Integrated analysis of DNA methylome and transcriptome revealing epigenetic regulation of CRIR1-promoted cold tolerance.

BACKGROUND: DNA methylation contributes to the epigenetic regulation of nuclear gene expression, and is associated with plant growth, development, and stress responses. Compelling evidence has emerged that long non-coding RNA (lncRNA) regulates DNA methylation. Previous genetic and physiological evidence indicates that lncRNA-CRIR1 plays a positive role in the responses of cassava plants to cold stress. However, it is unclear whether global DNA methylation changes with CRIR1-promoted cold tolerance. RESULTS: In this study, a comprehensive comparative analysis of DNA methylation and transcriptome profiles was performed to reveal the gene expression and epigenetic dynamics after CRIR1 overexpression. Compared with the wild-type plants, CRIR1-overexpressing plants present gained DNA methylation in over 37,000 genomic regions and lost DNA methylation in about 16,000 genomic regions, indicating a global decrease in DNA methylation after CRIR1 overexpression. Declining DNA methylation is not correlated with decreased/increased expression of the DNA methylase/demethylase genes, but is associated with increased transcripts of a few transcription factors, chlorophyll metabolism and photosynthesis-related genes, which could contribute to the CRIR1-promoted cold tolerance. CONCLUSIONS: In summary, a first set of transcriptome and epigenome data was integrated in this study to reveal the gene expression and epigenetic dynamics after CRIR1 overexpression, with the identification of several TFs, chlorophyll metabolism and photosynthesis-related genes that may be involved in CRIR1-promoted cold tolerance. Therefore, our study has provided valuable data for the systematic study of molecular insights for plant cold stress response.

Genome-Wide Survey of the RWP-RK Gene Family in Cassava (Manihot esculenta Crantz) and Functional Analysis.

The plant-specific RWP-RK transcription factor family plays a central role in the regulation of nitrogen response and gametophyte development. However, little information is available regarding the evolutionary relationships and characteristics of the RWP-RK family genes in cassava, an important tropical crop. Herein, 13 RWP-RK proteins identified in cassava were unevenly distributed across 9 of the 18 chromosomes (Chr), and these proteins were divided into two clusters based on their phylogenetic distance. The NLP subfamily contained seven cassava proteins including GAF, RWP-RK, and PB1 domains; the RKD subfamily contained six cassava proteins including the RWP-RK domain. Genes of the NLP subfamily had a longer sequence and more introns than the RKD subfamily. A large number of hormone- and stress-related cis-acting elements were found in the analysis of RWP-RK promoters. Real-time quantitative PCR revealed that all MeNLP1-7 and MeRKD1/3/5 genes responded to different abiotic stressors (water deficit, cold temperature, mannitol, polyethylene glycol, NaCl, and H2O2), hormonal treatments (abscisic acid and methyl jasmonate), and nitrogen starvation. MeNLP3/4/5/6/7 and MeRKD3/5, which can quickly and efficiently respond to different stresses, were found to be important candidate genes for further functional assays in cassava. The MeRKD5 and MeNLP6 proteins were localized to the cell nucleus in tobacco leaf. Five and one candidate proteins interacting with MeRKD5 and MeNLP6, respectively, were screened from the cassava nitrogen starvation library, including agamous-like mads-box protein AGL14, metallothionein 2, Zine finger FYVE domain containing protein, glyceraldehyde-3-phosphate dehydrogenase, E3 Ubiquitin-protein ligase HUWE1, and PPR repeat family protein. These results provided a solid basis to understand abiotic stress responses and signal transduction mediated by RWP-RK genes in cassava.

Biological function research of Fusarium oxysporum f. sp. cubense inducible banana long noncoding RNA Malnc2310 in Arabidopsis.

