Search

Diatoms adapt to various aquatic light environments and play major roles in the global carbon cycle using their unique light-harvesting system, i.e., fucoxanthin chlorophyll a/c binding proteins (FCPs). Structural analyses of photosystem II (PSII)-FCPII and photosystem I (PSI)-FCPI complexes from the diatom Chaetoceros gracilis have revealed the localization and interactions of many FCPs; however, the entire set of FCPs has not been characterized. Here, we identify 46 FCPs in the newly assembled genome and transcriptome of C. gracilis. Phylogenetic analyses suggest that these FCPs can be classified into five subfamilies: Lhcr, Lhcf, Lhcx, Lhcz, and the novel Lhcq, in addition to a distinct type of Lhcr, CgLhcr9. The FCPs in Lhcr, including CgLhcr9 and some Lhcqs, have orthologous proteins in other diatoms, particularly those found in the PSI-FCPI structure. By contrast, the Lhcf subfamily, some of which were found in the PSII-FCPII complex, seems to be diversified in each diatom species, and the number of Lhcqs differs among species, indicating that their diversification may contribute to species-specific adaptations to light. Further phylogenetic analyses of FCPs/light-harvesting complex (LHC) proteins using genome data and assembled transcriptomes of other diatoms and microalgae in public databases suggest that our proposed classification of FCPs is common among various red-lineage algae derived from secondary endosymbiosis of red algae, including Haptophyta. These results provide insights into the loss and gain of FCP/LHC subfamilies during the evolutionary history of the red algal lineage. This article is protected by copyright. All rights reserved.

Diatoms adapt to various aquatic light environments and play major roles in the global carbon cycle using their unique light-harvesting system, i.e., fucoxanthin chlorophyll a/c binding proteins (FCPs). Structural analyses of photosystem II (PSII)-FCPII and photosystem I (PSI)-FCPI complexes from the diatom Chaetoceros gracilis have revealed the localization and interactions of many FCPs; however, the entire set of FCPs has not been characterized. Here, we identified 46 FCPs in the newly assembled genome and transcriptome of C. gracilis. Phylogenetic analyses suggested that these FCPs could be classified into five subfamilies: Lhcr, Lhcf, Lhcx, Lhcz, and novel Lhcq, in addition to a distinct type of Lhcr, CgLhcr9. The FCPs in Lhcr, including CgLhcr9 and some Lhcqs, had orthologous proteins in other diatoms, particularly those found in the PSI-FCPI structure. By contrast, the Lhcf subfamily, some of which were found in the PSII-FCPII complex, seemed to be diversified in each diatom species, and the number of Lhcqs differed among species, indicating that their diversification may contribute to species-specific adaptations to light. Further phylogenetic analyses of FCPs/light-harvesting complex (LHC) proteins using genome data and assembled transcriptomes of other diatoms and microalgae in public databases suggest that our proposed classification of FCPs was common among various red-lineage algae derived from secondary endosymbiosis of red algae, including Haptophyta. These results provided insights into the loss and gain of FCP/LHC subfamilies during the evolutionary history of the red algal lineage.One sentence summaryPhylogenetic analysis of fucoxanthin-chlorophyll a/c proteins in C. gracilis revealed five major subfamilies and one minor subfamily, providing insights into the diversification of light-harvesting systems in red algae.

We present here the second complete genome of anaerobic ammonium oxidation (anammox) bacterium, Candidatus (Ca.) Brocadia pituitae, along with those of a nitrite oxidizer and two incomplete denitrifiers from the anammox bacterial community (ABC) metagenome. Although NO2− reduction to NO is considered to be the first step in anammox, Ca. B. pituitae lacks nitrite reductase genes (nirK and nirS) responsible for this reaction. Comparative genomics of Ca. B. pituitae with Ca. Kuenenia stuttgartiensis and six other anammox bacteria with nearly complete genomes revealed that their core genome structure contains 1,152 syntenic orthologues. But nitrite reductase genes were absent from the core, whereas two other Brocadia species possess nirK and these genes were horizontally acquired from multiple lineages. In contrast, at least five paralogous hydroxylamine oxidoreductase genes containing candidate ones (hao2 and hao3) encoding another nitrite reductase were observed in the core. Indeed, these two genes were also significantly expressed in Ca. B. pituitae as in other anammox bacteria. Because many nirS and nirK genes have been detected in the ABC metagenome, Ca. B. pituitae presumably utilises not only NO supplied by the ABC members but also NO and/or NH2OH by self-production for anammox metabolism.

