An EST-derived SNP and SSR genetic linkage map of cassava (Manihot esculenta Crantz)

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

KEY MESSAGE

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.

KASPar SNP genetic map of cassava for QTL discovery of productivity traits in moderate drought stress environment in Africa

Cassava is an important staple in Sub-Sahara Africa. While its production has rapidly expanded to the dry savannahs of the continent, productivity is low in this ecology due to drought by farmers, extending the growth cycle from 12 months to 18, and sometimes 24 months to ensure better harvests. Yield is a complex trait and often difficult to manipulate for genetic gain in conventional breeding. Unfortunately, the dearth of molecular tools for decades has hampered molecular breeding (MB) to improve cassava productivity. This study was conducted to explore KASpar SNPs to generate more molecular tools to enhance genetic dissection of elite African germplasm for improved cassava productivity in dry environments of Africa where molecular resources are highly limited for crop improvement. To aid molecular genetic analysis of traits, a linkage map covering 1582.8 cM with an average resolution of 3.69 cM was constructed using 505 polymorphic SNP markers distributed over 21 linkage groups. Composite interval mapping using 267 F1 progeny in initial QTL mapping identified 27 QTLs for productivity traits in the dry savannah of Nigeria. The availability of KASPar SNPs are anticipated to improve the implementation of MB for the development of high performing drought-tolerant cassava varieties in Africa.

Technological Innovations for Improving Cassava Production in Sub-Saharan Africa

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

Whole Genome Resequencing of Capsicum baccatum and Capsicum annuum to Discover Single Nucleotide Polymorphism Related to Powdery Mildew Resistance

The present study deals with genome wide identification of single-nucleotide polymorphism (SNP) markers related to powdery mildew (PM) resistance in two pepper varieties. Capsicum baccatum (PRH1- a PM resistant line) and Capsicum annuum (Saengryeg- a PM susceptible line), were resequenced to develop SNP markers. A total of 6,213,009 and 6,840,889 SNPs for PRH1 and Saengryeg respectively have been discovered. Among the SNPs, majority were classified as homozygous type SNPs, particularly in the resistant line. Moreover, the SNPs were differentially distributed among the chromosomes in both the resistant and susceptible lines. In total, 4,887,031 polymorphic SNP loci were identified between the two lines and 306,871 high-resolution melting (HRM) marker primer sets were designed. In order to understand the SNPs associated with the vital genes involved in diseases resistance and stress associated processes, chromosome-wise gene ontology analysis was performed. The results revealed the occurrence that SNPs related to diseases resistance genes were predominantly distributed in chromosome 4. In addition, 6281 SNPs associated with 46 resistance genes were identified. Among the lines, PRH1 consisted of maximum number of polymorphic SNPs related to NBS-LRR genes. The SNP markers were validated using HRM assay in 45 F4 populations and correlated with the phenotypic disease index.

QTL associated with resistance to cassava brown streak and cassava mosaic diseases in a bi-parental cross of two Tanzanian farmer varieties, Namikonga and Albert

KEY MESSAGE

QTL consistent across seasons were detected for resistance to cassava brown streak disease induced root necrosis and foliar symptoms. The CMD2 locus was detected in an East African landrace, and comprised two QTL. Cassava production in Africa is compromised by cassava brown streak disease (CBSD) and cassava mosaic disease (CMD). To reduce costs and increase the precision of resistance breeding, a QTL study was conducted to identify molecular markers linked to resistance against these diseases. A bi-parental F1 mapping population was developed from a cross between the Tanzanian farmer varieties, Namikonga and Albert. A one-step genetic linkage map comprising 943 SNP markers and 18 linkage groups spanning 1776.2 cM was generated. Phenotypic data from 240 F1 progeny were obtained from two disease hotspots in Tanzania, over two successive seasons, 2013 and 2014. Two consistent QTLs linked to resistance to CBSD-induced root necrosis were identified in Namikonga on chromosomes II (qCBSDRNFc2Nm) and XI (qCBSDRNc11Nm) and a putative QTL on chromosome XVIII (qCBSDRNc18Nm). qCBSDRNFc2Nm was identified at Naliendele in both seasons. The same QTL was also associated with CBSD foliar resistance. qCBSDRNc11Nm was identified at Chambezi in both seasons, and was characterized by three peaks, spanning a distance of 253 kb. Twenty-seven genes were identified within this region including two LRR proteins and a signal recognition particle. In addition, two highly significant CMD resistance QTL (qCMDc12.1A and qCMDc12.2A) were detected in Albert, on chromosome 12. Both qCMDc12.1A and qCMDc12.2A lay within the range of markers reported earlier, defining the CMD2 locus. This is the first time that two loci have been identified within the CMD2 QTL, and in germplasm of apparent East African origin. Additional QTLs with minor effects on CBSD and CMD resistance were also identified.

