Dynamics of CRISPR loci in microevolutionary process of Yersinia pestis strains

CRISPR genotyping methods: Tracing the evolution from spoligotyping to machine learning

CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) systems provide prokaryotes with adaptive immunity defenses against foreign genetic invaders. The identification of CRISPR-Cas function is among the most impactful discoveries of recent decades that have shaped the development of genome editing in various organisms paving the way for a plethora of promising applications in biotechnology and health. Even before the discovery of CRISPR-Cas biological role, the particular structure of CRISPR loci has been explored for epidemiological genotyping of bacterial pathogens. CRISPR-Cas loci are arranged in CRISPR arrays of mostly identical direct repeats intercalated with invader-derived spacers and an operon of cas genes encoding the Cas protein components. Each small CRISPR RNA (crRNA) encoded within the CRISPR array constitutes a key functional unit of this RNA-based CRISPR-Cas defense system guiding the Cas effector proteins toward the foreign nucleic acids for their destruction. The information acquired from prior invader encounters and stored within CRISPR arrays turns out to be extremely valuable in tracing the microevolution and epidemiology of major bacterial pathogens. We review here the history of CRISPR-based typing strategies highlighting the first PCR-based methods that have set the stage for recent developments of high-throughput sequencing and machine learning-based approaches. A great amount of whole genome sequencing and metagenomic data accumulated in recent years opens up new avenues for combining experimental and computational approaches of high-resolution CRISPR-based typing.

Analysis of CRISPR/Cas Genetic Structure, Spacer Content and Molecular Epidemiology in Brazilian Acinetobacter baumannii Clinical Isolates

CRISPR/Cas is a molecular mechanism to prevent predatory viruses from invading bacteria via the insertion of small viral sequences (spacers) in its repetitive locus. The nature of spacer incorporation and the viral origins of spacers provide an overview of the genetic evolution of bacteria, their natural viral predators, and the mechanisms that prokaryotes may use to protect themselves, or to acquire mobile genetic elements such as plasmids. Here, we report on the CRISPR/Cas genetic structure, its spacer content, and strain epidemiology through MLST and CRISPR typing in Acinetobacter baumannii, an opportunistic pathogen intimately related to hospital infections and antimicrobial resistance. Results show distinct genetic characteristics, such as polymorphisms specific to ancestor direct repeats, a well-defined degenerate repeat, and a conserved leader sequence, as well as showing most spacers as targeting bacteriophages, and several self-targeting spacers, directed at prophages. There was a particular relationship between CRISPR/Cas and CC113 in the study of Brazilian isolates, and CRISPR-related typing techniques are interesting for subtyping strains with the same MLST profile. We want to emphasize the significance of descriptive genetic research on CRISPR loci, and we argue that spacer or CRISPR typing are helpful for small-scale investigations, preferably in conjunction with other molecular typing techniques such as MLST.

CRISPR genotyping as complementary tool for epidemiological surveillance of Erwinia amylovora outbreaks

Fire blight is a destructive plant disease caused by Erwinia amylovora affecting pome fruit trees, and responsible for large yield declines, long phytosanitary confinements, and high economic losses. In Portugal, the first major fire blight outbreaks occurred in 2010 and 2011, and although later considered eradicated, the emergence of other outbreaks in recent years stressed the need to characterize the E. amylovora populations associated with these outbreaks. In this regard, CRISPR genotyping, assessment of three virulence markers, and semi-quantitative virulence bioassays, were carried out to determine the genotype, and assess the virulence of thirty-six E. amylovora isolates associated with outbreaks occurring between 2010 and 2017 and affecting apple and pear orchards located in the country central-west, known as the main producing region of pome fruits in Portugal. The data gathered reveal that 35 E. amylovora isolates belong to one of the widely-distributed CRISPR genotypes (5-24-38 / D-a-α) regardless the host species, year and region. Ea 680 was the single isolate revealing a new CRISPR genotype due to a novel CR2 spacer located closer to the leader sequence and therefore thought to be recently acquired. Regarding pathogenicity, although dot-blot hybridization assays showed the presence of key virulence factors, namely hrpL (T3SS), hrpN (T3E) and amsG from the amylovoran biosynthesis operon in all E. amylovora isolates studied, pathogenicity bioassays on immature pear slices allowed to distinguish four virulence levels, with most of the isolates revealing an intermediate to severe virulence phenotype. Regardless the clonal population structure of the E. amylovora associated to the outbreaks occurring in Portugal between 2010 and 2017, the different virulence phenotypes, suggests that E. amylovora may have been introduced at different instances into the country. This is the first study regarding E. amylovora in Portugal, and it discloses a novel CRISPR genotype for this bacterium.

