Cassava (Manihot esculenta) provides calories and nutrition for more than half a billion people. It was domesticated by native Amazonian peoples through cultivation of the wild progenitor M. esculenta ssp. flabellifolia and is now grown in tropical regions worldwide. Here we provide a high-quality genome assembly for cassava with improved contiguity, linkage, and completeness; almost 97% of genes are anchored to chromosomes. We find that paleotetraploidy in cassava is shared with the related rubber tree Hevea, providing a resource for comparative studies. We also sequence a global collection of 58 Manihot accessions, including cultivated and wild cassava accessions and related species such as Ceará or India rubber (M. glaziovii), and genotype 268 African cassava varieties. We find widespread interspecific admixture, and detect the genetic signature of past cassava breeding programs. As a clonally propagated crop, cassava is especially vulnerable to pathogens and abiotic stresses. This genomic resource will inform future genome-enabled breeding efforts to improve this staple crop. 13 International Institute of Tropical Agriculture (IITA), Nairobi, Kenya. 14 Dow AgroSciences, Indianapolis, Indiana, USA. 15 Molecular Genetics Unit, Okinawa Institute of Science and Technology Graduate University, Onna, Japan. 16 In this report we use "cassava" to refer to cultivated and/or domesticated varieties of M. esculenta, and the shorthand M. esc. flabellifolia for wild accessions 3 . We also shotgun-sequenced five Manihot accessions related to cassava, including three from the wild species M. glaziovii Muell. Arg., one named M. pseudoglaziovii Pax & K. Hoffman, and "tree" cassava, a suspected hybrid sometimes called M. catingea Ule 12,18 . The Ceará or India rubber tree species M. glaziovii, also domesticated in South America, was imported to East Africa in the early twentieth century. It is interfertile with cassava and has been used in African breeding programs to exploit the natural resistance of M. glaziovii to cassava pathogens 18 . To analyze genetic variation present in African varieties, we also characterized 268 cultivars of cassava using reduced representation genotypingby-sequencing (GBS) 19 (Table 2).
RESULTS
Chromosome structureTo produce a high-quality chromosome-scale reference genome for cassava, we augmented our earlier draft sequence 20 of the reference genotype AM560-2 with additional whole genome shotgun sequencing and mate pair data, fosmid-end sequences, and a paired-end library developed using proximity ligation of in vitro reconstituted chromatin 21 (Methods and Supplementary Note 1). AM560-2 is an S3 line bred at Centro Internacional de Agricultura Tropical (CIAT) from MCOL1505 (also known as Manihoica P-12 (ref. 22). Compared with the previous draft 23 , the contiguity of our new shotgun assembly has more than doubled (N50 length 27.7 kb vs. 11.5 kb), and an additional 135 Mb is anchored to chromosomes 23 (Supplementary Note 1). To organize the sequence into chromosomes we integrated the shotgun ...