Next Article in Journal
Assessing Soil Key Fertility Attributes Using a Portable X-ray Fluorescence: A Simple Method to Overcome Matrix Effect
Next Article in Special Issue
Impact of Ionic Liquids on Induction of Wheat Microspore Embryogenesis and Plant Regeneration
Previous Article in Journal
Role of the Fungicide Seed Dressing in Controlling Seed-Borne Fusarium spp. Infection and in Enhancing the Early Development and Grain Yield of Maize
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Agronomy
Volume 10
Issue 6
10.3390/agronomy10060785
agronomy-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Colchicine Mutagenesis from Long-term Cultured Adventitious Roots Increases Biomass and Ginsenoside Production in Wild Ginseng (Panax ginseng Mayer)

by
Kim-Cuong Le
1,2,
Thanh-Tam Ho
3,4,
Jong-Du Lee
5,
Kee-Yoeup Paek
1,6 and
So-Young Park
1,*
1
Department of Horticulture, Division of Animal, Horticultural and Food Sciences, Chungbuk National University, Cheongju 28644, Korea
2
Department of Forest Genetics and Plant Physiology, Swedish University of Agricultural Sciences, 901 83 Umea, Sweden
3
Institute for Global Health Innovations, Duy Tan University, Da Nang 550000, Vietnam
4
Faculty of Pharmacy, Duy Tan University, Da Nang 550000, Vietnam
5
Biodiversity Research Institute, Shinyedong-ro 338, Seogwipo-si, Jeju 63608, Korea
6
WellGreen Co., Chungdae-ro 1, Cheongju 28644, Korea
*
Author to whom correspondence should be addressed.
Agronomy 2020, 10(6), 785; https://doi.org/10.3390/agronomy10060785
Submission received: 5 May 2020 / Revised: 22 May 2020 / Accepted: 29 May 2020 / Published: 31 May 2020
(This article belongs to the Special Issue Role of Plant Tissue Culture in Agricultural Research and Production)
Browse Figures
Review Reports Versions Notes

Abstract

:
Panax ginseng Mayer is a perennial herb that has been used as a medicinal plant in Eastern Asia for thousands of years. The aim of this study was to enhance root biomass and ginsenoside content in cultured adventitious roots by colchicine mutagenesis. Adventitious P. ginseng roots were treated with colchicine at different concentrations (100, 200, and 300 mg·L−1) and for different durations (1, 2, and 3 days). Genetic variability of mutant lines was assessed using random amplification of polymorphic DNA (RAPD) analysis. Ginsenoside biosynthesis gene expression, ginsenoside content, enzyme activities, and performance in bioreactor culture were assessed in four mutant lines (100–1-2, 100–1-18, 300–1-16, and 300–2-8). The results showed that ginsenoside productivity was enhanced in all mutant lines, with mutant 100–1-18 exhibiting the most pronounced increase (4.8-fold higher than the control). Expression of some ginsenoside biosynthetic enzymes was elevated in mutant lines. Enzyme activities varied among lines, and lipid peroxidation activity correlated with root biomass. All four lines were suitable for bioreactor cultivation, with mutant 100–1-18 exhibiting the highest biomass after culture scale-up. The results indicated that colchicine mutagenesis of P. ginseng roots increased biomass and ginsenosides production. This technique, and the root lines produced in this study, may be used to increase industrial yields of P. ginseng biomass and ginsenosides.

1. Introduction

Cultivation of plant cell tissues and organs is one of the most powerful tools available for research into plant propagation, crop development, production of valuable phytocompounds, and preservation of endangered plant species. However, in vitro plant materials accumulate genetic and epigenetic changes at high frequencies during long-term culture. Long periods of continuous subculture also enhance somaclonal variation [1], which leads to changes in phenotypes such as plant size, yield, and disease and insect tolerance in ornamentals [2,3], and changes to secondary metabolite production in medicinal plants [2,4,5,6,7]. Somaclonal variation can lead to critical economic losses when plant cell and organ culture is used for biomass propagation and uniform plant germination [1].
Panax ginseng Mayer (Araliaceae), known as ginseng, is one of the most important oriental herbs worldwide. Cell suspension, adventitious root, and hairy root culture systems were established for ginseng at Chungbuk National University around 2000 [8,9]. However, production of the primary bioactive compounds, ginsenosides (triterpenoid saponins), decreased after long-term continuous subculture of P. ginseng. Kiselev et al [7] showed that ginsenosides comprised only approximately 0.024% dry weight of a 20 years cell subculture. Similarly, Li et al [6] reported that paclitaxel was undetectable after five years of Taxus chinesesis cell subculture, and anthocyanin and alkaloid contents in Vitis vinifera and Catharanthus roseus, respectively, were gradually lost during long-term culture [4,5].
Several biotechnological technologies, such as primary culture, induction of polyploidy, and mutagenesis, can be used to overcome the difficulties encountered during long-term subculture. However, ginseng plants collected from their natural environment can be expensive, depending on their age, and ginseng, particularly wild ginseng, is rare in nature. Primary culture methods are thus inappropriate for such endangered or perennial plant materials. Colchicine is an alkaloid derivative from amino acids phenylalanine and tyrosine of Colchicum autumnale. The DNA breakage and clastogenicity of alkylating agents in plants has been studied for over 40 years [10]. Colchicine, which has been used for a long time to induce polyploidy in higher plants by inhibition of spindle fiber formation [11,12,13], is also an effective mutagen in many plant species [14,15,16,17,18,19,20]. Mutagenesis induction techniques combined with in vitro culture and molecular marker systems have facilitated the development of genetic variability platforms to produce enhanced reproductive and vegetative organs [17,21,22]. Significant mutagenic effects of colchicine result in the physiological disturbances and chromosome aberrations by alterations of the signaling pathway [17,20,23]. In medicinal plants, the selection of an efficient mutagenic agent is critical for generation of profitable and desired mutations at high frequency. Mutant induction using colchicine was previously used successfully in a range of plants, such as Hordeum vulgare, Helianthus annuus, Sesame indicum, Calendula officinalis, Peniesetum purpureum, Glycine max and Bacopa monnieri [15,16,17,18,19,20,24]. However, very few studies report the enhancement of secondary metabolite production in P. ginseng after mutagenesis with colchicine.
In the present study, colchicine mutagenesis of wild ginseng produced mutant adventitious roots with enhanced root biomass and ginsenoside content. DNA polymorphisms in the mutant lines were analyzed using random amplification of polymorphic DNA (RAPD) markers, and expression of genes critical for ginsenoside biosynthesis (PgSS, PgSE2, PPDS, and PPTS) was assessed. Mutant lines were also assessed for their culture characteristics in flask and bioreactor cultivation. Finally, for the first time a correlation between Malondialdehyde (MDA) activity/biomass and adaptation capability in mutant roots was established.

