GenomeViz: visualizing microbial genomes

EDGAR3.0: comparative genomics and phylogenomics on a scalable infrastructure

The EDGAR platform, a web server providing databases of precomputed orthology data for thousands of microbial genomes, is one of the most established tools in the field of comparative genomics and phylogenomics. Based on precomputed gene alignments, EDGAR allows quick identification of the differential gene content, i.e. the pan genome, the core genome, or singleton genes. Furthermore, EDGAR features a wide range of analyses and visualizations like Venn diagrams, synteny plots, phylogenetic trees, as well as Amino Acid Identity (AAI) and Average Nucleotide Identity (ANI) matrices. During the last few years, the average number of genomes analyzed in an EDGAR project increased by two orders of magnitude. To handle this massive increase, a completely new technical backend infrastructure for the EDGAR platform was designed and launched as EDGAR3.0. For the calculation of new EDGAR3.0 projects, we are now using a scalable Kubernetes cluster running in a cloud environment. A new storage infrastructure was developed using a file-based high-performance storage backend which ensures timely data handling and efficient access. The new data backend guarantees a memory efficient calculation of orthologs, and parallelization has led to drastically reduced processing times. Based on the advanced technical infrastructure new analysis features could be implemented including POCP and FastANI genomes similarity indices, UpSet intersecting set visualization, and circular genome plots. Also the public database section of EDGAR was largely updated and now offers access to 24,317 genomes in 749 free-to-use projects. In summary, EDGAR 3.0 provides a new, scalable infrastructure for comprehensive microbial comparative gene content analysis. The web server is accessible at http://edgar3.computational.bio.

Discovery of in vivo Virulence Genes of Obligatory Intracellular Bacteria by Random Mutagenesis

Ehrlichia spp. are emerging tick-borne obligatory intracellular bacteria that cause febrile and sometimes fatal diseases with abnormal blood cell counts and signs of hepatitis. Ehrlichia HF strain provides an excellent mouse disease model of fatal human ehrlichiosis. We recently obtained and established stable culture of Ehrlichia HF strain in DH82 canine macrophage cell line, and obtained its whole genome sequence and annotation. To identify genes required for in vivo virulence of Ehrlichia, we constructed random insertional HF strain mutants by using Himar1 transposon-based mutagenesis procedure. Of total 158 insertional mutants isolated via antibiotic selection in DH82 cells, 74 insertions were in the coding regions of 55 distinct protein-coding genes, including TRP120 and multi-copy genes, such as p28/omp-1, virB2, and virB6. Among 84 insertions mapped within the non-coding regions, seven are located in the putative promoter region since they were within 50 bp upstream of the seven distinct genes. Using limited dilution methods, nine stable clonal mutants that had no apparent defect for multiplication in DH82 cells, were obtained. Mouse virulence of seven mutant clones was similar to that of wild-type HF strain, whereas two mutant clones showed significantly retarded growth in blood, livers, and spleens, and the mice inoculated with them lived longer than mice inoculated with wild-type. The two clones contained mutations in genes encoding a conserved hypothetical protein and a staphylococcal superantigen-like domain protein, respectively, and both genes are conserved among Ehrlichia spp., but lack homology to other bacterial genes. Inflammatory cytokine mRNA levels in the liver of mice infected with the two mutants were significantly diminished than those infected with HF strain wild-type, except IL-1β and IL-12 p40 in one clone. Thus, we identified two Ehrlichia virulence genes responsible for in vivo infection, but not for infection and growth in macrophages.

Comparative genomics analysis of Mycoplasma capricolum subsp. capripneumoniae 87001

Mycoplasma capricolum subsp. capripneumoniae (Mccp), belongs to Mycoplasma mycoides cluster and is a causal pathogen of contagious caprine pleuropneumonia (CCPP). This paper presents the complete annotated genome sequence of Mccp Strain 87001-a strain that was isolated from pneumonia affected goats on a farm in China, and comparative genomics analysis of five Mccp genomes in addition to comparative genomics within Mycoplasma mycoides cluster. The Mccp strain 87001 genome consists of a single circular chromosome 1017333 bp in length and encodes 898 open reading frames (orfs) averaging 944 bp in length. Fifty eight potential virulence genes were identified, including variable surface lipoproteins, hemolysin A, and P60 surface lipoprotein. Comparative genomic analysis revealed eight virulence genes and four extracellular genes which remained unchanged in five Mccp genomes for forty years, which can be used as potential target for drug development and vaccine design. We revealed 183 Mccp unique genes as markers to distinguish Mccp with other mycoplasma strains from goats, and different virulence factors contributing to host specificity and different syndrome of bovine pathogens and caprine pathogens.

Virulent Epidemic Pneumonia in Sheep Caused by the Human Pathogen Acinetobacter baumannii

The human pathogen Acinetobacter baumannii has emerged as a frequent cause of hospital-acquired infections, but infection of animals has rarely been observed. Here we analyzed an outbreak of epidemic pneumonia killing hundreds of sheep on a farm in Pakistan and identified A. baumannii as the infecting agent. A pure culture of strain AbPK1 isolated from lungs of sick animals was inoculated into healthy sheep, which subsequently developed similar disease symptoms. Bacteria re-isolated from the infected animals were shown to be identical to the inoculum, fulfilling Koch’s postulates. Comparison of the AbPK1 genome against 2283 A. baumannii genomes from the NCBI database revealed that AbPK1 carries genes for unusual surface structures, including a unique composition of iron acquisition genes, genes for O-antigen synthesis and sialic acid-specific acetylases of cell-surface carbohydrates that could enable immune evasion. Several of these unusual and otherwise rarely present genes were also identified in genomes of phylogenetically unrelated A. baumannii isolates from combat-wounded US military from Afghanistan indicating a common gene pool in this geographical region. Based on core genome MLST this virulent isolate represents a newly emerging lineage of Global Clone 2, suggesting a human source for this disease outbreak. The observed epidemic, direct transmission from sheep to sheep, which is highly unusual for A. baumannii, has important consequences for human and animal health. First, direct animal-to-animal transmission facilitates fast spread of pathogen and disease in the flock. Second, it may establish a stable ecological niche and subsequent spread in a new host. And third, it constitutes a serious risk of transmission of this hyper-virulent clone from sheep back to humans, which may result in emergence of contagious disease amongst humans.

Comparative Genomics of a Bovine Mycobacterium tuberculosis Isolate and Other Strains Reveals Its Potential Mechanism of Bovine Adaptation

The Mycobacterium tuberculosis complex causes tuberculosis (TB) in humans and other animal species, but Mycobacterium tuberculosis has a distinct host preference to humans. The present study aimed to determine whether a bovine M. tb strain 1458 has evolved some genetic properties in their genome that might be associated with their bovine adaptation. The genome of the M. tb strain 1458 was sequenced and subjected to an extensive comparative genomic analysis. A phylogenetic analysis showed that strain 1458 is most closely related to a Chinese M. tb strain, CCDC5079, of the same Beijing family. Compared with three human M. tb Beijing family strains, the strain 1458 has the fewest unique genes. However, there are most (21) IS6110 insertion sequences in the strain 1458 genome at either intragenic or intergenic sites, resulting in the interruption of 11 genes including three PPE family-encoding genes (PPE16, PPE38, and PPE59). Only the strain 1458 genome has the upstream insertion in esxS and phoP genes. PCR confirmed four upstream insertions and qPCR determined that transcription of esxS, phoP, dnaN, and ctpD genes differed significantly between M. tb strain 1458 and H37Rv or M. bovis. A Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis revealed that the genes affected by non-synonymous SNPs are enriched in RNA polymerase. Moreover, 127 of the 133 unique SNPs in strain 1458 are either different to those in the M. bovis genome. In conclusion, some critical genes responsible for bacterial virulence and immunogenicity were interrupted in the genome of bovine M. tb strain 1458 by IS insertions and non-synonymous SNPs, which might contribute to its bovine adaptation, and the modification of its virulence and immunogenicity in cattle.

Acquisition and loss of virulence-associated factors during genome evolution and speciation in three clades of Bordetella species

Background

