Next Article in Journal
Bithyniid Abundance in the South of Western Siberia Water-Courses and Water Reservoirs (Russia)
Next Article in Special Issue
Two New Species of the Family Canuellidae Lang, 1944 (Copepoda: Polyarthra), from Korea, with a Key to Species of the Genus Scottolana Huys, 2009
Previous Article in Journal
New Records of the Psychrophilic Tetracladium Species Isolated from Freshwater Environments in Korea
Previous Article in Special Issue
Eight New Records of Siphonophores (Cnidaria: Hydrozoa) in Korean Waters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Morphological and Molecular Studies of Three New Diatom Species from Mountain Streams in South Korea

1
Department of Environmental Science, Hanyang University, Seoul 04763, Korea
2
Department of Life Science and Research Institute for Natural Sciences, Hanyang University, Seoul 04763, Korea
*
Author to whom correspondence should be addressed.
Diversity 2022, 14(10), 790; https://doi.org/10.3390/d14100790
Submission received: 27 August 2022 / Revised: 17 September 2022 / Accepted: 19 September 2022 / Published: 23 September 2022
(This article belongs to the Collection Collection of Experts’ Researches on Aquatic Life (CEREAL))

Abstract

:
In January 2019, epilithic diatoms were collected from two streams on Mount Gumdan and Mount Yongma near Lake Paldang in central South Korea. A total of 16 diatoms were isolated and classified by molecular and morphological analysis. Morphology was studied by LM and SEM, while the molecular study focused on small subunit (SSU) rRNA and ribulose bisphosphate carboxylase (rbcL) genes. Molecular analysis showed that the three species had clear differences in phylogenetic distance. Based on these findings, we studied the ultrastructure of three species. Among the morphological characteristics, Hannaea librata is longer but narrower and always has conical spines, while the similar species H. pamirensis has bifurcated spines in the central region and conical spines near the pole. Gomphonema seminulum is wider in the axial–central area than G. pumilium. Nitzschia inclinata has a bended valve apex, while N. oligotraphenta has a straight apex.

1. Introduction

Diatoms are one of the most diverse groups of microalgae in aquatic habitats, with over 30,000 described and many more undescribed species [1]. They inhabit most of the world and have characteristics that are regulated by environmental conditions such as preferred water temperature, nutrients, conductivity, pH, and salinity [2,3,4,5], so diatoms are useful as bioindicators for interpreting the aquatic environment [6]. They have siliceous cell walls with intricate ornamentation and perforations [7,8], and their taxonomy and systematics are based on the complex patterns of their cell walls [7,9]. Therefore, reliable identification provides a consistent basis for further advances in ecological research [10].
Until the 19th century, most studies of diatom morphology were conducted using light microscopy (LM) at 1000× magnification [11,12], but LM has limitations for understanding small species and ultrastructures. Therefore, recent protocols include observing ultrastructures using scanning electron microscopy (SEM). In addition, accurate identification is achieved through molecular analysis. Among the genes used in this analysis, the small subunit (SSU) rRNA coding gene is the most widely used and suitable for inferring phylogenetic relationships [13,14], and the ribulose-1,5-bisphosphate carboxylase/oxygenase large subunit (rbcL) gene appears more suitable for evolutionary study [14].
Lake Paldang, the main source of water for citizens of Seoul and the metropolitan area, is one element of the Han River reservoir cooperation system [15,16]. Human activity is accelerating water pollution, and it is necessary to monitor the water quality of rivers and mountain streams in this area [12]. Mountain streams have harsh environmental conditions, such as low water temperature, poor nutrient status, and high turbidity, unlike rivers and lakes. In mountain streams, diatoms are the most abundant and dominant producers and are important taxa that can dominate the entire aquatic ecosystem [17]. As a result, the species composition of diatoms that inhabit the diverse environments of mountain streams varies, but many species have not yet been studied.
Therefore, this study conducted a morphological and molecular analysis of diatoms collected from streams on Mt. Gumdan and Yongma, which flow into Lake Paldang. We identified three new diatom species.

2. Materials and Methods

Sample collection: In January 2019, the epilithic diatoms were collected from streams on Mt. Gumdan (37°30′35″ N 127°16′24″ E) and Mt. Yongma (37°29′46″ N 127°16′55″ E) near Lake Paldang in central Korea (Figure 1). There was thick ice on the water surface at the time of collection, and the sediments consisted mostly of bedrock and some gravel. Each sample was collected by scraping 25 cm2 of the stone surface with a toothbrush.
Isolation: Under an inverted microscope (Olympus, Tokyo, Japan), diatoms were isolated by the capillary method [18] using Pasteur pipettes (Hilgenberg GmbH, Malsfeld, Germany). To minimize contamination by bacteria, fungi, and other algae, sterilized diatom medium (DM) [19] was dropped on the glass slide, and each cell was washed and then isolated on 96-well plates containing 160 µL of DM in each well.
Culture: At 10–14 days after isolation, when cells had reached the exponential growth stage, they were transferred to 24-well plates containing 1 mL of DM. After 1 week, cells were transferred to 25 cm3 flasks containing 20 mL DM. To maintain the health of the diatoms, the strains were subcultured at 30–50-day intervals. Samples were incubated at 20 °C under cool white fluorescent lamps with 12:12 h light:dark cycles and light intensity set to 30–50 µmol m−2 s−1 [12,20,21].
Diatom preparation for permanent slides: Morphological analysis—organic material was removed to observe morphological features. Each harvested subculture sample and acid were mixed at a 2:1 ratio (acid was a 1:1 mixture of sulfuric and nitric acid) and boiled for 2–3 minutes. To remove acid, we added distilled water, allowed samples to settle for 1 day, and carefully removed the supernatant with a pipette. This process was repeated 4–5 times per sample.
Light microscope (LM) observation: A few drops of pretreated sample were dropped onto cover glass and completely evaporated. Permanent slides were made using Mountain media (Wako Pure Chemical Industries, Ltd., Osaka, Japan) with a refraction index of greater than 1.5 for LM observation (Nikon E600, Tokyo, Japan) and photography (AmScope ToupView 3.7, Irvine, CA, USA).
Scanning electron microscope (SEM) observation: The pretreated sample was filtered using GTTP Millipore filter membrane (Millipore Filter Corporation, Cork, Ireland). The membrane was placed on an SEM stub with attached carbon tape (Shintron Enterprise Co., Ltd., Kaohsiung, Taiwan) and dried at room temperature for 24 h. Then, platinum coating was applied for 120 s using coater (Hyun corporation, Seoul, Korea), and images were acquired by SEM (Thermo fisher scientific, Waltham, MA, USA).
DNA extraction and PCR amplification: We harvested cells from the cultured flask and transferred them to 1.5 mL microtubes for centrifugation at 4000 rpm for 10 min. DNA was extracted using the DNeasy Plant mini kit (Qiagen, Valencia, CA, USA). PCR reaction mixtures of 40 μL contained 23.8 μL of distilled water, 4 μL of 10x Ex PCR buffer (TaKaRa, Tokyo, Japan), 4 μL of dNTP (TakaRa), 0.2 μL of Ex Taq polymerase (TaKaRa), 4 µL of DNA extraction template, and 2 μL of each primer. PCR was used to amplify the small subunit ribosomal RNA (SSU rRNA) and ribulose-1,5-bisphosphate carboxylase/oxygenase (rbcL) genes, as shown in Table 1. PCR assays were conducted in a Bio-Rad iCycler (Bio-Rad, Hercules, CA, USA) as follows: predenaturation at 94 °C for 4 min, 37 cycles at 94 °C for 20 s, 56 °C for 30 s, and 72 °C for 50 s, and final extension at 72 °C for 5 min. To confirm the PCR results, electrophoresis (ADVANCE, Tokyo, Japan) was conducted at 120 v using 1% agarose gel with 1% staining solution (Genetics, Dueren, Germany).
Phylogenetic analysis: DNA sequences were assembled using BioEdit v. 7.0.5.3 (Sequence Alignment Editor, Carlsbad, CA, USA, Hall 1999), and SSU rRNA and rbcL sequences were deposited in Genbank [22,23,24]. SSU and rbcL sequences of all species used in the phylogenetic analysis were taken from the National Center for Biotechnology Information (NCBI) and compared. Clustal W multiple alignment [25] was conducted in BioEdit v. 7.0.5.3 to match the sequence lengths of our and related species. MEGA version 7.0 [26] was used to calculate and represent phylogenetic relationships among the species. Phylogenetic trees were estimated by maximum likelihood based on the Kimura 2-parameter model [26]. Among the 24 models in MEGA 7.0, GTR + G + I was selected as the most appropriate. Bootstrap support was calculated with 1000 replicates for each branch of the phylogenetic tree. To calculate the similarity score and genetic distance (p-distance), 1000 bootstrap replicates and the Kimura 2-parameter model were also used in BioEdit v. 7.0.5.3 and MEGA 7.0.

