Genome resilience and prevalence of segmental duplications following fast neutron irradiation of soybean

Comparative Transcriptomics of Soybean Genotypes with Partial Resistance Toward Phytophthora sojae, Conrad, and M92-220 to Moderately Susceptible Fast Neutron Mutant Soybeans and Sloan

The breeding of disease-resistant soybeans cultivars to manage Phytophthora root and stem rot caused by the pathogen Phytophthora sojae involves combining quantitative disease resistance (QDR) and Rps gene-mediated resistance. To identify and confirm potential mechanisms of QDR toward P. sojae, we conducted a time course study comparing changes in gene expression among Conrad and M92-220 with high QDR to susceptible genotypes, Sloan, and three mutants derived from fast neutron irradiation of M92-220. Differentially expressed genes from Conrad and M92-220 indicated several shared defense-related pathways at the transcriptomic level but also defense pathways unique to each cultivar, such as stilbenoid, diarylheptanoid, and gingerol biosynthesis and monobactam biosynthesis. Gene Ontology pathway analysis showed that the susceptible fast neutron mutants lacked enrichment of three terpenoid-related pathways and two cell wall-related pathways at either one or both time points, in contrast to M92-220. The susceptible mutants also lacked enrichment of potentially important Kyoto Encyclopedia of Genes and Genomes pathways at either one or both time points, including sesquiterpenoid and triterpenoid biosynthesis; thiamine metabolism; arachidonic acid; stilbenoid, diarylheptanoid, and gingerol biosynthesis; and monobactam biosynthesis. Additionally, 31 genes that were differentially expressed in M92-220 following P. sojae infection were not expressed in the mutants. These 31 genes have annotations related to unknown proteins; valine, leucine, and isoleucine biosynthesis; and protein and lipid metabolic processes. The results of this study confirm previously proposed mechanisms of QDR, provide evidence for potential novel QDR pathways in M92-220, and further our understanding of the complex network associated with QDR mechanisms in soybean toward P. sojae.

Identification and Candidate Gene Evaluation of a Large Fast Neutron-Induced Deletion Associated with a High-Oil Phenotype in Soybean Seeds

Since the dawn of agriculture, crops have been genetically altered for desirable characteristics. This has included the selection of natural and induced mutants. Increasing the production of plant oils such as soybean (Glycine max) oil as a renewable resource for food and fuel is valuable. Successful breeding for higher oil levels in soybeans, however, usually results in reduced seed protein. A soybean fast neutron population was screened for oil content, and three high oil mutants with minimal reductions in protein levels were found. Three backcross F2 populations derived from these mutants exhibited segregation for seed oil content. DNA was pooled from the high-oil and normal-oil plants within each population and assessed by comparative genomic hybridization. A deletion encompassing 20 gene models on chromosome 14 was found to co-segregate with the high-oil trait in two of the three populations. Eighteen genes in the deleted region have known functions that appear unrelated to oil biosynthesis and accumulation pathways, while one of the unknown genes (Glyma.14G101900) may contribute to the regulation of lipid droplet formation. This high-oil trait can facilitate the breeding of high-oil soybeans without protein reduction, resulting in higher meal protein levels.

Scientific opinion on the ANSES analysis of Annex I of the EC proposal COM (2023) 411 (EFSA-Q-2024-00178)

EFSA was asked by the European Parliament to provide a scientific opinion on the analysis by the French Agency for Food, Environmental and Occupational Health & Safety (ANSES) of Annex I of the European Commission proposal for a regulation 'on plants obtained by certain new genomic techniques (NGTs) and their food and feed, and amending regulation (EU) 2017/625'. The Panel on genetically modified organisms (GMO) assessed the opinion published by ANSES, which focuses on (i) the need to clarify the definitions and scope, (ii) the scientific basis for the equivalence criteria and (iii) the need to take potential risks from category 1 NGT plants into account. The EFSA GMO Panel considered the ANSES analysis and comments on various terms used in the criteria in Annex I of the European Commission proposal and discussed definitions based on previous EFSA GMO Panel opinions. The EFSA GMO Panel concluded that the available scientific literature shows that plants containing the types and numbers of genetic modifications used as criteria to identify category 1 NGT plants in the European Commission proposal do exist as the result of spontaneous mutations or random mutagenesis. Therefore, it is scientifically justified to consider category 1 NGT plants as equivalent to conventionally bred plants with respect to the similarity of genetic modifications and the similarity of potential risks. The EFSA GMO Panel did not identify any additional hazards and risks associated with the use of NGTs compared to conventional breeding techniques in its previous Opinions.

Genomics of predictive radiation mutagenesis in oilseed rape: modifying seed oil composition

Rapeseed is a crop of global importance but there is a need to broaden the genetic diversity available to address breeding objectives. Radiation mutagenesis, supported by genomics, has the potential to supersede genome editing for both gene knockout and copy number increase, but detailed knowledge of the molecular outcomes of radiation treatment is lacking. To address this, we produced a genome re-sequenced panel of 1133 M2 generation rapeseed plants and analysed large-scale deletions, single nucleotide variants and small insertion-deletion variants affecting gene open reading frames. We show that high radiation doses (2000 Gy) are tolerated, gamma radiation and fast neutron radiation have similar impacts and that segments deleted from the genomes of some plants are inherited as additional copies by their siblings, enabling gene dosage decrease. Of relevance for species with larger genomes, we showed that these large-scale impacts can also be detected using transcriptome re-sequencing. To test the utility of the approach for predictive alteration of oil fatty acid composition, we produced lines with both decreased and increased copy numbers of Bna.FAE1 and confirmed the anticipated impacts on erucic acid content. We detected and tested a 21-base deletion expected to abolish function of Bna.FAD2.A5, for which we confirmed the predicted reduction in seed oil polyunsaturated fatty acid content. Our improved understanding of the molecular effects of radiation mutagenesis will underpin genomics-led approaches to more efficient introduction of novel genetic variation into the breeding of this crop and provides an exemplar for the predictive improvement of other crops.

