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.

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.

Influence of Drought Stress on the Rhizosphere Bacterial Community Structure of Cassava (Manihot esculenta Crantz).

Drought presents a significant abiotic stress that threatens crop productivity worldwide. Rhizosphere bacteria play pivotal roles in modulating plant growth and resilience to environmental stresses. Despite this, the extent to which rhizosphere bacteria are instrumental in plant responses to drought, and whether distinct cassava (Manihot esculenta Crantz) varieties harbor specific rhizosphere bacterial assemblages, remains unclear. In this study, we measured the growth and physiological characteristics, as well as the physical and chemical properties of the rhizosphere soil of drought-tolerant (SC124) and drought-sensitive (SC8) cassava varieties under conditions of both well-watered and drought stress. Employing 16S rDNA high-throughput sequencing, we analyzed the composition and dynamics of the rhizosphere bacterial community. Under drought stress, biomass, plant height, stem diameter, quantum efficiency of photosystem II (Fv/Fm), and soluble sugar of cassava decreased for both SC8 and SC124. The two varieties' rhizosphere bacterial communities' overall taxonomic structure was highly similar, but there were slight differences in relative abundance. SC124 mainly relied on Gamma-proteobacteria and Acidobacteriae in response to drought stress, and the abundance of this class was positively correlated with soil acid phosphatase. SC8 mainly relied on Actinobacteria in response to drought stress, and the abundance of this class was positively correlated with soil urease and soil saccharase. Overall, this study confirmed the key role of drought-induced rhizosphere bacteria in improving the adaptation of cassava to drought stress and clarified that this process is significantly related to variety.

Understanding starch gelatinization and rheology modeling of tapioca starch- NaCl/CaCl2 blends: Thermodynamic properties and gelatinization reaction kinetics during pre- and post-ultrasonication.

The presence of salt can impact the fluid phase and gelatinization process of starch granules. The variation in viscosity and rheology models including the Herschel-Bulkley, the Casson model, and the power law, were determined by adding salts before and after starch ultrasonication. Non-isothermal kinetics can be utilized for the mathematical modeling of the gelatinization process and the evolution of the reaction. Unlike Na+ ions, Ca+2 ions notably elevate viscosity. The Casson model accurately predicts viscosity data. Results indicate that the addition of Na+ ions decreases yield stress by up to 60.4 %, while Ca+2 ions increase by up to 100.8 %. Adding Na+ ions decreases the required thermal energy by as much as 49.6 %, while the presence of Ca+2 ions can lead to a substantial increase of up to 337.1 % compared to control samples. The positive ∆G indicates a non-spontaneous gelatinization process. The addition of NaCl promotes a spontaneous reaction, while the addition of CaCl2 increases the Gibbs energy. The changes in entropy are minimal, implying minimal changes in starches' disorder structure.

Batch and semi-continuous treatment of cassava wastewater using microbial fuel cells and metataxonomic analysis.

The treatment of agroindustrial wastewater using microbial fuel cells (MFCs) is a technological strategy to harness its chemical energy while simultaneously purifying the water. This manuscript investigates the organic load effect as chemical oxygen demand (COD) on the production of electricity during the treatment of cassava wastewater by means of a dual-chamber microbial fuel cell in batch mode. Additionally, specific conditions were selected to evaluate the semi-continuous operational mode. The dynamics of microbial communities on the graphite anode were also investigated. The maximum power density delivered by the batch MFC (656.4 µW m - 2 ) was achieved at the highest evaluated organic load (6.8 g COD L - 1 ). Similarly, the largest COD removal efficiency (61.9%) was reached at the lowest organic load (1.17 g COD L - 1 ). Cyanide degradation percentages (50-70%) were achieved across treatments. The semi-continuous operation of the MFC for 2 months revealed that the voltage across the cell is dependent on the supply or suspension of the organic load feed. The electrode polarization resistance was observed to decreases over time, possibly due to the enrichment of the anode with electrogenic microbial communities. A metataxonomic analysis revealed a significant increase in bacteria from the phylum Firmicutes, primarily of the genus Enterococcus.

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.

Continuous rotary kiln pyrolysis of cassava plant shoot system and wide speciation of oxygenated and nitrogen-containing compounds in bio-oil by HESI and APPI-Orbitrap MS.

This work proposes the pyrolysis of the cassava plant shoot system biomass and a comprehensive chemical characterization of the resulting bio-oil. The highest yields of liquid products were obtained at 600 °C, with 12.6 % bio-oil (organic fraction), which presented the lowest total acid number of 65.7 mg KOH g-1. The bio-oil produced at 500 °C exhibited the highest total phenolic content of approximately 41 % GAE, confirmed by GC/MS analysis (33.8 % of the total area). FT-Orbitrap MS analysis found hundreds of oxygenated constituents in the bio-oils, belonging to the O2-7 classes, as well as nitrogen compounds from the Ny and OxNy classes. Higher pyrolysis temperatures resulted in more oxygenated phenolics (O4-7) undergoing secondary degradation and deoxygenation reactions, generating O2-3 compounds. Additional classes affected were O3-5N2-3, while O1-2N1 presented more stable compounds. These findings show that cassava bio-oils are promising sources of renewable chemicals.

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.

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.

Tapioca-starch-based bionanocomposites with fructose and titanium dioxide for food packaging and fertilization applications.

Bionanocomposites offer a promising solution to the plastic waste crisis. Although tapioca starch shows potential as a bioplastic material, it is characterized by low mechanical properties, poor thermal stability, and high water absorption owing to its hydrophilic nature. To increase the flexibility of the material and reduce the transmission rate of oxygen and water vapor, additives such as fructose and titanium dioxide (TiO2) can be incorporated into the material. TiO2 nanoparticles are commonly utilized in agriculture to enhance nutrient release and promote plant growth. In this study, X-ray diffraction analysis revealed that TiO2 reduced crystal size while increasing the crystallinity of bionanocomposites. Fourier-transform infrared spectroscopy analysis revealed an absorption peak at 3397 cm-1, indicating hydrogen bonding between TiO2 and starch-OH groups, and a peak at 773 cm-1, indicating an increase in the intensity of Ti-O-Ti stretching vibrations with the incorporation of TiO2. Water absorption rate results confirmed that TiO2 addition enhanced bionanocomposite resistance to water vapor and moisture, evidenced by increased tensile strength from 0.11 to 0.49 MPa and Young's modulus from 2.48 to 5.26 MPa, as well as decreased elongation at break from 21.46 % to 2.36 % in bionanocomposites with TiO2. Furthermore, with TiO2 addition, the biodegradation rate of the bionanocomposites decreased, which is beneficial for enhancing plant nutrient content.

Heavy metal accumulation in root and shoot tapioca plant biomass grown in agriculture land situated around the magnesite mine tailings.

Heavy metal pollution in soil has emerged as a major environmental concern. This can be attributed to human activities such as mining, modern agriculture, and industrialization. This study was conducted to determine how heavy metals spread from mine tailings to surrounding farmland. Metal absorption and accumulation were also investigated in the root and shoot biomass of tapioca crops grown in those farmlands. Metal concentrations in MTAS1 were 85.3 ± 1.2, 45.8 ± 1.5, 134.8 ± 1.7, 92.4 ± 2.2, and 78.95 ± 1.4 mg kg-1, respectively. Heavy metal concentrations in MTAS2 and MTAS3 were found to be 79.62 ± 1.6, 75.4 ± 1.5, 41.31 ± 1.1, 47.8 ± 1.6, 142.5 ± 2.1, 128.4 ± 1.4, 86.2 ± 1.9, 79.5 ± 1.3, and 83.4 ± 1.2 mg kg-1, respectively. Tapioca crop shoot and root biomass grown at these metal polluted sites absorbed and accumulated significant amounts of Cd, Cu, Zn, Pb, Ni, and Mn. Notably, the metal content of the tapioca crop's root and shoot biomass exceeded national standards.

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.

Current and future scenarios of suitability and expansion of cassava brown streak disease, Bemisia tabaci species complex, and cassava planting in Africa.

Cassava (Manihot esculenta) is among the most important staple crops globally, with an imperative role in supporting the Sustainable Development Goal of 'Zero hunger'. In sub-Saharan Africa, it is cultivated mainly by millions of subsistence farmers who depend directly on it for their socio-economic welfare. However, its yield in some regions has been threatened by several diseases, especially the cassava brown streak disease (CBSD). Changes in climatic conditions enhance the risk of the disease spreading to other planting regions. Here, we characterise the current and future distribution of cassava, CBSD and whitefly Bemisia tabaci species complex in Africa, using an ensemble of four species distribution models (SDMs): boosted regression trees, maximum entropy, generalised additive model, and multivariate adaptive regression splines, together with 28 environmental covariates. We collected 1,422 and 1,169 occurrence records for cassava and Bemisia tabaci species complex from the Global Biodiversity Information Facility and 750 CBSD occurrence records from published literature and systematic surveys in East Africa. Our results identified isothermality as having the highest contribution to the current distribution of cassava, while elevation was the top predictor of the current distribution of Bemisia tabaci species complex. Cassava harvested area and precipitation of the driest month contributed the most to explain the current distribution of CBSD outbreaks. The geographic distributions of these target species are also expected to shift under climate projection scenarios for two mid-century periods (2041-2060 and 2061-2080). Our results indicate that major cassava producers, like Cameron, Ivory Coast, Ghana, and Nigeria, are at greater risk of invasion of CBSD. These results highlight the need for firmer agricultural management and climate-change mitigation actions in Africa to combat new outbreaks and to contain the spread of CBSD.

Overexpression of cassava melatonin receptor PMTR1 plays dual roles in development under light and dark conditions in Arabidopsis.

KEY MESSAGE: MePMTR1 is involved in plant development and production as well as photosynthesis in plant. Melatonin is widely involved in plant growth and development as well as stress responses. Compared with the extending studies of melatonin in stress responses, the direct link between melatonin and plant development in the whole stages remains unclear. With the identification of phytomelatonin receptor PMTR1 in plants, melatonin signalling is becoming much clearer. However, the function of MePMTR1 in tropical crop cassava remains elusive. In this study, we found that overexpression of MePMTR1 showed larger biomass than wild type (WT), including higher number and area of leaves, weight, and accompanying with higher photosynthetic efficiency. Consistently, exogenous melatonin accelerated photosynthetic rate in Arabidopsis. In addition, MePMTR1-overexpressed plants exhibited more resistance to dark-induced senescence compared with WT, demonstrated by higher chlorophyll, lower hydrogen peroxide and superoxide content. In summary, this study illustrated that melatonin and its receptor regulate growth, development and senescence in plants, highlighting the potential application of melatonin and its receptor in improving crop yield and photosynthesis.

Single-tube colorimetric loop-mediated isothermal amplification (LAMP) assay for high-sensitivity detection of SLCMV in cassava from southern India.

Sri Lankan cassava mosaic virus (SLCMV) is a major cause for mosaic infections in cassava leaves, resulting in significant economic losses in southern India. SLCMV leads to growth retardation, leaf curl, and chlorosis in the host, with rapid transmission through whitefly insect vectors. Detecting SLCMV promptly is crucial, and the study introduces a novel and efficient colorimetric Loop-mediated isothermal amplification (LAMP) assay for successful detection in 60 min. Three primer sets were designed to target the conserved region of the SLCMV genome, specifically the coat protein gene, making the assay highly specific. The LAMP assay offers rapid and sensitive detection, completing within 60 min in a temperature-controlled water bath or thermal cycler. Compared to PCR techniques, it demonstrates 100 times superior sensitivity. The visual inspection of LAMP tube results using a nucleic acid dye and observing ladder-like pattern bands in a 2 % agarose gel confirms the presence of SLCMV. The assay is specific to SLCMV, showing no false positives or contaminations when tested against other virus. The standardized SLCMV LAMP assay proves technically efficient, providing a rapid, specific, simple, and low-cost solution, streamlining the detection and management of SLCMV.

Cassava brown streak virus evolves with a nucleotide-substitution rate that is typical for the family Potyviridae.

The ipomoviruses (family Potyviridae) that cause cassava brown streak disease (cassava brown streak virus [CBSV] and Uganda cassava brown streak virus [UCBSV]) are damaging plant pathogens that affect the sustainability of cassava production in East and Central Africa. However, little is known about the rate at which the viruses evolve and when they emerged in Africa - which inform how easily these viruses can host shift and resist RNAi approaches for control. We present here the rates of evolution determined from the coat protein gene (CP) of CBSV (Temporal signal in a UCBSV dataset was not sufficient for comparable analysis). Our BEAST analysis estimated the CBSV CP evolves at a mean rate of 1.43 × 10-3 nucleotide substitutions per site per year, with the most recent common ancestor of sampled CBSV isolates existing in 1944 (95% HPD, between years 1922 - 1963). We compared the published measured and estimated rates of evolution of CPs from ten families of plant viruses and showed that CBSV is an average-evolving potyvirid, but that members of Potyviridae evolve more quickly than members of Virgaviridae and the single representatives of Betaflexiviridae, Bunyaviridae, Caulimoviridae and Closteroviridae.

Key fuel characteristics and techno-economic aspects of torrefied rubberwood biomass pellets produced by incorporating various cassava-based binders at varied doses.

Improving energy content and hydrophobic nature of woody biomass can be pursued through torrefaction. This gives torrefied biomass with a low bulk density, potentially increasing storage and transport costs. To overcome this issue, densifying the torrefied biomass is necessary. However, poor binding of particles makes densification challenging without using a binder. Therefore, the aim of this study was to investigate the physicochemical characteristics and techno-economic aspects of torrefied rubberwood biomass (TRWB) when pelletized using various cassava-based binders at different blending ratios. The selected binders included cassava starch (CS), cassava pulp (CP), and cassava chip (CC). Each binder at 5%, 10%, or 15% (wt.) was mixed with TRWB and water before pelletizing using a flat die machine. The results revealed that pelletizing TRWB with different cassava-based binders at various blending ratios influenced the physicochemical characteristics of the TRWB pellets, particularly dimensions, bulk density, fuel and atomic ratios, and energy content. The TRWB pellets demonstrated energy densities in the range of 7.95-11.39 GJ/m3, and their mechanical durability and fine content fell within acceptable ranges. The TRWB pellets maintained their shape during 120 min of water soaking, with water absorption levels varying by binder dose. The pelletizing ability, material, and energy costs of TRWB pellets depend on binder type and dose. CP can be applied as a binder for pelletizing torrefied rubberwood biomass. However, the mechanical durability of the product needs to be above the user requirement or standard.

Functional Characterization of the MeSSIII-1 Gene and Its Promoter from Cassava.

Soluble starch synthases (SSs) play important roles in the synthesis of cassava starch. However, the expression characteristics of the cassava SSs genes have not been elucidated. In this study, the MeSSIII-1 gene and its promoter, from SC8 cassava cultivars, were respectively isolated by PCR amplification. MeSSIII-1 protein was localized to the chloroplasts. qRT-PCR analysis revealed that the MeSSIII-1 gene was expressed in almost all tissues tested, and the expression in mature leaves was 18.9 times more than that in tuber roots. MeSSIII-1 expression was induced by methyljasmonate (MeJA), abscisic acid (ABA), and ethylene (ET) hormones in cassava. MeSSIII-1 expression patterns were further confirmed in proMeSSIII-1 transgenic cassava. The promoter deletion analysis showed that the -264 bp to -1 bp MeSSIII-1 promoter has basal activity. The range from -1228 bp to -987 bp and -488 bp to -264 bp significantly enhance promoter activity. The regions from -987 bp to -747 bp and -747 bp to -488 bp have repressive activity. These findings will provide an important reference for research on the potential function and transcriptional regulation mechanisms of the MeSSIII-1 gene and for further in-depth exploration of the regulatory network of its internal functional elements.

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.