Polymorphism of Type I-F CRISPR/Cas system in Escherichia coli of phylogenetic group B2 and its application in genotyping

E. coli of phylogenetic group B2 is responsible for many extraintestinal infections, posing a great threat to health. The relatively polymorphic nature of CRISPR in phylogenetically related E. coli strains makes them potential markers for bacterial typing and evolutionary studies. In the current work, we investigated the occurrence and diversity of CRISPR/Cas system and explored its potential for genotyping. Type I-F CRISPR/Cas systems were found in 413 of 1190 strains of E. coli and exhibited the clustering within certain CCs and STs. And CRISPR spacer contents correlated well with MLST types. The divergence analysis of CRISPR showed stronger discriminatory power than MLST, and CRISPR polymorphism was instrumental for differentiating highly closely related strains. The timeline of spacer acquisition and deletion provided important information for inferring the evolution model between distinct serotypes. Identical spacer sequences were shared by strains with the same H-antigen type but not strains with the same O-antigen type. The homology between spacers and antibiotic-resistant plasmids demonstrated the role of Type I-F system in limiting the acquisition of antimicrobial resistance. Collectively, our data presents the dynamic nature of Type I-F CRISPR in E. coli of phylogenetic group B2 and provides new insights into the application of CRISPR-based typing in the species.

Analysis of direct repeats and spacers of CRISPR/Cas systems type I-F in Brazilian clinical strains of Pseudomonas aeruginosa

CRISPR/Cas is an adaptive immune system found in prokaryotes, with the main function of protecting these cells from invasion and possible death by mobile genetic elements. Pseudomonas aeruginosa is considered a model for type I-F CRISPR/Cas system studies. However, its CRISPR loci characteristics have not yet been thoroughly described, and its function has not yet been fully unraveled. The aims of this study were to find the frequency of the system in Brazilian clinical isolates; to identify the loci sequence, its spacer diversity and its origins; as well as to propose a unified spacer library to aid in future structural studies of the CRISPR loci of P. aeruginosa. We investigated types I-F and I-E gene markers to establish CRISPR/Cas typing, and observed two strains harboring both systems simultaneously, a very rare feature. Through amplification and sequencing of CRISPR loci related to type I-F system, we describe polymorphisms in DRs and 350 spacers, of which 97 are new. The spacers that were identified had their possible organisms or proteins of origin identified. Spacer arrays were grouped in five different CRISPR patterns and the plasticity was inferred by rearrangements in spacer arrays. Here, we perform the first detailed and focused description of CRISPR/Cas elements in Brazilian clinical strains of P. aeruginosa. Our findings reflect active and highly diverse CRISPR loci, and we suggest that CRISPR/Cas may also pose as a transcriptional regulatory mechanism. The structural and diversity features described here can provide insights into the function of CRISPR/Cas in this pathogen and help guide the development of new therapeutic strategies.

Independent Microevolution Mediated by Mobile Genetic Elements of Individual Clostridium difficile Isolates from Clade 4 Revealed by Whole-Genome Sequencing

Mobile genetic elements play a key role in the continuing evolution of Clostridium difficile, resulting in the emergence of new phenotypes for individual isolates. On the basis of whole-genome sequencing analysis, we comprehensively explored transposons, CRISPR, prophage, and genetic sites for drug resistance within clade 4 C. difficile isolates with different sequence types. Great diversity in MGEs and a high rate of multidrug resistance were found within this clade, including new transposons, Tn4453a/b with aac(6′) aph(2′′) instead of catD, and a relatively high rate of prophage-carried CRISPR arrays. These findings provide important new insights into the mechanism of genome remodeling within clade 4 and offer a new method for typing and tracing the origins of closely related isolates.

