Skip to main content

Physiological and comparative transcriptome analyses reveal the mechanisms underlying waterlogging tolerance in a rapeseed anthocyanin-more mutant

Abstract

Background

Rapeseed (Brassica napus) is the second largest oil crop worldwide. It is widely used in food, energy production and the chemical industry, as well as being an ornamental. Consequently, it has a large economic value and developmental potential. Waterlogging is an important abiotic stress that restricts plant growth and development. However, little is known about the molecular mechanisms underlying waterlogging tolerance in B. napus.

Results

In the present study, the physiological changes and transcriptomes of germination-stage rapeseed in response to waterlogging stress were investigated in the B. napus cultivar ‘Zhongshuang 11’ (ZS11) and its anthocyanin-more (am) mutant, which was identified in our previous study. The mutant showed stronger waterlogging tolerance compared with ZS11, and waterlogging stress significantly increased anthocyanin, soluble sugar and malondialdehyde contents and decreased chlorophyll contents in the mutant after 12 days of waterlogging. An RNA-seq analysis identified 1370 and 2336 differently expressed genes (DEGs) responding to waterlogging stress in ZS11 and am, respectively. An enrichment analysis revealed that the DEGs in ZS11 were predominately involved in carbohydrate metabolism, whereas those in the am mutant were particularly enriched in plant hormone signal transduction and response to endogenous stimulation. In total, 299 DEGs were identified as anthocyanin biosynthesis-related structural genes (24) and regulatory genes encoding transcription factors (275), which may explain the increased anthocyanin content in the am mutant. A total of 110 genes clustered in the plant hormone signal transduction pathway were also identified as DEGs, including 70 involved in auxin and ethylene signal transduction that were significantly changed in the mutant. Furthermore, the expression levels of 16 DEGs with putative roles in anthocyanin accumulation and biotic/abiotic stress responses were validated by quantitative real-time PCR as being consistent with the transcriptome profiles.

Conclusion

This study provides new insights into the molecular mechanisms of increased anthocyanin contents in rapeseed in response to waterlogging stress, which should be useful for reducing the damage caused by waterlogging stress and for further breeding new rapeseed varieties with high waterlogging tolerance.

Background

Rapeseed (Brassica napus) is a main oil crop worldwide and has the largest cultivated area in China. Rapeseed oil contains high oleic acid, low saturated fatty acid, moderate linoleic acid, linolenic acid and many vitamins and sterols, which are important nutrient resources in edible vegetable oil [1, 2]. In addition, the applications of rapeseed oil in chemical, pharmaceutical and bioenergy fields are increasing [3, 4].

During growth and development, rapeseed encounters various adverse stresses, such as low temperature, salt, drought, waterlogging, heavy metal ions and pathogens, which can hinder the growth and reduce the yield by affecting gene expression and cell physiological metabolism. Among these stresses, rapeseed is especially sensitive to waterlogging stress owing to the lack of aerenchymal tissue [5]. Similar to wheat and corn, waterlogging stress also results in decreases in the photosynthetic rate, chlorophyll content, root absorption function and antioxidant enzyme activity of rapeseed, hindering the development of pods and eventually leading to yield reductions [6,7,8,9,10,11]. The Yangtze River Basin is the main rapeseed planting area in China, and heavy rainfalls during the growing season in this area increase the soil moisture and humidity, which also accelerates the propagation and spread of Sclerotinia sclerotiorum and other pathogens [8, 12]. Although measures adopted in production, such as ditching and drainage, plowing the land and increasing organic fertilizer, alleviate the harm to rapeseed caused by waterlogging stress to a certain extent, they undoubtedly increase the labor costs [13]. Therefore, breeding new rapeseed varieties with high waterlogging tolerance is still the most economical and effective way to reduce the damage, and in-depth analyses of waterlogging-tolerance mechanism are of great significance to achieve this breeding goal.

During waterlogging stress, plants usually form aerenchyma and radial oxygen loss barriers in adventitious roots and lateral roots to increase the oxygen transport capacity [14,15,16] or inhibit root extension to reduce energy consumption under hypoxic conditions [17]. However, owing to the lack of aerenchyma and the sensitivity to water, rapeseed can form numerous adventitious roots to maintain a high energy reserve, which improves its adaptability under stress conditions [5, 18]. In addition to morphological changes, plants also induce some positive physiological and molecular responses to waterlogging stress [19]. For example, many metabolic pathways, including glycolytic, fermentation and respiratory metabolic pathways, are changed accordingly. Moreover, some endogenous plant hormones, such as abscisic acid (ABA), ethylene (ET), auxin and brassinolide (BR) are very sensitive to waterlogging stress and can effectively mitigate the adverse effects caused by waterlogging [8, 20,21,22,23,24]. Waterlogging-tolerant rapeseed varieties can also significantly increase the contents of osmoregulation organic matter, such as soluble sugar, soluble protein and proline [7, 25]. The application of nitrogen can also effectively improve the physiological indexes of rapeseed, such as root development, chlorophyll content and antioxidant enzyme activity levels, and reduce the degree of membrane lipid peroxidation to alleviate the damage caused by waterlogging stress [11].

Anthocyanins, as one of the most important water-soluble pigments in plants, have a wide range of biological functions, such as plant abiotic stress adaptation, defense against pathogen invasion, antioxidation and other health-related functions [26,27,28]. Plant anthocyanin production can be affected by environmental factors, plant hormones, miRNAs (such as miR156, miR165/166 and miR778) and nutritional conditions [29,30,31,32,33]. In the natural state, there are few free anthocyanins, they are predominantly in the form of glycosides combined by glycosylation [26, 34]. Anthocyanins are flavonoids, and their biosynthesis and accumulation are parts of the flavonoid biosynthetic pathway, which is mainly controlled by structural and regulatory genes [35, 36]. Structural genes encode various enzymes involved in anthocyanin synthesis, such as chalcone synthase (CHS), dihydroflavonol 4-reductase (DFR) and anthocyanidin synthase (ANS) [37], whereas regulatory genes are related to plant color and regulate the expression levels of structural genes through a conserved MBW (MYB-bHLH-WD40) complex, which is composed of R2R3-MYB transcription factors (TFs) PAP1/MYB75, PAP2/MYB90, MYB113 and MYB114, basic helix–loop–helix (bHLH) TFs TT8, GL3 and EGL3, and the WD40-repeat protein TTG1 in Arabidopsis [38,39,40]. Similar MBW complexes have been identified in different plant species, such as R2R3-MYB TFs MYB6, ROS1 and SlAN2-like [41,42,43,44], and bHLH TFs, such as PPLS1 and bHLH33, which positively regulate anthocyanin biosynthesis [45, 46], whereas MYBL2, MYBC2L2 and SlMYBATV serve as transcriptional repressors of anthocyanin biosynthesis [31, 43, 47]. Owing to the lack of anthocyanin-related mutants, structural and regulatory genes that are important in controlling anthocyanin production in B. napus are still not well understood and need further identification and study.

In our previous study, we identified an anthocyanin-more (am) mutant of a B. napus cultivar that produced more anthocyanins compared with the wild type (WT) during growth and development [48]. The am mutant was more tolerant to waterlogging stress during B. napus seed germination and seedling growth and developmental stages. To systematically analyze the genetic and biochemical bases for the different anthocyanin production capabilities between the am mutant and WT, and the molecular mechanisms of the different waterlogging responses with respect to the accumulation of anthocyanin, an RNA sequencing (RNA-seq) analysis of leaf tissues was carried out to analyze transcriptional responses to waterlogging stress over 24 h. This study will expand our understanding of the transcriptional regulation of anthocyanin production and the molecular mechanisms of anthocyanin-induced waterlogging resistance in a rapeseed am mutant, and this lays the foundation for further candidate gene identification and the molecular breeding of waterlogging-tolerant rapeseed.

Results

Differences in seed germination and seedling growth between the WT and am mutant under waterlogging-stress conditions

The am mutant is a natural mutant of B. napus ‘Zhongshuang 11’ (ZS11) identified in our previous study. Compared with WT, the mutant presents purple tissues during different growth and developmental stages of plant [48]. During the early germination (36 h) stages, the root length, stem length and fresh weight of the mutant were significantly higher than those of ZS11 after various time periods (12, 24 and 36 h) under waterlogging stress (Fig. 1a, b). Following an 8-day recovery period, the seedlings rate, growth quantity and growth vigor of the am mutant were all superior to those of the WT under the same waterlogging treatment conditions (Additional file 1: Fig. S1). Specifically, no obvious differences in the relative root lengths and fresh weights were found between the two materials after 12 h of waterlogging. After 24 h of treatment and 8 days of recovery, the relative seedling rates of the two materials decreased, but the decline range of am was less than that of ZS11. The relative vigor indexes of the two materials also showed significant differences, with ZS11 being only 0.18, whereas that of the am mutant was 0.65, indicating a greater tolerance to waterlogging compared with WT. After a 36-h treatment and recovery period, most ZS11 seedlings did not grow normally, and cotyledons could not completely break away from the husks and extend outwards, whereas 46.1% of the am mutant germinated and grew normally after waterlogging (Table 1, Additional file 1: Fig. S1). Moreover, the mutant was prone to produce anthocyanins during germination under waterlogging stress, and the stems, petioles and veins appeared purple, whereas the WT plant showed no obvious changes (Fig. 1).

Fig. 1
figure 1

Phenotypic differences between the wild type (WT) and am mutant at different periods after waterlogging. a The seeds germinated for 36 h were waterlogged for 0, 12, 24 and 36 h. b Effects of waterlogging on root lengths, stem lengths and fresh weights of the WT and am mutant. Error bars represent the SDs from three replicated experiments. Significant differences between the am mutant and WT at different time points are indicated (Student’s t-test) as follows: ***P < 0.001; **P < 0.01; *P < 0.05. c The young seedlings of the WT ZS11 and am mutant were waterlogged for 5, 7 and 12 d. d The morphological changes in the ZS11 and am mutant seedlings after waterlogging for 12 d. Bars = 1.0 cm

Table 1 Effects of waterlogging stress on the seed germination characteristics of ZS11 and the am mutant

At the seedling stage after waterlogging stress, the am mutant had purple stems and leaves prior to the WT, and to a certain extent, as the waterlogging time increased, the purple in these parts gradually deepened (Fig. 1b). When the waterlogging time reached 12 days, the mutant am seedlings showed purple-colored veins, petioles and leaf edge areas, and they had well developed root systems, compared with WT seedlings (Fig. 1c). Thus, the am mutants were more likely to produce anthocyanins under waterlogging-stress conditions during both the germination and seedling stages, and this may be related to waterlogging tolerance.

Physiochemical changes in WT and the am mutant after waterlogging treatments

Under normal growth conditions, the chlorophyll a, chlorophyll b and total chlorophyll contents of the am seedlings leaves were lower than those of the WT. After 12 days of waterlogging stress, the contents of chlorophyll a, chlorophyll b and total chlorophyll in the leaves of WT seedlings decreased by 7.46%, 7.82% and 11.55%, respectively. The changes were more marked in the mutant, with the chlorophyll a, chlorophyll b and total chlorophyll contents decreasing by 14.55%, 14.52% and 25.31%, respectively, which was approximately twice as much as in WT seedlings (Fig. 2a–c). Anthocyanins play important roles in alleviating plant damage caused by environmental stresses. Under both normal water supply and waterlogging-stress conditions, the anthocyanin contents in leaves of the am mutant were greater than in WT. Accumulations of anthocyanin in both materials increased by different degrees after waterlogging. However, the amplitude of the anthocyanin increase in am seedlings was greater than that of WT seedlings, being fivefold in the mutant and twofold in the WT (Fig. 2d).

Fig. 2
figure 2

Waterlogging stress had different effects on the physiological indicators of WT and the am mutant. af The chlorophyll a (a), chlorophyll b (b), total chlorophyll (c), anthocyanin (d), MDA (e) and soluble sugar (f) contents were determined in the WT and am mutant after 12 d of waterlogging stress. Error bars represent the SDs from three replicated experiments. Significant differences between the am mutant and WT with and without (CK) waterlogging stress are indicated (Student’s t-test) as follows: ***P < 0.001; **P < 0.01; *P < 0.05

Plants trigger or intensify membrane lipid peroxidation in response to environmental stresses. Malondialdehyde (MDA) is one of the products, and it is used as an indicator of lipid peroxidation, indicating the degree of cell membrane peroxidation and the intensity of plant responses to stress [49]. As shown in Fig. 2e, the MDA contents in seedling leaves remained stable in non-waterlogged rapeseed plants, and there was no significant difference between the WT and the mutant. Waterlogging stress treatments triggered a marked increase in the MDA content, which increased by 4.5-fold in WT and eightfold in the mutant. Furthermore, when a plant encounters adversity, soluble sugars accumulate in vivo to regulate the osmotic capacity of the plant and improve its adaptability to stress. After 12 days of waterlogging stress, the soluble sugar levels in seedling leaves of the mutant increased more significantly than in the WT (Fig. 2f). These results indicated that the physiological metabolism of rapeseed seedlings was seriously affected by waterlogging stress, and the mutants may alleviate the corresponding damage by regulating the chlorophyll, anthocyanins, MDA and osmotic adjustment substance levels.

Transcriptome analysis of WT and the am mutant in response to waterlogging stress as assessed by RNA-seq

To further understand the molecular mechanisms underlying the different responses of the two materials to waterlogging stress, transcriptomes of seeding leaves of WT ZS11 and the am mutant after a 24-h waterlogging period followed by 8 days of recovery growth were investigated. Corresponding seedlings without waterlogging at the germination stage were used for controls in the screening of differentially expressed genes (DEGs) by RNA-seq. After data quality controls and data filtering, there were 26,187,420 (ZS11CK), 26,777,456 (ZS11WL), 27,536,936 (amCK) and 22,117,491 (amWL) 150-bp paired-end reads, with an average of 25,654,825 reads per sample. The average clean base reached 7.69 GB per sample, and approximately 83% of the reads from the cDNA libraries were mapped to the B. napus reference genome (Table 2). The Q30 percentage was greater than 95%, and the randomness of the sample sequencing was good, which ensures the sequencing quality (Additional file 1: Fig. S2). With an increase in sequencing quantity (number of reads), the number of genes detected also increases. We found that the numbers of genes detected in the treatment and control groups of the two varieties were close to saturation, within a certain Reads Per Kilo bases per Million reads (RPKM) interval, indicating that the sequencing quantity met the requirements (Additional file 1: Fig. S3). A correlation analysis showed that the correlation coefficients of gene expression levels between control and treatment groups of each material were close to 1, whereas the gene expression patterns between WT and mutant samples were different, suggesting that the reproducibility among samples was high and that sequencing data were reliable for further analyses (Additional file 1: Fig. S4).

Table 2 Statistics for sequencing reads after filtering and genomes alignment

A DEG analysis was conducted to detect global expression changes in the two materials during the waterlogging process. To visually display the distributions of false discovery rate (FDR) and fold change (FC) values of all the genes between control and treatment groups of each variety, MA (minus versus add) and volcano plots of both Z11CK-vs-Z11WL and amCK-vs-amWL were constructed. As shown in Additional file 1: Fig. S5, most of the genes concentrated in the black central region and showed no significant differences in expression. The DEGs (Red dots) were more concentrated in amCK-vs-amWL, indicating that there were more genes with significantly different expression levels in the mutant. Ultimately, 1370 (856 upregulated and 514 downregulated) genes and 2,336 (1145 upregulated and 1221 downregulated) DEGs were found in the WT and mutant, respectively (Fig. 3a, Additional file 2: Tables S1, S2). Among them, there were 208 (5.95%) common DEGs in the two varieties, and the others (94.05%) were DEGS specific to only one variety. Moreover, the number of specific DEGs in the mutant was approximately twice that in the WT (Fig. 3b). To reflect the gene expression levels and patterns of multiple samples, a hierarchical cluster analysis of the DEGs in the two varieties was performed based on the fragments per kilobase of transcript per million bases (FPKM) values, which provided an overall understanding of the transcriptional changes (Fig. 3c). Using co-expression clustering, 2296 known DEGs identified in the mutant were grouped into five subclusters and 1321 known DEGs identified in the WT were grouped into six subclusters on the basis of log10 (FPKM + 1) values (Additional file 1: Fig. S6, Additional file 2: Tables S3, S4). These results indicated that the global gene expression profiles of the two varieties varied greatly, and there were more transcriptomic variations in the mutant than in WT after exposure to waterlogging stress.

Fig. 3
figure 3

The expression profiles of waterlogging-regulated DEGs in the WT and am mutant. a Column diagram representing the numbers of DEGs in the WT and am mutant after the waterlogging treatment. b Venn diagrams showing DEGs and the overlaps of sets obtained across comparisons. c Heatmap clustering showing expression patterns of waterlogging-regulated DEGs based on FPKM values in the two different rapeseed cultivars

Functional classification and enrichment analysis of DEGs

To acknowledge the putative functions and pathways associated with the waterlogging tolerance phenotype of the am mutant, Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs were performed. For the GO enrichment analysis, the identified DEGs annotated with GO terms were assigned to three main functional categories, cellular component, biological process and molecular function, and then were divided into 44 and 47 sub-groups in ZS11CK-vs-ZS11WL and in amCK-vs-amWL, respectively (Fig. 4a, b). Many DEGs were assigned to more than one sub-group, and the first three enriched sub-groups in each of the three main categories were the same in ZS11CK-vs-ZS11WL and amCK-vs-amWL. Most of the common DEGs in the WT ZS11 and the am mutant were assigned to metabolism process, binding, catalytic activity, cellular process, cell, cell part and membrane. A further analysis showed that DEGs in the mutant were more enriched in response to endogenous stimulation and hormone, oxidoreductase activity, ion channel inhibitor activity and ion channel regulator activity (Additional file 2: Table S5), whereas the WT DEGs were mainly concentrated in isocitrate lyase activity, prochlorophyllate reductase activity, peptidase regulator activity, lipid transport, 2-acylglycerol o-acyltransferase activity, glyoxylic acid cycle and lipid localization (Additional file 2: Table S6).

Fig. 4
figure 4

Functional classifications of DEGs based on the a, b GO annotation and c, d KEGG pathway enrichments in the WT and am mutant

To further determine the functional distributions of these DEGs and supplemental GO classifications, we also conducted a gene functional annotation using the KEGG database. The KEGG enrichment analysis showed that the DEGs were divided into five major categories, cellular processes, environmental information processing, genetic information processing, metabolism and organismal systems. More DEGs were highlighted in ‘metabolism’ than in the other functional categories. A total of 2483 DEGs were mapped to 33 KEGG annotation pathways. As shown in Fig. 4c, the most significantly enriched pathway in WT (ZS11CK-vs-ZS11WL) was carbohydrate metabolism in the metabolic pathway branch, which contained 104 DEGs, followed by pathways related to signal transduction (83) in the environmental information processing branch and pathways related to transport and catabolism (79) in the cellular process branch (Additional file 2: Table S7). In contrast, the DEGs of amCK-vs-amWL were most significantly enriched in signal transduction (176) in the environmental information processing branch, followed by translation (115) in the genetic information processing branch and lipid metabolism (107), carbohydrate metabolism (106) and amino acid metabolism (95) in the metabolic pathway branch (Fig. 4d, Additional file 2: Table S8).

Detection of DEGs responsible for anthocyanin biosynthesis under waterlogging-stress conditions

To elucidate the differences in waterlogging stress between the two varieties at the transcriptional level, anthocyanin biosynthetic pathways were reconstructed, and structural genes including the early and late biosynthesis genes involved in anthocyanin synthesis were indicated [37] (Fig. 5a). In this study, DEGs associated with anthocyanin production were identified and presented in Fig. 5b using a log2FC(WL/CK) relative transcript changes. A total of 9 genes (7 upregulated and 2 downregulated) in WT and 15 genes (9 upregulated and 6 downregulated) in the mutant were identified as anthocyanin biosynthesis-related structural genes (Fig. 5b). The expression changes in these DEGs might be related to the different colors of the two varieties under waterlogging stress. Most synthetic genes were different members of the same gene family and showed different expression patterns between the two varieties. Additionally, some family members exhibited similar expression patterns between the two varieties, but the range of the expression changes was more marked in the mutant than in the WT after waterlogging. For example, phenylalanine is a key amino acid in anthocyanin biosynthesis, and 21 DEGs involved in phenylalanine biosynthesis and related metabolic pathways, belonging to nine gene families, were identified, and most of them were upregulated in the mutant compared with in the WT (Additional file 2: Table S9). One phenylalanine ammonia-lyase (PAL) gene (BnaC04T0337700ZS) displayed similar upregulated expression profiles in the two varieties, but the increased expression level was more pronounced in the mutant than in the WT. Dihydroflavonols are common substrates of flavonol [via flavonol synthase (FLS)] and anthocyanin (via DFR) synthesis, and flavonols and anthocyanins accumulation require competition for dihydroflavonols. In our study, one FLS (BnaC02T0021300ZS) gene was significantly downregulated only in the mutant. Meanwhile, two DFR genes (BnaC09T0215200ZS and BnaA09T0187400ZS) displayed similar upregulation expression profiles in the mutant. Additionally, one anthocyanidin reductase (ANR) gene (BnaA03T0496700ZS) that functions in proanthocyanin production and contributes to the normal green color of plants was only upregulated in the WT. These results suggested that the am mutant tended to preferentially synthesize anthocyanins compared with flavonols and proanthocyanins under waterlogging-stress conditions. Leucoanthocyanidins catalyzed by DFR are further catalyzed by leucoanthocyanidin dioxygenase (LDOX) to form anthocyanidins. Although no LDOXs were detected among the DEGs, UDP-flavonoid glucosyltransferases (UFGTs), which convert anthocyanidins into anthocyanins, were significantly upregulated in the mutant. This may contribute to the increased anthocyanin content and strong waterlogging resistance in the mutant compared with the WT (Fig. 5b).

Fig. 5
figure 5

Identification of DEGs involved in anthocyanin biosynthesis in the two different rapeseed cultivars. a Schematic diagram of anthocyanin biosynthetic pathways and key enzymes involved in the pathways. Anthocyanin early biosynthesis-related structural genes includes genes encoding phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), 4-coumaroyl:CoA ligase (4CL), chalcone synthase (CHS), chalcone isomerase (CHI), flavanone 3-hydroxylase (F3H) and flavonoid 3′-hydroxylase (F3′H). Anthocyanin late biosynthesis-related structural genes includes dihydroflavonol 4-reductase (DFR), leucoanthocyanidin dioxygenase (LDOX), UDP-flavonoid 3-O-glucosyltransferase (UF3GT) and anthocyanin acyltransferase (AAT). b Expression profiles of DEGs related to anthocyanin biosynthetic structural genes. Transcript level changes were normalized using log2 FC transformed counts. The blue and yellow represent up- and downregulated expression, respectively. The red and green words represent DEGs and their FCs in the mutant and WT, respectively

Detection of differentially expressed transcription factors under waterlogging-stress conditions

The MBW complex is a well-known transcriptional regulator of anthocyanin synthesis in plants, but the other transcription factors (TFs) involved in anthocyanin synthesis are still unclear. In our study, in addition to TF members belonging to the MYB-type and bHLH families, we identified many common TF families that showed significantly changed expression levels between the WT and the mutant, such as WRKY, NAC, C2H2 and bZIP. They play an important roles in responses to biotic and abiotic stresses. Among the 3,528 DEGs, 79 genes (54 upregulated and 25 downregulated) in WT and 196 genes (110 upregulated and 86 downregulated) in the mutant were identified as TFs under waterlogging-stress conditions (Additional file 2: Tables S10, S11). They all belonged to 43 TF families (Fig. 6). In the WT, most of the upregulated TFs belonged to the bHLH (11), AP2-EREBP (10), C2C2-GATA (5) and G2-like (5) families, and the downregulated TFs mostly belonged to AP2/ERFBP (3), ABI3VP1 (2), G2-like (3), WRKY (2) and some MYB families (2) (Fig. 6a). In the am mutant, most of the upregulated TFs were classified into MYB (12), AP2/EREBP (12), zf-HD (8), ABI3VP1 (6) and bHLH (6) families, whereas the downregulated TFs belonged to the NAC (11), bHLH (8) and WRKY (8) families (Fig. 6b). Among them, ABI3VP1 is a family of plant-specific TFs that contains a highly conserved B3 functional domain and has a potential role in regulating seed development to germination [50]. Moreover, some TFs showed line-specific expression changes in one of the two materials, and there were more specific TFs in the mutant (52) than in the WT (3). It is noticeable that ten TF families, including Alfin-like, C2C2-Dof, C2C2-YABBY, GeBP, PLATZ, Sigma 70-like, SRS, TCP, TIG and Trihelix, showed line-specific upregulated expression changes in the mutant. These families mainly function in biological processes of plant morphogenesis, plant growth and development, photosynthesis, hormone response and stress response, whereas three TF families, one upregulated TF family VOZ and two downregulated TF families C2C2-CO-like and CSD, showed line-specific expression changes in WT. Among them, the VOZ TFs positively regulate abiotic stress tolerance [51, 52], plant defense response [53] and flowering [54, 55], and they negatively regulate PhyB-mediated seed germination in Arabidopsis [56].

Fig. 6
figure 6

The distributions of differentially expressed TFs under waterlogging stress in a WT and b the mutant. The numbers in the histograms represent the number of up- or downregulated TFs in each TF family

Detection of DEGs associated with plant hormone signal transduction

Plant endogenous hormones, such as auxin, cytokinin, ABA, jasmonic acid (JA) and BR are not only growth regulators in higher plants, but they also participate in a wide range of physiological functions, such as plant responses to stress and anthocyanin synthesis [57]. There were 87 genes enriched in the term of “plant hormone signal transduction” in the mutant, including 24 upregulated DEGs and 63 downregulated DEGs (Fig. 7, Additional file 2: Table S12). Most of the upregulated DEGs belonged to the ET and auxin signaling pathways. Among them, ten DEGs involved in ET regulation, mainly including ET-responsive TFs (ERFs) and ET-response sensors, were upregulated in the mutant. A total of seven auxin signaling-related DEGs, which encoded three key proteins, auxin-responsive proteins, indole-3-acetic acid-amido synthetases and TRANSPORT INHIBITOR RESPONSE 1, were also upregulated in the mutant. Therefore, these two pathways might play important roles in anthocyanin accumulation and B. napus responses to waterlogging stress. Moreover, JA (3), SA (2) and BR (1) pathway-related genes were also identified among the upregulated DEGs in the mutant (Fig. 7a). Most downregulated genes in the mutant are homologs of auxin-responsive proteins SAUR (small auxin upregulated RNA) (41) and IAA (indole-3-acetic acid) (9), and auxin-induced proteins (3) (Fig. 7b). Only 23 DEGs in the WT belonged to “plant signal transduction”, including 13 upregulated DEGs and ten downregulated DEGs (Additional file 2: Table S12). Unlike in the mutant, the upregulated DEGs in the WT mainly encoded ABA receptor (4) and auxin transporter (5) homologs (Fig. 7c). Auxin (4), ABA (4), JA (1) and ET (1) signaling pathway-related DEGs were downregulated in the WT (Fig. 7d).

Fig. 7
figure 7

The distributions of DEGs related to plant hormone signal transduction pathways under waterlogging-stress conditions in a, b the mutant and c, d WT

Quantitative real-time PCR (qRT-PCR) validation of DEG expression levels

To confirm the differential expression levels of DEGs as assessed by RNA-seq, eight key TF genes associated with the regulation of anthocyanin accumulation, five anthocyanin biosynthetic genes and three genes related to hormone signal transduction were validated by qRT-PCR (Fig. 8). The expression levels of one TF, GRF3 (BnaA04T0234500ZS), and two structural genes, TCMO (BnaA03T0147300ZS) and DFR (BnaC09T0215200ZS), were all clearly upregulated in the mutant, whereas three key TFs, FLC (BnaC02T0039100ZS), NA102 (BnaA02T0405700ZS) and WRKY31 (BnaA01T0121300ZS), two structural genes, FLS (BnaC02T0021300ZS) and CAMT1 (BnaC06T0343500ZS), and two hormone-related genes JAR1 (BnaC03T0263700ZS) and A10A5 (BnaA01T0002400ZS), were all significantly downregulated in the mutant. The expression patterns of these genes were in good agreement with the RNA-seq trends. The remaining six genes, COL8 (BnaC05T0309900ZS), NAC13 (BnaC09T0521700ZS), CSP4 (BnaA09T0601800ZS), PHLD (Bnascaffold0163T0000200ZS), ANR (BnaA03T0496700ZS) and PYL4 (BnaA04T0246400ZS) were differentially expressed in the WT, and also showed the same trends as the transcriptome analysis. Therefore, the qRT-PCR results were consistent with the RNA-seq data, suggesting the reliability of the transcriptome sequencing.

Fig. 8
figure 8

The qRT-PCR analysis of selected DEGs in B. napus leaves of the two different rapeseed cultivars during the waterlogging treatment. Data are presented as means ± SDs from three independent experiments. Student’s t-test was used for the statistical analysis between treated and untreated (control) samples in each rapeseed cultivar (*P < 0.05; **P < 0.01; ***P < 0.001)

Discussion

Waterlogging stress affects phenotypic and physiological changes in the am mutant

As important water-soluble natural pigments in plants, anthocyanins not only play key roles in plant pollination and seed transmission, but also function in resisting abiotic and biotic stresses. There are some natural mutants that allow plants without anthocyanins to have the ability to synthesize and accumulate anthocyanins. Studying these mutants can increase our understanding of the regulatory mechanisms related to anthocyanin synthesis and metabolism. Waterlogging stress affects rapeseed germination, emergence and seedling growth, causing not only phenotypic changes, but also affecting the physiological and metabolic activities. In this study, the relative seedling survival rate and root lengths of the am mutant were both greater than those of WT at different times after waterlogging stress (12, 24 and 36 h) (Fig. 1). Moreover, the MDA and soluble sugar contents increased in the WT and mutant after waterlogging stress, whereas significant increases/change ranges in the mutant were observed compared with the WT (Fig. 2). MDA can reflect the degree of membrane damage caused by stress, and the accumulations of osmoregulatory substances, such as soluble sugars, can improve plant tolerance to stress by increasing osmotic potential and maintaining membrane stability [4, 58]. Our results showed significant increases in the MDA and soluble sugar contents of the am mutant, which were in accordance with previous studies concerning waterlogging-resistance in waterlogging-tolerant species, such as Medicago sativa [59], Zea mays L. [10] and Arachis hypogaea L. [60]. Thus, the mutants may reduce the damage caused by waterlogging stress by regulating the anthocyanins, MDA and osmotic adjustment substance contents. Phenotypic changes were mainly reflected by the leaves, stems and petioles of the am mutant changing to purple earlier than in the WT (Fig. 1c).

Comparative transcriptome analysis of the two B. napus lines using RNA-seq

Because the waterlogging tolerance of rapeseed is controlled by multiple genes, omics technologies can be used to more comprehensively explore waterlogging tolerance-related proteins/genes and analyze the molecular mechanisms of waterlogging sensitivity. For example, a comparative proteomic approach revealed that the response proteins are mainly enriched in the process of redox processes, protein phosphorylation, signal transduction and metabolism, and the gene expression changes in tolerant rapeseed cultivar occur faster than in the sensitive cultivar after exposure to waterlogging treatments [8]. Transcriptome analysis showed that gene expression profiles, representing various pathways involving multiple levels of regulation, of waterlogging-tolerant rapeseed cultivar changed in response to waterlogging stress [61]. These findings provide important clues for understanding the waterlogging-tolerance mechanisms in this important economic crop. However, to date, research in this field is quite limited, which has hindered the identification of key genes that function in waterlogging tolerance and in the molecular breeding of rapeseed having a high and stable yield.

In this study, transcriptome sequencing technology was used to analyze the response of WT ZS11 and the am mutant to waterlogging stress at the germination stage. Many genes related to germination and growth under waterlogging stress were identified by analyzing DEGs, and these increase our understanding of the response mechanisms of rapeseed to waterlogging stress. The GO enrichment analysis revealed that the DEGs were mainly involved in metabolic process, cellular process, single biological process, biological regulation, regulation of biological process and response to stimulation (Fig. 4a). The pathway enrichment analysis showed that the DEGs were involved in signal transduction, translation process, lipid metabolism, carbohydrate metabolism and amino acid metabolism, which were the main characteristics of waterlogging germination (Fig. 4b). The transcriptome data also showed that in the am mutant, the most DEGs were related to signal transduction, followed by metabolism-related genes (Fig. 4a), suggesting that the mutants were more likely to use limited resources to initiate defense responses under waterlogging-stress conditions.

DEGs related to the anthocyanin accumulation

In Arabidopsis, anthocyanin biosynthesis-related structural genes are divided into early and late biosynthesis genes [37]. After waterlogging stress, the mutant was more sensitive to environmental changes, and the anthocyanin content increased by approximately fivefold, which may explain the difference in seedling color between the two lines. Moreover, the changes in seedling color and anthocyanin content were related to the differential expression of structural genes involved in anthocyanin synthesis, including the early biosynthesis genes PAL, C4H, 4CL, CHS and F3H, as well as the late biosynthesis genes DFR and UFGT (Fig. 5b), as determined by the transcriptome data. Under natural conditions and waterlogging stress, the expression levels of these key structural genes in the am mutant were higher than in WT, which indicated that the anthocyanin biosynthetic pathway was involved in the responses of plants to waterlogging stress, and the mutated genes in the am mutant may function by affecting anthocyanin anabolic metabolism. Both dihydrokaempferol and dihydroquercetin are dihydroflavonols, which represent the intermediate products of anthocyanin synthesis and can be catalyzed by FLS to form flavonoids. There are two ways to increase anthocyanin synthesis in plants. One involves increasing the accumulation of dihydroflavonols to promote anthocyanin accumulation, and the other is to decrease the FLS activity and increase the DFR activity, thereby reducing flavonoid accumulation and increasing anthocyanin the synthesis [62, 63] (Fig. 5a). In our study, the significant downregulation of FLS and significant upregulation of DFR were detected in mutants. Therefore, the mechanism leading to different anthocyanin accumulations between the two varieties may be the latter. Similarly, some varieties prioritize flavonol biosynthesis by enhancing FLS activity and strongly reduce anthocyanin accumulation by down regulating the activities of enzymes such as DFR and UFGT [64, 65]. Moreover, UFGT can convert anthocyanidins into stable anthocyanins through glycosylation. The differential expression of UFGT in the two varieties may also be related to the high accumulation of purple anthocyanins in the mutant.

DEGs related to the regulation of anthocyanin synthesis and accumulation

The synthesis and regulatory network of plant anthocyanins is very complex, and many biosynthetic regulators of anthocyanins have been identified. The expression of genes involved in anthocyanin biosynthetic pathway is regulated by MBW complexes, which mainly regulate the expression levels of different structural genes at the transcriptional and post transcriptional levels [66, 67]. According to the results of transcriptome analysis, members of the MYB and bHLH TF families were upregulated in the two varieties, but we did not detect any WD40-repeat proteins differentially expressed between the two varieties under waterlogging conditions. In some species, such as apple and turnip, the ternary MBW complex is not necessary for the regulation of anthocyanin synthetic genes [68, 69]. Therefore, we speculated that the MYB plus bHLH dimeric complex is enough to regulate the expression of anthocyanin synthetic genes in B. napus. Among the differentially expressed TF genes, some belonged to the AP2/ERF TF family, which contains four subfamilies, AP2, RAV, ERF and DREB, showed specifically upregulated expression changes in the mutant. Based on our transcriptome data, only members of the AP2 and ERF subfamilies were identified as being differentially expressed, suggesting that they play a key role in the responses of rapeseed to waterlogging stress. Many studies have demonstrated that members of the AP2 and ERF subfamilies are not only key factors in stress signaling pathway [70, 71], but are also involved in plant seed development, as well as high salt and waterlogging stress responses [72,73,74]. The TFs that showed line-specific expression changes in the mutant, Alfin-like gene family, C2C2-Dof, C2C2-YABBY, GeBP, LOB, PLATZ, SBP, TCP, Tify, WRKY and zf-HD, which are unique to plants, played important roles in regulating plant growth and development and adversity stress responses. ABI3VP1 is also a family of plant-specific TFs that contain a highly conserved B3 functional domain, and it has a potential role in regulating seed development to germination [50, 75]. Other common TF families, such as WRKY, C2H2, NAC, ARF and bZIP, involved in plant responses to adversity, were also identified in our study. Notably, some members of TF families, including WRKY, showed inconsistent expression trends in the same material, indicating that different family members may play different roles after waterlogging stress. The WRKY TF family is involved in plant seed germination, hormone signal transduction, and the regulation of physiological and biochemical processes under stress conditions [76, 77]. Moreover, the expression patterns of these TFs genes were validated by performing qRT-PCR, which confirmed the reliability of transcriptome sequencing (Fig. 8). These TFs can also be used as targets for future functional identification.

Additionally, plant hormones are also important factors modulating anthocyanin accumulation. Several genetic mutants have been isolated, and they are mainly used to clarify the effects of hormones on anthocyanin biosynthesis. For example, the high level of ABA in ein2 mutants may be the reason for the high anthocyanin accumulation levels in these plants [78]. Plant hormones, such as cytokinin, ABA, JA and BR play positive regulatory roles [42, 68, 79, 80], and auxin, GA and melatonin play negative roles in the synthesis of anthocyanins [81,82,83]. GA signals inhibit anthocyanin biosynthesis through DELLA proteins, which mediate anthocyanin biosynthesis in plant responses to abiotic stress [82], whereas JA signals participate in anthocyanin biosynthesis regulated by DELLA and restore some anthocyanin accumulation inhibited by GA [31]. In agreement with previous research results, in addition to the GA pathway, DEGs involved in auxin ABA, BR, ET and JA signaling pathways were identified to be upregulated in the mutant, and DEGs associated with the auxin and ET pathways accounted for the greatest proportions, indicating that these two pathways played key roles in response to waterlogging stress and anthocyanin accumulation in rapeseed. In addition, we also detected DEGs encoding “pathogenesis-related protein 1” and “regulatory protein NPR5” belonging to SA signaling, which showed especially upregulated expression in the mutant. Spraying SA at the early stage of grape ripening can greatly improve the anthocyanin content in fruit [84]. However, the detailed mechanisms of SA regulation in anthocyanin accumulation are still unclear and need further study.

Screening of new transcripts after matching reference genome

After combining the genome mapping results of all the samples and removing the redundancy, the transcriptome was assembled by using the software of cufflinks, and then, the assembled transcripts were compared with the known gene annotations using cuffcompare software to obtain the classification of a reconstructed transcriptome and predict the new transcripts. An analysis of these new transcripts showed that they were mainly related to plant anthocyanin-induced response, kinase activity, signal transduction, hypersensitivity response, defense response and energy metabolism (Additional file 2: Table S13). The production of these new transcripts is an adaptive manner of gene evolution in rapeseed. Zou et al. [85] found that 11 transcripts containing early termination codons in waterlogging-tolerant varieties have specific regulatory functions under waterlogging-stress conditions. A more in-depth analysis of DEGs and the newly generated transcripts may help determine why the reason the mutants are tolerant to waterlogging stress.

Conclusions

A natural mutant that easily induces anthocyanin production was screened from strain ZS11 of B. napus. Compared with the WT, the mutant showed increased tolerance to waterlogging stress during germination. A transcriptome sequencing analysis showed that the DEGs in the mutants were mainly related to signal transduction, translation and metabolism. Among them, some key genes involved in anthocyanin accumulation and hormone signal transduction, and TF genes involved in stress were differentially expressed in mutants. In the future, we will establish a mapping population for gene localization and genetic analyses, and reduce the range of candidate genes in combination with the transcriptome data. The identification and cloning of the mutant gene will contribute to the breeding of waterlogging-tolerant varieties of rapeseed.

Materials and methods

Plant materials

Two B. napus varieties WT ZS11 and its am mutant, having different seedling colors, from the seed bank of the School of Life Sciences at Jiangsu University were used as plant materials. The plants were grown in a plant growth room in conditions described previously [27]. Fresh leaves were harvested after a 24-h waterlogging treatment followed by an 8-day recovery period [13, 61], and corresponding untreated seedling leaves were collected for the RNA-seq and qRT-PCR experiments. Each sample had three biological replicates.

Waterlogging treatment and phenotypic characteristics

The WT and mutant seeds were disinfected with 0.5% sodium hypochlorite solution, washed with ddH2O, evenly placed in a culture dish, having a diameter of 9 cm, and covered with three layers of wet filter paper. They were then germinated in an incubator for 36 h at 22 °C and 70% relative humidity with a 16-h light time. The germinated seeds (radicle length approximately 2–5 mm) of the treatment group were placed in a 10-mL centrifuge tube and subjected to waterlogging for 12 h, 24 h and 36 h. After washing with distilled water three times, they were transferred to a medium prepared with Hoagland’s nutrient solution for germination and growth for 8 days. The control group was not treated and was cultured normally. After 8 days, the survival number and seedling rate were determined. Then, ten seedlings were randomly selected and the root lengths, stem lengths and fresh weights measured. The experiment was repeated three times. The relative seedling rate, relative root length, relative stem length, relative fresh weight and relative vigor index of each seedling in both the treatment and control groups were calculated, and the relative vigor index was used as the main judgment index.

Physiological measurement

The WT and mutant seeds germinated in the petri dish, and then the emerged germinated seeds were transplanted into a plastic flowerpot containing nutrient soil to continue growing. When the plants grew to the six-leaf one-heart stage, WT and mutant seedlings exhibiting consistent growth were selected for waterlogging treatments. The flowerpots for planting rapeseed plants were placed in a plastic box containing water to keep the water surface in the plastic box from overflowing the soil layer, and the control group was provided a normal water supply. After waterlogging for 12 days, the leaves were harvested and the related physiological indexes were measured. The chlorophyll content was determined using ethanol extraction colorimetry as described previously [48]. Under different pH conditions, anthocyanins present different colors, and the color depth is proportional to the anthocyanin content. Consequently, the anthocyanin content was determined using colorimetry in accordance with Cao et al. [86]. The MDA contents were determined using the thiobarbituric acid method at 532 nm [87]. The soluble sugar content was determined using the sulfuric acid anthrone colorimetric method at 624 nm [88]. The Coomassie Brilliant Blue method was used to quantify the soluble protein contents in plant crude extracts. The experiment was repeated three times.

RNA extraction, library preparation, RNA sequencing and data analyses

The normally germinated WT and mutant seeds were divided into two groups, respectively. The treatment groups were subjected to waterlogging for 24 h and were recovery growth for 8 days on Hoagland medium. The control group was not treated and was cultured normally. The total RNA was extracted from the young leaves tissues of WT Z11 and the am mutant in control and treatment groups using an RNA extraction kit (Omega, Norcross, CA, USA) following the manufacturer’s instructions. After the RNA stock solution was diluted, the quality of each RNA sample was detected using an Agilent Technologies 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The OD260/280 of each RNA sample was between 1.9 and 2.1, indicating that the purity levels of the extracted RNA samples were high. After extraction, electrophoretic detection was performed to observe the electrophoretic band types. Equal amounts of RNA samples from three different plants per variety were pooled. The total RNA sample quality was determined by the Wuhan FISA Gene Information Co., Ltd., and then, the cDNA library was constructed and sequenced on an Illumina Hiseq X Ten platform (Illumina Inc., San Diego, CA, USA). The raw RNA-seq data were submitted to the National Center for Biotechnology Information (NCBI) Sequence Read Archive with a Bioproject ID: PRJNA821348.

