Genome-wide association study of cassava brown streak disease resistance in cassava germplasm conserved in South America.

Cassava (Manihot esculenta Crantz) is a vital carbohydrate source for over 800 million people globally, yet its production in East Africa is severely affected by cassava brown streak disease (CBSD). Genebanks, through ex-situ conservation, play a pivotal role in preserving crop diversity, providing crucial resources for breeding resilient and disease-resistant crops. This study genotyped 234 South American cassava accessions conserved at the CIAT genebank, previously phenotyped for CBSD resistance by an independent group, to perform a genome-wide association analysis (GWAS) to identify genetic variants associated with CBSD resistance. Our GWAS identified 35 single nucleotide polymorphism (SNP) markers distributed across various chromosomes, associated with disease severity or the presence/absence of viral infection. Markers were annotated within or near genes previously identified with functions related to pathogen recognition and immune response activation. Using the SNP candidates, we screened the world's largest cassava collection for accessions with a higher frequency of favorable genotypes, proposing 35 accessions with potential resistance to CBSD. Our results provide insights into the genetics of CBSD resistance and highlight the importance of genetic resources to equip breeders with the raw materials needed to develop new crop varieties resistant to pests and diseases.

Genotype by environment interaction effect and fresh root yield stability of cassava genotypes under contrasting nitrogen regimes.

Nitrogen (N) is an important nutrient element needed by cassava for optimum yield and it is a vital component of nucleotides (nucleic acids), enzymes, amino acids (proteins), chlorophyll molecules and hormones, among other essential compounds required for growth and development of cassava. Nitrogen stress is a major cassava production constraint, the study aimed to examine genotype by environment interaction (GEI) effects and fresh root yield stability of 203 diverse cassava clones to identify genotypes with stable performance under low and optimum nitrogen regimes across environments using AMMI and GGE biplot analysis. Experiments were conducted using an augmented block design with three replications for two years in three locations in Nigeria. There were significant differences (p < 0.001) in the genotype's mean performances as well as significant differences (p < 0.001) in the environment's mean performances for all the traits measured in both nitrogen regimes. The AMMI analysis of variance showed significant effects (p < 0.001) for genotypes, environments and the interactions for fresh root yield in both nitrogen regimes. The biplot analysis showed that for fresh root yield in the optimum nitrogen regime, the principal component accounted for 81.54% of the G + GE (Genotype plus and Genotype by Environment) variation. The G + GE for fresh root yield in the low nitrogen regime accounted for a total of 71.64% of the variation. Ten genotypes were identified as the best genotypes under the optimum nitrogen regime, while eleven genotypes were the best under the low nitrogen regime. Three genotypes under optimum nitrogen regimes were high-yielding. Still, they were unstable in their fresh root yield performance across the environments and can be recommended as specifically adapted to the environments they performed best. Three other genotypes were high-yielding genotypes under low nitrogen but were highly unstable in their fresh root yield mean performance across the environments. The environments Otobi_YR1, Igbariam_YR2, and Umudike_YR1 were identified as the most discriminatory among the test environments. The environments Umudike_YR2 and Igbariam_YR1 were identified as the most representative of the test environments and can represent a mega-environment. The best 21 genotypes that performed above the grand mean for fresh root yield in both nitrogen regimes can be further evaluated on the farmer's field for possible advancement.

Editing of the MeSWEET10a promoter yields bacterial blight resistance in cassava cultivar SC8.

Cassava starch is a widely used raw material for industrial production and food source for people. However, cassava bacterial blight (CBB) caused by Xanthomonas axonopodis pv. manihotis (Xam) results in severe yield losses and is the most destructive bacterial disease in all worldwide cassava-growing regions. Xam11 is a highly pathogenic subspecies from China that infects the Chinese local cassava South China No. 8 (SC8) cultivar with marked symptoms. This study showed that the transcription activator-like effector TALE20Xam11 of Xam11 strain regulates the expression of disease-susceptibility gene MeSWEET10a by binding to the EBETALE20 region of the MeSWEET10a promoter in cassava cultivar SC8. CRISPR/Cas9-generated mutations of the EBETALE20 region resulted in a significant reduction in MeSWEET10a expression after infection by Xam11, correlating with reduced disease symptoms, smaller lesion sizes and decreased bacterial proliferation compared with the wild type. Importantly, the edited plants maintained normal growth, development and yield characteristics under greenhouse conditions. The results lay a research foundation for breeding resistant cassava cultivar SC8 to bacterial blight.

Molecular characterization of cassava zinc finger-homeodomain (ZF-HD) transcription factors reveals their role in disease resistance.

