Natural variation in Fatty Acid 9 is a determinant of fatty acid and protein content.

Soybean is one of the most economically important crops worldwide and an important source of unsaturated fatty acids and protein for the human diet. Consumer demand for healthy fats and oils is increasing, and the global demand for vegetable oil is expected to double by 2050. Identification of key genes that regulate seed fatty acid content can facilitate molecular breeding of high-quality soybean varieties with enhanced fatty acid profiles. Here, we analysed the genetic architecture underlying variations in soybean seed fatty acid content using 547 accessions, including mainly landraces and cultivars from northeastern China. Through fatty acid profiling, genome re-sequencing, population genomics analyses, and GWAS, we identified a SEIPIN homologue at the FA9 locus as an important contributor to seed fatty acid content. Transgenic and multiomics analyses confirmed that FA9 was a key regulator of seed fatty acid content with pleiotropic effects on seed protein and seed size. We identified two major FA9 haplotypes in 1295 resequenced soybean accessions and assessed their phenotypic effects in a field planting of 424 accessions. Soybean accessions carrying FA9H2 had significantly higher total fatty acid contents and lower protein contents than those carrying FA9H1 . FA9H2 was absent in wild soybeans but present in 13% of landraces and 26% of cultivars, suggesting that it may have been selected during soybean post-domestication improvement. FA9 therefore represents a useful genetic resource for molecular breeding of high-quality soybean varieties with specific seed storage profiles.

Multi-Omics Analysis of a Chromosome Segment Substitution Line Reveals a New Regulation Network for Soybean Seed Storage Profile.

Soybean, a major source of oil and protein, has seen an annual increase in consumption when used in soybean-derived products and the broadening of its cultivation range. The demand for soybean necessitates a better understanding of the regulatory networks driving storage protein accumulation and oil biosynthesis to broaden its positive impact on human health. In this study, we selected a chromosome segment substitution line (CSSL) with high protein and low oil contents to investigate the underlying effect of donor introgression on seed storage through multi-omics analysis. In total, 1479 differentially expressed genes (DEGs), 82 differentially expressed proteins (DEPs), and 34 differentially expressed metabolites (DEMs) were identified in the CSSL compared to the recurrent parent. Based on Gene Ontology (GO) term analysis and the Kyoto Encyclopedia of Genes and Genomes enrichment (KEGG), integrated analysis indicated that 31 DEGs, 24 DEPs, and 13 DEMs were related to seed storage functionality. Integrated analysis further showed a significant decrease in the contents of the seed storage lipids LysoPG 16:0 and LysoPC 18:4 as well as an increase in the contents of organic acids such as L-malic acid. Taken together, these results offer new insights into the molecular mechanisms of seed storage and provide guidance for the molecular breeding of new favorable soybean varieties.

The impact of GmTSA and GmALS on soybean salt tolerance: uncovering the molecular landscape of amino acid and secondary metabolism pathways.

KEY MESSAGE: GmTSA and GmALS were screened out for salt stress in soybean and explore the poteintial amino acid secondary metabolism pathways. Soybean (Glycine max L.) is an oil and protein crop of global importance, and salinity has significant effects on soybean growth. Here, a population of soybean chromosome segment substitution lines was screened to identify highly salt-tolerant lines. In total, 24 quantitative trait loci (QTLs) on seven chromosomes were associated with salt tolerance, and CSSL_R71 was selected for further analysis. Although numerous genes were differentially expressed in CSSL_R71 in response to salt statically no differently, transcript levels of classical salt-response genes, including those of the salt overly sensitive pathway. Rather, salt tolerance in CSSL_R71 was associated with changes in amino acid and lipid metabolism. In particular, changes in p-coumaric acid, shikimic acid, and pyrrole-2-carboxylic acid levels accompanied salt tolerance in CSSL_R71. Eleven differentially expressed genes (DEGs) related to amino acid and secondary metabolism were identified as candidate genes on the substituted chromosome fragment. Six of these showed differences in coding sequence between the parental genotypes. Crucially, overexpression of GmTSA (Glyma.03G158400, tryptophan synthase) significantly enhanced salt tolerance in soybean hairy roots, whereas overexpression of GmALS (Glyma.13G241000, acetolactate synthase) decreased salt tolerance. Two KASP markers were developed for GmALS and used to genotype salt-tolerant and salt-sensitive lines in the CSSL population. Non-synonymous mutations were directly associated with salt tolerance. Taken together, these data provide evidence that changes in amino acid and secondary metabolism have the potential to confer salt tolerance in soybean.

