Regional Heritability Mapping Provides Insights into Dry Matter Content in African White and Yellow Cassava Populations

New protocol for rapid cassava multiplication in field conditions: a perspective on speed breeding

Despite the economic and social importance, high-yielding cassava cultivars are only released after extensive research, mainly due to the low multiplication rate. This study aimed to assess the impact of using smaller-sized seed cuttings treated with agrochemicals (8MP) compared to the conventional planting size (16 cm) on genetic parameters, agronomic performance, and the ranking of cassava clones based on yield and growth attributes. The evaluation was carried out in clonal evaluation trial (CET), preliminary yield trial (PYT), and uniform yield trials (UYT). Additionally, a new selection scheme for cassava breeding programs was proposed. A total of 169 clones were evaluated, including 154 improved clones at different stages of selection and 15 local varieties used as checks. Field trials were conducted using both sizes of propagative material (8MP and 16 cm) in each phase of the breeding program. The data were analyzed using mixed models, considering the random effects of genotype and genotype-environment interaction (G×E) to determine variances and heritabilities. Bland-Altman concordance and correlation analysis of selection indices were employed to examine the consistency in the ranking of cassava clones using different seed cutting sizes. The distribution of variance components, heritabilities, means, and range of the 8MP and 16 cm trials in different phases of the cassava breeding program exhibited remarkable similarity, thereby enabling a comparative assessment of similar genetic effects. With a selection intensity of 30%, the concordance in clone ranking was 0.41, 0.57, and 0.85 in CET, PYT, and UYT trials, respectively, when comparing the selection based on 8MP and 16 cm trials. It is worth noting that the ranking of the top 15% remained largely unchanged. Based on the findings, proposed changes in the cassava selection scheme involve increasing the number of trials starting from the CET phase, early incorporation of G×E interaction, elimination of the PYT trial, reduction of the breeding cycle from 5 to 3 years, and a decrease in the time required for variety development from 11 to 9 years. These modifications are expected to lead to cost reduction and enhance the effectiveness of cassava breeding programs.

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.

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.

Regional Heritability Mapping of Quantitative Trait Loci Controlling Traits Related to Growth and Productivity in Popcorn (Zea mays L.)

The method of regional heritability mapping (RHM) has become an important tool in the identification of quantitative trait loci (QTLs) controlling traits of interest in plants. Here, RHM was first applied in a breeding population of popcorn, to identify the QTLs and candidate genes involved in grain yield, plant height, kernel popping expansion, and first ear height, as well as determining the heritability of each significant genomic region. The study population consisted of 98 S1 families derived from the 9th recurrent selection cycle (C-9) of the open-pollinated variety UENF-14, which were genetically evaluated in two environments (ENV1 and ENV2). Seventeen and five genomic regions were mapped by the RHM method in ENV1 and ENV2, respectively. Subsequent genome-wide analysis based on the reference genome B73 revealed associations with forty-six candidate genes within these genomic regions, some of them are considered to be biologically important due to the proteins that they encode. The results obtained by the RHM method have the potential to contribute to knowledge on the genetic architecture of the growth and yield traits of popcorn, which might be used for marker-assisted selection in breeding programs.

A Low Resolution Epistasis Mapping Approach To Identify Chromosome Arm Interactions in Allohexaploid Wheat

Epistasis is an important contributor to genetic variance. In inbred populations, pairwise epistasis is present as additive by additive interactions. Testing for epistasis presents a multiple testing problem as the pairwise search space for modest numbers of markers is large. Single markers do not necessarily track functional units of interacting chromatin as well as haplotype based methods do. To harness the power of multiple markers while minimizing the number of tests conducted, we present a low resolution test for epistatic interactions across whole chromosome arms. Epistasis covariance matrices were constructed from the additive covariances of individual chromosome arms. These covariances were subsequently used to estimate an epistatic variance parameter while correcting for background additive and epistatic effects. We find significant epistasis for 2% of the interactions tested for four agronomic traits in a winter wheat breeding population. Interactions across homeologous chromosome arms were identified, but were less abundant than other chromosome arm pair interactions. The homeologous chromosome arm pair 4BL/4DL showed a strong negative relationship between additive and interaction effects that may be indicative of functional redundancy. Several chromosome arms appeared to act as hubs in an interaction network, suggesting that they may contain important regulatory factors. The differential patterns of epistasis across different traits demonstrate that detection of epistatic interactions is robust when correcting for background additive and epistatic effects in the population. The low resolution epistasis mapping method presented here identifies important epistatic interactions with a limited number of statistical tests at the cost of low precision.

Genome-Wide Association and Regional Heritability Mapping of Plant Architecture, Lodging and Productivity in Phaseolus vulgaris

The availability of high-density molecular markers in common bean has allowed to explore the genetic basis of important complex agronomic traits with increased resolution. Genome-Wide Association Studies (GWAS) and Regional Heritability Mapping (RHM) are two analytical approaches for the detection of genetic variants. We carried out GWAS and RHM for plant architecture, lodging and productivity across two important growing environments in Brazil in a germplasm of 188 common bean varieties using DArTseq genotyping strategies. The coefficient of determination of G × E interaction (c2int ) was equal to 17, 21 and 41%, respectively for the traits architecture, lodging, and productivity. Trait heritabilities were estimated at 0.81 (architecture), 0.79 (lodging) and 0.43 (productivity), and total genomic heritability accounted for large proportions (72% to ≈100%) of trait heritability. At the same probability threshold, three marker-trait associations were detected using GWAS, while RHM detected eight QTL encompassing 145 markers along five chromosomes. The proportion of genomic heritability explained by RHM was considerably higher (35.48 to 58.02) than that explained by GWAS (28.39 to 30.37). In general, RHM accounted for larger fractions of the additive genetic variance being captured by markers effects inside the defined regions. Nevertheless, a considerable proportion of the heritability is still missing (∼42% to ∼64%), probably due to LD between markers and genes and/or rare allele variants not sampled. RHM in autogamous species had the potential to identify larger-effect QTL combining allelic variants that could be effectively incorporated into whole-genome prediction models and tracked through breeding generations using marker-assisted selection.