Platanus_B: an accurate de novo assembler for bacterial genomes using an iterative error-removal process

Mobile genetic element-driven genomic changes in a community-associated methicillin-resistant Staphylococcus aureus clone during its transmission in a regional community outbreak in Japan

Community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) infections are now a public health concern in both community and healthcare settings worldwide. We previously identified a suspected case of a maternity clinic-centred outbreak of CA-MRSA skin infection in a regional community in Japan by PFGE-based analysis. In this study, we performed genome sequence-based analyses of 151 CA-MRSA isolates, which included not only outbreak-related isolates that we previously defined based on identical or similar PFGE patterns but also other isolates obtained during the same period in the same region. Our analysis accurately defined 133 isolates as outbreak-related isolates, collectively called the TDC clone. They belonged to a CA-MRSA lineage in clonal complex (CC) 30, known as the South West Pacific (SWP) clone. A high-resolution phylogenetic analysis of these isolates combined with their epidemiological data revealed that the TDC clone was already present and circulating in the region before the outbreak was recognized, and only the isolates belonging to two sublineages (named SL4 and SL5) were directly involved in the outbreak. Long persistence in patients/carriers and frequent intrahousehold transmission of the TDC clone were also revealed by this analysis. Moreover, by systematic analyses of the genome changes that occurred in this CA-MRSA clone during transmission in the community, we revealed that most variations were associated with mobile genetic elements (MGEs). Variant PFGE types were generated by alterations of prophages and genomic islands or insertion sequence (IS)-mediated insertion of a plasmid or a sequence of unknown origin. Dynamic changes in plasmid content, which were linked to changes in antimicrobial resistance profiles in specific isolates, were generated by frequent gain and loss of plasmids, most of which were self-transmissible or mobilizable. The introduction of IS256 by a plasmid (named pTDC02) into sublineage SL5 led to SL5-specific amplification of IS256, and amplified IS256 copies were involved in some of the structural changes of chromosomes and plasmids and generated variations in the repertoire of virulence-related genes in limited isolates. These data revealed how CA-MRSA genomes change during transmission in the community and how MGEs are involved in this process.

Complete genome sequences of 11 Streptococcus pneumoniae strains isolated from acute infectious purpura fulminans, sepsis, and pneumonia

The complete genome sequences of 11 Japanese Streptococcus pneumoniae isolates were determined by hybrid assembly of long and short reads, including two strains isolated from patients with acute infectious purpura fulminans, six strains from patients with sepsis, and three strains from patients with pneumonia.

Description of three new bat-associated species of hard ticks (Acari, Ixodidae) from Japan

In Eurasia, the geographically most widespread ixodid tick species of the bat families Rhinolophidae Gray, Vespertilionidae Gray, and Miniopteridae Dobson were considered to belong to four species, Ixodesvespertilionis Koch, I.collaris Hornok, I.ariadnae Hornok, and I.simplex Neumann. Previous data attest that bat-associated tick species from Eastern Asia show remarkable genetic difference from the above four tick species, but in the absence of detailed morphological comparison these were regarded as conspecific. In this study we compensate for this lack of data on three bat-associated tick species, reporting their morphological comparison, as well as molecular and phylogenetic relationships. According to the results we describe the females of three tick species new to science, i.e., I.nipponrhinolophi Hornok & Takano, sp. nov., I.fuliginosus Hornok & Takano, sp. nov., and I.fujitai Hornok & Takano, sp. nov. In case of all three new tick species the cytochrome c oxidase subunit (coxI) gene showed remarkably high sequence differences from the species that they previously were thought to belong to, well exceeding the average limit delineating ixodid tick species. This, as well as observed morphological differences fully justify their taxonomical status as new species.

Core single nucleotide polymorphism analysis reveals transmission of Mycobacterium marinum between animal and environmental sources in two aquaria

