ORIGINAL RESEARCH article

Front. Plant Sci., 29 November 2022

Sec. Crop and Product Physiology

Volume 13 - 2022 | https://doi.org/10.3389/fpls.2022.1034238

DNA barcoding and biomass accumulation rates of native Iranian duckweed species for biotechnological applications

  • 1. National Institute of Genetic Engineering and Biotechnology (NIGEB), Department of Agricultural Biotechnology, Tehran, Iran

  • 2. University of Greifswald, Institute of Botany and Landscape Ecology, Greifswald, Germany

Article metrics

View details

5

Citations

2,8k

Views

1,1k

Downloads

Abstract

The Lemnaceae family (duckweed) consists of at least three recognized genera with six reported species in Iran that are distributed in wetlands. Duckweeds are the simplest and smallest flowering aquatic monocots with free-floating fronds that can reproduce asexually every 2–3 days. Duckweed could be a major source of balanced amino acids and high protein content, which is increasingly promising for biotechnological applications. For molecular classification and species identification of the collected samples, DNA barcoding was performed using two standard chloroplast markers, the spacer region between the ATP synthase subunits F and H (atpF-atpH) and the intron region of the ribosomal protein S16 (rps16). The results confirm the presence of four species belonging to the two genera Lemna and Spirodela. In addition, L. turionifera was detected for the first time in Iran. Due to the high growth rates of duckweed, measurement of biomass accumulation and doubling time are important factors in determining growth potential, especially for native species. The relative growth rates (RGR), doubling times (DT), biomass accumulation, and relative weekly yields (RY) of 40 distinct duckweed clones were determined under standard cultivation conditions. The dry weight–based RGR ranged from 0.149 to more than 0.600 per day, DT from 1.12 to 9 days, and RY from 7 to 108.9 per week. All values are comparable with previous studies. RGR and RY of selected clones are higher than the growth potential for a wide range of wild plants and common crops. These data support that native duckweed has high productivity value and should be further investigated as a potentially rich protein source for alternative human food, livestock feed, and recombinant protein production.

Introduction

The cosmopolitan monocotyledonous family Lemnaceae (duckweed), which includes the smallest and fastest growing angiosperms known, comprises 36 species in five genera: Spirodela Schleid., Landoltia Les & Crawford, Lemna L., Wolffia Horkel ex Schleid., and Wolffiella Hegelm. (Bog et al., 2020). However, some authors claim that duckweeds should be reclassified as subfamily Lemnoideae in the Araceae family based on a close phylogenetic relationship. However, the inclusion of duckweed in the Araceae would remove a useful and well-defined taxonomic category of an angiosperm family that has been used by duckweed biologists for many years (Tippery et al., 2021). These tiny aquatic plants grow on or below the surface of slow-flowing, nutrient-enriched water bodies. Their morphology is highly reduced to simple leaf-like structures also known as fronds that appear genus-specific with or without roots. Duckweeds can create genetically uniform populations by their rapid vegetative propagation. These properties, their small size, rapid and high yield growth, and additionally relatively small genome sizes make duckweeds an ideal experimental material and a powerful platform for various biotechnological applications (Appenroth et al., 2013; Heenatigala et al., 2020; Tippery et al., 2021). Under optimized growth conditions, duckweed contains a high protein content of up to 45% with high-quality and easily digestible amino acids close to the recommendations of the World Health Organization (WHO), which is an important nutritional index (Mes et al., 2022a; Pagliuso et al., 2022). Therefore, there is increased interest in the use of duckweed species from the genera Wolffia and Lemna as a good protein source, particularly for use in human food and animal nutrition (Edelman and Colt, 2016; Appenroth et al., 2017; Appenroth et al., 2018). Moreover, the high fiber content (∼25% of dry weight) and polyunsaturated fatty acids (more than 60% of total fat) are shown to be a unique nutrient composition of duckweed (Appenroth et al., 2018). While major crops such as rice, maize, and wheat often have an imbalance in nutrient composition, the amino acids, vitamins, mineral profiles, and fatty acid fractions of duckweed have a high-quality composition (Edelman and Colt, 2016). The increase in world population, climate change, and decrease in food supply have increased pressure on food systems. In addition, excessive land use and agricultural activities have led to soil erosion, resulting in a 0.4% per year decline in global crop yields (Pennock, 2019). Therefore, there are increasing demands for the development of a sustainable food and feed safety system (Pagliuso et al., 2022). For example, the European Food and Safety Authority considers all genera of duckweed as novel foods (Mes et al., 2022b). Most notably, the ease of cultivation of these tiny aquatic plants in multilayered vertical farming systems and their high tolerance to a wide range of environmental conditions around the world may reduce competition with terrestrial crops (Coughlan et al., 2022; Pagliuso et al., 2022).

In addition, duckweed can be used in phytoremediation, water quality measurement, and wastewater treatment. Duckweeds have the ability to accumulate macronutrients and micronutrients hundreds of times compared with the mineral concentration of the water in which they proliferate (Chakrabarti et al., 2018; Walsh et al., 2021). The high biomass accumulation with a typically high protein content between 24% and 45% makes duckweed a suitable supplement for animal feed (Iatrou et al., 2015; Pagliuso et al., 2022). Further, duckweed can be used for bioethanol production due to its high starch accumulation under stress conditions (Sree and Appenroth, 2014; Liu et al., 2019).

The potential application of duckweed in plant bioreactors has attracted increasing attention due to its rapid doubling time (DT) and proliferation of uniform clones with nearly exponential growth (Edelman et al., 2020). Their growth rate is nearly 28 times faster than conventional crops used for human nutrition (Pagliuso et al., 2022; Coughlan et al., 2022). Undoubtedly, this higher reproductive rate will significantly shorten the production cycle of duckweed in bioreactors, leading to a maximum biomass accumulation of up to 100 tons of dry matter per hectare per year (Cao et al., 2018). Therefore, the study of biomass production and growth factors in duckweed under different cultivation conditions could be interesting. There is a number of pilot studies on duckweed biomass production under environmental conditions where wastewater or enriched medium is used to grow duckweed. The reports describe high-yield biomass production of 8 t dw/ha/y for Wolffia arrhiza (L.) Horkel ex Wimm. (Fujita et al., 1999) and 36 t dw/ha/y for Spirodela polyrhiza (L.) Schleid. (Xu et al., 2012) and up to 104 t/ha/year for Lemna minor L. (Frederic et al., 2006), which is comparable to the average yields of major land crops reported by the Food and Agriculture Organization of the United Nations (FAO, 2013) and the U.S. Department of Agriculture (USDA) (Ziegler et al., 2015).

The quality value of duckweed as a food source depends on the content and composition of constituents, particularly amino acid profiles, protein content, and high potential for rapid growth. Therefore, cultivation conditions and especially the genetic background of ecotypes are assumed to play a crucial role (Appenroth et al., 2018; Chakrabarti et al., 2018). Thus, assuming that different geographic isolates (ecotypes) of duckweeds have genetically differentiated due to adaptation to specific environmental conditions, the resulting clones may exhibit different physiological and growth behavior (Ziegler et al., 2015; Chakrabarti et al., 2018; Walsh et al., 2021). Studying a wide range of these ecotypes around the world may lead to the identification of superior clones in terms of growth characteristics that can be eligible for other studies in different fields, such as biochemical analysis to introduce an alternative food supply.

To date, only a limited number of extensive studies have been conducted to investigate growth factors and biomass production in duckweed ecotypes from different parts of the world. Bergmann et al. (2000) studied 41 geographic isolates of 12 species from Landolt’s worldwide stock collection at ETH Zurich (now hosted by the Istituto di Biologia e Biotecnologia Agraria in Milano, Italy) under in vitro conditions in a synthetic medium. To assess growth, they reported only wet weight gain and the percentage dry weight during the 11-day growth period for selection of superior geographic isolates (Bergmann et al., 2000). Among others, the relative growth rate (RGR) is an important growth factor that reflects the growth potential, especially in duckweeds. A comprehensive study of the RGR of duckweed was already carried out by Landolt (1957), who investigated 71 clones of 13 species. Ziegler et al. (2015) presented two other growth factors in addition to RGR to provide comparable data with other reports, e.g., on terrestrial crops. RGR, DT, and relative weekly yield (RY) of 39 ecotypes from 13 duckweed species were determined under standard cultivation conditions using a modified Schenk–Hildebrand medium for 7 days. Here, the mean RGR was 0.304 per day for Spirodela and 0.396 per day for Lemna. In general, RGR ranged from 0.153 to 0.519 per day, DT from 1.34 to 4.54 days, and RY from 2.9 to 37.8 per week for the duckweed species studied (Ziegler et al., 2015). Sree et al. (2015) investigated the RGR of 25 clones representing all 11 species of the genus Wolffia, the genus commonly used for human nutrition. They present a clone of Wolffia microscopica (Griff.) Kurz with a doubling time of 29.3 h, which is the fastest growing flowering plant (Sree et al., 2015). Other reports determining the RGR of duckweed species have reported an RGR of 0.31 per day for Lemna minor and 0.30 and 0.42 per day for Lemna gibba (Lasfar et al., 2007).

