Identification and characterization of variable-number tandem repeats in the Yersinia pestis genome

Review of genotyping methods for Yersinia pestis in Madagascar

BACKGROUND

Plague, a zoonotic disease caused by Yersinia pestis, was responsible for 3 historical human pandemics that killed millions of people. It remains endemic in rodent populations in Africa, Asia, North America, and South America but human plague is rare in most of these locations. However, human plague is still highly prevalent in Madagascar, which typically records a significant part of all annual global cases. This has afforded an opportunity to study contemporary human plague in detail using various typing methods for Y. pestis.

AIM

This review aims to summarize the methods that have been used to type Y. pestis in Madagascar along with the major discoveries that have been made using these approaches.

METHODS

Pubmed and Google Scholar were used to search for the keywords: "typing Yersinia pestis Madagascar," "evolution Yersinia pestis Madagascar," and "diversity Yersinia pestis Madagascar." Eleven publications were relevant to our topic and further information was retrieved from references cited in those publications.

RESULTS

The history of Y. pestis typing in Madagascar can be divided in 2 periods: the pre-genomics and genomics eras. During the pre-genomics era, ribotyping, direct observation of plasmid content and plasmid restriction fragment length polymorphisms (RFLP) were employed but only revealed a limited amount of diversity among Malagasy Y. pestis strains. Extensive diversity only started to be revealed in the genomics era with the use of clustered regularly interspaced palindromic repeats (CRISPR), multiple-locus variable number tandem repeats (VNTR) analysis (MLVA), and single-nucleotide polymorphisms (SNPs) discovered from whole genome sequences. These higher-resolution genotyping methods have made it possible to highlight the distribution and persistence of genotypes in the different plague foci of Madagascar (Mahajanga and the Central and Northern Highlands) by genotyping strains from the same locations across years, to detect transfers between foci, to date the emergence of genotypes, and even to document the transmission of antimicrobial resistant (AMR) strains during a pneumonic plague outbreak. Despite these discoveries, there still remain topics that deserve to be explored, such as the contribution of horizontal gene transfer to the evolution of Malagasy Y. pestis strains and the evolutionary history of Y. pestis in Madagascar.

CONCLUSIONS

Genotyping of Y. pestis has yielded important insights on plague in Madagascar, particularly since the advent of whole-genome sequencing (WGS). These include a better understanding of plague persistence in the environment, antimicrobial AMR and multi-drug resistance in Y. pestis, and the person-to-person spread of pneumonic plague. Considering that human plague is still a significant public health threat in Madagascar, these insights can be useful for controlling and preventing human plague in Madagascar and elsewhere, and also are relevant for understanding the historical pandemics and the possible use of Y. pestis as a biological weapon.

<i>Yersinia pestis</i> ssp. <i>pestis</i> Spatial MLVA25 Genotypic Structure in the Transboundary Saylyugem Natural Plague Focus

Advanced molecular-genetic methods for the diagnosis and typing of Yersinia pestis ssp. pestis in the field and clinical material are used for epidemiological surveillance of plague in the Saylyugem natural focus. The aim of the work was to study the spatial genotypic structure of Y. pestis ssp. pestis in the transboundary Saylyugem natural plague focus using MLVA25 typing. Materials and methods. The MLVA25 typing of 160 strains of Y. pestis ssp. Pestis isolated in the Saylyugem natural plague focus in 2012–2021 was carried out. Phylogenetic tree construction was performed with the help of UPGMA and MST methods. Results and discussion. The Y. pestis ssp. pestis strains isolated from the Saylyugem natural plague focus were differentiated into 15 MLVA types by the 25 VNTR loci cluster analysis. The studied strains form a homogeneous complex of MLVA25 types without marked geographical distribution across seven spatial groups. The analysis of the frequency of occurrence of the tandem repeats number for three variable loci of Y. pestis ssp. pestis strains shows the significant differences between the samples from the Mongolian and Russian parts of the Saylyugem natural plague focus. The most pronounced differences in spatial genotypic structure are traced through the yp4280ms62 locus.

An Asymmetric Alignment Algorithm for Estimating Ancestor-Descendant Edit Distance for Tandem Repeats

Tandem repeats are repetitive structures present in some DNA sequences, consisting of many repeated copies of a single motif. They can serve as important markers for phylogenetic and population genetic studies, due to the high polymorphism in the number of motif copies as well as variations in the motif. The first step in using tandem repeats for phylogenetic studies is to estimate the evolutionary distance between a pair D₁ and D₂ of tandem repeat sequences with homologous motifs. This problem can be broken into two sub-problems: 1) Construct the most recent common ancestor of the sequences. 2) Calculate the evolutionary distance between each sequence and the hypothesised common ancestor. We present an algorithm that estimates the solution to the second problem. This takes the form of an asymmetric alignment algorithm to estimate the evolutionary distance between two tandem repeat sequences A and D, where D is assumed to have descended from A, under a model that allows block duplication, deletion, and variant substitution. The algorithm is asymmetric in the sense that the two input sequences A and D play different roles in the calculations, reflecting the assumption that D descends from A. Our model assumes static motif boundaries, meaning that motif duplication and deletion events must respect the motif boundaries. The algorithm may also be applied without modification to more complex repetitive structures with two or more motifs, such as nested tandem repeats.

Small Insertions and Deletions Drive Genomic Plasticity during Adaptive Evolution of Yersinia pestis

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 sspA, rpoS, rnd, and YPO0624, having enriched mutations in Brazilian isolates, which might be advantageous for Y. pestis to cope with fluctuating environments. The diversified selection signals observed here demonstrates that indels are important contributors to the adaptive evolution of Y. pestis. Meanwhile, we provide potential targets for further exploration, as some genes/pseudogenes with indels we focus on remain uncharacterized. IMPORTANCE Yersinia pestis, the causative agent of plague, is a highly pathogenic clone of Yersinia pseudotuberculosis. Previous genome-wide SNP analysis provided few adaptive signatures during its evolution. Here by investigating 365 public genomes of Y. pestis, we give a comprehensive overview of general features of genome-wide indels on the core genome and their roles in Y. pestis evolution. Detection of "indel pockets," with average complexity scores declining around indel positions, in both Mycobacterium tuberculosis and Y. pestis, gives us a clue that this phenomenon might appear in other bacterial genomes. Importantly, the identification of four different forms of selection signals in indels would improve our understanding on adaptive evolution of Y. pestis, and provide targets for further physiological mechanism researches of this pathogen. As evolutionary research based on genome-wide indels is still rare in bacteria, our study would be a helpful reference in deciphering the role of indels in other species.

Genome-wide analyses of the Nodulin-like gene family in bread wheat revealed its potential roles during arbuscular mycorrhizal symbiosis

