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

Quantitative Trait Locus Mapping for Plant Height and Branch Number in CCRI70 Recombinant Inbred Line Population of Upland Cotton (Gossypium hirsutum)

Upland cotton accounts for a high percentage (95%) of the world's cotton production. Plant height (PH) and branch number (BN) are two important agronomic traits that have an impact on improving the level of cotton mechanical harvesting and cotton yield. In this research, a recombinant inbred line (RIL) population with 250 lines developed from the variety CCRI70 was used for constructing a high-density genetic map and identification of quantitative trait locus (QTL). The results showed that the map harbored 8298 single nucleotide polymorphism (SNP) markers, spanning a total distance of 4876.70 centimorgans (cMs). A total of 69 QTLs for PH (9 stable) and 63 for BN (11 stable) were identified and only one for PH was reported in previous studies. The QTLs for PH and BN harbored 495 and 446 genes, respectively. Combining the annotation information, expression patterns and previous studies of these genes, six genes could be considered as potential candidate genes for PH and BN. The results could be helpful for cotton researchers to better understand the genetic mechanism of PH and BN development, as well as provide valuable genetic resources for cotton breeders to manipulate cotton plant architecture to meet future demands.

Comprehensive analysis of the aldehyde dehydrogenase gene family in Phaseolus vulgaris L. and their response to saline-alkali stress

Background

Aldehyde dehydrogenase (ALDH) scavenges toxic aldehyde molecules by catalyzing the oxidation of aldehydes to carboxylic acids. Although ALDH gene family members in various plants have been extensively studied and were found to regulate plant response to abiotic stress, reports on ALDH genes in the common bean (Phaseolus vulgaris L.) are limited. In this study, we aimed to investigate the effects of neutral (NS) and basic alkaline (AS) stresses on growth, physiological and biochemical indices, and ALDH activity, ALDH gene expression of common bean. In addition, We used bioinformatics techniques to analyze the physical and chemical properties, phylogenetic relationships, gene replication, collinearity, cis-acting elements, gene structure, motifs, and protein structural characteristics of PvALDH family members.

Results

We found that both NS and AS stresses weakened the photosynthetic performance of the leaves, induced oxidative stress, inhibited common bean growth, and enhanced the antioxidative system to scavenge reactive oxygen species. Furthermore, we our findings revealed that ALDH in the common bean actively responds to NS or AS stress by inducing the expression of PvALDH genes. In addition, using the established classification criteria and phylogenetic analysis, 27 PvALDHs were identified in the common bean genome, belonging to 10 ALDH families. The primary expansion mode of PvALDH genes was segmental duplication. Cis-acting elemental analysis showed that PvALDHs were associated with abiotic stress and phytohormonal responses. Gene expression analysis revealed that the PvALDH gene expression was tissue-specific. For instance, PvALDH3F1 and PvALDH3H1 were highly expressed in flower buds and flowers, respectively, whereas PvALDH3H2 and PvALDH2B4 were highly expressed in green mature pods and young pods, respectively. PvALDH22A1 and PvALDH11A2 were highly expressed in leaves and young trifoliates, respectively; PvALDH18B2 and PvALDH18B3 were highly expressed in stems and nodules, respectively; and PvALDH2C2 and PvALDH2C3 were highly expressed in the roots. PvALDHs expression in the roots responded positively to NS-AS stress, and PvALDH2C3, PvALDH5F1, and PvALDH10A1 were significantly (P < 0.05) upregulated in the roots.

Conclusion

These results indicate that AS stress causes higher levels of oxidative damage than NS stress, resulting in weaker photosynthetic performance and more significant inhibition of common bean growth. The influence of PvALDHs potentially modulates abiotic stress response, particularly in the context of saline-alkali stress. These findings establish a basis for future research into the potential roles of ALDHs in the common bean.

Trehalose metabolism coordinates transcriptional regulatory control and metabolic requirements to trigger the onset of cassava storage root initiation

Cassava storage roots (SR) are an important source of food energy and raw material for a wide range of applications. Understanding SR initiation and the associated regulation is critical to boosting tuber yield in cassava. Decades of transcriptome studies have identified key regulators relevant to SR formation, transcriptional regulation and sugar metabolism. However, there remain uncertainties over the roles of the regulators in modulating the onset of SR development owing to the limitation of the widely applied differential gene expression analysis. Here, we aimed to investigate the regulation underlying the transition from fibrous (FR) to SR based on Dynamic Network Biomarker (DNB) analysis. Gene expression analysis during cassava root initiation showed the transition period to SR happened in FR during 8 weeks after planting (FR8). Ninety-nine DNB genes associated with SR initiation and development were identified. Interestingly, the role of trehalose metabolism, especially trehalase1 (TRE1), in modulating metabolites abundance and coordinating regulatory signaling and carbon substrate availability via the connection of transcriptional regulation and sugar metabolism was highlighted. The results agree with the associated DNB characters of TRE1 reported in other transcriptome studies of cassava SR initiation and Attre1 loss of function in literature. The findings help fill the knowledge gap regarding the regulation underlying cassava SR initiation.

Integrating Full-Length and Second-Generation Transcriptomics to Reveal Differentially Expressed Genes Associated with the Development of Corydalis yanhusuo Tuber

Corydalis yanhusuo is a medicinal herb in China that has been widely used to treat various kinds of pain. The tuber is the main organ of C. yanhusuo used for medicinal purposes, but changes in related genes during the development of the tuber have rarely been reported. To identify the differentially expressed genes during tuber development, C. yanhusuo full-length transcriptomic sequencing was performed using single-molecule real-time technology, and tubers at three development stages were selected for comparative transcriptome analysis. A total of 90,496 full-length non-chimeric transcripts were obtained, and 19,341 transcripts were annotated in at least one public database. A total of 9221 differentially expressed genes were identified during the swelling process of C. yanhusuo tuber. A Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis revealed that differentially expressed genes associated with a "starch and sucrose metabolism pathway", "phenylpropanoid biosynthesis pathway", "isoquinoline alkaloid biosynthesis pathway", "zeatin biosynthesis pathway", and "brassinosteroid biosynthesis pathway" were predominantly enriched. In addition, the genes involved in cell wall metabolism were potentially associated with tuber swelling. These processes regulated and were involved in C. yanhusuo tuber development. The results provide a foundation for further research on tuber formation in medicinal plants.

