Published papers

© 2019 The Author(s). Background: The rapid identification of lineage remains a challenge in the genotyping of clinical isolates of recombinogenic pathogens. The chromosome of Mycobacterium avium subsp. hominissuis (MAH), an agent of Mycobacterium avium complex (MAC) lung disease, is often mosaic and is composed of chromosomal segments originating from different lineages. This makes it difficult to infer the MAH lineage in a simple experimental set-up. To overcome this difficulty, we sought to identify chromosomal marker genes containing lineage-specific alleles by genome data mining. Results: We conducted genetic population structure analysis, phylogenetic analysis, and a survey of historical recombination using data from 125 global MAH isolates. Six MAH lineages (EA1, EA2, SC1, SC2, SC3, and SC4) were identified in the current dataset. One P-450 gene (locus-tag MAH-0788/MAV-0940) in the recombination-cold region was found to have multiple alleles that could discriminate five lineages. By combining the information about allele type from one additional gene, the six MAH lineages as well as other M. avium subspecies were distinguishable. A recombination-cold region of 116 kb contains an insertion hotspot and is flanked by a mammalian cell-entry protein operon where allelic variants have previously been reported to occur. Hence, we speculate that the acquisition of lineage-or strain-specific insertions has introduced homology breaks in the chromosome, thereby reducing the chance of interlineage recombination. Conclusions: The allele types of the newly identified marker genes can be used to predict major lineages of M. avium. The single nucleotide polymorphism typing approach targeting multiallelic loci in recombination-cold regions will facilitate the epidemiological study of MAC, and may also be useful for equivalent studies of other nontuberculous mycobacteria potentially carrying mosaic genomes.

The relationship between Legionella and protist hosts has a huge impact when considering the infectious risk in humans because it facilitates the long-term replication and survival of Legionella in the environment. The ciliate Paramecium is considered to be a protist host for Legionella in natural environments, but the details of their endosymbiosis are largely unknown. In this study, we determined candidate Legionella pneumophila genes that are likely to be involved in the establishment of endosymbiosis in Paramecium caudatum by comparing the genomes of Legionella spp. and Holospora spp. that are obligate endosymbiotic bacteria in Paramecium spp. Among the candidate genes, each single deletion mutant for five genes (lpg0492, lpg0522, lpg0523, lpg2141 and lpg2398) failed to establish endosymbiosis in P. caudatum despite showing intracellular growth in human macrophages. The mutants exhibited no characteristic changes in terms of their morphology, multiplication rate or capacity for modulating the phagosomes in which they were contained, but their resistance to lysozyme decreased significantly. This study provides insights into novel factors required by L. pneumophila for endosymbiosis in P. caudatum, and suggests that endosymbiotic organisms within conspecific hosts may have shared genes related to effective endosymbiosis establishment.

© 2015 Loftie-Eaton et al. We report the draft genome sequence of Pseudomonas sp. nov. H2, isolated from creek sediment in Moscow, ID, USA. The strain is most closely related to Pseudomonas putida. However, it has a slightly smaller genome that appears to have been impacted by horizontal gene transfer and poorly maintains IncP-1 plasmids.

© 2014 Ehara et al. Here we report the draft genome sequence of strain Q-1, an iodide (I-)-oxidizing heterotrophic bacterium in the class Alphaproteobacteria isolated from natural gas brine water. The genome sequence contained a multicopper oxidase gene probably responsible for iodide oxidation. A photosynthetic gene cluster was found but genes for carbon-fixation were absent.

© 2016 Useh et al. Here, we report the availability of draft genomes of several Salmonella serotypes, isolated from poultry sources from Nigeria. These genomes will help to further understand the biological diversity of S. enterica and will serve as references in microbial trace-back studies to improve food safety.

© 2015 Watanabe et al. Legionella pneumophila is the causative agent of legionellosis. Here, we report the draft genome sequences of five L. pneumophila strains, Bnt314, Ofk308, Twr292, Ymg289, and Ymt294, isolated from environmental water samples. Comparative analyses of these genomes may reveal the survival mechanisms and virulence of L. pneumophila in the natural environment.

© 2014 Federation of European Microbiological Societies. We present draft genome sequences of three Holospora species, hosted by the ciliate Paramecium caudatum; that is, the macronucleus-specific H. obtusa and the micronucleus-specific H. undulata and H. elegans. We investigate functions of orthologous core genes conserved across the three Holospora species, which may be essential for the infection and survival in the host nucleus.

Defining bacterial species and understanding the relative cohesiveness of different components of their genomes remains a fundamental problem in microbiology. Bacterial species tend to be comprised of both a set of core and dispensable genes, with the sum of these two components forming the species pan-genome. The role of the core and dispensable genes in defining bacterial species and the question of whether pan-genomes are finite or infinite remain unclear. Here we demonstrate, through the analysis of 96 genome sequences derived from two closely related sympatric sister species of pathogenic bacteria (Campylobacter coli and C. jejuni), that their pan-genome is indeed finite and that there are unique and cohesive features to each of their genomes defining their genomic identity. The two species have a similar pan-genome size; however, C. coli has acquired a larger core genome and each species has evolved a number of species-specific core genes, possibly reflecting different adaptive strategies. Genome-wide assessment of the level of lateral gene transfer within and between the two sister species, as well as within the core and non-core genes, demonstrates a resistance to interspecies recombination in the core genome of the two species and therefore provides persuasive support for the core genome hypothesis for bacterial species. © The Author(s) 2010.

Rates of nonsynonymous substitution (dN) significantly higher than rates of synonymous substitution (dS) have been used as evidence of positive selection for the fixation of advantageous point mutations. It has been suggested that positive selection contributes to the evolution of virulence factors and certain functional categories in bacterial pathogens. The genus Streptococcus contains a number of important human and agricultural pathogens. Here we assessed positive selection across 13 Streptococcus species, and their relationship with virulence factors and functional categories. We found that known virulence genes were subject to positive selection pressure as much as other genes. After false discovery rate correction for multiple comparisons, no functional categories were significantly over- or under-represented in positively selected genes relative to other genes. Our results suggest that within the genus Streptococcus positive selection based on dN/dS ratios is not distributed with bias across biological functions. © 2012.

Copyright © 2014 Crohn's & Colitis Foundation of America, Inc. Background: Perturbations of the intestinal microbiome, termed dysbiosis, are linked to intestinal inflammation. Isolation of adherent-invasive Escherichia coli (AIEC) from intestines of patients with Crohn's disease (CD), dogs with granulomatous colitis, and mice with acute ileitis suggests these bacteria share pathoadaptive virulence factors that promote inflammation. Methods: To identify genes associated with AIEC, we sequenced the genomes of phylogenetically diverse AIEC strains isolated from people with CD (4), dogs with granulomatous colitis (2), and mice with ileitis (2) and 1 non-AIEC strain from CD ileum and compared them with 38 genome sequences of E. coli and Shigella. We then determined the prevalence of AIEC-associated genes in 49 E. coli strains from patients with CD and controls and correlated genotype with invasion of intestinal epithelial cells, persistence within macrophages, AIEC pathotype, and growth in standardized conditions. Results: Genes encoding propanediol utilization (pdu operon) and iron acquisition (yersiniabactin, chu operon) were overrepresented in AIEC relative to nonpathogenic E. coli. PduC (propanediol dehydratase) was enriched in CD-derived AIEC, correlated with increased cellular invasion, and persistence in vitro and was increasingly expressed in fucose-containing media. Growth of AIEC required iron, and the presence of chuA (heme acquisition) correlated with persistence in macrophages. CD-associated AIEC with lpfA154 (long polar fimbriae) demonstrated increased invasion of epithelial cells and translocation across M cells. Conclusions: Our findings provide novel insights into the genetic basis of the AIEC pathotype, supporting the concept that AIEC are equipped to exploit and promote intestinal inflammation and reveal potential targets for intervention against AIEC and inflammation-associated dysbiosis.

Background: In many bacteria, intragenomic diversity in synonymous codon usage among genes has been reported. However, no quantitative attempt has been made to compare the diversity levels among different genomes. Here, we introduce a mean dissimilarity-based index (Dmean) for quantifying the level of diversity in synonymous codon usage among all genes within a genome. Results: The application of Dmean to 268 bacterial genomes shows that in bacteria with extremely biased genomic G+C compositions there is little diversity in synonymous codon usage among genes. Furthermore, our findings contradict previous reports. For example, a low level of diversity in codon usage among genes has been reported for Helicobacter pylori, but based on Dmean, the diversity level of this species is higher than those of more than half of bacteria tested here. The discrepancies between our findings and previous reports are probably due to differences in the methods used for measuring codon usage diversity. Conclusion: We recommend that Dmean be used to measure the diversity level of codon usage among genes. This measure can be applied to other compositional features such as amino acid usage and dinucleotide relative abundance as a genomic signature. © 2009 Suzuki et al; licensee BioMed Central Ltd.

Bacterial plasmids are extra-chromosomal genetic elements that code for a wide variety of phenotypes in their bacterial hosts and are maintained in bacterial communities through both vertical and horizontal transfer. Current mathematical models of plasmid-bacteria dynamics, based almost exclusively on mass-action differential equations that describe these interactions in completely mixed environments, fail to adequately explain phenomena such as the long-term persistence of plasmids in natural and clinical bacterial communities. This failure is, at least in part, due to the absence of any spatial structure in these models, whereas most bacterial populations are spatially structured in microcolonies and biofilms. To help bridge the gap between theoretical predictions and observed patterns of plasmid spread and persistence, an individual-based lattice model (interacting particle system) that provides a predictive framework for understanding the dynamics of plasmid-bacteria interactions in spatially structured populations is presented here. To assess the accuracy and flexibility of the model, a series of experiments that monitored plasmid loss and horizontal transfer of the IncP-1β plasmid pB10:: rfp in Escherichia coli K12 and other bacterial populations grown on agar surfaces were performed. The model-based visual patterns of plasmid loss and spread, as well as quantitative predictions of the effects of different initial parental strain densities and incubation time on densities of transconjugants formed on a 2D grid, were in agreement with this and previously published empirical data. These results include features of spatially structured populations that are not predicted by mass-action differential equation models. © 2007 SGM.

The self-transmissible, broad-host-range (BHR) plasmid pMOL98 was previously isolated from polluted soil using a triparental plasmid capture approach and shown to possess a replicon similar to that of the BHR plasmids pSB102 and pIPO2. Here, complete sequence analysis and comparative genomics reveal that the 55.5 kb nucleotide sequence of pMOL98 shows extensive sequence similarity and synteny with the BHR plasmid family that now includes pIPO2, pSB102, pTER331, and pMRAD02. They share a plasmid backbone comprising replication, partitioning and conjugative transfer functions. Comparison of the variable accessory regions of these plasmids shows that the majority of natural transposons, as well as the mini-transposon used to mark the plasmids, are inserted in the parA locus. The transposon unique to pMOL98 appears to have inserted from the chromosome of the recipient strain used in the plasmid capture procedure. This demonstrates the necessity for careful screening of plasmids and host chromosomes to avoid mis-interpretation of plasmid genome content. The presence of very similar BHR plasmids with different accessory genes in geographically distinct locations suggests an important role in horizontal gene exchange and bacterial adaptation for this recently defined plasmid group, which we propose to name "PromA". © Springer Science+Business Media B.V. 2009.

Despite the important contribution of self-transmissible plasmids to bacterial evolution, little is understood about the range of hosts in which these plasmids have evolved. Our goal was to infer this so-called evolutionary host range. The nucleotide composition, or genomic signature, of plasmids is often similar to that of the chromosome of their current host, suggesting that plasmids acquire their hosts' signature over time. Therefore, we examined whether the evolutionary host range of plasmids could be inferred by comparing their trinucleotide composition to that of all completely sequenced bacterial chromosomes. The diversity of candidate hosts was determined using taxonomic classification and genetic distance. The method was first tested using plasmids from six incompatibility (Inc) groups whose host ranges are generally thought to be narrow (IncF, IncH, and IncI) or broad (IncN, IncP, and IncW) and then applied to other plasmid groups. The evolutionary host range was found to be broad for IncP plasmids, narrow for IncF and IncI plasmids, and intermediate for IncH and IncN plasmids, which corresponds with their known host range. The IncW plasmids as well as several plasmids from the IncA/C, IncP, IncQ, IncU, and PromA groups have signatures that were not similar to any of the chromosomal signatures, raising the hypothesis that these plasmids have not been ameliorated in any host due to their promiscuous nature. The inferred evolutionary host range of IncA/C, IncP-9, and IncL/M plasmids requires further investigation. In this era of high-throughput sequencing, this genomic signature method is a useful tool for predicting the host range of novel mobile elements. Copyright © 2010, American Society for Microbiology. All Rights Reserved.

Proteins play a critical role in complex biological systems, yet about half of the proteins in publicly available databases are annotated as functionally unknown. Proteome-wide functional classification using bioinformatics approaches thus is becoming an important method for revealing unknown protein functions. Using the hyperthermophilic archaeon Pyrococcus furiosus as a model species, we used the support vector machine (SVM) method to discriminate DNA/RNA-binding proteins from proteins with other functions, using amino acid composition and periodicities as feature vectors. We defined this value as the composition score (CO) and periodicity score (PD). The P. furiosus proteins were classified into three classes (I-III) on the basis of the two-dimensional correlation analysis of CO score and PD score. As a result, approximately 87 of the functionally known proteins categorized as class I proteins (CO score + PD score > 0.6) were found to be DNA/RNA-binding proteins. Applying the two-dimensional correlation analysis to the 994 hypothetical proteins in P. furiosus, a total of 151 proteins were predicted to be novel DNA/RNA-binding protein candidates. DNA/RNA-binding activities of randomly chosen hypothetical proteins were experimentally verified. Six out of seven candidate proteins in class I possessed DNA/RNA-binding activities, supporting the efficacy of our method. © The Author 2007. Kazusa DNA Research Institute.

Background: Due to their bi-directional replication machinery starting from a single finite origin, bacterial genomes show characteristic nucleotide compositional bias between the two replichores, which can be visualised through GC skew or (C-G)/(C+G). Although this polarisation is used for computational prediction of replication origins in many bacterial genomes, the degree of GC skew visibility varies widely among different species, necessitating a quantitative measurement of GC skew strength in order to provide confidence measures for GC skew-based predictions of replication origins.Results: Here we discuss a quantitative index for the measurement of GC skew strength, named the generalised GC skew index (gGCSI), which is applicable to genomes of any length, including bacterial chromosomes and plasmids. We demonstrate that gGCSI is independent of the window size and can thus be used to compare genomes with different sizes, such as bacterial chromosomes and plasmids. It can suggest the existence of different replication mechanisms in archaea and of rolling-circle replication in plasmids. Correlation of gGCSI values between plasmids and their corresponding host chromosomes suggests that within the same strain, these replicons have reproduced using the same replication machinery and thus exhibit similar strengths of replication strand skew.Conclusions: gGCSI can be applied to genomes of any length and thus allows comparative study of replication-related mutation and selection pressures in genomes of different lengths such as bacterial chromosomes and plasmids. Using gGCSI, we showed that replication-related mutation or selection pressure is similar for replicons with similar machinery. © 2009 Arakawa et al; licensee BioMed Central Ltd.

The overall architecture of IncP-1 plasmids is very conserved in that the accessory genes are typically located in one or two specific regions: between oriV and trfA and between the tra and trb operons. Various hypotheses have been formulated to explain this, but none have been tested experimentally. We investigated whether this structural similarity is due to region-specific transposition alone or also is reliant on selection for plasmids with insertions limited to these two regions. We first examined the transposition of Tn21Km into IncP-1β plasmid pBP136 and found that most Tn21Km insertions (67%) were located around oriV. A similar experiment using the oriV region of IncP-1β plasmid pUO1 confirmed these results. We then tested the transferability, stability, and fitness cost of different pBP136 derivatives to determine if impairment of these key plasmid characters explained the conserved plasmid architecture. Most of the pBP136 derivatives with insertions in transfer genes were no longer transferable. The plasmids with insertions in the oriV-trfA and tra-trb regions were more stable than other plasmid variants, and one of these also showed a significantly lower fitness cost. In addition, our detailed sequence analysis of IncP-1 plasmids showed that Tn402/5053-like transposons are situated predominantly between the tra and trb operons and close to the putative resolution site for the ParA resolvase, a potential hot spot for those transposons. Our study presents the first empirical evidence that region-specific insertion of transposons in combination with selection for transferable and stable plasmids explains the structural similarity of IncP-1 plasmids. Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Shannon entropy from information theory has been applied to estimate the degree of deviation from equal usage of synonymous codons; however, previous attempts have failed to take into account all three aspects of amino acid usage, i.e. (i) the number of distinct amino acids, (ii) their relative frequencies, and (iii) their degree of codon degeneracy. A new index taking into account all of these aspects is proposed. The index, designated as the 'weighted sum of relative entropy' (Ew), is defined as the sum of the relative entropy of each amino acid weighted by its relative frequency in the sequence. In this paper, we demonstrate that Ew allows us to avoid some amino acid usage biases and can yield results contradictory to those obtained by previous methods. © 2004 Elsevier B.V. All rights reserved.

To determine whether different racial groups shared common types of vaginal microbiota, we characterized the composition and structure of vaginal bacterial communities in asymptomatic and apparently healthy Japanese women in Tokyo, Japan, and compared them with those of White and Black women from North America. The composition of vaginal communities was compared based on community profiles of terminal restriction fragments of 16S rRNA genes and phylogenetic analysis of cloned 16S rRNA gene sequences of the numerically dominant bacterial populations. The types of vaginal communities found in Japanese women were similar to those of Black and White women. As with White and Black women, most vaginal communities were dominated by lactobacilli, and only four species of Lactobacillus (Lactobacillus iners, Lactobacillus crispatus, Lactobacillus jensenii, and Lactobacillus gasseri) were commonly found. Communities dominated by multiple species of lactobacilli were common in Japanese and White women, but rare in Black women. The incidence, in Japanese women, of vaginal communities with several non-Lactobacillus species at moderately high frequencies was intermediate between Black women and White women. The limited number of community types found among women in different ethnic groups suggests that host genetic factors, including the innate and adaptive immune systems, may be more important in determining the species composition of vaginal bacterial communities than are cultural and behavioral differences. © 2009 Federation of European Microbiological Societies. Published by Blackwell Publishing Ltd. All rights reserved.