Controlling elemental composition is critical for plant growth and development as well as the nutrition of humans who utilize plants for food. Uncovering the genetic architecture underlying mineral ion homeostasis in plants is a critical first step towards understanding the biochemical networks that regulate a plant's elemental composition (ionome). Natural accessions of Arabidopsis thaliana provide a rich source of genetic diversity that leads to phenotypic differences. We analyzed the concentrations of 17 different elements in 12 A. thaliana accessions and three recombinant inbred line (RIL) populations grown in several different environments using high-throughput inductively coupled plasma- mass spectroscopy (ICP-MS). Significant differences were detected between the accessions for most elements and we identified over a hundred QTLs for elemental accumulation in the RIL populations. Altering the environment the plants were grown in had a strong effect on the correlations between different elements and the QTLs controlling elemental accumulation. All ionomic data presented is publicly available at www.ionomicshub.org.

Whole genome sequencing has allowed rapid progress in the application of forward genetics in model species. In this study, we demonstrated an application of next-generation sequencing for forward genetics in a complex crop genome. We sequenced an ethyl methanesulfonate-induced mutant of Sorghum bicolor defective in hydrogen cyanide release and identified the causal mutation. A workflow identified the causal polymorphism relative to the reference BTx623 genome by integrating data from single nucleotide polymorphism identification, prior information about candidate gene(s) implicated in cyanogenesis, mutation spectra, and polymorphisms likely to affect phenotypic changes. A point mutation resulting in a premature stop codon in the coding sequence of dhurrinase2, which encodes a protein involved in the dhurrin catabolic pathway, was responsible for the acyanogenic phenotype. Cyanogenic glucosides are not cyanogenic compounds but their cyanohydrins derivatives do release cyanide. The mutant accumulated the glucoside, dhurrin, but failed to efficiently release cyanide upon tissue disruption. Thus, we tested the effects of cyanide release on insect herbivory in a genetic background in which accumulation of cyanogenic glucoside is unchanged. Insect preference choice experiments and herbivory measurements demonstrate a deterrent effect of cyanide release capacity, even in the presence of wild-type levels of cyanogenic glucoside accumulation. Our gene cloning method substantiates the value of (1) a sequenced genome, (2) a strongly penetrant and easily measurable phenotype, and (3) a workflow to pinpoint a causal mutation in crop genomes and accelerate in the discovery of gene function in the postgenomic era.

In order to leverage novel sequencing techniques for cloning genes in eukaryotic organisms with complex genomes, the false positive rate of variant discovery must be controlled for by experimental design and informatics. We sequenced five lines from three pedigrees of ethyl methanesulfonate (EMS)-mutagenized Sorghum bicolor, including a pedigree segregating a recessive dwarf mutant. Comparing the sequences of the lines, we were able to identify and eliminate error-prone positions. One genomic region contained EMS mutant alleles in dwarfs that were homozygous reference sequences in wild-type siblings and heterozygous in segregating families. This region contained a single nonsynonymous change that cosegregated with dwarfism in a validation population and caused a premature stop codon in the Sorghum ortholog encoding the gibberellic acid (GA) biosynthetic enzyme ent-kaurene oxidase. Application of exogenous GA rescued the mutant phenotype. Our method for mapping did not require outcrossing and introduced no segregation variance. This enables work when line crossing is complicated by life history, permitting gene discovery outside of genetic models. This inverts the historical approach of first using recombination to define a locus and then sequencing genes. Our formally identical approach first sequences all the genes and then seeks cosegregation with the trait. Mutagenized lines lacking obvious phenotypic alterations are available for an extension of this approach: mapping with a known marker set in a line that is phenotypically identical to starting material for EMS mutant generation.

Leaf architecture determines plant structural integrity, light harvesting, and economic considerations such as plant density. Ligules, junctions at the leaf sheath and blade in grasses, protect stalks from environmental stresses and, in conjunction with auricles, controls leaf angle. Previous studies in mutants have recessive liguleless mutants (lg1 and lg2) and dominant mutations in knotted1-like homeobox genes (Lg3-O, Lg4, and Kn1) involved in ligule development. Recently, a new semidominant liguleless mutant, Liguleless narrow (Lgn-R), has been characterized in maize that affects ligule and auricle development and results in a narrow leaf phenotype. We show that quantitative genetic variation affects penetrance of Lgn-R. To examine the genetic architecture underlying Lgn-R expressivity, crosses between Lgn-R/+ mutants in a B73 background and intermated B73 x Mo17 recombinant inbred lines were evaluated in multiple years and locations. A single main-effect quantitative trait locus (QTL) on chromosome 1 (sympathy for the ligule; sol) was discovered with a Mo17-contributed allele that suppressed Lgn-R mutant phenotypes. This QTL has a genetic-interaction with a locus on chromosome 7 (lucifer; lcf) for which the B73-contributed allele increases the ability of the solMo17 allele to suppress Lgn-R. Neither of the genetic intervals likely to contain sol or lcf overlap with any current liguleless genes nor with previously identified genome-wide association QTL connected to leaf architecture. Analysis of phenotypes across environments further identified a genotype by enviroment interaction determining the strength of the sol x lcf interaction.

