Ming Peng

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.

Analysis of banana transcriptome and global gene expression profiles in banana roots in response to infection by race 1 and tropical race 4 of Fusarium oxysporum f. sp. cubense

BackgroundCavendish, the most widely grown banana cultivar, is relatively resistant to Race 1 of Fusarium oxysporum f. sp. cubense (Foc1) which caused widespread Panama disease during the first half of the 20th century but is susceptible to Tropical Race 4 of Foc (Foc TR4) which is threatening world banana production. The genome of the diploid species Musa acuminata which is the ancestor of a majority of triploid banana cultivars has recently been sequenced. Availability of banana transcriptomes will be highly useful for improving banana genome annotation and for biological research. The knowledge of global gene expression patterns influenced by infection of different Foc races will help to understand the host responses to the infection.ResultsRNA samples from different organs of the Cavendish cultivar were pooled for deep sequencing using the Illumina technology. Analysis of the banana transcriptome led to identification of over 842 genes that were not annotated by the Musa genome project. A large number of simple nucleotide polymorphisms (SNPs) and short insertions and deletion (indels) were identified from the transcriptome data. GFP-expressing Foc1 and Foc TR4 were used to monitor the infection process. Both Foc1 and Foc TR4 were found to be able to invade banana roots and spread to root vascular tissues in the first two days following inoculation. Digital gene expression (DGE) profiling analysis reveal that the infection by Foc1 and Foc TR4 caused very similar changes in the global gene expression profiles in the banana roots during the first two days of infection. The Foc infection led to induction of many well-known defense-related genes. Two genes encoding the ethylene biosynthetic enzyme ACC oxidase and several ethylene-responsive transcription factors (ERF) were among the strongly induced genes by both Foc1 and Foc TR4.ConclusionsBoth Foc1 and Foc TR4 are able to spread into the vascular system of banana roots during the early infection process and their infection led to similar gene expression profiles in banana roots. The transcriptome profiling analysis indicates that the ethylene synthetic and signalling pathways were activated in response to the Foc infection.

Genome-wide characterization and analysis of bZIP transcription factor gene family related to abiotic stress in cassava.

The basic leucine zipper (bZIP) transcription factor family plays crucial roles in various aspects of biological processes. Currently, no information is available regarding the bZIP family in the important tropical crop cassava. Herein, 77 bZIP genes were identified from cassava. Evolutionary analysis indicated that MebZIPs could be divided into 10 subfamilies, which was further supported by conserved motif and gene structure analyses. Global expression analysis suggested that MebZIPs showed similar or distinct expression patterns in different tissues between cultivated variety and wild subspecies. Transcriptome analysis of three cassava genotypes revealed that many MebZIP genes were activated by drought in the root of W14 subspecies, indicating the involvement of these genes in the strong resistance of cassava to drought. Expression analysis of selected MebZIP genes in response to osmotic, salt, cold, ABA, and H2O2 suggested that they might participate in distinct signaling pathways. Our systematic analysis of MebZIPs reveals constitutive, tissue-specific and abiotic stress-responsive candidate MebZIP genes for further functional characterization in planta, yields new insights into transcriptional regulation of MebZIP genes, and lays a foundation for understanding of bZIP-mediated abiotic stress response.

Genome-wide gene phylogeny of CIPK family in cassava and expression analysis of partial drought-induced genes.

Cassava is an important food and potential biofuel crop that is tolerant to multiple abiotic stressors. The mechanisms underlying these tolerances are currently less known. CBL-interacting protein kinases (CIPKs) have been shown to play crucial roles in plant developmental processes, hormone signaling transduction, and in the response to abiotic stress. However, no data is currently available about the CPK family in cassava. In this study, a total of 25 CIPK genes were identified from cassava genome based on our previous genome sequencing data. Phylogenetic analysis suggested that 25 MeCIPKs could be classified into four subfamilies, which was supported by exon-intron organizations and the architectures of conserved protein motifs. Transcriptomic analysis of a wild subspecies and two cultivated varieties showed that most MeCIPKs had different expression patterns between wild subspecies and cultivatars in different tissues or in response to drought stress. Some orthologous genes involved in CIPK interaction networks were identified between Arabidopsis and cassava. The interaction networks and co-expression patterns of these orthologous genes revealed that the crucial pathways controlled by CIPK networks may be involved in the differential response to drought stress in different accessions of cassava. Nine MeCIPK genes were selected to investigate their transcriptional response to various stimuli and the results showed the comprehensive response of the tested MeCIPK genes to osmotic, salt, cold, oxidative stressors, and ABA signaling. The identification and expression analysis of CIPK family suggested that CIPK genes are important components of development and multiple signal transduction pathways in cassava. The findings of this study will help lay a foundation for the functional characterization of the CIPK gene family and provide an improved understanding of abiotic stress responses and signaling transduction in cassava.

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.

Genome-Wide Identification and Expression Analysis of the NAC Transcription Factor Family in Cassava

NAC [no apical meristem (NAM), Arabidopsis transcription activation factor [ATAF1/2] and cup-shaped cotyledon (CUC2)] proteins is one of the largest groups of plant specific transcription factors and plays a crucial role in plant growth, development, and adaption to the environment. Currently, no information is known about the NAC family in cassava. In this study, 96 NAC genes (MeNACs) were identified from the cassava genome. Phylogenetic analysis of the NACs from cassava and Arabidopsis showed that MeNAC proteins can be clustered into 16 subgroups. Gene structure analysis found that the number of introns of MeNAC genes varied from 0 to 5, with the majority of MeNAC genes containing two introns, indicating a small gene structure diversity of cassava NAC genes. Conserved motif analysis revealed that all of the identified MeNACs had the conserved NAC domain and/or NAM domain. Global expression analysis suggested that MeNAC genes exhibited different expression profiles in different tissues between wild subspecies and cultivated varieties, indicating their involvement in the functional diversity of different accessions. Transcriptome analysis demonstrated that MeNACs had a widely transcriptional response to drought stress and that they had differential expression profiles in different accessions, implying their contribution to drought stress resistance in cassava. Finally, the expression of twelve MeNAC genes was analyzed under osmotic, salt, cold, ABA, and H2O2 treatments, indicating that cassava NACs may represent convergence points of different signaling pathways. Taken together, this work found some excellent tissue-specific and abiotic stress-responsive candidate MeNAC genes, which would provide a solid foundation for functional investigation of the NAC family, crop improvement and improved understanding of signal transduction in plants. These data bring new insight on the complexity of the transcriptional control of MeNAC genes and support the hypothesis that NACs play an important role in plant growth, development, and adaption of 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.

Comparative Physiological and Transcriptomic Analyses Reveal the Actions of Melatonin in the Delay of Postharvest Physiological Deterioration of Cassava

Melatonin plays important roles in various aspects of biological processes. However, it is less known on the effects and mechanism of melatonin on the postharvest physiological deterioration (PPD) process of cassava, which largely restricts the potential of cassava as a food and industrial crop. In this study, we found that exogenous application of melatonin significantly delayed PPD of cassava tuberous roots by reducing H2O2 content and improving activities of catalase and peroxidase. Moreover, 3425 differentially expressed genes by melatonin during the PPD process were identified by transcriptomic analysis. Several pathways were markedly affected by melatonin treatments, including metabolic-, ion homeostasis-, and enzyme activity-related processes. Further detailed analysis revealed that melatonin acted through activation of ROS-scavenging and ROS signal transduction pathways, including antioxidant enzymes, calcium signaling, MAPK cascades, and transcription factors at early stages. Notably, the starch degradation pathway was also activated at early stages, whereas it was repressed by melatonin at middle and late stages, thereby indicating its regulatory role in starch metabolism during PPD. Taken together, this study yields new insights into the effect and underlying mechanism of melatonin on the delay of PPD and provides a good strategy for extending shelf life and improvement of cassava tuberous roots.

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.