Published papers

Plasmids are ubiquitous mobile elements that serve as a pool of many host beneficial traits such as antibiotic resistance in bacterial communities. To understand the importance of plasmids in horizontal gene transfer, we need to gain insight into the 'evolutionary history' of these plasmids, i.e. the range of hosts in which they have evolved. Since extensive data support the proposal that foreign DNA acquires the host's nucleotide composition during long-term residence, comparison of nucleotide composition of plasmids and chromosomes could shed light on a plasmid's evolutionary history. The average absolute dinucleotide relative abundance difference, termed δ-distance, has been commonly used to measure differences in dinucleotide composition, or 'genomic signature', between bacterial chromosomes and plasmids. Here, we introduce the Mahalanobis distance, which takes into account the variance-covariance structure of the chromosome signatures. We demonstrate that the Mahalanobis distance is better than the δ-distance at measuring genomic signature differences between plasmids and chromosomes of potential hosts. We illustrate the usefulness of this metric for proposing candidate long-term hosts for plasmids, focusing on the virulence plasmids pXO1 from Bacillus anthracis, and pO157 from Escherichia coli O157:H7, as well as the broad host range multi-drug resistance plasmid pB10 from an unknown host. © 2008 The Author(s).

G + C composition at the third codon position (GC3) is widely reported to be correlated with synonymous codon usage bias. However, no quantitative attempt has been made to compare the extent of this correlation among different genomes. Here, we applied Shannon entropy from information theory to measure the degree of GC3 bias and that of synonymous codon usage bias of each gene. The strength of the correlation of GC3 with synonymous codon usage bias, quantified by a correlation coefficient, varied widely among bacterial genomes, ranging from -0.07 to 0.95. Previous analyses suggesting that the relationship between GC3 and synonymous codon usage bias is independent of species are thus inconsistent with the more detailed analyses obtained here for individual species.

We have developed a comprehensive software suite for bioinformatics research of cDNAs; it is aimed at rapid characterization of the features of genes and the proteins they code. Methods implemented include the detection of translation initiation and termination signals, statistical analysis of codon usage, comparative study of amino acid composition, comparative modeling of the structures of product proteins, prediction of alternative splice forms, and metabolic pathway reconstruction.

Multivariate analyses are often used to identify major trends of variation in synonymous codon usage among genes. These analyses need to be performed on properly normalized codon usage data to avoid biases masking this synonymous variation, i.e., gene length, amino acid usage, and codon degeneracy; however, previous studies have failed to do so. In this paper, we demonstrate that the use of alternative normalized data (called 'relative adaptiveness' in the literature) can avoid all these biases and furthermore, can identify more trends of variation among genes, including GC-ending codon usage, GT-ending codon usage, and gene expression level. © 2005 Published by Elsevier B.V. on behalf of the Federation of European Biochemical Societies.

© 2015 Federation of European Biochemical Societies. The ciliate Paramecium bursaria harbors the green-alga Chlorella symbionts. We reassembled the P. bursaria transcriptome to minimize falsely fused transcripts, and investigated amino acid and codon usage using the transcriptome data. Surface proteins preferentially use smaller amino acid residues like cysteine. Unusual synonymous codon and amino acid usage in highly expressed genes can reflect a balance between translational selection and other factors. A correlation of gene expression level with synonymous codon or amino acid usage is emphasized in genes down-regulated in symbiont-bearing cells compared to symbiont-free cells. Our results imply that the selection is associated with P. bursaria-Chlorella symbiosis.

Background: The white shark (Carcharodon carcharias) is a globally distributed,apex predator possessing physical, physiological, and behavioral traits thathave garnered it significant public attention. In addition to interest inthe genetic basis of its form and function, as a representative of theoldest extant jawed vertebrate lineage, white sharks are also ofconservation concern due to their small population size and threat fromoverfishing. Despite this, surprisingly little is known about the biology ofwhite sharks, and genomic resources are unavailable. To address thisdeficit, we combined Roche-454 and Illumina sequencing technologies tocharacterize the first transciptome of any tissue for this species. Results: From white shark heart cDNA we generated 665,399 Roche 454 reads (medianlength 387-bp) that were assembled into 141,626 contigs (mean length503-bp). We also generated 78,566,588 Illumina reads, which we aligned tothe 454 contigs producing 105,014 454/Illumina consensus sequences. Tothese, we added 3,432 non-singleton 454 contigs. By comparing thesesequences to the UniProtKB/Swiss-Prot database we were able to annotate21,019 translated open reading frames (ORFs) of ≥ 20 amino acids. Ofthese, 19,277 were additionally assigned Gene Ontology (GO) functionalannotations. While acknowledging the limitations of our single tissuetranscriptome, Fisher tests showed the white shark transcriptome to besignificantly enriched for numerous metabolic GO terms compared to the zebrafish and human transcriptomes, with white shark showing more similarity tohuman than to zebra fish (i.e. fewer terms were significantly different). Wealso compared the transcriptome to other available elasmobranch sequences,for signatures of positive selection and identified several genes ofputative adaptive significance on the white shark lineage. The white sharktranscriptome also contained 8,404 microsatellites (dinucleotide,trinucleotide, or tetranucleotide motifs ≥ five perfect repeats).Detailed characterization of these microsatellites showed that ORFs withtrinucleotide repeats, were significantly enriched for transcriptionregulatory roles and that trinucleotide frequency within ORFs was lower thanfor a wide range of taxonomic groups including other vertebrates.Conclusion: The white shark heart transcriptome represents a valuable resource for futureelasmobranch functional and comparative genomic studies, as well as forpopulation and other biological studies vital for effective conservation ofthis globally vulnerable species. © 2013 Richards et al.; licensee BioMed Central Ltd.

Lactobacilli are found in a wide variety of habitats. Four species, Lactobacillus crispatus, L. gasseri, L. iners, and L. jensenii, are common and abundant in the human vagina and absent from other habitats. These may be adapted to the vagina and possess characteristics enabling them to thrive in that environment. Furthermore, stable codominance of multiple Lactobacillus species in a single community is infrequently observed. Thus, it is possible that individual vaginal Lactobacillus species possess unique characteristics that confer to them host-specific competitive advantages. We performed comparative functional genomic analyses of representatives of 25 species of Lactobacillus, searching for habitat-specific traits in the genomes of the vaginal lactobacilli. We found that the genomes of the vaginal species were significantly smaller and had significantly lower GC content than those of the nonvaginal species. No protein families were found to be specific to the vaginal species analyzed, but some were either overor underrepresented relative to nonvaginal species. We also found that within the vaginal species, each genome coded for species- specific protein families. Our results suggest that even though the vaginal species show no general signatures of adaptation to the vaginal environment, each species has specific and perhaps unique ways of interacting with its environment, be it the host or other microbes in the community. These findings will serve as a foundation for further exploring the role of lactobacilli in the ecological dynamics of vaginal microbial communities and their ultimate impact on host health. © 2014, American Society for Microbiology.

© 2017 The Author(s). Background: Clostridium difficile is a Gram-positive anaerobe and the leading cause of antibiotic-associated diarrhea worldwide. The emergence of ribotype 027 (RT027) strains is associated with increased incidence of infection and mortality. To further understand the relationship between C. difficile NCKUH-21, a RT027 strain isolated from a patient in Taiwan, and other RT027 strains, we performed whole-genome shotgun sequencing on NCKUH-21 and comparative genomic analyses. Results: The genome size, G+C content, and gene number for the NCKUH-21 strain were determined to be similar to those for other C. difficile strains. The core genome phylogeny indicated that the five RT027 strains R20291, CD196, NCKUH-21, BI1, and 2007855 formed a clade. A pathogenicity locus, tcdR-tcdB-tcdE-orf-tcdA-tcdC, was conserved in the genome. A genomic region highly similar to the Clostridium phage $$\upvarphi$$ φ CD38-2 was present in the NCKUH-21 strain but absent in the other RT027 strains and designated as the prophage $$\upvarphi$$ φ NCKUH-21. The prophage $$\upvarphi$$ φ NCKUH-21 genes were significantly higher in G+C content than the other genes in the NCKUH-21 genome, indicating that the prophage does not match the base composition of the host genome. Conclusions: This is the first whole-genome analysis of a RT027 C. difficile strain isolated from Taiwan. Due to the high identity with $$\upvarphi$$ φ CD38-2, the prophage identified in the NCKUH-21 genome has the potential to regulate toxin production. These results provide important information for understanding the pathogenicity of RT027 C. difficile in Taiwan.

Background: Staphylococcus belongs to the Gram-positive low G + C content group of the Firmicutes division of bacteria. Staphylococcus aureus is an important human and veterinary pathogen that causes a broad spectrum of diseases, and has developed important multidrug resistant forms such as methicillin-resistant S. aureus (MRSA). Staphylococcus simiae was isolated from South American squirrel monkeys in 2000, and is a coagulase-negative bacterium, closely related, and possibly the sister group, to S. aureus. Comparative genomic analyses of closely related bacteria with different phenotypes can provide information relevant to understanding adaptation to host environment and mechanisms of pathogenicity.Results: We determined a Roche/454 draft genome sequence for S. simiae and included it in comparative genomic analyses with 11 other Staphylococcus species including S. aureus. A genome based phylogeny of the genus confirms that S. simiae is the sister group to S. aureus and indicates that the most basal Staphylococcus lineage is Staphylococcus pseudintermedius, followed by Staphylococcus carnosus. Given the primary niche of these two latter taxa, compared to the other species in the genus, this phylogeny suggests that human adaptation evolved after the split of S. carnosus. The two coagulase-positive species (S. aureus and S. pseudintermedius) are not phylogenetically closest but share many virulence factors exclusively, suggesting that these genes were acquired by horizontal transfer. Enrichment in genes related to mobile elements such as prophage in S. aureus relative to S. simiae suggests that pathogenesis in the S. aureus group has developed by gene gain through horizontal transfer, after the split of S. aureus and S. simiae from their common ancestor.Conclusions: Comparative genomic analyses across 12 Staphylococcus species provide hypotheses about lineages in which human adaptation has taken place and contributions of horizontal transfer in pathogenesis. © 2012 Suzuki et al; licensee BioMed Central Ltd.

Comparative genomics of closely related bacterial species with different pathogenesis and host preference can provide a means of identifying the specifics of adaptive differences. Streptococcus dysgalactiae (SD) is comprised of two subspecies: S. dysgalactiae subsp. equisimilis is both a human commensal organism and a human pathogen, and S. dysgalactiae subsp. dysgalactiae is strictly an animal pathogen. Here, we present complete genome sequences for both taxa, with analyses involving other species of Streptococcus but focusing on adaptation in the SD species group. We found little evidence for enrichment in biochemical categories of genes carried by each SD strain, however, differences in the virulence gene repertoire were apparent. Some of the differences could be ascribed to prophage and integrative conjugative elements. We identified approximately 9% of the nonrecombinant core genome to be under positive selection, some of which involved known virulence factors in other bacteria. Analyses of proteomes by pooling data across genes, by biochemical category, clade, or branch, provided evidence for increased rates of evolution in several gene categories, as well as external branches of the tree. Promoters were primarily evolving under purifying selection but with certain categories of genes evolving faster. Many of these fast-evolving categories were the same as those associated with rapid evolution in proteins. Overall, these results suggest that adaptation to changing environments and new hosts in the SD species group has involved the acquisition of key virulence genes along with selection of orthologous protein-coding loci and operon promoters. © 2011 The Author(s).

© 2015 Scaria et al. Background: Clostridium difficile and C. sordellii are two anaerobic, spore forming, gram positive pathogens with a broad host range and the ability to cause lethal infections. Despite strong similarities between the two Clostridial strains, differences in their host tissue preference place C. difficile infections in the gastrointestinal tract and C. sordellii infections in soft tissues. Results: In this study, to improve our understanding of C. sordellii and C. difficile virulence and pathogenesis, we have performed a comparative genomic and phenomic analysis of the two. The global phenomes of C. difficile and C. sordellii were compared using Biolog Phenotype microarrays. When compared to C. difficile, C. sordellii was found to better utilize more complex sources of carbon and nitrogen, including peptides. Phenotype microarray comparison also revealed that C. sordellii was better able to grow in acidic pH conditions. Using next generation sequencing technology, we determined the draft genome of C. sordellii strain 8483 and performed comparative genome analysis with C. difficile and other Clostridial genomes. Comparative genome analysis revealed the presence of several enzymes, including the urease gene cluster, specific to the C. sordellii genome that confer the ability of expanded peptide utilization and survival in acidic pH. Conclusions: The identified phenotypes of C. sordellii might be important in causing wound and vaginal infections respectively. Proteins involved in the metabolic differences between C. sordellii and C. difficile should be targets for further studies aimed at understanding C. difficile and C. sordellii infection site specificity and pathogenesis.

Plasmids of the incompatibility group IncP-1 are important agents of horizontal gene transfer and contribute to the spread of antibiotic resistance and xenobiotic degradation within bacterial communities. Even though some prototype plasmids have been studied in much detail, the diversity of this plasmid group was still greatly underestimated until recently, as only two of the five currently known divergent sub-groups had been described. To further improve our insight into the diversity and evolutionary history of this family of broad-host-range plasmids, we compared the complete nucleotide sequence of a new IncP-1δ plasmid pAKD4 to the genomes of other IncP-1 plasmids. Plasmid pAKD4 was previously isolated by exogenous plasmid isolation from an agricultural soil in Norway. Its 56,803. bp nucleotide sequence shows high similarity in gene sequence and gene order to both plasmids pEST4011 and pIJB1, the only other IncP-1δ plasmids sequenced so far. While all three plasmids have a typical IncP-1 backbone comprising replication, transfer, and stable inheritance/control genes, the low sequence similarity in some regions and presence/absence of some backbone genes compared to other IncP-1 plasmids cluster them in a divergent sub-group. Therefore this study validates the presence of a real IncP-1δ clade with multiple plasmids. Moreover, since both pEST4011 and pIJB1 are missing a portion of their transfer genes, pAKD4 represents the first completely sequenced self-transferable plasmid with a complete IncP-1δ backbone. We therefore propose it to be the prototype IncP-1δ plasmid. © 2009 Elsevier Inc.

