ABSTRACT
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.
ABSTRACT
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.