A Multi-Year, Multi-Cultivar Approach to Differential Expression Analysis of High- and Low-Protein Soybean (Glycine max).

Soybean (Glycine max (L.) Merr.) is among the most valuable crops based on its nutritious seed protein and oil. Protein quality, evaluated as the ratio of glycinin (11S) to ß-conglycinin (7S), can play a role in food and feed quality. To help uncover the underlying differences between high and low protein soybean varieties, we performed differential expression analysis on high and low total protein soybean varieties and high and low 11S soybean varieties grown in four locations across Eastern and Western Canada over three years (2018-2020). Simultaneously, ten individual differential expression datasets for high vs. low total protein soybeans and ten individual differential expression datasets for high vs. low 11S soybeans were assessed, for a total of 20 datasets. The top 15 most upregulated and the 15 most downregulated genes were extracted from each differential expression dataset and cross-examination was conducted to create shortlists of the most consistently differentially expressed genes. Shortlisted genes were assessed for gene ontology to gain a global appreciation of the commonly differentially expressed genes. Genes with roles in the lipid metabolic pathway and carbohydrate metabolic pathway were differentially expressed in high total protein and high 11S soybeans in comparison to their low total protein and low 11S counterparts. Expression differences were consistent between East and West locations with the exception of one, Glyma.03G054100. These data are important for uncovering the genes and biological pathways responsible for the difference in seed protein between high and low total protein or 11S cultivars.

Human-Soybean Allergies: Elucidation of the Seed Proteome and Comprehensive Protein-Protein Interaction Prediction.

The soybean crop, Glycine max (L.) Merr., is consumed by humans, Homo sapiens, worldwide. While the respective bodies of literature and -omics data for each of these organisms are extensive, comparatively few studies investigate the molecular biological processes occurring between the two. We are interested in elucidating the network of protein-protein interactions (PPIs) involved in human-soybean allergies. To this end, we leverage state-of-the-art sequence-based PPI predictors amenable to predicting the enormous comprehensive interactome between human and soybean. A network-based analytical approach is proposed, leveraging similar interaction profiles to identify candidate allergens and proteins involved in the allergy response. Interestingly, the predicted interactome can be explored from two complementary perspectives: which soybean proteins are predicted to interact with specific human proteins and which human proteins are predicted to interact with specific soybean proteins. A total of eight proteins (six specific to the human proteome and two to the soy proteome) have been identified and supported by the literature to be involved in human health, specifically related to immunological and neurological pathways. This study, beyond generating the most comprehensive human-soybean interactome to date, elucidated a soybean seed interactome and identified several proteins putatively consequential to human health.

Soybean (Glycine max) Haplotype Map (GmHapMap): a universal resource for soybean translational and functional genomics.

Here, we describe a worldwide haplotype map for soybean (GmHapMap) constructed using whole-genome sequence data for 1007 Glycine max accessions and yielding 14.9 million variants as well as 4.3 M tag single-nucleotide polymorphisms (SNPs). When sampling random subsets of these accessions, the number of variants and tag SNPs plateaued beyond approximately 800 and 600 accessions, respectively. This suggests extensive coverage of diversity within the cultivated soybean. GmHapMap variants were imputed onto 21 618 previously genotyped accessions with up to 96% success for common alleles. A local association analysis was performed with the imputed data using markers located in a 1-Mb region known to contribute to seed oil content and enabled us to identify a candidate causal SNP residing in the NPC1 gene. We determined gene-centric haplotypes (407 867 GCHs) for the 55 589 genes and showed that such haplotypes can help to identify alleles that differ in the resulting phenotype. Finally, we predicted 18 031 putative loss-of-function (LOF) mutations in 10 662 genes and illustrated how such a resource can be used to explore gene function. The GmHapMap provides a unique worldwide resource for applied soybean genomics and breeding.

Identifying new variation at the J locus, previously identified as e6, in long juvenile 'Paranagoiana' soybean.

KEY MESSAGE: A previously identified soybean maturity locus, E6, is discovered to be J, with the long juvenile allele in Paranagoiana now deemed j-x. Soybean grown at latitudes of ~20° or lower can produce lower grain yields due to the short days. This limitation can be overcome by using the long juvenile trait (LJ) which delays flowering under short day conditions. Two LJ loci have been mapped to the same location on Gm04, J and E6. The objective of this research was to investigate the e6 allele in 'Paranagoiana' and determine if E6 and J are the same locus or linked loci. KASP markers showed that e6 lines did not have the j-1 allele of LJ PI 159925. A population fixed for E1 but segregating for E6, with e6 introgressed from Paranagoiana, showed single gene control for flowering and maturity under short days. Sequencing Glyma.04G050200, the J gene, with long amplification Taq found that the e6 line 'Paranagoiana' contains a Ty1-copia retrotransposon of ~10,000 bp, inserted within exon 4. PCR amplification of the cDNA of Glyma.04G050200 also showed differences between the mRNA sequences (presence of insertion in j-x). Hence, we conclude that the loci E6 and J are one locus and deem this new variation found in Paranagoiana as j-x.

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.

Effect of change from a determinate to a semi-determinate growth habit on the yield and lodging resistance of soybeans in the northeast region of Japan.

Although an indeterminate growth habit is attractive to develop high-yield soybean varieties with higher number of pods (Glycine max (L). Merr.), lodging in indeterminate varieties remains a problem in Japan. As the semi-determinate varieties have shorter main stem length than the indeterminate varieties, this trait can be useful to improve varieties with high yield and low lodging risk. We introduced the genes Dt1 and Dt2, which regulate stem growth habit, into three determinate varieties by backcrossing and evaluated the resulting effects on yield and lodging tendency under four different growing environments. The yield and lodging degree of the semi-determinate and indeterminate lines were higher and more severe than those of the determinate lines. Despite the lower overall lodging score, the semi-determinate lines had marginally lower overall yield than that of the indeterminate lines. However, the effect of introduction of semi-determinate traits on yield and lodging degree was different in the three backgrounds, with the yield of semi-determinate lines being the highest and the difference in lodging degree between the semi-determinate and determinate lines being under 1.0 in one background. Therefore, semi-determinate growth habit has potential to develop high yielding varieties with low lodging risk.

Comprehensive description of genomewide nucleotide and structural variation in short-season soya bean.

Next-generation sequencing (NGS) and bioinformatics tools have greatly facilitated the characterization of nucleotide variation; nonetheless, an exhaustive description of both SNP haplotype diversity and of structural variation remains elusive in most species. In this study, we sequenced a representative set of 102 short-season soya beans and achieved an extensive coverage of both nucleotide diversity and structural variation (SV). We called close to 5M sequence variants (SNPs, MNPs and indels) and noticed that the number of unique haplotypes had plateaued within this set of germplasm (1.7M tag SNPs). This data set proved highly accurate (98.6%) based on a comparison of called genotypes at loci shared with a SNP array. We used this catalogue of SNPs as a reference panel to impute missing genotypes at untyped loci in data sets derived from lower density genotyping tools (150 K GBS-derived SNPs/530 samples). After imputation, 96.4% of the missing genotypes imputed in this fashion proved to be accurate. Using a combination of three bioinformatics pipelines, we uncovered ~92 K SVs (deletions, insertions, inversions, duplications, CNVs and translocations) and estimated that over 90% of these were accurate. Finally, we noticed that the duplication of certain genomic regions explained much of the residual heterozygosity at SNP loci in otherwise highly inbred soya bean accessions. This is the first time that a comprehensive description of both SNP haplotype diversity and SV has been achieved within a regionally relevant subset of a major crop.

The α' subunit of ß-conglycinin and various glycinin subunits of soy are not required to modulate hepatic lipid metabolism in rats.

