Gene expression atlas for the food security crop cassava

A foundational large language model for edible plant genomes

Significant progress has been made in the field of plant genomics, as demonstrated by the increased use of high-throughput methodologies that enable the characterization of multiple genome-wide molecular phenotypes. These findings have provided valuable insights into plant traits and their underlying genetic mechanisms, particularly in model plant species. Nonetheless, effectively leveraging them to make accurate predictions represents a critical step in crop genomic improvement. We present AgroNT, a foundational large language model trained on genomes from 48 plant species with a predominant focus on crop species. We show that AgroNT can obtain state-of-the-art predictions for regulatory annotations, promoter/terminator strength, tissue-specific gene expression, and prioritize functional variants. We conduct a large-scale in silico saturation mutagenesis analysis on cassava to evaluate the regulatory impact of over 10 million mutations and provide their predicted effects as a resource for variant characterization. Finally, we propose the use of the diverse datasets compiled here as the Plants Genomic Benchmark (PGB), providing a comprehensive benchmark for deep learning-based methods in plant genomic research. The pre-trained AgroNT model is publicly available on HuggingFace at https://huggingface.co/InstaDeepAI/agro-nucleotide-transformer-1b  for future research purposes.

Carbon usage in yellow-fleshed Manihot esculenta storage roots shifts from starch biosynthesis to cell wall and raffinose biosynthesis via the myo-inositol pathway

Cassava is a crucial staple crop for smallholder farmers in tropical Asia and Sub-Saharan Africa. Although high yield remains the top priority for farmers, the significance of nutritional values has increased in cassava breeding programs. A notable negative correlation between provitamin A and starch accumulation poses a significant challenge for breeding efforts. The negative correlation between starch and carotenoid levels in conventional and genetically modified cassava plants implies the absence of a direct genomic connection between the two traits. The competition among various carbon pathways seems to account for this relationship. In this study, we conducted a thorough analysis of 49 African cassava genotypes with varying levels of starch and provitamin A. Our goal was to identify factors contributing to differential starch accumulation. Considering carotenoid levels as a confounding factor in starch production, we found that yellow- and white-fleshed storage roots did not differ significantly in most measured components of starch or de novo fatty acid biosynthesis. However, genes and metabolites associated with myo-inositol synthesis and cell wall polymer production were substantially enriched in high provitamin A genotypes. These results indicate that yellow-fleshed cultivars, in comparison to their white-fleshed counterparts, direct more carbon toward the synthesis of raffinose and cell wall components. This finding is underlined by a significant rise in cell wall components measured within the 20 most contrasting genotypes for carotenoid levels. Our findings enhance the comprehension of the biosynthesis of starch and carotenoids in the storage roots of cassava.

MeGT2.6 increases cellulose synthesis and active gibberellin content to promote cell enlargement in cassava

Cassava, a pivotal tropical crop, exhibits rapid growth and possesses a substantial biomass. Its stem is rich in cellulose and serves as a crucial carbohydrate storage organ. The height and strength of stems restrict the mechanised operation and propagation of cassava. In this study, the triple helix transcription factor MeGT2.6 was identified through yeast one-hybrid assay using MeCesA1pro as bait, which is critical for cellulose synthesis. Over-expression and loss-of-function lines were generated, and results revealed that MeGT2.6 could promote a significant increase in the plant height, stem diameter, cell size and thickness of SCW of cassava plant. Specifically, MeGT2.6 upregulated the transcription activity of MeGA20ox1 and downregulated the expression level of MeGA2ox1, thereby enhancing the content of active GA3, resulting in a large cell size, high plant height and long stem diameter in cassava. Moreover, MeGT2.6 upregulated the transcription activity of MeCesA1, which promoted the synthesis of cellulose and hemicellulose and produced a thick secondary cell wall. Finally, MeGT2.6 could help supply additional substrates for the synthesis of cellulose and hemicellulose by upregulating the invertase genes (MeNINV1/6). Thus, MeGT2.6 was found to be a multiple regulator; it was involved in GA metabolism and sucrose decomposition and the synthesis of cellulose and hemicellulose.

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.

Genome-wide identification and expression analysis of heat shock protein gene family in cassava

Heat shock proteins are important molecular chaperones that are involved in plant growth and stress responses. However, members of the Hsp family have been poorly studied in cassava. In this study, 225 MeHsp genes were identified in the cassava genome, and their genetic structures exhibited relatively conserved features within each subfamily. The 225 MeHsp genes showed random chromosomal distribution, and at least 74 pairs of segmentally duplicated MeHsp genes. Eleven tandemly duplicated MeHsp genes were identified. Cis-element analysis revealed the importance of MeHsps in plant adaptations to the environment. The prediction of protein interactions suggested that MeHsp70-20 may play a critical regulatory role in the interactive network. Furthermore, the expression profiles of MeHsps in different tissues and cell subsets were analyzed using bulk transcriptomics and single-cell transcriptomic data. Several subfamily genes exhibited unique expression patterns in the transcriptome and were selected for detailed analysis of the single-cell transcriptome. Quantitative real-time polymerase chain reaction (qRT-PCR) revealed the expression patterns of these genes under temperature stress, further supporting the prediction of cis-acting elements. This study provides valuable information for understanding the functional characteristics of MeHsp genes and the evolutionary relationships between MeHsps.

LESION SIMULATING DISEASE 3 regulates disease resistance via fine-tuning histone acetylation in cassava

Bacterial blight seriously affects the growth and production of cassava (Manihot esculenta Crantz), but disease resistance genes and the underlying molecular mechanism remain unknown. In this study, we found that LESION SIMULATING DISEASE 3 (MeLSD3) is essential for disease resistance in cassava. MeLSD3 physically interacts with SIRTUIN 1 (MeSRT1), inhibiting MeSRT1-mediated deacetylation modification at the acetylation of histone 3 at K9 (H3K9Ac). This leads to increased H3K9Ac levels and transcriptional activation of SUPPRESSOR OF BIR1 (SOBIR1) and FLAGELLIN-SENSITIVE2 (FLS2) in pattern-triggered immunity, resulting in immune responses in cassava. When MeLSD3 was silenced, the release of MeSRT1 directly decreased H3K9Ac levels and inhibited the transcription of SOBIR1 and FLS2, leading to decreased disease resistance. Notably, DELLA protein GIBBERELLIC ACID INSENSITIVE 1 (MeGAI1) also interacted with MeLSD3, which enhanced the interaction between MeLSD3 and MeSRT1 and further strengthened the inhibition of MeSRT1-mediated deacetylation modification at H3K9Ac of defense genes. In summary, this study illustrates the mechanism by which MeLSD3 interacts with MeSRT1 and MeGAI1, thereby mediating the level of H3K9Ac and the transcription of defense genes and immune responses in cassava.

Genome-wide identification and expression analysis of Hsp70 family genes in Cassava (Manihot esculenta Crantz)

Hsp70 proteins function as molecular chaperones, regulating various cellular processes in plants. In this study, a genome-wide analysis led to the identification of 22 Hsp70 (MeHsp70) genes in cassava. Phylogenetic relationship studies with other Malpighiales genomes (Populus trichocarpa, Ricinus communis and Salix purpurea) classified MeHsp70 proteins into eight groups (Ia, Ib, Ic, Id, Ie, If, IIa and IIb). Promoter analysis of MeHsp70 genes revealed the presence of tissue-specific, light, biotic and abiotic stress-responsive cis-regulatory elements showing their functional importance in cassava. Meta-analysis of publically available RNA-seq transcriptome datasets showed constitutive, tissue-specific, biotic and abiotic stress-specific expression patterns among MeHsp70s in cassava. Among 22 Hsp70, six MeHsp70s viz., MecHsp70-3, MecHsp70-6, MeBiP-1, MeBiP-2, MeBiP-3 and MecpHsp70-2 displayed constitutive expression, while three MecHsp70s were induced under both drought and cold stress conditions. Five MeHsp70s, MecHsp70-7, MecHsp70-11, MecHsp70-12, MecHsp70-13, and MecHsp70-14 were induced under drought stress conditions. We predicted that 19 MeHsp70 genes are under the regulation of 24 miRNAs. This comprehensive genome-wide analysis of the Hsp70 gene family in cassava provided valuable insights into their functional roles and identified various potential Hsp70 genes associated with stress tolerance and adaptation to environmental stimuli.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03760-3.

Genome-wide identification and characterization of 14-3-3 gene family related to negative regulation of starch accumulation in storage root of Manihot esculenta

The 14-3-3 protein family is a highly conservative member of the acid protein family and plays an important role in regulating a series of important biological activities and various signal transduction pathways. The role of 14-3-3 proteins in regulating starch accumulation still remains largely unknown. To investigate the properties of 14-3-3 proteins, the structures and functions involved in starch accumulation in storage roots were analyzed, and consequently, 16 Me14-3-3 genes were identified. Phylogenetic analysis revealed that Me14-3-3 family proteins are split into two groups (ε and non-ε). All Me14-3-3 proteins contain nine antiparallel α-helices. Me14-3-3s-GFP fusion protein was targeted exclusively to the nuclei and cytoplasm. In the early stage of starch accumulation in the storage root, Me14-3-3 genes were highly expressed in high-starch cultivars, while in the late stage of starch accumulation, Me14-3-3 genes were highly expressed in low-starch cultivars. Me14-3-3 I, II, V, and XVI had relatively high expression levels in the storage roots. The transgenic evidence from Me14-3-3II overexpression in Arabidopsis thaliana and the virus-induced gene silencing (VIGS) in cassava leaves and storage roots suggest that Me14-3-3II is involved in the negative regulation of starch accumulation. This study provides a new insight to understand the molecular mechanisms of starch accumulation linked with Me14-3-3 genes during cassava storage root development.

Integrated Characterization of Cassava (Manihot esculenta) Pectin Methylesterase (MePME) Genes to Filter Candidate Gene Responses to Multiple Abiotic Stresses

Plant pectin methylesterases (PMEs) play crucial roles in regulating cell wall modification and response to various stresses. Members of the PME family have been found in several crops, but there is a lack of research into their presence in cassava (Manihot esculent), which is an important crop for world food security. In this research, 89 MePME genes were identified in cassava that were separated into two types (type-Ⅰ and type-Ⅱ) according to the existence or absence of a pro-region (PMEI domain). The MePME gene members were unevenly located on 17 chromosomes, with 19 gene pairs being identified that most likely arose via duplication events. The MePMEs could be divided into ten sub-groups in type-Ⅰ and five sub-groups in type-Ⅱ. The motif analysis revealed 11 conserved motifs in type-Ⅰ and 8 in type-Ⅱ MePMEs. The number of introns in the CDS region of type-Ⅰ MePMEs ranged between one and two, and the number of introns in type-Ⅱ MePMEs ranged between one and nine. There were 21 type-Ⅰ and 31 type-Ⅱ MePMEs that contained signal peptides. Most of the type-Ⅰ MePMEs had two conserved "RK/RLL" and one "FPSWVS" domain between the pro-region and the PME domain. Multiple stress-, hormone- and tissue-specific-related cis-acting regulatory elements were identified in the promoter regions of MePME genes. A total of five co-expressed genes (MePME1, MePME2, MePME27, MePME65 and MePME82) were filtered from different abiotic stresses via the use of UpSet Venn diagrams. The gene expression pattern analysis revealed that the expression of MePME1 was positively correlated with the degree of cassava postharvest physiological deterioration (PPD). The expression of this gene was also significantly upregulated by 7% PEG and 14 °C low-temperature stress, but slightly downregulated by ABA treatment. The tissue-specific expression analysis revealed that MePME1 and MePME65 generally displayed higher expression levels in most tissues than the other co-expressed genes. In this study, we obtain an in-depth understanding of the cassava PME gene family, suggesting that MePME1 could be a candidate gene associated with multiple abiotic tolerance.

Magnesium chelatase subunit D is not only required for chlorophyll biosynthesis and photosynthesis, but also affecting starch accumulation in Manihot esculenta Crantz

BACKGROUND

Magnesium chelatase plays an important role in photosynthesis, but only a few subunits have been functionally characterized in cassava.

RESULTS

Herein, MeChlD was successfully cloned and characterized. MeChlD encodes a magnesium chelatase subunit D, which has ATPase and vWA conservative domains. MeChlD was highly expressed in the leaves. Subcellular localization suggested that MeChlD:GFP was a chloroplast-localized protein. Furthermore, the yeast two-hybrid system and BiFC analysis indicated that MeChlD interacts with MeChlM and MePrxQ, respectively. VIGS-induce silencing of MeChlD resulted in significantly decreased chlorophyll content and reduction the expression of photosynthesis-related nuclear genes. Furthermore, the storage root numbers, fresh weight and the total starch content in cassava storage roots of VIGS-MeChlD plants was significantly reduced.

CONCLUSION

Taken together, MeChlD located at the chloroplast is not only required for chlorophyll biosynthesis and photosynthesis, but also affecting the starch accumulation in cassava. This study expands our understanding of the biological functions of ChlD proteins.

Genome-wide analysis of the metallothionein gene family in cassava reveals its role in response to physiological stress through the regulation of reactive oxygen species

BACKGROUND

Cassava (Manihot esculenta Crantz) is widely planted in tropical and several subtropical regions in which drought, high temperatures, and other abiotic stresses occur. Metallothionein (MT) is a group of conjugated proteins with small molecular weight and rich in cysteine. These proteins play a substantial role in response to physiological stress through the regulation of reactive oxygen species (ROS). However, the biological functions of MT genes in cassava are unknown.

RESULTS

A total of 10 MeMT genes were identified in the cassava genome. The MeMTs were divided into 3 groups (Types 2-4) based on the contents and distribution of Cys residues. The MeMTs exhibited tissue-specific expression and located on 7 chromosomes. The MeMT promoters contain some hormones regulatory and stresses responsiveness elements. MeMTs were upregulated under hydrogen peroxide (H₂O₂) treatment and in respond to post-harvest physiological deterioration (PPD). The results were consistent with defense-responsive cis-acting elements in the MeMT promoters. Further, four of MeMTs were selected and silenced by using the virus-induced gene silencing (VIGS) method to evaluate their functional characterization. The results of gene-silenced cassava suggest that MeMTs are involved in oxidative stress resistance, as ROS scavengers.

CONCLUSION