The production of cassava (Manihot esculenta Crantz) is constantly threatened by cassava bacterial blight (CBB), caused by Xanthomonas phaseoli pv. manihotis (Xpm). Zinc finger-homeodomain (ZF-HD) belongs to a family of homozygous heterotypic cassette genes widely implicated in various developmental and physiological processes in plants. Despite their importance, a comprehensive analysis of ZF-HD genes, particularly those involved in disease resistance, has not been performed for cassava. In the present study, we utilized bioinformatics methods to identify 21 ZF-HD genes distributed across 11 chromosomes of cassava genome, with the majority exhibiting gene structure without introns. Phylogenetic analysis categorized these genes into two major groups (MIF and ZHD) with five subgroups. We observed fourteen pairs of duplicated genes, suggesting that segmental duplication has likely facilitated the expansion of the cassava ZF-HD gene family. Comparative orthologous analyses between cassava and other plant species shed light on the evolutionary trajectory of this gene family. Promoter analyses revealed multiple hormone- and stress-related elements, indicative of a functional role in stress responses. Expression profiling through RNA-seq and quantitative real-time PCR (qRT-PCR) demonstrated that certain cassava ZF-HD genes are up-regulated in response to Xpm infection, suggesting their involvement in defense mechanisms. Notably, MeZHD7 gene was identified via virus induced gene silencing (VIGS) as potentially crucial in conferring resistance against CBB. Results from subcellular localization experiments indicated that MeZHD7 was localized in the nucleus. The Luciferase reporter assay demonstrated an interaction between MeZHD7 and MeMIF5. These findings may lay the foundation for further cloning and functional analyses of cassava ZF-HD genes, particularly those associated with pathogen resistance.

Tuber pathogenicity of Macrophomina phaseolina strain 3 a isolated from rotten cassava tuber from farm lands in Osogbo, Osun State, Nigeria, its virulence genes and ADMET properties.

BACKGROUND: Macrophomina phaseolina is a pathogen that causes an opportunistic disease that spreads by soil and seeds and affects more than 500 different plant species, like fruits, trees, and row crops. Mycotoxins, such as phaseolinic acid, and phaseolinone, are produced by M. phaseolina isolates in previous investigations; however, the production of these mycotoxins seems to vary depending on the host and the region. METHODS AND RESULTS: In this study, Macrophomina phaseolina strain 3 A was isolated from rotten cassava tuber and identified using the analysis of the sequences of the internal transcribed spacer region. The isolate was inoculated on a fresh healthy cassava tuber at 25 °C and tuber-rotting potential was monitored for 4 weeks. Virulence genes MPH_06603, MPH_06955, and MPH_01521 were determined with designed primers, and secondary metabolites were characterized by FTIR and GCMS. The rotten tuber effect was observed from the 2nd week of the experiment with severe tuber rot and weight reduction. The PCR showed the presence of MPH_06603 virulence gene. The GCMS showed N-Methylpivalamide (115.0 m/z), Butane, 1,4-dimethoxy- (119.0 m/z), and 5-Hydroxymethylfurfural (126.0 m/z) were the predominant metabolites produced by the pathogen. The compounds in the metabolites inhibit CYP3A4 enzymes, cause eye irritation, and Human Ether-a-go-go-related gene inhibition. CONCLUSION: This study revealed that M. phaseolina was responsible for the cassava tuber rot which leads to a lower yield of farm produce. The metabolites produced are toxic and unsafe for human consumption. It is suggested that farmers should destroy any cassava affected by this pathogen to prevent its toxic effects on humans and animals.

Promoter of Cassava MeAHL31 Responds to Diverse Abiotic Stresses and Hormone Signals in Transgenic Arabidopsis.

The AT-hook motif nuclear-localized (AHL) family is pivotal for the abiotic stress response in plants. However, the function of the cassava AHL genes has not been elucidated. Promoters, as important regulatory elements of gene expression, play a crucial role in stress resistance. In this study, the promoter of the cassava MeAHL31 gene was cloned. The MeAHL31 protein was localized to the cytoplasm and the nucleus. qRT-PCR analysis revealed that the MeAHL31 gene was expressed in almost all tissues tested, and the expression in tuber roots was 321.3 times higher than that in petioles. Promoter analysis showed that the MeAHL31 promoter contains drought, methyl jasmonate (MeJA), abscisic acid (ABA), and gibberellin (GA) cis-acting elements. Expression analysis indicated that the MeAHL31 gene is dramatically affected by treatments with salt, drought, MeJA, ABA, and GA3. Histochemical staining in the proMeAHL31-GUS transgenic Arabidopsis corroborated that the GUS staining was found in most tissues and organs, excluding seeds. Beta-glucuronidase (GUS) activity assays showed that the activities in the proMeAHL31-GUS transgenic Arabidopsis were enhanced by different concentrations of NaCl, mannitol (for simulating drought), and MeJA treatments. The integrated findings suggest that the MeAHL31 promoter responds to the abiotic stresses of salt and drought, and its activity is regulated by the MeJA hormone signal.

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.

