Role of genomic typing in taxonomy, evolutionary genetics, and microbial epidemiology

Prediction of viral families and hosts of single-stranded RNA viruses based on K-Mer coding from phylogenetic gene sequences

There are billions of virus species worldwide, and viruses, the smallest parasitic entities, pose a serious threat. Therefore, fighting associated disorders requires an understanding of the genetic structure of viruses. Considering the wide diversity and rapid evolution of viruses, there is a critical need to quickly and accurately classify viral species and their potential hosts to better understand transmission dynamics, facilitating the development of targeted therapies. Recognizing this, this study has investigated the classes of RNA viruses based on their genomic sequences using Machine Learning (ML) and Deep Learning (DL) models. The PhyVirus dataset, consisting of pathogenic Single-stranded RNA viruses of Baltimore group four (+ssRNA) and five (-ssRNA) with different hosts and species, was analyzed. The dataset containing viral gene sequences was analyzed using the K-Mer coding technique, which is based on base words of various lengths. The study used classical ML algorithms (Random Forest, Gradient Boosting and Extra Trees) and the Fully Connected Deep Neural Network, a Deep Learning algorithm, to predict viral families and hosts. Detailed analyses were performed on the classifier performance in scenarios with different train-test ratios and different word lengths (k-values) for K-Mer. The observed results show that Fully Connected Deep Neural Network has a high success rate of 99.60 % in predicting virus families. In predicting virus hosts, the Extra Trees classifier achieved the highest success rate of 81.53 %. This study is considered to be the first classification study in the literature on this dataset, which has a very large family and host diversity consisting of gene sequences of Single-stranded RNA viruses. Our detailed investigations on how varying word lengths based on K-Mer coding in gene sequences affect the classification into viral families and hosts make this study particularly valuable. This study shows that ML and DL methods have the potential to produce valuable results in phylogenetic studies. In addition, the results and high-performance values show that these methods can be successfully used in regenerative applications of gene sequences or in studies such as the elimination of losses in gene sequences.

Identification and characterisation of emerging fish pathogens Aeromonas veronii and Aeromonas hydrophila isolated from naturally infected Channa punctata

Naturally infected Channa punctata exhibiting bacterial septicemic syndrome including ulcerations along with mortality records were collected from a fish farm in Assam during winter season (early November 2020 to early January 2021). The moribund fishes were subjected for bacterial isolation followed by identification of the bacteria. Two dominant emerging bacterial pathogens were identified as Aeromonas veronii (isolate ZooGURD-01) and Aeromonas hydrophila (isolate ZooGURD-05) by standard biochemical characterization and 16S rRNA and rpo B gene amplification. Re-infection experiments of both the bacterial isolates in healthy disease-free C. punctata showed similar symptoms to that of natural infection thus confirming their virulence. The LD50 calculated during challenge test for both the isolates ZooGURD-01 and ZooGURD-05 found to be pathogenic at 2.6 × 104 and 1.6 × 104 CFU/fish respectively. Further PCR amplification of specific virulent genes (aerolysin, hemolysin and enterotoxin) confirmed pathogenicity for both isolates. Histopathological examinations of liver and kidney in re-infection experiments showed prominent changes supporting bacterial septicaemia. Antibiotic sensitivity pattern showed that the isolates ZooGURD-01 and ZooGURD-05 were sensitive to 22 and 19 out of 25 antimicrobials respectively. The present study was the first report on the mortality of farmed C. punctata associated with natural infection caused by A. veronii and A. hydrophila with no record of pathogenicity of A. veronii in C. punctata.

Phylogenomic analysis of Citrobacter sp. strain AAK_AS5 and its metabolic capabilities to support nitrogen removal behavior

Despite the ubiquity of the genus Citrobacter in clinical, industrial, and environmental scenarios, a large number of Citrobacter strains have not been explored at the genome-scale level. In this study, accurate taxonomic assignment of strain AAK_AS5 isolated from activated sludge was achieved by in-silico genomic comparison using Overall Genome-based Relatedness Indices (ANI(OAT): 97.55%, ANIb:97.28%, and ANIm: 97.83%) that indicated its closest identity to the related strain Citrobacter portucalensis A60^T . Results were consistent with a digital DNA-DNA hybridization value of 80% with C. portucalensis A60^T which was greater than the species boundary value >70% for delineating closely related bacterial species. Gene mining through Kyoto Encyclopedia of Genes and Genomes (KEGG), and annotation using rapid annotation subsystem technology (RAST) revealed the notable gene contents for nitrogen metabolism and other pathways associated with nitrate/nitrite ammonification (28 genes), ammonia assimilation (22 genes), and denitrification pathways (14 genes). Furthermore, the strain AAK_AS5 also exhibited a high soluble chemical oxygen demand (sCOD), NH₄ ⁺ -N, and NO₃ ^- -N removal efficiency of 91.4%, 90%, and 93.6%, respectively thus validating its genetic capability for utilizing both (NH₄ )₂ SO₄ and KNO₃ as the nitrogen source. The study provided deeper insights into the phylogenomics and the genetic potential of Citrobacter, sp. strain AAK AS5 associated with nitrogen metabolism thus signifying the potential application of the isolate for treating nitrogen-rich wastewaters.

A minisatellite-based MLVA for deciphering the global epidemiology of the bacterial cassava pathogen Xanthomonas phaseoli pv. manihotis

Cassava Bacterial Blight (CBB) is a destructive disease widely distributed in the different areas where this crop is grown. Populations studies have been performed at local and national scales revealing a geographical genetic structure with temporal variations. A global epidemiology analysis of its causal agent Xanthomonas phaseoli pv. manihotis (Xpm) is needed to better understand the expansion of the disease for improving the monitoring of CBB. We targeted new tandem repeat (TR) loci with large repeat units, i.e. minisatellites, that we multiplexed in a scheme of Multi-Locus Variable number of TR Analysis (MLVA-8). This genotyping scheme separated 31 multilocus haplotypes in three clusters of single-locus variants and a singleton within a worldwide collection of 93 Xpm strains isolated over a period of fifty years. The major MLVA-8 cluster 1 grouped strains originating from all countries, except the unique Chinese strain. On the contrary, all the Xpm strains genotyped using the previously developed MLVA-14 microsatellite scheme were separated as unique haplotypes. We further propose an MLVA-12 scheme which takes advantage of combining TR loci with different mutation rates: the eight minisatellites and four faster evolving microsatellite markers, for global epidemiological surveillance. This MLVA-12 scheme identified 78 haplotypes and separated most of the strains in groups of double-locus variants (DLV) supporting some phylogenetic relationships. DLV groups were subdivided into closely related clusters of strains most often sharing the same geographical origin and isolated over a short period, supporting epidemiological relationships. The main MLVA-12 DLV group#1 was composed by strains from South America and all the African strains. The MLVA-12 scheme combining both minisatellite and microsatellite loci with different discriminatory power is expected to increase the accuracy of the phylogenetic signal and to minimize the homoplasy effects. Further investigation of the global epidemiology of Xpm will be helpful for a better control of CBB worldwide.

The Dissemination and Molecular Characterization of Clonal Complex 361 (CC361) Methicillin-Resistant Staphylococcus aureus (MRSA) in Kuwait Hospitals

