An ordered EST catalogue and gene expression profiles of cassava (Manihot esculenta) at key growth stages

Integrated Metabolomic and Transcriptomic Analyses Reveal Novel Insights of Anthocyanin Biosynthesis on Color Formation in Cassava Tuberous Roots

Yellow roots are of higher nutritional quality and better appearance than white roots in cassava, a crucial tropical and subtropical root crop. In this work, two varieties with yellow and white cassava roots were selected to explore the mechanisms of color formation by using comparative metabolome and transcriptome analyses during seven developmental stages. Compared with the white-rooted cassava, anthocyanins, catechin derivatives, coumarin derivatives, and phenolic acids accumulated at higher levels in yellow-rooted cassava. Anthocyanins were particularly enriched and displayed different accumulation patterns during tuberous root development. This was confirmed by metabolic comparisons between five yellow-rooted and five white-rooted cassava accessions. The integrative metabolomic and transcriptomic analysis further revealed a coordinate regulation of 16 metabolites and 11 co-expression genes participating in anthocyanin biosynthesis, suggesting a vital role of anthocyanin biosynthesis in yellow pigmentation in cassava tuberous roots. In addition, two transcriptional factors, i.e., MeMYB5 and MeMYB42, were also identified to co-express with these anthocyanin biosynthesis genes. These findings expand our knowledge on the role of anthocyanin biosynthesis in cassava root color formation, and offer useful information for the genetic breeding of yellow-rooted cassava in the future.

Highly dynamic, coordinated, and stage-specific profiles are revealed by a multi-omics integrative analysis during tuberous root development in cassava

Cassava (Manihot esculenta) is an important starchy root crop that provides food for millions of people worldwide, but little is known about the regulation of the development of its tuberous root at the multi-omics level. In this study, the transcriptome, proteome, and metabolome were examined in parallel at seven time-points during the development of the tuberous root from the early to late stages of its growth. Overall, highly dynamic and stage-specific changes in the expression of genes/proteins were observed during development. Cell wall and auxin genes, which were regulated exclusively at the transcriptomic level, mainly functioned during the early stages. Starch biosynthesis, which was controlled at both the transcriptomic and proteomic levels, was mainly activated in the early stages and was greatly restricted during the late stages. Two main branches of lignin biosynthesis, coniferyl alcohol and sinapyl alcohol, also functioned during the early stages of development at both the transcriptomic and proteomic levels. Metabolomic analysis further supported the stage-specific roles of particular genes/proteins. Metabolites related to lignin and flavonoid biosynthesis showed high abundance during the early stages, those related to lipids exhibited high abundance at both the early and middle stages, while those related to amino acids were highly accumulated during the late stages. Our findings provide a comprehensive resource for broadening our understanding of tuberous root development and will facilitate future genetic improvement of cassava.

Exploring dynamic protein-protein interactions in cassava through the integrative interactome network

Protein-protein interactions (PPIs) play an essential role in cellular regulatory processes. Despite, in-depth studies to uncover the mystery of PPI-mediated regulations are still lacking. Here, an integrative interactome network (MePPI-Ux) was obtained by incorporating expression data into the improved genome-scale interactome network of cassava (MePPI-U). The MePPI-U, constructed by both interolog- and domain-based approaches, contained 3,638,916 interactions and 24,590 proteins (59% of proteins in the cassava AM560 genome version 6). After incorporating expression data as information of state, the MePPI-U rewired to represent condition-dependent PPIs (MePPI-Ux), enabling us to envisage dynamic PPIs (DPINs) that occur at specific conditions. The MePPI-Ux was exploited to demonstrate timely PPIs of cassava under various conditions, namely drought stress, brown streak virus (CBSV) infection, and starch biosynthesis in leaf/root tissues. MePPI-Ux_drought and MePPI-Ux_CBSV suggested involved PPIs in response to stress. MePPI-Ux_SB,leaf and MePPI-Ux_SB,root suggested the involvement of interactions among transcription factor proteins in modulating how leaf or root starch is synthesized. These findings deepened our knowledge of the regulatory roles of PPIs in cassava and would undeniably assist targeted breeding efforts to improve starch quality and quantity.

