Expression profiling of cassava storage roots reveals an active process of glycolysis/gluconeogenesis

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.

Temporal patterns of gene expression associated with tuberous root formation and development in sweetpotato (Ipomoea batatas)

BACKGROUND

The tuberous root of sweetpotato is undisputedly an important organ from agronomic and biological perspectives. Little is known regarding the regulatory networks programming tuberous root formation and development.

RESULTS

Here, as a first step toward understanding these networks, we analyzed and characterized the genome-wide transcriptional profiling and dynamics of sweetpotato root in seven distinct developmental stages using a customized microarray containing 39,724 genes. Analysis of these genes identified temporal programs of gene expression, including hundreds of transcription factor (TF) genes. We found that most genes active in roots were shared across all developmental stages, although significant quantitative changes in gene abundance were observed for 5,368 (including 435 TFs) genes. Clustering analysis of these differentially expressed genes pointed out six distinct expression patterns during root development. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that genes involved in different processes were enriched at specific stages of root development. In contrast with the large number of shared expressed genes in root development, each stage or period of root development has only a small number of specific genes. In total, 712 (including 27 TFs) and 1,840 (including 115 TFs) genes were identified as root-stage and root-period specific, respectively at the level of microarray. Several of the specific TF genes are known regulators of root development, including DA1-related protein, SHORT-ROOT and BEL1-like. The remaining TFs with unknown roles would also play critical regulatory roles during sweetpotato tuberous root formation and development.

CONCLUSIONS

The results generated in this study provided spatiotemporal patterns of root gene expression in support of future efforts for understanding the underlying molecular mechanism that control sweetpotato yield and quality.

Transition from somatic embryo to friable embryogenic callus in cassava: dynamic changes in cellular structure, physiological status, and gene expression profiles

Friable embryogenic callus (FEC) is considered as the most suitable material for efficient genetic transformation of cassava. Heavy genotype dependence of FEC induction and amenability to somaclonal variation limits the production and maintenance of reliable FEC. Identifying key elements involved in biological processes from somatic embryos (SEs) to FEC at different stages provides critical insights for FEC improvement. Cytological observation showed a dramatic change of subcellular structures among SEs, fresh FEC (FFEC), and old FEC (OFEC). Decrease of sucrose and increase of fructose and glucose were detected in OFEC. A total of 6871 differentially expressed genes (DEGs) were identified from SEs, FFEC, and OFEC by RNA-seq. Analysis of the DEGs showed that FEC induction was accompanied by the process of dedifferentiation, whereas the epigenetics modification occurred during the continuous subculturing process. The cell structure was reconstructed, mainly including the GO terms of "cell periphery" and "external encapsulating structure"; in parallel, the internal mechanisms changed correspondingly, including the biological process of glycolysis and metabolisms of alanine, aspartate, and glutamate. The significant reduction of genomic DNA methylation in OFEC indicated altered gene expression via chromatin modification. These results indicate that the induction and long-term subculture of FEC is a complicated biological process involving changes of genome modification, gene expression, and subcellular reconstruction. The findings will be useful for improving FEC induction and maintenance from farmer-preferred cassava cultivars recalcitrant to genetic transformation, hence improving cassava through genetic engineering.

EVALUACIÓN DE LA EXPRESIÓN DE GENES IMPLICADOS EN LA BIOSÍNTESIS DE ALMIDÓN EN DIFERENTES VARIEDADES DE YUCA

