Cassava genome from a wild ancestor to cultivated varieties

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

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

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.

Distinct heat response molecular mechanisms emerge in cassava vasculature compared to leaf mesophyll tissue under high temperature stress

With growing concerns over global warming, cultivating heat-tolerant crops has become paramount to prepare for the anticipated warmer climate. Cassava (Manihot esculenta Crantz), a vital tropical crop, demonstrates exceptional growth and productivity under high-temperature (HT) conditions. Yet, studies elucidating HT resistance mechanisms in cassava, particularly within vascular tissues, are rare. We dissected the leaf mid-vein from leaf, and did the comparative transcriptome profiling between mid-vein and leaf to figure out the cassava vasculature HT resistance molecular mechanism. Anatomical microscopy revealed that cassava leaf veins predominantly consisted of vasculature. A thermal imaging analysis indicated that cassava experienced elevated temperatures, coinciding with a reduction in photosynthesis. Transcriptome sequencing produced clean reads in total of 89.17G. Using Venn enrichment, there were 65 differentially expressed genes (DEGs) and 93 DEGs had been found highly specifically expressed in leaf and mid-vein. Further investigation disclosed that leaves enhanced pyruvate synthesis as a strategy to withstand high temperatures, while mid-veins fortified themselves by bolstering lignin synthesis by comprehensive GO and KEGG analysis of DEGs. The identified genes in these metabolic pathways were corroborated through quantity PCR (QPCR), with results aligning with the transcriptomic data. To verify the expression localization of DEGs, we used in situ hybridization experiments to identify the expression of MeCCoAMT(caffeoyl-coenzyme A-3-O-methyltransferase) in the lignin synthesis pathway in cassava leaf veins xylem. These findings unravel the disparate thermotolerance mechanisms exhibited by cassava leaves and mid-veins, offering insights that could potentially inform strategies for enhancing thermotolerance in other crops.

Microbiome and plant cell transformation trigger insect gall induction in cassava

Several specialised insects can manipulate normal plant development to induce a highly organised structure known as a gall, which represents one of the most complex interactions between insects and plants. Thus far, the mechanism for insect-induced plant galls has remained elusive. To study the induction mechanism of insect galls, we selected the gall induced by Iatrophobia brasiliensis (Diptera: Cecidomyiidae) in cassava (Euphorbiaceae: Manihot esculenta Crantz) as our model. PCR-based molecular markers and deep metagenomic sequencing data were employed to analyse the gall microbiome and to test the hypothesis that gall cells are genetically transformed by insect vectored bacteria. A shotgun sequencing discrimination approach was implemented to selectively discriminate between foreign DNA and the reference host plant genome. Several known candidate insertion sequences were identified, the most significant being DNA sequences found in bacterial genes related to the transcription regulatory factor CadR, cadmium-transporting ATPase encoded by the cadA gene, nitrate transport permease protein (nrtB gene), and arsenical pump ATPase (arsA gene). In addition, a DNA fragment associated with ubiquitin-like gene E2 was identified as a potential accessory genetic element involved in gall induction mechanism. Furthermore, our results suggest that the increased quality and rapid development of gall tissue are mostly driven by microbiome enrichment and the acquisition of critical endophytes. An initial gall-like structure was experimentally obtained in M. esculenta cultured tissues through inoculation assays using a Rhodococcus bacterial strain that originated from the inducing insect, which we related to the gall induction process. We provide evidence that the modification of the endophytic microbiome and the genetic transformation of plant cells in M. esculenta are two essential requirements for insect-induced gall formation. Based on these findings and having observed the same potential DNA marker in galls from other plant species (ubiquitin-like gene E2), we speculate that bacterially mediated genetic transformation of plant cells may represent a more widespread gall induction mechanism found in nature.

The role of orphan crops in the transition to nutritional quality-oriented crop improvement

Micronutrient malnutrition is a persisting problem threatening global human health. Biofortification via metabolic engineering has been proposed as a cost-effective and short-term means to alleviate this burden. There has been a recent rise in the recognition of potential that underutilized, orphan crops can hold in decreasing malnutrition concerns. Here, we illustrate how orphan crops can serve as a medium to provide micronutrients to populations in need, whilst promoting and maintaining dietary diversity. We provide a roadmap, illustrating which aspects to be taken into consideration when evaluating orphan crops. Recent developments have shown successful biofortification via metabolic engineering in staple crops. This review provides guidance in the implementation of these successes to relevant orphan crop species, with a specific focus on the relevant micronutrients iron, zinc, provitamin A and folates.

Isolation and Identification of Endophytic Bacteria Bacillus sp. ME9 That Exhibits Biocontrol Activity against Xanthomonas phaseoli pv. manihotis

In recent years, the bacterial blight of cassava has caused substantial economic losses to the Chinese cassava industry. Chemical control methods have become the primary approach to control this disease; however, their widespread usage and harmful residues have raised concerns about environmental pollution. In order to avoid this, it is urgent to seek a green ecological method to prevent and control it. Biological control through the utilization of microorganisms not only effectively inhibits the disease, but also gives consideration to environmental friendliness. Therefore, investigating an endophytic biological control method for cassava bacterial blight is of great importance. In this study, cassava leaf tissues were used as test specimens in order to isolate endophytic bacteria by using dilution and separation methods. Bacillus ME9, derived from cassava endophytic bacteria, exhibits good antagonism against a diverse range of pathogens, including Xpm11. Its genome consists of a series of genes encoding antibacterial lipopeptides, which may be directly related to its antibacterial capabilities. Furthermore, inoculation resulted in a substantial change in the diversity of the endophytic bacterial community, characterized by improved diversity, and displayed an obvious inhibition of pathogenic bacterial growth, demonstrating successful colonization within plants. The results laid a foundation and provided theoretical support for the development and utilization of cassava endophytic bacterial diversity and endogenous disease control strategies.

Phylogenetic Analysis of the PR-4 Gene Family in Euphorbiaceae and Its Expression Profiles in Tung Tree (Vernicia fordii)

Pathogenesis-related protein-4 (PR-4) is generally believed to be involved in physiological processes. However, a comprehensive investigation of this protein in tung tree (Vernicia fordii) has yet to be conducted. In this study, we identified 30 PR-4 genes in the genomes of Euphorbiaceae species and investigated their domain organization, evolution, promoter cis-elements, expression profiles, and expression profiles in the tung tree. Sequence and structural analyses indicated that VF16136 and VF16135 in the tung tree could be classified as belonging to Class II and I, respectively. Phylogenetic and Ka/Ks analyses revealed that Hevea brasiliensis exhibited a significantly expanded number of PR-4 genes. Additionally, the analysis of promoter cis-elements suggested that two VfPR-4 genes may play a role in the response to hormones and biotic and abiotic stress of tung trees. Furthermore, the expression patterns of VfPR-4 genes and their responses to 6-BA, salicylic acid, and silver nitrate in inflorescence buds of tung trees were evaluated using qRT-PCR. Notably, the expression of two VfPR-4 genes was found to be particularly high in leaves and early stages of tung seeds. These results suggest that VF16136 and VF16135 may have significant roles in the development of leaves and seeds in tung trees. Furthermore, these genes were found to be responsive to 6-BA, salicylic acid, and silver nitrate in the development of inflorescence buds. This research provides valuable insights for future investigation into the functions of PR-4 genes in tung trees.

Development of cassava common mosaic virus-based vector for protein expression and gene editing in cassava

BACKGROUND

Plant virus vectors designed for virus-mediated protein overexpression (VOX), virus-induced gene silencing (VIGS), and genome editing (VIGE) provide rapid and cost-effective tools for functional genomics studies, biotechnology applications and genome modification in plants. We previously reported that a cassava common mosaic virus (CsCMV, genus Potexvirus)-based VIGS vector was used for rapid gene function analysis in cassava. However, there are no VOX and VIGE vectors available in cassava.