One of the first important steps is the precise identification of plant species. Several methods, such as morphological, biochemical, and molecular comparisons, can be used for correct identification. Among the mentioned marker types, comparison of molecular data with references is the most reliable (Dogan et al., 2014; Hasanbegovic et al., 2021; Saran et al., 2021). Specifically for duckweed, identification using only morphological characteristics is nearly impossible even for experts because the morphological structure of duckweed is greatly reduced. Using molecular methods, such as DNA barcoding, which is based on DNA markers, it is possible to reproducibly and reliably identify most duckweed species (Borisjuk et al., 2015; Bog et al., 2019). To the best of our knowledge, the native duckweeds of Iran have not been studied at the molecular level or in terms of growth rate. This is the first report on DNA barcoding of the only duckweed collection in Iran.

In the present study, RGR, DT, and RY for duckweed species native to Iran are investigated for the first time. Most of the clones studied are from the north of Iran, where most of the duckweed habitats are located. The investigation of 40 Iranian clones, which can be assigned to four species within the two genera Lemna and Spirodela, aims to determine accumulation of biomass yields and to identify the superior ecotypes with high growth potential under laboratory conditions for future biotechnological applications and food safety research.

Material and methods

Plant material

Duckweed samples were collected from different natural ponds in the north of Iran (Mazandaran and Gilan provinces) and a region in the west of Iran (Kermanshah province). The geographical distribution map of sampling can be found in Supplementary Figure 1. Fronds were rinsed in clean tap water and sterilized using 2.5% sodium hypochlorite solution for 1 min and subsequently washed three times with sterile distilled water. Sterilized fronds of all Lemna species were then cultivated in modified Hoagland medium (Khvatkov et al., 2019) except L. gibba, which was cultivated in NF medium (Muranaka et al., 2015). Schenk–Hildebrand medium (Schenk and Hildebrandt, 1972) was used for the Spirodela species. All media were supplemented with 1% sucrose. All clones were cultivated under standard cultivation conditions (ISO 20079, 2005) with a 16/8 h light/dark photoperiod with 80 μmol/m2/s light intensity from fluorescent light tubes (40W) (Pars shahab, Tehran, Iran) at 25°C.

Morphological identification and molecular analysis

