The Genome of the Amoeba Symbiont “Candidatus Amoebophilus asiaticus” Reveals Common Mechanisms for Host Cell Interaction among Amoeba-Associated Bacteria

Protozoa play host for many intracellular bacteria and are important for the adaptation of pathogenic bacteria to eukaryotic cells. We analyzed the genome sequence of “Candidatus Amoebophilus asiaticus,” an obligate intracellular amoeba symbiont belonging to the Bacteroidetes. The genome has a size of 1.89 Mbp, encodes 1,557 proteins, and shows massive proliferation of IS elements (24% of all genes), although the genome seems to be evolutionarily relatively stable. The genome does not encode pathways for de novo biosynthesis of cofactors, nucleotides, and almost all amino acids. “Ca. Amoebophilus asiaticus” encodes a variety of proteins with predicted importance for host cell interaction; in particular, an arsenal of proteins with eukaryotic domains, including ankyrin-, TPR/SEL1-, and leucine-rich repeats, which is hitherto unmatched among prokaryotes, is remarkable. Unexpectedly, 26 proteins that can interfere with the host ubiquitin system were identified in the genome. These proteins include F- and U-box domain proteins and two ubiquitin-specific proteases of the CA clan C19 family, representing the first prokaryotic members of this protein family. Consequently, interference with the host ubiquitin system is an important host cell interaction mechanism of “Ca. Amoebophilus asiaticus”. More generally, we show that the eukaryotic domains identified in “Ca. Amoebophilus asiaticus” are also significantly enriched in the genomes of other amoeba-associated bacteria (including chlamydiae, Legionella pneumophila, Rickettsia bellii, Francisella tularensis, and Mycobacterium avium). This indicates that phylogenetically and ecologically diverse bacteria which thrive inside amoebae exploit common mechanisms for interaction with their hosts, and it provides further evidence for the role of amoebae as training grounds for bacterial pathogens of humans.Free-living amoebae, such as Acanthamoeba spp., are ubiquitous protozoa which can be found in such diverse habitats as soil, marine water, and freshwater and in many engineered environments (62, 100). They are important predators of prokaryotic and eukaryotic microorganisms, thereby having great influence on microbial community composition, soil mineralization, plant growth, and nutrient cycles (14, 100). Interestingly, many well-known pathogens of humans are able to infect, survive, and multiply within amoebae (39, 51). These protozoa can thus serve as reservoirs and vectors for the transmission of pathogenic bacteria to humans, as demonstrated for L. pneumophila and Mycobacterium avium (2, 115). It is also increasingly being recognized that protozoa are important for the adaptation of (pathogenic) bacteria to higher eukaryotic cells as a niche for growth (2, 24, 42, 78, 89).In addition to the many recognized transient associations between amoeba and pathogens, stable and obligate relationships between bacteria and amoebae also were described for members of the Alphaproteobacteria (11, 34, 48), the Betaproteobacteria (49), the Bacteroidetes (50), and the Chlamydiae (4, 12, 35, 52). These obligate amoeba symbionts show a worldwide distribution, since phylogenetically highly similar strains were found in amoeba isolates from geographically distant sources (51, 107). The phylogenetic diversity and the different lifestyles of these obligate intracellular bacteria—some are located directly in the host cell cytoplasm (11, 34, 48-50, 52), while others are enclosed in host-derived vacuoles (4, 35, 44)—suggest fundamentally different mechanisms of host cell interaction. However, with the exception of chlamydia-related amoeba symbionts (37, 46, 47), our knowledge of the biology of obligate intracellular symbionts of amoebae is still scarce.Comparative genomics has been extremely helpful for the analysis of intracellular bacteria. Numerous genome sequences from the Alpha- and Gammaproteobacteria and the Chlamydiae are available and have contributed significantly to our understanding of genome evolution, the biology of intracellular bacteria, and the interactions with their host cells (24, 26, 46, 79, 82). In this study, we determined and analyzed the complete genome sequence of “Candidatus Amoebophilus asiaticus” strain 5a2 in order to gain novel insights into its biology. “Ca. Amoebophilus asiaticus” is a Gram-negative, obligate intracellular amoeba symbiont belonging to the Bacteroidetes which has been discovered within an amoeba isolated from lake sediment (107). “Ca. Amoebophilus asiaticus” shows highest 16S rRNA similarity to “Candidatus Cardinium hertigii,” an obligate intracellular parasite of arthropods able to manipulate the reproduction of its hosts (131). According to 16S rRNA trees, both organisms are members of a monophyletic group within the phylogenetically diverse phylum Bacteroidetes, consisting only of symbionts and sequences which were directly retrieved from corals (113). Among members of the Bacteroidetes, the genome sequences of only three symbionts, which are only distantly related (75 to 80% 16S rRNA sequence similarity) to “Ca. Amoebophilus asiaticus,” have been determined to date: two strains of “Candidatus Sulcia muelleri,” a symbiont of sharpshooters, and “Azobacteroides pseudotrichonymphae,” a symbiont of an anaerobic termite gut ciliate (45, 72, 74, 127).The genome of “Ca. Amoebophilus asiaticus” is only moderately reduced in size compared to those of many other obligate intracellular bacteria (75, 123), but nevertheless, its biosynthetic capabilities are extremely limited. A large fraction of the genome consists of IS elements and an unparalleled high number of proteins with eukaryotic domains, such as ankyrin repeats, TPR/SEL1 repeats, leucine-rich repeats, and domains from the eukaryotic ubiquitin system, all of them most likely important for host cell interaction. Feature enrichment analysis across a nonredundant data set of all bacterial genomes showed that these domains are enriched in the genomes of bacteria (including several pathogens of humans) known to be able to infect amoebae, providing further evidence for an important role of amoebae in the evolution of mechanisms for host cell interaction in intracellular bacteria.     

Widely Conserved Recombination Patterns among Single-Stranded DNA Viruses

