text:
Anaplasma (formerly Ehrlichia) phagocytophilum, Ehrlichia chaffeensis, and Neorickettsia (formerly Ehrlichia) sennetsu are intracellular vector-borne pathogens that cause human ehrlichiosis, an emerging infectious disease. We present the complete genome sequences of these organisms along with comparisons to other organisms in the Rickettsiales order. Ehrlichia spp. and Anaplasma spp. display a unique large expansion of immunodominant outer membrane proteins facilitating antigenic variation. All Rickettsiales have a diminished ability to synthesize amino acids compared to their closest free-living relatives. Unlike members of the Rickettsiaceae family, these pathogenic Anaplasmataceae are capable of making all major vitamins, cofactors, and nucleotides, which could confer a beneficial role in the invertebrate vector or the vertebrate host. Further analysis identified proteins potentially involved in vacuole confinement of the Anaplasmataceae, a life cycle involving a hematophagous vector, vertebrate pathogenesis, human pathogenesis, and lack of transovarial transmission. These discoveries provide significant insights into the biology of these obligate intracellular pathogens.
Ehrlichiosis is an acute disease that triggers flu-like symptoms in both humans and animals. It is caused by a range of bacteria transmitted by ticks or flukes. Because these bacteria are difficult to culture, however, the organisms are poorly understood. The genomes of three emerging human pathogens causing ehrlichiosis were sequenced. A database was designed to allow the comparison of these three genomes to sixteen other bacteria with similar lifestyles. Analysis from this database reveals new species-specific and disease-specific genes indicating niche adaptations, pathogenic traits, and other features. In particular, one of the organisms contains more than 100 copies of a single gene involved in interactions with the host(s). These comparisons also enabled a reconstruction of the metabolic potential of five representative genomes from these bacteria and their close relatives. With this work, scientists can study these emerging pathogens in earnest.
Citation: Dunning Hotopp JC, Lin M, Madupu R, Crabtree J, Angiuoli SV, Eisen J, et al. (2006) Comparative Genomics of Emerging Human Ehrlichiosis Agents. PLoS Genet 2(2): e21. https://doi.org/10.1371/journal.pgen.0020021
Editor: Paul M. Richardson, The US DoE Joint Genome Institute, United States of America
Received: October 20, 2005; Accepted: January 9, 2006; Published: February 17, 2006
Copyright: © 2006 Dunning Hotopp. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This project was supported by National Institutes of Health grant R01 AI47885 to YR. In addition, JCDH was partially supported by National Science Foundation grant EF-0328363.
Competing interests: The authors have declared that no competing interests exist.
Abbreviations: HGA, human granulocytic anaplasmosis; HGE, human granulocytic ehrlichiosis; HME, human monocytic ehrlichiosis; NER, nucleotide excision repair; OMP, outer membrane protein; ORF, open reading frame
Anaplasma phagocytophilum, Ehrlichia chaffeensis, and Neorickettsia sennetsu are small (approximately 0.4–1.5 μm), pleomorphic α-Proteobacteria. These bacteria are human pathogens that replicate in membrane-bound compartments inside host granulocytes (A. phagocytophilum) or monocytes/macrophages (E. chaffeensis and N. sennetsu) [1–3]. They are obligate intracellular pathogens with a life cycle that involves both vertebrate and invertebrate hosts. A. phagocytophilum and E. chaffeensis depend on hematophagous ticks as vectors and wild mammals as reservoir hosts (Table 1) [2,4]. Unknown trematodes are suspected to be the vector and reservoir of N. sennetsu . No vaccine exists for any of these human pathogens.
Biological Characteristics of the Rickettsiales
A. phagocytophilum is the causative agent of human granulocytic anaplasmosis (HGA), formerly recognized as human granulocytic ehrlichiosis (HGE) [5,6]. Infection with A. phagocytophilum causes fever, headache, myalgia, anorexia, and chills . Prior to 1994, only ruminant and equine ehrlichiosis were known to be caused by this organism . A. phagocytophilum is transmitted by Ixodes spp. Cases of HGA correspond to the distribution of Ixodes spp. being identified in New England, the mid-Atlantic region, the upper Midwest, and northern California in the United States, as well as in parts of Europe. A. phagocytophilum is one of the leading causes of ehrlichiosis in the world. Recent serological data suggest that as much as 15%–36% of the population in endemic areas has been infected . Far fewer individuals are diagnosed with a symptomatic infection that varies in severity from fever to death . Half of all symptomatic patients require hospitalization, and 5%–7% require intensive care .
Human monocytic ehrlichiosis (HME), caused by E. chaffeensis, was discovered in 1986 [9–11]. HME is a systemic disease indistinguishable from HGA . E. chaffeensis has been most commonly identified in the Lone Star tick (Amblyomma americanum), with white-tailed deer considered to be the major reservoir. Over 500 cases of HME were diagnosed from 1986 to 1997, predominantly in the south-central and southeastern United States . The recognition and increased prevalence of the disease has been proposed to be related to changes in the host-vector ecology . As with all emerging diseases, it is likely outbreaks occurred in the preceding decades. Notably, 1,000 troops training in Texas contracted an unexplained disease with similar symptoms after exposure to the vector from 1942 to 1943 .
N. sennetsu is a monocytotropic species that causes sennetsu ehrlichiosis, an infectious mononucleosis-like disease with fever, fatigue, general malaise, and lymphadenopathy [1,13]. Less is known about the distribution of N. sennetsu when compared to Anaplasma and Ehrlichia. However, sequencing of its genome allows for interesting comparisons, since tissue tropism and clinical symptoms are similar but the vector (unknown trematodes) is different. Additionally, in the United States and Canada, domestic animals infected with the closely related N. risticii develop Potomac horse fever, an acute febrile disease accompanied by diarrhea with high morbidity and mortality [14,15]. The related N. helminthoeca causes acute and highly fatal salmon-poisoning disease of domestic and wild canines [14,16].
Along with Wolbachia, these bacteria are members of the Anaplasmataceae family (Figure 1) . Wolbachia infect arthropods and filarial nematodes, but have not been shown to infect vertebrates directly.
The protein sequences of select conserved genes were concatenated and aligned, and a phylogenetic tree was inferred of all sequenced α-Proteobacteria (see Materials and Methods). The Anaplasmataceae (purple) and the Rickettsiaceae (yellow) are highlighted.
Together with the Rickettsiaceae, the Anaplasmataceae are members of the order Rickettsiales (Figure 1) . The Rickettsiaceae include the obligate intracellular Rickettsia spp. Like the Anaplasmataceae, the Rickettsiaceae are obligate intracellular pathogens with a life cycle that involves both vertebrate and invertebrate hosts, but they replicate directly in the cytosol of endothelial cells. All organisms in the order Rickettsiales have relatively small genomes (0.8–1.5 Mb) that have arisen through reductive evolution as they developed dependence on the host cell for necessary functions . The Rickettsiales and other α-Proteobacteria also have an unresolved evolutionary relationship with the progenitor of the mitochondria [18,19].
Three Rickettsiaceae genomes have been published: Rickettsia prowazekii, R. conorii, and R. typhi [18,20,21]. Four Anaplasmataceae genomes have been published: the insect parasite W. pipientis wMel, the filarial nematode endosymbiont Wolbachia sp. wBm, the bovine pathogen Anaplasma marginale, and the bovine pathogen Ehrlichia ruminantium [19,22–24].
We present here a comparison of the previously completed Rickettsiales genomes to the first complete genomes of three representative Anaplasmataceae human pathogens: A. phagocytophilum, E. chaffeensis, and N. sennetsu. The complete genome sequence of these human pathogens will enhance the opportunities for investigation of virulence factors, pathogenesis, immune modulation, and novel targets for antimicrobial therapy and vaccines.
A. phagocytophilum, E. chaffeensis, and N. sennetsu each have a single circular chromosome (Figure S1). Most genomic features are typical of the sequenced Rickettsiales (Table 2). W. pipientis wMel, Ehrlichia spp., and Anaplasma spp., which are most closely related, all have numerous repeats in their genomes. In contrast, N. sennetsu and R. prowazekii have only six repeats in their respective genomes (Table 2). The repetitive nature of the Ehrlichia and Anaplasma genomes is exemplified by the expansion of outer membrane proteins of the OMP-1/P44/Msp2 family (discussed below). In addition numerous other functionally important genes are duplicated including those involved in type IV secretion and vitamin/cofactor biosynthesis.
The origin of replication was not experimentally determined in any of the genomes. As with other Rickettsiales , genes typically clustered near the origin (dnaA, gyrA, gyrB, rpmH, dnaN, parA, and parB) were dispersed throughout the genomes. For E. chaffeensis and N. sennetsu, a clear shift in GC-skew occurs near parA and parB (Figure 2). Therefore, basepair 1 was set in the intergenic region between the two genes. In A. phagocytophilum, none of these genes were found near the GC-skew transition. Therefore, basepair 1 was set in the intergenic region near polA. For E. chaffeensis and A. phagocytophilum, these predictions coincide with the predictions for E. ruminantium  and A. marginale .
Regions of conserved synteny were identified using the ortholog clusters (see Materials and Methods) and visualized with Sybil. The genes along each ordered chromosome were colored on a gradient from yellow to blue. The ortholog clusters for each query genome were then plotted relative to the order of the reference genome. Regions of synteny are then seen as continuous gradients across large regions of the genome. Above the synteny gradient display is the atypical nucleotide composition. Below the gradient display are the predicted coding regions on the plus strand and the minus strand, and the GC-skew. Representatives of all the Rickettsiales (A) and representative Ehrlichia spp. and Anaplasma spp. (B) were compared separately. AMA, A. marginale St. Maries; APH, A. phagocytophilum HZ; ECH, E. chaffeensis Arkansas; ERU, E. ruminantium Welgevonden; NES, N. sennetsu Miyayama; RPR, R. prowazekii Madrid E; WOL, W. pipientis wMel.
Only three islands of synteny over 10 kb in length are conserved among all the sequenced Anaplasmataceae, and these islands are shared among all the Rickettsiales (Figure 2). They include two operons of ribosomal proteins and one operon of proteins encoding portions of the type IV secretion system. Similar to the other Rickettsiales sequenced, all three genomes have the equivalent of a single rRNA operon with the 16S rRNA separated from the 23S-5S gene pair, as previously described for this order of bacteria .
Of genes typically clustered near the origin, parA and parB were not identified in A. phagocytophilum. Likewise, parA and parB are truncated in the Wolbachia sp. wBm. In various mutational studies in free-living prokaryotes, the effects of inactivation of parA or parB range from lethality to production of anucleated cells at low copy number [25,26]. Without parA and parB, A. phagocytophilum and the Wolbachia sp. wBm may have random chromosome partitioning, may require an alternate partitioning factor, or may have inefficient chromosome partitioning.
Of all the sequenced Anaplasmataceae, only the Anaplasma spp. and Ehrlichia spp. share conserved gene order (synteny) across their chromosome (Figure 2). E. ruminantium and E. chaffeensis have a single symmetrical inversion near two duplicate Rho termination factors (Figure 3). Symmetrical inversions around the origin are the most common large-scale rearrangements in microbial genomes . Genomic rearrangements between these Rho termination factors are also apparent in A. marginale. The presence of the same break in both the Anaplasma and Ehrlichia lineages suggests that the duplicate Rho termination factors allow for repeated inversions across this region of the genome.
Anaplasma spp. and Ehrlichia spp. share conserved gene order (synteny) across their chromosomes. E. ruminantium and E. chaffeensis have a single symmetrical inversion near two duplicate Rho termination factors (approximate positions shown in pink). Genomic rearrangements between these Rho termination factors are also apparent in A. marginale (pink). In addition to the synteny breaks near the Rho termination factors, A. marginale has rearrangements located near the msp2- and msp3- expression locus and pseudogenes (approximate positions shown in light blue). Likewise, in A. phagocytophilum, numerous changes in genome arrangement are located near the homologous p44 expression locus and silent genes (approximate positions shown in lavender).
In addition to the synteny breaks near the Rho termination factors, A. marginale has rearrangements located near the msp2 and msp3 expression loci and their corresponding pseudogenes (Figure 3). Likewise, numerous boundaries of genome rearrangements are located near the homologous p44 expression locus (p44ES/APH_1221) and silent genes. In both Anaplasma spp., the silent p44 and msp2 genes stored in reserve in the genome can recombine into the corresponding expression locus to generate antigenic variation in the immunodominant surface protein (discussed in detail below). These exact, repeated sequences throughout the genome facilitate recombination for antigenic variation and may also provide sites where chromosomal inversions occur.
In order to compare the genomic content of the Rickettsiales to that of other intracellular bacteria, ortholog clusters were delineated for 19 representatives of obligate and facultative intracellular pathogens and endosymbionts (see Materials and Methods). Such comparisons show conservation of 176 ortholog clusters across these intracellular bacteria (Table S1), most of which correspond to housekeeping functions.
Eleven ortholog clusters present in all the Rickettsiales distinguish the Rickettsiales from other intracellular bacteria examined (Table S2). These include a type I secretion system ATPase, a pyridine nucleotide-disulfide oxidoreductase family protein, a putative transporter, and type IV secretion system proteins VirB9 and VirB8. Thirteen ortholog clusters composed of 12 conserved hypothetical proteins and a GNAT family acetyltransferase distinguish all the Anaplasmataceae from the Rickettsiales (Table S3).
Five genera in the Rickettsiales order have at least one representative sequenced. In order to compare these five genera, the following genomes were compared: R. prowazekii, N. sennetsu, W. pipientis, A. phagocytophilum, and E. chaffeensis. This comparison shows conservation of 423 ortholog clusters (Table S4) generally associated with housekeeping functions. Most genes in the five compared genomes are either conserved among all genomes or unique to a given genome. Indeed, 60% of the two-, three-, and four-way comparisons shared fewer than ten ortholog clusters (Figure 4). In the three-way comparisons, the BDE (A. phagocytophilum, E. chaffeensis, and N. sennetsu) and CDE (A. phagocytophilum, E. chaffeensis, and W. pipientis) intersections harbor more than 20 ortholog clusters (Figure 4). The BDE intersection includes the organisms sequenced here and represents the human pathogens with very similar disease outcomes. Ortholog clusters conserved between these organisms include those for vitamin and cofactor biosynthesis enzymes, a monovalent cation/proton antiporter, a dicarboxylate transporter, and a DNA-binding protein (Table S5). Vitamin and cofactor biosynthesis is specific to the human ehrlichiosis agents, suggesting a niche adaptation or pathogenic trait. The CDE intersection is composed of the most closely related organisms. These ortholog clusters include genes for amino acid, fatty acid and nucleotide biosynthesis, an M48 family peptidase, a cytochrome c-type biogenesis protein, and the type IV secretion system protein VirB4 (Table S6).
The composition of ortholog clusters (see Materials and Methods) of representative Rickettsiales (A), Ehrlichia spp. (B), and Anaplasma spp. (C) were compared. Numbers within the intersections of different ovals indicate ortholog clusters shared by 2, 3, 4, or 5 organisms. Species compared are indicated in diagram intersections as follows. A, R. prowazekii; B, N. sennetsu; C, W. pipientis; D, A. phagocytophilum; E, E. chaffeensis; F, A. marginale; G, E. ruminantium Gardel; and H, E. ruminantium Welgevonden.
In two-way comparisons, the AC (R. prowazekii and W. pipientis) and DE (A. phagocytophilum and E. chaffeensis) intersections contain more than twenty ortholog clusters. Genes shared only by R. prowazekii and W. pipientis include those for cell wall biosynthesis, subunits of cytochrome D ubiquinol oxidase, a biotin transporter, a dinucleoside polyphosphate hydrolase, and an amino acid permease (Table S7). The presence of genes for cell wall biosynthesis in only R. prowazekii and W. pipientis likely reflects differences in the cell surface; A. phagocytophilum, E. chaffeensis, and N. sennetsu do not synthesize peptidoglycan . The peptidoglycan biosynthesis genes are also found in A. marginale, which suggests that if these genes are expressed, A. marginale may have a peptidoglycan layer . Since the peptidoglycan genes are present in A. marginale and W. pipientis but not in the other Anaplasmataceae, these genes have either been horizontally acquired in these organisms or have been lost numerous times in the Anaplasmataceae. Peptidoglycan binding to the Toll-like receptor 2 activates leukocytes. Neither A. marginale nor W. pipientis infects the immune cells of a vertebrate host. The peptidoglycan layer may have been lost to allow the organism to successfully infect vertebrate immune cells.
Genes shared only by A. phagocytophilum and E. chaffeensis include those encoding thiamine biosynthetic proteins, a potassium transporter, a peptide deformylase, and an ankyrin repeat protein (Table S8). Thiamine biosynthesis is distinctly absent from N. sennetsu, suggesting a possible trematode niche-specific adaptation.
A. phagocytophilum, E. chaffeensis, and N. sennetsu have 462, 312, and 303 open reading frames (ORFs) or paralog clusters that are unique with respect to the five-organism ortholog cluster analysis, respectively. The vast majority of these unique genes encode hypothetical, conserved hypothetical, and conserved domain proteins, as well as uncharacterized membrane proteins and lipoproteins. Other A. phagocytophilum-specific genes include those encoding the P44 outer membrane proteins and the HGE-14 and HGE-2 antigenic proteins (Table S9). E. chaffeensis-specific genes include those for the OMP-1 family of proteins, arginine biosynthesis, a major facilitator family transporter, and a variable-length PCR target protein (Table S10). N. sennetsu-specific genes include those for an F-type ATPase beta subunit, a cyclophilin-type peptidyl-prolyl cis-trans isomerase, a branched-chain amino acid transporter, a sensor histidine kinase, a strain-specific surface antigen, thioredoxin, and the type IV secretion system proteins VirB2 and VirB4 (Table S11).
Of the organism-specific genes detected in this five-way comparison, over half were hypothetical proteins, many of which formed genomic islands of hypothetical proteins (Figure 2). The majority of the genes identified as unique were not just unique to the genus, but to the species. Of the 462 A. phagocytophilum-unique genes in this comparison, 448 are also unique when compared with A. marginale. The 21 ortholog clusters shared only between Anaplasma spp. include conserved hypothetical proteins, OMP-1 proteins, membrane proteins, and HGE-2 (Table S12). Likewise, of the 312 E. chaffeensis-unique ORFs or paralog clusters in the five-way comparison, 267 are unique upon comparison with either E. ruminantium strain. The 52 ortholog clusters shared only between the Ehrlichia spp. include OMP-1 proteins, arginine biosynthetic proteins, a pyrroline-5-carboxylate reductase, a major facilitator protein, conserved hypothetical proteins, membrane proteins, and lipoproteins (Table S13).
Only one ortholog cluster containing conserved hypothetical proteins is shared between the animal pathogens E. ruminantium (Erum1840, ERGA_CDS_01780) and A. marginale (AM279) and are absent from the human pathogens E. chaffeensis, A. phagocytophilum, and N. sennetsu. In addition, a homolog of these proteins is present in the Ehrlichia canis Jake publicly available shotgun sequence. Since A. phagocytophilum and E. chaffeensis are maintained in animal reservoirs, presence of this gene is not associated with animal infection. Instead, loss of this protein could be required to establish infection in humans. These conserved hypothetical proteins have some homology to the eukaryotic patatin family of phospholipases. Patatin has been characterized to have phospholipase A-like activity .
Except for N. sennetsu, all of the sequenced pathogenic Anaplasmataceae require an arthropod-vector that feeds on blood (Table 1). Three ortholog clusters, including one for bacterioferritin and two for conserved hypothetical proteins, are absent in all of the tick-, flea-, and louse-borne Rickettsiales, but are present in Wolbachia spp. and N. sennetsu (Table S14). The proteins in these ortholog clusters may be correlated to the lack of a blood-sucking arthropod in the life cycles of these organisms.
The tick-borne Anaplasmataceae (Ehrlichia spp. and Anaplasma spp.) are the only Rickettsiales that are not transmitted transovarially in the invertebrate host. One ortholog cluster containing a class II aldolase/adducing domain protein (NSE_0849, RC0678, RP493, RT0479, WD0208) is absent only from Ehrlichia spp. and Anaplasma spp. Lack of this aldolase/adducing domain protein may prevent transovarial transmission in the arthropod vector.
Four ortholog clusters of conserved hypothetical proteins are present in all the pathogenic Rickettsiales but none of the endosymbionts. These proteins, which remain to be characterized, may be essential for pathogenesis or survival in the vertebrate host (Table S15).
A. phagocytophilum Strain Comparison
As an initial effort to use these genome sequences to identify the conserved genomic content of unsequenced members of these species, we conducted microarray-based comparative genome hybridization analyses with two A. phagocytophilum strains. Except for four p44 hypervariable regions (discussed below), the genomic content across all three strains is conserved (ratio < 3). Although A. phagocytophilum and A. marginale have very different complements of unique genes, the genomic content within the strains of A. phagocytophilum is highly conserved. Conservation of the gene content of the strains may explain the similarity of clinical signs of HGA from two geographic regions (New York, Minnesota) and equine ehrlichiosis in California .
Free-Living and Obligate-Intracellular α-Proteobacteria
In order to understand the differences between these obligate intracellular pathogens and a closely related free-living organism, the number of genes in each role category was compared between representative Anaplasmataceae and Caulobacter crescentus (Table 3). C. crescentus is a closely related and sequenced free-living α-Proteobacteria to the Rickettsiales . The scope of this comparison was limited to only these five α-Proteobacteria, as only these organisms had role categories assigned in an identical manner.
Comparison of Role Category Breakdown
All of the Anaplasmataceae examined have significantly higher percentages of their genomes involved in nucleotide biosynthesis, cofactor and vitamin biosynthesis, and protein synthesis. Enzymes in these biosynthetic pathways are likely to play an important role in interactions with their hosts and intracellular survival, as discussed below. The protein synthesis category includes many essential genes such as those encoding ribosomal proteins, tRNA synthetases, RNA modification enzymes, and translation factors. These genes are essential and cannot be sacrificed as the genome reduces. Therefore, as the genome size decreases, the proportion of genes involved in protein synthesis increases.
All of the Anaplasmataceae examined have a significantly lower coding capacity for central intermediary metabolism, transport, and regulatory functions. The decrease in central intermediary metabolism and transport reflects the differences in acquiring nutrients and energy. Since intracellular bacteria are exposed to a relatively restricted complement of nutrients and energy sources, they have evolved to be specialists in acquiring specific compounds from their hosts.
Likewise, these intracellular bacteria live in a homeostatic environment and have fewer regulatory genes. ORFs encoding σ70 and σ32 were identified (rpoD and rpoH, respectively), but σ24 and σ54 were not detected (rpoE and rpoN, respectively). Several two-component regulatory systems are retained and may be employed as these bacteria transition between their vertebrate and invertebrate hosts. Despite being identified in Rickettsia spp. , stringent response (a global regulatory response) may not be expected in the Anaplasmataceae, since neither RelA nor SpoT proteins were identified.
There are several role categories in which only specific organisms have significant differences from, or similarities to, C. crescentus. All the bacteria except E. chaffeensis have a statistically significant decrease in amino acid biosynthesis. The difference between Ehrlichia spp. and the other Anaplasmataceae is due to the presence of lysine and arginine biosynthesis pathways in Ehrlichia spp., as discussed below. A. phagocytophilum has a significant increase in the percentage of genes dedicated to the cell envelope due to expansion of the OMP-1 family in Anaplasma spp. (discussed below). W. pipientis has a significantly higher percentage of its genome involved in mobile and extrachromosomal functions due to the unique presence of phage and transposons in its genome . E. chaffeensis, A. phagocytophilum, and N. sennetsu have a significant decrease in mobile elements, as they have no intact prophage, no transposable elements, and only a few phage core components (HK97-like portal, major capsid, and prohead protease) scattered throughout their genomes. Lastly, A. phagocytophilum and W. pipientis both have an increased number of disrupted reading frames.
Based on comparisons of the intracellular and free-living α-Proteobacteria, the only overall theme that emerges is the conservation of housekeeping genes and the shuffling of the genomes resulting in the loss of many operon structures.
Little is known about the genetic determinants required for the Rickettsiales to invade a host and cause disease. Putative pathogenesis genes were identified, including enzymes to neutralize reactive oxygen species, outer membrane proteins, and protein secretion systems.
Oxidative stress response.
Reactive oxygen species have been implicated in both host defense to infection and host cell injury [31–33]. All of the Rickettsiales contain sodB, an iron superoxide dismutase. This superoxide dismutase may have an important role in pathogenesis since sodB is cotranscribed with components of the type IV secretion system in E. chaffeensis and A. phagocytophilum .
Further examination of conserved genes without functional annotation (e.g., conserved hypothetical proteins, conserved domain proteins) shows two other ortholog clusters of proteins that may be involved in response to oxidative stress—a putative heme copper oxidase and a putative flavohemoglobin. In both cases, there is no significant similarity to a protein of known function, but several conserved domains were identified. From a particular combination of domains and conservation of metal/cofactor ligands, a function of response to oxidative stress can be proposed for these proteins .
Indeed, ECH_1079, NSE_0121, and APH_1205 each contain the 12 transmembrane segments and six conserved histidine residues consistent with members of the heme-copper oxidase family. Members of this protein family include cytochrome oxidase subunit I, FixN for nitrogen fixation, and NorB for nitric oxide reduction . Each of these organisms is unlikely to be fixing nitrogen and already has a functional subunit I of cytochrome oxidase (ECH_1003, NSE_0622, and APH_1085), so these orthologs may be nitric oxide reductases. Alternatively, there may be another, as yet to be identified, role for this oxidase, which was identified in all the Rickettsiales genomes except the Wolbachia sp. wBm where it is truncated (an ORF that was not annotated but has genomic coordinates from 536343 to 536534).
APH_0545, NSE_0661, and ECH_0778 encode proteins with three functional motifs similar to flavohemoglobins—a heme binding site, an FAD binding domain, and an NAD binding domain. The biological function of the Escherichia coli flavohemoglobin has not been elucidated, but it has been shown to be an efficient alkylhydroperoxide reductase  and a nitric oxide reductase . This putative flavohemoglobin is conserved among the Anaplasmataceae, but Wolbachia spp. are missing the NAD oxidoreductase domain, and R. prowazekii is missing the heme ligands. Although the speculation of a role for these genes in pathogenicity is intriguing, the precise function of each of these proteins will need to be elucidated experimentally.
The OMP-1/MSP2/P44 protein superfamily.
The Anaplasmataceae all have a diverse complement of outer membrane proteins. Many of these outer membrane proteins (OMPs) are members of Pfam PF01617  and constitute the OMP-1/MSP2/P44 family. Anaplasma, Ehrlichia, and Wolbachia have each undergone variable levels of expansion of their omp-1/msp2 gene families (Figure S2). The N. sennetsu genome has only one uncharacterized protein from this family (NSE_0875). W. pipientis wMel and the Wolbachia sp. wBm have the smallest expansion with three wsp genes scattered throughout each genome. The largest expansion of this family is in Ehrlichia spp. and Anaplasma spp. These organisms cannot be transovarially inherited in their arthropod hosts. Instead, ticks acquire Ehrlichia or Anaplasma by feeding on an infected vertebrate reservoir animal. The expansion of this family may allow persistence in the vertebrate reservoir by providing antigenic variation, thus allowing for effective tick transmission.
E. chaffeensis, E. canis, and E. ruminantium have 17–22 paralogous tandemly arranged genes from this family that are flanked by a transcription regulator (tr1) and a preprotein translocase (secA) [40–42]. These genes all have signal peptides and are likely to be secreted across the cytoplasmic membrane by SecA . They encode immunodominant major outer membrane proteins that are differentially expressed in ticks and experimentally infected animals .
A. marginale St. Maries is reported to have 56 genes that have been placed into this superfamily, including eight msp2, eight msp3, one msp4, three opag, 15 omp-1, 12 orfX, seven orfY, and two msp3 remnants . These genes are scattered throughout the genome with a bias in location toward the origin of replication. MSP2 and MSP3 are the immunodominant proteins . The msp2 and msp3 gene subsets each include one full-length expression locus and seven reserve/silent sequences that are thought to recombine into the expression locus to generate antigenic variation .
The A. phagocytophilum genome has three omp-1, one msp2, two msp2 homologs, one msp4, and 113 p44 loci belonging to the OMP-1/MSP2/P44 superfamily. Although both Anaplasma spp. msp2 genes are members of PF01617 and the OMP1/MSP2/P44 superfamily, the A. marginale msp2 gene is distinct from the A. phagocytophilum msp2 gene. In addition, the previously identified omp-1N is not a member of this Pfam, but is homologous to E. chaffeensis omp-1N and the msp2 operon-associated gene 3 of A. marginale .
The largest expansion of this family is that of p44 genes in A. phagocytophilum. Only 36 copies of p44 are in this Pfam, but many smaller regions were identified, resulting in a total of 113 annotated p44 loci (Table S16). The p44s consist of a central hypervariable region of approximately 280 bp containing a signature of four conserved amino acid regions (C, C, WP, A) and conserved flanking sequences longer than 50 bp. Diverse p44 paralogs (p44–1 to p44–65) are expressed in mammals and ticks and confer antigenic environmental adaptation, especially during tick transmission [46–49]. The genomic loci of all 65 previously described p44 genes were determined in the present study (Figure S3). Twenty-three novel p44 genes (p44-66 to p44–88) were identified by genome sequencing, but have not yet been experimentally identified as being expressed.
The p44s were annotated as full-length, silent/reserve, truncated, and fragments (Figure 5). There are 22 full-length p44s identified that have ORFs longer than 1.0 kb with conserved start and stop codons. By locating highly conserved 5′ and 3′ flanking sequences and signature sequences within the hypervariable region, 64 shorter p44s were identified. These ORFs lack a translational start codon and likely serve as reserve/silent p44s that can be expressed after recombining into the previously described p44-expression locus (p44ES/APH_1221) [45,50]. The full-length and silent/reserve p44 genes are preferentially located near the replication origin (Figure S3) and symmetrically located around the p44 expression locus. Localization near the origin, where multiple replication forks coexist, may facilitate recombination between the expression locus and the reserve/silent p44 genes.
Full-length p44 genes contain conserved start and stop codons, an ORF longer than 1,000 bp, and a central hypervariable region of approximately 280 bp containing a signature of four conserved amino acid regions (C, C, WP, A). These genes can be expressed at their respective current genome location or can recombine into the expression locus (p44ES/APH_1221). A silent/reserve p44 is less than 1,000 bp. It may have either the conserved or alternative start and/or stop codons. A silent/reserve p44 is not likely to be expressed at its current genome location, but can recombine into the expression locus (p44ES/APH_1221). Truncated p44s carry the complete hypervariable region, or a portion thereof, but only one of the two conserved regions. Fragments of p44 have only a conserved region and no hypervariable region. Each annotated p44 is longer than 60 bp. It should be noted that smaller fragments can be identified throughout the genome. These, as well as p44 truncations and fragments, are likely to be nonfunctional remnants of previous recombination events.
In addition to the full-length and silent/reserve p44 genes, 21 5′ and 3′ fragments and six truncations of p44 genes larger than 60 nucleotides have been identified in the genome. Truncations include portions of a hypervariable region; fragments did not. The p44s annotated as truncated and fragments do not contain both conserved regions flanking the hypervariable region. These p44s are not expected to recombine through the homologous recombination model deduced by previous analyses of recombined p44s [49–52].
Microarray-based comparative genomic hybridization reveals that expansion of the p44 family is a common feature in A. phagocytophilum strains. All but four of the p44 unique hypervariable sequences used as targets on the microarray are present in the human isolate A. phagocytophilum MN and the horse isolate A. phagocytophilum California MRK (Figure S3; Table 4). The p44-12 and p44-9 unique regions are either absent or divergent only in strain MN. The p44-4 and p44-1 unique regions are absent or divergent in strains MN and MRK. This confirms previous results demonstrating that the p44–1 unique region is absent/divergent in MN and MRK .
Other important outer membrane proteins.
N. sennetsu has a single p51 gene (NSE_0242) encoding its immunodominant P51 major outer membrane protein . The p51 gene is highly conserved among N. risticii, N. sennetsu, and the Stellantchasmus falcatus fluke agent, but not in N. helminthoeca, the agent causing an acute, highly fatal salmon-poisoning disease of domestic and wild canines . Although a full-length, highly conserved homolog for P51 was not found in the Rickettsiales genome sequences, P51 was placed in an ortholog cluster of genes conserved among all the Rickettsiales due to short regions of similarity, particularly in a C-terminal region that may include a secretion peptide motif.
Other outer membrane proteins have been reported in A. marginale, including msp5, msp1a, and msp1b. The msp5 gene (a SCO1/SenC family protein) is found in all the Rickettsiales, whereas msp1a and msp1b are unique to A. marginale.
Only E. chaffeensis and E. canis encode a 120-kDa immunodominant surface protein (ECH_0039) . The variable-length PCR target useful in distinguishing various strains of E. chaffeensis  is present only in the genome of E. chaffeensis Arkansas (ECH_0170).
Protein secretion systems.
All of the strains sequenced here contain both a Sec-dependent and Sec-independent protein export pathway for secretion of proteins across the inner membrane. The Sec-independent pathway (Tat pathway) has been implicated in the transport of phospholipases in Pseudomonas aeruginosa . All of the strains sequenced here also contain two components of a putative type I secretion system, potentially for transporting toxins or proteases carrying a C-terminal secretion signal.
All of the Rickettsiales have a type IVa secretion system that uses a complex of transmembrane proteins and a pilus to deliver effector macromolecules from prokaryotic to eukaryotic cells. The reference Type IVa secretion system is that of Agrobacterium tumefaciens, which contains 11 genes in the virB locus and one gene in the virD locus. Several components of the A. tumefaciens type IVa secretion system are conserved in A. phagocytophilum, E. chaffeensis, and N. sennetsu. Like R. prowazekii and W. pipientis, the three organisms sequenced here are lacking virB1, virB5, and virB7. All but N. sennetsu lack virB2.
The virB3, virB4, and virB6 homologs are contiguous at one locus (Figure S4). Neighboring this locus in all of these organisms are three or four virB6 homologs. Contiguous at a second locus are virB8, virB9, virB10, virB11, and virD4. The type IV secretion system is one of the few sets of genes syntenic between all of the Rickettsiales sequenced, suggesting that tight coordination of expression of these genes is critical.
In A. tumefaciens, translocated type IV effector proteins have the consensus sequence R-X7-R-X-R-X-R-X-Xn, where lysine can substitute for arginine with no noticeable effect . In addition, effector molecules are often localized to a region of the chromosome near the type IV secretion apparatus. Examination of the regions around the type IV operons in A. phagocytophilum revealed numerous genes encoding HGE-14, which contain C-terminal sequences similar, but not identical, to this motif (Table S17), suggesting that it may be an excreted effector molecule. Subsequent searches of the Anaplasmataceae genomes with motifs like that found in HGE-14 did not reveal other potential effector molecules.
The metabolic potentials of A. phagocytophilum, E. chaffeensis, and N. sennetsu were compared to that of R. prowazekii and W. pipientis [18,19]. Overall, the Anaplasmataceae have very similar metabolic pathways but are quite distinct from those of R. prowazekii (Figure 6). W. pipientis differs from the other Anaplasmataceae in its inability to synthesize some cofactors.
Metabolic pathways of E. chaffeensis (magenta arrows), A. phagocytophilum (green arrows), N. sennetsu (gold arrows), W. pipientis (lavender arrows), and R. prowazekii (cyan arrows) were reconstructed and compared. The networks of some of the more important pathways are shown with metabolites color coded: red and purple, central and intermediary metabolites; blue, cofactors; green, amino acids; and black, cell structures. Transporters are shown in the membrane and are grouped by predicted substrate specificity: green, inorganic cations; magenta, inorganic anions; red, carbohydrates and carboxylates; blue, amino acids/peptides/amines; yellow, nucleotides/nucleosides; and black, drug/polysaccharide efflux or unknown.
Nucleotide and cofactor biosynthesis.
E. chaffeensis, A. phagocytophilum, N. sennetsu, and W. pipientis have the ability to synthesize all nucleotides. This differs from R. prowazekii, which cannot make purines or pyrimidines, and therefore must rely on nucleotide translocases and interconversion of the bases to obtain the full complement of nucleotides . E. chaffeensis, A. phagocytophilum, and N. sennetsu are able to synthesize most vitamins and cofactors. In contrast to the other Anaplasmataceae, W. pipientis has lost some of its ability to synthesize cofactors, and it has completely lost the biosynthetic pathways for biotin, thiamine, and NAD. In addition, it may be in the process of losing the ability to synthesize folate. R. prowazekii has also lost the ability to synthesize these cofactors as well as FAD, pantothenate, and pyridoxine-phosphate.
Biotin is one of the essential cofactors only synthesized by the vertebrate-infecting Anaplasmataceae. In most organisms, biotin is required for many carboxylation reactions, but is not synthesized by many multicellular eukaryotes. RT-PCR analysis showed that all four genes in the biotin biosynthesis pathway (BioA/B/D/F) were expressed by E. chaffeensis and A. phagocytophilum in THP-1 and HL-60 cells, respectively, at both 2 d and 3 d post infection (Figure S5).
The presence of nucleotide, vitamin, and cofactor biosynthetic pathway in E. chaffeensis, A. phagocytophilum, and N. sennetsu suggests that they do not need to compete with the host cell for, and may even supply host cells with, essential vitamins and nucleotides. It has been previously proposed that Wigglesworthia glossinidia supplies its host with vitamins that are rare in the blood meal of its arthropod host (tsetse fly) . Interestingly, Ehrlichia spp. and Anaplasma spp., the two tick-borne intracellular pathogens sequenced, both have a complement of pathways for cofactor and amino acid biosynthesis similar to W. glossinidia (Table 5). This raises the possibility that these pathogens may currently be, or historically have been, able to provide a benefit to their tick hosts by providing necessary cofactors.
Amino acid biosynthesis.
The Rickettsiales have a very limited ability to synthesize amino acids and must rely on transporting them from the host (Figure 6). All four of the Anaplasmataceae sequenced have the ability to make glycine, glutamine, glutamate, and aspartate. Additionally, E. chaffeensis is predicted to be able to synthesize arginine and lysine like E. ruminantium . One possible role for arginine biosynthesis may be to recover an intracellular arginine pool after exposure to inducible host nitric oxide. Nitric oxide is synthesized by nitric oxide synthases that convert arginine to citrulline and nitric oxide . The production of nitric oxide is likely to deplete the intracellular pool of arginine, further hampering intracellular growth. The presence of an arginine biosynthesis pathway and putative nitric oxide reductase(s) may allow Ehrlichia spp. to recover more rapidly and subvert the host immune response. This would be similar to the proposed retention of select tryptophan biosynthetic genes in Chlamydia spp. in order to replenish tryptophan pools after host enzymatic degradation of tryptophan in response to IFN-γ .
Glycolysis, tricarboxylic acid cycle, pentose phosphate, and respiration.
A complete pyruvate dehydrogenase, tricarboxylic acid cycle, F0F1-ATPase, and electron transport chain were found in all of the organisms. All five organisms are likely to use host-derived carboxylates and amino acids, but none of these organisms can obtain carbon or energy from fatty acids or actively carry out glycolysis. The glycolysis enzymes present are limited to those that produce glyceraldehyde-3-phosphate and dihydroxyacetone phosphate from phosphoenolpyruvate (Figure 6). The glyceraldehyde-3-phosphate produced in this manner is used in the nonoxidative pentose phosphate pathway, resulting in the production of pentoses needed for cofactor and nucleotide biosynthesis. Consistent with this role for the glycolytic enzymes, R. prowazekii and R. conorii retain neither the glycolytic enzymes nor the enzymes needed for the biosynthesis of nucleotides or cofactors from pentose. Similarly, dihydroxyacetone phosphate from these glycolytic enzymes can be converted to glycerol-3-phosphate for phospholipid biosynthesis in the Anaplasmataceae. Without the glycolytic enzymes, Rickettsia spp. must obtain glycerol-3-phosphate from the host via a glycerol-3-phosphate transporter.
Evolution and DNA Repair
A genome-scale phylogenetic analysis using a concatenated alignment of core proteins is consistent with rRNA studies and current taxonomic assignments. This indicates that Anaplasma and Ehrlichia are sister genera that share a common ancestor with Wolbachia (Figure 1). Neorickettsia is the deepest-branching lineage in the group.
The branch lengths on the whole genome tree can be used to get an indication of the relative rates of evolution of these organisms. In general, the branch lengths for these intracellular organisms are longer than those of their free-living relatives. This may be due to either differences in DNA repair or population genetic and selection-related force. For example, many intracellular organisms go through more stringent population bottlenecks, which in turn increase the amount of genetic drift and possibly the rate of accumulation of deleterious mutations.
Analysis of the genome of W. pipientis wMel revealed that it had a longer branch length than the closely related Rickettsia; the Rickettsia have higher rates of evolution than free-living organisms . Wu et al.  ascribed this increase to features of Wolbachia biology. However, there appears to be a general increase in the rate for all of the Anaplasmataceae (Figure 1). Thus, the increase reported for Wolbachia  is not likely due to the specific biology of Wolbachia, but instead to some feature shared by all Anaplasmataceae.
Examination of the putative DNA-repair capabilities of the different species does not reveal any significant differences between the Anaplasmataceae and the Rickettsia spp. (Table S18). Interestingly, within the Anaplasmataceae, N. sennetsu appears to have the longest branch length and the most limited suite of DNA repair genes within the group. For example, N. sennetsu is missing various glycosylases and exonucleases that contribute to repair, including uvrABC, which is involved in nucleotide excision repair. It is possible that the faster rate of evolution in this organism is related to the absence of some of these repair pathways.
The absence of uvrABC in N. sennetsu and the absence of uvrBC in the Ehrlichia spp. suggest that these species do not have nucleotide excision repair (NER). NER is used by other organisms, including bacteria, archaea, and eukaryotes, as a general repair process to remove sections of DNA with gross abnormalities. One important role of NER is in the repair of UV-induced DNA damage, and defects in NER in other species lead to great increases in UV sensitivity. It appears that Neorickettsia has compensated for this by acquiring a gene homologous to DNA photolyases, an alternative mechanism for repairing UV damage. The Neorickettsia photolyase is not particularly closely related to known photolyases from α-Proteobacteria but is instead most closely related to a photolyase from Coxiella burnetii, a γ-Proteobacteria. The Ehrlichia spp., however, do not encode a photolyase homolog, and thus these species may be highly UV-sensitive.
The dual existence of members of Anaplasma spp. and Ehrlichia spp. as invertebrate symbionts or commensals and effective human and animal pathogen requires flexibility, a fact reflected in the genome. Both organisms display an expansive inventory of paralogous genes encoding diverse functions that promote survival and success in different environments when compared to Neorickettsia spp. and Wolbachia spp., which do not require a mammalian host. This capacity is evident from the large repertoire of outer membrane proteins, and partial duplication of some of the virulence determinants (e.g., components of the type IV protein secretion system).
The large number of paralogous genes encoding immunodominant outer membrane proteins in Anaplasma spp. and Ehrlichia spp. has important implications for the study of pathogenesis and in the development of vaccination strategies. Adaptability in the human host may underlie significant disease manifestations. Genomic-level characterization of the full complement of variable antigens will facilitate the future development of more specific and sensitive diagnostic targets. In light of the growing recognition of the increased global burden of ehrlichiosis, development of such diagnostic targets will impact public health.
Between pairwise comparisons of different species within a single genus, there are hundreds of genes that are not shared. Often these gene differences are immunodominant outer membrane proteins, but the vast majority are genes that are not functionally characterized in any organism. Some are likely to be involved in zoonosis or specific disease characteristics. For instance, A. phagocytophilum is the only sequenced Rickettsiale that infects neutrophils. Therefore, some of the A. phagocytophilum-unique genes (e.g., genes encoding P44 and HGE-14) may be involved in neutrophil invasion.
Many pathogens are obligate intracellular bacteria. But since they are difficult or impossible to culture and tools for genetic manipulation are limited, they are less well characterized than the facultative intracellular bacteria or extracellular pathogens. The analysis of the genome sequences provides critical insights into the biology of these intracellular pathogens and will facilitate manipulation of the emerging human ehrlichiosis agents and leukocytotropic pathogens.
Materials and Methods
Intracellular bacteria purification and DNA preparation.
Organisms (infecting ~1 × 109 host cells; 50–100 175-cm2 flasks) were cultured in synchrony in respective host cells (E. chaffeensis in DH82 cells, A. phagocytophilum in HL-60 cells, and N. sennetsu in P388D1 cells). Bacterial cells were liberated from the infected host cells using Dounce homogenization, differential centrifugation, and Percoll density gradient centrifugation . Any specimens with host nuclei contamination were excluded. From these isolated bacteria, phenol extraction was used to purify DNA that was minimally fragmented and free of host-cell DNA. Levels of host DNA contamination were verified to be less than 0.001% by PCR using host G3PDH-specific primers. This method was highly successful, with only 14 sequencing reads identified as being of human origin from a total of over 57,000 good sequencing reads.
Sequencing and annotation.
The complete genome sequences were determined using the whole-genome shotgun sequencing approach , sequences were assembled into contigs using the Celera Assembler , and all gaps were closed . ORFs from each genome were predicted and annotated using a suite of automated tools that combine Glimmer gene prediction [64,65], ORF and non-ORF feature identification (e.g., protein motifs), and assignment of database matches and functional role categories to genes . Frameshifts and point mutations were detected and corrected where appropriate; those remaining were annotated as “authentic frameshift” or “authentic point mutation.” Repeats were identified using RepeatFinder [66,67] and were manually curated. The complete genome sequences for A. phagocytophilum HZ, E. chaffeensis Arkansas, and N. sennetsu Miyayama have been deposited in GenBank.
Annotation of the p44 genes.
Full-length p44s were defined as having ORFs greater than 1,000 bp with conserved start codon and stop codons. For shorter silent/reserve p44s, the ORFs were initially identified by locating highly conserved 5′ and 3′ sequences and signature sequences within the hypervariable region. Since these silent/reserve p44s lack a start and stop codon, the 5′ and 3′ ends were annotated on the basis of conserved genome features found in full-length p44 genes [50,68]. The annotated p44 fragments are at least 60 nucleotides in length, have either 5′ or 3′ conserved sequences, and may contain a partial hypervariable region (Figure 5).
Ortholog clusters were delineated for R. prowazekii Madrid E , R. typhi Wilmington , R. conorii Malish 7 , N. sennetsu Miyayama, Wolbachia sp. wBm , W. pipientis wMel , E. chaffeensis Arkansas, E. ruminantium Gardel (GenBank CR925677.1), E. ruminantium Welgevonden , A. marginale St. Maries , A. phagocytophilum HZ, Brucella suis 1330 , Bartonella henselae Houston-1 , Coxiella burnetii RSA 493 , Tropheryma whipplei Twist , Blochmannia floridanus , Buchnera sp. APS , Chlamydia pneumoniae AR39 , and W. glossinidia brevipalpis . For Wolbachia sp. wBm, the ORFs used in these comparisons were uncurated ORFs predicted and annotated using a suite of automated tools that combine Glimmer gene prediction [64,65], ORF and non-ORF feature identification (e.g., protein motifs), and assignment of database matches and functional role categories to genes . Upon release of the annotated genome , these uncurated ORFs were paired with the corresponding curated ones where possible, with exceptions noted in the text.
Paralog clusters within each of the genomes were identified using the Jaccard algorithm with the following parameters: 80% or greater identity and Jaccard coefficient 0.6 or higher [77, Text S1]. Members of paralog clusters were then organized into ortholog clusters by allowing any member of a paralog cluster to contribute to the reciprocal best matches used to construct the ortholog clusters. The conservation of ortholog clusters across the various genomes analyzed was determined using Sybil, a web-based software package for comparative genomics developed at TIGR (http://sybil.sourceforge.net). The database of these clusters and corresponding tools can be accessed through TIGR (http://www.tigr.org/sybil/rcd). Metabolic pathways and transporters were compared across genomes using (1) these calculated ortholog clusters, (2) Genome Properties , (3) TransportDB , and (4) Biocyc .
Significant differences in the role category composition was determined using χ2 calculated using the Yates continuity correction. A p-value less than 0.01 was considered significant.
GC-skew and origin prediction.
The GC-skew was calculated as (C − G)/(C + G) in windows of 1,000 bp along the chromosome . The origin of replication was not experimentally determined in any of the genomes. For E. chaffeensis and N. sennetsu, a clear shift in GC-skew occurs near parA and parB. Therefore basepair 1 was set in the intergenic region between the two genes. In A. phagocytophilum, a GC-skew transition occurs near polA. Therefore, basepair 1 was set in the intergenic region near polA.
Atypical nucleotide composition.
Regions of atypical nucleotide composition were identified by the χ2 analysis: the distribution of all 64 trinucleotides was computed for the complete genome in all six reading frames, followed by the trinucleotide distribution in 5,000-bp windows overlapping by 500 bp. For each window, the χ2 statistic was computed based on the difference between the trinucleotide content in that window and that of the whole genome. Peaks indicate regions of atypical nucleotide composition.
Genome tree construction.
Protein sequences of 31 housekeeping genes (frr, infC, nusA, pgk, pyrG, rplA, rplB, rplC, rplD, rplE, rplF, rplK, rplL, rplM, rplN, rplP, rplS, rplT, rpmA, rpoB, rpsB, rpsC, rpsE, rpsI, rpsJ, rpsK, rpsM, rpsS, smpB, and tsf) from complete α-Proteobacteria genomes were aligned to predefined HMM models and ambiguous regions were autotrimmed according to an embedded mask. Concatenated alignments were then used to build a maximum likelihood tree with bootstrapping using PHYML . The γ-Proteobacteria E. coli and the β-Proteobacteria Neisseria meningitidis were used as outgroups to root the tree.
Array construction and hybridizations.
Oligonucleotides (70-mer) were designed from the unique ORFs of each of the three genomes. The oligonucleotides (Illumina, San Diego, California, United States) were diluted to 25 μM in DMSO and spotted in quadruplicate onto UltraGap slides (Corning, Acton, Massachusetts, United States). Cy3 and Cy5 probes were synthesized from genomic DNA as previously described . In order to obtain enough DNA for microarray analysis, small amounts of DNA were prepared in the manner described above for genome sequencing. This DNA was then quantitatively amplified using GenomiPhi (Amersham, Piscataway, New Jersey, United States).
Appropriately labeled query and reference probes were hybridized overnight, washed, and scanned using an Axon GenePix 4000B scanner (Axon Instruments, Union City, California, United States). The corresponding images were analyzed with TIGR Spotfinder . Log mode centering was used to normalize the data alleviating the bias of expression microarray normalization methods, which expect a normal distribution of data. Briefly, a Perl script was designed to construct the histogram of the log2 of the ratio and adjust the histogram mode to zero. The data presented are the geometric means of the normalized ratios from at least two slides with different reference Cy dyes and with oligonucleotides printed in quadruplicate.
Transcript analysis of biotin biosynthetic genes.
Total RNA was extracted from E. chaffeensis or A. phagocytophilum-infected THP-1 or HL-60 cells at 2 d or 3 d postinfection using RNeasy (Qiagen, Valencia, California, United States). RNA was DNase I treated (Invitrogen, Carlsbad, California, United States) in the presence of 40 U of RNaseOUT (Invitrogen) for 15 min at room temperature, followed by inactivation at 65 °C in the presence of 2.5 mM EDTA for 10 min. For cDNA synthesis, total RNA (0.5 μg) was reverse-transcribed at 42 °C for 1 h in 50 mM Tris-HCl (pH 8.3), 75 mM KCl, 3 mM MgCl2, 0.5 mM of each dNTP, 1 U of RNase inhibitor (Invitrogen), 1.5 μM random hexamers (Invitrogen), and 10 U of Superscript II reverse transcriptase (Invitrogen). The reaction was terminated by heat inactivation at 70 °C for 15 min. To ensure the absence of DNA contamination in the RNA preparations, the assay was duplicated without reverse transcriptase. The subsequent amplification was conducted with standard conditions for 25 cycles of 95 °C for 45 s, 54 °C for 45 s, and 72 °C for 1 min and with the PCR primer pair (Table S19).
Figure S1. Linear Representation of E. chaffeensis Arkansas, A. phagocytophilum HZ, and N. sennetsu Miyayama Genomes
ORFs are oriented along the molecule and are color-coded by role category: violet, amino acid biosynthesis; light blue, biosynthesis of cofactors, prosthetic groups, and carriers; light green, cell envelope; red, cellular processes; brown, central intermediary metabolism; yellow, DNA metabolism; light gray, energy metabolism; black, mobile/extrachromosomal functions and truncated ORFs; magenta, fatty acid and phospholipid metabolism; pink, protein synthesis and fate; orange, purines, pyrimidines, nucleosides, and nucleotides; olive, regulatory functions and signal transduction; dark green, transcription; teal, transport and binding proteins; gray, unknown function; crosshatched, conserved hypothetical proteins; white, hypothetical proteins; salmon, other categories. The tRNA genes and rRNA genes are represented by their secondary structure. Repeats are shown by dual-arrowed line segments; paralogs are represented by arrows.
(10.4 MB PDF)
Figure S2. Phylogenetic Tree of OMP1 Proteins
The protein sequences of all the members of PFAM01617 were aligned and a phylogenetic tree inferred. The divergence of the OMP1/MSP2/P44 proteins in this superfamily did not permit robust inferences about the evolution of these proteins, but allowed classification of the proteins into superfamilies as reflected in their annotation. Particular families within this superfamily are highlighted, including the P44 proteins (pink), the OMP-1s (blue), and the Wsp (yellow).
(876 KB PDF)
Figure S3. Distribution of the p44 Genes in A. phagocytophilum Strains
Rims 1 and 2: predicted coding regions on the plus and minus strands, respectively, color-coded by role categories. Rim 3: atypical nucleotide composition. Rim 4: distribution of p44 silent genes (green), expression locus (cyan), full-length genes (magenta), fragments (brown), and truncations (blue). Rims 5 and 6: microarray-based comparative genome hybridization results for A. phagocytophilum strain HZ against strains MRK and MN. The ratios [(HZ normalized intensity)/(query normalized intensity)] were divided into three categories: ratio > 10 (red; absent); ratio 3–10 (blue; absent/divergent); and ratio < 3 (not plotted; present).
(1.0 MB PDF)
Figure S4. Type IV Secretion Systems in Rickettsiales
Genes encoding the type IV secretion system components can be found at two distinct regions of the Rickettsiales genome. At the larger of these regions, virB3, virB4, and virB6 show a typical arrangement. These are followed by a series of genes in the virB6 family that have been shown to be cotranscribed in W. pipientis wMel. Each of these regions is presented with ortholog clusters (see Materials and Methods) and color coded: cyan, virB3; orange, virB4; green, virB6; and purple, a virB6 family of genes. Orthologs conserved in location are connected with gray bars. The virB3 gene is not always annotated, due to its small size, but it is present in all Rickettsiales genomes examined.
(616 KB PDF)
Figure S5. Transcript Analysis of Biotin Biosynthetic Genes
DNase-treated total RNA was reverse-transcribed and subsequently PCR amplified using primers specific to each biotin biosynthesis gene. RT-PCR analysis showed that all four genes in the biotin biosynthesis pathway (BioA/B/D/F) were expressed by E. chaffeensis and A. phagocytophilum in THP-1 and HL-60 cells, respectively, at 2 d (unpublished data) and 3 d postinfection.
(61 KB JPG)
Table S1. Ortholog Clusters Conserved across All Representative Obligate and Facultative Intracellular Pathogens and Endosymbionts Presented in a Tab-Delimited Format
(640 KB DOC)
Table S2. Ortholog Clusters Present in All the Rickettsiales but Not in Any Other Intracellular Bacterium Examined
(49 KB DOC)
Table S3. Ortholog Clusters Present in All Anaplasmataceae but Not the Rickettsiales or Other Intracellular Bacterium Examined
(47 KB DOC)
Table S4. Ortholog Clusters Present in All of Five Representatives of the Genera in the Rickettsiales That Have at Least One Representative Sequenced
The following genomes were compared: Rickettsia prowazekii, Neorickettsia sennetsu, Wolbachia pipientis, Anaplasma phagocytophilum, and Ehrlichia chaffeensis.
(424 KB DOC)
Table S5. Ortholog Clusters Present in Neorickettsia sennetsu, Anaplasma phagocytophilum, and Ehrlichia chaffeensis, but Not in Wolbachia pipientis and Rickettsia prowazekii
(44 KB DOC)
Table S6. Ortholog Clusters Present in Wolbachia pipientis, Anaplasma phagocytophilum, and Ehrlichia chaffeensis, but Not in Neorickettsia sennetsu and Rickettsia prowazekii
(54 KB DOC)
Table S7. Ortholog Clusters Present in Rickettsia prowazekii and Wolbachia pipientis, but Not in Neorickettsia sennetsu, Anaplasma phagocytophilum, or Ehrlichia chaffeensis
(41 KB DOC)
Table S8. Ortholog Clusters Present in Anaplasma phagocytophilum and Ehrlichia chaffeensis, but Not in Rickettsia prowazekii, Wolbachia pipientis, or Neorickettsia sennetsu
(46 KB DOC)
Table S9. Individual Genes (Based on Ortholog Cluster Analysis) Present in Anaplasma phagocytophilum, but Not in Ehrlichia chaffeensis, Rickettsia prowazekii, Wolbachia pipientis, or Neorickettsia sennetsu
(394 KB DOC)
Table S10. Individual Genes (Based on Ortholog Cluster Analysis) Present in Ehrlichia chaffeensis but Not in Anaplasma phagocytophilum, Rickettsia prowazekii, Wolbachia pipientis, or Neorickettsia sennetsu
(244 KB DOC)
Table S11. Individual Genes (Based on Ortholog Cluster Analysis) Present in Neorickettsia sennetsu but Not in Anaplasma phagocytophilum, Rickettsia prowazekii, Wolbachia pipientis, or Ehrlichia chaffeensis
(223 KB DOC)
Table S12. Ortholog Clusters Present Only in Anaplasma phagocytophilum and Anaplasma marginale
(57 KB DOC)
Table S13. Ortholog Clusters Present Only in Ehrlichia chaffeensis, Ehrlichia ruminantium Welgevonden, and Ehrlichia ruminantium Gardel
(82 KB DOC)
Table S14. Ortholog Clusters That Are Absent in All of the Tick-, Flea-, and Louse-Borne Rickettsiales, but Are Present in Wolbachia spp. and Neorickettsia sennetsu
(27 KB DOC)
Table S15. Ortholog Clusters Present in All the Pathogenic Rickettsiales but None of the Endosymbionts
(30 KB DOC)
Table S16. Anaplasma phagocytophilum p44 Genes
(246 KB DOC)
Table S17. Putative Anaplasma phagocytophilum Type IV Effector Motifs in HGE-14 Proteins
(27 KB DOC)
Table S18. DNA Repair Genes
(110 KB DOC)
Table S19. Sequences of Oligonucleotides Used for RT-PCR
(31 KB DOC)
Text S1. PDF copy of Jaccard et al. reference
(4578 KB PDF)
The GenBank (http://www.ncbi.nlm.nih.gov) accession number for the Ehrlichia canis Jake shotgun sequence is ZP_00210380; the complete GenBank genome sequences for A. phagocytophilum HZ, E. chaffeensis Arkansas, and N. sennetsu Miyayama are CP000235, CP000236, and CP00237, respectively.
The ArrayExpress database (http://www.ebi.ac.uk/arrayexpress) accession numbers for the microarray slide type and study are E-TIGR-125 and A-TIGR-21.
We acknowledge Jessie Goodman, University of Minnesota, for providing the stock culture of A. phagocytophilum MN; Robert F. Massung at the Centers for Disease Control and Prevention, Atlanta for providing the stock culture of A. phagocytophilum MRK; Derrick Fouts for examination of the genomes for prophage elements; Robert DeBoy for examination of the genome for transposons; Leka Papazisi for assistance with statistical analysis; Karen Nelson, Ian Paulsen, and Emmanuel Mongodin for helpful discussions; Chun-Ha Wan for assistance depositing data in ArrayExpress; Robert Munson at Children's Hospital, Ohio State University for assistance with the grant proposal; David Dyer at University of Oklahoma for initial sequencing efforts; and our reviewers for their helpful comments and suggestions.
H. Tettelin coordinated genome sequencing. Y. Rikihisa performed general project coordination and strategy. J. C. Dunning Hotopp, M. Lin, S. V. Angiuoli, J. Crabtree, J. Sundaram, Y. Rikihisa, and H. Tettelin performed sequence analysis and comparative genomics. C. Zhang, N. Ohashi, H. Niu, N. Zhi, and M. Lin cultured bacteria, purified bacteria, and purified DNA. R. Madupu, L. M. Brinkac, Q. Lin, R. J. Dodson, M. J. Rosovitz, W. Nelson, S. C. Daugherty, A. S. Durkin, M. Gwinn, D. H. Haft, J. D. Selengut, S. A. Sullivan, L. Zhou, and O. White annotated the genomes. S. V. Angiuoli, J. Crabtree, J. Sundaram, M. Lin, T. Davidsen, N. Zafar, and O. White managed data. T. R. Utterback, S. Smith, M. Lewis, H. Khouri, F. Benahmed, H. Forberger, R Halpin, S. Mulligan, and J. Robinson sequenced and finished the genomes. J. Eisen, R. Seshadri, Q. Ren, and M. Wu performed specific analyses. M. Lin and J. C. Dunning Hotopp performed RT-PCR and CGH experiments, respectively. J. C. Dunning Hotopp, Y. Rikihisa, M. Lin, and H. Tettelin wrote the paper.
- 1. Rikihisa Y (1991) The tribe Ehrlichieae and ehrlichial diseases. Clin Microbiol Rev 4: 286–308.
- 2. Rikihisa Y (2003) Mechanisms to create a safe haven by members of the family Anaplasmataceae. In: Hechemy KE, Zvsic-Zupanc T, Childs JE, Rault DA, editors. New York: New York Academy of Sciences pp. 548–555.
- 3. Dumler JS, Barbet AF, Bekker CP, Dasch GA, Palmer GH, et al. (2001) Reorganization of genera in the families Rickettsiaceae and Anaplasmataceae in the order Rickettsiales: Unification of some species of Ehrlichia with Anaplasma, Cowdria with Ehrlichia and Ehrlichia with Neorickettsia, descriptions of six new species combinations and designation of Ehrlichia equi and “HGE agent” as subjective synonyms of Ehrlichia phagocytophila. Int J Syst Evol Microbiol 51: 2145–2165.
- 4. Telford SR 3rd, Dawson JE, Katavolos P, Warner CK, Kolbert CP, et al. (1996) Perpetuation of the agent of human granulocytic ehrlichiosis in a deer tick-rodent cycle. Proc Natl Acad Sci U S A 93: 6209–6214.
- 5. Chen SM, Dumler JS, Bakken JS, Walker DH (1994) Identification of a granulocytotropic Ehrlichia species as the etiologic agent of human disease. J Clin Microbiol 32: 589–595.
- 6. Goodman JL, Nelson C, Vitale B, Madigan JE, Dumler JS, et al. (1996) Direct cultivation of the causative agent of human granulocytic ehrlichiosis. N Engl J Med 334: 209–215.
- 7. Dumler JS, Bakken JS (1998) Human ehrlichioses: Newly recognized infections transmitted by ticks. Annu Rev Med 49: 201–213.
- 8. Dumler JS, Choi K-S, Garcia-Garcia JC, Barat NS, Scorpio DG, et al. (2005) Human granulocytic anaplasmosis and Anaplasma phagocytophilum. Emerg Infect Dis 11: 1828–1834.
- 9. Anderson BE, Dawson JE, Jones DC, Wilson KH (1991) Ehrlichia chaffeensis, a new species associated with human ehrlichiosis. J Clin Microbiol 29: 2838–2842.
- 10. Dawson JE, Anderson BE, Fishbein DB, Sanchez JL, Goldsmith CS, et al. (1991) Isolation and characterization of an Ehrlichia sp. from a patient diagnosed with human ehrlichiosis. J Clin Microbiol 29: 2741–2745.
- 11. Maeda K, Markowitz N, Hawley RC, Ristic M, Cox D, et al. (1987) Human infection with Ehrlichia canis, a leukocytic rickettsia. N Engl J Med 316: 853–856.
- 12. Paddock CD, Childs JE (2003) Ehrlichia chaffeensis: A prototypical emerging pathogen. Clin Microbiol Rev 16: 37–64.
- 13. Tachibana N, Shishime E, Okayama A, Ishizaki J, Murai K, et al. (1987) [Two cases of spotted fever rickettsiosis in Kyushu]. Kansenshogaku Zasshi 61: 1166–1172.
- 14. Rikihisa Y, Zhang C, Kanter M, Cheng Z, Ohashi N, et al. (2004) Analysis of p51, groESL, and the major antigen P51 in various species of Neorickettsia, an obligatory intracellular bacterium that infects trematodes and mammals. J Clin Microbiol 42: 3823–3826.
- 15. Gibson KE, Rikihisa Y, Zhang C, Martin C (2005) Neorickettsia risticii is vertically transmitted in the trematode Acanthatrium oregonense and horizontally transmitted to bats. Environ Microbiol 7: 203–212.
- 16. Rikihisa Y, Zimmerman G (1995) Salmon poisoning disease. In: Kirk RW, Bonagura JD, editors. Current veterinary therapy XII small animal practice. Philadelphia: W. B. Saunders. pp. 297–300.
- 17. Andersson SG, Kurland CG (1998) Reductive evolution of resident genomes. Trends Microbiol 6: 263–268.
- 18. Andersson SG, Zomorodipour A, Andersson JO, Sicheritz-Ponten T, Alsmark UC, et al. (1998) The genome sequence of Rickettsia prowazekii and the origin of mitochondria. Nature 396: 133–140.
- 19. Wu M, Sun LV, Vamathevan J, Riegler M, Deboy R, et al. (2004) Phylogenomics of the reproductive parasite Wolbachia pipientis wMel: A streamlined genome overrun by mobile genetic elements. PLoS Biol 2: e69.
- 20. Ogata H, Audic S, Renesto-Audiffren P, Fournier PE, Barbe V, et al. (2001) Mechanisms of evolution in Rickettsia conorii and R. prowazekii. Science 293: 2093–2098.
- 21. McLeod MP, Qin X, Karpathy SE, Gioia J, Highlander SK, et al. (2004) Complete genome sequence of Rickettsia typhi and comparison with sequences of other Rickettsiae. J Bacteriol 186: 5842–5855.
- 22. Foster J, Ganatra M, Kamal I, Ware J, Makarova K, et al. (2005) The Wolbachia genome of Brugia malayi: Endosymbiont evolution within a human pathogenic nematode. PLoS Biol 3: e121.
- 23. Brayton KA, Kappmeyer LS, Herndon DR, Dark MJ, Tibbals DL, et al. (2005) Complete genome sequencing of Anaplasma marginale reveals that the surface is skewed to two superfamilies of outer membrane proteins. Proc Natl Acad Sci U S A 102: 844–849.
- 24. Collins NE, Liebenberg J, de Villiers EP, Brayton KA, Louw E, et al. (2005) The genome of the heartwater agent Ehrlichia ruminantium contains multiple tandem repeats of actively variable copy number. Proc Natl Acad Sci U S A 102: 838–843.
- 25. Draper GC, Gober JW (2002) Bacterial chromosome segregation. Annu Rev Microbiol 56: 567–597.
- 26. Lewis RA, Bignell CR, Zeng W, Jones AC, Thomas CM (2002) Chromosome loss from par mutants of Pseudomonas putida depends on growth medium and phase of growth. Microbiology 148: 537–548.
- 27. Eisen JA, Heidelberg JF, White O, Salzberg SL (2000) Evidence for symmetric chromosomal inversions around the replication origin in bacteria. Genome Biol 1: RESEARCH0011.
- 28. Lin M, Rikihisa Y (2003) Ehrlichia chaffeensis and Anaplasma phagocytophilum lack genes for lipid A biosynthesis and incorporate cholesterol for their survival. Infect Immun 71: 5324–5331.
- 29. Hirschberg HJ, Simons JW, Dekker N, Egmond MR (2001) Cloning, expression, purification and characterization of patatin, a novel phospholipase A. Eur J Biochem 268: 5037–5044.
- 30. Nierman WC, Feldblyum TV, Laub MT, Paulsen IT, Nelson KE, et al. (2001) Complete genome sequence of Caulobacter crescentus. Proc Natl Acad Sci U S A 98: 4136–4141.
- 31. Carlyon JA, Abdel-Latif D, Pypaert M, Lacy P, Fikrig E (2004) Anaplasma phagocytophilum utilizes multiple host evasion mechanisms to thwart NADPH oxidase-mediated killing during neutrophil infection. Infect Immun 72: 4772–4783.
- 32. Park J, Rikihisa Y (1992) L-arginine-dependent killing of intracellular Ehrlichia risticii by macrophages treated with gamma interferon. Infect Immun 60: 3504–3508.
- 33. Feng HM, Whitworth T, Popov V, Walker DH (2004) Effect of antibody on the Rickettsia-host cell interaction. Infect Immun 72: 3524–3530.
- 34. Ohashi N, Zhi N, Lin Q, Rikihisa Y (2002) Characterization and transcriptional analysis of gene clusters for a type IV secretion machinery in human granulocytic and monocytic ehrlichiosis agents. Infect Immun 70: 2128–2138.
- 35. Namboori S, Mhatre N, Sujatha S, Srinivasan N, Pandit SB (2004) Enhanced functional and structural domain assignments using remote similarity detection procedures for proteins encoded in the genome of Mycobacterium tuberculosis H37Rv. J Biosci 29: 245–259.
- 36. Saraste M, Castresana J (1994) Cytochrome oxidase evolved by tinkering with denitrification enzymes. FEBS Lett 341: 1–4.
- 37. Poole RK, Hughes MN (2000) New functions for the ancient globin family: Bacterial responses to nitric oxide and nitrosative stress. Mol Microbiol 36: 775–783.
- 38. Bonamore A, Gentili P, Ilari A, Schinina ME, Boffi A (2003) Escherichia coli flavohemoglobin is an efficient alkylhydroperoxide reductase. J Biol Chem 278: 22272–22277.
- 39. Bateman A, Coin L, Durbin R, Finn RD, Hollich V, et al. (2004) The Pfam protein families database. Nucleic Acids Res 32: D138–141.
- 40. Yu X, McBride JW, Zhang X, Walker DH (2000) Characterization of the complete transcriptionally active Ehrlichia chaffeensis 28 kDa outer membrane protein multigene family. Gene 248: 59–68.
- 41. van Heerden H, Collins NE, Brayton KA, Rademeyer C, Allsopp BA (2004) Characterization of a major outer membrane protein multigene family in Ehrlichia ruminantium. Gene 330: 159–168.
- 42. Ohashi N, Rikihisa Y, Unver A (2001) Analysis of transcriptionally active gene clusters of major outer membrane protein multigene family in Ehrlichia canis and E. chaffeensis. Infect Immun 69: 2083–2091.
- 43. Unver A, Rikihisa Y, Stich RW, Ohashi N, Felek S (2002) The omp-1 major outer membrane multigene family of Ehrlichia chaffeensis is differentially expressed in canine and tick hosts. Infect Immun 70: 4701–4704.
- 44. Vidotto MC, McGuire TC, McElwain TF, Palmer GH, Knowles DP Jr. (1994) Intermolecular relationships of major surface proteins of Anaplasma marginale. Infect Immun 62: 2940–2946.
- 45. Barbet AF, Meeus PF, Belanger M, Bowie MV, Yi J, et al. (2003) Expression of multiple outer membrane protein sequence variants from a single genomic locus of Anaplasma phagocytophilum. Infect Immun 71: 1706–1718.
- 46. Felek S, Telford S 3rd, Falco RC, Rikihisa Y (2004) Sequence analysis of p44 homologs expressed by Anaplasma phagocytophilum in infected ticks feeding on naive hosts and in mice infected by tick attachment. Infect Immun 72: 659–666.
- 47. Lin Q, Zhi N, Ohashi N, Horowitz HW, Aguero-Rosenfeld ME, et al. (2002) Analysis of sequences and loci of p44 homologs expressed by Anaplasma phagocytophila in acutely infected patients. J Clin Microbiol 40: 2981–2988.
- 48. Wang X, Rikihisa Y, Lai TH, Kumagai Y, Zhi N, et al. (2004) Rapid sequential changeover of expressed p44 genes during the acute phase of Anaplasma phagocytophilum infection in horses. Infect Immun 72: 6852–6859.
- 49. Zhi N, Ohashi N, Tajima T, Mott J, Stich RW, et al. (2002) Transcript heterogeneity of the p44 multigene family in a human granulocytic ehrlichiosis agent transmitted by ticks. Infect Immun 70: 1175–1184.
- 50. Lin Q, Rikihisa Y, Ohashi N, Zhi N (2003) Mechanisms of variable p44 expression by Anaplasma phagocytophilum. Infect Immun 71: 5650–5661.
- 51. Lin Q, Rikihisa Y (2005) Establishment of cloned Anaplasma phagocytophilum and analysis of p44 gene conversion within an infected horse and infected SCID mice. Infect Immun 73: 5106–5114.
- 52. Lin Q, Rikihisa Y, Massung RF, Woldehiwet Z, Falco RC (2004) Polymorphism and transcription at the p44–1/p44–18 genomic locus in Anaplasma phagocytophilum strains from diverse geographic regions. Infect Immun 72: 5574–5581.
- 53. Yu XJ, McBride JW, Diaz CM, Walker DH (2000) Molecular cloning and characterization of the 120-kilodalton protein gene of Ehrlichia canis and application of the recombinant 120-kilodalton protein for serodiagnosis of canine ehrlichiosis. J Clin Microbiol 38: 369–374.
- 54. Sumner JW, Childs JE, Paddock CD (1999) Molecular cloning and characterization of the Ehrlichia chaffeensis variable-length PCR target: An antigen-expressing gene that exhibits interstrain variation. J Clin Microbiol 37: 1447–1453.
- 55. Voulhoux R, Ball G, Ize B, Vasil ML, Lazdunski A, et al. (2001) Involvement of the twin-arginine translocation system in protein secretion via the type II pathway. EMBO J 20: 6735–6741.
- 56. Vergunst AC, van Lier MC, den Dulk-Ras A, Stuve TA, Ouwehand A, et al. (2005) Positive charge is an important feature of the C-terminal transport signal of the VirB/D4-translocated proteins of Agrobacterium. Proc Natl Acad Sci U S A 102: 832–837.
- 57. Zientz E, Dandekar T, Gross R (2004) Metabolic interdependence of obligate intracellular bacteria and their insect hosts. Microbiol Mol Biol Rev 68: 745–770.
- 58. Chakravortty D, Hensel M (2003) Inducible nitric oxide synthase and control of intracellular bacterial pathogens. Microbes Infect 5: 621–627.
- 59. Wood H, Roshick C, McClarty G (2004) Tryptophan recycling is responsible for the interferon-gamma resistance of Chlamydia psittaci GPIC in indoleamine dioxygenase-expressing host cells. Mol Microbiol 52: 903–916.
- 60. Ohashi N, Zhi N, Zhang Y, Rikihisa Y (1998) Immunodominant major outer membrane proteins of Ehrlichia chaffeensis are encoded by a polymorphic multigene family. Infect Immun 66: 132–139.
- 61. Fleischmann RD, Adams MD, White O, Clayton RA, Kirkness EF, et al. (1995) Whole-genome random sequencing and assembly of Haemophilus influenzae Rd. Science 269: 496–512.
- 62. Myers EW, Sutton GG, Delcher AL, Dew IM, Fasulo DP, et al. (2000) A whole-genome assembly of Drosophila. Science 287: 2196–2204.
- 63. Tettelin H, Feldblyum TV (2004) Genome sequencing and analysis. In: Grandi G, editor. Genomics, proteomics and vaccines. London: John Wiley and Sons. pp. 45–73. pp.
- 64. Salzberg SL, Delcher AL, Kasif S, White O (1998) Microbial gene identification using interpolated Markov models. Nucleic Acids Res 26: 544–548.
- 65. Delcher AL, Harmon D, Kasif S, White O, Salzberg SL (1999) Improved microbial gene identification with GLIMMER. Nucleic Acids Res 27: 4636–4641.
- 66. Volfovsky N, Haas BJ, Salzberg SL (2001) A clustering method for repeat analysis in DNA sequences. Genome Biol 2: RESEARCH0027.
- 67. Kurtz S, Schleiermacher C (1999) REPuter: Fast computation of maximal repeats in complete genomes. Bioinformatics 15: 426–427.
- 68. Zhi N, Ohashi N, Rikihisa Y (1999) Multiple p44 genes encoding major outer membrane proteins are expressed in the human granulocytic ehrlichiosis agent. J Biol Chem 274: 17828–17836.
- 69. Paulsen IT, Seshadri R, Nelson KE, Eisen JA, Heidelberg JF, et al. (2002) The Brucella suis genome reveals fundamental similarities between animal and plant pathogens and symbionts. Proc Natl Acad Sci U S A 99: 13148–13153.
- 70. Alsmark CM, Frank AC, Karlberg EO, Legault BA, Ardell DH, et al. (2004) The louse-borne human pathogen Bartonella quintana is a genomic derivative of the zoonotic agent Bartonella henselae. Proc Natl Acad Sci U S A 101: 9716–9721.
- 71. Seshadri R, Paulsen IT, Eisen JA, Read TD, Nelson KE, et al. (2003) Complete genome sequence of the Q-fever pathogen Coxiella burnetii. Proc Natl Acad Sci U S A 100: 5455–5460.
- 72. Raoult D, Ogata H, Audic S, Robert C, Suhre K, et al. (2003) Tropheryma whipplei Twist: A human pathogenic Actinobacteria with a reduced genome. Genome Res 13: 1800–1809.
- 73. Gil R, Silva FJ, Zientz E, Delmotte F, Gonzalez-Candelas F, et al. (2003) The genome sequence of Blochmannia floridanus: Comparative analysis of reduced genomes. Proc Natl Acad Sci U S A 100: 9388–9393.
- 74. Shigenobu S, Watanabe H, Hattori M, Sakaki Y, Ishikawa H (2000) Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS. Nature 407: 81–86.
- 75. Read TD, Brunham RC, Shen C, Gill SR, Heidelberg JF, et al. (2000) Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39. Nucleic Acids Res 28: 1397–1406.
- 76. Akman L, Yamashita A, Watanabe H, Oshima K, Shiba T, et al. (2002) Genome sequence of the endocellular obligate symbiont of tsetse flies, Wigglesworthia glossinidia. Nat Genet 32: 402–407.
- 77. Jaccard P (1908) Nouvelles recherches sur la distribution florale. Bull Soc Vaud Sci Nat 44: 223–270.
- 78. Haft DH, Selengut JD, Brinkac LM, Zafar N, White O (2005) Genome properties: A system for the investigation of prokaryotic genetic content for microbiology, genome annotation and comparative genomics. Bioinformatics 21: 293–306.
- 79. Ren Q, Kang KH, Paulsen IT (2004) TransportDB: A relational database of cellular membrane transport systems. Nucleic Acids Res 32: D284–288.
- 80. Krieger CJ, Zhang P, Mueller LA, Wang A, Paley S, et al. (2004) MetaCyc: A multiorganism database of metabolic pathways and enzymes. Nucleic Acids Res 32: D438–442.
- 81. Lobry JR (1996) Asymmetric substitution patterns in the two DNA strands of bacteria. Mol Biol Evol 13: 660–665.
- 82. Guindon S, Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52: 696–704.
- 83. Tettelin H, Nelson KE, Paulsen IT, Eisen JA, Read TD, et al. (2001) Complete genome sequence of a virulent isolate of Streptococcus pneumoniae. Science 293: 498–506.
- 84. Saeed AI, Sharov V, White J, Li J, Liang W, et al. (2003) TM4: A free, open-source system for microarray data management and analysis. Biotechniques 34: 374–378.


text:
Fierce African aesthetics and cutting edge design collide as the fashion industry welcomes a brand new label into its sphere this week. Imagined by Zayna Gummi, Élan Bijoux retails exclusive handcrafted pieces, each as unique as its wearer. Every item is crafted with luxurious materials such as stingray leather, gemstones and precious metals. The result is a resplendent collection of distinctive jewellery that’s guaranteed to turn heads.
This season the Élan Bijoux aura is all about bold, statement pieces that can be rocked alone or teamed up to create a daring set. Think gemstone encrusted cuffs, handmade stacked bangles, stingray leather spiral bracelets in a cornucopia of colours, statement spherical earrings and more.
Born and raised in Nigeria, Gummi grew up with an insatiable infatuation with fashion. From the age of 11 she was sketching out designs for her sisters, a hobby that went on to become a professional passion. In 2012 Gummi moved to the UK to pursue a career in law. She has since graduated from the University of Buckingham and while she retains a passion for combating serious issues such as terrorism, money laundering and corruption, she balances it out with self-built luxury jewellery brand, Élan Bijoux. Despite the fact that Nigerian culture often frowns upon girls pursuing ‘non-professional’ careers such as fashion design or make-up artistry, Gummi received unwavering support from her family.
Gummi says, “Why make a choice? I don’t have any plans to work exclusively in a single field, and I adore the way I go from arguing about serious issues one day to discussing stone choices and metal finishes the next. It really is the best of best worlds!”
Élan Bijoux was borne from a bespoke jewellery set Gummi designed for her sister’s wedding. When manufacturers refused to create just a single set she went all out and designed a debut collection that has gone on to dazzle women across the globe. Inspired by Nigeria’s rich regal heritage and colourful culture, as well as the excitement and innovation of contemporary trends, Élan Bijoux is a reflection of Gummi’s past, present and future. She designs for fearless women who are inherently modern, infectiously confident, have an eye for quality and are unafraid to step out in something new.
One of the standout lines is the unique stingray leather range manufactured in Thailand. In South East Asia stingrays are not only considered a valuable source of protein but are also renowned for their extremely durable skin. Élan Bijoux then builds on the natural strength of the skin with luxurious finishes to create a lavish range with a glossy shine and roughened texture.
Gummi already has her sights set on the future, with a British made fine jewellery collection set to be released later this year. Expansion is also on the cards, with luxurious leather handbags and shoes designs currently in the works. Élan Bijoux is ever evolving and fans are encouraged to regularly check the e‑boutique for new additions and seasonal pieces.
All orders over £100 are eligible for free shipping to UK addresses.
For a limited time only new customers are invited to redeem 10% off their first order by signing up for the Élan Bijoux newsletter.
To find out more about Élan Bijoux, and to shop the range of luxury jewellery available to purchase from the e‑boutique, go to: www.elanbijoux.com

domain: artistic


text:
An intruder stabbed and wounded five people at a rabbi’s house in New York during a party to celebrate the Jewish festival of Hanukkah late Saturday, officials said.
Rabbi’s house is a synagogue located in the Forshay neighbourhood in Monsey, Rockland County, New York, United States.
The victims, all Hasidic members of the Jewish faith, were transported to local hospitals — two in critical condition — the Orthodox Jewish Public Affairs Council (OJPAC) tweeted after receiving a call at 9.50 pm.
According to Jerusalem Post, the intruder pulled out a machete, which he used to stab people
A suspect has been taken into custody and a vehicle safeguarded, a spokesman for the NYPD told AFP.
Governor of New York, Andrew Cuomo, said in a statement that he was “horrified” by the “despicable and cowardly act”, and had directed the State Police hate crimes task force to…
Source: Vanguard Newspaper

domain: news


text: It\u2019s hard to be a millennial these days.\r\n\r\nEverywhere we turn, we see stories on how record levels of student loan debt are killing millennial\u2019s chances at homeownership. Meanwhile, money woes for millennials mean many are putting off marriage, quite possibly for good. And according to government figures, the unemployment rate for millennials still hovers near 15 percent \u2013 a statistic The Atlantic used last year to claim that millennials \u201chave it worse than any generation in the last 50 years.\u201d\r\n\r\nYes, millennials have it rough. They\u2019re living in their parent\u2019s basements, struggling to keep the cabinets stocked with Ramen, and wondering which move to make next.\r\n\r\nIt\u2019s all tragic and indisputable. Like it or not, being under thirty these days isn\u2019t easy \u2013 at least in an economic sense.\r\nThe Millennial Advantage No One Wants to Talk About\r\nBut here\u2019s the thing: Millennials have one huge advantage \u2013 one huge gift \u2013 that they almost refuse to recognize. This gift is something so valuable that it transcends every other item of value on Earth, and it\u2019s something that is finite, at least on an individual level.\r\n\r\nThat\u2019s right; millennials have time. Time to save. Time to pay off debt. Time to cultivate their passions, skillsets, and careers. Time to invest in the future. And yes, time to invest in a financial sense too. Because as we all know, the magic of compound interest is far easier to harness when you\u2019re young.\r\nSo, Why Aren\u2019t Millennials Investing?\r\nUnfortunately, a recent survey from Bankrate.com revealed that millennials aren\u2019t taking advantage of the one thing working in their favor \u2013 their youth. According to the survey, which polled 1,001 adults to gauge millennial attitudes on investing, only 26 percent of individuals under the age of 30 owned stock in early 2015.\r\n\r\nSurvey participants listed several main reasons why they choose not to invest, none of which are surprising. More than half of those who don\u2019t invest (53%) use a lack of financial resources as their excuse, while another 21 percent said they didn\u2019t know enough about stocks to take the plunge. The rest cited a general distrust of stockbrokers and worry over paying too much in fees. Related: See how much you are paying in fees with our Fee Analyzer.\r\n\r\nThose are all good reasons in theory, but disastrous when applied to the real world. It makes you wonder - what will happen to all the young people who fail to save and invest while they\u2019re young?\r\n\r\n\u201cYoung people who don't invest in equities early are set to have a lot less money later on," John Salter, associate professor of financial planning at Texas Tech University, told Bankrate.com after analyzing their survey results.\r\n\r\n\u201cThey won't benefit from those extra years of letting their portfolio grow, eventually earning the \u2018interest on interest\u2019 that generates real wealth,\u201d he said, adding that "it's just the simple value of compounding.\u201d\r\n\r\nIn other words, millennials who fail to invest won\u2019t have enough money to retire. Instead, they\u2019ll need to continue working indefinitely, subsist on government assistance and social security, or shack up with their own kids in old age.\r\nRetirement: It\u2019s the End Game\r\nIf you\u2019re a millennial, you might find those ideas preposterous, or even hilarious. But don\u2019t start laughing yet. Here are a few facts from the National Council on Aging that you might want to chew on:\r\n\r\nAccording to the NCOA and various government agencies:\r\n\u2022 23 million Americans over 60 are already living in poverty\r\n\u2022 75 percent of those 65+ rely on social security alone\r\n\u2022 One-third of senior households lives paycheck to paycheck\r\n\u2022 In 2012, the average credit card debt among adults aged 65+ was $9,283\r\n\u2022 14% of adults aged 65+ face retirement with negative net worth\r\n\r\nWhen you look at what actually happens to people who don\u2019t save or invest enough to enjoy a comfortable retirement, it isn\u2019t so funny. In fact, the idea of living in your son\u2019s basement at 70 stops being hilarious the moment you realize that it\u2019s not only possible, but actually likely.\r\n\r\nAnd maybe it\u2019s just me, but spending retirement in poverty sounds a whole lot worse than being 25 and having student loan debt and a job you hate. Because, no matter what, you still have time. You have time to stop and think \u2013 time to make a change. Time to tackle your student loan debt with fervor. Time to get a second or third job. Time to go back to school, finish your degree, write a book, or improve yourself in a multitude of ways.\r\n\r\nYour future self doesn\u2019t have that privilege. When you\u2019re 75 and subsisting on frozen Hungry Mans, it\u2019s the end game. You\u2019re stuck with whatever future you, yourself designed. And if you happen to look around and think \u201cthis sucks\u201d in old age, you\u2019ll only have yourself to blame.\r\nThe Investing Advice Every Millennial Needs to Hear\r\nAnd that\u2019s why it\u2019s more important than ever to start investing as soon as you can. You may think you can\u2019t afford to, but the truth is, you can\u2019t afford not to. Your future depends on what you do here and now, and your \u201cfuture self\u201d is begging \u2013 pleading \u2013 for you to take action not just now, but yesterday.\r\n\r\nThe internet is a treasure trove for those who want to learn investing basics. Meanwhile, Personal Capital offers free tools that can help you monitor your spending, net worth, and growing portfolio of investments. You may not have a lot to start with, but you have to start somewhere. Whether you choose to start by socking away $50 or $100 a month at first, investing in your company 401(k) up to the match, or opening a traditional or Roth IRA with a brokerage firm, it\u2019s essential that you take that first step. Because, no matter what, no one can take it for you.\r\n\r\nWhen you\u2019re in your twenties, you think the world is something that happens to you. You wake up each day with the goal of doing the best you can with what you have; you react to your surroundings. You shrink at opportunity. You see maybe one or two years out, or perhaps five years or a decade, but you just can\u2019t picture yourself as a little old lady, no matter how hard you try.\r\n\r\nBut the years will creep by faster the older you get, and soon you\u2019ll look back and wonder where the time went. And if you think investing is hard now, just wait a while. The longer you sit idle, more impossible the idea of retirement will become.\r\n\r\nSo here it is, the investing advice every millennial needs to hear:\r\n\r\nDon\u2019t use student loans or anything you hear on the 24-hour news cycle as an excuse not to save or invest for your future. Sixty years from now, nothing will matter except for what you did with the resources you had.\r\n\r\nYou have more power and potential than you realize, but time goes by much faster than you think.

domain: economics


text:
MY NAME IS ASHER LEV Celebrates 100 Performances at The Westside Theatre
My Name Is Asher Lev, the hit off-Broadway play based on the best selling novel by Chaim Potok, celebrated 100 performances at The Westside Theatre on Sunday, February 3rd. The show was recently extended and tickets are currently on sale through May 26, 2013.
"My producing partners and I are just thrilled that theatre audiences have embraced My Name is Asher Lev," said producer Darren Bagert. "After 100 performances and with our extension, it is wonderful that Chaim Potok's universal message of following your dreams will have the chance to touch even more people. We are truly humbled by the response."
Directed by Gordon Edelstein (The Glass Menagerie), the production at The Westside Theatre (upstairs) features Ari Brand (Broadway Bound, directed by David Cromer), who plays Asher Lev in this breakout performance. The New York Post praises him "With his dark eyes burning, Ari Brand brings a sincere intensity, in one of the most affecting performances of the year". Mark Nelson (nine Broadway productions including After The Fall, The Invention of Love, A Few Good Men; Obie Award winner for Picasso at the Lapin Agile) is "excellent" writes Associated Press and Jenny Bacon (Heartless, Omnium Gatherum) "appears lit from within" - Associated Press.
Set in post-war Brooklyn, My Name is Asher Lev tells the coming- of-age story of a boy prodigy who must be an artist at any cost -- against the will of family, community and tradition. It is the powerful telling of a young artist with a universal theme to follow our passion - against everything.
The distinguished creative team for My Name is Asher Lev includes set design by Eugene Lee (Wicked); lighting design by James F. Ingalls (Glengarry Glen Ross); costume design by Ilona Somogyi (Clybourne Park) and original music/sound design by John Gromada (Gore Vidal's The Best Man, Proof).

domain: entertainment


text:
When looking for a retirement destination, couples or individuals seek destinations that fit their lifestyle, hobbies and retirement ideals. There are a few things for pre-retirees to consider in [...]
Securities offered through Valmark Securities, INC., Member FINRA, SIPC. 130 Springside Dr. #300, Akron, OH 44333. 330576.1234. Cirrus Wealth Management and related entities are separate entities from ValMark Securities, Inc.

domain: finance


text:
Great news! All four of Elk River's outdoor ice rinks are open and ready to be enjoyed! Crews have been working overnight to flood the rinks in order to build the ice back up after recent unseasonably warm temperatures and rain.
Please remember, due to pandemic restrictions, the Handke warming house is not open, however there are additional benches and seating near the skating areas.
Rinks are located at the following locations in Elk River:
- Handke Pit (1133 4th Street NW)
- Spectrum Athletic Fields (10129 181st Avenue NW)
- Mississippi Oaks Park (17998 Elgin Street)
- Riverplace Park (13818 186th Lane)
For additional updates and information on the city’s outdoor skating rinks, please visit ElkRiverMN.gov/OutdoorSkating.

domain: geography