A population based expression atlas provides insights into disease resistance and other physiological traits in cassava (Manihot esculenta Crantz)

A data-driven genome annotation approach for cassava

Genome annotation files play a critical role in dictating the quality of downstream analyses by providing essential predictions for gene positions and structures. These files are pivotal in decoding the complex information encoded within DNA sequences. Here, we generated experimental data resolving RNA 5'- and 3'-ends as well as full-length RNAs for cassava TME12 sticklings in ambient temperature and cold. We used these data to generate genome annotation files using the TranscriptomeReconstructoR (TR) tool. A careful comparison to high-quality genome annotations suggests that our new TR genome annotations identified additional genes, resolved the transcript boundaries more accurately and identified additional RNA isoforms. We enhanced existing cassava genome annotation files with the information from TR that maintained the different transcript models as RNA isoforms. The resultant merged annotation was subsequently utilized for comprehensive analysis. To examine the effects of genome annotation files on gene expression studies, we compared the detection of differentially expressed genes during cold using the same RNA-seq data but alternative genome annotation files. We found that our merged genome annotation that included cold-specific TR gene models identified about twice as many cold-induced genes. These data indicate that environmentally induced genes may be missing in off-the-shelf genome annotation files. In conclusion, TR offers the opportunity to enhance crop genome annotations with implications for the discovery of differentially expressed candidate genes during plant-environment interactions.

Identification of MeC3HDZ1/MeCNA as a potential regulator of cassava storage root development

The storage root (SR) of cassava is the main staple food in sub-Saharan Africa, where it feeds over 500 million people. However, little is known about the genetic and molecular regulation underlying its development. Unraveling such regulation would pave the way for biotechnology approaches aimed at enhancing cassava productivity. Anatomical studies indicate that SR development relies on the massive accumulation of xylem parenchyma, a cell-type derived from the vascular cambium. The C3HDZ family of transcription factors regulate cambial cells proliferation and xylem differentiation in Arabidopsis and other species. We thus aimed at identifying C3HDZ proteins in cassava and determining whether any of them shows preferential activity in the SR cambium and/or xylem. Using phylogeny and synteny studies, we identified eight C3HDZ proteins in cassava, namely MeCH3DZ1-8. We observed that MeC3HDZ1 is the MeC3HDZ gene displaying the highest expression in SR and that, within that organ, the gene also shows high expression in cambium and xylem. In-silico analyses revealed the existence of a number of potential C3HDZ targets displaying significant preferential expression in the SR. Subsequent Y1H analyses proved that MeC3HDZ1 can bind canonical C3HDZ binding sites, present in the promoters of these targets. Transactivation assays demonstrated that MeC3HDZ1 can regulate the expression of genes downstream of promoters harboring such binding sites, thereby demonstrating that MeC3HDZ1 has C3HDZ transcription factor activity. We conclude that MeC3HDZ1 may be a key factor for the regulation of storage root development in cassava, holding thus great promise for future biotechnology applications.