Genome sequence of the deep-rooted Yersinia pestis strain Angola reveals new insights into the evolution and pangenome of the plague bacterium

Fibrinolytic and coagulative activities of Yersinia pestis

The outer membrane protease Pla belongs to the omptin protease family spread by horizontal gene transfer into Gram-negative bacteria that infect animals or plants. Pla has adapted to support the life style of the plague bacterium Yersinia pestis. Pla has a β-barrel fold with 10 membrane-spanning β strands and five surface loops, and the barrel surface contains bound lipopolysaccharide (LPS) that is critical for the conformation and the activity of Pla. The biological activity of Pla is influenced by the structure of the surface loops around the active site groove and by temperature-induced LPS modifications. Several of the putative virulence-related functions documented for Pla in vitro address control of the human hemostatic system, i.e., coagulation and fibrinolysis. Pla activates human plasminogen to the serine protease plasmin and activates the physiological plasminogen activator urokinase. Pla also inactivates the protease inhibitors alpha-2-antiplasmin and plasminogen activator inhibitor 1 (PAI-1) and prevents the activation of thrombin-activatable fibrinolysis inhibitor (TAFI). These functions enhance uncontrolled fibrinolysis which is thought to improve Y. pestis dissemination and survival in the mammalian host, and lowered fibrin(ogen) deposition has indeed been observed in mice infected with Pla-positive Y. pestis. However, Pla also inactivates an anticoagulant, the tissue factor (TF) pathway inhibitor, which should increase fibrin formation and clotting. Thus, Pla and Y. pestis have complex interactions with the hemostatic system. Y. pestis modifies its LPS upon transfer to the mammalian host and we hypothesize that the contrasting biological activities of Pla in coagulation and fibrinolysis are influenced by LPS changes during infection.

Short and long-term genome stability analysis of prokaryotic genomes

Background

Gene organization dynamics is actively studied because it provides useful evolutionary information, makes functional annotation easier and often enables to characterize pathogens. There is therefore a strong interest in understanding the variability of this trait and the possible correlations with life-style. Two kinds of events affect genome organization: on one hand translocations and recombinations change the relative position of genes shared by two genomes (i.e. the backbone gene order); on the other, insertions and deletions leave the backbone gene order unchanged but they alter the gene neighborhoods by breaking the syntenic regions. A complete picture about genome organization evolution therefore requires to account for both kinds of events.

Results

We developed an approach where we model chromosomes as graphs on which we compute different stability estimators; we consider genome rearrangements as well as the effect of gene insertions and deletions. In a first part of the paper, we fit a measure of backbone gene order conservation (hereinafter called backbone stability) against phylogenetic distance for over 3000 genome comparisons, improving existing models for the divergence in time of backbone stability. Intra- and inter-specific comparisons were treated separately to focus on different time-scales. The use of multiple genomes of a same species allowed to identify genomes with diverging gene order with respect to their conspecific. The inter-species analysis indicates that pathogens are more often unstable with respect to non-pathogens. In a second part of the text, we show that in pathogens, gene content dynamics (insertions and deletions) have a much more dramatic effect on genome organization stability than backbone rearrangements.

Conclusion

In this work, we studied genome organization divergence taking into account the contribution of both genome order rearrangements and genome content dynamics. By studying species with multiple sequenced genomes available, we were able to explore genome organization stability at different time-scales and to find significant differences for pathogen and non-pathogen species. The output of our framework also allows to identify the conserved gene clusters and/or partial occurrences thereof, making possible to explore how gene clusters assembled during evolution.

Whole-Genome Draft Sequences of 26 Enterohemorrhagic Escherichia coli O157:H7 Strains

First identified in 1982, Escherichia coli O157:H7 is the dominant enterohemorrhagic serotype underlying food-borne human infections in North America. Here, we report the genomes of twenty-six strains derived from patients and the bovine reservoir. These resources enable detailed whole-genome comparisons and permit investigations of genotypic and phenotypic plasticity.