Horizontal gene transfer of mobile genetic elements (MGEs) accounts for the mosaic genome of Clostridium difficile, leading to acquisition of new phenotypes, including drug resistance and reconstruction of the genomes. MGEs were analyzed according to the whole-genome sequences of 37 C. difficile isolates with a variety of sequence types (STs) within clade 4 from China. Great diversity was found in each transposon even within isolates with the same ST. Two novel transposons were identified in isolates ZR9 and ZR18, of which approximately one third to half of the genes showed heterogenous origins compared with the usual intestinal bacterial genes. Most importantly, catD, known to be harbored by Tn4453a/b, was replaced by aac(6′) aph(2′′) in isolates 2, 7, and 28. This phenomenon illustrated the frequent occurrence of gene exchanges between C. difficile and other enterobacteria with individual heterogeneity. Numerous prophages and CRISPR arrays were identified in C. difficile isolates of clade 4. Approximately 20% of spacers were located in prophage-carried CRISPR arrays, providing a new method for typing and tracing the origins of closely related isolates, as well as in-depth studies of the mechanism underlying genome remodeling. The rates of drug resistance were obviously higher than those reported previously around the world, although all isolates retained high sensitivity to vancomycin and metronidazole. The increasing number of C. difficile isolates resistant to all antibiotics tested here suggests the ease with which resistance is acquired in vivo. This study gives insights into the genetic mechanism of microevolution within clade 4.

IMPORTANCE Mobile genetic elements play a key role in the continuing evolution of Clostridium difficile, resulting in the emergence of new phenotypes for individual isolates. On the basis of whole-genome sequencing analysis, we comprehensively explored transposons, CRISPR, prophage, and genetic sites for drug resistance within clade 4 C. difficile isolates with different sequence types. Great diversity in MGEs and a high rate of multidrug resistance were found within this clade, including new transposons, Tn4453a/b with aac(6′) aph(2′′) instead of catD, and a relatively high rate of prophage-carried CRISPR arrays. These findings provide important new insights into the mechanism of genome remodeling within clade 4 and offer a new method for typing and tracing the origins of closely related isolates.

A single introduction of Yersinia pestis to Brazil during the 3rd plague pandemic

Yersinia pestis was introduced to Brazil during the third plague pandemic and currently exists in several recognized foci. There is currently limited available phylogeographic data regarding Y. pestis in Brazil. We generated whole genome sequences for 411 Y. pestis strains from six Brazilian foci to investigate the phylogeography of Y. pestis in Brazil; these strains were isolated from 1966 to 1997. All 411 strains were assigned to a single monophyletic clade within the 1.ORI population, indicating a single Y. pestis introduction was responsible for the successful establishment of endemic foci in Brazil. There was a moderate level of genomic diversity but little population structure among the 411 Brazilian Y. pestis strains, consistent with a radial expansion wherein Y. pestis spread rapidly from the coast to the interior of Brazil and became ecologically established. Overall, there were no strong spatial or temporal patterns among the Brazilian strains. However, strains from the same focus tended to be more closely related and strains isolated from foci closer to the coast tended to fall in more basal positions in the whole genome phylogeny than strains from more interior foci. Overall, the patterns observed in Brazil are similar to other locations affected during the 3rd plague pandemic such as in North America and Madagascar.

The CRISPR-Cas system in Enterobacteriaceae

In nature, microorganisms are constantly exposed to multiple viral infections and thus have developed many strategies to survive phage attack and invasion by foreign DNA. One of such strategies is the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated proteins (Cas) bacterial immunological system. This defense mechanism is widespread in prokaryotes including several families such as Enterobacteriaceae. Much knowledge about the CRISPR-Cas system has been generated, including its biological functions, transcriptional regulation, distribution, utility as a molecular marker and as a tool for specific genome editing. This review focuses on these aspects and describes the state of the art of the CRISPR-Cas system in the Enterobacteriaceae bacterial family.

New Insights Into Functions and Possible Applications of Clostridium difficile CRISPR-Cas System

Over the last decades the enteric bacterium Clostridium difficile (novel name Clostridioides difficile) - has emerged as an important human nosocomial pathogen. It is a leading cause of hospital-acquired diarrhea and represents a major challenge for healthcare providers. Many aspects of C. difficile pathogenesis and its evolution remain poorly understood. Efficient defense systems against phages and other genetic elements could have contributed to the success of this enteropathogen in the phage-rich gut communities. Recent studies demonstrated the presence of an active CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) subtype I-B system in C. difficile. In this mini-review, we will discuss the recent advances in characterization of original features of the C. difficile CRISPR-Cas system in laboratory and clinical strains, as well as interesting perspectives for our understanding of this defense system function and regulation in this important enteropathogen. This knowledge will pave the way for the development of promising biotechnological and therapeutic tools in the future. Possible applications for the C. difficile strain monitoring and genotyping, as well as for CRISPR-based genome editing and antimicrobials are also discussed.

Reemergence of human plague in Yunnan, China in 2016

The third plague pandemic originated from Yunnan Province, China in the middle of the 19th century. The last human plague epidemic in Yunnan occurred from 1986-2005. On June 6, 2016, a case of human plague was reported in the Xishuangbanna Prefecture, Yunnan. The patient suffered from primary septicemic plague after exposure to a dead house rat (Rattus flavipectus), which has been identified as the main plague reservoir in the local epizootic area. Moreover, a retrospective investigation identified another bubonic plague case in this area. Based on these data, human plague reemerged after a silent period of ten years. In this study, three molecular typing methods, including a clustered regularly interspaced short palindromic repeats (CRISPR) analysis, different region analysis (DFR), and multiple-locus variable number of tandem repeats analysis (MLVA), were used to illustrate the molecular characteristics of Yersinia pestis (Y. pestis) strains isolated in Yunnan. The DFR profiles of the strains isolated in Yunnan in 2016 were the same as the strains that had previously been isolated in this Rattus flavipectus plague focus. The c3 spacer present in the previously isolated strains was absent in the spacer arrays of the Ypc CRISPR loci of the strains isolated in 2016. The MLVA analysis using MLVA (14+12) showed that the strains isolated from the human plague case and host animal plague infection in 2016 in Yunnan displayed different molecular patterns than the strains that had previously been isolated from Yunnan and adjacent provinces.

Genome rearrangements and phylogeny reconstruction in Yersinia pestis

Genome rearrangements have played an important role in the evolution of Yersinia pestis from its progenitor Yersinia pseudotuberculosis. Traditional phylogenetic trees for Y. pestis based on sequence comparison have short internal branches and low bootstrap supports as only a small number of nucleotide substitutions have occurred. On the other hand, even a small number of genome rearrangements may resolve topological ambiguities in a phylogenetic tree. We reconstructed phylogenetic trees based on genome rearrangements using several popular approaches such as Maximum likelihood for Gene Order and the Bayesian model of genome rearrangements by inversions. We also reconciled phylogenetic trees for each of the three CRISPR loci to obtain an integrated scenario of the CRISPR cassette evolution. Analysis of contradictions between the obtained evolutionary trees yielded numerous parallel inversions and gain/loss events. Our data indicate that an integrated analysis of sequence-based and inversion-based trees enhances the resolution of phylogenetic reconstruction. In contrast, reconstructions of strain relationships based on solely CRISPR loci may not be reliable, as the history is obscured by large deletions, obliterating the order of spacer gains. Similarly, numerous parallel gene losses preclude reconstruction of phylogeny based on gene content.

Spontaneous CRISPR loci generation in vivo by non-canonical spacer integration

The adaptation phase of CRISPR-Cas immunity depends on the precise integration of short segments of foreign DNA (spacers) into a specific genomic location within the CRISPR locus by the Cas1-Cas2 integration complex. Although off-target spacer integration outside of canonical CRISPR arrays has been described in vitro, no evidence of non-specific integration activity has been found in vivo. Here, we show that non-canonical off-target integrations can occur within bacterial chromosomes at locations that resemble the native CRISPR locus by characterizing hundreds of off-target integration locations within Escherichia coli. Considering whether such promiscuous Cas1-Cas2 activity could have an evolutionary role through the genesis of neo-CRISPR loci, we combed existing CRISPR databases and available genomes for evidence of off-target integration activity. This search uncovered several putative instances of naturally occurring off-target spacer integration events within the genomes of Yersinia pestis and Sulfolobus islandicus. These results are important in understanding alternative routes to CRISPR array genesis and evolution, as well as in the use of spacer acquisition in technological applications.

