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.

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.

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.

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.

Inactivated E. coli transformed with plasmids that produce dsRNA against infectious salmon anemia virus hemagglutinin show antiviral activity when added to infected ASK cells.

Infectious salmon anemia virus (ISAV) has caused great losses to the Chilean salmon industry, and the success of prevention and treatment strategies is uncertain. The use of RNA interference (RNAi) is a promising approach because during the replication cycle, the ISAV genome must be transcribed to mRNA in the cytoplasm. We explored the capacity of E. coli transformed with plasmids that produce double-stranded RNA (dsRNA) to induce antiviral activity when added to infected ASK cells. We transformed the non-pathogenic Escherichia coli HT115 (DE3) with plasmids that expressed highly conserved regions of the ISAV genes encoding the nucleoprotein (NP), fusion (F), hemagglutinin (HE), and matrix (M) proteins as dsRNA, which is the precursor of the RNAi mechanism. The inactivated transformed bacteria carrying dsRNA were tested for their capacity to silence the target ISAV genes, and the dsRNA that were able to inhibit gene expression were subsequently tested for their ability to attenuate the cytopathic effect (CPE) and reduce the viral load. Of the four target genes tested, inactivated E. coli transformed with plasmids producing dsRNA targeting HE showed antiviral activity when added to infected ASK cells.

Genome diversification within a clonal population of pandemic Vibrio parahaemolyticus seems to depend on the life circumstances of each individual bacteria.

BACKGROUND: New strains of Vibrio parahaemolyticus that cause diarrhea in humans by seafood ingestion periodically emerge through continuous evolution in the ocean. Influx and expansion in the Southern Chilean ocean of a highly clonal V. parahaemolyticus (serotype O3:K6) population from South East Asia caused one of the largest seafood-related diarrhea outbreaks in the world. Here, genomics analyses of isolates from this rapidly expanding clonal population offered an opportunity to observe the molecular evolutionary changes often obscured in more diverse populations. RESULTS: Whole genome sequence comparison of eight independent isolates of this population from mussels or clinical cases (from different years) was performed. Differences of 1366 to 217,729 bp genome length and 13 to 164 bp single nucleotide variants (SNVs) were found. Most genomic differences corresponded to the presence of regions unique to only one or two isolates, and were probably acquired by horizontal gene transfer (HGT). Some DNA gain was chromosomal but most was in plasmids. One isolate had a large region (8,644 bp) missing, which was probably caused by excision of a prophage. Genome innovation by the presence of unique DNA, attributable to HGT from related bacteria, varied greatly among the isolates, with values of 1,366 (ten times the number of highest number of SNVs) to 217,729 (a thousand times more than the number of highest number of SNVs). CONCLUSIONS: The evolutionary forces (SNVs, HGT) acting on each isolate of the same population were found to differ to an extent that probably depended on the ecological scenario and life circumstances of each bacterium.

Complete Genome Sequence of Vibrio anguillarum Phage CHOED Successfully Used for Phage Therapy in Aquaculture.

Vibrio anguillarum phage CHOED was isolated from Chilean mussels. It is a virulent phage showing effective inhibition of V. anguillarum. CHOED has potential in phage therapy, because it can protect fish from vibriosis in fish farms. Here, we announce the completely sequenced genome of V. anguillarum phage CHOED.

Rise and fall of pandemic Vibrio parahaemolyticus serotype O3:K6 in southern Chile.

Seafood consumption-related diarrhoea increased drastically in Chile when the pandemic strain of Vibrio parahaemolyticus serotype O3:K6 reached Region de Los Lagos, where most of Chile's seafood is produced. Outbreaks peaked in 2005 with 3725 clinical cases in this region and gradually decreased to fewer than 10 cases in 2010 and 2011. We show here that the pandemic strain concurrently vanished from mussels; we also report further environmental data. Integration of the 2010/2011 data with those obtained since 2004 suggests that after its arrival in southern Chile, the pandemic strain grew in mussels, likely facilitated by a minor rise in surface seawater temperature and by warming of the mussels in the intertidal region due to frequent sunny days. However, since these environmental parameters probably equally affected the pandemic strain and more than 30 V. parahaemolyticus DNA restriction clusters that inhabit local shellfish, a selective effect of bacteriophages is proposed. Lytic bacteriophage VP93 may have favoured the growth of the pandemic strain versus similar phage-sensitive strains, as shown here in a particular case. However, the pandemic strain's decline may have been promoted by temperate phage VP58.5, which kills the pandemic strain and increases the UV sensitivity of lysogenized phage-resistant cells.

Microevolution of pandemic Vibrio parahaemolyticus assessed by the number of repeat units in short sequence tandem repeat regions.

