A high-affinity potassium transporter (MeHKT1) from cassava (Manihot esculenta) negatively regulates the response of transgenic Arabidopsis to salt stress

BACKGROUND

High-affinity potassium transporters (HKTs) are crucial in facilitating potassium uptake by plants. Many types of HKTs confer salt tolerance to plants through regulating K+ and Na+ homeostasis under salinity stress. However, their specific functions in cassava (Manihot esculenta) remain unclear.

RESULTS

Herein, an HKT gene (MeHKT1) was cloned from cassava, and its expression is triggered by exposure to salt stress. The expression of a plasma membrane-bound protein functions as transporter to rescue a low potassium (K+) sensitivity of yeast mutant strain, but the complementation of MeHKT1 is inhibited by NaCl treatment. Under low K+ stress, transgenic Arabidopsis with MeHKT1 exhibits improved growth due to increasing shoot K+ content. In contrast, transgenic Arabidopsis accumulates more Na+ under salt stress than wild-type (WT) plants. Nevertheless, the differences in K+ content between transgenic and WT plants are not significant. Additionally, Arabidopsis expressing MeHKT1 displayed a stronger salt-sensitive phenotype.

CONCLUSION

These results suggest that under low K+ condition, MeHKT1 functions as a potassium transporter. In contrast, MeHKT1 mainly transports Na+ into cells under salt stress condition and negatively regulates the response of transgenic Arabidopsis to salt stress. Our results provide a reference for further research on the function of MeHKT1, and provide a basis for further application of MeHKT1 in cassava by molecular biological means.

Evaluation of potential increase in photosynthetic efficiency of cassava (Manihot esculenta Crantz) plants exposed to elevated carbon dioxide

Cassava (Manihot esculenta Crantz), an important tropical crop, is affected by extreme climatic events, including rising CO2 levels. We evaluated the short-term effect of elevated CO2 concentration (ECO2 ) (600, 800 and 1000ppm) on the photosynthetic efficiency of 14 cassava genotypes. ECO2 significantly altered gaseous exchange parameters (net photosynthetic rate (P n ), stomatal conductance (g s ), intercellular CO2 (C i ) and transpiration (E )) in cassava leaves. There were significant but varying interactive effects between ECO2 and varieties on these physiological characteristics. ECO2 at 600 and 800ppm increased the P n rate in the range of 13-24% in comparison to 400ppm (ambient CO2 ), followed by acclimation at the highest concentration of 1000ppm. A similar trend was observed in g s and E . Conversely, C i increased significantly and linearly across increasing CO2 concentration. Along with C i , a steady increase in water use efficiency [WUEintrinsic (P n /g s ) and WUEinstantaneous (P n /E )] across various CO2 concentrations corresponded with the central role of restricted stomatal activity, a common response under ECO2 . Furthermore, P n had a significant quadratic relationship with the ECO2 (R 2 =0.489) and a significant and linear relationship with C i (R 2 =0.227). Relative humidity and vapour pressure deficit during the time of measurements remained at 70-85% and ~0.9-1.31kPa, respectively, at 26±2°C leaf temperature. Notably, not a single variety exhibited constant performance for any of the parameters across CO2 concentrations. Our results indicate that the potential photosynthesis can be increased up to 800ppm cassava varieties with high sink capacity can be cultivated under protected cultivation to attain higher productivity.

Exogenous methyl jasmonate induced cassava defense response and enhanced resistance to Tetranychus urticae

Exogenous application of methyl jasmonate (MeJA) could activate plant defense response against the two-spotted spider mite (TSSM), Tetranychus urticae Koch, in different plants. However, whether MeJA can also serve as an elicitor in cassava (Manihot esculenta Crantz) remains unknown. In this study, induced defense responses were investigated in TSSM-resistant cassava variety C1115 and TSSM-susceptible cassava variety KU50 when applied with MeJA. The performance of TSSM feeding on cassava plants that were pre-treated with various concentrations of MeJA was first evaluated. Subsequently, the activities of antioxidative enzymes (superoxide dismutase and catalase), detoxification enzymes (glutathione S-transferase, cytochrome P450 and carboxylesterase) and digestive enzymes (protease, amylase and invertase) in TSSM were analyzed at days 1, 2, 4 and 8 post-feeding. The results showed that MeJA treatment can induce cassava defense responses to TSSM in terms of reducing egg production and adult longevity as well as slowing development and prolonging the egg stage. Noticeably, C1115 exhibited stronger inhibition of TSSM development and reproduction than KU50. In addition, the activities of all the tested enzymes were induced in both C1115 and KU50, the most in C1115. We conclude that exogenous methyl jasmonate can induce cassava defense responses and enhance resistance to TSSM.

Image-based phenotyping of cassava roots for diversity studies and carotenoids prediction

Phenotyping to quantify the total carotenoids content (TCC) is sensitive, time-consuming, tedious, and costly. The development of high-throughput phenotyping tools is essential for screening hundreds of cassava genotypes in a short period of time in the biofortification program. This study aimed to (i) use digital images to extract information on the pulp color of cassava roots and estimate correlations with TCC, and (ii) select predictive models for TCC using colorimetric indices. Red, green and blue images were captured in root samples from 228 biofortified genotypes and the difference in color was analyzed using L*, a*, b*, hue and chroma indices from the International Commission on Illumination (CIELAB) color system and lightness. Colorimetric data were used for principal component analysis (PCA), correlation and for developing prediction models for TCC based on regression and machine learning. A high positive correlation between TCC and the variables b* (r = 0.90) and chroma (r = 0.89) was identified, while the other correlations were median and negative, and the L* parameter did not present a significant correlation with TCC. In general, the accuracy of most prediction models (with all variables and only the most important ones) was high (R2 ranging from 0.81 to 0.94). However, the artificial neural network prediction model presented the best predictive ability (R2 = 0.94), associated with the smallest error in the TCC estimates (root-mean-square error of 0.24). The structure of the studied population revealed five groups and high genetic variability based on PCA regarding colorimetric indices and TCC. Our results demonstrated that the use of data obtained from digital image analysis is an economical, fast, and effective alternative for the development of TCC phenotyping tools in cassava roots with high predictive ability.

Correlation Between Relative Gene Expression Patterns of Two Flowering locus T (MeFT1 and MeFT2) in Cassava Leaf and Flowering Traits Under Different Flowering Induction Conditions

<b>Background and Objective:</b> <i>Flowering locus T</i> (<i>FT</i>) genes are involved in the flower induction mechanism in plants as florigen signals. The objective of this study was to study the relationship between the expression of <i>Flowering locus T</i> genes (<i>MeFTs</i>) in cassava and flowering traits under the different flowering induction conditions. <b>Materials and Methods:</b> The experimental design for flowering induction was RCBD for 4 replications. There were 5 treatment factors which were control, red light set from 5 pm to 7 am, 0.5 mM 6-benzyladenine (BA) with 2 mM silver thio-sulfate (STS), paclobutrazol for 6 g/plant and potassium chlorate (KCIO<sub>3</sub>) for 250 g/plant. The number of plants with flower bunches and the average number of bunches per plant in two cassava varieties were collected each month from 5-9 months after planting (MAP). The leaf samples were collected from HB80 and R9 varieties at 5-7 MAP for RNA extraction to study <i>MeFT1</i> and <i>MeFT2</i> expression. <b>Results:</b> The results show that <i>MeFT1</i> expression level positively correlated with flowering traits in the same month, while <i>MeFT2</i> expression level positively correlated with flowering traits in the following months. <b>Conclusion:</b> Therefore, expression of <i>MeFT2</i> can be used for the prediction of cassava flowering in the following month which will assist the breeder for the crossing management.

Character changes and Transcriptomic analysis of a cassava sexual Tetraploid

BACKGROUND

Cassava (Manihot esculenta Crantz) is an important food crop known for its high starch content. Polyploid breeding is effective in its genetic improvement, and use of 2n gametes in sexual polyploid breeding is one of the potential methods for cassava breeding and improvement. In our study, the cassava sexual tetraploid (ST), which carries numerous valuable traits, was successfully generated by hybridizing 2n female gametes SC5 (♀) and 2n male gametes SC10 (♂). However, the molecular mechanisms remain unclear. To understand these underlying molecular mechanisms behind the phenotypic alterations and heterosis in ST plants, we investigated the differences in gene expression between polyploids and diploids by determining the transcriptomes of the ST plant and its parents during the tuber root enlargement period. We also compared the characters and transcriptomes of the ST plant with its parents.

RESULTS

The ST plant was superior in plant height, stem diameter, leaf area, petiole length, plant weight, and root weight than the parent plants, except the leaf number, which was lower. The number of starch granules was higher in the roots of ST plants than those in the parent plants after five months (tuber root enlargement period), which could be due to a higher leaf net photosynthetic rate leading to early filling of starch granules. Based on transcriptome analysis, we identified 2934 and 3171 differentially expressed genes (DEGs) in the ST plant as compared to its female and male parents, respectively. Pathway enrichment analyses revealed that flavonoid biosynthesis and glycolysis/gluconeogenesis were significantly enriched in the ST plants, which might contribute to the colors of petiole (purple-red), root epidermis (dark brown), and tuber starch accumulation, respectively.

