Soybean Yield Formation Physiology - A Foundation for Precision Breeding Based Improvement

QTL mapping and genomic selection of stem and branch diameter in soybean (Glycine max L.)

Introduction

Soybean stem diameter (SD) and branch diameter (BD) are closely related traits, and genetic clarification of SD and BD is crucial for soybean breeding.

Methods

SD and BD were genetically analyzed by a population of 363 RIL derived from the cross between Zhongdou41 (ZD41) and ZYD02878 using restricted two-stage multi-locus genome-wide association, inclusive composite interval mapping, and three-variance component multi-locus random SNP effect mixed linear modeling. Then candidate genes of major QTLs were selected and genetic selection model of SD and BD were constructed respectively.

Results and discussion

The results showed that SD and BD were significantly correlated (r = 0.74, P < 0.001). A total of 93 and 84 unique quantitative trait loci (QTL) were detected for SD and BD, respectively by three different methods. There were two and ten major QTLs for SD and BD, respectively, with phenotypic variance explained (PVE) by more than 10%. Within these loci, seven genes involved in the regulation of phytohormones (IAA and GA) and cell proliferation and showing extensive expression of shoot apical meristematic genes were selected as candidate genes. Genomic selection (GS) analysis showed that the trait-associated markers identified in this study reached 0.47-0.73 in terms of prediction accuracy, which was enhanced by 6.56-23.69% compared with genome-wide markers. These results clarify the genetic basis of SD and BD, which laid solid foundation in regulation gene cloning, and GS models constructed could be potentially applied in future breeding programs.

Optimization of soybean physiochemical, agronomic, and genetic responses under varying regimes of day and night temperatures

Soybean is an important oilseed crop worldwide; however, it has a high sensitivity to temperature variation, particularly at the vegetative stage to the pod-filling stage. Temperature change affects physiochemical and genetic traits regulating the soybean agronomic yield. In this regard, the current study aimed to comparatively evaluate the effects of varying regimes of day and night temperatures (T1 = 20°C/12°C, T2 = 25°C/17°C, T3 = 30°C/22°C, T4 = 35°C/27°C, and T5 = 40°C/32°C) on physiological (chlorophyll, photosynthesis, stomatal conductance, transpiration, and membrane damage) biochemical (proline and antioxidant enzymes), genetic (GmDNJ1, GmDREB1G;1, GmHSF-34, GmPYL21, GmPIF4b, GmPIP1;6, GmGBP1, GmHsp90A2, GmTIP2;6, and GmEF8), and agronomic traits (pods per plant, seeds per plant, pod weight per plant, and seed yield per plant) of soybean cultivars (Swat-84 and NARC-1). The experiment was performed in soil plant atmosphere research (SPAR) units using two factorial arrangements with cultivars as one factor and temperature treatments as another factor. A significant increase in physiological, biochemical, and agronomic traits with increased gene expression was observed in both soybean cultivars at T4 (35°C/27°C) as compared to below and above regimes of temperatures. Additionally, it was established by correlation, principal component analysis (PCA), and heatmap analysis that the nature of soybean cultivars and the type of temperature treatments have a significant impact on the paired association of agronomic and biochemical traits, which in turn affects agronomic productivity. Furthermore, at corresponding temperature regimes, the expression of the genes matched the expression of physiochemical traits. The current study has demonstrated through extensive physiochemical, genetic, and biochemical analyses that the ideal day and night temperature for soybeans is T4 (35°C/27°C), with a small variation having a significant impact on productivity from the vegetative stage to the grain-filling stage.

Disturbed electron transport beyond PSI changes metabolome and transcriptome in Zn-deficient soybean