We identified the 10 MeMT genes, and explore their evolutionary relationship, conserved motif, and tissue-specific expression. The expression profiles of MeMTs under three kinds of abiotic stresses (wounding, low-temperature, and H₂O₂) and during PPD were analyzed. The tissue-specific expression and the response to abiotic stresses revealed the role of MT in plant growth and development. Furthermore, silenced expression of MeMTs in cassava leaves decreased its tolerance to ROS, consistent with its predicted role as ROS scavengers. In summary, our results suggest an important role of MeMTs in response to physiological stress as well as species adaptation via the regulation of ROS homeostasis.

Grapevine scion gene expression is driven by rootstock and environment interaction

BACKGROUND

Grafting is a horticultural practice used widely across woody perennial crop species to fuse together the root and shoot system of two distinct genotypes, the rootstock and the scion, combining beneficial traits from both. In grapevine, grafting is used in nearly 80% of all commercial vines to optimize fruit quality, regulate vine vigor, and enhance biotic and abiotic stress-tolerance. Rootstocks have been shown to modulate elemental composition, metabolomic profiles, and the shape of leaves in the scion, among other traits. However, it is currently unclear how rootstock genotypes influence shoot system gene expression as previous work has reported complex and often contradictory findings.

RESULTS

In the present study, we examine the influence of grafting on scion gene expression in leaves and reproductive tissues of grapevines growing under field conditions for three years. We show that the influence from the rootstock genotype is highly tissue and time dependent, manifesting only in leaves, primarily during a single year of our three-year study. Further, the degree of rootstock influence on scion gene expression is driven by interactions with the local environment.

CONCLUSIONS

Our results demonstrate that the role of rootstock genotype in modulating scion gene expression is not a consistent, unchanging effect, but rather an effect that varies over time in relation to local environmental conditions.

Genome-Wide Identification and Expression Analysis of the SHI-Related Sequence Family in Cassava

The SHORT INTERNODES (SHI)-related sequences (SRS) are plant-specific transcription factors that have been quantitatively characterized during plant growth, regeneration, and stress responses. However, the genome-wide discovery of SRS family genes and their involvement in abiotic stress-related activities in cassava have not been documented. A genome-wide search strategy was used to identify eight family members of the SRS gene family in cassava (Manihot esculenta Crantz). Based on their evolutionary linkages, all MeSRS genes featured homologous RING-like zinc finger and IXGH domains. Genetic architecture and conserved motif analysis validated the categorization of MeSRS genes into four groups. Eight pairs of segmental duplications were detected, resulting in an increase in the number of MeSRS genes. Orthologous studies of SRS genes among cassava and three different plant species (Arabidopsis thaliana, Oryza sativa, and Populus trichocarpa) provided important insights into the probable history of the MeSRS gene family. The functionality of MeSRS genes was elucidated through the prediction of protein–protein interaction networks and cis-acting domains. RNA-seq data demonstrated tissue/organ expression selectivity and preference of the MeSRS genes. Furthermore, qRT-PCR investigation of MeSRS gene expression after exposure to salicylic acid (SA) and methyl jasmonate (MeJA) hormone treatments, as well as salt (NaCl) and osmotic (polyethylene glycol, PEG) stresses, showed their stress-responsive patterns. This genome-wide characterization and identification of the evolutionary relationships and expression profiles of the cassava MeSRS family genes will be helpful for further research into this gene family and its function in stress response. It may also assist future agricultural efforts to increase the stress tolerance of cassava.

CRISPR-Cas9-mediated knockout of CYP79D1 and CYP79D2 in cassava attenuates toxic cyanogen production

Cassava (Manihot esculenta) is a starchy root crop that supports over a billion people in tropical and subtropical regions of the world. This staple, however, produces the neurotoxin cyanide and requires processing for safe consumption. Excessive consumption of insufficiently processed cassava, in combination with protein-poor diets, can have neurodegenerative impacts. This problem is further exacerbated by drought conditions which increase this toxin in the plant. To reduce cyanide levels in cassava, we used CRISPR-mediated mutagenesis to disrupt the cytochrome P450 genes CYP79D1 and CYP79D2 whose protein products catalyze the first step in cyanogenic glucoside biosynthesis. Knockout of both genes eliminated cyanide in leaves and storage roots of cassava accession 60444; the West African, farmer-preferred cultivar TME 419; and the improved variety TMS 91/02324. Although knockout of CYP79D2 alone resulted in significant reduction of cyanide, mutagenesis of CYP79D1 did not, indicating these paralogs have diverged in their function. The congruence of results across accessions indicates that our approach could readily be extended to other preferred or improved cultivars. This work demonstrates cassava genome editing for enhanced food safety and reduced processing burden, against the backdrop of a changing climate.

Reactive oxygen species turnover, phenolics metabolism, and some key gene expressions modulate postharvest physiological deterioration in cassava tubers

Rapid postharvest physiological deterioration (PPD) in cassava (Manihot esculenta Crantz) tuber is a significant concern during storage. The freshly harvested tubers start spoiling within 24 to 72 h. Accumulation of H2O2 is one of the earliest biochemical events that occurred during PPD, which was detected using the 3,3 diaminobenzidine (DAB) in two contrast cassava genotypes, MNP Local A (29-57 μg g-1) and Sree Prakash (64-141 μg g-1). Accumulating the fluorescence hydroxycoumarin compounds emitted by the cassava tubers observed under an ultraviolet (UV) lamp showed significant variations at 0, 3, 6, 9, 12, and 15 days of storage. The total phenolics and carotenoids significantly and negatively correlated with PPD progression; however, the anthocyanin and flavonoids positively correlated with the PPD-anchored ROS accumulation. The primary compound, Phthalic acid, di(2-propylpentyl) ester, was identified in both the cassava tubers, Sree Prakash (57.21 and 35.21%), and MNP Local A (75.58 and 60.21%) at 0, and 72 h of PPD, respectively. The expression of PPD-associated genes APX-2, APX-3, PAL, and AP was higher at 6-12 days of PPD, which signified the synthesis of ROS turnover and phenylpropanoid biosynthesis. A significant, strong, and positive correlation was established between the secondary metabolites and PPD signaling gene expression, which was inversely correlated with hydroxycoumarin and H2O2 accumulation. MNP Local A tubers exhibited longer storage life of 15 days with a low PPD score, higher metabolites synthesis, and gene expression. The PPD-resistant lines may be used to augment cassava breeding strategies for large-scale commercial and industrial use.

Toward kingdom-wide analyses of gene expression

Gene expression data for Archaeplastida are accumulating exponentially, with more than 300 000 RNA-sequencing (RNA-seq) experiments available for hundreds of species. The gene expression data stem from thousands of experiments that capture gene expression in various organs, tissues, cell types, (a)biotic perturbations, and genotypes. Advances in software tools make it possible to process all these data in a matter of weeks on modern office computers, giving us the possibility to study gene expression in a kingdom-wide manner for the first time. We discuss how the expression data can be accessed and processed and outline analyses that take advantage of cross-species analyses, allowing us to generate powerful and robust hypotheses about gene function and evolution.

Genome-wide identification, expression profiling, and functional analysis of ammonium transporter 2 (AMT2) gene family in cassava (Manihot esculenta crantz)

