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.