CC-type glutaredoxin, MeGRXC3, associates with catalases and negatively regulates drought tolerance in cassava (Manihot esculenta Crantz).

Glutaredoxins (GRXs) are essential for reactive oxygen species (ROS) homeostasis in responses of plants to environment changes. We previously identified several drought-responsive CC-type GRXs in cassava, an important tropical crop. However, how CC-type GRX regulates ROS homeostasis of cassava under drought stress remained largely unknown. Here, we report that a drought-responsive CC-type GRX, namely MeGRXC3, was associated with activity of catalase in the leaves of 100 cultivars (or unique unnamed genotypes) of cassava under drought stress. MeGRXC3 negatively regulated drought tolerance by modulating drought- and abscisic acid-induced stomatal closure in transgenic cassava. It antagonistically regulated hydrogen peroxide (H2 O2 ) accumulation in epidermal cells and guard cells. Moreover, MeGRXC3 interacted with two catalases of cassava, MeCAT1 and MeCAT2, and regulated their activity in vivo. Additionally, MeGRXC3 interacts with a cassava TGA transcription factor, MeTGA2, in the nucleus, and regulates the expression of MeCAT7 through a MeTGA2-MeMYB63 pathway. Overall, we demonstrated the roles of MeGRXC3 in regulating activity of catalase at both transcriptional and post-translational levels, therefore involving in ROS homeostasis and stomatal movement in responses of cassava to drought stress. Our study provides the first insights into how MeGRXC3 may be used in molecular breeding of cassava crops.

A Transformation and Genome Editing System for Cassava Cultivar SC8.

Cassava starch is a widely used raw material for industrial production. South Chinese cassava cultivar 8 (Manihot esculenta Crantz cv. SC8) is one of the main locally planted cultivars. In this study, an efficient transformation system for cassava SC8 mediated with Agrobacterium strain LBA4404 was presented for the first time. Cassava friable embryogenic calli (FECs) were transformed through the binary vector pCAMBIA1304 harboring GUS- and GFP-fused genes driven by the CaMV35S promoter. The transformation efficiency was increased in the conditions of Agrobacterium strain cell infection density (OD600 = 0.65), 250 µM acetosyringone induction, and agro-cultivation with wet FECs for 3 days in dark. Based on the optimized transformation protocol, approximately 120-140 independent transgenic lines per mL settled cell volume (SCV) of FECs were created by gene transformation in approximately 5 months, and 45.83% homozygous mono-allelic mutations of the MePDS gene with a YAO promoter-driven CRISPR/Cas9 system were generated. This study will open a more functional avenue for the genetic improvement of cassava SC8.

MeNINV1: An Alkaline/Neutral Invertase Gene of Manihot esculenta, Enhanced Sucrose Catabolism and Promoted Plant Vegetative Growth in Transgenic Arabidopsis.

Alkaline/neutral invertase (A/N-INV) is an invertase that irreversibly decomposes sucrose into fructose as well as glucose and plays a role in plant growth and development, starch synthesis, abiotic stress, and other plant-life activities. Cassava is an economically important starch crop in tropical regions. During the development of cassava tuber roots, A/N-INV activity is relatively high, which indicates that it may participate in sucrose metabolism and starch synthesis. In this study, MeNINV1 was confirmed to function as invertase to catalyze sucrose decomposition in yeast. The optimal enzymatic properties of MeNINV1 were a pH of 6.5, a reaction temperature of 40 °C, and sucrose as its specific catalytic substrate. VB6, Zn2+, and Pb2+ at low concentrations as well as EDTA, DTT, Tris, Mg2+, and fructose inhibited A/N-INV enzymic activity. In cassava, the MeNINV1 gene was mainly expressed in the fibrous roots and the tuber root phloem, and its expression decreased as the tuber root grew. MeNINV1 was confirmed to localize in chloroplasts. In Arabidopsis, MeNINV1-overexpressing Arabidopsis had higher A/N-INV activity, and the increased glucose, fructose, and starch content in the leaves promoted plant growth and delayed flowering time but did not change its resistance to abiotic stress. Our results provide new insights into the biological function of MeNINV1.

