Non-CG DNA hypomethylation promotes photosynthesis and nitrogen fixation in soybean.

Non-CG DNA methylation, a plant-specific epigenetic mark mainly regulated by chromomethylase (CMT), is known to play important roles in Arabidopsis thaliana. However, whether and to what extent non-CG DNA methylation modulates agronomic traits in crops remain to be explored. Here, we describe the consequences of non-CG DNA hypomethylation on development, seed composition, and yield in soybean (Glycine max). We created a Gmcmt mutant line lacking function of all four CMT genes. This line exhibited substantial hypomethylation of non-CG (CHG and CHH) sites. Non-CG hypomethylation enhanced chromatin accessibility and promoted or repressed the expression of hundreds of functionally relevant genes, including upregulation of GOLDEN-LIKE 10 (GmGLK10), which led to enhanced photosynthesis and, unexpectedly, improved nitrogen fixation efficiency. The Gmcmt line produced larger seeds with increased protein content. This study provides insights into the mechanisms of non-CG methylation-based epigenetic regulation of soybean development and suggests viable epigenetic strategies for improving soybean yield and nutritional value.

Mutations in the genes responsible for the synthesis of furan fatty acids resolve the light-induced off-odor in soybean oil.

Soybean oil is the second most produced edible vegetable oil and is used for many edible and industrial materials. Unfortunately, it has the disadvantage of 'reversion flavor' under photooxidative conditions, which produces an off-odor and decreases the quality of edible oil. Reversion flavor and off-odor are caused by minor fatty acids in the triacylglycerol of soybean oil known as furan fatty acids, which produce 3-methyl-2,4-nonanedione (3-MND) upon photooxidation. As a solution to this problem, a reduction in furan fatty acids leads to a decrease in 3-MND, resulting in a reduction in the off-odor induced by light exposure. However, there are no reports on the genes related to the biosynthesis of furan fatty acids in soybean oil. In this study, four mutant lines showing low or no furan fatty acid levels in soybean seeds were isolated from a soybean mutant library. Positional cloning experiments and homology search analysis identified two genes responsible for furan fatty acid biosynthesis in soybean: Glyma.20G201400 and Glyma.04G054100. Ectopic expression of both genes produced furan fatty acids in transgenic soybean hairy roots. The structure of these genes is different from that of the furan fatty acid biosynthetic genes in photosynthetic bacteria. Homologs of these two group of genes are widely conserved in the plant kingdom. The purified oil from the furan fatty acid mutant lines had lower amounts of 3-MND and reduced off-odor after light exposure, compared with oil from the wild-type.

The impact of multifactorial stress combination on reproductive tissues and grain yield of a crop plant.

Global warming, climate change, and industrial pollution are altering our environment subjecting plants, microbiomes, and ecosystems to an increasing number and complexity of abiotic stress conditions, concurrently or sequentially. These conditions, termed, "multifactorial stress combination" (MFSC), can cause a significant decline in plant growth and survival. However, the impacts of MFSC on reproductive tissues and yield of major crop plants are largely unknown. We subjected soybean (Glycine max) plants to a MFSC of up to five different stresses (water deficit, salinity, low phosphate, acidity, and cadmium), in an increasing level of complexity, and conducted integrative transcriptomic-phenotypic analysis of their reproductive and vegetative tissues. We reveal that MFSC has a negative cumulative effect on soybean yield, that each set of MFSC condition elicits a unique transcriptomic response (that is different between flowers and leaves), and that selected genes expressed in leaves or flowers of soybean are linked to the effects of MFSC on different vegetative, physiological, and/or reproductive parameters. Our study identified networks and pathways associated with reactive oxygen species, ascorbic acid and aldarate, and iron/copper signaling/metabolism as promising targets for future biotechnological efforts to augment the resilience of reproductive tissues of major crop plants to MFSC. In addition, we provide unique phenotypic and transcriptomic datasets for dissecting the mechanistic effects of MFSC on the vegetative, physiological, and reproductive processes of a crop plant.

Assembly, comparative analysis, and utilization of a single haplotype reference genome for soybean.

