DDBJ update in 2023: the MetaboBank for metabolomics data and associated metadata

The Bioinformation and DNA Data Bank of Japan (DDBJ) Center (https://www.ddbj.nig.ac.jp) provides database archives that cover a wide range of fields in life sciences. As a founding member of the International Nucleotide Sequence Database Collaboration (INSDC), DDBJ accepts and distributes nucleotide sequence data as well as their study and sample information along with the National Center for Biotechnology Information in the United States and the European Bioinformatics Institute (EBI). Besides INSDC databases, the DDBJ Center provides databases for functional genomics (GEA: Genomic Expression Archive), metabolomics (MetaboBank) and human genetic and phenotypic data (JGA: Japanese Genotype-phenotype Archive). These database systems have been built on the National Institute of Genetics (NIG) supercomputer, which is also open for domestic life science researchers to analyze large-scale sequence data. This paper reports recent updates on the archival databases and the services of the DDBJ Center, highlighting the newly redesigned MetaboBank. MetaboBank uses BioProject and BioSample in its metadata description making it suitable for multi-omics large studies. Its collaboration with MetaboLights at EBI brings synergy in locating and reusing public data.

A practical assembly guideline for genomes with various levels of heterozygosity

Although current long-read sequencing technologies have a long-read length that facilitates assembly for genome reconstruction, they have high sequence errors. While various assemblers with different perspectives have been developed, no systematic evaluation of assemblers with long reads for diploid genomes with varying heterozygosity has been performed. Here, we evaluated a series of processes, including the estimation of genome characteristics such as genome size and heterozygosity, de novo assembly, polishing, and removal of allelic contigs, using six genomes with various heterozygosity levels. We evaluated five long-read-only assemblers (Canu, Flye, miniasm, NextDenovo and Redbean) and five hybrid assemblers that combine short and long reads (HASLR, MaSuRCA, Platanus-allee, SPAdes and WENGAN) and proposed a concrete guideline for the construction of haplotype representation according to the degree of heterozygosity, followed by polishing and purging haplotigs, using stable and high-performance assemblers: Redbean, Flye and MaSuRCA.

Large-scale epidemiological study on feline autosomal dominant polycystic kidney disease and identification of novel PKD1 gene variants

OBJECTIVES

Autosomal dominant polycystic kidney disease (ADPKD) is a common inherited disease in cats. In most cases, the responsible abnormality is a nonsense single nucleotide polymorphism in exon 29 of the PKD1 gene (chrE3:g.42858112C>A, the conventional PKD1 variant). The aim of this study was to conduct a large-scale epidemiological study of ADPKD caused by the conventional PKD1 variant in Japan and to search for novel polymorphisms by targeted resequencing of the PKD1 using a next-generation sequencer.

METHODS

A total of 1281 cats visiting the Veterinary Medical Center of the University of Tokyo were included in this study. DNA was extracted from the blood of each cat. We established a novel TaqMan real-time PCR genotyping assay for the conventional PKD1 variant, and all cases were examined for the presence of this variant. Targeted resequencing of all exons of the PKD1 was performed on the DNA of 23 cats with the conventional PKD1 variant, six cats diagnosed with cystic kidneys but without this variant, and 61 wild-type normal cats.

RESULTS

Among the 1281 cats examined in this study, 23 (1.8%) harboured the conventional PKD1 variant. The odds of having the conventional PKD1 variant were significantly higher in Persian cats, Scottish Folds and Exotic Shorthairs than in the other breeds, although the number of cases in each breed was small. Furthermore, we identified four variants unique to cats with cystic kidneys that were not found in wild-type normal cats, all of which were in exon 15. In particular, two (chrE:g.42848725delC, pGly1641fs and chrE:g.42850283C>T, pArg2162Trp) were candidate variants.

CONCLUSIONS AND RELEVANCE

This study revealed that the conventional PKD1 variant was prevalent in Scottish Fold, Persian and Exotic Shorthair breeds in Japan, and variants in exon 15 of PKD1, in addition to the conventional variant in exon 29, would be key factors in the pathogenesis of ADPKD in cats.

Real-world treatment patterns and outcomes of abemaciclib for the treatment of HR + , HER2- metastatic breast cancer patients in Japan

INTRODUCTION

This study described, in routine clinical practice in Japan, the patient characteristics, treatment patterns, and outcomes of female patients with HR + /HER2- metastatic breast cancer (MBC) who started abemaciclib treatment.

METHODS

Clinical charts were reviewed for patients starting abemaciclib in 12/2018-08/2021 with a minimum of 3 months follow-up data post-abemaciclib initiation regardless of abemaciclib discontinuation. Patient characteristics, treatment patterns, and tumor response were descriptively summarized. Kaplan-Meier curves estimated progression-free survival (PFS).

RESULTS

200 patients from 14 institutions were included. At abemaciclib initiation, median age was 59 years, and the Eastern Cooperative Oncology Group performance status score was 0/1/2 for 102/68/5 patients (58.3/38.9/2.9%), respectively. Most had an abemaciclib starting dose of 150 mg (92.5%). The percentage of patients receiving abemaciclib as 1st, 2nd, or 3rd line treatment was 31.5%, 25.8%, and 25.2%, respectively. The most frequent endocrine therapy drugs used with abemaciclib were fulvestrant (59%) and aromatase inhibitors (40%). Evaluation of tumor response was available for 171 patients, 30.4% of whom had complete/partial response. Median PFS was 13.0 months (95% CI 10.1-15.8 months).

CONCLUSIONS

In a routine clinical practice setting in Japan, patients with HR + , HER2- MBC appear to benefit from abemaciclib treatment in terms of treatment response and median PFS, with the results broadly reflecting the evidence demonstrated in clinical trials.

Clostridium folliculivorans sp. nov., isolated from soil samples of an organic paddy in Japan

Two Gram-stain-negative, terminal endospore-forming, rod-shaped and aerotolerant bacterial strains designated D1-1^T and B3 were isolated from soil samples of an organic paddy in Japan. Strain D1-1^T grew at 15-37 °C, pH 5.0-7.3, and with up to 0.5 % (w/v) NaCl. Phylogenetic analysis of the 16S rRNA gene revealed that strain D1-1^T belonged to the genus Clostridium and was closely related to Clostridium zeae CSC2^T (99.7 % sequence similarity), Clostridium fungisolvens TW1^T (99.7 %) and Clostridium manihotivorum CT4^T (99.3 %). Strains D1-1^T and B3 were whole-genome sequenced and indistinguishable, with an average nucleotide identity value of 99.7 %. The average nucleotide identity (below 91.1 %) and digital DNA-DNA hybridization (below 43.6 %) values between the two novel isolates and their corresponding relatives showed that strains D1-1^T and B3 could be readily distinguished from their closely related species. A novel Clostridium species, Clostridium folliculivorans sp. nov., with type strain D1-1^T (=MAFF 212477^T=DSM 113523^T), is proposed based on genotypic and phenotypic data.

A leaf-emanated signal orchestrates grain size and number in response to maternal resources

In plants, variations in seed size and number are outcomes of different reproductive strategies. Both traits are often environmentally influenced, suggesting that a mechanism exists to coordinate these phenotypes in response to available maternal resources. Yet, how maternal resources are sensed and influence seed size and number is largely unknown. Here, we report a mechanism that senses maternal resources and coordinates grain size and number in the wild rice Oryza rufipogon, a wild progenitor of Asian cultivated rice. We showed that FT-like 9 (FTL9) regulates both grain size and number and that maternal photosynthetic assimilates induce FTL9 expression in leaves to act as a long-range signal that increases grain number and reduces size. Our findings highlight a strategy that benefits wild plants to survive in a fluctuating environment. In this strategy, when maternal resources are sufficient, wild plants increase their offspring number while preventing an increase in offspring size by the action of FTL9, which helps expand their habitats. In addition, we found that a loss-of-function allele (ftl9) is prevalent among wild and cultivated populations, offering a new scenario in the history of rice domestication.

Genome Sequence and Analysis of Nicotiana benthamiana, the Model Plant for Interactions between Organisms

Nicotiana benthamiana is widely used as a model plant for dicotyledonous angiosperms. In fact, the strains used in research are highly susceptible to a wide range of viruses. Accordingly, these strains are subject to plant pathology and plant-microbe interactions. In terms of plant-plant interactions, N. benthamiana is one of the plants that exhibit grafting affinity with plants from different families. Thus, N. benthamiana is a good model for plant biology and has been the subject of genome sequencing analyses for many years. However, N. benthamiana has a complex allopolyploid genome, and its previous reference genome is fragmented into 141,000 scaffolds. As a result, molecular genetic analysis is difficult to perform. To improve this effort, de novo whole-genome assembly was performed in N. benthamiana with Hifi reads, and 1,668 contigs were generated with a total length of 3.1 Gb. The 21 longest scaffolds, regarded as pseudomolecules, contained a 2.8-Gb sequence, occupying 95.6% of the assembled genome. A total of 57,583 high-confidence gene sequences were predicted. Based on a comparison of the genome structures between N. benthamiana and N. tabacum, N. benthamiana was found to have more complex chromosomal rearrangements, reflecting the age of interspecific hybridization. To verify the accuracy of the annotations, the cell wall modification genes involved in grafting were analyzed, which revealed not only the previously indeterminate untranslated region, intron and open reading frame sequences but also the genomic locations of their family genes. Owing to improved genome assembly and annotation, N. benthamiana would increasingly be more widely accessible.

Best Practices for Comprehensive Annotation of Neuropeptides of Gryllus bimaculatus

Genome annotation is critically important data that can support research. Draft genome annotations cover representative genes; however, they often do not include genes that are expressed only in limited tissues and stages, or genes with low expression levels. Neuropeptides are responsible for regulation of various physiological and biological processes. A recent study disclosed the genome draft of the two-spotted cricket Gryllus bimaculatus, which was utilized to understand the intriguing physiology and biology of crickets. Thus far, only two of the nine reported neuropeptides in G. bimaculatus were annotated in the draft genome. Even though de novo assembly using transcriptomic analyses can comprehensively identify neuropeptides, this method does not follow those annotations on the genome locus. In this study, we performed the annotations based on the reference mapping, de novo transcriptome assembly, and manual curation. Consequently, we identified 41 neuropeptides out of 43 neuropeptides, which were reported in the insects. Further, 32 of the identified neuropeptides on the genomic loci in G. bimaculatus were annotated. The present annotation methods can be applicable for the neuropeptide annotation of other insects. Furthermore, the methods will help to generate useful infrastructures for studies relevant to neuropeptides.

DNA Data Bank of Japan (DDBJ) update report 2022

The Bioinformation and DNA Data Bank of Japan (DDBJ) Center (https://www.ddbj.nig.ac.jp) maintains database archives that cover a wide range of fields in life sciences. As a founding member of the International Nucleotide Sequence Database Collaboration (INSDC), our primary mission is to collect and distribute nucleotide sequence data, as well as their study and sample information, in collaboration with the National Center for Biotechnology Information in the United States and the European Bioinformatics Institute. In addition to INSDC resources, the Center operates databases for functional genomics (GEA: Genomic Expression Archive), metabolomics (MetaboBank), and human genetic and phenotypic data (JGA: Japanese Genotype-Phenotype Archive). These databases are built on the supercomputer of the National Institute of Genetics, whose remaining computational capacity is actively utilized by domestic researchers for large-scale biological data analyses. Here, we report our recent updates and the activities of our services.

MarpolBase Expression: A Web-Based, Comprehensive Platform for Visualization and Analysis of Transcriptomes in the Liverwort Marchantia polymorpha

The liverwort Marchantia polymorpha is equipped with a wide range of molecular and genetic tools and resources that have led to its wide use to explore the evo-devo aspects of land plants. Although its diverse transcriptome data are rapidly accumulating, there is no extensive yet user-friendly tool to exploit such a compilation of data and to summarize results with the latest annotations. Here, we have developed a web-based suite of tools, MarpolBase Expression (MBEX, https://marchantia.info/mbex/), where users can visualize gene expression profiles, identify differentially expressed genes, perform co-expression and functional enrichment analyses and summarize their comprehensive output in various portable formats. Using oil body biogenesis as an example, we demonstrated that the results generated by MBEX were consistent with the published experimental evidence and also revealed a novel transcriptional network in this process. MBEX should facilitate the exploration and discovery of the genetic and functional networks behind various biological processes in M. polymorpha and promote our understanding of the evolution of land plants.

A highly contiguous genome assembly of red perilla (Perilla frutescens) domesticated in Japan

Perilla frutescens (Lamiaceae) is an important herbal plant with hundreds of bioactive chemicals, among which perillaldehyde and rosmarinic acid are the two major bioactive compounds in the plant. The leaves of red perilla are used as traditional Kampo medicine or food ingredients. However, the medicinal and nutritional uses of this plant could be improved by enhancing the production of valuable metabolites through the manipulation of key enzymes or regulatory genes using genome editing technology. Here, we generated a high-quality genome assembly of red perilla domesticated in Japan. A near-complete chromosome-level assembly of P. frutescens was generated contigs with N50 of 41.5 Mb from PacBio HiFi reads. 99.2% of the assembly was anchored into 20 pseudochromosomes, among which seven pseudochromosomes consisted of one contig, while the rest consisted of less than six contigs. Gene annotation and prediction of the sequences successfully predicted 86,258 gene models, including 76,825 protein-coding genes. Further analysis showed that potential targets of genome editing for the engineering of anthocyanin pathways in P. frutescens are located on the late-stage pathways. Overall, our genome assembly could serve as a valuable reference for selecting target genes for genome editing of P. frutescens.

A novel species of lactic acid bacteria, Ligilactobacillus pabuli sp. nov., isolated from alfalfa silage

In this study, we isolated a novel strain of lactic acid bacteria, AF129T, from alfalfa silage prepared locally in Morioka, Iwate, Japan. Polyphasic taxonomy was used to characterize the bacterial strain. The bacterium was rod-shaped, Gram-stain-positive, non-spore-forming and catalase-negative. The strain grew at various temperatures (15–40°C) and pH levels (4.0–8.0). The optimum growth conditions were a temperature of 30°C and a pH of 6.0. AF129T exhibited growth at salt (NaCl) concentrations of up to 6.5 % (w/v). The G+C content of the strain’s genomic DNA was 41.5 %. The major fatty acids were C16 : 0, C18 : 1ω9c, C19 : 0cyclo ω8c and summed feature 8. 16S rRNA gene sequencing revealed that AF129T represents a member of the genus Ligilactobacillus and it has higher sequence similarities with Ligilactobacillus pobuzihii (98.4 %), Ligilactobacillus acidipiscis (97.5 %) and Ligilactobacillus salitolerans (97.4 %). The digital DNA–DNA hybridization values for AF129T and phylogenetically related species of the genus Ligilactobacillus ranged from 19.8% to 24.1%. The average nucleotide identity of the strain with its closely related taxa was lower than the threshold (95 %–96 %) used for species differentiation. In the light of the above-mentioned physiological, genotypic, chemotaxonomic and phylogenetic evidence, we confirm that AF129T represents a member of the genus Ligilactobacillus and constitutes a novel species; we propose the name Ligilactobacillus pabuli sp. nov. for this species. The type strain is AF129T =MAFF 518002T =JCM 34518T=BCRC 81335T.

The renaissance and enlightenment of Marchantia as a model system

The liverwort Marchantia polymorpha has been utilized as a model for biological studies since the 18th century. In the past few decades, there has been a Renaissance in its utilization in genomic and genetic approaches to investigating physiological, developmental, and evolutionary aspects of land plant biology. The reasons for its adoption are similar to those of other genetic models, e.g. simple cultivation, ready access via its worldwide distribution, ease of crossing, facile genetics, and more recently, efficient transformation, genome editing, and genomic resources. The haploid gametophyte dominant life cycle of M. polymorpha is conducive to forward genetic approaches. The lack of ancient whole-genome duplications within liverworts facilitates reverse genetic approaches, and possibly related to this genomic stability, liverworts possess sex chromosomes that evolved in the ancestral liverwort. As a representative of one of the three bryophyte lineages, its phylogenetic position allows comparative approaches to provide insights into ancestral land plants. Given the karyotype and genome stability within liverworts, the resources developed for M. polymorpha have facilitated the development of related species as models for biological processes lacking in M. polymorpha.

Genome evolution of a nonparasitic secondary heterotroph, the diatom Nitzschia putrida

Secondary loss of photosynthesis is observed across almost all plastid-bearing branches of the eukaryotic tree of life. However, genome-based insights into the transition from a phototroph into a secondary heterotroph have so far only been revealed for parasitic species. Free-living organisms can yield unique insights into the evolutionary consequence of the loss of photosynthesis, as the parasitic lifestyle requires specific adaptations to host environments. Here, we report on the diploid genome of the free-living diatom Nitzschia putrida (35 Mbp), a nonphotosynthetic osmotroph whose photosynthetic relatives contribute ca. 40% of net oceanic primary production. Comparative analyses with photosynthetic diatoms and heterotrophic algae with parasitic lifestyle revealed that a combination of gene loss, the accumulation of genes involved in organic carbon degradation, a unique secretome, and the rapid divergence of conserved gene families involved in cell wall and extracellular metabolism appear to have facilitated the lifestyle of a free-living secondary heterotroph.

Identification of the sex-determining factor in the liverwort Marchantia polymorpha reveals unique evolution of sex chromosomes in a haploid system

Sex determination is a central process for sexual reproduction and is often regulated by a sex determinant encoded on a sex chromosome. Rules that govern the evolution of sex chromosomes via specialization and degeneration following the evolution of a sex determinant have been well studied in diploid organisms. However, distinct predictions apply to sex chromosomes in organisms where sex is determined in the haploid phase of the life cycle: both sex chromosomes, female U and male V, are expected to maintain their gene functions, even though both are non-recombining. This is in contrast to the X-Y (or Z-W) asymmetry and Y (W) chromosome degeneration in XY (ZW) systems of diploids. Here, we provide evidence that sex chromosomes diverged early during the evolution of haploid liverworts and identify the sex determinant on the Marchantia polymorpha U chromosome. This gene, Feminizer, encodes a member of the plant-specific BASIC PENTACYSTEINE transcription factor family. It triggers female differentiation via regulation of the autosomal sex-determining locus of FEMALE GAMETOPHYTE MYB and SUPPRESSOR OF FEMINIZATION. Phylogenetic analyses of Feminizer and other sex chromosome genes indicate dimorphic sex chromosomes had already been established 430 mya in the ancestral liverwort. Feminizer also plays a role in reproductive induction that is shared with its gametolog on the V chromosome, suggesting an ancestral function, distinct from sex determination, was retained by the gametologs. This implies ancestral functions can be preserved after the acquisition of a sex determination mechanism during the evolution of a dominant haploid sex chromosome system.

Lentilactobacillus fungorum sp. nov., isolated from spent mushroom substrates

In Japan, during a screening of lactic acid bacteria in spent mushroom substrates, an unknown bacterium was isolated and could not be assigned to any known species. Strain YK48G^T is Gram-stain-positive, rod-shaped, non-motile, non-spore-forming and catalase-negative. The isolate grew in 0-4 % (w/v) NaCl, at 15-37 °C (optimum, 30 °C) and at pH 4.0-8.0 (optimum, pH 6.0). The genomic DNA G+C content of strain YK48G^T was 42.5 mol%. Based on its 16S rRNA gene sequence, strain YK48G^T represented a member of the genus Lentilactobacillus and showed the highest pairwise similarity to Lentilactobacillus rapi DSM 19907^T (97.86 %). Phylogenetic analyses based on amino acid sequences of 466 shared protein-encoding genes also revealed that the strain was phylogenetically positioned in the genus Lentilactobacillus but did not suggest an affiliation with previously described species. The average nucleotide identity and digital DNA-DNA hybridization values between strain YK48G^T and the type strains of phylogenetically related species were 72.2-76.6% and 19.0-21.2 %, respectively, indicating that strain YK48G^T represents a novel species within the genus Lentilactobacillus. Phenotypic data further confirmed the differentiation of strain YK48G^T from other members of the genus Lentilactobacillus. According to the results of the polyphasic characterization presented in this study, strain YK48G^T represents a novel species of the genus Lentilactobacillus, for which the name Lentilactobacillus fungorum sp. nov. is proposed. The type strain is YK48G^T (=JCM 32598^T=DSM 107968^T).

Clostridium zeae sp. nov., isolated from corn silage

A rod-shaped, Gram-stain-negative, strictly anaerobic, catalase-negative and endospore-forming bacterial strain CSC2^T was isolated from corn silage preserved in Tochigi, Japan. The strain CSC2^T grew at 15-40 °C, at pH 5.0-7.7 and with up to 0.5 % (w/v) NaCl. The main cellular fatty acids were C_14 : 0, C_16 : 0 and C_16 : 0 dimethyl acetal. The cellular polar lipids detected were diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, phosphatidic acid, lysophosphatidylethanolamine, phosphatidylserine, lysophosphatidylcholine and two unidentified polar lipids. Phylogenetic analysis of the 16S rRNA gene showed that strain CSC2^T was a member of the genus Clostridium and closely related to Clostridium polyendosporum DSM 57272^T (95.6 % gene sequence similarity) and Clostridium fallax ATCC 19400^T (95.3 %). The genomic DNA G+C content of strain CSC2^T was 31.1 mol% (whole genome analysis). The average nucleotide identity based on blast and digital DNA-DNA hybridization values between strain CSC2^T and the type strains of phylogenetically related species were below 71 and 24 %, respectively. On the basis of the genotypic, phenotypic and chemotaxonomic characteristics, it is proposed to designate strain CSC2^T as representing Clostridium zeae sp. nov. The type strain is CSC2^T (=MAFF212476^T=JCM 33766^T=DSM 111242^T).

Taxonomic status of the species Clostridium methoxybenzovorans Mechichi et al. 1999

The taxonomic status of the species Clostridium methoxybenzovorans was assessed. The 16S rRNA gene sequence, whole-genome sequence and phenotypic characterizations suggested that the type strain deposited in the American Type Culture Collection (C. methoxybenzovorans ATCC 700855^T) is a member of the species Eubacterium callanderi. Hence, C. methoxybenzovorans ATCC 700855^T cannot be used as a reference for taxonomic study. The type strain deposited in the German Collection of Microorganism and Cell Cultures GmbH (DSM 12182^T) is no longer listed in its online catalogue. Also, both the 16S rRNA gene and the whole-genome sequences of the original strain SR3^T showed high sequence identity with those of Lacrimispora indolis (recently reclassified from Clostridium indolis) as the most closely related species. Analysis of the two genomes showed average nucleotide identity based on blast and digital DNA-DNA hybridization values of 98.3 and 87.9 %, respectively. Based on these results, C. methoxybenzovorans SR3^T was considered to be a member of L. indolis.

Lactobacillus corticis sp. nov., isolated from hardwood bark

During a study on the biodiversity of bacteria that inhabit woody biomass, we isolated a strain coded B40^T from hardwood bark used as a compost ingredient in Japan. The strain, characterized as B40^T, is a Gram-stain-positive, rod-shaped, non-motile, non-spore-forming and catalase-negative bacterium. This novel isolate showed growth at 30-50 °C, at pH 3.5-7.5 and in the presence of up to 4 % (w/v) NaCl. Its major fatty acids include C_16:0, C_18:1 ω9c and summed feature 8. The genomic DNA G+C content of strain B40^T is 42.2 mol%. Results of 16S rRNA gene sequence-based phylogenetic analysis indicated that strain B40^T belongs to the genus Lactobacillus and the closest neighbours of strain B40^T are Lactobacillus gigeriorum 202^T (95.7 %), Lactobacillus pasteurii CRBIP 24.76^T (95.6 %), Lactobacillus psittaci DSM 15354^T (95.4 %), Lactobacillus fornicalis TV1018^T (95.4 %) and Lactobacillus jensenii ATCC 25258^T (95.2 %). The amino acid sequence-based phylogenetic analyses of 489 shared protein-encoding genes showed that the strain forms a phylogenetically independent lineage in the genus Lactobacillus but could not be assigned to any known species. Strain B40^T has an average nucleotide identify of <70.2 % and a digital DNA-DNA hybridization value of 19.2 % compared with the strains of other closely related Lactobacillus species. Differential genomic, phenotypic and chemotaxonomic properties, in addition to phylogenetic analyses, indicated that strain B40^T represents a novel species of the genus Lactobacillus, for which the name Lactobacillus corticis sp. nov. is proposed. The strain type is B40^T (=JCM 32597^T=DSM 107967^T).

OryzaGenome2.1: Database of Diverse Genotypes in Wild Oryza Species

BACKGROUND

OryzaGenome ( http://viewer.shigen.info/oryzagenome21detail/index.xhtml ), a feature within Oryzabase ( https://shigen.nig.ac.jp/rice/oryzabase/ ), is a genomic database for wild Oryza species that provides comparative and evolutionary genomics approaches for the rice research community.

RESULTS

Here we release OryzaGenome2.1, the first major update of OryzaGenome. The main feature in this version is the inclusion of newly sequenced genotypes and their meta-information, giving a total of 217 accessions of 19 wild Oryza species (O. rufipogon, O. barthii, O. longistaminata, O. meridionalis, O. glumaepatula, O. punctata, O. minuta, O. officinalis, O. rhizomatis, O. eichingeri, O. latifolia, O. alta, O. grandiglumis, O. australiensis, O. brachyantha, O. granulata, O. meyeriana, O. ridleyi, and O. longiglumis). These 19 wild species belong to 9 genome types (AA, BB, CC, BBCC, CCDD, EE, FF, GG, and HHJJ), representing wide genomic diversity in the genus. Using the genotype information, we analyzed the genome diversity of Oryza species. Other features of OryzaGenome facilitate the use of information on single nucleotide polymorphisms (SNPs) between O. sativa and its wild progenitor O. rufipogon in rice research, including breeding as well as basic science. For example, we provide Variant Call Format (VCF) files for genome-wide SNPs of 33 O. rufipogon accessions against the O. sativa reference genome, IRGSP1.0. In addition, we provide a new SNP Effect Table function, allowing users to identify SNPs or small insertion/deletion polymorphisms in the 33 O. rufipogon accessions and to search for the effect of these polymorphisms on protein function if they reside in the coding region (e.g., are missense or nonsense mutations). Furthermore, the SNP Viewer for 446 O. rufipogon accessions was updated by implementing new tracks for possible selective sweep regions and highly mutated regions that were potentially exposed to selective pressures during the process of domestication.

CONCLUSION

OryzaGenome2.1 focuses on comparative genomic analysis of diverse wild Oryza accessions collected around the world and on the development of resources to speed up the identification of critical trait-related genes, especially from O. rufipogon. It aims to promote the use of genotype information from wild accessions in rice breeding and potential future crop improvements. Diverse genotypes will be a key resource for evolutionary studies in Oryza, including polyploid biology.

Unique niche-specific adaptation of fructophilic lactic acid bacteria and proposal of three Apilactobacillus species as novel members of the group

BACKGROUND

Fructophilic lactic acid bacteria (FLAB) found in D-fructose rich niches prefer D-fructose over D-glucose as a growth substrate. They need electron acceptors for growth on D-glucose. The organisms share carbohydrate metabolic properties. Fructobacillus spp., Apilactobacillus kunkeei, and Apilactobacillus apinorum are members of this unique group. Here we studied the fructophilic characteristics of recently described species Apilactobacillus micheneri, Apilactobacillus quenuiae, and Apilactobacillus timberlakei.

RESULTS

The three species prefer D-fructose over D-glucose and only metabolize D-glucose in the presence of electron acceptors. The genomic characteristics of the three species, i.e. small genomes and thus a low number of coding DNA sequences, few genes involved in carbohydrate transport and metabolism, and partial deletion of adhE gene, are characteristic of FLAB. The three species thus are novel members of FLAB. Reduction of genes involved in carbohydrate transport and metabolism in accordance with reduction of genome size were the common characteristics of the family Lactobacillaceae, but FLAB markedly reduced the gene numbers more than other species in the family. Pan-genome analysis of genes involved in metabolism displayed a lack of specific carbohydrate metabolic pathways in FLAB, leading to a unique cluster separation.

CONCLUSIONS

The present study expanded FLAB group. Fructose-rich environments have induced similar evolution in phylogenetically distant FLAB species. These are examples of convergent evolution of LAB.

Whole-exome sequencing of 79 xenografts as a potential approach for the identification of genetic variants associated with sensitivity to cytotoxic anticancer drugs

Chemotherapy response remains unpredictable in most patients with cancer. In this study, we performed whole-exome sequencing of 79 cancer xenografts derived from human cancer tissues to identify genetic predictors of chemosensitivity to nine cytotoxic anticancer drugs. Xenografts were harvested from 12 organs with cancer and implanted into nude mice. The mice were exposed to one of nine cytotoxic anticancer drugs (5-fluorouracil, nimustine, adriamycin, cyclophosphamide, cisplatin, mitomycin C, methotrexate, vincristine, and vinblastine) to assess the correlation between chemosensitivity response and variant allele frequency. We found 162 candidate variants that were possibly associated with chemosensitivity to one or more of the nine anticancer drugs (P < 0.01). In a subgroup analysis of breast and gastric cancer xenografts, 78 and 67 variants, respectively, were possibly associated with chemosensitivity. This approach may help to contribute to the development of personalized treatments that may allow for the prescription of optimal chemotherapy regimens among patients with cancer.

Lactobacillus buchneri subsp. silagei subsp. nov., isolated from rice grain silage

Two Gram-stain-positive, rod-shaped, non-motile, non-spore-forming, catalase-negative bacteria, designated strains SG162^T and NK01, were isolated from Japanese rice grain silage and total mixed ration silage, respectively. They were initially identified as Lactobacillus buchneri based on the 16S rRNA gene sequence similarities. However, the two strains were separated into a distinct clade from L. buchneri DSM 20057^T (=JCM 1115^T) through whole-genome sequence-based characterization, forming an infraspecific subgroup together with strains CD034 and S42, whose genomic sequences were available in the public sequence database. Strains within the subgroup shared 99.4-99.7 % average nucleotide identity (ANI) and 97.5-99.0 % digital DNA-DNA hybridization (dDDH) with each other, albeit 96.9-97.0 % ANI and 76.0-76.6 % dDDH against DSM 20057^T. Strains SG162^T and NK01 could utilize more substrates as sole carbon sources than DSM 20057^T, potentially owing to the abundance of genes involved in carbon metabolism, especially the Entner-Doudoroff pathway. The inability of γ-aminobutyric acid (GABA) production was evidenced by the lack of glutamate decarboxylase and glutamate/GABA antiporter genes in the new subgroup strains. Strain SG162^T grew at 10-45 °C (optimum, 30 °C), pH 3.5-8.0, and 0-8 % (w/v) NaCl. Its genomic DNA G+C content was 44.1 mol%. The predominant fatty acids were C_16 : 0, C_19 : 0 cyclo ω8c, and summed feature 8. On the basis of the polyphasic characterization findings, strains SG162^T and NK01 represent a novel subspecies of L. buchneri, for which the name Lactobacillus buchneri subsp. silagei subsp. nov. is proposed. The type strain is SG162^T (=JCM 32599^T=DSM 107969^T), and strains CD034 and S42 are also transferred to L. buchneri subsp. silagei.

Reclassification of Clostridium diolis Biebl and Spröer 2003 as a later heterotypic synonym of Clostridium beijerinckii Donker 1926 (Approved Lists 1980) emend. Keis et al. 2001

Clostridium diolis shares high similarity based on 16S rRNA gene sequences and fatty acid composition with Clostridium beijerinckii . In this study, the taxonomic status of C. diolis was clarified using genomic and phenotypic approaches. High similarity was detected among C. diolis DSM 15410T, C. beijerinckii DSM 791T and NCTC 13035T, showing average nucleotide identity on blast and in silico DNA–DNA hybridization values over 97 and 85 %, respectively. Results of investigations for substrate utilization and enzyme activity displayed no striking differences between C. diolis DSM 15410T and C. beijerinckii JCM 1390T. Based on the results, we propose the reclassification of Clostridium diolis as a later heterotypic synonym of Clostridium beijerinckii . The type strain is ATCC 25752T (=CIP 104308T=DSM 791T=JCM 1390T=LMG 5716T=NCTC 13035T).

Chromatin Organization in Early Land Plants Reveals an Ancestral Association between H3K27me3, Transposons, and Constitutive Heterochromatin

Genome packaging by nucleosomes is a hallmark of eukaryotes. Histones and the pathways that deposit, remove, and read histone modifications are deeply conserved. Yet, we lack information regarding chromatin landscapes in extant representatives of ancestors of the main groups of eukaryotes, and our knowledge of the evolution of chromatin-related processes is limited. We used the bryophyte Marchantia polymorpha, which diverged from vascular plants circa 400 mya, to obtain a whole chromosome genome assembly and explore the chromatin landscape and three-dimensional genome organization in an early diverging land plant lineage. Based on genomic profiles of ten chromatin marks, we conclude that the relationship between active marks and gene expression is conserved across land plants. In contrast, we observed distinctive features of transposons and other repetitive sequences in Marchantia compared with flowering plants. Silenced transposons and repeats did not accumulate around centromeres. Although a large fraction of constitutive heterochromatin was marked by H3K9 methylation as in flowering plants, a significant proportion of transposons were marked by H3K27me3, which is otherwise dedicated to the transcriptional repression of protein-coding genes in flowering plants. Chromatin compartmentalization analyses of Hi-C data revealed that repressed B compartments were densely decorated with H3K27me3 but not H3K9 or DNA methylation as reported in flowering plants. We conclude that, in early plants, H3K27me3 played an essential role in heterochromatin function, suggesting an ancestral role of this mark in transposon silencing.

Large-scale genome analysis of bovine commensal Escherichia coli reveals that bovine-adapted E. coli lineages are serving as evolutionary sources of the emergence of human intestinal pathogenic strains

How pathogens evolve their virulence to humans in nature is a scientific issue of great medical and biological importance. Shiga toxin (Stx)–producing Escherichia coli (STEC) and enteropathogenic E. coli (EPEC) are the major foodborne pathogens that can cause hemolytic uremic syndrome and infantile diarrhea, respectively. The locus of enterocyte effacement (LEE)–encoded type 3 secretion system (T3SS) is the major virulence determinant of EPEC and is also possessed by major STEC lineages. Cattle are thought to be the primary reservoir of STEC and EPEC. However, genome sequences of bovine commensal E. coli are limited, and the emerging process of STEC and EPEC is largely unknown. Here, we performed a large-scale genomic comparison of bovine commensal E. coli with human commensal and clinical strains, including EPEC and STEC, at a global level. The analyses identified two distinct lineages, in which bovine and human commensal strains are enriched, respectively, and revealed that STEC and EPEC strains have emerged in multiple sublineages of the bovine-associated lineage. In addition to the bovine-associated lineage-specific genes, including fimbriae, capsule, and nutrition utilization genes, specific virulence gene communities have been accumulated in stx- and LEE-positive strains, respectively, with notable overlaps of community members. Functional associations of these genes probably confer benefits to these E. coli strains in inhabiting and/or adapting to the bovine intestinal environment and drive their evolution to highly virulent human pathogens under the bovine-adapted genetic background. Our data highlight the importance of large-scale genome sequencing of animal strains in the studies of zoonotic pathogens.

Significance of Granulocyte Colony-Stimulating Factor-Combined High-Dose Cytarabine, Cyclophosphamide, and Total Body Irradiation in Allogeneic Hematopoietic Cell Transplantation for Myeloid Malignant Neoplasms

Allogeneic hematopoietic cell transplant (HCT) is a curative procedure for myeloid malignant neoplasms, but relapse after HCT remains critical. A conditioning regimen involving granulocyte colony-stimulating factor-combined high-dose cytarabine, cyclophosphamide, and total body irradiation (G-CSF-combined high-dose cytarabine/cyclophosphamide/total-body irradiation [HDCA/CY/TBI]) was reported to improve outcomes after cord blood transplant (CBT) for myeloid malignant neoplasms, but this regimen was not previously evaluated among patients undergoing bone marrow transplant (BMT) or peripheral blood stem cell transplant (PBSCT).

METHODS

We retrospectively analyzed 28 patients who underwent allogeneic HCT including BMT from a related (1 patient) or unrelated donor (9 patients), PBSCT from a related donor (7 patients), or single-unit CBT from an unrelated donor (11 patients) after a G-CSF-combined HDCA/CY/TBI regimen.

RESULTS

All patients achieved neutrophil and platelet engraftment, which were significantly more rapid in the BMT/PBSCT group than in the CBT group. Eighteen patients were alive at a median follow-up of 54.3 months. The 3-year relapse and nonrelapse mortality rates were 28.6% and 7.1%, respectively, which were similar between the BMT/PBSCT and CBT groups. Overall survival and disease-free survival at 5 years after HCT were 62.6% and 64.3%, respectively, which were also similar between the BMT/PBSCT and CBT groups. Only disease status at HCT had a significant impact on overall survival and disease-free survival (86.7% with standard risk vs 38.5% with high risk and 86.7% with standard risk vs 38.5% with high risk, respectively).

CONCLUSION

A G-CSF-combined HDCA/CY/TBI regimen is a promising conditioning in patients with myeloid malignant neoplasms who undergo not only CBT but also BMT or PBSCT.

Lactobacillus salitolerans sp. nov., a novel lactic acid bacterium isolated from spent mushroom substrates

A taxonomic study of a Gram-stain-positive, rod-shaped, non-motile, non-spore-forming, catalase-negative bacterium, strain YK43T, isolated from spent mushroom substrates stored in Nagano, Japan was performed. Growth was detected at 15–45 °C, pH 5.0–8.5, and 0–10 % (w/v) NaCl. The genomic DNA G+C content of strain YK43T was 43.6 mol%. The predominant fatty acids were C16 : 0, C18 : 1 ω9c and summed feature 8. Based on 16S rRNA gene sequence analysis, the type strains of Lactobacillus acidipiscis (sequence similarity, 97.6 %) and Lactobacillus pobuzihii (97.4 %) were most closely related to YK43T. The average nucleotide identities were 74.1 % between strain YK43T and L. acidipiscis DSM 15836T and 74.0 % between YK43T and L. pobuzihii E100301T. Based on a multilocus sequence analysis, comparative genomic analysis and a range of phenotypic and chemotaxonomic characteristics, strain YK43T represents a novel species of the genus Lactobacillus , for which the name Lactobacillus salitolerans sp. nov. is proposed. The type strain is YK43T (=JCM 31331T = DSM 103433T).

Isolation and Identification of Lactic Acid Bacteria from Environmental Samples

Isolation of lactic acid bacteria (LAB) is the first and crucial step to study possible roles of LAB in the environment, especially in food fermentation. This is also important to use the organisms for further application. LAB are diverse bacterial group and have diverse growth characteristics. Culture condition of LAB is thus varied, and selection of a suitable culture medium is essential for the purposes. Identification is also an important step, since certain desirable and undesirable characteristics are shared within species. Identification was classically carried out by phenotypic characteristics but is usually performed by DNA sequence-based approaches. 16S rRNA gene sequencing is generally used for identification, and sequencing of housekeeping genes is used when needed. In addition, identification based on whole-genome sequence similarities is becoming common. Here we describe isolation and identification of LAB briefly.

DFAST: a flexible prokaryotic genome annotation pipeline for faster genome publication

Summary

We developed a prokaryotic genome annotation pipeline, DFAST, that also supports genome submission to public sequence databases. DFAST was originally started as an on-line annotation server, and to date, over 7000 jobs have been processed since its first launch in 2016. Here, we present a newly implemented background annotation engine for DFAST, which is also available as a standalone command-line program. The new engine can annotate a typical-sized bacterial genome within 10 min, with rich information such as pseudogenes, translation exceptions and orthologous gene assignment between given reference genomes. In addition, the modular framework of DFAST allows users to customize the annotation workflow easily and will also facilitate extensions for new functions and incorporation of new tools in the future.

Availability and implementation

The software is implemented in Python 3 and runs in both Python 2.7 and 3.4—on Macintosh and Linux systems. It is freely available at https://github.com/nigyta/dfast_core/under the GPLv3 license with external binaries bundled in the software distribution. An on-line version is also available at https://dfast.nig.ac.jp/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Generating Publication-Ready Prokaryotic Genome Annotations with DFAST

DDBJ Fast Annotation and Submission Tool (DFAST) is a genome annotation pipeline for prokaryotes, which also assists data submission to the public sequence database. It is available both as a web service and as a stand-alone tool that runs on local machines. DFAST can annotate a typical-sized bacterial genome within 5 min. The default annotation workflow contains a gene prediction phase for protein coding sequence, rRNA, tRNA, and CRISPR, and a functional annotation phase to infer protein functions. DFAST generates result files in standard annotation formats and data files for submission to DNA Data Bank of Japan (DDBJ). In this chapter, the annotation workflow and applications of DFAST are introduced.

Fructophilic Lactic Acid Bacteria, a Unique Group of Fructose-Fermenting Microbes

Fructophilic lactic acid bacteria (FLAB) are a recently discovered group, consisting of a few Fructobacillus and Lactobacillus species. Because of their unique characteristics, including poor growth on glucose and preference of oxygen, they are regarded as "unconventional" lactic acid bacteria (LAB). Their unusual growth characteristics are due to an incomplete gene encoding a bifunctional alcohol/acetaldehyde dehydrogenase (adhE). This results in the imbalance of NAD/NADH and the requirement of additional electron acceptors to metabolize glucose. Oxygen, fructose, and pyruvate are used as electron acceptors. FLAB have significantly fewer genes for carbohydrate metabolism than other LAB, especially due to the lack of complete phosphotransferase system (PTS) transporters. They have been isolated from fructose-rich environments, including flowers, fruits, fermented fruits, and the guts of insects that feed on plants rich in fructose, and are separated into two groups on the basis of their habitats. One group is associated with flowers, grapes, wines, and insects, and the second group is associated with ripe fruits and fruit fermentations. Species associated with insects may play a role in the health of their host and are regarded as suitable vectors for paratransgenesis in honey bees. Besides their impact on insect health, FLAB may be promising candidates for the promotion of human health. Further studies are required to explore their beneficial properties in animals and humans and their applications in the food industry.

Factors predicting early postpartum glucose intolerance in Japanese women with gestational diabetes mellitus: decision-curve analysis

AIMS

To identify factors predicting early postpartum glucose intolerance in Japanese women with gestational diabetes mellitus, using decision-curve analysis.

METHODS

A retrospective cohort study was performed. The participants were 123 Japanese women with gestational diabetes who underwent 75-g oral glucose tolerance tests at 8-12 weeks after delivery. They were divided into a glucose intolerance and a normal glucose tolerance group based on postpartum oral glucose tolerance test results. Analysis of the pregnancy oral glucose tolerance test results showed predictive factors for postpartum glucose intolerance. We also evaluated the clinical usefulness of the prediction model based on decision-curve analysis.

RESULTS

Of 123 women, 78 (63.4%) had normoglycaemia and 45 (36.6%) had glucose intolerance. Multivariable logistic regression analysis showed insulinogenic index/fasting immunoreactive insulin and summation of glucose levels, assessed during pregnancy oral glucose tolerance tests (total glucose), to be independent risk factors for postpartum glucose intolerance. Evaluating the regression models, the best discrimination (area under the curve 0.725) was obtained using the basic model (i.e. age, family history of diabetes, BMI ≥25 kg/m² and use of insulin during pregnancy) plus insulinogenic index/fasting immunoreactive insulin <1.1. Decision-curve analysis showed that combining insulinogenic index/fasting immunoreactive insulin <1.1 with basic clinical information resulted in superior net benefits for prediction of postpartum glucose intolerance.

CONCLUSIONS

Insulinogenic index/fasting immunoreactive insulin calculated using oral glucose tolerance test results during pregnancy is potentially useful for predicting early postpartum glucose intolerance in Japanese women with gestational diabetes.

Comparative analysis of probiotic bacteria based on a new definition of core genome

The commensal genus Bifidobacterium has probiotic properties. We prepared a public library of the gene functions of the genus Bifidobacterium for its online annotation. Orthologous gene cluster analysis showed that the pan genomes of Bifidobacterium and Lactobacillus exhibit striking similarities when mapped to the Clusters of Orthologous Group (COG) database of proteins. When the core genes in each genus were selected based on our statistical definition of "core genome", core genes were present in at least 92% of 52 Bifidobacterium and in 97% of 178 Lactobacillus genomes. Functional comparison of the core genes of the two genera revealed a significant difference in the categories "amino acid transport and metabolism" representing their difference in niche specificity. Over-represented Bifidobacterium protein families were primarily involved in host interactions, the complex compound metabolism, and in stress responses. These findings coincide with the published information and validate our bias-resilient definition of the core genome.

Draft Sequencing of the Heterozygous Diploid Genome of Satsuma (Citrus unshiu Marc.) Using a Hybrid Assembly Approach

Satsuma (Citrus unshiu Marc.) is one of the most abundantly produced mandarin varieties of citrus, known for its seedless fruit production and as a breeding parent of citrus. De novo assembly of the heterozygous diploid genome of Satsuma ("Miyagawa Wase") was conducted by a hybrid assembly approach using short-read sequences, three mate-pair libraries, and a long-read sequence of PacBio by the PLATANUS assembler. The assembled sequence, with a total size of 359.7 Mb at the N₅₀ length of 386,404 bp, consisted of 20,876 scaffolds. Pseudomolecules of Satsuma constructed by aligning the scaffolds to three genetic maps showed genome-wide synteny to the genomes of Clementine, pummelo, and sweet orange. Gene prediction by modeling with MAKER-P proposed 29,024 genes and 37,970 mRNA; additionally, gene prediction analysis found candidates for novel genes in several biosynthesis pathways for gibberellin and violaxanthin catabolism. BUSCO scores for the assembled scaffold and predicted transcripts, and another analysis by BAC end sequence mapping indicated the assembled genome consistency was close to those of the haploid Clementine, pummel, and sweet orange genomes. The number of repeat elements and long terminal repeat retrotransposon were comparable to those of the seven citrus genomes; this suggested no significant failure in the assembly at the repeat region. A resequencing application using the assembled sequence confirmed that both kunenbo-A and Satsuma are offsprings of Kishu, and Satsuma is a back-crossed offspring of Kishu. These results illustrated the performance of the hybrid assembly approach and its ability to construct an accurate heterozygous diploid genome.

Revealing the genomic differences between two subgroups in Lactobacillus gasseri

Being an autochthonous species in humans, Lactobacillus gasseri is widely used as a probiotic for fermented products. We thoroughly compared the gene contents of 75 L. gasseri genomes and identified two intraspecific groups by the average nucleotide identity (ANI) threshold of 94%. Group I, with 48 strains, possessed 53 group-specific genes including the gassericin T cluster (9 genes) and N-acyl homoserine lactone lactonase. Group II, with 27 strains, including the type strain ATCC 33323, possessed group-specific genes with plasmid- or phage-related annotations. The genomic differences provide evidences for demarcating a new probiotic group within L. gasseri.

Visualization of consensus genome structure without using a reference genome

BACKGROUND

Standard graphical tools for whole genome comparison require a reference genome. However, any reference is also subject to annotation biases and rearrangements, and may not serve as the standard except for those of extensively studied model species. To fully exploit the rapidly accumulating sequence data from the recent sequencing technologies, genome comparison without any reference has been anticipated.

RESULTS

We introduce a circular genome visualizer to compare complete genomes of closely related species. This tool visualizes the position of orthologous gene clusters rather than actual sequences or their features, thereby achieving the comparative view without using a single reference genome. The essential information is the matrix of orthologous gene clusters whose positions (not sequences) are color-coded in circular graphics. As a demonstration, comparison of 14 Lactobacillus paracasei strains and one L. casei strain revealed not only large-scale rearrangements but also genomic islands that are strain-specific. Comparison of 73 Helicobacter pylori strains confirmed their genetic consistency and also revealed the three general patterns of large-scale genome inversions.

CONCLUSIONS

From the ample sequence information in the GenBank/ENA/DDBJ repository, we can reconstruct a genomic consensus for particular species. By visualizing multiple strains at a glance, we can identify conserved as well as strain-specific regions in multiply sequenced genomes. Positional consistency for orthologous genes provides information orthogonal to major sequence features such as the GC content or sequence similarity of marker genes. The positional comparison is therefore useful for identifying large-scale genome rearrangements or gene transfers.

Genomic characterization reconfirms the taxonomic status of Lactobacillus parakefiri

Whole-genome sequencing was performed for Lactobacillus parakefiri JCM 8573^T to confirm its hitherto controversial taxonomic position. Here, we report its first reliable reference genome. Genome-wide metrics, such as average nucleotide identity and digital DNA-DNA hybridization, and phylogenomic analysis based on multiple genes supported its taxonomic status as a distinct species in the genus Lactobacillus. The availability of a reliable genome sequence will aid future investigations on the industrial applications of L. parakefiri in functional foods such as kefir grains.

DNApod: DNA polymorphism annotation database from next-generation sequence read archives

With the rapid advances in next-generation sequencing (NGS), datasets for DNA polymorphisms among various species and strains have been produced, stored, and distributed. However, reliability varies among these datasets because the experimental and analytical conditions used differ among assays. Furthermore, such datasets have been frequently distributed from the websites of individual sequencing projects. It is desirable to integrate DNA polymorphism data into one database featuring uniform quality control that is distributed from a single platform at a single place. DNA polymorphism annotation database (DNApod; http://tga.nig.ac.jp/dnapod/) is an integrated database that stores genome-wide DNA polymorphism datasets acquired under uniform analytical conditions, and this includes uniformity in the quality of the raw data, the reference genome version, and evaluation algorithms. DNApod genotypic data are re-analyzed whole-genome shotgun datasets extracted from sequence read archives, and DNApod distributes genome-wide DNA polymorphism datasets and known-gene annotations for each DNA polymorphism. This new database was developed for storing genome-wide DNA polymorphism datasets of plants, with crops being the first priority. Here, we describe our analyzed data for 679, 404, and 66 strains of rice, maize, and sorghum, respectively. The analytical methods are available as a DNApod workflow in an NGS annotation system of the DNA Data Bank of Japan and a virtual machine image. Furthermore, DNApod provides tables of links of identifiers between DNApod genotypic data and public phenotypic data. To advance the sharing of organism knowledge, DNApod offers basic and ubiquitous functions for multiple alignment and phylogenetic tree construction by using orthologous gene information.

Genomic characterization of a fructophilic bee symbiont Lactobacillus kunkeei reveals its niche-specific adaptation

Lactobacillus kunkeei is classified as a sole obligate fructophilic lactic acid bacterium that is found in fructose-rich niches, including the guts of honeybees. The species is differentiated from other lactobacilli based on its poor growth with glucose, enhanced growth in the presence of oxygen and other electron acceptors, and production of high concentrations of acetate from the metabolism of glucose. These characteristics are similar to phylogenetically distant Fructobacillus spp. In the present study, the genomic structure of L. kunkeei was characterized by using 16 different strains, and it had significantly less genes and smaller genomes when compared with other lactobacilli. Functional gene classification revealed that L. kunkeei had lost genes specifically involved in carbohydrate transport and metabolism. The species also lacked most of the genes for respiration, although growth was enhanced in the presence of oxygen. The adhE gene of L. kunkeei, encoding a bifunctional alcohol dehydrogenase (ADH)/aldehyde dehydrogenase (ALDH) protein, lacked the part encoding the ADH domain, which is reported here for the first time in lactic acid bacteria. The deletion resulted in the lack of ADH activity, implying a requirement for electron acceptors in glucose assimilation. These results clearly indicated that L. kunkeei had undergone a specific reductive evolution in order to adapt to fructose-rich environments. The reduction characteristics were similar to those of Fructobacillus spp., but distinct from other lactobacilli with small genomes, such as Lactobacillus gasseri and Lactobacillus vaginalis. Fructose-richness thus induced an environment-specific gene reduction in phylogenetically distant microorganisms.