QTL analysis of seed-flooding tolerance in soybean (Glycine max [L.] Merr.)

Subcellular Proteomics to Elucidate Soybean Response to Abiotic Stress

Climate change jeopardizes soybean production by declining seed yield and quality. In this review, the morphophysiological alterations of soybean in response to abiotic stress are summarized, followed by illustrations of cellular metabolisms and regulatory mechanisms to organellar stress based on subcellular proteomics. This highlights the communications associated with reactive oxygen species scavenging, molecular chaperones, and phytohormone signals among subcellular compartments. Given the complexity of climate change and the limitations of plants in coping with multiple abiotic stresses, a generic response to environmental constraints is proposed between calcium and abscisic acid signals in subcellular organelles. This review summarizes the findings of subcellular proteomics in stressed soybean and discusses the future prospects of subcellular proteomics for promoting the improvement of climate-tolerant crops.

Genetic Dissection of Extreme Seed-Flooding Tolerance in a Wild Soybean PI342618B by Linkage Mapping and Candidate Gene Analysis

Seed-flooding stress is one of major abiotic constraints that adversely affects soybean production worldwide. Identifying tolerant germplasms and revealing the genetic basis of seed-flooding tolerance are imperative goals for soybean breeding. In the present study, high-density linkage maps of two inter-specific recombinant inbred line (RIL) populations, named NJIRNP and NJIR4P, were utilized to identify major quantitative trait loci (QTLs) for seed-flooding tolerance using three parameters viz., germination rate (GR), normal seedling rate (NSR), and electrical conductivity (EC). A total of 25 and 18 QTLs were detected by composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively, and 12 common QTLs were identified through both methods. All favorable alleles for the tolerance are notably from the wild soybean parent. Moreover, four digenic epistatic QTL pairs were identified, and three of them showed no main effects. In addition, the pigmented soybean genotypes exhibited high seed-flooding tolerance compared with yellow seed coat genotypes in both populations. Moreover, out of five identified QTLs, one major region containing multiple QTLs associated with all three traits was identified on Chromosome 8, and most of the QTLs within this hotspot were major loci (R² > 10) and detectable in both populations and multiple environments. Based on the gene expression and functional annotation information, 10 candidate genes from QTL "hotspot 8-2" were screened for further analysis. Furthermore, the results of qRT-PCR and sequence analysis revealed that only one gene, GmDREB2 (Glyma.08G137600), was significantly induced under flooding stress and displayed a TTC tribasic insertion mutation of the nucleotide sequence in the tolerant wild parent (PI342618B). GmDREB2 encodes an ERF transcription factor, and the subcellular localization analysis using green fluorescent protein (GFP) revealed that GmDREB2 protein was localized in the nucleus and plasma membrane. Furthermore, overexpression of GmDREB2 significantly promoted the growth of soybean hairy roots, which might indicate its critical role in seed-flooding stress. Thus, GmDREB2 was considered as the most possible candidate gene for seed-flooding tolerance.

The intervention of classical and molecular breeding approaches to enhance flooding stress tolerance in soybean - An review

Abiotic stresses and climate changes cause severe loss of yield and quality of crops and reduce the production area worldwide. Flooding stress curtails soybean growth, yield, and quality and ultimately threatens the global food supply chain. Flooding tolerance is a multigenic trait. Tremendous research in molecular breeding explored the potential genomic regions governing flood tolerance in soybean. The most robust way to develop flooding tolerance in soybean is by using molecular methods, including quantitative trait loci (QTL) mapping, identification of transcriptomes, transcription factor analysis, CRISPR/Cas9, and to some extent, genome-wide association studies (GWAS), and multi-omics techniques. These powerful molecular tools have deepened our knowledge about the molecular mechanism of flooding stress tolerance. Besides all this, using conventional breeding methods (hybridization, introduction, and backcrossing) and other agronomic practices is also helpful in combating the rising flooding threats to the soybean crop. The current review aims to summarize recent advancements in breeding flood-tolerant soybean, mainly by using molecular and conventional tools and their prospects. This updated picture will be a treasure trove for future researchers to comprehend the foundation of flooding tolerance in soybean and cover the given research gaps to develop tolerant soybean cultivars able to sustain growth under extreme climatic changes.

Identification of Functional Genetic Variations Underlying Flooding Tolerance in Brazilian Soybean Genotypes

Flooding is a frequent environmental stress that reduces soybean (Glycine max) growth and grain yield in many producing areas in the world, such as, e.g., in the United States, Southeast Asia and Southern Brazil. In these regions, soybean is frequently cultivated in lowland areas by rotating with rice (Oryza sativa), which provides numerous technical, economic and environmental benefits. Given these realities, this work aimed to characterize physiological responses, identify genes differentially expressed under flooding stress in Brazilian soybean genotypes with contrasting flooding tolerance, and select SNPs with potential use for marker-assisted selection. Soybean cultivars TECIRGA 6070 (flooding tolerant) and FUNDACEP 62 (flooding sensitive) were grown up to the V6 growth stage and then flooding stress was imposed. Total RNA was extracted from leaves 24 h after the stress was imposed and sequenced. In total, 421 induced and 291 repressed genes were identified in both genotypes. TECIRGA 6070 presented 284 and 460 genes up- and down-regulated, respectively, under flooding conditions. Of those, 100 and 148 genes were exclusively up- and down-regulated, respectively, in the tolerant genotype. Based on the RNA sequencing data, SNPs in differentially expressed genes in response to flooding stress were identified. Finally, 38 SNPs, located in genes with functional annotation for response to abiotic stresses, were found in TECIRGA 6070 and absent in FUNDACEP 62. To validate them, 22 SNPs were selected for designing KASP assays that were used to genotype a panel of 11 contrasting genotypes with known phenotypes. In addition, the phenotypic and grain yield impacts were analyzed in four field experiments using a panel of 166 Brazilian soybean genotypes. Five SNPs possibly related to flooding tolerance in Brazilian soybean genotypes were identified. The information generated from this research will be useful to develop soybean genotypes adapted to poorly drained soils or areas subject to flooding.

