Conservation of Small Regulatory RNAs in Vibrio parahaemolyticus: Possible role of RNA-OUT Encoded by the Pathogenicity Island (VPaI-7) of Pandemic Strains.

Small regulatory RNAs (sRNAs) are molecules that play an important role in the regulation of gene expression. sRNAs in bacteria can affect important processes, such as metabolism and virulence. Previous studies showed a significant role of sRNAs in the Vibrio species, but knowledge about Vibrio parahaemolyticus is limited. Here, we examined the conservation of sRNAs between V. parahaemolyticus and other human Vibrio species, in addition to investigating the conservation between V. parahaemolyticus strains differing in pandemic origin. Our results showed that only 7% of sRNAs were conserved between V. parahaemolyticus and other species, but 88% of sRNAs were highly conserved within species. Nonetheless, two sRNAs coding to RNA-OUT, a component of the Tn10/IS10 system, were exclusively present in pandemic strains. Subsequent analysis showed that both RNA-OUT were located in pathogenicity island-7 and would interact with transposase VPA1379, according to the model of pairing of IS10-encoded antisense RNAs. According to the location of RNA-OUT/VPA1379, we also investigated if they were expressed during infection. We observed that the transcriptional level of VPA1379 was significantly increased, while RNA-OUT was decreased at three hours post-infection. We suggest that IS10 transcription increases in pandemic strains during infection, probably to favor IS10 transposition and improve their fitness when they are facing adverse conditions.

Gene expression of Vibrio parahaemolyticus growing in laboratory isolation conditions compared to those common in its natural ocean environment.

BACKGROUND: Vibrio parahaemolyticus is an autochthonous marine bacterial species comprising strains able to grow in broth containing bile salts at 37 °C, a condition seldom found in the ocean. However, this condition is used for isolation in the laboratory because it is considered a necessary property for pathogenesis. In this context, revealing how gene expression enables V. parahaemolyticus to adapt to this particular condition -common to almost all V. parahaemolyticus isolates- will improve our understanding of the biology of this important pathogen. To determine the genes of V. parahaemolyticus differentially expressed when growing in isolation condition (37 °C, 0.9% NaCl, and 0.04% bile salts) referred to those at the temperature and salt concentration prevailing in ocean south of Chile (marine-like condition; 12 °C, 3% NaCl, and absence of bile salts) we used high-throughput sequencing of RNA. RESULTS: Our results showed that in the isolation condition, among the 5034 genes annotated in the V. parahaemolyticus RIMD2210633 genome, 344 were upregulated and 433 downregulated referred to the marine-like condition, managing an adjusted P-value (Padj) < E-5. Between the 50 more highly expressed genes, among the small RNAs (sRNA), the three carbon storage regulators B (CsrB) were up four to six times, while RyhB, related to iron metabolism besides motility control, was down about eight times. Among proteins, BfdA, a hemolysin-co-regulated protein (Hcp1) secreted by T6SS1, one of the most highly expressed genes, was about 140 times downregulated in isolation condition. The highest changes in relative expression were found among neighboring genes coding for proteins related to respiration, which were about 40 times upregulated. CONCLUSIONS: When V. parahaemolyticus is grown in conditions used for laboratory isolation 777 genes are up- or downregulated referred to conditions prevailing in the marine-like condition; the most significantly overrepresented categories among upregulated processes were those related to transport and localization, while secretion and pathogenesis were overrepresented among downregulated genes. Genes with the highest differential expression included the sRNAs CsrB and RhyB and the mRNAs related with secretion, nutritional upshift, respiration and rapid growing.

Contribution of the Type III Secretion System (T3SS2) of Vibrio parahaemolyticus in Mitochondrial Stress in Human Intestinal Cells.

