Development and utilization of a new chemically-induced soybean library with a high mutation density

Generating homozygous mutant populations of barley microspores by ethyl methanesulfonate treatment

Induced mutations are important for genetic research and breeding. Mutations induced by physical or chemical mutagenesis are usually heterozygous during the early generations. However, mutations must be fixed prior to phenotyping or field trials, which requires additional rounds of self-pollination. Microspore culture is an effective method to produce double-haploid (DH) plants that are fixed homozygotes. In this study, we conducted ethyl methanesulfonate (EMS)-induced mutagenesis of microspore cultures of barley (Hordeum vulgare) cultivar 'Hua30' and landrace 'HTX'. The EMS concentrations were negatively correlated with the efficiency of callus induction and the frequency of mutant plant regeneration. The two genotypes showed different regeneration efficiencies. The phenotypic variation of the regenerated M1 plants and the presence of genome-wide nucleotide mutations, revealed by whole-genome sequencing, highlight the utility of EMS-induced mutagenesis of isolated microspore cultures for developing DH mutants. Genome-wide analysis of the mutation frequency in the regenerated plants revealed that a considerable proportion of mutations resulted from microspore culture (somaclonal variation) rather than EMS-induced mutagenesis. In addition to producing a population of 1972 homozygous mutant lines that are available for future field trials, this study lays the foundation for optimizing the regeneration efficiency of DH plants and the richness of mutations (mainly by fine-tuning the mutagen dosage).

Variations in Total Protein and Amino Acids in the Sequenced Sorghum Mutant Library

Sorghum (Sorghum bicolor) is the fifth most important cereal crop worldwide; however, its utilization in food products can be limited due to reduced nutritional quality related to amino acid composition and protein digestibility in cooked products. Low essential amino acid levels and digestibility are influenced by the composition of the sorghum seed storage proteins, kafirins. In this study, we report a core collection of 206 sorghum mutant lines with altered seed storage proteins. Wet lab chemistry analysis was conducted to evaluate the total protein content and 23 amino acids, including 19 protein-bound and 4 non-protein amino acids. We identified mutant lines with diverse compositions of essential and non-essential amino acids. The highest total protein content in these lines was almost double that of the wild-type (BTx623). The mutants identified in this study can be used as a genetic resource to improve the sorghum grain quality and determine the molecular mechanisms underlying the biosynthesis of storage protein and starch in sorghum seeds.

Construction of a density mutant collection in bitter gourd via new germplasms innovation and gene functional study

Although a few studies have elucidated the creation of bitter gourd mutants, the suitable concentration and duration of ethyl methanesulfonate (EMS) mutagenesis have not been determined. In this study, mutant collection was conducted to create new germplasms and widen genetic diversity. By employing the seeds of the inbred line Y52 as the mutagenic material, EMS as the mutagen, and the suitable mutagenic conditions for bitter gourd seeds (EMS concentration 0.2%, mutagenic time 10 h), we mutated 10,000 seeds and acquired 3223 independent M₁ lines. For the randomly selected 1000 M₂ lines, 199 M₂ lines with visible phenotypes were found, and 167 M₂ lines were mutants of fruit shape, size, and tubercles. Furthermore, fourteen dwarf, eleven leaf color, five leaf shape, and eight meristem defect mutants were discovered in this mutant collection. In addition, three lines of 1253, 2284, and 3269 represented recessive mutants crossed with Y52. Furthermore, the yellow leaf lines of 2284 and 3269 were not mutated at the same gene locus. This study constructed a mutant collection through innovative new germplasms and provided valuable resources for bitter gourd breeding and functional gene research.

Genetic dissection of the soybean dwarf mutant dm with integrated genomic, transcriptomic and methylomic analyses

Plant height affects crop production and breeding practices, while genetic control of dwarfism draws a broad interest of researchers. Dwarfism in soybean (Glycine max) is mainly unexplored. Here, we characterized a dwarf mutant dm screened from ethyl methanesulfonate (EMS) mutated seeds of the soybean cultivar Zhongpin 661(ZP). Phenotypically, dm showed shorter and thinner stems, smaller leaves, and more nodes than ZP under greenhouse conditions. Genetically, whole-genome sequencing and comparison revealed that 210K variants of SNPs and InDel in ZP relative to the soybean reference genome Williams82, and EMS mutagenesis affected 636 genes with variants predicted to have a large impact on protein function in dm. Whole-genome methylation sequencing found 704 differentially methylated regions in dm. Further whole-genome RNA-Seq based transcriptomic comparison between ZP and dm leaves revealed 687 differentially expressed genes (DEGs), including 263 up-regulated and 424 down-regulated genes. Integrated omics analyses revealed 11 genes with both differential expressions and DNA variants, one gene with differential expression and differential methylation, and three genes with differential methylation and sequence variation, worthy of future investigation. Genes in cellulose, fatty acids, and energy-associated processes could be the key candidate genes for the dwarf phenotype. This study provides genetic clues for further understanding of the genetic control of dwarfism in soybean. The genetic resources could help to inbreed new cultivars with a desirable dwarf characteristic.

