ORIGINAL RESEARCH article
Novel Bioengineered Cassava Expressing an Archaeal Starch Degradation System and a Bacterial ADP-Glucose Pyrophosphorylase for Starch Self-Digestibility and Yield Increase
- 1College of Agriculture and Related Sciences, Delaware State University, Dover, DE, United States
- 2Department of Biochemistry and Molecular Biology, University of Georgia, Athens, GA, United States
- 3Department of Microbiology and Molecular Genetics, Michigan State University, East Lansing, MI, United States
- 4Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI, United States
To address national and global low-carbon fuel targets, there is great interest in alternative plant species such as cassava (Manihot esculenta), which are high-yielding, resilient, and are easily converted to fuels using the existing technology. In this study the genes encoding hyperthermophilic archaeal starch-hydrolyzing enzymes, α-amylase and amylopullulanase from Pyrococcus furiosus and glucoamylase from Sulfolobus solfataricus, together with the gene encoding a modified ADP-glucose pyrophosphorylase (glgC) from Escherichia coli, were simultaneously expressed in cassava roots to enhance starch accumulation and its subsequent hydrolysis to sugar. A total of 13 multigene expressing transgenic lines were generated and characterized phenotypically and genotypically. Gene expression analysis using quantitative RT-PCR showed that the microbial genes are expressed in the transgenic roots. Multigene-expressing transgenic lines produced up to 60% more storage root yield than the non-transgenic control, likely due to glgC expression. Total protein extracted from the transgenic roots showed up to 10-fold higher starch-degrading activity in vitro than the protein extracted from the non-transgenic control. Interestingly, transgenic tubers released threefold more glucose than the non-transgenic control when incubated at 85°C for 21-h without exogenous application of thermostable enzymes, suggesting that the archaeal enzymes produced in planta maintain their activity and thermostability.
Global energy consumption is predicted to increase by 48% from 580 quadrillion KJ in 2012 to 860 quadrillion KJ in 2040 (U.S. EIA, 2016). However, given the increasing concerns of environmental pollution and climate change due to production and use of petroleum-based fuels and chemicals (Mei and Rakhmatov, 2014), there is a need to develop alternative and renewable energy sources to meet the demands. Bioethanol as biofuel has been widely produced from sugarcane and corn in Brazil and United States, respectively, (Mei and Rakhmatov, 2014; Lopes et al., 2017). In the 2010/2011 crop season alone, grain ethanol refineries consumed 40% of the total United States corn production (Du and McPhail, 2012). This large diversion of grain to ethanol has led to growing concern over crop-based bioethanol due to its competition with food and feed supplies. Hence, there is a strong interest in second and third generation biofuels (e.g., lignocellulosic biomass and algae) whose technologies are still expensive due to the cost of feedstock production and processing (Carriquiry et al., 2011; Balan, 2014). Therefore, there is a need to identify alternative plant species that can be economically more viable in the near term, are compatible with existing starch-based biofuel technologies, and are high yielding, resilient, and multipurpose (Sanchez and Cardona, 2008; Cooke and Traore, 2010). Starch-rich feedstocks such as cassava (Manihot esculenta Crantz) offer great potential to address these needs. Given the ever-increasing demand for food, feed, and energy, as well as continuous changes in our societal demands, there is a need for cost-effective and sustainable starch-processing routes. Recent advances in genetic engineering have enabled the transfer of starch hydrolytic enzymes from the processing tank to starchy feedstock by directly producing enzymes in the biomass (Hebelstrup et al., 2015). This transfer requires the use of thermostable and thermoactive enzymes isolated from hyperthermophiles, which have very low activity at ambient temperatures, and thus do not compromise normal plant growth and development.
Cassava is one of the most important food crops, feeding more than a billion people in about 105 countries in Africa, Asia, and Latin America. It can be grown by resource-poor farmers due to its outstanding adaptation to and productivity on marginal lands (El-Sharkawy, 2004; Chetty et al., 2013). Increasing demands for food, feed, and fuel throughout the tropics has led to a 24% increase in land committed to cassava cultivation since 1990, surpassing the Solanum genus potato (Taylor et al., 2012). Cassava is an important source of starch in tropical and subtropical countries, ranking second to corn (De Souza et al., 2016; Karlstrom et al., 2016). According to the Food and Agriculture Organization of the United Nations, about 8.5 million tons of cassava starch and flour were exported globally in 2015 (FAO, 2015), of which Thailand contributed the largest share (7.9 million tons). The demand from ethanol sectors in Asia is driving rapid growth in world cassava production by an estimated 100 million tons (∼35% of total production) since 2000 (FAO, 2013), with at least 780 million liters of cassava ethanol produced annually in China alone (Xuan et al., 2015). Cassava is inherently able to produce very high harvestable yields (>070 tons/ha) of starchy roots in optimum environments, compared to corn (20 tons/ha), sorghum (13 tons/ha), rice (26 tons/ha), sweet potato (65 tons/ha), and wheat (12 tons/ha) (El-Sharkawy, 2004). Moreover, its remarkable tolerance to adverse conditions, such as poor soils, drought, and high temperature (El-Sharkawy and Cock, 1984; Angelov et al., 1993; El-Sharkawy, 1993; Fermont et al., 2009), where most other staple food crops would fail to produce reasonable yields (Jarvis et al., 2012), makes cassava one of the most promising crops in the face of climate change.
Improvement of cassava through traditional breeding has been challenging due to high heterozygosity, low natural fertility, erratic flowering and the lack of agronomically important genes in sexually compatible germplasm (Puonti-Kaerlas et al., 1999; Ihemere et al., 2006; Que et al., 2010; Ceballos et al., 2015; Chavarriaga-Aguirre et al., 2016). Therefore, genetic engineering is a preferred approach for trait improvement. Several laboratories have genetically engineered cassava for agronomically important traits such as vitamin content (Welsch et al., 2010; Li et al., 2015), virus resistance (Vanderschuren et al., 2009, 2012; Beyene et al., 2017), extended shelf-life (Zidenga et al., 2012; Xu et al., 2013), reduced toxic cyanogenic glycoside content (Siritunga and Sayre, 2004; Jorgensen et al., 2005), improved starch quality (Zhao et al., 2011), and increased cold-tolerance (Xu et al., 2014). Likewise, enhanced starch content has been reported in transgenic cassava expressing a mutated Escherichia coli ADP-glucose pyrophosphorylase (AGPase; glgCGly336Asp), which is less sensitive to its activator (fructose-1, 6-bisphosphate) and inhibitors (AMP and inorganic phosphate) (Stark et al., 1992; Meyer et al., 1998; Ballicora et al., 2003; Ihemere et al., 2006). AGPase catalyzes the first dedicated and rate-limiting step in starch biosynthesis (Preiss, 1988; Tiessen et al., 2002). Given that cassava is inherently able to produce high yields (El-Sharkawy, 2004), and has the highest annual bioethanol yield compared to other starchy feedstocks such as sweet sorghum, rice, maize, and wheat (Wang, 2002; Ziska et al., 2009), there is a need to develop novel cultivars with increased starch yields and starch self-processing ability for enhanced bioethanol production.
Hyperthermophilic bacteria and archaea growing optimally at temperatures above 80°C and have been a source of thermostable industrial enzymes (Richardson et al., 2002; Egorova and Antranikian, 2005). Complete starch breakdown to glucose requires the synergistic action of thermostable amylases, pullulanases, and glucoamylases (Egorova and Antranikian, 2005). Archaeal amylases (EC 126.96.36.199) have been isolated and characterized from Pyrococcus furiosus, P. woesei, Thermococcus profundus, and T. hydrothermalis (Dong et al., 1997a; Jorgensen et al., 1997; Bertoldo and Antranikian, 2001; Mehta and Satyanarayana, 2013). Similarly, pullulanases (EC 188.8.131.52) have been isolated and characterized from the archaea P. furiosus, P. woesei, and Desulfurococcus mucosus (Rudiger et al., 1995; Dong et al., 1997b; Duffner et al., 2000; Lee et al., 2006; Hii et al., 2012). Glucoamylases (E.C. 184.108.40.206), which convert partially processed starch/dextrin to glucose, have also been isolated from archaea such as Sulfolobus solfataricus, Picrophilus torridus, and Thermoplasma acidophilum (Bertoldo and Antranikian, 2001; Kim et al., 2004; Kumar and Satyanarayana, 2009) and characterized.
Several laboratories have successfully expressed thermostable bacterial α-amylase/amyllopullulanase in plants, including corn (Lanahan et al., 2006), rice (Chiang et al., 2005; Xu et al., 2008) and sweet potato (Santa-Maria et al., 2011) to induce starch self-processing ability. However, it remains unclear as to whether increased starch self-processing ability could be conferred by heterologous expression of a suite of hyperthermophilic archaeal enzymes in planta. Moreover, expressing multigene starch-degrading genes in cassava from a single T-DNA has never been reported before.
In this study, mutated E. coli glgCG336D as well as the genes encoding three archaeal starch-hydrolyzing enzymes, α-amylase and amylopullulanase from P. furiosus and glucoamylase from S. solfataricus, were simultaneously expressed in cassava to enhance yield and starch self-processing ability. Transgenic lines were phenotypically and genotypically characterized. Our findings reveal that simultaneous expression of the archaeal and bacterial genes improves starch self-processing ability and storage root yield, respectively. To our knowledge, this is the first report of in planta production and activation of multiple archaeal enzymes that led to increased starch hydrolysis directly in the biomass. This approach has a vast potential to improve cassava and other starch-rich feedstock for the production of bioethanol and industrial products.
Materials and Methods
The modified P. furiosus α-amylase gene containing the catalytic and Ca-binding domains (amylaseΔ1-37 of accession number WP_014835153, designated herein as α-amylase), and amylopullulanase gene (accession number NC_018092, Bridger et al., 2012), and the S. solfataricus glucoamylase gene (accession number NC_017274) were amplified from the chromosomal DNA of each archaea using Phusion Hot Start II High-Fidelity DNA Polymerase (Fisher Scientific, Pittsburgh, PA, United States). The amylase sequence starts with the first methionine and lacks 37 residues from the amino-terminus as compared to the sequence reported earlier (Dong et al., 1997a; Jorgensen et al., 1997). The amylase sequence was expressed in E. coli strain Rosetta2-DE3-pLyS and the activity of recombinant protein in degrading starch was validated (Supplementary Methods and Supplementary Figure S1A). The amylopullulanase (NC_018092) used in this study encodes the entire coding region (1355 amino acids). Prior to using this sequence in the multigene assembly its activity, was also validated using a recombinant protein produced in E. coli and pullulan as a substrate (Supplementary Methods and Supplementary Figure S1B). The coding region of E. coli glgC (Acc#S58224) was amplified from pO12 plasmid obtained from Dr. Tony Romeo (University of Florida). A 57-amino acid pea chloroplast transit peptide was fused to the N-terminus of glgC to target the protein to the amyloplast, a non-pigmented organelle responsible for starch synthesis and storage. Since single amino acid substitution Gly336Asp has been shown to reduce sensitivity of the enzyme to inhibitors and activators (Stark et al., 1992; Meyer et al., 1998; Ihemere et al., 2006), the mutation was introduced by site-directed mutagenesis. The root-specific promoters of the potato class I Patatin and granule-bound starch synthase (GBSS) were amplified from pUC19 and TA vectors, respectively, obtained from Dr. William Belknap (USDA-ARS, Albany, CA, United States). All primers used to amplify promoters or genes used in this study are listed in Supplementary Table S1. The sequence of each primer includes 5′ extra nucleotides, a restriction enzyme site, and a sequence specific to the target promoter or gene.
Generation of Expression Cassettes and Multigene Assembly
Expression cassettes were generated in pSAT modular vectors (pSAT1, pSAT4, pSAT5, and pSAT6) (Chung et al., 2005; Tzfira et al., 2005) under the control of Patatin or GBSS promoters. The Patatin promoter was cloned into pSAT1 and pSAT5, while GBSS was cloned into pSAT4 and pSAT6. Both promoters were cloned upstream of the multicloning sites by replacing the enhanced CaMV promoter (2x35s) in pSAT1, pSAT4, and pSAT6 and the mannopine synthase promoter (mas-P) in pSAT5 (Figure 1A). Expression cassettes Patatin-Amy-35sT (2.57 kb) for α-amylase, GBSS-Pull-35sT (5.68 kb) for amylopullulanase, Patatin-glgC-masT (2.72 kb) for glgC and GBSS-Gluco-35sT (3.48 kb) for glucoamylase were generated in pSAT1 (Figure 1A-i), pSAT4 (Figure 1A-ii), pSAT5 (Figure 1A-iii), and pSAT6 (Figure 1A-iv), respectively. The cassettes were then excised from the pSAT vector using AscI and homing endonucleases I-SceI, I-CeuI, and PI-PspI, respectively, and assembled into pPZP-RCS2-ocs-nptII binary vector (8.48 kb) that contains spectinomycin and neomycin phosphotransferase (NPTII, kanamycin) resistance genes for selection in bacterial and plant cells, respectively (Goderis et al., 2002). Intactness of the expression cassettes in the resulting ∼25.6 kb binary vector (Figure 1B) was validated using AscI and the homing endonucleases (Figure 1C) prior to introduction into Agrobacterium strain LBA4404 competent cells. Agrobacterium harboring the binary vectors was used to transform cassava friable embryogenic calli (FEC).
FIGURE 1. Multigenetic constructs used for cassava transformation. (A) Expression cassettes of Pyrococcus furiosus α-amylase (pSAT1-Amy, A-i), P. furiosus amylopullulanase (pSAT4-Pull, A-ii), Escherichia coli glgC (pSAT5-glgC, A-iii), and Sulfolobus solfataricus glucoamylase (pSAT6-Gluco, A-iv) constructed in pSAT modular vectors under the control of Patatin or GBSS promoters. (B) The four expression cassettes from (A) are assembled into pPZP-RCS2 containing the neomycin phosphotransferase selectable marker (NPTII), yielding ∼a 25.61 kb plasmid. (C) Validation of the expression cassettes with rare-cutting restriction enzymes (AscI, I-SceI, I-CeuI, and PI-PspI). Lower bands in lanes 1, 2, 3, and 4 represent expression cassettes of Amy, Pull, Gluco, and glgC, respectively, as indicated by the arrows. Lane 5 is undigested vector containing multigene expression cassettes (B).
Production of Friable Embryogenic Calli (FEC)
Cassava cultivar TMS60444 obtained from Danforth Plant Science Center was maintained in controlled growth chamber at 27°C, 80% relative humidity and 16-h/8-h day/night cycles. Stem segments with at least one shoot bud were cultured on MS medium (Murashige and Skoog, 1962) with supplements [CBM, 4.43 g/L MS salt with vitamins, 2% (w/v) sucrose, 2 mM CuSO4, Supplementary Table S2]. Friable embryogenic calli were generated according to Taylor et al. (2001, 2012) The FEC were continuously subcultured on GD (Greshoff and Doy, 1972) 2-50P (Supplementary Table S2) media until homogenous clusters were obtained. Homogenous calli were maintained on GD6-50P (supplemented with 6% sucrose) before they were used for transformation or before the calli were used to initiate embryogenic suspension (ES). Embryogenic suspension of the organized calli was initiated and maintained in Schenk and Hildebrandt SH6 (Schenk and Hildebrandt, 1972) modified medium [3.2 g/L SH salt, 6% (w/v) sucrose, 42 μM picloram, Supplementary Table S2]. Highly proliferating FEC or ES were used for transformation.
Agrobacterium-Mediated Transformation of Target Tissue
Agrobacterium tumefaciens strain LBA4404 harboring the control binary vector pPZP-RCS2-ocs-nptII and the vectors containing expression cassette of each starch-hydrolyzing enzymes (α-amylase, amylopullulanase, and glucoamylase genes) with or without glgC were used for transformation. Single colonies of LBA4404 carrying each construct were cultured in 2 mL LB medium containing appropriate antibiotics for about 10-h while shaking at 220 rpm at 28°C. About 1–1.5 mL of the suspension was used to initiate 25 mL of YM medium containing appropriate antibiotics and 100 μM acetosyringone. The culture was allowed to grow overnight to reach an OD600 of 0.5–0.75. The bacteria suspension was prepared according to Taylor et al. (2012). The non-ionic surfactant Pluronic F-68 was added to the bacterial suspension at a concentration of 0.03%. Embryogenic tissues obtained from GD2-50P agar media (Supplementary Table S1) or suspension culture was inoculated with the bacterial suspension in 6-well plates (Figure 2). At least two independent transformations were done for each construct. The tissues were vacuum-infiltrated for a total of 15 min with two pauses every 5 min, and incubated at room temperature for 30 min. Removal of excess Agrobacterium and coculturing was performed according to Taylor et al. (2012). After 5 days of-coculture, the tissues were transferred to Stage I selection medium GD2-50P containing 25 μM paromomycin and 500 mg/L carbenicillin. Putative transgenic tissues were further selected on GD2-50P containing 30 μM (Stage II) and 40 μM (Stage III) paromomycin each for 3 weeks. Transgene insertion was validated by PCR from genomic DNA samples isolated from proliferating putative transgenic calli.
FIGURE 2. Transformation of cassava friable embryogenic calli (FEC) and generation of transgenic lines. (A) FECs were generated as described in the Section “Materials and Methods.” (B) Inoculation of FEC with Agrobacterium strain LBA4404 containing expression cassettes of starch hydrolyzing enzymes and glgC as in Figure 1. (C) Putative transgenic calli selected under increasing concentrations of paromomycin. (D) Shoot regeneration from transgenic lines. (E) Rooting of transgenic lines. (F) One-month-old non-transformed control and representative transgenic cassava lines. (G) Integration of multiple genes into cassava genome as confirmed by PCR using genomic DNA as template. Primers specific to NPTII, Amy, Pull, Gluco, and glgC were used. When applicable, the bands in lanes 1, 2, 3, 4, and 5 represent NPTII, Amy, Pull, Gluco, and glgC, respectively.
Maturation and Regeneration of Putative Tissues
After selection under increasing paromomycin concentration (25, 30, and 40 μM), resistant calli were transferred to embryo maturation media (MSN, 4.43 g/L MS salt with vitamins, 2% (w/v) sucrose, 1 mg/L NAA, pH5.8) (Bull et al., 2009) as in Supplementary Table S1. Maturing cotyledonary stage embryos were transferred to embryo germination media [CEM, 4.43 g/L MS salt with vitamins, 2% (w/v) sucrose, 0.4 mg/L BAP %, Supplementary Table S2]. Germinating shoots were transferred to CBM without antibiotics for rooting.
Transplanting in Vitro Propagated Plants to Soil
Eighteen-day-old uniform micropropagated plantlets (at least nine per construct) were transplanted in Kord 3.0 square pots filled with 230 mL Fafard 51 mix1 that has been shown to promote tuber development, and treated with Gnatrol to control fungal gnats (Zidenga et al., 2012). The plants were established in the greenhouse under natural light and fertilized every 3 weeks with NPK (9-45-15) fertilizer. After 4 months when initiation of mini tubers was observed, plants were transferred to 2-L pots and grown for another 3 months before harvesting.
Quantitative Real-Time RT-PCR
To determine transcript levels of the archaeal genes in cassava tubers, quantitative RT-PCR was performed using an ABI 7500 real-time PCR system and SYBR Green Kit (Applied Biosystems, Grand Island, NY, United States). Total RNA isolated from the flesh (most inner tissue) of tubers using Spectrum Plant Total RNA Kit (Sigma) was treated with DNase I (Invitrogen, Carlsbad, CA, United States) to remove contaminating genomic DNA prior to first-strand complementary DNA (cDNA) synthesis. First-strand cDNA was synthesized using Superscript III (Life Technologies, Grand Island, NY, United States). One microgram of the first-strand cDNA was used for each real-time RT-PCR reaction containing 2xPower SYBR Green Master Mix and 0.15 μM primers in a final volume of 25 μl. Gene-specific sense and antisense primers (Supplementary Table S3) were used for the amplification, and cassava α-tubulin (TC30055) was used as an internal control. Relative expression levels were calculated using the ΔΔCT method available on SDS software (Applied Biosystems).
Total Protein Extraction from Fibrous Roots and Tubers
Total protein from fibrous or tuberous roots was extracted according to Santa-Maria et al. (2011) with minor modifications. The roots or tubers obtained from four independent plants were ground under liquid N2, and the ground tissue was homogenized for 10 min using 1.5 mL/g of cold extraction buffer containing 50 mM Tris-HCl (pH 8.0), 0.1 M KCl, 5% glycerol, 1.5% polyvinylpolypyrrolidone, 5% dimethylsulfoxide, 5 mM dithiothreitol, and 1x protease inhibitor cocktail (Sigma–Aldrich, St. Louis, MO, United States). The homogenate was centrifuged at 1,500 × g for 10 min. The supernatant was filtered through 250 micron mesh, heated at 80°C for 25 min to eliminate mesophilic enzymes which may interfere with the thermophilic enzyme activity assay, and centrifuged at 4°C for 15 min at 19000 × g. For assaying AGPase activity, total protein was extracted from fibrous roots of soil grown plants. Roots (100 mg) were transferred to 2 mL microcentrifuge tubes containing 1 mL extraction buffer containing 100 mM HEPES (pH 7.5), 8 mM MgCl2, 2 mM EDTA, 1 mM DTT, 12.5% (v/v) glycerol and 5% polyvinylpyrrolidone (Kulichikhin et al., 2016). Tissue was homogenized for two 15-s slow-speed periods followed by two 15-s fast-speed blending periods according to Ihemere et al. (2006). The suspension was filtered through miracloth and the extract was centrifuged for 20 min at 21, 000 × g. The supernatant was recovered and immediately frozen until it was used for AGPase assay. Total protein was estimated using Coomassie (Bradford) protein assay (Thermo Scientific Pierce, Waltham, MA, United States) using bovine serum albumin as the standard.
Assay of Enzymatic Activity
Archaeal enzymatic activity was assayed in 1.5 mL of 50 mM phosphate buffer (pH 6.0) containing 15 μg of heat-purified extract and 15 mg pure potato starch (Sigma). For total protein extracted from mini tubers, the activity of enzymes α-amylase, amylopullulanase, and glucoamylase was assayed in reactions containing 1% potato starch and pullulan (Sigma), and 2% maltose (Sigma)/maltotriose (Fisher Scientific), respectively. The reaction was conducted in phosphate buffer at pH 7.0 for starch and pullulan, and pH 5.6 for maltose/maltotriose. The reactions were incubated at 85°C for 30 min and centrifuged at 10,000 × g for 2 min. Total glucose in 500 μL of the hydrolysate was determined by the modified 3, 5-dinitrosalicylic acid (DNS) method (Miller, 1959). Increase in absorbance due to the release of glucose was measured at OD540. AGPase activity was assayed according to Kulichikhin et al. (2016) based on coupling enzymatic reactions where the product of the initial reaction is used as a substrate for subsequent reactions generating NADPH as the end product. The increase in absorbance due to NADPH was measured at OD340 using a control reaction without the enzyme as a reference.
Assay of Tuber Self-Processing Ability
Saccharification ability of the transgenic tubers was analyzed alongside with the wild-type control. Tubers were harvested and stored at -80°C until used. Frozen tubers were thawed and the bark and peels were removed. The inner starchy tissue (flesh) from same plant was pulled, and ground under liquid N2. Tubers obtained from four individual plants were analyzed for each treatment. To determine free sugar content, 1 g of the ground tissue was mixed thoroughly with 2 mL of water. Samples were centrifuged at 10,000 × g for 5 min and the supernatant was recovered to determine sugar content. To determine tuber starch self-processing ability, samples were incubated at 85°C, a temperature typically used for starch hydrolysis (Nahampun et al., 2013), for 21-h to activate the thermophilic enzymes. The reaction was cooled to room temperature and the samples were centrifuged for 2 min at 10,000 × g. The supernatant was recovered and used to quantify total reducing sugar using the DNS method as mentioned above.
Experiments were conducted in a complete randomized design in at least three replicates, and experiments were repeated twice. Data were analyzed using one-way ANOVA using the PROG GLM procedure (Westfall et al., 1996). After significant F-tests, the Tukey multiple comparison procedure was used to separate the means (P < 0.05).
Enzyme Cocktails from Fibrous Roots of Transgenic Cassava Hydrolyze Starch in Vitro
Starch conversion to glucose requires high temperature treatment and the action of thermostable glycosyl hydrolases that are commercially produced from thermophilic microorganisms. To study whether the archaeal enzymes are produced in the roots, the genes were targeted to the roots using root-specific promoters. A total of 11 lines expressing the three starch-degrading enzymes and the mutated bacterial AGPase (designated as 4-G-lines), two lines expressing the starch-degrading enzymes (3-G-lines), and five empty vector control (NPTII) lines were obtained. During in vitro growth and 3 weeks after transplanting to soil, the performance of the transgenic lines expressing multiple genes was not different from the non-transgenic control plants (Figure 2F). These results may suggest that expression of multiple hyperthermophilic archaeal enzymes does not interfere with normal plant growth and development due to these enzymes’ extremely low catalytic activities at ambient temperature.
Initial characterization of the transgenic lines was based on starch-degrading activity of enzymes extracted from fibrous roots. Plants were propagated in vitro on agar media and allowed to grow for at least 1 month until sufficient fibrous roots were obtained for protein extraction. Total protein from the roots was extracted and purified by heating, and used for saccharification of purified potato starch in vitro. The amount of glucose produced due to starch-hydrolyzing activity of the recombinant enzyme cocktails was determined using the DNS colorimetric assay.
As shown in Figure 3, statistically significant difference (p < 0.0001) was observed between the starch hydrolyzing activities of heat-purified (80°C for 25 min) protein extracts from fibrous roots of transgenic and that of control lines (WT and NPTII). Protein extracts from the majority of transgenic lines released more glucose from starch than those from wild-type and empty vector control lines. The enzyme cocktail from 4-G-L62 coexpressing the three starch-hydrolyzing enzymes (α-amylase, amylopullulanase, and glucoamylase) and E. coli AGPase showed the highest starch-degrading activity, releasing about sixfold more glucose than the wild-type control. Likewise, 3-G-L3, which expresses only the starch-hydrolyzing enzymes, released about fourfold more glucose than the wild-type control. This finding suggests that hyperthermophilic archaeal enzymes maintain their activity and thermostability when heterologously expressed in cassava roots.
FIGURE 3. Starch-degradation assay as determined by the DNS method. (A) Increase in color intensity of the DNS reagent with enzyme cocktail from wild-type (WT, non-transgenic), empty vector (NPTII) and transgenic lines (3-G, expressing the starch-degrading enzymes; 4-G, expressing the starch-degrading enzymes and AGPase). (B) Hydrolysis of purified potato starch by enzyme cocktail extracted from fibrous roots of control and transgenic cassava lines. The resulting reducing sugar/glucose was quantified by DNS colorimetric assay method using a maltose solution as the standard. Color intensity of DNS reagent was measured at 540 nm. Bars represent mean ± standard errors (SE) based on three biological replicates. Bars bearing the same letter are not statistically different (p < 0.05).
Characterization of Greenhouse- Established Transgenic Plants
To further characterize the multigene expressing lines, in vitro propagated transgenic and non-transgenic control plants were transplanted to soil according to Zidenga et al. (2012) and transferred to the greenhouse. Plants were grown under natural light from April to December.
While the transgenic lines developed normal shoots and fibrous roots on MS media in controled growth chambers and during the first 3 weeks after transplanting to soil (Figure 2F), we observed variations after greenhouse establishment in a later stage for most transgenic lines (Figure 4 and Table 1). The growth of 4-G-lines (4-G-L2 and 4-G-L53) and 3-G-L11 was severely stunted while 4-G-L34 and 4-G-L43 produced creeping shoots with weak stems (Figures 4A–D). This phenotypic alteration could be due to somaclonal variation as previously reported in cassava (Raemakers et al., 2001; Taylor et al., 2012). Since the FEC used in this study was about 4-months old, it is likely that the observed growth abnormality of the plants could be due to the age of the FEC. This is consistent with Taylor et al. (2012) who suggested that transformation of FEC at an optimum stage (about 9 weeks after induction), and timely subculturing of transgenic tissues could reduce the formation of abnormal plants. The observed growth abnormality could be due to genetic or epigenetic variation (Machczynska et al., 2015).
FIGURE 4. Characterization of transgenic lines grown in the greenhouse. (A–D) Transgenic lines exhibiting abnormal growth characteristics as compared to the non-transgenic control. (E) Greenhouse-established 6-month old plants. (F) Comparison of WT, 3-G-L3, and 4-G-L62 biomass. (G) Tubers harvested from 7-month-old WT and 4-G-L62. (H) Tubers pulverized under liquid N2 for protein or RNA extraction, or analysis of starch self-processing ability.
Initiation of tuberous root was observed after 4-month of growth in the soil. The frequency and time of tuber initiation varied among the greenhouse-established lines. The plants were allowed to grow for 7 months. Line 3-G-L3 (expressing only the archaeal genes) produced smaller shoot and tuber biomass than both the wild-type and 4-G-lines after 7 months (Table 1). Line 4-G-L62 whose root enzyme extract produced the highest amount of glucose (Figure 3) also produced more shoot and tuber biomass (Figures 4F,G and Table 1) than the control, 3-G-L3, and 4-G-L29. Interestingly, the average number of tubers per plant, which is positively correlated with yield (Edet et al., 2015), was also significantly higher in 4-G-L62 than the other transgenic lines and in the non-transformed control (Figure 4G and Table 1). This yield increase in 4-G-L62 is likely due to coexpression of modified glgC (G336D) (Table 2) in line with increased AGPase activity in 4-G-L62 (Figure 5E). However, the increase in storage root observed in transgenic TMS60444 in this study is lower than what was reported in glgC expressing TMS 71173 (Ihemere et al., 2006), which could be due to genotypic difference or growth conditions, for example, facilities for controlling light, temperature and humidity. Cassava growth and photosynthetic assimilate partitioning are largely influenced by environmental factors such as temperature and photoperiod (El-Sharkawy, 1993; Rao et al., 2016).
TABLE 2. Expression of archaeal (Amy, Pull, and Gluco) and bacterial (glgC) genes in tubers of wild-type (WT) and multigene-expressing lines as determined by qPCR.
FIGURE 5. Carbohydrate-hydrolyzing activity of enzyme cocktails extracted from tubers of WT or transgenic (4-G-L62) cassava as determined by the glucose oxidase–peroxidase assay. Degradation of starch (A), pullulan (B), maltose (C) and maltotriose (D) by heat-purified enzyme cocktail from WT (blue bar) and 4-G-L62 (red bar). (i) Increase in color intensity due to the release of glucose from each carbohydrate. (ii) Amount of glucose released from each polysaccharide as determined from a change in color intensity measured at 540 nm using glucose as a standard. (E) AGPase activity assay in roots of WT (blue bar) and 4-G-L62 (red bar) based on an increase in absorbance at OD340 due to NADPH formation generation. Bars represent mean ± SE based on four replicates. Bars bearing the same letter are not statistically different (p < 0.05).
Archaeal and Bacterial Genes Are Expressed in Transgenic Cassava Tubers
Given that the enzyme cocktail extracted from fibrous roots of transgenic lines significantly increased saccharification of starch in vitro as compared to the non-transgenic control (Figure 3), we studied transcript accumulation in the tubers to confirm expression of the archaeal and bacterial genes in transgenic storage roots. Total RNA was isolated from the inner part of the tuber (flesh, the part of the storage root used for starch extraction) of WT, 4-G-L29, and 4-G-L62. Transcript accumulation was quantified by quantitative RT-PCR.
As shown in Table 2, a significant increase in α-amylase, amylopullulanase, glucoamylase, and glgC transcript levels was observed in tubers of 4-G-L62 as compared to both WT and 4-G-L29. Amylopullulanase and glucoamylase transcript levels were also higher in 4-G-L29 than in the control, while no increase in the transcripts of α-amylase, and only a slight increase in that of glgC were observed in 4-G-L29. The observed increase in transcript accumulation of the archaeal genes is consistent with the higher starch degradation activity of the enzyme cocktail extracted from fibrous roots of transgenic lines as compared to the control (Figure 3). This correlation is more evident for amylase and amylopullulanase since expression of the two enzymes was higher in line 4-G-L62, which released the highest amount of glucose. Likewise, increased glgC transcript level in 4-G-L62 correlates well with enhanced tuber and overall biomass yield (Table 1 and Figure 4G).
A Cocktail of Enzymes from Transgenic Tubers Degrade Carbohydrates in Vitro
To determine the activity of the individual archaeal enzymes in the root-derived enzyme cocktail, we analyzed the activity of the enzyme cocktails extracted from tubers of 4-G-L62 and the WT on purified starch, pullulan, and maltotriose/maltose, chosen as substrates for amylase, amylopullulanase, and glucoamylase, respectively.
The enzyme cocktail extracted from 4-G-L62 showed increased saccharification of all the tested carbohydrates (Figure 5). The amount of glucose released from starch by the 4-G-L62 enzyme cocktail was at least 10-fold higher than the enzyme cocktail from WT control (80 mg/L vs. 8 mg/L) (Figure 5A). Likewise, purified proteins from 4-G-L62 released significantly more glucose from pullulan, which is a suitable substrate for amylopullulanase, than the proteins extracted from the WT. The activity of enzymes from 4-G-L62 on maltose and maltotriose (substrates of glucoamylase) was also higher than the WT (Figure 5B), suggesting that the archaeal enzymes are expressed in their active forms. The basal saccharification level observed for WT as compared to the control without a substrate could be due to substrate hydrolysis by heat.
Production of Hyperthermophilic Starch-Degrading Enzymes in Transgenic Tubers Enhance Starch Self-Processing Ability
Given that the primary goal of this study was to develop cassava with starch self-processing ability by directly producing thermostable starch-hydrolyzing enzymes in the tuber, we compared the self-processing ability of tubers from 4-G-lines to that of the non-transformed control. Starch self-processing ability was determined based on total reducing sugar content before and after heat incubation.
As shown in Figure 6, initial total reducing sugar content in the tubers was not significantly different between the transgenic lines and the wild-type control, suggesting that the thermophilic enzymes do not degrade starch in the intact tubers. After incubation of the tuber biomass for 21-h at 85°C, multigene expressing 4-G-L62 released significantly more reducing sugar than the WT, while 4-G-L29 was not different from the WT (Figures 6A,B). Lack of self-processing ability in 4-G-L29 could be due to a low level of enzyme accumulation, consistent with a low level of gene expression as compared to 4-G-L62. In particular, the α-amylase transcript did not accumulate in the tubers of 4-G-L29 beyond the background level (Table 2). The reducing sugar yield was three times higher in 4-G-L62 than in 4-G-L29 and the WT (Figures 6A,C). These findings establish that accumulation of starch hydrolyzing enzymes induce tuber self-processing ability.
FIGURE 6. Starch self-processing ability of WT and transgenic (4-G-L29 and 4-G-L62) tubers. (A) Increase in color intensity due to the release of glucose from storage roots after incubation at 85°C for 21-h. (B) Reducing sugar content determined from a change in color intensity of DNS reagent before (blue bars) and after (red bars) incubation at 85°C measured at 540 nm using maltose as the standard. (C) Fold increase in reducing sugar after 21-h incubation at 85°C as compared to initial sugar content. Bars represent mean ± SE based on four biological replicates. Bars bearing the same letter are not statistically different (p < 0.05).
Recombinant Archaeal Enzymes Maintain Their Catalytic Activity and Thermostability
In this study, expression of a suite of archaeal starch-degrading enzymes was targeted to the roots of cassava. To understand whether the enzymes are accumulated in their active forms, total protein was extracted from storage roots was used for saccharification of carbohydrates including starch, pullulan, maltose/maltotriose as substrates for amylase, amylopullulanase and glucoamylase, respectively.
The findings revealed that recombinant enzyme cocktails extracted from storage roots of the transgenic line showed significantly more carbohydrate-hydrolyzing activity, as compared to the cocktail extracted from the non-transgenic control (Figures 3, 5). This increase in carbohydrate saccharification establishes that the archaeal enzymes maintain their activity and thermostability when produced in cassava roots. The amount of glucose released from starch and pullulan by the 4-G-L62 enzyme cocktail was over 10-fold more than the enzyme cocktail from WT. Higher activity of the cocktail on starch and pullulan is likely due to higher expression of α-amylase and amylopullulanase (Table 2) as compared to glucoamylase, since amylopullulanase is a bifunctional enzyme possessing both pullulanase and amylase activities (Kim et al., 1993; Chiang et al., 2005). The amount of glucose released from both maltose and maltotriose was much lower than those released from starch or pullulan, which is likely due to low level of gene expression. Taken together, our findings establish that active forms of a suite of thermophilic carbohydrate-degrading enzymes can be produced in planta. Note, though, that this conclusion is based on the cumulative activity of the enzyme mix, and that the specific activity of each recombinant enzyme needs to be determined separately in the future. Furthermore, since the assays were performed at 85°C, a temperature typically used for starch hydrolysis, assay conditions such as temperature, pH, buffer, and ions need to be optimized for the three enzymes in the cocktail.
The observed increase in saccharification of several carbohydrates due to targeted accumulation of the archaeal enzymes in the tubers suggests a potential of the tubers for commercial-scale production of high-value enzymes because tubers have specialized molecular environment for product stability as previously reported in potato (Sparrow et al., 2007). Targeting of the proteins to the apoplast may further enhance enzyme production and improve saccharification efficiency. For example, targeting expression of hydroxynitrile lyase to the apoplast under the control of Patatin promoter has been shown to enhance protein accumulation by threefold in cassava (Narayanan et al., 2011). Likewise, accumulation of the enzymes may also be targeted to the endoplasmic reticulum (ER), which has several advantages including being a suitable environment for correct folding and disulfide bridge formation by the ER chaperones and enzymes (Schillberg et al., 1999; Wilkinson et al., 2005), and low levels of proteolytic activity in the ER lumen as compared to the cytosol.
Enhanced Self-Processing Ability of Transgenic Tubers
To validate starch hydrolyzing activity of the recombinant hyperthermophilic enzymes directly in the tubers, the tubers were pulverized, incubated at 85°C overnight to activate the enzymes, and the total reducing sugar yield in the hydrolysate was quantified. The result shows that transgenic tubers release three-times more reducing sugars as compared to the wild-type control, suggesting that the recombinant enzymes can hydrolyze starch in situ. A number of thermostable starch-hydrolyzing genes from bacteria have been expressed in transgenic plants. For example, coexpression of Bacillus stearothermophilus bifunctional α-amylase and glucose isomerase in potato tuber (Solanum tuberosum) (Beaujean et al., 2000), B. stearothermophilusα-amylase in transgenic rice seeds (Xu et al., 2008), Thermotoga maritima α-amylase in sweet potato (Ipomoea batatas) (Santa-Maria et al., 2011), and T. thermohydrosulfuricus truncated amylopullulanase in corn (Nahampun et al., 2013) induced starch saccharification at high temperatures. Interestingly, in all these studies the transgenes have not been shown to affect plant development and metabolism, and harvestable yield. However, cases in which thermostable starch-degrading enzymes modulated plant composition have also been reported, for example, increasing the expression level of thermostable Thermoanaerobacter ethanolicus amylopullulanase in rice modified the starch composition by reducing the amylose content (Chiang et al., 2005).
Recombinant hyperthermophilic archaeal glycosyl hydrolases have been produced in bacteria and their carbohydrate degrading activity has been established (Dong et al., 1997a,b; Kim et al., 2004; Lee et al., 2006). Only a few archaeal proteins have been successfully expressed in planta, though, primarily for studying the feasibility of recombinant protein production (Montalvo-Rodriguez et al., 2000; Conley et al., 2011; Klose et al., 2012) or stress tolerance (Im et al., 2009). To our knowledge, this study reports for the first time on heterologous expression of archaeal genes in plants for an in situ activation to enhance starch self-processing ability. This novel approach should lead to a sustainable route for starch conversion to sugars for biofuels and other industrial products. Given that cassava is the highest yielding starchy feedstock (Ziska et al., 2009; El-Sharkawy, 2012), and the most resilient crop to environmental variability (El-Sharkawy, 2012; Jarvis et al., 2012), transgenic cassava with improved traits will attract tremendous commercial interest.
Potential Applications of This Study
In this study we expressed multiple transgenes in cassava tubers to simultaneously induce starch autohydrolysis and increase starch yield. Trait improvement through traditional breeding is challenging in crops such as cassava, sugarcane (Saccharum spp.), and the biofuel grass miscanthus (Miscanthus spp.) (McGarry et al., 2017). The vector system used in this study offers a platform for multiple gene insertion into these crops to engineer traits that require a suite of enzymes, for example, polysaccharide conversion to sugars for biofuels or other industrial purposes (as in this study), disease resistance (Tricoli et al., 1995; Beyene et al., 2017), pest resistance (Santamaria et al., 2012), enhanced vitamin content (Ye et al., 2000; Lu et al., 2013; Pons et al., 2014), molecular farming of industrial and pharmaceutical products (De Martinis et al., 2016), and simultaneous improvement of two or more traits.
Studies have shown that cassava-derived bioethanol has economic advantages over other starchy feedstocks including corn (Dai et al., 2006; Ziska et al., 2009). The annual bioethanol yield of cassava is three-times that of corn (Wang, 2002; Ziska et al., 2009). Therefore, introducing starch self-processing trait in cassava will further increase its demand as bioethanol feedstock, and will have great economic and environmental benefits. Starch processing to sugar for bioethanol production or in food processing industries is a multistep process that requires the application of commercial enzymes such as α-amylase and glucoamylase. Considering the cost of large-scale enzyme production and applications to bioprocessing, the current approach whereby starch-degrading enzymes are directly produced and accumulated in the plant starchy biomass, and are heat-activated in situ to break down starch, would tremendously cut the cost of starchy feedstock processing. This technology has the potential to attract commercial interest, as observed for the genetically modified corn developed by Syngenta (Enogen) that expresses microbial α-amylase for ethanol production (Johnson et al., 2009). Furthermore, the current approach reduces energy/labor cost and increases process efficiency in many ways; (1) the archaeal hyperthermophilic enzymes with temperature optima above 90°C (Dong et al., 1997a,b; Jorgensen et al., 1997; Kim et al., 2004) will circumvent the need for the expensive cooling step normally required before saccharification step in the case of fungal glucoamylases with lower temperature optima (50–70°C) (Norouzian et al., 2006), (2) transgenic storage roots may be processed directly and require less water for grinding, omitting further processing such as drying which requires space, and is labor intensive, (3) provided that the recombinant enzymes are accumulated at higher levels, the transgenic biomass may be blended with non-transgenic biomass to reduce enzyme cost and the need to allocate more land for genetically engineered plants. Moreover, glucose produced by this potentially cheap route can serve as a feedstock for a plethora of industrial products. A number of products including succinic acid, itaconic acid, adipic acid, 3-hydroxypropanoic acid/aldehyde, isoprene/farmesene, glutamic acid and aspartic acid can be produced from glucose through biological fermentation (Sanchez et al., 2005; de Jong et al., 2010). Likewise, glucose can be chemically transformed to useful products such as sorbitol, furfural, glucaric acid, hydroxymethyl furfural and levulinic acid (Kusserow et al., 2003; de Jong et al., 2010).
In this study, we report on the production of a suite of hyperthermophilic archaeal enzymes that lead to enhanced starch self-processing ability directly in the transgenic biomass. However, to realize the full potential of this technology, recombinant enzyme production should be enhanced, for example, by codon-optimizing the native gene sequences, by using stronger promoters, translation enhancers, and/or by targeting the proteins to cellular compartments, such as apoplast and ER. All of these approaches have been shown to increase protein expression and storage levels elsewhere. To compensate for any potential harvestable yield penalty due to expression of multiple genes, there may be a need to co-express genes that enhance biomass/starch yield. Furthermore, future research will characterize the transgenic lines under well-controlled growth conditions, and establish optimum reaction conditions for the enzyme cocktail, including processing media, temperature, pH and cofactor additions. Moreover, the fermentability of the saccharification product by yeast needs to be investigated. Overall, this study reports on a novel approach for sustainable and inexpensive route to produce fermentable sugars from cassava, which can easily be transferred to other low-input tuber/root crops such as sweet potato, Solanum genus potato, and other tropical root crops including yam, taro, and arrowroot.
BH, AL-O, CV, and MA conceived the study and secured federal grant to support research activities. AL-O, MA, FP, CV, and BH designed all experiments. AL-O, BH, FP, JJ, SC, C-HW, and ED conducted experiments and collected data. BH, SC, FP, and C-HW cloned the starch hydrolyzing archaeal genes from Pyrococcus furiosus. AL-O cloned glucoamylase from Sulfolobus solfataricus and E. coli glgC, generated all the transgenic cassava plants, performed all biochemical characterization of these lines, and wrote the manuscript. AL-O and BH analyzed the data, designed and constructed and characterized all the expression vectors utilized to recover transgenic cassava lines. AL-O, JJ, and ED performed expression studies of Pyrococcus furiosus amylopullulanase and α-amylase in E. coli. CV, BH, MA, SC, FP, and C-HW edited the manuscript. All authors read and approved the manuscript.
Cassava-related research was sponsored by USDA-AFRI and USDA NIFA CBG Grants (2008-35504-04557 and 2011-38821-30974 to BH). Other projects in the laboratory and AL-O were supported by USDA-NIFA CBG (# 2014-38821-22417) and Evans-Allen Grant (#6011-220017-21). The research concerning the hyperthermophilic archaea was supported by the Department of Energy (DE-FG05-95ER20175 to MA).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Mr. Ifunanya Nwogbaga for technical assistance, and members of the USDA-ARS Lab including Drs. David Kingsley and Garry Richards for constructive suggestions and allowing AL-O to use their lab facilities as needed; Dr. Vitaly Citovsky Lab (Stony Brook University) for pSATs and pPZP-RCS2 vectors; David R. Rockhold and William R. Belknap (USDA-ARS) for kindly providing Patatin I and GBSS promoters, Dr. Tony Romeo (University of Florida) for providing pO12 plasmid containing the E. coli glgC sequence, Dr. Mayavan Subramani (Delaware State University) and Dr. Getu Beyene (Donald Danforth Plant Science) for technical advice, Donald Danforth Plant Science center for providing wild-type cassava cultivar TMS60444, and Dr. Dyremple Marsh (Delaware State University) for supporting our work with essential facilities and supplies as needed.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.00192/full#supplementary-material
FIGURE S1 | Validation of enzyme activity of recombinant Pyrococcus furiosus α-amylase and amylopullulanase. (A) The release of glucose from starch by α-amylase. (B) The release of glucose from pullulan by amylopullulanse. Values represent means ± SE based on three replicates.
TABLE S1 | Primers used for cloning root specific promoters Patatin I and GBSS, archaeal starch hydrolyzing enzymes and E. coli Glgc into pSAT modular vectors for subsequent assembly into pPZP-NPTII binary vector. The sequence of each primer includes 5′ extra nucleotides (italicized), restriction enzyme site (upper case letters and underlined) and a sequence specific to the target promoter or gene (lower case non-italic).
TABLE S2 | Types of media and composition used for cassava tissue culture. MS (MS with vitamins; Murashige and Skoog, 1962), GD (Greshoff and Doy, 1972) composition is based on phytotechnology laboratories (https://phytotechlab.com/), SH (Schenk and Hildebrandt, 1972), BAP (6-Benzylaminopurine), NAA (Naphthalene acetic acid). Media were solidified by phytagel (0.22%) for CBM and Noble Agar (0.8%) for all other media. pH (5.8).
TABLE S3 | Primers used for qPCR to study expression of archaeal starch hydrolyzing enzymes and E. coli glgc from control and transgenic cassava tubers. Cassava α-tubulin gene (TC3055) was used as an internal control.
Angelov, M. N., Sun, J., Byrd, G. T., Brown, R. H., and Black, C. C. (1993). Novel characteristics of cassava, Manihot esculenta Crantz, a reputed C3-C4 intermediate photosynthesis species. Photosynth. Res. 38, 61–72. doi: 10.1007/BF00015062
Ballicora, M. A., Iglesias, A. A., and Preiss, J. (2003). ADP-glucose pyrophosphorylase, a regulatory enzyme for bacterial glycogen synthesis. Microbiol. Mol. Biol. Rev. 67, 213–225. doi: 10.1128/MMBR.67.2.213-225.2003
Beaujean, A., Ducrocq-Assaf, C., Sangwan, R. S., Lilius, G., Bülow, L., and Sangwan-Norreel, B. S. (2000). Engineering direct fructose production in processed potato tubers by expressing a bifunctional alpha-amylase/glucose isomerase gene complex. Biotechnol. Bioeng. 70, 9–16. doi: 10.1002/1097-0290(20001005)70:1<9::AID-BIT2>3.0.CO;2-7
Beyene, G., Chauhan, R. D., Ilyas, M., Wagaba, H., Fauquet, C. M., Miano, D., et al. (2017). A virus-derived stacked RNAi construct confers robust resistance to cassava brown streak disease. Front. Plant Sci. 7:2052. doi: 10.3389/fpls.2016.02052
Bridger, S. L., Lancaster, W. A., Poole, F. L. II, Schut, G. J., and Adams, M. W. (2012). Genome sequencing of a genetically tractable Pyrococcus furiosus strain reveals a highly dynamic genome. J. Bacteriol. 194, 4097–4106. doi: 10.1128/JB.00439-12
Bull, S. E., Owiti, J. A., Niklaus, M., Beeching, J. R., Gruissem, W., and Vanderschuren, H. (2009). Agrobacterium-mediated transformation of friable embryogenic calli and regeneration of transgenic cassava. Nat. Protoc. 4, 1845–1854. doi: 10.1038/nprot.2009.208
Ceballos, H., Kawuki, R. S., Gracen, V. E., Yencho, G. C., and Hershey, C. H. (2015). Conventional breeding, marker-assisted selection, genomic selection and inbreeding in clonally propagated crops: a case study for cassava. Theor. Appl. Genet. 128, 1647–1667. doi: 10.1007/s00122-015-2555-4
Chavarriaga-Aguirre, P., Brand, A., Medina, A., Prías, M., Escobar, R., Martinez, J., et al. (2016). The potential of using biotechnology to improve cassava: a review. In Vitro Cell Dev. Biol. Plant. 52, 461–478. doi: 10.1007/s11627-016-9776-3
Chetty, C. C., Rossin, C. B., Gruissem, W., Vanderschuren, H., and Rey, M. E. (2013). Empowering biotechnology in southern Africa: establishment of a robust transformation platform for the production of transgenic industry-preferred cassava. New Biotechnol. 30, 136–143. doi: 10.1016/j.nbt.2012.04.006
Chiang, C. M., Yeh, F. S., Huang, L. F., Tseng, T. H., Chung, M. C., Wang, C. S., et al. (2005). Expression of a bifunctional and thermostable amylopullulanase in transgenic rice seeds leads to autohydrolysis and altered composition of starch. Mol. Breed. 15, 125–143. doi: 10.1007/s11032-004-3919-7
Conley, A. J., Zhu, H., Le, L. C., Jevnikar, A. M., Lee, B. H., Brandle, J. E., et al. (2011). Recombinant protein production in a variety of Nicotiana hosts: a comparative analysis. Plant Biotechnol. J. 9, 434–444. doi: 10.1111/j.1467-7652.2010.00563.x
Cooke, R., and Traore, G. (2010). “Global consultation on cassava as a potential bioenergy crop,” in International Fund for Agricultural Development, eds E. Kueneman, V. Raswant, N. Lutaladio, and R.Cooke (Rome: IFAD).
Dai, D., Hu, Z., Pu, G., Li, H., and Wang, C. (2006). Energy efficiency and potentials of cassava fuel ethanol in Guangxi region of China. Energy Convers. Manag. 47, 1686–1699. doi: 10.1016/j.enconman.2005.10.019
De Martinis, D., Rybicki, E. P., Fujiyama, K., Franconi, R., and Benvenuto, E. (2016). Editorial: plant molecular farming: fast, scalable, cheap, sustainable. Front. Plant Sci. 7:1148. doi: 10.3389/fpls.2016.01148
De Souza, A. P., Massenburg, L. N., Jaiswal, D., Cheng, S., Shekar, R., and Long, S. P. (2016). Rooting for cassava: insights into photosynthesis and associated physiology as a route to improve yield potential. New Phytol. 213, 50–65. doi: 10.1111/nph.14250
Dong, G., Vieille, C., Savchenko, A., and Zeikus, J. G. (1997a). Cloning, sequencing, and expression of the gene encoding extracellular alpha-amylase from Pyrococcus furiosus and biochemical characterization of the recombinant enzyme. Appl. Environ. Microbiol. 63, 3569–3576.
Dong, G., Vieille, C., and Zeikus, J. G. (1997b). Cloning, sequencing, and expression of the gene encoding amylopullulanase from Pyrococcus furiosus and biochemical characterization of the recombinant enzyme. Appl. Environ. Microbiol. 63, 3577–3584.
Duffner, F., Bertoldo, C., Andersen, J. T., Wagner, K., and Antranikian, G. (2000). A new thermoactive pullulanase from Desulfurococcus mucosus: cloning, sequencing, purification, and characterization of the recombinant enzyme after expression in Bacillus subtilis. J. Bacteriol. 182, 6331–6338. doi: 10.1128/JB.182.22.6331-6338.2000
Edet, M. A., Tijani-Eniola, H., Lagoke, S. T. O., and Tarawali, G. (2015). Relationship of cassava growth parameters with yield, yield related components and harvest time in ibadan, Southwestern Nigeria. J. Nat. Sci. Res. 5, 87-93.
El-Sharkawy, M. A. (1993). Drought-tolerant cassava for Africa, Asia, and Latin America: breeding projects work to stabilize productivity without increasing pressures on limited natural resources. Bioscience 43, 441–451. doi: 10.2307/1311903
El-Sharkawy, M. A., and Cock, J. H. (1984). Water use efficiency of cassava. I. Effects of air humidity and water stress on stomatal conductance and gas exchange. Crop Sci. 24, 497–502. doi: 10.2135/cropsci1984.0011183X002400030017x
Fermont, A. M., van Asten, P. J. A., Tittonell, P., van Wijk, M. T., and Giller, K. E. (2009). Closing the cassava yield gap: an analysis from smallholder farms in East Africa. Field Crops Res. 112, 24–36. doi: 10.1016/j.fcr.2009.01.009
Goderis, I. J., De Bolle, M. F., François, I. E., Wouters, P. F., Broekaert, W. F., Cammue, B. P., et al. (2002). A set of modular plant transformation vectors allowing flexible insertion of up to six expression units. Plant Mol. Biol. 50, 17–27. doi: 10.1023/A:1016052416053
Ihemere, U., Arias-Garzon, D., Lawrence, S., and Sayre, T. (2006). Genetic modification of cassava for enhanced starch production. Plant Biotechnol. J. 4, 453–465. doi: 10.1111/j.1467-7652.2006.00195.x
Im, Y. J., Ji, M., Lee, A., Killens, R., Grunden, A. M., and Boss, W. F. (2009). Expression of Pyrococcus furiosus superoxide reductase in Arabidopsis enhances heat tolerance. Plant Physiol. 151, 893–904. doi: 10.1104/pp.109.145409
Jorgensen, K., Bak, S., and Busk, P. K. (2005). Cassava plants with a depleted cyanogenic glucoside content in leaves and tubers. Distribution of cyanogenic glucosides, their site of synthesis and transport, and blockage of the biosynthesis by RNA interference technology. Plant Physiol. 139, 363–374. doi: 10.1104/pp.105.065904
Jorgensen, S., Vorgias, C. E., and Antranikian, G. (1997). Cloning, sequencing, characterization, and expression of an extracellular alpha-amylase from the hyperthermophilic archaeon Pyrococcus furiosus in Escherichia coli and Bacillus subtilis. J. Biol. Chem. 272, 16335–16342. doi: 10.1074/jbc.272.26.16335
Karlstrom, A., Calle, F., Salazar, S., Morante, N., Dufour, D., and Ceballos, H. (2016). Biological implications in cassava for the production of amylose-free starch: impact on root yield and related traits. Front. Plant Sci. 7:604. doi: 10.3389/fpls.2016.00604
Kim, C. H., Choi, H. I., and Lee, D. S. (1993). Pullulanases of alkaline and broad pH range from a newly isolated alkalophilic Bacillus sp. S-1 and a Micrococcus sp. Y-1. J. Ind. Microbiol. 12, 48–57. doi: 10.1007/BF01570128
Kim, M. S., Park, J. T., Kim, Y. W., Lee, H. S., Nyawira, R., and Shin, H. S. et al. (2004). Properties of a novel thermostable glucoamylase from the hyperthermophilic archaeon Sulfolobus solfataricus in relation to starch processing. Appl. Environ. Microbiol. 70, 3933–3940. doi: 10.1128/AEM.70.7.3933-3940.2004
Klose, H., Roeder, J., Girfoglio, M., Fischer, R., and Commandeur, U. (2012). Hyperthermophilic endoglucanase for in planta lignocellulose conversion. Biotechnol. Biofuels 5:63. doi: 10.1186/1754-6834-5-63
Kulichikhin, K., Mukherjee, S., and Ayele, B. (2016). Extraction and assays of ADP-glucose pyrophosphorylase, soluble starch synthase and granule bound starch synthase from wheat (Triticum aestivum L.) Grains. Bio Protocol 6:e1929. doi: 10.21769/BioProtoc.1929
Lee, H. S., Shockley, K. R., Schut, G. J., Conners, S. B., Montero, C. I., Johnson, M. R., et al. (2006). Transcriptional and biochemical analysis of starch metabolism in the hyperthermophilic archaeon Pyrococcus furiosus. J. Bacteriol. 188, 2115–2125. doi: 10.1128/JB.188.6.2115-2125.2006
Li, K. T., Moulin, M., Mangel, N., Albersen, M., Verhoeven-Duif, N. M., Ma, Q., et al. (2015). Increased bioavailable vitamin B-6 in field-grown transgenic cassava for dietary sufficiency. Nat. Biotechnol. 33, 1029–1032. doi: 10.1038/nbt.3318
Lopes, M. L., Paulillo, S. C., Godoy, A., Cherubin, R. A., Lorenzi, M. S., Giometti, F. H., et al. (2017). Ethanol production in Brazil: a bridge between science and industry. Braz. J. Microbiol. 47, 64–76. doi: 10.1016/j.bjm.2016.10.003
Lu, Y., Rijzaani, H., Karcher, D., Ruf, S., and Bock, R. (2013). Efficient metabolic pathway engineering in transgenic tobacco and tomato plastids with synthetic multigene operons. Proc. Natl. Acad. Sci. U.S.A. 110, E623–E632. doi: 10.1073/pnas.1216898110
Machczynska, J., Zimny, J., and Bednarek, P. T. (2015). Tissue culture-induced genetic and epigenetic variation in triticale (Triticosecale spp. Wittmack ex A. Camus 1927) regenerants. Plant Mol. Biol. 89, 279–292. doi: 10.1007/s11103-015-0368-0
McGarry, R. C., Klocko, A. L., Pang, M., Strauss, S. H., and Ayre, B. G. (2017). Virus-induced flowering: an application of reproductive biology to benefit plant research and breeding. Plant Physiol. 173, 47–55. doi: 10.1104/pp.16.01336
Mehta, D., and Satyanarayana, T. (2013). Dimerization mediates thermo-adaptation, substrate affinity and transglycosylation in a highly thermostable maltogenic amylase of Geobacillus thermoleovorans. PLoS One 8:e73612. doi: 10.1371/journal.pone.0073612
Mei, C., and Rakhmatov, C. (2014). “Advances in the genetic manipulation of cellulosic bioenergy crops for bioethanol production,” in Biological Conversion of Biomass for Fuels and Chemicals. Explorations from Natural Utilization Systems, eds J. Z. Sun, S. Y. Ding, and J. Doran-Peterson (Cambridge: Royal Society of Chemistry) 10, 53–82.
Meyer, C. R., Bork, J. A., Nadler, S., Yirsa, J., and Preiss, J. (1998). Site-directed mutagenesis of a regulatory site of Escherichia coli ADP glucose pyrophophorylase: the role of residue 336 in allosteric behavior. Arch. Biochem. Biophys. 353, 152–159. doi: 10.1006/abbi.1998.0648
Montalvo-Rodriguez, R., Haseltine, C., Huess-LaRossa, K., Clemente, T., Soto, J., and Staswick, P. (2000). Autohydrolysis of plant polysaccharides using transgenic hyperthermophilic enzymes. Biotechnol. Bioeng. 70, 151–159. doi: 10.1002/1097-0290(20001020)70:2<151::AID-BIT4>3.0.CO;2-D
Nahampun, H. N., Lee, C. J., Jane, J.-L., and Wang, K. (2013). Ectopic expression of bacterial amylopullulanase enhances bioethanol production from maize grain. Plant Cell Rep. 32, 1393–1405. doi: 10.1007/s00299-013-1453-y
Narayanan, N. N., Ihemere, U., Ellery, C., and Sayre, R. T. (2011). Overexpression of hydroxynitrile lyase in cassava roots elevates protein and free amino acids while reducing residual cyanogen levels. PLoS One 6:e21996. doi: 10.1371/journal.pone.0021996
Pons, E., Alqueìzar, B., Rodríguez, A., Ramón, D., Rodrigo, M. J., Zacarías, L., et al. (2014). Metabolic engineering of β-carotene in orange fruit increases its in vivo antioxidant properties. Plant Biotechnol. J. 12, 17–27. doi: 10.1111/pbi.12112
Que, Q., Chilton, M. D. M., de Fontes, C. M., He, C., Nuccio, M., Zhu, T., et al. (2010). Trait stacking in transgenic crops: challenges and opportunities. GM Crops 4, 220–229. doi: 10.4161/gmcr.1.4.13439
Rao, N. K. S., Laxman, R. H., and Shivashankara, K. S. (2016). “Physiological and morphological responses of horticultural crops to abiotic stress,” in Abiotic Stress Physiology of Horticultural Crops, eds N. K. Srinivasa Rao, K. S. Shivashankara, and R. H. Laxman (New Delhi: Springer). doi: 10.1007/978-81-322-2725-0
Richardson, T. H., Tan, X., Frey, G., Callen, W., Cabell, M., Lam, D., et al. (2002). A novel, high performance enzyme for starch liquefaction. Discovery and optimization of a low pH, thermostable alpha-amylase. J. Biol. Chem. 277, 26501–26507. doi: 10.1074/jbc.M203183200
Rudiger, A., Jorgensen P. L., and Antranikian, G. (1995). Isolation and characterization of a heat-stable pullulanase from the hyperthermophilic archaeon Pyrococcus woesei after cloning and expression of its gene in Escherichia coli. Appl. Environ. Microbiol. 61, 567–575.
Sanchez, A. M., Bennett, G. N., and San, K. Y. (2005). Efficient succinic acid production from glucose through overexpression of pyruvate carboxylase in an Escherichia coli alcohol dehydrogenase and lactate dehydrogenase mutant. Biotechnol. Prog. 21:358–365. doi: 10.1021/bp049676e
Santa-Maria, M. C., Yencho, C. G., Haigler, C. H., Thompson, W. F., Kelly, R. M., and Sosinski, B. (2011). Starch self-processing in transgenic sweet potato roots expressing a hyperthermophilic alpha-amylase. Biotechnol. Prog. 27, 351–359. doi: 10.1002/btpr.573
Santamaria, M. E., Cambra, I., Martinez, M., Pozancos, C., González-Melendi, P., and Grbi, V. (2012). Gene pyramiding of peptidase inhibitors enhances plant resistance to the spider mite Tetranychus urticae. PLoS One 7:e43011. doi: 10.1371/journal.pone.0043011
Schillberg, S., Zimmermann, S., Voss, A., and Fischer, R. (1999). Apoplastic and cytosolic expression of full size antibodies and antibody fragments in Nicotiana tabacum. Transgenic Res. 8, 255–263. doi: 10.1023/A:1008937011213
Sparrow, P. A., Irwin, J. A., Dale, P. J., Twyman, R. M., and Ma, J. K. (2007). Pharma-planta: road testing the developing regulatory guidelines for plant-made pharmaceuticals. Transgenic Res. 16, 147–161. doi: 10.1007/s11248-007-9074-2
Stark, D. M., Timmerman, K. P., Barry, G. F., Preiss, J., and Kishore, G. M. (1992). Regulation of the amount of starch in plant tissues by ADP-glucose pyrophosphorylase. Science 258, 287–292. doi: 10.1126/science.258.5080.287
Taylor, N., Gaitán-Solís, E., Moll, T., Trauterman, B., Jones, T., Pranjal, A., et al. (2012). A high-throughput platform for the production and analysis of transgenic cassava (Manihot esculenta) plants. Trop. Plant Biol. 5, 127–139. doi: 10.1007/s12042-012-9099-4
Taylor, N. J., Masona, M. V., Carcamo, R., Ho, T., Schopke, C., and Fauquet, C. M. (2001). Production of embryogenic tissues and regeneration of transgenic plants in cassava (Manihot esculenta Crantz). Euphytica 2001, 25–34. doi: 10.1023/A:1017526814789
Tiessen, A., Hendriks, J. H., Stitt, M., Branscheid, A., Gibon, Y., Farré, E. M., et al. (2002). Starch synthesis in potato tubers is regulated by post-transcriptional redox modification of ADP-glucose pyrophosphorylase: a novel regulatory mechanism linking starch synthesis to the sucrose supply. Plant Cell 14, 2191–2213. doi: 10.1105/tpc.003640
Tricoli, D. M., Careny, K. J., Russell, P. F., McMaster, J. R., Groff, D. W., Hadden, K. C., et al. (1995). Field evaluation of transgenic squash containing single or multiple virus coat protein gene constructs for resistance to cucumber mosaic virus, watermelon mosaic virus 2, and zucchini yellow mosaic virus. Nat. Biotechnol. 13, 1458–1464. doi: 10.1038/nbt1295-1458
Tzfira, T., Tian, G.-W., and Lacroix, B. (2005). pSAT vectors: a modular series of plasmids for autofluorescent protein tagging and expression of multiple genes in plants. Plant Mol. Biol. 57, 503–516. doi: 10.1007/s11103-005-0340-5
Vanderschuren, H., Alder, A., Zhang, P., and Gruissem, W. (2009). Dose-dependent RNAi-mediated geminivirus resistance in the tropical root crop cassava. Plant Mol. Biol. 70, 365–372. doi: 10.1007/s11103-009-9472-3
Vanderschuren, H., Moreno, I., Anjanappa, R. B., Zainuddin, I. M., and Gruissem, W. (2012). Exploiting the combination of natural and genetically engineered resistance to cassava mosaic and cassava brown streak viruses impacting cassava production in Africa. PLoS One 7:e45277. doi: 10.1371/journal.pone.0045277
Wang, W. (2002). “Cassava production for industrial utilization in the PRC-present and future perspective,” in Proceedings of the 7th Regional Cassava Workshop: Cassava Research and Development in Asia: Exploring New Opportunities for an Ancient Crop, Bangkok, 33–38.
Welsch, R., Arango, J., and Bar, C. (2010). Provitamin A accumulation in cassava (Manihot esculenta) roots driven by a single nucleotide polymorphism in a phytoene synthase gene. Plant Cell 22, 3348–3356. doi: 10.1105/tpc.110.077560
Wilkinson, B., Xiao, R., and Gilbert, H. F. (2005). A structural disulfide of yeast protein-disulfide isomerase destabilizes the active site disulfide of the N-terminal thioredoxin domain. J. Biol. Chem. 280, 11483–11487. doi: 10.1074/jbc.M414203200
Xu, J., Duan, X., Yang, J., Beeching, J. R., and Zhang, P. (2013). Enhanced reactive oxygen species scavenging by overproduction of superoxide dismutase and catalase delays postharvest physiological deterioration of cassava storage roots. Plant Physiol. 161, 1517–1528. doi: 10.1104/pp.112.212803
Xu, J., Yang, J., Duan, X., Jiang, Y., and Zhang, P. (2014). Increased expression of native cytosolic Cu/Zn superoxide dismutase and ascorbate peroxidase improves tolerance to oxidative and chilling stresses in cassava (Manihot esculenta Crantz). BMC Plant Biol. 14:208. doi: 10.1186/s12870-014-0208-4
Xu, X., Fang, J., Wang, W., Guo, J., Chen, P., Cheng, J., et al. (2008). Expression of a bacterial alpha-amylase gene in transgenic rice seeds. Transgenic Res. 17, 645–650. doi: 10.1007/s11248-007-9144-5
Xuan, T. D., Phuong, N. T., Khang, D. T., and Khan, T. D. (2015). Influence of sowing times, densities, and soils to biomass and ethanol yield of sweet sorghum. Sustainability 7, 11657–11678. doi: 10.3390/su70911657
Ye, X., Al-Babili, S., Kloti, A., Zhang, J., Lucca, P., and Beyer, P. I. (2000). Engineering the provitamin A (b-carotene) biosynthetic pathway into (carotenoid-free) rice endosperm. Science 287, 303–305. doi: 10.1126/science.287.5451.303
Zhao, S. S., Dufour, D., Sanchez, T., Ceballos, H., and Zhang, P. (2011). Development of waxy cassava with different biological and physico-chemical characteristics of starches for industrial applications. Biotechnol. Bioeng. 108, 1925–1935. doi: 10.1002/bit.23120
Zidenga, T., Leyva-Guerrero, E., Moon, H., Siritunga, D., and Sayre, R. (2012). Extending cassava root shelf life via reduction of reactive oxygen species production. Plant Physiol. 154, 1396–1407. doi: 10.1104/pp.112.200345
Ziska, L. H., Runion, G. B., Tomecek, M., Prior, S. A., Torbet, H. A., and Sicher, R. (2009). An evaluation of cassava, sweet potato and field corn as potential carbohydrate sources for bioethanol production in Alabama and Maryland. Biomass Bioeng. 33, 1503–1508. doi: 10.1016/j.biombioe.2009.07.014
Keywords: cassava, hyperthermophilic archaeal, glgC, starch self-processing, bioethanol, multigene-expression, starch-hydrolyzing enzymes
Citation: Ligaba-Osena A, Jones J, Donkor E, Chandrayan S, Pole F, Wu C-H, Vieille C, Adams MWW and Hankoua BB (2018) Novel Bioengineered Cassava Expressing an Archaeal Starch Degradation System and a Bacterial ADP-Glucose Pyrophosphorylase for Starch Self-Digestibility and Yield Increase. Front. Plant Sci. 9:192. doi: 10.3389/fpls.2018.00192
Received: 29 August 2017; Accepted: 01 February 2018;
Published: 26 February 2018.
Edited by:Zeng-Yu Wang, Noble Research Institute, LLC, United States
Reviewed by:Peng Zhang, Institute of Plant Physiology and Ecology, Shanghai Institutes for Biological Sciences (CAS), China
Alberto A. Iglesias, National University of the Littoral, Argentina
Copyright © 2018 Ligaba-Osena, Jones, Donkor, Chandrayan, Pole, Wu, Vieille, Adams and Hankoua. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bertrand B. Hankoua, email@example.com