Transcriptomic profiling suggests candidate molecular responses to waterlogging in cassava.

Owing to climate change impacts, waterlogging is a serious abiotic stress that affects crops, resulting in stunted growth and loss of productivity. Cassava (Manihot esculenta Grantz) is usually grown in areas that experience high amounts of rainfall; however, little research has been done on the waterlogging tolerance mechanism of this species. Therefore, we investigated the physiological responses of cassava plants to waterlogging stress and analyzed global gene transcription responses in the leaves and roots of waterlogged cassava plants. The results showed that waterlogging stress significantly decreased the leaf chlorophyll content, caused premature senescence, and increased the activities of superoxide dismutase (SOD), catalase (CAT) and peroxidase (POD) in the leaves and roots. In total, 2538 differentially expressed genes (DEGs) were detected in the leaves and 13364 in the roots, with 1523 genes shared between the two tissues. Comparative analysis revealed that the DEGs were related mainly to photosynthesis, amino metabolism, RNA transport and degradation. We also summarized the functions of the pathways that respond to waterlogging and are involved in photosynthesis, glycolysis and galactose metabolism. Additionally, many transcription factors (TFs), such as MYBs, AP2/ERFs, WRKYs and NACs, were identified, suggesting that they potentially function in the waterlogging response in cassava. The expression of 12 randomly selected genes evaluated via both quantitative real-time PCR (qRT-PCR) and RNA sequencing (RNA-seq) was highly correlated (R2 = 0.9077), validating the reliability of the RNA-seq results. The potential waterlogging stress-related transcripts identified in this study are representatives of candidate genes and molecular resources for further understanding the molecular mechanisms underlying the waterlogging response in cassava.

Genome-wide identification of cassava MeRboh genes and functional analysis in Arabidopsis.

Plant respiratory burst oxidase homolog (Rboh) gene family encodes NADPH oxidases, and plays important roles in the production of reactive oxygen species (ROS), plant signaling, growth and stress responses. Cassava is an important starchy crops in tropical region. Environmental stresses, such as drought, pathogen, have caused great yield loss. The mechanisms of stress response are little known in MeRBOH family of cassava. Investigation of Rboh genes response to disease may provide a clue for clarification the disease resistance mechanisms. In this study, eight MeRboh genes were identified from the cassava genome. Comparisons of gene structure, protein motifs, and a phylogenetic tree showed conservation of Rboh gene families in cassava, Arabidopsis and rice. Transcript levels of most MeRboh genes increased following treatment with a pathogen, Xanthomonas axonopodis pv. manihotis, or with phytohormones salicylic acid or jasmonic acid. Analysis of cis-acting elements also indicated that MeRboh genes could response to light, hormone, abiotic and biotic stress. Prediction of miRNA target and post-translation modification sites of MeRboh suggested possible regulations of miRNA and protein phosphorylation; and transient expression of MeRboh in cassava protoplasts confirmed their localization on plasma membrane. Expression of MeRbohB, MeRbohF partially complemented PAMP responses in Arabidopsis rboh mutants, including the expression of PTI marker FRK1, ROS production, peroxide accumulation and callose deposition. It suggesting that MeRbohB and MeRbohF may participate in the PTI pathway and contributed to ROS production triggered by pathogens. Moreover, overexpression of MeRbohB and MeRbohF enhanced the resistance of Arabidopsis against Pseudomonas syringae pv. tomato DC3000. Together, these results suggest the evolutionary conservation of MeRboh gene family and their important role in the immune response and in regulating the plant disease resistance, providing a foundation for revealing molecular mechanisms of cassava disease resistance.

Differential expression of the Hsf family in cassava under biotic and abiotic stresses.