Multi-Omics Techniques for Soybean Molecular Breeding.

Soybean is a major crop that provides essential protein and oil for food and feed. Since its origin in China over 5000 years ago, soybean has spread throughout the world, becoming the second most important vegetable oil crop and the primary source of plant protein for global consumption. From early domestication and artificial selection through hybridization and ultimately molecular breeding, the history of soybean breeding parallels major advances in plant science throughout the centuries. Now, rapid progress in plant omics is ushering in a new era of precision design breeding, exemplified by the engineering of elite soybean varieties with specific oil compositions to meet various end-use targets. The assembly of soybean reference genomes, made possible by the development of genome sequencing technology and bioinformatics over the past 20 years, was a great step forward in soybean research. It facilitated advances in soybean transcriptomics, proteomics, metabolomics, and phenomics, all of which paved the way for an integrated approach to molecular breeding in soybean. In this review, we summarize the latest progress in omics research, highlight novel findings made possible by omics techniques, note current drawbacks and areas for further research, and suggest that an efficient multi-omics approach may accelerate soybean breeding in the future. This review will be of interest not only to soybean breeders but also to researchers interested in the use of cutting-edge omics technologies for crop research and improvement.

SoySNP618K array: A high-resolution single nucleotide polymorphism platform as a valuable genomic resource for soybean genetics and breeding.

Innovations in genomics have enabled the development of low-cost, high-resolution, single nucleotide polymorphism (SNP) genotyping arrays that accelerate breeding progress and support basic research in crop science. Here, we developed and validated the SoySNP618K array (618,888 SNPs) for the important crop soybean. The SNPs were selected from whole-genome resequencing data containing 2,214 diverse soybean accessions; 29.34% of the SNPs mapped to genic regions representing 86.85% of the 56,044 annotated high-confidence genes. Identity-by-state analyses of 318 soybeans revealed 17 redundant accessions, highlighting the potential of the SoySNP618K array in supporting gene bank management. The patterns of population stratification and genomic regions enriched through domestication were highly consistent with previous findings based on resequencing data, suggesting that the ascertainment bias in the SoySNP618K array was largely compensated for. Genome-wide association mapping in combination with reported quantitative trait loci enabled fine-mapping of genes known to influence flowering time, E2 and GmPRR3b, and of a new candidate gene, GmVIP5. Moreover, genomic prediction of flowering and maturity time in 502 recombinant inbred lines was highly accurate (>0.65). Thus, the SoySNP618K array is a valuable genomic tool that can be used to address many questions in applied breeding, germplasm management, and basic crop research.

Fine mapping of a QTL locus (QNFSP07-1) and analysis of candidate genes for four-seeded pods in soybean.

