Phylogenetic understanding of clonal populations in an era of whole genome sequencing

Population sequencing for diversity and transmission analyses

Genomic diversity in a pathogen population is the foundation for evolution and adaptations in virulence, drug resistance, pathogenesis, and immune evasion. Characterizing, analyzing, and understanding population-level diversity is also essential for epidemiological and forensic tracking of sources and revealing detailed pathways of transmission and spread. For bacteria, culturing, isolating, and sequencing the large number of individual colonies required to adequately sample diversity can be prohibitively time-consuming and expensive. While sequencing directly from a mixed population will show variants among reads, they cannot be linked to reveal allele combinations associated with particular traits or phylogenetic inheritance patterns. Here, we describe the theory and method of how population sequencing directly from a mixed sample can be used in conjunction with sequencing a very small number of colonies to describe the phylogenetic diversity of a population without haplotype reconstruction. To demonstrate the utility of population sequencing in capturing phylogenetic diversity, we compared isogenic clones to population sequences of Burkholderia pseudomallei from the sputum of a single patient. We also analyzed population sequences of Staphylococcus aureus derived from different people and different body sites. Sequencing results confirm our ability to capture and characterize phylogenetic diversity in our samples. Our analyses of B. pseudomallei populations led to the surprising discovery that the pathogen population is highly structured in sputum, suggesting that for some pathogens, sputum sampling may preserve structuring in the lungs and thus present a non-invasive alternative to understanding colonization, movement, and pathogen/host interactions. Our analyses of S. aureus samples show how comparing phylogenetic diversity across populations can reveal directionality of transmission between hosts and across body sites, demonstrating the power and utility for characterizing the spread of disease and identification of reservoirs at the finest levels. We anticipate that population sequencing and analysis can be broadly applied to accelerate research in a broad range of fields reliant on a foundational understanding of population diversity.

In vitro and genomic mining studies of anti-Clostridium perfringens Compounds Derived from Bacillus amyloliquefaciens

Clostridium perfringens is an important opportunistic microorganism in commercial poultry production that is implicated in necrotic enteritis (NE) outbreaks. This disease poses a severe financial burden on the global poultry industry, causing estimated annual losses of $6 billion globally. The ban on in-feed antibiotic growth promoters has spurred investigations into approaches of alternatives to antibiotics, among which Bacillus probiotics have demonstrated varying degrees of effectiveness against NE. However, the precise mechanisms underlying Bacillus-mediated beneficial effects on host responses in NE remain to be further elucidated. In this manuscript, we conducted in vitro and genomic mining analysis to investigate anti-C. perfringens activity observed in the supernatants derived from 2 Bacillus amyloliquefaciens strains (FS1092 and BaD747). Both strains demonstrated potent anti-C. perfringens activities in in vitro studies. An analysis of genomes from 15 B. amyloliquefaciens, 11 B. velezensis, and 2 B. subtilis strains has revealed an intriguing clustering pattern among strains known to possess anti-C. perfringens activities. Furthermore, our investigation has identified 7 potential antimicrobial compounds, predicted as secondary metabolites through antiSMASH genomic mining within the published genomes of B. amyloliquefaciens species. Based on in vitro analysis, BaD747 may have the potential as a probiotic in the control of NE. These findings not only enhance our understanding of B. amyloliquefaciens's action against C. perfringens but also provide a scientific rationale for the development of novel antimicrobial therapeutic agents against NE.

Whole-Genome sequencing in routine Mycobacterium bovis epidemiology - scoping the potential

Mycobacterium bovis the main agent of bovine tuberculosis (bTB), presents as a series of spatially-localised micro-epidemics across landscapes. Classical molecular typing methods applied to these micro-epidemics, based on genotyping a few variable loci, have significantly improved our understanding of potential epidemiological links between outbreaks. However, they have limited utility owing to low resolution. Conversely, whole-genome sequencing (WGS) provides the highest resolution data available for molecular epidemiology, producing richer outbreak tracing, insights into phylogeography and epidemic evolutionary history. We illustrate these advantages by focusing on a common single lineage of M. bovis (1.140) from Northern Ireland. Specifically, we investigate the spatial sub-structure of 20 years of herd-level multi locus VNTR analysis (MLVA) surveillance data and WGS data from a down sampled subset of isolates of this MLVA type over the same time frame. We mapped 2108 isolate locations of MLVA type 1.140 over the years 2000-2022. We also mapped the locations of 148 contemporary WGS isolates from this lineage, over a similar geographic range, stratifying by single nucleotide polymorphism (SNP) relatedness cut-offs of 15 SNPs. We determined a putative core range for the 1.140 MLVA type and SNP-defined sequence clusters using a 50 % kernel density estimate, using cattle movement data to inform on likely sources of WGS isolates found outside of core ranges. Finally, we applied Bayesian phylogenetic methods to investigate past population history and reproductive number of the 1.140 M. bovis lineage. We demonstrate that WGS SNP-defined clusters exhibit smaller core ranges than the established MLVA type - facilitating superior disease tracing. We also demonstrate the superior functionality of WGS data in determining how this lineage was disseminated across the landscape, likely via cattle movement and to infer how its effective population size and reproductive number has been in flux since its emergence. These initial findings highlight the potential of WGS data for routine monitoring of bTB outbreaks.

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Spatial and phylogenetic patterns reveal hidden infection sources of Bacillus anthracis in an anthrax outbreak in Son La province, Vietnam

Bacillus anthracis, the bacterial cause of anthrax, is a zoonosis affecting livestock and wildlife often spilling over into humans. In Vietnam, anthrax has been nationally reportable since 2015 with cases occurring annually, mostly in the northern provinces. In April 2022, an outbreak was reported in Son La province following the butchering of a water buffalo, Bubalus bubalis. A total of 137 humans from three villages were likely exposed to contaminated meat from the animal. Early epidemiological investigations suggested a single animal was involved in all exposures. Five B. anthracis isolates were recovered from human clinical cases along with one from the buffalo hide, another from associated maggots, and one from soil at the carcass site. The isolates were whole genome sequenced, allowing global, regional, and local molecular epidemiological analyses of the outbreak strains. All recovered B. anthracis belong to the A.Br.001/002 lineage based on canonical single nucleotide polymorphism analysis (canSNP). Although not previously identified in Vietnam, this lineage has been identified in the nearby countries of China, India, Indonesia, Thailand, as well as Australia. A twenty-five marker multi-locus variable number tandem repeat analysis (MLVA-25) was used to investigate the relationship between human, soil, and buffalo strains. Locally, four MLVA-25 genotypes were identified from the eight isolates. This level of genetic diversity is unusual for the limited geography and timing of cases and differs from past literature using MLVA-25. The coupled spatial and phylogenetic data suggest this outbreak originated from multiple, likely undetected, animal sources. These findings were further supported by local news reports that identified at least two additional buffalo deaths beyond the initial animal sampled in response to the human cases. Future outbreak response should include intensive surveillance for additional animal cases and additional molecular epidemiological traceback to identify pathogen sources.

Geo-epidemiology of animal tuberculosis and Mycobacterium bovis genotypes in livestock in a small, high-incidence area in Sicily, Italy

Introduction

The persistence of animal tuberculosis (TB) in livestock is a major concern in Sicily, Italy. The objective of this study was to elucidate the transmission dynamics of M. bovis infection in a highly circumscribed, and at the same time geographically diverse, high-risk area of the island through an in-depth geo-epidemiological investigation of TB in cattle and black pigs raised in small-scale extensive farms across the district of Caronia.

Methods

We used genotype analysis coupled with geographic information system (GIS) technology and phylogenetic inference to characterize the spatial distribution of TB and M. bovis genotypes in livestock and the genetic relationships between M. bovis isolates. A total of 589 M. bovis isolates collected from slaughtered cattle (n = 527) and Sicilian black pigs (n = 62) over a 5-year period (2014–2018) were included in the study.

Results

TB was widespread throughout the district and was most frequent in the north-central area of the district, especially along one of the district’s streams. We identified a total of 62 M. bovis genotypes. Identical genetic profiles were isolated from both neighboring and non-neighburing herds. The 10 most frequent genotypes, accounting for 82% of M. bovis isolates, showed geographic specificities in that they tended to cluster in specific spatial niches. The landscape structure of these niches—i.e. steep slopes, rocky ridges, meadows and streams—is likely to have had a significant influence on the distribution of TB among livestock in Caronia. Higher concentrations of TB were observed along streams and in open meadows, while rocky ridges and slopes appeared to have hampered the spread of TB.

