Spectrum and Density of Gamma and X-ray Induced Mutations in a Non-Model Rice Cultivar.

Physical mutagens are a powerful tool used for genetic research and breeding for over eight decades. Yet, when compared to chemical mutagens, data sets on the effect of different mutagens and dosages on the spectrum and density of induced mutations remain lacking. To address this, we investigated the landscape of mutations induced by gamma and X-ray radiation in the most widely cultivated crop species: rice. A mutant population of a tropical upland rice, Oryza sativa L., was generated and propagated via self-fertilization for seven generations. Five dosages ranging from 75 Gy to 600 Gy in both X-ray and gamma-irradiated material were applied. In the process of a forward genetic screens, 11 unique rice mutant lines showing phenotypic variation were selected for mutation analysis via whole-genome sequencing. Thousands of candidate mutations were recovered in each mutant with single base substitutions being the most common, followed by small indels and structural variants. Higher dosages resulted in a higher accumulation of mutations in gamma-irradiated material, but not in X-ray-treated plants. The in vivo role of all annotated rice genes is yet to be directly investigated. The ability to induce a high density of single nucleotide and structural variants through mutagenesis will likely remain an important approach for functional genomics and breeding.

Fragmentation of Pooled PCR Products for Highly Multiplexed TILLING.

Improvements to massively parallel sequencing have allowed the routine recovery of natural and induced sequence variants. A broad range of biological disciplines have benefited from this, ranging from plant breeding to cancer research. The need for high sequence coverage to accurately recover single nucleotide variants and small insertions and deletions limits the applicability of whole genome approaches. This is especially true in organisms with a large genome size or for applications requiring the screening of thousands of individuals, such as the reverse-genetic technique known as TILLING. Using PCR to target and sequence chosen genomic regions provides an attractive alternative as the vast reduction in interrogated bases means that sample size can be dramatically increased through amplicon multiplexing and multi-dimensional sample pooling while maintaining suitable coverage for recovery of small mutations. Direct sequencing of PCR products is limited, however, due to limitations in read lengths of many next generation sequencers. In the present study we show the optimization and use of ultrasonication for the simultaneous fragmentation of multiplexed PCR amplicons for TILLING highly pooled samples. Sequencing performance was evaluated in a total of 32 pooled PCR products produced from 4096 chemically mutagenized Hordeum vulgare DNAs pooled in three dimensions. Evaluation of read coverage and base quality across amplicons suggests this approach is suitable for high-throughput TILLING and other applications employing highly pooled complex sampling schemes. Induced mutations previously identified in a traditional TILLING screen were recovered in this dataset further supporting the efficacy of the approach.

Induction and recovery of copy number variation in banana through gamma irradiation and low-coverage whole-genome sequencing.

Traditional breeding methods are hindered in bananas due to the fact that major cultivars are sterile, parthenocarpic, triploid and thus clonally propagated. This has resulted in a narrow genetic base and limited resilience to biotic and abiotic stresses. Mutagenesis of in vitro propagated bananas is one method to introduce novel alleles and broaden genetic diversity. We previously established a method for the induction and recovery of single nucleotide mutations generated with the chemical mutagen EMS. However, officially released mutant banana varieties have been created using gamma rays, a mutagen that can produce large genomic insertions and deletions (indels). Such dosage mutations may be important for generating observable phenotypes in polyploids. In this study, we establish a low-coverage whole-genome sequencing approach in triploid bananas to recover large genomic indels caused by treatment with gamma irradiation. We first evaluated the commercially released mutant cultivar 'Novaria' and found that it harbours multiple predicted deletions, ranging from 0.3 to 3.8 million base pairs (Mbp). In total, predicted deletions span 189 coding regions. To evaluate the feasibility of generating and maintaining new mutations, we developed a pipeline for mutagenesis and screening for copy number variation in Cavendish bananas using the cultivar 'Williams'. Putative mutations were recovered in 70% of lines treated with 20 Gy and 60% of the lines treated with 40 Gy. While deletion events predominate, insertions were identified in 20 Gy-treated material. Based on these results, we believe this approach can be scaled up to support large breeding projects.