Nodulin-like (NL) genes are involved in transporting of various substances and may play key roles during the establishment of symbiosis in legumes plants. However, basic biological information of NL genes in the wheat genome is still largely unknown. Here, we identified and characterized NL genes in wheat via integrating genomic information, collinearity analysis, co-expression network analysis (WGCNA) and transcriptome analysis. In addition, we analyzed the polymorphisms and the roles of NL genes during arbuscular mycorrhizal (AM) symbiosis using a large wheat panel consists of 259 wheat genotypes. We identified 181 NL genes in the wheat genome, which were classified into SWEET, Early Nodulin-Like (ENODL), Major Facilitator Superfamily-Nodulin (MFS), Vacuolar Iron Transporter (VIT) and Early nodulin 93 (ENOD93) subfamily. The expansion of NL genes was mainly driven by segmental duplication. The bHLH genes are potential unrecognized transcription factors regulating NL genes. Moreover, two NL genes were more sensitive than other NL genes to AM colonization. The polymorphisms of NL genes are mainly due to random drift, and the natural mutation of NL genes led to significant differences in the mycorrhizal dependence of wheat in phosphorus uptake. The results concluded that NL genes potentially play important roles during AM symbiosis with wheat.

Molecular epidemiology of clinically high-risk Pseudomonas aeruginosa strains: Practical overview

In recent years, numerous outbreaks of multidrug-resistant Pseudomonas aeruginosa have been reported across the world. Once an outbreak occurs, besides routinely testing isolates for susceptibility to antimicrobials, it is required to check their virulence genotypes and clonality profiles. Replacing pulsed-field gel electrophoresis DNA fingerprinting are faster, easier-to-use, and less expensive polymerase chain reaction (PCR)-based methods for characterizing hospital isolates. P. aeruginosa possesses a mosaic genome structure and a highly conserved core genome displaying low sequence diversity and a highly variable accessory genome that communicates with other Pseudomonas species via horizontal gene transfer. Multiple-locus variable-number tandem-repeat analysis and multilocus sequence typing methods allow for phylogenetic analysis of isolates by PCR amplification of target genes with the support of Internet-based services. The target genes located in the core genome regions usually contain low-frequency mutations, allowing the resulting phylogenetic trees to infer evolutionary processes. The multiplex PCR-based open reading frame typing (POT) method, integron PCR, and exoenzyme genotyping can determine a genotype by PCR amplifying a specific insertion gene in the accessory genome region using a single or a multiple primer set. Thus, analyzing P. aeruginosa isolates for their clonality, virulence factors, and resistance characteristics is achievable by combining the clonality evaluation of the core genome based on multiple-locus targeting methods with other methods that can identify specific virulence and antimicrobial genes. Software packages such as eBURST, R, and Dendroscope, which are powerful tools for phylogenetic analyses, enable researchers and clinicians to visualize clonality associations in clinical isolates.

Pulsed-field gel electrophoresis (PFGE): A review of the "gold standard" for bacteria typing and current alternatives

Pulsed-field gel electrophoresis (PFGE) is considered the "gold standard" for bacteria typing. The method involves enzyme restriction of bacteria DNA, separation of the restricted DNA bands using a pulsed-field electrophoresis chamber, followed by clonal assignment of bacteria based on PFGE banding patterns. Various PFGE protocols have been developed for typing different bacteria, leading it to be one of the most widely used methods for phylogenetic studies, food safety surveillance, infection control and outbreak investigations. On the other hand, as PFGE is lengthy and labourious, several PCR-based typing methods can be used as alternatives for research purposes. Recently, matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF MS) and whole genome sequencing (WGS) have also been proposed for bacteria typing. In fact, as WGS provides more information, such as antimicrobial resistance and virulence of the tested bacteria in comparison to PFGE, more and more laboratories are currently transitioning from PFGE to WGS for bacteria typing. Nevertheless, PFGE will remain an affordable and relevant technique for small laboratories and hospitals in years to come.

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.

Ten years of surveillance of the Yulong plague focus in China and the molecular typing and source tracing of the isolates

Plague, caused by Yersinia pestis, was classified as a reemerging infectious disease by the World Health Organization. The five human pneumonic plague cases in Yulong County in 2005 gave rise to the discovery of a Yulong plague focus in Yunnan province, China. Thereafter, continuous wild rodent plague (sylvatic plague) was identified as the main plague reservoir of this focus. In this study, the epizootics in Yulong focus were described, and three molecular typing methods, including the different region (DFR) analysis, clustered regularly interspaced short palindromic repeats (CRISPRs), and the multiple-locus variable number of tandem repeats (VNTR) analysis (MLVA) (14+12), were used for the molecular typing and source tracing of Y. pestis isolates in the Yulong plague focus. Simultaneously, several isolates from the vicinity of Yunnan were used as controls. The results showed that during the 10-year period from 2006 to 2016, an animal plague epidemic occurred in 6 of those years, and 5 villages underwent an animal plague epidemic within a 30-km² area of the Yulong plague focus. Searching for dead mice was the most effective monitoring method in this plague focus. No positive sample has been found in 6937 captured live rodents thus far, suggesting that the virulence of strains in the Yulong plague focus is stronger and the survival time of mice is shorter after infection. Strains from Lijiang, Sichuan and Tibet were of the same complex based on a typing analysis of DFR and CRISPR. The genetic relationship of Y. pestis illustrated by MLVA “14+12” demonstrates that Tibet and Sichuan strains evolved from the strains 1.IN2 (Qinghai, 1970 and Tibet, 1976), and Lijiang strains are closer to Batang strains (Batang County in Sichuan province, 2011, Himalaya marmot plague foci) in terms of genetic or phylogenic relationships. In conclusion, we have a deeper understanding of this new plague focus throughout this study, which provides a basis for effective prevention and control.

Plague is a type of zoonosis that is highly lethal to humans. The surveillance of animal hosts is critical for the prevention and control of plague. The Yulong plague focus is a newly discovered plague focus in China in recent years. The plague outbreak had attracted widespread attention because 5 people were infected in 2005, 2 of whom died. We have monitored the plague focus for a decade, and isolated strains and DNAs of Yersinia pestis were studied. The structure, origin and evolutionary trend of the Yulong plague focus were clarified, which provides a scientific basis for the effective prevention and control of human plague. This article also provides a set of paradigms for the systematic study of new plague foci, which is a perfect combination of traditional monitoring methods and modern research methods.

Genome-wide analysis of structural variants reveals genetic differences in Chinese pigs

Pigs have experienced long-term selections, resulting in dramatic phenotypic changes. Structural variants (SVs) are reported to exert extensive impacts on phenotypic changes. We built a high resolution and informative SV map based on high-depth sequencing data from 66 Chinese domestic and wild pigs. We inferred the SV formation mechanisms in the pig genome and used SVs as materials to perform a population-level analysis. We detected the selection signals on chromosome X for northern Chinese domestic pigs, as well as the differentiated loci across the whole genome. Analysis showed that these loci differ between southern and northern Chinese domestic pigs. Our results based on SVs provide new insights into genetic differences in Chinese pigs.

Population and evolutionary dynamics of Shiga-toxin producing Escherichia coli O157 in a beef herd: A longitudinal study

Shiga toxin producing Escherichia coli O157:H7 (STEC O157) is naturally found in the gastrointestinal tract of cattle and can cause severe disease in humans. There is limited understanding of the population dynamics and microevolution of STEC O157 at herd level. In this study, isolates from a closed beef herd of 23 cows were used to examine the population turnover in the herd. Of the nine STEC O157 clades previously described, clade 7 was found in 162 of the 169 isolates typed. Multiple locus variable number tandem repeat analysis (MLVA) differentiated 169 isolates into 33 unique MLVA types. Five predominant MLVA types were evident with most of the remaining types containing only a single isolate. MLVA data suggest that over time clonal replacement occurred within the herd. Genome sequencing of 18 selected isolates found that the isolates were divided into four lineages, representing four different 'clones' in the herd. Genome data confirmed clonal replacement over time and provided evidence of cross transmission of strains between cows. The findings enhanced our understanding of the population dynamics of STEC O157 in its natural host that will help developing effective control measures to prevent the spread of the pathogen to the human population.

Evolution of Variable Number Tandem Repeats and Its Relationship with Genomic Diversity in Salmonella Typhimurium

Salmonella enterica serovar Typhimurium is the most common Salmonella serovar causing human infections in Australia and many other countries. A total of 12,112 S. Typhimurium isolates from New South Wales were analyzed by multi-locus variable number of tandem repeat (VNTR) analysis (MLVA) using five VNTRs from 2007 to 2014. We found that mid ranges of repeat units of 8–14 in VNTR locus STTR5, 6–13 in STTR6, and 9–12 in STTR10 were always predominant in the population (>50%). In vitro passaging experiments using MLVA type carrying extreme length alleles found that the majority of long length alleles mutated to short ones and short length alleles mutated to longer ones. Both data suggest directional mutability of VNTRs toward mid-range repeats. Sequencing of 28 isolates from a newly emerged MLVA type and its five single locus variants revealed that single nucleotide variation between isolates with up to two MLVA differences ranged from 0 to 12 single nucleotide polymorphisms (SNPs). However, there was no relationship between SNP and VNTR differences. A population genetic model of the joint distribution of VNTRs and SNPs variations was used to estimate the mutation rates of the two markers, yielding a ratio of 1 VNTR change to 6.9 SNP changes. When only one VNTR repeat difference was considered, the majority of pairwise SNP difference between isolates were 4 SNPs or fewer. Based on this observation and our previous findings of SNP differences of outbreak isolates, we suggest that investigation of S. Typhimurium community outbreaks should include cases of 1 repeat difference to increase sensitivity. This study offers new insights into the short-term VNTR evolution of S. Typhimurium and its application for epidemiological typing.

Multiple Locus Variable number of tandem repeat Analysis: A molecular genotyping tool for Paenibacillus larvae

American Foulbrood, caused by Paenibacillus larvae, is the most severe bacterial disease of honey bees (Apis mellifera). To perform genotyping of P. larvae in an epidemiological context, there is a need of a fast and cheap method with a high resolution. Here, we propose Multiple Locus Variable number of tandem repeat Analysis (MLVA). MLVA has been used for typing a collection of 209 P. larvae strains from which 23 different MLVA types could be identified. Moreover, the developed methodology not only permits the identification of the four Enterobacterial Repetitive Intergenic Consensus (ERIC) genotypes, but allows also a discriminatory subdivision of the most dominant ERIC type I and ERIC type II genotypes. A biogeographical study has been conducted showing a significant correlation between MLVA genotype and the geographical region where it was isolated.

VNTR diversity in Yersinia pestis isolates from an animal challenge study reveals the potential for in vitro mutations during laboratory cultivation

Underlying mutation rates and other evolutionary forces shape the population structure of bacteria in nature. Although easily overlooked, similar forces are at work in the laboratory and may influence observed mutations. Here, we investigated tissue samples and Yersinia pestis isolates from a rodent laboratory challenge with strain CO92 using whole genome sequencing and multi-locus variable-number tandem repeat (VNTR) analysis (MLVA). We identified six VNTR mutations that were found to have occurred in vitro during laboratory cultivation rather than in vivo during the rodent challenge. In contrast, no single nucleotide polymorphism (SNP) mutations were observed, either in vivo or in vitro. These results were consistent with previously published mutation rates and the calculated number of Y. pestis generations that occurred during the in vitro versus the in vivo portions of the experiment. When genotyping disease outbreaks, the potential for in vitro mutations should be considered, particularly when highly variable genetic markers such as VNTRs are used.

Draft Genome of the Wheat Rust Pathogen (Puccinia triticina) Unravels Genome-Wide Structural Variations during Evolution

Leaf rust is one of the most important diseases of wheat and is caused by Puccinia triticina, a highly variable rust pathogen prevalent worldwide. Decoding the genome of this pathogen will help in unraveling the molecular basis of its evolution and in the identification of genes responsible for its various biological functions. We generated high quality draft genome sequences (approximately 100- 106 Mb) of two races of P. triticina; the variable and virulent Race77 and the old, avirulent Race106. The genomes of races 77 and 106 had 33X and 27X coverage, respectively. We predicted 27678 and 26384 genes, with average lengths of 1,129 and 1,086 bases in races 77 and 106, respectively and found that the genomes consisted of 37.49% and 39.99% repetitive sequences. Genome wide comparative analysis revealed that Race77 differs substantially from Race106 with regard to segmental duplication (SD), repeat element, and SNP/InDel characteristics. Comparative analyses showed that Race 77 is a recent, highly variable and adapted Race compared with Race106. Further sequence analyses of 13 additional pathotypes of Race77 clearly differentiated the recent, active and virulent, from the older pathotypes. Average densities of 2.4 SNPs and 0.32 InDels per kb were obtained for all P. triticina pathotypes. Secretome analysis demonstrated that Race77 has more virulence factors than Race 106, which may be responsible for the greater degree of adaptation of this pathogen. We also found that genes under greater selection pressure were conserved in the genomes of both races, and may affect functions crucial for the higher levels of virulence factors in Race77. This study provides insights into the genome structure, genome organization, molecular basis of variation, and pathogenicity of P. triticina The genome sequence data generated in this study have been submitted to public domain databases and will be an important resource for comparative genomics studies of the more than 4000 existing Puccinia species.