Low Zn availability in soils is a problem in many parts of the world, with tremendous consequences for food and feed production because Zn deficiency affects the yield and quality of plants. In this study we investigated the consequences of Zn-limitation in hydroponically cultivated soybean (Glycine max L.) plants. Parameters of photosynthesis biophysics were determined by spatially and spectrally resolved Kautsky and OJIP fluorescence kinetics and oxygen production at two time points (V4 stage, after five weeks, and pod development stage, R5-R6, after 8-10 weeks). Lower NPQ at 730 nm and lower quantum yield of electron transport flux until PSI acceptors were observed, indicating an inhibition of the PSI acceptor side. Metalloproteomics showed that down-regulation of Cu/Zn-superoxide dismutase (CuZnSOD) and Zn‑carbonic anhydrase (CA) were primary consequences of Zn-limitation. This explained the effects on photosynthesis in terms of decreased use of excitons, which consequently led to oxidative stress. Indeed, untargeted metabolomics revealed an accumulation of lipid oxidation products in the Zn-deficient leaves. Further response to Zn deficiency included up-regulation of gene expression of cell wall metabolism, response to (a)biotic stressors and antioxidant activity, which correlated with accumulation of antioxidants, Vit B6, (iso)flavonoids and phytoalexins.

A global-temporal analysis on Phytophthora sojae resistance-gene efficacy

Plant disease resistance genes are widely used in agriculture to reduce disease outbreaks and epidemics and ensure global food security. In soybean, Rps (Resistance to Phytophthora sojae) genes are used to manage Phytophthora sojae, a major oomycete pathogen that causes Phytophthora stem and root rot (PRR) worldwide. This study aims to identify temporal changes in P. sojae pathotype complexity, diversity, and Rps gene efficacy. Pathotype data was collected from 5121 isolates of P. sojae, derived from 29 surveys conducted between 1990 and 2019 across the United States, Argentina, Canada, and China. This systematic review shows a loss of efficacy of specific Rps genes utilized for disease management and a significant increase in the pathotype diversity of isolates over time. This study finds that the most widely deployed Rps genes used to manage PRR globally, Rps1a, Rps1c and Rps1k, are no longer effective for PRR management in the United States, Argentina, and Canada. This systematic review emphasizes the need to widely introduce new sources of resistance to P. sojae, such as Rps3a, Rps6, or Rps11, into commercial cultivars to effectively manage PRR going forward.

Polyploidy and diploidization in soybean

Polyploidy is widespread and particularly common in angiosperms. The prevalence of polyploidy in the plant suggests it as a crucial driver of diversification and speciation. The paleopolyploid soybean (Glycine max) is one of the most important crops of plant protein and oil for humans and livestock. Soybean experienced two rounds of whole genome duplication around 13 and 59 million years ago. Due to the relatively slow process of post-polyploid diploidization, most genes are present in multiple copies across the soybean genome. Growing evidence suggests that polyploidization and diploidization could cause rapid and dramatic changes in genomic structure and epigenetic modifications, including gene loss, transposon amplification, and reorganization of chromatin architecture. This review is focused on recent progresses about genetic and epigenetic changes during polyploidization and diploidization of soybean and represents the challenges and potentials for application of polyploidy in soybean breeding.

Evaluating the Effectiveness of Calcium Silicate in Enhancing Soybean Growth and Yield

The application of silicon (Si) fertilizer positively impacts crop health, yield, and seed quality worldwide. Si is a "quasi-essential" element that is crucial for plant nutrition and stress response but is less associated with growth. This study aimed to investigate the effect of Si on the yield of cultivated soybean (Glycine max L). Two locations, Gyeongsan and Gunwi, in the Republic of Korea were selected, and a land suitability analysis was performed using QGIS version 3.28.1. The experiments at both locations consisted of three treatments: the control, Si fertilizer application at 2.3 kg per plot (9 m × 9 m) (T1), and Si fertilizer application at 4.6 kg per plot (9 m × 9 m) (T2). The agronomic, root, and yield traits, as well as vegetative indices, were analyzed to evaluate the overall impact of Si. The results demonstrated that Si had consistently significant effects on most root and shoot parameters in the two experimental fields, which led to significantly increased crop yield when compared with the control, with T2 (22.8% and 25.6%, representing an output of 2.19 and 2.24 t ha^-1 at Gyeongsan and Gunwi, respectively) showing a higher yield than T1 (11% and 14.2%, representing 1.98 and 2.04 t ha^-1 at Gyeongsan and Gunwi, respectively). These results demonstrate the positive impact of exogenous Si application on the overall growth, morphological and physiological traits, and yield output of soybeans. However, the application of the optimal concentration of Si according to the crop requirement, soil status, and environmental conditions requires further studies.

Shade avoidance syndrome in soybean and ideotype toward shade tolerance