The genus Bordetella consists of nine species that include important respiratory pathogens such as the ‘classical’ species B. bronchiseptica, B. pertussis and B. parapertussis and six more distantly related and less extensively studied species. Here we analyze sequence diversity and gene content of 128 genome sequences from all nine species with focus on the evolution of virulence-associated factors.

Results

Both genome-wide sequence-based and gene content-based phylogenetic trees divide the genus into three species clades. The phylogenies are congruent between species suggesting genus-wide co-evolution of sequence diversity and gene content, but less correlated within species, mainly because of strain-specific presence of many different prophages. We compared the genomes with focus on virulence-associated genes and identified multiple clade-specific, species-specific and strain-specific events of gene acquisition and gene loss, including genes encoding O-antigens, protein secretion systems and bacterial toxins. Gene loss was more frequent than gene gain throughout the evolution, and loss of hundreds of genes was associated with the origin of several species, including the recently evolved human-restricted B. pertussis and B. holmesii, B. parapertussis and the avian pathogen B. avium.

Conclusions

Acquisition and loss of multiple genes drive the evolution and speciation in the genus Bordetella, including large scale gene loss associated with the origin of several species. Recent loss and functional inactivation of genes, including those encoding pertussis vaccine components and bacterial toxins, in individual strains emphasize ongoing evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3112-5) contains supplementary material, which is available to authorized users.

A systematic study of the whole genome sequence of Amycolatopsis methanolica strain 239T provides an insight into its physiological and taxonomic properties which correlate with its position in the genus

The complete genome of methanol-utilizing Amycolatopsis methanolica strain 239^T was generated, revealing a single 7,237,391 nucleotide circular chromosome with 7074 annotated protein-coding sequences (CDSs). Comparative analyses against the complete genome sequences of Amycolatopsis japonica strain MG417-CF17^T, Amycolatopsis mediterranei strain U32 and Amycolatopsis orientalis strain HCCB10007 revealed a broad spectrum of genomic structures, including various genome sizes, core/quasi-core/non-core configurations and different kinds of episomes. Although polyketide synthase gene clusters were absent from the A. methanolica genome, 12 gene clusters related to the biosynthesis of other specialized (secondary) metabolites were identified. Complete pathways attributable to the facultative methylotrophic physiology of A. methanolica strain 239^T, including both the mdo/mscR encoded methanol oxidation and the hps/hpi encoded formaldehyde assimilation via the ribulose monophosphate cycle, were identified together with evidence that the latter might be the result of horizontal gene transfer. Phylogenetic analyses based on 16S rDNA or orthologues of AMETH_3452, a novel actinobacterial class-specific conserved gene against 62 or 18 Amycolatopsis type strains, respectively, revealed three major phyletic lineages, namely the mesophilic or moderately thermophilic A. orientalis subclade (AOS), the mesophilic Amycolatopsis taiwanensis subclade (ATS) and the thermophilic A. methanolica subclade (AMS). The distinct growth temperatures of members of the subclades correlated with corresponding genetic variations in their encoded compatible solutes. This study shows the value of integrating conventional taxonomic with whole genome sequence data.

Changes in pathogenicity and immunogenicity of Mycoplasma mycoides subsp. mycoides strains revealed by comparative genomics analysis

Mycoplasma mycoides subsp. mycoides is the causative agent of contagious bovine pleuropneumonia. A pathogenic strain BEN-1 was isolated from bovine lung and underwent continuous passages in rabbits for 468 generations. During this process, the strain's strong virulence became weak and, gradually, it lost the ability to confer protective immunity in cattle but developed virulence in rabbits. In order to gain insight into the mechanisms behind the reduction in virulence and the loss of immunogenicity, we sequenced five representative strains of the BEN series, including the original strain (BEN-1), the strain generation that first acquired virulence in rabbits (BEN-50), the two vaccine strain generations (BEN-181 and BEN-326), and the strain generation showing the greatest loss of immunogenicity (BEN-468). The gene mutation rate in the four different propagation stages varied greatly, and over half of variations observed in each generation were removed during the propagation process. However, the variation maintained in the BEN-468 generation might contribute to its changes in virulence and immunogenicity. We thus identified 18 genes associated with host adaptation, six genes contributing to virulence in cattle, and 35 genes participating in conferring immunity in cattle. These findings might help us optimize the vaccine to obtain more effective immunization results.

Complete genome sequences and comparative genome analysis of Lactobacillus plantarum strain 5-2 isolated from fermented soybean

Lactobacillus plantarum is an important probiotic and is mostly isolated from fermented foods. We sequenced the genome of L. plantarum strain 5-2, which was derived from fermented soybean isolated from Yunnan province, China. The strain was determined to contain 3114 genes. Fourteen complete insertion sequence (IS) elements were found in 5-2 chromosome. There were 24 DNA replication proteins and 76 DNA repair proteins in the 5-2 genome. Consistent with the classification of L. plantarum as a facultative heterofermentative lactobacillus, the 5-2 genome encodes key enzymes required for the EMP (Embden-Meyerhof-Parnas) and phosphoketolase (PK) pathways. Several components of the secretion machinery are found in the 5-2 genome, which was compared with L. plantarum ST-III, JDM1 and WCFS1. Most of the specific proteins in the four genomes appeared to be related to their prophage elements.

Genome analysis and in vivo virulence of porcine extraintestinal pathogenic Escherichia coli strain PCN033

Background

Strains of extraintestinal pathogenic Escherichia coli (ExPEC) can invade and colonize extraintestinal sites and cause a wide range of infections. Genomic analysis of ExPEC has mainly focused on isolates of human and avian origins, with porcine ExPEC isolates yet to be sequenced. To better understand the genomic attributes underlying the pathogenicity of porcine ExPEC, we isolated two E. coli strains PCN033 and PCN061 from pigs, assessed their in vivo virulence, and completed and compared their genomes.

Results

Animal experiments demonstrated that strain PCN033, but not PCN061, was pathogenic in a pig model. The chromosome of PCN033 was 384 kb larger than that of PCN061. Among the PCN033-specific sequences, genes encoding adhesins, unique lipopolysaccharide, unique capsular polysaccharide, iron acquisition and transport systems, and metabolism were identified. Additionally, a large plasmid PCN033p3 harboring many typical ExPEC virulence factors was identified in PCN033. Based on the genetic variation between PCN033 and PCN061, corresponding phenotypic differences in flagellum-dependent swarming motility and metabolism were verified. Furthermore, the comparative genomic analyses showed that the PCN033 genome shared many similarities with genomic sequences of human ExPEC strains. Additionally, comparison of PCN033 genome with other nine characteristic E. coli genomes revealed 425 PCN033-special coding sequences. Genes of this subset included those encoding type I restriction-modification (R-M) system, type VI secretion system (T6SS) and membrane-associated proteins.

Conclusions

The genetic and phenotypic differences between PCN033 and PCN061 could partially explain their differences in virulence, and also provide insight towards the molecular mechanisms of porcine ExPEC infections. Additionally, the similarities between the genomes of PCN033 and human ExPEC strains suggest that some connections between porcine and human ExPEC strains exist. The first completed genomic sequence for porcine ExPEC and the genomic differences identified by comparative analyses provide a baseline understanding of porcine ExPEC genetics and lay the foundation for their further study.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1890-9) contains supplementary material, which is available to authorized users.

A mutation burst during the acute phase of Helicobacter pylori infection in humans and rhesus macaques

The evolution rate and genetic changes that occur during chronic infection with Helicobacter pylori have been analysed, but little is known about the genomic changes during the initial, acute bacterial infection phase. Here we analyse the rate and pattern of genome evolution in H. pylori from the genomes of two input strains isolated from human volunteers with asymptomatic infection, and the genomes of two output strains collected 20 and 44 days after re-infection. Similarly, we analyse genome evolution in bacteria from the genome sequences of input and output strains sequentially taken after experimental infection of a rhesus macaque. The estimated mutation rate reveals a mutation burst during the acute infection phase that is over 10 times faster than the mutation rate during chronic infection, and orders of magnitude faster than mutation rates in any other bacteria. The elevated frequency of mutations in outer membrane protein genes suggests that the mutation burst facilitates rapid host adaptation of the bacteria.

Helicobacter pylori genomic microevolution during naturally occurring transmission between adults