The emergence of the pandemic strain Vibrio parahaemolyticus O3:K6 in 1996 caused a large increase of diarrhea outbreaks related to seafood consumption in Southeast Asia, and later worldwide. Isolates of this strain constitutes a clonal complex, and their effectual differentiation is possible by comparison of their variable number tandem repeats (VNTRs). The differentiation of the isolates by the differences in VNTRs will allow inferring the population dynamics and microevolution of this strain but this requires knowing the rate and mechanism of VNTRs' variation. Our study of mutants obtained after serial cultivation of clones showed that mutation rates of the six VNTRs examined are on the order of 10(-4) mutant per generation and that difference increases by stepwise addition of single mutations. The single stepwise mutation (SSM) was deduced because mutants with 1, 2, 3, or more repeat unit deletions or insertions follow a geometric distribution. Plausible phylogenetic trees are obtained when, according to SSM, the genetic distance between clusters with different number of repeats is assessed by the absolute differences in repeats. Using this approach, mutants originated from different isolates of pandemic V. parahaemolyticus after serial cultivation are clustered with their parental isolates. Additionally, isolates of pandemic V. parahaemolyticus from Southeast Asia, Tokyo, and northern and southern Chile are clustered according their geographical origin. The deepest split in these four populations is observed between the Tokyo and southern Chile populations. We conclude that proper phylogenetic relations and successful tracing of pandemic V. parahaemolyticus requires measuring the differences between isolates by the absolute number of repeats in the VNTRs considered.

Origins and colonization history of pandemic Vibrio parahaemolyticus in South America.

The dynamics of dissemination of the environmental human pathogen Vibrio parahaemolyticus are uncertain. The O3:K6 clone was restricted to Asia until its detection along the Peruvian coasts and in northern Chile in 1997 in phase with the arrival of El Niño waters. A subsequent emergence of O3:K6 strains was detected in austral Chile in 2004. The origin of these 1997 and 2004 population radiations has not yet been conclusively determined. Multiple loci VNTR analysis using seven polymorphic loci was carried out with a number of representative strains from Asia, Peru and Chile to determine their genetic characteristics and population structure. Asian and Chilean subpopulations were the most genetically distant groups with an intermediate subpopulation in Peru. Population structure inferred from a minimum-spanning tree and Bayesian analysis divided the populations into two genetically distinct groups, consistent with the epidemic dynamics of the O3:K6 clone in South America. One group comprised strains from the original Asiatic population and strains arriving in Peru and Chile in 1997. The second group included the remaining Peruvian Strains and Chilean strains obtained from Puerto Montt in 2004. The analysis of the arrival of the O3:K6 clone at the Pacific coasts of South America has provided novel insights linking the origin of the invasion in 1997 to Asian populations and describing the successful establishment of the O3:K6 populations, first in Peru and subsequently in the South of Chile owing to a possible radiation of Peruvian populations.

A new group of cosmopolitan bacteriophages induce a carrier state in the pandemic strain of Vibrio parahaemolyticus.

A clonal population of pathogenic Vibrio parahaemolyticus O3 : K6 serovar has spread in coastal waters, causing outbreaks worldwide since 1996. Bacteriophage infection is one of the main factors affecting bacterial strain concentration in the ocean. We studied the occurrence and properties of phages infecting this V. parahaemolyticus pandemic strain in coastal waters. Analysing 143 samples, phages were found in 13. All isolates clustered in a closely related group of podophages with at least 90% nucleotide sequence identity in three essential genes, despite distant geographical origins. These bacteriophages were able to multiply on the V. parahaemolyticus pandemic strain, but the impact on host concentration and subsequent growth was negligible. Infected bacteria continued producing the phage but were not lysogenized. The phage genome of prototype strain VP93 is 43 931 nucleotides and contains 337 bp direct terminal repeats at both ends. VP93 is the first non-Pseudomonas phage related to the PhiKMV-like subgroup of the T7 supergroup. The lack of a major effect on host growth suggests that these phages exert little control on the propagation of the pandemic strain in the environment. This form of phage growth can be modelled if phage-sensitive and -resistant cells that convert to each other with a high frequency are present in clonal cultures of pandemic V. parahaemolyticus.

Dynamics of clinical and environmental Vibrio parahaemolyticus strains during seafood-related summer diarrhea outbreaks in southern Chile.

