Next Article in Journal
Exploitation of Drought Tolerance-Related Genes for Crop Improvement
Previous Article in Journal
The Role of Extracellular Vesicles from Human Macrophages on Host-Pathogen Interaction
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 19
10.3390/ijms221910264
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 Identification and Characterization of the Cystatin Gene Family in Bread Wheat (Triticum aestivum L.)

by
Long He
†,
Xuan Chen
†,
Miaoze Xu
,
Tingting Liu
,
Tianye Zhang
,
Juan Li
,
Jian Yang
,
Jianping Chen
* and
Kaili Zhong
*
State Key Laboratory for Quality and Safety of Agro-Products, Institute of Plant Virology, Ningbo University, Ningbo 315211, China
*
Authors to whom correspondence should be addressed.
†
These authors contribute equally to this study.
Int. J. Mol. Sci. 2021, 22(19), 10264; https://doi.org/10.3390/ijms221910264
Submission received: 3 September 2021 / Revised: 18 September 2021 / Accepted: 20 September 2021 / Published: 23 September 2021
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Cystatins, as reversible inhibitors of papain-like and legumain proteases, have been identified in several plant species. Although the cystatin family plays crucial roles in plant development and defense responses to various stresses, this family in wheat (Triticum aestivum L.) is still poorly understood. In this study, 55 wheat cystatins (TaCystatins) were identified. All TaCystatins were divided into three groups and both the conserved gene structures and peptide motifs were relatively conserved within each group. Homoeolog analysis suggested that both homoeolog retention percentage and gene duplications contributed to the abundance of the TaCystatin family. Analysis of duplication events confirmed that segmental duplications played an important role in the duplication patterns. The results of codon usage pattern analysis showed that TaCystatins had evident codon usage bias, which was mainly affected by mutation pressure. TaCystatins may be regulated by cis-acting elements, especially abscisic acid and methyl jasmonate responsive elements. In addition, the expression of all selected TaCystatins was significantly changed following viral infection and cold stress, suggesting potential roles in response to biotic and abiotic challenges. Overall, our work provides new insights into TaCystatins during wheat evolution and will help further research to decipher the roles of TaCystatins under diverse stress conditions.

1. Introduction

Protein hydrolysis in eukaryotic cells is a complex and sophisticated process that is regulated by a series of endogenous or exogenous proteases [1,2]. Proteases can be divided into different families according to their amino acid residues at their reaction sites [3]. Among these proteases, the family C1A proteases, namely papain-like cysteine proteases [PLCPs], and the family C13 proteases, namely legumain-like cysteine proteases [LLCPs], play important roles in various physiological processes [4,5,6,7]. The protease activity of PLCPs and LLCPs is affected by a group of small proteins called cystatins [8]. Cystatins have been reported to be tight and reversible inhibitors of C1A and C13 proteinases in multiple plants [9].
The cystatin family has been identified and characterized in pests, mammals, and plants [10,11,12]. All identified cystatins have three typical conserved motifs, which include a QxVxG motif in the reaction site, glycine residues in the N-terminus, and one tryptophan residue in the C-terminus [13]. These three conserved motifs can directly bind to the active center of the cysteine protease, resulting in the inhibition of catalytic activity [6,14]. In addition, a consensus sequence, namely (LVI)-(AGT)-(RKE)-(FY)-(AS)-(VI)-x-(EDQV), was found to be general for all cystatins identified in plants, which was related to a predicted secondary α-helix structure [15]. Compared to animal cystatins, plant cystatins are a group of proteins with a molecular weight (MW) ranging from 12 to 16 kDa, lacking glycosylation sites and disulfide bonds [6,16]. Several plant cystatins are usually thought to be special inhibitors of LLCPs due to an extended C-terminus which allows their MW to reach 23 kDa [15,17].
The primary function of plant cystatins is the regulation of cysteine proteases, which are involved in various physiological processes including plant growth and development, senescence, seed development and germination, nitrogen fixation, sexual reproduction, embryogenesis, and programmed cell death (PCD) [14,18,19,20,21,22]. More importantly, plant cystatins have been reported to participate in the regulation of plant defense against abiotic or biotic stress, including pathogens, pest attack, heat stress, and exogenous hormone treatments. For example, cystatins have been described to prevent attacks by mites and pathogenic fungi [23,24]. Application of exogenous methyl jasmonate (MeJA) positively modulates plant defense through the induction of cystatin expression against Tilletia indica infection, which causes serious losses to grain yields [25]. Karina et al. cloned a cystatin gene in maize called CC9, which inhibited the plant host immunity response by affecting apoplastic cysteine proteases [22]. A previous report has shown that the application of exogenous abscisic acid (ABA) or heat stress treatments can lead to the induction of AtCYS5 expression and further investigation revealed that the transgenic Arabidopsis lines overexpressing AtCYS5 display enhanced resistance against heat threats [26]. Tan et al. observed that overexpression of MpCYS2, a cystatin gene cloned from Malus prunifolia (Willd.) Borkh in A. thaliana, dramatically enhances drought tolerance [27].
Bread wheat (Triticum aestivum) is the most widely cultivated grain worldwide and supplies approximately one-fifth of the total calories consumed by humans [28]. With the growth of the world’s population, it is predicted that wheat agricultural production needs to be increased by 38% to satisfy the growing demand for food [29]. The cystatin family has been well characterized in a number of plant species such as A. thaliana, soybean, apple, and a variety of cereal crops, including rice, Brachypodium distachyon, and sorghum [10,11,12,30,31,32,33]. However, little is known about the cystatin family in bread wheat. As a complex allohexaploid with a large number of repetitive and transposable elements, bread wheat has one of the largest crop plant genomes (16 Gb genome size; AABBDD genomes), which makes working with bread wheat challenging from a genetics, genomics, and breeding perspective [34]. Fortunately, with the rapid development of genome sequencing technology, a high-quality complete genome assembly and annotation of wheat organized by the International Wheat Genome Sequencing Consortium have been completed, which provides us with a good opportunity to identify and characterize cystatin family members in wheat [35].
In the present study, a genome-wide investigation of the cystatin family in wheat was performed based on the recently released genome of T. aestivum [35]. We identified and characterized 55 members of the cystatin gene family in wheat bread. Furthermore, an overview of gene structures, evolutionary relationships, expansion, and expression levels of the wheat genes from the cystatin gene family is provided. In summary, our work provides a novel viewpoint for the subsequent research into cystatin genes in bread wheat and may contribute to further functional studies of cystatin genes to enhance the resistance of bread wheat against various stresses.

2. Results

2.1. Genome-Wide Identification of the TaCystatin Family

Cystatins were identified and characterized in Arabidopsis thaliana and rice (Oryza sativa) (10,30), and their locus IDs and sequences are listed in Supplementary Table S1. To extract all cystatin members in bread wheat (Triticum aestivum), we performed a genome-wide analysis through local BLSATP using A. thaliana and rice cystatin protein sequences as queries. All candidates were further filtered using Protein family database (Pfam) search and NCBI Batch CDD for function annotations. In this study, we identified 55 cystatins in wheat. The gene ID, location and open reading frame (ORF) length, amino acid length, MW, isoelectric points (PIs), gravy, and subcellular location are listed in Table 1. The ORF of TaCystatins ranged from 318 (TraesCS3D02G416600.1) to 546 (TraesCS1A02G2564000.1) bp in length. The lengths of TaCystatin proteins ranged from 100 (TraesCS2A02G576200.1) to 243 (TraesCS3B02G215400.1) amino acid residues, MW ranged from 11.5 to 26.7 kDa, and the PI ranged from 5.01 (TraesCS5D02G502100.1) to 10.23 (TraesCS2D02G274900.1). In addition, analysis of the subcellular localization showed that more than half of the identified TaCystatins were localized in the extracellular space, while only 10.9%, 26.7%, and 20.0% of the identified TaCystatins were predicted to be localized in the chloroplast, mitochondria, and cytoplasm, respectively.

2.2. Phylogenetic and Classification Analysis of Cystatins

To explore the phylogenetic relationships between TaCystatins and other known plant cystatins, seven A. thaliana, 11 rice, 13 barley, and 55 wheat cystatin protein sequences were used to construct a maximum likelihood phylogenetic tree using MEGA X (Figure 1). Similarly, we also constructed a maximum likelihood phylogenetic tree using only 55 wheat cystatin protein sequences (Figure 2a). Using the NCBI Batch CDD tools, we confirmed that all the identified TaCystatin sequences had one or two conserved cystatin domains (Figure 2b). All identified TaCystatins were clustered into three groups: Group A, Group B, and Group C, which was consistent with previous reports [6,11]. In addition, we found that wheat shows both a greater number of cystatin genes than the other analyzed plants and a greater abundance of all clades (Table 2).

2.3. Duplication Events Analysis of the TaCystatin Family

The cystatin family has been identified and characterized in several flowering plants (12, 15, 30). We found that the number of cystatin genes in bread wheat was the highest among the analyzed plants (Table 2). To better understand why the number of cystatin genes in wheat was quite high, we analyzed the homoeologous groups for the TaCystatin family. In this study, the results showed that 43.65% of the TaCystatin genes displayed in homoeologous groups (1:1:1), representing that three TaCystatins localized on the B, D, and A sub-genome shared high homology, which were also defined as triads (Table 3). Previous studies have reported that 35.8% of wheat genes are present in triads [36,37]. The proportion of homoeologous triads of the cystatin family was higher than that in the whole wheat genome (43.65% vs. 35.8%; Table 3). Hence, the higher proportion suggests that the expansion of TaCystatin may result from wheat polyploidization. However, even with respect to the ploidy level caused by the fact that wheat exhibits complex hexaploidy, the number of cystatin genes in wheat was still much higher than that in rice (Figure S1), which were the model species of monocots. In addition, the proportion of homoeologous groups with gene duplications in one sub-genome (1:1: N, N:1:1, or 1: N:1) was dramatically higher than that in the whole wheat genome (14.55% vs. 5.7%; Table 3), while the proportion of homoeologous groups in which there was one homoeolog (1:1:0, 0:1:1, or 1:0:1) was quite close to that of the whole wheat genome (14.54% vs. 13.2%; Table 3). Thus, our results suggest that the higher number of TaCystatin genes is partly due to the homoeolog retention percentage and gene duplications.
To better investigate the features of the TaCystatin family, we conducted chromosomal location analysis. The wheat A, B, and D sub-genomes contained 16, 18, and 21 cystatin genes, respectively (Figure S2). The wheat chromosome 2 contained the largest number of cystatin genes, with seven cystatin genes on chromosome 2A, six cystatin genes on chromosome 2B, and six cystatin genes on chromosome 2D. No cystatin gene family members were identified on chromosomes 6A, 6B, 6D, and 5B. Chromosomal location analysis revealed that the TaCystatins were unevenly distributed on chromosomes. We further analyzed the distribution of TaCystatin genes within each chromosome. Detailed information is provided in Table S2 and is illustrated in Figure 3a. The results showed that the proportion of TaCystatin distribution in the distal telomeric regions of the chromosomes and the proportion of TaCystatin distribution in the more central regions of the chromosomes were similar (47.45% vs. 54.55%; Table S2). Interestingly, we found that the proportion of the TaCystatin distribution in distal and proximal chromosomal regions varied greatly among Group A (53.33% vs. 46.67%), Group B (0 vs. 100.0%), and Group C (58.62% vs. 41.38%), which corresponded to the number of members of their phylogenetic group. Overall, the TaCystatins belonging to the greater phylogenetic group were inclined to be located in distal telomeric regions of the chromosomes, consistent with previous findings [37]. Tandem and segmental duplications have been recognized as two main factors in the expansion of gene families in plants [38]. Here, we identified eight tandem duplication clusters in the TaCystatin family (Figure 3b). Three tandem duplications occurred on chromosome 3B and the remaining tandem duplications occurred on chromosomes 2B, 2D, 3D, 5A, and 5D. In addition, 45 collinear cystatin gene pairs (Figure 3b) were characterized, which indicated that segmental duplications occupied a significant position in the expansion of the TaCystatin family.

2.4. Analysis of Cystatin Paralogs and Orthologs

To further investigate the evolutionary relationships of the cystatin family, we conducted syntenic analysis using the McscanX software to identify cystatin paralog gene pairs in wheat (T. aestivum) and cystatin orthologous gene pairs between wheat and dicotyledonous plants (Arabidopsis and G. max), monocotyledonous plants (O. sativa and Sorghum bicolor), and wheat relatives (Aegilops tauschii and Triticum dicoccoides). In this study, we identified 45 paralogues in bread wheat (Ta-Ta) (Table S3). No cystatin orthologous gene pairs were observed between Arabidopsis and wheat (Ta-At), and only one cystatin orthologous gene pair was found between G. max and wheat (Ta-Gm) (Figure 4 and Table S4). We also found that 10 and 15 TaCystatin genes were collinear with cystatin genes in rice (O. sativa) and Sorghum bicolor (Ta-Os and Ta-Sb), respectively (Figure 4 and Table S4). Eleven and twenty-five orthologous cystatin gene pairs were identified in wheat with Aegilops tauschii and Triticum dicoccoides (Ta-Aet and Ta-Td), respectively. These results suggest that cystatin genes in wheat are distantly related to those in dicotyledonous species and are most closely associated with those in T. dicoccoides, which might be due to the fact that T. dicoccoides are widely considered to be A-genome and B-genome donors of wheat [39]. The Ka/Ks ratio indicated the selective pressure on plant genes, which can be used to diagnose the evolutional form of the sequence [37]. The divergence time (T) was calculated according to Ks values. We estimated the Ka, Ks, Ka/Ks, and T values for each cystatin gene pair to further decipher the evolutionary trends of the cystatin family. In general, Ka/Ks > 1 represents accelerated evolution with advantageous selection, Ka/Ks = 1 represents neutral selection, and Ka/Ks < 1 represents purifying selection [38]. We found that the Ka/Ks ratios of most paralogous genes (Ta-Ta) and all orthologous genes (Ta-Gm, Ta-Os, Ta-Sb, Ta-Aet, and Ta-Td) were less than 1 (Tables S3 and S4), suggesting that purifying selection plays a more significant role during the evolution of the cystatin family. In addition, the results showed that the divergence time of 45 paralogues (Ta-Ta) ranged from 0.691149 to 39.22863 Mya. The divergence time of orthologues (Ta-Gm, Ta-Os, Ta-Sb, Ta-Aet, and Ta-Td) varied greatly depending on species and the divergence time of orthologues (Ta-Gm) was the longest, while Ta-Td was the shortest (Table S4).

2.5. Codon Usage Pattern Analysis of Cystatin Genes in Plants

Triplet codons play an important role in biological information transmission. Each codon encodes one amino acid, while an amino acid corresponds to at least one and up to six codons, among which the multiple codons encoding the same amino acid are called synonymous codons [40,41]. However, the frequency of synonymous codon usage varies greatly among different species or genes [42]. The codon preference of synonymous codon usage is an important parameter of evolutionary processes. CUB is also important for the level of gene expression, as it affects the translation efficiency and functional differentiation of genes [43,44]. In general, the GC content of the third codon position (GC3) is widely considered to be related to codon usage patterns [45]. Codon usage patterns of cystatin genes in seven species were analyzed using each cystatin CDS sequence and detailed information on codon usage indicators is provided in Table S5. We observed that the average GC proportion of cystatin genes was lower in dicot species than in monocot species (Figure 5a). The results also showed that the average proportion of A/T-terminated codons of cystatin genes in dicot species was relatively higher, while G/C-terminated codons were more common in monocot species (Figure 5a), which was consistent with previous reports [46,47]. The average GC3s and GC content of cystatin genes in rice (O. sativa) were higher than those in other species (Table 4). Compared with the cystatin genes in other species, the average effective number of codons (ENC) in rice was the lowest (Table 4), followed by B. distachyon and bread wheat (T. aestivum). This suggests that the CUB of cystatin genes in rice was the strongest, followed by B. distachyon and bread wheat. It is widely accepted that relative synonymous codon usage (RSCU) can intuitively indicate CUB [48]. RSCU > 1 represents more used codons, while RSCU < 1 means that codons are used less frequently than expected, and RSCU = 1 indicates that codons have no preference [49]. Thus, we performed relative (RSCU) analysis to better understand the codon usage patterns of cystatin genes in seven plant species. We found that the RSCU values of cystatin genes were close in two dicotyledons, namely A. thaliana and G.max, and the RSCU values of that in five monocots were also relatively similar (Figure 5d). Among the five monocotyledon plants, the RSCU values of cystatin genes in common wheat and barley (H. vulgare) were clustered in one group, while those in rice and Sorghum bicolor were clustered in the other (Figure 5d), which might be related to evolutionary relations between these species.
We performed parity rule 2 (PR2) analysis to investigate the bias of the AT and GC composition. We found that there were differences between T, A, G, and C proportions in most cystatin genes (Figure 5c and Figure S3). The results showed that C and A were used more frequently than G and T in common wheat. C and T were used more frequently than G and A in other species except for G3s/(G3s+C3s) in G. max, which suggests that mutation pressure might work in the nucleotide composition. We performed neutral plot analysis to further confirm the presence of other factors for codon preference, except for the base composition. A positive and significant correlation between GC3s and GC12 was detected in the cystatin gene CDS sequences of common wheat (r = 0.8925, p < 0.01) (Figure 5b). The GC3s of cystatin gene CDS sequences in common wheat ranged from 0.588 to 0.992. The slope of the regression line for cystatin genes in common wheat was 0.5941, which suggests that the effect of mutation pressure on the codon preference of cystatin genes in wheat was 59.41%. Neutral plot analysis was also performed for cystatin genes in the six other species. As shown in Figure S4, a positive and significant correlation between GC3s and GC12 was detected in cystatin gene CDS sequences of barley (r = 0.5800, p < 0.01), B. distachyon (r = 0.6667, p < 0.01), S. bicolor (r = 0.6244, p < 0.01), and G. max (r = 0.8032, p < 0.01). In addition, we found that the slope values of the regression line for cystatin genes in A. thaliana and G. max were lower than those in the rest of the analyzed species, which might be due to the fact that both A. thaliana and G. max are dicotyledons. In summary, our results suggest that the codon usage pattern of cystatin genes is relatively conserved within dicotyledonous or monocotyledonous plants and the CUB of cystatin genes in monocots was relatively stronger than that in dicots, consistent with previous findings [46].

2.6. Conserved Motifs and Gene Structure Analysis of TaCystatin

The comparison of the gene exon–intron structure provides novel insights into the evolution and function of gene family members [12]. To investigate the structural features of TaCysatin genes, we analyzed the exon–intron distribution of TaCysatins using TBtools software [50]. The results showed that the number of introns of TaCysatin genes ranged from one to four (Figure 6c). Most members of the TaCysatin genes from Group C or B had only one intron, while all members from Group A had the least number of two introns (Figure 6c). Overall, the exon–intron structures of most genes within the same group were relatively conserved. To characterize the conserved motifs of the identified TaCystatins, we submitted the amino acid sequences of all TaCystatin proteins to the MEME online website. Twenty conserved motifs were predicted, as illustrated in Figure 6b. The results revealed that the motif number of TaCystatins varied from three to nine. Several motifs were prevalent for most TaCystains, while others only existed in certain groups. For example, 98% of TaCystain members contained motif 1, while 87% of TaCystain members contained both motifs 3 and 5. Motif 2 was unique to Group C. Motifs 6 and 11 only existed in Group B. In summary, a number of TaCystatins within the same group shared similar motif structures and the genes of the same motifs’ composition might have similar functions. Detailed information on the 20 conserved motifs is provided in Figure S5.

2.7. Prediction of Cis-Acting Regulatory Elements in Promoter Regions of TaCystatins

