Soybean is one of the most important crops, providing large amounts of dietary proteins and edible oil, and is also an excellent model for studying evolution of duplicated genes. However, relative to the model plants Arabidopsis and rice, the present knowledge about soybean transcriptome is quite limited. In this study, we employed RNA-seq to investigate transcriptomes of 11 soybean tissues, for genome-wide discovery of truly expressed genes, and novel and alternative transcripts, as well as analyses of conservation and divergence of duplicated genes and their functional implications. We detected a total of 54,132 high-confidence expressed genes, and identified 6,718 novel transcriptional regions with a mean length of 372 bp. We also provided strong evidence for alternative splicing (AS) events for ~15.9% of the genes with two or more exons. Among them, 1,834 genes exhibited stage-dependent AS, and 202 genes had tissue-biased exon-skipping events. We further defined the conservation and divergence in expression patterns between duplicated gene pairs from recent whole genome duplications (WGDs); differentially expressed genes, tissue preferentially expressed genes, transcription factors and specific gene family members were identified for shoot apical meristem and flower development. Our results significantly improved soybean gene annotation, and also provide valuable resources for functional genomics and studies of the evolution of duplicated genes from WGDs in soybean.

One of the overarching goals of soybean breeding is to develop lines that combine increased yield with improved quality characteristics. High-density-marker QTL mapping can serve as an effective strategy to identify novel genomic information to facilitate crop improvement. In this study, we genotyped a recombinant inbred line (RIL) population (Zhonghuang 24 × Huaxia 3) using a restriction-site associated DNA sequencing (RAD-seq) approach. A high-density soybean genetic map was constructed and used to identify several QTLs that were shown to influence six yield-related and two quality traits. A total of 47,472 single-nucleotide polymorphisms (SNPs) were detected for the RILs that were integrated into 2639 recombination bin units, with an average distance of 1.00 cM between adjacent markers. Forty seven QTLs for yield-related traits and 13 QTLs for grain quality traits were found to be distributed on 16 chromosomes in the 2 year studies. Among them, 18 QTLs were stable, and were identified in both analyses. Twenty six QTLs were identified for the first time, with a single QTL (qNN19a) in a 56 kb region explaining 32.56% of phenotypic variation, and an additional 10 of these were novel, stable QTLs. Moreover, 8 QTL hotpots on four different chromosomes were identified for the correlated traits. With RAD-sequencing, some novel QTLs and important QTL clusters for both yield-related and quality traits were identified based on a new, high-density bin linkage map. Three predicted genes were selected as candidates that likely have a direct or indirect influence on both yield and quality in soybean. Our findings will be helpful for understanding common genetic control mechanisms of co-localized traits and to select cultivars for further analysis to predictably modulate soybean yield and quality simultaneously.

Multidrug and toxic compound extrusion (MATE) transporters, which exist widely in plants, function as crucial regulators in plant resistance to aluminum (Al) toxicity by inducing citrate efflux. However, the functions of most MATE family members in soybean (Glycine soja) remain to be elucidated. Expression pattern analysis showed that GsMATE was constitutively expressed in different soybean organs, with the highest level in root compared with those in stem, leaf and cotyledon. In addition, Al stress induced expression of GsMATE in soybean. Temporal analysis indicated that GsMATE expression was greatly enhanced by increasing concentrations of aluminum [Al3+] after short exposure, reaching the high levels detected in the BW69 (Al-resistant) and the JW81 (Al-sensitive) lines of Glycine soja of wild soybean at 6 h and 8 h, respectively. Furthermore, transient GsMATE expression in Arabidopsis protoplasts showed that GsMATE protein localized to the plasma membrane. Overexpression of GsMATE on an Arabidopsis columbia-0 (Col-0) background resulted in increased Al tolerance in transgenic plants. Analysis of hematoxylin staining showed that the roots of GsMATE transgenic lines were stained less intensely than those of the wild-type exposured to the same AlCl3 concentrations. Therefore, GsMATE enhanced the resistance of transgenic plants to Al toxicity by reducing Al accumulation in Arabidopsis roots. In summary, our results indicate that GsMATE is responsive to aluminum stress and may participate in the regulation of sensitivity to Al toxicity in Arabidopsis. In addition, the GsMATE protein is an Al-induced citrate transporter of the MATE family and exerts an essential role in Al tolerance in Glycine soja.

Acyl activating enzyme 3 (AAE3) was identified as being involved in the acetylation pathway of oxalate degradation, which regulates the responses to biotic and abiotic stresses in various higher plants. Here, we investigated the role of Glycine sojaAAE3 (GsAAE3) in Cadmium (Cd) and Aluminum (Al) tolerances. The recombinant GsAAE3 protein showed high activity toward oxalate, with a K_m of 105.10 ± 12.30 μM and V_max of 12.64 ± 0.34 μmol min⁻¹ mg⁻¹ protein, suggesting that it functions as an oxalyl–CoA synthetase. The expression of a GsAAE3–green fluorescent protein (GFP) fusion protein in tobacco leaves did not reveal a specific subcellular localization pattern of GsAAE3. An analysis of the GsAAE3 expression pattern revealed an increase in GsAAE3 expression in response to Cd and Al stresses, and it is mainly expressed in root tips. Furthermore, oxalate accumulation induced by Cd and Al contributes to the inhibition of root growth in wild soybean. Importantly, GsAAE3 overexpression increases Cd and Al tolerances in A. thaliana and soybean hairy roots, which is associated with a decrease in oxalate accumulation. Taken together, our data provide evidence that the GsAAE3-encoded protein plays an important role in coping with Cd and Al stresses.

Salt stress has a detrimental impact on crop yield, quality, and profitability. The tau-like glutathione transferases (GSTs) represent a significant group of enzymes that play a crucial role in plant stress responses, including salt stress. In this study, we identified a tau-like glutathione transferase family gene from soybean named GmGSTU23. Expression pattern analysis revealed that GmGSTU23 was predominantly expressed in the roots and flowers and exhibited a concentration–time-specific pattern in response to salt stress. Transgenic lines were generated and subjected to phenotypic characterization under salt stress. The transgenic lines exhibited increased salt tolerance, root length, and fresh weight compared to the wild type. Antioxidant enzyme activity and malondialdehyde content were subsequently measured, and the data revealed no significant differences between the transgenic and wild-type plants in the absence of salt stress. However, under salt stress, the wild-type plants exhibited significantly lower activities of SOD, POD, and CAT than the three transgenic lines, whereas the activity of APX and the content of MDA showed the opposite trend. We identified changes in glutathione pools and associated enzyme activity to gain insights into the underlying mechanisms of the observed phenotypic differences. Notably, under salt stress, the transgenic Arabidopsis’s GST activity, GR activity, and GSH content were significantly higher than those of the wild type. In summary, our findings suggest that GmGSTU23 mediates the scavenging of reactive oxygen species and glutathione by enhancing the activity of glutathione transferase, thereby conferring enhanced tolerance to salt stress in plants.