Next Article in Journal
Glutamine Addiction and Therapeutic Strategies in Lung Cancer
Next Article in Special Issue
Evaluation of Genomic Prediction for Pasmo Resistance in Flax
Previous Article in Journal
Dysregulation of Dopaminergic Regulatory Factors TH, Nurr1, and Pitx3 in the Ventral Tegmental Area Associated with Neuronal Injury Induced by Chronic Morphine Dependence
Previous Article in Special Issue
Identification and Characterization of the EXO70 Gene Family in Polyploid Wheat and Related Species
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 2
10.3390/ijms20020251
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Characterization and Identification of Trihelix Transcription Factor and Expression Profiling in Response to Abiotic Stresses in Rice (Oryza sativa L.)

by
Jiaming Li
1,
Minghui Zhang
2,
Jian Sun
1,
Xinrui Mao
1,
Jing Wang
3,
Jingguo Wang
1,
Hualong Liu
1,
Hongliang Zheng
1,
Zhen Zhen
2,
Hongwei Zhao
1 and
Detang Zou
1,*
1
College of Agriculture, Northeast Agricultural University, Harbin 150030, China
2
College of Life Science, Northeast Agricultural University, Harbin 150030, China
3
Agriculture Technology and Popularization Center, Jixi 158100, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(2), 251; https://doi.org/10.3390/ijms20020251
Submission received: 29 October 2018 / Revised: 21 December 2018 / Accepted: 6 January 2019 / Published: 10 January 2019
(This article belongs to the Special Issue Plant Genomics)
Browse Figures
Versions Notes

Abstract

:
Trihelix transcription factors play a role in plant growth, development and various stress responses. Here, we identified 41 trihelix family genes in the rice genome. These OsMSLs (Myb/SANT-LIKE) were located on twelve chromosomes. Synteny analysis indicated only six duplicated gene pairs in the rice trihelix family. Phylogenetic analysis of these OsMSLs and the trihelix genes from other species divided them into five clusters. OsMSLs from different groups significantly diverged in terms of gene structure and conserved functional domains. However, all OsMSLs contained the same five cis-elements. Some of these were responsive to light and dehydration stress. All OsMSLs expressed in four tissues and six developmental stages of rice but with different expression patterns. Quantitative real-time PCR analysis revealed that the OsMSLs responded to abiotic stresses including drought and high salt stress and stress signal molecule including ABA (abscisic acid), hydrogen peroxide. OsMSL39 were simultaneously expressed under all treatments, while OsMSL28 showed high expression under hydrogen peroxide, drought, and high salt treatments. Moreover, OsMSL16/27/33 displayed significant expression under ABA and drought treatments. Nevertheless, their responses were regulated by light. The expression levels of the 12 chosen OsMSLs differed between light and dark conditions. In conclusion, our results helped elucidate the biological functions of rice trihelix genes and provided a theoretical basis for further characterizing their biological roles in responding to abiotic stresses.

1. Introduction

Transcription factors are ubiquitous in plants. They play crucial roles in various growth and development processes and respond to abiotic stresses [1]. Previous studies reported more than 60 transcription factor families in plants [2,3]. However, little is known about several important transcription factor families. Trihelix transcription factors occur only in plants. They were first identified and isolated from pea (Pisum sativum) in the 1990s. They bind to the core sequence of 5’-G-Pu-(T/A)-A-A-(T/A)-3’ of the promoter region of rbcS-3A gene to regulate light-dependent expression [4]. They were initially called GT factors because they bind to light-responsive GT elements. The DNA-binding domain of the GT factors has a typical tandem trihelix (helix-loop-helix-loop-helix) structure which was later renamed the trihelix transcription factor. Subsequent research revealed that the trihelix structure of the GT factors resembles the solution structure of the Myb/SANT-LIKE DNA-binding domain [5]. GT factors evolved from Myb/SANT-LIKE proteins in plants. Gaps between helix pairs created different recognition sequences between GT factors and Myb/SANT-LIKE proteins [5,6]. According to databases like Pfam, the Myb/SANT-LIKE domain represents the trihelix conserved domain.
Trihelix is a family of transcription factors that have only recently received attention. However, the trihelix genes have been systematically studied mainly in dicotyledonous plants such as Arabidopsis, tomato and chrysanthemum, while almost no research has been carefully carried out in a monocotyledonous plant. In Arabidopsis, 30 GT family members were identified and divided into the GT-1, GT-2, GTγ, SH4, and SIP1 subfamilies named after their founding members [7]. The 96 trihelix proteins of tomato (Solanum lycopersicum) were classified into six subfamilies (clades GT-1, GT-2, SH4, SIP1, GTγ, and GTδ). The GTδ subfamily is apparently missing in Arabidopsis [8]. Most of the trihelix gene subfamily structures vary substantially, especially at the C-terminus. The exceptions are GT1 and GT2.
Earlier studies identified the trihelix family genes as a class of light regulators. Nevertheless, the roles of GT factors in light regulation must be systematically established. In Arabidopsis, the GT1 subfamily genes may participate in salt stress and pathogen responses and their expression was induced by light in 3-d seedlings [9]. In contrast, the rice GT-1 gene RML1 (OsMSL21 in the present study) was repressed by light in etiolated seedlings [10]. The trihelix transcription factors in soybean, GmGT-2A and GmGT-2B, were induced by ABA (abscisic acid), drought, high salt levels, and cold in soybean seedlings [11]. Loss-of-function analysis of GTL1 revealed that gtl1 mutants had fewer stomata than wild type plants. In this way, the former had comparatively lower water loss and higher drought tolerance than the latter [12]. The expression of the rice GTγ clade gene OsGTγ-1 increased 2.5 to 10 times in response to salt stress and was also upregulated by ABA treatment [13]. On the other hand, the expression of several trihelix genes in Chrysanthemum was downregulated by ABA [14]. Trihelix transcription factors are also associated with plant morphogenesis. The trihelix transcription factor PETAL LOSS (PTL) determines the number of petals per flower and sepal fusion in Arabidopsis. The rice SH4 clade gene (SH4) promotes the abscission layer development and function in mature seed peduncles [15]. However, the function of the SH4 clade has not yet been investigated. The Arabidopsis SIP1 genes ASIL1 and ASIL2 downregulated the LEA (Late Embryogenesis Abundant) genes in Arabidopsis seedlings [7]. The trihelix genes also have multiple functions throughout plant development. The molecular mechanisms of their stress responses and their involvement in the signaling pathway require elucidation.
Rice (Oryza sativa L.) is both a major global cereal crop and an important tool in plant research. In this study, we identified 41 rice trihelix genes by the Myb/SANT-LIKE domain using HMM-search in silico. We analyzed their chromosomal distributions, gene synteny, phylogenetic analysis, gene structures, motif compositions, cis-elements, and expression patterns in different tissues, developmental stages, and environmental stress responses. The aim of this study was to analyze the structure and function of rice trihelix genes and phylogenetic relationship between rice trihelix proteins and other species including dicotyledonous and monocotyledonous plant. To establish the role of the trihelix genes’ response to stress, we evaluated their response to abiotic stress factors including drought and high salt, and to stress signal molecules, such as abscisic acid and hydrogen peroxide. Our results provide a theoretical basis for the functional analysis of the rice trihelix family genes especially in abiotic stress responses.

2. Results

2.1. Identification of Trihelix Genes in Rice

The HMM (Hidden Markov Model) for the Myb/SANT-LIKE domain identified 117 gene candidates and a rice-specific Myb/SANT-LIKE domain was built using them. The HMM profile search was performed on the whole rice genome with the rice-specific Myb/SANT-LIKE domain and 79 new candidate genes were found. Only genes with E-value < 0.01 were classified in the trihelix family. Putative genes were verified in the Pfam and InterPro databases to confirm the existence of the complete Myb/SANT-LIKE domains. Finally, 41 trihelix genes were identified.
All trihelix genes mapped onto the rice chromosomes, they were named OsMSL01-OsMSL41 according to the gene distribution order on the chromosomes. OsMSL25 and OsMSL34 have two alternative splicing. “MSL” stands for “Myb/SANT-LIKE”. The characteristics of OsMSLs including the gene MSU_Locus ID, the chromosomes locations, the lengths of the CDS (coding sequence) and amino acid sequences, the number of exons, the protein sizes, and the isoelectric points are summarized in Table 1. OsMSL19 was the smallest protein with 266 amino acids, whereas OsMSL12 was the largest with 882 amino acids. The protein MW (Molecular Weight) ranged from 28.62 kDa to 97.37 kDa. Their predicted isoelectric points varied from 4.45 (OsMSL09) to 11.38 (OsMSL17). Twenty-nine of the trihelix transcription factors were localized in the nucleus, ten in the chloroplast, and two in the peroxisome (OsMSL04 and OsMSL29).

2.2. Chromosomal Distributions and Synteny Analysis of Rice Trihelix Genes

