Analyzing patterns of microbial evolution using the mauve genome alignment system

Unveiling the pathogenic and multidrug-resistant profiles of Vibrio alfacsensis: A potential identified threat in turbot (Scophthalmus maximus) aquaculture

Vibrio alfacsensis is traditionally seen as an environmental symbiont within its genus, with no detailedly documented pathogenicity in marine aquaculture to date. This study delves into the largely unexplored pathogenic potential and emerging antibiotic resistance of V. alfacsensis. The VA-1 strain, isolated from recirculating aquaculture system (RAS) effluent of cultured turbot (Scophthalmus maximus), underwent comprehensive analysis including biochemical identification, antibiotic susceptibility testing and reinfection trials. The results confirmed VA-1's pathogenicity and significant multiple antibiotic resistance. VA-1 could induce systemic infection in turbot, with symptoms like kidney enlargement, exhibiting virulence comparable to known Vibrio pathogens, with an LD50 around 2.36 × 106 CFU/fish. VA-1's remarkable resistance phenotype (14/22) suggested potential for genetic exchange and resistance factor acquisition in aquaculture environments. Phylogenetic analysis based on 16S rDNA sequences and whole-genome sequencing has firmly placed VA-1 within the V. alfacsensis clade, while genome-wide analysis highlights its similarity and diversity in relation to strains from across the globe. VA-1 contained numerous replicons, indicating the possibility for the spread of resistance and virulence genes. This study suggests V. alfacsensis may acquire and transfer pathogenic and resistant traits through horizontal gene transfer, a likelihood intensified by changing environmental and aquaculture conditions, highlighting the need for vigilant pathogen monitoring and new non-antibiotic treatments.

Comprehensive genomic analysis of Brevibacillus brevis BF19 reveals its biocontrol potential against bitter gourd wilt

Bitter gourd wilt, a severe vascular disease triggered by the soilborne pathogen Fusarium oxysporum f. sp. momordicae (FOM), markedly constrains bitter gourd yield. In this study, a novel strain BF19 of Brevibacillus brevis was isolated and identified, exhibiting strong antimicrobial activity against FOM through in vivo and in vitro experiments. To comprehensively assess the biocontrol potential of strain BF19, we conducted phenotypic, phylogenetic, and comparative genomics analyses. Phenotypic analysis revealed that BF19 exhibited 53.33% biocontrol efficacy and significantly increased the average plant height, root fresh weight, and dry weight. Whole-genome sequencing and comparative genomic analysis revealed numerous potential genes associated with biocontrol mechanisms in BF19. Importantly, the integration of metabolic cluster prediction with liquid chromatography‒tandem mass spectrometry (LC‒MS/MS) revealed the presence of a macrobrevin antibiotic, a product of polyketide synthases (PKSs), predominantly in BF19 fermentation products. The effectiveness of the Br. brevis strain BF19 and its crude extract against bitter gourd wilt has also been confirmed. This study provides a genetic framework for future investigations on PKSs and establishes a scientific basis for optimizing field applications of microbial biopesticides derived from Br. brevis BF19.

Genomic Analysis of Human-infecting Leptospira borgpetersenii isolates in Sri Lanka: expanded PF07598 gene family repertoire, less overall genome reduction than bovine isolates

Leptospira borgpetersenii commonly causes human leptospirosis, including severe disease. The first published analysis of L. borgpetersenii, performed on two strains of serovar Hardjo (L550 and JB197), concluded that the L. borgpetersenii genome is in the process of genome decay with functional consequences leading to a more obligately host-dependent life cycle. Yet whole genome analysis has only been carried out on few strains of L. borgpetersenii, with limited closed genomes and comprehensive analysis. Herein we report the complete, circularized genomes of seven non-Hardjo Leptospira borgpetersenii isolates from human leptospirosis patients in Sri Lanka. These isolates (all ST144) were found to be nearly identical by whole genome analysis; serotyping showed they are a novel serovar. We show that the L. borgpetersenii isolated from humans in Sri Lanka are less genomically decayed than previously reported isolates: fewer pseudogenes (N=141) and Insertion Sequence (IS) elements (N=46) compared to N=248, N=270, and N=400 pseudogenes, and N=121 and N=116 IS elements in published L. borgpetersenii Hardjo genomes (L550, JB197 and TC112). Compared to previously published L. borgpetersenii whole genome analyses showing two to three VM proteins in L. borgpetersenii isolates from cattle, rats and humans, we found that all of the human L. borgpetersenii isolates from Sri Lanka, including previously reported serovar Piyasena, have 4 encoded VM proteins, one ortholog of L. interrogans Copenhageni LIC12339 and 3 orthologs of LIC12844. Our findings of fewer pseudogenes, IS elements and expansion of the LIC12844 homologs of the PF07598 family in these human isolates suggests that this newly identified L. borgpetersenii serovar from Sri Lanka has unique pathogenicity. Comparative genome analysis and experimental studies of these L. borgpetersenii isolates will enable deeper insights into the molecular and cellular mechanisms of leptospirosis pathogenesis.

Halotolerant Bacillus sp. strain RA coordinates myo-inositol metabolism to confer salt tolerance to tomato

Soil salinity is a worldwide problem threatening crop yields. Some plant growth-promoting rhizobacteria (PGPR) could survive in high salt environment and assist plant adaptation to stress. Nevertheless, the genomic and metabolic features, as well as the regulatory mechanisms promoting salt tolerance in plants by these bacteria remain largely unknown. In the current work, a novel halotolerant PGPR strain, namely, Bacillus sp. strain RA can enhance tomato tolerance to salt stress. Comparative genomic analysis of strain RA with its closely related species indicated a high level of evolutionary plasticity exhibited by strain-specific genes and evolutionary constraints driven by purifying selection, which facilitated its genomic adaptation to salt-affected soils. The transcriptome further showed that strain RA could tolerate salt stress by balancing energy metabolism via the reprogramming of biosynthetic pathways. Plants exude a plethora of metabolites that can strongly influence plant fitness. The accumulation of myo-inositol in leaves under salt stress was observed, leading to the promotion of plant growth triggered by Bacillus sp. strain RA. Importantly, myo-inositol serves as a selective force in the assembly of the phyllosphere microbiome and the recruitment of plant-beneficial species. It promotes destabilizing properties in phyllosphere bacterial co-occurrence networks, but not in fungal networks. Furthermore, interdomain interactions between bacteria and fungi were strengthened by myo-inositol in response to salt stress. This work highlights the genetic adaptation of RA to salt-affected soils and its ability to impact phyllosphere microorganisms through the adjustment of myo-inositol metabolites, thereby imparting enduring resistance against salt stress in tomato.

Evolutionary analysis and biological characterization of a novel alphabaculovirus isolated from Mythimna separata

Baculoviruses are insect-specific pathogens. Novel baculovirus isolates provide new options for the biological control of pests. Therefore, research into the biological characteristics of newly isolated baculoviruses, including accurate classification and nomenclature, is important. In this study, a baculovirus was isolated from Mythimna separata and its complete genome sequence was determined by next-generation sequencing. The double-stranded DNA genome was 153 882 bp in length, encoding 163 open reading frames. The virus was identified as a variant of Mamestra brassicae multiple nucleopolyhedrovirus (MbMNPV) and designated Mamestra brassicae multiple nucleopolyhedrovirus CHN1 (MbMNPV-CHN1) according to ultrastructural analysis, genome comparison and phylogenetic analysis. Phylogenetic inference placed MbMNPV-CHN1 in a clade containing isolates of MacoNPV-A, MacoNPV-B and MbMNPV, which we have designated the Mb-McNPV group. The genomes of isolates in the Mb-McNPV group exhibited a high degree of collinearity with relatively minor differences in the content of annotated open reading frames. The development of codon usage bias in the Mb-McNPV group was affected mainly by natural selection. MbMNPV-CHN1 shows high infectivity against seven species of Lepidoptera. The yield of MbMNPV-CHN1 in the fourth- and fifth-instar M. separata larvae was 6.25×109-1.23×1010 OBs/cadaver. Our data provide insights into the classification, host range and virulence differences among baculoviruses of the Mb-McNPV group, as well as a promising potential new baculoviral insecticide.

Characterization and Genomic Analysis of the Naphthalene-Degrading Delftia tsuruhatensis ULwDis3 Isolated from Seawater

Strains of the genus Delftia are poorly studied microorganisms. In this work, the complete genome of the naphthalene-degrading Delftia tsuruhatensis strain ULwDis3 isolated from seawater of the Gulf of Finland of the Baltic Sea was assembled. For the first time, genes encoding naphthalene cleavage pathways via salicylate and gentisate were identified in a strain of the genus Delftia. The genes are part of one operon (nag genes). Three open reading frames (ORFs) were found in the genome of D. tsuruhatensis strain ULwDis3 that encode gentisate 1.2-dioxygenase. One of the ORFs is part of the nag operon. The physiological and biochemical characteristics of the strain ULwDis3 when cultured in mineral medium with naphthalene as the sole source of carbon and energy were also studied. It was found that after 22 h of growth, the strain stopped consuming naphthalene, and at the same time, naphthalene 1.2-dioxygenase and salicylate 5-hydroxylase activities were not detected. Later, a decrease in the number of living cells and the death of the culture were observed. Gentisate 1.2-dioxygenase activity was detected from the time of gentisate formation until culture death.