Methicillin-resistant Staphylococcus aureus (MRSA) belonging to clonal complex 361 (CC361-MRSA) is rare among patients’ populations globally. However, CC361-MRSA has been isolated with an increasing trend among patients in Kuwait hospitals since 2010. This study investigated the molecular characteristics of CC361-MRSA isolated from patients in Kuwait hospitals in 2016–2018 to understand their genetic relatedness and virulence determinants. Of 5,223 MRSA isolates investigated by DNA microarray, 182 (3.4%) isolates obtained in 2016 (N = 55), 2017 (N = 56), and 2018 (N = 71) were identified as CC361-MRSA. The CC361-MRSA isolates were analyzed further using antibiogram, spa typing and multi locus sequence typing (MLST). Most of the isolates were resistant to fusidic acid (64.8%), kanamycin (43.4%), erythromycin (36.3%), and clindamycin (14.3%) encoded by fusC, aphA3, and erm(B)/erm(C) respectively. Nine isolates (4.9%) were resistant to linezolid mediated by cfr. The isolates belonged to 22 spa types with t3841 (N = 113), t315 (N = 16), t1309 (N = 14), and t3175 (N = 5) constituting 81.3% of the spa types, four genotypes (strain types), CC361-MRSA-[V/VT + fus] (N = 112), CC361-MRSA-IV, WA MRSA-29 (N = 36), CC361-MRSA-V, WA MRSA-70/110 (N = 33) and CC361-MRSA-[V + fus] variant (N = 1). MLST conducted on 69 representative isolates yielded two sequence types: ST361 (11/69) and ST672 (58/69). All CC361-MRSA isolates were positive for cap8, agr1, and the enterotoxin egc gene cluster (seg, sei, selm, seln, selo, and selu). The tst1 was detected in 19 isolates. The immune evasion cluster (IEC) genes type B (scn, chp, and sak) and type E (scn and sak) were detected in 20 and 152 isolates, respectively. The CC361-MRSA circulating in Kuwait hospitals consisted of two closely related sequence types, ST361 and ST672 with ST672-MRSA [V/VT + fus] as the dominant genotype. The dissemination of these newly emerged clones and the emergence of linezolid resistance limits therapeutic options, as well as present significant challenges for the control of MRSA infections in Kuwait hospitals.

Characterisation, pathogenicity and hydrolytic enzyme profiling of selected Fusarium species and their inhibition by novel coumarins

Three Fusarium species isolated locally were characterised by the amplification of their rDNA ITS region, host specificity, and hydrolytic enzyme production. The strains were identified as Fusarium pseudoanthophilum, which is being reported for the first time in South Africa, as well as F. foetens and F. fujikuroi. All the three strains were capable of infecting vegetables such as tomatoes, bell and cayenne peppers, belonging to the Solanaceae family. The Fusarium strains also showed significant production of cell wall degrading enzymes in vitro, such as amylase, cellulase, xylanase, and polygalacturonase, thus highlighting the possibilities of these enzymes as pathogenic factors. Subsequently, the strains were discovered to be susceptible to three halogenated coumarins. The most effective of the tested coumarins, 6-bromo3-2,2-dibromoacetyl-2H-chromen-2-one, showed MIC values of 0.125, 0.0625 and 0.125 mg/ml against F. foetens, F. pseudoanthophilum and F. fujikuroi, respectively. The antifungal potentials of the halogenated coumarins were confirmed in silico through PASS analysis, toxicity prediction and docking studies. Findings from this study demonstrate the use of these coumarins as potential control agents against the Fusarium species and other pathogenic fungi in general.

Genomic investigation of clinically significant coagulase-negative staphylococci

Introduction. Coagulase-negative staphylococci have been recognized both as emerging pathogens and contaminants of clinical samples. High-resolution genomic investigation may provide insights into their clinical significance.Aims. To review the literature regarding coagulase-negative staphylococcal infection and the utility of genomic methods to aid diagnosis and management, and to identify promising areas for future research.Methodology. We searched Google Scholar with the terms (Staphylococcus) AND (sequencing OR (infection)). We prioritized papers that addressed coagulase-negative staphylococci, genomic analysis, or infection.Results. A number of studies have investigated specimen-related, phenotypic and genetic factors associated with colonization, infection and virulence, but diagnosis remains problematic.Conclusion. Genomic investigation provides insights into the genetic diversity and natural history of colonization and infection. Such information allows the development of new methodologies to identify and compare relatedness and predict antimicrobial resistance. Future clinical studies that employ suitable sampling frames coupled with the application of high-resolution whole-genome sequencing may aid the development of more discriminatory diagnostic approaches to coagulase-staphylococcal infection.

The bacterial communities of surface soils from desert sites in the eastern Utah (USA) portion of the Colorado Plateau

Desert-like areas located in the eastern portion of the state of Utah (USA) have geographic features that can resemble the surface of the planet Mars, characterized by red-colored hills, soils and sandstones. We examined the bacterial biodiversity of surface soil samples from several sites from the Colorado Plateau Desert in eastern Utah using pyrosequencing of PCR amplified bacterial 16S rRNA genes from total extracted soil DNA. The sample sites cover the Great Basin, Goblin Valley State Park and nearby regions on the Colorado Plateau. We also examined several physicochemical parameters of the soil samples to investigate any possible correlations between bacterial community structure and environmental drivers. The predominant bacterial phyla present in the samples were found to belong to members of the Proteobacteria, Actinobacteria, Bacteroidetes, and Gemmatimonadetes. The most abundant genera in our samples were found to belong to the Cesiribacter, Lysobacter, Adhaeribacter, Microvirga and Pontibacter genera. We found that the relative abundance of Proteobacteria and Gemmatimonadetes were significantly correlated with soil pH and a low concentration of organic matter, suggesting that, in these relatively high-altitude desert soils, these two parameters may be of primary importance to influence bacterial community composition.

Acute post-streptococcal glomerulonephritis: analysis of the pathogenesis

Increasing evidence supports a central role of the immune system in acute post streptococcal glomerulonephritis (APSGN), but the current view of how streptococcal biology affects immunity, and vice versa, remains to be clarified. Renal glomerular immune complex deposition is critical in the initiation of APSGN; however, mechanisms previous to immune complex formation could modulate the initiation and the progression of the disease. Initial and late renal events involved in the nephritis can also be related to host factors and streptococcal factors. In this review we describe the mechanisms reported for the APSGN pathogenesis, the interactions of streptococcal products with renal cells and leukocytes, the possible effects of different nephritogenic antigens in the renal environment and the possibility that APSGN is not just due to a single streptococcal antigen and its antibody; instead, kidney damage may be the result of different factors acting at the same time related to both streptococcus and host factors. Addressing these points should help us to better understand APSGN physiopathology.

Identification and characterization of virulent Aeromonas hydrophila Ah17 from infected Channa striata in river Cauvery and in vitro evaluation of shrimp chitosan

Aeromonas hydrophila, an inhabitant in the aquatic ecosystem is considered as an important foodborne bacterial zoonotic pathogen in aquaculture. The present study aimed to identify virulent A. hydrophila from naturally infected Channa striata in river Cauvery and in vitro evaluation of shrimp chitosan. Rimler Shotts (RS) and blood agar medium identified the presence of pathogenic Aeromonas sp. from the infected C. striata. A. hydrophila Ah17 was identified using 16S rRNA nucleotide sequence. Extracellular enzymes such as amylase, lipase, and protease were screened in A. hydrophila Ah17. Antibiotic susceptibility tests showed A. hydrophila Ah17 was highly resistant against β-lactam, glycopeptide, macrolides, phosphonic, fucidin, and oxazolidinone classes of antibiotics. Virulent genes such as hemolysin (aer and hly), heat-labile enterotoxin (act), cytotonic heat-stable enterotoxin (ast), elastase (ahyB), and lipase (lip) were identified. Growth and the viable cell population of virulent A. hydrophila Ah17 were significantly reduced in a dose-dependent manner against shrimp chitosan (CHS) from Penaeus indicus (P. indicus). Thus, the present study isolated virulent A. hydrophila Ah17 from the environmental source and characterized in vitro with shrimp chitosan.

Within-species contamination of bacterial whole-genome sequence data has a greater influence on clustering analyses than between-species contamination

Although it is assumed that contamination in bacterial whole-genome sequencing causes errors, the influences of contamination on clustering analyses, such as single-nucleotide polymorphism discovery, phylogenetics, and multi-locus sequencing typing, have not been quantified. By developing and analyzing 720 Listeria monocytogenes, Salmonella enterica, and Escherichia coli short-read datasets, we demonstrate that within-species contamination causes errors that confound clustering analyses, while between-species contamination generally does not. Contaminant reads mapping to references or becoming incorporated into chimeric sequences during assembly are the sources of those errors. Contamination sufficient to influence clustering analyses is present in public sequence databases.

Variation in Mycobacterium bovis genetic richness suggests that inwards cattle movements are a more important source of infection in beef herds than in dairy herds

Background

We used genetic Multi-Locus VNTR Analysis (MLVA) data gathered from surveillance efforts to better understand the ongoing bovine tuberculosis (bTB) epidemic in Northern Irish cattle herds. We modelled the factors associated with Mycobacterium bovis MLVA genotype richness at three analytical scales; breakdown level, herd level, and patch level, and compared the results between dairy and non-dairy production types.

Results

In 83% of breakdowns and in 63% of herds, a single MLVA genotype was isolated. Five or more MLVA genotypes were found in less than 3 % of herds. Herd size and the total number of reactors were important explanatory variables, suggesting that increasing MLVA genotype richness was positively related to increases in the number of host animals. Despite their smaller relative size, however, the highest MLVA genotype richness values were observed in non-dairy herds. Increasing inwards cattle movements were important positive predictors of MLVA genotype richness, but mainly in non-dairy settings.

Conclusions

The principal finding is that low MLVA genotype richness indicates that small-scale epidemics, e.g. wildlife, contiguous farms, and within-herd recrudescence, are important routes of M. bovis infection in cattle herds. We hypothesise that these mechanisms will maintain, but may not explicitly increase, MLVA genotype richness. The presence of elevated MLVA richness is relatively rare and likely indicates beef fattening enterprises, which purchase cattle from over long distances. Cattle movements were furthermore an important predictor of MLVA genotype richness in non-dairy herds, but not in dairy herds; this may represent reduced cattle purchasing levels in dairy enterprises, compared to beef. These observations allude to the relative contribution of different routes of bTB infection between production types; we posit that infection associated with local factors may be more evident in dairy herds than beef herds, however in beef herds, inwards movements offer additional opportunities for introducing M. bovis into the herd.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1530-7) contains supplementary material, which is available to authorized users.

Coagulase gene polymorphisms of Staphylococcus aureus isolates from patients at Kosti Teaching Hospital, Sudan

Background

Staphylococcus aureus is a common cause of nosocomial infections. Epidemiological typing of S. aureus enables control of its spread. The objective of this study was to investigate coagulase gene polymorphisms of S. aureus isolated from patients at Kosti Hospital in Sudan.

Methods

In total, 110 S. aureus isolates were recovered from 110 patients who were enrolled in the study. S. aureus strains were isolated on blood agar and MacConkey agar and then identified by conventional tests. Resistance to methicillin was determined by detection of the mecA gene. Polymorphism in the coagulase gene (coa) was investigated using PCR followed by AluI RFLP analysis.

Results

Methicillin-resistant S. aureus accounted for 62/110 (56 %) of the isolates. PCR of the coa gene showed two different amplicons, one of 680 bp detected in 83/110 (75.5 %) of the isolates and one of 790 bp detected in 27/110 (24.5 %). When digested with the AluI enzyme, the 790 bp amplicon resulted in three restriction fragments of 500, 210 and 80 bp (coa1). Restriction of the 680 bp amplicon gave two patterns; the first (coa2) was found in 22/110 (20 %) of the isolates and showed four fragments of 210, 210, 180 and 80 bp, and the second (coa3) was found in 61/110 (55.5 %) and revealed three fragments of 390, 210 and 80 bp. Most of the coa3 isolates (75.4%) were methicillin-resistant.

Conclusion

Three polymorphic genotypes of S. aureus were identified in patients at Kosti Hospital. The coa3 genotype was the predominant one and was mostly detected in methicillin-resistant isolates.

Molecular epidemiology of methicillin-resistant Staphylococcus aureus isolates in New South Wales, Australia, 2012-2017

BACKGROUND

To better understand the molecular epidemiology of MRSA and to assess the utility of 19-target binary typing we undertook large-scale epidemiological surveillance of MRSA from invasive and non-invasive clinical specimens, and screening swabs.

METHODS

Binary typing was performed on clinical MRSA isolates collected in New South Wales (NSW), Australia between 01/01/2012 - 31/12/2017. Binary type (BT) predicted multilocus sequence type (ST) and spa types based on results from isolates which had been characterised by both methods.

RESULTS

7624 MRSA isolates were analysed of which 3581 (47%) were wounds or skin & softtissue isolates (W/SSTI), 2436 (32%) screening swabs, 469 (6%) blood cultures (BC), 780 (10%) others, and 358 (5%) unknown. We identified 731 BTs, 54 spa types, and 31 STs. ST239 was the commonest MRSA clone in 2012 (30%), but it decreased to 7% in 2017 (p <0.001). In contrast, <0.5% of MRSA were ST45 in 2012 compared to 14% in 2017 (p<0.001). An emergence of PVL-positive ST22 was also noted. Of all isolates, 28% (2122/7624) were lukS/PVL positive; the proportion, among prospectively collected isolates increased from 24% (1406/5858) to 33% (1933/5858) between 2012 and 2017 (p <0.0001). 43% (1534/3581) W/SSTI, 20% (95/469) BC and 10% (239/2436) screening swabs were PVL-positive.

CONCLUSIONS

A major change in the epidemiology of MRSA was noted with a decline of ST239, an emergence of ST45 and PVL-positive ST22, and a significant increase in PVL-positive isolates. Binary typing can be a useful routine laboratory test for prospective molecular surveillance of MRSA colonisation and infection.

Toward an Alignment-Free Method for Feature Extraction and Accurate Classification of Viral Sequences

The classification of pathogens in emerging and re-emerging viruses represents major interests in taxonomic studies, functional genomics, host-pathogen interplay, prevention, and disease treatments. It consists of assigning a given sequence to its related group of known sequences sharing similar characteristics and traits. The challenges to such classification could be associated with several virus properties including recombination, mutation rate, multiplicity of motifs, and diversity. In domains such as pathogen monitoring and surveillance, it is important to detect and quantify known and novel taxa without exploiting the full and accurate alignments or virus family profiles. In this study, we propose an alignment-free method, CASTOR-KRFE, to detect discriminating subsequences within known pathogen sequences to classify accurately unknown pathogen sequences. This method includes three major steps: (1) vectorization of known viral genomic sequences based on k-mers to constitute the potential features, (2) efficient way of pattern extraction and evaluation maximizing classification performance, and (3) prediction of the minimal set of features fitting a given criterion (threshold of performance metric and maximum number of features). We assessed this method through a jackknife data partitioning on a dozen of various virus data sets, covering the seven major virus groups and including influenza virus, Ebola virus, human immunodeficiency virus 1, hepatitis C virus, hepatitis B virus, and human papillomavirus. CASTOR-KRFE provides a weighted average F-measure >0.96 over a wide range of viruses. Our method also shows better performance on complex virus data sets than multiple subsequences extractor for classification (MISSEL), a subsequence extraction method, and the Discriminative mode of MEME patterns extraction tool.

An Optimized Microsatellite Scheme for Assessing Populations of Xanthomonas phaseoli pv. manihotis

Diverse molecular markers have been used to analyze the genetic diversity of plant pathogens. Compared with traditional fingerprinting methods, multiple loci variable number of tandem repeat analyses (MLVAs) have gained importance recently due to their reproducibility, high discriminatory power, ease of performance, low cost, and throughput potential. These characteristics are desirable for continuous pathogen monitoring, especially for pathogens with relatively low genetic diversity, and for disease epidemiology studies. Genetic diversity studies of Xanthomonas phaseoli pv. manihotis, which is the causal agent of cassava bacterial blight, have shown variability and changes in the bacterial population over time. Thus, an easy and fast method needs to be developed to type populations of this pathogen in different countries of the world, especially on small scales. In this study, we developed an MLVA scheme to analyze X. phaseoli pv. manihotis variability on a local scale. The MLVA-15 scheme comprises 15 variable number of tandem repeat loci grouped into four multiplex polymerase chain reaction pools. We showed that the MLVA-15 scheme had slightly higher discriminatory ability at the locality level when compared with amplified fragment length polymorphisms. The MLVA-15 scheme allowed for an accurate determination of the number of genotypes in the sample and showed reproducibility and portability. Additionally, this scheme could be used to analyze numerous strains in a reasonable timeframe. The MLVA-15 scheme was highly specific to X. phaseoli but up to eight tandem repeat loci could be amplified from other Xanthomonas spp. Finally, we assessed the utility of the scheme for analyses of X. phaseoli pv. manihotis genetic variability in the Colombian Caribbean region. MLVA-15 distinguished 88.9% of the haplotypes in our sample. Strains originating from the same field and isolated at the same time could be discriminated. In this study, the advantages of the MLVA-15 scheme targeting 6- or 7-bp repeats were demonstrated. Moreover, this scheme was a fast method that was appropriate for routine monitoring of X. phaseoli pv. manihotis populations on a local scale and, thus, was useful for addressing epidemiological questions.

Rapid and Simple Universal Escherichia coli Genotyping Method Based on Multiple-Locus Variable-Number Tandem-Repeat Analysis Using Single-Tube Multiplex PCR and Standard Gel Electrophoresis

We developed a multiplex PCR method based on multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) that was designed for the rapid typing of Escherichia coli and Shigella isolates. The method amplifies seven VNTRs and does not require a sequencing capillary or fluorescent dyes. The amplification products are simply loaded on a standard agarose gel for electrophoresis, and the banding patterns are analyzed visually. We evaluated the method on 220 strains belonging to different collections: the E. coli reference (ECOR) collection (n = 72), O1:K1 isolates causing neonatal meningitis (n = 38), extended-spectrum beta-lactamase-producing fecal isolates belonging to the worldwide sequence type 131 (ST131) clone (n = 38), Shiga toxin-producing E. coli (STEC) isolates of serogroups O157:H7 (n = 21) and O26 (n = 16, 8 of which belonged to an outbreak), 27 Shigella isolates (22 Shigella sonnei isolates, including 5 epidemic strains), and 8 reference strains. The performances were compared to those of multilocus sequence typing (MLST), the DiversiLab automated repetitive element palindromic PCR (REP-PCR), pulsed-field gel electrophoresis (PFGE), and whole-genome sequencing (WGS). We found 66 different profiles among the isolates in the ECOR collection. Among the clonal group O1:K1 isolates, 14 different profiles were identified. For the 37 STEC isolates, we found 23 profiles, with 1 corresponding to the 8 epidemic strains. We found 19 profiles among the 27 Shigella isolates, with 1 corresponding to the epidemic strain. The method was able to recognize strains of the ST131 clone and to distinguish the O16 and O25b serogroups and identified 15 different MLVA types among them. This method allows the simple, fast, and inexpensive typing of E. coli/Shigella isolates that can be carried out in any laboratory equipped for molecular biology and has a discriminatory power superior to that of MLST and DiversiLab REP-PCR but slightly lower than that of PFGE.IMPORTANCE Fast typing methods that can easily and accurately distinguish clonal groups and unrelated isolates are of particular interest for microbiologists confronted with outbreaks or performing epidemiological studies. Highly discriminatory universal methods, like PFGE, optical mapping, or WGS, are expensive and/or time-consuming. MLST is useful for phylogeny but is less discriminatory and requires sequencing facilities. PCR methods, which are fast and easy to perform, also have drawbacks. Random PCRs and REP-PCR are universal but lack reproducibility. Other PCR methods may lack the discriminatory power to differentiate isolates during outbreaks. MLVA combines the advantages of PCR methods with a high discriminatory power but in its standard form requires sequencing capillary electrophoresis. The method that we have developed combines the advantages of standard PCR (simple, fast, and inexpensive) with the high discriminatory power of MLVA and permits the typing of all E. coli isolates (either intestinal or extraintestinal pathogenic isolates as well as commensal isolates).

Selection of marker genes for genetic barcoding of microorganisms and binning of metagenomic reads by Barcoder software tools

Background

Metagenomic approaches have revealed the complexity of environmental microbiomes with the advancement in whole genome sequencing displaying a significant level of genetic heterogeneity on the species level. It has become apparent that patterns of superior bioactivity of bacteria applicable in biotechnology as well as the enhanced virulence of pathogens often requires distinguishing between closely related species or sub-species. Current methods for binning of metagenomic reads usually do not allow for identification below the genus level and generally stops at the family level.

Results