The human gastric pathogen Helicobacter pylori is usually acquired during childhood and, in the absence of treatment, chronic infection persists through most of the host's life. However, the frequency and importance of H. pylori transmission between adults is underestimated. Here we sequenced the complete genomes of H. pylori strains that were transmitted between spouses and analysed the genomic changes. Similar to H. pylori from chronic infection, a significantly high proportion of the determined 31 SNPs and 10 recombinant DNA fragments affected genes of the hop family of outer membrane proteins, some of which are known to be adhesins. In addition, changes in a fucosyltransferase gene modified the LPS component of the bacterial cell surface, suggesting strong diversifying selection. In contrast, virulence factor genes were not affected by the genomic changes. We propose a model of the genomic changes that are associated with the transmission and adaptation of H. pylori to a new human host.

Comparative analysis of temperature-dependent transcriptome of Pseudomonas aeruginosa strains from rhizosphere and human habitats

In this study, we investigated the effects of a change in growth temperature on the transcriptome of two strains of Pseudomonas aeruginosa. The chosen P. aeruginosa strains were M18 and PAO1, which are adapted to two different niches, rhizosphere and human, respectively. To assess the changes induced by a change in temperature, we used a newly designed microarray covering the complete genome of four P. aeruginosa strains: PAO1, M18, PA14 and LESB58, which proved informative and reliable for the transcriptome study. Using the microarray, we analysed the transcriptome profile changes of two P. aeruginosa strains of M18 and PAO1 at their originating and non-originating temperatures: 28 °C for the rhizosphere and 37 °C for the human. The transcriptome profiles showed significant temperature-dependent differences (64.8 % in M18 and 66.8 % in PAO1) compared with the genome structure (6 % in M18 and 4.1 % in PAO1). Furthermore, we found that the specific induced genes at the non-originating growth temperature of the each strain (207 genes in M18 and 229 genes in PAO1) were evidently more than those induced at the originating growth temperature (158 genes in M18 and 169 genes in PAO1). The functional analysis of several newly found specific regulated operons (such as phh, liu, hmg) in the two strains indicated possible strategies implemented to respond to the non-originating temperature. This study provides new insight into how P. aeruginosa species responds to temperature change and a microarray platform covering the complete genomes of four widely studied P. aeruginosa strains.

Comparative genome analysis of the high pathogenicity Salmonella Typhimurium strain UK-1

Salmonella enterica serovar Typhimurium, a gram-negative facultative rod-shaped bacterium causing salmonellosis and foodborne disease, is one of the most common isolated Salmonella serovars in both developed and developing nations. Several S. Typhimurium genomes have been completed and many more genome-sequencing projects are underway. Comparative genome analysis of the multiple strains leads to a better understanding of the evolution of S. Typhimurium and its pathogenesis. S. Typhimurium strain UK-1 (belongs to phage type 1) is highly virulent when orally administered to mice and chickens and efficiently colonizes lymphoid tissues of these species. These characteristics make this strain a good choice for use in vaccine development. In fact, UK-1 has been used as the parent strain for a number of nonrecombinant and recombinant vaccine strains, including several commercial vaccines for poultry. In this study, we conducted a thorough comparative genome analysis of the UK-1 strain with other S. Typhimurium strains and examined the phenotypic impact of several genomic differences. Whole genomic comparison highlights an extremely close relationship between the UK-1 strain and other S. Typhimurium strains; however, many interesting genetic and genomic variations specific to UK-1 were explored. In particular, the deletion of a UK-1-specific gene that is highly similar to the gene encoding the T3SS effector protein NleC exhibited a significant decrease in oral virulence in BALB/c mice. The complete genetic complements in UK-1, especially those elements that contribute to virulence or aid in determining the diversity within bacterial species, provide key information in evaluating the functional characterization of important genetic determinants and for development of vaccines.

Genomic determinants of organohalide-respiration in Geobacter lovleyi, an unusual member of the Geobacteraceae

BACKGROUND

Geobacter lovleyi is a unique member of the Geobacteraceae because strains of this species share the ability to couple tetrachloroethene (PCE) reductive dechlorination to cis-1,2-dichloroethene (cis-DCE) with energy conservation and growth (i.e., organohalide respiration). Strain SZ also reduces U(VI) to U(IV) and contributes to uranium immobilization, making G. lovleyi relevant for bioremediation at sites impacted with chlorinated ethenes and radionuclides. G. lovleyi is the only fully sequenced representative of this distinct Geobacter clade, and comparative genome analyses identified genetic elements associated with organohalide respiration and elucidated genome features that distinguish strain SZ from other members of the Geobacteraceae.

RESULTS

Sequencing the G. lovleyi strain SZ genome revealed a 3.9 Mbp chromosome with 54.7% GC content (i.e., the percent of the total guanines (Gs) and cytosines (Cs) among the four bases within the genome), and average amino acid identities of 53-56% compared to other sequenced Geobacter spp. Sequencing also revealed the presence of a 77 kbp plasmid, pSZ77 (53.0% GC), with nearly half of its encoded genes corresponding to chromosomal homologs in other Geobacteraceae genomes. Among these chromosome-derived features, pSZ77 encodes 15 out of the 24 genes required for de novo cobalamin biosynthesis, a required cofactor for organohalide respiration. A plasmid with 99% sequence identity to pSZ77 was subsequently detected in the PCE-dechlorinating G. lovleyi strain KB-1 present in the PCE-to-ethene-dechlorinating consortium KB-1. Additional PCE-to-cis-DCE-dechlorinating G. lovleyi strains obtained from the PCE-contaminated Fort Lewis, WA, site did not carry a plasmid indicating that pSZ77 is not a requirement (marker) for PCE respiration within this species. Chromosomal genomic islands found within the G. lovleyi strain SZ genome encode two reductive dehalogenase (RDase) homologs and a putative conjugative pilus system. Despite the loss of many c-type cytochrome and oxidative-stress-responsive genes, strain SZ retained the majority of Geobacter core metabolic capabilities, including U(VI) respiration.

CONCLUSIONS

Gene acquisitions have expanded strain SZ's respiratory capabilities to include PCE and TCE as electron acceptors. Respiratory processes core to the Geobacter genus, such as metal reduction, were retained despite a substantially reduced number of c-type cytochrome genes. pSZ77 is stably maintained within its host strains SZ and KB-1, likely because the replicon carries essential genes including genes involved in cobalamin biosynthesis and possibly corrinoid transport. Lateral acquisition of the plasmid replicon and the RDase genomic island represent unique genome features of the PCE-respiring G. lovleyi strains SZ and KB-1, and at least the latter signifies adaptation to PCE contamination.

Edwardsiella comparative phylogenomics reveal the new intra/inter-species taxonomic relationships, virulence evolution and niche adaptation mechanisms

Edwardsiella bacteria are leading fish pathogens causing huge losses to aquaculture industries worldwide. E. tarda is a broad-host range pathogen that infects more than 20 species of fish and other animals including humans while E. ictaluri is host-adapted to channel catfish causing enteric septicemia of catfish (ESC). Thus, these two species consist of a useful comparative system for studying the intricacies of pathogen evolution. Here we present for the first time the phylogenomic comparisons of 8 genomes of E. tarda and E. ictaluri isolates. Genome-based phylogenetic analysis revealed that E. tarda could be separate into two kinds of genotypes (genotype I, EdwGI and genotype II, EdwGII) based on the sequence similarity. E. tarda strains of EdwGI were clustered together with the E. ictaluri lineage and showed low sequence conservation to E. tarda strains of EdwGII. Multilocus sequence analysis (MLSA) of 48 distinct Edwardsiella strains also supports the new taxonomic relationship of the lineages. We identified the type III and VI secretion systems (T3SS and T6SS) as well as iron scavenging related genes that fulfilled the criteria of a key evolutionary factor likely facilitating the virulence evolution and adaptation to a broad range of hosts in EdwGI E. tarda. The surface structure-related genes may underlie the adaptive evolution of E. ictaluri in the host specification processes. Virulence and competition assays of the null mutants of the representative genes experimentally confirmed their contributive roles in the evolution/niche adaptive processes. We also reconstructed the hypothetical evolutionary pathway to highlight the virulence evolution and niche adaptation mechanisms of Edwardsiella. This study may facilitate the development of diagnostics, vaccines, and therapeutics for this under-studied pathogen.

