Genetic analysis of cassava brown streak disease root necrosis using image analysis and genome-wide association studies.

Cassava brown streak disease (CBSD) poses a substantial threat to food security. To address this challenge, we used PlantCV to extract CBSD root necrosis image traits from 320 clones, with an aim of identifying genomic regions through genome-wide association studies (GWAS) and candidate genes. Results revealed strong correlations among certain root necrosis image traits, such as necrotic area fraction and necrotic width fraction, as well as between the convex hull area of root necrosis and the percentage of necrosis. Low correlations were observed between CBSD scores obtained from the 1-5 scoring method and all root necrosis traits. Broad-sense heritability estimates of root necrosis image traits ranged from low to moderate, with the highest estimate of 0.42 observed for the percentage of necrosis, while narrow-sense heritability consistently remained low, ranging from 0.03 to 0.22. Leveraging data from 30,750 SNPs obtained through DArT genotyping, eight SNPs on chromosomes 1, 7, and 11 were identified and associated with both the ellipse eccentricity of root necrosis and the percentage of necrosis through GWAS. Candidate gene analysis in the 172.2kb region on the chromosome 1 revealed 24 potential genes with diverse functions, including ubiquitin-protein ligase, DNA-binding transcription factors, and RNA metabolism protein, among others. Despite our initial expectation that image analysis objectivity would yield better heritability estimates and stronger genomic associations than the 1-5 scoring method, the results were unexpectedly lower. Further research is needed to comprehensively understand the genetic basis of these traits and their relevance to cassava breeding and disease management.

Genetic complexity of cassava brown streak disease: insights from qPCR-based viral titer analysis and genome-wide association studies.

Cassava, a vital global food source, faces a threat from Cassava Brown Streak Disease (CBSD). CBSD results from two viruses: Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV). These viruses frequently pose challenges to the traditional symptom-based 1-5 phenotyping method due to its limitations in terms of accuracy and objectivity. Quantitative polymerase chain reaction (qPCR) offers precise virus quantification, although high costs hinder its widespread adoption. In this research, we utilized qPCR to measure the viral titer/load of CBSV and UCBSV. The objectives were to evaluate titer variability within the Cycle 2 (C2) population in two different environments, establish connections between viral titers and CBSD severity scores from the 1-5 scoring method, perform Genome-Wide Association Studies (GWAS) to identify genomic regions associated with CBSV and UCBSV titers, and investigate the functional annotated genes. The results demonstrated a significantly higher prevalence of CBSV (50.2%) in clones compared to UCBSV (12.9%) with mixed infections in some cases. Genotypic effects, particularly concerning UCBSV, were significant, with genotype-by-environment effects primarily influencing CBSV titer. GWAS Studies identified genomic regions associated with CBSV and UCBSV titers. Twenty-one SNP markers on chromosomes 10, 13, 17, and 18 exhibited significant associations with CBSV titer, collectively explaining 43.14% of the phenotypic variation. Additionally, 25 SNP markers on chromosomes 1, 2, 4, 5, 8, 11, 12, 13, 16, and 18 were associated with UCBSV titer, and explained 70.71% of the phenotypic variation. No shared genomic regions were identified between CBSV and UCBSV viral titers. Gene ontology analysis also revealed diverse gene functions, especially in transport and catalytic activities. These findings enhance our understanding of virus prevalence, genetics, and molecular functions in cassava plants, offering valuable insights for targeted breeding strategies.

Genomic prediction of seed nutritional traits in biparental families of oat (Avena sativa).

Selection for more nutritious crop plants is an important goal of plant breeding to improve food quality and contribute to human health outcomes. While there are efforts to integrate genomic prediction to accelerate breeding progress, an ongoing challenge is identifying strategies to improve accuracy when predicting within biparental populations in breeding programs. We tested multiple genomic prediction methods for 12 seed fatty acid content traits in oat (Avena sativa L.), as unsaturated fatty acids are a key nutritional trait in oat. Using two well-characterized oat germplasm panels and other biparental families as training populations, we predicted family mean and individual values within families. Genomic prediction of family mean exceeded a mean accuracy of 0.40 and 0.80 using an unrelated and related germplasm panel, respectively, where the related germplasm panel outperformed prediction based on phenotypic means (0.54). Within family prediction accuracy was more variable: training on the related germplasm had higher accuracy than the unrelated panel (0.14-0.16 and 0.05-0.07, respectively), but variability between families was not easily predicted by parent relatedness, segregation of a locus detected by a genome-wide association study in the panel, or other characteristics. When using other families as training populations, prediction accuracies were comparable to the related germplasm panel (0.11-0.23), and families that had half-sib families in the training set had higher prediction accuracy than those that did not. Overall, this work provides an example of genomic prediction of family means and within biparental families for an important nutritional trait and suggests that using related germplasm panels as training populations can be effective.