CONCLUSIONS

After sexual polyploidization, the phenotype of ST has changed significantly in comparison to their diploid parents, mainly manifest as enlarged biomass, yield, early starch filling, deep colored petiole and root epidermis. The tetraploid plants were also mature early due to early starch grain filling. Owing to enriched flavonoid biosynthesis and glycolysis/gluconeogenesis, they are possibly resistant to adversity stresses and provide better yield, respectively.

Identification of cassava alternative splicing-related genes and functional characterization of MeSCL30 involvement in drought stress

Alternative splicing (AS) is an important post-transcriptional regulation strategy that can increase the proteome diversity and regulate mRNA level in eukaryote. Multi-exon genes can be alternative spliced to generate two or more transcripts, thereby increasing the adaptation to the external stress conditions in planta. However, AS-related proteins were less explored in cassava which is an important staple crop in the tropical area. A total of 365 genes encoding AS-related proteins were identified and renamed in the cassava genome, and the transcriptional and splicing changes of 15 randomly selected genes were systematically investigated in the tissues under diverse abiotic stress conditions. 13 out of 15 genes undergo AS in the tissues and under diverse environmental stress condition. Importantly, the greatest changes of splicing patterns were found in the leaf or in response to temperature stress, indicating that AS-related proteins had their tissue-specific regulation patterns and might be participated in the plant adaptation to temperature stress. We then found that overexpression of MeSCL30 in Arabidopsis enhanced the tolerance to drought stress through maintaining reactive oxygen species (ROS) homeostasis and increasing the expression of drought-responsive genes. Therefore, these findings refined the AS-related protein-coding genes and provided novel insights for manipulation of AS-related genes in order to enhance the resistance to environmental stress in plant.

The metabotyping of an East African cassava diversity panel: A core collection for developing biotic stress tolerance in cassava

Cassava will have a vital role to play, if food security is to be achieved in Sub-Saharan Africa, especially Central and East Africa. The whitefly Bemisia tabaci poses a major threat to cassava production by small holder farmers in part due to their role as a vector of cassava mosaic begomoviruses (CMBs) and cassava brown streak ipomoviruses (CBSIs). In the present study untargeted metabolomics has been used as a tool to assess natural variation, similarities and attempts to identify trait differentiators among an East African cassava diversity panel that displayed tolerance/resistance to the effects of Bemisia tabaci infestation. The metabolome captured, was represented by 1529 unique chemical features per accession. Principal component analysis (PCA) identified a 23% variation across the panel, with geographical origin/adaption the most influential classification factors. Separation based on resistance and susceptible traits to Bemisia tabaci could also be observed within the data and was corroborated by genotyping data. Thus the metabolomics pipeline represented an effective metabotyping approach. Agglomerative Hierarchical Clustering Analysis (HCA) of both the metabolomics and genotyping data was performed and revealed a high level of similarity between accessions. Specific differentiating features/metabolites were identified, including those potentially conferring vigour to whitefly tolerance on a constitutive manner. The implications of using these cassava varieties as parental breeding material and the future potential of incorporating more exotic donor material is discussed.

Genetic variation and evolutionary history of a mycorrhizal fungus regulate the currency of exchange in symbiosis with the food security crop cassava

Most land plants form symbioses with arbuscular mycorrhizal fungi (AMF). Diversity of AMF increases plant community productivity and plant diversity. For decades, it was known that plants trade carbohydrates for phosphate with their fungal symbionts. However, recent studies show that plant-derived lipids probably represent the most essential currency of exchange. Understanding the regulation of plant genes involved in the currency of exchange is crucial to understanding stability of this mutualism. Plants encounter many different AMF genotypes that vary greatly in the benefit they confer to plants. Yet the role that fungal genetic variation plays in the regulation of this currency has not received much attention. We used a high-resolution phylogeny of one AMF species (Rhizophagus irregularis) to show that fungal genetic variation drives the regulation of the plant fatty acid pathway in cassava (Manihot esculenta); a pathway regulating one of the essential currencies of trade in the symbiosis. The regulation of this pathway was explained by clearly defined patterns of fungal genome-wide variation representing the precise fungal evolutionary history. This represents the first demonstrated link between the genetics of AMF and reprogramming of an essential plant pathway regulating the currency of exchange in the symbiosis. The transcription factor RAM1 was also revealed as the dominant gene in the fatty acid plant gene co-expression network. Our study highlights the crucial role of variation in fungal genomes in the trade of resources in this important symbiosis and also opens the door to discovering characteristics of AMF genomes responsible for interactions between AMF and cassava that will lead to optimal cassava growth.

