Evolutionary perspective on the origin of Haitian cholera outbreak strain

Cholera epidemic has not been reported in Haiti for at least 100 years, although cholera has been present in Latin America since 1991. Surprisingly, the recent cholera epidemic in Haiti (October 2010) recorded more than 250,000 cases and 4000 deaths in the first 6 months and became one of the most explosive and deadly cholera outbreak in recent history. In the present study, we conducted genomic analyses of pathogenicity islands of three Haitian Vibrio cholerae strains and compared them with nine different V. cholerae O1 El Tor genomes. Although CIRS101 is evolutionarily most similar to the Haitian strains, our study also provides some important differences in the genetic organization of pathogenicity islands of Haitian strains with CIRS101. Evolutionary analysis suggests that unusual functional constraints have been imposed on the Haitian strains and we hypothesize that amino acid substitution is more deleterious in Haitian strains than in nonHaitian strains.     

High-Frequency Rugose Exopolysaccharide Production by Vibrio cholerae Strains Isolated in Haiti

In October, 2010, epidemic cholera was reported for the first time in Haiti in over 100 years. Establishment of cholera endemicity in Haiti will be dependent in large part on the continued presence of toxigenic V. cholerae O1 in aquatic reservoirs. The rugose phenotype of V. cholerae, characterized by exopolysaccharide production that confers resistance to environmental stress, is a potential contributor to environmental persistence. Using a microbiologic medium promoting high-frequency conversion of smooth to rugose (S–R) phenotype, 80 (46.5%) of 172 V. cholerae strains isolated from clinical and environmental sources in Haiti were able to convert to a rugose phenotype. Toxigenic V. cholerae O1 strains isolated at the beginning of the epidemic (2010) were significantly less likely to shift to a rugose phenotype than clinical strains isolated in 2012/2013, or environmental strains. Frequency of rugose conversion was influenced by incubation temperature and time. Appearance of the biofilm produced by a Haitian clinical rugose strain (altered biotype El Tor HC16R) differed from that of a typical El Tor rugose strain (N16961R) by confocal microscopy. On whole-genome SNP analysis, there was no phylogenetic clustering of strains showing an ability to shift to a rugose phenotype. Our data confirm the ability of Haitian clinical (and environmental) strains to shift to a protective rugose phenotype, and suggest that factors such as temperature influence the frequency of transition to this phenotype.     

Characterization of VPI pathogenicity island and CTXphi prophage in environmental strains of Vibrio cholerae

Environmental isolates of Vibrio cholerae of eight randomly amplified polymorphic DNA (RAPD) fingerprint types from Calcutta, India, that were unusual in containing toxin-coregulated pilus or cholera toxin genes but not O1 or O139 antigens of epidemic strains were studied by PCR and sequencing to gain insights into V. cholerae evolution. We found that each isolate contained a variant form of the VPI pathogenicity island. Distinguishing features included (i) four new alleles of tcpF (which encodes secreted virulence protein; its exact function is unknown), 20 to 70% divergent (at the protein level) from each other and canonical tcpF; (ii) a new allele of toxT (virulence regulatory gene), 36% divergent (at the protein level) in its 5' half and nearly identical in its 3' half to canonical toxT; (iii) a new tcpA (pilin) gene; and (iv) four variant forms of a regulatory sequence upstream of toxT. Also found were transpositions of an IS903-related element and function-unknown genes to sites in VPI. Cholera toxin (ctx) genes were found in isolates of two RAPD types, in each case embedded in CTXphi-like prophages. Fragments that are inferred to contain only putative repressor, replication, and integration genes were present in two other RAPD types. New possible prophage repressor and replication genes were also identified. Our results show marked genetic diversity in the virulence-associated gene clusters found in some nonepidemic V. cholerae strains, suggest that some of these genes contribute to fitness in nature, and emphasize the potential importance of interstrain gene exchange in the evolution of this species.     

Recombination shapes the structure of an environmental Vibrio cholerae population