Pan- and core- gene association networks: Integrative approaches to understanding biological regulation

The rapid increase in transcriptome data provides an opportunity to access the complex regulatory mechanisms in cellular systems through gene association network (GAN). Nonetheless, GANs derived from single datasets generally allow us to envisage only one side of the regulatory network, even under the particular condition of study. The circumstance is well demonstrated by inconsistent GANs of individual datasets proposed for similar experimental conditions, which always leads to ambiguous interpretation. Here, pan- and core-gene association networks (pan- and core-GANs), analogous to the pan- and core-genome concepts, are proposed to increase the power of inference through the integration of multiple, diverse datasets. The core-GAN represents the consensus associations of genes that were inferred from all individual networks. On the other hand, the pan-GAN represents the extensive gene-gene associations that occurred in each individual network. The pan- and core-GANs prospects were demonstrated based on three time series microarray datasets in leaves of Arabidopsis thaliana grown under diurnal conditions. We showed the overall performance of pan- and core-GANs was more robust to the number of data points in gene expression data compared to the GANs inferred from individual datasets. In addition, the incorporation of multiple data broadened our understanding of the biological regulatory system. While the pan-GAN enabled us to observe the landscape of gene association system, core-GAN highlighted the basic gene-associations in essence of the regulation regulating starch metabolism in leaves of Arabidopsis.

Unlocking conserved and diverged metabolic characteristics in cassava carbon assimilation via comparative genomics approach

Globally, cassava is an important source of starch, which is synthesized through carbon assimilation in cellular metabolism whereby harvested atmospheric carbon is assimilated into macromolecules. Although the carbon assimilation pathway is highly conserved across species, metabolic phenotypes could differ in composition, type, and quantity. To unravel the metabolic complexity and advantage of cassava over other starch crops, in terms of starch production, we investigated the carbon assimilation mechanisms in cassava through genome-based pathway reconstruction and comparative network analysis. First, MeRecon - the carbon assimilation pathway of cassava was reconstructed based upon six plant templates: Arabidopsis, rice, maize, castor bean, potato, and turnip. MeRecon, available at http://bml.sbi.kmutt.ac.th/MeRecon, comprises 259 reactions (199 EC numbers), 1,052 proteins (870 genes) and 259 metabolites in eight sub-metabolisms. Analysis of MeRecon and the carbon assimilation pathways of the plant templates revealed the overall topology is highly conserved, but variations at sub metabolism level were found in relation to complexity underlying each biochemical reaction, such as numbers of responsible enzymatic proteins and their evolved functions, which likely explain the distinct metabolic phenotype. Thus, this study provides insights into the network characteristics and mechanisms that regulate the synthesis of metabolic phenotypes of cassava.

Prediction of cassava protein interactome based on interolog method