Cultivar Williams 82 has served as the reference genome for the soybean research community since 2008, but is known to have areas of genomic heterogeneity among different sub-lines. This work provides an updated assembly (version Wm82.a6) derived from a specific sub-line known as Wm82-ISU-01 (seeds available under USDA accession PI 704477). The genome was assembled using Pacific BioSciences HiFi reads and integrated into chromosomes using HiC. The 20 soybean chromosomes assembled into a genome of 1.01Gb, consisting of 36 contigs. The genome annotation identified 48 387 gene models, named in accordance with previous assembly versions Wm82.a2 and Wm82.a4. Comparisons of Wm82.a6 with other near-gapless assemblies of Williams 82 reveal large regions of genomic heterogeneity, including regions of differential introgression from the cultivar Kingwa within approximately 30 Mb and 25 Mb segments on chromosomes 03 and 07, respectively. Additionally, our analysis revealed a previously unknown large (>20 Mb) heterogeneous region in the pericentromeric region of chromosome 12, where Wm82.a6 matches the 'Williams' haplotype while the other two near-gapless assemblies do not match the haplotype of either parent of Williams 82. In addition to the Wm82.a6 assembly, we also assembled the genome of 'Fiskeby III,' a rich resource for abiotic stress resistance genes. A genome comparison of Wm82.a6 with Fiskeby III revealed the nucleotide and structural polymorphisms between the two genomes within a QTL region for iron deficiency chlorosis resistance. The Wm82.a6 and Fiskeby III genomes described here will enhance comparative and functional genomics capacities and applications in the soybean community.

Effects of demographic history on recombination hotspots in soybean.

Recombination is the primary mechanism underlying genetic improvement in populations and allows plant breeders to create new allelic combinations for agronomic improvement. Soybean [Glycine max (L.) Merr.] has gone through multiple genetic bottlenecks that have significantly affected its genetic diversity, linkage disequilibrium, and altered allele frequencies. To investigate the impact of genetic bottlenecks on recombination hotspots in soybeans, historical recombination was studied in three soybean populations. The populations were wild soybean [Glycine soja (Sieb. and Zucc.)], landraces, and North American elite soybean cultivars that have been genotyped with the SoySNP50K BeadChip. While each population after a genetic bottleneck had an increased average haplotype block size, they did not have a significant difference in the number of hotspots between each population. Instead, the increase in observed haplotype block size is likely due to an elimination of individuals that contained historical recombination at hotspots which decreased the observed rate of recombination for the hotspot after each genetic bottleneck. Conversely, heterochromatic DNA which has an increased haplotype block size compared to euchromatic DNA had a significantly different number of hotspots but not a significant difference in the average hotspot recombination rate. Previously identified genomic motifs associated with hotspots were also associated with hotspots found in the historical populations suggesting a common mechanism. This characterization of historical recombination hotspots in soybeans provides further insights into the effect genetic bottlenecks and selection have on recombination hotspots.

Knockdown of ß-conglycinin α' and α subunits alters seed protein composition and improves salt tolerance in soybean.

Soybean is an important plant source of protein worldwide. Increasing demands for soybean can be met by improving the quality of its seed protein. In this study, GmCG-1, which encodes the ß-conglycinin α' subunit, was identified via combined genome-wide association study and transcriptome analysis. We subsequently knocked down GmCG-1 and its paralogues GmCG-2 and GmCG-3 with CRISPR-Cas9 technology and generated two stable multigene knockdown mutants. As a result, the ß-conglycinin content decreased, whereas the 11S/7S ratio, total protein content and sulfur-containing amino acid content significantly increased. Surprisingly, the globulin mutant exhibited salt tolerance in both the germination and seedling stages. Little is known about the relationship between seed protein composition and the salt stress response in soybean. Metabonomics and RNA-seq analysis indicated that compared with the WT, the mutant was formed through a pathway that was more similar to that of active salicylic acid biosynthesis; however, the synthesis of cytokinin exhibited greater defects, which could lead to increased expression of plant dehydrin-related salt tolerance proteins and cell membrane ion transporters. Population evolution analysis suggested that GmCG-1, GmCG-2, and GmCG-3 were selected during soybean domestication. The soybean accessions harboring GmCG-1Hap1 presented relatively high 11S/7S ratios and relatively high salt tolerance. In conclusion, knockdown of the ß-conglycinin α and α' subunits can improve the nutritional quality of soybean seeds and increase the salt tolerance of soybean plants, providing a strategy for designing soybean varieties with high nutritional value and high salt tolerance.

