Genome-Wide Association and Prediction Reveals Genetic Architecture of Cassava Mosaic Disease Resistance and Prospects for Rapid Genetic Improvement

Identification of Genomic Regions for Traits Associated with Flowering in Cassava (Manihot esculenta Crantz)

Flowering in cassava (Manihot esculenta Crantz) is crucial for the generation of botanical seed for breeding. However, genotypes preferred by most farmers are erect and poor at flowering or never flower. To elucidate the genetic basis of flowering, 293 diverse cassava accessions were evaluated for flowering-associated traits at two locations and seasons in Uganda. Genotyping using the Diversity Array Technology Pty Ltd. (DArTseq) platform identified 24,040 single-nucleotide polymorphisms (SNPs) distributed on the 18 cassava chromosomes. Population structure analysis using principal components (PCs) and kinships showed three clusters; the first five PCs accounted for 49.2% of the observed genetic variation. Linkage disequilibrium (LD) estimation averaged 0.32 at a distance of ~2850 kb (kilo base pairs). Polymorphism information content (PIC) and minor allele frequency (MAF) were 0.25 and 0.23, respectively. A genome-wide association study (GWAS) analysis uncovered 53 significant marker-trait associations (MTAs) with flowering-associated traits involving 27 loci. Two loci, SNPs S5_29309724 and S15_11747301, were associated with all the traits. Using five of the 27 SNPs with a Phenotype_Variance_Explained (PVE) ≥ 5%, 44 candidate genes were identified in the peak SNP sites located within 50 kb upstream or downstream, with most associated with branching traits. Eight of the genes, orthologous to Arabidopsis and other plant species, had known functional annotations related to flowering, e.g., eukaryotic translation initiation factor and myb family transcription factor. This study identified genomic regions associated with flowering-associated traits in cassava, and the identified SNPs can be useful in marker-assisted selection to overcome hybridization challenges, like unsynchronized flowering, and candidate gene validation.

Cassava molecular genetics and genomics for enhanced resistance to diseases and pests

Cassava (Manihot esculenta) is one of the most important sources of dietary calories in the tropics, playing a central role in food and economic security for smallholder farmers. Cassava production is highly constrained by several pests and diseases, mostly cassava mosaic disease (CMD) and cassava brown streak disease (CBSD). These diseases cause significant yield losses, affecting food security and the livelihoods of smallholder farmers. Developing resistant varieties is a good way of increasing cassava productivity. Although some levels of resistance have been developed for some of these diseases, there is observed breakdown in resistance for some diseases, such as CMD. A frequent re-evaluation of existing disease resistance traits is required to make sure they are still able to withstand the pressure associated with pest and pathogen evolution. Modern breeding approaches such as genomic-assisted selection in addition to biotechnology techniques like classical genetic engineering or genome editing can accelerate the development of pest- and disease-resistant cassava varieties. This article summarizes current developments and discusses the potential of using molecular genetics and genomics to produce cassava varieties resistant to diseases and pests.

The Development of Thematic Core Collections in Cassava Based on Yield, Disease Resistance, and Root Quality Traits

Thematic collections (TCs), which are composed of genotypes with superior agronomic traits and reduced size, offer valuable opportunities for parental selection in plant breeding programs. Three TCs were created to focus on crucial attributes: root yield (CC_Yield), pest and disease resistance (CC_Disease), and root quality traits (CC_Root_quality). The genotypes were ranked using the best linear unbiased predictors (BLUP) method, and a truncated selection was implemented for each collection based on specific traits. The TCs exhibited minimal overlap, with each collection comprising 72 genotypes (CC_Disease), 63 genotypes (CC_Root_quality), and 64 genotypes (CC_Yield), representing 4%, 3.5%, and 3.5% of the total individuals in the entire collection, respectively. The Shannon-Weaver Diversity Index values generally varied but remained below 10% when compared to the entire collection. Most TCs exhibited observed heterozygosity, genetic diversity, and the inbreeding coefficient that closely resembled those of the entire collection, effectively retaining 90.76%, 88.10%, and 88.99% of the alleles present in the entire collection (CC_Disease, CC_Root_quality, and CC_Disease, respectively). A PCA of molecular and agro-morphological data revealed well-distributed and dispersed genotypes, while a discriminant analysis of principal components (DAPC) displayed a high discrimination capacity among the accessions within each collection. The strategies employed in this study hold significant potential for advancing crop improvement efforts.

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.

Increasing cassava root yield: Additive-dominant genetic models for selection of parents and clones

Genomic selection has been promising in situations where phenotypic assessments are expensive, laborious, and/or inefficient. This work evaluated the efficiency of genomic prediction methods combined with genetic models in clone and parent selection with the goal of increasing fresh root yield, dry root yield, as well as dry matter content in cassava roots. The bias and predictive ability of the combinations of prediction methods Genomic Best Linear Unbiased Prediction (G-BLUP), Bayes B, Bayes Cπ, and Reproducing Kernel Hilbert Spaces with additive and additive-dominant genetic models were estimated. Fresh and dry root yield exhibited predominantly dominant heritability, while dry matter content exhibited predominantly additive heritability. The combination of prediction methods and genetic models did not show significant differences in the predictive ability for dry matter content. On the other hand, the prediction methods with additive-dominant genetic models had significantly higher predictive ability than the additive genetic models for fresh and dry root yield, allowing higher genetic gains in clone selection. However, higher predictive ability for genotypic values did not result in differences in breeding value predictions between additive and additive-dominant genetic models. G-BLUP with the classical additive-dominant genetic model had the best predictive ability and bias estimates for fresh and dry root yield. For dry matter content, the highest predictive ability was obtained by G-BLUP with the additive genetic model. Dry matter content exhibited the highest heritability, predictive ability, and bias estimates compared with other traits. The prediction methods showed similar selection gains with approximately 67% of the phenotypic selection gain. By shortening the breeding cycle time by 40%, genomic selection may overcome phenotypic selection by 10%, 13%, and 18% for fresh root yield, dry root yield, and dry matter content, respectively, with a selection proportion of 15%. The most suitable genetic model for each trait allows for genomic selection optimization in cassava with high selection gains, thereby accelerating the release of new varieties.

Analysis of proteomic changes in cassava cv. Kasetsart 50 caused by Sri Lankan cassava mosaic virus infection

BACKGROUND

Sri Lankan cassava mosaic virus (SLCMV) is a plant virus causing significant economic losses throughout Southeast Asia. While proteomics has the potential to identify molecular markers that could assist the breeding of virus resistant cultivars, the effects of SLCMV infection in cassava have not been previously explored in detail.

RESULTS

Liquid Chromatography-Tandem Mass Spectrometry (LC/MS-MS) was used to identify differentially expressed proteins in SLCMV infected leaves, and qPCR was used to confirm changes at mRNA levels. LC/MS-MS identified 1,813 proteins, including 479 and 408 proteins that were upregulated in SLCMV-infected and healthy cassava plants respectively, while 109 proteins were detected in both samples. Most of the identified proteins were involved in biosynthetic processes (29.8%), cellular processes (20.9%), and metabolism (18.4%). Transport proteins, stress response molecules, and proteins involved in signal transduction, plant defense responses, photosynthesis, and cellular respiration, although present, only represented a relatively small subset of the detected differences. RT-qPCR confirmed the upregulation of WRKY 77 (A0A140H8T1), WRKY 83 (A0A140H8T7), NAC 6 (A0A0M4G3M4), NAC 35 (A0A0M5JAB4), NAC 22 (A0A0M5J8Q6), NAC 54 (A0A0M4FSG8), NAC 70 (A0A0M4FEU9), MYB (A0A2C9VER9 and A0A2C9VME6), bHLH (A0A2C9UNL9 and A0A2C9WBZ1) transcription factors. Additional upregulated transcripts included receptors, such as receptor-like serine/threonine-protein kinase (RSTK) (A0A2C9UPE4), Toll/interleukin-1 receptor (TIR) (A0A2C9V5Q3), leucine rich repeat N-terminal domain (LRRNT_2) (A0A2C9VHG8), and cupin (A0A199UBY6). These molecules participate in innate immunity, plant defense mechanisms, and responses to biotic stress and to phytohormones.

CONCLUSIONS

We detected 1,813 differentially expressed proteins infected cassava plants, of which 479 were selectively upregulated. These could be classified into three main biological functional groups, with roles in gene regulation, plant defense mechanisms, and stress responses. These results will help identify key proteins affected by SLCMV infection in cassava plants.

Genome-wide association studies reveal genetic basis of ionomic variation in cassava

As one of the most important food crops, cassava (Manihot esculenta) is the main dietary source of micronutrients for about 1 billion people. However, the ionomic variation in cassava and the underlying genetic mechanisms remain unclear so far. Herein, genome-wide association studies were performed to reveal the specific single nucleotide polymorphisms (SNPs) that affect the ionomic variation in cassava. We identified 164 SNPs with P-values lower than the threshold located in 88 loci associated with divergent ionomic variations. Among them, 13 SNPs are related to both calcium (Ca) and magnesium (Mg), and many loci for different ionomic traits seem to be clustered on specific chromosome regions. Moreover, we identified the peak SNPs in the promoter regions of Sc10g003170 (encoding methionyl-tRNA synthetase [MetRS]) and Sc18g015190 (encoding the transcriptional regulatory protein AlgP) for nitrogen (N) and phosphorus (P) accumulation, respectively. Notably, these two SNPs (chr10_32807962 and chr18_31343738) were directly correlated with the transcript levels of Sc10g003170 (MetRS) and Sc18g015190 (AlgP), which positively modulated N accumulation and P concentration in cassava, respectively. Taken together, this study provides important insight into the genetic basis of cassava natural ionomic variation, which will promote genetic breeding to improve nutrient use and accumulation of elements in cassava.

Validation of KASP markers associated with cassava mosaic disease resistance, storage root dry matter and provitamin A carotenoid contents in Ugandan cassava germplasm

INTRODUCTION

The intrinsic high heterozygosity of cassava makes conventional breeding ineffective for rapid genetic improvement. However, recent advances in next generation sequencing technologies have enabled the use of high-density markers for genome-wide association studies, aimed at identifying single nucleotide polymorphisms (SNPs) linked to major traits such as cassava mosaic disease (CMD) resistance, dry matter content (DMC) and total carotenoids content (TCC). A number of these trait-linked SNPs have been converted to Kompetitive allele-specific polymerase chain reaction (KASP) markers for downstream application of marker assisted selection.

METHODS

We assayed 13 KASP markers to evaluate their effectiveness in selecting for CMD, DMC and TCC in 1,677 diverse cassava genotypes representing two independent breeding populations in Uganda.

RESULTS

Five KASP markers had significant co-segregation with phenotypes; CMD resistance (2), DMC (1) and TCC (2), with each marker accounting for at least 30% of the phenotypic variation. Markers located within the chromosomal regions for which strong marker-trait association loci have been characterised (chromosome 12 markers for CMD, chromosome 1 markers for DMC and TCC) had consistently superior ability to discriminate the respective phenotypes.

DISCUSSION

The results indicate varying discriminatory abilities of the KASP markers assayed and the need for their context-based use for MAS, with PSY2_572 particularly effective in selecting for high TCC. Availing the effective KASP markers on cost-effective genotyping platforms could facilitate practical implementation of marker-assisted cassava breeding for accelerated genetic gains for CMD, DMC and provitamin A carotenoids.

Utility of Ugandan genomic selection cassava breeding populations for prediction of cassava viral disease resistance and yield in West African clones

Cassava (Manihot esculenta Crantz) is a staple crop for ~800 million people in sub-Saharan Africa. Its production and productivity are being heavily affected by the two viral diseases: cassava brown streak disease (CBSD) and cassava mosaic disease (CMD), impacting greatly on edible root yield. CBSD is currently endemic to central, eastern and southern Africa, if not contained could spread to West Africa the largest cassava producer and consumer in the continent. Genomic selection (GS) has been implemented in Ugandan cassava breeding for accelerated development of virus resistant and high yielding clones. This study leveraged available GS training data in Uganda for pre-emptive CBSD breeding in W. Africa alongside CMD and fresh root yield (FRW). First, we tracked genetic gain through the current three cycles of GS in Uganda. The mean genomic estimated breeding values (GEBVs), indicated general progress from initial cycle zero (C0) to cycle one (C1) and cycle two (C2) for CBSD traits and yield except for CMD. Secondly, we used foliar data of both CBSD and CMD, as well as harvest root necrosis and yield data to perform cross-validation predictions. Cross-validation prediction accuracies of five GS models were tested for each of the three GS cycles and West African (WA) germplasm as a test set. In all cases, cross-validation prediction accuracies were low to moderate, ranging from -0.16 to 0.68 for CBSD traits, -0.27 to 0.57 for CMD and -0.22 to 0.41 for fresh root weight (FRW). Overall, the highest prediction accuracies were recorded in C0 for all traits tested across models and the best performing model in cross-validation was G-BLUP. Lastly, we tested the predictive ability of the Ugandan training sets to predict CBSD in W. African clones. In general, the Ugandan training sets had low prediction accuracies for all traits across models in West African germplasm, varying from -0.18 to 0.1. Based on the findings of this study, the cassava breeding program in Uganda has made progress through application of GS for most target traits, but the utility of the training population for pre-emptive breeding in WA is limiting. In this case, efforts should be devoted to sharing Ugandan germplasm that possess resistance with the W. African breeding programs for hybridization to fully enable deployment of genomic selection as a pre-emptive CBSD breeding strategy in W. Africa.

Validation of KASP-SNP markers in cassava germplasm for marker-assisted selection of increased carotenoid content and dry matter content

Provitamin A biofortification and increased dry matter content are important breeding targets in cassava improvement programs worldwide. Biofortified varieties contribute to the alleviation of provitamin A deficiency, a leading cause of preventable blindness common among pre-school children and pregnant women in developing countries particularly Africa. Dry matter content is a major component of dry yield and thus underlies overall variety performance and acceptability by growers, processors, and consumers. Single nucleotide polymorphism (SNP) markers linked to these traits have recently been discovered through several genome-wide association studies but have not been deployed for routine marker-assisted selection (MAS). This is due to the lack of useful information on markers' performances in diverse genetic backgrounds. To overcome this bottleneck, technical and biological validation of the loci associated with increased carotenoid content and dry matter content were carried out using populations independent of the marker discovery population. In the present study, seven previously identified markers for these traits were converted to a robust set of uniplex allele-specific polymerase chain reaction (PCR) assays and validated in two independent pre-breeding and breeding populations. These assays were efficient in discriminating marker genotypic classes and had an average call rate greater than 98%. A high correlation was observed between the predicted and observed carotenoid content as inferred by root yellowness intensity in the breeding (r = 0.92) and pre-breeding (r = 0.95) populations. On the other hand, dry matter content-markers had moderately low predictive accuracy in both populations (r< 0.40) due to the more quantitative nature of the trait. This work confirmed the markers' effectiveness in multiple backgrounds, therefore, further strengthening their value in cassava biofortification to ensure nutritional security as well as dry matter content productivity. Our study provides a framework to guide future marker validation, thus leading to the more routine use of markers in MAS in cassava improvement programs.

CRISPR for accelerating genetic gains in under-utilized crops of the drylands: Progress and prospects

Technologies and innovations are critical for addressing the future food system needs where genetic resources are an essential component of the change process. Advanced breeding tools like "genome editing" are vital for modernizing crop breeding to provide game-changing solutions to some of the "must needed" traits in agriculture. CRISPR/Cas-based tools have been rapidly repurposed for editing applications based on their improved efficiency, specificity and reduced off-target effects. Additionally, precise gene-editing tools such as base editing, prime editing, and multiplexing provide precision in stacking of multiple traits in an elite variety, and facilitating specific and targeted crop improvement. This has helped in advancing research and delivery of products in a short time span, thereby enhancing the rate of genetic gains. A special focus has been on food security in the drylands through crops including millets, teff, fonio, quinoa, Bambara groundnut, pigeonpea and cassava. While these crops contribute significantly to the agricultural economy and resilience of the dryland, improvement of several traits including increased stress tolerance, nutritional value, and yields are urgently required. Although CRISPR has potential to deliver disruptive innovations, prioritization of traits should consider breeding product profiles and market segments for designing and accelerating delivery of locally adapted and preferred crop varieties for the drylands. In this context, the scope of regulatory environment has been stated, implying the dire impacts of unreasonable scrutiny of genome-edited plants on the evolution and progress of much-needed technological advances.

Advances in integrated genomic selection for rapid genetic gain in crop improvement: a review

Genomic selection and its importance in crop breeding. Integration of GS with new breeding tools and developing SOP for GS to achieve maximum genetic gain with low cost and time. The success of conventional breeding approaches is not sufficient to meet the demand of a growing population for nutritious food and other plant-based products. Whereas, marker assisted selection (MAS) is not efficient in capturing all the favorable alleles responsible for economic traits in the process of crop improvement. Genomic selection (GS) developed in livestock breeding and then adapted to plant breeding promised to overcome the drawbacks of MAS and significantly improve complicated traits controlled by gene/QTL with small effects. Large-scale deployment of GS in important crops, as well as simulation studies in a variety of contexts, addressed G × E interaction effects and non-additive effects, as well as lowering breeding costs and time. The current study provides a complete overview of genomic selection, its process, and importance in modern plant breeding, along with insights into its application. GS has been implemented in the improvement of complex traits including tolerance to biotic and abiotic stresses. Furthermore, this review hypothesises that using GS in conjunction with other crop improvement platforms accelerates the breeding process to increase genetic gain. The objective of this review is to highlight the development of an appropriate GS model, the global open source network for GS, and trans-disciplinary approaches for effective accelerated crop improvement. The current study focused on the application of data science, including machine learning and deep learning tools, to enhance the accuracy of prediction models. Present study emphasizes on developing plant breeding strategies centered on GS combined with routine conventional breeding principles by developing GS-SOP to achieve enhanced genetic gain.

PP2C1 fine-tunes melatonin biosynthesis and phytomelatonin receptor PMTR1 binding to melatonin in cassava

Melatonin is an important molecule in both animals and plants, regulating circadian rhythms and stress responses. Therefore, the improvement of melatonin accumulation not only strengthens the function of melatonin but also improves stress resistance in crops. Although melatonin biosynthetic enzymes have been identified through reverse genetics previously, an investigation of melatonin level-related genes through forward genetics in plants has yet to be performed. In this study, a genome-wide association study using cassava natural population of 298 genetic resources identified melatonin accumulation 1 (MA1), which regulates the natural variation of melatonin levels in cassava. We found that MA1 encodes type 2C protein phosphatase 1 (PP2C1), which serves as a negative regulator of melatonin levels in cassava. MePP2C1 physically interacts with MeRAV1/2 and MeWRKY20 and dephosphorylates them at serine (S) 35 residue, S34 residue, and S176 residue, respectively, thereby hindering their transcriptional activation on downstream melatonin biosynthetic genes. Notably, MePP2C1 interacts with phytomelatonin receptor MePMTR1 and dephosphorylates it at S11 residue, repressing its binding to melatonin. In summary, this study demonstrates that MePP2C1 as MA1 plays dual roles in negatively regulating both melatonin accumulation and signaling, extending the understanding of the molecular mechanism underlying melatonin accumulation and signaling through forward genetics in plants.

Mutations in DNA polymerase δ subunit 1 co-segregate with CMD2-type resistance to Cassava Mosaic Geminiviruses

Cassava mosaic disease (CMD) suppresses cassava yields across the tropics. The dominant CMD2 locus confers resistance to cassava mosaic geminiviruses. It has been reported that CMD2-type landraces lose resistance after regeneration through de novo morphogenesis. As full genome bisulfite sequencing failed to uncover an epigenetic mechanism for this loss of resistance, whole genome sequencing and genetic variant analysis was performed and the CMD2 locus was fine-mapped to a 190 kilobase interval. Collectively, these data indicate that CMD2-type resistance is caused by a nonsynonymous, single nucleotide polymorphism in DNA polymerase δ subunit 1 (MePOLD1) located within this region. Virus-induced gene silencing of MePOLD1 in a CMD-susceptible cassava variety produced a recovery phenotype typical of CMD2-type resistance. Analysis of other CMD2-type cassava varieties identified additional candidate resistance alleles within MePOLD1. Genetic variation of MePOLD1, therefore, could represent an important genetic resource for resistance breeding and/or genome editing, and elucidating mechanisms of resistance to geminiviruses.

Current status and impending progress for cassava structural genomics

We demystify recent advances in genome assemblies for the heterozygous staple crop cassava (Manihot esculenta), and highlight key cassava genomic resources. Cassava, Manihot esculenta Crantz, is a crop of societal and agricultural importance in tropical regions around the world. Genomics provides a platform for accelerated improvement of cassava's nutritional and agronomic traits, as well as for illuminating aspects of cassava's history including its path towards domestication. The highly heterozygous nature of the cassava genome is widely recognized. However, the full extent and context of this heterozygosity has been difficult to reveal because of technological limitations within genome sequencing. Only recently, with several new long-read sequencing technologies coming online, has the genomics community been able to tackle some similarly difficult genomes. In light of these recent advances, we provide this review to document the current status of the cassava genome and genomic resources and provide a perspective on what to look forward to in the coming years.

Genome-wide association analysis reveals new insights into the genetic architecture of defensive, agro-morphological and quality-related traits in cassava

More than 40 QTLs associated with 14 stress-related, quality and agro-morphological traits were identified. A catalogue of favourable SNP markers for MAS and a list of candidate genes are provided. Cassava (Manihot esculenta) is one of the most important starchy root crops in the tropics due to its adaptation to marginal environments. Genetic progress in this clonally propagated crop can be accelerated through the discovery of markers and candidate genes that could be used in cassava breeding programs. We carried out a genome-wide association study (GWAS) using a panel of 5130 clones developed at the International Institute of Tropical Agriculture-Nigeria. The population was genotyped at more than 100,000 SNP markers via genotyping-by-sequencing (GBS). Genomic regions underlying genetic variation for 14 traits classified broadly into four categories: biotic stress (cassava mosaic disease and cassava green mite severity); quality (dry matter content and carotenoid content) and plant agronomy (harvest index and plant type) were investigated. We also included several agro-morphological traits related to leaves, stems and roots with high heritability. In total, 41 significant associations were uncovered. While some of the identified loci matched with those previously reported, we present additional association signals for the traits. We provide a catalogue of favourable alleles at the most significant SNP for each trait-locus combination and candidate genes occurring within the GWAS hits. These resources provide a foundation for the development of markers that could be used in cassava breeding programs and candidate genes for functional validation.

Genome-wide association study of cassava starch paste properties

An understanding of cassava starch paste properties (CSPP) can contribute to the selection of clones with differentiated starches. This study aimed to identify genomic regions associated with CSPP using different genome-wide association study (GWAS) methods (MLM, MLMM, and Farm-CPU). The GWAS was performed using 23,078 single-nucleotide polymorphisms (SNPs). The rapid viscoanalyzer (RVA) parameters were pasting temperature (PastTemp), peak viscosity (PeakVisc), hot-paste viscosity (Hot-PVisc), cool-paste viscosity (Cold-PVisc), final viscosity (FinalVis), breakdown (BreDow), and setback (Setback). Broad phenotypic and molecular diversity was identified based on the genomic kinship matrix. The broad-sense heritability estimates (h²) ranged from moderate to high magnitudes (0.66 to 0.76). The linkage disequilibrium (LD) declined to between 0.3 and 2.0 Mb (r² <0.1) for most chromosomes, except chromosome 17, which exhibited an extensive LD. Thirteen SNPs were found to be significantly associated with CSPP, on chromosomes 3, 8, 17, and 18. Only the BreDow trait had no associated SNPs. The regional marker-trait associations on chromosome 18 indicate a LD block between 2907312 and 3567816 bp and that SNP S18_3081635 was associated with SetBack, FinalVis, and Cold-PVisc (all three GWAS methods) and with Hot-PVisc (MLM), indicating that this SNP can track these four traits simultaneously. The variance explained by the SNPs ranged from 0.13 to 0.18 for SetBack, FinalVis, and Cold-PVisc and from 0.06 to 0.09 for PeakVisc and Hot-PVisc. The results indicated additive effects of the genetic control of Cold-PVisc, FinalVis, Hot-PVisc, and SetBack, especially on the large LD block on chromosome 18. One transcript encoding the glycosyl hydrolase family 35 enzymes on chromosome 17 and one encoding the mannose-p-dolichol utilization defect 1 protein on chromosome 18 were the most likely candidate genes for the regulation of CSPP. These results underline the potential for the assisted selection of high-value starches to improve cassava root quality through breeding programs.

Two decades of association mapping: Insights on disease resistance in major crops

Climate change across the globe has an impact on the occurrence, prevalence, and severity of plant diseases. About 30% of yield losses in major crops are due to plant diseases; emerging diseases are likely to worsen the sustainable production in the coming years. Plant diseases have led to increased hunger and mass migration of human populations in the past, thus a serious threat to global food security. Equipping the modern varieties/hybrids with enhanced genetic resistance is the most economic, sustainable and environmentally friendly solution. Plant geneticists have done tremendous work in identifying stable resistance in primary genepools and many times other than primary genepools to breed resistant varieties in different major crops. Over the last two decades, the availability of crop and pathogen genomes due to advances in next generation sequencing technologies improved our understanding of trait genetics using different approaches. Genome-wide association studies have been effectively used to identify candidate genes and map loci associated with different diseases in crop plants. In this review, we highlight successful examples for the discovery of resistance genes to many important diseases. In addition, major developments in association studies, statistical models and bioinformatic tools that improve the power, resolution and the efficiency of identifying marker-trait associations. Overall this review provides comprehensive insights into the two decades of advances in GWAS studies and discusses the challenges and opportunities this research area provides for breeding resistant varieties.

Genome-wide association mapping and genomic prediction of yield-related traits and starch pasting properties in cassava

GWAS identified eight yield-related, peak starch type of waxy and wild-type starch and 21 starch pasting property-related traits (QTLs). Prediction ability of eight GS models resulted in low to high predictability, depending on trait, heritability, and genetic architecture. Cassava is both a food and an industrial crop in Africa, South America, and Asia, but knowledge of the genes that control yield and starch pasting properties remains limited. We carried out a genome-wide association study to clarify the molecular mechanisms underlying these traits and to explore marker-based breeding approaches. We estimated the predictive ability of genomic selection (GS) using parametric, semi-parametric, and nonparametric GS models with a panel of 276 cassava genotypes from Thai Tapioca Development Institute, International Center for Tropical Agriculture, International Institute of Tropical Agriculture, and other breeding programs. The cassava panel was genotyped via genotyping-by-sequencing, and 89,934 single-nucleotide polymorphism (SNP) markers were identified. A total of 31 SNPs associated with yield, starch type, and starch properties traits were detected by the fixed and random model circulating probability unification (FarmCPU), Bayesian-information and linkage-disequilibrium iteratively nested keyway and compressed mixed linear model, respectively. GS models were developed, and forward predictabilities using all the prediction methods resulted in values of - 0.001-0.71 for the four yield-related traits and 0.33-0.82 for the seven starch pasting property traits. This study provides additional insight into the genetic architecture of these important traits for the development of markers that could be used in cassava breeding programs.

Genome-Wide Association Study of Root Mealiness and Other Texture-Associated Traits in Cassava

Cassava breeders have made significant progress in developing new genotypes with improved agronomic characteristics such as improved root yield and resistance against biotic and abiotic stresses. However, these new and improved cassava (Manihot esculenta Crantz) varieties in cultivation in Nigeria have undergone little or no improvement in their culinary qualities; hence, there is a paucity of genetic information regarding the texture of boiled cassava, particularly with respect to its mealiness, the principal sensory quality attribute of boiled cassava roots. The current study aimed at identifying genomic regions and polymorphisms associated with natural variation for root mealiness and other texture-related attributes of boiled cassava roots, which includes fibre, adhesiveness (ADH), taste, aroma, colour, and firmness. We performed a genome-wide association (GWAS) analysis using phenotypic data from a panel of 142 accessions obtained from the National Root Crops Research Institute (NRCRI), Umudike, Nigeria, and a set of 59,792 high-quality single nucleotide polymorphisms (SNPs) distributed across the cassava genome. Through genome-wide association mapping, we identified 80 SNPs that were significantly associated with root mealiness, fibre, adhesiveness, taste, aroma, colour and firmness on chromosomes 1, 4, 5, 6, 10, 13, 17 and 18. We also identified relevant candidate genes that are co-located with peak SNPs linked to these traits in M. esculenta. A survey of the cassava reference genome v6.1 positioned the SNPs on chromosome 13 in the vicinity of Manes.13G026900, a gene recognized as being responsible for cell adhesion and for the mealiness or crispness of vegetables and fruits, and also known to play an important role in cooked potato texture. This study provides the first insights into understanding the underlying genetic basis of boiled cassava root texture. After validation, the markers and candidate genes identified in this novel work could provide important genomic resources for use in marker-assisted selection (MAS) and genomic selection (GS) to accelerate genetic improvement of root mealiness and other culinary qualities in cassava breeding programmes in West Africa, especially in Nigeria, where the consumption of boiled and pounded cassava is low.

Can Cross-Country Genomic Predictions Be a Reasonable Strategy to Support Germplasm Exchange? - A Case Study With Hydrogen Cyanide in Cassava

Genomic prediction (GP) offers great opportunities for accelerated genetic gains by optimizing the breeding pipeline. One of the key factors to be considered is how the training populations (TP) are composed in terms of genetic improvement, kinship/origin, and their impacts on GP. Hydrogen cyanide content (HCN) is a determinant trait to guide cassava's products usage and processing. This work aimed to achieve the following objectives: (i) evaluate the feasibility of using cross-country (CC) GP between germplasm's of Embrapa Mandioca e Fruticultura (Embrapa, Brazil) and The International Institute of Tropical Agriculture (IITA, Nigeria) for HCN; (ii) provide an assessment of population structure for the joint dataset; (iii) estimate the genetic parameters based on single nucleotide polymorphisms (SNPs) and a haplotype-approach. Datasets of HCN from Embrapa and IITA breeding programs were analyzed, separately and jointly, with 1,230, 590, and 1,820 clones, respectively. After quality control, ∼14K SNPs were used for GP. The genomic estimated breeding values (GEBVs) were predicted based on SNP effects from analyses with TP composed of the following: (i) Embrapa genotypic and phenotypic data, (ii) IITA genotypic and phenotypic data, and (iii) the joint datasets. Comparisons on GEBVs' estimation were made considering the hypothetical situation of not having the phenotypic characterization for a set of clones for a certain research institute/country and might need to use the markers' effects that were trained with data from other research institutes/country's germplasm to estimate their clones' GEBV. Fixation index (F_ST) among the genetic groups identified within the joint dataset ranged from 0.002 to 0.091. The joint dataset provided an improved accuracy (0.8-0.85) compared to the prediction accuracy of either germplasm's sources individually (0.51-0.67). CC GP proved to have potential use under the present study's scenario, the correlation between GEBVs predicted with TP from Embrapa and IITA was 0.55 for Embrapa's germplasm, whereas for IITA's it was 0.1. This seems to be among the first attempts to evaluate the CC GP in plants. As such, a lot of useful new information was provided on the subject, which can guide new research on this very important and emerging field.

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.

Deciphering the Genetic Architecture of Plant Virus Resistance by GWAS, State of the Art and Potential Advances

Growing virus resistant varieties is a highly effective means to avoid yield loss due to infection by many types of virus. The challenge is to be able to detect resistance donors within plant species diversity and then quickly introduce alleles conferring resistance into elite genetic backgrounds. Until now, mainly monogenic forms of resistance with major effects have been introduced in crops. Polygenic resistance is harder to map and introduce in susceptible genetic backgrounds, but it is likely more durable. Genome wide association studies (GWAS) offer an opportunity to accelerate mapping of both monogenic and polygenic resistance, but have seldom been implemented and described in the plant-virus interaction context. Yet, all of the 48 plant-virus GWAS published so far have successfully mapped QTLs involved in plant virus resistance. In this review, we analyzed general and specific GWAS issues regarding plant virus resistance. We have identified and described several key steps throughout the GWAS pipeline, from diversity panel assembly to GWAS result analyses. Based on the 48 published articles, we analyzed the impact of each key step on the GWAS power and showcase several GWAS methods tailored to all types of viruses.

Genomic mating in outbred species: predicting cross usefulness with additive and total genetic covariance matrices

Diverse crops are both outbred and clonally propagated. Breeders typically use truncation selection of parents and invest significant time, land, and money evaluating the progeny of crosses to find exceptional genotypes. We developed and tested genomic mate selection criteria suitable for organisms of arbitrary homozygosity level where the full-sibling progeny are of direct interest as future parents and/or cultivars. We extended cross variance and covariance variance prediction to include dominance effects and predicted the multivariate selection index genetic variance of crosses based on haplotypes of proposed parents, marker effects, and recombination frequencies. We combined the predicted mean and variance into usefulness criteria for parent and variety development. We present an empirical study of cassava (Manihot esculenta), a staple tropical root crop. We assessed the potential to predict the multivariate genetic distribution (means, variances, and trait covariances) of 462 cassava families in terms of additive and total value using cross-validation. Most variance (89%) and covariance (70%) prediction accuracy estimates were greater than zero. The usefulness of crosses was accurately predicted with good correspondence between the predicted and the actual mean performance of family members breeders selected for advancement as new parents and candidate varieties. We also used a directional dominance model to quantify significant inbreeding depression for most traits. We predicted 47,083 possible crosses of 306 parents and contrasted them to those previously tested to show how mate selection can reveal the new potential within the germplasm. We enable breeders to consider the potential of crosses to produce future parents (progeny with top breeding values) and varieties (progeny with top own performance).

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.

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.

RNA polymerase mapping in plants identifies intergenic regulatory elements enriched in causal variants

Control of gene expression is fundamental at every level of cell function. Promoter-proximal pausing and divergent transcription at promoters and enhancers, which are prominent features in animals, have only been studied in a handful of research experiments in plants. PRO-Seq analysis in cassava (Manihot esculenta) identified peaks of transcriptionally engaged RNA polymerase at both the 5' and 3' end of genes, consistent with paused or slowly moving Polymerase. In addition, we identified divergent transcription at intergenic sites. A full genome search for bi-directional transcription using an algorithm for enhancer detection developed in mammals (dREG) identified many intergenic regulatory element (IRE) candidates. These sites showed distinct patterns of methylation and nucleotide conservation based on genomic evolutionary rate profiling (GERP). SNPs within these IRE candidates explained significantly more variation in fitness and root composition than SNPs in chromosomal segments randomly ascertained from the same intergenic distribution, strongly suggesting a functional importance of these sites. Maize GRO-Seq data showed RNA polymerase occupancy at IREs consistent with patterns in cassava. Furthermore, these IREs in maize significantly overlapped with sites previously identified on the basis of open chromatin, histone marks, and methylation, and were enriched for reported eQTL. Our results suggest that bidirectional transcription can identify intergenic genomic regions in plants that play an important role in transcription regulation and whose identification has the potential to aid crop improvement.

Fifty years of a public cassava breeding program: evolution of breeding objectives, methods, and decision-making processes

This paper reviews and analyzes key features from cassava breeding at the International Center for Tropical Agriculture (CIAT) over 50 years and draws lessons for public breeding efforts broadly. The breeding team, jointly with national program partners and the private processing sector, defined breeding objectives and guiding business plans. These have evolved through the decades and currently focus on four global product profiles. The recurrent selection method also evolved and included innovations such as estimation of phenotypic breeding values, increasing the number of locations in the first stage of agronomic evaluations, gradual reduction of the duration of breeding cycles (including rapid cycling for high-heritability traits), the development of protocols for the induction of flowering, and the introduction of genome-wide predictions. The impact of cassava breeding depends significantly on the type of target markets. When roots are used for large processing facilities for starch, animal feeding or ethanol production (such as in SE Asia), the adoption of improved varieties is nearly universal and productivity at the regional scale increases significantly. When markets and relevant infrastructure are weak or considerable proportion of the production goes for local artisanal processing and on-farm consumption, the impact has been lower. The potential of novel breeding tools needs to be properly assessed for the most effective allocation of resources. Finally, a brief summary of challenges and opportunities for the future of cassava breeding is presented. The paper describes multiple ways that public and private sector breeding programs can learn from each other to optimize success.

Outlook of Cassava Brown Streak Disease Assessment: Perspectives of the Screening Methods of Breeders and Pathologists

Cassava production and productivity in Eastern, Central, and Southern Africa are ravaged by cassava brown streak disease (CBSD), causing yield losses of up to 100% when susceptible varieties are grown. Efforts to develop CBSD-resistant clones are underway. However, the methods for screening CBSD resistance currently vary between breeders and pathologists, with the limited empirical data to support their choices. In this study, we used the empirical CBSD foliar and root necrosis data from two breeding populations, termed cycle zero (C₀) and cycle one (C₁), to assess and compare the effectiveness of the CBSD screening methods of breeders vs. pathologists. On the one hand, the estimates of broad-sense heritability (H ²) for the CBSD root necrosis assessment of breeder ranged from 0.15 to 0.87, while for the assessment method of pathologists, H ² varied from 0.00 to 0.71 in C₀ clones. On the other hand, the marker-based heritability estimates (h ²) for C₀ ranged from 0.00 to 0.70 for the assessment method of breeders and from 0.00 to 0.63 for the assessment method of pathologists. For cycle one (C₁) population, where both foliar and root necrosis data were analyzed for clones assessed at clonal evaluation trials (CETs) and advanced yield trials (AYTs), H ² varied from 0.10 to 0.59 for the assessment method of breeders, while the H ² values ranged from 0.09 to 0.35 for the CBSD computation method of pathologists. In general, higher correlations were recorded for foliar severity from the assessment method of breeders (r = 0.4, p ≤ 0.01 for CBSD3s and r = 0.37, p ≤ 0.01 for CBSD6s) in C₁ clones evaluated at both clonal and advanced breeding stages than from the approach of pathologists. Ranking of top 10 C₁ clones by their indexed best linear unbiased predictors (BLUPs) for CBSD foliar and root necrosis showed four overlapping clones between clonal and advanced selection stages for the method of breeders; meanwhile, only a clone featured in both clonal and advanced selection stages from the CBSD assessment method of pathologists. Overall, the CBSD assessment method of breeders was more effective than the assessment method of pathologists, and thus, it justifies its continued use in CBSD resistance breeding.

Genetic Gains for Yield and Virus Disease Resistance of Cassava Varieties Developed Over the Last Eight Decades in Uganda

Achieving food security for an ever-increasing human population requires faster development of improved varieties. To this end, assessment of genetic gain for key traits is important to inform breeding processes. Despite the improvements made to increase production and productivity of cassava in Uganda at research level, there has been limited effort to quantify associated genetic gains. Accordingly, a study was conducted in Uganda to assess whether or not genetic improvement was evident in selected cassava traits using cassava varieties that were released from 1940 to 2019. Thirty-two varieties developed during this period, were evaluated simultaneously in three major cassava production zones; central (Namulonge), eastern (Serere), and northern (Loro). Best linear unbiased predictors (BLUPs) of the genotypic value for each clone were obtained across environments and regressed on order of release year to estimate annual genetic gains. We observed that genetic trends were mostly quadratic. On average, cassava mosaic disease (CMD) resistance increased by 1.9% per year, while annual genetic improvements in harvest index (0.0%) and fresh root yield (-5 kg per ha or -0.03% per ha) were non-substantial. For cassava brown streak disease (CBSD) resistance breeding which was only initiated in 2003, average annual genetic gains for CBSD foliar and CBSD root necrosis resistances were 2.3% and 1.5%, respectively. It's evident that cassava breeding has largely focused on protecting yield against diseases. This underpins the need for simultaneous improvement of cassava for disease resistance and high yield for the crop to meet its current and futuristic demands for food and industry.

Comprehensive genotyping of a Brazilian cassava (Manihot esculenta Crantz) germplasm bank: insights into diversification and domestication

KEY MESSAGE

Brazilian cassava diversity was characterized through population genetics and clustering approaches, highlighting contrasted genetic groups and spatial genetic differentiation. Cassava (Manihot esculenta Crantz) is a major staple root crop of the tropics, originating from the Amazonian region. In this study, 3354 cassava landraces and modern breeding lines from the Embrapa Cassava Germplasm Bank (CGB) were characterized. All individuals were subjected to genotyping-by-sequencing (GBS), identifying 27,045 single-nucleotide polymorphisms (SNPs). Identity-by-state and population structure analyses revealed a unique set of 1536 individuals and 10 distinct genetic groups with heterogeneous linkage disequilibrium (LD). On this basis, a density of 1300-4700 SNP markers were selected for large-effect quantitative trait loci (QTL) detection. Identified genetic groups were further characterized for population genetics parameters including minor allele frequency (MAF), observed heterozygosity [Formula: see text], effective population size estimate [Formula: see text]) and polymorphism information content (PIC). Selection footprints and introgressions of M. glaziovii were detected. Spatial population structure analysis revealed five ancestral populations related to distinct Brazilian ecoregions. Estimation of historical relationships among identified populations suggests an early population split from Amazonian to Atlantic forest and Caatinga ecoregions and active gene flows. This study provides a thorough genetic characterization of ex situ germplasm resources from cassava's center of origin, South America, with results shedding light on Brazilian cassava characteristics and its biogeographical landscape. These findings support and facilitate the use of genetic resources in modern breeding programs including implementation of association mapping and genomic selection strategies.

Comprehensive genotyping of a Brazilian cassava (Manihot esculenta Crantz) germplasm bank: insights into diversification and domestication

Brazilian cassava diversity was characterized through population genetics and clustering approaches, highlighting contrasted genetic groups and spatial genetic differentiation. Cassava (Manihot esculenta Crantz) is a major staple root crop of the tropics, originating from the Amazonian region. In this study, 3354 cassava landraces and modern breeding lines from the Embrapa Cassava Germplasm Bank (CGB) were characterized. All individuals were subjected to genotyping-by-sequencing (GBS), identifying 27,045 single-nucleotide polymorphisms (SNPs). Identity-by-state and population structure analyses revealed a unique set of 1536 individuals and 10 distinct genetic groups with heterogeneous linkage disequilibrium (LD). On this basis, a density of 1300-4700 SNP markers were selected for large-effect quantitative trait loci (QTL) detection. Identified genetic groups were further characterized for population genetics parameters including minor allele frequency (MAF), observed heterozygosity [Formula: see text], effective population size estimate [Formula: see text]) and polymorphism information content (PIC). Selection footprints and introgressions of M. glaziovii were detected. Spatial population structure analysis revealed five ancestral populations related to distinct Brazilian ecoregions. Estimation of historical relationships among identified populations suggests an early population split from Amazonian to Atlantic forest and Caatinga ecoregions and active gene flows. This study provides a thorough genetic characterization of ex situ germplasm resources from cassava's center of origin, South America, with results shedding light on Brazilian cassava characteristics and its biogeographical landscape. These findings support and facilitate the use of genetic resources in modern breeding programs including implementation of association mapping and genomic selection strategies.

Recent Large-Scale Genotyping and Phenotyping of Plant Genetic Resources of Vegetatively Propagated Crops

Several recent national and international projects have focused on large-scale genotyping of plant genetic resources in vegetatively propagated crops like fruit and berries, potatoes and woody ornamentals. The primary goal is usually to identify true-to-type plant material, detect possible synonyms, and investigate genetic diversity and relatedness among accessions. A secondary goal may be to create sustainable databases that can be utilized in research and breeding for several years ahead. Commonly applied DNA markers (like microsatellite DNA and SNPs) and next-generation sequencing each have their pros and cons for these purposes. Methods for large-scale phenotyping have lagged behind, which is unfortunate since many commercially important traits (yield, growth habit, storability, and disease resistance) are difficult to score. Nevertheless, the analysis of gene action and development of robust DNA markers depends on environmentally controlled screening of very large sets of plant material. Although more time-consuming, co-operative projects with broad-scale data collection are likely to produce more reliable results. In this review, we will describe some of the approaches taken in genotyping and/or phenotyping projects concerning a wide variety of vegetatively propagated crops.

Large-scale genome-wide association study, using historical data, identifies conserved genetic architecture of cyanogenic glucoside content in cassava (Manihot esculenta Crantz) root

Manihot esculenta (cassava) is a root crop originating from South America that is a major staple in the tropics, including in marginal environments. This study focused on South American and African germplasm and investigated the genetic architecture of hydrogen cyanide (HCN), a major component of root quality. HCN, representing total cyanogenic glucosides, is a plant defense component against herbivory but is also toxic for human consumption. We genotyped 3354 landraces and modern breeding lines originating from 26 Brazilian states and 1389 individuals were phenotypically characterized across multi-year trials for HCN. All plant material was subjected to high-density genotyping using genotyping by sequencing. We performed genome-wide association mapping to characterize the genetic architecture and gene mapping of HCN. Field experiments revealed strong broad- and narrow-sense trait heritability (0.82 and 0.41, respectively). Two major loci were identified, encoding for an ATPase and a MATE protein, and contributing up to 7 and 30% of the HCN concentration in roots, respectively. We developed diagnostic markers for breeding applications, validated trait architecture consistency in African germplasm and investigated further evidence for the domestication of sweet and bitter cassava. Fine genomic characterization revealed: (i) the major role played by vacuolar transporters in regulating HCN content; (ii) the co-domestication of sweet and bitter cassava major alleles are dependent upon geographical zone; and (iii) the major loci allele for high HCN in M. esculenta Crantz seems to originate from its ancestor, M. esculenta subsp. flabellifolia. Taken together, these findings expand our insights into cyanogenic glucosides in cassava roots and its glycosylated derivatives in plants.

Technological Innovations for Improving Cassava Production in Sub-Saharan Africa

Cassava is crucial for food security of millions of people in sub-Saharan Africa. The crop has great potential to contribute to African development and is increasing its income-earning potential for small-scale farmers and related value chains on the continent. Therefore, it is critical to increase cassava production, as well as its quality attributes. Technological innovations offer great potential to drive this envisioned change. This paper highlights genomic tools and resources available in cassava. The paper also provides a glimpse of how these resources have been used to screen and understand the pattern of cassava genetic diversity on the continent. Here, we reviewed the approaches currently used for phenotyping cassava traits, highlighting the methodologies used to link genotypic and phenotypic information, dissect the genetics architecture of key cassava traits, and identify quantitative trait loci/markers significantly associated with those traits. Additionally, we examined how knowledge acquired is utilized to contribute to crop improvement. We explored major approaches applied in the field of molecular breeding for cassava, their promises, and limitations. We also examined the role of national agricultural research systems as key partners for sustainable cassava production.

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.

Genetic analysis and QTL mapping for multiple biotic stress resistance in cassava

In Sub-Saharan Africa cassava (Manihot esculenta Crantz) is one of the most important food crops where more than 40% of the population relies on it as their staple carbohydrate source. Biotic constraints such as viral diseases, mainly Cassava Mosaic Disease (CMD) and Cassava Brown Streak Disease (CBSD), and arthropod pests, particularly Cassava Green Mite (CGM), are major constraints to the realization of cassava’s full production potential in Africa. To address these problems, we aimed to map the quantitative trait loci (QTL) associated with resistance to CBSD foliar and root necrosis symptoms, foliar CMD and CGM symptoms in a full-sib mapping population derived from the genotypes AR40-6 and Albert. A high-density linkage map was constructed with 2,125 SNP markers using a genotyping-by-sequencing approach. For phenotyping, clonal evaluation trials were conducted with 120 F₁ individuals for two consecutive field seasons using an alpha-lattice design at Chambezi and Naliendele, Tanzania. Previously identified QTL for resistance to CBSD foliar symptoms were corroborated, and a new putative QTL for CBSD root necrosis identified (qCBSDRNc14AR) from AR40-6. Two QTL were identified within the region of the previously recognized CMD2 locus from this population in which both parents are thought to possess the CMD2 locus. Interestingly, a minor but consistent QTL, qCGM18AR, for CGM resistance at 3 months after planting stage was also detected and co-localized with a previously identified SSR marker, NS346, linked with CGM resistance. Markers underlying these QTL may be used to increase efficiencies in cassava breeding programs.

Genome-Wide Association Analysis for Tuber Dry Matter and Oxidative Browning in Water Yam (Dioscorea alata L.)

Yam (Dioscorea spp.) is a nutritional and medicinal staple tuber crop grown in the tropics and sub-tropics. Among the food yam species, water yam (Dioscorea alata L.) is the most widely distributed and cultivated species worldwide. Tuber dry matter content (DMC) and oxidative browning (OxB) are important quality attributes that determine cultivar acceptability in water yam. This study used a single nucleotide polymorphism (SNP) assay from a diversity arrays technology (DArT) platform for a genome-wide association study (GWAS) of the two quality traits in a panel of 100 water yam clones grown in three environments. The marker–trait association analysis identified significant SNPs associated with tuber DMC on chromosomes 6 and 19 and with OxB on chromosome 5. The significant SNPs cumulatively explained 45.87 and 12.74% of the total phenotypic variation for the tuber DMC and OxB, respectively. Gene annotation for the significant SNP loci identified important genes associated in the process of the proteolytic modification of carbohydrates in the dry matter accumulation pathway as well as fatty acid β-oxidation in peroxisome for enzymatic oxidation. Additional putative genes were also identified in the peak SNP sites for both tuber dry matter and enzymatic oxidation with unknown functions. The results of this study provide valuable insight for further dissection of the genetic architecture of tuber dry matter and enzymatic oxidation in water yam. They also highlight SNP variants and genes useful for genomics-informed selection decisions in the breeding process for improving food quality traits in water yam.

Improving root characterisation for genomic prediction in cassava

Cassava is cultivated due to its drought tolerance and high carbohydrate-containing storage roots. The lack of uniformity and irregular shape of storage roots poses constraints on harvesting and post-harvest processing. Here, we phenotyped the Genetic gain and offspring (C1) populations from the International Institute of Tropical Agriculture (IITA) breeding program using image analysis of storage root photographs taken in the field. In the genome-wide association analysis (GWAS), we detected for most shape and size-related traits, QTL on chromosomes 1 and 12. In a previous study, we found the QTL on chromosome 12 to be associated with cassava mosaic disease (CMD) resistance. Because the root uniformity is important for breeding, we calculated the standard deviation (SD) of individual root measurements per clone. With SD measurements we identified new significant QTL for Perimeter, Feret and Aspect Ratio on chromosomes 6, 9 and 16. Predictive accuracies of root size and shape image-extracted traits were mostly higher than yield trait prediction accuracies. This study aimed to evaluate the feasibility of the image phenotyping protocol and assess GWAS and genomic prediction for size and shape image-extracted traits. The methodology described and the results are promising and open up the opportunity to apply high-throughput methods in cassava.

Mapping crown rust resistance at multiple time points in elite oat germplasm

Crown rust, caused by Puccinia coronata f. sp. avenae Erikss., is the most important disease impacting cultivated oat (Avena sativa L.). Genetic resistance is the most desirable management strategy. The genetic architecture of crown rust resistance is not fully understood, and previous mapping investigations have mostly ignored temporal variation. A collection of elite oat lines sourced from oat breeding programs in the American Upper Midwest and Canada was genotyped using a high-density genotyping-by-sequencing system and evaluated for crown rust disease severity at multiple time points throughout the growing season in three disease nursery environments. Genome-wide association mapping was conducted for disease severity on each observation date of each trial, area under the disease progress curve for each trial, heading date for each trial, and area under the disease progress curve in a multi-environment model. Crown rust resistance quantitative trait loci (QTL) were detected on linkage groups Mrg05, Mrg12, Mrg15, Mrg18, Mrg20, and Mrg33. None of these QTL were coincident with a days-to-heading QTL detected on Mrg02. Only the QTL detected on Mrg15 was detected in multiple mapping models. The QTL on Mrg05, Mrg12, Mrg18, Mrg20, and Mrg33 were detected on only a single observation date and were not detected on observations just days before and after. This result uncovers the importance of temporal variation in mapping experiments which is usually ignored. It is possible that high density temporal data could be used to more precisely characterize the nature of plant resistance in other systems.

Genome wide association study of 5 agronomic traits in olive (Olea europaea L.)

Olive (Olea europaea L.) is one of the most economically and historically important fruit crops worldwide. Genetic progress for valuable agronomic traits has been slow in olive despite its importance and benefits. Advances in next generation sequencing technologies provide inexpensive and highly reproducible genotyping approaches such as Genotyping by Sequencing, enabling genome wide association study (GWAS). Here we present the first comprehensive GWAS study on olive using GBS. A total of 183 accessions (FULL panel) were genotyped using GBS, 94 from the Turkish Olive GenBank Resource (TOGR panel) and 89 from the USDA-ARS National Clonal Germplasm Repository (NCGR panel) in the USA. After filtering low quality and redundant markers, GWAS was conducted using 24,977 SNPs in FULL, TOGR and NCGR panels. In total, 52 significant associations were detected for leaf length, fruit weight, stone weight and fruit flesh to pit ratio using the MLM_K. Significant GWAS hits were mapped to their positions and 19 candidate genes were identified within a 10-kb distance of the most significant SNP. Our findings provide a framework for the development of markers and identification of candidate genes that could be used in olive breeding programs.

Historical Introgressions from a Wild Relative of Modern Cassava Improved Important Traits and May Be Under Balancing Selection

Introgression of alleles from wild relatives has often been adaptive in plant breeding. However, the significance of historical hybridization events in modern breeding is often not clear. Cassava (Manihot esculenta) is among the most important staple foods in the world, sustaining hundreds of millions of people in the tropics, especially in sub-Saharan Africa. Widespread genotyping makes cassava a model for clonally propagated root and tuber crops in the developing world, and provides an opportunity to study the modern benefits and consequences of historical introgression. We detected large introgressed Manihot glaziovii genome-segments in a collection of 2742 modern cassava landraces and elite germplasm, the legacy of a 1930s era breeding to combat disease epidemics. African landraces and improved varieties were, on average, 3.8% (max 13.6%) introgressed. Introgressions accounted for a significant (mean 20%, max 56%) portion of the heritability of tested traits. M. glaziovii alleles on the distal 10 Mb of chr. 1 increased dry matter and root number. On chr. 4, introgressions in a 20 Mb region improved harvest index and brown streak disease tolerance. We observed the introgression frequency on chr. 1 double over three cycles of selection, and that later stage trials selectively excluded homozygotes from consideration as varieties. This indicates a heterozygous advantage of introgressions. However, we also found that maintaining large recombination-suppressed introgressions in the heterozygous state allowed the accumulation of deleterious mutations. We conclude that targeted recombination of introgressions would increase the efficiency of cassava breeding by allowing simultaneous fixation of beneficial alleles and purging of genetic load.

Cassava yield traits predicted by genomic selection methods

Genomic selection (GS) has been used to optimize genetic gains when phenotypic selection is considered costly and difficult to measure. The objective of this work was to evaluate the efficiency and consistency of GS prediction for cassava yield traits (Manihot esculenta Crantz) using different methods, taking into account the effect of population structure. BLUPs and deregressed BLUPs were obtained for 888 cassava accessions and evaluated for fresh root yield, dry root yield and dry matter content in roots in 21 trials conducted from 2011 to 2016. The deregressed BLUPs obtained for the accessions from a 48K single nucleotide polymorphism dataset were used for genomic predictions based on the BayesB, BLASSO, RR-BLUP, G-BLUP and RKHS methods. The accessions’ BLUPs were used in the validation step using four cross-validation strategies, taking into account population structure and different GS methods. Similar estimates of predictive ability and bias were identified for the different genomic selection methods in the first cross-validation strategy. Lower predictive ability was observed for fresh root yield (0.4569 –RR-BLUP to 0.4756—RKHS) and dry root yield (0.4689 –G-BLUP to 0.4818—RKHS) in comparison with dry matter content (0.5655 –BLASSO to 0.5670 –RKHS). However, the RKHS method exhibited higher efficiency and consistency in most of the validation scenarios in terms of prediction ability for fresh root yield and dry root yield. The correlations of the genomic estimated breeding values between the genomic selection methods were quite high (0.99–1.00), resulting in high coincidence of clone selection regardless of the genomic selection method. The deviance analyses within and between the validation clusters formed by the discriminant analysis of principal components were significant for all traits. Therefore, this study indicated that i) the prediction of dry matter content was more accurate compared to that of yield traits, possibly as a result of the smaller influence of non-additive genetic effects; ii) the RKHS method resulted in high and stable prediction ability in most of the validation scenarios; and iii) some kinship between the validation and training populations is desirable in order for genomic selection to succeed due to the significant effect of population structure on genomic selection predictions.

Engineering Disease-Resistant Cassava

Manihot esculenta Crantz (cassava) is a food crop originating from South America grown primarily for its starchy storage roots. Today, cassava is grown in the tropics of South America, Africa, and Asia with an estimated 800 million people relying on it as a staple source of calories. In parts of sub-Saharan Africa, cassava is particularly crucial for food security. Cassava root starch also has use in the pharmaceutical, textile, paper, and biofuel industries. Cassava has seen strong demand since 2000 and production has increased consistently year-over-year, but potential yields are hampered by susceptibility to biotic and abiotic stresses. In particular, bacterial and viral diseases can cause severe yield losses. Of note are cassava bacterial blight (CBB), cassava mosaic disease (CMD), and cassava brown streak disease (CBSD), all of which can cause catastrophic losses for growers. In this article, we provide an overview of the major microbial diseases of cassava, discuss current and potential future efforts to engineer new sources of resistance, and conclude with a discussion of the regulatory hurdles that face biotechnology.

RNA silencing of South African cassava mosaic virus in transgenic cassava expressing AC1/AC4 hp- RNA induces tolerance

•

A hairpin-RNA construct targeting South African cassava mosaic virus AC1/AC4 overlapping genome region confers tolerance in cassava.

•

Viral load and symptom severity were significantly lower in transgenic lines transformed with AC1/AC4 construct compared to wildtype cv.60444.

•

SACMV transgenic lines infected with SACMV displayed a recovery phenotype similar to cassava mosaic disease-tolerant cassava landrace TME3.

Cassava mosaic disease (CMD), caused by geminiviruses, is a major hurdle to cassava production. Due to the heterozygous nature of cassava, breeding for virus resistance is difficult, but cassava has been shown to be a good candidate for genetic engineering using RNA interference (RNAi). T This study reports on the ability of a transgene-derived RNA hairpin, homologous to an overlapping region of the SACMV replication associated protein and putative virus suppressor of silencing protein (AC1/AC4), to confer tolerance in the CMD-susceptible model cassava cultivar 60444. Three of the fourteen transgenic lines expressing SACMV AC1/AC4 hairpin-derived siRNAs showed decreased symptoms and viral loads compared to untransformed control plants. Expression of SACMV AC1/AC4 homologous siRNAs showed that this tolerance is most likely associated with post-transcriptional gene silencing of the virus. This is the first report of targeting the overlapping AC1 and AC4 genes of SACMV conferring CMD tolerance in cassava.

Haplotype-resolved genomes of geminivirus-resistant and geminivirus-susceptible African cassava cultivars

Background

Cassava is an important food crop in tropical and sub-tropical regions worldwide. In Africa, cassava production is widely affected by cassava mosaic disease (CMD), which is caused by the African cassava mosaic geminivirus that is transmitted by whiteflies. Cassava breeders often use a single locus, CMD2, for introducing CMD resistance into susceptible cultivars. The CMD2 locus has been genetically mapped to a 10-Mbp region, but its organization and genes as well as their functions are unknown.

Results

We report haplotype-resolved de novo assemblies and annotations of the genomes for the African cassava cultivar TME (tropical Manihot esculenta), which is the origin of CMD2, and the CMD-susceptible cultivar 60444. The assemblies provide phased haplotype information for over 80% of the genomes. Haplotype comparison identified novel features previously hidden in collapsed and fragmented cassava genomes, including thousands of allelic variants, inter-haplotype diversity in coding regions, and patterns of diversification through allele-specific expression. Reconstruction of the CMD2 locus revealed a highly complex region with nearly identical gene sets but limited microsynteny between the two cultivars.

Conclusions

The genome maps of the CMD2 locus in both 60444 and TME3, together with the newly annotated genes, will help the identification of the causal genetic basis of CMD2 resistance to geminiviruses. Our de novo cassava genome assemblies will also facilitate genetic mapping approaches to narrow the large CMD2 region to a few candidate genes for better informed strategies to develop robust geminivirus resistance in susceptible cassava cultivars.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0697-6) contains supplementary material, which is available to authorized users.

Genetic Variation and Trait Correlations in an East African Cassava Breeding Population for Genomic Selection

Cassava (Manihot esculenta Crantz) is a major source of dietary carbohydrates for >700 million people globally. However, its long breeding cycle has slowed the rate of genetic gain for target traits. This study aimed to asses genetic variation, the level of inbreeding, and trait correlations in genomic selection breeding cycles. We used phenotypic and genotypic data from the National Crops Resources Research Institute (NaCRRI) foundation population (Cycle 0, C₀) and the progeny (Cycle 1, C₁) derived from crosses of 100 selected C0 clones as progenitors, both to evaluate and optimize genomic selection. The highest broad-sense heritability (H ² = 0.95) and narrow-sense heritability (h ² = 0.81) were recorded for cassava mosaic disease severity and the lowest for root weight per plot (H ² = 0.06 and h ² = 0.00). We observed the highest genetic correlation (r _g= 0.80) between cassava brown streak disease root incidence measured at seedling and clonal stages of evaluation, suggesting the usefulness of seedling data in predicting clonal performance for cassava brown streak root necrosis. Similarly, high genetic correlations were observed between cassava brown streak disease severity (r _g= 0.83) scored at 3 and 6 mo after planting (MAP) and cassava mosaic disease, scored at 3 and 6 MAP (r _g= 0.95), indicating that data obtained on these two diseases at 6 MAP would suffice. Population differentiation between C₀ and C₁ was not well defined, implying that the 100 selected progenitors of C₁ captured the diversity in the C₀. Overall, genetic gain for most traits were observed from C₀ to C₁.