Wide-ranging transcriptome remodelling mediated by alternative polyadenylation in response to abiotic stresses in Sorghum

Full-length transcriptome sequencing provides insights into alternative splicing under cold stress in peanut

Introduction

Peanut (Arachis hypogaea L.), also called groundnut is an important oil and cash crop grown widely in the world. The annual global production of groundnuts has increased to approximately 50 million tons, which provides a rich source of vegetable oils and proteins for humans. Low temperature (non-freezing) is one of the major factors restricting peanut growth, yield, and geographic distribution. Since the complexity of cold-resistance trait, the molecular mechanism of cold tolerance and related gene networks were largely unknown in peanut.

Methods

In this study, comparative transcriptomic analysis of two peanut cultivars (SLH vs. ZH12) with differential cold tolerance under low temperature (10°C) was performed using Oxford Nanopore Technology (ONT) platform.

Results and discussion

As a result, we identified 8,949 novel gene loci and 95,291 new/novel isoforms compared with the reference database. More differentially expressed genes (DEGs) were discovered in cold-sensitive cultivar (ZH12) than cold-tolerant cultivar (SLH), while more alternative splicing events were found in SLH compared to ZH12. Gene Ontology (GO) analyses of the common DEGs showed that the "response to stress", "chloroplast part", and "transcription factor activity" were the most enriched GO terms, indicating that photosynthesis process and transcription factors play crucial roles in cold stress response in peanut. We also detected a total of 708 differential alternative splicing genes (DASGs) under cold stress compared to normal condition. Intron retention (IR) and exon skipping (ES) were the most prevalent alternative splicing (AS) events. In total, 4,993 transcription factors and 292 splicing factors were detected, many of them had differential expression levels and/or underwent AS events in response to cold stress. Overexpression of two candidate genes (encoding trehalose-6-phosphatephosphatases, AhTPPs) in yeast improves cold tolerance. This study not only provides valuable resources for the study of cold resistance in peanut but also lay a foundation for genetic modification of cold regulators to enhance stress tolerance in crops.

Alternative polyadenylation profiles of susceptible and resistant rice (Oryza sativa L.) in response to bacterial leaf blight using RNA-seq

BACKGROUND

Alternative polyadenylation (APA) is an important pattern of post-transcriptional regulation of genes widely existing in eukaryotes, involving plant physiological and pathological processes. However, there is a dearth of studies investigating the role of APA profile in rice leaf blight.

RESULTS

In this study, we compared the APA profile of leaf blight-susceptible varieties (CT 9737-613P-M) and resistant varieties (NSIC RC154) following bacterial blight infection. Through gene enrichment analysis, we found that the genes of two varieties typically exhibited distal poly(A) (PA) sites that play different roles in two kinds of rice, indicating differential APA regulatory mechanisms. In this process, many disease-resistance genes displayed multiple transcripts via APA. Moreover, we also found five polyadenylation factors of similar expression patterns of rice, highlighting the critical roles of these five factors in rice response to leaf blight about PA locus diversity.

CONCLUSION

Notably, the present study provides the first dynamic changes of APA in rice in early response to biotic stresses and proposes a possible functional conjecture of APA in plant immune response, which lays the theoretical foundation for in-depth determination of the role of APA events in plant stress response and other life processes.

Natural variation in the plant polyadenylation complex

Messenger RNA polyadenylation, the process wherein the primary RNA polymerase II transcript is cleaved and a poly(A) tract added, is a key step in the expression of genes in plants. Moreover, it is a point at which gene expression may be regulated by determining the functionality of the mature mRNA. Polyadenylation is mediated by a complex (the polyadenylation complex, or PAC) that consists of between 15 and 20 subunits. While the general functioning of these subunits may be inferred by extending paradigms established in well-developed eukaryotic models, much remains to be learned about the roles of individual subunits in the regulation of polyadenylation in plants. To gain further insight into this, we conducted a survey of variability in the plant PAC. For this, we drew upon a database of naturally-occurring variation in numerous geographic isolates of Arabidopsis thaliana. For a subset of genes encoding PAC subunits, the patterns of variability included the occurrence of premature stop codons in some Arabidopsis accessions. These and other observations lead us to conclude that some genes purported to encode PAC subunits in Arabidopsis are actually pseudogenes, and that others may encode proteins with dispensable functions in the plant. Many subunits of the PAC showed patterns of variability that were consistent with their roles as essential proteins in the cell. Several other PAC subunits exhibit patterns of variability consistent with selection for new or altered function. We propose that these latter subunits participate in regulatory interactions important for differential usage of poly(A) sites.

Smart reprograming of plants against salinity stress using modern biotechnological tools

Climate change gives rise to numerous environmental stresses, including soil salinity. Salinity/salt stress is the second biggest abiotic factor affecting agricultural productivity worldwide by damaging numerous physiological, biochemical, and molecular processes. In particular, salinity affects plant growth, development, and productivity. Salinity responses include modulation of ion homeostasis, antioxidant defense system induction, and biosynthesis of numerous phytohormones and osmoprotectants to protect plants from osmotic stress by decreasing ion toxicity and augmented reactive oxygen species scavenging. As most crop plants are sensitive to salinity, improving salt tolerance is crucial in sustaining global agricultural productivity. In response to salinity, plants trigger stress-related genes, proteins, and the accumulation of metabolites to cope with the adverse consequence of salinity. Therefore, this review presents an overview of salinity stress in crop plants. We highlight advances in modern biotechnological tools, such as omics (genomics, transcriptomics, proteomics, and metabolomics) approaches and different genome editing tools (ZFN, TALEN, and CRISPR/Cas system) for improving salinity tolerance in plants and accomplish the goal of "zero hunger," a worldwide sustainable development goal proposed by the FAO.

Differential alternative polyadenylation of homoeologous genes of allohexaploid wheat ABD subgenomes during drought stress response

Because allohexaploid wheat genome contains ABD subgenomes, how the expression of homoeologous genes is coordinated remains largely unknown, particularly at the co-transcriptional level. Alternative polyadenylation (APA) is an important part of co-transcriptional regulation, which is crucial in developmental processes and stress responses. Drought stress is a major threat to the stable yield of wheat. Focusing on APA, we used poly(A) tag sequencing to track poly(A) site dynamics in wheat under drought stress. The results showed that drought stress led to extensive APA involving 37-47% of differentially expressed genes in wheat. Significant poly(A) site switching was found in stress-responsive genes. Interestingly, homoeologous genes exhibit unequal numbers of poly(A) sites, divergent APA patterns with tissue specificity and time-course dynamics, and distinct 3'-UTR length changes. Moreover, differentially expressed transcripts in leaves and roots used different poly(A) signals, the up- and downregulated isoforms had distinct preferences for non-canonical poly(A) sites. Genes that encode key polyadenylation factors showed differential expression patterns under drought stress. In summary, poly(A) signals and the changes in core poly(A) factors may widely affect the selection of poly(A) sites and gene expression levels during the response to drought stress, and divergent APA patterns among homoeologous genes add extensive plasticity to this responsive network. These results not only reveal the significant role of APA in drought stress response, but also provide a fresh perspective on how homoeologous genes contribute to adaptability through transcriptome diversity. In addition, this work provides information about the ends of transcripts for a better annotation of the wheat genome.

Alternative Polyadenylation Is a Novel Strategy for the Regulation of Gene Expression in Response to Stresses in Plants

Precursor message RNA requires processing to generate mature RNA. Cleavage and polyadenylation at the 3'-end in the maturation of mRNA is one of key processing steps in eukaryotes. The polyadenylation (poly(A)) tail of mRNA is an essential feature that is required to mediate its nuclear export, stability, translation efficiency, and subcellular localization. Most genes have at least two mRNA isoforms via alternative splicing (AS) or alternative polyadenylation (APA), which increases the diversity of transcriptome and proteome. However, most previous studies have focused on the role of alternative splicing on the regulation of gene expression. In this review, we summarize the recent advances concerning APA in the regulation of gene expression and in response to stresses in plants. We also discuss the mechanisms for the regulation of APA for plants in the adaptation to stress responses, and suggest that APA is a novel strategy for the adaptation to environmental changes and response to stresses in plants.

Review: Mechanisms underlying alternative polyadenylation in plants - looking in the right places