In this work, an attempt was made to improve metagenomic binning resolution by creating genome specific barcodes based on the core and accessory genomes. This protocol was implemented in novel software tools available for use and download from http://bargene.bi.up.ac.za/. The most abundant barcode genes from the core genomes were found to encode for ribosomal proteins, certain central metabolic genes and ABC transporters. Performance of metabarcode sequences created by this package was evaluated using artificially generated and publically available metagenomic datasets. Furthermore, a program (Barcoding 2.0) was developed to align reads against barcode sequences and thereafter calculate various parameters to score the alignments and the individual barcodes. Taxonomic units were identified in metagenomic samples by comparison of the calculated barcode scores to set cut-off values. In this study, it was found that varying sample sizes, i.e. number of reads in a metagenome and metabarcode lengths, had no significant effect on the sensitivity and specificity of the algorithm. Receiver operating characteristics (ROC) curves were calculated for different taxonomic groups based on the results of identification of the corresponding genomes in artificial metagenomic datasets. The reliability of distinguishing between species of the same genus or family by the program was nearly perfect.

Conclusions

The results showed that the novel online tool BarcodeGenerator (http://bargene.bi.up.ac.za/) is an efficient approach for generating barcode sequences from a set of complete genomes provided by users. Another program, Barcoder 2.0 is available from the same resource to enable an efficient and practical use of metabarcodes for visualization of the distribution of organisms of interest in environmental and clinical samples.

Current Bioinformatics resources in combating infectious diseases

Bioinformatics tools and techniques analyzing next-generation sequencing (NGS) data are increasingly used for the diagnosis and monitoring of infectious diseases. It is of interest to review the application of bioinformatics tools, commonly used databases and NGS data in clinical microbiology, focusing on molecular identification, genotypic, microbiome research, antimicrobial resistance analysis and detection of unknown disease-associated pathogens in clinical specimens. This review documents available bioinformatics resources and databases that are used by medical microbiology scientists and physicians to control emerging infectious pathogens.

Genetic studies to re-affiliate Edwardsiella tarda fish isolates to Edwardsiella piscicida and Edwardsiella anguillarum species

Until 2012, the genus Edwardsiella was composed by three species Edwardsiella tarda, Edwardsiella hoshinae and Edwardsiella ictaluri. In 2013, Edwardsiella piscicida, compiling fish pathogenic strains previously identified as E. tarda was described, and more recently a new species isolated from diseased eel was reported, namely Edwardsiella anguillarum. The incorporation of these species into the genus makes necessary a revision of the taxonomic position of the isolates previously identified as E. tarda. Using AFLP technique, MLSA studies and in silico DNA-DNA hybridization, 46 of 49 E. tarda isolates were re-assigned as E. piscicida and 2 as E. anguillarum, whereas it was confirmed previous classification of the Edwardsiella types and reference strains used. The study of the taxonomic resolution of the genes 16S rRNA, adk, atpD, dnaJ, glnA, hsp60, tuf as well as the possible combinations among housekeeping genes, showed that the gene dnaJ was the more resolutive. In conclusion, the use of molecular techniques is necessary to accurately identify Edwardsiella isolates, especially when differentiating new species from E. tarda.

Development and evaluation of a bioinformatics approach for designing molecular assays for viral detection

BACKGROUND

Viruses belonging to the Flaviviridae and Bunyaviridae families show considerable genetic diversity. However, this diversity is not necessarily taken into account when developing diagnostic assays, which are often based on the pairwise alignment of a limited number of sequences. Our objective was to develop and evaluate a bioinformatics workflow addressing two recurrent issues of molecular assay design: (i) the high intraspecies genetic diversity in viruses and (ii) the potential for cross-reactivity with close relatives.

METHODOLOGY

The workflow developed herein was based on two consecutive BLASTn steps; the first was utilized to select highly conserved regions among the viral taxon of interest, and the second was employed to assess the degree of similarity of these highly-conserved regions to close relatives. Subsequently, the workflow was tested on a set of eight viral species, including various strains from the Flaviviridae and Bunyaviridae families.

PRINCIPAL FINDINGS

The genetic diversity ranges from as low as 0.45% variable sites over the complete genome of the Japanese encephalitis virus to more than 16% of variable sites on segment L of the Crimean-Congo hemorrhagic fever virus. Our proposed bioinformatics workflow allowed the selection-based on computing scores-of the best target for a diagnostic molecular assay for the eight viral species investigated.

CONCLUSIONS/SIGNIFICANCE

Our bioinformatics workflow allowed rapid selection of highly conserved and specific genomic fragments among the investigated viruses, while considering up to several hundred complete genomic sequences. The pertinence of this workflow will increase in parallel to the number of sequences made publicly available. We hypothesize that our workflow might be utilized to select diagnostic molecular markers for higher organisms with more complex genomes, provided the sequences are made available.

A machine learning approach for viral genome classification

Background

Advances in cloning and sequencing technology are yielding a massive number of viral genomes. The classification and annotation of these genomes constitute important assets in the discovery of genomic variability, taxonomic characteristics and disease mechanisms. Existing classification methods are often designed for specific well-studied family of viruses. Thus, the viral comparative genomic studies could benefit from more generic, fast and accurate tools for classifying and typing newly sequenced strains of diverse virus families.

Results

Here, we introduce a virus classification platform, CASTOR, based on machine learning methods. CASTOR is inspired by a well-known technique in molecular biology: restriction fragment length polymorphism (RFLP). It simulates, in silico, the restriction digestion of genomic material by different enzymes into fragments. It uses two metrics to construct feature vectors for machine learning algorithms in the classification step. We benchmark CASTOR for the classification of distinct datasets of human papillomaviruses (HPV), hepatitis B viruses (HBV) and human immunodeficiency viruses type 1 (HIV-1). Results reveal true positive rates of 99%, 99% and 98% for HPV Alpha species, HBV genotyping and HIV-1 M subtyping, respectively. Furthermore, CASTOR shows a competitive performance compared to well-known HIV-1 specific classifiers (REGA and COMET) on whole genomes and pol fragments.

Conclusion

The performance of CASTOR, its genericity and robustness could permit to perform novel and accurate large scale virus studies. The CASTOR web platform provides an open access, collaborative and reproducible machine learning classifiers. CASTOR can be accessed at http://castor.bioinfo.uqam.ca.

Electronic supplementary material

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

Pangenome and immuno-proteomics analysis of Acinetobacter baumannii strains revealed the core peptide vaccine targets

BACKGROUND

Acinetobacter baumannii has emerged as a significant nosocomial pathogen during the last few years, exhibiting resistance to almost all major classes of antibiotics. Alternative treatment options such as vaccines tend to be most promising and cost effective approaches against this resistant pathogen. In the current study, we have explored the pan-genome of A. baumannii followed by immune-proteomics and reverse vaccinology approaches to identify potential core vaccine targets.

RESULTS

The pan-genome of all available A. baumannii strains (30 complete genomes) is estimated to contain 7,606 gene families and the core genome consists of 2,445 gene families (~32 % of the pan-genome). Phylogenetic tree, comparative genomic and proteomic analysis revealed both intra- and inter genomic similarities and evolutionary relationships. Among the conserved core genome, thirteen proteins, including P pilus assembly protein, pili assembly chaperone, AdeK, PonA, OmpA, general secretion pathway protein D, FhuE receptor, Type VI secretion system OmpA/MotB, TonB dependent siderophore receptor, general secretion pathway protein D, outer membrane protein, peptidoglycan associated lipoprotein and peptidyl-prolyl cis-trans isomerase are identified as highly antigenic. Epitope mapping of the target proteins revealed the presence of antigenic surface exposed 9-mer T-cell epitopes. Protein-protein interaction and functional annotation have shown their involvement in significant biological and molecular processes. The pipeline is validated by predicting already known immunogenic targets against Gram negative pathogen Helicobacter pylori as a positive control.

CONCLUSION

The study, based upon combinatorial approach of pan-genomics, core genomics, proteomics and reverse vaccinology led us to find out potential vaccine candidates against A. baumannii. The comprehensive analysis of all the completely sequenced genomes revealed thirteen putative antigens which could elicit substantial immune response. The integration of computational vaccinology strategies would facilitate in tackling the rapid dissemination of resistant A.baumannii strains. The scarcity of effective antibiotics and the global expansion of sequencing data making this approach desirable in the development of effective vaccines against A. baumannii and other bacterial pathogens.

Whole Genome Sequencing Investigation of a Tuberculosis Outbreak in Port-au-Prince, Haiti Caused by a Strain with a "Low-Level" rpoB Mutation L511P - Insights into a Mechanism of Resistance Escalation

The World Health Organization recommends diagnosing Multidrug-Resistant Tuberculosis (MDR-TB) in high burden countries by detection of mutations in Rifampin (RIF) Resistance Determining Region of Mycobacterium tuberculosis rpoB gene with rapid molecular tests GeneXpert MTB/RIF and Hain MTBDRplus. Such mutations are found in >95% of Mycobacterium tuberculosis strains resistant to RIF by conventional culture-based drug susceptibility testing (DST). However routine diagnostic screening with molecular tests uncovered specific “low level” rpoB mutations conferring resistance to RIF below the critical concentration of 1 μg/ml in some phenotypically susceptible strains. Cases with discrepant phenotypic (susceptible) and genotypic (resistant) results for resistance to RIF account for at least 10% of resistant diagnoses by molecular tests and urgently require new guidelines to inform therapeutic decision making. Eight strains with a “low level” rpoB mutation L511P were isolated by GHESKIO laboratory between 2008 and 2012 from 6 HIV-negative and 2 HIV-positive patients during routine molecular testing. Five isolates with a single L511P mutation and two isolates with double mutation L511P&M515T had MICs for RIF between 0.125 and 0.5 μg/ml and tested susceptible in culture-based DST. The eighth isolate carried a double mutation L511P&D516C and was phenotypically resistant to RIF. All eight strains shared the same spoligotype SIT 53 commonly found in Haiti but classic epidemiological investigation failed to uncover direct contacts between the patients. Whole Genome Sequencing (WGS) revealed that L511P cluster isolates resulted from a clonal expansion of an ancestral strain resistant to Isoniazid and to a very low level of RIF. Under the selective pressure of RIF-based therapy the strain acquired mutation in the M306 codon of embB followed by secondary mutations in rpoB and escalation of resistance level. This scenario highlights the importance of subcritical resistance to RIF for both clinical management of patients and public health and provides support for introducing rpoB mutations as proxy for MICs into laboratory diagnosis of RIF resistance. This study illustrates that WGS is a promising multi-purpose genotyping tool for high-burden settings as it provides both “gold standard” sequencing results for prediction of drug susceptibility and a high-resolution data for epidemiological investigation in a single assay.

Practical benefits of knowing the enemy: modern molecular tools for diagnosing the etiology of bacterial diseases and understanding the taxonomy and diversity of plant-pathogenic bacteria

Knowing the identity of bacterial plant pathogens is essential to strategic and sustainable disease management in agricultural systems. This knowledge is critical for growers, diagnosticians, extension agents, and others dealing with crops. However, such identifications are linked to bacterial taxonomy, a complicated and changing discipline that depends on methods and information that are often not used by those who are diagnosing field problems. Modern molecular tools for fingerprinting and sequencing allow for pathogen identification in the absence of distinguishing or conveniently tested phenotypic characteristics. These methods are also useful in studying the etiology and epidemiology of phytopathogenic bacteria from epidemics, as was done in numerous studies conducted in California's Salinas Valley. Multilocus and whole-genome sequence analyses are becoming the cornerstones of studies of microbial diversity and bacterial taxonomy. Whole-genome sequence analysis needs to become adequately accessible, automated, and affordable in order to be used routinely for identification and epidemiology. The power of molecular tools in accurately identifying bacterial pathogenesis is therefore of value to the farmer, diagnostician, phytobacteriologist, and taxonomist.

Multilocus sequence analysis of Anaplasma phagocytophilum reveals three distinct lineages with different host ranges in clinically ill French cattle

Molecular epidemiology represents a powerful approach to elucidate the complex epidemiological cycles of multi-host pathogens, such as Anaplasma phagocytophilum. A. phagocytophilum is a tick-borne bacterium that affects a wide range of wild and domesticated animals. Here, we characterized its genetic diversity in populations of French cattle; we then compared the observed genotypes with those found in horses, dogs, and roe deer to determine whether genotypes of A. phagocytophilum are shared among different hosts. We sampled 120 domesticated animals (104 cattle, 13 horses, and 3 dogs) and 40 wild animals (roe deer) and used multilocus sequence analysis on nine loci (ankA, msp4, groESL, typA, pled, gyrA, recG, polA, and an intergenic region) to characterize the genotypes of A. phagocytophilum present. Phylogenic analysis revealed three genetic clusters of bacterial variants in domesticated animals. The two principal clusters included 98% of the bacterial genotypes found in cattle, which were only distantly related to those in roe deer. One cluster comprised only cattle genotypes, while the second contained genotypes from cattle, horses, and dogs. The third contained all roe deer genotypes and three cattle genotypes. Geographical factors could not explain this clustering pattern. These results suggest that roe deer do not contribute to the spread of A. phagocytophilum in cattle in France. Further studies should explore if these different clusters are associated with differing disease severity in domesticated hosts. Additionally, it remains to be seen if the three clusters of A. phagocytophilum genotypes in cattle correspond to distinct epidemiological cycles, potentially involving different reservoir hosts.

A MLVA genotyping scheme for global surveillance of the citrus pathogen Xanthomonas citri pv. citri suggests a worldwide geographical expansion of a single genetic lineage

MultiLocus Variable number of tandem repeat Analysis (MLVA) has been extensively used to examine epidemiological and evolutionary issues on monomorphic human pathogenic bacteria, but not on bacterial plant pathogens of agricultural importance albeit such tools would improve our understanding of their epidemiology, as well as of the history of epidemics on a global scale. Xanthomonas citri pv. citri is a quarantine organism in several countries and a major threat for the citrus industry worldwide. We screened the genomes of Xanthomonas citri pv. citri strain IAPAR 306 and of phylogenetically related xanthomonads for tandem repeats. From these in silico data, an optimized MLVA scheme was developed to assess the global diversity of this monomorphic bacterium. Thirty-one minisatellite loci (MLVA-31) were selected to assess the genetic structure of 129 strains representative of the worldwide pathological and genetic diversity of X. citri pv. citri. Based on Discriminant Analysis of Principal Components (DAPC), four pathotype-specific clusters were defined. DAPC cluster 1 comprised strains that were implicated in the major geographical expansion of X. citri pv. citri during the 20th century. A subset of 12 loci (MLVA-12) resolved 89% of the total diversity and matched the genetic structure revealed by MLVA-31. MLVA-12 is proposed for routine epidemiological identification of X. citri pv. citri, whereas MLVA-31 is proposed for phylogenetic and population genetics studies. MLVA-31 represents an opportunity for international X. citri pv. citri genotyping and data sharing. The MLVA-31 data generated in this study was deposited in the Xanthomonas citri genotyping database (http://www.biopred.net/MLVA/).

Polymorphisms and ambiguous sites present in DNA sequences of Leishmania clones: looking closer

In genetic studies of Leishmania parasites, co-dominant markers are chosen for their ability to detect heterozygous polymorphisms, to infer the occurrence of inbreeding and to resolve genetic variability. The majority of DNA sequence based reports perform conventional dye terminator cycle sequencing where perfectly ambiguous sites or double peaks in the chromatogram are interpreted as heterozygous strains. However, molecular peculiarities of the parasite such as aneuploidy, mixed populations and homologous recombination advise that data from regular DNA sequence analysis should be carefully evaluated. We report here a closer look at ambiguous sites observed in 6pgd DNA sequences obtained for a multilocus sequence analysis project on Leishmania (Viannia) strains. After comparing 286 DNA sequences from biological and molecular clones of six L. (Viannia) strains we could distinguish events that contribute to genetic variation in Leishmania (recombination, mutation, chromosomal mosaics). Also, the results suggest how diversity might not be completely revealed through regular DNA sequence analysis and demonstrate the importance for molecular epidemiology research to be aware of such possibilities while choosing samples for studies.

Clustering acinetobacter strains by optical mapping

Optical mapping is a technique that produces an ordered restriction map of a bacterial or eukaryotic chromosome. We have developed a new method, the BOP method, to compare experimental optical maps with in silico optical maps of complete genomes to infer the presence/absence of short DNA sequences (bops) in each genome. The BOP method, as implemented by the Optical Mapping suite of four programs, circumvents the necessity of whole-genome multiple alignments and permits reliable strain typing and clustering on the basis of optical maps. We have applied the Optical Mapping Suite to 125 strains of Acinetobacter sp., including 11 completely sequenced genomes and 114 Acinetobacter complex from three US military hospitals. We found that optical mapping completely resolves all 125 strains. Signal to noise analysis showed that when the 125 strains were considered together almost 1/3 of the experimental fragments were misidentified. We found that the set of 125 genomes could be divided into three clusters, two of which included sequenced genomes. Signal to noise analysis after clustering showed that only 3.5% of the experimental restriction fragments were misidentified. Minimum spanning trees of the two clusters that included sequenced genomes are presented. The programs we have developed provide a more rigorous approach for analyzing optical map data than previously existed.

Methicillin-resistant Staphylococcus aureus sequence type 239-III, Ohio, USA, 2007-2009

Identification of virulent strains emphasizes the need for molecular surveillance.

High-throughput single-nucleotide polymorphism-based typing of shared Pseudomonas aeruginosa strains in cystic fibrosis patients using the Sequenom iPLEX platform

Shared strains of Pseudomonas aeruginosa are now well recognized in people with cystic fibrosis (CF), and suitable P. aeruginosa laboratory typing tools are pivotal to understanding their clinical significance and guiding infection control policies in CF clinics. We therefore compared a single-nucleotide polymorphism (SNP)-based typing method using Sequenom iPLEX matrix-assisted laser desorption ionization with time-of-flight mass spectrometry (MALDI-TOF MS) with typing methods used routinely by our laboratory. We analysed 617 P. aeruginosa isolates that included 561 isolates from CF patients collected between 2001 and 2009 in two Brisbane CF clinics and typed previously by enterobacterial repetitive intergenic consensus (ERIC)-PCR, as well as 56 isolates from non-CF patients analysed previously by multilocus sequence typing (MLST). The isolates were tested using a P. aeruginosa Sequenom iPLEX MALDI-TOF (PA iPLEX) method comprising two multiplex reactions, a 13-plex and an 8-plex, to characterize 20 SNPs from the P. aeruginosa housekeeping genes acsA, aroE, guaA, mutL, nuoD, ppsA and trpE. These 20 SNPs were employed previously in a real-time format involving 20 separate assays in our laboratory. The SNP analysis revealed 121 different SNP profiles for the 561 CF isolates. Overall, there was at least 96% agreement between the ERIC-PCR and SNP analyses for all predominant shared strains among patients attending our CF clinics: AUST-01, AUST-02 and AUST-06. For the less frequently encountered shared strain AUST-07, 6/25 (24%) ERIC-PCR profiles were misidentified initially as AUST-02 or as unique, illustrating the difficulty of gel-based analyses. SNP results for the 56 non-CF isolates were consistent with previous MLST data. Thus, the PA iPLEX format provides an attractive high-throughput alternative to ERIC-PCR for large-scale investigations of shared P. aeruginosa strains.

Automated extraction of typing information for bacterial pathogens from whole genome sequence data: Neisseria meningitidis as an exemplar

Whole genome sequence (WGS) data are increasingly used to characterise bacterial pathogens. These data provide detailed information on the genotypes and likely phenotypes of aetiological agents, enabling the relationships of samples from potential disease outbreaks to be established precisely. However, the generation of increasing quantities of sequence data does not, in itself, resolve the problems that many microbiological typing methods have addressed over the last 100 years or so; indeed, providing large volumes of unstructured data can confuse rather than resolve these issues. Here we review the nascent field of storage of WGS data for clinical application and show how curated sequence-based typing schemes on websites have generated an infrastructure that can exploit WGS for bacterial typing efficiently. We review the tools that have been implemented within the PubMLST website to extract clinically useful, strain-characterisation information that can be provided to physicians and public health professionals in a timely, concise and understandable way. These data can be used to inform medical decisions such as how to treat a patient, whether to instigate public health action, and what action might be appropriate. The information is compatible both with previous sequence-based typing data and also with data obtained in the absence of WGS, providing a flexible infrastructure for WGS-based clinical microbiology.

Molecular characterization of Escherichia coli O157:H7 strains isolated from different sources and geographic regions

Escherichia (E.) coli serotype O157:H7 is a globally distributed human enteropathogen and is comprised of microorganisms with closely related genotypes. The main reservoir for this group is bovine bowels, and infection mainly occurs after ingestion of contaminated water and food. Virulence genetic markers of 28 O157:H7 strains were investigated and multilocus enzyme electrophoresis (MLEE) was used to evaluate the clonal structure. O157:H7 strains from several countries were isolated from food, human and bovine feces. According to MLEE, O157:H7 strains clustered into two main clonal groups designated A and B. Subcluster A1 included 82% of the O157:H7 strains exhibiting identical MLEE pattern. Most enterohemorrhagic E. coli (EHEC) O157:H7 strains from Brazil and Argentina were in the same MLEE subgroup. Bovine and food strains carried virulence genes associated with EHEC pathogenicity in humans.

PHYLOViZ: phylogenetic inference and data visualization for sequence based typing methods

Background

With the decrease of DNA sequencing costs, sequence-based typing methods are rapidly becoming the gold standard for epidemiological surveillance. These methods provide reproducible and comparable results needed for a global scale bacterial population analysis, while retaining their usefulness for local epidemiological surveys. Online databases that collect the generated allelic profiles and associated epidemiological data are available but this wealth of data remains underused and are frequently poorly annotated since no user-friendly tool exists to analyze and explore it.

Results

PHYLOViZ is platform independent Java software that allows the integrated analysis of sequence-based typing methods, including SNP data generated from whole genome sequence approaches, and associated epidemiological data. goeBURST and its Minimum Spanning Tree expansion are used for visualizing the possible evolutionary relationships between isolates. The results can be displayed as an annotated graph overlaying the query results of any other epidemiological data available.

Conclusions

PHYLOViZ is a user-friendly software that allows the combined analysis of multiple data sources for microbial epidemiological and population studies. It is freely available at http://www.phyloviz.net.

Suspicion of Mycobacterium avium subsp. paratuberculosis transmission between cattle and wild-living red deer (Cervus elaphus) by multitarget genotyping

Multitarget genotyping of the etiologic agent Mycobacterium avium subsp. paratuberculosis is necessary for epidemiological tracing of paratuberculosis (Johne's disease). The study was undertaken to assess the informative value of different typing techniques and individual genome markers by investigation of M. avium subsp. paratuberculosis transmission between wild-living red deer and farmed cattle with known shared habitats. Fifty-three M. avium subsp. paratuberculosis type II isolates were differentiated by short sequence repeat analysis (SSR; 4 loci), mycobacterial interspersed repetitive-unit-variable-number tandem-repeat analysis (MIRU-VNTR; 8 loci), and restriction fragment length polymorphism analysis based on IS900 (IS900-RFLP) using BstEII and PstI digestion. Isolates originated from free-living red deer (Cervus elaphus) from Eifel National Park (n = 13), six cattle herds living in the area of this park (n = 23), and five cattle herds without any contact with these red deer (n = 17). Data based on individual herds and genotypes verified that SSR G2 repeats did not exhibit sufficient stability for epidemiological studies. Two common SSR profiles (without G2 repeats), nine MIRU-VNTR patterns, and nine IS900-RFLP patterns were detected, resulting in 17 genotypes when combined. A high genetic variability was found for red deer and cattle isolates within and outside Eifel National Park, but it was revealed only by combination of different typing techniques. Results imply that within this restricted area, wild-living and farmed animals maintain a reservoir for specific M. avium subsp. paratuberculosis genotypes. No host relation of genotypes was obtained. Results suggested that four genotypes had been transmitted between and within species and that one genotype had been transmitted between cattle herds only. Use of multitarget genotyping for M. avium subsp. paratuberculosis type II strains and sufficiently stable genetic markers is essential for reliable interpretations of epidemiological studies on paratuberculosis.

Adjusted Wallace coefficient as a measure of congruence between typing methods

We propose a new coefficient, the adjusted Wallace coefficient (AW), and corresponding confidence intervals (CI) as quantitative measures of congruence between typing methods. The performance of the derived CI was evaluated using simulated data. Published microbial typing data were used to demonstrate the advantages of AW over the Wallace coefficient.

Comparative sequence analysis of recA gene among Vibrio cholerae isolates from Iran with globally reported sequences

AIMS

To study the genetic relatedness between V. cholerae isolates from Iran and other countries based on housekeeping gene recA sequence analysis.

METHODS AND RESULTS

A 995-bp region of the recA gene from 24 V. cholerae isolates obtained from human and surface water origins in Iran over a 5-year period was sequenced and compared with the sequence data from the isolates belonging to other places. Cluster analysis of the constructed dendrogram based on recA sequence divergence for our clinical isolates showed one sequence type (ST), whereas environmental isolates revealed eight STs. Interestingly, one of our environmental isolates was intermixed with clinical isolates in the largest cluster containing the epidemic strains. Our 24 isolates plus 198 global isolates available in the GenBank showed 77 sequence types (STs) with at least one nucleotide difference.

CONCLUSIONS

Our result suggested that recA sequencing is a reliable analysis method for understanding the relatedness of the local isolates with the isolates obtained elsewhere.

SIGNIFICANCE AND IMPACT OF THE STUDY

Understanding the genetic relatedness between V. cholerae isolates could give insights into the health care system for better control and prevention of the cholera.

Horizontal gene transfer and assortative recombination within the Acinetobacter baumannii clinical population provide genetic diversity at the single carO gene, encoding a major outer membrane protein channel

We described previously the presence in Acinetobacter baumannii of a novel outer membrane (OM) protein, CarO, which functions as an L-ornithine OM channel and whose loss was concomitant with increased carbapenem resistance among clonally related nosocomial isolates of this opportunistic pathogen. Here, we describe the existence of extensive genetic diversity at the carO gene within the A. baumannii clinical population. The systematic analysis of carO sequences from A. baumannii isolates obtained from public hospitals in Argentina revealed the existence of four highly polymorphic carO variants among them. Sequence polymorphism between the different A. baumannii CarO variants was concentrated in three well-defined protein regions that superimposed mostly to predicted surface-exposed loops. Polymorphism among A. baumannii CarO variants was manifested in differential electrophoretic mobilities, antigenic properties, abilities to form stable oligomeric structures, and l-ornithine influx abilities through the A. baumannii OM under in vivo conditions. Incongruence between the phylogenies of the clinical A. baumannii isolates analyzed and those of the carO variants they harbor suggests the existence of assortative (entire-gene) carO recombinational exchange within the A. baumannii population. Exchange of carO variants possessing differential characteristics mediated by horizontal gene transfer may constitute an A. baumannii population strategy to survive radically changing environmental conditions, such as the leap from inanimate sources to human hosts and vice versa, persistence in a compromised host, and/or survival in health care facilities.

Insertion sequence- and tandem repeat-based genotyping techniques for Xanthomonas citri pv. mangiferaeindicae

Molecular fingerprinting techniques that have the potential to identify or subtype bacteria at the strain level are needed for improving diagnosis and understanding of the epidemiology of pathogens such as Xanthomonas citri pv. mangiferaeindicae, which causes mango bacterial canker disease. We developed a ligation-mediated polymerase chain reaction targeting the IS1595 insertion sequence as a means to differentiate pv. mangiferaeindicae from the closely related pv. anacardii (responsible for cashew bacterial spot), which has the potential to infect mango but not to cause significant disease. This technique produced weakly polymorphic fingerprints composed of ≈70 amplified fragments per strain for a worldwide collection of X. citri pv. mangiferaeindicae but produced no or very weak amplification for pv. anacardii strains. Together, 12 tandem repeat markers were able to subtype X. citri pv. mangiferaeindicae at the strain level, distinguishing 231 haplotypes from a worldwide collection of 299 strains. Multilocus variable number of tandem repeats analysis (MLVA), IS1595-ligation-mediated polymerase chain reaction, and amplified fragment length polymorphism showed differences in discriminatory power and were congruent in describing the diversity of this strain collection, suggesting low levels of recombination. The potential of the MLVA scheme for molecular epidemiology studies of X. citri pv. mangiferaeindicae is discussed.