Genetic architecture of protein and oil content in soybean seed and meal.

Soybean is grown primarily for the protein and oil extracted from its seed and its value is influenced by these components. The objective of this study was to map marker-trait associations (MTAs) for the concentration of seed protein, oil, and meal protein using the soybean nested association mapping (SoyNAM) population. The composition traits were evaluated on seed harvested from over 5000 inbred lines of the SoyNAM population grown in 10 field locations across 3 years. Estimated heritabilities were at least 0.85 for all three traits. The genotyping of lines with single nucleotide polymorphism markers resulted in the identification of 107 MTAs for the three traits. When MTAs for the three traits that mapped within 5 cM intervals were binned together, the MTAs were mapped to 64 intervals on 19 of the 20 soybean chromosomes. The majority of the MTA effects were small and of the 107 MTAs, 37 were for protein content, 39 for meal protein, and 31 for oil content. For cases where a protein and oil MTAs mapped to the same interval, most (94%) significant effects were opposite for the two traits, consistent with the negative correlation between these traits. A coexpression analysis identified candidate genes linked to MTAs and 18 candidate genes were identified. The large number of small effect MTAs for the composition traits suggest that genomic prediction would be more effective in improving these traits than marker-assisted selection.

Genomic selection performs as effectively as phenotypic selection for increasing seed yield in soybean.

Increasing the rate of genetic gain for seed yield remains the primary breeding objective in both public and private soybean [Glycine max (L.) Merr.] breeding programs. Genomic selection (GS) has the potential to accelerate the rate of genetic gain for soybean seed yield. Limited studies to date have validated GS accuracy and directly compared GS with phenotypic selection (PS), and none have been reported in soybean. This study conducted the first empirical validation of GS for increasing seed yield using over 1,500 lines and over 7 yr (2010-2016) of replicated experiments in the University of Nebraska-Lincoln soybean breeding program. The study was designed to capture the varying genetic relatedness of the training population to three validation sets: two large biparental populations (TBP-1 and TBP-2) and a large validation set comprised of 457 preselected advanced lines derived from 45 biparental populations (TMP). We found that prediction accuracy (.54) realized in our validation experiments was comparable with what we obtained from a series of cross-validation experiments (.64). Both GS and PS were more effective for increasing the population mean performance compared with random selection (RS). We found a selection advantage of GS over PS, where higher genetic gain and identification of top-performing lines was maximized at 10-20% selected proportion. Genomic selection led to small increases in genetic similarity when compared with PS and RS presumably because of a significant shift on allelic frequencies toward the extremes, suggesting that it could erode genetic diversity more quickly. Overall, we found that GS can perform as effectively as PS but that measures should be considered to protect against loss of genetic variance when using GS.

Decoupling between leaf nitrogen and radiation use efficiency in vegetative and early reproductive stages in high-yielding soybean.

Ontogenic changes in soybean radiation use efficiency (RUE) have been attributed to variation in specific leaf nitrogen (SLN) based only on data collected during seed filling. We evaluated this hypothesis using data on leaf area, absorbed radiation (ARAD), aboveground dry matter (ADM), and plant nitrogen (N) concentration collected during the entire crop season from seven field experiments conducted in a stress-free environment. Each experiment included a full-N treatment that received ample N fertilizer and a zero-N treatment that relied on N fixation and soil N mineralization. We estimated RUE based on changes in ADM between sampling times and associated ARAD, accounting for changes in biomass composition. The RUE and SLN exhibited different seasonal patterns: a bell-shaped pattern with a peak around the beginning of seed filling, and a convex pattern followed by an abrupt decline during late seed filling, respectively. Changes in SLN explained the decline in RUE during seed filling but failed to predict changes in RUE in earlier stages and underestimated the maximum RUE observed during pod setting. Comparison between observed and simulated RUE using a process-based crop simulation model revealed similar discrepancies. The decoupling between RUE and SLN during early crop stages suggests that leaf N is above that needed to maximize crop growth but may play a role in storing N that can be used in later reproductive stages to meet the large seed N demand associated with high-yielding crops.

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.

An Assessment of the Factors Influencing the Prediction Accuracy of Genomic Prediction Models Across Multiple Environments.

