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.

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.

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.

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.

Lacrimispora brassicae sp. nov. isolated from fermented cabbage, and proposal of Clostridium indicum Gundawar et al. 2019 and Clostridium methoxybenzovorans Mechichi et al. 1999 as heterotypic synonyms of Lacrimispora amygdalina (Parshina et al. 2003) Haas and Blanchard 2020 and Lacrimispora indolis (McClung and McCoy 1957) Haas and Blanchard 2020, respectively.

A Gram-stain-negative, endospore-forming, rod-shaped, indole-producing bacterial strain, designated YZC6T, was isolated from fermented cabbage. Strain YZC6T grew at 10-37  °C, pH 5.5-8.5, and with up to 2  % (w/v) NaCl. The major cellular fatty acids were C16 : 0 and C18 : 1 cis 11 dimethyl acetal. Phylogenetic analysis of the 16S rRNA gene revealed that strain YZC6T belonged to the genus Lacrimispora and was closely related to Lacrimispora aerotolerans DSM 5434T (98.3  % sequence similarity), Lacrimispora saccharolytica WM1T (98.1  %), and Lacrimispora algidixylanolytica SPL73T (98.1  %). The average nucleotide identity based on blast (below 87.8  %) and digital DNA-DNA hybridization (below 36.1 %) values between the novel isolate and its corresponding relatives showed that strain YZC6T could be readily distinguished from its closely related species. Based on genotypic, phenotypic, and chemotaxonomic data, a novel Lacrimispora species, Lacrimispora brassicae sp. nov., was proposed, with YZC6T as the type strain (=MAFF 212518T=JCM 32810T=DSM 112100T). This study also proposed Clostridium indicum Gundawar et al. 2019 as a later heterotypic synonym of Lacrimispora amygdalina (Parshina et al. 2003) Haas and Blanchard 2020 and Clostridium methoxybenzovorans Mechichi et al. 1999 as a later heterotypic synonym of Lacrimispora indolis (McClung and McCpy 1957) Haas and Blanchard 2020.

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.

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-1T and B3 were isolated from soil samples of an organic paddy in Japan. Strain D1-1T 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-1T belonged to the genus Clostridium and was closely related to Clostridium zeae CSC2T (99.7 % sequence similarity), Clostridium fungisolvens TW1T (99.7 %) and Clostridium manihotivorum CT4T (99.3 %). Strains D1-1T 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-1T and B3 could be readily distinguished from their closely related species. A novel Clostridium species, Clostridium folliculivorans sp. nov., with type strain D1-1T (=MAFF 212477T=DSM 113523T), is proposed based on genotypic and phenotypic data.

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.

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.

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.

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.

Clostridium zeae sp. nov., isolated from corn silage.

A rod-shaped, Gram-stain-negative, strictly anaerobic, catalase-negative and endospore-forming bacterial strain CSC2T was isolated from corn silage preserved in Tochigi, Japan. The strain CSC2T 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 C14 : 0, C16 : 0 and C16 : 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 CSC2T was a member of the genus Clostridium and closely related to Clostridium polyendosporum DSM 57272T (95.6 % gene sequence similarity) and Clostridium fallax ATCC 19400T (95.3 %). The genomic DNA G+C content of strain CSC2T was 31.1 mol% (whole genome analysis). The average nucleotide identity based on blast and digital DNA-DNA hybridization values between strain CSC2T 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 CSC2T as representing Clostridium zeae sp. nov. The type strain is CSC2T (=MAFF212476T=JCM 33766T=DSM 111242T).

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 B40T from hardwood bark used as a compost ingredient in Japan. The strain, characterized as B40T, 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 C16:0, C18:1 ω9c and summed feature 8. The genomic DNA G+C content of strain B40T is 42.2 mol%. Results of 16S rRNA gene sequence-based phylogenetic analysis indicated that strain B40T belongs to the genus Lactobacillus and the closest neighbours of strain B40T are Lactobacillus gigeriorum 202T (95.7 %), Lactobacillus pasteurii CRBIP 24.76T (95.6 %), Lactobacillus psittaci DSM 15354T (95.4 %), Lactobacillus fornicalis TV1018T (95.4 %) and Lactobacillus jensenii ATCC 25258T (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 B40T 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 B40T represents a novel species of the genus Lactobacillus, for which the name Lactobacillus corticis sp. nov. is proposed. The strain type is B40T (=JCM 32597T=DSM 107967T).

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 YK48GT 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 YK48GT was 42.5 mol%. Based on its 16S rRNA gene sequence, strain YK48GT represented a member of the genus Lentilactobacillus and showed the highest pairwise similarity to Lentilactobacillus rapi DSM 19907T (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 YK48GT and the type strains of phylogenetically related species were 72.2-76.6% and 19.0-21.2 %, respectively, indicating that strain YK48GT represents a novel species within the genus Lentilactobacillus. Phenotypic data further confirmed the differentiation of strain YK48GT from other members of the genus Lentilactobacillus. According to the results of the polyphasic characterization presented in this study, strain YK48GT represents a novel species of the genus Lentilactobacillus, for which the name Lentilactobacillus fungorum sp. nov. is proposed. The type strain is YK48GT (=JCM 32598T=DSM 107968T).

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 700855T) is a member of the species Eubacterium callanderi. Hence, C. methoxybenzovorans ATCC 700855T cannot be used as a reference for taxonomic study. The type strain deposited in the German Collection of Microorganism and Cell Cultures GmbH (DSM 12182T) is no longer listed in its online catalogue. Also, both the 16S rRNA gene and the whole-genome sequences of the original strain SR3T 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 SR3T was considered to be a member of L. indolis.

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

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 SG162T 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 20057T (=JCM 1115T) 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 20057T. Strains SG162T and NK01 could utilize more substrates as sole carbon sources than DSM 20057T, 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 SG162T 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 C16 : 0, C19 : 0 cyclo ω8c, and summed feature 8. On the basis of the polyphasic characterization findings, strains SG162T 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 SG162T (=JCM 32599T=DSM 107969T), and strains CD034 and S42 are also transferred to L. buchneri subsp. silagei.

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/. Contact: yn@nig.ac.jp. Supplementary information: Supplementary data are available at Bioinformatics online.