Efficient CRISPR/Cas9 Genome Editing of Phytoene desaturase in Cassava.

CRISPR/Cas9 has become a powerful genome-editing tool for introducing genetic changes into crop species. In order to develop capacity for CRISPR/Cas9 technology in the tropical staple cassava (Manihot esculenta), the Phytoene desaturase (MePDS) gene was targeted in two cultivars using constructs carrying gRNAs targeting two sequences within MePDS exon 13. After Agrobacterium-mediated delivery of CRISPR/Cas9 reagents into cassava cells, both constructs induced visible albino phenotypes within cotyledon-stage somatic embryos regenerating on selection medium and the plants regenerated therefrom. A total of 58 (cv. 60444) and 25 (cv. TME 204) plant lines were recovered, of which 38 plant lines (19 from each cultivar) were analyzed for mutagenesis. The frequency of plant lines showing albino phenotype was high, ranging from 90 to 100% in cv. TME 204. Observed albino phenotypes were comprised of full albinos devoid of green tissue and chimeras containing a mixture of white and green tissues. Sequence analysis revealed that 38/38 (100%) of the plant lines examined carried mutations at the targeted MePDS site, with insertions, deletions, and substitutions recorded. One putatively mono-allelic homozygous line (1/19) was found from cv. 60444, while 1 (1/19) and 4 (4/19) putatively bi-allelic homozygous lines were found in 60444 and TME204, respectively. The remaining plant lines, comprised mostly of the chimeras, were found to be putatively heterozygous. We observed minor (1 bp) nucleotide substitutions and or deletions upstream of the 5' and or downstream of the 3' targeted MePDS region. The data reported demonstrates that CRISPR/Cas9-mediated genome editing of cassava is highly efficient and relatively simple, generating multi-allelic mutations in both cultivars studied. Modification of MePDS described here generates visually detectable mutated events in a relatively short time frame of 6-8 weeks, and does not require sequencing to confirm editing at the target. It therefore provides a valuable platform to facilitate rapid assessment and optimization of CRISPR/Cas9 and other genome-editing technologies in cassava.

Gene-based microsatellites for cassava (Manihot esculenta Crantz): prevalence, polymorphisms, and cross-taxa utility.

BACKGROUND: Cassava (Manihot esculenta Crantz), a starchy root crop grown in tropical and subtropical climates, is the sixth most important crop in the world after wheat, rice, maize, potato and barley. The repertoire of simple sequence repeat (SSR) markers for cassava is limited and warrants a need for a larger number of polymorphic SSRs for germplasm characterization and breeding applications. RESULTS: A total of 846 putative microsatellites were identified in silico from an 8,577 cassava unigene set with an average density of one SSR every 7 kb. One hundred and ninety-two candidate SSRs were screened for polymorphism among a panel of cassava cultivars from Africa, Latin America and Asia, four wild Manihot species as well as two other important taxa in the Euphorbiaceae, leafy spurge (Euphorbia esula) and castor bean (Ricinus communis). Of 168 markers with clean amplification products, 124 (73.8%) displayed polymorphism based on high resolution agarose gels. Of 85 EST-SSR markers screened, 80 (94.1%) amplified alleles from one or more wild species (M epruinosa, M glaziovii, M brachyandra, M tripartita) whereas 13 (15.3%) amplified alleles from castor bean and 9 (10.6%) amplified alleles from leafy spurge; hence nearly all markers were transferable to wild relatives of M esculenta while only a fraction was transferable to the more distantly related taxa. In a subset of 20 EST-SSRs assessed by fluorescence-based genotyping the number of alleles per locus ranged from 2 to 10 with an average of 4.55 per locus. These markers had a polymorphism information content (PIC) from 0.19 to 0.75 with an average value of 0.55 and showed genetic relationships consistent with existing information on these genotypes. CONCLUSION: A set of 124 new, unique polymorphic EST-SSRs was developed and characterized which extends the repertoire of SSR markers for cultivated cassava and its wild relatives. The markers show high PIC values and therefore will be useful for cultivar identification, taxonomic studies, and genetic mapping. The study further shows that mining ESTs is a highly efficient strategy for polymorphism detection within the cultivated cassava gene pool.