Multi-omics-driven advances in the understanding of triacylglycerol biosynthesis in oil seeds.

Vegetable oils are rich sources of polyunsaturated fatty acids and energy as well as valuable sources of human food, animal feed, and bioenergy. Triacylglycerols, which are comprised of three fatty acids attached to a glycerol backbone, are the main component of vegetable oils. Here, we review the development and application of multiple-level omics in major oilseeds and emphasize the progress in the analysis of the biological roles of key genes underlying seed oil content and quality in major oilseeds. Finally, we discuss future research directions in functional genomics research based on current omics and oil metabolic engineering strategies that aim to enhance seed oil content and quality, and specific fatty acids components according to either human health needs or industrial requirements.

Rhizobia cystathionine γ-lyase-derived H2S delays nodule senescence in soybean.

Hydrogen sulphide (H2S) is required for optimal establishment of soybean (Glycine max)-Sinorhizobium fredii symbiotic interaction, yet its role in regulating the nitrogen fixation-senescence transition remains poorly understood. A S. fredii cystathionine γ-lyase (CSE) mutant deficient in H2S synthesis showed early nodule senescence characterized by reduced nitrogenase activity, structural changes in nodule cells, and accelerated bacteroid death. In parallel, the CSE mutant facilitated the generation of reactive oxygen species (ROS) and elicited antioxidant responses. We observed that H2S-mediated persulfidation of cysteine C31/C80 in ascorbate peroxidase (APX) and C32 in APX2 modulated enzyme activity, thereby participating in hydrogen peroxide (H2O2) detoxification and delaying nodule senescence. Comparative transcriptomic analysis revealed a significant up-regulation of GmMYB128, an MYB transcription factor (TF), in the CSE mutant nodules. Functional analysis through overexpression and RNAi lines of GmMYB128 demonstrated its role as a positive regulator in nodule senescence. MYB128-OE inoculated with the CSE mutant strain exhibited a reduction in nitrogenase activity and a significant increase in DD15 expression, both of which were mitigated by NaHS addition. Changes at the protein level encompassed the activation of plant defenses alongside turnover in carbohydrates and amino acids. Our results suggest that H2S plays an important role in maintaining efficient symbiosis and preventing premature senescence of soybean nodules.

An NBS-LRR protein in the Rpp1 locus negates the dominance of Rpp1-mediated resistance against Phakopsora pachyrhizi in soybean.

The soybean Rpp1 locus confers resistance to Phakopsora pachyrhizi, causal agent of rust, and resistance is usually dominant over susceptibility. However, dominance of Rpp1-mediated resistance is lost when a resistant genotype (Rpp1 or Rpp1b) is crossed with susceptible line TMG06_0011, and the mechanism of this dominant susceptibility (DS) is unknown. Sequencing the Rpp1 region reveals that the TMG06_0011 Rpp1 locus has a single nucleotide-binding site leucine-rich repeat (NBS-LRR) gene (DS-R), whereas resistant PI 594760B (Rpp1b) is similar to PI 200492 (Rpp1) and has three NBS-LRR resistance gene candidates. Evidence that DS-R is the cause of DS was reflected in virus-induced gene silencing of DS-R in Rpp1b/DS-R or Rpp1/DS-R heterozygous plants with resistance partially restored. In heterozygous Rpp1b/DS-R plants, expression of Rpp1b candidate genes was not significantly altered, indicating no effect of DS-R on transcription. Physical interaction of the DS-R protein with candidate Rpp1b resistance proteins was supported by yeast two-hybrid studies and in silico modeling. Thus, we conclude that suppression of resistance most likely does not occur at the transcript level, but instead probably at the protein level, possibly with Rpp1 function inhibited by binding to the DS-R protein. The DS-R gene was found in other soybean lines, with an estimated allele frequency of 6% in a diverse population, and also found in wild soybean (Glycine soja). The identification of a dominant susceptible NBS-LRR gene provides insight into the behavior of NBS-LRR proteins and serves as a reminder to breeders that the dominance of an R gene can be influenced by a susceptibility allele.

Root growth of monocotyledons and dicotyledons is limited by different tissues.