Recent years have seen an explosion of interest in the subject of alternative polyadenylation in plants. Connections between the polyadenylation complex and numerous developmental and stress responses are well-established. However, those that link stimuli with the functioning of the polyadenylation complex are less well understood. To this end, it is imperative to clearly delineate the roles of the polyadenylation complex in both plant growth AND alternative polyadenylation. It is also necessary to understand the ways by which other molecular processes may contribute to alternative polyadenylation. This review discusses these issues, with a focus on instances that reveal mechanisms by which mRNA polyadenylation may be regulated. Insights from from characterizations of mutants affected in the polyadenylation complex are discussed, as are the limitations of such characterizations when it comes to teasing out cause and effect. These limitations encourage explorations to other processes that are beyond the core polyadenylation complex. Two such processes that sculpt the plant transcriptome - transcription termination and the epigenetic control of transposon activity - also contribute to regulated poly(A) site choice. These subjects define "the right places" - molecular mechanisms that contribute to the wide-ranging control of gene expression via mRNA polyadenylation.

Drought induces epitranscriptome and proteome changes in stem-differentiating xylem of Populus trichocarpa

Understanding gene expression and regulation requires insights into RNA transcription, processing, modification, and translation. However, the relationship between the epitranscriptome and the proteome under drought stress remains undetermined in poplar (Populus trichocarpa). In this study, we used Nanopore direct RNA sequencing and tandem mass tag-based proteomic analysis to examine epitranscriptomic and proteomic regulation induced by drought treatment in stem-differentiating xylem (SDX). Our results revealed a decreased full-length read ratio under drought treatment and, especially, a decreased association between transcriptome and proteome changes in response to drought. Epitranscriptome analysis of cellulose- and lignin-related genes revealed an increased N6-Methyladenosine (m6A) ratio, which was accompanied by decreased RNA abundance and translation, under drought stress. Interestingly, usage of the distal poly(A) site increased during drought stress. Finally, we found that transcripts of highly expressed genes tend to have shorter poly(A) tail length (PAL), and drought stress increased the percentage of transcripts with long PAL. These findings provide insights into the interplay among m6A, polyadenylation, PAL, and translation under drought stress in P. trichocarpa SDX.

Phased small RNA-mediated systemic signaling in plants

Systemic acquired resistance (SAR) involves the generation of systemically transported signal that arms distal plant parts against secondary infections. We show that two phased 21-nucleotide (nt) trans-acting small interfering RNA3a RNAs (tasi-RNA) derived from TAS3a and synthesized within 3 hours of pathogen infection are the early mobile signal in SAR. TAS3a undergoes alternate polyadenylation, resulting in the generation of 555- and 367-nt transcripts. The 555-nt transcripts likely serves as the sole precursor for tasi-RNAs D7 and D8, which cleave Auxin response factors (ARF) 2, 3, and 4 to induce SAR. Conversely, increased expression of ARF3 represses SAR. Knockout mutations in TAS3a or RNA silencing components required for tasi-RNA biogenesis compromise SAR without altering levels of known SAR-inducing chemicals. Both tasi-ARFs and the 367-nt transcripts are mobile and transported via plasmodesmata. Together, we show that tasi-ARFs are the early mobile signal in SAR.

Systematical characterization of GRF gene family in sorghum, and their potential functions in aphid resistance

Sorghum (Sorghum bicolor) is the fifth important cereal and an industrial energy crop in the world. Growth Regulation Factors (GRFs) play an important role in response to environmental stress, however, the knowledge of GRFs relating to the pest resistance is lacking. Here, we identified 8 GRF genes harboring the typical QLQ (glutamine, leucine, glutamine) and WRC (tryptophan, arginine, cysteine) domains in Sorghum, which could be classified into 4 clades through phylogenetic analysis. The SbGRF genes express in most tissues, while more than half of them express at the highest level in inflorescence. To further investigate their possible role in stress response, we analyzed the transcriptomics data. The results showed that SbGRFs could respond to the abiotic stresses including heat, salt and drought stress. Furthermore, combined the data with qRT-PCR, SbGRF1, 2, 4 and 7 were identified as dominant genes response to the aphid-induced stress. SSR markers close to these genes were also searched. Above all, we summarized the SbGRFs and provided their potential roles in aphid response.

Full-Length Transcriptional Analysis of the Same Soybean Genotype With Compatible and Incompatible Reactions to Heterodera glycines Reveals Nematode Infection Activating Plant Defense Response

