A time series transcriptome analysis of cassava (Manihot esculenta Crantz) varieties challenged with Ugandan cassava brown streak virus

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

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

Genetic analysis of cassava brown streak disease root necrosis using image analysis and genome-wide association studies

Cassava brown streak disease (CBSD) poses a substantial threat to food security. To address this challenge, we used PlantCV to extract CBSD root necrosis image traits from 320 clones, with an aim of identifying genomic regions through genome-wide association studies (GWAS) and candidate genes. Results revealed strong correlations among certain root necrosis image traits, such as necrotic area fraction and necrotic width fraction, as well as between the convex hull area of root necrosis and the percentage of necrosis. Low correlations were observed between CBSD scores obtained from the 1-5 scoring method and all root necrosis traits. Broad-sense heritability estimates of root necrosis image traits ranged from low to moderate, with the highest estimate of 0.42 observed for the percentage of necrosis, while narrow-sense heritability consistently remained low, ranging from 0.03 to 0.22. Leveraging data from 30,750 SNPs obtained through DArT genotyping, eight SNPs on chromosomes 1, 7, and 11 were identified and associated with both the ellipse eccentricity of root necrosis and the percentage of necrosis through GWAS. Candidate gene analysis in the 172.2kb region on the chromosome 1 revealed 24 potential genes with diverse functions, including ubiquitin-protein ligase, DNA-binding transcription factors, and RNA metabolism protein, among others. Despite our initial expectation that image analysis objectivity would yield better heritability estimates and stronger genomic associations than the 1-5 scoring method, the results were unexpectedly lower. Further research is needed to comprehensively understand the genetic basis of these traits and their relevance to cassava breeding and disease management.

Candidate genes for field resistance to cassava brown streak disease revealed through the analysis of multiple data sources

Cassava (Manihot esculenta Crantz) is a food and industrial storage root crop with substantial potential to contribute to managing risk associated with climate change due to its inherent resilience and in providing a biodegradable option in manufacturing. In Africa, cassava production is challenged by two viral diseases, cassava brown streak disease (CBSD) and cassava mosaic disease. Here we detect quantitative trait loci (QTL) associated with CBSD in a biparental mapping population of a Tanzanian landrace, Nachinyaya and AR37-80, phenotyped in two locations over three years. The purpose was to use the information to ultimately facilitate either marker-assisted selection or adjust weightings in genomic selection to increase the efficiency of breeding. Results from this study were considered in relation to those from four other biparental populations, of similar genetic backgrounds, that were phenotyped and genotyped simultaneously. Further, we investigated the co-localization of QTL for CBSD resistance across populations and the genetic relationships of parents based on whole genome sequence information. Two QTL on chromosome 4 for resistance to CBSD foliar symptoms and one on each of chromosomes 11 and 18 for root necrosis were of interest. Of significance within the candidate genes underlying the QTL on chromosome 4 are Phenylalanine ammonia-lyase (PAL) and Cinnamoyl-CoA reductase (CCR) genes and three PEPR1-related kinases associated with the lignin pathway. In addition, a CCR gene was also underlying the root necrosis-resistant QTL on chromosome 11. Upregulation of key genes in the cassava lignification pathway from an earlier transcriptome study, including PAL and CCR, in a CBSD-resistant landrace compared to a susceptible landrace suggests a higher level of basal lignin deposition in the CBSD-resistant landrace. Earlier RNAscope® in situ hybridisation imaging experiments demonstrate that cassava brown streak virus (CBSV) is restricted to phloem vessels in CBSV-resistant varieties, and phloem unloading for replication in mesophyll cells is prevented. The results provide evidence for the involvement of the lignin pathway. In addition, five eukaryotic initiation factor (eIF) genes associated with plant virus resistance were found within the priority QTL regions.

Comparative analysis of salicylic acid levels and gene expression in resistant, tolerant, and susceptible cassava varieties following whitefly-mediated SLCMV infection

Sri Lankan cassava mosaic virus (SLCMV), the primary pathogen responsible for cassava mosaic disease in cassava plantations, is transmitted via infected cutting stems and the whitefly vector, Bemisia tabaci. To obtain better insights into the defense mechanism of cassava against SLCMV, whiteflies were used to induce SLCMV infection for activating the salicylic acid (SA) signaling pathway, which triggers the innate immune system. The study aimed to investigate the specific interactions between viruliferous whiteflies and SA accumulation in resistant (C33), tolerant (Kasetsart 50; KU50), and susceptible (Rayong 11) cassava cultivars by infecting with SLCMV. Leaf samples were collected at various time points, from 1 to 7 days after inoculation (dai). The SA levels were quantified by gas chromatography-mass spectrometry and validated by quantitative reverse transcription polymerase chain reaction. The SA levels increased in KU50 and C33 plants at 2 and 3 dai, respectively, but remained undetected in Rayong11 plants. The expression of PR-9e, PR-7f5, SPS1, SYP121, Hsf8, and HSP90 increased in infected C33 plants at 4 dai, whereas that of KU50 plants decreased immediately at 2 dai, and that of Rayong11 plants increased at 1 dai but gradually decreased thereafter. These findings strongly indicate that SA plays a crucial role in regulating antiviral defense mechanisms, especially in SLCMV-resistant plants. Altogether, the findings provide valuable insights into the mechanisms underlying the activation of SA-mediated anti-SLCMV defense pathways, and the resistance, tolerance, and susceptibility of cassava, which can aid future breeding programs aimed at enhancing SLCMV resistance.

