The Cassava Genome: Current Progress, Future Directions

A Coiled-Coil Nucleotide-Binding Domain Leucine-Rich Repeat Receptor Gene MeRPPL1 Plays a Role in the Replication of a Geminivirus in Cassava

Disease resistance gene (R gene)-encoded nucleotide-binding leucine-rich repeat proteins (NLRs) are critical players in plant host defence mechanisms because of their role as receptors that recognise pathogen effectors and trigger plant effector-triggered immunity (ETI). This study aimed to determine the putative role of a cassava coiled-coil (CC)-NLR (CNL) gene MeRPPL1 (Manes.12G091600) (single allele) located on chromosome 12 in the tolerance or susceptibility to South African cassava mosaic virus (SACMV), one of the causal agents of cassava mosaic disease (CMD). A transient protoplast system was used to knock down the expression of MeRPPL1 by clustered regularly interspaced short palindromic repeats-CRISPR-associated protein 9 (CRISPR-Cas9). The MeRPPL1-targeting CRISPR vectors and/or SACMV DNA A and DNA B infectious clones were used to transfect protoplasts isolated from leaf mesophyll cells from the SACMV-tolerant cassava (Manihot esculenta) cultivar TME3. The CRISPR/Cas9 silencing vector significantly reduced MeRPPL1 expression in protoplasts whether with or without SACMV co-infection. Notably, SACMV DNA A replication was higher in protoplasts with lower MeRPPL1 expression levels than in non-silenced protoplasts. Mutagenesis studies revealed that protoplast co-transfection with CRISPR-MeRPPL1 silencing vector + SACMV and transfection with only SACMV induced nucleotide substitution mutations that led to altered amino acids in the highly conserved MHD motif of the MeRPPL1-translated polypeptide. This may abolish or alter the regulatory role of the MHD motif in controlling R protein activity and could contribute to the increase in SACMV-DNA A accumulation observed in MeRPPL1-silenced protoplasts. The results herein demonstrate for the first time a role for a CNL gene in tolerance to a geminivirus in TME3.

Comparative analysis of infected cassava root transcriptomics reveals candidate genes for root rot disease resistance

Cassava root-rot incited by soil-borne pathogens is one of the major diseases that reduces root yield. Although the use of resistant cultivars is the most effective method of management, the genetic basis for root-rot resistance remains poorly understood. Therefore, our work analyzed the transcriptome of two contrasting genotypes (BRS Kiriris/resistant and BGM-1345/susceptible) using RNA-Seq to understand the molecular response and identify candidate genes for resistance. Cassava seedlings (resistant and susceptible to root-rot) were both planted in infested and sterilized soil and samples from Initial-time and Final-time periods, pooled. Two controls were used: (i) seedlings collected before planting in infested soil (absolute control) and, (ii) plants grown in sterilized soil (mock treatments). For the differentially expressed genes (DEGs) analysis 23.912 were expressed in the resistant genotype, where 10.307 were differentially expressed in the control treatment, 15 DEGs in the Initial Time-period and 366 DEGs in the Final Time-period. Eighteen candidate genes from the resistant genotype were related to plant defense, such as the MLP-like protein 31 and the peroxidase A2-like gene. This is the first model of resistance at the transcriptional level proposed for the cassava × root-rot pathosystem. Gene validation will contribute to screening for resistance of germplasm, segregating populations and/or use in gene editing in the pursuit to develop most promising cassava clones with resistance to root-rot.

Whole-genome sequencing and comparative genomics reveal candidate genes associated with quality traits in Dioscorea alata

BACKGROUND

Quality traits are essential determinants of consumer preferences. Dioscorea alata (Greater Yam), is a starchy tuber crop in tropical regions. However, a comprehensive understanding of the genetic basis underlying yam tuber quality remains elusive. To address this knowledge gap, we employed population genomics and candidate gene association approaches to unravel the genetic factors influencing the quality attributes of boiled yam.

METHODS AND RESULTS

Comparative genomics analysis of 45 plant species revealed numerous novel genes absent in the existing D. alata gene annotation. This approach, adding 48% more genes, significantly enhanced the functional annotation of three crucial metabolic pathways associated with boiled yam quality traits: pentose and glucuronate interconversions, starch and sucrose metabolism, and flavonoid biosynthesis. In addition, the whole-genome sequencing of 127 genotypes identified 27 genes under selection and 22 genes linked to texture, starch content, and color through a candidate gene association analysis. Notably, five genes involved in starch content and cell wall composition, including 1,3-beta Glucan synthase, β-amylase, and Pectin methyl esterase, were common to both approaches and their expression levels were assessed by transcriptomic data.

CONCLUSIONS

The analysis of the whole-genome of 127 genotypes of D. alata and the study of three specific pathways allowed the identification of important genes for tuber quality. Our findings provide insights into the genetic basis of yam quality traits and will help the enhancement of yam tuber quality through breeding programs.

Genomic and Expression Analysis of Cassava (Manihot esculenta Crantz) Chalcone Synthase Genes in Defense against Tetranychus cinnabarinus Infestation

Cassava is susceptible to mites, especially Tetranychus cinnabarinus. Secondary metabolism products such as flavonoids play an important role as antimicrobial metabolites protecting plants against biotic stressors including fungal, pathogen, bacterial, and pest defense. The chalcone synthase (CHS) is the initial step of the phenylpropanoid pathway for producing flavonoids and is the gatekeeper of the pathway. Until recently, the CHS genes family has not been systematically studied in cassava. Thirty-nine CHS genes were identified from the cassava genome database. Based on phylogenetic and sequence composition analysis, these CHSs were divided into 3 subfamilies. Within the same subfamily, the gene structure and motif compositions of these CHS genes were found to be quite conserved. Duplication events, particularly segmental duplication of the cassava CHS genes, were identified as one of the main driving force of its expansion. Various cis-elements contained in the promoter might regulate the gene expression patterns of MeCHS. Protein-protein interaction (PPI) network analysis showed that MeCHS1 and MeCHS10 protein are more closely related to other family members. The expression of MeCHS genes in young leaves was higher than that in other tissues, and their expression varies even within the same tissue. Coincidentally, these CHS genes of most LAP subclasses were highly expressed in young leaves. The verified MeCHS genes showed consistent with the real-time reverse transcription quantitative PCR (RT-qPCR) and proteomic expression in protected and affected leaves respectively, indicating that these MeCHS genes play crucial roles in the response to T. cinnabarinus. This study is the first to comprehensively expatiate the information on MeCHS family members. These data will further enhance our understanding both the molecular mechanisms and the effects of CHS genes. In addition, the results will help to further clarify the effects on T. cinnabarinus and provide a theoretical basis for the potential functions of the specific CHS gene in resistance to mites and other biotic stress.