Full-length transcriptome sequencing with long reads is a powerful tool to analyze transcriptional and post-transcriptional events; however, it has not been applied on soybean (Glycine max). Here, a comparative full-length transcriptome analysis was performed on soybean genotype 09-138 infected with soybean cyst nematode (SCN, Heterodera glycines) race 4 (SCN4, incompatible reaction) and race 5 (SCN5, compatible reaction) using Oxford Nanopore Technology. Each of 9 full-length samples collected 8 days post inoculation with/without nematodes generated an average of 6.1 GB of clean data and a total of 65,038 transcript sequences. After redundant transcripts were removed, 1,117 novel genes and 41,096 novel transcripts were identified. By analyzing the sequence structure of the novel transcripts, a total of 28,759 complete open reading frame (ORF) sequences, 5,337 transcription factors, 288 long non-coding RNAs, and 40,090 novel transcripts with function annotation were predicted. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of differentially expressed genes (DEGs) revealed that growth hormone, auxin-activated signaling pathway and multidimensional cell growth, and phenylpropanoid biosynthesis pathway were enriched by infection with both nematode races. More DEGs associated with stress response elements, plant-hormone signaling transduction pathway, and plant-pathogen interaction pathway with more upregulation were found in the incompatible reaction with SCN4 infection, and more DEGs with more upregulation involved in cell wall modification and carbohydrate bioprocess were detected in the compatible reaction with SCN5 infection when compared with each other. Among them, overlapping DEGs with a quantitative difference was triggered. The combination of protein-protein interaction with DEGs for the first time indicated that nematode infection activated the interactions between transcription factor WRKY and VQ (valine-glutamine motif) to contribute to soybean defense. The knowledge of the SCN-soybean interaction mechanism as a model will present more understanding of other plant-nematode interactions.

Translational Landscape of a C4 Plant, Sorghum bicolor, Under Normal and Sulfur-Deficient Conditions

Recent accumulation of genomic and transcriptomic information has facilitated genetic studies. Increasing evidence has demonstrated that translation is an important regulatory step, and the transcriptome does not necessarily reflect the profile of functional protein production. Deep sequencing of ribosome-protected mRNA fragments (ribosome profiling or Ribo-seq) has enabled genome-wide analysis of translation. Sorghum is a C4 cereal important not only as food but also as forage and a bioenergy resource. Its resistance to harsh environments has made it an agriculturally important research subject. Yet genome-wide translational profiles in sorghum are still missing. In this study, we took advantage of Ribo-seq and identified actively translated reading frames throughout the genome. We detected translation of 4,843 main open reading frames (ORFs) annotated in the sorghum reference genome version 3.1 and revealed a number of unannotated translational events. A comparison of the transcriptome and translatome between sorghums grown under normal and sulfur-deficient conditions revealed that gene expression is modulated independently at transcript and translation levels. Our study revealed the translational landscape of sorghum's response to sulfur and provides datasets that could serve as a fundamental resource to extend genetic research on sorghum, including studies on translational regulation.

Biological Parts for Engineering Abiotic Stress Tolerance in Plants

It is vital to ramp up crop production dramatically by 2050 due to the increasing global population and demand for food. However, with the climate change projections showing that droughts and heatwaves becoming common in much of the globe, there is a severe threat of a sharp decline in crop yields. Thus, developing crop varieties with inbuilt genetic tolerance to environmental stresses is urgently needed. Selective breeding based on genetic diversity is not keeping up with the growing demand for food and feed. However, the emergence of contemporary plant genetic engineering, genome-editing, and synthetic biology offer precise tools for developing crops that can sustain productivity under stress conditions. Here, we summarize the systems biology-level understanding of regulatory pathways involved in perception, signalling, and protective processes activated in response to unfavourable environmental conditions. The potential role of noncoding RNAs in the regulation of abiotic stress responses has also been highlighted. Further, examples of imparting abiotic stress tolerance by genetic engineering are discussed. Additionally, we provide perspectives on the rational design of abiotic stress tolerance through synthetic biology and list various bioparts that can be used to design synthetic gene circuits whose stress-protective functions can be switched on/off in response to environmental cues.

Divergence in the Regulation of the Salt Tolerant Response Between Arabidopsis thaliana and Its Halophytic Relative Eutrema salsugineum by mRNA Alternative Polyadenylation

