Comprehensive Genome Analysis of Neisseria meningitidis from South America Reveals a Distinctive Pathogenicity-Related Prophage Repertoire.

Neisseria meningitidis, a bacterium that colonizes in the human nasopharynx, occasionally causes invasive meningococcal disease leading to meningitis or septicemia. Different serogroups and lineages (clonal complexes) are related to the occurrence and epidemiology of N. meningitidis. Despite vaccines for most serogroups, N. meningitidis lineages causing unusual clinical manifestations and a higher fatality rate compared to other lineages have been reported in South America. The present study focused on exploring the diversity of N. meningitidis prophages from South America and their relationship with the epidemiological variables of these strains. We found a high diversity of prophages among the different clonal complexes. By comparing them with previously described N. meningitidis phages and prophages, we revealed groups of prophages sharing similar compositions, which could be useful for prophage comparison in N. meningitidis. Furthermore, we observed a high correlation between the prophage content and epidemiological features, e.g., pathogenicity or clonal complex. Additionally, a distinctive filamentous prophage named here as IMSAR-11 (Invasive Meningococci from South America Related to cc11) was identified. Interestingly, two versions of IMSAR-11, circular and chromosomally integrated, were found. Overall, this study reinforces the importance of the genomic characterization of circulating N. meningitidis lineages to generate new targets for lineage monitoring, diagnosis, or appropriateness of vaccine development. Further studies are necessary to understand the role of these prophages in the persistence, dispersal, and virulence of N. meningitidis in the world.

Continuation and replacement of Vibrio cholerae non-O1 clonal genomic groups isolated from Plecoglossus altivelis fish in freshwaters.

The dissemination and abundances of Vibrio species in aquatic environments are of interest, as some species cause emerging diseases in humans and in aquatic organisms like fish. It is suggested that Vibrio cholerae non-O1 infections of Plecoglossus altivelis ('ayu') were spread to various parts of Japan through the annual transplantation of juvenile fish. To investigate this, we used genome-aided tracing of 17 V. cholerae strains isolated from ayu between the 1970s and 1990s in different Japanese freshwater systems. The strains formed a genomic clade distinct from all known clades, which we designate as the Ayu clade. Two clonal genomic groups identified within the clade, Ayu-1 and Ayu-2, persisted for a few years (between 1977 to 1979 and 1987 to 1990, respectively), and clonal replacement of Ayu-1 by Ayu-2 took place over an 8-year period. Despite the high similarity between Ayu-1 and Ayu-2 (> 99.9% identity and > 97% fraction of genomes shared), differences in their gene repertoires were found, raising the possibility that they are phenotypically distinct. These results highlight the importance of genome-based studies for understanding the long-term dynamics of populations over the timescale of years.

Multi-gene phylogenetic analysis reveals that shochu-fermenting Saccharomyces cerevisiae strains form a distinct sub-clade of the Japanese sake cluster.

Shochu is a traditional Japanese distilled spirit. The formation of the distinguishing flavour of shochu produced in individual distilleries is attributed to putative indigenous yeast strains. In this study, we performed the first (to our knowledge) phylogenetic classification of shochu strains based on nucleotide gene sequences. We performed phylogenetic classification of 21 putative indigenous shochu yeast strains isolated from 11 distilleries. All of these strains were shown or confirmed to be Saccharomyces cerevisiae, sharing species identification with 34 known S. cerevisiae strains (including commonly used shochu, sake, ale, whisky, bakery, bioethanol and laboratory yeast strains and clinical isolate) that were tested in parallel. Our analysis used five genes that reflect genome-level phylogeny for the strain-level classification. In a first step, we demonstrated that partial regions of the ZAP1, THI7, PXL1, YRR1 and GLG1 genes were sufficient to reproduce previous sub-species classifications. In a second step, these five analysed regions from each of 25 strains (four commonly used shochu strains and the 21 putative indigenous shochu strains) were concatenated and used to generate a phylogenetic tree. Further analysis revealed that the putative indigenous shochu yeast strains form a monophyletic group that includes both the shochu yeasts and a subset of the sake group strains; this cluster is a sister group to other sake yeast strains, together comprising a sake-shochu group. Differences among shochu strains were small, suggesting that it may be possible to correlate subtle phenotypic differences among shochu flavours with specific differences in genome sequences. Copyright © 2017 John Wiley & Sons, Ltd.

Birth and death of genes linked to chromosomal inversion.

The birth and death of genes is central to adaptive evolution, yet the underlying genome dynamics remain elusive. The availability of closely related complete genome sequences helps to follow changes in gene contents and clarify their relationship to overall genome organization. Helicobacter pylori, bacteria in our stomach, are known for their extreme genome plasticity through mutation and recombination and will make a good target for such an analysis. In comparing their complete genome sequences, we found that gain and loss of genes (loci) for outer membrane proteins, which mediate host interaction, occurred at breakpoints of chromosomal inversions. Sequence comparison there revealed a unique mechanism of DNA duplication: DNA duplication associated with inversion. In this process, a DNA segment at one chromosomal locus is copied and inserted, in an inverted orientation, into a distant locus on the same chromosome, while the entire region between these two loci is also inverted. Recognition of this and three more inversion modes, which occur through reciprocal recombination between long or short sequence similarity or adjacent to a mobile element, allowed reconstruction of synteny evolution through inversion events in this species. These results will guide the interpretation of extensive DNA sequencing results for understanding long- and short-term genome evolution in various organisms and in cancer cells.

Unique genomic sequences in a novel Mycobacterium avium subsp. hominissuis lineage enable fine scale transmission route tracing during pig movement.

Mycobacterium avium subsp. hominissuis (MAH) is one of the most prevalent mycobacteria causing non-tuberculous mycobacterial disease in humans and animals. Of note, MAH is a major cause of mycobacterial granulomatous mesenteric lymphadenitis outbreaks in pig populations. To determine the precise source of infection of MAH in a pig farm and to clarify the epidemiological relationship among pig, human and environmental MAH lineages, we collected 50 MAH isolates from pigs reared in Japan and determined draft genome sequences of 30 isolates. A variable number of tandem repeat analysis revealed that most pig MAH isolates in Japan were closely related to North American, European and Russian human isolates but not to those from East Asian human and their residential environments. Historical recombination analysis revealed that most pig isolates could be classified into SC2/4 and SC3, which contain MAH isolated from pig, European human and environmental isolates. Half of the isolates in SC2/4 had many recombination events with MAH lineages isolated from humans in East Asia. To our surprise, four isolates belonged to a new lineage (SC5) in the global MAH population. Members of SC5 had few footprints of inter-lineage recombination in the genome, and carried 80 unique genes, most of which were located on lineage specific-genomic islands. Using unique genetic features, we were able to trace the putative transmission route via their host pigs. Together, we clarify the possibility of species-specificity of MAH in addition to local adaptation. Our results highlight two transmission routes of MAH, one exposure on pig farms from the environment and the other via pig movement. Moreover, our study also warns that the evolution of MAH in pigs is influenced by MAH from patients and their residential environments, even if the MAH are genetically distinct.

MBGD update 2010: toward a comprehensive resource for exploring microbial genome diversity.

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/.

A putative mobile genetic element carrying a novel type IIF restriction-modification system (PluTI).

Genome comparison and genome context analysis were used to find a putative mobile element in the genome of Photorhabdus luminescens, an entomopathogenic bacterium. The element is composed of 16-bp direct repeats in the terminal regions, which are identical to a part of insertion sequences (ISs), a DNA methyltransferase gene homolog, two genes of unknown functions and an open reading frame (ORF) (plu0599) encoding a protein with no detectable sequence similarity to any known protein. The ORF (plu0599) product showed DNA endonuclease activity, when expressed in a cell-free expression system. Subsequently, the protein, named R.PluTI, was expressed in vivo, purified and found to be a novel type IIF restriction enzyme that recognizes 5'-GGCGC/C-3' (/ indicates position of cleavage). R.PluTI cleaves a two-site supercoiled substrate at both the sites faster than a one-site supercoiled substrate. The modification enzyme homolog encoded by plu0600, named M.PluTI, was expressed in Escherichia coli and shown to protect DNA from R.PluTI cleavage in vitro, and to suppress the lethal effects of R.PluTI expression in vivo. These results suggested that they constitute a restriction-modification system, present on the putative mobile element. Our approach thus allowed detection of a previously uncharacterized family of DNA-interacting proteins.