<p>Las raíces almacenadoras de yuca representan una fuente importante de almidón. La ruta metabólica del almidón ha sido reconstruida recientemente en yuca gracias a la liberación de la secuencia completa de su genoma. En este estudio se evaluó la expresión de los genes que codifican para las enzimas Pululanasa, Isoamilasa, α-amilasa, Enzima Desproporcionante, ADP-glucosa pirofoforilasa, Almidón sintasa unida al gránulo, Enzima ramificante del almidón y Sintasa soluble del almidón, en las raíces almacenadoras de plantas de 5 y 11 meses de edad, en un grupo de cinco variedades de yuca. Se evidenciaron diferencias importantes en la expresión de estos genes entre las variedades evaluadas y entre los dos tiempos. Las variedades CM523-7 y SM1219-2 presentaron uno de los niveles más altos de expresión para los genes ADP-glucosa pirofoforilasa y Almidón sintasa unida al gránulo mientras que el gen para α-amilasa fue el más bajo en estas dos variedades. Aunque  la variedad TMS60444 presentó niveles de expresión similares en genes implicados en la síntesis de almidón, fue la que presentó el mayor nivel de expresión de la α-amilasa. Estos datos se pueden correlacionar con el relativo bajo contenido de materia seca en esta variedad. Los datos de expresión génica presentados en este trabajo permitirán complementar información sobre actividad enzimática con miras a identificar los elementos más importantes en la acumulación diferencial de almidón entre variedades de yuca.</p><p><strong>ABSTRACT</strong></p><p>Cassava storage roots represent an important starch source. Recently, the starch metabolic pathway in cassava has been reconstructed thanks to the full release of its genome. In this study gene expression was evaluated for genes coding Pullulanase, Isoamylase, α-amylase, Deproportionating enzyme, ADP-glucose pyrophosphorylase, Granule bound starch synthase, Starch branching enzyme and Soluble starch synthase, in cassava storage roots 5 and 11 months old, in 5 cassava varieties. Important gene expression differences were detected both at the variety and time level. CM523-7 and SM1219-2 showed one of the highest expression levels for AGPase and GBSS genes, while α-amylase showed the lowest level in these two varieties. TMS60444 variety showed similar expression levels in starch biosynthesis-related genes, but conversely also showed the highest α-amylase expression. This correlates with the relative low dry-matter content in TMS60444. Gene expression data reported here will allow complementing actual information on enzymatic activity, in order to identify the most relevant factors in differential starch accumulation between cassava varieties.</p><br /><p> </p>

Cassava genome from a wild ancestor to cultivated varieties

Cassava is a major tropical food crop in the Euphorbiaceae family that has high carbohydrate production potential and adaptability to diverse environments. Here we present the draft genome sequences of a wild ancestor and a domesticated variety of cassava and comparative analyses with a partial inbred line. We identify 1,584 and 1,678 gene models specific to the wild and domesticated varieties, respectively, and discover high heterozygosity and millions of single-nucleotide variations. Our analyses reveal that genes involved in photosynthesis, starch accumulation and abiotic stresses have been positively selected, whereas those involved in cell wall biosynthesis and secondary metabolism, including cyanogenic glucoside formation, have been negatively selected in the cultivated varieties, reflecting the result of natural selection and domestication. Differences in microRNA genes and retrotransposon regulation could partly explain an increased carbon flux towards starch accumulation and reduced cyanogenic glucoside accumulation in domesticated cassava. These results may contribute to genetic improvement of cassava through better understanding of its biology.

Why mosaic? Gene expression profiling of African cassava mosaic virus-infected cassava reveals the effect of chlorophyll degradation on symptom development

Cassava mosaic disease, caused by cassava begomoviruses, is the most serious disease for cassava in Africa. However, the pathogenesis of this disease is poorly understood. We employed high throughput digital gene expression profiling based on the Illumina Solexa sequencing technology to investigate the global transcriptional response of cassava to African cassava mosaic virus infection. We found that 3,210 genes were differentially expressed in virus-infected cassava leaves. Gene ontology term and Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that genes implicated in photosynthesis were most affected, consistent with the chlorotic symptoms observed in infected leaves. The upregulation of chlorophyll degradation genes, including the genes encoding chlorophyllase, pheophytinase, and pheophorbide a oxygenase, and downregulation of genes encoding the major apoproteins in light-harvesting complex II were confirmed by qRT-PCR. These findings, together with the reduction of chlorophyll b content and fewer grana stacks in the infected leaf cells, reveal that the degradation of chlorophyll plays an important role in African cassava mosaic virus symptom development. This study will provide a road map for future investigations into viral pathogenesis.

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.

Starch biosynthesis in cassava: a genome-based pathway reconstruction and its exploitation in data integration

Background

Cassava is a well-known starchy root crop utilized for food, feed and biofuel production. However, the comprehension underlying the process of starch production in cassava is not yet available.

Results