Comparative genomics and transcriptomics of lineages I, II, and III strains of Listeria monocytogenes

Background

Listeria monocytogenes is a food-borne pathogen that causes infections with a high-mortality rate and has served as an invaluable model for intracellular parasitism. Here, we report complete genome sequences for two L. monocytogenes strains belonging to serotype 4a (L99) and 4b (CLIP80459), and transcriptomes of representative strains from lineages I, II, and III, thereby permitting in-depth comparison of genome- and transcriptome -based data from three lineages of L. monocytogenes. Lineage III, represented by the 4a L99 genome is known to contain strains less virulent for humans.

Results

The genome analysis of the weakly pathogenic L99 serotype 4a provides extensive evidence of virulence gene decay, including loss of several important surface proteins. The 4b CLIP80459 genome, unlike the previously sequenced 4b F2365 genome harbours an intact inlB invasion gene. These lineage I strains are characterized by the lack of prophage genes, as they share only a single prophage locus with other L. monocytogenes genomes 1/2a EGD-e and 4a L99. Comparative transcriptome analysis during intracellular growth uncovered adaptive expression level differences in lineages I, II and III of Listeria, notable amongst which was a strong intracellular induction of flagellar genes in strain 4a L99 compared to the other lineages. Furthermore, extensive differences between strains are manifest at levels of metabolic flux control and phosphorylated sugar uptake. Intriguingly, prophage gene expression was found to be a hallmark of intracellular gene expression. Deletion mutants in the single shared prophage locus of lineage II strain EGD-e 1/2a, the lma operon, revealed severe attenuation of virulence in a murine infection model.

Conclusion

Comparative genomics and transcriptome analysis of L. monocytogenes strains from three lineages implicate prophage genes in intracellular adaptation and indicate that gene loss and decay may have led to the emergence of attenuated lineages.

Safety assessment of Bifidobacterium longum JDM301 based on complete genome sequences

AIM

To assess the safety of Bifidobacterium longum (B. longum) JDM301 based on complete genome sequences.

METHODS

The complete genome sequences of JDM301 were determined using the GS 20 system. Putative virulence factors, putative antibiotic resistance genes and genes encoding enzymes responsible for harmful metabolites were identified by blast with virulence factors database, antibiotic resistance genes database and genes associated with harmful metabolites in previous reports. Minimum inhibitory concentration of 16 common antimicrobial agents was evaluated by E-test.

RESULTS

JDM301 was shown to contain 36 genes associated with antibiotic resistance, 5 enzymes related to harmful metabolites and 162 nonspecific virulence factors mainly associated with transcriptional regulation, adhesion, sugar and amino acid transport. B. longum JDM301 was intrinsically resistant to ciprofloxacin, amikacin, gentamicin and streptomycin and susceptible to vancomycin, amoxicillin, cephalothin, chloramphenicol, erythromycin, ampicillin, cefotaxime, rifampicin, imipenem and trimethoprim-sulphamethoxazol. JDM301 was moderately resistant to bacitracin, while an earlier study showed that bifidobacteria were susceptible to this antibiotic. A tetracycline resistance gene with the risk of transfer was found in JDM301, which needs to be experimentally validated.

CONCLUSION

The safety assessment of JDM301 using information derived from complete bacterial genome will contribute to a wider and deeper insight into the safety of probiotic bacteria.

Phylogenomic evidence for the presence of a flagellum and cbb(3) oxidase in the free-living mitochondrial ancestor

The initiation of the intracellular symbiosis that would give rise to mitochondria and eukaryotes was a major event in the history of life on earth. Hypotheses to explain eukaryogenesis fall into two broad and competing categories: those proposing that the host was a phagocytotic proto-eukaryote that preyed upon the free-living mitochondrial ancestor (hereafter FMA), and those proposing that the host was an archaebacterium that engaged in syntrophy with the FMA. Of key importance to these hypotheses are whether the FMA was motile or nonmotile, and the atmospheric conditions under which the FMA thrived. Reconstructions of the FMA based on genome content of Rickettsiales representatives-generally considered to be the closest living relatives of mitochondria-indicate that it was nonmotile and aerobic. We have sequenced the genome of Candidatus Midichloria mitochondrii, a novel and phylogenetically divergent member of the Rickettsiales. We found that it possesses unique gene sets found in no other Rickettsiales, including 26 genes associated with flagellar assembly, and a cbb(3)-type cytochrome oxidase. Phylogenomic analyses show that these genes were inherited in a vertical fashion from an ancestral α-proteobacterium, and indicate that the FMA possessed a flagellum, and could undergo oxidative phosphorylation under both aerobic and microoxic conditions. These results indicate that the FMA played a more active and potentially parasitic role in eukaryogenesis than currently appreciated and provide an explanation for how the symbiosis could have evolved under low levels of oxygen.

Genomic analysis and temperature-dependent transcriptome profiles of the rhizosphere originating strain Pseudomonas aeruginosa M18

Background

Our previously published reports have described an effective biocontrol agent named Pseudomonas sp. M18 as its 16S rDNA sequence and several regulator genes share homologous sequences with those of P. aeruginosa, but there are several unusual phenotypic features. This study aims to explore its strain specific genomic features and gene expression patterns at different temperatures.

Results

The complete M18 genome is composed of a single chromosome of 6,327,754 base pairs containing 5684 open reading frames. Seven genomic islands, including two novel prophages and five specific non-phage islands were identified besides the conserved P. aeruginosa core genome. Each prophage contains a putative chitinase coding gene, and the prophage II contains a capB gene encoding a putative cold stress protein. The non-phage genomic islands contain genes responsible for pyoluteorin biosynthesis, environmental substance degradation and type I and III restriction-modification systems. Compared with other P. aeruginosa strains, the fewest number (3) of insertion sequences and the most number (3) of clustered regularly interspaced short palindromic repeats in M18 genome may contribute to the relative genome stability. Although the M18 genome is most closely related to that of P. aeruginosa strain LESB58, the strain M18 is more susceptible to several antimicrobial agents and easier to be erased in a mouse acute lung infection model than the strain LESB58. The whole M18 transcriptomic analysis indicated that 10.6% of the expressed genes are temperature-dependent, with 22 genes up-regulated at 28°C in three non-phage genomic islands and one prophage but none at 37°C.

Conclusions

The P. aeruginosa strain M18 has evolved its specific genomic structures and temperature dependent expression patterns to meet the requirement of its fitness and competitiveness under selective pressures imposed on the strain in rhizosphere niche.

Complete genome and proteome of Acholeplasma laidlawii

We present the complete genome sequence and proteogenomic map for Acholeplasma laidlawii PG-8A (class Mollicutes, order Acholeplasmatales, family Acholeplasmataceae). The genome of A. laidlawii is represented by a single 1,496,992-bp circular chromosome with an average G+C content of 31 mol%. This is the longest genome among the Mollicutes with a known nucleotide sequence. It contains genes of polymerase type I, SOS response, and signal transduction systems, as well as RNA regulatory elements, riboswitches, and T boxes. This demonstrates a significant capability for the regulation of gene expression and mutagenic response to stress. Acholeplasma laidlawii and phytoplasmas are the only Mollicutes known to use the universal genetic code, in which UGA is a stop codon. Within the Mollicutes group, only the sterol-nonrequiring Acholeplasma has the capacity to synthesize saturated fatty acids de novo. Proteomic data were used in the primary annotation of the genome, validating expression of many predicted proteins. We also detected posttranslational modifications of A. laidlawii proteins: phosphorylation and acylation. Seventy-four candidate phosphorylated proteins were found: 16 candidates are proteins unique to A. laidlawii, and 11 of them are surface-anchored or integral membrane proteins, which implies the presence of active signaling pathways. Among 20 acylated proteins, 14 contained palmitic chains, and six contained stearic chains. No residue of linoleic or oleic acid was observed. Acylated proteins were components of mainly sugar and inorganic ion transport systems and were surface-anchored proteins with unknown functions.

Transcriptomic response of Listeria monocytogenes during the transition to the long-term-survival phase