Vibrio cholerae consists of pathogenic strains that cause sporadic gastrointestinal illness or epidemic cholera disease and nonpathogenic strains that grow and persist in coastal aquatic ecosystems. Previous studies of disease-causing strains have shown V. cholerae to be a primarily clonal bacterial species, but isolates analyzed have been strongly biased toward pathogenic genotypes, while representing only a small sample of the vast diversity in environmental strains. In this study, we characterized homologous recombination and structure among 152 environmental V. cholerae isolates and 13 other putative Vibrio isolates from coastal waters and sediments in central California, as well as four clinical V. cholerae isolates, using multilocus sequence analysis of seven housekeeping genes. Recombinant regions were identified by at least three detection methods in 72% of our V. cholerae isolates. Despite frequent recombination, significant linkage disequilibrium was still detected among the V. cholerae sequence types. Incongruent but nonrandom associations were observed for maximum likelihood topologies from the individual loci. Overall, our estimated recombination rate in V. cholerae of 6.5 times the mutation rate is similar to those of other sexual bacteria and appears frequently enough to restrict selection from purging much of the neutral intraspecies diversity. These data suggest that frequent recombination among V. cholerae may hinder the identification of ecotypes in this bacterioplankton population.     

Identifying Likely Transmission Pathways within a 10-Year Community Outbreak of Tuberculosis by High-Depth Whole Genome Sequencing

Background

Improved tuberculosis control and the need to contain the spread of drug-resistant strains provide a strong rationale for exploring tuberculosis transmission dynamics at the population level. Whole-genome sequencing provides optimal strain resolution, facilitating detailed mapping of potential transmission pathways.

Methods

We sequenced 22 isolates from a Mycobacterium tuberculosis cluster in New South Wales, Australia, identified during routine 24-locus mycobacterial interspersed repetitive unit typing. Following high-depth paired-end sequencing using the Illumina HiSeq 2000 platform, two independent pipelines were employed for analysis, both employing read mapping onto reference genomes as well as de novo assembly, to control biases in variant detection. In addition to single-nucleotide polymorphisms, the analyses also sought to identify insertions, deletions and structural variants.

Results

Isolates were highly similar, with a distance of 13 variants between the most distant members of the cluster. The most sensitive analysis classified the 22 isolates into 18 groups. Four of the isolates did not appear to share a recent common ancestor with the largest clade; another four isolates had an uncertain ancestral relationship with the largest clade.

Conclusion

Whole genome sequencing, with analysis of single-nucleotide polymorphisms, insertions, deletions, structural variants and subpopulations, enabled the highest possible level of discrimination between cluster members, clarifying likely transmission pathways and exposing the complexity of strain origin. The analysis provides a basis for targeted public health intervention and enhanced classification of future isolates linked to the cluster.     

Molecular typing of epidemic and nonepidemic Vibrio cholerae isolates and differentiation of V. cholerae and V. mimicus isolates by PCR-single-strand conformation polymorphism analysis

AIMS: To examine the utility of polymerase chain reaction (PCR)-single-strand conformation polymorphism (SSCP) analysis to differentiate epidemic and nonepidemic Vibrio cholerae isolates as well as to differentiate V. cholerae and Vibrio mimicus isolates. METHODS AND RESULTS: By both PCR-restriction fragment length polymorphism (RFLP) and PCR-SSCP analysis of groEL-I on chromosome 1 and groEL-II on chromosome 2, V. cholerae isolates gave distinct profiles compared with V. mimicus isolates. In addition, PCR-SSCP analysis of groEL-I and groEL-II could differentiate between V. cholerae epidemic and nonepidemic isolates. Interestingly, the relationships among strains based on groEL-I from chromosome 1 and groEL-II from chromosome 2 were congruent with each other, highlighting the conserved evolutionary history of both chromosomes in this species. CONCLUSIONS: PCR-SSCP is a powerful typing technique, which has the ability to differentiate V. cholerae and V. mimicus isolates. The epidemic V. cholerae O1/O139 serogroup isolates represent a clonal complex distinct from non-O1/non-O139 isolates that can be identified by PCR-SSCP analysis. SIGNIFICANCE AND IMPACT OF THE STUDY: This study highlights the effectiveness of using reliable molecular typing methods and in particular PCR-SSCP, to identify genetic variation among V. cholerae and V. mimicus isolates.     

Comparative sequence analysis of recA gene among Vibrio cholerae isolates from Iran with globally reported sequences