Genome-wide characterization of aldehyde dehydrogenase gene family members in groundnut (Arachis hypogaea) and the analysis under saline-alkali stress

Groundnut or peanut (Arachis hypogaea) is a legume crop. Its seeds are rich in protein and oil. Aldehyde dehydrogenase (ALDH, EC: 1.2.1.3) is an important enzyme involved in detoxification of aldehyde and cellular reactive oxygen species, as well as in attenuation of lipid peroxidation-meditated cellular toxicity under stress conditions. However, few studies have been identified and analyzed about ALDH members in Arachis hypogaea. In the present study, 71 members of the ALDH superfamily (AhALDH) were identified using the reference genome obtained from the Phytozome database. A systematic analysis of the evolutionary relationship, motif, gene structure, cis-acting elements, collinearity, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, and expression patterns was conducted to understand the structure and function of AhALDHs. AhALDHs exhibited tissue-specific expression, and quantitative real-time PCR identified significant differences in the expression levels of AhALDH members under saline-alkali stress. The results revealed that some AhALDHs members could be involved in response to abiotic stress. Our findings on AhALDHs provide insights for further study.

Integrative omics approaches revealed a crosstalk among phytohormones during tuberous root development in cassava

Integrative omics approaches revealed a crosstalk among phytohormones during tuberous root development in cassava. Tuberous root formation is a complex process consisting of phase changes as well as cell division and elongation for radial growth. We performed an integrated analysis to clarify the relationships among metabolites, phytohormones, and gene transcription during tuberous root formation in cassava (Manihot esculenta Crantz). We also confirmed the effects of the auxin (AUX), cytokinin (CK), abscisic acid (ABA), jasmonic acid (JA), gibberellin (GA), brassinosteroid (BR), salicylic acid, and indole-3-acetic acid conjugated with aspartic acid on tuberous root development. An integrated analysis of metabolites and gene expression indicated the expression levels of several genes encoding enzymes involved in starch biosynthesis and sucrose metabolism are up-regulated during tuberous root development, which is consistent with the accumulation of starch, sugar phosphates, and nucleotides. An integrated analysis of phytohormones and gene transcripts revealed a relationship among AUX signaling, CK signaling, and BR signaling, with AUX, CK, and BR inducing tuberous root development. In contrast, ABA and JA inhibited tuberous root development. These phenomena might represent the differences between stem tubers (e.g., potato) and root tubers (e.g., cassava). On the basis of these results, a phytohormonal regulatory model for tuberous root development was constructed. This model may be useful for future phytohormonal studies involving cassava.

Tuber and Tuberous Root Development

Root and tuber crops have been an important part of human nutrition since the early days of humanity, providing us with essential carbohydrates, proteins, and vitamins. Today, they are especially important in tropical and subtropical regions of the world, where they help to feed an ever-growing population. Early induction and storage organ size are important agricultural traits, as they determine yield over time. During potato tuberization, environmental and metabolic status are sensed, ensuring proper timing of tuberization mediated by phloem-mobile signals. Coordinated cellular restructuring and expansion growth, as well as controlled storage metabolism in the tuber, are executed. This review summarizes our current understanding of potato tuber development and highlights similarities and differences to important tuberous root crop species like sweetpotato and cassava. Finally, we point out knowledge gaps that need to be filled before a complete picture of storage organ development can emerge.

Transcriptome integrated metabolic modeling of carbon assimilation underlying storage root development in cassava

The existing genome-scale metabolic model of carbon metabolism in cassava storage roots, rMeCBM, has proven particularly resourceful in exploring the metabolic basis for the phenotypic differences between high and low-yield cassava cultivars. However, experimental validation of predicted metabolic fluxes by carbon labeling is quite challenging. Here, we incorporated gene expression data of developing storage roots into the basic flux-balance model to minimize infeasible metabolic fluxes, denoted as rMeCBMx, thereby improving the plausibility of the simulation and predictive power. Three different conceptual algorithms, GIMME, E-Flux, and HPCOF were evaluated. The rMeCBMx-HPCOF model outperformed others in predicting carbon fluxes in the metabolism of storage roots and, in particular, was highly consistent with transcriptome of high-yield cultivars. The flux prediction was improved through the oxidative pentose phosphate pathway in cytosol, as has been reported in various studies on root metabolism, but hardly captured by simple FBA models. Moreover, the presence of fluxes through cytosolic glycolysis and alanine biosynthesis pathways were predicted with high consistency with gene expression levels. This study sheds light on the importance of prediction power in the modeling of complex plant metabolism. Integration of multi-omics data would further help mitigate the ill-posed problem of constraint-based modeling, allowing more realistic simulation.

Comparison of carbohydrate partitioning and expression patterns of some genes involved in carbohydrate biosynthesis pathways in annual and biennial species of Cichorium spp

Variation in metabolism and partitioning of carbohydrates, particularly fructans, between annual and perennial Cichorium species remains a challenging topic. To address this problem, an annual (endive, Cichorium endive L. var. Crispum; Asteraceae) and a biennial species (chicory, Cichorium intybus L. var. Witloof; Asteraceae) were compared with in terms of variability in carbohydrate accumulation and expression patterns of fructan-active enzyme genes, as well as sucrose metabolism at various growth and developmental stages. In general, constituents such as 1-kestose, nystose, and inulin were detected only in the root of chicory and were not present in any of the endive tissues. For both species, flower tissue contained maximum levels of both fructose and glucose, while for sucrose, more fluctuations were observed. On the other hand, all the genes under study exhibited variation, not only between the two species but also among different tissues at different sampling times. In endive root compared to endive leaf, the expression of cell wall invertase genes and sucrose accumulation decreased simultaneously, indicating the limited capacity of its roots to absorb sucrose, a precursor to inulin production. In addition, low expression of fructan: fructan fructosyltransferase in endive root compared to chicory root confirmed the inability of endive to inulin synthesis. Overall, annual and biennial species were different in the production of inulin, transport, remobilization, and unloading of sucrose.

De novo transcriptome in roots of switchgrass (Panicum virgatum L.) reveals gene expression dynamic and act network under alkaline salt stress

Background

Soil salinization is a major limiting factor for crop cultivation. Switchgrass is a perennial rhizomatous bunchgrass that is considered an ideal plant for marginal lands, including sites with saline soil. Here we investigated the physiological responses and transcriptome changes in the roots of Alamo (alkaline-tolerant genotype) and AM-314/MS-155 (alkaline-sensitive genotype) under alkaline salt stress.

Results

Alkaline salt stress significantly affected the membrane, osmotic adjustment and antioxidant systems in switchgrass roots, and the ASTTI values between Alamo and AM-314/MS-155 were divergent at different time points. A total of 108,319 unigenes were obtained after reassembly, including 73,636 unigenes in AM-314/MS-155 and 65,492 unigenes in Alamo. A total of 10,219 DEGs were identified, and the number of upregulated genes in Alamo was much greater than that in AM-314/MS-155 in both the early and late stages of alkaline salt stress. The DEGs in AM-314/MS-155 were mainly concentrated in the early stage, while Alamo showed greater advantages in the late stage. These DEGs were mainly enriched in plant-pathogen interactions, ubiquitin-mediated proteolysis and glycolysis/gluconeogenesis pathways. We characterized 1480 TF genes into 64 TF families, and the most abundant TF family was the C2H2 family, followed by the bZIP and bHLH families. A total of 1718 PKs were predicted, including CaMK, CDPK, MAPK and RLK. WGCNA revealed that the DEGs in the blue, brown, dark magenta and light steel blue 1 modules were associated with the physiological changes in roots of switchgrass under alkaline salt stress. The consistency between the qRT-PCR and RNA-Seq results confirmed the reliability of the RNA-seq sequencing data. A molecular regulatory network of the switchgrass response to alkaline salt stress was preliminarily constructed on the basis of transcriptional regulation and functional genes.

Conclusions

Alkaline salt tolerance of switchgrass may be achieved by the regulation of ion homeostasis, transport proteins, detoxification, heat shock proteins, dehydration and sugar metabolism. These findings provide a comprehensive analysis of gene expression dynamic and act network induced by alkaline salt stress in two switchgrass genotypes and contribute to the understanding of the alkaline salt tolerance mechanism of switchgrass and the improvement of switchgrass germplasm.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07368-w.

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.

Cd exposure-induced growth retardation involves in energy metabolism disorder of midgut tissues in the gypsy moth larvae

Cadmium, a common environmental contaminant in both terrestrial and aquatic ecosystems, presented a serious hazard to growth and development of phytophagous insects. For better understanding the toxicology of Cd exposure on phytophagous insects, the physiological and molecular mechanisms underlying the energy metabolism disorder in midgut tissue of gypsy moth larvae fed on Cd-amended artificial diets (3.248 or 44.473 mg Cd/kg fresh food) were investigated. Our results showed that compared with control, Cd exposure at both two levels triggered detriment effects on growth indexes, and with the increase of exposure concentrations, the adverse effects were significantly exacerbated. Larval growth and nutritional indexes (except approximate digestibility) showed a strong positive correlation, indicating that growth retardation in the gypsy moth larvae under Cd stress was tightly related to the food utilization. The key genes at mRNA level in glycolysis/gluconeogenesis, citrate cycle pathway and starch/sucrose metabolism pathway also presented a significant and positive correlation with growth indexes, once again demonstrating that energy metabolism was the key factor that controls the growth and development of the gypsy moth larvae under Cd stress. Antioxidant system collapse and oxidative damage, a chief cause of histopathological alterations in midgut tissue, consist of the physiological basis of energy metabolism disorder in Cd-treated gypsy moth larvae. Together, these results suggest that histopathological alterations or oxidative damage of tissue structure significant disturbed physiological functions of midgut tissue in gypsy moth larvae exposed to Cd stress, as reflected via food utilization or energy metabolism disorder, and eventually resulted in larval growth retardation.

The Cassava Source-Sink project: opportunities and challenges for crop improvement by metabolic engineering

Cassava (Manihot esculenta Crantz) is one of the important staple foods in Sub-Saharan Africa. It produces starchy storage roots that provide food and income for several hundred million people, mainly in tropical agriculture zones. Increasing cassava storage root and starch yield is one of the major breeding targets with respect to securing the future food supply for the growing population of Sub-Saharan Africa. The Cassava Source-Sink (CASS) project aims to increase cassava storage root and starch yield by strategically integrating approaches from different disciplines. We present our perspective and progress on cassava as an applied research organism and provide insight into the CASS strategy, which can serve as a blueprint for the improvement of other root and tuber crops. Extensive profiling of different field-grown cassava genotypes generates information for leaf, phloem, and root metabolic and physiological processes that are relevant for biotechnological improvements. A multi-national pipeline for genetic engineering of cassava plants covers all steps from gene discovery, cloning, transformation, molecular and biochemical characterization, confined field trials, and phenotyping of the seasonal dynamics of shoot traits under field conditions. Together, the CASS project generates comprehensive data to facilitate conventional breeding strategies for high-yielding cassava genotypes. It also builds the foundation for genome-scale metabolic modelling aiming to predict targets and bottlenecks in metabolic pathways. This information is used to engineer cassava genotypes with improved source-sink relations and increased yield potential.

Genome-wide identification and evolution of Dof transcription factor family in cultivated and ancestral cotton species

The DNA-binding with one finger (Dof) proteins are transcription factors involved in many biological processes in plants. To predict the evolutionary pattern, a genome-wide in-silico analysis of Dof TFs family in diploid (Gossypium arboreum and Gossypium raimondii) and allotetraploid (Gossypium hirsutum and Gossypium barbadense) cotton species were carried out. In G. arboreum, we have identified 58 non-redundant genes encoding Dof proteins renamed as GaDof (G. arboreum Dof), 55 Dof genes were identified in G. raimondii (GrDof), 89 were predicted ffrom G. hirsutum (GhDof) and the highest, 110 Dof genes were identified in G. barbadense (GbDof). The phylogenetic analysis, physical location, gene structure, conserved domain analyses were also investigated for G. arboreum, G. raimondii, and G. hirsutum. The gene expression pattern in G. hirsutum, at different growth stages, revealing the probable involvement of some GhDof genes in growth and development. These genes may improve seed germination and growth in cotton.

Gene Regulatory Network Guided Investigations and Engineering of Storage Root Development in Root Crops

The plasticity of plant development relies on its ability to balance growth and stress resistance. To do this, plants have established highly coordinated gene regulatory networks (GRNs) of the transcription factors and signaling components involved in developmental processes and stress responses. In root crops, yields of storage roots are mainly determined by secondary growth driven by the vascular cambium. In relation to this, a dynamic yet intricate GRN should operate in the vascular cambium, in coordination with environmental changes. Despite the significance of root crops as food sources, GRNs wired to mediate secondary growth in the storage root have just begun to emerge, specifically with the study of the radish. Gene expression data available with regard to other important root crops are not detailed enough for us directly to infer underlying molecular mechanisms. Thus, in this review, we provide a general overview of the regulatory programs governing the development and functions of the vascular cambium in model systems, and the role of the vascular cambium on the growth and yield potential of the storage roots in root crops. We then undertake a reanalysis of recent gene expression data generated for major root crops and discuss common GRNs involved in the vascular cambium-driven secondary growth in storage roots using the wealth of information available in Arabidopsis. Finally, we propose future engineering schemes for improving root crop yields by modifying potential key nodes in GRNs.

Production of very-high-amylose cassava by post-transcriptional silencing of branching enzyme genes

High amylose starch can be produced by plants deficient in the function of branching enzymes (BEs). Here we report the production of transgenic cassava (Manihot esculenta Crantz) with starches containing up to 50% amylose due to the constitutive expression of hair-pin dsRNAs targeting the BE1 or BE2 genes. All BE1-RNAi plant lines (BE1i) and BE2-RNAi plant lines (BE2i) were grown up in the field, but with reduced total biomass production. Considerably high amylose content in the storage roots of BE2i plant lines was achieved. Storage starch granules of BE1i and BE2i plants had similar morphology as wild type (WT), however, the size of BE1i starch granules were bigger than that of WT. Comparisons of amylograms and thermograms of all three sources of storage starches revealed dramatic changes to the pasting properties and a higher melting temperature for BE2i starches. Glucan chain length distribution analysis showed a slight increase in chains of DP>36 in BE1i lines and a dramatic increase in glucan chains between DP 10-20 and DP>40 in BE2i lines. Furthermore, BE2i starches displayed a B-type X-ray diffraction pattern instead of the A-type pattern found in BE1i and WT starches. Therefore, cassava BE1 and BE2 function differently in storage root starch biosynthesis.

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.

Genetic characterization of cassava (Manihot esculenta Crantz) genotypes using agro-morphological and single nucleotide polymorphism markers

Dearth of information on extent of genetic variability in cassava limits the genetic improvement of cassava genotypes in Sierra Leone. The aim of this study was to assess the genetic diversity and relationships within 102 cassava genotypes using agro-morphological and single nucleotide polymorphism markers. Morphological classification based on qualitative traits categorized the germplasm into five different groups, whereas the quantitative trait set had four groups. The SNP markers classified the germplasm into three main cluster groups. A total of seven principal components (PCs) in the qualitative and four PCs in the quantitative trait sets accounted for 79.03% and 72.30% of the total genetic variation, respectively. Significant and positive correlations were observed between average yield per plant and harvest index (r = 0.76***), number of storage roots per plant and harvest index (r = 0.33*), height at first branching and harvest index (0.26*), number of storage roots per plant and average yield per plant (r = 0.58*), height at first branching and average yield per plant (r = 0.24*), length of leaf lobe and petiole length (r = 0.38*), number of leaf lobe and petiole length (r = 0.31*), width of leaf lobe and length of leaf lobe (r = 0.36*), number of leaf lobe and length of leaf lobe (r = 0.43*), starch content and dry matter content (r = 0.99***), number of leaf lobe and root dry matter (r = 0.30*), number of leaf lobe and starch content (r = 0.28*), and height at first branching and plant height (r = 0.45**). Findings are useful for conservation, management, short term recommendation for release and genetic improvement of the crop.

Melatonin as an inducer of arecoline and their coordinated roles in anti-oxidative activity and immune responses

Melatonin and it induced arecoline in arecoline play coordinated roles in immune responses.

Cell Wall Invertase 3 Affects Cassava Productivity via Regulating Sugar Allocation From Source to Sink

Storage roots are the main sink for photo-assimilate accumulation and reflect cassava yield and productivity. Regulation of sugar partitioning from leaves to storage roots has not been elucidated. Cell wall invertases are involved in the hydrolysis of sugar during phloem unloading of vascular plants to control plant development and sink strength but have rarely been studied in root crops like cassava. MeCWINV3 encodes a typical cell wall invertase in cassava and is mainly expressed in vascular bundles. The gene is highly expressed in leaves, especially mature leaves, in response to diurnal rhythm. When MeCWINV3 was overexpressed in cassava, sugar export from leaves to storage roots was largely inhibited and sucrose hydrolysis in leaves was accelerated, leading to increased transient starch accumulation by blocking starch degradation and reduced overall plant growth. The progress of leaf senescence was promoted in the MeCWINV3 over-expressed cassava plants with increased expression of senescence-related genes. Storage root development was also delayed because of dramatically reduced sugar allocation from leaves. As a result, the transcriptional expression of starch biosynthetic genes such as small subunit ADP-glucose pyrophosphorylase, granule-bound starch synthase I, and starch branching enzyme I was reduced in accordance with insufficient sugar supply in the storage roots of the transgenic plants. These results show that MeCWINV3 regulates sugar allocation from source to sink and maintains sugar balance in cassava, thus affecting yield of cassava storage roots.

Proteomics and microstructure profiling of goat milk protein after homogenization

This study investigated the protein changes in goat milk during the homogenization process using label-free quantification. We quantified 310 and 315 proteins in the control group (CG) and homogenized group (HG), respectively, and 16 proteins were significantly different between the 2 groups. For HG, the goat milk protein particle sizes were smaller and more evenly distributed and exhibited an increase in the regular arrangement of the secondary structures. Proteomics analysis verified that xanthine dehydrogenase and asparaginase-like 1 expression in CG were higher than in HG, whereas the opposite was observed for fructose-bisphosphate aldolase, κ-casein, and β-casein. Significant changes were found in the homogenization-treated goat milk proteome that were related to goat milk glycolysis/gluconeogenesis metabolism. This work provides updated information on the current proteome characteristics of homogenized goat milk, which may be important for applying the protein component of goat milk to human nutrition and health.

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.

Gibberellin Promotes Sweetpotato Root Vascular Lignification and Reduces Storage-Root Formation

Sweetpotato yield depends on a change in the developmental fate of adventitious roots into storage-roots. The mechanisms underlying this developmental switch are still unclear. We examined the hypothesis claiming that regulation of root lignification determines storage-root formation. We show that application of the plant hormone gibberellin increased stem elongation and root gibberellin levels, while having inhibitory effects on root system parameters, decreasing lateral root number and length, and significantly reducing storage-root number and diameter. Furthermore, gibberellin enhanced root xylem development, caused increased lignin deposition, and, at the same time, decreased root starch accumulation. In accordance with these developmental effects, gibberellin application upregulated expression levels of sweetpotato orthologues of Arabidopsis vascular development regulators (IbNA075, IbVND7, and IbSND2) and of lignin biosynthesis genes (IbPAL, IbC4H, Ib4CL, IbCCoAOMT, and IbCAD), while downregulating starch biosynthesis genes (IbAGPase and IbGBSS) in the roots. Interestingly, gibberellin downregulated root expression levels of orthologues of the Arabidopsis BREVIPEDICELLUS transcription factor (IbKN2 and IbKN3), regulator of meristem maintenance. The results substantiate our hypothesis and mark gibberellin as an important player in regulation of sweetpotato root development, suggesting that increased fiber formation and lignification inhibit storage-root formation and yield. Taken together, our findings provide insight into the mechanisms underlying sweetpotato storage-root formation and provide a valuable database of genes for further research.

Sugars en route to the roots. Transport, metabolism and storage within plant roots and towards microorganisms of the rhizosphere

In plants, the root is a typical sink organ that relies exclusively on the import of sugar from the aerial parts. Sucrose is delivered by the phloem to the most distant root tips and, en route to the tip, is used by the different root tissues for metabolism and storage. Besides, a certain portion of this carbon is exuded in the rhizosphere, supplied to beneficial microorganisms and diverted by parasitic microbes. The transport of sugars toward these numerous sinks either occurs symplastically through cell connections (plasmodesmata) or is apoplastically mediated through membrane transporters (MST, mononsaccharide tranporters, SUT/SUC, H+/sucrose transporters and SWEET, Sugar will eventually be exported transporters) that control monosaccharide and sucrose fluxes. Here, we review recent progresses on carbon partitioning within and outside roots, discussing membrane transporters involved in plant responses to biotic and abiotic factors.

Plant Vascular Tissues-Connecting Tissue Comes in All Shapes

For centuries, humans have grown and used structures based on vascular tissues in plants. One could imagine that life would have developed differently without wood as a resource for building material, paper, heating energy, or fuel and without edible tubers as a food source. In this review, we will summarise the status of research on Arabidopsis thaliana vascular development and subsequently focus on how this knowledge has been applied and expanded in research on the wood of trees and storage organs of crop plants. We will conclude with an outlook on interesting open questions and exciting new research opportunities in this growing and important field.

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.

Expression and Interaction Analysis among Saffron ALDHs and Crocetin Dialdehyde

In saffron, the cleavage of zeaxanthin by means of CCD2 generates crocetin dialdehyde, which is then converted by an unknown aldehyde dehydrogenase to crocetin. A proteome from saffron stigma was released recently and, based on the expression pattern and correlation analyses, five aldehyde dehydrogenases (ALDHs) were suggested as possible candidates to generate crocetin from crocetin dialdehydes. We selected four of the suggested ALDHs and analyzed their expression in different tissues, determined their activity over crocetin dialdehyde, and performed structure modeling and docking calculation to find their specificity. All the ALDHs were able to convert crocetin dialdehyde to crocetin, but two of them were stigma tissue-specific. Structure modeling and docking analyses revealed that, in all cases, there was a high coverage of residues in the models. All of them showed a very close conformation, indicated by the low root-mean-square deviation (RMSD) values of backbone atoms, which indicate a high similarity among them. However, low affinity between the enzymes and the crocetin dialdehyde were observed. Phylogenetic analysis and binding affinities calculations, including some ALDHs from Gardenia jasmonoides, Crocus sieberi, and Buddleja species that accumulate crocetin and Bixa orellana synthetizing the apocarotenoid bixin selected on their expression pattern matching with the accumulation of either crocins or bixin, pointed out that family 2 C4 members might be involved in the conversion of crocetin dialdehyde to crocetin with high specificity.

The root transcriptome of Achyranthes bidentata and the identification of the genes involved in the replanting benefit

The transcriptome profiling in replanting roots revealed that expression pattern changes of key genes promoted important metabolism pathways, antioxidant and pathogen defense systems, adjusted phytohormone signaling and inhibited lignin biosynthesis. The yield of the medicinal plant Achyranthes bidentata could be significantly increased when replanted into a field cultivated previously for the same crop, but the biological basis of this so-called "replanting benefit" is unknown. Here, the RNA-seq technique was used to identify candidate genes responsible for the benefit. The analysis of RNA-seq libraries prepared from mRNA extracted from the roots of first year planting (normal growth, NG) and second year replanting (consecutive monoculture, CM) yielded about 40.22 GB sequencing data. After de novo assembly, 87,256 unigenes were generated with an average length of 1060 bp. Among these unigenes, 55,604 were annotated with public databases, and 52,346 encoding sequences and 2881 transcription factors were identified. A contrast between the NG and CM libraries resulted in a set of 3899 differentially transcribed genes (DTGs). The DTGs related to the replanting benefit and their expression profiles were further analyzed by bioinformatics and qRT-PCR approaches. The major differences between the NG and CM transcriptomes included genes encoding products involved in glycolysis/gluconeogenesis, glutathione metabolism and antioxidant defense, in aspects of the plant/pathogen interaction, phytohormone signaling and phenylpropanoid biosynthesis. The indication was that replanting material enjoyed a stronger level of defense systems, a balance regulation of hormone signals and a suppression of lignin formation, thereby promoting root growth and development. The study provides considerable significant insights for a better understanding of the molecular mechanism of the replanting benefit and suggests their possible application in developing methods to reinforce the effects in medicinal plants.

Genome-Wide Identification, Expression, and Functional Analysis of the Sugar Transporter Gene Family in Cassava (Manihot esculenta)

