Transcriptional response of virus-infected cassava and identification of putative sources of resistance for cassava brown streak disease

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.

Genome-wide identification and expression analysis of Hsp70 family genes in Cassava (Manihot esculenta Crantz)

Hsp70 proteins function as molecular chaperones, regulating various cellular processes in plants. In this study, a genome-wide analysis led to the identification of 22 Hsp70 (MeHsp70) genes in cassava. Phylogenetic relationship studies with other Malpighiales genomes (Populus trichocarpa, Ricinus communis and Salix purpurea) classified MeHsp70 proteins into eight groups (Ia, Ib, Ic, Id, Ie, If, IIa and IIb). Promoter analysis of MeHsp70 genes revealed the presence of tissue-specific, light, biotic and abiotic stress-responsive cis-regulatory elements showing their functional importance in cassava. Meta-analysis of publically available RNA-seq transcriptome datasets showed constitutive, tissue-specific, biotic and abiotic stress-specific expression patterns among MeHsp70s in cassava. Among 22 Hsp70, six MeHsp70s viz., MecHsp70-3, MecHsp70-6, MeBiP-1, MeBiP-2, MeBiP-3 and MecpHsp70-2 displayed constitutive expression, while three MecHsp70s were induced under both drought and cold stress conditions. Five MeHsp70s, MecHsp70-7, MecHsp70-11, MecHsp70-12, MecHsp70-13, and MecHsp70-14 were induced under drought stress conditions. We predicted that 19 MeHsp70 genes are under the regulation of 24 miRNAs. This comprehensive genome-wide analysis of the Hsp70 gene family in cassava provided valuable insights into their functional roles and identified various potential Hsp70 genes associated with stress tolerance and adaptation to environmental stimuli.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03760-3.

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.

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.

Enhanced Convolutional Neural Network Model for Cassava Leaf Disease Identification and Classification

Cassava is a crucial food and nutrition security crop cultivated by small-scale farmers because it can survive in a brutal environment. It is a significant source of carbohydrates in African countries. Sometimes, Cassava crops can be infected by leaf diseases, affecting the overall production and reducing farmers’ income. The existing Cassava disease research encounters several challenges, such as poor detection rate, higher processing time, and poor accuracy. This research provides a comprehensive learning strategy for real-time Cassava leaf disease identification based on enhanced CNN models (ECNN). The existing Standard CNN model utilizes extensive data processing features, increasing the computational overhead. A depth-wise separable convolution layer is utilized to resolve CNN issues in the proposed ECNN model. This feature minimizes the feature count and computational overhead. The proposed ECNN model utilizes a distinct block processing feature to process the imbalanced images. To resolve the color segregation issue, the proposed ECNN model uses a Gamma correction feature. To decrease the variable selection process and increase the computational efficiency, the proposed ECNN model uses global average election polling with batch normalization. An experimental analysis is performed over an online Cassava image dataset containing 6256 images of Cassava leaves with five disease classes. The dataset classes are as follows: class 0: “Cassava Bacterial Blight (CBB)”; class 1: “Cassava Brown Streak Disease (CBSD)”; class 2: “Cassava Green Mottle (CGM)”; class 3: “Cassava Mosaic Disease (CMD)”; and class 4: “Healthy”. Various performance measuring parameters, i.e., precision, recall, measure, and accuracy, are calculated for existing Standard CNN and the proposed ECNN model. The proposed ECNN classifier significantly outperforms and achieves 99.3% accuracy for the balanced dataset. The test findings prove that applying a balanced database of images improves classification performance.

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.

Differential Tropism in Roots and Shoots of Resistant and Susceptible Cassava (Manihot esculenta Crantz) Infected by Cassava Brown Streak Viruses

Cassava brown streak disease (CBSD) is a destructive disease of cassava in Eastern and Central Africa. Because there was no source of resistance in African varieties to provide complete protection against the viruses causing the disease, we searched in South American germplasm and identified cassava lines that did not become infected with the cassava brown streak viruses. These findings motivated further investigations into the mechanism of virus resistance. We used RNAscope^® in situ hybridization to localize cassava brown streak virus in cassava germplasm lines that were highly resistant (DSC 167, immune) or that restricted virus infections to stems and roots only (DSC 260). We show that the resistance in those lines is not a restriction of long-distance movement but due to preventing virus unloading from the phloem into parenchyma cells for replication, thus restricting the virus to the phloem cells only. When DSC 167 and DSC 260 were compared for virus invasion, only a low CBSV signal was found in phloem tissue of DSC 167, indicating that there is no replication in this host, while the presence of intense hybridization signals in the phloem of DSC 260 provided evidence for virus replication in companion cells. In neither of the two lines studied was there evidence of virus replication outside the phloem tissues. Thus, we conclude that in resistant cassava lines, CBSV is confined to the phloem tissues only, in which virus replication can still take place or is arrested.

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.

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.

Computational models to improve surveillance for cassava brown streak disease and minimize yield loss

Cassava brown streak disease (CBSD) is a rapidly spreading viral disease that affects a major food security crop in sub-Saharan Africa. Currently, there are several proposed management interventions to minimize loss in infected fields. Field-scale data comparing the effectiveness of these interventions individually and in combination are limited and expensive to collect. Using a stochastic epidemiological model for the spread and management of CBSD in individual fields, we simulate the effectiveness of a range of management interventions. Specifically we compare the removal of diseased plants by roguing, preferential selection of planting material, deployment of virus-free ‘clean seed’ and pesticide on crop yield and disease status of individual fields with varying levels of whitefly density crops under low and high disease pressure. We examine management interventions for sustainable production of planting material in clean seed systems and how to improve survey protocols to identify the presence of CBSD in a field or quantify the within-field prevalence of CBSD. We also propose guidelines for practical, actionable recommendations for the deployment of management strategies in regions of sub-Saharan Africa under different disease and whitefly pressure.

Cassava is the second largest source of calories in sub-Saharan Africa and is particularly important for the poorest farmers in the region. Cassava brown streak disease is a viral disease that causes cassava tubers to rot, rendering the roots inedible. Recently, the disease has begun to spread towards major cassava growing regions in West Africa from East Africa, where it continues to cause significant yield losses. Improved approaches for disease control are needed to enable small-holder farmers to prepare for and minimize the impact of the disease when their fields become infected. Using a combination of computational methods and mathematical models enables us to screen a much larger range of potential treatments for their likely effectiveness in managing disease and reducing crop loss than would be possible in conventional field trials, which are expensive and logistically difficult to conduct. Our results indicate that regularly planting part of the field with virus-free cassava greatly improves the yield. Removing visibly infected plants and replanting using visibly uninfected plants also improves yield, even when some of these plants may be infected but not yet showing symptoms. We also show how the survey protocol can be optimized to improve estimates of disease severity leading to more effective tailored advice to farmers in regions with different disease pressures.

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.

Assessing the Degeneration of Cassava Under High-Virus Inoculum Conditions in Coastal Tanzania

Cassava brown streak disease (CBSD), caused by cassava brown streak ipomoviruses (CBSIs), has become the most debilitating biotic stress to cassava production in East and Central Africa. Lack of CBSD-resistant varieties has necessitated the search for alternative control measures. Most smallholder farmers reuse stems from previous crops for planting in the new season. Recycling planting material in this way can lead to "degeneration" owing to the compounding effects of disease. In this study, degeneration was defined as the increase in CBSD incidence and reduction in marketable root yield over time. An experiment was established to study the rates of degeneration in selected cassava varieties Chereko, KBH2002_135, Kipusa, Kizimbani, and Mkuranga1 and cultivars Kiroba and Kikombe under high-CBSD inoculum conditions in Bagamoyo, Tanzania from 2013 to 2017. The experiment was replicated across two seasons: the first planted during the long rains (Masika) between March and June and the second planted during the short rains (Vuli) between October and December. Mean abundance of the whitefly vector (Bemisia tabaci) was much greater during the Vuli season (>19 insects per plant) than the Masika season (<2 insects per plant). CBSD shoot symptoms occurred naturally and were observed only on Kikombe, Kiroba, and Kipusa. New materials had overall lower CBSD shoot incidences (1.5%) compared with recycled materials (6.9%) in Masika, although no significant differences were obvious in Vuli. However, Masika (8.7%) had an overall lower CBSD shoot incidence than Vuli (16.5%) in the varieties that had shoot symptoms. CBSD root incidences were higher in Vuli (10.3%) than Masika (4.4%), and root yields in Masika (29.4 t/ha) were significantly greater than those in Vuli (22.5 t/ha). The highest percentage of roots rendered unusable owing to CBSD was observed in Vuli. There was significantly higher unusable root incidence in recycled materials (3.7%) than in new materials (1.4%) in Masika but not in Vuli. Overall root yield was similar between recycled and new materials in either season. Significant reductions in root yield over the course of the experiment were observed both in Masika and Vuli, whereas changes in marketable yield were significant only in Masika. Differences in the response of varieties to degeneration led to the identification of four degeneration patterns, namely "strong," "moderate," "mild," and "delayed" degeneration. The strongest effects of degeneration were most obvious in the susceptible cultivar (Kikombe), which also had the lowest marketable yield in either season. Seasonal differences were a key driver of degeneration, because its effects were much greater in Vuli than Masika. To the best of our knowledge, this work reports the first study of degeneration caused by cassava viruses.[Formula: see text] Copyright © 2019 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

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.

Resistance Against Cassava Brown Streak Viruses From Africa in Cassava Germplasm From South America

Cassava brown streak disease (CBSD) is a severe virus disease of cassava and prevalent in the eastern regions of Africa. The disease is characterized by distinct vein chlorosis and streak symptoms on leaves and stems and necrosis of storage roots. This necrosis can encompass large areas of the root, rendering it inedible so that the entire cassava harvest can be lost. African cassava varieties are susceptible to either of the two viruses causing the disease, cassava brown streak virus (CBSV) and Uganda cassava brown streak virus, and while there are less sensitive varieties, all cassava eventually succumb to the disease. The lack of CBSD resistance in African cassava varieties prompted this search for new sources of virus resistance in the diversity of South American cassava germplasm held in the collection at International Center for Tropical Agriculture, Columbia. Our search for CBSD resistance in South American cassava germplasm accessions revealed that most of the 238 South American cassava lines infected with CBSV established systemic virus infections with moderate to severe disease symptoms on leaves and stems. Fifteen cassava accessions did not become virus infected, remained free of symptoms, and CBSV was undetected by qRT-PCR. When tuberous roots of those lines were examined, necrotic tissue was found in eight lines and CBSV was detected. The remaining seven cassava accessions remained clear of symptoms on all tissues and organs and were virus free. A broad spectrum of virus resistance also including other virus isolates was confirmed for the breeding lines DSC167 and DSC118. While detailed infection experiments with other cassava lines selected for resistance are still ongoing, this indicates that the resistance identified may also hold against a broader diversity of CBSVs. Taken together, we present the results of a comprehensive study on CBSV resistance and susceptibility in cassava germplasm accessions from South America. The virus resistance in cassava germplasm identified provides compelling evidence for the invaluable contribution of germplasm collections to supply the genetic resources for the improvement of our crops.

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.

Simultaneous CRISPR/Cas9-mediated editing of cassava eIF4E isoforms nCBP-1 and nCBP-2 reduces cassava brown streak disease symptom severity and incidence

Cassava brown streak disease (CBSD) is a major constraint on cassava yields in East and Central Africa and threatens production in West Africa. CBSD is caused by two species of positive-sense RNA viruses belonging to the family Potyviridae, genus Ipomovirus: Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV). Diseases caused by the family Potyviridae require the interaction of viral genome-linked protein (VPg) and host eukaryotic translation initiation factor 4E (eIF4E) isoforms. Cassava encodes five eIF4E proteins: eIF4E, eIF(iso)4E-1, eIF(iso)4E-2, novel cap-binding protein-1 (nCBP-1), and nCBP-2. Protein-protein interaction experiments consistently found that VPg proteins associate with cassava nCBPs. CRISPR/Cas9-mediated genome editing was employed to generate ncbp-1, ncbp-2, and ncbp-1/ncbp-2 mutants in cassava cultivar 60444. Challenge with CBSV showed that ncbp-1/ncbp-2 mutants displayed delayed and attenuated CBSD aerial symptoms, as well as reduced severity and incidence of storage root necrosis. Suppressed disease symptoms were correlated with reduced virus titre in storage roots relative to wild-type controls. Our results demonstrate the ability to modify multiple genes simultaneously in cassava to achieve tolerance to CBSD. Future studies will investigate the contribution of remaining eIF4E isoforms on CBSD and translate this knowledge into an optimized strategy for protecting cassava from disease.

Resistance of advanced cassava breeding clones to infection by major viruses in Uganda

Cassava brown streak disease (CBSD) and cassava mosaic disease (CMD) are two viral diseases that cause severe yield losses in cassava of up to 100%, thereby persistently threatening food and income security in sub-Saharan Africa. For effective management of these diseases, there is a critical need to develop and deploy varieties with dual resistance to CBSD and CMD. In this study, we determined the response of advanced breeding lines to field infection by cassava brown streak viruses (CBSVs) and cassava mosaic begomoviruses (CMBs). This aim helped in identifying superior clones for downstream breeding. In total, 220 cassava clones, three in uniform yield trials (UYTs) and 217 in a crossing block trial (CBT), were evaluated for virus and disease resistance. Field data were collected on disease incidence and severity. To detect and quantify CBSVs, 448 and 128 leaf samples from CBSD symptomatic and symptomless plants were analyzed by reverse transcription PCR and real-time quantitative PCR, respectively. In addition, 93 leaf samples from CMD symptomatic plants in the CBT were analyzed by conventional PCR using CMB species-specific primers. In the CBT, 124 (57%) cassava clones did not express CMD symptoms. Of the affected plants, 44 (55%) had single African cassava mosaic virus infection. Single Cassava brown streak virus (CBSV) infections were more prevalent (81.6%) in CBT clones than single Ugandan cassava brown streak virus (UCBSV) infection (3.2%). Of the three advanced clones in the UYT, NAROCASS 1 and NAROCASS 2 had significantly lower (P < 0.05) CBSD severity, incidence, and CBSV load than MH04/0300. In the UYT, only 22% of samples tested had CBSVs, and all showed a negative result for CMBs. The low disease incidence, severity, and viral load associated with NAROCASS 1 and NAROCASS 2 is evidence of their tolerance to both CBSD and CMD. Therefore, these two cassava clones should be utilized in CBSD and CMD management in Uganda, including their utilization as progenitors in further virus resistance breeding.

A method for generating virus-free cassava plants to combat viral disease epidemics in Africa

Here, we report a method to clean cassava plants from viral infections that cause cassava mosaic and brown streak diseases in Africa. Infected plants of resistant or tolerant varieties from Malawi, Mozambique, Kenya, Tanzania and Uganda were cleaned in the UK using a combination of tissue culture, chemotherapy and thermotherapy. In the first cycle of our virus-indexing procedure, we successfully cleaned 27 of the 31 varieties (87%), and after an additional three cleaning cycles, all plants were virus-free. Virus-free tissue-cultured plants were shipped back to Africa for distribution to farmers. This first cross-boundary effort provides important lessons for mitigating the two-major cassava viral diseases.

Immunohistochemical localization of Cassava brown streak virus and its morphological effect on cassava leaves

The localization of Cassava brown streak virus (CBSV) in cassava (Manihot esculenta) leaf tissues was determined and cellular morphological changes in CBSV-infected tissues were evaluated. CBSV-symptomatic leaves were screened with CBSV-specific primers using reverse-transcriptase polymerase chain reaction. Immunohistochemical reactions showed precipitation in CBSV-infected but not CBSV-free tissues, demonstrating successful localization of CBSV. Microscopic inspection showed significantly larger (P < 0.001) midribs in CBSV-infected compared with control (uninfected) leaves. Viral accumulation occurred in middle and lower but rarely in young upper leaves. This immunohistochemical method for virus localization will be invaluable for efficient screening of CBSV and for breeding resistant cassava.

Overexpression of MinE gene affects the plastid division in cassava

The MinE protein plays an important role in plastid division. In this study, the MinE gene was isolated from the cassava (Manihot esculenta Crantz) genome. We isolated high quality and quantity protoplasts and succeed in performing the transient expression of the GFP-fused Manihot esculenta MinE (MeMinE) protein in cassava mesophyll protoplasts. The transient expression of MeMinE-GFP in cassava protoplasts showed that the MeMinE protein was located in the chloroplast. Due to the abnormal division of chloroplasts, overexpression of MeMinE proteins in cassava mesophyll protoplasts could result in fewer and smaller chloroplasts. Overexpression of MeMinE proteins also showed abnormal cell division characteristics and minicell occurrence in Escherichia coli caused by aberrant septation events in the cell poles.

Transcriptome analysis provides insights into the delayed sticky disease symptoms in Carica papaya

Global gene expression analysis indicates host stress responses, mainly those mediated by SA, associated to the tolerance to sticky disease symptoms at pre-flowering stage in Carica papaya. Carica papaya plants develop the papaya sticky disease (PSD) as a result of the combined infection of papaya meleira virus (PMeV) and papaya meleira virus 2 (PMeV2), or PMeV complex. PSD symptoms appear only after C. papaya flowers. To understand the mechanisms involved in this phenomenon, the global gene expression patterns of PMeV complex-infected C. papaya at pre-and post-flowering stages were assessed by RNA-Seq. The result was 633 and 88 differentially expressed genes at pre- and post-flowering stages, respectively. At pre-flowering stage, genes related to stress and transport were up-regulated while metabolism-related genes were down-regulated. It was observed that induction of several salicylic acid (SA)-activated genes, including PR1, PR2, PR5, WRKY transcription factors, ROS and callose genes, suggesting SA signaling involvement in the delayed symptoms. In fact, pre-flowering C. papaya treated with exogenous SA showed a tendency to decrease the PMeV and PMeV2 loads when compared to control plants. However, pre-flowering C. papaya also accumulated transcripts encoding a NPR1-inhibitor (NPR1-I/NIM1-I) candidate, genes coding for UDP-glucosyltransferases (UGTs) and several genes involved with ethylene pathway, known to be negative regulators of SA signaling. At post-flowering, when PSD symptoms appeared, the down-regulation of PR-1 encoding gene and the induction of BSMT1 and JA metabolism-related genes were observed. Hence, SA signaling likely operates at the pre-flowering stage of PMeV complex-infected C. papaya inhibiting the development of PSD symptoms, but the induction of its negative regulators prevents the full-scale and long-lasting tolerance.

Cassava brown streak disease: historical timeline, current knowledge and future prospects

Cassava is the second most important staple food crop in terms of per capita calories consumed in Africa and holds potential for climate change adaptation. Unfortunately, productivity in East and Central Africa is severely constrained by two viral diseases: cassava mosaic disease (CMD) and cassava brown streak disease (CBSD). CBSD was first reported in 1936 from northeast Tanzania. For approximately 70 years, CBSD was restricted to coastal East Africa and so had a relatively low impact on food security compared with CMD. However, at the turn of the 21st century, CBSD re-emerged further inland, in areas around Lake Victoria, and it has since spread through many East and Central African countries, causing high yield losses and jeopardizing the food security of subsistence farmers. This recent re-emergence has attracted intense scientific interest, with studies shedding light on CBSD viral epidemiology, sequence diversity, host interactions and potential sources of resistance within the cassava genome. This review reflects on 80 years of CBSD research history (1936-2016) with a timeline of key events. We provide insights into current CBSD knowledge, management efforts and future prospects for improved understanding needed to underpin effective control and mitigation of impacts on food security.

Molecular insights into Cassava brown streak virus susceptibility and resistance by profiling of the early host response

Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) are responsible for significant cassava yield losses in eastern sub-Saharan Africa. To study the possible mechanisms of plant resistance to CBSVs, we inoculated CBSV-susceptible and CBSV-resistant cassava varieties with a mixed infection of CBSVs using top-cleft grafting. Transcriptome profiling of the two cassava varieties was performed at the earliest time point of full infection (28 days after grafting) in the susceptible scions. The expression of genes encoding proteins in RNA silencing, salicylic acid pathways and callose deposition was altered in the susceptible cassava variety, but transcriptional changes were limited in the resistant variety. In total, the expression of 585 genes was altered in the resistant variety and 1292 in the susceptible variety. Transcriptional changes led to the activation of β-1,3-glucanase enzymatic activity and a reduction in callose deposition in the susceptible cassava variety. Time course analysis also showed that CBSV replication in susceptible cassava induced a strong up-regulation of RDR1, a gene previously reported to be a susceptibility factor in other potyvirus-host pathosystems. The differences in the transcriptional responses to CBSV infection indicated that susceptibility involves the restriction of callose deposition at plasmodesmata. Aniline blue staining of callose deposits also indicated that the resistant variety displays a moderate, but significant, increase in callose deposition at the plasmodesmata. Transcriptome data suggested that resistance does not involve typical antiviral defence responses (i.e. RNA silencing and salicylic acid). A meta-analysis of the current RNA-sequencing (RNA-seq) dataset and selected potyvirus-host and virus-cassava RNA-seq datasets revealed that the conservation of the host response across pathosystems is restricted to genes involved in developmental processes.

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.

Cassava geminivirus agroclones for virus-induced gene silencing in cassava leaves and roots

AIM

We report the construction of a Virus-Induced Gene Silencing (VIGS) vector and an agroinoculation protocol for gene silencing in cassava (Manihot esculenta Crantz) leaves and roots. The African cassava mosaic virus isolate from Nigeria (ACMV-[NOg]), which was initially cloned in a binary vector for agroinoculation assays, was modified for application as VIGS vector. The functionality of the VIGS vector was validated in Nicotiana benthamiana and subsequently applied in wild-type and transgenic cassava plants expressing the uidA gene under the control of the CaMV 35S promoter in order to facilitate the visualization of gene silencing in root tissues. VIGS vectors were targeted to the Mg2+-chelatase gene in wild type plants and both the coding and promoter sequences of the 35S::uidA transgene in transgenic plants to induce silencing. We established an efficient agro-inoculation method with the hyper-virulent Agrobacterium tumefaciens strain AGL1, which allows high virus infection rates. The method can be used as a low-cost and rapid high-throughput evaluation of gene function in cassava leaves, fibrous roots and storage roots.

BACKGROUND

VIGS is a powerful tool to trigger transient sequence-specific gene silencing in planta. Gene silencing in different organs of cassava plants, including leaves, fibrous and storage roots, is useful for the analysis of gene function.

RESULTS

We developed an African cassava mosaic virus-based VIGS vector as well as a rapid and efficient agro-inoculation protocol to inoculate cassava plants. The VIGS vector was validated by targeting endogenous genes from Nicotiana benthamiana and cassava as well as the uidA marker gene in transgenic cassava for visualization of gene silencing in cassava leaves and roots.

CONCLUSIONS

The African cassava mosaic virus-based VIGS vector allows efficient and cost-effective inoculation of cassava for high-throughput analysis of gene function in cassava leaves and roots.

The role of the whitefly, Bemisia tabaci (Gennadius), and farmer practices in the spread of cassava brown streak ipomoviruses

Cassava brown streak disease (CBSD) is arguably the most dangerous current threat to cassava, which is Africa's most important food security crop. CBSD is caused by two RNA viruses: Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV). The roles of the whitefly Bemisia tabaci (Gennadius) and farmer practices in the spread of CBSD were investigated in a set of field and laboratory experiments. The virus was acquired and transmitted by B. tabaci within a short time (5-10 min each for virus acquisition and inoculation), and was retained for up to 48 hr. Highest virus transmission (60%) was achieved using 20-25 suspected viruliferous whiteflies per plant that were given acquisition and inoculation periods of 24 and 48 hr, respectively. Experiments mimicking the agronomic practices of cassava leaf picking or the use of contaminated tools for making cassava stem cuttings did not show the transmission of CBSV or UCBSV. Screenhouse and field experiments in Tanzania showed that the spread of CBSD next to spreader rows was high, and that the rate of spread decreased with increasing distance from the source of inoculum. The disease spread in the field up to a maximum of 17 m in a cropping season. These results collectively confirm that CBSV and UCBSV are transmitted by B. tabaci semipersistently, but for only short distances in the field. This implies that spread over longer distances is due to movements of infected stem cuttings used for planting material. These findings have important implications for developing appropriate management strategies for CBSD.

QTL associated with resistance to cassava brown streak and cassava mosaic diseases in a bi-parental cross of two Tanzanian farmer varieties, Namikonga and Albert

KEY MESSAGE

QTL consistent across seasons were detected for resistance to cassava brown streak disease induced root necrosis and foliar symptoms. The CMD2 locus was detected in an East African landrace, and comprised two QTL. Cassava production in Africa is compromised by cassava brown streak disease (CBSD) and cassava mosaic disease (CMD). To reduce costs and increase the precision of resistance breeding, a QTL study was conducted to identify molecular markers linked to resistance against these diseases. A bi-parental F1 mapping population was developed from a cross between the Tanzanian farmer varieties, Namikonga and Albert. A one-step genetic linkage map comprising 943 SNP markers and 18 linkage groups spanning 1776.2 cM was generated. Phenotypic data from 240 F1 progeny were obtained from two disease hotspots in Tanzania, over two successive seasons, 2013 and 2014. Two consistent QTLs linked to resistance to CBSD-induced root necrosis were identified in Namikonga on chromosomes II (qCBSDRNFc2Nm) and XI (qCBSDRNc11Nm) and a putative QTL on chromosome XVIII (qCBSDRNc18Nm). qCBSDRNFc2Nm was identified at Naliendele in both seasons. The same QTL was also associated with CBSD foliar resistance. qCBSDRNc11Nm was identified at Chambezi in both seasons, and was characterized by three peaks, spanning a distance of 253 kb. Twenty-seven genes were identified within this region including two LRR proteins and a signal recognition particle. In addition, two highly significant CMD resistance QTL (qCMDc12.1A and qCMDc12.2A) were detected in Albert, on chromosome 12. Both qCMDc12.1A and qCMDc12.2A lay within the range of markers reported earlier, defining the CMD2 locus. This is the first time that two loci have been identified within the CMD2 QTL, and in germplasm of apparent East African origin. Additional QTLs with minor effects on CBSD and CMD resistance were also identified.

Community phytosanitation to manage cassava brown streak disease

Cassava viruses are the major biotic constraint to cassava production in Africa. Community-wide action to manage them has not been attempted since a successful cassava mosaic disease control programme in the 1930s/40s in Uganda. A pilot initiative to investigate the effectiveness of community phytosanitation for managing cassava brown streak disease (CBSD) was implemented from 2013 to 2016 in two communities in coastal (Mkuranga) and north-western (Chato) Tanzania. CBSD incidence in local varieties at the outset was >90%, which was typical of severely affected regions of Tanzania. Following sensitization and monitoring by locally-recruited taskforces, there was effective community-wide compliance with the initial requirement to replace local CBSD-infected material with newly-introduced disease-free planting material of improved varieties. The transition was also supported by the free provision of additional seed sources, including maize, sweet potato, beans and cowpeas. Progress of the initiative was followed in randomly-selected monitoring fields in each of the two locations. Community phytosanitation in both target areas produced an area-wide reduction in CBSD incidence, which was sustained over the duration of the programme. In Chato, maximum CBSD incidence was 39.1% in the third season, in comparison with an incidence of >60% after a single season in a control community where disease-free planting material was introduced in the absence of community phytosanitation. Kriging and geospatial analysis demonstrated that inoculum pressure, which was a function of vector abundance and the number of CBSD-infected plants surrounding monitored fields, was a strong determinant of the pattern of CBSD development in monitored fields. In the first year, farmers achieved yield increases with the new varieties relative to the local variety baseline of 94% in Chato (north-west) and 124% in Mkuranga (coast). Yield benefits of the new material were retained up to the final season in each location. The new variety (Mkombozi) introduced under community phytosanitation conditions in Chato yielded 86% more than the same variety from the same source planted in the no-phytosanitation control location. Although there was an 81% reduction in CBSD incidence in the new variety Kiroba introduced under community phytosanitation compared to control conditions in Mkuranga, there was no concomitant yield increase. Variety Kiroba is known to be tolerant to the effects of CBSD, and tuberous roots of infected plants are frequently asymptomatic. Community phytosanitation has the potential to deliver area-wide and sustained reductions in the incidence of CBSD, which also provide significant productivity gains for growers, particularly where introduced varieties do not have high levels of resistant/tolerance to CBSD. The approach should therefore be considered as a potential component for integrated cassava virus management programmes, particularly where new cassava plantations are being established in areas severely affected by CBSD.

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

A time-course transcriptome analysis of two cassava varieties that are either resistant or susceptible to cassava brown streak disease (CBSD) was conducted using RNASeq, after graft inoculation with Ugandan cassava brown streak virus (UCBSV). From approximately 1.92 billion short reads, the largest number of differentially expressed genes (DEGs) was obtained in the resistant (Namikonga) variety at 2 days after grafting (dag) (3887 DEGs) and 5 dag (4911 DEGs). At the same time points, several defense response genes (encoding LRR-containing, NBARC-containing, pathogenesis-related, late embryogenesis abundant, selected transcription factors, chaperones, and heat shock proteins) were highly expressed in Namikonga. Also, defense-related GO terms of 'translational elongation', 'translation factor activity', 'ribosomal subunit' and 'phosphorelay signal transduction', were overrepresented in Namikonga at these time points. More reads corresponding to UCBSV sequences were recovered from the susceptible variety (Albert) (733 and 1660 read counts per million (cpm)) at 45 dag and 54 dag compared to Namikonga (10 and 117 cpm respectively). These findings suggest that Namikonga's resistance involves restriction of multiplication of UCBSV within the host. These findings can be used with other sources of evidence to identify candidate genes and biomarkers that would contribute substantially to knowledge-based resistance breeding.

Variations of five eIF4E genes across cassava accessions exhibiting tolerant and susceptible responses to cassava brown streak disease

Cassava (Manihot esculenta) is an important tropical subsistence crop that is severely affected by cassava brown streak disease (CBSD) in East Africa. The disease is caused by Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV). Both have a (+)-sense single-stranded RNA genome with a 5' covalently-linked viral protein, which functionally resembles the cap structure of mRNA, binds to eukaryotic translation initiation factor 4E (eIF4E) or its analogues, and then enable the translation of viral genomic RNA in host cells. To characterize cassava eIF4Es and their potential role in CBSD tolerance and susceptibility, we cloned five eIF4E transcripts from cassava (accession TMS60444). Sequence analysis indicated that the cassava eIF4E family of proteins consisted of one eIF4E, two eIF(iso)4E, and two divergent copies of novel cap-binding proteins (nCBPs). Our data demonstrated experimentally the coding of these five genes as annotated in the published cassava genome and provided additional evidence for refined annotations. Illumina resequencing data of the five eIF4E genes were analyzed from 14 cassava lines tolerant or susceptible to CBSD. Abundant single nucleotide polymorphisms (SNP) and biallelic variations were observed in the eIF4E genes; however, most of the SNPs were located in the introns and non-coding regions of the exons. Association studies of non-synonymous SNPs revealed no significant association between any SNP of the five eIF4E genes and the tolerance or susceptibility to CBSD. However, two SNPs in two genes were weakly associated with the CBSD responses but had no direct causal-effect relationship. SNPs in an intergenic region upstream of eIF4E_me showed a surprising strong association with CBSD responses. Digital expression profile analysis showed differential expression of different eIF4E genes but no significant difference in gene expression was found between susceptible and tolerant cassava accessions despite the association of the intergenic SNPs with CBSD responses.

Absolute quantification of cassava brown streak virus mRNA by real-time qPCR

Cassava brown streak disease (CBSD) is the most important virus disease of cassava and a major food security threat in Africa. Yearly economic losses of up to $100 million USD have been attributed to CBSD. The lack of information on plant-virus interactions has restricted progress in breeding for CBSD resistance. Virus quantification is becoming a major tool for the quick and reliable assessment of plant host resistance. Therefore, a protocol for specific absolute quantification of Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) was developed. CBSV and UCBSV coat protein (CP) specific standard templates: CBSV (pFer2, 826bp) and UCBSV (pUF1-R1-1, 732) respectively were generated and maintained in a TA cloning vector. These were used to construct standard curves using a TaqMan qPCR assay. Standard curves with acceptable amplification efficiencies (90-105%) and coefficients of determination (R2) greater than 0.99 were obtained. Infected cassava plants were sampled from a screenhouse and the field and used to validate this assay. Results obtained by testing several screenhouse and field samples revealed consistent absolute quantification assays for different CBSV and UCBSV isolates. This study presents the first protocol for absolute quantification of CBSVs and is expected to accelerate screening for CBSD resistance and hence breeding for CBSD resistance. The use of the method presented here should improve the clarity of virus quantification data as the results obtained are not influenced by varietal, host, seasonal or environmental conditions. Screening efficiency will also be greatly improved as there is no need for the use of reference genes consequently allowing for a larger number of samples to be analyzed. This will increase experimental precision in a timely and cost effective manner.

Cassava Mosaic and Brown Streak Diseases: Current Perspectives and Beyond

Cassava is the fourth largest source of calories in the world but is subject to economically important yield losses due to viral diseases, including cassava brown streak disease and cassava mosaic disease. Cassava mosaic disease occurs in sub-Saharan Africa and the Asian subcontinent and is associated with nine begomovirus species, whereas cassava brown streak disease has to date been reported only in sub-Saharan Africa and is caused by two distinct ipomovirus species. We present an overview of key milestones and their significance in the understanding and characterization of these two major diseases as well as their associated viruses and whitefly vector. New biotechnologies offer a wide range of opportunities to reduce virus-associated yield losses in cassava for farmers and can additionally enable the exploitation of this valuable crop for industrial purposes. This review explores established and new technologies for genetic manipulation to achieve desired traits such as virus resistance.

Structure, Expression, and Functional Analysis of the Hexokinase Gene Family in Cassava

Hexokinase (HXK) proteins play important roles in catalyzing hexose phosphorylation and sugar sensing and signaling. To investigate the roles of HXKs in cassava tuber root development, seven HXK genes (MeHXK1-7) were isolated and analyzed. A phylogenetic analysis revealed that the MeHXK family can be divided into five subfamilies of plant HXKs. MeHXKs were clearly divided into type A (MeHXK1) and type B (MeHXK2-7) based on their N-terminal sequences. MeHXK1-5 all had typical conserved regions and similar protein structures to the HXKs of other plants; while MeHXK6-7 lacked some of the conserved regions. An expression analysis of the MeHXK genes in cassava organs or tissues demonstrated that MeHXK2 is the dominant HXK in all the examined tissues (leaves, stems, fruits, tuber phloems, and tuber xylems). Notably, the expression of MeHXK2 and the enzymatic activity of HXK were higher at the initial and expanding tuber stages, and lower at the mature tuber stage. Furthermore, the HXK activity of MeHXK2 was identified by functional complementation of the HXK-deficient yeast strain YSH7.4-3C (hxk1, hxk2, glk1). The gene expression and enzymatic activity of MeHXK2 suggest that it might be the main enzyme for hexose phosphorylation during cassava tuber root development, which is involved in sucrose metabolism to regulate the accumulation of starch.

Eleven years of breeding efforts to combat cassava brown streak disease

Cassava (Manihot esculenta Crantz) production is currently under threat from cassava brown streak disease (CBSD), a disease that is among the seven most serious obstacles to world's food security. Three issues are of significance for CBSD. Firstly, the virus associated with CBSD, has co-evolved with cassava outside its center of origin for at least 90 years. Secondly, that for the last 74 years, CBSD was only limited to the low lands. Thirdly, that most research has largely focused on CBSD epidemiology and virus diversity. Accordingly, this paper focuses on CBSD genetics and/or breeding and hence, presents empirical data generated in the past 11 years of cassava breeding in Uganda. Specifically, this paper provides: 1) empirical data on CBSD resistance screening efforts to identify sources of resistance and/or tolerance; 2) an update on CBSD resistance population development comprising of full-sibs, half-sibs and S1 families and their respective field performances; and 3) insights into chromosomal regions and genes involved in CBSD resistance based on genome wide association analysis. It is expected that this information will provide a foundation for harmonizing on-going CBSD breeding efforts and consequently, inform the future breeding interventions aimed at combating CBSD.

Transcriptome Analysis of Capsicum Chlorosis Virus-Induced Hypersensitive Resistance Response in Bell Capsicum

BACKGROUND

Capsicum chlorosis virus (CaCV) is an emerging pathogen of capsicum, tomato and peanut crops in Australia and South-East Asia. Commercial capsicum cultivars with CaCV resistance are not yet available, but CaCV resistance identified in Capsicum chinense is being introgressed into commercial Bell capsicum. However, our knowledge of the molecular mechanisms leading to the resistance response to CaCV infection is limited. Therefore, transcriptome and expression profiling data provide an important resource to better understand CaCV resistance mechanisms.

METHODOLOGY/PRINCIPAL FINDINGS

We assembled capsicum transcriptomes and analysed gene expression using Illumina HiSeq platform combined with a tag-based digital gene expression system. Total RNA extracted from CaCV/mock inoculated CaCV resistant (R) and susceptible (S) capsicum at the time point when R line showed a strong hypersensitive response to CaCV infection was used in transcriptome assembly. Gene expression profiles of R and S capsicum in CaCV- and buffer-inoculated conditions were compared. None of the genes were differentially expressed (DE) between R and S cultivars when mock-inoculated, while 2484 genes were DE when inoculated with CaCV. Functional classification revealed that the most highly up-regulated DE genes in R capsicum included pathogenesis-related genes, cell death-associated genes, genes associated with hormone-mediated signalling pathways and genes encoding enzymes involved in synthesis of defense-related secondary metabolites. We selected 15 genes to confirm DE expression levels by real-time quantitative PCR.

CONCLUSION/SIGNIFICANCE

DE transcript profiling data provided comprehensive gene expression information to gain an understanding of the underlying CaCV resistance mechanisms. Further, we identified candidate CaCV resistance genes in the CaCV-resistant C. annuum x C. chinense breeding line. This knowledge will be useful in future for fine mapping of the CaCV resistance locus and potential genetic engineering of resistance into CaCV-susceptible crops.

Analysis of Interactions Between Cassava Brown Streak Disease Symptom Types Facilitates the Determination of Varietal Responses and Yield Losses

Cassava brown streak disease (CBSD), caused by cassava brown streak viruses, is recognized as one of the most important plant disease threats to African food security. This study describes the incidence and severity of the different symptom types caused by CBSD, derived from extensive surveys in the country most severely affected by the disease: Tanzania. Total plant incidence and mean root severity of CBSD, recorded from 341 farmers' fields, were both greater in the Coast Zone (49.5% and 3.05), than in the Lake Zone (32.7% and 2.57). Overall, the differing incidences recorded declined in the following order: total plant incidence (39.1%), plant shoot incidence (33.4%), plant root incidence (19.3%), root incidence (10.5%), and unusable root incidence (5.4%). The much lower-than-anticipated loss due to the root necrosis that is characteristic of CBSD was offset by large reductions of root number in plants expressing foliar symptoms of CBSD (15.7% in the Coast Zone and 5.5% in the Lake Zone). These data suggest that the effects of CBSD on the growth of affected plants are greater than those due to root spoilage. Based on these two factors, annual losses due to CBSD in the parts of Tanzania surveyed were estimated at >860,000 t, equivalent to more than U.S.$51 million. A novel approach to using farm-derived data on the responses of the most frequently cultivated varieties to CBSD infection allowed the grouping of the varieties into four categories, based on their relative resistance or tolerance to infection. This tool should be of value to breeders in identifying and selecting for sources of resistance or tolerance in both local and exotic germplasm, and should ultimately contribute to enhancing the management of one of Africa's most damaging crop diseases.

Artificial microRNA-derived resistance to Cassava brown streak disease

Artificial miRNAs (amiRNA) were generated targeting conserved sequences within the genomes of the two causal agents of Cassava brown streak disease (CBSD): Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV). Transient expression studies on ten amiRNAs targeting 21nt conserved sequences of P1(CBSV and UCBSV), P3(CBSV and UCBSV), CI(UCBSV), NIb(CBSV and UCBSV), CP(UCBSV) and the un-translated region (3'-UTR) were tested in Nicotiana benthamiana. Four out of the ten amiRNAs expressed the corresponding amiRNA at high levels. Transgenic N. benthamiana plants were developed for the four amiRNAs targeting the P1 and NIb genes of CBSV and the P1 and CP genes of UCBSV and shown to accumulate miRNA products. Transgenic plants challenged with CBSV and UCBSV isolates showed resistance levels that ranged between ∼20-60% against CBSV and UCBSV and correlated with expression levels of the transgenically derived miRNAs. MicroRNAs targeting P1 and NIb of CBSV showed protection against CBSV and UCBSV, while amiRNAs targeting the P1 and CP of UCBSV showed protection against UCBSV but were less efficient against CBSV. These results indicate a potential application of amiRNAs for engineering resistance to CBSD-causing viruses in cassava.

High-throughput multiplex real-time PCR assay for the simultaneous quantification of DNA and RNA viruses infecting cassava plants

AIMS

To develop a multiplex TaqMan-based real-time PCR assay (qPCR) for the simultaneous detection and quantification of both RNA and DNA viruses affecting cassava (Manihot esculenta) in eastern Africa.

METHODS AND RESULTS

The diagnostic assay was developed for two RNA viruses; Cassava brown streak virus (CBSV) and Uganda cassava brown streak virus (UCBSV) and two predominant DNA viruses; African cassava mosaic virus (ACMV) and East African cassava mosaic virus (EACMV), which cause the economically important cassava brown streak disease (CBSD) and cassava mosaic disease (CMD) respectively. Our method, developed by analysing PCR products of viruses, was highly sensitive to detect target viruses from very low quantities of 4-10 femtograms. Multiplexing did not diminish sensitivity or accuracy compared to uniplex alternatives. The assay reliably detected and quantified four cassava viruses in field samples where CBSV and UCBSV synergy was observed in majority of mixed-infected varieties.

CONCLUSIONS

We have developed a high-throughput qPCR diagnostic assay capable of specific and sensitive quantification of predominant DNA and RNA viruses of cassava in eastern Africa.

SIGNIFICANCE AND IMPACT OF THE STUDY

The qPCR methods are a great improvement on the existing methods and can be used for monitoring virus spread as well as for accurate evaluation of the cassava varieties for virus resistance.

Comparative analysis of virus-derived small RNAs within cassava (Manihot esculenta Crantz) infected with cassava brown streak viruses

•

The 21-nt virus-derived small RNAs were predominant, followed by the 22-nt class.

•

Susceptible cassava genotypes accumulated higher CBSV- than UCBSV-derived small RNA.

•

Tolerant cassava genotype accumulated high CBSV- and low UCBSV-derived small RNAs.

•

AGO2, DCL2 and DCL4 were differentially regulated in CBSV/UCBSV-infected plants.

•

CBSV and UCBSV interact differently in the same host genetic background

Infection of plant cells by viral pathogens triggers RNA silencing, an innate antiviral defense mechanism. In response to infection, small RNAs (sRNAs) are produced that associate with Argonaute (AGO)-containing silencing complexes which act to inactivate viral genomes by posttranscriptional gene silencing (PTGS). Deep sequencing was used to compare virus-derived small RNAs (vsRNAs) in cassava genotypes NASE 3, TME 204 and 60444 infected with the positive sense single-stranded RNA (+ssRNA) viruses cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV), the causal agents of cassava brown streak disease (CBSD). An abundance of 21–24 nt vsRNAs was detected and mapped, covering the entire CBSV and UCBSV genomes. The 21 nt vsRNAs were most predominant, followed by the 22 nt class with a slight bias toward sense compared to antisense polarity, and a bias for adenine and uracil bases present at the 5′-terminus. Distribution and frequency of vsRNAs differed between cassava genotypes and viral genomes. In susceptible genotypes TME 204 and 60444, CBSV-derived sRNAs were seen in greater abundance than UCBSV-derived sRNAs. NASE 3, known to be resistant to UCBSV, accumulated negligible UCBSV-derived sRNAs but high populations of CBSV-derived sRNAs. Transcript levels of cassava homologues of AGO2, DCL2 and DCL4, which are central to the gene-silencing complex, were found to be differentially regulated in CBSV- and UCBSV-infected plants across genotypes, suggesting these proteins play a role in antiviral defense. Irrespective of genotype or viral pathogen, maximum populations of vsRNAs mapped to the cytoplasmic inclusion, P1 and P3 protein-encoding regions. Our results indicate disparity between CBSV and UCBSV host-virus interaction mechanisms, and provide insight into the role of virus-induced gene silencing as a mechanism of resistance to CBSD.

Risk assessment of gene flow from genetically engineered virus resistant cassava to wild relatives in Africa: an expert panel report

The probability and consequences of gene flow to wild relatives is typically considered in the environmental risk assessment of genetically engineered crops. This is a report from a discussion by a group of experts who used a problem formulation approach to consider existing information for risk assessment of gene flow from cassava (Manihot esculenta) genetically engineered for virus resistance to the 'wild' (naturalized) relative M. glaziovii in East Africa. Two environmental harms were considered in this case: (1) loss of genetic diversity in the germplasm pool, and (2) loss of valued species, ecosystem resources, or crop yield and quality due to weediness or invasiveness of wild relatives. Based on existing information, it was concluded that gene flow will occur, but it is not likely that this will reduce the genetic diversity in the germplasm pool. There is little existing information about the impact of the virus in natural populations that could be used to inform a prediction about whether virus resistance would lead to an increase in reproduction or survival, hence abundance of M. glaziovii. However, an increase in the abundance of M. glaziovii should be manageable, and would not necessarily lead to the identified environmental harms.

Resistance gene analogs involved in tolerant cassava--geminivirus interaction that shows a recovery phenotype

The current literature describes recovery from virus-induced symptoms as a RNA silencing defense, but immunity-related genes, including the structurally specific resistance gene analogs (RGAs) that may play a key role in tolerance and recovery is not yet reported. In this study, the transcriptome data of tolerant cassava TME3 (which exhibits a recovery phenotype) and susceptible cassava T200 infected with South African cassava mosaic virus were explored for RGAs. Putative resistance protein analogs (RPAs) with amide-like indole-3-acetic acid-Ile-Leu-Arg (IAA-ILR) and leucine-rich repeat (LRR)-kinase conserved domains were unique to TME3. Common responsive RPAs in TME3 and T200 were the dirigent-like protein, coil-coil nucleotide-binding site (NBS) and toll-interleukin-resistance, disease resistance zinc finger chromosome condensation-like protein (DZC), and NBS-apoptosis repressor with caspase recruitment (ARC)-LRR domains. Mutations in RPAs in the MHD motif of the NBS-ARC2 subdomain associated with the recovery phase in TME3 were observed. Additionally, a cohort of 25 RGAs mined solely during the recovery process in TME3 was identified. Phylogenetic and expression analyses support that diverse RGAs are differentially expressed during tolerance and recovery. This study reveals that in cassava, a perennial crop, RGAs participate in tolerance and differentially accumulate during recovery as a complementary defense mechanism to natural occurring RNA silencing to impair viral replication.

Cassava brown streak disease: a threat to food security in Africa

Cassava brown streak disease (CBSD) has emerged as the most important viral disease of cassava (Manihot esculenta) in Africa and is a major threat to food security. CBSD is caused by two distinct species of ipomoviruses, Cassava brown streak virus and Ugandan cassava brown streak virus, belonging to the family Potyviridae. Previously, CBSD was reported only from the coastal lowlands of East Africa, but recently it has begun to spread as an epidemic throughout the Great Lakes region of East and Central Africa. This new spread represents a major threat to the cassava-growing regions of West Africa. CBSD-resistant cassava cultivars are being developed through breeding, and transgenic RNA interference-derived field resistance to CBSD has also been demonstrated. This review aims to provide a summary of the most important studies on the aetiology, epidemiology and control of CBSD and to highlight key research areas that need prioritization.

Conventional breeding, marker-assisted selection, genomic selection and inbreeding in clonally propagated crops: a case study for cassava

Consolidates relevant molecular and phenotypic information on cassava to demonstrate relevance of heterosis, and alternatives to exploit it by integrating different tools. Ideas are useful to other asexually reproduced crops. Asexually propagated crops offer the advantage that all genetic effects can be exploited in farmers' production fields. However, non-additive effects complicate selection because, while influencing the performance of the materials under evaluation, they cannot be transmitted efficiently to the following cycle of selection. Cassava can be used as a model crop for asexually propagated crops because of its diploid nature and the absence of (known) incompatibility effects. New technologies such as genomic selection (GS), use of inbred progenitors based on doubled haploids and induction of flowering can be employed for accelerating genetic gains in cassava. Available information suggests that heterosis, non-additive genetic effects and within-family variation are relatively large for complex traits such as fresh root yield, moderate for dry matter or starch content in the roots, and low for defensive traits (pest and disease resistance) and plant architecture. The present article considers the potential impact of different technologies for maximizing gains for key traits in cassava, and highlights the advantages of integrating them. Exploiting heterosis would be optimized through the implementation of reciprocal recurrent selection. The advantages of using inbred progenitors would allow shifting the current cassava phenotypic recurrent selection method into line improvement, which in turn would allow designing outstanding hybrids rather than finding them by trial and error.

Identification and distribution of the NBS-LRR gene family in the Cassava genome

BACKGROUND

Plant resistance genes (R genes) exist in large families and usually contain both a nucleotide-binding site domain and a leucine-rich repeat domain, denoted NBS-LRR. The genome sequence of cassava (Manihot esculenta) is a valuable resource for analysing the genomic organization of resistance genes in this crop.

RESULTS

With searches for Pfam domains and manual curation of the cassava gene annotations, we identified 228 NBS-LRR type genes and 99 partial NBS genes. These represent almost 1% of the total predicted genes and show high sequence similarity to proteins from other plant species. Furthermore, 34 contained an N-terminal toll/interleukin (TIR)-like domain, and 128 contained an N-terminal coiled-coil (CC) domain. 63% of the 327 R genes occurred in 39 clusters on the chromosomes. These clusters are mostly homogeneous, containing NBS-LRRs derived from a recent common ancestor.

CONCLUSIONS

This study provides insight into the evolution of NBS-LRR genes in the cassava genome; the phylogenetic and mapping information may aid efforts to further characterize the function of these predicted R genes.