Quantitative Trait Locus Mapping of Seed Vigor in Soybean under -20 °C Storage and Accelerated Aging Conditions via RAD Sequencing

Due to its fast deterioration, soybean (Glycine max L.) has an inherently poor seed vigor. Vigor loss occurring during storage is one of the main obstacles to soybean production in the tropics. To analyze the genetic background of seed vigor, soybean seeds of a recombinant inbred line (RIL) population derived from the cross between Zhonghuang24 (ZH24, low vigor cultivar) and Huaxia3hao (HX3, vigorous cultivar) were utilized to identify the quantitative trait loci (QTLs) underlying the seed vigor under -20 °C conservation and accelerated aging conditions. According to the linkage analysis, multiple seed vigor-related QTLs were identified under both -20 °C and accelerated aging storage. Two major QTLs and eight QTL hotspots localized on chromosomes 3, 6, 9, 11, 15, 16, 17, and 19 were detected that were associated with seed vigor across two storage conditions. The indicators of seed vigor did not correlate well between the two aging treatments, and no common QTLs were detected in RIL populations stored in two conditions. These results indicated that deterioration under accelerated aging conditions was not reflective of natural aging at -20 °C. Additionally, we suggest 15 promising candidate genes that could possibly determine the seed vigor in soybeans, which would help explore the mechanisms responsible for maintaining high seed vigor.

Genome-wide association study uncovers major genetic loci associated with seed flooding tolerance in soybean

BACKGROUND

Seed flooding stress is one of the threatening environmental stressors that adversely limits soybean at the germination stage across the globe. The knowledge on the genetic basis underlying seed-flooding tolerance is limited. Therefore, we performed a genome-wide association study (GWAS) using 34,718 single nucleotide polymorphism (SNPs) in a panel of 243 worldwide soybean collections to identify genetic loci linked to soybean seed flooding tolerance at the germination stage.

RESULTS

In the present study, GWAS was performed with two contrasting models, Mixed Linear Model (MLM) and Multi-Locus Random-SNP-Effect Mixed Linear Model (mrMLM) to identify significant SNPs associated with electrical conductivity (EC), germination rate (GR), shoot length (ShL), and root length (RL) traits at germination stage in soybean. With MLM, a total of 20, 40, 4, and 9 SNPs associated with EC, GR, ShL and RL, respectively, whereas in the same order mrMLM detected 27, 17, 13, and 18 SNPs. Among these SNPs, two major SNPs, Gm_08_11971416, and Gm_08_46239716 were found to be consistently connected with seed-flooding tolerance related traits, namely EC and GR across two environments. We also detected two SNPs, Gm_05_1000479 and Gm_01_53535790 linked to ShL and RL, respectively. Based on Gene Ontology enrichment analysis, gene functional annotations, and protein-protein interaction network analysis, we predicted eight candidate genes and three hub genes within the regions of the four SNPs with Cis-elements in promoter regions which may be involved in seed-flooding tolerance in soybeans and these warrant further screening and functional validation.

CONCLUSIONS

Our findings demonstrate that GWAS based on high-density SNP markers is an efficient approach to dissect the genetic basis of complex traits and identify candidate genes in soybean. The trait associated SNPs could be used for genetic improvement in soybean breeding programs. The candidate genes could help researchers better understand the molecular mechanisms underlying seed-flooding stress tolerance in soybean.

Tolerance of four grain legume species to waterlogging, hypoxia and anoxia at germination and recovery

Legume seeds, when relay sown following rice, may suffer from soil waterlogging and the associated hypoxia or even anoxia. This study evaluated the tolerance of grain legume species, grass pea (three genotypes), lentil (two genotypes), faba bean (two genotypes) and field pea (one genotype), to soil waterlogging in a glasshouse, to anoxia and hypoxia in temperature-controlled room at germination and seedling stages. Changes in oxygen in the surface layers of soil, with time after waterlogging, were measured by microelectrode profiling. The soil profiling showed that soil oxygen declined and then stabilized by the fourth day after waterlogging and oxygen was not detected at 8 mm below the soil surface. Germination of seeds under waterlogging for up to 12 days and seedling survival after the soil was drained for up to 36 days were measured in pot experiments. Seed germination and/or survival in anoxia (N₂-flushed solutions) and hypoxia (1.0 and 2.5 kPa oxygen) were evaluated, and so were post-anoxia or post-hypoxia recoveries, all in comparison with aerated controls. Lentil had higher seedling emergence (55 %) than the other species during soil waterlogging. However, lentil had lower seedling survival (9 %) than grass pea (28 %) during recovery following soil drainage. Grass pea seeds were more tolerant of anoxia and of hypoxia than the seeds of the three other species. In conclusion, grass pea, with higher percent germination and seedling survival during recovery, is more tolerant to waterlogging and subsequent soil drainage than the three other grain legume species. Grass pea was also more tolerant of hypoxia and of anoxia at the seed germination stage. These findings demonstrate the superior waterlogging tolerance of grass pea in relay sowing, as compared with the other grain legumes.

Hydropriming Applied on Fast Germinating Solanum villosum Miller Seeds: Impact on Pre-germinative Metabolism