Critical role for uricase and xanthine dehydrogenase in soybean nitrogen fixation and nodule development

De novo purine biosynthesis is required for the incorporation of fixed nitrogen in ureide exporting nodules, as formed on soybean [Glycine max (L.) Merr.] roots. However, in many cases, the enzymes involved in this pathway have been deduced strictly from genome annotations with little direct genetic evidence, such as mutant studies, to confirm their biochemical function or importance to nodule development. While efforts to develop large mutant collections of soybean are underway, research on this plant is still hampered by the inability to obtain mutations in any specific gene of interest. Using a forward genetic approach, as well as CRISPR/Cas9 gene editing via Agrobacterium rhizogenes-mediated hairy root transformation, we identified and characterized the role of GmUOX (Uricase) and GmXDH (Xanthine Dehydrogenase) in nitrogen fixation and nodule development in soybean. The gmuox knockout soybean mutants displayed nitrogen deficiency chlorosis and early nodule senescence, as exemplified by the reduced nitrogenase (acetylene reduction) activity in nodules, the internal greenish-white internal appearance of nodules, and diminished leghemoglobin production. In addition, gmuox1 nodules showed collapsed infected cells with degraded cytoplasm, aggregated bacteroids with no discernable symbiosome membranes, and increased formation of poly-β-hydroxybutyrate granules. Similarly, knockout gmxdh mutant nodules, generated with the CRISPR/Cas9 system, also exhibited early nodule senescence. These genetic studies confirm the critical role of the de novo purine metabolisms pathway not only in the incorporation of fixed nitrogen but also in the successful development of a functional, nitrogen-fixing nodule. Furthermore, these studies demonstrate the great utility of the CRISPR/Cas9 system for studying root-associated gene traits when coupled with hairy root transformation.

Comparison and Characterization of Phenotypic and Genomic Mutations Induced by a Carbon-Ion Beam and Gamma-ray Irradiation in Soybean (Glycine max (L.) Merr.)

Soybean (Glycine max (L.) Merr.) is a nutritious crop that can provide both oil and protein. A variety of mutagenesis methods have been proposed to obtain better soybean germplasm resources. Among the different types of physical mutagens, carbon-ion beams are considered to be highly efficient with high linear energy transfer (LET), and gamma rays have also been widely used for mutation breeding. However, systematic knowledge of the mutagenic effects of these two mutagens during development and on phenotypic and genomic mutations has not yet been elucidated in soybean. To this end, dry seeds of Williams 82 soybean were irradiated with a carbon-ion beam and gamma rays. The biological effects of the M₁ generation included changes in survival rate, yield and fertility. Compared with gamma rays, the relative biological effectiveness (RBE) of the carbon-ion beams was between 2.5 and 3.0. Furthermore, the optimal dose for soybean was determined to be 101 Gy to 115 Gy when using the carbon-ion beam, and it was 263 Gy to 343 Gy when using gamma rays. A total of 325 screened mutant families were detected from out of 2000 M₂ families using the carbon-ion beam, and 336 screened mutant families were found using gamma rays. Regarding the screened phenotypic M₂ mutations, the proportion of low-frequency phenotypic mutations was 23.4% when using a carbon ion beam, and the proportion was 9.8% when using gamma rays. Low-frequency phenotypic mutations were easily obtained with the carbon-ion beam. After screening the mutations from the M₂ generation, their stability was verified, and the genome mutation spectrum of M₃ was systemically profiled. A variety of mutations, including single-base substitutions (SBSs), insertion-deletion mutations (INDELs), multinucleotide variants (MNVs) and structural variants (SVs) were detected with both carbon-ion beam irradiation and gamma-ray irradiation. Overall, 1988 homozygous mutations and 9695 homozygous + heterozygous genotype mutations were detected when using the carbon-ion beam. Additionally, 5279 homozygous mutations and 14,243 homozygous + heterozygous genotype mutations were detected when using gamma rays. The carbon-ion beam, which resulted in low levels of background mutations, has the potential to alleviate the problems caused by linkage drag in soybean mutation breeding. Regarding the genomic mutations, when using the carbon-ion beam, the proportion of homozygous-genotype SVs was 0.45%, and that of homozygous + heterozygous-genotype SVs was 6.27%; meanwhile, the proportions were 0.04% and 4.04% when using gamma rays. A higher proportion of SVs were detected when using the carbon ion beam. The gene effects of missense mutations were greater under carbon-ion beam irradiation, and the gene effects of nonsense mutations were greater under gamma-ray irradiation, which meant that the changes in the amino acid sequences were different between the carbon-ion beam and gamma rays. Taken together, our results demonstrate that both carbon-ion beam and gamma rays are effective techniques for rapid mutation breeding in soybean. If one would like to obtain mutations with a low-frequency phenotype, low levels of background genomic mutations and mutations with a higher proportion of SVs, carbon-ion beams are the best choice.

A chromosome 16 deletion conferring a high sucrose phenotype in soybean

