GET_HOMOLOGUES, a versatile software package for scalable and robust microbial pangenome analysis

Genomic characterization reconfirms the taxonomic status of Lactobacillus parakefiri

Whole-genome sequencing was performed for Lactobacillus parakefiri JCM 8573^T to confirm its hitherto controversial taxonomic position. Here, we report its first reliable reference genome. Genome-wide metrics, such as average nucleotide identity and digital DNA-DNA hybridization, and phylogenomic analysis based on multiple genes supported its taxonomic status as a distinct species in the genus Lactobacillus. The availability of a reliable genome sequence will aid future investigations on the industrial applications of L. parakefiri in functional foods such as kefir grains.

Ecotype diversification of an abundant Roseobacter lineage

The Roseobacter DC5-80-3 cluster (also known as the RCA clade) is among the most abundant bacterial lineages in temperate and polar oceans. Previous studies revealed two phylotypes within this cluster that are distinctly distributed in the Antarctic and other ocean provinces. Here, we report a nearly complete genome co-assembly of three closely related single cells co-occurring in the Antarctic, and compare it to the available genomes of the other phylotype from ocean regions where iron is more accessible but phosphorus and nitrogen are less. The Antarctic phylotype exclusively contains an operon structure consisting of a dicitrate transporter fecBCDE and an upstream regulator likely for iron uptake, whereas the other phylotype consistently carry a high-affinity phosphate pst transporter and the phoB-phoR regulatory system, a high-affinity ammonium amtB transporter, urea and taurine utilization systems. Moreover, the Antarctic phylotype uses proteorhodopsin to acquire light, whereas the other uses bacteriochlorophyll-a and the sulfur-oxidizing sox cluster for energy acquisition. This is potentially an iron-saving strategy for the Antarctic phylotype because only the latter two pathways have iron-requiring cytochromes. Therefore, the two DC5-80-3 phylotypes, while diverging by only 1.1% in their 16S rRNA genes, have evolved systematic differences in metabolism to support their distinct ecologies.

Patterns and processes of Mycobacterium bovis evolution revealed by phylogenomic analyses

Mycobacterium bovis is an important animal pathogen worldwide that parasitizes wild and domesticated vertebrate livestock as well as humans. A comparison of the five M. bovis complete genomes from the United Kingdom, South Korea, Brazil, and the United States revealed four novel large-scale structural variations of at least 2,000 bp. A comparative phylogenomic study including 2,483 core genes of 38 taxa from eight countries showed conflicting phylogenetic signal among sites. By minimizing this effect, we obtained a tree that better agrees with sampling locality. Results supported a relatively basal position of African strains (all isolated from Homo sapiens ), confirming that Africa was an important region for early diversification and that humans were one of the earliest hosts. Selection analyses revealed that functional categories such as “Lipid transport and metabolism,” “Cell cycle control, cell division, chromosome partitioning” and “Cell motility” were significant for the evolution of the group, besides other categories previously described, showing importance of genes associated with virulence and cholesterol metabolism in the evolution of M. bovis . PE/PPE genes, many of which are known to be associated with virulence, were major targets for large-scale polymorphisms, homologous recombination, and positive selection, evincing for the first time a plethora of evolutionary forces possibly contributing to differential adaptability in M. bovis . By assuming different priors, US strains originated and started to diversify around 150–5,210 ya. By further analyzing the largest set of US genomes to date (76 in total), obtained from 14 host species, we detected that hosts were not clustered in clades (except for a few cases), with some faster-evolving strains being detected, suggesting fast and ongoing reinfections across host species, and therefore, the possibility of new bovine tuberculosis outbreaks.

The first complete genome sequences of the acI lineage, the most abundant freshwater Actinobacteria, obtained by whole-genome-amplification of dilution-to-extinction cultures

The acI lineage of the phylum Actinobacteria is the most abundant bacterial group in most freshwater lakes. However, due to difficulties in laboratory cultivation, only two mixed cultures and some incomplete single-amplified or metagenome-derived genomes have been reported for the lineage. Here, we report the initial cultivation and complete genome sequences of four novel strains of the acI lineage from the tribes acI-A1, -A4, -A7, and -C1. The acI strains, initially isolated by dilution-to-extinction culturing, eventually failed to be maintained as axenic cultures. However, the first complete genomes of the acI lineage were successfully obtained from these initial cultures through whole genome amplification applied to more than hundreds of cultured acI cells. The genome sequences exhibited features of genome streamlining and showed that the strains are aerobic chemoheterotrophs sharing central metabolic pathways, with some differences among tribes that may underlie niche diversification within the acI lineage. Actinorhodopsin was found in all strains, but retinal biosynthesis was complete in only A1 and A4 tribes.

Prophages and adaptation of Staphylococcus aureus ST398 to the human clinic

BACKGROUND

It has been suggested that prophages in the ST398 S. aureus clone are responsible for expanding ST398's spectrum of action and increasing its ability to cause human infections. We carried out the first characterization of the various prophages carried by 76 ST398 bloodstream infection (BSI) isolates obtained over 9 years of observation.

RESULTS

Whole-genome sequencing of 22 representative isolates showed (1) the presence of the φ3-prophage and diverse genetic features typical of animal-associated isolates (i.e., SCCmec XI element, Tn916 transposon and non φ3-prophages) in a majority of BSI isolates, (2) one BSI isolate devoid of the φ3-prophage but otherwise similar to an animal-infecting isolate, (3) 35 prophages carrying numerous genes previously associated with virulence or immune evasion in animal models of staphylococcal infections. The analysis of prophage content in all 76 BSI isolates showed an increasing prevalence of polylysogeny over time. Overall, over the course of the last 10 years, the BSI isolates appear to have acquired increasing numbers of genetic features previously shown to contribute to bacterial adaptation and virulence in animal models of staphylococcal infections.

CONCLUSIONS

We hypothesize that lysogeny has played a significant role in increasing the ability of the ST398 clone to cause infections in humans. Our findings highlight the risk that the ST398 lineage will increase its threat to public health by continuing to acquire virulence and/or multiple antibiotic-resistance genes from hospital-associated clones of Staphylococcus aureus.

Genomic Characterization of Dairy Associated Leuconostoc Species and Diversity of Leuconostocs in Undefined Mixed Mesophilic Starter Cultures

Undefined mesophilic mixed (DL-type) starter cultures are composed of predominantly Lactococcus lactis subspecies and 1-10% Leuconostoc spp. The composition of the Leuconostoc population in the starter culture ultimately affects the characteristics and the quality of the final product. The scientific basis for the taxonomy of dairy relevant leuconostocs can be traced back 50 years, and no documentation on the genomic diversity of leuconostocs in starter cultures exists. We present data on the Leuconostoc population in five DL-type starter cultures commonly used by the dairy industry. The analyses were performed using traditional cultivation methods, and further augmented by next-generation DNA sequencing methods. Bacterial counts for starter cultures cultivated on two different media, MRS and MPCA, revealed large differences in the relative abundance of leuconostocs. Most of the leuconostocs in two of the starter cultures were unable to grow on MRS, emphasizing the limitations of culture-based methods and the importance of careful media selection or use of culture independent methods. Pan-genomic analysis of 59 Leuconostoc genomes enabled differentiation into twelve robust lineages. The genomic analyses show that the dairy-associated leuconostocs are highly adapted to their environment, characterized by the acquisition of genotype traits, such as the ability to metabolize citrate. In particular, Leuconostoc mesenteroides subsp. cremoris display telltale signs of a degenerative evolution, likely resulting from a long period of growth in milk in association with lactococci. Great differences in the metabolic potential between Leuconostoc species and subspecies were revealed. Using targeted amplicon sequencing, the composition of the Leuconostoc population in the five commercial starter cultures was shown to be significantly different. Three of the cultures were dominated by Ln. mesenteroides subspecies cremoris. Leuconostoc pseudomesenteroides dominated in two of the cultures while Leuconostoc lactis, reported to be a major constituent in fermented dairy products, was only present in low amounts in one of the cultures. This is the first in-depth study of Leuconostoc genomics and diversity in dairy starter cultures. The results and the techniques presented may be of great value for the dairy industry.

Genomic Epidemiology of Penicillin-Nonsusceptible Pneumococci with Nonvaccine Serotypes Causing Invasive Disease in the United States

