Forward Genetics by Sequencing EMS Variation-Induced Inbred Lines

Augmenting tomato functional genomics with a genome-wide induced genetic variation resource

Induced mutations accelerate crop improvement by providing novel disease resistance and yield alleles. However, the alleles with no perceptible phenotype but have an altered function remain hidden in mutagenized plants. The whole-genome sequencing (WGS) of mutagenized individuals uncovers the complete spectrum of mutations in the genome. Genome-wide induced mutation resources can improve the targeted breeding of tomatoes and facilitate functional genomics. In this study, we sequenced 132 doubly ethyl methanesulfonate (EMS)-mutagenized lines of tomato and detected approximately 41 million novel mutations and 5.5 million short InDels not present in the parental cultivar. Approximately 97% of the genome had mutations, including the genes, promoters, UTRs, and introns. More than one-third of genes in the mutagenized population had one or more deleterious mutations predicted by Sorting Intolerant From Tolerant (SIFT). Nearly one-fourth of deleterious genes mapped on tomato metabolic pathways modulate multiple pathway steps. In addition to the reported GC>AT transition bias for EMS, our population also had a substantial number of AT>GC transitions. Comparing mutation frequency among synonymous codons revealed that the most preferred codon is the least mutagenic toward EMS. The validation of a potato leaf-like mutation, reduction in carotenoids in ζ-carotene isomerase mutant fruits, and chloroplast relocation loss in phototropin1 mutant validated the mutation discovery pipeline. Our database makes a large repertoire of mutations accessible to functional genomics studies and breeding of tomatoes.

Identification and Analysis of Candidate Genes Associated with Yield Structure Traits and Maize Yield Using Next-Generation Sequencing Technology

The main challenge of agriculture in the 21st century is the continuous increase in food production. In addition to ensuring food security, the goal of modern agriculture is the continued development and production of plant-derived biomaterials. Conventional plant breeding methods do not allow breeders to achieve satisfactory results in obtaining new varieties in a short time. Currently, advanced molecular biology tools play a significant role worldwide, markedly contributing to biological progress. The aim of this study was to identify new markers linked to candidate genes determining grain yield. Next-generation sequencing, gene association, and physical mapping were used to identify markers. An additional goal was to also optimize diagnostic procedures to identify molecular markers on reference materials. As a result of the conducted research, 19 SNP markers significantly associated with yield structure traits in maize were identified. Five of these markers (28629, 28625, 28640, 28649, and 29294) are located within genes that can be considered candidate genes associated with yield traits. For two markers (28639 and 29294), different amplification products were obtained on the electrophorograms. For marker 28629, a specific product of 189 bp was observed for genotypes 1, 4, and 10. For marker 29294, a specific product of 189 bp was observed for genotypes 1 and 10. Both markers can be used for the preliminary selection of well-yielding genotypes.

A cryptic natural variant allele of BYPASS2 suppresses the bypass1 mutant phenotype

The Arabidopsis (Arabidopsis thaliana) BYPASS1 (BPS1) gene encodes a protein with no functionally characterized domains, and loss-of-function mutants (e.g. bps1-2 in Col-0) present a severe growth arrest phenotype that is evoked by a root-derived graft-transmissible small molecule that we call dalekin. The root-to-shoot nature of dalekin signaling suggests it could be an endogenous signaling molecule. Here, we report a natural variant screen that allowed us to identify enhancers and suppressors of the bps1-2 mutant phenotype (in Col-0). We identified a strong semi-dominant suppressor in the Apost-1 accession that largely restored shoot development in bps1 and yet continued to overproduce dalekin. Using bulked segregant analysis and allele-specific transgenic complementation, we showed that the suppressor is the Apost-1 allele of a BPS1 paralog, BYPASS2 (BPS2). BPS2 is one of four members of the BPS gene family in Arabidopsis, and phylogenetic analysis demonstrated that the BPS family is conserved in land plants and the four Arabidopsis paralogs are retained duplicates from whole genome duplications. The strong conservation of BPS1 and paralogous proteins throughout land plants, and the similar functions of paralogs in Arabidopsis, suggests that dalekin signaling might be retained across land plants.

Systematic prediction of EMS-induced mutations in a sorghum mutant population

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