Aims: To study the genetic relatedness between V. cholerae isolates from Iran and other countries based on housekeeping gene recA sequence analysis. Methods and Results: A 995‐bp region of the recA gene from 24 V. cholerae isolates obtained from human and surface water origins in Iran over a 5‐year period was sequenced and compared with the sequence data from the isolates belonging to other places. Cluster analysis of the constructed dendrogram based on recA sequence divergence for our clinical isolates showed one sequence type (ST), whereas environmental isolates revealed eight STs. Interestingly, one of our environmental isolates was intermixed with clinical isolates in the largest cluster containing the epidemic strains. Our 24 isolates plus 198 global isolates available in the GenBank showed 77 sequence types (STs) with at least one nucleotide difference. Conclusions: Our result suggested that recA sequencing is a reliable analysis method for understanding the relatedness of the local isolates with the isolates obtained elsewhere. Significance and Impact of the Study: Understanding the genetic relatedness between V. cholerae isolates could give insights into the health care system for better control and prevention of the cholera.     

Evolutionary perspective on the origin of Haitian cholera outbreak strain

Cholera epidemic has not been reported in Haiti for at least 100?years, although cholera has been present in Latin America since 1991. Surprisingly, the recent cholera epidemic in Haiti (October 2010) recorded more than 250,000 cases and 4000 deaths in the first 6?months and became one of the most explosive and deadly cholera outbreak in recent history. In the present study, we conducted genomic analyses of pathogenicity islands of three Haitian Vibrio cholerae strains and compared them with nine different V. cholerae O1 El Tor genomes. Although CIRS101 is evolutionarily most similar to the Haitian strains, our study also provides some important differences in the genetic organization of pathogenicity islands of Haitian strains with CIRS101. Evolutionary analysis suggests that unusual functional constraints have been imposed on the Haitian strains and we hypothesize that amino acid substitution is more deleterious in Haitian strains than in nonHaitian strains.     

Molecular analyses of a putative CTXphi precursor and evidence for independent acquisition of distinct CTX(phi)s by toxigenic Vibrio cholerae

The genes encoding cholera toxin (ctxA and ctxB) are encoded in the genome of CTXphi, a filamentous phage that infects Vibrio cholerae. To study the evolutionary history of CTXphi, we examined genome diversity in CTX(phi)s derived from a variety of epidemic and nonepidemic Vibrio sp. natural isolates. Among these were three V. cholerae strains that contained CTX prophage sequences but not the ctxA and ctxB genes. These prophages each gave rise to a plasmid form whose genomic organization was very similar to that of the CTXphi replicative form, with the exception of missing ctxAB. Sequence analysis of these three plasmids revealed that they lacked the upstream control region normally found 5' of ctxA, as well as the ctxAB promoter region and coding sequences. These findings are consistent with the hypothesis that a CTXphi precursor that lacked ctxAB simultaneously acquired the toxin genes and their regulatory sequences. To assess the evolutionary relationships among additional CTX(phi)s, two CTXphi-encoded genes, orfU and zot, were sequenced from 13 V. cholerae and 4 V. mimicus isolates. Comparative nucleotide sequence analyses revealed that the CTX(phi)s derived from classical and El Tor V. cholerae isolates comprise two distinct lineages within otherwise nearly identical chromosomal backgrounds (based on mdh sequences). These findings suggest that nontoxigenic precursors of the two V. cholerae O1 biotypes independently acquired distinct CTX(phi)s.     

966. CLAVIJA DOMINGENSIS

Clavija domingensis Urb. & Ekman was one of the many Haitian endemics that were described based on collections made by the great Swedish botanist Leonard Ekman between 1924 and 1928. The species is Critically Endangered sensu IUCN (criteria c2a(i); D) and it is currently the focus of conservation initiatives in Jardin Botanique des Cayes (Haiti), Jardín Botánico Nacional Dr. Rafael M. Moscoso (Dominican Republic), and Fairchild Tropical Botanic Garden (U.S.A.). Now known from only six localities from southern Haiti, each locality only represents a single individual. The species is illustrated based on plants grown in Fairchild Tropical Botanic Garden.     

Toward genomic prediction from whole-genome sequence data: impact of sequencing design on genotype imputation and accuracy of predictions

Genomic prediction from whole-genome sequence data is attractive, as the accuracy of genomic prediction is no longer bounded by extent of linkage disequilibrium between DNA markers and causal mutations affecting the trait, given the causal mutations are in the data set. A cost-effective strategy could be to sequence a small proportion of the population, and impute sequence data to the rest of the reference population. Here, we describe strategies for selecting individuals for sequencing, based on either pedigree relationships or haplotype diversity. Performance of these strategies (number of variants detected and accuracy of imputation) were evaluated in sequence data simulated through a real Belgian Blue cattle pedigree. A strategy (AHAP), which selected a subset of individuals for sequencing that maximized the number of unique haplotypes (from single-nucleotide polymorphism panel data) sequenced gave good performance across a range of variant minor allele frequencies. We then investigated the optimum number of individuals to sequence by fold coverage given a maximum total sequencing effort. At 600 total fold coverage (x 600), the optimum strategy was to sequence 75 individuals at eightfold coverage. Finally, we investigated the accuracy of genomic predictions that could be achieved. The advantage of using imputed sequence data compared with dense SNP array genotypes was highly dependent on the allele frequency spectrum of the causative mutations affecting the trait. When this followed a neutral distribution, the advantage of the imputed sequence data was small; however, when the causal mutations all had low minor allele frequencies, using the sequence data improved the accuracy of genomic prediction by up to 30%.     

Genome sequencing reveals unique mutations in characteristic metabolic pathways and the transfer of virulence genes between V. mimicus and V. cholerae

Vibrio mimicus, the species most similar to V. cholerae, is a microbe present in the natural environmental and sometimes causes diarrhea and internal infections in humans. It shows similar phenotypes to V. cholerae but differs in some biochemical characteristics. The molecular mechanisms underlying the differences in biochemical metabolism between V. mimicus and V. cholerae are currently unclear. Several V. mimicus isolates have been found that carry cholera toxin genes (ctxAB) and cause cholera-like diarrhea in humans. Here, the genome of the V. mimicus isolate SX-4, which carries an intact CTX element, was sequenced and annotated. Analysis of its genome, together with those of other Vibrio species, revealed extensive differences within the Vibrionaceae. Common mutations in gene clusters involved in three biochemical metabolism pathways that are used for discrimination between V. mimicus and V. cholerae were found in V. mimicus strains. We also constructed detailed genomic structures and evolution maps for the general types of genomic drift associated with pathogenic characters in polysaccharides, CTX elements and toxin co-regulated pilus (TCP) gene clusters. Overall, the whole-genome sequencing of the V. mimicus strain carrying the cholera toxin gene provides detailed information for understanding genomic differences among Vibrio spp. V. mimicus has a large number of diverse gene and nucleotide differences from its nearest neighbor, V. cholerae. The observed mutations in the characteristic metabolism pathways may indicate different adaptations to different niches for these species and may be caused by ancient events in evolution before the divergence of V. cholerae and V. mimicus. Horizontal transfers of virulence-related genes from an uncommon clone of V. cholerae, rather than the seventh pandemic strains, have generated the pathogenic V. mimicus strain carrying cholera toxin genes.     

Survey for potential soft tick (Acari: Argasidae) vectors of African swine fever on the Island of Hispaniola

A survey of the occurrence ofOrnithodoros ticks in animal burrows was conducted in Haiti and the Dominican Republic. Vacuum sampling techniques were used. Four of 8 sample sites in Haiti and 10 of 70 sample sites in the Dominican Republic were positive forO. puertoricensis Fox. Positive sample sites in Haiti were usually near swine. Sites in the Dominican Republic were in drier regions of the country and were not directly associated with previous swine locations. The determination thatO. puertoricensis, a potential vector and reservoir of African swine fever (ASF), is present in this region may pose a serious problem for eradication of ASF from the island of Hispaniola.Florida Agricultural Experiment Station Journal Series no. 4585.     

A study of the prevalence of regulatory genes controlling virulence gene expression among Vibrio choleraeeltor biovariant strains varying in their pandemic potential

