Identification and Comparative Analysis of Differential Gene Expression in Soybean Leaf Tissue under Drought and Flooding Stress Revealed by RNA-Seq

Characterization of fatty acid desaturase gene family in Glycine max and their expression patterns in seeds after Fusarium fujikuroi infection

Background

The family of membrane-bound fatty acid desaturase (FAD) genes play a vital role in plant growth, development, and stress responses. The seed-borne pathogen Fusarium fujikuroi causes seed decay disease during pre-harvest and post-harvest stages of soybean, leading to a significant reduction in yield and quality. Therefore, it is very meaningful to characterize the diversity and function of the GmFAD gene family in soybean and to elucidate their roles in seed resistance to F. fujikuroi.

Results

In this study, 30 full-length GmFAD genes were identified from the soybean genome. A range of analysis was conducted to characterize gene and protein structures, chromosomal locations, conserved motif and conserved structural domains, and results showed that GmFAD genes were clustered into seven subfamilies (FAB2, ADS, SLD, DES, FAD6, FAD2, FAD3/7/8), which is also supported by phylogenetic analysis. The diversity and expansion of the GmFAD gene family were mainly caused by segmental duplication, and their encoding proteins were observed to locate in chloroplast or endoplasmic reticulum. The promoters of GmFAD genes contained a set of cis-acting elements in response to plant hormone, defense and stress, light, and plant growth and development, indicating these genes have the complex expression regulation and diverse functions. Gene ontology (GO) and KEGG enrichment pathway analyses showed that GmFAD genes were closely related to the biosynthesis and metabolism of lipid and unsaturated fatty acids (UFAs). In addition, the expression of GmFADs was significantly changed in soybean seeds when challenged by the seed decay pathogen F. fujikuroi. Specifically, GmFAB2.1/2.2, GmFAD3.3/3-2B/7-1//8-2, and GmFAD2.3/2.5 genes displayed distinct temporal expression patterns in the resistant ND25 and susceptible CX12, highlighting their potential roles in soybean resistance against F. fujikuroi infection.

Conclusion

Our findings contribute to a deeper understanding of the GmFAD gene family and their intricate roles in soybean resistance against the seed-borne pathogen F. fujikuroi. Moreover, several distinct genes provide valuable candidates for further application in soybean resistant breeding.

Enhancing soybean germination and vigor under water stress: the efficacy of bio-priming with sodium carboxymethyl cellulose and gum arabic

Seed priming can significantly enhance the tolerance of soybean against different environmental stresses by improving seed water uptake and modulating stress-response mechanisms. In particular, seed priming with sodium carboxymethylcellulose (SCMC) and gum Arabic (GA) can support seeds to withstand extreme conditions better, promoting more consistent germination and robust seedling establishment, which is crucial for achieving stable agricultural yields. The present study investigated the effects of seed priming using a combination of SCMC and GA (10% CG) on the germination, growth, and biochemical responses of six soybean varieties under drought and flooding stress conditions. The results revealed significant differences among varieties and applied treatments on germination, vigor, and physiological traits. Under drought stress, seed priming with 10% CG significantly improved germination percentage, germination rate, shoot length, root length, and biomass compared to unprimed seeds. Notable reductions in malondialdehyde (MDA) content and enhanced antioxidant enzyme activities, including superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD), suggest that 10% CG priming mitigates oxidative damage through enhanced antioxidant defense mechanisms. Moreover, 10% CG seed priming improved germination and growth parameters under flooding stress, but the advantages were less significant. In addition, the priming treatment significantly reduced electrolyte conductivity (EC) across all varieties compared to unprimed seeds, indicating improved membrane stability. Overall, 10% CG seed priming was more effective under drought and flooding conditions, demonstrating a potential strategy for enhancing stress tolerance in soybean varieties.

Could flooding undermine progress in building climate-resilient crops?

Flooding threatens crop productivity, agricultural sustainability, and global food security. In this article I review the effects of flooding on plants and highlight three important gaps in our understanding: (i) effects of flooding on ecological interactions mediated by plants both below (changing root metabolites and exudates) and aboveground (changing plant quality and metabolites, and weakening the plant immune system), (ii) flooding impacts on soil health and microorganisms that underpin plant and ecosystems health, and (iii) the legacy impacts of flooding. Failure to address these overlooked aspects could derail and undermine the monumental progress made in building climate-resilient crops and soil-microbe-assisted plant resilience. Addressing the outlined knowledge gaps will enhance solutions developed to mitigate flooding and preserve gains made to date.

Metabolomics combined with transcriptomics and physiology reveals the regulatory responses of soybean plants to drought stress

Drought, a prevalent environmental stressor, has had significant consequences on soybean (Glycine max L.), notably impeding its growth and production. Therefore, it is crucial to gain insight into the regulatory responses of soybean plants exposed to drought stress during soybean flowering in the field. In this study, the cultivar 'Liaodou 15' was performed light drought (LD, 24.3% soil moisture content), moderate drought (MD, 20.6% soil moisture content) and severe drought (SD, 16.9% soil moisture content) treatments at flowering stages of soybean and then rehydrated (30% soil moisture content) until harvest. The yield-related indicators were measured and revealed that MD and SD treatments significantly reduced 6.3% and 10.8% of the 100-grain weight. Soybean plants subjected to three drought stresses showed that net photosynthetic rates were 20.8%, 51.5% and 71.8% lower in LD, MD and SD than that of CK. The WUE increased by 31.8%, 31.5% and 18.8% under three drought stress treatments compared to CK. In addition, proline content was 25.94%, 41.01% and 65.43% greater than that of CK under three drought stress treatments. The trend of the MDA content was consistent with that of the proline content. SOD activity was significantly increasing by 10.86%, 46.73% and 14.54% under three drought stress treatments. The activity of CAT in the SD treatment increased by 49.28%. All the indices recovered after rehydration. Furthermore, 54,78 and 51 different expressed metabolomics (DEMs) were identified in the LDCK/LD, MDCK/MD and SDCK/SD groups, respectively. There were 1,211, 1,265 and 1,288 different expressed genes (DEGs) were upregulated and 1,003, 1,819 and 1,747 DEGs were downregulated. Finally, combined transcriptomic and metabolomic analysis suggested that 437 DEGs and 24 DEMs of LDCK/LD group, 741 DEGs and 35 DEMs of MDCK/MD group, 633 DEGs and 23 DEMs of SDCK/SD group, were highly positively correlated in soybean plants under drought stress. Drought stress induced the expression of the PAO1, PAO4, PAO5 and P5CS genes to promote the accumulation of spermidine and proline. Our study elucidates the responses of drought-stressed soybean plants in the field and provides a genetic basis for the breeding of drought-tolerant soybean plants.

Transcriptome analysis and characteristics of drought resistance related genes in four varieties of foxtail millet [Setaria italica]

Foxtail millet [Setaria italica] plays a crucial role as a multigrain crop in agricultural production. However, due to future extreme weather conditions, drought remains the main abiotic stress that limits foxtail millet yield, it is highly significant to screen for drought-tolerant varieties throughout the entire growth period and identify the regulatory genes associated with drought resistance in foxtail millet breeding. We identified 217 foxtail millet seed resources for drought resistance during the maturity stage in the field, and subsequently categorized them into different levels of drought resistance. Two cultivars with extremely strong drought resistance during the maturity stage in the field, JKH4 (Chi 5422) and JKH6 (Chigu 26), as well as two cultivars with extremely weak drought resistance during the maturity stage in the field, JRK3 (17M1309) and JRK6 (Canggu 9), were selected for physiological comparison and transcriptome sequencing before and after drought treatment. Transcriptome analysis at the seedling stage revealed that JRK3 and JRK6 cultivar primarily regulated phenylpropanoid biosynthesis, MAPK signaling pathogen-plant, and plant hormone signal transduction pathway in response to drought stress. On the other hand, the fatty acid elongation pathway of JKH4 and JKH6 variety was found to be more significant. Furthermore, 22 drought resistance related genes were screened through transcriptome analysis of four foxtail millet varieties. These findings could offer valuable theoretical guidance for breeding foxtail millet with enhanced drought resistance and potentially facilitate the development of genetically engineered drought-resistant foxtail millet varieties.

Assessment of waterlogging-induced changes in enzymatic antioxidants and carbohydrate metabolism in peanuts genotypes

