Next Article in Journal
Biofumigation for Fighting Replant Disease- A Review
Next Article in Special Issue
Legume Biofortification and the Role of Plant Growth-Promoting Bacteria in a Sustainable Agricultural Era
Previous Article in Journal
Comparison of Greenhouse Energy Requirements for Rose Cultivation in Europe and North Africa
Previous Article in Special Issue
Soil Type and Zinc Doses in Agronomic Biofortification of Lettuce Genotypes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gene Expression and Metabolite Profiling of Thirteen Nigerian Cassava Landraces to Elucidate Starch and Carotenoid Composition

1
Swedish University of Agricultural Sciences, Sundsvägen 10, SE-230 53 Alnarp, Sweden
2
International Institute of Tropical Agriculture, PMB 5320, Oyo Road, Ibadan 200001, Oyo State, Nigeria
3
Leibniz Institute of Vegetable and Ornamental Crops, Plant Quality and Food Security, Theodor-Echtermeyer-Weg 1, 14979 Großbeeren, Germany
4
Department of Food Chemistry, Institute of Nutritional Science, University of Potsdam, Arthur-Scheunert-Allee 114-116, 14558 Nuthetal, Germany
5
Institute for Sustainable Plant Protection, CNR, via Amendola 122, I-70126 Bari, Italy
*
Author to whom correspondence should be addressed.
Agronomy 2020, 10(3), 424; https://doi.org/10.3390/agronomy10030424
Submission received: 5 February 2020 / Revised: 5 March 2020 / Accepted: 10 March 2020 / Published: 20 March 2020
(This article belongs to the Special Issue Biofortification of Crops)

Abstract

:
The prevalence of vitamin A deficiency in sub-Saharan Africa necessitates effective approaches to improve provitamin A content of major staple crops. Cassava holds much promise for food security in sub-Saharan Africa, but a negative correlation between β-carotene, a provitamin A carotenoid, and dry matter content has been reported, which poses a challenge to cassava biofortification by conventional breeding. To identify suitable material for genetic transformation in tissue culture with the overall aim to increase β-carotene and maintain starch content as well as better understand carotenoid composition, root and leaf tissues from thirteen field-grown cassava landraces were analyzed for agronomic traits, carotenoid, chlorophyll, and starch content. The expression of five genes related to carotenoid biosynthesis were determined in selected landraces. Analysis revealed a weak negative correlation between starch and β-carotene content, whereas there was a strong positive correlation between root yield and many carotenoids including β-carotene. Carotenoid synthesis genes were expressed in both white and yellow cassava roots, but phytoene synthase 2 (PSY2), lycopene-ε-cyclase (LCYε), and β-carotenoid hydroxylase (CHYβ) expression were generally higher in yellow roots. This study identified lines with reasonably high content of starch and β-carotene that could be candidates for biofortification by further breeding or plant biotechnological means.

1. Introduction

Cassava (Manihot esculenta Crantz) is an important starchy root crop in the tropics. It is a common staple food in Africa, Asia, and Latin America [1], and particularly of high importance in sub-Saharan African countries. Cassava is cultivated for its tolerance to stress and ability to grow under drought conditions [2]. As a result, cassava serves as an important component for food security for low income farmers who rely on it as a major source of dietary energy. While cassava is rich in starch and serves as a good energy source, it is extremely low in protein and important micronutrients, such as iron, zinc, and provitamin A carotenoids [3]. Therefore, a diet reliant on cassava predisposes one to malnutrition, especially provitamin A deficiency.
Carotenoids are C40 lipophilic isoprenoids synthesized in plants, algae, bacteria, and fungi [4]. Carotenoids play essential roles in photosynthesis and photoprotection and provide precursors for phytohormones abscisic acid (ABA) and strigolactones (SLs) [5]. They are also precursors of provitamin A, an essential component of human diet due to its roles in human health [6].
Carotenoids are synthesized in plants via the methylerythritol 4-phosphate (MEP) pathway, localized in the plastids [4] (Figure 1). Carotenoid synthesis starts with the formation of the C5 prenyl phosphates; isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP) are formed when isomerized by isopentenyl diphosphate isomerase (IPI) [6]. IPP and DMAPP then undergo reaction condensation by geranylgeranyl diphosphate synthase (GGPPS) to form geranylgeranyl diphosphate (GGPP) [4]. Condensation of two molecules of GGPP by phytoene synthase (PSY) leads to the formation of phytoene, the first carotenoid product and known rate-limiting metabolite of carotenoid biosynthesis [4,5]. Plant tissues exhibiting low levels of carotenoid production are thought to have a low expression level of PSY; additionally, PSY expression driven to high levels has resulted in accumulation of β-carotene crystals in Arabidopsis and carrot [3]. An overexpression of PSY in nonphotosynthetic storage roots such as cassava (Manihot esculenta) and some nongreen plants like tomato fruits (Solanum lycopersicum), canola seeds (Brassica napus), and rice endosperm (Oryza sativa) has also resulted in an increased flux through the biosynthetic pathway, further confirming that high carotenoid levels depend on PSY expression [3,7]. Phytoene is converted into lycopene by a series of desaturation and isomerization reactions catalyzed by phytoene desaturase (PDS), ζ-carotene desaturase (ZDS), ζ-carotene isomerase (ZISO), and carotenoid isomerase (CRTISO) [4,5,6]. Lycopene undergoes a cyclization reaction by lycopene ε-cyclase (LCYε) and/or lycopene β-cyclase (LCYβ) to produce α-carotene and β-carotene, representing the α- and β-branch of the pathway, respectively [5]. Subsequently, α-carotene and β-carotene undergo hydroxylation by β-carotenoid hydroxylases (CHYβ) to form lutein in the α-branch and zeaxanthin in the β-branch of the pathway [4,5,6]. Epoxidation and de-epoxidation of zeaxanthin by zeaxanthin epoxidase (ZEP) and violaxanthin de-epoxidase (VDE) form the xanthophyll cycle, a mechanism to protect plants against photodamage [5]. The pathway concludes by the conversion of violaxanthin into neoxanthin by neoxanthin synthase (NXS). Furthermore, oxidative cleavage of carotenoids by carotenoid cleavage dioxygenases (CCDs) such as 9-cis-epoxycarotenoid dioxygenase (NCED), produces apocarotenoids, phytohormones, signaling molecules and volatiles, thus maintaining a steady-state level of carotenoids and curbing excessive accumulation [5]. Figure 1 shows a simplified carotenoid pathway.
Food high in β-carotene is related to several health benefits because of its contribution to vitamin A in the body. Other carotenoids, such as β-cryptoxanthin and α-carotene can also be converted to vitamin A following cleavage by β-carotene oxygenase 1 (BCO1) [8]. Vitamin A maintains good eyesight, keeps the skin and epithelial linings of the internal organs healthy, and maintains immunity against diseases [8]. It is also involved in growth, development, and reproduction [9]. Vitamin A deficiency (VAD) is a global health problem affecting mostly vulnerable groups, such as preschool children and pregnant women [10], and health consequences include a weakened immune system, poor eyesight, and eventual blindness. Globally, 1–2.5 million child deaths are attributed to nutritional deficiencies; an estimated 761,000 of these deaths are caused by vitamin A, iron, and zinc deficiency [11]. Due to these consequences, efforts are being made to improve the provitamin A carotenoid content of major staple crops such as cassava. These efforts include fortification, dietary diversification, supplementation, and biofortification of major staple crops such as rice and cassava. Of all these efforts, biofortification proves to be most sustainable [12].
Biofortification is the use of breeding or genetic engineering to enhance nutritional content of crops and subsequently improve human nutritional value [7]. Genetic engineering of cassava has the advantage in that conventional cassava breeding requires a lengthy process to breed for new traits, usually about 8–13 years [13]. Studies have reported a negative correlation between dry matter and carotenoid content due to genetic linkage in African cassava germplasm [2,13,14]. The observed negative correlation may pose a challenge to biofortification efforts for cassava.
Carotenoid content and expression of key transcript in carotenoid biosynthesis were studied previously in cassava landraces of storage roots with varying color. Analysis of 19 carotenoid types by high performance liquid chromatography with diode array detection (HPLC–DAD) detected a complex carotenoid composition as well as a significant association between lycopene β-cyclase (LYCβ) expression and total β-carotene content [15].
Our study compares the variations in carotenoid profiles of thirteen Nigerian cassava landraces as well as the expression levels of carotenoid biosynthesis genes in a select subset of these landraces. These landraces are screened in an effort to introduce more cassava lines to tissue culture, thereby widening the genetic base with desirable traits such as high dry matter, high starch and high micronutrient content, for further improvement by genetic transformation or gene editing. The aim of this study is also to better understand the relationship between agronomic traits, carotenoids, and expression of genes in the carotenoid biosynthesis pathway in cassava.