Photosynthetic light-harvesting complexes (LHCs) play a pivotal role in collecting solar energy for photochemical reactions in photosynthesis. One of the major LHCs are fucoxanthin chlorophyll a/c-binding proteins (FCPs) present in diatoms, a group of organisms having important contribution to the global carbon cycle. Here, we report a 2.40-Å resolution structure of the diatom photosystem I (PSI)-FCPI supercomplex by cryo-electron microscopy. The supercomplex is composed of 16 different FCPI subunits surrounding a monomeric PSI core. Each FCPI subunit showed different protein structures with different pigment contents and binding sites, and they form a complicated pigment–protein network together with the PSI core to harvest and transfer the light energy efficiently. In addition, two unique, previously unidentified subunits were found in the PSI core. The structure provides numerous insights into not only the light-harvesting strategy in diatom PSI-FCPI but also evolutionary dynamics of light harvesters among oxyphototrophs., 珪藻の光化学系I-集光性色素タンパク質複合体の立体構造を解明 --進化すると色素タンパク質が増える--. 京都大学プレスリリース. 2020-05-26.

Founder animals carrying high proportions of somatic mutation induced by CRISPR-Cas9 enable a rapid and scalable strategy for the functional screening of numerous target genes in vivo. In this functional screening, genotyping using pooled amplicons with next-generation sequencing is the most suitable approach for large-scale management of multiple samples and accurate evaluation of the efficiency of Cas9-induced somatic mutations at target sites. Here, we present a simple workflow for genotyping of multiple CRISPR-Cas9-based knockout founders by pooled amplicon sequencing. Using custom barcoded primers, pooled amplicons from multiple individuals can be run in a single-indexed library on the Illumina MiSeq platform. Additionally, a user-friendly web tool, CLiCKAR, is available to simultaneously perform demultiplexing of pooled sequence data and evaluation of somatic mutation in each phenotype. CLiCKAR provides users with practical reports regarding the positions of insertions/deletions, as well as the frameshift ratio and tables containing mutation sequences, and read counts of each phenotype, with just a few clicks by the implementation of demultiplexing for pooled sample data and calculation of the frameshift ratio. This genotyping workflow can be harnessed to evaluate genotype-phenotype correlations in CRISPR-Cas9-based loss-of-function screening of numerous target genes in various organisms.

We report on the activities of the 2015 edition of the BioHackathon, an annual event that brings together researchers and developers from around the world to develop tools and technologies that promote the reusability of biological data. We discuss issues surrounding the representation, publication, integration, mining and reuse of biological data and metadata across a wide range of biomedical data types of relevance for the life sciences, including chemistry, genotypes and phenotypes, orthology and phylogeny, proteomics, genomics, glycomics, and metabolomics. We describe our progress to address ongoing challenges to the reusability and reproducibility of research results, and identify outstanding issues that continue to impede the progress of bioinformatics research. We share our perspective on the state of the art, continued challenges, and goals for future research and development for the life sciences Semantic Web.

The molecular mechanisms of cell reprogramming and differentiation involve various signaling factors. Small molecule compounds have been identified to artificially influence these factors through interacting cellular proteins. Although such small molecule compounds are useful to enhance reprogramming and differentiation and to show the mechanisms that underlie these events, the screening usually requires a large number of compounds to identify only a very small number of hits (e.g., one hit among several tens of thousands of compounds). Here, we show a proof of concept that xenospecific gene products can affect the efficiency of cell reprogramming to pluripotency. Thirty genes specific for the bacterium Wolbachia pipientis were forcibly expressed individually along with reprogramming factors (Oct4, Sox2, Klf4 and c-Myc) that can generate induced pluripotent stem cells in mammalian cells, and eight were found to affect the reprogramming efficiency either positively or negatively (hit rate 26.7%). Mechanistic analysis suggested one of these proteins interacted with cytoskeleton to promote reprogramming. Our results raise the possibility that xenospecific gene products provide an alternative way to study the regulatory mechanism of cell identity.

Nitrification, the microbial oxidation of ammonia to nitrate via nitrite, occurs in a wide range of acidic soils. However, the ammonia-oxidizing bacteria (AOB) that have been isolated from soil to date are acid-sensitive. Here we report the isolation and characterization of an acid-adapted AOB from an acidic agricultural soil. The isolated AOB, strain TAO100, is classified within the Gammaproteobacteria based on phylogenetic characteristics. TAO100 can grow in the pH range of 5-7.5 and survive in highly acidic conditions until pH 2 by forming cell aggregates. Whereas all known gammaproteobacterial AOB (γ-AOB) species, which have been isolated from marine and saline aquatic environments, are halophiles, TAO100 is not phenotypically halophilic. Thus, TAO100 represents the first soil-originated and non-halophilic γ-AOB. The TAO100 genome is considerably smaller than those of other γ-AOB and lacks several genes associated with salt tolerance which are unnecessary for survival in soil. The ammonia monooxygenase subunit A gene of TAO100 and its transcript are higher in abundance than those of ammonia-oxidizing archaea and betaproteobacterial AOB in the strongly acidic soil. These results indicate that TAO100 plays an important role in the nitrification of acidic soils. Based on these results, we propose TAO100 as a novel species of a new genus, Candidatus Nitrosoglobus terrae.