Cis-acting regulatory elements located in the promoter region can regulate the expression levels of target genes by binding to transcription factors [51,52]. Cis-acting elements have been reported to be involved in various plant responses to abiotic or biotic stress [53,54,55,56]. To explore the expression pattern of TaCystatin genes, we submitted the 2.0 kb promoter region sequences of TaCystatins to the PlantCARE database website [57]. In summary, 3647 cis-acting regulatory elements were identified. All identified cis-acting regulatory elements could be classified into several categories, including development-related, environmental stress-related, hormone-responsive, light-responsive, promoter and enhancer, site-binding-related, and transcription-related elements (Figure 7a), which suggested that cis-acting elements of TaCystatins play a significant role during wheat growth and development. We identified 278 environmental stress-related elements (Figure 7c). These predicted environmental stress-related elements were involved in temperature, drought, and pathogen responses. In addition, a total of 651 hormone-responsive elements were identified. ABREs responding to ABA and CGTCA-motifs responding to MeJA accounted for the majority of these predicted hormone-responsive elements. In fact, ABA-responsive cis-acting and MeJA-responsive cis-acting elements were characterized in the promoter regions of all the identified TaCystatin genes (Figure 7b).
Cis-acting regulatory element analysis showed that TaCystatin genes were mostly characterized by ABA and MeJA. To investigate the association between the cystatin family in wheat and ABA or MeJA, we randomly selected two members from each phylogenetic group of the TaCystatin family as representatives to measure the expression profiles of six selected TaCystatins upon ABA or MeJA exogenous treatments by quantitative reverse transcription polymerase chain reaction (RT-qPCR). The results showed that all selected TaCystatin genes were sensitive to ABA or MeJA application (Figure 8), which, to some extent, indicated a close relationship between the selected TaCystatin genes’ regulation and ABA together with MeJA. Both ABA and MeJA have been shown to play important roles in plant stress biology [58,59,60,61]. Our results suggest that the expression of TaCystatin genes may be involved in several different stresses.

2.8. Tissue-Specific Expression Analysis of TaCystatins

To comprehensively decipher the functions of TaCystatin genes, we calculated the expression levels of six selected TaCystatin genes in five different tissues of bread wheat by RT-qPCR. Five different organs from which we collected samples contained roots (RO), stem (ST), bottom leaf (BL), middle leaf (ML), and top leaf (TL) in three-leaf-stage bread wheat. The expression of the selected genes in RO were regarded as the mock control. As illustrated in Figure 9, all selected TaCystatin genes were expressed in at least one organ. Most of the selected TaCystatin genes showed distinct expression patterns between different phylogenetic groups. However, their expression pattern was relatively similar within the same phylogenetic group. For instance, the expression levels of all the selected genes were relatively higher in the top leaves and relatively lower in the roots (Figure 9), except for TraesCS2D02G274900.1.
These results suggest that various TaCystatins may be involved in the development of different tissues during various stages.

2.9. Expression Analysis of TaCystatins under Abiotic and Biotic Stress

Previous studies have found that plant cystatins are widely involved in plant growth and development, senescence, and PCD [20,21,22,61,62]. The expression levels of several plant cystatin genes have been reported to be affected by various conditions, including pathogens and cold stress [63,64]. Climate temperature is one of the most important factors affecting wheat production; the Chinese wheat mosaic virus (CWMV), as well as the wheat yellow mosaic virus (WYMV), pose a severe threat to winter wheat production in China [65,66]. Thus, to explore the potential roles of TaCystatins in response to biotic or abiotic stress, we analyzed the effects of viral inoculation (biotic) and gradient temperature treatment (abiotic) on the expression levels of selected TaCystatin genes. We found that the expression levels of all the selected genes were dramatically up-regulated at seven days post inoculation (dpi) with CWMV or WYMV (Figure 10b). The expression levels of TraesCS1B02G322100.1, TraesCS3B02G456800.1, and TraesCS3B02G77600.1 did not change significantly at 10 days post CWMV infection, while significant changes in the expression of TraesCS3B02G215400.1, TraesCS2D02G274900.1, and TraesCS4D02G066000.1 were observed at 10 days post CWMV infection. Regarding 10 days post WYMV infection, TraesCS1B02G322100.1 and TraesCS3B02G77600.1 displayed important expression changes. Meanwhile, there were no significant changes observed in the expression of TraesCS3B02G456800.1, TraesCS3B02G215400.1, TraesCS2D02G274900.1, and TraesCS4D02G066000.1. We found that the expression of TraesCS3B02G456800.1 was dramatically upregulated at 13 days post CWMV inoculation, while its expression levels did not change significantly at 13 days post WYMV inoculation. The expression of TraesCS1B02G322100.1 was significantly increased by WYMV but not by CWMV (Figure 10b). In the final stage of plant responses for inoculation with the virus, the expression of TraesCS2D02G274900.1 and TraesCS4D02G066000.1 were both up-regulated by CWMV infection, while down-regulated by WYMV infection. Moreover, we found that the expression of TraesCS1B02G322100.1 and TraesCS3B02G215400.1 did not respond to WYMV infection, and TraesCS3B02G456800.1 as well as TraesCS3B02G77600.1 did not respond to CWMV infection. Most TaCystatin genes were highly expressed at 8 °C on the seventh day (Figure 10a) compared to that under 15, 20, or 25 °C. Interestingly, the situation became quite complex and diverse on the 10th day. TraesCS2D02G274900.1 showed high expression levels at 25 and 15 °C on the 13th day, while its expression was relatively weakened at 20 °C (Figure 10a). On the 16th day, we observed that the expression levels of some members including TraesCS1B02G322100.1, TraesCS3B02G456800.1, TraesCS3B02G77600.1, and TraesCS4D02G066000.1 were upregulated at high temperatures (20 and 25 °C) (Figure 10a). Additionally, TraesCS3B02G215400.1 showed relatively low expression levels at all time points at high temperatures. Interestingly, we also observed that the expression of TraesCS2D02G274900.1 was up-regulated to a six-fold change on the 13th day under 25 °C and was decreased by a six-fold change between the 13th and 16th day. In addition, the expression of TraesCS3B02G456800.1 also responded to the elevated temperature and its expression was specifically increased on the 13th day at 20 °C. The expression of TraesCS2D02G274900.1 and TraesCS3B02G456800.1 in responding to the elevated temperature might have been due to the fact that they were both in the same phylogenetic group. Overall, although our work suggests that the relative expression levels of TaCystatin genes change greatly under stress conditions, their expression patterns under stress challenges were complex and varied.

3. Discussion

Cystatins in plants, as an intrinsic small protein, have been reported to play important roles in multiple stress-signaling pathways and are widely involved in the response to environmental stress [67,68,69]. Although several previous studies have identified and characterized cystatin members in Arabidopsis, rice, sorghum, barley, Glycine max, and Brachypodium distachyon [10,12,31,32], knowledge of the cystatin family of bread wheat has still been limited until now. As bread wheat (Triticum aestivum) occupies a significant position in the supply of food crops for humans, here, we identified and characterized 55 TaCystatin family members using the latest completion of the wheat genome [35]. According to the phylogenetic and gene structure analysis (Figure 1 and Figure 2a), 55 TaCystatins could be categorized into three groups (Groups A, B, and C), consistent with previous findings [6,30]. We found that the structural domains, gene structures, and motif compositions of the TaCystatin family were relatively conserved within each group (Figure 2b and Figure 6b,c). In addition, a majority of TaCystatins were predicted to be localized in the extracellular components (Table 1), which might be attributed to the fact that cystatins are reversible inhibitors of C1A and C13 proteinases, and most of the C1A and C13 proteinases in plants were localized in the extracellular components [9,70].
The number of cystatins in wheat was the highest among the several plant species (Table 2). This might be due to the fact that wheat had undergone two whole genome duplications and wheat (16 Gb genome size; genomes BBAADD) both is a complex allohexaploid and has a large genome (Table 2). To investigate the reasons for the high number of cystatin genes in wheat, the homoeologous groups for the TaCystatin family were analyzed. Approximately 43.65% of TaCystatin genes could be assigned to 1:1:1 homoeologous groups (Table 3), which was above the average homoeologous retention proportion in wheat (43.65% vs. 35.8%; Table 3) [35]. Moreover, we identified 11 and 25 orthologous gene pairs between TaCystatins and cystatins in A. tauschii and T. dicoccoides (Figure 4c and Table S4), respectively. Previous studies have demonstrated that Aegilops tauschii (genomes DD) is the natural source of D sub-genomes of wheat (genomes BBAADD) and that Triticum dicoccoides (BBAA) is the natural foundation of B and A sub-genomes of wheat (genomes BBAADD) [71]. Together, these results indicate that two whole genome duplications resulting from hybridization might partly be responsible for the abundance of cystatin members in wheat. However, even considering the ploidy level, we found that the number of cystatin genes in wheat was still more than three times that in rice (Figure S1). We found that in wheat, the proportion of cystatin genes with the homeology (1:1:N, N:1:1, or 1:N:1) was (14.55% vs. 5.7%; Table 3) higher than for all other wheat genes. This suggests that gene duplications play a vital role in the expansion of the cystatin family in wheat. Chromosomal locations and synteny analyses were then performed to explore the relationships within the cystatin genes family in wheat more comprehensively. The results showed that 55 TaCystatins were irregularly distributed on chromosome 1, 2, 3, 4, 5, and 7 (Figure S2), and the TaCystatin family members in the larger group were more likely to be in distal telomeric regions (Figure 3a and Table S2), which is also consistent with previous reports [36]. Tandem and segmental duplications have been reported to be the two main causes of gene duplication in plants [37]. In this study, eight tandem duplication clusters and forty-five collinear cystatin gene pairs were identified (Figure 3b), which suggests that both tandem and segmental duplication events were necessary for the expansion of the cystatin family in wheat, while segmental duplications appeared to be more advantageous in duplication patterns.
It is quite common for CUB to occur in the genome, which indicates that genes encoding proteins are not uniformly used. CUB is also considered to be important in gene regulation and molecular evolution [72,73]. To analyze the CUB of the cystatin family in plants, we calculated several representative parameters including CBI, Fop, ENC, GC3s, and GC contents in this study. Among them, G/C in the third base of the cystatin family was more preferable in monocotyledonous plants, whereas T/A in the third base of the cystatin family was more predominant in dicotyledonous plants (Figure 5a), which is consistent with most previous findings [46]. We also deciphered the relative synonymous codon usage (RSCU) of genes from the cystatin family in seven species. The results showed that the RSCU was relatively conserved with monocots and dicots, and the average RSCU of the cystatin family in monocots was higher than that in dicots (Figure 5d). In addition, the results of PR2 together with the neutral plot analysis showed that CUB of cystatin genes in both monocotyledonous and dicotyledonous plants were affected by mature pressure (Figure 5b,c, and Figures S3 and S4). Thus, the results showed that CUB of the genes from the cystatin family was relatively conserved in dicotyledonous or monocotyledonous plants, and the cystatin genes in monocotyledonous plants had enhanced codon preference compared to dicotyledonous plants.
Cystatin genes have shown divergent expression patterns in several plant species. For instance, the expression levels of SbCys15 and SbCys7 in Sorghum bicolor were considerably higher in vegetative tissues than in reproductive tissues, while others were more expressed in reproductive tissues [12]. In this study, we found that most of the selected TaCystatin genes were highly expressed in the top leaves (Figure 9) compared to that in RO, suggesting that they may participate in plant growth and development. We predicted cis-acting regulatory elements to analyze the putative biological functions of TaCystatins in the view of the promoter structure. The results showed that MeJA-responsive cis-acting and ABA-responsive cis-acting elements accounted for most among the predicted hormone-responsive elements. Additionally, both ABA and MeJA-responsive elements of predicted hormone-responsive elements were not absent in the promoter region of all members of the TaCystatin family (Figure 7). Thus, we analyzed the expression of TaCystatins upon ABA and MeJA application. The results suggested that the expression of TaCystatins was significantly changed under ABA or MeJA treatments (Figure 8). As MeJA and ABA have been reported to play important roles in response to different kinds of stress [58], we hypothesized that the TaCystatin family might be involved in the response to stress biology. Moreover, cystatins are thought to be widely used to regulate endogenous processes that respond to different kinds of abiotic or biotic stresses [68,74]. Therefore, we analyzed the expression of the TaCystatin family under biotic stress (viral infection) and abiotic stress (cold treatment). The results showed that the expression levels of the TaCystatin family changed significantly under biotic stress (viral infection) or abiotic stress (cold treatment) (Figure 10) and their relative expression levels varied significantly, which indicates that the TaCystatins might play specific roles under both viral infection and cold stress.