Background: Nitrogen (N), absorbed primarily as ammonium (NH₄ ⁺) from soil by plant, is a necessary macronutrient in plant growth and development. Ammonium transporter (AMT) plays a vital role in the absorption and transport of ammonium (NH₄ ⁺). Cassava (Manihot esculenta Crantz) has a strong adaptability to nitrogen deprivation. However, little is known about the functions of ammonium transporter AMT2 in cassava. Methods: The cassava AMT2-type genes were identified and their characteristics were analyzed using bioinformatic techniques. The spatial expression patterns were analyzed based on the public RNA-seq data and their expression profiles under low ammonium treatment were studied using Real-time quantitative PCR (RT-qPCR) method. The cassava AMT2 genes were transformed into yeast mutant strain TM31019b by PEG/LiAc method to investigate their functions. Results: Seven AMT2-type genes (MeAMT2.1-2.7) were identified in cassava and they were distributed on 6 chromosomes and included two segmental duplication events (MeAMT2.2/MeAMT2.4 and MeAMT2.3/MeAMT2.5). Based on their amino acid sequences, seven MeAMT2 were further divided into four subgroups, and each subgroup contained similar motif constitution and protein structure. Synteny analysis showed that two and four MeAMT2 genes in cassava were collinear with those in the Arabidopsis and soybean genomes, respectively. Sixteen types of cis-elements were identified in the MeAMT2 promoters, and they were related to light-, hormone-, stress-, and plant growth and development-responsive elements, respectively. Most of the MeAMT2 genes displayed tissue-specific expression patterns according to the RNA-seq data, of them, three MeAMT2 (MeAMT2.3, MeAMT2.5, and MeATM2.6) expressions were up-regulated under ammonium deficiency. Complementation experiments showed that yeast mutant strain TM31019b transformed with MeAMT2.3, MeAMT2.5, or MeATM2.6 grew better than untransgenic yeast cells under ammonium deficiency, suggesting that MeAMT2.3, MeAMT2.5, and MeATM2.6 might be the main contributors in response to ammonium deficiency in cassava. Conclusion: This study provides a basis for further study of nitrogen efficient utilization in cassava.

Comparative Transcriptome Profiling of Cassava Tuberous Roots in Response to Postharvest Physiological Deterioration

Cassava is one of the most versatile tuberous-root crops on Earth. However, the postharvest storage properties of cassava tuberous root mean that it is perishable through a process known as postharvest physiological deterioration (PPD), which seriously affects its starch quality. Therefore, a comprehensive understanding of the transcriptional regulatory activity of cassava against the PPD response is necessary in order to extract key molecular mechanisms related to PPD tolerance. In this study, we found that RYG1 tuberous roots showed delayed PPD compared to those of SC8. In addition, RYG1 roots maintained a more stable cell wall structure after storage than those of SC8. The transcriptome changes in tuberous roots were analyzed for both RYG1 and SC8 after 21 days of storage (SR and SS) compared to fresh (FR and FS) by the RNA-Seq method. The total number of differentially expressed genes (DEGs) in the various comparisons of these four samples ranged from 68 to 3847. Of these, a total of 2008 co-DEGs in SR vs. SS were shared by either SR vs. FR or SS vs. FS. GO and KEGG enrichment analysis revealed that upregulated co-DEGs in SR vs. SS were mainly enriched in photosynthesis, protein processing, hormone and cutin, suberine and wax biosynthesis. By contrast, the downregulated co-DEGs were mainly related to cell wall organization, starch and sucrose metabolism, galactose metabolism, phenylpropanoid biosynthesis, diterpenoid biosynthesis, cysteine and methionine metabolism and flavonoid biosynthesis. The protein-protein interaction (PPI) networks of the co-DEGs showed a complex interaction of genes in different pathways, and 16 hub genes were characterized to have a degree in excess of 15, among which eight genes were associated with photosynthesis. These results provide new information for the study of cassava resistance to PPD and lay a foundation for the further molecular breeding of storage-tolerant cassava varieties.

A promoter toolbox for tissue-specific expression supporting translational research in cassava (Manihot esculenta)

There is an urgent need to stimulate agricultural output in many tropical and subtropical countries of the world to combat hunger and malnutrition. The starchy crop cassava (Manihot esculenta), growing even under sub-optimal conditions, is a key staple food in these regions, providing millions of people with food. Cassava biotechnology is an important technique benefiting agricultural progress, but successful implementation of many biotechnological concepts depends on the availability of the right spatiotemporal expression tools. Yet, well-characterized cassava promoters are scarce in the public domain. In this study, we investigate the promoter activity and tissue specificity of 24 different promoter elements in stably transformed cassava plants. We show that many of the investigated promoters, especially from other species, have surprisingly low activity and/or tissue specificity, but feature several promoter sequences that can drive tissue-specific expression in either autotrophic-, transport- or storage tissues. We especially highlight pAtCAB1, pMePsbR, and pSlRBCS2 as strong and specific source promoters, pAtSUC2, pMeSWEET1-like, and pMeSUS1 as valuable tools for phloem and phloem parenchyma expression, and pStB33, pMeGPT, pStGBSS1, as well as pStPatatin Class I, as strong and specific promoters for heterotrophic storage tissues. We hope that the provided information and sequences prove valuable to the cassava community by contributing to the successful implementation of biotechnological concepts aimed at the improvement of cassava nutritional value and productivity.

The Effect of Climate Change on Food Crop Production in Lagos State

Climate change is set to be particularly disruptive in poor agricultural communities. This study examines the effects of, and farmer's perceptions of, climate change on farming practices for cassava and maize in Lagos, Nigeria. Analysis of weather data from 1998 to 2018 (the most recent available) reveals little impact on cassava yield but a significant impact on maize yield. Furthermore, survey results indicate that farmers in this area are currently implementing techniques to adapt to changes in climate based on the type of crop grown. Agriculture in Lagos, Nigeria, is largely rain-fed and climate change negatively impacts crop productivity by decreasing crop yield and soil fertility, limiting the availability of soil water, increasing soil erosion, and contributing to the spread of pests. A decline in crop production due to climate change may be further exasperated by a lack of access to farming technology that reduces over-reliance on the rain-fed farming system and subsistence agriculture. This study indicates that there is a need for initiatives to motivate young and older farmers through access to credits, irrigation facilities, and innovative climate change adaptive strategies.

Genome-wide identification and expression of TIFY family in cassava (Manihot esculenta Crantz)

Plant-specific TIFY [TIF(F/Y)XG] proteins serve important roles in the regulation of plant stress responses. This family encodes four subfamilies of proteins, JAZ (JASMONATE ZIM-domain), PPD (PEAPOD), ZML (Zinc-finger Inflorescence-like), and TIFY. In this work, a total of 16 JAZ, 3 PPD, 7 ZML, and 2 TIFY genes were found in cassava (Manihot esculenta Crantz) at the genome-wide level. The phylogenetics, exon-intron structure, motif organization, and conserved domains of these genes were analyzed to characterize the members of the JAZ, PPD, and ZML subfamilies. Chromosome location and synteny analyses revealed that 26 JAZ, PPD, and ZML genes were irregularly distributed across 14 of the 18 chromosomes, and 18 gene pairs were implicated in large-scale interchromosomal segmental duplication events. In addition, JAZ, PPD, and ZML gene synteny comparisons between cassava and three other plant species (Arabidopsis, Populus trichocarpa, and rice) uncovered vital information about their likely evolution. The prediction of protein interaction network and cis-acting elements reveal the function of JAZ, PPD, and ZML genes. Subsequently, expression patterns of JAZ, PPD, and ZML genes were validated by qRT-PCR as being expressed in response to osmotic, salt, and cadmium stress. Moreover, almost all JAZ subfamily genes were responsive to jasmonic acid (JA) treatment. In particular, MeJAZ1, MeJAZ13, and MeJAZ14, were highly up-regulated by three treatments, and these genes may deserve further study. This comprehensive study lays the groundwork for future research into TIFY family genes in cassava and may be valuable for genetic improvement of cassava and other related species.

Mutations in DNA polymerase δ subunit 1 co-segregate with CMD2-type resistance to Cassava Mosaic Geminiviruses

Cassava mosaic disease (CMD) suppresses cassava yields across the tropics. The dominant CMD2 locus confers resistance to cassava mosaic geminiviruses. It has been reported that CMD2-type landraces lose resistance after regeneration through de novo morphogenesis. As full genome bisulfite sequencing failed to uncover an epigenetic mechanism for this loss of resistance, whole genome sequencing and genetic variant analysis was performed and the CMD2 locus was fine-mapped to a 190 kilobase interval. Collectively, these data indicate that CMD2-type resistance is caused by a nonsynonymous, single nucleotide polymorphism in DNA polymerase δ subunit 1 (MePOLD1) located within this region. Virus-induced gene silencing of MePOLD1 in a CMD-susceptible cassava variety produced a recovery phenotype typical of CMD2-type resistance. Analysis of other CMD2-type cassava varieties identified additional candidate resistance alleles within MePOLD1. Genetic variation of MePOLD1, therefore, could represent an important genetic resource for resistance breeding and/or genome editing, and elucidating mechanisms of resistance to geminiviruses.

Current status and impending progress for cassava structural genomics

We demystify recent advances in genome assemblies for the heterozygous staple crop cassava (Manihot esculenta), and highlight key cassava genomic resources. Cassava, Manihot esculenta Crantz, is a crop of societal and agricultural importance in tropical regions around the world. Genomics provides a platform for accelerated improvement of cassava's nutritional and agronomic traits, as well as for illuminating aspects of cassava's history including its path towards domestication. The highly heterozygous nature of the cassava genome is widely recognized. However, the full extent and context of this heterozygosity has been difficult to reveal because of technological limitations within genome sequencing. Only recently, with several new long-read sequencing technologies coming online, has the genomics community been able to tackle some similarly difficult genomes. In light of these recent advances, we provide this review to document the current status of the cassava genome and genomic resources and provide a perspective on what to look forward to in the coming years.

Cassava shrunken-2 homolog MeAPL3 determines storage root starch and dry matter content and modulates storage root postharvest physiological deterioration

Among the five cassava isoforms (MeAPL1-MeAPL5), MeAPL3 is responsible for determining storage root starch content. Degree of storage root postharvest physiological deterioration (PPD) is directly correlated with starch content. AGPase is heterotetramer composed of two small and two large subunits each coded by small gene families in higher plants. Studies in cassava (Manihot esculenta) identified and characterized five isoforms of Manihot esculenta ADP-glucose pyrophosphorylase large subunit (MeAPL1-MeAPL5) and employed virus induced gene silencing (VIGS) to show that MeAPL3 is the key isoform responsible for starch and dry matter accumulation in cassava storage roots. Silencing of MeAPL3 in cassava through stable transgenic lines resulted in plants displaying significant reduction in storage root starch and dry matter content (DMC) and induced a distinct phenotype associated with increased petiole/stem angle, resulting in a droopy leaf phenotype. Plants with reduced starch and DMC also displayed significantly reduced or no postharvest physiological deterioration (PPD) compared to controls and lines with high DMC and starch content. This provides strong evidence for direct relationships between starch/dry matter content and its role in PPD and canopy architecture traits in cassava.

Genome-Wide Identification of Cassava Glyoxalase I Genes and the Potential Function of MeGLYⅠ-13 in Iron Toxicity Tolerance

Glyoxalase I (GLYI) is a key enzyme in the pathway of the glyoxalase system that degrades the toxic substance methylglyoxal, which plays a crucial part in plant growth, development, and stress response. A total of 19 GLYI genes were identified from the cassava genome, which distributed randomly on 11 chromosomes. These genes were named MeGLYI-1–19 and were systematically characterized. Transcriptome data analysis showed that MeGLYIs gene expression is tissue-specific, and MeGLYI-13 is the dominant gene expressed in young tissues, while MeGLYI-19 is the dominant gene expressed in mature tissues and organs. qRT-PCR analysis showed that MeGLYI-13 is upregulated under 2 h excess iron stress, but downregulated under 6, 12, and 20 h iron stress. Overexpression of MeGLYI-13 enhanced the growth ability of transgenic yeast under iron stress. The root growth of transgenic Arabidopsis seedlings was less inhibited by iron toxicity than that of the wild type (WT). Potted transgenic Arabidopsis blossomed and podded under iron stress, but flowering of the WT was significantly delayed. The GLYI activity in transgenic Arabidopsis was improved under both non-iron stress and iron stress conditions compared to the WT. The SOD activity in transgenic plants was increased under iron stress, while the POD and CAT activity and MDA content were decreased compared to that in the WT. These results provide a basis for the selection of candidate genes for iron toxicity tolerance in cassava, and lay a theoretical foundation for further studies on the functions of these MeGLYI genes.

The haplotype-resolved chromosome pairs of a heterozygous diploid African cassava cultivar reveal novel pan-genome and allele-specific transcriptome features

BACKGROUND

Cassava (Manihot esculenta) is an important clonally propagated food crop in tropical and subtropical regions worldwide. Genetic gain by molecular breeding has been limited, partially because cassava is a highly heterozygous crop with a repetitive and difficult-to-assemble genome.

FINDINGS

Here we demonstrate that Pacific Biosciences high-fidelity (HiFi) sequencing reads, in combination with the assembler hifiasm, produced genome assemblies at near complete haplotype resolution with higher continuity and accuracy compared to conventional long sequencing reads. We present 2 chromosome-scale haploid genomes phased with Hi-C technology for the diploid African cassava variety TME204. With consensus accuracy >QV46, contig N50 >18 Mb, BUSCO completeness of 99%, and 35k phased gene loci, it is the most accurate, continuous, complete, and haplotype-resolved cassava genome assembly so far. Ab initio gene prediction with RNA-seq data and Iso-Seq transcripts identified abundant novel gene loci, with enriched functionality related to chromatin organization, meristem development, and cell responses. During tissue development, differentially expressed transcripts of different haplotype origins were enriched for different functionality. In each tissue, 20-30% of transcripts showed allele-specific expression (ASE) differences. ASE bias was often tissue specific and inconsistent across different tissues. Direction-shifting was observed in <2% of the ASE transcripts. Despite high gene synteny, the HiFi genome assembly revealed extensive chromosome rearrangements and abundant intra-genomic and inter-genomic divergent sequences, with large structural variations mostly related to LTR retrotransposons. We use the reference-quality assemblies to build a cassava pan-genome and demonstrate its importance in representing the genetic diversity of cassava for downstream reference-guided omics analysis and breeding.

CONCLUSIONS

The phased and annotated chromosome pairs allow a systematic view of the heterozygous diploid genome organization in cassava with improved accuracy, completeness, and haplotype resolution. They will be a valuable resource for cassava breeding and research. Our study may also provide insights into developing cost-effective and efficient strategies for resolving complex genomes with high resolution, accuracy, and continuity.

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