Soybean [Glycine max (L.) Merr.] is an important grain and oil crop in the world, and it is the main source of high-quality protein. The number of four-seeded pods is a quantitative trait in soybean and is closely related to yield in terms of breeding. Therefore, it is of great significance to study the inheritance of four-seed pods and to excavate related genes for improving soybean yield. In this study, individuals with high ratio of four-seed pods which from chromosome segment substitution lines (CSSLs) that can be stably inherited were selected as the parent, and Suinong 14 (SN14) was used as recurrent parent to construct secondary mapping population via marker-assisted selection. From 2006 to 2017, QTL analysis was performed using secondary mapping populations, and the initial QTL mapping interval was 0.67 Mb and was located on Gm07. Based on the initial QTL mapping results, individuals that were heterozygous at the interval (36,116,118-37,399,738 bp) were screened in 2018, and the heterozygous individuals were subjected to inbreeding to obtain 13 F3 populations, with a target interval of 321 kb. Gene annotation was performed on the fine mapping interval, and 27 genes were obtained. Among 27 genes, Glyma.07G200900 and Glyma.07G201200 were identified as candidate genes. qRT-PCR was used to measure the expression of the 2 candidate genes at different developmental stages of soybean, and the expression levels of the 2 candidate genes in terms of cell division (axillary buds, COTs, EMs) were higher than those in terms of cell expansion (MM, LM), and these genes play a positive regulatory role in the formation of four-seeded pods. Haplotype analysis of 2 candidate genes which shows that Glyma.07G201200 has two excellent haplotypes, and the significance level between the two excellent haplotypes at p < 0.05. Those results provide the information for gene map-based cloning and molecular marker-assisted breeding of the number of four-seeded pod in soybean. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-021-01265-6.

Enhanced production of seed oil with improved fatty acid composition by overexpressing NAD+ -dependent glycerol-3-phosphate dehydrogenase in soybean.

There is growing interest in expanding the production of soybean oils (mainly triacylglycerol, or TAG) to meet rising feed demand and address global energy concerns. We report that a plastid-localized glycerol-3-phosphate dehydrogenase (GPDH), encoded by GmGPDHp1 gene, catalyzes the formation of glycerol-3-phosphate (G3P), an obligate substrate required for TAG biosynthesis. Overexpression of GmGPDHp1 increases soybean seed oil content with high levels of unsaturated fatty acids (FAs), especially oleic acid (C18:1), without detectably affecting growth or seed protein content or seed weight. Based on the lipidomic analyses, we found that the increase in G3P content led to an elevated diacylglycerol (DAG) pool, in which the Kennedy pathway-derived DAG was mostly increased, followed by PC-derived DAG, thereby promoting the synthesis of TAG containing relatively high proportion of C18:1. The increased G3P levels induced several transcriptional alterations of genes involved in the glycerolipid pathways. In particular, genes encoding the enzymes responsible for de novo glycerolipid synthesis were largely upregulated in the transgenic lines, in-line with the identified biochemical phenotype. These results reveal a key role for GmGPDHp1-mediated G3P metabolism in enhancing TAG synthesis and demonstrate a strategy to modify the FA compositions of soybean oils for improved nutrition and biofuel.

RNA Sequencing-Associated Study Identifies GmDRR1 as Positively Regulating the Establishment of Symbiosis in Soybean.

In soybean (Glycine max)-rhizobium interactions, the type III secretion system (T3SS) of rhizobium plays a key role in regulating host specificity. However, the lack of information on the role of T3SS in signaling networks limits our understanding of symbiosis. Here, we conducted an RNA sequencing analysis of three soybean chromosome segment substituted lines, one female parent and two derived lines with different chromosome-substituted segments of wild soybean and opposite nodulation patterns. By analyzing chromosome-linked differentially expressed genes in the substituted segments and quantitative trait loci (QTL)-assisted selection in the substituted-segment region, genes that may respond to type III effectors to mediate plant immunity-related signaling were identified. To narrow down the number of candidate genes, QTL assistant was used to identify the candidate region consistent with the substituted segments. Furthermore, one candidate gene, GmDRR1, was identified in the substituted segment. To investigate the role of GmDRR1 in symbiosis establishment, GmDRR1-overexpression and RNA interference soybean lines were constructed. The nodule number increased in the former compared with wild-type soybean. Additionally, the T3SS-regulated effectors appeared to interact with the GmDDR1 signaling pathway. This finding will allow the detection of T3SS-regulated effectors involved in legume-rhizobium interactions.

QTL analysis of nodule traits and the identification of loci interacting with the type III secretion system in soybean.