2. Materials and Methods

2.1. Plant Materials

Adventitious roots derived from calluses of 100-year-old wild ginseng (P. ginseng Mayer) roots were cultured in MS medium [25] supplemented with indole-3-butyric acid (IBA) and 3% (w/v) sucrose. The cultures were maintained at 24 ± 1 °C in the dark.

2.2. Mutant Root Induction with Colchicine

MS liquid medium supplemented with different concentrations (0, 100, 200, and 300 mg·L−1) of filtered-sterilized colchicine (Sigma-Aldrich, St. Louis, MO, USA) was utilized for mutant root induction. Ten adventitious roots (1.5 cm length) were utilized per replicate, with four replicates per treatment. Adventitious roots were placed in 100 mL flasks containing 20 mL of MS liquid medium and colchicine. The flasks were incubated on a rotary shaker (100 rpm) for 1, 2, and 3 days in darkness. Root tips, 1.5 cm long, were washed three times with MS medium and then transferred to fresh solid medium supplement with 5.0 mg·L−1 IBA, 3.0% sucrose and cultured in the dark at 24 ± 1 °C. Lateral root induction, the number of lateral roots per explant and lateral root length were recorded after eight weeks. The lateral roots of treated roots from each treatment were detached and cultured on the same medium describe above. The putative mutant roots were firstly screened for morphology and biomass after six weeks of culture. The four putative lines with the highest biomass were selected and screened again for ginsenoside content.
Mutant roots were transferred to a flask containing 100 mL of MS liquid medium containing IBA and 3% sucrose and assessed in a scaled-up 3 L bioreactor containing 2 L of MS liquid medium. The inoculation density of root in flask and bioreactor culture was 5.0 mg·L−1. Lateral root number, root length, and biomass (fresh and dry weight) were measured in the flask and bioreactor cultures. Mutant root morphology was also examined.

2.3. DNA Content Analysis by Flow Cytometry

DNA content of colchicine-treated adventitious roots was analyzed by flow cytometry (Partec PA, Münster, Germany) after six weeks of culture. Three to four lateral roots were chopped using a razor blade in 0.4 mL nuclei extraction buffer (CyStain UV Precise P, Partec, Münster, Germany) for 30 s. The sample was incubated for 5 min and filtered through a Partec 50 µm celltrics® disposable filter. Staining buffer (1.6 mL) was added to the sample tube, incubated for 30 s and analyzed in a flow cytometer. At least 3000–6000 nuclei were measured per sample and processed using FloMax software (Partec, Münster, Germany). The peak of nuclei insulated from a tetraploid (control) lateral root of P. ginseng was adjusted at channel 100.

2.4. RAPD Analysis

2.4.1. DNA Isolation

Fresh 6-week-old adventitious roots (0.1 g) from selected putative mutant lines and from three tetraploid control lines were homogenized using a Tissuelyser II (Qiagen, Hilden Germany). Genomic DNA was isolated following the cetyl trimethyl ammonium bromide (CTAB) extraction method [26]. CTAB buffer solution and 0.2% β-mercaptoethanol were added to the homogenized adventitious roots samples. Samples were incubated at 65 °C for 30 min, chloroform:isoamyl alcohol extraction solution (24:1) was added, and DNA was precipitated with isopropanol. Pellets were washed with 70% ethanol and resuspended in 50 µL of sterilized water. DNA concentration was measured using a DS-11 + spectrophotometer (Denovix, Inc., Wilmington, DE, USA).

2.4.2. PCR Amplification

Amplification was performed using a C1000TM Thermal Cycler (CFX96 Touch™ Real-time PCR Detection System, Bio-Rad Laboratories, Inc., Singapore) as described by Williams et al. (1990). Twenty oligonucleotides from Operon Inc (USA) were used for initial PCR amplifications, and six primers (A-05, A-08, A-09, A-10, A-13, and A-15) were subsequently selected for further amplifications (Table 1). PCR was conducted in a total volume of 20 µL containing DNA template (10 ng·µL−1), Primer Tag Premix (2×) (Genetbio Inc., Daejeon, Korea), primers, and sterilized water. RAPD amplification was performed over 35 cycles as follows: 94 °C for 1 min, 37 °C for 30 s, 72 °C for 1 min, and 72 °C for 5 min, after which the reactions were held at 4 °C. PCR products were analyzed using 2% agarose gel electrophoresis in 1× TEB buffer. Electrophoresed gels were examined with UV transillumination.

2.5. RNA Isolation and Quantitative Real-Time PCR Analysis of Gene Expression