Soil flooding, manifesting as submergence or waterlogging stress, significantly impacts plant species composition and agricultural productivity, particularly in regions with low rainfall. This study investigates the biochemical responses of two peanut (Arachis hypogaea L.) genotypes, DH-86 and GJG-32, under waterlogging stress. The experiment involved in-vivo pot trials where peanut plants were subjected to continuous waterlogging for 12 days at the flowering stage. Biochemical analyses of leaves conducted and revealed significant alterations in enzyme activities and metabolite concentrations. Key findings include variations in superoxide dismutase (SOD), catalase (CAT), guaiacol peroxidase (GPOD), α-amylase, invertase, acid phosphomonoesterase activities, and changes in starch, proline, reducing sugars, and chlorophyll content. SOD, CAT, and GPOD activities exhibited differential responses between genotypes, highlighting DH-86's quicker recovery post-waterlogging. Notably, DH-86 demonstrated higher resilience, reflected in its rapid normalization of biochemical parameters, while GJG-32 showed prolonged stress effects. These findings underscore the importance of antioxidative enzyme systems in mitigating oxidative damage induced by waterlogging. This study enhances our understanding of the biochemical adaptations of peanut genotypes to waterlogging stress, offering valuable insights for breeding programs focused on improving flood tolerance in crops.

Identification of Drought-Resistant Response in Proso Millet (Panicum miliaceum L.) Root through Physiological and Transcriptomic Analysis

Proso millet (Panicum miliaceum L.) is resilient to abiotic stress, especially to drought. However, the mechanisms by which its roots adapt and tolerate salt stress are obscure. In this study, to clarify the molecular mechanism of proso millet in response to drought stress, the physiological indexes and transcriptome in the root of seedlings of the proso millet cultivar 'Yumi 2' were analyzed at 0, 0.5, 1.0, 1.5, and 3.0 h of stimulated drought stress by using 20% PEG-6000 and after 24 h of rehydration. The results showed that the SOD activity, POD activity, soluble protein content, MDA, and O2-· content of 'Yumi 2' increased with the time of drought stress, but rapidly decreased after rehydration. Here, 130.46 Gb of clean data from 18 samples were obtained, and the Q30 value of each sample exceeded 92%. Compared with 0 h, the number of differentially expressed genes (DEGs) reached the maximum of 16,105 after 3 h of drought, including 9153 upregulated DEGs and 6952 downregulated DEGs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed that upregulated DEGs were mainly involved in ATP binding, nucleus, protein serine/threonine phosphatase activity, MAPK signaling pathway-plant, plant-pathogen interactions, and plant hormone signal transduction under drought stress, while downregulated DEGs were mainly involved in metal ion binding, transmembrane transporter activity, and phenylpropanoid biosynthesis. Additionally, 1441 TFs screened from DEGs were clustered into 64 TF families, such as AP2/ERF-ERF, bHLH, WRKY, NAC, MYB, and bZIP TF families. Genes related to physiological traits were closely related to starch and sucrose metabolism, phenylpropanoid biosynthesis, glutathione metabolism, and plant hormone signal transduction. In conclusion, the active oxygen metabolism system and the soluble protein of proso millet root could be regulated by the activity of protein serine/threonine phosphatase. AP2/ERF-ERF, bHLH, WRKY, NAC, MYB, and bZIP TF families were found to be closely associated with drought tolerance in proso millet root. This study will provide data to support a subsequent study on the function of the drought tolerance gene in proso millet.

Integrated meta-analysis and transcriptomics pinpoint genomic loci and novel candidate genes associated with submergence tolerance in rice

BACKGROUND

Due to rising costs, water shortages, and labour shortages, farmers across the globe now prefer a direct seeding approach. However, submergence stress remains a major bottleneck limiting the success of this approach in rice cultivation. The merger of accumulated rice genetic resources provides an opportunity to detect key genomic loci and candidate genes that influence the flooding tolerance of rice.

RESULTS

In the present study, a whole-genome meta-analysis was conducted on 120 quantitative trait loci (QTL) obtained from 16 independent QTL studies reported from 2004 to 2023. These QTL were confined to 18 meta-QTL (MQTL), and ten MQTL were successfully validated by independent genome-wide association studies from diverse natural populations. The mean confidence interval (CI) of the identified MQTL was 3.44 times narrower than the mean CI of the initial QTL. Moreover, four core MQTL loci with genetic distance less than 2 cM were obtained. By combining differentially expressed genes (DEG) from two transcriptome datasets with 858 candidate genes identified in the core MQTL regions, we found 38 common differentially expressed candidate genes (DECGs). In silico expression analysis of these DECGs led to the identification of 21 genes with high expression in embryo and coleoptile under submerged conditions. These DECGs encode proteins with known functions involved in submergence tolerance including WRKY, F-box, zinc fingers, glycosyltransferase, protein kinase, cytochrome P450, PP2C, hypoxia-responsive family, and DUF domain. By haplotype analysis, the 21 DECGs demonstrated distinct genetic differentiation and substantial genetic distance mainly between indica and japonica subspecies. Further, the MQTL7.1 was successfully validated using flanked marker S2329 on a set of genotypes with phenotypic variation.

CONCLUSION

This study provides a new perspective on understanding the genetic basis of submergence tolerance in rice. The identified MQTL and novel candidate genes lay the foundation for marker-assisted breeding/engineering of flooding-tolerant cultivars conducive to direct seeding.

Physiological and transcriptomic analyses of leaves from Gardenia jasminoides Ellis under waterlogging stress

Abstract Gardenia jasminoides Ellis is a Chinese herbal medicine with medicinal and economic value, but its mechanism of response to waterlogging stress remains unclear. In this study, the “double pots method” was used to simulate the waterlogging stress of Gardenia jasminoides Ellis to explore its physiological and transcriptomic response mechanism. We found no significant damage to Gardenia jasminoides Ellis membrane lipid during stress. POD played a vital antioxidant role, KEGG enrichment showed that secondary metabolites such as flavonoids might also play an antioxidant role, and PRO played a significant osmotic adjustment. Endogenous hormones regulate the Gardenia jasminoides Ellis's growth and development and play a role in signal transduction. Among them, light waterlogging stress is delayed. At the same time, there were 19631, 23693, and 15045 differentially expressed genes on the 5th, 10d, and 15d of Gardenia jasminoides Ellis under waterlogging stress. These genes were closely associated with the proteasome, endopeptidase, ribosome, MAPK signal transduction, and endogenous hormone signal transduction, plant-pathogen interaction and phenylpropanoid biosynthesis and other physiological and metabolic pathways, which regulate the turnover and transportation of protein, the reinforcement and adhesion of cell walls, the induction of stomatal closure, allergic reactions, defense reactions, leaf movements and others. It also can absorb ultraviolet rays to reduce the generation of oxygen free radicals, change the way of energy utilization and adjust the osmotic pressure of plant cells.

Biochemical Analysis to Understand the Flooding Tolerance of Mutant Soybean Irradiated with Gamma Rays

Flooding stress, which reduces plant growth and seed yield, is a serious problem for soybean. To improve the productivity of flooded soybean, flooding-tolerant soybean was produced by gamma-ray irradiation. Three-day-old wild-type and mutant-line plants were flooded for 2 days. Protein, RNA, and genomic DNA were then analyzed based on oppositely changed proteins between the wild type and the mutant line under flooding stress. They were associated with cell organization, RNA metabolism, and protein degradation according to proteomic analysis. Immunoblot analysis confirmed that the accumulation of beta-tubulin/beta-actin increased in the wild type under flooding stress and recovered to the control level in the mutant line; however, alpha-tubulin increased in both the wild type and the mutant line under stress. Ubiquitin was accumulated and genomic DNA was degraded by flooding stress in the wild type; however, they were almost the same as control levels in the mutant line. On the other hand, the gene expression level of RNase H and 60S ribosomal protein did not change in either the wild type or the mutant line under flooding stress. Furthermore, chlorophyll a/b decreased and increased in the wild type and the mutant line, respectively, under flooding stress. These results suggest that the regulation of cell organization and protein degradation might be an important factor in the acquisition of flooding tolerance in soybean.

GmWRKY17-mediated transcriptional regulation of GmDREB1D and GmABA2 controls drought tolerance in soybean

Drought affects soybean growth and ultimately led to yield reduction. WRKY transcription factors involve in the regulation of abiotic stress. Few functions of WRKY transcription factors underlying drought tolerance in soybean are clear. Here, we reported a WRKY transcription factor named GmWRKY17 that positively regulates soybean drought tolerance by regulating drought-induced genes and ABA-related genes. Transcriptome sequencing (RNA-Seq) and yeast one hybrid analysis identified downstream genes regulated by GmWRKY17. ChIP-qPCR, EMSA and dual-luciferase reporter assay showed that GmWRKY17 directly bound to the promoters of the GmDREB1D and GmABA2, and activated their expression under drought stress. Overexpression of GmDREB1D gene enhanced drought tolerance of soybean. Taken together, our study revealed a regulatory mechanism that GmWRKY17 transcription factor may improve soybean drought tolerance by mediating ABA synthesis and DREB signaling pathway.

Genome-wide association mapping for yield-related traits in soybean (Glycine max) under well-watered and drought-stressed conditions

Soybean (Glycine max) productivity is significantly reduced by drought stress. Breeders are aiming to improve soybean grain yields both under well-watered (WW) and drought-stressed (DS) conditions, however, little is known about the genetic architecture of yield-related traits. Here, a panel of 188 soybean germplasm was used in a genome wide association study (GWAS) to identify single nucleotide polymorphism (SNP) markers linked to yield-related traits including pod number per plant (PN), biomass per plant (BM) and seed weight per plant (SW). The SLAF-seq genotyping was conducted on the population and three phenotype traits were examined in WW and DS conditions in four environments. Based on best linear unbiased prediction (BLUP) data and individual environmental analyses, 39 SNPs were significantly associated with three soybean traits under two conditions, which were tagged to 26 genomic regions by linkage disequilibrium (LD) analysis. Of these, six QTLs qPN-WW19.1, qPN-DS8.8, qBM-WW1, qBM-DS17.4, qSW-WW4 and qSW-DS8 were identified controlling PN, BM and SW of soybean. There were larger proportions of favorable haplotypes for locus qPN-WW19.1 and qSW-WW4 rather than qBM-WW1, qBM-DS17.4, qPN-DS8.8 and qSW-DS8 in both landraces and improved cultivars. In addition, several putative candidate genes such as Glyma.19G211300, Glyma.17G057100 and Glyma.04G124800, encoding E3 ubiquitin-protein ligase BAH1, WRKY transcription factor 11 and protein zinc induced facilitator-like 1, respectively, were predicted. We propose that the further exploration of these locus will facilitate accelerating breeding for high-yield soybean cultivars.

The relative impact of parental and current environment on plant transcriptomes depends on type of stress and genotype

Through developmental plasticity, an individual organism integrates influences from its immediate environment with those due to the environment of its parents. While both effects on phenotypes are well documented, their relative impact has been little studied in natural systems, especially at the level of gene expression. We examined this issue in four genotypes of the annual plant Persicaria maculosa by varying two key resources-light and soil moisture-in both generations. Transcriptomic analyses showed that the relative effects of parent and offspring environment on gene expression (i.e. the number of differentially expressed transcripts, DETs) varied both for the two types of resource stress and among genotypes. For light, immediate environment induced more DETs than parental environment for all genotypes, although the precise proportion of parental versus immediate DETs varied among genotypes. By contrast, the relative effect of soil moisture varied dramatically among genotypes, from 8-fold more DETs due to parental than immediate conditions to 10-fold fewer. These findings provide evidence at the transcriptomic level that the relative impacts of parental and immediate environment on the developing organism may depend on the environmental factor and vary strongly among genotypes, providing potential for the interplay of these developmental influences to evolve.

Phosphate starvation response precedes abscisic acid response under progressive mild drought in plants

Drought severely damages crop production, even under conditions so mild that the leaves show no signs of wilting. However, it is unclear how field-grown plants respond to mild drought. Here, we show through six years of field trials that ridges are a useful experimental tool to mimic mild drought stress in the field. Mild drought reduces inorganic phosphate levels in the leaves to activate the phosphate starvation response (PSR) in soybean plants in the field. Using Arabidopsis thaliana and its mutant plants grown in pots under controlled environments, we demonstrate that PSR occurs before abscisic acid response under progressive mild drought and that PSR plays a crucial role in plant growth under mild drought. Our observations in the field and laboratory using model crop and experimental plants provide insight into the molecular response to mild drought in field-grown plants and the relationship between nutrition and drought stress response.

Elucidating the molecular responses to waterlogging stress in onion (Allium cepa L.) leaf by comparative transcriptome profiling

Introduction

Waterlogging is a major stress that severely affects onion cultivation worldwide, and developing stress-tolerant varieties could be a valuable measure for overcoming its adverse effects. Gathering information regarding the molecular mechanisms and gene expression patterns of waterlogging-tolerant and sensitive genotypes is an effective method for improving stress tolerance in onions. To date, the waterlogging tolerance-governing molecular mechanism in onions is unknown.

Methods

This study identified the differentially expressed genes (DEGs) through transcriptome analysis in leaf tissue of two onion genotypes (Acc. 1666; tolerant and W-344; sensitive) presenting contrasting responses to waterlogging stress.

Results

Differential gene expression analysis revealed that in Acc. 1666, 1629 and 3271 genes were upregulated and downregulated, respectively. In W-344, 2134 and 1909 genes were upregulated and downregulated, respectively, under waterlogging stress. The proteins coded by these DEGs regulate several key biological processes to overcome waterlogging stress such as phytohormone production, antioxidant enzymes, programmed cell death, and energy production. The clusters of orthologous group pathway analysis revealed that DEGs contributed to the post-translational modification, energy production, and carbohydrate metabolism-related pathways under waterlogging stress. The enzyme assay demonstrated higher activity of antioxidant enzymes in Acc. 1666 than in W-344. The differential expression of waterlogging tolerance related genes, such as those related to antioxidant enzymes, phytohormone biosynthesis, carbohydrate metabolism, and transcriptional factors, suggested that significant fine reprogramming of gene expression occurs in response to waterlogging stress in onion. A few genes such as ADH, PDC, PEP carboxylase, WRKY22, and Respiratory burst oxidase D were exclusively upregulated in Acc. 1666.

Discussion

The molecular information about DEGs identified in the present study would be valuable for improving stress tolerance and for developing waterlogging tolerant onion varieties.

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.

Physiological and Transcriptome Analyses of Photosynthesis in Three Mulberry Cultivars within Two Propagation Methods (Cutting and Grafting) under Waterlogging Stress

Mulberry is a valuable woody plant with significant economic importance. It can be propagated through two main methods: cutting and grafting. Waterlogging can have a major impact on mulberry growth and can significantly reduce production. In this study, we examined gene expression patterns and photosynthetic responses in three waterlogged mulberry cultivars propagated through cutting and grafting. Compared to the control group, waterlogging treatments reduced levels of chlorophyll, soluble protein, soluble sugars, proline, and malondialdehyde (MDA). Additionally, the treatments significantly decreased the activities of ascorbate peroxidase (APX), peroxidase (POD), and catalase (CAT) in all three cultivars, except for superoxide dismutase (SOD). Waterlogging treatments also affected the rate of photosynthesis (Pn), stomatal conductance (Gs), and transpiration rate (Tr) in all three cultivars. However, no significant difference in physiological response was observed between the cutting and grafting groups. Gene expression patterns in the mulberry changed dramatically after waterlogging stress and varied between the two propagation methods. A total of 10,394 genes showed significant changes in expression levels, with the number of differentially expressed genes (DEGs) varying between comparison groups. GO and KEGG analysis revealed important DEGs, including photosynthesis-related genes that were significantly downregulated after waterlogging treatment. Notably, these genes were upregulated at day 10 in the cutting group compared to the grafting group. In particular, genes involved in carbon fixation were significantly upregulated in the cutting group. Finally, cutting propagation methods displayed better recovery capacity from waterlogging stress than grafting. This study provides valuable information for improving mulberry genetics in breeding programs.

Heterologous expression of Zygophyllum xanthoxylon zinc finger protein gene (ZxZF) enhances the tolerance of poplar photosynthetic function to drought stress

The ZxZF transcription factor (TF) of Zygophyllum xanthoxylon (Bunge) Maxim, an extremely drought-resistant woody plant, is a C₂H₂ zinc finger protein. Studies have shown that C₂H₂ zinc finger proteins play important roles in activating stress-related genes and enhancing plant resistance. However, their function in regulating plant photosynthesis under drought stress is not well understood. Since poplar is an important greening and afforestation tree species, it is particularly important to cultivate excellent drought-tolerant varieties. The ZxZF transcription factor (TF) was heterogeneously expressed in Euroamerican poplar (Populus × euroameracana cl.'Bofengl') by genetic transformation. Based on the mechanism and potential function of poplar photosynthesis regulated by ZxZF under drought stress, transcriptomic and physiological determinations were used to reveal the important role of this gene in improving the drought resistance of poplar. The results showed that the overexpression of ZxZF TF in transgenic poplars could improve the inhibition of Calvin cycle by regulating stomatal opening and increasing the concentration of intercellular CO₂. The chlorophyll content, photosynthetic performance index, and photochemical efficiency of transgenic lines under drought stress were significantly higher than those of the wild type (WT). The overexpression of ZxZF TFs could alleviate the degree of photoinhibition of photosystems II and I under drought stress and maintain the efficiency of light energy capture and the photosynthetic electron transport chain. The transcriptomic data also showed that differentially expressed genes between the transgenic poplar and WT under drought stress were primarily enriched in metabolic pathways related to photosynthesis, such as photosynthesis, photosynthesis-antenna protein, porphyrin and chlorophyll metabolism, and photosynthetic carbon fixation, and the downregulation of genes related to chlorophyll synthesis, photosynthetic electron transport and Calvin cycle were alleviated. In addition, the overexpression of ZxZF TF can alleviate the inhibition of NADH dehydrogenase-like (NDH) cyclic electron flow of the poplar NDH pathway under drought stress, which plays an important role in reducing the excess pressure of electrons on the photosynthetic electron transport chain and maintaining the normal photosynthetic electron transport. In summary, the overexpression of ZxZF TFs can effectively alleviate the inhibition of drought on the assimilation of carbon in poplar and have a positive impact on light energy capture, the orderly transport of photosynthetic electron transport chain and the integrity of the photosystem, which is highly significant to acheivean in-depth understanding of the function of ZxZF TFs. This also provides an important basis for the breeding of new transgenic poplar varieties.