Two Distinct Yersinia pestis Populations Causing Plague among Humans in the West Nile Region of Uganda

Background

Plague is a life-threatening disease caused by the bacterium, Yersinia pestis. Since the 1990s, Africa has accounted for the majority of reported human cases. In Uganda, plague cases occur in the West Nile region, near the border with Democratic Republic of Congo. Despite the ongoing risk of contracting plague in this region, little is known about Y. pestis genotypes causing human disease.

Methodology/Principal Findings

During January 2004–December 2012, 1,092 suspect human plague cases were recorded in the West Nile region of Uganda. Sixty-one cases were culture-confirmed. Recovered Y. pestis isolates were analyzed using three typing methods, single nucleotide polymorphisms (SNPs), pulsed field gel electrophoresis (PFGE), and multiple variable number of tandem repeat analysis (MLVA) and subpopulations analyzed in the context of associated geographic, temporal, and clinical data for source patients. All three methods separated the 61 isolates into two distinct 1.ANT lineages, which persisted throughout the 9 year period and were associated with differences in elevation and geographic distribution.

Conclusions/Significance

We demonstrate that human cases of plague in the West Nile region of Uganda are caused by two distinct 1.ANT genetic subpopulations. Notably, all three typing methods used, SNPs, PFGE, and MLVA, identified the two genetic subpopulations, despite recognizing different mutation types in the Y. pestis genome. The geographic and elevation differences between the two subpopulations is suggestive of their maintenance in highly localized enzootic cycles, potentially with differing vector-host community composition. This improved understanding of Y. pestis subpopulations in the West Nile region will be useful for identifying ecologic and environmental factors associated with elevated plague risk.

Plague, a severe and often fatal zoonotic disease, is caused by the bacterium Yersinia pestis. Currently, the majority of human cases have been reported from resource limited areas of Africa, where the proximity to commensal rats and other small mammals increases the likelihood for human contact with infected animals or their fleas. Over a 9 year time period, >1000 suspect cases were recorded in the West Nile region of Uganda within the districts of Arua and Zombo. Culture-confirmed cases were shown by three independent typing methods to be due to two distinct 1.ANT genetic subpopulations of Y. pestis. The two genetic subpopulations persisted throughout the 9 year time period, consistent with their ongoing maintenance in local enzootic cycles. Additionally, the two subpopulations were found to differ with respect to geographic location and elevation, with SNP Group 1 strains being found further north and at lower elevations as compared to SNP Group 2. The relative independence of the two Y. pestis subpopulations is suggestive of their maintenance in distinct foci involving enzootic cycles with differing vector-host community composition.

Structural Variant Detection by Large-scale Sequencing Reveals New Evolutionary Evidence on Breed Divergence between Chinese and European Pigs

In this study, we performed a genome-wide SV detection among the genomes of thirteen pigs from diverse Chinese and European originated breeds by next genetation sequencing, and constrcuted a single-nucleotide resolution map involving 56,930 putative SVs. We firstly identified a SV hotspot spanning 35 Mb region on the X chromosome specifically in the genomes of Chinese originated individuals. Further scrutinizing this region by large-scale sequencing data of extra 111 individuals, we obtained the confirmatory evidence on our initial finding. Moreover, thirty five SV-related genes within the hotspot region, being of importance for reproduction ability, rendered significant different evolution rates between Chinese and European originated breeds. The SV hotspot identified herein offers a novel evidence for assessing phylogenetic relationships, as well as likely explains the genetic difference of corresponding phenotypes and features, among Chinese and European pig breeds. Furthermore, we employed various SVs to infer genetic structure of individuls surveyed. We found SVs can clearly detect the difference of genetic background among individuals. This clues us that genome-wide SVs can capture majority of geneic variation and be applied into cladistic analyses. Characterizing whole genome SVs demonstrated that SVs are significantly enriched/depleted with various genomic features.

Taxonomy of Yersinia pestis

This chapter summarized the taxonomy and typing works of Yersinia pestis since it's firstly identified in Hong Kong in 1894. Phenotyping methods that based on phenotypic characteristics, including biotyping, serotyping, antibiogram analysis, bacteriocin typing, phage typing, and plasmid typing, were firstly applied in classification of Y. pestis in subspecies level. And then, with the advancement of molecular biological technology, the methods based on outer membrane protein profiles, fatty acid composition, and bacterial mass fingerprinting were also used to identify the populations within Y. pestis. However, Y. pestis is a highly homogenous species; therefore, the above typing methods could only provide low resolution, e.g., only one serotype and one phage type were observed for the whole species. Since the 1990s, molecular typing based on DNA variations, including single-nucleotide polymorphism, gene gain/loss, variable-number tandem repeats, clustered regularly interspaced short palindromic repeat, etc., was introduced and improved the resolution and robust of typing result. Especially in recent years, genotyping-based whole-genome-wide variations were successfully employed in Y. pestis, which built the "gold standard" of typing scheme of the species and could distinguish the samples under the strain level. The taxonomy and typing works leaved us enormous polymorphism data; therefore, a comprehensive fingerprint database of Y. pestis was needed to collect and standardize these data, for facilitating future works on evolution, plague surveillance and control, anti-bioterrorism, and microbial forensic researches.

Whole genome multilocus sequence typing as an epidemiologic tool for Yersinia pestis

Human plague is a severe and often fatal zoonotic disease caused by Yersinia pestis. For public health investigations of human cases, nonintensive whole genome molecular typing tools, capable of defining epidemiologic relationships, are advantageous. Whole genome multilocus sequence typing (wgMLST) is a recently developed methodology that simplifies genomic analyses by transforming millions of base pairs of sequence into character data for each gene. We sequenced 13 US Y. pestis isolates with known epidemiologic relationships. Sequences were assembled de novo, and multilocus sequence typing alleles were assigned by comparison against 3979 open reading frames from the reference strain CO92. Allele-based cluster analysis accurately grouped the 13 isolates, as well as 9 publicly available Y. pestis isolates, by their epidemiologic relationships. Our findings indicate wgMLST is a simplified, sensitive, and scalable tool for epidemiologic analysis of Y. pestis strains.

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.

Development of Multiple-Locus Variable-Number Tandem-Repeat Analysis for Molecular Subtyping of Campylobacter jejuni by Using Capillary Electrophoresis