Discussion

The geographical distribution of TB cases among livestock in Caronia is consistent with several epidemiological scenarios (e.g., high density of infected herds along the streams or in hilly plateau where livestock share pastures). Landscape structure is likely to play an important role in the transmission and persistence of M. bovis infection across the district. Additional potential risk factors, such as livestock trading and extensive breeding methods, are also discussed. Our results will contribute to the improvement of surveillance, control and eradication activities of TB in Sicily by the implementation of ad hoc TB control measures, especially in farms located along streams, sharing common pastures or with mixed animal species.

Whole Genome Sequencing Refines Knowledge on the Population Structure of Mycobacterium bovis from a Multi-Host Tuberculosis System

Classical molecular analyses of Mycobacterium bovis based on spoligotyping and Variable Number Tandem Repeat (MIRU-VNTR) brought the first insights into the epidemiology of animal tuberculosis (TB) in Portugal, showing high genotypic diversity of circulating strains that mostly cluster within the European 2 clonal complex. Previous surveillance provided valuable information on the prevalence and spatial occurrence of TB and highlighted prevalent genotypes in areas where livestock and wild ungulates are sympatric. However, links at the wildlife-livestock interfaces were established mainly via classical genotype associations. Here, we apply whole genome sequencing (WGS) to cattle, red deer and wild boar isolates to reconstruct the M. bovis population structure in a multi-host, multi-region disease system and to explore links at a fine genomic scale between M. bovis from wildlife hosts and cattle. Whole genome sequences of 44 representative M. bovis isolates, obtained between 2003 and 2015 from three TB hotspots, were compared through single nucleotide polymorphism (SNP) variant calling analyses. Consistent with previous results combining classical genotyping with Bayesian population admixture modelling, SNP-based phylogenies support the branching of this M. bovis population into five genetic clades, three with apparent geographic specificities, as well as the establishment of an SNP catalogue specific to each clade, which may be explored in the future as phylogenetic markers. The core genome alignment of SNPs was integrated within a spatiotemporal metadata framework to further structure this M. bovis population by host species and TB hotspots, providing a baseline for network analyses in different epidemiological and disease control contexts. WGS of M. bovis isolates from Portugal is reported for the first time in this pilot study, refining the spatiotemporal context of TB at the wildlife-livestock interface and providing further support to the key role of red deer and wild boar on disease maintenance. The SNP diversity observed within this dataset supports the natural circulation of M. bovis for a long time period, as well as multiple introduction events of the pathogen in this Iberian multi-host system.

Emerging diversity and ongoing expansion of the genus Brucella

Remarkable genetic diversity and breadth of host species has been uncovered in the Brucella genus over the past decade, fundamentally changing our concept of what it means to be a Brucella. From ocean fishes and marine mammals, to pond dwelling amphibians, forest foxes, desert rodents, and cave-dwelling bats, Brucella have revealed a variety of previously unknown niches. Classical microbiological techniques have been able to help us classify many of these new strains but at times have limited our ability to see the true relationships among or within species. The closest relatives of Brucella are soil bacteria and the adaptations of Brucella spp. to live intracellularly suggest that the genus has evolved to live in vertebrate hosts. Several recently discovered species appear to have phenotypes that are intermediate between soil bacteria and core Brucella, suggesting that they may represent ancestral traits that were subsequently lost in the traditional species. Remarkably, the broad relationships among Brucella species using a variety of sequence and fragment-based approaches have been upheld when using comparative genomics with whole genomes. Nonetheless, genomes are required for fine-scale resolution of many of the relationships and for understanding the evolutionary history of the genus. We expect that the coming decades will reveal many more hosts and previously unknown diversity in a wide range of environments.

Pathogen to commensal? Longitudinal within-host population dynamics, evolution, and adaptation during a chronic >16-year Burkholderia pseudomallei infection

Although acute melioidosis is the most common outcome of Burkholderia pseudomallei infection, we have documented a case, P314, where disease severity lessened with time, and the pathogen evolved towards a commensal relationship with the host. In the current study, we used whole-genome sequencing to monitor this long-term symbiotic relationship to better understand B. pseudomallei persistence in P314’s sputum despite intensive initial therapeutic regimens. We collected and sequenced 118 B. pseudomallei isolates from P314’s airways over a >16-year period, and also sampled the patient’s home environment, recovering six closely related B. pseudomallei isolates from the household water system. Using comparative genomics, we identified 126 SNPs in the core genome of the 124 isolates or 162 SNPs/indels when the accessory genome was included. The core SNPs were used to construct a phylogenetic tree, which demonstrated a close relationship between environmental and clinical isolates and detailed within-host evolutionary patterns. The phylogeny had little homoplasy, consistent with a strictly clonal mode of genetic inheritance. Repeated sampling revealed evidence of genetic diversification, but frequent extinctions left only one successful lineage through the first four years and two lineages after that. Overall, the evolution of this population is nonadaptive and best explained by genetic drift. However, some genetic and phenotypic changes are consistent with in situ adaptation. Using a mouse model, P314 isolates caused greatly reduced morbidity and mortality compared to the environmental isolates. Additionally, potentially adaptive phenotypes emerged and included differences in the O-antigen, capsular polysaccharide, motility, and colony morphology. The >13-year co-existence of two long-lived lineages presents interesting hypotheses that can be tested in future studies to provide additional insights into selective pressures, niche differentiation, and microbial adaptation. This unusual melioidosis case presents a rare example of the evolutionary progression towards commensalism by a highly virulent pathogen within a single human host.

Pathogens frequently jump between different hosts, and associated adaptation may lead to the emergence of new infectious agents. Such host-jumping evolution is witnessed through endpoint analyses but these cannot capture genetic changes in lineages that have gone extinct. In this study, we have identified and monitored an example of the evolution of a bacterium often deadly to its mammalian host, in an unprecedented case whereby disease lessened through time and the pathogen became a part of the commensal human flora. We used genomic analyses to characterize more than 16 years of this evolutionary process and the stepwise mutations that control pathogen interactions with the patient. Soon after infection, mutational changes occurred that allowed the bacterium to remain in the airways without causing disease. This shift towards avirulence was determined based on clinical data and virulence testing in an animal model. In addition, mutations occurred that contributed to the persistence of the bacteria in the patient's lungs. Finally, we found evidence for the evolutionary emergence and persistence of two distinct lineages of the bacterium over the last 13 years, presenting interesting questions about niche utilization. Bacteria are ubiquitous in the human body and almost all are beneficial or benign. In this study, we document the evolutionary conversion of a normally deadly bacterium towards a commensal.

Epidemiological and comparative genomic analysis of Bacillus anthracis isolated from northern Vietnam

To understand the epidemiological and genetic background of anthrax cases occurring in Vietnam from 2011 to 2015, we surveilled and genetically analyzed Bacillus anthracis isolated in the north of the country. Epidemiological surveillance showed that most human cutaneous anthrax cases occurred in association with animal dissection. Whole-genome sequences were obtained from six B. anthracis strains from human patients with cutaneous anthrax in the endemic area. Comparative genomic analysis showed that the genetic homogeneity among Vietnamese B. anthracis strains was very high. All Vietnamese B. anthracis strains belonged to the canSNP lineage of A.Br.011/009, which mostly consists of strains of the trans-Eurasian (TEA) group, including the most closely related strain, Carbosap. To clarify the genetic diversity of Vietnamese strains and strains belonging to A.Br.011/009 and A.Br.008/011 canSNP lineages, we applied a reference genome-based single-nucleotide polymorphism (SNP) and gene-by-gene genomic analysis (whole-genome MLST) strategy. The phylogeny from core genome SNPs revealed that the Vietnamese strains were positioned close to each other; moreover, several SNPs specific to Vietnamese B. anthracis were identified. Whole-genome MLST analysis revealed the differences in the number of SNPs between Vietnamese strains, which could enable discrimination at the strain level.

Genomic epizootiology of a Brucella abortus outbreak in Northern Ireland (1997-2012)

BACKGROUND

In the recent past (1997-2012), Northern Ireland in the United Kingdom suffered an outbreak of Brucella abortus, which at its height affected over 200 cattle herds. Initially, isolates were characterized using multi-locus variable number tandem repeats analysis (MLVA). While informative in this setting, hyper-variability in some loci limited the resolution necessary to infer fine-scale disease transmission networks. Consequently, we applied whole-genome sequencing to isolates from this outbreak to evaluate higher resolution markers for disease epizootiology.