Symbiotic nitrogen fixation is the main source of nitrogen for soybean growth. Since the genotypes of rhizobia and soybean germplasms vary, the nitrogen-fixing ability of soybean after inoculation also varies. A few studies have reported that quantitative trait loci (QTLs) control biological nitrogen fixation traits, even soybean which is an important crop. The present study reported that the Sinorhizobium fredii HH103 gene rhcJ belongs to the tts (type III secretion) cluster and that the mutant HH103ΩrhcJ can clearly decrease the number of nodules in American soybeans. However, few QTLs of nodule traits have been identified. This study used a soybean (Glycine max (L.) Merr.) 'Charleston' × 'Dongnong 594' (C × D, n = 150) recombinant inbred line (RIL). Nodule traits were analysed in the RIL population after inoculation with S. fredii HH103 and the mutant HH103ΩrhcJ. Plants were grown in a greenhouse with a 16-h light cycle at 26 °C and an 8-h dark cycle at 18 °C. Then, 4 weeks after inoculation, plants were harvested for evaluation of nodule traits. Through QTL mapping, 16 QTLs were detected on 8 chromosomes. Quantitative PCR (qRT-PCR) and RNA-seq analysis determined that the genes Glyma.04g060600, Glyma.18g159800 and Glyma.13g252600 might interact with rhcJ.

Meta-analysis and transcriptome profiling reveal hub genes for soybean seed storage composition during seed development.

Soybean is an important crop providing edible oil and protein source. Soybean oil and protein contents are quantitatively inherited and significantly affected by environmental factors. In this study, meta-analysis was conducted based on soybean physical maps to integrate quantitative trait loci (QTLs) from multiple experiments in different environments. Meta-QTLs for seed oil, fatty acid composition, and protein were identified. Of them, 11 meta-QTLs were located on hot regions for both seed oil and protein. Next, we selected 4 chromosome segment substitution lines with different seed oil and protein contents to characterize their 3 years of phenotype selection in the field. Using strand-specific RNA-sequencing analysis, we profile the time-course transcriptome patterns of soybean seeds at early maturity, middle maturity, and dry seed stages. Pairwise comparison and K-means clustering analysis revealed 7,482 differentially expressed genes and 45 expression patterns clusters. Weighted gene coexpression network analysis uncovered 46 modules of gene expression patterns. The 2 most significant coexpression networks were visualized, and 7 hub genes were identified that were involved in soybean oil and seed storage protein accumulation processes. Our results provided a transcriptome dataset for soybean seed development, and the candidate hub genes represent a foundation for further research.

Identification of Soybean Genes Whose Expression is Affected by the Ensifer fredii HH103 Effector Protein NopP.

In some legume⁻rhizobium symbioses, host specificity is influenced by rhizobial nodulation outer proteins (Nops). However, the genes encoding host proteins that interact with Nops remain unknown. We generated an Ensifer fredii HH103 NopP mutant (HH103ΩNopP), and analyzed the nodule number (NN) and nodule dry weight (NDW) of 10 soybean germplasms inoculated with the wild-type E. fredii HH103 or the mutant strain. An analysis of recombinant inbred lines (RILs) revealed the quantitative trait loci (QTLs) associated with NopP interactions. A soybean genomic region containing two overlapping QTLs was analyzed in greater detail. A transcriptome analysis and qRT-PCR assay were used to identify candidate genes encoding proteins that interact with NopP. In some germplasms, NopP positively and negatively affected the NN and NDW, while NopP had different effects on NN and NDW in other germplasms. The QTL region in chromosome 12 was further analyzed. The expression patterns of candidate genes Glyma.12g031200 and Glyma.12g073000 were determined by qRT-PCR, and were confirmed to be influenced by NopP.

Oceanobacillus picturae alleviates cadmium stress and promotes growth in soybean seedlings.

Cadmium (Cd) is a heavy metal that significantly impacts human health and the environment. Microorganisms play a crucial role in reducing heavy metal stress in plants; however, the mechanisms by which microorganisms enhance plant tolerance to Cd stress and the interplay between plants and microorganisms under such stress remain unclear. In this study, Oceanobacillus picturae (O. picturae) was isolated for interaction with soybean seedlings under Cd stress. Results indicated that Cd treatment alone markedly inhibited soybean seedling growth. Conversely, inoculation with O. picturae significantly improved growth indices such as plant height, root length, and fresh weight, while also promoting recovery in soil physiological indicators and pH. Metabolomic and transcriptomic analyses identified 157 genes related to aspartic acid, cysteine, and flavonoid biosynthesis pathways. Sixty-three microbial species were significantly associated with metabolites in these pathways, including pathogenic, adversity-resistant, and bioconductive bacteria. This research experimentally demonstrates, for the first time, the growth-promoting effect of the O. picturae strain on soybean seedlings under non-stress conditions. It also highlights its role in enhancing root growth and reducing Cd accumulation in the roots under Cd stress. Additionally, through the utilization of untargeted metabolomics, metagenomics, and transcriptomics for a multi-omics analysis, we investigated the impact of O. picturae on the soil microbiome and its correlation with differential gene expression in plants. This innovative approach unveils the molecular mechanisms underlying O. picturae's promotion of root growth and adaptation to Cd stress.

Integrated Transcriptomic and Proteomic Characterization of a Chromosome Segment Substitution Line Reveals the Regulatory Mechanism Controlling the Seed Weight in Soybean.

Soybean is the major global source of edible oils and vegetable proteins. Seed size and weight are crucial traits determining the soybean yield. Understanding the molecular regulatory mechanism underlying the seed weight and size is helpful for improving soybean genetic breeding. The molecular regulatory pathways controlling the seed weight and size were investigated in this study. The 100-seed weight, seed length, seed width, and seed weight per plant of a chromosome segment substitution line (CSSL) R217 increased compared with those of its recurrent parent 'Suinong14' (SN14). Transcriptomic and proteomic analyses of R217 and SN14 were performed at the seed developmental stages S15 and S20. In total, 2643 differentially expressed genes (DEGs) and 208 differentially accumulated proteins (DAPs) were detected at S15, and 1943 DEGs and 1248 DAPs were detected at S20. Furthermore, integrated transcriptomic and proteomic analyses revealed that mitogen-activated protein kinase signaling and cell wall biosynthesis and modification were potential pathways associated with seed weight and size control. Finally, 59 candidate genes that might control seed weight and size were identified. Among them, 25 genes were located on the substituted segments of R217. Two critical pathways controlling seed weight were uncovered in our work. These findings provided new insights into the seed weight-related regulatory network in soybean.

Wild rodents seed choice is relevant for sustainable agriculture.

Mitigating pre-harvest sprouting (PHS) and post-harvest food loss (PHFL) is essential for enhancing food securrity. To reduce food loss, the use of plant derived specialized metabolites can represent a good approach to develop a more eco-friendly agriculture. Here, we have discovered that soybean seeds hidden underground during winter by Tscherskia triton and Apodemus agrarius during winter possess a higher concentration of volatile organic compounds (VOCs) compared to those remaining exposed in fields. This selection by rodents suggests that among the identified volatiles, 3-FurAldehyde (Fur) and (E)-2-Heptenal (eHep) effectively inhibit the growth of plant pathogens such as Aspergillus flavus, Alternaria alternata, Fusarium solani and Pseudomonas syringae. Additionally, compounds such as Camphene (Cam), 3-FurAldehyde, and (E)-2-Heptenal, suppress the germination of seeds in crops including soybean, rice, maize, and wheat. Importantly, some of these VOCs also prevent rice seeds from pre-harvest sprouting. Consequently, our findings offer straightforward and practical approaches to seed protection and the reduction of PHS and PHFL, indicating potential new pathways for breeding, and reducing both PHS and pesticide usage in agriculture.

The type III effector NopL interacts with GmREM1a and GmNFR5 to promote symbiosis in soybean.