The extraction of the chromosomal information of the OsMSLs identified their chromosomal locations. As shown in Figure 1, all OsMSLs had precise positions in the chromosomes. Each rice chromosome contains ≥1 OsMSL. The OsMSLs are unevenly and non-randomly distributed on 12 chromosomes. Chr2 (chromosomal 2) contains the largest number of OsMSLs (eight) whereas Chr6 contains only one. The first four chromosomes contain 23 trihelix genes while chromosomes 5–12 have on average only 2–3 genes per chromosome. Therefore, OsMSLs are distributed mainly on the first four rice chromosomes. Although Chr2 is relatively short, it contains the most OsMSLs. Chr1 is the longest in rice and also contains numerous OsMSLs. Chr10, the shortest chromosome, contains two OsMSLs. In contrast, Chr6 is longer than Chr10 but contains only one OsMSL. There is no apparent correlation between the chromosome length and OsMSL gene distribution. Moreover, only OsMSL12 and OsMSL13 form gene clusters on Chr2.
Synteny was also used to analyze rice trihelix gene duplication. Chromosomal region within 200 kb containing two or more genes is defined as a tandem duplication event [16]. As shown in Figure 1, four rice trihelix genes (OsMSL09/10 and OsMSL12/13) were clustered into two tandem duplication event regions on rice chromosomal 2. Besides the tandem duplication events, segmental duplications were also investigated by BLASTP and MCScanX methods [17]. Four segmental duplication events with eight rice trihelix genes were also identified, which are located on duplicated segments on chromosomes 1, 2, 4, 5, and 11 (Figure 2). This finding is consistent with the highly divergent, non-conservative evolution of OsMSLs.
To further understand the gene duplication mechanisms of the rice trihelix family, we constructed four comparative syntenic maps of rice associated with four representative species, including one dicots (Arabidopsis) (Figure 3A) and three monocots (Brachypodium distachyon, wheat and maize) (Figure 3B). A total of 23 rice trihelix genes showed a syntenic relationship with those in maize, followed by wheat (21), Brachypodium distachyon (19) and Arabidopsis (2), indicating that in comparison with monocotyledonous plants, rice trihelix genes show a high evolution divergence with dicotyledonous plants. Congruously, previous research reported that 14 pairs of orthologous trihelix genes were found between tomato and Arabidopsis [8]. Some OsMSLs were found to be associated with at least three syntenic gene pairs, such as OsMSL14, OsMSL17, and OsMSL21. These genes may have played a crucial role in the trihelix gene family during evolution. To better understand the evolutionary constraints acting on the trihelix gene family, the Ka/Ks ratios of the trihelix gene pairs were calculated (Tables S1–S5). All segmental and tandem duplicated OsMSL gene pairs, and the majority of orthologous trihelix gene pairs had Ka/Ks < 1, suggesting that the rice trihelix gene family might have experienced a strong purifying selective pressure during evolution.

2.3. Phylogenetic Analysis, Gene Structure, and Motif Composition of Trihelix Genes

To better understand the phylogenetic relationships of trihelix genes, a maximum likelihood phylogenetic tree was built based on the multiple sequence alignment of Myb/SANT-LIKE domains among rice and other species which include dicotyledonous plants such as Arabidopsis, soybean, tomato, chrysanthemum and monocotyledonous plant such as maize, wheat, wild rice, Brachypodium distachyon. As shown in Figure 4, OsMSLs were divided into five subfamilies named SIP1, GTγ, GT, SH4, and GTδ according to the characteristics of their trihelix DNA binding domains. Some genes that have been classified previously such as SIGT-4/7/12/18/36 in tomato [8], CmTH2/6/12/17/19/20 in chrysanthemum [14], GmGT-2A and GmGT-2B in soybean [11] was as a classified marker. The GT clade was the largest subfamily, containing 28 trihelix genes, whereas the SH4 clade was the smallest, consisting of 13 members, indicating that trihelix genes were distributed unevenly in the different clades. All clades consisted of genes both from dicot and monocot species. There is a similar classification in rice which was previously named GTδ in tomato and two tomato trihelix genes SIGT-4 and SIGT-12 have been found in this subfamily. To demonstrate the evolutionary relationships among OsMSLs, we constructed an unrooted phylogenetic tree using the full-length amino acid sequences of the OsMSLs. Of the 43 transcripts of the 41 rice trihelix genes, nine belonged to SIP1, 10 belonged to GTγ, 11 belonged to GT, five belonged to SH4, and eight belonged to GTδ (Figure 5A). Most of the duplicated genes were present in the GTδ classification. The phylogenetic tree of the all MSLs between rice and Arabidopsis was constructed and is shown in Figure S1. However, we found that the GTδ subfamily does not contain Arabidopsis trihelix genes.
To identify the differences between the rice trihelix family genes, we analyzed the OsMSL gene structure by comparing each coding sequence with its corresponding genomic sequence. As shown in Figure 5B, the number of OsMSL exons is discontinuously distributed from 1 through 18. Combining the gene structure with the phylogenetic tree, we found that the OsMSL exon-intron distribution is related to its classification. Closely related genes usually have homologs. Therefore, their gene structures are similar. For example, the OsMSL genome sequences in the SIP1 subfamily have no introns and only one exon. Therefore, the evolution of this gene subfamily is relatively conservative. The genes in the GTδ subfamily have no UTR region and only exons and introns except for OsMSL16 and OsMSL41. In contrast, the structures of the various genes in the GTγ, GT, and SH4 subfamilies are relatively different. These results indicate that although the OsMSLs are subdivided into five families, their genes are relatively conservative.
To determine the functions of the trihelix family genes, the OsMSL motif composition was analyzed by amino acid sequence in the MEME program. Ten motifs with E < 1.8×10−45 were identified. These resemble the MSLs in chrysanthemum. The genes for each subfamily were classified [14]. As shown in Figure 6, except for OsMSL01 and OsMSL09, most trihelix family genes contain motif 1 (Myb-type DNA-binding domain) located at the N-terminus of the amino acid sequence. Motifs 2, 6, and 8 are various trihelix DNA binding domains (WWW, WWF, and WWI). These determine OsMSL classification, structure, and function [18]. As the gene structure analysis indicated, the gene motifs and distribution patterns are closely related to their subfamilies. SIP1 contains motif 8, GTγ contains motif 6, and only OsMSL06 contains an extra motif 8. Both GT and SH4 contain motif 2 but that in SH4 is longer than that in GT. OsMSL09 and OsMSL12 in SH4 also contain an additional motif 8. OsMSL02, OsMSL03, OsMSL04, OsMSL26, OsMSL29, and OsMSL38 in the GT subfamily also contain motif 6. Motif 2 with other functional domains and conservative sequences are contained in the rice-specific GTδ subfamily. Although their functions have yet to be elucidated, they may indicate that the GTδ gene in rice has multiple functions.

2.4. Cis-Element Analysis of Rice Trihelix Genes

To understand the genetic functions, metabolic networks, and regulatory mechanisms of rice trihelix genes, the shared cis-elements in the promoter regions of the OsMSLs were analyzed. The 1500-bp upstream OsMSL sequence was obtained and identified as a hypothetical promoter. The potential shared OsMSL cis-element was scanned and screened out and its distribution and function were analyzed. Two dehydration-responsive- and three light-responsive cis-elements common to all OsMSLs were identified and labeled by different colors in the promoter sequence (Figure 7).
As shown in Figure 7, A (ACGTATERD1) and M (GT1GMSCAM4) element are two dehydration-responsive elements and M is a core element. Therefore, OsMSLs probably participate in dehydration (including drought and salt) stress responses. GB (GATABOX element), GC (GT1CONSENSUS), and I (INRNTPSADB) are three light-responsive elements. They indicate that the OsMSLs family potentially consists of light-inducible/repressible genes. Light responsiveness is typical of the GT factor (now known as the trihelix family gene) and was confirmed in our cis-element study. To verify whether OsMSLs are regulated by light under both normal- and stress conditions, a dark treatment was added to the OsMSL expression analysis.

2.5. Expression Profiles of Trihelix Genes in Rice Tissues and Developmental Stages

