Genetic variability of Yersinia pestis isolates as predicted by PCR-based IS100 genotyping and analysis of structural genes encoding glycerol-3-phosphate dehydrogenase (glpD)

Identifying sources of human exposure to plague

Yersinia pestis, the etiologic agent of plague, has shaped the course of human history, killing millions of people in three major pandemics. This bacterium is still endemic in parts of Asia, Africa, and the Americas, where it poses a natural disease threat to human populations. Y. pestis has also recently received attention as a possible bioterrorism agent. Thus, rapid methods to distinguish between bioterrorism and naturally occurring plague infections are of major importance. Our study is the first to demonstrate that variable-number tandem repeats (VNTRs) in the Y. pestis genome can link human case isolates to those obtained from suspected environmental sources of infection. We demonstrate the valuable utility of VNTR markers in epidemiological investigations of naturally occurring plague and the forensic analysis of possible bioterrorism events.

Natural history of plague: perspectives from more than a century of research

For more than a century, scientists have investigated the natural history of plague, a highly fatal disease caused by infection with the gram-negative bacterium Yersinia pestis. Among their most important discoveries were the zoonotic nature of the disease and that plague exists in natural cycles involving transmission between rodent hosts and flea vectors. Other significant findings include those on the evolution of Y. pestis; geographic variation among plague strains; the dynamics and maintenance of transmission cycles; mechanisms by which fleas transmit Y. pestis; resistance and susceptibility among plague hosts; the structure and typology of natural foci; and how landscape features influence the focality, maintenance, and spread of the disease. The knowledge gained from these studies is essential for the development of effective prevention and control strategies.

Microevolution and history of the plague bacillus, Yersinia pestis

The association of historical plague pandemics with Yersinia pestis remains controversial, partly because the evolutionary history of this largely monomorphic bacterium was unknown. The microevolution of Y. pestis was therefore investigated by three different multilocus molecular methods, targeting genomewide synonymous SNPs, variation in number of tandem repeats, and insertion of IS100 insertion elements. Eight populations were recognized by the three methods, and we propose an evolutionary tree for these populations, rooted on Yersinia pseudotuberculosis. The tree invokes microevolution over millennia, during which enzootic pestoides isolates evolved. This initial phase was followed by a binary split 6,500 years ago, which led to populations that are more frequently associated with human disease. These populations do not correspond directly to classical biovars that are based on phenotypic properties. Thus, we recommend that henceforth groupings should be based on molecular signatures. The age of Y. pestis inferred here is compatible with the dates of historical pandemic plague. However, it is premature to infer an association between any modern molecular grouping and a particular pandemic wave that occurred before the 20th century.

Use of highly variable intergenic spacer sequences for multispacer typing of Rickettsia conorii strains

By use of the nearly perfectly colinear genomes of Rickettsia conorii and Rickettsia prowazekii, we compared the usefulness of three types of sequences for typing of R. conorii isolates: (i) 5 variable coding genes comprising the 16S ribosomal DNA, gltA, ompB, and sca4 (gene D) genes, which are present in both genomes, and the ompA gene, which is degraded in R. prowazekii; (ii) 28 genes degraded in R. conorii but intact in R. prowazekii, including 23 split and 5 remnant genes; and (iii) 27 conserved and 25 variable intergenic spacers. The 4 conserved and 23 split genes as well as the 27 conserved intergenic spacers each had identical sequences in 34 human and 5 tick isolates of R. conorii. Analysis of the ompA sequences identified three genotypes of R. conorii. The variable intergenic spacers were significantly more variable than conserved genes, split genes, remnant genes, and conserved spacers (P < 10(-2) in all cases). Four of the variable intergenic spacers (dksA-xerC, mppA-purC, rpmE-tRNA(fMet), and tRNA(Gly)-tRNA(Tyr)) had highly variable sequences; when they were combined for typing, multispacer typing (MST) identified 27 different genotypes in the 39 R. conorii isolates. Two batches from the same R. conorii strain, Malish (Seven), with different culture passage histories were found to exhibit the same MST type. MST was more discriminatory for strain genotyping than multiple gene sequencing (P < 10(-2)). Phylogenetic analysis based on MST sequences was concordant with the geographic origins of R. conorii isolates. Our study supports the usefulness of MST for strain genotyping. This tool may be useful for tracing a strain and identifying its source during outbreaks, including those resulting from bioterrorism.

Insights into the evolution of Yersinia pestis through whole-genome comparison with Yersinia pseudotuberculosis

Yersinia pestis, the causative agent of plague, is a highly uniform clone that diverged recently from the enteric pathogen Yersinia pseudotuberculosis. Despite their close genetic relationship, they differ radically in their pathogenicity and transmission. Here, we report the complete genomic sequence of Y. pseudotuberculosis IP32953 and its use for detailed genome comparisons with available Y. pestis sequences. Analyses of identified differences across a panel of Yersinia isolates from around the world reveal 32 Y. pestis chromosomal genes that, together with the two Y. pestis-specific plasmids, to our knowledge, represent the only new genetic material in Y. pestis acquired since the the divergence from Y. pseudotuberculosis. In contrast, 149 other pseudogenes (doubling the previous estimate) and 317 genes absent from Y. pestis were detected, indicating that as many as 13% of Y. pseudotuberculosis genes no longer function in Y. pestis. Extensive insertion sequence-mediated genome rearrangements and reductive evolution through massive gene loss, resulting in elimination and modification of preexisting gene expression pathways, appear to be more important than acquisition of genes in the evolution of Y. pestis. These results provide a sobering example of how a highly virulent epidemic clone can suddenly emerge from a less virulent, closely related progenitor.

Genotyping, Orientalis-like Yersinia pestis, and plague pandemics

Three pandemics have been attributed to plague in the last 1,500 years. Yersinia pestis caused the third, and its DNA was found in human remains from the second. The Antiqua biovar of Y. pestis may have caused the first pandemic; the other two biovars, Medievalis and Orientalis, may have caused the second and third pandemics, respectively. To test this hypothesis, we designed an original genotyping system based on intergenic spacer sequencing called multiple spacer typing (MST). We found that MST differentiated every biovar in a collection of 36 Y. pestis isolates representative of the three biovars. When MST was applied to dental pulp collected from remains of eight persons who likely died in the first and second pandemics, this system identified original sequences that matched those of Y. pestis Orientalis. These data indicate that Y. pestis caused cases of Justinian plague. The two historical plague pandemics were likely caused by Orientalis-like strains.

Genetics of metabolic variations between Yersinia pestis biovars and the proposal of a new biovar, microtus

Yersinia pestis has been historically divided into three biovars: antiqua, mediaevalis, and orientalis. On the basis of this study, strains from Microtus-related plague foci are proposed to constitute a new biovar, microtus. Based on the ability to ferment glycerol and arabinose and to reduce nitrate, Y. pestis strains can be assigned to one of four biovars: antiqua (glycerol positive, arabinose positive, and nitrate positive), mediaevalis (glycerol positive, arabinose positive, and nitrate negative), orientalis (glycerol negative, arabinose positive, and nitrate positive), and microtus (glycerol positive, arabinose negative, and nitrate negative). A 93-bp in-frame deletion in glpD gene results in the glycerol-negative characteristic of biovar orientalis strains. Two kinds of point mutations in the napA gene may cause the nitrate reduction-negative characteristic in biovars mediaevalis and microtus, respectively. A 122-bp frameshift deletion in the araC gene may lead to the arabinose-negative phenotype of biovar microtus strains. Biovar microtus strains have a unique genomic profile of gene loss and pseudogene distribution, which most likely accounts for the human attenuation of this new biovar. Focused, hypothesis-based investigations on these specific genes will help delineate the determinants that enable this deadly pathogen to be virulent to humans and give insight into the evolution of Y. pestis and plague pathogenesis. Moreover, there may be the implications for development of biovar microtus strains as a potential vaccine.

DNA microarray analysis of genome dynamics in Yersinia pestis: insights into bacterial genome microevolution and niche adaptation

Genomics research provides an unprecedented opportunity for us to probe into the pathogenicity and evolution of the world's most deadly pathogenic bacterium, Yersinia pestis, in minute detail. In our present work, extensive microarray analysis in conjunction with PCR validation revealed that there are considerable genome dynamics, due to gene acquisition and loss, in natural populations of Y. pestis. We established a genomotyping system to group homologous isolates of Y. pestis, based on profiling or gene acquisition and loss in their genomes, and then drew an outline of parallel microevolution of the Y. pestis genome. The acquisition of a number of genomic islands and plasmids most likely induced Y. pestis to evolve rapidly from Yersinia pseudotuberculosis to a new, deadly pathogen. Horizontal gene acquisition also plays a key role in the dramatic evolutionary segregation of Y. pestis lineages (biovars and genomovars). In contrast to selective genome expansion by gene acquisition, genome reduction occurs in Y. pestis through the loss of DNA regions. We also theorized about the links between niche adaptation and genome microevolution. The transmission, colonization, and expansion of Y. pestis in the natural foci of endemic plague are parallel and directional and involve gradual adaptation to the complex of interactions between the environment, the hosts, and the pathogen itself. These adaptations are based on the natural selections against the accumulation of genetic changes within genome. Our data strongly support that the modern plague originated from Yunnan Province in China, due to the arising of biovar orientalis from biovar antiqua rather than mediaevalis.