Seed priming can circumvent poor germination rate and uniformity, frequently reported in eggplant (Solanum melongena L.) and its crop wild relatives (CWRs). However, there is still a gap of knowledge on how these treatments impact the pre-germinative metabolism in a genotype- and/or species-dependent manner. The CWR Solanum villosum Miller (hairy nightshade) investigated in this study showed a quite unique profile of fast germination. Although this accelerated germination profile would not apparently require further improvement, we wanted to test whether priming would still be able to impact the pre-germinative metabolism, eventually disclosing the predominant contribution of specific antioxidant components. Hydropriming followed by dry-back resulted in synchronized germination, as revealed by the lowest MGR (Mean Germination Rate) and U (Uncertainty) values, compared to unprimed seeds. No significant changes in ROS (reactive oxygen species) were observed throughout the treatment. Increased tocopherols levels were detected at 2 h of hydropriming whereas, overall, a low lipid peroxidation was evidenced by the malondialdehyde (MDA) assay. Hydropriming resulted in enhanced accumulation of the naturally occurring antioxidant phenolic compounds chlorogenic acid and iso-orientin, found in the dry seeds and ex novo accumulation of rutin. The dynamic changes of the pre-germinative metabolism induced by hydropriming are discussed in view of future applications that might boost the use of eggplant CWRs for breeding, upon upgrade mediated by seed technology.

Prioritization and Evaluation of Flooding Tolerance Genes in Soybean [Glycine max (L.) Merr.]

Soybean [Glycine max (L.) Merr.] is one of the most important legume crops abundant in edible protein and oil in the world. In recent years there has been increasingly more drastic weather caused by climate change, with flooding, drought, and unevenly distributed rainfall gradually increasing in terms of the frequency and intensity worldwide. Severe flooding has caused extensive losses to soybean production and there is an urgent need to breed strong soybean seeds with high flooding tolerance. The present study demonstrates bioinformatics big data mining and integration, meta-analysis, gene mapping, gene prioritization, and systems biology for identifying prioritized genes of flooding tolerance in soybean. A total of 83 flooding tolerance genes (FTgenes), according to the appropriate cut-off point, were prioritized from 36,705 test genes collected from multidimensional genomic features linking to soybean flooding tolerance. Several validation results using independent samples from SoyNet, genome-wide association study, SoyBase, GO database, and transcriptome databases all exhibited excellent agreement, suggesting these 83 FTgenes were significantly superior to others. These results provide valuable information and contribution to research on the varieties selection of soybean.

Review: Proteomic Techniques for the Development of Flood-Tolerant Soybean

Soybean, which is rich in protein and oil as well as phytochemicals, is cultivated in several climatic zones. However, its growth is markedly decreased by flooding stress, which is caused by climate change. Proteomic techniques were used for understanding the flood-response and -tolerant mechanisms in soybean. Subcellular proteomics has potential to elucidate localized cellular responses and investigate communications among subcellular components during plant growth and under stress stimuli. Furthermore, post-translational modifications play important roles in stress response and tolerance to flooding stress. Although many flood-response mechanisms have been reported, flood-tolerant mechanisms have not been fully clarified for soybean because of limitations in germplasm with flooding tolerance. This review provides an update on current biochemical and molecular networks involved in soybean tolerance against flooding stress, as well as recent developments in the area of functional genomics in terms of developing flood-tolerant soybeans. This work will expedite marker-assisted genetic enhancement studies in crops for developing high-yielding stress-tolerant lines or varieties under abiotic stress.

Quantitative proteomic analyses of two soybean low phytic acid mutants to identify the genes associated with seed field emergence

BACKGROUND

Seed germination is essential to crop growth and development, and ultimately affects its harvest. It is difficult to breed soybeans low in phytic acid with a higher seed field emergence. Although additional management and selection could overcome the phytate reduction, the mechanisms of seed germination remain unknown.

RESULTS

A comparative proteomic analysis was conducted between two low phytic acid (LPA) soybean mutants (TW-1-M and TW-1), both of which had a deletion of 2 bp in the GmMIPS1 gene. However, the TW-1 seeds showed a significantly lower field emergence compared to the TW-1-M. There were 282 differentially accumulated proteins (DAPs) identified between two mutants at the three stages. Among these DAPs, 80 were down-accumulated and 202 were up-accumulated. Bioinformatic analysis showed that the identified proteins were related to functional categories of oxidation reduction, response to stimulus and stress, dormancy and germination processes and catalytic activity. KEGG analysis showed that these DAPs were mainly involved in energy metabolism and anti-stress pathways. Based upon the conjoint analysis of DAPs with the differentially expressed genes (DEGs) previously published among three germination stages in two LPA mutants, 30 shared DAPs/DEGs were identified with different patterns, including plant seed protein, beta-amylase, protein disulfide-isomerase, disease resistance protein, pyrophosphate-fructose 6-phosphate 1-phosphotransferase, cysteine proteinase inhibitor, non-specific lipid-transfer protein, phosphoenolpyruvate carboxylase and acyl-coenzyme A oxidase.

CONCLUSIONS

Seed germination is a very complex process in LPA soybean mutants. The TW-1-M and TW-1 showed many DAPs involved in seed germination. The differential accumulation of these proteins could result in the difference of seed field emergence between the two mutants. The high germination rate in the TW-1-M might be strongly attributed to reactive oxygen species-related and plant hormone-related genes. All these findings would help us further explore the germination mechanisms in LPA crops.

Identification of QTN and Candidate Gene for Seed-flooding Tolerance in Soybean [Glycine max (L.) Merr.] using Genome-Wide Association Study (GWAS)

Seed-flooding stress is one of the major abiotic constraints severely affecting soybean yield and quality. Understanding the molecular mechanism and genetic basis underlying seed-flooding tolerance will be of greatly importance in soybean breeding. However, very limited information is available about the genetic basis of seed-flooding tolerance in soybean. The present study performed Genome-Wide Association Study (GWAS) to identify the quantitative trait nucleotides (QTNs) associated with three seed-flooding tolerance related traits, viz., germination rate (GR), normal seedling rate (NSR) and electric conductivity (EC), using a panel of 347 soybean lines and the genotypic data of 60,109 SNPs with MAF > 0.05. A total of 25 and 21 QTNs associated with all three traits were identified via mixed linear model (MLM) and multi-locus random-SNP-effect mixed linear model (mrMLM) in three different environments (JP14, HY15, and Combined). Among these QTNs, three major QTNs, viz., QTN13, qNSR-10 and qEC-7-2, were identified through both methods MLM and mrMLM. Interestingly, QTN13 located on Chr.13 has been consistently identified to be associated with all three studied traits in both methods and multiple environments. Within the 1.0 Mb physical interval surrounding the QTN13, nine candidate genes were screened for their involvement in seed-flooding tolerance based on gene annotation information and available literature. Based on the qRT-PCR and sequence analysis, only one gene designated as GmSFT (Glyma.13g248000) displayed significantly higher expression level in all tolerant genotypes compared to sensitive ones under flooding treatment, as well as revealed nonsynonymous mutation in tolerant genotypes, leading to amino acid change in the protein. Additionally, subcellular localization showed that GmSFT was localized in the nucleus and cell membrane. Hence, GmSFT was considered as the most likely candidate gene for seed-flooding tolerance in soybean. In conclusion, the findings of the present study not only increase our knowledge of the genetic control of seed-flooding tolerance in soybean, but will also be of great utility in marker-assisted selection and gene cloning to elucidate the mechanisms of seed-flooding tolerance.

Identification of Genes/Proteins Related to Submergence Tolerance by Transcriptome and Proteome Analyses in Soybean

Flooding can lead to yield reduction of soybean. Therefore, identification of flooding tolerance genes has great significance in production practice. In this study, Qihuang 34, a highly-resistant variety to flooding stress, was selected for submergence treatments. Transcriptome and proteome analyses were conducted, by which twenty-two up-regulated differentially expressed genes (DEGs)/differentially expressed proteins (DEPs) associated with five KEGG pathways were isolated. The number of the DEGs/DEPs enriched in glycolysis/gluconeogenesis was the highest. Four of these genes were confirmed by RT-qPCR, suggesting that glycolysis/gluconeogenesis may be activated to generate energy for plant survival under anaerobic conditions. Thirty-eight down-regulated DEGs/DEPs associated with six KEGG pathways were identified under submergence stress. Eight DEGs/DEPs enriched in phenylpropanoid biosynthesis were assigned to peroxidase, which catalyzes the conversion of coumaryl alcohol to hydroxy-phenyl lignin in the final step of lignin biosynthesis. Three of these genes were confirmed by RT-qPCR. The decreased expression of these genes led to the inhibition of lignin biosynthesis, which may be the cause of plant softening under submergence stress for a long period of time. This study revealed a number of up-/down-regulated pathways and the corresponding DEGs/DEPs, by which, a better understanding of the mechanisms of submergence tolerance in soybean may be achieved.

Waterlogging Tolerance at Germination in Field Pea: Variability, Genetic Control, and Indirect Selection

In the Eastern Gangetic Plain of South Asia field pea (Pisum sativum L.) is often grown as a relay crop where soil waterlogging (WL) causes germination failure. To assess if selection for WL tolerance is feasible, we studied the response to WL stress at germination stage in a recombinant inbred line (RIL) population from a bi-parental cross between WL-contrasting parents and in a diversity panel to identify extreme phenotypes, understand the genetics of WL tolerance and find traits for possible use in indirect selection. The RIL population and the diversity panel were screened to test the ability of germination under both waterlogged and drained soils. A total of 50, most WL tolerant and sensitive, genotypes from each of both the RIL and the diversity panel were further evaluated to assay testa integrity/leakage in CaSO₄ solution. Morphological characterization of both populations was undertaken. A wide range of variation in the ability to germination in waterlogged soil was observed in the RIL population (6-93%) and the diversity panel (5-100%) with a high broad-sense heritability (H ² > 85%). The variation was continuously distributed indicating polygenic control. Most genotypes with a dark colored testa (90%) were WL tolerant, whereas those with a light colored testa were all WL sensitive in both the RIL population and diversity panel. Testa integrity, measured by electrical conductivity of the leakage solute, was strongly associated with WL tolerance in the RIL population (r _G = -1.00) and the diversity panel (r _G = -0.90). Therefore, testa integrity can be effectively used in indirect selection for WL tolerance. Response to selection for WL tolerance at germination is confidently predicted enabling the adaptation of the ancient model pea to extreme precipitation events at germination.

Improving Flooding Tolerance of Crop Plants

A major problem of climate change is the increasing duration and frequency of heavy rainfall events. This leads to soil flooding that negatively affects plant growth, eventually leading to death of plants if the flooding persists for several days. Most crop plants are very sensitive to flooding, and dramatic yield losses occur due to flooding each year. This review summarizes recent progress and approaches to enhance crop resistance to flooding. Most experiments have been done on maize, barley, and soybean. Work on other crops such as wheat and rape has only started. The most promising traits that might enhance crop flooding tolerance are anatomical adaptations such as aerenchyma formation, the formation of a barrier against radial oxygen loss, and the growth of adventitious roots. Metabolic adaptations might be able to improve waterlogging tolerance as well, but more studies are needed in this direction. Reasonable approaches for future studies are quantitative trait locus (QTL) analyses or genome-wide association (GWA) studies in combination with specific tolerance traits that can be easily assessed. The usage of flooding-tolerant relatives or ancestral cultivars of the crop of interest in these experiments might enhance the chances of finding useful tolerance traits to be used in breeding.

Genetic Analysis of Flooding Tolerance in an Andean Diversity Panel of Dry Bean (Phaseolus vulgaris L.)

Climate change models predict temporal and spatial shifts in precipitation resulting in more frequent incidents of flooding, particularly in regions with poor soil drainage. In these flooding conditions, crop losses are inevitable due to exposure of plants to hypoxia and the spread of root rot diseases. Improving the tolerance of bean cultivars to flooding is crucial to minimize crop losses. In this experiment, we evaluated the phenotypic responses of 277 genotypes from the Andean Diversity Panel to flooding at germination and seedling stages. A randomized complete block design, with a split plot arrangement, was employed for phenotyping germination rate, total weight, shoot weight, root weight, hypocotyl length, SPAD index, adventitious root rate, and survival score. A subset of genotypes (n = 20) were further evaluated under field conditions to assess correlations between field and greenhouse data and to identify the most tolerant genotypes. A genome-wide association study (GWAS) was performed using ~203 K SNP markers to understand the genetic architecture of flooding tolerance in this panel. Survival scores between field and greenhouse data were significantly correlated (r = 0.55, P = 0.01). Subsequently, a subset of the most tolerant and susceptible genotypes were evaluated under pathogenic Pythium spp. pressure. This experiment revealed a potential link between flooding tolerance and Pythium spp. resistance. Several tolerant genotypes were identified that could be used as donor parents in breeding pipelines, especially ADP-429 and ADP-604. Based on the population structure analysis, a subpopulation consisting of 20 genotypes from the Middle American gene pool was detected that also possessed the highest root weight, hypocotyl length, and adventitious root development under flooding conditions. Genomic regions associated with flooding tolerance were identified including a region on Pv08/3.2 Mb, which is associated with germination rate and resides in vicinity of SnRK1.1, a central gene involved in response of plants to hypoxia. Furthermore, a QTL at Pv07/4.7 Mb was detected that controls survival score of seedlings under flooding conditions. The association of these QTL with the survivability traits including germination rate and survival score, indicates that these loci can be used in marker-assisted selection breeding to improve flooding tolerance in the Andean germplasm.

QTLs underlying the genetic interrelationship between efficient compatibility of Bradyrhizobium strains with soybean and genistein secretion by soybean roots

Soybean plants establish symbiotic relationships with soil rhizobia which form nodules on the plant roots. Nodule formation starts when the plant roots exudate isoflavonoids that induce nod gene expression of a specific Bradyrhizobium. We examined the specific indigenous rhizobia that form nodules with the soybean cultivars Peking and Tamahomare in different soils. PCR-RFLP analysis targeted to the 16S-23S rRNA gene internal transcribed spacer (ITS) region of the bacterial type of each root nodule showed that Bradyrhizobium japonicum (USDA110-type) and Bradyrhizobium elkanii (USDA94-type) had high compatibility with the Tamahomare and Peking cultivars, respectively. We grew 93 recombinant inbred lines (RIL) of soybean seeds derived from the cross between Peking and Tamahomare in three different field soils and identified the indigenous rhizobia nodulating each line using the same PCR-RFLP analysis. QTL analysis identified one QTL region in chromosome-18 with a highly significant additive effect that controls compatibility with both B. japonicum USDA110 and B. elkanii USDA94. We also measured the amount of daidzein and genistein secretion from roots of the 93 RILs by HPLC analysis. QTL analysis showed one QTL region in chromosome-18 controlling genistein secretion from roots and coinciding with that regulating compatibility of specific indigenous rhizobia with soybean. The amount of genistein may be a major regulatory factor in soybean-rhizobium compatibility.

Development of a high-density linkage map and chromosome segment substitution lines for Japanese soybean cultivar Enrei

Using progeny of a cross between Japanese soybean Enrei and Chinese soybean Peking, we developed a high-density linkage map and chromosomal segment substitution lines (CSSLs). The map consists of 2,177 markers with polymorphism information for 32 accessions and provides a detailed genetic framework for these markers. The marker order on the linkage map revealed close agreement with that on the chromosome-scale assembly, Wm82.a2.v1. The differences, especially on Chr. 5 and Chr. 11, in the present map provides information to identify regions in the genome assembly where additional information is required to resolve marker order and assign remaining scaffolds. To cover the entire soybean genome, we used 999 BC3F2 backcross plants and selected 103 CSSLs carrying chromosomal segments from Peking in the genetic background of Enrei. Using these low-genetic-complexity resources, we dissected variation in traits related to flowering, maturity and yield into approximately 50 reproducible quantitative trait loci (QTLs) and evaluated QTLs with small genetic effects as single genetic factors in a uniform genetic background. CSSLs developed in this study may be good starting material for removing the unfavourable characteristics of Peking during pre-breeding and for isolation of genes conferring disease and stress resistance that have not yet been characterized.

Comprehensive analysis of response and tolerant mechanisms in early-stage soybean at initial-flooding stress

Soybean is one of the most cultivated crops in the world; however, it is very sensitive to flooding stress, which markedly reduces its growth and yield. Morphological and biochemical changes such as an increase of fresh weight and a decrease of ATP content happen in early-stage soybean at initial-flooding stress, indicating that soybean responses to flooding stress are keys for its survival and seedling growth. Phosphoproteomics and nuclear proteomics are useful tools to detect protein-phosphorylation status and to identify transcriptional factors. In the review, the effect of flooding on soybean response to initial flooding stress is discussed based on recent results of proteomic, phosphoproteomic, nuclear proteomic, and nuclear phosphoproteomic studies. In addition, soybean survival under flooding stress, which is defined as tolerance mechanism, is discussed with the results of comprehensive analysis in flooding-tolerant mutant line and abscisic acid-treated soybean.

BIOLOGICAL SIGNIFICANCE

Soybean is one of the most cultivated crops in the world; however, it is very sensitive to flooding stress, especially soybean responses to initial flooding stress is key for its survival and seedling growth. Recently, proteomic techniques are applied to investigate the response and tolerant mechanisms of soybean at initial flooding condition. In this review, the progress in proteomic, phosphoproteomic, nuclear proteomic, and nuclear phosphoproteomic studies about the initial-flooding response mechanism in early-stage soybean is presented. In addition, the tolerant mechanism in soybean is discussed with the results of comprehensive analysis in flooding-tolerant mutant line and abscisic acid-treated soybean. Through this review, the key proteins and genes involved in initial flooding response and tolerance at early stage soybean are summarized and they contribute greatly to uncover response and tolerance mechanism at early stage under stressful environmental conditions in soybean.

Genetic Architecture of Flooding Tolerance in the Dry Bean Middle-American Diversity Panel

Flooding is a devastating abiotic stress that endangers crop production in the twenty-first century. Because of the severe susceptibility of common bean (Phaseolus vulgaris L.) to flooding, an understanding of the genetic architecture and physiological responses of this crop will set the stage for further improvement. However, challenging phenotyping methods hinder a large-scale genetic study of flooding tolerance in common bean and other economically important crops. A greenhouse phenotyping protocol was developed to evaluate the flooding conditions at early stages. The Middle-American diversity panel (n = 272) of common bean was developed to capture most of the diversity exits in North American germplasm. This panel was evaluated for seven traits under both flooded and non-flooded conditions at two early developmental stages. A subset of contrasting genotypes was further evaluated in the field to assess the relationship between greenhouse and field data under flooding condition. A genome-wide association study using ~150 K SNPs was performed to discover genomic regions associated with multiple physiological responses. The results indicate a significant strong correlation (r > 0.77) between greenhouse and field data, highlighting the reliability of greenhouse phenotyping method. Black and small red beans were the least affected by excess water at germination stage. At the seedling stage, pinto and great northern genotypes were the most tolerant. Root weight reduction due to flooding was greatest in pink and small red cultivars. Flooding reduced the chlorophyll content to the greatest extent in the navy bean cultivars compared with other market classes. Races of Durango/Jalisco and Mesoamerica were separated by both genotypic and phenotypic data indicating the potential effect of eco-geographical variations. Furthermore, several loci were identified that potentially represent the antagonistic pleiotropy. The GWAS analysis revealed peaks at Pv08/1.6 Mb and Pv02/41 Mb that are associated with root weight and germination rate, respectively. These regions are syntenic with two QTL reported in soybean (Glycine max L.) that contribute to flooding tolerance, suggesting a conserved evolutionary pathway involved in flooding tolerance for these related legumes.

Mapping quantitative trait loci for root development under hypoxia conditions in soybean (Glycine max L. Merr.)

KEY MESSAGE

Greatest potential, QTLs for hypoxia and waterlogging tolerance in soybean roots were detected using a new phenotypic evaluation method. Waterlogging is a major environmental stress limiting soybean yield in wet parts of the world. Root development is an important indicator of hypoxia tolerance in soybean. However, little is known about the genetic control of root development under hypoxia. This study was conducted to identify quantitative trait loci (QTLs) responsible for root development under hypoxia. Recombinant inbred lines (RILs) developed from a cross between a hypoxia-sensitive cultivar, Tachinagaha, and a tolerant landrace, Iyodaizu, were used. Seedlings were subjected to hypoxia, and root development was evaluated with the value change in root traits between after and before treatments. We found 230 polymorphic markers spanning 2519.2 cM distributed on all 20 chromosomes (Chrs.). Using these, we found 11 QTLs for root length (RL), root length development (RLD), root surface area (RSA), root surface area development (RSAD), root diameter (RD), and change in average root diameter (CARD) on Chrs. 11, 12, 13 and 14, and 7 QTLs for hypoxia tolerance of these root traits. These included QTLs for RLD and RSAD between markers Satt052 and Satt302 on Chr. 12, which are important markers of hypoxia tolerance in soybean; those QTLs were stable between 2 years. To validate the QTLs, we developed a near-isogenic line with the QTL region derived from Iyodaizu. The line performed well under both hypoxia and waterlogging, suggesting that the region contains one or more genes with large effects on root development. These findings may be useful for fine mapping and positional cloning of gene responsible for root development under hypoxia.

Identification of environmentally stable QTLs controlling Saponin content in Glycine max

Saponins are secondary metabolites that are widely distributed in plants. There are two major saponin precursors in soybean: soyasapogenol A, contributing to the undesirable taste, and soyasapogenol B, some of which have health benefits. It is important to control the ratio and content of the two major saponin groups to enhance the appeal of soybean as a health food. The structural diversity of saponin in the sugar chain composition makes it hard to quantify the saponin content. We measured the saponin content in soybean by removing the sugar chain from the saponin using acidic hydrolysis and detected novel quantitative trait loci (QTLs) for saponin content. Major QTLs in the hypocotyl were identified on chromosome 5 near the SSR marker, Satt 384, while those in the cotyledon were on chromosome 6 near Sat_312, which is linked to the T and E1 loci. Our results suggest that saponin contents in the hypocotyl and cotyledon are controlled by different genes and that it is difficult to increase the beneficial group B saponin and to decrease the undesirable group A saponin at the same time.

Comparison of Genetic Diversity between Chinese and American Soybean (Glycine max (L.)) Accessions Revealed by High-Density SNPs

Soybean is one of the most important economic crops for both China and the United States (US). The exchange of germplasm between these two countries has long been active. In order to investigate genetic relationships between Chinese and US soybean germplasm, 277 Chinese soybean accessions and 300 US soybean accessions from geographically diverse regions were analyzed using 5,361 SNP markers. The genetic diversity and the polymorphism information content (PIC) of the Chinese accessions was higher than that of the US accessions. Population structure analysis, principal component analysis, and cluster analysis all showed that the genetic basis of Chinese soybeans is distinct from that of the USA. The groupings observed in clustering analysis reflected the geographical origins of the accessions; this conclusion was validated with both genetic distance analysis and relative kinship analysis. F_ST-based and EigenGWAS statistical analysis revealed high genetic variation between the two subpopulations. Analysis of the 10 loci with the strongest selection signals showed that many loci were located in chromosome regions that have previously been identified as quantitative trait loci (QTL) associated with environmental-adaptation-related and yield-related traits. The pattern of diversity among the American and Chinese accessions should help breeders to select appropriate parental accessions to enhance the performance of future soybean cultivars.

Association mapping of soybean seed germination under salt stress

Soil salinity is a serious threat to agriculture sustainability worldwide. Seed germination is a critical phase that ensures the successful establishment and productivity of soybeans in saline soils. However, little information is available regarding soybean salt tolerance at the germination stage. The objective of this study was to identify the genetic mechanisms of soybean seed germination under salt stress. One natural population consisting of 191 soybean landraces was used in this study. Soybean seeds produced in four environments were used to evaluate the salt tolerance at their germination stage. Using 1142 single-nucleotide polymorphisms (SNPs), the molecular markers associated with salt tolerance were detected by genome-wide association analysis. Eight SNP-trait associations and 13 suggestive SNP-trait associations were identified using a mixed linear model and the TASSEL 4.0 software. Eight SNPs or suggestive SNPs were co-associated with two salt tolerance indices, namely (1) the ratio of the germination index under salt conditions to the germination index under no-salt conditions (ST-GI) and (2) the ratio of the germination rate under salt conditions to the germination rate under no-salt conditions (ST-GR). One SNP (BARC-021347-04042) was significantly associated with these two traits (ST-GI and ST-GR). In addition, nine possible candidate genes were located in or near the genetic region where the above markers were mapped. Of these, five genes, Glyma08g12400.1, Glyma08g09730.1, Glyma18g47140.1, Glyma09g00460.1, and Glyma09g00490.3, were verified in response to salt stress at the germination stage. The SNPs detected could facilitate a better understanding of the genetic basis of soybean salt tolerance at the germination stage, and the marker BARC-021347-04042 could contribute to future breeding for soybean salt tolerance by marker-assisted selection.

Core clock, SUB1, and ABAR genes mediate flooding and drought responses via alternative splicing in soybean

Circadian clocks are a great evolutionary innovation and provide competitive advantage during the day/night cycle and under changing environmental conditions. The circadian clock mediates expression of a large proportion of genes in plants, achieving a harmonious relationship between energy metabolism, photosynthesis, and biotic and abiotic stress responses. Here it is shown that multiple paralogues of clock genes are present in soybean (Glycine max) and mediate flooding and drought responses. Differential expression of many clock and SUB1 genes was found under flooding and drought conditions. Furthermore, natural variation in the amplitude and phase shifts in PRR7 and TOC1 genes was also discovered under drought and flooding conditions, respectively. PRR3 exhibited flooding- and drought-specific splicing patterns and may work in concert with PRR7 and TOC1 to achieve energy homeostasis under flooding and drought conditions. Higher expression of TOC1 also coincides with elevated levels of abscisic acid (ABA) and variation in glucose levels in the morning and afternoon, indicating that this response to abiotic stress is mediated by ABA, endogenous sugar levels, and the circadian clock to fine-tune photosynthesis and energy utilization under stress conditions. It is proposed that the presence of multiple clock gene paralogues with variation in DNA sequence, phase, and period could be used to screen exotic germplasm to find sources for drought and flooding tolerance. Furthermore, fine tuning of multiple clock gene paralogues (via a genetic engineering approach) should also facilitate the development of flooding- and drought-tolerant soybean varieties.

Identification of novel QTL governing root architectural traits in an interspecific soybean population

Cultivated soybean (Glycine max L.) cv. Dunbar (PI 552538) and wild G. soja (PI 326582A) exhibited significant differences in root architecture and root-related traits. In this study, phenotypic variability for root traits among 251 BC2F5 backcross inbred lines (BILs) developed from the cross Dunbar/PI 326582A were identified. The root systems of the parents and BILs were evaluated in controlled environmental conditions using a cone system at seedling stage. The G. max parent Dunbar contributed phenotypically favorable alleles at a major quantitative trait locus on chromosome 8 (Satt315-I locus) that governed root traits (tap root length and lateral root number) and shoot length. This QTL accounted for >10% of the phenotypic variation of both tap root and shoot length. This QTL region was found to control various shoot- and root-related traits across soybean genetic backgrounds. Within the confidence interval of this region, eleven transcription factors (TFs) were identified. Based on RNA sequencing and Affymetrix expression data, key TFs including MYB, AP2-EREBP and bZIP TFs were identified in this QTL interval with high expression in roots and nodules. The backcross inbred lines with different parental allelic combination showed different expression pattern for six transcription factors selected based on their expression pattern in root tissues. It appears that the marker interval Satt315-I locus on chromosome 8 contain an essential QTL contributing to early root and shoot growth in soybean.