Comparative Genomics of Legionella pneumophila Isolates from the West Bank and Germany Support Molecular Epidemiology of Legionnaires' Disease

Legionella pneumophila is an environmental bacterium and clinical pathogen that causes many life-threating outbreaks of an atypical pneumonia called Legionnaires' disease (LD). Studies of this pathogen have focused mainly on Europe and the United States. A shortage in L. pneumophila data is clearly observed for developing countries. To reduce this knowledge gap, L. pneumophila isolates were studied in two widely different geographical areas, i.e., the West Bank and Germany. For this study, we sequenced and compared the whole genome of 38 clinical and environmental isolates of L. pneumophila covering different MLVA-8(12) genotypes in the two areas. Sequencing was conducted using the Illumina HiSeq 2500 platform. In addition, two isolates (A194 and H3) were sequenced using a Pacific Biosciences (PacBio) RSII platform to generate complete reference genomes from each of the geographical areas. Genome sequences from 55 L. pneumophila strains, including 17 reference strains, were aligned with the genome sequence of the closest strain (L. pneumophila strain Alcoy). A whole genome phylogeny based on single nucleotide polymorphisms (SNPs) was created using the ParSNP software v 1.0. The reference genomes obtained for isolates A194 and H3 consisted of circular chromosomes of 3,467,904 bp and 3,691,263 bp, respectively. An average of 36,418 SNPs (min. 8569, max. 70,708 SNPs) against our reference strain L. pneumophila str. Alcoy, and 2367 core-genes were identified among the fifty-five strains. An analysis of the genomic population structure by SNP comparison divided the fifty-five L. pneumophila strains into six branches. Individual isolates in sub-lineages in these branches differed by less than 120 SNPs if they had the same MLVA genotype and were isolated from the same location. A bioinformatics analysis identified the genomic islands (GIs) for horizontal gene transfer and mobile genetic elements, demonstrating that L. pneumophila showed high genome plasticity. Four L. pneumophila isolates (H3, A29, A129 and L10-091) contained well-defined plasmids. On average, only about half of the plasmid genes could be matched to proteins in databases. In silico phage findings suggested that 43 strains contained at least one phage. However, none of them were found to be complete. BLASTp analysis of proteins from the type IV secretion Dot/Icm system showed those proteins highly conserved, with less than 25% structural differences in the new L. pneumophila isolates. Overall, we demonstrated that whole genome sequencing provides a molecular surveillance tool for L. pneumophila at the highest conceivable discriminatory level, i.e., two to eight SNPs were observed for isolates from the same location but several years apart.

Pseudomonas veronii strain 7-41 degrading medium-chain n-alkanes and polycyclic aromatic hydrocarbons

Pollution of the environment by crude oil and oil products (represented by various types of compounds, mainly aliphatic, mono- and polyaromatic hydrocarbons) poses a global problem. The strain Pseudomonas veronii 7-41 can grow on medium-chain n-alkanes (C₈-C₁₂) and polycyclic aromatic hydrocarbons such as naphthalene. We performed a genetic analysis and physiological/biochemical characterization of strain 7-41 cultivated in a mineral medium with decane, naphthalene or a mixture of the hydrocarbons. The genes responsible for the degradation of alkanes and PAHs are on the IncP-7 conjugative plasmid and are organized into the alk and nah operons typical of pseudomonads. A natural plasmid carrying functional operons for the degradation of two different classes of hydrocarbons was first described. In monosubstrate systems, 28.4% and 68.8% of decane and naphthalene, respectively, were biodegraded by the late stationary growth phase. In a bisubstrate system, these parameters were 25.4% and 20.8% by the end of the exponential growth phase. Then the biodegradation stopped, and the bacterial culture started dying due to the accumulation of salicylate (naphthalene-degradation metabolite), which is toxic in high concentrations. The activity of the salicylate oxidation enzymes was below the detection limit. These results indicate that the presence of decane and a high concentration of salicylate lead to impairment of hydrocarbon degradation by the strain.

Genomic Analysis of Halotolerant Bacterial Strains Martelella soudanensis NC18T and NC20

Two novel, halotolerant strains of Martelella soudanensis, NC18^T and NC20, were isolated from deep subsurface sediment, deeply sequenced, and comparatively analyzed with related strains. Based on a phylogenetic analysis using 16S rRNA gene sequences, the two strains grouped with members of the genus Martelella. Here, we sequenced the complete genomes of NC18^T and NC20 to understand the mechanisms of their halotolerance. The genome sizes and G+C content of the strains were 6.1 Mb and 61.8 mol%, respectively. Moreover, NC18^T and NC20 were predicted to contain 5,849 and 5,830 genes, and 5,502 and 5,585 protein-coding genes, respectively. Both strains contain the identically predicted 6 rRNAs and 48 tRNAs. The harboring of halotolerant-associated genes revealed that strains NC18^T and NC20 might tolerate high salinity through the accumulation of potassium ions in a "salt-in" strategy induced by K⁺ uptake protein (kup) and the K⁺ transport system (trkAH and kdpFABC). These two strains also use the ectoine transport system (dctPQM), the glycine betaine transport system (proVWX), and glycine betaine uptake protein (opu) to accumulate "compatible solutes," such as ectoine and glycine betaine, to protect cells from salt stress. This study reveals the halotolerance mechanism of strains NC18^T and NC20 in high salt environments and suggests potential applications for these halotolerant and halophilic strains in environmental biotechnology.

Complete Genome and Molecular Characterization of a New Cyprinid Herpesvirus 2 (CyHV-2) SH-01 Strain Isolated from Cultured Crucian Carp

Cyprinid herpesvirus 2 (CyHV-2) is a causative factor of herpesviral hematopoietic necrosis (HVHN) in farmed crucian carp (Carassius carassius) and goldfish (Carassius auratus). In this study, we analyzed the genomic characteristics of a new strain, CyHV-2 SH-01, isolated during outbreaks in crucian carp at a local fish farm near Shanghai, China. CyHV-2 SH-01 exhibited a high sensitivity to goldfish and crucian carp in our previous research. The complete genome of SH-01 is 290,428 bp with 154 potential open reading frames (ORFs) and terminal repeat (TR) regions at both ends. Compared to the sequenced genomes of other CyHVs, Carassius auratus herpesvirus (CaHV) and Anguillid herpesvirus 1 (AngHV-1), several variations were found in SH-01, including nucleotide mutations, deletions, and insertions, as well as gene duplications, rearrangements, and horizontal transfers. Overall, the genome of SH-01 shares 99.60% of its identity with that of ST-J1. Genomic collinearity analysis showed that SH-01 has a high degree of collinearity with another three CyHV-2 isolates, and it is generally closely related to CaHV, CyHV-1, and CyHV-3, although it contains many differences in locally collinear blocks (LCBs). The lowest degree of collinearity was found with AngHV-1, despite some homologous LCBs, indicating that they are evolutionarily the most distantly related. The results provide new clues to better understand the CyHV-2 genome through sequencing and sequence mining.

Metabacillus dongyingensis sp. nov. Is Represented by the Plant Growth-Promoting Bacterium BY2G20 Isolated from Saline-Alkaline Soil and Enhances the Growth of Zea mays L. under Salt Stress

A novel plant growth-promoting rhizobacterium (PGPR), which was designated strain BY2G20, was isolated from saline-alkaline soil in Dongying, China. Strain BY2G20 can grow at a NaCl range from 0 to 7% and a pH range from 7 to 9 and can prevent the growth of the phytopathogen Ralstonia solanacearum. Based on its phenotypic and genomic characteristics and phylogenetic analysis, strain BY2G20 represents a novel species of the genus Metabacillus, for which the name Metabacillus dongyingensis sp. nov. is proposed. Comparative genomic analysis of strain BY2G20 with its closely related species exhibited a high level of evolutionary plasticity derived by horizontal gene transfer, which facilitated adaptative evolution. Different evolutionary constraints have operated on the diverse functions of BY2G20, with the gene adapted to saline-alkaline ecosystems experiencing functional constraints. We determined the genetic properties of saline-alkaline tolerance and plant growth promotion, such as cation-proton antiporters, cation transporters, osmoprotectant synthesis and transport, H⁺-transporting F₁F₀-ATPase, indole-3-acetic acid production, and secondary metabolite synthesis. We also evaluated the effects of strain BY2G20 on the growth of Zea mays L. (maize) under salt stress. The physiological parameters of maize such as plant height, stem diameter, dry biomass, and fresh biomass were significantly higher after inoculating strain BY2G20 under salt stress, indicating that inoculation with BY2G20 enhanced the growth of maize in saline areas. This study demonstrates that M. dongyingensis sp. nov. BY2G20 is a potential candidate for organic agriculture biofertilizers in saline-alkaline areas.

