Genome-wide association study reveals GmFulb as candidate gene for maturity time and reproductive length in soybeans (Glycine max).

The ability of soybean [Glycine max (L.) Merr.] to adapt to different latitudes is attributed to genetic variation in major E genes and quantitative trait loci (QTLs) determining flowering time (R1), maturity (R8), and reproductive length (RL). Fully revealing the genetic basis of R1, R8, and RL in soybeans is necessary to enhance genetic gains in soybean yield improvement. Here, we performed a genome-wide association analysis (GWA) with 31,689 single nucleotide polymorphisms (SNPs) to detect novel loci for R1, R8, and RL using a soybean panel of 329 accessions with the same genotype for three major E genes (e1-as/E2/E3). The studied accessions were grown in nine environments and observed for R1, R8 and RL in all environments. This study identified two stable peaks on Chr 4, simultaneously controlling R8 and RL. In addition, we identified a third peak on Chr 10 controlling R1. Association peaks overlap with previously reported QTLs for R1, R8, and RL. Considering the alternative alleles, significant SNPs caused RL to be two days shorter, R1 two days later and R8 two days earlier, respectively. We identified association peaks acting independently over R1 and R8, suggesting that trait-specific minor effect loci are also involved in controlling R1 and R8. From the 111 genes highly associated with the three peaks detected in this study, we selected six candidate genes as the most likely cause of R1, R8, and RL variation. High correspondence was observed between a modifying variant SNP at position 04:39294836 in GmFulb and an association peak on Chr 4. Further studies using map-based cloning and fine mapping are necessary to elucidate the role of the candidates we identified for soybean maturity and adaptation to different latitudes and to be effectively used in the marker-assisted breeding of cultivars with optimal yield-related traits.

Cataloging SCN resistance loci in North American public soybean breeding programs.

Soybean cyst nematode (SCN) is a destructive pathogen of soybeans responsible for annual yield loss exceeding $1.5 billion in the United States. Here, we conducted a series of genome-wide association studies (GWASs) to understand the genetic landscape of SCN resistance in the University of Missouri soybean breeding programs (Missouri panel), as well as germplasm and cultivars within the United States Department of Agriculture (USDA) Uniform Soybean Tests-Northern Region (NUST). For the Missouri panel, we evaluated the resistance of breeding lines to SCN populations HG 2.5.7 (Race 1), HG 1.2.5.7 (Race 2), HG 0 (Race 3), HG 2.5.7 (Race 5), and HG 1.3.6.7 (Race 14) and identified seven quantitative trait nucleotides (QTNs) associated with SCN resistance on chromosomes 2, 8, 11, 14, 17, and 18. Additionally, we evaluated breeding lines in the NUST panel for resistance to SCN populations HG 2.5.7 (Race 1) and HG 0 (Race 3), and we found three SCN resistance-associated QTNs on chromosomes 7 and 18. Through these analyses, we were able to decipher the impact of seven major genetic loci, including three novel loci, on resistance to several SCN populations and identified candidate genes within each locus. Further, we identified favorable allelic combinations for resistance to individual SCN HG types and provided a list of available germplasm for integration of these unique alleles into soybean breeding programs. Overall, this study offers valuable insight into the landscape of SCN resistance loci in U.S. public soybean breeding programs and provides a framework to develop new and improved soybean cultivars with diverse plant genetic modes of SCN resistance.

Loss-of-function of an α-SNAP gene confers resistance to soybean cyst nematode.

Plant-parasitic nematodes are one of the most economically impactful pests in agriculture resulting in billions of dollars in realized annual losses worldwide. Soybean cyst nematode (SCN) is the number one biotic constraint on soybean production making it a priority for the discovery, validation and functional characterization of native plant resistance genes and genetic modes of action that can be deployed to improve soybean yield across the globe. Here, we present the discovery and functional characterization of a soybean resistance gene, GmSNAP02. We use unique bi-parental populations to fine-map the precise genomic location, and a combination of whole genome resequencing and gene fragment PCR amplifications to identify and confirm causal haplotypes. Lastly, we validate our candidate gene using CRISPR-Cas9 genome editing and observe a gain of resistance in edited plants. This demonstrates that the GmSNAP02 gene confers a unique mode of resistance to SCN through loss-of-function mutations that implicate GmSNAP02 as a nematode virulence target. We highlight the immediate impact of utilizing GmSNAP02 as a genome-editing-amenable target to diversify nematode resistance in commercially available cultivars.

Genomic Variations Explorer (GenVarX): a toolset for annotating promoter and CNV regions using genotypic and phenotypic differences.

The rapid growth of sequencing technology and its increasing popularity in biology-related research over the years has made whole genome re-sequencing (WGRS) data become widely available. A large amount of WGRS data can unlock the knowledge gap between genomics and phenomics through gaining an understanding of the genomic variations that can lead to phenotype changes. These genomic variations are usually comprised of allele and structural changes in DNA, and these changes can affect the regulatory mechanisms causing changes in gene expression and altering the phenotypes of organisms. In this research work, we created the GenVarX toolset, that is backed by transcription factor binding sequence data in promoter regions, the copy number variations data, SNPs and Indels data, and phenotypes data which can potentially provide insights about phenotypic differences and solve compelling questions in plant research. Analytics-wise, we have developed strategies to better utilize the WGRS data and mine the data using efficient data processing scripts, libraries, tools, and frameworks to create the interactive and visualization-enhanced GenVarX toolset that encompasses both promoter regions and copy number variation analysis components. The main capabilities of the GenVarX toolset are to provide easy-to-use interfaces for users to perform queries, visualize data, and interact with the data. Based on different input windows on the user interface, users can provide inputs corresponding to each field and submit the information as a query. The data returned on the results page is usually displayed in a tabular fashion. In addition, interactive figures are also included in the toolset to facilitate the visualization of statistical results or tool outputs. Currently, the GenVarX toolset supports soybean, rice, and Arabidopsis. The researchers can access the soybean GenVarX toolset from SoyKB via https://soykb.org/SoybeanGenVarX/, rice GenVarX toolset, and Arabidopsis GenVarX toolset from KBCommons web portal with links https://kbcommons.org/system/tools/GenVarX/Osativa and https://kbcommons.org/system/tools/GenVarX/Athaliana, respectively.

The Allele Catalog Tool: a web-based interactive tool for allele discovery and analysis.

BACKGROUND: The advancement of sequencing technologies today has made a plethora of whole-genome re-sequenced (WGRS) data publicly available. However, research utilizing the WGRS data without further configuration is nearly impossible. To solve this problem, our research group has developed an interactive Allele Catalog Tool to enable researchers to explore the coding region allelic variation present in over 1,000 re-sequenced accessions each for soybean, Arabidopsis, and maize. RESULTS: The Allele Catalog Tool was designed originally with soybean genomic data and resources. The Allele Catalog datasets were generated using our variant calling pipeline (SnakyVC) and the Allele Catalog pipeline (AlleleCatalog). The variant calling pipeline is developed to parallelly process raw sequencing reads to generate the Variant Call Format (VCF) files, and the Allele Catalog pipeline takes VCF files to perform imputations, functional effect predictions, and assemble alleles for each gene to generate curated Allele Catalog datasets. Both pipelines were utilized to generate the data panels (VCF files and Allele Catalog files) in which the accessions of the WGRS datasets were collected from various sources, currently representing over 1,000 diverse accessions for soybean, Arabidopsis, and maize individually. The main features of the Allele Catalog Tool include data query, visualization of results, categorical filtering, and download functions. Queries are performed from user input, and results are a tabular format of summary results by categorical description and genotype results of the alleles for each gene. The categorical information is specific to each species; additionally, available detailed meta-information is provided in modal popups. The genotypic information contains the variant positions, reference or alternate genotypes, the functional effect classes, and the amino-acid changes of each accession. Besides that, the results can also be downloaded for other research purposes. CONCLUSIONS: The Allele Catalog Tool is a web-based tool that currently supports three species: soybean, Arabidopsis, and maize. The Soybean Allele Catalog Tool is hosted on the SoyKB website ( https://soykb.org/SoybeanAlleleCatalogTool/ ), while the Allele Catalog Tool for Arabidopsis and maize is hosted on the KBCommons website ( https://kbcommons.org/system/tools/AlleleCatalogTool/Zmays and https://kbcommons.org/system/tools/AlleleCatalogTool/Athaliana ). Researchers can use this tool to connect variant alleles of genes with meta-information of species.

AccuCalc: A Python Package for Accuracy Calculation in GWAS.

The genome-wide association study (GWAS) is a popular genomic approach that identifies genomic regions associated with a phenotype and, thus, aims to discover causative mutations (CM) in the genes underlying the phenotype. However, GWAS discoveries are limited by many factors and typically identify associated genomic regions without the further ability to compare the viability of candidate genes and actual CMs. Therefore, the current methodology is limited to CM identification. In our recent work, we presented a novel approach to an empowered "GWAS to Genes" strategy that we named Synthetic phenotype to causative mutation (SP2CM). We established this strategy to identify CMs in soybean genes and developed a web-based tool for accuracy calculation (AccuTool) for a reference panel of soybean accessions. Here, we describe our further development of the tool that extends its utilization for other species and named it AccuCalc. We enhanced the tool for the analysis of datasets with a low-frequency distribution of a rare phenotype by automated formatting of a synthetic phenotype and added another accuracy-based GWAS evaluation criterion to the accuracy calculation. We designed AccuCalc as a Python package for GWAS data analysis for any user-defined species-independent variant calling format (vcf) or HapMap format (hmp) as input data. AccuCalc saves analysis outputs in user-friendly tab-delimited formats and also offers visualization of the GWAS results as Manhattan plots accentuated by accuracy. Under the hood of Python, AccuCalc is publicly available and, thus, can be used conveniently for the SP2CM strategy utilization for every species.

Natural and artificial selection of multiple alleles revealed through genomic analyses.

Genome-to-phenome research in agriculture aims to improve crops through in silico predictions. Genome-wide association study (GWAS) is potent in identifying genomic loci that underlie important traits. As a statistical method, increasing the sample quantity, data quality, or diversity of the GWAS dataset positively impacts GWAS power. For more precise breeding, concrete candidate genes with exact functional variants must be discovered. Many post-GWAS methods have been developed to narrow down the associated genomic regions and, ideally, to predict candidate genes and causative mutations (CMs). Historical natural selection and breeding-related artificial selection both act to change the frequencies of different alleles of genes that control phenotypes. With higher diversity and more extensive GWAS datasets, there is an increased chance of multiple alleles with independent CMs in a single causal gene. This can be caused by the presence of samples from geographically isolated regions that arose during natural or artificial selection. This simple fact is a complicating factor in GWAS-driven discoveries. Currently, none of the existing association methods address this issue and need to identify multiple alleles and, more specifically, the actual CMs. Therefore, we developed a tool that computes a score for a combination of variant positions in a single candidate gene and, based on the highest score, identifies the best number and combination of CMs. The tool is publicly available as a Python package on GitHub, and we further created a web-based Multiple Alleles discovery (MADis) tool that supports soybean and is hosted in SoyKB (https://soykb.org/SoybeanMADisTool/). We tested and validated the algorithm and presented the utilization of MADis in a pod pigmentation L1 gene case study with multiple CMs from natural or artificial selection. Finally, we identified a candidate gene for the pod color L2 locus and predicted the existence of multiple alleles that potentially cause loss of pod pigmentation. In this work, we show how a genomic analysis can be employed to explore the natural and artificial selection of multiple alleles and, thus, improve and accelerate crop breeding in agriculture.

A novel Synthetic phenotype association study approach reveals the landscape of association for genomic variants and phenotypes.

INTRODUCTION: Genome-Wide Association Studies (GWAS) identify tagging variants in the genome that are statistically associated with the phenotype because of their linkage disequilibrium (LD) relationship with the causative mutation (CM). When both low-density genotyped accession panels with phenotypes and resequenced data accession panels are available, tagging variants can assist with post-GWAS challenges in CM discovery. OBJECTIVES: Our objective was to identify additional GWAS evaluation criteria to assess correspondence between genomic variants and phenotypes, as well as enable deeper analysis of the localized landscape of association. METHODS: We used genomic variant positions as Synthetic phenotypes in GWAS that we named "Synthetic phenotype association study" (SPAS). The extreme case of SPAS is what we call an "Inverse GWAS" where we used CM positions of cloned soybean genes. We developed and validated the Accuracy concept as a measure of the correspondence between variant positions and phenotypes. RESULTS: The SPAS approach demonstrated that the genotype status of an associated variant used as a Synthetic phenotype enabled us to explore the relationships between tagging variants and CMs, and further, that utilizing CMs as Synthetic phenotypes in Inverse GWAS illuminated the landscape of association. We implemented the Accuracy calculation for a curated accession panel to an online Accuracy calculation tool (AccuTool) as a resource for gene identification in soybean. We demonstrated our concepts on three examples of soybean cloned genes. As a result of our findings, we devised an enhanced "GWAS to Genes" analysis (Synthetic phenotype to CM strategy, SP2CM). Using SP2CM, we identified a CM for a novel gene. CONCLUSION: The SP2CM strategy utilizing Synthetic phenotypes and the Accuracy calculation of correspondence provides crucial information to assist researchers in CM discovery. The impact of this work is a more effective evaluation of landscapes of GWAS associations.

Utilization of Plant Architecture Genes in Soybean to Positively Impact Adaptation to High Yield Environments.

Optimization of plant architecture by modifying stem termination and timing of flowering and maturity of soybean is a promising strategy to improve its adaptability to specific production environments. Therefore, it is important to choose a proper stem termination type and to understand morphological differences between each stem termination type under various environmental conditions. Variations in abruptness of stem termination have been generally classified into three classical genetic types, indeterminate (Dt1), determinate (dt1), and semi-determinate (Dt2). However, an additional stem termination type, termed tall determinate, and its genetic symbol, dt1-t, were introduced about 25 years ago. The tall determinate soybean lines show delayed cessation of apical stem growth and about 50% taller plant heights than the typical determinate soybeans, even though the genetic control of the tall determinate phenotype was found to be allelic to dt1. Despite the potential agronomic merits of the alternative stem termination type, knowledge about the tall determinate soybean remains limited. We clarified the molecular basis of the tall determinate stem termination type and examined potential agronomic merits of the alternative stem type under three different production environments in the US. Sequence analysis of the classical tall determinate soybean lines revealed that the dt1-t allele responsible for tall determinate stem architecture is caused by two of the identified independent missense alleles of dt1, dt1-t1 (R130K), and dt1-t2 (R62S). Also, from the comparison among soybean accessions belonging to each of the genotype categories for stem termination types, soybean accessions with tall determinate alleles were found to have a high discrepancy rate in phenotyping. Newly developed tall determinate late-maturing soybean germplasm lines had taller plant heights and a greater number of nodes with a similar stem diameter and similar pod density at the apical stem compared to typical determinate soybeans having dt1 (R166W) alleles in Southern environments in the US. The phenotype of increased pod-bearing nodes with lodging resistance has the potential to improve yield, especially grown in high yield environments. This study suggests an alternative strategy to remodel the shape of soybean plants, which can possibly lead to yield improvement through the modification of soybean plant architecture.

Candidate Genes Modulating Reproductive Timing in Elite US Soybean Lines Identified in Soybean Alleles of Arabidopsis Flowering Orthologs With Divergent Latitude Distribution.

Adaptation of soybean cultivars to the photoperiod in which they are grown is critical for optimizing plant yield. However, despite its importance, only the major loci conferring variation in flowering time and maturity of US soybean have been isolated. By contrast, over 200 genes contributing to floral induction in the model organism Arabidopsis thaliana have been described. In this work, putative alleles of a library of soybean orthologs of these Arabidopsis flowering genes were tested for their latitudinal distribution among elite US soybean lines developed in the United States. Furthermore, variants comprising the alleles of genes with significant differences in latitudinal distribution were assessed for amino acid conservation across disparate genera to infer their impact on gene function. From these efforts, several candidate genes from various biological pathways were identified that are likely being exploited toward adaptation of US soybean to various maturity groups.

Geographic distribution of the E1 family of genes and their effects on reproductive timing in soybean.

BACKGROUND: Soybean is an economically important crop which flowers predominantly in response to photoperiod. Several major loci controlling the quantitative trait for reproductive timing have been identified, of which allelic combinations at three of these loci, E1, E2, and E3, are the dominant factors driving time to flower and reproductive period. However, functional genomics studies have identified additional loci which affect reproductive timing, many of which are less understood. A better characterization of these genes will enable fine-tuning of adaptation to various production environments. Two such genes, E1La and E1Lb, have been implicated in flowering by previous studies, but their effects have yet to be assessed under natural photoperiod regimes. RESULTS: Natural and induced variants of E1La and E1Lb were identified and introgressed into lines harboring either E1 or its early flowering variant, e1-as. Lines were evaluated for days to flower and maturity in a Maturity Group (MG) III production environment. These results revealed that variation in E1La and E1Lb promoted earlier flowering and maturity, with stronger effects in e1-as background than in an E1 background. The geographic distribution of E1La alleles among wild and cultivated soybean revealed that natural variation in E1La likely contributed to northern expansion of wild soybean, while breeding programs in North America exploited e1-as to develop cultivars adapted to northern latitudes. CONCLUSION: This research identified novel alleles of the E1 paralogues, E1La and E1Lb, which promote flowering and maturity under natural photoperiods. These loci represent sources of genetic variation which have been under-utilized in North American breeding programs to control reproductive timing, and which can be valuable additions to a breeder's molecular toolbox.

Temporal changes in metabolism late in seed development affect biomass composition.

The negative association between protein and oil production in soybean (Glycine max) seed is well-documented. However, this inverse relationship is based primarily on the composition of mature seed, which reflects the cumulative result of events over the course of soybean seed development and therefore does not convey information specific to metabolic fluctuations during developmental growth regimes. In this study, we assessed maternal nutrient supply via measurement of seed coat exudates and metabolite levels within the cotyledon throughout development to identify trends in the accumulation of central carbon and nitrogen metabolic intermediates. Active metabolic activity during late seed development was probed through transient labeling with 13C substrates. The results indicated: (1) a drop in lipid contents during seed maturation with a concomitant increase in carbohydrates, (2) a transition from seed filling to maturation phases characterized by quantitatively balanced changes in carbon use and CO2 release, (3) changes in measured carbon and nitrogen resources supplied maternally throughout development, (4) 13C metabolite production through gluconeogenic steps for sustained carbohydrate accumulation as the maternal nutrient supply diminishes, and (5) oligosaccharide biosynthesis within the seed coat during the maturation phase. These results highlight temporal engineering targets for altering final biomass composition to increase the value of soybeans and a path to breaking the inverse correlation between seed protein and oil content.

Molecular-assisted breeding for soybean with high oleic/low linolenic acid and elevated vitamin E in the seed oil.

The uses of vegetable oils are determined by functional properties arising from their chemical composition. Soybean oil was previously used in margarines and baked foods after partial hydrogenation to achieve heat and oxidative stability. This process, however, generates trans fats that are now excluded from food use because of cardiovascular health risks. Also present in soybean oil are the anti-oxidant tocopherols, with α-tocopherol (vitamin E) typically present as a minor component compared to γ-tocopherol. Genetic improvement of the fatty acid profile and tocopherol profile is an attractive solution to increase the functional and health qualities of soybean oil. The objective of this research was to develop resources to directly select with molecular markers for the elevated vitamin E trait in soybean oil and to use a molecular breeding approach to combine elevated vitamin E with the high oleic/low linolenic acid seed oil trait that improves oil functionality and nutrition. New soybean germplasm was developed from the molecular breeding strategy that selected for alleles of six targeted genes. Seed oil from the novel soybean germplasm was confirmed to contain increased vitamin E α-tocopherol along with a high oleic acid/low linolenic acid profile.

The effects and interaction of soybean maturity gene alleles controlling flowering time, maturity, and adaptation in tropical environments.

BACKGROUND: Soybean is native to the temperate zones of East Asia. Poor yields of soybean in West African countries may be partially attributed to inadequate adaptation of soybean to tropical environments. Adaptation will require knowledge of the effects of allelic combinations of major maturity genes (E1, E2, and E3) and stem architecture. The long juvenile trait (J) influences soybean flowering time in short, ~ 12 h days, which characterize tropical latitudes. Soybean plant architecture includes determinate or indeterminate stem phenotypes controlled by the Dt1 gene. Understanding the influence of these genetic components on plant development and adaptation is key to optimize phenology and improve soybean yield potential in tropical environments. RESULTS: Soybean lines from five recombinant inbred populations were developed that varied in their combinations of targeted genes. The soybean lines were field tested in multiple environments and characterized for days to flowering (DTF), days to maturity (DTM), and plant height in locations throughout northern Ghana, and allelic combinations were determined for each line for associating genotype with phenotype. The results revealed significant differences based on genotype for DTF and DTM and allowed the comparison of different variant alleles of those genes. The mutant alleles of J and E1 had significant impact on DTF and DTM, and alleles of those genes interacted with each other for DTF but not DTM. The Dt1 gene significantly influenced plant height but not DTF or DTM. CONCLUSIONS: This research identified major and minor effect alleles of soybean genes that can be combined to control DTF, DTM, and plant height in short day tropical environments in Ghana. These phenotypes contribute to adaptation to a low latitude environment that can be optimized in a soybean breeding program with targeted selection of desired allele combinations. The knowledge of the genetic control of these traits will enhance molecular breeding to produce optimally adapted soybean varieties targeted to tropical environments.

Molecular-assisted breeding for improved carbohydrate profiles in soybean seed.

KEY MESSAGE: Two independent variant raffinose synthase 3 (RS3) alleles produced an equivalent phenotype and implicated the gene as a key contributor to soybean seed carbohydrate phenotype. Soybean is an important crop because the processed seed is utilized as a vegetable oil and a high protein meal typically used in livestock feeds. Raffinose and stachyose, the raffinose family of oligosaccharides (RFO) carbohydrate components of the seed, are synthesized in developing soybean seeds from sucrose and galactinol. Sucrose is considered positive for metabolizable energy, while RFO are anti-nutritional factors in diets of monogastric animals such as humans, poultry, and swine. To increase metabolizable energy available in soybean seed meal, prior research has been successful in deploying variant alleles of key soybean raffinose synthase (RS) genes leading to reductions or near elimination of seed RFO, with significant increases in seed sucrose. The objective of this research was to investigate the specific role of variants of the RS3 gene in a genomic context and improve molecular marker-assisted selection for the ultra-low (UL) RFO phenotype in soybean seeds. The results revealed a new variant of the RS3 allele (rs3 snp5, rs3 snp6) contributed to the UL RFO phenotype when mutant alleles of RS2 were present. The variant RS3 allele identified was present in about 15% of a small set of soybean cultivars released in North America. A missense allele of the RS3 gene (rs3 G75E) also produced the UL RFO phenotype when combined with mutant alleles of RS2. The discoveries reported here enable direct marker-assisted selection for an improved soybean meal trait that has the potential to add value to soybean by improving the metabolizable energy of the meal.

The development and use of a molecular model for soybean maturity groups.

BACKGROUND: Achieving appropriate maturity in a target environment is essential to maximizing crop yield potential. In soybean [Glycine max (L.) Merr.], the time to maturity is largely dependent on developmental response to dark periods. Once the critical photoperiod is reached, flowering is initiated and reproductive development proceeds. Therefore, soybean adaptation has been attributed to genetic changes and natural or artificial selection to optimize plant development in specific, narrow latitudinal ranges. In North America, these regions have been classified into twelve maturity groups (MG), with lower MG being shorter season than higher MG. Growing soybean lines not adapted to a particular environment typically results in poor growth and significant yield reductions. The objective of this study was to develop a molecular model for soybean maturity based on the alleles underlying the major maturity loci: E1, E2, and E3. RESULTS: We determined the allelic variation and diversity of the E maturity genes in a large collection of soybean landraces, North American ancestors, Chinese cultivars, North American cultivars or expired Plant Variety Protection lines, and private-company lines. The E gene status of accessions in the USDA Soybean Germplasm Collection with SoySNP50K Beadchip data was also predicted. We determined the E allelic combinations needed to adapt soybean to different MGs in the United States (US) and discovered a strong signal of selection for E genotypes released in North America, particularly the US and Canada. CONCLUSIONS: The E gene maturity model proposed will enable plant breeders to more effectively transfer traits into different MGs and increase the overall efficiency of soybean breeding in the US and Canada. The powerful yet simple selection strategy for increasing soybean breeding efficiency can be used alone or to directly enhance genomic prediction/selection schemes. The results also revealed previously unrecognized aspects of artificial selection in soybean imposed by soybean breeders based on geography that highlights the need for plant breeding that is optimized for specific environments.

Environmental Stability of Seed Carbohydrate Profiles in Soybeans Containing Different Alleles of the Raffinose Synthase 2 (RS2) Gene.

Soybean [Glycine max (L.) Merr.] is important for the high protein meal used for livestock feed formulations. Carbohydrates contribute positively or negatively to the potential metabolizable energy in soybean meal. The positive carbohydrate present in soybean meal consists primarily of sucrose, whereas the negative carbohydrate components are the raffinose family of oligosaccharides (RFOs), raffinose and stachyose. Increasing sucrose and decreasing raffinose and stachyose are critical targets to improve soybean. In three recently characterized lines, variant alleles of the soybean raffinose synthase 2 (RS2) gene were associated with increased sucrose and decreased RFOs. The objective of this research was to compare the environmental stability of seed carbohydrates in soybean lines containing wild-type or variant alleles of RS2 utilizing a field location study and a date of planting study. The results define the carbohydrate variation in distinct regional and temporal environments using soybean lines with different alleles of the RS2 gene.

Analysis of the siRNA-Mediated Gene Silencing Process Targeting Three Homologous Genes Controlling Soybean Seed Oil Quality.

In the past decade, RNA silencing has gained significant attention because of its success in genomic scale research and also in the genetic improvement of crop plants. However, little is known about the molecular basis of siRNA processing in association with its target transcript. To reveal this process for improving hpRNA-mediated gene silencing in crop plants, the soybean GmFAD3 gene family was chosen as a test model. We analyzed RNAi mutant soybean lines in which three members of the GmFAD3 gene family were silenced. The silencing levels of FAD3A, FAD3B and FAD3C were correlated with the degrees of sequence homology between the inverted repeat of hpRNA and the GmFAD3 transcripts in the RNAi lines. Strikingly, transgenes in two of the three RNAi lines were heavily methylated, leading to a dramatic reduction of hpRNA-derived siRNAs. Small RNAs corresponding to the loop portion of the hairpin transcript were detected while much lower levels of siRNAs were found outside of the target region. siRNAs generated from the 318-bp inverted repeat were found to be diced much more frequently at stem sequences close to the loop and associated with the inferred cleavage sites on the target transcripts, manifesting "hot spots". The top candidate hpRNA-derived siRNA share certain sequence features with mature miRNA. This is the first comprehensive and detailed study revealing the siRNA-mediated gene silencing mechanism in crop plants using gene family GmFAD3 as a test model.

Detection of novel QTLs for foxglove aphid resistance in soybean.

KEY MESSAGE: The Raso2 , novel QTL for Korea biotype foxglove aphid resistance in soybean from PI 366121 was identified on chromosome 7 using GoldenGate SNP microarray. Foxglove aphid, Aulacorthum solani (Kaltenbach), is a hemipteran insect that infects a wide variety of plants worldwide and causes serious yield losses in crops. The objective of this study was to identify the putative QTL for foxglove aphid resistance in wild soybean, PI 366121, (Glycine soja Sieb. and Zucc.). One hundred and forty-one F4-derived F8 recombinant inbred lines developed from a cross of susceptible Williams 82 and PI 366121 were used. The phenotyping of antibiosis and antixenosis resistance was done through choice and no-choice tests with total plant damage and primary infestation leaf damage; a genome-wide molecular linkage map was constructed with 504 single-nucleotide polymorphism markers utilizing a GoldenGate assay. Using inclusive composite interval mapping analysis for foxglove aphid resistance, one major candidate QTL on chromosome 7 and three minor QTL regions on chromosomes 3, 6 and 18 were identified. The major QTL on chromosome 7 showed both antixenosis and antibiosis resistance responses. However, the minor QTLs showed only antixenosis resistance response. The major QTL mapped to a different chromosome than the previously identified foxglove aphid resistance QTL, Raso1, from the cultivar Adams. Also, the responses to the Korea biotype foxglove aphid were different for Raso1, and the gene from PI 366121 against the Korea biotype foxglove aphid was different. Thus, the foxglove aphid resistance gene from PI 366121 was determined to be an independent gene from Raso1 and was designated as Raso2. This result could be useful in breeding for new foxglove aphid-resistant soybean cultivars.

Phytic Acid and Inorganic Phosphate Composition in Soybean Lines with Independent IPK1 Mutations.

Soybean seeds contain a large amount of P, which is stored as phytic acid (PA). Phytic acid is indigestible by nonruminant livestock and considered an antinutritional factor in soybean meal. Several low PA soybean lines have been discovered, but many of these lines have either minor reductions in PA or inadequate germination and emergence. The reduced PA phenotype of soybean line Gm-lpa-ZC-2 was previously shown to be the result of a mutation in a gene encoding an inositol pentakisphosphate 2-kinase on chromosome 14 (14IPK1). While the 14IPK1 mutation was shown to have no impact on germination and emergence, the reduction in PA was modest (up to 50%). Our objective was to determine the effect on seed P partitioning for a novel mutation of an independent IPK1 gene on chromosome six (06IPK1) on its own and in combination with mutant alleles of the 14IPK1. We developed soybean populations and conducted genotype and phenotype association analyses based on the genotype of the 06IPK1 and 14IPK1 genes and the seed P partitioning profile. The lines with both mutant IPK1 genes had very low PA levels, moderate accumulation of inorganic phosphate (Pi), and accumulation of high amounts of P in lower inositols. The developed lines did not have significant reductions in germination or field emergence. In addition, characterization of the lower inositols produced in the mutant lines suggests that IPK1 is a polyphosphate kinase and provides some insight into the PA biosynthesis pathway in soybean seeds.