Improving root characterisation for genomic prediction in cassava

Cassava is cultivated due to its drought tolerance and high carbohydrate-containing storage roots. The lack of uniformity and irregular shape of storage roots poses constraints on harvesting and post-harvest processing. Here, we phenotyped the Genetic gain and offspring (C1) populations from the International Institute of Tropical Agriculture (IITA) breeding program using image analysis of storage root photographs taken in the field. In the genome-wide association analysis (GWAS), we detected for most shape and size-related traits, QTL on chromosomes 1 and 12. In a previous study, we found the QTL on chromosome 12 to be associated with cassava mosaic disease (CMD) resistance. Because the root uniformity is important for breeding, we calculated the standard deviation (SD) of individual root measurements per clone. With SD measurements we identified new significant QTL for Perimeter, Feret and Aspect Ratio on chromosomes 6, 9 and 16. Predictive accuracies of root size and shape image-extracted traits were mostly higher than yield trait prediction accuracies. This study aimed to evaluate the feasibility of the image phenotyping protocol and assess GWAS and genomic prediction for size and shape image-extracted traits. The methodology described and the results are promising and open up the opportunity to apply high-throughput methods in cassava.

The chaperone MeHSP90 recruits MeWRKY20 and MeCatalase1 to regulate drought stress resistance in cassava

The 90 kDa heat shock protein (HSP90) is widely involved in various developmental processes and stress responses in plants. However, the molecular chaperone HSP90-constructed protein complex and its function in cassava remain elusive. In this study, we report that HSP90 is essential for drought stress resistance in cassava by regulating abscisic acid (ABA) and hydrogen peroxide (H₂ O₂ ) using two specific protein inhibitors of HSP90 (geldanamycin (GDA) and radicicol (RAD)). Among 10 MeHSP90s, the transcript of MeHSP90.9 is largely induced during drought stress. Further investigation identifies MeWRKY20 and MeCatalase1 as MeHSP90.9-interacting proteins. MeHSP90.9-, MeWRKY20-, or MeCatalase1-silenced plants through virus-induced gene silencing display drought sensitivity in cassava, indicating that they are important to drought stress response. MeHSP90.9 can promote the direct transcriptional activation of MeWRKY20 on the W-box element of MeNCED5 promoter, encoding a key enzyme in ABA biosynthesis. Moreover, MeHSP90.9 positively regulates the activity of MeCatalase1, and MeHSP90.9-silenced cassava leaves accumulate more H₂ O₂ under drought stress. Taken together, we demonstrate that the MeHSP90.9 chaperone complex is a regulator of drought stress resistance in cassava.

Genomic and Transcriptomic Analysis Identified Novel Putative Cassava lncRNAs Involved in Cold and Drought Stress

Long non-coding RNAs (lncRNAs) play important roles in the regulation of complex cellular processes, including transcriptional and post-transcriptional regulation of gene expression relevant for development and stress response, among others. Compared to other important crops, there is limited knowledge of cassava lncRNAs and their roles in abiotic stress adaptation. In this study, we performed a genome-wide study of ncRNAs in cassava, integrating genomics- and transcriptomics-based approaches. In total, 56,840 putative ncRNAs were identified, and approximately half the number were verified using expression data or previously known ncRNAs. Among these were 2229 potential novel lncRNA transcripts with unmatched sequences, 250 of which were differentially expressed in cold or drought conditions, relative to controls. We showed that lncRNAs might be involved in post-transcriptional regulation of stress-induced transcription factors (TFs) such as zinc-finger, WRKY, and nuclear factor Y gene families. These findings deepened our knowledge of cassava lncRNAs and shed light on their stress-responsive roles.

Genome-Wide Identification of Cassava Serine/Arginine-Rich Proteins: Insights into Alternative Splicing of Pre-mRNAs and Response to Abiotic Stress

Serine/arginine-rich (SR) proteins have an essential role in the splicing of pre-messenger RNA (pre-mRNA) in eukaryote. Pre-mRNA with introns can be alternatively spliced to generate multiple transcripts, thereby increasing adaptation to the external stress conditions in planta. However, pre-mRNA of SR proteins can also be alternatively spliced in different plant tissues and in response to diverse stress treatments, indicating that SR proteins might be involved in regulating plant development and adaptation to environmental changes. We identified and named 18 SR proteins in cassava and systematically studied their splicing and transcriptional changes under tissue-specific and abiotic stress conditions. Fifteen out of 18 SR genes showed alternative splicing in the tissues. 45 transcripts were found from 18 SR genes under normal conditions, whereas 55 transcripts were identified, and 21 transcripts were alternate spliced in some SR genes under salt stress, suggesting that SR proteins might participate in the plant adaptation to salt stress. We then found that overexpression of MeSR34 in Arabidopsis enhanced the tolerance to salt stress through maintaining reactive oxygen species homeostasis and increasing the expression of calcineurin B-like proteins (CBL)-CBL-interacting protein kinases and osmotic stress-related genes. Therefore, our findings highlight the critical role of cassava SR proteins as regulators of RNA splicing and salt tolerance in planta.

Symplasmic phloem unloading and radial post-phloem transport via vascular rays in tuberous roots of Manihot esculenta

Cassava (Manihot esculenta) is one of the most important staple food crops worldwide. Its starchy tuberous roots supply over 800 million people with carbohydrates. Yet, surprisingly little is known about the processes involved in filling of those vital storage organs. A better understanding of cassava carbohydrate allocation and starch storage is key to improving storage root yield. Here, we studied cassava morphology and phloem sap flow from source to sink using transgenic pAtSUC2::GFP plants, the phloem tracers esculin and 5(6)-carboxyfluorescein diacetate, as well as several staining techniques. We show that cassava performs apoplasmic phloem loading in source leaves and symplasmic unloading into phloem parenchyma cells of tuberous roots. We demonstrate that vascular rays play an important role in radial transport from the phloem to xylem parenchyma cells in tuberous roots. Furthermore, enzymatic and proteomic measurements of storage root tissues confirmed high abundance and activity of enzymes involved in the sucrose synthase-mediated pathway and indicated that starch is stored most efficiently in the outer xylem layers of tuberous roots. Our findings form the basis for biotechnological approaches aimed at improved phloem loading and enhanced carbohydrate allocation and storage in order to increase tuberous root yield of cassava.

The chloroplast proteome response to drought stress in cassava leaves

Cassava is an important tropical crop with strong resistance to drought stress. The chloroplast, the site of photosynthesis, is sensitive to stress, and the drought-response proteins in cassava chloroplasts are worthy of investigation. In this study, cassava leaves were collected for ultra-structure observation from plants subjected to different drought stress conditions. Our results showed that drought stress can promote starch accumulation in cassava chloroplasts. To evaluate changes in chloroplast proteins under different drought conditions, two-dimensional electrophoresis was performed using purified chloroplasts, which resulted in the identification of 26 unique chloroplast proteins responsive to drought stress. These drought-responsive proteins are predominantly related to photosynthesis, carbon and nitrogen metabolism, and amino acid metabolism. Among them, most photosynthesis-related proteins are downregulated, with decreases in photosynthetic parameters upon drought stress. Several proteins associated with carbon and nitrogen metabolism, including rubisco and carbonic anhydrase, were upregulated, which might promote drought tolerance in cassava by enhancing the carbohydrate conversion efficiency and protecting the plant from oxidative stress. Our proteomic data not only provide insight into the complement of proteins in cassava chloroplasts but also further our overall understanding of drought-responsive proteins in cassava chloroplasts.

Transcriptomic and proteomic response of Manihot esculenta to Tetranychus urticae infestation at different densities

Tetranychus urticae (Acari: Tetranychidae) is an extremely serious cassava (Manihot esculenta) pest. Building a genomic resource to investigate the molecular mechanisms of cassava responses to T. urticae is vital for characterizing cassava resistance to mites. Based on the tolerance of cassava varieties to mite infestation (focusing on mite development rate, fecundity and physiology), cassava variety SC8 was selected to analyze transcriptomic and proteomic changes after 5 days of T. urticae feeding. Transcriptomic analysis revealed 698 and 2140 genes with significant expression changes under low and high mite infestation, respectively. More defense-related genes were found in the enrichment pathways at high mite density than at low density. In addition, iTRAQ-labeled proteomic analysis revealed 191 proteins with significant expression changes under low mite infestation. Differentially expressed genes and proteins were mainly found in the following defense-related pathways: flavonoid biosynthesis, phenylpropanoid biosynthesis, and glutathione metabolism under low-density mite feeding and plant hormone signal transduction and plant-pathogen interaction pathways under high-density mite feeding. The plant hormone signal transduction network, involving ethylene, jasmonic acid, and salicylic acid transduction pathways, was explored in relation to the M. esculenta response to T. urticae. Correlation analysis of the transcriptome and proteome generated a Pearson correlation coefficients of R = 0.2953 (P < 0.01), which might have been due to post-transcriptional or post-translational regulation resulting in many genes being inconsistently expressed at both the transcript and protein levels. In summary, the M. esculenta transcriptome and proteome changed in response to T. urticae, providing insight into the general activation of plant defense pathways in response to mite infestation.