Vibrio parahaemolyticus is an important human pathogen that is currently the leading cause of shellfish-borne gastroenteritis in the world. Particularly, the pandemic strain has the capacity to induce cytotoxicity and enterotoxicity through its Type 3 Secretion System (T3SS2) that leads to massive cell death. However, the specific mechanism by which the T3SS2 induces cell death remains unclear and its contribution to mitochondrial stress is not fully understood. In this work, we evaluated the contribution of the T3SS2 of V. parahaemolyticus in generating mitochondrial stress during infection in human intestinal HT-29 cells. To evaluate the contribution of the T3SS2 of V. parahaemolyticus in mitochondrial stress, infection assays were carried out to evaluate mitochondrial transition pore opening, mitochondrial fragmentation, ATP quantification, and cell viability during infection. Our results showed that the Δvscn1 (T3SS2+) mutant strain contributes to generating the sustained opening of the mitochondrial transition pore. Furthermore, it generates perturbations in the ATP production in infected cells, leading to a significant decrease in cell viability and loss of membrane integrity. Our results suggest that the T3SS2 from V. parahaemolyticus plays a role in generating mitochondrial stress that leads to cell death in human intestinal HT-29 cells. It is important to highlight that this study represents the first report indicating the possible role of the V. parahaemolyticus T3SS2 and its effector proteins involvement in generating mitochondrial stress, its impact on the mitochondrial pore, and its effect on ATP production in human cells.

Inclusion Depth for Contour Ensembles.

Ensembles of contours arise in various applications like simulation, computer-aided design, and semantic segmentation. Uncovering ensemble patterns and analyzing individual members is a challenging task that suffers from clutter. Ensemble statistical summarization can alleviate this issue by permitting analyzing ensembles' distributional components like the mean and median, confidence intervals, and outliers. Contour boxplots, powered by Contour Band Depth (CBD), are a popular non-parametric ensemble summarization method that benefits from CBD's generality, robustness, and theoretical properties. In this work, we introduce Inclusion Depth (ID), a new notion of contour depth with three defining characteristics. First, ID is a generalization of functional Half-Region Depth, which offers several theoretical guarantees. Second, ID relies on a simple principle: the inside/outside relationships between contours. This facilitates implementing ID and understanding its results. Third, the computational complexity of ID scales quadratically in the number of members of the ensemble, improving CBD's cubic complexity. This also in practice speeds up the computation enabling the use of ID for exploring large contour ensembles or in contexts requiring multiple depth evaluations like clustering. In a series of experiments on synthetic data and case studies with meteorological and segmentation data, we evaluate ID's performance and demonstrate its capabilities for the visual analysis of contour ensembles.

Accelerating Hyperbolic t-SNE.

The need to understand the structure of hierarchical or high-dimensional data is present in a variety of fields. Hyperbolic spaces have proven to be an important tool for embedding computations and analysis tasks as their non-linear nature lends itself well to tree or graph data. Subsequently, they have also been used in the visualization of high-dimensional data, where they exhibit increased embedding performance. However, none of the existing dimensionality reduction methods for embedding into hyperbolic spaces scale well with the size of the input data. That is because the embeddings are computed via iterative optimization schemes and the computation cost of every iteration is quadratic in the size of the input. Furthermore, due to the non-linear nature of hyperbolic spaces, euclidean acceleration structures cannot directly be translated to the hyperbolic setting. This article introduces the first acceleration structure for hyperbolic embeddings, building upon a polar quadtree. We compare our approach with existing methods and demonstrate that it computes embeddings of similar quality in significantly less time. Implementation and scripts for the experiments can be found at https://graphics.tudelft.nl/accelerating-hyperbolic-tsne.

Presence of Zonula Occludens Toxin-Coding Genes among Vibrio parahaemolyticus Isolates of Clinical and Environmental Origin.

In recent studies, emphasis has been placed on the zonula occludens toxin (Zot) from the non-toxigenic Vibrio parahaemolyticus strain PMC53.7 as an agent inducing alterations in the actin cytoskeleton of infected Caco-2 cells and which appears as a relevant virulence factor. Universal zot primers were designed by the alignment of different types of zot gene and identification of conserved sequences to investigate the presence in diverse environmental and clinical V. parahaemolyticus isolates, in co-occurrence with virulence factors, such a hemolysins and secretion systems. The study screened a total of 390 isolates from environmental sources from Chile and Italy and 95 Chilean clinical isolates. The results revealed that around 37.2% of Chilean environmental strains and 25.9% of Italian strains, and 24.2% of clinical isolates carried the zot gene. The Zot-C2 cluster was present in 71.4% of Chilean environmental strains but absent in clinical isolates, while the Zot-C4 cluster was identified in 28.6% of environmental and 100% of clinical isolates. Understanding the role of zot in V. parahaemolyticus virulence is crucial, especially considering the risk associated with consuming diverse isolates from bivalves and the co-occurrence with virulence factors such as TDH, TRH or T3SS2.

Large-scale dose evaluation of deep learning organ contours in head-and-neck radiotherapy by leveraging existing plans.

Background and purpose: Retrospective dose evaluation for organ-at-risk auto-contours has previously used small cohorts due to additional manual effort required for treatment planning on auto-contours. We aimed to do this at large scale, by a) proposing and assessing an automated plan optimization workflow that used existing clinical plan parameters and b) using it for head-and-neck auto-contour dose evaluation. Materials and methods: Our automated workflow emulated our clinic's treatment planning protocol and reused existing clinical plan optimization parameters. This workflow recreated the original clinical plan (POG) with manual contours (PMC) and evaluated the dose effect (POG-PMC) on 70 photon and 30 proton plans of head-and-neck patients. As a use-case, the same workflow (and parameters) created a plan using auto-contours (PAC) of eight head-and-neck organs-at-risk from a commercial tool and evaluated their dose effect (PMC-PAC). Results: For plan recreation (POG-PMC), our workflow had a median impact of 1.0% and 1.5% across dose metrics of auto-contours, for photon and proton respectively. Computer time of automated planning was 25% (photon) and 42% (proton) of manual planning time. For auto-contour evaluation (PMC-PAC), we noticed an impact of 2.0% and 2.6% for photon and proton radiotherapy. All evaluations had a median ΔNTCP (Normal Tissue Complication Probability) less than 0.3%. Conclusions: The plan replication capability of our automated program provides a blueprint for other clinics to perform auto-contour dose evaluation with large patient cohorts. Finally, despite geometric differences, auto-contours had a minimal median dose impact, hence inspiring confidence in their utility and facilitating their clinical adoption.

Vibrio type III secretion system 2 is not restricted to the Vibrionaceae and encodes differentially distributed repertoires of effector proteins.

Vibrio parahaemolyticus is the leading cause of seafood-borne gastroenteritis worldwide. A distinctive feature of the O3:K6 pandemic clone, and its derivatives, is the presence of a second, phylogenetically distinct, type III secretion system (T3SS2) encoded within the genomic island VPaI-7. The T3SS2 allows the delivery of effector proteins directly into the cytosol of infected eukaryotic cells to subvert key host-cell processes, critical for V. parahaemolyticus to colonize and cause disease. Furthermore, the T3SS2 also increases the environmental fitness of V. parahaemolyticus in its interaction with bacterivorous protists; hence, it has been proposed that it contributed to the global oceanic spread of the pandemic clone. Several reports have identified T3SS2-related genes in Vibrio and non-Vibrio species, suggesting that the T3SS2 gene cluster is not restricted to the Vibrionaceae and can mobilize through horizontal gene transfer events. In this work, we performed a large-scale genomic analysis to determine the phylogenetic distribution of the T3SS2 gene cluster and its repertoire of effector proteins. We identified putative T3SS2 gene clusters in 1130 bacterial genomes from 8 bacterial genera, 5 bacterial families and 47 bacterial species. A hierarchical clustering analysis allowed us to define six T3SS2 subgroups (I-VI) with different repertoires of effector proteins, redefining the concepts of T3SS2 core and accessory effector proteins. Finally, we identified a subset of the T3SS2 gene clusters (subgroup VI) that lacks most T3SS2 effector proteins described to date and provided a list of 10 novel effector candidates for this subgroup through bioinformatic analysis. Collectively, our findings indicate that the T3SS2 extends beyond the family Vibrionaceae and suggest that different effector protein repertories could have a differential impact on the pathogenic potential and environmental fitness of each bacterium that has acquired the Vibrio T3SS2 gene cluster.

Analysis of the Zonula occludens Toxin Found in the Genome of the Chilean Non-toxigenic Vibrio parahaemolyticus Strain PMC53.7.

Vibrio parahaemolyticus non-toxigenic strains are responsible for about 10% of acute gastroenteritis associated with this species, suggesting they harbor unique virulence factors. Zonula occludens toxin (Zot), firstly described in Vibrio cholerae, is a secreted toxin that increases intestinal permeability. Recently, we identified Zot-encoding genes in the genomes of highly cytotoxic Chilean V. parahaemolyticus strains, including the non-toxigenic clinical strain PMC53.7. To gain insights into a possible role of Zot in V. parahaemolyticus, we analyzed whether it could be responsible for cytotoxicity. However, we observed a barely positive correlation between Caco-2 cell membrane damage and Zot mRNA expression during PMC53.7 infection and non-cytotoxicity induction in response to purified PMC53.7-Zot. Unusually, we observed a particular actin disturbance on cells infected with PMC53.7. Based on this observation, we decided to compare the sequence of PMC53.7-Zot with Zot of human pathogenic species such as V. cholerae, Campylobacter concisus, Neisseria meningitidis, and other V. parahaemolyticus strains, using computational tools. The PMC53.7-Zot was compared with other toxins and identified as an endotoxin with conserved motifs in the N-terminus and a variable C-terminal region and without FCIGRL peptide. Notably, the C-terminal diversity among Zots meant that not all of them could be identified as toxins. Structurally, PMC53.7-Zot was modeled as a transmembrane protein. Our results suggested that it has partial 3D structure similarity with V. cholerae-Zot. Probably, the PMC53.7-Zot would affect the actin cytoskeletal, but, in the absence of FCIGRL, the mechanisms of actions must be elucidated.