Campylobacter jejuni is a common cause of the frequently reported food-borne diseases in developed and developing nations. This study describes the development of multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) using capillary electrophoresis as a novel typing method for microbial source tracking and epidemiological investigation of C. jejuni. Among 36 tandem repeat loci detected by the Tandem Repeat Finder program, 7 VNTR loci were selected and used for characterizing 60 isolates recovered from chicken meat samples from retail shops, samples from chicken meat processing factory, and stool samples. The discrimination ability of MLVA was compared with that of multilocus sequence typing (MLST). MLVA (diversity index of 0.97 with 31 MLVA types) provided slightly higher discrimination than MLST (diversity index of 0.95 with 25 MLST types). The overall concordance between MLVA and MLST was estimated at 63% by adjusted Rand coefficient. MLVA predicted MLST type better than MLST predicted MLVA type, as reflected by Wallace coefficient (Wallace coefficient for MLVA to MLST versus MLST to MLVA, 86% versus 51%). MLVA is a useful tool and can be used for effective monitoring of C. jejuni and investigation of epidemics caused by C. jejuni.

Sequence Analysis of lip R: A Good Method for Molecular Epidemiology of Clinical Isolates of Mycobacterium tuberculosis

Advances in DNA sequencing have greatly enhanced the molecular epidemiology studies. In order to assess evolutionary and phylogenetic relation of Mycobacterium tuberculosis isolates several gene targets were evaluated. In this study, appropriate fragments of 5 highly variable genes (rpsL, mprA, lipR, katG, and fgd1 genes) were sequenced. The sequence data were analyzed with neighbor-joining method using mega and Geneious software. The phylogenetic trees analyzes revealed that the discriminatory power of lipR is much stronger than that observed in the other genes. lipR could distinguish between more clinical isolates. Therefore, lipR is a promising target for sequence analyzes of M. tuberculosis.

Single-Nucleotide Polymorphisms Reveal Spatial Diversity Among Clones of Yersinia pestis During Plague Outbreaks in Colorado and the Western United States

BACKGROUND

In western North America, plague epizootics caused by Yersinia pestis appear to sweep across landscapes, primarily infecting and killing rodents, especially ground squirrels and prairie dogs. During these epizootics, the risk of Y. pestis transmission to humans is highest. While empirical models that include climatic conditions and densities of rodent hosts and fleas can predict when epizootics are triggered, bacterial transmission patterns across landscapes, and the scale at which Y. pestis is maintained in nature during inter-epizootic periods, are poorly defined. Elucidating the spatial extent of Y. pestis clones during epizootics can determine whether bacteria are propagated across landscapes or arise independently from local inter-epizootic maintenance reservoirs.

MATERIAL AND METHODS

We used DNA microarray technology to identify single-nucleotide polymorphisms (SNPs) in 34 Y. pestis isolates collected in the western United States from 1980 to 2006, 21 of which were collected during plague epizootics in Colorado. Phylogenetic comparisons were used to elucidate the hypothesized spread of Y. pestis between the mountainous Front Range and the eastern plains of northern Colorado during epizootics. Isolates collected from across the western United States were included for regional comparisons.

RESULTS

By identifying SNPs that mark individual clones, our results strongly suggest that Y. pestis is maintained locally and that widespread epizootic activity is caused by multiple clones arising independently at small geographic scales. This is in contrast to propagation of individual clones being transported widely across landscapes. Regionally, our data are consistent with the notion that Y. pestis diversifies at relatively local scales following long-range translocation events. We recommend that surveillance and prediction by public health and wildlife management professionals focus more on models of local or regional weather patterns and ecological factors that may increase risk of widespread epizootics, rather than predicting or attempting to explain epizootics on the basis of movement of host species that may transport plague.

Microbial Forensics☆

Biothreats are a high priority concern for public safety and national security. The field of microbial forensics was developed to analyze evidence associated with biological crimes in which microbes or their toxins are used as weapons. Microbial forensics is the scientific discipline dedicated to analyzing evidence from a bioterrorism act, biocrime, hoax, or inadvertent microorganism/toxin release for attribution purposes. Microbial forensics combines the practices of epidemiology with the characterization of microbial and microbial-related evidence to assist in determining the specific source of the sample, as individualizing as possible, and/or the methods, means, processes and locations involved to determine the identity of the perpetrator(s) of an attack.

Insights to genetic characterization tools for epidemiological tracking of Francisella tularensis in Sweden

Tularaemia, caused by the bacterium Francisella tularensis, is endemic in Sweden and is poorly understood. The aim of this study was to evaluate the effectiveness of three different genetic typing systems to link a genetic type to the source and place of tularemia infection in Sweden. Canonical single nucleotide polymorphisms (canSNPs), MLVA including five variable number of tandem repeat loci and PmeI-PFGE were tested on 127 F. tularensis positive specimens collected from Swedish case-patients. All three typing methods identified two major genetic groups with near-perfect agreement. Higher genetic resolution was obtained with canSNP and MLVA compared to PFGE; F. tularensis samples were first assigned into ten phylogroups based on canSNPs followed by 33 unique MLVA types. Phylogroups were geographically analysed to reveal complex phylogeographic patterns in Sweden. The extensive phylogenetic diversity found within individual counties posed a challenge to linking specific genetic types with specific geographic locations. Despite this, a single phylogroup (B.22), defined by a SNP marker specific to a lone Swedish sequenced strain, did link genetic type with a likely geographic place. This result suggests that SNP markers, highly specific to a particular reference genome, may be found most frequently among samples recovered from the same location where the reference genome originated. This insight compels us to consider whole-genome sequencing (WGS) as the appropriate tool for effectively linking specific genetic type to geography. Comparing the WGS of an unknown sample to WGS databases of archived Swedish strains maximizes the likelihood of revealing those rare geographically informative SNPs.

Limitations to estimating bacterial cross-species transmission using genetic and genomic markers: inferences from simulation modeling

Cross-species transmission (CST) of bacterial pathogens has major implications for human health, livestock, and wildlife management because it determines whether control actions in one species may have subsequent effects on other potential host species. The study of bacterial transmission has benefitted from methods measuring two types of genetic variation: variable number of tandem repeats (VNTRs) and single nucleotide polymorphisms (SNPs). However, it is unclear whether these data can distinguish between different epidemiological scenarios. We used a simulation model with two host species and known transmission rates (within and between species) to evaluate the utility of these markers for inferring CST. We found that CST estimates are biased for a wide range of parameters when based on VNTRs and a most parsimonious reconstructed phylogeny. However, estimations of CST rates lower than 5% can be achieved with relatively low bias using as low as 250 SNPs. CST estimates are sensitive to several parameters, including the number of mutations accumulated since introduction, stochasticity, the genetic difference of strains introduced, and the sampling effort. Our results suggest that, even with whole-genome sequences, unbiased estimates of CST will be difficult when sampling is limited, mutation rates are low, or for pathogens that were recently introduced.

Phylogeography and population structure of the biologically invasive phytopathogen Erwinia amylovora inferred using minisatellites

Erwinia amylovora causes a major disease of pome fruit trees worldwide, and is regulated as a quarantine organism in many countries. While some diversity of isolates has been observed, molecular epidemiology of this bacterium is hindered by a lack of simple molecular typing techniques with sufficiently high resolution. We report a molecular typing system of E. amylovora based on variable number of tandem repeats (VNTR) analysis. Repeats in the E. amylovora genome were identified with comparative genomic tools, and VNTR markers were developed and validated. A Multiple-Locus VNTR Analysis (MLVA) was applied to E. amylovora isolates from bacterial collections representing global and regional distribution of the pathogen. Based on six repeats, MLVA allowed the distinction of 227 haplotypes among a collection of 833 isolates of worldwide origin. Three geographically separated groups were recognized among global isolates using Bayesian clustering methods. Analysis of regional outbreaks confirmed presence of diverse haplotypes but also high representation of certain haplotypes during outbreaks. MLVA analysis is a practical method for epidemiological studies of E. amylovora, identifying previously unresolved population structure within outbreaks. Knowledge of such structure can increase our understanding on how plant diseases emerge and spread over a given geographical region.

Dynamic of mutational events in variable number tandem repeats of Escherichia coli O157:H7

VNTRs regions have been successfully used for bacterial subtyping; however, the hypervariability in VNTR loci is problematic when trying to predict the relationships among isolates. Since few studies have examined the mutation rate of these markers, our aim was to estimate mutation rates of VNTRs specific for verotoxigenic E. coli O157:H7. The knowledge of VNTR mutational rates and the factors affecting them would make MLVA more effective for epidemiological or microbial forensic investigations. For this purpose, we analyzed nine loci performing parallel, serial passage experiments (PSPEs) on 9 O157:H7 strains. The combined 9 PSPE population rates for the 8 mutating loci ranged from 4.4 × 10⁻⁰⁵ to 1.8 × 10⁻⁰³ mutations/generation, and the combined 8-loci mutation rate was of 2.5 × 10⁻⁰³ mutations/generation. Mutations involved complete repeat units, with only one point mutation detected. A similar proportion between single and multiple repeat changes was detected. Of the 56 repeat mutations, 59% were insertions and 41% were deletions, and 72% of the mutation events corresponded to O157-10 locus. For alleles with up to 13 UR, a constant and low mutation rate was observed; meanwhile longer alleles were associated with higher and variable mutation rates. Our results are useful to interpret data from microevolution and population epidemiology studies and particularly point out that the inclusion or not of O157-10 locus or, alternatively, a differential weighting data according to the mutation rates of loci must be evaluated in relation with the objectives of the proposed study.

Multiplex variable number of tandem repeats for Oenococcus oeni and applications

Oenococcus oeni is responsible for the malolactic fermentation of wine. Genomic diversity has already been established in this species. In addition, winemakers usually report varying starter-culture efficiency. It is essential to monitor indigenous and selected strains in order to understand strain survival and development during the winemaking process. A previous article described a variable number of tandem repeats (VNTR) scheme, based on five polymorphic loci of the genome. VNTR typing of O. oeni was highly discriminating, faster, and more reliable than the PFGE or MLST methods. The objective of this study was to set up a faster protocol by multiplexing, taking advantage of the high performance of multicolor capillary electrophoresis. The primers were labeled with multiple fluorescent dyes. PCR conditions were adapted by multiplexing amplifications in two separate PCR mixtures for the five loci, both at the same annealing temperature. The resulting assay proved to be robust, accurate, fast and easy to perform. Thanks to this new protocol, all O. oeni strains used in the study were typed using the five tandem repeats (TR). As expected, the primers for the five TR loci were specific to O. oeni. The method was improved to analyze isolated and mixed colonies, as well as bacteria harvested from wine using fast technology for analysis of nucleic acids (FTA(®)) technology. Finally, predictive models were constructed, to predict phylogenetic relationships and associate bacterial strain resistance to freeze-drying with fragment length analysis (FLA) profiles and genotypic and phenotypic characters.

Features of Variable Number of Tandem Repeats in Yersinia pestis and the Development of a Hierarchical Genotyping Scheme

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.

Multilocus variable-number-tandem-repeats analysis (MLVA) distinguishes a clonal complex of Clavibacter michiganensis subsp. michiganensis strains isolated from recent outbreaks of bacterial wilt and canker in Belgium

BACKGROUND

Clavibacter michiganensis subsp. michiganensis (Cmm) causes bacterial wilt and canker in tomato. Cmm is present nearly in all European countries. During the last three years several local outbreaks were detected in Belgium. The lack of a convenient high-resolution strain-typing method has hampered the study of the routes of transmission of Cmm and epidemiology in tomato cultivation. In this study the genetic relatedness among a worldwide collection of Cmm strains and their relatives was approached by gyrB and dnaA gene sequencing. Further, we developed and applied a multilocus variable number of tandem repeats analysis (MLVA) scheme to discriminate among Cmm strains.

RESULTS

A phylogenetic analysis of gyrB and dnaA gene sequences of 56 Cmm strains demonstrated that Belgian Cmm strains from recent outbreaks of 2010-2012 form a genetically uniform group within the Cmm clade, and Cmm is phylogenetically distinct from other Clavibacter subspecies and from non-pathogenic Clavibacter-like strains. MLVA conducted with eight minisatellite loci detected 25 haplotypes within Cmm. All strains from Belgian outbreaks, isolated between 2010 and 2012, together with two French strains from 2010 seem to form one monomorphic group. Regardless of the isolation year, location or tomato cultivar, Belgian strains from recent outbreaks belonged to the same haplotype. On the contrary, strains from diverse geographical locations or isolated over longer periods of time formed mostly singletons.

CONCLUSIONS

We hypothesise that the introduction might have originated from one lot of seeds or contaminated tomato seedlings that was the source of the outbreak in 2010 and that these Cmm strains persisted and induced infection in 2011 and 2012. Our results demonstrate that MLVA is a promising typing technique for a local surveillance and outbreaks investigation in epidemiological studies of Cmm.

Development of variable number of tandem repeats typing schemes for Ralstonia solanacearum, the agent of bacterial wilt, banana Moko disease and potato brown rot