Probabilistic models for CRISPR spacer content evolution

BACKGROUND

The CRISPR/Cas system is known to act as an adaptive and heritable immune system in Eubacteria and Archaea. Immunity is encoded in an array of spacer sequences. Each spacer can provide specific immunity to invasive elements that carry the same or a similar sequence. Even in closely related strains, spacer content is very dynamic and evolves quickly. Standard models of nucleotide evolution cannot be applied to quantify its rate of change since processes other than single nucleotide changes determine its evolution.

METHODS

We present probabilistic models that are specific for spacer content evolution. They account for the different processes of insertion and deletion. Insertions can be constrained to occur on one end only or are allowed to occur throughout the array. One deletion event can affect one spacer or a whole fragment of adjacent spacers. Parameters of the underlying models are estimated for a pair of arrays by maximum likelihood using explicit ancestor enumeration.

RESULTS

Simulations show that parameters are well estimated on average under the models presented here. There is a bias in the rate estimation when including fragment deletions. The models also estimate times between pairs of strains. But with increasing time, spacer overlap goes to zero, and thus there is an upper bound on the distance that can be estimated. Spacer content similarities are displayed in a distance based phylogeny using the estimated times.We use the presented models to analyze different Yersinia pestis data sets and find that the results among them are largely congruent. The models also capture the variation in diversity of spacers among the data sets. A comparison of spacer-based phylogenies and Cas gene phylogenies shows that they resolve very different time scales for this data set.

CONCLUSIONS

The simulations and data analyses show that the presented models are useful for quantifying spacer content evolution and for displaying spacer content similarities of closely related strains in a phylogeny. This allows for comparisons of different CRISPR arrays or for comparisons between CRISPR arrays and nucleotide substitution rates.

Historical variations in mutation rate in an epidemic pathogen, Yersinia pestis

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

Evolution and virulence contributions of the autotransporter proteins YapJ and YapK of Yersinia pestis CO92 and their homologs in Y. pseudotuberculosis IP32953

Yersinia pestis, the causative agent of plague, evolved from the gastrointestinal pathogen Yersinia pseudotuberculosis. Both species have numerous type Va autotransporters, most of which appear to be highly conserved. In Y. pestis CO92, the autotransporter genes yapK and yapJ share a high level of sequence identity. By comparing yapK and yapJ to three homologous genes in Y. pseudotuberculosis IP32953 (YPTB0365, YPTB3285, and YPTB3286), we show that yapK is conserved in Y. pseudotuberculosis, while yapJ is unique to Y. pestis. All of these autotransporters exhibit >96% identity in the C terminus of the protein and identities ranging from 58 to 72% in their N termini. By extending this analysis to include homologous sequences from numerous Y. pestis and Y. pseudotuberculosis strains, we determined that these autotransporters cluster into a YapK (YPTB3285) class and a YapJ (YPTB3286) class. The YPTB3286-like gene of most Y. pestis strains appears to be inactivated, perhaps in favor of maintaining yapJ. Since autotransporters are important for virulence in many bacterial pathogens, including Y. pestis, any change in autotransporter content should be considered for its impact on virulence. Using established mouse models of Y. pestis infection, we demonstrated that despite the high level of sequence identity, yapK is distinct from yapJ in its contribution to disseminated Y. pestis infection. In addition, a mutant lacking both of these genes exhibits an additive attenuation, suggesting nonredundant roles for yapJ and yapK in systemic Y. pestis infection. However, the deletion of the homologous genes in Y. pseudotuberculosis does not seem to impact the virulence of this organism in orogastric or systemic infection models.

Whole-genome sequences of Bacillus subtilis and close relatives