Unraveling novel and rare mutations for alpha-amylase and key transcription factors in EMS-induced wheat mutants for amylose by TILLING

BACKGROUND

TILLING (Targeting Induced Local Lesions in Genomes) is a reverse-genetic strategy that is used to locate an allelic series of induced point mutations in genes of interest. High-throughput TILLING allows the rapid and cost-effective detection of induced point mutations in populations of chemically mutagenized individuals. Grain amylose content is the major constraints for its nutritional quality and have drawn research interest. Identification of allelic variations in genes involved in starch biosynthesis in wheat endosperm is pre-requisite to amenable for nutritional quality improvement.

METHODS AND RESULTS

In this study, 44 EMS-induced (M4 generation) mutant lines having variation for amylose content were used for TILLING sequencing. Overall 2098.08 kb of the sequence was analyzed, and the average mutation density was 1/65.56 kb. In analysis, at the high depth score a total of 32 variations were identified including three natural variations, 76% transitions, 10% transversions, and 14% InDels respectively. The substitutions led to intronic variants, UTRs and up-downstream gene variants in Alpha-amylase, TabZIP77.1, TabZIP1 and Myb respectively. In the Myb transcription factor two missense mutations recorded namely Myb_7B c.680G > A and c.1358 T > C led to p.Gly227Asp and p.Met453Thr and c.1390G > A one substitution in Myb_7D led to p.Val464Ile.

CONCLUSION

The identified missense substitutions were predicted to affect the protein function; hence they may have a probable role in context to the amylose content in mutants. The mutations ascertained in the current study will help in gene discovery in wheat and identified mutants can be used as genetic resources to improve nutritional quality of wheat.

Frequency and Spectrum of Mutations Induced by Gamma Rays Revealed by Phenotype Screening and Whole-Genome Re-Sequencing in Arabidopsis thaliana

Genetic variations are an important source of germplasm diversity, as it provides an allele resource that contributes to the development of new traits for plant breeding. Gamma rays have been widely used as a physical agent for mutation creation in plants, and their mutagenic effect has attracted extensive attention. However, few studies are available on the comprehensive mutation profile at both the large-scale phenotype mutation screening and whole-genome mutation scanning. In this study, biological effects on M₁ generation, large-scale phenotype screening in M₂ generation, as well as whole-genome re-sequencing of seven M₃ phenotype-visible lines were carried out to comprehensively evaluate the mutagenic effects of gamma rays on Arabidopsis thaliana. A total of 417 plants with visible mutated phenotypes were isolated from 20,502 M₂ plants, and the phenotypic mutation frequency of gamma rays was 2.03% in Arabidopsis thaliana. On average, there were 21.57 single-base substitutions (SBSs) and 11.57 small insertions and deletions (InDels) in each line. Single-base InDels accounts for 66.7% of the small InDels. The genomic mutation frequency was 2.78 × 10⁻¹⁰/bp/Gy. The ratio of transition/transversion was 1.60, and 64.28% of the C > T events exhibited the pyrimidine dinucleotide sequence; 69.14% of the small InDels were located in the sequence with 1 to 4 bp terminal microhomology that was used for DNA end rejoining, while SBSs were less dependent on terminal microhomology. Nine genes, on average, were predicted to suffer from functional alteration in each re-sequenced line. This indicated that a suitable mutation gene density was an advantage of gamma rays when trying to improve elite materials for one certain or a few traits. These results will aid the full understanding of the mutagenic effects and mechanisms of gamma rays and provide a basis for suitable mutagen selection and parameter design, which can further facilitate the development of more controlled mutagenesis methods for plant mutation breeding.

Sorghum genetic, genomic, and breeding resources

Sorghum research has entered an exciting and fruitful era due to the genetic, genomic, and breeding resources that are now available to researchers and plant breeders. As the world faces the challenges of a rising population and a changing global climate, new agricultural solutions will need to be developed to address the food and fiber needs of the future. To that end, sorghum will be an invaluable crop species as it is a stress-resistant C₄ plant that is well adapted for semi-arid and arid regions. Sorghum has already remained as a staple food crop in many parts of Africa and Asia and is critically important for animal feed and niche culinary applications in other regions, such as the United States. In addition, sorghum has begun to be developed into a promising feedstock for forage and bioenergy production. Due to this increasing demand for sorghum and its potential to address these needs, the continuous development of powerful community resources is required. These resources include vast collections of sorghum germplasm, high-quality reference genome sequences, sorghum association panels for genome-wide association studies of traits involved in food and bioenergy production, mutant populations for rapid discovery of causative genes for phenotypes relevant to sorghum improvement, gene expression atlas, and online databases that integrate all resources and provide the sorghum community with tools that can be used in breeding and genomic studies. Used in tandem, these valuable resources will ensure that the rate, quality, and collaborative potential of ongoing sorghum improvement efforts is able to rival that of other major crops.

Metabolomics and complementary techniques to investigate the plant phytochemical cosmos

Covering: up to 2021Plants and their associated microbial communities are known to produce millions of metabolites, a majority of which are still not characterized and are speculated to possess novel bioactive properties. In addition to their role in plant physiology, these metabolites are also relevant as existing and next-generation medicine candidates. Elucidation of the plant metabolite diversity is thus valuable for the successful exploitation of natural resources for humankind. Herein, we present a comprehensive review on recent metabolomics approaches to illuminate molecular networks in plants, including chemical isolation and enzymatic production as well as the modern metabolomics approaches such as stable isotope labeling, ultrahigh-resolution mass spectrometry, metabolome imaging (spatial metabolomics), single-cell analysis, cheminformatics, and computational mass spectrometry. Mass spectrometry-based strategies to characterize plant metabolomes through metabolite identification and annotation are described in detail. We also highlight the use of phytochemical genomics to mine genes associated with specialized metabolites' biosynthesis. Understanding the metabolic diversity through biotechnological advances is fundamental to elucidate the functions of the plant-derived specialized metabolome.

slim shady is a novel allele of PHYTOCHROME B present in the T-DNA line SALK_015201

Auxin is a hormone that is required for hypocotyl elongation during seedling development. In response to auxin, rapid changes in transcript and protein abundance occur in hypocotyls, and some auxin responsive gene expression is linked to hypocotyl growth. To functionally validate proteomic studies, a reverse genetics screen was performed on mutants in auxin-regulated proteins to identify novel regulators of plant growth. This uncovered a long hypocotyl mutant, which we called slim shady, in an annotated insertion line in IMMUNOREGULATORY RNA-BINDING PROTEIN (IRR). Overexpression of the IRR gene failed to rescue the slim shady phenotype and characterization of a second T-DNA allele of IRR found that it had a wild-type (WT) hypocotyl length. The slim shady mutant has an elevated expression of numerous genes associated with the brassinosteroid-auxin-phytochrome (BAP) regulatory module compared to WT, including transcription factors that regulate brassinosteroid, auxin, and phytochrome pathways. Additionally, slim shady seedlings fail to exhibit a strong transcriptional response to auxin. Using whole genome sequence data and genetic complementation analysis with SALK_015201C, we determined that a novel single nucleotide polymorphism in PHYTOCHROME B was responsible for the slim shady phenotype. This is predicted to induce a frameshift and premature stop codon at leucine 1125, within the histidine kinase-related domain of the carboxy terminus of PHYB, which is required for phytochrome signaling and function. Genetic complementation analyses with phyb-9 confirmed that slim shady is a mutant allele of PHYB. This study advances our understanding of the molecular mechanisms in seedling development, by furthering our understanding of how light signaling is linked to auxin-dependent cell elongation. Furthermore, this study highlights the importance of confirming the genetic identity of research material before attributing phenotypes to known mutations sourced from T-DNA stocks.

BSAseq: an interactive and integrated web-based workflow for identification of causal mutations in bulked F2 populations

SUMMARY

With the advance of next-generation sequencing technologies and reductions in the costs of these techniques, bulked segregant analysis (BSA) has become not only a powerful tool for mapping quantitative trait loci but also a useful way to identify causal gene mutations underlying phenotypes of interest. However, due to the presence of background mutations and errors in sequencing, genotyping, and reference assembly, it is often difficult to distinguish true causal mutations from background mutations. In this study, we developed the BSAseq workflow, which includes an automated bioinformatics analysis pipeline with a probabilistic model for estimating the linked region (the region linked to the causal mutation) and an interactive Shiny web application for visualizing the results. We deeply sequenced a sorghum male-sterile parental line (ms8) to capture the majority of background mutations in our bulked F2 data. We applied the workflow to 11 bulked sorghum F2 populations and 1 rice F2 population and identified the true causal mutation in each population. The workflow is intuitive and straightforward, facilitating its adoption by users without bioinformatics analysis skills. We anticipate that the BSAseq workflow will be broadly applicable to the identification of causal mutations for many phenotypes of interest.

AVAILABILITY AND IMPLEMENTATION

BSAseq is freely available on https://www.sciapps.org/page/bsa.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mapping-by-Sequencing via MutMap Identifies a Mutation in ZmCLE7 Underlying Fasciation in a Newly Developed EMS Mutant Population in an Elite Tropical Maize Inbred

Induced point mutations are important genetic resources for their ability to create hypo- and hypermorphic alleles that are useful for understanding gene functions and breeding. However, such mutant populations have only been developed for a few temperate maize varieties, mainly B73 and W22, yet no tropical maize inbred lines have been mutagenized and made available to the public to date. We developed a novel Ethyl Methanesulfonate (EMS) induced mutation resource in maize comprising 2050 independent M2 mutant families in the elite tropical maize inbred ML10. By phenotypic screening, we showed that this population is of comparable quality with other mutagenized populations in maize. To illustrate the usefulness of this population for gene discovery, we performed rapid mapping-by-sequencing to clone a fasciated-ear mutant and identify a causal promoter deletion in ZmCLE7 (CLE7). Our mapping procedure does not require crossing to an unrelated parent, thus is suitable for mapping subtle traits and ones affected by heterosis. This first EMS population in tropical maize is expected to be very useful for the maize research community. Also, the EMS mutagenesis and rapid mapping-by-sequencing pipeline described here illustrate the power of performing forward genetics in diverse maize germplasms of choice, which can lead to novel gene discovery due to divergent genetic backgrounds.

Development and characterization of an EMS-mutagenized wheat population and identification of salt-tolerant wheat lines

BACKGROUND

Triticum aestivum (wheat) is one of the world's oldest crops and has been used for >8000 years as a food crop in North Africa, West Asia and Europe. Today, wheat is one of the most important sources of grain for humans, and is cultivated on greater areas of land than any other crop. As the human population increases and soil salinity becomes more prevalent, there is increased pressure on wheat breeders to develop salt-tolerant varieties in order to meet growing demands for yield and grain quality. Here we developed a mutant wheat population using the moderately salt-tolerant Bangladeshi variety BARI Gom-25, with the primary goal of further increasing salt tolerance.

RESULTS

After titrating the optimal ethyl methanesulfonate (EMS) concentration, ca 30,000 seeds were treated with 1% EMS, and 1676 lines, all originating from single seeds, survived through the first four generations. Most mutagenized lines showed a similar phenotype to BARI Gom-25, although visual differences such as dwarfing, giant plants, early and late flowering and altered leaf morphology were seen in some lines. By developing an assay for salt tolerance, and by screening the mutagenized population, we identified 70 lines exhibiting increased salt tolerance. The selected lines typically showed a 70% germination rate on filter paper soaked in 200 mM NaCl, compared to 0-30% for BARI Gom-25. From two of the salt-tolerant OlsAro lines (OA42 and OA70), genomic DNA was sequenced to 15x times coverage. A comparative analysis against the BARI Gom-25 genomic sequence identified a total of 683,201 (OA42), and 768,954 (OA70) SNPs distributed throughout the three sub-genomes (A, B and D). The mutation frequency was determined to be approximately one per 20,000 bp. All the 70 selected salt-tolerant lines were tested for root growth in the laboratory, and under saline field conditions in Bangladesh. The results showed that all the lines selected for tolerance showed a better salt tolerance phenotype than both BARI Gom-25 and other local wheat varieties tested.

CONCLUSION

The mutant wheat population developed here will be a valuable resource in the development of novel salt-tolerant varieties for the benefit of saline farming.

Mutagenesis of Puccinia graminis f. sp. tritici and Selection of Gain-of-Virulence Mutants

Wheat stem rust caused by the fungus Puccinia graminis f. sp. tritici (Pgt), is regaining prominence due to the recent emergence of virulent isolates and epidemics in Africa, Europe and Central Asia. The development and deployment of wheat cultivars with multiple stem rust resistance (Sr) genes stacked together will provide durable resistance. However, certain disease resistance genes can suppress each other or fail in particular genetic backgrounds. Therefore, the function of each Sr gene must be confirmed after incorporation into an Sr-gene stack. This is difficult when using pathogen disease assays due to epistasis from recognition of multiple avirulence (Avr) effectors. Heterologous delivery of single Avr effectors can circumvent this limitation, but this strategy is currently limited by the paucity of cloned Pgt Avrs. To accelerate Avr gene cloning, we outline a procedure to develop a mutant population of Pgt spores and select for gain-of-virulence mutants. We used ethyl methanesulphonate (EMS) to mutagenize urediniospores and create a library of > 10,000 independent mutant isolates that were combined into 16 bulks of ~658 pustules each. We sequenced random mutants and determined the average mutation density to be 1 single nucleotide variant (SNV) per 258 kb. From this, we calculated that a minimum of three independently derived gain-of-virulence mutants is required to identify a given Avr gene. We inoculated the mutant library onto plants containing Sr43, Sr44, or Sr45 and obtained 9, 4, and 14 mutants with virulence toward Sr43, Sr44, or Sr45, respectively. However, only mutants identified on Sr43 and Sr45 maintained their virulence when reinolculated onto the lines from which they were identified. We further characterized 8 mutants with virulence toward Sr43. These also maintained their virulence profile on the stem rust international differential set containing 20 Sr genes, indicating that they were most likely not accidental contaminants. In conclusion, our method allows selecting for virulent mutants toward targeted resistance (R) genes. The development of a mutant library from as little as 320 mg spores creates a resource that enables screening against several R genes without the need for multiple rounds of spore multiplication and mutagenesis.

Common metabolic networks contribute to carbon sink strength of sorghum internodes: implications for bioenergy improvement

BACKGROUND

Sorghum bicolor (L.) is an important bioenergy source. The stems of sweet sorghum function as carbon sinks and accumulate large amounts of sugars and lignocellulosic biomass and considerable amounts of starch, therefore providing a model of carbon allocation and accumulation for other bioenergy crops. While omics data sets for sugar accumulation have been reported in different genotypes, the common features of primary metabolism in sweet genotypes remain unclear. To obtain a cohesive and comparative picture of carbohydrate metabolism between sorghum genotypes, we compared the phenotypes and transcriptome dynamics of sugar-accumulating internodes among three different sweet genotypes (Della, Rio, and SIL-05) and two non-sweet genotypes (BTx406 and R9188).

RESULTS

Field experiments showed that Della and Rio had similar dynamics and internode patterns of sugar concentration, albeit distinct other phenotypes. Interestingly, cellulose synthases for primary cell wall and key genes in starch synthesis and degradation were coordinately upregulated in sweet genotypes. Sweet sorghums maintained active monolignol biosynthesis compared to the non-sweet genotypes. Comparative RNA-seq results support the role of candidate Tonoplast Sugar Transporter gene (TST), but not the Sugars Will Eventually be Exported Transporter genes (SWEETs) in the different sugar accumulations between sweet and non-sweet genotypes.

CONCLUSIONS

Comparisons of the expression dynamics of carbon metabolic genes across the RNA-seq data sets identify several candidate genes with contrasting expression patterns between sweet and non-sweet sorghum lines, including genes required for cellulose and monolignol synthesis (CesA, PTAL, and CCR), starch metabolism (AGPase, SS, SBE, and G6P-translocator SbGPT2), and sucrose metabolism and transport (TPP and TST2). The common transcriptome features of primary metabolism identified here suggest the metabolic networks contributing to carbon sink strength in sorghum internodes, prioritize the candidate genes for manipulating carbon allocation with bioenergy purposes, and provide a comparative and cohesive picture of the complexity of carbon sink strength in sorghum stem.

Combining next-generation sequencing and progeny testing for rapid identification of induced recessive and dominant mutations in maize M2 individuals