Identification of MeC3HDZ1/MeCNA as a potential regulator of cassava storage root development

The storage root (SR) of cassava is the main staple food in sub-Saharan Africa, where it feeds over 500 million people. However, little is known about the genetic and molecular regulation underlying its development. Unraveling such regulation would pave the way for biotechnology approaches aimed at enhancing cassava productivity. Anatomical studies indicate that SR development relies on the massive accumulation of xylem parenchyma, a cell-type derived from the vascular cambium. The C3HDZ family of transcription factors regulate cambial cells proliferation and xylem differentiation in Arabidopsis and other species. We thus aimed at identifying C3HDZ proteins in cassava and determining whether any of them shows preferential activity in the SR cambium and/or xylem. Using phylogeny and synteny studies, we identified eight C3HDZ proteins in cassava, namely MeCH3DZ1-8. We observed that MeC3HDZ1 is the MeC3HDZ gene displaying the highest expression in SR and that, within that organ, the gene also shows high expression in cambium and xylem. In-silico analyses revealed the existence of a number of potential C3HDZ targets displaying significant preferential expression in the SR. Subsequent Y1H analyses proved that MeC3HDZ1 can bind canonical C3HDZ binding sites, present in the promoters of these targets. Transactivation assays demonstrated that MeC3HDZ1 can regulate the expression of genes downstream of promoters harboring such binding sites, thereby demonstrating that MeC3HDZ1 has C3HDZ transcription factor activity. We conclude that MeC3HDZ1 may be a key factor for the regulation of storage root development in cassava, holding thus great promise for future biotechnology applications.

Integrating genomics and genome editing for orphan crop improvement: a bridge between orphan crops and modern agriculture system

Domestication of orphan crops could be explored by editing their genomes. Genome editing has a lot of promise for enhancing agricultural output, and there is a lot of interest in furthering breeding in orphan crops, which are sometimes plagued with unwanted traits that resemble wild cousins. Consequently, applying model crop knowledge to orphan crops allows for the rapid generation of targeted allelic diversity and innovative breeding germplasm. We explain how plant breeders could employ genome editing as a novel platform to accelerate the domestication of semi-domesticated or wild plants, resulting in a more diversified base for future food and fodder supplies. This review emphasizes both the practicality of the strategy and the need to invest in research that advances our understanding of plant genomes, genes, and cellular systems. Planting more of these abandoned orphan crops could help alleviate food scarcities in the challenge of future climate crises.

Haplotype-resolved genome of heterozygous African cassava cultivar TMEB117 (Manihot esculenta)

Cassava (Manihot esculenta Crantz) is a vital tropical root crop providing essential dietary energy to over 800 million people in tropical and subtropical regions. As a climate-resilient crop, its significance grows as the human population expands. However, yield improvement faces challenges from biotic and abiotic stress and limited breeding. Advanced sequencing and assembly techniques enabled the generation of a highly accurate, nearly complete, haplotype-resolved genome of the African cassava cultivar TMEB117. It is the most accurate cassava genome sequence to date with a base-level accuracy of QV > 64, N50 > 35 Mbp, and 98.9% BUSCO completeness. Over 60% of the genome comprises repetitive elements. We predicted over 45,000 gene models for both haplotypes. This achievement offers valuable insights into the heterozygosity genome organization of the cassava genome, with improved accuracy, completeness, and phased genomes. Due to its high susceptibility to African Cassava Mosaic Virus (ACMV) infections compared to other cassava varieties, TMEB117 provides an ideal reference for studying virus resistance mechanisms, including epigenetic variations and smallRNA expressions.

Superoxide dismutase gene family in cassava revealed their involvement in environmental stress via genome-wide analysis

Superoxide dismutase (SOD) is a crucial metal-containing enzyme that plays a vital role in catalyzing the dismutation of superoxide anions, converting them into molecular oxygen and hydrogen peroxide, essential for enhancing plant stress tolerance. We identified 8 SOD genes (4 CSODs, 2 FSODs, and 2 MSODs) in cassava. Bioinformatics analyses provided insights into chromosomal location, phylogenetic relationships, gene structure, conserved motifs, and gene ontology annotations. MeSOD genes were classified into two groups through phylogenetic analysis, revealing evolutionary connections. Promoters of these genes harbored stress-related cis-elements. Duplication analysis indicated the functional significance of MeCSOD2/MeCSOD4 and MeMSOD1/MeMSOD2. Through qRT-PCR, MeCSOD2 responded to salt stress, MeMSOD2 to drought, and cassava bacterial blight. Silencing MeMSOD2 increased XpmCHN11 virulence, indicating MeMSOD2 is essential for cassava's defense against XpmCHN11 infection. These findings enhance our understanding of the SOD gene family's role in cassava and contribute to strategies for stress tolerance improvement.

Allele mining through TILLING and EcoTILLING approaches in vegetable crops

MAIN CONCLUSION

The present review illustrates a comprehensive overview of the allele mining for genetic improvement in vegetable crops, and allele exploration methods and their utilization in various applications related to pre-breeding of economically important traits in vegetable crops. Vegetable crops have numerous wild descendants, ancestors and terrestrial races that could be exploited to develop high-yielding and climate-resilient varieties resistant/tolerant to biotic and abiotic stresses. To further boost the genetic potential of economic traits, the available genomic tools must be targeted and re-opened for exploitation of novel alleles from genetic stocks by the discovery of beneficial alleles from wild relatives and their introgression to cultivated types. This capability would be useful for giving plant breeders direct access to critical alleles that confer higher production, improve bioactive compounds, increase water and nutrient productivity as well as biotic and abiotic stress resilience. Allele mining is a new sophisticated technique for dissecting naturally occurring allelic variants in candidate genes that influence important traits which could be used for genetic improvement of vegetable crops. Target-induced local lesions in genomes (TILLINGs) is a sensitive mutation detection avenue in functional genomics, particularly wherein genome sequence information is limited or not available. Population exposure to chemical mutagens and the absence of selectivity lead to TILLING and EcoTILLING. EcoTILLING may lead to natural induction of SNPs and InDels. It is anticipated that as TILLING is used for vegetable crops improvement in the near future, indirect benefits will become apparent. Therefore, in this review we have highlighted the up-to-date information on allele mining for genetic enhancement in vegetable crops and methods of allele exploration and their use in pre-breeding for improvement of economic traits.

Genome assemblies and comparison of two Neotropical spiral gingers: Costus pulverulentus and C. lasius

The spiral gingers (Costus L.) are a pantropical genus of herbaceous perennial monocots; the Neotropical clade of Costus radiated rapidly in the past few million years into over 60 species. The Neotropical spiral gingers have a rich history of evolutionary and ecological research that can motivate and inform modern genetic investigations. Here, we present the first 2 chromosome-level genome assemblies in the genus, for C. pulverulentus and C. lasius, and briefly compare their synteny. We assembled the C. pulverulentus genome from a combination of short-read data, Chicago and Dovetail Hi-C chromatin-proximity sequencing, and alignment with a linkage map. We annotated the genome by mapping a C. pulverulentus transcriptome and querying mapped transcripts against a protein database. We assembled the C. lasius genome with Pacific Biosciences HiFi long reads and alignment to the C. pulverulentus genome. These 2 assemblies are the first published genomes for non-cultivated tropical plants. These genomes solidify the spiral gingers as a model system and will facilitate research on the poorly understood genetic basis of tropical plant diversification.

Fully sequencing the cassava full-length cDNA library reveals unannotated transcript structures and alternative splicing events in regions with a high density of single nucleotide variations, insertions-deletions, and heterozygous sequences

The primary transcript structure provides critical insights into protein diversity, transcriptional modification, and functions. Cassava transcript structures are highly diverse because of alternative splicing (AS) events and high heterozygosity. To precisely determine and characterize transcript structures, fully sequencing cloned transcripts is the most reliable method. However, cassava annotations were mainly determined according to fragmentation-based sequencing analyses (e.g., EST and short-read RNA-seq). In this study, we sequenced the cassava full-length cDNA library, which included rare transcripts. We obtained 8,628 non-redundant fully sequenced transcripts and detected 615 unannotated AS events and 421 unannotated loci. The different protein sequences resulting from the unannotated AS events tended to have diverse functional domains, implying that unannotated AS contributes to the truncation of functional domains. The unannotated loci tended to be derived from orphan genes, implying that the loci may be associated with cassava-specific traits. Unexpectedly, individual cassava transcripts were more likely to have multiple AS events than Arabidopsis transcripts, suggestive of the regulated interactions between cassava splicing-related complexes. We also observed that the unannotated loci and/or AS events were commonly in regions with abundant single nucleotide variations, insertions-deletions, and heterozygous sequences. These findings reflect the utility of completely sequenced FLcDNA clones for overcoming cassava-specific annotation-related problems to elucidate transcript structures. Our work provides researchers with transcript structural details that are useful for annotating highly diverse and unique transcripts and alternative splicing events.

Comparative Transcriptome Profiling of Cassava Tuberous Roots in Response to Postharvest Physiological Deterioration

Cassava is one of the most versatile tuberous-root crops on Earth. However, the postharvest storage properties of cassava tuberous root mean that it is perishable through a process known as postharvest physiological deterioration (PPD), which seriously affects its starch quality. Therefore, a comprehensive understanding of the transcriptional regulatory activity of cassava against the PPD response is necessary in order to extract key molecular mechanisms related to PPD tolerance. In this study, we found that RYG1 tuberous roots showed delayed PPD compared to those of SC8. In addition, RYG1 roots maintained a more stable cell wall structure after storage than those of SC8. The transcriptome changes in tuberous roots were analyzed for both RYG1 and SC8 after 21 days of storage (SR and SS) compared to fresh (FR and FS) by the RNA-Seq method. The total number of differentially expressed genes (DEGs) in the various comparisons of these four samples ranged from 68 to 3847. Of these, a total of 2008 co-DEGs in SR vs. SS were shared by either SR vs. FR or SS vs. FS. GO and KEGG enrichment analysis revealed that upregulated co-DEGs in SR vs. SS were mainly enriched in photosynthesis, protein processing, hormone and cutin, suberine and wax biosynthesis. By contrast, the downregulated co-DEGs were mainly related to cell wall organization, starch and sucrose metabolism, galactose metabolism, phenylpropanoid biosynthesis, diterpenoid biosynthesis, cysteine and methionine metabolism and flavonoid biosynthesis. The protein-protein interaction (PPI) networks of the co-DEGs showed a complex interaction of genes in different pathways, and 16 hub genes were characterized to have a degree in excess of 15, among which eight genes were associated with photosynthesis. These results provide new information for the study of cassava resistance to PPD and lay a foundation for the further molecular breeding of storage-tolerant cassava varieties.

On-Farm Multi-Environment Evaluation of Selected Cassava (Manihot esculenta Crantz) Cultivars in South Africa

Cassava is an important starchy root crop grown globally in tropical and subtropical regions. The ability of cassava to withstand difficult growing conditions and long-term storability underground makes it a resilient crop, contributing to food and nutrient security. This study was conducted to evaluate the performance and adaptability of exotic cassava cultivars across different environments in South Africa and to recommend genotypes for cultivation. A total of 11 cassava cultivars were evaluated at six on-farm sites, using a randomized complete block design with three replications. There were highly significant (p < 0.001) variations between genotypes, environments, and their interaction for all yield and yield-related traits studied. This indicates the need to test the genotypes in multiple environments before effective selection and commercialization can be undertaken. MSAF2 and UKF4 showed the overall best performances for most of the traits, whilst UKF9 (49.5%) and P1/19 (48.5%) had the highest dry matter yield. UKF4 (102.7 t ha−1) had the highest yield and greatest root yield stability across environments. MSAF2 did not perform consistently across environments because it was highly susceptible to cassava mosaic disease (CMD). MSAF2 could be used as a donor parent to generate novel clones with large numbers of marketable roots, and high fresh root yields, if the other parent can provide effective resistance to CMD. Based on genotype and environmental mean, Mabuyeni (KwaZulu-Natal), Mandlakazi (Limpopo), and Shatale (Mpumalanga) were found to be better environments for cassava cultivation and testing. This study is a pioneer in cassava research using multiple environments in South Africa. It provides baseline information on the performance of currently available cassava clones, their adaptation to multiple sites, the identification of suitable test sites, and information on current genetic resources for a future breeding program.

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.

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.

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.