How Central Carbon Metabolites of Mexican Mint (Plectranthus amboinicus) Plants Are Impacted under Different Watering Regimes

Plants are sessile, and their ability to reprogram their metabolism to adapt to fluctuations in soil water level is crucial but not clearly understood. A study was performed to determine alterations in intermediate metabolites involved in central carbon metabolism (CCM) following exposure of Mexican mint (Plectranthus amboinicus) to varying watering regimes. The water treatments were regular watering (RW), drought (DR), flooding (FL), and resumption of regular watering after flooding (DHFL) or after drought (RH). Leaf cluster formation and leaf greening were swift following the resumption of regular watering. A total of 68 key metabolites from the CCM routes were found to be significantly (p < 0.01) impacted by water stress. Calvin cycle metabolites in FL plants, glycolytic metabolites in DR plants, total tricarboxylic acid (TCA) cycle metabolites in DR and DHFL plants, and nucleotide biosynthetic molecules in FL and RH plants were significantly (p < 0.05) increased. Pentose phosphate pathway (PPP) metabolites were equally high in all the plants except DR plants. Total Calvin cycle metabolites had a significantly (p < 0.001) strong positive association with TCA cycle (r = 0.81) and PPP (r = 0.75) metabolites. Total PPP metabolites had a moderately positive association with total TCA cycle metabolites (r = 0.68; p < 0.01) and a negative correlation with total glycolytic metabolites (r = -0.70; p < 0.005). In conclusion, the metabolic alterations of Mexican mint plants under different watering regimes were revealed. Future studies will use transcriptomic and proteomic approaches to identify genes and proteins that regulate the CCM route.

Genome-Wide Identification and Analysis of the Maize Serine Peptidase S8 Family Genes in Response to Drought at Seedling Stage

Subtilisin-like proteases (subtilases) are found in almost all plant species and are involved in regulating various biotic and abiotic stresses. Although the literature on subtilases in different plant species is vast, the gene function of the serine peptidase S8 family and its maize subfamily is still unknown. Here, a bioinformatics analysis of this gene family was conducted by describing gene structure, conserved motifs, phylogenetic relationships, chromosomal distributions, gene duplications, and promoter cis-elements. In total, we identified 18 ZmSPS8 genes in maize, distributed on 7 chromosomes, and half of them were hydrophilic. Most of these proteins were located at the cell wall and had similar secondary and tertiary structures. Prediction of cis-regulatory elements in promoters illustrated that they were mainly associated with hormones and abiotic stress. Maize inbred lines B73, Zheng58, and Qi319 were used to analyze the spatial-temporal expression patterns of ZmSPS8 genes under drought treatment. Seedling drought results showed that Qi319 had the highest percent survival after 14 d of withholding irrigation, while B73 was the lowest. Leaf relative water content (LRWC) declined more rapidly in B73 and to lower values, and the nitrotetrazolium blue chloride (NBT) contents of leaves were higher in Qi319 than in the other inbreds. The qPCR results indicated that 6 serine peptidase S8 family genes were positively or negatively correlated with plant tolerance to drought stress. Our study provides a detailed analysis of the ZmSPS8s in the maize genome and finds a link between drought tolerance and the family gene expression, which was established by using different maize inbred lines.

Physiological and comparative transcriptome analysis of the response and adaptation mechanism of the photosynthetic function of mulberry (Morus alba L.) leaves to flooding stress

Flooding has become one of the major abiotic stresses that seriously affects plant growth and development owing to changes in the global precipitation pattern. Mulberry (Morus alba L.) is a desirable tree spePhysocarpus amurensis Maxim andcies with high ecological and economic benefits. To reveal the response and adaptive mechanisms of the photosynthetic functions of mulberry leaves to flooding stress, chlorophyll synthesis, photosynthetic electron transfer and the Calvin cycle were investigated by physiological studies combined with an analysis of the transcriptome. Flooding stress inhibited the synthesis of chlorophyll (Chl) and decreased its content in mulberry leaves. The sensitivity of Chl a to flooding stress was higher than that of Chl b owing to the changes of CHLG (LOC21385082) and CAO (LOC21408165) that encode genes during chlorophyll synthesis. The levels of expression of Chl b reductase NYC (LOC112094996) and NYC (LOC21385774), which are involved in Chl b degradation, were upregulated on the fifteenth day of flooding, which accelerated the transformation of Chl b to Chl a, and upregulated the expression of PPH (LOC21385040) and PAO (LOC21395013). This accelerated the degradation of chlorophyll. Flooding stress significantly inhibited the photosynthetic function of mulberry leaves. A Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of differentially expressed genes under different days of flooding stress indicated significant enrichment in Photosynthesis-antenna proteins (map00196), Photosynthesis (map00195) and Carbon fixation in photosynthetic organisms (map00710). On the fifth day of flooding, 7 and 5 genes that encode antenna proteins were identified on LHCII and LHCI, respectively. They were significantly downregulated, and the degree of downregulation increased as the trees were flooded longer. Therefore, the power of the leaves to capture solar energy and transfer this energy to the reaction center was reduced. The chlorophyll fluorescence parameters and related changes in the expression of genes in the transcriptome indicated that the PSII and PSI of mulberry leaves were damaged, and their activities decreased under flooding stress. On the fifth day of flooding, electron transfer on the PSII acceptor side of mulberry leaves was blocked, and the oxygen-evolving complex (OEC) on the donor side was damaged. On the tenth day of flooding, the thylakoid membranes of mulberry leaves were damaged. Five of the six coding genes that mapped to the OEC were significantly downregulated. Simultaneously, other coding genes located at the PSII reaction center and those located at the PSI reaction center, including Cytb6/f, PC, Fd, FNR and ATP, were also significantly downregulated. In addition, the gas exchange parameters (P_n, G_s, T_r, and C_i) of the leaves decreased after 10 days of flooding stress primarily owing to the stomatal factor. However, on the fifteenth day of flooding, the value for the intracellular concentration of CO₂ was significantly higher than that on the tenth day of flooding. In addition, the differentially expressed genes identified in the Calvin cycle were significantly downregulated, suggesting that in addition to stomatal factors, non-stomatal factors were also important factors that mediated the decrease in the photosynthetic capacity of mulberry leaves. In conclusion, the inhibition of growth of mulberry plants caused by flooding stress was primarily related to the inhibition of chlorophyll synthesis, antenna proteins, photosynthetic electron transfer and the Calvin cycle. The results of this study provide a theoretical basis for the response and mechanism of adaptation of the photosynthetic function of mulberry to flooding stress.

Pinpointing Genomic Regions and Candidate Genes Associated with Seed Oil and Protein Content in Soybean through an Integrative Transcriptomic and QTL Meta-Analysis

Soybean with enriched nutrients has emerged as a prominent source of edible oil and protein. In the present study, a meta-analysis was performed by integrating quantitative trait loci (QTLs) information, region-specific association and transcriptomic analysis. Analysis of about a thousand QTLs previously identified in soybean helped to pinpoint 14 meta-QTLs for oil and 16 meta-QTLs for protein content. Similarly, region-specific association analysis using whole genome re-sequenced data was performed for the most promising meta-QTL on chromosomes 6 and 20. Only 94 out of 468 genes related to fatty acid and protein metabolic pathways identified within the meta-QTL region were found to be expressed in seeds. Allele mining and haplotyping of these selected genes were performed using whole genome resequencing data. Interestingly, a significant haplotypic association of some genes with oil and protein content was observed, for instance, in the case of FAD2-1B gene, an average seed oil content of 20.22% for haplotype 1 compared to 15.52% for haplotype 5 was observed. In addition, the mutation S86F in the FAD2-1B gene produces a destabilizing effect of (ΔΔG Stability) -0.31 kcal/mol. Transcriptomic analysis revealed the tissue-specific expression of candidate genes. Based on their higher expression in seed developmental stages, genes such as sugar transporter, fatty acid desaturase (FAD), lipid transporter, major facilitator protein and amino acid transporter can be targeted for functional validation. The approach and information generated in the present study will be helpful in the map-based cloning of regulatory genes, as well as for marker-assisted breeding in soybean.

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.