Selection has been a central issue in biology in eukaryotes as well as prokaryotes. Inference of selection in recombining bacterial species, compared with clonal ones, has been a challenge. It is not known how codons under diversifying selection are distributed along the chromosome or among functional categories or how frequently such codons are subject to mutual homologous recombination. Here, we explored these questions by analysing genes present in >90% among 29 genomes of Helicobacter pylori, one of the bacterial species with the highest mutation and recombination rates. By a method for recombining sequences, we identified codons under diversifying selection (dN/dS > 1), which were widely distributed and accounted for ∼0.2% of all the codons of the genome. The codons were enriched in genes of host interaction/cell surface and genome maintenance (DNA replication, recombination, repair, and restriction modification system). The encoded amino acid residues were sometimes found adjacent to critical catalytic/binding residues in protein structures. Furthermore, by estimating the intensity of homologous recombination at a single nucleotide level, we found that these codons appear to be more frequently subject to recombination. We expect that the present study provides a new approach to population genomics of selection in recombining prokaryotes.

Urodele amphibian newts have unique biological properties, notably including prominent regeneration ability. Iberian ribbed newt,Pleurodeles waltl, is a promising model newt along with the successful development of the easy breeding system and efficient transgenic and genome editing methods. However, genetic information ofP. waltlwas limited. In the present study, we conducted an intensive transcriptome analysis ofP. waltlusing RNA-sequencing to build gene models and annotate them. We generated 1.2 billion Illumina reads from a wide variety of samples across 11 different tissues and 9 time points during embryogenesis. They were assembled into 202,788 non-redundant contigs that appear to cover nearly complete (~98%)P. waltlprotein-coding genes. Using the gene set as a reference, our gene network analysis identified regeneration-, developmental-stage-, and tissue-specific co-expressed gene modules. Ortholog analyses with other vertebrates revealed the gene repertoire evolution of amphibians which includes urodele-specific loss ofbmp4and duplications ofwnt11b. Our transcriptome resource will enhance future research employing this emerging model animal for regeneration research as well as other areas such as developmental biology, stem cell biology, cancer research, ethology and toxico-genomics. These data are available via our portal website, iNewt (http://www.nibb.ac.jp/imori/main/).

The Microbial Genome Database for Comparative Analysis (MBGD) is a database for comparative genomics based on comprehensive orthology analysis of bacteria, archaea and unicellular eukaryotes. MBGD now contains 6318 genomes. To utilize the database for both closely related and distantly related genomes, MBGD previously provided two types of ortholog tables: the standard ortholog table containing one representative genome from each genus covering the entire taxonomic range and the taxon specific ortholog tables for each taxon. However, this approach has a drawback in that the standard ortholog table contains only genes that are conserved in the representative genomes. To address this problem, we developed a stepwise procedure to construct ortholog tables hierarchically in a bottom-up manner. By using this approach, the new standard ortholog table now covers the entire gene repertoire stored in MBGD. In addition, we have enhanced several functionalities, including rapid and flexible keyword searching, profile-based sequence searching for orthology assignment to a user query sequence, and displaying a phylogenetic tree of each taxon based on the concatenated core gene sequences. For integrative database searching, the core data in MBGD are represented in Resource Description Framework (RDF) and a SPARQL interface is provided to search them. MBGD is available at http://mbgd.genome.ad.jp/.

Summary Spore formation is considered to be one of the microbial strategies for long-term survival in subseafloor sedimentary habitats. However, our knowledge of the genetic and physiological characteristics of subseafloor microbes is limited. Here, we studied the distribution and frequency of genes that are related to endospore formation in 10 subseafloor sedimentary metagenomes from Site C9001 off Japan and Site 1229 off Peru. None or very low frequencies of endospore-specific genes (e.g. dpaA, dpaB, sspA, spo0A, spoIIGA, spoIIM, spoIIIAB, spoIVA, spoIVB, yabP, yunB, spoVM) were observed in the subseafloor metagenomes. Based on the number of universally conserved single copy genes, the frequency ratio of putative endospore-formers was estimated to be

Comparative genomics is becoming an essential approach for identification of genes associated with a specific function or phenotype. Here, we introduce the microbial genome database for comparative analysis (MBGD), which is a comprehensive ortholog database among the microbial genomes available so far. MBGD contains several precomputed ortholog tables including the standard ortholog table covering the entire taxonomic range and taxon-specific ortholog tables for various major taxa. In addition, MBGD allows the users to create an ortholog table within any specified set of genomes through dynamic calculations. In particular, MBGD has a "My MBGD" mode where users can upload their original genome sequences and incorporate them into orthology analysis. The created ortholog table can serve as the basis for various comparative analyses. Here, we describe the use of MBGD and briefly explain how to utilize the orthology information during comparative genome analysis in combination with the stand-alone comparative genomics software RECOG, focusing on the application to comparison of closely related microbial genomes.

Comparative genomics is becoming an essential approach for identification of genes associated with a specific function or phenotype. Here, we introduce the microbial genome database for comparative analysis (MBGD), which is a comprehensive ortholog database among the microbial genomes available so far. MBGD contains several precomputed ortholog tables including the standard ortholog table covering the entire taxonomic range and taxon-specific ortholog tables for various major taxa. In addition, MBGD allows the users to create an ortholog table within any specified set of genomes through dynamic calculations. In particular, MBGD has a "My MBGD" mode where users can upload their original genome sequences and incorporate them into orthology analysis. The created ortholog table can serve as the basis for various comparative analyses. Here, we describe the use of MBGD and briefly explain how to utilize the orthology information during comparative genome analysis in combination with the stand-alone comparative genomics software RECOG, focusing on the application to comparison of closely related microbial genomes.

The microbial genome database for comparative analysis (MBGD) (available at http://mbgd.genome.ad.jp/) is a comprehensive ortholog database for flexible comparative analysis of microbial genomes, where the users are allowed to create an ortholog table among any specified set of organisms. Because of the rapid increase in microbial genome data owing to the next-generation sequencing technology, it becomes increasingly challenging to maintain high-quality orthology relationships while allowing the users to incorporate the latest genomic data available into an analysis. Because many of the recently accumulating genomic data are draft genome sequences for which some complete genome sequences of the same or closely related species are available, MBGD now stores draft genome data and allows the users to incorporate them into a user-specific ortholog database using the MyMBGD functionality. In this function, draft genome data are incorporated into an existing ortholog table created only from the complete genome data in an incremental manner to prevent low-quality draft data from affecting clustering results. In addition, to provide high-quality orthology relationships, the standard ortholog table containing all the representative genomes, which is first created by the rapid classification program DomClust, is now refined using DomRefine, a recently developed program for improving domain-level clustering using multiple sequence alignment information.

Background: Although the infection rate of Helicobacter suis is significantly lower than that of Helicobacter pylori, the H. suis infection is associated with a high rate of gastric mucosa-associated lymphoid tissue (MALT) lymphoma. In addition, in vitro cultivation of H. suis remains difficult, and some H. suis-infected patients show negative results on the urea breath test (UBT). Materials and Methods: Female C57BL/6J mice were orally inoculated with mouse gastric mucosal homogenates containing H. suis strains TKY or SNTW101 isolated from a cynomolgus monkey or a patient suffering from nodular gastritis, respectively. The high-purity chromosomal DNA samples of H. suis strains TKY and SNTW101 were prepared from the infected mouse gastric mucosa. The SOLiD sequencing of two H. suis genomes enabled comparative genomics of 20 Helicobacter and 11 Campylobacter strains for the identification of the H. suis-specific nucleotide sequences. Results: Oral inoculation with mouse gastric mucosal homogenates containing H. suis strains TKY and SNTW101 induced gastric MALT lymphoma and the formation of gastric lymphoid follicles, respectively, in C57BL/6J mice. Two conserved nucleotide sequences among six H. suis strains were identified and were used to design diagnostic PCR primers for the detection of H. suis. Conclusions: There was a strong association between the H. suis infection and gastric diseases in the C57BL/6 mouse model. PCR diagnosis using an H. suis-specific primer pair is a valuable method for detecting H. suis in gastric biopsy specimens.

Publishing databases in the Resource Description Framework (RDF) model is becoming widely accepted to maximize the syntactic and semantic interoperability of open data in life sciences. Here we report advancements made in the 6th and 7th annual BioHackathons which were held in Tokyo and Miyagi respectively. This review consists of two major sections covering: 1) improvement and utilization of RDF data in various domains of the life sciences and 2) meta-data about these RDF data, the resources that store them, and the service quality of SPARQL Protocol and RDF Query Language (SPARQL) endpoints. The first section describes how we developed RDF data, ontologies and tools in genomics, proteomics, metabolomics, glycomics and by literature text mining. The second section describes how we defined descriptions of datasets, the provenance of data, and quality assessment of services and service discovery. By enhancing the harmonization of these two layers of machine-readable data and knowledge, we improve the way community wide resources are developed and published. Moreover, we outline best practices for the future, and prepare ourselves for an exciting and unanticipatable variety of real world applications in coming years.