Identification of FT family genes that respond to photoperiod, temperature and genotype in relation to flowering in cassava (Manihot esculenta, Crantz)

Cassava is a starch-storing root crop that is an important source of dietary energy in tropical regions of the world. Genetic improvement of cassava by breeding is hindered by late flowering and sparse flower production in lines that are needed as parents. To advance understanding of regulatory mechanisms in cassava, this work sought to identify and characterize homologs of the FLOWERING LOCUS T (FT) gene. Ten members of the phosphatidylethanolamine-binding protein gene family, to which FT belongs, were obtained from the cassava genome database. Phylogenetic and sequence analysis of these proteins was used to identify two putative FT homologs which had amino acid sequences at key positions in accordance with those predicted for functional FTs. Expression of these ten genes was determined in mature leaves, immature leaves, flower buds, fibrous roots, storage roots and stem. The FT transcripts were expressed in mature leaves, as expected for their possible role in leaf-to-apical meristem signaling. In growth chamber studies, plants flowered earlier in long-day photoperiod than in short-day photoperiod. Expression studies indicated that while MeFT1 was expressed in leaves without a clear-cut photoperiod response, MeFT2 was expressed in a photoperiod-dependent manner, consistent with its involvement in photoperiodic control of flowering. In growth chambers that subjected plants to a range of temperatures from 22 to 34 °C, flowering was delayed by warmer temperatures although MeFT1 and MeFT2 expression declined in only one genotype, indicating other factors regulate this response. The earliest flowering genotype, IBA980002, had high levels of MeFT1 and MeFT2 expression, suggesting that both homologs contribute to earliness of this genotype.

Extensive Post-Transcriptional Regulation Revealed by Transcriptomic and Proteomic Integrative Analysis in Cassava under Drought

Cassava is a major tropical/subtropical food crop and its yield is greatly restrained by drought; however, the mechanism underlying the drought stress remains largely unknown. In this study, totally 1242 and 715 differentially expressed genes (DEGs), together with 237 and 307 differentially expressed proteins (DEPs), were respectively identified in cassava leaves and roots through RNA-seq and iTRAQ techniques. The majority of DEGs and DEPs were exclusively regulated at the mRNA and protein level, respectively, whereas only a few were commonly regulated, indicating the major involvement of post-transcriptional regulation under drought. Subsequently, the functions of these specifically or commonly regulated DEGs and DEPs were analyzed, and the post-transcriptional regulation of genes involved in heat shock protein, secondary metabolism biosynthesis, and hormone biosynthesis was extensively discussed. This is the first report on an integration of transcriptomic and proteomic analysis in cassava, and it provides new insights into the post-transcriptional regulation of cassava drought stress.

Identification and characterization of drought-responsive CC-type glutaredoxins from cassava cultivars reveals their involvement in ABA signalling

BACKGROUND

CC-type glutaredoxins (GRXs) are plant-specific glutaredoxin, play regulatory roles in response of biotic and abiotic stress. However, it is not clear whether the CC-type GRXs are involve in drought response in cassava (Manihot esculenta), an important tropical tuber root crop.

RESULTS

Herein, genome-wide analysis identified 18 CC-type GRXs in the cassava genome, of which six (namely MeGRXC3, C4, C7, C14, C15, and C18) were induced by drought stress in leaves of two cassava cultivars Argentina 7 (Arg7) and South China 124 (SC124). Exogenous abscisic acid (ABA) application induced the expression of all the six CC-type GRXs in leaves of both Arg7 and SC124 plants. Overexpression of MeGRXC15 in Arabidopsis (Col-0) increases tolerance of ABA on the sealed agar plates, but results in drought hypersensitivity in soil-grown plants. The results of microarray assays show that MeGRXC15 overexpression affected the expression of a set of transcription factors which involve in stress response, ABA, and JA/ET signalling pathway. The results of protein interaction analysis show that MeGRXC15 can interact with TGA5 from Arabidopsis and MeTGA074 from cassava.

CONCLUSIONS

