BackgroundDNA methylation is an important feature of plant epigenomes, involved in the formation of heterochromatin and affecting gene expression. Extensive variation of DNA methylation patterns within a species has been uncovered from studies of natural variation. However, the extent to which DNA methylation varies between flowering plant species is still unclear. To understand the variation in genomic patterning of DNA methylation across flowering plant species, we compared single base resolution DNA methylomes of 34 diverse angiosperm species.ResultsBy analyzing whole-genome bisulfite sequencing data in a phylogenetic context, it becomes clear that there is extensive variation throughout angiosperms in gene body DNA methylation, euchromatic silencing of transposons and repeats, as well as silencing of heterochromatic transposons. The Brassicaceae have reduced CHG methylation levels and also reduced or loss of CG gene body methylation. The Poaceae are characterized by a lack or reduction of heterochromatic CHH methylation and enrichment of CHH methylation in genic regions. Furthermore, low levels of CHH methylation are observed in a number of species, especially in clonally propagated species.ConclusionsThese results reveal the extent of variation in DNA methylation in angiosperms and show that DNA methylation patterns are broadly a reflection of the evolutionary and life histories of plant species.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-1059-0) contains supplementary material, which is available to authorized users.

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

To understand the variation in genomic patterning of DNA methylation we compared methylomes of 34 diverse angiosperm species. By analyzing whole-genome bisulfite sequencing data in a phylogenetic context it becomes clear that there is extensive variation throughout angiosperms in gene body DNA methylation, euchromatic silencing of transposons and repeats, as well as silencing of heterochromatic transposons. The Brassicaceae have reduced CHG methylation levels and also reduced or loss of CG gene body methylation. The Poaceae are characterized by a lack or reduction of heterochromatic CHH methylation and enrichment of CHH methylation in genic regions. Reduced CHH methylation levels are found in clonally propagated species, suggesting that these methods of propagation may alter the epigenomic landscape over time. These results show that DNA methylation patterns are broadly a reflection of the evolutionary and life histories of plant species.

Cassava (Manihot esculenta Crantz) is an important staple food crop in Africa and South America; however, ubiquitous deleterious mutations may severely decrease its fitness. To evaluate these deleterious mutations, we constructed a cassava haplotype map through deep sequencing 241 diverse accessions and identified >28 million segregating variants. We found that (i) although domestication has modified starch and ketone metabolism pathways to allow for human consumption, the concomitant bottleneck and clonal propagation have resulted in a large proportion of fixed deleterious amino acid changes, increased the number of deleterious alleles by 26%, and shifted the mutational burden toward common variants; (ii) deleterious mutations have been ineffectively purged, owing to limited recombination in the cassava genome; (iii) recent breeding efforts have maintained yield by masking the most damaging recessive mutations in the heterozygous state but have been unable to purge the mutation burden; such purging should be a key target in future cassava breeding.For millions of people in the tropics, cassava is the third most consumed carbohydrate source, after rice and maize 1 . Even though cassava was domesticated in Latin America 2,3 , it has spread widely and has become a major staple crop in Africa. Although its wild progenitor, M. esculenta sp. falbellifolia, reproduces by seed 4 , cultivated cassava is notably almost exclusively clonally propagated via stem cutting 5 . The limited number of recombination events in such vegetatively propagated crops may result in an accumulation of deleterious alleles throughout the genome 6 . Thus, mutation burden in cassava is expected to be more severe than that in sexually propagated species. Deleterious mutations are considered to be at the heart of inbreeding depression 7 . Even in elite cassava accessions, inbreeding depression is extremely severe, and a single generation of inbreeding may result in a >60% decrease in fresh root yield 8,9 . In this study, we sought to identify deleterious mutations in cassava populations, with the goal of accelerating cassava breeding by allowing breeders to purge deleterious mutations more efficiently.We conducted a comprehensive characterization of genetic variation by whole-genome sequencing (WGS) of 241 cassava accessions ( Fig. 1, Supplementary Fig. 1 and Supplementary Table 1). On average, more than 30× coverage was generated for each accession. To ensure high-quality variant discovery, variants from low-copynumber regions of the cassava genome 10,11 were identified to develop the cassava haplotype map II (HapMapII) (Supplementary Fig. 2), containing 25.9 million SNPs and 1.9 million insertions/deletions (indels) (Supplementary Table 2), of which nearly 50% were rare (minor-allele frequency <0.05) (Supplementary Fig. 3). The error rate of variant calling, i.e., the proportion of segregating sites in the reference accession, was 0.01%. The correlation between read depth and the proportion of SNP heterozygosity was extremely low (r 2 = 6 × 10...

More than 250 million Africans rely on the starchy root crop cassava (Manihot esculenta) as their staple source of calories. A typical cassava-based diet, however, provides less than 30% of the minimum daily requirement for protein and only 10%-20% of that for iron, zinc, and vitamin A. The BioCassava Plus (BC+) program has employed modern biotechnologies intended to improve the health of Africans through the development and delivery of genetically engineered cassava with increased nutrient (zinc, iron, protein, and vitamin A) levels. Additional traits addressed by BioCassava Plus include increased shelf life, reductions in toxic cyanogenic glycosides to safe levels, and resistance to viral disease. The program also provides incentives for the adoption of biofortified cassava. Proof of concept was achieved for each of the target traits. Results from field trials in Puerto Rico, the first confined field trials in Nigeria to use genetically engineered organisms, and ex ante impact analyses support the efficacy of using transgenic strategies for the biofortification of cassava.

Cassava (Manihot esculenta Crantz) is a clonally propagated staple food crop in the tropics. Genomic selection (GS) has been implemented at three breeding institutions in Africa to reduce cycle times. Initial studies provided promising estimates of predictive abilities. Here, we expand on previous analyses by assessing the accuracy of seven prediction models for seven traits in three prediction scenarios: cross-validation within populations, cross-population prediction and cross-generation prediction. We also evaluated the impact of increasing the training population (TP) size by phenotyping progenies selected either at random or with a genetic algorithm. Cross-validation results were mostly consistent across programs, with nonadditive models predicting of 10% better on average. Cross-population accuracy was generally low (mean = 0.18) but prediction of cassava mosaic disease increased up to 57% in one Nigerian population when data from another related population were combined. Accuracy across generations was poorer than within-generation accuracy, as expected, but accuracy for dry matter content and mosaic disease severity should be sufficient for rapid-cycling GS. Selection of a prediction model made some difference across generations, but increasing TP size was more important. With a genetic algorithm, selection of one-third of progeny could achieve an accuracy equivalent to phenotyping all progeny. We are in the early stages of GS for this crop but the results are promising for some traits. General guidelines that are emerging are that TPs need to continue to grow but phenotyping can be done on a cleverly selected subset of individuals, reducing the overall phenotyping burden.

Cassava (Manihot esculenta Crantz) roots spoil 2 to 3 d after harvest because of postharvest physiological deterioration (PPD), which has remained an unsolved problem. Roots from different sources of germplasm were evaluated 5, 10, 20, and 40 d after harvest and some were found to be tolerant to PPD. Three genotypes showed zero levels of PPD even 40 d after harvest. Tolerance to PPD in roots with high carotenoid levels may be explained by their antioxidant properties. Irradiation of seeds in mutagenized populations may have silenced one of the genes involved in the expression of PPD. The tolerance to PPD found in other sources cannot be properly explained. The identification of several sources of resistance (and at least two different modes of action) suggests that now there are alternatives available for solving this problem, benefiting millions of resource‐limited farmers worldwide.

Cassava (Manihot esculenta Crantz) is a crucial, under-researched crop feeding millions worldwide, especially in Africa. Cassava mosaic disease (CMD) has plagued production in Africa for over a century. Biparental mapping studies suggest primarily a single major gene mediates resistance. To investigate this genetic architecture, we conducted the first genome-wide association mapping study in cassava with up to 6128 genotyping-by-sequenced African breeding lines and 42,113 reference genome-mapped single-nucleotide polymorphism (SNP) markers. We found a single region on chromosome 8 that accounts for 30 to 66% of genetic resistance in the African cassava germplasm. Thirteen additional regions with small effects were also identified. Further dissection of the major quantitative trait locus (QTL) on chromosome 8 revealed the presence of two possibly epistatic loci and/or multiple resistance alleles, which may account for the difference between moderate and strong disease resistances in the germplasm. Search of potential candidate genes in the major QTL region identified two peroxidases and one thioredoxin. Finally, we found genomic prediction accuracy of 0.53 to 0.58 suggesting that genomic selection (GS) will be effective both for improving resistance in breeding populations and identifying highly resistant clones as varieties.