The effects of climate change create formidable challenges for breeders striving to produce sufficient food quantities in rapidly changing environments. It is therefore critical to investigate the ability of multi-environment genomic prediction (GP) models to predict genomic estimated breeding values (GEBVs) in extreme environments. Exploration of the impact of training set composition on the accuracy of such GEBVs is also essential. Accordingly, we examined the influence of the number of training environments and the use of environmental covariates (ECs) in GS models on four subsets of n = 500 lines of the soybean nested association mapping (SoyNAM) panel grown in nine environments in the US-North Central Region. The ensuing analyses provided insights into the influence of both of these factors for predicting grain yield in the most and the least extreme of these environments. We found that only a subset of the available environments was needed to obtain the highest observed prediction accuracies. The inclusion of ECs in the GP model did not substantially increase prediction accuracies relative to competing models, and instead more often resulted in negative prediction accuracies. Combined with the overall low prediction accuracies for grain yield in the most extreme environment, our findings highlight weaknesses in current GP approaches for prediction in extreme environments, and point to specific areas on which to focus future research efforts.

Comparing a Mixed Model Approach to Traditional Stability Estimators for Mapping Genotype by Environment Interactions and Yield Stability in Soybean [Glycine max (L.) Merr.].

Identifying genetic loci associated with yield stability has helped plant breeders and geneticists begin to understand the role and influence of genotype by environment (GxE) interactions in soybean [Glycine max (L.) Merr.] productivity, as well as other crops. Quantifying a genotype's range of performance across testing locations has been developed over decades with dozens of methodologies available. This includes directly modeling GxE interactions as part of an overall model for yield, as well as methods which generate overall yield "stability" values from multi-environment trial data. Correspondence between these methods as it pertains to the outcomes of genome wide association studies (GWAS) has not been well defined. In this study, the GWAS results for yield and yield stability were compared in 213 soybean lines across 11 environments to determine their utility and potential intersection. Both univariate and multivariate conventional stability estimates were considered alongside a mixed model for yield that fit marker by environment interactions as a random effect. One-hundred and six total QTL were discovered across all mapping results, however, genetic loci that were significant in the mixed model for grain yield that fit marker by environment interactions were completely distinct from those that were significant when mapping using traditional stability measures as a phenotype. Furthermore, 73.21% of QTL discovered in the mixed model were determined to cause a crossover interaction effect which cause genotype rank changes between environments. Overall, the QTL discovered via explicitly mapping GxE interactions also explained more yield variance that those QTL associated with differences in traditional stability estimates making their theoretical impact on selection greater. A lack of intersecting results between mapping approaches highlights the importance of examining stability in multiple contexts when attempting to manipulate GxE interactions in soybean.

Early Prediction of Soybean Traits through Color and Texture Features of Canopy RGB Imagery.

Global crop production is facing the challenge of a high projected demand, while the yields of major crops are not increasing at sufficient speeds. Crop breeding is an important way to boost crop productivity, however its improvement rate is partially hindered by the long crop generation cycles. If end-season crop traits such as yield can be predicted through early-season phenotypic measurements, crop selection can potentially be made before a full crop generation cycle finishes. This study explored the possibility of predicting soybean end-season traits through the color and texture features of early-season canopy images. Six thousand three hundred and eighty-three images were captured at V4/V5 growth stage over 6039 soybean plots growing at four locations. One hundred and forty color features and 315 gray-level co-occurrence matrix-based texture features were derived from each image. Another two variables were also introduced to account for location and timing differences between the images. Five regression and five classification techniques were explored. Best results were obtained using all 457 predictor variables, with Cubist as the regression technique and Random Forests as the classification technique. Yield (RMSE = 9.82, R2 = 0.68), Maturity (RMSE = 3.70, R2 = 0.76) and Seed Size (RMSE = 1.63, R2 = 0.53) were identified as potential soybean traits that might be early predictable.

Generating High Density, Low Cost Genotype Data in Soybean [Glycine max (L.) Merr.].

Obtaining genome-wide genotype information for millions of SNPs in soybean [Glycine max (L.) Merr.] often involves completely resequencing a line at 5X or greater coverage. Currently, hundreds of soybean lines have been resequenced at high depth levels with their data deposited in the NCBI Short Read Archive. This publicly available dataset may be leveraged as an imputation reference panel in combination with skim (low coverage) sequencing of new soybean genotypes to economically obtain high-density SNP information. Ninety-nine soybean lines resequenced at an average of 17.1X were used to generate a reference panel, with over 10 million SNPs called using GATK's Haplotype Caller tool. Whole genome resequencing at approximately 1X depth was performed on 114 previously ungenotyped experimental soybean lines. Coverages down to 0.1X were analyzed by randomly subsetting raw reads from the original 1X sequence data. SNPs discovered in the reference panel were genotyped in the experimental lines after aligning to the soybean reference genome, and missing markers imputed using Beagle 4.1. Sequencing depth of the experimental lines could be reduced to 0.3X while still retaining an accuracy of 97.8%. Accuracy was inversely related to minor allele frequency, and highly correlated with marker linkage disequilibrium. The high accuracy of skim sequencing combined with imputation provides a low cost method for obtaining dense genotypic information that can be used for various genomics applications in soybean.

Response Surface Analysis of Genomic Prediction Accuracy Values Using Quality Control Covariates in Soybean.

An important and broadly used tool for selection purposes and to increase yield and genetic gain in plant breeding programs is genomic prediction (GP). Genomic prediction is a technique where molecular marker information and phenotypic data are used to predict the phenotype (eg, yield) of individuals for which only marker data are available. Higher prediction accuracy can be achieved not only by using efficient models but also by using quality molecular marker and phenotypic data. The steps of a typical quality control (QC) of marker data include the elimination of markers with certain level of minor allele frequency (MAF) and missing marker values and the imputation of missing marker values. In this article, we evaluated how the prediction accuracy is influenced by the combination of 12 MAF values, 27 different percentages of missing marker values, and 2 imputation techniques (IT; naïve and Random Forest (RF)). We constructed a response surface of prediction accuracy values for the two ITs as a function of MAF and percentage of missing marker values using soybean data from the University of Nebraska-Lincoln Soybean Breeding Program. We found that both the genetic architecture of the trait and the IT affect the prediction accuracy implying that we have to be careful how we perform QC on the marker data. For the corresponding combinations MAF-percentage of missing values we observed that implementing the RF imputation increased the number of markers by 2 to 5 times than the simple naïve imputation method that is based on the mean allele dosage of the non-missing values at each loci. We conclude that there is not a unique strategy (combination of the QCs and imputation method) that outperforms the results of the others for all traits.

Genetic Architecture of Soybean Yield and Agronomic Traits.

Soybean is the world's leading source of vegetable protein and demand for its seed continues to grow. Breeders have successfully increased soybean yield, but the genetic architecture of yield and key agronomic traits is poorly understood. We developed a 40-mating soybean nested association mapping (NAM) population of 5,600 inbred lines that were characterized by single nucleotide polymorphism (SNP) markers and six agronomic traits in field trials in 22 environments. Analysis of the yield, agronomic, and SNP data revealed 23 significant marker-trait associations for yield, 19 for maturity, 15 for plant height, 17 for plant lodging, and 29 for seed mass. A higher frequency of estimated positive yield alleles was evident from elite founder parents than from exotic founders, although unique desirable alleles from the exotic group were identified, demonstrating the value of expanding the genetic base of US soybean breeding.

Field-Based Scoring of Soybean Iron Deficiency Chlorosis Using RGB Imaging and Statistical Learning.

Iron deficiency chlorosis (IDC) is an abiotic stress in soybean that can cause significant biomass and yield reduction. IDC is characterized by stunted growth and yellowing and interveinal chlorosis of early trifoliate leaves. Scoring IDC severity in the field is conventionally done by visual assessment. The goal of this study was to investigate the usefulness of Red Green Blue (RGB) images of soybean plots captured under the field condition for IDC scoring. A total of 64 soybean lines with four replicates were planted in 6 fields over 2 years. Visual scoring (referred to as Field Score, or FS) was conducted at V3-V4 growth stage; and concurrently RGB images of the field plots were recorded with a high-throughput field phenotyping platform. A second set of IDC scores was done on the plot images (displayed on a computer screen) consistently by one person in the office (referred to as Office Score, or OS). Plot images were then processed to remove weeds and extract six color features, which were used to train computer-based IDC scoring models (referred to as Computer Score, or CS) using linear discriminant analysis (LDA) and support vector machine (SVM). The results showed that, in the fields where severe IDC symptoms were present, FS and OS were strongly positively correlated with each other, and both of them were strongly negatively correlated with yield. CS could satisfactorily predict IDC scores when evaluated using FS and OS as the reference (overall classification accuracy > 81%). SVM models appeared to outperform LDA models; and the SVM model trained to predict IDC OS gave the highest prediction accuracy. It was anticipated that coupling RGB imaging from the high-throughput field phenotyping platform with real-time image processing and IDC CS models would lead to a more rapid, cost-effective, and objective scoring pipeline for soybean IDC field screening and breeding.

An epigenetic breeding system in soybean for increased yield and stability.

Epigenetic variation has been associated with a wide range of adaptive phenotypes in plants, but there exist few direct means for exploiting this variation. RNAi suppression of the plant-specific gene, MutS HOMOLOG1 (MSH1), in multiple plant species produces a range of developmental changes accompanied by modulation of defence, phytohormone and abiotic stress response pathways along with methylome repatterning. This msh1-conditioned developmental reprogramming is retained independent of transgene segregation, giving rise to transgene-null 'memory' effects. An isogenic memory line crossed to wild type produces progeny families displaying increased variation in adaptive traits that respond to selection. This study investigates amenability of the MSH1 system for inducing agronomically valuable epigenetic variation in soybean. We developed MSH1 epi-populations by crossing with msh1-acquired soybean memory lines. Derived soybean epi-lines showed increase in variance for multiple yield-related traits including pods per plant, seed weight and maturity time in both glasshouse and field trials. Selected epi-F2:4 and epi-F2:5 lines showed an increase in seed yield over wild type. By epi-F2:6, we observed a return of MSH1-derived enhanced growth back to wild-type levels. Epi-populations also showed evidence of reduced epitype-by-environment (e × E) interaction, indicating higher yield stability. Transcript profiling of epi-lines identified putative signatures of enhanced growth behaviour across generations. Genes related to cell cycle, abscisic acid biosynthesis and auxin response, particularly SMALL AUXIN UP RNAs (SAURs), were differentially expressed in epi-F2:4 lines that showed increased yield when compared to epi-F2:6 . These data support the potential of MSH1-derived epigenetic variation in plant breeding for enhanced yield and yield stability.

Gene Expression Profiling of Iron Deficiency Chlorosis Sensitive and Tolerant Soybean Indicates Key Roles for Phenylpropanoids under Alkalinity Stress.

Alkaline soils comprise 30% of the earth and have low plant-available iron (Fe) concentration, and can cause iron deficiency chlorosis (IDC). IDC causes soybean yield losses of $260 million annually. However, it is not known whether molecular responses to IDC are equivalent to responses to low iron supply. IDC tolerant and sensitive soybean lines provide a contrast to identify specific factors associated with IDC. We used RNA-seq to compare gene expression under combinations of normal pH (5.7) or alkaline pH (7.7, imposed by 2.5 mM bicarbonate, or pH 8.2 imposed by 5 mM bicarbonate) and normal (25 µM) or low (1 µM) iron conditions from roots of these lines. Thus, we were able to treat pH and Fe supply as separate variables. We also noted differential gene expression between IDC sensitive and tolerant genotypes in each condition. Classical iron uptake genes, including ferric-chelate reductase (FCR) and ferrous transporters, were upregulated by both Fe deficiency and alkaline stress, however, their gene products did not function well at alkaline pH. In addition, genes in the phenylpropanoid synthesis pathway were upregulated in both alkaline and low Fe conditions. These genes lead to the production of fluorescent root exudate (FluRE) compounds, such as coumarins. Fluorescence of nutrient solution increased with alkaline treatment, and was higher in the IDC tolerant line. Some of these genes also localized to previously identified QTL regions associated with IDC. We hypothesize that FluRE become essential at alkaline pH where the classical iron uptake system does not function well. This work could result in new strategies to screen for IDC tolerance, and provide breeding targets to improve crop alkaline stress tolerance.

Genome-Wide Analysis of Grain Yield Stability and Environmental Interactions in a Multiparental Soybean Population.

Genetic improvement toward optimized and stable agronomic performance of soybean genotypes is desirable for food security. Understanding how genotypes perform in different environmental conditions helps breeders develop sustainable cultivars adapted to target regions. Complex traits of importance are known to be controlled by a large number of genomic regions with small effects whose magnitude and direction are modulated by environmental factors. Knowledge of the constraints and undesirable effects resulting from genotype by environmental interactions is a key objective in improving selection procedures in soybean breeding programs. In this study, the genetic basis of soybean grain yield responsiveness to environmental factors was examined in a large soybean nested association population. For this, a genome-wide association to performance stability estimates generated from a Finlay-Wilkinson analysis and the inclusion of the interaction between marker genotypes and environmental factors was implemented. Genomic footprints were investigated by analysis and meta-analysis using a recently published multiparent model. Results indicated that specific soybean genomic regions were associated with stability, and that multiplicative interactions were present between environments and genetic background. Seven genomic regions in six chromosomes were identified as being associated with genotype-by-environment interactions. This study provides insight into genomic assisted breeding aimed at achieving a more stable agronomic performance of soybean, and documented opportunities to exploit genomic regions that were specifically associated with interactions involving environments and subpopulations.

Dissecting the Genetic Basis of Local Adaptation in Soybean.

Soybean (Glycine max) is the most widely grown oilseed in the world and is an important source of protein for both humans and livestock. Soybean is widely adapted to both temperate and tropical regions, but a changing climate demands a better understanding of adaptation to specific environmental conditions. Here, we explore genetic variation in a collection of 3,012 georeferenced, locally adapted landraces from a broad geographical range to help elucidate the genetic basis of local adaptation. We used geographic origin, environmental data and dense genome-wide SNP data to perform an environmental association analysis and discover loci displaying steep gradients in allele frequency across geographical distance and between landrace and modern cultivars. Our combined application of methods in environmental association mapping and detection of selection targets provide a better understanding of how geography and selection may have shaped genetic variation among soybean landraces. Moreover, we identified several important candidate genes related to drought and heat stress, and revealed important genomic regions possibly involved in the geographic divergence of soybean.

Genome-wide Association Mapping of Qualitatively Inherited Traits in a Germplasm Collection.

Genome-wide association (GWA) has been used as a tool for dissecting the genetic architecture of quantitatively inherited traits. We demonstrate here that GWA can also be highly useful for detecting many major genes governing categorically defined phenotype variants that exist for qualitatively inherited traits in a germplasm collection. Genome-wide association mapping was applied to categorical phenotypic data available for 10 descriptive traits in a collection of ∼13,000 soybean [ (L.) Merr.] accessions that had been genotyped with a 50,000 single nucleotide polymorphism (SNP) chip. A GWA on a panel of accessions of this magnitude can offer substantial statistical power and mapping resolution, and we found that GWA mapping resulted in the identification of strong SNP signals for 24 classical genes as well as several heretofore unknown genes controlling the phenotypic variants in those traits. Because some of these genes had been cloned, we were able to show that the narrow GWA mapping SNP signal regions that we detected for the phenotypic variants had chromosomal bp spans that, with just one exception, overlapped the bp region of the cloned genes, despite local variation in SNP number and nonuniform SNP distribution in the chip set.

Multi-Population Selective Genotyping to Identify Soybean [Glycine max (L.) Merr.] Seed Protein and Oil QTLs.

Plant breeders continually generate ever-higher yielding cultivars, but also want to improve seed constituent value, which is mainly protein and oil, in soybean [Glycine max (L.) Merr.]. Identification of genetic loci governing those two traits would facilitate that effort. Though genome-wide association offers one such approach, selective genotyping of multiple biparental populations offers a complementary alternative, and was evaluated here, using 48 F2:3 populations (n = âˆ¼224 plants) created by mating 48 high protein germplasm accessions to cultivars of similar maturity, but with normal seed protein content. All F2:3 progeny were phenotyped for seed protein and oil, but only 22 high and 22 low extreme progeny in each F2:3 phenotypic distribution were genotyped with a 1536-SNP chip (ca 450 bimorphic SNPs detected per mating). A significant quantitative trait locus (QTL) on one or more chromosomes was detected for protein in 35 (73%), and for oil in 25 (52%), of the 48 matings, and these QTL exhibited additive effects of ≥ 4 g kg(-1) and R(2) values of 0.07 or more. These results demonstrated that a multiple-population selective genotyping strategy, when focused on matings between parental phenotype extremes, can be used successfully to identify germplasm accessions possessing large-effect QTL alleles. Such accessions would be of interest to breeders to serve as parental donors of those alleles in cultivar development programs, though 17 of the 48 accessions were not unique in terms of SNP genotype, indicating that diversity among high protein accessions in the germplasm collection is less than what might ordinarily be assumed.

Soy protein is beneficial but high-fat diet and voluntary running are detrimental to bone structure in mice.

Physical activity and soy protein isolate (SPI) augmentation have been reported to be beneficial for bone health. We hypothesized that combining voluntary running and SPI intake would alleviate detrimental changes in bone induced by a high-fat diet. A 2 × 2 × 2 experiment was designed with diets containing 16% or 45% of energy as corn oil and 20% SPI or casein fed to sedentary or running male C57BL/6 mice for 14 weeks. Distal femurs were assessed for microstructural changes. The high-fat diet significantly decreased trabecular number (Tb.N) and bone mineral density (BMD) and increased trabecular separation (Tb.Sp). Soy protein instead of casein, regardless of fat content, in the diet significantly increased bone volume fraction, Tb.N, connectivity density, and BMD and decreased Tb.Sp. Voluntary running, regardless of fat content, significantly decreased bone volume fraction, Tb.N, connectivity density, and BMD and increased Tb.Sp. The high-fat diet significantly decreased osteocalcin and increased tartrate-resistant acid phosphatase 5b (TRAP 5b) concentrations in plasma. Plasma concentrations of osteocalcin were increased by both SPI and running. Running alleviated the increase in TRAP 5b induced by the high-fat diet. These findings demonstrate that a high-fat diet is deleterious, and SPI is beneficial to trabecular bone properties. The deleterious effect of voluntary running on trabecular structural characteristics indicates that there may be a maximal threshold of running beyond which beneficial effects cease and detrimental effects occur. Increases in plasma osteocalcin and decreases in plasma TRAP 5b in running mice suggest that a compensatory response occurs to counteract the detrimental effects of excessive running.

Genotyping by sequencing for genomic prediction in a soybean breeding population.

BACKGROUND: Advances in genotyping technology, such as genotyping by sequencing (GBS), are making genomic prediction more attractive to reduce breeding cycle times and costs associated with phenotyping. Genomic prediction and selection has been studied in several crop species, but no reports exist in soybean. The objectives of this study were (i) evaluate prospects for genomic selection using GBS in a typical soybean breeding program and (ii) evaluate the effect of GBS marker selection and imputation on genomic prediction accuracy. To achieve these objectives, a set of soybean lines sampled from the University of Nebraska Soybean Breeding Program were genotyped using GBS and evaluated for yield and other agronomic traits at multiple Nebraska locations. RESULTS: Genotyping by sequencing scored 16,502 single nucleotide polymorphisms (SNPs) with minor-allele frequency (MAF) > 0.05 and percentage of missing values ≤ 5% on 301 elite soybean breeding lines. When SNPs with up to 80% missing values were included, 52,349 SNPs were scored. Prediction accuracy for grain yield, assessed using cross validation, was estimated to be 0.64, indicating good potential for using genomic selection for grain yield in soybean. Filtering SNPs based on missing data percentage had little to no effect on prediction accuracy, especially when random forest imputation was used to impute missing values. The highest accuracies were observed when random forest imputation was used on all SNPs, but differences were not significant. A standard additive G-BLUP model was robust; modeling additive-by-additive epistasis did not provide any improvement in prediction accuracy. The effect of training population size on accuracy began to plateau around 100, but accuracy steadily climbed until the largest possible size was used in this analysis. Including only SNPs with MAF > 0.30 provided higher accuracies when training populations were smaller. CONCLUSIONS: Using GBS for genomic prediction in soybean holds good potential to expedite genetic gain. Our results suggest that standard additive G-BLUP models can be used on unfiltered, imputed GBS data without loss in accuracy.

Stacking of a stearoyl-ACP thioesterase with a dual-silenced palmitoyl-ACP thioesterase and ∆12 fatty acid desaturase in transgenic soybean.

Soybean (Glycine max (L.) Merr) is valued for both its protein and oil, whose seed is composed of 40% and 20% of each component, respectively. Given its high percentage of polyunsaturated fatty acids, linoleic acid and linolenic acid, soybean oil oxidative stability is relatively poor. Historically food processors have employed a partial hydrogenation process to soybean oil as a means to improve both the oxidative stability and functionality in end-use applications. However, the hydrogenation process leads to the formation of trans-fats, which are associated with negative cardiovascular health. As a means to circumvent the need for the hydrogenation process, genetic approaches are being pursued to improve oil quality in oilseeds. In this regard, we report here on the introduction of the mangosteen (Garcinia mangostana) stearoyl-ACP thioesterase into soybean and the subsequent stacking with an event that is dual-silenced in palmitoyl-ACP thioesterase and ∆12 fatty acid desaturase expression in a seed-specific fashion. Phenotypic analyses on transgenic soybean expressing the mangosteen stearoyl-ACP thioesterase revealed increases in seed stearic acid levels up to 17%. The subsequent stacked with a soybean event silenced in both palmitoyl-ACP thioesterase and ∆12 fatty acid desaturase activity, resulted in a seed lipid phenotype of approximately 11%-19% stearate and approximately 70% oleate. The oil profile created by the stack was maintained for four generations under greenhouse conditions and a fifth generation under a field environment. However, in generation six and seven under field conditions, the oleate levels decreased to 30%-40%, while the stearic level remained elevated.