RESULTS

Phylogenetic analysis revealed that the B. abortus outbreak in Northern Ireland was caused by two distinct pathogen lineages. One contained isolates consistent with the 1997-2012 outbreak being linked to a previous endemic infection thought eradicated. The dominant second lineage exhibited little genetic diversity throughout the recrudescent outbreak, with limited population sub-structure evident. This finding was inconsistent with prior MLVA molecular characterizations that suggested the presence of seven clonal complexes. Spatio-temporal modeling revealed a significant association of pairwise SNP differences between isolates and geographic distances. However, effect sizes were very small due to reduced pathogen diversity.

CONCLUSIONS

Genome sequence data suggested that hyper-variability in some MLVA loci contributed to an overestimate of pathogen diversity in the most recent outbreak. The low diversity observed in our genomic dataset made it inappropriate to apply phylodynamic methods to these data. We conclude that maintaining data repositories of genome sequence data will be invaluable for source attribution/epizootiological inference should recrudescence ever re-occur. However genomic epizootiological methods may have limited utility in some settings, such as when applied to recrudescent/re-emergent infections of slowly-evolving bacterial pathogens.

Phylogenomic structure of Bacillus anthracis isolates in the Northern Cape Province, South Africa revealed novel single nucleotide polymorphisms

Bacillus anthracis, the aetiological agent of anthrax, is regarded as a highly monomorphic pathogen that presents a low genetic diversity using standard molecular techniques. Whole genome sequencing and single nucleotide polymorphisms (SNPs) are definitive signatures for subtyping of B. anthracis. Here we employed whole genome single nucleotide polymorphism (wgSNP) analysis to investigate the genetic diversity of B. anthracis in the historically endemic region of Northern Cape Province (NCP), South Africa. Twenty-six isolates from anthrax outbreaks that occurred between 1998 and 2008/9 in NCP as well as from Namibia-South Africa Transfontier Conservation area and Botswana were compared to global B. anthracis genomes. Most NCP B. anthracis strains (n = 22) clustered in the A.Br.003/004 (A.Br.101) branch and are closely related to the Zimbabwe and Mozambique strains (A.Br.102 branch). A total of 4923 parsimony informative-SNPs accurately established the A.Br.003/004 phylogenetic relationships of the NCP isolates into two distinct sub-clades and SNP markers designated as A.Br.172 and A.Br.173 were developed. Other NCP strains (n = 2) grouped in the A.Br.001/002 (Sterne) branch while strains (n = 2) from the Namibia-South Africa Transfontier Conservation area and Botswana clustered in A.Br.005/006 (Ancient A) branch. The sequenced B. anthracis strains (A0094, A0096 and A0097) that clustered in the A.Br.064 (V770) clade were isolated from Vaalbos National Park and similar strains have not been isolated. The B. anthracis A0088 strain cluster with the NCP strains in the A.Br.003/004 (A.Br.172) SNP branch which has been isolated in NCP, South Africa. This study highlights the phylogenetic structure of NCP B. anthracis strains with distinctive SNP branches important for forensic tracing and novel SNP discovery purposes. The sequenced strains will serve as a means to further trace the dissemination of B. anthracis outbreaks in NCP, South Africa, and on the continent, as well as for forensic tracking on a global scale.

Phylogeography and epidemiology of Brucella suis biovar 2 in wildlife and domestic swine

Swine brucellosis due to Brucella suis biovar 2 (bv2) is enzootic in wild boar and hare in continental Europe and may cause major economic losses to the pig industry, mainly in free-ranged pig farms. The high nucleotide identity found among the B. suis biovar 2 isolates has long hindered the full understanding of the epidemiology and the phylogeography of the disease. Here, we used multilocus variable-number tandem-repeat (VNTR) analysis (MLVA) and whole-genome analysis to identify single-nucleotide polymorphisms (SNPs) in order to gain insights from the largest B. suis bv2 dataset analyzed so far composed of domestic pigs and wildlife isolates collected throughout Europe since the 1970s. We found four major clades with a specific phylogeographic pattern. The Iberian clade contains isolates exclusively from the Iberian Peninsula. The Central European clade includes most isolates from France, Northern Italy, Switzerland and an important proportion of those of Northern Spain. The Eastern European clade clustered isolates from Croatia and Hungary mainly but also from areas of France, Germany, Italy and Poland. Finally, a separated Sardinian clade grouped three isolates from this island. At fine scale, MLVA demonstrated an endemic status of the infection in Europe and it allowed tracking a large outbreak formed by different farms from Spain linked to the same infection source. The whole genome SNP analysis showed that the strains form genetically distinct clades, shared between wild boar and pigs, in agreement with the MLVA clades. Interestingly, all hare isolates clustered together within two groups composed exclusively of wildlife isolates. Our results support the hypothesis that maintenance and spread of B. suis bv2 in Europe is a dynamic process linked to the natural expansion of wild boar as the main wild reservoir of the infection, while spread over long distances is found largely dependent on anthropogenic activities.

Comparative genomics and phylogenetic analysis of Bacillus anthracis strains isolated from domestic animals in Japan

Anthrax, caused by Bacillus anthracis, is a severe zoonosis with a great impact on both human and animal health. In the present study, we identified the phylogenetic relationships among 16 Japanese strains of B. anthracis, including eight bovine strains, two equine strains, five swine strains, and one former vaccine strain, using in silico canonical single nucleotide polymorphism (canSNP) and core genome SNP analyses. The results of our in silico canSNP analysis suggest that these 16 Japanese strains could be divided into four lineages: i) one equine strain in A.Br.Ames, ii) one equine and six bovine strains in A.Br.001/002, iii) five swine and one bovine strain in A.Br.Aust94, and iv) one bovine and one vaccine strain in A.Br.008/011. A comparison with non-Japanese B. anthracis strains revealed a total of 3787 SNPs identified from the whole genome sequences of the Japanese strains; these SNP data were subjected to a phylogenetic analysis using the maximum parsimony (MP) method. Our core genome SNP analysis was also able to detect differences of a few chromosomal SNPs across clonal strains from the same cases that had different storage and passage histories. Additionally, our whole genome SNP analysis clearly indicated that the Japanese swine anthrax cases of 1982 were caused by at least three independent strains; however, their phylogeny revealed no clear relationship with swine strains from other countries. The bovine strain belonging to the A.Br.008/011 lineage differed from a former Japanese vaccine strain by only 12 SNPs. Together with the phylogenic results and epidemiological circumstances, the diversity of strains reveals that the B. anthracis available in Japan probably resulted from multiple relatively recent import events, rather than reflecting the persistence of a more ancient ecologically established group.

Microbial sequence typing in the genomic era

Next-generation sequencing (NGS), also known as high-throughput sequencing, is changing the field of microbial genomics research. NGS allows for a more comprehensive analysis of the diversity, structure and composition of microbial genes and genomes compared to the traditional automated Sanger capillary sequencing at a lower cost. NGS strategies have expanded the versatility of standard and widely used typing approaches based on nucleotide variation in several hundred DNA sequences and a few gene fragments (MLST, MLVA, rMLST and cgMLST). NGS can now accommodate variation in thousands or millions of sequences from selected amplicons to full genomes (WGS, NGMLST and HiMLST). To extract signals from high-dimensional NGS data and make valid statistical inferences, novel analytic and statistical techniques are needed. In this review, we describe standard and new approaches for microbial sequence typing at gene and genome levels and guidelines for subsequent analysis, including methods and computational frameworks. We also present several applications of these approaches to some disciplines, namely genotyping, phylogenetics and molecular epidemiology.

Comparative pan-genomic analyses of Orientia tsutsugamushi reveal an exceptional model of bacterial evolution driving genomic diversity