MINI BODY1, encoding a MATE/DTX family transporter, affects plant architecture in mungbean (Vigna radiata L.)

It has been shown that multidrug and toxic compound extrusion/detoxification (MATE/DTX) family transporters are involved in the regulation of plant development and stress response. Here, we characterized the mini body1 (mib1) mutants in mungbean, which gave rise to increased branches, pentafoliate compound leaves, and shortened pods. Map-based cloning revealed that MIB1 encoded a MATE/DTX family protein in mungbean. qRT-PCR analysis showed that MIB1 was expressed in all tissues of mungbean, with the highest expression level in the young inflorescence. Complementation assays in Escherichia coli revealed that MIB1 potentially acted as a MATE/DTX transporter in mungbean. It was found that overexpression of the MIB1 gene partially rescued the shortened pod phenotype of the Arabidopsis dtx54 mutant. Transcriptomic analysis of the shoot buds and young pods revealed that the expression levels of several genes involved in the phytohormone pathway and developmental regulators were altered in the mib1 mutants. Our results suggested that MIB1 plays a key role in the control of plant architecture establishment in mungbean.

Soybean Protein and Oil Variants Identified through a Forward Genetic Screen for Seed Composition

Mutagenesis remains an important tool in soybean biology. In classical plant mutation breeding, mutagenesis has been a trusted approach for decades, creating stable non-transgenic variation, and many mutations have been incorporated into germplasm for several crops, especially to introduce favorable seed composition traits. We performed a genetic screen for aberrant oil or protein composition of soybean seeds, and as a result isolated over 100 mutant lines for seed composition phenotypes, with particular interest in high protein or high oil phenotypes. These lines were followed for multiple seasons and generations to select the most stable traits for further characterization. Through backcrossing and outcrossing experiments, we determined that a subset of the lines showed recessive inheritance, while others showed a dominant inheritance pattern that suggests the involvement of multiple loci and genetic mechanisms. These lines can be used as a resource for future studies of the genetic control of seed protein and oil content in soybean.

An enhancing effect attributed to a nonsynonymous mutation in SOYBEAN SEED SIZE 1, a SPINDLY-like gene, is exploited in soybean domestication and improvement

Soybean (Glycine max) was domesticated from its wild relative Glycine soja. One-hundred-seed weight is one of the most important domesticated traits determining soybean yield; however, its underlying genetic basis remains elusive. We characterized a soybean seed size 1 (sss1) mutant featuring large seeds compared to its wild-type background. Positional cloning revealed that the candidate gene GmSSS1 encoded a SPINDLY homolog and was co-located in a well-identified quantitative trait locus (QTL)-rich region on chromosome 19. Knocking out GmSSS1 resulted in small seeds, while overexpressing GmSSS1/Gmsss1 induced large seeds. Modulating GmSSS1/Gmsss1 in transgenic plants can positively influence cell expansion and cell division. Relative to GmSSS1, one mutation leading to an E to Q substitution at the 182^nd residue in Gmsss1 conferred an enhancing effect on seed weight. GmSSS1 underwent diversification in wild-type and cultivated soybean, and the alleles encoding the Gmsss1-type substitution of 182^nd -Q, which originated along the central and downstream parts of the Yellow River, were selected and expanded during soybean domestication and improvement. We cloned the causative gene for the sss1 mutant, which is linked with a seed weight QTL, identified an elite allele of this gene for increasing seed weight, and provided new insights into soybean domestication and breeding.

GmIAA27 Encodes an AUX/IAA Protein Involved in Dwarfing and Multi-Branching in Soybean