Cassava, a food security crop in Africa, is grown throughout the tropics and subtropics. Although cassava can provide high productivity in suboptimal conditions, the yield in Africa is substantially lower than in other geographies. The yield gap is attributable to many challenges faced by cassava in Africa, including susceptibility to diseases and poor soil conditions. In this study, we carried out 3’RNA sequencing on 150 accessions from the National Crops Resources Research Institute, Uganda for 5 tissue types, providing population-based transcriptomics resources to the research community in a web-based queryable cassava expression atlas. Differential expression and weighted gene co-expression network analysis were performed to detect 8820 significantly differentially expressed genes (DEGs), revealing similarity in expression patterns between tissue types and the clustering of detected DEGs into 18 gene modules. As a confirmation of data quality, differential expression and pathway analysis targeting cassava mosaic disease (CMD) identified 27 genes observed in the plant–pathogen interaction pathway, several previously identified CMD resistance genes, and two peroxidase family proteins different from the CMD2 gene. Present research work represents a novel resource towards understanding complex traits at expression and molecular levels for the development of resistant and high-yielding cassava varieties, as exemplified with CMD.

Large structural variations in the haplotype-resolved African cassava genome

Cassava (Manihot esculenta Crantz, 2n = 36) is a global food security crop. It has a highly heterozygous genome, high genetic load, and genotype-dependent asynchronous flowering. It is typically propagated by stem cuttings and any genetic variation between haplotypes, including large structural variations, is preserved by such clonal propagation. Traditional genome assembly approaches generate a collapsed haplotype representation of the genome. In highly heterozygous plants, this results in artifacts and an oversimplification of heterozygous regions. We used a combination of Pacific Biosciences (PacBio), Illumina, and Hi-C to resolve each haplotype of the genome of a farmer-preferred cassava line, TME7 (Oko-iyawo). PacBio reads were assembled using the FALCON suite. Phase switch errors were corrected using FALCON-Phase and Hi-C read data. The ultralong-range information from Hi-C sequencing was also used for scaffolding. Comparison of the two phases revealed >5000 large haplotype-specific structural variants affecting over 8 Mb, including insertions and deletions spanning thousands of base pairs. The potential of these variants to affect allele-specific expression was further explored. RNA-sequencing data from 11 different tissue types were mapped against the scaffolded haploid assembly and gene expression data are incorporated into our existing easy-to-use web-based interface to facilitate use by the broader plant science community. These two assemblies provide an excellent means to study the effects of heterozygosity, haplotype-specific structural variation, gene hemizygosity, and allele-specific gene expression contributing to important agricultural traits and further our understanding of the genetics and domestication of cassava.

Comprehensive analysis and expression profiles of cassava UDP-glycosyltransferases (UGT) family reveal their involvement in development and stress responses in cassava

UDP-glycosyltransferases (UGTs) are widely involved in plant growth and stress responses. However, UGT family are not well understood in cassava. Here, we identified 121 MeUGT genes and classified them into 14 subfamilies by phylogenetic analysis. All MeUGT proteins have typical feature of the UGTs family. Tandem duplications are the crucial driving force for the expansion of MeUGT family. Cis-Acting elements analysis uncovered those 14 kinds of cis-elements associated with biotic and abiotic stress responses. Transcriptomic and qRT-PCR analyses indicated that MeUGT genes participate in postharvest physiological deterioration of storage root and the responses of biotic and abiotic stresses. Of which, MeUGT-14/41 were significantly induced after Xam treatment. Silencing of MeUGT-14 or MeUGT-41 reduced cassava resistance to Xam, verifying the accuracy of transcriptomic data for function prediction. Together, this study characterized the MeUGTs family and revealed their potential functions, which build a solid foundation for MeUGTs associated genetic improvement of cassava.

Autophagy-related genes serve as heat shock protein 90 co-chaperones in disease resistance against cassava bacterial blight

Heat shock protein 90 (HSP90) is involved in plant growth and various stress responses via regulating protein homeostasis. Autophagy keeps cellular homeostasis by recycling the components of cellular cytoplasmic constituents. Although they have similar effects on cellular protein homeostasis, the direct association between HSP90 and autophagy signaling remains unclear in plants, especially in tropical crops. In this study, the correlation between HSP90 and autophagy signaling was systematically analyzed by protein-protein interaction in cassava, one of the most important economy fruit in tropic. In addition, their effects on plant disease response and underlying mechanisms in cassava were investigated by functional genomics and genetic phenotype assay. The potential MeHSP90.9-MeSGT1-MeRAR1 chaperone complex interacts with MeATGs and subsequently triggers autophagy signaling, conferring improved disease resistance to cassava bacterial blight (CBB). On the contrary, HSP90 inhibitor and autophagy inhibitor decreased disease resistance against CBB in cassava, and autophagy may be involved in the potential MeHSP90.9-MeSGT1-MeRAR1 chaperone complex-mediated multiple immune responses. This study highlights the precise modulation of autophagy signaling by potential MeHSP90.9-MeSGT1-MeRAR1 chaperone complex in autophagy-mediated disease resistance to CBB.

Metabolomic and transcriptomic profiling reveals distinct nutritional properties of cassavas with different flesh colors

•

A total of 508 metabolites were identified in three cassava cultivars.

•

White-fleshed cassava had the highest contents of amino acids and organic acids.

•

Yellow-fleshed cassava was enriched in metabolites related to specific pathways.

•

Several pathways were found to be regulated at the transcriptional level.

Cassava is a significant food security crop in several developing countries. Metabolites in cassava roots provide numerous nutrients essential for human health. Exploiting the diversity of nutritional ingredients present in cassavas is vital for improving its nutritional value. To address this problem, root metabolomes of three cassava cultivars with white-flesh, light-yellow-flesh and yellow-flesh were comprehensively measured, respectively. A total of 508 metabolites were detected in cassava roots, including 300 primary metabolites and 185 secondary metabolites. There were 22.6% to 34.1% metabolites exhibiting significant variations among the three cassava cultivars. The light-yellow-flesh cassava contained higher contents of secondary metabolites, especially flavone, phenylpropanoids and alkaloids, and lower contents of primary metabolites except lipids, alcohols, vitamins and derivatives. Compared with light-yellow-flesh cassava, the yellow-flesh cassava contained higher contents of amino acid and derivatives, but lower contents of phenylpropanoids, nucleotide and derivates. White-flesh cassava contained higher contents of primary metabolites, especially amino acid and derivatives, but lower contents of secondary metabolites except flavonoid and indole derivatives. Transcriptome analyses were parallelly performed to decipher the potential mechanisms regulating the accumulations of related metabolites. Several pathways were both enriched by differentially expressed genes and differentially accumulated metabolites, supporting that metabolisms of these metabolites were regulated at transcriptional level. These results expand the knowledge on metabolite compositions in cassava roots and provide substantial information for genetic improvement of cassavas with high nutritional values.

Overexpression of Cassava MeAnn2 Enhances the Salt and IAA Tolerance of Transgenic Arabidopsis