Sucrose in soybean seeds is desirable for many end-uses. Increased sucrose contents were discovered to associate with a chromosome 16 deletion resulting from fast neutron irradiation. Soybean is one of the most economically important crops in the United States. A primary end-use of soybean is for livestock feed. Therefore, genetic improvement of seed composition is one of the most important goals in soybean breeding programs. Sucrose is desired in animal feed due to its role as an easily digestible energy source. An elite soybean line was irradiated with fast neutrons and the seed from plants were screened for altered seed composition with near-infrared spectroscopy (NIR). One mutant line, G15FN-54, was found to have higher sucrose content (8-9%) than the parental line (5-6%). Comparative genomic hybridization (CGH) revealed three large deletions on chromosomes (Chrs) 10, 13, and 16 in the mutant, which were confirmed through whole genome sequencing (WGS). A bi-parental population derived from the mutant G15FN-54 and the cultivar Benning was developed to conduct a bulked segregant analysis (BSA) with SoySNP50K BeadChips, revealing that the deletion on Chr 16 might be responsible for the altered phenotype. The mapping result using the bi-parental population confirmed that the deletion on Chr 16 conferred elevated sucrose content and a total of 21 genes are located within this Chr 16 deletion. NIR and high-pressure liquid chromatography (HPLC) were used to confirm the stability of the phenotype across generations in the bi-parental population. The mutation will be useful to understand the genetic control of soybean seed sucrose content.

Genetic mapping and functional genomics of soybean seed protein

Soybean is an utterly important crop for high-quality meal protein and vegetative oil. Soybean seed protein content has become a key factor in nutrients for livestock feed as well as human dietary consumption. Genetic improvement of soybean seed protein is highly desired to meet the demands of rapidly growing world population. Molecular mapping and genomic analysis in soybean have identified many quantitative trait loci (QTL) underlying seed protein content control. Exploring the mechanisms of seed storage protein regulation will be helpful to achieve the improvement of protein content. However, the practice of breeding higher protein soybean is challenging because soybean seed protein is negatively correlated with seed oil content and yield. To overcome the limitation of such inverse relationship, deeper insights into the property and genetic control of seed protein are required. Recent advances of soybean genomics have strongly enhanced the understandings for molecular mechanisms of soybean with better seed quality. Here, we review the research progress in the genetic characteristics of soybean storage protein, and up-to-date advances of molecular mappings and genomics of soybean protein. The key factors underlying the mechanisms of the negative correlation between protein and oil in soybean seeds are elaborated. We also briefly discuss the future prospects of breaking the bottleneck of the negative correlation to develop high protein soybean without penalty of oil and yield.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01373-5.

Flagellin and mannitol modulate callose biosynthesis and deposition in soybean seedlings

Callose is a polymer deposited on the cell wall and is necessary for plant growth and development. Callose is synthesized by genes from the glucan synthase-like family (GSL) and dynamically responds to various types of stress. Callose can inhibit pathogenic infection, in the case of biotic stresses, and maintain cell turgor and stiffen the plant cell wall in abiotic stresses. Here, we report the identification of 23 GSL genes (GmGSL) in the soybean genome. We performed phylogenetic analyses, gene structure prediction, duplication patterns, and expression profiles on several RNA-Seq libraries. Our analyses show that WGD/Segmental duplication contributed to expanding this gene family in soybean. Next, we analyzed the callose responses in soybean under abiotic and biotic stresses. The data show that callose is induced by both osmotic stress and flagellin 22 (flg22) and is related to the activity of β-1,3-glucanases. By using RT-qPCR, we evaluated the expression of GSL genes during the treatment of soybean roots with mannitol and flg22. The GmGSL23 gene was upregulated in seedlings treated with osmotic stress or flg22, showing the essential role of this gene in the soybean defense response to pathogenic organisms and osmotic stress. Our results provide an important understanding of the role of callose deposition and regulation of GSL genes in response to osmotic stress and flg22 infection in soybean seedlings.

Irradiation with carbon ion beams affects soybean nutritional quality in early generations

As people's demand for healthy diet increases, improving soybean seed nutritional quality is becoming as important as yield. Carbon ion beam radiation (CIBR) is an effective method to create soybean mutants, and thus breeding cultivars with better seed nutritional quality. In this study, the high-yield soybean line 'Dongsheng 28' was used, and three CIBR doses (100, 120, and 140 Gy) were used to explore the characteristics of quality separation and variation in the offspring of early mutant populations. Eleven quality traits, including protein, oil, sucrose, soluble sugar, iron (Fe), manganese (Mn), zinc (Zn), cupper (Cu), daidzin, glycitin, and genistin concentrations were analyzed in the M₂ and M₃ generations. The results revealed that the range of protein and oil concentration of all three CIBR doses changed by 38.5-42.9% and 18.8-23.8% in the M₂ and M₃ generations, respectively, while soluble sugar and sucrose concentrations changed by 48.1-123.4 and 22.7-74.7 mg/g, with significant effects by 140 Gy across the two generations. Therefore, around the optimum range, a higher CIBR dose is better for high protein, oil, and sugar varieties selection. In general, irradiation raised isoflavone concentrations, but 140 Gy had an inhibitory effect on isoflavone concentrations in the M₃ generation. Although a variety could not be released in the M₂ or M₃ generation, the results of this study have important guiding significance for the targeted cultivation of specific nutritional quality materials. For instance, a lower irradiation dose is preferable when breeding targets are higher isoflavones and Mn concentrations. It is essential to increase the irradiation dose if the breeding targets contain high levels of protein, oil, sucrose, soluble sugars, Fe, Zn, and Cu.

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.

Fast neutron mutagenesis in soybean enriches for small indels and creates frameshift mutations

The mutagenic effects of ionizing radiation have been used for decades to create novel variants in experimental populations. Fast neutron (FN) bombardment as a mutagen has been especially widespread in plants, with extensive reports describing the induction of large structural variants, i.e., deletions, insertions, inversions, and translocations. However, the full spectrum of FN-induced mutations is poorly understood. We contrast small insertions and deletions (indels) observed in 27 soybean lines subject to FN irradiation with the standing indels identified in 107 diverse soybean lines. We use the same populations to contrast the nature and context (bases flanking a nucleotide change) of single-nucleotide variants. The accumulation of new single-nucleotide changes in FN lines is marginally higher than expected based on spontaneous mutation. In FN-treated lines and in standing variation, C→T transitions and the corresponding reverse complement G→A transitions are the most abundant and occur most frequently in a CpG local context. These data indicate that most SNPs identified in FN lines are likely derived from spontaneous de novo processes in generations following mutagenesis rather than from the FN irradiation mutagen. However, small indels in FN lines differ from standing variants. Short insertions, from 1 to 6 bp, are less abundant than in standing variation. Short deletions are more abundant and prone to induce frameshift mutations that should disrupt the structure and function of encoded proteins. These findings indicate that FN irradiation generates numerous small indels, increasing the abundance of loss-of-function mutations that impact single genes.

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

Plant Genome Editing and the Relevance of Off-Target Changes

Site-directed nucleases (SDNs) used for targeted genome editing are powerful new tools to introduce precise genetic changes into plants. Like traditional approaches, such as conventional crossing and induced mutagenesis, genome editing aims to improve crop yield and nutrition. Next-generation sequencing studies demonstrate that across their genomes, populations of crop species typically carry millions of single nucleotide polymorphisms and many copy number and structural variants. Spontaneous mutations occur at rates of ∼10-8 to 10-9 per site per generation, while variation induced by chemical treatment or ionizing radiation results in higher mutation rates. In the context of SDNs, an off-target change or edit is an unintended, nonspecific mutation occurring at a site with sequence similarity to the targeted edit region. SDN-mediated off-target changes can contribute to a small number of additional genetic variants compared to those that occur naturally in breeding populations or are introduced by induced-mutagenesis methods. Recent studies show that using computational algorithms to design genome editing reagents can mitigate off-target edits in plants. Finally, crops are subject to strong selection to eliminate off-type plants through well-established multigenerational breeding, selection, and commercial variety development practices. Within this context, off-target edits in crops present no new safety concerns compared to other breeding practices. The current generation of genome editing technologies is already proving useful to develop new plant varieties with consumer and farmer benefits. Genome editing will likely undergo improved editing specificity along with new developments in SDN delivery and increasing genomic characterization, further improving reagent design and application.

Transcriptome sequencing reveals hotspot mutation regions and dwarfing mechanisms in wheat mutants induced by γ-ray irradiation and EMS

Induced mutation is an important approach for creating novel plant germplasms. The introduction of dwarf or semi-dwarf genes into wheat has led to great advancements in yield improvement. In this study, four elite dwarf wheat mutants, named dm1-dm4, induced from γ-ray irradiation or ethyl methanesulfonate (EMS) mutagenesis, were used to identify transcriptome variations and dwarfing mechanisms. The results showed that the hotspot regions of mutations distributed on the chromosomes were consistent among the four mutant lines and these regions were mainly located around the 50, 360 and 400 Mb positions of chromosome 1A and the distal regions of chromosomes 2A and 2BL. Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that 'protein processing in endoplasmic reticulum' was the most common significantly enriched pathway based on the differentially expressed genes (DEGs) between wildtype (WT) and the mutants. Notably, 18 out of 20 genes involved in this process encode heat shock proteins (HSPs). The results implied that HSPs might participate in wheat dwarfism response and function in the dwarfism process through protein folding and/or degradation. Moreover, seven genes in dm4 involved in modulating auxin levels were down-regulated and dm4 was more sensitive to auxin treatment compared with WT, indicating the important roles of auxin in regulation of dwarf phenotype in dm4. This study not only identified transcriptome sequence variation induced by physical and chemical mutagenesis but also revealed potential dwarfing mechanisms in the wheat mutant lines.

On the Inverse Correlation of Protein and Oil: Examining the Effects of Altered Central Carbon Metabolism on Seed Composition Using Soybean Fast Neutron Mutants

Protein and oil levels measured at maturity are inversely correlated across soybean lines; however, carbon is in limited supply during maturation resulting in tradeoffs for the production of other reserves including oligosaccharides. During the late stages of seed development, the allocation of carbon for storage reserves changes. Lipid and protein levels decline while concentrations of indigestible raffinose family oligosaccharides (RFOs) increase, leading to a decreased crop value. Since the maternal source of carbon is diminished during seed maturation stages of development, carbon supplied to RFO synthesis likely comes from an internal, turned-over source and may contribute to the reduction in protein and lipid content in mature seeds. In this study, fast neutron (FN) mutagenized soybean populations with deletions in central carbon metabolic genes were examined for trends in oil, protein, sugar, and RFO accumulation leading to an altered final composition. Two lines with concurrent increases in oil and protein, by combined 10%, were identified. A delayed switch in carbon allocation towards RFO biosynthesis resulted in extended lipid accumulation and without compromising protein. Strategies for future soybean improvement using FN resources are described.

Identification and characterization of a fast-neutron-induced mutant with elevated seed protein content in soybean

KEY MESSAGE

Protein content of soybean is critical for utility of soybean meal. A fast-neutron-induced deletion on chromosome 12 was found to be associated with increased protein content. Soybean seed composition affects the utility of soybean, and improving seed composition is an essential breeding goal. Fast neutron radiation introduces genomic mutations resulting in novel variation for traits of interest. Two elite soybean lines were irradiated with fast neutrons and screened for altered seed composition. Twenty-three lines with altered protein, oil, or sucrose content were selected based on near-infrared spectroscopy data from five environments and yield tested at five locations. Mutants with significantly increased protein averaged 19.1-36.8 g kg-1 more protein than the parents across 10 environments. Comparative genomic hybridization (CGH) identified putative mutations in a mutant, G15FN-12, that has 36.8 g kg-1 higher protein than the parent genotype, and whole genome sequencing (WGS) of the mutant has confirmed these mutations. An F2:3 population was developed from G15FN-12 to determine association between genomic changes and increased protein content. Bulked segregant analysis of the population using the SoySNP50K BeadChip identified a CGH- and WGS-confirmed deletion on chromosome 12 to be responsible for elevated protein content. The population was genotyped using a KASP marker designed at the mutation region, and significant association (P < 0.0001) between the deletion on chromosome 12 and elevated protein content was observed and confirmed in the F3:4 generation. The F2 segregants homozygous for the deletion averaged 27 g kg-1 higher seed protein and 8 g kg-1 lower oil than homozygous wild-type segregants. Mutants with altered seed composition are a new resource for gene function studies and provide elite materials for genetic improvement of seed composition.