3. Results and Discussion

3.1. Morphological Characteristics of Hannaea librata sp. nov.

Class Bacillariophyceae
Subclass Fragilariophycidae
Order Licmophorales
Family Ulnariaceae
Genus Hannaea
Hannaea librata E.A. Hwang and B.H. Kim (Figure 2, LM; Figure 3, SEM)
Description: Valves are linear with strongly rostrate apices, slightly arcuate, and slight swelling at the central area of the curved (ventral) side. Length 44–99 μm, width 5–5.5 μm, 12–16 striae in 10 μm, and 70–80 areolae in 10 μm. Axial area is narrow, linear, and slightly bent in the central region. Striae are alternate, uniseriate, perpendicular to the axial area, and parallel to each other, but one or two striae at the apical end were parallel to the axis. Apical pore field is located at each pole of the valve. Valve face is undulating; valve with striae is slightly sunken, and without striae is swollen. In central area, absent striae but the undulating valve face forms ghost striae and forms a wide U-shape extending into both striae. Valve mantle is flat regardless of striae. Areolae are poroid type, elongated oval in shape, and near the axial and margin are round in shape. Single rimoportula with slit-like opening located at valve apices. Girdle band is open ring type with single row of round areolae and scalloped advalvar edge. Spines present along the valve face–mantle junction until apical, shape conical in all parts of the valve, and tips radiate from the central area. In girdle view, frustules are rectangular. Cells form linear colonies.
Remarks: To compare the morphological characteristics of Hannaea librata with similar species, its structure was observed using LM, and its ultrastructure was closely observed using SEM. The morphological characteristics of H. librata are different from those of related taxa (Table 2). The valve of H. pamirensis [27] is shorter and wider than that of H. librata, and the density of striae is higher. The central area of H. librata is wide and U-shaped, extending into both striae and wider than H. pamirensis. Moreover, the spine of H. librata is conical, but H. pamirensis has conical spines near each pole and bifurcated thorn-shaped spines near the center of the valve, showing a distinct difference. H. hattoriana [28] is another species with similar morphological characteristics and has a more lanceolate valve shape than H. librata.; H hattoriana has a narrower central area than H. librata. The valve shape of H. recta [29,30] is more lanceolate than H. librata. The valve mantle of H. librata is flat, while that of H. recta is waved. Lastly, the shape of the spine of H. librata is constant, whereas the shape of the H. recta changes from the center valve and each pole. Hannaea shares morphological similarities with both Fragilaria and Synedra [31] and is most closely related to Fragilaria [32]. In comparison with F. capucina, which has the most similar morphology, the shape spine in the center of the valve is conical in H. librata and triangular in F. capucina. In the center of the valve, H. librata has unilateral swelling, but F. capucina remains flat from both sides [33,34,35]. Hannaea librata is distinguished from similar related taxa in the shape of the valve, central area, and spine.
Holotype: permanent slides were deposited with the Korea Collection for Type Cultures KCTC under the deposit number AG9001.
Isotype: permanent slides were deposited with the Korea Collection for Type Cultures KCTC under the deposit numbers AG9002, AG9003, AG9006, and AG9007.
Habitat: free-living on rocky substrates. Water quality: WT (water temperature), 1.0 °C; DO (dissolved oxygen), 15.6 mg/L; pH, 6.4; EC (electrical conductivity), 42 μS/cm; turbidity, 8.5 NTU. Major species (more than 5% of the total abundance) recorded from the same location: Achnanthes convergens, Gomphonema parvulum, Achnanthes minutissima, Fragilaria gracilis, and Gomphonema gracile.
Type locality: Republic of Korea, Gyeonggi-do Province, Mount Gumdan, 37°30′35″ N 127°16′24″ E, 3 January 2019.
Molecular characterization: SSU, rbcL accession no. The nucleotide sequences of the SSU rRNA and rbcL genes of the strain were deposited in GenBank (NCBI) under the SSU accession numbers ON040635, ON040636, ON040640, ON040641, and ON040642, and rbcL accession numbers OP137173, OP183209, OP183212, OP183213, and OP183214.
Etymology:librata” is derived from the Latin “librátĭo”, which means that spines at the center and at each apex are constant in shape.

3.2. Phylogenetic Characteristics of H. librata sp. nov.

We selected several species that were similar molecularly to Hannaea using GenBank and inferred phylogenetic characteristics using 28 aligned sequences of SSU and 23 aligned sequences of rbcL (Figure 4 and Figure 5). Hannea has been separated from Ceratoneis, and there are 20 species worldwide. However, there are few molecular phylogenetic studies, and the sequence is unknown. In the SSU and rbcL tree, H. librata collected in different places belong to the same clade. Its sister groups include Fragilaria, Ulnaria, and Synedra, which are distinguished from Hannaea. Genetic distance and similarity analysis results for SSU and rbcL indicated that the genetic distance between H. librata collected from different places was the same for all five individuals (Table 3 and Table 4). H. librata showed a genetic distance of 0.015 or less from Fragilaria in SSU genes and 0.081 or less for the rbcL gene and appeared molecularly similar to F. capucina, with a genetic distance for the SSU gene of 0.011 and 0.012 for the rbcL gene, making it clearly distinct from H. librata. Therefore, the phylogenetic features of H. librata differentiate it from Fragilaria, occupying distinctly different clades, and are clearly distinguished from those of F. capucina, which is molecularly and phylogenetically similar.

3.3. Morphological Characteristics of Gomphonema Seminulum sp. nov.

Class Bacillariophyceae
Subclass Bacillariophycidae
Order Cymbellales
Family Gomphonemataceae
Genus Gomphonema
Gomphonema seminulum E.A. Hwang and B.H. Kim (Figure 6, LM; Figure 7, SEM)
Description: Valves are narrow lanceolate, slightly asymmetric about the axis, valve apices are bluntly rounded, and central and axial area is wide lanceolate shape. Length 15–20 μm, width 4–5 μm, 13–15 striae in 10 μm, and 40–50 areolae in 10 μm. Striae are uniseriate and parallel to slightly radiate toward the apices. Areola is broadly C-shaped, occluded by raised flaps; 2–4 lines from the axis are flapped open toward the center, and the other lines are open in the opposite direction. The large openings of the alveolus are formed inside the valve, and the areola is located outside the valve along this ultrastructure. Single stigma is circular in shape, located at the center of valve. Raphe extends along the entire valvar face and is slightly undulate. External raphe ending, with central endings terminating in elliptical central pores, and polar endings deflecting in the opposite direction of the stigma. Internal raphe ending, with central endings bent to stigma and hook-shaped toward each pole, and polar endings forming helictoglossa. The apical pore field is located at the foot pole of the valve and is bifurcated by the raphe. Pseudosepta is formed at each apex of the internal valve. In girdle view, frustules are wedge-shaped.
Remarks: In this study, we conducted LM and SEM to study the morphology of G. seminulum and confirmed new species belonging to the genus Gomphonema. The morphology of G. seminulum differs from closely related taxa (Table 5); among them, G. bourbonense [36], which lacks descriptions of ultrastructure, was further analyzed by examining type material (specimen ID TCC441, TCC458, TCC930 in the Thonon Culture Collection, France). The G. seminulum valve is narrower and more lanceolate than that of G. bourbonense. In addition, the valve of G. seminulum is asymmetric with a very wide central area, whereas G. bourbonense is symmetric with a narrower central area. Finally, C-shaped areolae in G. seminulum are wider than in G. bourbonense. It shows distinct differences from G. pumilum, another species that is morphologically similar to G. seminulum. G pumilum striae slightly radiate to strongly radiate toward the apices, but G. seminulum is almost parallel in the central area, slightly radiate near the head pole, and strongly radiate toward the foot pole. The valve of the head pole in G. pumilum is strongly tapered compared with G. seminulum. G angustum has a significantly lower central area density of striae than G. seminulum. Lastly, in G. clevei, the valve is rhombic or elliptical clavate, whereas G. seminulum is narrow lanceolate, and the external polar raphe ending of G. clevei is curved toward the stigma but deflected in the opposite direction in G. seminulum. Based on these observations, G. seminulum is clearly distinguished from closely related taxa.
Holotype: permanent slides were deposited with the Korea Collection for Type Cultures KCTC under the deposit number AG9005.
Habitat: free-living on rocky substrates. Water quality: WT, 1.3 °C; DO, 14.6 mg/L; pH, 5.5; EC, 49 μS/cm; turbidity, 10.7 NTU. Major species (more than 5% of the total abundance) recorded from the same location: Gomphoneis quadripunctata, Meridion constrictum, Achnanthidium minutissimum, and Achnanthes convergens.
Type locality: Republic of Korea, Gyeonggi-do Province, Mount Yongma, 37°29′46″ N 127°16′55″ E, 3 January 2019.
Molecular characterization: SSU, rbcL accession no. The nucleotide sequences of the SSU rRNA and rbcL genes of the strain were deposited in GenBank (NCBI) under the accession numbers ON040637 and OP183210.
Etymology:seminulu” is derived from the Latin “seminúdus”, which means that the central area of G. seminulum is empty and simple.