Annexins are a superfamily of soluble calcium-dependent phospholipid-binding proteins that have considerable regulatory effects in plants, especially in response to adversity and stress. The Arabidopsis thaliana AtAnn1 gene has been reported to play a significant role in various abiotic stress responses. In our study, the cDNA of an annexin gene highly similar to AtAnn1 was isolated from the cassava genome and named MeAnn2. It contains domains specific to annexins, including four annexin repeat sequences (I–IV), a Ca²⁺-binding sequence, Ca²⁺-independent membrane-binding-related tryptophan residues, and a salt bridge-related domain. MeAnn2 is localized in the cell membrane and cytoplasm, and it was found to be preferentially expressed in the storage roots of cassava. The overexpression of MeAnn2 reduced the sensitivity of transgenic Arabidopsis to various Ca²⁺, NaCl, and indole-3-acetic acid (IAA) concentrations. The expression of the stress resistance-related gene AtRD29B and auxin signaling pathway-related genes AtIAA4 and AtLBD18 in transgenic Arabidopsis was significantly increased under salt stress, while the Malondialdehyde (MDA) content was significantly lower than that of the control. These results indicate that the MeAnn2 gene may increase the salt tolerance of transgenic Arabidopsis via the IAA signaling pathway.

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.

Gene tagging via CRISPR-mediated homology-directed repair in cassava

Research on a few model plant-pathogen systems has benefitted from years of tool and resource development. This is not the case for the vast majority of economically and nutritionally important plants, creating a crop improvement bottleneck. Cassava bacterial blight (CBB), caused by Xanthomonas axonopodis pv. manihotis (Xam), is an important disease in all regions where cassava (Manihot esculenta Crantz) is grown. Here, we describe the development of cassava that can be used to visualize one of the initial steps of CBB infection in vivo. Using CRISPR-mediated homology-directed repair (HDR), we generated plants containing scarless insertion of GFP at the 3' end of CBB susceptibility (S) gene MeSWEET10a. Activation of MeSWEET10a-GFP by the transcription activator-like (TAL) effector TAL20 was subsequently visualized at transcriptional and translational levels. To our knowledge, this is the first such demonstration of HDR via gene editing in cassava.

Genome-Wide Analysis of Tubulin Gene Family in Cassava and Expression of Family Member FtsZ2-1 during Various Stress

Filamentous temperature-sensitive protein Z (Tubulin/FtsZ) family is a group of conserved GTP-binding (guanine nucleotide-binding) proteins, which are closely related to plant tissue development and organ formation as the major component of the cytoskeleton. According to the published genome sequence information of cassava (Manihot esculenta Crantz), 23 tubulin genes (MeTubulins) were identified, which were divided into four main groups based on their type and phylogenetic characteristics. The same grouping generally has the same or similar motif composition and exon–intron structure. Collinear analysis showed that fragment repetition event is the main factor in amplification of cassava tubulin superfamily gene. The expression profiles of MeTubulin genes in various tissue were analyzed, and it was found that MeTubulins were mainly expressed in leaf, petiole, and stem, while FtsZ2-1 was highly expressed in storage root. The qRT-PCR results of the FtsZ2-1 gene under hormone and abiotic stresses showed that indole-3-acetic acid (IAA) and gibberellin A3 (GA3) stresses could significantly increase the expression of the FtsZ2-1 gene, thereby revealing the potential role of FtsZ2-1 in IAA and GA3 stress-induced responses.

Viral inosine triphosphatase: A mysterious enzyme with typical activity, but an atypical function

Plant viruses typically have highly condensed genomes, yet the plant-pathogenic viruses Cassava brown streak virus, Ugandan cassava brown streak virus, and Euphorbia ringspot virus are unusual in encoding an enzyme not yet found in any other virus, the "house-cleaning" enzyme inosine triphosphatase. Inosine triphosphatases (ITPases) are highly conserved enzymes that occur in all kingdoms of life and perform a house-cleaning function by hydrolysing the noncanonical nucleotide inosine triphosphate to inosine monophosphate. The ITPases encoded by cassava brown streak virus and Ugandan cassava brown streak virus have been characterized biochemically and are shown to have typical ITPase activity. However, their biological role in virus infection has yet to be elucidated. Here we review what is known of viral-encoded ITPases and speculate on potential roles in infection with the aim of generating a greater understanding of cassava brown streak viruses, a group of the world's most devastating viruses.

Heat shock protein 90 co-chaperone modules fine-tune the antagonistic interaction between salicylic acid and auxin biosynthesis in cassava

Heat shock protein 90 (HSP90) is an important molecular chaperone in plants. However, HSP90-mediated plant immune response remains elusive in cassava. In this study, cassava bacterial blight (CBB) induces the expression of MeHsf8, which directly targets MeHSP90.9 to activate its expression and immune response. Further identification of SHI-related sequence 1 (MeSRS1) and MeWRKY20 as MeHSP90.9 co-chaperones revealed the underlying mechanism of MeHSP90.9-mediated immune response. MeHSP90.9 interacts with MeSRS1 and MeWRKY20 to promote their transcriptional activation of salicylic acid (SA) biosynthetic gene avrPphB Susceptible 3 (MePBS3) and tryptophan metabolic gene N-acetylserotonin O-methyltransferase 2 (MeASMT2), respectively, so as to activate SA biosynthesis but inhibit tryptophan-derived auxin biosynthesis. Notably, genetic experiments confirmed that overexpressing MePBS3 and MeASMT2 could rescue the effects of silencing MeHsf8-MeHSP90.9 on disease resistance. This study highlights the dual regulation of SA and auxin biosynthesis by MeHSP90.9, providing the mechanistic understanding of MeHSP90.9 client partners in plant immunity.

Large-scale genome-wide association study, using historical data, identifies conserved genetic architecture of cyanogenic glucoside content in cassava (Manihot esculenta Crantz) root

Manihot esculenta (cassava) is a root crop originating from South America that is a major staple in the tropics, including in marginal environments. This study focused on South American and African germplasm and investigated the genetic architecture of hydrogen cyanide (HCN), a major component of root quality. HCN, representing total cyanogenic glucosides, is a plant defense component against herbivory but is also toxic for human consumption. We genotyped 3354 landraces and modern breeding lines originating from 26 Brazilian states and 1389 individuals were phenotypically characterized across multi-year trials for HCN. All plant material was subjected to high-density genotyping using genotyping by sequencing. We performed genome-wide association mapping to characterize the genetic architecture and gene mapping of HCN. Field experiments revealed strong broad- and narrow-sense trait heritability (0.82 and 0.41, respectively). Two major loci were identified, encoding for an ATPase and a MATE protein, and contributing up to 7 and 30% of the HCN concentration in roots, respectively. We developed diagnostic markers for breeding applications, validated trait architecture consistency in African germplasm and investigated further evidence for the domestication of sweet and bitter cassava. Fine genomic characterization revealed: (i) the major role played by vacuolar transporters in regulating HCN content; (ii) the co-domestication of sweet and bitter cassava major alleles are dependent upon geographical zone; and (iii) the major loci allele for high HCN in M. esculenta Crantz seems to originate from its ancestor, M. esculenta subsp. flabellifolia. Taken together, these findings expand our insights into cyanogenic glucosides in cassava roots and its glycosylated derivatives in plants.

CAMTA3 negatively regulates disease resistance through modulating immune response and extensive transcriptional reprogramming in cassava