2. Materials and Methods

2.1. Plant Material

Thirteen cassava landraces were selected from the collection of the International Institute of Tropical Agriculture (IITA), Nigeria, based on variance in starch and β-carotenoid accumulation.

2.2. Field Trial and Layout

Trials were conducted over two seasons between July 2017 and August 2019, each lasting for 13 months. Trials for the two seasons were carried out at the same location at IITA Ibadan, Nigeria (N 7.490250,E 3.884143). Temperature and rainfall data were collected from the IITA weather station over the two-season period to monitor the potential impact of seasonal changes of weather conditions on detected trait variations. For the agronomic measurements, trials were laid out in a completely randomized block design with two replications. There were twenty plants per plot, each plot had four rows with five plants in each row; two replications provided forty plants per genotype. Plants from both replication plots were harvested both seasons to confirm agronomic parameters and performance of all genotypes in the field.
Planting was done with a 1 m spacing between each plant and adjacent plots were separated by 1 m alleys. All the stakes used for planting were generated from fresh portions of mature stems. Hand weeding was done when needed. Since both field trials started in July (during the rainy season), plants received adequate water from rainfall during the initial growth stage, and there was no need for irrigation.

2.3. Sample Collection for Gene Expression and Metabolite Analyses

Source leaves and storage roots were collected for gene expression and metabolite analyses. Leaves and roots were immediately snap frozen in liquid nitrogen. Root samples were collected and prepared under a canopy tent to shield them from sunlight as much as possible. For molecular and metabolic analysis, three roots were randomly selected from the three plants harvested from the same plot to ensure homogeneity and minimize differences in soil characteristics. Selected roots were peeled, cleaned, and portions of the central part were chopped into small cubes and mixed together. The roots were then placed in aluminum foil and teabags, snap frozen in liquid nitrogen and stored at −80 °C.

2.4. Agronomic Data and Dry Matter Content (DMC)

For each trial, dry matter content (DMC) and fresh root weight (FRW) were measured at 4, 8, and 13 months after planting. Three plants per genotype were uprooted at random from the same plot at every sampling time to measure traits. iCheck Carotene by BioAnalyt (Teltow, Germany) a portable spectrophotometer, was used to measure total carotenoids. This tool determines the total carotenoid concentration based on absorbance at 450 and 525 nm [14]. At harvest, biomass from plants was used to estimate yield by weighing fresh roots and foliage separately. Harvest index (HI) was computed following the formula by Esuma et al. [2], for the calculation of harvest index: HI = FRW/(FRW + FSW), where HI is harvest index, FRW is fresh root weight, and FSW is fresh shoot weight. To measure dry matter content, roots were peeled, cleaned, and chopped. Three hundred grams of homogeneously mixed roots were collected into paper bags and oven dried at 60−80 °C for a minimum of 48 h. Samples were weighed after drying and the percentage DMC calculated as follows: DMC = (DSW/FSW) × 100 (where DSW is dry sample weight and FSW is fresh sample weight).

2.5. Metabolite Analysis

At eight months after planting, samples were collected for metabolite analysis. Thirty grams of root cubes were snap frozen and lyophilized prior to analysis. Lyophilized root samples were analyzed for starch content using a Total Starch Assay kit by Megazyme (Bray, Ireland). The standard protocol, which excludes D-glucose rinsing and measurement of the resistant starch of samples, was followed. Using a glucose standard, absorbance was measured at 510 nm and starch values were calculated using the weight and the absorbance value. Carotenoid content and composition were measured for the same samples by ultra-high performance liquid chromatography coupled with diode array detection and time-of-flight mass spectrometry (UHPLC–DAD–ToF–MS) at the Leibniz Institute of Vegetable and Ornamental Crops, Plant Quality and Food Security, Germany. Extraction and analysis of carotenoids and chlorophyll were conducted according to Frede et al. [16] using an Agilent Technologies 1290 Infinity II UHPLC (Agilent Technologies Sales and Services GmbH and Co. KG, Waldbronn, Germany) coupled to Agilent Technologies 6230 TOF LC/MS.

2.6. DNA and RNA Extraction and Amplification

Total nucleic acids (TNA) was extracted from cassava roots using a modified cetyltrimethylammonium bromide (CTAB) extraction protocol [16], followed by a clean-up step using a DNA clean and concentrator kit (ZYMO Research, Irvine, CA USA). Total RNA was extracted in biological and technical triplicates from TNA using RNA clean and concentrator with DNase I (ZYMO research, Irvine, CA USA) treatment to eliminate genomic DNA. The quantity of RNA was assessed using a NANODROP 8000 spectrophotometer (Thermo Scientific, Waltham, MA USA) and quality was further assessed by agarose gel electrophoresis. RNA was reverse-transcribed into cDNA using M-MuLV Reverse Transcriptase (New England Biolabs, Ipswich, MA United States) according to manufacturer’s instructions. The cDNA products were then used in quantitative real-time polymerase chain reaction (qPCR) experiment.
Using the Applied Biosystems Veriti Thermal Cycler, cDNA and gDNA samples were amplified with Phusion proofreading DNA polymerase (Thermo Scientific) and primers complementary to the selected carotenoid genes, CHYβ, LCYε, and NCED1 under the following thermal cycling conditions: 94 °C for 30 s (1 cycle), 94 °C for 10 s, 60 °C for 30 s, 72 °C for 1.45 min (35 cycles), 72 °C for 10 min. The qPCR products were run on 2% agarose gel with 1× Tris base, Acetic acid and EDTA (TAE) buffer (Sigma-Aldrich, St. Louis, MO, USA), at 100 V for 2 h to confirm PCR products. To confirm product sizes match those of the selected genes, a 100 bp DNA ladder (New England Biolabs) was used.

2.7. Primer Design and RT–qPCR Conditions

The expression levels of carotenoid genes, phytoene synthase 1 and 2 (PSY1 and 2), β-carotenoid hydroxylase (CHYβ), lycopene-ε-cyclase (LCYε), and 9-cis-epoxycarotenoid dioxygenase (NCED1) were investigated in roots of selected genotypes by qPCR. Using a standard curve, four reference genes: Actin, protein phosphatase 2 (PP2A), ubiquitin (UBQ), and G/T binding protein (GTBP), were tested across all samples. Since Actin gave a higher and relatively stable expression than the other reference genes, it was selected as the reference gene for the qPCR experiment. The primers for RT–PCR were designed using Integrated DNA Technologies (IDT) primer design tool (PrimerQuest® Tool, Integrated DNA Technologies, Coralville, Iowa, USA), using the following parameters: annealing temperature of 60–65 °C, GC content of 40–60%, and amplicon size of 50–150 bp (Table S1). The qPCR was performed in three biological and three technical replicates in a Roche LightCycler 480 II Instrument (Roche Diagnostics, Indianapolis, IN, USA) using the Luna Universal qPCR Master Mix according to the manufacturer’s instructions (New England Biolabs). Primer specificity was confirmed by the presence of a single melting peak in the dissociation curves. Each PCR reaction mix was prepared according to manufacturer’s instructions. PCR cycling was performed as follows: 5 min at 94 °C followed by 40 rounds of 15 s at 94 °C, 10 s at 60 °C, 15 s at 72 °C, and finally 1 round of 35 s at 60 °C. Melting curve cycling consisted of: 15 s at 95 °C, 1 min at 60 °C, 30 s at 95° C, and 15 s at 60 °C. All primers used in the RT–qPCR reaction had an efficiency greater than 90%. The comparative ΔCt method was used to determine the standard curve Actin as reference gene for relative quantification. The relative expression data were calculated according to the Livak (2-ΔΔCt) method [17]. Primer sequences are included in Table S1.

2.8. Statistical Analysis

To distinguish groups of common genotypes, ANOVA tests were performed on β-carotene measurements from the UHPLC–DAD–ToF–MS analysis of the root tissues, percent starch measurements, and dry matter content from both 2018 and 2019 seasons using R. Tukey groupings were then generated using the agricolae library [18].
To further understand the results presented by UHPLC–DAD–ToF–MS analyses, principal component analysis was done to show the patterns of variation of metabolites. Principal component analyses (PCAs) were constructed for leaf and root data. For each test, data were first normalized in the following way: (1) for each sample, the total sum of all metabolites was found, (2) the average of all the samples’ sums was calculated, (3) a scaling ratio for each sample was determined by dividing the average of all sums by each sample’s total sum, and (4) all of the metabolites in each sample were multiplied by that sample’s scaling ratio. This normalization method scaled measurements up if the total amount of metabolites found in the sample was less than average, and vice versa. These datasets were then input into the prcomp function in R. Visualizations were created using the ggbiplot library [19].
Spearman correlations were found for pairs of elements between datasets. Elements are defined as the following: various agronomic factors from the agronomic data, amount of metabolites from the metabolite data, and relative expression from the qPCR data. This determined if any elements correlated or anticorrelated across all genotypes. For each comparison, calculations were performed using only shared genotypes. Comparisons between agronomic factors and amount of metabolites were visualized in a correlation network using Cytoscape [20]. Only comparisons with a Spearman index greater than 0.5 and less than −0.5 were included in the network visualization.

3. Results

3.1. Starch, Dry Matter Content and Agronomic Data

In the cassava landraces analyzed, starch ranged from 50% to 80% dry weight. The genotype, white cassava Kaleso, had the highest starch content and IBA083565, yellow cassava, has the lowest starch content. The starch assay also showed that several genotypes with yellow storage roots have a high percentage of starch (Figure 2A,B). The agronomic field data can be found in Table S2.
Dry matter content (DMC) was measured in grams dry weight and calculated as percentages. We observed that the white variety Kaleso had the highest DMC for both seasons, while IBA130799 had the lowest DMC in 2018 and IBA102103 had the lowest DMC in 2019. Comparing the two seasons, we observed a lower DMC in 2018 for all genotypes except IBA102103. Based on the meteorological data collected, the 2018 season was much drier compared to the 2019 season, which is likely to explain the variance in DMC between seasons. The dry matter content could not be measured for two of the cultivars, IBA141092 and IBA141102, in 2018 and are thus not included in Figure 2B. Our results show that starch and dry matter content are positively correlated (Figure 2).