We sequenced four strains of Bacillus subtilis and the type strains for two closely related species, Bacillus vallismortis and Bacillus mojavensis. We report the high-quality Sanger genome sequences of B. subtilis subspecies subtilis RO-NN-1 and AUSI98, B. subtilis subspecies spizizenii TU-B-10(T) and DV1-B-1, Bacillus mojavensis RO-H-1(T), and Bacillus vallismortis DV1-F-3(T).

Plague in the genomic area

With plague being not only a subject of interest for historians, but still a disease of public health concern in several countries, mainly in Africa, there were hopes that analyses of the Yersinia pestis genomes would put an end to this deadly epidemic pathogen. Genomics revealed that Y. pestis isolates evolved from Yersinia pseudotuberculosis in Central Asia some millennia ago, after the acquisition of two Y. pestis-specific plasmids balanced genomic reduction parallel with the expansion of insertion sequences, illustrating the modern concept that, except for the acquisition of plasmid-borne toxin-encoding genes, the increased virulence of Y. pestis resulted from gene loss rather than gene acquisition. The telluric persistence of Y. pestis reminds us of this close relationship, and matters in terms of plague epidemiology. Whereas biotype Orientalis isolates spread worldwide, the Antiqua and Medievalis isolates showed more limited expansion. In addition to animal ectoparasites, human ectoparasites such as the body louse may have participated in this expansion and in devastating historical epidemics. The recent analysis of a Black Death genome indicated that it was more closely related to the Orientalis branch than to the Medievalis branch. Modern Y. pestis isolates grossly exhibit the same gene content, but still undergo micro-evolution in geographically limited areas by differing in the genome architecture, owing to inversions near insertion sequences and the stabilization of the YpfPhi prophage in Orientalis biotype isolates. Genomics have provided several new molecular tools for the genotyping and phylogeographical tracing of isolates and description of plague foci. However, genomics and post-genomics approaches have not yet provided new tools for the prevention, diagnosis and management of plague patients and the plague epidemics still raging in some sub-Saharan countries.

Averting Behavior Framework for Perceived Risk of Yersinia enterocolitica Infections

The focus of this research is to present a theoretical model of averting actions that households take to avoid exposure to Yersinia enterocolitica in contaminated food. The cost of illness approach only takes into account the value of a cure, while the averting behavior approach can estimate the value of preventing the illness. The household, rather than the individual, is the unit of analysis in this model, where one household member is primarily responsible for procuring uncontaminated food for their family. Since children are particularly susceptible and live with parents who are primary decision makers for sustenance, the designated household head makes the choices that are investigated in this paper. This model uses constrained optimization to characterize activities that may offer protection from exposure to Yersinia enterocolitica contaminated food. A representative household decision maker is assumed to allocate family resources to maximize utility of an altruistic parent, an assumption used in most research involving economics of the family.

Novel plasmids and resistance phenotypes in Yersinia pestis: unique plasmid inventory of strain Java 9 mediates high levels of arsenic resistance

Growing evidence suggests that the plasmid repertoire of Yersinia pestis is not restricted to the three classical virulence plasmids. The Java 9 strain of Y. pestis is a biovar Orientalis isolate obtained from a rat in Indonesia. Although it lacks the Y. pestis-specific plasmid pMT, which encodes the F1 capsule, it retains virulence in mouse and non-human primate animal models. While comparing diverse Y. pestis strains using subtractive hybridization, we identified sequences in Java 9 that were homologous to a Y. enterocolitica strain carrying the transposon Tn2502, which is known to encode arsenic resistance. Here we demonstrate that Java 9 exhibits high levels of arsenic and arsenite resistance mediated by a novel promiscuous class II transposon, named Tn2503. Arsenic resistance was self-transmissible from Java 9 to other Y. pestis strains via conjugation. Genomic analysis of the atypical plasmid inventory of Java 9 identified pCD and pPCP plasmids of atypical size and two previously uncharacterized cryptic plasmids. Unlike the Tn2502-mediated arsenic resistance encoded on the Y. enterocolitica virulence plasmid; the resistance loci in Java 9 are found on all four indigenous plasmids, including the two novel cryptic plasmids. This unique mobilome introduces more than 105 genes into the species gene pool. The majority of these are encoded by the two entirely novel self-transmissible plasmids, which show partial homology and synteny to other enterics. In contrast to the reductive evolution in Y. pestis, this study underlines the major impact of a dynamic mobilome and lateral acquisition in the genome evolution of the plague bacterium.