Conjugate vaccination against seven pneumococcal serotypes (PCV7) reduced disease prevalence due to antibiotic-resistant strains throughout the 2000s. However, diseases caused by resistant nonvaccine type (NVT) strains increased. Some of these emerging strains were derived from vaccine types (VT) that had changed their capsule by recombination. The introduction of a vaccine targeting 13 serotypes (PCV13) in 2010 has led to concern that this scenario will repeat itself. We generated high-quality draft genomes from 265 isolates of NVT pneumococci not susceptible to penicillin (PNSP) in 2009 and compared them with the genomes of 581 isolates from 2012 to 2013 collected by the Active Bacterial Core surveillance (ABCs) of the Centers for Disease Control and Prevention (CDC). Of the seven sequence clusters (SCs) identified, three SCs fell into a single lineage associated with serogroup 23, which had an origin in 1908 as dated by coalescent analysis and included isolates with a divergent 23B capsule locus. Three other SCs represented relatively deep-branching lineages associated with serotypes 35B, 15A, and 15BC. In all cases, the resistant clones originated prior to 2010, indicating that PNSP are at present dominated by descendants of NVT clones present before vaccination. With one exception (15BC/ST3280), these SCs were related to clones identified by the Pneumococcal Molecular Epidemiology Network (PMEN). We conclude that postvaccine diversity in NVT PNSP between 2009 and 2013 was driven mainly by the persistence of preexisting strains rather than through de novo adaptation, with few cases of serotype switching. Future surveillance is essential for documenting the long-term dynamics and resistance of NVT PNSP.

Genome Analysis of Clostridium difficile PCR Ribotype 014 Lineage in Australian Pigs and Humans Reveals a Diverse Genetic Repertoire and Signatures of Long-Range Interspecies Transmission

Clostridium difficile PCR ribotype (RT) 014 is well-established in both human and porcine populations in Australia, raising the possibility that C. difficile infection (CDI) may have a zoonotic or foodborne etiology. Here, whole genome sequencing and high-resolution core genome phylogenetics were performed on a contemporaneous collection of 40 Australian RT014 isolates of human and porcine origin. Phylogenies based on MLST (7 loci, STs 2, 13, and 49) and core orthologous genes (1260 loci) showed clustering of human and porcine strains indicative of very recent shared ancestry. Core genome single nucleotide variant (SNV) analysis found 42% of human strains showed a clonal relationship (separated by ≤2 SNVs in their core genome) with one or more porcine strains, consistent with recent inter-host transmission. Clones were spread over a vast geographic area with 50% of the human cases occurring without recent healthcare exposure. These findings suggest a persistent community reservoir with long-range dissemination, potentially due to agricultural recycling of piggery effluent. We also provide the first pan-genome analysis for this lineage, characterizing its resistome, prophage content, and in silico virulence potential. The RT014 is defined by a large "open" pan-genome (7587 genes) comprising a core genome of 2296 genes (30.3% of the total gene repertoire) and an accessory genome of 5291 genes. Antimicrobial resistance genotypes and phenotypes varied across host populations and ST lineages and were characterized by resistance to tetracycline [tetM, tetA(P), tetB(P) and tetW], clindamycin/erythromycin (ermB), and aminoglycosides (aph3-III-Sat4A-ant6-Ia). Resistance was mediated by clinically important mobile genetic elements, most notably Tn6194 (harboring ermB) and a novel variant of Tn5397 (harboring tetM). Numerous clinically important prophages (Siphoviridae and Myoviridae) were identified as well as an uncommon accessory gene regulator locus (agr3). Conservation in the pathogenicity locus and S-layer correlated with ST affiliation, further extending the concept of clonal C. difficile lineages. This study provides novel insights on the genetic variability and strain relatedness of C. difficile RT014, a lineage of emerging One Health importance. Ongoing molecular and genomic surveillance of strains in humans, animals, food, and the environment is imperative to identify opportunities to reduce the overall CDI burden.

The driving force of prophages and CRISPR-Cas system in the evolution of Cronobacter sakazakii

Cronobacter sakazakii is an important foodborne pathogens causing rare but life-threatening diseases in neonates and infants. CRISPR-Cas system is a new prokaryotic defense system that provides adaptive immunity against phages, latter play an vital role on the evolution and pathogenicity of host bacteria. In this study, we found that genome sizes of C. sakazakii strains had a significant positive correlation with total genome sizes of prophages. Prophages contributed to 16.57% of the genetic diversity (pan genome) of C. sakazakii, some of which maybe the potential virulence factors. Subtype I-E CRISPR-Cas system and five types of CRISPR arrays were found in the conserved site of C. sakazakii strains. CRISPR1 and CRISPR2 loci with high variable spacers were active and showed potential protection against phage attacks. The number of spacers from two active CRISPR loci in clinical strains was significant less than that of foodborne strains, it maybe a reason why clinical strains were found to have more prophages than foodborne strains. The frequently gain/loss of prophages and spacers in CRISPR loci is likely to drive the quick evolution of C. sakazakii. Our study provides a new insight into the co-evolution of phages and C. sakazakii.

Analysis of Plant Pan-Genomes and Transcriptomes with GET_HOMOLOGUES-EST, a Clustering Solution for Sequences of the Same Species

The pan-genome of a species is defined as the union of all the genes and non-coding sequences found in all its individuals. However, constructing a pan-genome for plants with large genomes is daunting both in sequencing cost and the scale of the required computational analysis. A more affordable alternative is to focus on the genic repertoire by using transcriptomic data. Here, the software GET_HOMOLOGUES-EST was benchmarked with genomic and RNA-seq data of 19 Arabidopsis thaliana ecotypes and then applied to the analysis of transcripts from 16 Hordeum vulgare genotypes. The goal was to sample their pan-genomes and classify sequences as core, if detected in all accessions, or accessory, when absent in some of them. The resulting sequence clusters were used to simulate pan-genome growth, and to compile Average Nucleotide Identity matrices that summarize intra-species variation. Although transcripts were found to under-estimate pan-genome size by at least 10%, we concluded that clusters of expressed sequences can recapitulate phylogeny and reproduce two properties observed in A. thaliana gene models: accessory loci show lower expression and higher non-synonymous substitution rates than core genes. Finally, accessory sequences were observed to preferentially encode transposon components in both species, plus disease resistance genes in cultivated barleys, and a variety of protein domains from other families that appear frequently associated with presence/absence variation in the literature. These results demonstrate that pan-genome analyses are useful to explore germplasm diversity.

SEC-Translocon Dependent Extracytoplasmic Proteins of Candidatus Liberibacter asiaticus

Citrus Huanglongbing (HLB) is the most destructive citrus disease worldwide. HLB is associated with three species of the phloem-limited, gram-negative, fastidious α-proteobacteria: Candidatus Liberibacter asiaticus (Las), Ca. L. americanus (Lam), and Ca. L. africanus (Laf) with Las being the most widespread species. Las has not been cultured in artificial media, which has greatly hampered our efforts to understand its virulence mechanisms. Las contains a complete Sec-translocon, which has been suggested to transport Las proteins including virulence factors into the extracytoplasmic milieu. In this study, we characterized the Sec-translocon dependent, signal peptide containing extracytoplasmic proteins of Las. A total of 166 proteins of Las-psy62 strain were predicted to contain signal peptides targeting them out of the cell cytoplasm via the Sec-translocon using LipoP, SigalP 3.0, SignalP 4.1, and Phobius. We also predicated SP containing extracytoplasmic proteins for Las-gxpsy and Las-Ishi-1, Lam, Laf, Ca. L. solanacearum (Lso), and L. crescens (Lcr). For experimental validation of the predicted extracytoplasmic proteins, Escherichia coli based alkaline phosphatase (PhoA) gene fusion assays were conducted. A total of 86 out of the 166 predicted Las proteins were experimentally validated to contain signal peptides. Additionally, Las-psy62 lepB (CLIBASIA_04190), the gene encodes signal peptidase I, was able to partially complement the amber mutant of lepB of E. coli. This work will contribute to the identification of Sec-translocon dependent effector proteins of Las, which might be involved in virulence of Las.

Characterisation of Commensal Escherichia coli Isolated from Apparently Healthy Cattle and Their Attendants in Tanzania