Cassava is a starchy root crop whose role in food security becomes more significant nowadays. Together with the industrial uses for versatile purposes, demand for cassava starch is continuously growing. However, in-depth study to uncover the mystery of cellular regulation, especially the interaction between proteins, is lacking. To reduce the knowledge gap in protein-protein interaction (PPI), genome-scale PPI network of cassava was constructed using interolog-based method (MePPI-In, available at http://bml.sbi.kmutt.ac.th/ppi). The network was constructed from the information of seven template plants. The MePPI-In included 90,173 interactions from 7,209 proteins. At least, 39 percent of the total predictions were found with supports from gene/protein expression data, while further co-expression analysis yielded 16 highly promising PPIs. In addition, domain-domain interaction information was employed to increase reliability of the network and guide the search for more groups of promising PPIs. Moreover, the topology and functional content of MePPI-In was similar to the networks of Arabidopsis and rice. The potential contribution of MePPI-In for various applications, such as protein-complex formation and prediction of protein function, was discussed and exemplified. The insights provided by our MePPI-In would hopefully enable us to pursue precise trait improvement in cassava.

Rooting for cassava: insights into photosynthesis and associated physiology as a route to improve yield potential

Contents 50 I. 50 II. 52 III. 54 IV. 55 V. 57 VI. 57 VII. 59 60 References 61 SUMMARY: As a consequence of an increase in world population, food demand is expected to grow by up to 110% in the next 30-35 yr. The population of sub-Saharan Africa is projected to increase by > 120%. In this region, cassava (Manihot esculenta) is the second most important source of calories and contributes c. 30% of the daily calorie requirements per person. Despite its importance, the average yield of cassava in Africa has not increased significantly since 1961. An evaluation of modern cultivars of cassava showed that the interception efficiency (ɛ_i ) of photosynthetically active radiation (PAR) and the efficiency of conversion of that intercepted PAR (ɛ_c ) are major opportunities for genetic improvement of the yield potential. This review examines what is known of the physiological processes underlying productivity in cassava and seeks to provide some strategies and directions toward yield improvement through genetic alterations to physiology to increase ɛ_i and ɛ_c . Possible physiological limitations, as well as environmental constraints, are discussed.

Transcriptome response of cassava leaves under natural shade

Cassava is an important staple crop in tropical and sub-tropical areas. As a common farming practice, cassava is usually cultivated intercropping with other crops and subjected to various degrees of shading, which causes reduced productivity. Herein, a comparative transcriptomic analysis was performed on a series of developmental cassava leaves under both full sunlight and natural shade conditions. Gene expression profiles of these two conditions exhibited similar developmental transitions, e.g. genes related to cell wall and basic cellular metabolism were highly expressed in immature leaves, genes involved in lipid metabolism and tetrapyrrole synthesis were highly expressed during the transition stages, and genes related to photosynthesis and carbohydrates metabolism were highly expressed in mature leaves. Compared with the control, shade significantly induced the expression of genes involved in light reaction of photosynthesis, light signaling and DNA synthesis/chromatin structure; however, the genes related to anthocyanins biosynthesis, heat shock, calvin cycle, glycolysis, TCA cycle, mitochondrial electron transport, and starch and sucrose metabolisms were dramatically depressed. Moreover, the shade also influenced the expression of hormone-related genes and transcriptional factors. The findings would improve our understanding of molecular mechanisms of shade response, and shed light on pathways associated with shade-avoidance syndrome for cassava improvement.

Characters related to higher starch accumulation in cassava storage roots

Cassava (Manihot esculenta) is valued mainly for high content starch in its roots. Our understanding of mechanisms promoting high starch accumulation in the roots is, however, still very limited. Two field-grown cassava cultivars, Huanan 124(H124) with low root starch and Fuxuan 01(F01) with high root starch, were characterised comparatively at four main growth stages. Changes in key sugars in the leaves, stems and roots seemed not to be strongly associated with the final amount of starch accumulated in the roots. However, when compared with H124, F01 exhibited a more compact arrangement of xylem vascular bundles in the leaf axils, much less callose around the phloem sieve plates in the stems, higher starch synthesis-related enzymatic activity but lower amylase activity in the roots, more significantly up-regulated expression of related genes, and a much higher stem flow rate (SFR). In conclusion, higher starch accumulation in the roots results from the concurrent effects of powerful stem transport capacity highlighted by higher SFR, high starch synthesis but low starch degradation in the roots, and high expression of sugar transporter genes in the stems. A model of high starch accumulation in cassava roots was therefore proposed and discussed.

Proteomics Profiling Reveals Carbohydrate Metabolic Enzymes and 14-3-3 Proteins Play Important Roles for Starch Accumulation during Cassava Root Tuberization

Cassava is one of the most important root crops as a reliable source of food and carbohydrates. Carbohydrate metabolism and starch accumulation in cassava storage root is a cascade process that includes large amounts of proteins and cofactors. Here, comparative proteomics were conducted in cassava root at nine developmental stages. A total of 154 identified proteins were found to be differentially expressed during starch accumulation and root tuberization. Many enzymes involved in starch and sucrose metabolism were significantly up-regulated, and functional classification of the differentially expressed proteins demonstrated that the majority were binding-related enzymes. Many proteins were took part in carbohydrate metabolism to produce energy. Among them, three 14-3-3 isoforms were induced to be clearly phosphorylated during storage root enlargement. Overexpression of a cassava 14-3-3 gene in Arabidopsis thaliana confirmed that the older leaves of these transgenic plants contained higher sugar and starch contents than the wild-type leaves. The 14-3-3 proteins and their binding enzymes may play important roles in carbohydrate metabolism and starch accumulation during cassava root tuberization. These results not only deepened our understanding of the tuberous root proteome, but also uncovered new insights into carbohydrate metabolism and starch accumulation during cassava root enlargement.

Genome-wide analysis reveals phytohormone action during cassava storage root initiation

Development of storage roots is a process associated with a phase change from cell division and elongation to radial growth and accumulation of massive amounts of reserve substances such as starch. Knowledge of the regulation of cassava storage root formation has accumulated over time; however, gene regulation during the initiation and early stage of storage root development is still poorly understood. In this study, transcription profiling of fibrous, intermediate and storage roots at eight weeks old were investigated using a 60-mer-oligo microarray. Transcription and gene expression were found to be the key regulating processes during the transition stage from fibrous to intermediate roots, while homeostasis and signal transduction influenced regulation from intermediate roots to storage roots. Clustering analysis of significant genes and transcription factors (TF) indicated that a number of phytohormone-related TF were differentially expressed; therefore, phytohormone-related genes were assembled into a network of correlative nodes. We propose a model showing the relationship between KNOX1 and phytohormones during storage root initiation. Exogeneous treatment of phytohormones N (6) -benzylaminopurine and 1-Naphthaleneacetic acid were used to induce the storage root initiation stage and to investigate expression patterns of the genes involved in storage root initiation. The results support the hypothesis that phytohormones are acting in concert to regulate the onset of cassava storage root development. Moreover, MeAGL20 is a factor that might play an important role at the onset of storage root initiation when the root tip becomes swollen.

Transcriptional analysis of South African cassava mosaic virus-infected susceptible and tolerant landraces of cassava highlights differences in resistance, basal defense and cell wall associated genes during infection

BACKGROUND

Cassava mosaic disease is caused by several distinct geminivirus species, including South African cassava mosaic virus-[South Africa:99] (SACMV). To date, there is limited gene regulation information on viral stress responses in cassava, and global transcriptome profiling in SACMV-infected cassava represents an important step towards understanding natural host responses to plant geminiviruses.

RESULTS

A RNA-seq time course (12, 32 and 67 dpi) study, monitoring gene expression in SACMV-challenged susceptible (T200) and tolerant (TME3) cassava landraces, was performed using the Applied Biosystems (ABI) SOLiD next-generation sequencing platform. The multiplexed paired end sequencing run produced a total of 523 MB and 693 MB of paired-end reads for SACMV-infected susceptible and tolerant cDNA libraries, respectively. Of these, approximately 50.7% of the T200 reads and 55.06% of TME3 reads mapped to the cassava reference genome available in phytozome. Using a log2 fold cut-off (p<0.05), comparative analysis between the six normalized cDNA libraries showed that 4181 and 1008 transcripts in total were differentially expressed in T200 and TME3, respectively, across 12, 32 and 67 days post infection, compared to mock-inoculated. The number of responsive transcripts increased dramatically from 12 to 32 dpi in both cultivars, but in contrast, in T200 the levels did not change significantly at 67 dpi, while in TME3 they declined. GOslim functional groups illustrated that differentially expressed genes in T200 and TME3 were overrepresented in the cellular component category for stress-related genes, plasma membrane and nucleus. Alterations in the expression of other interesting genes such as transcription factors, resistance (R) genes, and histone/DNA methylation-associated genes, were observed. KEGG pathway analysis uncovered important altered metabolic pathways, including phenylpropanoid biosynthesis, sucrose and starch metabolism, and plant hormone signalling.

CONCLUSIONS

Molecular mechanisms for TME3 tolerance are proposed, and differences in patterns and levels of transcriptome profiling between T200 and TME3 with susceptible and tolerant phenotypes, respectively, support the hypothesis that viruses rearrange their molecular interactions in adapting to hosts with different genetic backgrounds.

Increased expression of native cytosolic Cu/Zn superoxide dismutase and ascorbate peroxidase improves tolerance to oxidative and chilling stresses in cassava (Manihot esculenta Crantz)

BACKGROUND

Cassava (Manihot esculenta Crantz) is a tropical root crop, and is therefore, extremely sensitive to low temperature; its antioxidative response is pivotal for its survival under stress. Timely turnover of reactive oxygen species (ROS) in plant cells generated by chilling-induced oxidative damages, and scavenging can be achieved by non-enzymatic and enzymatic reactions in order to maintain ROS homeostasis.

RESULTS

Transgenic cassava plants that co-express cytosolic superoxide dismutase (SOD), MeCu/ZnSOD, and ascorbate peroxidase (APX), MeAPX2, were produced and tested for tolerance against oxidative and chilling stresses. The up-regulation of MeCu/ZnSOD and MeAPX2 expression was confirmed by the quantitative reverse transcriptase-polymerase chain reaction, and enzymatic activity analyses in the leaves of transgenic cassava plant lines with a single-transgene integration site. Upon exposure to ROS-generating agents, 100 μM ROS-generating reagent methyl viologen and 0.5 M H₂O₂, higher levels of enzymatic activities of SOD and APX were detected in transgenic plants than the wild type. Consequently, the oxidative stress parameters, such as lipid peroxidation, chlorophyll degradation and H₂O₂ synthesis, were lower in the transgenic lines than the wild type. Tolerance to chilling stress at 4°C for 2 d was greater in transgenic cassava, as observed by the higher levels of SOD, catalase, and ascorbate-glutathione cycle enzymes (e.g., APX, monodehydroascorbate reductase, dehydroascorbate reducatase and glutathione reductase) and lower levels of malondialdehyde content.

CONCLUSIONS

These results suggest that the expression of native cytosolic SOD and APX simultaneously activated the antioxidative defense mechanisms via cyclic ROS scavenging, thereby improving its tolerance to cold stress.

Next-generation sequencing-based mRNA and microRNA expression profiling analysis revealed pathways involved in the rapid growth of developing culms in Moso bamboo

BACKGROUND

As one of the fastest-growing lignocellulose-abundant plants on Earth, bamboos can reach their final height quickly due to the expansion of individual internodes already present in the buds; however, the molecular processes underlying this phenomenon remain unclear. Moso bamboo (Phyllostachys heterocycla cv. Pubescens) internodes from four different developmental stages and three different internodes within the same stage were used in our study to investigate the molecular processes at the transcriptome and post-transcriptome level.

RESULTS

Our anatomical observations indicated the development of culms was dominated by cell division in the initial stages and by cell elongation in the middle and late stages. The four major endogenous hormones appeared to actively promote culm development. Using next-generation sequencing-based RNA-Seq, mRNA and microRNA expression profiling technology, we produced a transcriptome and post-transcriptome in possession of a large fraction of annotated Moso bamboo genes, and provided a molecular basis underlying the phenomenon of sequentially elongated internodes from the base to the top. Several key pathways such as environmental adaptation, signal transduction, translation, transport and many metabolisms were identified as involved in the rapid elongation of bamboo culms.

CONCLUSIONS

This is the first report on the temporal and spatial transcriptome and gene expression and microRNA profiling in a developing bamboo culms. In addition to gaining more insight into the unique growth characteristics of bamboo, we provide a good case study to analyze gene, microRNA expression and profiling of non-model plant species using high-throughput short-read sequencing. Also, we demonstrate that the integrated analysis of our multi-omics data, including transcriptome, post-transcriptome, proteome, yield more complete representations and additional biological insights, especially the complex dynamic processes occurring in Moso bamboo culms.

Physiological and molecular characterization of drought responses and identification of candidate tolerance genes in cassava

Cassava is an important root crop to resource-poor farmers in marginal areas, where its production faces drought stress constraints. Given the difficulties associated with cassava breeding, a molecular understanding of drought tolerance in cassava will help in the identification of markers for use in marker-assisted selection and genes for transgenic improvement of drought tolerance. This study was carried out to identify candidate drought-tolerance genes and expression-based markers of drought stress in cassava. One drought-tolerant (improved variety) and one drought-susceptible (farmer-preferred) cassava landrace were grown in the glasshouse under well-watered and water-stressed conditions. Their morphological, physiological and molecular responses to drought were characterized. Morphological and physiological measurements indicate that the tolerance of the improved variety is based on drought avoidance, through reduction of water loss via partial stomatal closure. Ten genes that have previously been biologically validated as conferring or being associated with drought tolerance in other plant species were confirmed as being drought responsive in cassava. Four genes (MeALDH, MeZFP, MeMSD and MeRD28) were identified as candidate cassava drought-tolerance genes, as they were exclusively up-regulated in the drought-tolerant genotype to comparable levels known to confer drought tolerance in other species. Based on these genes, we hypothesize that the basis of the tolerance at the cellular level is probably through mitigation of the oxidative burst and osmotic adjustment. This study provides an initial characterization of the molecular response of cassava to drought stress resembling field conditions. The drought-responsive genes can now be used as expression-based markers of drought stress tolerance in cassava, and the candidate tolerance genes tested in the context of breeding (as possible quantitative trait loci) and engineering drought tolerance in transgenics.

Transcriptome analysis using a high-density oligomicroarray under drought stress in various genotypes of cassava: an important tropical crop

Cassava is an important crop that provides food security and income generation in many tropical countries and is known for its adaptability to various environmental conditions. Despite its global importance, the development of cassava microarray tools has not been well established. Here, we describe the development of a 60-mer oligonucleotide Agilent microarray representing ∼20,000 cassava genes and how it can be applied to expression profiling under drought stress using three cassava genotypes (MTAI16, MECU72 and MPER417-003). Our results identified about 1300 drought stress up-regulated genes in cassava and indicated that cassava has similar mechanisms for drought stress response and tolerance as other plant species. These results demonstrate that our microarray is a useful tool for analysing the cassava transcriptome and that it is applicable for various cassava genotypes.

Transcriptome profiling of low temperature-treated cassava apical shoots showed dynamic responses of tropical plant to cold stress

BACKGROUND

Cassava is an important tropical root crop adapted to a wide range of environmental stimuli such as drought and acid soils. Nevertheless, it is an extremely cold-sensitive tropical species. Thus far, there is limited information about gene regulation and signalling pathways related to the cold stress response in cassava. The development of microarray technology has accelerated the study of global transcription profiling under certain conditions.

RESULTS

A 60-mer oligonucleotide microarray representing 20,840 genes was used to perform transcriptome profiling in apical shoots of cassava subjected to cold at 7°C for 0, 4 and 9 h. A total of 508 transcripts were identified as early cold-responsive genes in which 319 sequences had functional descriptions when aligned with Arabidopsis proteins. Gene ontology annotation analysis identified many cold-relevant categories, including 'Response to abiotic and biotic stimulus', 'Response to stress', 'Transcription factor activity', and 'Chloroplast'. Various stress-associated genes with a wide range of biological functions were found, such as signal transduction components (e.g., MAP kinase 4), transcription factors (TFs, e.g., RAP2.11), and reactive oxygen species (ROS) scavenging enzymes (e.g., catalase 2), as well as photosynthesis-related genes (e.g., PsaL). Seventeen major TF families including many well-studied members (e.g., AP2-EREBP) were also involved in the early response to cold stress. Meanwhile, KEGG pathway analysis uncovered many important pathways, such as 'Plant hormone signal transduction' and 'Starch and sucrose metabolism'. Furthermore, the expression changes of 32 genes under cold and other abiotic stress conditions were validated by real-time RT-PCR. Importantly, most of the tested stress-responsive genes were primarily expressed in mature leaves, stem cambia, and fibrous roots rather than apical buds and young leaves. As a response to cold stress in cassava, an increase in transcripts and enzyme activities of ROS scavenging genes and the accumulation of total soluble sugars (including sucrose and glucose) were also detected.

CONCLUSIONS

The dynamic expression changes reflect the integrative controlling and transcriptome regulation of the networks in the cold stress response of cassava. The biological processes involved in the signal perception and physiological response might shed light on the molecular mechanisms related to cold tolerance in tropical plants and provide useful candidate genes for genetic improvement.