# North by Southwest: Screening the Naturally Isolated Microalgal Strains from Different Habitats of Iran for Various Pharmaceutical and Biotechnology Applications

Authors: Zahra Dinpazhooh; Seyyed Vahid Niknezhad; Fardin Fadaei; Saeedeh Shaker; Ghasem Najafpour; Younes Ghasemi; Pegah Mousavi; Mohammad Hossein Morowvat
Journal: International Journal of Microbiology
 (2022)

---

## Abstract

Background and Aims. Microalgae are known as a promising source for food, pharmaceutical, and biofuel production while providing environmental advantages. The present study evaluates some newly isolated microalgal strains from north and southwest of Iran as a potential source for high-value products. Methods. Primitive screening was carried out regarding growth parameters. The molecular and morphological identifications of the selected strains were performed using 18S rRNA gene sequencing. After phylogenic and evolutionary studies, the selected microalgal strains were characterized in terms of protein and pigment content, in addition to the fatty acid profile content. Besides, the CO2 fixation rate was determined to assess capability for various environmental applications. Results. All of the selected strains were predominantly belonging to Scenedesmus sp. and Desmodesmus sp. The isolated Scenedesmus sp. VN 009 possessed the highest productivity content and CO2 fixation rate of 0.054 g·L−1d−1 and 0.1 g·L−1d−1, respectively. Moreover, data from GC/MS analysis demonstrated the high robustness of this strain to produce several valuable fatty acids including α-linolenic acid and linoleic acid in 45% and 20% of total fatty acids. Conclusions. The identified strains have a great but different potential for SCP, β-carotene, and ω-3 production, as well as CO2 fixation for environmental purposes. In this study, considering the wide range of microalgal strains in different habitats of Iran, the potential applications of native microalgae for various pharmaceutical, food, and biotechnology purposes were investigated.


---

## Body

## 1. Introduction

Nowadays, there is an increased affinity to use microalgae for a variety of industrial purposes including food, pharmaceutical, and biofuel production. Thanks to its renewability and being environmentally friendly, microalgae are considered as a robust candidate for CO2 bio-fixation, and also bioremediation of wastewater streams [1–4]. Biofuel from microalgae gained significant attention as an alternative fuel source, however, biofuel production from microalgae is faced with some major challenges since they are not economically sustainable. On the other hand, several successful attempts have been conducted toward exploiting microalgae as a food and feed source or other pharmaceutical biotechnology applications [5–7].Microalgae possess a range of advantages including high lipid content, high growth rate; the ability to grow in poor quality waters and on non-arable lands while; at the same time, providing environmental benefits [8–10].Microalgae have demonstrated diverse adaptation in natural or man-made ecosystems such as fresh, saline water, wastewater, and soil. Microalgae widely occurred in freshwater and marine habitats such as lakes, rivers, ponds, and streams [11]. Algae have been estimated over 1 million species around the world according to AlgaeBase. Currently, about 165,000 species have been processed by AlgaeBase (https://www.algaebase.org) to date [12]. Therefore, exploiting indigenous strains that possess high dominance and adaptability to local environmental conditions should be a rational strategy to obtain novel isolates with potential applications to produce bioactive compounds.A vast range of energy and non-energy product can be obtained from microalgae based on their potential [13]. Hence, the only biodiesel or single production line approach cannot be commercially viable compared to the matrix approach. To date, the system of biomass Biorefining led to the sustainable processing of biomass into a wide spectrum of marketable products [14, 15]. The various co-products derived from the algal biorefinery like lipids, proteins, carbohydrates, polyunsaturated fatty acids (PUFAs), pigments, natural colorants, antioxidants, as well as pharmaceutical active ingredients can be produced from microalgal biomass. In addition to the biorefinery approach, nowadays, combinatorial genetic engineering based on increasing the number of the gene of interest is being addressed toward more efficient biosynthesis and more industrial reality [13].Microalgae is an emerging alternative protein source that could meet expected global protein requirements in food and feed applications. Environmental problems that are associated with the high ecological footprint of conventional agriculture practices resuscitated the interest in microalgae as a source of sustainable protein [16–18]. Freshwater strains can produce high amounts of protein up to 70% DW, which is composed of essential amino acids that are equal to or even superior compared to conventional plant sources. High quality amino acids, environmentally friendly feathers like CO2 consumption, and not needing arable land or drinking water have made microalgae as alternative food source. Nowadays, Spirulina and Chlorella have been known as rich sources of protein content (about 60%) as well as mineral, carotenoid, antioxidant, and vitamin. Spirulina is being incorporated into many foods’ formulations including bread, pasta, and dairy product, most of which, use as coloring agents or as a marketing strategy. Both Spirulina and Chlorella show very similar benefits in terms of the bioactive compounds that affect our physiology. Chlorella is higher in fat and calories when compared to spirulina while containing higher amount of minerals and vitamins. While Spirulina may have 10 percent more protein. This algal biomass has positioned firmly in the food and beverage, nutraceutical, and pharmaceutical market [18–21]. However, the industry has faced some critical challenges including undesirable flavor, color, and eventually consumer acceptance. To date, the major microalgal market can be found in the health food market in the form of encapsulated algal powder. The suitability of these microalgae as a feed supplement, mainly for Poultry and aquaculture is demonstrated by a large number of nutritional and toxicological evaluations [19, 20].Currently, single-cell protein (SCP) products from microalgae and microbial sources are actively being commercialized as a food ingredient in aquaculture [22, 23].Microalgal species are also recognized as the richest source of pigments, mainly including chlorophylla, b, and c; ß-carotene, xanthophylls, phycocyanin, and phycoerythrin. Carotenoid pigments are widely used in food, pharmaceuticals, and cosmetics. Microalgal carotenoid has been indicated various health benefits; it is converted to vitamin A in the human body. which is assist body’s immune system and also has antioxidant property that eliminates the harmful effects of free radicals which are implied several diseases such as premature aging, various forms of cancer, coronary heart disease, and arthritis [21]. Nowadays, ß-carotene has been widely commercially applied as natural food and cosmetic colorant because of an increasing global trend and consumer preference toward natural products [24]. Commercially viable carotenoids are extracted from Dunaliella salina yielding up to 400 mg ß-carotene m−2 of cultivation area under ideal conditions [25].Algal lipid content could be increased up to 75% DW by nitrogen depletion [26–28], genetic modifications [29], and exploring novel microalgae strains. Although microalgal lipid has mainly been considered as a feed for biodiesel production, still not economically competitive compared to terrestrial plants and lacks an evaluation of its full potential [30]. Recently obtaining high-value compounds from algal lipids, especially polyunsaturated fatty acids (PUFA), has aroused much more attention in the food and pharmaceutical industries [31–33]. The polyunsaturated fatty acid comprises some major fatty acids, such as essential fatty acids ω-3, ω-6, and conjugated fatty acids [34]. Docosahexaenoic acid (DHA; 22 : 6 n − 3) and eicosapentaenoic acid (EPA; 20 : 5 n − 3) along with α-linolenic acid (ALA; 18 : 3 n − 3) and docosapentaenoic acid (22 : 5 n − 3) known as essential fatty acids (EFAs), which are imperative for human health [35]. DHA possesses a cardioprotective function; as a consequence, it is commonly used as a nutraceutical agent in nutritional products such as infant formula, dairy, and certain other food categories [21, 36]. There are extensive studies using the strain of Schizochytrium sp. and Crypthecodinium cohnii for scale-up production of DHA [37–39].Owing to its climate diversity as well as a variety of natural resources, Iran has been considered as one of the excellent sources to find a large number of unstudied microalgal strains with special metabolic abilities [34, 40]. However, there has not been a comprehensive study on the screening of naturally isolated microalgae to assess high-value products, a possibility for large-scale production. Hence, this study was conducted to investigate naturally isolated microalgal strains from some habitats close to the Caspian Sea (north of Iran) and Kohgiluyeh and Boyer-Ahmad provinces (southwest of Iran), as precious reservoirs for different microalgal strains, to identify the best strain regarding biomass productivity, biochemical compositions along with providing the environmental benefits, to assess information for industrial application of high-value products.


## 2. Materials and Methods

### 2.1. Isolation and Cultivation

The microalgal strains used in this study were aseptically isolated, from various habitats including rivers, headwaters, paddy fields, and a sea, which are located in Kohgiluyeh and Boyer-Ahmad provinces (30 40′ 12″ N, 51 36′ 0″ E) and Mazandaran province (36° 42′ 9″ N, 52° 39′ 27″ E), the southwest and north of Iran, respectively. Water and wet soil samples were taken from the water surface and up to 5 cm of the top of the soil.Preliminary cultivation and isolation of the microalgal strains were performed in Pharmaceutical Sciences Research Center, Shiraz University of Medical Sciences. The collected samples were enriched using BG-11 liquid medium described in Table1 [41]. They were incubated under constant illumination at 60 mol·m−2·s−1 intensity with 12 : 12 h dark: light photoperiod of white fluorescent lamps at 25 ± 2°C for three weeks. After a visible growth occurred, the purified colonies were obtained by repeated sub-culturing cells across BG-11 solid agar plates and regular microscopic observation (Figure 1). The pure unialgal colonies were inoculated into a new 20 mL flask containing liquid BG-11 medium, then incubated in the abovementioned culture room.Table 1
BG-11 culture medium composition.
CompoundAmount (g/L)NaNO31.5K2HPO40.04MgSO4.7H2O0.075CaCl2.2H2O0.036Citric acid0.006Ferric ammonium citrate0.006EDTA0.001Na2CO30.002Trace-elements solution∗1mL/L∗H3BO3, 2.86 g/L; MnCl2.4H20, 1.81 g/L; ZnSO4.7H20, 0.222 g/L; Na2MoO4.2H2O, 0.39 g/L; CuSO4.5H2O, 0.079 g/L; Co (NO3)2 .6H2O, 0.0494 g/L.Figure 1
An agar plate streaked showing small isolated colonies arising from a single cell.


### 2.2. Morphological Identification

Preliminary identification of isolated microalgae was performed by direct samples microscopy of BG-11 culture (10μL) between 10 to 14 days. Among forty isolated microalgae (Figure 2), a total of twenty-nine single-cell strains were identified by evaluation of morphologic traits (cell shape, size, chloroplast, pyrenoid, flagellum, and cells arrangement) using an Olympus BX41 optical microscope equipped with an Olympus DP80 CCD camera, based on standard morphological feature keys [42, 43].Figure 2
Light microscopy photograph of all isolated microalgal strains using a 40X microspore lens.


### 2.3. Molecular Identification of Selected Microalgae

#### 2.3.1. Primers and PCR Condition

A polymerase chain reaction (PCR) targeting 18S rRNA was used to identify the isolated microalgae. Genomic DNA was extracted using a DNP™ DNA extraction kit (SinaClon, Tehran, Iran) according to the manufacturer’s instructions and adjusted to the final concentration of 20 ng·μL−1.The partial 18S rRNA gene was PCR amplified using a universal primer set; 5′-GTCAGAGGTGAAATTCTTGGATTTA-3′ (Forward), and 5′-AGGGCAGGGACGTAATCAACG-3′ (Reverse) [40] in a total reaction volume of 80 μL containing: 40 μL Taq DNA polymerase master mix red 2X (Ampliqon, Odense, Denmark), 3 μL of each primer (20 μM), 6 μL template DNA and 28 μL sterile deionized distilled water. The PCR conditions were as follows: initial denaturation at 95°C for 5 min, followed by 30 cycles of denaturation at 95°C for 1 min, annealing at 58°C for 1 min, elongation at 72°C for 90 s, and the final extension step at 72°C for 90 s, according to the method used by Radha et al. [44] with slight modifications.


#### 2.3.2. Electrophoresis and Gel Purification

The 1.8% Tris-Borate-EDTA (TBE) agarose gel and U: Genius3 GelDoc system (SYNGENE, Cambridge, United Kingdom) was used for electrophoresis and visualization of PCR products. The desired amplicons (∼750 bp) were extracted and purified using a Gel Purification Kit (Bioneer, Daejeon, South Korea) and sent for sequencing (Macrogen, Seoul, South Korea).


### 2.4. Analyzing Sequences

#### 2.4.1. Multiple Sequence Alignment

The BioEdit (version 7.1.9) and Chromas Pro (version 2.1.3) programs were utilized to modify the obtained sequences. The final resulting 18S rRNA gene sequence was analyzed and used the BLASTn search tool for finding similarity searches against previously published 18S rRNA gene sequences in the NCBI databases. The representative sequences were selected and aligned using CLC Genomics workbench V20.0.


#### 2.4.2. Molecular Phylogeny

A phylogenetic tree was constructed with the Neighbor-Joining algorithm using CLC genomics software version 20.0 with a bootstrap analysis of 500 replicates. The 18S rRNA gene sequence obtained in this study was deposited in GenBank databases.


### 2.5. Screening of Isolated Microalgae for their Growth Potential

#### 2.5.1. Cell Dry Weight Measurement

A total of twenty-nine isolated strains were inoculated in 20 mL flasks (10% (v/v)) containing 12 mL of BG-11 medium and were kept on a rotary shaker at 150 rpm and 20°С under a 12 : 12 h LD cycle of 100μE·m−2·S−1 illumination.The experiment was performed for 18 days. The sampling procedure was exploited every two days. The biomass concentrations were assessed gravimetrically. Biomass was harvested at 10000 rpm, for 5 min at 4°C and washed twice with distilled water followed by drying at 60°C in an oven for 48 h. Eventually, ten candidate microalgae were selected from the above step based on their high biomass productivity for further studies.


#### 2.5.2. Growth Measurement and Kinetics

The biomass productivity (BP, gL−1·day−1), maximum specific growth rate (μmax day−1), and doubling time (h) were determined using equations (1), (2), and (3), respectively:(1)BP=X2−X1t2−t1,(2)μmax=lnX2−lnX1t2−t1,(3)td=ln2µmax,where X1 and X2 represent the dry cell weight (DCW, gL−1) and (t2 − t1) is required time for increasing cell concentration from X1 to X2 in the exponential phase [45]. All the experiments were done in triplicates to calculate standard deviations.


### 2.6. Biochemical Analysis of Selected Strains

#### 2.6.1. Protein Content Analysis

Total protein contents were evaluated by Bradford method. A calibration curve was assessed using bovine serum albumin. Freeze-dried biomass of each sample (5 mg) that was harvested on the 18th day, was assayed for alkaline hydrolysis in 1 mL NaOH (1 N) for 1 h at 100°C using a hot plate (WiseTherm HB-R) and then centrifuged at 10000 rpm for 3 min. Protein content of each sample was calculated by measuring absorbance and comparing their absorbance at 595 nm with the bovine serum albumin standard curve [40]. This analysis was performed in triplicates and the average values were used for statistical evaluation of the results.


#### 2.6.2. Pigments Analysis

(1) Determination of Chlorophyll a, Chlorophyll b, and Total Chlorophyll. Chlorophyll concentration determination was defined according to the method described by Eijckelhoff and Dekker [46]. Three mL of the microalgae suspension was taken and centrifuged at 2500 rpm for 10 min and rinsed with acetone (80%) and centrifuged again. The definition of the amount of chlorophyll extracted in the mentioned method was determined by colorimetric assay at 668.2 and 646.2 nm using an ELISA plate reader with the following equation: (4)ChlorophyllaμgmL=12.25×A668.2−2.79×A646.2.Five mL of the microalgae suspension was taken and added to the 3 mL acetone (80%), and the definition of the amount of total chlorophyll (a + b) extracted was specified using the following equation:(5)Total Chlorophyllaa+bμgmL=7.93×A664−19.53×A647.(2) β-Carotene Extraction and Assay. β-carotene extraction from the cell pellets was performed using the ethanol/n-hexane (2 : 1, v/v). The microalgal suspension (1 mL) was collected, centrifuged at 3000 rpm for 5 min, and rinsed with ethanol/n-hexane (2/1) and distilled water. β-carotene concentration is solved in the overlapping phase with n-hexane was determined by colorimetric assay at 450 nm using ELISA plate reader [47]. The definition of the amount of β-carotene extracted in n-hexane was specified using the following equation:(6)β−caroteneμgmL=25.2×A450.


### 2.7. CO2 Fixation Rate

By removing the carbon source from BG-11 media (Na2CO3), CO2 fixation rate for each selected strain was calculated in terms of CO2 consumption by each selected strain according to the following equation:(7)RCO2mgL−1d−1=P×Ccarbon×MCO2MC,where RCO2 is the rate of CO2 fixation (mg·L−1·d−1), P is the biomass productivity in 18 days (mg·L−1·d−1), MCO2 and MC represent the molecular weight of carbon dioxide and atomic weight of carbon, respectively. Ccarbon is the carbon content of the biomass was measured by a CHNS-O elemental analyzer equipped with gas chromatography and TCD detector (Costech, ECS 4010) [48].


### 2.8. FAME Analysis

#### 2.8.1. Lipid Extraction

Fatty acid extraction was performed according to the Bligh and Dyer method with some modifications [49]. Briefly, lyophilized microalgal biomass (1 g) was boiled in 5 mL of isopropanol for 2 min and dried by vacuum rotary. A mixture of and 5 mL of chloroform was added, followed by 5 mL of potassium chloride solution (0.88%), to give a final solvent ratio of chloroform: methanol: water of 1 : 1 : 0.8. The mixture chloroform-methanol (1 : 2 v/v) was added to dried pellet in the presence of 0.015 g butylated hydroxytoluene (BHT) as an antioxidant agent, vortexed for 30 s and centrifuged (2000g, 5 min). To the above-obtained supernatant, 0.8 mL distilled water was shaken for 5 min using a separating funnel and allowed for 30 min to be separated. The total lipid was collected by concentrating the solvent phase under nitrogen gas.


#### 2.8.2. Fatty Acid Esterification and GC/MS Analysis

The acid-catalyzed esterification was carried out in the Dien–Stark apparatus and was heated to reflux for 24 h at 75°C. The reaction mixture consisted of 0.5 g of extracted fatty acid dissolved in 3 mL methanol in the presence of 0.3 mL sulfuric acid. The mixture was washed with 4 mL saturated sodium hydrogen carbonate and then dried over anhydrous sodium sulfate. To assess an oily substance ready for GC/MS analysis, the solvent was removed under nitrogen gas [40].The GC/MS analysis was performed using an Agilent 7890A GC 5977B MSD equipped with an HP-5MS capillary column (30 m × 250μm × 0.25 μm, Agilent 19091S-433). The oven temperature was programmed between 70°C (5 min) to 270°C for 10 min and the rate was 7°C·min−1. The carrier gas was helium with a flow rate of 1 mL·min−1. The identification of fatty acid peaks was carried out by comparing the obtained mass spectra with the Wiley 7n.1 library.


## 2.1. Isolation and Cultivation

The microalgal strains used in this study were aseptically isolated, from various habitats including rivers, headwaters, paddy fields, and a sea, which are located in Kohgiluyeh and Boyer-Ahmad provinces (30 40′ 12″ N, 51 36′ 0″ E) and Mazandaran province (36° 42′ 9″ N, 52° 39′ 27″ E), the southwest and north of Iran, respectively. Water and wet soil samples were taken from the water surface and up to 5 cm of the top of the soil.Preliminary cultivation and isolation of the microalgal strains were performed in Pharmaceutical Sciences Research Center, Shiraz University of Medical Sciences. The collected samples were enriched using BG-11 liquid medium described in Table1 [41]. They were incubated under constant illumination at 60 mol·m−2·s−1 intensity with 12 : 12 h dark: light photoperiod of white fluorescent lamps at 25 ± 2°C for three weeks. After a visible growth occurred, the purified colonies were obtained by repeated sub-culturing cells across BG-11 solid agar plates and regular microscopic observation (Figure 1). The pure unialgal colonies were inoculated into a new 20 mL flask containing liquid BG-11 medium, then incubated in the abovementioned culture room.Table 1
BG-11 culture medium composition.
CompoundAmount (g/L)NaNO31.5K2HPO40.04MgSO4.7H2O0.075CaCl2.2H2O0.036Citric acid0.006Ferric ammonium citrate0.006EDTA0.001Na2CO30.002Trace-elements solution∗1mL/L∗H3BO3, 2.86 g/L; MnCl2.4H20, 1.81 g/L; ZnSO4.7H20, 0.222 g/L; Na2MoO4.2H2O, 0.39 g/L; CuSO4.5H2O, 0.079 g/L; Co (NO3)2 .6H2O, 0.0494 g/L.Figure 1
An agar plate streaked showing small isolated colonies arising from a single cell.


## 2.2. Morphological Identification

Preliminary identification of isolated microalgae was performed by direct samples microscopy of BG-11 culture (10μL) between 10 to 14 days. Among forty isolated microalgae (Figure 2), a total of twenty-nine single-cell strains were identified by evaluation of morphologic traits (cell shape, size, chloroplast, pyrenoid, flagellum, and cells arrangement) using an Olympus BX41 optical microscope equipped with an Olympus DP80 CCD camera, based on standard morphological feature keys [42, 43].Figure 2
Light microscopy photograph of all isolated microalgal strains using a 40X microspore lens.


## 2.3. Molecular Identification of Selected Microalgae

### 2.3.1. Primers and PCR Condition

A polymerase chain reaction (PCR) targeting 18S rRNA was used to identify the isolated microalgae. Genomic DNA was extracted using a DNP™ DNA extraction kit (SinaClon, Tehran, Iran) according to the manufacturer’s instructions and adjusted to the final concentration of 20 ng·μL−1.The partial 18S rRNA gene was PCR amplified using a universal primer set; 5′-GTCAGAGGTGAAATTCTTGGATTTA-3′ (Forward), and 5′-AGGGCAGGGACGTAATCAACG-3′ (Reverse) [40] in a total reaction volume of 80 μL containing: 40 μL Taq DNA polymerase master mix red 2X (Ampliqon, Odense, Denmark), 3 μL of each primer (20 μM), 6 μL template DNA and 28 μL sterile deionized distilled water. The PCR conditions were as follows: initial denaturation at 95°C for 5 min, followed by 30 cycles of denaturation at 95°C for 1 min, annealing at 58°C for 1 min, elongation at 72°C for 90 s, and the final extension step at 72°C for 90 s, according to the method used by Radha et al. [44] with slight modifications.


### 2.3.2. Electrophoresis and Gel Purification

The 1.8% Tris-Borate-EDTA (TBE) agarose gel and U: Genius3 GelDoc system (SYNGENE, Cambridge, United Kingdom) was used for electrophoresis and visualization of PCR products. The desired amplicons (∼750 bp) were extracted and purified using a Gel Purification Kit (Bioneer, Daejeon, South Korea) and sent for sequencing (Macrogen, Seoul, South Korea).


## 2.3.1. Primers and PCR Condition

A polymerase chain reaction (PCR) targeting 18S rRNA was used to identify the isolated microalgae. Genomic DNA was extracted using a DNP™ DNA extraction kit (SinaClon, Tehran, Iran) according to the manufacturer’s instructions and adjusted to the final concentration of 20 ng·μL−1.The partial 18S rRNA gene was PCR amplified using a universal primer set; 5′-GTCAGAGGTGAAATTCTTGGATTTA-3′ (Forward), and 5′-AGGGCAGGGACGTAATCAACG-3′ (Reverse) [40] in a total reaction volume of 80 μL containing: 40 μL Taq DNA polymerase master mix red 2X (Ampliqon, Odense, Denmark), 3 μL of each primer (20 μM), 6 μL template DNA and 28 μL sterile deionized distilled water. The PCR conditions were as follows: initial denaturation at 95°C for 5 min, followed by 30 cycles of denaturation at 95°C for 1 min, annealing at 58°C for 1 min, elongation at 72°C for 90 s, and the final extension step at 72°C for 90 s, according to the method used by Radha et al. [44] with slight modifications.


## 2.3.2. Electrophoresis and Gel Purification

The 1.8% Tris-Borate-EDTA (TBE) agarose gel and U: Genius3 GelDoc system (SYNGENE, Cambridge, United Kingdom) was used for electrophoresis and visualization of PCR products. The desired amplicons (∼750 bp) were extracted and purified using a Gel Purification Kit (Bioneer, Daejeon, South Korea) and sent for sequencing (Macrogen, Seoul, South Korea).


## 2.4. Analyzing Sequences

### 2.4.1. Multiple Sequence Alignment

The BioEdit (version 7.1.9) and Chromas Pro (version 2.1.3) programs were utilized to modify the obtained sequences. The final resulting 18S rRNA gene sequence was analyzed and used the BLASTn search tool for finding similarity searches against previously published 18S rRNA gene sequences in the NCBI databases. The representative sequences were selected and aligned using CLC Genomics workbench V20.0.


### 2.4.2. Molecular Phylogeny

A phylogenetic tree was constructed with the Neighbor-Joining algorithm using CLC genomics software version 20.0 with a bootstrap analysis of 500 replicates. The 18S rRNA gene sequence obtained in this study was deposited in GenBank databases.


## 2.4.1. Multiple Sequence Alignment

The BioEdit (version 7.1.9) and Chromas Pro (version 2.1.3) programs were utilized to modify the obtained sequences. The final resulting 18S rRNA gene sequence was analyzed and used the BLASTn search tool for finding similarity searches against previously published 18S rRNA gene sequences in the NCBI databases. The representative sequences were selected and aligned using CLC Genomics workbench V20.0.


## 2.4.2. Molecular Phylogeny

A phylogenetic tree was constructed with the Neighbor-Joining algorithm using CLC genomics software version 20.0 with a bootstrap analysis of 500 replicates. The 18S rRNA gene sequence obtained in this study was deposited in GenBank databases.


## 2.5. Screening of Isolated Microalgae for their Growth Potential

### 2.5.1. Cell Dry Weight Measurement

A total of twenty-nine isolated strains were inoculated in 20 mL flasks (10% (v/v)) containing 12 mL of BG-11 medium and were kept on a rotary shaker at 150 rpm and 20°С under a 12 : 12 h LD cycle of 100μE·m−2·S−1 illumination.The experiment was performed for 18 days. The sampling procedure was exploited every two days. The biomass concentrations were assessed gravimetrically. Biomass was harvested at 10000 rpm, for 5 min at 4°C and washed twice with distilled water followed by drying at 60°C in an oven for 48 h. Eventually, ten candidate microalgae were selected from the above step based on their high biomass productivity for further studies.


### 2.5.2. Growth Measurement and Kinetics

The biomass productivity (BP, gL−1·day−1), maximum specific growth rate (μmax day−1), and doubling time (h) were determined using equations (1), (2), and (3), respectively:(1)BP=X2−X1t2−t1,(2)μmax=lnX2−lnX1t2−t1,(3)td=ln2µmax,where X1 and X2 represent the dry cell weight (DCW, gL−1) and (t2 − t1) is required time for increasing cell concentration from X1 to X2 in the exponential phase [45]. All the experiments were done in triplicates to calculate standard deviations.


## 2.5.1. Cell Dry Weight Measurement

A total of twenty-nine isolated strains were inoculated in 20 mL flasks (10% (v/v)) containing 12 mL of BG-11 medium and were kept on a rotary shaker at 150 rpm and 20°С under a 12 : 12 h LD cycle of 100μE·m−2·S−1 illumination.The experiment was performed for 18 days. The sampling procedure was exploited every two days. The biomass concentrations were assessed gravimetrically. Biomass was harvested at 10000 rpm, for 5 min at 4°C and washed twice with distilled water followed by drying at 60°C in an oven for 48 h. Eventually, ten candidate microalgae were selected from the above step based on their high biomass productivity for further studies.


## 2.5.2. Growth Measurement and Kinetics

The biomass productivity (BP, gL−1·day−1), maximum specific growth rate (μmax day−1), and doubling time (h) were determined using equations (1), (2), and (3), respectively:(1)BP=X2−X1t2−t1,(2)μmax=lnX2−lnX1t2−t1,(3)td=ln2µmax,where X1 and X2 represent the dry cell weight (DCW, gL−1) and (t2 − t1) is required time for increasing cell concentration from X1 to X2 in the exponential phase [45]. All the experiments were done in triplicates to calculate standard deviations.


## 2.6. Biochemical Analysis of Selected Strains

### 2.6.1. Protein Content Analysis

Total protein contents were evaluated by Bradford method. A calibration curve was assessed using bovine serum albumin. Freeze-dried biomass of each sample (5 mg) that was harvested on the 18th day, was assayed for alkaline hydrolysis in 1 mL NaOH (1 N) for 1 h at 100°C using a hot plate (WiseTherm HB-R) and then centrifuged at 10000 rpm for 3 min. Protein content of each sample was calculated by measuring absorbance and comparing their absorbance at 595 nm with the bovine serum albumin standard curve [40]. This analysis was performed in triplicates and the average values were used for statistical evaluation of the results.


### 2.6.2. Pigments Analysis

(1) Determination of Chlorophyll a, Chlorophyll b, and Total Chlorophyll. Chlorophyll concentration determination was defined according to the method described by Eijckelhoff and Dekker [46]. Three mL of the microalgae suspension was taken and centrifuged at 2500 rpm for 10 min and rinsed with acetone (80%) and centrifuged again. The definition of the amount of chlorophyll extracted in the mentioned method was determined by colorimetric assay at 668.2 and 646.2 nm using an ELISA plate reader with the following equation: (4)ChlorophyllaμgmL=12.25×A668.2−2.79×A646.2.Five mL of the microalgae suspension was taken and added to the 3 mL acetone (80%), and the definition of the amount of total chlorophyll (a + b) extracted was specified using the following equation:(5)Total Chlorophyllaa+bμgmL=7.93×A664−19.53×A647.(2) β-Carotene Extraction and Assay. β-carotene extraction from the cell pellets was performed using the ethanol/n-hexane (2 : 1, v/v). The microalgal suspension (1 mL) was collected, centrifuged at 3000 rpm for 5 min, and rinsed with ethanol/n-hexane (2/1) and distilled water. β-carotene concentration is solved in the overlapping phase with n-hexane was determined by colorimetric assay at 450 nm using ELISA plate reader [47]. The definition of the amount of β-carotene extracted in n-hexane was specified using the following equation:(6)β−caroteneμgmL=25.2×A450.


## 2.6.1. Protein Content Analysis

Total protein contents were evaluated by Bradford method. A calibration curve was assessed using bovine serum albumin. Freeze-dried biomass of each sample (5 mg) that was harvested on the 18th day, was assayed for alkaline hydrolysis in 1 mL NaOH (1 N) for 1 h at 100°C using a hot plate (WiseTherm HB-R) and then centrifuged at 10000 rpm for 3 min. Protein content of each sample was calculated by measuring absorbance and comparing their absorbance at 595 nm with the bovine serum albumin standard curve [40]. This analysis was performed in triplicates and the average values were used for statistical evaluation of the results.


## 2.6.2. Pigments Analysis

(1) Determination of Chlorophyll a, Chlorophyll b, and Total Chlorophyll. Chlorophyll concentration determination was defined according to the method described by Eijckelhoff and Dekker [46]. Three mL of the microalgae suspension was taken and centrifuged at 2500 rpm for 10 min and rinsed with acetone (80%) and centrifuged again. The definition of the amount of chlorophyll extracted in the mentioned method was determined by colorimetric assay at 668.2 and 646.2 nm using an ELISA plate reader with the following equation: (4)ChlorophyllaμgmL=12.25×A668.2−2.79×A646.2.Five mL of the microalgae suspension was taken and added to the 3 mL acetone (80%), and the definition of the amount of total chlorophyll (a + b) extracted was specified using the following equation:(5)Total Chlorophyllaa+bμgmL=7.93×A664−19.53×A647.(2) β-Carotene Extraction and Assay. β-carotene extraction from the cell pellets was performed using the ethanol/n-hexane (2 : 1, v/v). The microalgal suspension (1 mL) was collected, centrifuged at 3000 rpm for 5 min, and rinsed with ethanol/n-hexane (2/1) and distilled water. β-carotene concentration is solved in the overlapping phase with n-hexane was determined by colorimetric assay at 450 nm using ELISA plate reader [47]. The definition of the amount of β-carotene extracted in n-hexane was specified using the following equation:(6)β−caroteneμgmL=25.2×A450.


## 2.7. CO2 Fixation Rate

By removing the carbon source from BG-11 media (Na2CO3), CO2 fixation rate for each selected strain was calculated in terms of CO2 consumption by each selected strain according to the following equation:(7)RCO2mgL−1d−1=P×Ccarbon×MCO2MC,where RCO2 is the rate of CO2 fixation (mg·L−1·d−1), P is the biomass productivity in 18 days (mg·L−1·d−1), MCO2 and MC represent the molecular weight of carbon dioxide and atomic weight of carbon, respectively. Ccarbon is the carbon content of the biomass was measured by a CHNS-O elemental analyzer equipped with gas chromatography and TCD detector (Costech, ECS 4010) [48].


## 2.8. FAME Analysis

### 2.8.1. Lipid Extraction

Fatty acid extraction was performed according to the Bligh and Dyer method with some modifications [49]. Briefly, lyophilized microalgal biomass (1 g) was boiled in 5 mL of isopropanol for 2 min and dried by vacuum rotary. A mixture of and 5 mL of chloroform was added, followed by 5 mL of potassium chloride solution (0.88%), to give a final solvent ratio of chloroform: methanol: water of 1 : 1 : 0.8. The mixture chloroform-methanol (1 : 2 v/v) was added to dried pellet in the presence of 0.015 g butylated hydroxytoluene (BHT) as an antioxidant agent, vortexed for 30 s and centrifuged (2000g, 5 min). To the above-obtained supernatant, 0.8 mL distilled water was shaken for 5 min using a separating funnel and allowed for 30 min to be separated. The total lipid was collected by concentrating the solvent phase under nitrogen gas.


### 2.8.2. Fatty Acid Esterification and GC/MS Analysis

The acid-catalyzed esterification was carried out in the Dien–Stark apparatus and was heated to reflux for 24 h at 75°C. The reaction mixture consisted of 0.5 g of extracted fatty acid dissolved in 3 mL methanol in the presence of 0.3 mL sulfuric acid. The mixture was washed with 4 mL saturated sodium hydrogen carbonate and then dried over anhydrous sodium sulfate. To assess an oily substance ready for GC/MS analysis, the solvent was removed under nitrogen gas [40].The GC/MS analysis was performed using an Agilent 7890A GC 5977B MSD equipped with an HP-5MS capillary column (30 m × 250μm × 0.25 μm, Agilent 19091S-433). The oven temperature was programmed between 70°C (5 min) to 270°C for 10 min and the rate was 7°C·min−1. The carrier gas was helium with a flow rate of 1 mL·min−1. The identification of fatty acid peaks was carried out by comparing the obtained mass spectra with the Wiley 7n.1 library.


## 2.8.1. Lipid Extraction

Fatty acid extraction was performed according to the Bligh and Dyer method with some modifications [49]. Briefly, lyophilized microalgal biomass (1 g) was boiled in 5 mL of isopropanol for 2 min and dried by vacuum rotary. A mixture of and 5 mL of chloroform was added, followed by 5 mL of potassium chloride solution (0.88%), to give a final solvent ratio of chloroform: methanol: water of 1 : 1 : 0.8. The mixture chloroform-methanol (1 : 2 v/v) was added to dried pellet in the presence of 0.015 g butylated hydroxytoluene (BHT) as an antioxidant agent, vortexed for 30 s and centrifuged (2000g, 5 min). To the above-obtained supernatant, 0.8 mL distilled water was shaken for 5 min using a separating funnel and allowed for 30 min to be separated. The total lipid was collected by concentrating the solvent phase under nitrogen gas.


## 2.8.2. Fatty Acid Esterification and GC/MS Analysis

The acid-catalyzed esterification was carried out in the Dien–Stark apparatus and was heated to reflux for 24 h at 75°C. The reaction mixture consisted of 0.5 g of extracted fatty acid dissolved in 3 mL methanol in the presence of 0.3 mL sulfuric acid. The mixture was washed with 4 mL saturated sodium hydrogen carbonate and then dried over anhydrous sodium sulfate. To assess an oily substance ready for GC/MS analysis, the solvent was removed under nitrogen gas [40].The GC/MS analysis was performed using an Agilent 7890A GC 5977B MSD equipped with an HP-5MS capillary column (30 m × 250μm × 0.25 μm, Agilent 19091S-433). The oven temperature was programmed between 70°C (5 min) to 270°C for 10 min and the rate was 7°C·min−1. The carrier gas was helium with a flow rate of 1 mL·min−1. The identification of fatty acid peaks was carried out by comparing the obtained mass spectra with the Wiley 7n.1 library.


## 3. Results and Discussion

### 3.1. Isolation and Identification of Microalgal Strains

A total of forty microalgal strains were isolated from thirteen samples collected along the southern coast of the Caspian Sea, Mazandaran province, and different habitats of Kohgiluyeh and Boyer-Ahmad provinces, the southwest of Iran. To obtain the biodiversity of mentioned regions in the south and north of the country, samples were collected from different habitats of these locations.Microscopic observation of these isolated microalgae indicated their colonial existence and purities which includes unicellular, colonial, and filamentous forms of Chlorophyceae and Cyanophyceae families. Photomicrograph of all isolated strains is presented in Figure2. The bio-diversity of strains is quite likely that reflects the vast adaptation ability of these strains in the freshwaters.All isolated microalgae were represented with specific codes and maintained in Microalgal Culture Collection of Shiraz University of Medical Sciences. Twenty-nine single-cell isolated strains were screened. Among them, ten strains were selected based on their high biomass productivity, and subjected to more characterizations. Preliminary morphology-based identification showed that the most morphological features of these algal isolates resemble the genus ofScenedesmus, which belongs to the family of Scenedesmaceae within the class of Chlorophyceae. Photomicrographs and morphological features of ten selected algal isolates are shown in supplementary data, Figure S1 and Table S1, respectively.Metzger and Largeau noted that for the same strain and within each chemical race, the morphological feature could vary regarding age and culture conditions [50]; therefore, PCR-amplification of 18S rRNA was used to confirm our morphological identification.


### 3.2. PCR Amplification and Phylogenetic Analysis

Molecular identification of ten selected microalgae, which displayed high productivity, was performed based on the partial 18S rRNA gene by comparing them with the reference sequences retrieved from the GenBank database. All of the selected strains belonged to the Chlorophyceae including, the generaTetradesmus, Desmodesmus, and Scenedesmus as the most abundant genus among them. Accession numbers of these sequences from GenBank and also the relationships among partial 18S rRNA that was inferred using phylogenic analysis by Neighbor-Joining algorithm are indicated in Table 2 and Figure 3(a), respectively. In parallel, the phylogenetic analysis was represented and compared with Acutodesmus sp., Coelastrum sp., and Tetradesmus sp. as similar species (Figure 3(b)).Table 2
Accession numbers for 18S rRNA gene sequence, strain, family, length in base pairs, and location for ten microalgal isolates.
CodeAccession numberStrainFamilyLength (bp)LocationMCCS∗1MG653564.1Scenedesmus sp. VN 001Chlorophyceae567Yasuj, 25 Km SEMCCS4MG653565.1Scenedesmus sp. VN 002Chlorophyceae596Yasuj, 30 Km SEMCCS17MG653566.1Scenedesmus sp. VN 003Chlorophyceae687Yasuj, 12 Km NWMCCS20MG653567.1Desmodesmus sp. VN 004Chlorophyceae671Yasuj, 25 Km NWMCCS21MG653568.1Scenedesmus sp. VN 005Chlorophyceae623Yasuj, 25 Km NWMCCS25MG653569.1Scenedesmus sp. VN 006Chlorophyceae701Yasuj, 30 Km SEMCCS28MG653570.1Desmodesmus sp. VN 007Chlorophyceae224Yasuj, Beshar river SEMCCS32MG653571.1Tetradesmus sp. VN 008Chlorophyceae613Yasuj, 25 Km NWMCCS35MG653572.1Scenedesmus sp. VN 009Chlorophyceae686Babol, 70 Km·SMCCS41MG653573.1Scenedesmus sp. VN 010Chlorophyceae688Yasuj, 125 Km SW∗Microalgal culture collection of Shiraz university of medical sciences (MCCS).Figure 3
(a) Phylogenic analysis of ten microalgal isolates with higher biomass productivities using MEGA-X analyzed by maximum likelihood (distance measure = Jukes-Cantor model, bootstrap = 500 replicates) with nearest neighbor interchange. (b) Phylogenic tree of three promising strains (Scenedesmus sp. VN009, Scenedesmus sp. VN006, and Tetradesmus sp. VN008) by compared with other closely related species using MEGA-X analyzed by maximum likelihood (distance measure = Jukes-Cantor model, bootstrap = 500 replicates) with nearest neighbor interchange.
(a)(b)


### 3.3. Kinetics and Growth Parameters

A total of twenty-nine single-cell strains were cultured in BG-11 under uniform conditions as described in Section2.5.1 section. During this period, the growth curve of one typical isolated strain was portrayed in Figure 4. All cultures were harvested at the end of the exponential phase of the typical isolate (18th day).Figure 4
Growth curve of typical isolated microalgal strain (MCCS 37) from Kohgiluyeh and Boyer-Ahmad provinces.Screening of the isolated microalgae has been performed primarily by growth characterization analyses, TableS2. After the screening program, ten strains were chosen where the biomass productivity was the main selection parameter and identified afterward (Table 3).Table 3
Top ten biomass productivity values of the isolated microalgae.
StrainSpecific growth rate (d−1)Doubling time (day)Biomass yield (gL−1)Biomass productivity (gL−1·d−1)Scenedesmus sp. VN 0090.0858.120.97 ± 0.100.054 ± 0.005Scenedesmus sp. VN 0060.0739.450.79 ± 0.090.044 ± 0.005Tetradesmus sp. VN 0080.0897.750.80 ± 0.100.044 ± 0.005Scenedesmus sp. VN 0100.0848.210.74 ± 0.080.041 ± 0.004Scenedesmus sp. VN 0010.0798.730.72 ± 0.030.040 ± 0.001Scenedesmus sp. VN 0020.088.620.74 ± 0.090.040 ± 0.005Scenedesmus sp. VN 0030.0967.190.70 ± 0.070.039 ± 0.004Scenedesmus sp. VN 0050.0858.120.68 ± 0.070.038 ± 0.004Desmodesmus sp. VN 0040.0729.580.63 ± 0.070.035 ± 0.004Desmodesmus sp. VN 0070.06410.780.58 ± 0.090.032 ± 0.005The observed biomass productivities were revealed to be higher than average values for previously reported concentrations, 0.03 gL−1·d−1 [11], 0.016 gL−1·d−1 [51]. The highest productivity among these selected strains belongs to Scenedesmus sp. VN009 (MCCS35), reaching 0.054 gL−1·d−1, followed by Scenedesmus sp. VN006 (MCCS25) and Tetradesmus sp. VN008 (MCCS32) at 0.044 gL−1·d−1. The lowest biomass productivity was found in two strains of Desmodesmus (Desmodesmus sp. VN007: 0.032 gL−1·d−1; Desmodesmus sp. VN004: 0.035 gL−1·d−1).These findings warrant the high potential of our selected strains to be exploited for scale-up studies for SCP production, or food and biofuel.


### 3.4. Protein and Pigment Content of the Selected Isolated Microalgae

The quantitative analysis of protein and pigment of the ten chosen strains is shown in Table4. The isolated Desmodesmus sp. VN 004 and Scenedesmus sp. VN 010 strains showed the highest protein content among the ten strains assayed (Desmodesmus sp. VN 004: 21.50% DW; Scenedesmus sp. VN 010: 19.30% DW). Scenedesmus sp. VN 006 protein content (18.58% DW) was similar to that of Scenedesmus sp. VN010 (Table 3). The least observed amount of protein was comprehended in Scenedesmus sp. VN 001 with 12.10% DW. It should be mentioned that the provided data in Table 4 only account for a small fraction of the weight percent as protein. The other identified parts were lipids, carbohydrates, nucleic acid fractions, other impurities, and possible errors.Table 4
Biomass composition analysis of selected microalgal strains.
StrainProtein (% w/w)Totalβ-carotene (% w/w)Totalβ-carotene (μg·mL−1)Chlorophylla(μ moles)Chlorophyllb(μ moles)Total chlorophyll (μ moles)Scenedesmus sp. VN 00112.10 ± 1.500.186 ± 0.131.77 ± 0.1460.08 ± 0.29269.13 ± 0.52343.82 ± 0.94Scenedesmus sp. VN 00214.35 ± 1.900.036 ± 0.090.35 ± 0.07237.53 ± 0.11172.27 ± 0.39423.1 ± 0.88Scenedesmus sp. VN 00316.40 ± 1.100.142 ± 0.121.21 ± 0.24139.25 ± 0.43247.60 ± 0.70401.95 ± 0.73Desmodesmus sp. VN 00421.50 ± 3.120.04 ± 0.010.35 ± 0.107.54 ± 0.11278.04 ± 0.81299.58 ± 0.89Scenedesmus sp. VN 00517.12 ± 0.510.052 ± 0.0100.45 ± 0.0613.2 ± 0.08231.62 ± 0.74256.62 ± 0.92Scenedesmus sp. VN 00618.58 ± 5.080.075 ± 0.0210.81 ± 0159.26 ± 0.29337.72 ± 0.97414.98 ± 1.07Desmodesmus sp. VN 00712.75 ± 1.360.047 ± 0.0050.4 ± 0.019.79 ± 0.37228.99 ± 0.73260.58 ± 0.99Tetradesmus sp. VN 00812.75 ± 3.00N. d∗N. d67.41 ± 0.22250.15 ± 0.66331.36 ± 1.21Scenedesmus sp. VN 00915.23 ± 2.120.057 ± 0.0090.71 ± 0.03165.43 ± 0.41331.51 ± 0.85516.74 ± 1.44Scenedesmus sp. VN 01019.30 ± 3.130.372 ± 0.0563.53 ± 0.08176.65 ± 0.25136.32 ± 0.63323.28 ± 1.01∗Not determined.The typical amount of protein in different microalgae during physiologic conditions has been reported to be about 10% to 71% DW [20]. However, the composition of culture medium and seasonal change might have influenced biomass content and compositions. Besides, the natural specifications of each strain are crucial as well.Theß-carotene content of all the ten selected microalgae indicated that the maximum value of ß-carotene content (0.372% DW) was achieved by Scenedesmus sp. VN 010, followed by Scenedesmus sp. VN 001 (0.186% DW) and Scenedesmus sp. VN 003 (0.142% DW). The ß-carotene content in Scenedesmus sp. VN 010 was found to be significantly more than the second highest strains (about twice the amount). ß-carotene possesses anti-oxidant activity; therefore, it plays an important role in human health. It is usually present in the range of 0.1–0.2% of the total dry weight of microalgae and up to 14% DW in industrialized species [52]. Different approaches of bioprocess and genetic engineering seem to be useful for enhancing the ß-carotene content of isolated microalgae for different industrial applications.Results from Table4 also show the highest chlorophyll content is owing to Scenedesmus sp. VN 009 (516.74 μmoles).


### 3.5. CO2 Fixation Rate

To avoid the role of the carbon source of BG-11 culture medium in CO2 fixation by selected microalgae, the main source of culture medium (Na2CO3) was removed, and the ability to stabilize CO2 in the air was investigated. Concerning the CO2 fixation rate, the data shown in Table 5 reveal that Scenedesmus sp. VN 009 had the best performance of 0.101 gL−1·d−1. The rate of carbon dioxide photosynthetically fixed by the rest of the strains was almost similar ranging from 0.067 to 0.078 gL−1·d−1, except for Desmodesmus sp. VN 007 that showed the lowest CO2 fixation rate of 0.042 gL−1·d−1. The maximum content of the CO2 fixation rate belonged to the strain, which was isolated from around Babol with −2 up to 50 m elevations, in contrast, those strains have been isolated from Yasuj with an average height of 1810 m height. It would seem that the strain isolated from the region with less height shows more CO2 fixation adaptability.Table 5
CO2 fixation rate was observed in each studied microalgal strain.
StrainCarbon concentration in biomassBiomass productivity (gL−1·d−1)CO2 fixation rate (gL−1·d−1)Scenedesmus sp. VN 00950.900.054 ± 0.0050.101 ± 0.009Scenedesmus sp. VN 00643.630.044 ± 0.0050.071 ± 0.008Tetradesmus sp. VN 00848.180.044 ± 0.0050.078 ± 0.008Scenedesmus sp. VN 01048.430.041 ± 0.0040.073 ± 0.007Scenedesmus sp. VN 00153.170.040 ± 0.0010.078 ± 0.002Scenedesmus sp. VN 00251.150.040 ± 0.0050.075 ± 0.009Scenedesmus sp. VN 00351.900.039 ± 0.0040.074 ± 0.008Scenedesmus sp. VN 00550.190.038 ± 0.0040.070 ± 0.007Desmodesmus sp. VN 00451.940.035 ± 0.0040.067 ± 0.008Desmodesmus sp. VN 00750.550.032 ± 0.0050.042 ± 0.009Previous reports indicate thatScenedesmus spp. Strains are drawn extensive attention for their appropriate CO2 fixation ability. Similarly, Isolated strains in this study have exhibited excellent performance in CO2 fixation, and the results are quite comparable to previous research, which have CO2 fixation rate in the range of 0.03 to 0.06 gL−1·d−1 with regard to their control conditions with 0.03% CO2 (v/v) [53–55].Considering the recent environmental issues such as global warming, industrialization, and population growth, having renewable energy and discovering potential strains for CO2 fixation are vitally important.According to obtained results for CO2 fixation, it could be suggested that isolated microalgae in this study are quite appropriate for CO2 fixation purposes.


### 3.6. Fatty Acid Composition

The Fatty acid profile gives vital information to choose the strain with the targeted application. Hence, Table6 is prepared to show the fatty acid profile of the top three isolated microalgae based on high biomass productivity values.Table 6
Fatty acid composition profile for microalgal lipids from three isolates exhibiting the highest biomass productivities.
Fatty acidsScenedesmus sp. VN 006Tetradesmus sp. VN 008Scenedesmus sp. VN 009Myristic acid (C14 : 0)—0.690.639Palmitic acid (C16 : 0)47.239.0521.13Palmitoleic acid (C16 : 1) (ω-7)—0.74—7,10-Hexadecadienoic acid (C16 : 2)—0.72—Stearic acid (C18 : 0)9.578.523.13Oleic acid (18 : 1) (ω-9)——5.48Linoleic acid (C18 : 2) (ω-6)8.0953.3220.7α-linolenic acid (C18 : 3) (ω-3)35.0926.9345.43γ-linolenic acid (C18 : 3) (ω-6)——1.21Arachidic acid (C20 : 0)——0.63Behenic acid (C22 : 0)——0.67Lignoceric acid (C24 : 0)——0.94PUFA43.1880.9767.34The results of GC/MS analysis comprehended that the three selected microalgal strains contain valuable saturated fatty acids (SFA). Palmitic acid (C 16 : 0), Stearic acid (C 18 : 0), and Arachidic acid (C 20 : 0) were regarded as the major identified SFAs. It can be observed that the algal lipid comprised high content of PUFA, ranging from 43% to 80%, which is desirable for pharmaceuticals and food industries.Linoleic acid (ω-6), α-linolenic acid (ω-3), and ɤ-linolenic acid (ω-6) (only in Scenedesmus sp. VN 009) are considered as the valuable PUFAs. Scenedesmus sp. VN 009 exhibited 45% ω-3 production, which was higher than previously reported [11, 56, 57]. This indicates a great promising candidate for ω-3 production. Also, Scenedesmus sp. VN 006 α-linolenic acid content was almost similar to that of Scenedesmus sp. VN009 (35% of total fatty acid). These polyunsaturated fatty acids have been known as essential precursors of longer n-3 PUFAs, especially DHA, that eventually could possess several beneficial impacts on human health, from influencing fetal growth and development to reducing body fat [58], prevention of cardiac arrhythmias, sudden cardiac death [59], Alzheimer’s disease [33], alleviate coronary heart disease [60], and also diabetes [61].Furthermore, there are some studies indicate thatα-linolenic acid (ALA) and linoleic acid have exerted therapeutic effects as they could reduce ischemic brain damage, and alleviate functional stroke recovery and psychiatric disorders [62]. Furthermore, details are presented in Table 6.


## 3.1. Isolation and Identification of Microalgal Strains

A total of forty microalgal strains were isolated from thirteen samples collected along the southern coast of the Caspian Sea, Mazandaran province, and different habitats of Kohgiluyeh and Boyer-Ahmad provinces, the southwest of Iran. To obtain the biodiversity of mentioned regions in the south and north of the country, samples were collected from different habitats of these locations.Microscopic observation of these isolated microalgae indicated their colonial existence and purities which includes unicellular, colonial, and filamentous forms of Chlorophyceae and Cyanophyceae families. Photomicrograph of all isolated strains is presented in Figure2. The bio-diversity of strains is quite likely that reflects the vast adaptation ability of these strains in the freshwaters.All isolated microalgae were represented with specific codes and maintained in Microalgal Culture Collection of Shiraz University of Medical Sciences. Twenty-nine single-cell isolated strains were screened. Among them, ten strains were selected based on their high biomass productivity, and subjected to more characterizations. Preliminary morphology-based identification showed that the most morphological features of these algal isolates resemble the genus ofScenedesmus, which belongs to the family of Scenedesmaceae within the class of Chlorophyceae. Photomicrographs and morphological features of ten selected algal isolates are shown in supplementary data, Figure S1 and Table S1, respectively.Metzger and Largeau noted that for the same strain and within each chemical race, the morphological feature could vary regarding age and culture conditions [50]; therefore, PCR-amplification of 18S rRNA was used to confirm our morphological identification.


## 3.2. PCR Amplification and Phylogenetic Analysis

Molecular identification of ten selected microalgae, which displayed high productivity, was performed based on the partial 18S rRNA gene by comparing them with the reference sequences retrieved from the GenBank database. All of the selected strains belonged to the Chlorophyceae including, the generaTetradesmus, Desmodesmus, and Scenedesmus as the most abundant genus among them. Accession numbers of these sequences from GenBank and also the relationships among partial 18S rRNA that was inferred using phylogenic analysis by Neighbor-Joining algorithm are indicated in Table 2 and Figure 3(a), respectively. In parallel, the phylogenetic analysis was represented and compared with Acutodesmus sp., Coelastrum sp., and Tetradesmus sp. as similar species (Figure 3(b)).Table 2
Accession numbers for 18S rRNA gene sequence, strain, family, length in base pairs, and location for ten microalgal isolates.
CodeAccession numberStrainFamilyLength (bp)LocationMCCS∗1MG653564.1Scenedesmus sp. VN 001Chlorophyceae567Yasuj, 25 Km SEMCCS4MG653565.1Scenedesmus sp. VN 002Chlorophyceae596Yasuj, 30 Km SEMCCS17MG653566.1Scenedesmus sp. VN 003Chlorophyceae687Yasuj, 12 Km NWMCCS20MG653567.1Desmodesmus sp. VN 004Chlorophyceae671Yasuj, 25 Km NWMCCS21MG653568.1Scenedesmus sp. VN 005Chlorophyceae623Yasuj, 25 Km NWMCCS25MG653569.1Scenedesmus sp. VN 006Chlorophyceae701Yasuj, 30 Km SEMCCS28MG653570.1Desmodesmus sp. VN 007Chlorophyceae224Yasuj, Beshar river SEMCCS32MG653571.1Tetradesmus sp. VN 008Chlorophyceae613Yasuj, 25 Km NWMCCS35MG653572.1Scenedesmus sp. VN 009Chlorophyceae686Babol, 70 Km·SMCCS41MG653573.1Scenedesmus sp. VN 010Chlorophyceae688Yasuj, 125 Km SW∗Microalgal culture collection of Shiraz university of medical sciences (MCCS).Figure 3
(a) Phylogenic analysis of ten microalgal isolates with higher biomass productivities using MEGA-X analyzed by maximum likelihood (distance measure = Jukes-Cantor model, bootstrap = 500 replicates) with nearest neighbor interchange. (b) Phylogenic tree of three promising strains (Scenedesmus sp. VN009, Scenedesmus sp. VN006, and Tetradesmus sp. VN008) by compared with other closely related species using MEGA-X analyzed by maximum likelihood (distance measure = Jukes-Cantor model, bootstrap = 500 replicates) with nearest neighbor interchange.
(a)(b)


## 3.3. Kinetics and Growth Parameters

A total of twenty-nine single-cell strains were cultured in BG-11 under uniform conditions as described in Section2.5.1 section. During this period, the growth curve of one typical isolated strain was portrayed in Figure 4. All cultures were harvested at the end of the exponential phase of the typical isolate (18th day).Figure 4
Growth curve of typical isolated microalgal strain (MCCS 37) from Kohgiluyeh and Boyer-Ahmad provinces.Screening of the isolated microalgae has been performed primarily by growth characterization analyses, TableS2. After the screening program, ten strains were chosen where the biomass productivity was the main selection parameter and identified afterward (Table 3).Table 3
Top ten biomass productivity values of the isolated microalgae.
StrainSpecific growth rate (d−1)Doubling time (day)Biomass yield (gL−1)Biomass productivity (gL−1·d−1)Scenedesmus sp. VN 0090.0858.120.97 ± 0.100.054 ± 0.005Scenedesmus sp. VN 0060.0739.450.79 ± 0.090.044 ± 0.005Tetradesmus sp. VN 0080.0897.750.80 ± 0.100.044 ± 0.005Scenedesmus sp. VN 0100.0848.210.74 ± 0.080.041 ± 0.004Scenedesmus sp. VN 0010.0798.730.72 ± 0.030.040 ± 0.001Scenedesmus sp. VN 0020.088.620.74 ± 0.090.040 ± 0.005Scenedesmus sp. VN 0030.0967.190.70 ± 0.070.039 ± 0.004Scenedesmus sp. VN 0050.0858.120.68 ± 0.070.038 ± 0.004Desmodesmus sp. VN 0040.0729.580.63 ± 0.070.035 ± 0.004Desmodesmus sp. VN 0070.06410.780.58 ± 0.090.032 ± 0.005The observed biomass productivities were revealed to be higher than average values for previously reported concentrations, 0.03 gL−1·d−1 [11], 0.016 gL−1·d−1 [51]. The highest productivity among these selected strains belongs to Scenedesmus sp. VN009 (MCCS35), reaching 0.054 gL−1·d−1, followed by Scenedesmus sp. VN006 (MCCS25) and Tetradesmus sp. VN008 (MCCS32) at 0.044 gL−1·d−1. The lowest biomass productivity was found in two strains of Desmodesmus (Desmodesmus sp. VN007: 0.032 gL−1·d−1; Desmodesmus sp. VN004: 0.035 gL−1·d−1).These findings warrant the high potential of our selected strains to be exploited for scale-up studies for SCP production, or food and biofuel.


## 3.4. Protein and Pigment Content of the Selected Isolated Microalgae

The quantitative analysis of protein and pigment of the ten chosen strains is shown in Table4. The isolated Desmodesmus sp. VN 004 and Scenedesmus sp. VN 010 strains showed the highest protein content among the ten strains assayed (Desmodesmus sp. VN 004: 21.50% DW; Scenedesmus sp. VN 010: 19.30% DW). Scenedesmus sp. VN 006 protein content (18.58% DW) was similar to that of Scenedesmus sp. VN010 (Table 3). The least observed amount of protein was comprehended in Scenedesmus sp. VN 001 with 12.10% DW. It should be mentioned that the provided data in Table 4 only account for a small fraction of the weight percent as protein. The other identified parts were lipids, carbohydrates, nucleic acid fractions, other impurities, and possible errors.Table 4
Biomass composition analysis of selected microalgal strains.
StrainProtein (% w/w)Totalβ-carotene (% w/w)Totalβ-carotene (μg·mL−1)Chlorophylla(μ moles)Chlorophyllb(μ moles)Total chlorophyll (μ moles)Scenedesmus sp. VN 00112.10 ± 1.500.186 ± 0.131.77 ± 0.1460.08 ± 0.29269.13 ± 0.52343.82 ± 0.94Scenedesmus sp. VN 00214.35 ± 1.900.036 ± 0.090.35 ± 0.07237.53 ± 0.11172.27 ± 0.39423.1 ± 0.88Scenedesmus sp. VN 00316.40 ± 1.100.142 ± 0.121.21 ± 0.24139.25 ± 0.43247.60 ± 0.70401.95 ± 0.73Desmodesmus sp. VN 00421.50 ± 3.120.04 ± 0.010.35 ± 0.107.54 ± 0.11278.04 ± 0.81299.58 ± 0.89Scenedesmus sp. VN 00517.12 ± 0.510.052 ± 0.0100.45 ± 0.0613.2 ± 0.08231.62 ± 0.74256.62 ± 0.92Scenedesmus sp. VN 00618.58 ± 5.080.075 ± 0.0210.81 ± 0159.26 ± 0.29337.72 ± 0.97414.98 ± 1.07Desmodesmus sp. VN 00712.75 ± 1.360.047 ± 0.0050.4 ± 0.019.79 ± 0.37228.99 ± 0.73260.58 ± 0.99Tetradesmus sp. VN 00812.75 ± 3.00N. d∗N. d67.41 ± 0.22250.15 ± 0.66331.36 ± 1.21Scenedesmus sp. VN 00915.23 ± 2.120.057 ± 0.0090.71 ± 0.03165.43 ± 0.41331.51 ± 0.85516.74 ± 1.44Scenedesmus sp. VN 01019.30 ± 3.130.372 ± 0.0563.53 ± 0.08176.65 ± 0.25136.32 ± 0.63323.28 ± 1.01∗Not determined.The typical amount of protein in different microalgae during physiologic conditions has been reported to be about 10% to 71% DW [20]. However, the composition of culture medium and seasonal change might have influenced biomass content and compositions. Besides, the natural specifications of each strain are crucial as well.Theß-carotene content of all the ten selected microalgae indicated that the maximum value of ß-carotene content (0.372% DW) was achieved by Scenedesmus sp. VN 010, followed by Scenedesmus sp. VN 001 (0.186% DW) and Scenedesmus sp. VN 003 (0.142% DW). The ß-carotene content in Scenedesmus sp. VN 010 was found to be significantly more than the second highest strains (about twice the amount). ß-carotene possesses anti-oxidant activity; therefore, it plays an important role in human health. It is usually present in the range of 0.1–0.2% of the total dry weight of microalgae and up to 14% DW in industrialized species [52]. Different approaches of bioprocess and genetic engineering seem to be useful for enhancing the ß-carotene content of isolated microalgae for different industrial applications.Results from Table4 also show the highest chlorophyll content is owing to Scenedesmus sp. VN 009 (516.74 μmoles).


## 3.5. CO2 Fixation Rate

To avoid the role of the carbon source of BG-11 culture medium in CO2 fixation by selected microalgae, the main source of culture medium (Na2CO3) was removed, and the ability to stabilize CO2 in the air was investigated. Concerning the CO2 fixation rate, the data shown in Table 5 reveal that Scenedesmus sp. VN 009 had the best performance of 0.101 gL−1·d−1. The rate of carbon dioxide photosynthetically fixed by the rest of the strains was almost similar ranging from 0.067 to 0.078 gL−1·d−1, except for Desmodesmus sp. VN 007 that showed the lowest CO2 fixation rate of 0.042 gL−1·d−1. The maximum content of the CO2 fixation rate belonged to the strain, which was isolated from around Babol with −2 up to 50 m elevations, in contrast, those strains have been isolated from Yasuj with an average height of 1810 m height. It would seem that the strain isolated from the region with less height shows more CO2 fixation adaptability.Table 5
CO2 fixation rate was observed in each studied microalgal strain.
StrainCarbon concentration in biomassBiomass productivity (gL−1·d−1)CO2 fixation rate (gL−1·d−1)Scenedesmus sp. VN 00950.900.054 ± 0.0050.101 ± 0.009Scenedesmus sp. VN 00643.630.044 ± 0.0050.071 ± 0.008Tetradesmus sp. VN 00848.180.044 ± 0.0050.078 ± 0.008Scenedesmus sp. VN 01048.430.041 ± 0.0040.073 ± 0.007Scenedesmus sp. VN 00153.170.040 ± 0.0010.078 ± 0.002Scenedesmus sp. VN 00251.150.040 ± 0.0050.075 ± 0.009Scenedesmus sp. VN 00351.900.039 ± 0.0040.074 ± 0.008Scenedesmus sp. VN 00550.190.038 ± 0.0040.070 ± 0.007Desmodesmus sp. VN 00451.940.035 ± 0.0040.067 ± 0.008Desmodesmus sp. VN 00750.550.032 ± 0.0050.042 ± 0.009Previous reports indicate thatScenedesmus spp. Strains are drawn extensive attention for their appropriate CO2 fixation ability. Similarly, Isolated strains in this study have exhibited excellent performance in CO2 fixation, and the results are quite comparable to previous research, which have CO2 fixation rate in the range of 0.03 to 0.06 gL−1·d−1 with regard to their control conditions with 0.03% CO2 (v/v) [53–55].Considering the recent environmental issues such as global warming, industrialization, and population growth, having renewable energy and discovering potential strains for CO2 fixation are vitally important.According to obtained results for CO2 fixation, it could be suggested that isolated microalgae in this study are quite appropriate for CO2 fixation purposes.


## 3.6. Fatty Acid Composition

The Fatty acid profile gives vital information to choose the strain with the targeted application. Hence, Table6 is prepared to show the fatty acid profile of the top three isolated microalgae based on high biomass productivity values.Table 6
Fatty acid composition profile for microalgal lipids from three isolates exhibiting the highest biomass productivities.
Fatty acidsScenedesmus sp. VN 006Tetradesmus sp. VN 008Scenedesmus sp. VN 009Myristic acid (C14 : 0)—0.690.639Palmitic acid (C16 : 0)47.239.0521.13Palmitoleic acid (C16 : 1) (ω-7)—0.74—7,10-Hexadecadienoic acid (C16 : 2)—0.72—Stearic acid (C18 : 0)9.578.523.13Oleic acid (18 : 1) (ω-9)——5.48Linoleic acid (C18 : 2) (ω-6)8.0953.3220.7α-linolenic acid (C18 : 3) (ω-3)35.0926.9345.43γ-linolenic acid (C18 : 3) (ω-6)——1.21Arachidic acid (C20 : 0)——0.63Behenic acid (C22 : 0)——0.67Lignoceric acid (C24 : 0)——0.94PUFA43.1880.9767.34The results of GC/MS analysis comprehended that the three selected microalgal strains contain valuable saturated fatty acids (SFA). Palmitic acid (C 16 : 0), Stearic acid (C 18 : 0), and Arachidic acid (C 20 : 0) were regarded as the major identified SFAs. It can be observed that the algal lipid comprised high content of PUFA, ranging from 43% to 80%, which is desirable for pharmaceuticals and food industries.Linoleic acid (ω-6), α-linolenic acid (ω-3), and ɤ-linolenic acid (ω-6) (only in Scenedesmus sp. VN 009) are considered as the valuable PUFAs. Scenedesmus sp. VN 009 exhibited 45% ω-3 production, which was higher than previously reported [11, 56, 57]. This indicates a great promising candidate for ω-3 production. Also, Scenedesmus sp. VN 006 α-linolenic acid content was almost similar to that of Scenedesmus sp. VN009 (35% of total fatty acid). These polyunsaturated fatty acids have been known as essential precursors of longer n-3 PUFAs, especially DHA, that eventually could possess several beneficial impacts on human health, from influencing fetal growth and development to reducing body fat [58], prevention of cardiac arrhythmias, sudden cardiac death [59], Alzheimer’s disease [33], alleviate coronary heart disease [60], and also diabetes [61].Furthermore, there are some studies indicate thatα-linolenic acid (ALA) and linoleic acid have exerted therapeutic effects as they could reduce ischemic brain damage, and alleviate functional stroke recovery and psychiatric disorders [62]. Furthermore, details are presented in Table 6.


## 4. Conclusions

To sum up, we have isolated twenty-nine microalgal strains from different locations in Iran, selected strains were further characterized to find the most promising strains for food along with CO2 fixation ability. Scenedesmus sp. and Desmodesmus sp. were the most abundant microalgae in the Kohgiluyeh and Boyer-Ahmad provinces (southwest of Iran). Some strains showed high biomass productivity and CO2 fixation rate. Fatty acid analyzed in this study shows the ability of Scenedesmus sp. VN009 and Scenedesmus sp. VN 006 to accumulate a high level of α-linolenic acid (ω-3). As far as we know, the observed amounts for ω-3 productions were higher than any other available reports in the literature, thereby these strains could be targeted for large-scale production.



---

*Source: 4386268-2022-08-12.xml*