Genetic dissection of cassava brown streak disease in a genomic selection population

Introduction

Cassava brown streak disease (CBSD) is a major threat to food security in East and central Africa. Breeding for resistance against CBSD is the most economical and sustainable way of addressing this challenge.

Methods

This study seeks to assess the (1) performance of CBSD incidence and severity; (2) identify genomic regions associated with CBSD traits and (3) candidate genes in the regions of interest, in the Cycle 2 population of the National Crops Resources Research Institute.

Results

A total of 302 diverse clones were screened, revealing that CBSD incidence across growing seasons was 44%. Severity scores for both foliar and root symptoms ranged from 1.28 to 1.99 and 1.75 to 2.28, respectively across seasons. Broad sense heritability ranged from low to high (0.15 - 0.96), while narrow sense heritability ranged from low to moderate (0.03 - 0.61). Five QTLs, explaining approximately 19% phenotypic variation were identified for CBSD severity at 3 months after planting on chromosomes 1, 13, and 18 in the univariate GWAS analysis. Multivariate GWAS analysis identified 17 QTLs that were consistent with the univariate analysis including additional QTLs on chromosome 6. Seventy-seven genes were identified in these regions with functions such as catalytic activity, ATP-dependent activity, binding, response to stimulus, translation regulator activity, transporter activity among others.

Discussion

These results suggest variation in virulence in the C2 population, largely due to genetics and annotated genes in these QTLs regions may play critical roles in virus initiation and replication, thus increasing susceptibility to CBSD.

Analysis of proteomic changes in cassava cv. Kasetsart 50 caused by Sri Lankan cassava mosaic virus infection

BACKGROUND

Sri Lankan cassava mosaic virus (SLCMV) is a plant virus causing significant economic losses throughout Southeast Asia. While proteomics has the potential to identify molecular markers that could assist the breeding of virus resistant cultivars, the effects of SLCMV infection in cassava have not been previously explored in detail.

RESULTS

Liquid Chromatography-Tandem Mass Spectrometry (LC/MS-MS) was used to identify differentially expressed proteins in SLCMV infected leaves, and qPCR was used to confirm changes at mRNA levels. LC/MS-MS identified 1,813 proteins, including 479 and 408 proteins that were upregulated in SLCMV-infected and healthy cassava plants respectively, while 109 proteins were detected in both samples. Most of the identified proteins were involved in biosynthetic processes (29.8%), cellular processes (20.9%), and metabolism (18.4%). Transport proteins, stress response molecules, and proteins involved in signal transduction, plant defense responses, photosynthesis, and cellular respiration, although present, only represented a relatively small subset of the detected differences. RT-qPCR confirmed the upregulation of WRKY 77 (A0A140H8T1), WRKY 83 (A0A140H8T7), NAC 6 (A0A0M4G3M4), NAC 35 (A0A0M5JAB4), NAC 22 (A0A0M5J8Q6), NAC 54 (A0A0M4FSG8), NAC 70 (A0A0M4FEU9), MYB (A0A2C9VER9 and A0A2C9VME6), bHLH (A0A2C9UNL9 and A0A2C9WBZ1) transcription factors. Additional upregulated transcripts included receptors, such as receptor-like serine/threonine-protein kinase (RSTK) (A0A2C9UPE4), Toll/interleukin-1 receptor (TIR) (A0A2C9V5Q3), leucine rich repeat N-terminal domain (LRRNT_2) (A0A2C9VHG8), and cupin (A0A199UBY6). These molecules participate in innate immunity, plant defense mechanisms, and responses to biotic stress and to phytohormones.

CONCLUSIONS

We detected 1,813 differentially expressed proteins infected cassava plants, of which 479 were selectively upregulated. These could be classified into three main biological functional groups, with roles in gene regulation, plant defense mechanisms, and stress responses. These results will help identify key proteins affected by SLCMV infection in cassava plants.

Current status and impending progress for cassava structural genomics

KEY MESSAGE

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.

Cassava Brown Streak Viruses express second 6-kilodalton (6K2) protein with varied polarity and three dimensional (3D) structures: Basis for trait discrepancy between the virus species

Cassava Brown Streak Virus (CBSV) and Ugandan Cassava Brown Streak Virus (UCBSV) are the two among six virus species speculated to cause the most catastrophic Brown Streak Disease of Cassava (CBSD) in Africa and Asia. Cassava Brown Streak Virus (CBSV) is hard to breed resistance for compared to Ugandan Cassava Brown Streak Virus (UCBSV) species. This is exemplified by incidences of CBSV species rather than UCBSV species in elite breeding line, KBH 2006/0026 at Bagamoyo, Tanzania. It is not yet understood as to why CBSV species could breakdown CBSD-resistance in the KBH 2006/0026 unlike the UCBSV species. This marks the first in silico study conducted to understand molecular basis for the trait discrepancy between CBSV and UCBSV species from structural biology view point. Following ab initio modelling and analysis of physical-chemical properties of second 6-kilodalton (6K2) protein encoded by CBSV and UCBSV species, using ROBETTA server and Protein Parameters tool, respectively we report that; three dimensional (3D) structures and polarity of the protein differs significantly between the two virus species. (95% and 5%) and (85% and 15%) strains of 20 CBSV and 20 UCBSV species respectively, expressed the protein in homo-trimeric and homo-tetrameric forms, correspondingly. 95% and 85% of studied strain population of the two virus species expressed hydrophilic and hydrophobic 6K2, respectively. Based on findings of the curent study, we hypothesize that; (i) The hydrophilic 6K2 expressed by the CBSV species, favour its faster systemic movement via vascular tissues of cassava host and hence result into higher tissue titres than the UCBSV species encoding hydrophobic form of the protein. t and (ii) The hydrophilic 6K2 expressed byCBSV species have additional interaction advantage with Nuclear Inclusion b protease domain (NIb) and Viral genome-linked protein (VPg), components of Virus Replication Complex (VRC) and hence contributing to faster replication of viral genome than the hydrophobic 6K2 expressed by the UCBSV species. Experimental studies are needed to resolve the 3D structures of the 6K2, VPg and NIb and comprehend complex molecular interactions between them. We suggest that, the 6K2 gene should be targeted for improvement of RNA interference (RNAi)-directed transgenesis of virus-resistant cassava as a more effective way to control the CBSD besides breeding.

The Cassava NBS-LRR Genes Confer Resistance to Cassava Bacterial Blight

Cassava bacterial blight (CBB) caused by Xanthomonas axonopodis pv. manihotis (Xam) seriously affects cassava yield. Genes encoding nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains are among the most important disease resistance genes in plants that are specifically involved in the response to diverse pathogens. However, the in vivo roles of NBS-LRR remain unclear in cassava (Manihot esculenta). In this study, we isolated four MeLRR genes and assessed their expression under salicylic acid (SA) treatment and Xam inoculation. Four MeLRR genes positively regulate cassava disease general resistance against Xam via virus-induced gene silencing (VIGS) and transient overexpression. During cassava-Xam interaction, MeLRRs positively regulated endogenous SA and reactive oxygen species (ROS) accumulation and pathogenesis-related gene 1 (PR1) transcripts. Additionally, we revealed that MeLRRs positively regulated disease resistance in Arabidopsis. These pathogenic microorganisms include Pseudomonas syringae pv. tomato, Alternaria brassicicola, and Botrytis cinerea. Our findings shed light on the molecular mechanism underlying the regulation of cassava resistance against Xam inoculation.

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.

WRKY Transcription Factors in Cassava Contribute to Regulation of Tolerance and Susceptibility to Cassava Mosaic Disease through Stress Responses

Among the numerous biological constraints that hinder cassava (Manihot esculenta Crantz) production, foremost is cassava mosaic disease (CMD) caused by virus members of the family Geminiviridae, genus Begomovirus. The mechanisms of CMD tolerance and susceptibility are not fully understood; however, CMD susceptible T200 and tolerant TME3 cassava landraces have been shown to exhibit different large-scale transcriptional reprogramming in response to South African cassava mosaic virus (SACMV). Recent identification of 85 MeWRKY transcription factors in cassava demonstrated high orthology with those in Arabidopsis, however, little is known about their roles in virus responses in this non-model crop. Significant differences in MeWRKY expression and regulatory networks between the T200 and TME3 landraces were demonstrated. Overall, WRKY expression and associated hormone and enriched biological processes in both landraces reflect oxidative and other biotic stress responses to SACMV. Notably, MeWRKY11 and MeWRKY81 were uniquely up and downregulated at 12 and 67 days post infection (dpi) respectively in TME3, implicating a role in tolerance and symptom recovery. AtWRKY28 and AtWRKY40 homologs of MeWRKY81 and MeWRKY11, respectively, have been shown to be involved in regulation of jasmonic and salicylic acid signaling in Arabidopsis. AtWRKY28 is an interactor in the RPW8-NBS resistance (R) protein network and downregulation of its homolog MeWRKY81 at 67 dpi in TME3 suggests a negative role for this WRKY in SACMV tolerance. In contrast, in T200, nine MeWRKYs were differentially expressed from early (12 dpi), middle (32 dpi) to late (67 dpi) infection. MeWRKY27 (homolog AtWRKY33) and MeWRKY55 (homolog AtWRKY53) were uniquely up-regulated at 12, 32 and 67 dpi in T200. AtWRKY33 and AtWRKY53 are positive regulators of leaf senescence and oxidative stress in Arabidopsis, suggesting MeWRKY55 and 27 contribute to susceptibility in T200.

PAL1 gene of the phenylpropanoid pathway increases resistance to the Cassava brown streak virus in cassava

BACKGROUND

The phenylalanine ammonia lyase genes play crucial role in plant response to biotic and abiotic stresses. In this study, we characterized the role of PAL genes in increasing resistance to the Cassava brown streak virus that causes the economically important cassava brown streak disease (CBSD) on cassava in Africa.

METHODS

The whole transcriptomes of eight cassava varieties differing in resistance to CBSD were obtained at 1, 5 and 8 weeks after CBSV infection.

RESULTS

Analysis of RNA-Seq data identified the overexpression of PAL1, PAL2, cinnamic acid and two chalcone synthase genes in CBSD-resistant cassava varieties, which was subsequently confirmed by RT-qPCR. The exogenous application of Acibenzolar-S-Methyl induced PAL1 gene expression to enhance resistance in the susceptible var. Kalawe. In contrast, the silencing of PAL1 by RNA interference led to increased susceptibility of the resistant var. Kaleso to CBSD.

CONCLUSIONS

PAL1 gene of the phenylpropanoid pathway has a major role in inducing resistance to CBSD in cassava plants and its early induction is key for CBSD resistance.

Why has permanent control of cassava brown streak disease in Sub-Saharan Africa remained a dream since the 1930s?

Effective control of ipomoviruses that cause cassava brown streak disease (CBSD) in Africa has remained problematic despite eight remarkable decades (1930-2021) of research efforts. Molecular mechanisms underlying resistance breakdown in genetically improved cassava are still unknown. The vast genetic diversity of cassava brown streak viruses, which is crucial for the improvement of routine reverse transcription polymerase chain reaction (RT-qPCR) assays in CBSD-endemic regions of Africa, is controversial and underrepresented. From a molecular epidemiology viewpoint, this review discusses the reasons for why permanent control of CBSD is difficult in the modern era, even with the presence of diverse in silico and omics tools, recombinant DNA, and high throughput next-generation sequencing technologies. Following an extensive nucleotide data search in the National Centre for Biotechnology Information (NCBI) database and a literature review in PubMed and Scopus, we report that genomic data of 87.62% (474/541) strains of cassava brown streak virus are missing due to poor sequencing capacity in Africa. The evolution dynamics of viral virulence and pathogenicity has not yet been fully explored from the available 67 (12.38%) genomic sequences, owing to poor bioinformatics capacity. Tanzania and Zambia have the highest and lowest disease inoculum pressure, correspondingly. Knowledge gaps in molecular biology and the overall molecular pathogenesis of CBSD viruses impede effective disease control in Africa. Recommendations for possible solutions to the research questions, controversies, and hypotheses raised in this study serve as a roadmap for the invention of more effective CBSD control methods.

Multiple cellular compartments engagement in Nicotiana benthamiana-peanut stunt virus-satRNA interactions revealed by systems biology approach

KEY MESSAGE

PSV infection changed the abundance of host plant's transcripts and proteins associated with various cellular compartments, including ribosomes, chloroplasts, mitochondria, the nucleus and cytosol, affecting photosynthesis, translation, transcription, and splicing. Virus infection is a process resulting in numerous molecular, cellular, and physiological changes, a wide range of which can be analyzed due to development of many high-throughput techniques. Plant RNA viruses are known to replicate in the cytoplasm; however, the roles of chloroplasts and other cellular structures in the viral replication cycle and in plant antiviral defense have been recently emphasized. Therefore, the aim of this study was to analyze the small RNAs, transcripts, proteins, and phosphoproteins affected during peanut stunt virus strain P (PSV-P)-Nicotiana benthamiana interactions with or without satellite RNA (satRNA) in the context of their cellular localization or functional connections with particular cellular compartments to elucidate the compartments most affected during pathogenesis at the early stages of infection. Moreover, the processes associated with particular cell compartments were determined. The 'omic' results were subjected to comparative data analyses. Transcriptomic and small RNA (sRNA)-seq data were obtained to provide new insights into PSV-P-satRNA-plant interactions, whereas previously obtained proteomic and phosphoproteomic data were used to broaden the analysis to terms associated with cellular compartments affected by virus infection. Based on the collected results, infection with PSV-P contributed to changes in the abundance of transcripts and proteins associated with various cellular compartments, including ribosomes, chloroplasts, mitochondria, the nucleus and the cytosol, and the most affected processes were photosynthesis, translation, transcription, and mRNA splicing. Furthermore, sRNA-seq and phosphoproteomic analyses indicated that kinase regulation resulted in decreases in phosphorylation levels. The kinases were associated with the membrane, cytoplasm, and nucleus components.

Topological clustering of regulatory genes confers pathogenic tolerance to cassava brown streak virus (CBSV) in cassava

Robustness, a naïve property of biological systems, enables organisms to maintain functions during perturbation and is crucial for improving the resilience of crops to prevailing stress conditions and diseases, guaranteeing food security. Most studies of robustness in crops have focused on genetic superiority based upon individual genes, overlooking the collaborative actions of multiple responsive genes and the regulatory network topology. This research aims to uncover patterns of gene cooperation leading to organismal robustness by studying the topology of gene co-expression networks (GCNs) of both CBSV virus resistant and susceptible cassava cultivars. The resulting GCNs show higher topological clustering of cooperative genes in the resistant cultivar, suggesting that the network architecture is central to attaining robustness. Despite a reduction in the number of hub genes in the resistant cultivar following the perturbation, essential biological functions contained in the network were maintained through neighboring genes that withstood the shock. The susceptible cultivar seemingly coped by inducing more gene actions in the network but could not maintain the functions required for plant growth. These findings underscore the importance of regulatory network architecture in ensuring phenotypic robustness and deepen our understanding of transcriptional regulation.

Viral inosine triphosphatase: A mysterious enzyme with typical activity, but an atypical function

Plant viruses typically have highly condensed genomes, yet the plant-pathogenic viruses Cassava brown streak virus, Ugandan cassava brown streak virus, and Euphorbia ringspot virus are unusual in encoding an enzyme not yet found in any other virus, the "house-cleaning" enzyme inosine triphosphatase. Inosine triphosphatases (ITPases) are highly conserved enzymes that occur in all kingdoms of life and perform a house-cleaning function by hydrolysing the noncanonical nucleotide inosine triphosphate to inosine monophosphate. The ITPases encoded by cassava brown streak virus and Ugandan cassava brown streak virus have been characterized biochemically and are shown to have typical ITPase activity. However, their biological role in virus infection has yet to be elucidated. Here we review what is known of viral-encoded ITPases and speculate on potential roles in infection with the aim of generating a greater understanding of cassava brown streak viruses, a group of the world's most devastating viruses.

A cassava protoplast system for screening genes associated with the response to South African cassava mosaic virus

Background

The study of transient gene expression in cassava plants during virus infection using existing protocols is laborious and may take approximately fifteen weeks due to cassava’s recalcitrance to transformation. The combination of a protoplast system with CRISPR-mediated gene editing promises to shorten the turnaround time from plant tissue culture to high-throughput gene expression screening for candidate genes. Here, we detail a protocol for screening genes associated with the response to South African cassava mosaic virus (SACMV) in cassava protoplasts, with reference to the ubiquitin E3 ligase gene, MeE3L.

Methods

Cassava protoplasts of model, and SACMV-susceptible and -tolerant genotypes, were transformed with SACMV infectious clones and/or a CRISPR-editing construct targeting the MeE3L using PEG4000-mediated transfection. DNA and RNA were extracted from transformed protoplasts at 24 h post-transfection. Relative SACMV DNA accumulation was determined via qPCR using DpnI-digested total DNA, MeE3L relative expression was determined via reverse transcriptase qPCR, and results were analysed using one-way ANOVA, Tukey’s HSD test and the 2^−ΔΔCTstatistical method. The MeE3L exonic region was sequenced on the ABI 3500XL Genetic Analyzer platform; and sequences were analysed for mutations using MAFTT and MEGA-X software. Construction of a phylogenetic tree was done using the Maximum Likelihood method and Jones-Taylor-Thornton (JTT) matrix-based model.

Results

The differential expression of unedited and mutant MeE3L during SACMV infection of model, susceptible and tolerant cassava protoplasts was determined within 7 weeks after commencement of tissue culture. The study also revealed that SACMV DNA accumulation in cassava protoplasts is genotype-dependent and induces multiple mutations in the tolerant landrace MeE3L homolog. Notably, the susceptible cassava landrace encodes a RINGless MeE3Lwhich is silenced by SACMV-induced mutations. SACMV also induces mutations which silence the MeE3L RING domain in protoplasts from and tolerant cassava landraces.

Conclusions

This protocol presented here halves the turnaround time for high-throughput screening of genes associated with the host response to SACMV. It provides evidence that a cassava E3 ligase is associated with the response to SACMV and forms a basis for validation of these findings by in planta functional and interaction studies.

Genotype by environment cultivar evaluation for cassava brown streak disease resistance in Tanzania

Cassava brown steak disease (CBSD), caused by Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV), is the most important biotic constraint to cassava production in East and Central Africa. Concerted efforts are required to prevent further spread into West Africa as well as to reduce losses in areas already affected. The study reported here was part of a five-country (Kenya, Malawi, Mozambique, Tanzania and Uganda) programme that aimed to identify superior cassava cultivars resistant to CBSD and to disseminate them widely in the region. Seventeen tissue-cultured and virus-tested cultivars were evaluated in Tanzania across nine sites with diverse CBSD inoculum conditions. Experiments were planted using an alpha-lattice design and assessments were made of surrounding inoculum pressure, CBSD foliar and root incidence and root yield at harvest. There were large differences in CBSD infection between sites, with greatest spread recorded from the north-western Lake (Victoria) zone. Differences were driven by Bemisia tabaci whitefly vector abundance and CBSD inoculum pressure. Both CBSV and UCBSV were almost equally represented in cassava fields surrounding experimental plots, although CBSV predominated in the north-west whilst UCBSV was more frequent in coastal and southern sites. However, the incidence of CBSV was much greater than that of UCBSV in initially virus-free experimental plots, suggesting that CBSV is more virulent. Cultivars could be categorised into three groups based on the degree of CBSD symptom expression in shoots and roots. The seven cultivars (F10_30R2, Eyope, Mkumba, Mkuranga1, Narocass1, Nase3 and Orera) in the most resistant category each had shoot and root incidences of less than 20%. Fresh root yield differed between sites and cultivars, but there was no genotype by environment interaction for this trait, probably attributable to the large fertility and soil moisture differences between sites. Susceptible cultivars and the local check performed well in the absence of CBSD pressure, highlighting the importance of exploiting quality and yield traits of local landraces in breeding programmes. Overall, our results emphasized the importance of applying a balanced strategy for CBSD management. This should use both improved and local germplasm resources to generate high yielding cultivars for specific end-user traits, and combine the deployment of improved cultivars with phytosanitary control measures including the use of healthy planting material and planting during periods of reduced CBSD infection.

Enhancing of anthracnose disease resistance indicates a potential role of antimicrobial peptide genes in cassava

Cassava (Manihot esculenta Crantz.) is an important economic crop in tropical countries. Demands for using cassava in food, feed and biofuel industries have been increasing worldwide. Cassava anthracnose disease, caused by Colletotrichum gloeosporioides f.sp. manihotis (CAD), is considered a major problem in cassava production. To minimize the effects of such disease, this study investigated the response of cassava to attack by CAD and how the plants defend themselves against this threat. Genome-wide identification of antimicrobial peptide genes (AMPs) and their expression in response to fungal infection was performed in the resistant cassava cultivar (Huay Bong 60; HB60) in comparison with the highly susceptible cultivar (Hanatee; HN). A total of 114 gene members of AMP were identified in the cassava genome database. Fifty-six gene members were selected for phylogenetic tree construction and analysis of putative cis-acting elements in their promoter regions. Differential expression profiles of six candidate genes were observed in response to CAD infection of both cassava cultivars. Upregulation of snakins, MeSN1 and MeSN2 was found in HB60, whereas MeHEL, Me-AMP-D2 and MeLTP2 were highly induced in HN. The MeLTP1 gene was not expressed in either cultivar. HB60 showed a reduced severity rating in comparison to HN after CAD infection. The biomembrane permeability test of fungal CAD was strongly affected after treatment with protein extract derived from CAD-infected HB60. Altogether, these findings suggest that snakins have a potential function in the CAD defense response in cassava. These results could be useful for cassava improvement programs to fight fungal pathogen.

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.

Genomic and Transcriptomic Analysis Identified Novel Putative Cassava lncRNAs Involved in Cold and Drought Stress

Long non-coding RNAs (lncRNAs) play important roles in the regulation of complex cellular processes, including transcriptional and post-transcriptional regulation of gene expression relevant for development and stress response, among others. Compared to other important crops, there is limited knowledge of cassava lncRNAs and their roles in abiotic stress adaptation. In this study, we performed a genome-wide study of ncRNAs in cassava, integrating genomics- and transcriptomics-based approaches. In total, 56,840 putative ncRNAs were identified, and approximately half the number were verified using expression data or previously known ncRNAs. Among these were 2229 potential novel lncRNA transcripts with unmatched sequences, 250 of which were differentially expressed in cold or drought conditions, relative to controls. We showed that lncRNAs might be involved in post-transcriptional regulation of stress-induced transcription factors (TFs) such as zinc-finger, WRKY, and nuclear factor Y gene families. These findings deepened our knowledge of cassava lncRNAs and shed light on their stress-responsive roles.

Genome-wide analyses of cassava Pathogenesis-related (PR) gene families reveal core transcriptome responses to whitefly infestation, salicylic acid and jasmonic acid

Background

Whiteflies are a threat to cassava (Manihot esculenta), an important staple food in many tropical/subtropical regions. Understanding the molecular mechanisms regulating cassava’s responses against this pest is crucial for developing control strategies. Pathogenesis-related (PR) protein families are an integral part of plant immunity. With the availability of whole genome sequences, the annotation and expression programs of the full complement of PR genes in an organism can now be achieved. An understanding of the responses of the entire complement of PR genes during biotic stress and to the defense hormones, salicylic acid (SA) and jasmonic acid (JA), is lacking. Here, we analyze the responses of cassava PR genes to whiteflies, SA, JA, and other biotic aggressors.

Results

The cassava genome possesses 14 of the 17 plant PR families, with a total of 447 PR genes. A cassava PR gene nomenclature is proposed. Phylogenetic relatedness of cassava PR proteins to each other and to homologs in poplar, rice and Arabidopsis identified cassava-specific PR gene family expansions. The temporal programs of PR gene expression in response to the whitefly (Aleurotrachelus socialis) in four whitefly-susceptible cassava genotypes showed that 167 of the 447 PR genes were regulated after whitefly infestation. While the timing of PR gene expression varied, over 37% of whitefly-regulated PR genes were downregulated in all four genotypes. Notably, whitefly-responsive PR genes were largely coordinately regulated by SA and JA. The analysis of cassava PR gene expression in response to five other biotic stresses revealed a strong positive correlation between whitefly and Xanthomonas axonopodis and Cassava Brown Streak Virus responses and negative correlations between whitefly and Cassava Mosaic Virus responses. Finally, certain associations between PR genes in cassava expansions and response to biotic stresses were observed among PR families.

Conclusions

This study represents the first genome-wide characterization of PR genes in cassava. PR gene responses to six biotic stresses and to SA and JA are demonstrably different to other angiosperms. We propose that our approach could be applied in other species to fully understand PR gene regulation by pathogens, pests and the canonical defense hormones SA and JA.

Haplotype-resolved genomes of geminivirus-resistant and geminivirus-susceptible African cassava cultivars

Background

Cassava is an important food crop in tropical and sub-tropical regions worldwide. In Africa, cassava production is widely affected by cassava mosaic disease (CMD), which is caused by the African cassava mosaic geminivirus that is transmitted by whiteflies. Cassava breeders often use a single locus, CMD2, for introducing CMD resistance into susceptible cultivars. The CMD2 locus has been genetically mapped to a 10-Mbp region, but its organization and genes as well as their functions are unknown.

Results

We report haplotype-resolved de novo assemblies and annotations of the genomes for the African cassava cultivar TME (tropical Manihot esculenta), which is the origin of CMD2, and the CMD-susceptible cultivar 60444. The assemblies provide phased haplotype information for over 80% of the genomes. Haplotype comparison identified novel features previously hidden in collapsed and fragmented cassava genomes, including thousands of allelic variants, inter-haplotype diversity in coding regions, and patterns of diversification through allele-specific expression. Reconstruction of the CMD2 locus revealed a highly complex region with nearly identical gene sets but limited microsynteny between the two cultivars.

Conclusions

The genome maps of the CMD2 locus in both 60444 and TME3, together with the newly annotated genes, will help the identification of the causal genetic basis of CMD2 resistance to geminiviruses. Our de novo cassava genome assemblies will also facilitate genetic mapping approaches to narrow the large CMD2 region to a few candidate genes for better informed strategies to develop robust geminivirus resistance in susceptible cassava cultivars.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0697-6) contains supplementary material, which is available to authorized users.

Using a transcriptome sequencing approach to explore candidate resistance genes against stemphylium blight in the wild lentil species Lens ervoides

BACKGROUND

Stemphylium blight (SB), caused by Stemphylium botryosum, is a devastating disease in lentil production. Although it is known that accessions of Lens ervoides possess superior SB resistance at much higher frequency than the cultivated lentil species, very little is known about the molecular basis regulating SB resistance in L. ervoides. Therefore, a comprehensive molecular study of SB resistance in L. ervoides was needed to exploit this wild resource available at genebanks for use by plant breeders in resistance breeding.

RESULTS

Microscopic and qPCR quantification of fungal growth revealed that 48, 96, and 144 h post-inoculation (hpi) were interesting time points for disease development in L. ervoides recombinant inbred lines (RILs) LR-66-637 (resistant to SB) and LR-66-577 (susceptible to SB). Results of transcriptome sequencing at 0, 48, 96 and 144 hpi showed that 8810 genes were disease-responsive genes after challenge by S. botryosum. Among them, 7526 genes displayed a similar expression trend in both RILs, and some of them were likely involved in non-host resistance. The remaining 1284 genes were differentially expressed genes (DEGs) between RILs. Of those, 712 DEGs upregulated in LR-66-637 were mostly enriched in 'carbohydrate metabolic process', 'cell wall organization or biogenesis', and 'polysaccharide metabolic process'. In contrast, there were another 572 DEGs that were upregulated in LR-66-577, and some of them were enriched in 'oxidation-reduction process', 'asparagine metabolic process' and 'asparagine biosynthetic process'. After comparing DEGs to genes identified in previously described quantitative trait loci (QTLs) for resistance to SB, nine genes were common and three of them showed differential gene expression between a resistant and a susceptible bulk consisting of five RILs each. Results showed that two genes encoding calcium-transporting ATPase and glutamate receptor3.2 were candidate resistance genes, whereas one gene with unknown function was a candidate susceptibility gene.

CONCLUSION

This study provides new insights into the mechanisms of resistance and susceptibility in L. ervoides RILs responding to S. botryosum infection. Furthermore, we identified candidate resistance or susceptibility genes which warrant further gene function analyses, and which could be valuable for resistance breeding, if their role in resistance or susceptibility can be confirmed.

Next generation sequencing and proteomics in plant virology: how is Colombia doing?

Crop production and trade are two of the most economically important activities in Colombia, and viral diseases cause a high negative impact to agricultural sector. Therefore, the detection, diagnosis, control, and management of viral diseases are crucial. Currently, Next-Generation Sequencing (NGS) and ‘Omic’ technologies constitute a right-hand tool for the discovery of novel viruses and for studying virus-plant interactions. This knowledge allows the development of new viral diagnostic methods and the discovery of key components of infectious processes, which could be used to generate plants resistant to viral infections. Globally, crop sciences are advancing in this direction. In this review, advancements in ‘omic’ technologies and their different applications in plant virology in Colombia are discussed. In addition, bioinformatics pipelines and resources for omics data analyses are presented. Due to their decreasing prices, NGS technologies are becoming an affordable and promising means to explore many phytopathologies affecting a wide variety of Colombian crops so as to improve their trade potential.

Differential expression of the Hsf family in cassava under biotic and abiotic stresses

Heat shock transcription factors (Hsfs) are important regulators of biotic and abiotic stress responses in plants. Currently, the Hsf gene family is not well understood in cassava, an important tropical crop. In the present study, 32 MeHsf genes were identified from the cassava genome database, which were divided into three types based on functional domain and motif distribution analyses. Analysis of the differential expression of the genes belonging to the Hsf family in cassava was carried out based on published cassava transcriptome data from tissues/organs (leaf blade, leaf midvein, lateral buds, organized embryogenic structures, friable embryogenic callus, fibrous roots, storage roots, stem, petiole, shoot apical meristem, and root apical meristem) under abiotic stress (cold, drought) or biotic stress (mealybugs. cassava brown streak disease, cassava bacterial blight). The results show the expression diversity of cassava Hsfs genes in various tissues/organs. The transcript levels of MeHsfB3a, MeHsfA6a, MeHsfA2a, and MeHsfA9b were upregulated by abiotic and biotic stresses, such as cold, drought, cassava bacterial blight, cassava brown streak disease, and mealybugs, indicating their potential roles in mediating the response of cassava plants to environment stresses. Further interaction network and co-expression analyses suggests that Hsf genes may interact with Hsp70 family members to resist environmental stresses in cassava. These results provide valuable information for future studies of the functional characterization of the MeHsf gene family.

An optimized isolation protocol yields high-quality RNA from cassava tissues (Manihot esculenta Crantz)

We developed and modified a precise, rapid, and reproducible protocol isolating high-quality RNA from tissues of multiple varieties of cassava plants (Manihot esculenta Crantz). The resulting method is suitable for use in mini, midi, and maxi preparations and rapidly achieves high total RNA yields (170-600 μg·g^-1) using low-cost chemicals and consumables and with minimal contamination from polysaccharides, polyphenols, proteins, and other secondary metabolites. In particular, A ₂₆₀ : A ₂₈₀ ratios were > 2.0 for RNA from various tissues, and all of the present RNA samples yielded ribosomal integrity number values of greater than six. The resulting high purity and quality of isolated RNA will facilitate downstream applications (quantitative reverse transcriptase-polymerase chain reaction or RNA sequencing) in cassava molecular breeding.

Genome-wide association mapping and genomic prediction for CBSD resistance in Manihot esculenta

Cassava (Manihot esculenta Crantz) is an important security crop that faces severe yield loses due to cassava brown streak disease (CBSD). Motivated by the slow progress of conventional breeding, genetic improvement of cassava is undergoing rapid change due to the implementation of quantitative trait loci mapping, Genome-wide association mapping (GWAS), and genomic selection (GS). In this study, two breeding panels were genotyped for SNP markers using genotyping by sequencing and phenotyped for foliar and CBSD root symptoms at five locations in Uganda. Our GWAS study found two regions associated to CBSD, one on chromosome 4 which co-localizes with a Manihot glaziovii introgression segment and one on chromosome 11, which contains a cluster of nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes. We evaluated the potential of GS to improve CBSD resistance by assessing the accuracy of seven prediction models. Predictive accuracy values varied between CBSD foliar severity traits at 3 months after planting (MAP) (0.27-0.32), 6 MAP (0.40-0.42) and root severity (0.31-0.42). For all traits, Random Forest and reproducing kernel Hilbert spaces regression showed the highest predictive accuracies. Our results provide an insight into the genetics of CBSD resistance to guide CBSD marker-assisted breeding and highlight the potential of GS to improve cassava breeding.