Listeria monocytogenes can change its cellular morphology from bacilli to cocci during the transition to the long-term-survival (LTS) phase. The LTS cells demonstrated increased baro- and thermotolerance compared to their vegetative counterparts. So far, the underlying mechanisms that trigger this morphological and physiological transition remain largely unknown. In this study, we compared the transcriptomic profiles of L. monocytogenes serotype 4b strain F2365 at different growth stages in tryptic soy broth with yeast extract (TSBYE) using a whole-genome DNA chip approach. We identified a total of 225 differentially expressed genes (≥4-fold; P < 0.05) during the transition to the LTS phase in TSBYE. Genes related to cell envelope structure, energy metabolism, and transport were most significantly upregulated in the LTS phase. The upregulation of compatible solute transporters may lead to the accumulation of cellular solutes, lowering intracellular water activity and thus increasing bacterial stress resistance during the transition to the LTS phase. The downregulation of genes associated with protein synthesis may indicate a status of metabolic dormancy of the LTS cells. The transcriptomic profiles of resuscitated LTS cells in fresh TSBYE resembled those of log-phase cells (r=0.94), as the LTS cells rapidly resume metabolic activities and transit back to log phase with decreased baro- and thermotolerance.

The complete genome sequence of Mycoplasma bovis strain Hubei-1

Infection by Mycoplasma bovis (M. bovis) can induce diseases, such as pneumonia and otitis media in young calves and mastitis and arthritis in older animals. Here, we report the finished and annotated genome sequence of M. bovis strain Hubei-1, a strain isolated in 2008 that caused calf pneumonia on a Chinese farm. The genome of M. bovis strain Hubei-1 contains a single circular chromosome of 953,114 bp with a 29.37% GC content. We identified 803 open reading frames (ORFs) that occupy 89.5% of the genome. While 34 ORFs were Hubei-1 specific, 662 ORFs had orthologs in the M. bovis type strain PG45 genome. Genome analysis validated lateral gene transfer between M. bovis and the Mycoplasma mycoides subspecies mycoides, while phylogenetic analysis found that the closest M. bovis neighbor is Mycoplasma agalactiae. Glycerol may be the main carbon and energy source of M. bovis, and most of the biosynthesis pathways were incomplete. We report that 47 lipoproteins, 12 extracellular proteins and 18 transmembrane proteins are phase-variable and may help M. bovis escape the immune response. Besides lipoproteins and phase-variable proteins, genomic analysis found two possible pathogenicity islands, which consist of four genes and 11 genes each, and several other virulence factors including hemolysin, lipoate protein ligase, dihydrolipoamide dehydrogenase, extracellular cysteine protease and 5′-nucleotidase.

CompaGB: An open framework for genome browsers comparison

Background

Tools to visualize and explore genomes hold a central place in genomics and the diversity of genome browsers has increased dramatically over the last few years. It often turns out to be a daunting task to compare and choose a well-adapted genome browser, as multidisciplinary knowledge is required to carry out this task and the number of tools, functionalities and features are overwhelming.

Findings

To assist in this task, we propose a community-based framework based on two cornerstones: (i) the implementation of industry promoted software qualification method (QSOS) adapted for genome browser evaluations, and (ii) a web resource providing numerous facilities either for visualizing comparisons or performing new evaluations. We formulated 60 criteria specifically for genome browsers, and incorporated another 65 directly from QSOS's generic section. Those criteria aim to answer versatile needs, ranging from a biologist whose interest primarily lies into user-friendly and informative functionalities, a bioinformatician who wants to integrate the genome browser into a wider framework, or a computer scientist who might choose a software according to more technical features. We developed a dedicated web application to enrich the existing QSOS functionalities (weighting of criteria, user profile) with features of interest to a community-based framework: easy management of evolving data, user comments...

Conclusions

The framework is available at http://genome.jouy.inra.fr/CompaGB. It is open to anyone who wishes to participate in the evaluations. It helps the scientific community to (1) choose a genome browser that would better fit their particular project, (2) visualize features comparatively with easily accessible formats, such as tables or radar plots and (3) perform their own evaluation against the defined criteria. To illustrate the CompaGB functionalities, we have evaluated seven genome browsers according to the implemented methodology. A summary of the features of the compared genome browsers is presented and discussed.

Highly specific and quick detection of Mycobacterium avium subsp. paratuberculosis in feces and gut tissue of cattle and humans by multiple real-time PCR assays

Mycobacterium avium subsp. paratuberculosis is the causative agent of Johne's disease (JD) in cattle and may be associated with Crohn's disease (CD) in humans. It is the slowest growing of the cultivable mycobacteria, and culture from clinical, veterinary, food, or environmental specimens can take 4 months or even longer. Currently, the insertion element IS900 is used to detect M. avium subsp. paratuberculosis DNA. However, closely related IS900 elements are also present in other mycobacteria, thus limiting its specificity as a target. Here we describe the use of novel primer sets derived from the sequences of two highly specific single copy genes, MAP2765c and MAP0865, for the quantitative detection of M. avium subsp. paratuberculosis within 6 h by using real-time PCR. Specificity of the target was established using 40 M. avium subsp. paratuberculosis isolates, 67 different bacterial species, and two intestinal parasites. Using the probes and methods described, we detected 27 (2.09%) M. avium subsp. paratuberculosis-positive stool specimens from 1,293 individual stool samples by the use of either IS900 or probes deriving from the MAP2765c and MAP0865 genes described here. In general, bacterial load due to M. avium subsp. paratuberculosis was uniformly low in these samples and we estimated 500 to 5,000 M. avium subsp. paratuberculosis bacteria per gram of stool in assay-positive samples. Thus, the methods described here are useful for rapid and specific detection of M. avium subsp. paratuberculosis in clinical samples.

Comparative genomic and transcriptomic analysis revealed genetic characteristics related to solvent formation and xylose utilization in Clostridium acetobutylicum EA 2018

Background

Clostridium acetobutylicum, a gram-positive and spore-forming anaerobe, is a major strain for the fermentative production of acetone, butanol and ethanol. But a previously isolated hyper-butanol producing strain C. acetobutylicum EA 2018 does not produce spores and has greater capability of solvent production, especially for butanol, than the type strain C. acetobutylicum ATCC 824.

Results

Complete genome of C. acetobutylicum EA 2018 was sequenced using Roche 454 pyrosequencing. Genomic comparison with ATCC 824 identified many variations which may contribute to the hyper-butanol producing characteristics in the EA 2018 strain, including a total of 46 deletion sites and 26 insertion sites. In addition, transcriptomic profiling of gene expression in EA 2018 relative to that of ATCC824 revealed expression-level changes of several key genes related to solvent formation. For example, spo0A and adhEII have higher expression level, and most of the acid formation related genes have lower expression level in EA 2018. Interestingly, the results also showed that the variation in CEA_G2622 (CAC2613 in ATCC 824), a putative transcriptional regulator involved in xylose utilization, might accelerate utilization of substrate xylose.

Conclusions

Comparative analysis of C. acetobutylicum hyper-butanol producing strain EA 2018 and type strain ATCC 824 at both genomic and transcriptomic levels, for the first time, provides molecular-level understanding of non-sporulation, higher solvent production and enhanced xylose utilization in the mutant EA 2018. The information could be valuable for further genetic modification of C. acetobutylicum for more effective butanol production.

The intracellular sRNA transcriptome of Listeria monocytogenes during growth in macrophages

Small non-coding RNAs (sRNAs) are widespread effectors of post-transcriptional gene regulation in bacteria. Currently extensive information exists on the sRNAs of Listeria monocytogenes expressed during growth in extracellular environments. We used deep sequencing of cDNAs obtained from fractioned RNA (<500 nt) isolated from extracellularly growing bacteria and from L. monocytogenes infected macrophages to catalog the sRNA repertoire during intracellular bacterial growth. Here, we report on the discovery of 150 putative regulatory RNAs of which 71 have not been previously described. A total of 29 regulatory RNAs, including small non-coding antisense RNAs, are specifically expressed intracellularly. We validated highly expressed sRNAs by northern blotting and demonstrated by the construction and characterization of isogenic mutants of rli31, rli33-1 and rli50* for intracellular expressed sRNA candidates, that their expression is required for efficient growth of bacteria in macrophages. All three mutants were attenuated when assessed for growth in mouse and insect models of infection. Comparative genomic analysis revealed the presence of lineage specific sRNA candidates and the absence of sRNA loci in genomes of naturally occurring infection-attenuated bacteria, with additional loss in non-pathogenic listerial genomes. Our analyses reveal extensive sRNA expression as an important feature of bacterial regulation during intracellular growth.