Ralstonia solanacearum is an important soil borne bacterial plant pathogen causing bacterial wilt on many important crops. To better monitor epidemics, efficient tools that can identify and discriminate populations are needed. In this study, we assessed variable number of tandem repeats (VNTR) genotyping as a new tool for epidemiological surveillance of R. solanacearum phylotypes, and more specifically for the monitoring of the monomorphic ecotypes "Moko" (banana-pathogenic) and "brown rot" (potato-pathogenic under cool conditions). Screening of six R. solanacearum genome sequences lead to select 36 VNTR loci that were preliminarily amplified on 24 strains. From this step, 26 single-locus primer pairs were multiplexed, and applied to a worldwide collection of 337 strains encompassing the whole phylogenetic diversity, with revelation on a capillary-electrophoresis genotype. Four loci were monomorphic within all phylotypes and were not retained; the other loci were highly polymorphic but displayed a clear phylotype-specificity. Phylotype-specific MLVA schemes were thus defined, based on 13 loci for phylotype I, 12 loci for phylotype II, 11 loci for phylotype III and 6 for phylotype IV. MLVA typing was significantly more discriminative than egl-based sequevar typing, particularly on monomorphic "brown rot" ecotype (phylotype IIB/sequevar 1) and "Moko disease" clade 4 (Phylotype IIB/sequevar 4). Our results raise promising prospects for studies of population genetic structures and epidemiological monitoring.

Low numbers of repeat units in variable number of tandem repeats (VNTR) regions of white spot syndrome virus are correlated with disease outbreaks

White spot syndrome virus (WSSV) is the most important pathogen in shrimp farming systems worldwide including the Mekong Delta, Vietnam. The genome of WSSV is characterized by the presence of two major 'indel regions' found at ORF14/15 and ORF23/24 (WSSV-Thailand) and three regions with variable number tandem repeats (VNTR) located in ORF75, ORF94 and ORF125. In the current study, we investigated whether or not the number of repeat units in the VNTRs correlates with virus outbreak status and/or shrimp farming practice. We analysed 662 WSSV samples from individual WSSV-infected Penaeus monodon shrimp from 104 ponds collected from two important shrimp farming regions of the Mekong Delta: Ca Mau and Bac Lieu. Using this large data set and statistical analysis, we found that for ORF94 and ORF125, the mean number of repeat units (RUs) in VNTRs was significantly lower in disease outbreak ponds than in non-outbreak ponds. Although a higher mean RU number was observed in the improved-extensive system than in the rice-shrimp or semi-intensive systems, these differences were not significant. VNTR sequences are thus not only useful markers for studying WSSV genotypes and populations, but specific VNTR variants also correlate with disease outbreaks in shrimp farming systems.

Genotyping of present-day and historical Geobacillus species isolates from milk powders by high-resolution melt analysis of multiple variable-number tandem-repeat loci

Spores of thermophilic Geobacillus species are a common contaminant of milk powder worldwide due to their ability to form biofilms within processing plants. Genotyping methods can provide information regarding the source and monitoring of contamination. A new genotyping method was developed based on multilocus variable-number tandem-repeat (VNTR) analysis (MLVA) in conjunction with high-resolution melt analysis (MLV-HRMA) and compared to the currently used method, randomized amplified polymorphic DNA PCR (RAPD-PCR). Four VNTR loci were identified and used to genotype 46 Geobacillus isolates obtained from retailed powder and samples from 2 different milk powder processing plants. These 46 isolates were differentiated into 16 different groups using MLV-HRMA (D = 0.89). In contrast, only 13 RAPD-PCR genotypes were identified among the 46 isolates (D = 0.79). This new method was then used to analyze 35 isolates obtained from powders with high spore counts (>10(4) spores · g(-1)) from a single processing plant together with 27 historical isolates obtained from powder samples processed in the same region of Australia 17 years ago. Results showed that three genotypes can coexist in a single processing run, while the same genotypes observed 17 years ago are present today. While certain genotypes could be responsible for powders with high spore counts, there was no correlation to specific genotypes being present in powder plants and retailed samples. In conclusion, the MLV-HRMA method is useful for genotyping Geobacillus spp. to provide insight into the prevalence and persistence of certain genotypes within milk powder processing plants.

Highly discriminatory typing method for Listeria monocytogenes using polymorphic tandem repeat regions

Tandem repeats (TR), which are repetitive nucleotide sequences in DNA, are polymorphic both in repeat number and sequence. In this study, we developed a new typing method, multilocus TR sequence analysis (MLTSA), for the foodborne pathogen Listeria monocytogenes using sequence polymorphisms in three tandem repeat regions. The obtained dendrogram clustered L. monocytogenes strains of lineage I and lineage II separately, and formed three groups within the lineage I cluster, each of which included one of the three major L. monocytogenes epidemic clones (ECI, ECIa, and ECII). These results were consistent with a previously established virulence-gene-based MLST method. In comparison, our method grouped some epidemiologically related isolates together, which virulence-gene-based MLST did not. Moreover, our method, using three tandem repeat regions, showed a higher discriminatory power than the MLST method, which uses six virulence gene regions. This MLTSA approach using sequence polymorphisms in TR regions could be a useful tool in the epidemiological study of L. monocytogenes.

Microfluidic-chip-based multiple-locus variable-number tandem-repeat fingerprinting with new primer sets for methicillin-resistant Staphylococcus aureus

The detection of outbreaks of methicillin-resistant Staphylococcus aureus (MRSA) infections and a rapid and accurate identification of sources and routes of transmission should be conducted in hospital settings as early and swiftly as possible. In this study, we investigated the application potential of a new approach based on multiple-locus variable-number tandem-repeat fingerprinting (MLVF) and microfluidics technology for a rapid discrimination of MRSA lineages in outbreak settings. A total of 206 nonrepetitive MRSA isolates recovered from infected patients at the University Medical Center Groningen between 2000 and 2010 were tested. The results obtained by MLVF using microcapillary electrophoresis with newly designed primers were compared to those obtained by spa typing and multiple-locus variable-number tandem-repeat analysis (MLVA). The discriminatory power was 0.980 (107 patterns), 0.969 (85 allelic profiles), and 0.959 (66 types) for MLVF, MLVA, and spa typing, respectively. All methods tested showed a good concordance of results calculated by the adjusted Rand's coefficient method. Comparisons of data obtained by the three approaches allowed us to propose an 88% cutoff value for the similarity between any two MLVF patterns, which can be used in S. aureus epidemiological studies, including analyses of outbreaks and strain transmission events. Of the three tested methods, MLVF is the cheapest, fastest, and easiest to perform. MLVF applied to microfluidic polymer chips is a rapid, cheap, reproducible, and highly discriminating tool to determine the clonality of MRSA isolates and to trace the spread of MRSA strains over periods of many years. Although spa typing should be used due to its portability of data, MLVF has a high added value because it is more discriminatory.

Global Wolbachia prevalence, titer fluctuations and their potential of causing cytoplasmic incompatibilities in tsetse flies and hybrids of Glossina morsitans subgroup species

