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 ...

Symbiotic root nodules are beneficial to leguminous host plants; however, excessive nodulation damages the host because it interferes with the distribution of nutrients in the plant. To keep a steady balance, the nodulation programme is regulated systemically in leguminous hosts. Leguminous mutants that have lost this ability display a hypernodulating phenotype. Through the use of reciprocal and self-grafting studies using Lotus japonicus hypernodulating mutants, har1 (also known as sym78), we show that the shoot genotype is responsible for the negative regulation of nodule development. A map-based cloning strategy revealed that HAR1 encodes a protein with a relative molecular mass of 108,000, which contains 21 leucine-rich repeats, a single transmembrane domain and serine/threonine kinase domains. The har1 mutant phenotype was rescued by transfection of the HAR1 gene. In a comparison of Arabidopsis receptor-like kinases, HAR1 showed the highest level of similarity with CLAVATA1 (CLV1). CLV1 negatively regulates formation of the shoot and floral meristems through cell-cell communication involving the CLV3 peptide. Identification of hypernodulation genes thus indicates that genes in leguminous plants bearing a close resemblance to CLV1 regulate nodule development systemically, by means of organ-organ communication.

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.