Soybean plant height and branching affect plant architecture and yield potential in soybean. In this study, the mutant dmbn was obtained by treating the cultivar Zhongpin 661 with ethylmethane sulfonate. The dmbn mutant plants were shorter and more branched than the wild type. The genetic analysis showed that the mutant trait was controlled by a semi-dominant gene. The candidate gene was fine-mapped to a 91 kb interval on Chromosome 9 by combining BSA-seq and linkage analysis. Sequence analysis revealed that Glyma.09g193000 encoding an Aux/IAA protein (GmIAA27) was mutated from C to T in the second exon of the coding region, resulting to amino acid substitution of proline to leucine. Overexpression of the mutant type of this gene in Arabidopsis thaliana inhibited apical dominance and promoted lateral branch development. Expression analysis of GmIAA27 and auxin response genes revealed that some GH3 genes were induced. GmIAA27 relies on auxin to interact with TIR1, whereas Gmiaa27 cannot interact with TIR1 owing to the mutation in the degron motif. Identification of this unique gene that controls soybean plant height and branch development provides a basis for investigating the mechanisms regulating soybean plant architecture development.

Protein profiling of fast neutron soybean mutant seeds reveals differential accumulation of seed and iron storage proteins

A fast neutron (FN) radiated mutant soybean (Glycine max (L.) Merr., Fabaceae) displaying large duplications exhibited an increase in total seed protein content. A tandem mass tag (TMT) based protein profiling of matured seeds resulted in the identification of 4338 proteins. Gene duplication resulted in a significant increase in several seed storage proteins and protease inhibitors. Among the storage proteins, basic 7 S globulin, glycinin G4, and beta-conglycinin showed higher abundance in matured FN mutant seeds in addition to protease inhibitors. A significantly higher abundance of L-ascorbate peroxidases, acid phosphatases, and iron storage proteins was also observed. A higher amount of albumin, sucrose synthase, iron storage, and ascorbate family proteins in the mutant seeds was observed at the mid-stage of seed filling. We anticipate that the duplicated genes might have a cascading effect on the genome constituents, thus, resulting in increased storage and iron-containing protein content in the mutant seeds.

Comparison between Germinated Seed and Isolated Microspore EMS Mutagenesis in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Mutagenesis is an important tool for breeding and genomic research. In this study, the germinated seeds and isolated microspores of a double haploid line ‘FT’ were treated with EMS, respectively, with the aim of comparing the effects of the two approaches on generating mutants in Chinese cabbage. For microspore EMS mutagenesis, the isolated microspores were treated with 0.12% EMS for 20 min, a total of 1268 plantlets were obtained, and 15 M1 mutants were screened with a mutation frequency of 1.2%. For seed EMS mutagenesis, 7800 germinated seeds were treated with 0.8% EMS for 12 h, and a total of 701 M2 mutants were screened, with a mutation frequency of 18.78%. In total, 716 mutants with heritable morphological variation including leaf color, leaf shape, leafy head, bolting, and fertility, were obtained from the EMS mutagenesis experiments. Homozygous mutant plants could be screened from M1 lines by microspore mutagenesis, and M2 lines by seed mutagenesis. The mutation frequency was higher in seed mutagenesis than in microspore mutagenesis. Based on these results, we propose that seed EMS mutagenesis is more suitable to generate a large-scale mutant library, and the microspore EMS mutagenesis is conducive to rapidly obtaining homozygous mutants.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Proteomic Profiling of Fast Neutron-Induced Soybean Mutant Unveiled Pathways Associated with Increased Seed Protein Content

