Draft Genome Sequence of Naganishia liquefaciens Strain N6, Isolated from the Japan Trench.

概要

The draft genome sequence of the deep-sea yeast Naganishia liquefa- ciens strain N6, isolated from the Japan Trench, is reported here. This strain was pre- viously classified into a Cryptococcus clade. Phylogenetic analysis using the presented sequence suggests that strain N6 is in the clade of the genus Naganishia.

　Naganishia liquefaciens (phylum Basidiomycota, class Tremellomycetes) strain N6 was isolated from the Japan Trench (6,500 m) (1, 2). Sediment samples were spread onto yeast extract-peptone-dextrose (YPD) plates (3) containing penicillin and strepto- mycin sulfate, and colonies were isolated. Although initially named Cryptococcus liquefa- ciens based on a comparison of the 18S rDNA with that of Cryptococcus albidus (1, 2), Cryptococcus albidus was subsequently assigned to the genus Naganishia; Cryptococcus liquefaciens was therefore renamed Naganishia liquefaciens (4, 5). Strain N6 is tolerant to heavy metals (1, 2, 6). Many deep-sea microorganisms have evolved to survive under extreme conditions, and their characterization is potentially crucial for the pro-duction of useful biomolecules.

　Strain N6 was cultivated on YPD medium (3) at 30°C, and genomic DNA was pre- pared (Dr. GenTLE kit, TaKaRa Bio). Whole-genome sequencing was performed using the Illumina MiSeq platform. Three paired-end (kit, TruSeq; 34,806,922 reads; read length, 150 or 300 bp; insert sizes, 400 to 550 bp) and three mate pair (kit, Nextera mate pair; 25,211,788 reads; read length, 250 bp; insert sizes, 4,000 to 10,000 bp) libra- ries were generated (total, 14.7 Gbp). Default parameters were used except where oth- erwise noted. The reads were trimmed using Platanus_trim v1.0.2 (http://platanus.bio.titech.ac.jp/pltanus_trim). De novo assembly was performed by Platanus v1.2.1 (7), inputting all but the 10-kbp mate pair libraries. Misassemblies were corrected based on the physical coverage of the 10-kbp mate pairs. Some gaps were filled by additional Sanger sequencing of the PCR products of gap-flanking primers (44 reads; SRA acces- sion number DRR244395; alignment tool, BLASTN). The complete mitochondrial ge- nome (GenBank accession number BLZA01000059.1) was constructed using Platanus (“assemble -n 200” command). The numbers of resulting scaffolds and contigs, the total length, the scaffold N50 value, the contig N50 value, the gap rate, and the GC con- tent were 59, 87, 19.44 Mbp, 1.03 Mbp, 0.62 Mbp, 0.10%, and 53.38%, respectively.

　To determine strain N6’s gene structure, total RNA from cells grown in YPD plus adenine (YPAD) or YPAD containing 10 mM CuSO4 was prepared (Nucleospin RNA kit, Macherey-Nagel). Four transcriptome sequencing (RNA-seq) libraries (two repli- cates for each condition) were prepared [TruSeq kit with poly(A) selection], and RNA-seq was performed using an Illumina MiSeq instrument (read length, 300 bp), resulting in 37,287,236 reads (total 8.10 Gbp). The protein-coding gene structure on the scaffolds was predicted using FunGAP v1.0.0 (8) and proteins from Cryptococcus neoformans and Cryptococcus gattii (GenBank accession numbers GCA_000091045.1 and GCA_000185945.1, respectively). A total of 6,883 genes were predicted. The av- erage lengths of transcripts and coding DNA sequences (CDS) were 1,999 bp and 1,621 bp, respectively. In addition, 38,748 introns were identified, with 93.7% of genes containing at least one intron.

　We downloaded 113 genomes in the class Tremellomycetes and the genome of Ustilago maydis JCM 2005 as an outgroup from the GenBank database. Single-copy orthologs were identified using BUSCO v4.0.6 (9) with the basidiomycota_odb10 data set (1,764 orthologs). In the strain N6 genome, 1,613 (91.3%) benchmarking universal single-copy ortholog (BUSCO) complete genes were detected. For each ortholog group, protein sequences were aligned using MAFFT v7.455 (10), and gaps were removed. From the concatenated alignment, a maximum-likelihood tree was recon- structed using IQ-TREE v1.6.12 (11) with the LG+I+G4+F model and 1,000 boot- strap replicates. The resulting tree (Fig. 1) suggests that strain N6 belongs to the genus Naganishia.

　Data availability. The raw reads and the draft genome have been deposited in DDBJ/ ENA/GenBank under BioProject accession number PRJDB10172 and the whole-genome shotgun project number BLZA00000000.1, respectively. The SRA accession numbers for paired-end reads are DRR237062, DRR237063, and DRR237064. The SRA accession num- bers for mate pair reads are DRR237065, DRR237066, and DRR237067. The SRA accession numbers for RNA-seq reads are DRR237068, DRR237069, DRR237070, and DRR237071. The SRA accession number for Sanger reads is DRR244395.