Resequencing of 388 cassava accessions identifies valuable loci and selection for variation in heterozygosity

Transcription factor AHL17 fine-tunes reactive oxygen species accumulation and cassava disease resistance

Zinc finger protein 67 interacts with a transcription factor to fine-tune reactive oxygen species accumulation and immune responses during pathogen infection in cassava.

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

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

TCOD: an integrated resource for tropical crops

Tropical crops are vital for tropical agriculture, with resource scarcity, functional diversity and extensive market demand, providing considerable economic benefits for the world's tropical agriculture-producing countries. The rapid development of sequencing technology has promoted a milestone in tropical crop research, resulting in the generation of massive amount of data, which urgently needs an effective platform for data integration and sharing. However, the existing databases cannot fully satisfy researchers' requirements due to the relatively limited integration level and untimely update. Here, we present the Tropical Crop Omics Database (TCOD, https://ngdc.cncb.ac.cn/tcod), a comprehensive multi-omics data platform for tropical crops. TCOD integrates diverse omics data from 15 species, encompassing 34 chromosome-level de novo assemblies, 1 255 004 genes with functional annotations, 282 436 992 unique variants from 2048 WGS samples, 88 transcriptomic profiles from 1997 RNA-Seq samples and 13 381 germplasm items. Additionally, TCOD not only employs genes as a bridge to interconnect multi-omics data, enabling cross-species comparisons based on homology relationships, but also offers user-friendly online tools for efficient data mining and visualization. In short, TCOD integrates multi-species, multi-omics data and online tools, which will facilitate the research on genomic selective breeding and trait biology of tropical crops.

MePP2C24, a cassava (Manihot esculenta) gene encoding protein phosphatase 2C, negatively regulates drought stress and abscisic acid responses in transgenic Arabidopsis thaliana

Abscisic acid (ABA) signaling plays a crucial role in plant development and response to abiotic/biotic stress. However, the function and regulation of protein phosphatase 2C (PP2C), a key component of abscisic acid signaling, under abiotic stress are still unknown in cassava, a drought-tolerant crop. In this study, a cassava PP2C gene (MePP2C24) was cloned and characterized. The MePP2C24 transcripts increased in response to mannitol, NaCl, and ABA. Overexpression of MePP2C24 in Arabidopsis resulted in increased sensitivity to drought stress and decreased sensitivity to exogenous ABA. This was demonstrated by transgenic lines having higher levels of malondialdehyde (MDA), ion leakage (IL), and reactive oxygen species (ROS), lower activities of catalase (CAT) and peroxidase (POD), and lower proline content than wild type (WT) under drought stress. Moreover, MePP2C24 overexpression caused decrease in expression of drought-responsive genes related to ABA signaling pathway. In addition, MePP2C24 was localized in the cell nucleus and showed self-activation. Furthermore, many MePYLs (MePYL1, MePYL4, MePYL7-9, and MePYL11-13) could interact with MePP2C24 in the presence of ABA, and MePYL1 interacted with MePP2C24 in both the presence and absence of ABA. Additionally, MebZIP11 interacted with the promoter of MePP2C24 and exerted a suppressive effect. Taken together, our results suggest that MePP2C24 acts as a negative regulator of drought tolerance and ABA response.

Metabolic GWAS-based dissection of genetic basis underlying nutrient quality variation and domestication of cassava storage root

BACKGROUND

Metabolites play critical roles in regulating nutritional qualities of plants, thereby influencing their consumption and human health. However, the genetic basis underlying the metabolite-based nutrient quality and domestication of root and tuber crops remain largely unknown.

RESULTS

We report a comprehensive study combining metabolic and phenotypic genome-wide association studies to dissect the genetic basis of metabolites in the storage root (SR) of cassava. We quantify 2,980 metabolic features in 299 cultivated cassava accessions. We detect 18,218 significant marker-metabolite associations via metabolic genome-wide association mapping and identify 12 candidate genes responsible for the levels of metabolites that are of potential nutritional importance. Me3GT, MeMYB4, and UGT85K4/UGT85K5, which are involved in flavone, anthocyanin, and cyanogenic glucoside metabolism, respectively, are functionally validated through in vitro enzyme assays and in vivo gene silencing analyses. We identify a cluster of cyanogenic glucoside biosynthesis genes, among which CYP79D1, CYP71E7b, and UGT85K5 are highly co-expressed and their allelic combination contributes to low linamarin content. We find MeMYB4 is responsible for variations in cyanidin 3-O-glucoside and delphinidin 3-O-rutinoside contents, thus controlling SR endothelium color. We find human selection affects quercetin 3-O-glucoside content and SR weight per plant. The candidate gene MeFLS1 is subject to selection during cassava domestication, leading to decreased quercetin 3-O-glucoside content and thus increased SR weight per plant.

CONCLUSIONS

These findings reveal the genetic basis of cassava SR metabolome variation, establish a linkage between metabolites and agronomic traits, and offer useful resources for genetically improving the nutrition of cassava and other root crops.

Plant genome resequencing and population genomics: Current status and future prospects

Advances in DNA sequencing technology have sparked a genomics revolution, driving breakthroughs in plant genetics and crop breeding. Recently, the focus has shifted from cataloging genetic diversity in plants to exploring their functional significance and delivering beneficial alleles for crop improvement. This transformation has been facilitated by the increasing adoption of whole-genome resequencing. In this review, we summarize the current progress of population-based genome resequencing studies and how these studies affect crop breeding. A total of 187 land plants from 163 countries have been resequenced, comprising 54 413 accessions. As part of resequencing efforts 367 traits have been surveyed and 86 genome-wide association studies have been conducted. Economically important crops, particularly cereals, vegetables, and legumes, have dominated the resequencing efforts, leaving a gap in 49 orders, including Lycopodiales, Liliales, Acorales, Austrobaileyales, and Commelinales. The resequenced germplasm is distributed across diverse geographic locations, providing a global perspective on plant genomics. We highlight genes that have been selected during domestication, or associated with agronomic traits, and form a repository of candidate genes for future research and application. Despite the opportunities for cross-species comparative genomics, many population genomic datasets are not accessible, impeding secondary analyses. We call for a more open and collaborative approach to population genomics that promotes data sharing and encourages contribution-based credit policy. The number of plant genome resequencing studies will continue to rise with the decreasing DNA sequencing costs, coupled with advances in analysis and computational technologies. This expansion, in terms of both scale and quality, holds promise for deeper insights into plant trait genetics and breeding design.

Genomic analysis unveils reduced genetic variability but increased proportion of heterozygotic genotypes of the intensively managed mezcal agave, Agave angustifolia

PREMISE

The central Oaxaca Basin has a century-long history of agave cultivation and is hypothesized to be the region of origin of other cultivated crops. Widely cultivated for mezcal production, the perennial crop known as "espadín" is putatively derived from wild Agave angustifolia. Nevertheless, little is known about its genetic relationship to the wild A. angustifolia or how the decades-long clonal propagation has affected its genetics.

METHODS

Using restriction-site-associated DNA sequencing and over 8000 single-nucleotide polymorphisms, we studied aspects of the population genomics of wild and cultivated A. angustifolia in Puebla and Oaxaca, Mexico. We assessed patterns of genetic diversity, inbreeding, distribution of genetic variation, and differentiation among and within wild populations and plantations.

RESULTS

Genetic differentiation between wild and cultivated plants was strong, and both gene pools harbored multiple unique alleles. Nevertheless, we found several cultivated individuals with high genetic affinity with wild samples. Higher heterozygosity was observed in the cultivated individuals, while in total, they harbored considerably fewer alleles and presented higher linkage disequilibrium compared to the wild plants. Independently of geographic distance among sampled plantations, the genetic relatedness of the cultivated plants was high, suggesting a common origin and prevalent role of clonal propagation.

CONCLUSIONS

The considerable heterozygosity found in espadín is contained within a network of highly related individuals, displaying high linkage disequilibrium generated by decades of clonal propagation and possibly by the accumulation of somatic mutations. Wild A. angustifolia, on the other hand, represents a significant genetic diversity reservoir that should be carefully studied and conserved.

Analysis of the molecular and biochemical mechanisms involved in the symbiotic relationship between Arbuscular mycorrhiza fungi and Manihot esculenta Crantz

INTRODUCTION

Plants and arbuscular mycorrhizal fungi (AMF) mutualistic interactions are essential for sustainable agriculture production. Although it is shown that AMF inoculation improves cassava physiological performances and yield traits, the molecular mechanisms involved in AM symbiosis remain largely unknown. Herein, we integrated metabolomics and transcriptomics analyses of symbiotic (Ri) and asymbiotic (CK) cassava roots and explored AM-induced biochemical and transcriptional changes.

RESULTS

Three weeks (3w) after AMF inoculations, proliferating fungal hyphae were observable, and plant height and root length were significantly increased. In total, we identified 1,016 metabolites, of which 25 were differentially accumulated (DAMs) at 3w. The most highly induced metabolites were 5-aminolevulinic acid, L-glutamic acid, and lysoPC 18:2. Transcriptome analysis identified 693 and 6,481 differentially expressed genes (DEGs) in the comparison between CK (3w) against Ri at 3w and 6w, respectively. Functional enrichment analyses of DAMs and DEGs unveiled transport, amino acids and sugar metabolisms, biosynthesis of secondary metabolites, plant hormone signal transduction, phenylpropanoid biosynthesis, and plant-pathogen interactions as the most differentially regulated pathways. Potential candidate genes, including nitrogen and phosphate transporters, transcription factors, phytohormone, sugar metabolism-related, and SYM (symbiosis) signaling pathway-related, were identified for future functional studies.

DISCUSSION

Our results provide molecular insights into AM symbiosis and valuable resources for improving cassava production.

In Silico Identification of Cassava Genome-Encoded MicroRNAs with Predicted Potential for Targeting the ICMV-Kerala Begomoviral Pathogen of Cassava

Cassava mosaic disease (CMD) is caused by several divergent species belonging to the genus Begomovirus (Geminiviridae) transmitted by the whitefly Bemisia tabaci cryptic species group. In India and other parts of Asia, the Indian cassava mosaic virus-Kerala (ICMV-Ker) is an emergent begomovirus of cassava causing damage that results in reduced yield loss and tuber quality. Double-stranded RNA-mediated interference (RNAi) is an evolutionary conserved mechanism in eukaryotes and highly effective, innate defense system to inhibit plant viral replication and/or translation. The objective of this study was to identify and characterize cassava genome-encoded microRNAs (mes-miRNA) that are predicted to target ICMV-Ker ssDNA-encoded mRNAs, based on four in silico algorithms: miRanda, RNA22, Tapirhybrid, and psRNA. The goal is to deploy the predicted miRNAs to trigger RNAi and develop cassava plants with resistance to ICMV-Ker. Experimentally validated mature cassava miRNA sequences (n = 175) were downloaded from the miRBase biological database and aligned with the ICMV-Ker genome. The miRNAs were evaluated for base-pairing with the cassava miRNA seed regions and to complementary binding sites within target viral mRNAs. Among the 175 locus-derived mes-miRNAs evaluated, one cassava miRNA homolog, mes-miR1446a, was identified to have a predicted miRNA target binding site, at position 2053 of the ICMV-Ker genome. To predict whether the cassava miRNA might bind predicted ICMV-Ker mRNA target(s) that could disrupt viral infection of cassava plants, a cassava locus-derived miRNA-mRNA regulatory network was constructed using Circos software. The in silico-predicted cassava locus-derived mes-miRNA-mRNA network corroborated interactions between cassava mature miRNAs and the ICMV-Ker genome that warrant in vivo analysis, which could lead to the development of ICMV-Ker resistant cassava plants.

Genome-wide association studies reveal genetic basis of ionomic variation in cassava

As one of the most important food crops, cassava (Manihot esculenta) is the main dietary source of micronutrients for about 1 billion people. However, the ionomic variation in cassava and the underlying genetic mechanisms remain unclear so far. Herein, genome-wide association studies were performed to reveal the specific single nucleotide polymorphisms (SNPs) that affect the ionomic variation in cassava. We identified 164 SNPs with P-values lower than the threshold located in 88 loci associated with divergent ionomic variations. Among them, 13 SNPs are related to both calcium (Ca) and magnesium (Mg), and many loci for different ionomic traits seem to be clustered on specific chromosome regions. Moreover, we identified the peak SNPs in the promoter regions of Sc10g003170 (encoding methionyl-tRNA synthetase [MetRS]) and Sc18g015190 (encoding the transcriptional regulatory protein AlgP) for nitrogen (N) and phosphorus (P) accumulation, respectively. Notably, these two SNPs (chr10_32807962 and chr18_31343738) were directly correlated with the transcript levels of Sc10g003170 (MetRS) and Sc18g015190 (AlgP), which positively modulated N accumulation and P concentration in cassava, respectively. Taken together, this study provides important insight into the genetic basis of cassava natural ionomic variation, which will promote genetic breeding to improve nutrient use and accumulation of elements in cassava.

PP2C1 fine-tunes melatonin biosynthesis and phytomelatonin receptor PMTR1 binding to melatonin in cassava

Melatonin is an important molecule in both animals and plants, regulating circadian rhythms and stress responses. Therefore, the improvement of melatonin accumulation not only strengthens the function of melatonin but also improves stress resistance in crops. Although melatonin biosynthetic enzymes have been identified through reverse genetics previously, an investigation of melatonin level-related genes through forward genetics in plants has yet to be performed. In this study, a genome-wide association study using cassava natural population of 298 genetic resources identified melatonin accumulation 1 (MA1), which regulates the natural variation of melatonin levels in cassava. We found that MA1 encodes type 2C protein phosphatase 1 (PP2C1), which serves as a negative regulator of melatonin levels in cassava. MePP2C1 physically interacts with MeRAV1/2 and MeWRKY20 and dephosphorylates them at serine (S) 35 residue, S34 residue, and S176 residue, respectively, thereby hindering their transcriptional activation on downstream melatonin biosynthetic genes. Notably, MePP2C1 interacts with phytomelatonin receptor MePMTR1 and dephosphorylates it at S11 residue, repressing its binding to melatonin. In summary, this study demonstrates that MePP2C1 as MA1 plays dual roles in negatively regulating both melatonin accumulation and signaling, extending the understanding of the molecular mechanism underlying melatonin accumulation and signaling through forward genetics in plants.