While pathogenic types of Escherichia coli are well characterized, relatively little is known about the commensal E. coli flora. In the current study, antimicrobial resistance in commensal E. coli and distribution of ERIC-PCR genotypes among isolates of such bacteria from cattle and cattle attendants on cattle farms in Tanzania were investigated. Seventeen E. coli genomes representing different ERIC-PCR types of commensal E. coli were sequenced in order to determine their possible importance as a reservoir for both antimicrobial resistance genes and virulence factors. Both human and cattle isolates were highly resistant to tetracycline (40.8% and 33.1%), sulphamethazole-trimethoprim (49.0% and 8.8%) and ampicillin (44.9% and 21.3%). However, higher proportion of resistant E. coli and higher frequency of resistance to more than two antimicrobials was found in isolates from cattle attendants than isolates from cattle. Sixteen out of 66 ERIC-PCR genotypes were shared between the two hosts, and among these ones, seven types contained isolates from cattle and cattle attendants from the same farm, suggesting transfer of strains between hosts. Genome-wide analysis showed that the majority of the sequenced cattle isolates were assigned to phylogroups B1, while human isolates represented phylogroups A, C, D and E. In general, in silico resistome and virulence factor identification did not reveal differences between hosts or phylogroups, except for lpfA and iss found to be cattle and B1 phylogroup specific. The most frequent plasmids replicon genes found in strains from both hosts were of IncF type, which are commonly associated with carriage of antimicrobial and virulence genes. Commensal E. coli from cattle and attendants were found to share same genotypes and to carry antimicrobial resistance and virulence genes associated with both intra and extraintestinal E. coli pathotypes.

Evolution of a zoonotic pathogen: investigating prophage diversity in enterohaemorrhagic Escherichia coli O157 by long-read sequencing

Enterohaemorrhagic Escherichia coli (EHEC) O157 is a zoonotic pathogen for which colonization of cattle and virulence in humans is associated with multiple horizontally acquired genes, the majority present in active or cryptic prophages. Our understanding of the evolution and phylogeny of EHEC O157 continues to develop primarily based on core genome analyses; however, such short-read sequences have limited value for the analysis of prophage content and its chromosomal location. In this study, we applied Single Molecule Real Time (SMRT) sequencing, using the Pacific Biosciences long-read sequencing platform, to isolates selected from the main sub-clusters of this clonal group. Prophage regions were extracted from these sequences and from published reference strains. Genome position and prophage diversity were analysed along with genetic content. Prophages could be assigned to clusters, with smaller prophages generally exhibiting less diversity and preferential loss of structural genes. Prophages encoding Shiga toxin (Stx) 2a and Stx1a were the most diverse, and more variable compared to prophages encoding Stx2c, further supporting the hypothesis that Stx2c-prophage integration was ancestral to acquisition of other Stx types. The concept that phage type (PT) 21/28 (Stx2a+, Stx2c+) strains evolved from PT32 (Stx2c+) was supported by analysis of strains with excised Stx-encoding prophages. Insertion sequence elements were over-represented in prophage sequences compared to the rest of the genome, showing integration in key genes such as stx and an excisionase, the latter potentially acting to capture the bacteriophage into the genome. Prophage profiling should allow more accurate prediction of the pathogenic potential of isolates.

Comparative genomics of Beauveria bassiana: uncovering signatures of virulence against mosquitoes

BACKGROUND

Entomopathogenic fungi such as Beauveria bassiana are promising biological agents for control of malaria mosquitoes. Indeed, infection with B. bassiana reduces the lifespan of mosquitoes in the laboratory and in the field. Natural isolates of B. bassiana show up to 10-fold differences in virulence between the most and the least virulent isolate. In this study, we sequenced the genomes of five isolates representing the extremes of low/high virulence and three RNA libraries, and applied a genome comparison approach to uncover genetic mechanisms underpinning virulence.

RESULTS

A high-quality, near-complete genome assembly was achieved for the highly virulent isolate Bb8028, which was compared to the assemblies of the four other isolates. Whole genome analysis showed a high level of genetic diversity between the five isolates (2.85-16.8 SNPs/kb), which grouped into two distinct phylogenetic clusters. Mating type gene analysis revealed the presence of either the MAT1-1-1 or the MAT1-2-1 gene. Moreover, a putative new MAT gene (MAT1-2-8) was detected in the MAT1-2 locus. Comparative genome analysis revealed that Bb8028 contains 163 genes exclusive for this isolate. These unique genes have a tendency to cluster in the genome and to be often located near the telomeres. Among the genes unique to Bb8028 are a Non-Ribosomal Peptide Synthetase (NRPS) secondary metabolite gene cluster, a polyketide synthase (PKS) gene, and five genes with homology to bacterial toxins. A survey of candidate virulence genes for B. bassiana is presented.

CONCLUSIONS

Our results indicate several genes and molecular processes that may underpin virulence towards mosquitoes. Thus, the genome sequences of five isolates of B. bassiana provide a better understanding of the natural variation in virulence and will offer a major resource for future research on this important biological control agent.

Overexpressed Proteins in Hypervirulent Clade 8 and Clade 6 Strains of Escherichia coli O157:H7 Compared to E. coli O157:H7 EDL933 Clade 3 Strain

Escherichia coli O157:H7 is responsible for severe diarrhea and hemolytic uremic syndrome (HUS), and predominantly affects children under 5 years. The major virulence traits are Shiga toxins, necessary to develop HUS and the Type III Secretion System (T3SS) through which bacteria translocate effector proteins directly into the host cell. By SNPs typing, E. coli O157:H7 was separated into nine different clades. Clade 8 and clade 6 strains were more frequently associated with severe disease and HUS. In this study, we aimed to identify differentially expressed proteins in two strains of E. coli O157:H7 (clade 8 and clade 6), obtained from cattle and compared them with the well characterized reference EDL933 strain (clade 3). Clade 8 and clade 6 strains show enhanced pathogenicity in a mouse model and virulence-related properties. Proteins were extracted and analyzed using the TMT-6plex labeling strategy associated with two dimensional liquid chromatography and mass spectrometry in tandem. We detected 2241 proteins in the cell extract and 1787 proteins in the culture supernatants. Attention was focused on the proteins related to virulence, overexpressed in clade 6 and 8 strains compared to EDL933 strain. The proteins relevant overexpressed in clade 8 strain were the curli protein CsgC, a transcriptional activator (PchE), phage proteins, Stx2, FlgM and FlgD, a dienelactone hydrolase, CheW and CheY, and the SPATE protease EspP. For clade 6 strain, a high overexpression of phage proteins was detected, mostly from Stx2 encoding phage, including Stx2, flagellin and the protease TagA, EDL933_p0016, dienelactone hydrolase, and Haemolysin A, amongst others with unknown function. Some of these proteins were analyzed by RT-qPCR to corroborate the proteomic data. Clade 6 and clade 8 strains showed enhanced transcription of 10 out of 12 genes compared to EDL933. These results may provide new insights in E. coli O157:H7 mechanisms of pathogenesis.

Genome of the carbapenemase-producing clinical isolate Elizabethkingia miricola EM_CHUV and comparative genomics with Elizabethkingia meningoseptica and Elizabethkingia anophelis: evidence for intrinsic multidrug resistance trait of emerging pathogens

Elizabethkingia miricola is a Gram-negative non-fermenting rod emerging as a life-threatening human pathogen. The multidrug-resistant (MDR) carbapenemase-producing clinical isolate E. miricola EM_CHUV was recovered in the setting of severe nosocomial pneumonia. In this study, the genome of E. miricola EM_CHUV was sequenced and a functional analysis was performed, including a comparative genomic study with Elizabethkingia meningoseptica and Elizabethkingia anophelis. The resistome of EM_CHUV revealed the presence of a high number of resistance genes, including the presence of the bla_GOB-13 and bla_B-9 carbapenemase-encoding genes. Twelve mobility genes, with only two of them located in the proximity of resistance genes, and four potential genomic islands were identified in the genome of EM_CHUV, but no prophages or CRISPR sequences. Ten restriction-modification system (RMS) genes were also identified. In addition, we report the presence of a putative conjugative plasmid (pEM_CHUV) that does not encode any antibiotic resistance genes. Altogether, these findings point towards a limited number of DNA exchanges with other bacteria and suggest that multidrug resistance is an intrinsic trait of E. miricola owing to the presence of a high number of resistance genes within the bacterial core genome.

Kaumoebavirus, a New Virus That Clusters with Faustoviruses and Asfarviridae

In this study, we report the isolation of a new giant virus found in sewage water from the southern area of Jeddah (Saudi Arabia), with morphological and genomic resemblance to Faustoviruses. This new giant virus, named Kaumoebavirus, was obtained from co-culture with Vermamoeba vermiformis, an amoeboid protozoa considered to be of special interest to human health and the environment. This new virus has ~250 nm icosahedral capsids and a 350,731 bp DNA genome length. The genome of Kaumoebavirus has a coding density of 86%, corresponding to 465 genes. Most of these genes (59%) are closely related to genes from members of the proposed order Megavirales, and the best matches to its proteins with other members of the Megavirales are Faustoviruses (43%) and Asfarviruses (23%). Unsurprisingly, phylogenetic reconstruction places Kaumoebavirus as a distant relative of Faustoviruses and Asfarviruses.

Non-toxigenic environmental Vibrio cholerae O1 strain from Haiti provides evidence of pre-pandemic cholera in Hispaniola

Vibrio cholerae is ubiquitous in aquatic environments, with environmental toxigenic V. cholerae O1 strains serving as a source for recurrent cholera epidemics and pandemic disease. However, a number of questions remain about long-term survival and evolution of V. cholerae strains within these aquatic environmental reservoirs. Through monitoring of the Haitian aquatic environment following the 2010 cholera epidemic, we isolated two novel non-toxigenic (ctxA/B-negative) Vibrio cholerae O1. These two isolates underwent whole-genome sequencing and were investigated through comparative genomics and Bayesian coalescent analysis. These isolates cluster in the evolutionary tree with strains responsible for clinical cholera, possessing genomic components of 6^th and 7^th pandemic lineages, and diverge from "modern" cholera strains around 1548 C.E. [95% HPD: 1532-1555]. Vibrio Pathogenicity Island (VPI)-1 was present; however, SXT/R391-family ICE and VPI-2 were absent. Rugose phenotype conversion and vibriophage resistance evidenced adaption for persistence in aquatic environments. The identification of V. cholerae O1 strains in the Haitian environment, which predate the first reported cholera pandemic in 1817, broadens our understanding of the history of pandemics. It also raises the possibility that these and similar environmental strains could acquire virulence genes from the 2010 Haitian epidemic clone, including the cholera toxin producing CTXϕ.

Pangenome and Phylogenomic Analysis of the Pathogenic Actinobacterium Rhodococcus equi

We report a comparative study of 29 representative genomes of the animal pathogen Rhodococcus equi The analyses showed that R. equi is genetically homogeneous and clonal, with a large core genome accounting for ≈80% of an isolates' gene content. An open pangenome, even distribution of accessory genes among the isolates, and absence of significant core-genome recombination, indicated that gene gain/loss is a main driver of R. equi genome evolution. Traits previously predicted to be important in R. equi physiology, virulence and niche adaptation were part of the core genome. This included the lack of a phosphoenolpyruvate:carbohydrate transport system (PTS), unique among the rhodococci except for the closely related Rhodococcus defluvii, reflecting selective PTS gene loss in the R. equi-R. defluvii sublineage. Thought to be asaccharolytic, rbsCB and glcP non-PTS sugar permease homologues were identified in the core genome and, albeit inefficiently, R. equi utilized their putative substrates, ribose and (irregularly) glucose. There was no correlation between R. equi whole-genome phylogeny and host or geographical source, with evidence of global spread of genomovars. The distribution of host-associated virulence plasmid types was consistent with the exchange of the plasmids (and corresponding host shifts) across the R. equi population, and human infection being zoonotically acquired. Phylogenomic analyses demonstrated that R. equi occupies a central position in the Rhodococcus phylogeny, not supporting the recently proposed transfer of the species to a new genus.

Genomic characterization of a fructophilic bee symbiont Lactobacillus kunkeei reveals its niche-specific adaptation

Lactobacillus kunkeei is classified as a sole obligate fructophilic lactic acid bacterium that is found in fructose-rich niches, including the guts of honeybees. The species is differentiated from other lactobacilli based on its poor growth with glucose, enhanced growth in the presence of oxygen and other electron acceptors, and production of high concentrations of acetate from the metabolism of glucose. These characteristics are similar to phylogenetically distant Fructobacillus spp. In the present study, the genomic structure of L. kunkeei was characterized by using 16 different strains, and it had significantly less genes and smaller genomes when compared with other lactobacilli. Functional gene classification revealed that L. kunkeei had lost genes specifically involved in carbohydrate transport and metabolism. The species also lacked most of the genes for respiration, although growth was enhanced in the presence of oxygen. The adhE gene of L. kunkeei, encoding a bifunctional alcohol dehydrogenase (ADH)/aldehyde dehydrogenase (ALDH) protein, lacked the part encoding the ADH domain, which is reported here for the first time in lactic acid bacteria. The deletion resulted in the lack of ADH activity, implying a requirement for electron acceptors in glucose assimilation. These results clearly indicated that L. kunkeei had undergone a specific reductive evolution in order to adapt to fructose-rich environments. The reduction characteristics were similar to those of Fructobacillus spp., but distinct from other lactobacilli with small genomes, such as Lactobacillus gasseri and Lactobacillus vaginalis. Fructose-richness thus induced an environment-specific gene reduction in phylogenetically distant microorganisms.

Effector Diversification Contributes to Xanthomonas oryzae pv. oryzae Phenotypic Adaptation in a Semi-Isolated Environment

Understanding the processes that shaped contemporary pathogen populations in agricultural landscapes is quite important to define appropriate management strategies and to support crop improvement efforts. Here, we took advantage of an historical record to examine the adaptation pathway of the rice pathogen Xanthomonas oryzae pv. oryzae (Xoo) in a semi-isolated environment represented in the Philippine archipelago. By comparing genomes of key Xoo groups we showed that modern populations derived from three Asian lineages. We also showed that diversification of virulence factors occurred within each lineage, most likely driven by host adaptation, and it was essential to shape contemporary pathogen races. This finding is particularly important because it expands our understanding of pathogen adaptation to modern agriculture.

The Genetic Analysis of an Acinetobacter johnsonii Clinical Strain Evidenced the Presence of Horizontal Genetic Transfer

Acinetobacter johnsonii rarely causes human infections. While most A. johnsonii isolates are susceptible to virtually all antibiotics, strains harboring a variety of β-lactamases have recently been described. An A. johnsonii Aj2199 clinical strain recovered from a hospital in Buenos Aires produces PER-2 and OXA-58. We decided to delve into its genome by obtaining the whole genome sequence of the Aj2199 strain. Genome comparison studies on Aj2199 revealed 240 unique genes and a close relation to strain WJ10621, isolated from the urine of a patient in China. Genomic analysis showed evidence of horizontal genetic transfer (HGT) events. Forty-five insertion sequences and two intact prophages were found in addition to several resistance determinants such as blaPER-2, blaOXA-58, blaTEM-1, strA, strB, ereA, sul1, aacC2 and a new variant of blaOXA-211, called blaOXA-498. In particular, blaPER-2 and blaTEM-1 are present within the typical contexts previously described in the Enterobacteriaceae family. These results suggest that A. johnsonii actively acquires exogenous DNA from other bacterial species and concomitantly becomes a reservoir of resistance genes.

Comparative Genomics Analysis of Streptococcus tigurinus Strains Identifies Genetic Elements Specifically and Uniquely Present in Highly Virulent Strains

Streptococcus tigurinus is responsible for severe invasive infections such as infective endocarditis, spondylodiscitis and meningitis. As described, S. tigurinus isolates AZ_3aT and AZ_14 were highly virulent (HV phenotype) in an experimental model of infective endocarditis and showed enhanced adherence and invasion of human endothelial cells when compared to low virulent S. tigurinus isolate AZ_8 (LV phenotype). Here, we sought whether genetic determinants could explain the higher virulence of AZ_3aT and AZ_14 isolates. Several genetic determinants specific to the HV strains were identified through extensive comparative genomics amongst which some were thought to be highly relevant for the observed HV phenotype. These included i) an iron uptake and metabolism operon, ii) an ascorbate assimilation operon, iii) a newly acquired PI-2-like pilus islets described for the first time in S. tigurinus, iv) a hyaluronate metabolism operon, v) an Entner-Doudoroff pathway of carbohydrates metabolism, and vi) an alternate pathways for indole biosynthesis. We believe that the identified genomic features could largely explain the phenotype of high infectivity of the two HV S. tigurinus strains. Indeed, these features include determinants that could be involved at different stages of the disease such as survival of S. tigurinus in blood (iron uptake and ascorbate metabolism operons), initial attachment of bacterial pathogen to the damaged cardiac tissue and/or vegetation that formed on site (PI-2-like pilus islets), tissue invasion (hyaluronate operon and Entner-Doudoroff pathway) and regulation of pathogenicity (indole biosynthesis pathway).