Evolution in an oncogenic bacterial species with extreme genome plasticity: Helicobacter pylori East Asian genomes.

BACKGROUND: The genome of Helicobacter pylori, an oncogenic bacterium in the human stomach, rapidly evolves and shows wide geographical divergence. The high incidence of stomach cancer in East Asia might be related to bacterial genotype. We used newly developed comparative methods to follow the evolution of East Asian H. pylori genomes using 20 complete genome sequences from Japanese, Korean, Amerind, European, and West African strains. RESULTS: A phylogenetic tree of concatenated well-defined core genes supported divergence of the East Asian lineage (hspEAsia; Japanese and Korean) from the European lineage ancestor, and then from the Amerind lineage ancestor. Phylogenetic profiling revealed a large difference in the repertoire of outer membrane proteins (including oipA, hopMN, babABC, sabAB and vacA-2) through gene loss, gain, and mutation. All known functions associated with molybdenum, a rare element essential to nearly all organisms that catalyzes two-electron-transfer oxidation-reduction reactions, appeared to be inactivated. Two pathways linking acetyl~CoA and acetate appeared intact in some Japanese strains. Phylogenetic analysis revealed greater divergence between the East Asian (hspEAsia) and the European (hpEurope) genomes in proteins in host interaction, specifically virulence factors (tipα), outer membrane proteins, and lipopolysaccharide synthesis (human Lewis antigen mimicry) enzymes. Divergence was also seen in proteins in electron transfer and translation fidelity (miaA, tilS), a DNA recombinase/exonuclease that recognizes genome identity (addA), and DNA/RNA hybrid nucleases (rnhAB). Positively selected amino acid changes between hspEAsia and hpEurope were mapped to products of cagA, vacA, homC (outer membrane protein), sotB (sugar transport), and a translation fidelity factor (miaA). Large divergence was seen in genes related to antibiotics: frxA (metronidazole resistance), def (peptide deformylase, drug target), and ftsA (actin-like, drug target). CONCLUSIONS: These results demonstrate dramatic genome evolution within a species, especially in likely host interaction genes. The East Asian strains appear to differ greatly from the European strains in electron transfer and redox reactions. These findings also suggest a model of adaptive evolution through proteome diversification and selection through modulation of translational fidelity. The results define H. pylori East Asian lineages and provide essential information for understanding their pathogenesis and designing drugs and therapies that target them.

Capsular polysaccharide switching in Streptococcus suis modulates host cell interactions and virulence.

The capsular polysaccharide (CPS) of Streptococcus suis defines various serotypes based on its composition and structure. Though serotype switching has been suggested to occur between S. suis strains, its impact on pathogenicity and virulence remains unknown. Herein, we experimentally generated S. suis serotype-switched mutants from a serotype 2 strain that express the serotype 3, 4, 7, 8, 9, or 14 CPS. The effects of serotype switching were then investigated with regards to classical properties conferred by presence of the serotype 2 CPS, including adhesion to/invasion of epithelial cells, resistance to phagocytosis by macrophages, killing by whole blood, dendritic cell-derived pro-inflammatory mediator production and virulence using mouse and porcine infection models. Results demonstrated that these properties on host cell interactions were differentially modulated depending on the switched serotypes, although some different mutations other than loci of CPS-related genes were found in each the serotype-switched mutant. Among the serotype-switched mutants, the mutant expressing the serotype 8 CPS was hyper-virulent, whereas mutants expressing the serotype 3 or 4 CPSs had reduced virulence. By contrast, switching to serotype 7, 9, or 14 CPSs had little to no effect. These findings suggest that serotype switching can drastically alter S. suis virulence and host cell interactions.

A Standardized Procedure for Monitoring Harmful Algal Blooms in Chile by Metabarcoding Analysis.

Harmful algae blooms (HABs) monitoring has been implemented worldwide, and Chile, a country famous for its fisheries and aquaculture, has intensively used microscopic and toxin analyses for decades for this purpose. Molecular biological methods, such as high-throughput DNA sequencing and bacterial assemblage-based approaches, are just beginning to be introduced in Chilean HAB monitoring, and the procedures have not yet been standardized. Here, 16S rRNA and 18S rRNA metabarcoding analyses for monitoring Chilean HABs are introduced stepwise. According to a recent hypothesis, algal-bacterial mutualistic association plays a critical synergetic or antagonistic relationship accounting for bloom initiation, maintenance, and regression. Thus, monitoring HAB from algal-bacterial perspectives may provide a broader understanding of HAB mechanisms and the basis for early warning. Metabarcoding analysis is one of the best suited molecular-based tools for this purpose because it can detect massive algal-bacterial taxonomic information in a sample. The visual procedures of sampling to metabarcoding analysis herein provide specific instructions, aiming to reduce errors and collection of reliable data.

Genomic features of Mycobacterium avium subsp. hominissuis isolated from pigs in Japan.

Mycobacterium avium subsp. hominissuis (MAH) is one of the most important agents causing non-tuberculosis mycobacterial infection in humans and pigs. There have been advances in genome analysis of MAH from human isolates, but studies of isolates from pigs are limited despite its potential source of infection to human. Here, we obtained 30 draft genome sequences of MAH from pigs reared in Japan. The 30 draft genomes were 4,848,678-5,620,788 bp in length, comprising 4652-5388 coding genes and 46-75 (median: 47) tRNAs. All isolates had restriction modification-associated genes and 185-222 predicted virulence genes. Two isolates had tRNA arrays and one isolate had a clustered regularly interspaced short palindromic repeat (CRISPR) region. Our results will be useful for evaluation of the ecology of MAH by providing a foundation for genome-based epidemiological studies.

Single-Cell Genomics of Novel Actinobacteria With the Wood-Ljungdahl Pathway Discovered in a Serpentinizing System.

Serpentinite-hosted systems represent modern-day analogs of early Earth environments. In these systems, water-rock interactions generate highly alkaline and reducing fluids that can contain hydrogen, methane, and low-molecular-weight hydrocarbons-potent reductants capable of fueling microbial metabolism. In this study, we investigated the microbiota of Hakuba Happo hot springs (∼50°C; pH∼10.5-11), located in Nagano (Japan), which are impacted by the serpentinization process. Analysis of the 16S rRNA gene amplicon sequences revealed that the bacterial community comprises Nitrospirae (47%), "Parcubacteria" (19%), Deinococcus-Thermus (16%), and Actinobacteria (9%), among others. Notably, only 57 amplicon sequence variants (ASV) were detected, and fifteen of these accounted for 90% of the amplicons. Among the abundant ASVs, an early-branching, uncultivated actinobacterial clade identified as RBG-16-55-12 in the SILVA database was detected. Ten single-cell genomes (average pairwise nucleotide identity: 0.98-1.00; estimated completeness: 33-93%; estimated genome size: ∼2.3 Mb) that affiliated with this clade were obtained. Taxonomic classification using single copy genes indicates that the genomes belong to the actinobacterial class-level clade UBA1414 in the Genome Taxonomy Database. Based on metabolic pathway predictions, these actinobacteria are anaerobes, capable of glycolysis, dissimilatory nitrate reduction and CO2 fixation via the Wood-Ljungdahl (WL) pathway. Several other genomes within UBA1414 and two related class-level clades also encode the WL pathway, which has not yet been reported for the Actinobacteria phylum. For the Hakuba actinobacterium, the energy metabolism related to the WL pathway is likely supported by a combination of the Rnf complex, group 3b and 3d [NiFe]-hydrogenases, [FeFe]-hydrogenases, and V-type (H+/Na+ pump) ATPase. The genomes also harbor a form IV ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) complex, also known as a RubisCO-like protein, and contain signatures of interactions with viruses, including clustered regularly interspaced short palindromic repeat (CRISPR) regions and several phage integrases. This is the first report and detailed genome analysis of a bacterium within the Actinobacteria phylum capable of utilizing the WL pathway. The Hakuba actinobacterium is a member of the clade UBA1414/RBG-16-55-12, formerly within the group "OPB41." We propose to name this bacterium 'Candidatus Hakubanella thermoalkaliphilus.'