Genomic anatomy of Escherichia coli O157:H7 outbreaks

The rapid emergence of Escherichia coli O157:H7 from an unknown strain in 1982 to the dominant hemorrhagic E. coli serotype in the United States and the cause of widespread outbreaks of human food-borne illness highlights a need to evaluate critically the extent to which genomic plasticity of this important enteric pathogen contributes to its pathogenic potential and its evolution as well as its adaptation in different ecological niches. Aimed at a better understanding of the evolution of the E. coli O157:H7 pathogenome, the present study presents the high-quality sequencing and comparative phylogenomic analysis of a comprehensive panel of 25 E. coli O157:H7 strains associated with three nearly simultaneous food-borne outbreaks of human disease in the United States. Here we present a population genetic analysis of more than 200 related strains recovered from patients, contaminated produce, and zoonotic sources. High-resolution phylogenomic approaches allow the dynamics of pathogenome evolution to be followed at a high level of phylogenetic accuracy and resolution. SNP discovery and study of genome architecture and prophage content identified numerous biomarkers to assess the extent of genetic diversity within a set of clinical and environmental strains. A total of 1,225 SNPs were identified in the present study and are now available for typing of the E. coli O157:H7 lineage. These data should prove useful for the development of a refined phylogenomic framework for forensic, diagnostic, and epidemiological studies to define better risk in response to novel and emerging E. coli O157:H7 resistance and virulence phenotypes.

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.

Everything at once: comparative analysis of the genomes of bacterial pathogens

The sum of unique genes in all genomes of a bacterial species is referred to as the pan-genome and is comprised of variably absent or present accessory genes and universally present core genes. The accessory genome is an important source of genetic variability in bacterial populations, allowing sub-populations of bacteria to better adapt to specific niches. Such subgroups may themselves have a relatively stable core genome that may influence host preference, virulence, or an association with specific disease syndromes. The core genome provides a useful means of phylogenetic reconstruction as well as contributing to phenotypic heterogeneity. Variation within the pan-genome forms the basis of comparative genotyping techniques, which have evolved alongside technology. Current high-throughput sequencing platforms have created an unprecedented opportunity for comparisons among multiple, closely related genomes. The computer algorithms and software for such comparisons continue to evolve and promise exciting advances in the world of bacterial comparative genomics. We review genotyping techniques based upon phenotypic traits, both core and accessory genomes, and look at some of the software programs currently available to perform whole-genome comparative analyses.

Genome signatures of Escherichia coli O157:H7 isolates from the bovine host reservoir

Cattle comprise a main reservoir of Shiga toxin-producing Escherichia coli O157:H7 (STEC). The significant differences in host prevalence, transmissibility, and virulence phenotypes among strains from bovine and human sources are of major interest to the public health community and livestock industry. Genomic analysis revealed divergence into three lineages: lineage I and lineage I/II strains are commonly associated with human disease, while lineage II strains are overrepresented in the asymptomatic bovine host reservoir. Growing evidence suggests that genotypic differences between these lineages, such as polymorphisms in Shiga toxin subtypes and synergistically acting virulence factors, are correlated with phenotypic differences in virulence, host ecology, and epidemiology. To assess the genomic plasticity on a genome-wide scale, we have sequenced the whole genome of strain EC869, a bovine-associated E. coli O157:H7 isolate. Comparative phylogenomic analysis of this key isolate enabled us to place accurately bovine lineage II strains within the genetically homogenous E. coli O157:H7 clade. Identification of polymorphic loci that are anchored both in the chromosomal backbone and horizontally acquired regions allowed us to associate bovine genotypes with altered virulence phenotypes and host prevalence. This study catalogued numerous novel lineage II-specific genome signatures, some of which appear to be associated intimately with the altered pathogenic potential and niche adaptation within the bovine rumen. The presented extended list of polymorphic markers is valuable in the development of a robust typing system critical for forensic, diagnostic, and epidemiological studies of this emerging human pathogen.

Population structure and evolution of pathogenicity of Yersinia pseudotuberculosis

Yersinia pseudotuberculosis is an enteric human pathogen but is widespread in the environment. Pathogenicity is determined by a number of virulence factors, including the virulence plasmid pYV, the high-pathogenicity island (HPI), and the Y. pseudotuberculosis-derived mitogen (YPM), a superantigen. The presence of the 3 virulence factors varies among Y. pseudotuberculosis isolates. We developed a multilocus sequence typing (MLST) scheme to address the population structure of Y. pseudotuberculosis and the evolution of its pathogenicity. The seven housekeeping genes selected for MLST were mdh, recA, sucA, fumC, aroC, pgi, and gyrB. An MLST analysis of 83 isolates of Y. pseudotuberculosis, representing 19 different serotypes and six different genetic groups, identified 61 sequence types (STs) and 12 clonal complexes. Out of 26 allelic changes that occurred in the 12 clonal complexes, 13 were mutational events while 13 were recombinational events, indicating that recombination and mutation contributed equally to the diversification of the clonal complexes. The isolates were separated into 2 distinctive clusters, A and B. Cluster A is the major cluster, with 53 STs (including Y. pestis strains), and is distributed worldwide, while cluster B is restricted to the Far East. The YPM gene is widely distributed on the phylogenetic tree, with ypmA in cluster A and ypmB in cluster B. pYV is present in cluster A only but is sporadically absent in some cluster A isolates. In contrast, an HPI is present only in a limited number of lineages and must be gained by lateral transfer. Three STs carry all 3 virulence factors and can be regarded as high-pathogenicity clones. Isolates from the same ST may not carry all 3 virulence factors, indicating frequent gain or loss of these factors. The differences in pathogenicity among Y. pseudotuberculosis strains are likely due to the variable presence and instability of the virulence factors.

Characterization of pPCP1 Plasmids in Yersinia pestis Strains Isolated from the Former Soviet Union

Complete sequences of 9.5-kb pPCP1 plasmids in three Yersinia pestis strains from the former Soviet Union (FSU) were determined and compared with those of pPCP1 plasmids in three well-characterized, non-FSU Y. pestis strains (KIM, CO92, and 91001). Two of the FSU plasmids were from strains C2614 and C2944, isolated from plague foci in Russia, and one plasmid was from strain C790 from Kyrgyzstan. Sequence analyses identified four sequence types among the six plasmids. The pPCP1 plasmids in the FSU strains were most genetically related to the pPCP1 plasmid in the KIM strain and least related to the pPCP1 plasmid in Y. pestis 91001. The FSU strains generally had larger pPCP1 plasmid copy numbers compared to strain CO92. Expression of the plasmid's pla gene was significantly (P ≤ .05) higher in strain C2944 than in strain CO92. Given pla's role in Y. pestis virulence, this difference may have important implications for the strain's virulence.

Yersinia pestis genome sequencing identifies patterns of global phylogenetic diversity

Pandemic infectious diseases have accompanied humans since their origins1, and have shaped the form of civilizations2. Of these, plague is possibly historically the most dramatic. We reconstructed historical patterns of plague transmission through sequence variation in 17 complete genome sequences and 933 single nucleotide polymorphisms (SNPs) within a global collection of 286 Yersinia pestis isolates. Y. pestis evolved in or near China, and has been transmitted via multiple epidemics that followed various routes, probably including transmissions to West Asia via the Silk Road and to Africa by Chinese marine voyages. In 1894, Y. pestis spread to India and radiated to diverse parts of the globe, leading to country-specific lineages that can be traced by lineage-specific SNPs. All 626 current isolates from the U.S.A. reflect one radiation and 82 isolates from Madagascar represent a second. Subsequent local microevolution of Y. pestis is marked by sequential, geographically-specific SNPs.

Delineation and analysis of chromosomal regions specifying Yersinia pestis

Yersinia pestis, the causative agent of plague, has recently diverged from the less virulent enteropathogen Yersinia pseudotuberculosis. Its emergence has been characterized by massive genetic loss and inactivation and limited gene acquisition. The acquired genes include two plasmids, a filamentous phage, and a few chromosomal loci. The aim of this study was to characterize the chromosomal regions acquired by Y. pestis. Following in silico comparative analysis and PCR screening of 98 strains of Y. pseudotuberculosis and Y. pestis, we found that eight chromosomal loci (six regions [R1pe to R6pe] and two coding sequences [CDS1pe and CDS2pe]) specified Y. pestis. Signatures of integration by site specific or homologous recombination were identified for most of them. These acquisitions and the loss of ancestral DNA sequences were concentrated in a chromosomal region opposite to the origin of replication. The specific regions were acquired very early during Y. pestis evolution and were retained during its microevolution, suggesting that they might bring some selective advantages. Only one region (R3pe), predicted to carry a lambdoid prophage, is most likely no longer functional because of mutations. With the exception of R1pe and R2pe, which have the potential to encode a restriction/modification and a sugar transport system, respectively, no functions could be predicted for the other Y. pestis-specific loci. To determine the role of the eight chromosomal loci in the physiology and pathogenicity of the plague bacillus, each of them was individually deleted from the bacterial chromosome. None of the deletants exhibited defects during growth in vitro. Using the Xenopsylla cheopis flea model, all deletants retained the capacity to produce a stable and persistent infection and to block fleas. Similarly, none of the deletants caused any acute flea toxicity. In the mouse model of infection, all deletants were fully virulent upon subcutaneous or aerosol infections. Therefore, our results suggest that acquisition of new chromosomal materials has not been of major importance in the dramatic change of life cycle that has accompanied the emergence of Y. pestis.

The importance of the small RNA chaperone Hfq for growth of epidemic Yersinia pestis, but not Yersinia pseudotuberculosis, with implications for plague biology

Yersinia pestis, the etiologic agent of plague, has only recently evolved from Yersinia pseudotuberculosis. hfq deletion caused severe growth restriction at 37 degrees C in Y. pestis but not in Y. pseudotuberculosis. Strains from all epidemic plague biovars were similarly affected, implicating Hfq, and likely small RNAs (sRNAs), in the unique biology of the plague bacillus.

Ergatis: a web interface and scalable software system for bioinformatics workflows

MOTIVATION

The growth of sequence data has been accompanied by an increasing need to analyze data on distributed computer clusters. The use of these systems for routine analysis requires scalable and robust software for data management of large datasets. Software is also needed to simplify data management and make large-scale bioinformatics analysis accessible and reproducible to a wide class of target users.

RESULTS

We have developed a workflow management system named Ergatis that enables users to build, execute and monitor pipelines for computational analysis of genomics data. Ergatis contains preconfigured components and template pipelines for a number of common bioinformatics tasks such as prokaryotic genome annotation and genome comparisons. Outputs from many of these components can be loaded into a Chado relational database. Ergatis was designed to be accessible to a broad class of users and provides a user friendly, web-based interface. Ergatis supports high-throughput batch processing on distributed compute clusters and has been used for data management in a number of genome annotation and comparative genomics projects.

AVAILABILITY

Ergatis is an open-source project and is freely available at http://ergatis.sourceforge.net.