RESULTS

In this study, we developed an efficient VOX vector (CsCMV2-NC) for cassava by modifying the CsCMV-based VIGS vector. Specifically, the length of the duplicated putative subgenomic promoter (SGP1) of the CsCMV CP gene was increased to improve heterologous protein expression in cassava plants. The modified CsCMV2-NC-based VOX vector was engineered to express genes encoding green fluorescent protein (GFP), bacterial phytoene synthase (crtB), and Xanthomonas axonopodis pv. manihotis (Xam) type III effector XopAO1 for viral infection tracking, carotenoid biofortification and Xam virulence effector identification in cassava. In addition, we used CsCMV2-NC to deliver single guide RNAs (gMePDS1/2) targeting two loci of the cassava phytoene desaturase gene (MePDS) in Cas9-overexpressing transgenic cassava lines. The CsCMV-gMePDS1/2 efficiently induced deletion mutations of the targeted MePDS with the albino phenotypes in systemically infected cassava leaves.

CONCLUSIONS

Our results provide a useful tool for rapid and efficient heterologous protein expression and guide RNA delivery in cassava. This expands the potential applications of CsCMV-based vector in gene function studies, biotechnology research, and precision breeding for cassava.

Fully sequencing the cassava full-length cDNA library reveals unannotated transcript structures and alternative splicing events in regions with a high density of single nucleotide variations, insertions-deletions, and heterozygous sequences

The primary transcript structure provides critical insights into protein diversity, transcriptional modification, and functions. Cassava transcript structures are highly diverse because of alternative splicing (AS) events and high heterozygosity. To precisely determine and characterize transcript structures, fully sequencing cloned transcripts is the most reliable method. However, cassava annotations were mainly determined according to fragmentation-based sequencing analyses (e.g., EST and short-read RNA-seq). In this study, we sequenced the cassava full-length cDNA library, which included rare transcripts. We obtained 8,628 non-redundant fully sequenced transcripts and detected 615 unannotated AS events and 421 unannotated loci. The different protein sequences resulting from the unannotated AS events tended to have diverse functional domains, implying that unannotated AS contributes to the truncation of functional domains. The unannotated loci tended to be derived from orphan genes, implying that the loci may be associated with cassava-specific traits. Unexpectedly, individual cassava transcripts were more likely to have multiple AS events than Arabidopsis transcripts, suggestive of the regulated interactions between cassava splicing-related complexes. We also observed that the unannotated loci and/or AS events were commonly in regions with abundant single nucleotide variations, insertions-deletions, and heterozygous sequences. These findings reflect the utility of completely sequenced FLcDNA clones for overcoming cassava-specific annotation-related problems to elucidate transcript structures. Our work provides researchers with transcript structural details that are useful for annotating highly diverse and unique transcripts and alternative splicing events.

Genome-Wide Identification and Expression Analysis of the SHI-Related Sequence Family in Cassava

The SHORT INTERNODES (SHI)-related sequences (SRS) are plant-specific transcription factors that have been quantitatively characterized during plant growth, regeneration, and stress responses. However, the genome-wide discovery of SRS family genes and their involvement in abiotic stress-related activities in cassava have not been documented. A genome-wide search strategy was used to identify eight family members of the SRS gene family in cassava (Manihot esculenta Crantz). Based on their evolutionary linkages, all MeSRS genes featured homologous RING-like zinc finger and IXGH domains. Genetic architecture and conserved motif analysis validated the categorization of MeSRS genes into four groups. Eight pairs of segmental duplications were detected, resulting in an increase in the number of MeSRS genes. Orthologous studies of SRS genes among cassava and three different plant species (Arabidopsis thaliana, Oryza sativa, and Populus trichocarpa) provided important insights into the probable history of the MeSRS gene family. The functionality of MeSRS genes was elucidated through the prediction of protein–protein interaction networks and cis-acting domains. RNA-seq data demonstrated tissue/organ expression selectivity and preference of the MeSRS genes. Furthermore, qRT-PCR investigation of MeSRS gene expression after exposure to salicylic acid (SA) and methyl jasmonate (MeJA) hormone treatments, as well as salt (NaCl) and osmotic (polyethylene glycol, PEG) stresses, showed their stress-responsive patterns. This genome-wide characterization and identification of the evolutionary relationships and expression profiles of the cassava MeSRS family genes will be helpful for further research into this gene family and its function in stress response. It may also assist future agricultural efforts to increase the stress tolerance of cassava.

NetMoST: A network-based machine learning approach for subtyping schizophrenia using polygenic SNP allele biomarkers

Subtyping neuropsychiatric disorders like schizophrenia is essential for improving the diagnosis and treatment of complex diseases. Subtyping schizophrenia is challenging because it is polygenic and genetically heterogeneous, rendering the standard symptom-based diagnosis often unreliable and unrepeatable. We developed a novel network-based machine-learning approach, netMoST, to subtyping psychiatric disorders. NetMoST identifies polygenic risk SNP-allele modules from genome-wide genotyping data as polygenic haplotype biomarkers (PHBs) for disease subtyping. We applied netMoST to subtype a cohort of schizophrenia subjects into three distinct biotypes with differentiable genetic, neuroimaging and functional characteristics. The PHBs of the first biotype (36.9% of all patients) were related to neurodevelopment and cognition, the PHBs of the second biotype (28.4%) were enriched for neuroimmune functions, and the PHBs of the third biotype (34.7%) were associated with the transport of calcium ions and neurotransmitters. Neuroimaging patterns provided additional support to the new biotypes, with unique regional homogeneity (ReHo) patterns observed in the brains of each biotype compared with healthy controls. Our findings demonstrated netMoST's capability for uncovering novel biotypes of complex diseases such as schizophrenia. The results also showed the power of exploring polygenic allelic patterns that transcend the conventional GWAS approaches.

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.

The transcription factor MebHLH18 in cassava functions in decreasing low temperature-induced leaf abscission to promote low-temperature tolerance

The reactive oxygen species (ROS) signal regulates stress-induced leaf abscission in cassava. The relationship between the function of the cassava transcription factor bHLH gene and low temperature-induced leaf abscission is still unclear. Here, we report that MebHLH18, a transcription factor, involved in regulating low temperature-induced leaf abscission in cassava. The expression of the MebHLH18 gene was significantly related to low temperature-induced leaf abscission and POD level. Under low temperatures, the levels of ROS scavengers in different cassava genotypes were significantly different in the low temperature-induced leaf abscission process. Cassava gene transformation showed that MebHLH18 overexpression significantly decreased the low temperature-induced leaf abscission rate. Simultaneously, interference expression increased the rate of leaf abscission under the same conditions. ROS analysis showed a connection between the decrease in the low temperature-induced leaf abscission rate caused by MebHLH18 expression and the increase in antioxidant activity. A Genome-wide association studies analysis showed a relationship between the natural variation of the promoter region of MebHLH18 and low temperature-induced leaf abscission. Furthermore, studies showed that the change in MebHLH18 expression was caused by a single nucleotide polymorphism variation in the promoter region upstream of the gene. The high expression of MebHLH18 led to a significant increase in POD activity. The increased POD activity decreased the accumulation of ROS at low temperatures and the rate of leaf abscission. It indicates that the natural variation in the promoter region of MebHLH18 increases antioxidant levels under low temperatures and slows down low temperature-induced leaf abscission.

MeGATAs, functional generalists in interactions between cassava growth and development, and abiotic stresses