Daphnia magna has been used extensively to evaluate organism- and population-level responses to pollutants in acute toxicity and reproductive toxicity tests. We have previously reported that exposure to juvenile hormone (JH) agonists results in a reduction of reproductive function and production of male offspring in a cyclic parthenogenesis, D. magna. Recent advances in molecular techniques have provided tools to understand better the responses to pollutants in aquatic organisms, including D. magna. DNA microarray was used to evaluate gene expression profiles of neonatal daphnids exposed to JH agonists: methoprene (125, 250 and 500 ppb), fenoxycarb (0.5, 1 and 2 ppb) and epofenonane (50, 100 and 200 ppb). Exposure to these JH analogs resulted in chemical-specific patterns of gene expression. The heat map analyses based on hierarchical clustering revealed a similar pattern between treatments with a high dose of methoprene and with epofenonane. In contrast, treatment with low to middle doses of methoprene resulted in similar profiles to fenoxycarb treatments. Hemoglobin and JH epoxide hydrolase genes were clustered as JH-responsive genes. These data suggest that fenoxycarb has high activity as a JH agonist, methoprene shows high toxicity and epofenonane works through a different mechanism compared with other JH analogs, agreeing with data of previously reported toxicity tests. In conclusion, D. magna DNA microarray is useful for the classification of JH analogs and identification of JH-responsive genes. Copyright © 2013 John Wiley & Sons, Ltd.