Blueprint for a minimal photoautotrophic cell: conserved and variable genes in Synechococcus elongatus PCC 7942

BACKGROUND

Simpler biological systems should be easier to understand and to engineer towards pre-defined goals. One way to achieve biological simplicity is through genome minimization. Here we looked for genomic islands in the fresh water cyanobacteria Synechococcus elongatus PCC 7942 (genome size 2.7 Mb) that could be used as targets for deletion. We also looked for conserved genes that might be essential for cell survival.

RESULTS

By using a combination of methods we identified 170 xenologs, 136 ORFans and 1401 core genes in the genome of S. elongatus PCC 7942. These represent 6.5%, 5.2% and 53.6% of the annotated genes respectively. We considered that genes in genomic islands could be found if they showed a combination of: a) unusual G+C content; b) unusual phylogenetic similarity; and/or c) a small number of the highly iterated palindrome 1 (HIP1) motif plus an unusual codon usage. The origin of the largest genomic island by horizontal gene transfer (HGT) could be corroborated by lack of coverage among metagenomic sequences from a fresh water microbialite. Evidence is also presented that xenologous genes tend to cluster in operons. Interestingly, most genes coding for proteins with a diguanylate cyclase domain are predicted to be xenologs, suggesting a role for horizontal gene transfer in the evolution of Synechococcus sensory systems.

CONCLUSIONS

Our estimates of genomic islands in PCC 7942 are larger than those predicted by other published methods like SIGI-HMM. Our results set a guide to non-essential genes in S. elongatus PCC 7942 indicating a path towards the engineering of a model photoautotrophic bacterial cell.

A linear molecule with two large inverted repeats: the mitochondrial genome of the stramenopile Proteromonas lacertae

Mitochondrial evolution has given rise to a complex array of organelles, ranging from classical aerobic mitochondria to mitochondrial remnants known as hydrogenosomes and mitosomes. The latter are found in anaerobic eukaryotes, and these highly derived organelles often retain only scant evidence of their mitochondrial origins. Intermediate evolutionary stages have also been reported as facultatively or even strictly anaerobic mitochondria, and hydrogenosomes that still retain some mitochondrial features. However, the diversity among these organelles with transitional features remains rather unclear and barely studied. Here, we report the sequence, structure, and gene content of the mitochondrial DNA of the anaerobic stramenopile Proteromonas lacertae. It has a linear genome with a unique central region flanked by two identical large inverted repeats containing numerous genes and “telomeres” with short inverted repeats. Comparison with the organelle genome of the strictly anaerobic human parasite Blastocystis reveals that, despite the close similarity of the sequences, features such as the genome structure display striking differences. It remains unclear whether the virtually identical gene repertoires are the result of convergence or descent.

ChIP on Chip: surprising results are often artifacts

Background

The method of chromatin immunoprecipitation combined with microarrays (ChIP-Chip) is a powerful tool for genome-wide analysis of protein binding. However, a high background signal is a common phenomenon.

Results

Reinvestigation of the chromatin immunoprecipitation procedure led us to discover four causes of high background: i) non-unique sequences, ii) incomplete reversion of crosslinks, iii) retention of protein in spin-columns and iv) insufficient RNase treatment. The chromatin immunoprecipitation method was modified and applied to analyze genome-wide binding of SeqA and σ³² in Escherichia coli.

Conclusions

False positive findings originating from these shortcomings of the method could explain surprising and contradictory findings in published ChIP-Chip studies. We present a modified chromatin immunoprecipitation method greatly reducing the background signal.

The complete genome sequence of Cupriavidus metallidurans strain CH34, a master survivalist in harsh and anthropogenic environments

Many bacteria in the environment have adapted to the presence of toxic heavy metals. Over the last 30 years, this heavy metal tolerance was the subject of extensive research. The bacterium Cupriavidus metallidurans strain CH34, originally isolated by us in 1976 from a metal processing factory, is considered a major model organism in this field because it withstands milli-molar range concentrations of over 20 different heavy metal ions. This tolerance is mostly achieved by rapid ion efflux but also by metal-complexation and -reduction. We present here the full genome sequence of strain CH34 and the manual annotation of all its genes. The genome of C. metallidurans CH34 is composed of two large circular chromosomes CHR1 and CHR2 of, respectively, 3,928,089 bp and 2,580,084 bp, and two megaplasmids pMOL28 and pMOL30 of, respectively, 171,459 bp and 233,720 bp in size. At least 25 loci for heavy-metal resistance (HMR) are distributed over the four replicons. Approximately 67% of the 6,717 coding sequences (CDSs) present in the CH34 genome could be assigned a putative function, and 9.1% (611 genes) appear to be unique to this strain. One out of five proteins is associated with either transport or transcription while the relay of environmental stimuli is governed by more than 600 signal transduction systems. The CH34 genome is most similar to the genomes of other Cupriavidus strains by correspondence between the respective CHR1 replicons but also displays similarity to the genomes of more distantly related species as a result of gene transfer and through the presence of large genomic islands. The presence of at least 57 IS elements and 19 transposons and the ability to take in and express foreign genes indicates a very dynamic and complex genome shaped by evolutionary forces. The genome data show that C. metallidurans CH34 is particularly well equipped to live in extreme conditions and anthropogenic environments that are rich in metals.

Identification of protective and broadly conserved vaccine antigens from the genome of extraintestinal pathogenic Escherichia coli

Extraintestinal pathogenic Escherichia coli (ExPEC) are a common cause of disease in both mammals and birds. A vaccine to prevent such infections would be desirable given the increasing antibiotic resistance of these bacteria. We have determined the genome sequence of ExPEC IHE3034 (ST95) isolated from a case of neonatal meningitis and compared this to available genome sequences of other ExPEC strains and a few nonpathogenic E. coli. We found 19 genomic islands present in the genome of IHE3034, which are absent in the nonpathogenic E. coli isolates. By using subtractive reverse vaccinology we identified 230 antigens present in ExPEC but absent (or present with low similarity) in nonpathogenic strains. Nine antigens were protective in a mouse challenge model. Some of them were also present in other pathogenic non-ExPEC strains, suggesting that a broadly protective E. coli vaccine may be possible. The gene encoding the most protective antigen was detected in most of the E. coli isolates, highly conserved in sequence and found to be exported by a type II secretion system which seems to be nonfunctional in nonpathogenic strains.

Legionella pneumophila pangenome reveals strain-specific virulence factors

BACKGROUND

Legionella pneumophila subsp. pneumophila is a gram-negative gamma-Proteobacterium and the causative agent of Legionnaires' disease, a form of epidemic pneumonia. It has a water-related life cycle. In industrialized cities L. pneumophila is commonly encountered in refrigeration towers and water pipes. Infection is always via infected aerosols to humans. Although many efforts have been made to eradicate Legionella from buildings, it still contaminates the water systems. The town of Alcoy (Valencian Region, Spain) has had recurrent outbreaks since 1999. The strain "Alcoy 2300/99" is a particularly persistent and recurrent strain that was isolated during one of the most significant outbreaks between the years 1999-2000.

RESULTS

We have sequenced the genome of the particularly persistent L. pneumophila strain Alcoy 2300/99 and have compared it with four previously sequenced strains known as Philadelphia (USA), Lens (France), Paris (France) and Corby (England).Pangenome analysis facilitated the identification of strain-specific features, as well as some that are shared by two or more strains. We identified: (1) three islands related to anti-drug resistance systems; (2) a system for transport and secretion of heavy metals; (3) three systems related to DNA transfer; (4) two CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) systems, known to provide resistance against phage infections, one similar in the Lens and Alcoy strains, and another specific to the Paris strain; and (5) seven islands of phage-related proteins, five of which seem to be strain-specific and two shared.

