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.

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.

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.

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.

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.

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.

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

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

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.

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.

Oligosaccharide Metabolism and Lipoteichoic Acid Production in Lactobacillus gasseri and Lactobacillus paragasseri.

Lactobacillus gasseri and Lactobacillus paragasseri are human commensal lactobacilli that are candidates for probiotic application. Knowledge of their oligosaccharide metabolic properties is valuable for synbiotic application. The present study characterized oligosaccharide metabolic systems and their impact on lipoteichoic acid (LTA) production in the two organisms, i.e., L. gasseri JCM 1131T and L. paragasseri JCM 11657. The two strains grew well in medium with glucose but poorly in medium with raffinose, and growth rates in medium with kestose differed between the strains. Oligosaccharide metabolism markedly influenced their LTA production, and apparent molecular size of LTA in electrophoresis recovered from cells cultured with glucose and kestose differed from that from cells cultured with raffinose in the strains. On the other hand, more than 15-fold more LTA was observed in the L. gasseri cells cultured with raffinose when compared with glucose or kestose after incubation for 15 h. Transcriptome analysis identified glycoside hydrolase family 32 enzyme as a potential kestose hydrolysis enzyme in the two strains. Transcriptomic levels of multiple genes in the dlt operon, involved in D-alanine substitution of LTA, were lower in cells cultured with raffinose than in those cultured with kestose or glucose. This suggested that the different sizes of LTA observed among the carbohydrates tested were partly due to different levels of alanylation of LTA. The present study indicates that available oligosaccharide has the impact on the LTA production of the industrially important lactobacilli, which might influence their probiotic properties.

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