Evolutionary architecture of the infant-adapted group of Bifidobacterium species associated with the probiotic function

Bifidobacteria, often associated with the gastrointestinal tract of animals, are well known for their roles as probiotics. Among the dozens of Bifidobacterium species, Bifidobacterium bifidum, B. breve, and B. longum are the ones most frequently isolated from the feces of infants and known to help the digestion of human milk oligosaccharides. To investigate the correlation between the metabolic properties of bifidobacteria and their phylogeny, we performed a phylogenomic analysis based on 452 core genes of forty-four completely sequenced Bifidobacterium species. Results show that a major evolutionary event leading to the clade of the infant-adapted species is linked to carbohydrate metabolism, but it is not the only factor responsible for the adaptation of bifidobacteria to the gut. The genome of B. longum subsp. infantis, a typical bifidobacterium in the gut of breast-fed infants, encodes proteins associated with several kinds of species-specific metabolic pathways, including urea metabolism and biosynthesis of riboflavin and lantibiotics. Our results demonstrate that these metabolic features, which are associated with the probiotic function of bifidobacteria, are species-specific and highly correlate with their phylogeny.

Identification of source and sink populations for the emergence and global spread of the East-Asia clone of community-associated MRSA

BACKGROUND

Our understanding of the factors influencing the emergence, dissemination and global distribution of epidemic clones of bacteria is limited. ST59 is a major epidemic clone of community-associated MRSA in East Asia, responsible for extensive morbidity and mortality, but has a much lower prevalence in other parts of the world. The geographic origin of ST59 and its international routes of dissemination are unclear and disputed in the literature.

RESULTS

To investigate the origin and spread of the ST59 clone, we obtained whole genome sequences of isolates from four continents, sampled over more than a decade, and carried out a time-scaled phylogeographic analysis. We discover that two distinct ST59 clades emerged concurrently, in East Asia and the USA, but underwent clonal expansion at different times. The East Asia clade was strongly enriched for gene determinants associated with antibiotic resistance, consistent with regional differences in antibiotic usage. Both clones spread independently to Australia and Europe, and we found evidence of the persistence of multi-drug resistance following export from East Asia. Direct transfer of strains between Taiwan and the USA was not observed in either direction, consistent with geographic niche exclusion.

CONCLUSIONS

Our results resolve a longstanding controversy regarding the origin of the ST59 clone, revealing the major global source and sink populations and routes for the spread of multi-drug resistant clones. Additionally, our findings indicate that diversification of the accessory genome of epidemic clones partly reflects region-specific patterns of antibiotic usage, which may influence bacterial fitness after transmission to different geographic locations.

Genome Evolution and Nitrogen Fixation in Bacterial Ectosymbionts of a Protist Inhabiting Wood-Feeding Cockroaches

By combining genomics and isotope imaging analysis using high-resolution secondary ion mass spectrometry (NanoSIMS), we examined the function and evolution of Bacteroidales ectosymbionts of the protist Barbulanympha from the hindguts of the wood-eating cockroach Cryptocercus punctulatus In particular, we investigated the structure of ectosymbiont genomes, which, in contrast to those of endosymbionts, has been little studied to date, and tested the hypothesis that these ectosymbionts fix nitrogen. Unlike with most obligate endosymbionts, genome reduction has not played a major role in the evolution of the Barbulanympha ectosymbionts. Instead, interaction with the external environment has remained important for this symbiont as genes for synthesis of transporters, outer membrane proteins, lipopolysaccharides, and lipoproteins have been retained. The ectosymbiont genome carried two complete operons for nitrogen fixation, a urea transporter, and a urease, indicating the availability of nitrogen as a driving force behind the symbiosis. NanoSIMS analysis of C. punctulatus hindgut symbionts exposed in vivo to (15)N2 supports the hypothesis that Barbulanympha ectosymbionts are capable of nitrogen fixation. This genomic and in vivo functional investigation of protist ectosymbionts highlights the diversity of evolutionary forces and trajectories that shape symbiotic interactions.

IMPORTANCE

The ecological and evolutionary importance of symbioses is increasingly clear, but the overall diversity of symbiotic interactions remains poorly explored. In this study, we investigated the evolution and nitrogen fixation capabilities of ectosymbionts attached to the protist Barbulanympha from the hindgut of the wood-eating cockroach Cryptocercus punctulatus In addressing genome evolution of protist ectosymbionts, our data suggest that the ecological pressures influencing the evolution of extracellular symbionts clearly differ from intracellular symbionts and organelles. Using NanoSIMS analysis, we also obtained direct imaging evidence of a specific hindgut microbe playing a role in nitrogen fixation. These results demonstrate the power of combining NanoSIMS and genomics tools for investigating the biology of uncultivable microbes. This investigation paves the way for a more precise understanding of microbial interactions in the hindguts of wood-eating insects and further exploration of the diversity and ecological significance of symbiosis between microbes.

DFAST and DAGA: web-based integrated genome annotation tools and resources

Quality assurance and correct taxonomic affiliation of data submitted to public sequence databases have been an everlasting problem. The DDBJ Fast Annotation and Submission Tool (DFAST) is a newly developed genome annotation pipeline with quality and taxonomy assessment tools. To enable annotation of ready-to-submit quality, we also constructed curated reference protein databases tailored for lactic acid bacteria. DFAST was developed so that all the procedures required for DDBJ submission could be done seamlessly online. The online workspace would be especially useful for users not familiar with bioinformatics skills. In addition, we have developed a genome repository, DFAST Archive of Genome Annotation (DAGA), which currently includes 1,421 genomes covering 179 species and 18 subspecies of two genera, Lactobacillus and Pediococcus, obtained from both DDBJ/ENA/GenBank and Sequence Read Archive (SRA). All the genomes deposited in DAGA were annotated consistently and assessed using DFAST. To assess the taxonomic position based on genomic sequence information, we used the average nucleotide identity (ANI), which showed high discriminative power to determine whether two given genomes belong to the same species. We corrected mislabeled or misidentified genomes in the public database and deposited the curated information in DAGA. The repository will improve the accessibility and reusability of genome resources for lactic acid bacteria. By exploiting the data deposited in DAGA, we found intraspecific subgroups in Lactobacillus gasseri and Lactobacillus jensenii, whose variation between subgroups is larger than the well-accepted ANI threshold of 95% to differentiate species. DFAST and DAGA are freely accessible at https://dfast.nig.ac.jp.

Isolation and Genomic Characterization of 'Desulfuromonas soudanensis WTL', a Metal- and Electrode-Respiring Bacterium from Anoxic Deep Subsurface Brine

Reaching a depth of 713 m below the surface, the Soudan Underground Iron Mine (Soudan, MN, USA) transects a massive Archaean (2.7 Ga) banded iron formation, providing a remarkably accessible window into the terrestrial deep biosphere. Despite organic carbon limitation, metal-reducing microbial communities are present in potentially ancient anoxic brines continuously emanating from exploratory boreholes on Level 27. Using graphite electrodes deposited in situ as bait, we electrochemically enriched and isolated a novel halophilic iron-reducing Deltaproteobacterium, ‘Desulfuromonas soudanensis’ strain WTL, from an acetate-fed three-electrode bioreactor poised at +0.24 V (vs. standard hydrogen electrode). Cyclic voltammetry revealed that ‘D. soudanensis’ releases electrons at redox potentials approximately 100 mV more positive than the model freshwater surface isolate Geobacter sulfurreducens, suggesting that its extracellular respiration is tuned for higher potential electron acceptors. ‘D. soudanensis’ contains a 3,958,620-bp circular genome, assembled to completion using single-molecule real-time (SMRT) sequencing reads, which encodes a complete TCA cycle, 38 putative multiheme c-type cytochromes, one of which contains 69 heme-binding motifs, and a LuxI/LuxR quorum sensing cassette that produces an unidentified N-acyl homoserine lactone. Another cytochrome is predicted to lie within a putative prophage, suggesting that horizontal gene transfer plays a role in respiratory flexibility among metal reducers. Isolation of ‘D. soudanensis’ underscores the utility of electrode-based approaches for enriching rare metal reducers from a wide range of habitats.

Complete genome sequence of Leuconostoc gelidum subsp. gasicomitatum KG16-1, isolated from vacuum-packaged vegetable sausages