PURPOSE: This study examined the effect of soy proteins with depletion of different subunits of the two major storage proteins, ß-conglycinin and glycinin, on hepatic lipids and proteins involved in lipid metabolism in rats, since the bioactive component of soy responsible for lipid-lowering is unclear. METHODS: Weanling Sprague Dawley rats were fed diets containing either 20% casein protein in the absence (casein) or presence (casein + ISF) of isoflavones or 20% alcohol-washed soy protein isolate (SPI) or 20% soy protein concentrates derived from a conventional (Haro) or 2 soybean lines lacking the α' subunit of ß-conglycinin and the A1-3 (1TF) or A1-5 (1a) subunits of glycinin. After 8 weeks, the rats were necropsied and liver proteins and lipids were extracted and analysed. RESULTS: The results showed that soy protein diets reduced lipid droplet accumulation and content in the liver compared to casein diets. The soy protein diets also decreased the level of hepatic mature SREBP-1 and FAS in males, with significant decreases in diets 1TF and 1a compared to the casein diets. The effect of the soy protein diets on female hepatic mature SREBP-1, FAS, and HMGCR was confounded since casein + ISF decreased these levels compared to casein alone perhaps muting the decrease by soy protein. A reduction in both phosphorylated and total STAT3 in female livers by ISF may account for the gender difference in mechanism in the regulation and protein expression of the lipid modulators. CONCLUSIONS: Overall, soy protein deficient in the α' subunit of ß-conglycinin and A1-5 subunits of glycinin maintain similar hypolipidemic function compared to the conventional soy protein. The exact bioactive component(s) warrant identification.

Mapping and identification of a potential candidate gene for a novel maturity locus, E10, in soybean.

KEY MESSAGE: E10 is a new maturity locus in soybean and FT4 is the predicted/potential functional gene underlying the locus. Flowering and maturity time traits play crucial roles in economic soybean production. Early maturity is critical for north and west expansion of soybean in Canada. To date, 11 genes/loci have been identified which control time to flowering and maturity; however, the molecular bases of almost half of them are not yet clear. We have identified a new maturity locus called "E10" located at the end of chromosome Gm08. The gene symbol E10e10 has been approved by the Soybean Genetics Committee. The e10e10 genotype results in 5-10 days earlier maturity than E10E10. A set of presumed E10E10 and e10e10 genotypes was used to identify contrasting SSR and SNP haplotypes. These haplotypes, and their association with maturity, were maintained through five backcross generations. A functional genomics approach using a predicted protein-protein interaction (PPI) approach (Protein-protein Interaction Prediction Engine, PIPE) was used to investigate approximately 75 genes located in the genomic region that SSR and SNP analyses identified as the location of the E10 locus. The PPI analysis identified FT4 as the most likely candidate gene underlying the E10 locus. Sequence analysis of the two FT4 alleles identified three SNPs, in the 5'UTR, 3'UTR and fourth exon in the coding region, which result in differential mRNA structures. Allele-specific markers were developed for this locus and are available for soybean breeders to efficiently develop earlier maturing cultivars using molecular marker assisted breeding.

Identification of loci governing eight agronomic traits using a GBS-GWAS approach and validation by QTL mapping in soya bean.

Soya bean is a major source of edible oil and protein for human consumption as well as animal feed. Understanding the genetic basis of different traits in soya bean will provide important insights for improving breeding strategies for this crop. A genome-wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of agronomic traits in soya bean. A genotyping-by-sequencing (GBS) approach was used to provide dense genome-wide marker coverage (>47,000 SNPs) for a panel of 304 short-season soya bean lines. A subset of 139 lines, representative of the diversity among these, was characterized phenotypically for eight traits under six environments (3 sites × 2 years). Marker coverage proved sufficient to ensure highly significant associations between the genes known to control simple traits (flower, hilum and pubescence colour) and flanking SNPs. Between one and eight genomic loci associated with more complex traits (maturity, plant height, seed weight, seed oil and protein) were also identified. Importantly, most of these GWAS loci were located within genomic regions identified by previously reported quantitative trait locus (QTL) for these traits. In some cases, the reported QTLs were also successfully validated by additional QTL mapping in a biparental population. This study demonstrates that integrating GBS and GWAS can be used as a powerful complementary approach to classical biparental mapping for dissecting complex traits in soya bean.

The α' subunit of ß-conglycinin and the A1-5 subunits of glycinin are not essential for many hypolipidemic actions of dietary soy proteins in rats.