Identifying population structure forms an important basis for genetic and evolutionary studies. Most current methods to identify population structure have limitations in analyzing haplotypes and recombination across the genome. Recently, a method of chromosome painting in silico has been developed to overcome these shortcomings and has been applied to multiple human genome sequences. This method detects the genome-wide transfer of DNA sequence chunks through homologous recombination. Here, we apply it to the frequently recombining bacterial species Helicobacter pylori that has infected Homo sapiens since their birth in Africa and shows wide phylogeographic divergence. Multiple complete genome sequences were analyzed including sequences from Okinawa, Japan, that we recently sequenced. The newer method revealed a finer population structure than revealed by a previous method that examines only MLST housekeeping genes or a phylogenetic network analysis of the core genome. Novel subgroups were found in Europe, Amerind, and East Asia groups. Examination of genetic flux showed some singleton strains to be hybrids of subgroups and revealed evident signs of population admixture in Africa, Europe, and parts of Asia. We expect this approach to further our understanding of intraspecific bacterial evolution by revealing population structure at a finer scale.

Ipomoea is the largest genus in the family Convolvulaceae. Ipomoea nil (Japanese morning glory) has been utilized as a model plant to study the genetic basis of floricultural traits, with over 1,500 mutant lines. In the present study, we have utilized second- and third-generation-sequencing platforms, and have reported a draft genome of I. nil with a scaffold N50 of 2.88 Mb (contig N50 of 1.87 Mb), covering 98% of the 750 Mb genome. Scaffolds covering 91.42% of the assembly are anchored to 15 pseudo-chromosomes. The draft genome has enabled the identification and cataloguing of the Tpn1 family transposons, known as the major mutagen of I. nil, and analysing the dwarf gene, CONTRACTED, located on the genetic map published in 1956. Comparative genomics has suggested that a whole genome duplication in Convolvulaceae, distinct from the recent Solanaceae event, has occurred after the divergence of the two sister families., Japanese morning glory (Ipomoea nil) has diverse flowering traits. Here, the authors describe the reference genome sequence of I. nil, annotations of genes and transposons, and compare evolution of the I. nil genome to other Convolvulaceae and Solanales genomes.

We present here the draft whole-genome shotgun sequence of an uncultivated strain SNTW101 of Helicobacter suis , which has been maintained in the stomachs of mice. This strain was originally isolated from gastric biopsy specimens of a urea breath test-negative Japanese patient suffering from nodular gastritis.

Aurantimicrobium minutum type strain KNC T is a planktonic ultramicrobacterium isolated from river water in western Japan. Strain KNC T has an extremely small, streamlined genome of 1,622,386 bp comprising 1,575 protein-coding sequences. The genome annotation suggests that strain KNC T has an actinorhodopsin-based photometabolism.

Genomes of a given bacterial species can show great variation in gene content and thus systematic analysis of the entire gene repertoire, termed the pan-genome, is important for understanding bacterial intra-species diversity, population genetics, and evolution. Here, we analyzed the pan-genome from 30 completely sequenced strains of the human gastric pathogen Helicobacter pylori belonging to various phylogeographic groups, focusing on 991 accessory (not fully conserved) orthologous groups (OGs). We developed a method to evaluate the mobility of genes within a genome, using the gene order in the syntenically conserved regions as a reference, and classified the 991 accessory OGs into five classes: Core, Stable, Intermediate, Mobile, and Unique. Phylogenetic networks based on the gene content of Core and Stable classes are highly congruent with that created from the concatenated alignment of fully conserved core genes, in contrast to those of Intermediate and Mobile classes, which show quite different topologies. By clustering the accessory OGs on the basis of phylogenetic pattern similarity and chromosomal proximity, we identified 60 co-occurring gene clusters (CGCs). In addition to known genomic islands, including cag pathogenicity island, bacteriophages, and integrating conjugative elements, we identified some novel ones. One island encodes TerY-phosphorylation triad, which includes the eukaryote-type protein kinase/phosphatase gene pair, and components of type VII secretion system. Another one contains a reverse-transcriptase homolog, which may be involved in the defense against phage infection through altruistic suicide. Many of the CGCs contained restriction-modification (RM) genes. Different RM systems sometimes occupied the same (orthologous) locus in the strains. We anticipate that our method will facilitate pan-genome studies in general and help identify novel genomic islands in various bacterial species.