CONCLUSIONS

The dispensable genome disclosed by the pangenomic analysis seems to be a reservoir of new traits that have mainly been acquired by horizontal gene transfer and could confer evolutionary advantages over strains lacking them.

3D genome tuner: compare multiple circular genomes in a 3D context

Circular genomes, being the largest proportion of sequenced genomes, play an important role in genome analysis. However, traditional 2D circular map only provides an overview and annotations of genome but does not offer feature-based comparison. For remedying these shortcomings, we developed 3D Genome Tuner, a hybrid of circular map and comparative map tools. Its capability of viewing comparisons between multiple circular maps in a 3D space offers great benefits to the study of comparative genomics. The program is freely available (under an LGPL licence) at http://sourceforge.net/projects/dgenometuner.

Identification of potential therapeutic targets for Burkholderia cenocepacia by comparative transcriptomics

Background

Burkholderia cenocepacia is an endemic soil dweller and emerging opportunistic pathogen in patients with cystic fibrosis (CF). The identification of virulence factors and potential therapeutic targets has been hampered by the genomic diversity within the species as many factors are not shared among the pathogenic members of the species.

Methodology/Principal Findings

In this study, global identification of putative virulence factors was performed by analyzing the transcriptome of two related strains of B. cenocepacia (one clinical, one environmental) under conditions mimicking cystic fibrosis sputum versus soil. Soil is a natural reservoir for this species; hence, genes induced under CF conditions relative to soil may represent adaptations that have occurred in clinical strains. Under CF conditions, several genes encoding proteins thought to be involved in virulence were induced and many new ones were identified. Our analysis, in combination with previous studies, reveals 458 strain-specific genes, 126 clinical-isolate-specific, and at least four species-specific genes that are induced under CF conditions. The chromosomal distribution of the induced genes was disproportionate to the size of the chromosome as genes expressed under soil conditions by both strains were more frequent on the second chromosome and those differentially regulated between strains were more frequent on the third chromosome. Conservation of these induced genes was established using the 11 available Bcc genome sequences to indicate whether potential therapeutic targets would be species-wide.

Conclusions/Significance

Comparative transcriptomics is a useful way to identify new potential virulence factors and therapeutic targets for pathogenic bacteria. We identified eight genes induced under CF conditions that were also conserved in the Bcc and may constitute particularly attractive therapeutic targets due to their signal sequence, predicted cellular location, and homology to known therapeutic targets.

Genomic insights into the convergence and pathogenicity factors of Campylobacter jejuni and Campylobacter coli species

Whether or not bacteria form coherent evolutionary groups via means of genetic exchange and, hence, elicit distinct species boundaries remains an unsettled issue. A recent report implied that not only may the former be true but also, in fact, the clearly distinct Campylobacter jejuni and Campylobacter coli species may be converging as a consequence of increased interspecies gene flow fostered, presumably, by the recent invasion of an overlapping ecological niche (S. K. Sheppard, N. D. McCarthy, D. Falush, and M. C. Maiden, Science 320:237-239, 2008). We have reanalyzed the Campylobacter multilocus sequence typing database used in the previous study and found that the number of interspecies gene transfer events may actually be too infrequent to account, unequivocally, for species convergence. For instance, only 1 to 2% of the 4,507 Campylobacter isolates examined appeared to have imported gene alleles from another Campylobacter species. Furthermore, by analyzing the available Campylobacter genomic sequences, we show that although there seems to be a slightly higher number of exchanged genes between C. jejuni and C. coli relative to other comparable species ( approximately 10% versus 2 to 3% of the total genes in the genome, respectively), the function and spatial distribution in the genome of the exchanged genes are far from random, and hence, inconsistent with the species convergence hypothesis. In fact, the exchanged genes appear to be limited to a few environmentally selected cellular functions. Accordingly, these genes may represent important pathogenic determinants of pathogenic Campylobacter, and convergence of (any) two bacterial species remains to be seen.

Mapping the Burkholderia cenocepacia niche response via high-throughput sequencing

Determining how an organism responds to its environment by altering gene expression is key to understanding its ecology. Here, we used RNA-seq to comprehensively and quantitatively assess the transcriptional response of the bacterial opportunistic cystic fibrosis (CF) pathogen and endemic soil dweller, Burkholderia cenocepacia, in conditions mimicking these 2 environments. By sequencing 762 million bases of cDNA from 2 closely related B. cenocepacia strains (one isolated from a CF patient and one from soil), we identified a number of potential virulence factors expressed under CF-like conditions, whereas genes whose protein products are involved in nitrogen scavenging and 2-component sensing were among those induced under soil-like conditions. Interestingly, 13 new putative noncoding RNAs were discovered using this technique, 12 of which are preferentially induced in the soil environment, suggesting that ncRNAs play an important role in survival in the soil. In addition, we detected a surprisingly large number of regulatory differences between the 2 strains, which may represent specific adaptations to the niches from which each strain was isolated, despite their high degree of DNA sequence similarity. Compared with the CF strain, the soil strain shows a stronger global gene expression response to its environment, which is consistent with the need for a more dynamic reaction to the heterogeneous conditions of soil.

The virophage as a unique parasite of the giant mimivirus

Viruses are obligate parasites of Eukarya, Archaea and Bacteria. Acanthamoeba polyphaga mimivirus (APMV) is the largest known virus; it grows only in amoeba and is visible under the optical microscope. Mimivirus possesses a 1,185-kilobase double-stranded linear chromosome whose coding capacity is greater than that of numerous bacteria and archaea1, 2, 3. Here we describe an icosahedral small virus, Sputnik, 50 nm in size, found associated with a new strain of APMV. Sputnik cannot multiply in Acanthamoeba castellanii but grows rapidly, after an eclipse phase, in the giant virus factory found in amoebae co-infected with APMV4. Sputnik growth is deleterious to APMV and results in the production of abortive forms and abnormal capsid assembly of the host virus. The Sputnik genome is an 18.343-kilobase circular double-stranded DNA and contains genes that are linked to viruses infecting each of the three domains of life Eukarya, Archaea and Bacteria. Of the 21 predicted protein-coding genes, eight encode proteins with detectable homologues, including three proteins apparently derived from APMV, a homologue of an archaeal virus integrase, a predicted primase-helicase, a packaging ATPase with homologues in bacteriophages and eukaryotic viruses, a distant homologue of bacterial insertion sequence transposase DNA-binding subunit, and a Zn-ribbon protein. The closest homologues of the last four of these proteins were detected in the Global Ocean Survey environmental data set5, suggesting that Sputnik represents a currently unknown family of viruses. Considering its functional analogy with bacteriophages, we classify this virus as a virophage. The virophage could be a vehicle mediating lateral gene transfer between giant viruses.

Analysis of two genomes from the mitochondrion-like organelle of the intestinal parasite Blastocystis: complete sequences, gene content, and genome organization

Acquisition of mitochondria by the ancestor of all living eukaryotes represented a crucial milestone in the evolution of the eukaryotic cell. Nevertheless, a number of anaerobic unicellular eukaryotes have secondarily discarded certain mitochondrial features, leading to modified organelles such as hydrogenosomes and mitosomes via degenerative evolution. These mitochondrion-derived organelles have lost many of the typical characteristics of aerobic mitochondria, including certain metabolic pathways, morphological traits, and, in most cases, the organellar genome. So far, the evolutionary pathway leading from aerobic mitochondria to anaerobic degenerate organelles has remained unclear due to the lack of examples representing intermediate stages. The human parasitic stramenopile Blastocystis is a rare example of an anaerobic eukaryote with organelles that have retained some mitochondrial characteristics, including a genome, whereas they lack others, such as cytochromes. Here we report the sequence and comparative analysis of the organellar genome from two different Blastocystis isolates as well as a comparison to other genomes from stramenopile mitochondria. Analysis of the characteristics displayed by the unique Blastocystis organelle genome gives us an insight into the initial evolutionary steps that may have led from mitochondria to hydrogenosomes and mitosomes.

Genome sequence of a Lancefield group C Streptococcus zooepidemicus strain causing epidemic nephritis: new information about an old disease