Protocols for Monitoring Harmful Algal Blooms for Sustainable Aquaculture and Coastal Fisheries in Chile.

Harmful algae blooms (HABs) cause acute effects on marine ecosystems due to their production of endogenous toxins or their enormous biomass, leading to significant impacts on local economies and public health. Although HAB monitoring has been intensively performed at spatiotemporal scales in coastal areas of the world over the last decades, procedures have not yet been standardized. HAB monitoring procedures are complicated and consist of many methodologies, including physical, chemical, and biological water sample measurements. Each monitoring program currently uses different combinations of methodologies depending on site specific purposes, and many prior programs refer to the procedures in quotations. HAB monitoring programs in Chile have adopted the traditional microscopic and toxin analyses but not molecular biology and bacterial assemblage approaches. Here we select and optimize the HAB monitoring methodologies suitable for Chilean geography, emphasizing on metabarcoding analyses accompanied by the classical tools with considerations including cost, materials and instrument availability, and easiness and efficiency of performance. We present results from a pilot study using the standardized stepwise protocols, demonstrating feasibility and plausibility for sampling and analysis for the HAB monitoring. Such specific instructions in the standardized protocol are critical obtaining quality data under various research environments involving multiple stations, different analysts, various time-points, and long HAB monitoring duration.

Draft genome sequences of Mycolicibacterium peregrinum isolated from a pig with lymphadenitis and from soil on the same Japanese pig farm.

OBJECTIVES: Mycolicibacterium peregrinum, a rapidly growing mycobacterial species, can opportunistically infect humans and other animals. Although M. peregrinum infections in animals have been reported, the infection sources are unknown, as is information on its virulence and drug resistant genes, which limits our current understanding of this bacterium. To address this knowledge gap, we obtained draft genome sequences for two M. peregrinum isolates; one from a case of pig lymphadenitis and one from the pig farm's soil. DATA DESCRIPTION: We report here the draft genome sequences of M. peregrinum isolates 131_1 and 138 (6,451,733-bp and 6,479,047-bp). They were isolated from a pig with mesenteric lymph node lymphadenitis and from soil on the Japanese farm where the pig was reared. A sequence alignment identity score of 100% was obtained by in silico DNA-DNA hybridization of the two isolates, while 98.28% (isolate 131_1) and 98.27% (isolate 138) scores were recorded for hybridization with a human isolate. Both isolates carry arr-1, AAC(2')-Ib, RbpA, mtrA and tap drug-resistance genes. Isolates 131_1 and 138 carry 234 and 236 putative virulence genes, respectively. Therefore, environment M. peregrinum is potentially drug resistant and can cause swine lymphadenitis. Our data provides valuable new information for future studies on nontuberculous mycobacteria.

Type II restriction endonuclease R.Hpy188I belongs to the GIY-YIG nuclease superfamily, but exhibits an unusual active site.