Heat shock transcription factors (Hsfs) are important regulators of biotic and abiotic stress responses in plants. Currently, the Hsf gene family is not well understood in cassava, an important tropical crop. In the present study, 32 MeHsf genes were identified from the cassava genome database, which were divided into three types based on functional domain and motif distribution analyses. Analysis of the differential expression of the genes belonging to the Hsf family in cassava was carried out based on published cassava transcriptome data from tissues/organs (leaf blade, leaf midvein, lateral buds, organized embryogenic structures, friable embryogenic callus, fibrous roots, storage roots, stem, petiole, shoot apical meristem, and root apical meristem) under abiotic stress (cold, drought) or biotic stress (mealybugs. cassava brown streak disease, cassava bacterial blight). The results show the expression diversity of cassava Hsfs genes in various tissues/organs. The transcript levels of MeHsfB3a, MeHsfA6a, MeHsfA2a, and MeHsfA9b were upregulated by abiotic and biotic stresses, such as cold, drought, cassava bacterial blight, cassava brown streak disease, and mealybugs, indicating their potential roles in mediating the response of cassava plants to environment stresses. Further interaction network and co-expression analyses suggests that Hsf genes may interact with Hsp70 family members to resist environmental stresses in cassava. These results provide valuable information for future studies of the functional characterization of the MeHsf gene family.

Identification and expression analysis of MinD gene involved in plastid division in cassava.

Cassava is a tropical crop known for its starchy root and excellent properties. Considering that starch biosynthesis in the amyloplast is affected by its division, it appears conceivable that the regulation of plastid division plays an important role in starch accumulation. As a member of the Min system genes, MinD participated in the spatial regulation of the position of the plastid division site.In our studies, sequence analysis and phylogenetic analysis showed that MeMinD has been highly conserved during the evolutionary process. Subcellular localisation indicated that MeMinD carries a chloroplast transit peptide and was localised in the chloroplast. Overexpression of MeMinD resulted in division site misplacement and filamentous formation in E. coli, indicating that MeMinD protein was functional across species. MeMinD exhibited different spatial and temporal expression patterns which was highly expressed in the source compared to that in the sink organ.

Overexpression of MinE gene affects the plastid division in cassava.

The MinE protein plays an important role in plastid division. In this study, the MinE gene was isolated from the cassava (Manihot esculenta Crantz) genome. We isolated high quality and quantity protoplasts and succeed in performing the transient expression of the GFP-fused Manihot esculenta MinE (MeMinE) protein in cassava mesophyll protoplasts. The transient expression of MeMinE-GFP in cassava protoplasts showed that the MeMinE protein was located in the chloroplast. Due to the abnormal division of chloroplasts, overexpression of MeMinE proteins in cassava mesophyll protoplasts could result in fewer and smaller chloroplasts. Overexpression of MeMinE proteins also showed abnormal cell division characteristics and minicell occurrence in Escherichia coli caused by aberrant septation events in the cell poles.

Isolation and Characterization of Ftsz Genes in Cassava.

The filamenting temperature-sensitive Z proteins (FtsZs) play an important role in plastid division. In this study, three FtsZ genes were isolated from the cassava genome, and named MeFtsZ1, MeFtsZ2-1, and MeFtsZ2-2, respectively. Based on phylogeny, the MeFtsZs were classified into two groups (FtsZ1 and FtsZ2). MeFtsZ1 with a putative signal peptide at N-terminal, has six exons, and is classed to FtsZ1 clade. MeFtsZ2-1 and MeFtsZ2-2 without a putative signal peptide, have seven exons, and are classed to FtsZ2 clade. Subcellular localization found that all the three MeFtsZs could locate in chloroplasts and form a ring in chloroplastids. Structure analysis found that all MeFtsZ proteins contain a conserved guanosine triphosphatase (GTPase) domain in favor of generate contractile force for cassava plastid division. The expression profiles of MeFtsZ genes by quantitative reverse transcription-PCR (qRT-PCR) analysis in photosynthetic and non-photosynthetic tissues found that all of the MeFtsZ genes had higher expression levels in photosynthetic tissues, especially in younger leaves, and lower expression levels in the non-photosynthetic tissues. During cassava storage root development, the expressions of MeFtsZ2-1 and MeFtsZ2-2 were comparatively higher than MeFtsZ1. The transformed Arabidopsis of MeFtsZ2-1 and MeFtsZ2-2 contained abnormally shape, fewer number, and larger volume chloroplasts. Phytohormones were involved in regulating the expressions of MeFtsZ genes. Therefore, we deduced that all of the MeFtsZs play an important role in chloroplast division, and that MeFtsZ2 (2-1, 2-2) might be involved in amyloplast division and regulated by phytohormones during cassava storage root development.