Seafood consumption-related diarrhea became prevalent in Chile when the pandemic strain of Vibrio parahaemolyticus serotype O3:K6 reached a region in the south of Chile (Region de los Lagos) where approximately 80% of the country's seafood is produced. In spite of the large outbreaks of clinical infection, the load of V. parahaemolyticus in shellfish of this region is relatively low. The pandemic strain constitutes a small but relatively stable group of a diverse V. parahaemolyticus population, composed of at least 28 genetic groups. Outbreaks in Region de los Lagos began in 2004 and reached a peak in 2005 with 3,725 clinical cases, all associated with the pandemic strain. After 2005, reported cases steadily decreased to a total of 477 cases in 2007. At that time, 40% of the clinical cases were associated with a pandemic strain of a different serotype (O3:K59), and 27% were related to V. parahaemolyticus isolates unrelated to the pandemic strain. In the results published here, we report that in the summer of 2008, when reported cases unexpectedly increased from 477 to 1,143, 98% of the clinical cases were associated with the pandemic strain serotype O3:K6, a change from 2007. Nevertheless, in 2009, when clinical cases decreased to 441, only 64% were related to the pandemic strain; the remaining cases were related to a nonpandemic tdh- and trh-negative strain first identified in shellfish in 2006. Overall, our observations indicate that the pandemic strain has become a relatively stable subpopulation and that when the number of diarrhea cases related to the pandemic strain is low, previously undetected V. parahaemolyticus pathogenic strains become evident.

The linear plasmid prophage Vp58.5 of Vibrio parahaemolyticus is closely related to the integrating phage VHML and constitutes a new incompatibility group of telomere phages.

Vibrio parahaemolyticus O3:K6 pandemic strains recovered in Chile frequently possess a 42-kb plasmid which is the prophage of a myovirus. We studied the prototype phage VP58.5 and show that it does not integrate into the host cell chromosome but replicates as a linear plasmid (Vp58.5) with covalently closed ends (telomeres). The Vp58.5 replicon coexists with other plasmid prophages (N15, PY54, and PhiKO2) in the same cell and thus belongs to a new incompatibility group of telomere phages. We determined the complete nucleotide sequence (42,612 nucleotides) of the VP58.5 phage DNA and compared it with that of the plasmid prophage. The two molecules share the same nucleotide sequence but are 35% circularly permuted to each other. In contrast to the hairpin ends of the plasmid, VP58.5 phage DNA contains 5'-protruding ends. The VP58.5 sequence is 92% identical to the sequence of phage VHML, which was reported to integrate into the host chromosome. However, the gene order and termini of the phage DNAs are different. The VHML genome exhibits the same gene order as does the Vp58.5 plasmid. VHML phage DNA has been reported to contain terminal inverted repeats. This repetitive sequence is similar to the telomere resolution site (telRL) of VP58.5 which, after processing by the phage protelomerase, forms the hairpin ends of the Vp58.5 prophage. It is discussed why these closely related phages may be so different in terms of their genome ends and their lifestyle.

Multiple-locus variable-number tandem-repeat analysis for clonal identification of Vibrio parahaemolyticus isolates by using capillary electrophoresis.

Epidemics of Vibrio parahaemolyticus in Chile have occurred since 1998. Direct genome restriction enzyme analysis (DGREA) using conventional gel electrophoresis permitted discrimination of different V. parahaemolyticus isolates obtained from these outbreaks and showed that this species consists of a highly diverse population. A multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) approach was developed and applied to 22 clinical and 91 environmental V. parahaemolyticus isolates from Chile to understand their clonal structures. To this end, an advanced molecular technique was developed by applying multiplex PCR, fluorescent primers, and capillary electrophoresis, resulting in a high-resolution and high-throughput (HRHT) genotyping method. The genomic basis of this HRHT method was eight VNTR loci described previously by Kimura et al. (J. Microbiol. Methods 72:313-320, 2008) and two new loci which were identified by a detailed molecular study of 24 potential VNTR loci on both chromosomes. The isolates of V. parahaemolyticus belonging to the same DGREA pattern were distinguishable by the size variations in the indicative 10 VNTRs. This assay showed that these 10 VNTR loci were useful for distinguishing isolates of V. parahaemolyticus that had different DGREA patterns and also isolates that belong to the same group. Isolates that differed in their DGREA patterns showed polymorphism in their VNTR profiles. A total of 81 isolates was associated with 59 MLVA groups, providing fine-scale differentiation, even among very closely related isolates. The developed approach enables rapid and high-resolution analysis of V. parahaemolyticus with pandemic potential and provides a new surveillance tool for food-borne pathogens.

Epidemiology of Vibrio parahaemolyticus outbreaks, southern Chile.

Disease outbreaks caused by Vibrio parahaemolyticus in Puerto Montt, Chile, began in 2004 and reached a peak in 2005 at 3,600 clinical cases. Until 2006, every analyzed case was caused by the serovar O3:K6 pandemic strain. In the summer of 2007, only 475 cases were reported; 73% corresponded to the pandemic strain. This decrease was associated with a change in serotype of many pandemic isolates to O3:K59 and the emergence of new clinical strains. One of these strains, associated with 11% of the cases, was genotypically different from the pandemic strain but contained genes that were identical to those found on its pathogenicity island. These findings suggest that pathogenicity-related genes were laterally transferred from the pandemic strain to one of the different V. parahaemolyticus groups comprising the diverse and shifting bacterial population in shellfish in this region.

Enhancement of UV light sensitivity of a Vibrio parahaemolyticus O3:K6 pandemic strain due to natural lysogenization by a telomeric phage.

The Vibrio parahaemolyticus O3:K6 pandemic clonal strain was first observed in southern Chile in 2004 and has since caused approximately 8,000 seafood-related diarrhea cases in this region. The massive proliferation of the original clonal population offers a unique opportunity to study the evolution of a bacterial pathogen in its natural environment by detection and characterization of emerging bacterial variants. Here, we describe a group of pandemic variants characterized by the presence of a 42-kb extrachromosomal DNA that can be recovered by alkaline extraction. Upon treatment with mitomycin C, these variants lyse with production of a myovirus containing DNA of equal size to the plasmid but which cannot be recovered by alkaline extraction. Plasmid and phage DNAs show similar restriction patterns corresponding to enzyme sites in a circular permutation. Sequenced regions showed 81 to 99% nucleotide similarity to bacteriophage VHML of Vibrio harveyi. Altogether these observations indicate that the 42-kb plasmid corresponds to a prophage, consisting of a linear DNA with terminal hairpins of a telomeric temperate phage with a linear genome. Bacteria containing the prophage were 7 to 15 times more sensitive to UV radiation, likely due to phage induction by UV irradiation as plasmid curing restored the original sensitivity. The enhanced UV sensitivity could have a significant role in reducing the survival and propagation capability of the V. parahaemolyticus pandemic strain in the ocean.

Determination of molecular phylogenetics of Vibrio parahaemolyticus strains by multilocus sequence typing.

Vibrio parahaemolyticus is an important human pathogen whose transmission is associated with the consumption of contaminated seafood. There is a growing public health concern due to the emergence of a pandemic strain causing severe outbreaks worldwide. Many questions remain unanswered regarding the evolution and population structure of V. parahaemolyticus. In this work, we describe a multilocus sequence typing (MLST) scheme for V. parahaemolyticus based on the internal fragment sequences of seven housekeeping genes. This MLST scheme was applied to 100 V. parahaemolyticus strains isolated from geographically diverse clinical (n = 37) and environmental (n = 63) sources. The sequences obtained from this work were deposited and are available in a public database (http://pubmlst.org/vparahaemolyticus). Sixty-two unique sequence types were identified, and most (50) were represented by a single isolate, suggesting a high level of genetic diversity. Three major clonal complexes were identified by eBURST analysis. Separate clonal complexes were observed for V. parahaemolyticus isolates originating from the Pacific and Gulf coasts of the United States, while a third clonal complex consisted of strains belonging to the pandemic clonal complex with worldwide distribution. The data reported in this study indicate that V. parahaemolyticus is genetically diverse with a semiclonal population structure and an epidemic structure similar to that of Vibrio cholerae. Genetic diversity in V. parahaemolyticus appears to be driven primarily by frequent recombination rather than mutation, with recombination ratios estimated at 2.5:1 and 8.8:1 by allele and site, respectively. Application of this MLST scheme to more V. parahaemolyticus strains and by different laboratories will facilitate production of a global picture of the epidemiology and evolution of this pathogen.

Intragenomic heterogeneity and intergenomic recombination among Vibrio parahaemolyticus 16S rRNA genes.

Vibrio parahaemolyticus is a marine bacterium bearing 11 copies of ribosomal operons. In some strains, such as RIMD2210633, the genome includes identical copies of 16S rRNA genes (rrs). However, it is known that other strains of the species, such as strains ATCC 17802 and RIMD2210856, show conspicuous intragenomic rrs heterogeneity. The extent and diversity of the rrs heterogeneity in V. parahaemolyticus were studied in further detail by characterization of the rrs copies in environmental isolates belonging to 21 different genotype groups. Thirteen of these groups showed intragenomic heterogeneity, containing altogether 16 sequences differing within a 25 bp segment of their rrs. These sequences grouped into four clusters differing in at least four nucleotide sites. Some isolates contained rrs alleles from up to three different clusters. Each segment sequence conserved the stem-loop characteristic of the 16S rRNA structure of this 25 bp sequence. The double-stranded stem sequence was quite variable, but almost every variation had a compensatory change to maintain seven to eight paired bases. Conversely, the single-strand loop sequence was conserved. The results may be explained as a consequence of recombination among rrs evolving in different bacteria. The results suggest that intergenomic rrs recombination is very high in V. parahaemolyticus and that it occurs solely among Vibrio species. This high rrs homologous intergenomic recombination could be an effective mechanism to maintain intragenomic rrs cohesion, mediating the dispersal of the most abundant rrs version among the 11 intragenomic loci.