Changying Zeng

Conservation and divergence of microRNAs and their functions in Euphorbiaceous plants

MicroRNAs (miRNAs) are ∼21 nt non-coding RNAs which regulate post-transcriptional gene expression. miRNAs are key regulators of nearly all essential biological processes. Aiming at understanding miRNA’s functions in Euphorbiaceae, a large flowering plant family, we performed a genome-scale systematic study of miRNAs in Euphorbiaceae, by combining computational prediction and experimental analysis to overcome the difficulty of lack of genomes for most Euphorbiaceous species. Specifically, we predicted 85 conserved miRNAs in 23 families in the Castor bean (Ricinus communis), and experimentally verified and characterized 58 (68.2%) of the 85 miRNAs in at least one of four Euphorbiaceous species, the Castor bean, the Cassava (Manihot esculenta), the Rubber tree (Hevea brasiliensis) and the Jatropha (Jatropha curcas) during normal seedling development. To elucidate their function in stress response, we verified and profiled 48 (56.5%) of the 85 miRNAs under cold and drought stresses as well as during the processes of stress recovery. The results revealed some species- and condition-specific miRNA expression patterns. Finally, we predicted 258 miRNA:target partners, and identified the cleavage sites of six out of ten miRNA targets by a modified 5′ RACE. This study produced the first collection of miRNAs and their targets in Euphorbiaceae. Our results revealed wide conservation of many miRNAs and diverse functions in Euphorbiaceous plants during seedling growth and in response to abiotic stresses.

Cassava genome from a wild ancestor to cultivated varieties

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

Chilling acclimation provides immunity to stress by altering regulatory networks and inducing genes with protective functions in Cassava

BackgroundStress acclimation is an effective mechanism that plants acquired for adaption to dynamic environment. Even though generally considered to be sensitive to low temperature, Cassava, a major tropical crop, can be tolerant to much lower temperature after chilling acclimation. Improvement to chilling resistance could be beneficial to breeding. However, the underlying mechanism and the effects of chilling acclimation on chilling tolerance remain largely unexplored.ResultsIn order to understand the mechanism of chilling acclimation, we profiled and analyzed the transcriptome and microRNAome of Cassava, using high-throughput deep sequencing, across the normal condition, a moderate chilling stress (14°C), a harsh stress (4°C) after chilling acclimation (14°C), and a chilling shock from 24°C to 4°C. The results revealed that moderate stress and chilling shock triggered comparable degrees of transcriptional perturbation, and more importantly, about two thirds of differentially expressed genes reversed their expression from up-regulation to down-regulation or vice versa in response to hash stress after experiencing moderate stress. In addition, microRNAs played important roles in the process of this massive genetic circuitry rewiring. Furthermore, function analysis revealed that chilling acclimation helped the plant develop immunity to further harsh stress by exclusively inducing genes with function for nutrient reservation therefore providing protection, whereas chilling shock induced genes with function for viral reproduction therefore causing damage.ConclusionsOur study revealed, for the first time, the molecular basis of chilling acclimation, and showed potential regulation role of microRNA in chilling response and acclimation in Euphorbia.

Potential functions of microRNAs in starch metabolism and development revealed by miRNA transcriptome profiling of cassava cultivars and their wild progenitor

BackgroundMicroRNAs (miRNAs) are small (approximately 21 nucleotide) non-coding RNAs that are key post-transcriptional gene regulators in eukaryotic organisms. More than 100 cassava miRNAs have been identified in a conservation analysis and a repertoire of cassava miRNAs have also been characterised by next-generation sequencing (NGS) in recent studies. Here, using NGS, we profiled small non-coding RNAs and mRNA genes in two cassava cultivars and their wild progenitor to identify and characterise miRNAs that are potentially involved in plant growth and starch biosynthesis.ResultsSix small RNA and six mRNA libraries from leaves and roots of the two cultivars, KU50 and Arg7, and their wild progenitor, W14, were subjected to NGS. Analysis of the sequencing data revealed 29 conserved miRNA families and 33 new miRNA families. Together, these miRNAs potentially targeted a total of 360 putative target genes. Whereas 16 miRNA families were highly expressed in cultivar leaves, another 13 miRNA families were highly expressed in storage roots of cultivars. Co-expression analysis revealed that the expression level of some targets had negative relationship with their corresponding miRNAs in storage roots and leaves; these targets included MYB33, ARF10, GRF1, RD19, APL2, NF-YA3 and SPL2, which are known to be involved in plant development, starch biosynthesis and response to environmental stimuli.ConclusionThe identified miRNAs, target mRNAs and target gene ontology annotation all shed light on the possible functions of miRNAs in Manihot species. The differential expression of miRNAs between cultivars and their wild progenitor, together with our analysis of GO annotation and confirmation of miRNA: target pairs, might provide insight into know the differences between wild progenitor and cultivated cassava.

Endogenous small-noncoding RNAs and their roles in chilling response and stress acclimation in Cassava.

BackgroundSmall noncoding RNA (sncRNA), including microRNAs (miRNAs) and endogenous small-interfering RNAs (endo-siRNAs) are key gene regulators in eukaryotes, playing critical roles in plant development and stress tolerance. Trans-acting siRNAs (ta-siRNAs), which are secondary siRNAs triggered by miRNAs, and siRNAs from natural antisense transcripts (nat-siRNAs) are two well-studied classes of endo-siRNAs.ResultsIn order to understand sncRNAs’ roles in plant chilling response and stress acclimation, we performed a comprehensive study of miRNAs and endo-siRNAs in Cassava (Manihot esculenta), a major source of food for the world populations in tropical regions. Combining Next-Generation sequencing and computational and experimental analyses, we profiled and characterized sncRNA species and mRNA genes from the plants that experienced severe and moderate chilling stresses, that underwent further severe chilling stress after chilling acclimation at moderate stress, and that grew under the normal condition. We also included castor bean (Ricinus communis) in our study to understand conservation of sncRNAs. In addition to known miRNAs, we identified 32 (22 and 10) novel miRNAs as well as 47 (26 and 21) putative secondary siRNA-yielding and 8 (7 and 1) nat-siRNA-yielding candidate loci in Cassava and castor bean, respectively. Among the expressed sncRNAs, 114 miRNAs, 12 ta-siRNAs and 2 nat-siRNAs showed significant expression changes under chilling stresses.ConclusionSystematic and computational analysis of microRNAome and experimental validation collectively showed that miRNAs, ta-siRNAs, and possibly nat-siRNAs play important roles in chilling response and chilling acclimation in Cassava by regulating stress-related pathways, e.g. Auxin signal transduction. The conservation of these sncRNA might shed lights on the role of sncRNA-mediated pathways affected by chilling stress and stress acclimation in Euphorbiaceous plants.

MeSAUR1, Encoded by a Small Auxin-Up RNA Gene, Acts as a Transcription Regulator to Positively Regulate ADP-Glucose Pyrophosphorylase Small Subunit1a Gene in Cassava

Cassava, being one of the top three tuberous crops, features highly efficient starch accumulation in the storage root to adapt the tropical resources and environments. The molecular mechanism for the process, however, is still unclear. ADP-glucose pyrophosphorylase, the first and rate-limited enzyme in starch biosynthesis pathway, is a heterotetramer comprised of two small/catalytic and two large/modulatory subunits. To understand the regulation of MeAGPase, the promoter of a highly expressed small subunit, MeAGPs1a, was used as bait for a yeast one-hybrid assay to screen storage root cDNA library. One cDNA, coding for a small auxin-up RNA protein, named MeSAUR1, was isolated from cassava. MeSAUR1 could bind to the promoter of MeAGPS1a in yeast one-hybrid test and in vitro, and was located in cell nucleus. MeSAUR1 displayed a higher transcript level in cassava root cortex, and its expression was induced by indole-3-acetic acid, gibberellin and ethylene, but repressed by abscisic acid. A dual-luciferase interaction test further convinced that MeSAUR1 could bind to the promoter of MeAGPS1a, and positively regulate the transcription of MeAGPS1a in cassava.

Phylogeny and expression pattern analysis of TCP transcription factors in cassava seedlings exposed to cold and/or drought stress