We demonstrate the high applicability of a novel VNTR-based (Variable-Number-Tandem-Repeat) molecular screening tool for fingerprinting Wolbachia-infections in tsetse flies. The VNTR-141 locus provides reliable and concise differentiation between Wolbachia strains deriving from Glossina morsitans morsitans, Glossina morsitans centralis, and Glossina brevipalpis. Moreover, we show that certain Wolbachia-infections in Glossina spp. are capable of escaping standard PCR screening methods by 'hiding' as low-titer infections below the detection threshold. By applying a highly sensitive PCR-blot technique to our Glossina specimen, we were able to enhance the symbiont detection limit substantially and, consequently, trace unequivocally Wolbachia-infections at high prevalence in laboratory-reared G. swynnertoni individuals. To our knowledge, Wolbachia-persistence was reported exclusively for field-collected samples, and at low prevalence only. Finally, we highlight the substantially higher Wolbachia titer levels found in hybrid Glossina compared to non-hybrid hosts and the possible impact of these titers on hybrid host fitness that potentially trigger incipient speciation in tsetse flies.

Yersinia pestis lineages in Mongolia

Background

Whole genome sequencing allowed the development of a number of high resolution sequence based typing tools for Yersinia (Y.) pestis. The application of these methods on isolates from most known foci worldwide and in particular from China and the Former Soviet Union has dramatically improved our understanding of the population structure of this species. In the current view, Y. pestis including the non or moderate human pathogen Y. pestis subspecies microtus emerged from Yersinia pseudotuberculosis about 2,600 to 28,600 years ago in central Asia. The majority of central Asia natural foci have been investigated. However these investigations included only few strains from Mongolia.

Methodology/Principal Findings

Clustered Regularly Interspaced Short Prokaryotic Repeats (CRISPR) analysis and Multiple-locus variable number of tandem repeats (VNTR) analysis (MLVA) with 25 loci was performed on 100 Y. pestis strains, isolated from 37 sampling areas in Mongolia. The resulting data were compared with previously published data from more than 500 plague strains, 130 of which had also been previously genotyped by single nucleotide polymorphism (SNP) analysis. The comparison revealed six main clusters including the three microtus biovars Ulegeica, Altaica, and Xilingolensis. The largest cluster comprises 78 isolates, with unique and new genotypes seen so far in Mongolia only. Typing of selected isolates by key SNPs was used to robustly assign the corresponding clusters to previously defined SNP branches.

Conclusions/Significance

We show that Mongolia hosts the most recent microtus clade (Ulegeica). Interestingly no representatives of the ancestral Y. pestis subspecies pestis nodes previously identified in North-western China were identified in this study. This observation suggests that the subsequent evolution steps within Y. pestis pestis did not occur in Mongolia. Rather, Mongolia was most likely re-colonized by more recent clades coming back from China contemporary of the black death pandemic, or more recently in the past 600 years.

Comparative genomics of 2009 seasonal plague (Yersinia pestis) in New Mexico

Plague disease caused by the Gram-negative bacterium Yersinia pestis routinely affects animals and occasionally humans, in the western United States. The strains native to the North American continent are thought to be derived from a single introduction in the late 19^th century. The degree to which these isolates have diverged genetically since their introduction is not clear, and new genomic markers to assay the diversity of North American plague are highly desired. To assay genetic diversity of plague isolates within confined geographic areas, draft genome sequences were generated by 454 pyrosequencing from nine environmental and clinical plague isolates. In silico assemblies of Variable Number Tandem Repeat (VNTR) loci were compared to laboratory-generated profiles for seven markers. High-confidence SNPs and small Insertion/Deletions (Indels) were compared to previously sequenced Y. pestis isolates. The resulting panel of mutations allowed clustering of the strains and tracing of the most likely evolutionary trajectory of the plague strains. The sequences also allowed the identification of new putative SNPs that differentiate the 2009 isolates from previously sequenced plague strains and from each other. In addition, new insertion points for the abundant insertion sequences (IS) of Y. pestis are present that allow additional discrimination of strains; several of these new insertions potentially inactivate genes implicated in virulence. These sequences enable whole-genome phylogenetic analysis and allow the unbiased comparison of closely related isolates of a genetically monomorphic pathogen.

Phylogeography and molecular epidemiology of Yersinia pestis in Madagascar

Background

Plague was introduced to Madagascar in 1898 and continues to be a significant human health problem. It exists mainly in the central highlands, but in the 1990s was reintroduced to the port city of Mahajanga, where it caused extensive human outbreaks. Despite its prevalence, the phylogeography and molecular epidemiology of Y. pestis in Madagascar has been difficult to study due to the great genetic similarity among isolates. We examine island-wide geographic-genetic patterns based upon whole-genome discovery of SNPs, SNP genotyping and hypervariable variable-number tandem repeat (VNTR) loci to gain insight into the maintenance and spread of Y. pestis in Madagascar.

Methodology/Principal Findings

We analyzed a set of 262 Malagasy isolates using a set of 56 SNPs and a 43-locus multi-locus VNTR analysis (MLVA) system. We then analyzed the geographic distribution of the subclades and identified patterns related to the maintenance and spread of plague in Madagascar. We find relatively high levels of VNTR diversity in addition to several SNP differences. We identify two major groups, Groups I and II, which are subsequently divided into 11 and 4 subclades, respectively. Y. pestis appears to be maintained in several geographically separate subpopulations. There is also evidence for multiple long distance transfers of Y. pestis, likely human mediated. Such transfers have resulted in the reintroduction and establishment of plague in the port city of Mahajanga, where there is evidence for multiple transfers both from and to the central highlands.

Conclusions/Significance

The maintenance and spread of Y. pestis in Madagascar is a dynamic and highly active process that relies on the natural cycle between the primary host, the black rat, and its flea vectors as well as human activity.

Plague, caused by the bacterium Yersinia pestis, has been a problem in Madagascar since it was introduced in 1898. It mainly affects the central highlands, but also has caused several large outbreaks in the port city of Mahajanga, after it was reintroduced there in the 1990s. Despite its prevalence, the genetic diversity and related geographic distribution of different genetic groups of Y. pestis in Madagascar has been difficult to study due to the great genetic similarity among isolates. We subtyped a set of Malagasy isolates and identified two major genetic groups that were subsequently divided into 11 and 4 subgroups, respectively. Y. pestis appears to be maintained in several geographically separate subpopulations. There is also evidence for multiple long distance transfers of Y. pestis, likely human mediated. Such transfers have resulted in the reintroduction and establishment of plague in the port city of Mahajanga where there is evidence for multiple transfers both from and to the central highlands. The maintenance and spread of Y. pestis in Madagascar is a dynamic and highly active process that relies on the natural cycle between the primary host, the black rat, and its flea vectors as well as human activity.