IMPORTANCE Plant growth and yield are adversely affected by soil salinity. PGPRs can promote plant growth and enhance plant tolerance to salt stress. In this study, a saline-alkaline tolerant PGPR strain BY2G20 was isolated from the rhizosphere of Ulmus pumila in Dongying, China. Strain BY2G20 represents a novel species within the genus Metabacillus based on phenotypic, genomic, and phylogenetic analysis. Genomic components have undergone different functional constraints, and the disparity in the evolutionary rate may be associated with the adaptation to a specific niche. Genomic analysis revealed numerous adaptive features of strain BY2G20 to a saline-alkaline environment and rhizosphere, especially genes related to salt tolerance, pH adaptability, and plant growth promotion. Our work also exhibited that inoculation of strain BY2G20 enhanced the growth of maize under salt stress. This study demonstrates that PGPRs play an important role in stimulating salt tolerance in plants and can be used as biofertilizers to enhance the growth of crops in saline-alkaline areas.

Pan-Genome Analysis of Delftia tsuruhatensis Reveals Important Traits Concerning the Genetic Diversity, Pathogenicity, and Biotechnological Properties of the Species

Delftia tsuruhatensis strains have long been known to promote plant growth and biological control. Recently, it has become an emerging opportunistic pathogen in humans. However, the genomic characteristics of the genetic diversity, pathogenicity, and biotechnological properties have not yet been comprehensively investigated. Here, a comparative pan-genome analysis was constructed. The open pan-genome with a large and flexible gene repertoire exhibited a high degree of genetic diversity. The purifying selection was the main force to drive pan-genome evolution. Significant differences were observed in the evolutionary relationship, functional enrichment, and degree of selective pressure between the different components of the pan-genome. A high degree of genetic plasticity was characterized by the determinations of diverse mobile genetic elements (MGEs), massive genomic rearrangement, and horizontal genes. Horizontal gene transfer (HGT) plays an important role in the genetic diversity of this bacterium and the formation of genomic traits. Our results revealed the occurrence of diverse virulence-related elements associated with macromolecular secretion systems, virulence factors associated with multiple nosocomial infections, and antimicrobial resistance, indicating the pathogenic potential. Lateral flagellum, T1SS, T2SS, T6SS, Tad pilus, type IV pilus, and a part of virulence-related genes exhibited general properties, whereas polar flagellum, T4SS, a part of virulence-related genes, and resistance genes presented heterogeneous properties. The pan-genome also harbors abundant genetic traits related to secondary metabolism, carbohydrate active enzymes (CAZymes), and phosphate transporter, indicating rhizosphere adaptation, plant growth promotion, and great potential uses in agriculture and biological control. This study provides comprehensive insights into this uncommon species from the genomic perspective.

IMPORTANCED. tsuruhatensis is considered a plant growth-promoting rhizobacterium (PGPR), an organic pollutant degradation strain, and an emerging opportunistic pathogen to the human. However, the genetic diversity, the evolutionary dynamics, and the genetic basis of these remarkable traits are still little known. We constructed a pan-genome analysis for D. tsuruhatensis and revealed extensive genetic diversity and genetic plasticity exhibited by open pan-genome, diverse mobile genetic elements (MGEs), genomic rearrangement, and horizontal genes. Our results highlight that horizontal gene transfer (HGT) and purifying selection are important forces in D. tsuruhatensis genetic evolution. The abundant virulence-related elements associated with macromolecular secretion systems, virulence factors, and antimicrobial resistance could contribute to the pathogenicity of this bacterium. Therefore, clinical microbiologists need to be aware of D. tsuruhatensis as an opportunistic pathogen. The genetic profiles of secondary metabolism, carbohydrate active enzymes (CAZymes), and phosphate transporter could provide insight into the genetic armory of potential applications for agriculture and biological control of D. tsuruhatensis in general.

The Evolutionary Trend and Genomic Features of an Emerging Lineage of Elizabethkingia anophelis Strains in Taiwan

The incidence of Elizabethkingia anophelis bacteremia increased significantly in a tertiary hospital, Changhua Christian Hospital (CCH) since 2013. The infection density was 1.3 and 8.1 cases per 100,000 patient-days between 2005 and 2012 and 2013 and 2020, respectively (P < 0.05). During an outbreak investigation, a specific lineage of E. anophelis strains was identified by the pulsed-field gel electrophoresis analysis. To evaluate the evolution of the specific E. anophelis lineage, whole-genome sequencing was performed, and unique genomic features (GRs) were determined by comparative genomic analysis. The specific E. anophelis lineage was novel compared to worldwide strains ever reported by cg-MLST phylogenic and whole-genome comparative analysis. Multiplex PCR using primers designed from unique GRs were performed for prevalence screening among isolates from the CCH and nationwide isolates from the Taiwan surveillance of Antimicrobial Resistance (TSAR) Program. The proportion of the specific E. anophelis lineage increased from 7.9% (3/38) during 2005-2012 to 89.2% (223/250) during 2013-2020 (P < 0.05). Although E. anophelis usually confers resistance to multiple antibiotics with limited therapeutic options, the E. anophelis strains in the specific lineage had higher ciprofloxacin resistance (100% [226/226] versus 27.4% [17/62], P < 0.05) and was associated with a higher 14-day mortality rates (33.2% [37/226] versus 16.1% [10/62], P < 0.05) than other strains at CCH. A similarly increasing trend was also found in the national TSAR program during 2002-2018 (p for trend <0.05). We concluded that a novel lineage of E. anophelis strains has emerged dominantly in Taiwan. The genomic features are important for further investigations of epidemiology, resistance, virulence, and appropriate treatment.

IMPORTANCEElizabethkingia anophelis is an emerging multidrug resistant pathogen caused several global outbreaks recently. E. anophelis was frequently misidentified as E. meningoseptica in the past by conventional culture methods; therefore, the prevalence was often underestimated. Through revised identification, an increasing trend of E. anophelis infection was noted in a tertiary hospital and a dominant lineage of strains was recognized by genotyping. To our best knowledge, the dominant lineage of E. anophelis is novel in comparison to other worldwide strains by whole-genome comparative analysis and several unique genomic regions were found. The whole-genome sequencing data also demonstrated multiple putative virulence factors and genes associated with multidrug resistance. In our study, we identified a specially evolved E. anophelis in Taiwan with increasing nationwide dominance. This study will assist in further epidemiology surveillance and developing corresponsive infection control policies to restrain it potential of global dissemination.

Changing Molecular Epidemiology of Vibrio cholerae Outbreaks in Shanghai, China

The 7th cholera pandemic began in 1961 in Sulawesi, Indonesia, and then spread around the world in at least three waves. However, the lack of genome sequences for Vibrio cholerae strains under long-term surveillance in East Asia, especially in China, has restricted our understanding of the dynamics of the intracountry and intercountry evolution and transmission of the 7th-pandemic clones. In this study, we obtained the genome sequences of 60 V. cholerae strains isolated in Shanghai, the largest port in the world and the largest city in China, from 1961 to 2011. Our whole-genome-based phylogeny of 7th-pandemic strains revealed that all but one fell into five "stages," most of which are single clades and share independent ancestors. Each stage dominated in succession for a period, with little overlap between them. In addition, two near-identical Shanghai strains belonging to a pre-7th-pandemic precursor and 4 nontoxigenic O1/O139 strains attributed to independent recombination events at the O-antigen loci were present. The major lineages of the 7th pandemic in Shanghai appeared to be closely related to V. cholerae strains isolated from South or Southeast Asia. Stage succession was consistently related to changes in society and human activity, implying that human-caused niche change may play a vital role in the cholera dynamics in Shanghai.IMPORTANCE V. cholerae is the causative agent of cholera, a life-threatening disease characterized by severe, watery diarrhea. The 7th pandemic started in Indonesia in 1961 and spread globally, currently infecting 1.3 million to 4 million people annually. Here, we applied whole-genome sequencing to analyze a long-term collection of V. cholerae clinical strains to reveal the phylogenetic background and evolutionary dynamics of the 7th pandemic in Shanghai, which had undergone breathtakingly rapid development in the last half-century. All but one of the Shanghai 7th-pandemic strains fell into five "stages" that were dominant in Shanghai and appeared to be closely related to 7th-pandemic strains of South or Southeast Asia. Our findings extended the understanding of the dynamics of the evolution and transmission of the 7th-pandemic clones in East Asia and the relationship between social changes and cholera epidemiology.

Comparative genomics of Salmonella enterica serovar Montevideo reveals lineage-specific gene differences that may influence ecological niche association

Salmonella enterica serovar Montevideo has been linked to recent foodborne illness outbreaks resulting from contamination of products such as fruits, vegetables, seeds and spices. Studies have shown that Montevideo also is frequently associated with healthy cattle and can be isolated from ground beef, yet human salmonellosis outbreaks of Montevideo associated with ground beef contamination are rare. This disparity fuelled our interest in characterizing the genomic differences between Montevideo strains isolated from healthy cattle and beef products, and those isolated from human patients and outbreak sources. To that end, we sequenced 13 Montevideo strains to completion, producing high-quality genome assemblies of isolates from human patients (n=8) or from healthy cattle at slaughter (n=5). Comparative analysis of sequence data from this study and publicly available sequences (n=72) shows that Montevideo falls into four previously established clades, differentially occupied by cattle and human strains. The results of these analyses reveal differences in metabolic islands, environmental adhesion determinants and virulence factors within each clade, and suggest explanations for the infrequent association between bovine isolates and human illnesses.

A Tad pilus promotes the establishment and resistance of Vibrio vulnificus biofilms to mechanical clearance

Vibrio vulnificus is autochthonous to estuaries and warm coastal waters. Infection occurs via open wounds or ingestion, where its asymptomatic colonization of seafood, most infamously oysters, provides a gateway into the human food chain. Colonization begins with initial surface contact, which is often mediated by bacterial surface appendages called pili. Type IV Tad pili are widely distributed in the Vibrionaceae, but evidence for a physiological role for these structures is scant. The V. vulnificus genome codes for three distinct tad loci. Recently, a positive correlation was demonstrated between the expression of tad-3 and the phenotypes of a V. vulnificus descendent (NT) that exhibited increased biofilm formation, auto-aggregation, and oyster colonization relative to its parent. However, the mechanism by which tad pilus expression promoted these phenotypes was not determined. Here, we show that deletion of the tad pilin gene (flp) altered the near-surface motility profile of NT cells from high curvature, orbital retracing patterns characteristic of cells actively probing the surface to low curvature traces indicative of wandering and diminished bacteria-surface interactions. The NT flp pilin mutant also exhibited decreased initial surface attachment, attenuated auto-aggregation and formed fragile biofilms that disintegrated under hydrodynamic flow. Thus, the tad-3 locus, designated iam, promoted initial surface attachment, auto-aggregation and resistance to mechanical clearance of V. vulnificus biofilms. The prevalence of tad loci in the Vibrionaceae suggests that they may play equally important roles in other family members.

A Novel Aerobic Degradation Pathway for Thiobencarb Is Initiated by the TmoAB Two-Component Flavin Mononucleotide-Dependent Monooxygenase System in Acidovorax sp. Strain T1

Thiobencarb is a thiocarbamate herbicide used in rice paddies worldwide. Microbial degradation plays a crucial role in the dissipation of thiobencarb in the environment. However, the physiological and genetic mechanisms underlying thiobencarb degradation remain unknown. In this study, a novel thiobencarb degradation pathway was proposed in Acidovorax sp. strain T1. Thiobencarb was oxidized and cleaved at the C-S bond, generating diethylcarbamothioic S-acid and 4-chlorobenzaldehyde (4CDA). 4CDA was then oxidized to 4-chlorobenzoic acid (4CBA) and hydrolytically dechlorinated to 4-hydroxybenzoic acid (4HBA). The identification of catabolic genes suggested further hydroxylation to protocatechuic acid (PCA) and finally degradation through the protocatechuate 4,5-dioxygenase pathway. A novel two-component monooxygenase system identified in the strain, TmoAB, was responsible for the initial catabolic reaction. TmoA shared 28 to 32% identity with the oxygenase components of pyrimidine monooxygenase from Agrobacterium fabrum, alkanesulfonate monooxygenase from Pseudomonas savastanoi, and dibenzothiophene monooxygenase from Rhodococcus sp. TmoB shared 25 to 37% identity with reported flavin reductases and oxidized NADH but not NADPH. TmoAB is a flavin mononucleotide (FMN)-dependent monooxygenase and catalyzed the C-S bond cleavage of thiobencarb. Introduction of tmoAB into cells of the thiobencarb degradation-deficient mutant T1m restored its ability to degrade and utilize thiobencarb. A dehydrogenase gene, tmoC, was located 7,129 bp downstream of tmoAB, and its transcription was clearly induced by thiobencarb. The purified TmoC catalyzed the dehydrogenation of 4CDA to 4CBA using NAD⁺ as a cofactor. A gene cluster responsible for the complete 4CBA metabolic pathway was also cloned, and its involvement in thiobencarb degradation was preliminarily verified by transcriptional analysis.IMPORTANCE Microbial degradation is the main factor in thiobencarb dissipation in soil. In previous studies, thiobencarb was degraded initially via N-deethylation, sulfoxidation, hydroxylation, and dechlorination. However, enzymes and genes involved in the microbial degradation of thiobencarb have not been studied. This study revealed a new thiobencarb degradation pathway in Acidovorax sp. strain T1 and identified a novel two-component FMN-dependent monooxygenase system, TmoAB. Under TmoAB-mediated catalysis, thiobencarb was cleaved at the C-S bond, producing diethylcarbamothioic S-acid and 4CDA. Furthermore, the downstream degradation pathway of thiobencarb was proposed. Our study provides the physiological, biochemical, and genetic foundation of thiobencarb degradation in this microorganism.

PacBio But Not Illumina Technology Can Achieve Fast, Accurate and Complete Closure of the High GC, Complex Burkholderia pseudomallei Two-Chromosome Genome

Although PacBio third-generation sequencers have improved the read lengths of genome sequencing which facilitates the assembly of complete genomes, no study has reported success in using PacBio data alone to completely sequence a two-chromosome bacterial genome from a single library in a single run. Previous studies using earlier versions of sequencing chemistries have at most been able to finish bacterial genomes containing only one chromosome with de novo assembly. In this study, we compared the robustness of PacBio RS II, using one SMRT cell and the latest P6-C4 chemistry, with Illumina HiSeq 1500 in sequencing the genome of Burkholderia pseudomallei, a bacterium which contains two large circular chromosomes, very high G+C content of 68–69%, highly repetitive regions and substantial genomic diversity, and represents one of the largest and most complex bacterial genomes sequenced, using a reference genome generated by hybrid assembly using PacBio and Illumina datasets with subsequent manual validation. Results showed that PacBio data with de novo assembly, but not Illumina, was able to completely sequence the B. pseudomallei genome without any gaps or mis-assemblies. The two large contigs of the PacBio assembly aligned unambiguously to the reference genome, sharing >99.9% nucleotide identities. Conversely, Illumina data assembled using three different assemblers resulted in fragmented assemblies (201–366 contigs), sharing only 92.2–100% and 92.0–100% nucleotide identities to chromosomes I and II reference sequences, respectively, with no indication that the B. pseudomallei genome consisted of two chromosomes with four copies of ribosomal operons. Among all assemblies, the PacBio assembly recovered the highest number of core and virulence proteins, and housekeeping genes based on whole-genome multilocus sequence typing (wgMLST). Most notably, assembly solely based on PacBio outperformed even hybrid assembly using both PacBio and Illumina datasets. Hybrid approach generated only 74 contigs, while the PacBio data alone with de novo assembly achieved complete closure of the two-chromosome B. pseudomallei genome without additional costly bench work and further sequencing. PacBio RS II using P6-C4 chemistry is highly robust and cost-effective and should be the platform of choice in sequencing bacterial genomes, particularly for those that are well-known to be difficult-to-sequence.

Draft genome sequences of strains Salinicola socius SMB35T, Salinicola sp. MH3R3-1 and Chromohalobacter sp. SMB17 from the Verkhnekamsk potash mining region of Russia

Halomonads are moderately halophilic bacteria that are studied as models of prokaryotic osmoadaptation and sources of enzymes and chemicals for biotechnological applications. Despite the progress in understanding the diversity of these organisms, our ability to explain ecological, metabolic, and biochemical traits of halomonads at the genomic sequence level remains limited. This study addresses this gap by presenting draft genomes of Salinicola socius SMB35^T, Salinicola sp. MH3R3–1 and Chromohalobacter sp. SMB17, which were isolated from potash mine tailings in the Verkhnekamsk salt deposit area of Russia. The analysis of these genomes confirmed the importance of ectoines and quaternary amines to the capacity of halomonads to tolerate osmotic stress and adapt to hypersaline environments. The study also revealed that Chromohalobacter and Salinicola share 75–90% of the predicted proteome, but also harbor a set of genus-specific genes, which in Salinicola amounted to approximately 0.5 Mbp. These genus-specific genome segments may contribute to the phenotypic diversity of the Halomonadaceae and the ability of these organisms to adapt to changing environmental conditions and colonize new ecological niches.

Origins of the current seventh cholera pandemic

Vibrio cholerae has caused seven cholera pandemics since 1817, imposing terror on much of the world, but bacterial strains are currently only available for the sixth and seventh pandemics. The El Tor biotype seventh pandemic began in 1961 in Indonesia, but did not originate directly from the classical biotype sixth-pandemic strain. Previous studies focused mainly on the spread of the seventh pandemic after 1970. Here, we analyze in unprecedented detail the origin, evolution, and transition to pandemicity of the seventh-pandemic strain. We used high-resolution comparative genomic analysis of strains collected from 1930 to 1964, covering the evolution from the first available El Tor biotype strain to the start of the seventh pandemic. We define six stages leading to the pandemic strain and reveal all key events. The seventh pandemic originated from a nonpathogenic strain in the Middle East, first observed in 1897. It subsequently underwent explosive diversification, including the spawning of the pandemic lineage. This rapid diversification suggests that, when first observed, the strain had only recently arrived in the Middle East, possibly from the Asian homeland of cholera. The lineage migrated to Makassar, Indonesia, where it gained the important virulence-associated elements Vibrio seventh pandemic island I (VSP-I), VSP-II, and El Tor type cholera toxin prophage by 1954, and it then became pandemic in 1961 after only 12 additional mutations. Our data indicate that specific niches in the Middle East and Makassar were important in generating the pandemic strain by providing gene sources and the driving forces for genetic events.

Comparative genomic analysis of toxin-negative strains of Clostridium difficile from humans and animals with symptoms of gastrointestinal disease

Background

Clostridium difficile infections (CDI) are a significant health problem to humans and food animals. Clostridial toxins ToxA and ToxB encoded by genes tcdA and tcdB are located on a pathogenicity locus known as the PaLoc and are the major virulence factors of C. difficile. While toxin-negative strains of C. difficile are often isolated from faeces of animals and patients suffering from CDI, they are not considered to play a role in disease. Toxin-negative strains of C. difficile have been used successfully to treat recurring CDI but their propensity to acquire the PaLoc via lateral gene transfer and express clinically relevant levels of toxins has reinforced the need to characterise them genetically. In addition, further studies that examine the pathogenic potential of toxin-negative strains of C. difficile and the frequency by which toxin-negative strains may acquire the PaLoc are needed.

Results

We undertook a comparative genomic analysis of five Australian toxin-negative isolates of C. difficile that lack tcdA, tcdB and both binary toxin genes cdtA and cdtB that were recovered from humans and farm animals with symptoms of gastrointestinal disease. Our analyses show that the five C. difficile isolates cluster closely with virulent toxigenic strains of C. difficile belonging to the same sequence type (ST) and have virulence gene profiles akin to those in toxigenic strains. Furthermore, phage acquisition appears to have played a key role in the evolution of C. difficile.

Conclusions

Our results are consistent with the C. difficile global population structure comprising six clades each containing both toxin-positive and toxin-negative strains. Our data also suggests that toxin-negative strains of C. difficile encode a repertoire of putative virulence factors that are similar to those found in toxigenic strains of C. difficile, raising the possibility that acquisition of PaLoc by toxin-negative strains poses a threat to human health. Studies in appropriate animal models are needed to examine the pathogenic potential of toxin-negative strains of C. difficile and to determine the frequency by which toxin-negative strains may acquire the PaLoc.

Electronic supplementary material

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

Recombination Does Not Hinder Formation or Detection of Ecological Species of Synechococcus Inhabiting a Hot Spring Cyanobacterial Mat