Analyses of flooding tolerance of soybean varieties at emergence and varietal differences in their proteomes

Flooding of fields due to heavy and/or continuous rainfall influences soybean production. To identify soybean varieties with flooding tolerance at the seedling emergence stage, 128 soybean varieties were evaluated using a flooding tolerance index, which is based on plant survival rates, the lack of apparent damage and lateral root development, and post-flooding radicle elongation rate. The soybean varieties were ranked according to their flooding tolerance index, and it was found that the tolerance levels of soybean varieties exhibit a continuum of differences between varieties. Subsequently, tolerant, moderately tolerant and sensitive varieties were selected and subjected to comparative proteomic analysis to clarify the tolerance mechanism. Proteomic analysis of the radicles, combined with correlation analysis, showed that the ratios of RNA binding/processing related proteins and flooding stress indicator proteins were significantly correlated with flooding tolerance index. The RNA binding/processing related proteins were positively correlated in untreated soybeans, whereas flooding stress indicator proteins were negatively correlated in flooded soybeans. These results suggest that flooding tolerance is regulated by mechanisms through multiple factors and is associated with abundance levels of the identified proteins.

Differential expression and elution behavior of basic 7S globulin among cultivars under hot water treatment of soybean seeds

Basic 7S globulin (Bg7S), which accumulates in mature soybean (Glycine max) seeds, is an extracellular matrix protein. A large amount of Bg7S is synthesized de novo and is eluted from soybean seeds when immersed in 50-60°C water (hot water treatment, HWT). However, the Bg7S elution mechanism remains unclear. Under HWT, the seeds probably undergo heat stress and flooding stress. To obtain fundamental knowledge related to how Bg7S is eluted from hot-water-treated seeds, this study compared Bg7S elution among soybean cultivars having different flooding tolerance during pre-germination. The amounts of Bg7S eluted from seeds varied significantly among cultivars. Elution was suppressed by seed coats regarded as preventing the leakage of seed contents by rapid water imbibition. Furthermore, Bg7S expression levels differed among cultivars, although the difference did not result from any variation in Bg7S promoter sequences. However, the expression levels of Bg7S under HWT were not associated with the flooding tolerance level. Immunoelectron microscopy revealed that the Bg7S accumulated in the intercellular space of hot-water-treated seeds. Plasma membrane shrinkage was observed. The main proteins eluted from seeds under HWT were located in the extracellular space. This study clarified the mechanism of Bg7S elution from seeds under HWT.

Single- and double-SSR primer combined analyses in rice

Polymerase chain reaction (PCR) is the foundation of SSR molecular marker technology. We used sib rice varieties J518, XD1 and SD23 as experimental materials, selecting 30 pairs of SSR primers, including RM127, RM337 and RM5172, covering the rice genome, and performed single- and double-SSR primer combined analyses. We found that under the same PCR system and conditions, a single primer of the SSR primer pairs could amplify the same fragments as double primers do. The sequencing results demonstrated that some amplified fragments that we previously believed to come from double primers were actually produced by a single primer. The use of this kind of primer, such as the RM127 primer pair, for marker-assisted breeding will therefore be misleading. Additionally, using the same PCR system and conditions, some single primers that are part of SSR primer pairs can amplify many more specific fragments than double-SSR primers. For instance, in the case of the RM5172 primer pair, a single primer P1 amplified approximately three times the number of fragments as the double primer. This information can contribute to research on genetic diversity of species, understanding of genetic relationships and identification of germplasm resources. Accordingly, combined analyses of single- and double-primer amplification products not only can remove single-primer amplification fragments and false-positives from double-primer amplification products in order to improve test accuracy, but also can facilitate research on genetic diversity, exploration of phylogenetic relationships and identification of germplasm resources. We define this method as "single- and double-SSR primer combined analyses".

Problems with and a system to eliminate single-primer PCR product contamination in simple sequence repeat molecular marker-assisted selection in soybean

Polymerase chain reaction (PCR) provides a foundation for simple sequence repeat molecular marker-assisted selection (SSR MAS) in soybean. This PCR system and its various conditions have been optimized by many researchers. However, current research on the optimization of the PCR system focuses on double-primer PCR products. We compared single- and double-SSR primer PCR products from 50 soybean samples and found that the use of single-PCR primers in the reaction system can lead to amplified fragments of portions of the SSR primers in the PCR process, resulting in both false-positives and fragment impurity of double-primer PCR amplification, inconvenient for subsequent analysis. We used "single-primer PCR correction" to eliminate interference caused by single-primer nonspecific PCR amplification and improve PCR quality. Using this method, the precision and success rates of SSR MAS in soybean can be increased.

Functional genomics of soybean for improvement of productivity in adverse conditions

Global soybean production is frequently impacted by various stresses, including both abiotic and biotic stresses. To develop soybean plants with enhanced tolerance to different stressors, functional genomics of soybean and a comprehensive understanding of available biotechnological resources and approaches are essential. In this review, we will discuss recent advances in soybean functional genomics which provide unprecedented opportunities to understand global patterns of gene expression, gene regulatory networks, various physiological, biochemical, and metabolic pathways as well as their association with the development of specific phenotypes. Soybean functional genomics, therefore, will ultimately enable us to develop new soybean varieties with improved productivity under adverse conditions by genetic engineering.