BACKGROUND: Catalytic domains of Type II restriction endonucleases (REases) belong to a few unrelated three-dimensional folds. While the PD-(D/E)XK fold is most common among these enzymes, crystal structures have been also determined for single representatives of two other folds: PLD (R.BfiI) and half-pipe (R.PabI). Bioinformatics analyses supported by mutagenesis experiments suggested that some REases belong to the HNH fold (e.g. R.KpnI), and that a small group represented by R.Eco29kI belongs to the GIY-YIG fold. However, for a large fraction of REases with known sequences, the three-dimensional fold and the architecture of the active site remain unknown, mostly due to extreme sequence divergence that hampers detection of homology to enzymes with known folds. RESULTS: R.Hpy188I is a Type II REase with unknown structure. PSI-BLAST searches of the non-redundant protein sequence database reveal only 1 homolog (R.HpyF17I, with nearly identical amino acid sequence and the same DNA sequence specificity). Standard application of state-of-the-art protein fold-recognition methods failed to predict the relationship of R.Hpy188I to proteins with known structure or to other protein families. In order to increase the amount of evolutionary information in the multiple sequence alignment, we have expanded our sequence database searches to include sequences from metagenomics projects. This search resulted in identification of 23 further members of R.Hpy188I family, both from metagenomics and the non-redundant database. Moreover, fold-recognition analysis of the extended R.Hpy188I family revealed its relationship to the GIY-YIG domain and allowed for computational modeling of the R.Hpy188I structure. Analysis of the R.Hpy188I model in the light of sequence conservation among its homologs revealed an unusual variant of the active site, in which the typical Tyr residue of the YIG half-motif had been substituted by a Lys residue. Moreover, some of its homologs have the otherwise invariant Arg residue in a non-homologous position in sequence that nonetheless allows for spatial conservation of the guanidino group potentially involved in phosphate binding. CONCLUSION: The present study eliminates a significant "white spot" on the structural map of REases. It also provides important insight into sequence-structure-function relationships in the GIY-YIG nuclease superfamily. Our results reveal that in the case of proteins with no or few detectable homologs in the standard "non-redundant" database, it is useful to expand this database by adding the metagenomic sequences, which may provide evolutionary linkage to detect more remote homologs.

Draft Genome Sequences of Mycolicibacter senuensis Isolate GF74 and Mycobacterium colombiense Isolates GF28 and GF76 from a Swine Farm in Japan.

Several nontuberculous mycobacteria (NTM) occasionally infect humans and animals. Here, we report the draft genome sequences of Mycolicibacter senuensis isolate GF74 (4,792,997 bp) and Mycobacterium colombiense isolates GF28 and GF76 (5,473,554 bp and 5,426,852 bp, respectively) isolated from a swine farm in Japan. These sequences provide further information on NTM research.

Complete genome sequence and expression profile of the commercial lytic enzyme producer Lysobacter enzymogenes M497-1.

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.

How genomes rearrange: genome comparison within bacteria Neisseria suggests roles for mobile elements in formation of complex genome polymorphisms.

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.

Low frequency of endospore-specific genes in subseafloor sedimentary metagenomes.

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 < 10%, which is consistent with the frequency of Clostridia-derived genomes (2-4%) but is lower than previous estimates based on the concentration of dipicolinic acid. Conceivable explanations for this discrepancy are as follows: the efficiency of lysis and DNA extraction of subseafloor endospore cells may have been lower than those of vegetative cells, conversion factor of dipicolinic acid content per cell may differ, and/or sporulation-related genes and other functional strategies for long-term survival in the deep subseafloor biosphere are evolutionarily distinct from known spore-forming gene repertoires.

Complete genome and gene expression analyses of Asaia bogorensis reveal unique responses to culture with mammalian cells as a potential opportunistic human pathogen.

Asaia bogorensis, a member of acetic acid bacteria (AAB), is an aerobic bacterium isolated from flowers and fruits, as well as an opportunistic pathogen that causes human peritonitis and bacteraemia. Here, we determined the complete genomic sequence of the As. bogorensis type strain NBRC 16594, and conducted comparative analyses of gene expression under different conditions of co-culture with mammalian cells and standard AAB culture. The genome of As. bogorensis contained 2,758 protein-coding genes within a circular chromosome of 3,198,265 bp. There were two complete operons encoding cytochrome bo3-type ubiquinol terminal oxidases: cyoABCD-1 and cyoABCD-2. The cyoABCD-1 operon was phylogenetically common to AAB genomes, whereas the cyoABCD-2 operon belonged to a lineage distinctive from the cyoABCD-1 operon. Interestingly, cyoABCD-1 was less expressed under co-culture conditions than under the AAB culture conditions, whereas the converse was true for cyoABCD-2. Asaia bogorensis shared pathogenesis-related genes with another pathogenic AAB, Granulibacter bethesdensis, including a gene coding pathogen-specific large bacterial adhesin and additional genes for the inhibition of oxidation and antibiotic resistance. Expression alteration of the respiratory chain and unique hypothetical genes may be key traits that enable the bacterium to survive under the co-culture conditions.