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

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.

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

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

Whole-genome resequencing analysis of the medicinal plant Gardenia jasminoides

Background

Gardenia jasminoides is a species of Chinese medicinal plant, which has high medicinal and economic value and rich genetic diversity, but the study on its genetic diversity is far not enough.

Methods

In this study, one wild and one cultivated gardenia materials were resequenced using IlluminaHiSeq sequencing platform and the data were evaluated to understand the genomic characteristics of G. jasminoides.

Results

After data analysis, the results showed that clean data of 11.77G, Q30 reached 90.96%. The average comparison rate between the sample and reference genome was 96.08%, the average coverage depth was 15X, and the genome coverage was 85.93%. The SNPs of FD and YP1 were identified, and 3,087,176 and 3,241,416 SNPs were developed, respectively. In addition, SNP non-synonymous mutation, InDel mutation, SV mutation and CNV mutation were also detected between the sample and the reference genome, and KEGG, GO and COG database annotations were made for genes with DNA level variation. The structural gene variation in the biosynthetic pathway of crocin and gardenia, the main medicinal substance of G. jasminoides was further explored, which provided basic data for molecular breeding and genetic diversity of G. jasminoides in the future.

Unravelling the molecular mechanism underlying drought stress response in chickpea via integrated multi-omics analysis

Drought stress affects growth and productivity significantly in chickpea. An integrated multi-omics analysis can provide a better molecular-level understanding of drought stress tolerance. In the present study, comparative transcriptome, proteome and metabolome analyses of two chickpea genotypes with contrasting responses to drought stress, ICC 4958 (drought-tolerant, DT) and ICC 1882 (drought-sensitive, DS), was performed to gain insights into the molecular mechanisms underlying drought stress response/tolerance. Pathway enrichment analysis of differentially abundant transcripts and proteins suggested the involvement of glycolysis/gluconeogenesis, galactose metabolism, and starch and sucrose metabolism in the DT genotype. An integrated multi-omics analysis of transcriptome, proteome and metabolome data revealed co-expressed genes, proteins and metabolites involved in phosphatidylinositol signaling, glutathione metabolism and glycolysis/gluconeogenesis pathways, specifically in the DT genotype under drought. These stress-responsive pathways were coordinately regulated by the differentially abundant transcripts, proteins and metabolites to circumvent the drought stress response/tolerance in the DT genotype. The QTL-hotspot associated genes, proteins and transcription factors may further contribute to improved drought tolerance in the DT genotype. Altogether, the multi-omics approach provided an in-depth understanding of stress-responsive pathways and candidate genes involved in drought tolerance in chickpea.

Transcriptome Analysis Reveals Genes Associated with Flooding Tolerance in Mulberry Plants

Mulberry (Morus alba), a widely distributed economic plant, can withstand long-term flooding stress. However, the regulatory gene network underlying this tolerance is unknown. In the present study, mulberry plants were subjected to submergence stress. Subsequently, mulberry leaves were collected to perform quantitative reverse-transcription PCR (qRT-PCR) and transcriptome analysis. Genes encoding ascorbate peroxidase and glutathione S-transferase were significantly upregulated after submergence stress, indicating that they could protect the mulberry plant from flood damage by mediating ROS homeostasis. Genes that regulate starch and sucrose metabolism; genes encoding pyruvate kinase, alcohol dehydrogenase, and pyruvate decarboxylase (enzymes involved in glycolysis and ethanol fermentation); and genes encoding malate dehydrogenase and ATPase (enzymes involved in the TCA cycle) were also obviously upregulated. Hence, these genes likely played a key role in mitigating energy shortage during flooding stress. In addition, genes associated with ethylene, cytokinin, abscisic acid, and MAPK signaling; genes involved in phenylpropanoid biosynthesis; and transcription factor genes also showed upregulation under flooding stress in mulberry plants. These results provide further insights into the adaptation mechanisms and genetics of submergence tolerance in mulberry plants and could aid in the molecular breeding of these plants.

Integrative pathway and network analysis provide insights on flooding-tolerance genes in soybean

Soybean is highly sensitive to flooding and extreme rainfall. The phenotypic variation of flooding tolerance is a complex quantitative trait controlled by many genes and their interaction with environmental factors. We previously constructed a gene-pool relevant to soybean flooding-tolerant responses from integrated multiple omics and non-omics databases, and selected 144 prioritized flooding tolerance genes (FTgenes). In this study, we proposed a comprehensive framework at the systems level, using competitive (hypergeometric test) and self-contained (sum-statistic, sum-square-statistic) pathway-based approaches to identify biologically enriched pathways through evaluating the joint effects of the FTgenes within annotated pathways. These FTgenes were significantly enriched in 36 pathways in the Gene Ontology database. These pathways were related to plant hormones, defense-related, primary metabolic process, and system development pathways, which plays key roles in soybean flooding-induced responses. We further identified nine key FTgenes from important subnetworks extracted from several gene networks of enriched pathways. The nine key FTgenes were significantly expressed in soybean root under flooding stress in a qRT-PCR analysis. We demonstrated that this systems biology framework is promising to uncover important key genes underlying the molecular mechanisms of flooding-tolerant responses in soybean. This result supplied a good foundation for gene function analysis in further work.

Metabolomic and transcriptomic responses of Adiantum (Adiantum nelumboides) leaves under drought, half-waterlogging, and rewater conditions

Introduction: Adiantum nelumboides (Adiantum) is an endangered fern with a narrow distribution along the Yangtze River in China. Due to its cliff-dwelling habit, it experiences water stress conditions, which further endangers its survival. However, no information is available about its molecular responses to drought and half-waterlogging conditions. Methods: Here, we applied five and ten days of half-waterlogging stress, five days of drought stress, and rewatering after five days of drought stress, and studied the resulting metabolome profiles and transcriptome signatures of Adiantum leaves. Results and Discussion: The metabolome profiling detected 864 metabolites. The drought and half-waterlogging stress induced up-accumulation of primary and secondary metabolites including amino acids and derivatives, nucleotides and derivatives, flavonoids, alkaloids, and phenolic acid accumulation in Adiantum leaves. Whereas, rewatering the drought-stressed seedlings reversed most of these metabolic changes. Transcriptome sequencing confirmed the differential metabolite profiles, where the genes enriched in pathways associated with these metabolites showed similar expression patterns. Overall, the half-waterlogging stress for 10 days induced large-scale metabolic and transcriptomic changes compared to half-waterlogging stress for 05 days, drought stress for 05 days or rewatering for 05 days. Conclusion: This pioneering attempt provides a detailed understanding of molecular responses of Adiantum leaves to drought and half-waterlogging stresses and rewater conditions. This study also provides useful clues for the genetic improvement of Adiantum for drought/half-waterlogging stress 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.

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.

Morphological, Biochemical, and Proteomic Analyses to Understand the Promotive Effects of Plant-Derived Smoke Solution on Wheat Growth under Flooding Stress

Wheat is an important staple food crop for one-third of the global population; however, its growth is reduced by flooding. On the other hand, a plant-derived smoke solution enhances plant growth; however, its mechanism is not fully understood. To reveal the effects of the plant-derived smoke solution on wheat under flooding, morphological, biochemical, and proteomic analyses were conducted. The plant-derived smoke solution improved wheat-leaf growth, even under flooding. According to the functional categorization of proteomic results, oppositely changed proteins were correlated with photosynthesis, glycolysis, biotic stress, and amino-acid metabolism with or without the plant-derived smoke solution under flooding. Immunoblot analysis confirmed that RuBisCO activase and RuBisCO large/small subunits, which decreased under flooding, were recovered by the application of the plant-derived smoke solution. Furthermore, the contents of chlorophylls a and b significantly decreased by flooding stress; however, they were recovered by the application of the plant-derived smoke solution. In glycolysis, fructose-bisphosphate aldolase and glyceraldehyde-3-phosphate dehydrogenase decreased with the application of the plant-derived smoke solution under flooding as compared with flooding alone. Additionally, glutamine, glutamic acid, aspartic acid, and serine decreased under flooding; however, they were recovered by the plant-derived smoke solution. These results suggest that the application of the plant-derived smoke solution improves the recovery of wheat growth through the regulation of photosynthesis and glycolysis even under flooding conditions. Furthermore, the plant-derived smoke solution might promote wheat tolerance against flooding stress through the regulation of amino-acid metabolism.

Plant proteomic research for improvement of food crops under stresses: a review