CC-type glutaredoxins play regulatory roles in cassava response to drought possibly through ABA signalling pathway.

Comparison of leaf transcriptomes of cassava "Xinxuan 048" diploid and autotetraploid plants

Polyploidy breeding of cassava has been used to improve cassava traits over the past years. We previously reported in vitro induction of tetraploids in the cassava variety "Xinxuan 048" using colchicine. Significant differences in morphology and anatomy were found between the diploid and tetraploid plants. However, very little is known about the transcriptome difference between them. In this study, morphological and physiological characteristics including leaf thickness, plant height, internode length, chlorophyll content, and photosynthetic capacity were measured. Further, we investigated and validated the difference in gene expression patterns between cassava "Xinxuan 048" tetraploid genotype and its diploid plants using RNA sequencing (RNAseq) and quantitative reverse transcription-polymerase chain reaction (qRT-PCR). Significant differences in morphology and physiology were observed during tetraploidization. A comparison revealed that tetraploidy induced very limited changes in the leaf transcriptomes of cassava "Xinxuan 048" diploid and autotetraploid plants. However, the differentially expressed genes (DEGs) between 2× and 4× plants, especially those upregulated in 4× plants, were strongly associated with hormonal and stress responses. Large changes in morphology and physiology between the diploid cassava "Xinxuan 048" and its autotetraploid were not associated with large changes in their leaf transcriptomes. Moreover, the differently expressed genes related to the regulation of gibberellin and brassinosteroids potentially explained why the plant height and internode length of 4× plants became shorter. Collectively, our results suggest that 4× cassava is potentially valuable for breeding strains with improved stress resistance.

Molecular identification of GAPDHs in cassava highlights the antagonism of MeGAPCs and MeATG8s in plant disease resistance against cassava bacterial blight

MeGAPCs were identified as negative regulators of plant disease resistance, and the interaction of MeGAPCs and MeATG8s was highlighted in plant defense response. As an important enzyme of glycolysis metabolic pathway, glyceraldehyde-3-P dehydrogenase (GAPDH) plays important roles in plant development, abiotic stress and immune responses. Cassava (Manihot esculenta) is most important tropical crop and one of the major food crops, however, no information is available about GAPDH gene family in cassava. In this study, 14 MeGAPDHs including 6 cytosol GAPDHs (MeGAPCs) were identified from cassava, and the transcripts of 14 MeGAPDHs in response to Xanthomonas axonopodis pv manihotis (Xam) indicated their possible involvement in immune responses. Further investigation showed that MeGAPCs are negative regulators of disease resistance against Xam. Through transient expression in Nicotiana benthamiana, we found that overexpression of MeGAPCs led to decreased disease resistance against Xam. On the contrary, MeGAPCs-silenced cassava plants through virus-induced gene silencing (VIGS) conferred improved disease resistance. Notably, MeGAPCs physically interacted with autophagy-related protein 8b (MeATG8b) and MeATG8e and inhibited autophagic activity. Moreover, MeATG8b and MeATG8e negatively regulated the activities of NAD-dependent MeGAPDHs, and are involved in MeGAPCs-mediated disease resistance. Taken together, this study highlights the involvement of MeGAPCs in plant disease resistance, through interacting with MeATG8b and MeATG8e.

Species- and density-dependent induction of volatile organic compounds by three mite species in cassava and their role in the attraction of a natural enemy

Upon damage by herbivores, plants induce an array of volatile organic compounds (VOCs) that mediate ecological interactions involving communication with organisms of the second and third trophic levels. VOC-mediated tritrophic interactions have largely been studied in various systems, including cassava (Manihot esculenta), but little is known about the chemical nature of herbivore-induced VOCs in this crop and the response they evoke in natural enemies. Several tetranychid and predatory mites are associated with cassava. Here, VOC emissions from uninfested plants and plants infested with 200 or 400 Mononychellus tanajoa, a specialist herbivore on cassava, and the generalists Tetranychus urticae and T. gloveri were measured. Dual-choice experiments were also conducted to assess the preference of inexperienced (reared on prey-infested bean plants) and experienced (adapted on prey-infested cassava plants) predatory mites, Neoseiulus idaeus (Phytoseiidae), between odors of uninfested plants versus odors of plants infested with M. tanajoa, T. urticae or T. gloveri. Two hundred individuals significantly increased the emissions of (Z)-3-hexen-1-ol, (E)-β-ocimene, β-caryophyllene, alloaromadendrene and (E)-geranyl acetone in T. urticae-infested plants, and (E)-β-ocimene and methyl salicylate (MeSA) in T. gloveri-infested plants. Four hundred individuals significantly increased the emissions of (Z)-3-hexen-1-ol, MeSA, α-pinene and D-limonene in M. tanajoa-infested plants. In addition, T. urticae at this density induced (E)-β-ocimene, D-limonene, (E)-geranyl acetone and six compounds that were not detected in other treatments. Tetranychus gloveri-infested plants induced the emissions of (E)-2-hexenal and D-limonene. Regardless of the infesting species, inexperienced N. idaeus did not discriminate between uninfested or infested plants. Upon experience, they discriminated between the odors of uninfested and T. urticae-damaged plants. Our findings reveal that mite infestations in cassava result in density-dependent and species-specific emission of VOCs, and that N. idaeus relies on associative learning to forage for its prey.