The combinatorial nature of genetic recombination can potentially provide organisms with immediate access to many more positions in sequence space than can be reached by mutation alone. Recombination features particularly prominently in the evolution of a diverse range of viruses. Despite rapid progress having been made in the characterization of discrete recombination events for many species, little is currently known about either gross patterns of recombination across related virus families or the underlying processes that determine genome-wide recombination breakpoint distributions observable in nature. It has been hypothesized that the networks of coevolved molecular interactions that define the epistatic architectures of virus genomes might be damaged by recombination and therefore that selection strongly influences observable recombination patterns. For recombinants to thrive in nature, it is probably important that the portions of their genomes that they have inherited from different parents work well together. Here we describe a comparative analysis of recombination breakpoint distributions within the genomes of diverse single-stranded DNA (ssDNA) virus families. We show that whereas nonrandom breakpoint distributions in ssDNA virus genomes are partially attributable to mechanistic aspects of the recombination process, there is also a significant tendency for recombination breakpoints to fall either outside or on the peripheries of genes. In particular, we found significantly fewer recombination breakpoints within structural protein genes than within other gene types. Collectively, these results imply that natural selection acting against viruses expressing recombinant proteins is a major determinant of nonrandom recombination breakpoint distributions observable in most ssDNA virus families.Genetic recombination is a ubiquitous biological process that is both central to DNA repair pathways (10, 57) and an important evolutionary mechanism. By generating novel combinations of preexisting nucleotide polymorphisms, recombination can potentially accelerate evolution by increasing the population-wide genetic diversity upon which adaptive selection relies. Recombination can paradoxically also prevent the progressive accumulation of harmful mutations within individual genomes (18, 35, 53). Whereas its ability to defend high-fitness genomes from mutational decay possibly underlies the evolutionary value of sexuality in higher organisms, in many microbial species where pseudosexual genetic exchange is permissible among even highly divergent genomes, recombination can enable access to evolutionary innovations that would otherwise be inaccessible by mutation alone.Such interspecies recombination is fairly common in many virus families (8, 17, 27, 44, 82). It is becoming clear, however, that as with mutation events, most recombination events between distantly related genomes are maladaptive (5, 13, 38, 50, 63, 80). As genetic distances between parental genomes increase, so too does the probability of fitness defects in their recombinant offspring (16, 51). The viability of recombinants is apparently largely dependent on how severely recombination disrupts coevolved intragenome interaction networks (16, 32, 51). These networks include interacting nucleotide sequences that form secondary structures, sequence-specific protein-DNA interactions, interprotein interactions, and amino acid-amino acid interactions within protein three-dimensional folds.One virus family where such interaction networks appear to have a large impact on patterns of natural interspecies recombination are the single-stranded DNA (ssDNA) geminiviruses. As with other ssDNA viruses, recombination is very common among the species of this family (62, 84). Partially conserved recombination hot and cold spots have been detected in different genera (39, 81) and are apparently caused by both differential mechanistic predispositions of genome regions to recombination and natural selection disfavoring the survival of recombinants with disrupted intragenome interaction networks (38, 51).Genome organization and rolling circle replication (RCR)—the mechanism by which geminiviruses and many other ssDNA viruses replicate (9, 67, 79; see reference 24 for a review)—seem to have a large influence on basal recombination rates in different parts of geminivirus genomes (20, 33, 39, 61, 81). To initiate RCR, virion-strand ssDNA molecules are converted by host-mediated pathways into double-stranded “replicative-form” (RF) DNAs (34, 67). Initiated by a virus-encoded replication-associated protein (Rep) at a well-defined virion-strand replication origin (v-ori), new virion strands are synthesized on the complementary strand of RF DNAs (28, 73, 74) by host DNA polymerases. Virion-strand replication is concomitant with the displacement of old virion strands, which, once complete, yields covalently closed ssDNA molecules which are either encapsidated or converted into additional RF DNAs. Genome-wide basal recombination rates in ssDNA viruses are probably strongly influenced by the specific characteristics of host DNA polymerases that enable RCR. Interruption of RCR has been implicated directly in geminivirus recombination (40) and is most likely responsible for increased basal recombination rates both within genes transcribed in the opposite direction from that of virion-strand replication (40, 71) and at the v-ori (1, 9, 20, 69, 74).Whereas most ssDNA virus families replicate via either a rolling circle mechanism (the Nanoviridae, Microviridae, and Geminiviridae) (3, 23, 24, 31, 59, 67, 74) or a related rolling hairpin mechanism (the Parvoviridae) (25, 76), among the Circoviridae only the Circovirus genus is known to use RCR (45). Although the Gyrovirus genus (the other member of the Circoviridae) and the anelloviruses (a currently unclassified ssDNA virus group) might also use RCR, it is currently unknown whether they do or not (78). Additionally, some members of the Begomovirus genus of the Geminiviridae either have a second genome component, called DNA-B, or are associated with satellite ssDNA molecules called DNA-1 and DNA-Beta, all of which also replicate by RCR (1, 47, 68).Recombination is known to occur in the parvoviruses (19, 43, 70), microviruses (66), anelloviruses (40, 46), circoviruses (11, 26, 60), nanoviruses (30), geminivirus DNA-B components, and geminivirus satellite molecules (2, 62). Given that most, if not all, of these ssDNA replicons are evolutionarily related to and share many biological features with the geminiviruses (22, 31, 36), it is of interest to determine whether conserved recombination patterns observed in the geminiviruses (61, 81) are evident in these other groups. To date, no comparative analyses have ever been performed with different ssDNA virus families to identify, for example, possible influences of genome organization on recombination breakpoint distributions found in these viruses.Here we compare recombination frequencies and recombination breakpoint distributions in most currently described ssDNA viruses and satellite molecules and identify a number of sequence exchange patterns that are broadly conserved across this entire group.     

A New Borrelia Species Defined by Multilocus Sequence Analysis of Housekeeping Genes