Synonymous codon usage varies both between organisms and among genes within a genome, and arises due to differences in G + C content, replication strand skew, or gene expression levels. Correspondence analysis (CA) is widely used to identify major sources of variation in synonymous codon usage among genes and provides a way to identify horizontally transferred or highly expressed genes. Four methods of CA have been developed based on three kinds of input data: absolute codon frequency, relative codon frequency, and relative synonymous codon usage (RSCU) as well as within-group CA (WCA). Although different CA methods have been used in the past, no comprehensive comparative study has been performed to evaluate their effectiveness. Here, the four CA methods were evaluated by applying them to 241 bacterial genome sequences. The results indicate that WCA is more effective than the other three methods in generating axes that reflect variations in synonymous codon usage. Furthermore, WCA reveals sources that were previously unnoticed in some genomes; e.g. synonymous codon usage related to replication strand skew was detected in Rickettsia prowazekii. Though CA based on RSCU is widely used, our evaluation indicates that this method does not perform as well as WCA.

Background: The ciliate Paramecium bursaria harbors several hundred cells of the green-alga Chlorella sp. in their cytoplasm. Irrespective of the mutual relation between P. bursaria and the symbiotic algae, both cells retain the ability to grow without the partner. They can easily reestablish endosymbiosis when put in contact with each other. Consequently, P. bursaria is an excellent model for studying cell-cell interaction and the evolution of eukaryotic cells through secondary endosymbiosis between different protists. Despite the importance of this organism, no genomic resources have been identified for P. bursaria to date. This investigation compared gene expressions through RNA-Seq analysis and de novo transcriptome assembly of symbiont-free and symbiont-bearing host cells.Results: To expedite the process of gene discovery related to the endosymbiosis, we have undertaken Illumina deep sequencing of mRNAs prepared from symbiont-bearing and symbiont-free P. bursaria cells. We assembled the reads de novo to build the transcriptome. Sequencing using Illumina HiSeq2000 platform yielded 232.3 million paired-end sequence reads. Clean reads filtered from the raw reads were assembled into 68,175 contig sequences. Of these, 10,557 representative sequences were retained after removing Chlorella sequences and lowly expressed sequences. Nearly 90% of these transcript sequences were annotated by similarity search against protein databases. We identified differentially expressed genes in the symbiont-bearing P. bursaria cells relative to the symbiont-free cells, including heat shock 70 kDa protein and glutathione S-transferase.Conclusions: This is the first reported comprehensive sequence resource of Paramecium - Chlorella endosymbiosis. Results provide some keys for the elucidation of secondary endosymbiosis in P. bursaria. We identified P. bursaria genes that are differentially expressed in symbiont-bearing and symbiont-free conditions. © 2014 Kodama et al.; licensee BioMed Central Ltd.

© 2013 Lin et al. Reported herein is the complete genome sequence of Brachyspira pilosicoli strain P43/6/78T, isolated from a pig with clinical disease. This sequence will aid in the study of genome-wide comparison among Brachyspira species.

© 2015 Tonomura et al. Here, we report a draft genome sequence of Anaeromyxobacter sp. strain PSR-1, an arsenate-respiring bacterium isolated from arsenic-contaminated soil. It contained three distinct arsenic resistance gene clusters (ars operons), while no respiratory arsenate reductase gene (arr) was identified.

© 2016 Ito et al. Arenibacter sp. strain C-21, isolated from surface marine sediment of Japan, accumulates iodine in the presence of glucose and iodide (I-). We report here the draft genome sequence of this strain to provide insight into the molecular mechanism underlying its iodine-accumulating ability.

© 2015 Suzuki et al. Caedibacter varicaedens is a kappa killer endosymbiont bacterium of the ciliate Paramecium biaurelia. Here, we present the draft genome sequence of C. varicaedens.

© 2015 Ehara et al. Here, we report a draft genome sequence of Geobacter sp. strain OR-1, an arsenate-respiring bacterium isolated from Japanese paddy soil. It contained two distinct arsenic islands, one including genes for a respiratory arsenate reductase (Arr) as well as for arsenic resistance (arsD-arsA-acr3-arsR-arrA-arrB) and the second containing only genes for arsenic resistance.

Gluconobacter thailandicus strain NBRC 3257, isolated from downy cherry (Prunus tomentosa), is a strict aerobic rod-shaped Gram-negative bacterium. Here, we report the features of this organism, together with the draft genome sequence and annotation. The draft genome sequence is composed of 107 contigs for 3,446,046 bp with 56.17% G+C content and contains 3,360 protein-coding genes and 54 RNA genes.