Two novel genes identified by large-scale transcriptomic analysis are essential for biofilm and rugose colony development of Vibrio vulnificus.

Many pathogenic bacteria form biofilms to survive under environmental stresses and host immune defenses. Differential expression (DE) analysis of the genes in biofilm and planktonic cells under a single condition, however, has limitations to identify the genes essential for biofilm formation. Independent component analysis (ICA), a machine learning algorithm, was adopted to comprehensively identify the biofilm genes of Vibrio vulnificus, a fulminating human pathogen, in this study. ICA analyzed the large-scale transcriptome data of V. vulnificus cells under various biofilm and planktonic conditions and then identified a total of 72 sets of independently co-regulated genes, iModulons. Among the three iModulons specifically activated in biofilm cells, BrpT-iModulon mainly consisted of known genes of the regulon of BrpT, a transcriptional regulator controlling biofilm formation of V. vulnificus. Interestingly, the BrpT-iModulon additionally contained two novel genes, VV1_3061 and VV2_1694, designated as cabH and brpN, respectively. cabH and brpN were shared in other Vibrio species and not yet identified by DE analyses. Genetic and biochemical analyses revealed that cabH and brpN are directly up-regulated by BrpT. The deletion of cabH and brpN impaired the robust biofilm and rugose colony formation. CabH, structurally similar to the previously known calcium-binding matrix protein CabA, was essential for attachment to the surface. BrpN, carrying an acyltransferase-3 domain as observed in BrpL, played an important role in exopolysaccharide production. Altogether, ICA identified two novel genes, cabH and brpN, which are regulated by BrpT and essential for the development of robust biofilms and rugose colonies of V. vulnificus.

Independent Component Analysis Identifies the Modulons Expanding the Transcriptional Regulatory Networks of Enterohemorrhagic Escherichia coli.

The elucidation of the transcriptional regulatory networks (TRNs) of enterohemorrhagic Escherichia coli (EHEC) is critical to understand its pathogenesis and survival in the host. However, the analyses of current TRNs are still limited to comprehensively understand their target genes generally co-regulated under various conditions regardless of the genetic backgrounds. In this study, independent component analysis (ICA), a machine learning-based decomposition method, was used to decompose the large-scale transcriptome data of EHEC into the modulons, which contain the target genes of several TRNs. The locus of enterocyte effacement (LEE) and the Shiga toxin (Stx) modulons mainly consisted of the Ler regulon and the Stx prophage genes, respectively, confirming that ICA properly grouped the co-regulated major virulence genes of EHEC. Further investigation revealed that the LEE modulon contained the hypothetical Z0395 gene as a novel member of the Ler regulon, and the Stx modulon contained the thi and cus locus genes in addition to the Stx prophage genes. Correspondingly, the Stx prophage genes were also regulated by thiamine and copper ions known to control the thi and cus locus genes, respectively. The modulons effectively clustered the genes co-regulated regardless of the growth conditions and the genetic backgrounds of EHEC. The changed activities of the individual modulons successfully explained the differential expressions of the virulence and survival genes during the course of infection in bovines. Altogether, these results suggested that ICA of the large-scale transcriptome data can expand and enhance the current understanding of the TRNs of EHEC.

A Master Regulator BrpR Coordinates the Expression of Multiple Loci for Robust Biofilm and Rugose Colony Development in Vibrio vulnificus.

Vibrio vulnificus, a fulminating human pathogen, forms biofilms to enhance its survival in nature and pathogenicity during host infection. BrpR is the transcriptional regulator governing robust biofilm and rugose colony formation in V. vulnificus, but little is known about both the direct regulon of BrpR and the role of BrpR in regulation of downstream genes. In this study, transcript analyses revealed that BrpR is highly expressed and thus strongly regulates the downstream gene in the stationary and elevated cyclic di-GMP conditions. Transcriptome analyses discovered the genes, whose expression is affected by BrpR but not by the downstream regulator BrpT. Two unnamed adjacent genes (VV2_1626-1627) were newly identified among the BrpR regulon and designated as brpL and brpG in this study. Genetic analyses showed that the deletion of brpL and brpG impairs the biofilm and rugose colony formation, indicating that brpLG plays a crucial role in the development of BrpR-regulated biofilm phenotypes. Comparison of the colony morphology and exopolysaccharide (EPS) production suggested that although the genetic location and regulation of brpLG are distinct from the brp locus, brpABCDFHIJK (VV2_1574-1582), brpLG is also responsible for the robust EPS production together with the brp locus genes. Electrophoretic mobility shift assays and DNase I protection assays demonstrated that BrpR regulates the expression of downstream genes in distinct loci by directly binding to their upstream regions, revealing a palindromic binding sequence. Altogether, this study suggests that BrpR is a master regulator coordinating the expression of multiple loci responsible for EPS production and thus, contributing to the robust biofilm and rugose colony formation of V. vulnificus.

A Nitric Oxide-Responsive Transcriptional Regulator NsrR Cooperates With Lrp and CRP to Tightly Control the hmpA Gene in Vibrio vulnificus.

Nitric oxide (NO) is an important antimicrobial effector produced by the host innate immune system to counteract invading pathogens. To survive and establish a successful infection, a fulminating human pathogen Vibrio vulnificus expresses the hmpA gene encoding an NO dioxygenase in an NO-responsive manner. In this study, we identified an Rrf2-family transcriptional regulator NsrR that is predicted to contain the Fe-S cluster coordinated by three cysteine residues. Transcriptome analysis showed that NsrR controls the expression of multiple genes potentially involved in nitrosative stress responses. Particularly, NsrR acts as a strong repressor of hmpA transcription and relieves the repression of hmpA upon exposure to NO. Notably, nsrR and hmpA are transcribed divergently, and their promoter regions overlap with each other. Molecular biological analyses revealed that NsrR directly binds to this overlapping promoter region, which is alleviated by loss of the Fe-S cluster, leading to the subsequent derepression of hmpA under nitrosative stress. We further found that a leucine-responsive regulatory protein (Lrp) negatively regulates hmpA in an NsrR-dependent manner by directly binding to the promoter region, presumably resulting in a DNA conformation change to support the repression by NsrR. Meanwhile, a cyclic AMP receptor protein (CRP) positively regulates hmpA probably through repression of nsrR and lrp by directly binding to each promoter region in a sequential cascade. Altogether, this collaborative regulation of NsrR along with Lrp and CRP enables an elaborate control of hmpA transcription, contributing to survival under host-derived nitrosative stress and thereby the pathogenesis of V. vulnificus.

Pathogenic potential assessment of the Shiga toxin-producing Escherichia coli by a source attribution-considered machine learning model.

Instead of conventional serotyping and virulence gene combination methods, methods have been developed to evaluate the pathogenic potential of newly emerging pathogens. Among them, the machine learning (ML)-based method using whole-genome sequencing (WGS) data are getting attention because of the recent advances in ML algorithms and sequencing technologies. Here, we developed various ML models to predict the pathogenicity of Shiga toxin-producing Escherichia coli (STEC) isolates using their WGS data. The input dataset for the ML models was generated using distinct gene repertoires from positive (pathogenic) and negative (nonpathogenic) control groups in which each STEC isolate was designated based on the source attribution, the relative risk potential of the isolation sources. Among the various ML models examined, a model using the support vector machine (SVM) algorithm, the SVM model, discriminated between the two control groups most accurately. The SVM model successfully predicted the pathogenicity of the isolates from the major sources of STEC outbreaks, the isolates with the history of outbreaks, and the isolates that cannot be assessed by conventional methods. Furthermore, the SVM model effectively differentiated the pathogenic potentials of the isolates at a finer resolution. Permutation importance analyses of the input dataset further revealed the genes important for the estimation, proposing the genes potentially essential for the pathogenicity of STEC. Altogether, these results suggest that the SVM model is a more reliable and broadly applicable method to evaluate the pathogenic potential of STEC isolates compared with conventional methods.

The transcriptional regulator IscR integrates host-derived nitrosative stress and iron starvation in activation of the vvhBA operon in Vibrio vulnificus.

For successful infection of their hosts, pathogenic bacteria recognize host-derived signals that induce the expression of virulence factors in a spatiotemporal manner. The fulminating food-borne pathogen Vibrio vulnificus produces a cytolysin/hemolysin protein encoded by the vvhBA operon, which is a virulence factor preferentially expressed upon exposure to murine blood and macrophages. The Fe-S cluster containing transcriptional regulator IscR activates the vvhBA operon in response to nitrosative stress and iron starvation, during which the cellular IscR protein level increases. Here, electrophoretic mobility shift and DNase I protection assays revealed that IscR directly binds downstream of the vvhBA promoter P vvhBA , which is unusual for a positive regulator. We found that in addition to IscR, the transcriptional regulator HlyU activates vvhBA transcription by directly binding upstream of P vvhBA , whereas the histone-like nucleoid-structuring protein (H-NS) represses vvhBA by extensively binding to both downstream and upstream regions of its promoter. Of note, the binding sites of IscR and HlyU overlapped with those of H-NS. We further substantiated that IscR and HlyU outcompete H-NS for binding to the P vvhBA regulatory region, resulting in the release of H-NS repression and vvhBA induction. We conclude that concurrent antirepression by IscR and HlyU at regions both downstream and upstream of P vvhBA provides V. vulnificus with the means of integrating host-derived signal(s) such as nitrosative stress and iron starvation for precise regulation of vvhBA transcription, thereby enabling successful host infection.