Front. Microbiol.Frontiers in MicrobiologyFront. Microbiol.1664-302XFrontiers Media S.A.10.3389/fmicb.2025.1517865MicrobiologyOriginal ResearchPhylogeny and fatty acid profiles of new Coccomyxa (Chlorophyta) species from soils of VietnamMaltsevYevhen1*†KezlyaElena1†MaltsevaSvetlana1†KrivovaZinaida1†ÐinhCù Nguyên2†KulikovskiyMaxim1†1K.A. Timiryazev Institute of Plant Physiology RAS, IPP RAS, Moscow, Russia2Southern Branch of Joint Vietnam-Russian Tropical Science and Technology Research Center, Ho Chi Minh City, Vietnam
Edited by: Jiangxin Wang, Shenzhen University, China
Reviewed by: Zixi Chen, Shenzhen University, China
Zheng Yuan, China Academy of Chinese Medical Sciences, China
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Introduction
A new algal species of Coccomyxa, C. cattiensis sp. nov., C. fusiformis sp. nov., and C. tropica sp. nov. are recorded in the tropical forest soil from the Cát Tiên National Park, Vietnam.
Methods
The analysis is based on morphological characters, evolutionary distance, 18S rDNA phylogeny, and ITS2 secondary structure.
Results
New species differed from other species of the genus by the size and shape of vegetative cells, and habitat type. The evolutionary distance matrix based on the 18S rRNA gene shared 97.9–100% similarities with other Coccomyxa sequences. The phylogeny inferred by maximum likelihood and Bayesian inference placed new species in the independent lineages close to clades C. subellipsoidea, C. polymorpha, and C. parasitica. Predicted secondary structures of the ITS2 of new species differed from the Coccomyxa species by compensatory base changes, deletions, and single bases. Fatty acid composition analysis was performed on Coccomyxa strains from tropical habitats for the first time. α-linolenic (26.8–47.0%), palmitic (13.6–31.2%), and linoleic (up to 21.7%) acids were the dominant fatty acids in the algae cultured on the BBM medium.
Conclusion
High concentrations of omega-3 polyunsaturated fatty acids (up to 485 mg L–1) position the novel Coccomyxa strains as a promising feedstock for the food and pharmaceutical, agricultural, and aquaculture industries.
Cát Tiên National Parkgreen algaeITS secondary structurenew speciesphylogeny18S rDNA23-74-10081124052200012-7Russian Science Foundation10.13039/501100006769Ministry of Science and Higher Education of the Russian Federation10.13039/501100012190section-at-acceptanceTerrestrial Microbiology1 Introduction
The diversity of soil microalgae remains poorly understood despite the widespread use of molecular methods in their study (Lawley et al., 2004; Bates et al., 2013; Hodač, 2016). One of the most extensive studies on the molecular diversity of soil green algae by culture was based on Hodač (2016) in Germany. This work selected 188 monoclonal cultures of green microalgae from 57 sample areas covering deciduous and coniferous forests, meadows, and pastures. Based on an analysis of 18S rDNA, the author distinguished 73 different phylotypes of soil green algae. The species belong mainly to the green algae: Chlorophyta (59 spp.) and Streptophyta (2 spp.), also to Stramenopiles: Xanthophyceae (11 spp.) and Eustigmatophyceae (1 sp.). We compared results obtained with a culture-dependent approach (morphological observations using light microscopy of mixed or monoclonal algal cultures) with those obtained with a culture-independent approach, NGS method based on the 18S rDNA marker (Kezlya et al., 2023). Almost 90% of cultivated species were shown to correspond to clones obtained from environmental samples, with an 18S rDNA similarity threshold of ≤ 0.01. The high proportion of matches confirms the relevance of monoclonal algal cultures studying for replenishing algal collections and amplifying other markers besides 18S rDNA. Seven new species have also been described, and it has been suggested that algae diversity in Central Europe’s soils cannot be fully explored without cultivation (Hodač, 2016).
Tropics are unique regions with high richness of eukaryotic algae and cyanobacteria. The studies of the last two decades have mainly focused on the composition and dynamics of species diversity in cultivated soils or rice fields in India (Vijayan and Ray, 2015), Egypt (Hameed, 2006), and Laos (Fujita and Ohtsuka, 2005). However, these works were carried out without using molecular approaches. The most common direction in tropical areas is the study of aerophyte algae using molecular methods. Many new genera and species of green algae have been described precisely from subaerial habitats: Kalinella (Neustupa et al., 2009) from the bark of Gigantochloa sp. in a park in Singapore, Hylodesmus from decaying bare wood in a tropical forest in Singapore (Eliaš et al., 2010), Xylochloris from tropical trees in Singapore (Neustupa et al., 2011), Elliptochloris bilobata var. corticola from Cibodas Botanical Garden West Java in Indonesia (Eliáš et al., 2008), Symbiochloris tropica from the bark sample in the secondary lowland rain forest in Malaysia (Škaloud et al., 2016). Data on the phylogeny and distribution of green aerophytic algae in the tropical rainforest in Yunnan, China, were also renewed, and a new order Watanabenales was proposed (Li et al., 2021). New species from the genera Calidiella, Jaagichlorella, Massjukichlorella, Watanabea were described by Li et al. (2020; 2021).
The genus Coccomyxa (Trebouxiophyceae, Chlorophyta) was first described by Schmidle (1901). The genus included green coccoid small size microalgae with a single parietal chloroplast lacking pyrenoids and the absence of flagellate stages and with mucilaginous sheath (Cao et al., 2018a,2018b). Several algal strains morphologically similar to Coccomyxa have been redefined and even described as new species. For instance, Choricystis sp. CAUP H 5107, Choricystis sp. GSE4G, Monodus sp. CR2-4, Paradoxia multiseta UTEX LB 2460, and Pseudococcomyxa simplex CCAP 812/5 were transferred to Coccomyxa based on their phylogeny (Darienko et al., 2015; Malavasi et al., 2016). The first revision of the genus using molecular methods was carried out by Darienko et al. (2015). As a result of a phylogenetic analysis of 41 strains, it was proposed to distinguish seven species. Later, Malavasi et al. (2016) performed a similar study by adding 18 more 18S rDNA + ITS sequences with the strains’ ecological features (life forms and habitats). Coccomyxa is widely distributed in various niches: marine, freshwater, and soil, free-living strains with epiphytic, endophytic, parasitic, and symbiotic lifestyles (Darienko et al., 2015; Malavasi et al., 2016; Gustavs et al., 2017; Cao et al., 2018b; Sciuto et al., 2019). Some species are found in extreme habitats: C. onubensis and C. silvae-gabretae in acidic waters (Garbayo et al., 2012; Fuentes et al., 2016; Barcytë and Nedbalová, 2017), C. melkonianii in the river with high content of heavy metals (Malavasi et al., 2016). The following species have been recorded in the tropical region: C. actinabiotis, strain KN-2011-T4 on the Bambussa sp., C. polymorpha, strain KN-2011-T2 on the bark of Ehretia javanica, C. solorinae, strain KN-2011-T5 on the tree bark in Singapore, strain Coccomyxa sp. KN-2011-T1 on the tree bark in Indonesia, strain Coccomyxa sp. KN-2011-T3 on the bark of Syzygium nervosum (Malavasi et al., 2016). Coccomyxa was found in Cuba as a symbiont of lichen-forming fungi of the genus Sticta (Ascomycota, Peltigeraceae) (Lindgren et al., 2020).
Coccoid green algae are one of the promising sources of metabolites that can be utilized for bioremediation and biomass production, with subsequent processing into biofuel, fertilizer, or feed (Lamminen et al., 2019; Ma et al., 2019; Maltsev and Maltseva, 2021). A group with a C16–C18 chain length usually dominates green algae’s fatty acid profiles. The most common fatty acids (FAs) are 16:0 palmitic, 16:3n-3 roughanic, 18:1n-9 oleic, 18:3n-3 α-linolenic acids (Lang et al., 2011; Maltsev and Maltseva, 2021). Some green algae can produce large amounts of unsaturated fatty acids (UFA), which is more than 97.0% of the total fatty acid content (Øezanka et al., 2008; Lang et al., 2011). The maximum values are set for Chlamydomonas, Chloromonas and Choricystis strains. Only in some Chlamydomonas and Didymogenes strains is the omega-3 fatty acid content higher than 80% (Maltsev and Maltseva, 2021). At the same time, Coccomyxa strains represent a source of different fatty acid groups, which can reach high concentrations depending on culture conditions. For example, saturated fatty acids (SFAs) marked in quantity up to 45.50% of total FA content (Soru et al., 2019), monounsaturated fatty acids (MUFAs) with 60.09% (Maltsev et al., 2019), polyunsaturated fatty acids (PUFAs) with 75.30% or omega-3 FAs with 56.39% (Lang et al., 2011). The high content of PUFAs and omega-3 FAs, along with the ability to rapidly increase biomass in photobioreactors (Abe et al., 2014; Soru et al., 2019), enhances the biotechnological potential of Coccomyxa strains for agricultural feed production.
The authors first started the study of soil algae in Cát Tiên National Park using an integrative approach in 2019. Soil samples were collected from flooded and non-flooded forest areas, as disturbed and undisturbed fields. Previous papers present new species of the different groups of algae: diatoms (Kezlya et al., 2020; Kezlya et al., 2021; Glushchenko et al., 2022; Kezlya et al., 2022), as well as cryptophytes (Martynenko et al., 2022) and cyanobacteria (Maltseva et al., 2022). During a study of the tropical soils of the Cát Tiên National Park, we isolated five strains with a Coccomyxa-like morphology but with several unique morphological and phylogenetic characteristics. Thus, we describe three new species from the genus Coccomyxa based on morphology, 18S rDNA phylogeny, ITS2 secondary structures and fatty acid profiles.
2 Materials and methods2.1 Isolation and cultivation
Tropical forest soil samples (Figure 1) were collected in the Cát Tiên National Park, Vietnam, which has a monsoon tropical climate with summer rains. Table 1 presents a list of all samples, along with their geographic locations and measured ecological parameters.
Sampling sites in Cát Tiên National Park, Đồng Nai Province, Vietnam. (A) Forest (KT55); (B) forest (KT30); (C) dry swamp (KT61); (D) forest (KT23).
List of sampling locations with descriptive information.
Sample name
Geographic coordinates, local name
Habitat type
pH
Humidity, %
Relief and vegetation
KT23*
11°26.906′ N, 107°26.490′ E “Dipterocarpus”
Forest
3.8
41.1
Top of the shale ridge. The soil is red-yellow tropical shallow underdeveloped clayey on metamorphic shale. Vegetation: Dipterocarpus turbinatus, Shorea roxburghii, Swintonia floribunda, herbaceous vegetation with a protective cover of 50%, background: Taenitis blechnoides. The site is sufficiently lightened. The Dipterocarpaceae leaves are dominate in litterfall.
KT30*
11°25.725′ N,72107°25.650′ E “Lagerstroemia”
Forest
5.5
39.24
The drive-side slope was very gently sloping and flooded during the wet season. Clayey slates. Lagerstroemia calyculata dominates the forest stands with an admixture of Tetrameles nudiflora, no grass. Litter is composed mainly of the bark and leaves of L. calyculata. The dark clayey tropical soils developed from the derivatives of basaltic rocks.
KT55
11°23.917′ N, 107°22.523′ E
Forest
5.55
44.8
Level surface with slightly pronounced microtopography/restored forest, lightened area, no flooding during the wet season. There is no exact description of vegetation and soil features.
KT61
11°30.999′ N, 107s°23.269′ E
Dry swamp
5.11
49.02
Grassy wetland. Usually, it floods during the wet season and dries up during the dry season. There was no water at the time of sampling. There is no exact description of vegetation and soil features.
*Morphological and chemical properties, bulk elemental composition and forms of iron compounds in the Cát Tiên National Park soil plot presented in Khokhlova et al. (2017).
A brief description of vegetation and soil type for samples KT23 and KT30 is given using data from Anichkin (2011) and Khokhlova et al. (2017). The authors describe the samples KT55 and KT61. Immediately after sampling, the absolute humidity was determined by the “hot drying” method (Vadjunina and Korchagina, 1986) in the laboratory room, then brought to the air-dry state and packaged. To measure pH, 30 g of soil was weighed, to which 150 mL of distilled water was added (Arinushkina, 1970). The suspension was poured into a clean glass and measured using the Hanna HI 98112 Piccolo 2 (Hanna Instruments, Inc., Woonsocket, RI, United States).
The novel strains VP336, VP339, VP449, VP451, and VP521 were isolated by micropipetting from algae enrichment cultures using an inverted Zeiss Axio Vert A1 microscope (Zeiss, Germany). Small amounts of soil samples were taken to obtain enrichment cultures, placed in Petri dishes, and moisturized. The strains were deposited in the Culture and Barcode Collection of Microalgae and Cyanobacteria “Algabank” (WDCM 1318) at K.A. Timiryazev Institute of Plant Physiology RAS (Moscow, Russia) as perpetually-transferred pure cultures.
Light microscopic observations were performed with a Zeiss Axio Scope A1 (Carl Zeiss Microscopy GmbH, Göttingen, Germany) microscope equipped with an oil immersion objective (Plan-apochromatic × 100/n.a.1.4, Nomarski differential interference contrast, DIC) and a Zeiss Axio Cam ERc 5s camera (Carl Zeiss NTS Ltd., Oberkochen, Germany). Cells were stained with 0.1% (w/v) methylene blue solution and 1.0% (w/v) ink solution to determine the mucilage structure. Observation of the strain lasted from 24 h to 6 months. The culture was maintained on the BBM medium (Bischoff and Bold, 1963). Starvation experiments using a 3N-BBM + V medium, as described by Darienko et al. (2015), were conducted to induce mucilage formation.
For biochemical analysis, the cultures were maintained in 250 mL Erlenmeyer glass flasks containing 150 mL of standard BBM medium under constant orbital shaking (150 rpm in ELMI Sky Line Shaker S-3 L, ELMI Ltd., Riga, Latvia) for 15 days. The cultures grew at standardized conditions with 24°C under light intensity 70 μmol photons m–2 s–1 with color temperature 4,000 K and a 16:8 h light/dark photoperiod. The light intensity and color temperature were measured using the Sekonic C-800 spectrometer (Sekonic Corporation, Tokyo, Japan). During strain growth, biomass increase was monitored by optical density changes using IMPLEN Nanophotometer P300 (Implen GmbH, Germany) at λ = 720 nm (OD720).
2.2 Molecular analysis
The DNA of the investigated strains VP336, VP339, VP449, VP451, and VP521 were extracted using Chelex 100 Chelating Resin, molecular biology grade (Bio-Rad Laboratories, Hercules, CA, USA), according to the manufacturer’s protocol (Bio-Rad, 2020). Amplification of the 18S rRNA gene (1,741–1,791 bp) was performed with a pair of primers 18S–FA2 (5′–ACC TGG TTG ATC CTG CCA GTA–3′) and 18S–RB2 (5′–GAT CCT TCT GCA GGT TCA CCT ACG–3′) (Nakada et al., 2010). Amplification conditions for the 18S rRNA gene were as follows: initial denaturation for 5 min at 95°C followed by 32 cycles of 30 s denaturation at 94°C, 40 s annealing at 64°C, and 90 s extension at 72°C, with the final extension for 5 min at 72°C. Amplification of the site of ITS1–5.8S rDNA–ITS2 region (576–744 bp) was performed with a pair of primers ITS1 (5′–TCC GTA GGT GAA CCT GCG G–3′) and ITS4 (5′–TCC TCC GCT TAT TGA TAT GC–3′) (White et al., 1990). Amplification conditions for the ITS1–5.8S rDNA–ITS2 region were as follows: initial denaturation for 5 min at 95°C followed by 35 cycles of 30 s denaturation at 94°C, 30 s annealing at 60°C, and 60 s extension at 72°C, with the final extension for 5 min at 72°C.
The PCR products were visualized by horizontal electrophoresis in 1.0% agarose gel stained with SYBR™ Safe (Life Technologies, Carlsbad, CA, United States). The products were purified with a mixture of FastAP, 10 × FastAP Buffer, Exonuclease I (Thermo Fisher Scientific, Waltham, MA, United States), and water. Sequencing was performed using a Genetic Analyzer 3500 instrument (Applied Biosystems, Waltham, MA, United States). We used internal primers provided by Maltsev et al. (2018) for sequencing.
Editing and assembling of the consensus sequences were carried out by processing the direct and reverse chromatograms in Ridom TraceEdit (ver. 1.1.0) (Ridom GmbH, Münster, Germany) and Mega (ver. 7.0.26) software. The 18S rDNA and ITS1–5.8S rDNA–ITS2 sequences of the novel strains were included in the alignments of 75 Coccomyxa sequences from the National Center for Biotechnology Information1 (taxa names and Accession Numbers are given in Figure 2). Elliptochloris bilobata and Hemichloris antarctica were chosen as the outgroups. The 18S rDNA and ITS1–5.8S rRNA–ITS2 sequences of all strains (including outgroups) were aligned in Mega (ver. 7.0.26) software according to their secondary structure by means of comparing the structure presented for the strain NIES 2166 of C. subellipsoidea and SAG 216-3b of C. simplex (Darienko et al., 2015). The resulting alignments had lengths of 2,706 characters.
Phylogenetic position of new Coccomyxa species (indicated in bold) based on Bayesian inference from an alignment of 82 sequences and 2,706 characters (18S rRNA gene and ITS1–5.8S rDNA–ITS2 region). Values above the horizontal lines are bootstrap support from ML analyses (values below 50 are not shown). Values under the horizontal lines (or to the right of the slash) are Bayesian posterior probabilities (values below 0.9 are not shown). Strain numbers (if available) and NCBI database accession numbers are indicated for all sequences. Colored boxes represent clades according to the most recent classification or OTUs described in this work or previous papers. The classification sensu Darienko et al. (2015) marked to the right of the vertical line.
The data set was analyzed using the Bayesian inference (BI) method implemented in Beast ver. 1.10.1 software (BEAST Developers, Auckland, New Zealand) (Drummond and Rambaut, 2007) to construct a phylogeny. The most appropriate substitution model for the alignment partition, shape parameter α and a proportion of invariable sites (pinvar) were estimated using the Bayesian information criterion (BIC) as implemented in jModelTest ver. 2.1.10 (Vigo, Spain) (Darriba et al., 2012). This BIC-based model selection procedure selected the TrN + I + G model, shape parameter α = 0.3770 and a proportion of invariable sites (pinvar) = 0.5500.
We used the HKY model of nucleotide substitution instead of TrN, given that it was the best matching model available for BI. A Yule process tree prior was used as a speciation model. The analysis ran for 5 million generations with chain sampling every 1,000 generations. The parameters-estimated convergence, effective sample size (ESS), and burn-in period were checked using the Tracer ver. 1.7.1 software (MCMC Trace Analysis Tool, Edinburgh, United Kingdom) (Drummond and Rambaut, 2007). The initial 25% of the trees were removed, and the rest were retained to reconstruct a final phylogeny. The phylogenetic tree and posterior probabilities of its branching were obtained based on the remaining trees, having stable estimates of the parameter models of nucleotide substitutions and likelihood. The Maximum Likelihood (ML) analysis was performed using RAxML software (Stamatakis et al., 2008). The non-parametric bootstrap analysis with 1,000 replicas was used. FigTree ver. 1.4.4 (University of Edinburgh, Edinburgh, United Kingdom) and Adobe Photoshop CC ver. 19.0 software (Adobe, San Jose, CA, United States) were used for viewing and editing the trees.
The Mfold version 2.5 software was used to model the secondary structure of ITS2 (Zuker, 2003). The presence of UU pyrimidine–pyrimidine unpaired section in the second hairpin and the conservative GGUAG motive in the 5′ side of helix III (Coleman, 2007), as well as the length and nucleotide composition of the spacers in the central loop determining the helix boundaries, were taken into consideration when constructing the final ITS2 model (Caisová et al., 2013). Also, we used a concept of ITS2 secondary structure provided by Darienko et al. (2015) for Coccomyxa to avoid ambiguity in the interpretation of results. The construction of a hybrid stem with the 5.8S rDNA terminal site and the complementary 28S rDNA start site determined the beginning and the end of ITS2. The PseudoViewer3 program (Byun and Han, 2009) and Adobe Photoshop CC ver. 19.0 software were used to visualize the resulting secondary structure.
The 18S rRNA gene or whole 18S rDNA–ITS1–5.8S rDNA region (if available) was also used to estimate the degree of similarity between gene sequences of different Coccomyxa strains. Using Mega (ver. 7.0.26) software, the p-distances were determined to calculate the sequence similarity with the formula (1–p) × 100.
2.3 Fatty acid composition analysis
Biomass preparation for the fatty acid methyl ester (FAME) profiling was performed according to Maltseva et al. (2022). The VP336, VP339, VP449 and VP521 cultural suspensions with 150 ± 6.0 mg algal wet weight were transferred to 50 mL tubes. The cells were pelleted at room temperature for 3 min at 3,600 g. The supernatants were removed, and the pelleted cells were resuspended in a 10 mL volume of distilled water, quantitatively transferred to 15 mL centrifuge tubes, and pelleted again by centrifugation. The supernatants were removed, and all the quantities of biomass were transferred to a 50-mL round-bottom flask. To avoid the oxidation of unsaturated FAs, all samples were dried with an argon flow. Heptadecanoic acid (Sigma-Aldrich, United States) was used as an internal standard for determining fatty acid composition (Maltseva et al., 2022). Ten milliliter of 1M KOH in 80% aqueous ethanol were added to the dry residue, the flask was sealed with a reflux condenser and kept for 60 min at the mixture’s boiling point (∼80°C). After the time-lapse, the solvents were evaporated in vacuo. The resulting volume of ∼3 mL was quantitatively transferred to a 50-mL centrifuge tube with distilled water added to a total volume of 25 mL. Unsaponifiable components were extracted with three 10 mL changes of n-hexane (Himmed, Russia); the tube was centrifuged at room temperature for 5 min at 2,022 g to accelerate phase separation. After that, the aqueous phase was acidified with a few drops of 20% sulfuric acid (Himmed, Russia) to a slightly acidic pH (checked with indicator paper). Free FAs were extracted with 15 mL of n-hexane. The extraction was repeated three times. The hexane solution of free FAs was transferred to a dry 50-mL round-bottom flask, and the solvent was evaporated to dryness on a rotary evaporator (IKA RV-10, Germany), after which 10 mL of absolute methanol (Sigma-Aldrich, United States) and 1 mL of acetyl chloride (Sigma-Aldrich, United States) were added to the dry residue. The flask, closed with a reflux condenser, was kept for 1 h at 70°C, then the solvents were evaporated to dryness, a few drops of distilled water were added to the dry residue, and FAMEs were extracted with n-hexane.
The obtained FAMEs were analyzed using an Agilent 7890A GC gas chromatograph (Agilent Technologies, United States) equipped with an Agilent 5975Ñ mass spectrometric detector and G2591A inert ion source. A DB-23 capillary column (Ser. no. US8897617H, B&W, United States), 60 m long and 0.25 mm in diameter and containing a grafted (50% cyanopropyl) methylpolysiloxane polar liquid phase as a 0.25-mm-thick film, was used for the analysis. The remaining conditions of the analysis were as follows: pressure in the injector 191 kPa, carrier gas helium at a flow rate of 1 mL min–1, injected sample volume 1 μL, 1:5 flow split ratio, and evaporation temperature of 260°C. The temperature gradient program was as follows: 130–170°C at 6.5°C min–1 steps; 170–215°C at 2.5°C min–1 increments, hold at 215°C for 25 min, 215–240°C at 40°C min–1 increments, and the final hold at 240°C for 50 min; the operating temperature of the mass spectrometric detector was 240°C and ionization energy was 70 eV.
2.4 Data analysis
All analyses were performed in triplicate. Tables show the mean values and standard errors. The significance of differences between the groups was evaluated by one-way analysis of variance (ANOVA) with a Tukey’s post-hoc test performed using Statgraphics Centurion ver. 18 software. A difference between the two groups was declared significant at P < 0.05.
3 Results
The studies performed light microscopy, 18S rDNA phylogeny, ITS2 secondary structures, and fatty acid profiles of five new Coccomyxa strains. Molecular and morphological investigations showed that strains VP339, VP449, VP451, and VP521 are unknown taxa. Therefore, they are described herein as new species: Coccomyxa fusiformis Maltsev et Kezlya, Coccomyxa tropica Maltsev et Kezlya, Coccomyxa cattiensis Maltsev et Kezlya.
3.1 Description of new species and strains3.1.1 <italic>Coccomyxa subellipsoidea</italic> E. Acton3.1.1.1 Strain <italic>Coccomyxa subellipsoidea</italic> VP336
Morphological description (Figure 3): Mature vegetative cells solitary, ellipsoidal, pointed on one side, sometimes slightly curved, asymmetrical, without mucilaginous sheath, cell size 5.0–7.1 μm in length and 1.6–3.9 μm in width. Chloroplast parietal, trough-shaped, thin, straight or with small concavity in the middle, covering about half of the cell volume. Pyrenoid is absent. Cell wall is thin. Reproduction by 2–4 autospores. Protoplast division is oblique. Liberation of autospores through the rupture of mother cell wall at one end. The released autospores have the same morphology as mature vegetative cells.
Nomarski interference micrographs of Coccomyxa subellipsoidea, strain VP336 in culture. Scale bar = 10 μM. (A–F) Variable shapes of culture cells, age 2 weeks; (G,H) Asexual reproduction through cell division with the formation of 2–4 autospores, age 4 weeks.
Culture information: The culture was deposited and maintained as an active culture under the designation VP336 in the Culture and Barcode Collection of Microalgae and Cyanobacteria “Algabank” (WDCM 1318) at K.A. Timiryazev Institute of Plant Physiology RAS.
Habitat: Forest soil, Cát Tiên National Park, Đồng Nai Province, Vietnam (11°23.917′ N, 107°22.523′ E): soil sample KT55 collected on 25 June 2019.
Sequence data: Accession number PV748123 in the National Center for Biotechnology Information (see text footnote 1) for the 18S rRNA gene and ITS1–5.8S rDNA–ITS2 region sequence.
Remarks: In size and form of cells similar to C. onubensis, C. vinatzeri, C. melkonianii, C. fottii and C. cattiensis, C. tropica from this study (Supplementary material 1). The new strain has a thinner chloroplast, covering no more than half of the cell volume, whereas in the above species, a massive chloroplast occupies more than half of the cell volume. Also, C. fottii, C. onubensis, and C. tropica lipid droplets were present in the cytosol, whereas in C. subellipsoidea they were not observed.
ITS2 analysis: Comparative analysis of ITS2 sequences for the strain VP336 vs C. subellipsoidea SAG 216-13 (type culture), SAG 216-7 and “Choricystis sp.” CAUP H5105 strains revealed several evolutionary events in the clade C. subellipsoidea (Figures 2, 4). The VP336 and SAG 216-13, SAG 216-7 and “Choricystis sp.” CAUP H5105 sequences differed by one A→C transversion within the barcode region in the helix III; one compensatory base change (CBC) U–A → C–G and one hemi-compensatory base change (hCBC) U–G → C–G in the helix IV (outside the barcode region). Also, one hCBC in the helix III (G–C → G–U) and one U→C transition in the 5.5S/LSU stem differed VP336 from CAUP H5105.
Predicted secondary structure of the ITS2 sequence of the Coccomyxa subellipsoidea, strain VP336. Base numbering is indicated every 10 bases. The four helices are numbered with Roman numerals. The barcode region of the ITS2 sensu Darienko et al. (2015) is shown by lines inside the secondary structure. hCBC within the barcode region for “Choricystis sp.” CAUP H5105 is shown by a triangle. Single bases within the barcode region for “Choricystis sp.” CAUP H5105 are shown by squares with dotted sides. Single base within the barcode region for C. subellipsoidea SAG 216-13, SAG 216-7 and “Choricystis sp.” CAUP H5105 is shown by a circle. CBC outside the barcode region for C. subellipsoidea SAG 216-13, SAG 216-7 and “Choricystis sp.” CAUP H5105 is shown by rectangle and hCBC is shown by pentagon. Deletions and single bases outside the barcode region are not shown.
3.1.2 <italic>Coccomyxa fusiformis</italic> Maltsev et Kezlya
Diagnosis (Figure 5): Mature vegetative cells solitary, elongated-oval, very often pointed on both sides, spindle-shaped, slightly curved, asymmetric, without mucilaginous sheath, cell size 8.0–15.0 (16.2) μm in length and 1.9–3.0 (3.4) μm in width. Chloroplast parietal, thin, trough-shaped or ribbon-like, with waved margin in large cells, covering about half of the cell volume. Pyrenoid is absent. Lipid droplets present in the cytosol. Cell wall is thin. Reproduction by 2–4 autospores. Protoplast division is oblique. Liberation of autospores through the rupture of mother cell wall at one end. The released autospores have the same morphology as mature vegetative cells.
Nomarski interference micrographs of Coccomyxa fusiformis, strain VP339 in culture. Scale bar = 10 μM. (A–F) Variable shapes of culture cells, age 2 weeks; (G–J) asexual reproduction through cell division with the formation of 2–4 autospores, age 4 weeks.
Holotype (designated here): Material from the authentic strain VP339 is stored in the IPPAS collection of microalgae and cyanobacteria (WDCM 596), Moscow, Russian Federation (metabolically inactive cryopreserved culture). Preserved specimen with unfixed, dried cells (in a metabolically inactive state) from a large drop of a batch of strain VP339 added to watercolor paper (Figure 5), deposited at the MHA (Herbarium, Main Botanical Garden, Botanicheskaya Str. 4, Moscow, 127276, Russia) under the designation Coccomyxa fusiformis Vietnam Kezlya strain VP339. The culture was also deposited and maintained as an active culture under the designation VP339 in the Culture and Barcode Collection of Microalgae and Cyanobacteria “Algabank” (WDCM 1318) at K.A. Timiryazev Institute of Plant Physiology RAS. Also, the strain is conserved as a formaldehyde-fixed sample.
Reference strain: VP339.
Type locality: Forest soil, Cát Tiên National Park, Đồng Nai Province, Vietnam (11°30.999′ N, 107°23.269′ E): soil sample KT61 from dry swamp collected on 25 June 2019.
Sequence data: Accession number PV748126 in the National Center for Biotechnology Information (see text footnote 1) for the 18S rRNA gene and ITS1–5.8S rDNA–ITS2 region sequence.
Etymology: The species epithet “fusiformis” refers to the spindle-like cell shape.
Remarks: The new species differs in cell size and shape from both those studied in this study and the phylogenetically closely C. subellipsoidea (Darienko et al., 2015) as it has relatively large elongated cells (8.0–16.2 μm in length, 1.9–3.4 μm in width, Supplementary material 1). Morphologically similar to C. polymorpha (8.0–13.0 μm in length, 2.0–3.5 μm in width) but cells not gathered in star-like structures. Exact identification is possible only by using 18S rDNA and ITS2 sequence analysis.
Distribution: As yet known only from the type locality.
ITS2 analysis: Several modifications within the barcode region have been shown in the result of the comparison of the secondary structure of ITS2 between C. fusiformis vs C. subellipsoidea SAG 216-13, SAG 216-7 and “Choricystis sp.” CAUP H5105: one hCBC in the helix I (U–G → C–G) and one G→A transition in the helix III (Figure 6). Additionally, one hCBC (G–U → G–C) and one G→A transition in helix III distinguished VP339 from SAG 216-13 and SAG 216-7. One hCBC (G–C → U–C) in the helix IV (outside the barcode region) differed C. fusiformis VP339 from other taxa in the clade C. subellipsoidea (Figure 2). Thus, the most significant differences in nucleotide sequences of ITS2 of the strain VP339 were located inside and outside the barcode region.
Predicted secondary structure of the ITS2 sequence of the Coccomyxa fusiformis, strain VP339. Base numbering is indicated every 10 bases. The four helices are numbered with Roman numerals. The barcode region of the ITS2 sensu Darienko et al. (2015) is shown by lines inside the secondary structure. hCBC within the barcode region for C. subellipsoidea SAG 216-13, SAG 216-7 and “Choricystis sp.” CAUP H5105 is shown by a rectangle. hCBC within the barcode region for C. subellipsoidea SAG 216-13, SAG 216-7 is shown by a triangle. Single base within the barcode region for C. subellipsoidea SAG 216-13, SAG 216-7 and “Choricystis sp.” CAUP H5105 is shown by a circle. A single base within the barcode region for C. subellipsoidea SAG 216-13, and SAG 216-7 is shown by hexagon. Single base within the barcode region for “Choricystis sp.” CAUP H5105 is shown by a square with dotted sides. hCBC outside the barcode region for C. subellipsoidea SAG 216-13, SAG 216-7 and “Choricystis sp.” CAUP H5105 is shown by the pentagon. Deletions and single bases outside the barcode region are not shown.
3.1.3 <italic>Coccomyxa tropica</italic> Maltsev et Kezlya
Diagnosis (Figure 7): Mature vegetative cells solitary, ellipsoidal, sometimes ovoid, without mucilaginous sheath, cell size 3.3–6.9 μm in length and 1.5–2.7 μm in width. Chloroplast parietal, trough-shaped, thick, covering about half of the cell volume. Pyrenoid is absent. Lipid droplets present in cytosol. Cell wall is thin. Reproduction by 2–4 autospores. Protoplast division is oblique, rarely transverse. Liberation of autospores through the rupture of mother cell wall at one end. The released autospores have the same morphology as mature vegetative cells.
Nomarski interference micrographs of Coccomyxa tropica, strain VP521 in culture. Scale bar = 10 μM. (A,B) Variable shapes of culture cells, age 2 weeks; (C,D) asexual reproduction through cell division with the formation of 2–4 autospores, age 4 weeks.
Holotype (designated here): Material from the authentic strain VP521 is stored in the IPPAS collection of microalgae and cyanobacteria (WDCM 596), Moscow, Russian Federation (metabolically inactive cryopreserved culture). Preserved specimen with unfixed, dried cells (in a metabolically inactive state) from a large drop of a batch of strain VP521 added to watercolor paper (Figure 7), deposited at the MHA (Herbarium, Main Botanical Garden, Botanicheskaya Str. 4, Moscow, 127276, Russia) under the designation Coccomyxa tropica Vietnam Kezlya strain VP521. The culture was also deposited and maintained as an active culture under the designation VP521 in the Culture and Barcode Collection of Microalgae and Cyanobacteria “Algabank” (WDCM 1318) at K.A. Timiryazev Institute of Plant Physiology RAS. Also, the strain is conserved as a formaldehyde-fixed sample.
Reference strain: VP521.
Type locality: Forest soil, Cát Tiên National Park, Đồng Nai Province, Vietnam (11°26.906′ N, 107°26.490′ E): soil sample KT23 collected on 05 June 2019.
Sequence data: Accession number PV748127 in the National Center for Biotechnology Information (see text footnote 1) for the 18S rRNA gene and ITS1–5.8S rDNA–ITS2 region sequence.
Etymology: The species epithet ‘tropica’ refers to the climate region in which the strain was isolated.
Remarks: C. tropica is similar in shape and cell size to the most closely phylogenetically related C. onubensis, C. actinabiotis and some other species: C. cimbrica, C. greatwallensis (Supplementary material 1). Exact identification is possible only by using 18S rDNA and ITS2 sequence analysis.
Distribution: Known from the type locality and the bark of Syzygium nervosum in Indonesia (Java) for the strain KN-2011-T3.
ITS2 analysis: We analyzed the ITS2 secondary structure of C. tropica and strains from the clade C. polymorpha: Coccomyxa sp. KN-2011-T3, “C. onubensis” ACCV1 and “C. actinabiotis” CCAP 216/25 (Figures 2, 8). The VP521 and KN-2011-T3 differed only by deletions and single bases outside the barcode region. One hCBC (G–U → G–C) in the helix IV (outside the barcode region) differed C. tropica from “C. onubensis” ACCV1. Comparison of the secondary structure of ITS2 for C. tropica and “C. actinabiotis” CCAP 216/25 showed several distinctive features within the barcode region: two CBCs (U–G → A–U; C–G → U–A) and one hCBC (G–C → G–U) in the helix III; total six single bases in the 5.5S/LSU stem, helices II and III. These differences in ITS2 between the strain VP521 and the CCAP 216/25 can be considered interspecific variations.
Predicted secondary structure of the ITS2 sequence of the Coccomyxa tropica, strain VP521. Base numbering is indicated every 10 bases. The four helices are numbered with Roman numerals. The barcode region of the ITS2 sensu Darienko et al. (2015) is shown by lines inside the secondary structure. Deletions and single bases for Coccomyxa sp. KN-2011-T3 are shown by a square with dotted sides. hCBC outside the barcode region for “C. onubensis” ACCV1 is shown by a triangle. CBCs within the barcode region for “C. actinabiotis” CCAP 216/25 are shown by a rectangle. hCBC within the barcode region for “C. actinabiotis” CCAP 216/25 is shown by a pentagon. Deletions and single bases within the barcode region for “C. actinabiotis” CCAP 216/25 are shown by a circle. Deletions and single bases outside the barcode region for “C. actinabiotis” CCAP 216/25 and “C. onubensis” ACCV1 are not shown.
3.1.4 <italic>Coccomyxa cattiensis</italic> Maltsev et Kezlya
Diagnosis (Figure 9): Mature vegetative cells solitary, ovoid, ellipsoidal to long ellipsoidal, sometimes slightly curved, asymmetrical, pointed on one side, without mucilaginous sheath, cell size 3.6–5.7 (7.1) μm in length 1.5–2.6 μm in width. Chloroplast parietal, trough-shaped, thin, straight or with small concavity in the middle, covering half or 2/3 of the cell volume. Pyrenoid is absent. Cell wall is thin. Reproduction by 2–4 autospores. Protoplast division is transverse or oblique. Liberation of autospores through the rupture of mother cell wall at one end. The released autospores have the same morphology as mature vegetative cells.
Nomarski interference micrographs of Coccomyxa cattiensis, strain VP451 in culture. Scale bar = 10 μM. (A–D) Variable shapes of culture cells, age 2 weeks; (E–H) asexual reproduction through cell division with the formation of 2–4 autospores, age 4 weeks.
Holotype (designated here): Material from the authentic strain VP451 is stored in the IPPAS collection of microalgae and cyanobacteria (WDCM 596), Moscow, Russian Federation (metabolically inactive cryopreserved culture). Preserved specimen with unfixed, dried cells (in a metabolically inactive state) from a large drop of a batch of strain VP451 added to watercolor paper (Figure 9), deposited at the MHA (Herbarium, Main Botanical Garden, Botanicheskaya Str. 4, Moscow, 127276, Russia) under the designation Coccomyxa cattiensis Vietnam Kezlya strain VP451. The culture was also deposited and maintained as an active culture under the designation VP451 in the Culture and Barcode Collection of Microalgae and Cyanobacteria “Algabank” (WDCM 1318) at K.A. Timiryazev Institute of Plant Physiology RAS. Also, the strain is conserved as a formaldehyde-fixed sample.
Reference strain: VP451, paratype strain: VP449.
Type locality: Forest soil, Cát Tiên National Park, Đồng Nai Province, Vietnam (11°23.917′ N, 107°22.523′ E): soil sample KT55 collected on 07 June 2019.
Sequence data: Accession numbers PV748125 for the strain VP451, and PV748124 for the strain VP449 in the National Center for Biotechnology Information (see text footnote 1) for the 18S rRNA gene and ITS1–5.8S rDNA–ITS2 region sequences.
Etymology: The species epithet “cattiensis” refers to the name of Cát Tiên National Park, where this species was observed.
Remarks: C. cattiensis differs in size and shape of cells from a closely related in the phylogeny C. antarctica. The cells of the latter are larger: 8.0–12.0 μm in length, 4.0–7.0 in μm width, ovoid to ellipsoidal in the form vs. 3.6–7.1 μm in length, 1.5–2.6 μm in width in C. cattiensis ovoid to long ellipsoidal in the form. For sister taxon in a phylogenetic lineage, C. parasitica. South showed a very variable in shape from cells: spherical to elliptical, oblong or sickle-shaped, frequently with an attenuated hyaline tip, whereas spherical and sickle-shaped not observed in new species. Several species are very similar in morphology: C. vinatzeri and C. fottii (Supplementary material 1). There are no clear differences in the size and shape of cells or the structure of the chloroplast in these species. Exact identification is possible only by using 18S rDNA and ITS2 sequence analysis.
Distribution: As yet known only from the type locality. Strain VP449 isolated from forest soil, sample KT30, Cát Tiên National Park, Đồng Nai Province, Vietnam (11°25.725′ N, 107°25.650′ E) (Table 1).
ITS2 analysis: The present study detected several distinctive features between C. cattiensis vs C. antarctica Ua6 (FACHB-2140) inside the barcode region: one CBC in the helix I (C–G → U–A) and one CBC in the helix III (U–A → C–G), one hCBC (C–G → U–G) in the helix II, and several deletions and single bases in the 5.5S/LSU stem, helices II and III (Figure 10). Thus, the described species C. cattiensis has high evolutionary isolation from C. antarctica. Strains VP451 and VP449 differed by one deletion outside the barcode region and one U→C transition in helix II.
Predicted secondary structure of the ITS2 sequence of the Coccomyxa cattiensis, strain VP451. Base numbering is indicated every 10 bases. The four helices are numbered with Roman numerals. Differences with the strain VP449 are shown by circle and ellipse. The barcode region of the ITS2 sensu Darienko et al. (2015) is shown by lines inside the secondary structure. Differences in characteristics within the barcode region for Coccomyxa antarctica Ua6 (FACHB-2140) are shown by nucleotides outside the secondary structure. CBCs are shown by rectangles, and hCBC is shown by triangles. Deletion and single bases are shown by squares with dotted sides. Deletions and single bases outside the barcode region are not shown.
3.2 Growth parameters
The novel Coccomyxa strains were used in cultivation experiments to evaluate growth dynamics. The initial cultures had OD720 of 0.08. The growth curves showed typical sigmoid character for all cultures (Figure 11). The lag phase was observed during the initial 3 days of cultivation for all strains. After day 3, all cultures entered the exponential phase with a steady increase in optical density till day 11 for VP339 and VP521 strains or day 13 for VP336 and VP449 strains. The cultures subsequently slowed their growth and entered a stationary phase, during which OD720 and cell counts remained almost constant. Starting from day 7, OD720 was significantly lower for VP339 and VP521 than VP336 and VP449 (Figure 11). The observed character of biomass accumulation by Coccomyxa strains suggests that the strain C. cattiensis VP449 had a higher optical density of culture.
Growth kinetics of Coccomyxa strains displaying changes in an optical density (OD), arithmetic means ± standard errors, n = 3.
3.3 Fatty acid profiles
Analysis of the biomass at the stationary phase revealed the domination of PUFAs for all strains, notably 18:2n-6 linoleic acid (within the range of 4.54–21.69% of total FAs) and 18:3n-3 α-linolenic (26.78–47.02%) (Tables 2, 3; Supplementary material 2). The highest content of PUFAs (75.11% of the total FA content corresponding to the total yield of 655.55 mg L–1) was observed for C. fusiformis VP339. In addition, C. subellipsoidea VP336, C. tropica VP521, and C. cattiensis VP449 accumulated high amounts of saturated 16:0 palmitic (22.01–31.16%) and monounsaturated 16:3n-3 roughanic (7.70–16.12%) acids. Small amounts of monounsaturated 16:1n-9 hypogeic (≤ 0.65%) and 16:1n-7 palmitoleic (≤ 3.18%) acids were determined in all strains, whereas C. subellipsoidea VP336, C. tropica VP521 and C. cattiensis VP449 also produced small amounts of the long-chain 20:4n-6 arachidonic acid (≤ 1.90%). We determined the long-chain omega-3 20:5n-3 eicosapentaenoic acid in small concentrations (0.85%) exclusively in C. subellipsoidea VP336 biomass. In general, novel strains were marked by a high content of omega-3 FAs. The greatest concentration of omega-3 acids (62.64%), including 18:3n-3 α-linolenic acid, was observed in the strain C. subellipsoidea VP336. At the same time, the VP339 culture contained the highest total yield of omega-3 acids, 485.28 mg L–1 (Table 3). New species also had specific features in FA profiles. C. fusiformis VP339 differed from the other species by the highest concentrations of 16:2n-6 hexadecadienoic (F3,8 = 15.83, P = 0.001), 16:4n-3 hexadecatetraenoic (F1,4 = 74.14, P = 0.001) and 18:1n-9 oleic (F3,8 = 15.83, P = 0.001) acids and low content of 16:0 palmitic (F3,8 = 15.83, P = 0.001) and 16:3n-3 roughanic acids (F3,8 = 15.83, P = 0.001) (Table 2). C. cattiensis VP449 had the highest content of 16:0 palmitic, 18:2n-6 linoleic and omega-6 FAs (F3,8 = 15.83, P = 0.001 for all values). C. subellipsoidea VP336 was characterized by maximum value of omega-3 FAs (F3,8 = 6.57, P = 0.015 with VP339; F3,8 = 5.79, P = 0.021 with VP521; F3,8 = 15.83, P = 0.001 with VP449) and the highest ratio of omega-3 and omega-6 FAs (9.4:1 in the total yield).
Fatty acid composition of new Coccomyxa strains.
Fatty acid
C. subellipsoidea VP336
C. fusiformis VP339
C. tropica VP521
C. cattiensis VP449
14:0 Myristic
1.22ABC ± 0.07
0.16A ± 0.01
0.16B ± 0.02
0.24C ± 0.01
16:0 Palmitic
22.53A ± 0.37
13.59AB ± 0.31
22.01B ± 0.59
31.16AB ± 1.03
16:1n-9 Hypogeic
0.16Aa ± 0.01
0.29aB ± 0.02
0.65AB ± 0.04
0.11B ± 0.02
16:1n-7 Palmitoleic
3.18AB ± 0.09
2.49Aa ± 0.09
2.10aB ± 0.09
1.29AB ± 0.04
16:2n-6 7,10-Hexadecadienoic
0.26A ± 0.02
4.95A ± 0.14
2.90A ± 0.09
3.83A ± 0.11
16:3n-3 Roughanic
10.77A ± 0.33
1.44A ± 0.06
16.12A ± 0.39
7.70A ± 0.21
16:3n-6 4,7,10-Hexadecatrienoic
4.96 ± 0.14
16:4n-3 4,7,10,13-Hexadecatetraenoic
3.95A ± 0.14
17.16A ± 0.57
18:0 Stearic
1.74A ± 0.07
0.31A ± 0.03
1.04A ± 0.04
2.91A ± 0.10
18:1n-9 Oleic
1.73A ± 0.05
8.02A ± 0.24
5.60A ± 0.14
3.03A ± 0.07
18:1n-7 Vaccenic
0.48a ± 0.04
0.68a ± 0.04
18:2n-6 Linoleic
4.54AB ± 0.14
9.48A ± 0.12
8.88B ± 0.24
21.69AB ± 0.34
18:3n-6 γ-Linolenic
0.11Aa ± 0.02
0.05aB ± 0.01
0.41AB ± 0.02
18:3n-3 α-Linolenic
47.02A ± 0.76
32.83A ± 0.64
39.92A ± 1.14
26.78A ± 0.62
18:4n-3 Stearidonic
0.05A ± 0.01
4.18A ± 0.12
20:0 Arachidic
0.10A ± 0.01
0.03A ± 0.01
20:4n-6 Arachidonic
1.90AB ± 0.06
0.09A ± 0.01
0.17B ± 0.02
20:5n-3 Eicosapentaenoic
0.85 ± 0.20
SFA
25.59A ± 0.36
14.09AB ± 0.29
23.21B ± 0.60
34.31AB ± 1.0
MUFA
5.07A ± 0.08
10.80AB ± 0.23
8.83AB ± 0.17
5.11B ± 0.09
PUFA
69.34A ± 1.35
75.11aB ± 1.39
67.96ab ± 1.75
60.58ABb ± 0.81
Omega-3
62.64Aab ± 1.23
55.61aB ± 1.26
56.04bC ± 1.53
34.48ABC ± 0.74
Omega-6
6.70A ± 0.12
19.50A ± 0.13
11.92A ± 0.23
26.10A ± 0.29
Omega-3/omega-6
9.3
2.6
4.7
1.3
The data are reported as the mean (% of total fatty acids) ± standard error from three independent biological replicates.
A,B,CMean with unlike superscript in row differ significantly (P ≤ 0.01). a,bMean with unlike superscript in row differ significantly (P ≤ 0.05).
Fatty acid composition of new Coccomyxa strains.
Fatty acid
C. subellipsoidea VP336
C. fusiformis VP339
C. tropica VP521
C. cattiensis VP449
14:0 Myristic
4.65AB ± 0.16
2.78AB ± 0.11
0.97B ± 0.03
0.57A ± 0.01
16:0 Palmitic
86.0A ± 3.06
237.0AB ± 6.66
131.0AB ± 4.51
74.57B ± 1.78
16:1n-9 Hypogeic
0.61Aa ± 0.02
5.11AB ± 0.13
3.89AB ± 0.05
0.27aB ± 0.01
16:1n-7 Palmitoleic
12.14A ± 0.27
43.51AB ± 0.71
12.50B ± 0.28
3.08AB ± 0.03
16:2n-6 7,10-Hexadecadienoic
1.01A ± 0.04
86.61A ± 1.50
17.21A ± 0.98
9.19A ± 0.53
16:3n-3 Roughanic
41.10A ± 0.96
25.26A ± 0.31
95.96A ± 1.9
18.44A ± 0.44
16:3n-6 4,7,10-Hexadecatrienoic
86.61 ± 0.94
16:4n-3 4,7,10,13-Hexadecatetraenoic
15.06A ± 0.53
299.49A ± 9.77
18:0 Stearic
6.67A ± 0.20
5.34AB ± 0.09
6.19B ± 0.08
6.96B ± 0.05
18:1n-9 Oleic
6.61A ± 0.11
139.87AB ± 2.18
33.30AB ± 0.83
7.25B ± 0.15
18:1n-7 Vaccenic
2.88A ± 0.03
1.63A ± 0.05
18:2n-6 Linoleic
17.32AB ± 0.33
165.40AB ± 3.99
52.87A ± 0.72
51.90B ± 0.64
18:3n-6 γ-Linolenic
1.89A ± 0.03
0.30A ± 0.01
0.98A ± 0.01
18:3n-3 α-Linolenic
179.44A ± 3.51
572.85A ± 11.66
237.65A ± 5.17
64.09A ± 0.84
18:4n-3 Stearidonic
0.19A ± 0.01
72.94A ± 0.65
20:0 Arachidic
0.37A ± 0.01
0.56A ± 0.02
20:4n-6 Arachidonic
7.24AB ± 0.36
0.55A ± 0.01
0.42B ± 0.01
20:5n-3 Eicosapentaenoic
3.24 ± 0.05
SFA
97.69A ± 3.42
245.68AB ± 6.86
138.16AB ± 4.60
82.10B ± 1.84
MUFA
19.36A ± 0.40
188.49AB ± 2.96
52.57AB ± 1.18
12.23B ± 0.24
PUFA
264.60A ± 5.72
1311.05A ± 28.72
404.54A ± 8.73
145.02A ± 2.42
Omega-3
239.03A ± 5.01
970.54A ± 22.37
333.61A ± 7.03
82.53A ± 1.27
Omega-6
25.57AB ± 0.72
340.51AB ± 6.36
70.93A ± 1.71
62.49B ± 1.19
Omega-3/omega-6
9.4
2.9
4.7
1.3
The data are reported as the mean (mg L–1) ± standard error from three independent biological replicates.
A,BMean with unlike superscript in row differ significantly (P ≤ 0.01). aMean with unlike superscript in row differ significantly (P ≤ 0.05).
4 Discussion
The results of recent studies of Coccomyxa using molecular methods have shown that identification based solely on morphological characteristics is possible up to the genus level (Darienko et al., 2015; Malavasi et al., 2016; Barcytë and Nedbalová, 2017; Cao et al., 2018b). Species identification is possible only with an integrative approach (Darienko et al., 2015). The experiments of Darienko et al. (2015) demonstrated changes in morphological characteristics under different physiological conditions. Coccomyxa cells changed size and shape depending on the concentration of NaCl in the medium. Morphological plasticity has also been shown for C. silvae-gabretae from Plešné Lake in the Czech Republic (Barcytë and Nedbalová, 2017). Under natural conditions, the species had elongated spindle-shaped cells, while under laboratory conditions, these features disappeared, and the cells acquired an ellipsoidal shape. Therefore, morphometric characteristics cannot be considered diagnostic features for Coccomyxa species identification.
The formation of a one-side mucilage cap was previously considered a distinctive feature of the genus (Andreyeva, 1998). Later, as a result of the phylogenetic analysis of Darienko et al. (2015), it was shown that two closely related genera, Coccomyxa and Pseudococcomyxa, represent a monophyletic group. Coccomyxa and Pseudococcomyxa species were previously separated according to the ability of cells to form mucilage (Andreyeva, 1998). In this regard, the ability of mucilage formation in Coccomyxa-like strains under starvation conditions was tested (Darienko et al., 2015). The authors investigated nine strains representing all the phylogenetic lineages of Coccomyxa and showed that after 6 weeks of starvation, only 2 out of 9 strains were observed to produce mucilaginous sheath. Thus, it has been proven that mucilage formation is not characteristic of all Coccomyxa species. In recent years, the absence of mucilage has been shown for several new species, such as C. antarctica (Cao et al., 2018a), C. greatwallensis (Cao et al., 2018b), C. fottii, C. silvae-gabretae (Barcytë and Nedbalová, 2017).
Determining species boundaries is based on DNA sequence analysis, considering morphology and ecology (integrative approach). The first paper on a detailed study of Coccomyxa was published by Darienko et al. (2015). Based on the analysis of the molecular structure of 18S rRNA and ITS2 markers of 41 strains from public and working collections, it was proposed to distinguish 7 Coccomyxa species. Later, as additional criteria for the species, Malavasi et al. (2016) proposed to evaluate environmental data. As a result of the analysis of 61 sequences (18S rDNA and ITS1–5.8S rDNA–ITS2) of Coccomyxa, it was shown that the living stage and habitat ecological characteristics of species are consistent with the phylogenetic lineages (Malavasi et al., 2016). Other authors also support this approach (Barcytë and Nedbalová, 2017; Sciuto et al., 2019).
All studied strains were isolated from the upper 5 cm layer of soil on four different sampling sites, of which 3 were forests and 1 was a dry swamp (Figure 1 and Table 1). Only C. cattiensis was found at two points; the remaining species were recorded once. The described species have a typical morphology for Coccomyxa: small cells with ellipsoidal or ovoid shape, parietal chloroplast shape without a pyrenoid, and oblique cell division. A specific characteristic of C. fusiformis is the presence of cells with large sizes and spindle-like shape, reminiscent of the crescent-shaped cells of Monoraphidium (Figure 5). We investigated all strains under starvation conditions within 10 weeks with cultivation in a 1:20 diluted 3N-BBM + V medium containing 1% glucose. The mucilage production in the culture cells was not observed.
Our phylogenetic analysis includes 18S rDNA and ITS1–5.8S rDNA–ITS2 sequences of 80 Coccomyxa strains. The tree topology is consistent with previous studies (Darienko et al., 2015; Malavasi et al., 2016; Barcytë and Nedbalová, 2017; Sciuto et al., 2019). Coccomyxa species form a monophyletic group, and previously designated clades are supported (Figure 2). New species described in this study form sister lineages to clades G and C. subellipsoidea (C. fusiformis), clade I (C. tropica) sensu Malavasi et al. (2016) and clade “C. parasitica” (C. cattiensis) sensu Sciuto et al. (2019) (Figure 2).
Phylogenetic positions in the tree are consistent with the ecological differentiation of Coccomyxa lineages, as suggested by Malavasi et al. (2016). Most strains in the clade C. subellipsoidea are photobionts of soil lichens, rarely epiphytic. Terrestrial and epilithic taxa represent neighboring clades E, G, and H (Malavasi et al., 2016). New strain C. subellipsoidea VP336 forms the lineage next to the terrestrial taxa of the clade E: Coccomyxa sp. NIES 2252 (unknown habitat), Coccomyxa sp. IB-GF-12 (soil habitat). New species C. fusiformis forms the lineage next to the epilithic strain CAUP H5105 from the clade G in Malavasi et al. (2016). The evolutionary distance matrix based on 18S rDNA–ITS1–5.8S rDNA region showed that C. subellipsoidea shared 99.0–99.3% similarities with other Coccomyxa strains from the clade C. subellipsoidea sensu Darienko et al. (2015) (Figure 12) and 94.9–98.7% similarities with Coccomyxa species from other clades. The p-distance matrix showed that C. fusiformis shared 98.8–99.1% similarities with strains from the clade C. subellipsoidea, and 94.5–98.9% similarities with strains from clades G, C. simplex and C. polymorpha (Figure 12).
Pairwise nucleotide identity (p-distance) heatmap based on the 18S rDNA–ITS1–5.8S rDNA region (2,203 bp) showing 10 Coccomyxa strains including VP336 and VP339. Red indicates a high nucleotide identity, and blue indicates a lower nucleotide identity value. The scale bar represents color codes of nucleotide identity values, showing sequence identity percentage. See Supplementary material 3, Table 1 for the exact p-distance values.
C. tropica complements clade I with high statistical support (92ML/0.94BI), most closely related to the tropical epiphytic strain Coccomyxa sp. KN-2011-T3. Pairwise comparisons of 18S rDNA–ITS1–5.8S rDNA region with taxa from the clade C. polymorpha sensu Darienko et al. (2015) showed that C. tropica was 98.0–100% similar Coccomyxa strains (Figure 13). The highest sequence similarity (100%) was recorded with strain Coccomyxa sp. KN-2011-T3. Additionally, the absence of CBC, hCBC, and single bases within the barcode region of the ITS2 suggests that strain KN-2011-T3 belongs to C. tropica as C. cf. tropica. Another strain in this clade, Coccomyxa sp. LH08AW1017 was isolated from soil in Germany. Thus, the clade I was formed by Coccomyxa from terrestrial habitat.
Pairwise nucleotide identity (p-distance) heatmap based on the 18S rDNA–ITS1–5.8S rDNA region (2,227 bp) showing 11 Coccomyxa strains including VP521. Red indicates a high nucleotide identity, and blue indicates a lower nucleotide identity value. The scale bar represents color codes of nucleotide identity values, showing sequence identity percentage. See Supplementary material 3, Table 2 for the exact p-distance values.
C. cattiensis is represented by VP449 and VP451, which are closely related to the free-living acid-tolerant strain Coccomyxa sp. CPCC 508, isolated from the metal-contaminated area in Boomerang Lake, near Red Lake (ON, Canada) (Verma et al., 2009). The highest 18S rRNA sequence similarity (100%) between C. cattiensis and Coccomyxa sp. CPCC 508 (Figure 14) suggests this strain belongs to the described species as C. cf. cattiensis. Neighboring sequences from the clade C. parasitica sensu Sciuto et al. (2019) belong to strains infesting in mussels and probably represent separate species from C. cattiensis. Strain CP-01 was isolated from horse mussel Modiolus kurilensis (Sokolnikova et al., 2016). The two Flensburg fjord strains (1 and 2) were also attributed to the species C. parasitica and characterized as a facultative parasite blue mussel Mytilus edulis (Rodríguez et al., 2008). C. cattiensis showed sequence similarity of 99.6–99.7% to “C. parasitica” Flensburg fjord (1 and 2) and 99.4% to C. antarctica FACHB-2140 (Figure 14). Given the ecological and phylogenetic characteristics, the new species C. cattiensis is separated from other members of the clade C. parasitica sensu (Sciuto et al., 2019).
Pairwise nucleotide identity (p-distance) heatmap based on the 18S rRNA gene (1,788 bp) showing 13 Coccomyxa strains including VP449 and V451. Red indicates a high nucleotide identity, and blue indicates a lower nucleotide identity value. The scale bar represents color codes of nucleotide identity values, showing sequence identity percentage. See Supplementary material 3, Table 3 for the exact p-distance values.
Obtained FA profiles were compared with the published evidence on FA profiles of other Coccomyxa strains in Supplementary material 2. The analysis compared the FA profiles of Coccomyxa strains representing different habitats and varying culture times. The study included strains C. elongata SAG 216-3a, C. mucigena SAG 216-4, C. simplex SAG 216-8, SAG 216-9a, C. solorinae SAG 216-5, and C. subellipsoidea SAG 216-7, which maintained in culture since 1954, C. subellipsoidea SAG 69.80 since 1971, C. dispar SAG 49.84 since 1983 and C. elongata MZ–Ch64 since 2015. Coccomyxa strains of lichen photobionts from the lineage C. solorinae revealed a high content of PUFAs in the range of 55.60–75.30% of total fatty acids and a low concentration of MUFAs (Lang et al., 2011). Free living strains from different habitats and clades (freshwater C. elongata MZ–Ch64, epiphyte C. avernensis SAG 216-1 and Coccomyxa sp. SAG 2040) showed enrichment with MUFAs: 60.09% (Maltsev et al., 2019), 46.02 and 42.94% (Lang et al., 2011), respectively. The leader in SFA accumulation was the strain C. melkonianii SCCA048 (45.50%) from a river highly polluted by heavy metals (Soru et al., 2019). Conversely, this FA group’s minimum values (8.83–9.64%) are shown for Coccomyxa strains of lichen photobionts from the lineage C. subellipsoidea (Lang et al., 2011). Within PUFAs special attention is paid to omega-3 and omega-6 FAs since they are characterized by different features and play an essential role in providing regular activity for humans and animals (Horrocks and Yeo, 1999). New species from the tropical forest soil (C. tropica) together with strains C. subellipsoidea SAG 216-7 and SAG 69.80 accumulated omega-3 FAs in the range of 54.75–62.64% with maximum value for strain C. subellipsoidea VP336. Compared to oils obtained from other productive green algae strains (Maltsev and Maltseva, 2021), oil produced by the strain VP336 has a similar content of omega-3 FAs with Chlamydomonas reinhardtii strains SAG 11-32a, SAG 11-32b, and SAG 73.72 (62.80–63.81% of omega-3). On the other hand, new Coccomyxa strains are characterized by a more minor concentration of omega-3 FAs in contrast to Chlamydomonas hydra strain SAG 11-6b and SAG 11-6c (76.63–85.26%, respectively). Coccomyxa strains (SAG 216-14, SAG 2040) from the clade C. viridis revealed a high content of omega-6 FAs (40.0%, 36.54%) and the lowest content of omega-3 FAs (Supplementary material 2). To sum it up, new Coccomyxa strains accumulate remarkably high amounts of PUFAs (up to 655.55 mg L–1) compared to Coccomyxa onubensis with 161.2 mg L–1 (Bermejo et al., 2018) and such green algae strain as Coelastrella multistriata MZ–Ch23 with 300.1 mg L–1 (Maltsev et al., 2021). The ratio of omega-3 and omega-6 FAs is essential for the biotechnological use of algal biomass. Most notably, it should have high omega-3 FA concentrations accompanied by low omega-6 FA content. Different values of this ratio characterize Coccomyxa strains: extra low content omega-3 over omega-6 with the factor of 0.148–0.969 (strains from the clade C. viridis), slight predominance omega-3 over omega-6 by 1.551–3.1-fold (strains from the clade C. simplex, except MZ–Ch64, SAG 216-2, SCCA048), significant predominance omega-3 over omega-6 by 2.637–9.342-fold (C. fusiformis, C. tropica and strains from the clade C. subellipsoidea, except NIES 2166 from Antarctica). We found the maximum value of the omega-3/omega-6 ratio in strain C. subellipsoidea VP336 (Table 2; Supplementary material 2). Generally, the analysis shows equitable trends in fatty acid composition among strains from similar habitats and close phylogenetic clades. At the same time, the long-term culture time of Coccomyxa is not accompanied by a loss of specific biochemical characteristics. For example, strains C. subellipsoidea SAG 216-7 and SAG 69.80 are lichen photobionts and belong to the same phylogenetic subclade (Figure 2), but were isolated from geographically distant areas (Finland and Germany) and maintained as cultures for different times (since 1954 and 1970). However, these strains are characterized by similar values of omega-3 FA content (56.39% and 54.75%) and omega-3/omega-6 ratio (5.3 and 4.4). Among other trebouxiophycean species, the highest omega-3 to omega-6 ratio was observed in Neglectella solitaria SAG 83.80 (30.6:1), Oocystis parva SAG 82.80 (28.7:1) and Koliella longiseta SAG 470-1 (9.1:1). The opposite excess of omega-6 over omega-3 fatty acids has been reported to Parachlorella kessleri SAG 10.80 (14.4:1), Micractinium pusillum SAG 13.81 (10.3:1) and Myrmecia bisecta SAG 2043 (9.3:1) (Lang et al., 2011; Maltsev and Maltseva, 2021).
A well-balanced proportion of omega-3 and omega-6 FAs and a high content of 18:3n-3 α-linolenic acid determine the possibility of using biomass of new Coccomyxa species from forest soils of Vietnam (especially strains VP336 and VP521) as a component of highly effective food additives or forage in aquaculture and farming. Enrichment of feed with fatty acids is a new strategy for agriculture (Ma et al., 2019). The high content of SFAs and PUFAs (up to 94.89%) positions the novel strain C. cattiensis VP449 as a promising feed additive that promotes methane mitigation from ruminants (Beauchemin et al., 2022). The favorable effect of the algal biomass inclusion in the diet was observed for growing various birds and farm animals (Lamminen et al., 2019; Petrolli et al., 2019). Fortification of the cows’ diet with microalgae-based feed resulted in the successful enrichment of milk with omega-3 docosahexaenoic acid up to 4.5–6.4 mg/100 mL milk (Moran et al., 2019). In the future, it will also be necessary to carry out cultivation experiments to evaluate growth dynamics and the impact of stress factors on fatty acid profiles and lipid productivity.
Thus, phylogenetic analysis and considering habitat ecological characteristics confirm the description of three new species of Coccomyxa. The recognition and description of these novel species suggest that the algal flora of Vietnam’s soil habitats has yet to be fully documented.
5 Conclusion
Investigations of soil green algae using molecular methods are still minimal. Using an integrative approach, we have described three new Coccomyxa species from an understudied region. Phylogenetic analysis based on the 18S rRNA gene indicates that new species were from separate lineages. The fatty acid composition of the studied soil Coccomyxa strains showed an increased content of palmitic, α-linolenic, and linoleic acids. Given the fatty acid content of the strain’s biomass, they are a promising natural source for feed production, promoting methane mitigation from ruminants and the nutraceutical and pharmaceutical industries. In future studies, adjustments to lighting intensity, the composition of the medium, and two-stage cultivation will be examined to optimize cultivation conditions and improve the productivity of Coccomyxa strains.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
YM: Conceptualization, Formal Analysis, Methodology, Writing – original draft. EK: Conceptualization, Formal Analysis, Investigation, Methodology, Resources, Visualization, Writing – original draft. SM: Resources, Validation, Visualization, Writing – original draft. ZK: Investigation, Writing – original draft. CÐ: Writing – review and editing. MK: Conceptualization, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This publication was based on research carried out with financial support by Russian Science Foundation (project number 23-74-10081), percentage contribution was 80%. Algal strain’s cultivation was obtained with financial support by state assignment of the Ministry of Science and Higher Education of the Russian Federation (theme 124052200012-7 No. FFES-2024-0001), percentage contribution was 20%.
We thank A.A. Kotov (A.N. Severtsov Institute of Ecology and Evolution) for collecting samples. The expedition was organized and permitted by the Joint Russian-Vietnam Tropical Center, Ecolan 1.2 theme.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1517865/full#supplementary-material
https://www.ncbi.nlm.nih.gov
ReferencesAbeK.IshiwatariT.WakamatsuM.AburaiN. (2014). Fatty acid content and profile of the aerial microalga Coccomyxa sp. isolated from dry environments.Appl. Biochem. Biotechnol.1741724–1735. 10.1007/s12010-014-1181-y25146196AndreyevaV. (1998). Terrestial and Aerophytic Green Algae (Chlorophyta: Tetrasporales, Chlorococcales, Chlorosarcinales).St. Peterburg: Nauka.AnichkinA. (2011). “Soil macrofauna: Structure and seasonal dynamics,” in Structure and Functions of soil Communities of a Monsoon Tropical Forest (Cat Tien National Park, Southern Vietnam), ed.TiunovA. V. (Moscow: KMK Scientific Press), 45–75.ArinushkinaE. V. (1970). Handbook for Chemical Soil Analysis.Moscow: Moscow State University.BarcytëD.NedbalováL. (2017). Coccomyxa: A dominant planktic alga in two acid lakes of different origin.Extremophiles21245–257. 10.1007/s00792-016-0899-627942983BatesS.ClementeJ.FloresG.WaltersW.ParfreyL.KnightR. (2013). Global biogeography of highly diverse protistan communities in soil.ISME J.7652–659. 10.1038/ismej.2012.14723235291BeaucheminK. A.UngerfeldE. M.AbdallaA. L.AlvarezC.ArndtC.BecquetP. (2022). Invited review: Current enteric methane mitigation options.J. Dairy Sci.1059297–9326. 10.3168/jds.2022-2209136270879BermejoE.Ruiz-DomínguezM. C.CuaresmaM.VaqueroI.Ramos-MerchanteA.VegaJ. M. (2018). Production of lutein, and polyunsaturated fatty acids by the acidophilic eukaryotic microalga Coccomyxa onubensis under abiotic stress by salt or ultraviolet light.J. Biosci. Bioengin.125669–675. 10.1016/j.jbiosc.2017.12.02529370982Bio-RadL. (2020). Chelex§-100 and Chelex§-20 Chelating Ion Exchange Resin Instruction Manual.Hercules, CA: Bio-Rad Laboratories.BischoffH. W.BoldH. C. (1963). Phycological Studies IV. Some Soil Algae From Enchanted Rock and Related Algal Species.Austin, TX: University of Texas Publication.ByunY.HanK. (2009). PseudoViewer3: Generating planar drawings of large–scale RNA structures with pseudoknots.Bioinformatics251435–1437. 10.1093/bioinformatics/btp25219369500CaisováL.MarinB.MelkonianM. (2013). A consensus secondary structure of ITS2 in the Chlorophyta identified by phylogenetic reconstruction.Protist164482–496. 10.1016/j.protis.2013.04.00523770573CaoS.ZhangF.ZhengH.LiuC.PengF.ZhouQ. (2018a). Coccomyxa antarctica sp. nov. from the Antarctic lichen Usnea aurantiacoatra.PhytoKeys98107–115. 10.3897/phytokeys.98.2536029849474CaoS.ZhangF.ZhengH.PengF.LiuC.ZhouQ. (2018b). Coccomyxa greatwallensis sp. nov. (Trebouxiophyceae, Chlorophyta), a lichen epiphytic alga from Fildes Peninsula, Antarctica.PhytoKeys110:39. 10.3897/phytokeys.110.2696130473613ColemanA. W. (2007). Pan–eukaryote ITS2 homologies revealed by RNA secondary structure.Nucl. Acids Res.353322–3329. 10.1093/nar/gkm23317459886DarienkoT.GustavsL.EggertA.WolfW.PröscholdT. (2015). Evaluating the species boundaries of green microalgae (Coccomyxa, Trebouxiophyceae, Chlorophyta) using integrative taxonomy and DNA barcoding with further implications for the species identification in environmental samples.PLoS One10:e0127838. 10.1371/journal.pone.012783826080086DarribaD.TaboadaG. L.DoalloR.PosadaD. (2012). jModelTest 2: More models, new heuristics and parallel computing.Nat. Meth.9:772. 10.1038/nmeth.210922847109DrummondA. J.RambautA. (2007). BEAST: Bayesian evolutionary analysis by sampling trees.BMC Evol. Biol.7:214. 10.1186/1471-2148-7-21417996036EliašM.NìmcováY.ŠkaloudP.NeustupaJ.KaufnerovaV.ŠejnohováL. (2010). Hylodesmus singaporensis gen. et sp. nov., a new autosporic subaerial green alga (Scenedesmaceae, Chlorophyta) from Singapore.Int. J. Syst. Evol. Microbiol.601224–1235. 10.1099/ijs.0.012963-019666799EliášM.NeustupaJ.ŠkaloudP. (2008). Elliptochloris bilobata var. corticola var. nov. (Trebouxiophyceae, Chlorophyta), a novel subaerial coccal green alga.Biologia63791–798. 10.2478/s11756-008-0100-5FuentesJ.HussV.MonteroZ.TorronterasR.CuaresmaM.GarbayoI. (2016). Phylogenetic characterization and morphological and physiological aspects of a novel acidotolerant and halotolerant microalga Coccomyxa onubensis sp. nov. (Chlorophyta, Trebouxiophyceae).J. Appl. Phycol.283269–3279. 10.1007/s10811-016-0887-3FujitaY.OhtsukaT. (2005). Diatoms from paddy fields in northern Laos.Diatom2171–89. 10.11464/DIATOM1985.21.0_71GarbayoI.TorronterasR.ForjánE.CuaresmaM.CasalC.MogedasB. (2012). Identification and physiological aspects of a novel carotenoid-enriched, metal-resistant microalga isolated from an acidic river in Huelva (Spain).J. Phycol.48607–614. 10.1111/j.1529-8817.2012.01160.x27011076GlushchenkoA.KezlyaE.MaltsevY.GenkalS.KociolekJ. P.KulikovskiyM. (2022). Description of the soil diatom Sellaphora terrestris sp. nov. (Bacillariophyceae, Sellaphoraceae) from Vietnam, with remarks on the phylogeny and taxonomy of Sellaphora and systematic position of Microcostatus.Plants11:2148. 10.3390/plants1116214836015452GustavsL.SchiefelbeinU.DarienkoT. (2017). “Symbioses of the green algal genera Coccomyxa and Elliptochloris (Trebouxiophyceae, Chlorophyta),” in Algal and Cyanobacterial Symbioses, edsGrubeM.SeckbachJ.MuggiaL. (London: World Scientific Publishing Europe Ltd), 169–208.HameedM. S. A. (2006). Survey of soil algal flora of some cultivated soils in Beni Suef.Egypt. Egypt. J. Phycol.71–15. 10.21608/egyjs.2006.114130HodačL. (2016). Green Algae in Soil: Assessing their Biodiversity and Biogeography with Molecular-Phylogenetic Methods Based on Cultures. Ph.D. Thesis. Göttingen: Georg August University of Göttingen.HorrocksL. A.YeoY. K. (1999). Health benefits of docosahexaenoic acid (DHA).Pharmacol. Res.40211–225. 10.1006/phrs.1999.049510479465KezlyaE.GlushchenkoA.KociolekJ. P.MaltsevY.MartynenkoN.GenkalS. (2020). Mayamaea vietnamica sp. nov.: A new, terrestrial diatom (Bacillariophyceae) species from Vietnam.Algae35325–335. 10.4490/algae.2020.35.11.23KezlyaE.GlushchenkoA.MaltsevY.GusevE.GenkalS.KociolekJ. P. (2021). Three new species of Placoneis Mereschkowsky (Bacillariophyceae: Cymbellales) with comments on cryptic diversity in the P. elginensis–group.Water13:3276. 10.3390/w13223276KezlyaE.MaltsevY.GenkalS.KrivovaZ.KulikovskiyM. (2022). Phylogeny and fatty acid profiles of new Pinnularia (Bacillariophyta) species from soils of Vietnam.Cells11:2446. 10.3390/cells1115244635954290KezlyaE.TseplikN.KulikovskiyM. (2023). Genetic markers for metabarcoding of freshwater microalgae: Review.Biology12:1038. 10.3390/biology1207103837508467KhokhlovaO.MyakshinaT.KuznetsovA.GubinS. (2017). Morphogenetic features of soils in the Cat Tien national park, southern Vietnam.Eurasian Soil Sci.50158–175. 10.1134/S1064229316120085LamminenM.Halmemies-Beauchet-FilleauA.KokkonenT.JaakkolaS.VanhataloA. (2019). Different microalgae species as a substitutive protein feed for soya bean meal in grass silage based dairy cow diets.Anim. Feed Sci. Technol.247112–126. 10.1016/j.anifeedsci.2018.11.005LangI.HodacL.FriedlT.FeussnerI. (2011). Fatty acid profiles and their distribution patterns in microalgae: A comprehensive analysis of more than 2000 strains from the SAG culture collection.BMC Plant Biol.11:124. 10.1186/1471-2229-11-12421896160LawleyB.RipleyS.BridgeP.ConveyP. (2004). Molecular analysis of geographic patterns of eukaryotic diversity in Antarctic soils.Appl. Environm. Microbiol.705963–5672. 10.1128/AEM.70.10.5963-5972.200415466539LiS.SunH.HuY.LiuB.ZhuH.HuZ. (2020). Four new members of foliicolous green algae within the Watanabea clade (Trebouxiophyceae, Chlorophyta) from China.J. Eukaryot. Microbiol.67369–382. 10.1111/jeu.1278731958181LiS.TanH.LiuB.ZhuH.HuZ.LiuG. (2021). Watanabeales ord. nov. and twelve novel species of Trebouxiophyceae (Chlorophyta).J. Phycol.571167–1186. 10.1111/jpy.1316533713360LindgrenH.MoncadaB.LückingR.MagainN.SimonA.GoffinetB. (2020). Cophylogenetic patterns in algal symbionts correlate with repeated symbiont switches during diversification and geographic expansion of lichen-forming fungi in the genus Sticta (Ascomycota, Peltigeraceae).Molec. Phylogen. Evol.150:106860. 10.1016/j.ympev.2020.10686032473336MaY.ZhouW.ChenP.UrriolaP. E.ShursonG. C.RuanR. (2019). Metabolomic evaluation of Scenedesmus sp. as a feed ingredient revealed dose-dependent effects on redox balance, intermediary and microbial metabolism in a mouse model.Nutrients11:1971. 10.3390/nu1109197131438641MalavasiV.ŠkaloudP.RindiF.TempestaS.PaolettiM.PasqualettiM. (2016). DNA-based taxonomy in ecologically versatile microalgae: A re-evaluation of the species concept within the coccoid green algal genus Coccomyxa (Trebouxiophyceae, Chlorophyta).PLoS One11:e0151137. 10.1371/journal.pone.015113727028195MaltsevY.MaltsevaK. (2021). Fatty acids of microalgae: Diversity and applications.Rev. Environ. Sci. Biotechnol.20515–547. 10.1007/s11157-021-09571-3MaltsevY.GusevE.MaltsevaI.KulikovskiyM.NamsaraevZ.PetrushkinaM. (2018). Description of a new species of soil algae, Parietochloris grandis sp. nov., and study of its fatty acid profiles under different culturing conditions.Algal Res.33358–368. 10.1016/j.algal.2018.06.008MaltsevY.KrivovaZ.MaltsevaS.MaltsevaK.GorshkovaE.KulikovskiyM. (2021). Lipid accumulation by Coelastrella multistriata (Scenedesmaceae, Sphaeropleales) during nitrogen and phosphorus starvation.Sci. Rep.11:19818. 10.1038/s41598-021-99376-934615938MaltsevY.MaltsevaI.MaltsevaS.KociolekJ. P.KulikovskiyM. (2019). Fatty acid content and profile of the novel strain of Coccomyxa elongata (Trebouxiophyceae, Chlorophyta) cultivated at reduced nitrogen and phosphorus concentrations.J. Phycol.551154–1165. 10.1111/jpy.1290331318981MaltsevaS.KezlyaE.KrivovaZ.GusevE.KulikovskiyM.MaltsevY. (2022). Phylogeny and fatty acid profiles of Aliinostoc vietnamicum sp. nov. (Cyanobacteria) from the soils of Vietnam.J. Phycol.58789–803. 10.1111/jpy.1328336054495MartynenkoN.KezlyaE.GusevE. (2022). Description of a new species of the genus Cryptomonas (Cryptophyceae: Cryptomonadales), isolated from soils in a tropical forest.Diversity14:1001. 10.3390/d14111001MoranC. A.MorlacchiniM.KeeganJ. D.WarrenH.FusconiG. (2019). Dietary supplementation of dairy cows with a docosahexaenoic acid-rich thraustochytrid, Aurantiochytrium limacinum: Effects on milk quality, fatty acid composition and cheese making properties.J. Anim. Feed Sci.283–14. 10.22358/jafs/105105/2019NakadaT.ShinkawaH.ItoT.TomitaM. (2010). Recharacterization of Chlamydomonas reinhardtii and its relatives with new isolates from Japan.J. Pl. Res.12367–78. 10.1007/s10265-009-0266-019882207NeustupaJ.EliášM.ŠkaloudP.NìmcováY.ŠejnohováL. (2011). Xylochloris irregularis gen. et sp. nov. (Trebouxiophyceae, Chlorophyta), a novel subaerial coccoid green alga.Phycologia5057–66. 10.2216/08-64.1NeustupaJ.NìmcováY.EliášM.ŠkaloudP. (2009). Kalinella bambusicola gen. et sp. nov. (Trebouxiophyceae, Chlorophyta), a novel coccoid Chlorella-like subaerial alga from Southeast Asia.Phycol. Res.57159–169. 10.1111/j.1440-1835.2009.00534.xPetrolliT. G.PetrolliO. J.PereiraA. S. C.ZottiC. A.RomaniJ.VillaniR. (2019). Effects of the dietary supplementation with a microalga extract on broiler performance and fatty-acid meat profile.Braz. J. Poult. Sci.21eRBCA–2018–0958. 10.1590/1806-9061-2018-0958ØezankaT.NedbalováL.SiglerK. (2008). Unusual medium-chain polyunsaturated fatty acids from the snow alga Chloromonas brevispina.Microbiol. Res.163373–379. 10.1016/j.micres.2006.11.02117403601RodríguezF.FeistS.GuillouL.HarkestadL.BatemanK.RenaultT. (2008). Phylogenetic and morphological characterisation of the green algae infesting blue mussel Mytilus edulis in the North and South Atlantic oceans.Dis. Aquat. Org.81231–240. 10.3354/dao0195618998587SchmidleW. (1901). Über drei Algengenera.Berichte der Deutschen Botanischen Gessellschaft1910–24.SciutoK.BaldanB.MarcatoS.MoroI. (2019). Coccomyxa cimbrica sp. nov., a green microalga found in association with carnivorous plants of the genus Drosera L.Eur. J. Phycol.54531–547. 10.1080/09670262.2019.1618920ŠkaloudP.FriedlT.HallmannC.BeckA.Dal GrandeF. (2016). Taxonomic revision and species delimitation of coccoid green algae currently assigned to the genus Dictyochloropsis (Trebouxiophyceae, Chlorophyta).J. Phycol.52599–617. 10.1111/jpy.1242227135898SokolnikovaY.MagarlamovT.StenkovaA.KumeikoV. (2016). Permanent culture and parasitic impact of the microalga Coccomyxa parasitica, isolated from horse mussel Modiolus kurilensis.J. Invert. Pathol.14025–34. 10.1016/j.jip.2016.07.01227469189SoruS.MalavasiV.CaboniP.ConcasA.CaoG. (2019). Behavior of the extremophile green alga Coccomyxa melkonianii SCCA 048 in terms of lipids production and morphology at different pH values.Extremophiles2379–89. 10.1007/s00792-018-1062-330341564StamatakisA.HooverP.RougemontJ. (2008). A rapid bootstrap algorithm for the RAxML web–servers.System. Biol.57758–771. 10.1080/1063515080242964218853362VadjuninaA. F.KorchaginaZ. A. (1986). Methods of Studying the Physical Properties of Soils.Moscow: Agropromizdat.VermaV.BhattiS.HussV. A.ColmanB. (2009). Photosynthetic inorganic carbon acquisition in an acid-tolerant, free-living species of Coccomyxa (Chlorophyta).J. Phycol.45847–854. 10.1111/j.1529-8817.2009.00718.x27034214VijayanD.RayJ. G. (2015). Green algae of a unique tropical wetland, Kuttanadu, Kerala, India, in relation to soil regions, seasons, and paddy growth stages.Int. J. Sci. Environm. Technol.4770–803. 10.4236/ajps.2015.618288WhiteT. J.BrunsT.LeeS.TaylorJ. W. (1990). “Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics,” in PCR Protocols: A Guide to Methods and Applications, edsInnisM. A.GelfandD. H.SninskyJ. J.WhiteT. J. (New York: Academic Press Inc), 315–322.ZukerM. (2003). Mfold web server for nucleic acid folding and hybridization prediction.Nucl. Acids Res.313406–3415. 10.1093/nar/gkg59512824337