The data analysis after sequencing mainly included sequencing data quality control filtering, sequencing quality evaluation, reference sequence alignment, overall RNA-seq quality evaluation, sample correlation analysis and differential expression analysis. For the original sequencing data, reads with adaptor sequences, N (uncertain base) content > 10% and low-quality base (Q ≤ 20) content > 50% were filtered to get clean data using trimmatic software (v0.33) [89]. The statistics of sequencing data showed that Q20 > 95% and Q30 > 89%, indicating that the sequencing quality was good and could be used for subsequent comparison and data analysis. After quality control, calling by Tophat2 (v2.1.1, parameter: -library-type fr-unstranded -G) [90] and bowtie2 (v2.2.2, parameter: default) [91] was used to compare the sequence of each sample to NCBI EST database of B. napus and ZS11 B. napus reference genome (http://cbi.hzau.edu.cn/cgi-bin/bnapus//blast), and then, the alignment results were statistically analyzed by using the bowtie2 function in RSEM (v1.3.0, parameter: default) [92]. The FPKM values were used to calculate the number of aligned reads and normalize the expression levels of genes and transcripts. EdgeR (v3.6.8) package method [93] was used for screening DEGs using the thresholds of false discovery rate ≤ 0.05 and log2FC (fold change (WL / CK) for a gene) > 1 or log2FC < − 1 between untreated or waterlogging-treated groups in each variety.

Functional and enrichment analysis of DEGs

The functions of DEGs were identified using GO (http://www.geneontology.org/) and KEGG (http://www.genome.jp/kegg/) databases. The GO enrichment was performed using GOseq, and DEGs with a corrected p value ≤ 0.05 were considered as significantly enriched GO terms [94]. The KEGG pathway enrichment analysis was perform using KOBAS3.0, and pathways with a corrected p value ≤ 0.05 were defined as significantly enriched pathways [95].

Quantitative real-time PCR

Sub-samples (total RNA of 500 ng) of RNA-seq were treated with gDNA wiper Mix to remove genomic DNA, and then, cDNA was synthesized using Hiscript®II QRT SuperMix Kit (Vazyme, Nanjing, China) in accordance with the manufacturer’s protocols. The qRT-PCR was conducted using AceQ®qPCR SYBR®Green Master Mix (Vazyme, Nanjing, China) in accordance with the manufacturer’s instructions and run on an ABI 7300 Real Time PCR System (Applied Biosystems, Carlsbad, CA, USA). The PCR program was as follows: 95 °C for 5 min and 40 cycles of 95 °C for 10 s, followed by 60 °C (± 5 °C) for 31 s. Primer 5.0 was used to design primers, the specificity of the primer sequences was tested using a Primer-BLAST of the NCBI database, and the reasonably designed primer sequences were sent to Sangong Bioengineering (Shanghai, China) Co., Ltd. for synthesis and purification (Additional file 2: Table S14). The ACTIN gene was used as the internal control. The application of the 2[–ΔΔC(t)] method converts the threshold cycle value output by the instrument into the relative gene expression level [96]. Each reaction included two biological replications with three technical replications.

Availability of data and materials

All data generated or analyzed during this study are included in the manuscript and additional files.

Abbreviations

DEGs:

Differently expressed genes

TF:

Transcription factor

WT:

Wild type

ABA:

Abscisic acid

ET:

Ethylene

BR:

Brassinolide

CHS:

Chalcone synthase

DFR:

Dihydroflavonol 4-reductase

ANS:

Anthocyanidin synthase

bHLH:

Basic helix–loop–helix

MDA:

Malondialdehyde

FDR:

False discovery rate

FC:

Fold change

FPKM:

Fragments per kilobase of transcript per million bases

GO:

Gene ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

PAL:

Phenylalanine ammonia-lyase

C4H:

Cinnamate 4-hydroxylase

4CL:

4-Coumaroyl:CoA ligase

CHI:

Chalcone isomerase

F3H:

Flavanone 3-hydroxylase

F3′H:

Flavonoid 3′-hydroxylase

LDOX:

Leucoanthocyanidin dioxygenase

UF3GT:

UDP-flavonoid 3-O-glucosyltransferase

AAT:

Anthocyanin acyltransferase

qRT-PCR:

Quantitative real-time PCR

References

  1. Kaur H, Wang L, Stawniak N, Sloan R, van Erp H, Eastmond P, Bancroft I. The impact of reducing fatty acid desaturation on the composition and thermal stability of rapeseed oil. Plant Biotechnol J. 2020;18:983–91.

    Article  CAS  PubMed  Google Scholar 

  2. Wang X, Zheng M, Liu H, Zhang L, Chen F, Zhang W, et al. Fine-mapping and transcriptome analysis of a candidate gene controlling plant height in Brassica napus L. Biotechnol Biofuels. 2020;13:42.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Pullen J, Saeed K. Investigation of the factors affecting the progress of base-catalyzed transesterification of rapeseed oil to biodiesel FAME. Fuel Process Technol. 2015;130:127–35.

    Article  CAS  Google Scholar 

  4. Ding LN, Li T, Guo XJ, Li M, Liu XY, Cao J, et al. Sclerotinia stem rot resistance in rapeseed: recent progress and future prospects. J Agric Food Chem. 2021;69:2965–78.

    Article  CAS  PubMed  Google Scholar 

  5. Ploschuk RA, Miralles DJ, Colmer TD, Ploschuk EL, Striker GG. Waterlogging of winter crops at early and late stages: impacts on leaf physiology, growth and yield. Front Plant Sci. 2018;9:1863.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Leul M, Zhou WJ. Alleviation of waterlogging damage in winter rape by uniconazole application: effects on enzyme activity, lipid peroxidation, and membrane integrity. J Plant Growth Regul. 1999;18:9–14.

    Article  CAS  PubMed  Google Scholar 

  7. Lee YH, Kim KS, Jang YS, Hwang JH, Lee DH, Choi IH. Global gene expression responses to waterlogging in leaves of rape seedlings. Plant Cell Rep. 2014;33:289–99.

    Article  CAS  PubMed  Google Scholar 

  8. Xu J, Qiao X, Tian Z, Zhang X, Zou X, Cheng Y, et al. Proteomic analysis of rapeseed root response to waterlogging stress. Plants. 2018;7:71.

    Article  CAS  PubMed Central  Google Scholar 

  9. Herzog M, Striker GG, Colmer TD, Pedersen O. Mechanisms of waterlogging tolerance in wheat—a review of root and shoot physiology. Plant Cell Environ. 2016;39:1068–86.

    Article  CAS  PubMed  Google Scholar 

  10. Ren B, Zhang J, Dong S, Liu P, Zhao B. Effects of waterlogging on leaf mesophyll cell ultrastructure and photosynthetic characteristics of summer maize. PLoS ONE. 2016;11: e0161424.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Men S, Chen H, Chen S, Zheng S, Shen X, Wang C, et al. Effects of supplemental nitrogen application on physiological characteristics, dry matter and nitrogen accumulation of winter rapeseed (Brassica napus L.) under waterlogging stress. Sci Rep. 2020;10:10201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Yan LY, Kang YP, Lei Y, Huang JQ, Wan LY, Liao BS. First report of Sclerotinia sclerotiorum causing sclerotinia blight on peanut (Arachis hypogaea) in Northeastern China. Plant Dis. 2014;98:156.

    Article  CAS  PubMed  Google Scholar 

  13. Zou X, Hu C, Zeng L, Cheng Y, Xu M, Zhang X. A comparison of screening methods to identify waterlogging tolerance in the field in Brassica napus L. during plant ontogeny. PLoS ONE. 2014;9: e89731.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Abiko T, Kotula L, Shiono K, Malik AI, Colmer TD, Nakazono M. Enhanced formation of aerenchyma and induction of a barrier to radial oxygen loss in adventitious roots of Zea nicaraguensis contribute to its waterlogging tolerance as compared with maize (Zea mays ssp. mays). Plant Cell Environ. 2012;35:1618–30.

    Article  CAS  PubMed  Google Scholar 

  15. Yamauchi T, Abe F, Tsutsumi N, Nakazono M. Root cortex provides a venue for gas-space formation and is essential for plant adaptation to waterlogging. Front Plant Sci. 2019;10:259.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Pedersen O, Nakayama Y, Yasue H, Kurokawa Y, Takahashi H, Heidi Floytrup A, et al. Lateral roots, in addition to adventitious roots, form a barrier to radial oxygen loss in Zea nicaraguensis and a chromosome segment introgression line in maize. New Phytol. 2021;229:94–105.

    Article  CAS  PubMed  Google Scholar 

  17. van Dongen JT, Fröhlich A, Ramírez-Aguilar SJ, Schauer N, Fernie AR, Erban A, et al. Transcript and metabolite profiling of the adaptive response to mild decreases in oxygen concentration in the roots of Arabidopsis plants. Ann Bot. 2009;103:269–80.

    Article  PubMed  Google Scholar 

  18. Xie YY, Yu CB, Hou JJ, Hu XJ, Li YS, Shen H, et al. Effect of low oxygen stress on adventitious roots of rape seedlings. Chin J Oil Crop Sci. 2013;35:284–9.

    Google Scholar 

  19. Dennis ES, Dolferus R, Ellis M, Rahman M, Wu Y, Hoeren FU, et al. Molecular strategies for improving waterlogging tolerance in plants. J Exp Bot. 2000;51:89–97.

    Article  CAS  PubMed  Google Scholar 

  20. Yang SH, Choi D. Characterization of genes encoding ABA 8′-hydroxylase in ethylene-induced stem growth of deepwater rice (Oryza sativa L.). Biochem Biophys Res Commun. 2006;350:685–90.

    Article  CAS  PubMed  Google Scholar 

  21. Bashar KK. Hormone dependent survival mechanisms of plants during post-waterlogging stress. Plant Signal Behav. 2018;13: e1529522.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Hu J, Ren B, Dong S, Liu P, Zhao B, Zhang J. Comparative proteomic analysis reveals that exogenous 6-benzyladenine (6-BA) improves the defense system activity of waterlogged summer maize. BMC Plant Biol. 2020;20:44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Qi X, Li Q, Ma X, Qian C, Wang H, Ren N, et al. Waterlogging-induced adventitious root formation in cucumber is regulated by ethylene and auxin through reactive oxygen species signalling. Plant Cell Environ. 2019;42:1458–70.

    Article  CAS  PubMed  Google Scholar 

  24. Qi X, Li Q, Shen J, Qian C, Xu X, Xu Q, et al. Sugar enhances waterlogging-induced adventitious root formation in cucumber by promoting auxin transport and signalling. Plant Cell Environ. 2020;43:1545–57.

    Article  CAS  PubMed  Google Scholar 

  25. Zhang XK, Fan QX, Chen J, Li JN, Wang HZ. Physiological reaction differences of different waterlogging tolerant genotype rapeseed (Brassica napus L.) to anoxia. Sci Agric Sin. 2007;40:485–91.

    CAS  Google Scholar 

  26. Saigo T, Wang T, Watanabe M, Tohge T. Diversity of anthocyanin and proanthocyanin biosynthesis in land plants. Curr Opin Plant Biol. 2020;55:93–9.

    Article  CAS  PubMed  Google Scholar 

  27. Ding LN, Li M, Guo XJ, Tang MQ, Cao J, Wang Z, et al. Arabidopsis GDSL1 overexpression enhances rapeseed Sclerotinia sclerotiorum resistance and the functional identification of its homolog in Brassica napus. Plant Biotechnol J. 2020;18:1255–70.

    Article  CAS  PubMed  Google Scholar 

  28. Qin L, Sun L, Wei L, Yuan J, Kong F, Zhang Y, Miao X, Xia G, Liu S. Maize SRO1e represses anthocyanin synthesis through regulating the MBW complex in response to abiotic stress. Plant J. 2021;105:1010–25.

    Article  CAS  PubMed  Google Scholar 

  29. Gou JY, Felippes FF, Liu CJ, Weigel D, Wang JW. Negative regulation of anthocyanin biosynthesis in Arabidopsis by a miR156-targeted SPL transcription factor. Plant Cell. 2011;23:1512–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Soubeyrand E, Basteau C, Hilbert G, van Leeuwen C, Delrot S, Gomès E. Nitrogen supply affects anthocyanin biosynthetic and regulatory genes in grapevine cv. Cabernet-Sauvignon berries. Phytochemistry. 2014;103:38–49.

    Article  CAS  PubMed  Google Scholar 

  31. Xie Y, Tan H, Ma Z, Huang J. DELLA proteins promote anthocyanin biosynthesis via sequestering MYBL2 and JAZ suppressors of the MYB/bHLH/WD40 complex in Arabidopsis thaliana. Mol Plant. 2016;9:711–21.

    Article  CAS  PubMed  Google Scholar 

  32. An JP, Zhang XW, Bi SQ, You CX, Wang XF, Hao YJ. The ERF transcription factor MdERF38 promotes drought stress-induced anthocyanin biosynthesis in apple. Plant J. 2020;101:573–89.

    Article  CAS  PubMed  Google Scholar 

  33. Ni J, Zhao Y, Tao R, Yin L, Gao L, Strid Å, et al. Ethylene mediates the branching of the jasmonate-induced flavonoid biosynthesis pathway by suppressing anthocyanin biosynthesis in red Chinese pear fruits. Plant Biotechnol J. 2020;18:1223–40.

    Article  CAS  PubMed  Google Scholar 

  34. Cheng J, Wei G, Zhou H, Gu C, Vimolmangkang S, Liao L, Han Y. Unraveling the mechanism underlying the glycosylation and methylation of anthocyanins in peach. Plant Physiol. 2014;166:1044–58.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Naing AH, Park DY, Park KI, Kim CK. Differential expression of anthocyanin structural genes and transcription factors determines coloration patterns in gerbera flowers. 3 Biotech. 2018;8:393.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Wang N, Liu W, Zhang T, Jiang S, Xu H, Wang Y, et al. Transcriptomic analysis of red-fleshed apples reveals the novel role of MdWRKY11 in flavonoid and anthocyanin biosynthesis. J Agric Food Chem. 2018;66:7076–86.

    Article  CAS  PubMed  Google Scholar 

  37. Tanaka Y, Sasaki N, Ohmiya A. Biosynthesis of plant pigments: anthocyanins, betalains and carotenoids. Plant J. 2008;54:733–49.

    Article  CAS  PubMed  Google Scholar 

  38. Walker AR, Davison PA, Bolognesi-Winfield AC, James CM, Srinivasan N, Blundell TL, et al. The TRANSPARENT TESTA GLABRA1 locus, which regulates trichome differentiation and anthocyanin biosynthesis in Arabidopsis, encodes a WD40 repeat protein. Plant Cell. 1999;11:1337–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Gonzalez A, Zhao M, Leavitt JM, Lloyd AM. Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings. Plant J. 2008;53:814–27.

    Article  CAS  PubMed  Google Scholar 

  40. Shin DH, Cho M, Choi MG, Das PK, Lee SK, Choi SB, et al. Identification of genes that may regulate the expression of the transcription factor production of anthocyanin pigment 1 (PAP1)/MYB75 involved in Arabidopsis anthocyanin biosynthesis. Plant Cell Rep. 2015;34:805–15.

    Article  CAS  PubMed  Google Scholar 

  41. Bedoya LC, Martínez F, Orzáez D, Daròs JA. Visual tracking of plant virus infection and movement using a reporter MYB transcription factor that activates anthocyanin biosynthesis. Plant Physiol. 2012;158:1130–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Wang L, Lu W, Ran L, Dou L, Yao S, Hu J, et al. R2R3-MYB transcription factor MYB6 promotes anthocyanin and proanthocyanidin biosynthesis but inhibits secondary cell wall formation in Populus tomentosa. Plant J. 2019;99:733–51.

    Article  CAS  PubMed  Google Scholar 

  43. Sun C, Deng L, Du M, Zhao J, Chen Q, Huang T, et al. A transcriptional network promotes anthocyanin biosynthesis in tomato flesh. Mol Plant. 2020;13:42–58.

    Article  CAS  PubMed  Google Scholar 

  44. Yan S, Chen N, Huang Z, Li D, Zhi J, Yu B, et al. Anthocyanin fruit encodes an R2R3-MYB transcription factor, SlAN2-like, activating the transcription of SlMYBATV to fine-tune anthocyanin content in tomato fruit. New Phytol. 2020;225:2048–63.

    Article  CAS  PubMed  Google Scholar 

  45. Bai H, Song Z, Zhang Y, Li Z, Wang Y, Liu X, et al. The bHLH transcription factor PPLS1 regulates the color of pulvinus and leaf sheath in foxtail millet (Setaria italica). Theor Appl Genet. 2020;133:1911–26.

    Article  CAS  PubMed  Google Scholar 

  46. An JP, Wang XF, Zhang XW, Xu HF, Bi SQ, You CX, et al. An apple MYB transcription factor regulates cold tolerance and anthocyanin accumulation and undergoes MIEL1-mediated degradation. Plant Biotechnol J. 2020;18:337–53.

    Article  CAS  PubMed  Google Scholar 

  47. Zhu Z, Li G, Liu L, Zhang Q, Han Z, Chen X, et al. A R2R3-MYB transcription factor, VvMYBC2L2, functions as a transcriptional repressor of anthocyanin biosynthesis in grapevine (Vitis vinifera L.). Molecules. 2018;24:92.

    Article  PubMed Central  Google Scholar 

  48. Liu R, Ding LN, Li M, Cao W, Wang YK, Wang WJ, et al. Characterization of a rapeseed anthocyanin-more mutant with enhanced resistance to Sclerotinia sclerotiorum. J Plant Growth Regul. 2019;39:703–16.

    Article  Google Scholar 

  49. Tsikas D. Assessment of lipid peroxidation by measuring malondialdehyde (MDA) and relatives in biological samples: analytical and biological challenges. Anal Biochem. 2017;524:13–30.

    Article  CAS  PubMed  Google Scholar 

  50. Feng CZ, Chen Y, Wang C, Kong YH, Wu WH, Chen YF. Arabidopsis RAV1 transcription factor, phosphorylated by SnRK2 kinases, regulates the expression of ABI3, ABI4, and ABI5 during seed germination and early seedling development. Plant J. 2014;80:654–68.

    Article  CAS  PubMed  Google Scholar 

  51. Prasad KVSK, Xing D, Reddy ASN. Vascular plant one-zinc-finger (VOZ) transcription factors are positive regulators of salt tolerance in Arabidopsis. Int J Mol Sci. 2018;19:3731.

    Article  PubMed Central  Google Scholar 

  52. Li B, Zheng JC, Wang TT, Min DH, Wei WL, Chen J, et al. Expression analyses of soybean VOZ transcription factors and the role of GmVOZ1G in drought and salt stress tolerance. Int J Mol Sci. 2020;21:2177.

    Article  CAS  PubMed Central  Google Scholar 

  53. Wang J, Wang R, Fang H, Zhang C, Zhang F, Hao Z, et al. Two VOZ transcription factors link an E3 ligase and an NLR immune receptor to modulate immunity in rice. Mol Plant. 2021;14:253–66.

    Article  PubMed  Google Scholar 

  54. Celesnik H, Ali GS, Robison FM, Reddy AS. Arabidopsis thaliana VOZ (vascular plant one-zinc finger) transcription factors are required for proper regulation of flowering time. Biol Open. 2013;2:424–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Kumar S, Choudhary P, Gupta M, Nath U. Vascular plant one-zinc finger 1 (VOZ1) and VOZ2 interact with CONSTANS and promote photoperiodic flowering transition. Plant Physiol. 2018;176:2917–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Luo D, Qu L, Zhong M, Li X, Wang H, Miao J, et al. Vascular plant one-zinc finger 1 (VOZ1) and VOZ2 negatively regulate phytochrome B-mediated seed germination in Arabidopsis. Biosci Biotechnol Biochem. 2020;84:1384–93.

    Article  CAS  PubMed  Google Scholar 

  57. Gao HN, Jiang H, Cui JY, You CX, Li YY. Review: The effects of hormones and environmental factors on anthocyanin biosynthesis in apple. Plant Sci. 2021;312: 111024.

    Article  CAS  PubMed  Google Scholar 

  58. Yang H, Wang T, Yu X, Yang Y, Wang C, Yang Q, et al. Enhanced sugar accumulation and regulated plant hormone signalling genes contribute to cold tolerance in hypoploid Saccharum spontaneum. BMC Genom. 2020;21:507.

    Article  CAS  Google Scholar 

  59. Zhang Q, Liu X, Zhang Z, Liu N, Li D, Hu L. Melatonin improved waterlogging tolerance in Alfalfa (Medicago sativa) by reprogramming polyamine and ethylene metabolism. Front Plant Sci. 2019;10:44.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Zeng R, Chen T, Wang X, Cao J, Li X, Xu X, et al. Physiological and expressional regulation on photosynthesis, starch and sucrose metabolism response to waterlogging stress in peanut. Front Plant Sci. 2021;12: 601771.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Zou X, Tan X, Hu C, Zeng L, Lu G, Fu G, et al. The transcriptome of Brassica napus L. roots under waterlogging at the seedling stage. Int J Mol Sci. 2013;14:2637–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Bai YC, Li CL, Zhang JW, Li SJ, Luo XP, Yao HP, et al. Characterization of two tartary buckwheat R2R3-MYB transcription factors and their regulation of proanthocyanidin biosynthesis. Physiol Plant. 2014;152:431–40.

    Article  CAS  PubMed  Google Scholar 

  63. Feng X, Zhang Y, Wang H, Tian Z, Xin S, Zhu P. The dihydroflavonol 4-reductase BoDFR1 drives anthocyanin accumulation in pink-leaved ornamental kale. Theor Appl Genet. 2021;134:159–69.

    Article  CAS  PubMed  Google Scholar 

  64. Matsui K, Oshima Y, Mitsuda N, Sakamoto S, Nishiba Y, Walker AR, et al. Buckwheat R2R3 MYB transcription factor FeMYBF1 regulates flavonol biosynthesis. Plant Sci. 2018;274:466–75.

    Article  CAS  PubMed  Google Scholar 

  65. Yu Z, Dong W, Teixeira da Silva JA, He C, Si C, Duan J. Ectopic expression of DoFLS1 from Dendrobium officinale enhances flavonol accumulation and abiotic stress tolerance in Arabidopsis thaliana. Protoplasma. 2021;258:803–15.

    Article  CAS  PubMed  Google Scholar 

  66. Kodama M, Brinch-Pedersen H, Sharma S, Holme IB, Joernsgaard B, Dzhanfezova T, et al. Identification of transcription factor genes involved in anthocyanin biosynthesis in carrot (Daucus carota L.) using RNA-Seq. BMC Genom. 2018;19:811.

    Article  CAS  Google Scholar 

  67. Zhao M, Li J, Zhu L, Chang P, Li L, Zhang L. Identification and characterization of MYB-bHLH-WD40 regulatory complex members controlling anthocyanidin biosynthesis in blueberry fruits development. Genes. 2019;10:496.

    Article  CAS  PubMed Central  Google Scholar 

  68. An JP, Zhang XW, Liu YJ, Wang XF, You CX, Hao YJ. ABI5 regulates ABA-induced anthocyanin biosynthesis by modulating the MYB1-bHLH3 complex in apple. J Exp Bot. 2021;72:1460–72.

    Article  CAS  PubMed  Google Scholar 

  69. Zhuang H, Lou Q, Liu H, Han H, Wang Q, Tang Z, et al. Differential regulation of anthocyanins in green and purple turnips revealed by combined de novo transcriptome and metabolome analysis. Int J Mol Sci. 2019;20:4387.

    Article  CAS  PubMed Central  Google Scholar 

  70. Chen HY, Hsieh EJ, Cheng MC, Chen CY, Hwang SY, Lin TP. ORA47 (octadecanoid-responsive AP2/ERF-domain transcription factor 47) regulates jasmonic acid and abscisic acid biosynthesis and signaling through binding to a novel cis-element. New Phytol. 2016;211:599–613.

    Article  CAS  PubMed  Google Scholar 

  71. Zhou X, Zhang ZL, Park J, Tyler L, Yusuke J, Qiu K, et al. The ERF11 transcription factor promotes internode elongation by activating gibberellin biosynthesis and signaling. Plant Physiol. 2016;171:2760–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Fan W, Yang Y, Wang Z, Yin Y, Yu C, Shi Q, Guo J, Xuan L, Hua J. Molecular cloning and expression analysis of three ThERFs involved in the response to waterlogging stress of Taxodium ‘Zhongshanshan406’, and subcellular localization of the gene products. PeerJ. 2018;6: e4434.

    Article  PubMed  PubMed Central  Google Scholar 

  73. Jiang W, Zhang X, Song X, Yang J, Pang Y. Genome-wide identification and characterization of APETALA2/ethylene-responsive element binding factor superfamily genes in soybean seed development. Front Plant Sci. 2020;11: 566647.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Qu Y, Nong Q, Jian S, Lu H, Zhang M, Xia K. An AP2/ERF gene, HuERF1, from pitaya (Hylocereus undatus) positively regulates salt tolerance. Int J Mol Sci. 2020;21:4586.

    Article  CAS  PubMed Central  Google Scholar 

  75. Sakata Y, Nakamura I, Taji T, Tanaka S, Quatrano RS. Regulation of the ABA-responsive Em promoter by ABI3 in the moss Physcomitrella patens: role of the ABA response element and the RY element. Plant Signal Behav. 2010;5:1061–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Amorim LLB, da Fonseca Dos Santos R, Neto JPB, Guida-Santos M, Crovella S, Benko-Iseppon AM. Transcription factors involved in plant resistance to pathogens. Curr Protein Pept Sci. 2017;18:335–51.

    Article  CAS  PubMed  Google Scholar 

  77. Baillo EH, Kimotho RN, Zhang Z, Xu P. Transcription factors associated with abiotic and biotic stress tolerance and their potential for crops improvement. Genes. 2019;10:771.

    Article  CAS  PubMed Central  Google Scholar 

  78. Wang Y, Liu C, Li K, Sun F, Hu H, Li X, et al. Arabidopsis EIN2 modulates stress response through abscisic acid response pathway. Plant Mol Biol. 2007;64:633–44.

    Article  CAS  PubMed  Google Scholar 

  79. Qi T, Song S, Ren Q, Wu D, Huang H, Chen Y, et al. The Jasmonate-ZIM-domain proteins interact with the WD-repeat/bHLH/MYB complexes to regulate jasmonate-mediated anthocyanin accumulation and trichome initiation in Arabidopsis thaliana. Plant Cell. 2011;23(5):1795–814.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Peng Z, Han C, Yuan L, Zhang K, Huang H, Ren C. Brassinosteroid enhances jasmonate-induced anthocyanin accumulation in Arabidopsis seedlings. J Integr Plant Biol. 2011;53:632–40.

    Article  CAS  PubMed  Google Scholar 

  81. Wang CK, Han PL, Zhao YW, Ji XL, Yu JQ, You CX, et al. Auxin regulates anthocyanin biosynthesis through the auxin repressor protein MdIAA26. Biochem Biophys Res Commun. 2020;533:717–22.

    Article  CAS  PubMed  Google Scholar 

  82. Zhang Y, Liu Z, Liu J, Lin S, Wang J, Lin W, et al. GA-DELLA pathway is involved in regulation of nitrogen deficiency-induced anthocyanin accumulation. Plant Cell Rep. 2017;36:557–69.

    Article  CAS  PubMed  Google Scholar 

  83. Li D, Guo Y, Zhang D, He S, Gong J, Ma H, et al. Melatonin represses oil and anthocyanin accumulation in seeds. Plant Physiol. 2020;183:898–914.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Oraei M, Panahirad S, Zaare-Nahandi F, Gohari G. Pre-veraison treatment of salicylic acid to enhance anthocyanin content of grape (Vitis vinifera L.) berries. J Sci Food Agric. 2019;99:5946–52.

    Article  CAS  PubMed  Google Scholar 

  85. Zou X, Jiang Y, Zheng Y, Zhang M, Zhang Z. Prolyl 4-hydroxylase genes are subjected to alternative splicing in roots of maize seedlings under waterlogging. Ann Bot. 2011;108:1323–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Cao H, Wang J, Dong X, Han Y, Ma Q, Ding Y, et al. Carotenoid accumulation affects redox status, starch metabolism, and flavonoid/anthocyanin accumulation in citrus. BMC Plant Biol. 2015;15:27.

    Article  PubMed  PubMed Central  Google Scholar 

  87. Chen T, Zhang B. Measurements of proline and malondialdehyde contents and antioxidant enzyme activities in leaves of drought stressed cotton. Bio-Protoc. 2016;6: e1913.

    Article  Google Scholar 

  88. Zhang Z, Huang R. Analysis of malondialdehyde, chlorophyll proline, soluble sugar, and glutathione content in Arabidopsis seedling. Bio-Protoc. 2013;3: e817.

    Article  Google Scholar 

  89. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  92. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011;12:323.

    Article  CAS  Google Scholar 

  93. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:129–40.

    Article  Google Scholar 

  94. Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11:R14.

    Article  PubMed  PubMed Central  Google Scholar 

  95. Mao X, Cai T, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification using the KEGG orthology (KO) as a controlled vocabulary. Bioinformatics. 2005;21:3787–93.

    Article  CAS  PubMed  Google Scholar 

  96. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25:402–8.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Ye Tao of Shanghai Lingen Biotechnology Co., Ltd for assistance with reference sequence alignment. We thank Shuo Xu of Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, for his help in raw RNA-seq data submission.

Funding

Financial support was provided by the National Natural Science Foundation of China (32172060 and 31471527), the Natural Science Foundation of Jiangsu Province (BK20201416), Jiangsu Agriculture science and technology innovation fund (CX(21)3112), and the National Key R & D Program of China (2016YFD0101904 and 2016YFD0100305).

Author information

Authors and Affiliations

Authors

Contributions

XLT conceived and designed research. LND, RL, TL and XYL conducted the experiments and interpreted the data. LND, JC, and XLT contributed new reagents or analytical tools. ML, WJW and YKY contributed to data analysis. LND wrote the manuscript. TL, JC and XLT revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xiao-Li Tan.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors consent to the publication of this manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

The phenotypes of seedlings after 8 day of recovery growth following waterlogging for 0, 12, 24 and 36 h. Bars = 1 cm. Figure S2. Sequencing randomness analysis of the waterlogging-treatment and control groups in both the WT and am mutant. Figure S3. Schematic diagram of the expression level saturation distribution curve. Figure S4. Schematic diagram of the correlations between samples. Figure S5. The MA and volcano plots of gene expression in (a) amCK-vs-amWL and (b) Z11CK-vs-Z11WL. “Z11CK-vs-Z11WL” indicates waterlogging-treated WT compared with the untreated WT, and “amCK-vs-amWL” indicates waterlogging-treated am mutant compared with the untreated am mutant. Figure S6. Co-expression clustering showing the expression profile of DEGs in WT (a) and am mutant (b). The X-axis represents with or without waterlogging treatment. The Y-axis represents the value of the relative expression level [log10 (FPKM + 1)].

Additional file 2: Table S1.

List of DEGs detected in the WT ZS11after the waterlogging treatment. Table S2. List of DEGs detected in the am mutant after the waterlogging treatment. Table S3. Differential expression of subclusters 1–5 in the am mutant. Table S4. Differential expression of subclusters 1–6 in the WT. Table S5. All the up- and downregulated DEGs with GO annotations in the am mutant. Table S6. All the up- and downregulated DEGs with GO annotations in the WT. Table S7. KEGG pathway enrichment of the DEGs in the WT. Table S8. KEGG pathway enrichment of the DEGs in the am mutant. Table S9. The DEG enrichment in the phenylalanine biosynthetic and metabolic pathways. Table S10. A list of differentially expressed TFs in the WT. Table S11. A list of differentially expressed TFs in the am mutant. Table S12. The DEG enrichment in plant hormone signal transduction pathways. Table S13. New transcript predictions as assessed by RNA-seq. Table S14. Primers used for qPCR validation of the selected DEGs.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ding, LN., Liu, R., Li, T. et al. Physiological and comparative transcriptome analyses reveal the mechanisms underlying waterlogging tolerance in a rapeseed anthocyanin-more mutant. Biotechnol Biofuels 15, 55 (2022). https://doi.org/10.1186/s13068-022-02155-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-022-02155-5

Keywords