Long noncoding RNAs (lncRNAs) participate in plant biological processes under biotic and abiotic stresses. However, little is known about the function and regulation mechanism of lncRNAs related to the pathogen at a molecular level. A banana lncRNA, Malnc2310, is a Fusarium oxysporum f. sp. cubense inducible lncRNA in roots. In this study, we demonstrate the nuclear localization of Malnc2310 by fluorescence in situ hybridization and it can bind to several proteins that are related to flavonoid pathway, pathogen response and programmed cell death. Overexpression of Malnc2310 increases susceptibility to Fusarium crude extract (Fu), salinity, and cold in transgenic Arabidopsis. In addition, Malnc2310 transgenic Arabidopsis accumulated more anthocyanins under Fusarium crude extract and cold treatments that are related to upregulation of these genes involved in anthocyanin biosynthesis. Based on our findings, we propose that Malnc2310 may participate in flavonoid metabolism in plants under stress. Furthermore, phenylalanine ammonia lyase (PAL) protein expression was enhanced in Malnc2310 overexpressed transgenic Arabidopsis, and Malnc2310 may participate in PAL regulation by binding to it. This study provides new insights into the role of Malnc2310 in mediating plant stress adaptation.

An Oxalate Transporter Gene, AtOT, Enhances Aluminum Tolerance in Arabidopsis thaliana by Regulating Oxalate Efflux.

Secretion and efflux of oxalic acid from roots is an important aluminum detoxification mechanism for various plants; however, how this process is completed remains unclear. In this study, the candidate oxalate transporter gene AtOT, encoding 287 amino acids, was cloned and identified from Arabidopsis thaliana. AtOT was upregulated in response to aluminum stress at the transcriptional level, which was closely related to aluminum treatment concentration and time. The root growth of Arabidopsis was inhibited after knocking out AtOT, and this effect was amplified by aluminum stress. Yeast cells expressing AtOT enhanced oxalic acid resistance and aluminum tolerance, which was closely correlated with the secretion of oxalic acid by membrane vesicle transport. Collectively, these results underline an external exclusion mechanism of oxalate involving AtOT to enhance oxalic acid resistance and aluminum tolerance.

Turnip crinkle virus-encoded suppressor of RNA silencing interacts with Arabidopsis SGS3 to enhance virus infection.

Most plant viruses encode suppressors of RNA silencing (VSRs) to protect themselves from antiviral RNA silencing in host plants. The capsid protein (CP) of Turnip crinkle virus (TCV) is a well-characterized VSR, whereas SUPPRESSOR OF GENE SILENCING 3 (SGS3) is an important plant-encoded component of the RNA silencing pathways. Whether the VSR activity of TCV CP requires it to engage SGS3 in plant cells has yet to be investigated. Here, we report that TCV CP interacts with SGS3 of Arabidopsis in both yeast and plant cells. The interaction was identified with the yeast two-hybrid system, and corroborated with bimolecular fluorescence complementation and intracellular co-localization assays in Nicotiana benthamiana cells. While multiple partial TCV CP fragments could independently interact with SGS3, its hinge domain connecting the surface and protruding domains appears to be essential for this interaction. Conversely, SGS3 enlists its N-terminal domain and the XS rice gene X and SGS3 (XS) domain as the primary CP-interacting sites. Interestingly, SGS3 appears to stimulate TCV accumulation because viral RNA levels of a TCV mutant with low VSR activities decreased in the sgs3 knockout mutants, but increased in the SGS3-overexpressing transgenic plants. Transgenic Arabidopsis plants overexpressing TCV CP exhibited developmental abnormalities that resembled sgs3 knockout mutants and caused similar defects in the biogenesis of trans-acting small interfering RNAs. Our data suggest that TCV CP interacts with multiple RNA silencing pathway components that include SGS3, as well as previously reported DRB4 (dsRNA-binding protein 4) and AGO2 (ARGONAUTE protein 2), to achieve efficient suppression of RNA silencing-mediated antiviral defence.

The cassava (Manihot-esculenta Crantz)'s nitrate transporter NPF4.5, expressed in seedling roots, involved in nitrate flux and osmotic stress.