The evolution of the genome of the pathogenic agent of the seventh cholera pandemia Vibrio cholerae eltor biovariant was thought to occur by acquiring not only structural genes of virulence but also regulatory systems as a result of horizontal transfer events. The polymerase chain reaction revealed the presence of the following regulatory genes that control the virulence gene expression in the chromosome of pre-pandemic and pandemic strains of cholera vibrios eltor: toxR, toxT, tcpP, tcpH, luxS, luxO, crp, vicH, pepA. The avirulent V. cholerae strain ATCC14033 isolated in 1910 (hypothetical predecessor of the cholera eltor agent) was shown to be lacking the regulatory genes toxT, tcpP, tcpHlocalized in the pathogenicity island VPI-1, and to be capable of realizing positive control over the expression of the virulence genes involved in the ToxR regulon. The virulent strains isolated from cholera patients during the local cholera outbreak in Indonesia in 1937 did not differ from the strains that caused cholera eltor pandemic in 1961. The strains had identical content of the regulatory genes tested. Only one strain of the four isolates studied contained no tcpPgene. Two key regulatory genes, toxR and toxT, were sequenced in all the isolates. The toxR nucleotide sequence of three pre-pandemic strains was shown to be indistinguishable from that of the pandemic isolates. On the other hand, the clinical strain MAK757 isolated prior to the emergence of the epidemic demonstrated an altered nucleotide sequence in its toxR gene. Experiments with the intra-intestinal challenge of suckling rabbits were indicative of similar virulence levels for the pre-pandemic and pandemic clinical strains. These results may serve as the evidence of the in vivo activity of the pre-pandemic strains of the toxT, tcpH, and tcpP positive regulatory genes that acquired in V. cholerae during the evolutionary process.     

Epidemiological study of Vibrio cholerae using variable number of tandem repeats

By conventional genetic methods, including pulse-field gel electrophoresis and multilocus sequence typing, most pathogenic, cholera toxin-positive O1 and O139 isolates of Vibrio cholerae cannot be distinguished. We evaluated relationships among 173 V. cholerae isolates collected between 1992 and 2007 from different geographic areas in India by analyzing five variable number of tandem repeat (VNTR) loci. Each VNTR locus was highly variable, with between 5 and 19 alleles. eburst analysis revealed four large groups of genetically related isolates. Two groups contained genotypes of isolates with the O139 serogroup (which emerged for the first time in epidemic form in 1992), with the other two groups containing O1 strains. In subsequent analysis, it was possible to track the spread of specific genotypes across time and space. Our data highlight the utility of the methodology as an epidemiologic tool for assessing spread of isolates in both epidemic and endemic settings.     

Enterobacterial repetitive intergenic consensus sequences and the PCR to generate fingerprints of genomic DNAs from Vibrio cholerae O1, O139, and non-O1 strains.

Enterobacterial repetitive intergenic consensus (ERIC) sequence polymorphism was studied in Vibrio Cholerae strains isolated before and after the cholera epidemic in Brazil (in 1991), along with epidemic strains from Peru, Mexico, and India, by PCR. A total of 17 fingerprint patterns (FPs) were detected in the V. cholerae strains examined; 96.7% of the toxigenic V. cholerae O1 strains and 100% of the O139 serogroup strains were found to belong to the same FP group comprising four fragments (FP1). The nontoxigenic V. cholerae O1 also yielded four fragments but constituted a different FP group (FP2). A total of 15 different patterns were observed among the V. cholerae non-O1 strains. Two patterns were observed most frequently for V. cholerae non-01 strains, 25% of which have FP3, with five fragments, and 16.7% of which have FP4, with two fragments. Three fragments, 1.75, 0.79, and 0.5 kb, were found to be common to both toxigenic and nontoxigenic V. cholerae O1 strains as well as to group FP3, containing V. cholerae non-O1 strains. Two fragments of group FP3, 1.3 and 1.0 kb, were present in FP1 and FP2 respectively. The 0.5-kb fragment was common to all strains and serogroups of V. cholerae analyzed. It is concluded from the results of this study, based on DNA FPs of environmental isolates, that it is possible to detect an emerging virulent strain in a cholera-endemic region. ERIC-PCR constitutes a powerful tool for determination of the virulence potential of V. cholerae O1 strains isolated in surveillance programs and for molecular epidemiological investigations.     

霍乱弧菌和副溶血弧菌分离株的gyrB基因系统发育分析