Flooding and herbivory: the effect of concurrent stress factors on plant volatile emissions and gene expression in two heirloom tomato varieties

BACKGROUND

In nature and in cultivated fields, plants encounter multiple stress factors. Nonetheless, our understanding of how plants actively respond to combinatorial stress remains limited. Among the least studied stress combination is that of flooding and herbivory, despite the growing importance of these stressors in the context of climate change. We investigated plant chemistry and gene expression changes in two heirloom tomato varieties: Cherokee Purple (CP) and Striped German (SG) in response to flooding, herbivory by Spodoptera exigua, and their combination.

RESULTS

Volatile organic compounds (VOCs) identified in tomato plants subjected to flooding and/or herbivory included several mono- and sesquiterpenes. Flooding was the main factor altering VOCs emission rates, and impacting plant biomass accumulation, while different varieties had quantitative differences in their VOC emissions. At the gene expression levels, there were 335 differentially expressed genes between the two tomato plant varieties, these included genes encoding for phenylalanine ammonia-lyase (PAL), cinnamoyl-CoA-reductase-like, and phytoene synthase (Psy1). Flooding and variety effects together influenced abscisic acid (ABA) signaling genes with the SG variety showing higher levels of ABA production and ABA-dependent signaling upon flooding. Flooding downregulated genes associated with cytokinin catabolism and general defense response and upregulated genes associated with ethylene biosynthesis, anthocyanin biosynthesis, and gibberellin biosynthesis. Combining flooding and herbivory induced the upregulation of genes including chalcone synthase (CHS), PAL, and genes encoding BAHD acyltransferase and UDP-glucose iridoid glucosyltransferase-like genes in one of the tomato varieties (CP) and a disproportionate number of heat-shock proteins in SG. Only the SG variety had measurable changes in gene expression due to herbivory alone, upregulating zeatin, and O-glucosyltransferase and thioredoxin among others.

CONCLUSION

Our results suggest that both heirloom tomato plant varieties differ in their production of secondary metabolites including phenylpropanoids and terpenoids and their regulation and activation of ABA signaling upon stress associated with flooding. Herbivory and flooding together had interacting effects that were evident at the level of plant chemistry (VOCs production), gene expression and biomass markers. Results from our study highlight the complex nature of plant responses to combinatorial stresses and point at specific genes and pathways that are affected by flooding and herbivory combined.

PacBio full-length sequencing integrated with RNA-seq reveals the molecular mechanism of waterlogging and its recovery in Paeonia ostii

Paeonia ostii, a widely cultivated tree peony species in China, is a resourceful plant with medicinal, ornamental and oil value. However, fleshy roots lead to a low tolerance to waterlogging in P. ostii. In this study, P. ostii roots were sequenced using a hybrid approach combining single-molecule real-time and next-generation sequencing platforms to understand the molecular mechanism underlying the response to this sequentially waterlogging stress, the normal growth, waterlogging treatment (WT), and waterlogging recovery treatment (WRT). Our results indicated that the strategy of P. ostii, in response to WT, was a hypoxic resting syndrome, wherein the glycolysis and fermentation processes were accelerated to maintain energy levels and the tricarboxylic acid cycle was inhibited. P. ostii enhanced waterlogging tolerance by reducing the uptake of nitrate and water from the soil. Moreover, transcription factors, such as AP2/EREBP, WRKY, MYB, and NAC, played essential roles in response to WT and WRT. They were all induced in response to the WT condition, while the decreasing expression levels were observed under the WRT condition. Our results contribute to understanding the defense mechanisms against waterlogging stress in P. ostii.

Comparative physiological and coexpression network analyses reveal the potential drought tolerance mechanism of peanut

BACKGROUND

Drought stress has negative effects on plant growth and productivity. In this study, a comprehensive analysis of physiological responses and gene expression was performed. The responses and expressions were compared between drought-tolerant (DT) and drought-sensitive (DS) peanut varieties to investigate the regulatory mechanisms and hub genes involved in the impact of drought stress on culture.

RESULTS

The drought-tolerant variety had robust antioxidative capacities with higher total antioxidant capacity and flavonoid contents, and it enhanced osmotic adjustment substance accumulation to adapt to drought conditions. KEGG analysis of differentially expressed genes demonstrated that photosynthesis was strongly affected by drought stress, especially in the drought-sensitive variety, which was consistent with the more severe suppression of photosynthesis. The hub genes in the key modules related to the drought response, including genes encoding protein kinase, E3 ubiquitin-protein ligase, potassium transporter, pentatricopeptide repeat-containing protein, and aspartic proteinase, were identified through a comprehensive combined analysis of genes and physiological traits using weighted gene co-expression network analysis. There were notably differentially expressed genes between the two varieties, suggesting the positive roles of these genes in peanut drought tolerance.

CONCLUSION

A comprehensive analysis of physiological traits and relevant genes was conducted on peanuts with different drought tolerances. The findings revealed diverse drought-response mechanisms and identified candidate genes for further research.

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.

Drought resistance of tobacco overexpressing the AfNAC1 gene of Amorpha fruticosa Linn

Plants are often adversely affected by abiotic stresses such as drought, low temperature, and salinity during growth, and plant NAC-like transcription factors are involved in regulating growth and developmental processes in response to stresses such as drought and salinity. In this study, to investigate the function of AfNAC1, a co-expression network of AfNAC1 genes was constructed using gene expression data from the Chinese legume deciduous shrub, Amorpha fruticosa Linn. A 576 bp NAC transcription factor (AfNAC1 gene, MN180266) encoding 191 amino acids was isolated from Amorpha fruticosa seedlings by RT-PCR. qRT-PCR showed that the AfNAC1 gene was expressed in the roots, stems, leaves, and flowers of Amorpha fruticosa. However, drought stress significantly increased root expression, and the AfNAC1 protein was localized in the nucleus by green fluorescence detection. This study analyzed the drought resistance of overexpressing tobacco in depth. Under natural drought stress, the chlorophyll and antioxidant enzyme activities of overexpressing plants were significantly higher than those of wild-type (WT) plants, but the MDA content was lower than that of WT; after rehydration the Fv/Fm values of AfNAC1-overexpressing tobacco recovered faster than those of wild-type tobacco and rapidly reached the control levels; AfNAC1 may be involved in the regulation of the photosystem and indirectly in the regulation of the plant in response to drought stress.

The Choice of Soybean Cultivar Alters the Underyielding of Protein and Oil under Drought Conditions in Central Poland

The popularization of soybean cultivation in Central Poland is progressing due to the European Soy Declaration signed by 13 member states in Brussels on 17 July 2017. Hence, this research was initiated under the European Innovation Partnership on phenotyping soybean cultivars in two regions, i.e., Kuyavian-Pomeranian and Greater Poland for integrated systems. The aim of this study was to determine soybean potential in the agrotechnical and agroclimatic conditions, with the selection of the most suitable cultivars for tillage and no-tillage cultivation in the region of Central Poland. Strict field research was carried out in six locations from 2018–2020 on 20 cultivars selected in terms of their earliness of maturation to the climatic conditions. On the basis of meteorological data, it was found that half of the plantations suffered from drought stress, as evidenced by the hydrothermal coefficients (K) for the growing season of soybean K < 1.2. The significant multivariate relationship (R = 0.67; p < 0.001) between the hydrothermal coefficient K and the vegetation period days (VPD) as the predictors of the soybean yield was determined by the regression equation Y(yield) = 21.8 + 12.0 X(K) + 0.20 X(VPD). A significant multivariate relationship (R = 0.43; p < 0.01) was also found between the seed yield, hydrothermal coefficient K and the protein content, quantified with the regression equation Y(protein) = 32.6 + 0.25 X(Yield) + 0.28 X(K). In the no-tillage system of cultivation, the seed yield of soybean constantly increased with increasing K (+32.3% between dry/relatively dry and optimal periods, +22.4% between optimal and humid periods), while in the tillage system, the yield increased by 22.1% only when K rose from optimal to humid. In the seasons with optimum and humid conditions, no-tilled soybean produced more oil, in comparison to the season of dry conditions, and the opposite trend was found with the greater oil content in the tillage system. Meanwhile, a higher protein content was observed in the tillage system under humid conditions. All soybean cultivars were grouped according to the cluster analysis (k-means) with ANOVA in terms of vegetation period in days, seed yield, oil, and protein content for tillage and no-tillage cultivation. The study confirmed that soybean might be considered in the future as a profitable crop in Central Poland and create perspective for a low-input source of protein and oil.

The AP2/ERF GmERF113 Positively Regulates the Drought Response by Activating GmPR10-1 in Soybean

Ethylene response factors (ERFs) are involved in biotic and abiotic stress; however, the drought resistance mechanisms of many ERFs in soybeans have not been resolved. Previously, we proved that GmERF113 enhances resistance to the pathogen Phytophthora sojae in soybean. Here, we determined that GmERF113 is induced by 20% PEG-6000. Compared to the wild-type plants, soybean plants overexpressing GmERF113 (GmERF113-OE) displayed increased drought tolerance which was characterized by milder leaf wilting, less water loss from detached leaves, smaller stomatal aperture, lower Malondialdehyde (MDA) content, increased proline accumulation, and higher Superoxide dismutase (SOD) and Peroxidase (POD) activities under drought stress, whereas plants with GmERF113 silenced through RNA interference were the opposite. Chromatin immunoprecipitation and dual effector-reporter assays showed that GmERF113 binds to the GCC-box in the GmPR10-1 promoter, activating GmPR10-1 expression directly. Overexpressing GmPR10-1 improved drought resistance in the composite soybean plants with transgenic hairy roots. RNA-seq analysis revealed that GmERF113 downregulates abscisic acid 8′-hydroxylase 3 (GmABA8’-OH 3) and upregulates various drought-related genes. Overexpressing GmERF113 and GmPR10-1 increased the abscisic acid (ABA) content and reduced the expression of GmABA8’-OH3 in transgenic soybean plants and hairy roots, respectively. These results reveal that the GmERF113-GmPR10-1 pathway improves drought resistance and affects the ABA content in soybean, providing a theoretical basis for the molecular breeding of drought-tolerant soybean.

The Hidden Potential of High-Throughput RNA-Seq Re-Analysis, a Case Study for DHDPS, Key Enzyme of the Aspartate-Derived Lysine Biosynthesis Pathway and Its Role in Abiotic and Biotic Stress Responses in Soybean

DHDPS is a key enzyme in the aspartate-derived lysine biosynthesis pathway and an evident object of study for biofortification strategies in plants. DHDPS isoforms with novel regulatory properties in Medicago truncatula were demonstrated earlier and hypothesized to be involved in abiotic and biotic stress responses. Here, we present a phylogenetic analysis of the DHPDS gene family in land plants which establishes the existence of a legume-specific class of DHDPS, termed DHDPS B-type, distinguishable from the DHDPS A-type commonly present in all land plants. The G. max genome comprises two A-type DHDPS genes (Gm.DHDPS-A1; Glyma.09G268200, Gm.DHDPS-A2; Glyma.18G221700) and one B-type (Gm.DHDPS-B; Glyma.03G022300). To further investigate the expression pattern of the G. max DHDPS isozymes in different plant tissues and under various stress conditions, 461 RNA-seq experiments were exploited and re-analyzed covering two expression atlases, 13 abiotic and 5 biotic stress studies. Gm.DHDPS-B is seen almost exclusively expressed in roots and nodules in addition to old cotyledons or senescent leaves while both DHDPS A-types are expressed constitutively in all tissues analyzed with the highest expression in mature seeds. Furthermore, Gm.DHDPS-B expression is significantly upregulated in some but not all stress responses including salt stress, flooding, ethylene or infection with Phytophthora sojae and coincides with downregulation of DHDPS A-types. In conclusion, we demonstrate the potential of an in-depth RNA-seq re-analysis for the guidance of future experiments and to expand on current knowledge.

Integrated Transcriptomic and Proteomic Analyses Uncover the Regulatory Mechanisms of Myricaria laxiflora Under Flooding Stress

Flooding is one of the major environmental stresses that severely influence plant survival and development. However, the regulatory mechanisms underlying flooding stress remain largely unknown in Myricaria laxiflora, an endangered plant mainly distributed in the flood zone of the Yangtze River, China. In this work, transcriptome and proteome were performed in parallel in roots of M. laxiflora during nine time-points under the flooding and post-flooding recovery treatments. Overall, highly dynamic and stage-specific expression profiles of genes/proteins were observed during flooding and post-flooding recovery treatment. Genes related to auxin, cell wall, calcium signaling, and MAP kinase signaling were greatly down-regulated exclusively at the transcriptomic level during the early stages of flooding. Glycolysis and major CHO metabolism genes, which were regulated at the transcriptomic and/or proteomic levels with low expression correlations, mainly functioned during the late stages of flooding. Genes involved in reactive oxygen species (ROS) scavenging, mitochondrial metabolism, and development were also regulated exclusively at the transcriptomic level, but their expression levels were highly up-regulated upon post-flooding recovery. Moreover, the comprehensive expression profiles of genes/proteins related to redox, hormones, and transcriptional factors were also investigated. Finally, the regulatory networks of M. laxiflora in response to flooding and post-flooding recovery were discussed. The findings deepen our understanding of the molecular mechanisms of flooding stress and shed light on the genes and pathways for the preservation of M. laxiflora and other endangered plants in the flood zone.

Carbon fluxes and environmental interactions during legume development, with a specific focus on Pisum sativum

Grain legumes are major food crops cultivated worldwide for their seeds with high nutritional content. To answer the growing concern about food safety and protein autonomy, legume cultivation must increase in the coming years. In parallel, current agricultural practices are facing environmental challenges, including global temperature increase and more frequent and severe episodes of drought stress. Crop yield directly relies on carbon allocation and is particularly affected by these global changes. We review the current knowledge on source-sink relationships and carbon resource allocation at all developmental stages, from germination to vegetative growth and seed production in grain legumes, focusing on pea (Pisum sativum). We also discuss how these source-sink relationships and carbon fluxes are influenced by biotic and abiotic factors. Major agronomic traits, including seed yield and quality, are particularly impacted by drought, temperatures, salinity, waterlogging, or pathogens and can be improved through the promotion of beneficial soil microorganisms or through optimized plant carbon resource allocation. Altogether, our review highlights the need for a better understanding of the cellular and molecular mechanisms regulating carbon fluxes from source leaves to sink organs, roots, and seeds. These advancements will further improve our understanding of yield stability and stress tolerance and contribute to the selection of climate-resilient crops.

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.

Full-Length Transcriptome and RNA-Seq Analyses Reveal the Mechanisms Underlying Waterlogging Tolerance in Kiwifruit (Actinidia valvata)

Actinidia valvata possesses waterlogging tolerance; however, the mechanisms underlying this trait are poorly characterized. Here, we performed a transcriptome analysis by combining single-molecule real-time (SMRT) sequencing and Illumina RNA sequencing and investigated the physiological responses of the roots of KR5 (A. valvata, a tolerant genotype) after 0, 12, 24 and 72 h of waterlogging stress. KR5 roots responded to waterlogging stress mainly via carbohydrate and free amino acids metabolism and reactive oxygen species (ROS) scavenging pathways. Trehalose-6-phosphate synthase (TPS) activity, alcohol dehydrogenase (ADH) activity and the total free amino acid content increased significantly under waterlogging stress. The nicotinamide adenine dinucleotide-dependent glutamate synthase/alanine aminotransferase (NADH-GOGAT/AlaAT) cycle was correlated with alanine accumulation. Levels of genes encoding peroxidase (POD) and catalase (CAT) decreased and enzyme activity increased under waterlogging stress. Members of the LATERAL ORGAN BOUNDARIES (LOB), AP2/ERF-ERF, Trihelix and C3H transcription factor families were identified as potential regulators of the transcriptional response. Several hub genes were identified as key factors in the response to waterlogging stress by a weighted gene co-expression network analysis (WGCNA). Our results provide insights into the factors contributing to waterlogging tolerance in kiwifruit, providing a basis for further studies of interspecific differences in an important plant trait and for molecular breeding.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Variety-specific transcript accumulation during reproductive stage in drought-stressed rice

The divergence of natural stress tolerance mechanisms between species is an intriguing phenomenon. To study it in rice, a comparative transcriptome analysis was carried out in 'heading' stage tissue (flag leaf, panicles and roots) of Nagina 22 (N22; drought-tolerant) and IR64 (drought-sensitive) plants subjected to field drought. Interestingly, N22 showed almost double the number of differentially expressed genes (DEGs) than IR64. Many DEGs colocalized within drought-related QTLs responsible for grain yield and drought tolerance and also associated with drought tolerance and critical drought-related plant traits such as leaf rolling, trehalose content, sucrose and cellulose content. Besides, co-expression analysis of the DEGs revealed several 'hub' genes known to actively regulate drought stress response. Strikingly, 1366 DEGs, including 21 'hub' genes, showed a distinct opposite regulation in the two rice varieties under similar drought conditions. Annotation of these variety-specific DEGs (VS-DEGs) revealed that they are distributed in various biological pathways. Furthermore, 103 VS-DEGs were found to physically interact with over 1300 genes, including 32 that physically interact with other VS-DEGs as well. The promoter region of these genes has sequence variations among the two rice varieties, which might be in part responsible for their unique expression pattern.

Characterization of Root System Architecture Traits in Diverse Soybean Genotypes Using a Semi-Hydroponic System

Phenotypic variation and correlations among root traits form the basis for selecting and breeding soybean varieties with efficient access to water and nutrients and better adaptation to abiotic stresses. Therefore, it is important to develop a simple and consistent system to study root traits in soybean. In this study, we adopted the semi-hydroponic system to investigate the variability in root morphological traits of 171 soybean genotypes popularized in the Yangtze and Huaihe River regions, eastern China. Highly diverse phenotypes were observed: shoot height (18.7-86.7 cm per plant with a median of 52.3 cm); total root length (208-1663 cm per plant with a median of 885 cm); and root mass (dry weight) (19.4-251 mg per plant with a median of 124 mg). Both total root length and root mass exhibited significant positive correlation with shoot mass (p ≤ 0.05), indicating their relationship with plant growth and adaptation strategies. The nine selected traits contributed to one of the two principal components (eigenvalues > 1), accounting for 78.9% of the total genotypic variation. Agglomerative hierarchical clustering analysis separated the 171 genotypes into five major groups based on these root traits. Three selected genotypes with contrasting root systems were validated in soil-filled rhizoboxes (1.5 m deep) until maturity. Consistent ranking of the genotypes in some important root traits at various growth stages between the two experiments indicates the reliability of the semi-hydroponic system in phenotyping root trait variability at the early growth stage in soybean germplasms.

Physiological, Biochemical and Transcriptomic Analysis of the Aerial Parts (Leaf-Blade and Petiole) of Asarum sieboldii Responding to Drought Stress

Asarum sieboldii Miq. is a leading economic crop and a traditional medicinal herb in China. Leaf-blade and petiole are the only aerial tissues of A. sieboldii during the vegetative growth, playing a vital role in the accumulation and transportation of biomass energy. They also act as critical indicators of drought in agricultural management, especially for crops having underground stems. During drought, variations in the morphology and gene expression of the leaves and petioles are used to control agricultural irrigation and production. Besides, such stress can also alter the differential gene expression in these tissues. However, little is known about the drought-tolerant character of the aerial parts of A. sieboldii. In this study, we examined the physiological, biochemical and transcriptomic responses to the drought stress in the leaf blades and petioles of A. sieboldii. The molecular mechanism, involving in drought stress response, was elucidated by constructing the cDNA libraries and performing transcriptomic sequencing. Under drought stress, a total of 2912 and 2887 unigenes were differentially expressed in the leaf blade and petiole, respectively. The detection of many transcription factors and functional genes demonstrated that multiple regulatory pathways were involved in drought tolerance. In response to drought, the leaf blade and petiole displayed a general physiological character, a higher SOD and POD activity, a higher MDA content and lower chlorophyll content. Three unigenes encoding POD were up-regulated, which can improve POD activity. Essential oil in petiole was extracted. The relative contents of methyleugenol and safrole in essential oil were increased from 0.01% to 0.05%, and 3.89% to 16.97%, respectively, while myristicin slightly reduced from 24.87% to 21.52%. Additionally, an IGS unigene, involved in eugenol biobiosynthesis, was found up-regulated under drought stress, which was predicated to be responsible for the accumulation of methyleugenol and safrole. Simple sequence repeats (SSRs) were characterized in of A. sieboldii, and a total of 5466 SSRs were identified. Among them, mono-nucleotides were the most abundant repeat units, accounting for 44.09% followed by tri-, tetra-, penta and hexa-nucleotide repeats. Overall, the present work provides a valuable resource for the population genetics studies of A. sieboldii. Besides, it provides much genomic information for the functional dissection of the drought-resistance in A. sieboldii, which will be useful to understand the bio-regulatory mechanisms linked with drought-tolerance to enhance its yield.

Identifying Traits Associated With Terminal Drought Tolerance in Sesame (Sesamum indicum L.) Genotypes

Sesame is predominantly cultivated in rainfed and low fertile lands and is frequently exposed to terminal drought. Sesamum species inhabiting dryland ecosystems adaptively diverge from those inhabiting rainfed habitats, and drought-specific traits have a genetic basis. In sesame, traits associated with drought conditions have not been explored to date, yet studies of these traits are needed given that drought is predicted to become more frequent and severe in many parts of the world because of climate change. Here, 76 accessions from the available Indian core set were used to quantify variation in several traits under irrigated (WW) and terminal drought stress (WS) conditions as well as their association with seed yield over two consecutive years. The range of trait variation among the studied genotypes under WW and WS was significant. Furthermore, the traits associated with seed yield under WW and WS differed. The per se performance of the accessions indicated that the expression of most traits was reduced under WS. The correlation analysis revealed that the number of branches, leaf area (LA), leaves dry weight (LDW), number of capsules plant^-1, and harvest index (HI) were positively correlated with seed yield under WW and WS, and total dry matter (TDM), plant stem weight, and canopy temperature (CT) were negatively correlated with seed yield under WW and WS, indicating that smaller and cooler canopy genotypes had higher yields. The genotypes IC-131936, IC-204045, IC-204861, IC-205363, IC-205311, and IC-73576 with the highest seed yields were characterized by low canopy temperature, high relative water content, and high harvest index under WS. Phenotypic and molecular diversity analysis was conducted on genotypes along with checks. Phenotypic diversity was assessed using multivariate analysis, whereas molecular diversity was estimated using simple sequence repeat (SSR) loci to facilitate the use of sesame in breeding and genetic mapping. SSRs showed low allelic variation, as indicated by a low average number of alleles (2.31) per locus, gene diversity (0.25), and polymorphism information content (0.22). Cluster analysis (CA) [neighbor-joining (NJ) tree] revealed three major genotypic groups and structure analysis showed 4 populations. The diverse genotypes identified with promising morpho-physiological traits can be used in breeding programs to develop new varieties.

Aerenchyma formation in the root of leaf-vegetable sweet potato: Programmed cell death initiated by ethylene-mediated H2 O2 accumulation

Sweet potato, commonly planted in Southeast Asia and South America with abundant rainfall, often suffers from waterlogging. The aerenchyma formation in roots is an effective way for plants to facilitate gas exchange. In the present study, tolerant and sensitive varieties, respectively, designated NC1 and C211, were evaluated under water oxygen content at 2.0 mg·L^-1 (hypoxia treatment) and 8.0 mg·L^-1 (control). The results showed that NC1 variety has a relatively higher root growth rate under low oxygen condition. In NC1 plants, aerenchyma was observed in the mid-section of the main adventitious root and spread to the proximal and distal ends, forming a complete channel in the cortex. However, in C211 plants, the aerenchyma occurred relatively later and could not turn into a whole channel. Ethylene synthesis-related (ACS1, ACS4, ACS5, etc.) and signal transduction-related (ETR1, ERS1, EIN2, etc.) genes were upregulated in the NC1 plants and led to changes in the reactive oxygen species-related genes (RBOHA, SOD, CAT, etc.) and enzyme activities. It was found that programmed cell death was induced by H₂ O₂ accumulation. A regulatory model of lysigenous aerenchyma formation in the root of sweet potato was constructed. Our study enriches the understanding of the mechanisms of the aerenchyma formation in plants.

Physiological and Expressional Regulation on Photosynthesis, Starch and Sucrose Metabolism Response to Waterlogging Stress in Peanut

Waterlogging has negative effects on crop yield. Physiological and transcriptome data of two peanut cultivars [Zhongkaihua 1 (ZKH 1) and Huayu 39 (HY 39)] were studied under normal water supply and waterlogging stress for 5 or 10 days at the flowering stage. The results showed that the main stem height, the number of lateral branches, lateral branch length, and the stem diameter increased under waterlogging stress, followed by an increase in dry matter accumulation, which was correlated with the increase in the soil and plant analysis development (SPAD) and net photosynthetic rate (Pn) and the upregulation of genes related to porphyrin and chlorophyll metabolism and photosynthesis. However, the imbalance of the source-sink relationship under waterlogging was the main cause of yield loss, and waterlogging caused an increase in the sucrose and soluble sugar contents and a decrease in the starch content; it also decreased the activities of sucrose synthetase (SS) and sucrose phosphate synthetase (SPS), which may be due to the changes in the expression of genes related to starch and sucrose metabolism. However, the imbalance of the source-sink relationship led to the accumulation of photosynthate in the stems and leaves, which resulted in the decrease of the ratio of pod dry weight to total dry weight (PDW/TDW) and yield. Compared with ZKH 1, the PDW of HY 39 decreased more probably because more photosynthate accumulated in the stem and leaves of HY 39 and could not be effectively transported to the pod.

Based on the whole genome clarified the evolution and expression process of fatty acid desaturase genes in three soybeans

Soybean is an important oil crop cultivated worldwide. With the increasing global population crossed with growing challenging cultivation conditions, improving soybean breeding by selecting important traits is urgent needed. Genes coding for plant fatty acid desaturases (FADs) genes are major candidates for that, because they are involving in controlling fatty acid composition and holding membrane fluidity under abiotic stress. Here, 75 FADs were found in three soybean genomes, which were further classified into four sub-groups. Phylogenetic tree, gene structure, motif and promoter analysis showed that the FAD gene family was conserved in the three soybeans. In addition, the numbers of omega desaturase from Chinese cultivated varieties were significantly higher than those in Chinese wild soybean and ancient polyploid soybean, respectively. However, it was the opposite for the sphingolipid subfamily. These results indicated that each subfamily was subjected to different selection pressures during cultivation and domestication. As the extra genes of the subfamily were very close to other family members' positions on chromosomes, they should be produced by duplication. The cis-element analysis of FAD promoter sequences revealed that upstream sequences of FAD contained abundant light, hormone and abiotic stress responsive cis-elements, suggesting that the quality of soybean could be improved by regulating these stresses. Expression analysis of Chinese wild soybean under salt stress showed that GsDES1.1, GsDES1.2, GsFAD2.1 and GsSLD1 in leaves and GsSLD2, GsSLD5 and GsSLD6 in roots were not closely related to salt stress response. Therefore, we explored the significant role of conserved, duplicated and neofunctionalized FAD in the domestication of soybean, which contributes to the importance of soybean as a global oil crop.

Overlapping and stress-specific transcriptomic and hormonal responses to flooding and drought in soybean

Flooding and drought are serious constraints that reduce crop productivity worldwide. Previous studies identified genes conferring tolerance to both water extremes in various plants. However, overlapping responses to flooding and drought at the genome-scale remain obscure. Here, we defined overlapping and stress-specific transcriptomic and hormonal responses to submergence, drought and recovery from these stresses in soybean (Glycine max). We performed comparative RNA-sequencing and hormone profiling, identifying genes, hormones and biological processes that are differentially regulated in an overlapping or stress-specific manner. Overlapping responses included positive regulation of trehalose and sucrose metabolism and negative regulation of cellulose, tubulin, photosystem II and I, and chlorophyll biosynthesis, facilitating the economization of energy reserves under both submergence and drought. Additional energy-consuming pathways were restricted in a stress-specific manner. Downregulation of distinct pathways for energy saving under each stress suggests energy-consuming processes that are relatively unnecessary for each stress adaptation are turned down. Our newly developed transcriptomic-response analysis revealed that abscisic acid and ethylene responses were activated in common under both stresses, whereas stimulated auxin response was submergence-specific. The energy-saving strategy is the key overlapping mechanism that underpins adaptation to both submergence and drought in soybean. Abscisic acid and ethylene are candidate hormones that coordinate transcriptomic energy-saving processes under both stresses. Auxin may be a signaling component that distinguishes submergence-specific regulation of the stress response.

Molecular Characterisation of Soybean Osmotins and Their Involvement in Drought Stress Response

Osmotins are multifunctional proteins belonging to the thaumatin-like family related to plant stress responses. To better understand the functions of soybean osmotins in drought stress response, the current study presents the characterisation of four previously described proteins and a novel putative soybean osmotin (GmOLPa-like). Gene and protein structure as well as gene expression analyses were conducted on different tissues and developmental stages of two soybean cultivars with varying dehydration sensitivities (BR16 and EMB48 are highly and slightly sensitive, respectively). The analysed osmotin sequences share the conserved amino acid signature and 3D structure of the thaumatin-like family. Some differences were observed in the conserved regions of protein sequences and in the electrostatic surface potential. P21-like present the most similar electrostatic potential to osmotins previously characterised as promoters of drought tolerance in Nicotiana tabacum and Solanum nigrum. Gene expression analysis indicated that soybean osmotins were differentially expressed in different organs (leaves and roots), developmental stages (R1 and V3), and cultivars in response to dehydration. In addition, under dehydration conditions, the highest level of gene expression was detected for GmOLPa-like and P21-like osmotins in the leaves and roots, respectively, of the less drought sensitive cultivar. Altogether, the results suggest an involvement of these genes in drought stress tolerance.

Comprehensive RNA-seq analysis revealed molecular pathways and genes associated with drought tolerance in wild soybean (Glycine soja Sieb. and Zucc.)

Drought stress at the germination stage is an important environmental stress limiting crop yield. Hence, our study investigated comparative root transcriptome profiles of four contrasting soybean genotypes viz., drought-tolerant (PI342618B/DTP and A214/DTL) and drought-sensitive (NN86-4/DSP and A195/DSL) under drought stress using RNA-Seq approach. A total of 4850 and 6272 differentially expressed genes (DEGs) were identified in tolerant (DTP and DTL) and sensitive (DSP and DSL) genotypes, respectively. Principle component analysis (PCA) and correlation analysis revealed higher correlation between DTP and DTL. Both gene ontology (GO) and MapMan analyses showed that the drought response was enriched in DEGs associated with water and auxin transport, cell wall/membrane, antioxidant activity, catalytic activity, secondary metabolism, signaling and transcription factor (TF) activities. Out of 981 DEGs screened from above terms, only 547 showed consistent opposite expression between contrasting genotypes. Twenty-eight DEGs of 547 were located on Chr.08 rich in QTLs and "Hotspot regions" associated with drought stress, and eight of them showed non-synonymous single nucleotide polymorphism. Hence, 10 genes (including above eight genes plus two hub genes) were predicated as possible candidates regulating drought tolerance, which needs further functional validation. Overall, the transcriptome profiling provided in-depth understanding about the genetic mechanism and candidate genes underlying drought tolerance in soybean.

Soybean transporter database: A comprehensive database for identification and exploration of natural variants in soybean transporter genes

Transporters, a class of membrane proteins that facilitate exchange of solutes including diverse molecules and ions across the cellular membrane, are vital component for the survival of all organisms. Understanding plant transporters is important to get insight of the basic cellular processes, physiology, and molecular mechanisms including nutrient uptake, signaling, response to external stress, and many more. In this regard, extensive analysis of transporters predicted in soybean and other plant species was performed. In addition, an integrated database for soybean transporter protein, SoyTD, was developed that will facilitate the identification, classification, and extensive characterization of transporter proteins by integrating expression, gene ontology, conserved domain and motifs, gene structure organization, and chromosomal distribution features. A comprehensive analysis was performed to identify highly confident transporters by integrating various prediction tools. Initially, 7541 transmembrane (TM) proteins were predicted in the soybean genome; out of these, 3306 non-redundant transporter genes carrying two or more transmembrane domains were selected for further analysis. The identified transporter genes were classified according to a standard transporter classification (TC) system. Comparative analysis of transporter genes among 47 plant genomes provided insights into expansion and duplication of transporter genes in land plants. The whole genome resequencing (WGRS) and tissue-specific transcriptome datasets of soybean were integrated to investigate the natural variants and expression profile associated with transporter(s) of interest. Overall, SoyTD provides a comprehensive interface to study genetic and molecular function of soybean transporters. SoyTD is publicly available at http://artemis.cyverse.org/soykb_dev/SoyTD/.

VvWRKY22 transcription factor interacts with VvSnRK1.1/VvSnRK1.2 and regulates sugar accumulation in grape

WRKY transcription factors are one of the largest transcriptional regulator families, involved in various signaling networks in plants. However, only limited functional exploration of the sugar signaling of Vitis vinifera WRKY22 transcription factor (VvWRKY22) has been conducted. In this study, the roles played by VvWRKY22 in sugar accumulation in grapes were investigated. VvWRKY22 was co-expressed with 16 sugar-related genes, and the expression of VvWRKY22 in grape suspension cells was inhibited by sucrose, and induced by fructose and abscisic acid (ABA). Results showed that over-expression of VvWRKY22 decreased the sucrose, glucose and fructose content, and regulated the expression levels of sugar and ABA-related genes. Moreover, it was found that VvWRKY22 interacted with VvSnRK1.1 or VvSnRK1.2 proteins (Sucrose non-fermenting-1-related protein kinase 1), which are important kinases related to sugar metabolism. These results, thus, provide new genetic evidences to support the view that VvWRKY22 functions in regulating sugar metabolism in grapes.

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.