Identification, Expression, and Functional Analysis of the Fructokinase Gene Family in Cassava.

Fructokinase (FRK) proteins play important roles in catalyzing fructose phosphorylation and participate in the carbohydrate metabolism of storage organs in plants. To investigate the roles of FRKs in cassava tuber root development, seven FRK genes (MeFRK1-7) were identified, and MeFRK1-6 were isolated. Phylogenetic analysis revealed that the MeFRK family genes can be divided into α (MeFRK1, 2, 6, 7) and ß (MeFRK3, 4, 5) groups. All the MeFRK proteins have typical conserved regions and substrate binding residues similar to those of the FRKs. The overall predicted three-dimensional structures of MeFRK1-6 were similar, folding into a catalytic domain and a ß-sheet ''lid" region, forming a substrate binding cleft, which contains many residues involved in the binding to fructose. The gene and the predicted three-dimensional structures of MeFRK3 and MeFRK4 were the most similar. MeFRK1-6 displayed different expression patterns across different tissues, including leaves, stems, tuber roots, flowers, and fruits. In tuber roots, the expressions of MeFRK3 and MeFRK4 were much higher compared to those of the other genes. Notably, the expression of MeFRK3 and MeFRK4 as well as the enzymatic activity of FRK were higher at the initial and early expanding tuber stages and were lower at the later expanding and mature tuber stages. The FRK activity of MeFRK3 and MeFRK4 was identified by the functional complementation of triple mutant yeast cells that were unable to phosphorylate either glucose or fructose. The gene expression and enzymatic activity of MeFRK3 and MeFRK4 suggest that they might be the main enzymes in fructose phosphorylation for regulating the formation of tuber roots and starch accumulation at the tuber root initial and expanding stages.

Structure, Expression, and Functional Analysis of the Hexokinase Gene Family in Cassava.

Hexokinase (HXK) proteins play important roles in catalyzing hexose phosphorylation and sugar sensing and signaling. To investigate the roles of HXKs in cassava tuber root development, seven HXK genes (MeHXK1-7) were isolated and analyzed. A phylogenetic analysis revealed that the MeHXK family can be divided into five subfamilies of plant HXKs. MeHXKs were clearly divided into type A (MeHXK1) and type B (MeHXK2-7) based on their N-terminal sequences. MeHXK1-5 all had typical conserved regions and similar protein structures to the HXKs of other plants; while MeHXK6-7 lacked some of the conserved regions. An expression analysis of the MeHXK genes in cassava organs or tissues demonstrated that MeHXK2 is the dominant HXK in all the examined tissues (leaves, stems, fruits, tuber phloems, and tuber xylems). Notably, the expression of MeHXK2 and the enzymatic activity of HXK were higher at the initial and expanding tuber stages, and lower at the mature tuber stage. Furthermore, the HXK activity of MeHXK2 was identified by functional complementation of the HXK-deficient yeast strain YSH7.4-3C (hxk1, hxk2, glk1). The gene expression and enzymatic activity of MeHXK2 suggest that it might be the main enzyme for hexose phosphorylation during cassava tuber root development, which is involved in sucrose metabolism to regulate the accumulation of starch.

Effect of HbDHN1 and HbDHN2 Genes on Abiotic Stress Responses in Arabidopsis.

Dehydrin is a type of late embryogenesis abundant (LEA) protein. The dehydrin genes, HbDHN1 and HbDHN2, in Hevea brasiliensis were previously found to be induced at the wounding site of epicormic shoots, with local tissue dehydration identified as the key signal for laticifer differentiation. However, the exact role of the HbDHNs remains unknown. In this study, HbDHN1 and HbDHN2 expression was examined under multiple abiotic stresses; namely, cold, salt, drought, wounding, abscisic acid (ABA), ethylene (ET), and jasmonic acid (JA) treatment. Although, both HbDHNs were defined as SK2-type dehydrin, they showed different cellular localizations. Overexpression of the HbDHNs in Arabidopsis thaliana further revealed a significant increase in tolerance to salt, drought and osmotic stresses. Increased accumulation of proline and a reduction in electrolyte leakage were also observed under salt and drought stress, and a higher water content was indicated under osmotic stress. The transgenic plants also showed higher activity levels of ascorbate peroxidase (APX), superoxide dismutase (SOD) and catalase, and accumulated less hydrogen peroxide (H2O2) and superoxide ([Formula: see text]). Given that reactive oxygen species (ROS) are thought to be a key signal for laticifer differentiation, these findings suggest that HbDHNs act as ROS scavengers, directly or indirectly affecting laticifer differentiation. Both HbDHNs therefore influence physiological processes, improving plant tolerance to multiple abiotic stresses.

Cloning, 3D modeling and expression analysis of three vacuolar invertase genes from cassava (Manihot Esculenta Crantz).

Vacuolar invertase is one of the key enzymes in sucrose metabolism that irreversibly catalyzes the hydrolysis of sucrose to glucose and fructose in plants. In this research, three vacuolar invertase genes, named MeVINV1-3, and with 653, 660 and 639 amino acids, respectively, were cloned from cassava. The motifs of NDPNG (ß-fructosidase motif), RDP and WECVD, which are conserved and essential for catalytic activity in the vacuolar invertase family, were found in MeVINV1 and MeVINV2. Meanwhile, in MeVINV3, instead of NDPNG we found the motif NGPDG, in which the three amino acids GPD are different from those in other vacuolar invertases (DPN) that might result in MeVINV3 being an inactivated protein. The N-terminal leader sequence of MeVINVs contains a signal anchor, which is associated with the sorting of vacuolar invertase to vacuole. The overall predicted 3D structure of the MeVINVs consists of a five bladed ß-propeller module at N-terminus domain, and forms a ß-sandwich module at the C-terminus domain. The active site of the protein is situated in the ß-propeller module. MeVINVs are classified in two subfamilies, α and ß groups, in which α group members of MeVINV1 and 2 are highly expressed in reproductive organs and tuber roots (considered as sink organs), while ß group members of MeVINV3 are highly expressed in leaves (source organs). All MeVINVs are highly expressed in leaves, while only MeVINV1 and 2 are highly expressed in tubers at cassava tuber maturity stage. Thus, MeVINV1 and 2 play an important role in sucrose unloading and starch accumulation, as well in buffering the pools of sucrose, hexoses and sugar phosphates in leaves, specifically at later stages of plant development.

Genome-wide identification, 3D modeling, expression and enzymatic activity analysis of cell wall invertase gene family from cassava (Manihot esculenta Crantz).

The cell wall invertases play a crucial role on the sucrose metabolism in plant source and sink organs. In this research, six cell wall invertase genes (MeCWINV1-6) were cloned from cassava. All the MeCWINVs contain a putative signal peptide with a predicted extracellular location. The overall predicted structures of the MeCWINV1-6 are similar to AtcwINV1. Their N-terminus domain forms a ß-propeller module and three conserved sequence domains (NDPNG, RDP and WECP(V)D), in which the catalytic residues are situated in these domains; while the C-terminus domain consists of a ß-sandwich module. The predicted structure of Pro residue from the WECPD (MeCWINV1, 2, 5, and 6), and Val residue from the WECVD (MeCWINV3 and 4) are different. The activity of MeCWINV1 and 3 were higher than other MeCWINVs in leaves and tubers, which suggested that sucrose was mainly catalyzed by the MeCWINV1 and 3 in the apoplastic space of cassava source and sink organs. The transcriptional levels of all the MeCWINVs and their enzymatic activity were lower in tubers than in leaves at all the stages during the cassava tuber development. It suggested that the major role of the MeCWINVs was on the regulation of carbon exportation from source leaves, and the ratio of sucrose to hexose in the apoplasts; the role of these enzymes on the sucrose unloading to tuber was weaker.