Salt tolerance is an important mechanism by which plants can adapt to a saline environment. To understand the process of salt tolerance, we performed global analyses of mRNA alternative polyadenylation (APA), an important regulatory mechanism during eukaryotic gene expression, in Arabidopsis thaliana and its halophytic relative Eutrema salsugineum with regard to their responses to salt stress. Analyses showed that while APA occurs commonly in both Arabidopsis and Eutrema, Eutrema possesses fewer APA genes than Arabidopsis (47% vs. 54%). However, the proportion of APA genes was significantly increased in Arabidopsis under salt stress but not in Eutrema. This indicated that Arabidopsis is more sensitive to salt stress and that Eutrema exhibits an innate response to such conditions. Both species utilized distal poly(A) sites under salt stress; however, only eight genes were found to overlap when their 3' untranslated region (UTR) lengthen genes were compared, thus revealing their distinct responses to salt stress. In Arabidopsis, genes that use distal poly(A) sites were enriched in response to salt stress. However, in Eutrema, the use of poly(A) sites was less affected and fewer genes were enriched. The transcripts with upregulated poly(A) sites in Arabidopsis showed enriched pathways in plant hormone signal transduction, starch and sucrose metabolism, and fatty acid elongation; in Eutrema, biosynthetic pathways (stilbenoid, diarylheptanoid, and gingerol) and metabolic pathways (arginine and proline) showed enrichment. APA was associated with 42% and 29% of the differentially expressed genes (DE genes) in Arabidopsis and Eutrema experiencing salt stress, respectively. Salt specific poly(A) sites and salt-inducible APA events were identified in both species; notably, some salt tolerance-related genes and transcription factor genes exhibited differential APA patterns, such as CIPK21 and LEA4-5. Our results suggest that adapted species exhibit more orderly response at the RNA maturation step under salt stress, while more salt-specific poly(A) sites were activated in Arabidopsis to cope with salinity conditions. Collectively, our findings not only highlight the importance of APA in the regulation of gene expression in response to salt stress, but also provide a new perspective on how salt-sensitive and salt-tolerant species perform differently under stress conditions through transcriptome diversity.

Abiotic stress responses in plants

Plants cannot move, so they must endure abiotic stresses such as drought, salinity and extreme temperatures. These stressors greatly limit the distribution of plants, alter their growth and development, and reduce crop productivity. Recent progress in our understanding of the molecular mechanisms underlying the responses of plants to abiotic stresses emphasizes their multilevel nature; multiple processes are involved, including sensing, signalling, transcription, transcript processing, translation and post-translational protein modifications. This improved knowledge can be used to boost crop productivity and agricultural sustainability through genetic, chemical and microbial approaches.

Decoding co-/post-transcriptional complexities of plant transcriptomes and epitranscriptome using next-generation sequencing technologies

Next-generation sequencing (NGS) technologies - Illumina RNA-seq, Pacific Biosciences isoform sequencing (PacBio Iso-seq), and Oxford Nanopore direct RNA sequencing (DRS) - have revealed the complexity of plant transcriptomes and their regulation at the co-/post-transcriptional level. Global analysis of mature mRNAs, transcripts from nuclear run-on assays, and nascent chromatin-bound mRNAs using short as well as full-length and single-molecule DRS reads have uncovered potential roles of different forms of RNA polymerase II during the transcription process, and the extent of co-transcriptional pre-mRNA splicing and polyadenylation. These tools have also allowed mapping of transcriptome-wide start sites in cap-containing RNAs, poly(A) site choice, poly(A) tail length, and RNA base modifications. The emerging theme from recent studies is that reprogramming of gene expression in response to developmental cues and stresses at the co-/post-transcriptional level likely plays a crucial role in eliciting appropriate responses for optimal growth and plant survival under adverse conditions. Although the mechanisms by which developmental cues and different stresses regulate co-/post-transcriptional splicing are largely unknown, a few recent studies indicate that the external cues target spliceosomal and splicing regulatory proteins to modulate alternative splicing. In this review, we provide an overview of recent discoveries on the dynamics and complexities of plant transcriptomes, mechanistic insights into splicing regulation, and discuss critical gaps in co-/post-transcriptional research that need to be addressed using diverse genomic and biochemical approaches.

Alternative Polyadenylation in response to temperature stress contributes to gene regulation in Populus trichocarpa

BACKGROUND

