Next-Generation Sequencing of an 88-Year-Old Specimen of the Poorly Known Species Liagora japonica (Nemaliales, Rhodophyta) Supports the Recognition of Otohimella gen. nov.

Liagora japonica is a red algal species distributed in temperate regions of Japan. This species has not been collected from its type locality on the Pacific coast of Japan since 1927 and seems to have become extinct in this area. For molecular characterization of L. japonica, we extracted DNA from the topotype material of L. japonica collected in 1927, analyzed seven genes using Illumina next-generation sequencing, and compared these data with sequences from modern samples of similar red algae collected from the Japan Sea coast of Japan. Both morphological and molecular data from modern samples and historical specimens (including the lectotype and topotype) suggest that the specimens from the Pacific and Japan Sea coasts of Japan should be treated as a single species, and that L. japonica is phylogenetically separated from the genus Liagora. Based on the phylogenetic results and examination of reproductive structures, we propose Otohimella japonica gen. et comb. nov., characterized morphologically by diffuse carposporophytes, undivided carposporangia, and involucral filaments initiated only from the cortical cell on the supporting cell.


Introduction
Taxonomic studies are based on type specimens, which are permanently attached to taxonomic names. The type specimens of red algal seaweeds are basically specimens housed in herbaria.
To determine the correct names of modern collections, type specimens have been used for comparisons, including habit and anatomy. In many cases, morphological information from type specimens is not sufficient for taxonomic studies, because many type materials are fragmented and lack reproductive structures required for comparison with modern collections. Instead of type materials, detailed morphological observations and DNA sequencing from topotypes newly collected from their type locality have been used for taxonomic studies [1,2]. However, in some cases the coastal environment of type localities has changed and the species are missing. The most reliable method to determine the correct name of modern collections involves sequencing of type specimens and comparing them to related sequences from fieldcollected material [3,4]. However, DNA fragments of less than 300 base pairs can be obtained typically [3,4], and thus complete gene sequences required for phylogenetic analyses are not available. Many type specimens were established more than 50 years ago, and the DNA has become highly fragmented due to deterioration over time [5]; therefore, Sanger sequencing is not reliable with old DNA samples. Recently, next-generation sequencing (NGS) has been used to examine DNA sequences from old samples [6]. For red algal seaweeds, Hughey et al. [7] performed NGS to determine the complete plastid and mitochondrial genomes from 140-year-old type specimens of the bangiophycean species of Pyropia using the published complete organelle genome data from Pyropia spp. for reference mapping and sequence assembly. However, in the case of Florideophyceae, the largest group of red algae, limited genomic data are available as references in NGS, compared with bangiophycean red algae, and thus old DNA sequences have not been determined previously by NGS for molecular phylogenetic analyses.
The red algal genus Liagora is the largest genus of Liagoraceae and is widely distributed in warm temperate to tropical regions of the Atlantic, Indian, and Pacific Oceans [8]. Recent molecular analyses suggest that Liagora is polyphyletic, and the generic concept of Liagora sensu stricto has been revised substantially in recent years [9][10][11]. Although three new genera, Macrocarpus, Neoizziella, and Titanophycus, have since been separated from Liagora sensu lato [10,12], the generic positions of some species of Liagora remain unresolved. topotype material of L. japonica collected in 1927 and sequenced nuclear-encoded 28S ribosomal RNA (28S rRNA), plastid-encoded photosystem I P700 chlorophyll a apoprotein A1 (psaA), plastid-encoded ribulose-1,5-bisphosphate carboxylase/oxygenase (rbcL) genes, and the universal mitochondrial barcode marker cytochrome oxidase 1 (COI) using Illumina NGS and determined the phylogenetic position of the type material of L. japonica. We also examined herbarium specimens of lectotype and topotype as well as formalin-preserved samples used by Chihara and Yosihzaki [17] for anatomical observations. Both morphological characteristics and molecular analyses indicated that the specimens from the Pacific and Japan Sea coasts are the same species and can therefore be referred to as L. japonica. However, our molecular analyses indicated that L. japonica is incorrectly placed in Liagora and requires assignment to a new genus. We address the generic placement of L. japonica and propose a new genus, Otohimella gen. nov.

Ethics statement
We collected Liagora japonica from Sado Island and Oki Island in Japan. Collection locations and details are shown in Table B in S2 File. The collection site in Sado Island is not a protected area, while the coast of Oki Island lies within the ordinary zone of Daisen-Oki National Park. According to the National Park Act in Japan, there are no restrictions collecting marine algae from the ordinary zones of national parks. In addition, no specific permission was required for the locations on Sado Island and Oki Island.
The specimens of L. japonica housed in herbaria were subjected to morphological and molecular investigations. Official permission for this study was obtained from the Herbarium of the Graduate School of Science, Hokkaido University (SAP), and the Department of Botany, National Museum of Nature and Science, Japan (TNS).

Morphological observations
Specimens were preserved in Silica gel for DNA extraction, or 10% formalin/seawater for anatomical observations. Voucher herbarium specimens were deposited at the Department of Botany, National Museum of Nature and Science, Japan (TNS). Lectotype and topotype specimen of L. japonica, and formalin preserved sample collected from Tsushima used for Chihara and Yoshizaki [17] were also studied. Specimens on herbarium sheets, or preserved in formalin seawater were examined under an Olympus BX50 microscope (Olympus, Japan) after the material had been rehydrated and decalcified (in some cases stained with 1% aniline blue, acidified with 1% HCl and mounted in 50% aqueous Karo 1 syrup (Englewood Cliffs, NJ, USA) with 3% formaldehyde to prevent microbial growth. Drawings were made with the aid of a camera lucida.

DNA extraction and sequencing procedures for modern specimens
For phylogenetic analyses, partial 28S rRNA, rbcL genes and COI were sequenced with Sanger sequencing from two modern samples (Table B in S2 File). DNA extractions and sequencing procedures were performed in S1 Text.
DNA extraction from the historical herbarium specimens DNA was extracted from each of the herbarium specimens, ranging in age from 88 years (collected in 1927) to 57 years (collected in 1958), using a class 100 NK System Clean Bench VSF-1600RA (NK Systems, Tokyo, Japan) to prevent microbial and human contamination.
Approximately a 1 × 5 mm 2 section of material was ruptured manually using a Handy Pestle 1 (Toyobo, Osaka, Japan) in liquid nitrogen, and DNA was subsequently extracted using a QIA-GEN 1 Genomic-tip 20/G according to the manufacturer's protocol (Qiagen, Valencia, CA, USA). The DNA concentration was determined using a Quant-iT dsDNA HS assay kit with a Qubit fluorometer (Life Technologies, Carlsbad, CA, USA).
Aliquot of DNA (40-80ng) from herbarium specimens were sheared to a target peak size of 500 bp using the Covaris S220 Focused-Ultrasonicator system (Covaris, Woburn, MA, USA) according to the manufacturer's recommendations. To generate DNA sequencing libraries for high-throughput DNA sequencing, the NEBNext Ultra DNA Library Prep kit for Illumina (New England Biolabs) was used, according to the manufacturer's instructions with the exception of performing the 10-12 cycles of PCR. The amplified library products (size range, 250-600 bp) were isolated on agarose gels and purified using the NucleoSpin Gel and PCR Cleanup kit (TaKaRa, Kyoto, Japan).
The paired-end reads were generated on the Illumina MiSeq platform using the MiSeq Reagent Kit version 2 (Illumina, San Diego, CA, USA). FASTQ files were generated using the MiSeq Reporter software version 2.3.32 (Illumina). Raw sequence reads, 10,452,985 base pairs (bp) for the 57-year-old Wakayama prefecture specimen and 9,335,303 bp for the 88-year-old Kanagawa prefecture specimen, were generated.

Quality filtering of Illumina sequence data
We discarded the Illumina MiSeq reads that contained ambiguous nucleotides or were mapped to PhiX genomic sequences using Bowtie 2 version 2.2.3 with the default parameters [18]. Subsequently, we removed the adapter sequences from the reads, using Cutadapt version 1.2.135, and low-quality regions with a Phred-like quality score <17within the 3' end of the reads. In addition, we discarded reads that were < 20 bp in length or were associated with an average Phred-like quality score < 25. The quality filtering yielded high-quality reads: 10,057,049 bp for the 57-years-old specimen and 8,913,992 bp for the 88-year-old specimen.

Identification of Liagora japonica sequences
The sequences derived from the L. japonica genomes were identified based on two different methods to avoid sequencing error. The first was a read-based method and the other a scaffold-based method.
Read-based method. The high-quality MiSeq reads derived from L. japonica genomes were identified as follows. (1) An in-house nucleotide sequence database consisting of the phylogenetic marker gene sequences of Nemaliales including Liagoraceae and its relatives (designated as in-house Nemaliales database) was constructed by combining previously reported sequence data (Table C in S2 File). (2) All of the high-quality reads were subjected to BLASTN 2.2.27 searches [19] against the in-house Nemaliales database with an E-value < 0.01. (3) The reads that matched the sequences in the in-house Nemaliales database were subjected to BLASTN 2.2.27 searches against the GenBank nucleotide database (September 2014) with an E-value < 0.01. (4) The reads that matched the sequences from Rhodophyta in the GenBank nucleotide database were regarded as genome fragments of L. japonica. The target read coverages of L. japonica are shown in Tables D and E in S2 File.
MiSeq reads identified as sequences belonging to nuclear genes or organelle genes of L. japonica were used for reconstruction of each gene sequence as follows. (1) In each sample, the identified reads were assembled independently for each gene using CAP3 [20]. (2) CLUS-TALW was used with the default parameters to generate multiple alignments of the assembled contigs and singletons with reference gene sequences in the in-house Nemaliales database [21].
(3) Partially matched false-positive reads were removed manually by checking the multiple alignment. (4) In each sample, each gene sequence of L. japonica was reconstructed manually from the multiple alignment.
Scaffold-based method. The high-quality MiSeq reads of each sample were assembled using IDBA-UD version 1.1.0 with the following parameters:-mink 20-maxk 120-step 5) (Table F in S2 File) [22]. Assembled scaffolds were subjected to BLASTN 2.2.27 searches against the in-house Nemaliales database with an E-value < 0.0001. The scaffolds that matched the sequences in the in-house Nemaliales database were subjected to BLASTN 2.2.27 searches against the GenBank nucleotide database with an E-value < 0.001. The scaffolds that matched the sequences from Rhodophyta in the GenBank nucleotide database were regarded as genome fragments of Liagora japonica. The sequences similar to the reference nuclear gene sequences or putative organelle gene sequences were obtained from the scaffolds by (1) identifying a similar region within the reference gene sequence from the BLASTN results and (2) extracting a ±500 bp region of the BLASTN-aligned region. CLUSTAL W was used with default parameters to generate multiple alignments of the extracted sequences with reference gene sequences from the in-house Nemaliales database. Partially matched false-positive sequences were removed manually by checking the multiple alignment. In each sample, each gene sequence of L. japonica was manually reconstructed from the multiple alignment. The results of CLUSTALW and manual refinement of phylogenetic marker genes of two samples were checked, and the read coverage was calculated by Bowtie2 (version 2.2.6) mapping of the corresponding sample reads. The read alignments (BAM files) are available from the web server (http://liagora. paleogenome.jp).

Phylogenetic analysis
We sequenced two psaA, four rbcL, four 28S rRNA, and four COI gene sequences from modern samples and historical herbarium specimens of L. japonica (Table B in S2 File). PsaA, rbcL, and 28S rRNA genes were selected to infer the phylogeny of Liagoraceae, and the COI gene was selected to assess the effectiveness of DNA barcoding. The sequence data of Liagoraceae, available from GenBank, were compiled. PsaA sequences for 55 taxa, rbcL sequences for 99 taxa, and COI sequences for 22 taxa were aligned using CLUSTAL W [21]. The 28S rRNA gene sequences for 51 taxa of Nemaliales were aligned using CLUSTAL W [21] and were refined based on published secondary structures of the 28S rRNA gene of Palmaria palmata (Linnaeus) Kuntze [23] using SeaView 4.1 [24]. The ambiguous regions of the alignments were removed. Samples with identical nucleotide sequences were treated as a single operational taxonomic unit (OTU). As the Liagoraceae has been resolved previously as a monophyletic group [25], four to five other families belonging to the Nemaliales were designated as outgroups for psaA, rbcL, combined psaA and rbcL, and 28S rRNA gene analyses. Based on the results of psaA, rbcL, combined psaA and rbcL, and 28S rRNA gene analyses, the genera Cumagloia, Hommersandiophycus, and Nemalion were designated as outgroups for COI analysis. For confirmation of the identification of 18S rRNA, psaB, and psbA genes derived from the historical specimens of L. japonica, phylogenetic analyses based on those genes were performed. The alignments used for the present phylogenetic analyses are available from TreeBASE at https:// treebase.org/treebase-web/search/study/summary.html?id=18829 (matrix accession number S18829).
Phylogenetic analyses of the aligned sequences from each dataset, were subjected to Bayesian inference (BI) and maximum likelihood (ML) analysis. The substitution models applied to BI analyses are shown in Table G in S2 File. BI analysis was performed using MrBayes 3.2.1 [26], as described previously [27]. Four chains of Markov chain Monte Carlo (MCMC) iterations were carried out for 2,000,000 or 3,000,000 generations, keeping one tree every 500 generations. Convergence of log-likelihood and parameter values was assessed in Tracer ver-sion1.4. [28]. A burn-in sample of 5,000-7,500 trees was removed before constructing the majority rule consensus tree, and the remaining trees were used to calculate a 50% majorityrule tree and to determine the posterior probabilities (PP) of the individual branches. The Bayesian analyses are summarized in Table G in S2 File. ML analysis was performed using RAxML version7.0.4 software [29]. The GTR+I+Γ model was applied to each dataset in the analysis. Bootstrap values (BP) for ML analysis were calculated based on 1000 pseudoreplicates. The p distances and K2P genetic distances for each pair of liagoracean species were calculated using PAUP 4.0b10 [30].

Nucleotide sequence accession numbers
The Sanger sequence datasets have been submitted to DDBJ under accession numbers LC066217 to LC066223, and LC066521 to LC066532 and LC093491 to LC093498, and the Illumina sequence datasets have been submitted to the DDBJ Short Read Archive under accession number DRA003813. IDBA-UD assemblies of two samples can be accessed under BCQK01000001-BCQK01275014 (suzuki-1) and BCQL01000001-BCQL01381344 (suzuki-2), respectively.

Nomenclature acts
The electronic version of this article in Portable Document Format (PDF) in a work with an ISSN or ISBN will represent a published work according to the International Code of Nomenclature for algae, fungi, and plants, and hence the new names contained in the electronic publication of a PLOS ONE article are effectively published under that Code from the electronic edition alone, so there is no longer any need to provide printed copies.
In addition, new names contained in this work have been submitted to IPNI, from where they will be made available to the Global Names Index. The IPNI LSIDs can be resolved and the associated information viewed through any standard web browser by appending the LSID contained in this publication to the prefix http://ipni.org/. The online version of this work is archived and available from the following digital repositories: PubMed Central, LOCKSS.

Sequences determined from historical specimens
Based on limited publicly available sequence data from the Nemaliales, including Liagoraceae (Table C in S2 File), reliable sequences from Liagora japonica were determined for four genes: 28S rRNA, psaA, rbcL, and COI (Fig 1, Figures B-E in S1 File). Another three genes, nuclearencoded 18S rRNA, putative plastid-encoded photosystem I P700 chlorophyll a apoprotein A2 (psaB), and putative plastid-encoded photosystem II core 32 kDa protein (psbA) genes of potentially nemalialean affinity, were obtained from the historical specimens (Figure F in S1 File, Table B in S2 File). Both the read-based and scaffold based methods, the determined sequences of the seven genes were the same. The psaA, psaB, psbA, rbcL, and COI gene sequences were complete, while only partial sequences of the 28S rRNA and 18S rRNA genes were available. The length of the 28S and 18S rRNA gene sequences were 82.8% and 91.5%, respectively, of the published complete 28S and 18S rRNA gene sequence of P. palmata [23].

Phylogenetic markers
Analysis of the rbcL and COI sequences identified few polymorphisms (Figures B and E in S1 File) among populations of L. japonica from Misaki, type locality, Wakayama, Oki Island, and Sado Island. The rbcL gene of L. japonica showed 0-4 bp variation. The p distance among the species of the genus Liagora sensu stricto used for rbcL was 5.5-10.9%, whereas that among populations of L. japonica was < 0.3%. The COI gene for L. japonica showed 0-9 bp variation. The p distance among the species of the genus Liagora sensu stricto used for COI was 7.8-13.5%, whereas that among populations of L. japonica was < 1.5%. The p distances and K2P genetic distances among the liagoracean species used in COI analysis are shown in Tables H and I in S2 File. Both the p distances and K2P genetic distances generated similar results.
To determine whether the third nucleotides of the codons of the psaA and rbcL genes exhibit saturation of substitutions, p distances among liagoracean species based on the third nucleotide of codons were compared with those of the first and second nucleotides of the same codons ( Figure G in S1 File). The data indicated that these third nucleotide positions have been saturated with substitutions in liagoracean species, but not in L. japonica or four related genera: Akalaphycus, Macrocarpus, Neoizziella, and Titanophycus.
The topologies of individual psaA, rbcL, and combined psaA and rbcL trees, except the position of L. japonica, were basically similar to previous phylogenetic analyses for Liagoraceae [9-11, 25, 31, 32]. The phylogenetic trees generated by BI and ML showed the same topologies; therefore, we present only the BI tree topology (Fig 1, Figures B and C in S1 File). BI topologies of the individual psaA and rbcL trees (Figures B and C in S1 File) were similar to the combined psaA and rbcL data (Fig 1), except for the position of Neoizziella, but with weaker statistical support. In the individual rbcL tree, L. japonica and Macrocarpus were sisters to each other, with low support (< 0.95 PP and 59% BP), while in the individual psaA and combined psaA and rbcL trees, L. japonica and Neoizziella were sisters to each other, with moderate to high support (0.98-0.99 PP and 79-81% BP). In the combined psaA and rbcL tree, the monophyly of each genus belonging to Liagoraceae was highly supported (1.00 PP and > 95% BP), excluding Izziella. Liagora japonica, Macrocarpus, and Neoizziella formed a monophyletic clade with full support (1.00 PP and 100% BP), whereas L. japonica was clearly separated from Liagora sensu stricto.
The topology of the 28S rRNA gene tree, except the position of L. japonica was similar to previous phylogenetic analyses for Liagoraceae [10,12,33]. In the tree based on the 28S rRNA gene, Liagora japonica was separated from Liagora sensu stricto; however, the boundaries of the genera belonging to Liagoraceae were not resolved ( Figure D in S1 File).

Morphological observations of modern samples
Thalli were found on rocks at a depth of approximately 1-3 m (Fig 2A). The thalli were erect, 4-7 cm in height, and composed of 2-3 terete axes branched subdichotomously to 5-7 orders, arising from a small discoid holdfast, loosely calcified, pinkish or reddish brown in color (Fig  2A and 2B).
The thalli were multiaxial and composed of assimilatory filaments and medullary filaments (Fig 2C). The assimilatory filaments were subdichotomously branched 3-4 times. The upper parts of the assimilatory filaments were composed of ellipsoidal cells, which were 10-18 μm long and 5-8 μm wide. The lower parts were composed of elongated cells, which were 30-60 μm long and 4-6 μm wide. Gametophytes were monoecious. Spermatangial parent cells derived from apical cells and spermatangia were cut off terminally (Figs 2D and 3A). One or, rarely, two spermatangia were cut off from each spermatangial parental cell.
Carpogonial branches were curved slightly, 4-or 5-celled, and borne on the middle part of assimilatory filaments (Figs 2E and 2F, 3B and 3C). After presumed fertilization, the carpogonium divided transversely to produce a gonimoblast initial, after which the gonimoblast initial cut off primary gonimoblast cells obliquely (Fig 2G). Meanwhile, the involucral filaments were initiated mostly from cortical cells on the supporting cell (Figs 2G, 3D and 3E). Gonimoblast initials cut off second and third gonimoblast cells, and the gonimoblast developed radially (Figs 2H, 3D and 3E). At an early stage of gonimoblast development, growth of the involucral filaments and gonimoblast cells became dominant (Fig 2I). At maturity, the gonimoblast cells were embedded in and intermingled with the involucral filaments, and carposporophytes were diffuse (Figs 2J and 3F). The cells of the carpogonial branch were not fused through carposporophyte development (Figs 2G-2J and 3D-3F). Carposporangia were elliptical to oblong, 10-15 μm long and 5-6 μm wide (Fig 2K).

Morphological observations of historical materials
Based on habit and vegetative anatomies, the lectotype and topotype specimens examined corresponded to the descriptions of Yamada [13], while samples collected from Tsushima fit those of Chihara and Yoshizaki [17]. Thalli were erect and bushy, consisting of 3-5 main axes, subdichotomously branched, with 3-9 orders of branching, 7-18 cm in height, arising from a discoid holdfast with moderately calcified branches (Fig 4A and 4B). After fertilization, cells of the carpogonial branch were not fused. At an early stage of gonimoblast development, gonimoblast cells were loosely elongated, and involucral filaments intermingled with gonimoblast cells ( Fig  4C). Mature carposporophytes were not observed in the lectotype or topotype, while those of Tsushima were diffuse (Fig 4D).

Discussion
The rbcL and COI have been used for species delineation of Nemaliales, including Liagoraceae [9,11,34]. Analyses of rbcL and COI sequences suggest that modern and historical specimens collected from four localities of both the Japan Sea and Pacific sides of Japan, including topotypes, are very closely related. The typical intraspecific divergence of rbcL of the Liagoraceae based on p distances is < 0.4% [31], while that of COI of the Nemaliales is < 1.0% [34]. Liagora japonica showed < 0.3% intraspecific divergence of rbcL, which was within the range of that of Liagoraceae [31], while the intraspecific divergence of COI was < 1.5%, which was higher than typical intraspecific divergence of the Nemaliales [34]. The typical minimum interspecific divergence of COI of the Liagoraceae is 5.6% [32], and thus the intraspecific divergence of L. japonica is much lower than the interspecific divergence. Thus, all specimens should be identified as a single species: L. japonica.
However, the structures of carposporophytes from modern samples differed from those originally described by Yamada [13]. Yamada [13] presented a drawing of compact carposporophytes, whereas carposporophytes of the modern samples were diffuse. Furthermore, Chihara and Yoshizaki [17] observed L. japonica collected from Tsushima, Japan, and also presented a drawing of compact carposporophytes, similar to Yamada [13]. However, the present examination of lectotype and topotype specimens, as well as specimens used by Chihara and Yoshizaki [17], showed elongated gonimoblast cells and involucral filaments intermingled with gonimoblast cells, similar to the carposporophyte development of diffuse carposporophytes ( Table 1). The habit and vegetative structures of the modern samples collected from the Japan Sea correspond to the original description of Yamada [13]. Thus, Yamada [13] and Chihara and Yoshizaki [17] did not recognize diffuse carposporophytes, because those of the genus Liagora sensu lato were not recognized until the 1990s [35,36].
Yamada [13] Chihara and Yoshizaki [17] OJ5 and involucral filaments produced from the cortical cell on the supporting cell ( Table 2). The critical features of carposporophyte development of L. japonica are most similar to those of Yoshizakia. Yoshizakia is characterized by involucral filaments, which are not intermingled with gonimoblast filaments. Liagora japonica, however, has involucral filaments that intermingle with gonimoblast filaments. Individual psaA and rbcL, and combined psaA and rbcL analyses suggested that L. japonica, Macrocarpus, and Neoizziella form a monophyletic clade separated from Liagora sensu stricto. However, L. japonica was not included in the clade of Macrocarpus or Neoizziella. Unfortunately, psaA sequence data for the species of Yoshizakia are not yet available, but individual rbcL analyses suggest that L. japonica is clearly separated from Yoshizakia. Both morphological and molecular results suggest that L. japonica is a distinct genus in the Liagoraceae, and we propose a new genus, Otohimella Mas. Suzuki, to accommodate this species.

Taxonomic Treatment
Otohimella Mas. Suzuki gen. nov. Description: Thalli are moderately calcified, and 5 -18 cm in height, arising from a discoid holdfast with a short stipe. Thalli are subdichotomously branched to 5 or 6 orders. The cells of assimilatory filaments are ovoid to ellipsoidal and borne on colorless medullary filaments. Etymology: Named for Otohime, a sea goddess in Japanese mythology.

Conclusion
NGS can be used to determine DNA sequences from historical herbarium specimens containing highly fragmented DNA molecules. The liagoracean species are non-model organisms; thus, there are no useful complete genome sequences in the GenBank database for reference mapping and sequence assembly of NGS. However, we were able to determine seven gene sequences from historical herbarium specimens of L. japonica using NGS with limited reference sequence data. Four of the seven genes can be used for phylogenetic analyses and species identification. This study showed that sequencing ofhistorical specimens using NGS is a powerful tool for systematics and identification of not only model, but also non-model, organisms. We addressed the identity of L. japonica based on both morphological and molecular data, including those of the lectotype and topotype specimens. The species had been considered extinct on the Pacific Ocean side of Japan. However, we showed that L. japonica survives on the Japan Sea side of Japan. Further, we propose Otohimella japonica gen. et comb. nov. based on this species.
Supporting Information S1 File. Figure  The thick branches represent highest statistic supports (1.00 PP and 100% BP). Note that the trifurcation (asterisk) represents lack of bifurcation with 0.50 or more PP values in BI. Figure G. Comparson of p distances among the liagoracean species based on the first and second nucleotides of codons and based on the third nucleotide of codons in the combined psaA and rbcL dataset used for the present phylogenetic analyses (Fig 1). Red diamonds indicate p distances among "Liagora" japonica, Akalaphycus, Macrocarpus, Neoizziella, and Titanophycus. (PDF) S2 File. Table A. Collection and herbarium information for specimens of Liagora japonica used in the morphological analyses that have no molecular data. Table B. Collection locations and details, and GenBank accession numbers of samples used in the psaA, psaB, psbA, rbcL, COI, 18S rRNA, and 28S rRNA genes analyses. Table C. GenBank accession numbers of species used in the identification of Liagora japonica sequences from old herbarium specimens.