3.2. Carotenoid Measurement by UHPLC–DAD–ToF–MS and iCheck Carotene

Based on root β-carotenoid content, genotypes were classified into low, intermediate/average, and high carotenoid content lines (Table 1). Analysis of leaf and storage root samples by UHPLC–DAD–ToF–MS identified chlorophyll A and B and nine carotenoid types in leaf tissues (Table 2 and Table 3), and three carotenoid types in roots. Most carotenoids measured were below detection level in roots; conversely, in many cases, high levels were identified in leaves. Carotenoids detected in storage roots include all trans-β-carotene, 9-cis-β-carotene, and 15-cis-β-carotene. Total β-carotenoid content varied, with genotype IBA141092 having the highest content of all-trans β-carotene and its isomers in roots (Figure 3). It is important to note that β-carotenoids were also identified in cassava genotypes with white storage roots. Variation in root carotenoid component ranged from 0.3 to 36.1 µg g−1 dry weight for total β-carotenes and 0.2 to 22.7 µg g−1 dry weight for all-trans β-carotene, the highest component across all genotypes. While only three carotenoids were identified in root samples, nine carotenoids were identified in leaf samples, including β-carotenoids, and at much higher levels than in storage roots (Table 2).
Comparison between starch and β-carotenoids showed some cultivars that are high in β-carotenoids were also high in starch content (Figure 4). Pearson and Spearman correlations between the two traits showed a week negative correlation of −0.3 and −0.2, respectively.
iCheck Carotene by BioAnalyt, a portable spectrophotometer, was used to measure total carotenoids based on absorbance. This method identified the genotypes IBA141092, IBA141102, and IBA130799 as high carotenoid cassava genotypes. Of these three cultivars, IBA141092 had the highest carotenoid content, while Kaleso, a cultivar with a white storage root, had a carotenoid content below detection level.
Both iCheck and the the UHPLC–DAD–Tof–MS analyses identified IBA141092 and Kaleso as having the highest and lowest carotene content, respectively. However, LC–MS analysis gave higher total carotene content compared to the iCheck analysis. It is important to note that the samples analyzed by UHPLC–MS were younger than samples analyzed by iCHECK at harvest, with plants being 8 months and 13 months, respectively. High performance liquid chromatography (HPLC) methods have the advantage to quantify individual carotenoids and are reported to give higher carotenoid values than spectrophotometric methods [9,21].

3.3. Network Analysis and Correlations of Agronomic and Metabolite Data

PCA showed the opposite relation between zeaxanthin and violaxanthin, as can be expected in the xantophyll cycle (Figure S1). Similarly, there was an opposite relation between zeaxanthin and chlorophyll A. There was no clear clustering depending on the root color, based on the carotenoid and chlorophyll content in either leaves or roots.
For a deeper analysis of the data, we designed a correlation network based on Spearman correlations based on all agronomic and metabolite data using Cytoscape [20] (Figure 5). This analysis showed strong positive correlation between root yield and total carotenoid content irrespective of whether it was determined by UHPLC–MS or colorimetric (iCheck) methods, and a positive correlation was observed between root yield and β-carotene composition. Weak positive correlations were observed between storage root mass (SRM), storage root count (SRC), harvest index (HI), and carotenoid content. A strong negative correlation was observed between DMC and carotenoid content in both planting seasons. As could be expected, a strong positive correlation was observed between dry matter content and percent starch in planting season 1 (S2); however, there was a much weaker positive correlation observed with planting season 1 (S1).

3.4. Gene Expression

Gene expression levels were determined in the storage roots of six cultivars: Kaleso, IBA011231, IBA083724, IBA011663, IBA141092, and TMS60444 (Figure 6 and Table S3). These genotypes have low to high carotenoid levels based on the LC–MS data. Because this study was partly aimed at selecting suitable genotypes for biotechnological biofortification, we included the cultivar, TMS60444, which is not a Nigerian landrace. While this genotype is susceptible to cassava mosaic virus infection and therefore not ideal for cultivation, it is easy to handle in tissue culture and relatively easy to transform, and can thus be used for proof-of-concept for various genetic modifications.
The following five genes encoding enzymes in the carotenoid biosynthesis and cleavage pathways were analyzed by qPCR: phytoene synthases 1 and 2 (PSY1 and PSY2), which catalyzes the condensation of two molecules of geranylgeranyl pyrophosphate to form phytoene; β-carotenoid hydroxylase (CHYβ), which catalyzes the hydroxylation of β-carotene to β-cryptoxanthin and zeaxanthin; lycopene-epsilon cyclase (LCYε), which catalyzes the cyclization of lycopene to α-carotene; and lastly, 9-cis-epoxycarotenoid dioxygenase (NCED1), which catalyzes the formation of abscisic acid.
The genotype IBA011663 (D), had the highest expression of PSY2, CHYβ and LCYε, while IBA141092 (E) showed the highest expression of PSY1 and relatively high expression of PSY2 compared to other cultivars. Kaleso (A), a white root cultivar, had the lowest expression for all five genes and genotype IBA011231 (B) had the highest expression level of NCED1. The differences in carotenoid gene expression levels and carotenoid composition between the genotypes show the complexity in the regulation and accumulation of carotenoids.
We compared gene expression, carotenoid composition, and starch in five selected genotypes and observed a negative correlation between PSY2, LCYε, and NCED1, and starch (Figure 6).

4. Discussion

4.1. Agronomic Traits and Metabolite Profiling