Plant growth and morphogenesis are determined by the mechanical properties of its cell walls. Using atomic force microscopy, we have characterized the dynamics of cell wall elasticity in different tissues in developing roots of several plant species. The elongation growth zone of roots of all species studied was distinguished by a reduced modulus of elasticity of most cell walls compared to the meristem or late elongation zone. Within the individual developmental zones of roots, there were also significant differences in the elasticity of the cell walls of the different tissues, thus identifying the tissues that limit root growth in the different species. In cereals, this is mainly the inner cortex, whereas in dicotyledons this function is performed by the outer tissues-rhizodermis and cortex. These differences result in a different behaviour of the roots of these species during longitudinal dissection. Modelling of longitudinal root dissection using measured properties confirmed the difference shown. Thus, the morphogenesis of monocotyledonous and dicotyledonous roots relies on different tissues as growth limiting, which should be taken into account when analyzing the localization of associated molecular events. At the same time, no matrix polysaccharide was found whose immunolabelling in type I or type II cell walls would predict their mechanical properties. However, assessment of the degree of anisotropy of cortical microtubules showed a striking correlation with the elasticity of the corresponding cell walls in all species studied.

A chromosome-level genome assembly of the soybean pod borer: insights into larval transcriptional response to transgenic soybean expressing the pesticidal Cry1Ac protein.

BACKGROUND: Genetically modified (GM) crop plants with transgenic expression of Bacillus thuringiensis (Bt) pesticidal proteins are used to manage feeding damage by pest insects. The durability of this technology is threatened by the selection for resistance in pest populations. The molecular mechanism(s) involved in insect physiological response or evolution of resistance to Bt is not fully understood. RESULTS: To investigate the response of a susceptible target insect to Bt, the soybean pod borer, Leguminivora glycinivorella (Lepidoptera: Tortricidae), was exposed to soybean, Glycine max, expressing Cry1Ac pesticidal protein or the non-transgenic parental cultivar. Assessment of larval changes in gene expression was facilitated by a third-generation sequenced and scaffolded chromosome-level assembly of the L. glycinivorella genome (657.4 Mb; 27 autosomes + Z chromosome), and subsequent structural annotation of 18,197 RefSeq gene models encoding 23,735 putative mRNA transcripts. Exposure of L. glycinivorella larvae to transgenic Cry1Ac G. max resulted in prediction of significant differential gene expression for 204 gene models (64 up- and 140 down-regulated) and differential splicing among isoforms for 10 genes compared to unexposed cohorts. Differentially expressed genes (DEGs) included putative peritrophic membrane constituents, orthologs of Bt receptor-encoding genes previously linked or associated with Bt resistance, and those involved in stress responses. Putative functional Gene Ontology (GO) annotations assigned to DEGs were significantly enriched for 36 categories at GO level 2, respectively. Most significantly enriched cellular component (CC), biological process (BP), and molecular function (MF) categories corresponded to vacuolar and microbody, transport and metabolic processes, and binding and reductase activities. The DEGs in enriched GO categories were biased for those that were down-regulated (≥ 0.783), with only MF categories GTPase and iron binding activities were bias for up-regulation genes. CONCLUSIONS: This study provides insights into pathways and processes involved larval response to Bt intoxication, which may inform future unbiased investigations into mechanisms of resistance that show no evidence of alteration in midgut receptors.

Analysis of Ion Transport Properties of Glycine max HKT Transporters and Identifying a Regulation of GmHKT1;1 by the Non-Functional GmHKT1;4.

High-affinity potassium transporters (HKTs) play an important role in plants responding to salt stress, but the transport properties of the soybean HKT transporters at the molecular level are still unclear. Here, using Xenopus oocyte as a heterologous expression system and two-electrode voltage-clamp technique, we identified four HKT transporters, GmHKT1;1, GmHKT1;2, GmHKT1;3 and GmHKT1;4, all of which belong to type I subfamily, but have distinct ion transport properties. While GmHKT1;1, GmHKT1;2 and GmHKT1;3 function as Na+ transporters, GmHKT1;1 is less selective against K+ than the two other transporters. Astonishingly, GmHKT1;4, which lacks transmembrane segments and has no ion permeability, is significantly expressed, and its gene expression pattern is different from the other three GmHKTs under salt stress. Interestingly, GmHKT1;4 reduced the Na+/K+ currents mediated by GmHKT1;1. Further study showed that the transport ability of GmHKT1;1 regulated by GmHKT1;4 was related to the structural differences in the first intracellular domain and the fourth repeat domain. Overall, we have identified one unique GmHKT member, GmHKT1;4, which modulates the Na+ and K+ transport ability of GmHKT1;1 via direct interaction. Thus, we have revealed a new type of HKT interaction model for altering their ion transport properties.

Variation in shoot architecture traits and their relationship to canopy coverage and light interception in soybean (Glycine max).

BACKGROUND: In soybeans, faster canopy coverage (CC) is a highly desirable trait but a fully covered canopy is unfavorable to light interception at lower levels in the canopy with most of the incident radiation intercepted at the top of the canopy. Shoot architecture that influences CC is well studied in crops such as maize and wheat, and altering architectural traits has resulted in enhanced yield. However, in soybeans the study of shoot architecture has not been as extensive. RESULTS: This study revealed significant differences in CC among the selected soybean accessions. The rate of CC was found to decrease at the beginning of the reproductive stage (R1) followed by an increase during the R2-R3 stages. Most of the accessions in the study achieved maximum rate of CC between R2-R3 stages. We measured Light interception (LI), defined here as the ratio of Photosynthetically Active Radiation (PAR) transmitted through the canopy to the incoming PAR or the radiation above the canopy. LI was found to be significantly correlated with CC parameters, highlighting the relationship between canopy structure and light interception. The study also explored the impact of plant shape on LI and CO2 assimilation. Plant shape was characterized into distinct quantifiable parameters and by modeling the impact of plant shape on LI and CO2 assimilation, we found that plants with broad and flat shapes at the top maybe more photosynthetically efficient at low light levels, while conical shapes were likely more advantageous when light was abundant. Shoot architecture of plants in this study was described in terms of whole plant, branching and leaf-related traits. There was significant variation for the shoot architecture traits between different accessions, displaying high reliability. We found that that several shoot architecture traits such as plant height, and leaf and internode-related traits strongly influenced CC and LI. CONCLUSION: In conclusion, this study provides insight into the relationship between soybean shoot architecture, canopy coverage, and light interception. It demonstrates that novel shoot architecture traits we have defined here are genetically variable, impact CC and LI and contribute to our understanding of soybean morphology. Correlations between different architecture traits, CC and LI suggest that it is possible to optimize soybean growth without compromising on light transmission within the soybean canopy. In addition, the study underscores the utility of integrating low-cost 2D phenotyping as a practical and cost-effective alternative to more time-intensive 3D or high-tech low-throughput methods. This approach offers a feasible means of studying basic shoot architecture traits at the field level, facilitating a broader and efficient assessment of plant morphology.

The application potential of mepiquat chloride in soybean: improvement of yield characteristics and drought resistance.

BACKGROUND: Drought can result in yield losses, the application of plant growth regulators is an effective measure to improve drought resistance and yield. The objective of the study was to explore the application potential of mepiquat chloride (MC) in regulating soybean yield and drought resistance. METHODS: In this study, a three-year field experiment was designed and combined with drought experiments to measure the yield of popularized varieties during 2021-2022 and drought-resistant and drought-sensitive varieties were selected, and planted in the field in 2023. RESULTS: MC increased the yield of HN84 and HN87 for two consecutive years from 2021 to 2022 and improved their physiological characteristics under field conditions. Under M200 treatment, the yield of HN84 increased by 6.93% and 9.46%, and HN87 increased by 11.11% and 15.72%. Different concentrations of MC have different effects on soybeans. The maximum increase of SOD, POD and proline in HN84 under M400 treatment reached 71.92%, 63.26% and 71.54%, respectively; the maximum increase of SOD, POD and proline in HN87 under M200 treatment reached 21.96%, 93.49% and 40.45%, respectively. In 2023, the foliar application of MC improved the physiological characteristics of HN44 and HN65 under drought-stress conditions. On the eighth day of drought treatment, compared to the drought treatment, the leaf and root dry weight of HN44 under M100 treatment increased by 17.91% and 32.76%, respectively; the dry weight of leaves and roots of HN65 increased by 20.74% and 29.29% under M200 treatment, respectively. MC also reduced malondialdehyde (MDA) content, decreased antioxidant enzyme activity and proline content. In addition, different concentrations of MC increased the chlorophyll fluorescence parameters (Fs, Fv/Fm, YII, and SPAD). In the field, the plant height of the two varieties decreased significantly, the yield increased, the number of two-grain and three-grain pods increased, and the stem length at the bottom and middle decreased with MC induction. CONCLUSIONS: The application of 100-200 mg/L MC effectively improved drought resistance and increased yield. This study provided support for the rational application of MC in soybean production.

Lipid mediated plant immunity in susceptible and tolerant soybean cultivars in response to Phytophthora sojae colonization and infection.

BACKGROUND: Soybean is one of the most cultivated crops globally and a staple food for much of the world's population. The annual global crop losses due to infection by Phytophthora sojae is currently estimated at $20B USD, yet we have limited understanding of the role of lipid mediators in the adaptative strategies used by the host plant to limit infection. Since root is the initial site of this infection, we examined the infection process in soybean root infected with Phytophthora sojae using scanning electron microscopy to observe the changes in root morphology and a multi-modal lipidomics approach to investigate how soybean cultivars remodel their lipid mediators to successfully limit infection by Phytophthora sojae. RESULTS: The results reveal the presence of elevated biogenic crystals and more severe damaged cells in the root morphology of the infected susceptible cultivar compared to the infected tolerant cultivars. Furthermore, induced accumulation of stigmasterol was observed in the susceptible cultivar whereas, induced accumulation of phospholipids and glycerolipids occurred in tolerant cultivar. CONCLUSION: The altered lipidome reported in this study suggest diacylglycerol and phosphatidic acid mediated lipid signalling impacting phytosterol anabolism appears to be a strategy used by tolerant soybean cultivars to successfully limit infection and colonization by Phytophthora sojae.

Soybean steroids improve crop abiotic stress tolerance and increase yield.

Sterols have long been associated with diverse fields, such as cancer treatment, drug development, and plant growth; however, their underlying mechanisms and functions remain enigmatic. Here, we unveil a critical role played by a GmNF-YC9-mediated CCAAT-box transcription complex in modulating the steroid metabolism pathway within soybeans. Specifically, this complex directly activates squalene monooxygenase (GmSQE1), which is a rate-limiting enzyme in steroid synthesis. Our findings demonstrate that overexpression of either GmNF-YC9 or GmSQE1 significantly enhances soybean stress tolerance, while the inhibition of SQE weakens this tolerance. Field experiments conducted over two seasons further reveal increased yields per plant in both GmNF-YC9 and GmSQE1 overexpressing plants under drought stress conditions. This enhanced stress tolerance is attributed to the reduction of abiotic stress-induced cell oxidative damage. Transcriptome and metabolome analyses shed light on the upregulation of multiple sterol compounds, including fucosterol and soyasaponin II, in GmNF-YC9 and GmSQE1 overexpressing soybean plants under stress conditions. Intriguingly, the application of soybean steroids, including fucosterol and soyasaponin II, significantly improves drought tolerance in soybean, wheat, foxtail millet, and maize. These findings underscore the pivotal role of soybean steroids in countering oxidative stress in plants and offer a new research strategy for enhancing crop stress tolerance and quality from gene regulation to chemical intervention.

Complete genome sequence of Bacillus velezensis strain Ag109, a biocontrol agent against plant-parasitic nematodes and Sclerotinia sclerotiorum.

Soybean is the main oilseed cultivated worldwide. Even though Brazil is the world's largest producer and exporter of soybean, its production is severely limited by biotic factors. Soil borne diseases are the most damaging biotic stressors since they significantly reduce yield and are challenging to manage. In this context, the present study aimed to evaluate the potential of a bacterial strain (Ag109) as a biocontrol agent for different soil pathogens (nematodes and fungi) of soybean. In addition, the genome of Ag109 was wholly sequenced and genes related to secondary metabolite production and plant growth promotion were mined. Ag109 showed nematode control in soybean and controlled 69 and 45% of the populations of Meloidogyne javanica and Pratylenchus brachyurus, respectively. Regarding antifungal activity, these strains showed activity against Macrophomia phaseolina, Rhizoctonia solani, and Sclerotinia sclerotiorum. For S. sclerotiorum, this strain increased the number of healthy plants and root dry mass compared to the control (with inoculation). Based on the average nucleotide identity and digital DNA-DNA hybridization, this strain was identified as Bacillus velezensis. Diverse clusters of specific genes related to secondary metabolite biosynthesis and root growth promotion were identified, highlighting the potential of this strain to be used as a multifunctional microbial inoculant that acts as a biological control agent while promoting plant growth in soybean.

The soybean sugar transporter GmSWEET6 participates in sucrose transport towards fungi during arbuscular mycorrhizal symbiosis.

In arbuscular mycorrhizal (AM) symbiosis, sugars in root cortical cells could be exported as glucose or sucrose into peri-arbuscular space for use by AM fungi. However, no sugar transporter has been identified to be involved in sucrose export. An AM-inducible SWEET transporter, GmSWEET6, was functionally characterised in soybean, and its role in AM symbiosis was investigated via transgenic plants. The expression of GmSWEET6 was enhanced by inoculation with the cooperative fungal strain in both leaves and roots. Heterologous expression in a yeast mutant showed that GmSWEET6 mainly transported sucrose. Transgenic plants overexpressing GmSWEET6 increased sucrose concentration in root exudates. Overexpression or knockdown of GmSWEET6 decreased plant dry weight, P content, and sugar concentrations in non-mycorrhizal plants, which were partly recovered in mycorrhizal plants. Intriguingly, overexpression of GmSWEET6 increased root P content and decreased the percentage of degraded arbuscules, while knockdown of GmSWEET6 increased root sugar concentrations in RNAi2 plants and the percentage of degraded arbuscules in RNAi1 plants compared with wild-type plants when inoculated with AM fungi. These results in combination with subcellular localisation of GmSWEET6 to peri-arbuscular membranes strongly suggest that GmSWEET6 is required for AM symbiosis by mediating sucrose efflux towards fungi.

Dynamic Landscapes of Long Noncoding RNAs During Early Root Development and Differentiation in Glycine max and Glycine soja.

Soybean (Glycine max) is an important crop for its nutritional value. Its wild relative, Glycine soja, provides a valuable genetic resource for improving soybean productivity. Root development and differentiation are essential for soybean plants to take up water and nutrients, store energy and anchor themselves. Long noncoding RNAs (lncRNAs) have been reported to play critical roles in various biological processes. However, the spatiotemporal landscape of lncRNAs during early root development and differentiation in soybeans is scarcely characterized. Using RNA sequencing and transcriptome assembly, we identified 1578 lncRNAs in G. max and 1454 in G. soja, spanning various root portions and time points. Differential expression analysis revealed 82 and 69 lncRNAs exhibiting spatiotemporally differential expression patterns in G. max and G. soja, respectively, indicating their involvement in the early stage of root architecture formation. By elucidating multiple competitive endogenous RNA (ceRNA) networks involving lncRNAs, microRNAs and protein-coding RNAs, we unveiled intricate regulatory mechanisms of lncRNA in early root development and differentiation. Our efforts significantly expand the transcriptome annotations of soybeans, unravel the dynamic landscapes of lncRNAs during early root development and differentiation, and provide valuable resources into the field of soybean root research.

Contextualized Metabolic Modelling Revealed Factors Affecting Isoflavone Accumulation in Soybean Seeds.

Isoflavones, secondary metabolites with numerous health benefits, are predominantly found in legume seeds, especially soybean; however, their contents in domesticated soybean seeds are highly variable. Wild soybeans are known for higher seed isoflavone contents than cultivars. Here we used experimental and modelling approaches on wild soybean (W05) and cultivated soybean (C08) to delineate factors influencing isoflavone accumulation. We found imported nutrients were converted into storage compounds, with isoflavone accumulation in W05 seeds being faster than in C08 ones. The isoflavone accumulation during seed development was simulated using context-specific cotyledon metabolic models of four developmental stages on cultivar C08, and the metabolic burden imposed by increasing biomass was evaluated. Trade-off analyses between biomass and isoflavone suggest that high biomass requirement in cultivars could limit the reallocation of resources for secondary metabolite production. Isoflavone production in mature seeds was also influenced by biomass compositions. Seeds with higher carbohydrate contents favour isoflavone production, while those with highest protein and oil contents had lowest isoflavone contents. Although seeds could synthesize isoflavones on their own, the predicted fluxes from biosynthesis alone were lower than the empirical levels. Shadow price analyses indicated that isoflavone accumulation depended on both intrinsic biosynthesis and direct contribution from the plant.