Total RNA (100 mg) was isolated from mutant roots after 6 weeks of culture using NucleoZol reagent (MACHEREY-NAGEL GmbH& Co. KG, Düren, Germany). Total RNA was incubated at 65 °C for 5 min and reverse transcribed using ReverTra Ace® qPCR RT Master (TOYOBO CO, LTD, Osaka, Japan) in 10 µL reaction volumes containing 1 µg of total RNA and 5× RT Master (ReverTra Ace®, RNase inhibitor, oligo dT primer, random primer, MgCl2 and dNTPs). Reverse transcription reactions were performed at 37 °C for 15 min, and then 50 °C for 5 min. The mixture was then heated at 98 °C for 5 min followed by cooling at 4 °C to terminate the reaction. For quantitative real-time PCR (qPCR), specific primers were designed using the Primer3Plus program (http://carbon.bioneer.co.kr/primer3plus/). The primers used were P. ginseng ACT1 (actin, GenBank accession no. KF699319), P. ginseng SS2 (squalene synthase, GenBank accession no. GQ468527), P. ginseng SE2 (squalene epoxidase, GenBank accession no. FJ393274), P. ginseng PPDS (cytochrome P450 CYP716A47, GenBank accession no. JN604537), and P. ginseng PPTS (cytochrome P450 CYP716A53v2, GenBank accession no. JX036031) (Table S1). Gene expression levels of PgSS, PgSE2, PPDS, and PPTS in P. ginseng mutant roots were qPCR quantified using a C1000TM Thermal Cycler (CFX96 Touch™ Real-time PCR Detection System, Bio-Rad Laboratories, Inc., Singapore, Singapore) with the following cycle: 95 °C for 30 s; then 40 cycles of 95 °C for 5 s, 55 °C for 10 s, and 72 °C for 15 s; and a final 10 min extension at 72 °C.

2.6. Preparation and Extraction of Mutant Lines

The 6-week-old root samples (0.5 g dry weight) of control and selected mutant roots were placed in 100 mL flasks with 50 mL of 80% (v/v) ethanol. Flasks were sonicated in an ultra-sonication bath (SD-D250H by Mujigae Co., Seoul, Korea) for 1 h at room temperature. Thereafter, extracts were filtered through Whatman filter paper (No. 1002 110) and collected in a round-bottom flask. The solvent was evaporated to dryness using a rotary evaporator (N-1000, Eyela, Tokyo, Japan) at 40 °C, and the residue was dissolved in 50 mL of distilled water. The aqueous solution was washed twice with 50 mL of ethyl ether, and the aqueous layer was then extracted twice with 50 mL of water-saturated n-butanol. The n-butanol fraction was evaporated using a rotary evaporator at 50 °C. The sample solution was dissolved in 2 mL of methanol and filtered through a 0.2 µm Whatman syringe filter before analysis by High Performance Liquid Chromatography (HPLC).

2.7. HPLC Analysis of Ginsenosides

Extracted samples were analyzed by HPLC using a Waters 2695 separation module with a 2996 photodiode array detector on a Fortis C18 column (φ 5 µm, 150 × 4.6 mm). The mobile phase comprised acetonitrile (solvent A) and water (solvent B). The following gradient elution formula was used: 0 min, 18% A, 82% B; 0–42 min, 24% A, 76% B; 42–46 min, 29% A, 71% B; 46–75 min, 40% A, 60% B; 75–100 min, 65% A, 35% B; 100–135 min, 85% A, 15% B; and 135–150 min, 85% A, 15% B. A flow rate of 0.6 mL·min-1 was maintained throughout. Ginsenosides were detected at 203 nm. The sample injection volume was 20 µL, and the temperature of the column was controlled at 35 °C. Ginsenoside standards were purchased from ChromaDex (USA). Protopanaxadiol-type (Rb1, Rb2, Rb3, Rc, Rd, Rg3, Rh2) and protopanaxatriol-type (Re, Rf, Rg1, Rh1) standard solutions were used. Ginsenoside content was calculated as follows: Ginsenoside content (mg·g−1 DW) = ginsenoside concentration in sample (mg·L−1) × sample volume (L)/root dry weight (g). Total ginsenosides content was calculated as the sum of the ginsenoside fractions. Ginsenoside productivity was calculated as follows: Ginsenoside productivity (mg·L−1) = Total ginsenoside content (mg·g−1 DW) × harvested root dry weight (g)/volume of culture medium (L).

2.8. Lipid Peroxidation

Lipid peroxidation contents were measured using a modified malondialdehyde (MDA) method [27]. The 6-week-old root tissue (0.4 g) of control and selected mutant roots were disrupted in liquid nitrogen with a mortar and pestle and homogenized in 0.1% trichloroacetic acid (TCA). The homogenate was centrifuged at 5000 rpm for 5 min. The supernatant was then blended with 0.6% TBA and placed in a water bath at 90 °C for 30 min. Absorbance was measured at 532 and 600 nm using a spectrophotometer (Optizen POP, Mecasys Co., Ltd, Daejeon, Korea).

2.9. Catalase and Peroxidase Activity

Fresh 6-week-old adventitious root tissue (1.0 g) of control and mutant root lines (100–1-2, 100–1-18, 300–1-16, and 300–2-8) were disrupted in liquid nitrogen and homogenized in 4 mL of extraction buffer (50 mM potassium phosphate buffer, 2.0% polyvinylpolypyrrolidone, and 1.0 mM phenylmethylsulfonylfluoride). The extracted sample was centrifuged at 13,000 rpm for 10 min at 4 °C (Smart R17, Hanil science CO., LTD, Gimpo, Korea). The supernatant was stored at 2 °C and used for enzyme assays within 4 h.
Catalase, CAT (EC 1.11.1.6) activity was estimated in a reaction mixture containing 500 µmol H2O2 in 10 mL of 100 mM phosphate buffer (pH 7.0). CAT activity was determined by monitoring H2O2 consumption at 240 nm for 3 min. The results were expressed as CAT unit mg−1 protein (1 mM of H2O2 reduction min−1·mg−1 protein) [28].
Peroxidase, POD (EC 1.11.1.7) activity was measured using a modified method described by Bisht et al. [28]. The reaction mixture contained 1.5 mL of phosphate buffer (pH 7.0), 1% guaiacol, and 1% H2O2. POD activity was assessed as the change in absorption at 470 nm as a result of guaiacol oxidation. The results were expressed as unit min-1·mg-1 protein.

2.10. Statistics

RAPD data were scored for the absence “0” or presence “1” of electrophoresis band. Scoring data were entered into a binary matrix and analyzed using NTSYSpc 2.1. Similarities were calculated using Jaccard’s coefficient, and a dendrogram was established using Unweighted Pair Group Method with Arithmetic Mean (UPGMA). Data were analyzed using SPSS version 16.0 (SPSS Inc., Chicago, IL, USA). Where a significant difference (p < 0.05) was observed for a measured parameter, means were separated using Duncan’s multiple range test at the 5% level.

3. Results and Discussion

3.1. Effect of Colchicine Treatment on Lateral Root Induction and Growth

The induction of secondary roots after exposure to different concentrations of colchicine (100, 200, and 300 mg·L−1) for different durations (1, 2, and 3 days) was estimated eight weeks after treatment (Table 2; Figure 1). All the colchicine-treated adventitious root explants survived until eight weeks after treatment (data not shown). Colchicine treatment appreciably affected the number of explants with induced lateral roots (Table 2). Treatment slowed the induction and branching of lateral roots during the initial 2–4 weeks of culture (data not shown). However, there was no significant difference in lateral root formation between control and treated roots after eight weeks of cultivation (Table 2). It was possible that the colchicine-treated explants exhibited slowed root induction as a result of a physiological disturbance that reduced the cell division rate and caused initial growth retardation [29]. Colchicine-treated adventitious roots produced fewer, shorter lateral roots than untreated roots (Table 2). This may have been due to the loss of microtubules and the occurrence of sticky supercoiled chromosomes and c-mitoses. Loss of cortical microtubules likely resulted in cell expansion rather than cell elongation, producing shorter and thicker lateral roots after colchicine treatment (Figure 1C,D,G,H) [29,30]. Obute et al. [23], and Hewawasam et al. [31] observed similar suppression of growth in Vigna unguiculata, Cucumeropsis mannii, and Crossandra infundibuiformis. The presence of IBA in culture medium might have counteracted the initial harmful effects of colchicine [29]. New adventitious roots were detached, cultured, and screened by biomass. Four biomass-enhanced root lines were selected, and the histogram of nuclei (Figure. S1) and DNA index (Table 3) were assessed using flow cytometry, which showed that all four lines and controls were tetraploid. The capacity of colchicine to infiltrate inside cell of living organisms to interact with the DNA produces the typical toxic effects related to colchicine properties. Thereby, the mutagenic effects are principally caused by the direct interaction between the mutagen and DNA molecules [10,17].
The four putative lines with the highest biomass among the total of 254 lines after screened by morphology and biomass (data not showed) were selected and screened again for ginsenoside content. Putative mutant roots were propagated in 250 mL flasks containing 100 mL of MS medium supplemented with 5 mg·L−1 IBA and 3% sucrose, in the dark, on a rotary shaker (100 rpm) (Table 3). The putative mutant roots were initially screened for high root biomass after six weeks of culture. Further experiments using DNA marker techniques were used to elucidate the genetic backgrounds of the four selected lines.
Four mutant root lines (100–1-2, 100–1-18, 300–1-16, and 300–2-8) exhibited higher growth ratios and more rapid growth than the control. Lateral root number, root length, and fresh and dry biomass were also recorded in the four lines after six weeks of culture (Table 3). The mutant roots (Figure 1E,F,G,H) were thicker and formed more abundant and longer lateral roots than the control roots (Figure 1A,B,C,D). The biomass obtained after six weeks of culture was approximately 1.5-fold higher for mutant roots than for control roots (Table 3).

3.2. RAPD Analysis of Mutant Roots

RAPD analysis was used to assess genetic differences between the control and four putative mutant lines because it is a rapid detect method that does not require prior information regarding the nucleotide sequence, is inexpensive, and is suitable for detecting DNA alternations after exposure with mutagenic agents [32,33]. Six of eighteen RAPD primers (Table 1) produced different band patterns between control and mutant lines (data not shown). The principal differences seen in the RAPD profiles were the presence or absence of varying bands. Amplification data from the differential primers were scored and used to assess genetic distance (Figure 2A). Genetic distances between the control and the four mutant lines were 0.46–1.0, as assessed using Jaccard’s coefficient matrix (Figure 2B). UPGMA analysis of RAPD data assigned the control to one cluster and the four mutant lines to a separate cluster. Of the four mutant lines, 100–1-18 differed substantially from the other three lines. Although RAPD analysis is most commonly used for phytogenetic, taxonomic, and genetic mapping studies, RAPD can also be used for detection of DNA damage and mutation [34]. In the present study, RAPD markers were effective for analysis of P. ginseng mutant adventitious root lines.

3.3. Expression Analysis of Ginsenoside Biosynthetic Genes

The expression levels of four important genes involved in ginsenoside synthesis (PgSS, PgSE2, PPDS, and PPTS) were analyzed to better understand the molecular characteristics of the four mutant lines (Figure 3). Triterpene ginsenosides are principally biosynthesized through the mevalonic acid pathway in the cytoplasm and the methylerythriol phosphate pathway in the chloroplast [35,36]. Squalene synthase (PgSS) catalyzes the biosynthesis of triterpenes [37], and squalene epoxidase (PgSE2) is involved in the production of 2,3-oxidosqualene [38]. Cytochrome P450 enzymes are involved in a further stage of ginsenoside synthesis. PPDS (CYP716A47) catalyzes the formation of protopanaxadiol (PPD) from drammarenediol-II, and PPTS (CYP716A53v2) catalyzes the formation of protopanaxatriol (PPT) from PPD. Differences in PgSE2, PPDS, and PPTS expression were observed among mutant and control lines. PgSS transcription was elevated in the four mutants compared with the control. PgSS, PgSE2, PPDS, and PPTS were highly expressed in mutant 100–1-18 (Figure 3). The PgSE2 gene was also highly expressed in mutant 300–2-8, whereas the PgSE2, PPDS, and PPTS genes were minimally expressed in mutants 100–1-2 and 300–1-16 (Figure 3). The results clearly demonstrated that ginsenoside biosynthesis genes were expressed at higher levels in the mutant 100–1-18 than in the other mutant lines.

3.4. Ginsenoside Content Analysis by HPLC

Ginsenoside content was determined in the four mutant lines after six weeks of culture. The HPLC chromatogram with the retention time values of the determined ginsenosides and their standards are shown at Table S2 and Figure S2. The total contents of 11 ginsenosides increased in all mutant lines compared with the control (Table 4). Mutant 100–1-18 produced the highest amounts of total PPD and PPT ginsenosides (2.5-, 2.12-, and 2.68-fold higher than the control, respectively) (Table 4). In the PPD group, the Rb1, Rb2, and Rc contents of 100–1-18 were 3.7-, 4.1-, and 2.3-fold higher than the control, respectively (Table 4). In the PPT group, the Re, Rf, Rg1, and Rg2 contents of 100–1-18 were 2.87-, 4.7-, 1.7-, and 3.2-fold higher than the control, respectively (Table 4). The PPD/PPT ratio differed among mutant and control roots: the 100–1-2 roots exhibited the highest PPD/PPT ratio, and the 300–2-8 roots had the lowest ratio. Overall, mutant 100–1-18 had the highest ginsenoside productivity (186.6 mg·L−1), which was 4.85-fold higher than that of the control. Ginsenoside contents differed among the mutant and control lines as a result of disparate expression of ginsenoside synthesis genes. The high accumulation of ginsenosides in mutant 100–1-18 corresponded with elevated expression of PgSE2, PPDS, and PPTS (Figure 3; Table 4). The increased accumulation of ginsenoside content in mutant lines indicated that colchicine had a pronounced effect on ginsenoside synthesis. Mutants with elevated biomass and increased production of effective compounds are particularly valuable for medicinal plants [20,39]. Kharde et al [20] reported that bacoside production increased 4-fold in colchicine-treated in vitro Brahmi plants. In this study, mutant roots exhibited higher ginsenoside contents than control roots, demonstrating that this phytochemical characteristic was strongly affected by mutagenesis (Table 4). Mutant roots exhibited thicker and longer growth than control roots, which may have contributed to the elevated ginsenoside levels. Further studies are needed to fully understand the molecular mechanisms governing high ginsenoside content and productivity in P. ginseng. The ginsenoside content of mutant line 100–1-18 produced the highest PPD (Rb1, Rb2, Rg3, and Rh2) and PPT (Re, Rg1, and Rg2) levels. The 100–1-18 line is therefore of possible commercial pharmacological value.

3.5. MDA, CAT, and POD Activity of Mutant Roots

Antioxidase (POD and CAT) activities and MDA content were analyzed for clarify the distinctive characteristic in tetraploid and mutant roots. This allows a more detailed understanding of the molecular mechanism through which the mutagenesis effects on bioactive compounds stay obscure. MDA content was higher in mutant root lines than in the control line, and MDA contents were directly proportional to biomass which showed through the linear regression coefficients (r2) in the four mutant root lines (Figure 4). Mutant root lines 100–1-2 and 300–2-8 exhibited the lowest POD activities, whereas 100–1-18 and 300–1-16 activities were similar to those of the control. Mutant 100–1-2 also exhibited the lowest CAT value, whereas the other mutants resembled the control (Figure 5). No regression relationship among the POD and CAT actives and biomass was observed (data not shown). The colchicine-induced genetic alterations in the mutant lines may explain the differences in concentration and activities of some enzymes [40,41]. The difference in antioxidant enzyme activities and MDA content may be due to the upregulated expression of the corresponding genes.

3.6. Bioreactor Culture of Mutant Lines

The selected mutant roots were scaled-up in a bioreactor system. All four mutant lines grew successfully in the bioreactor system during the six-week culture duration (Table 3; Figure 6). Mutant 300–2-8, which showed the highest fresh biomass, exhibited particularly rapid adaptation capacity in the bioreactor system. After six weeks of culture, lateral root number and length were highest in mutant 100–1-18 (Table 3). Overall, all four mutant lines were suitable for scaled-up cultivation in the bioreactor system, suggesting their amenability to commercial production at industrial levels. Mutagenesis led to increased vegetative and reproductive organ production in P. ginseng. Similarly, Nura et al. [21] showed that leaf and seed yield were improved in sesame mutants. Colchicine induction of mutagenesis in sesame caused genomic alterations that enhanced cell division and meristematic expansion [21].

4. Conclusions

Mutagenesis is a powerful strategy for inducing desirable, valuable features such as secondary metabolite production in plants and is particularly applicable to ginseng. Here, in vitro mutant induction by colchicine in long-term adventitious root cultures of P. ginseng is demonstrated. Mutant induction of wild ginseng adventitious roots is essential for improving the biomass and ginsenoside content in the ginseng in vitro culture systems. Mutant roots characterized by rapid growth and high production of secondary metabolites may improve the efficiency of industrial production systems for foods, pharmaceutical products, and cosmetic compounds, among others.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4395/10/6/785/s1. Table S1: List of primers related to ginsenoside synthetic genes of Panax ginseng. Table S2. The retention time values of the ginsenoside standard of Panax ginseng by HPLC analysis. Figure S1: Flow cytometric histogram of nuclei of Panax ginseng.

Author Contributions

K.-C.L., T.-T.H. and J.-D.L. contributed to data acquisition and writing of the manuscript. K.-Y.P. participated in interpreting the data. S.-Y.P. made substantial contributions to the conception and design of this study, and of revising important intellectual content. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Korea Institute of Planning and Evaluation for Technology in Food, Agriculture, Forestry and Fisheries (IPET) through Advanced Production Technology Development Program, funded by Ministry of Agriculture, Food and Rural Affairs (Grant number 315013-4).

Acknowledgments

Authors were supported by Brain Korea (BK) 21 Plus Program through the National Research Foundation (NRF) of Korea.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bairu, M.W.; Kane, M.E. Physiological and developmental problems encountered by in vitro cultured plants. Plant Growth Regul. 2011, 63, 101–103. [Google Scholar] [CrossRef] [Green Version]
  2. Bairu, M.W.; Aremu, A.O.; Van Staden, J. Somaclonal variation in plants: Causes and detection methods. Plant Growth Regul. 2011, 63, 147–173. [Google Scholar] [CrossRef]
  3. Lee, H.J.; Kim, Y.E.; Yoon, Y.J.; Jeong, C.S.; Lian, M.L.; Paek, K.Y.; Park, S.Y. Highly endoreduplicated floral organs of somaclonal variants in clonally propagated Phalaenopsis “Spring Dancer”. Plant Cell Tissue Organ Cult. 2016, 126, 67–77. [Google Scholar] [CrossRef]
  4. Whitmer, S.; Canel, C.; van der Heijden, R.; Verpoorte, R. Long-term instability of alkaloid production by stably transformed cell lines of Catharanthus roseus. Plant Cell Tissue Organ Cult. 2003, 74, 73–80. [Google Scholar] [CrossRef]
  5. Qu, J.; Zhang, W.; Yu, X.; Jin, M. Instability of anthocyanin accumulation in Vitis vinifera L. var. Gamay Fréaux suspension cultures. Biotechnol. Bioprocess. Eng. 2005, 10, 155. [Google Scholar] [CrossRef]
  6. Li, L.Q.; Fu, C.H.; Zhao, C.F.; Xia, J.; Wu, W.J.; Yu, L.J. Efficient extraction of RNA and analysis of gene expression in a long-term Taxus cell culture using real-time RT-PCR. Z. Naturforsch C 2009, 64, 125–130. [Google Scholar] [CrossRef]
  7. Kiselev, K.V.; Dubrovina, A.S.; Shumakova, O.A. DNA mutagenesis in 2-and 20-yr-old Panax ginseng cell cultures. In Vitro Cell Dev. Biol. 2013, 49, 175–182. [Google Scholar] [CrossRef]
  8. Yu, K.W.; Gao, W.Y.; Son, S.H.; Paek, K.Y. Improvement of ginsenoside production by jasmonic acid and some other elicitors in hairy root culture of ginseng (Panax ginseng CA Meyer). In Vitro Cell Dev. Biol. 2000, 36, 424–428. [Google Scholar] [CrossRef]
  9. Lian, M.L.; Chakrabarty, D.; Paek, K.Y. Effect of plant growth regulators and medium composition on cell growth and saponin production during cell-suspension culture of mountain ginseng (Panax ginseng C. A. mayer). J. Plant Biol. 2002, 45, 201–206. [Google Scholar] [CrossRef]
  10. Leitão, J.M. Chemical mutagenesis. In Plant Mutation Breeding and Biotechnology; Shu, Q., Forster, B., Nakagawa, H., Eds.; CABI: Wallingford, UK, 2012; pp. 135–158. [Google Scholar]
  11. Davidson, D.; Pertens, E.; Zhao, J.P. Chromosome distribution between two restitution nuclei in a cell following colchicine treatment. Can. J. Genet. Cytol. 1983, 25, 437–445. [Google Scholar] [CrossRef]
  12. Ade, R.; Rai, M.K. Colchicine, current advances and future prospects. Nusant. Biosci. 2010, 2, 90–96. [Google Scholar]
  13. Naranjo, T. Contribution of structural chromosome mutants to the study of meiosis in plants. Cytogenet Genome Res. 2015, 147, 55–69. [Google Scholar] [CrossRef] [PubMed]
  14. Bragdo, M. Production of polyploids by colchicine. Euphytica 1955, 4, 76–82. [Google Scholar]
  15. Balkanjieva, J. Influence of genotype on mutagenic variability in barley (H. vulgare L.) following colchicine treatment. Barley Genet. Newsl. 1980, 10, 7–10. [Google Scholar]
  16. Downes, R.W.; Marshall, D.R. Colchicine-induced variants in sunflower. Euphytica 1983, 32, 757–766. [Google Scholar] [CrossRef]
  17. Mensah, J.K.; Obadoni, B.O.; Akomeah, P.A.; Ikhajiagbe, B.; Ajibolu, J. The effects of sodium azide and colchicine treatments on morphological and yield traits of sesame seed (Sesame indicum L.). Afr. J. Biotechnol. 2007, 6, 534–538. [Google Scholar]
  18. El-Nashar, Y.I.; Ammar, M.H. Mutagenic influences of colchicine on phenological and molecular diversity of Calendula officinalis L. Genet. Mol. Res. GMR 2016, 15. [Google Scholar] [CrossRef]
  19. Kamwean, P.; Chaisan, T.; Thobunluep, P.; Phumichai, C.; Bredemeier, M. Changing of morphological characteristic and biomass properties in Pennisetum purpureum by colchicine treatment. J. Agron. 2016, 16, 23–31. [Google Scholar] [CrossRef] [Green Version]
  20. Kharde, A.V.; Chavan, N.S.; Chandre, M.A.; Autade, R.H.; Khetmalas, M.B. In vitro enhancement of bacoside in brahmi (Bacopa monnieri) using colchicine. J. Plant Biochem. Physiol. 2017, 5, 1–6. [Google Scholar]
  21. Nura, S.; Adamu, A.K.; Mu’Azu, S.; Dangora, D.B.; Fagwalawa, L.D. Morphological characterization of colchicine-induced mutants in sesame (Sesamum indicum L.). J. Biol. Sci. 2013, 13, 277–282. [Google Scholar] [CrossRef] [Green Version]
  22. Ahloowalia, B.S.; Maluszynski, M. Induced mutations–A new paradigm in plant breeding. Euphytica 2001, 118, 167–173. [Google Scholar] [CrossRef]
  23. Obute, G.C.; Ndukwu, B.C.; Chukwu, O.F. Targeted mutagenesis in Vigna unguiculata (L.) Walp. and Cucumeropsis mannii (NAUD) in Nigeria. Afr. J. Biotechnol. 2007, 6, 2467–2472. [Google Scholar] [CrossRef] [Green Version]
  24. Shehu, A.S.; Adelanwa, M.A.; Alonge, S.O. Effects of gamma radiation and colchicine on crude protein content of soy-bean [Glycine max (L.) Merrill] varieties. Int. J. Res. 2016, 3, 1186–1193. [Google Scholar]
  25. Murashige, T.; Skoog, F. A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol. Plant 1962, 15, 473–497. [Google Scholar] [CrossRef]
  26. Doyle, J.J. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem. Bull. Bot. Soc. Am. 1987, 19, 11–15. [Google Scholar]
  27. Heath, R.L.; Packer, L. Photoperoxidation in isolated chloroplasts: I. Kinetics and stoichiometry of fatty acid peroxidation. Arch. Biochem. Biophys. 1968, 125, 189–198. [Google Scholar] [CrossRef]
  28. Bisht, S.S.; Sharma, A.; Chaturvedi, K. Certain metabolic lesions of chromium toxicity in radish. Indian J. Agric. Biochem. 1989, 2, 109–115. [Google Scholar]
  29. Roy, A.; Leggett, G.; Koutoulis, A. In vitro tetraploid induction and generation of tetraploids from mixoploids in hop (Humulus lupulus L.). Plant Cell Rep. 2001, 20, 489–495. [Google Scholar] [CrossRef]
  30. Dhooghe, E.; Van Laere, K.; Eeckhaut, T.; Leus, L.; Van Huylenbroeck, J. Mitotic chromosome doubling of plant tissues in vitro. Plant Cell Tissue Organ Cult. 2011, 104, 359–373. [Google Scholar] [CrossRef]
  31. Hewawasam, W.; Bandara, D.C.; Aberathne, W.M. New phenotypes of Crossandra infundibuliformis var. Danica through in vitro culture and induced mutations. Trop. Agric. Res. 2004, 16, 253–270. [Google Scholar]
  32. Atak, Ç.; Celik, O.; Acik, L. Genetic analysis of Rhododendron mutants using random amplified polymorphic DNA (RAPD). Pak. J. Bot. 2011, 43, 1173–1182. [Google Scholar]
  33. Zainudin, A.; Fitriani, H. Analysis of genetic diversity on mutants Jatropha curcas using RAPD. Energy Procedia 2014, 47, 1–6. [Google Scholar] [CrossRef] [Green Version]
  34. Williams, J.G.K.; Kubelik, A.R.; Livak, K.J.; Rafalski, J.A.; Tingey, S.V. DNA polymorphisms amplified by arbitrary primers are useful as genetic markers. Nucleic Acids Res. 1990, 18, 6531–6535. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Murthy, H.N.; Georgiev, M.I.; Kim, Y.S.; Jeong, C.S.; Kim, S.J.; Park, S.Y.; Paek, K.Y. Ginsenosides: Prospective for sustainable biotechnological production. Appl. Microbiol. Biotechnol. 2014, 98, 6243–6254. [Google Scholar] [CrossRef]
  36. Kim, Y.J.; Zhang, D.; Yang, D.C. Biosynthesis and biotechnological production of ginsenosides. Biotechnol. Adv. 2015, 33, 717–735. [Google Scholar] [CrossRef]
  37. Abe, I.; Rohmer, M.; Prestwich, G.D. Enzymatic cyclization of squalene and oxidosqualene to sterols and triterpenes. Chem. Rev. 1993, 93, 2189–2206. [Google Scholar] [CrossRef]
  38. Vincken, J.P.; Heng, L.; de Groot, A.; Gruppen, H. Saponins, classification and occurrence in the plant kingdom. Phytochemistry 2007, 68, 275–297. [Google Scholar] [CrossRef]
  39. Mehreen, K. Assessment of biochemical attributes of Praecitrullus fistulosus treated with mutagens. J. Plant Breed. Genet. 2015, 3, 39–47. [Google Scholar]
  40. Caruso, I.; Lepore, L.; De Tommasi, N.; Dal Piaz, F.; Frusciante, L.; Aversano, R.; Garramone, R.; Carputo, D. Secondary metabolite profile in induced tetraploids of wild Solanum commersonii Dun. Chem. Biodivers. 2011, 8, 2226–2237. [Google Scholar] [CrossRef]
  41. Xu, C.; Tang, T.; Chen, R.; Liang, C.; Liu, X.; Wu, C.; Yang, Y.S.; Yang, D.P.; Wu, H. A comparative study of bioactive secondary metabolite production in diploid and tetraploid Echinacea purpurea (L.) Moench. Plant Cell Tissue Organ Cult. 2014, 116, 323–332. [Google Scholar] [CrossRef]
Agronomy 10 00785 g001 550
Figure 1. Morphology (A,B,E,F) and histology (C,D, G,H) of control and mutant root (100–1-18) in Panax ginseng after six weeks of culture. Adventitious root of control (A) and mutant 100–1-18 (E). Close-up adventitious root of control (B) and mutant 100–1-18 (F). Cross section of root (C,G), and longitudinal section of primordia (D,H) of control (C,D) and mutant root (100–1-18) (G,H) in Panax ginseng.
Figure 1. Morphology (A,B,E,F) and histology (C,D, G,H) of control and mutant root (100–1-18) in Panax ginseng after six weeks of culture. Adventitious root of control (A) and mutant 100–1-18 (E). Close-up adventitious root of control (B) and mutant 100–1-18 (F). Cross section of root (C,G), and longitudinal section of primordia (D,H) of control (C,D) and mutant root (100–1-18) (G,H) in Panax ginseng.
Agronomy 10 00785 g001
Agronomy 10 00785 g002 550
Figure 2. RAPD amplification (A) and similarity coefficient dendrogram (B) of mutant adventitious roots of Panax ginseng. Primers A1: A-05, A2:A10, and A3:A-15),. M, marker; lane 1, control; lane 2, 100–1-2; lane 3, 100–1-18; lane 4, 300–1-16; and lane 5, 300–2-8.
Figure 2. RAPD amplification (A) and similarity coefficient dendrogram (B) of mutant adventitious roots of Panax ginseng. Primers A1: A-05, A2:A10, and A3:A-15),. M, marker; lane 1, control; lane 2, 100–1-2; lane 3, 100–1-18; lane 4, 300–1-16; and lane 5, 300–2-8.
Agronomy 10 00785 g002
Agronomy 10 00785 g003 550
Figure 3. Gene expression levels of four mutant Panax ginseng root lines. SS, squalene synthase; SE, squalene epoxidase; PPDS, cytochrome P450 CYP716A47; and PPTS, cytochrome P450 CYP716A53v2.
Figure 3. Gene expression levels of four mutant Panax ginseng root lines. SS, squalene synthase; SE, squalene epoxidase; PPDS, cytochrome P450 CYP716A47; and PPTS, cytochrome P450 CYP716A53v2.
Agronomy 10 00785 g003
Agronomy 10 00785 g004 550
Figure 4. Malondialdehyde (MDA) content (A) and correlation between MDA content and fresh weight (B) in four mutant lines of Panax ginseng after six weeks of culture. CT, control. Data present the mean ± SE of three replications. z Different letters indicate significant difference at p < 0.05 according to Duncan’s multiple range test (n = 3).
Figure 4. Malondialdehyde (MDA) content (A) and correlation between MDA content and fresh weight (B) in four mutant lines of Panax ginseng after six weeks of culture. CT, control. Data present the mean ± SE of three replications. z Different letters indicate significant difference at p < 0.05 according to Duncan’s multiple range test (n = 3).
Agronomy 10 00785 g004
Agronomy 10 00785 g005 550
Figure 5. Peroxidase (POD) (A) and catalase (CAT) (B) activities in four mutant lines of Panax ginseng after six weeks of culture. CT, control. Data present the mean ± SE of three replications. z Different letters indicate significant difference at p < 0.05 according to Duncan’s multiple range test (n = 3).
Figure 5. Peroxidase (POD) (A) and catalase (CAT) (B) activities in four mutant lines of Panax ginseng after six weeks of culture. CT, control. Data present the mean ± SE of three replications. z Different letters indicate significant difference at p < 0.05 according to Duncan’s multiple range test (n = 3).
Agronomy 10 00785 g005
Agronomy 10 00785 g006 550
Figure 6. Scaled-up cultivation of four mutant lines in the bioreactor system.
Figure 6. Scaled-up cultivation of four mutant lines in the bioreactor system.
Agronomy 10 00785 g006
Table
Table 1. Nucleotide sequences of RAPD primers that detected polymorphism. RAPD: random amplification of polymorphic DNA.
Table 1. Nucleotide sequences of RAPD primers that detected polymorphism. RAPD: random amplification of polymorphic DNA.
PrimersSequence
A-055ʹ-AGGGGTCTTG-3ʹ
A-085ʹ-GTGACGTAGG-3ʹ
A-095ʹ-GGGTAACGCC-3ʹ
A-105ʹ-GTGATCGCAG-3ʹ
A-135ʹ-CAGCACCCAC-3ʹ
A-155ʹ-TTCCGAACCC-3ʹ
Table
Table 2. Effect of colchicine treatments on DNA content and growth in adventitious root of Panax ginseng.
Table 2. Effect of colchicine treatments on DNA content and growth in adventitious root of Panax ginseng.
Colchicine Concentration
(mg·L−1)		Exposure
Time
(day)		Lateral Root Formation
(%)		No. of Lateral Roots/ExplantLateral Root Length
(cm)
00100.0a z27.75a2.51a
1001100.0a25.75b2.24c
2100.0a20.13c1.79e
3100.0a12.25g1.50g
2001100.0a18.25d2.31b
2100.0a15.00e2.01d
397.5ab11.13h1.60f
3001100.0a17.50d2.34b
297.5ab13.38f2.05d
395.0b10.75h2.01d
z Different letters within a column indicate significant difference at p < 0.05 according to Duncan’s multiple range test (n = 4).
Table
Table 3. Characteristics of mutant roots of Panax ginseng after six weeks of flask and bioreactor culture.
Table 3. Characteristics of mutant roots of Panax ginseng after six weeks of flask and bioreactor culture.
Culture ScaleLinesDNA Index zNo. of Lateral Roots/ExplantFresh Mass
(g·L−1)		Dry Mass
(g·L−1)		Growth
Ratio y
Flask (250 mL)Control1.00a24.00c x70.08b4.94b6.01b
100–1-21.08a32.00bc111.32a8.69a10.14a
100–1-180.99a38.25ab116.27a9.95a10.63a
300–1-161.00a28.50bc118.60a8.14ab10.86a
300–2-80.96a48.00a110.71a8.45a10.07a
Bioreactor (3 L)Control1.00a25.83c73.36d5.39c13.91c
100–1-21.08a39.83b84.31c6.25bc16.10b
100–1-180.99a50.50a94.87b7.19ab17.98a
300–1-161.00a39.50b93.28b7.60ab17.66a
300–2-80.96a38.50b99.84a6.87ab18.98a
z DNA index = (DNA value of mutant line G0/G1)/(DNA value of control G0/G1). y Growth ratio = [Harvested dry weight (g) − Inoculated dry weight (g)]/Inoculated dry weight (g). x Different letters within a column indicate significant difference at p < 0.05 according to Duncan’s multiple range test (n = 4).
Table
Table 4. Ginsenoside profiling of mutant Panax ginseng roots by HPLC analysis.
Table 4. Ginsenoside profiling of mutant Panax ginseng roots by HPLC analysis.
LinePPD (mg·g−1 DW)PPT (mg·g−1 DW)PPD (mg·g−1 DW)PPT (mg·g−1 DW)PPD/PPTTotal Ginsenoside (mg·g−1 DW) zGinsenoside Productivity (mg·L−1)
Rb1Rb2Rb3RcRdRg3Rh2ReRfRg1Rg2
Control1.08d y0.29c1.23a0.39b0.24a0.06a0.52a1.64c0.07b0.62b0.22b3.71b2.54c1.46b6.25b38.45c
100–1-22.98b0.51b0.15b0.59b0.46a0.06ab0.53a1.13c0.09b0.53b0.18b5.08b1.93c2.63a7.01b77.21b
100–1-184.01a1.19a0.28b0.88a0.50a0.05a0.48a4.72a0.33a1.06a0.71a7.85a6.81a1.15c14.66a186.60a
300–1-162.18bc0.52b0.14b0.35b0.30a0.05ab0.44a1.60c0.05b0.64b0.25b4.22b2.54c1.66b6.76b72.60b
300–2-81.74cd0.45bc0.10b0.29b0.31a0.10b0.44a2.29b0.12b1.16a0.30b3.56b3.87b0.92d7.43b71.83b
z Total content= (Rb1+Rb2+Rb3+Rc+Rd+Rg3+Rh2+Re+Rf+Rg1+Rg2). y Different letters within a column indicate significant difference at p < 0.05 according to Duncan’s multiple range test (n = 3).

Share and Cite

MDPI and ACS Style

Le, K.-C.; Ho, T.-T.; Lee, J.-D.; Paek, K.-Y.; Park, S.-Y. Colchicine Mutagenesis from Long-term Cultured Adventitious Roots Increases Biomass and Ginsenoside Production in Wild Ginseng (Panax ginseng Mayer). Agronomy 2020, 10, 785. https://doi.org/10.3390/agronomy10060785

AMA Style

Le K-C, Ho T-T, Lee J-D, Paek K-Y, Park S-Y. Colchicine Mutagenesis from Long-term Cultured Adventitious Roots Increases Biomass and Ginsenoside Production in Wild Ginseng (Panax ginseng Mayer). Agronomy. 2020; 10(6):785. https://doi.org/10.3390/agronomy10060785

Chicago/Turabian Style

Le, Kim-Cuong, Thanh-Tam Ho, Jong-Du Lee, Kee-Yoeup Paek, and So-Young Park. 2020. "Colchicine Mutagenesis from Long-term Cultured Adventitious Roots Increases Biomass and Ginsenoside Production in Wild Ginseng (Panax ginseng Mayer)" Agronomy 10, no. 6: 785. https://doi.org/10.3390/agronomy10060785

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

Article Metrics

Back to TopTop