Isolation and Identification of Salinity-Tolerant Rhizobia and Nodulation Phenotype Analysis in Different Soybean Germplasms.

Increasing the soybean-planting area and increasing the soybean yield per unit area are two effective solutions to improve the overall soybean yield. Northeast China has a large saline soil area, and if soybeans could be grown there with the help of isolated saline-tolerant rhizobia, the soybean cultivation area in China could be effectively expanded. In this study, soybeans were planted in soils at different latitudes in China, and four strains of rhizobia were isolated and identified from the soybean nodules. According to the latitudes of the soil-sampling sites from high to low, the four isolated strains were identified as HLNEAU1, HLNEAU2, HLNEAU3, and HLNEAU4. In this study, the isolated strains were identified for their resistances, and their acid and saline tolerances and nitrogen fixation capacities were preliminarily identified. Ten representative soybean germplasm resources in Northeast China were inoculated with these four strains, and the compatibilities of these four rhizobium strains with the soybean germplasm resources were analyzed. All four isolates were able to establish different extents of compatibility with 10 soybean resources. Hefeng 50 had good compatibility with the four isolated strains, while Suinong 14 showed the best compatibility with HLNEAU2. The isolated rhizobacteria could successfully establish symbiosis with the soybeans, but host specificity was also present. This study was a preliminary exploration of the use of salinity-tolerant rhizobacteria to help the soybean nitrogen fixation in saline soils in order to increase the soybean acreage, and it provides a valuable theoretical basis for the application of saline-tolerant rhizobia.

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.

QTL Mapping and Data Mining to Identify Genes Associated with Soybean Epicotyl Length Using Cultivated Soybean and Wild Soybean.

Soybean (Glycine max) plants first emerged in China, and they have since been established as an economically important oil crop and a major source of daily protein for individuals throughout the world. Seed emergence height is the first factor that ensures seedling adaptability to field management practices, and it is closely related to epicotyl length. In the present study, the Suinong 14 and ZYD00006 soybean lines were used as parents to construct chromosome segment substitution lines (CSSLs) for quantitative trait loci (QTL) identification. Seven QTLs were identified using two years of epicotyl length measurement data. The insertion region of the ZYD00006 fragment was identified through whole genome resequencing, with candidate gene screening and validation being performed through RNA-Seq and qPCR, and Glyma.08G142400 was ultimately selected as an epicotyl length-related gene. Through combined analyses of phenotypic data from the study population, Glyma.08G142400 expression was found to be elevated in those varieties exhibiting longer epicotyl length. Haplotype data analyses revealed that epicotyl data were consistent with haplotype typing. In summary, the QTLs found to be associated with the epicotyl length identified herein provide a valuable foundation for future molecular marker-assisted breeding efforts aimed at improving soybean emergence height in the field, with the Glyma.08G142400 gene serving as a regulator of epicotyl length, offering new insight into the mechanisms that govern epicotyl development.

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.

Detection of type III effector-induced transcription factors that regulate phytohormone content during symbiosis establishment in soybean.

Soybean is a pivotal protein and oil crop that utilizes atmospheric nitrogen via symbiosis with rhizobium soil bacteria. Rhizobial type III effectors (T3Es) are essential regulators during symbiosis establishment. However, how the transcription factors involved in the interaction between phytohormone synthesis and type III effectors are connected is unclear. To detect the responses of phytohormone and transcription factor genes to rhizobial type III effector NopAA and type III secretion system, the candidate genes underlying soybean symbiosis were identified using RNA sequencing (RNA-seq) and phytohormone content analysis of soybean roots infected with wild-type Rhizobium and its derived T3E mutant. Via RNA-seq analysis the WRKY and ERF transcription factor families were identified as the most differentially expressed factors in the T3E mutant compared with the wild-type. Next, qRT-PCR was used to confirm the candidate genes Glyma.09g282900, Glyma.08g018300, Glyma.18g238200, Glyma.03g116300, Glyma.07g246600, Glyma.16g172400 induced by S. fredii HH103, S. fredii HH103ΩNopAA, and S. fredii HH103ΩRhcN. Since the WRKY and ERF families may regulate abscisic acid (ABA) content and underlying nodule formation, we performed phytohormone content analysis at 0.5 and 24 h post-inoculation (hpi). A significant change in ABA content was found between wild Rhizobium and type III effector mutant. Our results support that NopAA can promote the establishment of symbiosis by affecting the ABA signaling pathways by regulating WRKY and ERF which regulate the phytohormone signaling pathway. Specifically, our work provides insights into a signaling interaction of prokaryotic effector-induced phytohormone response involved in host signaling that regulates the establishment of symbiosis and increases nitrogen utilization efficiency in soybean plants.