In order to leverage novel sequencing techniques for cloning genes in eukaryotic organisms with complex genomes, the false positive rate of variant discovery must be controlled for by experimental design and informatics. We sequenced five lines from three pedigrees of EMS mutagenized Sorghum bicolor, including a pedigree segregating a recessive dwarf mutant. Comparing the sequences of the lines, we were able to identify and eliminate error prone positions. One genomic region contained EMS mutant alleles in dwarfs that were homozygous reference sequence in wild-type siblings and heterozygous in segregating families. This region contained a single non-synonymous change that cosegregated with dwarfism in a validation population and caused a premature stop codon in the sorghum ortholog encoding the giberellic acid biosynthetic enzyme entkaurene oxidase. Application of exogenous giberillic acid rescued the mutant phenotype.Our method for mapping did not require outcrossing and introduced no segregation variance. This enables work when line crossing is complicated by life history, permitting gene discovery outside of genetic models.This inverts the historical approach of first using recombination to define a locus and then sequencing genes. Our formally identical approach first sequences all the genes and then seeks co-segregation with the trait.Mutagenized lines lacking obvious phenotypic alterations are available for an extention of this approach: mapping with a known marker set in a line that is phenotypically identical to starting material for EMS mutant generation.

Centromere positions on 7 maize chromosomes were compared on the basis of data from 4 to 6 mapping techniques per chromosome. Centromere positions were first located relative to molecular markers by means of radiation hybrid lines and centric fission lines recovered from oat-maize chromosome addition lines. These centromere positions were then compared with new data from centric fission lines recovered from maize plants, half-tetrad mapping, and fluorescence in situ hybridizations and to data from earlier studies. Surprisingly, the choice of mapping technique was not the critical determining factor. Instead, on 4 chromosomes, results from all techniques were consistent with a single centromere position. On chromosomes 1, 3, and 6, centromere positions were not consistent even in studies using the same technique. The conflicting centromere map positions on chromosomes 1, 3, and 6 could be explained by pericentric inversions or alternative centromere positions on these chromosomes.

The precise detection of causal DNA mutations (deoxyribonucleic acid) is very crucial for forward genetic studies. Several sources of errors contribute to false‐positive detections by current variant‐calling algorithms, which impact associating phenotypes with genotypes. To improve the accuracy of mutation detection, we implemented a binning method for the accurate detection of likely ethyl methanesulfonate (EMS)‐induced mutations in a sequenced mutant population. We also implemented a clustering algorithm for detecting likely false negatives with high accuracy. Sorghum bicolor is a very valuable crop species with tremendous potential for uncovering novel gene functions associated with highly desirable agronomical traits. We demonstrate the precision of the described approach in the detection of likely EMS‐induced mutations from the publicly available low‐cost sequencing of the M3 generation from 600 sorghum BTx623 mutants. The approach detected 3,274,606 single nucleotide polymorphisms (SNPs), of which 96% (3,141,908) were G/C to A/T DNA substitutions, as expected by EMS‐mutagenesis mode of action. We demonstrated the general applicability of the described method and showed a high concordance, 94% (3,074,759) SNPs overlap between SAMtools‐based and GATK‐based variant‐calling algorithms. Our clustering algorithm uncovered evidence for an additional 223,048 likely false‐negative shared EMS‐induced mutations. The final 3,497,654 SNPs represent an 87% increase in SNPs detected from the previous analysis of the mutant population, with an average of one SNP per 125 kb in the sorghum genome. Annotation of the final SNPs revealed 10,263 high‐impact and 136,639 moderate‐impact SNPs, including 7217 stop‐gained mutations, which averages 12 stop‐gained mutations per mutant, and four high‐ or medium‐impact SNPs per sorghum gene. We have implemented a public search database for this new genetic resource of 30,285 distinct sorghum genes containing medium‐ or high‐impact EMS‐induced mutations. Seedstock for a select 486 of the 600 described mutants are publicly available in the Germplasm Resources Information Network (GRIN) database.

The precise detection of causal DNA mutations is very crucial for forward genetic studies. Several sources of errors contribute to false-positive detections by current variant-calling algorithms, and these impact associating phenotypes with genotypes. To improve the accuracy of mutation detection we propose and implemented a high-resolution binning method for the accurate detection of likely EMS-induced mutations in a sequenced mutant population. The approach also incorporates a novel clustering algorithm for detecting likely false negatives with high accuracy. Sorghum bicolor is a very valuable crop species with tremendous potential for uncovering novel gene functions associated with highly desirable agronomical traits. We demonstrate the precision of the proposed method in the detection of likely EMS-induced mutations in the publicly available low-cost sequencing of the M3 generation from 600 sorghum BTx623 mutants. The method detected 3,274,606 single nucleotide polymorphisms (SNPs) of which 96% (3,141,908) were G/C to A/T DNA substitutions, as expected by EMS-mutagenesis action. We demonstrated the general applicability of the method, and showed a high concordance, 94% (3,074,759) SNPs overlap between SAMtools-based and GATK-based variant-calling algorithms. We also implemented a novel clustering algorithm which uncovered evidence for an additional 223,048 likely false-negative shared EMS-induced mutations. The final 3,497,654 SNPs represents an 87% increase in SNPs detected in the previous analysis of the sorghum mutant population. Annotation of the final SNPs revealed 10,263 high impact and 136,639 moderate impact SNPs, including 7,217 stop-gained mutations, and an average of 12 stop-gained mutations per mutant. We have implemented a public search database for this new genetic resource of 30,285 distinct sorghum genes containing medium or high impact EMS-induced mutations. Seedstock for a select 486 of the 600 described mutants are publicly available in the Germplasm Resources Information Network (GRIN) database.