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

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

Cassava Brown Streak Disease Response and Association With Agronomic Traits in Elite Nigerian Cassava Cultivars

Cassava mosaic geminiviruses (CMGs) and cassava brown streak viruses (CBSVs) cause the highest yield losses in cassava production in Africa. In particular, cassava brown streak disease (CBSD) is and continues to be a significant constraint to optimal cassava production in Eastern and Southern Africa. While CBSD has not been reported in West Africa, its recent rapid spread and damage to cassava productivity in Eastern, and Southern Africa is alarming. The aim of this study was to evaluate Nigerian cassava genotypes in order to determine their responses to CBSD, in the event that it invades Nigeria, the world’s largest cassava producer. The study gathered information on whether useful CBSD resistance alleles are present in the elite Nigerian cassava accessions. A total of 1,980 full-sib cassava seedlings from 106 families were assessed in the field at the seedling stage for a year. A subset of 569 clones were selected and assessed for another year at the clonal stage in Namulonge, central Uganda, a known hotspot for CBSD screening. Results indicated that foliar and root incidences and severities varied significantly (p ≤ 0.01, p ≤ 0.001) except for CBSD foliar incidence at 6 months (CBSD6i). Highest and lowest plot-based heritability estimates for CBSD were registered for CBSD root severity (CBSDrs) (0.71) and CBSD6i (0.5). Positive and highly significant correlations were noted between CBSD root incidence (CBSDri) and CBSDrs (r = 0.90***). Significant positive correlations were also noted between CBSD foliar severity at 3 months (CBSD3s) and CBSD foliar incidence at 6 months (CBSD6i) (r = 0.77***), CBSD3s and CBSDrs (r = 0.35***). Fresh root weight (FreshRW) negatively correlated with CBSDri and CBSDrs, respectively (r = −0.21*** and r = −0.22***). Similarly, CBSD3s correlated negatively with cassava mosaic disease severity at 3 (CMD3s) and 6 months (CMD6s), respectively (r = −0.25*** and r = −0.21***). Fifteen clones were selected using a non-weighted summation selection index for further screening. In conclusion, results revealed that the elite Nigerian accessions exhibited significant susceptibility to CBSD within 2 years of evaluation period. It is expected that this information will aid future breeding decisions for the improvement of CBSD resistance among the Nigerian cassava varieties.

Evaluation and Characterization of High Yielding Cassava Mosaic Resistant Variety YTP2

The cassava variety YTP2 (Me 681) has been developed through selection from Thondamuthur type at Tapioca and Castor Research Station, TNAU, Yethapur. The performance of YTP2 in the Adaptive Research Trial (ART) and On Farm Trial (OFT) in the farmer’s field inferred that this new variety is well adapted to cassava growing districts of Tamil Nadu. In addition to the above, YTP2 was found to be resistant to cassava mosaic disease incidence (CMD). Plants are erect, medium growing and non-branching type and suitable for growing under irrigated and rainfed conditions. The internodal length is shorter and the leaf size is medium with sufficient canopy. The leaves of the plants droop down to reduce the transpiration loss which is more advantageous to overcome or escape from drought and heat stress during summer season. It is a dual purpose variety wherein the tubers contain high starch content which is much favourable for the manufacture of starch, sago and also suited for table purpose. The overall performance of this variety showed higher tuber yield (42.20 t ha-1) and starch content (28.40%) which is 15.94% and 18.20% increase over the check varieties YTP1 and H226 respectively. The results of DNA fingerprint data involving SSR markers (SSRY235, NS169 and NS928) showed that it is genetically distinct from the existing commercial varieties viz., YTP1, H226 and Sree Athulya.

Proteome Mapping of South African Cassava Mosaic Virus-Infected Susceptible and Tolerant Landraces of Cassava

The production of cassava is threatened by the geminivirus South African cassava mosaic virus (SACMV), which causes cassava mosaic disease. Cassava landrace TME3 shows tolerance to SACMV, while T200 is highly susceptible. This study aimed to identify the leaf proteome involved in anti-viral defence. Liquid chromatography mass spectrometry (LC-MS) identified 2682 (54 differentially expressed) and 2817 (206 differentially expressed) proteins in both landraces at systemic infection (32 days post infection) and symptom recovery (67 days post infection), respectively. Differences in the number of differentially expressed proteins (DEPs) between the two landraces were observed. Gene ontology analysis showed that defence-associated pathways such as the chloroplast, proteasome, and ribosome were overrepresented at 67 days post infection (dpi) in SACMV-tolerant TME3. At 67 dpi, a high percentage (56%) of over-expressed proteins were localized in the chloroplast in TME3 compared to T200 (31% under-expressed), proposing that chloroplast proteins play a role in tolerance in TME3. Ribosomal_L7Ae domain-containing protein (Manes.12G139100) was over-expressed uniquely in TME3 at 67 dpi and interacts with the ribosomal protein Sac52 (RPL10). RPL10 is a known key player in the NIK1-mediated effector triggered immunity (ETI) response to geminivirus infection, indicating a possible role for Sac52 in SACMV recovery in TME3. In conclusion, differential protein expression responses in TME3 and T200 may be key to unravel tolerance to CMD.

Survey and molecular detection of Sri Lankan cassava mosaic virus in Thailand

Cassava plantations in an area of 458 hectares spanning five provinces along the Thailand–Cambodia border were surveyed from October 2018 to July 2019 to determine the prevalence of cassava mosaic disease (CMD) caused by Sri Lankan cassava mosaic virus (SLCMV) in the region. CMD prevalence was 40% in the whole area and 80% in Prachinburi, 43% in Sakaeo, 37% in Burium, 25% in Surin, and 19% in Sisaket provinces. Disease incidence of CMD was highest 43.08% in Sakaeo, followed by 26.78% in Prachinburi, 7% in Burium, 2.58% in Surin, and 1.25% in Sisaket provinces. Disease severity of CMD symptoms was mild chlorosis to moderate mosaic (2–3). The greatest disease severity was recorded in Prachinburi and Sakaeo provinces. Asymptomatic plants were identified in Surin (12%), Prachinburi (5%), Sakaeo (0.2%), and Buriram (0.1%) by PCR analysis. Cassava cultivars CMR-89 and Huai Bong 80 were susceptible to CMD. In 95% of cases, the infection was transmitted by whiteflies (Bemisia tabaci), which were abundant in Sakaeo, Buriram, and Prachinburi but were sparse in Surin; their densities were highest in May and June 2019. Nucleotide sequencing of the mitochondrial cytochrome oxidase 1 (mtCO1) gene of whiteflies in Thailand revealed that it was similar to the mtCO1 gene of Asia II 1 whitefly. Furthermore, the AV1 gene of SLCMV—which encodes the capsid protein—showed 90% nucleotide identity with SLCMV. Phylogenetic analysis of completed nucleotide sequences of DNA-A and DNA-B components of the SLCMV genome determined by rolling circle amplification (RCA) indicated that they were similar to the nucleotide sequence of SLCMV isolates from Thailand, Vietnam, and Cambodia. These results provide important insights into the distribution, impact, and spread of CMD and SLCMV in Thailand.

Distribution and Molecular Diversity of Whitefly Species Colonizing Cassava in Kenya

The whitefly, Bemisia tabaci (Gennadium, Hemiptera) has been reported to transmit viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD) in many parts of sub-Saharan Africa (SSA). Currently, there is limited information on the distribution, species and haplotype composition of the whitefly populations colonizing cassava in Kenya. A study was conducted in the major cassava growing regions of Kenya to address this gap. Analyses of mitochondrial DNA cytochrome oxidase 1 (mtCO1) sequences revealed the presence of four distinct whitefly species: Bemisia tabaci, Bemisia afer, Aleurodicus dispersus and Paraleyrodes bondari in Kenya. The B. tabaci haplotypes were further resolved into SSA1, SSA2 and Indian Ocean (IO) putative species. The SSA1 population had three haplogroups of SSA1-SG1, SSA-SG2 and SSA1-SG3. Application of KASP genotyping grouped the Bemisia tabaci into two haplogroups namely sub-Saharan Africa East and Southern Africa (SSA-ESA) and sub-Saharan Africa East and Central Africa (SSA-ECA). The study presents the first report of P. bondari (Bondar’s nesting whitefly) on cassava in Kenya. Bemisia tabaci was widely distributed in all the major cassava growing regions in Kenya. The increased detection of different whitefly species on cassava and genetically diverse B. tabaci mitotypes indicates a significant influence on the dynamics of cassava virus epidemics in the field. The study highlights the need for continuous monitoring of invasive whitefly species population on cassava for timely application of management practices to reduce the impact of cassava viral diseases and prevent potential yield losses.

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, Ugandan 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 8,820 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.

Genetic Analysis of Yield and Quality Characteristics in Provitamin A Hybrid Cassava Families in Ghana

Cassava is the most widely cultivated and consumed crop in Ghana. Malnutrition is endemic in cassava-producing regions of Africa, partly due to the low micronutrient content of this crop. The aim of this study was to generate genetic information on characteristics such as total carotenoid content, dry matter content, root weight and number, and cassava mosaic disease (CMD), and their possible combination in cassava clones, using a North Carolina II breeding scheme. Five genetically diverse yellow-fleshed clones at advanced selection stages, with CMD resistance, were used as females and two high dry matter content white-fleshed clones, selected from farmers’ fields in Ghana, were used as males. Ten F1 families were generated, and evaluated at two locations in Ghana. General combining ability (GCA) mean squares were larger than specific combining ability (SCA) mean squares for harvest index, CMD, and carotenoid content, indicating additive genetic effects. The positive significant correlations that were observed between pulp color and carotenoid content; carotenoid content and CMD; pulp color and CMD; and pulp color and cortex color, make the screening of large numbers of progenies possible in the cassava breeding program. This could allow breeders to combine carotenoid content and CMD at the early breeding stages through visual assessment of pulp color and CMD symptoms. Large numbers of genotypes can be evaluated and a few can be selected to be quantified for carotenoid content at later stages of the breeding program, to save costs. One of the parents (P6), showed positive GCA effects for carotenoid content, dry matter content, CMD and storage root weight, hence could be used as a parent to generate clones that combine carotenoid content and dry matter content.

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.