AtNPF4.5/AIT2, which was predicted to be a low-affinity transporter capable for nitrate uptake, was screened by ABA receptor complex in Arabidopsis ten years ago. However, the molecular and biochemical characterizations of AtNPF4.5 in plants remained largely unclear. In this study, the function of a plasma-membrane-localized and root-specifically-expressed gene MeNPF4.5 (Manihot-esculenta NITRATE TRANSPORTER 1 PTR FAMILY4.5), an ortholog of the Arabidopsis thaliana NPF4.5, was investigated in cassava roots as a nitrate efflux transporter on low nitrate medium and an influx transporter following exposure to high concentration of external nitrates. Moreover, RNA interference (RNAi) of MeNPF4.5 reduced the nitrate efflux capacity but the overexpressing cassava seedlings increased the ability of efflux from the elongation to the mature zone of root under low nitrate treatments. Besides, MeNPF4.5-RNAi expression reduced the nitrate influx capacity but enhanced nitrate absorption in parts of overexpressing plants from the meristem, elongation to mature zone of roots under high nitrate conditions. Furthermore, MeNPF4.5-RNAi seedlings survived owing to roots that could grow normally, but the MeNPF4.5-over-expressors showed adverse growth under 7% PEG6000 stress, suggesting that MeNPF4.5 negatively regulated the osmotic stress and was involved in nitrate flux through cassava seedlings.

Cloning, Expression Analysis, and Functional Characterization of Candidate Oxalate Transporter Genes of HbOT1 and HbOT2 from Rubber Tree (Hevea brasiliensis).

Secretion of oxalic acid from roots is an important aluminum detoxification mechanism for many plants such as Hevea brasiliensis (rubber tree). However, the underlying molecular mechanism and oxalate transporter genes in plants have not yet been reported. In this study, the oxalate transporter candidate genes HbOT1 and HbOT2 from the rubber tree were cloned and preliminarily identified. It was found that HbOT1 had a full length of 1163 bp with CDS size of 792 bp, encoding 263 amino acids, and HbOT2 had a full length of 1647 bp with a CDS region length of 840 bp, encoding 279 amino acid residues. HbOT1 and HbOT2 were both stable hydrophobic proteins with transmembrane structure and SNARE_assoc domains, possibly belonging to the SNARE_assoc subfamily proteins of the SNARE superfamily. qRT-PCR assays revealed that HbOT1 and HbOT2 were constitutively expressed in different tissues, with HbOT1 highly expressed in roots, stems, barks, and latex, while HbOT2 was highly expressed in latex. In addition, the expressions of HbOT1 and HbOT2 were up-regulated in response to aluminum stress, and they were inducible by metals, such as copper and manganese. Heterologous expression of HbOT1 and HbOT2 in the yeast mutant AD12345678 enhanced the tolerance to oxalic acid and high concentration aluminum stress, which was closely correlated with the secretion of oxalic acid. This study is the first report on oxalate transporter genes in plants, which provides a theoretical reference for the study on the molecular mechanism of oxalic acid secretion to relieve aluminum toxicity and on aluminum-tolerance genetic engineering breeding.

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

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

Role of cassava CC-type glutaredoxin MeGRXC3 in regulating sensitivity to mannitol-induced osmotic stress dependent on its nuclear activity.

BACKGROUND: We previously identified six drought-inducible CC-type glutaredoxins in cassava cultivars, however, less is known about their potential role in the molecular mechanism by which cassava adapted to abiotic stress. RESULTS: Herein, we investigate one of cassava drought-responsive CC-type glutaredoxins, namely MeGRXC3, that involved in regulation of mannitol-induced inhibition on seed germination and seedling growth in transgenic Arabidopsis. MeGRXC3 overexpression up-regulates several stress-related transcription factor genes, such as PDF1.2, ERF6, ORA59, DREB2A, WRKY40, and WRKY53 in Arabidopsis. Protein interaction assays show that MeGRXC3 interacts with Arabidopsis TGA2 and TGA5 in the nucleus. Eliminated nuclear localization of MeGRXC3 failed to result mannitol-induced inhibition of seed germination and seedling growth in transgenic Arabidopsis. Mutation analysis of MeGRXC3 indicates the importance of conserved motifs for its transactivation activity in yeast. Additionally, these motifs are also indispensable for its functionality in regulating mannitol-induced inhibition of seed germination and enhancement of the stress-related transcription factors in transgenic Arabidopsis. CONCLUSIONS: MeGRXC3 overexpression confers mannitol sensitivity in transgenic Arabidopsis possibly through interaction with TGA2/5 in the nucleus, and nuclear activity of MeGRXC3 is required for its function.