Genome-Wide Identification and Comparative Analysis of WOX Genes in Four Euphorbiaceae Species and Their Expression Patterns in Jatropha curcas

The WUSCHEL-related homeobox (WOX) proteins are widely distributed in plants and play important regulatory roles in growth and development processes such as embryonic development and organ development. Here, series of bioinformatics methods were utilized to unravel the structural basis and genetic hierarchy of WOX genes, followed by regulation of the WOX genes in four Euphorbiaceae species. A genome-wide survey identified 59 WOX genes in Hevea brasiliensis (H. brasiliensis: 20 genes), Jatropha curcas (J. curcas: 10 genes), Manihot esculenta (M. esculenta: 18 genes), and Ricinus communis (R. communis: 11 genes). The phylogenetic analysis revealed that these WOX members could be clustered into three close proximal clades, such as namely ancient, intermediate and modern/WUS clades. In addition, gene structures and conserved motif analyses further validated that the WOX genes were conserved within each phylogenetic clade. These results suggested the relationships among WOX members in the four Euphorbiaceae species. We found that WOX genes in H. brasiliensis and M. esculenta exhibit close genetic relationship with J. curcas and R. communis. Additionally, the presence of various cis-acting regulatory elements in the promoter of J. curcas WOX genes (JcWOXs) reflected distinct functions. These speculations were further validated with the differential expression profiles of various JcWOXs in seeds, reflecting the importance of two JcWOX genes (JcWOX6 and JcWOX13) during plant growth and development. Our quantitative real-time PCR (qRT-PCR) analysis demonstrated that the JcWOX11 gene plays an indispensable role in regulating plant callus. Taken together, the present study reports the comprehensive characteristics and relationships of WOX genes in four Euphorbiaceae species, providing new insights into their characterization.

Advances in Genetic Analysis and Breeding of Cassava (Manihot esculenta Crantz): A Review

Cassava (Manihot esculenta Crantz) is the sixth most important food crop and consumed by 800 million people worldwide. In Africa, cassava is the second most important food crop after maize and Africa is the worlds' largest producer. Though cassava is not one of the main commodity crops in South Africa, it is becoming a popular crop among farming communities in frost-free areas, due to its climate-resilient nature. This necessitated the establishment of a multi-disciplinary research program at the Agricultural Research Council of South Africa. The objective of this review is to highlight progress made in cassava breeding and genetic analysis. This review highlights the progress of cassava research worldwide and discusses research findings on yield, quality, and adaptability traits in cassava. It also discusses the limitations and the prospects of the cassava R&D program towards development of the cassava industry in South Africa.

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.

The haplotype-resolved chromosome pairs of a heterozygous diploid African cassava cultivar reveal novel pan-genome and allele-specific transcriptome features

BACKGROUND

Cassava (Manihot esculenta) is an important clonally propagated food crop in tropical and subtropical regions worldwide. Genetic gain by molecular breeding has been limited, partially because cassava is a highly heterozygous crop with a repetitive and difficult-to-assemble genome.

FINDINGS

Here we demonstrate that Pacific Biosciences high-fidelity (HiFi) sequencing reads, in combination with the assembler hifiasm, produced genome assemblies at near complete haplotype resolution with higher continuity and accuracy compared to conventional long sequencing reads. We present 2 chromosome-scale haploid genomes phased with Hi-C technology for the diploid African cassava variety TME204. With consensus accuracy >QV46, contig N50 >18 Mb, BUSCO completeness of 99%, and 35k phased gene loci, it is the most accurate, continuous, complete, and haplotype-resolved cassava genome assembly so far. Ab initio gene prediction with RNA-seq data and Iso-Seq transcripts identified abundant novel gene loci, with enriched functionality related to chromatin organization, meristem development, and cell responses. During tissue development, differentially expressed transcripts of different haplotype origins were enriched for different functionality. In each tissue, 20-30% of transcripts showed allele-specific expression (ASE) differences. ASE bias was often tissue specific and inconsistent across different tissues. Direction-shifting was observed in <2% of the ASE transcripts. Despite high gene synteny, the HiFi genome assembly revealed extensive chromosome rearrangements and abundant intra-genomic and inter-genomic divergent sequences, with large structural variations mostly related to LTR retrotransposons. We use the reference-quality assemblies to build a cassava pan-genome and demonstrate its importance in representing the genetic diversity of cassava for downstream reference-guided omics analysis and breeding.

CONCLUSIONS

The phased and annotated chromosome pairs allow a systematic view of the heterozygous diploid genome organization in cassava with improved accuracy, completeness, and haplotype resolution. They will be a valuable resource for cassava breeding and research. Our study may also provide insights into developing cost-effective and efficient strategies for resolving complex genomes with high resolution, accuracy, and continuity.

High-density genetic map and genome-wide association studies of aesthetic traits in Phalaenopsis orchids

Phalaenopsis spp. represent the most popular orchids worldwide. Both P. equestris and P. aphrodite are the two important breeding parents with the whole genome sequence available. However, marker–trait association is rarely used for floral traits in Phalaenopsis breeding. Here, we analyzed markers associated with aesthetic traits of Phalaenopsis orchids by using genome-wide association study (GWAS) with the F1 population P. Intermedia of 117 progenies derived from the cross between P. aphrodite and P. equestris. A total of 113,517 single nucleotide polymorphisms (SNPs) were identified in P. Intermedia by using genotyping-by-sequencing with the combination of two different restriction enzyme pairs, Hinp1 I/Hae III and Apek I/Hae III. The size-related traits from flowers were negatively related to the color-related traits. The 1191 SNPs from Hinp1 I/ Hae III and 23 simple sequence repeats were used to establish a high-density genetic map of 19 homolog groups for P. equestris. In addition, 10 quantitative trait loci were highly associated with four color-related traits on chromosomes 2, 5 and 9. According to the sequence within the linkage disequilibrium regions, 35 candidate genes were identified and related to anthocyanin biosynthesis. In conclusion, we performed marker-assisted gene identification of aesthetic traits with GWAS in Phalaenopsis orchids.

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.