MeWRKY20 and its interacting and activating autophagy-related protein 8 (MeATG8) regulate plant disease resistance in cassava

As a highly conserved mechanism, autophagy is responsible for the transport of cytoplasmic constituents in the vacuoles or lysosomes. Moreover, autophagy is essential for plant development and various stress responses. In this study, 34 MeATGs were systematically identified in cassava, and their transcripts were commonly regulated by Xanthomonas axonopodis pv manihotis (Xam). Through transient expression in Nicotiana benthamiana, the subcellular locations of 4 MeATG8s were revealed. Notably, MeWRKY20 was identified as physical interacting protein of MeATG8a/8f/8h and upstream transcriptional activator of MeATG8a. Through virus-induced gene silencing (VIGS) in cassava, we found that MeATG8-silenced and MeWRKY20-silenced plants resulted in disease sensitive, with less callose depositions and lower autophagic activity. This study may facilitate our understanding of the upstream MeWRKY20 and underlying target as well as interacting proteins of MeATG8s in immune response. Taken together, MeWRKY20 and MeATG8a/8f/8h are essential for disease resistance against bacterial blight by forming various transcriptional modules and interacting complex in cassava.

Prospects for Genomic Selection in Cassava Breeding

Cassava ( Crantz) is a clonally propagated staple food crop in the tropics. Genomic selection (GS) has been implemented at three breeding institutions in Africa to reduce cycle times. Initial studies provided promising estimates of predictive abilities. Here, we expand on previous analyses by assessing the accuracy of seven prediction models for seven traits in three prediction scenarios: cross-validation within populations, cross-population prediction and cross-generation prediction. We also evaluated the impact of increasing the training population (TP) size by phenotyping progenies selected either at random or with a genetic algorithm. Cross-validation results were mostly consistent across programs, with nonadditive models predicting of 10% better on average. Cross-population accuracy was generally low (mean = 0.18) but prediction of cassava mosaic disease increased up to 57% in one Nigerian population when data from another related population were combined. Accuracy across generations was poorer than within-generation accuracy, as expected, but accuracy for dry matter content and mosaic disease severity should be sufficient for rapid-cycling GS. Selection of a prediction model made some difference across generations, but increasing TP size was more important. With a genetic algorithm, selection of one-third of progeny could achieve an accuracy equivalent to phenotyping all progeny. We are in the early stages of GS for this crop but the results are promising for some traits. General guidelines that are emerging are that TPs need to continue to grow but phenotyping can be done on a cleverly selected subset of individuals, reducing the overall phenotyping burden.

A Risk Assessment Framework for Seed Degeneration: Informing an Integrated Seed Health Strategy for Vegetatively Propagated Crops

Pathogen buildup in vegetative planting material, termed seed degeneration, is a major problem in many low-income countries. When smallholder farmers use seed produced on-farm or acquired outside certified programs, it is often infected. We introduce a risk assessment framework for seed degeneration, evaluating the relative performance of individual and combined components of an integrated seed health strategy. The frequency distribution of management performance outcomes was evaluated for models incorporating biological and environmental heterogeneity, with the following results. (1) On-farm seed selection can perform as well as certified seed, if the rate of success in selecting healthy plants for seed production is high; (2) when choosing among within-season management strategies, external inoculum can determine the relative usefulness of 'incidence-altering management' (affecting the proportion of diseased plants/seeds) and 'rate-altering management' (affecting the rate of disease transmission in the field); (3) under severe disease scenarios, where it is difficult to implement management components at high levels of effectiveness, combining management components can be synergistic and keep seed degeneration below a threshold; (4) combining management components can also close the yield gap between average and worst-case scenarios. We also illustrate the potential for expert elicitation to provide parameter estimates when empirical data are unavailable. [Formula: see text] Copyright © 2017 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .