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.

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.

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.

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.

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.

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.

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.

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.

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.

Bio-synthesis and characterization of silver nanoparticles from Trichoderma species against cassava root rot disease.

Cassava root rot disease caused by the fungal pathogens Fusarium solani and Lasiodiplodia theobromae produces severe damages on cassava production. This research was conducted to produce and assess silver nanoparticles (AgNPs) synthesized by Trichoderma harzianum for reducing root rot disease. The results revealed that using the supernatants of T. harzianum on a silver nitrate solution changed it to reddish color at 48 h, indicating the formation of AgNPs. Further characterization was identified using dynamic light scattering (DLS) and scanning electron microscope (SEM). DLS supported that the Z-average size is at 39.79 nm and the mean zeta potential is at - 36.5 mV. SEM revealed the formation of monodispersed spherical shape with a diameter between 60-75 nm. The antibacterial action of AgNPs as an antifungal agent was demonstrated by an observed decrease in the size of the fungal colonies using an increasing concentration of AgNPs until the complete inhibition growth of L. theobromae and F. solani at > 58 µg mL-1 and at ≥ 50 µg mL-1, respectively. At in vitro conditions, the applied AgNPs caused a decrease in the percentage of healthy aerial hyphae of L. theobromae (32.5%) and of F. solani (70.0%) compared to control (100%). The SR-FTIR spectra showed the highest peaks in the first region (3000-2800 cm-1) associated with lipids and fatty acids located at 2962, 2927, and 2854 cm-1 in the AgNPs treated samples. The second region (1700-1450 cm-1) consisting of proteins and peptides revealed the highest peaks at 1658, 1641, and 1548 cm-1 in the AgNPs treated samples. The third region (1300-900 cm-1), which involves nucleic acid, phospholipids, polysaccharides, and carbohydrates, revealed the highest peaks at 1155, 1079, and 1027 cm-1 in the readings from the untreated samples. Finally, the observed root rot severity on cassava roots treated with AgNPs (1.75 ± 0.50) was significantly lower than the control samples (5.00 ± 0.00).

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.

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.

Real-time reverse transcription recombinase polymerase amplification (RT-RPA) assay for detection of cassava brown streak viruses.

Cassava brown streak disease (CBSD) caused by Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) is the most economically important viral disease of cassava. As cassava is a vegetatively propagated crop, the development of rapid and sensitive diagnostics would aid in the identification of virus-free planting material and development of effective management strategies. In this study, a rapid, specific and sensitive real-time reverse transcription recombinase polymerase amplification (RT-RPA) assay was developed for real-time detection of CBSV and UCBSV. The RT-RPA was able to detect as little as 2 pg/µl of purified RNA obtained from infected cassava leaves, a sensitivity equivalent to that obtained by quantitative real-time reverse transcription PCR (qRT-PCR), within 20 min at 37 °C. Further, the RT-RPA detected each target virus directly from crude leaf and stem extracts, avoiding the tedious and costly isolation of high-quality RNA. The developed RT-RPA assay provides a valuable diagnostic tool that can be adopted by cassava seed certification and virus resistance breeding programs to ensure distribution of virus-free cassava planting materials to farmers. This is the first report on the development and validation of crude sap-based RT-RPA assay for the detection of cassava brown streak viruses (UCBSV and CBSV) infection in cassava plants.

Proline inhibits postharvest physiological deterioration of cassava by improving antioxidant capacity.

Cassava (Manihot esculenta Crantz), a crucial global tuber crop, encounters significant economic losses attributed to postharvest physiological deterioration (PPD). The PPD phenomenon in cassava is closely related to the accumulation of reactive oxygen species (ROS), and amino acids play a pivotal role in regulating signaling pathways and eliminating ROS. In this study, the storage performance of eight cassava varieties were conducted. Cassava cultivar SC5 showed the best storage performance among the eight cassava varieties, but the edible cassava cultivar SC9 performed much worse. Comparative analysis of free amino acids was conducted in eight cassava varieties, revealing changes in proline, aspartic acid, histidine, glutamic acid, threonine, and serine. Exogenous supplementation of these six amino acids was performed to inhibit PPD of SC9. Proline was confirmed as the key amino acid for inhibiting PPD. Treatment with optimal exogenous proline of 5 g/L resulted in a 17.9% decrease in the deterioration rate compared to untreated cassava. Accompanied by a decrease in H2O2 content and an increase in catalase, superoxide dismutase and ascorbate peroxidase activity. Proline treatment proved to be an effective approach to alleviate cell oxidative damage, inhibit PPD in cassava, and prolong shelf life.

Micronization induced gelatinization of tapioca starch and its effects on starch physicochemical and structural properties.

The vibrating superfine mill (VSM) is a machine that belongs to the micronization technique. In this study, VSM was employed to produce micronized tapioca starch by varying micronization times (15, 30, 45, and 60 min). The structural and physicochemical properties of the micronized starch were then examined. Scanning electron microscopy studies revealed that micronized starch was partially gelatinized, and the granule size dramatically increased when micronization time increased. X-ray diffraction patterns showed that the relative crystallinity was decreased from 24.67% (native) to 4.13% after micronization treatment for 15 min and slightly decreased after that. The solubility of micronized starch significantly increased as the micronization time increased, which was associated with the destruction of the starch crystalline structure. Differential scanning calorimetry investigations confirmed that micronized starch was "partly gelatinized," and the degree of gelatinization increased to 81.27% when the micronization time was 60 min. The weight-average molar mass was reduced by 15.0% (15 min), 30.9% (30 min), 55.7% (45 min), and 70.5% (60 min), respectively, indicating that the molecular structure was seriously degraded. The results demonstrated that the physicochemical changes of micronized starch granules were related to the destruction of the starch structure. These observations would provide details on micronized starch and its potential applications. PRACTICAL APPLICATION: These observations would provide details on micronized starch and its potential applications. Moreover, we believe that when the structures of starches were known, it is probable that the effect of VSM on the structural and physicochemical properties change of other starches might be predicted by adjusting the processing time.

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.