Mycobacterium marinum is a slow-growing, photochromogenic nontuberculous mycobacterium, which can cause mycobacteriosis in various animals, including humans. Several cases of fish mycobacteriosis have been reported to date. Mycobacterium marinum has also been isolated from aquatic environmental sources such as water, sand, biofilms, and plants in the natural environments. Hence, we hypothesized that a wide variety of sources could be involved in the transmission of M. marinum. In this study, we tested this hypothesis by isolating M. marinum from various sources such as fish, invertebrates, seagrass, periphytons, biofilms, sand, and/or water in two aquaria in Japan and conducting a phylogenetic analysis based on single-nucleotide polymorphisms (SNPs) using whole-genome sequences of the isolated strains. The analysis revealed that the strains from animal and environmental sources belonged to the same clusters. This molecular-based study epidemiologically confirmed that various sources, including fish, invertebrates, and environmental sources, could be involved in transmission of M. marinum in a closed-rearing environment. This is the first report where M. marinum was isolated from different sources, and various transmission routes were confirmed in actual cases, which provided essential information to improve the epidemiology of M. marinum.

Escherichia cryptic clade I is an emerging source of human intestinal pathogens

BACKGROUND

Within the genus Escherichia, several monophyletic clades other than the traditionally defined species have been identified. Of these, cryptic clade I (C-I) appears to represent a subspecies of E. coli, but due to the difficulty in distinguishing it from E. coli sensu stricto, the population structure and virulence potential of C-I are unclear.

RESULTS

We defined a set of true C-I strains (n = 465), including a Shiga toxin 2a (Stx2a)-producing isolate from a patient with bloody diarrhoea identified by the retrospective analyses using a C-I-specific detection system. Through genomic analysis of 804 isolates from the cryptic clades, including these C-I strains, we revealed their global population structures and the marked accumulation of virulence genes and antimicrobial resistance genes in C-I. In particular, half of the C-I strains contained hallmark virulence genes of Stx-producing E. coli (STEC) and/or enterotoxigenic E. coli (ETEC). We also found the host-specific distributions of virulence genes, which suggests bovines as the potential source of human infections caused by STEC- and STEC/ETEC hybrid-type C-I strains, as is known in STEC.

CONCLUSIONS

Our findings demonstrate the emergence of human intestinal pathogens in C-I lineage. To better understand the features of C-I strains and their infections, extensive surveillance and larger population studies of C-I strains are needed. The C-I-specific detection system developed in this study will be a powerful tool for screening and identifying C-I strains.

Genomic and phylogenetic characterization of Elizabethkingia anophelis strains: The first two cases of life-threatening infection in Japan

OBJECTIVE

Elizabethkingia anophelis causes meningitis, bloodstream infections, and respiratory infections in immunocompromised individuals. We examined two E. anophelis strains isolated from the first life-threatening cases caused by this species in Japan to determine the phylogenetic origin and genomic features of them.

METHODS

We performed whole genome-based analysis to clarify the genetic relationship for the two strains (EK0004 and EK0079) and Elizabethkingia sp. strains isolated from worldwide and to characterize the genomic features such as the prevalence of virulence- and antimicrobial resistance (AMR)-related genes.

PATIENTS

A 29-year-old man with hepatosplenic T-cell lymphoma and a 52-year-old man with systemic lupus erythematosus developed fatal bacteremia and meningitis due to E. anophelis, respectively.

RESULTS

Two strains, EK0004 and EK0079, were genetically different but most closely related to the strains isolated from the largest outbreak in Wisconsin, USA from 2015 to 2016, and the strain isolated from cerebrospinal fluid of a patient in Florida, USA in 1982, respectively. The two strains contained AMR-related genes such as those encoding for an extended-spectrum β-lactamase and multiple metallo-β-lactamases and several virulence-related genes such as capsular polysaccharide synthesis gene clusters.

CONCLUSIONS

Although further functional analyses are required to understand the virulence of these clones, these finding suggests that enough caution of E. anophelis infection in immunocompromised patients is required since the number of infections by this species is increasing outside Japan.

Discovery of diphenyl ether-degrading Streptomyces strains by direct screening based on ether bond-cleaving activity