Genome-wide change of polyadenylation (polyA) sites (also known as alternative polyadenylation, APA) is emerging as an important strategy of gene regulation in response to stress in plants. But little is known in woody perennials that are persistently dealing with multiple abiotic stresses.

RESULTS

Here, we performed a genome-wide profiling of polyadenylation sites under heat and cold treatments in Populus trichocarpa. Through a comprehensive analysis of polyA tail sequences, we identified 25,919 polyA-site clusters (PACs), and revealed 3429 and 3139 genes shifted polyA sites under heat and cold stresses respectively. We found that a small proportion of genes possessed APA that affected the open reading frames; and some shifts were commonly identified. Functional analysis of genes displaying shifted polyA tails suggested that pathways related to RNA metabolism were linked to regulate the APA events under both heat and cold stresses. Interestingly, we found that the heat stress induced a significantly more antisense PACs comparing to cold and control conditions. Furthermore, we showed that a unique cis-element (AAAAAA) was predominately enriched downstream of PACs in P. trichocarpa genes; and this sequence signal was only absent in shifted PACs under the heat condition, indicating a distinct APA mechanism responsive to heat tolerance.

CONCLUSIONS

This work provides a comprehensive picture of global polyadenylation patterns in response to temperatures stresses in trees. We show that the frequent change of polyA tail is a potential mechanism of gene regulation responsive to stress, which are associated with distinctive sequence signatures.

PSDX: A Comprehensive Multi-Omics Association Database of Populus trichocarpa With a Focus on the Secondary Growth in Response to Stresses

Populus trichocarpa (P. trichocarpa) is a model tree for the investigation of wood formation. In recent years, researchers have generated a large number of high-throughput sequencing data in P. trichocarpa. However, no comprehensive database that provides multi-omics associations for the investigation of secondary growth in response to diverse stresses has been reported. Therefore, we developed a public repository that presents comprehensive measurements of gene expression and post-transcriptional regulation by integrating 144 RNA-Seq, 33 ChIP-seq, and six single-molecule real-time (SMRT) isoform sequencing (Iso-seq) libraries prepared from tissues subjected to different stresses. All the samples from different studies were analyzed to obtain gene expression, co-expression network, and differentially expressed genes (DEG) using unified parameters, which allowed comparison of results from different studies and treatments. In addition to gene expression, we also identified and deposited pre-processed data about alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI). The post-transcriptional regulation, differential expression, and co-expression network datasets were integrated into a new P. trichocarpa Stem Differentiating Xylem (PSDX) database (http://forestry.fafu.edu.cn/db/SDX), which further highlights gene families of RNA-binding proteins and stress-related genes. The PSDX also provides tools for data query, visualization, a genome browser, and the BLAST option for sequence-based query. Much of the data is also available for bulk download. The availability of PSDX contributes to the research related to the secondary growth in response to stresses in P. trichocarpa, which will provide new insights that can be useful for the improvement of stress tolerance in woody plants.

Transcriptome Analysis of Drought-Resistant and Drought-Sensitive Sorghum (Sorghum bicolor) Genotypes in Response to PEG-Induced Drought Stress

Drought is a major limiting factor of crop yields. In response to drought, plants reprogram their gene expression, which ultimately regulates a multitude of biochemical and physiological processes. The timing of this reprogramming and the nature of the drought-regulated genes in different genotypes are thought to confer differential tolerance to drought stress. Sorghum is a highly drought-tolerant crop and has been increasingly used as a model cereal to identify genes that confer tolerance. Also, there is considerable natural variation in resistance to drought in different sorghum genotypes. Here, we evaluated drought resistance in four genotypes to polyethylene glycol (PEG)-induced drought stress at the seedling stage and performed transcriptome analysis in seedlings of sorghum genotypes that are either drought-resistant or drought-sensitive to identify drought-regulated changes in gene expression that are unique to drought-resistant genotypes of sorghum. Our analysis revealed that about 180 genes are differentially regulated in response to drought stress only in drought-resistant genotypes and most of these (over 70%) are up-regulated in response to drought. Among these, about 70 genes are novel with no known function and the remaining are transcription factors, signaling and stress-related proteins implicated in drought tolerance in other crops. This study revealed a set of drought-regulated genes, including many genes encoding uncharacterized proteins that are associated with drought tolerance at the seedling stage.