Tapioca: a platform for predicting de novo protein-protein interactions in dynamic contexts.

Protein-protein interactions (PPIs) drive cellular processes and responses to environmental cues, reflecting the cellular state. Here we develop Tapioca, an ensemble machine learning framework for studying global PPIs in dynamic contexts. Tapioca predicts de novo interactions by integrating mass spectrometry interactome data from thermal/ion denaturation or cofractionation workflows with protein properties and tissue-specific functional networks. Focusing on the thermal proximity coaggregation method, we improved the experimental workflow. Finely tuned thermal denaturation afforded increased throughput, while cell lysis optimization enhanced protein detection from different subcellular compartments. The Tapioca workflow was next leveraged to investigate viral infection dynamics. Temporal PPIs were characterized during the reactivation from latency of the oncogenic Kaposi's sarcoma-associated herpesvirus. Together with functional assays, NUCKS was identified as a proviral hub protein, and a broader role was uncovered by integrating PPI networks from alpha- and betaherpesvirus infections. Altogether, Tapioca provides a web-accessible platform for predicting PPIs in dynamic contexts.

Early Holocene crop cultivation and landscape modification in Amazonia.

The onset of plant cultivation is one of the most important cultural transitions in human history1-4. Southwestern Amazonia has previously been proposed as an early centre of plant domestication, on the basis of molecular markers that show genetic similarities between domesticated plants and wild relatives4-6. However, the nature of the early human occupation of southwestern Amazonia, and the history of plant cultivation in this region, are poorly understood. Here we document the cultivation of squash (Cucurbita sp.) at about 10,250 calibrated years before present (cal. yr BP), manioc (Manihot sp.) at about 10,350 cal. yr BP and maize (Zea mays) at about 6,850 cal. yr BP, in the Llanos de Moxos (Bolivia). We show that, starting at around 10,850 cal. yr BP, inhabitants of this region began to create a landscape that ultimately comprised approximately 4,700 artificial forest islands within a treeless, seasonally flooded savannah. Our results confirm that the Llanos de Moxos is a hotspot for early plant cultivation and demonstrate that-ever since their arrival in Amazonia-humans have markedly altered the landscape, with lasting repercussions for habitat heterogeneity and species conservation.

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.

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.

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.

Cassava phosphatase PP2C1 modulates thermotolerance via fine-tuning dephosphorylation of antioxidant enzymes.

Global warming is an adverse environmental factor that threatens crop yields and food security. 2C-type protein phosphatases (PP2Cs), as core protein phosphatase components, play important roles in plant hormone signaling to cope with various environmental stresses. However, the function and underlying mechanism of PP2Cs in the heat stress response remain elusive in tropical crops. Here, we report that MePP2C1 negatively regulated thermotolerance in cassava (Manihot esculenta Crantz), accompanied by the modulation of reactive oxygen species (ROS) accumulation and the underlying antioxidant enzyme activities of catalase (CAT) and ascorbate peroxidase (APX). Further investigation found that MePP2C1 directly interacted with and dephosphorylated MeCAT1 and MeAPX2 at serine (S) 112 and S160 residues, respectively. Moreover, in vitro and in vivo assays showed that protein phosphorylation of MeCAT1S112 and MeAPX2S160 was essential for their enzyme activities, and MePP2C1 negatively regulated thermotolerance and redox homeostasis by dephosphorylating MeCAT1S112 and MeAPX2S160. Taken together, this study illustrates the direct relationship between MePP2C1-mediated protein dephosphorylation of MeCAT1 and MeAPX2 and ROS accumulation in thermotolerance to provide insights for adapting to global warming via fine-tuning thermotolerance of the tropical crop cassava.

RNA splicing modulates the postharvest physiological deterioration of cassava storage root.