Characterization of cassava ORANGE proteins and their capability to increase provitamin A carotenoids accumulation

Cassava (Manihot esculenta Crantz) biofortification with provitamin A carotenoids is an ongoing process that aims to alleviate vitamin A deficiency. The moderate content of provitamin A carotenoids achieved so far limits the contribution to providing adequate dietary vitamin A levels. Strategies to increase carotenoid content focused on genes from the carotenoids biosynthesis pathway. In recent years, special emphasis was given to ORANGE protein (OR), which promotes the accumulation of carotenoids and their stability in several plants. The aim of this work was to identify, characterize and investigate the role of OR in the biosynthesis and stabilization of carotenoids in cassava and its relationship with phytoene synthase (PSY), the rate-limiting enzyme of the carotenoids biosynthesis pathway. Gene and protein characterization of OR, expression levels, protein amounts and carotenoids levels were evaluated in roots of one white (60444) and two yellow cassava cultivars (GM5309-57 and GM3736-37). Four OR variants were found in yellow cassava roots. Although comparable expression was found for three variants, significantly higher OR protein amounts were observed in the yellow varieties. In contrast, cassava PSY1 expression was significantly higher in the yellow cultivars, but PSY protein amount did not vary. Furthermore, we evaluated whether expression of one of the variants, MeOR_X1, affected carotenoid accumulation in cassava Friable Embryogenic Callus (FEC). Overexpression of maize PSY1 alone resulted in carotenoids accumulation and induced crystal formation. Co-expression with MeOR_X1 led to greatly increase of carotenoids although PSY1 expression was high in the co-expressed FEC. Our data suggest that posttranslational mechanisms controlling OR and PSY protein stability contribute to higher carotenoid levels in yellow cassava. Moreover, we showed that cassava FEC can be used to study the efficiency of single and combinatorial gene expression in increasing the carotenoid content prior to its application for the generation of biofortified cassava with enhanced carotenoids levels.

Gramene: A Resource for Comparative Analysis of Plants Genomes and Pathways

Gramene is an integrated bioinformatics resource for accessing, visualizing, and comparing plant genomes and biological pathways. Originally targeting grasses, Gramene has grown to host annotations for over 90 plant genomes including agronomically important cereals (e.g., maize, sorghum, wheat, teff), fruits and vegetables (e.g., apple, watermelon, clementine, tomato, cassava), specialty crops (e.g., coffee, olive tree, pistachio, almond), and plants of special or emerging interest (e.g., cotton, tobacco, cannabis, or hemp). For some species, the resource includes multiple varieties of the same species, which has paved the road for the creation of species-specific pan-genome browsers. The resource also features plant research models, including Arabidopsis and C4 warm-season grasses and brassicas, as well as other species that fill phylogenetic gaps for plant evolution studies. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. This chapter outlines system requirements for end-users and database hosting, data types and basic navigation within Gramene, and provides examples of how to (1) explore Gramene's search results, (2) explore gene-centric comparative genomics data visualizations in Gramene, and (3) explore genetic variation associated with a gene locus. This is the first publication describing in detail Gramene's integrated search interface-intended to provide a simplified entry portal for the resource's main data categories (genomic location, phylogeny, gene expression, pathways, and external references) to the most complete and up-to-date set of plant genome and pathway annotations.

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.

Utilizing CRISPR-Cas in Tropical Crop Improvement: A Decision Process for Fitting Genome Engineering to Your Species

Adopting modern gene-editing technologies for trait improvement in agriculture requires important workflow developments, yet these developments are not often discussed. Using tropical crop systems as a case study, we describe a workflow broken down into discrete processes with specific steps and decision points that allow for the practical application of the CRISPR-Cas gene editing platform in a crop of interest. While we present the steps of developing genome-edited plants as sequential, in practice parts can be done in parallel, which are discussed in this perspective. The main processes include 1) understanding the genetic basis of the trait along with having the crop’s genome sequence, 2) testing and optimization of the editing reagents, development of efficient 3) tissue culture and 4) transformation methods, and 5) screening methods to identify edited events with commercial potential. Our goal in this perspective is to help any lab that wishes to implement this powerful, easy-to-use tool in their pipeline, thus aiming to democratize the technology.

Large structural variations in the haplotype-resolved African cassava genome

Cassava (Manihot esculenta Crantz, 2n = 36) is a global food security crop. It has a highly heterozygous genome, high genetic load, and genotype-dependent asynchronous flowering. It is typically propagated by stem cuttings and any genetic variation between haplotypes, including large structural variations, is preserved by such clonal propagation. Traditional genome assembly approaches generate a collapsed haplotype representation of the genome. In highly heterozygous plants, this results in artifacts and an oversimplification of heterozygous regions. We used a combination of Pacific Biosciences (PacBio), Illumina, and Hi-C to resolve each haplotype of the genome of a farmer-preferred cassava line, TME7 (Oko-iyawo). PacBio reads were assembled using the FALCON suite. Phase switch errors were corrected using FALCON-Phase and Hi-C read data. The ultralong-range information from Hi-C sequencing was also used for scaffolding. Comparison of the two phases revealed >5000 large haplotype-specific structural variants affecting over 8 Mb, including insertions and deletions spanning thousands of base pairs. The potential of these variants to affect allele-specific expression was further explored. RNA-sequencing data from 11 different tissue types were mapped against the scaffolded haploid assembly and gene expression data are incorporated into our existing easy-to-use web-based interface to facilitate use by the broader plant science community. These two assemblies provide an excellent means to study the effects of heterozygosity, haplotype-specific structural variation, gene hemizygosity, and allele-specific gene expression contributing to important agricultural traits and further our understanding of the genetics and domestication of cassava.

Cassava diseases caused by Xanthomonas phaseoli pv. manihotis and Xanthomonas cassavae

Xanthomonas phaseoli pv. manihotis (Xpm) and X. cassavae (Xc) are two bacterial pathogens attacking cassava. Cassava bacterial blight (CBB) is a systemic disease caused by Xpm, which might have dramatic effects on plant growth and crop production. Cassava bacterial necrosis is a nonvascular disease caused by Xc with foliar symptoms similar to CBB, but its impacts on the plant vigour and the crop are limited. In this review, we describe the epidemiology and ecology of the two pathogens, the impacts and management of the diseases, and the main research achievements for each pathosystem. Because Xc data are sparse, our main focus is on Xpm and CBB.

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.

Adoption and Promotion of Resilient Crops for Climate Risk Mitigation and Import Substitution: A Case Analysis of Cassava for South African Agriculture

Cassava is an important starchy root crop grown globally in tropical and subtropical regions. The ability of cassava to withstand difficult growing conditions and long-term storability underground makes it a resilient crop, contributing to food security. Historically, small-scale farmers have grown cassava as a minor crop in the far north-eastern part of the country. However, there is an initiative to scale up cassava production, with two discrete areas of interest: large-scale production for industrial starch, and expanding its footprint as a food security crop for small-scale farmers, especially in the context of climate change. In this scoping study, production, processing and marketing data for cassava were accessed from the FAO and US Commercial trade databases. Other domestic market and demand analysis case studies were also explored. There is no cassava data available for South Africa. The study indicated that South Africa imports more than 66,000 tons of starch annually, of which 33% is cassava starch, showing the availability of a local market. The potential of cassava for the South African economy is discussed. Significant industrial opportunities exist for the production and use of cassava in South Africa. However, the realization of these opportunities will depend on the reliable supply of good quality cassava roots. However, the lack of a well-established cassava research program, and a lack of an existing value chain for the industrial scale cassava production and processing are barriers to the development of cassava industry in South Africa. As the initial step to the development of a successful cassava industry, high potential germplasm is imported, characterized and bred for local conditions to ensure the sustainable primary production of cassava. Subsequently, industrial value chains will need to be developed as the optimization of the breeding and agronomy of the crop are completed, and yield potentials are quantified in the different regions of the country.

Citrus Tristeza Virus Genotype Detection Using High-Throughput Sequencing

The application of high-throughput sequencing (HTS) has successfully been used for virus discovery to resolve disease etiology in many agricultural crops. The greatest advantage of HTS is that it can provide a complete viral status of a plant, including information on mixed infections of viral species or virus variants. This provides insight into the virus population structure, ecology, or evolution and can be used to differentiate among virus variants that may contribute differently toward disease etiology. In this study, the use of HTS for citrus tristeza virus (CTV) genotype detection was evaluated. A bioinformatic pipeline for CTV genotype detection was constructed and evaluated using simulated and real data sets to determine the parameters to discriminate between false positive read mappings and true genotype-specific genome coverage. A 50% genome coverage cut-off was identified for non-target read mappings. HTS with the associated bioinformatic pipeline was validated and proposed as a CTV genotyping assay.

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.

The Cassava Source-Sink project: opportunities and challenges for crop improvement by metabolic engineering

Cassava (Manihot esculenta Crantz) is one of the important staple foods in Sub-Saharan Africa. It produces starchy storage roots that provide food and income for several hundred million people, mainly in tropical agriculture zones. Increasing cassava storage root and starch yield is one of the major breeding targets with respect to securing the future food supply for the growing population of Sub-Saharan Africa. The Cassava Source-Sink (CASS) project aims to increase cassava storage root and starch yield by strategically integrating approaches from different disciplines. We present our perspective and progress on cassava as an applied research organism and provide insight into the CASS strategy, which can serve as a blueprint for the improvement of other root and tuber crops. Extensive profiling of different field-grown cassava genotypes generates information for leaf, phloem, and root metabolic and physiological processes that are relevant for biotechnological improvements. A multi-national pipeline for genetic engineering of cassava plants covers all steps from gene discovery, cloning, transformation, molecular and biochemical characterization, confined field trials, and phenotyping of the seasonal dynamics of shoot traits under field conditions. Together, the CASS project generates comprehensive data to facilitate conventional breeding strategies for high-yielding cassava genotypes. It also builds the foundation for genome-scale metabolic modelling aiming to predict targets and bottlenecks in metabolic pathways. This information is used to engineer cassava genotypes with improved source-sink relations and increased yield potential.

Production of very-high-amylose cassava by post-transcriptional silencing of branching enzyme genes

High amylose starch can be produced by plants deficient in the function of branching enzymes (BEs). Here we report the production of transgenic cassava (Manihot esculenta Crantz) with starches containing up to 50% amylose due to the constitutive expression of hair-pin dsRNAs targeting the BE1 or BE2 genes. All BE1-RNAi plant lines (BE1i) and BE2-RNAi plant lines (BE2i) were grown up in the field, but with reduced total biomass production. Considerably high amylose content in the storage roots of BE2i plant lines was achieved. Storage starch granules of BE1i and BE2i plants had similar morphology as wild type (WT), however, the size of BE1i starch granules were bigger than that of WT. Comparisons of amylograms and thermograms of all three sources of storage starches revealed dramatic changes to the pasting properties and a higher melting temperature for BE2i starches. Glucan chain length distribution analysis showed a slight increase in chains of DP>36 in BE1i lines and a dramatic increase in glucan chains between DP 10-20 and DP>40 in BE2i lines. Furthermore, BE2i starches displayed a B-type X-ray diffraction pattern instead of the A-type pattern found in BE1i and WT starches. Therefore, cassava BE1 and BE2 function differently in storage root starch biosynthesis.

The chromosome-scale genome reveals the evolution and diversification after the recent tetraploidization event in tea plant

Tea is one of the most popular nonalcoholic beverages due to its characteristic secondary metabolites with numerous health benefits. Although two draft genomes of tea plant (Camellia sinensis) have been published recently, the lack of chromosome-scale assembly hampers the understanding of the fundamental genomic architecture of tea plant and potential improvement. Here, we performed a genome-wide chromosome conformation capture technique (Hi-C) to obtain a chromosome-scale assembly based on the draft genome of C. sinensis var. sinensis and successfully ordered 2984.7 Mb (94.7%) scaffolds into 15 chromosomes. The scaffold N50 of the improved genome was 218.1 Mb, ~157-fold higher than that of the draft genome. Collinearity comparison of genome sequences and two genetic maps validated the high contiguity and accuracy of the chromosome-scale assembly. We clarified that only one Camellia recent tetraploidization event (CRT, 58.9-61.7 million years ago (Mya)) occurred after the core-eudicot common hexaploidization event (146.6-152.7 Mya). Meanwhile, 9243 genes (28.6%) occurred in tandem duplication, and most of these expanded after the CRT event. These gene duplicates increased functionally divergent genes that play important roles in tea-specific biosynthesis or stress response. Sixty-four catechin- and caffeine-related quantitative trait loci (QTLs) were anchored to chromosome assembly. Of these, two catechin-related QTL hotspots were derived from the CRT event, which illustrated that polyploidy has played a dramatic role in the diversification of tea germplasms. The availability of a chromosome-scale genome of tea plant holds great promise for the understanding of genome evolution and the discovery of novel genes contributing to agronomically beneficial traits in future breeding programs.

Cassava breeding and agronomy in Asia: 50 years of history and future directions

In Asia, cassava (Manihot esculenta) is cultivated by more than 8 million farmers, driving the rural economy of many countries. The International Center for Tropical Agriculture (CIAT), in partnership with national agricultural research institutes (NARIs), instigated breeding and agronomic research in Asia, 1983. The breeding program has successfully released high-yielding cultivars resulting in an average yield increase from 13.0 t ha-1 in 1996 to 21.3 t ha-1 in 2016, with significant economic benefits. Following the success in increasing yields, cassava breeding has turned its focus to higher-value traits, such as waxy cassava, to reach new market niches. More recently, building resistance to invasive pests and diseases has become a top priority due to the emergent threat of cassava mosaic disease (CMD). The agronomic research involves driving profitability with advanced technologies focusing on better agronomic management practices thereby maintaining sustainable production systems. Remote sensing technologies are being tested for trait discovery and large-scale field evaluation of cassava. In summary, cassava breeding in Asia is driven by a combination of food and market demand with technological innovations to increase the productivity. Further, exploration in the potential of data-driven agriculture is needed to empower researchers and producers for sustainable advancement.

A population genomics appraisal suggests independent dispersals for bitter and sweet manioc in Brazilian Amazonia

Amazonia is a major world centre of plant domestication, but the genetics of domestication remains unclear for most Amazonian crops. Manioc (Manihot esculenta) is the most important staple food crop that originated in this region. Although manioc is relatively well-studied, little is known about the diversification of bitter and sweet landraces and how they were dispersed across Amazonia. We evaluated single nucleotide polymorphisms (SNPs) in wild and cultivated manioc to identify outlier SNPs putatively under selection and to assess the neutral genetic structure of landraces to make inferences about the evolution of the crop in Amazonia. Some outlier SNPs were in putative manioc genes possibly related to plant architecture, transcriptional regulation and responses to stress. The neutral SNPs revealed contrasting genetic structuring for bitter and sweet landraces. The outlier SNPs may be signatures of the genomic changes resulting from domestication, while the neutral genetic structure suggests independent dispersals for sweet and bitter manioc, possibly related to the earlier domestication and diversification of the former. Our results highlight the role of ancient peoples and current smallholders in the management and conservation of manioc genetic diversity, including putative genes and specific genetic resources with adaptive potential in the context of climate change.

Genome-wide identification, phylogeny, and expression analysis of the SBP-box gene family in Euphorbiaceae

Background

Euphorbiaceae is one of the largest families of flowering plants. Due to its exceptional growth form diversity and near-cosmopolitan distribution, it has attracted much interest since ancient times. SBP-box (SBP) genes encode plant-specific transcription factors that play critical roles in numerous biological processes, especially flower development. We performed genome-wide identification and characterization of SBP genes from four economically important Euphorbiaceae species.

Results

In total, 77 SBP genes were identified in four Euphorbiaceae genomes. The SBP proteins were divided into three length ranges and 10 groups. Group-6 was absent in Arabidopsis thaliana but conserved in Euphorbiaceae. Segmental duplication played the most important role in the expansion processes of Euphorbiaceae SBP genes, and all the duplicated genes were subjected to purify selection. In addition, about two-thirds of the Euphorbiaceae SBP genes are potential targets of miR156, and some miR-regulated SBP genes exhibited high intensity expression and differential expression in different tissues. The expression profiles related to different stress treatments demonstrated broad involvement of Euphorbiaceae SBP genes in response to various abiotic factors and hormonal treatments.

Conclusions

In this study, 77 SBP genes were identified in four Euphorbiaceae species, and their phylogenetic relationships, protein physicochemical characteristics, duplication, tissue and stress response expression, and potential roles in Euphorbiaceae development were studied. This study lays a foundation for further studies of Euphorbiaceae SBP genes, providing valuable information for future functional exploration of Euphorbiaceae SBP genes.

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.

Hybrid de novo transcriptome assembly of poinsettia (Euphorbia pulcherrima Willd. Ex Klotsch) bracts

BACKGROUND

Poinsettia is a popular and important ornamental crop, mostly during the Christmas season. Its bract coloration ranges from pink/red to creamy/white shades. Despite its ornamental value, there is a lack of knowledge about the genetics and molecular biology of poinsettia, especially on the mechanisms of color formation. We performed an RNA-Seq analysis in order to shed light on the transcriptome of poinsettia bracts. Moreover, we analyzed the transcriptome differences of red- and white-bracted poinsettia varieties during bract development and coloration. For the assembly of a bract transcriptome, two paired-end cDNA libraries from a red and white poinsettia pair were sequenced with the Illumina technology, and one library from a red-bracted variety was used for PacBio sequencing. Both short and long reads were assembled using a hybrid de novo strategy. Samples of red- and white-bracted poinsettias were sequenced and comparatively analyzed in three color developmental stages in order to understand the mechanisms of color formation and accumulation in the species.

RESULTS

The final transcriptome contains 288,524 contigs, with 33% showing confident protein annotation against the TAIR10 database. The BUSCO pipeline, which is based on near-universal orthologous gene groups, was applied to assess the transcriptome completeness. From a total of 1440 BUSCO groups searched, 77% were categorized as complete (41% as single-copy and 36% as duplicated), 10% as fragmented and 13% as missing BUSCOs. The gene expression comparison between red and white varieties of poinsettia showed a differential regulation of the flavonoid biosynthesis pathway only at particular stages of bract development. An initial impairment of the flavonoid pathway early in the color accumulation process for the white poinsettia variety was observed, but these differences were no longer present in the subsequent stages of bract development. Nonetheless, GSTF11 and UGT79B10 showed a lower expression in the last stage of bract development for the white variety and, therefore, are potential candidates for further studies on poinsettia coloration.

CONCLUSIONS

In summary, this transcriptome analysis provides a valuable foundation for further studies on poinsettia, such as plant breeding and genetics, and highlights crucial information on the molecular mechanism of color formation.

Density threshold-based acaricide application for the two-spotted spider mite Tetranychus urticae on cassava: from laboratory to the field

BACKGROUD

The red spider mite Tetranychus urticae is known as a worldwide pest in cassava production which can caused serious economic losses. Because no threshold levels are established for T. urticae in cassava cropping, many growers in China are applying excessive chemical applications per cassava cropping period. This is not only expensive, but also results in lower yields because of severe leaf damage and development of resistance. This study aims to develop an immediate threshold level for T. urticae control which could be determined in the laboratory and fit the field application requirements.

RESULTS

The mite density of 25 mites/leaf was the threshold which caused significant decrease of photosynthetic pigment contents and protective enzyme activities in cassava leaves in the laboratory. Moreover, 25 mites/leaf was also the threshold density which resulted in significantly lower leaf damage and higher yield levels compared with a higher level of mites/leaf where calendar sprays were used.

CONCLUSION

The mite density threshold that brought about significant physiological and biochemical changes in the laboratory basically coincided with the threshold that resulted in significant yield loss in the field. It is therefore concluded that the optimum threshold acaricide spray for T. urticae on cassava is 25 mites/leaf. © 2019 Society of Chemical Industry.

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.

Genomic selection for productive traits in biparental cassava breeding populations

Cassava improvement using traditional breeding strategies is slow due to the species’ long breeding cycle. However, the use of genomic selection can lead to a shorter breeding cycle. This study aimed to estimate genetic parameters for productive traits based on pedigree (pedigree and phenotypic information) and genomic (markers and phenotypic information) analyses using biparental crosses at different stages of selection. A total of 290 clones were genotyped and phenotyped for fresh root yield (FRY), dry matter content (DMC), dry yield (DY), fresh shoot yield (FSY) and harvest index (HI). The clones were evaluated in clonal evaluation trials (CET), preliminary yield trials (PYT), advanced yield trials (AYT) and uniform yield trials (UYT), from 2013 to 2018 in ten locations. The breeding stages were analyzed as follows: one stage (CET), two stages (CET and PYT), three stages (CET, PYT and AYT) and four stages (CET, PYT, AYT and UYT). The genomic predictions were analyzed via k-fold cross-validation based on the genomic best linear unbiased prediction (GBLUP) considering a model with genetic additive effects and genotype × location interactions. Genomic and pedigree accuracies were moderate to high (0.56–0.72 and 0.62–0.78, respectively) for important starch-related traits such as DY and FRY; when considering one breeding stage (CET) with the aim of early selection, the genomic accuracies ranged from 0.60 (DMC) to 0.71 (HI). Moreover, the correlations between the genomic estimation breeding values of one-stage genomic analysis and the estimated breeding values of the four-stage (full data set) pedigree analysis were high for all traits as well as for a selection index including all traits. The results indicate great possibilities for genomic selection in cassava, especially for selection early in the breeding cycle (saving time and effort).

Reproductive biology in cassava: stigma receptivity and pollen tube growth

Knowledge on the reproductive biology of cassava, relevant to breeders and molecular geneticists, is still limited. Therefore, different studies were carried out to determine the duration of stigma receptivity and the rate of pollen tube growth. Inflorescences were covered for up to 3 days after the first opening of the bracts (e.g. anthesis day) to prevent open pollination. Results indicate that fruit and seed set are drastically reduced when flowers were covered for 2 or 3 days. However, fruits and seeds were obtained even from flowers that had been covered for 3 days after anthesis, although at low frequency. The rate of pollen tube growth was assessed in many combinations of female and male progenitors crossed through controlled pollinations and collecting the pistils at varying hours after pollination (HAP). Pollen tube growth is fast during the first 6 HAP reaching the tip of the nucellar beak. The growth slows down thereafter, taking 10 additional hours to reach the end of the beak. The growth of pollen tubes slows down even further until they enter the embryo sac. Only 10% of samples showed pollen tubes entering the embryo sac between 48 and 66 HAP. Although several tubes may reach the nucellar beak, only one was observed entering the embryo sac. Results, across the different experiments, were highly variable suggesting that the timeline of fertilization is influenced both by genotypic and environmental factors as well as the manual manipulation of inflorescences and cyathia.

Phylogenetic approaches to identifying fragments of the same gene, with application to the wheat genome

MOTIVATION

As the time and cost of sequencing decrease, the number of available genomes and transcriptomes rapidly increases. Yet the quality of the assemblies and the gene annotations varies considerably and often remains poor, affecting downstream analyses. This is particularly true when fragments of the same gene are annotated as distinct genes, which may cause them to be mistaken as paralogs.

RESULTS

In this study, we introduce two novel phylogenetic tests to infer non-overlapping or partially overlapping genes that are in fact parts of the same gene. One approach collapses branches with low bootstrap support and the other computes a likelihood ratio test. We extensively validated these methods by (i) introducing and recovering fragmentation on the bread wheat, Triticum aestivum cv. Chinese Spring, chromosome 3B; (ii) by applying the methods to the low-quality 3B assembly and validating predictions against the high-quality 3B assembly; and (iii) by comparing the performance of the proposed methods to the performance of existing methods, namely Ensembl Compara and ESPRIT. Application of this combination to a draft shotgun assembly of the entire bread wheat genome revealed 1221 pairs of genes that are highly likely to be fragments of the same gene. Our approach demonstrates the power of fine-grained evolutionary inferences across multiple species to improving genome assemblies and annotations.

AVAILABILITY AND IMPLEMENTATION

An open source software tool is available at https://github.com/DessimozLab/esprit2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.