Molecular identification of mutant alleles responsible for certain phenotypic alterations is a central goal of genetic analyses. In this study we describe a rapid procedure suitable for the identification of induced recessive and dominant mutations applied to two Zea mays mutants expressing a dwarf and a pale green phenotype, respectively, which were obtained through pollen ethyl methanesulfonate (EMS) mutagenesis. First, without prior backcrossing, induced mutations (single nucleotide polymorphisms, SNPs) segregating in a (M₂ ) family derived from a heterozygous (M₁ ) parent were identified using whole-genome shotgun (WGS) sequencing of a small number of (M₂ ) individuals with mutant and wild-type phenotypes. Second, the state of zygosity of the mutation causing the phenotype was determined for each sequenced individual by phenotypic segregation analysis of the self-pollinated (M₃ ) offspring. Finally, we filtered for segregating EMS-induced SNPs whose state of zygosity matched the determined state of zygosity of the mutant locus in each sequenced (M₂ ) individuals. Through this procedure, combining sequencing of individuals and Mendelian inheritance, three and four SNPs in linkage passed our zygosity filter for the homozygous dwarf and heterozygous pale green mutation, respectively. The dwarf mutation was found to be allelic to the an1 locus and caused by an insertion in the largest exon of the AN1 gene. The pale green mutation affected the nuclear W2 gene and was caused by a non-synonymous amino acid exchange in encoded chloroplast DNA polymerase with a predicted deleterious effect. This coincided with lower cpDNA levels in pale green plants.

Bulked-Segregant Analysis Coupled to Whole Genome Sequencing (BSA-Seq) for Rapid Gene Cloning in Maize

Forward genetics remains a powerful method for revealing the genes underpinning organismal form and function, and for revealing how these genes are tied together in gene networks. In maize, forward genetics has been tremendously successful, but the size and complexity of the maize genome made identifying mutant genes an often arduous process with traditional methods. The next generation sequencing revolution has allowed for the gene cloning process to be significantly accelerated in many organisms, even when genomes are large and complex. Here, we describe a bulked-segregant analysis sequencing (BSA-Seq) protocol for cloning mutant genes in maize. Our simple strategy can be used to quickly identify a mapping interval and candidate single nucleotide polymorphisms (SNPs) from whole genome sequencing of pooled F2 individuals. We employed this strategy to identify narrow odd dwarf as an enhancer of teosinte branched1, and to identify a new allele of defective kernel1 Our method provides a quick, simple way to clone genes in maize.

Whole-Genome Sequence Accuracy Is Improved by Replication in a Population of Mutagenized Sorghum

The accurate detection of induced mutations is critical for both forward and reverse genetics studies. Experimental chemical mutagenesis induces relatively few single base changes per individual. In a complex eukaryotic genome, false positive detection of mutations can occur at or above this mutagenesis rate. We demonstrate here, using a population of ethyl methanesulfonate (EMS)-treated Sorghum bicolor BTx623 individuals, that using replication to detect false positive-induced variants in next-generation sequencing (NGS) data permits higher throughput variant detection with greater accuracy. We used a lower sequence coverage depth (average of 7×) from 586 independently mutagenized individuals and detected 5,399,493 homozygous single nucleotide polymorphisms (SNPs). Of these, 76% originated from only 57,872 genomic positions prone to false positive variant calling. These positions are characterized by high copy number paralogs where the error-prone SNP positions are at copies containing a variant at the SNP position. The ability of short stretches of homology to generate these error-prone positions suggests that incompletely assembled or poorly mapped repeated sequences are one driver of these error-prone positions. Removal of these false positives left 1,275,872 homozygous and 477,531 heterozygous EMS-induced SNPs, which, congruent with the mutagenic mechanism of EMS, were >98% G:C to A:T transitions. Through this analysis, we generated a collection of sequence indexed mutants of sorghum. This collection contains 4035 high-impact homozygous mutations in 3637 genes and 56,514 homozygous missense mutations in 23,227 genes. Each line contains, on average, 2177 annotated homozygous SNPs per genome, including seven likely gene knockouts and 96 missense mutations. The number of mutations in a transcript was linearly correlated with the transcript length and also the G+C count, but not with the GC/AT ratio. Analysis of the detected mutagenized positions identified CG-rich patches, and flanking sequences strongly influenced EMS-induced mutation rates. This method for detecting false positive-induced mutations is generally applicable to any organism, is independent of the choice of in silico variant-calling algorithm, and is most valuable when the true mutation rate is likely to be low, such as in laboratory-induced mutations or somatic mutation detection in medicine.