The microbial genome database for comparative analysis (MBGD, available at http://mbgd.genome.ad.jp/) is a platform for microbial genome comparison based on orthology analysis. As its unique feature, MBGD allows users to conduct orthology analysis among any specified set of organisms; this flexibility allows MBGD to adapt to a variety of microbial genomic study. Reflecting the huge diversity of microbial world, the number of microbial genome projects now becomes several thousands. To efficiently explore the diversity of the entire microbial genomic data, MBGD now provides summary pages for pre-calculated ortholog tables among various taxonomic groups. For some closely related taxa, MBGD also provides the conserved synteny information (core genome alignment) pre-calculated using the CoreAligner program. In addition, efficient incremental updating procedure can create extended ortholog table by adding additional genomes to the default ortholog table generated from the representative set of genomes. Combining with the functionalities of the dynamic orthology calculation of any specified set of organisms, MBGD is an efficient and flexible tool for exploring the microbial genome diversity.

"Helicobacter heilmannii" ("H. heilmannii"), which belongs to the genus Helicobacter, is a group of bacterial species that display a long spiral-shaped morphology. Recent studies have demonstrated that "H. heilmannii" type 1 is actually H. suis, which mainly colonizes the stomachs of various animals and humans. However, the influence of H. suis on gastric diseases remains to be fully elucidated. In this report, we revealed the relationship between natural H. suis infection and follicular gastritis in the pig stomachs. From sequence analysis of the 16S rRNA, urease A, and urease B genes, the presence of H. suis was confirmed in pig gastric lymphoid follicles, and this bacterium was named H. suis KB1. In addition, H. suis KB1 was inoculated into C57BL/6J mice, and the following mouse model of the pathogenesis of follicular gastritis by H. suis infection was established: H. suis KB1 colonizes the mouse stomach, and moreover, induces the development of lymphoid follicles and acquired immune responses characterized by the activation of B cells and CD4 positive cells. These results may lead to better understanding of the relationship between H. suis and gastric diseases, especially follicular gastritis; and furthermore, our findings emphasize the zoonotic aspects of animal-human infection by H. suis.

The microbial genome database (MBGD) for comparative analysis is a platform for microbial comparative genomics based on automated ortholog group identification. A prominent feature of MBGD is that it allows users to create ortholog groups using a specified subgroup of organisms. The database is constantly updated and now contains almost 1000 genomes. To utilize the MBGD database as a comprehensive resource for investigating microbial genome diversity, we have developed the following advanced functionalities: (i) enhanced assignment of functional annotation, including external database links to each orthologous group, (ii) interface for choosing a set of genomes to compare based on phenotypic properties, (iii) the addition of more eukaryotic microbial genomes (fungi and protists) and some higher eukaryotes as references and (iv) enhancement of the MyMBGD mode, which allows users to add their own genomes to MBGD and now accepts raw genomic sequences without any annotation (in such a case, it runs a gene-finding procedure before identifying the orthologs). Some analysis functions, such as the function to find orthologs with similar phylogenetic patterns, have also been improved. MBGD is accessible at http://mbgd.genome.ad.jp/.

We isolated the methicillin-resistantMacrococcus caseolyticusstrain JCSC5402 from animal meat in a supermarket and determined its whole-genome nucleotide sequence. This is the first report on the genome analysis of a macrococcal species that is evolutionarily closely related to the human pathogensStaphylococcus aureusandBacillus anthracis. The essential biological pathways ofM. caseolyticusare similar to those of staphylococci. However, the species has a small chromosome (2.1 MB) and lacks many sugar and amino acid metabolism pathways and a plethora of virulence genes that are present inS. aureus. On the other hand,M. caseolyticuspossesses a series of oxidative phosphorylation machineries that are closely related to those in the familyBacillaceae. We also discovered a probable primordial form of aMacrococcusmethicillin resistance gene complex,mecIRAm, on one of the eight plasmids harbored by theM. caseolyticusstrain. This is the first finding of a plasmid-encoding methicillin resistance gene.Macrococcusis considered to reflect the genome of ancestral bacteria before the speciation of staphylococcal species and may be closely associated with the origin of the methicillin resistance gene complex of the notorious human pathogen methicillin-resistantS. aureus.

Scrub typhus ('Tsutsugamushi' disease in Japanese) is a mite-borne infectious disease. The causative agent is Orientia tsutsugamushi, an obligate intracellular bacterium belonging to the family Rickettsiaceae of the subdivision alpha-Proteobacteria. In this study, we determined the complete genome sequence of O. tsutsugamushi strain Ikeda, which comprises a single chromosome of 2 008 987 bp and contains 1967 protein coding sequences (CDSs). The chromosome is much larger than those of other members of Rickettsiaceae, and 46.7% of the sequence was occupied by repetitive sequences derived from an integrative and conjugative element, 10 types of transposable elements, and seven types of short repeats of unknown origins. The massive amplification and degradation of these elements have generated a huge number of repeated genes (1196 CDSs, categorized into 85 families), many of which are pseudogenes (766 CDSs), and also induced intensive genome shuffling. By comparing the gene content with those of other family members of Rickettsiacea, we identified the core gene set of the family Rickettsiaceae and found that, while much more extensive gene loss has taken place among the housekeeping genes of Orientia than those of Rickettsia, O. tsutsugamushi has acquired a large number of foreign genes. The O. tsutsugamushi genome sequence is thus a prominent example of the high plasticity of bacterial genomes, and provides the genetic basis for a better understanding of the biology of O. tsutsugamushi and the pathogenesis of 'Tsutsugamushi' disease.

Comparison of closely related genome sequences can provide a clue as to how macroscopic genome polymorphisms were formed through various events of recombination. However, this approach has been limited to relatively simple polymorphisms such as insertion, deletion and inversion. In the present study, we tried to extend this approach to more complex genome polymorphisms that were observed when four genome sequences of bacterial genus Neisseria were compared. The first polymorphism was an apparent translocation (ab-cd to cd-ba; a region 'ab' was translocated). The second one was a re-ordering of adjacent regions (ab-cd-ef-gh to ef-cd-ab-gh; ab, cd and ef were in reverse order). The third one was a translocation of two adjacent regions with permutation of their order (ab-cd to cd-ab elsewhere in the genome). The fourth one was a genome-wide inversion associated with a genome-specific insertion into the joints (-ab-cd- to -y-ba-x-cd-). We were able to explain their formation by only a few steps of plausible events of recombination that involved linked IS copies and prophages. Our approach would help to reconstruct a history of apparently complex genome polymorphisms in any forms of organisms and to understand genome rearrangements in the natural environments in non-model organisms.

The microbial genome database for comparative analysis (MBGD) is a comprehensive platform for microbial comparative genomics. The central function of MBGD is to create orthologous groups among multiple genomes from precomputed all-against-all similarity relationships using the DomClust algorithm. The database now contains >300 published genomes and the number continues to grow. For researchers who are interested in ongoing genome projects, we have now started a new service called 'My MBGD,' which allows users to add their own genome sequences to MBGD for the purpose of identifying orthologs among both the new and the existing genomes. Furthermore, in order to make available the rapidly accumulating information on closely related genome sequences, we enhanced the interface for pairwise genome comparisons using the CGAT interface, which allows users to see nucleotide sequence alignments of non-coding as well as coding regions. MBGD is available at http://mbgd.genome.ad.jp/.

Analysis of evolution of paralogous genes in a genome is central to our understanding of genome evolution. Comparison of closely related bacterial genomes, which has provided clues as to how genome sequences evolve under natural conditions, would help in such an analysis. With species Staphylococcus aureus, whole-genome sequences have been decoded for seven strains. We compared their DNA sequences to detect large genome polymorphisms and to deduce mechanisms of genome rearrangements that have formed each of them. We first compared strains N315 and Mu50, which make one of the most closely related strain pairs, at the single-nucleotide resolution to catalogue all the middle-sized (more than 10 bp) to large genome polymorphisms such as indels and substitutions. These polymorphisms include two paralogous gene sets, one in a tandem paralogue gene cluster for toxins in a genomic island and the other in a ribosomal RNA operon. We also focused on two other tandem paralogue gene clusters and type I restriction-modification (RM) genes on the genomic islands. Then we reconstructed rearrangement events responsible for these polymorphisms, in the paralogous genes and the others, with reference to the other five genomes. For the tandem paralogue gene clusters, we were able to infer sequences for homologous recombination generating the change in the repeat number. These sequences were conserved among the repeated paralogous units likely because of their functional importance. The sequence specificity (S) subunit of type I RM systems showed recombination, likely at the homology of a conserved region, between the two variable regions for sequence specificity. We also noticed novel alleles in the ribosomal RNA operons and suggested a role for illegitimate recombination in their formation. These results revealed importance of recombination involving long conserved sequence in the evolution of paralogous genes in the genome.

Ortholog identification is a crucial first step in comparative genomics. Here, we present a rapid method of ortholog grouping which is effective enough to allow the comparison of many genomes simultaneously. The method takes as input all-against-all similarity data and classifies genes based on the traditional hierarchical clustering algorithm UPGMA. In the course of clustering, the method detects domain fusion or fission events, and splits clusters into domains if required. The subsequent procedure splits the resulting trees such that intra-species paralogous genes are divided into different groups so as to create plausible orthologous groups. As a result, the procedure can split genes into the domains minimally required for ortholog grouping. The procedure, named DomClust, was tested using the COG database as a reference. When comparing several clustering algorithms combined with the conventional bidirectional best-hit (BBH) criterion, we found that our method generally showed better agreement with the COG classification. By comparing the clustering results generated from datasets of different releases, we also found that our method showed relatively good stability in comparison to the BBH-based methods.

We have identified filamentous prophages, Nf (Neisserial filamentous phages), during an in silico genome comparison in Neisseria. Comparison of three genomes of Neisseria meningitidis and one of Neisseria gonorrhoeae revealed four subtypes of Nf. Eleven intact copies are located at different loci in the four genomes. Each intact copy of Nf is flanked by duplication of 5'-CT and, at its right end, carries a transposase homologue (pivNM/irg) of RNaseH/Retroviral integrase superfamily. The phylogeny of these putative transposases and that of phage-related proteins on Nfs are congruent. Following circularization of Nfs, a promoter-like sequence forms. The sequence at the junction of these predicted circular forms (5'-atCTtatat) was found in a related plasmid (pMU1) at a corresponding locus. Several structural variants of Nfs--partially inverted, internally deleted and truncated--were also identified. The partial inversion seems to be a product of site-specific recombination between two 5'-CTtat sequences that are in inverse orientation, one at its end and the other upstream of pivNM/irg. Formation of internally deleted variants probably proceeded through replicative transposition that also involved two 5'-CTtat sequences. We concluded that the PivNM/Irg transposase on Nfs integrated their circular forms into the chromosomal 5'-CT-containing sequences and probably mediated the above rearrangements.

MBGD is a workbench system for comparative analysis of completely sequenced microbial genomes. The central function of MBGD is to create an orthologous gene classification table using precomputed all-against-all similarity relationships among genes in multiple genomes. In MBGD, an automated classification algorithm has been implemented so that users can create their own classification table by specifying a set of organisms and parameters. This feature is especially useful when the user's interest is focused on some taxonomically related organisms. The created classification table is stored into the database and can be explored combining with the data of individual genomes as well as similarity relationships among genomes. Using these data, users can carry out comparative analyses from various points of view, such as phylogenetic pattern analysis, gene order comparison and detailed gene structure comparison. MBGD is accessible at http://mbgd.genome.ad.jp/.

Lysobacter enzymogenes M497-1 is a producer of commercialized achromopeptidase and is expected to harbour genes encoding various other antimicrobial enzymes. Here, we present the complete sequence of the genome of M497-1 and the expression profiles of the genes for various antimicrobial enzymes. Of the 117 peptidase-encoding genes found in the 6.1-Mb genome of M497-1, 15 genes (aside from the gene encoding the achromopeptidase) were expressed at a level higher than that of the average ribosomal protein genes in the 24-h culture. Thus, the strain was found more valuable than hitherto considered. In addition, M497-1 harbours 98 genes involved in the biosynthesis of various natural products, 16 of which are M497-1-specific across 4 Lysobacter species. A gene cluster starting at LEN_2603 through LEN_2673 among the 98 genes closely resembled the lysobactin biosynthesis gene cluster of Lysobacter sp. ATCC 53042. It is likely that M497-1 may produce lysobactin or related antibacterial compounds. Furthermore, comparative genomic analysis of M497-1 and four other Lysobacter species revealed that their core genome structure comprises 3,737 orthologous groups. Our findings are expected to advance further biotechnological application of Lysobacter spp. as a promising source of natural bioactive compounds.