Rapid postharvest physiological deterioration (PPD) of cassava (Manihot esculenta Crantz) storage roots is a major constraint that limits the potential of this plant as a food and industrial crop. Extensive studies have been performed to explore the regulatory mechanisms underlying the PPD processes in cassava to understand their molecular and physiological responses. However, the exceptional functional versatility of alternative splicing (AS) remains to be explored during the PPD process in cassava. Here, we identified several aberrantly spliced genes during the early PPD stage. An in-depth analysis of AS revealed that the abscisic acid (ABA) biosynthesis pathway might serve as an additional molecular layer in attenuating the onset of PPD. Exogenous ABA application alleviated PPD symptoms through maintaining ROS generation and scavenging. Interestingly, the intron retention transcript of MeABA1 (ABA DEFICIENT 1) was highly correlated with PPD symptoms in cassava storage roots. RNA yeast 3-hybrid and RNA immunoprecipitation (RIP) assays showed that the serine/arginine-rich protein MeSCL33 (SC35-like splicing factor 33) binds to the precursor mRNA of MeABA1. Importantly, overexpressing MeSCL33 in cassava conferred improved PPD resistance by manipulating the AS and expression levels of MeABA1 and then modulating the endogenous ABA levels in cassava storage roots. Our results uncovered the pivotal role of the ABA biosynthesis pathway and RNA splicing in regulating cassava PPD resistance and proposed the essential roles of MeSCL33 for conferring PPD resistance, broadening our understanding of SR proteins in cassava development and stress responses.

Coat protein of cassava common mosaic virus targets RAV1 and RAV2 transcription factors to subvert immunity in cassava.

Cassava common mosaic virus (CsCMV, genus Potexvirus) is a prevalent virus associated with cassava mosaic disease, so it is essential to elucidate the underlying molecular mechanisms of the coevolutionary arms race between viral pathogenesis and the cassava (Manihot esculenta Crantz) defense response. However, the molecular mechanism underlying CsCMV infection is largely unclear. Here, we revealed that coat protein (CP) acts as a major pathogenicity determinant of CsCMV via a mutant infectious clone. Moreover, we identified the target proteins of CP-related to abscisic acid insensitive3 (ABI3)/viviparous1 (VP1) (MeRAV1) and MeRAV2 transcription factors, which positively regulated disease resistance against CsCMV via transcriptional activation of melatonin biosynthetic genes (tryptophan decarboxylase 2 (MeTDC2), tryptamine 5-hydroxylase (MeT5H), N-aceylserotonin O-methyltransferase 1 (MeASMT1)) and MeCatalase6 (MeCAT6) and MeCAT7. Notably, the interaction between CP, MeRAV1, and MeRAV2 interfered with the protein phosphorylation of MeRAV1 and MeRAV2 individually at Ser45 and Ser44 by the protein kinase, thereby weakening the transcriptional activation activity of MeRAV1 and MeRAV2 on melatonin biosynthetic genes, MeCAT6 and MeCAT7 dependent on the protein phosphorylation of MeRAV1 and MeRAV2. Taken together, the identification of the CP-MeRAV1 and CP-MeRAV2 interaction module not only illustrates a molecular mechanism by which CsCMV orchestrates the host defense system to benefit its infection and development but also provides a gene network with potential value for the genetic improvement of cassava disease resistance.

Symplasmic phloem loading and subcellular transport in storage roots are key factors for carbon allocation in cassava.

Cassava (Manihot esculenta) is a deciduous woody perennial shrub that stores large amounts of carbon and water in its storage roots. Previous studies have shown that assimilating unloading into storage roots happens symplasmically once secondary anatomy is established. However, mechanisms controlling phloem loading and overall carbon partitioning to different cassava tissues remain unclear. Here, we used a combination of histological, transcriptional, and biochemical analyses on different cassava tissues and at different timepoints to better understand source-sink carbon allocation. We found that cassava likely utilizes a predominantly passive symplasmic phloem loading strategy, indicated by the lack of expression of genes coding for key players of sucrose transport, the existence of branched plasmodesmata in the companion cell/bundle sheath interface of minor leaf veins, and very high leaf sucrose concentrations. Furthermore, we showed that tissue-specific changes in anatomy and non-structural carbohydrate contents are associated with tissue-specific modification in gene expression for sucrose cleavage/synthesis, as well as subcellular compartmentalization of sugars. Overall, our data suggest that carbon allocation during storage root filling is mostly facilitated symplasmically and is likely mostly regulated by local tissue demand and subcellular compartmentalization.