QTL Mapping for Pest and Disease Resistance in Cassava and Coincidence of Some QTL with Introgression Regions Derived from Manihot glaziovii

Genetic mapping of quantitative trait loci (QTL) for resistance to cassava brown streak disease (CBSD), cassava mosaic disease (CMD), and cassava green mite (CGM) was performed using an F1 cross developed between the Tanzanian landrace, Kiroba, and a breeding line, AR37-80. The population was evaluated for two consecutive years in two sites in Tanzania. A genetic linkage map was derived from 106 F1 progeny and 1,974 SNP markers and spanned 18 chromosomes covering a distance of 1,698 cM. Fifteen significant QTL were identified; two are associated with CBSD root necrosis only, and were detected on chromosomes V and XII, while seven were associated with CBSD foliar symptoms only and were detected on chromosomes IV, VI, XVII, and XVIII. QTL on chromosomes 11 and 15 were associated with both CBSD foliar and root necrosis symptoms. Two QTL were found to be associated with CMD and were detected on chromosomes XII and XIV, while two were associated with CGM and were identified on chromosomes V and X. There are large Manihot glaziovii introgression regions in Kiroba on chromosomes I, XVII, and XVIII. The introgression segments on chromosomes XVII and XVIII overlap with QTL associated with CBSD foliar symptoms. The introgression region on chromosome I is of a different haplotype to the characteristic "Amani haplotype" found in the landrace Namikonga and others, and unlike some other genotypes, Kiroba does not have a large introgression block on chromosome IV. Kiroba is closely related to a sampled Tanzanian "tree cassava." This supports the observation that some of the QTL associated with CBSD resistance in Kiroba are different to those observed in another variety, Namikonga.

Molecular diversity analysis, drought related marker-traits association mapping and discovery of excellent alleles for 100-day old plants by EST-SSRs in cassava germplasms (Manihot esculenta Cranz)

Cassava is the third largest food crop of the world and has strong ability of drought tolerance. In order to evaluate the molecular diversity and to discover novel alleles for drought tolerance in cassava germplasms, we examined a total of 107 abiotic stress related expressed sequence tags-simple sequence repeat (EST-SSR) markers in 134 cassava genotypes coming from planting regions worldwide and performed drought related marker-traits association mapping. As results, we successfully amplified 98 of 107 markers in 97 polymorphic loci and 279 alleles, with 2.87 alleles per locus, gene diversity of 0.48 and polymorphic information content (PIC) of 0.41 on average. The genetic coefficient between every two lines was 0.37 on average, ranging from 0.21 to 0.82. According to our population structure analysis, these samples could be divided into three sub-populations showing obvious gene flow between them. We also performed water stress experiments using 100-day old cassava plants in two years and calculated the drought tolerance coefficients (DTCs) and used them as phenotypes for marker-trait association mapping. We found that 53 markers were significantly associated with these drought-related traits, with a contribution rate for trait variation of 8.60% on average, ranging between 2.66 and 28.09%. Twenty-four of these 53 associated genes showed differential transcription or protein levels which were confirmed by qRT-PCR under drought stress when compared to the control conditions in cassava. Twelve of twenty-four genes were the same differential expression patterns in omics data and results of qRT-PCR. Out of 33 marker-traits combinations on 24 loci, 34 were positive and 53 negative alleles according to their phenotypic effects and we also obtained the typical materials which carried these elite alleles. We also found 23 positive average allele effects while 10 loci were negative according to their allele effects (AAEs). Our results on molecular diversity, locus association and differential expression under drought can prove beneficial to select excellent materials through marker assisted selection and for functional genes research in the future.

Cassava brown streak disease: a threat to food security in Africa