Mutagenesis through fast neutron (FN) radiation of soybean resulted in a mutant with a 15% increase in seed protein content. A comparative genomic hybridization analysis confirmed that the mutant is lacking 24 genes located at chromosomes 5 and 10. A tandem mass tag-based proteomic profiling of the wild type and the FN mutant revealed 3,502 proteins, of which 206 proteins exhibited increased abundance and 214 proteins showed decreased abundance. Among the abundant proteins, basic 7S globulin increased fourfold, followed by vacuolar-sorting receptor and protein transporters. The differentially expressed proteins were mapped on the global metabolic pathways. It was observed that there was an enrichment of 29 ribosomal proteins, 16 endoplasmic reticular proteins, and several proteins in export metabolic pathways. The deletion of the sequence-specific DNA binding transcription factor along with 23 other genes may have altered the negative regulation of protein syntheses processes, resulting in an increase in the overall protein content of the mutant seed. This mutant is a valuable resource for researchers to understand the metabolic pathways that may affect an increase in seed protein content (the mass spectrometry data files were submitted to massive.ucsd.edu # MassIVE MSV000084228).

Identification of a cytochrome P450 hydroxylase, CYP81E22, as a causative gene for the high sensitivity of soybean to herbicide bentazon

KEY MESSAGE

A frame shift invoked by a single-base deletion in the gene encoding a cytochrome P450 hydroxylase, CYP81E22, causes the loss of bentazon detoxification function in soybean. Bentazon is an effective herbicide in soybean cultivation applied at post-emergence stages for control of several broadleaf weeds. However, some soybean cultivars are highly sensitive to bentazon and are killed upon application. In this study, the gene related to the high sensitivity of soybean cultivars to bentazon was mapped to chromosome 16, and its location was narrowed down to a 257-kb region where three cytochrome P450 genes were located. In these genes, a single-base deletion of cytosine was detected in the coding region of Glyma.16G149300, CYP81E22, at + 1465 bp downstream from the translation start codon, leading to a frame shift in the open reading frame and creating a premature stop codon. This stop codon resulted in the loss of more than half of the P450, and consequently, the remaining molecule failed to form a functioning protein. This single-base deletion was common among the highly sensitive cultivars screened from the soybean mini-core collection and other previously reported highly sensitive cultivars. Furthermore, we screened plant lines from the targeting-induced local lesions in genomes library of the soybean cultivar Enrei based on a modelled 3D structure of CYP81E22. The lines with mutations in Glyma.16G149300 were highly sensitive to bentazon, which provides strong evidence that Glyma.16G149300 is the gene responsible for high sensitivity to bentazon.

Flavonoids and Isoflavonoids Biosynthesis in the Model Legume Lotus japonicus; Connections to Nitrogen Metabolism and Photorespiration

Phenylpropanoid metabolism represents an important metabolic pathway from which originates a wide number of secondary metabolites derived from phenylalanine or tyrosine, such as flavonoids and isoflavonoids, crucial molecules in plants implicated in a large number of biological processes. Therefore, various types of interconnection exist between different aspects of nitrogen metabolism and the biosynthesis of these compounds. For legumes, flavonoids and isoflavonoids are postulated to play pivotal roles in adaptation to their biological environments, both as defensive compounds (phytoalexins) and as chemical signals in symbiotic nitrogen fixation with rhizobia. In this paper, we summarize the recent progress made in the characterization of flavonoid and isoflavonoid biosynthetic pathways in the model legume Lotus japonicus (Regel) Larsen under different abiotic stress situations, such as drought, the impairment of photorespiration and UV-B irradiation. Emphasis is placed on results obtained using photorespiratory mutants deficient in glutamine synthetase. The results provide different types of evidence showing that an enhancement of isoflavonoid compared to standard flavonol metabolism frequently occurs in Lotus under abiotic stress conditions. The advance produced in the analysis of isoflavonoid regulatory proteins by the use of co-expression networks, particularly MYB transcription factors, is also described. The results obtained in Lotus japonicus plants can be also extrapolated to other cultivated legume species, such as soybean, of extraordinary agronomic importance with a high impact in feeding, oil production and human health.

Single-base deletion in GmCHR5 increases the genistein-to-daidzein ratio in soybean seed

Novel mutant alleles related to isoflavone content are useful for breeding programs to improve the disease resistance and nutritional content of soybean. However, identification of mutant alleles from high-density mutant libraries is expensive and time-consuming because soybean has a large, complicated genome. Here, we identified the gene responsible for increased genistein-to-daidzein ratio in seed of the mutant line F333ES017D9. For this purpose, we used a time- and cost-effective approach based on selective genotyping of a small number of F₂ plants showing the mutant phenotype with nearest-neighboring-nucleotide substitution-high-resolution melting analysis markers, followed by alignment of short reads obtained by next-generation sequencing analysis with the identified locus. In the mutant line, GmCHR5 harbored a single-base deletion that caused a change in the substrate flow in the isoflavone biosynthetic pathway towards genistein. Mutated GmCHR5 was expressed at a lower level during seed development than wild-type GmCHR5. Ectopic overexpression of GmCHR5 increased the production of daidzein derivatives in both the wild-type and mutant plants. The present strategy will be useful for accelerating identification of mutant alleles responsible for traits of interest in agronomically important crops.

Toward a "Green Revolution" for Soybean

Soybean (Glycine max), as an economically important food and oilseedcrop, is a major source of plant proteins and oils. Although considerable progress has been made in increasing the yields of rice, wheat, and maize through the "Green Revolution", little improvements have been made for soybean. With the increasing demand of soybean production and the rapid development of crop breeding technologies, time has come for this important crop to undergo a Green Revolution. Here, we briefly summarize the history of crop breeding and Green Revolution in other crops. We then discuss the possible directions and potential approaches toward achieving a Green Revolution for soybean. We provide our views and perspectives on how to breed new soybean varieties with improved yield.

Generation of a multiplex mutagenesis population via pooled CRISPR-Cas9 in soya bean

The output of genetic mutant screenings in soya bean [Glycine max (L.) Merr.] has been limited by its paleopolypoid genome. CRISPR-Cas9 can generate multiplex mutants in crops with complex genomes. Nevertheless, the transformation efficiency of soya bean remains low and, hence, remains the major obstacle in the application of CRISPR-Cas9 as a mutant screening tool. Here, we report a pooled CRISPR-Cas9 platform to generate soya bean multiplex mutagenesis populations. We optimized the key steps in the screening protocol, including vector construction, sgRNA assessment, pooled transformation, sgRNA identification and gene editing verification. We constructed 70 CRISPR-Cas9 vectors to target 102 candidate genes and their paralogs which were subjected to pooled transformation in 16 batches. A population consisting of 407 T0 lines was obtained containing all sgRNAs at an average mutagenesis frequency of 59.2%, including 35.6% lines carrying multiplex mutations. The mutation frequency in the T1 progeny could be increased further despite obtaining a transgenic chimera. In this population, we characterized gmric1/gmric2 double mutants with increased nodule numbers and gmrdn1-1/1-2/1-3 triple mutant lines with decreased nodulation. Our study provides an advanced strategy for the generation of a targeted multiplex mutant population to overcome the gene redundancy problem in soya bean as well as in other major crops.

An efficient TILLING platform for cultivated tobacco

Targeting-induced local lesions in genomes (TILLING) is a powerful reverse-genetics tool that enables high-throughput screening of genomic variations in plants. Although TILLING has been developed for many diploid plants, the technology has been used in very few polyploid species due to their genomic complexity. Here, we established an efficient capillary electrophoresis-based TILLING platform for allotetraploid cultivated tobacco (Nicotiana tabacum L.) using an ethyl methanesulfonate (EMS)-mutagenized population of 1,536 individuals. We optimized the procedures for endonuclease preparation, leaf tissue sampling, DNA extraction, normalization, pooling, PCR amplification, heteroduplex formation, and capillary electrophoresis. In a test screen using seven target genes with eight PCR fragments, we obtained 118 mutants. The mutation density was estimated to be approximately one mutation per 106 kb on average. Phenotypic analyses showed that mutations in two heavy metal transporter genes, HMA2S and HMA4T, led to reduced accumulation of cadmium and zinc, which was confirmed independently using CRISPR/Cas9 to generate knockout mutants. Our results demonstrate that this powerful TILLING platform (available at http://www.croptilling.org) can be used in tobacco to facilitate functional genomics applications.

Assessment of Phenotypic Variations and Correlation among Seed Composition Traits in Mutagenized Soybean Populations

Soybean [Glycine max (L.) Merr.] seed is a valuable source of protein and oil worldwide. Traditionally, the natural variations were heavily used in conventional soybean breeding programs to select desired traits. However, traditional plant breeding is encumbered with low frequencies of spontaneous mutations. In mutation breeding, genetic variations from induced mutations provide abundant sources of alterations in important soybean traits; this facilitated the development of soybean germplasm with modified seed composition traits to meet the different needs of end users. In this study, a total of 2366 'Forrest'-derived M2 families were developed for both forward and reverse genetic studies. A subset of 881 M3 families was forward genetically screened to measure the contents of protein, oil, carbohydrates, and fatty acids. A total of 14 mutants were identified to have stable seed composition phenotypes observed in both M3 and M4 generations. Correlation analyses have been conducted among ten seed composition traits and compared to a collection of 103 soybean germplasms. Mainly, ethyl methanesulfonate (EMS) mutagenesis had a strong impact on the seed-composition correlation that was observed among the 103 soybean germplasms, which offers multiple benefits for the soybean farmers and industry to breed for desired multiple seed phenotypes.

Legume genomics and transcriptomics: From classic breeding to modern technologies

Legumes are essential and play a significant role in maintaining food standards and augmenting physiochemical soil properties through the biological nitrogen fixation process. Biotic and abiotic factors are the main factors limiting legume production. Classical breeding methodologies have been explored extensively about the problem of truncated yield in legumes but have not succeeded at the desired rate. Conventional breeding improved legume genotypes but with more resources and time. Recently, the invention of next-generation sequencing (NGS) and high-throughput methods for genotyping have opened new avenues for research and developments in legume studies. During the last decade, genome sequencing for many legume crops documented. Sequencing and re-sequencing of important legume species have made structural variation and functional genomics conceivable. NGS and other molecular techniques such as the development of markers; genotyping; high density genetic linkage maps; quantitative trait loci (QTLs) identification, expressed sequence tags (ESTs), single nucleotide polymorphisms (SNPs); and transcription factors incorporated into existing breeding technologies have made possible the accurate and accelerated delivery of information for researchers. The application of genome sequencing, RNA sequencing (transcriptome sequencing), and DNA sequencing (re-sequencing) provide considerable insights for legume development and improvement programs.

Moreover, RNA-Seq helps to characterize genes, including differentially expressed genes, and can be applied for functional genomics studies, especially when there is limited information available for the studied genomes. Genome-based crop development studies and the availability of genomics data as well as decision-making gears look be specific for breeding programs. This review mainly presents an overview of the path from classical breeding to new emerging genomics tools, which will trigger and accelerate genomics-assisted breeding for recognition of novel genes for yield and quality characters for sustainable legume crop production.

Genomic changes and biochemical alterations of seed protein and oil content in a subset of fast neutron induced soybean mutants

BACKGROUND

Soybean is subjected to genetic manipulation by breeding, mutation, and transgenic approaches to produce value-added quality traits. Among those genetic approaches, mutagenesis through fast neutrons radiation is intriguing because it yields a variety of mutations, including single/multiple gene deletions and/or duplications. Characterizing the seed composition of the fast neutron mutants and its relationship with gene mutation is useful towards understanding oil and protein traits in soybean.

RESULTS

From a large population of fast neutron mutagenized plants, we selected ten mutants based on a screening of total oil and protein content using near infra-red spectroscopy. These ten mutants were regrown, and the seeds were analyzed for oil by GC-MS, protein profiling by SDS-PAGE and gene mapping by comparative genomic hybridization. The mutant 2R29C14Cladecr233cMN15 (nicknamed in this study as L10) showed higher protein and lower oil content compared to the wild type, followed by three other lines (nicknamed in this study as L03, L05, and L06). We characterized the fatty acid methyl esters profile of the trans-esterified oil and found the presence of five major fatty acids (palmitic, stearic, oleic, linoleic, and linolenic acids) at varying proportions among the mutants. Protein profile using SDS-PAGE of the ten mutants did exhibit discernable variation between storage (glycinin and β-conglycinin) and anti-nutritional factor (trypsin inhibitor) proteins. In addition, we physically mapped the position of the gene deletions or duplications in each mutant using comparative genomic hybridization.

CONCLUSION

Characterization of oil and protein profile in soybean fast neutron mutants will assist scientist and breeders to develop new value-added soybeans with improved protein and oil quality traits.

Effective identification of soybean candidate genes involved in resistance to soybean cyst nematode via direct whole genome re-sequencing of two segregating mutants

KEY MESSAGE

Three soybean candidate genes involved in resistance to soybean cyst nematode race 4 were identified via direct whole genome re-sequencing of two segregating mutants. The genes conferring resistance to soybean cyst nematode (SCN) race 4 (Hg type 1.2.3.5.7) in soybean (Glycine max L. Merr.) remains unknown. Next generation sequencing-based methods identify a wide range of targets, it is difficult to identify genes underlying traits. Use of the MutMap and QTL-seq methods to identify trait candidate genes needs backcrossing and is very time-consuming. Here we report a simple method to effectively identify candidate genes involved in resistance to SCN race 4. Two ethane methylsulfonate mutagenized mutants of soybean 'PI 437654', whose SCN race 4-infection phenotype altered, were selected. Six relevant whole genomes were re-sequenced, and then calling of genomic variants (SNPs and InDels) was conducted and compared to 'Williams 82'. The comparison eliminated many genomic variants from the mutant lines that overlapped two non-phenotypic but mutant progeny plants, wild-type PI 437654 and 'Zhonghuang 13'. Finally, only 27 mutations were found among 10 genes. Of these 10 genes, 3 genes, Glyma.09g054000, Glyma.16g065700 and Glyma.18g192200 were overlapped between two phenotypic mutant progeny plants. Therefore, the three genes may be the candidate genes involved in resistance of PI 437654 to soybean cyst nematode race 4. This method simplifies the effective identification of candidate genes.

High-Throughput Screening and Characterization of a High-Density Soybean Mutant Library Elucidate the Biosynthesis Pathway of Triterpenoid Saponins

Triterpenes (C30) constitute one of the diverse class of natural products with potential applications in food, cosmetic and pharmaceutical industries. Soyasaponins are oleanane-type triterpenoids widespread among legumes and particularly abundant in soybean seeds. They have associated with various pharmacological implications and undesirable taste properties of soybean-based food products. Uncovering the biosynthetic genes of soyasaponins will provide new opportunities to control the pathway for human benefits. However, the pathway of soyasaponin biosynthesis has not been fully elucidated in part because of a paucity of natural mutants. Here, we applied a structured high-density soybean mutant library for the forward genetic screening of triterpenoid biosynthesis. The seed soyasaponin polymorphism in the mutant library was evaluated using a high-throughput thin-layer chromatography and liquid chromatography tandem mass spectrometry analysis. This screening identified 35 mutants (3.85% of 909 mutant lines) with seven unusual soyasaponin phenotypes (Categories 1-7), which was greater than the number of natural mutants reported previously (22 mutants, 0.18% of ∼12,428 accessions). Nine unique intermediates of soyasaponin biosynthesis were identified and their chemical structures were estimated based on their MS/MS fragment patterns. Based on published information, 19 mutants could be associated with loss of function of four individual soyasaponin biosynthesis genes identified through expressed sequence tag mining or positional cloning, whereas the remaining 16 mutants were novel and may facilitate discovery of the unknown biosynthetic genes of soyasaponins. Our approach and library may help to identify new phenotype materials and causative genes associated with specialized metabolite production and other traits.

Regulation of compound leaf development in mungbean (Vigna radiata L.) by CUP-SHAPED COTYLEDON/NO APICAL MERISTEM (CUC/NAM) gene

The results provide a significant verification of functional redundancy and diversity of CUC/NAM genes in legumes. The CUP-SHAPED COTYLEDON/NO APICAL MERISTEM (CUC/NAM) orthologs play key roles for plant organ boundary formation and organ development. Here, we performed a forward screen of the gamma irradiation mutagenesis population in mungbean and characterised a mutant, reduced rachis and fused leaflets (rrf1), which gave rise to the formation of compound leaves with reduced rachis and fused leaflets. Map-based cloning revealed that RRF1 encoded a CUC/NAM protein in mungbean. Phylogenetic analysis indicated that legume CUC1/CUC2 genes were classified as belonging to two subclades, and there are different copies of CUC1/CUC2 genes in legumes. Transcriptomic analysis showed that expression levels of a set of developmental regulators, including class I KNOTTED-LIKE HOMEOBOXI (KNOXI) gene and LATERAL ORGAN BOUNDARIES DOMAIN (LBD) gene, were altered in rrf1 mutants compared to the wild-type plants. Furthermore, rrf1 genetically interacted with heptafoliate leaflets1 (hel1), a mutant displaying a seven-leaflet compound leaf, to regulate leaf development in mungbean. Our results suggest functional redundancy and diversity of two subclades of CUC1/CUC2 genes in legumes, following the duplication of an ancestral gene.

Generation and Characterization of a Foxtail Millet (Setaria italica) Mutant Library

Foxtail millet (Setaria italica) is attractive to plant scientists as a model plant because of several distinct characteristics, such as its short stature, rapid life cycle, sufficient seed production per plant, self-compatibility, true diploid nature, high photosynthetic efficiency, small genome size, and tolerance to abiotic and biotic stress. However, the study on the genetic resources of foxtail millet largely lag behind those of the other model plants such as Arabidopsis, rice and maize. Mutagenized populations cannot only create new germplasm resources, but also provide materials for gene function research. In this manuscript, an ethyl methanesulfonate (EMS)-induced foxtail millet population comprising ∼15,000 individual M₁ lines was established. Total 1353 independent lines with diverse abnormal phenotypes of leaf color, plant morphologies and panicle shapes were identified in M₂. Resequencing of sixteen randomly selected M₂ plants showed an average estimated mutation density of 1 loci/213 kb. Moreover, we provided an example for rapid cloning of the WP1 gene by a map-based cloning method. A white panicle mutant, named as wp1.a, exhibited significantly reduced chlorophyll (Chl) and carotenoid contents in leaf and panicle. Map-based cloning results showed an eight-base pair deletion located at the sixth exon of wp1.a in LOC101786849, which caused the premature termination. WP1 encoded phytoene synthase. Moreover, the sequencing analysis and cross test verified that a white panicle mutant wp1.b was an allelic mutant of wp1.a. The filed phenotypic observation and gene cloning example showed that our foxtail millet EMS-induced mutant population would be used as an important resource for functional genomics studies of foxtail millet.

Genetic control of compound leaf development in the mungbean (Vigna radiata L.)

Many studies suggest that there are distinct regulatory processes controlling compound leaf development in different clades of legumes. Loss of function of the LEAFY (LFY) orthologs results in a reduction of leaf complexity to different degrees in inverted repeat-lacking clade (IRLC) and non-IRLC species. To further understand the role of LFY orthologs and the molecular mechanism in compound leaf development in non-IRLC plants, we studied leaf development in unifoliate leaf (un) mutant, a classical mutant of mungbean (Vigna radiata L.), which showed a complete conversion of compound leaves into simple leaves. Our analysis revealed that UN encoded the mungbean LFY ortholog (VrLFY) and played a significant role in leaf development. In situ RNA hybridization results showed that STM-like KNOXI genes were expressed in compound leaf primordia in mungbean. Furthermore, increased leaflet number in heptafoliate leaflets1 (hel1) mutants was demonstrated to depend on the function of VrLFY and KNOXI genes in mungbean. Our results suggested that HEL1 is a key factor coordinating distinct processes in the control of compound leaf development in mungbean and its related non-IRLC legumes.

Identification of the dwarf gene GmDW1 in soybean (Glycine max L.) by combining mapping-by-sequencing and linkage analysis

KEY MESSAGE

GmDW1 encodes an ent-kaurene synthase (KS) acting at the early step of the biosynthesis pathway for gibberellins (GAs) and regulates the development of plant height in soybean. Plant height is an important component of plant architecture, and significantly affects crop breeding practices and yield. Here, we report the characterization of an EMS-induced dwarf mutant (dw) of the soybean cultivar Zhongpin 661 (ZDD23893). The dw mutant displayed reduced plant height and shortened internodes, both of which were mainly attributed to the longitudinally decreased cell length. The bioactive GA1 (gibberellin A1) and GA4 (gibberellin A4) were not detectable in the stem of dw, and the dwarf phenotype could be rescued by treatment with exogenous GA3. Genetic analysis showed that the dwarf trait of dw was controlled by a recessive nuclear gene. By combining linkage analysis and mapping-by-sequencing, we mapped the GmDW1 gene to an approximately 460-kb region on chromosome (Chr.) 8, containing 36 annotated genes in the reference Willliams 82 genome. Of these genes, we identified two nonsynonymous single nucleotide polymorphisms (SNPs) that are present in the encoding regions of Gmdw1 and Glyma.08G165100 in dw, respectively. However, only the SNP mutation (T>A) at nucleotide 1224 in Gmdw1 cosegregated with the dwarf phenotype. GmDW1 encodes an ent-kaurene synthase, and was expressed in various tissues including root, stem, and leaf. Further phenotypic analysis of the allelic variations in soybean accessions strongly indicated that GmDW1 is responsible for the dwarf phenotype in dw. Our results provide important information for improving our understanding of the genetics of soybean plant height and crop breeding.

Next-Generation Sequencing from Bulked-Segregant Analysis Accelerates the Simultaneous Identification of Two Qualitative Genes in Soybean

Next-generation sequencing (NGS)-based bulked-segregant analysis (BSA) approaches have been proven successful for rapidly mapping genes in plant species. However, most such methods are based on mutants and usually only one gene controlling the mutant phenotype is identified. In this study, NGS-based BSA was employed to map simultaneously two qualitative genes controlling cotyledon color of seed in soybean. Yellow-cotyledon (YC) and green-cotyledon (GC) bulks from progenies of a biparental population (Zhonghuang 30 × Jiyu 102) were sequenced. The SNP-index of each SNP locus in YC and GC bulks was calculated and two genomic regions on chromosomes 1 and 11 harboring, respectively, loci qCC1 and qCC2 were identified by Δ(SNP-index) analysis. These two BSA-seq-derived loci were further validated with SSR markers and fine-mapped. qCC1 was mapped to a 30.7-kb region containing four annotated genes and qCC2 was mapped to a 67.7-kb region with nine genes. These two regions contained, respectively, genes D1 and D2, which had previously been identified by homology-based cloning as being associated with cotyledon color. Sequence analysis of the NGS data also identified a frameshift deletion in the coding region of D1. These results suggested that BSA-seq could accelerate the mapping of loci controlling qualitative traits, even if a trait is controlled by more than one locus.