Unified Workflow for the Rapid and In-Depth Characterization of Bacterial Proteomes.

Bacteria are the most abundant and diverse organisms among the kingdoms of life. Due to this excessive variance, finding a unified, comprehensive, and safe workflow for quantitative bacterial proteomics is challenging. In this study, we have systematically evaluated and optimized sample preparation, mass spectrometric data acquisition, and data analysis strategies in bacterial proteomics. We investigated workflow performances on six representative species with highly different physiologic properties to mimic bacterial diversity. The best sample preparation strategy was a cell lysis protocol in 100% trifluoroacetic acid followed by an in-solution digest. Peptides were separated on a 30-min linear microflow liquid chromatography gradient and analyzed in data-independent acquisition mode. Data analysis was performed with DIA-NN using a predicted spectral library. Performance was evaluated according to the number of identified proteins, quantitative precision, throughput, costs, and biological safety. With this rapid workflow, over 40% of all encoded genes were detected per bacterial species. We demonstrated the general applicability of our workflow on a set of 23 taxonomically and physiologically diverse bacterial species. We could confidently identify over 45,000 proteins in the combined dataset, of which 30,000 have not been experimentally validated before. Our work thereby provides a valuable resource for the microbial scientific community. Finally, we grew Escherichia coli and Bacillus cereus in replicates under 12 different cultivation conditions to demonstrate the high-throughput suitability of the workflow. The proteomic workflow we present in this manuscript does not require any specialized equipment or commercial software and can be easily applied by other laboratories to support and accelerate the proteomic exploration of the bacterial kingdom.

Non-canonical transcriptional start sites in E. coli O157:H7 EDL933 are regulated and appear in surprisingly high numbers.

Analysis of genome wide transcription start sites (TSSs) revealed an unexpected complexity since not only canonical TSS of annotated genes are recognized by RNA polymerase. Non-canonical TSS were detected antisense to, or within, annotated genes as well new intergenic (orphan) TSS, not associated with known genes. Previously, it was hypothesized that many such signals represent noise or pervasive transcription, not associated with a biological function. Here, a modified Cappable-seq protocol allows determining the primary transcriptome of the enterohemorrhagic E. coli O157:H7 EDL933 (EHEC). We used four different growth media, both in exponential and stationary growth phase, replicated each thrice. This yielded 19,975 EHEC canonical and non-canonical TSS, which reproducibly occurring in three biological replicates. This questions the hypothesis of experimental noise or pervasive transcription. Accordingly, conserved promoter motifs were found upstream indicating proper TSSs. More than 50% of 5,567 canonical and between 32% and 47% of 10,355 non-canonical TSS were differentially expressed in different media and growth phases, providing evidence for a potential biological function also of non-canonical TSS. Thus, reproducible and environmentally regulated expression suggests that a substantial number of the non-canonical TSSs may be of unknown function rather than being the result of noise or pervasive transcription.

Biological factors in the synthetic construction of overlapping genes.

BACKGROUND: Overlapping genes (OLGs) with long protein-coding overlapping sequences are disallowed by standard genome annotation programs, outside of viruses. Recently however they have been discovered in Archaea, diverse Bacteria, and Mammals. The biological factors underlying life's ability to create overlapping genes require more study, and may have important applications in understanding evolution and in biotechnology. A previous study claimed that protein domains from viruses were much better suited to forming overlaps than those from other cellular organisms - in this study we assessed this claim, in order to discover what might underlie taxonomic differences in the creation of gene overlaps. RESULTS: After overlapping arbitrary Pfam domain pairs and evaluating them with Hidden Markov Models we find OLG construction to be much less constrained than expected. For instance, close to 10% of the constructed sequences cannot be distinguished from typical sequences in their protein family. Most are also indistinguishable from natural protein sequences regarding identity and secondary structure. Surprisingly, contrary to a previous study, virus domains were much less suitable for designing OLGs than bacterial or eukaryotic domains were. In general, the amount of amino acid change required to force a domain to overlap is approximately equal to the variation observed within a typical domain family. The resulting high similarity between natural sequences and those altered so as to overlap is mostly due to the combination of high redundancy in the genetic code and the evolutionary exchangeability of many amino acids. CONCLUSIONS: Synthetic overlapping genes which closely resemble natural gene sequences, as measured by HMM profiles, are remarkably easy to construct, and most arbitrary domain pairs can be altered so as to overlap while retaining high similarity to the original sequences. Future work however will need to assess important factors not considered such as intragenic interactions which affect protein folding. While the analysis here is not sufficient to guarantee functional folding proteins, further analysis of constructed OLGs will improve our understanding of the origin of these remarkable genetic elements across life and opens up exciting possibilities for synthetic biology.

Facklamia lactis sp. nov., isolated from raw milk.

Two strains of a Gram-staining-positive species were isolated from German bulk tank milk. On the basis of their 16S rRNA sequences they were affiliated to the genus Facklamia but could not be assigned to any species with a validly published name. Facklamia miroungae ATCC BAA-466T (97.3 % 16S rRNA sequence similarity), Facklamia languida CCUG 37842T (96.9 %), and Facklamia hominis CCUG 36813T (96.6 %) are the closest relatives. In the 16S rRNA phylogeny and in the core-genome phylogeny strains WS 5301T and WS 5302 form a well-supported, separate lineage. Pairwise average nucleotide identity calculated using MUMmer (ANIm) between WS 5301T and type strains of other Facklamia species is well below the species cut-off (95 %) and ranges from 83.4 to 87.7 %. The DNA G+C content of the type strain is 36.4 mol% and the assembly size of the genome is 2.2 Mb. Cells of WS 5301T are non-motile, non-endospore-forming, oxidase-negative, catalase-negative and facultatively anaerobic cocci. The fastidious species grows at 10-40 °C and with up to 7.0 % (w/v) NaCl in BHI supplemented with 5 g l-1 yeast extract. Major polar lipids are phosphatidylglycerol, diphosphatidylglycerol and two glycolipids. Predominant fatty acids are C16 : 1ω9c and C18 : 1ω9c. On the basis of their genomic, physiological and chemotaxonomic characteristics the strains examined in this study represent the same, hitherto unknown species. We propose the name Facklamia lactis sp. nov. for which WS 5301T (=DSM 111018T=LMG 31861T) is the type strain and WS 5302 (=DSM 111019=LMG 31862) is an additional strain of this novel species.

Pseudomonas cremoris sp. nov., a novel proteolytic species isolated from cream.

During a study investigating the microbiota of raw milk and its semi-finished products, strains WS 5106T and WS 5096 were isolated from cream and skimmed milk concentrate. They could be assigned to the genus Pseudomonas by their 16S rRNA sequences, but not to any validly named species. In this work, a polyphasic approach was used to characterize the novel strains and to investigate their taxonomic status. Examinations based on the topology of core genome phylogenomy as well as average nucleotide identity (ANIm) comparisons suggested a novel Pseudomonas species within the Pseudomonas fluorescens subgroup. With pairwise ANIm values of 90.1 and 89.8 %, WS 5106T was most closely related to Pseudomonas nabeulensis CECT 9765T and Pseudomonas kairouanensis CECT 9766T. The G+C content of strain WS 5106T was 60.1 mol%. Morphologic analyses revealed Gram-stain-negative, aerobic, catalase and oxidase positive, rod-shaped and motile cells. Proteolysis on skimmed milk agar as well as lipolysis on tributyrin agar occurred at both 28 and 6 °C. Tolerated growth conditions were temperatures between 4 and 34 °C, pH values between 6.0 and 8.0, and salt concentrations of up to 5 %. Fatty acid profiles showed a pattern typical for Pseudomonas, with C16 : 0 as the dominant component. The major cellular polar lipids were phosphatidylethanolamine, phosphatidylglycerol and diphosphatidylglycerol and the dominating quinone was Q-9. Based on these results, it is proposed to classify the strains as a novel species, Pseudomonas cremoris sp. nov., with WS 5106T (=DSM 111143T=LMG 31863T) as type strain and WS 5096 (=DSM 111129=LMG 31864) as an additional strain.

Simultaneous quantification of the most common and proteolytic Pseudomonas species in raw milk by multiplex qPCR.

The heat-stable peptidase AprX, secreted by psychrotolerant Pseudomonas species in raw milk, is a major cause of destabilization and premature spoilage of ultra-high temperature (UHT) milk and milk products. To enable rapid detection and quantification of seven frequent and proteolytic Pseudomonas species (P. proteolytica, P. gessardii, P. lactis, P. fluorescens, P. protegens, P. lundensis, and P. fragi) in raw milk, we developed two triplex qPCR assays taking into account species-dependent differences in AprX activity. Besides five species-specific hydrolysis probes, targeting the aprX gene, a universal rpoB probe was included in the assay to determine the total Pseudomonas counts. For all six probes, linear regression lines between Cq value and target DNA concentration were obtained in singleplex as well as in multiplex approaches, yielding R2 values of > 0.975 and amplification efficiencies of 85-97%. Moreover, high specificity was determined using genomic DNA of 75 Pseudomonas strains, assigned to 57 species, and 40 other bacterial species as templates in the qPCR. Quantification of the target species and total Pseudomonas counts resulted in linear detection ranges of approx. 103-107 cfu/ml, which correspond well to common Pseudomonas counts in raw milk. Application of the assay using 60 raw milk samples from different dairies showed good agreement of total Pseudomonas counts calculated by qPCR with cell counts derived from cultivation. Furthermore, a remarkably high variability regarding the species composition was observed for each milk sample, whereby P. lundensis and P. proteolytica/P. gessardii were the predominant species detected. KEY POINTS: • Multiplex qPCR for quantification of seven proteolytic Pseudomonas species and total Pseudomonas counts in raw milk • High specificity and sensitivity via hydrolysis probes against aprX and rpoB • Rapid method to determine Pseudomonas contamination in raw milk and predict spoilage potential.

Amplicon-sequencing of raw milk microbiota: impact of DNA extraction and library-PCR.

The highly complex raw milk matrix challenges the sample preparation for amplicon-sequencing due to low bacterial counts and high amounts of eukaryotic DNA originating from the cow. In this study, we optimized the extraction of bacterial DNA from raw milk for microbiome analysis and evaluated the impact of cycle numbers in the library-PCR. The selective lysis of eukaryotic cells by proteinase K and digestion of released DNA before bacterial lysis resulted in a high reduction of mostly eukaryotic DNA and increased the proportion of bacterial DNA. Comparative microbiome analysis showed that a combined enzymatic and mechanical lysis procedure using the DNeasy® PowerFood® Microbial Kit with a modified protocol was best suitable to achieve high DNA quantities after library-PCR and broad coverage of detected bacterial biodiversity. Increasing cycle numbers during library-PCR systematically altered results for species and beta-diversity with a tendency to overrepresentation or underrepresentation of particular taxa. To limit PCR bias, high cycle numbers should thus be avoided. An optimized DNA extraction yielding sufficient bacterial DNA and enabling higher PCR efficiency is fundamental for successful library preparation. We suggest that a protocol using ethylenediaminetetraacetic acid (EDTA) to resolve casein micelles, selective lysis of somatic cells, extraction of bacterial DNA with a combination of mechanical and enzymatic lysis, and restriction of PCR cycles for analysis of raw milk microbiomes is optimal even for samples with low bacterial numbers. KEY POINTS: • Sample preparation for high-throughput 16S rRNA gene sequencing of raw milk microbiota. • Reduction of eukaryotic DNA by enzymatic digestion. • Shift of detected microbiome caused by high cycle numbers in library-PCR.

Fundicoccus ignavus gen. nov., sp. nov., a novel genus of the family Aerococcaceae isolated from bulk tank milk.