3.4. Phylogenetic Characteristics of Gomphonema Seminulum sp. nov.

We used 20 aligned sequences of SSU rRNA and 18 aligned sequences of rbcL from GenBank and inferred the phylogenetic characteristics of species (Figure 8 and Figure 9). In the SSU rRNA and rbcL trees, G. seminulum was placed in the Gomphonema clade and was distant from other similar genera. In the Gomphonema clade, G. seminulum was closer to the clade in which G. bourbonense, G. pumilum, G. truncatum, G. subclavatum, and G. acuminatum are placed, and its sister taxon was G. bourbonense. However, they were not close to G. parvulum, G. affine and G. clevei. Among the SSU rRNA alignment sequences (Table 6), G. bourbonense showed the smallest genetic distance from G. seminulum at 0.007. However, among morphologically similar species, G. pumilum was 0.008, and G. clevei was 0.031, which were greater distances. In rbcL-aligned sequences (Table 7), G. acuminatum showed the smallest genetic distance from G. seminulum of 0.037, while G. bourbonense, located in a nearby clade, was 0.044, and G. pumilum was 0.049. In the SSU rRNA and rbcL phylogenetic analysis, G. seminulum and G. bourbonense were close genetically and located in the neighboring clades, but the molecular phylogenetic characteristics of the two species were distinct. In addition, within the Gomponema clade, morphologically similar species such as G. clevei were clearly distinguished in a different clade. As a result, morphological and phylogenetic characteristics of G. seminulum clearly distinguish it from similar species.

3.5. Morphological Characteristics of Nitzschia inclinata sp. nov.

Class Bacillariophyceae
Subclass Bacillariophycidae D.G.
Order Bacillariales
Family Bacillariaceae
Genus Nitzschia
Nitzschia inclinata E.A. Hwang and B.H. Kim (Figure 10, LM; Figure 11, SEM)
Description: Valves are narrow linear, and valve apices are capitate and slightly bent to one side. Length 32–34.5 μm, width 3.5–5 μm, 7–10 fibula in 10 μm, 50 striae in 10 μm, and 60–80 areolae in 10 μm. Striae are uniseriate, parallel, and not observed under LM. Areola on the valve is loculate in type and circular in shape. The terminal fissure is located at each valve apex with a broad sweeping curve, and a helictoglossa is formed inside the valve. Fibula is to one side of the valve face and is irregular in width. Conopeum is positioned on the external valve face along the fibula. On the face of the conopeum, circular areolae form one line, but for the valve apices, the number of areolae increases to 1–3. Raphe is eccentric from the valve face, centered on the conopeum. Support points connecting the valve face to the conopeum are observed. In girdle view, frustules are rectangular. Girdle bands with double rows of small punctuate are associated with a valve.
Remarks: The morphology, including the ultrastructure, of Nitzschia inclinata was observed and compared with morphologically similar species belonging to the same genus (Table 8). The N. dissipata [41] valve is more lanceolate than N. inclinata.; N. dissipata formed subrostrate apices but N. inclinata has capitate apices. Another similar species, N. oligotraphenta [42], has linear lanceolate valves, but N. inclinata is narrow and linear in shape. A more obvious difference is that the valve apices of N. oligotraphenta are straight, whereas those of N. inclinata are slightly curved to one side of the valve. N. sigmoidea [43] is longer and larger than N. inclinata. In addition, in girdle view, N. sigmoidea is sigmoid, and N. inclinata is rectangular. Finally, N. sigmoidea fibula is along the margin of one side of the valve, but N. inclinata fibula is located on the valve face. The valve of N. angularis [44] is rhombic in shape, so the center of the valve is swollen and sharply narrowed toward the end of the apices. N. inclinata is clearly distinct from similar species in morphology and microstructure.
Holotype: permanent slides were deposited with the Korea Collection for Type Cultures KCTC under the deposit number AG9004.
Habitat: free-living on rocky substrates. Water quality: WT, 1.3 °C; DO, 14.6 mg/L; pH, 5.5; EC, 49 μS/cm; turbidity, 10.7 NTU. Major species (more than 5% of the total abundance) recorded from the same location: Gomphoneis quadripunctata, Meridion constrictum, Achnanthidium minutissimum, and Achnanthes convergens.
Type locality: Republic of Korea, Gyeonggi-do Province, Mount Yongma, 37°29′46″ N 127°16′55″ E, 3 January 2019.
Molecular characterization: SSU, rbcL accession no. The nucleotide sequences of the SSU rRNA and rbcL genes of the strain were deposited in GenBank (NCBI) under the accession numbers ON040638 and OP183211.
Etymology:inclinata” is derived from the Latin “incliato”, which means that the valve apices of N. inclinata are inclined to one side.

3.6. Phylogenetic Characteristics of Nitzschia inclinata sp. nov.

We performed phylogenetic analysis by selecting 20 aligned sequences of SSU rRNA and 18 aligned sequences of rbcL from the molecularly similar species to N. inclinata, and Cyclotella meneghiniana was used as an outgroup (Figure 12 and Figure 13). In the SSU rRNA and rbcL phylogenetic trees, N. inclinata was located in the Nitzschia clade and was distinct from Navicula and Amphora. In SSU rRNA tree, N. inclinata was related to N. dissipata and N. sigmoidea, with strong support (ML bootstrap = 100%). The genetic distance between N. inclinata and these species was 0.016 for N. dissipata and 0.019 for N. simoidea, located within the same clade but distinct (Table 9). In the rbcL tree, N. inclinata was closely related to N. sigmoidea and N. dissipata and was the sister taxon of N. sigmoidea with 64% support. In rbcL-aligned sequences, the genetic distance between N. inclinata and N. sigmodea was 0.027, closer than to other species (Table 10). In the SSU rRNA and rbcL phylogenetic trees, N. inclinata was located on a branch close to N. dissipata and N. sigmodea, but its individual branch had high support values. Although they were close in genetic distance, they showed clear differences. Therefore, we are confident that N. inclinata is a new species with different molecular characteristics from similar species.

Author Contributions

E.-A.H. and B.-H.K. conceived the research; E.-A.H., H.-K.K., I.-H.C., and B.-H.K. performed fieldwork; E.-A.H., B.-H.K. and C.Y. analyzed the data; E.-A.H. and B.-H.K. wrote and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a grant from the Nakdonggang National Institute of Biological Resources (NNIBR) funded by the Ministry of Environment (MOE) of the Republic of Korea (NNIBR202201205).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank anonymous reviewers for their valuable and constructive comments, and special thanks to Kim, Y.-H. for kind assistance.

Conflicts of Interest

The authors report no potential conflict of interest.