Orientia tsutsugamushi, formerly Rickettsia tsutsugamushi, is an obligate intracellular pathogen that causes scrub typhus, an underdiagnosed acute febrile disease with high morbidity. Scrub typhus is transmitted by the larval stage (chigger) of Leptotrombidium mites and is irregularly distributed across endemic regions of Asia, Australia and islands of the western Pacific Ocean. Previous work to understand population genetics in O. tsutsugamushi has been based on sub-genomic sampling methods and whole-genome characterization of two genomes. In this study, we compared 40 genomes from geographically dispersed areas and confirmed patterns of extensive homologous recombination likely driven by transposons, conjugative elements and repetitive sequences. High rates of lateral gene transfer (LGT) among O. tsutsugamushi genomes appear to have effectively eliminated a detectable clonal frame, but not our ability to infer evolutionary relationships and phylogeographical clustering. Pan-genomic comparisons using 31 082 high-quality bacterial genomes from 253 species suggests that genomic duplication in O. tsutsugamushi is almost unparalleled. Unlike other highly recombinant species where the uptake of exogenous DNA largely drives genomic diversity, the pan-genome of O. tsutsugamushi is driven by duplication and divergence. Extensive gene innovation by duplication is most commonly attributed to plants and animals and, in contrast with LGT, is thought to be only a minor evolutionary mechanism for bacteria. The near unprecedented evolutionary characteristics of O. tsutsugamushi, coupled with extensive intra-specific LGT, expand our present understanding of rapid bacterial evolutionary adaptive mechanisms.

The reporting of a Bacillus anthracis B-clade strain in South Africa after more than 20 years

OBJECTIVES

Anthrax is a disease with an age old history in Africa caused by the Gram-positive endospore forming soil bacterium Bacillus anthracis. Epizootics of wild ungulates occur annually in the enzootic region of Pafuri, Kruger National Park (KNP) in the Limpopo Province of South Africa. Rigorous routine surveillance and diagnostics in KNP, has not revealed these rare isolates since the 1990s, despite unabated annual outbreaks. In 2011 a cheetah was diagnosed as anthrax positive from a private game reserve in Limpopo Province and reported to State Veterinary Services for further investigation. Isolation, molecular diagnostics, whole genome sequencing and comparative genomics were carried out for B. anthracis KC2011.

RESULTS

Bacteriological and molecular diagnostics confirmed the isolate as B. anthracis. Subsequent typing and whole genome single nucleotide polymorphisms analysis indicated it clustered alongside B. anthracis SA A0091 in the B.Br.010 SNP branch. Unlike B. anthracis KrugerB strain, KC2011 strain has unique SNPs and represents a new branch in the B-clade. The isolation and genotypic characterisation of KC2011 demonstrates a gap in the reporting of anthrax outbreaks in the greater Limpopo province area. The identification of vulnerable and susceptible cheetah mortalities due to this strain has implications for conservation measures and disease control.

Phylogenomics

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. The abundance of genomic data for an enormous variety of organisms has enabled phylogenomic inference of many groups, and this has motivated the development of many computer programs implementing the associated methods. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Phylogenetics of a Fungal Invasion: Origins and Widespread Dispersal of White-Nose Syndrome

Globalization has facilitated the worldwide movement and introduction of pathogens, but epizoological reconstructions of these invasions are often hindered by limited sampling and insufficient genetic resolution among isolates. Pseudogymnoascus destructans, a fungal pathogen causing the epizootic of white-nose syndrome in North American bats, has exhibited few genetic polymorphisms in previous studies, presenting challenges for both epizoological tracking of the spread of this fungus and for determining its evolutionary history. We used single nucleotide polymorphisms (SNPs) from whole-genome sequencing and microsatellites to construct high-resolution phylogenies of P. destructans. Shallow genetic diversity and the lack of geographic structuring among North American isolates support a recent introduction followed by expansion via clonal reproduction across the epizootic zone. Moreover, the genetic relationships of isolates within North America suggest widespread mixing and long-distance movement of the fungus. Genetic diversity among isolates of P. destructans from Europe was substantially higher than in those from North America. However, genetic distance between the North American isolates and any given European isolate was similar to the distance between the individual European isolates. In contrast, the isolates we examined from Asia were highly divergent from both European and North American isolates. Although the definitive source for introduction of the North American population has not been conclusively identified, our data support the origin of the North American invasion by P. destructans from Europe rather than Asia.

This phylogenetic study of the bat white-nose syndrome agent, P. destructans, uses genomics to elucidate evolutionary relationships among populations of the fungal pathogen to understand the epizoology of this biological invasion. We analyze hypervariable and abundant genetic characters (microsatellites and genomic SNPs, respectively) to reveal previously uncharacterized diversity among populations of the pathogen from North America and Eurasia. We present new evidence supporting recent introduction of the fungus to North America from a diverse Eurasian population, with limited increase in genetic variation in North America since that introduction.

Knowledge gaps that hamper prevention and control of Mycobacterium avium subspecies paratuberculosis infection

In the last decades, many regional and country-wide control programmes for Johne's disease (JD) were developed due to associated economic losses, or because of a possible association with Crohn's disease. These control programmes were often not successful, partly because management protocols were not followed, including the introduction of infected replacement cattle, because tests to identify infected animals were unreliable, and uptake by farmers was not high enough because of a perceived low return on investment. In the absence of a cure or effective commercial vaccines, control of JD is currently primarily based on herd management strategies to avoid infection of cattle and restrict within-farm and farm-to-farm transmission. Although JD control programmes have been implemented in most developed countries, lessons learned from JD prevention and control programmes are underreported. Also, JD control programmes are typically evaluated in a limited number of herds and the duration of the study is less than 5 year, making it difficult to adequately assess the efficacy of control programmes. In this manuscript, we identify the most important gaps in knowledge hampering JD prevention and control programmes, including vaccination and diagnostics. Secondly, we discuss directions that research should take to address those knowledge gaps.

Advances in Understanding Bacterial Pathogenesis Gained from Whole-Genome Sequencing and Phylogenetics

The development of next-generation sequencing as a cost-effective technology has facilitated the analysis of bacterial population structure at a whole-genome level and at scale. From these data, phylogenic trees have been constructed that define population structures at a local, national, and global level, providing a framework for genetic analysis. Although still at an early stage, these approaches have yielded progress in several areas, including pathogen transmission mapping, the genetics of niche colonization and host adaptation, as well as gene-to-phenotype association studies. Antibiotic resistance has proven to be a major challenge in the early 21(st) century, and phylogenetic analyses have uncovered the dramatic effect that the use of antibiotics has had on shaping bacterial population structures. An update on insights into bacterial evolution from comparative genomics is provided in this review.

Temporal phylogeography of Yersinia pestis in Madagascar: Insights into the long-term maintenance of plague

BACKGROUND

Yersinia pestis appears to be maintained in multiple, geographically separate, and phylogenetically distinct subpopulations within the highlands of Madagascar. However, the dynamics of these locally differentiated subpopulations through time are mostly unknown. To address that gap and further inform our understanding of plague epidemiology, we investigated the phylogeography of Y. pestis in Madagascar over an 18 year period.

METHODOLOGY/PRINCIPAL FINDINGS

We generated whole genome sequences for 31 strains and discovered new SNPs that we used in conjunction with previously identified SNPs and variable-number tandem repeats (VNTRs) to genotype 773 Malagasy Y. pestis samples from 1995 to 2012. We mapped the locations where samples were obtained on a fine geographic scale to examine phylogeographic patterns through time. We identified 18 geographically separate and phylogenetically distinct subpopulations that display spatial and temporal stability, persisting in the same locations over a period of almost two decades. We found that geographic areas with higher levels of topographical relief are associated with greater levels of phylogenetic diversity and that sampling frequency can vary considerably among subpopulations and from year to year. We also found evidence of various Y. pestis dispersal events, including over long distances, but no evidence that any dispersal events resulted in successful establishment of a transferred genotype in a new location during the examined time period.

CONCLUSIONS/SIGNIFICANCE

Our analysis suggests that persistent endemic cycles of Y. pestis transmission within local areas are responsible for the long term maintenance of plague in Madagascar, rather than repeated episodes of wide scale epidemic spread. Landscape likely plays a role in maintaining Y. pestis subpopulations in Madagascar, with increased topographical relief associated with increased levels of localized differentiation. Local ecological factors likely affect the dynamics of individual subpopulations and the associated likelihood of observing human plague cases in a given year in a particular location.

A single amino acid polymorphism in the glycosyltransferase CpsK defines four Streptococcus suis serotypes

The capsular polysaccharide (CPS) is the major virulence factor of the emerging zoonotic pathogen Streptococcus suis. CPS differences are also the basis for serological differentiation of the species into 29 serotypes. Serotypes 2 and 1/2, which possess identical gene content in their cps loci, express CPSs that differ only by substitution of galactose (Gal) by N-acetylgalactosamine (GalNAc) in the CPS side chain. The same sugar substitution differentiates the CPS of serotypes 14 and 1, whose cps loci are also identical in gene content. Here, using mutagenesis, CPS structural analysis, and protein structure modeling, we report that a single amino acid polymorphism in the glycosyltransferase CpsK defines the enzyme substrate predilection for Gal or GalNAc and therefore determines CPS composition, structure, and strain serotype. We also show that the different CPS structures have similar antiphagocytic properties and that serotype switching has limited impact on the virulence of S. suis.