依据gyrB基因部分编码序列构建系统发育树以分类和鉴别霍乱弧菌和副溶血弧菌,并探讨其种系发生关系。扩增并测序13株霍乱弧菌、8株副溶血弧菌、2株嗜水气单胞菌及1株类志贺邻单胞菌的gyrB基因(编码DNA促旋酶B亚单位)序列,并采用距离法与最大似然法构建系统发育树。两种方法所构建的树结构完全一致,霍乱弧菌、副溶血弧菌、嗜水气单胞菌及类志贺邻单胞菌各自形成一个独立的簇。其中,霍乱肠毒素基因(ctxA)阳性的霍乱弧菌(8株O139群与2株O1群ElTor型)聚类成一分枝;3株副溶血弧菌临床株(1株2002年流行株,2株2004年分离株)与1日本菌株及2001年1株自环境分离的毒力株聚类。系统发育分析靶分子gyrB基因可以良好区分上述4种常见病原菌。产毒O139群霍乱弧菌与产毒O1群ElTor型霍乱弧菌关系密切。副溶血弧菌环境毒力株与本地区临床主要流行株在系统发育关系上较为接近,可能是潜在的致病菌。     

A physical map of the highly heterozygous Populus genome: integration with the genome sequence and genetic map and analysis of haplotype variation

As part of a larger project to sequence the Populus genome and generate genomic resources for this emerging model tree, we constructed a physical map of the Populus genome, representing one of the few such maps of an undomesticated, highly heterozygous plant species. The physical map, consisting of 2802 contigs, was constructed from fingerprinted bacterial artificial chromosome (BAC) clones. The map represents approximately 9.4-fold coverage of the Populus genome, which has been estimated from the genome sequence assembly to be 485 ± 10 Mb in size. BAC ends were sequenced to assist long-range assembly of whole-genome shotgun sequence scaffolds and to anchor the physical map to the genome sequence. Simple sequence repeat-based markers were derived from the end sequences and used to initiate integration of the BAC and genetic maps. A total of 2411 physical map contigs, representing 97% of all clones assigned to contigs, were aligned to the sequence assembly (JGI Populus trichocarpa , version 1.0). These alignments represent a total coverage of 384 Mb (79%) of the entire poplar sequence assembly and 295 Mb (96%) of linkage group sequence assemblies. A striking result of the physical map contig alignments to the sequence assembly was the co-localization of multiple contigs across numerous regions of the 19 linkage groups. Targeted sequencing of BAC clones and genetic analysis in a small number of representative regions showed that these co-aligning contigs represent distinct haplotypes in the heterozygous individual sequenced, and revealed the nature of these haplotype sequence differences.     

The Population Structure of Vibrio cholerae from the Chandigarh Region of Northern India

Background

Cholera infection continues to be a threat to global public health. The current cholera pandemic associated with Vibrio cholerae El Tor has now been ongoing for over half a century.

Methodology/Principal Findings

Thirty-eight V. cholerae El Tor isolates associated with a cholera outbreak in 2009 from the Chandigarh region of India were characterised by a combination of microbiology, molecular typing and whole-genome sequencing. The genomic analysis indicated that two clones of V. cholera circulated in the region and caused disease during this time. These clones fell into two distinct sub-clades that map independently onto wave 3 of the phylogenetic tree of seventh pandemic V. cholerae El Tor. Sequence analyses of the cholera toxin gene, the Vibrio seventh Pandemic Island II (VSPII) and SXT element correlated with this phylogenetic position of the two clades on the El Tor tree. The clade 2 isolates, characterized by a drug-resistant profile and the expression of a distinct cholera toxin, are closely related to the recent V. cholerae isolated elsewhere, including Haiti, but fell on a distinct branch of the tree, showing they were independent outbreaks. Multi-Locus Sequence Typing (MLST) distinguishes two sequence types among the 38 isolates, that did not correspond to the clades defined by whole-genome sequencing. Multi-Locus Variable-length tandem-nucleotide repeat Analysis (MLVA) identified 16 distinct clusters.

Conclusions/Significance

The use of whole-genome sequencing enabled the identification of two clones of V. cholerae that circulated during the 2009 Chandigarh outbreak. These clones harboured a similar structure of ICEVchHai1 but differed mainly in the structure of CTX phage and VSPII. The limited capacity of MLST and MLVA to discriminate between the clones that circulated in the 2009 Chandigarh outbreak highlights the value of whole-genome sequencing as a route to the identification of further genetic markers to subtype V. cholerae isolates.