Biochemical characterisation of a cassava (Manihot esculenta crantz) diversity panel for post-harvest physiological deterioration; metabolite involvement and environmental influence.

Cassava (Manihot esculenta Crantz) produces edible roots, a major carbohydrate source feeding more than 800 million people in Africa, Latin America, Oceania and Asia. Post-harvest physiological deterioration (PPD) renders harvested cassava roots unpalatable and unmarketable. Decades of research on PPD have elucidated several genetic, enzymatic and metabolic processes involved. Breeding populations were established to enable verification of robust biomarkers for PPD resistance. For comparison, these PPD populations have been cultivated concurrently with diversity population for carotenoid (ß-carotene) content. Results highlighted a significant variation of the chemotypes due to environmental factors. Less than 3% of the detected molecular features showed consistent trends between the two harvest years and were putatively identified as phenylpropanoid derived compounds (e.g. caffeoyl rutinoside). The data corroborated that ∼20 µg ß-carotene/g DW can reduced the PPD response of the cassava roots to a score of ∼1. Correlation analysis showed a significant correlation of ß-carotene content at harvest to PPD response (R2 -0.55). However, the decrease of ß-carotene over storage was not significantly correlated to initial content or PPD response. Volatile analysis observed changes of apocarotenoids derived from ß-carotene, lipid oxidation products (alkanes, alcohols and carbonyls and esters) and terpenes. The majority of these volatiles (>90%) showed no significant correlation to ß-carotene or PPD. Observed data indicated an increase (∼2-fold) of alkanes in varieties with ß-carotene >10 µg/g DW and a decrease (∼60%) in varieties with less ß-carotene. Fatty acid methyl esters with a chain length > C9 were detected solely after storage and show lower levels in varieties with higher ß-carotene content. In combination with correlation values to PPD (R2 ∼0.3; P-value >0.05), the data indicated a more efficient ROS quenching mechanism in PPD resistant varieties.

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.

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.

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.

Three cassava A20/AN1 family genes, Metip3 (5, and 7), can bestow on tolerance of plants to multiple abiotic stresses but show functional convergence and divergence.

A20/AN1 zinc-finger domain-containing genes are very promising candidates in improving plant tolerance to abiotic stresses, but considerably less is known about functions and mechanisms for many of them. In this study, Metip3 (5, and 7), cassava (Manihot esculenta) A20/AN1 genes carrying one A20 domain and one AN1 domain, were functionally characterized at different layers. Metip3 (5, and 7) proteins were all located in the nucleus. No interactions were found between these three proteins. Metip3 (5, and 7)-expressing Arabidopsis was more tolerant to multiple abiotic stresses by Na, Cd, Mn, Al, drought, high temperature, and low temperature. Metip3- and Metip5-expressing Arabidopsis was sensitive to Cu stress, while Metip7-expressing Arabidopsis was insensitive. The H2O2 production significantly decreased in all transgenic Arabidopsis, however, O2·- production significantly decreased in Metip3- and Metip5-expressing Arabidopsis but did not significantly changed in Metip7-expressing Arabidopsis under drought. Metip3 (5, and 7) expression-silenced cassava showed the decreased tolerance to drought and NaCl, presented significant decreases in superoxide dismutase and catalase activities and proline content, and displayed a significant increase in malondialdehyde content under drought. Taken together with transcriptome sequencing analysis, it is suggested that Metip5 gene can not only affect signal transduction related to plant hormone, mitogen activated protein kinases, and starch and sucrose metabolism, DRE-binding transcription factors, and antioxidants, conferring the drought tolerance, but also might deliver the signals from DREB2A INTERACTING PROTEIN1, E3 ubiquitin-protein ligases to proteasome, leading to the drought intolerance. The results are informative not only for further study on evolution of A20/AN1 genes but also for development of climate resilient crops.

Genome-wide association study for yield and quality of granulated cassava processed product.

The starchy storage roots of cassava are commonly processed into a variety of products, including cassava granulated processed products (gari). The commercial value of cassava roots depends on the yield and quality of processed products, directly influencing the acceptance of new varieties by farmers, processors, and consumers. This study aims to estimate genetic advance through phenotypic selection and identify genomic regions associated and candidate genes linked with gari yield and quality. Higher single nucleotide polymorphism (SNP)-based heritability estimates compared to broad-sense heritability estimates were observed for most traits highlighting the influence of genetic factors on observed variation. Using genome-wide association analysis of 188 clones, genotyped using 53,150 genome-wide SNPs, nine SNPs located on seven chromosomes were significantly associated with peel loss, gari yield, color parameters for gari and eba, bulk density, swelling index, and textural properties of eba. Future research will focus on validating and understanding the functions of identified genes and their influence on gari yield and quality traits.

Q&A: Methods for estimating genetic gain in sub-Saharan Africa and achieving improved gains.

Regular measurement of realized genetic gain allows plant breeders to assess and review the effectiveness of their strategies, allocate resources efficiently, and make informed decisions throughout the breeding process. Realized genetic gain estimation requires separating genetic trends from nongenetic trends using the linear mixed model (LMM) on historical multi-environment trial data. The LMM, accounting for the year effect, experimental designs, and heterogeneous residual variances, estimates best linear unbiased estimators of genotypes and regresses them on their years of origin. An illustrative example of estimating realized genetic gain was provided by analyzing historical data on fresh cassava (Manihot esculenta Crantz) yield in West Africa (https://github.com/Biometrics-IITA/Estimating-Realized-Genetic-Gain). This approach can serve as a model applicable to other crops and regions. Modernization of breeding programs is necessary to maximize the rate of genetic gain. This can be achieved by adopting genomics to enable faster breeding, accurate selection, and improved traits through genomic selection and gene editing. Tracking operational costs, establishing robust, digitalized data management and analytics systems, and developing effective varietal selection processes based on customer insights are also crucial for success. Capacity building and collaboration of breeding programs and institutions also play a significant role in accelerating genetic gains.

A Coiled-Coil Nucleotide-Binding Domain Leucine-Rich Repeat Receptor Gene MeRPPL1 Plays a Role in the Replication of a Geminivirus in Cassava.

Disease resistance gene (R gene)-encoded nucleotide-binding leucine-rich repeat proteins (NLRs) are critical players in plant host defence mechanisms because of their role as receptors that recognise pathogen effectors and trigger plant effector-triggered immunity (ETI). This study aimed to determine the putative role of a cassava coiled-coil (CC)-NLR (CNL) gene MeRPPL1 (Manes.12G091600) (single allele) located on chromosome 12 in the tolerance or susceptibility to South African cassava mosaic virus (SACMV), one of the causal agents of cassava mosaic disease (CMD). A transient protoplast system was used to knock down the expression of MeRPPL1 by clustered regularly interspaced short palindromic repeats-CRISPR-associated protein 9 (CRISPR-Cas9). The MeRPPL1-targeting CRISPR vectors and/or SACMV DNA A and DNA B infectious clones were used to transfect protoplasts isolated from leaf mesophyll cells from the SACMV-tolerant cassava (Manihot esculenta) cultivar TME3. The CRISPR/Cas9 silencing vector significantly reduced MeRPPL1 expression in protoplasts whether with or without SACMV co-infection. Notably, SACMV DNA A replication was higher in protoplasts with lower MeRPPL1 expression levels than in non-silenced protoplasts. Mutagenesis studies revealed that protoplast co-transfection with CRISPR-MeRPPL1 silencing vector + SACMV and transfection with only SACMV induced nucleotide substitution mutations that led to altered amino acids in the highly conserved MHD motif of the MeRPPL1-translated polypeptide. This may abolish or alter the regulatory role of the MHD motif in controlling R protein activity and could contribute to the increase in SACMV-DNA A accumulation observed in MeRPPL1-silenced protoplasts. The results herein demonstrate for the first time a role for a CNL gene in tolerance to a geminivirus in TME3.

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.

Comparative analysis of infected cassava root transcriptomics reveals candidate genes for root rot disease resistance.

Cassava root-rot incited by soil-borne pathogens is one of the major diseases that reduces root yield. Although the use of resistant cultivars is the most effective method of management, the genetic basis for root-rot resistance remains poorly understood. Therefore, our work analyzed the transcriptome of two contrasting genotypes (BRS Kiriris/resistant and BGM-1345/susceptible) using RNA-Seq to understand the molecular response and identify candidate genes for resistance. Cassava seedlings (resistant and susceptible to root-rot) were both planted in infested and sterilized soil and samples from Initial-time and Final-time periods, pooled. Two controls were used: (i) seedlings collected before planting in infested soil (absolute control) and, (ii) plants grown in sterilized soil (mock treatments). For the differentially expressed genes (DEGs) analysis 23.912 were expressed in the resistant genotype, where 10.307 were differentially expressed in the control treatment, 15 DEGs in the Initial Time-period and 366 DEGs in the Final Time-period. Eighteen candidate genes from the resistant genotype were related to plant defense, such as the MLP-like protein 31 and the peroxidase A2-like gene. This is the first model of resistance at the transcriptional level proposed for the cassava × root-rot pathosystem. Gene validation will contribute to screening for resistance of germplasm, segregating populations and/or use in gene editing in the pursuit to develop most promising cassava clones with resistance to root-rot.

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.

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.