Global identification of full-length cassava lncRNAs unveils the role of cold-responsive intergenic lncRNA 1 in cold stress response.

Long noncoding RNAs (lncRNAs) have been considered to be important regulators of gene expression in a range of biological processes in plants. A large number of lncRNAs have been identified in plants. However, most of their biological functions still remain to be determined. Here, we identified a total of 3004 lncRNAs in cassava under normal or cold-treated conditions from Iso-seq data. We further characterized a cold-responsive intergenic lncRNA 1 (CRIR1) as a novel positive regulator of the plant response to cold stress. CRIR1 can be significantly induced by cold treatment. Ectopic expression of CRIR1 in cassava enhanced the cold tolerance of transgenic plants. Transcriptome analysis demonstrated that CRIR1 regulated a range of cold stress-related genes in a CBF-independent pathway. We further found that CRIR1 RNA can interact with cassava cold shock protein 5 (MeCSP5), which acts as an RNA chaperone, indicating that CRIR1 may recruit MeCSP5 to improve the translation efficiency of messenger RNA. In summary, our study extends the repertoire of lncRNAs in plants as well as their role in cold stress responses. Moreover, it reveals a mechanism by which CRIR1 affected cold stress response by modulating the expression of stress-responsive genes and increasing their translational yield.

MeSPL9 attenuates drought resistance by regulating JA signaling and protectant metabolite contents in cassava.

KEY MESSAGE: Analysis of drought-related genes in cassava shows the involvement of MeSPL9 in drought stress tolerance and overexpression of a dominant-negative form of this gene demonstrates its negative roles in drought stress resistance. Drought stress severely impairs crop yield and is considered a primary threat to food security worldwide. Although the SQUAMOSA promoter binding protein-like 9 (SPL9) gene participates extensively in numerous developmental processes and in plant response to abiotic stimuli, its role and regulatory pathway in cassava (Manihot esculenta) response to the drought condition remain elusive. In the current study, we show that cassava SPL9 (MeSPL9) plays negative roles in drought stress resistance. MeSPL9 expression was strongly repressed by drought treatment. Overexpression of a dominant-negative form of miR156-resistant MeSPL9, rMeSPL9-SRDX, in which a 12-amino acid repressor sequence was fused to rMeSPL9 at the C terminus, conferred drought tolerance without penalizing overall growth. rMeSPL9-SRDX-overexpressing lines not only exhibited increased osmoprotectant metabolites including proline and anthocyanin, but also accumulated more endogenous jasmonic acid (JA) and soluble sugars. Transcriptomic and real-time PCR analysis suggested that differentially expressed genes were involved in sugar or JA biosynthesis, signaling, and metabolism in transgenic cassava under drought conditions. Exogenous application of JA further confirmed that JA conferred improved drought resistance and promoted stomatal closure in cassava leaves. Taken together, our findings suggest that MeSPL9 affects drought resistance by modulating protectant metabolite levels and JA signaling, which have substantial implications for engineering drought tolerant crops.

Genome-wide survey of the phosphofructokinase family in cassava and functional characterization in response to oxygen-deficient stress.