Cassava brown streak disease (CBSD) has emerged as the most important viral disease of cassava (Manihot esculenta) in Africa and is a major threat to food security. CBSD is caused by two distinct species of ipomoviruses, Cassava brown streak virus and Ugandan cassava brown streak virus, belonging to the family Potyviridae. Previously, CBSD was reported only from the coastal lowlands of East Africa, but recently it has begun to spread as an epidemic throughout the Great Lakes region of East and Central Africa. This new spread represents a major threat to the cassava-growing regions of West Africa. CBSD-resistant cassava cultivars are being developed through breeding, and transgenic RNA interference-derived field resistance to CBSD has also been demonstrated. This review aims to provide a summary of the most important studies on the aetiology, epidemiology and control of CBSD and to highlight key research areas that need prioritization.

Conventional breeding, marker-assisted selection, genomic selection and inbreeding in clonally propagated crops: a case study for cassava

Consolidates relevant molecular and phenotypic information on cassava to demonstrate relevance of heterosis, and alternatives to exploit it by integrating different tools. Ideas are useful to other asexually reproduced crops. Asexually propagated crops offer the advantage that all genetic effects can be exploited in farmers' production fields. However, non-additive effects complicate selection because, while influencing the performance of the materials under evaluation, they cannot be transmitted efficiently to the following cycle of selection. Cassava can be used as a model crop for asexually propagated crops because of its diploid nature and the absence of (known) incompatibility effects. New technologies such as genomic selection (GS), use of inbred progenitors based on doubled haploids and induction of flowering can be employed for accelerating genetic gains in cassava. Available information suggests that heterosis, non-additive genetic effects and within-family variation are relatively large for complex traits such as fresh root yield, moderate for dry matter or starch content in the roots, and low for defensive traits (pest and disease resistance) and plant architecture. The present article considers the potential impact of different technologies for maximizing gains for key traits in cassava, and highlights the advantages of integrating them. Exploiting heterosis would be optimized through the implementation of reciprocal recurrent selection. The advantages of using inbred progenitors would allow shifting the current cassava phenotypic recurrent selection method into line improvement, which in turn would allow designing outstanding hybrids rather than finding them by trial and error.

A genetic map of cassava (Manihot esculenta Crantz) with integrated physical mapping of immunity-related genes

BACKGROUND

Cassava, Manihot esculenta Crantz, is one of the most important crops world-wide representing the staple security for more than one billion of people. The development of dense genetic and physical maps, as the basis for implementing genetic and molecular approaches to accelerate the rate of genetic gains in breeding program represents a significant challenge. A reference genome sequence for cassava has been made recently available and community efforts are underway for improving its quality. Cassava is threatened by several pathogens, but the mechanisms of defense are far from being understood. Besides, there has been a lack of information about the number of genes related to immunity as well as their distribution and genomic organization in the cassava genome.

RESULTS

A high dense genetic map of cassava containing 2,141 SNPs has been constructed. Eighteen linkage groups were resolved with an overall size of 2,571 cM and an average distance of 1.26 cM between markers. More than half of mapped SNPs (57.4%) are located in coding sequences. Physical mapping of scaffolds of cassava whole genome sequence draft using the mapped markers as anchors resulted in the orientation of 687 scaffolds covering 45.6% of the genome. One hundred eighty nine new scaffolds are anchored to the genetic cassava map leading to an extension of the present cassava physical map with 30.7 Mb. Comparative analysis using anchor markers showed strong co-linearity to previously reported cassava genetic and physical maps. In silico based searching for conserved domains allowed the annotation of a repertory of 1,061 cassava genes coding for immunity-related proteins (IRPs). Based on physical map of the corresponding sequencing scaffolds, unambiguous genetic localization was possible for 569 IRPs.

CONCLUSIONS

This is the first study reported so far of an integrated high density genetic map using SNPs with integrated genetic and physical localization of newly annotated immunity related genes in cassava. These data build a solid basis for future studies to map and associate markers with single loci or quantitative trait loci for agronomical important traits. The enrichment of the physical map with novel scaffolds is in line with the efforts of the cassava genome sequencing consortium.

Large-scale SNP discovery through RNA sequencing and SNP genotyping by targeted enrichment sequencing in cassava (Manihot esculenta Crantz)

Cassava (Manihot esculenta Crantz) is one of the most important crop species being the main source of dietary energy in several countries. Marker-assisted selection has become an essential tool in plant breeding. Single nucleotide polymorphism (SNP) discovery via transcriptome sequencing is an attractive strategy for genome complexity reduction in organisms with large genomes. We sequenced the transcriptome of 16 cassava accessions using the Illumina HiSeq platform and identified 675,559 EST-derived SNP markers. A subset of those markers was subsequently genotyped by capture-based targeted enrichment sequencing in 100 F1 progeny segregating for starch viscosity phenotypes. A total of 2,110 non-redundant SNP markers were used to construct a genetic map. This map encompasses 1,785 cM and consists of 19 linkage groups. A major quantitative trait locus (QTL) controlling starch pasting properties was identified and shown to coincide with the QTL previously reported for this trait. With a high-density SNP-based linkage map presented here, we also uncovered a novel QTL associated with starch pasting time on LG 10.

High-resolution linkage map and chromosome-scale genome assembly for cassava (Manihot esculenta Crantz) from 10 populations

Cassava (Manihot esculenta Crantz) is a major staple crop in Africa, Asia, and South America, and its starchy roots provide nourishment for 800 million people worldwide. Although native to South America, cassava was brought to Africa 400-500 years ago and is now widely cultivated across sub-Saharan Africa, but it is subject to biotic and abiotic stresses. To assist in the rapid identification of markers for pathogen resistance and crop traits, and to accelerate breeding programs, we generated a framework map for M. esculenta Crantz from reduced representation sequencing [genotyping-by-sequencing (GBS)]. The composite 2412-cM map integrates 10 biparental maps (comprising 3480 meioses) and organizes 22,403 genetic markers on 18 chromosomes, in agreement with the observed karyotype. We used the map to anchor 71.9% of the draft genome assembly and 90.7% of the predicted protein-coding genes. The chromosome-anchored genome sequence will be useful for breeding improvement by assisting in the rapid identification of markers linked to important traits, and in providing a framework for genomic selection-enhanced breeding of this important crop.

AFSM sequencing approach: a simple and rapid method for genome-wide SNP and methylation site discovery and genetic mapping

We describe methods for the assessment of amplified-fragment single nucleotide polymorphism and methylation (AFSM) sites using a quick and simple molecular marker-assisted breeding strategy based on the use of two restriction enzyme pairs (EcoRI-MspI and EcoRI-HpaII) and a next-generation sequencing platform. Two sets of 85 adapter pairs were developed to concurrently identify SNPs, indels and methylation sites for 85 lines of cassava population in this study. In addition to SNPs and indels, the simplicity of the AFSM protocol makes it particularly suitable for high-throughput full methylation and hemi-methylation analyses. To further demonstrate the ease of this approach, a cassava genetic linkage map was constructed. This approach should be widely applicable for genetic mapping in a variety of organisms and will improve the application of crop genomics in assisted breeding.

Transcriptional analysis of South African cassava mosaic virus-infected susceptible and tolerant landraces of cassava highlights differences in resistance, basal defense and cell wall associated genes during infection

BACKGROUND

Cassava mosaic disease is caused by several distinct geminivirus species, including South African cassava mosaic virus-[South Africa:99] (SACMV). To date, there is limited gene regulation information on viral stress responses in cassava, and global transcriptome profiling in SACMV-infected cassava represents an important step towards understanding natural host responses to plant geminiviruses.

RESULTS

A RNA-seq time course (12, 32 and 67 dpi) study, monitoring gene expression in SACMV-challenged susceptible (T200) and tolerant (TME3) cassava landraces, was performed using the Applied Biosystems (ABI) SOLiD next-generation sequencing platform. The multiplexed paired end sequencing run produced a total of 523 MB and 693 MB of paired-end reads for SACMV-infected susceptible and tolerant cDNA libraries, respectively. Of these, approximately 50.7% of the T200 reads and 55.06% of TME3 reads mapped to the cassava reference genome available in phytozome. Using a log2 fold cut-off (p<0.05), comparative analysis between the six normalized cDNA libraries showed that 4181 and 1008 transcripts in total were differentially expressed in T200 and TME3, respectively, across 12, 32 and 67 days post infection, compared to mock-inoculated. The number of responsive transcripts increased dramatically from 12 to 32 dpi in both cultivars, but in contrast, in T200 the levels did not change significantly at 67 dpi, while in TME3 they declined. GOslim functional groups illustrated that differentially expressed genes in T200 and TME3 were overrepresented in the cellular component category for stress-related genes, plasma membrane and nucleus. Alterations in the expression of other interesting genes such as transcription factors, resistance (R) genes, and histone/DNA methylation-associated genes, were observed. KEGG pathway analysis uncovered important altered metabolic pathways, including phenylpropanoid biosynthesis, sucrose and starch metabolism, and plant hormone signalling.

CONCLUSIONS

Molecular mechanisms for TME3 tolerance are proposed, and differences in patterns and levels of transcriptome profiling between T200 and TME3 with susceptible and tolerant phenotypes, respectively, support the hypothesis that viruses rearrange their molecular interactions in adapting to hosts with different genetic backgrounds.

Cassava genome from a wild ancestor to cultivated varieties

Cassava is a major tropical food crop in the Euphorbiaceae family that has high carbohydrate production potential and adaptability to diverse environments. Here we present the draft genome sequences of a wild ancestor and a domesticated variety of cassava and comparative analyses with a partial inbred line. We identify 1,584 and 1,678 gene models specific to the wild and domesticated varieties, respectively, and discover high heterozygosity and millions of single-nucleotide variations. Our analyses reveal that genes involved in photosynthesis, starch accumulation and abiotic stresses have been positively selected, whereas those involved in cell wall biosynthesis and secondary metabolism, including cyanogenic glucoside formation, have been negatively selected in the cultivated varieties, reflecting the result of natural selection and domestication. Differences in microRNA genes and retrotransposon regulation could partly explain an increased carbon flux towards starch accumulation and reduced cyanogenic glucoside accumulation in domesticated cassava. These results may contribute to genetic improvement of cassava through better understanding of its biology.

Construction of a SSR-based genetic map and identification of QTLs for catechins content in tea plant (Camellia sinensis)

Catechins are the most important bioactive compounds in tea, and have been demonstrated to possess a wide variety of pharmacological activities. To characterize quantitative trait loci (QTLs) for catechins content in the tender shoots of tea plant, we constructed a moderately saturated genetic map using 406 simple sequence repeat (SSR) markers, based on a pseudo-testcross population of 183 individuals derived from an intraspecific cross of two Camellia sinensis varieties with diverse catechins composition. The map consisted of fifteen linkage groups (LGs), corresponding to the haploid chromosome number of tea plant (2n = 2x = 30). The total map length was 1,143.5 cM, with an average locus spacing of 2.9 cM. A total of 25 QTLs associated with catechins content were identified over two measurement years. Of these, nine stable QTLs were validated across years, and clustered into four main chromosome regions on LG03, LG11, LG12 and LG15. The population variability explained by each QTL was predominantly at moderate-to-high levels and ranged from 2.4% to 71.0%, with an average of 17.7%. The total number of QTL for each trait varied from four to eight, while the total population variability explained by all QTLs for a trait ranged between 38.4% and 79.7%. This is the first report on the identification of QTL for catechins content in tea plant. The results of this study provide a foundation for further cloning and functional characterization of catechin QTLs for utilization in improvement of tea plant.

Genome-wide discovery and information resource development of DNA polymorphisms in cassava

Cassava (Manihot esculenta Crantz) is an important crop that provides food security and income generation in many tropical countries, and is known for its adaptability to various environmental conditions. Its draft genome sequence and many expressed sequence tags are now publicly available, allowing the development of cassava polymorphism information. Here, we describe the genome-wide discovery of cassava DNA polymorphisms. Using the alignment of predicted transcribed sequences from the cassava draft genome sequence and ESTs from GenBank, we discovered 10,546 single-nucleotide polymorphisms and 647 insertions and deletions. To facilitate molecular marker development for cassava, we designed 9,316 PCR primer pairs to amplify the genomic region around each DNA polymorphism. Of the discovered SNPs, 62.7% occurred in protein-coding regions. Disease-resistance genes were found to have a significantly higher ratio of nonsynonymous-to-synonymous substitutions. We identified 24 read-through (changes of a stop codon to a coding codon) and 38 premature stop (changes of a coding codon to a stop codon) single-nucleotide polymorphisms, and found that the 5 gene ontology terms in biological process were significantly different in genes with read-through single-nucleotide polymorphisms compared with all cassava genes. All data on the discovered DNA polymorphisms were organized into the Cassava Online Archive database, which is available at http://cassava.psc.riken.jp/.