Diphenyl ethers (DEs), which are widely used in the agricultural and chemical industries, have become hazardous contaminants in the environment. Although several DE-degrading bacteria have been reported, discovering new types of such microorganisms could enhance understanding of the degradation mechanism in the environment. In this study, we used a direct screening method based on detection of ether bond-cleaving activity to screen for microorganisms that degrade 4,4'-dihydroxydiphenyl ether (DHDE) as a model DE. Microorganisms isolated from soil samples were incubated with DHDE, and strains producing hydroquinone via ether bond cleavage were selected using hydroquinone-sensitive Rhodanine reagent. This screening procedure resulted in the isolation of 3 bacteria and 2 fungi that transform DHDE. Interestingly, all of the isolated bacteria belonged to one genus, Streptomyces. To our knowledge, these are the first microorganisms of the genus Streptomyces shown to degrade a DE. Streptomyces sp. TUS-ST3 exhibited high and stable DHDE-degrading activity. HPLC, LC-MS, and GC-MS analyses revealed that strain TUS-ST3 converts DHDE to its hydroxylated analogue and generates hydroquinone as an ether bond-cleavage product. Strain TUS-ST3 also transformed DEs other than DHDE. In addition, glucose-grown TUS-ST3 cells began to transform DHDE after incubation with this compound for 12 h, and produced 75 μM hydroquinone in 72 h. These activities of streptomycetes may play an important role in DE degradation in the environment. We also report the whole genome sequence of strain TUS-ST3.

Draft Genome Sequence of a Mycobacterium chelonae subsp. bovis Strain Isolated from a Baikal Seal (Pusa sibirica) in Captivity

Mycobacterium chelonae is a nontuberculous mycobacterium that causes infections in various animals, including humans. In this study, we report the draft genome sequence of M. chelonae subsp. bovis strain NJB1701, which was isolated from a Baikal seal (Pusa sibirica) in captivity in Japan.

Draft Genome Sequences of Endophytic Pseudomonas Strains, Isolated from the Interior of Brassicaceae Plants

Members of the genus Pseudomonas have often been studied as agricultural biocontrol agents to activate plant immune responses. Here, we report the draft genome sequences of six Pseudomonas strains that were isolated from the interior of Brassicaceae plants.

Alteration of a Shiga toxin-encoding phage associated with a change in toxin production level and disease severity in Escherichia coli

Among the nine clades of Shiga toxin (Stx)-producing Escherichia coli O157:H7, clade 8 is thought to be highly pathogenic, as it causes severe disease more often than other clades. Two subclades have been proposed, but there are conflicting reports on intersubclade differences in Stx2 levels, although Stx2 production is a risk factor for severe disease development. The global population structure of clade 8 has also yet to be fully elucidated. Here, we present genome analyses of a global clade 8 strain set (n=510), including 147 Japanese strains sequenced in this study. The complete genome sequences of 18 of the 147 strains were determined to perform detailed clade-wide genome analyses together with 17 publicly available closed genomes. Intraclade variations in Stx2 production level and disease severity were also re-evaluated within the phylogenetic context. Based on phylogenomic analysis, clade 8 was divided into four lineages corresponding to the previously proposed SNP genotypes (SGs): SG8_30, SG8_31A, SG8_31B and SG8_32. SG8_30 and the common ancestor of the other SGs were first separated, with SG8_31A and SG8_31B emerging from the latter and SG8_32 emerging from SG8_31B. Comparison of 35 closed genomes revealed the overall structure of chromosomes and pO157 virulence plasmids and the prophage contents to be well conserved. However, Stx2a phages exhibit notable genomic diversity, even though all are integrated into the argW locus, indicating that subtype changes in Stx2a phage occurred from the γ subtype to its variant (γ_v1) in SG8_31A and from γ to δ in SG8_31B and SG8_32 via replacement of parts or almost entire phage genomes, respectively. We further show that SG8_30 strains (all carrying γ Stx2a phages) produce significantly higher levels of Stx2 and cause severe disease more frequently than SG8_32 strains (all carrying δ Stx2a phages). Clear conclusions on SG8_31A and SG8_31B cannot be made due to the small number of strains available, but as SG8_31A (carrying γ_v1 Stx2a phages) contains strains that produce much more Stx2 than SG8_30 strains, attention should also be paid to this SG.

Draft Genome Sequences of Eight Mycobacterium montefiorense Strains Isolated from Salamanders in Captivity

Mycobacterium montefiorense is a nontuberculous mycobacterium that causes infections in fish and salamanders. Here, we report annotated draft genome sequences of eight strains that were isolated in 2014 and 2018 from salamanders reared in an aquarium in Japan.

Draft Genome Sequences of 25 Mycobacterium marinum Strains Isolated from Animals and Environmental Components in Aquaria and an Aquaculture Farm