The proteins with DNA-binding preference to the consensus DNA sequence (A/T) GATA (A/G) belong to a GATA transcription factor family, with a wide array of biological processes in plants. Cassava (Manihot esculenta) is an important food crop with high production of starch in storage roots. Little was however known about cassava GATA domain-containing genes (MeGATAs). Thirty-six MeGATAs, MeGATA1 to MeGATA36, were found in this study. Some MeGATAs showed a collinear relationship with orthologous genes of Arabidopsis, poplar and potato, rice, maize and sorghum. Eight MeGATA-encoded proteins (MeGATAs) analysed were all localized in the nucleus. Some MeGATAs had potentials of binding ligands and/or enzyme activity. One pair of tandem-duplicated MeGATA17-MeGATA18 and 30 pairs of whole genome-duplicated MeGATAs were found. Fourteen MeGATAs showed low or no expression in the tissues. Nine analysed MeGATAs showed expression responses to abiotic stresses and exogenous phytohormones. Three groups of MeGATA protein interactions were found. Fifty-three miRNAs which can target 18 MeGATAs were identified. Eight MeGATAs were found to target other 292 cassava genes, which were directed to radial pattern formation and phyllome development by gene ontology enrichment, and autophagy by Kyoto Encyclopaedia of Genes and Genomes enrichment. These data suggest that MeGATAs are functional generalists in interactions between cassava growth and development, abiotic stresses and starch metabolism.

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.

Conservation genomics of Agave tequilana Weber var. azul: low genetic differentiation and heterozygote excess in the tequila agave from Jalisco, Mexico

Background

Genetic diversity is fundamental for the survival of species. In particular, in a climate change scenario, it is crucial that populations maintain genetic diversity so they can adapt to novel environmental conditions. Genetic diversity in wild agaves is usually high, with low genetic differentiation among populations, in part maintained by the agave pollinators such as the nectarivorous bats. In cultivated agaves, patterns of genetic diversity vary according to the intensity of use, management, and domestication stage. In Agave tequilana Weber var. azul (A. tequilana thereafter), the plant used for tequila production, clonal propagation has been strongly encouraged. These practices may lead to a reduction in genetic diversity.

Methods

We studied the diversity patterns with genome-wide SNPs, using restriction site associated DNA sequencing in cultivated samples of A. tequilana from three sites of Jalisco, Mexico. For one locality, seeds were collected and germinated in a greenhouse. We compared the genomic diversity, levels of inbreeding, genetic differentiation, and connectivity among studied sites and between adults and juvenile plants.

Results

Agave tequilana presented a genomic diversity of H_T = 0.12. The observed heterozygosity was higher than the expected heterozygosity. Adults were more heterozygous than juveniles. This could be a consequence of heterosis or hybrid vigor. We found a shallow genetic structure (average paired F_ST = 0.0044). In the analysis of recent gene flow, we estimated an average migration rate among the different populations of m = 0.25. In particular, we found a population that was the primary source of gene flow and had greater genomic diversity (H_E and H_O ), so we propose that this population should continue to be monitored as a potential genetic reservoir.

Discussion

Our results may be the consequence of more traditional management in the studied specific region of Jalisco. Also, the exchange of seeds or propagules by producers and the existence of gene flow due to occasional sexual reproduction may play an important role in maintaining diversity in A. tequilana. For populations to resist pests, to continue evolving and reduce their risk of extinction under a climate change scenario, it is necessary to maintain genetic diversity. Under this premise we encourage to continue acting in conservation programs for this species and its pollinators.

Divergence of flowering-related genes to control flowering in five Euphorbiaceae genomes

Reproductive growth and vegetative growth are a pair of main contradictions in the process of plant growth. Flowering, as part of reproductive growth, is a key switch in the life cycle of higher plants, which affects the yield and economic benefits of plants to a certain extent. The Euphorbiaceae species, including castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), cassava (Manihot esculenta), and rubber tree (Hevea brasiliensis), have important economic values because they are raw materials for the production of biodiesel, rubber, etc. The flowering mechanisms are still excluded in the Euphorbiaceae species. The flowering-related genes of Arabidopsis thaliana (Arabidopsis) were used as a reference to determine the orthologs of these genes in Euphorbiaceae genomes. The result showed that 146, 144, 114, 114, and 149 of 207 A. thaliana genes were respectively matched to R. communis, V. fordii, J. curcas, H. brasiliensis, and M. esculenta. These identified genes were clustered into seven pathways including gibberellins, floral meristem identity (FMI), vernalization, photoperiod, floral pathway integrators (FPIs), and autonomous pathways. Then, some key numbers of flowering-related genes are widely conserved in the Euphorbiaceae genomes including but not limited to FPI genes LFY, SOC1, FT, and FMI genes AG, CAL, and FUL. However, some genes, including FRI, FLC, and GO, were missing in several or all five Euphorbiaceae species. In this study, we proposed the putative mechanisms of flowering-related genes to control flowering and provided new candidate flowering genes for using marker-assisted breeding to improve variety quality.

Genome-wide identification and expression of TIFY family in cassava (Manihot esculenta Crantz)

Plant-specific TIFY [TIF(F/Y)XG] proteins serve important roles in the regulation of plant stress responses. This family encodes four subfamilies of proteins, JAZ (JASMONATE ZIM-domain), PPD (PEAPOD), ZML (Zinc-finger Inflorescence-like), and TIFY. In this work, a total of 16 JAZ, 3 PPD, 7 ZML, and 2 TIFY genes were found in cassava (Manihot esculenta Crantz) at the genome-wide level. The phylogenetics, exon-intron structure, motif organization, and conserved domains of these genes were analyzed to characterize the members of the JAZ, PPD, and ZML subfamilies. Chromosome location and synteny analyses revealed that 26 JAZ, PPD, and ZML genes were irregularly distributed across 14 of the 18 chromosomes, and 18 gene pairs were implicated in large-scale interchromosomal segmental duplication events. In addition, JAZ, PPD, and ZML gene synteny comparisons between cassava and three other plant species (Arabidopsis, Populus trichocarpa, and rice) uncovered vital information about their likely evolution. The prediction of protein interaction network and cis-acting elements reveal the function of JAZ, PPD, and ZML genes. Subsequently, expression patterns of JAZ, PPD, and ZML genes were validated by qRT-PCR as being expressed in response to osmotic, salt, and cadmium stress. Moreover, almost all JAZ subfamily genes were responsive to jasmonic acid (JA) treatment. In particular, MeJAZ1, MeJAZ13, and MeJAZ14, were highly up-regulated by three treatments, and these genes may deserve further study. This comprehensive study lays the groundwork for future research into TIFY family genes in cassava and may be valuable for genetic improvement of cassava and other related species.

Regulation of cyanogenic glucosides in wild and domesticated Eusorghum taxa

Domesticated sorghum (Sorghum bicolor [L.] Moench subsp. bicolor) diverts significant amounts of nitrogen away from primary metabolism to the synthesis of cyanogenic glucosides (CNglc) - specialized metabolites that release toxic hydrogen cyanide (HCN). Our aim was to identify the point in the genus Sorghum Moench at which plants gained the ability to maintain hazardous concentrations of cyanogenic glucosides in their leaves into maturity (HCN potential >0.4 mg g^-1 ). This ability occurs in domesticated sorghum (in the subgenus Eusorghum), but not in wild taxa in other Sorghum subgenera. Eight accessions from the subgenus Eusorghum were grown in a common garden: an improved sorghum line, five sorghum landraces, the crop's wild progenitor (S. bicolor subsp. verticilliflorum [Steud.] de Wet ex Wiersema & J. Dahlb.) and wild Sorghum propinquum (Kunth) Hitchc. HCN potential was measured in plants (n = 80) at the three-leaf stage and at 6 weeks old. All study accessions, including the wild taxa, had hazardous CNglc concentrations in the leaves at both the three-leaf stage (mean HCN potentials > = 2.5 mg g^-1 ) and at 6 weeks old (mean HCN potentials > = 0.68 mg g^-1 ), greatly contrasting the much lower mature leaf HCN potentials previously found in wild Sorghum taxa outside subgenus Eusorghum (generally <= 0.01 mg g^-1 ). Our results suggest that the ability to maintain hazardous leaf HCN potentials into maturity might have arisen during the divergence of Eusorghum from other Sorghum subgenera, rather than during the speciation or domestication of S. bicolor, and highlights the value of utilizing Sorghum taxa outside Eusorghum in efforts to improve the crop safety of sorghum.

Pre-Columbian cultivation of vegetatively propagated and fruit tree tropical crops in the Atacama Desert

South America is a megadiverse continent that witnessed the domestication, translocation and cultivation of various plant species from seemingly contrasting ecosystems. It was the recipient and supplier of crops brought to and from Mesoamerica (such as maize and cacao, respectively), and Polynesia to where the key staple crop sweet potato was exported. Not every instance of the trans-ecological expansion of cultivated plants (both domesticated and wild), however, resulted in successful farming. Here, we review the transregional circulation and introduction of five food tropical crops originated in the tropical and humid valleys of the eastern Andes—achira, cassava, ahipa, sweet potato, and pacay—to the hyper-arid coastal valleys of the Atacama Desert of northern Chile, where they have been found in early archeological sites. By means of an evaluation of the contexts of their deposition and supported by direct radiocarbon dating, stable isotopes analyses, and starch grain analysis, we evaluate different hypotheses for explaining their introduction and adaptation to the hyper-arid soils of northern Chile, by societal groups that after the introduction of cultigens still retained a strong dependence on marine hunting, gathering and fishing ways of life based on wide variety of marine coast resources. Many of the studied plants were part of a broader package of introduced goods and technological devices and procedures, linked to food, therapeutic medicine, social and ritual purposes that transformed previous hunter-gatherer social, economic, and ideological institutions. Based on archeological data, we discuss some of the possible socio-ecological processes involved in the development of agricultural landscapes including the adoption of tropical crops originated several hundred kilometers away from the Atacama Desert during the Late Holocene.

Advances in Genetic Analysis and Breeding of Cassava (Manihot esculenta Crantz): A Review

Cassava (Manihot esculenta Crantz) is the sixth most important food crop and consumed by 800 million people worldwide. In Africa, cassava is the second most important food crop after maize and Africa is the worlds' largest producer. Though cassava is not one of the main commodity crops in South Africa, it is becoming a popular crop among farming communities in frost-free areas, due to its climate-resilient nature. This necessitated the establishment of a multi-disciplinary research program at the Agricultural Research Council of South Africa. The objective of this review is to highlight progress made in cassava breeding and genetic analysis. This review highlights the progress of cassava research worldwide and discusses research findings on yield, quality, and adaptability traits in cassava. It also discusses the limitations and the prospects of the cassava R&D program towards development of the cassava industry in South Africa.

DNA fingerprinting reveals varietal composition of Vietnamese cassava germplasm (Manihot esculenta Crantz) from farmers' field and genebank collections

A molecular analysis using informative SNP markers in 1570 clones of cassava from Vietnam reveals varietal composition from farmers' field and genebank collections Cassava is the most important smallholder cash crops in Southeast Asia and is especially used in industrial products. Yet, systematic genetic studies on molecular markers from Vietnamese germplasm have not been considered for breeding and conservation programs. We conducted a molecular analysis of 1570 clones of cassava germplasm from farms across six agro-ecological zones using informative SNP markers. We unraveled the genetic diversity and population structure and provided insights into the value of breeding and conservation programs. Duplicated genotypes comprised 98% of the total sample of the Central Highlands region. Ninety-six SNPs were amplified Central Highlands and South East provinces had the highest allelic richness, covering up to 83% of alleles. The average observed heterozygosity (Ho = 0.43) was slightly higher than expected (He = 0.40) across SNP markers, suggesting an excess of heterozygotes plants. Diversity indexes indicated that cassava populations from North West and Eastern Vietnam are genetically diverse (mean He = 0.40). Genetic parentage tests identified 85 unique genetic groups within the varieties KM94, KM419, BRA1305, KM101, KM140, PER262, KM60, KM57 and two unidentified varieties, which accounted for 82% of the frequency distribution. KM94 is the most dominant variety in Vietnamese farms surveyed (38%), reflecting its superior quality and productivity. Discriminant analysis of principal components (DAPC) revealed four main subgroups, which were partially corroborated by neighbor joining (NJ) analyses. After removing duplicates, 31 unique genotypes were distributed across five of the agro-ecological zones. These were well distributed in the subgroups revealed via DAPC and NJ analyses. The genetic groups identified herein could be used to select unique accessions that should ideally conform with ex situ germplasm collections and identify areas where on-farm conservation programs should be targeted. Newly identified genotypes may also contribute as genetic breeding resources that could be used to adapt cassava to future changes and farmers' needs.

Current status and impending progress for cassava structural genomics

We demystify recent advances in genome assemblies for the heterozygous staple crop cassava (Manihot esculenta), and highlight key cassava genomic resources. Cassava, Manihot esculenta Crantz, is a crop of societal and agricultural importance in tropical regions around the world. Genomics provides a platform for accelerated improvement of cassava's nutritional and agronomic traits, as well as for illuminating aspects of cassava's history including its path towards domestication. The highly heterozygous nature of the cassava genome is widely recognized. However, the full extent and context of this heterozygosity has been difficult to reveal because of technological limitations within genome sequencing. Only recently, with several new long-read sequencing technologies coming online, has the genomics community been able to tackle some similarly difficult genomes. In light of these recent advances, we provide this review to document the current status of the cassava genome and genomic resources and provide a perspective on what to look forward to in the coming years.

Beyond a reference genome: pangenomes and population genomics of underutilized and orphan crops for future food and nutrition security

Underutilized crops are, by definition, under-researched compared to staple crops yet come with traits that may be especially important given climate change and the need to feed a globally increasing population. These crops are often stress-tolerant, and this combined with unique and beneficial nutritional profiles. Whilst progress is being made by generating reference genome sequences, in this Tansley Review, we show how this is only the very first step. We advocate that going 'beyond a reference genome' should be a priority, as it is only at this stage one can identify the specific genes and the adaptive alleles that underpin the valuable traits. We sum up how population genomic and pangenomic approaches have led to the identification of stress- and disease-tolerant alleles in staple crops and compare this to the small number of examples from underutilized crops. We also demonstrate how previously underutilized crops have benefitted from genomic advances and that many breeding targets in underutilized crops are often well studied in staple crops. This cross-crop population-level resequencing could lead to an understanding of the genetic basis of adaptive traits in underutilized crops. This level of investment may be crucial for fully understanding the value of these crops before they are lost.

The haplotype-resolved chromosome pairs of a heterozygous diploid African cassava cultivar reveal novel pan-genome and allele-specific transcriptome features

BACKGROUND

Cassava (Manihot esculenta) is an important clonally propagated food crop in tropical and subtropical regions worldwide. Genetic gain by molecular breeding has been limited, partially because cassava is a highly heterozygous crop with a repetitive and difficult-to-assemble genome.

FINDINGS

Here we demonstrate that Pacific Biosciences high-fidelity (HiFi) sequencing reads, in combination with the assembler hifiasm, produced genome assemblies at near complete haplotype resolution with higher continuity and accuracy compared to conventional long sequencing reads. We present 2 chromosome-scale haploid genomes phased with Hi-C technology for the diploid African cassava variety TME204. With consensus accuracy >QV46, contig N50 >18 Mb, BUSCO completeness of 99%, and 35k phased gene loci, it is the most accurate, continuous, complete, and haplotype-resolved cassava genome assembly so far. Ab initio gene prediction with RNA-seq data and Iso-Seq transcripts identified abundant novel gene loci, with enriched functionality related to chromatin organization, meristem development, and cell responses. During tissue development, differentially expressed transcripts of different haplotype origins were enriched for different functionality. In each tissue, 20-30% of transcripts showed allele-specific expression (ASE) differences. ASE bias was often tissue specific and inconsistent across different tissues. Direction-shifting was observed in <2% of the ASE transcripts. Despite high gene synteny, the HiFi genome assembly revealed extensive chromosome rearrangements and abundant intra-genomic and inter-genomic divergent sequences, with large structural variations mostly related to LTR retrotransposons. We use the reference-quality assemblies to build a cassava pan-genome and demonstrate its importance in representing the genetic diversity of cassava for downstream reference-guided omics analysis and breeding.

CONCLUSIONS

The phased and annotated chromosome pairs allow a systematic view of the heterozygous diploid genome organization in cassava with improved accuracy, completeness, and haplotype resolution. They will be a valuable resource for cassava breeding and research. Our study may also provide insights into developing cost-effective and efficient strategies for resolving complex genomes with high resolution, accuracy, and continuity.

Genomic insights into lineage-specific evolution of the oleosin family in Euphorbiaceae

BACKGROUND

Lipid droplets (LDs) present in land plants serve as an essential energy and carbon reserve for seed germination and seedling development. Oleosins, the most abundant structural proteins of LDs, comprise a small family involved in LD formation, stabilization and degradation. Despite their importance, our knowledge on oleosins is still poor in Euphorbiaceae, a large plant family that contains several important oil-bearing species.

RESULTS

To uncover lineage-specific evolution of oleosin genes in Euphorbiaceae, in this study, we performed a genome-wide identification and comprehensive comparison of the oleosin family in Euphorbiaceae species with available genome sequences, i.e. castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), Mercurialis annua, cassava (Manihot esculenta) and rubber tree (Hevea brasiliensis), and a number of five, five, five, five, eight and eight members were found, respectively. Synteny analysis revealed one-to-one collinear relationship of oleosin genes between the former four (i.e. castor bean, physic nut, tung tree and M. annua) as well as latter two species (i.e. cassava and rubber tree), whereas one-to-one and one-to-two collinear relationships were observed between physic nut and cassava, reflecting the occurrence of one recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber tree. The presence of five ortholog groups representing three previously defined clades (i.e. U, SL and SH) dates back at least to the Malpighiales ancestor, because they are also conserved in poplar (Populus trichocarpa), a tree having experienced one Salicaceae-specific recent WGD. As observed in poplar, WGD was shown to be the main driver for the family expansion in both cassava and rubber tree. Nevertheless, same retention patterns of WGD-derived duplicates observed in cassava and rubber tree are somewhat different from that of poplar, though certain homologous fragments are still present in rubber tree. Further transcriptional profiling revealed an apparent seed-predominant expression pattern of oleosin genes in physic nut, castor bean and rubber tree. Moreover, structure and expression divergence of paralogous pairs were also observed in both cassava and rubber tree.

CONCLUSION

Comparative genomics analysis of oleosin genes reported in this study improved our knowledge on lineage-specific family evolution in Euphorbiaceae, which also provides valuable information for further functional analysis and utilization of key members and their promoters.

Utilizing evolutionary conservation to detect deleterious mutations and improve genomic prediction in cassava

Introduction

Cassava (Manihot esculenta) is an annual root crop which provides the major source of calories for over half a billion people around the world. Since its domestication ~10,000 years ago, cassava has been largely clonally propagated through stem cuttings. Minimal sexual recombination has led to an accumulation of deleterious mutations made evident by heavy inbreeding depression.

Methods

To locate and characterize these deleterious mutations, and to measure selection pressure across the cassava genome, we aligned 52 related Euphorbiaceae and other related species representing millions of years of evolution. With single base-pair resolution of genetic conservation, we used protein structure models, amino acid impact, and evolutionary conservation across the Euphorbiaceae to estimate evolutionary constraint. With known deleterious mutations, we aimed to improve genomic evaluations of plant performance through genomic prediction. We first tested this hypothesis through simulation utilizing multi-kernel GBLUP to predict simulated phenotypes across separate populations of cassava.

Results

Simulations showed a sizable increase of prediction accuracy when incorporating functional variants in the model when the trait was determined by<100 quantitative trait loci (QTL). Utilizing deleterious mutations and functional weights informed through evolutionary conservation, we saw improvements in genomic prediction accuracy that were dependent on trait and prediction.

Conclusion

We showed the potential for using evolutionary information to track functional variation across the genome, in order to improve whole genome trait prediction. We anticipate that continued work to improve genotype accuracy and deleterious mutation assessment will lead to improved genomic assessments of cassava clones.

Phosphorylation of RAV1/2 by KIN10 is essential for transcriptional activation of CAT6/7, which underlies oxidative stress response in cassava

Related to ABI3/VP1 (RAV) transcription factors have important roles in plant stress responses; however, it is unclear whether RAVs regulates oxidative stress response in cassava (Manihot esculenta). In this study, we report that MeRAV1/2 positively regulate oxidative stress resistance and catalase (CAT) activity in cassava. Consistently, RNA sequencing (RNA-seq) identifies three MeCATs that are differentially expressed in MeRAV1/2-silenced cassava leaves. Interestingly, MeCAT6 and MeCAT7 are identified as direct transcriptional targets of MeRAV1/2 via binding to their promoters. In addition, protein kinase MeKIN10 directly interacts with MeRAV1/2 to phosphorylate them at Ser45 and Ser44 residues, respectively, to promote their direct transcriptional activation on MeCAT6 and MeCAT7. Site mutation of MeRAV1S45A or MeRAV2S44A has no significant effect on the activities of MeCAT6 and MeCAT7 promoters or on oxidative stress resistance. In summary, this study demonstrates that the phosphorylation of MeRAV1/2 by MeKIN10 is essential for its direct transcriptional activation of MeCAT6/7 in response to oxidative stress.

A population based expression atlas provides insights into disease resistance and other physiological traits in cassava (Manihot esculenta Crantz)

Cassava, a food security crop in Africa, is grown throughout the tropics and subtropics. Although cassava can provide high productivity in suboptimal conditions, the yield in Africa is substantially lower than in other geographies. The yield gap is attributable to many challenges faced by cassava in Africa, including susceptibility to diseases and poor soil conditions. In this study, we carried out 3’RNA sequencing on 150 accessions from the National Crops Resources Research Institute, Uganda for 5 tissue types, providing population-based transcriptomics resources to the research community in a web-based queryable cassava expression atlas. Differential expression and weighted gene co-expression network analysis were performed to detect 8820 significantly differentially expressed genes (DEGs), revealing similarity in expression patterns between tissue types and the clustering of detected DEGs into 18 gene modules. As a confirmation of data quality, differential expression and pathway analysis targeting cassava mosaic disease (CMD) identified 27 genes observed in the plant–pathogen interaction pathway, several previously identified CMD resistance genes, and two peroxidase family proteins different from the CMD2 gene. Present research work represents a novel resource towards understanding complex traits at expression and molecular levels for the development of resistant and high-yielding cassava varieties, as exemplified with CMD.

Cassava diseases caused by Xanthomonas phaseoli pv. manihotis and Xanthomonas cassavae

Xanthomonas phaseoli pv. manihotis (Xpm) and X. cassavae (Xc) are two bacterial pathogens attacking cassava. Cassava bacterial blight (CBB) is a systemic disease caused by Xpm, which might have dramatic effects on plant growth and crop production. Cassava bacterial necrosis is a nonvascular disease caused by Xc with foliar symptoms similar to CBB, but its impacts on the plant vigour and the crop are limited. In this review, we describe the epidemiology and ecology of the two pathogens, the impacts and management of the diseases, and the main research achievements for each pathosystem. Because Xc data are sparse, our main focus is on Xpm and CBB.

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

Background

Heterozygous genomes are widespread in outcrossing and clonally propagated crops. However, the variation in heterozygosity underlying key agronomic traits and crop domestication remains largely unknown. Cassava is a staple crop in Africa and other tropical regions and has a highly heterozygous genome.

Results

We describe a genomic variation map from 388 resequenced genomes of cassava cultivars and wild accessions. We identify 52 loci for 23 agronomic traits through a genome-wide association study. Eighteen allelic variations in heterozygosity for nine candidate genes are significantly associated with seven key agronomic traits. We detect 81 selective sweeps with decreasing heterozygosity and nucleotide diversity, harboring 548 genes, which are enriched in multiple biological processes including growth, development, hormone metabolisms and responses, and immune-related processes. Artificial selection for decreased heterozygosity has contributed to the domestication of the large starchy storage root of cassava. Selection for homozygous GG allele in MeTIR1 during domestication contributes to increased starch content. Selection of homozygous AA allele in MeAHL17 is associated with increased storage root weight and cassava bacterial blight (CBB) susceptibility. We have verified the positive roles of MeTIR1 in increasing starch content and MeAHL17 in resistance to CBB by transient overexpression and silencing analysis. The allelic combinations in MeTIR1 and MeAHL17 may result in high starch content and resistance to CBB.

Conclusions

This study provides insights into allelic variation in heterozygosity associated with key agronomic traits and cassava domestication. It also offers valuable resources for the improvement of cassava and other highly heterozygous crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02524-7.

Artificial selection optimizes clonality in chaya (Cnidoscolus aconitifolius)

The clonal propagation of crops offers several advantages to growers, such as skipping the juvenile phase, faster growth, and reduced mortality. However, it is not known if the wild ancestors of most clonal crops have a similar ability to reproduce clonally. Therefore, it is unclear whether clonality was an ancestral condition, or if it evolved during domestication in the majority of these crops. Here, I assessed some traits that are relevant to clonal propagation using stem cuttings from chaya (Cnidoscolus aconitifolius) and compared these traits to those of its wild ancestor. Chaya is highly relevant crop to food security in its domestication center (Yucatan Peninsula) and is now cultivated in several countries. Chaya is also an excellent model for assessing the effect of domestication on clonality because wild relatives and selection targets are known. Specifically, I compared resistance to desiccation, water and resource storage, as well as the production of new organs (shoots and leaves) by the stems of wild and domesticated plants. I also compared their performance in root development and clone survival. I found that, relative to their wild ancestors, the stem cuttings of domesticated chaya had 1.1 times greater storage capacity for water and starch. Additionally, the stems of domesticated plants produced 1.25 times more roots, 2.69 times more shoots and 1.94 more leaves, and their clones lived 1.87 times longer than their wild relatives. In conclusion, the results suggest that artificial selection has optimized water and starch storage by stems in chaya. Because these traits also confer greater fitness (i.e. increased fecundity and survival of clones), they can be considered adaptations to clonal propagation in the agroecosystems where this crop is cultivated.

Genome-level diversification of eight ancient tea populations in the Guizhou and Yunnan regions identifies candidate genes for core agronomic traits

The ancient tea plant, as a precious natural resource and source of tea plant genetic diversity, is of great value for studying the evolutionary mechanism, diversification, and domestication of plants. The overall genetic diversity among ancient tea plants and the genetic changes that occurred during natural selection remain poorly understood. Here, we report the genome resequencing of eight different groups consisting of 120 ancient tea plants: six groups from Guizhou Province and two groups from Yunnan Province. Based on the 8,082,370 identified high-quality SNPs, we constructed phylogenetic relationships, assessed population structure, and performed genome-wide association studies (GWAS). Our phylogenetic analysis showed that the 120 ancient tea plants were mainly clustered into three groups and five single branches, which is consistent with the results of principal component analysis (PCA). Ancient tea plants were further divided into seven subpopulations based on genetic structure analysis. Moreover, it was found that the variation in ancient tea plants was not reduced by pressure from the external natural environment or artificial breeding (nonsynonymous/synonymous = 1.05). By integrating GWAS, selection signals, and gene function prediction, four candidate genes were significantly associated with three leaf traits, and two candidate genes were significantly associated with plant type. These candidate genes can be used for further functional characterization and genetic improvement of tea plants.

Autophagy-related genes serve as heat shock protein 90 co-chaperones in disease resistance against cassava bacterial blight

Heat shock protein 90 (HSP90) is involved in plant growth and various stress responses via regulating protein homeostasis. Autophagy keeps cellular homeostasis by recycling the components of cellular cytoplasmic constituents. Although they have similar effects on cellular protein homeostasis, the direct association between HSP90 and autophagy signaling remains unclear in plants, especially in tropical crops. In this study, the correlation between HSP90 and autophagy signaling was systematically analyzed by protein-protein interaction in cassava, one of the most important economy fruit in tropic. In addition, their effects on plant disease response and underlying mechanisms in cassava were investigated by functional genomics and genetic phenotype assay. The potential MeHSP90.9-MeSGT1-MeRAR1 chaperone complex interacts with MeATGs and subsequently triggers autophagy signaling, conferring improved disease resistance to cassava bacterial blight (CBB). On the contrary, HSP90 inhibitor and autophagy inhibitor decreased disease resistance against CBB in cassava, and autophagy may be involved in the potential MeHSP90.9-MeSGT1-MeRAR1 chaperone complex-mediated multiple immune responses. This study highlights the precise modulation of autophagy signaling by potential MeHSP90.9-MeSGT1-MeRAR1 chaperone complex in autophagy-mediated disease resistance to CBB.

Integrating Omics and Gene Editing Tools for Rapid Improvement of Traditional Food Plants for Diversified and Sustainable Food Security

Indigenous communities across the globe, especially in rural areas, consume locally available plants known as Traditional Food Plants (TFPs) for their nutritional and health-related needs. Recent research shows that many TFPs are highly nutritious as they contain health beneficial metabolites, vitamins, mineral elements and other nutrients. Excessive reliance on the mainstream staple crops has its own disadvantages. Traditional food plants are nowadays considered important crops of the future and can act as supplementary foods for the burgeoning global population. They can also act as emergency foods in situations such as COVID-19 and in times of other pandemics. The current situation necessitates locally available alternative nutritious TFPs for sustainable food production. To increase the cultivation or improve the traits in TFPs, it is essential to understand the molecular basis of the genes that regulate some important traits such as nutritional components and resilience to biotic and abiotic stresses. The integrated use of modern omics and gene editing technologies provide great opportunities to better understand the genetic and molecular basis of superior nutrient content, climate-resilient traits and adaptation to local agroclimatic zones. Recently, realizing the importance and benefits of TFPs, scientists have shown interest in the prospection and sequencing of TFPs for their improvements, cultivation and mainstreaming. Integrated omics such as genomics, transcriptomics, proteomics, metabolomics and ionomics are successfully used in plants and have provided a comprehensive understanding of gene-protein-metabolite networks. Combined use of omics and editing tools has led to successful editing of beneficial traits in several TFPs. This suggests that there is ample scope for improvement of TFPs for sustainable food production. In this article, we highlight the importance, scope and progress towards improvement of TFPs for valuable traits by integrated use of omics and gene editing techniques.

A cassava common mosaic virus vector for virus-induced gene silencing in cassava

BACKGROUND

Cassava is an important crop for food security and industry in the least-developed and developing countries. The completion of the cassava genome sequence and identification of large numbers of candidate genes by next-generation sequencing provide extensive resources for cassava molecular breeding and increase the need for rapid and efficient gene function analysis systems in cassava. Several plant virus-induced gene silencing (VIGS) systems have been developed as reverse genetic tools for rapid gene function analysis in cassava. However, these VIGS vectors could cause severe viral symptoms or inefficient gene silencing.

RESULTS

In this study, we constructed agroinfection-compatible infectious cDNA clones of cassava common mosaic virus isolate CM (CsCMV-CM, genus Potexvirus, family Alphaflexiviridae) that causes systemic infection with mild symptoms in cassava. CsCMV-CM was then modified to a viral vector carrying the Nimble cloning frame, which facilitates the rapid and high-throughput cloning of silencing fragments into the viral genome. The CsCMV-based vector successfully silenced phytoene desaturase (PDS) and magnesium chelatase subunit I (ChlI) in different cassava varieties and Nicotiana benthamiana. The silencing of the ChlI gene could persist for more than two months.

CONCLUSIONS

This CsCMV-based VIGS system provides a new tool for rapid and efficient gene function studies in cassava.

Using precision phenotyping to inform de novo domestication

An update on the use of precision phenotyping to assess the potential of lesser cultivated species as candidates for de novo domestication or similar development for future agriculture.

Phospholipases C and D and Their Role in Biotic and Abiotic Stresses

Plants, as sessile organisms, have adapted a fine sensing system to monitor environmental changes, therefore allowing the regulation of their responses. As the interaction between plants and environmental changes begins at the surface, these changes are detected by components in the plasma membrane, where a molecule receptor generates a lipid signaling cascade via enzymes, such as phospholipases (PLs). Phospholipids are the key structural components of plasma membranes and signaling cascades. They exist in a wide range of species and in different proportions, with conversion processes that involve hydrophilic enzymes, such as phospholipase-C (PLC), phospholipase-D (PLD), and phospholipase-A (PLA). Hence, it is suggested that PLC and PLD are highly conserved, compared to their homologous genes, and have formed clusters during their adaptive history. Additionally, they generate responses to different functions in accordance with their protein structure, which should be reflected in specific signal transduction responses to environmental stress conditions, including innate immune responses. This review summarizes the phospholipid systems associated with signaling pathways and the innate immune response.

Transcriptome integrated metabolic modeling of carbon assimilation underlying storage root development in cassava

The existing genome-scale metabolic model of carbon metabolism in cassava storage roots, rMeCBM, has proven particularly resourceful in exploring the metabolic basis for the phenotypic differences between high and low-yield cassava cultivars. However, experimental validation of predicted metabolic fluxes by carbon labeling is quite challenging. Here, we incorporated gene expression data of developing storage roots into the basic flux-balance model to minimize infeasible metabolic fluxes, denoted as rMeCBMx, thereby improving the plausibility of the simulation and predictive power. Three different conceptual algorithms, GIMME, E-Flux, and HPCOF were evaluated. The rMeCBMx-HPCOF model outperformed others in predicting carbon fluxes in the metabolism of storage roots and, in particular, was highly consistent with transcriptome of high-yield cultivars. The flux prediction was improved through the oxidative pentose phosphate pathway in cytosol, as has been reported in various studies on root metabolism, but hardly captured by simple FBA models. Moreover, the presence of fluxes through cytosolic glycolysis and alanine biosynthesis pathways were predicted with high consistency with gene expression levels. This study sheds light on the importance of prediction power in the modeling of complex plant metabolism. Integration of multi-omics data would further help mitigate the ill-posed problem of constraint-based modeling, allowing more realistic simulation.

Variation in production of cyanogenic glucosides during early plant development: A comparison of wild and domesticated sorghum

Domestication has narrowed the genetic diversity found in crop wild relatives, potentially reducing plasticity to cope with a changing climate. The tissues of domesticated sorghum (Sorghum bicolor), especially in younger plants, are cyanogenic and potentially toxic. Species of wild sorghum produce lower levels of the cyanogenic glucoside (CNglc) dhurrin than S. bicolor at maturity, but it is not known if this is also the case during germination and early growth. CNglcs play multiple roles in primary and specialised metabolism in domesticated sorghum and other crop plants. In this study, the temporal and spatial distribution of dhurrin in wild and domesticated sorghum at different growth stages was monitored in leaf, sheath and root tissues up to 35 days post germination using S. bicolor and the wild species S. brachypodum and S. macrospermum as the experimental systems. Growth parameters were also measured and allocation of plant total nitrogen (N%) to both dhurrin and nitrate (NO₃^-) was calculated. Negligible amounts of dhurrin were produced in the leaves of the two wild species compared to S. bicolor. The morphology of the two wild sorghums also differed from S. bicolor, with the greatest differences observed for the more distantly related S. brachypodum. S. bicolor had the highest leaf N% whilst the wild species had significantly higher root N%. Allocation of nitrogen to dhurrin in aboveground tissue was significantly higher in S. bicolor compared to the wild species but did not differ in the roots across the three species. The differences in plant morphology, dhurrin content and re-mobilisation, and nitrate/nitrogen allocation suggest that domestication has affected the functional roles of dhurrin in sorghum.

Microbiome-wide association studies reveal correlations between the structure and metabolism of the rhizosphere microbiome and disease resistance in cassava

Cassava is one of the most important staple food crops in tropical regions. To date, an understanding of the relationship between microbial communities and disease resistance in cassava has remained elusive. In order to explore the relationship among microbiome and phenotypes for further targeted design of microbial community, 16S rRNA and ITS of microbiome of ten cassava varieties were analysed, and a distinctive microbial community in the rhizosphere showed significant interdependence with disease resistance. Shotgun metagenome sequencing was performed to elucidate the structure of microbiomes of cassava rhizosphere. Comprehensive microbiome studies were performed to assess the correlation between the rhizosphere microbiome and disease resistance. Subsequently, the metagenome of rhizosphere microbiome was annotated to obtain taxonomic information at species level and identify metabolic pathways that were significantly associated with cassava disease resistance. Notably, cassava disease resistance was significantly associated with Lactococcus sp., which specifically produces nisin. To definitively explain the role of nisin and underlying mechanism, analysis of nisin biosynthesis-associated genes together with in vitro and in vivo experiments highlighted the effect of nisin on inhibiting the growth of Xanthomonas axonopodis pv. manihotis (Xam) and activating immune response in cassava. The new insights between cassava rhizosphere microbiome especially Lactococcus sp. and disease resistance provide valuable information into further control of cassava disease.

Technological Innovations for Improving Cassava Production in Sub-Saharan Africa

Cassava is crucial for food security of millions of people in sub-Saharan Africa. The crop has great potential to contribute to African development and is increasing its income-earning potential for small-scale farmers and related value chains on the continent. Therefore, it is critical to increase cassava production, as well as its quality attributes. Technological innovations offer great potential to drive this envisioned change. This paper highlights genomic tools and resources available in cassava. The paper also provides a glimpse of how these resources have been used to screen and understand the pattern of cassava genetic diversity on the continent. Here, we reviewed the approaches currently used for phenotyping cassava traits, highlighting the methodologies used to link genotypic and phenotypic information, dissect the genetics architecture of key cassava traits, and identify quantitative trait loci/markers significantly associated with those traits. Additionally, we examined how knowledge acquired is utilized to contribute to crop improvement. We explored major approaches applied in the field of molecular breeding for cassava, their promises, and limitations. We also examined the role of national agricultural research systems as key partners for sustainable cassava production.

Orphan Crops and their Wild Relatives in the Genomic Era

More than half of the calories consumed by humans are provided by three major cereal crops (rice, maize, and wheat). Orphan crops are usually well adapted to low-input agricultural conditions, and they not only play vital roles in local areas but can also contribute to food and nutritional needs worldwide. Interestingly, many wild relatives of orphan crops are important weeds of major crops. Although orphan crops and their wild relatives have received little attentions from researchers for many years, genomic studies have recently been performed on these plants. Here, we provide an overview of genomic studies on orphan crops, with a focus on orphan cereals and their wild relatives. The genomes of at least 12 orphan cereals and/or their wild relatives have been sequenced. In addition to genomic benefits for orphan crop breeding, we discuss the potential ways for mutual utilization of genomic data from major crops, orphan crops, and their wild relatives (including weeds) and provide perspectives on genetic improvement of both orphan and major crops (including de novo domestication of orphan crops) in the coming genomic era.