The sugar transporter (STP) gene family encodes monosaccharide transporters that contain 12 transmembrane domains and belong to the major facilitator superfamily. STP genes play critical roles in monosaccharide distribution and participate in diverse plant metabolic processes. To investigate the potential roles of STPs in cassava (Manihot esculenta) tuber root growth, genome-wide identification and expression and functional analyses of the STP gene family were performed in this study. A total of 20 MeSTP genes (MeSTP1–20) containing the Sugar_tr conserved motifs were identified from the cassava genome, which could be further classified into four distinct groups in the phylogenetic tree. The expression profiles of the MeSTP genes explored using RNA-seq data showed that most of the MeSTP genes exhibited tissue-specific expression, and 15 out of 20 MeSTP genes were mainly expressed in the early storage root of cassava. qRT-PCR analysis further confirmed that most of the MeSTPs displayed higher expression in roots after 30 and 40 days of growth, suggesting that these genes may be involved in the early growth of tuber roots. Although all the MeSTP proteins exhibited plasma membrane localization, variations in monosaccharide transport activity were found through a complementation analysis in a yeast (Saccharomyces cerevisiae) mutant, defective in monosaccharide uptake. Among them, MeSTP2, MeSTP15, and MeSTP19 were able to efficiently complement the uptake of five monosaccharides in the yeast mutant, while MeSTP3 and MeSTP16 only grew on medium containing galactose, suggesting that these two MeSTP proteins are transporters specific for galactose. This study provides significant insights into the potential functions of MeSTPs in early tuber root growth, which possibly involves the regulation of monosaccharide distribution.

Regional Heritability Mapping Provides Insights into Dry Matter Content in African White and Yellow Cassava Populations

The HarvestPlus program for cassava ( Crantz) fortifies cassava with β-carotene by breeding for carotene-rich tubers (yellow cassava). However, a negative correlation between yellowness and dry matter (DM) content has been identified. We investigated the genetic control of DM in white and yellow cassava. We used regional heritability mapping (RHM) to associate DM with genomic segments in both subpopulations. Significant segments were subjected to candidate gene analysis and candidates were validated with prediction accuracies. The RHM procedure was validated via a simulation approach and revealed significant hits for white cassava on chromosomes 1, 4, 5, 10, 17, and 18, whereas hits for the yellow were on chromosome 1. Candidate gene analysis revealed genes in the carbohydrate biosynthesis pathway including plant serine-threonine protein kinases (SnRKs), UDP (uridine diphosphate)-glycosyltransferases, UDP-sugar transporters, invertases, pectinases, and regulons. Validation using 1252 unique identifiers from the SnRK gene family genome-wide recovered 50% of the predictive accuracy of whole-genome single nucleotide polymorphisms for DM, whereas validation using 53 likely genes (extracted from the literature) from significant segments recovered 32%. Genes including an acid invertase, a neutral or alkaline invertase, and a glucose-6-phosphate isomerase were validated on the basis of an a priori list for the cassava starch pathway, and also a fructose-biphosphate aldolase from the Calvin cycle pathway. The power of the RHM procedure was estimated as 47% when the causal quantitative trait loci generated 10% of the phenotypic variance (sample size = 451). Cassava DM genetics are complex and RHM may be useful for complex traits.

Proteomic analysis reveals the distinct energy and protein metabolism characteristics involved in myofiber type conversion and resistance of atrophy in the extensor digitorum longus muscle of hibernating Daurian ground squirrels

Previous hibernation studies demonstrated that such a natural model of skeletal muscle disuse causes limited muscle atrophy and a significant fast-to-slow fiber type shift. However, the underlying mechanism as defined in a large-scale analysis remains unclarified. Isobaric tags for relative and absolute quantification (iTRAQ) based quantitative analysis were used to examine proteomic changes in the fast extensor digitorum longus muscles (EDL) of Daurian ground squirrels (Spermophilus dauricus). Although the wet weights and fiber cross-sectional area of the EDL muscle showed no significant decrease, the percentage of slow type fiber was 61% greater (P < 0.01) in the hibernation group. Proteomics analysis identified 264 proteins that were significantly changed (ratio < 0.83 or >1.2-fold and P < 0.05) in the hibernation group, of which 23 proteins were categorized into energy production and conversion and translation and 22 proteins were categorized into ribosomal structure and biogenesis. Along with the validation by western blot, MAPKAP kinase 2, ATP5D, ACADSB, calcineurin, CSTB and EIF2S were up-regulated in the hibernation group, whereas PDK4, COX II and EIF3C were down-regulated in the hibernation group. MAPKAP kinase 2 and PDK4 were associated with glycolysis, COX II and ATP5D were associated with oxidative phosphorylation, ACADSB was associated with fatty acid metabolism, calcineurin and CSTB were associated with catabolism, and EIF2S and EIF3C were associated with anabolism. Moreover, the total proteolysis rate of EDL in the hibernation group was significantly inhibited compared with that in the pre-hibernation group. These distinct energy and protein metabolism characteristics may be involved in myofiber type conversion and resistance to atrophy in the EDL of hibernating Daurian ground squirrels.

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.

Genome-Wide Association Mapping of Correlated Traits in Cassava: Dry Matter and Total Carotenoid Content

Cassava is a starchy root crop cultivated in the tropics for fresh consumption and commercial processing. Primary selection objectives in cassava breeding include dry matter content and micronutrient density, particularly provitamin A carotenoids. These traits are negatively correlated in the African germplasm. This study aimed at identifying genetic markers associated with these traits and uncovering whether linkage and/or pleiotropy were responsible for observed negative correlation. A genome-wide association mapping using 672 clones genotyped at 72,279 single nucleotide polymorphism (SNP) loci was performed. Root yellowness was used indirectly to assess variation in carotenoid content. Two major loci for root yellowness were identified on chromosome 1 at positions 24.1 and 30.5 Mbp. A single locus for dry matter content that colocated with the 24.1 Mbp peak for carotenoids was identified. Haplotypes at these loci explained 70 and 37% of the phenotypic variability for root yellowness and dry matter content, respectively. Evidence of megabase-scale linkage disequilibrium (LD) around the major loci of the two traits and detection of the major dry matter locus in independent analysis for the white- and yellow-root subpopulations suggests that physical linkage rather that pleiotropy is more likely to be the cause of the negative correlation between the target traits. Moreover, candidate genes for carotenoid () and starch biosynthesis ( and ) occurred in the vicinity of the identified locus at 24.1 Mbp. These findings elucidate the genetic architecture of carotenoids and dry matter in cassava and provide an opportunity to accelerate breeding of these traits.