PURPOSE: This study examined the effects of dietary soy protein (SP) lacking different storage protein subunits and isoflavones (ISF) on the abdominal fat, blood lipids, thyroid hormones, and enzymatic activities in rats. METHODS: Weanling Sprague-Dawley rats (8 males and 8 females/group) were fed diets containing either 20 % casein without or with supplemental isoflavones or alcohol-washed SP isolate or SP concentrates (SPC) prepared from 6 different soy bean lines for 8 weeks. RESULTS: Feeding of diets containing SPC regardless of their subunit compositions significantly lowered relative liver weights, blood total, free, and LDL cholesterol in both genders (P < 0.05) and also reduced serum free fatty acids (FFA) and abdominal fat in females (P < 0.05) compared to the casein or casein + ISF diets. Dietary SPC significantly elevated the plasma free triiodothyronine (T3) in both genders and total T3 in females compared to the casein diet (P < 0.05). The SPC lacking ß-conglycinin α' and either the glycinin A1-3 or A1-5 subunits increased total T3 in males and reduced plasma enzymatic activities of creatine kinase and lactate dehydrogenase compared to casein or casein + ISF diet (P < 0.05). CONCLUSIONS: Soy isoflavones were mainly responsible for the hypocholesterolemic effects and increased plasma free T3, whereas reduction in FFA, abdominal fat, liver weight and increased plasma total T3 were the effects of the soy proteins. Neither the α' subunit of ß-conglycinin nor the A1-5 subunits of glycinin are essential for the hypolipidemic properties of soy proteins.

Using the candidate gene approach for detecting genes underlying seed oil concentration and yield in soybean.

Increasing the oil concentration in soybean seeds has been given more attention in recent years because of demand for both edible oil and biodiesel production. Oil concentration in soybean is a complex quantitative trait regulated by many genes as well as environmental conditions. To identify genes governing seed oil concentration in soybean, 16 putative candidate genes of three important gene families (GPAT: acyl-CoA:sn-glycerol-3-phosphate acyltransferase, DGAT: acyl-CoA:diacylglycerol acyltransferase, and PDAT: phospholipid:diacylglycerol acyltransferase) involved in triacylglycerol (TAG) biosynthesis pathways were selected and their sequences retrieved from the soybean database ( http://www.phytozome.net/soybean ). Three sequence mutations were discovered in either coding or noncoding regions of three DGAT soybean isoforms when comparing the parents of a 203 recombinant inbreed line (RIL) population; OAC Wallace and OAC Glencoe. The RIL population was used to study the effects of these mutations on seed oil concentration and other important agronomic and seed composition traits, including seed yield and protein concentration across three field locations in Ontario, Canada, in 2009 and 2010. An insertion/deletion (indel) mutation in the GmDGAT2B gene in OAC Wallace was significantly associated with reduced seed oil concentration across three environments and reduced seed yield at Woodstock in 2010. A mutation in the 3' untranslated (3'UTR) region of GmDGAT2C was associated with seed yield at Woodstock in 2009. A mutation in the intronic region of GmDGAR1B was associated with seed yield and protein concentration at Ottawa in 2010. The genes identified in this study had minor effects on either seed yield or oil concentration, which was in agreement with the quantitative nature of the traits. However, the novel gene-specific markers designed in the present study can be used in soybean breeding for marker-assisted selection aimed at increasing seed yield and oil concentration with no significant impact on seed protein concentration.

Genetic control of soybean seed oil: II. QTL and genes that increase oil concentration without decreasing protein or with increased seed yield.

Soybean [Glycine max (L.) Merrill] seed oil is the primary global source of edible oil and a major renewable and sustainable feedstock for biodiesel production. Therefore, increasing the relative oil concentration in soybean is desirable; however, that goal is complex due to the quantitative nature of the oil concentration trait and possible effects on major agronomic traits such as seed yield or protein concentration. The objectives of the present study were to study the relationship between seed oil concentration and important agronomic and seed quality traits, including seed yield, 100-seed weight, protein concentration, plant height, and days to maturity, and to identify oil quantitative trait loci (QTL) that are co-localized with the traits evaluated. A population of 203 F4:6 recombinant inbred lines, derived from a cross between moderately high oil soybean genotypes OAC Wallace and OAC Glencoe, was developed and grown across multiple environments in Ontario, Canada, in 2009 and 2010. Among the 11 QTL associated with seed oil concentration in the population, which were detected using either single-factor ANOVA or multiple QTL mapping methods, the number of QTL that were co-localized with other important traits QTL were six for protein concentration, four for seed yield, two for 100-seed weight, one for days to maturity, and one for plant height. The oil-beneficial allele of the QTL tagged by marker Sat_020 was positively associated with seed protein concentration. The oil favorable alleles of markers Satt001 and GmDGAT2B were positively correlated with seed yield. In addition, significant two-way epistatic interactions, where one of the interacting markers was solely associated with seed oil concentration, were identified for the selected traits in this study. The number of significant epistatic interactions was seven for yield, four for days to maturity, two for 100-seed weight, one for protein concentration, and one for plant height. The identified molecular markers associated with oil-related QTL in this study, which also have positive effects on other important traits such as seed yield and protein concentration, could be used in the soybean marker breeding programs aimed at developing either higher seed yield and oil concentration or higher seed protein and oil concentration per hectare. Alternatively, selecting complementary parents with greater breeding values due to positive epistatic interactions could lead to the development of higher oil soybean cultivars.

Genetic control of soybean seed oil: I. QTL and genes associated with seed oil concentration in RIL populations derived from crossing moderately high-oil parents.

Soybean seed is a major source of oil for human consumption worldwide and the main renewable feedstock for biodiesel production in North America. Increasing seed oil concentration in soybean [Glycine max (L.) Merrill] with no or minimal impact on protein concentration could be accelerated by exploiting quantitative trait loci (QTL) or gene-specific markers. Oil concentration in soybean is a polygenic trait regulated by many genes with mostly small effects and which is negatively associated with protein concentration. The objectives of this study were to discover and validate oil QTL in two recombinant inbred line (RIL) populations derived from crosses between three moderately high-oil soybean cultivars, OAC Wallace, OAC Glencoe, and RCAT Angora. The RIL populations were grown across several environments over 2 years in Ontario, Canada. In a population of 203 F(3:6) RILs from a cross of OAC Wallace and OAC Glencoe, a total of 11 genomic regions on nine different chromosomes were identified as associated with oil concentration using multiple QTL mapping and single-factor ANOVA. The percentage of the phenotypic variation accounted for by each QTL ranged from 4 to 11 %. Of the five QTL that were tested in a population of 211 F(3:5) RILs from the cross RCAT Angora × OAC Wallace, a "trait-based" bidirectional selective genotyping analysis validated four QTL (80 %). In addition, a total of seven two-way epistatic interactions were identified for oil concentration in this study. The QTL and epistatic interactions identified in this study could be used in marker-assisted introgression aimed at pyramiding high-oil alleles in soybean cultivars to increase oil concentration for biodiesel as well as edible oil applications.

P-TarPmiR accurately predicts plant-specific miRNA targets.

microRNAs (miRNAs) are small non-coding ribonucleic acids that post-transcriptionally regulate gene expression through the targeting of messenger RNA (mRNAs). Most miRNA target predictors have focused on animal species and prediction performance drops substantially when applied to plant species. Several rule-based miRNA target predictors have been developed in plant species, but they often fail to discover new miRNA targets with non-canonical miRNA-mRNA binding. Here, the recently published TarDB database of plant miRNA-mRNA data is leveraged to retrain the TarPmiR miRNA target predictor for application on plant species. Rigorous experiment design across four plant test species demonstrates that animal-trained predictors fail to sustain performance on plant species, and that the use of plant-specific training data improves accuracy depending on the quantity of plant training data used. Surprisingly, our results indicate that the complete exclusion of animal training data leads to the most accurate plant-specific miRNA target predictor indicating that animal-based data may detract from miRNA target prediction in plants. Our final plant-specific miRNA prediction method, dubbed P-TarPmiR, is freely available for use at http://ptarpmir.cu-bic.ca . The final P-TarPmiR method is used to predict targets for all miRNA within the soybean genome. Those ranked predictions, together with GO term enrichment, are shared with the research community.

Species-specific microRNA discovery and target prediction in the soybean cyst nematode.

The soybean cyst nematode (SCN) is a devastating pathogen for economic and food security considerations. Although the SCN genome has recently been sequenced, the presence of any miRNA has not been systematically explored and reported. This paper describes the development of a species-specific SCN miRNA discovery pipeline and its application to the SCN genome. Experiments on well-documented model nematodes (Caenorhabditis elegans and Pristionchus pacificus) are used to tune the pipeline's hyperparameters and confirm its recall and precision. Application to the SCN genome identifies 3342 high-confidence putative SCN miRNA. Prediction specificity within SCN is confirmed by applying the pipeline to RNA hairpins from known exonic regions of the SCN genome (i.e., sequences known to not be miRNA). Prediction recall is confirmed by building a positive control set of SCN miRNA, based on a limited deep sequencing experiment. Interestingly, a number of novel miRNA are predicted to be encoded within the intronic regions of effector genes, known to be involved in SCN parasitism, suggesting that these miRNA may also be involved in the infection process or virulence. Beyond miRNA discovery, gene targets within SCN are predicted for all high-confidence novel miRNA using a miRNA:mRNA target prediction system. Lastly, cross-kingdom miRNA targeting is investigated, where putative soybean mRNA targets are identified for novel SCN miRNA. All predicted miRNA and gene targets are made available in appendix and through a Borealis DataVerse open repository ( https://borealisdata.ca/dataset.xhtml?persistentId=doi:10.5683/SP3/30DEXA ).

Large-scale data mining pipeline for identifying novel soybean genes involved in resistance against the soybean cyst nematode.

The soybean cyst nematode (SCN) [Heterodera glycines Ichinohe] is a devastating pathogen of soybean [Glycine max (L.) Merr.] that is rapidly becoming a global economic issue. Two loci conferring SCN resistance have been identified in soybean, Rhg1 and Rhg4; however, they offer declining protection. Therefore, it is imperative that we identify additional mechanisms for SCN resistance. In this paper, we develop a bioinformatics pipeline to identify protein-protein interactions related to SCN resistance by data mining massive-scale datasets. The pipeline combines two leading sequence-based protein-protein interaction predictors, the Protein-protein Interaction Prediction Engine (PIPE), PIPE4, and Scoring PRotein INTeractions (SPRINT) to predict high-confidence interactomes. First, we predicted the top soy interacting protein partners of the Rhg1 and Rhg4 proteins. Both PIPE4 and SPRINT overlap in their predictions with 58 soybean interacting partners, 19 of which had GO terms related to defense. Beginning with the top predicted interactors of Rhg1 and Rhg4, we implement a "guilt by association" in silico proteome-wide approach to identify novel soybean genes that may be involved in SCN resistance. This pipeline identified 1,082 candidate genes whose local interactomes overlap significantly with the Rhg1 and Rhg4 interactomes. Using GO enrichment tools, we highlighted many important genes including five genes with GO terms related to response to the nematode (GO:0009624), namely, Glyma.18G029000, Glyma.11G228300, Glyma.08G120500, Glyma.17G152300, and Glyma.08G265700. This study is the first of its kind to predict interacting partners of known resistance proteins Rhg1 and Rhg4, forming an analysis pipeline that enables researchers to focus their search on high-confidence targets to identify novel SCN resistance genes in soybean.

Genome-wide analysis of cold imbibition stress in soybean, Glycine max.

In Canada, the length of the frost-free season necessitates planting crops as early as possible to ensure that the plants have enough time to reach full maturity before they are harvested. Early planting carries inherent risks of cold water imbibition (specifically less than 4°C) affecting seed germination. A marker dataset developed for a previously identified Canadian soybean GWAS panel was leveraged to investigate the effect of cold water imbibition on germination. Seed from a panel of 137 soybean elite cultivars, grown in the field at Ottawa, ON, over three years, were placed on filter paper in petri dishes and allowed to imbibe water for 16 hours at either 4°C or 20°C prior to being transferred to a constant 20°C. Observations on seed germination, defined as the presence of a 1 cm radicle, were done from day two to seven. A three-parameter exponential rise to a maximum equation (3PERM) was fitted to estimate germination, time to the one-half maximum germination, and germination uniformity for each cultivar. Genotype-by-sequencing was used to identify SNPs in 137 soybean lines, and using genome-wide association studies (GWAS - rMVP R package, with GLM, MLM, and FarmCPU as methods), haplotype block analysis, and assumed linkage blocks of ±100 kbp, a threshold for significance was established using the qvalue package in R, and five significant SNPs were identified on chromosomes 1, 3, 4, 6, and 13 for maximum germination after cold water imbibition. Percent of phenotypic variance explained (PVE) and allele substitution effect (ASE) eliminated two of the five candidate SNPs, leaving three QTL regions on chromosomes 3, 6, and 13 (Chr3-3419152, Chr6-5098454, and Chr13-29649544). Based on the gene ontology (GO) enrichment analysis, 14 candidate genes whose function is predicted to include germination and cold tolerance related pathways were identified as candidate genes. The identified QTLs can be used to select future soybean cultivars tolerant to cold water imbibition and mitigate risks associated with early soybean planting.

Differential gene expression provides leads to environmentally regulated soybean seed protein content.

Soybean is an important global source of plant-based protein. A persistent trend has been observed over the past two decades that soybeans grown in western Canada have lower seed protein content than soybeans grown in eastern Canada. In this study, 10 soybean genotypes ranging in average seed protein content were grown in an eastern location (control) and three western locations (experimental) in Canada. Seed protein and oil contents were measured for all lines in each location. RNA-sequencing and differential gene expression analysis were used to identify differentially expressed genes that may account for relatively low protein content in western-grown soybeans. Differentially expressed genes were enriched for ontologies and pathways that included amino acid biosynthesis, circadian rhythm, starch metabolism, and lipid biosynthesis. Gene ontology, pathway mapping, and quantitative trait locus (QTL) mapping collectively provide a close inspection of mechanisms influencing nitrogen assimilation and amino acid biosynthesis between soybeans grown in the East and West. It was found that western-grown soybeans had persistent upregulation of asparaginase (an asparagine hydrolase) and persistent downregulation of asparagine synthetase across 30 individual differential expression datasets. This specific difference in asparagine metabolism between growing environments is almost certainly related to the observed differences in seed protein content because of the positive correlation between seed protein content at maturity and free asparagine in the developing seed. These results provided pointed information on seed protein-related genes influenced by environment. This information is valuable for breeding programs and genetic engineering of geographically optimized soybeans.

Relationship between asparagine metabolism and protein concentration in soybean seed.

The relationship between asparagine metabolism and protein concentration was investigated in soybean seed. Phenotyping of a population of recombinant inbred lines adapted to Illinois confirmed a positive correlation between free asparagine levels in developing seeds and protein concentration at maturity. Analysis of a second population of recombinant inbred lines adapted to Ontario associated the elevated free asparagine trait with two of four quantitative trait loci determining population variation for protein concentration, including a major one on chromosome 20 (linkage group I) which has been reported in multiple populations. In the seed coat, levels of asparagine synthetase were high at 50 mg and progressively declined until 150 mg seed weight, suggesting that nitrogenous assimilates are pre-conditioned at early developmental stages to enable a high concentration of asparagine in the embryo. The levels of asparaginase B1 showed an opposite pattern, being low at 50 mg and progressively increased until 150 mg, coinciding with an active phase of storage reserve accumulation. In a pair of genetically related cultivars, ∼2-fold higher levels of asparaginase B1 protein and activity in seed coat, were associated with high protein concentration, reflecting enhanced flux of nitrogen. Transcript expression analyses attributed this difference to a specific asparaginase gene, ASPGB1a. These results contribute to our understanding of the processes determining protein concentration in soybean seed.