References

  1. Mann, D.G.; Vanormelingen, P. An inordinate fondness? The number, distributions, and origins of diatom species. J. Eukaryot. Microbiol. 2013, 60, 414–420. [Google Scholar] [CrossRef] [PubMed]
  2. Kovács, C.; Kahlert, M.; Padisák, J. Benthic diatom communities along pH and TP gradients in Hungarian and Swedish streams. J. Appl. Phycol. 2006, 18, 105–117. [Google Scholar] [CrossRef]
  3. Pestryakova, L.A.; Herzschuh, U.; Gorodnichev, R.; Wetterich, S. The sensitivity of diatom taxa from Yakutian lakes (north-eastern Siberia) to electrical conductivity and other environmental variables. Polar Res. 2018, 37, 1485625. [Google Scholar] [CrossRef]
  4. Kim, H.K.; Cho, I.H.; Hwang, E.A.; Kim, Y.J.; Kim, B.H. Benthic diatom communities in Korean estuaries: Species appearances in relation to environmental variables. Int. J. Environ. Health Res. 2019, 16, 2681. [Google Scholar] [CrossRef] [PubMed]
  5. Pan, Y.; Lowe, R.L. Independent and interactive effects of nutrients and grazers on benthic algal community structure. Hydrobiologia 1994, 291, 201–209. [Google Scholar] [CrossRef]
  6. Pan, Y.; Stevenson, R.J.; Hill, B.H.; Herlihy, A.T.; Collins, G.B. Using diatoms as indicators of ecological conditions in lotic systems: A regional assessment. J. N. Am. Benthol. Soc. 1996, 15, 481–495. [Google Scholar] [CrossRef]
  7. Round, F.E.; Crawford, R.M.; Mann, D.G. The Diatoms: Biology and Morphology of the Genera; Cambridge University Press: Cambridge, UK, 1990; 747p. [Google Scholar]
  8. Green, J.C.; Leadbeater, B.S.C.; Diver, W.L. (Eds.) The origins of the diatom and its life cycle. In The Chromophyte Algae: Problems and Perspectives; Clarendon Press: Oxford, UK, 1989; pp. 305–321. [Google Scholar]
  9. Witkowski, A.; Sieminska, J. (Eds.) A review of the evolution of the diatoms—A total approach using molecules, morphology and geology. In The Origin and Early Evolution of the Diatoms: Fossil, Molecular and Biogeographical Approaches; Polish Academy of Sciences: Krakow, Poland, 2000; pp. 13–35. [Google Scholar]
  10. Cox, E.J. The identity and typification of some naviculoid diatoms (Bacillariophyta) from freshwater or brackish habitats. Phycologia 1998, 37, 162–175. [Google Scholar] [CrossRef]
  11. Morales, E.A.; Siver, P.A.; Trainor, F.R. Identification of diatoms (Bacillariophyceae) during ecological assessments: Comparison between light microscopy and scanning electron microscopy techniques. Proc. Acad. Nat. Sci. Phila. 2001, 151, 95–103. [Google Scholar] [CrossRef]
  12. Tan, L.; Wang, P.; Cho, I.H.; Hwang, E.A.; Lee, H.; Kim, B.H. Morphology and phylogenetic position of three new raphid diatoms (Bacillariophyceae) from Hangang River, South Korea. Phytotaxa 2020, 442, 153–182. [Google Scholar] [CrossRef]
  13. Auwera, G.V.D.; Wachter, R.D. Structure of the large subunit rDNA from a diatom, and comparison between small and large subunit ribosomal RNA for studying stramenopile evolution. J. Eukaryot. Microbiol. 1998, 45, 521–527. [Google Scholar] [CrossRef]
  14. Mann, D.G.; Simpson, G.E.; Sluiman, H.J.; Möller, M. rbcL gene tree of diatoms: A second large data-set for phylogenetic reconstruction. Phycologia 2001, 40, 1–2. [Google Scholar]
  15. Kim, H.P. A Study on Ecological Model to Estimate Water Quality in the Lake Pal’tang; Hanyang University: Seoul, Korea, 1999. [Google Scholar]
  16. Choi, H.J.; Cho, Y.C.; Yu, S.; Song, Y.S.; Ryu, I. Analysis of Water Circulation Characteristics for Hydraulic and Water Temperature Investigation in Paldang Reservoir. Ecol. Resilient Infrastruct. 2019, 6, 12–22. [Google Scholar]
  17. Tsoi, W.Y.; Hadwen, W.L.; Sheldon, F. Identifying diatom indicator species of nutrient enrichment: An in situ nutrient enrichment experiment in subtropical upland streams. Ecol. Indic. 2020, 119, 106744. [Google Scholar] [CrossRef]
  18. Andersen, R.A. (Ed.) Algal Culturing Techniques; Elsevier: San Diego, CA, USA, 2005; p. 578. [Google Scholar]
  19. Beakes, G.W.; Canter, H.M.; Jaworski, G.H.M. Zoospore ultrastructure of Zygorhizidium affluens and Z. planktonicum, two chytrids parasitizing the diatom Asterionella formosa. Can. J. Bot. 1988, 66, 1054–1067. [Google Scholar] [CrossRef]
  20. Miao, M.; Li, Z.; Hwang, E.A.; Kim, H.K.; Lee, H.; Kim, B.H. Two New Benthic Diatoms of the Genus Achnanthidium (Bacillariophyceae) from the Hangang River, Korea. Diversity 2020, 12, 285. [Google Scholar] [CrossRef]
  21. Shi, Y.; Wang, P.; Kim, H.K.; Lee, H.; Han, M.S.; Kim, B.H. Lemnicola hungarica (Bacillariophyceae) and the new monoraphid diatom Lemnicola uniseriata sp. nov. (Bacillariophyceae) from South Korea. Diatom Res. 2018, 33, 69–87. [Google Scholar] [CrossRef]
  22. Ki, J.S.; Han, M.S. Molecular analysis of complete ssu to lsu rdna sequence in the harmful dinoflagellate Alexandrium tamarense (Korean isolate, HY970328M). Ocean Sci. 2005, 40, 43–54. [Google Scholar] [CrossRef]
  23. Medlin, L.; Elwood, H.J.; Stickel, S.; Sogin, M.L. The characterization of enzymatically amplified eukaryotic 16S-like rRNA-coding regions. Gene 1988, 71, 491–499. [Google Scholar] [CrossRef]
  24. Bruder, K.; Medlin, L. Molecular assessment of phylogenetic relationships in selected species/genera in the naviculoid diatoms (Bacillariophyta). I. The genus Placoneis. Nova Hedwig. 2007, 85, 331. [Google Scholar] [CrossRef]
  25. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar] [CrossRef]
  26. Kumar, S.; Stecher, G.; Tamura, K. MEGA7 (2016): Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  27. Liu, Q.; Glushchenko, A.; Kulikovskiy, M.; Maltsev, Y.; Kociolek, J.P. New Hannaea Patrick (Fragilariaceae, Bacillariophyta) species from Asia, with comments on the biogeography of the genus. Cryptogam. Algol. 2019, 40, 41–61. [Google Scholar] [CrossRef]
  28. Meister, F. Beiträge zur Bacillariaceenflora Japan. Arch. Hydrobiol. 1914, 9, 226–232. [Google Scholar]
  29. Skvortzow, B.W. A contribution to the diatoms of Baikal Lake. Proc. R. Soc. B 1928, 1, 3. [Google Scholar]
  30. Kobayasi, H.; Idei, M.; Mayama, S.; Nagumon, T.; Osada, K. (Eds.) Kobayasi’s Atlas of Japanese Diatoms Based on Electron Microscopy; Unchida Rokakuho: Tokyo, Japan, 2006; p. 531. [Google Scholar]
  31. Bixby, R.J. Morphology, Phytogeography, and Systematics of the Diatom Genus Hannaea (Bacillariophyceae); University of Michigan: Ann Arbor, MI, USA, 2001. [Google Scholar]
  32. Williams, D.M. Tabulariopsis, a new genus of marine araphid diatom, with notes on the taxonomy of Tabularia (Kütz.) Williams et Round. Nova Hedwig. 1988, 47, 247–254. [Google Scholar]
  33. Lange-Bertalot, H. Zur systematischen Bewertung der bandformigen Kolonien von Navicula und Fragilaria. Kriterien fur die Vereinigung von Synedra (subgen. Synedra) Ehrenberg mit Fragilaria Lyngbye. Nova Hedwig. 1980, 33, 723–787. [Google Scholar]
  34. Tuji, A.; Williams, D.M. Typification of Conferva pectinalis OF Müll. (Bacillariophyceae) and the identity of the type of an alleged synonym, Fragilaria capucina Desm. Taxon 2006, 55, 193–199. [Google Scholar] [CrossRef]
  35. Heudre, D.; Wetzel, C.E.; Moreau, L.; Van de Vijver, B.; Ector, L. On the identity of the rare Fragilaria subconstricta (Fragilariaceae), with Fragilaria species forming ribbon-like colonies shortly reconsidered. Plant Ecol. Evol. 2019, 152, 327–339. [Google Scholar] [CrossRef]
  36. Reichardt, E. Taxonomische revision des artenkomplexes um Gomphonema pumilum (Bacillariophyceae). Nova Hedwig. 1997, 99–129. [Google Scholar] [CrossRef]
  37. Reichardt, E.; Lange-Bertalot, H. Taxonomische Revision des Artenkomplexes um Gomphonoma angustumG. dichotomumG. intricatumG. vibrio und ähnliche Taxa (Bacillariophyceae). Nova Hedwig. 1991, 53, 519–544. [Google Scholar]
  38. Wojtal, A. Diatoms of the genus Gomphonema Ehr. (Bacillariophyceae) from a karstic stream in the Krakowsko-Czestochowska Upland. Acta Soc. Bot. Pol. 2003, 72, 213–220. [Google Scholar] [CrossRef]
  39. Agardh, C.A. (Ed.) Conspectus Criticus Diatomacearum. Part 3; Litteris Berlingianis: Lundae, Sweden, 1831; pp. 33–48. [Google Scholar]
  40. Fricke, F. (Ed.) Verzeichniss der in A. Schmidt’s Atlas der Diatomaceenkunde Tafel 1–240 (Serie IV) Abgebildeten und Bekannten Formen/Herausgegeben von Friedr; O.R. Reisland: Leipzig, Germany, 1902. [Google Scholar]
  41. Rabenhorst, L. (Ed.) Die Algen Sachsens. Resp. Mittel-Europa’s Gesammelt und Herausgegeben von Dr. L.; Exsiccata: Dresden, Germany, 1860. [Google Scholar]
  42. Lange-Bertalot, H. (Ed.) Indicators of Oligotrophy. 800 Taxa Representative of Three Ecologically Distinct Lake Types; Koeltz Scientific Books: Königstein, Germany, 1996; 390p. [Google Scholar]
  43. Knattrup, A.; Yde, M.; Lundholm, N.; Ellegaard, M. A detailed description of a Danish strain of Nitzschia sigmoidea, the type species of Nitzschia, providing a reference for future morphological and phylogenetic studies of the genus. Diatom Res. 2007, 22, 105–116. [Google Scholar] [CrossRef]
  44. Smith, W. (Ed.) A Synopsis of the British Diatomaceae; with Remarks on Their Structure, Function and Distribution; and Instructions for Collecting and Preserving Specimens; Paternoster Row: London, UK, 1853; Volume 1, pp. 1–89. [Google Scholar]
  45. Joh, G. (Ed.) Chrysophyta: Bacillariophyceae: Pennales: Raphidineae: Naviculaceae: Cymbella, Cymbopleura, Encyonema, Encyonopsis, Reimeria, Gomphonema. Freshwater diatoms IV. In Algal Flora of Korea; National Institute of Biological Resources: Incheon, Korea, 2011; Volume 3, pp. 1–70. [Google Scholar]
  46. Lange-Bertalot, H. (Ed.) Diatomeen im Süßwasser—Benthos von Mitteleuropa. Bestimmungsflora Kieselalgen für die Ökologische Praxis. Über 700 der Häufigsten Arten und Ihre Ökologie, Königstein; Koeltz Scientific Books: Königstein, Germany, 2013; pp. 1–908. [Google Scholar]
  47. In Diatoms of North America. Available online: https://diatoms.org/species/nitzschia_dissipata (accessed on 15 August 2022).
Figure 1. Map of sampling sites in mountain streams from Korea.
Figure 1. Map of sampling sites in mountain streams from Korea.
Diversity 14 00790 g001
Figure 2. Light microscope photo of Hannaea librata sp. nov. (AL) valve view. (AG) taken from isotype population (AG9002); (HL) from holotype population (AG001); scale bar = 10 µm.
Figure 2. Light microscope photo of Hannaea librata sp. nov. (AL) valve view. (AG) taken from isotype population (AG9002); (HL) from holotype population (AG001); scale bar = 10 µm.
Diversity 14 00790 g002
Figure 3. SEM of Hannaea librata sp. nov. taken from holotype population (AG9001), scale bar = 2 μm: (A) external view of apex with rimoportula (arrow); (B) external view of central area, ghost strain (arrow); (C) external view without rimoportula; (D) external view of one side showing valve apices with spines, conical shape (solid arrow), girdle band (dotted arrow); (E) internal view of apex with rimoportula (arrow); (F) internal view of central area; (G) internal view of apex without rimoportula, vertical striae (arrow); (H) external view of central area with spines, conical shape (arrow).
Figure 3. SEM of Hannaea librata sp. nov. taken from holotype population (AG9001), scale bar = 2 μm: (A) external view of apex with rimoportula (arrow); (B) external view of central area, ghost strain (arrow); (C) external view without rimoportula; (D) external view of one side showing valve apices with spines, conical shape (solid arrow), girdle band (dotted arrow); (E) internal view of apex with rimoportula (arrow); (F) internal view of central area; (G) internal view of apex without rimoportula, vertical striae (arrow); (H) external view of central area with spines, conical shape (arrow).
Diversity 14 00790 g003
Figure 4. SSU rRNA maximum likelihood phylogenetic tree of Hannaea librata sp. nov. and other species. The tree with the highest log likelihood (−3434.22) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis included 28 nucleotide sequences and 1478 positions in the final data set. Cyclotella meneghiniana was used as an outgroup.
Figure 4. SSU rRNA maximum likelihood phylogenetic tree of Hannaea librata sp. nov. and other species. The tree with the highest log likelihood (−3434.22) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis included 28 nucleotide sequences and 1478 positions in the final data set. Cyclotella meneghiniana was used as an outgroup.
Diversity 14 00790 g004
Figure 5. rbcL maximum likelihood phylogenetic tree of Hannaea librata sp. nov. and other species. The tree with the highest log likelihood (−1639.58) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis involved 22 nucleotide sequences and 596 positions in the final data set. Cyclotella meneghiniana was used as the outgroup.
Figure 5. rbcL maximum likelihood phylogenetic tree of Hannaea librata sp. nov. and other species. The tree with the highest log likelihood (−1639.58) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis involved 22 nucleotide sequences and 596 positions in the final data set. Cyclotella meneghiniana was used as the outgroup.
Diversity 14 00790 g005
Figure 6. Light microscope photo of Gomphonema seminulum sp. nov. taken from holotype population (AG9005): (A) girdle view; (BK) valve view; scale bar = 10 µm.
Figure 6. Light microscope photo of Gomphonema seminulum sp. nov. taken from holotype population (AG9005): (A) girdle view; (BK) valve view; scale bar = 10 µm.
Diversity 14 00790 g006
Figure 7. SEM of Gomphonema seminulum sp. nov. taken from holotype population (AG9005): (A,B) scale bar = 5 μm; (CH) scale bar = 2 μm; (A) external view of whole valve; (B) internal view of whole valve; (C) external view of head pole; (D) external view of valve middle, showing the stigma on central area of valve face (arrow); (E) internal view of foot pole, showing the apical pore field (arrow); (F) internal view of head pole, showing the helictoglossa (arrow), pseudosepta (arrow P); (G) internal view of valve middle, showing the alveoli (arrow); (H) internal view of foot pole, showing the pseudosepta (arrow).
Figure 7. SEM of Gomphonema seminulum sp. nov. taken from holotype population (AG9005): (A,B) scale bar = 5 μm; (CH) scale bar = 2 μm; (A) external view of whole valve; (B) internal view of whole valve; (C) external view of head pole; (D) external view of valve middle, showing the stigma on central area of valve face (arrow); (E) internal view of foot pole, showing the apical pore field (arrow); (F) internal view of head pole, showing the helictoglossa (arrow), pseudosepta (arrow P); (G) internal view of valve middle, showing the alveoli (arrow); (H) internal view of foot pole, showing the pseudosepta (arrow).
Diversity 14 00790 g007
Figure 8. SSU rRNA phylogenetic tree by maximum likelihood method of Gomphonema seminulum sp. nov. and other species molecular position. The tree with the highest log likelihood (−3766.04) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis involved 20 nucleotide sequences and 1475 positions in the final data set. Cyclotella meneghiniana was used as the outgroup.
Figure 8. SSU rRNA phylogenetic tree by maximum likelihood method of Gomphonema seminulum sp. nov. and other species molecular position. The tree with the highest log likelihood (−3766.04) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis involved 20 nucleotide sequences and 1475 positions in the final data set. Cyclotella meneghiniana was used as the outgroup.
Diversity 14 00790 g008
Figure 9. rbcL phylogenetic tree by maximum likelihood method of Gomphonema seminulum sp. nov. and other species molecular position. The tree with the highest log likelihood (−1964.35) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis involved 18 nucleotide sequences and 596 positions in the final data set. Cyclotella meneghiniana was used as the outgroup.
Figure 9. rbcL phylogenetic tree by maximum likelihood method of Gomphonema seminulum sp. nov. and other species molecular position. The tree with the highest log likelihood (−1964.35) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis involved 18 nucleotide sequences and 596 positions in the final data set. Cyclotella meneghiniana was used as the outgroup.
Diversity 14 00790 g009
Figure 10. Light microscope photo of Nitzschia inclinata sp. nov. taken from holotype population (AG9004): (A) girdle view; (BN) valve view; scale bar = 10 µm.
Figure 10. Light microscope photo of Nitzschia inclinata sp. nov. taken from holotype population (AG9004): (A) girdle view; (BN) valve view; scale bar = 10 µm.
Diversity 14 00790 g010
Figure 11. SEM of Nitzschia inclinata sp. nov. taken from holotype population (AG9004): (A,B) scale bar = 10 μm; (CG) scale bar = 2 μm; (H) scale bar = 200 nm; (A) external view of whole valve; (B) internal view of whole valve; (C) external view of valve apex, terminal fissure (arrow); (D) external view of valve middle; (E) internal view of valve apex with girdle band (arrow); (F) internal view of valve apex, helictoglossa (arrow); (G) internal view of valve middle, fibulae (arrow); (H) external view of areolae.
Figure 11. SEM of Nitzschia inclinata sp. nov. taken from holotype population (AG9004): (A,B) scale bar = 10 μm; (CG) scale bar = 2 μm; (H) scale bar = 200 nm; (A) external view of whole valve; (B) internal view of whole valve; (C) external view of valve apex, terminal fissure (arrow); (D) external view of valve middle; (E) internal view of valve apex with girdle band (arrow); (F) internal view of valve apex, helictoglossa (arrow); (G) internal view of valve middle, fibulae (arrow); (H) external view of areolae.
Diversity 14 00790 g011
Figure 12. SSU rRNA phylogenetic tree by maximum likelihood method of Nitzschia inclinata sp. nov. and other species molecular position. The tree with the highest log likelihood (−4374.77) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis involved 20 nucleotide sequences and 1392 positions in the final data set. Cyclotella meneghiniana was used as the outgroup.
Figure 12. SSU rRNA phylogenetic tree by maximum likelihood method of Nitzschia inclinata sp. nov. and other species molecular position. The tree with the highest log likelihood (−4374.77) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis involved 20 nucleotide sequences and 1392 positions in the final data set. Cyclotella meneghiniana was used as the outgroup.
Diversity 14 00790 g012
Figure 13. rbcL phylogenetic tree by maximum likelihood method of Nitzschia inclinata sp. nov. and other species molecular position. The tree with the highest log likelihood (−2490.51) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis involved 18 nucleotide sequences and 557 positions in the final data set. Cyclotella meneghiniana was used as the outgroup.
Figure 13. rbcL phylogenetic tree by maximum likelihood method of Nitzschia inclinata sp. nov. and other species molecular position. The tree with the highest log likelihood (−2490.51) is shown. The percentages of trees in which the associated taxa clustered together are shown next to the branches. The analysis involved 18 nucleotide sequences and 557 positions in the final data set. Cyclotella meneghiniana was used as the outgroup.
Diversity 14 00790 g013
Table 1. Primers used for amplification and sequencing of the nuclear SSU rRNA and rbcL.
Table 1. Primers used for amplification and sequencing of the nuclear SSU rRNA and rbcL.
GenePrimerNucleotide Sequence (5′ to 3′)References
SSU
rRNA
AT18F01YAC CTG GTT GAT CCT GCC AGT AG[21]
AT18R02GTTTCAGCCTTGCGACCATACTCC[21]
AT18F02AGA ACG AAA GTT AAG GGA TCG AAG ACG[21]
AT18R01GCTTGATCCTTCTGCAGGTTCACC[21]
EulAAAC CTG GTT GAT CCT GCC AGT[22]
EukBGAT CCT TCT GCA GGT TCA CCT AC[22]
rbcLF3GCT TAC CGT GTA GAT CCA GTT CC[23]
R3CCT TCT AAT TTA CCA ACA ACT G[23]
Table 2. Comparison of morphological characteristics for Hannaea librata nov. and closely related species.
Table 2. Comparison of morphological characteristics for Hannaea librata nov. and closely related species.
Hannaea librataH. pamirensisH. hattorianaH. rectaFragilaria
capucina
Length (μm)44–9942–4540–8529–7128–47
Width (μm)5–5.55.5–7.06–76–73.3–4.2
No. striae (/10 μm)12–1617–1813–1412–1414–17
No. areolae (/10 μm)70–8075–80n.d.n.d.n.d.
ValvesLinearWeakly arcuate to almost linearLanceolate,
slightly arcuate
to almost linear
LanceolateLinear
Valve apicesStrongly rostrateCapitate to rostrate apicesCapitate to rostrate apicesRostrateWeakly rostrate
StriationParallelParallelParallelParallel
Ghost strainCentral areaCentral areaCentral areaCentral arean.d.
Valve faceWavedWavedn.d.WavedAlternate, parallel to
slightly radiate toward the apices
Valve mantleFlatFlatn.d.Wavedn.d.
Valves of central areaExpanded to unilateralExpanded to unilateralExpanded to unilateralExpanded to unilateralFlat
Central areaWide U-shape extending into both striaeUnilaterally tumid on the concave marginUnilaterally tumid on the concave marginHorseshoe shapedRectangular to rhombic
Girdle bandSingle row of small punctuaten.d.n.d.n.d.n.d.
Apical pore fieldsRectangularRectangularn.d.Rectangularn.d.
Spine shapeConicalBifurcated thorn (center)/
conical (near the pole)
n.d.Bifurcated thorn (center)/
conical (near the pole)
Conical near the apex to
triangular in the middle
RimoportulaOne per valveOne per valven.d.One per valveTwo per valve
ReferencesThis study[27][27,28][29,30][31,32,33,34,35]
Table 3. Similarity scores and genetic distances among 15 aligned sequences (1478 bp) based on SSU rRNA gene.
Table 3. Similarity scores and genetic distances among 15 aligned sequences (1478 bp) based on SSU rRNA gene.
SpeciesAccession123456789101112131415
Similarity
1 Hannaea librataON040635 1.0001.0001.0001.0000.9890.9890.9890.9870.9870.9850.9740.9740.9690.969
2 Hannaea librataON0406360.000 1.0001.0001.0000.9890.9890.9890.9870.9870.9850.9740.9740.9690.969
3 Hannaea librataON0406400.0000.000 1.0001.0000.9890.9890.9890.9870.9870.9850.9740.9740.9690.969
4 Hannaea librataON0406410.0000.0000.000 1.0000.9890.9890.9890.9870.9870.9850.9740.9740.9690.969
5 Hannaea librataON0406420.0000.0000.0000.000 0.9890.9890.9890.9870.9870.9850.9740.9740.9690.969
6 Fragilaria capucinaMH3562570.0110.0110.0110.0110.011 0.9991.0000.9960.9960.9940.9790.9790.9650.965
7 F. capucina var. mesoleptaMH9978450.0110.0110.0110.0110.0110.001 0.9990.9960.9960.9920.9790.9790.9660.966
8 Fragilaria bidensAB4305990.0110.0110.0110.0110.0110.0000.001 0.9960.9960.9940.9790.9790.9650.965
9 Fragilaria crotonensisAM7126160.0130.0130.0130.0130.0130.0040.0040.004 1.0000.9970.9770.9770.9650.965
10 Fragilaria vaucheriaeAM4977330.0130.0130.0130.0130.0130.0040.0040.0040.000 0.9970.9770.9770.9650.965
11 Fragilaria vaucheriaeEU2604690.0150.0150.0150.0150.0150.0060.0070.0060.0030.003 0.9740.9740.9630.963
12 Synedra minusculaEF4234150.0270.0270.0270.0270.0270.0220.0220.0220.0240.0240.026 1.0000.9590.959
13 Synedra hyperboreaAY4854640.0270.0270.0270.0270.0270.0220.0220.0220.0240.0240.0260.001 0.9590.959
14 Ulnaria acusKF9596590.0300.0300.0300.0300.0300.0340.0330.0340.0340.0340.0360.0410.041 0.999
15 Ulnaria ulnaMG6843610.0300.0300.0300.0300.0300.0340.0330.0340.0340.0340.0360.0410.0410.001
p-distance
Accession; Gene Bank accession number.
Table 4. Similarity scores and genetic distances among 17 aligned sequences (596 bp) based on rbcL gene.
Table 4. Similarity scores and genetic distances among 17 aligned sequences (596 bp) based on rbcL gene.
SpeciesAccession1234567891011121314151617
Similarity
1 Hannaea librataOP137173 1.0001.0001.0001.0000.9880.9800.9860.9820.9840.9840.9820.9140.9140.9470.9450.947
2 Hannaea librataOP1832090.000 1.0001.0001.0000.9880.9800.9860.9820.9840.9840.9820.9140.9140.9470.9450.947
3 Hannaea librataOP1832120.0000.000 1.0001.0000.9880.9800.9860.9820.9840.9840.9820.9140.9140.9470.9450.947
4 Hannaea librataOP1832130.0000.0000.000 1.0000.9880.9800.9860.9820.9840.9840.9820.9140.9140.9470.9450.947
5 Hannaea librataOP1832140.0000.0000.0000.000 0.9880.9800.9860.9820.9840.9840.9820.9140.9140.9470.9450.947
6 Fragilaria capucinaKT0729280.0120.0120.0120.0120.012 0.9860.9980.9940.9960.9900.9880.9210.9210.9450.9430.945
7 Fragilaria capucinaKC7365940.0200.0200.0200.0200.0200.014 0.9880.9840.9860.9880.9860.9210.9210.9450.9430.945
8 Fragilaria bidensAB4306760.0140.0140.0140.0140.0140.0020.012 0.9920.9940.9920.9900.9230.9230.9470.9450.947
9 Fragilaria crotonensisKF9596400.0180.0180.0180.0180.0180.0060.0160.008 0.9980.9880.9860.9170.9170.9410.9380.941
10 Fragilaria crotonensisKT0729030.0160.0160.0160.0160.0160.0040.0140.0060.002 0.9900.9880.9190.9190.9430.9410.943
11 Fragilaria vaucheriaeHYUD0100.0160.0160.0160.0160.0160.0100.0120.0080.0120.010 0.9940.9230.9230.9470.9450.947
12 Fragilaria perminutaKF9596500.0180.0180.0180.0180.0180.0120.0140.0100.0140.0120.006 0.9210.9210.9450.9430.945
13 Synedra minusculaJN1628250.0810.0810.0810.0810.0810.0750.0750.0730.0790.0770.0730.075 1.0000.9340.9360.934
14 Synedra hyperboreaHQ9124850.0810.0810.0810.0810.0810.0750.0750.0730.0790.0770.0730.0750.000 0.9340.9360.934
15 Ulnaria acusKF9596450.0510.0510.0510.0510.0510.0530.0530.0510.0570.0550.0510.0530.0630.063 0.9981.000
16 Ulnaria ulnaMG6843320.0530.0530.0530.0530.0530.0550.0550.0530.0590.0570.0530.0550.0610.0610.002 0.998
17 Ulnaria ulnaKT0729420.0510.0510.0510.0510.0510.0530.0530.0510.0570.0550.0510.0530.0630.0630.0000.002
p-distance
Accession: GenBank accession number.
Table 5. Comparison of morphological characteristics for Gomphonema seminulum sp. nov. and closely related species.
Table 5. Comparison of morphological characteristics for Gomphonema seminulum sp. nov. and closely related species.
Gomphonema seminulumG. bourbonenseG. pumilumG. angustumG. clevei
Length (μm)15–209.4–2817–3713–13011.5–34
Width (μm)4–53.3–4.75–83–123.5–7
No. striae
(/10 μm)
13–1510.5–1311–1411–1510–15
No. areolae
(/1 μm)
4–5n.d.n.d.n.d.3–4
StriaeParallel to slightly radiate near the head pole,
strongly radiate toward the foot pole
Slightly radiate, almost parallelSlight radiate to strongly radiate toward the apicesSlightly radiate,
parallel, low density of central area
Slightly radiate
AreolaBroadly C-shaped, occluded by raised flaps, 2–4 lines from the axis are flap open toward the center, and the other lines are open in the opposite directionOccluded by raised flaps,
C-shaped, slit-like
C-shaped, 1–3 lines from the axis are flap open toward the center, and the other lines are open in the opposite directionn.d.C-shaped,
occluded by raised flaps, 2–3 lines from the axis are flap open toward the center, and the other lines are open in the opposite direction
Axial areaWide lanceolateAnguste lanceolateSmallWide rectangleWide lanceolate
StigmaOne per valve, circularOne per valve, circularOne per valveOne per valveOne per valve, circular with a thickened margin
ValvesNarrow lanceolate,
slightly asymmetric about the axis
Linear and slightly ellipticalNarrow elliptic lanceolate,
tapering more strongly
toward the head pole
Narrow lanceolate–ovataRhombic clavate,
elliptical clavate
Valve apicesBluntly roundedObtuse, wide circularbluntly rounded,
slightly protracted apex
Wide circularObtuse
(head pole is wider than foot pole)
RapheExtends along the entire valvar face,
slightly undulate
ParallelExtends along the entire valvar faceExtends along the entire valvar faceUndulate (external),
relatively linear (internal)
Central raphe endingElliptical (external), bent to stigma, hooks in opposite directions (internal)Elliptical(external), bent to stigma, hooks in opposite directions (internal)Elliptical (external),
sharp hook shape (internal)
Hook shapeElliptical (external),
bending to stigma,
hook shape (internal)
Polar raphe endingdeflects in the opposite
direction of the
stigma (external),
formed helictoglossa (internal)
n.d.Deflects in the opposite
direction of the
stigma (external),
formed helictoglossa (internal)
n.d.Curved toward the
stigma (external), formed
small helictoglossa (internal)
ReferenceThis study[36][36,37,38][37,38,39][30,40]
Table 6. Similarity scores and genetic distances among 15 aligned sequences (1475 bp) of Gomphonema species based on SSU rRNA gene.
Table 6. Similarity scores and genetic distances among 15 aligned sequences (1475 bp) of Gomphonema species based on SSU rRNA gene.
SpeciesAccession123456789101112131415
Similarity
1 Gomphonema seminulumON040637 0.9920.9920.9920.9910.9920.9920.9920.9910.9910.9760.9760.9740.9740.969
2 Gomphonema bourbonenseKC7366210.007 1.0001.0000.9990.9900.9920.9920.9910.9920.9740.9740.9720.9720.968
3 Gomphonema bourbonenseKC7366200.0070.000 1.0000.9990.9900.9920.9920.9910.9920.9740.9740.9720.9720.968
4 Gomphonema bourbonenseJN7902810.0070.0000.000 0.9990.9900.9920.9920.9910.9920.9740.9740.9720.9720.968
5 Gomphonema bourbonenseKT0729570.0090.0010.0010.001 0.9890.9900.9900.9900.9910.9720.9720.9720.9720.967
6 Gomphonema pumilumKC7366290.0080.0100.0100.0100.012 0.9920.9920.9920.9920.9730.9730.9720.9720.967
7 Gomphonema acuminatumJN7902800.0080.0080.0080.0080.0090.008 1.0000.9990.9990.9780.9780.9770.9770.969
8 Gomphonema acuminatumAM5020190.0080.0080.0080.0080.0090.0080.000 0.9990.9990.9780.9780.9770.9770.969
9 Gomphonema subclavatumKJ0116660.0090.0090.0090.0090.0100.0090.0010.001 0.9990.9770.9770.9760.9760.969
10 Gomphonema truncatumAM5019560.0090.0070.0070.0070.0090.0090.0010.0010.001 0.9770.9770.9760.9760.970
11 Gomphonema parvulumKC7366250.0240.0260.0260.0260.0270.0270.0220.0220.0220.022 1.0000.9970.9980.985
12 Gomphonema parvulumKT0729610.0240.0260.0260.0260.0270.0270.0220.0220.0220.0220.000 0.9970.9980.985
13 Gomphonema dichotomumKJ0116550.0250.0270.0270.0270.0280.0280.0230.0230.0240.0240.0030.003 0.9980.985
14 Gomphonema affineMH9978430.0250.0270.0270.0270.0280.0280.0230.0230.0240.0240.0020.0020.002 0.987
15 Gomphonema cleveiKT0729590.0310.0310.0310.0310.0320.0330.0300.0300.0310.0290.0150.0150.0150.013
p-distance
Accession; Gene Bank accession number.
Table 7. Similarity scores and genetic distances among 12 aligned sequences (596 bp) of Gomphonema species based on rbcL gene.
Table 7. Similarity scores and genetic distances among 12 aligned sequences (596 bp) of Gomphonema species based on rbcL gene.
SpeciesAccession123456789101112
Similarity
1 Gomphonema seminulumOP183210 0.9550.9490.9620.9620.9600.9530.9420.9420.9480.9480.927
2 Gomphonema bourbonenseKC7365950.044 0.9470.9620.9580.9530.9570.9480.9480.9420.9420.926
3 Gomphonema pumilumKC7365990.0490.050 0.9570.9570.9570.9570.9460.9460.9440.9440.933
4 Gomphonema acuminatumLR7425500.0370.0420.037 0.9970.9900.9830.9530.9530.9530.9530.933
5 Gomphonema acuminatumLR7425530.0370.0420.0400.003 0.9930.9860.9530.9530.9530.9530.933
6 Gomphonema subclavatumLR7425560.0390.0420.0450.0100.007 0.9860.9500.9500.9500.9500.929
7 Gomphonema truncatumLR7425470.0450.0420.0420.0170.0130.013 0.9500.9500.9490.9490.933
8 Gomphonema parvulumKY6615710.0550.0520.0500.0450.0450.0490.049 1.0000.9850.9850.964
9 Gomphonema parvulumKY6615540.0550.0520.0500.0450.0450.0490.0490.000 0.9850.9850.964
10 Gomphonema affineAM7104690.0500.0540.0550.0450.0450.0490.0490.0150.015 1.0000.964
11 Gomphonema affineMK5760420.0500.0540.0550.0450.0450.0490.0490.0150.0150.000 0.964
12 Gomphonema cleveiJQ3546820.0690.0640.0700.0640.0640.0670.0640.0350.0350.0350.035
p-distance
Accession; Gene Bank accession number.
Table 8. Comparison of morphological characteristics for Nitzschia inclinata sp. nov. and closely related species.
Table 8. Comparison of morphological characteristics for Nitzschia inclinata sp. nov. and closely related species.
Nitzschia inclinataN. dissipataN. oligotraphentaN. sigmoideaN. angularis
Length (μm)32–34.512.5–8530–45346–35960–200
Width (μm)3.5–53.5–73–49–136–15
No. striae (/10 μm)5039–5046–4824–2631–32
No. areolae (/10 μm)60–8069n.d.40–50n.d.
No. fibulae (/10 μm)7–105–11.58.5–11.56–92.5–5
FibulaeOne-sided valve, irregularly distributedOne side of valve,
irregularly distributed
One-sided valveOne side of valve margin, irregularly distributedCenter of valve longitudinally striated
StriationParallelParallelParallelParallelParallel
HelictoglossaVery apex of the valve
(raphe ends internally formed)
n.d.n.d.Very apex of the valve
(raphe ends internally formed)
Very apex of the valve
ValvesNarrow linearNarrow lanceolateLinear lanceolateLinearRhomboidal
Valve apicesCapitate, bendedsubrostrateDistinctly capitulumTaperingn.d.
RapheEccentricSlightly eccentricEccentricHighly eccentric, one side of valve marginAlmost central of valve face
Terminal fissureBroad sweeping curve over the valve apexHook shape, sometimes widened or bifurcaten.d.Hook-shapedHook-shaped
Girdle viewRectangularn.d.n.d.Sigmoidn.d.
Girdle bandDouble rows of small punctuateDouble rows of small punctuaten.d.Double rows of small punctuaten.d.
ReferenceThis study[41,45,46,47][42,46,47][43][44]
Table 9. Similarity scores and genetic distances among 14 aligned sequences (1392 bp) of Nitzschia species based on SSU rRNA gene.
Table 9. Similarity scores and genetic distances among 14 aligned sequences (1392 bp) of Nitzschia species based on SSU rRNA gene.
SpeciesAccession1234567891011121314
Similarity
1 Nitzschia inclinataON040638 0.9840.9810.9780.9760.9760.9750.9750.9740.9750.9750.9630.9530.937
2 Nitzschia dissipataKY3203940.016 0.9880.9770.9760.9760.9740.9740.9760.9760.9740.9640.9580.943
3 Nitzschia sigmoideaMN7504230.0190.012 0.9750.9720.9720.9740.9740.9710.9720.9700.9590.9530.938
4 Nitzschia pusillaKY3203900.0220.0220.025 0.9930.9930.9940.9940.9890.9890.9930.9810.9600.948
5 Nitzschia paleaMH3584610.0240.0240.0270.006 1.0000.9950.9950.9910.9900.9950.9820.9620.951
6 Nitzschia paleaKY8634740.0240.0240.0270.0060.000 0.9950.9950.9910.9900.9950.9820.9620.951
7 Nitzschia laevisKF1777750.0250.0250.0260.0060.0050.005 1.0000.9900.9910.9930.9820.9600.948
8 Nitzschia thermalisAY4854580.0250.0250.0260.0060.0050.0050.000 0.9900.9910.9930.9820.9600.948
9 Nitzschia epithemoidesFR8655010.0250.0240.0280.0110.0090.0090.0100.010 0.9930.9900.9790.9580.949
10 Nitzschia paleaeformisKY3203830.0250.0240.0270.0110.0100.0100.0090.0090.007 0.9880.9790.9610.948
11 Nitzschia capitellataKT0729780.0250.0250.0290.0060.0050.0050.0060.0060.0100.012 0.9800.9610.948
12 Nitzschia frustulumAJ5351640.0380.0370.0410.0200.0190.0190.0190.0190.0220.0220.021 0.9490.936
13 Nitzschia hantaniiMK5678930.0450.0410.0450.0390.0370.0370.0390.0390.0400.0380.0380.051 0.968
14 Nitzschia inconspicuaMN7504630.0610.0550.0590.0510.0480.0480.0510.0510.0490.0500.0500.0630.031
p-distance
Accession: GenBank accession number.
Table 10. Similarity scores and genetic distances among 12 aligned sequences (557 bp) of Nitzschia species based on rbcL gene.
Table 10. Similarity scores and genetic distances among 12 aligned sequences (557 bp) of Nitzschia species based on rbcL gene.
SpeciesAccession123456789101112
Similarity
1 Nitzschia inclinataOP183211 0.9560.9730.9670.8900.8860.8820.8820.9000.8860.8840.873
2 Nitzschia dissipataKY3203320.043 0.9480.9560.8670.8710.8710.8710.9000.8730.8710.860
3 Nitzschia sigmoideaFN5570330.0270.050 0.9650.8920.8820.8770.8770.9020.8860.8790.875
4 Nitzschia cf. sigmoideaKM9991130.0320.0430.034 0.8900.8900.8810.8810.9120.8960.8880.871
5 Nitzschia capitellataKT0729240.1020.1220.1010.102 0.9760.9650.9650.9120.9260.9280.916
6 Nitzschia paleaKC7366090.1060.1180.1100.1020.023 0.9890.9890.9140.9120.9250.914
7 Nitzschia paleaKF9596390.1100.1180.1130.1100.0340.011 1.0000.9060.9040.9170.906
8 Nitzschia paleaKJ5424640.1100.1180.1130.1100.0340.0110.000 0.9060.9040.9170.906
9 Nitzschia paleaeformisKY3203220.0930.0930.0920.0830.0830.0810.0880.088 0.9300.9400.928
10 Nitzschia hantaniiMK5760400.1060.1170.1060.0970.0700.0830.0900.0900.066 0.9610.950
11 Nitzschia inconspicuaKC7366070.1080.1180.1110.1040.0680.0720.0790.0790.0570.038 0.969
12 Nitzschia frustulumHF6750700.1170.1270.1150.1180.0790.0810.0880.0880.0680.0480.031
p-distance
Accession; Gene Bank accession number.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hwang, E.-A.; Kim, H.-K.; Cho, I.-H.; Yi, C.; Kim, B.-H. Morphological and Molecular Studies of Three New Diatom Species from Mountain Streams in South Korea. Diversity 2022, 14, 790. https://doi.org/10.3390/d14100790

AMA Style

Hwang E-A, Kim H-K, Cho I-H, Yi C, Kim B-H. Morphological and Molecular Studies of Three New Diatom Species from Mountain Streams in South Korea. Diversity. 2022; 14(10):790. https://doi.org/10.3390/d14100790

Chicago/Turabian Style

Hwang, Eun-A, Ha-Kyung Kim, In-Hwan Cho, Chen Yi, and Baik-Ho Kim. 2022. "Morphological and Molecular Studies of Three New Diatom Species from Mountain Streams in South Korea" Diversity 14, no. 10: 790. https://doi.org/10.3390/d14100790

APA Style

Hwang, E. -A., Kim, H. -K., Cho, I. -H., Yi, C., & Kim, B. -H. (2022). Morphological and Molecular Studies of Three New Diatom Species from Mountain Streams in South Korea. Diversity, 14(10), 790. https://doi.org/10.3390/d14100790

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop