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

Validation of SNP Markers for Diversity Analysis, Quality Control, and Trait Selection in a Biofortified Cassava Population

A validated marker system is crucial to running an effective genomics-assisted breeding program. We used 36 Kompetitive Allele-Specific PCR (KASP) markers to genotype 376 clones from the biofortified cassava pipeline, and fingerprinted 93 of these clones with DArTseq markers to characterize breeding materials and evaluate their relationships. The discriminating ability of the 36-quality control (QC) KASP and 6602 DArTseq markers was assessed using 92 clones genotyped in both assays. In addition, trait-specific markers were used to determine the presence or absence of target genomic regions. Hierarchical clustering identified two major groups, and the clusters were consistent with the breeding program origins. There was moderate genetic differentiation and a low degree of variation between the identified groups. The general structure of the population was similar using both assays. Nevertheless, KASP markers had poor resolution when it came to differentiating the genotypes by seed sources and overestimated the prevalence of duplicates. The trait-linked markers did not achieve optimal performance as all markers displayed variable levels of false positive and/or false negative. These findings represent the initial step in the application of genomics-assisted breeding for the biofortified cassava pipeline, and will guide the use of genomic selection in the future.

Genetic diversity and population structure of Uganda cassava germplasm

Cassava (Manihot esculenta Crantz) holds significant economic importance globally. Evaluating a diverse range of germplasm based on molecular characteristics not only enhances its preservation but also supports its utilization in breeding programs. In this study, we assessed genetic diversity and population structure among 155 cassava genotypes from Uganda using 5247 single nucleotide polymorphism (SNP) markers. Genotyping by sequencing (GBS) was employed for SNP discovery and to evaluate genetic diversity and population structure using the ADMIXTURE software. The cassava accessions comprised two populations: 49 accessions from Ugandan lines and 106 accessions resulting from crosses between South American and Ugandan lines. The average call rate of 96% was utilized to assess marker polymorphism. Polymorphic information content values of the markers ranged from 0.1 to 0.5 with an average of 0.4 which was moderately high. The principal component analysis (PCA) showed that the first two components captured ~ 24.2% of the genetic variation. The average genetic diversity was 0.3. The analysis of molecular variance (AMOVA) indicated that 66.02% and 33.98% of the total genetic variation occurred within accessions and between sub-populations, respectively. Five sub-populations were identified based on ADMIXTURE structure analysis (K = 5). Neighbor-joining tree and hierarchical clustering tree revealed the presence of three different groups which were primarily based on the source of the genotypes. The results suggested that there was considerable genetic variation among the cassava genotypes which is useful in cassava improvement and conservation efforts.

Cassava Breeding and Cultivation Challenges in Thailand: Past, Present, and Future Perspectives

Cassava (Manihot esculenta Crantz) was introduced to Southeast Asia in the 16th-17th centuries and has since flourished as an industrial crop. Since the 1980s, Thailand has emerged as the leading producer and exporter of cassava products. This growth coincided with the initiation of cassava breeding programs in collaboration with the International Center for Tropical Agriculture (CIAT), focusing on root yield and starch production. The success of Thai cassava breeding programs can be attributed to the incorporation of valuable genetic diversity from international germplasm resources to cross with the local landraces, which has become the genetic foundation of many Thai commercial varieties. Effective evaluation under diverse environmental conditions has led to the release of varieties with high yield stability. A notable success is the development of Kasetsart 50. However, extreme climate change poses significant challenges, including abiotic and biotic stresses that threaten cassava root yield and starch content, leading to a potential decline in starch-based industries. Future directions for cassava breeding must include hybrid development, marker-assisted recurrent breeding, and gene editing, along with high-throughput phenotyping and flower induction. These strategies are essential to achieve breeding objectives focused on drought tolerance and disease resistance, especially for CMD and CBSD.

Phenotypic Variability in Resistance to Anthracnose, White, Brown, and Blight Leaf Spot in Cassava Germplasm

Despite fungal diseases affecting the aerial parts of cassava (Manihot esculenta Crantz) and causing significant yield losses, there is a lack of comprehensive studies assessing resistance in the species' germplasm. This study aimed to evaluate the phenotypic diversity for resistance to anthracnose disease (CAD), blight leaf spot (BliLS), brown leaf spot (BLS), and white leaf spot (WLS) in cassava germplasm and to identify genotypes suitable for breeding purposes. A total of 837 genotypes were evaluated under field conditions across two production cycles (2021 and 2022). Artificial inoculations were carried out in the field, and data on yield and disease severity were collected using a standardized rating scale. The top 25 cassava genotypes were selected based on a selection index for disease resistance and agronomic traits. High environmental variability resulted in low heritabilities (h2) for CAD, WLS, and BLS (h2 = 0.42, 0.34, 0.29, respectively) and moderate heritability for BliLS (h2 = 0.51). While the range of data for disease resistance was narrow, it was considerably wider for yield traits. Cluster analysis revealed that increased yield traits and disease severity were associated with higher scores of the first and second discriminant functions, respectively. Thus, most clusters comprised genotypes with hybrid characteristics for both traits. Overall, there was a strong correlation among aerial diseases, particularly between BLS and BliLS (r = 0.96), while the correlation between CAD and other diseases ranged from r = 0.53 to 0.58. Yield traits showed no significant correlations with disease resistance. Although the mean selection differential for disease resistance was modest (between -2.31% and -3.61%), selection based on yield traits showed promising results, particularly for fresh root yield (82%), dry root yield (39%), shoot yield (49%), and plant vigor (26%). This study contributes to enhancing genetic gains for resistance to major aerial part diseases and improving yield traits in cassava breeding programs.