NopC/T/L Signal Crosstalk Gene GmPHT1-4.

Symbiotic nodulation between leguminous plants and rhizobia is a critical biological interaction. The type III secretion system (T3SS) employed by rhizobia manipulates the host's nodulation signaling, analogous to mechanisms used by certain bacterial pathogens for effector protein delivery into host cells. This investigation explores the interactive signaling among type III effectors HH103ΩNopC, HH103ΩNopT, and HH103ΩNopL from SinoRhizobium fredii HH103. Experimental results revealed that these effectors positively regulate nodule formation. Transcriptomic analysis pinpointed GmPHT1-4 as the key gene facilitating this effector-mediated signaling. Overexpression of GmPHT1-4 enhances nodulation, indicating a dual function in nodulation and phosphorus homeostasis. This research elucidates the intricate regulatory network governing Rhizobium-soybean (Glycine max (L.) Merr) interactions and the complex interplay between type III effectors.

Identifying the Soybean microRNAs Related to Phytophthora sojae Based on RNA Sequencing and Bioinformatics Analysis.

Phytophthora root rot in soybeans is caused by a pathogen called Phytophthora sojae (P. sojae), which results in a significant decrease in soybean production within affected regions. MicroRNAs (miRNAs) are a class of small non-coding RNA molecules that play a key post-transcriptional regulatory role in eukaryotes. In this paper, the miRNAs that respond to P. sojae were analyzed from the gene level to complement the study of molecular resistance mechanisms in soybean. The study utilized high-throughput sequencing of soybean data to predict miRNAs that respond to P. sojae, analyze their specific functions, and verify regulatory relationships using qRT-PCR. The results showed that the miRNAs in soybean respond to P. sojae infection. MiRNAs can be transcribed independently, suggesting the presence of transcription factor binding sites in the promoter regions. Additionally, we performed an evolutionary analysis on conserved miRNAs that respond to P. sojae. Finally, we investigated the regulatory relationships among miRNAs, genes, and transcription factors, and identified five regulatory patterns. These findings lay the groundwork for future studies on the evolution of miRNAs responsive to P. sojae.

Identifications of QTLs and Candidate Genes Associated with Pseudomonas syringae Responses in Cultivated Soybean (Glycine max) and Wild Soybean (Glycine soja).

Soybeans (Glycine max) are a key food crop, serving as a valuable source of both oil and plant-derived protein. Pseudomonas syringae pv. glycinea (Psg) is among the most aggressive and prevalent pathogens affecting soybean production, causing a form of bacterial spot disease that impacts soybean leaves and thereby reduces crop yields. In this study, 310 natural soybean varieties were screened for Psg resistance and susceptibility. The identified susceptible and resistant varieties were then used for linkage mapping, BSA-seq, and whole genome sequencing (WGS) analyses aimed at identifying key QTLs associated with Psg responses. Candidate Psg-related genes were further confirmed through WGS and qPCR analyses. Candidate gene haplotype analyses were used to explore the associations between haplotypes and soybean Psg resistance. In addition, landrace and wild soybean plants were found to exhibit a higher degree of Psg resistance as compared to cultivated soybean varieties. In total, 10 QTLs were identified using chromosome segment substitution lines derived from Suinong14 (cultivated soybean) and ZYD00006 (wild soybean). Glyma.10g230200 was found to be induced in response to Psg, with the Glyma.10g230200 haplotype corresponding to soybean disease resistance. The QTLs identified herein can be leveraged to guide the marker-assisted breeding of soybean cultivars that exhibit partial resistance to Psg. Moreover, further functional and molecular studies of Glyma.10g230200 have the potential to offer insight into the mechanistic basis for soybean Psg resistance.