Crop improvement approaches have been changed due to technological advancements in traditional plant-breeding methods. Abiotic and biotic stresses limit plant growth and development, which ultimately lead to reduced crop yield. Proteins encoded by genomes have a considerable role in the endurance and adaptation of plants to different environmental conditions. Biotechnological applications in plant breeding depend upon the information generated from proteomic studies. Proteomics has a specific advantage to contemplate post-translational modifications, which indicate the functional effects of protein modifications on crop production. Subcellular proteomics helps in exploring the precise cellular responses and investigating the networking among subcellular compartments during plant development and biotic/abiotic stress responses. Large-scale mass spectrometry-based plant proteomic studies with a more comprehensive overview are now possible due to dramatic improvements in mass spectrometry, sample preparation procedures, analytical software, and strengthened availability of genomes for numerous plant species. Development of stress-tolerant or resilient crops is essential to improve crop productivity and growth. Use of high throughput techniques with advanced instrumentation giving efficient results made this possible. In this review, the role of proteomic studies in identifying the stress-response processes in different crops is summarized. Advanced techniques and their possible utilization on plants are discussed in detail. Proteomic studies accelerate marker-assisted genetic augmentation studies on crops for developing high yielding stress-tolerant lines or varieties under stresses.

Proteomic and Biochemical Analyses of the Mechanism of Tolerance in Mutant Soybean Responding to Flooding Stress

To investigate the mechanism of flooding tolerance of soybean, flooding-tolerant mutants derived from gamma-ray irradiated soybean were crossed with parent cultivar Enrei for removal of other factors besides the genes related to flooding tolerance in primary generated mutant soybean. Although the growth of the wild type was significantly suppressed by flooding compared with the non-flooding condition, that of the mutant lines was better than that of the wild type even if it was treated with flooding. A two-day-old mutant line was subjected to flooding for 2 days and proteins were analyzed using a gel-free/label-free proteomic technique. Oppositely changed proteins in abundance between the wild type and mutant line under flooding stress were associated in endoplasmic reticulum according to gene-ontology categorization. Immunoblot analysis confirmed that calnexin accumulation increased in both the wild type and mutant line; however, calreticulin accumulated in only the mutant line under flooding stress. Furthermore, although glycoproteins in the wild type decreased by flooding compared with the non-flooding condition, those in the mutant line increased even if it was under flooding stress. Alcohol dehydrogenase accumulated in the wild type and mutant line; however, this enzyme activity significantly increased and mildly increased in the wild type and mutant line, respectively, under flooding stress compared with the non-flooding condition. Cell death increased and decreased in the wild type and mutant line, respectively, by flooding stress. These results suggest that the regulation of cell death through the fermentation system and glycoprotein folding might be an important factor for the acquisition of flooding tolerance in mutant soybean.

Recombinant Expression and Characterization of Lemon (Citrus limon) Peroxidase

BACKGROUND

Class III plant peroxidases play important role in a number of physiological processes in plants such as lignin biosynthesis, suberization, cell wall biosynthesis, reactive oxygen species metabolism and plant defense against pathogens. Peroxidases are also of significance in several industrial applications. In view of this, the production and identification of novel peroxidases having resistance towards temperature, pH, salts is desirable.

OBJECTIVE

The objective of the present work was to clone and characterize a novel plant peroxidase suitable for industrial application.

METHODS

A full length cDNA clone of lemon peroxidase was isolated using PCR and RACE approaches, characterized and heterologously expressed in Escherichia coli using standard protocols. The expressed peroxidase was purified using Ni-NTA agarose column and biochemically characterized using standard protocols. The peroxidase was also in-silico characterized at nucleotide as well as protein levels using standard protocols.

RESULTS

A full length cDNA clone of lemon peroxidase was isolated and expressed heterologously in E. coli. The expressed recombinant lemon peroxidase (LPRX) was activated by in-vitro refolding and purified. The purified LPRX exhibited pH and temperature optima of pH 7.0 and 50°C, respectively. The LPRX was found to be activated by metal ions (Na+, Ca2+, Mg2+ and Mn2+) at lower concentration. The expressional analysis of the transcripts suggested involvement of lemon peroxidase in plant defense. The lemon peroxidase was in silico modelled and docked with the substrates guaiacol, and pyrogallol and shown the favourability of pyrogallol over guaiacol, which is in agreement with the in-vitro findings. The protein function annotation analyses suggested the involvement of lemon peroxidase in the phenylpropanoid biosynthesis pathway and plant defense mechanisms.

CONCLUSION

Based on the biochemical characterization, the purified peroxidase was found to be resistant towards the salts and thus, might be a good candidate for industrial exploitation. The in-silico protein function annotation and transcript analyses highlighted the possible involvement of the lemon peroxidase in plant defense response.

Long-Term Waterlogging as Factor Contributing to Hypoxia Stress Tolerance Enhancement in Cucumber: Comparative Transcriptome Analysis of Waterlogging Sensitive and Tolerant Accessions

Waterlogging (WL), excess water in the soil, is a phenomenon often occurring during plant cultivation causing low oxygen levels (hypoxia) in the soil. The aim of this study was to identify candidate genes involved in long-term waterlogging tolerance in cucumber using RNA sequencing. Here, we also determined how waterlogging pre-treatment (priming) influenced long-term memory in WL tolerant (WL-T) and WL sensitive (WL-S) i.e., DH2 and DH4 accessions, respectively. This work uncovered various differentially expressed genes (DEGs) activated in the long-term recovery in both accessions. De novo assembly generated 36,712 transcripts with an average length of 2236 bp. The results revealed that long-term waterlogging had divergent impacts on gene expression in WL-T DH2 and WL-S DH4 cucumber accessions: after 7 days of waterlogging, more DEGs in comparison to control conditions were identified in WL-S DH4 (8927) than in WL-T DH2 (5957). Additionally, 11,619 and 5007 DEGs were identified after a second waterlogging treatment in the WL-S and WL-T accessions, respectively. We identified genes associated with WL in cucumber that were especially related to enhanced glycolysis, adventitious roots development, and amino acid metabolism. qRT-PCR assay for hypoxia marker genes i.e., alcohol dehydrogenase (adh), 1-aminocyclopropane-1-carboxylate oxidase (aco) and long chain acyl-CoA synthetase 6 (lacs6) confirmed differences in response to waterlogging stress between sensitive and tolerant cucumbers and effectiveness of priming to enhance stress tolerance.

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.

Genome-wide characterization and expression analysis of soybean trihelix gene family

Trihelix transcription factors play multiple roles in plant growth, development and various stress responses. In this study, we identified 71 trihelix family genes in the soybean genome. These trihelix genes were located at 19 out of 20 soybean chromosomes unevenly and were classified into six distinct subfamilies: GT-1, GT-2, GTγ, SIP1, SH4 and GTδ. The gene structure and conserved functional domain of these trihelix genes were similar in the same subfamily but diverged between different subfamilies. Thirteen segmental duplicated gene pairs were identified and all of them experienced a strong purifying selective pressure during evolution. Various stress-responsive cis-elements presented in the promoters of soybean trihelix genes, suggesting that the trihelix genes might respond to the environmental stresses in soybean. The expression analysis suggests that trihelix genes are involved in diverse functions during soybean development, flood or salinity tolerance, and plant immunity. Our results provide genomic information of the soybean trihelix genes and a basis for further characterizing their roles in response to environmental stresses.

A Comprehensive Transcriptomics Analysis Reveals Long Non-Coding RNA to be Involved in the Key Metabolic Pathway in Response to Waterlogging Stress in Maize

Waterlogging stress (WS) in a dynamic environment seriously limits plant growth, development, and yield. The regulatory mechanism underlying WS conditions at an early stage in maize seedlings is largely unknown. In the present study, the primary root tips of B73 seedlings were sampled before (0 h) and after (2 h, 4 h, 6 h, 8 h, 10 h, and 12 h) WS and then subjected to transcriptome sequencing, resulting in the identification of differentially expressed protein-coding genes (DEpcGs) and long non-coding RNAs (DElncRs) in response to WS. These DEpcGs were classified into nine clusters, which were significantly enriched in several metabolic pathways, such as glycolysis and methionine metabolism. Several transcription factor families, including AP2-EREBP, bZIP, NAC, bHLH, and MYB, were also significantly enriched. In total, 6099 lncRNAs were identified, of which 3190 were DElncRs. A co-expression analysis revealed lncRNAs to be involved in 11 transcription modules, 10 of which were significantly associated with WS. The DEpcGs in the four modules were enriched in the hypoxia response pathways, including phenylpropanoid biosynthesis, MAPK signaling, and carotenoid biosynthesis, in which 137 DElncRs were also co-expressed. Most of the co-expressed DElncRs were co-localized with previously identified quantitative trait loci associated with waterlogging tolerance. A quantitative reverse transcription-polymerase chain reaction analysis of DEpcG and DElncR expression among the 32 maize genotypes after 4 h of WS verified significant expression correlations between them as well as significant correlation with the phenotype of waterlogging tolerance. Moreover, the high proportion of hypoxia response elements in the promoter region increased the reliability of the DElncRs identified in this study. These results provide a comprehensive transcriptome in response to WS at an early stage of maize seedlings and expand our understanding of the regulatory network involved in hypoxia in plants.