BACKGROUND: Glycolytic pathway is common in all plant organs, especially in oxygen-deficient tissues. Phosphofructokinase (PFK) is a rate-limiting enzyme in the glycolytic pathway and catalyses the phosphorylation of fructose-6-phosphate to fructose-1,6-bisphosphate. Cassava (M. esculenta) root is a huge storage organ with low amount of oxygen. However, less is known about the functions of PFK from M. esculenta (MePFK). We conducted a systematic analysis of MePFK genes to explore the function of the MePFK gene family under hypoxic stress. RESULTS: We identified 13 MePFK genes and characterised their sequence structure. The phylogenetic tree divided the 13 genes into two groups: nine were MePFKs and four were pyrophosphate-fructose-6-phosphate phosphotransferase (MePFPs). We confirmed by green fluorescent protein fusion protein expression that MePFK03 and MePFPA1 were localised in the chloroplast and cytoplasm, respectively. The expression profiles of the 13 MePFKs detected by quantitative reverse transcription polymerase chain reaction revealed that MePFK02, MePFK03, MePFPA1, MePFPB1 displayed higher expression in leaves, root and flower. The expression of MePFK03, MePFPA1 and MePFPB1 in tuber root increased gradually with plant growth. We confirmed that hypoxia occurred in the cassava root, and the concentration of oxygen was sharply decreasing from the outside to the inside root. The expression of MePFK03, MePFPA1 and MePFPB1 decreased with the decrease in the oxygen concentration in cassava root. Waterlogging stress treatment showed that the transcript level of PPi-dependent MePFP and MeSuSy were up-regulated remarkably and PPi-dependent glycolysis bypass was promoted. CONCLUSION: A systematic survey of phylogenetic relation, molecular characterisation, chromosomal and subcellular localisation and cis-element prediction of MePFKs were performed in cassava. The expression profiles of MePFKs in different development stages, organs and under waterlogging stress showed that MePFPA1 plays an important role during the growth and development of cassava. Combined with the transcriptional level of MeSuSy, we found that pyrophosphate (PPi)-dependent glycolysis bypass was promoted when cassava was under waterlogging stress. The results would provide insights for further studying the function of MePFKs under hypoxic stress.

Molecular diversity analysis, drought related marker-traits association mapping and discovery of excellent alleles for 100-day old plants by EST-SSRs in cassava germplasms (Manihot esculenta Cranz).

Cassava is the third largest food crop of the world and has strong ability of drought tolerance. In order to evaluate the molecular diversity and to discover novel alleles for drought tolerance in cassava germplasms, we examined a total of 107 abiotic stress related expressed sequence tags-simple sequence repeat (EST-SSR) markers in 134 cassava genotypes coming from planting regions worldwide and performed drought related marker-traits association mapping. As results, we successfully amplified 98 of 107 markers in 97 polymorphic loci and 279 alleles, with 2.87 alleles per locus, gene diversity of 0.48 and polymorphic information content (PIC) of 0.41 on average. The genetic coefficient between every two lines was 0.37 on average, ranging from 0.21 to 0.82. According to our population structure analysis, these samples could be divided into three sub-populations showing obvious gene flow between them. We also performed water stress experiments using 100-day old cassava plants in two years and calculated the drought tolerance coefficients (DTCs) and used them as phenotypes for marker-trait association mapping. We found that 53 markers were significantly associated with these drought-related traits, with a contribution rate for trait variation of 8.60% on average, ranging between 2.66 and 28.09%. Twenty-four of these 53 associated genes showed differential transcription or protein levels which were confirmed by qRT-PCR under drought stress when compared to the control conditions in cassava. Twelve of twenty-four genes were the same differential expression patterns in omics data and results of qRT-PCR. Out of 33 marker-traits combinations on 24 loci, 34 were positive and 53 negative alleles according to their phenotypic effects and we also obtained the typical materials which carried these elite alleles. We also found 23 positive average allele effects while 10 loci were negative according to their allele effects (AAEs). Our results on molecular diversity, locus association and differential expression under drought can prove beneficial to select excellent materials through marker assisted selection and for functional genes research in the future.

Genome-wide identification and functional prediction of cold and/or drought-responsive lncRNAs in cassava.