Alpha-Glucan, Water Dikinase 1 Affects Starch Metabolism and Storage Root Growth in Cassava (Manihot esculenta Crantz)

Regulation of storage root development by source strength remains largely unknown. The cassava storage root delay (srd) T-DNA mutant postpones storage root development but manifests normal foliage growth as wild-type plants. The SRD gene was identified as an orthologue of α-glucan, water dikinase 1 (GWD1), whose expression is regulated under conditions of light/dark cycles in leaves and is associated with storage root development. The GWD1-RNAi cassava plants showed both retarded plant and storage root growth, as a result of starch excess phenotypes with reduced photosynthetic capacity and decreased levels of soluble saccharides in their leaves. These leaves contained starch granules having greatly increased amylose content and type C semi-crystalline structures with increased short chains that suggested storage starch. In storage roots of GWD1-RNAi lines, maltose content was dramatically decreased and starches with much lower phosphorylation levels showed a drastically reduced β-amylolytic rate. These results suggested that GWD1 regulates transient starch morphogenesis and storage root growth by decreasing photo-assimilation partitioning from the source to the sink and by starch mobilization in root crops.

The industrial applications of cassava: current status, opportunities and prospects

Cassava (Manihot esculenta Crantz) is a drought-tolerant, staple food crop that is grown in tropical and subtropical areas. As an important raw material, cassava is a valuable food source in developing countries and is also extensively employed for producing starch, bioethanol and other bio-based products (e.g. feed, medicine, cosmetics and biopolymers). These cassava-based industries also generate large quantities of wastes/residues rich in organic matter and suspended solids, providing great potential for conversion into value-added products through biorefinery. However, the community of cassava researchers is relatively small and there is very limited information on cassava. Therefore this review summarizes current knowledge on the system biology, economic value, nutritional quality and industrial applications of cassava and its wastes in an attempt to accelerate understanding of the basic biology of cassava. The review also discusses future perspectives with respect to integrating and utilizing cassava information resources for increasing the economic and environmental sustainability of cassava industries. © 2017 Society of Chemical Industry.

Gene expression atlas for the food security crop cassava

Cassava (Manihot esculenta) feeds c. 800 million people world-wide. Although this crop displays high productivity under drought and poor soil conditions, it is susceptible to disease, postharvest deterioration and the roots contain low nutritional content. Here, we provide molecular identities for 11 cassava tissue/organ types through RNA-sequencing and develop an open access, web-based interface for further interrogation of the data. Through this dataset, we consider the physiology of cassava. Specifically, we focus on identification of the transcriptional signatures that define the massive, underground storage roots used as a food source and the favored target tissue for transgene integration and genome editing, friable embryogenic callus (FEC). Further, we identify promoters able to drive strong expression in multiple tissue/organs. The information gained from this study is of value for both conventional and biotechnological improvement programs.

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.

Transcriptome Analysis of Storage Roots and Fibrous Roots of the Traditional Medicinal Herb Callerya speciosa (Champ.) ScHot

Callerya speciosa (Champ.) ScHot is a woody perennial plant in Fabaceae, the roots of which are used medicinally. The storage roots of C. speciosa are derived from fibrous roots, but not all fibrous roots can develop into storage roots. To detect key genes involved in storage roots formation, we performed Illumina sequencing of the C. speciosa storage roots and fibrous roots. De novo assembly resulted in 161,926 unigenes, which were subsequently annotated by BLAST, GO and KEGG analyses. After expression profiling, 4538 differentially expressed genes were identified. The KEGG pathway enrichment analysis revealed changes in the biosynthesis of cytokinin, phenylpropanoid, starch, sucrose, flavone and other secondary metabolites. Transcription factor-related differentially expressed genes (DEGs) were also identified, including such gene families as GRAS, COL, MIKC, ERF, LBD, and NAC. The DEGs related to light signaling, starch, sugar, photohormones and cell wall-loosening might be involved in the formation of storage roots. This study provides the first transcriptome profiling of C. speciosa roots, data that will facilitate future research of root development and metabolites with medicinal value as well as the breeding of C. speciosa.

Cassava (Manihot esculenta) transcriptome analysis in response to infection by the fungus Colletotrichum gloeosporioides using an oligonucleotide-DNA microarray

Cassava anthracnose disease (CAD), caused by the fungus Colletotrichum gloeosporioides f. sp. Manihotis, is a serious disease of cassava (Manihot esculenta) worldwide. In this study, we established a cassava oligonucleotide-DNA microarray representing 59,079 probes corresponding to approximately 30,000 genes based on original expressed sequence tags and RNA-seq information from cassava, and applied it to investigate the molecular mechanisms of resistance to fungal infection using two cassava cultivars, Huay Bong 60 (HB60, resistant to CAD) and Hanatee (HN, sensitive to CAD). Based on quantitative real-time reverse transcription PCR and expression profiling by the microarray, we showed that the expressions of various plant defense-related genes, such as pathogenesis-related (PR) genes, cell wall-related genes, detoxification enzyme, genes related to the response to bacterium, mitogen-activated protein kinase (MAPK), genes related to salicylic acid, jasmonic acid and ethylene pathways were higher in HB60 compared with HN. Our results indicated that the induction of PR genes in HB60 by fungal infection and the higher expressions of defense response-related genes in HB60 compared with HN are likely responsible for the fungal resistance in HB60. We also showed that the use of our cassava oligo microarray could improve our understanding of cassava molecular mechanisms related to environmental responses and development, and advance the molecular breeding of useful cassava plants.

Transcriptomic analysis of Pseudostellariae Radix from different fields using RNA-seq

Pseudostellariae Radix is an important traditional Chinese medicine (TCM), which is consumed commonly for its positive health effects. However, a lack of transcriptomic and genomic information hinders research on Pseudostellariae Radix. Here, high-throughput RNA sequencing (RNA-seq) was employed for the de novo assembly to analyze the transcriptome in Pseudostellariae Radix, finding significantly differentially expressed genes in this TCM from different fields based on RNA-seq and bioinformatic analysis. A total of 146,408,539 paired-end reads were generated and assembled into 89,857 unigenes with an average length of 862bp. All of the assembly unigenes were annotated by running BLASTx and BLASTn similarity searches on the Non-redundant nucleotide database (NT), the Non-redundant protein database (NR), Swiss-Prot, Cluster of Orthologous Groups (COG), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Interpro. On the basis of bioinformatic analysis and the expression profiles for Pseudostellariae Radix, 29 significantly differentially expressed genes were identified, which provides the basic information for exploring the molecular mechanisms that determine the quality of Pseudostellariae Radix from different fields. The expression levels of 29 genes were validated by real-time quantitative PCR (RT-qPCR). This is the first study to sample Pseudostellariae Radix, which provides an invaluable resource for understanding the genome of this herb.

Physiological Investigation and Transcriptome Analysis of Polyethylene Glycol (PEG)-Induced Dehydration Stress in Cassava

Cassava is an important tropical and sub-tropical root crop that is adapted to drought environment. However, severe drought stress significantly influences biomass accumulation and starchy root production. The mechanism underlying drought-tolerance remains obscure in cassava. In this study, changes of physiological characters and gene transcriptome profiles were investigated under dehydration stress simulated by polyethylene glycol (PEG) treatments. Five traits, including peroxidase (POD) activity, proline content, malondialdehyde (MDA), soluble sugar and soluble protein, were all dramatically induced in response to PEG treatment. RNA-seq analysis revealed a gradient decrease of differentially expressed (DE) gene number in tissues from bottom to top of a plant, suggesting that cassava root has a quicker response and more induced/depressed DE genes than leaves in response to drought. Overall, dynamic changes of gene expression profiles in cassava root and leaves were uncovered: genes related to glycolysis, abscisic acid and ethylene biosynthesis, lipid metabolism, protein degradation, and second metabolism of flavonoids were significantly induced, while genes associated with cell cycle/organization, cell wall synthesis and degradation, DNA synthesis and chromatin structure, protein synthesis, light reaction of photosynthesis, gibberelin pathways and abiotic stress were greatly depressed. Finally, novel pathways in ABA-dependent and ABA-independent regulatory networks underlying PEG-induced dehydration response in cassava were detected, and the RNA-Seq results of a subset of fifteen genes were confirmed by real-time PCR. The findings will improve our understanding of the mechanism related to dehydration stress-tolerance in cassava and will provide useful candidate genes for breeding of cassava varieties better adapted to drought environment.

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.

Altered Phenylpropanoid Metabolism in the Maize Lc-Expressed Sweet Potato (Ipomoea batatas) Affects Storage Root Development

There is no direct evidence of the effect of lignin metabolism on early storage root development in sweet potato. In this study, we found that heterologous expression of the maize leaf color (Lc) gene in sweet potato increased anthocyanin pigment accumulation in the whole plant and resulted in reduced size with an increased length/width ratio, low yield and less starch content in the early storage roots. RT-PCR analysis revealed dramatic up-regulation of the genes involved in the lignin biosynthesis pathway in developing storage roots, leading to greater lignin content in the Lc transgenic lines, compared to the wild type. This was also evidenced by the enhanced lignification of vascular cells in the early storage roots. Furthermore, increased expression of the β-amylase gene in leaves and storage roots also accelerated starch degradation and increased the sugar use efficiency, providing more energy and carbohydrate sources for lignin biosynthesis in the Lc transgenic sweet potato. Lesser starch accumulation was observed in the developing storage roots at the initiation stage in the Lc plants. Our study provides experimental evidence of the basic carbohydrate metabolism underlying the development of storage roots, which is the transformation of lignin biosynthesis to starch biosynthesis.

Cassava root membrane proteome reveals activities during storage root maturation

Cassava (Manihot esculenta Crantz) is one of the most important crops of Thailand. Its storage roots are used as food, feed, starch production, and be the important source for biofuel and biodegradable plastic production. Despite the importance of cassava storage roots, little is known about the mechanisms involved in their formation. This present study has focused on comparison of the expression profiles of cassava root proteome at various developmental stages using two-dimensional gel electrophoresis and LC-MS/MS. Based on an anatomical study using Toluidine Blue, the secondary growth was confirmed to be essential during the development of cassava storage root. To investigate biochemical processes occurring during storage root maturation, soluble and membrane proteins were isolated from storage roots harvested from 3-, 6-, 9-, and 12-month-old cassava plants. The proteins with differential expression pattern were analysed and identified to be associated with 8 functional groups: protein folding and degradation, energy, metabolism, secondary metabolism, stress response, transport facilitation, cytoskeleton, and unclassified function. The expression profiling of membrane proteins revealed the proteins involved in protein folding and degradation, energy, and cell structure were highly expressed during early stages of development. Integration of these data along with the information available in genome and transcriptome databases is critical to expand knowledge obtained solely from the field of proteomics. Possible role of identified proteins were discussed in relation with the activities during storage root maturation in cassava.

Integrative Analysis of mRNA and miRNA Expression Profiles of the Tuberous Root Development at Seedling Stages in Turnips

The tuberous root of Brassica rapa L. (turnip) is an important modified organ for nutrition storage. A better understanding of the molecular mechanisms involved in the process of tuberous root development is of great value in both economic and biological context. In this study, we analyzed the expression profiles of both mRNAs and miRNAs in tuberous roots at an early stage before cortex splitting (ES), cortex splitting stage (CSS), and secondary root thickening stage (RTS) in turnip based on high-throughput sequencing technology. A large number of differentially expressed genes (DEGs) and several differentially expressed miRNAs (DEMs) were identified. Based on the DEG analysis, we propose that metabolism is the dominant pathway in both tuberous root initiation and secondary thickening process. The plant hormone signal transduction pathway may play a predominant role in regulating tuberous root initiation, while the starch and sucrose metabolism may be more important for the secondary thickening process. These hypotheses were partially supported by sequential DEM analyses. Of all DEMs, miR156a, miR157a, and miR172a exhibited relatively high expression levels, and were differentially expressed in both tuberous root initiation and the secondary thickening process with the expression profiles negatively correlated with those of their target genes. Our results suggest that these miRNAs play important roles in tuberous root development in turnips.

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.