4. Materials and Methods

4.1. Identification of Cystatin Family in Wheat

The cystatin protein sequences of Arabidopsis thaliana and rice (Oryza sativa) were obtained from the Ensemble Plants database (http://plants.ensembl.org/index.html; updated on 12 February 2021) as previously described [10]. The cystatin protein sequences of Glycine max, barley (Hordeum vulgare), Sorghum bicolor, and Brachypodium distachyon were obtained from the Phytozome database (https://phytozome-next.jgi.doe.gov/info; accessed on 12 February 2021) as previously described [12,31,32]. The newly released reference genome of bread wheat (Triticum aestivum) used in this study was downloaded from the Ensemble Plants database (http://plants.ensembl.org/Triticum_aestivum/Info/Index; accessed on 12 February 2021). These cystatin sequences in A. thaliana and rice were used as queries to conduct local BlastP against the latest bread wheat genome (threshold E-value < 1 × 10–10). The hidden Markov model (HMM) profile (PF00031) of the cystatin family was downloaded from the PFAM database (http://pfam.xfam.org/; accessed on 13 February 2021). The cystatin HMM profile (PF00031) was used for functional annotation filters using the HMMER software (version 3.0) [75]. Then, all candidate protein sequences were further filtered using the NCBI Batch Web CD-Search Tool (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/Bwrpsb.cgi; accessed on 13 February 2021) to confirm the structural integrity of the cystatin domain in each tag sequence. In summary, 55 TaCystatins were identified. Detailed information on the TaCystatin family including gene locations, gene length, ORF length, and size were collected from the Ensemble Plants database. The theoretical PI, molecular weight (MW), and grand average of hydropathy (gravy) of TaCystatins were analyzed using the ExPAsy tool (https://web.expasy.org/compute; accessed on 13 February 2021) [76]. The subcellular localization of TaCystatins was predicted using the Plant-mPLoc tool [77].

4.2. Multiple Sequence Alignment and Phylogenetic Analysis

We conducted sequence alignment analysis of the cystatins from bread wheat, rice, Arabidopsis, and barley using MUSCLE in the MEGA X software (Mega Limited, Auckland, New Zealand) with default parameters. Then, all the sequences imported into the MEGA X software were used to construct a maximum likelihood (ML) phylogenetic tree, with a set of 1000 bootstrap replications and the Poisson distribution mode. We also used a similar method to build an ML phylogenetic tree of the TaCystatin protein sequences.

4.3. Gene Duplication Analysis of TaCystatins

McscanX software was used to investigate tandem and segmental duplications in the TaCystatin family [78]. The synteny relationships between several members of the cystatin family in bread wheat and several other species were analyzed using McscanX. Segmental and tandem duplication relationships were virtualized using the Advanced Cicros function of the TBtools software [50]. The Ka/Ks ratios were calculated for tandem duplications using the Ka/Ks Calculator function of the TBtools software and divergence times (T) were estimated on the basis of T = Ks/(2 × 9.1 × 10−9) Mya [38].

4.4. Codon Usage Pattern Analysis

The cystatin CDS sequences longer than 300 bp of bread wheat, A. thaliana, rice, barley, sorghum, Glycine max, and Brachypodium distachyon were obtained to calculate the parameters of codon usage bias using CodonW 1.4.2 software (http://codonw.sourceforge.net/; accessed on 16 March 2021) [79]. These parameters included GC content, GC3s content, frequency of optimal codons, and the codon bias index. GC12 content, relative synonymous codon usage (RSCU), and ENC were calculated using the EMBOSS tool (https://www.bioinformatics.nl/emboss-explorer/; accessedd on 16 March 2021).

4.5. Gene Structure and Motif Analysis

The genome feature format (GTF) file for bread wheat was downloaded from the Ensemble Plants database [80]. We used the gene structure view (advanced) function of the TBtools software to analyze and visualize the gene structures of TaCystatins [50]. We used motif-based sequence analysis tools (MEME) (https://meme-suite.org/meme/; accessed on 18 February 2021) to predict the conserved motifs of TaCystatins, with a set maximum section of up to 20 motifs [81]. Finally, the results were visualized using TBtools software [50].

4.6. Promoter Analysis

We extracted the 2000 bp upstream sequences of the transcription start site of genes from thd TaCystatin family using the GTF file of the bread wheat genome. The obtained sequences were submitted to the PlantCARE website (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/; accessed on 21 March 2021) to predict the putative cis-acting regulatory elements.

4.7. Plant Cultivation and Viral Inoculation

Bread wheat (T. aestivum cv Yangmai 158) seedlings were cultivated in an artificial greenhouse at 23 °C with a 16/8 h (light/dark) photoperiod. Three-leaf-stage wheat seedlings were used for viral inoculation, temperature stress, and hormone treatments. CWMV and WYMV inoculations were performed by applying mechanical friction using in vitro synthesized transcripts, as previously described [82,83]. For temperature stress analysis, wheat seedlings were divided into four groups and laid up in greenhouses with different temperatures (8, 15, 20, and 25 °C), respectively. The plants placed under 8 °C were regarded as the mock control. For hormone treatments, methyl jasmonate (MeJA) and abscisic acid (ABA) were dissolved in 100% ethanol to suitable concentrations as stock solutions and then were diluted with sterile distilled water containing 0.1% Triton X-100. Wheat seedings (YangMai 158) were treated with 100 μmol L−1 MeJA and 100 μmol L−1 ABA, and 0.1% Triton X-100 was regarded as the mock control. Three biological replicates of leaf samples were collected at divergent times for RNA extraction and RT-qPCR analysis.

4.8. RNA Isolation and Real-Time Quantitative Polymerase Chain Reaction (RT-qPCR)

Total RNA was extracted using TRIzol Reagent (Invitrogen, Carlsbad, CA, USA). Then, strand cDNA was synthesized using a First Strand cDNA Synthesis Kit (TOYOBO, Kita-ku, Osaka, Japan). RT-qPCR was performed using an ABI QuantStudio5 Detection System (Applied Biosystems, Foster City, CA, USA) using the Hieff qPCR SYBR Green Master Mix (YEASEN, Shanghai, China). Each treatment was performed using at least three biological replicates, with three technical replicates. The relative expression levels of target genes were calculated in the 2–ΔΔC(t) manner as described in a previous study [84]. The CDC gene was used as the internal reference for each reaction. The primers used for RT-qPCR are provided in Table S6.

5. Conclusions

In this study, we identified and characterized 55 members of the cystatin family in bread wheat, which could be divided into three phylogenetic groups. TaCystatin genes’ structure and the composition of the amino acid motifs in proteins were conserved in each of the three clades of this gene family. The homoeolog retention rate and gene duplication partly explain the expansion of this family and segmental duplications played a predominant role in duplication patterns. Codon usage pattern analysis showed that the TaCystatin family had an obvious codon preference. The expression of selected TaCystatins was organ-specific and greatly changed due to viral infection or cold stress, with several exceptions. Our results will be helpful to attain a comprehensive understanding of the cystatin family in wheat and to investigate the relationships between TaCystatins and responses to biotic or abiotic stress.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/article/10.3390/ijms221910264/s1.

Author Contributions

L.H., K.Z. and J.Y. designed the experiment and prepared the manuscript; K.Z. and J.C. acquired funding; L.H. and X.C. performed experiments; L.H., M.X. and T.L. performed the bioinformatics analysis; T.Z. and J.L. offered valuable feedback to this study. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the China Agriculture Research System from the Ministry of Agriculture of the P.R. of China (CARS-03) and by K.C. Wong Magna Funding Ningbo University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Van der Hoorn, R.A.; Klemenčič, M. Plant proteases: From molecular mechanisms to functions in development and immunity. J. Exp. Bot. 2021, 72, 3337–3339. [Google Scholar] [CrossRef]
  2. Nduwimana, J.; Guenet, L.; Dorval, I.; Blayau, M.; Le Gall, J.Y.; Le Treut, A. Proteases. Ann. Biol. Clin-Paris. 1995, 53, 251–264. [Google Scholar]
  3. Rawlings, N.D. Peptidase inhibitors in the MEROPS database. Biochimie 2010, 92, 1463–1483. [Google Scholar] [CrossRef] [PubMed]
  4. Barrett, A.J. Classification of peptidases. Methods Enzymol. 1994, 244, 1. [Google Scholar]
  5. Misas-Villamil, J.C.; Van der Hoorn, R.A.; Doehlemann, G. Papain-like cysteine proteases as hubs in plant immunity. New Phytol. 2016, 212, 902–907. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Martinez, M.; Cambra, I.; Carrillo, L.; Diaz-Mendoza, M.; Diaz, I. Characterization of the entire cystatin gene family in barley and their target cathepsin L-like cysteine-proteases, partners in the hordein mobilization during seed germination. Plant. Physiol. 2009, 151, 1531–1545. [Google Scholar] [CrossRef] [Green Version]
  7. Li, N.; Zhang, D.-S.; Liu, H.-S.; Yin, C.-S.; Li, X.-X.; Liang, W.-Q.; Yuan, Z.; Xu, B.; Chu, H.-W.; Wang, J.; et al. The rice tapetum degeneration retardation gene is required for tapetum degradation and anther development. Plant. Cell 2006, 18, 2999–3014. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Arai, S.; Watanabe, H.; Kondo, H.; Emori, Y.; Abe, K. Papain-inhibitory activity of oryzacystatin, a rice seed cysteine proteinase inhibitor, depends on the central Gln-Val-Val-Ala-Gly region conserved among cystatin superfamily members. J. Biochem. 1991, 109, 294–298. [Google Scholar]
  9. Brown, W.M.; Dziegielewska, K.M. Friends and relations of the cystatin superfamily—New members and their evolution. Protein Sci. 1997, 6, 5–12. [Google Scholar] [CrossRef] [Green Version]
  10. Martínez, M.; Abraham, Z.; Carbonero, P.; Díaz, I. Comparative phylogenetic analysis of cystatin gene families from Arabidopsis, rice and barley. Mol. Genet. Genom. 2005, 273, 423–432. [Google Scholar] [CrossRef]
  11. Tan, Y.; Wang, S.; Liang, D.; Li, M.; Ma, F. Genome-wide identification and expression profiling of the cystatin gene family in apple (Malusx domestica Borkh.). Plant. Physiol Biochem. 2014, 79, 88–97. [Google Scholar] [CrossRef] [PubMed]
  12. Li, J.; Liu, X.; Wang, Q.; Sun, J.; He, D. Genome-wide identification and analysis of cystatin family genes in Sorghum (Sorghum bicolor L. Moench). PeerJ 2021, 9, e10617. [Google Scholar] [CrossRef]
  13. Stubbs, M.T.; Laber, B.; Bode, W.; Huber, R.; Jerala, R.; Lenarcic, B.; Turk, V. The refined 2.4 A X-ray crystal structure of recombinant human stefin B in complex with the cysteine proteinase papain: A novel type of proteinase inhibitor interaction. EMBO J. 1990, 9, 1939–1947. [Google Scholar] [CrossRef] [PubMed]
  14. Benchabane, M.; Schlueter, U.; Vorster, J.; Goulet, M.C.; Biochimie, D.M.J. Plant cystatins. Biochimie 2010, 92, 1657–1666. [Google Scholar] [CrossRef]
  15. Martinez, M.; Diaz, I. The origin and evolution of plant cystatins and their target cysteine proteinases indicate a complex functional relationship. BMC Evol. Biol. 2008, 8, 198–206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Misaka, T.; Kuroda, M.; Iwabuchi, K.; Abe, K.; Arai, S. Soyacystatin, a novel cysteine proteinase inhibitor in soybean, is distinct in protein structure and gene organization from other cystatins of animal and plant origin. Eur. J. Biochem. 1996, 240, 609–614. [Google Scholar] [CrossRef]
  17. Martinez, M.; Diaz-Mendoza, M.; Carrillo, L.; Diaz, I. Carboxy terminal extended phytocystatins are bifunctional inhibitors of papain and legumain cysteine proteinases. FEBS Lett. 2007, 581, 2914–2918. [Google Scholar] [CrossRef] [Green Version]
  18. Solomon, M.; Belenghi, B.; Delledonne, M.; Menachem, E.; Levine, A. The involvement of cysteine proteases and protease inhibitor genes in the regulation of programmed cell death in plants. Plant. Cell 1999, 11, 431–444. [Google Scholar] [CrossRef]
  19. Kiyosaki, T.; Matsumoto, I.; Asakura, T.; Funaki, J.; Kuroda, M.; Misaka, T.; Arai, S.; Abe, K. Gliadain, a gibberellin-inducible cysteine proteinase occurring in germinating seeds of wheat, Triticum aestivum L., specifically digests gliadin and is regulated by intrinsic cystatins. FEBS J. 2007, 274, 1908–1917. [Google Scholar] [CrossRef]
  20. Díaz-Mendoza, M.; Velasco-Arroyo, B.; González-Melendi, P.; Martínez, M.; Díaz, I. C1A cysteine protease-cystatin interactions in leaf senescence. J. Exp. Bot. 2014, 65, 3825–3833. [Google Scholar] [CrossRef] [Green Version]
  21. Zhao, P.; Zhou, X.M.; Zou, J.; Wang, W.; Wang, L.; Peng, X.B.; Sun, M.X. Comprehensive analysis of cystatin family genes suggests their putative functions in sexual reproduction, embryogenesis, and seed formation. J. Exp. Bot. 2014, 65, 5093–5107. [Google Scholar] [CrossRef] [Green Version]
  22. Van der Linde, K.; Hemetsberger, C.; Kastner, C.; Kaschani, F.; Van der Hoorn, R.A.; Kumlehn, J.; Doehlemann, G. A maize cystatin suppresses host immunity by inhibiting apoplastic cysteine proteases. Plant. Cell 2012, 24, 1285–1300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Carrillo, L.; Martinez, M.; Alvarez-Alfageme, F.; Castañera, P.; Smagghe, G.; Diaz, I.; Ortego, F. A barley cysteine-proteinase inhibitor reduces the performance of two aphid species in artificial diets and transgenic Arabidopsis plants. Transgenic Res. 2011, 20, 305–319. [Google Scholar] [CrossRef]
  24. Martínez, M.; López-Solanilla, E.; Rodríguez-Palenzuela, P.; Carbonero, P.; Díaz, I. Inhibition of plant-pathogenic fungi by the barley cystatin Hv-CPI is not associated with its cysteine-proteinase inhibitory properties. Mol. Plant. Microbe 2003, 16, 876–883. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Dutt, S.; Pandey, D.; Kumar, A. Jasmonate signal induced expression of cystatin genes for providing resistance against Karnal bunt in wheat. Plant. Signal. Behav. 2011, 6, 821–830. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Song, C.; Kim, T.; Chung, W.S.; Lim, C.O. The Arabidopsis phytocystatin AtCYS5 enhances seed germination and seedling growth under heat stress conditions. Mol. Cells 2017, 40, 577–586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Tan, Y.; Li, M.; Ma, F.J.P.C.T.; Culture, O. Overexpression of MpCYS2, a phytocystatin gene from Malus prunifolia (Willd.) Borkh., confers drought tolerance and protects against oxidative stress in Arabidopsis. Plant. Cell Tiss. Org. 2015, 123, 15–27. [Google Scholar] [CrossRef]
  28. Pfeifer, M.; Kugler, K.G.; Sandve, S.R.; Zhan, B.; Rudi, H.; Hvidsten, T.R.; Mayer, K.F.; Olsen, O.A. Genome interplay in the grain transcriptome of hexaploid bread wheat. Science 2014, 345, 1250091. [Google Scholar] [CrossRef] [Green Version]
  29. Ray, D.K.; Mueller, N.D.; West, P.C.; Foley, J.A. Yield trends are insufficient to double global crop production by 2050. PLoS ONE 2013, 8, e66428. [Google Scholar] [CrossRef] [Green Version]
  30. Wang, W.; Zhao, P.; Zhou, X.M.; Xiong, H.X.; Sun, M.X. Genome-wide identification and characterization of cystatin family genes in rice (Oryza sativa L.). Plant. Cell Rep. 2015, 34, 1579–1592. [Google Scholar] [CrossRef] [PubMed]
  31. Yuan, S.; Li, R.; Wang, L.; Chen, H.; Zhang, C.; Chen, L.; Hao, Q.; Shan, Z.; Zhang, X.; Chen, S.; et al. Search for nodulation and nodule development-related cystatin genes in the genome of soybean (Glycine max). Front. Plant. Sci. 2016, 7, 1595. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Subburaj, S.; Zhu, D.; Li, X.; Hu, Y.; Yan, Y. Molecular Characterization and Expression Profiling of Brachypodium distachyon L. Cystatin genes reveal high evolutionary conservation and functional divergence in response to abiotic stress. Front. Plant. Sci. 2017, 8, 743–751. [Google Scholar] [PubMed]
  33. Abraham, Z.; Martinez, M.; Carbonero, P.; Diaz, I. Structural and functional diversity within the cystatin gene family of Hordeum vulgare. J. Exp. Bot. 2006, 57, 4245–4255. [Google Scholar] [CrossRef] [PubMed]
  34. Borrill, P.; Adamski, N.; Uauy, C. Genomics as the key to unlocking the polyploid potential of wheat. New Phytol. 2015, 208, 1008–1022. [Google Scholar] [CrossRef] [PubMed]
  35. International Wheat Genome Sequencing (IWGSC). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 2018, 361, 7191. [Google Scholar] [CrossRef] [Green Version]
  36. Schilling, S.; Kennedy, A.; Pan, S.; Jermiin, L.S.; Melzer, R. Genome-wide analysis of MIKC-type MADS-box genes in wheat: Pervasive duplications, functional conservation and putative neofunctionalization. New Phytol 2020, 225, 511–529. [Google Scholar] [CrossRef] [Green Version]
  37. 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–17. [Google Scholar] [CrossRef] [Green Version]
  38. Hurst, L.D. The Ka/Ks ratio: Diagnosing the form of sequence evolution. Trends Genet. 2002, 18, 486. [Google Scholar] [CrossRef]
  39. Brenchley, R.; Spannagl, M.; Pfeifer, M.; Barker, G.L.; D’Amore, R.; Allen, A.M.; McKenzie, N.; Kramer, M.; Kerhornou, A.; Bolser, D.; et al. Analysis of the bread wheat genome using whole-genome shotgun sequencing. Nature 2012, 491, 705–710. [Google Scholar] [CrossRef] [Green Version]
  40. Hershberg, R.; Petrov, D.A. Selection on codon bias. Annu. Rev. Genet. 2008, 42, 287–299. [Google Scholar] [CrossRef] [Green Version]
  41. Liu, Y.; Li, Y.; Ji, K.K.; Zhu, J.; Ling, P.; Zhou, T.; Fan, L.Y.; Xie, S.Q. Genome-wide codon usage pattern analysis reveals the correlation between codon usage bias and gene expression in Cuscuta australis. Genomics 2020, 112, 2695–2702. [Google Scholar] [CrossRef]
  42. Roberts, R.J. Restriction and modification enzymes and their recognition sequences. Nucleic Acids Res. 1982, 10, 117–144. [Google Scholar] [CrossRef] [Green Version]
  43. Angellotti, M.C.; Bhuiyan, S.B.; Chen, G.; Wan, X.F. Codon usage bias analysis within and across genomes. Nucleic Acids Res. 2007, 35, 132–136. [Google Scholar] [CrossRef] [PubMed]
  44. LaBella, A.L.; Opulente, D.A.; Steenwyk, J.L.; Hittinger, C.T.; Rokas, A. Variation and selection on codon usage bias across an entire subphylum. PLoS Genet. 2019, 15, e1008304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Mazumdar, P.; Binti Othman, R.; Mebus, K.; Ramakrishnan, N.; Ann Harikrishna, J. Codon usage and codon pair patterns in non-grass monocot genomes. Ann. Bot. 2017, 120, 893–909. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Brown, C.M.; Stockwell, P.A.; Trotman, C.N.; Tate, W.P. Sequence analysis suggests that tetra-nucleotides signal the termination of protein synthesis in eukaryotes. Nucleic Acids Res. 1990, 18, 6339–6345. [Google Scholar] [CrossRef]
  47. Kawabe, A.; Miyashita, N.T. Patterns of codon usage bias in three dicot and four monocot plant species. Genes Genet. Syst. 2003, 78, 343–352. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Li, G.; Pan, Z.; Gao, S.; He, Y.; Xia, Q.; Jin, Y.; Yao, H. Analysis of synonymous codon usage of chloroplast genome in Porphyra umbilicalis. Genes Genom. 2019, 41, 1173–1181. [Google Scholar] [CrossRef]
  49. Tao, L.; Zhihua, N.; Mingshan, S.; Liang, C.J. High-frequency codon analysis and its application in codon analysis of tobacco. PLoS ONE 2002, 41, 551–554. [Google Scholar]
  50. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant. 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  51. Wang, X.; Yan, X.; Li, S.; Jing, Y.; Gu, L.; Zou, S.; Zhang, J.; Liu, B. Genome-wide identification, evolution and expression analysis of the aspartic protease gene family during rapid growth of moso bamboo shoots. BMC Genom. 2021, 22, 45. [Google Scholar]
  52. Hernandez-Garcia, C.M.; Finer, J.J. Identification and validation of promoters and cis-acting regulatory elements. Plant. Sci. 2014, 19, 217–218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Yamaguchi-Shinozaki, K.; Shinozaki, K. Organization of cis-acting regulatory elements in osmotic and cold-stress-responsive promoters. Trends Plant. Sci. 2005, 10, 88–94. [Google Scholar] [CrossRef]
  54. Kaur, A.; Pati, P.K.; Pati, A.M.; Nagpal, A.K. In-silico analysis of cis-acting regulatory elements of pathogenesis-related proteins of Arabidopsis thaliana and Oryza sativa. PLoS ONE 2017, 12, e0184523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Wanapu, C.; Shinmyo, A. Cis-regulatory elements of the peroxidase gene in Arabidopsis thaliana involved in root-specific expression and responsiveness to high-salt stress. Ann. N. Y. Acad. Sci. 1996, 782, 107–114. [Google Scholar] [CrossRef]
  56. Vishwakarma, K.; Upadhyay, N.; Kumar, N.; Yadav, G.; Singh, J.; Mishra, R.K.; Kumar, V.; Verma, R.; Pandey, M.; Sharma, S. Abscisic acid signaling and abiotic Stress tolerance in plants. Front. Plant. Sci. 2017, 8, 161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van de Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef] [PubMed]
  58. Lee, S.C.; Luan, S. ABA signal transduction at the crossroad of biotic and abiotic stress responses. Plant. Cell Environ. 2012, 35, 53–60. [Google Scholar] [CrossRef]
  59. Alazem, M.; Lin, N.S. Antiviral roles of abscisic acid in plants. Front. Plant. Sci. 2017, 8, 1760. [Google Scholar] [CrossRef] [Green Version]
  60. Santner, A.; Estelle, M. Recent advances and emerging trends in plant hormone signalling. Nature 2009, 459, 1071–1078. [Google Scholar] [CrossRef]
  61. He, L.; Chen, X.; Yang, J.; Zhang, T.; Li, J.; Zhang, S.; Zhong, K.; Zhang, H.; Chen, J.; Yang, J. Rice black-streaked dwarf virus-encoded P5-1 regulates the ubiquitination activity of SCF E3 ligases and inhibits jasmonate signaling to benefit its infection in rice. New Phytol. 2020, 225, 896–912. [Google Scholar] [CrossRef] [Green Version]
  62. Belenghi, B.; Acconcia, F.; Trovato, M.; Perazzolli, M.; Bocedi, A.; Polticelli, F.; Ascenzi, P.; Delledonne, M. AtCYS1, a cystatin from Arabidopsis thaliana, suppresses hypersensitive cell death. Eur. J. Biochem. 2003, 270, 2593–2604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Hwang, J.E.; Hong, J.K.; Lim, C.J.; Chen, H.; Je, J.; Yang, K.A.; Kim, D.Y.; Choi, Y.J.; Lee, S.Y.; Lim, C.O. Distinct expression patterns of two Arabidopsis phytocystatin genes, AtCYS1 and AtCYS2, during development and abiotic stresses. Plant. Cell Rep. 2010, 29, 905–915. [Google Scholar] [CrossRef] [Green Version]
  64. Chen, P.J.; Senthilkumar, R.; Jane, W.N.; He, Y.; Tian, Z.; Yeh, K.W. Transplastomic Nicotiana benthamiana plants expressing multiple defence genes encoding protease inhibitors and chitinase display broad-spectrum resistance against insects, pathogens and abiotic stresses. Plant. Biotechnol. J. 2014, 12, 503–515. [Google Scholar] [CrossRef] [PubMed]
  65. Rezaei, E.E.; Siebert, S.; Hüging, H.; Ewert, F. Climate change effect on wheat phenology depends on cultivar change. Sci. Rep. 2018, 8, 4891. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Yang, J.; Zhang, T.Y.; Liao, Q.S.; He, L.; Li, J.; Zhang, H.M.; Chen, X.; Li, J.; Yang, J.; Li, J.B.; et al. Chinese wheat mosaic virus-induced gene slencing in monocots and dicots at Low temperature. Front. Plant. Sci. 2018, 9, 1627. [Google Scholar] [CrossRef] [PubMed]
  67. Martinez, M.; Santamaria, M.E.; Diaz-Mendoza, M.; Arnaiz, A.; Carrillo, L.; Ortego, F.; Diaz, I. Phytocystatins: Defense proteins against phytophagous insects and acari. Int. J. Mol. Sci. 2016, 17, 1747. [Google Scholar] [CrossRef] [Green Version]
  68. Wang, Y.; Zhan, Y.; Wu, C.; Gong, S.; Zhu, N.; Chen, S.; Li, H. Cloning of a cystatin gene from sugar beet M14 that can enhance plant salt tolerance. Plant. Sci. 2012, 191, 93–99. [Google Scholar] [CrossRef]
  69. Haq, S.K.; Atif, S.M.; Khan, R.H. Protein proteinase inhibitor genes in combat against insects, pests, and pathogens: Natural and engineered phytoprotection. Arch. Biochem. Biophys. 2004, 431, 145–159. [Google Scholar] [CrossRef]
  70. Shindo, T.; Van der Hoorn, R.A. Papain-like cysteine proteases: Key players at molecular battlefields employed by both plants and their invaders. Mol. Plant. Pathol. 2008, 9, 119–125. [Google Scholar] [CrossRef]
  71. Dubcovsky, J.; Dvorak, J. Genome plasticity a key factor in the success of polyploid wheat under domestication. Science 2007, 316, 1862–1866. [Google Scholar] [CrossRef] [Green Version]
  72. Brule, C.E.; Grayhack, E.J. Synonymous codons: Choose wisely for expression. Trends Genet. 2017, 33, 283–297. [Google Scholar] [CrossRef] [Green Version]
  73. Qian, W.; Yang, J.R.; Pearson, N.M.; Maclean, C.; Zhang, J. Balanced codon usage optimizes eukaryotic translational efficiency. PLoS Genet. 2012, 8, e1002603. [Google Scholar] [CrossRef] [Green Version]
  74. Zhang, X.; Liu, S.; Takano, T. Two cysteine proteinase inhibitors from Arabidopsis thaliana, AtCYSa and AtCYSb, increasing the salt, drought, oxidation and cold tolerance. Plant. Mol. Biol. 2008, 68, 131–143. [Google Scholar] [CrossRef]
  75. Mistry, J.; Finn, R.D.; Eddy, S.R.; Bateman, A.; Punta, M. Challenges in homology search: HMMER3 and convergent evolution of coiled-coil regions. Nucleic Acids Res. 2013, 41, e121. [Google Scholar] [CrossRef] [Green Version]
  76. Wilkins, M.R.; Gasteiger, E.; Bairoch, A.; Sanchez, J.C.; Williams, K.L.; Appel, R.D.; Hochstrasser, D.F. Protein identification and analysis tools in the ExPASy server. Methods Mol. Biol. 1999, 112, 531–552. [Google Scholar] [PubMed]
  77. Chou, K.C.; Shen, H.B. Plant-mPLoc: A top-down strategy to augment the power for predicting plant protein subcellular localization. PLoS ONE 2010, 5, e11335. [Google Scholar] [CrossRef] [Green Version]
  78. Wang, Y.; Tang, H.; Debarry, J.D.; Tan, X.; Li, J.; Wang, X.; Lee, T.H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [Green Version]
  79. Peden, J.F. Analysis of codon usage. PeerJ 2000, 90, 73–74. [Google Scholar]
  80. Bolser, D.M.; Kerhornou, A.; Walts, B.; Kersey, P. Triticeae resources in Ensembl plants. Plant. Cell Physiol. 2015, 56, e3. [Google Scholar] [CrossRef] [Green Version]
  81. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, 202–208. [Google Scholar] [CrossRef]
  82. Zhang, F.; Liu, S.; Zhang, T.; Ye, Z.; Han, X.; Zhong, K.; Yang, J.; Chen, J.; Liu, P. Construction and biological characterization of an infectious full-length cDNA clone of a Chinese isolate of Wheat yellow mosaic virus. Virology 2021, 556, 101–109. [Google Scholar] [CrossRef]
  83. Yang, J.; Zhang, F.; Xie, L.; Song, X.J.; Li, J.; Chen, J.P.; Zhang, H.M. Functional identification of two minor capsid proteins from Chinese wheat mosaic virus using its infectious full-length cDNA clones. J. Gen. Virol. 2016, 97, 2441–2450. [Google Scholar] [CrossRef]
  84. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Ijms 22 10264 g001 550
Figure 1. Phylogenetic tree of cystatins. The maximum likelihood (ML) tree was constructed on the basis of the full-length amino acid sequences of bread wheat (55), Arabidopsis thaliana (7), rice (11), and barley (13) using MEGA X software, with a set of 1000 replications. All cystatins were divided into three phylogenetic groups, and each group is represented by a different color. Triangles, stars, circles, and squares correspondingly indicate bread wheat, rice, barley, and A. thaliana.
Figure 1. Phylogenetic tree of cystatins. The maximum likelihood (ML) tree was constructed on the basis of the full-length amino acid sequences of bread wheat (55), Arabidopsis thaliana (7), rice (11), and barley (13) using MEGA X software, with a set of 1000 replications. All cystatins were divided into three phylogenetic groups, and each group is represented by a different color. Triangles, stars, circles, and squares correspondingly indicate bread wheat, rice, barley, and A. thaliana.
Ijms 22 10264 g001
Ijms 22 10264 g002 550
Figure 2. Phylogenetic and conserved domain analysis of the TaCystatin family. (a) Phylogenetic tree of TaCystatins. The ML tree was constructed on the basis of the full-length amino acid sequences of TaCystatins by MEGA X, with a set of 1000 replications. All TaCystatins were divided into three phylogenetic groups. (b) Conserved domain of 55 TaCystatins.
Figure 2. Phylogenetic and conserved domain analysis of the TaCystatin family. (a) Phylogenetic tree of TaCystatins. The ML tree was constructed on the basis of the full-length amino acid sequences of TaCystatins by MEGA X, with a set of 1000 replications. All TaCystatins were divided into three phylogenetic groups. (b) Conserved domain of 55 TaCystatins.
Ijms 22 10264 g002
Ijms 22 10264 g003 550
Figure 3. Chromosomal distribution of the TaCystatin family. (a) The ratio between the TaCystatin genes in R2A, R2B together with C segments of every chromosome, and R1 together with R3 segments of every chromosome. (b) Genomic location and duplication events analysis of 55 TaCystatin genes. Light-gray lines in the background indicate the synteny blocks within the bread wheat genome. The duplication events are highlighted with different colored lines.
Figure 3. Chromosomal distribution of the TaCystatin family. (a) The ratio between the TaCystatin genes in R2A, R2B together with C segments of every chromosome, and R1 together with R3 segments of every chromosome. (b) Genomic location and duplication events analysis of 55 TaCystatin genes. Light-gray lines in the background indicate the synteny blocks within the bread wheat genome. The duplication events are highlighted with different colored lines.
Ijms 22 10264 g003
Ijms 22 10264 g004 550
Figure 4. Syntenic relationships of the cystatin genes in bread wheat and six other species. (a) Syntenic relationships of cystatins between wheat, Arabidopsis thaliana, and Glycine max. (b) Syntenic relationships of cystatins between wheat, Oryza sativa, and S.bicolor. (c) Syntenic relationships of cystatins between wheat, Aegilops tauschii, and Triticum dicoccoides. Gray lines in the background represent the synteny blocks of wheat and other plants, while the red lines highlight the syntenic cystatin gene pairs.
Figure 4. Syntenic relationships of the cystatin genes in bread wheat and six other species. (a) Syntenic relationships of cystatins between wheat, Arabidopsis thaliana, and Glycine max. (b) Syntenic relationships of cystatins between wheat, Oryza sativa, and S.bicolor. (c) Syntenic relationships of cystatins between wheat, Aegilops tauschii, and Triticum dicoccoides. Gray lines in the background represent the synteny blocks of wheat and other plants, while the red lines highlight the syntenic cystatin gene pairs.
Ijms 22 10264 g004
Ijms 22 10264 g005 550
Figure 5. Codon usage pattern analysis. (a) Vertical axis representing the contents of various bases at the third position of the codons in seven species. (b) Neutrality plot analysis of TaCystatin CDS sequences. (c) Parity rule 2 (PR2) analysis of cystatin CDS sequences in seven species. (d) Heatmap showing relative synonymous codon usage (RSCU) values of cystatin CDS sequences in seven species. Blue-to-red color indicates low-to-high RSCU values of codons.
Figure 5. Codon usage pattern analysis. (a) Vertical axis representing the contents of various bases at the third position of the codons in seven species. (b) Neutrality plot analysis of TaCystatin CDS sequences. (c) Parity rule 2 (PR2) analysis of cystatin CDS sequences in seven species. (d) Heatmap showing relative synonymous codon usage (RSCU) values of cystatin CDS sequences in seven species. Blue-to-red color indicates low-to-high RSCU values of codons.
Ijms 22 10264 g005
Ijms 22 10264 g006 550
Figure 6. Phylogenetic relationship, gene structure, and motif composition analysis of TaCystatins. (a) Phylogenetic tree of TaCystatins. The ML tree was built on the basis of the full-length amino acid sequences of TaCystatins using MEGA X, with a set of 1000 replications. All TaCystatins members were divided into three phylogenetic groups. (b) Architecture of the conserved protein motifs of TaCystatins. Different motifs numbered 1 to 20 are indicated by different colors. (c) Intron/exon structures analysis of TaCystatins. Untranslated regions (UTRs) are represented by green boxes, introns are represented by gray lines, and the coding sequences (CDS) of exons are indicated by yellow boxes.
Figure 6. Phylogenetic relationship, gene structure, and motif composition analysis of TaCystatins. (a) Phylogenetic tree of TaCystatins. The ML tree was built on the basis of the full-length amino acid sequences of TaCystatins using MEGA X, with a set of 1000 replications. All TaCystatins members were divided into three phylogenetic groups. (b) Architecture of the conserved protein motifs of TaCystatins. Different motifs numbered 1 to 20 are indicated by different colors. (c) Intron/exon structures analysis of TaCystatins. Untranslated regions (UTRs) are represented by green boxes, introns are represented by gray lines, and the coding sequences (CDS) of exons are indicated by yellow boxes.
Ijms 22 10264 g006
Ijms 22 10264 g007 550
Figure 7. Prediction analysis of cis-acting regulatory elements in TaCystatins. (a) The number of seven different kinds of cis-acting elements in the promoter region of TaCystatins. (b) The type and position of hormone-responsive elements in TaCystatins. (c) The type and position of environmental stress-related elements in TaCystatins.
Figure 7. Prediction analysis of cis-acting regulatory elements in TaCystatins. (a) The number of seven different kinds of cis-acting elements in the promoter region of TaCystatins. (b) The type and position of hormone-responsive elements in TaCystatins. (c) The type and position of environmental stress-related elements in TaCystatins.
Ijms 22 10264 g007
Ijms 22 10264 g008 550
Figure 8. Expression pattern analysis of TaCystatin genes upon hormone applications. (a) Relative expression levels of six selected TaCystatin genes at 0.5, 2, 4, 6, and 8 h after 100 μmol L−1 ABA treatments. Means ± standard deviations (SE) were deciphered from three biological replicates and each biological replicate had three technical replicates. Note: ns, not significant; **, p < 0.01; and *, p < 0.05 (Student’s t-test). (b) Relative expression levels of six selected TaCystatin genes at 0.5, 2, 4, 6, and 8 h after 100 μmol L−1 MeJA treatments. Means ± standard deviations (SE) were deciphered from three biological replicates and each biological replicate had three technical replicates. Note: ns, not significant; **, p < 0.01; and *, p < 0.05 (Student’s t-test). Plants treated with 0.1% Triton X-100 were used as the mock control.
Figure 8. Expression pattern analysis of TaCystatin genes upon hormone applications. (a) Relative expression levels of six selected TaCystatin genes at 0.5, 2, 4, 6, and 8 h after 100 μmol L−1 ABA treatments. Means ± standard deviations (SE) were deciphered from three biological replicates and each biological replicate had three technical replicates. Note: ns, not significant; **, p < 0.01; and *, p < 0.05 (Student’s t-test). (b) Relative expression levels of six selected TaCystatin genes at 0.5, 2, 4, 6, and 8 h after 100 μmol L−1 MeJA treatments. Means ± standard deviations (SE) were deciphered from three biological replicates and each biological replicate had three technical replicates. Note: ns, not significant; **, p < 0.01; and *, p < 0.05 (Student’s t-test). Plants treated with 0.1% Triton X-100 were used as the mock control.
Ijms 22 10264 g008
Ijms 22 10264 g009 550
Figure 9. Expression pattern analysis of selected TaCystatin genes among different organs. (a) Raw data indicating the relative expression levels of selected TaCystatin genes in distinct organs. Means ± standard deviations (SE) were deciphered from three biological replicates and each biological replicate had three technical replicates. Note: ns, not significant and **, p < 0.01 (Student’s t-test). (b) Heatmap indicating the relative expression of the selected TaCystatin genes in distinct organs. Color scale represents relative expression values, with the color from blue to red indicating low to high expression abundance. Abbreviations: RO, roots; ST, stem; BL, bottom leaf; ML, middle leaf; and TL, top leaf. The expression of RO was regarded as the controls.
Figure 9. Expression pattern analysis of selected TaCystatin genes among different organs. (a) Raw data indicating the relative expression levels of selected TaCystatin genes in distinct organs. Means ± standard deviations (SE) were deciphered from three biological replicates and each biological replicate had three technical replicates. Note: ns, not significant and **, p < 0.01 (Student’s t-test). (b) Heatmap indicating the relative expression of the selected TaCystatin genes in distinct organs. Color scale represents relative expression values, with the color from blue to red indicating low to high expression abundance. Abbreviations: RO, roots; ST, stem; BL, bottom leaf; ML, middle leaf; and TL, top leaf. The expression of RO was regarded as the controls.
Ijms 22 10264 g009
Ijms 22 10264 g010 550
Figure 10. Expression analysis of selected TaCystatins genes under biotic or abiotic stress. (a) Expression profiles of selected TaCystatins genes under different temperatures calculated by RT-qPCR. Means ± standard deviations (SE) were deciphered from three biological replicates and each biological replicate had three technical replicates. Note: ns, not significant; **, p < 0.01; and *, p < 0.05 (Student’s t-test). (b) Expression profiles of selected TaCystatin genes under viral infection analyzed by RT-qPCR. Means ± standard deviations (SE) were deciphered from three biological replicates and each biological replicate had three technical replicates. Note: ns, not significant; ** p < 0.01; and * p < 0.05 (Student’s t-test).
Figure 10. Expression analysis of selected TaCystatins genes under biotic or abiotic stress. (a) Expression profiles of selected TaCystatins genes under different temperatures calculated by RT-qPCR. Means ± standard deviations (SE) were deciphered from three biological replicates and each biological replicate had three technical replicates. Note: ns, not significant; **, p < 0.01; and *, p < 0.05 (Student’s t-test). (b) Expression profiles of selected TaCystatin genes under viral infection analyzed by RT-qPCR. Means ± standard deviations (SE) were deciphered from three biological replicates and each biological replicate had three technical replicates. Note: ns, not significant; ** p < 0.01; and * p < 0.05 (Student’s t-test).
Ijms 22 10264 g010
Table
Table 1. Detailed information of the putative cystatin proteins in bread wheat.
Table 1. Detailed information of the putative cystatin proteins in bread wheat.
Sequence IDGene LocationORF Length (bp)Size (aa)MW (KDa)pIGRAVYSplice VariantsSubcellular Location
TraesCS2A02G576200.177034268877034306730310011.516189.07−0.5511Cytoplasmic
TraesCS3D02G416600.152808447052808478731810511.424068.57−0.0161Cytoplasmic
TraesCS4B02G067000.1598663585986667832110611.719116.83−0.4441Mitochondrial
TraesCS1D02G310300.140655214040655287632410711.787466.74−0.3611Cytoplasmic
TraesCS1B02G322100.154653451454653527732410711.770569.13−0.2931Mitochondrial
TraesCS4D02G252400.142122074742122114834211312.320085.210.0441Extracellular
TraesCS7D02G337700.143112037543112073335111612.609579.40.0391Extracellular
TraesCS3B02G361700.157340358757340465735411712.740769.180.1131Extracellular
TraesCS7A02G341900.150192143650192178935411712.938089.30.1261Extracellular
TraesCS7B02G241000.144863952144863987435411712.836839.40.0141Extracellular
TraesCS4A02G052100.1430313034303188235411713.1815410.11−0.1121Mitochondrial
TraesCS5D02G502100.152980089952980176735711813.526465.01−0.4191Cytoplasmic
TraesCS5A02G557600.170845191870845250635711812.965059.61−0.0571Mitochondrial
TraesCS3D02G020500.16949696695028635711812.498369.980.1171Extracellular
TraesCS4D02G066000.1408244474082480636011913.149979.15−0.2671Extracellular
TraesCS4A02G247900.155867600855867636736011913.212159.15−0.211Mitochondrial
TraesCS2B02G148400.111443071911443108436612112.939729.59−0.061Mitochondrial
TraesCS4B02G396100.167045611567045648336912213.437649.76 1Mitochondrial
TraesCS2A02G126100.1745613767456232037212313.010769.50.0271Extracellular
TraesCS2D02G128900.1749855917498662637212313.0829210.12−0.0591Extracellular
TraesCS3B02G077600.1486679664866852838712814.544675.88−0.2031Cytoplasmic
TraesCS2D02G589700.164595864264595922138712814.584686.04−0.2411Cytoplasmic
TraesCS2A02G578500.177177253877177312438712814.55866.04−0.2821Cytoplasmic
TraesCS2B02G609300.178956395378956454038712814.471536.06−0.1781Cytoplasmic
TraesCS2A02G581500.177356288077356338538712814.590856.4−0.0541Extracellular
TraesCS2B02G611300.179062366479062426138712814.485656.83−0.1351Chloroplast
TraesCS2D02G586800.164492373164492441938712814.485697.78−0.1411Extracellular
TraesCS2B02G502800.169739823969739880238712814.586937.8−0.1631Cytoplasmic
TraesCS2B02G611200.179061284379061340538712814.44678.93−0.161Extracellular
TraesCS2A02G576100.177029035177029092238712814.454688.95−0.1691Extracellular
TraesCS2D02G586900.164497526964497565538712814.515769.14−0.2121Extracellular
TraesCS2A02G575600.177002312577002351138712814.399639.260.0051Extracellular
TraesCS2D02G477600.157951058957951116338712814.553059.33−0.1091Cytoplasmic
TraesCS3B02G456600.169800019569800063839313013.825747.670.0681Chloroplast
TraesCS3D02G325100.143796653843796795942914215.714276.14−0.0561Cytoplasmic
TraesCS5D02G502600.152984504352984627542914215.685238.8−0.1361Extracellular
TraesCS5A02G487700.165775683365775816142914216.081899.69−0.1421Mitochondrial
TraesCS5D02G502200.152981566352981700043514416.30717.95−0.0921Extracellular
TraesCS2A02G275900.145548109445548171243514415.0515310.010.1831Chloroplast
TraesCS2B02G293600.140940047340940117143514414.7881710.160.2021Chloroplast
TraesCS2D02G274900.134430305134430363743514414.9373610.230.1661Chloroplast
TraesCS3B02G456800.169819487569819531243814515.559558.55−0.241Extracellular
TraesCS3A02G421000.166252008766252086545615116.171317.66−0.1641Extracellular
TraesCS4D02G025600.1109997101100009338412713.307299.170.2271Extracellular
TraesCS3B02G456700.169812225769812271545915216.275377.66−0.1991Extracellular
TraesCS3D02G416700.152809485952809531745915216.102267.66−0.1141Extracellular
TraesCS3D02G416500.152773567852773613946215316.200318.43−0.151Extracellular
TraesCS3B02G456500.169798325369798371446215316.199288.45−0.1121Extracellular
TraesCS5A02G487600.165772516665772596650716819.18526.16−0.2831Extracellular
TraesCS1D02G255800.134806343634806724251317018.153719.42−0.1151Extracellular
TraesCS1B02G267100.146999934147000314754318019.016677.79−0.0181Extracellular
TraesCS1A02G256400.144892677544893074454618119.114798.62−0.0442Chloroplast
TraesCS3D02G189800.117748278517748600073224326.605376.37−0.2491Extracellular
TraesCS3A02G185800.121811745921812038873224326.605376.37−0.2491Extracellular
TraesCS3B02G215400.125635715125636009573224326.748516.38−0.2881Extracellular
TraesCS3B02G215400.125635715125636009573224326.748516.38−0.2881Extracellular
Table
Table 2. Number of cystatin proteins in seven different species.
Table 2. Number of cystatin proteins in seven different species.
LineageOrganismGenome SizeTotal Number of Cystatin Proteins
DicotsArabidopsis thaliana (2n)135 Mb7
Glycine max (2n)1.15 Gb20
MonocotsOryza sativa (2n)500 Mb11
Sorghum bicolor (2n)700 Mb22
Hordeum vulgare (6n)1.35 Gb13
Brachypodium distachyon (2n)300 Mb23
Triticum aestivum (6n)15.8 Gb55
Table
Table 3. Groups of homeologous cystatin genes in bread wheat.
Table 3. Groups of homeologous cystatin genes in bread wheat.
Homoeologous (A:B:D)All Wheat Genes 1Classes 2Number of GroupsNumber of Genes% of Total TaCystatins
ABC
1:1:135.8%21582443.65%
n: 1:1, 1: n:1, or 1:1: n &5.7% 112814.55%
1:1:0, 1:0:1, or 0:1:113.2%2 24814.55%
Orphans37.1% 1111.81%
Other ratios8%11241425.45%
& n > 1. 1 All wheat genes distributed among homeologous groups of the whole wheat genome according to IWGSC. 2 The number of cystatin family members within each phylogenetic group (A, B, and C).
Table
Table 4. Codon usage indicators of the cystatin family in seven different species. Abbreviations: CBI, codon bias index; Fop, frequency of optimal codons; ENC, effective number of codon; and GC3s, contents of G or C bases at the third position of the codons; and GC content, the contents of the G and C bases of the codons.
Table 4. Codon usage indicators of the cystatin family in seven different species. Abbreviations: CBI, codon bias index; Fop, frequency of optimal codons; ENC, effective number of codon; and GC3s, contents of G or C bases at the third position of the codons; and GC content, the contents of the G and C bases of the codons.
Species NameCBIFopENCGC3sGC Content
Triticum aestivum0.093 0.470 42.185 0.774 0.611
Arabidopsis thaliana0.010 0.423 52.233 0.509 0.473
Brachypodium distachyon0.109 0.479 40.929 0.810 0.626
Hordeum vulgare0.096 0.474 41.722 0.770 0.604
Glycine max0.028 0.433 50.735 0.595 0.508
Oryza sativa0.106 0.477 38.076 0.857 0.669
Sorghum bicolor0.094 0.468 42.536 0.794 0.634
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

He, L.; Chen, X.; Xu, M.; Liu, T.; Zhang, T.; Li, J.; Yang, J.; Chen, J.; Zhong, K. Genome-Wide Identification and Characterization of the Cystatin Gene Family in Bread Wheat (Triticum aestivum L.). Int. J. Mol. Sci. 2021, 22, 10264. https://doi.org/10.3390/ijms221910264

AMA Style

He L, Chen X, Xu M, Liu T, Zhang T, Li J, Yang J, Chen J, Zhong K. Genome-Wide Identification and Characterization of the Cystatin Gene Family in Bread Wheat (Triticum aestivum L.). International Journal of Molecular Sciences. 2021; 22(19):10264. https://doi.org/10.3390/ijms221910264

Chicago/Turabian Style

He, Long, Xuan Chen, Miaoze Xu, Tingting Liu, Tianye Zhang, Juan Li, Jian Yang, Jianping Chen, and Kaili Zhong. 2021. "Genome-Wide Identification and Characterization of the Cystatin Gene Family in Bread Wheat (Triticum aestivum L.)" International Journal of Molecular Sciences 22, no. 19: 10264. https://doi.org/10.3390/ijms221910264

APA Style

He, L., Chen, X., Xu, M., Liu, T., Zhang, T., Li, J., Yang, J., Chen, J., & Zhong, K. (2021). Genome-Wide Identification and Characterization of the Cystatin Gene Family in Bread Wheat (Triticum aestivum L.). International Journal of Molecular Sciences, 22(19), 10264. https://doi.org/10.3390/ijms221910264

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