In this study, we identify some cultivars with high β-carotene content and high starch and dry matter content (Figure 2 and Figure 3). This proves that some yellow root cultivars can be an adequate source of both starch and provitamin A carotenoids, and partly contradicts the generally reported negative correlation between DMC and total carotenoid content in cassava storage roots [13]. According to Ceballos et al. [9], dry matter content increased along with carotenoid content following years of recurrent selection, suggesting that cassava germplasms are able to combine high total carotenoid (TC) content with adequate DMC as observed in a genotype with 25.8 mg g–1 TC and an average dry matter content of 37.1%.
We observed a positive correlation between root yield and carotenoid content, coupled to a low level of root rot in cultivars with high carotenoid content (Figure 5). This may be as a result of the antioxidant properties of carotenoids, which are thought to delay the onset postharvest physiological deterioration (PPD) in cassava. Sánchez et al. [22] previously reported an inverse association between total carotenoid content in cassava roots and PPD, suggesting that the onset of PPD is reduced or delayed in cassava roots with high carotenoid content. We believe the higher root yield observed at time of harvest in the cassava cultivars with high carotenoid content can also be attributed to the antioxidant properties of carotenoids.
Only three of the eleven identified carotenoids were detected in storage roots: all-trans β-carotene and its cis isomers, 9-cis-, and 15-cis-β-carotene. A similar finding was published in a study by Ceballos et al. [9] in which most carotenoids detected by HPLC–DAD were in the form of all-trans β-carotene, accounting for 61.2% of total β-carotene (TBC). Other isomers detected in the study were 9-cis-, 13-cis-, and 15-cis-β-carotene, which accounted for 82.2% of total carotenoid. Contrary to our findings, the study of Ceballos et al., which also included a pink variety, identified xanthophylls such as violaxanthin, antheraxanthin, lutein, lycopene, and β-cryptoxanthin in several lines of yellow cassava root.
Chlorophylls are leaf photosynthetic pigments and harvesters of light in plants. In the present study, the chlorophyll A/B ratio ranged from 6.7 to 8. Even if this ratio is species-specific, this value is high compared to a previous report from two cassava cultivars [23]. The proportions of chlorophyll A and B varies according to leaf age and physiological stress [24]. For instance, a decrease in chlorophyll A/B ratio in rice seedlings has been linked to leaf senescence.

4.2. Association between Carotenoid Content and Gene Expression

Carotenoid levels and composition are regulated in many different ways [25]. In our study, PSY2, LCYε, and CHYβ expression were generally positively associated to β-carotene content (Figure 6). This is in line with the findings of Carvalho et al., which showed a correlation between LCYε expression and total β-carotene content, as well as correlation between PDS and total xanthophyll content [15]. A high expression of PSY is thought to drive the flux through the carotenoid synthesis pathway which results in the formation of lycopene [4]. LCYβ and LCYε both compete for lycopene [6]. Regulation of carotenoid accumulation is seen at the molecular level in the expression of the LCYε and CHYβ which serve to drive carotenoid synthesis further down the pathway. This leads to the formation of either lutein or zeaxanthin. Subsequent epoxidation and de-epoxidation of zeaxanthin by zeaxanthin epoxidase (ZEP) and violaxanthin de-epoxidase (VDE) forms the xanthophyll cycle. This may be the explanation for the high expression of CHYβ in the cultivar IBA011663. It could therefore be inferred for IBA011663 that both α- and β-branches of the carotenoid synthesis pathway were favored equally, and that higher expression of one gene produces a bias for its associated branch.
We see a positive correlation between CHYβ expression and total carotenoid content in cassava roots. Even if the number of lines tested for expression in this study is limited, this is in contrast to many transgenic studies. It has been shown a number of times, also in other species, that high CHYβ expression is linked to low β-carotene content. For example, RNAi-silenced CHYβ expression in white sweet potato increased the overall carotenoid contentt and also increased salt stress tolerance. In fact, the suppression of CHYβ expression in white sweet potato led to a doubling of β-carotene in storage roots [25]. Silenced CHYβ1 and CHYβ2 in potato also achieved a 38-fold increase in β-carotene [26]. However, during testing of South American landraces, a white cultivar with low β-carotene content was associated with high CHYβ expression [15]. These findings suggests that the level of expression of CHYβ influences the β-carotene content. In a similar manner, LCYε expression levels were also shown to be negatively correlated with β-carotene content. For example, in potato, an RNAi construct downregulating LCYε led to an increase in β-carotene levels [26].
Various other factors could affect the synthesis and accumulation of carotenoids in cassava. Plastids are the sites for carotenoid biosynthesis in plants, and the type and size of plastids influence carotenoid accumulation and stability [4]. While chloroplasts are the site of carotenoid synthesis, chromoplasts are the carotenoid-accumulating plastids and the structure of the chromoplasts determines the amount of carotenoids accumulated. Crystalline chromoplasts are known to overaccumulate lycopene and β-carotene [5]. Chromoplasts are thought to exert a great influence on carotenoid accumulation. For instance, in cauliflower, the transcript levels of carotenoid synthesis pathway genes are similar in both the white and orange curd of the Or mutant and no significant differences in the expression of carotenoid pathway genes were seen in orange fleshed and nonorange fleshed melon fruits [5]. In addition to this ongoing study, further investigation is needed to understand the underlying mechanism of chromoplast biogenesis in cassava.

5. Conclusions

Biofortification of staple crops is one possible way to alleviate vitamin A deficiency in sub-Saharan Africa. Several studies observed a negative correlation between carotenoid content and dry matter content in African cassava germplasm. However, this present study identified cassava genotypes with high carotenoid content and adequate amount of dry matter and starch content. Expression of carotenoid genes indicated a positive correlation between PSY2, LCYε, and CHYβ expression and total carotenoid content of cassava landraces. A major observation was that gene expression and regulation of carotenoid synthesis may differ between cassava genotypes and there is a high possibility that more than one regulatory mechanism is involved in accumulation and maintenance of carotenoid accumulation. This study identified a suitable cassava genotype, IBA011663, for in vitro transformation. Since this cultivar has an adequate amount of dry matter and starch content, we would like to see if it is able to retain these properties if the carotenoid content is increased by genetic engineering.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4395/10/3/424/s1, Table S1: qPCR primer sequences, Table S2: Agronomic data for average dry matter (DM) and average total carotenoid (TC) composition in storage roots at 13 months, Table S3: Relative expression determined by qPCR of carotenoid genes in six cassava genotypes, Figure S1: PCA loadings of leaf (A) and root (B) metabolites.