Leuconostoc gelidum subsp. gasicomitatum is a predominant lactic acid bacterium (LAB) in spoilage microbial communities of different kinds of modified-atmosphere packaged (MAP) food products. So far, only one genome sequence of a poultry-originating type strain of this bacterium (LMG 18811(T)) has been available. In the current study, we present the completely sequenced and functionally annotated genome of strain KG16-1 isolated from a vegetable-based product. In addition, six other vegetable-associated strains were sequenced to study possible "niche" specificity suggested by recent multilocus sequence typing. The genome of strain KG16-1 consisted of one circular chromosome and three plasmids, which together contained 2,035 CDSs. The chromosome carried at least three prophage regions and one of the plasmids encoded a galactan degradation cluster, which might provide a survival advantage in plant-related environments. The genome comparison with LMG 18811(T) and six other vegetable strains suggests no major differences between the meat- and vegetable-associated strains that would explain their "niche" specificity. Finally, the comparison with the genomes of other leuconostocs highlights the distribution of functionally interesting genes across the L. gelidum strains and the genus Leuconostoc.

Comparative Genomics of Herpesviridae Family to Look for Potential Signatures of Human Infecting Strains

Herpesviridae family is one of the significant viral families which comprises major pathogens of a wide range of hosts. This family includes at least eight species of viruses which are known to infect humans. This family has evolved 180–220 million years ago and the present study highlights that it is still evolving and more genes can be added to the repertoire of this family. In addition, its core-genome includes important viral proteins including glycoprotein B and helicase. Most of the infections caused by human herpesviruses have no definitive cure; thus, search for new therapeutic strategies is necessary. The present study finds core-genome of human herpesviruses that differs from that of Herpesviridae family and nonhuman herpes strains of this family and might be a putative target for vaccine development. The phylogenetic reconstruction based upon the protein sequences of core gene set of Herpesviridae family reveals the sharp splits of its different subfamilies and supports the hypothesis of coevolution of viruses with their hosts. In addition, data mining for cis-elements in the genomes of human herpesviruses results in the prediction of numerous regulatory elements which can be used for regulating the expression of viral based vectors implicated in gene therapies.

Genome analysis of a clinical isolate of Shewanella sp. uncovered an active hybrid integrative and conjugative element carrying an integron platform inserted in a novel genomic locus

Shewanella spp. are currently considered to be emerging pathogens that can code for a blaOXA carbapenemase in their chromosome. Complete genome analysis of the clinical isolate Shewanella sp. Sh95 revealed that this strain is a novel species, which shares a lineage with marine isolates. Characterization of its resistome showed that it codes for genes drfA15, qacH and blaOXA-48. We propose that Shewanella sp. Sh95 acts as reservoir of blaOXA-48. Moreover, analysis of mobilome showed that it contains a novel integrative and conjugative element (ICE), named ICESh95. Comparative analysis between the close relatives ICESpuPO1 from Shewanella sp. W3-18-1 and ICE SXTMO10 from Vibrio cholerae showed that ICESh95 encompassed two new regions, a type III restriction modification system and a multidrug resistance integron. The integron platform contained a novel arrangement formed by gene cassettes drfA15 and qacH, and a class C-attC group II intron. Furthermore, insertion of ICESh95 occurred at a unique target site, which correlated with the presence of a different xis/int module. Mobility of ICESh95 was assessed and demonstrated its ability to self-transfer with high efficiency to different species of bacteria. Our results show that ICESh95 is a self-transmissible, mobile element, which can contribute to the dissemination of antimicrobial resistance; this is clearly a threat when natural bacteria from water ecosystems, such as Shewanella, act as vectors in its propagation.

The Geobacillus Pan-Genome: Implications for the Evolution of the Genus

The genus Geobacillus is comprised of a diverse group of spore-forming Gram-positive thermophilic bacterial species and is well known for both its ecological diversity and as a source of novel thermostable enzymes. Although the mechanisms underlying the thermophilicity of the organism and the thermostability of its macromolecules are reasonably well understood, relatively little is known of the evolutionary mechanisms, which underlie the structural and functional properties of members of this genus. In this study, we have compared 29 Geobacillus genomes, with a specific focus on the elements, which comprise the conserved core and flexible genomes. Based on comparisons of conserved core and flexible genomes, we present evidence of habitat delineation with specific Geobacillus genomes linked to specific niches. Our analysis revealed that Geobacillus and Anoxybacillus share a high proportion of genes. Moreover, the results strongly suggest that horizontal gene transfer is a major factor deriving the evolution of Geobacillus from Bacillus, with genetic contributions from other phylogenetically distant taxa.

Metagenomic Signatures of Bacterial Adaptation to Life in the Phyllosphere of a Salt-Secreting Desert Tree

The leaves of Tamarix aphylla, a globally distributed, salt-secreting desert tree, are dotted with alkaline droplets of high salinity. To successfully inhabit these organic carbon-rich droplets, bacteria need to be adapted to multiple stress factors, including high salinity, high alkalinity, high UV radiation, and periodic desiccation. To identify genes that are important for survival in this harsh habitat, microbial community DNA was extracted from the leaf surfaces of 10 Tamarix aphylla trees along a 350-km longitudinal gradient. Shotgun metagenomic sequencing, contig assembly, and binning yielded 17 genome bins, six of which were >80% complete. These genomic bins, representing three phyla (Proteobacteria,Bacteroidetes, and Firmicutes), were closely related to halophilic and alkaliphilic taxa isolated from aquatic and soil environments. Comparison of these genomic bins to the genomes of their closest relatives revealed functional traits characteristic of bacterial populations inhabiting the Tamarix phyllosphere, independent of their taxonomic affiliation. These functions, most notably light-sensing genes, are postulated to represent important adaptations toward colonization of this habitat.

IMPORTANCE

Plant leaves are an extensive and diverse microbial habitat, forming the main interface between solar energy and the terrestrial biosphere. There are hundreds of thousands of plant species in the world, exhibiting a wide range of morphologies, leaf surface chemistries, and ecological ranges. In order to understand the core adaptations of microorganisms to this habitat, it is important to diversify the type of leaves that are studied. This study provides an analysis of the genomic content of the most abundant bacterial inhabitants of the globally distributed, salt-secreting desert tree Tamarix aphylla Draft genomes of these bacteria were assembled, using the culture-independent technique of assembly and binning of metagenomic data. Analysis of the genomes reveals traits that are important for survival in this habitat, most notably, light-sensing and light utilization genes.

Towards plant pangenomics

As an increasing number of genome sequences become available for a wide range of species, there is a growing understanding that the genome of a single individual is insufficient to represent the gene diversity within a whole species. Many studies examine the sequence diversity within genes, and this allelic variation is an important source of phenotypic variation which can be selected for by man or nature. However, the significant gene presence/absence variation that has been observed within species and the impact of this variation on traits is only now being studied in detail. The sum of the genes for a species is termed the pangenome, and the determination and characterization of the pangenome is a requirement to understand variation within a species. In this review, we explore the current progress in pangenomics as well as methods and approaches for the characterization of pangenomes for a wide range of plant species.

Comparative Genomic Analysis Reveals a Diverse Repertoire of Genes Involved in Prokaryote-Eukaryote Interactions within the Pseudovibrio Genus

Strains of the Pseudovibrio genus have been detected worldwide, mainly as part of bacterial communities associated with marine invertebrates, particularly sponges. This recurrent association has been considered as an indication of a symbiotic relationship between these microbes and their host. Until recently, the availability of only two genomes, belonging to closely related strains, has limited the knowledge on the genomic and physiological features of the genus to a single phylogenetic lineage. Here we present 10 newly sequenced genomes of Pseudovibrio strains isolated from marine sponges from the west coast of Ireland, and including the other two publicly available genomes we performed an extensive comparative genomic analysis. Homogeneity was apparent in terms of both the orthologous genes and the metabolic features shared amongst the 12 strains. At the genomic level, a key physiological difference observed amongst the isolates was the presence only in strain P. axinellae AD2 of genes encoding proteins involved in assimilatory nitrate reduction, which was then proved experimentally. We then focused on studying those systems known to be involved in the interactions with eukaryotic and prokaryotic cells. This analysis revealed that the genus harbors a large diversity of toxin-like proteins, secretion systems and their potential effectors. Their distribution in the genus was not always consistent with the phylogenetic relationship of the strains. Finally, our analyses identified new genomic islands encoding potential toxin-immunity systems, previously unknown in the genus. Our analyses shed new light on the Pseudovibrio genus, indicating a large diversity of both metabolic features and systems for interacting with the host. The diversity in both distribution and abundance of these systems amongst the strains underlines how metabolically and phylogenetically similar bacteria may use different strategies to interact with the host and find a niche within its microbiota. Our data suggest the presence of a sponge-specific lineage of Pseudovibrio. The reduction in genome size and the loss of some systems potentially used to successfully enter the host, leads to the hypothesis that P. axinellae strain AD2 may be a lineage that presents an ancient association with the host and that may be vertically transmitted to the progeny.

Inside the Pan-genome - Methods and Software Overview

The number of genomes that have been deposited in databases has increased exponentially after the advent of Next-Generation Sequencing (NGS), which produces high-throughput sequence data; this circumstance has demanded the development of new bioinformatics software and the creation of new areas, such as comparative genomics. In comparative genomics, the genetic content of an organism is compared against other organisms, which helps in the prediction of gene function and coding region sequences, identification of evolutionary events and determination of phylogenetic relationships. However, expanding comparative genomics to a large number of related bacteria, we can infer their lifestyles, gene repertoires and minimal genome size. In this context, a powerful approach called Pan-genome has been initiated and developed. This approach involves the genomic comparison of different strains of the same species, or even genus. Its main goal is to establish the total number of non-redundant genes that are present in a determined dataset. Pan-genome consists of three parts: core genome; accessory or dispensable genome; and species-specific or strain-specific genes. Furthermore, pan-genome is considered to be “open” as long as new genes are added significantly to the total repertoire for each new additional genome and “closed” when the newly added genomes cannot be inferred to significantly increase the total repertoire of the genes. To perform all of the required calculations, a substantial amount of software has been developed, based on orthologous and paralogous gene identification.

Genomic Insights into a New Citrobacter koseri Strain Revealed Gene Exchanges with the Virulence-Associated Yersinia pestis pPCP1 Plasmid

The history of infectious diseases raised the plague as one of the most devastating for human beings. Far too often considered an ancient disease, the frequent resurgence of the plague has led to consider it as a reemerging disease in Madagascar, Algeria, Libya, and Congo. The genetic factors associated with the pathogenicity of Yersinia pestis, the causative agent of the plague, involve the acquisition of the pPCP1 plasmid that promotes host invasion through the expression of the virulence factor Pla. The surveillance of plague foci after the 2003 outbreak in Algeria resulted in a positive detection of the specific pla gene of Y. pestis in rodents. However, the phenotypic characterization of the isolate identified a Citrobacter koseri. The comparative genomics of our sequenced C. koseri URMITE genome revealed a mosaic gene structure resulting from the lifestyle of our isolate and provided evidence for gene exchanges with different enteric bacteria. The most striking was the acquisition of a continuous 2 kb genomic fragment containing the virulence factor Pla of the Y. pestis pPCP1 plasmid; however, the subcutaneous injection of the CKU strain in mice did not produce any pathogenic effect. Our findings demonstrate that fast molecular detection of plague using solely the pla gene is unsuitable and should rather require Y. pestis gene marker combinations. We also suggest that the evolutionary force that might govern the expression of pathogenicity can occur through the acquisition of virulence genes but could also require the loss or the inactivation of resident genes such as antivirulence genes.

Comparative Genomic Analysis of Sulfurospirillum cavolei MES Reconstructed from the Metagenome of an Electrosynthetic Microbiome

Sulfurospirillum spp. play an important role in sulfur and nitrogen cycling, and contain metabolic versatility that enables reduction of a wide range of electron acceptors, including thiosulfate, tetrathionate, polysulfide, nitrate, and nitrite. Here we describe the assembly of a Sulfurospirillum genome obtained from the metagenome of an electrosynthetic microbiome. The ubiquity and persistence of this organism in microbial electrosynthesis systems suggest it plays an important role in reactor stability and performance. Understanding why this organism is present and elucidating its genetic repertoire provide a genomic and ecological foundation for future studies where Sulfurospirillum are found, especially in electrode-associated communities. Metabolic comparisons and in-depth analysis of unique genes revealed potential ecological niche-specific capabilities within the Sulfurospirillum genus. The functional similarities common to all genomes, i.e., core genome, and unique gene clusters found only in a single genome were identified. Based upon 16S rRNA gene phylogenetic analysis and average nucleotide identity, the Sulfurospirillum draft genome was found to be most closely related to Sulfurospirillum cavolei. Characterization of the draft genome described herein provides pathway-specific details of the metabolic significance of the newly described Sulfurospirillum cavolei MES and, importantly, yields insight to the ecology of the genus as a whole. Comparison of eleven sequenced Sulfurospirillum genomes revealed a total of 6246 gene clusters in the pan-genome. Of the total gene clusters, 18.5% were shared among all eleven genomes and 50% were unique to a single genome. While most Sulfurospirillum spp. reduce nitrate to ammonium, five of the eleven Sulfurospirillum strains encode for a nitrous oxide reductase (nos) cluster with an atypical nitrous-oxide reductase, suggesting a utility for this genus in reduction of the nitrous oxide, and as a potential sink for this potent greenhouse gas.

Ecological Genomics of the Uncultivated Marine Roseobacter Lineage CHAB-I-5

Members of the marine Roseobacter clade are major participants in global carbon and sulfur cycles. While roseobacters are well represented in cultures, several abundant pelagic lineages, including SAG-O19, DC5-80-3, and NAC11-7, remain largely uncultivated and show evidence of genome streamlining. Here, we analyzed the partial genomes of three single cells affiliated with CHAB-I-5, another abundant but exclusively uncultivated Roseobacter lineage. Members of this lineage encode several metabolic potentials that are absent in streamlined genomes. Examples are quorum sensing and type VI secretion systems, which enable them to effectively interact with host and other bacteria. Further analysis of the CHAB-I-5 single-cell amplified genomes (SAGs) predicted that this lineage comprises members with relatively large genomes (4.1 to 4.4 Mbp) and a high fraction of noncoding DNA (10 to 12%), which is similar to what is observed in many cultured, nonstreamlined Roseobacter lineages. The four uncultured lineages, while exhibiting highly variable geographic distributions, together represent >60% of the global pelagic roseobacters. They are consistently enriched in genes encoding the capabilities of light harvesting, oxidation of "energy-rich" reduced sulfur compounds and methylated amines, uptake and catabolism of various carbohydrates and osmolytes, and consumption of abundant exudates from phytoplankton. These traits may define the global prevalence of the four lineages among marine bacterioplankton.

Phylogenetic Analyses of Shigella and Enteroinvasive Escherichia coli for the Identification of Molecular Epidemiological Markers: Whole-Genome Comparative Analysis Does Not Support Distinct Genera Designation

As a leading cause of bacterial dysentery, Shigella represents a significant threat to public health and food safety. Related, but often overlooked, enteroinvasive Escherichia coli (EIEC) can also cause dysentery. Current typing methods have limited ability to identify and differentiate between these pathogens despite the need for rapid and accurate identification of pathogens for clinical treatment and outbreak response. We present a comprehensive phylogeny of Shigella and EIEC using whole genome sequencing of 169 samples, constituting unparalleled strain diversity, and observe a lack of monophyly between Shigella and EIEC and among Shigella taxonomic groups. The evolutionary relationships in the phylogeny are supported by analyses of population structure and hierarchical clustering patterns of translated gene homolog abundance. Lastly, we identified a panel of 254 single nucleotide polymorphism (SNP) markers specific to each phylogenetic cluster for more accurate identification of Shigella and EIEC. Our findings show that Shigella and EIEC are not distinct evolutionary groups within the E. coli genus and, thus, EIEC as a group is not the ancestor to Shigella. The multiple analyses presented provide evidence for reconsidering the taxonomic placement of Shigella. The SNP markers offer more discriminatory power to molecular epidemiological typing methods involving these bacterial pathogens.

Complete genome sequence of the potato pathogen Ralstonia solanacearum UY031