Exploiting DNA methylation in cassava under water deficit for crop improvement

DNA methylation plays a key role in the development and plant responses to biotic and abiotic stresses. This work aimed to evaluate the DNA methylation in contrasting cassava genotypes for water deficit tolerance. The varieties BRS Formosa (bitter) and BRS Dourada (sweet) were grown under greenhouse conditions for 50 days, and afterwards, irrigation was suspended. The stressed (water deficit) and non-stressed plants (negative control) consisted the treatments with five plants per variety. The DNA samples of each variety and treatment provided 12 MethylRAD-Seq libraries (two cassava varieties, two treatments, and three replicates). The sequenced data revealed methylated sites covering 18 to 21% of the Manihot esculenta Crantz genome, depending on the variety and the treatment. The CCGG methylated sites mapped mostly in intergenic regions, exons, and introns, while the CCNGG sites mapped mostly intergenic, upstream, introns, and exons regions. In both cases, methylated sites in UTRs were less detected. The differentially methylated sites analysis indicated distinct methylation profiles since only 12% of the sites (CCGG and CCNGG) were methylated in both varieties. Enriched gene ontology terms highlighted the immediate response of the bitter variety to stress, while the sweet variety appears to suffer more potential stress-damages. The predicted protein-protein interaction networks reinforced such profiles. Additionally, the genomes of the BRS varieties uncovered SNPs/INDELs events covering genes stood out by the interactomes. Our data can be useful in deciphering the roles of DNA methylation in cassava drought-tolerance responses and adaptation to abiotic stresses.

Individual-based landscape genomics for conservation: An analysis pipeline

Landscape genomics can harness environmental and genetic data to inform conservation decisions by providing essential insights into how landscapes shape biodiversity. The massive increase in genetic data afforded by the genomic era provides exceptional resolution for answering critical conservation genetics questions. The accessibility of genomic data for non-model systems has also enabled a shift away from population-based sampling to individual-based sampling, which now provides accurate and robust estimates of genetic variation that can be used to examine the spatial structure of genomic diversity, population connectivity and the nature of environmental adaptation. Nevertheless, the adoption of individual-based sampling in conservation genetics has been slowed due, in large part, to concerns over how to apply methods developed for population-based sampling to individual-based sampling schemes. Here, we discuss the benefits of individual-based sampling for conservation and describe how landscape genomic methods, paired with individual-based sampling, can answer fundamental conservation questions. We have curated key landscape genomic methods into a user-friendly, open-source workflow, which we provide as a new R package, A Landscape Genomics Analysis Toolkit in R (algatr). The algatr package includes novel added functionality for all of the included methods and extensive vignettes designed with the primary goal of making landscape genomic approaches more accessible and explicitly applicable to conservation biology.

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.

Selective signatures and high genome-wide diversity in traditional Brazilian manioc (Manihot esculenta Crantz) varieties

Knowledge about genetic diversity is essential to promote effective use and conservation of crops, because it enables farmers to adapt their crops to specific needs and is the raw material for breeding. Manioc (Manihot esculenta ssp. esculenta) is one of the world's major food crops and has the potential to help achieve food security in the context of on-going climate changes. We evaluated single nucleotide polymorphisms in traditional Brazilian manioc varieties conserved in the gene bank of the Luiz de Queiroz College of Agriculture, University of São Paulo. We assessed genome-wide diversity and identified selective signatures contrasting varieties from different biomes with samples of manioc's wild ancestor M. esculenta ssp. flabellifolia. We identified signatures of selection putatively associated with resistance genes, plant development and response to abiotic stresses that might have been important for the crop's domestication and diversification resulting from cultivation in different environments. Additionally, high neutral genetic diversity within groups of varieties from different biomes and low genetic divergence among biomes reflect the complexity of manioc's evolutionary dynamics under traditional cultivation. Our results exemplify how smallholder practices contribute to conserve manioc's genetic resources, maintaining variation of potential adaptive significance and high levels of neutral genetic diversity.

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.

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.

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 (FST) 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.