Draft Genome Sequence of Pseudomonas sp. Strain M7D1, Isolated from the Rhizosphere of Desert Bloom Plants.

We announce the draft genome sequence of Pseudomonas sp. strain M7D1, isolated from the rhizosphere of a plant in the Atacama Desert bloom event. The genome sequence had 6,170,633 bp with a G+C content of 59.9%. This draft genome sequence gives information about the presence of genes related to iron acquisition, alleviation of abiotic stress, and other essential traits of plant growth-promoting rhizobacteria.

Time-Frequency Approach Applied to Finding Interaction Regions in Pathogenic Proteins.

Protein-protein interactions govern all molecular processes for living organisms, even those involved in pathogen infection. Pathogens such as virus, bacteria, and parasites contain proteins that help the pathogen to attach, penetrate, and settle inside the target cell. Thus, it is necessary to know the regions in pathogenic proteins that interact with host cell receptors. Currently, powerful pathogen databases are available and many pathogenic proteins have been recognized, but many pathogenic proteins have not been characterized. This work developed a program in MATLAB environment based on the time-frequency analysis to recognize important sites in proteins. Our program highlights the highest energy patches in proteins from their time-frequency distribution and matches the corresponding frequency. We sought to know if this approach is able to recognize stretches residues related to interaction. Our approach was applied to five study cases from pathogenic co-crystallized structures that have been well characterized. We searched the frequencies that characterize interaction regions in pathogenic proteins and with this information tried to identify new interaction patches in either paralogs or orthologs. We found that our program generates a well-interpretable graphic under several descriptors that can show important regions in proteins even those related to interaction. We propose that this MATLAB program could be used as a tool to explore outstanding regions in uncharacterized proteins.

Multiple Ribosomal RNA Operons in Bacteria; Their Concerted Evolution and Potential Consequences on the Rate of Evolution of Their 16S rRNA.

Bacterial species differ greatly in the number and location of the rRNA operons which may be present in the bacterial chromosomes and plasmids. Most bacterial species contain more than one ribosomal RNA operon copy in their genomes, with some species containing up to 15 such copies. We review the number and location of the rRNA operons and discuss evolution of 16S rRNA (rrs) genes -which are considered as ultimate chronometers for phylogenetic classification- in bacteria with multiple copies of these genes. In these bacterial species, the rrs genes must evolve in concert and sequence changes generated by mutation or horizontal gene transfer must be either erased or spread to every gene copy to avoid divergence, as it occurs when they are present in different species. Analysis of polymorphic sites in intra-genomic rrs copies identifies putative conversion events and demonstrates that sequence conversion is patchy and occurs in small conversion tracts. Sequence conversion probably arises by a non-reciprocal transfer between two or more copies where one copy contributes only a small contiguous segment of DNA, whereas the other copy contributes the rest of the genome in a fairly well understood molecular process. Because concerted evolution implies that a mutation in any of the rrs copies is either eliminated or transferred to every rrs gene in the genome, this process should slow their evolution rate relative to that of single copy genes. However, available data on the rrs genes in bacterial genomes do not show a clear relationship between their evolution rates and the number of their copies in the genome.

Exploring the Genomic Traits of Non-toxigenic Vibrio parahaemolyticus Strains Isolated in Southern Chile.

Vibrio parahaemolyticus is the leading cause of seafood-borne gastroenteritis worldwide. As reported in other countries, after the rise and fall of the pandemic strain in Chile, other post-pandemic strains have been associated with clinical cases, including strains lacking the major toxins TDH and TRH. Since the presence or absence of tdh and trh genes has been used for diagnostic purposes and as a proxy of the virulence of V. parahaemolyticus isolates, the understanding of virulence in V. parahaemolyticus strains lacking toxins is essential to detect these strains present in water and marine products to avoid possible food-borne infection. In this study, we characterized the genome of four environmental and two clinical non-toxigenic strains (tdh-, trh-, and T3SS2-). Using whole-genome sequencing, phylogenetic, and comparative genome analysis, we identified the core and pan-genome of V. parahaemolyticus of strains of southern Chile. The phylogenetic tree based on the core genome showed low genetic diversity but the analysis of the pan-genome revealed that all strains harbored genomic islands carrying diverse virulence and fitness factors or prophage-like elements that encode toxins like Zot and RTX. Interestingly, the three strains carrying Zot-like toxin have a different sequence, although the alignment showed some conserved areas with the zot sequence found in V. cholerae. In addition, we identified an unexpected diversity in the genetic architecture of the T3SS1 gene cluster and the presence of the T3SS2 gene cluster in a non-pandemic environmental strain. Our study sheds light on the diversity of V. parahaemolyticus strains from the southern Pacific which increases our current knowledge regarding the global diversity of this organism.

Insight Into the Origin and Evolution of the Vibrio parahaemolyticus Pandemic Strain.

A strain of Vibrio parahaemolyticus that emerged in 1995 caused the first known pandemic involving this species. This strain comprises clonal autochthonous ocean-dwelling bacteria whose evolution has occurred in the ocean environment. The low sequence diversity in this population enabled the discovery of information on its origin and evolution that has been hidden in bacterial clones that have evolved over a long period. Multilocus sequencing and microarray analysis, together with phylogenetic analysis, of pandemic and pre-pandemic isolates has suggested that the founder clone was an O3:K6 non-pathogenic strain that initially acquired a toxRS/new region and subsequently acquired at least seven novel genomic islands. Sequencing and comparison of whole genomes later confirmed these early observations, and it confirmed that most of the genetic changes occurred via gene conversion involving horizontally transmitted DNA. The highly clonal population rapidly diversified, especially in terms of antigenicity, and 27 serotypes have already been reported. Comparisons of the core genomes derived from the founder clone indicate that there are only a few hundred single-nucleotide variations between isolates. However, when the whole genome is considered (the core plus non-core genome and from any clonal frame), the amount of DNA with a different clonal frame can reach up to 4.2% and the number of single-nucleotide variations can reach several hundred thousand. Altogether, these and previous observations based on multilocus sequence typing, microarray analysis, and whole-genome sequencing indicate the large contribution made by DNA with different clonal genealogy to genome diversification. The evidence also indicates that horizontal gene transfer (HGT) caused the emergence of new pathogens. Furthermore, the extent of HGT seems to depend on the vicissitudes of the life of each bacterium, as exemplified by differences in thousands of base pairs acquired by HGT among almost identical genetic isolates.

Role of Non-coding Regulatory RNA in the Virulence of Human Pathogenic Vibrios.

In recent decades, the identification of small non-coding RNAs in bacteria has revealed an important regulatory mechanism of gene expression involved in the response to environmental signals and to the control of virulence. In the family Vibrionaceae, which includes several human and animal pathogens, small non-coding RNAs (sRNAs) are closely related to important processes including metabolism, quorum sensing, virulence, and fitness. Studies conducted in silico and experiments using microarrays and high-throughput RNA sequencing have led to the discovery of an unexpected number of sRNAs in Vibrios. The present review discusses the most relevant reports regarding the mechanisms of action of sRNAs and their implications in the virulence of the main human pathogens in the family Vibrionaceae: Vibrio parahaemolyticus, V. vulnificus and V. cholerae.

Bacteriophages in the control of pathogenic vibrios

Vibrios are common inhabitants of marine and estuarine environments. Some of them can be pathogenic to humans and/or marine animals using a broad repertory of virulence factors. Lately, several reports have indicated that the incidence of Vibrio infections in humans is rising and also in animals constitute a continuing threat for aquaculture. Moreover, the continuous use of antibiotics has been accompanied by an emergence of antibiotic resistance in Vibrio species, implying a necessity for efficient treatments. One promising alternative that emerges is the use of lytic bacteriophages; however, there are some drawbacks that should be overcome to make phage therapy a widely accepted method. In this work, we discuss about the major pathogenic Vibrio species and the progress, benefits and disadvantages that have been detected during the experimental use of bacteriophages to their control.