The establishment of symbiotic interactions between leguminous plants and rhizobia requires complex cellular programming activated by Rhizobium Nod factors (NFs) as well as type III effector (T3E)-mediated symbiotic signaling. However, the mechanisms by which different signals jointly affect symbiosis are still unclear. Here we describe the mechanisms mediating the cross-talk between the broad host range rhizobia Sinorhizobium fredii HH103 T3E Nodulation Outer Protein L (NopL) effector and NF signaling in soybean. NopL physically interacts with the Glycine max Remorin 1a (GmREM1a) and the NFs receptor NFR5 (GmNFR5) and promotes GmNFR5 recruitment by GmREM1a. Furthermore, NopL and NF influence the expression of GmRINRK1, a receptor-like kinase (LRR-RLK) ortholog of the Lotus RINRK1, that mediates NF signaling. Taken together, our work indicates that S. fredii NopL can interact with the NF signaling cascade components to promote the symbiotic interaction in soybean.

Identification of QTLs for seed and pod traits in soybean and analysis for additive effects and epistatic effects of QTLs among multiple environments.

Soybean seed and pod traits are important yield components. Selection for high yield style in seed and pod along with agronomic traits is a goal of many soybean breeders. The intention of this study was to identify quantitative trait loci (QTL) underlying seed and pod traits in soybean among eleven environments in China. 147 recombinant inbred lines were advanced through single-seed-descent method. The population was derived from a cross between Charleston (an American high yield soybean cultivar) and DongNong594 (a Chinese high yield soybean cultivar). A total of 157 polymorphic simple sequence repeat markers were used to construct a genetic linkage map. The phenotypic data of seed and pod traits [number of one-seed pod, number of two-seed pod, number of three-seed pod, number of four-seed pod, number of (two plus three)-seed pod, number of (three plus four)-seed pod, seed weight per plant, number of pod per plant] were recorded in eleven environments. In the analysis of single environment, fourteen main effect QTLs were identified. In the conjoint analysis of multiple environments, twenty-four additive QTLs were identified, and additive QTLs by environments interactions (AE) were evaluated and analyzed at the same time among eleven environments; twenty-three pairs of epistatic QTLs were identified, and epistasis (additive by additive) by environments interactions (AAE) were also analyzed and evaluated among eleven environments. Comparing the results of identification between single environment mapping and multiple environments conjoint mapping, three main effect QTLs with positive additive values and another three main effect QTLs with negative additive values, had no interactions with all environments, supported that these QTLs could be used in molecular assistant breeding in the future. These different effect QTLs could supply a good foundation to the gene clone and molecular asisstant breeding of soybean seed and pod traits.

Meta-Analysis and Multiomics of a Chromosome Segment Substitution Line Reveal Candidate Genes Associated with Seed Hardness in Soybean.

Soybean seed hardness is a key trait that influences planting, nutritional quality, and postharvest processing, but its genetic and molecular mechanisms remain to be clarified. We used meta-analysis to detect 17 meta-quantitative trait locus (QTLs) for soybean seed hardness. We then identified a hard-seeded chromosome segment substitution line, R75, with fragments introduced from hard-seeded wild germplasm in four of the meta-QTL intervals. Observations of the seed coat ultrastructure revealed thicker palisade tissue in R75 than in its soft-seeded recurrent parent. Transcriptomics and proteomics of R75 and its recurrent parent revealed multiple candidate genes associated with seed hardness. Fifty-seven were located on homozygous introduced fragments, 26 in meta-QTL intervals, and one in both (Glyma.02G268600). Five initial candidates were selected for KASP marker development on the basis of their predicted functions and nonsynonymous SNPs. The selection efficiency of the markers was as high as 90% for nonhard lines and 43% for hard lines in the chromosome segment substitution line (CSSL) population.

An analysis of differentially expressed and differentially m6A-modified transcripts in soybean roots treated with lead.

