eISSN: 2353-9461
ISSN: 0860-7796
BioTechnologia
Current issue Archive About the journal Editorial board Abstracting and indexing Subscription Contact Instructions for authors Publication charge Ethical standards and procedures
Editorial System
Submit your Manuscript
SCImago Journal & Country Rank
1/2021
vol. 102
 
Share:
Share:
RESEARCH PAPERS

Integrating expression data and genomic sequences to investigate the transcriptional regulation in barley in response to abiotic stress

Azar Delavari
1
,
Zahra Zinati
1
,
Sima Sazegari
2
,
Ahmad Tahmasebi
2

  1. Department of Agroecology, College of Agriculture and Natural Resources of Darab, Shiraz University, Shiraz, Iran
  2. Institute of Biotechnology, College of Agriculture, Shiraz University, Shiraz, Iran
BioTechnologia vol. 102 (1) C pp. 21–32 C 2021
Online publish date: 2021/03/31
Get citation
 
PlumX metrics:
 

Introduction

Abiotic stress severely inhibits crop development and productivity worldwide, and hence, there is a great interest in understanding the various mechanisms by which plants cope with abiotic stress. To overcome abiotic stresses, plants have evolved adaptive developmental and physiological strategies through the regulation of transcriptional networks (Lee et al., 2017). Transcription factors (TFs) are one of the critical components of these networks. TFs perform their regulatory functions by binding to their specific binding sites of 5–25 bp length (Rani et al., 2007; Zhang et al., 2019). Extensive progress has been made in determining regulatory genes involved in stress responses. This could enable to develop plants tolerant to abiotic stress (Galiba et al., 2009).

Cis-acting regulatory elements are critical molecular switches that play a role in the transcriptional regulation of different biological processes (e.g., abiotic stress responses, hormone responses, and development process). It is possible to obtain an accurate knowledge of the regulatory systems in the stress-responsive genes by analyzing the cis-acting elements (Yamaguchi-Shinozaki and Shinozaki, 2005).

Understanding the intricate mechanisms that govern gene expression regulation is critical and challenging (Das et al., 2019). Determining the factors that regulate genes and identifying their related binding sites is the most common approach to achieve it (Wang et al., 2018). Cohen and Leach (2019) implemented in silico analysis of publicly available rice transcriptome data to identify genes and pathways regulated by abiotic stress, biotic stress, and both stress types.

Barley, as one of the oldest cereal crops, is an important food source for livestock as well as a source of food and drink for humans (Newton et al., 2011). It possesses natural tolerance to the most important abiotic stresses such as drought and salinity, which makes it an excellent model in research on abiotic stress (Munns et al., 2006; Sallam et al., 2019). Despite the genetic complexity of abiotic stress tolerance, there exist valuable genes in the genome of barley with potential for useing in biotechnological improvement programs for abiotic stresses. A genome-wide association study (GWAS) of salinity tolerance performed by Mwando et al. (2020) led to the detection of 19 loci consisting of 52 significant salt tolerant-associated markers. In another study, GWAS was performed in spring barley, and 12 genes that regulated traits under drought conditions were recognized (Thabet et al., 2020). Therefore, genomic and transcriptomic data on environmental stress responses of barley can provide valuable information on single and combined abiotic stress tolerance.

Several miRNAs have been confirmed to play a role in stress responses in many plants (Chuck et al., 2009). For instance, next-generation sequencing (NGS) profiling and northern blot analysis confirmed the upregulation and involvement of miR5655 and miR2933b in Arabidopsis under drought stress (Barciszewska-Pacak et al., 2015). The downregulation of mir419 and the upregulation of its target gene confirmed the regulatory role of this miRNA gene in barley (Deng et al., 2015). Transgenic rice lines in which OsmiR156 was overexpressed showed improved cell viability and growth under cold stress (Zhou and Tang, 2019). Moreover, Zare et al. (2019) reported the upregulation of miR414 and miR2102 in two barley genotypes under drought. However, no association was found between these two miRNAs and the two highest upregulated differential expressed genes HvsnLTP (non-specific lipid transport protein) and HvPiP1;4 (N-butyl-N-methylpiperidinium). Accordingly, the investigation of abiotic stress-associated miRNA families in barley provides substantial opportunities to enhance the tolerance of barley to abiotic stresses.

Most of the studies conducted thus far focus only on one type of stress. However, there is a network of pathways regulated by abiotic stresses (Zhu et al., 2013). Therefore, meta-analysis is a useful tool to identify differentially expressed genes (DEGs). An integrative meta-analysis of stress response studies using micro-array gene expression data and a system biology analysis would provide a potential approach for probing the regulatory gene sets. Gene sets such as TFs, transcription regulators (TRs), and protein kinases (PKs) govern the stress responses.

The present study was designed to analyze gene promoters to reveal conserved and enriched motifs with their possible roles in stress response. We used DEGs obtained from the meta-analysis of transcriptomic responses under multiple abiotic stresses in barley. Furthermore, miRNA identification was performed on the DEGs under abiotic stresses. Additionally, regulatory gene sets, TFs, TRs, and PKs families that orchestrate abiotic stress responses were identified through several bioinformatics tools such as PlantTFDB (Tian et al., 2019; Jin et al., 2017) and iTAK (Zheng et al., 2016). Finally, signaling networks were constructed using the TFs, TRs, and PKs, and the key regulators were detected.

Methods

Identification of DEGs

Microarray datasets for abiotic stresses such as drought, cold, salinity, and low-temperature stresses were extracted from Gene Expression Omnibus (www.ncbi.nlm.nih.gov/geo) – Table 1. All microarray gene expression data were implemented on a single platform (Barley Genome Array, GPL1340).

Table 1

Characteristics of data sets selected for meta-analysis

DatasetSamplesTreatmentPlatform
stresscontrol
GSE699393low temperature stressAffymetrix Barley Genome Array
GSE1033233cold stressAffymetrix Barley Genome Array
GSE15970269drought stressAffymetrix Barley Genome Array
GSE176691212drought stressAffymetrix Barley Genome Array
GSE2389666heat stressAffymetrix Barley Genome Array
GSE2782193cold stressAffymetrix Barley Genome Array
GSE2782233cold stressAffymetrix Barley Genome Array
GSE309733salinity stressAffymetrix Barley Genome Array
GSE317066drought stressAffymetrix Barley Genome Array
GSE560562salinity stressAffymetrix Barley Genome Array
GSE564373636drought stressAffymetrix Barley Genome Array
GSE6736766drought stressAffymetrix Barley Genome Array
GSE6990123drought stressAffymetrix Barley Genome Array
GSE63251212salinity stressAffymetrix Barley Genome Array

The probe-gene maps and probe annotation files were downloaded from the Affymetrix site (www.affymetrix.com). The robust multiple average (RMA) method was used to normalize the raw expression data of each array (Irizarry et al. 2003). For this purpose, the Affymetrix Expression Console software was used. When normalization was performed to adjust batch effects by applying the SVA R package (Leek et al., 2012), the empirical Bayes method was used (Johnson et al., 2007).

To find the probe sets that were expressed differentially based on the false discovery rate, a meta-analysis was performed using a nonparametric rank production method (Breitling et al., 2004). A probe set was considered to be differentially expressed if FDR was < 0.01. More analyses were performed on the 3721 founded probe sets. The analyses included determination of TFs and PKs, cis-element analysis, and construction of signaling networks. Overall, 1694 and 2027 probe sets showed significant upregulation and downregulation of expression, respectively. BLASTX from IPK Barley BLAST Server was used to identify the Gene ID accessions related to the probe sets.

Promoter motif analysis of genes corresponding to DEGs

The 1 kbp upstream anking regions of genes corresponding to DEGs were extracted from Ensembl Plants (http://plants.ensembl.org). MEME (Multiple EM for Motif Elicitation) is one of the most widely used tools for identifying novel signals in sets of biological sequences. It includes the discovery of new transcription factor binding sites (Bailey et al., 2006). MEME (meme. nbcr.net/meme/intro.html) was used to discover conserved motifs on the sequences with its default parameters, except for the maximum number of motifs, which was set to 20, with the threshold E-value of < 1e−4 (Bailey et al., 2009). The Tomtom v 5.0.1 tool (http://meme-suite.org/tools/tomtom) was used to eliminate redundant motifs and to de ne known CRE based on the JASPAR CORE 2018 database with the threshold E-value of 0.05 (Gupta et al., 2007; Khan et al., 2017). The GOMo tool (http://meme-suite.org/tools/gomo) was also applied to identify possible biological roles for motifs (Buske et al., 2010).

MAST (Motif Alignment and Search Tool) was used for searching biological sequence databases for sequences containing an occurrence of each motif in a given set of motifs. AME (Analysis of Motif Enrichment) of the MEME suite was used to identify known motifs that are relatively enriched in a set of sequences. To scan given sequences for individual matches to each of the motifs, FIMO (Find Individual Motif Occurrences) of the MEME suite was implemented (McLeay and Bailey, 2010).

Identification of TFs, TRs, and PKs

To identify TFs families and TRs, BLASTX search against the barley transcription factors (http://planttfdb.cbi.pku.edu) with a cutoff value of E ≤ 10−5 was performed on DEG sequences. Moreover, PKs in DEGs were identified from the BLASTX search against the iTAK database (http://bioinfo.bti.cornell.edu/cgi-bin/itak/index.cgi) with a cutoff value of E ≤ 10−5.

Construction of the signaling networks consisting of TFs, TRs, and PKs

The Pathway Studio software (Ariadne Genomics, 2010) was used to construct a network that is established based on the ResNet database in which direct connection is considered between TFs, TRs, and PKs and their neighbors. The downstream term was selected as the directionality of the relationship to detect upstream regulators. The software calculates the local connectivity (node degrees in the network built) and compares them with the connectivity (node degrees in the ResNet database).

Identification of candidate abiotic stress-associated miRNAs

To identify potential miRNAs that might be related to the downregulated DEGs, their sequences were used as queries for the psRNATarget server (http://plantgrn.noble.org/psRNATarget/) with default parameters, except for the maximum expectation, which was set to 3, and compared with the 71 published miRNAs related to Hordium vulgare. Figure 1 shows the schematic procedure of the main steps of this study.

Fig. 1

Schematic of the procedure for clarifying the aspects of response of barley to abiotic stress

/f/fulltexts/BTA/43334/BTA-102-1-43334-g001_min.jpg

Results

Promoter motif analysis of DEGs

The 1 kbp upstream flanking regions of genes corresponding to DEGs were analyzed to discover conserved motifs and consensus cis-regulatory elements (CREs).

In this regard, 1424 promoter regions corresponding to the upregulated probe sets and 2018 promoter regions corresponding to the downregulated probe sets were retrieved from Ensemble Plants. Twenty significant motifs were detected with lengths ranging from 15 to 50 bp in the promoters of DEGs by using MEME (Suppl. mater. 1). The GOMo analysis for the motifs demonstrated several biological functions. For the upregulated genes, GO motifs were related to kinase activity, monooxygenase activity, plasma membrane, chloroplast envelope, oxygen binding, and protein glycosylation. On the other hand, the downregulated motifs showed monooxygenase activity, tyrosine kinase signaling, nitrogen compound metabolic process, and protein heterodimerization activity.

As observed in motif analysis, the majority of them show a relationship with the APETALA2/ETHYLENE-RESPONSIVE FACTOR (AP2/ERF) domain and C2H2 zinc finger TF family (28% for each). Analysis of motif enrichment revealed that the highest number of TF members belong to the AP2/ERF class: more than 20% of the whole TFs (Table 2). We found 28 and 54 motifs that exist only in the upregulated and downregulated genes, respectively. Table 3 illustrates their functions in abiotic stresses in barley searched through the PLANT-CARE database.

Table 2

The classes of TFs that bind to the related motif

Number of TF members
upregulateddownregulated
AP2/ERF8181
bHLH6973
bZIP3027
C2H26362
Helix-turn-helix1213
MADS box factors1318
Tryptophan cluster factors812
Other classes that contain < 10 members6180
Total TF members337366
Table 3

Functions related to the detected motifs that were obtained by checking in the PLANTCARE database

Motif function
motif_alt_IDconsensusmotif classsite namefunction
EmBP-1ASACGTGGbZIPChs-unit 1part of a light-responsive element
CMTA3CCGCGTNNNCG-1 domainGTGGC-motifpart of a light-responsive element
EmBP-1ASACGTGGbZIPABREcis-acting element involved in the abscisic acid responsiveness
ARF5RCCGACAAB3 domainChs-unit 1part of a light-responsive element
RSC3CGCGCVNC6 zinc cluster factorsABREcis-acting element involved in abscisic acid responsiveness
GATA8NMGATCYRNC4 zinc finger-type factorsBox Ipart of a light-responsive element
Chs-CMA2bpart of a light-responsive element
SPL4NGTACGDNNSBP-type zinc fingerABREcis-acting element involved in abscisic acid responsiveness
Chs-CMA2bpart of a light-responsive element
SPL7CGTACGVCSBP-type zinc fingerABREcis-acting element involved in abscisic acid responsiveness
BHLH104GGCACGTGCCbHLHABREcis-acting element involved in abscisic acid responsiveness

Identification of the TFs, TRs, and PKs

DEG sequences were aligned against the PlntTFDB database to identify TFs. We found 35 TFs related to the upregulated genes and 44 TFs related to the downregulated genes, which, respectively, belong to 19 and 23 families of TFs. Four upregulated genes were categorized in the NAC and bZIP class, while in the bHLH family, seven members were downregulated, which was the highest number of downregulated TFs (Table 4). The C2H2 and NF-YA families had just three upregulated members. The MYB-related family had five out of seven and the WRKY family had five out of six downregulated members, which indicates that these TFs families are mostly downregulated in barley under abiotic stresses.

Table 4

TF families identified using BLASTX search against the barley TFs database

TF familyUpregulatedDownregulated
AP2/ERF-ERF32
C2H230
MYB-related25
Tify11
AP2/ERF-AP211
C3H11
NAC41
TUB23
bHLH17
GARP-G2-like22
NF-YA30
VOZ10
bZIP41
HSF10
NF-YB10
WRKY15
C2C2-Dof11
MYB20
PLATZ10
C2C2-GATA01
MADS-M-type01
C2C2-CO-like01
GRAS02
B301
HB-HD-ZIP02
zf-HD01
HB-other01
NF-YA02
FAR101
LIM01
Total3544

TRs family was one of the identified regulatory genes in our study, with 11 upregulated and 24 downregulated members. TRs related to the auxin/indole-3-acetic acid (AUX/IAA) family contains six members, out of which four members were downregulated. The GNAT (GCN5 related N terminal acetyltransferase) family had the highest number of downregulated TRs (Table 5).

Table 5

TR families identified using BLASTX search against the barley TFs database

TR familyUpregulatedDownregulated
AUX/IAA24
SWI/SNF-BAF60b11
HMG11
MBF110
Others58
Pseudo ARR-B10
GNAT04
SET01
LUG01
TAZ02
mTERF02
Total1124

Among the DEGs, 18 upregulated and 32 downregulated PK genes were categorized into six PK families. The largest upregulated PK family was CAMK with seven members. RLK-Pelle containing 18 members was the largest among the downregulated PKs (Table 6).

Table 6

PKs in DEGs were obtained from BLASTX search against the iTAK database

Kinase familyUpregulatedDownregulated
AGC01
CAMK74
CK110
CMGC33
RLK-Pelle618
TKL16
Total1832

Construction of the signaling networks consisting of the TFs, TRs, and PKs

The connection between TFs, TRs, and PKs and their neighbors was analyzed by network construction using the Pathway Studio software. Table 7 presents the top five genes that contain the highest number of relationships. AP2, a member of the AP2/ERF class of TFs, has 30, the highest number of relationships with other TFs and PKs genes (Table 7). Therefore, this gene can be introduced as one of the most critical TFs in barley involved in abiotic stresses. AP2 has 22 relationship nodes with TFs, one with PK, and one with TR gene. The role of AP2 has been reported previously. To help activate abscisic acid (ABA) and ethylene (ET) dependent and independent stress-responsive genes, many AP2/ERFs respond to ABA and ET (Xie et al., 2019; Hoang et al., 2017). The results showed that HY5 (bZIP TF family member), SHY2 (AUX/IAA TR family member), MYC2 (bHLH TF family member), and LHY (Helix-Turn-Helix TF family member) formed 21–24 nodes with other genes. These genes can be presented as the most important genes of the network. These genes have 6, 8, 6, and 5 nodes with TFs, respectively; 1, 4, 0, and 3 nodes with TRs, respectively; and no relationship with any of the PKs. The remaining genes in the network analysis had less than 17 nodes. The excel format and network image of the genes comprising all relationships and entities are provided as Supplementary materials 2.

Table 7

The calculated connectivity of TFs, TRs, and Pks with other genes obtained from the ResNet database

NameTypeConnectivityLocal connectivity
AP2 ( AT4G36920)protein11930
HY5 ( AT5G11260)protein10424
SHY2 (AT1G04240)protein6422
MYC2 (AT1G32640)protein7722
LHY ( AT1G01060)protein5721

Identification of candidate abiotic stress-associated miRNAs

One of the objectives of this study was to identify abiotic stress-associated miRNAs in barley. Hence, the potential miRNAs targeting the downregulated DEGs were predicted using the psRNATarget server. A total of 49 miRNAs belonging to 23 distinct conserved families were found (Fig. 2, Suppl. mater. 3). According to the output psRNATarget server, only five members of the candidate miRNAs participate in translation inhibition, and the remaining are involved in mRNA cleavage processes. The hvu-miR5049 was the largest family with 15 members, followed by the hvu-miR6214 family containing 5 members.

Fig. 2

Prediction of potential miRNAs targeting the downregulated DEGs using psRNATarget. A total of 49 miRNAs were found to belong to 23 distinct conserved families

/f/fulltexts/BTA/43334/BTA-102-1-43334-g002_min.jpg

Discussion

Transcriptional regulation plays the main role in the activation and suppression of gene expression. It is controlled mainly by gene promoters and their contributing cis-acting elements and motifs (Zou et al., 2011). According to the results obtained from FIMO, Motif3 was found in the promoters of the highest number (568) of upregulated genes (Suppl. mater. 1). As reported in the motif analysis part, the ERF family transcription factor binds to the Motif3 sequence. This shows that the ERF family transcription factor is the upstream component in a transcriptional cascade that functions as a mediator of the upregulation of gene expression caused by abiotic stress. Previously, it has been indicated that ERF TFs are involved in Arabidopsis (Dubois et al., 2013), rice (Schmidt et al., 2013), and wheat (Zhu et al., 2014) responses to drought, salt, and freezing stresses, respectively. Transgenic Arabidopsis expressing barley ERF TF showed also improved salt tolerance (Jung et al., 2007). Our analysis revealed the participation of the ERF TF family in abiotic stress response in barley.

TFs, TRs, and PKs among the DEGs of barley were identified in this study. It was previously shown that transcription factors such as dehydration-responsive element-binding factors (HvDREB1) and WRKYs (HvWRKY38) are involved in drought stress responses through regulation of the expression of stress-related genes (Gürel et al., 2016). NAC proteins are one of the largest families of TFs that play a role in plant growth and abiotic stress responses. Many studies have shown that overexpression of NAC TFs caused an improvement in stress tolerance such as salt tolerance in Arabidopsis (Jiang and Deyholos, 2006); drought, salinity, and low temperature tolerance in rice (Fang et al., 2008); and dehydration stress tolerance in soybean (Le et al., 2011). The NAC TF family with the highest members among the upregulated genes was also identified in our study. In addition to NAC, the Leucine Zipper (bZIP) family was also identified with the maximum number of members. It was reported that bZIP is one of the largest TFs families that is induced in response to different plant development stages and stress conditions (Pourabed et al., 2015). Both NAC and bZIP transcription factors are known as genes that are induced by the ABA hormone under abiotic stresses (Yoon et al., 2020). Overall, 114 bZIP proteins were reported in H. vulgare, some of which were upregulated under drought and cold conditions (Pourabed et al., 2015). bHLH TFs are the largest family of downregulated genes identified in our study. According to previous reports, nitrogen deficiency caused a decrease in the number of bHLH transcripts in barley (Comadira et al., 2015).

To date, many studies have indicated the role of WRKY TFs in response to cold and drought in barley (Mare et al., 2004), drought in common bean (Wu et al., 2017), and in other crops (reviewed by Bai et al., 2018). The WRKY TF family has 94 members in barley plants (H. vulgare ) (Liu et al., 2014). We identified six differentially expressed WRKY TF family members among barley DEGs, and most of them were downregulated. Although the MYB TF family is large and involved in abiotic stress responses in Arabidopsis (Lippold et al., 2009), rice (El-Kereamy et al., 2012), apple (Cao et al., 2013), and Triticum aestivum (Mao et al., 2011), there is a lack of information regarding the MYB TF family participation in abiotic stress in barley. Our studies revealed that five out of seven members of this family are downregulated.

AUX/IAA plays a role in the growth and development of plants such as wheat (Qiao et al., 2017), Arabidopsis (Kim et al., 2020), and barley (Shi et al., 2020). Our studies indicated that most of the AUX/IAA TR family members were downregulated during abiotic stresses. Nevertheless, the downregulation of ten members of the AUX/IAA family was also reported at the embryogenesis stage in barley following mannitol application (Muñoz Amatriaín et al., 2006).

According to our analysis, the RLK-Pelle PK family members were among the DEGs in barley. To the best of our knowledge, there is only one report on the role of receptor-like kinases (RLKs) in crops under abiotic stress. Yang et al. (2017a) reported RLKs as downregulated kinase genes with negative regulatory effects during drought stress in jute. CaMK is the largest upregulated PK family identified in our study. Several studies have been conducted on the role of calcium-dependent protein kinase (CDPK) in response to abiotic stresses (Yang et al., 2017b). However, there is no report on the role of other CAMK members in the response of barley to abiotic stresses.

The analysis of network connections of the TFs, TRs, and PKs performed in this study showed that HY5 (bZIP TF family member) and LHY (Helix-Turn-Helix TF family member) have a relationship with PHYB in the network, and PHYB has a relationship with short hypocotyl/suppressor of HY2 (SHY2) (AUX/IAA TR family). SHY2 encodes IAA3 and plays a role as a negative regulator of auxin signaling. It connects phytohormone signaling and stress signaling in root development (Li et al., 2020). It shows that HY5 and LHY are upstream genes of SHY2. PHYB belongs to the phytochrome family and plays a role in seed germination and de-etiolation, and it inhibits shade avoidance responses under a high ratio of red : far-red light (R : FR) in Arabidopsis (Reed et al., 1994; Legris et al., 2019).

MicroRNAs (miRNAs), a class of small noncoding RNAs, regulate the expression of genes involved in many different biological processes such as development, environmental adaptation, and stress tolerance by post-transcriptional regulation of gene expression (Alptekin et al., 2017; Parmar et al., 2020; Fan et al., 2020). Manipulation of abiotic stress-associated miRNAs and their targets may pave the way for improving crop performance under several abiotic stresses (Alptekin et al., 2017; Jiang et al., 2019). In the present study, 49 miRNAs belonging to 23 miRNA families were predicted in barley under abiotic stresses. The highest number of miRNAs belong to the conserved family hvu-miR5049 (Fig. 2). Previous studies have reported the role of miR5049, miR169, miR1120, and miR444 in abiotic stresses. For instance, hvu-miR5049 was significantly upregulated under drought conditions in H. vulgare (Hackenberg et al., 2015). miR5049b has been further reported to be upregulated in leaf and root tissues of a drought-tolerant wheat cultivar under drought stress (Akdogan et al., 2015). hvu-miR5049a was upregulated in leaves but downregulated in roots of barley under drought conditions (Hackenberg et al., 2015). The hvumiR169 family, with two members, was one of the abiotic stress-associated miRNAs identified in this study. The miR169 family is found in many different plant species and is considered to be the largest plant miRNA family. Nevertheless, its ubiquitous expression in various plant genera implies its important regulatory effects (Rao et al., 2020). miR169 was shown to be involved in several stresses. For instance, nuclear transcription factor Y (NF-Y), one of the targets of miR169 in several plants such as maize and Arabidopsis, was upregulated in both heat- and cold-stressed wheat (Ni et al., 2013; Gupta et al., 2014; Kumar et al., 2014; Sorin et al., 2014; Luan et al., 2015). It has been shown that the downregulation of miR169 conduces to the transcript accumulation of NFYA5, a critical positive regulator of drought stress resistance in Arabidopsis (Li et al., 2008). Similarly, the decreased expression of miR169 in Arabidopsis was observed by high-throughput sequencing and qRT-PCR analysis under both drought and salt stresses (Pegler et al., 2019). Furthermore, the regulatory role of miR169 and its targets in drought, salt, and cold stress responses in tomato have been recently reported (Rao et al., 2020).

The hvu-miR444 family with three members is another abiotic stress-associated miRNA identified in this study. hvu-miR444b was upregulated in leaves but downregulated in roots under drought conditions in barley (Hackenberg et al., 2015). The same miRNA was also found to be downregulated in response to salinity stress in barley (Deng et al., 2015). hvu-miR444 targets putative MADS-box transcription factor 27 and cryptochrome 1b (Kantar et al., 2010); thus, it controls gene expression related to plant development, morphology, and flowering time.

The miR1120 is yet another abiotic-responsive miRNA identified in this study among the downregulated DEGs in barley under abiotic stress conditions. In accordance with this result, Sinha et al. (2015) reported miR1120 downregulation under N-deficient conditions in wheat (Sinha et al., 2015).

In the present study, the 19 miRNA families involved in barley stress responses have not been reported previously. The members of these miRNA families can be investigated in future studies for revealing their roles in abiotic stresses.

Conclusions

Exploring the abiotic stress-responsive regulatory gene sets in barley in order to understand the mechanism of gene regulatory system under abiotic stress conditions is highly important. For this purpose, TFs, TRs, and PKs families that respond to abiotic stresses in barley were identified based on the analysis of DEGs by bioinformatics tools. Our results also identified new motifs involved in abiotic stress responses. This study sheds light on cis-elements that play a role in response to environmental stimuli in barley. Finally, critical genes and miRNAs were identified as potential candidate genes involved in responses to abiotic stresses.

References

1 

Akdogan G., Tufekci E.D., Uranbey S., Unvert T. (2016) miRNA-based drought regulation in wheat. Funct. Integr. Genomics 16: 221.

2 

Alptekin B., Langridge P., Budak H. (2017) Abiotic stress miRNomes in the Triticeae. Funct. Integr. Genomics 17: 145–170.

3 

Bai Y., Sunarti S., Kissoudis C., Visser R.G.F., van der Linden C.G. (2018) The role of tomato WRKY genes in plant responses to combined abiotic and biotic stresses. Front. Plant Sci. 9: 801.

4 

Bailey T.L., Boden M., Buske F.A., Frith M., Grant C.E., Clementi L., Ren J., Li W.W., Noble S. (2009) MEME Suite: tools for motif discovery and searching. Nucl. Acids Res. 37: W202–W208.

5 

Bailey T.L., Williams N., Misleh C., Li W.W. (2006) MEME: discovering and analyzing DNA and protein sequence motifs. Nucl. Acids Res. 34: W369–W373.

6 

Barciszewska-Pacak M., Milanowska K., Knop K., Bielewicz D., Nuc P., Plewka P., Pacak A.M., Vazquez F., Karlowski W., Jarmolowski A., Szweykowska-Kulinska Z. (2015) Arabidopsis microRNA expression regulation in a wide range of abiotic stress responses. Front. Plant Sci. 6: 410.

7 

Breitling R., Armengaud P., Amtmann A., Herzyk P. (2004) Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett. 573: 83–92.

8 

Buske F.A., Bodén M., Bauer D.C., Bailey T.L. (2010) Assigning roles to DNA regulatory motifs using comparative genomics. Bioinformatics 26: 860–866.

9 

Cao Z.H., Zhang S.Z., Wang R.K., Zhang R.F., Hao Y.J. (2013) Genome wide analysis of the apple MYB transcription factor family allows the identification of MdoMYB121 gene conferring abiotic stress tolerance in plants. PLoS One 8: e69955.

10 

Chuck G., Candela H., Hake S. (2009) Big impacts by small RNAs in plant development. Curr. Opin. Plant Biol. 12: 81–86.

11 

Cohen S.P., Leach J.E. (2019) Abiotic and biotic stresses induce a core transcriptome response in rice. Sci. Rep. 9: 6273.

12 

Comadira G., Rasool B., Karpinska B., Karpinska B., Morris J., Verrall S.R., Hedley P.E., Foyer C.H., Hancock R.D. (2015) Nitrogen deficiency in barley (Hordeum vulgare) seedlings induces molecular and metabolic adjustments that trigger aphid resistance. J. Exp. Bot. 66: 3639–3655.

13 

Das A., Pramanik K., Sharma R., Gantait S., Banerjee J. (2019) In-silico study of biotic and abiotic stress-related transcription factor binding sites in the promoter regions of rice germin-like protein genes. PloS one 14: e0211887.

14 

Deng P., Wang L., Cui L., Feng K., Liu F., Du X., Tong W., Nie X., Ji W., Weining S. (2015) Global identification of micro-RNAs and their targets in barley under salinity stress. PLoS One 10: e0137990.

15 

Dubois M., Skirycz A., Claeys H., Maleux K., Dhondt S., De Bodt S., Vanden Bossche R., De Milde L., Yoshizumi T., Matsui M., Inzé D. (2013) Ethylene response factor6 acts as a central regulator of leaf growth under water-limiting conditions in Arabidopsis. Plant Physiol. 162: 319–332.

16 

El-kereamy A., Bi Y.M., Ranathunge K., Beatty P.H., Good A.G., Rothestein S.J. (2012) The rice R2R3-MYB transcription factor OsMYB55 is involved in the tolerance to high temperature and modulates amino acid metabolism. PLoS One 7: e52030.

17 

Fan G., Liu Y., Du H., Kuang T., Zhang Y. (2020) Identification of drought-responsive miRNAs in Hippophae tibetana using high-throughput sequencing. 3 Biotech. 10(2): 53.

18 

Fang Y., You J., Xie K., Xie W., Xiong L. (2008) Systematic sequence analysis and identification of tissue-specific or stress-responsive genes of NAC transcription factor family in rice. Mol. Genet. Genomics 280: 547–563.

19 

Galiba G., Vágújfalvi A., Li C., Soltész A., Dubcovsky J. (2009) Regulatory genes involved in the determination of frost tolerance in temperate cereals. Plant Sci. 176: 12–19.

20 

Gupta S., Stamatoyannopoulos J.A., Bailey T.L., Noble W.S. (2007) Quantifying similarity between motifs. Genome Biol. 8: R24.

21 

Gupta O.P., Meena N.L., Sharma I., Sharma P. (2014) Differential regulation of microRNAs in response to osmotic, salt and cold stresses in wheat. Mol. Biol. Rep. 41: 4623–4629.

22 

Gürel F., Öztürk Z.N., Uçarl2, C., Rosellini D. (2016) Barley genes as tools to confer abiotic stress tolerance in crops. Front. Plant Sci. 7: 1137.

23 

Hackenberg M., Gustafson P.J., Langridge P., Shi B. (2015) Differential expression of microRNAs and other small RNAs in barley between water and drought conditions. Plant Biotechnol. J. 13: 2–13.

24 

Hoang X.L.T., Nhi D.N.H., Thu N.B.A., Thao N.P., Tran L.P. (2017) Transcription factors and their roles in signal transduction in plants under abiotic stresses. Curr. Genom. 18: 483–497.

25 

Irizarry R.A., Hobbs B., Collin F., Beazer-Barclay Y.D., Antonellis K.J., Scherf U., Speed T.P. (2003) Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostat. 4: 249–264.

26 

Jiang Y., Deyholos M.K. (2006) Comprehensive transcriptional profiling of NaCl-stressed Arabidopsis roots reveals novel classes of responsive genes. BMC Plant Biol. 6: 25.

27 

Jiang D., Zhou L., Chen W., Ye N., Xia J., Zhuang C. (2019) Overexpression of a microRNA-targeted NAC transcription factor improves drought and salt tolerance in Rice via ABA-mediated pathways. Rice (N.Y.) 12(1): 76.

28 

Jin J.P., Tian F., Yang D.C., Meng Y.Q., Kong L., Luo J.C., Gao G. (2017) PlantTFDB 4.0: toward a central hub for transcription factors and regulatory interactions in plants. Nucl. Acids Res. 45: D1040–D1045.

29 

Johnson W.E., Li C., Rabinovic A. (2007) Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostat. 8: 118–127.

30 

Jung J., Won S.Y., Suh S.C., Kim H., Wing R., Jeong Y., Hwang I., Kim M. (2007) The barley ERF-type transcription factor HvRAF confers enhanced pathogen resistance and salt tolerance in Arabidopsis. Planta 225: 575–588.

31 

Kantar M., Unver T., Budak H. (2010) Regulation of barley miRNAs upon dehydration stress correlated with target gene expression. Funct. Integ. Genomics 10: 493–507.

32 

Khan A., Fornes O., Stigliani A., Gheorghe M., Castro-Mondragon J.M., van der Lee R., Bessy A., Chèneby J., Kulkarni S.R., Tan G., et al. (2018) JASPAR 2018: update of the open-access database of transcription factor binding profiles and its web framework. Nucl. Acids Res. 46: D260–D266.

33 

Kim Y., Park C., Cha S., Han M., Ryu K.S., Suh J.Y. (2020) Determinants of PB1 domain interactions in auxin response factor ARF5 and repressor IAA17. J. Mol. Biol. 432(14): 4010–4022.

34 

Kumar R.R., Pathak H., Sharma S.K., Kala Y.K., Nirjal M.K., Singh G.P., Goswami S., Rai R.D. (2014) Novel and conserved heat-responsive microRNAs in wheat (Triticum aestivum L.). Funct. Integr. Genomics 15: 323–348.

35 

Le D.T., Nishiyama R., Watanabe Y., Mochida K., Yamaguchi-Shinozaki K., Shinozaki K., Tran L.S. (2011) Genome-wide survey and expression analysis of the plant specific NAC transcription factor family in soybean during development and dehydration stress. DNA Res. 18: 263–276.

36 

Lee M.C., Park J.C., Kim D.H., Kang S., Shin K.H., Park H.G., Han j., Lee J.S. (2017) Interrelationship of salinity shift with oxidative stress and lipid metabolism in the monogonont rotifer Brachionus koreanus. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 214: 79–84.

37 

Leek J.T., Johnson W.E., Parker H.S., Jaffe A.E., Storey J.D. (2012) The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics 28: 882–883.

38 

Legris M., Ince Y.Ç., Fankhauser C. (2019) Molecular mechanisms underlying phytochrome-controlled morphogenesis in plants. Nat. Commun. 10: 1–15.

39 

Li T., Kang X., Lei W., Yao X., Zou L., Zhang D., Lin H. (2020) SHY2 as a node in the regulation of root meristem development by auxin, brassinosteroids, and cytokinin. J. Integr. Plant Biol. 62(10): 1500–1517. doi.org 10.1111/jipb.12931

40 

Li W.X., Oono Y., Zhu J., He X.J., Wu J.M., Iida K., Lu X.Y., Cui X., Jin H., Zhu J.K. (2008) The Arabidopsis NFYA5 transcription factor is regulated transcriptionally and post-transcriptionally to promote drought resistance. Plant Cell 20: 2238–2251.

41 

Lippold F., Sanchez D.H., Musialak M., Schlereth A., Scheible W.R., Hincha D.K., Udvardi M.K. (2009) AtMyb41 regulates transcriptional and metabolic responses to osmotic stress in Arabidopsis. Plant Physiol. 149: 1761–1772.

42 

Liu D., Leib K., Zhao P., Kogel K.H., Langen G. (2014) Phylogenetic analysis of barley WRKY proteins and characterization of HvWRKY1 and-2 as repressors of the pathogen-inducible gene HvGER4c. Mol. Genet. Genomics 289: 1331–1345.

43 

Luan M., Xu M., Lu Y., Zhang L., Fan Y., Wang L. (2015) Expression of zma-miR169 miRNAs and their target ZmNFYA genes in response to abiotic stress in maize leaves. Gene 555: 178–185.

44 

Mao X., Jia D., Li A., Zhang H., Tian S., Zhang X., Jia J., Jing R. (2011) Transgenic expression of TaMYB2A confers enhanced tolerance to multiple abiotic stresses in Arabidopsis. Funct. Integr. Gen. 11: 445–465.

45 

Mare C., Mazzucotelli E., Crosatti C., Francia E., Cattivelli L. (2004) Hv-WRKY38: a new transcription factor involved in cold-and drought response in barley. Plant Mol. Biol. 55: 399–416.

46 

McLeay R.C., Bailey T.L. (2010) Motif enrichment analysis: a unified framework and an evaluation on ChIP data. BMC Bioinformatics 11: 165.

47 

Munns R., James R., Lauchli A. (2006) Approaches to increasing the salt tolerance of wheat and other cereals. J. Exp. Bot. 57: 1025–1043.

48 

Muñoz Amatriaín M., Svensson J.T., Castillo A.M., Cistue L., Close T.J., Valles M.P. (2006) Transcriptome analysis of barley anthers: effect of mannitol treatment on microspore embryogenesis. Physiol. Plant 127: 551–560.

49 

Mwando E., Han Y., Angessa T.T., Zhou G., Hill C.B., Zhang X.Q., Li C. (2020) Genome-wide association study of salinity tolerance during germination in barley (Hordeum vulgare L.). Front. Plant Sci. 11: 1–15.

50 

Newton A.C., Flavell A.J., George T.S., Leat P.M.K., Mullholland B., Ramsay L., Revoredo-Giha C., Russell J., Steffenson B.J., Swanston J.S. et al. (2011) Crops that feed the world 4. Barley: a resilient crop? Strengths and weaknesses in the context of food security. Food Sec. 3: 141.

51 

Ni Z., Hu Z., Jiang Q., Zhang H. (2013) GmNFYA3, a target gene of miR169, is a positive regulator of plant tolerance to drought stress. Plant Mol. Biol. 82: 113–129.

52 

Parmar S., Gharat S.A., Tagirasa R., Chandra T., Behera L., Dash S.K., Shaw B.P. (2020) Identification and expression analysis of miRNAs and elucidation of their role in salt tolerance in rice varieties susceptible and tolerant to salinity. PLoS One 15(4): e0230958.

53 

Pegler J.L., Oultram J.M., Grof C.P., Eamens A.L. (2019) Profiling the abiotic stress responsive microRNA landscape of Arabidopsis thaliana. Plants 8(3): 58.

54 

Pourabed E., Ghane Golmohamadi F., Soleymani Monfared P., Razavi S.M., Shobbar Z.S. (2015) Basic leucine zipper family in barley: genome-wide characterization of members and expression analysis. Mol. Biotechnol. 57: 12–26.

55 

Qiao L., Zhang L., Zhang X., Zhang L., Li X., Chang J., Zhan H., Guo H., Zheng J., Chang Z. (2017) Evolution of the Aux/IAA gene family in hexaploid wheat. J. Mol. Evol. 85(3–4): 107–119.

56 

Rani V. (2007) Computational methods to dissect cis-regulatory transcriptional networks. J. Biol. Sci. 32: 1325–1330.

57 

Rao S., Balyan S., Jha, S. Mathur S. (2020) Novel insights into expansion and functional diversification of MIR169 family in tomato. Planta 251(2): 55.

58 

Reed J.W., Nagatani A., Elich T.D., Fagan M., Chory J. (1994) Phytochrome A and phytochrome B have overlapping but distinct functions in Arabidopsis development. Plant Physiol. 104: 1139–1149.

59 

Sallam A., Alqudah A.M., Dawood M.F., Baenziger P.S., Börner A. (2019) Drought stress tolerance in wheat and barley: advances in physiology, breeding and genetics research. Int. J. Mol. 20: 3137–3143.

60 

Schmidt R., Mieulet D., Hubberten H.M., Obata T., Hoefgen R., Fernie A.R., Fisahn J., San Segundo B., Guiderdoni E., Schippers J.H. et al. (2013) Salt-responsive ERF1 regulates reactive oxygen species-dependent signaling during the initial response to salt stress in rice. Plant Cell 25: 2115–2131.

61 

Shi Q., Zhang Y., To V.T., Shi J., Zhang D., Cai W. (2020) Genome-wide characterization and expression analyses of the auxin/indole-3-acetic acid (Aux/IAA) gene family in barley (Hordeum vulgare L.). Sci. Rep. 10(1): 10242.

62 

Sinha S.K., Rani M., Bansal N., Gayatri Venkatesh K., Mandal P.K. (2015) Nitrate starvation induced changes in root system architecture, carbon: nitrogen metabolism, and miRNA expression in nitrogen-responsive wheat genotypes. Appl. Biochem. Biotechnol. 177: 1299–1312.

63 

Sorin C., Declerck M., Christ A., Blein T., Ma L., Lelandais-Brière C., Njo M.F., Beeckman T., Crespi M., Hartmann C. (2014) A miR169 isoform regulates specific NF-YA targets and root architecture in Arabidopsis. New Phytol. 202: 1197–1211.

64 

Thabet S.G., Moursi Y.S., Karam M.A., Börner A., Alqudah A.M. (2020) Natural variation uncovers candidate genes for barley spikelet number and grain yield under drought stress. Genes 11: 533.

65 

Tian F., Yang D.C., Meng Y.Q., Jin J.P., Gao G. (2020) Plant-RegMap: charting functional regulatory maps in plants. Nucl. Acids Res. 48: D1104–D1113.

66 

Wang G., Luo, X., Wang L., Wan J., Xia, S., Zhu, H., Qian, J., Wang, Y. (2018) MeDReaders: a database for transcription factors that bind to methylated DNA. Nucl. Acids Res. 46: D146–D151.

67 

Wu J., Chen J., Wang L., Wang S. (2017) Genome-wide investigation of WRKY transcription factors involved in terminal drought stress response in common bean. Front. Plant Sci. 8: 380.

68 

Xie Z., Nolan T.M., Jiang H., Yin Y. (2019) AP2/ERF transcription factor regulatory networks in hormone and abiotic stress responses in Arabidopsis. Front. Plant. Sci. 10: 228.

69 

Yamaguchi-Shinozaki K., Shinozaki K. (2005) Organization of cis-acting regulatory elements in osmotic- and cold-stress-responsive promoters. Trends Plant Sci. 10: 88–94.

70 

Yang Y., Wang Q., Chen Q., Yin X., Qian M., Sun X., Yang Y. (2017a) Genome-wide survey indicates diverse physiological roles of the barley (Hordeum vulgare L.) calcium-dependent protein kinase genes. Sci. Rep. 7: 5306.

71 

Yang Z., Dai Z., Lu R., Wu B., Tang Q., Xu Y., Cheng C., Su J. (2017b) Transcriptome analysis of two species of jute in response to polyethylene glycol (PEG)- induced drought stress. Sci. Rep. 7: 16565.

72 

Yoon Y., Seo D.H., Shin H., Kim H.J., Kim C.M., Jang G. (2020) The role of stress-responsive transcription factors in modulating abiotic stress tolerance in plants. Agronomy 10(6): 788.

73 

Zare S., Nazarian-Firouzabadi F., Ismailia A., Pakniyat H. (2019) Identification of miRNAs and evaluation of candidate genes expression profile associated with drought stress in barley. Plant Gene 20: 100205.

74 

Zhang L., Song Z., Li F., Li X., Ji H., Yang S. (2019) The specific MYB binding sites bound by TaMYB in the GAPCp2/3 promoters are involved in the drought stress response in wheat. BMC Plant Biol. 19: 366–370.

75 

Zheng Y., Jiao C., Sun H., Rosli H.G., Pombo M.A., Zhang P., Banf M., Dai X., Martin G.B., Giovannoni J.J., et al. (2016) iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. Mol. Plant 9: 1667–1670.

76 

Zhou M., Tang W. (2019) MicroRNA156 amplifies transcription factor-associated cold stress tolerance in plant cells. Mol. Genet. Genomics 294(2): 379–393.

77 

Zhu X., Qi L., Liu X., Cai S., Xu H., Huang R., Li J., Wei X., Zhang Z. (2014) The wheat ethylene response factor transcription factor pathogen-induced ERF1 mediates host responses to both the necrotrophic pathogen Rhizoctonia cerealis and freezing stresses. Plant Physiol. 164: 1499–1514.

78 

Zhu Y.N., Shi D.Q., Ruan M.B., Zhang, L.L., Meng Z.H., Liu J., Yang W.C. (2013) Transcriptome analysis reveals crosstalk of responsive genes to multiple abiotic stresses in cotton (Gossypium hirsutum L.). PLoS One 8: e80218.

79 

Zou C., Sun K., Mackaluso J.D., Seddon A.E., Jin R., Thomashow M.F., Shiu S.H. (2011) Cis-regulatory code of stress-responsive transcription in Arabidopsis thaliana. Proc. Natl. Acad. Sci. USA 108: 14992–14997.

Copyright: © 2021 Institute of Bioorganic Chemistry, Polish Academy of Sciences This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs (CC BY-NC-ND) (https://creativecommons.org/licenses/by-nc-nd/3.0/legalcode),.allowing third parties to download and share its works but not commercially purposes or to create derivative works.
 
Stosujemy się do standardu HONcode dla wiarygodnej informacji zdrowotnej This site complies with the HONcode standard for trustworthy health information: verify here