Recent studies of bacterial speciation have claimed to support the biological species concept—that reduced recombination is required for bacterial populations to diverge into species. This conclusion has been reached from the discovery that ecologically distinct clades show lower rates of recombination than that which occurs among closest relatives. However, these previous studies did not attempt to determine whether the more-rapidly recombining close relatives within the clades studied may also have diversified ecologically, without benefit of sexual isolation. Here we have measured the impact of recombination on ecological diversification within and between two ecologically distinct clades (A and B') of Synechococcus in a hot spring microbial mat in Yellowstone National Park, using a cultivation-free, multi-locus approach. Bacterial artificial chromosome (BAC) libraries were constructed from mat samples collected at 60°C and 65°C. Analysis of multiple linked loci near Synechococcus 16S rRNA genes showed little evidence of recombination between the A and B' lineages, but a record of recombination was apparent within each lineage. Recombination and mutation rates within each lineage were of similar magnitude, but recombination had a somewhat greater impact on sequence diversity than mutation, as also seen in many other bacteria and archaea. Despite recombination within the A and B' lineages, there was evidence of ecological diversification within each lineage. The algorithm Ecotype Simulation identified sequence clusters consistent with ecologically distinct populations (ecotypes), and several hypothesized ecotypes were distinct in their habitat associations and in their adaptations to different microenvironments. We conclude that sexual isolation is more likely to follow ecological divergence than to precede it. Thus, an ecology-based model of speciation appears more appropriate than the biological species concept for bacterial and archaeal diversification.

Whole-genome sequencing reveals the effect of vaccination on the evolution of Bordetella pertussis

Herd immunity can potentially induce a change of circulating viruses. However, it remains largely unknown that how bacterial pathogens adapt to vaccination. In this study, Bordetella pertussis, the causative agent of whooping cough, was selected as an example to explore possible effect of vaccination on the bacterial pathogen. We sequenced and analysed the complete genomes of 40 B. pertussis strains from Finland and China, as well as 11 previously sequenced strains from the Netherlands, where different vaccination strategies have been used over the past 50 years. The results showed that the molecular clock moved at different rates in these countries and in distinct periods, which suggested that evolution of the B. pertussis population was closely associated with the country vaccination coverage. Comparative whole-genome analyses indicated that evolution in this human-restricted pathogen was mainly characterised by ongoing genetic shift and gene loss. Furthermore, 116 SNPs were specifically detected in currently circulating ptxP3-containing strains. The finding might explain the successful emergence of this lineage and its spread worldwide. Collectively, our results suggest that the immune pressure of vaccination is one major driving force for the evolution of B. pertussis, which facilitates further exploration of the pathogenicity of B. pertussis.

The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes

Whole-genome sequences are now available for many microbial species and clades, however existing whole-genome alignment methods are limited in their ability to perform sequence comparisons of multiple sequences simultaneously. Here we present the Harvest suite of core-genome alignment and visualization tools for the rapid and simultaneous analysis of thousands of intraspecific microbial strains. Harvest includes Parsnp, a fast core-genome multi-aligner, and Gingr, a dynamic visual platform. Together they provide interactive core-genome alignments, variant calls, recombination detection, and phylogenetic trees. Using simulated and real data we demonstrate that our approach exhibits unrivaled speed while maintaining the accuracy of existing methods. The Harvest suite is open-source and freely available from: http://github.com/marbl/harvest.

Genome Sequence of the Heteropolysaccharide-Producing Strain Lactobacillus mucosae DPC 6426

Exopolysaccharide-synthesizing Lactobacillus mucosae DPC 6426 is a heterofermentative strain, which has demonstrated cholesterol-lowering properties in an animal model of lipid-driven atherosclerosis. The genome revealed a plethora of homologues linked to carbohydrate metabolism and mucin binding.

The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes

Whole-genome sequences are now available for many microbial species and clades, however existing whole-genome alignment methods are limited in their ability to perform sequence comparisons of multiple sequences simultaneously. Here we present the Harvest suite of core-genome alignment and visualization tools for the rapid and simultaneous analysis of thousands of intraspecific microbial strains. Harvest includes Parsnp, a fast core-genome multi-aligner, and Gingr, a dynamic visual platform. Together they provide interactive core-genome alignments, variant calls, recombination detection, and phylogenetic trees. Using simulated and real data we demonstrate that our approach exhibits unrivaled speed while maintaining the accuracy of existing methods. The Harvest suite is open-source and freely available from: http://github.com/marbl/harvest.

Whole genome phylogeny of Prochlorococcus marinus group of cyanobacteria: genome alignment and overlapping gene approach

Prochlorococcus is the smallest known oxygenic phototrophic marine cyanobacterium dominating the mid-latitude oceans. Physiologically and genetically distinct P. marinus isolates from many oceans in the world were assigned two different groups, a tightly clustered high-light (HL)-adapted and a divergent low-light (LL-) adapted clade. Phylogenetic analysis of this cyanobacterium on the basis of 16S rRNA and other conserved genes did not show consistency with its phenotypic behavior. We analyzed phylogeny of this genus on the basis of complete genome sequences through genome alignment, overlapping-gene content and gene-order approach. Phylogenetic tree of P. marinus obtained by comparing whole genome sequences in contrast to that based on 16S rRNA gene, corresponded well with the HL/LL ecotypic distinction of twelve strains and showed consistency with phenotypic classification of P. marinus. Evidence for the horizontal descent and acquisition of genes within and across the genus was observed. Many genes involved in metabolic functions were found to be conserved across these genomes and many were continuously gained by different strains as per their needs during the course of their evolution. Consistency in the physiological and genetic phylogeny based on whole genome sequence is established. These observations improve our understanding about the adaptation and diversification of these organisms under evolutionary pressure.

Less is more: Burkholderia pseudomallei and chronic melioidosis

The Gram-negative bacterium Burkholderia pseudomallei is the causative agent of melioidosis, a serious infectious disease of humans and animals. Once considered an esoteric tropical disease confined to Southeast Asia and northern Australia, research on B. pseudomallei has recently gained global prominence due to its classification as a potential bioterrorism agent by countries such as the United States and also by increasing numbers of case reports from regions where it is not endemic. An environmental bacterium typically found in soil and water, assessing the true global prevalence of melioidosis is challenged by the fact that clinical symptoms associated with B. pseudomallei infection are extremely varied and may be confused with diverse conditions such as lung cancer, tuberculosis, or Staphyloccocus aureus infection. These diagnostic challenges, coupled with lack of awareness among clinicians, have likely contributed to underdiagnosis and the high mortality rate of melioidosis, as initial treatment is often either inappropriate or delayed. Even after antibiotic treatment, relapses are frequent, and after resolution of acute symptoms, chronic melioidosis can also occur, and the symptoms can persist for months to years. In a recent article, Price et al. [mBio 4(4):e00388-13, 2013, doi:10.1128/mBio.00388-13] demonstrate how comparative genomic sequencing can reveal the repertoire of genetic changes incurred by B. pseudomallei during chronic human infection. Their results have significant clinical ramifications and highlight B. pseudomallei’s ability to survive in a wide range of potential niches within hosts, through the acquisition of genetic adaptations that optimize fitness and resource utilization.

Pathoadaptive mutations in Salmonella enterica isolated after serial passage in mice

How pathogenic bacteria adapt and evolve in the complex and variable environment of the host remains a largely unresolved question. Here we have used whole genome sequencing of Salmonella enterica serovar Typhimurium LT2 populations serially passaged in mice to identify mutations that adapt bacteria to systemic growth in mice. We found unique pathoadaptive mutations in two global regulators, phoQ and stpA, which increase the competitive indexes of the bacteria 3- to 5-fold. Also, all mouse-adapted lineages had changed the orientation of the hin invertable element, resulting in production of a FliC type of flagellum. Competition experiments in mice with locked flagellum mutants showed that strains expressing the FliC type of flagellum had a 5-fold increase in competitive index as compared to those expressing FljB type flagellum. Combination of the flagellum cassette inversion with the stpA mutation increased competitive indexes up to 20-fold. These experiments show that Salmonella can rapidly adapt to a mouse environment by acquiring a few mutations of moderate individual effect that when combined confer substantial increases in growth.

Impact of homologous and non-homologous recombination in the genomic evolution of Escherichia coli

Background

Escherichia coli is an important species of bacteria that can live as a harmless inhabitant of the guts of many animals, as a pathogen causing life-threatening conditions or freely in the non-host environment. This diversity of lifestyles has made it a particular focus of interest for studies of genetic variation, mainly with the aim to understand how a commensal can become a deadly pathogen. Many whole genomes of E. coli have been fully sequenced in the past few years, which offer helpful data to help understand how this important species evolved.

Results

We compared 27 whole genomes encompassing four phylogroups of Escherichia coli (A, B1, B2 and E). From the core-genome we established the clonal relationships between the isolates as well as the role played by homologous recombination during their evolution from a common ancestor. We found strong evidence for sexual isolation between three lineages (A+B1, B2, E), which could be explained by the ecological structuring of E. coli and may represent on-going speciation. We identified three hotspots of homologous recombination, one of which had not been previously described and contains the aroC gene, involved in the essential shikimate metabolic pathway. We also described the role played by non-homologous recombination in the pan-genome, and showed that this process was highly heterogeneous. Our analyses revealed in particular that the genomes of three enterohaemorrhagic (EHEC) strains within phylogroup B1 have converged from originally separate backgrounds as a result of both homologous and non-homologous recombination.

Conclusions

Recombination is an important force shaping the genomic evolution and diversification of E. coli, both by replacing fragments of genes with an homologous sequence and also by introducing new genes. In this study, several non-random patterns of these events were identified which correlated with important changes in the lifestyle of the bacteria, and therefore provide additional evidence to explain the relationship between genomic variation and ecological adaptation.

A comparative genomics perspective on the genetic content of the alkaliphilic haloarchaeon Natrialba magadii ATCC 43099T

BACKGROUND

Natrialba magadii is an aerobic chemoorganotrophic member of the Euryarchaeota and is a dual extremophile requiring alkaline conditions and hypersalinity for optimal growth. The genome sequence of Nab. magadii type strain ATCC 43099 was deciphered to obtain a comprehensive insight into the genetic content of this haloarchaeon and to understand the basis of some of the cellular functions necessary for its survival.

RESULTS

The genome of Nab. magadii consists of four replicons with a total sequence of 4,443,643 bp and encodes 4,212 putative proteins, some of which contain peptide repeats of various lengths. Comparative genome analyses facilitated the identification of genes encoding putative proteins involved in adaptation to hypersalinity, stress response, glycosylation, and polysaccharide biosynthesis. A proton-driven ATP synthase and a variety of putative cytochromes and other proteins supporting aerobic respiration and electron transfer were encoded by one or more of Nab. magadii replicons. The genome encodes a number of putative proteases/peptidases as well as protein secretion functions. Genes encoding putative transcriptional regulators, basal transcription factors, signal perception/transduction proteins, and chemotaxis/phototaxis proteins were abundant in the genome. Pathways for the biosynthesis of thiamine, riboflavin, heme, cobalamin, coenzyme F420 and other essential co-factors were deduced by in depth sequence analyses. However, approximately 36% of Nab. magadii protein coding genes could not be assigned a function based on Blast analysis and have been annotated as encoding hypothetical or conserved hypothetical proteins. Furthermore, despite extensive comparative genomic analyses, genes necessary for survival in alkaline conditions could not be identified in Nab. magadii.

CONCLUSIONS

Based on genomic analyses, Nab. magadii is predicted to be metabolically versatile and it could use different carbon and energy sources to sustain growth. Nab. magadii has the genetic potential to adapt to its milieu by intracellular accumulation of inorganic cations and/or neutral organic compounds. The identification of Nab. magadii genes involved in coenzyme biosynthesis is a necessary step toward further reconstruction of the metabolic pathways in halophilic archaea and other extremophiles. The knowledge gained from the genome sequence of this haloalkaliphilic archaeon is highly valuable in advancing the applications of extremophiles and their enzymes.

Bacteriophages φMR299-2 and φNH-4 can eliminate Pseudomonas aeruginosa in the murine lung and on cystic fibrosis lung airway cells

Pseudomonas aeruginosa is a common cause of infection in the lungs of patients with cystic fibrosis (CF). In addition, biofilm formation and antibiotic resistance of Pseudomonas are major problems that can complicate antibiotic therapy. We evaluated the efficacy of using bacteriophages to kill the pathogen in both biofilms and in the murine lung. We isolated and characterized two phages from a local wastewater treatment plant, a myovirus (ϕNH-4) and a podovirus (ϕMR299-2). Both phages were active against clinical isolates of P. aeruginosa. Together, the two phages killed all 9 clinical isolate strains tested, including both mucoid and nonmucoid strains. An equal mixture of the two phages was effective in killing P. aeruginosa NH57388A (mucoid) and P. aeruginosa MR299 (nonmucoid) strains when growing as a biofilm on a cystic fibrosis bronchial epithelial CFBE41o- cell line. Phage titers increased almost 100-fold over a 24-h period, confirming replication of the phage. Furthermore, the phage mix was also effective in killing the pathogen in murine lungs containing 1 × 10⁷ to 2 × 10⁷P. aeruginosa. Pseudomonas was effectively cleared (reduced by a magnitude of at least 3 to 4 log units) from murine lungs in 6 h. Our study demonstrates the efficacy of these two phages in killing clinical Pseudomonas isolates in the murine lung or as a biofilm on a pulmonary cell line and supports the growing interest in using phage therapy for the control and treatment of multidrug-resistant Pseudomonas lung infections in CF patients.

Given the rise in antibiotic resistance, nonantibiotic therapies are required for the treatment of infection. This is particularly true for the treatment of Pseudomonas infection in patients with cystic fibrosis. We have identified two bacterial viruses (bacteriophages) that can kill Pseudomonas growing on human lung cells and in an animal model of lung infection. The use of bacteriophages is particularly appropriate because the killing agent can replicate on the target cell, generating fresh copies of the bacteriophage. Thus, in the presence of a target, the killing agent multiplies. By using two bacteriophages we can reduce the risk of resistant colonies developing at the site of infection. Bacteriophage therapy is an exciting field, and this study represents an important demonstration of efficacy in validated infection models.

Genome sequence of Bifidobacterium breve DPC 6330, a strain isolated from the human intestine

The draft genome of Bifidobacterium breve DPC 6330, isolated from an elderly patient, was determined. B. breve DPC 6330 was previously identified to synthesize the beneficial metabolite conjugated linoleic acid from free linoleic acid. The sequence will allow identification and characterization of the genetic determinants of its putative beneficial properties.

Differences in genotype and virulence among four multidrug-resistant Streptococcus pneumoniae isolates belonging to the PMEN1 clone

We report on the comparative genomics and characterization of the virulence phenotypes of four S. pneumoniae strains that belong to the multidrug resistant clone PMEN1 (Spain^23F ST81). Strains SV35-T23 and SV36-T3 were recovered in 1996 from the nasopharynx of patients at an AIDS hospice in New York. Strain SV36-T3 expressed capsule type 3 which is unusual for this clone and represents the product of an in vivo capsular switch event. A third PMEN1 isolate – PN4595-T23 – was recovered in 1996 from the nasopharynx of a child attending day care in Portugal, and a fourth strain – ATCC700669 – was originally isolated from a patient with pneumococcal disease in Spain in 1984. We compared the genomes among four PMEN1 strains and 47 previously sequenced pneumococcal isolates for gene possession differences and allelic variations within core genes. In contrast to the 47 strains – representing a variety of clonal types – the four PMEN1 strains grouped closely together, demonstrating high genomic conservation within this lineage relative to the rest of the species. In the four PMEN1 strains allelic and gene possession differences were clustered into 18 genomic regions including the capsule, the blp bacteriocins, erythromycin resistance, the MM1-2008 prophage and multiple cell wall anchored proteins. In spite of their genomic similarity, the high resolution chinchilla model was able to detect variations in virulence properties of the PMEN1 strains highlighting how small genic or allelic variation can lead to significant changes in pathogenicity and making this set of strains ideal for the identification of novel virulence determinants.

Horizontal gene transfer in Histophilus somni and its role in the evolution of pathogenic strain 2336, as determined by comparative genomic analyses

BACKGROUND

Pneumonia and myocarditis are the most commonly reported diseases due to Histophilus somni, an opportunistic pathogen of the reproductive and respiratory tracts of cattle. Thus far only a few genes involved in metabolic and virulence functions have been identified and characterized in H. somni using traditional methods. Analyses of the genome sequences of several Pasteurellaceae species have provided insights into their biology and evolution. In view of the economic and ecological importance of H. somni, the genome sequence of pneumonia strain 2336 has been determined and compared to that of commensal strain 129Pt and other members of the Pasteurellaceae.

RESULTS

The chromosome of strain 2336 (2,263,857 bp) contained 1,980 protein coding genes, whereas the chromosome of strain 129Pt (2,007,700 bp) contained only 1,792 protein coding genes. Although the chromosomes of the two strains differ in size, their average GC content, gene density (total number of genes predicted on the chromosome), and percentage of sequence (number of genes) that encodes proteins were similar. The chromosomes of these strains also contained a number of discrete prophage regions and genomic islands. One of the genomic islands in strain 2336 contained genes putatively involved in copper, zinc, and tetracycline resistance. Using the genome sequence data and comparative analyses with other members of the Pasteurellaceae, several H. somni genes that may encode proteins involved in virulence (e.g., filamentous haemaggutinins, adhesins, and polysaccharide biosynthesis/modification enzymes) were identified. The two strains contained a total of 17 ORFs that encode putative glycosyltransferases and some of these ORFs had characteristic simple sequence repeats within them. Most of the genes/loci common to both the strains were located in different regions of the two chromosomes and occurred in opposite orientations, indicating genome rearrangement since their divergence from a common ancestor.

CONCLUSIONS

Since the genome of strain 129Pt was ~256,000 bp smaller than that of strain 2336, these genomes provide yet another paradigm for studying evolutionary gene loss and/or gain in regard to virulence repertoire and pathogenic ability. Analyses of the complete genome sequences revealed that bacteriophage- and transposon-mediated horizontal gene transfer had occurred at several loci in the chromosomes of strains 2336 and 129Pt. It appears that these mobile genetic elements have played a major role in creating genomic diversity and phenotypic variability among the two H. somni strains.

Rates of mutation and host transmission for an Escherichia coli clone over 3 years

Although over 50 complete Escherichia coli/Shigella genome sequences are available, it is only for closely related strains, for example the O55:H7 and O157:H7 clones of E. coli, that we can assign differences to individual evolutionary events along specific lineages. Here we sequence the genomes of 14 isolates of a uropathogenic E. coli clone that persisted for 3 years within a household, including a dog, causing a urinary tract infection (UTI) in the dog after 2 years. The 20 mutations observed fit a single tree that allows us to estimate the mutation rate to be about 1.1 per genome per year, with minimal evidence for adaptive change, including in relation to the UTI episode. The host data also imply at least 6 host transfer events over the 3 years, with 2 lineages present over much of that period. To our knowledge, these are the first direct measurements for a clone in a well-defined host community that includes rates of mutation and host transmission. There is a concentration of non-synonymous mutations associated with 2 transfers to the dog, suggesting some selection pressure from the change of host. However, there are no changes to which we can attribute the UTI event in the dog, which suggests that this occurrence after 2 years of the clone being in the household may have been due to chance, or some unknown change in the host or environment. The ability of a UTI strain to persist for 2 years and also to transfer readily within a household has implications for epidemiology, diagnosis, and clinical intervention.

Gnom(Cmp): a quantitative approach for comparative analysis of closely related genomes of bacterial pathogens

Comparative genome analysis is a powerful approach to understanding the biology of infectious bacterial pathogens. In this study, a quantitative approach, referred to as Gnom(Cmp), was developed to study the microevolution of bacterial pathogens. Although much more time-consuming than existing tools, this procedure provides a much higher resolution. Gnom(Cmp) accomplishes this by establishing genome-wide heterogeneity genotypes, which are then quantified and comparatively analyzed. The heterogeneity genotypes are defined as chromosomal base positions that have multiple variants within particular genomes, resulted from DNA duplications and subsequent mutations. To prove the concept, the procedure was applied on the genomes of 15 Staphylococcus aureus strains, focusing extensively on two pairs of hVISA/VISA strains. hVISA refers to heteroresistant vancomycin-intermediate S. aureus strains and VISA is their VISA mutants. hVISA/VISA displays some remarkable properties. hVISA is susceptible to vancomycin, but VISA mutants emerge soon after a short period of vancomycin therapy, therefore making the pathogen a great model organism for fast-evolving bacterial pathogens. The analysis indicated that Gnom(Cmp) could reveal variants within the genomes, which can be analyzed within the global genome context. Gnom(Cmp) discovered evolutionary hotspots and their dynamics among many closely related, even isogenic genomes. The analysis thus allows the exploration of the molecular mechanisms behind hVISA/VISA evolution, providing a working hypotheses for experimental testing and validation.

Host specific diversity in Lactobacillus johnsonii as evidenced by a major chromosomal inversion and phage resistance mechanisms

Genetic diversity and genomic rearrangements are a driving force in bacterial evolution and niche adaptation. We sequenced and annotated the genome of Lactobacillus johnsonii DPC6026, a strain isolated from the porcine intestinal tract. Although the genome of DPC6026 is similar in size (1.97mbp) and GC content (34.8%) to the sequenced human isolate L. johnsonii NCC 533, a large symmetrical inversion of approximately 750 kb differentiated the two strains. Comparative analysis among 12 other strains of L. johnsonii including 8 porcine, 3 human and 1 poultry isolate indicated that the genome architecture found in DPC6026 is more common within the species than that of NCC 533. Furthermore a number of unique features were annotated in DPC6026, some of which are likely to have been acquired by horizontal gene transfer (HGT) and contribute to protection against phage infection. A putative type III restriction-modification system was identified, as were novel Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) elements. Interestingly, these particular elements are not widely distributed among L. johnsonii strains. Taken together these data suggest intra-species genomic rearrangements and significant genetic diversity within the L. johnsonii species and indicate towards a host-specific divergence of L. johnsonii strains with respect to genome inversion and phage exposure.

