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.

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 α' 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.

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.

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.

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.

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

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.

Photoperiod Affects Node Appearance Rate and Flowering in Early Maturing Soybean.

The photoperiod plays a critical role in the control of flowering timing in soybean (Glycine max (L.) Merr.) with long days increasing the time to flowering. Early flowering cultivars have been developed from breeding programs for environments with long photoperiods; however, this effect is challenging to isolate in field experiments because of other environmental influences. Our experiment examined the effect of photoperiod on the node appearance rate and time to flower for 13 early maturing soybean cultivars ranging in maturity group (MG) between 000.9 and 1.3. Growth chambers were programmed to 14, 15, 16, and 17 h photoperiods and temperature was kept at 25 °C. The date of emergence and main stem node appearance were recorded until flowering. The node appearance rate was slowest for the first node and increased thereafter. All cultivars required more time to flowering in the longer photoperiod treatments and the later rated MG had the greatest sensitivity to photoperiod. A delay in time to flower from a longer photoperiod can delay maturity and expose the crop to fall frost that can reduce seed yield and quality. Understanding and documentation of soybean photoperiod sensitivity will help plant breeders develop suitable cultivars for environments with long photoperiods.

A Broad Review of Soybean Research on the Ongoing Race to Overcome Soybean Cyst Nematode.

Plant pathogens greatly impact food security of the ever-growing human population. Breeding resistant crops is one of the most sustainable strategies to overcome the negative effects of these biotic stressors. In order to efficiently breed for resistant plants, the specific plant-pathogen interactions should be understood. Soybean is a short-day legume that is a staple in human food and animal feed due to its high nutritional content. Soybean cyst nematode (SCN) is a major soybean stressor infecting soybean worldwide including in China, Brazil, Argentina, USA and Canada. There are many Quantitative Trait Loci (QTLs) conferring resistance to SCN that have been identified; however, only two are widely used: rhg1 and Rhg4. Overuse of cultivars containing these QTLs/genes can lead to SCN resistance breakdown, necessitating the use of additional strategies. In this manuscript, a literature review is conducted on research related to soybean resistance to SCN. The main goal is to provide a current understanding of the mechanisms of SCN resistance and list the areas of research that could be further explored.

Variability in Maturity, Oil and Protein Concentration, and Genetic Distinctness among Soybean Accessions Conserved at Plant Gene Resources of Canada.

Soybean (Glycine max (L.) Merr.) is one of the important crops in Canada and has the potential to expand its production further north into the Canadian Prairies. Such expansion, however, requires the search for adapted soybean germplasm useful for the development of productive cultivars with earlier maturity and increased protein concentration. We initiated several research activities to characterize 848 accessions of the soybean collection conserved at Plant Gene Resources of Canada (PGRC) for maturity, oil and protein concentration, and genetic distinctness. The characterization revealed a wide range of variations present in each assessed trait among the PGRC soybean accessions. The trait variabilities allowed for the identification of four core subsets of 35 PGRC soybean accessions, each specifically targeted for early maturity for growing in Saskatoon and Ottawa, and for high oil and protein concentration. The two early maturity core subsets for Saskatoon and Ottawa displayed days to maturity ranging from 103 to 126 days and 94 to 102 days, respectively. The two core subsets for high oil and protein concentration showed the highest oil and protein concentration from 25.0 to 22.7% and from 52.8 to 46.7%, respectively. However, these core subsets did not differ significantly in genetic distinctness (as measured with 19,898 SNP markers across 20 soybean chromosomes) from the whole PGRC soybean collection. These findings are useful, particularly for the management and utilization of the conserved soybean germplasm.

Novel QTL for Low Seed Cadmium Accumulation in Soybean.

Soybean is a valuable crop, used in animal feed and for human consumption. Selecting soybean cultivars with low seed cadmium (Cd) concentration is important for the purpose of minimizing the transfer of Cd into the human body. To ensure international trade, farmers need to produce soybean that meets the European Union (EU) Cd limit of 0.2 mg kg-1. In this study, we evaluated two populations of recombinant inbred lines (RILs), X5154 and X4050, for seed Cd accumulation. Linkage maps were constructed with 325 and 280 polymorphic simple sequence repeat (SSR) markers, respectively, and used to identify a novel minor quantitative trait locus (QTL) on chromosome 13 in the X4050 population between SSR markers Satt522 and Satt218. Based on a gene ontology search within the QTL region, seven genes were identified as candidates responsible for low seed Cd accumulation, including Glyma.13G308700 and Glyma.13G309100. In addition, we confirmed the known major gene, Cda1, in the X5154 population and developed KASP and CAPS/dCAPS allele-specific markers for efficient marker-assisted breeding for Cda1.

GmSWEET29 and Paralog GmSWEET34 Are Differentially Expressed between Soybeans Grown in Eastern and Western Canada.

Over the past two decades soybeans grown in western Canada have persistently had lower seed protein than those grown in eastern Canada. To understand the discrepancy in seed protein content between eastern- and western-grown soybeans, RNA-seq and differential expression analysis have been investigated. Ten soybean genotypes, ranging from low to high in seed protein content, were grown in four locations across eastern (Ottawa) and western (Morden, Brandon, and Saskatoon) Canada. Differential expression analysis revealed 34 differentially expressed genes encoding Glycine max Sugars Will Eventually be Exported Transporters (GmSWEETs), including paralogs GmSWEET29 and GmSWEET34 (AtSWEET2 homologs) that were consistently upregulated across all ten genotypes in each of the western locations over three years. GmSWEET29 and GmSWEET34 are likely candidates underlying the lower seed protein content of western soybeans. GmSWEET20 (AtSWEET12 homolog) was downregulated in the western locations and may also play a role in lower seed protein content. These findings are valuable for improving soybean agriculture in western growing regions, establishing more strategic and efficient agricultural practices.

Patterns of Genetic Variation in a Soybean Germplasm Collection as Characterized with Genotyping-by-Sequencing.

Genomic characterization is playing an increasing role in plant germplasm conservation and utilization, as it can provide higher resolution with genome-wide SNP markers than before to identify and analyze genetic variation. A genotyping-by-sequencing technique was applied to genotype 541 soybean accessions conserved at Plant Gene Resources of Canada and 30 soybean cultivars and breeding lines developed by the Ottawa soybean breeding program of Agriculture and Agri-Food Canada. The sequencing generated an average of 952,074 raw sequence reads per sample. SNP calling identified 43,891 SNPs across 20 soybean chromosomes and 69 scaffolds with variable levels of missing values. Based on 19,898 SNPs with up to 50% missing values, three distinct genetic groups were found in the assayed samples. These groups were a mixture of the samples that originated from different countries and the samples of known maturity groups. The samples that originated from Canada were clustered into all three distinct groups, but 30 Ottawa breeding lines fell into two groups only. Based on the average pairwise dissimilarity estimates, 40 samples with the most genetic distinctness were identified from three genetic groups with diverse sample origin and known maturity. Additionally, 40 samples with the highest genetic redundancy were detected and they consisted of different sample origins and maturity groups, largely from one genetic group. Moreover, some genetically duplicated samples were identified, but the overall level of genetic duplication was relatively low in the collection. These findings are useful for soybean germplasm management and utilization.