The shade avoidance syndrome (SAS) in soybean can have destructive effects on yield, as essential carbon resources reserved for yield are diverted to the petiole and stem for exaggerated elongation, resulting in lodging and susceptibility to disease. Despite numerous attempts to reduce the unfavorable impacts of SAS for the development of cultivars suitable for high-density planting or intercropping, the genetic bases and fundamental mechanisms of SAS remain largely unclear. The extensive research conducted in the model plant Arabidopsis provides a framework for understanding the SAS in soybean. Nevertheless, recent investigations suggest that the knowledge obtained from model Arabidopsis may not be applicable to all processes in soybean. Consequently, further efforts are required to identify the genetic regulators of SAS in soybean for molecular breeding of high-yield cultivars suitable for density farming. In this review, we present an overview of the recent developments in SAS studies in soybean and suggest an ideal planting architecture for shade-tolerant soybean intended for high-yield breeding.

Identification of genomic regions associated with the plasticity of carbon 13 ratio in soybean

Improving water use efficiency (WUE) for soybean [Glycine max (L.) Merr.] through selection for high carbon isotope (C13) ratio may increase drought tolerance, but increased WUE may limit growth in productive environments. An ideal genotype would be plastic for C13 ratio; that is, be able to alter C13 ratio in response to the environment. Our objective was to identify genomic regions associated with C13 ratio plasticity, C13 ratio stability, and overall C13 ratio in two panels of diverse Maturity Group IV soybean accessions. A second objective was to identify accessions that differed in their C13 ratio plasticity. Panel 1 (205 accessions) was evaluated in seven irrigated and four drought environments, and Panel 2 (373 accessions) was evaluated in four environments. Plasticity was quantified as the slope from regressing C13 ratio of individual genotypes against an environmental index calculated based on the mean within and across environments. The regression intercept was considered a measure of C13 ratio over all environments, and the root mean square error was considered a measure of stability. Combined over both panels, genome-wide association mapping (GWAM) identified 19 single nucleotide polymorphisms (SNPs) for plasticity, 39 SNPs for C13 ratio, and 16 SNPs for stability. Among these SNPs, 71 candidate genes had annotations associated with transpiration or water conservation and transport, root development, root hair elongation, and stomatal complex morphogenesis. The genomic regions associated with plasticity and stability identified in the current study will be a useful resource for implementing genomic selection for improving drought tolerance in soybean.

Effects of Co-Inoculating Saccharomyces spp. with Bradyrhizobium japonicum on Atmospheric Nitrogen Fixation in Soybeans (Glycine max (L.))

Crop production encounters challenges due to the dearth of nitrogen (N) and phosphorus (P), while excessive chemical fertilizer use causes environmental hazards. The use of N-fixing microbes and P-solubilizing microbes (PSMs) can be a sustainable strategy to overcome these problems. Here, we conducted a greenhouse pot experiment following a completely randomized blocked design to elucidate the influence of co-inoculating N-fixing bacteria (Bradyrhizobium japonicum) and PSMs (Saccharomyces cerevisiae and Saccharomyces exiguus) on atmospheric N₂-fixation, growth, and yield. The results indicate a significant influence of interaction on Indole-3-acetic acid production, P solubilization, seedling germination, and growth. It was also found that atmospheric N₂-fixation, nodule number per plant, nodule dry weight, straw, and root dry weight per plant at different growth stages were significantly increased under dual inoculation treatments relative to single inoculation or no inoculation treatment. Increased seed yield and N and P accumulation were also noticed under co-inoculation treatments. Soil available N was highest under sole bacterial inoculation and lowest under the control treatment, while soil available P was highest under co-inoculation treatments and lowest under the control treatment. We demonstrated that the co-inoculation of N-fixing bacteria and PSMs enhances P bioavailability and atmospheric N₂-fixation in soybeans leading to improved soil fertility, raising crop yields, and promoting sustainable agriculture.

Identification of genomic regions associated with soybean responses to off-target dicamba exposure

The widespread adoption of genetically modified (GM) dicamba-tolerant (DT) soybean was followed by numerous reports of off-target dicamba damage and yield losses across most soybean-producing states. In this study, a subset of the USDA Soybean Germplasm Collection consisting of 382 genetically diverse soybean accessions originating from 15 countries was used to identify genomic regions associated with soybean response to off-target dicamba exposure. Accessions were genotyped with the SoySNP50K BeadChip and visually screened for damage in environments with prolonged exposure to off-target dicamba. Two models were implemented to detect significant marker-trait associations: the Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK) and a model that allows the inclusion of population structure in interaction with the environment (G×E) to account for variable patterns of genotype responses in different environments. Most accessions (84%) showed a moderate response, either moderately tolerant or moderately susceptible, with approximately 8% showing tolerance and susceptibility. No differences in off-target dicamba damage were observed across maturity groups and centers of origin. Both models identified significant associations in regions of chromosomes 10 and 19. The BLINK model identified additional significant marker-trait associations on chromosomes 11, 14, and 18, while the G×E model identified another significant marker-trait association on chromosome 15. The significant SNPs identified by both models are located within candidate genes possessing annotated functions involving different phases of herbicide detoxification in plants. These results entertain the possibility of developing non-GM soybean cultivars with improved tolerance to off-target dicamba exposure and potentially other synthetic auxin herbicides. Identification of genetic sources of tolerance and genomic regions conferring higher tolerance to off-target dicamba may sustain and improve the production of other non-DT herbicide soybean production systems, including the growing niche markets of organic and conventional soybean.

Combined Transcriptome and Metabolome Profiling Provide Insights into Cold Responses in Rapeseed (Brassica napus L.) Genotypes with Contrasting Cold-Stress Sensitivity

Low temperature is a major environmental factor, which limits rapeseed (Brassica napus L.) growth, development, and productivity. So far, the physiological and molecular mechanisms of rapeseed responses to cold stress are not fully understood. Here, we explored the transcriptome and metabolome profiles of two rapeseed genotypes with contrasting cold responses, i.e., XY15 (cold-sensitive) and GX74 (cold-tolerant). The global metabolome profiling detected 545 metabolites in siliques of both genotypes before (CK) and after cold-stress treatment (LW). The contents of several sugar metabolites were affected by cold stress with the most accumulated saccharides being 3-dehydro-L-threonic acid, D-xylonic acid, inositol, D-mannose, D-fructose, D-glucose, and L-glucose. A total of 1943 and 5239 differentially expressed genes were identified from the transcriptome sequencing in XY15CK_vs_XY15LW and GX74CK_vs_GX74LW, respectively. We observed that genes enriched in sugar metabolism and biosynthesis-related pathways, photosynthesis, reactive oxygen species scavenging, phytohormone, and MAPK signaling were highly expressed in GX74LW. In addition, several genes associated with cold-tolerance-related pathways, e.g., the CBF-COR pathway and MAPK signaling, were specifically expressed in GX74LW. Contrarily, genes in the above-mentioned pathways were mostly downregulated in XY15LW. Thus, our results indicate the involvement of these pathways in the differential cold-stress responses in XY15 and GX74.

An APETALA2/ethylene responsive factor transcription factor GmCRF4a regulates plant height and auxin biosynthesis in soybean

Plant height is one of the key agronomic traits affecting soybean yield. The cytokinin response factors (CRFs), as a branch of the APETALA2/ethylene responsive factor (AP2/ERF) super gene family, have been reported to play important roles in regulating plant growth and development. However, their functions in soybean remain unknown. This study characterized a soybean CRF gene named GmCRF4a by comparing the performance of the homozygous Gmcrf4a-1 mutant, GmCRF4a overexpression (OX) and co-silencing (CS) lines. Phenotypic analysis showed that overexpression of GmCRF4a resulted in taller hypocotyls and epicotyls, more main stem nodes, and higher plant height. While down-regulation of GmCRF4a conferred shorter hypocotyls and epicotyls, as well as a reduction in plant height. The histological analysis results demonstrated that GmCRF4a promotes epicotyl elongation primarily by increasing cell length. Furthermore, GmCRF4a is required for the expression of GmYUCs genes to elevate endogenous auxin levels, which may subsequently enhance stem elongation. Taken together, these observations describe a novel regulatory mechanism in soybean, and provide the basis for elucidating the function of GmCRF4a in auxin biosynthesis pathway and plant heigh regulation in plants.

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

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