Analysis of leaf morphology, secondary metabolites and proteins related to the resistance to Tetranychus cinnabarinus in cassava (Manihot esculenta Crantz)

Constitutive resistance of plant can be divided into physical and chemical barriers. Cassava (Manihot esculenta Crantz) is susceptible to mites, especially Tetranychus cinnabarinus. Although significant differences in the resistance to T. cinnabarinus are observed in different cassava cultivars, limited research has been done on the mechanism accounting for the resistance. The aim of this study was to explore the mechanism of resistance to T. cinnabarinus by comparing morphology, secondary metabolites and proteins in different cassava cultivars. The anatomical structure of leaves showed that the cassava cultivar Xinxuan 048 (XX048), which showed a stronger resistance to T. cinnabarinus in both greenhouse testing and three years field evaluation tests (2016–2018), had thicker palisade tissue, spongy tissue, lower epidermis and leaf midrib tissue compared to cultivar Guire 4 (GR4). Greenhouse evaluation demonstrated that originally these cultivars were different, leading to differences in constitutive levels of metabolites. The proteomic analysis of protected leaves in XX048 and GR4 revealed that up-regulated differentially expressed proteins (DEPs) were highly enriched in secondary metabolic pathways, especially in the biosynthesis of flavonoids. This study not only provides a comprehensive data set for overall proteomic changes of leaves in resistant and susceptible cassava, but also sheds light on the morphological characteristics of cassava-mite interaction, secondary metabolite defense responses, and molecular breeding of mite-resistant cassava for effective pest control.

Exploring the chemotypes underlying important agronomic and consumer traits in cassava (Manihot esculenta crantz)

A broad diversity of phenotypes are available within the cassava germplasm collections. The phenotypes include improved nutritional, starch or culinary root quality as well as abiotic and biotic resistance properties. Some of these traits can be found naturally occurring in cassava landraces, whereas others are the result of targeted breeding efforts. For future breeding programmes it is important to know the underlying mechanisms of these desirable traits. Metabolomics can assist in the elucidation of these mechanisms by measuring the end products of the cellular processes conferring the traits of interest. The present study focused on the comparison of two or more variants of the same trait such as high and low culinary quality or resistance and susceptibility to thrips. Overall, eight different traits were assessed. Results showed that amino acids and umami compounds were associated with superior culinary attributes and the phenylpropanoid superpathway plays an important role in pest resistance. Furthermore, the data highlighted a low chemodiversity in African cassavas and that the source-sink relation was still active at the harvest stage.

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.

Gene Regulatory Network Guided Investigations and Engineering of Storage Root Development in Root Crops

The plasticity of plant development relies on its ability to balance growth and stress resistance. To do this, plants have established highly coordinated gene regulatory networks (GRNs) of the transcription factors and signaling components involved in developmental processes and stress responses. In root crops, yields of storage roots are mainly determined by secondary growth driven by the vascular cambium. In relation to this, a dynamic yet intricate GRN should operate in the vascular cambium, in coordination with environmental changes. Despite the significance of root crops as food sources, GRNs wired to mediate secondary growth in the storage root have just begun to emerge, specifically with the study of the radish. Gene expression data available with regard to other important root crops are not detailed enough for us directly to infer underlying molecular mechanisms. Thus, in this review, we provide a general overview of the regulatory programs governing the development and functions of the vascular cambium in model systems, and the role of the vascular cambium on the growth and yield potential of the storage roots in root crops. We then undertake a reanalysis of recent gene expression data generated for major root crops and discuss common GRNs involved in the vascular cambium-driven secondary growth in storage roots using the wealth of information available in Arabidopsis. Finally, we propose future engineering schemes for improving root crop yields by modifying potential key nodes in GRNs.

Production of very-high-amylose cassava by post-transcriptional silencing of branching enzyme genes

High amylose starch can be produced by plants deficient in the function of branching enzymes (BEs). Here we report the production of transgenic cassava (Manihot esculenta Crantz) with starches containing up to 50% amylose due to the constitutive expression of hair-pin dsRNAs targeting the BE1 or BE2 genes. All BE1-RNAi plant lines (BE1i) and BE2-RNAi plant lines (BE2i) were grown up in the field, but with reduced total biomass production. Considerably high amylose content in the storage roots of BE2i plant lines was achieved. Storage starch granules of BE1i and BE2i plants had similar morphology as wild type (WT), however, the size of BE1i starch granules were bigger than that of WT. Comparisons of amylograms and thermograms of all three sources of storage starches revealed dramatic changes to the pasting properties and a higher melting temperature for BE2i starches. Glucan chain length distribution analysis showed a slight increase in chains of DP>36 in BE1i lines and a dramatic increase in glucan chains between DP 10-20 and DP>40 in BE2i lines. Furthermore, BE2i starches displayed a B-type X-ray diffraction pattern instead of the A-type pattern found in BE1i and WT starches. Therefore, cassava BE1 and BE2 function differently in storage root starch biosynthesis.

Exploring dynamic protein-protein interactions in cassava through the integrative interactome network

Protein-protein interactions (PPIs) play an essential role in cellular regulatory processes. Despite, in-depth studies to uncover the mystery of PPI-mediated regulations are still lacking. Here, an integrative interactome network (MePPI-Ux) was obtained by incorporating expression data into the improved genome-scale interactome network of cassava (MePPI-U). The MePPI-U, constructed by both interolog- and domain-based approaches, contained 3,638,916 interactions and 24,590 proteins (59% of proteins in the cassava AM560 genome version 6). After incorporating expression data as information of state, the MePPI-U rewired to represent condition-dependent PPIs (MePPI-Ux), enabling us to envisage dynamic PPIs (DPINs) that occur at specific conditions. The MePPI-Ux was exploited to demonstrate timely PPIs of cassava under various conditions, namely drought stress, brown streak virus (CBSV) infection, and starch biosynthesis in leaf/root tissues. MePPI-Ux_drought and MePPI-Ux_CBSV suggested involved PPIs in response to stress. MePPI-Ux_SB,leaf and MePPI-Ux_SB,root suggested the involvement of interactions among transcription factor proteins in modulating how leaf or root starch is synthesized. These findings deepened our knowledge of the regulatory roles of PPIs in cassava and would undeniably assist targeted breeding efforts to improve starch quality and quantity.

Cassava breeding and agronomy in Asia: 50 years of history and future directions

In Asia, cassava (Manihot esculenta) is cultivated by more than 8 million farmers, driving the rural economy of many countries. The International Center for Tropical Agriculture (CIAT), in partnership with national agricultural research institutes (NARIs), instigated breeding and agronomic research in Asia, 1983. The breeding program has successfully released high-yielding cultivars resulting in an average yield increase from 13.0 t ha-1 in 1996 to 21.3 t ha-1 in 2016, with significant economic benefits. Following the success in increasing yields, cassava breeding has turned its focus to higher-value traits, such as waxy cassava, to reach new market niches. More recently, building resistance to invasive pests and diseases has become a top priority due to the emergent threat of cassava mosaic disease (CMD). The agronomic research involves driving profitability with advanced technologies focusing on better agronomic management practices thereby maintaining sustainable production systems. Remote sensing technologies are being tested for trait discovery and large-scale field evaluation of cassava. In summary, cassava breeding in Asia is driven by a combination of food and market demand with technological innovations to increase the productivity. Further, exploration in the potential of data-driven agriculture is needed to empower researchers and producers for sustainable advancement.