Author Contributions

P.O., L.S. (Livia Stavolone), and E.A. planned the experiments and wrote the first draft of the manuscript; HPLC and starch analyses were done by L.S. (Linnea Stridh); UHPLC–MS analysis was done by S.B.; P.O. conducted field sampling and performed qPCR; I.R. provided the agronomic field data; A.L. did the statistical analysis and correlation plots; P.O., L.S. (Livia Stavolone), and E.A. did the main interpretation of the results. All authors contributed to the writing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Vetenskapsrådet: 2016-05703 and Vetenskapsrådet: 2018-04660.

Acknowledgments

This project was supported by Swedish Research Council grants 2016-05703 and 2018-04660. The authors are grateful to Marit Lenman for providing guidance for qPCR data analysis; Elizabeth Parkes for providing high carotenoid cassava genotypes; Peter Illuebey for providing cassava stem cuttings, field labor, and management throughout the experiment; Anetor Omonuwa and Ayo Akinwumi for providing field support for planting, sampling, and harvesting. A.L. would like to acknowledge that this work was completed during her fellowship with the Fulbright student program.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. El-sharkawy, M.A. Drought-Tolerant Cassava for Africa, Asia, and Latin America. BioScience 1993, 7, 441–451. [Google Scholar] [CrossRef]
  2. Esuma, W.; Kawuki, R.S.; Herselman, L.; Labuschagne, M.T. Stability and genotype by environment interaction of provitamin A carotenoid and dry matter content in cassava in Uganda. Breed. Sci. 2016. [Google Scholar] [CrossRef] [Green Version]
  3. Welsch, R.; Arango, J.; Bar, C.; Salazar, B.; Al-Babili, S.; Beltran, J.; Chavarriaga, P.; Ceballos, H.; Tohme, J.; Beyer, P. Provitamin A Accumulation in Cassava (Manihot esculenta) Roots Driven by a Single Nucleotide Polymorphism in a Phytoene Synthase Gene. Plant Cell 2010, 22, 3348–3356. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Nisar, N.; Li, L.; Lu, S.; Khin, N.C.; Pogson, B.J. Carotenoid metabolism in plants. Mol. Plant 2015, 8, 68–82. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Sun, T.; Yuan, H.; Cao, H.; Yazdani, M.; Tadmor, Y.; Li, L. Carotenoid Metabolism in Plants: The Role of Plastids. Mol Plant 2018, 11, 58–74. [Google Scholar] [CrossRef] [Green Version]
  6. Giuliano, G. Provitamin A biofortification of crop plants: A gold rush with many miners. Curr. Opin. Biotechnol. 2017, 44, 169–180. [Google Scholar] [CrossRef]
  7. Rao, C.K.; Annadana, S. Nutrient Biofortification of Staple Food Crops: Technologies, Products and Prospects. In Phytonutritional Improvement of Crops; Wiley Blackwell: Hoboken, NJ, USA, 2017; pp. 113–183. [Google Scholar] [CrossRef]
  8. Bohn, T. Determinants and Determination of Carotenoid Bioavailability from Infant Food Formulas and Adult Nutritionals Including Liquid Dairy Products. J. AOAC Int. 2019, 102, 1044–1058. [Google Scholar] [CrossRef]
  9. Ceballos, H.; Morante, N.; Sanchez, T.; Ortiz, D.; Aragon, I.; Chavez, A.L.; Pizarro, M.; Calle, F.; Dufour, D. Rapid Cycling Recurrent Selection for Increased Carotenoids Content in Cassava Roots. Crop Sci. 2013, 53, 2342–2351. [Google Scholar] [CrossRef] [Green Version]
  10. Low, J.W.; Mwanga, R.O.M.; Andrade, M.; Carey, E.; Ball, A.M. Tackling vitamin A deficiency with biofortified sweetpotato in sub-Saharan Africa. Glob. Food Sec. 2017, 14, 23–30. [Google Scholar] [CrossRef]
  11. Sayre, R.; Beeching, J.R.; Cahoon, E.B.; Egesi, C.; Fauquet, C.; Fellman, J.; Fregene, M.; Gruissem, W.; Mallowa, S.; Manary, M.; et al. The BioCassava plus program: Biofortification of cassava for sub-Saharan Africa. Annu. Rev. Plant Biol. 2011, 62, 251–272. [Google Scholar] [CrossRef]
  12. Atwijukire, E.; Hawumba, J.F.; Baguma, Y.; Wembabazi, E.; Esuma, W.; Kawuki, R.S.; Nuwamanya, E. Starch quality traits of improved provitamin A cassava (Manihot esculenta Crantz). Heliyon 2019, 5, e01215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Njoku, D.N.; Gracen, V.E.; Offei, S.K.; Asante, I.K.; Egesi, C.N.; Kulakow, P.; Ceballos, H. Parent-Offspring regression analysis for total carotenoids and some agronomic traits in cassava. Euphytica 2015, 206, 657–666. [Google Scholar] [CrossRef]
  14. Rabbi, I.Y.; Udoh, L.I.; Wolfe, M.; Parkes, E.Y.; Gedil, M.A.; Dixon, A.; Ramu, P.; Jannink, J.L.; Kulakow, P. Genome-Wide Association Mapping of Correlated Traits in Cassava: Dry Matter and Total Carotenoid Content. Plant Genome 2017, 10. [Google Scholar] [CrossRef] [PubMed]
  15. Carvalho, L.J.C.B.; Agustini, M.A.V.; Anderson, J.V.; Vieira, E.A.; de Souza, C.R.B.; Chen, S.B.; Schaal, B.A.; Silva, J.P. Natural variation in expression of genes associated with carotenoid biosynthesis and accumulation in cassava (Manihot esculenta Crantz) storage root. BMC Plant Biol. 2016, 16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Frede, K.; Schreiner, M.; Baldermann, S. Light quality-Induced changes of carotenoid composition in pak choi Brassica rapa ssp. chinensis. J. Photochem. Photobiol. B 2019, 193, 18–30. [Google Scholar] [CrossRef] [PubMed]
  17. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-Time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  18. De Mendiburu, F. Agricolae: Statistical procedures for agricultural research. R Package Version 2014, 1. [Google Scholar]
  19. Vu, V.Q. ggbiplot: A ggplot2 based biplot. R Package 2011, 342. [Google Scholar]
  20. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  21. Jaramillo, A.M.; Londono, L.F.; Orozco, J.C.; Patino, G.; Belalcazar, J.; Davrieux, F.; Talsma, E.F. A comparison study of five different methods to measure carotenoids in biofortified yellow cassava (Manihot esculenta). PLoS ONE 2018, 13, e0209702. [Google Scholar] [CrossRef]
  22. Sánchez, T.; Chávez, A.L.; Ceballos, H.; Rodriguez-Amaya, D.B.; Nestel, P.; Ishitami, M. Reduction or delay of post-Harvest physiological deterioration in cassava roots with higher carotenoid content. J. Sci. Food sAgric. 2005. [Google Scholar] [CrossRef]
  23. Anand, M.H.; Byju, G. Chlorophyll meter and leaf colour chart to estimate chlorophyll content, leaf colour, and yield of cassava. Photosynthetica 2008, 46, 511–516. [Google Scholar] [CrossRef]
  24. Zhang, Y.; Huang, J.; Wang, F.; Blackburn, G.A.; Zhang, H.K.; Wang, X.; Wei, C.; Zhang, K.; Wei, C. An extended PROSPECT: Advance in the leaf optical properties model separating total chlorophylls into chlorophyll a and b. Sci. Rep. 2017, 7, 6429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kang, L.; Ji, C.Y.; Kim, S.H.; Ke, Q.; Park, S.C.; Kim, H.S.; Lee, H.U.; Lee, J.S.; Park, W.S.; Ahn, M.J.; et al. Suppression of the beta-Carotene hydroxylase gene increases beta-Carotene content and tolerance to abiotic stress in transgenic sweetpotato plants. Plant Physiol. Biochem. 2017, 117, 24–33. [Google Scholar] [CrossRef] [PubMed]
  26. Diretto, G.; Welsch, R.; Tavazza, R.; Mourgues, F.; Pizzichini, D.; Beyer, P.; Giuliano, G. Silencing of beta-Carotene hydroxylase increases total carotenoid and beta-Carotene levels in potato tubers. BMC Plant Biol. 2007, 7, 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Schematic carotenoid biosynthesis pathway showing metabolites and genes involved in their synthesis. Genes measured by quantitative real-time polymerase reaction (qPCR) are marked. Geranylgenranyl diphosphosphate synthase (GGPPS), Phytoene synthase (PSY), phytoene desaturase (PDS), ζ-carotene desaturase (ZDS), ζ-carotene isomerase (ZISO), carotenoid isomerase (CRTISO), lycopene ε-cyclase (LCYε), lycopene β-cyclase (LCYβ) β-carotenoid hydroxylases (CHYβ), zeaxanthin epoxidase (ZEP), violaxanthin de-epoxidase (VDE), 9-cis-expoxycarotenoid dioxygenase.