In this work, we exploited the recently released genome information and utilized the post-genomic approaches to reconstruct the metabolic pathway of starch biosynthesis in cassava using multiple plant templates. The quality of pathway reconstruction was assured by the employed parsimonious reconstruction framework and the collective validation steps. Our reconstructed pathway is presented in the form of an informative map, which describes all important information of the pathway, and an interactive map, which facilitates the integration of omics data into the metabolic pathway. Additionally, to demonstrate the advantage of the reconstructed pathways beyond just the schematic presentation, the pathway could be used for incorporating the gene expression data obtained from various developmental stages of cassava roots. Our results exhibited the distinct activities of the starch biosynthesis pathway in different stages of root development at the transcriptional level whereby the activity of the pathway is higher toward the development of mature storage roots.

Conclusions

To expand its applications, the interactive map of the reconstructed starch biosynthesis pathway is available for download at the SBI group’s website (http://sbi.pdti.kmutt.ac.th/?page_id=33). This work is considered a big step in the quantitative modeling pipeline aiming to investigate the dynamic regulation of starch biosynthesis in cassava roots.

Transcriptional profiling of sweetpotato (Ipomoea batatas) roots indicates down-regulation of lignin biosynthesis and up-regulation of starch biosynthesis at an early stage of storage root formation

Background

The number of fibrous roots that develop into storage roots determines sweetpotato yield. The aim of the present study was to identify the molecular mechanisms involved in the initiation of storage root formation, by performing a detailed transcriptomic analysis of initiating storage roots using next-generation sequencing platforms. A two-step approach was undertaken: (1) generating a database for the sweetpotato root transcriptome using 454-Roche sequencing of a cDNA library created from pooled samples of two root types: fibrous and initiating storage roots; (2) comparing the expression profiles of initiating storage roots and fibrous roots, using the Illumina Genome Analyzer to sequence cDNA libraries of the two root types and map the data onto the root transcriptome database.

Results

Use of the 454-Roche platform generated a total of 524,607 reads, 85.6% of which were clustered into 55,296 contigs that matched 40,278 known genes. The reads, generated by the Illumina Genome Analyzer, were found to map to 31,284 contigs out of the 55,296 contigs serving as the database. A total of 8,353 contigs were found to exhibit differential expression between the two root types (at least 2.5-fold change). The Illumina-based differential expression results were validated for nine putative genes using quantitative real-time PCR. The differential expression profiles indicated down-regulation of classical root functions, such as transport, as well as down-regulation of lignin biosynthesis in initiating storage roots, and up-regulation of carbohydrate metabolism and starch biosynthesis. In addition, data indicated delicate control of regulators of meristematic tissue identity and maintenance, associated with the initiation of storage root formation.

Conclusions

This study adds a valuable resource of sweetpotato root transcript sequences to available data, facilitating the identification of genes of interest. This resource enabled us to identify genes that are involved in the earliest stage of storage root formation, highlighting the reduction in carbon flow toward phenylpropanoid biosynthesis and its delivery into carbohydrate metabolism and starch biosynthesis, as major events involved in storage root initiation. The novel transcripts related to storage root initiation identified in this study provide a starting point for further investigation into the molecular mechanisms underlying this process.

Genome-wide analysis of plant-specific Dof transcription factor family in tomato

The Dof (DNA binding with One Finger) family encoding single zinc finger proteins has been known as a family of plant-specific transcription factors. These transcription factors are involved in a variety of functions of importance for different biological processes in plants. In the current study, we identified 34 Dof family genes in tomato, distributed on 11 chromosomes. A complete overview of SlDof genes in tomato is presented, including the gene structures, chromosome locations, phylogeny, protein motifs and evolution pattern. Phylogenetic analysis of 34 SlDof proteins resulted in four classes constituting six clusters. In addition, a comparative analysis between these genes in tomato, Arabidopsis and rice was also performed. The tomato Dof family expansion has been dated to recent duplication events, and segmental duplication is predominant for the SlDof genes. Furthermore, the SlDof genes displayed differential expression either in their transcript abundance or in their expression patterns under normal growth conditions. This is the first step towards genome-wide analyses of the Dof genes in tomato. Our study provides a very useful reference for cloning and functional analysis of the members of this gene family in tomato and other species.

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.

Aldehyde dehydrogenase (ALDH) superfamily in plants: gene nomenclature and comparative genomics

In recent years, there has been a significant increase in the number of completely sequenced plant genomes. The comparison of fully sequenced genomes allows for identification of new gene family members, as well as comprehensive analysis of gene family evolution. The aldehyde dehydrogenase (ALDH) gene superfamily comprises a group of enzymes involved in the NAD(+)- or NADP(+)-dependent conversion of various aldehydes to their corresponding carboxylic acids. ALDH enzymes are involved in processing many aldehydes that serve as biogenic intermediates in a wide range of metabolic pathways. In addition, many of these enzymes function as 'aldehyde scavengers' by removing reactive aldehydes generated during the oxidative degradation of lipid membranes, also known as lipid peroxidation. Plants and animals share many ALDH families, and many genes are highly conserved between these two evolutionarily distinct groups. Conversely, both plants and animals also contain unique ALDH genes and families. Herein we carried out genome-wide identification of ALDH genes in a number of plant species-including Arabidopsis thaliana (thale crest), Chlamydomonas reinhardtii (unicellular algae), Oryza sativa (rice), Physcomitrella patens (moss), Vitis vinifera (grapevine) and Zea mays (maize). These data were then combined with previous analysis of Populus trichocarpa (poplar tree), Selaginella moellindorffii (gemmiferous spikemoss), Sorghum bicolor (sorghum) and Volvox carteri (colonial algae) for a comprehensive evolutionary comparison of the plant ALDH superfamily. As a result, newly identified genes can be more easily analyzed and gene names can be assigned according to current nomenclature guidelines; our goal is to clarify previously confusing and conflicting names and classifications that might confound results and prevent accurate comparisons between studies.

Large-scale collection and annotation of gene models for date palm (Phoenix dactylifera, L.)

The date palm (Phoenix dactylifera L.), famed for its sugar-rich fruits (dates) and cultivated by humans since 4,000 B.C., is an economically important crop in the Middle East, Northern Africa, and increasingly other places where climates are suitable. Despite a long history of human cultivation, the understanding of P. dactylifera genetics and molecular biology are rather limited, hindered by lack of basic data in high quality from genomics and transcriptomics. Here we report a large-scale effort in generating gene models (assembled expressed sequence tags or ESTs and mapped to a genome assembly) for P. dactylifera, using the long-read pyrosequencing platform (Roche/454 GS FLX Titanium) in high coverage. We built fourteen cDNA libraries from different P. dactylifera tissues (cultivar Khalas) and acquired 15,778,993 raw sequencing reads-about one million sequencing reads per library-and the pooled sequences were assembled into 67,651 non-redundant contigs and 301,978 singletons. We annotated 52,725 contigs based on the plant databases and 45 contigs based on functional domains referencing to the Pfam database. From the annotated contigs, we assigned GO (Gene Ontology) terms to 36,086 contigs and KEGG pathways to 7,032 contigs. Our comparative analysis showed that 70.6 % (47,930), 69.4 % (47,089), 68.4 % (46,441), and 69.3 % (47,048) of the P. dactylifera gene models are shared with rice, sorghum, Arabidopsis, and grapevine, respectively. We also assigned our gene models into house-keeping and tissue-specific genes based on their tissue specificity.

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.

Cassava genetic transformation and its application in breeding

As a major source of food, cassava (Manihot esculenta Crantz) is an important root crop in the tropics and subtropics of Africa and Latin America, and serves as raw material for the production of starches and bioethanol in tropical Asia. Cassava improvement through genetic engineering not only overcomes the high heterozygosity and serious trait separation that occurs in its traditional breeding, but also quickly achieves improved target traits. Since the first report on genetic transformation in cassava in 1996, the technology has gradually matured over almost 15 years of development and has overcome cassava genotype constraints, changing from mode cultivars to farmer-preferred ones. Significant progress has been made in terms of an increased resistance to pests and diseases, biofortification, and improved starch quality, building on the fundamental knowledge and technologies related to planting, nutrition, and the processing of this important food crop that has often been neglected. Therefore, cassava has great potential in food security and bioenergy development worldwide.