If it ain't broke, don't fix it: evaluating the effect of increased recombination on response to selection for wheat breeding.

Meiotic recombination is a source of allelic diversity, but the low frequency and biased distribution of crossovers that occur during meiosis limits the genetic variation available to plant breeders. Simulation studies previously identified that increased recombination frequency can retain more genetic variation and drive greater genetic gains than wildtype recombination. Our study was motivated by the need to define desirable recombination intervals in regions of the genome with fewer crossovers. We hypothesized that deleterious variants, which can negatively impact phenotypes and occur at higher frequencies in low recombining regions where they are linked in repulsion with favorable loci, may offer a signal for positioning shifts of recombination distributions. Genomic selection breeding simulation models based on empirical wheat data were developed to evaluate increased recombination frequency and changing recombination distribution on response to selection. Comparing high and low values for a range of simulation parameters identified that few combinations retained greater genetic variation and fewer still achieved higher genetic gain than wildtype. More recombination was associated with loss of genomic prediction accuracy, which outweighed the benefits of disrupting repulsion linkages. Irrespective of recombination frequency or distribution and deleterious variant annotation, enhanced response to selection under increased recombination required polygenic trait architecture, high heritability, an initial scenario of more repulsion than coupling linkages, and greater than 6 cycles of genomic selection. Altogether, the outcomes of this research discourage a controlled recombination approach to genomic selection in wheat as a more efficient path to retaining genetic variation and increasing genetic gains compared with existing breeding methods.

Parsimonious genotype by environment interaction covariance models for cassava (Manihot esculenta).

The assessment of cassava clones across multiple environments is often carried out at the uniform yield trial, a late evaluation stage, before variety release. This is to assess the differential response of the varieties across the testing environments, a phenomenon referred to as genotype-by-environment interaction (GEI). This phenomenon is considered a critical challenge confronted by plant breeders in developing crop varieties. This study used the data from variety trials established as randomized complete block design (RCBD) in three replicates across 11 locations in different agro-ecological zones in Nigeria over four cropping seasons (2016-2017, 2017-2018, 2018-2019, and 2019-2020). We evaluated a total of 96 varieties, including five checks, across 48 trials. We exploited the intricate pattern of GEI by fitting variance-covariance structure models on fresh root yield. The goodness-of-fit statistics revealed that the factor analytic model of order 3 (FA3) is the most parsimonious model based on Akaike Information Criterion (AIC). The three-factor loadings from the FA3 model explained, on average across the 27 environments, 53.5% [FA (1)], 14.0% [FA (2)], and 11.5% [FA (3)] of the genetic effect, and altogether accounted for 79.0% of total genetic variability. The association of factor loadings with weather covariates using partial least squares regression (PLSR) revealed that minimum temperature, precipitation and relative humidity are weather conditions influencing the genotypic response across the testing environments in the southern region and maximum temperature, wind speed, and temperature range for those in the northern region of Nigeria. We conclude that the FA3 model identified the common latent factors to dissect and account for complex interaction in multi-environment field trials, and the PLSR is an effective approach for describing GEI variability in the context of multi-environment trials where external environmental covariables are included in modeling.

Skinny kelp (Saccharina angustissima) provides valuable genetics for the biomass improvement of farmed sugar kelp (Saccharina latissima).