Outbreaks of disease attributable to human error or natural causes can provide unique opportunities to gain new information about host-pathogen interactions and new leads for pathogenesis research. Poststreptococcal glomerulonephritis (PSGN), a sequela of infection with pathogenic streptococci, is a common cause of preventable kidney disease worldwide. Although PSGN usually occurs after infection with group A streptococci, organisms of Lancefield group C and G also can be responsible. Despite decades of study, the molecular pathogenesis of PSGN is poorly understood. As a first step toward gaining new information about PSGN pathogenesis, we sequenced the genome of Streptococcus equi subsp. zooepidemicus strain MGCS10565, a group C organism that caused a very large and unusually severe epidemic of nephritis in Brazil. The genome is a circular chromosome of 2,024,171 bp. The genome shares extensive gene content, including many virulence factors, with genetically related group A streptococci, but unexpectedly lacks prophages. The genome contains many apparently foreign genes interspersed around the chromosome, consistent with the presence of a full array of genes required for natural competence. An inordinately large family of genes encodes secreted extracellular collagen-like proteins with multiple integrin-binding motifs. The absence of a gene related to speB rules out the long-held belief that streptococcal pyrogenic exotoxin B or antibodies reacting with it singularly cause PSGN. Many proteins previously implicated in GAS PSGN, such as streptokinase, are either highly divergent in strain MGCS10565 or are not more closely related between these species than to orthologs present in other streptococci that do not commonly cause PSGN. Our analysis provides a comparative genomics framework for renewed appraisal of molecular events underlying APSGN pathogenesis.

Genetic basis of virulence attenuation revealed by comparative genomic analysis of Mycobacterium tuberculosis strain H37Ra versus H37Rv

Tuberculosis, caused by Mycobacterium tuberculosis, remains a leading infectious disease despite the availability of chemotherapy and BCG vaccine. The commonly used avirulent M. tuberculosis strain H37Ra was derived from virulent strain H37 in 1935 but the basis of virulence attenuation has remained obscure despite numerous studies. We determined the complete genomic sequence of H37Ra ATCC25177 and compared that with its virulent counterpart H37Rv and a clinical isolate CDC1551. The H37Ra genome is highly similar to that of H37Rv with respect to gene content and order but is 8,445 bp larger as a result of 53 insertions and 21 deletions in H37Ra relative to H37Rv. Variations in repetitive sequences such as IS6110 and PE/PPE/PE-PGRS family genes are responsible for most of the gross genetic changes. A total of 198 single nucleotide variations (SNVs) that are different between H37Ra and H37Rv were identified, yet 119 of them are identical between H37Ra and CDC1551 and 3 are due to H37Rv strain variation, leaving only 76 H37Ra-specific SNVs that affect only 32 genes. The biological impact of missense mutations in protein coding sequences was analyzed in silico while nucleotide variations in potential promoter regions of several important genes were verified by quantitative RT-PCR. Mutations affecting transcription factors and/or global metabolic regulations related to in vitro survival under aging stress, and mutations affecting cell envelope, primary metabolism, in vivo growth as well as variations in the PE/PPE/PE-PGRS family genes, may underlie the basis of virulence attenuation. These findings have implications not only for improved understanding of pathogenesis of M. tuberculosis but also for development of new vaccines and new therapeutic agents.

Comparative microbial genome visualization using GenomeViz

Recent years have brought a tremendous increase in the amount of sequence data from various bacterial genome sequencing projects, an increase that is projected to accelerate over the next years. Comparative genomics of microbial strains has provided us with unprecedented information to describe a bacterial species and examine for microbial diversity. This has allowed us to define core genomes based on genes commonly present in all strains of a species or genus and to identify dispensable regions in the genome harboring genus-, species-, and even strain-specific genes. Nevertheless, the task of organizing and summarizing the data to extract the most informative features remains a challenging yet critical endeavor. Visualization is an effective way of structuring and presenting such information effectively, in a concise and eloquent fashion. The large-scale views unveil commonalities and differences between the genomes that may shed light on their evolutionary relationships and define characteristics that are typical of pathogenicity or other ecological adaptations. We describe GenomeViz, a tool for comparative visualization of bacterial genomes that allows the user to actively create, modify and query a genome plot in a visually compact, user-friendly, and interactive manner.

Variation in virulence among clades of Escherichia coli O157:H7 associated with disease outbreaks

Escherichia coli O157:H7, a toxin-producing food and waterborne bacterial pathogen, has been linked to large outbreaks of gastrointestinal illness for more than two decades. E. coli O157 causes a wide range of clinical illness that varies by outbreak, although factors that contribute to variation in disease severity are poorly understood. Several recent outbreaks involving O157 contamination of fresh produce (e.g., spinach) were associated with more severe disease, as defined by higher hemolytic uremic syndrome and hospitalization frequencies, suggesting that increased virulence has evolved. To test this hypothesis, we developed a system that detects SNPs in 96 loci and applied it to >500 E. coli O157 clinical strains. Phylogenetic analyses identified 39 SNP genotypes that differ at 20% of SNP loci and are separated into nine distinct clades. Differences were observed between clades in the frequency and distribution of Shiga toxin genes and in the type of clinical disease reported. Patients with hemolytic uremic syndrome were significantly more likely to be infected with clade 8 strains, which have increased in frequency over the past 5 years. Genome sequencing of a spinach outbreak strain, a member of clade 8, also revealed substantial genomic differences. These findings suggest that an emergent subpopulation of the clade 8 lineage has acquired critical factors that contribute to more severe disease. The ability to detect and rapidly genotype O157 strains belonging to such lineages is important and will have a significant impact on both disease diagnosis and treatment guidelines.

Temporal transcriptomic analysis of the Listeria monocytogenes EGD-e sigmaB regulon

Background

The opportunistic food-borne gram-positive pathogen Listeria monocytogenes can exist as a free-living microorganism in the environment and grow in the cytoplasm of vertebrate and invertebrate cells following infection. The general stress response, controlled by the alternative sigma factor, σ^B, has an important role for bacterial survival both in the environment and during infection. We used quantitative real-time PCR analysis and immuno-blot analysis to examine σ^B expression during growth of L. monocytogenes EGD-e. Whole genome-based transcriptional profiling was used to identify σ^B-dependent genes at different growth phases.

Results

We detected 105 σ^B-positively regulated genes and 111 genes which appeared to be under negative control of σ^B and validated 36 σ^B-positively regulated genes in vivo using a reporter gene fusion system.

Conclusion

Genes comprising the σ^B regulon encode solute transporters, novel cell-wall proteins, universal stress proteins, transcriptional regulators and include those involved in osmoregulation, carbon metabolism, ribosome- and envelope-function, as well as virulence and niche-specific survival genes such as those involved in bile resistance and exclusion. Ten of the σ^B-positively regulated genes of L. monocytogenes are absent in L. innocua. A total of 75 σ^B-positively regulated listerial genes had homologs in B. subtilis, but only 33 have been previously described as being σ^B-regulated in B. subtilis even though both species share a highly conserved σ^B-dependent consensus sequence. A low overlap of genes may reflects adaptation of these bacteria to their respective environmental conditions.

Novel bacterial artificial chromosome vector pUvBBAC for use in studies of the functional genomics of Listeria spp

Bacterial artificial chromosome (BAC) vectors are important tools for microbial genome research. We constructed a novel BAC vector, pUvBBAC, for replication in both gram-negative and gram-positive bacterial hosts. The pUvBBAC vector was used to generate a BAC library for the facultative intracellular pathogen Listeria monocytogenes EGD-e. The library had insert sizes ranging from 68 to 178 kb. We identified two recombinant BACs from the L. monocytogenes pUvBBAC library that each contained the entire virulence gene cluster (vgc) of L. monocytogenes and transferred them to a nonpathogenic Listeria innocua strain. Recombinant L. innocua strains harboring pUvBBAC+vgc1 and pUvBBAC+vgc2 produced the vgc-specific listeriolysin (LLO) and actin assembly protein ActA and represent the first reported cloning of the vgc locus in its entirety. The use of the novel broad-host-range BAC vector pUvBBAC extends the versatility of this technology and provides a powerful platform for detailed functional genomics of gram-positive bacteria as well as its use in explorative functional metagenomics.