Origins and global context of Brucella abortus in Italy

Background

Brucellosis is a common and chronic disease of cattle and other bovids that often causes reproductive disorders. Natural infection in cattle is caused by Brucella abortus and transmission typically occurs during abortions, calving, or nursing. Brucellosis is also a major zoonotic disease due to contamination of dairy products or contact with the tissues of infected animals. Brucellosis has been eradicated from most of the developed world in the last 40 years but persists in many regions—the disease remains prevalent in portions of Africa, the Middle East, Asia, and Central and South America, as well as in the Mediterranean basin. In Italy, B. abortus has persisted in southern regions in both cattle and water buffalo. Previous attempts at analyzing the phylogenetics of B. abortus in Italy have been challenging due to limited genetic variability and unresolved global population genetic structure of this pathogen.

Results

We conducted genome-wide phylogenetic analyses on 11 representative strains of B. abortus from Italy, and compared these sequences to a worldwide collection of publically available genomes. Italian isolates belong to three clades that are basal to the main and global B. abortus lineage. Using six SNP-based assays designed to identify substructure within the Italian clades, we surveyed a collection of 261 isolates and found that one clade predominates throughout endemic districts in the country, while the other two clades are more geographically restricted to portions of southern Italy.

Conclusions

Although related strains exist worldwide, B. abortus isolates from Italy are substantially different than those found in much of the rest of Europe and North America, and are more closely related to strains from the Middle East and Asia. Our assays targeting genetic substructure within Italy allowed us to identify the major lineages quickly and inexpensively, without having to generate whole genome sequences for a large isolate collection. These findings highlight the importance of genetic studies to assess the status and the history of pathogens.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-017-0939-0) contains supplementary material, which is available to authorized users.

Is Predominant Clonal Evolution a Common Evolutionary Adaptation to Parasitism in Pathogenic Parasitic Protozoa, Fungi, Bacteria, and Viruses?

We propose that predominant clonal evolution (PCE) in microbial pathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure. The main features of PCE are (1) strong linkage disequilibrium, (2) the widespread occurrence of stable genetic clusters blurred by occasional bouts of genetic exchange ('near-clades'), (3) the existence of a "clonality threshold", beyond which recombination is efficiently countered by PCE, and near-clades irreversibly diverge. We hypothesize that the PCE features are not mainly due to natural selection but also chiefly originate from in-built genetic properties of pathogens. We show that the PCE model obtains even in microbes that have been considered as 'highly recombining', such as Neisseria meningitidis, and that some clonality features are observed even in Plasmodium, which has been long described as panmictic. Lastly, we provide evidence that PCE features are also observed in viruses, taking into account their extremely fast genetic turnover. The PCE model provides a convenient population genetic framework for any kind of micropathogen. It makes it possible to describe convenient units of analysis (clones and near-clades) for all applied studies. Due to PCE features, these units of analysis are stable in space and time, and clearly delimited. The PCE model opens up the possibility of revisiting the problem of species definition in these organisms. We hypothesize that PCE constitutes a major evolutionary strategy for protozoa, fungi, bacteria, and viruses to adapt to parasitism.

A Bacillus anthracis Genome Sequence from the Sverdlovsk 1979 Autopsy Specimens

Anthrax is a zoonotic disease that occurs naturally in wild and domestic animals but has been used by both state-sponsored programs and terrorists as a biological weapon. A Soviet industrial production facility in Sverdlovsk, USSR, proved deficient in 1979 when a plume of spores was accidentally released and resulted in one of the largest known human anthrax outbreaks. In order to understand this outbreak and others, we generated a Bacillus anthracis population genetic database based upon whole-genome analysis to identify all single-nucleotide polymorphisms (SNPs) across a reference genome. Phylogenetic analysis has defined three major clades (A, B, and C), B and C being relatively rare compared to A. The A clade has numerous subclades, including a major polytomy named the trans-Eurasian (TEA) group. The TEA radiation is a dominant evolutionary feature of B. anthracis, with many contemporary populations having resulted from a large spatial dispersal of spores from a single source. Two autopsy specimens from the Sverdlovsk outbreak were deep sequenced to produce draft B. anthracis genomes. This allowed the phylogenetic placement of the Sverdlovsk strain into a clade with two Asian live vaccine strains, including the Russian Tsiankovskii strain. The genome was examined for evidence of drug resistance manipulation or other genetic engineering, but none was found. The Soviet Sverdlovsk strain genome is consistent with a wild-type strain from Russia that had no evidence of genetic manipulation during its industrial production. This work provides insights into the world’s largest biological weapons program and provides an extensive B. anthracis phylogenetic reference.

The 1979 Russian anthrax outbreak resulted from an industrial accident at the Soviet anthrax spore production facility in the city of Sverdlovsk. Deep genomic sequencing of two autopsy specimens generated a draft genome and phylogenetic placement of the Soviet Sverdlovsk anthrax strain. While it is known that Soviet scientists had genetically manipulated Bacillus anthracis with the potential to evade vaccine prophylaxis and antibiotic therapeutics, there was no genomic evidence of this from the Sverdlovsk production strain genome. The whole-genome SNP genotype of the Sverdlovsk strain was used to precisely identify it and its close relatives in the context of an extensive global B. anthracis strain collection. This genomic identity can now be used for forensic tracking of this weapons material on a global scale and for future anthrax investigations.

Massive dispersal of Coxiella burnetii among cattle across the United States

Q-fever is an underreported disease caused by the bacterium Coxiella burnetii, which is highly infectious and has the ability to disperse great distances. It is a completely clonal pathogen with low genetic diversity and requires whole-genome analysis to identify discriminating features among closely related isolates. C. burnetii, and in particular one genotype (ST20), is commonly found in cow's milk across the entire dairy industry of the USA. This single genotype dominance is suggestive of host-specific adaptation, rapid dispersal and persistence within cattle. We used a comparative genomic approach to identify SNPs for high-resolution and high-throughput genotyping assays to better describe the dispersal of ST20 across the USA. We genotyped 507 ST20 cow milk samples and discovered three subgenotypes, all of which were present across the entire country and over the complete time period studied. Only one of these sub-genotypes was observed in a single dairy herd. The temporal and geographic distribution of these sub-genotypes is consistent with a model of large-scale, rapid, frequent and continuous dissemination on a continental scale. The distribution of subgenotypes is not consistent with wind-based dispersal alone, and it is likely that animal husbandry and transportation practices, including pooling of milk from multiple herds, have also shaped the patterns. On the scale of an entire country, there appear to be few barriers to rapid, frequent and large-scale dissemination of the ST20 subgenotypes.

Gall-ID: tools for genotyping gall-causing phytopathogenic bacteria

Understanding the population structure and genetic diversity of plant pathogens, as well as the effect of agricultural practices on pathogen evolution, is important for disease management. Developments in molecular methods have contributed to increase the resolution for accurate pathogen identification, but those based on analysis of DNA sequences can be less straightforward to use. To address this, we developed Gall-ID, a web-based platform that uses DNA sequence information from 16S rDNA, multilocus sequence analysis and whole genome sequences to group disease-associated bacteria to their taxonomic units. Gall-ID was developed with a particular focus on gall-forming bacteria belonging to Agrobacterium, Pseudomonas savastanoi, Pantoea agglomerans, and Rhodococcus. Members of these groups of bacteria cause growth deformation of plants, and some are capable of infecting many species of field, orchard, and nursery crops. Gall-ID also enables the use of high-throughput sequencing reads to search for evidence for homologs of characterized virulence genes, and provides downloadable software pipelines for automating multilocus sequence analysis, analyzing genome sequences for average nucleotide identity, and constructing core genome phylogenies. Lastly, additional databases were included in Gall-ID to help determine the identity of other plant pathogenic bacteria that may be in microbial communities associated with galls or causative agents in other diseased tissues of plants. The URL for Gall-ID is http://gall-id.cgrb.oregonstate.edu/.

Genome Sequences of 11 Brucella abortus Isolates from Persistently Infected Italian Regions