Figure 1. Schematic carotenoid biosynthesis pathway showing metabolites and genes involved in their synthesis. Genes measured by quantitative real-time polymerase reaction (qPCR) are marked. Geranylgenranyl diphosphosphate synthase (GGPPS), Phytoene synthase (PSY), phytoene desaturase (PDS), ζ-carotene desaturase (ZDS), ζ-carotene isomerase (ZISO), carotenoid isomerase (CRTISO), lycopene ε-cyclase (LCYε), lycopene β-cyclase (LCYβ) β-carotenoid hydroxylases (CHYβ), zeaxanthin epoxidase (ZEP), violaxanthin de-epoxidase (VDE), 9-cis-expoxycarotenoid dioxygenase.
Agronomy 10 00424 g001
Figure 2. (A) Starch in dried 8-month-old root tissue from 13 cassava genotypes using starch enzymatic assay. (B) Dry matter content (DMC, %). Error bars represent the standard deviation of the mean from three biological replicates with three technical repetitions per genotype. Groupings generated from ANOVA and Tukey ad hoc tests are above each bar.
Figure 2. (A) Starch in dried 8-month-old root tissue from 13 cassava genotypes using starch enzymatic assay. (B) Dry matter content (DMC, %). Error bars represent the standard deviation of the mean from three biological replicates with three technical repetitions per genotype. Groupings generated from ANOVA and Tukey ad hoc tests are above each bar.
Agronomy 10 00424 g002
Figure 3. (A) All-trans β-carotene concentration in µg g–1 dry weight of cassava. (B) Total β-carotene in µg g–1 dry weight of cassava. Error bars represent the standard deviation from the mean. Groupings generated from ANOVA and Tukey adhoc tests are above each bar.
Figure 3. (A) All-trans β-carotene concentration in µg g–1 dry weight of cassava. (B) Total β-carotene in µg g–1 dry weight of cassava. Error bars represent the standard deviation from the mean. Groupings generated from ANOVA and Tukey adhoc tests are above each bar.
Agronomy 10 00424 g003
Figure 4. Comparison of total β-carotene and starch content for cassava genotypes.
Figure 4. Comparison of total β-carotene and starch content for cassava genotypes.
Agronomy 10 00424 g004
Figure 5. Correlations between agronomic traits and carotenoids. Red circular nodes are agronomic traits, while purple diamond nodes are metabolites. Green lines represent positive correlations, while red lines represent negative correlations. The thickness of the lines reflects the strength of the correlation. Spearman correlation cutoff are values greater than 0.5 or less than −0.5 (Rt = roots, Lf = leaf, S1 = season 1, S2 = season 2).
Figure 5. Correlations between agronomic traits and carotenoids. Red circular nodes are agronomic traits, while purple diamond nodes are metabolites. Green lines represent positive correlations, while red lines represent negative correlations. The thickness of the lines reflects the strength of the correlation. Spearman correlation cutoff are values greater than 0.5 or less than −0.5 (Rt = roots, Lf = leaf, S1 = season 1, S2 = season 2).
Agronomy 10 00424 g005
Figure 6. Relative gene expression (RQ) and comparison between selected metabolites across five cassava genotypes. (A)–(E) show relative expression of PSY1, PSY2, LCYε, CHYβ and NCED1 and their expression correlation with starch and total β-carotene composition across selected genotypes.
Figure 6. Relative gene expression (RQ) and comparison between selected metabolites across five cassava genotypes. (A)–(E) show relative expression of PSY1, PSY2, LCYε, CHYβ and NCED1 and their expression correlation with starch and total β-carotene composition across selected genotypes.
Agronomy 10 00424 g006
Table 1. Genotypes classification as low, intermediate, and high carotenoid content in storage roots based on the UHPLC–DAD–ToF–MS analysis. Three plants per genotype were analyzed. W: white, PY: pale yellow, Y: yellow.
Table 1. Genotypes classification as low, intermediate, and high carotenoid content in storage roots based on the UHPLC–DAD–ToF–MS analysis. Three plants per genotype were analyzed. W: white, PY: pale yellow, Y: yellow.
GenotypeFlesh ColorTotal ß-CaroteneCategory
KALESOW0.3low
IBA050128PY6.2intermediate
IBA011231PY4.4intermediate
IBA083724Y11.2high
IBA090144W14.2high
Z010116W0.6low
IBA083565Y12.4high
MM990214W2.7low
IBA011663Y14.9high
IBA102103W15high
IBA141092Y36.1high
IBA141102Y14.7high
IBA130799Y12.4high
Table 2. Carotenoid composition of leaf carotenoids (µg g−1 leaf dry weight). Three plants per genotype were analyzed. Variation between samples is given as SD.
Table 2. Carotenoid composition of leaf carotenoids (µg g−1 leaf dry weight). Three plants per genotype were analyzed. Variation between samples is given as SD.
GenotypeLuteinβ-Carotene9-cis-β-Carotene15-cis-β-CaroteneZeaxanthinNeoxanthin1Neoxanthin2Violaxanthin1Violaxanthin2
KALESO393 ± 135229 ± 8835 ± 165 ± 2328 ± 2770 ± 2522 ± 755 ± 2446 ± 17
IBA050128529 ± 41326 ± 25648 ± 39 ± 0.8481 ± 9093 ± 1227 ± 5139 ± 4182 ± 20
IBA011231526 ± 71299 ± 7542 ± 107 ± 21849 ± 306106 ± 1617 ± 465 ± 2934 ± 11
IBA083724528 ± 38289 ± 3240 ± 57 ± 1641 ± 8287 ± 921 ± 475 ± 2256 ± 7
IBA090144570 ± 120268 ± 6233 ± 77 ± 2702 ± 3788 ± 1922 ± 662 ± 3333 ± 8
Z010116582 ± 96305 ± 5045 ± 98 ± 2523 ± 167113 ± 2226 ± 7102 ± 4161 ± 13
IBA083565721 ± 29395 ± 3854 ± 910 ± 0.9791 ± 131116 ± 529 ± 5123 ± 3667 ± 22
MM990214781 ± 188402 ± 8451 ± 712 ± 3915 ± 462137 ± 3120 ± 471 ± 452 ± 7
IBA011663575 ± 88329 ± 6144 ± 109 ± 1993 ± 28297 ± 1419 ± 655 ± 1340 ± 8
IBA102103841 ± 36500 ± 3168 ± 215 ± 0.91478 ± 211166 ± 1027 ± 6140 ± 3694 ± 12
IBA141092812 ± 74405 ± 3349 ± 211 ± 1734 ± 34125 ± 1420 ± 286 ± 1261 ± 10
IBA141102722 ± 48366 ± 3145 ± 710 ± 11201 ± 135118 ± 127 ± 587 ± 2843 ± 7
IBA130799716 ± 26385 ± 8748 ± 1411 ± 41165 ± 315119 ± 2223 ± 193 ± 550 ± 15
Table 3. Chlorophyll A and B composition in leaf samples (µg g−1). Three plants per genotype were analyzed.
Table 3. Chlorophyll A and B composition in leaf samples (µg g−1). Three plants per genotype were analyzed.
GenotypeChlorophyll AChlorophyll BChl A: Chl B
KALESO8006 ± 2688913 ± 3148.77
IBA0501289775 ± 3971111 ± 1068.8
IBA0112319025 ± 11851173 ± 1697.69
IBA0837248453 ± 9991096 ± 1097.71
IBA0901447856 ± 16891060 ± 2657.41
Z01011610,229 ± 14151321 ± 2217.75
IBA08356510,564 ± 8021397 ± 977.56
MM99021411,653 ± 20371609 ± 4037.24
IBA0116639087 ± 18611220 ± 2677.45
IBA10210313,231 ± 9471841 ± 757.19
IBA14109210,215 ± 8121380 ± 847.4
IBA14110210,281 ± 8111525 ± 1246.74
IBA13079910,238 ± 21261406 ± 3687.28