Mycobacterium marinum is a ubiquitous nontuberculous mycobacterium that causes infections in various animals. Here, we report the annotated draft genome sequences of 25 strains isolated from vertebrates, invertebrates, and environmental components in aquaria and an aquaculture farm in Japan, sampled between 2015 and 2020.

Helicobacter cinaedi is a human-adapted lineage in the Helicobacter cinaedi/canicola/'magdeburgensis' complex

Helicobacter cinaedi is an enterohepatic Helicobacter that causes bacteremia and other diseases in humans. While H. cinaedi-like strains are isolated from animals, including dog isolates belonging to a recently proposed H. canicola, little is known about the genetic differences between H. cinaedi and these animal isolates. Here, we sequenced 43 H. cinaedi- or H. canicola-like strains isolated from humans, hamsters, rats and dogs and collected 81 genome sequences of H. cinaedi, H. canicola and other enterohepatic Helicobacter strains from public databases. Genomic comparison of these strains identified four distinct clades (clades I–IV) in H. cinaedi/canicola/‘magderbugensis’ (HCCM) complex. Among these, clade I corresponds to H. cinaedi sensu stricto and represents a human-adapted lineage in the complex. We identified several genomic features unique to clade I. They include the accumulation of antimicrobial resistance-related mutations that reflects the human association of clade I and the larger genome size and the presence of a CRISPR-Cas system and multiple toxin-antitoxin and restriction-modification systems, both of which indicate the contribution of horizontal gene transfer to the evolution of clade I. In addition, nearly all clade I strains but only a few strains belonging to one minor clade contained a highly variable genomic region encoding a type VI secretion system (T6SS), which could play important roles in gut colonization by killing competitors or inhibiting their growth. We also developed a method to systematically search for H. cinaedi sequences in large metagenome data sets based on the results of genome comparison. Using this method, we successfully identified multiple HCCM complex-containing human faecal metagenome samples and obtained the sequence information covering almost the entire genome of each strain. Importantly, all were clade I strains, supporting our conclusion that H. cinaedi sensu stricto is a human-adapted lineage in the HCCM complex.

Genome-based, phenotypic and chemotaxonomic classification of Faecalibacterium strains: proposal of three novel species Faecalibacterium duncaniae sp. nov., Faecalibacterium hattorii sp. nov. and Faecalibacterium gallinarum sp. nov

Faecalibacterium prausnitzii is one of the most important butyrate-producing bacteria in the human gut. Previous studies have suggested the presence of several phylogenetic groups, with differences at the species level, in the species, and a taxonomic re-evaluation is thus essential for further understanding of ecology of the important human symbiont. Here we examine the phenotypic, physiological, chemotaxonomic and phylogenomic characteristics of six F. prausnitzii strains (BCRC 81047T=ATCC 27768T, A2-165T=JCM 31915T, APC918/95b=JCM 39207, APC942/30−2=JCM 39208, APC924/119=JCM 39209 and APC922/41−1T=JCM 39210T) deposited in public culture collections with two reference strains of Faecalibacterium butyricigenerans JCM 39212T and Faecalibacterium longum JCM 39211T. Faecalibacterium sp. JCM 17207T isolated from caecum of broiler chicken was also included. Three strains of F. prausnitzii (BCRC 81047T, JCM 39207 and JCM 39209) shared more than 96.6 % average nucleotide identity (ANI) and 69.6 % digital DNA–DNA hybridization (dDDH) values, indicating that the three strains are members of the same species. On the other hand, the remaining three strains of F. prausnitzii (JCM 31915T, JCM 39208 and JCM 39210T) were clearly separated from the above three strains based on the ANI and dDDH values. Rather, JCM 39208 showed ANI and dDDH values over the cut-off values of species discrimination (>70 % dDDH and >95–96 % ANI) with F. longum JCM 39211T, whereas JCM 31915T, JCM 39210T and JCM 17207T did not share dDDH and ANI values over the currently accepted cut-off values with any of the tested strains, including among them. Furthermore, the cellular fatty acid patterns of these strains were slightly different from other F. prausnitzii strains. Based on the collected data, F. prausnitzii JCM 31915T, F. prausnitzii JCM 39210T and Faecalibacterium sp. JCM 17207T represent three novel species of the genus Faecalibacterium , for which the names Faecalibacterium duncaniae sp. nov. (type strain JCM 31915T=DSM 17677T=A2-165T), Faecalibacterium hattorii sp. nov. (type strain JCM 39210T=DSM 107841T=APC922/41-1T) and Faecalibacterium gallinarum sp. nov. (type strain JCM 17207T=DSM 23680T=ic1379T) are proposed.

Global population structure of the Serratia marcescens complex and identification of hospital-adapted lineages in the complex

Serratia marcescens is an important nosocomial pathogen causing various opportunistic infections, such as urinary tract infections, bacteremia and sometimes even hospital outbreaks. The recent emergence and spread of multidrug-resistant (MDR) strains further pose serious threats to global public health. This bacterium is also ubiquitously found in natural environments, but the genomic differences between clinical and environmental isolates are not clear, including those between S. marcescens and its close relatives. In this study, we performed a large-scale genome analysis of S. marcescens and closely related species (referred to as the 'S. marcescens complex'), including more than 200 clinical and environmental strains newly sequenced here. Our analysis revealed their phylogenetic relationships and complex global population structure, comprising 14 clades, which were defined based on whole-genome average nucleotide identity. Clades 10, 11, 12 and 13 corresponded to S. nematodiphila, S. marcescens sensu stricto, S. ureilytica and S. surfactantfaciens, respectively. Several clades exhibited distinct genome sizes and GC contents and a negative correlation of these genomic parameters was observed in each clade, which was associated with the acquisition of mobile genetic elements (MGEs), but different types of MGEs, plasmids or prophages (and other integrative elements), were found to contribute to the generation of these genomic variations. Importantly, clades 1 and 2 mostly comprised clinical or hospital environment isolates and accumulated a wide range of antimicrobial resistance genes, including various extended-spectrum β-lactamase and carbapenemase genes, and fluoroquinolone target site mutations, leading to a high proportion of MDR strains. This finding suggests that clades 1 and 2 represent hospital-adapted lineages in the S. marcescens complex although their potential virulence is currently unknown. These data provide an important genomic basis for reconsidering the classification of this group of bacteria and reveal novel insights into their evolution, biology and differential importance in clinical settings.

MetaPlatanus: a metagenome assembler that combines long-range sequence links and species-specific features

De novo metagenome assembly is effective in assembling multiple draft genomes, including those of uncultured organisms. However, heterogeneity in the metagenome hinders assembly and introduces interspecies misassembly deleterious for downstream analysis. For this purpose, we developed a hybrid metagenome assembler, MetaPlatanus. First, as a characteristic function, it assembles the basic contigs from accurate short reads and then iteratively utilizes long-range sequence links, species-specific sequence compositions, and coverage depth. The binning information was also used to improve contiguity. Benchmarking using mock datasets consisting of known bacteria with long reads or mate pairs revealed the high contiguity MetaPlatanus with a few interspecies misassemblies. For published human gut data with nanopore reads from potable sequencers, MetaPlatanus assembled many biologically important elements, such as coding genes, gene clusters, viral sequences, and over-half bacterial genomes. In the benchmark with published human saliva data with high-throughput nanopore reads, the superiority of MetaPlatanus was considerably more evident. We found that some high-abundance bacterial genomes were assembled only by MetaPlatanus as near-complete. Furthermore, MetaPlatanus can circumvent the limitations of highly fragmented assemblies and frequent interspecies misassembles obtained by the other tools. Overall, the study demonstrates that MetaPlatanus could be an effective approach for exploring large-scale structures in metagenomes.

Vibrio Clade 3.0: New Vibrionaceae Evolutionary Units Using Genome-Based Approach

Currently, over 190 species in family Vibrionaceae, including not-yet-cultured taxa, have been described and classified into over nine genera, in which the number of species has doubled compared to the previous vibrio evolutionary update (Vibrio Clade 2.0) (Sawabe et al. 2014). In this study, "Vibrio Clade 3.0," the second update of the molecular phylogenetic analysis was performed based on nucleotide sequences of eight housekeeping genes (8-HKGs) retrieved from genome sequences, including 22 newly determined genomes. A total of 51 distinct clades were observed, of which 21 clades are newly described. We further evaluated the delineation powers of the clade classification based on nucleotide sequences of 34 single-copy genes and 11 ribosomal protein genes (11-RPGs) retrieved from core-genome sequences; however, the delineation power of 8-HKGs is still high and that gene set can be reliably used for the classification and identification of Vibrionaceae. Furthermore, the 11-RPGs set proved to be useful in identifying uncultured species among metagenome-assembled genome (MAG) and/or single-cell genome-assembled genome (SAG) pools. This study expands the awareness of the diversity and evolutionary history of the family Vibrionaceae and accelerates the taxonomic applications in classifying as not-yet-cultured taxa among MAGs and SAGs.

Oligosaccharide Metabolism and Lipoteichoic Acid Production in Lactobacillus gasseri and Lactobacillus paragasseri

Lactobacillus gasseri and Lactobacillus paragasseri are human commensal lactobacilli that are candidates for probiotic application. Knowledge of their oligosaccharide metabolic properties is valuable for synbiotic application. The present study characterized oligosaccharide metabolic systems and their impact on lipoteichoic acid (LTA) production in the two organisms, i.e., L. gasseri JCM 1131^T and L. paragasseri JCM 11657. The two strains grew well in medium with glucose but poorly in medium with raffinose, and growth rates in medium with kestose differed between the strains. Oligosaccharide metabolism markedly influenced their LTA production, and apparent molecular size of LTA in electrophoresis recovered from cells cultured with glucose and kestose differed from that from cells cultured with raffinose in the strains. On the other hand, more than 15-fold more LTA was observed in the L. gasseri cells cultured with raffinose when compared with glucose or kestose after incubation for 15 h. Transcriptome analysis identified glycoside hydrolase family 32 enzyme as a potential kestose hydrolysis enzyme in the two strains. Transcriptomic levels of multiple genes in the dlt operon, involved in D-alanine substitution of LTA, were lower in cells cultured with raffinose than in those cultured with kestose or glucose. This suggested that the different sizes of LTA observed among the carbohydrates tested were partly due to different levels of alanylation of LTA. The present study indicates that available oligosaccharide has the impact on the LTA production of the industrially important lactobacilli, which might influence their probiotic properties.

Draft Genome Sequence of Neisseria mucosa Strain HSUH001, Isolated from an Aggressive Periodontal Lesion

We report the draft genome sequence (143 contigs, with a total length of 2,424,805 bp and an N ₅₀ value of 36,066 bp) of a bacterium isolated from an aggressive periodontal lesion in a patient. We assigned strain HSUH001 to Neisseria mucosa through a multilocus sequence analysis.

Lactobacillus nasalidis sp. nov., isolated from the forestomach of a captive proboscis monkey (Nasalis larvatus)

Three strains (YZ01^T, YZ02 and YZ03) of Gram-stain-positive, facultatively anaerobic rods were isolated from the forestomach contents collected from a captive male proboscis monkey (Nasalis larvatus) at Yokohama Zoo in Japan. Phylogenetic analysis of the 16S rRNA gene sequences revealed that these strains belonged to the genus Lactobacillus. Based on the sequence similarity of the 16S rRNA gene, Lactobacillus delbrueckii subsp. indicus JCM 15610^T was the closest phylogenetic neighbour to YZ01^T. Sequence analyses of two partial concatenated housekeeping genes, the RNA polymerase alpha subunit (rpoA) and phenylalanyl-tRNA synthase alpha subunit (pheS) also indicated that the novel strains belonged to the genus Lactobacillus. The average nucleotide identity and digital DNA-DNA hybridization (dDDH) between L. delbrueckii subsp. indicus and YZ01^T were 85.9 and 31.4 %, respectively. The phylogenetic tree based on the whole genomic data of strains YZ01^T, YZ02 and YZ03 suggested that these three strains formed a single monophyletic cluster in the genus Lactobacillus, indicating that it belonged to a new species. The DNA G+C content of strain YZ01^T was 51.6 mol%. The major fatty acids were C_16 : 0 and C_18 : 1 ω9c. Therefore, based on phylogenetic, phenotypic and physiological evidence, strains YZ01^T, YZ02 and YZ03 represent a novel species of the genus Lactobacillus, for which the name Lactobacillus nasalidis sp. nov. is proposed with the type strain YZ01^T (=JCM 33769^T=DSM 110539^T).

Prophages integrating into prophages: A mechanism to accumulate type III secretion effector genes and duplicate Shiga toxin-encoding prophages in Escherichia coli

Bacteriophages (or phages) play major roles in the evolution of bacterial pathogens via horizontal gene transfer. Multiple phages are often integrated in a host chromosome as prophages, not only carrying various novel virulence-related genetic determinants into host bacteria but also providing various possibilities for prophage-prophage interactions in bacterial cells. In particular, Escherichia coli strains such as Shiga toxin (Stx)-producing E. coli (STEC) and enteropathogenic E. coli (EPEC) strains have acquired more than 10 prophages (up to 21 prophages), many of which encode type III secretion system (T3SS) effector gene clusters. In these strains, some prophages are present at a single locus in tandem, which is usually interpreted as the integration of phages that use the same attachment (att) sequence. Here, we present phages integrating into T3SS effector gene cluster-associated loci in prophages, which are widely distributed in STEC and EPEC. Some of the phages integrated into prophages are Stx-encoding phages (Stx phages) and have induced the duplication of Stx phages in a single cell. The identified attB sequences in prophage genomes are apparently derived from host chromosomes. In addition, two or three different attB sequences are present in some prophages, which results in the generation of prophage clusters in various complex configurations. These phages integrating into prophages represent a medically and biologically important type of inter-phage interaction that promotes the accumulation of T3SS effector genes in STEC and EPEC, the duplication of Stx phages in STEC, and the conversion of EPEC to STEC and that may be distributed in other types of E. coli strains as well as other prophage-rich bacterial species.

Multiple prophages are often integrated in a bacterial host chromosome and some are present at a single locus in tandem. The most striking examples are Shiga toxin (Stx)-producing and enteropathogenic Escherichia coli (STEC and EPEC) strains, which usually contain more than 10 prophages (up to 21). Many of them encode a cluster of type III secretion system (T3SS) effector genes, contributing the acquisition of a large number of effectors (>30) by STEC and EPEC. Here, we describe prophages integrating into T3SS effector gene cluster-associated loci in prophages, which are widely distributed in STEC and EPEC. Two or three different attachment sequences derived from host chromosomes are present in some prophages, generating prophage clusters in various complex configurations. Of note, some of such phages integrating into prophages are Stx-encoding phages (Stx phages) and have induced the duplication of Stx phages. Thus, “prophage-in-prophage” represents an important inter-phage interaction as they can promote not only the accumulation of T3SS effectors in STEC and EPEC but also the duplication of Stx phages and the conversion of EPEC to STEC.

Niche-specific adaptation of Lactobacillus helveticus strains isolated from malt whisky and dairy fermentations

Lactobacillus helveticus is a well characterized lactobacillus for dairy fermentations that is also found in malt whisky fermentations. The two environments contain considerable differences related to microbial growth, including the presence of different growth inhibitors and nutrients. The present study characterized L. helveticus strains originating from dairy fermentations (called milk strains hereafter) and malt whisky fermentations (called whisky strains hereafter) by in vitro phenotypic tests and comparative genomics. The whisky strains can tolerate ethanol more than the milk strains, whereas the milk strains can tolerate lysozyme and lactoferrin more than the whisky strains. Several plant-origin carbohydrates, including cellobiose, maltose, sucrose, fructooligosaccharide and salicin, were generally metabolized only by the whisky strains, whereas milk-derived carbohydrates, i.e. lactose and galactose, were metabolized only by the milk strains. Milk fermentation properties also distinguished the two groups. The general genomic characteristics, including genomic size, number of coding sequences and average nucleotide identity values, differentiated the two groups. The observed differences in carbohydrate metabolic properties between the two groups correlated with the presence of intact specific enzymes in glycoside hydrolase (GH) families GH1, GH4, GH13, GH32 and GH65. Several GHs in the milk strains were inactive due to the presence of stop codon(s) in genes encoding the GHs, and the inactivation patterns of the genes encoding specific enzymes assigned to GH1 in the milk strains suggested a possible diversification manner of L. helveticus strains. The present study has demonstrated how L. helveticus strains have adapted to their habitats.