Tandem repeats analysis for the high resolution phylogenetic analysis of Yersinia pestis

Background

Yersinia pestis, the agent of plague, is a young and highly monomorphic species. Three biovars, each one thought to be associated with the last three Y. pestis pandemics, have been defined based on biochemical assays. More recently, DNA based assays, including DNA sequencing, IS typing, DNA arrays, have significantly improved current knowledge on the origin and phylogenetic evolution of Y. pestis. However, these methods suffer either from a lack of resolution or from the difficulty to compare data. Variable number of tandem repeats (VNTRs) provides valuable polymorphic markers for genotyping and performing phylogenetic analyses in a growing number of pathogens and have given promising results for Y. pestis as well.

Results

In this study we have genotyped 180 Y. pestis isolates by multiple locus VNTR analysis (MLVA) using 25 markers. Sixty-one different genotypes were observed. The three biovars were distributed into three main branches, with some exceptions. In particular, the Medievalis phenotype is clearly heterogeneous, resulting from different mutation events in the napA gene. Antiqua strains from Asia appear to hold a central position compared to Antiqua strains from Africa. A subset of 7 markers is proposed for the quick comparison of a new strain with the collection typed here. This can be easily achieved using a Web-based facility, specifically set-up for running such identifications.

Conclusion

Tandem-repeat typing may prove to be a powerful complement to the existing phylogenetic tools for Y. pestis. Typing can be achieved quickly at a low cost in terms of consumables, technical expertise and equipment. The resulting data can be easily compared between different laboratories. The number and selection of markers will eventually depend upon the type and aim of investigations.

Intraspecific diversity of Yersinia pestis

Increased interest in the pathogenic potential of Yersinia pestis has emerged because of the potential threats from bioterrorism. Pathogenic potential is based on genetic factors present in a population of microbes, yet most studies evaluating the role of specific genes in virulence have used a limited number of strains. For Y. pestis this issue is complicated by the fact that most strains available for study in the Americas are clonally derived and thus genetically restricted, emanating from a strain of Y. pestis introduced into the United States in 1902 via marine shipping and subsequent spread of this strain throughout North and South America. In countries from the former Soviet Union (FSU), Mongolia, and China there are large areas of enzootic foci of Y. pestis infection containing genetically diverse strains that have been intensely studied by scientists in these countries. However, the results of these investigations are not generally known outside of these countries. Here we describe the variety of methods used in the FSU to classify Y. pestis strains based on genetic and phenotypic variation and show that there is a high level of diversity in these strains not reflected by ones obtained from sylvatic areas and patients in the Americas.

Application of DNA microarrays to study the evolutionary genomics of Yersinia pestis and Yersinia pseudotuberculosis

Yersinia pestis, the causative agent of plague, diverged from Yersinia pseudotuberculosis, an enteric pathogen, an estimated 1500-20,000 years ago. Genetic characterization of these closely related organisms represents a useful model to study the rapid emergence of bacterial pathogens that threaten mankind. To this end, we undertook genome-wide DNA microarray analysis of 22 strains of Y. pestis and 10 strains of Y. pseudotuberculosis of diverse origin. Eleven Y. pestis DNA loci were deemed absent or highly divergent in all strains of Y. pseudotuberculosis. Four were regions of phage origin, whereas the other seven included genes encoding a vitamin B12 receptor and the insect toxin sepC. Sixteen differences were identified between Y. pestis strains, with biovar Antiqua and Mediaevalis strains showing most divergence from the arrayed CO92 Orientalis strain. Fifty-eight Y. pestis regions were specific to a limited number of Y. pseudotuberculosis strains, including the high pathogenicity island, three putative autotransporters, and several possible insecticidal toxins and hemolysins. The O-antigen gene cluster and one of two possible flagellar operons had high levels of divergence between Y. pseudotuberculosis strains. This study reports chromosomal differences between species, biovars, serotypes, and strains of Y. pestis and Y. pseudotuberculosis that may relate to the evolution of these species in their respective niches.

Characterisation of IS901 integration sites in the Mycobacterium avium genome

Data are presented on the identification and characterisation of 17 chromosomal integration loci of the insertion element IS901 in the Mycobacterium avium (cervine strain JD88/118) genome. Thirteen of these integration loci have been mapped to their corresponding positions on the M. avium strain 104 (an IS901(-) strain) genome (The Institute for Genome Research (TIGR) unfinished genome-sequencing project). Sequence data for both upstream and downstream sequence flanking regions were obtained for 12 insertion loci, while upstream sequence was obtained for five others. A consensus IS901 insertion target sequence compiled from all 17 integration sites was in broad agreement with earlier reports that were based on only two such loci. Analysis of IS901 integration site flanking sequences revealed that, like IS900 in M. avium subspecies paratuberculosis, IS901 inserts preferentially between a putative ribosome-binding sequence (RBS) and the translational start codon of an open reading frame (ORF). In BLAST X and BLAST P searches of the GenBank database, these ORFs were shown to share significant homologies with a number of other prokaryotic genes.

Genomic signature tags (GSTs): a system for profiling genomic DNA

Genomic signature tags (GSTs) are the products of a method we have developed for identifying and quantitatively analyzing genomic DNAs. The DNA is initially fragmented with a type II restriction enzyme. An oligonucleotide adaptor containing a recognition site for MmeI, a type IIS restriction enzyme, is then used to release 21-bp tags from fixed positions in the DNA relative to the sites recognized by the fragmenting enzyme. These tags are PCR-amplified, purified, concatenated, and then cloned and sequenced. The tag sequences and abundances are used to create a high-resolution GST sequence profile of the genomic DNA. GSTs are shown to be long enough for use as oligonucleotide primers to amplify adjacent segments of the DNA, which can then be sequenced to provide additional nucleotide information or used as probes to identify specific clones in metagenomic libraries. GST analysis of the 4.7-Mb Yersinia pestis EV766 genome using BamHI as the fragmenting enzyme and NlaIII as the tagging enzyme validated the precision of our approach. The GST profile predicts that this strain has several changes relative to the archetype CO92 strain, including deletion of a 57-kb region of the chromosome known to be an unstable pathogenicity island.

Conserved filamentous prophage in Escherichia coli O18:K1:H7 and Yersinia pestis biovar orientalis

Microbial virulence is known to emerge by horizontal gene transfer mechanisms. Here we describe the discovery of a novel filamentous prophage, designated CUS-1, which is integrated into the chromosomal dif homologue of the high-virulence clone Escherichia coli O18:K1:H7. An homologous chromosomal element (CUS-2) in Yersinia pestis biovar orientalis is integrated at the same relative location as CUS-1; both lysogenic E. coli and Y. pestis strains produce particles with properties expected of single-stranded DNA virions. CUS(phi) is epidemiologically correlated with the emergence of K1 strains with increased virulence and with the Y. pestis biovar responsible for the current (third) plague pandemic.

Genome plasticity in Yersinia pestis

Yersinia pestis, the causative agent of bubonic plague, emerged recently (<20000 years ago) as a clone of Yersinia pseudotuberculosis. There is scant evidence of genome diversity in Y. pestis, although it is possible to differentiate three biovars (antiqua, mediaevalis or orientalis) based on two biochemical tests. There are a few examples of restriction fragment length polymorphisms (RFLPs) within Y. pestis; however, their genetic basis is poorly understood. In this study, six difference regions (DFRs) were identified in Y. pestis, by using subtractive hybridization, which ranged from 4.6 to 19 kb in size. Four of the DFRs are flanked by insertion sequences, and their sequences show similarity to bacterial genes encoding proteins for flagellar synthesis, ABC transport, insect toxicity and bacteriophage functions. The presence or absence of these DFRs (termed the DFR profile) was demonstrated in 78 geographically diverse strains of Y. pestis. Significant genome plasticity was observed among these strains and suggests the acquisition and deletion of these DNA regions during the recent evolution of Y. pestis. Y. pestis biovar orientalis possesses DFR profiles that are different from antiqua and mediaevalis biovars, reflecting the recent origins of this biovar. Whereas some DFR profiles are specific for antiqua and mediaevalis, some DFR profiles are shared by both biovars. Furthermore, the progenitor of Y. pestis, Y. pseudotuberculosis (an enteric pathogen), possesses its own DFR profile. The DFR profiles detailed here demonstrate genome plasticity within Y. pestis, and they imply evolutionary relationships among the three biovars of Y. pestis, as well as between Y. pestis and Y. pseudotuberculosis.