Phenotypic analysis of different barley varieties under the waterlogging treatment
The phenotypes of two genotypes (TX9425 and Franklin) after 21-day waterlogging treatment are shown in Fig. 1. Franklin leaves became more wilted and chlorotic than TX9425 under waterlogging treatment. Plant height, tillers, leaf area, shoot fresh weight and dry weight of Franklin significantly decreased. Comparing with a small decline values were detected in TX9425 leaves (Table 1). No significant difference between Franklin and TX9425 under control condition was observed in adventitious root parameters. However, after three-week waterlogging treatment, the adventitious root length, surface area, volume and number of TX9425 significantly increased and the fold change value ranged from 2.36 to 4.06 compared to the control, while there was no significant difference in Franklin except adventitious root number (Table 2). Therefore, the performance of Franklin and TX9425 displayed significant difference after three-week waterlogging treatment.
Table 1
The effect of waterlogging on agronomic traits between TX9425 and Franklin
Treatment
|
Leaf Age
|
Leaf Chlorosis
|
Plant Height
|
Tillers
|
Leaf Area
(cm2)
|
Shoot Fresh
|
Shoot Dry
|
|
|
(cm)
|
|
Weight (g)
|
Weight (g)
|
Franklin
|
|
|
|
|
|
|
|
Control
|
7.84 ± 0.44a
|
0.36 ± 0.45a
|
47.35 ± 2.39a
|
9.14 ± 3.25a
|
33.31 ± 4.32a
|
24 ± 1.13a
|
2.32 ± 0.15a
|
Waterlogging
|
7.97 ± 0.51a
|
4.96 ± 0.89b
|
28.38 ± 2.34b
|
6.75 ± 1.95b
|
16.67 ± 3.14b
|
9.97 ± 0.52b
|
1.35 ± 0.12b
|
TX9425
|
|
|
|
|
|
|
|
Control
|
7.67 ± 0.4a
|
0.64 ± 0.49a
|
52.05 ± 3.48a
|
6.66 ± 1.54a
|
41.19 ± 6.15a
|
23.99 ± 1.4a
|
2.4 ± 0.2a
|
Waterlogging
|
7.86 ± 0.35a
|
1 ± 0.41b
|
45.11 ± 3.67b
|
6.41 ± 1.86a
|
34.15 ± 7.82b
|
20.24 ± 1.39b
|
2.3 ± 0.2a
|
The different letters within a column for the same line represent significant difference between waterlogging treatment and control. |
Table 2
The effect of waterlogging on the adventitious roots between TX9425 and Franklin
Treatment
|
Total Adventitious
Length
|
Total Adventitious Surface
Area
|
Adventitious Average
Diameter
|
Total Adventitious
Root Volume
|
Adventitious Root Number
|
(cm)
|
(cm2)
|
(mm)
|
(cm3)
|
Franklin
|
|
|
|
|
|
Control
|
63.42 ± 12.36 a
|
14.8 ± 5.72a
|
0.75 ± 0.21a
|
0.29 ± 0.16a
|
6.17 ± 1.33a
|
Waterlogging
|
68.37 ± 12.31a
|
14.29 ± 3.04a
|
0.66 ± 0.04a
|
0.24 ± 0.06a
|
14.83 ± 2.33b
|
TX9425
|
|
|
|
|
|
Control
|
68.41 ± 12.07a
|
13.93 ± 2.96a
|
0.62 ± 0.09a
|
0.22 ± 0.07a
|
5.83 ± 1.33a
|
Waterlogging
|
168.85 ± 13.87b
|
33.24 ± 8.38b
|
0.63 ± 0.03a
|
0.52 ± 0.13b
|
23.67 ± 3.83b
|
The different letters within a column for the same line represent significant difference between waterlogging treatment and control. |
Physiological And Anatomical Analysis Of Different Barley Varieties Under The Waterlogging Treatment
As shown in Fig. 2A-C, a significant genotype difference of the activities of SOD, CAT and POD in leaf was found. The antioxidant enzyme activities in both varieties all decreased under waterlogging, while the decrease in the tolerant TX9425 was lower. Under waterlogging treatment, the increase of MDA content of Franklin was about 2.1-fold, but that of TX9425 only increased by 1.3-fold (Fig. 2D). The root SOD activity of Franklin and TX9425 increased by 1.2- and 1.5-fold, respectively (Fig. 2E). The CAT activity of Franklin root increased by 1.6-fold, and that of TX9425 increased by 2.1-fold change (Fig. 2G). Moreover, POD enzyme activity of TX9425 increased by 1.4-fold change, and no significant difference was observed in the root of Franklin between waterlogging treatment and control (Fig. 2F). On the contrary, MDA content of Franklin increased by 2.1-fold change compared with the control in root, but the change in TX9425 was not significant (Fig. 2H).
Barley leaf anatomy is a typical monocotyledonous type consisting of epidermis, mesophyll and vascular tissue. Intercellular spaces existed among the mesophyll cells in the control. Under waterlogging, mesophyll cells of Franklin were severely damaged; on the contrary, the leaves of TX9425 developed more lysigenous aerenchyma under waterlogging compared with that of control (Fig. 3A). The adventitious root of barley was composed of epidermis, cortex and cylinder of vascular tissues. Cortex parenchyma cells of adventitious root formed larger number of lysigenous aerenchyma under waterlogging condition, compared with small intercellular space under control condition. Remarkably, the proportion of TX9425 aerenchyma was significantly higher than Franklin after three-week treatment (Fig. 3B). Under waterlogging, the adventitious root was formed in the section of shoot base in both lines, and more adventitious root primordiums were observed in TX9425 than Franklin. Otherwise, in the absence of waterlogging, few adventitious roots were found in both accessions (Fig. 3C, Table 2).
Analysis Of Barley Root Transcriptome Under Waterlogging Stress
In order to reveal molecular mechanisms of barley in response to waterlogging stress, roots were collected from TX9425 and Franklin after 0h, 24 h and 72 h waterlogging treatments. Each sample was subjected to three replicate treatments, and a total of 18 libraries were constructed. A high-throughput Illumina sequencing platform was used to sequence the transcriptome of barley. After removing adaptor sequences, low quality reads, and the reads with more than 10% ambiguous “N” bases, 2.87–7.58 GB data were obtained from each sample. The Q20 values of all transcriptomes were all above 96.42%, and the Q30 values at least 92.41%, indicating the high-quality sequencing data in the RNA-seq experiments (Table S1). Averagely, more than 63% of the valid reads were mapped into the reference barley genome. A principal component analysis (PCA) was conducted on the RNA-seq data set of 18 samples. The control and treatment samples of two genotypes could be clearly separated by the first principal component (PC1), which accounted for 98.53% of the total variation (Fig. 4A).
Identification Of Degs In Two Barley Varieties In Response To Waterlogging Stress
We further compared the DEGs in the two barley varieties subjected to waterlogging stress. We found a total of 3064 DEGs in TX9425 and 2297 DEGs in Franklin after 24 h waterlogging stress compared to the control, by using the parameters of log2FC ≥ 1 and q value ≤ 0.05. 1335 DEGs were up-regulated and 1729 DEGs were down-regulated in TX9425, while there were 967 up-regulated genes and 1330 down-regulated genes in Franklin (Fig. 4B). By comparing the transcriptome profiles of TX9425 and Franklin, it was observed that total of 2183 DEGs were uniquely expressed in TX9425 only, whereas 1416 DEGs were distinctively found in Franklin under 24 h waterlogging treatment. In addition, 881 DEGs were common between the two genotypes (Fig. 4C).
5693 DEGs and 8462 DEGs were identified under waterlogging treatment (72 h) vs. control in TX9425 and Franklin, respectively. 2012 DEGs were up-regulated and 3681 DEGs were down-regulated in TX9425, while there were 3314 up-regulated genes and 5148 down-regulated genes in Franklin. There were more DEGs after 72 h waterlogging stress than 24 h waterlogging stress. The number of DEGs were significantly different between TX9425 and Franklin (Fig. 4B). 1664 DEGs were uniquely expressed in TX9425 only, whereas a total of 4083 DEGs were distinctively found in Franklin under 72 h waterlogging treatment. Besides, 4029 DEGs were common between the two genotypes (Fig. 4C).
In addition, to verify the reliability of the RNA-seq data, 10 DEGs were randomly selected for the qRT-PCR analysis. Significantly positive correlations were observed between qRT-PCR and RNA sequencing data (r2 = 0.82). These results suggested that the RNA-seq data was credible (Fig. 4D).
Functional Annotation Of Waterlogging-responsive Degs
Gene ontology (GO) functional classification analysis was performed to categorize the functions of DEGs during waterlogging stress. The DEGs could be classified into three main ontologies, namely Molecular function, Biological process, and Cellular component (Table S2). As determined through a GO enrichment analysis of these DEGs, the DEGs in TX9425 and Franklin under 24 h waterlogging stress mostly functioned in biological process, metabolic process, transferase activity and catalytic activity (Fig. 5AB). After 72 h waterlogging, the DEGs of TX9425 mainly functioned in metabolic process, biological process, organic cyclic compound binding, heterocyclic compound binding and catalytic activity. However, the DEGs in Franklin mostly showed localization, oxidation-reduction process, protein binding and catalytic activity (Fig. 5CD, Table S2).
For KEGG pathway enrichment analysis, these DEGs were significantly (p < 0.01) enriched into 27 KEGG pathways (Table S3). Under 24 h waterlogging stress, the DEGs of TX9425 were mostly enriched into metabolic pathways and biosynthesis of secondary metabolites. However, the DEGs in Franklin primarily associated with biosynthesis of secondary metabolites and phenylpropanoid biosynthesis. Under 72 h waterlogging stress, the DEGs of TX9425 were mostly enriched into biosynthesis of secondary metabolites, MAPK signaling pathway, toll-like receptor signaling pathway. However, the DEGs in Franklin primarily associated with biosynthesis of secondary metabolites, biosynthesis of antibiotics and phenylpropanoid biosynthesis.
Analysis Of Degs Related To The Energy Metabolism
Energy deprivation is one of the major factors affecting waterlogged plant survival. This is due to mitochondrial respiration in plant roots is limited by hypoxia, and the acceleration of carbohydrate metabolism play more important role for plant survival. The KEGG enrichment analysis showed that many DEGs were involved in the starch and sucrose metabolism, glycolysis/fermentation pathway. As expected, we found that several DEGs, such as sucrose synthase, pyruvate kinase family protein, ATP-dependent 6-phosphofructokinase, alpha-amylase/subtilisin inhibitor, fructose-bisphosphate aldolase 2, were significantly accumulated in both TX9425 and Franklin.
In addition, some DEGs involved in glycolysis/fermentation pathway such as alanine aminotransferase, glyceraldehyde-3-phosphate dehydrogenase C2, alcohol dehydrogenase 1, L-lactate dehydrogenase A, pyruvate decarboxylase-2, were also significantly induced by waterlogging stress in two genotypes. In this study, we found that some genes had different expression levels in two varieties. For example, the pyruvate kinase family protein (HORVU2Hr1G040570) and fructose-bisphosphate aldolase 2 (HORVU3Hr1G088500) were induced at higher levels in TX9425 than Franklin after 24 or 72 h waterlogging treatment. The expression level of ATP-dependent 6-phosphofructokinase (HORVU5Hr1G019030), Alpha-amylase/trypsin inhibitor (HORVU7Hr1G035020), alcohol dehydrogenase 1 (HORVU1Hr1G082250, HORVU4Hr1G016810) were first increased and then decreased in TX9425, while it continuously increased in Franklin. Consequently, TX9425 had a greater energy state than Franklin under waterlogging stress (Table S4).
Hormone Regulation Of Barley Occurrence Under Waterlogging Stress
Hormones play an important role in plant response to environment stress. Here, we identified some DEGs related to hormones which are mainly involved in the biosynthesis of ethylene and auxin. Ethylene is biosynthesized by the activation of 1-aminocyclopropane-1-carboxylicacid synthase (ACS) and ACC oxidase (ACO). 2 ACS and 6 ACO were identified in TX9425 and Franklin. Two ACO genes (HORVU5Hr1G067490, HORVU5Hr1G067530) were significantly accumulated in both varieties, but the genes inductions in TX9425 were greater. 31 DEGs involved in auxin metabolism were identified in two genotypes, including 23 down-regulated and 8 up-regulated genes. After 72 h of waterlogging treatment, the expressions levels of HORVU1Hr1G025670, HORVU3Hr1G064590 and HORVU3Hr1G084840 in TX9425 were significantly higher than that in Franklin (Table S4).
Degs Involved In Ros Scavenging And Cell Wall Modifying Enzymes
Reactive oxygen species (ROS), which are produced when plants experiencing adversity stresses, can damage the normal functions in plant cells. To survive, plants have evolved multiple strategies such as activating antioxidant system remove excess ROS. 124 DEGs involved in ROS scavenging were found in our study, and most of them were down-regulated. These DEGs involves in the synthesis of glutathione S-transferase, peroxidase, catalase, L-ascorbate oxidase, most of which (82 genes, 66.12% of 124) were related to peroxidase. 8 genes related to glutathione S-transferase and 8 genes related to peroxidase were up-regulated in both genotypes, and the fold changes of these genes in TX9425 were significantly higher than that in Franklin (Table S4).
To adapt to waterlogging stress, plants also have evolved many mechanisms, such as the formation of adventitious roots and aerenchyma. The formation of aerenchyma was related to cell wall biosynthesis and loosening. As expected, we found 34 DEGs were involved in cell wall modifying enzymes, such as xyloglucan galactosyltransferase, pectinesterase, respiratory burst oxidase homologue. 8 DEGs were significantly up-regulated in both genotypes. Under waterlogging stress, the genes of HORVU2Hr1G101150 and HORVU4Hr1G081670 in TX9425 had significantly higher expression levels than that in Franklin (Table S4).
Characterization of HvADH4 genes in barley
A total of 44 ADH genes were identified in the barley genome based on the BLAST program. These genes were named HvADH1- HvADH44 according to their order of distribution on the chromosomes (Table S5). In the HvADH genes family, the length of coding sequences ranged from 99 (HvADH17) to 1524 (HvADH37). The size of corresponding amino acids varied between 32 and 507.The theoretical isoelectric point (PI) of these genes ranged from 4.51 to 9.66, and the molecular weight (Mw) varied from 3.47 to 48.11 kDa.
In this study, 17 ADH genes were found to have differential expressions between waterlogging treatment and control (Fig. 6). The highest differential expression was found for the HvADH4 in TX9425, and there was about a 50- fold difference between 24 h and control. We thus performed a standard method to isolate the HvADH4 from TX9425. Sequencing of HvADH4 showed that the full-length gene was 1158 bp in length and encoded 385 amino acids. Multiple amino acids alignment showed that HvADH4 protein shared two highly conserved ADH GroES-like (amino acid 36–156) and Zinc-binding dehydrogenase domains (amino acid 205–336) (Fig. 7A). Phylogenetic tree indicated that the HvADH4 has relatively high homology with protein from Triticum turgidum, and relatively distant sequence homology with the protein from Setaria italica (Fig. 7B).
Overexpression of HvADH4 enhanced waterlogging tolerance
To further verify the function of barley HvADH4 (HORVU1Hr1G082250), transgenic Arabidopsis plants overexpressing the HvADH4 gene from TX9425 were generated. Five week-old plants of the WT and three homozygous T3 transgenic lines were selected for waterlogging stress experiments. As shown in Fig. 8, no obvious changes in morphological and developmental phenotypes appeared between the WT and transgenic lines under normal conditions, while the transgenic lines grew better than WT plants after two weeks of waterlogging (Fig. 8A). Under waterlogging conditions, plant height was reduced by 49.1% in the WT, and 31.2, 36.1 and 40.3% in the transgenic lines (Fig. 8B). Compared to the control, the SPAD value was 61.6% lower in the WT, and 41.4, 51.4, 48.8% lower in the transgenic lines (Fig. 8C). The shoot fresh weights of the transgenic lines were 29.2, 37.2 and 36.5%, respectively, which were lower than in the control, while 65.8% smaller than that in the WT (Fig. 8D). The shoot dry weight decreased by 51.0% in the WT, and by 29.7, 13.3and 22.9% in the transgenic lines (Fig. 8E). In addition, the root lengths of the WT plants reduced greater than that of the transgenic lines during waterlogging stress (Fig. 8F). Furthermore, the average survival rate of the transgenic lines after waterlogging was 81.8%, but that of the WT was only 37.4% (Fig. 8G). Taken together, these data indicate that the overexpression of HvADH4 in Arabidopsis significantly enhances plant waterlogging tolerance.
Overexpression of HvADH4 increased the ROS scavenging capacity
To investigate the difference in physiological response to the waterlogging stress between the WT plants and the transgenic plants, the activities of antioxidant enzymes (SOD, CAT, and POD), ADH activity and MDA content were examined under normal and waterlogging conditions. Transgenic lines showed higher ADH activity than the WT plants even when they were under control conditions, and its remained significantly higher at subsequent times (Fig. 9D). There were no significant differences in the activities of antioxidant enzymes between transgenic lines and WT under normal growth conditions. After waterlogging, the major antioxidant enzymes activities increased markedly in both WT and transgenic plants, reaching peak levels at 6 days of treatment and then decreasing after 9 days of treatment. However, the fold changes were significantly greater in the transgenic lines than in the WT (Fig. 9ABC). MDA content is an important indicator to measure the level of lipid peroxidation. As shown in Fig. 9E, the MDA content in WT plants was significantly higher than that in transgenic lines, and this difference was more pronounced in the 6 d samples. Therefore, these results suggest that the overexpression of HvADH4 enhanced the scavenging ability of ROS in the plants and reduced the oxidative damage of plants under waterlogging stress.