As one of the important crops in the world, cassava production is seriously threatened by Xanthomonas axonopodis pv. manihotis (Xam) all year round. Calmodulin-binding transcription activators (CAMTAs) play key roles in biotic stress and abiotic stress in plants, however, their roles in cassava remain elusive. In this study, six MeCAMTAs were identified, and MeCAMTA3 with the highest induction upon Xam infection was confirmed as a transcription factor that binds to the vCGCGb motif. MeCAMTA3 negatively regulates plant disease resistance against Xam. On the one hand, MeCAMTA3 negatively regulated endogenous salicylic acid and reactive oxygen species accumulation, pathogenesis-related genes MePRs' transcripts and callose deposition during cassava-Xam interaction but not under control conditions. On the other hand, RNA sequencing showed extensive transcriptional reprogramming by MeCAMTA3, especially 18 genes with a vCGCGb motif in the promoter region in hormone signaling, antioxidant signaling and other disease resistance signaling. Notably, chromatin immunoprecipitation-polymerase chain reaction showed that eight of these genes might be directly regulated by MeCAMTA3 through transcriptional repression. In summary, MeCAMTA3 negatively regulates plant disease resistance against cassava bacterial blight through modulation of multiple immune responses during cassava-Xam interaction and extensive transcriptional reprogramming.

Melatonin synthesis enzymes interact with ascorbate peroxidase to protect against oxidative stress in cassava

Melatonin is an important indole amine hormone in animals and plants. The enzymes that catalyse melatonin synthesis positively regulate plant stress responses through modulation of the accumulation of reactive oxygen species (ROS). However, the relationship between melatonin biosynthetic enzymes and ROS-scavenging enzymes has not been characterized. In this study, we demonstrate that two enzymes of the melatonin synthesis pathway in Manihot esculenta (MeTDC2 and MeASMT2) directly interact with ascorbate peroxidase (MeAPX2) in both in vitro and in vivo experiments. Notably, in the presence of MeTDC2 and MeASMT2, MeAPX2 showed significantly higher activity and antioxidant capacity than the purified MeAPX2 protein alone. These findings indicate that MeTDC2-MeAPX2 and MeASMT2-MeAPX2 interactions both activate APX activity and increase antioxidant capacity. In addition, the combination of MeTDC2, MeASMT2, and MeAPX2 conferred improved resistance to hydrogen peroxide in Escherichia coli. Moreover, this combination also positively regulates oxidative stress tolerance in cassava. Taken together, these findings not only reveal a direct interaction between MeTDC2, MeASMT2, and MeAPX2, but also highlight the importance of this interaction in regulating redox homoeostasis and stress tolerance in cassava.

Analysis of leaf morphology, secondary metabolites and proteins related to the resistance to Tetranychus cinnabarinus in cassava (Manihot esculenta Crantz)

Constitutive resistance of plant can be divided into physical and chemical barriers. Cassava (Manihot esculenta Crantz) is susceptible to mites, especially Tetranychus cinnabarinus. Although significant differences in the resistance to T. cinnabarinus are observed in different cassava cultivars, limited research has been done on the mechanism accounting for the resistance. The aim of this study was to explore the mechanism of resistance to T. cinnabarinus by comparing morphology, secondary metabolites and proteins in different cassava cultivars. The anatomical structure of leaves showed that the cassava cultivar Xinxuan 048 (XX048), which showed a stronger resistance to T. cinnabarinus in both greenhouse testing and three years field evaluation tests (2016–2018), had thicker palisade tissue, spongy tissue, lower epidermis and leaf midrib tissue compared to cultivar Guire 4 (GR4). Greenhouse evaluation demonstrated that originally these cultivars were different, leading to differences in constitutive levels of metabolites. The proteomic analysis of protected leaves in XX048 and GR4 revealed that up-regulated differentially expressed proteins (DEPs) were highly enriched in secondary metabolic pathways, especially in the biosynthesis of flavonoids. This study not only provides a comprehensive data set for overall proteomic changes of leaves in resistant and susceptible cassava, but also sheds light on the morphological characteristics of cassava-mite interaction, secondary metabolite defense responses, and molecular breeding of mite-resistant cassava for effective pest control.

Genomic analyses of heat stress transcription factors (HSFs) in simulated drought stress response and storage root deterioration after harvest in cassava

Heat shock factors (HSFs) play crucial roles in various plant stress responses. However, the current knowledge about HSFs in cassava, an important crop, is still insufficient. In this research, we identified 32 cassava HSF genes (MeHSFs) and clustered them into three groups (A, B, C) based on phylogenetic analysis and structural characteristics. Conserved motif analyses showed that MeHSFs display domains characteristic to HSF transcription factors. Gene structure analyses suggested that 29 MeHSFs contained only two exons. All identified 32 cassava MeHSFs were distributed on 13 chromosomes. Their expression profiles revealed that the different MeHSFs were expressed differentially in different tissues, most high expression genes belonged to group A. The similar MeHSFs were up-regulated after treatment with both PEG and abscisic acid (ABA), which implied that these MeHSFs may participate in resistance to simulated drought stress associated with the ABA signaling pathway. In addition, several MeHSFs were induced during postharvest physiological deterioration (PPD) in cassava. Our results provided basic but important knowledge for future gene function analysis of MeHSFs toward efforts in improving tolerance to abiotic stress and PPD in cassava.

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.

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.

Cassava breeding and agronomy in Asia: 50 years of history and future directions

In Asia, cassava (Manihot esculenta) is cultivated by more than 8 million farmers, driving the rural economy of many countries. The International Center for Tropical Agriculture (CIAT), in partnership with national agricultural research institutes (NARIs), instigated breeding and agronomic research in Asia, 1983. The breeding program has successfully released high-yielding cultivars resulting in an average yield increase from 13.0 t ha-1 in 1996 to 21.3 t ha-1 in 2016, with significant economic benefits. Following the success in increasing yields, cassava breeding has turned its focus to higher-value traits, such as waxy cassava, to reach new market niches. More recently, building resistance to invasive pests and diseases has become a top priority due to the emergent threat of cassava mosaic disease (CMD). The agronomic research involves driving profitability with advanced technologies focusing on better agronomic management practices thereby maintaining sustainable production systems. Remote sensing technologies are being tested for trait discovery and large-scale field evaluation of cassava. In summary, cassava breeding in Asia is driven by a combination of food and market demand with technological innovations to increase the productivity. Further, exploration in the potential of data-driven agriculture is needed to empower researchers and producers for sustainable advancement.

Genomic and Transcriptomic Analysis Identified Novel Putative Cassava lncRNAs Involved in Cold and Drought Stress

Long non-coding RNAs (lncRNAs) play important roles in the regulation of complex cellular processes, including transcriptional and post-transcriptional regulation of gene expression relevant for development and stress response, among others. Compared to other important crops, there is limited knowledge of cassava lncRNAs and their roles in abiotic stress adaptation. In this study, we performed a genome-wide study of ncRNAs in cassava, integrating genomics- and transcriptomics-based approaches. In total, 56,840 putative ncRNAs were identified, and approximately half the number were verified using expression data or previously known ncRNAs. Among these were 2229 potential novel lncRNA transcripts with unmatched sequences, 250 of which were differentially expressed in cold or drought conditions, relative to controls. We showed that lncRNAs might be involved in post-transcriptional regulation of stress-induced transcription factors (TFs) such as zinc-finger, WRKY, and nuclear factor Y gene families. These findings deepened our knowledge of cassava lncRNAs and shed light on their stress-responsive roles.