Functional analysis and development of a CRISPR/Cas9 allelic series for a CPR5 ortholog necessary for proper growth of soybean trichomes

Developments in genomic and genome editing technologies have facilitated the mapping, cloning, and validation of genetic variants underlying trait variation. This study combined bulked-segregant analysis, array comparative genomic hybridization, and CRISPR/Cas9 methodologies to identify a CPR5 ortholog essential for proper trichome growth in soybean (Glycine max). A fast neutron mutant line exhibited short trichomes with smaller trichome nuclei compared to its parent line. A fast neutron-induced deletion was identified within an interval on chromosome 6 that co-segregated with the trichome phenotype. The deletion encompassed six gene models including an ortholog of Arabidopsis thaliana CPR5. CRISPR/Cas9 was used to mutate the CPR5 ortholog, resulting in five plants harboring a total of four different putative knockout alleles and two in-frame alleles. Phenotypic analysis of the mutants validated the candidate gene, and included intermediate phenotypes that co-segregated with the in-frame alleles. These findings demonstrate that the CPR5 ortholog is essential for proper growth and development of soybean trichomes, similar to observations in A. thaliana. Furthermore, this work demonstrates the value of using CRISPR/Cas9 to generate an allelic series and intermediate phenotypes for functional analysis of candidate genes and/or the development of novel traits.

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.

DNA Copy Number Variations as Markers of Mutagenic Impact

DNA copy number variation (CNV) occurs due to deletion or duplication of DNA segments resulting in a different number of copies of a specific DNA-stretch on homologous chromosomes. Implications of CNVs in evolution and development of different diseases have been demonstrated although contribution of environmental factors, such as mutagens, in the origin of CNVs, is poorly understood. In this review, we summarize current knowledge about mutagen-induced CNVs in human, animal and plant cells. Differences in CNV frequencies induced by radiation and chemical mutagens, distribution of CNVs in the genome, as well as adaptive effects in plants, are discussed. Currently available information concerning impact of mutagens in induction of CNVs in germ cells is presented. Moreover, the potential of CNVs as a new endpoint in mutagenicity test-systems is discussed.

Screening populations for copy number variation using genotyping-by-sequencing: a proof of concept using soybean fast neutron mutants

BACKGROUND

The effective use of mutant populations for reverse genetic screens relies on the population-wide characterization of the induced mutations. Genome- and population-wide characterization of the mutations found in fast neutron populations has been hindered, however, by the wide range of mutations generated and the lack of affordable technologies to detect DNA sequence changes. In this study, we therefore aimed to test whether genotyping-by-sequencing (GBS) technology could be used to characterize copy number variation (CNV) induced by fast neutrons in a soybean mutant population.

RESULTS

We called CNVs from GBS data in 79 soybean mutants and assessed the sensitivity and precision of this approach by validating our results against array comparative genomic hybridization (aCGH) data for 19 of these mutants as well as targeted PCR and ddPCR assays for a representative subset of the smallest events detected by GBS. Our GBS pipeline detected 55 of the 96 events found by aCGH, with approximate detection thresholds of 60 kb, 500 kb and 1 Mb for homozygous deletions, hemizygous deletions and duplications, respectively. Among the whole set of 79 mutants, the GBS data revealed 105 homozygous deletions, 32 hemizygous deletions and 19 duplications. This included several extremely large events, exhibiting maximum sizes of ~ 11.2 Mb for a homozygous deletion, ~ 11.6 Mb for a hemizygous deletion, and ~ 50 Mb for a duplication.

CONCLUSIONS

This study provides a proof of concept that GBS can be used as an affordable high-throughput method for assessing CNVs in fast neutron mutants. The modularity of this GBS approach allows combining as many different libraries or sequencing runs as is necessary for reaching the goals of a particular study. This method should enable the low-cost genome-wide characterization of hundreds to thousands of individuals in fast neutron mutant populations or any population with large genomic deletions and duplications.

Whole-Genome Sequencing Identifies a Rice Grain Shape Mutant, gs9-1

BACKGROUND

Breeding for genes controlling key agronomic traits is an important goal of rice genetic improvement. To gain insight into genes controlling grain morphology, we screened M3 plants derived from 1,000 whole-genome sequenced (WGS) M2 Kitaake mutants to identify lines with altered grain size.

RESULTS

In this study, we isolated a mutant, named fast-neutron (FN) 60-4, which exhibits a significant reduction in grain size. We crossed FN60-4 with the parental line Kitaake and analyzed the resulting backcross population. Segregation analysis of 113 lines from the BC2F2 population revealed that the mutant phenotype is controlled by a single semi-dominant locus. Mutant FN60-4 is reduced 20% in plant height and 8.8% in 1000-grain weight compared with Kitaake. FN60-4 also exhibits an 8% reduction in cell number and a 9% reduction in cell length along the vertical axis of the glume. We carried out whole-genome sequencing of DNA pools extracted from segregants with long grains or short grains, and revealed that one gene, LOC_Os09g02650, cosegregated with the grain size phenotype in the BC1F2 and BC2F2 populations. This mutant allele was named grain shape 9-1 (gs9-1). gs9-1 carries a 3-bp deletion that affects two amino acids. This locus is a new allele of the BC12/GDD1/MTD1 gene that encodes a kinesin-like protein involved in cell-cycle progression, cellulose microfibril deposition and gibberellic acid (GA) biosynthesis. The GA biosynthesis-related gene KO2 is down-regulated in gs9-1. The dwarf phenotype of gs9-1 can be rescued by adding exogenous GA3. In contrast to the phenotypes for the other alleles, the gs9-1 is less severe, consistent with the nature of the mutation, which does not disrupt the open reading frame as observed for the other alleles.

CONCLUSIONS

In this study, we isolated a mutant, which exhibits altered grain shape and identified the mutated gene, gs9-1. Our study reveals that gs9-1 is a semi-dominant gene that carries a two-amino acid mutation. gs9-1 is allelic to the BC12/GDD1/MTD1 gene involved in GA biosynthesis. These results demonstrate the efficiency and convenience of cloning genes from the whole-genome sequenced Kitaake mutant population to advance investigations into genes controlling key agronomic traits in rice.

Expanding Avenue of Fast Neutron Mediated Mutagenesis for Crop Improvement

: Fast neutron (FN) radiation mediated mutagenesis is a unique approach among the several induced mutagenesis methods being used in plant science in terms of impacted mutations. The FN mutagenesis usually creates deletions from few bases to several million bases (Mb). A library of random deletion generated using FN mutagenesis lines can provide indispensable resources for the reverse genetic approaches. In this review, information from several efforts made using FN mutagenesis has been compiled to understand the type of induced mutations, frequency, and genetic stability. Concerns regarding the utilization of FN mutagenesis technique for a plant with different level of ploidy and genome complexity are discussed. We have highlighted the utility of next-generation sequencing techniques that can be efficiently utilized for the characterization of mutant lines as well as for the mapping of causal mutations. Pros and cons of mapping by mutation (MutMap), mutant chromosome sequencing (MutChromSeq), exon capture, whole genome sequencing, MutRen-Seq, and different tilling approaches that can be used for the detection of FN-induced mutation has also been discussed. Genomic resources developed using the FN mutagenesis have been catalogued wooing to meaningful utilization of the available resources. The information provided here will be helpful for the efficient exploration for the crop improvement programs and for better understanding of genetic regulations.

Molecular mapping and genomics of soybean seed protein: a review and perspective for the future

KEY MESSAGE

Genetic improvement of soybean protein meal is a complex process because of negative correlation with oil, yield, and temperature. This review describes the progress in mapping and genomics, identifies knowledge gaps, and highlights the need of integrated approaches. Meal protein derived from soybean [Glycine max (L) Merr.] seed is the primary source of protein in poultry and livestock feed. Protein is a key factor that determines the nutritional and economical value of soybean. Genetic improvement of soybean seed protein content is highly desirable, and major quantitative trait loci (QTL) for soybean protein have been detected and repeatedly mapped on chromosomes (Chr.) 20 (LG-I), and 15 (LG-E). However, practical breeding progress is challenging because of seed protein content's negative genetic correlation with seed yield, other seed components such as oil and sucrose, and interaction with environmental effects such as temperature during seed development. In this review, we discuss rate-limiting factors related to soybean protein content and nutritional quality, and potential control factors regulating seed storage protein. In addition, we describe advances in next-generation sequencing technologies for precise detection of natural variants and their integration with conventional and high-throughput genotyping technologies. A syntenic analysis of QTL on Chr. 15 and 20 was performed. Finally, we discuss comprehensive approaches for integrating protein and amino acid QTL, genome-wide association studies, whole-genome resequencing, and transcriptome data to accelerate identification of genomic hot spots for allele introgression and soybean meal protein improvement.

The Sequences of 1504 Mutants in the Model Rice Variety Kitaake Facilitate Rapid Functional Genomic Studies

The availability of a whole-genome sequenced mutant population and the cataloging of mutations of each line at a single-nucleotide resolution facilitate functional genomic analysis. To this end, we generated and sequenced a fast-neutron-induced mutant population in the model rice cultivar Kitaake (Oryza sativa ssp japonica), which completes its life cycle in 9 weeks. We sequenced 1504 mutant lines at 45-fold coverage and identified 91,513 mutations affecting 32,307 genes, i.e., 58% of all rice genes. We detected an average of 61 mutations per line. Mutation types include single-base substitutions, deletions, insertions, inversions, translocations, and tandem duplications. We observed a high proportion of loss-of-function mutations. We identified an inversion affecting a single gene as the causative mutation for the short-grain phenotype in one mutant line. This result reveals the usefulness of the resource for efficient, cost-effective identification of genes conferring specific phenotypes. To facilitate public access to this genetic resource, we established an open access database called KitBase that provides access to sequence data and seed stocks. This population complements other available mutant collections and gene-editing technologies. This work demonstrates how inexpensive next-generation sequencing can be applied to generate a high-density catalog of mutations.

An Induced Chromosomal Translocation in Soybean Disrupts a KASI Ortholog and Is Associated with a High-Sucrose and Low-Oil Seed Phenotype

Mutagenesis is a useful tool in many crop species to induce heritable genetic variability for trait improvement and gene discovery. In this study, forward screening of a soybean fast neutron (FN) mutant population identified an individual that produced seed with nearly twice the amount of sucrose (8.1% on dry matter basis) and less than half the amount of oil (8.5% on dry matter basis) as compared to wild type. Bulked segregant analysis (BSA), comparative genomic hybridization, and genome resequencing were used to associate the seed composition phenotype with a reciprocal translocation between chromosomes 8 and 13. In a backcross population, the translocation perfectly cosegregated with the seed composition phenotype and exhibited non-Mendelian segregation patterns. We hypothesize that the translocation is responsible for the altered seed composition by disrupting a β-ketoacyl-[acyl carrier protein] synthase 1 (KASI) ortholog. KASI is a core fatty acid synthesis enzyme that is involved in the conversion of sucrose into oil in developing seeds. This finding may lead to new research directions for developing soybean cultivars with modified carbohydrate and oil seed composition.

Identification of Substitutions and Small Insertion-Deletions Induced by Carbon-Ion Beam Irradiation in Arabidopsis thaliana

Heavy-ion beam irradiation is one of the principal methods used to create mutants in plants. Research on mutagenic effects and molecular mechanisms of radiation is an important subject that is multi-disciplinary. Here, we re-sequenced 11 mutagenesis progeny (M3) Arabidopsis thaliana lines derived from carbon-ion beam (CIB) irradiation, and subsequently focused on substitutions and small insertion-deletion (INDELs). We found that CIB induced more substitutions (320) than INDELs (124). Meanwhile, the single base INDELs were more prevalent than those in large size (≥2 bp). In details, the detected substitutions showed an obvious bias of C > T transitions, by activating the formation of covalent linkages between neighboring pyrimidine residues in the DNA sequence. An A and T bias was observed among the single base INDELs, in which most of these were induced by replication slippage at either the homopolymer or polynucleotide repeat regions. The mutation rate of 200-Gy CIB irradiation was estimated as 3.37 × 10^-7 per site. Different from previous researches which mainly focused on the phenotype, chromosome aberration, genetic polymorphism, or sequencing analysis of specific genes only, our study revealed genome-wide molecular profile and rate of mutations induced by CIB irradiation. We hope our data could provide valuable clues for explaining the potential mechanism of plant mutation breeding by CIB irradiation.

Identification of Homogentisate Dioxygenase as a Target for Vitamin E Biofortification in Oilseeds

Soybean (Glycine max) is a major plant source of protein and oil and produces important secondary metabolites beneficial for human health. As a tool for gene function discovery and improvement of this important crop, a mutant population was generated using fast neutron irradiation. Visual screening of mutagenized seeds identified a mutant line, designated MO12, which produced brown seeds as opposed to the yellow seeds produced by the unmodified Williams 82 parental cultivar. Using forward genetic methods combined with comparative genome hybridization analysis, we were able to establish that deletion of the GmHGO1 gene is the genetic basis of the brown seeded phenotype exhibited by the MO12 mutant line. GmHGO1 encodes a homogentisate dioxygenase (HGO), which catalyzes the committed enzymatic step in homogentisate catabolism. This report describes to our knowledge the first functional characterization of a plant HGO gene, defects of which are linked to the human genetic disease alkaptonuria. We show that reduced homogentisate catabolism in a soybean HGO mutant is an effective strategy for enhancing the production of lipid-soluble antioxidants such as vitamin E, as well as tolerance to herbicides that target pathways associated with homogentisate metabolism. Furthermore, this work demonstrates the utility of fast neutron mutagenesis in identifying novel genes that contribute to soybean agronomic traits.

Fast neutron-induced structural rearrangements at a soybean NAP1 locus result in gnarled trichomes

KEY MESSAGE

Three adjacent and distinct sequence rearrangements were identified at a NAP1 locus in a soybean mutant. Genetic dissection and validation revealed the function of this gene in soybean trichome development. A soybean (Glycine max (L.) Merr.) gnarled trichome mutant, exhibiting stunted trichomes compared to wild-type, was identified in a fast neutron mutant population. Genetic mapping using whole genome sequencing-based bulked segregant analysis identified a 26.6 megabase interval on chromosome 20 that co-segregated with the phenotype. Comparative genomic hybridization analysis of the mutant indicated that the chromosome 20 interval included a small structural variant within the coding region of a soybean ortholog (Glyma.20G019300) of Arabidopsis Nck-Associated Protein 1 (NAP1), a regulator of actin nucleation during trichome morphogenesis. Sequence analysis of the candidate allele revealed multiple rearrangements within the coding region, including two deletions (approximately 1-2 kb each), a translocation, and an inversion. Further analyses revealed that the mutant allele perfectly co-segregated with the phenotype, and a wild-type soybean NAP1 transgene functionally complemented an Arabidopsis nap1 mutant. In addition, mapping and exon sequencing of NAP1 in a spontaneous soybean gnarled trichome mutant (T31) identified a frame shift mutation resulting in a truncation of the coding region. These data indicate that the soybean NAP1 gene is essential for proper trichome development and show the utility of the soybean fast neutron population for forward genetic approaches for identifying genes.

Genome engineering and plant breeding: impact on trait discovery and development

KEY MESSAGE

New tools for the precise modification of crops genes are now available for the engineering of new ideotypes. A future challenge in this emerging field of genome engineering is to develop efficient methods for allele mining. Genome engineering tools are now available in plants, including major crops, to modify in a predictable manner a given gene. These new techniques have a tremendous potential for a spectacular acceleration of the plant breeding process. Here, we discuss how genetic diversity has always been the raw material for breeders and how they have always taken advantage of the best available science to use, and when possible, increase, this genetic diversity. We will present why the advent of these new techniques gives to the breeders extremely powerful tools for crop breeding, but also why this will require the breeders and researchers to characterize the genes underlying this genetic diversity more precisely. Tackling these challenges should permit the engineering of optimized alleles assortments in an unprecedented and controlled way.

Soybean (Glycine max) Mutant and Germplasm Resources: Current Status and Future Prospects

Genetic bottlenecks during domestication and modern breeding limited the genetic diversity of soybean (Glycine max (L.) Merr.). Therefore, expanding and diversifying soybean genetic resources is a major priority for the research community. These resources, consisting of natural and induced genetic variants, are valuable tools for improving soybean and furthering soybean biological knowledge. During the twentieth century, researchers gathered a wealth of genetic variation in the forms of landraces, Glycine soja accessions, Glycine tertiary germplasm, and the U.S. Department of Agriculture (USDA) Type and Isoline Collections. During the twenty-first century, soybean researchers have added several new genetic and genomic resources. These include the reference genome sequence, genotype data for the USDA soybean germplasm collection, next-generation mapping populations, new irradiation and transposon-based mutagenesis populations, and designer nuclease platforms for genome engineering. This paper briefly surveys the publicly accessible soybean genetic resources currently available or in development and provides recommendations for developing such genetic resources in the future. © 2016 by John Wiley & Sons, Inc.

Fighting Asian Soybean Rust

Phakopsora pachyrhizi is a biotrophic fungus provoking SBR disease. SBR poses a major threat to global soybean production. Though several R genes provided soybean immunity to certain P. pachyrhizi races, the pathogen swiftly overcame this resistance. Therefore, fungicides are the only current means to control SBR. However, insensitivity to fungicides is soaring in P. pachyrhizi and, therefore, alternative measures are needed for SBR control. In this article, we discuss the different approaches for fighting SBR and their potential, disadvantages, and advantages over other measures. These encompass conventional breeding for SBR resistance, transgenic approaches, exploitation of transcription factors, secondary metabolites, and antimicrobial peptides, RNAi/HIGS, and biocontrol strategies. It seems that an integrating approach exploiting different measures is likely to provide the best possible means for the effective control of SBR.

Genomic variation and DNA repair associated with soybean transgenesis: a comparison to cultivars and mutagenized plants

BACKGROUND

The safety of mutagenized and genetically transformed plants remains a subject of scrutiny. Data gathered and communicated on the phenotypic and molecular variation induced by gene transfer technologies will provide a scientific-based means to rationally address such concerns. In this study, genomic structural variation (e.g. large deletions and duplications) and single nucleotide polymorphism rates were assessed among a sample of soybean cultivars, fast neutron-derived mutants, and five genetically transformed plants developed through Agrobacterium based transformation methods.

RESULTS

On average, the number of genes affected by structural variations in transgenic plants was one order of magnitude less than that of fast neutron mutants and two orders of magnitude less than the rates observed between cultivars. Structural variants in transgenic plants, while rare, occurred adjacent to the transgenes, and at unlinked loci on different chromosomes. DNA repair junctions at both transgenic and unlinked sites were consistent with sequence microhomology across breakpoints. The single nucleotide substitution rates were modest in both fast neutron and transformed plants, exhibiting fewer than 100 substitutions genome-wide, while inter-cultivar comparisons identified over one-million single nucleotide polymorphisms.

CONCLUSIONS

Overall, these patterns provide a fresh perspective on the genomic variation associated with high-energy induced mutagenesis and genetically transformed plants. The genetic transformation process infrequently results in novel genetic variation and these rare events are analogous to genetic variants occurring spontaneously, already present in the existing germplasm, or induced through other types of mutagenesis. It remains unclear how broadly these results can be applied to other crops or transformation methods.

MicroRNA Maturation and MicroRNA Target Gene Expression Regulation Are Severely Disrupted in Soybean dicer-like1 Double Mutants

Small nonprotein-coding microRNAs (miRNAs) are present in most eukaryotes and are central effectors of RNA silencing-mediated mechanisms for gene expression regulation. In plants, DICER-LIKE1 (DCL1) is the founding member of a highly conserved family of RNase III-like endonucleases that function as core machinery proteins to process hairpin-like precursor transcripts into mature miRNAs, small regulatory RNAs, 21-22 nucleotides in length. Zinc finger nucleases (ZFNs) were used to generate single and double-mutants of putative soybean DCL1 homologs, DCL1a and DCL1b, to confirm their functional role(s) in the soybean miRNA pathway. Neither DCL1 single mutant, dcl1a or dcl1b plants, exhibited a pronounced morphological or molecular phenotype. However, the dcl1a/dcl1b double mutant expressed a strong morphological phenotype, characterized by reduced seed size and aborted seedling development, in addition to defective miRNA precursor transcript processing efficiency and deregulated miRNA target gene expression. Together, these findings indicate that the two soybean DCL1 paralogs, DCL1a and DCL1b, largely play functionally redundant roles in the miRNA pathway and are essential for normal plant development.

Functional analysis of duplicated Symbiosis Receptor Kinase (SymRK) genes during nodulation and mycorrhizal infection in soybean (Glycine max)

Association between legumes and rhizobia results in the formation of root nodules, where symbiotic nitrogen fixation occurs. The early stages of this association involve a complex of signalling events between the host and microsymbiont. Several genes dealing with early signal transduction have been cloned, and one of them encodes the leucine-rich repeat (LRR) receptor kinase (SymRK; also termed NORK). The Symbiosis Receptor Kinase gene is required by legumes to establish a root endosymbiosis with Rhizobium bacteria as well as mycorrhizal fungi. Using degenerate primer and BAC sequencing, we cloned duplicated SymRK homeologues in soybean called GmSymRKα and GmSymRKβ. These duplicated genes have high similarity of nucleotide (96%) and amino acid sequence (95%). Sequence analysis predicted a malectin-like domain within the extracellular domain of both genes. Several putative cis-acting elements were found in promoter regions of GmSymRKα and GmSymRKβ, suggesting a participation in lateral root development, cell division and peribacteroid membrane formation. The mutant of SymRK genes is not available in soybean; therefore, to know the functions of these genes, RNA interference (RNAi) of these duplicated genes was performed. For this purpose, RNAi construct of each gene was generated and introduced into the soybean genome by Agrobacterium rhizogenes-mediated hairy root transformation. RNAi of GmSymRKβ gene resulted in an increased reduction of nodulation and mycorrhizal infection than RNAi of GmSymRKα, suggesting it has the major activity of the duplicated gene pair. The results from the important crop legume soybean confirm the joint phenotypic action of GmSymRK genes in both mycorrhizal and rhizobial infection seen in model legumes.