Saccharina latissima (sugar kelp) is one of the most widely cultivated brown marine macroalgae species in the North Atlantic and the eastern North Pacific Oceans. To meet the expanding demands of the sugar kelp mariculture industry, selecting and breeding sugar kelp that is best suited to offshore farm environments is becoming necessary. To that end, a multi-year, multi-institutional breeding program was established by the U.S. Department of Energy's (DOE) Advanced Research Projects Agency-Energy (ARPA-E) Macroalgae Research Inspiring Novel Energy Resources (MARINER) program. Hybrid sporophytes were generated using 203 unique gametophyte cultures derived from wild-collected Saccharina spp. for two seasons of farm trials (2019-2020 and 2020-2021). The wild sporophytes were collected from 10 different locations within the Gulf of Maine (USA) region, including both sugar kelp (Saccharina latissima) and the skinny kelp species (Saccharina angustissima). We harvested 232 common farm plots during these two seasons with available data. We found that farmed kelp plots with skinny kelp as parents had an average increased yield over the mean (wet weight 2.48 ± 0.90 kg m-1 and dry weight 0.32 ± 0.10 kg m-1) in both growing seasons. We also found that blade length positively correlated with biomass in skinny kelp x sugar kelp crosses or pure sugar kelp crosses. The skinny x sugar progenies had significantly longer and narrower blades than the pure sugar kelp progenies in both seasons. Overall, these findings suggest that sugar x skinny kelp crosses provide improved yield compared to pure sugar kelp crosses. Supplementary Information: The online version contains supplementary material available at 10.1007/s10811-022-02811-1.

Exploring genotype by environment interaction on cassava yield and yield related traits using classical statistical methods.

Variety advancement decisions for root quality and yield-related traits in cassava are complex due to the variable patterns of genotype-by-environment interactions (GEI). Therefore, studies focused on the dissection of the existing patterns of GEI using linear-bilinear models such as Finlay-Wilkinson (FW), additive main effect and multiplicative interaction (AMMI), and genotype and genotype-by-environment (GGE) interaction models are critical in defining the target population of environments (TPEs) for future testing, selection, and advancement. This study assessed 36 elite cassava clones in 11 locations over three cropping seasons in the cassava breeding program of IITA based in Nigeria to quantify the GEI effects for root quality and yield-related traits. Genetic correlation coefficients and heritability estimates among environments found mostly intermediate to high values indicating high correlations with the major TPE. There was a differential clonal ranking among the environments indicating the existence of GEI as also revealed by the likelihood ratio test (LRT), which further confirmed the statistical model with the heterogeneity of error variances across the environments fit better. For all fitted models, we found the main effects of environment, genotype, and interaction significant for all observed traits except for dry matter content whose GEI sensitivity was marginally significant as found using the FW model. We identified TMS14F1297P0019 and TMEB419 as two topmost stable clones with a sensitivity values of 0.63 and 0.66 respectively using the FW model. However, GGE and AMMI stability value in conjunction with genotype selection index revealed that IITA-TMS-IBA000070 and TMS14F1036P0007 were the top-ranking clones combining both stability and yield performance measures. The AMMI-2 model clustered the testing environments into 6 mega-environments based on winning genotypes for fresh root yield. Alternatively, we identified 3 clusters of testing environments based on genotypic BLUPs derived from the random GEI component.

GrainGenes: a data-rich repository for small grains genetics and genomics.

As one of the US Department of Agriculture-Agricultural Research Service flagship databases, GrainGenes (https://wheat.pw.usda.gov) serves the data and community needs of globally distributed small grains researchers for the genetic improvement of the Triticeae family and Avena species that include wheat, barley, rye and oat. GrainGenes accomplishes its mission by continually enriching its cross-linked data content following the findable, accessible, interoperable and reusable principles, enhancing and maintaining an intuitive web interface, creating tools to enable easy data access and establishing data connections within and between GrainGenes and other biological databases to facilitate knowledge discovery. GrainGenes operates within the biological database community, collaborates with curators and genome sequencing groups and contributes to the AgBioData Consortium and the International Wheat Initiative through the Wheat Information System (WheatIS). Interactive and linked content is paramount for successful biological databases and GrainGenes now has 2917 manually curated gene records, including 289 genes and 254 alleles from the Wheat Gene Catalogue (WGC). There are >4.8 million gene models in 51 genome browser assemblies, 6273 quantitative trait loci and >1.4 million genetic loci on 4756 genetic and physical maps contained within 443 mapping sets, complete with standardized metadata. Most notably, 50 new genome browsers that include outputs from the Wheat and Barley PanGenome projects have been created. We provide an example of an expression quantitative trait loci track on the International Wheat Genome Sequencing Consortium Chinese Spring wheat browser to demonstrate how genome browser tracks can be adapted for different data types. To help users benefit more from its data, GrainGenes created four tutorials available on YouTube. GrainGenes is executing its vision of service by continuously responding to the needs of the global small grains community by creating a centralized, long-term, interconnected data repository. Database URL:https://wheat.pw.usda.gov.

Generalizable approaches for genomic prediction of metabolites in plants.

Plant metabolites are important traits for plant breeders seeking to improve nutrition and agronomic performance yet integrating selection for metabolomic traits can be limited by phenotyping expense and degree of genetic characterization, especially of uncommon metabolites. As such, developing generalizable genomic selection methods based on biochemical pathway biology for metabolites that are transferable across plant populations would benefit plant breeding programs. We tested genomic prediction accuracy for >600 metabolites measured by gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS) in oat (Avena sativa L.) seed. Using a discovery germplasm panel, we conducted metabolite genome-wide association study (mGWAS) and selected loci to use in multikernel models that encompassed metabolome-wide mGWAS results or mGWAS from specific metabolite structures or biosynthetic pathways. Metabolite kernels developed from LC-MS metabolites in the discovery panel improved prediction accuracy of LC-MS metabolite traits in the validation panel consisting of more advanced breeding lines. No approach, however, improved prediction accuracy for GC-MS metabolites. We ranked model performance by metabolite and found that metabolites with similar polarity had consistent rankings of models. Overall, testing biological rationales for developing kernels for genomic prediction across populations contributes to developing frameworks for plant breeding for metabolite traits.

Breedbase: a digital ecosystem for modern plant breeding.

Modern breeding methods integrate next-generation sequencing and phenomics to identify plants with the best characteristics and greatest genetic merit for use as parents in subsequent breeding cycles to ultimately create improved cultivars able to sustain high adoption rates by farmers. This data-driven approach hinges on strong foundations in data management, quality control, and analytics. Of crucial importance is a central database able to (1) track breeding materials, (2) store experimental evaluations, (3) record phenotypic measurements using consistent ontologies, (4) store genotypic information, and (5) implement algorithms for analysis, prediction, and selection decisions. Because of the complexity of the breeding process, breeding databases also tend to be complex, difficult, and expensive to implement and maintain. Here, we present a breeding database system, Breedbase (https://breedbase.org/, last accessed 4/18/2022). Originally initiated as Cassavabase (https://cassavabase.org/, last accessed 4/18/2022) with the NextGen Cassava project (https://www.nextgencassava.org/, last accessed 4/18/2022), and later developed into a crop-agnostic system, it is presently used by dozens of different crops and projects. The system is web based and is available as open source software. It is available on GitHub (https://github.com/solgenomics/, last accessed 4/18/2022) and packaged in a Docker image for deployment (https://hub.docker.com/u/breedbase, last accessed 4/18/2022). The Breedbase system enables breeding programs to better manage and leverage their data for decision making within a fully integrated digital ecosystem.

Simulation of sugar kelp (Saccharina latissima) breeding guided by practices to accelerate genetic gains.

Though Saccharina japonica cultivation has been established for many decades in East Asian countries, the domestication process of sugar kelp (Saccharina latissima) in the Northeast United States is still at its infancy. In this study, by using data from our breeding experience, we will demonstrate how obstacles for accelerated genetic gain can be assessed using simulation approaches that inform resource allocation decisions. Thus far, we have used 140 wild sporophytes that were sampled in 2018 from the northern Gulf of Maine to southern New England. From these sporophytes, we sampled gametophytes and made and evaluated over 600 progeny sporophytes from crosses among the gametophytes in 2019 and 2020. The biphasic life cycle of kelp gives a great advantage in selective breeding as we can potentially select both on the sporophytes and gametophytes. However, several obstacles exist, such as the amount of time it takes to complete a breeding cycle, the number of gametophytes that can be maintained in the laboratory, and whether positive selection can be conducted on farm-tested sporophytes. Using the Gulf of Maine population characteristics for heritability and effective population size, we simulated a founder population of 1,000 individuals and evaluated the impact of overcoming these obstacles on rate of genetic gain. Our results showed that key factors to improve current genetic gain rely mainly on our ability to induce reproduction of the best farm-tested sporophytes, and to accelerate the clonal vegetative growth of released gametophytes so that enough gametophyte biomass is ready for making crosses by the next growing season. Overcoming these challenges could improve rates of genetic gain more than 2-fold. Future research should focus on conditions favorable for inducing spring reproduction, and on increasing the amount of gametophyte tissue available in time to make fall crosses in the same year.

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

KEY MESSAGE: 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.

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.

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

KEY MESSAGE: 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.

Development of the Wheat Practical Haplotype Graph database as a resource for genotyping data storage and genotype imputation.

To improve the efficiency of high-density genotype data storage and imputation in bread wheat (Triticum aestivum L.), we applied the Practical Haplotype Graph (PHG) tool. The Wheat PHG database was built using whole-exome capture sequencing data from a diverse set of 65 wheat accessions. Population haplotypes were inferred for the reference genome intervals defined by the boundaries of the high-quality gene models. Missing genotypes in the inference panels, composed of wheat cultivars or recombinant inbred lines genotyped by exome capture, genotyping-by-sequencing (GBS), or whole-genome skim-seq sequencing approaches, were imputed using the Wheat PHG database. Though imputation accuracy varied depending on the method of sequencing and coverage depth, we found 92% imputation accuracy with 0.01× sequence coverage, which was slightly lower than the accuracy obtained using the 0.5× sequence coverage (96.6%). Compared to Beagle, on average, PHG imputation was ∼3.5% (P-value < 2 × 10-14) more accurate, and showed 27% higher accuracy at imputing a rare haplotype introgressed from a wild relative into wheat. We found reduced accuracy of imputation with independent 2× GBS data (88.6%), which increases to 89.2% with the inclusion of parental haplotypes in the database. The accuracy reduction with GBS is likely associated with the small overlap between GBS markers and the exome capture dataset, which was used for constructing PHG. The highest imputation accuracy was obtained with exome capture for the wheat D genome, which also showed the highest levels of linkage disequilibrium and proportion of identity-by-descent regions among accessions in the PHG database. We demonstrate that genetic mapping based on genotypes imputed using PHG identifies SNPs with a broader range of effect sizes that together explain a higher proportion of genetic variance for heading date and meiotic crossover rate compared to previous studies.

Selection for seed size has uneven effects on specialized metabolite abundance in oat (Avena sativa L.).

Plant breeding strategies to optimize metabolite profiles are necessary to develop health-promoting food crops. In oats (Avena sativa L.), seed metabolites are of interest for their antioxidant properties, yet have not been a direct target of selection in breeding. In a diverse oat germplasm panel spanning a century of breeding, we investigated the degree of variation of these specialized metabolites and how it has been molded by selection for other traits, like yield components. We also ask if these patterns of variation persist in modern breeding pools. Integrating genomic, transcriptomic, metabolomic, and phenotypic analyses for three types of seed specialized metabolites-avenanthramides, avenacins, and avenacosides-we found reduced heritable genetic variation in modern germplasm compared with diverse germplasm, in part due to increased seed size associated with more intensive breeding. Specifically, we found that abundance of avenanthramides increases with seed size, but additional variation is attributable to expression of biosynthetic enzymes. In contrast, avenacoside abundance decreases with seed size and plant breeding intensity. In addition, these different specialized metabolites do not share large-effect loci. Overall, we show that increased seed size associated with intensive plant breeding has uneven effects on the oat seed metabolome, but variation also exists independently of seed size to use in plant breeding. This work broadly contributes to our understanding of how plant breeding has influenced plant traits and tradeoffs between traits (like growth and defense) and the genetic bases of these shifts.

Sexual dimorphism and the effect of wild introgressions on recombination in cassava (Manihot esculenta Crantz) breeding germplasm.

Recombination has essential functions in meiosis, evolution, and breeding. The frequency and distribution of crossovers dictate the generation of new allele combinations and can vary across species and between sexes. Here, we examine recombination landscapes across the 18 chromosomes of cassava (Manihot esculenta Crantz) with respect to male and female meioses and known introgressions from the wild relative Manihot glaziovii. We used SHAPEIT2 and duoHMM to infer crossovers from genotyping-by-sequencing data and a validated multigenerational pedigree from the International Institute of Tropical Agriculture cassava breeding germplasm consisting of 7020 informative meioses. We then constructed new genetic maps and compared them to an existing map previously constructed by the International Cassava Genetic Map Consortium. We observed higher recombination rates in females compared to males, and lower recombination rates in M. glaziovii introgression segments on chromosomes 1 and 4, with suppressed recombination along the entire length of the chromosome in the case of the chromosome 4 introgression. Finally, we discuss hypothesized mechanisms underlying our observations of heterochiasmy and crossover suppression and discuss the broader implications for plant breeding.

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.

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.