Bovine brucellosis, typically caused by Brucella abortus, has been eradicated from much of the developed world. However, the disease remains prevalent in southern Italy, persisting as a public and livestock health concern. We report here the whole-genome sequences of 11 isolates from cattle (Bos taurus) and water buffalo (Bubalus bubalis) that are representative of the current genetic diversity of B. abortus lineages circulating in Italy.

Whole-genome sequencing investigation of animal-skin-drum-associated UK anthrax cases reveals evidence of mixed populations and relatedness to a US case

There have been two anthrax cases affecting people that played and/or made animal-skin drums in the UK during the last 10 years, with single fatal occurrences in Scotland in 2006 and London in 2008. Investigations by the Health Protection Agency (now Public Health England) employing multi-locus-variable number tandem repeat analysis had previously linked the clinical cases to spores associated with animal skins and drums the patients had been in contact with. In this study, whole-genome sequencing of 23 Bacillus anthracis isolates harvested during the investigations was performed. High-quality draft assemblies of these genomes provided greater characterization of the B. anthracis strains present and placed them all upon a new branch of the global phylogeny. Although closely related, the clinical isolates from the two events, and another isolated from a drum-skin-associated case in New York in 2006, were distinct from each other. Multiple distinct genotypes were found during both investigations, implying either multiple contamination events or a single heterogeneous contamination. One environmental isolate from the Scottish incident was more closely related to London isolates than to the other Scottish isolates. As B. anthracis of this subgroup was present at both geographically and temporally distinct events, it may be more widespread at the source of contamination. All isolates were distinct from currently characterized West African strains, despite this being the likely origin of the drums and hides, therefore adding to our knowledge of B. anthracis diversity in the region.

The population genetics of Trypanosoma cruzi revisited in the light of the predominant clonal evolution model

Comparing the population structure of Trypanosoma cruzi with that of other pathogens, including parasitic protozoa, fungi, bacteria and viruses, shows that the agent of Chagas disease shares typical traits with many other species, related to a predominant clonal evolution (PCE) pattern: statistically significant linkage disequilibrium, overrepresented multilocus genotypes, near-clades (genetic subdivisions somewhat blurred by occasional genetic exchange/hybridization) and "Russian doll" patterns (PCE is observed, not only at the level of the whole species, but also, within the near-clades). Moreover, T. cruzi population structure exhibits linkage with the diversity of several strongly selected genes, with gene expression profiles, and with some major phenotypic traits. We discuss the evolutionary significance of these results, and their implications in terms of applied research (molecular epidemiology/strain typing, analysis of genes of interest, vaccine and drug design, immunological diagnosis) and of experimental evolution. Lastly, we revisit the long-term debate of describing new species within the T. cruzi taxon.

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.

Whole Genome Analysis of Injectional Anthrax Identifies Two Disease Clusters Spanning More Than 13 Years

Background

Anthrax is a rare disease in humans but elicits great public fear because of its past use as an agent of bioterrorism. Injectional anthrax has been occurring sporadically for more than ten years in heroin consumers across multiple European countries and this outbreak has been difficult to trace back to a source.

Methods

We took a molecular epidemiological approach in understanding this disease outbreak, including whole genome sequencing of Bacillus anthracis isolates from the anthrax victims. We also screened two large strain repositories for closely related strains to provide context to the outbreak.

Findings

Analyzing 60 Bacillus anthracis isolates associated with injectional anthrax cases and closely related reference strains, we identified 1071 Single Nucleotide Polymorphisms (SNPs). The synapomorphic SNPs (350) were used to reconstruct phylogenetic relationships, infer likely epidemiological sources and explore the dynamics of evolving pathogen populations. Injectional anthrax genomes separated into two tight clusters: one group was exclusively associated with the 2009–10 outbreak and located primarily in Scotland, whereas the second comprised more recent (2012–13) cases but also a single Norwegian case from 2000.

Interpretation

Genome-based differentiation of injectional anthrax isolates argues for at least two separate disease events spanning > 12 years. The genomic similarity of the two clusters makes it likely that they are caused by separate contamination events originating from the same geographic region and perhaps the same site of drug manufacturing or processing. Pathogen diversity within single patients challenges assumptions concerning population dynamics of infecting B. anthracis and host defensive barriers for injectional anthrax.

Funding

This work was supported by the United States Department of Homeland Security grant no. HSHQDC-10-C-00,139 and via a binational cooperative agreement between the United States Government and the Government of Germany. This work was supported by funds from the German Ministry of Defense (Sonderforschungsprojekt 25Z1-S-431,214). Support for sequencing was also obtained from Illumina, Inc. These sources had no role in the data generation or interpretation, and had not role in the manuscript preparation.

Panel 1: Research in Context Systematic Review

We searched PubMed for any article published before Jun. 17, 2015, with the terms “Bacillus anthracis” and “heroin”, or “injectional anthrax”. Other than our previously published work (Price et al., 2012), we found no other relevant studies on elucidating the global phylogenetic relationships of B. anthracis strains associated with injectional anthrax caused by recreational heroin consumption of spore-contaminated drug. There were, however, publically available genome sequences of two strains involved (Price et al., 2012, Grunow et al., 2013) and the draft genome sequence of Bacillus anthracis UR-1, isolated from a German heroin user (Ruckert et al., 2012) with only limited information on the genotyping of closely related strains (Price et al., 2012, Grunow et al., 2013).

Lay Person Interpretation

Injectional anthrax has been plaguing heroin drug users across Europe for more than 10 years. In order to better understand this outbreak, we assessed genomic relationships of all available injectional anthrax strains from four countries spanning a > 12 year period. Very few differences were identified using genome-based analysis, but these differentiated the isolates into two distinct clusters. This strongly supports a hypothesis of at least two separate anthrax spore contamination events perhaps during the drug production processes. Identification of two events would not have been possible from standard epidemiological analysis. These comprehensive data will be invaluable for classifying future injectional anthrax isolates and for future geographic attribution.

•

Whole genome sequences of injectional anthrax B. anthracis isolates fall in two tight but distinct genomic clusters.

•

The distinct genomic clusters are consistent with two or more disease events that overlap in time and space.

•

Defining pathogen clusters will lead to better public health responses to difficult to track disease outbreaks.

Molecular Epidemiology of Brucella abortus in Northern Ireland-1991 to 2012

Background

Brucellosis is the most common bacterial zoonoses worldwide. Bovine brucellosis caused by Brucella abortus has far reaching animal health and economic impacts at both the local and national levels. Alongside traditional veterinary epidemiology, the use of molecular typing has recently been applied to inform on bacterial population structure and identify epidemiologically-linked cases of infection. Multi-locus variable number tandem repeat VNTR analysis (MLVA) was used to investigate the molecular epidemiology of a well-characterised Brucella abortus epidemic in Northern Ireland involving 387 herds between 1991 and 2012.

Results

MLVA identified 98 unique B. abortus genotypes from disclosing isolates in the 387 herds involved in the epidemic. Clustering algorithms revealed the relatedness of many of these genotypes. Combined with epidemiological information on chronology of infection and geographic location, these genotype data helped to identify 7 clonal complexes which underpinned the outbreak over the defined period. Hyper-variability of some VNTR loci both within herds and individual animals led to detection of multiple genotypes associated with single outbreaks. However with dense sampling, these genotypes could still be associated with specific clonal complexes thereby permitting inference of epidemiological links. MLVA- based epidemiological monitoring data were congruent with an independent classical veterinary epidemiology study carried out in the same territory.

Conclusions

MLVA is a useful tool in ongoing disease surveillance of B. abortus outbreaks, especially when combined with accurate epidemiological information on disease tracings, geographical clustering of cases and chronology of infection.

Best practices for evaluating single nucleotide variant calling methods for microbial genomics

Innovations in sequencing technologies have allowed biologists to make incredible advances in understanding biological systems. As experience grows, researchers increasingly recognize that analyzing the wealth of data provided by these new sequencing platforms requires careful attention to detail for robust results. Thus far, much of the scientific Communit’s focus for use in bacterial genomics has been on evaluating genome assembly algorithms and rigorously validating assembly program performance. Missing, however, is a focus on critical evaluation of variant callers for these genomes. Variant calling is essential for comparative genomics as it yields insights into nucleotide-level organismal differences. Variant calling is a multistep process with a host of potential error sources that may lead to incorrect variant calls. Identifying and resolving these incorrect calls is critical for bacterial genomics to advance. The goal of this review is to provide guidance on validating algorithms and pipelines used in variant calling for bacterial genomics. First, we will provide an overview of the variant calling procedures and the potential sources of error associated with the methods. We will then identify appropriate datasets for use in evaluating algorithms and describe statistical methods for evaluating algorithm performance. As variant calling moves from basic research to the applied setting, standardized methods for performance evaluation and reporting are required; it is our hope that this review provides the groundwork for the development of these standards.