Three strains of a Gram-stain-positive, catalase-negative, facultative anaerobic, and coccoid species were isolated from German bulk tank milk. Phylogenetic analyses based on the 16S rRNA gene sequences indicated that the three strains (WS4937T, WS4759 and WS5303) constitute an independent phylogenetic lineage within the family Aerococcaceae with Facklamia hominis CCUG 36813T (93.7-94.1 %) and Eremococcus coleocola M1831/95/2T (93.5 %) as most closely related type species. The unclassified strains demonstrated variable growth with 6.5 % (w/v) NaCl and tolerated pH 6.5-9.5. Growth was observed from 12 to 39 °C. Their cell-wall peptidoglycan belongs to the A1α type (l-Lys-direct) consisting of alanine, glutamic acid and lysine. The predominant fatty acids were C16 : 1 ω9c, C16 : 0 and C18 : 1 ω9c and in the polar lipids profile three glycolipids, a phospholipid, phosphatidylglycerol, phosphoglycolipid and diphosphatidylglycerol were found. The G+C content of strain WS4937T was 37.4 mol% with a genome size of ~3.0 Mb. Based on phylogenetic, phylogenomic and biochemical characterizations, the isolates can be demarcated from all other genera of the family Aerococcaceae and, therefore, the novel genus Fundicoccus gen. nov. is proposed. The type species of the novel genus is Fundicoccus ignavus gen. nov., sp. nov. WS4937T (=DSM 109652T=LMG 31441T).

Pseudomonas haemolytica sp. nov., isolated from raw milk and skimmed milk concentrate.

Two strains, WS 5063T and WS 5067, isolated from raw cow's milk and skimmed milk concentrate, could be affiliated as members of the same, hitherto unknown, Pseudomonas species by 16S rRNA and rpoD gene sequences. Multilocus sequence and average nucleotide identity (ANIm) analyses based on draft genome sequences confirmed the discovery of a novel Pseudomonas species. It was most closely related to Pseudomonas synxantha DSM 18928T with an ANIm of 91.4 %. The DNA G+C content of WS 5063T was 60.0 mol %. Phenotypic characterizations showed that the isolates are rod-shaped, motile, catalase- and oxidase-positive, and aerobic. Growth occurred at 4-34 °C and at pH values of pH 5.5-8.0. Both strains showed strong ß-haemolysis on blood agar. The major cellular polar lipids were phosphatidylethanolamine, phosphatidylglycerol and diphosphatidylglycerol. The dominant quinone was Q-9 (90 %), but noticeable amounts of Q-8 (9 %) and traces of Q-7 were also detected. Fatty acid profiles were typical for Pseudomonas species and exhibited C16 : 0 as a major component. Based on these results, we conclude that both strains belong to a novel species, for which the name Pseudomonas haemolytica sp. nov. is proposed. The type strain is WS 5063T (=DSM 108987T=LMG 31232T) and an additional strain is WS 5067 (=DSM 108988=LMG 31233).

Pseudomonas saxonica sp. nov., isolated from raw milk and skimmed milk concentrate.

A polyphasic approach was used to investigate the taxonomic status of two bacterial strains, WS 5072T and WS 5092, isolated from skimmed milk concentrate and raw cow's milk. The 16S rRNA and rpoD gene sequences affiliated the strains to the same, hitherto unknown, Pseudomonas species. Further examinations of the draft genomes based on multilocus sequence analysis and average nucleotide identity confirmed the presence of a novel Pseudomonas species. It was most closely related to Pseudomonas fragi DSM 3456T with 86.3 % ANIm. The DNA G+C content of strain WS 5072T was 56.3 mol%. Cells were aerobic, Gram-negative, catalase and oxidase positive, rod-shaped and motile. Growth occurred at 4-34 °C, pH 5.5-8.0 and with salt concentrations of up to 7 %. The major cellular polar lipids were phosphatidylethanolamine, phosphatidylglycerol, and diphosphatidylglycerol. The dominating quinone was Q-9 with 94 %, with noticeable amounts of Q-8 (5 %) and traces of Q-7 and Q-10. Fatty acid profiles showed a composition common for Pseudomonas with the major component C16 : 0. Based on these results, the novel species Pseudomonas saxonica sp. nov. is proposed, with the type strain WS 5072T (=DSM 108989T=LMG 31234T) and the additional strain WS 5092 (=DSM 108990=LMG 31235).

Brevilactibacter flavus gen. nov., sp. nov., a novel bacterium of the family Propionibacteriaceae isolated from raw milk and dairy products and reclassification of Propioniciclava sinopodophylli as Brevilactibacter sinopodophylli comb. nov.

Eight facultatively anaerobic rod-shaped bacteria were isolated from raw milk and two other dairy products. Results of phylogenetic analyses based on 16S rRNA gene sequences showed that the isolates are placed in a distinct lineage within the family Propionibacteriaceae with Propioniciclava sinopodophylli and Propioniciclava tarda as the closest relatives (94.6 and 93.5 % similarity, respectively). The cell-wall peptidoglycan contained meso-diaminopimelic acid, alanine and glutamic acid and was of the A1γ type (meso-DAP-direct). The major cellular fatty acid was anteiso-C15 : 0 and the major polar lipids were diphosphatidylglycerol, phosphatidyglycerol and three unidentified glycolipids. The quinone system contained predominantly menaquinone MK-9(H4). The G+C content of the genomic DNA of strain VG341T was 67.7 mol%. The whole-cell sugar pattern contained ribose, rhamnose, arabinose and galactose. On the basis of phenotypic and genetic data, eight strains (VG341T, WS4684, WS4769, WS 4882, WS4883, WS4901, WS4902 and WS4904) are proposed to be classified as members of a novel species in a new genus of the family Propionibacteriaceae, for which the name Brevilactibacter flavus gen. nov., sp. nov. is proposed. The type strain is VG341T (=WS4900T=DSM 100885T=LMG 29089T) and seven additional strains are WS4684, WS4769, WS4882, WS4883, WS4901, WS4902 and WS4904. Furthermore, we propose the reclassification of P. sinopodophylli as Brevilactibacter sinopodophylli comb. nov.

Resistance of thermophilic spore formers isolated from milk and whey products towards cleaning-in-place conditions: Influence of pH, temperature and milk residues.

The occurrence of thermophilic spore formers in dairy powders is a major concern for producers worldwide. This study aims to investigate the resistance of thermophilic endospores towards cleaning solutions typically used for cleaning-in-place in dairy manufacturing plants. From eleven tested strains, all were able to survive an alkaline treatment (NaOH) at 65 °C for 10 min (0.5%), whereas at concentrations of 2% eight strains withstood the treatment. Acid solutions were more sporicidal. At 0.5% of HNO3, only three strains survived the treatment. Milk impurities reduced the inactivation effect of the NaOH solutions towards thermophilic spore formers. For two selected strains, a detailed kinetic inactivation in NaOH and HNO3 solutions at different temperatures was performed and non-log-linear inactivation curves were observed. This study highlights the risk of reusing cleaning solutions in dairies.

A novel short L-arginine responsive protein-coding gene (laoB) antiparallel overlapping to a CadC-like transcriptional regulator in Escherichia coli O157:H7 Sakai originated by overprinting.

BACKGROUND: Due to the DNA triplet code, it is possible that the sequences of two or more protein-coding genes overlap to a large degree. However, such non-trivial overlaps are usually excluded by genome annotation pipelines and, thus, only a few overlapping gene pairs have been described in bacteria. In contrast, transcriptome and translatome sequencing reveals many signals originated from the antisense strand of annotated genes, of which we analyzed an example gene pair in more detail. RESULTS: A small open reading frame of Escherichia coli O157:H7 strain Sakai (EHEC), designated laoB (L-arginine responsive overlapping gene), is embedded in reading frame -2 in the antisense strand of ECs5115, encoding a CadC-like transcriptional regulator. This overlapping gene shows evidence of transcription and translation in Luria-Bertani (LB) and brain-heart infusion (BHI) medium based on RNA sequencing (RNAseq) and ribosomal-footprint sequencing (RIBOseq). The transcriptional start site is 289 base pairs (bp) upstream of the start codon and transcription termination is 155 bp downstream of the stop codon. Overexpression of LaoB fused to an enhanced green fluorescent protein (EGFP) reporter was possible. The sequence upstream of the transcriptional start site displayed strong promoter activity under different conditions, whereas promoter activity was significantly decreased in the presence of L-arginine. A strand-specific translationally arrested mutant of laoB provided a significant growth advantage in competitive growth experiments in the presence of L-arginine compared to the wild type, which returned to wild type level after complementation of laoB in trans. A phylostratigraphic analysis indicated that the novel gene is restricted to the Escherichia/Shigella clade and might have originated recently by overprinting leading to the expression of part of the antisense strand of ECs5115. CONCLUSIONS: Here, we present evidence of a novel small protein-coding gene laoB encoded in the antisense frame -2 of the annotated gene ECs5115. Clearly, laoB is evolutionarily young and it originated in the Escherichia/Shigella clade by overprinting, a process which may cause the de novo evolution of bacterial genes like laoB.

Reprogramming Human Siderocalin To Neutralize Petrobactin, the Essential Iron Scavenger of Anthrax Bacillus.

Bacillus anthracis owes its pronounced virulence-apart from specific toxins-to a twofold import mechanism for FeIII ions. This pathogenic bacterium secretes the siderophores bacillibactin (BB) and petrobactin (PB), of which only BB is neutralized by human siderocalin, an abundant lipocalin in plasma. We describe its reshaping via combinatorial protein design to bind PB⋅FeIII instead of BB⋅FeIII , and with even higher affinity (KD ≈20 pm). X-ray crystallographic analysis of the resulting "petrocalin" in complex with PB⋅GaIII reveals a positively charged ligand pocket while the extended butterfly-like conformation of the bound PB provides a rationale for the missing recognition by the natural siderocalin. In microbiological studies, a combination of petrocalin and siderocalin effectively suppressed the growth of a BB+ /PB+ strain of Bacillus cereus under iron-limiting culture conditions. Thus, our reprogrammed lipocalin may offer novel treatment options for devastating infections caused by B. anthracis.

Differentiation of ncRNAs from small mRNAs in Escherichia coli O157:H7 EDL933 (EHEC) by combined RNAseq and RIBOseq - ryhB encodes the regulatory RNA RyhB and a peptide, RyhP.

BACKGROUND: While NGS allows rapid global detection of transcripts, it remains difficult to distinguish ncRNAs from short mRNAs. To detect potentially translated RNAs, we developed an improved protocol for bacterial ribosomal footprinting (RIBOseq). This allowed distinguishing ncRNA from mRNA in EHEC. A high ratio of ribosomal footprints per transcript (ribosomal coverage value, RCV) is expected to indicate a translated RNA, while a low RCV should point to a non-translated RNA. RESULTS: Based on their low RCV, 150 novel non-translated EHEC transcripts were identified as putative ncRNAs, representing both antisense and intergenic transcripts, 74 of which had expressed homologs in E. coli MG1655. Bioinformatics analysis predicted statistically significant target regulons for 15 of the intergenic transcripts; experimental analysis revealed 4-fold or higher differential expression of 46 novel ncRNA in different growth media. Out of 329 annotated EHEC ncRNAs, 52 showed an RCV similar to protein-coding genes, of those, 16 had RIBOseq patterns matching annotated genes in other enterobacteriaceae, and 11 seem to possess a Shine-Dalgarno sequence, suggesting that such ncRNAs may encode small proteins instead of being solely non-coding. To support that the RIBOseq signals are reflecting translation, we tested the ribosomal-footprint covered ORF of ryhB and found a phenotype for the encoded peptide in iron-limiting condition. CONCLUSION: Determination of the RCV is a useful approach for a rapid first-step differentiation between bacterial ncRNAs and small mRNAs. Further, many known ncRNAs may encode proteins as well.

Pseudomonas lactis sp. nov. and Pseudomonas paralactis sp. nov., isolated from bovine raw milk.

Five strains, designated WS 4672T, WS 4998, WS 4992T, WS 4997 and WS 5000, isolated from bovine raw milk formed two individual groups in a phylogenetic analysis. The most similar species on the basis of 16S rRNA gene sequences were Pseudomonas azotoformans IAM 1603T, Pseudomonas gessardii CIP 105469T and Pseudomonas libanensis CIP 105460T showing 99.7-99.6 % similarity. Using rpoD gene sequences Pseudomonas veronii LMG 17761T (93.3 %) was most closely related to strain WS 4672T and Pseudomonas libanensis CIP 105460T to strain WS 4992T (93.3 %). The five strains could be differentiated from their closest relatives and from each other by phenotypic and chemotaxonomic characterization and ANIb values calculated from draft genome assemblies. ANIb values of strains WS 4992T and WS4671T to the closest relatives are lower than 90 %. The major cellular polar lipids of both strains are phosphatidylethanolamine, phosphatidylglycerol, a phospholipid and diphosphatidylglycerol, and their major quinone is Q-9. The DNA G+C content of strains WS 4992T and WS 4672T were 60.0  and 59.7  mol%, respectively. Based on these genotypic and phenotypic traits two novel species of the genus Pseudomonas are proposed: Pseudomonas lactis sp. nov. [with type strain WS 4992T (=DSM 29167T=LMG 28435T) and the additional strains WS 4997 and WS 5000], and Pseudomonasparalactis sp. nov. [with type strain WS 4672T (=DSM 29164T=LMG 28439T) and additional strain WS 4998].

Translatomics combined with transcriptomics and proteomics reveals novel functional, recently evolved orphan genes in Escherichia coli O157:H7 (EHEC).

BACKGROUND: Genomes of E. coli, including that of the human pathogen Escherichia coli O157:H7 (EHEC) EDL933, still harbor undetected protein-coding genes which, apparently, have escaped annotation due to their small size and non-essential function. To find such genes, global gene expression of EHEC EDL933 was examined, using strand-specific RNAseq (transcriptome), ribosomal footprinting (translatome) and mass spectrometry (proteome). RESULTS: Using the above methods, 72 short, non-annotated protein-coding genes were detected. All of these showed signals in the ribosomal footprinting assay indicating mRNA translation. Seven were verified by mass spectrometry. Fifty-seven genes are annotated in other enterobacteriaceae, mainly as hypothetical genes; the remaining 15 genes constitute novel discoveries. In addition, protein structure and function were predicted computationally and compared between EHEC-encoded proteins and 100-times randomly shuffled proteins. Based on this comparison, 61 of the 72 novel proteins exhibit predicted structural and functional features similar to those of annotated proteins. Many of the novel genes show differential transcription when grown under eleven diverse growth conditions suggesting environmental regulation. Three genes were found to confer a phenotype in previous studies, e.g., decreased cattle colonization. CONCLUSIONS: These findings demonstrate that ribosomal footprinting can be used to detect novel protein coding genes, contributing to the growing body of evidence that hypothetical genes are not annotation artifacts and opening an additional way to study their functionality. All 72 genes are taxonomically restricted and, therefore, appear to have evolved relatively recently de novo.

Permanent colonization of creek sediments, creek water and limnic water plants by four Listeria species in low population densities.

During a 1-year longitudinal study, water, sediment and water plants from two creeks and one pond were sampled monthly and analyzed for the presence of Listeria species. A total of 90 % of 30 sediment samples, 84 % of 31 water plant samples and 67 % of 36 water samples were tested positive. Generally, most probable number counts ranged between 1 and 40 g-1, only occasionally >110 cfu g-1 were detected. Species differentiation based on FT-IR spectroscopy and multiplex PCR of a total of 1220 isolates revealed L. innocua (46 %), L. seeligeri (27 %), L. monocytogenes (25 %) and L. ivanovii (2 %). Titers and species compositions were similar during all seasons. While the species distributions in sediments and associated Ranunculus fluitans plants appeared to be similar in both creeks, RAPD typing did not provide conclusive evidence that the populations of these environments were connected. It is concluded that (i) the fresh-water sediments and water plants are year-round populated by Listeria, (ii) no clear preference for growth in habitats as different as sediments and water plants was found and (iii) the RAPD-based intraspecific biodiversity is high compared to the low population density.

Massive horizontal gene transfer, strictly vertical inheritance and ancient duplications differentially shape the evolution of Bacillus cereus enterotoxin operons hbl, cytK and nhe.

BACKGROUND: Bacillus cereus sensu lato comprises eight closely related species including the human pathogens Bacillus anthracis and Bacillus cereus. Within B. cereus sensu lato, chromosomally and plasmid-encoded toxins exist. While plasmid-mediated horizontal gene transfer of the emetic toxin, anthrax and insecticidal toxins is known, evolution of enterotoxin genes within the group has not been studied. RESULTS: We report draft genome assemblies of 25 strains, a phylogenetic network of 142 strains based on ANI derived from genome sequences and a phylogeny based on whole-genome SNP analysis. The data clearly support subdivision of B. cereus sensu lato into seven phylogenetic groups. While group I, V and VII represent B. pseudomycoides, B. toyonensis and B. cytotoxicus, which are distinguishable at species level (ANI border ≥ 96 %), strains ascribed to the other five species do not match phylogenic groups. The chromosomal enterotoxin operons nheABC and hblCDAB are abundant within B. cereus both isolated from infections and from the environment. While the duplicated hbl variant hbl a is present in 22 % of all strains investigated, duplication of nheABC is extremely rare (0.02 %) and appears to be phylogenetically unstable. Distribution of toxin genes was matched to a master tree based on seven concatenated housekeeping genes, which depicts species relationships in B. cereus sensu lato as accurately as whole-genome comparisons. Comparison to the phylogeny of enterotoxin genes uncovered ample evidence for horizontal transfer of hbl, cytK and plcR, as well as frequent deletion of both toxins and duplication of hbl. No evidence for nhe deletion was found and stable horizontal transfer of nhe is rare. Therefore, evolution of B. cereus enterotoxin operons is shaped unexpectedly different for yet unknown reasons. CONCLUSIONS: Frequent exchange of the pathogenicity factors hbl, cytK and plcR in B. cereus sensu lato appears to be an important mechanism of B. cereus virulence evolution, including so-called probiotic or non-pathogenic species, which might have consequences for risk assessment procedures. In contrast, exclusively vertical inheritance of nhe was observed, and since nhe-negative strains appear to be extremely rare, we suggest that fitness loss may be associated with deletion or horizontal transfer of the nhe operon.

Evidence for the recent origin of a bacterial protein-coding, overlapping orphan gene by evolutionary overprinting.

BACKGROUND: Gene duplication is believed to be the classical way to form novel genes, but overprinting may be an important alternative. Overprinting allows entirely novel proteins to evolve de novo, i.e., formerly non-coding open reading frames within functional genes become expressed. Only three cases have been described for Escherichia coli. Here, a fourth example is presented. RESULTS: RNA sequencing revealed an open reading frame weakly transcribed in cow dung, coding for 101 residues and embedded completely in the -2 reading frame of citC in enterohemorrhagic E. coli. This gene is designated novel overlapping gene, nog1. The promoter region fused to gfp exhibits specific activities and 5' rapid amplification of cDNA ends indicated the transcriptional start 40-bp upstream of the start codon. nog1 was strand-specifically arrested in translation by a nonsense mutation silent in citC. This Nog1-mutant showed a phenotype in competitive growth against wild type in the presence of MgCl2. Small differences in metabolite concentrations were also found. Bioinformatic analyses propose Nog1 to be inner membrane-bound and to possess at least one membrane-spanning domain. A phylogenetic analysis suggests that the orphan gene nog1 arose by overprinting after Escherichia/Shigella separated from the other γ-proteobacteria. CONCLUSIONS: Since nog1 is of recent origin, non-essential, short, weakly expressed and only marginally involved in E. coli's central metabolism, we propose that this gene is in an initial stage of evolution. While we present specific experimental evidence for the existence of a fourth overlapping gene in enterohemorrhagic E. coli, we believe that this may be an initial finding only and overlapping genes in bacteria may be more common than is currently assumed by microbiologists.