Lead is one of the most common toxic heavy metal pollutants in nature, and exposure to lead can cause serious toxicity to many organisms. In this study, we collected root growth data from soybean plants exposed to lead for seven days and confirmed that lead significantly inhibited root growth. We performed a transcriptome-wide m6A methylation analysis to study the response of soybean RNA methylation groups to lead. The m6A modified regions were enriched near the 3'UTR region and stop codon, and m6A methylation was positively correlated with transcript abundance. In the presence of lead, the transcriptome range of m6A RNA methylation peaks increased, and we identified 1144 m6A modification peaks and 1094 differentially expressed genes. The integration of m6A methylation and transcriptomic results enabled us to identify 16 candidate genes whose transcripts were differentially methylated and differentially expressed under lead stress. Annotation results suggest that these genes may promote abiotic stress tolerance by impacting lead uptake, transport, and accumulation through ROS pathways, enzymes, transporters, and hormones. These results provide candidate genes for future studies of lead stress tolerance mechanisms in soybean roots and provide genetic resources for studying plant heavy metal stress in soybean breeding.

Changes in the m6A RNA methylome accompany the promotion of soybean root growth by rhizobia under cadmium stress.

Cadmium (Cd) is the most widely distributed heavy metal pollutant in soil and has significant negative effects on crop yields and human health. Rhizobia can enhance soybean growth in the presence of heavy metals, and the legume-rhizobia symbiosis has been used to promote heavy-metal phytoremediation, but much remains to be learned about the molecular networks that underlie these effects. Here, we demonstrated that soybean root growth was strongly suppressed after seven days of Cd exposure but that the presence of rhizobia largely eliminated this effect, even prior to nodule development. Moreover, rhizobia did not appear to promote root growth by limiting plant Cd uptake: seedlings with and without rhizobia had similar root Cd concentrations. Previous studies have demonstrated a role for m6A RNA methylation in the response of rice and barley to Cd stress. We therefore performed transcriptome-wide m6A methylation profiling to investigate changes in the soybean RNA methylome in response to Cd with and without rhizobia. Here, we provide some of the first data on transcriptome-wide m6a RNA methylation patterns in soybean; m6A modifications were concentrated at the 3' UTR of transcripts and showed a positive relationship with transcript abundance. Transcriptome-wide m6A RNA methylation peaks increased in the presence of Cd, and the integration of m6A methylome and transcriptome results enabled us to identify 154 genes whose transcripts were both differentially methylated and differentially expressed in response to Cd stress. Annotation results suggested that these genes were associated with Ca2+ homeostasis, ROS pathways, polyamine metabolism, MAPK signaling, hormones, and biotic stress responses. There were 176 differentially methylated and expressed transcripts under Cd stress in the presence of rhizobia. In contrast to the Cd-only gene set, they were also enriched in genes related to auxin, jasmonic acid, and brassinosteroids, as well as abiotic stress tolerance. They contained fewer genes related to Ca2+ homeostasis and also included candidates with known functions in the legume-rhizobia symbiosis. These findings offer new insights into how rhizobia promote soybean root growth under Cd stress; they provide candidate genes for research on plant heavy metal responses and for the use of legumes in phytoremediation.

Multi-environment mapping and meta-analysis of 100-seed weight in soybean.

100-Seed weight (100-SW) of soybean is an important but complicated quantitative trait to yield. This study was focus on the quantitative trait loci (QTLs) of soybean 100-SW from 2006 to 2010, using recombination inbred lines population that was derived from a cross between Charleston and Dongnong 594. A total of 23 QTLs for 100-SW were detected in the linkage group C2, D1a, F, G and O. Nine QTLs were identified by composite interval mapping including one QTL with the minimum confidence interval (CI) of 1.3 cM, while 14 QTLs by multiple interval mapping. Furthermore, 94 reported QTLs of 100-SW were integrated with our QTL mapping results using BioMercator. As a result, 15 consensus QTLs and their corresponding markers were identified. The minimum CI was reduced to 1.52 cM by the combination of meta-analysis. These findings may merit fine-mapping of these QTL in soybean.