Search

The life cycle of Yersinia pestis has changed a lot to adapt to flea-borne transmission since it evolved from an enteric pathogen, Yersinia pseudotuberculosis. Small insertions and deletions (indels), especially frameshift mutations, can have major effects on phenotypes and contribute to virulence and host adaptation through gene disruption and inactivation. Here, we analyzed 365 Y. pestis genomes and identified 2,092 genome-wide indels on the core genome. As recently reported in Mycobacterium tuberculosis, we also detected "indel pockets" in Y. pestis, with average complexity scores declining around indel positions, which we speculate might also exist in other prokaryotes. Phylogenic analysis showed that indel-based phylogenic tree could basically reflect the phylogenetic relationships of major phylogroups in Y. pestis, except some inconsistency around the Big Bang polytomy. We observed 83 indels arising in the trunk of the phylogeny, which played a role in accumulation of pseudogenes related to key metabolism and putatively pathogenicity. We also discovered 32 homoplasies at the level of phylogroups and 7 frameshift scars (i.e., disrupted reading frame being rescued by a second frameshift). Additionally, our analysis showed evidence of parallel evolution at the level of genes, with

The great gerbil (Rhombomys opimus) is a social rodent living in permanent, complex burrow systems distributed throughout Central Asia, where it serves as the main host of several important vector-borne infectious pathogens including the well-known plague bacterium (Yersinia pestis). Here, we present a continuous annotated genome assembly of the great gerbil, covering over 96% of the estimated 2.47-Gb genome. Taking advantage of the recent genome assemblies of the sand rat (Psammomys obesus) and the Mongolian gerbil (Meriones unguiculatus), comparative immunogenomic analyses reveal shared gene losses within TLR gene families (i.e., TLR8, TLR10, and the entire TLR11-subfamily) for Gerbillinae, accompanied with signs of diversifying selection of TLR7 and TLR9. Most notably, we find a great gerbil-specific duplication of the MHCII DRB locus. In silico analyses suggest that the duplicated gene provides high peptide binding affinity for Yersiniae epitopes as well as Leishmania and Leptospira epitopes, putatively leading to increased capability to withstand infections by these pathogens. Our study demonstrates the power of whole-genome sequencing combined with comparative genomic analyses to gain deeper insight into the immunogenomic landscape of the great gerbil and its close relatives.

In this study, a detailed genetic dissection of the huge and complex blaNDM-carrying genetic elements and their related mobile genetic elements was performed in Enterobacteriaceae. An extensive comparison was applied to 12 chromosomal genetic elements, including six sequenced in this study and the other six from GenBank. These 12 genetic elements were divided into five groups: a novel IME Tn6588; two related IMEs Tn6523 (SGI1) and Tn6589; four related ICEs Tn6512 (R391), Tn6575 (ICEPvuChnBC22), Tn6576, and Tn6577; Tn7 and its derivatives Tn6726 and 40.7-kb Tn7-related element; and two related IMEs Tn6591 (GIsul2) and Tn6590. At least 51 resistance genes, involved in resistance to 18 different categories of antibiotics and heavy metals, were found in these 12 genetic elements. Notably, Tn6576 carried another ICE Tn6582. In particular, the six blaNDM-carrying genetic elements Tn6588, Tn6589, Tn6575, Tn6576, Tn6726, and 40.7-kb Tn7-related element contained large accessory multidrug resistance (MDR) regions, each of which had a very complex mosaic structure that comprised intact or residual mobile genetic elements including insertion sequences, unit or composite transposons, integrons, and putative resistance units. Core blaNDM genetic environments manifested as four different Tn125 derivatives and, notably, two or more copies of relevant Tn125 derivatives were found in each of Tn6576, Tn6588, Tn6589, and 40.7-kb Tn7-related element. The huge and complex blaNDM-carrying genetic elements were assembled from complex transposition and homolog recombination. Firstly identified were eight novel mobile elements, including three ICEs Tn6576, Tn6577, and Tn6582, two IMEs, Tn6588 and Tn6589, two composite transposons Tn6580a and Tn6580b, and one integron In1718.

Pathogens may elicit a high selective pressure on hosts and can alter genetic diversity over short evolutionary timescales. Intraspecific variation in immune response can be observed as variable survivability from specific infections. The great gerbil (Rhombomys opimus) is a rodent plague host with a heterogenic but highly resistant phenotype. Here, we investigate if the most plague-resistant phenotypes are linked to genomic differences between survivors and susceptible individuals by exposure of wild-caught great gerbils from Northwest China to plague (Yersinia pestis). Whole genome sequencing of ten survivors and ten moribund individuals revealed a low genome-wide mean divergence, except for a subset of genomic regions that showed elevated differentiation. Gene ontology (GO) analysis of candidate genes within regions of increased differentiation, demonstrated enrichment of pathways involved in transcription and translation and their regulation), as well as genes directly involved in immune functions, cellular metabolism and the regulation of apoptosis. Differential RNA expression analysis revealed that the early activated great gerbil immune response to plague consisted of classical components of the innate immune system. Our approach combining challenge experiments with transcriptomics and population level sequencing, provides new insight into the genetic background of plague-resistance and confirms its complex nature, most likely involving multiple genes and pathways of both the immune system and regulation of basic cellular functions.

Yunnan Province, China is thought to be the original source of biovar Orientalis of Yersinia pestis, the causative agent of the third plague pandemic that has spread globally since the end of the 19th century. Although encompassing a large area of natural plague foci, Y. pestis strains have rarely been found in live rodents during surveillance in Yunnan, and most isolates are from rodent corpses and their fleas. In 2017, 10 Y. pestis strains were isolated from seven live rodents and three fleas in Heqing County of Yunnan. These strains were supposed to have low virulence to local rodents Eothenomys miletus and Apodemus chevrieri because the rodents were healthy and no dead animals were found in surrounding areas, as had occurred in previous epizootic disease. We performed microscopic and biochemical examinations of the isolates, and compared their whole-genome sequences and transcriptome with those of 10 high virulence Y. pestis strains that were isolated from nine rodents and one parasitic flea in adjacent city (Lijiang). We analyzed the phenotypic, genomic, and transcriptomic characteristics of live rodent isolates. The isolates formed a previously undefined monophyletic branch of Y. pestis that was named 1.IN5. Six SNPs, two indels, and one copy number variation were detected between live rodent isolates and the high virulence neighbors. No obvious functional consequence of these variations was found according to the known annotation information. Among genes which expression differential in the live rodent isolates compared to their high virulent neighbors, we found five iron transfer related ones that were significant up-regulated (| log2 (FC) | > 1, p.adjust < 0.05), indicating these genes may be related to the low-virulence phenotype. The novel genotype of Y. pestis reported here provides further insights into the evolution and spread of plague as well as clues that may help to decipher the virulence mechanism of this notorious pathogen.

To explore the genetic diversity of Yersinia pestis in Erenhot, China, and their relationship with Mongolia strains, we collected and sequenced three Y. pestis strains from Erenhot, China, in 2018. Here, we report the draft genome sequences of three Y. pestis bv. Medievalis strains belonging to the 2.MED phylogroup that were circulating in Meriones unguiculatus populations., To explore the genetic diversity of Yersinia pestis strains in Erenhot, China, and their relationship with Mongolian strains, we collected and sequenced three Y. pestis strains from Erenhot, China, in 2018. Here, we report the draft genome sequences of three Y. pestis bv. Medievalis strains belonging to the 2.MED phylogroup that were circulating in Meriones unguiculatus populations.

Here, we report the draft genome sequences of two Yersinia pestis bv. Antiqua strains, belonging to the 3.ANT phylogroup, that were isolated in Mongolia and were circulating in marmot populations.

This study aimed to genetically characterize two fully-sequenced novel IncFII-type multidrug resistant (MDR) plasmids, p0716-KPC and p12181-KPC, recovered from two different clinical Klebsiella pneumoniae isolates. p0716-KPC and p12181-KPC had a very similar genomic content. The backbones of p0716-KPC/p12181-KPC contained two different replicons (belonging to a novel IncFII subtype and the Rep_3 family), the IncFIIK and IncFIIY maintenance regions, and conjugal transfer gene sets from IncFIIK-type plasmids and unknown origins. p0716-KPC and p12181-KPC carried similar three accessory resistance regions, namely ΔTn6209, a MDR region, and the blaKPC-2 region. Resistance genes blaKPC-2, mph(A), strAB, aacC2, qacEΔ1, sul1, sul2, and dfrA25, which are associated with transposons, integrons, and insertion sequence-based mobile units, were located in these accessory regions. p0716-KPC carried two additional resistance genes: aphA1a and blaTEM-1. Together, our analyses showed that p0716-KPC and p12181-KPC belong to a novel IncFII subtype and display a complex chimeric nature, and that the carbapenem resistance gene blaKPC-2 coexists with a lot of additional resistance genes on these two plasmids.

Incompatibility group IncHI plasmids are important vectors of antibiotic resistance in Enterobacteriaceae. In this study, a scheme for typing IncHI into five separately clustering subgroups, including previously designated IncHI1–3 as well as IncHI4–5, was proposed based on sequenced IncHI plasmids. The complete nucleotide sequences of the IncHI2 plasmids pT5282-mphA and p112298-catA and the IncHI5 plasmid pYNKP001-dfrA from clinical Enterobacter cloacae , Citrobacter freundii and Raoultella ornithinolytica isolates, respectively, were determined and were compared with IncHI2 and IncHI5 reference plasmids. Considerable genetic conservation was observed within the backbone sequences of each of the IncHI2 and IncHI5 subgroups, but the backbone sequences of the two subgroups were dramatically different from each other. However, the conjugal transfer regions tra1 and tra2 as well as the tellurium resistance gene cluster ter were present in all five plasmids. A number of accessory regions associated with integrons, transposons and insertion sequence-based mobile elements have been inserted at various sites of the plasmid backbones, among which were several large regions harbouring genes conferring resistance to multiple classes of antibiotics. Data generated from this study provide us with a deeper understanding of the diversification of IncHI-type resistance plasmids.

Objectives To dissect genomic features of IncpRBL16 plasmids from Pseudomonas. Methods An extensive genomic comparison was applied to all 17 available sequenced IncpRBL16 plasmids, including 8 sequenced in this study and another 2 sequenced in two of our previous studies. Results Conserved IncpRBL16 backbone markers repAIncpRBL16 together with its iterons, parB2–parA, che, pil and ter were present in all 17 plasmids. At least 18 regions or sites across IncpRBL16 genomes exhibited major modular differences, including insertion of accessory modules, deletion of backbone regions surrounding insertion sites and substitution of multiple-gene backbone regions. Ten plasmids carried a sole IncpRBL16 replicon, while exogenous acquisition of an auxiliary replicon (located in an accessory module) besides the primary IncpRBL16 replicon was observed in each of the remaining seven plasmids. The 17 IncpRBL16 plasmids carried at least 71 different accessory modules, notably including Tn1403-related regions, Tn7-family transposons, Tn6571-family transposons, integrative and conjugative elements, and integrative and mobilizable elements. There were a total of 40 known resistance genes, which were involved in resistance to 15 categories of antibiotics and heavy metals, notably including blaIMP-9, blaIMP-45, blaVIM-2, blaDIM-2, blaOXA-246, blaPER-1, aphA and armA. Conclusions Different IncpRBL16 plasmids contain different profiles of accessory modules and thus diverse collections of resistance genes. To the best of our knowledge, this is the first report of fully sequenced blaOXA-246-carrying (p12939-PER) and blaPER-1-carrying (p12939-PER and pA681-IMP) IncpRBL16 plasmids and also that of 14 novel (first identified in this study) and additionally 31 newly named (first designated in this study, but with previously determined sequences) mobile elements.

Many genes in the bacterial pathogen Yersinia pestis, the causative agent of three plague pandemics, remain uncharacterized, greatly hampering the development of measures for plague prevention and control. Clustered regularly interspaced short palindromic repeat interference (CRISPRi) has been shown to be an effective tool for gene knockdown in model bacteria. In this system, a catalytically dead Cas9 (dCas9) and a small guide RNA (sgRNA) form a complex, binding to the specific DNA target through base pairing, thereby impeding RNA polymerase binding and causing target gene repression. Here, we introduce an optimized CRISPRi system using Streptococcus pyogenes Cas9-derived dCas9 for gene knockdown in Y. pestis. Multiple genes harbored on either the chromosome or plasmids of Y. pestis were efficiently knocked down (up to 380-fold) in a strictly anhydrotetracycline-inducible manner using this CRISPRi approach. Knockdown of hmsH (responsible for biofilm formation) or cspB (encoding a cold shock protein) resulted in greatly decreased biofilm formation or impaired cold tolerance in in vitro phenotypic assays. Furthermore, silencing of the virulence-associated genes yscB or ail using this CRISPRi system resulted in attenuation of virulence in HeLa cells and mice similar to that previously reported for yscB and ail null mutants. Taken together, our results confirm that this optimized CRISPRi system can reversibly and efficiently repress the expression of target genes in Y. pestis, providing an alternative to conventional gene knockdown techniques, as well as a strategy for high-throughput phenotypic screening of Y. pestis genes with unknown functions. IMPORTANCEYersiniapestis is a lethal pathogen responsible for millions of human deaths in history. It has also attracted much attention for potential uses as a bioweapon or bioterrorism agent, against which new vaccines are desperately needed. However, many Y. pestis genes remain uncharacterized, greatly hampering the development of measures for plague prevention and control. Clustered regularly interspaced short palindromic repeat interference (CRISPRi) has been successfully used in a variety of bacteria in functional genomic studies, but no such genetic tool has been reported in Y. pestis. Here, we systematically optimized the CRISPRi approach for use in Y. pestis, which ultimately repressed target gene expression with high efficiency in a reversible manner. Knockdown of functional genes using this method produced phenotypes that were readily detected by in vitro assays, cell infection assays, and mouse infection experiments. This is a report of a CRISPRi approach in Y. pestis and highlights the potential use of this approach in high-throughput functional genomics studies of this pathogen.

We employed a stepwise selection model for investigating the dynamics of antibiotic-resistant variants in Escherichia coli K-12 treated with increasing concentrations of ciprofloxacin (CIP). Firstly, we used Sanger sequencing to screen the variations in the fluoquinolone target genes, then, employed Illumina NGS sequencing for amplicons targeted regions with variations. The results demonstrated that variations G81C in gyrA and K276N and K277L in parC are standing resistance variations (SRVs), while S83A and S83L in gyrA and G78C in parC were emerging resistance variations (ERVs). The variants containing SRVs and/or ERVs were selected successively based on their sensitivities to CIP. Variant strain 1, containing substitution G81C in gyrA, was immediately selected following ciprofloxacin exposure, with obvious increases in the parC SRV, and parC and gyrA ERV allele frequencies. Variant strain 2, containing the SRVs, then dominated the population following a 20× increase in ciprofloxacin concentration, with other associated allele frequencies also elevated. Variant strains 3 and 4, containing ERVs in gyrA and parC, respectively, were then selected at 40× and 160× antibiotic concentrations. Two variants, strains 5 and 6, generated in the selection procedure, were lost because of higher fitness costs or a lower level of resistance compared with variants 3 and 4. For the second induction, all variations/indels were already present as SRVs and selected out step by step at different passages. Whatever the first induction or second induction, our results confirmed the soft selective sweep hypothesis and provided critical information for guiding clinical treatment of pathogens containing SRVs.

The great gerbil (Rhombomys opimus) is a social rodent living in permanent, complex burrow systems distributed throughout Central Asia, where it serves as the main host of several important vector-borne infectious diseases and is defined as a key reservoir species for plague (Yersinia pestis). Studies from the wild have shown that the great gerbil is largely resistant to plague but the genetic basis for resistance is yet to be determined. Here, we present a highly contiguous annotated genome assembly of great gerbil, covering over 96 % of the estimated 2.47 Gb genome. Comparative genomic analyses focusing on the immune gene repertoire, reveal shared gene losses withinTLRgene families (i.e.TLR8, TLR10and all members ofTLR11-subfamily) for the Gerbillinae lineage, accompanied with signs of diversifying selection ofTLR7andTLR9. Most notably, we find a great gerbil-specific duplication of theMHCII DRBlocus.In silicoanalyses suggest that the duplicated gene provides high peptide binding affinity forYersiniaeepitopes. The great gerbil genome provides new insights into the genomic landscape that confers immunological resistance towards plague. The high affinity forYersiniaepitopes could be key in our understanding of the high resistance in great gerbils, putatively conferring a faster initiation of the adaptive immune response leading to survival of the infection. Our study demonstrates the power of studying zoonosis in natural hosts through the generation of a genome resource for further comparative and experimental work on plague survival and evolution of host-pathogen interactions.

Yersinia pestis is transmitted from fleas to rodents when the bacterium develops an extensive biofilm in the foregut of a flea, starving it into a feeding frenzy, or, alternatively, during a brief period directly after feeding on a bacteremic host. These two transmission modes are in a trade-off regulated by the amount of biofilm produced by the bacterium. Here by investigating 446 global isolated Y. pestis genomes, including 78 newly sequenced isolates sampled over 40 years from a plague focus in China, we provide evidence for strong selection pressures on the RNA polymerase ω-subunit encoding gene rpoZ. We demonstrate that rpoZ variants have an increased rate of biofilm production in vitro, and that they evolve in the ecosystem during colder and drier periods. Our results support the notion that the bacterium is constantly adapting—through extended phenotype changes in the fleas—in response to climate-driven changes in the niche., Yersinia pestis, the causative agent of plague, can change its biofilm production to influence the dynamics of flea-borne transmission. Here, the authors sequence Y. pestis isolates sampled over 40 years in China and show evidence for climate-associated selection on rpoZ to increase biofilm production.

We investigated global patterns of variation in 157 whole-genome sequences of Vibrio parahaemolyticus, a free-living and seafood associated marine bacterium. Pandemic clones, responsible for recent outbreaks of gastroenteritis in humans, have spread globally. However, there are oceanic gene pools, one located in the oceans surrounding Asia and another in the Mexican Gulf. Frequent recombination means that most isolates have acquired the genetic profile of their current location. We investigated the genetic structure in the Asian gene pool by calculating the effective population size in two different ways. Under standard neutral models, the two estimates should give similar answers but we found a 27-fold difference. We propose that this discrepancy is caused by the subdivision of the species into a hundred or more ecotypes which are maintained stably in the population. To investigate the genetic factors involved, we used 51 unrelated isolates to conduct a genome-wide scan for epistatically interacting loci. We found a single example of strong epistasis between distant genome regions. A majority of strains had a type VI secretion system associated with bacterial killing. The remaining strains had genes associated with biofilm formation and regulated by cyclic dimeric GMP signaling. All strains had one or other of the two systems and none of isolate had complete complements of both systems, although several strains had remnants. Further "top down" analysis of patterns of linkage disequilibrium within frequently recombining species will allow a detailed understanding of how selection acts to structure the pattern of variation within natural bacterial populations.

This work presents the complete nucleotide sequences of p0801-IMP from Klebsiella pneumoniae, p7121-IMP from K. oxytoca, and p17285-IMP from Citrobacter freundii, which are recovered from three different cases of nosocomial infection. These three plasmids represent the first fully sequenced blaIMP-carrying IncN2 plasmids. Further comparative genomics analysis of all the five integron-carrying IncN2 plasmids p0801-IMP, p7121-IMP, p17285-IMP, pJIE137, and p34983-59.134kb indicates that they possess conserved IncN2 backbones with limited genetic variations with respect to gene content and organization. Four class 1 integrons (blaIMP-1-carrying In1223 in p0801-IMP/p7121-IMP, blaIMP-8-carrying In655 in p17285-IMP, In27 in pJIE137, and In1130 in p34983-59.134kb), two insertion sequence-based transposition units (ISEcp1-orfRA1-14 in p17285-IMP, and ISEcp1-blaCTX-M-62-Δorf477-orfRA1-14 in pJIE137), and a novel Tn1696-related transposon Tn6325 carrying In1130 in p34983-59.134kb are indentified in the plasmid accessory regions. In1223 and In655 represent ancestral Tn402-associated integrons, while In27 and In1130 belong to complex class 1 integrons. The relatively small IncN2 backbones are able to integrate different mobile elements which carry various resistance markers, promoting the accumulation and spread of antimicrobial resistance genes among enterobacterial species.

Plague, one of the most devastating infectious diseases in human history, is caused by the bacterial species Yersinia pestis. A live attenuated Y. pestis strain (EV76) has been widely used as a plague vaccine in various countries around the world. Here we compared the whole genome sequence of an EV76 strain used in China (EV76-CN) with the genomes of Y. pestis wild isolates to identify genetic variations specific to the EV76 lineage. We identified 6 SNPs and 6 Indels (insertions and deletions) differentiating EV76-CN from its counterparts. Then, we screened these polymorphic sites in 28 other strains of EV76 lineage that were stored in different countries. Based on the profiles of SNPs and Indels, we reconstructed the parsimonious dissemination history of EV76 lineage. This analysis revealed that there have been at least three independent imports of EV76 strains into China. Additionally, we observed that the pyrE gene is a mutation hotspot in EV76 lineages. The fine comparison results based on whole genome sequence in this study provide better understanding of the effects of laboratory passages on the accumulation of genetic polymorphisms in plague vaccine strains. These variations identified here will also be helpful in discriminating different EV76 derivatives.

Yersinia pestis is one of the most dangerous bacterial pathogens. PhoP and cyclic AMP receptor protein (CRP) are global regulators of Y. pestis, and they control two distinct regulons that contain multiple virulence-related genes. The PhoP regulator and its cognate sensor PhoQ constitute a two-component regulatory system. The regulatory activity of CRP is triggered only by binding to its cofactor cAMP, which is synthesized from ATP by adenylyl cyclase (encoded by cyaA). However, the association between the two regulatory systems PhoP/PhoQ and CRP-cAMP is still not understood for Y. pestis. In the present work, the four consecutive genes YPO1635, phoP, phoQ, and YPO1632 were found to constitute an operon, YPO1635-phoPQ-YPO1632, transcribed as a single primary RNA, whereas the last three genes comprised another operon, phoPQ-YPO1632, transcribed with two adjacent transcriptional starts. Through direct PhoP-target promoter association, the transcription of these two operons was stimulated and repressed by PhoP, respectively; thus, both positive autoregulation and negative autoregulation of PhoP/PhoQ were detected. In addition, PhoP acted as a direct transcriptional activator of crp and cyaA. The translational/transcriptional start sites, promoter −10 and −35 elements, PhoP sites, and PhoP box-like sequences were determined for these PhoP-dependent genes, providing a map of the PhoP-target promoter interaction. The CRP and PhoP regulons have evolved to merge into a single regulatory cascade in Y. pestis because of the direct regulatory association between PhoP/PhoQ and CRP-cAMP.

Yersinia pestis is among the most dangerous human pathogens, and systematic research of this pathogen is important in bacterial pathogenomics research. To fully interpret the biological functions, physiological characteristics, and pathogenesis of Y. pestis, a comprehensive annotation of its entire genome is necessary. The emergence of omics-based research has brought new opportunities to better annotate the genome of this pathogen. Here, the complete genome of Y. pestis strain 91001 was reannotated using genomics and proteogenomics data. One hundred and thirty-seven unreliable coding sequences were removed, and 41 homologous genes were relocated with their translational initiation sites, while the functions of seven pseudogenes and 392 hypothetical genes were revised. Moreover, annotations of noncoding RNAs, repeat sequences, and transposable elements have also been incorporated. The reannotated results are freely available at http://tody.bmi.ac.cn.

Here, we report the draft genome sequence of “ Paramesorhizobium deserti ” A-3-E T , a strain isolated from the Taklimakan Desert of Xinjiang, China, which is resistant to multiple β-lactam antibiotics and other antibiotics (kanamycin, erythromycin, streptomycin, etc.) as well.

The catalase or catalase-peroxidase activity commonly exists in many pathogens and plays an important role in resisting the oxidative burst of phagocytes helping the pathogen persistently colonize in the host. Yersinia pestis is a facultative pathogen and the causative agent of plague. KatY has been identified as a thermosensing antigen with modest catalase activity in this pathogen. Here Y. pestis KatA and KatY were experimentally confirmed as a monofunctional catalase and bifunctional catalase-peroxidase, respectively. Their expression induced by H2O2 was proven to be mediated by the oxidative regulator, OxyR. Expression of KatA changed with growth phases and was crucial to its traditional physiological role in protecting Y. pestis cells against toxicity of exogenous H2O2. KatY was regulated by temperature and H2O2, two major elements of phagolysosomal microenvironments. Consistent with the above results, gene expression of katY increased significantly during intracellular growth of Y. pestis compared with that in vitro growth. However, a DeltakatY mutant was fully virulent to mice, suggesting that KatY is not required for Y. pestis virulence.

Using DNA microarray analysis, mRNA levels from wild-type Yersinia pestis cells treated with the iron chelator 2,2′-dipyridyl were compared with those supplemented with excessive iron, and subsequent to this, gene expression in the fur mutant was compared with that in the wild-type strain under iron rich conditions. The microarray analysis revealed many iron transport or storage systems that had been induced in response to the iron starvation, which is mediated by the Fur protein, using the iron as a co-repressor. The iron–Fur complex also affected some genes involved in various non-iron functions (ribonucleoside-diphosphate reductase, membrane proteins, electron transport and oxidative defense, etc.). The Fur protein still participated in the regulation of genes involved in broad cellular processes (virulence factors, pesticin activity, haemin storage and many proteins with unknown functions) that were not affected by iron depletion conditions. In addition to its classical negative regulatory activities, the Fur protein activates gene transcription. Using bioinformatics tools, we were able to predict the Y. pestis Fur box sequence that was clearly the over-presented motif in the promoter regions of members of the iron–Fur modulon.

Yersinia pestis, the causative agent of bubonic and pneumonic plague, has been classified into four biovars: Antiqua, Mediaevalis, Orientalis and Microtus. Although the entire genome sequences of three Y. pestis strains, CO92, KIM and 91001, of biovar Orientalis, Mediaevalis and Microtus, respectively, have been decoded, the genome sequence of the biovar Antiqua strain is unknown. In an initial effort to find Antiqua-specific sequences, suppression subtractive hybridization (SSH) was performed and four different regions (DFRs) were identified. Among the four DFRs, only DFR4 was specific to the tester (strain 49006, biovar Antiqua). PCR demonstrated that DFR4 was present only in 57 of 60 Antiqua strains from the Marmota baibacina-Spermophilus undulates plague focus in the Tianshan Mountains (focus B) and in three strains of Y. pseudotuberculosis (serotypes I and II), showing that not all Antiqua strains had DFR4. Five DFR profiles were identified based on the presence or absence of these four DFRs in 636 strains of Y. pestis from 10 plague foci in China.

Plague is a natural focus-based disease, and for better understanding of this disease it is crucial to determine the molecular mechanisms of its pathogen, Yersinia pestis, for adapting to different foci. Gene inactivation, loss and acquisition are the main mechanisms that contribute to a pathogen's fitness. Determination of the whole-genome sequences of three Y. pestis strains, CO92, KIM and 91001, provided a good opportunity to probe into its genome in minute detail. Many genetic variations were found between the three strains. The present work focused on adaptive microevolutionary analysis of Y. pestis from different natural plague foci in China based on pseudogene profiles. Twenty-four mutations that led to inactivation in the corresponding genes were analysed, and a PCR-based screening method was employed to investigate the distribution of these mutations among Y. pestis isolates from different foci and also among seven strains of Yersinia pseudotuberculosis. It was found that Y. pestis isolates from the same focus had identical mutation profiles, and 260 isolates of Y. pestis were divided into eight genotypes, while Y. pseudotuberculosis harboured wild-type alleles for all the mutations. The isolates of three known biovars were grouped into distinct branches in the phylogenetic tree, which supports the proposition that biovars mediaevalis and orientalis directly arose from biovar antiqua individually. The constructed phylogenetic tree suggests that the isolates from focus B should be the oldest lineage of Y. pestis in China except for isolates from foci L and M, which might be a special lineage of Y. pestis and originated differently to the others.

A PCR-based screening method was used to study the genetic variations of the pgm locus among natural isolates of Yersinia pestis from China. Our results indicate that genetic variations in the pgm locus are well correlated with biovars of Y. pestis and plague foci, suggesting that the pgm locus plays a role in Y. pestis adaptation to its environment. The gene encoding two-component regulatory system sensor kinase became a pseudogene in all strains of biovar Orientalis due to a thymidine deletion, while it is intact in all the strains of the other biovars. Only strains from Foci H and L are the same as Yersinia pseudotuberculosis in that they have an intact transmembrane helix in the sensor kinase protein, which is lost in all the other strains because of the 18 bp in-frame deletion. The IS100 element that flanks the 39 terminus of the pgm locus was inserted into the chromosome during the within-species microevolution of Y. pestis, which is absent in strains from Foci G, H and L and also in Y. pseudotuberculosis. This fact indicates that the strains from these three foci are of an older lineage of Chinese Y. pestis. It is this IS100 element's absence that maintained high stability of the pgm locus in the Y. pestis strains from these three foci. The IS285 element insertion in the pigmentation segment and the IS100 element insertion in the downstream flanking region of the pgm locus are only present in strains from Foci H and L. The flanking region outside the 59 terminus of the upstream IS100 element is identical in the strains from these two foci, which is different in the other strains. All of these unique characteristics suggest that they are of a special lineage of Chinese Y. pestis.

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.

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.

In this study, 159 Legionella pneumophila strains isolated from various natural and artificial water sources in Guangzhou and Jiangmen, China, were subjected to genotyping by the sequence-based typing (SBT) scheme. These isolates were assigned into 53 sequence types (STs) (50 STs with seven loci data and three unidentified STs with incomplete loci profiles) with ST1 as the dominant one (14.5%), and the index of diversity (IOD) was 0.950. Eight new alleles and 34 new STs were reported here. Notably, most of the newly identified STs with seven loci data (24/34) contained no new allele, implying frequent recombination events in L. pneumophila. Five intragenic recombination events were identified in the concatenated sequences of seven loci. The diversity of STs in natural environmental isolates (41 STs, IOD = 0.956) is higher than that of artificial environmental ones (17 STs, IOD = 0.824). The ST patterns varied in isolates from these two sources: the most common STs from artificial water sources, ST1 and ST752 (39.2% and 13.7%), were only occasionally isolated from natural water sources (2.9% and 3.8%, respectively); while the predominant STs from natural water sources, ST1048, ST739 and ST1267 (15.2%, 6.7% and 6.7%), were less frequently seen in artificial environments (2.0%, 0% and 0%, respectively). We also found out that Legionnaires’ disease associated STs might be more frequently isolated in artificial environments than in natural ones. Our data revealed remarkable genetic diversity of L. pneumophila isolates from environmental water systems of Guangzhou and Jiangmen, and the different ST distribution patterns between natural water and artificial water sources as well. 2014 Published by Elsevier B.V.

Bacterial pathogens have evolved various mechanisms to modulate host immune responses for successful infection. In this study, RNA-sequencing technology was used to analyze the responses of human monocytes THP1 to Yersinia pestis infection. Over 6000 genes were differentially expressed over the 12 h infection. Kinetic responses of pathogen recognition receptor signaling pathways, apoptosis, antigen processing, and presentation pathway and coagulation system were analyzed in detail. Among them, RIG-I-like receptor (RLR) signaling pathway, which was established for antiviral defense, was significantly affected. Mice lacking MAVS, the adaptor of the RLR signaling pathway, were less sensitive to infection and exhibited lower IFN-β production, higher Th1-type cytokines IFN-γ and IL-12 production, and lower Th2-type cytokines IL-4 and IL-13 production in the serum compared with wild-type mice. Moreover, infection of pathogenic bacteria other than Y. pestis also altered the expression of the RLR pathway, suggesting that the response of RLR pathway to bacterial infection is a universal mechanism.

The genetic diversity of Yersinia pestis , the etiologic agent of plague, is extremely limited because of its recent origin coupled with a slow clock rate. Here we identified 2,326 SNPs from 133 genomes of Y. pestis strains that were isolated in China and elsewhere. These SNPs define the genealogy of Y. pestis since its most recent common ancestor. All but 28 of these SNPs represented mutations that happened only once within the genealogy, and they were distributed essentially at random among individual genes. Only seven genes contained a significant excess of nonsynonymous SNP, suggesting that the fixation of SNPs mainly arises via neutral processes, such as genetic drift, rather than Darwinian selection. However, the rate of fixation varies dramatically over the genealogy: the number of SNPs accumulated by different lineages was highly variable and the genealogy contains multiple polytomies, one of which resulted in four branches near the time of the Black Death. We suggest that demographic changes can affect the speed of evolution in epidemic pathogens even in the absence of natural selection, and hypothesize that neutral SNPs are fixed rapidly during intermittent epidemics and outbreaks.

The spread of the bla NDM-1 gene is gaining worldwide attentions. This gene is usually carried by large plasmids and has been discovered in diverse bacteria since it was originally found in Klebsiella pneumoniae . Here we report the complete sequences of a bla NDM-1 -bearing plasmid, pNDM-BJ01, and its variant, pNDM-BJ02, isolated from clinical Acinetobacter lwoffii strains. The plasmid pNDM-BJ01 is 47.3 kb in size and cannot be classified into any known plasmid incompatibility group, thus representing a novel plasmid with an unknown maintenance mechanism. This plasmid contains both a bla NDM-1 gene and a type IV secretion system (T4SS) gene cluster. The T4SS is assigned to the P-type T4SS group, which usually encode a short, rigid pilus, and the bla NDM-1 gene is located within a composite transposon flanked by two insertion elements of IS Aba125 . Plasmid pNDM-BJ02 is nearly identical to pNDM-BJ01 except that one copy of the IS Aba125 element is missing, and it is therefore regarded as a variant of pNDM-BJ01. Sequence alignment indicated that this bla NDM-1 -containing composite transposon, which can also be captured by other mobile elements, was probably a product of multiple recombination events and can move as a whole by transposition.

Here we report the draft genome sequence of one Acinetobacter genomic species 3 strain, D499, which harbors the bla NDM-1 gene. The total length of the assembled genome is 4,103,824 bp, and 3,896 coding sequences (CDSs) were predicted within the genome. A previously unreported bla NDM-1 -bearing plasmid was identified in this strain.

International audience; Plague is a pandemic human invasive disease caused by the bacterial agent Yersinia pestis. We here report a comparison of 17 whole genomes of Y. pestis isolates from global sources. We also screened a global collection of 286 Y. pestis isolates for 933 SNPs using Sequenom MassArray SNP typing. We conducted phylogenetic analyses on this sequence variation dataset, assigned isolates to populations based on maximum parsimony and, from these results, made inferences regarding historical transmission routes. Our phylogenetic analysis suggests that Y. pestis evolved in or near China and spread through multiple radiations to Europe, South America, Africa and Southeast Asia, leading to country-specific lineages that can be traced by lineage-specific SNPs. All 626 current isolates from the United States reflect one radiation, and 82 isolates from Madagascar represent a second radiation. Subsequent local microevolution of Y. pestis is marked by sequential, geographically specific SNPs.

To gain insights into the evolutionary origin, emergence, and pathogenicity of the etiologic agent of plague, we have sequenced the genomes of four Yersinia pestis strains isolated from the zoonotic rodent reservoir in foci of endemic plague in China. These resources enable in-depth studies of Y. pestis sequence variations and detailed whole-genome comparisons of very closely related genomes from the supposed site of the origin and the emergence of global pandemics of plague.

International audience; BACKGROUND: The species Yersinia pestis is commonly divided into three classical biovars, Antiqua, Medievalis, and Orientalis, belonging to subspecies pestis pathogenic for human and the (atypical) non-human pathogenic biovar Microtus (alias Pestoides) including several non-pestis subspecies. Recent progress in molecular typing methods enables large-scale investigations in the population structure of this species. It is now possible to test hypotheses about its evolution which were proposed decades ago. For instance the three classical biovars of different geographical distributions were suggested to originate from Central Asia. Most investigations so far have focused on the typical pestis subspecies representatives found outside of China, whereas the understanding of the emergence of this human pathogen requires the investigation of strains belonging to subspecies pestis from China and to the Microtus biovar. METHODOLOGY/PRINCIPAL FINDINGS: Multi-locus VNTR analysis (MLVA) with 25 loci was performed on a collection of Y. pestis isolates originating from the majority of the known foci worldwide and including typical rhamnose-negative subspecies pestis as well as rhamnose-positive subspecies pestis and biovar Microtus. More than 500 isolates from China, the Former Soviet Union (FSU), Mongolia and a number of other foci around the world were characterized and resolved into 350 different genotypes. The data revealed very close relationships existing between some isolates from widely separated foci as well as very high diversity which can conversely be observed between nearby foci. CONCLUSIONS/SIGNIFICANCE: The results obtained are in full agreement with the view that the Y. pestis subsp. pestis pathogenic for humans emerged in the Central Asia region between China, Kazakhstan, Russia and Mongolia, only three clones of which spread out of Central Asia. The relationships among the strains in China, Central Asia and the rest of the world based on the MLVA25 assay provide an unprecedented view on the expansion and microevolution of Y. pestis.

Background DFR (different region) analysis has been developed for typing Yesinia pestis in our previous study, and in this study, we extended this method by using 23 DFRs to investigate 909 Chinese Y. pestis strains for validating DFR-based genotyping method and better understanding adaptive microevolution of Y. pestis. Methodology/Principal Findings On the basis of PCR and Bionumerics data analysis, 909 Y. pestis strains were genotyped into 32 genomovars according to their DFR profiles. New terms, Major genomovar and Minor genomovar, were coined for illustrating evolutionary relationship between Y. pestis strains from different plague foci and different hosts. In silico DFR profiling of the completed or draft genomes shed lights on the evolutionary scenario of Y. pestis from Y. pseudotuberculosis. Notably, several sequenced Y. pestis strains share the same DFR profiles with Chinese strains, providing data for revealing the global plague foci expansion. Conclusions/significance Distribution of Y. pestis genomovars is plague focus-specific. Microevolution of biovar Orientalis was deduced according to DFR profiles. DFR analysis turns to be an efficient and inexpensive method to portrait the genome plasticity of Y. pestis based on horizontal gene transfer (HGT). DFR analysis can also be used as a tool in comparative and evolutionary genomic research for other bacteria with similar genome plasticity.

The ferric uptake regulator (Fur) is a predominant bacterial regulator controlling the iron assimilation functions in response to iron availability. Our previous microarray analysis on Yersinia pestis defined the iron-Fur modulon. In the present work, we reannotated the iron assimilation genes in Y. pestis , and the resulting genes in complementation with those disclosed by microarray constituted a total of 34 genome loci (putative operons) that represent the potential iron-responsive targets of Fur. The subsequent real-time reverse transcription-PCR (RT-PCR) in conjunction with the primer extension analysis showed that 32 of them were regulated by Fur in response to iron starvation. A previously predicted Fur box sequence was then used to search against the promoter regions of the 34 operons; the homologue of the above box could be predicted in each promoter tested. The subsequent electrophoretic mobility shift assay (EMSA) demonstrated that a purified His 6 tag-fused Fur protein was able to bind in vitro to each of these promoter regions. Therefore, Fur is a global regulator, both an activator and a repressor, and directly controls not only almost all of the iron assimilation functions but also a variety of genes involved in various non-iron functions for governing a complex regulatory cascade in Y. pestis . In addition, real-time RT-PCR, primer extension, EMSA, and DNase I footprinting assay were used to elucidate the Fur regulation of the ybt locus encoding a virulence-required iron uptake system. By combining the published data on the YbtA regulation of ybt , we constructed a concise Fur/YbtA regulatory network with a map of the Fur-promoter DNA interactions within the ybt locus. The data presented here give us an overview of the iron-responsive Fur regulon in Y. pestis .

BackgroundYersinia pestis, the pathogen of plague, has greatly influenced human history on a global scale. Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR), an element participating in immunity against phages' invasion, is composed of short repeated sequences separated by unique spacers and provides the basis of the spoligotyping technology. In the present research, three CRISPR loci were analyzed in 125 strains of Y. pestis from 26 natural plague foci of China, the former Soviet Union and Mongolia were analyzed, for validating CRISPR-based genotyping method and better understanding adaptive microevolution of Y. pestis.Methodology/principal findingsUsing PCR amplification, sequencing and online data processing, a high degree of genetic diversity was revealed in all three CRISPR elements. The distribution of spacers and their arrays in Y. pestis strains is strongly region and focus-specific, allowing the construction of a hypothetic evolutionary model of Y. pestis. This model suggests transmission route of microtus strains that encircled Takla Makan Desert and ZhunGer Basin. Starting from Tadjikistan, one branch passed through the Kunlun Mountains, and moved to the Qinghai-Tibet Plateau. Another branch went north via the Pamirs Plateau, the Tianshan Mountains, the Altai Mountains and the Inner Mongolian Plateau. Other Y. pestis lineages might be originated from certain areas along those routes.Conclusions/significanceCRISPR can provide important information for genotyping and evolutionary research of bacteria, which will help to trace the source of outbreaks. The resulting data will make possible the development of very low cost and high-resolution assays for the systematic typing of any new isolate.

BackgroundEnvironmental modulation of gene expression inYersinia pestisis critical for its life style and pathogenesis. Using cDNA microarray technology, we have analyzed the global gene expression of this deadly pathogen when grown under different stress conditionsin vitro.ResultsTo provide us with a comprehensive view of environmental modulation of global gene expression inY. pestis, we have analyzed the gene expression profiles of 25 different stress conditions. Almost all known virulence genes ofY. pestiswere differentially regulated under multiple environmental perturbations. Clustering enabled us to functionally classify co-expressed genes, including some uncharacterized genes. Collections of operons were predicted from the microarray data, and some of these were confirmed by reverse-transcription polymerase chain reaction (RT-PCR). Several regulatory DNA motifs, probably recognized by the regulatory protein Fur, PurR, or Fnr, were predicted from the clustered genes, and a Fur binding site in the corresponding promoter regions was verified by electrophoretic mobility shift assay (EMSA).ConclusionThe comparative transcriptomics analysis we present here not only benefits our understanding of the molecular determinants of pathogenesis and cellular regulatory circuits inY. pestis, it also serves as a basis for integrating increasing volumes of microarray data using existing methods.

The precise nature of the pathogen having caused early plague pandemics is uncertain. Although Yersinia pestis is a likely candidate for all three plague pandemics, the very rare direct evidence that can be deduced from ancient DNA (aDNA) analysis is controversial. Moreover, which of the three biovars, Antiqua, Medievalis or Orientalis, was associated with these pandemics is still debated. There is a need for phylogenetic analysis performed on Y. pestis strains isolated from countries from which plague probably arose and is still endemic. In addition there exist technical difficulties inherent to aDNA investigations and a lack of appropriate genetic targets. The recently described CRISPRs (clustered regularly interspaced short palindromic repeats) may represent such a target. CRISPR loci consist of a succession of highly conserved regions separated by specific "spacers" usually of viral origin. To be of use, data describing the mechanisms of evolution and diversity of CRISPRs in Y. pestis, its closest neighbors, and other species which might contaminate ancient DNA, are necessary. The investigation of closely related Y. pestis isolates has revealed recent mutation events in which elements constituting CRISPRs were acquired or lost, providing essential insight on their evolution. Rules deduced represent the basis for subsequent interpretation. In the present study, the CRISPR loci from representative Y. pestis and Yersinia pseudotuberculosis strains were investigated by PCR amplification and sequence analysis. The investigation of this wider panel of strains, including other subspecies or ecotypes within Y. pestis and also Y. pseudotuberculosis strains provides a database of the existing CRISPR spacers and helps predict the expected CRISPR structure of the Y. pestis ancestor. This knowledge will open the way to the development of a spoligotyping assay, in which spacers can be amplified even from highly degraded DNA samples. The data obtained show that CRISPR analysis can provide a very powerful typing tool, adapted to the systematic, large-scale genotyping of Y. pestis isolates, and the creation of international typing databases. In addition, CRISPRs do constitute a very promising new tool and genetic target to investigate ancient DNA. The corresponding genetic targets are small (

Plague is a deadly disease caused by Yersinia pestis. Human plague can be effectively controlled by timely antibiotic administration, chloramphenicol being a drug of choice. In this study, a DNA microarray was used to investigate the gene expression profile of Y. pestis in response to chloramphenicol. Seven hundred and fifty-five genes were differentially expressed on chloramphenicol treatment: 364 genes were up-regulated and 391 were down-regulated. In addition to a large number of genes encoding unknown or unassigned functions, genes encoding the components of the translation apparatus, cell envelope and transport/binding functions were strongly represented amongst the induced genes. Genes encoding proteins involved in energy metabolism and synthesis and modification of macromolecules were strongly represented amongst the down-regulated genes. A number of heat-shock genes were also repressed. These global transcriptional changes provide an insight into the mechanisms of action of chloramphenicol against Y. pestis.

Udgivelsesdato: 2006-Aug In order to further figure out the genetic differences between Yersinia pestis and Yersinia pseudotuberculosis, and to provide novel insights into the evolution of Y. pestis, we compared the genomes of Y. pseudotuberculosis serogroup I strain ATCC29833 and Y. pestis Antiqua strain 49006 using a combination of suppression subtractive hybridization (SSH) and comparative genomic hybridization with DNAs from a diverse panel of Y. pestis and Y. pseudotuberculosis strains. SSH followed by BLAST analysis revealed 112 SSH fragments specific to strain ATCC29833, compared to the genomic sequence data of Y. pestis strains CO92, KIM and 91001. We identified 17 SSH fragments that appeared to be newly determined genetic contents of Y. pseudotuberculosis. The combination of SSH and microarray analysis showed that the parallel loss of genes contributed greatly not only to the significant genomic divergence between Y. pestis and Y. pseudotuberculosis but also to the intra-species microevolution of both of species. The results confirmed our earlier hypothesis that Y. pestis Antiqua isolates from the natural plague focus B in China represented the most ancestral strains in China, hence phylogenetically the closest isolates to Y. pseudotuberculosis.

Yersinia pestis, the etiologic agent of plague, must acclimatize itself to temperature shifts between the temperature (26 C) for flea blockage and the body temperature (37 C) of warm-blooded hosts during its life cycle. Here a whole-genome DNA microarray was used to investigate transcriptional regulation upon the upshift of growth temperature from 26 to 37 C in a chemically defined medium. Four hundred and one genes were regulated differentially under the two temperatures. About 39% of these genes were up-regulated at 37 C, whereas 61% were down-regulated. Temperature-induced changes occurred at the level of transcription of genes encoding proven or predicted virulence factors, regulators, metabolism-associated proteins, prophages, and hypothetical proteins. Strikingly, many gene clusters displayed a co-transcription pattern in response to temperature upshift. Our data provided a genome-wide profile of gene transcription induced by temperature shift and should shed light on the pathogenicity and host-microbe interaction of this deadly pathogen.

Plague, caused by Y ersinia pestis, is one of the most dangerous diseases that impressed a horror onto human consciousness that persists to this day. Cases of plague can be normally controlled by timely antibiotic administration. Streptomycin is the first-line antibiotic for plague treatment. In this study, a DNA microarray was used to investigate the changes in the gene expression profile of Y. pestis upon exposure to streptomycin. A total of 345 genes were identified to be differentially regulated, 144 of which were up-regulated, and 201 down-regulated. Streptomycin-induced transcriptional changes occurred in genes responsible for heat shock response, drug/analogue sensitivity, biosynthesis of the branched-chain amino acids, chemotaxis and mobility and broad regulatory functions. A wide set of genes involved in energy metabolism, biosynthesis of small macromolecules, synthesis and modification of macromoclecules and degradation of small and macro molecules were among those down-regulated. The results reveal general changes in gene expression that are consistent with known mechanisms of action of streptomycin and many new genes that are likely to play important roles in the response to streptomycin treatment, providing useful candidates for investigating the specific mechanisms of streptomycin action.

DNA microarray was used as a tool to define the heat- and cold-shock stimulons in Yersinia pestis. Heat shock dramatically enhanced the transcription of genes encoding major heat-shock proteins (MHSPs) that are important for cell survival against the heat. Many other genes were also greatly up-regulated, but their roles in heat-shock response need to be elucidated. Meanwhile, heat shock retarded most of the metabolic processes, i.e. RNA transcription, protein translation, aerobic respiration, energy metabolism, small molecule metabolism, peptidoglycan biosynthesis, sulfate uptake and cysteine biosynthesis. In response to cold shock, Y. pestis has evolved complex adaptive mechanisms by elevating the transcription of a specific set of genes whose protein products are designed to prevent or eliminate cold-induced DNA or RNA structuring, to remodel cell membrane components for maintenance of normal functions, to elevate the energy generation for ensuring ATP-dependent responses during cold adaptation and to synthesize or transport compatible solutes such as cryoprotectants, and at the same time, by repressing the mRNA level of certain genes whose protein products are not needed for bacterial growth at low temperatures, such as the MHSPs. These results provide a set of new candidate genes for hypothesis-based investigations of their roles in stress response, host adaptation and pathogenicity of this deadly pathogen.

Polymerase chain reaction (PCR) amplification of DNA-based unique markers, the signature sequences, is ideal for rapid detection and identification of pathogens. We described the discovery of twenty-eight signature genes of Yersinia pestis by DNA microarray-based comparative genome hybridization in conjunction with PCR validation. Three pairs of Y. pestis-specific primers designed from signature genes were demonstrated to have the expected specificity to this target bacterium, without cross-reaction with the closely related Y. pseudotuberculosis or a large collection of genomic DNAs from other organisms.

Yersinia pestis whole-genome DNA microarrays were developed to perform genomic comparison of a collection of live plague vaccines. By using the genomic DNA to probe the DNA microarrays, we detected dozens of deletions and amplifications of the genomic regions in the 19 vaccine strains analyzed. The revealed genomic differences within the vaccine strains of different origins provide us an unprecedented opportunity to understand the molecular background of the variability of the immunogenic and protective potency of plague live vaccine. The whole-genome DNA microarray also provides an ideal tool to perform the pre-evaluation of a vaccine strain for its high throughput to determine the genomic features essential or unallowable for the live vaccines.

Background Variable number of tandem repeats (VNTRs) that are widely distributed in the genome of Yersinia pestis proved to be useful markers for the genotyping and source-tracing of this notorious pathogen. In this study, we probed into the features of VNTRs in the Y. pestis genome and developed a simple hierarchical genotyping system based on optimized VNTR loci. Methodology/Principal Findings Capillary electrophoresis was used in this study for multi-locus VNTR analysis (MLVA) in 956 Y. pestis strains. The general features and genetic diversities of 88 VNTR loci in Y. pestis were analyzed with BioNumerics, and a “14+12” loci-based hierarchical genotyping system, which is compatible with single nucleotide polymorphism-based phylogenic analysis, was established. Conclusions/Significance Appropriate selection of target loci reduces the impact of homoplasies caused by the rapid mutation rates of VNTR loci. The optimized “14+12” loci are highly discriminative in genotyping and source-tracing Y. pestis for molecular epidemiological or microbial forensic investigations with less time and lower cost. An MLVA genotyping datasets of representative strains will improve future research on the source-tracing and microevolution of Y. pestis.