Evidence for the horizontal transfer of an unusual capsular polysaccharide biosynthesis locus in marine bacteria

The most intensely studied of the Vibrio vulnificus virulence factors is the capsular polysaccharide (CPS). All virulent strains produce copious amounts of CPS. Acapsular strains are avirulent. The structure of the CPS from the clinical isolate ATCC 27562 is unusual. It is serine modified and contains, surprisingly, N-acetylmuramic acid. We identified the complete 25-kb CPS biosynthesis locus from ATCC 27562. It contained 21 open reading frames and was allelic to O-antigen biosynthesis loci. Two of the genes, murA(CPS) and murB(CPS), were paralogs of the murA(PG) and murB(PG) genes of the peptidoglycan biosynthesis pathway; only a single copy of these genes is present in the strain CMCP6 and YJ016 genomes. Although MurA(CPS) and MurB(CPS) were functional when expressed in Escherichia coli, lesions in either gene had no effect on CPS production, virulence, or growth in V. vulnificus; disruption of 8 other genes within the locus resulted in an acapsular phenotype and attenuated virulence. Thus, murA(CPS) and murB(CPS) were functional but redundant. Comparative genomic analysis revealed that while completely different CPS biosynthesis loci were found in the same chromosomal region in other V. vulnificus strains, most of the CPS locus of ATCC 27562 was conserved in another marine bacterium, Shewanella putrefaciens strain 200. However, the average GC content of the CPS locus was significantly lower than the average GC content of either genome. Furthermore, several of the encoded proteins appeared to be of Gram-positive and archaebacterial origin. These data indicate that the horizontal transfer of intact and partial CPS loci drives CPS diversity in marine bacteria.

Generation of genic diversity among Streptococcus pneumoniae strains via horizontal gene transfer during a chronic polyclonal pediatric infection

Although there is tremendous interest in understanding the evolutionary roles of horizontal gene transfer (HGT) processes that occur during chronic polyclonal infections, to date there have been few studies that directly address this topic. We have characterized multiple HGT events that most likely occurred during polyclonal infection among nasopharyngeal strains of Streptococcus pneumoniae recovered from a child suffering from chronic upper respiratory and middle-ear infections. Whole genome sequencing and comparative genomics were performed on six isolates collected during symptomatic episodes over a period of seven months. From these comparisons we determined that five of the isolates were genetically highly similar and likely represented a dominant lineage. We analyzed all genic and allelic differences among all six isolates and found that all differences tended to occur within contiguous genomic blocks, suggestive of strain evolution by homologous recombination. From these analyses we identified three strains (two of which were recovered on two different occasions) that appear to have been derived sequentially, one from the next, each by multiple recombination events. We also identified a fourth strain that contains many of the genomic segments that differentiate the three highly related strains from one another, and have hypothesized that this fourth strain may have served as a donor multiple times in the evolution of the dominant strain line. The variations among the parent, daughter, and grand-daughter recombinant strains collectively cover greater than seven percent of the genome and are grouped into 23 chromosomal clusters. While capturing in vivo HGT, these data support the distributed genome hypothesis and suggest that a single competence event in pneumococci can result in the replacement of DNA at multiple non-adjacent loci.

Bacterial infections have long been studied using Koch's postulates wherein the paradigm is that a single clone leads to a given infection. Over the past decade, it has become clear that chronic bacterial infections often do not fit this paradigm. Instead these are associated with the presence of multiple strains or species (polyclonal) of bacteria that are organized into highly structured communities, termed biofilms, which can persist in the body and are recalcitrant to antibiotic treatment. In addition, there is extensive evidence that bacteria can incorporate genes from neighboring bacteria into their own genomes. This process can produce new strains and is known as horizontal gene transfer. In this study, we investigated for the first time, the tempo and relevance of gene transfer among bacterial strains of Streptococcus pneumoniae during a naturally occurring chronic childhood infection. We identified extensive gene transfer among multiple infecting strains, by sequencing of isolates recovered sequentially over a seven-month period. This gene transfer may serve as a counterpoint to the host's adaptive immune response and help explain the phenomenon of bacterial persistence, since, as occurs with some chronic viral and parasitic infections, the immune system may become overwhelmed by a set of related strains.

A genomic survey of positive selection in Burkholderia pseudomallei provides insights into the evolution of accidental virulence

Certain environmental microorganisms can cause severe human infections, even in the absence of an obvious requirement for transition through an animal host for replication ("accidental virulence"). To understand this process, we compared eleven isolate genomes of Burkholderia pseudomallei (Bp), a tropical soil microbe and causative agent of the human and animal disease melioidosis. We found evidence for the existence of several new genes in the Bp reference genome, identifying 282 novel genes supported by at least two independent lines of supporting evidence (mRNA transcripts, database homologs, and presence of ribosomal binding sites) and 81 novel genes supported by all three lines. Within the Bp core genome, 211 genes exhibited significant levels of positive selection (4.5%), distributed across many cellular pathways including carbohydrate and secondary metabolism. Functional experiments revealed that certain positively selected genes might enhance mammalian virulence by interacting with host cellular pathways or utilizing host nutrients. Evolutionary modifications improving Bp environmental fitness may thus have indirectly facilitated the ability of Bp to colonize and survive in mammalian hosts. These findings improve our understanding of the pathogenesis of melioidosis, and establish Bp as a model system for studying the genetics of accidental virulence.

Derivation of Escherichia coli O157:H7 from its O55:H7 precursor

There are 29 E. coli genome sequences available, mostly related to studies of species diversity or mode of pathogenicity, including two genomes of the well-known O157:H7 clone. However, there have been no genome studies of closely related clones aimed at exposing the details of evolutionary change. Here we sequenced the genome of an O55:H7 strain, closely related to the major pathogenic O157:H7 clone, with published genome sequences, and undertook comparative genomic and proteomic analysis. We were able to allocate most differences between the genomes to individual mutations, recombination events, or lateral gene transfer events, in specific lineages. Major differences include a type II secretion system present only in the O55:H7 chromosome, fewer type III secretion system effectors in O55:H7, and 19 phage genomes or phagelike elements in O55:H7 compared to 23 in O157:H7, with only three common to both. Many other changes were found in both O55:H7 and O157:H7 lineages, but in general there has been more change in the O157:H7 lineages. For example, we found 50% more synonymous mutational substitutions in O157:H7 compared to O55:H7. The two strains also diverged at the proteomic level. Mutational synonymous SNPs were used to estimate a divergence time of 400 years using a new clock rate, in contrast to 14,000 to 70,000 years using the traditional clock rates. The same approaches were applied to three closely related extraintestinal pathogenic E. coli genomes, and similar levels of mutation and recombination were found. This study revealed for the first time the full range of events involved in the evolution of the O157:H7 clone from its O55:H7 ancestor, and suggested that O157:H7 arose quite recently. Our findings also suggest that E. coli has a much lower frequency of recombination relative to mutation than was observed in a comparable study of a Vibrio cholerae lineage.

Genome sequence of Aggregatibacter actinomycetemcomitans serotype c strain D11S-1

Aggregatibacter actinomycetemcomitans is a major etiological agent of periodontitis. Here we report the complete genome sequence of serotype c strain D11S-1, which was recovered from the subgingival plaque of a patient diagnosed with generalized aggressive periodontitis.

Tracing the evolution of competence in Haemophilus influenzae

Natural competence is the genetically encoded ability of some bacteria to take up DNA from the environment. Although most of the incoming DNA is degraded, occasionally intact homologous fragments can recombine with the chromosome, displacing one resident strand. This potential to use DNA as a source of both nutrients and genetic novelty has important implications for the ecology and evolution of competent bacteria. However, it is not known how frequently competence changes during evolution, or whether non-competent strains can persist for long periods of time. We have previously studied competence in H. influenzae and found that both the amount of DNA taken up and the amount recombined varies extensively between different strains. In addition, several strains are unable to become competent, suggesting that competence has been lost at least once. To investigate how many times competence has increased or decreased during the divergence of these strains, we inferred the evolutionary relationships of strains using the largest datasets currently available. However, despite the use of three datasets and multiple inference methods, few nodes were resolved with high support, perhaps due to extensive mixing by recombination. Tracing the evolution of competence in those clades that were well supported identified changes in DNA uptake and/or transformation in most strains. The recency of these events suggests that competence has changed frequently during evolution but the poor support of basal relationships precludes the determination of whether non-competent strains can persist for long periods of time. In some strains, changes in transformation have occurred that cannot be due to changes in DNA uptake, suggesting that selection can act on transformation independent of DNA uptake.

Genetics and virulence association of the Shigella flexneri sit iron transport system

The sit-encoded iron transport system is present within pathogenicity islands in all Shigella spp. and some pathogenic Escherichia coli strains. The islands contain numerous insertion elements and sequences with homology to bacteriophage genes. The Shigella flexneri sit genes can be lost as a result of deletion within the island. The formation of deletions was dependent upon RecA and occurred at relatively high frequency. This suggests that the sit region is inherently unstable, yet sit genes are maintained in all of the clinical isolates tested. Characterization of the sitABCD genes in S. flexneri indicates that they encode a ferrous iron transport system, although the genes are induced aerobically. The sit genes provide a competitive advantage to S. flexneri growing within epithelial cells, and a sitA mutant is outcompeted by the wild type in cultured epithelial cells. The Sit system is also required for virulence in a mouse lung model. The sitA mutant was able to infect the mice and induce a protective immune response but was avirulent compared to its wild-type parent strain.