Analysis of Lyme borreliosis (LB) spirochetes, using a novel multilocus sequence analysis scheme, revealed that OspA serotype 4 strains (a rodent-associated ecotype) of Borrelia garinii were sufficiently genetically distinct from bird-associated B. garinii strains to deserve species status. We suggest that OspA serotype 4 strains be raised to species status and named Borrelia bavariensis sp. nov. The rooted phylogenetic trees provide novel insights into the evolutionary history of LB spirochetes.Multilocus sequence typing (MLST) and multilocus sequence analysis (MLSA) have been shown to be powerful and pragmatic molecular methods for typing large numbers of microbial strains for population genetics studies, delineation of species, and assignment of strains to defined bacterial species (4, 13, 27, 40, 44). To date, MLST/MLSA schemes have been applied only to a few vector-borne microbial populations (1, 6, 30, 37, 40, 41, 47).Lyme borreliosis (LB) spirochetes comprise a diverse group of zoonotic bacteria which are transmitted among vertebrate hosts by ixodid (hard) ticks. The most common agents of human LB are Borrelia burgdorferi (sensu stricto), Borrelia afzelii, Borrelia garinii, Borrelia lusitaniae, and Borrelia spielmanii (7, 8, 12, 35). To date, 15 species have been named within the group of LB spirochetes (6, 31, 32, 37, 38, 41). While several of these LB species have been delineated using whole DNA-DNA hybridization (3, 20, 33), most ecological or epidemiological studies have been using single loci (5, 9-11, 29, 34, 36, 38, 42, 51, 53). Although some of these loci have been convenient for species assignment of strains or to address particular epidemiological questions, they may be unsuitable to resolve evolutionary relationships among LB species, because it is not possible to define any outgroup. For example, both the 5S-23S intergenic spacer (5S-23S IGS) and the gene encoding the outer surface protein A (ospA) are present only in LB spirochete genomes (36, 43). The advantage of using appropriate housekeeping genes of LB group spirochetes is that phylogenetic trees can be rooted with sequences of relapsing fever spirochetes. This renders the data amenable to detailed evolutionary studies of LB spirochetes.LB group spirochetes differ remarkably in their patterns and levels of host association, which are likely to affect their population structures (22, 24, 46, 48). Of the three main Eurasian Borrelia species, B. afzelii is adapted to rodents, whereas B. valaisiana and most strains of B. garinii are maintained by birds (12, 15, 16, 23, 26, 45). However, B. garinii OspA serotype 4 strains in Europe have been shown to be transmitted by rodents (17, 18) and, therefore, constitute a distinct ecotype within B. garinii. These strains have also been associated with high pathogenicity in humans, and their finer-scale geographical distribution seems highly focal (10, 34, 52, 53).In this study, we analyzed the intra- and interspecific phylogenetic relationships of B. burgdorferi, B. afzelii, B. garinii, B. valaisiana, B. lusitaniae, B. bissettii, and B. spielmanii by means of a novel MLSA scheme based on chromosomal housekeeping genes (30, 48).     

Coexistence of Different Genotypes in the Same Bat and Serological Characterization of Rousettus Bat Coronavirus HKU9 Belonging to a Novel Betacoronavirus Subgroup

Rousettus bat coronavirus HKU9 (Ro-BatCoV HKU9), a recently identified coronavirus of novel Betacoronavirus subgroup D, from Leschenault''s rousette, was previously found to display marked sequence polymorphism among genomes of four strains. Among 10 bats with complete RNA-dependent RNA polymerase (RdRp), spike (S), and nucleocapsid (N) genes sequenced, three and two sequence clades for all three genes were codetected in two and five bats, respectively, suggesting the coexistence of two or three distinct genotypes of Ro-BatCoV HKU9 in the same bat. Complete genome sequencing of the distinct genotypes from two bats, using degenerate/genome-specific primers with overlapping sequences confirmed by specific PCR, supported the coexistence of at least two distinct genomes in each bat. Recombination analysis using eight Ro-BatCoV HKU9 genomes showed possible recombination events between strains from different bat individuals, which may have allowed for the generation of different genotypes. Western blot assays using recombinant N proteins of Ro-BatCoV HKU9, Betacoronavirus subgroup A (HCoV-HKU1), subgroup B (SARSr-Rh-BatCoV), and subgroup C (Ty-BatCoV HKU4 and Pi-BatCoV HKU5) coronaviruses were subgroup specific, supporting their classification as separate subgroups under Betacoronavirus. Antibodies were detected in 75 (43%) of 175 and 224 (64%) of 350 tested serum samples from Leschenault''s rousette bats by Ro-BatCoV HKU9 N-protein-based Western blot and enzyme immunoassays, respectively. This is the first report describing coinfection of different coronavirus genotypes in bats and coronavirus genotypes of diverse nucleotide variation in the same host. Such unique phenomena, and the unusual instability of ORF7a, are likely due to recombination which may have been facilitated by the dense roosting behavior and long foraging range of Leschenault''s rousette.Coronaviruses infect a wide variety of animals in which they can cause respiratory, enteric, hepatic, and neurological diseases of various severities. Based on genotypic and serological characterization, coronaviruses were traditionally classified into three distinct groups, groups 1, 2, and 3 (3, 27, 59). Recently, the Coronavirus Study Group of the International Committee for Taxonomy of Viruses has renamed the traditional group 1, 2, and 3 coronaviruses as Alphacoronavirus, Betacoronavirus, and Gammacoronavirus, respectively (http://talk.ictvonline.org/media/p/1230.aspx). Coronaviruses are known to have a high frequency of recombination as a result of their unique mechanism of viral replication (27). Such tendency for recombination and high mutation rates may allow them to adapt to new hosts and ecological niches (24, 47, 52).The severe acute respiratory syndrome (SARS) epidemic has boosted interest in the study of coronaviruses in humans and animals (21, 34, 38, 41, 54). In the past few years, there has been a dramatic increase in the number of newly described human and animal coronaviruses (2, 4, 5, 8-10, 15-20, 23, 25, 28, 30, 32, 35, 36, 39, 43, 45, 50, 51, 53, 56, 58). Two novel human coronaviruses, human coronavirus NL63 (HCoV-NL63) and human coronavirus HKU1 (HCoV-HKU1), belonging to Alphacoronavirus and Betacoronavirus, respectively, have been discovered, in addition to the human coronavirus OC43 (HCoV-OC43), human coronavirus 229E (HCoV-229E), and SARS coronavirus (SARS-CoV) (17, 29, 45, 53, 55). We have also previously described the discovery of a diversity of novel coronaviruses in wild bats and birds in China, including SARSr-Rh-BatCoV, belonging to Betacoronavirus subgroup B, from Chinese horseshoe bats (30, 48, 56). Among these novel coronaviruses, three avian coronaviruses were found to belong to a novel subgroup of Gammacoronavirus (Gammacoronavirus subgroup C), while three bat coronaviruses were found to belong to two novel subgroups of Betacoronavirus (Betacoronavirus subgroups C and D) (48, 50). Based on the presence of the huge diversity of coronaviruses in Alphacoronavirus and Betacoronavirus among various bat species, bats are likely the reservoir for the ancestor of these two coronavirus genera (47).During our genome analysis of these novel coronaviruses, one of them, Rousettus bat coronavirus HKU9 (Ro-BatCoV HKU9), belonging to Betacoronavirus subgroup D, which was identified in Leschenault''s rousette bats, was found to display marked nucleotide and amino acid sequence polymorphism among the four strains with complete genome sequences (50). In our study on HCoV-HKU1, it has been shown that such sequence polymorphisms may indicate the presence of different genotypes (52). By complete genome sequence analysis of the potentially different genotypes of HCoV-HKU1, we have demonstrated for the first time natural recombination in a human coronavirus, resulting in the generation of at least three genotypes (52). We have also recently shown that recombination between different strains of SARSr-Rh-BatCoV from different regions of China may have given rise to the emergence of civet SARSr-CoV (31). To investigate the presence of different genotypes of Ro-BatCoV HKU9, the complete RNA-dependent RNA polymerase (RdRp) (corresponding to nsp12), spike (S), and nucleocapsid (N) gene sequences of Ro-BatCoV HKU9 from 10 additional bats were determined. Since sequence analysis showed the possible coexistence of different genotypes in seven bat individuals, complete genome sequencing of these distinct genotypes from two bats was carried out to investigate for possible recombination events among the different genotypes. In addition, serological characterization of Ro-BatCoV HKU9 was also performed by Western blot and enzyme immunoassays using recombinant Ro-BatCoV HKU9 nucleocapsid proteins and recombinant nucleocapsid proteins of Betacoronavirus subgroup A, B, and C coronaviruses to determine possible cross-reactivity among the different Betacoronavirus subgroups and the seroepidemiology of Ro-BatCoV HKU9 in Leschenault''s rousette bats.     

Optical Mapping Reveals a Large Genetic Inversion between Two Methicillin-Resistant Staphylococcus aureus Strains

Staphylococcus aureus is a highly versatile and evolving bacterium of great clinical importance. S. aureus can evolve by acquiring single nucleotide polymorphisms and mobile genetic elements and by recombination events. Identification and location of novel genomic elements in a bacterial genome are not straightforward, unless the whole genome is sequenced. Optical mapping is a new tool that creates a high-resolution, in situ ordered restriction map of a bacterial genome. These maps can be used to determine genomic organization and perform comparative genomics to identify genomic rearrangements, such as insertions, deletions, duplications, and inversions, compared to an in silico (virtual) restriction map of a known genome sequence. Using this technology, we report here the identification, approximate location, and characterization of a genetic inversion of ∼500 kb of a DNA element between the NRS387 (USA800) and FPR3757 (USA300) strains. The presence of the inversion and location of its junction sites were confirmed by site-specific PCR and sequencing. At both the left and right junction sites in NRS387, an IS1181 element and a 73-bp sequence were identified as inverted repeats, which could explain the possible mechanism of the inversion event.Staphylococcus aureus is a gram-positive bacterium of immense clinical importance. This opportunistic pathogen is capable of causing a wide range of diseases from skin and soft-tissue infections to life-threatening bacteremia, endocarditis, and osteomyelitis (14). Approximately 75% of the S. aureus genome is composed of a core genome that is common in all strains, and 25% of the genome is composed of variable regions which can differ between different strains (4, 16, 24-26). S. aureus evolves primarily by introducing single nucleotide polymorphisms in its core genome and by acquiring mobile genetic elements (MGEs) through horizontal gene transfer. These MGEs include pathogenicity/genomic islands, plasmids, transposons, and bacteriophages that become integrated in the chromosome (4, 11, 16, 31, 32). Despite being a heterogeneous organism, genetic recombination in S. aureus is proposed to be rather rare (20, 24, 29, 35). Its clones are more likely to evolve by point mutations than by recombination events (12). The MGEs contribute to the phenotypic and genotypic diversity seen with the S. aureus population. Acquisition of the staphylococcal cassette chromosome (SCCmec) elements through site-specific recombinases has led to the epidemic of methicillin-resistant S. aureus (MRSA) strains in hospitals and communities all over the world (6, 10, 15). In recent years, the integration of arginine catabolite mobile element in the USA300 lineage of MRSA has been proposed to give the pathogen its epidemiological advantage, including traits for surviving in low-pH conditions and oxygen tension environments (11). In addition, chromosomal replacements have been observed within lineages of sequence type 34 (ST34) and ST42 (34) and ST8 and ST30 (13).Genomic rearrangements, such as inversions, have been observed with genomes of enteric bacteria, such as Salmonella enterica, Shigella flexneri, Yersinia pestis KIM, Escherichia coli (K12 and O157:H7), and group A Streptococcus pyogenes (8, 9, 18, 27, 28, 30, 37). No genomic inversions in S. aureus have been reported to date. With the use of optical mapping, large genomic rearrangements, such as inversions, that would otherwise be missed with other comparative genotyping approaches, including microarray analysis, can be identified. Optical mapping uses high-resolution restriction maps (optical maps) of a bacterial genome to determine its genomic organization (5, 21, 23, 33, 36). These optical maps can be compared to an in silico (virtual) restriction map of a known genome sequence and can be used to identify gene rearrangements and their locations. Using optical mapping in conjunction with subsequent site-specific PCR and sequencing, we report the identification, approximate location, and partial characterization of an ∼500-kb DNA element in NRS387, a USA800 strain that was found to be inverted relative to USA300FPR3757. Identification of IS1181 elements and a novel 73-bp element at both ends of the ∼500-kb element in NRS387 could suggest the possibility of an inversion event in an ancestral strain of NRS387.     

Comparison of the Complete Genome Sequences of Bifidobacterium animalis subsp. lactis DSM 10140 and Bl-04

Bifidobacteria are important members of the human gut flora, especially in infants. Comparative genomic analysis of two Bifidobacterium animalis subsp. lactis strains revealed evolution by internal deletion of consecutive spacer-repeat units within a novel clustered regularly interspaced short palindromic repeat locus, which represented the largest differential content between the two genomes. Additionally, 47 single nucleotide polymorphisms were identified, consisting primarily of nonsynonymous mutations, indicating positive selection and/or recent divergence. A particular nonsynonymous mutation in a putative glucose transporter was linked to a negative phenotypic effect on the ability of the variant to catabolize glucose, consistent with a modification in the predicted protein transmembrane topology. Comparative genome sequence analysis of three Bifidobacterium species provided a core genome set of 1,117 orthologs complemented by a pan-genome of 2,445 genes. The genome sequences of the intestinal bacterium B. animalis subsp. lactis provide insights into rapid genome evolution and the genetic basis for adaptation to the human gut environment, notably with regard to catabolism of dietary carbohydrates, resistance to bile and acid, and interaction with the intestinal epithelium. The high degree of genome conservation observed between the two strains in terms of size, organization, and sequence is indicative of a genomically monomorphic subspecies and explains the inability to differentiate the strains by standard techniques such as pulsed-field gel electrophoresis.Actinobacteria, Firmicutes, Proteobacteria, and Bacteroidetes are dominant microbial phyla widely distributed in diverse ecosystems on the planet (10, 13, 20, 23, 33, 40, 51). Metagenomic analyses of the microbial landscape inhabiting various mammalian environments, notably the human gastrointestinal tract (GIT) and skin, have specifically identified Actinobacteria as an important and occasionally dominant phylum (18, 21, 33). Among the members of the large, diverse, and dynamic microbial community residing in the human GIT, Bifidobacterium is a dominant genus considered beneficial to humans and includes probiotic strains (live microorganisms which, when administered in adequate amounts, confer a health benefit on the host) (11). The population of bifidobacteria in the human intestine varies over time. Following vaginal delivery, the GIT of healthy newborns is typically colonized by bifidobacteria, especially in breast-fed infants, during the first few days of life (12). Interindividual variation, however, is remarkable in the human infant intestinal flora (41), and dominant genera are not always consistent across metagenomic analyses of the human gut flora (18, 30, 33, 41). Over time, the infant intestinal ecosystem becomes more complex as the diet becomes more diverse, with bifidobacteria typically remaining dominant until weaning (30).Bifidobacterium animalis subsp. lactis is a gram-positive lactic acid bacterium commonly found in the guts of healthy humans and has been identified in the infant gut biota, particularly in ileal, fecal, and mucosal samples (52, 56). Some strains of B. animalis subsp. lactis are able to survive in the GIT, to adhere to human epithelial cells in vitro, to modify fecal flora, to modulate the host immune response, or to prevent microbial gastroenteritis and colitis (4, 15, 20, 40, 52, 56). Additionally, B. animalis subsp. lactis has been reported to utilize nondigestible oligosaccharides, which may contribute to the organism''s ability to compete in the human gut. Carbohydrates resistant to enzymatic degradation and not absorbed in the upper intestinal tract are a primary source of energy for microbes residing in the large intestine. The benefits associated with probiotic strains of B. animalis subsp. lactis have resulted in their inclusion in the human diet via formulation into a large array of dietary supplements and foods, including dairy products such as yogurt. Deciphering the complete genome sequences of such microbes will provide additional insight into the genetic basis for survival and residence in the human gut, notably with regard to the ability to survive gastric passage and utilize available nutrients. Also, these genomes provide reference sequences for ongoing metagenomic analyses of the human environment, including the gut metagenome.Bifidobacterium animalis subsp. lactis is the most common bifidobacterium utilized as a probiotic in commercial dairy products in North America and Europe (22, 38). However, despite this commercial and probiotic significance, strain-level differentiation of B. animalis subsp. lactis strains has been hindered by the high genetic similarity of these organisms, as determined by pulsed-field gel electrophoresis and other nucleic acid-based techniques (6, 55, 56), and the lack of available genomic sequence information. The genome sequence of strain BB-12 (17) is not currently publicly available, and only a draft genome sequence in 28 contigs is available for strain HN019 (GenBank project 28807). The complete B. animalis subsp. lactis genome for strain AD011 (28) was only recently (2009) published. While this was an important first step, a single genome does not allow identification of unique targets for strain differentiation or comparative analyses within the subspecies.The objectives of this study were to determine the complete genome sequences of two B. animalis subsp. lactis strains, the type strain and a widely used commercial strain, to provide insights into the functionality of this species and into species identification and strain specialization.     

Host-Interactive Genes in Amerindian Helicobacter pylori Diverge from Their Old World Homologs and Mediate Inflammatory Responses

Helicobacter pylori is the dominant member of the gastric microbiota and has been associated with an increased risk of gastric cancer and peptic ulcers in adults. H. pylori populations have migrated and diverged with human populations, and health effects vary. Here, we describe the whole genome of the cag-positive strain V225d, cultured from a Venezuelan Piaroa Amerindian subject. To gain insight into the evolution and host adaptation of this bacterium, we undertook comparative H. pylori genomic analyses. A robust multiprotein phylogenetic tree reflects the major human migration out of Africa, across Europe, through Asia, and into the New World, placing Amerindian H. pylori as a particularly close sister group to East Asian H. pylori. In contrast, phylogenetic analysis of the host-interactive genes vacA and cagA shows substantial divergence of Amerindian from Old World forms and indicates new genotypes (e.g., VacA m3) involving these loci. Despite deletions in CagA EPIYA and CRPIA domains, V225d stimulates interleukin-8 secretion and the hummingbird phenotype in AGS cells. However, following a 33-week passage in the mouse stomach, these phenotypes were lost in isolate V225-RE, which had a 15-kb deletion in the cag pathogenicity island that truncated CagA and eliminated some of the type IV secretion system genes. Thus, the unusual V225d cag architecture was fully functional via conserved elements, but the natural deletion of 13 cag pathogenicity island genes and the truncation of CagA impaired the ability to induce inflammation.Helicobacter pylori is a microaerophilic bacterium of the Epsilonproteobacteria that has colonized the stomach since early in human evolution (45) and diverged with ancient human migrations (24, 45, 92). Thus, several major H. pylori populations, such as hpAfrica1, hpEurope, hspEAsia, and hspAmerind, whose names indicate their original geographic associations (45, 51), have been defined. In particular, similarities between the hspAmerind and hspEAsia populations suggest that the first colonizers of the New World brought H. pylori with them (24, 28). With recent mixing of human groups, H. pylori populations are also mixing and competing, with an apparent dominance by the hpEurope population at least in Latin America (19).H. pylori usually does not cause illness, but colonization with strains bearing the cag (cytotoxin-associated gene) pathogenicity island (cag PAI) (3, 7, 25, 52, 57, 61, 63) is associated with an increased risk of noncardia gastric adenocarcinoma and peptic ulcer disease (56, 64). Nonetheless, a high prevalence of cag-positive H. pylori strains occurs concurrently with low gastric cancer rates in Africa (40) and some regions in Latin America, such as the Venezuelan savannas and Amazonas (29, 53). Moreover, clinical and epidemiological data provide evidence for an inverse relationship between H. pylori colonization and the prevalence of certain metabolic disorders, esophageal diseases, asthma and allergic disorders, and acute infectious diseases, as well as a direct relationship with improved nutritional status of rural children (3, 14, 34, 37, 49, 68). That the host interaction with an indigenous gastric microbe provides some health benefits to the host is not unexpected given the well-established role of gastrointestinal microflora in maintaining gastroenteric homeostasis (8).The most thoroughly studied H. pylori proteins that interact with human cells are CagA and VacA. CagA is an effector protein injected into gastric epithelial cells by a type IV secretion system encoded by the cag PAI (10, 12, 15, 83). VacA is initially secreted from the bacterial cell by an autotransporter mechanism (16). Both proteins have multiple effects on host cells. Inside the host cell, phosphorylation of CagA on EPIYA repeats in the phosphotyrosine (PY) region (73) induces cellular elongation known as the hummingbird phenotype (72). CagA may also induce secretion of interleukin-8 (IL-8) (11), a process commonly attributed to NF-κB, and disrupt the barrier function of the tight junctions in polarized epithelial cells, leading to a loss of adhesion (1, 5). Other motifs in the PY region promote phosphorylation-independent effects (79). In addition, cagA may be considered an oncogene (60), since transgenic expression of cagA in mice leads to gastric epithelial hyperplasia through aberrant epithelial cell signaling and gastric carcinogenesis (60, 62). In contrast, VacA is a multifunctional protein with several activities in epithelial and immune cells (16). VacA induces cell vacuolation (43), alters mitochondrial membrane permeability (27, 41, 90), and increases epithelial monolayer permeability. VacA also activates several signal transduction pathways that are important in immune and epithelial cells, including the mitogen-activated protein (MAP) kinase and p38/ATF-2-mediated signal pathways (9, 55).Genomic analysis provides insights into the evolution of H. pylori strains and their relation with their human hosts and may be useful for the development of diagnostic tools and novel therapies. To date, there are six published complete H. pylori genomes, mostly from the hpEurope population (see Table SA1 in the supplemental material). Here, we report the whole genome of a newly characterized hspAmerind strain, V225d, and assess its genetic structure in comparison to those of Old World H. pylori strains through a comprehensive multiprotein phylogenetic analysis, as well as through single-gene examination of cagA and vacA, revealing clues to the evolution and migration of this strain into the New World and the implications for human health. We also present the results of functional and genomic studies using gastric epithelial cells demonstrating that V225d can induce an inflammatory host response, an effect that was lost following passage through the mouse stomach.     

Ecoepidemiology and Complete Genome Comparison of Different Strains of Severe Acute Respiratory Syndrome-Related Rhinolophus Bat Coronavirus in China Reveal Bats as a Reservoir for Acute,Self-Limiting Infection That Allows Recombination Events

Despite the identification of severe acute respiratory syndrome-related coronavirus (SARSr-CoV) in Rhinolophus Chinese horseshoe bats (SARSr-Rh-BatCoV) in China, the evolutionary and possible recombination origin of SARSr-CoV remains undetermined. We carried out the first study to investigate the migration pattern and SARSr-Rh-BatCoV genome epidemiology in Chinese horseshoe bats during a 4-year period. Of 1,401 Chinese horseshoe bats from Hong Kong and Guangdong, China, that were sampled, SARSr-Rh-BatCoV was detected in alimentary specimens from 130 (9.3%) bats, with peak activity during spring. A tagging exercise of 511 bats showed migration distances from 1.86 to 17 km. Bats carrying SARSr-Rh-BatCoV appeared healthy, with viral clearance occurring between 2 weeks and 4 months. However, lower body weights were observed in bats positive for SARSr-Rh-BatCoV, but not Rh-BatCoV HKU2. Complete genome sequencing of 10 SARSr-Rh-BatCoV strains showed frequent recombination between different strains. Moreover, recombination was detected between SARSr-Rh-BatCoV Rp3 from Guangxi, China, and Rf1 from Hubei, China, in the possible generation of civet SARSr-CoV SZ3, with a breakpoint at the nsp16/spike region. Molecular clock analysis showed that SARSr-CoVs were newly emerged viruses with the time of the most recent common ancestor (tMRCA) at 1972, which diverged between civet and bat strains in 1995. The present data suggest that SARSr-Rh-BatCoV causes acute, self-limiting infection in horseshoe bats, which serve as a reservoir for recombination between strains from different geographical locations within reachable foraging range. Civet SARSr-CoV is likely a recombinant virus arising from SARSr-CoV strains closely related to SARSr-Rh-BatCoV Rp3 and Rf1. Such frequent recombination, coupled with rapid evolution especially in ORF7b/ORF8 region, in these animals may have accounted for the cross-species transmission and emergence of SARS.Coronaviruses can infect a wide variety of animals, causing respiratory, enteric, hepatic, and neurological diseases with different degrees of severity. On the basis of genotypic and serological characteristics, coronaviruses were classified into three distinct groups (2, 20, 54). Among coronaviruses that infect humans, human coronavirus 229E (HCoV-229E) and human coronavirus NL63 (HCoV-NL63) belong to group 1 coronaviruses and human coronavirus OC43 (HCoV-OC43), and human coronavirus HKU1 (HCoV-HKU1) belong to group 2 coronaviruses, whereas severe acute respiratory syndrome-related coronavirus (SARSr-CoV) has been classified as a group 2b coronavirus, distantly related to group 2a, and the recently discovered group 2c and 2d coronaviruses (6, 8, 10, 18, 31, 38, 43, 46, 49, 50). Recently, the Coronavirus Study Group of the International Committee for Taxonomy of Viruses has proposed renaming the traditional group 1, 2, and 3 coronaviruses Alphacoronavirus, Betacoronavirus, and Gammacoronavirus, respectively (http://talk.ictvonline.org/media/p/1230.aspx).Among all coronaviruses, SARSr-CoV has caused the most severe disease in humans, with over 700 fatalities since the SARS epidemic in 2003. Although the identification of SARSr-CoV in Himalayan palm civets and raccoon dogs in live animal markets in southern China suggested that wild animals could be the origin of SARS (11), the presence of the virus in only market or farmed civets, but not civets in the wild, and the rapid evolution of SARSr-CoV genomes in market civets suggested that these caged animals were only intermediate hosts (24, 39, 42, 52). Since bats are commonly found and served in wild animal markets and restaurants in Guangdong, China (47), we have previously carried out a study of bats from the region and identified a SARSr-CoV in Rhinolophus Chinese horseshoe bats (SARSr-Rh-BatCoV) (21). Similar viruses have also been found in three other species of horseshoe bats in mainland China (25), supporting the hypothesis that horseshoe bats are a reservoir of SARSr-CoV. Recently, viruses closely related to SARSr-Rh-BatCoV in China were also reported in Chaerophon bats from Africa, although only partial RNA-dependent RNA polymerase (RdRp) sequences were available (41). In addition, more than 10 previously unrecognized coronaviruses of huge diversity have since been identified in bats from China and other countries (1, 3, 5, 9, 22, 27, 32, 33, 40, 46, 51), suggesting that bats play an important role in the ecology and evolution of coronaviruses.As a result of the unique mechanism of viral replication, coronaviruses have a high frequency of recombination (20). Such a high recombination rate, coupled with the infidelity of the polymerases of RNA viruses, may allow them to adapt to new hosts and ecological niches (12, 48). Recombination in coronaviruses was first recognized between different strains of murine hepatitis virus (MHV) and subsequently in other coronaviruses, such as infectious bronchitis virus, between MHV and bovine coronavirus, and between feline coronavirus type I and canine coronavirus generating feline coronavirus type II (12, 16, 17, 23). Recently, by complete genome analysis of 22 strains of HCoV-HKU1, we have also documented for the first time that natural recombination events in a human coronavirus can give rise to three different genotypes (48).Although previous studies have attempted to study the possible evolutionary and recombination origin of SARSr-CoV, no definite conclusion can be made on whether the viruses from bats are the direct ancestor of SARSr-CoV in civets and humans, given the paucity of available strains and genome sequences. To better define the epidemiology and evolution of SARSr-Rh-BatCoV in China and their role as a recombination origin of SARSr-CoV in civets, we carried out a 4-year study on coronaviruses in Chinese horseshoe bats in Hong Kong and Guangdong Province of southern China. Bat tagging was also performed to study the migration pattern of bats and viral persistence. The complete genomes of 10 strains of SARSr-Rh-BatCoV obtained at different time were sequenced and compared to previously sequenced genomes. With the availability of this larger set of genome sequences for more accurate analysis, recombination and molecular clock analyses were performed to elucidate the evolutionary origin and time of interspecies transmission of SARSr-CoV.     

Conserved Symbiotic Plasmid DNA Sequences in the Multireplicon Pangenomic Structure of Rhizobium etli

Strains of the same bacterial species often show considerable genomic variation. To examine the extent of such variation in Rhizobium etli, the complete genome sequence of R. etli CIAT652 and the partial genomic sequences of six additional R. etli strains having different geographical origins were determined. The sequences were compared with each other and with the previously reported genome sequence of R. etli CFN42. DNA sequences common to all strains constituted the greater part of these genomes and were localized in both the chromosome and large plasmids. About 700 to 1,000 kb of DNA that did not match sequences of the complete genomes of strains CIAT652 and CFN42 was unique to each R. etli strain. These sequences were distributed throughout the chromosome as individual genes or chromosomal islands and in plasmids, and they encoded accessory functions, such as transport of sugars and amino acids, or secondary metabolism; they also included mobile elements and hypothetical genes. Sequences corresponding to symbiotic plasmids showed high levels of nucleotide identity (about 98 to 99%), whereas chromosomal sequences and the sequences with matches to other plasmids showed lower levels of identity (on average, about 90 to 95%). We concluded that R. etli has a pangenomic structure with a core genome composed of both chromosomal and plasmid sequences, including a highly conserved symbiotic plasmid, despite the overall genomic divergence.It is becoming clear that bacterial genomes of strains of the same species vary widely both in size and in gene composition (39). An unexpected degree of genomic diversity has been found by comparing whole genomes (39). For instance, in Escherichia coli strains, differences of up to 1,400 kb account for some strain-specific pathogenic traits (5, 56). The extent of intraspecies genome diversity varies in different bacterial lineages. Some species have a wide range of variation; these species include E. coli (42), Streptococcus agalactiae (53), and Haloquadratum walsbyi (34). Other bacteria display only limited gene content diversity; an example is Ureaplasma urealyticum (1, 54). Tettelin and colleagues have suggested that bacterial species can be characterized by the presence of a pangenome consisting of a core genome containing genes present in all strains and a dispensable genome consisting of partially shared and strain-specific genes (53, 54). This concept is rooted in the earlier ideas of Reanney (43) and Campbell (7) concerning the structure of bacterial populations, and it indicates both that there is a pool of accessory genetic information in bacterial species and that strains of the same or even different species can obtain this information by horizontal transfer mechanisms (7, 43).Genome size and diversity are related to bacterial lifestyle. Small genomes are typical of strict pathogens such as Rickettsia prowazekii (2) and endosymbionts such as Buchnera aphidicola (44a). In contrast, free-living bacteria, such as Pseudomonas syringae and Streptomyces coelicolor, have large genomes (4, 6). The bacteria with the largest genomes are common inhabitants of heterogeneous environments, such as soil, where energy sources are limited but diverse (32). An increase in genome size is attributable mainly to expansion of functions such as secondary metabolism, transport of metabolites, and gene regulation. All these features are common to the nitrogen-fixing symbiotic bacteria of legumes, which are collectively known as rhizobia, and their close relative the plant pathogen Agrobacterium. The genomes of such bacterial species have diverse architectures with circular chromosomes that are different sizes or linear chromosomes, like that in Agrobacterium species, and the organisms contain variable numbers of large plasmids (31, 49). Comparative genomic studies have highlighted the conservation of gene content and order among the chromosomes of some species of rhizobia (22, 23, 25, 40). Furthermore, Guerrero and colleagues (25) observed that most essential genes occur in syntenic arrangements and display a higher level of sequence identity than nonsyntenic genes. In contrast, plasmids, including symbiotic plasmids and symbiotic chromosomal islands (like those in Mesorhizobium loti and Bradyrhizobium japonicum) are poorly conserved in terms of both gene content and gene order (21). It is not clear what evolutionary advantage, if any, is provided by multipartite genomes, but some authors have speculated that such genomes may allow further accumulation of genes independent of the chromosome. Recently, Slater and coworkers (46) proposed a model for the origin of secondary chromosomes. Their idea is based on the notion of intragenomic gene transfers that might occur from primary chromosomes to ancestral plasmids of the repABC type. Observations of conservation of clusters of genes in secondary chromosomes or in large plasmids that retain synteny with respect to the main chromosome support this hypothesis (46).We have been studying Rhizobium etli as a multipartite genome model species (23). This organism is a free-living soil bacterium that is able to form nodules and fix nitrogen in the roots of bean plants. The genome of R. etli is partitioned into several replicons, a circular chromosome, and several large plasmids. In the reference strain R. etli CFN42, the genome is composed of a circular chromosome consisting of about 4,381 kb and 6 large plasmids whose total size is 2,148 kb (23). A 371-kb plasmid, termed pSym or the symbiotic plasmid, contains most of the genes required for symbiosis (21). Previous studies have described the high level of genetic diversity among geographically different R. etli isolates (41). The strains are also variable with respect to the number and size of plasmids. Nevertheless, there has been no direct measurement of diversity at the genomic level, nor have comparative studies of shared and particular genomic features of R. etli strains been reported. Therefore, to assess the degrees of genomic difference and genomic similarity in R. etli, we obtained the complete genomic sequence of an additional R. etli strain and partial genomic sequences of six other R. etli strains isolated worldwide. Our results support the concept of a pangenomic structure at the multireplicon level and show that a highly conserved symbiotic plasmid is present in divergent R. etli isolates.     

Alignment of the c-Type Cytochrome OmcS along Pili of Geobacter sulfurreducens

Immunogold localization revealed that OmcS, a cytochrome that is required for Fe(III) oxide reduction by Geobacter sulfurreducens, was localized along the pili. The apparent spacing between OmcS molecules suggests that OmcS facilitates electron transfer from pili to Fe(III) oxides rather than promoting electron conduction along the length of the pili.There are multiple competing/complementary models for extracellular electron transfer in Fe(III)- and electrode-reducing microorganisms (8, 18, 20, 44). Which mechanisms prevail in different microorganisms or environmental conditions may greatly influence which microorganisms compete most successfully in sedimentary environments or on the surfaces of electrodes and can impact practical decisions on the best strategies to promote Fe(III) reduction for bioremediation applications (18, 19) or to enhance the power output of microbial fuel cells (18, 21).The three most commonly considered mechanisms for electron transfer to extracellular electron acceptors are (i) direct contact between redox-active proteins on the outer surfaces of the cells and the electron acceptor, (ii) electron transfer via soluble electron shuttling molecules, and (iii) the conduction of electrons along pili or other filamentous structures. Evidence for the first mechanism includes the necessity for direct cell-Fe(III) oxide contact in Geobacter species (34) and the finding that intensively studied Fe(III)- and electrode-reducing microorganisms, such as Geobacter sulfurreducens and Shewanella oneidensis MR-1, display redox-active proteins on their outer cell surfaces that could have access to extracellular electron acceptors (1, 2, 12, 15, 27, 28, 31-33). Deletion of the genes for these proteins often inhibits Fe(III) reduction (1, 4, 7, 15, 17, 28, 40) and electron transfer to electrodes (5, 7, 11, 33). In some instances, these proteins have been purified and shown to have the capacity to reduce Fe(III) and other potential electron acceptors in vitro (10, 13, 29, 38, 42, 43, 48, 49).Evidence for the second mechanism includes the ability of some microorganisms to reduce Fe(III) that they cannot directly contact, which can be associated with the accumulation of soluble substances that can promote electron shuttling (17, 22, 26, 35, 36, 47). In microbial fuel cell studies, an abundance of planktonic cells and/or the loss of current-producing capacity when the medium is replaced is consistent with the presence of an electron shuttle (3, 14, 26). Furthermore, a soluble electron shuttle is the most likely explanation for the electrochemical signatures of some microorganisms growing on an electrode surface (26, 46).Evidence for the third mechanism is more circumstantial (19). Filaments that have conductive properties have been identified in Shewanella (7) and Geobacter (41) species. To date, conductance has been measured only across the diameter of the filaments, not along the length. The evidence that the conductive filaments were involved in extracellular electron transfer in Shewanella was the finding that deletion of the genes for the c-type cytochromes OmcA and MtrC, which are necessary for extracellular electron transfer, resulted in nonconductive filaments, suggesting that the cytochromes were associated with the filaments (7). However, subsequent studies specifically designed to localize these cytochromes revealed that, although the cytochromes were extracellular, they were attached to the cells or in the exopolymeric matrix and not aligned along the pili (24, 25, 30, 40, 43). Subsequent reviews of electron transfer to Fe(III) in Shewanella oneidensis (44, 45) appear to have dropped the nanowire concept and focused on the first and second mechanisms.Geobacter sulfurreducens has a number of c-type cytochromes (15, 28) and multicopper proteins (12, 27) that have been demonstrated or proposed to be on the outer cell surface and are essential for extracellular electron transfer. Immunolocalization and proteolysis studies demonstrated that the cytochrome OmcB, which is essential for optimal Fe(III) reduction (15) and highly expressed during growth on electrodes (33), is embedded in the outer membrane (39), whereas the multicopper protein OmpB, which is also required for Fe(III) oxide reduction (27), is exposed on the outer cell surface (39).OmcS is one of the most abundant cytochromes that can readily be sheared from the outer surfaces of G. sulfurreducens cells (28). It is essential for the reduction of Fe(III) oxide (28) and for electron transfer to electrodes under some conditions (11). Therefore, the localization of this important protein was further investigated.