NopAA and NopD Signaling Association-Related Gene GmNAC27 Promotes Nodulation in Soybean (Glycine max).

Rhizobia secrete effectors that are essential for the effective establishment of their symbiotic interactions with leguminous host plants. However, the signaling pathways governing rhizobial type III effectors have yet to be sufficiently characterized. In the present study, the type III effectors, NopAA and NopD, which perhaps have signaling pathway crosstalk in the regulation of plant defense responses, have been studied together for the first time during nodulation. Initial qRT-PCR experiments were used to explore the impact of NopAA and NopD on marker genes associated with symbiosis and defense responses. The effects of these effectors on nodulation were then assessed by generating bacteria in which both NopAA and NopD were mutated. RNA-sequencing analyses of soybean roots were further utilized to assess signaling crosstalk between NopAA and NopD. NopAA mutant and NopD mutant were both found to repress GmPR1, GmPR2, and GmPR5 expression in these roots. The two mutants also significantly reduced nodules dry weight and the number of nodules and infection threads, although these changes were not significantly different from those observed following inoculation with double-mutant (HH103ΩNopAA&NopD). NopAA and NopD co-mutant inoculation was primarily found to impact the plant-pathogen interaction pathway. Common differentially expressed genes (DEGs) associated with both NopAA and NopD were enriched in the plant-pathogen interaction, plant hormone signal transduction, and MAPK signaling pathways, and no further changes in these common DEGs were noted in response to inoculation with HH103ΩNopAA&NopD. Glyma.13G279900 (GmNAC27) was ultimately identified as being significantly upregulated in the context of HH103ΩNopAA&NopD inoculation, serving as a positive regulator of nodulation. These results provide new insight into the synergistic impact that specific effectors can have on the establishment of symbiosis and the responses of host plant proteins.

GmTCP and GmNLP Underlying Nodulation Character in Soybean Depending on Nitrogen.

Soybean is a cereal crop with high protein and oil content which serves as the main source of plant-based protein and oil for human consumption. The symbiotic relationship between legumes and rhizobia contributes significantly to soybean yield and quality, but the underlying molecular mechanisms remain poorly understood, hindering efforts to improve soybean productivity. In this study, we conducted a transcriptome analysis and identified 22 differentially expressed genes (DEGs) from nodule-related quantitative trait loci (QTL) located in chromosomes 12 and 19. Subsequently, we performed functional characterisation and haplotype analysis to identify key candidate genes among the 22 DEGs that are responsive to nitrate. Our findings identified GmTCP (TEOSINTE-BRANCHED1/CYCLOIDEA/PCF) and GmNLP (NIN-LIKE PROTEIN) as the key candidate genes that regulate the soybean nodule phenotype in response to nitrogen concentration. We conducted homologous gene mutant analysis in Arabidopsis thaliana, which revealed that the homologous genes of GmTCP and GmNLP play a vital role in regulating root development in response to nitrogen concentration. We further performed overexpression and gene knockout of GmTCP and GmNLP through hairy root transformation in soybeans and analysed the effects of GmTCP and GmNLP on nodulation under different nitrogen concentrations using transgenic lines. Overexpressing GmTCP and GmNLP resulted in significant differences in soybean hairy root nodulation phenotypes, such as nodule number (NN) and nodule dry weight (NDW), under varying nitrate conditions. Our results demonstrate that GmTCP and GmNLP are involved in regulating soybean nodulation in response to nitrogen concentration, providing new insights into the mechanism of soybean symbiosis establishment underlying different nitrogen concentrations.