CRISPR/Cas9-generated mutations in a sugar transporter gene reduce cassava susceptibility to bacterial blight.

Bacteria from the genus Xanthomonas are prolific phytopathogens that elicit disease in over 400 plant species. Xanthomonads carry a repertoire of specialized proteins called transcription activator-like (TAL) effectors that promote disease and pathogen virulence by inducing the expression of host susceptibility (S) genes. Xanthomonas phaseoli pv. manihotis (Xpm) causes bacterial blight on the staple food crop cassava (Manihot esculenta Crantz). The Xpm effector TAL20 induces ectopic expression of the S gene Manihot esculenta Sugars Will Eventually be Exported Transporter 10a (MeSWEET10a), which encodes a sugar transporter that contributes to cassava bacterial blight (CBB) susceptibility. We used CRISPR/Cas9 to generate multiple cassava lines with edits to the MeSWEET10a TAL20 effector binding site and/or coding sequence. In several of the regenerated lines, MeSWEET10a expression was no longer induced by Xpm, and in these cases, we observed reduced CBB disease symptoms post Xpm infection. Because MeSWEET10a is expressed in cassava flowers, we further characterized the reproductive capability of the MeSWEET10a promoter and coding sequence mutants. Lines were crossed to themselves and to wild-type plants. The results indicated that expression of MeSWEET10a in female, but not male, flowers is critical to produce viable F1 seed. In the case of promoter mutations that left the coding sequence intact, viable F1 progeny were recovered. Taken together, these results demonstrate that blocking MeSWEET10a induction is a viable strategy for decreasing cassava susceptibility to CBB and that ideal lines will contain promoter mutations that block TAL effector binding while leaving endogenous expression of MeSWEET10a unaltered.

Dynamic DNA methylation modifications in the cold stress response of cassava.

Cassava, a crucial tropical crop, faces challenges from cold stress, necessitating an exploration of its molecular response. Here, we investigated the role of DNA methylation in moderating the response to moderate cold stress (10 °C) in cassava. Using whole-genome bisulfite sequencing, we examined DNA methylation patterns in leaf blades and petioles under control conditions, 5 h, and 48 h of cold stress. Tissue-specific responses were observed, with leaf blades exhibiting subtle changes, while petioles displayed a pronounced decrease in methylation levels under cold stress. We identified cold stress-induced differentially methylated regions (DMRs) that demonstrated both tissue and treatment specificity. Importantly, these DMRs were enriched in genes with altered expression, implying functional relevance. The cold-response transcription factor ERF105 associated with DMRs emerged as a significant and conserved regulator across tissues and treatments. Furthermore, we investigated DNA methylation dynamics in transposable elements, emphasizing the sensitivity of MITEs with bHLH binding motifs to cold stress. These findings provide insights into the epigenetic regulation of response to cold stress in cassava, contributing to an understanding of the molecular mechanisms underlying stress adaptation in this tropical plant.

Efficient sugar utilization and transition from oxidative to substrate-level phosphorylation in high starch storage roots of African cassava genotypes.

Cassava's storage roots represent one of the most important sources of nutritional carbohydrates worldwide. Particularly, smallholder farmers in sub-Saharan Africa depend on this crop plant, where resilient and yield-improved varieties are of vital importance to support steadily increasing populations. Aided by a growing understanding of the plant's metabolism and physiology, targeted improvement concepts already led to visible gains in recent years. To expand our knowledge and to contribute to these successes, we investigated storage roots of eight cassava genotypes with differential dry matter content from three successive field trials for their proteomic and metabolic profiles. At large, the metabolic focus in storage roots transitioned from cellular growth processes toward carbohydrate and nitrogen storage with increasing dry matter content. This is reflected in higher abundance of proteins related to nucleotide synthesis, protein turnover, and vacuolar energization in low starch genotypes, while proteins involved in sugar conversion and glycolysis were more prevalent in high dry matter genotypes. This shift in metabolic orientation was underlined by a clear transition from oxidative- to substrate-level phosphorylation in high dry matter genotypes. Our analyses highlight metabolic patterns that are consistently and quantitatively associated with high dry matter accumulation in cassava storage roots, providing fundamental understanding of cassava's metabolism as well as a data resource for targeted genetic improvement.

Comparative transcriptomics analysis reveals defense mechanisms of Manihot esculenta Crantz against Sri Lanka Cassava MosaicVirus.

BACKGROUND: Cassava mosaic disease (CMD), caused by Sri Lankan cassava mosaic virus (SLCMV) infection, has been identified as a major pernicious disease in Manihot esculenta Crantz (cassava) plantations. It is widespread in Southeast Asia, especially in Thailand, which is one of the main cassava supplier countries. With the aim of restricting the spread of SLCMV, we explored the gene expression of a tolerant cassava cultivar vs. a susceptible cassava cultivar from the perspective of transcriptional regulation and the mechanisms underlying plant immunity and adaptation. RESULTS: Transcriptomic analysis of SLCMV-infected tolerant (Kasetsart 50 [KU 50]) and susceptible (Rayong 11 [R 11]) cultivars at three infection stages-that is, at 21 days post-inoculation (dpi) (early/asymptomatic), 32 dpi (middle/recovery), and 67 dpi (late infection/late recovery)-identified 55,699 expressed genes. Differentially expressed genes (DEGs) between SLCMV-infected KU 50 and R 11 cultivars at (i) 21 dpi to 32 dpi (the early to middle stage), and (ii) 32 dpi to 67 dpi (the middle stage to late stage) were then identified and validated by real-time quantitative PCR (RT-qPCR). DEGs among different infection stages represent genes that respond to and regulate the viral infection during specific stages. The transcriptomic comparison between the tolerant and susceptible cultivars highlighted the role of gene expression regulation in tolerant and susceptible phenotypes. CONCLUSIONS: This study identified genes involved in epigenetic modification, transcription and transcription factor activities, plant defense and oxidative stress response, gene expression, hormone- and metabolite-related pathways, and translation and translational initiation activities, particularly in KU 50 which represented the tolerant cultivar in this study.

Morphological, anatomical, and transcriptomics analysis reveals the regulatory mechanisms of cassava plant height development.

BACKGROUND: Cassava is one of three major potato crops and the sixth most important food crop globally. Improving yield remains a primary aim in cassava breeding. Notably, plant height significantly impacts the yield and quality of crops; however, the mechanisms underlying cassava plant height development are yet to be elucidated. RESULTS: In this study, we investigated the mechanisms responsible for cassava plant height development using phenotypic, anatomical, and transcriptomic analyses. Phenotypic and anatomical analysis revealed that compared to the high-stem cassava cultivar, the dwarf-stem cassava cultivar exhibited a significant reduction in plant height and a notable increase in internode tissue xylem area. Meanwhile, physiological analysis demonstrated that the lignin content of dwarf cassava was significantly higher than that of high cassava. Notably, transcriptome analysis of internode tissues identified several differentially expressed genes involved in cell wall synthesis and expansion, plant hormone signal transduction, phenylpropanoid biosynthesis, and flavonoid biosynthesis between the two cassava cultivars. CONCLUSIONS: Our findings suggest that internode tissue cell division, secondary wall lignification, and hormone-related gene expression play important roles in cassava plant height development. Ultimately, this study provides new insights into the mechanisms of plant height morphogenesis in cassava and identifies candidate regulatory genes associated with plant height that can serve as valuable genetic resources for future crop dwarfing breeding.

Transcriptomic study of the role of MeFtsZ2-1 in pigment accumulation in cassava leaves.

MeFtsZ2-1 is a key gene for plant plastid division, but the mechanism by which MeFtsZ2-1 affects pigment accumulation in cassava (Manihot esculenta Crantz) through plastids remains unclear. We found that MeFtsZ2-1 overexpression in cassava (OE) exhibited darker colors of leaves, with increased levels of anthocyanins and carotenoids. Further observation via Transmission Electron Microscopy (TEM) revealed no apparent defects in chloroplast structure but an increase in the number of plastoglobule in OE leaves. RNA-seq results showed 1582 differentially expressed genes (DEGs) in leaves of OE. KEGG pathway analysis indicated that these DEGs were enriched in pathways related to flavonoid, anthocyanin, and carotenoid biosynthesis. This study reveals the role of MeFtsZ2-1 in cassava pigment accumulation from a physiological and transcriptomic perspective, providing a theoretical basis for improving cassava quality.

Identification of the 4CL family in cassava (Manihot esculenta Crantz) and expression pattern analysis of the Me4CL32 gene.

The 4-coumarate coenzyme A ligase (4CL) plays a critical role in the phenylpropane metabolic pathway and is a key enzyme in plant growth metabolism and stress responses. Using bioinformatics methods, 50 Me4CL gene were identified within the cassava genome u, and a comprehensive analysis of the cassava 4CL gene family was conducted. The results showed that these 50 4CL proteins are divided into four subfamilies, with members within the same subfamily sharing similar or identical gene structures. Co-linearity analysis revealed that cassava and rubber trees have the highest number of homologous genes, indicating a close homologous relationship between them. Analysis of 20 cis-acting elements in the promoter region of Me4CL32 revealed the presence of hormone-responsive elements such as gibberellin, auxin, abscisic acid, and as well as elements related to meristematic tissue regulation. results Quantitative real-time PCR (qRT-PCR) results showed that Me4CL32 gene expression changes in response to abiotic stressors (drought, salt, cold, heat) and hormonal stimuli(GA3 and ABA), indicating that Me4CL32 can respond to both environmental pressures and hormone signals. RNA-seq transcriptome and single-cell transcriptome analyses were used to examine the expression patterns of Me4CLs. Additionally, subcellular localization studies demonstrated that the Me4CL32 protein is confined to the chloroplasts of cassava leaves.Investigating the functionality of this gene family aids in comprehending plant growth, development, and stress resistance mechanisms. Furthermore, it furnishes a theoretical basis for future research on developing resilient cassava germplasm and the enhancing cassava's environmental tolerance.

A high-affinity potassium transporter (MeHKT1) from cassava (Manihot esculenta) negatively regulates the response of transgenic Arabidopsis to salt stress.

BACKGROUND: High-affinity potassium transporters (HKTs) are crucial in facilitating potassium uptake by plants. Many types of HKTs confer salt tolerance to plants through regulating K+ and Na+ homeostasis under salinity stress. However, their specific functions in cassava (Manihot esculenta) remain unclear. RESULTS: Herein, an HKT gene (MeHKT1) was cloned from cassava, and its expression is triggered by exposure to salt stress. The expression of a plasma membrane-bound protein functions as transporter to rescue a low potassium (K+) sensitivity of yeast mutant strain, but the complementation of MeHKT1 is inhibited by NaCl treatment. Under low K+ stress, transgenic Arabidopsis with MeHKT1 exhibits improved growth due to increasing shoot K+ content. In contrast, transgenic Arabidopsis accumulates more Na+ under salt stress than wild-type (WT) plants. Nevertheless, the differences in K+ content between transgenic and WT plants are not significant. Additionally, Arabidopsis expressing MeHKT1 displayed a stronger salt-sensitive phenotype. CONCLUSION: These results suggest that under low K+ condition, MeHKT1 functions as a potassium transporter. In contrast, MeHKT1 mainly transports Na+ into cells under salt stress condition and negatively regulates the response of transgenic Arabidopsis to salt stress. Our results provide a reference for further research on the function of MeHKT1, and provide a basis for further application of MeHKT1 in cassava by molecular biological means.

NF-YC15 transcription factor activates ethylene biosynthesis and improves cassava disease resistance.

The nuclear factor Y (NF-Y) transcription factors play important roles in plant development and physiological responses. However, the relationship between NF-Y, plant hormone and plant stress resistance in tropical crops remains unclear. In this study, we identified MeNF-YC15 gene in the NF-Y family that significantly responded to Xanthomonas axonopodis pv. manihotis (Xam) treatment. Using MeNF-YC15-silenced and -overexpressed cassava plants, we elucidated that MeNF-YC15 positively regulated disease resistance to cassava bacterial blight (CBB). Notably, we illustrated MeNF-YC15 downstream genes and revealed the direct genetic relationship between MeNF-YC15 and 1-aminocyclopropane-1-carboxylic acid (ACC) oxidase (MeACO1)-ethylene module in disease resistance, as evidenced by the rescued disease susceptibility of MeNF-YC15 silenced cassava plants with ethylene treatment or overexpressing MeACO1. In addition, the physical interaction between 2C-type protein phosphatase 1 (MePP2C1) and MeNF-YC15 inhibited the transcriptional activation of MeACO1 by MeNF-YC15. In summary, MePP2C1-MeNF-YC15 interaction modulates ethylene biosynthesis and cassava disease resistance, providing gene network for cassava genetic improvement.

KIN10-mediated HB16 protein phosphorylation and self-association improve cassava disease resistance by transcriptional activation of lignin biosynthesis genes.

Cassava bacterial blight significantly affects cassava yield worldwide, while major cassava cultivars are susceptible to this disease. Therefore, it is crucial to identify cassava disease resistance gene networks and defence molecules for the genetic improvement of cassava cultivars. In this study, we found that MeHB16 transcription factor as a differentially expressed gene in cassava cultivars with contrasting disease resistance, positively modulated disease resistance by modulating defence molecule lignin accumulation. Further investigation showed that MeHB16 physically interacted with itself via the leucine-Zippe domain (L-Zip), which was necessary for the transcriptional activation of downstream lignin biosynthesis genes. In addition, protein kinase MeKIN10 directly interacted with MeHB16 to promote its phosphorylation at Ser6, which in turn enhanced MeHB16 self-association and downstream lignin biosynthesis. In summary, this study revealed the molecular network of MeKIN10-mediated MeHB16 protein phosphorylation improved cassava bacterial blight resistance by fine-tuning lignin biosynthesis and provides candidate genes and the defence molecule for improving cassava disease resistance.

Maf/ham1-like pyrophosphatases of non-canonical nucleotides are host-specific partners of viral RNA-dependent RNA polymerases.

Cassava brown streak disease (CBSD), dubbed the "Ebola of plants", is a serious threat to food security in Africa caused by two viruses of the family Potyviridae: cassava brown streak virus (CBSV) and Ugandan (U)CBSV. Intriguingly, U/CBSV, along with another member of this family and one secoviridae, are the only known RNA viruses encoding a protein of the Maf/ham1-like family, a group of widespread pyrophosphatase of non-canonical nucleotides (ITPase) expressed by all living organisms. Despite the socio-economic impact of CDSD, the relevance and role of this atypical viral factor has not been yet established. Here, using an infectious cDNA clone and reverse genetics, we demonstrate that UCBSV requires the ITPase activity for infectivity in cassava, but not in the model plant Nicotiana benthamiana. HPLC-MS/MS experiments showed that, quite likely, this host-specific constraint is due to an unexpected high concentration of non-canonical nucleotides in cassava. Finally, protein analyses and experimental evolution of mutant viruses indicated that keeping a fraction of the yielded UCBSV ITPase covalently bound to the viral RNA-dependent RNA polymerase (RdRP) optimizes viral fitness, and this seems to be a feature shared by the other members of the Potyviridae family expressing Maf/ham1-like proteins. All in all, our work (i) reveals that the over-accumulation of non-canonical nucleotides in the host might have a key role in antiviral defense, and (ii) provides the first example of an RdRP-ITPase partnership, reinforcing the idea that RNA viruses are incredibly versatile at adaptation to different host setups.