Share and Cite

MDPI and ACS Style

Olayide, P.; Large, A.; Stridh, L.; Rabbi, I.; Baldermann, S.; Stavolone, L.; Alexandersson, E. Gene Expression and Metabolite Profiling of Thirteen Nigerian Cassava Landraces to Elucidate Starch and Carotenoid Composition. Agronomy 2020, 10, 424. https://doi.org/10.3390/agronomy10030424

AMA Style

Olayide P, Large A, Stridh L, Rabbi I, Baldermann S, Stavolone L, Alexandersson E. Gene Expression and Metabolite Profiling of Thirteen Nigerian Cassava Landraces to Elucidate Starch and Carotenoid Composition. Agronomy. 2020; 10(3):424. https://doi.org/10.3390/agronomy10030424

Chicago/Turabian Style

Olayide, Priscilla, Annabel Large, Linnea Stridh, Ismail Rabbi, Susanne Baldermann, Livia Stavolone, and Erik Alexandersson. 2020. "Gene Expression and Metabolite Profiling of Thirteen Nigerian Cassava Landraces to Elucidate Starch and Carotenoid Composition" Agronomy 10, no. 3: 424. https://doi.org/10.3390/agronomy10030424

APA Style

Olayide, P., Large, A., Stridh, L., Rabbi, I., Baldermann, S., Stavolone, L., & Alexandersson, E. (2020). Gene Expression and Metabolite Profiling of Thirteen Nigerian Cassava Landraces to Elucidate Starch and Carotenoid Composition. Agronomy, 10(3), 424. https://doi.org/10.3390/agronomy10030424

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