Diversification and Distribution of Ruminant Chlamydia abortus Clones Assessed by MLST and MLVA

Chlamydia abortus, an obligate intracellular bacterium, is the most common infectious cause of abortion in small ruminants worldwide and has zoonotic potential. We applied multilocus sequence typing (MLST) together with multiple-locus variable-number tandem repeat analysis (MLVA) to genotype 94 ruminant C. abortus strains, field isolates and samples collected from 1950 to 2011 in diverse geographic locations, with the aim of delineating C. abortus lineages and clones. MLST revealed the previously identified sequence types (STs) ST19, ST25, ST29 and ST30, plus ST86, a recently-assigned type on the Chlamydiales MLST website and ST87, a novel type harbouring the hemN_21 allele, whereas MLVA recognized seven types (MT1 to MT7). Minimum-spanning-tree analysis suggested that all STs but one (ST30) belonged to a single clonal complex, possibly reflecting the short evolutionary timescale over which the predicted ancestor (ST19) has diversified into three single-locus variants (ST86, ST87 and ST29) and further, through ST86 diversification, into one double-locus variant (ST25). ST descendants have probably arisen through a point mutation evolution mode. Interestingly, MLVA showed that in the ST19 population there was a greater genetic diversity than in other STs, most of which exhibited the same MT over time and geographical distribution. However, the evolutionary pathways of C. abortus STs seem to be diverse across geographic distances with individual STs restricted to particular geographic locations. The ST30 singleton clone displaying geographic specificity and represented by the Greek strains LLG and POS was effectively distinguished from the clonal complex lineage, supporting the notion that possibly two separate host adaptations and hence independent bottlenecks of C. abortus have occurred through time. The combination of MLST and MLVA assays provides an additional level of C. abortus discrimination and may prove useful for the investigation and surveillance of emergent C. abortus clonal populations.

Using multilocus sequence typing to study bacterial variation: prospects in the genomic era

Multilocus sequence typing (MLST) indexes the sequence variation present in a small number (usually seven) of housekeeping gene fragments located around the bacterial genome. Unique alleles at these loci are assigned arbitrary integer identifiers, which effectively summarizes the variation present in several thousand base pairs of genome sequence information as a series of numbers. Comparing bacterial isolates using allele-based methods efficiently corrects for the effects of lateral gene transfer present in many bacterial populations and is computationally efficient. This 'gene-by-gene' approach can be applied to larger collections of loci, such as the ribosomal protein genes used in ribosomal MLST (rMLST), up to and including the complete set of coding sequences present in a genome, whole-genome MLST (wgMLST), providing scalable, efficient and readily interpreted genome analysis.

Limitations of variable number of tandem repeat typing identified through whole genome sequencing of Mycobacterium avium subsp. paratuberculosis on a national and herd level

Background

Mycobacterium avium subsp. paratuberculosis (MAP), the causative bacterium of Johne’s disease in dairy cattle, is widespread in the Canadian dairy industry and has significant economic and animal welfare implications. An understanding of the population dynamics of MAP can be used to identify introduction events, improve control efforts and target transmission pathways, although this requires an adequate understanding of MAP diversity and distribution between herds and across the country. Whole genome sequencing (WGS) offers a detailed assessment of the SNP-level diversity and genetic relationship of isolates, whereas several molecular typing techniques used to investigate the molecular epidemiology of MAP, such as variable number of tandem repeat (VNTR) typing, target relatively unstable repetitive elements in the genome that may be too unpredictable to draw accurate conclusions. The objective of this study was to evaluate the diversity of bovine MAP isolates in Canadian dairy herds using WGS and then determine if VNTR typing can distinguish truly related and unrelated isolates.

Results

Phylogenetic analysis based on 3,039 SNPs identified through WGS of 124 MAP isolates identified eight genetically distinct subtypes in dairy herds from seven Canadian provinces, with the dominant type including over 80% of MAP isolates. VNTR typing of 527 MAP isolates identified 12 types, including “bison type” isolates, from seven different herds. At a national level, MAP isolates differed from each other by 1–2 to 239–240 SNPs, regardless of whether they belonged to the same or different VNTR types. A herd-level analysis of MAP isolates demonstrated that VNTR typing may both over-estimate and under-estimate the relatedness of MAP isolates found within a single herd.

Conclusions

The presence of multiple MAP subtypes in Canada suggests multiple introductions into the country including what has now become one dominant type, an important finding for Johne’s disease control. VNTR typing often failed to identify closely and distantly related isolates, limiting the applicability of using this typing scheme to study the molecular epidemiology of MAP at a national and herd-level.

Genomic epidemiology of the Haitian cholera outbreak: a single introduction followed by rapid, extensive, and continued spread characterized the onset of the epidemic

For centuries, cholera has been one of the most feared diseases. The causative agent Vibrio cholerae is a waterborne Gram-negative enteric pathogen eliciting a severe watery diarrheal disease. In October 2010, the seventh pandemic reached Haiti, a country that had not experienced cholera for more than a century. By using whole-genome sequence typing and mapping strategies of 116 serotype O1 strains from global sources, including 44 Haitian genomes, we present a detailed reconstructed evolutionary history of the seventh pandemic with a focus on the Haitian outbreak. We catalogued subtle genomic alterations at the nucleotide level in the genome core and architectural rearrangements from whole-genome map comparisons. Isolates closely related to the Haitian isolates caused several recent outbreaks in southern Asia. This study provides evidence for a single-source introduction of cholera from Nepal into Haiti followed by rapid, extensive, and continued clonal expansion. The phylogeographic patterns in both southern Asia and Haiti argue for the rapid dissemination of V. cholerae across the landscape necessitating real-time surveillance efforts to complement the whole-genome epidemiological analysis. As eradication efforts move forward, phylogeographic knowledge will be important for identifying persistent sources and monitoring success at regional levels. The results of molecular and epidemiological analyses of this outbreak suggest that an indigenous Haitian source of V. cholerae is unlikely and that an indigenous source has not contributed to the genomic evolution of this clade.

In this genomic epidemiology study, we have applied high-resolution whole-genome-based sequence typing methodologies on a comprehensive set of genome sequences that have become available in the aftermath of the Haitian cholera epidemic. These sequence resources enabled us to reassess the degree of genomic heterogeneity within the Vibrio cholerae O1 serotype and to refine boundaries and evolutionary relationships. The established phylogenomic framework showed how outbreak isolates fit into the global phylogeographic patterns compared to a comprehensive globally and temporally diverse strain collection and provides strong molecular evidence that points to a nonindigenous source of the 2010 Haitian cholera outbreak and refines epidemiological standards used in outbreak investigations for outbreak inclusion/exclusion following the concept of genomic epidemiology. The generated phylogenomic data have major public health relevance in translating sequence-based information to assist in future diagnostic, epidemiological, surveillance, and forensic studies of cholera.

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 of Bacillus anthracis in the country of Georgia shows evidence of population structuring and is dissimilar to other regional genotypes

Sequence analyses and subtyping of Bacillus anthracis strains from Georgia reveal a single distinct lineage (Aust94) that is ecologically established. Phylogeographic analysis and comparisons to a global collection reveals a clade that is mostly restricted to Georgia. Within this clade, many groups are found around the country, however at least one subclade is only found in the eastern part. This pattern suggests that dispersal into and out of Georgia has been rare and despite historical dispersion within the country, for at least for one lineage, current spread is limited.

Genetic diversity of Bacillus anthracis in Europe: genotyping methods in forensic and epidemiologic investigations