Spontaneous CRISPR loci generation in vivo by non-canonical spacer integration

The adaptation phase of CRISPR-Cas immunity depends on the precise integration of short segments of foreign DNA (spacers) into a specific genomic location within the CRISPR locus by the Cas1-Cas2 integration complex. Although off-target spacer integration outside of canonical CRISPR arrays has been described in vitro, no evidence of non-specific integration activity has been found in vivo. Here, we show that non-canonical off-target integrations can occur within bacterial chromosomes at locations that resemble the native CRISPR locus by characterizing hundreds of off-target integration locations within Escherichia coli. Considering whether such promiscuous Cas1-Cas2 activity could have an evolutionary role through the genesis of neo-CRISPR loci, we combed existing CRISPR databases and available genomes for evidence of off-target integration activity. This search uncovered several putative instances of naturally occurring off-target spacer integration events within the genomes of Yersinia pestis and Sulfolobus islandicus. These results are important in understanding alternative routes to CRISPR array genesis and evolution, as well as in the use of spacer acquisition in technological applications.

Diversity and microevolution of CRISPR loci in Helicobacter cinaedi

Helicobacter cinaedi is associated with nosocomial infections. The CRISPR-Cas system provides adaptive immunity against foreign genetic elements. We investigated the CRISPR-Cas system in H. cinaedi to assess the potential of the CRISPR-based microevolution of H. cinaedi strains. A genotyping method based on CRISPR spacer organization was carried out using 42 H. cinaedi strains. The results of sequence analysis showed that the H. cinaedi strains used in this study had two CRISPR loci (CRISPR1 and CRISPR2). The lengths of the consensus direct repeat sequences in CRISPR1 and CRISPR2 were both 36 bp-long, and 224 spacers were found in the 42 H. cinaedi strains. Analysis of the organization and sequence similarity of the spacers of the H. cinaedi strains showed that CRISPR arrays could be divided into 7 different genotypes. Each genotype had a different ancestral spacer, and spacer acquisition/deletion events occurred while isolates were spreading. Spacer polymorphisms of conserved arrays across the strains were instrumental for differentiating closely-related strains collected from the same hospital. MLST had little variability, while the CRISPR sequences showed remarkable diversity. Our data revealed the structural features of H. cinaedi CRISPR loci for the first time. CRISPR sequences constitute a valuable basis for genotyping, provide insights into the divergence and relatedness between closely-related strains, and reflect the microevolutionary process of H. cinaedi.

Characterization and Exploitation of CRISPR Loci in Bifidobacterium longum

Diverse CRISPR-Cas systems provide adaptive immunity in many bacteria and most archaea, via a DNA-encoded, RNA-mediated, nucleic-acid targeting mechanism. Over time, CRISPR loci expand via iterative uptake of invasive DNA sequences into the CRISPR array during the adaptation process. These genetic vaccination cards thus provide insights into the exposure of strains to phages and plasmids in space and time, revealing the historical predatory exposure of a strain. These genetic loci thus constitute a unique basis for genotyping of strains, with potential of resolution at the strain-level. Here, we investigate the occurrence and diversity of CRISPR-Cas systems in the genomes of various Bifidobacterium longum strains across three sub-species. Specifically, we analyzed the genomic content of 66 genomes belonging to B. longum subsp. longum, B. longum subsp. infantis and B. longum subsp. suis, and identified 25 strains that carry 29 total CRISPR-Cas systems. We identify various Type I and Type II CRISPR-Cas systems that are widespread in this species, notably I-C, I-E, and II-C. Noteworthy, Type I-C systems showed extended CRISPR arrays, with extensive spacer diversity. We show how these hypervariable loci can be used to gain insights into strain origin, evolution and phylogeny, and can provide discriminatory sequences to distinguish even clonal isolates. By investigating CRISPR spacer sequences, we reveal their origin and implicate phages and prophages as drivers of CRISPR immunity expansion in this species, with redundant targeting of select prophages. Analysis of CRISPR spacer origin also revealed novel PAM sequences. Our results suggest that CRISPR-Cas immune systems are instrumental in mounting diversified viral resistance in B. longum, and show that these sequences are useful for typing across three subspecies.

Whole genome analysis of Yersinia ruckeri isolated over 27 years in Australia and New Zealand reveals geographical endemism over multiple lineages and recent evolution under host selection

Yersinia ruckeri is a salmonid pathogen with widespread distribution in cool-temperate waters including Australia and New Zealand, two isolated environments with recently developed salmonid farming industries. Phylogenetic comparison of 58 isolates from Australia, New Zealand, USA, Chile, Finland and China based on non-recombinant core genome SNPs revealed multiple deep-branching lineages, with a most recent common ancestor estimated at 18 500 years BP (12 355–24 757 95% HPD) and evidence of Australasian endemism. Evolution within the Tasmanian Atlantic salmon serotype O1b lineage has been slow, with 63 SNPs describing the variance over 27 years. Isolates from the prevailing lineage are poorly/non-motile compared to a lineage pre-vaccination, introduced in 1997, which is highly motile but has not been isolated since from epizootics. A non-motile phenotype has arisen independently in Tasmania compared to Europe and USA through a frameshift in fliI, encoding the ATPase of the flagella cluster. We report for the first time lipopolysaccharide O-antigen serotype O2 isolates in Tasmania. This phenotype results from deletion of the O-antigen cluster and consequent loss of high-molecular-weight O-antigen. This phenomenon has occurred independently on three occasions on three continents (Australasia, North America and Asia) as O2 isolates from the USA, China and Tasmania share the O-antigen deletion but occupy distant lineages. Despite the European and North American origins of the Australasian salmonid stocks, the lineages of Y. ruckeri in Australia and New Zealand are distinct from those of the northern hemisphere, suggesting they are pre-existing ancient strains that have emerged and evolved with the introduction of susceptible hosts following European colonization.

CRISPR Diversity and Microevolution in Clostridium difficile

Virulent strains of Clostridium difficile have become a global health problem associated with morbidity and mortality. Traditional typing methods do not provide ideal resolution to track outbreak strains, ascertain genetic diversity between isolates, or monitor the phylogeny of this species on a global basis. Here, we investigate the occurrence and diversity of clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated genes (cas) in C. difficile to assess the potential of CRISPR-based phylogeny and high-resolution genotyping. A single Type-IB CRISPR-Cas system was identified in 217 analyzed genomes with cas gene clusters present at conserved chromosomal locations, suggesting vertical evolution of the system, assessing a total of 1,865 CRISPR arrays. The CRISPR arrays, markedly enriched (8.5 arrays/genome) compared with other species, occur both at conserved and variable locations across strains, and thus provide a basis for typing based on locus occurrence and spacer polymorphism. Clustering of strains by array composition correlated with sequence type (ST) analysis. Spacer content and polymorphism within conserved CRISPR arrays revealed phylogenetic relationship across clades and within ST. Spacer polymorphisms of conserved arrays were instrumental for differentiating closely related strains, e.g., ST1/RT027/B1 strains and pathogenicity locus encoding ST3/RT001 strains. CRISPR spacers showed sequence similarity to phage sequences, which is consistent with the native role of CRISPR-Cas as adaptive immune systems in bacteria. Overall, CRISPR-Cas sequences constitute a valuable basis for genotyping of C. difficile isolates, provide insights into the micro-evolutionary events that occur between closely related strains, and reflect the evolutionary trajectory of these genomes.

Exploring Diversity and Origins of Streptomycin-Resistant Erwinia amylovora Isolates in New York Through CRISPR Spacer Arrays

Streptomycin is the most effective and widely used chemical control in the eastern United States for blossom blight of apple caused by Erwinia amylovora; however, resistance to this antibiotic has been a concern in New York since 2002. From 2011 to 2014, statewide collections of E. amylovora were conducted resulting in the isolation of streptomycin-resistant (SmR) E. amylovora from several commercial orchards. Further genetic analysis of isolates was necessary to understand the origins and the diversity of these bacteria. Clustered regularly interspaced short palindromic repeat (CRISPR) spacer sequencing was employed to explore the diversity and possible origins of New York SmR E. amylovora isolates. The spacer array CR1, CR2, and CR3 regions of 27 SmR E. amylovora isolates and 76 streptomycin-sensitive (SmS) E. amylovora isolates were amplified and subsequently sequenced, revealing 19 distinct CRISPR spacer profiles for New York isolates. The majority of SmR E. amylovora isolates had the same CRISPR profile as SmR E. amylovora isolates discovered in 2002. This may infer that eradication efforts in 2002 failed and the bacterial populations continued to spread throughout the state. Several CRISPR profiles for SmR E. amylovora were identical to SmS E. amylovora collected from the same orchards, leading to the hypothesis that resistance may be developing within New York. Profiles not unique to New York were identical to many isolates from the Midwestern, eastern, and western United States, implying that streptomycin resistance may be due to the introduction of SmR E. amylovora from other regions of the United States. The increased understanding as to how SmR E. amylovora isolates are introduced, evolve, or have become established afforded by CRISPR profiling has been useful for disease management and restricting the movement of streptomycin resistance in New York.

Building International Genomics Collaboration for Global Health Security

Genome science and technologies are transforming life sciences globally in many ways and becoming a highly desirable area for international collaboration to strengthen global health. The Genome Science Program at the Los Alamos National Laboratory is leveraging a long history of expertise in genomics research to assist multiple partner nations in advancing their genomics and bioinformatics capabilities. The capability development objectives focus on providing a molecular genomics-based scientific approach for pathogen detection, characterization, and biosurveillance applications. The general approaches include introduction of basic principles in genomics technologies, training on laboratory methodologies and bioinformatic analysis of resulting data, procurement, and installation of next-generation sequencing instruments, establishing bioinformatics software capabilities, and exploring collaborative applications of the genomics capabilities in public health. Genome centers have been established with public health and research institutions in the Republic of Georgia, Kingdom of Jordan, Uganda, and Gabon; broader collaborations in genomics applications have also been developed with research institutions in many other countries.

Viability of Yersinia pestis subcultures in agar stabs

Since its identification as the causative agent of plague in 1894, thousands of Yersinia pestis strains have been isolated and stored. Here, we report the ability of Y. pestis to survive up to 47 years in agar stabs, in rubber-stoppered tubes, under refrigeration (+4 to +10°C), although overall subculture recovery rates were poor and inversely related to the length of time stored. Genetic characterization of virulence gene presence among these subcultures was suggestive of significant variation in the genomic stability of Y. pestis subcultures stored under these conditions. Specifically, we found variation in the presence of plasmid and chromosomal virulence markers (genes pla, lcrV, caf1 and irp2) among multiple subcultures of Y. pestis strains in the 'Collection of Yersinia pestis' (Fiocruz-CYP) maintained by the SRP of FIOCRUZ-PE in Brazil. This variation, together with all of the inherent temporal, geographic and other genetic variation represented by all of the recoverable strains in this historical collection was preserved in new frozen culture stocks stored at -70°C as a result of this study. These frozen culture stocks represent a valuable resource for future comparative studies of Y. pestis.

SIGNIFICANCE AND IMPACT OF THE STUDY

We report the ability of Yersinia pestis to survive up to 47 years in agar stabs, in rubber-stoppered tubes, under refrigeration (+4 to +10°C), although overall subculture recovery rates were poor and inversely related to the length of time stored. Genetic characterization of virulence gene presence among these subcultures was suggestive of significant variation in the genomic stability of Y. pestis subcultures stored under these conditions. This variation, together with all of the inherent temporal, geographic and other genetic variation represented by all of the recoverable strains in the historical 'Collection of Yersinia pestis' (Fiocruz-CYP) maintained by the SRP of FIOCRUZ-PE in Brazil was preserved in new frozen culture stocks stored at -70°C as a result of this study. These frozen culture stocks represent a valuable resource for future comparative studies of Y. pestis.

A review of methods for subtyping Yersinia pestis: From phenotypes to whole genome sequencing

Numerous subtyping methods have been applied to Yersinia pestis with varying success. Here, we review the various subtyping methods that have been applied to Y. pestis and their capacity for answering questions regarding the population genetics, phylogeography, and molecular epidemiology of this important human pathogen. Methods are evaluated in terms of expense, difficulty, transferability among laboratories, discriminatory power, usefulness for different study questions, and current applicability in light of the advent of whole genome sequencing.