The TCP transcription factors usually act as integrators of multiple growth regulatory and environmental stimuli. However, little is known about this gene family in the important tropical crop cassava (Manihot esculenta). In this study, 36 TCP genes were identified and renamed based on cassava whole-genome sequence and their sequence similarity with Arabidopsis TCPs. Typical TCP domains were detected in these proteins by multiple sequence alignment analysis. Evolutionary analysis indicated that MeTCPs could be divided into 8 subgroups, which was further supported by gene structure and conserved motif analyses. qRT-PCR analysis revealed tissue-specific and hormone-responsive expression patterns of MeTCP genes. Moreover, with global expression and promoter analysis, we found that MeTCPs showed similar or distinct expression patterns under cold and/or drought stress, suggesting that they might participate in distinct signaling pathways. Our study provides the first comprehensive analysis of TCP gene family in the cassava genome. The data will be useful for uncovering the potential functions of MeTCP genes, and their possible roles in mediating hormone and abiotic stress responses in cassava.

The Discrepant and Similar Responses of Genome-Wide Transcriptional Profiles between Drought and Cold Stresses in Cassava

Background: Cassava, an important tropical crop, has remarkable drought tolerance, but is very sensitive to cold. The growth, development, and root productivity of cassava are all adversely affected under cold and drought. Methods: To profile the transcriptional response to cold and drought stresses, cassava seedlings were respectively subjected to 0, 6, 24, and 48 h of cold stress and 0, 4, 6, and 10 days of drought stress. Their folded leaves, fully extended leaves, and roots were respectively investigated using RNA-seq. Results: Many genes specifically and commonly responsive to cold and drought were revealed: genes related to basic cellular metabolism, tetrapyrrole synthesis, and brassinosteroid metabolism exclusively responded to cold; genes related to abiotic stress and ethylene metabolism exclusively responded to drought; and genes related to cell wall, photosynthesis, and carbohydrate metabolism, DNA synthesis/chromatic structure, abscisic acid and salicylic acid metabolism, and calcium signaling commonly responded to both cold and drought. Discussion: Combined with cold- and/or drought-responsive transcription factors, the regulatory networks responding to cold and drought in cassava were constructed. All these findings will improve our understanding of the specific and common responses to cold and drought in cassava, and shed light on genetic improvement of cold and drought tolerance in cassava.

Ethylene Responsive Factor MeERF72 Negatively Regulates Sucrose synthase 1 Gene in Cassava

Cassava, an important food and industrial crop globally, is characterized by its powerful starch accumulation in its storage root. However, the underlying molecular mechanism for this feature remains unclear. Sucrose synthase initializes the conversion of sucrose to starch, and, to a certain extent, its enzyme activity can represent sink strength. To understand the modulation of MeSus gene family, the relatively high expressed member in storage root, MeSus1, its promoter was used as bait to screen cassava storage root full-length cDNA library through a yeast one-hybrid system. An ethylene responsive factor cDNA, designated as MeERF72 according to its homolog in Arabidopsis, was screened out. The transcript level of MeERF72 was induced by ethylene, drought, and salt treatments and repressed by abscisic acid, Auxin, gibberellin, salicylic acid, and low and high temperatures. The MeERF72 protein has a conserved APETALA2 domain in its N-terminus and an activated domain of 30 amino acids in its C-terminus, can bind to MeSus1 promoter in vitro and in vivo, and represses the promoter activity of MeSus1. MeERF72 is a transcription factor that can negatively regulate the expression level of MeSus1 in cassava.

MePHD1 as a PHD-Finger Protein Negatively Regulates ADP-Glucose Pyrophosphorylase Small Subunit1a Gene in Cassava

ADP-glucose pyrophosphorylase (AGPase) is an important enzyme in the starch synthesis pathway. Its enzyme activity can determine the efficiency of starch biosynthesis. Cassava (Manihot esculenta Crantz) is the main staple crop worldwide and has a high starch content in its storage root. However, the inner regulatory mechanism of AGPase gene family is unclear. MePHD1; a plant homeodomain transcription factor; was isolated through a yeast one-hybrid screening using the promoter of ADP-glucose pyrophosphorylase small subunit1a (MeAGPS1a) as bait, and cassava storage root cDNA library as prey. This factor could bind to the MeAGPS1a promoter in vitro and in vivo, and its predicted binding region ranged from −400 bp to −201 bp, at the translation initiation site. The transcript level of MePHD1 could be induced by five plant hormones, and a temperature of 42 °C. This was down-regulated during the maturation process of the storage root. MePHD1 protein could repress the promoter activity of MeAGPS1a gene by a dual-luciferase assay; which indicated that MePHD1 is a negative regulator for the transcript level of MeAGPS1a gene.