The expression profiles of the various rice tissues including the root, stem, leaf, and sheath of four-leaf rice seedlings were investigated (Figure 8). As shown in Figure 8A, all 43 transcripts expressed in all tissues but their expression levels varied greatly in each tissue. Particularly, there are more highly expressing trihelix genes in the leaves and sheaths than the other organs. There are almost no genes with low expression levels and most of the genes remained at high expression levels in the leaves. In contrast, comparatively fewer genes with extremely high expression exist in the stems and none of them express at an extremely high level in the roots.
The OsMSLs were clustered into groups I to V according to their expression characteristics in different tissues (Figure 8B). OsMSLs in groups I and II expressed at lower levels in all tissues, but members in group I displayed relatively higher expression in root tissue. For instance, although OsMSL19 maintained very low expression levels in most tissues, its expression displayed a distinct enrichment in root. In contrast, the genes in group II, especially OsMSL02, OsMSL08, and OsMSL12, remained low in the various tissues. On the other hand, in Arabidopsis, the At5g63420 gene, the orthology of OsMSL12, is highly expressed specifically in the seeds, suggesting that also OsMSL12 could show a dominant expression level in seeds rather than in any of these four tissues [19]. Contrary to the group I, most of the genes in group III expressed at higher levels in the stems or leaves and at lower levels in the roots. However, the transcription levels of OsMSL34a and OsMSL34b in group IV were high in all four tissues. Group V members also displayed tissue specific expression, but expressed at higher levels than group I and III. Besides, the genes in group V, especially OsMSL16, OsMSL27, OsMSL28, OsMSL35, and OsMSL39, expressed extremely highly in the leaves and sheaths, whereas their expression levels in the roots and stems were lower than those of genes in group I and III, respectively. In conclusion, OsMSLs displayed tissue expression specificity, indicating their potential roles in different mechanisms.
We also investigated the OsMSL expression profiles at different rice developmental stages including 7 (S1), 20 (S2), 40 (S3), 80 (S4), 100(S5), and 140(S6) days after sowing. As shown in Figure 9A, the numbers of genes with high- and extremely high expression levels are greater in S2 and S3 than in the other stages. In contrast, no genes expressed at high levels during S4, S5, or S6. The genes were subdivided into five groups according to their expression characteristics at different stages (Figure 9B). OsMSLs in group I, were at comparatively higher transcription levels during the S2 and S3 stages except for OsMSL28, which displayed high expression levels at S1 stage as well. In general, the expression levels of the genes in groups II and III were similar to those in group I, but their expression levels in S2 and S3 were lower than in group I. Notably, the group IV genes had lower expression levels than the other four clusters at almost all developmental stages. Contrary to groups II and III, high gene expression levels were observed mainly at S1 in group V. Besides, with the exception of OsMSL35 and OsMSL06, all other genes in group V expressed at significantly lower levels by S2. In conclusion, OsMSLs play a potential role at early developmental stages.

2.6. Quantitative Real-Time PCR Analysis of Rice Trihelix Genes in Responses to Different Treatments under Normal and Dark Conditions

A prediction of the cis-elements of the OsMSLs suggested that they may participate in rice dehydration stress tolerance and light-mediated signaling pathways. To verify this hypothesis, we subjected rice seedlings to ABA, hydrogen peroxide, drought, and high salt then carefully selected 12 genes expressing positively in the leaves (Figure 10 and Figure 11). On the whole, although trihelix genes were induced by multiple treatments, their expression levels under one treatment were much higher than other treatments. For instance, six OsMSLs (OsMSL25a, OsMSL25b, OsMSL28, OsMSL34a, OsMSL34b, and OsMSL35) were induced by multiple treatments, but the transcript levels significantly increased under hydrogen peroxide treatment compared to other treatments. Several genes were induced after being repressed, such as OsMSL41, which remained downregulated up until 12 h of ABA treatment. Some trihelix genes showed high transcript levels under multiple treatments. For example, OsMSL39 simultaneously responded to all treatments. OsMSL28 was significantly induced by three tested treatments except ABA treatment. In conclusion, all 12 genes are induced by various abiotic stress or stress signaling molecules, but the expression levels are different and there was no significant correlation between gene expression and its classification.
To determine whether light influences trihelix genes expression, the expression levels of 12 OsMSLs from each treatment group under light and dark conditions were investigated (Figure 10 and Figure 11). Overall, some OsMSLs were induced under normal conditions. For example, at 0 h, OsMSL01, OsMSL25b, OsMSL33, OsMSL39, and OsMSL41 expression levels significantly differed between the two conditions. Several genes began to be light-regulated expressed under different treatments. For instance, after ABA exposure, the expression levels of OsMSL16, OsMSL25a, OsMSL25b, OsMSL28, OsMSL34a, and OsMSL39 at different time points under normal conditions were significantly higher than those at the same intervals under dark conditions. In contrast, some genes were repressed by light at different points, such as the expression levels of OsMSL41 at 2 h under ABA treatment, OsMSL28 at 6 h under hydrogen peroxide treatment, OsMSL34b at 12 h under drought stress and OsMSL16 at 12 h under high salty stress were higher in the dark than under normal conditions. Therefore, the relationship between trihelix genes and light is complex and trihelix gene expression is not directly regulated by light but may be controlled by multiple regulatory mechanisms.

3. Discussion

The trihelix DNA binding domain is usually associated with that of Myb/SANT LIKE. There are three strongly conserved, regularly spaced tryptophan (W) residues in each repeating Myb α-helix. The residues of the Myb α-helix regions are also strongly conserved between the GT factors and the Myb/SANT-LIKE proteins. Individual helices are longer, so the trihelix domains of the GT factors and the Myb/SANT LIKE proteins are related [20]. In contrast, other amino acids at this location have longer Myb repeat sequences than the helix-turn-helix structure formed by the Myb-type DNA-binding domain. Therefore, the trihelix family genes have functions and target gene sequences differing from those of the MYB transcription factor family [20]. In the present study, 41 rice trihelix genes were identified using the Myb-type (Myb/SANT-LIKE) DNA-binding domain. However, in a previous study, only 31 rice trihelix genes were identified [21]. Because the repeated search method was performed in the present study, the trihelix gene could be identified in the rice genome more comprehensively. Based on the wide ranges of OsMSL protein MW, isoelectric point, and subcellular localization, we speculated that OsMSLs are not conservatively evolved. Besides, in Arabidopsis, the chimeric trihelix gene At4g17060 was misannotated in the genome [22]. Its ortholog LOC_Os10g41460 (OsMSL36), was identified in this study.
In general, gene families expand by tandem and segmental duplications [23]. Evolutionary conservatism increased with the number of duplicated genes in a gene family [16]. Chromosomal distribution and gene duplication analyses indicated that there are six pairs of duplicated genes in rice among total 41 trihelix gene, including two pairs of tandem genes and four pairs of segmental genes. However, in Arabidopsis, fifteen pairs of duplicated trihelix genes were previously detected among 34 trihelix genes [24]. These results suggested that OsMSLs are less conserved, and most genes may not originate from the same ancestor. On the other hand, these results demonstrate that the rice trihelix family has a high degree of evolutionary divergence and is non-conservative. These properties may account for the substantial differences among the rice trihelix proteins.
We conducted a phylogenetic analysis to elucidate the evolutionary relationships within the rice and other species trihelix gene family. In a previous study, no GTδ was designated for rice or Arabidopsis [18] because the Arabidopsis family Myb/SANT-LIKE DNA-binding domain or protein sequence was aligned with the rice genome in the attempt to identify the potential rice trihelix gene. However, there is a relatively long evolutionary distance between rice and Arabidopsis. Therefore, this method may overlook the specific trihelix genes in rice. Misclassification may have resulted in deviations because of the influence of the Arabidopsis trihelix genes. In the present study, a class of rice trihelix genes and some tomato trihelix genes previously assigned to this subfamily have been found in the GTδ clade. Subsequent investigation revealed that this subfamily has a high structure similarity. For this reason, the evolutionary relationships of its members are more conservative than those in other subgroups. In previous studies, the GT clade was divided into the GT1 and GT2 subfamilies cited [13,14,24]. According to our evolutionary analysis, there are two GT subfamily clusters (Figure 3A). The difference between GT1 and GT2 is smaller than those among SIP1, SH4, GTγ, and GTδ. The genes in GT1 each have one trihelix DNA binding domain with three conserved tryptophans. Those in GT2 each have an additional trihelix DNA binding domain with two tryptophans and one phenylalanine [25]. Consequently, we classified both the GT1 and GT2 clades in the GT subfamily.
Although the evolution of the trihelix family was not conservative, our gene structure and conservative functional domain analyses indicated that the genes within the same subfamilies (especially SIP1 and GTδ) were still relatively conserved. In fact, most duplicated genes occurred in GTδ. The MEME analysis revealed that the functional domain distribution of each OsMSL was related to its classification (Figure 4). Therefore, these conserved functional domains play central roles in group-specific functions. In contrast, the gene structures among the various groups differed greatly from the conserved functional domains. For this reason, they may have different downstream regulatory genes and participate in different signaling pathways.
The distribution and type of cis-elements on the gene promoter may determine OsMSL functions. In this study, we identified five cis-elements shared by all genes out of a large number among 41 OsMSLs. The results showed that OsMSLs were mainly involved in abiotic stress and light-induced responses. To date, little functional analysis has been conducted on the trihelix transcription factors in plants. Previous studies showed that they participate in responses to pathogens and abiotic stress, light induction, and nitrogen metabolism [18]. The light-induced process is a major feature of the trihelix genes. Light induces massive reprogramming of the plant transcriptome and upregulates or downregulates gene expression and its corresponding signaling pathway [26]. Light signaling coordinates the induction or repression of specific downstream genes like bHLH [27], bZIP [28], R2R3-MYB [29], FAR1 [30], and FHY3 [31]. Studies on light regulatory mechanisms in plants focused on the long-term effects of light exposure. However, little attention has been paid to the transient light-responsive processes of transcription factors in plant stress reactions. In the present study, the OsMSL expression profiles disclosed that their responses to light are transient and change with the processing time. Therefore, OsMSLs may be regulated by light in response to abiotic stress in the same way that phototropism, chloroplast movement, and stomatal opening participate in rapid light-responsive processes and are not under extensive transcriptional regulation. This mechanism substantially differs from that observed in relation to gene expression changes in response to the long-term effects of light on photoperiod.
Gene expression specificity in plant tissues and developmental stages may indicate possible gene functions. Previous studies showed that certain trihelix genes in tomato and chrysanthemum exhibited stable expression levels in all tissues [8,14]. However, most OsMSLs do not maintain stable expression levels in different tissues. OsMSL expression levels may vary substantially among tissues. For example, the expression levels of OsMSL16, OsMSL27, OsMSL28, OsMSL35, and OsMSL39 were extremely high in the leaves and sheaths but comparatively low in the roots and stems. Subcellular localization revealed that these genes are expressed in the chloroplasts, and these are present only in the rice leaves and sheaths. The expression patterns of OsMSL16, OsMSL27, OsMSL28, OsMSL35, and OsMSL39 in GTδ were similar. Therefore, the different OsMSLs within the same group explain the parallel functions of the rice trihelix family. OsMSL02, OsMSL08, and OsMSL14 in Group II were expressed at low levels in all four tissues. Either they are inducible or they are only upregulated under special conditions [32].
High OsMSL expression was observed mainly in 7–40 days seedlings (from germination to early tillering). Therefore, they indicate that OsMSLs contribute primarily to the early stages of rice growth, as do many other genes. For example, ZFP182 overexpression enhanced salt, drought, and cold tolerance in transgenic rice seedlings [33]. Loss of the ABA transporter OsPM1 in 35 days rice seedlings conferred greater drought sensitivity than that seen in the WT [34]. The expression patterns of OsMSL34a and OsMSL34b disclosed that they were, in fact, different variable splicing forms of the same gene with very different expression levels. Since OsMSL34b might be the primary variable splicing form of the gene, OsMSL34a was downregulated to some extent.
In the previous study on rice trihelix genes, their responses to various plant hormones were highlighted [21], and the present study focused on abiotic stress and stress signaling molecules. Other than plant growth and development, MSLs participate in stress responses [18]. Although no cis-elements related to the ABA signaling pathway were found in the OsMSL promoter region, ABA nonetheless, induced the rice trihelix genes. ABA accumulates when plants are subjected to a water deficit. It regulates the expression of drought stress-related genes and modulates the molecular, cellular, and physiological mechanisms for adaptation to environmental stress [35]. In chrysanthemum, however, ABA downregulated the trihelix genes but others were upregulated after prolonged ABA exposure [14]. In contrast, GmGT-2A and GmGT-2B in soybean were upregulated by ABA [11]. In the present study, twelve OsMSLs from all subfamilies have been induced by ABA. These results suggest that the signaling mechanisms of the trihelix family genes vary with species. Whether the ABA presence has a negative regulation on the trihelix gene requires further experimental verification.
ROS (reactive oxygen species) are produced in response to most environmental stress. Excessive ROS accumulation may irreversibly damage cells [36,37]. In previous studies, trihelix family genes had at least one response during plant osmotic stress defense [18]. It follows that trihelix family genes may also participate in ROS scavenging and enhance plant tolerance to various stresses. However, little is known about the mechanism of the peroxide reaction mediated by the trihelix family. We performed a quantitative PCR analysis on 12 OsMSLs subjected to hydrogen peroxide. All 12 OsMSLs responded to hydrogen peroxide stress. Therefore, they may help improve the permeation tolerance by increasing the ROS scavenging capacity in rice.
The trihelix transcription factors bind to GT elements on the light-regulating genes [25]. In darkness, photo regulatory genes are repressed and their associated trihelix family genes are also affected. Nevertheless, certain trihelix family genes are downregulated in response to light exposure, apparently because they must be repressed to be able to downregulate target whose expression is light-dependent. Certain constitutively expressed trihelix genes occur in Arabidopsis. Their expression is ubiquitous and indifferent to the light regime. They are concentrated mainly in GT1 and GT2 [38]. In the present study, however, the expression of the 12 rice trihelix genes changed direction at least once after light or dark treatment. Therefore, they may be inducible rather than constitutive. It remains to be determined whether OsMSLs are regulated by one wavelength or a wide light spectrum. Further investigation of the light path and its components is necessary. The results of our study have helped initiate the research of the rice trihelix transcription factors. In future research, the relationships among the OsMSLs, ABA-mediated dehydration stress tolerance, and ROS scavenging ability under light regulation should be explored.

4. Materials and Methods

4.1. Identification and Sequence Analysis of Trihelix Transcription Factor Family in Rice

The rice trihelix transcription factors were identified according to a previously described method with minor changes [39]. The Hidden Markov Model of Myb/SANT-LIKE domain (PF13837) was downloaded from the Pfam database (http://pfam.xfam.org/) [40]. The entire rice amino acid, genome, and CDS sequence assembly and corresponding annotation were downloaded from the EnsemblPlants database (http://plants.ensembl.org/index.html) [41]. The candidate proteins were sought by the HMMSEARCH program (https://www.ebi.ac.uk/Tools/hmmer/search/hmmsearch) base on the Bio-Linux system (Dr Tracey Timms-Wilson, Centre for Ecology & Hydrology (CEH), Oxfordshire, UK). The domain sequences of these candidate proteins were extracted and used to build a rice-specific Hidden Markov Model. All rice proteins were detected by the rice-specific Hidden Markov Model. Those with E-value < 0.01 were selected. The trihelix proteins were verified using the Pfam and InterPro databases (http://www.ebi.ac.uk/interpro/) [42]. Proteins obtained by the domain and database screening confirmation were considered trihelix family members. The corresponding CDS and gene sequences were extracted according to their protein identifications.
The MEME program (http://meme-suite.org/) identified conserved motifs of the trihelix family proteins with the following parameters: Any number of repetitions; minimum seven motifs; maximum 49 motifs; optimum 10–200 amino acids; expected E-value < 1 × 10−48. The trihelix family gene structures were displayed by comparing the coding and genomic sequences with the Gene Structure Display Server tools (http://gsds.cbi.pku.edu.cn/) [43]. The chromosomal locations of the trihelix family genes were mapped onto the rice linkage map with an online tool according to their TIGR numbers [44]. The isoelectric points and molecular weights of trihelix family proteins were estimated with ExPASy (http://expasy.org/) [45].

4.2. Phylogenetic Analysis

A multiple alignment was performed with the full-length amino acid sequences of the rice trihelix family proteins using MEGA v. 7.0 (https://www.megasoftware.net/) [46]. Unrooted trees were constructed by the maximum-likelihood (ML) method with the following parameters: Poisson correction; pairwise deletion; 1000 bootstrap replicates.

4.3. Gene Duplication and Ka/Ks Analysis

Synteny blocks of the rice genome were downloaded from the Plant Genome Duplication Database (PGDD, http://chibba.agtec.uga.edu/duplication/) [47]. Duplicated OsMSLs pairs were connected by solid lines.

4.4. Cis-Element Analysis of Trihelix Transcription Factor Family

Promoters of the trihelix family genes were downloaded from the Phytozome database (https://phytozome.jgi.doe.gov/pz/portal.html#) [48]. The PLACE database (https://sogo.dna.affrc.go.jp/) was used to analyze the cis-regulatory elements of the trihelix family gene promoters [49].

4.5. Plant Growth Conditions and Treatments

Nipponbare rice seeds (O. sativa L. ssp. japonica) were surface-sterilized with 10% sodium hypochlorite solution for 30 min then sown on 1/2 MS (Murashige & Skoog) solid medium and cultured in a light incubator. After 2 weeks, the seedlings were at the two true leaf stage. They were transplanted into Hoagland’s nutrient solution and cultured in an artificial climate chamber under controlled conditions (14 h light at 28 °C/10 h dark at 22 °C; relative humidity 70%). The rice seedlings were subjected to various stresses at the three-leaf stage (4 weeks) [50].
For the drought, salt, and hydrogen peroxide stress treatments, the rice seedlings were transferred to Hoagland’s nutrient solution containing 20% polyethylene glycol (PEG)-6000 (w/v), 150 mM NaCl, or 2% hydrogen peroxide (v/v), respectively. For the ABA treatment, the rice seedlings were cultured on 1/2 MS solid medium containing 10 μM ABA. The control group was maintained on normal nutrient solution or medium. All other culture conditions were the same as described above. Treated rice tissues were harvested at 0 h, 1 h, 2 h, 4 h, 6 h, 12 h, 24 h, and 48 h. The samples were immediately placed in liquid nitrogen and stored at −80 °C until use. Untreated material was used as a control. The experimental procedure was repeated at least three times.

4.6. Expression Analysis of Trihelix Transcription Factor Family

Total RNA was extracted from rice tissues by the TRIzol method (Thermo Fisher Scientific, Waltham, MA, USA) and treated with DNase to eliminate any DNA contamination. RNA quality was assessed by electrophoresis and stored at −80 °C until use. First-strand cDNA (10 μL) was synthesized according to the instructions for the PrimeScript™ RT Master Mix (Takara Biomedical Technology (Beijing) Co., Ltd., Beijing, China). Primers were designed with Primer Premier v. 5.0 (PREMIER Biosoft International, Palo Alto, CA, USA) and were based on the trihelix gene family transcript sequences. Gene specific primers for quantitative real-time PCR are listed in Table S6. Primer amplification specificity was verified in the rice genome database using Blast from NCBI (https://www.ncbi.nlm.nih.gov/) [51]. Rice β-actin was the internal reference gene. Quantitative real-time PCR was performed in the ABI 7300 Real Time PCR System (Applied Biosystems, Foster City, CA, USA) using SYBR Green chemistry and reaction mix consists of 10 μL SYBR qPCR Master Mix (Vazyme Biotech Co.,Ltd., Nanjing, China), 0.4 μL upstream and downstream primers respectively, 0.4 μL ROX, 2 μL cDNA (10 times dilution) and 6.8 μL ddH2O to 20 μL. The PCR reaction protocol was 95 °C for 5 min; 95 °C for 10 s; 60 °C for 20 s; 72 °C for 20 s; 45 cycles. Gene expression levels were calculated by the 2-ΔΔCT method: ΔΔCT = (CTtarget − CTactin) at time x − (CTtarget − CTactin) at time 0 [52]. The test was repeated three times. Expression data for the rice trihelix family genes were retrieved from the Expression Atlas database (https://www.ebi.ac.uk/gxa/home) [53]. Heatmaps were created in HemI v.1.0 (The CUCKOO Workgroup, Hubei, China) and based on the expression data [54].

5. Conclusions

Trihelix transcription factors participate in many plant biological processes but have not been systematically studied. Here, 41 rice trihelix transcription factors were identified by bioinformatics analysis. Gene synteny analysis showed that OsMSLs are less conserved, and most genes may not originate from the same ancestor. Phylogenetic analysis categorized them into five subfamilies. The gene structures and conserved functional domains of the rice trihelix transcription factors varied greatly but they shared five cis-elements governing light dehydration stress responses. Expression pattern analysis revealed that the trihelix transcription factors had the highest expression levels in the early rice growth stages and most of the strongly upregulated genes were localized in the leaves and sheaths. Real-time quantitative PCR analysis of the trihelix family genes subjected to various stressors or ABA revealed that they were induced in response to drought, high salt, hydrogen peroxide, and ABA. However, these responses were also regulated by light and individual OsMSL expressions differed with the presence or absence of light. Our study helped elucidate the biological functions of the trihelix transcription factors in rice.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/20/2/251/s1.

Author Contributions

Conceptualization, J.L.; data curation, X.M.; formal analysis, J.S.; funding acquisition, D.Z.; investigation, J.W. (Jing Wang) and J.W. (Jingguo Wang); project administration, H.Z. (Hongwei Zhao); resources, H.Z. (Hualong Liu); software, H.Z. (Hongliang Zheng); supervision, Z.Z.; writing—original draft, J.L.; Writing—review & editing, M.Z. All authors read and approved the final manuscript.

Funding

This research was funded by 1. National Science and Technology Major Project (2018ZX0800912B-002). 2. Youth Science Foundation of Heilongjiang Province (QC2017015). 3. National Natural Science Foundation (31701507). 4. Natural Science Foundation of Heilongjiang Province of China (No. C2015003).

Acknowledgments

We thank Editage for its linguistic assistance during the preparation of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lindemose, S.; O’Shea, C.; Jensen, M.K.; Skriver, K. Structure, Function and Networks of Transcription Factors Involved in Abiotic Stress Responses. Int. J. Mol. Sci. 2013, 14, 5842–5878. [Google Scholar] [CrossRef] [PubMed]
  2. Jin, J.P.; Tian, F.; Yang, D.C.; Meng, Y.Q.; Kong, L.; Luo, J.C.; Gao, G. PlantTFDB 4.0: Toward a central hub for transcription factors and regulatory interactions in plants. Nucleic. Acids Res. 2017, 45, D1040–D1045. [Google Scholar] [CrossRef] [PubMed]
  3. Riechmann, J.L.; Heard, J.; Martin, G.; Reuber, L.; Jiang, C.; Keddie, J.; Adam, L.; Pineda, O.; Ratcliffe, O.J.; Samaha, R.R. Arabidopsis transcription factors: Genome-wide comparative analysis among eukaryotes. Science 2000, 290, 2105. [Google Scholar] [CrossRef] [PubMed]
  4. Green, P.J.; Kay, S.A.; Chua, N.H. Sequence-specific interactions of a pea nuclear factor with light-responsive elements upstream of the rbcS-3A gene. EMBO J. 1987, 6, 2543–2549. [Google Scholar] [CrossRef] [PubMed]
  5. Nagano, Y. Several Features of the GT-Factor Trihelix Domain Resemble Those of the Myb DNA-Binding Domain. Plant Physiol. 2000, 124, 491–493. [Google Scholar] [CrossRef] [PubMed]
  6. Qin, Y.; Ma, X.; Yu, G.H.; Wang, Q.; Wang, L.; Kong, L.R.; Kim, W.; Wang, H.W. Evolutionary History of Trihelix Family and Their Functional Diversification. DNA Res. 2014, 21, 499–510. [Google Scholar] [CrossRef] [Green Version]
  7. Gao, M.J.; Lydiate, D.J.; Li, X.; Lui, H.; Gjetvaj, B.; Hegedus, D.D.; Rozwadowski, K. Repression of Seed Maturation Genes by a Trihelix Transcriptional Repressor in Arabidopsis Seedlings. Plant Cell 2009, 21, 54–71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Yu, C.Y.; Cal, X.F.; Ye, Z.B.; Li, H.X. Genome-wide identification and expression profiling analysis of trihelix gene family in tomato. Biochem. Biophys. Res. Commun. 2015, 468, 653–659. [Google Scholar] [CrossRef] [PubMed]
  9. Murata, J.; Takase, H.; Hiratsuka, K. Characterization of a Novel GT-box Binding Protein from Arabidopsis. Plant Biotechnol. 2002, 19, 103–112. [Google Scholar] [CrossRef]
  10. Wang, R.H.G.; Han, B. Transcript abundance of rml1, encoding a putative GT1-like factor in rice, is up-regulated by Magnaporthe grisea and down-regulated by light. Gene 2004, 324, 105–115. [Google Scholar] [CrossRef] [PubMed]
  11. Xie, Z.M.; Zou, H.F.; Lei, G.; Wei, W.; Zhou, Q.Y.; Niu, C.F.; Liao, Y.; Tian, A.G.; Ma, B.; Zhang, W.K. Soybean Trihelix Transcription Factors GmGT-2A and GmGT-2B Improve Plant Tolerance to Abiotic Stresses in Transgenic Arabidopsis. PLoS ONE 2009, 4, e6898. [Google Scholar] [CrossRef]
  12. Chan, Y.Y.; Pence, H.E.; Jing, B.J.; Miura, K.; Gosney, M.J.; Hasegawa, P.M.; Mickelbart, M.V. The Arabidopsis GTL1 Transcription Factor Regulates Water Use Efficiency and Drought Tolerance by Modulating Stomatal Density via Transrepression of SDD1. Plant Cell 2010, 22, 4128. [Google Scholar]
  13. Fang, Y.; Xie, K.; Xin, H.; Hu, H.; Xiong, L. Systematic analysis of GT factor family of rice reveals a novel subfamily involved in stress responses. Mol. Genet. Genom. 2010, 283, 157–169. [Google Scholar] [CrossRef] [PubMed]
  14. Song, A.P.; Wu, D.; Fan, Q.Q.; Tian, C.; Chen, S.M.; Guan, Z.Y.; Xin, J.J.; Zhao, K.K.; Chen, F.D. Transcriptome-Wide Identification and Expression Profiling Analysis of Chrysanthemum Trihelix Transcription Factors. Int. J. Mol. Sci. 2016, 17, 198. [Google Scholar] [CrossRef] [PubMed]
  15. Li, C.; Zhou, A.; Sang, T. Rice Domestication by Reducing Shattering. Science 2006, 311, 1936–1939. [Google Scholar] [CrossRef] [PubMed]
  16. Holub, E.B. The arms race is ancient history in Arabidopsis, the wildflower. Nat. Rev. Genet. 2001, 2, 516–527. [Google Scholar] [CrossRef]
  17. Wang, Y.; Tang, H.; DeBarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.-H.; Jin, H.; Marler, B.; Guo, H. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [PubMed]
  18. KaplanLevy, R.N.; Brewer, P.B.; Quon, T.; Smyth, D.R. The trihelix family of transcription factors—Light, stress and development. Trends Plant Sci. 2012, 17, 163–171. [Google Scholar] [CrossRef] [PubMed]
  19. Markus, S.; Davison, T.S.; Henz, S.R.; Pape, U.J.; Monika, D.; Martin, V.; Bernhard, S.L.; Detlef, W.; Lohmann, J.U. A gene expression map of Arabidopsis thaliana development. Nat. Genet. 2005, 37, 501–506. [Google Scholar]
  20. Nagata, T.; Niyada, E.; Fujimoto, N.; Nagasaki, Y.; Noto, K.; Miyanoiri, Y.; Murata, J.; Hiratsuka, K.; Katahira, M. Solution structures of the trihelix DNA-binding domains of the wild-type and a phosphomimetic mutant of Arabidopsis GT-1: Mechanism for an increase in DNA-binding affinity through phosphorylation. Proteins 2010, 78, 3033–3047. [Google Scholar] [CrossRef] [PubMed]
  21. Jianhui, J.; Yingjun, Z.; Hehe, W.; Liming, Y. Genome-wide analysis and functional prediction of the Trihelix transcription factor family in rice. Hereditas 2015, 37, 1228. (In Chinese) [Google Scholar]
  22. Geraldo, N.; Baurle, I.; Kidou, S.; Hu, X.; Dean, C. FRIGIDA Delays Flowering in Arabidopsis via a Cotranscriptional Mechanism Involving Direct Interaction with the Nuclear Cap-Binding Complex. Plant Physiol. 2009, 150, 1611. [Google Scholar] [CrossRef]
  23. Cannon, S.B.; Mitra, A.; Baumgarten, A.; Young, N.D.; May, G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 2004, 4, 10. [Google Scholar] [CrossRef] [PubMed]
  24. Ali, M.A.; Yasmeen, E.; Riaz, M.; Azeem, F.; Sultan, S.; Abbas, A.; Riaz, K. Genome-wide analysis of trihelix Transcription factor gene family in Arabidopsis thaliana. Pak. J. Agric. Sci. 2016, 53, 439–448. [Google Scholar]
  25. Lam, E. Domain analysis of the plant DNA-binding protein GT1a: Requirement of four putative alpha-helices for DNA binding and identification of a novel oligomerization region. Mol. Cell. Biol. 1995, 15, 1014–1020. [Google Scholar] [CrossRef] [PubMed]
  26. Jiao, Y.; Lau, O.S.; Deng, X.W. Light-regulated transcriptional networks in higher plants. Nat. Rev. Genet. 2007, 8, 217–230. [Google Scholar] [CrossRef] [PubMed]
  27. Toledoortiz, G.; Huq, E.; Quail, P.H. The Arabidopsis Basic/Helix-Loop-Helix Transcription Factor Family. Plant Cell 2003, 15, 1749–1770. [Google Scholar] [CrossRef] [Green Version]
  28. Jakoby, M.; Weisshaar, B.; Dröge-Laser, W.; Vicente-Carbajosa, J.; Tiedemann, J.; Kroj, T.; Parcy, F. bZIP transcription factors in Arabidopsis. Trends Plant Sci. 2002, 7, 106–111. [Google Scholar] [CrossRef]
  29. Ballesteros, M.L.; Bolle, C.; Lois, L.M.; Moore, J.M.; Viellecalzada, J.P.; Grossniklaus, U.; Chua, N.H. LAF1, a MYB transcription activator for phytochrome A signaling. Gene Dev. 2001, 15, 2613–2625. [Google Scholar] [CrossRef] [Green Version]
  30. Hudson, M.; Ringli, C.; Boylan, M.T.; Quail, P.H. The FAR1 locus encodes a novel nuclear protein specific to phytochrome A signaling. Gene Dev. 1999, 13, 2017–2027. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, H.; Deng, X.W. Arabidopsis FHY3 defines a key phytochrome A signaling component directly interacting with its homologous partner FAR1. EMBO J. 2002, 21, 1339–1349. [Google Scholar] [CrossRef]
  32. Tang, Y.; Qin, S.; Guo, Y.; Chen, Y.; Wu, P.; Chen, Y.; Li, M.; Jiang, H.; Wu, G. Genome-Wide Analysis of the AP2/ERF Gene Family in Physic Nut and Overexpression of the JcERF011 Gene in Rice Increased Its Sensitivity to Salinity Stress. PLoS ONE 2016, 11, e0150879. [Google Scholar] [CrossRef] [PubMed]
  33. Huang, J.; Sun, S.; Xu, D.; Lan, H.; Sun, H.; Wang, Z.; Bao, Y.; Wang, J.; Tang, H.; Zhang, H. A TFIIIA-type zinc finger protein confers multiple abiotic stress tolerances in transgenic rice (Oryza sativa L.). Plant Mol. Biol. 2012, 80, 337. [Google Scholar] [CrossRef] [PubMed]
  34. Yao, L.; Cheng, X.; Gu, Z.; Huang, W.; Li, S.; Wang, L.; Wang, Y.F.; Xu, P.; Ma, H.; Ge, X. The AWPM-19 Family Protein OsPM1 Mediates Abscisic Acid Influx and Drought Response in Rice. Plant Cell 2018, 30, 1258–1276. [Google Scholar] [CrossRef]
  35. Yamaguchi-Shinozaki, K.; Shinozaki, K. Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Annu. Rev. Plant Biol. 2006, 57, 781–803. [Google Scholar] [CrossRef] [PubMed]
  36. Apel, K.; Hirt, H. Reactive oxygen species: Metabolism, oxidative stress, and signal transduction. Annu. Rev. Plant Biol. 2004, 55, 373–399. [Google Scholar] [CrossRef] [PubMed]
  37. Miller, G.; Suzuki, N.; Ciftci-Yilmaz, S.; Mittler, R. Reactive oxygen species homeostasis and signalling during drought and salinity stresses. Plant Cell Environ. 2010, 33, 453–467. [Google Scholar] [CrossRef] [PubMed]
  38. Gilmartin, P.M.; Memelink, J.; Hiratsuka, K.; Kay, S.A.; Chua, N.H. Characterization of a gene encoding a DNA binding protein with specificity for a light-responsive element. Plant Cell 1992, 4, 839–849. [Google Scholar] [CrossRef]
  39. Lozano, R.; Hamblin, M.T.; Prochnik, S.; Jannink, J.L. Identification and distribution of the NBS-LRR gene family in the Cassava genome. BMC Genom. 2015, 16, 360. [Google Scholar] [CrossRef] [PubMed]
  40. El-Gebali, S.; Mistry, J.; Bateman, A.; Eddy, S.R.; Luciani, A.; Potter, S.C.; Qureshi, M.; Richardson, L.J.; Salazar, G.A.; Smart, A. The Pfam protein families database in 2019. Nucleic Acids Res. 2018, 28, 263–266. [Google Scholar] [CrossRef]
  41. Kersey, P.J.; Allen, J.E.; Allot, A.; Barba, M.; Boddu, S.; Bolt, B.J.; Carvalho-Silva, D.; Christensen, M.; Davis, P.; Grabmueller, C. Ensembl Genomes 2018: An integrated omics infrastructure for non-vertebrate species. Nucleic Acids Res. 2018, 46, D802–D808. [Google Scholar] [CrossRef] [PubMed]
  42. Mitchell, A.L.; Attwood, T.K.; Babbitt, P.C.; Blum, M.; Bork, P.; Bridge, A.; Brown, S.D.; Chang, H.-Y.; El-Gebali, S.; Fraser, M.I. InterPro in 2019: Improving coverage, classification and access to protein sequence annotations. Nucleic Acids Res. 2018. [Google Scholar] [CrossRef] [PubMed]
  43. Hu, B.; Jin, J.; Guo, A.Y.; Zhang, H.; Luo, J.; Gao, G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 2014, 31, 1296. [Google Scholar] [CrossRef] [PubMed]
  44. Kurata, N.; Yamazaki, Y. Oryzabase. An integrated biological and genome information database for rice. Plant Physiol. 2006, 140, 12–17. [Google Scholar] [CrossRef] [PubMed]
  45. Artimo, P.; Jonnalagedda, M.; Arnold, K.; Baratin, D.; Csardi, G.; De Castro, E.; Duvaud, S.; Flegel, V.; Fortier, A.; Gasteiger, E. ExPASy: SIB bioinformatics resource portal. Nucleic Acids Res. 2012, 40, W597–W603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  47. Lee, T.H.; Tang, H.; Wang, X.; Paterson, A.H. PGDD: A database of gene and genome duplication in plants. Nucleic Acids Res. 2013, 41, 1152–1158. [Google Scholar] [CrossRef] [PubMed]
  48. Goodstein, D.M.; Shu, S.; Howson, R.; Neupane, R.; Hayes, R.D.; Fazo, J.; Mitros, T.; Dirks, W.; Hellsten, U.; Putnam, N. Phytozome: A comparative platform for green plant genomics. Nucleic Acids Res. 2011, 40, D1178–D1186. [Google Scholar] [CrossRef]
  49. Higo, K.; Ugawa, Y.; Iwamoto, M.; Korenaga, T. Plant cis-acting regulatory DNA elements (PLACE) database: 1999. Nucleic Acids Res. 1999, 27, 297–300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Ren, J.; Gao, F.; Wu, X.; Lu, X.; Zeng, L.; Lv, J.; Su, X.; Luo, H.; Ren, G. Bph32, a novel gene encoding an unknown SCR domain-containing protein, confers resistance against the brown planthopper in rice. Sci. Rep. 2016, 6, 37645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Coordinators, N.R. Database resources of the national center for biotechnology information. Nucleic Acids Res. 2017, 45, D12. [Google Scholar]
  52. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  53. Papatheodorou, I.; Fonseca, N.A.; Keays, M.; Tang, Y.A.; Barrera, E.; Bazant, W.; Burke, M.; Füllgrabe, A.; Fuentes, A.M.-P.; George, N. Expression Atlas: Gene and protein expression across multiple studies and organisms. Nucleic Acids Res. 2017, 46, D246–D251. [Google Scholar] [CrossRef] [PubMed]
  54. Deng, W.; Wang, Y.; Liu, Z.; Cheng, H.; Xue, Y. HemI: A Toolkit for Illustrating Heatmaps. PLoS ONE 2014, 9, e111988. [Google Scholar] [CrossRef] [PubMed]
Ijms 20 00251 g001 550
Figure 1. Chromosomal locations of rice trihelix genes. Black bars represent the chromosomes. Chromosome numbers are shown at the tops of the bar. Trihelix genes are labeled at the right of the chromosomes. Scale bar on the left indicates the chromosome lengths (Mb).
Figure 1. Chromosomal locations of rice trihelix genes. Black bars represent the chromosomes. Chromosome numbers are shown at the tops of the bar. Trihelix genes are labeled at the right of the chromosomes. Scale bar on the left indicates the chromosome lengths (Mb).
Ijms 20 00251 g001
Ijms 20 00251 g002 550
Figure 2. Schematic representations of segmental duplications of rice trihelix genes. Different color lines indicate all synteny blocks in rice genome between each chromosome, and the thick red lines indicate duplicated trihelix gene pairs. The chromosome number is indicated at the bottom of each chromosome. Scale bar marked on the chromosome indicating chromosome lengths (Mb).
Figure 2. Schematic representations of segmental duplications of rice trihelix genes. Different color lines indicate all synteny blocks in rice genome between each chromosome, and the thick red lines indicate duplicated trihelix gene pairs. The chromosome number is indicated at the bottom of each chromosome. Scale bar marked on the chromosome indicating chromosome lengths (Mb).
Ijms 20 00251 g002
Ijms 20 00251 g003 550
Figure 3. Synteny analysis of trihelix genes between rice and (A) dicotyledonous plant Arabidopsi thaliana, (B) monocotyledonous plant Brachypodium distachyon, wheat and maize. Gray lines in the background indicate the collinear blocks within rice and other plant genomes, while the red lines highlight the syntenic trihelix gene pairs. The species names with the prefixes ‘A. thaliana’, ‘B. distachyon’, ‘T. aestivum’, ‘Z. mays’ and ‘O. sativa’ indicate Arabidopsi thaliana, Brachypodium distachyon, Triticum aestivum, Zea mays and Oryza sativa, respectively. Red or green bars represent the chromosomes. The chromosome number is labeled at the top or bottom of each chromosome.
Figure 3. Synteny analysis of trihelix genes between rice and (A) dicotyledonous plant Arabidopsi thaliana, (B) monocotyledonous plant Brachypodium distachyon, wheat and maize. Gray lines in the background indicate the collinear blocks within rice and other plant genomes, while the red lines highlight the syntenic trihelix gene pairs. The species names with the prefixes ‘A. thaliana’, ‘B. distachyon’, ‘T. aestivum’, ‘Z. mays’ and ‘O. sativa’ indicate Arabidopsi thaliana, Brachypodium distachyon, Triticum aestivum, Zea mays and Oryza sativa, respectively. Red or green bars represent the chromosomes. The chromosome number is labeled at the top or bottom of each chromosome.
Ijms 20 00251 g003
Ijms 20 00251 g004 550
Figure 4. Phylogenetic relationships among 105 trihelix proteins in rice, Arabidopsis, soybean, maize, tomato, wheat, chrysanthemum, wild rice and Brachypodium distachyon. The maximum likelihood tree was created using MEGA v. 7.0 (bootstrap value = 1000) and the bootstrap value of each branch is displayed. Forty-three OsMSL proteins are marked with black circles and other species are marked with white circles. The phylogenetic tree was clustered into SIP1, GTγ, GT, SH4, and GTδ.
Figure 4. Phylogenetic relationships among 105 trihelix proteins in rice, Arabidopsis, soybean, maize, tomato, wheat, chrysanthemum, wild rice and Brachypodium distachyon. The maximum likelihood tree was created using MEGA v. 7.0 (bootstrap value = 1000) and the bootstrap value of each branch is displayed. Forty-three OsMSL proteins are marked with black circles and other species are marked with white circles. The phylogenetic tree was clustered into SIP1, GTγ, GT, SH4, and GTδ.
Ijms 20 00251 g004
Ijms 20 00251 g005 550
Figure 5. Phylogenetic analysis and gene structure of the rice trihelix family. (A) Phylogenetic analysis of the rice trihelix family. The phylogenetic tree was constructed based on the full-length amino acid sequences of the rice trihelix proteins by using MEGA v. 7.0 with the maximum-likelihood method. Bootstrap = 1,000. SIP1, GTγ, GT, SH4, and GTδ are marked with different colors. (B) Gene structures of the rice trihelix family. These were analyzed by the Gene Structure Display Server (GSDS v. 2.0). Exons, introns, and untranslated regions are marked by round red rectangles, black lines, and blue rectangles, respectively. The scale bar at the bottom estimates the lengths of the exons, introns, and untranslated regions.
Figure 5. Phylogenetic analysis and gene structure of the rice trihelix family. (A) Phylogenetic analysis of the rice trihelix family. The phylogenetic tree was constructed based on the full-length amino acid sequences of the rice trihelix proteins by using MEGA v. 7.0 with the maximum-likelihood method. Bootstrap = 1,000. SIP1, GTγ, GT, SH4, and GTδ are marked with different colors. (B) Gene structures of the rice trihelix family. These were analyzed by the Gene Structure Display Server (GSDS v. 2.0). Exons, introns, and untranslated regions are marked by round red rectangles, black lines, and blue rectangles, respectively. The scale bar at the bottom estimates the lengths of the exons, introns, and untranslated regions.
Ijms 20 00251 g005
Ijms 20 00251 g006 550
Figure 6. Motif composition of rice trihelix proteins. Motif analysis was performed using the MEME program as described in the methods section. The trihelix proteins are listed on the left. Boxes of different colors represent the various motifs. Their location in each sequence is marked. Motif sequences are shown in Figure S2. The scale bar at the bottom indicates the lengths of the trihelix protein sequences.
Figure 6. Motif composition of rice trihelix proteins. Motif analysis was performed using the MEME program as described in the methods section. The trihelix proteins are listed on the left. Boxes of different colors represent the various motifs. Their location in each sequence is marked. Motif sequences are shown in Figure S2. The scale bar at the bottom indicates the lengths of the trihelix protein sequences.
Ijms 20 00251 g006
Ijms 20 00251 g007 550
Figure 7. Predicted cis-elements in the promoter regions of the rice trihelix genes. All promoter sequences (−1500 bp) were analyzed. The trihelix genes are shown on the left side of the figure. The scale bar at the bottom indicates the length of promoter sequence. Green bar (GB): GATABOX element; purple bar (A): ACGTATERD1 element; red bar (GC): GT1CONSENSUS element; gray bar (I): INRNTPSADB element; blue bar (M): GT1GMSCAM4 element.
Figure 7. Predicted cis-elements in the promoter regions of the rice trihelix genes. All promoter sequences (−1500 bp) were analyzed. The trihelix genes are shown on the left side of the figure. The scale bar at the bottom indicates the length of promoter sequence. Green bar (GB): GATABOX element; purple bar (A): ACGTATERD1 element; red bar (GC): GT1CONSENSUS element; gray bar (I): INRNTPSADB element; blue bar (M): GT1GMSCAM4 element.
Ijms 20 00251 g007
Ijms 20 00251 g008 550
Figure 8. Expression of the rice trihelix gene family in various tissues (R: Root; St: Stem; L: Leaf; Sh: Sheath). (A) Numbers of expressed genes in each tissue. Expression data of the rice trihelix gene family were retrieved from the Expression Atlas database. Extremely high: Expression value > 6, high: 6 ≥ expression value > 4, medium: 4 ≥ expression value > 2, low: 2 ≥ expression value > 0; (B) Expression patterns of the trihelix genes in various rice tissues. Heatmaps were created in HemI v.1.0 and based on the expression data. Expression levels are depicted by different colors on the scale. Green and red represent low and high expression levels, respectively.
Figure 8. Expression of the rice trihelix gene family in various tissues (R: Root; St: Stem; L: Leaf; Sh: Sheath). (A) Numbers of expressed genes in each tissue. Expression data of the rice trihelix gene family were retrieved from the Expression Atlas database. Extremely high: Expression value > 6, high: 6 ≥ expression value > 4, medium: 4 ≥ expression value > 2, low: 2 ≥ expression value > 0; (B) Expression patterns of the trihelix genes in various rice tissues. Heatmaps were created in HemI v.1.0 and based on the expression data. Expression levels are depicted by different colors on the scale. Green and red represent low and high expression levels, respectively.
Ijms 20 00251 g008
Ijms 20 00251 g009 550
Figure 9. Expression of the rice trihelix gene family at different developmental stages (S1: 7 days after sowing; S2: 20 days after sowing; S3: 40 days after sowing; S4: 80 days after sowing; S5: 100 days after sowing; S6: 140 days after sowing). (A) Numbers of expressed genes in various developmental stages. Expression data of the rice trihelix gene family genes were retrieved from the Expression Atlas database. Extremely high: Expression value > 6, high: 6 ≥ expression value > 4, medium: 4 ≥ expression value > 2, low: 2 ≥ expression value > 0; (B) Expression patterns of trihelix genes in various rice developmental stages. Heatmaps were created in HemI v.1.0 and based on the expression data. Expression levels are depicted by different colors on the scale. Blue and red represent low and high expression levels, respectively.
Figure 9. Expression of the rice trihelix gene family at different developmental stages (S1: 7 days after sowing; S2: 20 days after sowing; S3: 40 days after sowing; S4: 80 days after sowing; S5: 100 days after sowing; S6: 140 days after sowing). (A) Numbers of expressed genes in various developmental stages. Expression data of the rice trihelix gene family genes were retrieved from the Expression Atlas database. Extremely high: Expression value > 6, high: 6 ≥ expression value > 4, medium: 4 ≥ expression value > 2, low: 2 ≥ expression value > 0; (B) Expression patterns of trihelix genes in various rice developmental stages. Heatmaps were created in HemI v.1.0 and based on the expression data. Expression levels are depicted by different colors on the scale. Blue and red represent low and high expression levels, respectively.
Ijms 20 00251 g009
Ijms 20 00251 g010 550
Figure 10. Expression profiles of 12 selected OsMSLs in response to (A) ABA (abscisic acid), (B) hydrogen peroxide treatments. Data were normalized to the β-actin gene. Vertical bars indicate standard deviations. Asterisks indicate corresponding genes significantly upregulated or downregulated compared with the control. (* p < 0.05; ** p < 0.01; Student’s t-test).
Figure 10. Expression profiles of 12 selected OsMSLs in response to (A) ABA (abscisic acid), (B) hydrogen peroxide treatments. Data were normalized to the β-actin gene. Vertical bars indicate standard deviations. Asterisks indicate corresponding genes significantly upregulated or downregulated compared with the control. (* p < 0.05; ** p < 0.01; Student’s t-test).
Ijms 20 00251 g010
Ijms 20 00251 g011 550
Figure 11. Expression profiles of 12 selected OsMSLs in response to (A) drought, (B) high salt stress, treatments. Data were normalized to the β-actin gene. Vertical bars indicate standard deviations. Asterisks indicate corresponding genes significantly upregulated or downregulated compared with the control. (* p < 0.05; ** p < 0.01; Student’s t-test).
Figure 11. Expression profiles of 12 selected OsMSLs in response to (A) drought, (B) high salt stress, treatments. Data were normalized to the β-actin gene. Vertical bars indicate standard deviations. Asterisks indicate corresponding genes significantly upregulated or downregulated compared with the control. (* p < 0.05; ** p < 0.01; Student’s t-test).
Ijms 20 00251 g011
Table
Table 1. Detailed information of all trihelix family genes identified in the rice genome.
Table 1. Detailed information of all trihelix family genes identified in the rice genome.
Gene NameGene LocusChrORF (bp)Exon No.Length (aa)MW (kD)pILocalization
OsMSL01LOC_Os01g11200.1Chr1828327531.88.57Nucleus
OsMSL02LOC_Os01g27590.1Chr11131237641.088.6Nucleus
OsMSL03LOC_Os01g34400.1Chr11227340846.069.22Chloroplast
OsMSL04LOC_Os01g36850.1Chr11170238957.619.46Peroxisome
OsMSL05LOC_Os01g48320.1Chr1999133235.929.62Nucleus
OsMSL06LOC_Os01g52090.1Chr1969132235.765.54Nucleus
OsMSL07LOC_Os01g62410.1Chr11764758764.18.5Nucleus
OsMSL08LOC_Os02g01380.1Chr21113137040.55.33Nucleus
OsMSL09LOC_Os02g07800.1Chr21308243546.244.45Chloroplast
OsMSL10LOC_Os02g08450.1Chr21968465574.357.09Chloroplast
OsMSL11LOC_Os02g31160.1Chr21125237439.555.17Nucleus
OsMSL12LOC_Os02g33610.1Chr226491888297.378.97Chloroplast
OsMSL13LOC_Os02g33770.1Chr21233141046.676.28Nucleus
OsMSL14LOC_Os02g35690.1Chr21260141944.446.11Nucleus
OsMSL15LOC_Os02g43300.1Chr21887362867.744.87Nucleus
OsMSL16LOC_Os03g44130.1Chr31104336741.746.61Chloroplast
OsMSL17LOC_Os03g46350.1Chr31038234542.4911.38Nucleus
OsMSL18LOC_Os04g21860.1Chr41254141747.159.1Nucleus
OsMSL19LOC_Os04g30890.1Chr4801126628.629.64Nucleus
OsMSL20LOC_Os04g36790.1Chr41254141743.936.53Nucleus
OsMSL21LOC_Os04g40930.1Chr41158538541.935.82Nucleus
OsMSL22LOC_Os04g45750.1Chr41587252857.465.74Nucleus
OsMSL23LOC_Os04g57530.1Chr41173239041.449.01Nucleus
OsMSL24LOC_Os05g03740.1Chr51002133336.955.85Nucleus
OsMSL25LOC_Os05g48690.1Chr51041134637.449.92Nucleus
OsMSL26LOC_Os06g32944.1Chr6876329132.828.22Nucleus
OsMSL27LOC_Os07g02500.1Chr71104336741.746.61Chloroplast
OsMSL28LOC_Os07g10950.1Chr71437847854.239.06Chloroplast
OsMSL29LOC_Os08g08130.1Chr81170238944.286.2Peroxisome
OsMSL30LOC_Os08g12950.1Chr81254141747.158.92Nucleus
OsMSL31LOC_Os08g37810.1Chr8948131535.067.09Nucleus
OsMSL32LOC_Os08g44690.1Chr8912130334.839.29Nucleus
OsMSL33LOC_Os09g03570.1Chr91932764373.579.06Chloroplast
OsMSL34LOC_Os09g38570.1Chr91011133636.346.58Nucleus
OsMSL35LOC_Os10g33030.1Chr101104336741.76.61Chloroplast
OsMSL36LOC_Os10g41460.1Chr101011133635.648.9Nucleus
OsMSL37LOC_Os11g06410.1Chr111492248355.066.24Nucleus
OsMSL38LOC_Os11g17954.1Chr111545451458.18.3Nucleus
OsMSL39LOC_Os11g38660.1Chr111938564572.977.34Chloroplast
OsMSL40LOC_Os12g06640.1Chr121299143248.776.13Nucleus
OsMSL41LOC_Os12g10550.1Chr12888229533.57.74Nucleus

Share and Cite

MDPI and ACS Style

Li, J.; Zhang, M.; Sun, J.; Mao, X.; Wang, J.; Wang, J.; Liu, H.; Zheng, H.; Zhen, Z.; Zhao, H.; et al. Genome-Wide Characterization and Identification of Trihelix Transcription Factor and Expression Profiling in Response to Abiotic Stresses in Rice (Oryza sativa L.). Int. J. Mol. Sci. 2019, 20, 251. https://doi.org/10.3390/ijms20020251

AMA Style

Li J, Zhang M, Sun J, Mao X, Wang J, Wang J, Liu H, Zheng H, Zhen Z, Zhao H, et al. Genome-Wide Characterization and Identification of Trihelix Transcription Factor and Expression Profiling in Response to Abiotic Stresses in Rice (Oryza sativa L.). International Journal of Molecular Sciences. 2019; 20(2):251. https://doi.org/10.3390/ijms20020251

Chicago/Turabian Style

Li, Jiaming, Minghui Zhang, Jian Sun, Xinrui Mao, Jing Wang, Jingguo Wang, Hualong Liu, Hongliang Zheng, Zhen Zhen, Hongwei Zhao, and et al. 2019. "Genome-Wide Characterization and Identification of Trihelix Transcription Factor and Expression Profiling in Response to Abiotic Stresses in Rice (Oryza sativa L.)" International Journal of Molecular Sciences 20, no. 2: 251. https://doi.org/10.3390/ijms20020251

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop