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

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

Genetic complexity of cassava brown streak disease: insights from qPCR-based viral titer analysis and genome-wide association studies.

Cassava, a vital global food source, faces a threat from Cassava Brown Streak Disease (CBSD). CBSD results from two viruses: Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV). These viruses frequently pose challenges to the traditional symptom-based 1-5 phenotyping method due to its limitations in terms of accuracy and objectivity. Quantitative polymerase chain reaction (qPCR) offers precise virus quantification, although high costs hinder its widespread adoption. In this research, we utilized qPCR to measure the viral titer/load of CBSV and UCBSV. The objectives were to evaluate titer variability within the Cycle 2 (C2) population in two different environments, establish connections between viral titers and CBSD severity scores from the 1-5 scoring method, perform Genome-Wide Association Studies (GWAS) to identify genomic regions associated with CBSV and UCBSV titers, and investigate the functional annotated genes. The results demonstrated a significantly higher prevalence of CBSV (50.2%) in clones compared to UCBSV (12.9%) with mixed infections in some cases. Genotypic effects, particularly concerning UCBSV, were significant, with genotype-by-environment effects primarily influencing CBSV titer. GWAS Studies identified genomic regions associated with CBSV and UCBSV titers. Twenty-one SNP markers on chromosomes 10, 13, 17, and 18 exhibited significant associations with CBSV titer, collectively explaining 43.14% of the phenotypic variation. Additionally, 25 SNP markers on chromosomes 1, 2, 4, 5, 8, 11, 12, 13, 16, and 18 were associated with UCBSV titer, and explained 70.71% of the phenotypic variation. No shared genomic regions were identified between CBSV and UCBSV viral titers. Gene ontology analysis also revealed diverse gene functions, especially in transport and catalytic activities. These findings enhance our understanding of virus prevalence, genetics, and molecular functions in cassava plants, offering valuable insights for targeted breeding strategies.

Rapid analysis of hydrogen cyanide in fresh cassava roots using NIRSand machine learning algorithms: Meeting end user demand for low cyanogenic cassava.

This study focuses on meeting end-users' demand for cassava (Manihot esculenta Crantz) varieties with low cyanogenic potential (hydrogen cyanide potential [HCN]) by using near-infrared spectrometry (NIRS). This technology provides a fast, accurate, and reliable way to determine sample constituents with minimal sample preparation. The study aims to evaluate the effectiveness of machine learning (ML) algorithms such as logistic regression (LR), support vector machine (SVM), and partial least squares discriminant analysis (PLS-DA) in distinguishing between low and high HCN accessions. Low HCN accessions averagely scored 1-5.9, while high HCN accessions scored 6-9 on a 1-9 categorical scale. The researchers used 1164 root samples to test different NIRS prediction models and six spectral pretreatments. The wavelengths 961, 1165, 1403-1505, 1913-1981, and 2491 nm were influential in discrimination of low and high HCN accessions. Using selected wavelengths, LR achieved 100% classification accuracy and PLS-DA achieved 99% classification accuracy. Using the full spectrum, the best model for discriminating low and high HCN accessions was the PLS-DA combined with standard normal variate with second derivative, which produced an accuracy of 99.6%. The SVM and LR had moderate classification accuracies of 75% and 74%, respectively. This study demonstrates that NIRS coupled with ML algorithms can be used to identify low and high HCN accessions, which can help cassava breeding programs to select for low HCN accessions.

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.

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.

Alternative Approaches for Assessing Cassava Brown Streak Root Necrosis to Guide Resistance Breeding and Selection.

Cassava brown streak disease (CBSD) caused by the rapidly evolving cassava brown streak viruses (CBSVs), causes immense yield losses to the cassava value chain in eastern and southern Africa. Western Africa, another region that heavily depends on cassava is under eminent threat from CBSD. Resistance breeding is the best practical solution. However, complexities associated with CBSD resistance screening i.e., variable root sampling units, limit systematic attainment of genetic progress. Accordingly, we compared efficiency of five CBSD root necrosis assessment methods to guide selection: cassava brown streak disease root incidence (CBSDRi), cassava brown streak disease root severity (CBSDRs), cassava brown streak disease root severity computed as harmonic mean (CBSD-Harmonic), proportion-based root necrosis index (CBSD-proportion), and standardized root necrosis index (CBSD-standardized). The indexes (CBSD-proportion and CBSD-standardized) correct for variable sample size. We analyzed CBSD root necrosis data of 256 clones evaluated across 12 environments. Higher and variable standard errors were associated with root severity score 1 (no CBSD root necrosis). Lowest and highest plot-based heritability were respectively registered for CBSD-standardized (0.22) and CBSD-proportion (0.71). CBSDRs was only positively correlated with CBSDRi (r = 0.92) and CBSD-Harmonic (r = 0.97). Using best linear unbiased predictions (BLUPs), we ranked the top 15 CBSD resistant clones; only one clone (UG130014) featured in all the five assessment methods; two clones (UG130006 and UG120156) featured in four (CBSD-Harmonic, CBSDRi, CBSDRs, and CBSD-standardized); and five clones (UG120180, UG120063, UG130002, UG130033, and UG120183) featured in three methods (CBSD-Harmonic, CBSDRi, and CBSDRs). Influence of sample size was also quantified by sub-setting and analyzing CBSDRs data to have plots with at least 40 or 30 roots. Data stabilization was evident in plots with 30 roots. The significant influence of root sample sizes on overall ranking of clones, justifies the use of CBSD root necrosis indexes in early selection stages i.e., seedling and/or clonal trials, that are often characterized by high variations in roots assessed per plot. It is expected that this information will provide a foundation for harmonizing and/or optimizing on-going and future CBSD resistance breeding efforts.