Ralstonia solanacearum is the causative agent of bacterial wilt of potato. Ralstonia solanacearum strain UY031 belongs to the American phylotype IIB, sequevar 1, also classified as race 3 biovar 2. Here we report the completely sequenced genome of this strain, the first complete genome for phylotype IIB, sequevar 1, and the fourth for the R. solanacearum species complex. In addition to standard genome annotation, we have carried out a curated annotation of type III effector genes, an important pathogenicity-related class of genes for this organism. We identified 60 effector genes, and observed that this effector repertoire is distinct when compared to those from other phylotype IIB strains. Eleven of the effectors appear to be nonfunctional due to disruptive mutations. We also report a methylome analysis of this genome, the first for a R. solanacearum strain. This analysis helped us note the presence of a toxin gene within a region of probable phage origin, raising the hypothesis that this gene may play a role in this strain's virulence.

Comparative genomics of Fructobacillus spp. and Leuconostoc spp. reveals niche-specific evolution of Fructobacillus spp

Background

Fructobacillus spp. in fructose-rich niches belong to the family Leuconostocaceae. They were originally classified as Leuconostoc spp., but were later grouped into a novel genus, Fructobacillus, based on their phylogenetic position, morphology and specific biochemical characteristics. The unique characters, so called fructophilic characteristics, had not been reported in the group of lactic acid bacteria, suggesting unique evolution at the genome level. Here we studied four draft genome sequences of Fructobacillus spp. and compared their metabolic properties against those of Leuconostoc spp.

Results

Fructobacillus species possess significantly less protein coding sequences in their small genomes. The number of genes was significantly smaller in carbohydrate transport and metabolism. Several other metabolic pathways, including TCA cycle, ubiquinone and other terpenoid-quinone biosynthesis and phosphotransferase systems, were characterized as discriminative pathways between the two genera. The adhE gene for bifunctional acetaldehyde/alcohol dehydrogenase, and genes for subunits of the pyruvate dehydrogenase complex were absent in Fructobacillus spp. The two genera also show different levels of GC contents, which are mainly due to the different GC contents at the third codon position.

Conclusion

The present genome characteristics in Fructobacillus spp. suggest reductive evolution that took place to adapt to specific niches.

Electronic supplementary material

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

Comparative genomics of Pseudomonas fluorescens subclade III strains from human lungs

Background

While the taxonomy and genomics of environmental strains from the P. fluorescens species-complex has been reported, little is known about P. fluorescens strains from clinical samples. In this report, we provide the first genomic analysis of P. fluorescens strains in which human vs. environmental isolates are compared.

Results

Seven P. fluorescens strains were isolated from respiratory samples from cystic fibrosis (CF) patients. The clinical strains could grow at a higher temperature (>34 °C) than has been reported for environmental strains. Draft genomes were generated for all of the clinical strains, and multi-locus sequence analysis placed them within subclade III of the P. fluorescens species-complex. All strains encoded type- II, −III, −IV, and -VI secretion systems, as well as the widespread colonization island (WCI). This is the first description of a WCI in P. fluorescens strains. All strains also encoded a complete I2/PfiT locus and showed evidence of horizontal gene transfer. The clinical strains were found to differ from the environmental strains in the number of genes involved in metal resistance, which may be a possible adaptation to chronic antibiotic exposure in the CF lung.

Conclusions

This is the largest comparative genomics analysis of P. fluorescens subclade III strains to date and includes the first clinical isolates. At a global level, the clinical P. fluorescens subclade III strains were largely indistinguishable from environmental P. fluorescens subclade III strains, supporting the idea that identifying strains as ‘environmental’ vs ‘clinical’ is not a phenotypic trait. Rather, strains within P. fluorescens subclade III will colonize and persist in any niche that provides the requirements necessary for growth.

Electronic supplementary material

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

The pan-genome of Lactobacillus reuteri strains originating from the pig gastrointestinal tract

BACKGROUND

Lactobacillus reuteri is a gut symbiont of a wide variety of vertebrate species that has diversified into distinct phylogenetic clades which are to a large degree host-specific. Previous work demonstrated host specificity in mice and begun to determine the mechanisms by which gut colonisation and host restriction is achieved. However, how L. reuteri strains colonise the gastrointestinal (GI) tract of pigs is unknown.

RESULTS

To gain insight into the ecology of L. reuteri in the pig gut, the genome sequence of the porcine small intestinal isolate L. reuteri ATCC 53608 was completed and consisted of a chromosome of 1.94 Mbp and two plasmids of 138.5 kbp and 9.09 kbp, respectively. Furthermore, we generated draft genomes of four additional L. reuteri strains isolated from pig faeces or lower GI tract, lp167-67, pg-3b, 20-2 and 3c6, and subjected all five genomes to a comparative genomic analysis together with the previously completed genome of strain I5007. A phylogenetic analysis based on whole genomes showed that porcine L. reuteri strains fall into two distinct clades, as previously suggested by multi-locus sequence analysis. These six pig L. reuteri genomes contained a core set of 1364 orthologous gene clusters, as determined by OrthoMCL analysis, that contributed to a pan-genome totalling 3373 gene clusters. Genome comparisons of the six pig L. reuteri strains with 14 L. reuteri strains from other host origins gave a total pan-genome of 5225 gene clusters that included a core genome of 851 gene clusters but revealed that there were no pig-specific genes per se. However, genes specific for and conserved among strains of the two pig phylogenetic lineages were detected, some of which encoded cell surface proteins that could contribute to the diversification of the two lineages and their observed host specificity.

CONCLUSIONS

This study extends the phylogenetic analysis of L. reuteri strains at a genome-wide level, pointing to distinct evolutionary trajectories of porcine L. reuteri lineages, and providing new insights into the genomic events in L. reuteri that occurred during specialisation to their hosts. The occurrence of two distinct pig-derived clades may reflect differences in host genotype, environmental factors such as dietary components or to evolution from ancestral strains of human and rodent origin following contact with pig populations.

Divergence between the Highly Virulent Zoonotic Pathogen Helicobacter heilmannii and Its Closest Relative, the Low-Virulence "Helicobacter ailurogastricus" sp. nov

Helicobacter heilmannii naturally colonizes the stomachs of dogs and cats and has been associated with gastric disorders in humans. Nine feline Helicobacter strains, classified as H. heilmannii based on ureAB and 16S rRNA gene sequences, were divided into a highly virulent and a low-virulence group. The genomes of these strains were sequenced to investigate their phylogenetic relationships, to define their gene content and diversity, and to determine if the differences in pathogenicity were associated with the presence or absence of potential virulence genes. The capacities of these helicobacters to bind to the gastric mucosa were investigated as well. Our analyses revealed that the low-virulence strains do not belong to the species H. heilmannii but to a novel, closely related species for which we propose the name Helicobacter ailurogastricus. Several homologs of H. pylori virulence factors, such as IceA1, HrgA, and jhp0562-like glycosyltransferase, are present in H. heilmannii but absent in H. ailurogastricus. Both species contain a VacA-like autotransporter, for which the passenger domain is remarkably larger in H. ailurogastricus than in H. heilmannii. In addition, H. ailurogastricus shows clear differences in binding to the gastric mucosa compared to H. heilmannii. These findings highlight the low-virulence character of this novel Helicobacter species.

Positive selection is the main driving force for evolution of citrus canker-causing Xanthomonas

Understanding the evolutionary history and potential of bacterial pathogens is critical to prevent the emergence of new infectious bacterial diseases. Xanthomonas axonopodis subsp. citri (Xac) (synonym X. citri subsp. citri), which causes citrus canker, is one of the hardest-fought plant bacterial pathogens in US history. Here, we sequenced 21 Xac strains (14 XacA, 3 XacA* and 4 XacA(w)) with different host ranges from North America and Asia and conducted comparative genomic and evolutionary analyses. Our analyses suggest that acquisition of beneficial genes and loss of detrimental genes most likely allowed XacA to infect a broader range of hosts as compared with XacA(w) and XacA*. Recombination was found to have occurred frequently on the relative ancient branches, but rarely on the young branches of the clonal genealogy. The ratio of recombination/mutation ρ/θ was 0.0790±0.0005, implying that the Xac population was clonal in structure. Positive selection has affected 14% (395 out of 2822) of core genes of the citrus canker-causing Xanthomonas. The genes affected are enriched in 'carbohydrate transport and metabolism' and 'DNA replication, recombination and repair' genes (P<0.05). Many genes related to virulence, especially genes involved in the type III secretion system and effectors, are affected by positive selection, further highlighting the contribution of positive selection to the evolution of citrus canker-causing Xanthomonas. Our results suggest that both metabolism and virulence genes provide advantages to endow XacA with higher virulence and a wider host range. Our analysis advances our understanding of the genomic basis of specialization by positive selection in bacterial evolution.