Cold and drought stresses seriously affect cassava (Manihot esculenta) plant growth and yield. Recently, long noncoding RNAs (lncRNAs) have emerged as key regulators of diverse cellular processes in mammals and plants. To date, no systematic screening of lncRNAs under abiotic stress and their regulatory roles in cassava has been reported. In this study, we present the first reference catalog of 682 high-confidence lncRNAs based on analysis of strand-specific RNA-seq data from cassava shoot apices and young leaves under cold, drought stress and control conditions. Among them, 16 lncRNAs were identified as putative target mimics of cassava known miRNAs. Additionally, by comparing with small RNA-seq data, we found 42 lncNATs and sense gene pairs can generate nat-siRNAs. We identified 318 lncRNAs responsive to cold and/or drought stress, which were typically co-expressed concordantly or discordantly with their neighboring genes. Trans-regulatory network analysis suggested that many lncRNAs were associated with hormone signal transduction, secondary metabolites biosynthesis, and sucrose metabolism pathway. The study provides an opportunity for future computational and experimental studies to uncover the functions of lncRNAs in cassava.

Analysis of different strategies adapted by two cassava cultivars in response to drought stress: ensuring survival or continuing growth.

Cassava is one of the most drought-tolerant crops, however, the underlying mechanism for its ability to survive and produce under drought remains obscure. In this study, two cassava cultivars, SC124 and Arg7, were treated by gradually reducing the soil water content. Their responses to the drought stress were examined through their morphological and physiological traits and isobaric tags for relative and absolute quantitation (iTRAQ)-based proteomic analysis. SC124 plants adapted a 'survival' mode under mild drought stress as evidenced by early stomatal closure and a reduction in the levels of various photosynthetic proteins and photosynthetic capacity, resulting in early growth quiescence. In contrast, Arg7 plants underwent senescence of older leaves but continued to grow, although at a reduced rate, under mild drought. SC124 plants were more capable of surviving prolonged severe drought than Arg7. The iTRAQ analysis identified over 5000 cassava proteins. Among the drought-responsive proteins identified in the study were an aquaporin, myo-inositol 1-phosphate synthases, and a number of proteins involved in the antioxidant systems and secondary metabolism. Many proteins that might play a role in signalling or gene regulation were also identified as drought-responsive proteins, which included several protein kinases, two 14-3-3 proteins, several RNA-binding proteins and transcription factors, and two histone deacetylases. Our study also supports the notion that linamarin might play a role in nitrogen reallocation in cassava under drought.

Compositional equivalency of RNAi-mediated virus-resistant transgenic soybean and its nontransgenic counterpart.

RNA silencing or RNA interference (RNAi), which is triggered by double-stranded RNA (dsRNA), is an evolutionarily conserved process that is active in a wide variety of eukaryotic organisms. Engineering plants with hairpin construct in which the viral gene is arranged in inverted repeats (IR) renders plants resistant to plant virus infection. However, there is no report on whether biologically important changes occurred by the insertion of IR, which confer transgenic plants virus resistance. In the present study, the compositions of virus-resistant transgenic soybean seeds developed by insertion of three short IRs, each containing the specific, highly conserved sequences derived from one virus, were compared with those of nontransgenic counterparts by applying the principle of substantial equivalence to determine whether significant undesirable biological changes occurred by IR insertion. The results revealed that the nutrient components as well as antinutrient contents of these virus-resistant soybean lines are substantially equivalent to those of the nontransgenic counterparts, and the majority of the measured amounts of nutritional components and antinutrient contents are well within the range of values reported for other commercial soybean lines. The results imply that no biologically important changes occurred by the insertion of IRs in the RNAi-mediated virus-resistant transgenic soybeans. The results can serve as baseline information for developing RNAi-mediated transgenic soybean cultivars or other crops with broader spectrum virus resistance.

Recent advances in microbial raw starch degrading enzymes.

Raw starch degrading enzymes (RSDE) refer to enzymes that can directly degrade raw starch granules below the gelatinization temperature of starch. These promising enzymes can significantly reduce energy and simplify the process in starch industry. RSDE are ubiquitous and produced by plants, animals, and microorganisms. However, microbial sources are the most preferred one for large-scale production. During the past few decades, RSDE have been studied extensively. This paper reviews the recent development in the production, purification, properties, and application of microbial RSDE. This is the first review on microbial RSDE to date.