Bacillus anthracis, the etiological agent of anthrax, a zoonosis relatively common throughout the world, can be used as an agent of bioterrorism. In naturally occurring outbreaks and in criminal release of this pathogen, a fast and accurate diagnosis is crucial to an effective response. Microbiological forensics and epidemiologic investigations increasingly rely on molecular markers, such as polymorphisms in DNA sequence, to obtain reliable information regarding the identification or source of a suspicious strain. Over the past decade, significant research efforts have been undertaken to develop genotyping methods with increased power to differentiate B. anthracis strains. A growing number of DNA signatures have been identified and used to survey B. anthracis diversity in nature, leading to rapid advances in our understanding of the global population of this pathogen. This article provides an overview of the different phylogenetic subgroups distributed across the world, with a particular focus on Europe. Updated information on the anthrax situation in Europe is reported. A brief description of some of the work in progress in the work package 5.1 of the AniBioThreat project is also presented, including (1) the development of a robust typing tool based on a suspension array technology and multiplexed single nucleotide polymorphisms scoring and (2) the typing of a collection of DNA from European isolates exchanged between the partners of the project. The know-how acquired will contribute to improving the EU's ability to react rapidly when the identity and real origin of a strain need to be established.

Genetic variations of live attenuated plague vaccine strains (Yersinia pestis EV76 lineage) during laboratory passages in different countries

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

Genotyping of Coxiella burnetii from domestic ruminants and human in Hungary: indication of various genotypes

BACKGROUND

Information about the genotypic characteristic of Coxiella burnetii from Hungary is lacking. The aim of this study is to describe the genetic diversity of C. burnetii in Hungary and compare genotypes with those found elsewhere. A total of 12 samples: (cattle, n = 6, sheep, n = 5 and human, n = 1) collected from across Hungary were studied by a 10-loci multispacer sequence typing (MST) and 6-loci multiple-locus variable-number of tandem repeat analysis (MLVA). Phylogenetic relationships among MST genotypes show how these Hungarian samples are related to others collected around the world.

RESULTS

Three MST genotypes were identified: sequence type (ST) 20 has also been identified in ruminants from other European countries and the USA, ST28 was previously identified in Kazakhstan, and the proposed ST37 is novel. All MST genotypes yielded different MLVA genotypes and three different MLVA genotypes were identified within ST20 samples alone. Two novel MLVA types 0-9-5-5-6-2 (AG) and 0-8-4-5-6-2 (AF) (Ms23-Ms24-Ms27-Ms28-Ms33-Ms34) were defined in the ovine materials correlated with ST28 and ST37. Samples from different parts of the phylogenetic tree were associated with different hosts, suggesting host-specific adaptations.

CONCLUSIONS

Even with the limited number of samples analysed, this study revealed high genetic diversity among C. burnetii in Hungary. Understanding the background genetic diversity will be essential in identifying and controlling outbreaks.

High-throughput sequencing of Bacillus anthracis in France: investigating genome diversity and population structure using whole-genome SNP discovery

Background

Single nucleotide polymorphisms (SNPs) are ideal signatures for subtyping monomorphic pathogens such as Bacillus anthracis. Here we report the use of next-generation sequencing technology to investigate the historical, geographic and genetic diversity of Bacillus anthracis in France. 122 strains isolated over a 60-years period throughout the country were whole-genome sequenced and comparative analyses were carried out with a focus on SNPs discovery to discriminate regional sub-groups of strains.

Results

A total of 1581 chromosomal SNPs precisely establish the phylogenetic relationships existing between the French strains. Phylogeography patterns within the three canSNP sub-lineages present in France (i.e. B.Br.CNEVA, A.Br.011/009 and A.Br.001/002) were observed. One of the more remarkable findings was the identification of a variety of genotypes within the A.Br.011/009 sub-group that are persisting in the different regions of France. The 560 SNPs defining the A.Br.011/009- affiliated French strains split the Trans-Eurasian sub-group into six distinct branches without any intermediate nodes. Distinct sub-branches, with some geographic clustering, were resolved. The 345 SNPs defining the major B.Br CNEVA sub-lineage clustered three main phylogeographic clades, the Alps, the Pyrenees, and the Massif Central, with a small Saône-et-Loire sub-cluster nested within the latter group. The French strains affiliated to the minor A.Br.001/002 group were characterized by 226 SNPs. All recent isolates collected from the Doubs department were closely related. Identification of SNPs from whole-genome sequences facilitates high-resolution strain tracking and provides the level of discrimination required for outbreak investigations. Eight diagnostic SNPs, representative of the main French-specific phylogeographic clusters, were therefore selected and developed into high-resolution melting SNP discriminative assays.

Conclusions

This work has established one of the most accurate phylogenetic reconstruction of B. anthracis population structure in a country. An extensive next-generation sequencing (NGS) dataset of 122 French strains have been created that allowed the identification of novel diagnostic SNPs useful to rapidly determine the geographic origin of any strain found in France.

The epidemic of extended-spectrum-β-lactamase-producing Escherichia coli ST131 is driven by a single highly pathogenic subclone, H30-Rx

The Escherichia coli sequence type 131 (ST131) clone is notorious for extraintestinal infections, fluoroquinolone resistance, and extended-spectrum beta-lactamase (ESBL) production, attributable to a CTX-M-15-encoding mobile element. Here, we applied pulsed-field gel electrophoresis (PFGE) and whole-genome sequencing to reconstruct the evolutionary history of the ST131 clone. PFGE-based cluster analyses suggested that both fluoroquinolone resistance and ESBL production had been acquired by multiple ST131 sublineages through independent genetic events. In contrast, the more robust whole-genome-sequence-based phylogenomic analysis revealed that fluoroquinolone resistance was confined almost entirely to a single, rapidly expanding ST131 subclone, designated H30-R. Strikingly, 91% of the CTX-M-15-producing isolates also belonged to a single, well-defined clade nested within H30-R, which was named H30-Rx due to its more extensive resistance. Despite its tight clonal relationship with H30Rx, the CTX-M-15 mobile element was inserted variably in plasmid and chromosomal locations within the H30-Rx genome. Screening of a large collection of recent clinical E. coli isolates both confirmed the global clonal expansion of H30-Rx and revealed its disproportionate association with sepsis (relative risk, 7.5; P < 0.001). Together, these results suggest that the high prevalence of CTX-M-15 production among ST131 isolates is due primarily to the expansion of a single, highly virulent subclone, H30-Rx.

We applied an advanced genomic approach to study the recent evolutionary history of one of the most important Escherichia coli strains in circulation today. This strain, called sequence type 131 (ST131), causes multidrug-resistant bladder, kidney, and bloodstream infections around the world. The rising prevalence of antibiotic resistance in E. coli is making these infections more difficult to treat and is leading to increased mortality. Past studies suggested that many different ST131 strains gained resistance to extended-spectrum cephalosporins independently. In contrast, our research indicates that most extended-spectrum-cephalosporin-resistant ST131 strains belong to a single highly pathogenic subclone, called H30-Rx. The clonal nature of H30-Rx may provide opportunities for vaccine or transmission prevention-based control strategies, which could gain importance as H30-Rx and other extraintestinal pathogenic E. coli subclones become resistant to our best antibiotics.

Four genotyping schemes for phylogenetic analysis of Pseudomonas aeruginosa: comparison of their congruence with multi-locus sequence typing

Several molecular typing schemes have been proposed to differentiate among isolates and clonal groups, and hence establish epidemiological or phylogenetic links. It has been widely accepted that multi-locus sequence typing (MLST) is the gold standard for phylogenetic typing/long-term epidemiological surveillance, but other recently described methods may be easier to carry out, especially in settings with limited access to DNA sequencing. Comparing the performance of such techniques to MLST is therefore of relevance. A study was therefore carried out with a collection of P. aeruginosa strains (n = 133) typed by four typing schemes: MLST, multiple-locus variable number tandem repeat analysis (MLVA), pulsed-field gel electrophoresis (PFGE) and the commercial DiversiLab microbial typing system (DL). The aim of this study was to compare the results of each typing method with MLST. The Simpson's indices of diversity were 0.989, 0.980, 0.961 and 0.906 respectively for PFGE, MLVA, DL and MLST. The congruence between techniques was measured by the adjusted Wallace index (W): this coefficient indicates the probability that a pair of isolates which is assigned to the same type by one typing method is also typed as identical by the other. In this context, the congruence between techniques was recorded as follow: MLVA-type to predict MLST-type (93%), PFGE to MLST (92%), DL to MLST (64.2%), PFGE to MLVA (63.5%) and PFGE to DL (61.7%). Conversely, for all above combinations, prediction was very poor. The congruence was increased at the clonal complex (CC) level. MLST is regarded the gold standard for phylogenetic classification of bacteria, but is rather laborious to carry out in many settings. Our data suggest that MLVA can predict the MLST-type with high accuracy, and even higher when studying the clonal complex level. Of the studied three techniques MLVA was therefore the best surrogate method to predict MLST.