Effector-Dependent and -Independent Molecular Mechanisms of Soybean-Microbe Interaction.

Soybean is a pivotal staple crop worldwide, supplying the main food and feed plant proteins in some countries. In addition to interacting with mutualistic microbes, soybean also needs to protect itself against pathogens. However, to grow inside plant tissues, plant defense mechanisms ranging from passive barriers to induced defense reactions have to be overcome. Pathogenic but also symbiotic micro-organisms effectors can be delivered into the host cell by secretion systems and can interfere with the immunity system and disrupt cellular processes. This review summarizes the latest advances in our understanding of the interaction between secreted effectors and soybean feedback mechanism and uncovers the conserved and special signaling pathway induced by pathogenic soybean cyst nematode, Pseudomonas, Xanthomonas as well as by symbiotic rhizobium.

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.

GsRSS3L, a Candidate Gene Underlying Soybean Resistance to Seedcoat Mottling Derived from Wild Soybean (Glycine soja Sieb. and Zucc).

Soybeans are a major crop that produce the best vegetable oil and protein for use in food and beverage products worldwide. However, one of the most well-known viral infections affecting soybeans is the Soybean Mosaic Virus (SMV), a member of the Potyviridae family. A crucial method for preventing SMV damage is the breeding of resistant soybean cultivars. Adult resistance and resistance of seedcoat mottling are two types of resistance to SMV. Most studies have focused on adult-plant resistance but not on the resistance to seedcoat mottling. In this study, chromosome segment-substituted lines derived from a cross between Suinong14 (cultivated soybean) and ZYD00006 (wild soybean) were used to identify the chromosome region and candidate genes underlying soybean resistance to seed coat mottling. Herein, two quantitative trait loci (QTLs) were found on chromosome 17, and eighteen genes were found in the QTL region. RNA-seq was used to evaluate the differentially expressed genes (DEGs) among the eighteen genes located in the QTLs. According to the obtained data, variations were observed in the expression of five genes following SMV infection. Furthermore, Nicotiana benthamiana was subjected to an Agrobacterium-mediated transient expression assay to investigate the role of the five candidate genes in SMV resistance. It has also been revealed that Glyma.17g238900 encoding a RICE SALT SENSITIVE 3-like protein (RSS3L) can inhibit the multiplication of SMV in N.benthamiana. Moreover, two nonsynonymous single-nucleotide polymorphisms (SNPs) were found in the coding sequence of Glyma.17g238900 derived from the wild soybean ZYD00006 (GsRSS3L), and the two amino acid mutants may be associated with SMV resistance. Hence, it has been suggested that GsRSS3L confers seedcoat mottling resistance, shedding light on the mechanism of soybean resistance to SMV.

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.

Combined Linkage Mapping and BSA to Identify QTL and Candidate Genes for Plant Height and the Number of Nodes on the Main Stem in Soybean.

Soybean is one of the most important food and oil crops in the world. Plant height (PH) and the number of nodes on the main stem (NNMS) are quantitative traits closely related to soybean yield. In this study, we used 208 chromosome segment substitution lines (CSSL) populations constructed using "SN14" and "ZYD00006" for quantitative trait locus (QTL) mapping of PH and NNMS. Combined with bulked segregant analysis (BSA) by extreme materials, 8 consistent QTLs were identified. According to the gene annotation of the QTL interval, a total of 335 genes were obtained. Five of which were associated with PH and NNMS, potentially representing candidate genes. RT-qPCR of these 5 candidate genes revealed two genes with differential relative expression levels in the stems of different materials. Haplotype analysis showed that different single nucleotide polymorphisms (SNPs) between the excellent haplotypes in Glyma.04G251900 and Glyma.16G156700 may be the cause of changes in these traits. These results provide the basis for research on candidate genes and marker-assisted selection (MAS) in soybean breeding.

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.