Duckweed samples were primarily identified morphologically using the key from Riedl (1976) and the updated key from Bog et al. (2020) based on frond shape, frond size, and number of roots and veins. The identity of the clones that were chosen for the biomass accumulation test was confirmed by DNA barcoding. For molecular analysis, total DNA was extracted by a modified CTAB protocol (Murray and Thompson, 1980). The chloroplast marker from the noncoding spacer atpF-atpH was amplified using the primers atpF-atpH forward (5’ ACTCGCACACACTCCCTTTCC 3’) and atpF-atpH reverse (5’ GCTTTTATGGAAGCTTTAACAAT 3’) as described previously (Wang et al., 2010). The PCR conditions were predenaturation at 94°C for 2 min, followed by 35 cycles of 94°C, 15 s; 51°C, 15 s; 72°C, 40 s; and a final extension at 72°C for 5 min. The primer set of the second marker, rps16, was used to amplify the chloroplast ribosomal protein S16 gene intron with the degenerate primers rps16 F (5’ AAACGATGTGGTARAAAGCAAC 3’) and rps16 R (5’ AACATCWATTGCAASGATTCGATA 3’) as described previously (Shaw et al., 2005). The PCR conditions were predenaturation at 94°C for 5 min, followed by 45 cycles at 94°C, 30 s; 61°C, 50 s; 72°C, 80s; and a final extension at 72°C for 7 min. The PCR fragments were purified and further processed for sequencing by the Beijing Genomic Institute (BGI, Shenzhen, China). Sequences were deposited in GenBank (https://www.ncbi.nlm.nih.gov/). Accession numbers of the rps16 and atpF-atpH sequences are listed in Table 1.

Table 1

RowSpeciesStrainOriginAccession number
atpF-atpHrps 16
1Lemna minor1CTonekabon pond, IranMT891062MW308182
2Lemna minor2CTonekabon pond, IranMT891063MZ422535
3Lemna minor3aTonekabon pond, IranMT891064MW308183
4Lemna minor4BMMansoori pond, IranMT891065MZ422534
5Lemna minor5CTonekabon pond, IranMT891066MW308184
6Lemna minor6aTonekabon pond, IranMT891067MW308185
7Lemna minor7WTonekabon pond, IranMT891068MW308186
8Lemna minor8CTonekabon pond, IranMT891069MW308187
9Lemna minor9aTonekabon pond, IranMT891070MW308188
10Lemna minor10CTonekabon pond, IranMT891071MW308189
11Lemna minor11CTonekabon pond, IranMT891072MW308190
12Lemna minor12WTonekabon pond, IranMT891073MZ422533
13Lemna minor13CTonekabon pond, IranMT891074MW308191
14Lemna minor14BMMansoori pond, IranMT891075MZ422532
15Lemna minor15CTonekabon pond, IranMT891076MW308192
16Lemna minor16CTonekabon pond, IranMT891077MW308193
17Lemna turionifera17AMAmir kelaye international lagoon, IranMT891086MW308200
18Lemna minor18CTonekabon pond, IranMT891078MW308194
19Lemna minor19CTonekabon pond, IranMT891079MW308195
20Lemna turionifera20AMAmir kelaye international lagoon, IranMT891087MW308201
21Lemna minor21BMMansoori pond, IranMT891080MW308196
22Lemna minor22aTonekabon pond, IranMT891081MW308197
23Lemna minor23CTonekabon pond, IranMT891082MW308198
24Lemna minor24WTonekabon pond, IranMT891083MZ422531
25Lemna minor25CTonekabon pond, IranMT891084MZ422530
26Lemna minor26WTonekabon pond, IranMT891085MW308199
27Lemna gibba1LALangarud paddy,IranMT891088MW308202
28Lemna gibba2SSoostan Lagoon, IranMT891089MW308203
29Lemna gibba3gGovaver, IranMT891092MW308206
30Lemna gibba4LALangarud paddy, IranMT891093MW308207
31Lemna gibba5SSoostan Lagoon, IranMT891091MW308205
32Lemna gibba6LALangarud paddy, IranMT891090MW308204
33Lemna gibba7SSoostan Lagoon, IranMT891094MZ422528
34Lemna gibba8SSoostan Lagoon, IranMT891095MZ422529
35Spirodela polyrhiza1AMAmir kelaye international lagoon, IranMT891096MW308208
36Spirodela polyrhiza2AMAmir kelaye international lagoon, IranMT891097MW308209
37Spirodela polyrhiza3LALangarud lagoon, IranMT891098MW308210
38Spirodela polyrhiza4AMAmir kelaye international lagoon, IranMT891099MW308211
39Spirodela polyrhiza5BMMansoori pond, IranMT891100MW308212
40Spirodela polyrhiza6AMAmir kelaye international lagoon, IranMT891101MZ422536

Duckweed species investigated for this study with their NCBI accession numbers.

DNA barcoding analysis

To analyze the genetic diversity, atpF-atpH and rps16 sequences were checked using Chromas Lite 2.6.2 (Technelysium Pty Ltd, South Brisbane, Australia) and aligned using the BioEdit Sequence Alignment Editor 7.1.3.0 (Hall, 1999). Reference sequences of atpF-atpH and rps16 markers for DNA barcoding were prepared from already identified clones from the duckweed stock collection of the University Greifswald (Germany) or taken from GenBank (Table 2). They were chosen to represent a wide geographical distribution of the species. SeqState 1.4.1 (Müller, 2005) was used to recode insertion and deletion (indel) positions using the implemented Simmons and Ochoterena simple coding algorithm, leading to a final alignment length of 798 sites including 14 indel coded sites for rps16 and a final alignment length of 661 sites including 12 indel coded sites for atpF-atpH. The indel coded alignments can be found as Supplementary Material 1 and 2. Subsequently, TCS 1.23 (Clement et al., 2000) was used with default settings to build haplotype networks for each chloroplast marker. Based on the haplotype results, the alignments were collapsed to unique haplotypes for which a maximum-likelihood tree was built using iqTREE 2.1.3 (Minh et al., 2020) with 1000 bootstrap replicates. The implemented ModelFinder (Kalyaanamoorthy et al., 2017) chose F81+F (atpF-atpH) and K3Pu+F (rps16) as the best-fit models according to the Bayesian information criterion. Finally, DnaSP 6.12.03 (Rozas et al., 2017) was run to count polymorphic and parsimony informative sites and to estimate nucleotide and haplotype diversity.

Table 2

RowReference sequencesStrainOriginAccession number
atpF-atpHrps 16
1Lemna minor7123Canada, Saskachewan, SaskatoonMG000397*
2Lemna minor8292Iran, Mazanda, Ramsar, Ghassem Abbath**
3Lemna minor9441Germany, Marburg (clone St)**
4Lemna turionifera6573USA, Montana, Lincoln Co.MG775403*
5Lemna turionifera7683Korea, Kyonggi, SosaMG775404*
6Lemna turionifera9434Russia, Lake BaikalMG775405*
7Lemna gibba7589USA, California, Los Angeles Co., CovinaGU454219*
8Lemna gibba7741Italy, Sicilia, Siracusa (clone G3)KX212887*
9Lemna gibba8703Japan, Honshu AichiGU454222*
10Spirodela polyrhiza7373Egypt, Mahallet, El RahabeinHG938145HG938250
11Spirodela polyrhiza7498USA, North Carolina, Durham Co., DurhamGU454204HG938251
12Spirodela polyrhiza9500Germany, Jena, Porstendorf 1967 (clone SJ)GU454208HG938257

Reference sequences of atpF-atpH and rps16 markers for DNA barcoding.

*sequenced but no Genbank number yet.

Culture conditions for growth factor analysis

Forty geographic isolates representing four species from two genera were used for biomass accumulation and DT analysis as described in Table 1. According to the ISO 20079 protocol (ISO 20079, 2005), the sterilized duckweed samples were precultivated for 1 month to acclimatize the clones to the cultivation conditions. Nutrient media were replenished every week. A single clone with the same frond number from each species was used for initial inoculation of 50 ml nutrient medium in glass jars covered with plastic caps. All glasses were kept under axenic conditions at 25°C in a standard growth chamber. The investigated duckweed species showed optimum growth in different media (unpublished results). For this reason, we used a specific nutrient medium for each species instead of the Steinberg medium specified in the ISO 20079 protocol as mentioned. The growth factors were determined starting with a four-frond colony for each species, and the initial weight was determined. The main cultivation phase lasted 7 days, taking care that the fronds never completely covered the surface of the medium, which may limit growth.

Calculation of growth parameters

All growth parameters were determined at the onset of the experiment (t0) and 7 days later (t7). The number of fronds (FN0 and FN7), fresh weight (FW0 and FW7), and dry weight (DW0 and DW7) were measured. At the initiation of the experiment, the frond numbers were recorded. Then, an equal frond number and size was surface-dried by filter paper and weighed (FW0). These reference samples were dried at 37°C for 72 h to determine the dry weight of the preliminary inoculum (DW0). Frond number, fresh weight, and dry weight of fronds at t7 were determined as for the reference samples from t0. Three independent experiments were conducted with each of the clones.

RGR was calculated using Equation (1) (Naumann et al., 2007; Ziegler et al., 2015). This equation was simplified to Equation (2) for better interpretation of growth potentials. The values of measured parameters x (Frond number or fresh and dry weight) in two time points (t0 and t7) were placed in Equation (2).

The RGR unit is based on time (per day). DT (days) or biomass accumulation (per day), was calculated by Equation (3), when RGR is measured with frond number values or fresh weight and dry weight of fronds at the two time points, respectively.

The yield obtained from the initial inoculum of one frond (or 1 mg duckweed biomass) after 7 days of cultivation is known as RY. It was calculated using Equation (4):

RY is equal to lnxt7, and x is one of the growth parameters measured in the experiment, such as FN, FW, and DW at t0 (lnxt0) and at t7 (lnxt7). The RY of one frond or 1 mg duckweed initial inoculum after 7 days has the unit per week.

Data analysis

Statistical analysis was carried out in SPSS 16.0 (IBM, USA). The normality of the data was confirmed by the Kolmogorov–Smirnov test. Therefore, parametric methods were used to compare the means. The variation of means among groups was compared with one-way ANOVA using the Student–Newman–Keuls test (SNK), a post hoc test for analysis of the differences in means, at the level of P ≤.05.

Results

DNA barcoding of duckweed ecotypes based on rps16 and atpF-atpH sequences

A total of 40 duckweed clones were collected from the north of Iran (lakes of Tonekabon, Lahijan, and Langarud) and the Kermanshah Govaver River. The collected duckweed accessions were morphologically determined and validated by molecular methods, i.e., DNA barcoding (Figure 1).

Figure 1

Figure 1

Four investigated duckweed species native to Iran; (A)Spirodela polyrhiza.(B)Lemna minor.(C)L. gibba in dorsal view. (D)L. gibba in ventral view. (E)L. turionifera.

In summary, a total of 24 ecotypes of L. minor, eight ecotypes of L. gibba, six ecotypes of S. polyrhiza, and two ecotypes of L. turionifera, a rare species for Iran, were identified with the chloroplast fragments rps16 and atpF-atpH. All identified clones were successfully propagated to produce pure clones.

The species could be very well-distinguished by both chloroplast markers (rps16 and atpF-atpH) as represented in the maximum-likelihood phylogenetic trees (Supplementary Figures 2A, B). Comparison of the two sequence alignments for both markers separately shows that rps16 has a higher haplotype and nucleotide diversity than atpF-atpH although differences between different haplotypes within one species are most often caused by indels (Table 3, Figures 2, 3). For rps16, both clones identified as L. turionifera from Iran showed the same haplotype (LT1) as the reference sequence of clone 9434 from Lake Baikal, Russia. For L. minor the Iranian clones were identical to the haplotype (LM1) of 8292, a reference clone from Iran, too. Two further clones (25C – haplotype LM3 and 12W – haplotype LM4) showed one or two additional bases but are more similar to the main haplotype LM1 found for Iran than to the other two reference clones from Canada and Germany (haplotype LM2) (Figures 2A, B). For the marker atpF-atpH only L. gibba revealed different haplotypes, in which the Iranian clones differed by an additional stretch of three A’s (haplotype lg1) from the three reference clones, which showed the same haplotype (lg2) (Figures 3A, B).

Table 3

speciesnumberclonesrps16atpF-atpH
alignmentlength (bp)*indel codedsitesPSPIHnumHd ± SDл ± SDalignmentlength (bp)*indel codedsitesPSPIHnumHd ± SDл ± SD
S. polyrhiza9798140010.000±
0.000
0.0000±
0.0000
661120010.000±
0.000
0.0000±
0.0000
L. gibba110010.000±
0.000
0.0000±
0.0000
1120.436±
0.133
0.0008±
0.0002
L. turionifera52220.600±
0.175
0.0016±
0.0005
0010.000±
0.000
0.0000±
0.0000
L. minor273240.276±
0.109
0.0005±
0.0002
0010.000±
0.000
0.0000±
0.0000

Alignment characteristics for the two investigated chloroplast markers.

PS, polymorphic sites; PI, parsimony informative sites; Hnum, number of haplotypes; Hd, haplotype diversity; л, nucleotide diversity; SD, standard deviation.* including indel coded sites.

Figure 2

Figure 2

Molecular genetic results of the cp marker rps16. (A) Maximum-likelihood tree of unique haplotypes found for all investigated clones. Bootstrap values based on 1,000 replicates. Spirodela polyrhiza was set as outgroup. Scale indicates number of substitutions per site. (B) Identification of haplotypes and their differences for all investigated clones. Substitutions are given on the lines. - denotes deletion.

Figure 3

Figure 3

Molecular genetic results of the cp marker atpF-atpH. (A) Maximum-likelihood tree of unique haplotypes found for all investigated clones. Bootstrap values based on 1,000 replicates. Spirodela polyrhiza was set as outgroup. Scale indicates number of substitutions per site. (B) Identification of haplotypes and their differences for all investigated clones. Substitutions are given on the lines. - denotes deletion.

Biomass accumulation and doubling time

The growth potential measured as RGR and RY based on fresh and dry weight and DT based on frond numbers of 40 geographic duckweed isolates under axenic cultivation conditions are shown in Table 4. Overcrowding of populations was not observed during the experiment or at the end after 7 days. This ensured that the growth was not inhibited by intraspecific competition. In addition, axenic cultivation prevented inhibitory effects of undesirable microorganisms (Ziegler et al., 2015). During the 7 days of the experiment, the increase in frond number, fresh weight, and dry weight never deviated from an exponential progression.

Table 4

RowSpeciesstrainDTFNRGRFWRGRDWBAFWBADWRYFWRYDW
1Lemna minor1C3.77 ± 0.110.337 ± 0.0290.478 ± 0.012.10 ± 0.181.51 ± 0.0228.94 ± 5.553.71 ± 0.02
2Lemna minor2C3.44 ± 0.120.279 ± 0.0050.255 ± 0.012.49 ± 0.052.72 ± 0.0630.68 ± 1.142.14 ± 0.09
3Lemna minor3a5.97 ± 0.420.215 ± 0.0010.307 ± 0.013.10 ± 0.062.39 ± 0.0214.44 ± 0.252.05 ± 0.03
4Lemna minor4BM3.73 ± 0.0010.338 ± 0.0060.277 ± 0.012.05 ± 0.032.50 ± 0.0629.69 ± 1.141.75 ± 0.09
5Lemna minor5C3.67 ± 0.150.343 ± 0.0020.256 ± 0.012.02 ± 0.012.72 ± 0.1134.87 ± 0.432.05 ± 0.14
6Lemna minor6a7.05 ± 0.02 (-)0.301 ± 0.0250.287 ± 0.022.35 ± 0.202.44 ± 0.1428.21 ± 4.841.75 ± 0.2
7Lemna minor7W4.04 ± 0.010.215 ± 0.0010.332 ± 0.023.23 ± 0.022.18 ± 0.0328.06 ± 0.202.94 ± 0.09
8Lemna minor8C5.34 ± 0.160.202 ± 0.0050.156 ± 0.01 (-)3.44 ± 0.084.52 ± 0.37 (-)12.65 ± 0.431.30 ± 0.12
9Lemna minor9a3.74 ± 0.010.277 ± 0.0070.294 ± 0.0032.50 ± 0.072.36 ± 0.0622.78 ± 1.141.10 ± 0.06
10Lemna minor10C3.40 ± 0.040.260 ± 0.0020.273 ± 0.012.67 ± 0.012.68 ± 0.0115.86 ± 0.011.40 ± 0.01
11Lemna minor11C *2.96 ± 0.02 (+)0.373 ± 0.017 (+)0.563 ± 0.011.87 ± 0.08 (+)1.23 ± 0.0337.09 ± 4.27 (+)2.59 ± 0.23
12Lemna minor12W4.97 ± 0.010.293 ± 0.0140.302 ± 0.022.38 ± 0.122.31 ± 0.1220.31± 21.35 ± 0.14
13Lemna minor13C4.43 ± 0.010.306 ± 0.0070.191 ± 0.012.27 ± 0.053.63 ± 0.1423.77 ± 1.142.25 ± 0.11
14Lemna minor14BM4.04 ± 0.010.301 ± 0.0030.258 ± 0.012.23 ± 0.032.71 ± 0.1229.44 ± 0.713.24 ± 0.26
15Lemna minor15C4.97 ± 0.020.350 ± 0.0050.369 ± 0.021.98 ± 0.031.90 ± 0.1225.25 ± 0.852.74 ± 0.43
16Lemna minor16C4.43 ± 0.010.326 ± 0.0020.489 ± 0.012.12 ± 0.011.42 ± 0.0426.24 ± 0.281.55 ± 0.14
18Lemna minor18C4.68 ± 0.110.287 ± 0.0290.326 ± 0.032.40 ± 0.302.18 ± 0.1919.40 ± 0.241.35 ± 0.26
19Lemna minor19C4.61 ± 0.450.234 ± 0.0150.216 ± 0.023 ± 0.193.26 ± 0.2421.79 ± 2.281.85 ± 0.20
21Lemna minor21BM5.07 ± 0.270.256 ± 0.0060.327 ± 0.022.72 ± 0.072.14 ± 0.1026.73 ± 1.141.10 ± 0.12 (-)
22Lemna minor22a3.78 ± 0.260.317 ± 0.0190.512 ± 0.012.21 ± 0.131.36 ± 0.0435.60 ± 4.553.64 ± 0.37
23Lemna minor23C3.96 ± 0.190.260 ± 0.0090.612 ± 0.012.68 ± 0.091.15 ± 0.0133.14 ± 1.998.21 ± 0.85 (+)
24Lemna minor24W5.07 ± 0.270.229 ± 0.0070.654 ± 0.03 (+)3.03 ± 0.091.06 ± 0.01 (+)17.35 ± 0.861.95 ± 0.03
25Lemna minor25C2.97 ± 0.020.175 ± 0.0020.406 ± 0.073.95 ± 0.051.90 ± 0.3520.31 ± 0.294.59 ± 0.34
26Lemna minor26W5.07 ± 0.270.094 ± 0.010 (-)0.229 ± 0.017.62 ± 0.84 (-)3.04 ± 0.137.95 ± 0.57 (-)2.14 ± 0.14
Lemna minormean4.38 ± 0.050.274 ± 0.0020.349 ± 0.022.77 ± 0.042.31 ± 0.0324.61 ± 0.423.70 ± 0.05
17Lemna turionifera17AM9.57 ± 0.030.196 ± 0.0330.199 ± 0.033.36 ± 0.363.78 ± 0.617.20 ± 1.570.95 ± 0.20
20Lemna toriunifera20AM4.22 ± 0.300.433 ± 0.0030.628 ± 0.021.60 ± 0.011.11 ± 0.0414.38 ± 0.290.85 ± 0.14
Lemna turioniferamean6.90 ± 0.130.314 ± 0.0180.413 ± 0.022.48 ± 0.192.45 ± 0.2810.79 ± 0.930.90 ± 0.03
27Lemna gibba1LA2.50 ± 0.01 (-)0.274 ± 0.002 (-)0.329 ± 0.01 (-)2.53 ± 0.02 (-)2.11 ± 0.0260.49 ± 0.99 (-)3.98 ± 0.01 (-)
28Lemna gibba2S *2.16 ± 0.040.393 ± 0.011 (+)0.376 ± 0.011.77 ± 0.05 (+)1.85 ± 0.06108.9 ± 7.94 (+)6.96 ± 0.57 (+)
29Lemna gibba3g2.16 ± 0.070.379 ± 0.0110.391 ± 0.031.83 ± 0.051.76 ± 0.0669.10 ± 5.114.68 ± 0.40
30Lemna gibba4LA2.17 ± 0.030.355 ± 0.0020.338 ± 0.011.95 ± 0.012.05 ± 0.0297.13 ± 1.426.36 ± 0.11
31Lemna gibba5S2.35 ± 0.020.306 ± 0.0040.403 ± 0.012.26 ± 0.031.72 ± 0.02 (+)75.98 ± 2.275.02 ± 0.20
32Lemna gibba6LA2.14 ± 0.01 (+)0.321 ± 0.0020.333 ± 0.032.16 ± 0.012.15 ± 0.01 (-)93.20 ± 1.426.66 ± 0.06
33Lemna gibba7S2.48 ± 0.050.298 ± 0.0070.420 ± 0.02 (+)2.33 ± 0.061.80 ± 0.0160.73 ± 3.135.92 ± 0.03
34Lemna gibba8S2.16 ± 0.010.333 ± 0.0030.412 ± 0.012.08 ± 0.011.77 ± 0.0293.10 ± 0.016.17 ± 0.06
Lemna gibbamean2.26 ± 0.010.332 ± 0.0030.375 ± 0.012.12 ± 0.011.90 ± 0.0282.34 ± 1.655.72 ± 0.13
35Spirodela polyrhiza1AM7.22 ± 0.93 (-)0.323 ± 0.0110.425 ± 0.042.16 ± 0.081.63 ± 0.0228.70 ± 2.282.74 ± 0.09 (-)
36Spirodela polyrhiza2AM *2.84 ± 0.020.472 ± 0.003 (+)0.498 ± 0.051.47 ± 0.01 (+)1.43 ± 0.1459.25 ± 1.14 (+)10.85 ± 0.03
37Spirodela polyrhiza3LA3.86 ± 0.230.379 ± 0.0110.382 ± 0.011.83 ± 0.051.82 ± 0.0652.36 ± 3.984.08 ± 0.34
38Spirodela polyrhiza4AM *2.33 ± 0.05 (+)0.200 ± 0.0090.784 ± 0.02 (+)3.48 ± 0.160.88 ± 0.02 (+)40.04 ± 2.5650.88 ± 4.83 (+)
39Spirodela polyrhiza5BM2.93 ± 0.150.214 ± 0.0180.650 ± 0.023.30 ± 0.271.07 ± 0.0416.36 ± 2 (-)40.07 ± 0.04
40Spirodela polyrhiza6AM3.9 ± 0.210.117 ± 0.014 (-)0.149 ± 0.07 (-)6.16 ± 0.73 (-)13.55 ± 6.35 (-)17.84 ± 1.7134.01 ± 0.01
Spirodela polyrhizamean3.85 ± 0.210.284 ± 0.0020.481 ± 0.013.07 ± 0.083.40 ± 1.0535.76 ± 0.4723.78 ± 0.76

Fresh and dry weight biomass accumulation and doubling time of 40 geographical isolates of native Iranian duckweed representing four species.

DT, doubling time (days); BA, biomass accumulation (per day); RGR, relative growth rate (per day); RY, relative yield (per week). Mean values for each species are presented in high light rows. Values are mean ± SE.

(+) indicates the highest growth rate values and (-) indicates the lowest growth rate values.

*Strains with the highest growth rate.

The mean RGR for all 40 investigated ecotypes was 0.301 per day for fresh weight, ranging from 0.094 to 0.472 per day for individual clones belonging to L. minor 26W and S. polyrhiza 2AM, respectively. Additionally, the mean RGR for dry weight was 0.435 per day with a range from 0.149 to 0.784 per day for Spirodela clones (Table 4).

As shown in Figure 4A, the mean RGR based on fresh weight for four species belonging to two duckweed genera ranged from 0.284 per day for S. polyrhiza to 0.332 per day for L. gibba. The mean RGR for L. minor is 0.274 per day. The RGR for L. gibba was significantly higher than for any other species (P<.05). This was due to a higher weight gain and more fronds after 7 days of cultivation. There were no significant differences between the mean RGRFW of L. minor and S. polyrhiza due to a wide range of RGR values among L. minor clones.

Figure 4

Figure 4

The mean fresh and dry weight–based relative growth rate (biomass accumulation, relative weekly yield, and doubling time based on frond number) of the four investigated duckweed species. (A) The mean RGR for the species represented by several clones. The number of clones per species is given above. (B) The mean value of DTFN between species is significantly different. Lemna gibba by doubling its fronds every 2.26 days is faster than other species. (C) The mean BAFW of four species of two genera is shown. (D) The maximum rate of mean RYFW of four species belongs to L. gibba (82.3 per week). (E) The mean of BADW is not significantly different among species. (F) The mean dry weight (mg) produced after 1 week (RY) from primary inoculum among species is significantly different (P<.05). In all figures: The column height shows the mean growth parameters of the four investigated species. Error bar was indicated in the figures. Letters above the columns indicate significance according to ANOVA: means of columns marked with the same letter (either ‘a’ or ‘b’) do not differ to a statistically significant extent; differences are statistically significant when the columns are marked with single, different letters. The means of the columns marked with ‘ab’ do not differ significantly from means of columns marked with either ‘a’ or ‘b’.

The dry weight analysis gave similar results to fresh weight. Among the 24 L. minor clones, the highest RGRDW values were obtained by clones 11C, 23C, and 24W with 0.563, 0.612, and 0.654 per day, respectively. These values were not significantly different from each other (P<.05). The clone with the significantly highest RGRDW among the eight L. gibba clones was clone 7S with 0.420 per day. Among the six clones of S. polyrhiza, 4AM showed the highest value with RGRDW = 0.784 per day. Fresh and dry weight RGR within species are significantly different (P<.05), especially in L. minor and S. polyrhiza. Some of the ecotypes of L. minor and S. polyrhiza studied as well as one of the clones of L. turionifera had an RGR higher than 0.600 per day (Table 4). This is a significantly higher RGR value compared with other ecotypes.

DT, which is based on the number of fronds, and BA, based on the fresh or dry weight, reflect the RGR value but numerically in the opposite way (see Equation (3) above). It indicates how much time is needed to double the number of fronds or biomass. DT based on frond number was investigated for all ecotypes. The lowest DT (rapid growth) within species was measured for L. minor 11C (2.96 days), L. gibba 6LA (2.14 days), and S. polyrhiza 4AM (2.33 days), which doubled their frond number every 51 to 71 h (Table 4). A comparison of the mean values for DT of the investigated clones of the four species is shown in Figure 4B. The mean value of frond doubling time for L. gibba with 2.26 days is significantly lower than that of the other species (P<.05).

The mean biomass accumulation based on fresh weight (BAFW) of eight L. gibba clones was 2.12 days. Thus, this species had significantly higher productivity than the other species (P<.05; Figure 4C). On the other hand, the mean biomass accumulation based on dry weight (BADW) of L. gibba is consistent with BAFW. However, this is not significantly different between species (Figure 4E). The significantly fastest BAFW within L. minor was observed for strain 11C with 1.87 days (P<.05; Table 4). The two clones of the rare L. turionifera had BA rates of 1.60 and 3.78 days. To increase the accuracy of growth parameters in L. turionifera, it is necessary to continue the work with more ecotypes.

The RY of biomass accumulation based on fresh weight (RYFW) of all investigated ecotypes ranged from 108.92 (L. gibba 2S) to 7.2 per week (L. turionifera 17AM). The highest mean RYFW among the species belonged to L. gibba (82.34 per week) and the lowest mean RYFW to Spirodela (35.76 per week) and L. minor (24.61 per week). The differences among the four species were statistically significant (P<.05; Figure 4D). Analysis of the RY value showed that the results were consistent with BA and DT. In addition, intraspecies data for L. minor showed that the highest relative yield of 37.09 per week was obtained by L. minor 11C due to its better growth potential. Despite the results for mean fresh weight (Figures 4D, F), where L. gibba showed the highest RYFW, S. polyrhiza had the highest value in RYDW in the species comparison with 23.78 per week (Figure 4F) while L. gibba (5.72 per week) and L. minor (3.70 per week) were both in the lower range.

Discussion

In the present study, biomass production screening was performed based on growth potential analysis on native duckweed species from Iran coupled with DNA barcoding based on two standard chloroplast markers. Growth parameters were used to identify the most productive ecotypes of each of the four native duckweed species (Spirodela polyrhiza, Lemna minor, L. gibba, and L. turionifera) among the 40 clones studied. Lemna minor 11C, L. gibba 2S, and S. polyrhiza 4AM showed growth rates and a relative yield even higher than any terrestrial plants reported in previous studies (Ziegler et al., 2015; Koca and Erekul, 2016). These data demonstrate the importance of comprehensive studies on duckweed for biomass accumulation in biological production systems and food security. In addition, for the first time, the species L. turionifera was detected for Iran based on DNA barcoding analysis, and the two clones were included in the biomass screening.

Primary identification of duckweed based on morphological characters identified three species belonging to two genera: Spirodela polyrhiza, Lemna minor, and L. gibba. However, subsequent DNA barcoding analysis revealed a fourth species: L. turionifera, which is not easily distinguishable from L. minor due to the strong reduction in their morphology. After an appropriate literature research, L. turionifera was not listed for Iran in Landolt’s monograph (Landolt, 1986), nor in the online Flora of Iran (2022) (http://flora-iran.com/central-herbarium-of-tehran-university/plant-list/), nor in the two online platforms “Plants of the World Online (2022)” (https://powo.science.kew.org/) and “Global Biodiversity Information Facility – GBIF (2022)” (https://www.gbif.org/).

The variation of the studied sequences is rather low, as is known for several duckweed species (Borisjuk et al., 2015; Chen et al., 2022), but also to be expected for a DNA barcode as they should show lower genetic variation within than between species (Dasmahapatra and Mallet, 2006). Interestingly, with the exception of two accessions, all other accessions identified as L. minor show the same haplotype as the reference sequence of rps16 for clone 8292 from Iran, which was already mentioned in Landolt and Urbanska-Worytkiewicz (1980) and, thus, has been kept in culture for more than 40 years, which again could indicate a relatively constant haplotype pool. However, further studies are needed to reach similar conclusions as for S. polyrhiza, namely, that the mutation rate in this species is very low (Ho et al., 2019; Xu et al., 2019). Because the molecular markers used do not allow us to draw conclusions about hybridization events, the possibility that the L. minor clones identified here may be hybrids between L. minor and L. turionifera cannot be ruled out. Both species occur in the area, and hybridization between these two species has already been demonstrated by Braglia et al. (2021).

RGR is an important factor to show physiological responses of plants to light, temperature, CO2, and nutrients, but its interpretation is less intuitive (Buxbaum et al., 2022). Therefore, it is often mathematically transformed and reported as BA or RY, which are common parameters for large-scale screening of plant growth potential (Ziegler et al., 2015). Fast-growing species show higher RGR under standard cultivation conditions. Higher RGR is due to efficient nutrient and CO2 uptake in fronds (Naumann et al., 2007; Ziegler et al., 2015; Ghanem et al., 2019). High RGR causes duckweed to rapidly double its frond number and biomass, which increases its photosynthetically active surface area per unit area. One of the most important variations is the nutrient medium. Many studies demonstrate that duckweed exhibits optimal growth potential and accumulates more biomass in species-specific nutrient media (Kittiwongwattana and Vuttipongchaikij, 2013; Muranaka et al., 2015; Ghanem et al., 2019). This study was the first to use a nutrient medium optimized for growth for each species (data in preparation).

In this study, a relative linear relationship between RGR, DT or BA, and RY was investigated (Figure 4). In addition, growth factors based on dry weight, especially RYDW, are more reliable parameters to study the stability of biomass gain in the duckweed family. Approximately 92%–94% of the fresh weight of duckweed consist of water, which is lost after drying (Pagliuso et al., 2022). This is confirmed by the relatively low yield (RYDW) results for L. gibba shown in Figure 4F. Despite the highest mean RYFW value (Figure 4D, F) for L. gibba, almost 90% of the fresh weight was lost to drying. In contrast, for Spirodela, only 35% of the fresh weight was lost to drying, and on average, 65% of the weight was retained as dried biomass (Table 4 and Figure 4F).

Based on the analysis of growth potential of 40 ecotypes, ecotypes with better growth potential under standard growing conditions were identified. Spirodela polyrhiza and L. gibba showed higher average growth factors among the studied species in this experiment. On the other hand, L. minor, the most widespread and easily manipulated species, with ecotype 11C provided one of the most productive duckweed clones in this experiment for future research. Considering that growth potential is shown to have a wide range of values within genera and species and, thus, significant differences among clones and ecotypes, it is concluded that species or genus specificity is not a reliable method for screening growth potential. In other words, growth parameters determined for one species or genus cannot be generalized to all clones of a species. This is consistent with Bergmann et al. (2000), who suggest focusing on the geographic isolate (ecotypes) rather than the species level. Ziegler et al. (2015) also confirm that screening on the ecotype level is the most reliable method for screening growth factors in duckweed. We found three ecotypes with high growth potential such as S. polyrhiza 2AM with an RGRFW of 0.472 per day, which is higher than the value reported by Ziegler et al. (2015), ranging from 0.168 to 0.386 per day for seven S. polyrhiza clones. Other dependent parameters, such as BAFW (1.47 per day) and RYFW (59.25 per week), are consistent with those reported in the literature (Ziegler et al., 2015). In addition, L. gibba 2S was found to be the fastest ecotype in biomass accumulation with 1.77 per day. Similarly, the relative yield after 7 days of inoculation is the highest value (108.92 per week) for L. gibba 2S, which is consistent with those from previous reports (Ziegler et al., 2015). It is noteworthy that L. gibba already has a good potential to gain biomass in a short time. This is due to the rapid proliferation rate and gibbous fronds in this species. As mentioned earlier, under unsuitable culture conditions, L. gibba loses its gibbosity and looks like L. minor. NF is the best culture medium for L. gibba to form gibbous fronds and increase the rate of biomass accumulation (Muranaka et al., 2015). In addition, the widely distributed duckweed species L. minor with ecotype 11C has an RGRFW value of 0.373 per day. This is comparable to the value of RGRFW = 0.422 per day reported by Chakrabarti et al. (2018), where L. minor was cultivated with high-efficiency organic manure. The RGRDW of L. minor 11C was 0.563 per day, which was among the highest values in the study, while Petersen et al. (2021) obtained a maximum RGRDW of 0.23 per day for L. minor under nonsterile conditions and with high-yielding agricultural fertilizer treatment. Despite the different cultivation conditions, the highly concentrated nitrate-N medium and light intensity (270 μmol/m2/s higher than in the present study) were effective factors for optimal growth.

Duckweed is described as one of the fastest angiosperms due to its ability to double its biomass in a short time period. Demmig-Adams et al. (2022) suggest that the rapid growth of duckweed, even under limited light conditions, may be related to relatively thin photosynthetic organs (fronds) without complex structures on the water surface, allowing all chloroplasts to be involved in sugar production. The free-floating fronds with high availability of nutrient resources have higher photosynthetic yields. Terrestrial plants, on the other hand, use a significant amount of sugar to build the complex structures of their stems, leaves, and roots, which is a time-consuming process that involves the production of some organs that are unusable for human nutrition. The presented results shown in Table 4 confirm the superiority of growth rate of duckweeds by comparing their RGR with the RGR of crops. Of course, plant species differ greatly in their relative growth rate even when compared under similar environmental conditions (Tomlinson et al., 2014). However, duckweed is shown to have a higher RGRDW compared with many crops despite its reduced and leaf-like structure that makes it one of the lightweights among crops. Poorter (1989) reports RGRDW of eight herbaceous wild species under optimized growth conditions with the highest RGR value of 0.268 per day for Urtica dioica. Potter and Jones (1977) report RGRDW after 28 days for nine species of important crops with a value ranging from 0.202 per day for soybean to 0.391 per day for sorghum. In contrast, an RGR value of 0.255 per day was determined for maize under the best experimental conditions. The data presented in Table 4 and Figure 4 show that most of the duckweed clones that were studied had an RGR based on fresh and dry weight (and corresponding RY) that was higher than the 0.255 per day described for maize. Overall, the superior growth characteristics and short harvest time of duckweed compared with crops may lead to higher biomass and economic production.

As reported by the Food and Agriculture Organization of the United Nations (FAO), food security is one of the most important challenges in the world (https://www.fao.org/state-of-food-security-nutrition/2021). Due to the increase in world population, reduction in food resources, depletion of nutrients in soils, and global climate change leading to a reduction in crop production, it is imperative to pay more attention to alternatives with low water and soil requirements that are cost-efficient, have short harvesting times, and produce more biomass. Duckweed as an aquatic crop has several advantages over terrestrial crops: it absorbs nutrients directly from water, is easy to grow and harvest, has low water requirements (due to its short growing season), and does not compete with crops in agricultural land use, allowing for higher biomass production per hectare even in dry areas (Toyama et al., 2017; Tursi, 2019). To date, interest in large-scale cultivation of duckweed in greenhouses and under protective structures has grown through commercial companies, such as LENTEIN™ (https://www.parabel.com/) and Rubisco Foods (2022) (https://rubiscofoods.com/).

Conclusion

The results of this study show that three selected ecotypes of Lemna and Spirodela species can provide high yields of fresh and dry biomass under optimal growth conditions. The data confirm that most ecotypes of duckweed can grow faster than traditional crop plants. However, this high RGR obtained under the optimized culture conditions for the selected ecotypes may be different under real environmental conditions. However, the fact that duckweed as an emergent crop has a high relative yield and accumulates more biomass in a short time was also confirmed under agricultural cultivation conditions. The data in this study illustrate the numerous potentials of the selected duckweed ecotypes for commercial biomass production and biotechnological application. However, several strategies are needed to optimize duckweed growth with low cost, simplicity, and scalability.

Funding

This research was supported partially by grants from Iran National Science Foundation (INSF) (97024978), Center for International Scientific Studies and Collaboration (CISSC) (990160) and National Institute of Genetic Engineering and Biotechnology (NIGEB) (980301-II-714).

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.

Statements

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

AS is the executor of plan, ET is PhD student and project manager, MB contributed to DNA barcoding analysis, FF, SS, NR and MA cooperated in biomass measurements, MJ organized the database. All authors contributed to manuscript revision, read, and approved the submitted version.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.1034238/full#supplementary-material

References

  • 1

    AppenrothK. J.BorisjukN.LamE. (2013). Telling duckweed apart: genotyping technologies for the lemnaceae. Chin. J. Appl. Environ. Biol.19, 110. doi: 10.3724/SP.J.1145.2013.00001

  • 2

    AppenrothK. J.SreeK. S.BogM.EckerJ.SeeligerC.BöhmV.et al. (2018). Nutritional value of the duckweed species of the genus wolffia (Lemnaceae) as human food. Front. Chem.6, 483. doi: 10.3389/fchem.2018.00483

  • 3

    AppenrothK.-J.SreeK. S.BöhmV.HammannS.VetterW.LeitererM.et al. (2017). Nutritional value of duckweeds (Lemnaceae) as human food. Food Chem.217, 266273. doi: 10.1016/j.foodchem.2016.08.116

  • 4

    BergmannB. A.ChengJ.ClassenJ.StompA. M. (2000). In vitro selection of duckweed geographical isolates for potential use in swine lagoon effluent renovation. Bioresour. Technol.73 (1), 1320. doi: 10.1016/S0960-8524(99)00137-6

  • 5

    BogM.AppenrothK. J.SreeK. S. (2019). Duckweed (Lemnaceae): its molecular taxonomy. Front. Sustain. Food Syst.3, 117. doi: 10.3389/fsufs.2019.00117

  • 6

    BogM.AppenrothK. J.SreeK. S. (2020). Key to the determination of taxa of Lemnaceae: an update. Nord. J. Bot.38 (8), e02658. doi: 10.1111/njb.02658

  • 7

    BorisjukN.ChuP.GutierrezR.ZhangH.AcostaK.FriesenN.et al. (2015). Assessment, validation and deployment strategy of a two-barcode protocol for facile genotyping of duckweed species. Plant Biol.17, 4249. doi: 10.1111/plb.12229

  • 8

    BragliaL.LauriaM.AppenrothK. J.BogM.BreviarioD.GrassoA.et al. (2021). Duckweed species genotyping and interspecific hybrid discovery by tubulin-based polymorphism fingerprinting. Front. Plant Sci.12, 625670. doi: 10.3389/fpls.2021.625670

  • 9

    BuxbaumN.LiethJ. H.EarlesM. (2022). Non-destructive plant biomass monitoring with high spatio-temporal resolution via proximal RGB-d imagery and end-to-End deep learning. Front. Plant Sci.13, 758818. doi: 10.3389/fpls.2022.758818

  • 10

    CaoH. X.FourounjianP.WangW. (2018). “The importance and potential of duckweeds as a model and crop plant for biomass-based applications and beyond,” in Handbook of environmental materials management (Cham: Springer International Publishing), 116.

  • 11

    ChakrabartiR.ClarkW. D.SharmaJ. G.GoswamiR. K.ShrivastavA. K.TocherD. R. (2018). Mass production of Lemna minor and its amino acid and fatty acid profiles. Front. Chem.6, 479. doi: 10.3389/fchem.2018.00479

  • 12

    ChenG.StepanenkoA.LakhnekoO.ZhouY.KishchenkoO.PetersonA.et al. (2022). Biodiversity of duckweed (Lemnaceae) in water reservoirs of Ukraine and China assessed by chloroplast DNA barcoding. Plants11 (11), 1468. doi: 10.3390/plants11111468

  • 13

    ClementM.PosadaD. C. K. A.CrandallK. A. (2000). TCS: a computer program to estimate gene genealogies. Mol. Ecol.9 (10), 16571659. doi: 10.1046/j.1365-294x.2000.01020.x

  • 14

    CoughlanN. E.WalshÉ.BolgerP.BurnellG.O'LearyN.O'MahoneyM.et al. (2022). Duckweed bioreactors: Challenges and opportunities for large-scale indoor cultivation of lemnaceae. J. Clean. Prod.336, 130285. doi: 10.1016/j.jclepro.2021.130285

  • 15

    DasmahapatraK. K.MalletJ. (2006). DNA Barcodes: recent successes and future prospects. Heredity97 (4), 254255. doi: 10.1038/sj.hdy.6800858

  • 16

    Demmig-AdamsB.López-PozoM.PolutchkoS. K.FourounjianP.StewartJ. J.ZenirM. C.et al. (2022). Growth and nutritional quality of lemnaceae viewed comparatively in an ecological and evolutionary context. Plants11 (2), 145. doi: 10.3390/plants11020145

  • 17

    DoganH.ErcisliS.TemimE.HadziabulicA.TosunM.YilmazS. O.et al. (2014). Diversity of chemical content and biological activity in flower buds of a wide number of wild grown caper (Capparis ovate desf.) genotypes from Turkey. C R Acad. Bulg Sci.67 (11), 15931600.

  • 18

    EdelmanM.AppenrothK. J.SreeK. S. (2020). Duckweed: Biological chemistry and applications. Front. Sustain. Food Syst.4, 615135. doi: 10.3389/fsufs.2020.615135

  • 19

    EdelmanM.ColtM. (2016). Nutrient value of leaf vs. seed. front. Chem.4, 32. doi: 10.3389/fchem.2016.00032

  • 20

    FAO. (2013). (Food and Agriculture of the United Nations). Available at: https://www.faostat3.fao.org/faostat-gateway (accessed 30 May 2013).

  • 21

    Flora of Iran (2022) Plant-list. Available at: http://flora-iran.com/central-herbarium-of-tehran-university/plant-list/ (Accessed 30/05/2022).

  • 22

    Food and Agriculture Organization of the United Nations (2021) Food-security. Available at: https://www.fao.org/state-of-food-security-nutrition (Accessed 23/07/2022).

  • 23

    FredericM.SamirL.LouiseM.AbdelkrimA. (2006). Comprehensive modeling of mat density effect on duckweed (Lemna minor) growth under controlled eutrophication. Water Res.40, 29012910. doi: 10.1016/j.watres.2006.05.026

  • 24

    FujitaM.MoriK.KoderaT. (1999). Nutrient removal and starch production through cultivation of Wolffia arrhiza. J. Biosci. Bioeng.87, 194198. doi: 10.1016/S1389-1723(99)89012-4

  • 25

    GhanemH.HaddadA.BaydounS.Abou HamdanH.KorfaliS.ChalakL. (2019). In vitro proliferation of Lebanese Lemna minor and Lemna gibba on different nutrient media. J. Taibah Univ. Sci.13 (1), 497503. doi: 10.1080/16583655.2019.1597450

  • 26

    Global Biodiversity Information Facility – GBIF (2022). Available at: https://www.gbif.org/ (Accessed 30/05/2022).

  • 27

    HallT. A. (1999). “BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT,” in Nucleic acids symp. ser, vol. Vol. 41. (London: Information Retrieval Ltd), 9598.

  • 28

    HasanbegovicJ.HadziabulicS.KurtovicM.GasiF.LazovicB.DorbicB.et al. (2021). Genetic characterization of almond (Prunus amygdalus l) using microsatellite markers in the area of Adriatic Sea. Turk. J. Agric. For.45, 797806. doi: 10.3906/tar-2103-82

  • 29

    HeenatigalaP. P. M.SunZ.YangJ.ZhaoX.HouH. (2020). Expression of LamB vaccine antigen in Wolffia globosa (duck weed) against fish vibriosis. Front. Immunol.11, 1857. doi: 10.3389/fimmu.2020.01857

  • 30

    HoE. K.BartkowskaM.WrightS. I.AgrawalA. F. (2019). Population genomics of the facultatively asexual duckweed Spirodela polyrhiza. New Phytol.224 (3), 13611371. doi: 10.1111/nph.16056

  • 31

    IatrouE. I.StasinakisA. S.AloupiM. (2015). Cultivating duckweed Lemna minor in urine and treated domestic wastewater for simultaneous biomass production and removal of nutrients and antimicrobials. Ecol. Eng.84, 632639. doi: 10.1016/j.ecoleng.2015.09.071

  • 32

    ISO 20079 (2005). Water quality – determination of the toxic effect of water constituents and waste water on duckweed (Lemna minor) – duckweed growth inhibition test. ISO TC 147/SC 5/WG 5.

  • 33

    KalyaanamoorthyS.MinhB. Q.WongT. K.Von HaeselerA.JermiinL. S. (2017). ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods14 (6), 587589. doi: 10.1038/nmeth.4285

  • 34

    KhvatkovP.ChernobrovkinaM.OkunevaA.DolgovS. (2019). Creation of culture media for efficient duckweeds micropropagation (Wolffia arrhiza and Lemna minor) using artificial mathematical optimization models. PCTOC136 (1), 85100. doi: 10.1007/s11240-018-1494-6

  • 35

    KittiwongwattanaC.VuttipongchaikijS. (2013). Effects of nutrient media on vegetative growth of Lemna minor and Landoltia punctata during in vitro and ex vitro cultivation. Maejo Int. J. Sci. Technol.7 (1), 6069. doi: 10.14456/mijst.2013.5

  • 36

    KocaY. O.ErekulO. (2016). Changes of dry matter, biomass and relative growth rate with different phenological stages of corn. Agric. Agric. Sci. Procedia.10, 6775. doi: 10.1016/j.aaspro.2016.09.015

  • 37

    LandoltE. (1986). The Family of Lemnaceae - A Monographic Study. Vol. 1. Biosystematic Investigations in the Family of Duckweeds (Lemnaceae). Veröffentlichungen des Geobotanischen Institutes der Eidg. Technische Hochschule, Stiftung Rübel, Zurich.

  • 38

    LandoltE. (1957). “Physiologische und Ökologische untersuchungen an lemnaceen,” in Berichte der schweizerischen botanischen gesellschaft, (Bern, Switzerland: Druck und Kommissionsverlag von Büchler & Co.) 67, 271410.

  • 39

    LandoltE.Urbanska-WorytkiewiczK. (1980). “List of the studied lemnaceae samples: origin and chromosome numbers,” in Biosystematische untersuchungen in der familie der wasserlinsen (Lemnaceae). (Zürich, Switzerland: Veröff. Geobot. Inst. Rübel) 1, 205247.

  • 40

    LasfarS.MonetteF.MilletteL.AzzouzA. (2007). Intrinsic growth rate: a new approach to evaluate the effects of temperature, photoperiod and phosphorus–nitrogen concentrations on duckweed growth under controlled eutrophication. Water Res.41 (11), 23332340. doi: 10.1016/j.watres.2007.01.059

  • 41

    Lemnature AquaFarms - A Sustainable Aquaculture Company (2022). Available at: https://www.parabel.com/ (Accessed 06/05/2022).

  • 42

    LiuY.ChenX.WangX.FangY.ZhangY.HuangM.et al. (2019). The influence of different plant hormones on biomass and starch accumulation of duckweed: A renewable feedstock for bioethanol production. Renew. Energy.138, 659665. doi: 10.1016/j.renene.2019.01.128

  • 43

    MesJ. J.EsserD.OosterinkE.van den DoolR. T.EngelJ.de JongG. A.et al. (2022a). A controlled human intervention trial to study protein quality by amino acid uptake kinetics with the novel Lemna protein concentrate as case study. Int. J. Food. Sci. Nutr.73 (2), 251262. doi: 10.1080/09637486.2021.1960958

  • 44

    MesJ. J.EsserD.SomhorstD.OosterinkE.van der HaarS.UmmelsM.et al. (2022b). Daily intake of Lemna minor or spinach as vegetable does not show significant difference on health parameters and taste preference. Plant Foods Hum. Nutr.77 (1), 121127. doi: 10.1007/s11130-022-00952-9

  • 45

    MinhB. Q.SchmidtH. A.ChernomorO.SchrempfD.WoodhamsM. D.Von HaeselerA.et al. (2020). IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol.37 (5), 15301534. doi: 10.1093/molbev/msaa015

  • 46

    MüllerK. (2005). SeqState - primer design and sequence statistics for phylogenetic DNA data sets. Appl. Bioinf.4 (1), 6569. doi: 10.2165/00822942-200504010-00008

  • 47

    MuranakaT.OkadaM.YomoJ.KubotaS.OyamaT. (2015). Characterization of circadian rhythms of various duckweeds. Plant Biol.17, 6674. doi: 10.1111/plb.12202

  • 48

    MurrayM. G.ThompsonW. (1980). Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res.8 (19), 43214326. doi: 10.1093/nar/8.19.4321

  • 49

    NaumannB.EberiusM.AppenrothK. J. (2007). Growth rate based dose–response relationships and EC-values of ten heavy metals using the duckweed growth inhibition test (ISO 20079) with Lemna minor l. clone St. J. Plant Physiol.164 (12), 16561664. doi: 10.1016/j.jplph.2006.10.011

  • 50

    PagliusoD.GrandisA.FortirerJ. S.CamargoP.FlohE. I. S.BuckeridgeM. S. (2022). Duckweeds as promising food feedstocks globally. Agron.12 (4), 796. doi: 10.3390/agronomy12040796

  • 51

    PennockD. (2019). Soil erosion: The greatest challenge for sustainable soil management (Food and Agriculture Organization of the United Nations (FAO).

  • 52

    PetersenF.DemannJ.RestemeyerD.UlbrichA.OlfsH. W.WestendarpH.et al. (2021). Influence of the nitrate-n to ammonium-n ratio on relative growth rate and crude protein content in the duckweeds Lemna minor and Wolffiella hyalina. Plants10 (8), 1741. doi: 10.3390/plants10081741

  • 53

    Plants of the World Online (2022). Available at: https://powo.science.kew.org/ (Accessed 30/05/2022).

  • 54

    PoorterH. (1989). Plant growth analysis: towards a synthesis of the classical and the functional approach. Physiol. Plant75, 237244. doi: 10.1111/j.1399-3054.1989.tb06175.x

  • 55

    PotterJ. R.JonesJ. W. (1977). Leaf area partitioning as an important factor in growth. Plant Physiol.59 (1), 1014. doi: 10.1104/pp.59.1.10

  • 56

    RiedlH. (1976). Lemnaceae in K. H. Rechinger, Flora Iranica, no. 119. -Graz.

  • 57

    RozasJ.Ferrer-MataA.Sánchez-DelBarrioJ. C.Guirao-RicoS.LibradoP.Ramos-OnsinsS. E.et al. (2017). DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol.34 (12), 32993302. doi: 10.1093/molbev/msx248

  • 58

    Rubisco Foods (2022) Next generation ingredients. Available at: https://rubiscofoods.com/ (Accessed 06/05/2022).

  • 59

    SaranP. L.SinghS.SolankiV.ChoudharyR.ManivelP. (2021). Evaluation of asparagus adscendens accessions for root yield and shatavarin IV content in India. Turk J. Agric. For.45 (4), 475483. doi: 10.3906/tar-2006-42

  • 60

    SchenkR. U.HildebrandtA. C. (1972). Medium and techniques for induction and growth of monocotyledonous and dicotyledonous plant cell cultures. Can. J. Bot.50, 199204. doi: 10.1139/b72-026

  • 61

    ShawJ.LickeyE. B.BeckJ. T.FarmerS. B.LiuW.MillerJ.et al. (2005). The tortoise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. Am. J. Bot.92 (1), 142166. doi: 10.3732/ajb.92.1.142

  • 62

    SreeK. S.AppenrothK. J. (2014). Increase of starch accumulation in the duckweed Lemna minor under abiotic stress. Albanian J. Agric. Sci., 13, 1114.

  • 63

    SreeK. S.SudakaranS.AppenrothK. J. (2015). How fast can angiosperms grow? species and clonal diversity of growth rates in the genus Wolffia (Lemnaceae). Acta Physiol. Plant10 (37), 17. doi: 10.1007/s11738-015-1951-3

  • 64

    TipperyN. P.LesD. H.AppenrothK. J.SreeK. S.CrawfordD. J.BogM. (2021). Lemnaceae and orontiaceae are phylogenetically and morphologically distinct from araceae. Plants.10 (12), 2639. doi: 10.3390/plants10122639

  • 65

    TomlinsonK. W.PoorterL.BongersF.BorghettiF.JacobsL.van LangeveldeF. (2014). Relative growth rate variation of evergreen and deciduous savanna tree species is driven by different traits. Ann. Bot.114 (2), 315324. doi: 10.1093/aob/mcu107

  • 66

    ToyamaT.KurodaM.OgataY.HachiyaY.QuachA.TokuraK.et al. (2017). Enhanced biomass production of duckweeds by inoculating a plant growth-promoting bacterium, acinetobacter calcoaceticus P23, in sterile medium and non-sterile environmental waters. Water Sci. Technol.76 (6), 14181428. doi: 10.2166/wst.2017.296

  • 67

    TursiA. (2019). A review on biomass: importance, chemistry, classification, and conversion. Biofuel Res. J.6 (2), 962979. doi: 10.18331/BRJ2019.6.2.3

  • 68

    WalshÉ.KuehnholdH.O’BrienS.CoughlanN. E.JansenM. A. (2021). Light intensity alters the phytoremediation potential of Lemna minor. environ. Sci. pollut. Res.28 (13), 1639416407. doi: 10.1007/s11356-020-11792-y

  • 69

    WangW.WuY.YanY.ErmakovaM.KerstetterR.MessingJ. (2010). DNA Barcoding of the lemnaceae, a family of aquatic monocots. BMC Plant Biol.10 (1), 111. doi: 10.1186/1471-2229-10-205

  • 70

    XuJ.ChengJ. J.StompA. M. (2012). Growing Spirodela polyrrhiza in swine wastewater for the production of animal feed and fuel ethanol: a pilot study. Clean (Weinh).40 (7), 760765. doi: 10.1002/clen.201100108

  • 71

    XuS.StapleyJ.GablenzS.BoyerJ.AppenrothK. J.SreeK. S.et al. (2019). Low genetic variation is associated with low mutation rate in the giant duckweed. Nat. Commun.10 (1), 16. doi: 10.1038/s41467-019-09235-5

  • 72

    ZieglerP.AdelmannK.ZimmerS.SchmidtC.AppenrothK. J. (2015). Relative in vitro growth rates of duckweeds (Lemnaceae)–the most rapidly growing higher plants. Plant Biol.17 (s1), 3341. doi: 10.1111/plb.12184

Summary

Keywords

duckweed, Lemna, RGR, DNA barcoding, biomass accumulation, doubling time, biotechnology, food security

Citation

Taghipour E, Bog M, Frootan F, Shojaei S, Rad N, Arezoumandi M, Jafari M and Salmanian AH (2022) DNA barcoding and biomass accumulation rates of native Iranian duckweed species for biotechnological applications. Front. Plant Sci. 13:1034238. doi: 10.3389/fpls.2022.1034238

Received

01 September 2022

Accepted

01 November 2022

Published

29 November 2022

Volume

13 - 2022

Edited by

Beckley Ikhajiagbe, University of Benin, Nigeria

Reviewed by

Sezai Ercisli, Atatürk University, Turkey; David Dewez, Université du Québec à Montréal, Canada

Updates

Copyright

*Correspondence: Ali Hatef Salmanian,

This article was submitted to Crop and Product Physiology, a section of the journal Frontiers in Plant Science

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics