[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Interaction between Growth Regulators Controls In Vitro Shoot Multiplication in Paulownia and Selection of NaCl-Tolerant Variants
Next Article in Special Issue
Crosstalk between Ca2+ and Other Regulators Assists Plants in Responding to Abiotic Stress
Previous Article in Journal
Improvement of Selected Morphological, Physiological, and Biochemical Parameters of Roselle (Hibiscus sabdariffa L.) Grown under Different Salinity Levels Using Potassium Silicate and Aloe saponaria Extract
Previous Article in Special Issue
A Possible Mode of Action of Methyl Jasmonate to Induce the Secondary Abscission Zone in Stems of Bryophyllum calycinum: Relevance to Plant Hormone Dynamics
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 11
Issue 4
10.3390/plants11040496
plants-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
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Analysis and Characterization of the Proline-Rich Extensin-like Receptor Kinases (PERKs) Gene Family Reveals Their Role in Different Developmental Stages and Stress Conditions in Wheat (Triticum aestivum L.)

by
Mahipal Singh Kesawat
1,2,
Bhagwat Singh Kherawat
3,
Anupama Singh
1,
Prajjal Dey
1,
Snehasish Routray
4,
Chinmayee Mohapatra
4,
Debanjana Saha
5,
Chet Ram
6,
Kadambot H. M. Siddique
7,
Ajay Kumar
8,
Ravi Gupta
9,
Sang-Min Chung
10 and
Manu Kumar
10,*
1
Department of Genetics and Plant Breeding, Faculty of Agriculture, Sri Sri University, Cuttack 754006, Odisha, India
2
School of Biological Sciences and Institute for Molecular Biology and Genetics, Seoul National University, Seoul 08826, Korea
3
Krishi Vigyan Kendra, Bikaner II, Swami Keshwanand Rajasthan Agricultural University, Bikaner 334603, Rajasthan, India
4
Department of Entomology and Plant Pathology, Faculty of Agriculture, Sri Sri University, Cuttack 754006, Odisha, India
5
Department of Biotechnology, Centurion University of Technology and Management, Bhubaneshwar 752050, Odisha, India
6
ICAR-Central Institute for Arid Horticulture, Bikaner 334006, Rajasthan, India
7
The UWA Institute of Agriculture, The University of Western Australia, Perth, WA 6009, Australia
8
Agriculture Research Organization, Volcani Center, Department of Postharvest Science, Rishon Lezzion 50250, Israel
9
College of General Education, Kookmin University, Seoul 02707, Korea
10
Department of Life Science, Dongguk University, Dong-gu, Ilsan, Seoul 10326, Korea
*
Author to whom correspondence should be addressed.
Plants 2022, 11(4), 496; https://doi.org/10.3390/plants11040496
Submission received: 19 January 2022 / Revised: 8 February 2022 / Accepted: 9 February 2022 / Published: 11 February 2022
(This article belongs to the Special Issue What Makes the Life of Stressed Plants a Little Easier? Defense Mechanisms against Adverse Conditions)
Browse Figures
Figure 1
<p>Phylogenetic analysis of TaPERK proteins with Arabidopsis (15), rice (8), and soybean (16). The phylogenetic analysis was executed using the ClustalW program as well as MEGAX software by the neighbor-joining method and bootstrap values of 1000 replicates. The numbers on the nodes indicate the bootstrap values. Distinct groups are represented by the different colors.</p> "> Figure 2
<p>Genomic distribution of identified PERK genes on the 21 chromosomes of wheat and within the three sub-genomes. (<b>A</b>) Schematic representations of the chromosomal distribution of PERK genes on the 21 chromosomes of wheat and the name of the gene on the right. The colored circles on the chromosomes indicate the position of the PERK genes. The chromosome numbers of the three sub-genomes are indicated at the top of each bar. (<b>B</b>) Distribution of PERK genes in the three sub-genomes. (<b>C</b>) Distribution of PERK genes across 21 chromosomes, Un: unaligned contig.</p> "> Figure 2 Cont.
<p>Genomic distribution of identified PERK genes on the 21 chromosomes of wheat and within the three sub-genomes. (<b>A</b>) Schematic representations of the chromosomal distribution of PERK genes on the 21 chromosomes of wheat and the name of the gene on the right. The colored circles on the chromosomes indicate the position of the PERK genes. The chromosome numbers of the three sub-genomes are indicated at the top of each bar. (<b>B</b>) Distribution of PERK genes in the three sub-genomes. (<b>C</b>) Distribution of PERK genes across 21 chromosomes, Un: unaligned contig.</p> "> Figure 3
<p>Syntenic relationships of TaPERK genes between <span class="html-italic">Aegilops tauschii</span>, <span class="html-italic">Brachypodium distachyon</span>, and <span class="html-italic">Oryza sativa</span>. The gray lines in the background represent the collinear blocks within <span class="html-italic">Triticum aestivum</span> and other plant genomes, while the red lines highlight the syntenic PERK gene pairs.</p> "> Figure 4
<p>Diagrammatic representation of the exon–intron organization of the TaPERK genes. Yellow boxes represent exons, untranslated regions (UTRs) are indicated by blue boxes, and black lines represent introns. The lengths of the boxes and lines are scaled based on gene length. The exon and intron sizes can be estimated using the scale at the bottom.</p> "> Figure 5
<p>Conserved motifs of TaPERK genes elucidated by MEME. Up to 10 motifs were shown in different colors. (<b>A</b>) Colored boxes representing different conserved motifs with different sequences and sizes. (<b>B</b>) Sequence logo conserved motif of the wheat PERK proteins. The overall height of each stack represents the degree of conservation at this position, while the height of individual letters within each stack indicates the relative frequency of the corresponding amino acids. The sequence of each motif, combined <span class="html-italic">p</span>-value, and length are shown on the left side of the figure. MEME Parameters: number of repetitions, any; maximum number of motifs, 10; optimum motif width, between 6 and 50.</p> "> Figure 5 Cont.
<p>Conserved motifs of TaPERK genes elucidated by MEME. Up to 10 motifs were shown in different colors. (<b>A</b>) Colored boxes representing different conserved motifs with different sequences and sizes. (<b>B</b>) Sequence logo conserved motif of the wheat PERK proteins. The overall height of each stack represents the degree of conservation at this position, while the height of individual letters within each stack indicates the relative frequency of the corresponding amino acids. The sequence of each motif, combined <span class="html-italic">p</span>-value, and length are shown on the left side of the figure. MEME Parameters: number of repetitions, any; maximum number of motifs, 10; optimum motif width, between 6 and 50.</p> "> Figure 6
<p>Multiple sequence alignment of the TaPERK protein sequences. The conserved protein tyrosine kinase domain is boxed in red. Colored and shaded amino acids are chemically similar residues. Dashes indicate gaps introduced to maximize the alignment of the homologous region. * indicates positions which have a single, fully conserved residue.</p> "> Figure 6 Cont.
<p>Multiple sequence alignment of the TaPERK protein sequences. The conserved protein tyrosine kinase domain is boxed in red. Colored and shaded amino acids are chemically similar residues. Dashes indicate gaps introduced to maximize the alignment of the homologous region. * indicates positions which have a single, fully conserved residue.</p> "> Figure 7
<p>Cis-acting regulatory elements (CAREs) in the promoter region of the TaPERK genes family. The CAREs analysis was performed with a 2kb upstream region using PlantCARE online server. The different numbers of <span class="html-italic">cis</span>-regulatory elements represent different colors. (<b>A</b>) Hormone-responsive elements, stress-responsive elements, growth and development-related elements, light-responsive elements, and other elements with unknown functions are differentiated by color. (<b>B</b>) Most commonly occurring CAREs in TaPERKs.</p> "> Figure 8
<p>Heatmap representing expression profile of the TaPERK genes at various developmental stages. Columns represent genes, and rows represent different developmental stages. TPM values were used directly to create the heatmaps. The “z” nomenclature refers to Zadok’s growth stage.</p> "> Figure 9
<p>Quantitative real-time PCR analysis of selected TaPERK genes in response to drought stress (DS), heat stress (HS), and cold stress to verify RNA seq data. The wheat actin gene was used as the internal control to standardize the RNA samples for each reaction. Asterisks indicate significant differences compared with control. Bars represent results of Tukey’s HSD test at the &lt;0.05 and &lt;0.001 level (* <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> lies in between the values of 0.05 and 0.001, and *** <span class="html-italic">p</span> &lt; 0.001). Error bars show standard deviation. Data are mean ± SD (<span class="html-italic">n</span> = 3).</p> "> Figure 10
<p>Protein–protein interaction analysis of TaPERKs proteins. Protein–protein interaction network produced by STRINGV9.1, each node represents a protein, and each edge represents an interaction, colored by evidence type. The figure highlights the connections between differentially represented proteins.</p> ">
Versions Notes

Abstract

:
Proline-rich extensin-like receptor kinases (PERKs) are a class of receptor kinases implicated in multiple cellular processes in plants. However, there is a lack of information on the PERK gene family in wheat. Therefore, we identified 37 PERK genes in wheat to understand their role in various developmental processes and stress conditions. Phylogenetic analysis of PERK genes from Arabidopsis thaliana, Oryza sativa, Glycine max, and T. aestivum grouped them into eight well-defined classes. Furthermore, synteny analysis revealed 275 orthologous gene pairs in B. distachyon, Ae. tauschii, T. dicoccoides, O. sativa and A. thaliana. Ka/Ks values showed that most TaPERK genes, except TaPERK1, TaPERK2, TaPERK17, and TaPERK26, underwent strong purifying selection during evolutionary processes. Several cis-acting regulatory elements, essential for plant growth and development and the response to light, phytohormones, and diverse biotic and abiotic stresses, were predicted in the promoter regions of TaPERK genes. In addition, the expression profile of the TaPERK gene family revealed differential expression of TaPERK genes in various tissues and developmental stages. Furthermore, TaPERK gene expression was induced by various biotic and abiotic stresses. The RT-qPCR analysis also revealed similar results with slight variation. Therefore, this study’s outcome provides valuable information for elucidating the precise functions of TaPERK in developmental processes and diverse stress conditions in wheat.

1. Introduction

Protein phosphorylation is a key post-translational modification that regulates cell signaling networks and cellular processes in response to internal and external environmental stimulation through the reversible regulation of protein function by activation or deactivation, formation of protein complexes, and determination of the subcellular location of proteins [1,2,3]. Phosphorylation is the most common event in which the phosphoryl group transfers from adenosine triphosphate to hydroxyl residue of the protein substrate [4]. Although plants deploy receptor kinases at the cell surface to perceive, the signal, generated in the sudden changing environment, activates the various signaling pathways and regulates the growth, reproduction, and response against diverse stresses [2,5]. Receptor kinases are the most prominent gene family in crop species, including Arabidopsis, rice, maize, soybean, cotton, and sorghum [6,7,8]. For example, the Arabidopsis receptor kinase gene family encompasses ~610 members, and their homologs have been identified and characterized in many plant species [7,9,10,11]. The biological functions of these predicted receptor kinase genes remain to be elucidated. However, receptor kinases play a crucial roles in cell differentiation, pollen tube growth, pollen development, symbiosis, pathogen recognition, phytohormone response, signal transduction, self-incompatibility and response towards internal and external stimuli [8,9,12,13,14,15,16]. A few studies have identified and characterized the ligands that activate specific receptor kinases, as well as some signaling components [8,13]. Receptor kinases bind to different kinds of biomolecules such as polypeptides, steroids, carbohydrates, and cell wall components. These receptor kinases perceive and transduce the signals across the plasma membrane via diverse signaling complexes, which have been developed during the long course evolution of complex multicellular organisms [2,9].
Receptor kinases were divided into different groups based on the motifs structure in their extracellular domains [9,17]. For instance, the leucine-rich repeat receptor kinase family comprises the main class of receptor kinases in plants, including BAK1 and BRI1, which have been implicated in brassinosteroid signaling [18,19,20]. Proline-rich extensin-like receptor kinases (PERKs) are one of the main classes of receptor kinases. Fifteen PERK genes have been found in the Arabidopsis genome. However, their functions are poorly understood [6,12]. Most of the AtPERK gene family members are ubiquitously expressed, while few genes are specifically expressed [6]. For instance, AtPERK1 is broadly expressed, whereas the expression of AtPERK2 is mainly observed in rosette leaf veins, stems, and pollen [21]. In addition, the expression of AtPERK8 and AtPERK13 were detected in root hairs [6,22,23]. The expression of AtPERK5, AtPERK6, AtPERK7, AtPERK11, and AtPERK12 were highly elevated in pollen. However, unnoticeable expressions were observed in the sporophytic tissues [23,24,25]. Furthermore, AtPERK4 regulates the root growth function at an early stage of ABA signaling by perturbing calcium homeostasis in Arabidopsis [14]. A few studies have demonstrated that increased concentrations of calcium in the cells also enhances antioxidant enzyme activities and eventually regulates the lipid peroxidation of cell membranes and stomatal apertures [26,27,28]. The PERKs suppress the accumulation of reactive oxygen species (ROS) in the root, which is necessary for root hair growth [29,30]. MAPK cascade is an essential regulator of high light-induced Cu/Zn SODs and anti-PERK antibodies from animals, used to detect the presence of homologous proteins such as MPK3 and MPK6 in plants [30,31]. AtPERK5 and AtPERK12 are essential for the pollen tube growth in Arabidopsis [16]. Furthermore, AtPERK8, AtPERK9, and AtPERK10 negatively regulate root growth in Arabidopsis [23]. PERK1 rapidly induces early perception and response to a wound in Chinese cabbage [12]. Antisense suppression of BnPERK1 has exhibited various growth defects, such as amplified secondary branching, loss of apical dominance, and defects in floral organ formation. At the same time, the overexpression line showed increased lateral shoot production, seed set, and unusual deposition of callose and cellulose in Brassica napus [21]. A PERK-like receptor kinase specifically interacts with the nuclear shuttle protein (NSP), led viral infection, and positively regulates the NSP function in cabbage leaf curl virus and geminivirus [32].
Plants are sessile organisms, and constantly face fluctuating environmental conditions and various biotic and abiotic stresses during growth and development [33,34,35]. For example, wheat is an important cereal crop cultivated worldwide [36,37], and its quality and productivity are largely influenced by different biotic and abiotic factors [38,39]. With the recent advent of sequencing technology, a rapid increase in sequenced plant genomes has been accessed in the past few years [40]. However, identifying the genes in plant species’ genomes is now a great challenge, particularly in terms of their structure to functionally characterization [41,42]. For example, the wheat genome sequenced, completed and identified 124,201 genes [43]. Thus, this project’s completion has made it possible to complete genome-wide analysis and identification of the PERK gene family in wheat. We performed a comprehensive analysis of 37 PERK genes using several computational approaches in this work.
Furthermore, phylogenetic analysis, physical and biochemical properties, exon/intron, conserved motifs, chromosomal distribution, subcellular localization, gene duplication, Ka/Ks values, synteny analysis, and three-dimensional (3D) structure were also determined. In addition, tissue-specific expression profiles and responses to diverse stress conditions were also examined for the TaPERK genes. The outcome of the present study will be helpful in the detailed understanding of the TaPERK gene’s role in plant growth, development, and survivability under different stress conditions.

2. Results

2.1. Identification of TaPERK in Wheat

In this study, we identified 37 PERK genes in the wheat genome using various computational approaches (Table 1).
This number is relatively higher than the earlier reported PERK genes in Arabidopsis, soybean, rice, sorghum, maize, and cotton (Table 2).
This might be due to the higher chromosome number and big size of the wheat genome, which indicates that the PERK genes underwent a substantial expansion in wheat. In addition, wheat is derived from the hybridization of three progenitor genomes: A, B, and D. The TaPERK family had protein lengths ranging from 401–1052, and amino acid with molecular weight (MW) 44.37–113.73 kDa for TaPERK23 and TaPERK7, respectively. The isoelectric point (pI) ranged from 5.22 and 9.04 for TaPERK34 and TaPERK3, respectively. We also plotted the MW of TaPERK with their pI to understand the MW distribution of different TaPERK proteins (Figure S1). The plots showed that most of the TaPERKs had similar MW and pI. Hence, pI values ranged from acidic to basic, and the heaviest TaPERK was over twice the weight of the lightest. Furthermore, the grand average of hydropathy index values ranged from −0.09 to −0.645, indicating that TaPERK proteins are hydrophilic in nature. Moreover, the subcellular localization prediction of TaPERK proteins indicated that most of the TaPERKs were situated on the plasma membrane (Table 1).
To understand the origins and evolutionary dynamics between plant species PERKs, the phylogenetic tree was produced with TaPERKs, AtPERKs OsPERKs, and GmPERK proteins (Table S2). The phylogenetic analysis revealed that TaPERK proteins were classified into eight groups (Figure 1).
Group III was the biggest with 14 members, while Group I, II, IV, V, VI, VII, and VIII contained 10, 3, 0, 6, 0, 4, and 0 members, respectively (Figure S2).

2.2. Chromosomal Distribution, Gene Duplication, and Synteny Analysis

To map the chromosomal distribution of the identified TaPERK genes in wheat, corresponding to chromosomal locations of PERK genes were determined using the PhenGram online server. The TaPERK genes were found on 17 wheat chromosomes (Figure 2A and Table 1). TaPERK genes showed a higher presence on A sub-genomes (Figure 2B). Maximum TaPERK genes (Fourteen) were located on the chromosomes of the A sub-genome.
The B and D sub-genome had a minimum number of TaPERK genes (Eleven). Five TaPERKs were mapped on chromosomes 3A and 3B (Figure 2C). The lowest number of TaPERKs was detected on the chromosomes 1A, 1B, 4B, 5A, 5B, and 7B (single gene, respectively). On the contrary, none of the TaPERK genes were located on the chromosomes 4D, 5D, and 6D. In addition, one TaPERK was located on an unaligned contig. Thus, all the PERK family members were uniformly distributed on the wheat’s A, B, and D sub-genome.
To explore why the wheat was polyploidy with the largest genome, we further investigated the duplication events in the TaPERK gene family. The phylogenetic analysis of the TaPERK genes also revealed many duplication events (Figure S3). We observed that 26 PERK genes in wheat involved duplication events (Figure S4 and Table S3), indicating expanding the PERK gene family in wheat. Furthermore, to examine the selective pressure on the duplicated TaPERK genes, we analyzed the synonymous substitution (Ks), non-synonymous (Ka), and the Ka/Ks ratios for the 13 TaPERK genes pairs (Table S3). The value of Ka/Ks = 1 indicates that genes underwent a neutral selection; <1 denotes negative selection or purifying, and >1 suggests a positive selection [44]. The Ka/Ks values for all 11 gene pairs were <1, which indicates that TaPERK genes experienced a robust purifying selection pressure with slight alteration after duplication. However, 2 gene pairs, TaPERK1/TaPERK2 and TaPERK17/TaPERK26, had more than 1, which suggests that two pairs of TaPERK genes experienced a positive selection (Table S3). These findings showed the conserved evolution of TaPERKs.
To further elucidate the synteny relationships of TaPERK genes with wheat relatives and other model plants, including B. distachyon, Ae. tauschii, T. dicoccoides, O. sativa, and A. thaliana, multiple collinearity scan tools were run to identify the orthologous genes between genomes of these plant species (Figure 3 and Table S4).
We found 35, 30, 61, 66, and 83 orthologous gene pairs between TaPERKs with other PERK genes in B. distachyon, Ae. tauschii, T. dicoccoides, O. sativa, and A. thaliana, respectively. The results showed that 26, 24, 37, 42, and 66 TaPERK genes were collinear with PERK genes in B. distachyon, Ae. tauschii, T. dicoccoides, O. sativa, and A. thaliana, respectively. Some of the TaPERK genes had five pairs of orthologous genes, for example; TaPERK5, TaPERK6, TaPERK8, TaPERK9, and TaPERK11, while few of the TaPERK genes had four pairs of orthologous genes, TaPERK7, TaPERK10, TaPERK12, TaPERK15, TaPERK20, TaPERK28, TaPERK32 and TaPERK35 that might have played an essential role in the evolution of PERK genes. Thus, these results indicated that PERK genes in wheat-derived from a common ancestor.

2.3. Exon/Intron Structure and Motif analysis of TaPERK Genes

To elucidate the structural character of the TaPERK genes, the exon/intron organization and conserved motifs (Figure 4) of TaPERK genes were examined.
Exon–intron analysis showed that the TaPERK gene family greatly varied in terms of gene structure. For instance, most TaPERK genes contain 3–23 introns. Maximum twenty-three introns were detected in the TaPERK5, TaPERK6, TaPERK7, TaPERK8, TaPERK9, TaPERK11, TaPERK12, TaPERK32, and TaPERK35, while TaPERK31 had three introns (Figure S5). Furthermore, we also analyzed the conserved motif of TaPERK genes using the Multiple Em for Motif Elicitation (MEME) webserver. Eventually, ten well-preserved motifs were found in 37 TaPERK genes (Figure 5A,B).
Furthermore, the TaPERK gene family was detected by the presence of the tyrosine kinase domain (Pfam PF07714), and all TaPERKs consist of at least one tyrosine kinase domain (Table S5) involved in signal transduction. Furthermore, to understand the biological function of TaPERK genes in wheat, 3D protein models of all TaPERKs were produced using a phyre2 webserver. TaPERKs 3D protein structure had two distinct subdomains, a smaller N-terminal lobe and a bigger C-terminal lobe connected by a small hinge loop (Figure S6B). In addition, protein sequence alignment also showed that all TaPERK proteins consisted of a conserved tyrosine kinase domain (Figure 6 and Figure S6A).
This result will help understand and explain the substrate specificity and molecular function of TaPERK genes in activating the PERK signal transduction pathway.

2.4. Cis-Acting Regulatory Elements (CAREs) Analysis of TaPERK Genes

To further understand the function of TaPERK genes, upstream 2000 bp sequences from the transcription start site of TaPERKs were analyzed using the PlantCARE web server. This analysis revealed that the promoter region of TaPERKs gene families contained the multiple cis-elements related to phytohormones, developmental processes, and different stresses (Figure 7A and Table S6).
The TaPERKs gene family consists of five hormone response elements, including the auxin response element (AuxRE), gibberellin response element (GARE), methyl jasmonate response element (MeJARE), abscisic acid response element (ABRE), and salicylic acid response element (SARE). The response elements belong to light responses, MeJARE, defense, and stress response and ABRE were the most abundant CAREs in the TaPERK gene family (Figure 7B). This result indicates that TaPERKs play a crucial role in plant growth and development.
Furthermore, TaPERKs contain cis-elements related to zein metabolism, endosperm expression, circadian control, meristem expression, seed-specific, and cell cycle regulation. Thus, the CAREs found in the TaPERK gene family indicate that TaPERKs might be participating in a wide range of biological processes. Furthermore, various types of CAREs in the TaPERK genes suggest that these genes might be involved in diverse developmental processes. Therefore, these results provide valuable insights to understand the regulatory mechanism of the TaPERK gene family in response to phytohormone, defense, stress, and various developmental processes.

2.5. Gene Ontology (GO) Enrichment of TaPERK Genes

Gene ontology (GO) assists in understanding the biological function of any genes by comparing their sequence similarity with the known function of genes and gene products with other species [40,42]. All TaPERKs were successfully annotated and allotted GO terms using AgriGO, and further verified using eggNOG-Mapper (Figure S7; Table S7 and Table S8), giving almost the same results as AgriGO. In the biological process category, TaPERK genes were enriched in cell communication (GO:0007154), signaling (GO:0023052), cellular process (GO:0009987), and regulation of biological process (GO:0050789) categories (Figure S7A). In the cellular component category, TaPERK displayed enrichment in the cell (GO:0005623), cell junction (GO:0030054), and membrane (GO:0016020) (Figure S7B). Furthermore, subcellular localization prediction (Table 1) also provided indistinguishable results. In the molecular function category, molecular transducer activity (GO:0060089) and catalytic activity (GO:0003824) were the most prevalent category which was primarily involved in signal transduction (Figure S7C). Apart from cell communication and signaling, the GO term analysis also suggested a variety of roles of TaPERK genes, such as maintenance of dormancy, tissue development, organ formation, post-embryonic organ development, gametophyte development, seedling development, and regulation of developmental process and metabolism. Thus, these results demonstrate that TaPERK genes play a critical role in plant growth and development.

2.6. Expression Profiling of TaPERK Genes in Various Developmental Stages and under Diverse Stress Conditions

To investigate the precise function of TaPERK genes, the expression pattern of TaPERK genes was examined during different developmental and in diverse stress conditions. The TPM values of all TaPERKs were retrieved from the wheat gene expression database. These TPM values were directly used to generate the PCA and heatmaps (Figure S8A,B, Figure 8 and Figure 9).
To examine the expression pattern of TaPERKs, five tissues from three different developmental stages were taken in this work. The TaPERK genes displayed differential induction among the different tissues; for example, TaPERK2, TaPERK4, TaPERK13, TaPERK21, TaPERK23, and TaPERK25 exhibited induction at the spike z39 stage, while TaPERK15, TaPERK20, TaPERK27, TaPERK29, and TaPERK36 exhibited induction at spike z65 stage (Figure 8).
The expression of TaPERK5, TaPERK6, TaPERK7, TaPERK8, TaPERK9, TaPERK11, TaPERK12, TaPERK17, TaPERK22, TaPERK26, TaPERK28, TaPERK32, TaPERK33, TaPERK34, TaPERK35, and TaPERK37 were elevated in roots at z13 and z39 stage, respectively. TaPERK10, TaPERK30, and TaPERK31 were also up-regulated in the leaf at the z23 and z71 stages, respectively. In addition, TaPERK13, TaPERK14, TaPERK23, TaPERK24, and TaPERK31 showed induction at the grain z71 stage. TaPERK18 and TaPERK24 showed higher expression at the stem z30 stage, whereas TaPERK24, TaPERK30, and TaPERK31 expression was raised at the stem z65 stage (Figure 8). These results showed that the TaPERK gene family members might be involved in developing different tissues and stages.
Expression patterns of TaPERKs were also investigated under the different stress conditions, including septoria tritici blotch (STB), stripe rust, powdery mildew, drought, and heat stress. The expression of several members of the TaPERK gene family was elevated in biotic and abiotic stress (Figure S9). The expression of TaPERK9, TaPERK13, TaPERK15, TaPERK17, TaPERK20, TaPERK22, TaPERK23, TaPERK26, TaPERK28, TaPERK33, TaPERK35, and TaPERK36 were induced during the septoria tritici blotch, while the expression of TaPERK2, TaPERK6, TaPERK7, TaPERK8, TaPERK10, TaPERK11, TaPERK12 TaPERK27, TaPERK35, and TaPERK37 were significantly raised upon the powdery mildew infection. TaPERK1, TaPERK8, TaPERK16, TaPERK21, TaPERK24, TaPERK25, TaPERK29, TaPERK30, TaPERK31, and TaPERK34 were highly elevated during the stripe rust infection. In the case of abiotic stress, the expression profile indicates the expression of a few members of the TaPERK family, for instance, TaPERK3, TaPERK4, TaPERK18, and TaPERK32 were raised during the initial hours of heat stress. It seems that the TaPERK family does not participate in drought stress. However, only TaPERK4 and TaPERK18 genes were elevated during the combined drought and heat stress (Figure S9). The expression level of TaPERK14, TaPERK19, TaPERK24, and TaPERK29 was significantly raised during cold stress. Furthermore, the expression patterns of a few selected TaPERK genes were validated through RT-qPCR, and the results displayed nearly similar expression patterns (Figure 9). Overall, these results demonstrated that different TaPERK genes respond to diverse stress conditions.

2.7. Protein–Protein Network Analysis of the TaPERK Family Genes

A protein network was produced using the STRING online webserver to examine the interactions between TaPERKs and other T. aestivum proteins (Figure 10 and Table S9).
We found eighteen TaPERKs interacting with 10 different wheat proteins according to the STRING results. TaPERK29 can interact with seven other wheat proteins (Traes_3AL_394642923.1, Traes_3AL_5E8DEE3E8.1, Traes_3B_514AAB5F3.1, Traes_3B_9EBD47B52.1, Traes_4AS_E28B34320.1, Traes_6BL_DFDCD5B11.1 and Traes_4DL_05BF7F181.1), which were cGMP-dependent protein kinase/PKG II, protein of unknown function (DUF1645) and BRASSINOSTEROID INSENSITIVE 1, and play critical roles in the Brassinosteroids signaling. TaPERK2 and TaPERK4 can interact with four other wheat proteins (Traes_7AL_5E0DD589E.1, Traes_4DL_E447FD9FD.1, Traes_6BL_DFDCD5B11.1 and Traes_7DL_909EA97B3.1) which were cGMP-dependent protein kinase/PKG II and non-specific serine/threonine-protein kinase. cGMP-dependent protein-kinase is a phosphorylated diverse biologically important pathway [45,46,47]. PKG is activated by cGMP and has been implicated in the regulation of cell division, nucleic acid synthesis response to biotic stress, stomata closure during osmotic stress, and development of adventitious roots [45,46,47,48,49]. These results provide important insight for further elucidating the complex biological functions of TaPERK genes.

3. Discussion

PERKs are a class of receptor kinases that have been implicated during various stages of growth and developments in plants, including cell differentiation, pollen tube growth, pollen development, symbiosis, pathogen recognition, phytohormone response, signal transduction, self-incompatibility, and response to internal and external stimuli [8,9,12,13,14,15,16]. PERKS gene family members have also been identified in other plant species, such as 15 genes in Arabidopsis, 8 in O. sativa, 23 in Z. mays, 16 in G. max, 15 in S. bicolor, 15 in G. arboreum, 16 in G. raimondii, and 33 from G. hirsutum [6,7,12]. However, this is the first time we have identified the PERK gene family in the wheat genome. Many studies have reported PERKS genes in ancient land plants, which have expanded during evolutionary processes [6,7,15]. However, in this study, we identified 37 TaPERK genes in the wheat genome (Table 1), which, upon phylogenetic analysis, classified the TaPERK gene family into eight subfamilies or groups (Figure 1). Phylogenetic analysis revealed that groups III and VII were monocot-specific TaPERKs, while groups IV, VI, and VIII contained dicot-specific TaPERKs (Figure 1). The evolution of this type of gene indicates the monocot’s specific functions that might play an essential role in establishing physiological and morphological development [40,42,50]. Although, TaPERK genes were distributed into the well-known rice, Arabidopsis, and the soybean cluster, indicating that TaPERKs might be derived from a common ancestor. In addition, most of the TaPERKs showed orthologous relationships with rice, Arabidopsis, and soybean PERKs.
Furthermore, the phylogenetic tree also displayed that all subfamilies have an expanded number of members (Figure 1 and Figure S2), suggesting that the duplication of TaPERKs results from a long course of evolution. Similar results were reported in Arabidopsis, B. rapa, and cotton [6,7,15]. Collectively, these results demonstrated a lineage-specific expansion of TaPERKs via the partial alteration of the genome to adapt to internal and external environments during evolution [40,42,50,51].
The wheat PERKS gene family was widely expanded and had comparatively more PERKs than the previously reported PERKs in A. thaliana, O. sativa, G. max, S. bicolor, Z. mays, G. max, G. arboreum, G. raimondii, and G. hirsutum [6,7,12]. Many previous studies have demonstrated that polyploidy enabled numerous plant species to adapt to adverse environmental conditions [44,52,53]. Mostly, polyploidy is linked with gene duplication and, in our study, we also found that tandem, segmental and whole-genome duplication was the critical driving force responsible for the duplication of TaPERK genes. Segmental duplication is the fundamental drive factor, and occurs in numerous plant genomes during evolution consisting of several duplicated chromosomal blocks [54]. For instance, several Arabidopsis gene families experienced coherent evolutionary dynamics directed to expanding the gene family [55,56]. Moreover, several gene families, including cotton GRAS, RH2FE3, MADS-Box, MIKC-Type, YABBY, WOX, sesame heat shock proteins, and soybean WRKY, experienced segmental expansion and whole-genome duplication events [7,57,58,59,60,61,62]. The chromosomal map of TaPERK genes revealed that the 37 TaPERKs were unequally distributed throughout chromosomes, excluding chromosome 6 (Figure 2). The gene number on each chromosome varied from one to five: chromosomes 3A and 3B had five genes; chromosome 3D had four genes; chromosome 2A and 2D contained three genes; chromosomes 1D, 2B, 4A,7A and 7D had two genes; and 1A, 1B, 4B, 5A, 5B and 7B consisted of a single gene. Hence, uneven distribution of the TaPERK genes on the 17 chromosomes of wheat indicates probable gene addition or loss via whole genome or segmental duplication events and errors during genome sequencing and assembly. Gene duplication analysis showed 13 pairs of duplicated genes, which shared high sequence similarity at the nucleotide level. The duplicated pairs were TaPERK1:TaPERK2, TaPERK14:TaPERK24, TaPERK27:TaPERK29, TaPERK30:TaPERK31, TaPERK17:TaPERK26, TaPERK7:TaPERK12, TaPERK32:TaPERK35, TaPERK8:TaPERK10, TaPERK9:TaPERK11, TaPERK15:TaPERK20, TaPERK34:TaPERK37, TaPERK13:TaPERK23, and TaPERK16:TaPERK21. Furthermore, the Ka/Ks value of 11 gene pairs was <1, indicating that TaPERK genes experienced a robust purifying selection pressure (Figure S3 and Table S3). However, two gene pairs, TaPERK1:TaPERK2 and TaPERK17:TaPERK26, had more than 1, suggesting that two pairs of TaPERK genes underwent a positive selection. Therefore, these results indicate that TaPERK genes were not changed much in function after duplication and exhibited the conserved evolution of TaPERK genes. Qanmber and colleagues (2019) also reported similar results in cotton. Furthermore, ten gene pairs were the results of segmental duplications in cotton [7]. Furthermore, 146 out of 149 duplicated gene pairs had a Ka/Ks ratio of <1.0, and only three duplicated gene pairs displayed more than 1, which indicates the positive selection pressure. Similar type gene duplication events were also described in the BrPERKs genes [15]. Our gene duplication analysis also demonstrated that the TaPERK gene duplication events were similar, as previously reported in the cotton and Brassica rapa [7,15]. Thus, these results showed that segmental and whole-genome duplications might play a critical role in the evolution and expansion of the PERK genes in wheat.
To further elucidate the synteny relationships of TaPERK genes with wheat relatives and other model plants, we identified 35, 30, 61, 66, and 83 orthologous gene pairs between TaPERKs with other PERK genes in B. distachyon, Ae. tauschii, T. dicoccoides, O. sativa, and A. thaliana, respectively (Figure 3 and Table S4). Additionally, Ae. speltoides (BB, diploid) and Ae. tauschii (DD, diploid) were the foundation of B and D subgenomes of wheat. The synteny relationship displayed that nine orthologous gene pairs between Ae. tauschii with a wheat D subgenome were found on the same chromosomes with two on 1D, one on 2D, four on 3D, and two on 7D (Figure 3 and Table S4). Furthermore, twenty-two orthologous gene pairs between T. dicoccoides with a wheat AABB subgenome were detected on the same chromosomes with one on 1A, two on 2A, five on 3A, two on 4A, one on 5A, two on 7A, one on 1B, one on 2B, four on 3B, one on 4B, one on 5B, and one on 7B (Figure 3 and Table S4). These findings suggest that PERK genes might have come from Ae. tauschii and T. dicoccoides during natural hybridization events. Furthermore, more orthologous gene pairs were found in T. aestivum with A. thaliana and O. sativa, which exhibited that TaPERK and other PERKS genes might be derived from these orthologous genes during evolution.
The gene structure analysis of TaPERKs revealed that TaPERKs greatly varied in gene structure. The majority of the TaPERK genes contained more than five exons, except for TaPERK31 with four exons, while TaPERK30 had five exons (Figure 4). A maximum of twenty-four exons were detected in TaPERK5, TaPERK6, TaPERK7, TaPERK8, TaPERK9, TaPERK11, TaPERK12, TaPERK32, and TaPERK35. Furthermore, a maximum of twenty-three introns were found in TaPERK5, TaPERK6, TaPERK7, TaPERK8, TaPERK9, TaPERK11, TaPERK12, TaPERK32, and TaPERK35, while TaPERK31 had a minimum of three introns (Figure S5). The size of an intron is a critical player that affects the gene size; for example, a notable difference in gene size was found between the biggest gene TaPERK17 (4 kb) and the smallest gene TaPERK23 (2.1 kb), and this was mainly caused by the total intron length (4 kb vs. 1.1 kb). Many studies have shown the significance of introns in the evolution of numerous plant genes [63,64]. Several gene families had less, lack, or more introns in their gene families [7,59,65,66]. The exon and intron differences might be due to deletion/insertion events, which would predict the evolutionary processes [67]. All PERK gene family members in cotton had no introns, indicating that GhPERK genes might have evolved comparatively quickly [7].
Furthermore, it has been established that gene families containing larger or more introns can acquire new functions during evolution processes. There were more intron gains than losses in the plant lineages and chordates, while in arthropods and fungi, losses prevailed over gains [63,64,67]. In our study, almost all TaPERK genes had more and larger introns. Hence, we can speculate that PERK genes gained new functions during evolution in wheat. Furthermore, conserved motif analysis showed ten different types of motif compositions amidst the TaPERK proteins. We observed that five motifs were found in all the TaTERK proteins (Motif 1, 2, 3, 4 and 6), and proteins of the same subfamilies usually shared the same motifs and were more conservative. Thus, we hypothesize that proteins of the same subfamilies may have the same function.
Additionally, amino acid sequence alignment of TaPERK with other plant species PERK proteins also showed that all TaPERK proteins consisted of a conserved tyrosine kinase domain (Figure 6 and Figure S6). The amino acid residues of PERK were highly conserved in rice, Arabidopsis, soybean, and wheat, which might be helpful to find the pattern of PERK protein sequence conservation in different plant species. Yang and colleagues also found that YABBY and WOX gene families were evolutionarily conserved in cotton [57,58]. Furthermore, 3D protein structure analysis revealed that TaPERKs had two distinct subdomains, a smaller N-terminal lobe, and a more prominent C-terminal lobe connected by a small hinge loop (Figure S6A,B). These findings will be helpful to understand and explain the substrate specificity and molecular function of TaPERK genes in activating the PERK signal transduction pathway.
The cis-acting regulatory element in the promoter plays an important role in regulating and functioning genes [68]. The promoter region of TaPERKs gene families contains the multiple cis-acting elements related to plant hormones, growth, development, defense, and stress-related functions (Figure 7A and Table S6). We predicted more than eight CAREs in the promoter region of each TaPERK (Table S6). A total of 15 CAREs related to light response were detected, including AE-box, Box 4 and ATCT motif, chs-Unit 1 m1, TCT-motif, I-box, chs-CMA1a, chs-CMA2a, GA-motif, GATA-motif, LAMP-element and TCCC-motif, ACE and GT1-motif, Sp1 and 3-AF1 binding site [69,70]. We also detected the six CAREs related to growth and development, such as MSA-like (cell cycle regulation), GCN4-motif (endosperm expression), O2-site (zein metabolism regulation), CAT-box (meristem expression), RY-element (seed-specific regulation) and CAAAGATATC-motif (circadian control) [71,72]. In addition, we also found the CARE related to phytohormone response, for instance, CGTCA-motif (MeJA-responsive element), ABRE (abscisic acid-responsive element), TCA-element (salicylic acid responsiveness), TGA-motif (auxin-responsive element), P-box, and GARE-motif (gibberellin-responsive element). The MeJA-responsive element was predicted in most TaPERK genes except TaPERK8, TaPERK11, TaPERK23, and TaPERK28. Moreover, we also predicted that other cis-elements had been involved in different stress conditions, such as LTR (low-temperature responsiveness), MBS (drought inducibility), and TC-rich repeats (defense and stress responsiveness) in the TaPERK promoters. Several studies have demonstrated that light plays a crucial role in plant growth and development processes [73]. Several CAREs related to low temperature, fungal elicitors, stress and defense, auxins, MeJA, gibberellin, ethylene abscisic acid, and the salicylic acid-responsive element were also predicted in GhPERK and BrPERK gene promoter regions [7,15]. In this study, almost all TaPERK genes contained the multiple CAREs involved in plant growth and the responses to diverse stress. GhPERK8, GhPERK 9, GhPERK12, GhPERK23, GhPERK27, and GhPERK29 expression levels were elevated upon exposure to plant hormones such as indole-3-acetic acid, gibberellin, salicylic acid, and MeJA; however, the expression level of GhPERK5 declined [7]. PERK4 regulates the root growth function at an early stage of ABA signaling by perturbing calcium homeostasis in Arabidopsis [14]. PERK1 rapidly induced early perception and response to a wound stimulus in Chinese cabbage [12]. Antisense suppression of BnPERK1 exhibited various growth defects such as amplified secondary branching, loss of apical dominance, and defects in floral organ formation. At the same time, the overexpression line showed increased lateral shoot production, seed set, and unusual deposition of callose and cellulose in Brassica napus [21]. Collectively, these results showed that PERKS gene family members might regulate diverse biological processes, responses to phytohormones, and work against different biotic and abiotic stress. Of course, this needs to be established by experimental studies in the near future. Therefore, these data provide the valuable information to understand TaPERKs’ function in plant growth and development, response to phytohormones, and different stresses.
Receptor kinases play a critical role in different biological processes and responses to internal and external stimuli [8,9,12,13,14,15,16,21,25]. Different TaPERK genes displayed differential expressions in various tissues. For example, TaPERK2, TaPERK4, TaPERK13, TaPERK21, TaPERK23, and TaPERK25 exhibited induction at the spike z39 stage, while TaPERK15, TaPERK20, TaPERK27, TaPERK29, and TaPERK36 exhibited induction at spike z65 stage (Figure 8). The expression of TaPERK5, TaPERK6, TaPERK7, TaPERK8, TaPERK9, TaPERK11, TaPERK12, TaPERK17, TaPERK22, TaPERK26, TaPERK28, TaPERK32, TaPERK33, TaPERK34, TaPERK35, and TaPERK37 were elevated in the roots at the z13 and z39 stages, respectively. TaPERK10, TaPERK30, and TaPERK31 were also up-regulated in the leaf at the z23 and z71 stages, respectively. In addition, TaPERK13, TaPERK14, TaPERK23, TaPERK24, and TaPERK31 showed induction at the grain z71 stage. TaPERK18 and TaPERK24 showed higher expression at the stem z30 stage, whereas TaPERK24, TaPERK30, and TaPERK31 expressions were raised at the stem z65 stage. PERKs proteins have been involved in various developmental processes, including cell differentiation, pollen tube growth, pollen development, symbiosis, pathogen recognition, phytohormone response, signal transduction, and self-incompatibility [9,12,14,16,17,21]. AtPERK gene family members are ubiquitously expressed, while few genes are specifically expressed [6]. For instance, AtPERK1 is broadly expressed, whereas AtPERK2 is mainly expressed in rosette leaf veins, stems, and pollen [21,32]. AtPERK8 and AtPERK13 expression were found in the root hairs [6,22,23]. In addition, AtPERK5, AtPERK6, AtPERK7, AtPERK11, and AtPERK12 expression were up-regulated in the pollens [23,24,25]. Furthermore, PERK4 modulates the root tip growth at an early stage of ABA signaling via the disruption of calcium homeostasis in Arabidopsis [14]. Some researchers have shown that increased calcium concentration in the cells also enhances antioxidant enzyme activities to regulate the lipid peroxidation of cell membranes and stomatal aperture [26,27,28]. AtPERK5 and AtPERK12 play an essential role in pollen tube growth in Arabidopsis [16]. Furthermore, AtPERK8, AtPERK9, and AtPERK10 negatively regulate root growth in Arabidopsis [23]. The expression level of twelve GhPERK genes was significantly elevated in leaves and ovule development in cotton [7]. BrPERK genes were differentially expressed in various tissues of Chinese cabbage, but some BrPERK genes were specially expressed in reproductive organs [15]. Our GO analysis also indicated the critical roles of the TaPERK gene in the cell (Figure S7A–C). Thus, this spatial and temporal expression of the TaPERK genes suggests that these PERKs might have an essential function in different wheat tissue.
Receptor kinases are crucial in plant adaptations and responses to internal and external stimuli [8,9,12,13,14,15,74]. Our results also showed that several TaPERK gene family members’ expressions were elevated in different stress conditions (Figure S9). TaPERK9, TaPERK13, TaPERK15, TaPERK17, TaPERK20, TaPERK22, TaPERK23, TaPERK26, TaPERK28, TaPERK33, TaPERK35, and TaPERK36 were induced during the septoria tritici blotch, while the expressions of TaPERK2, TaPERK6, TaPERK7, TaPERK8, TaPERK10, TaPERK11, TaPERK12 TaPERK27, TaPERK35, and TaPERK37 were significantly raised after powdery mildew infection. TaPERK1, TaPERK8, TaPERK16, TaPERK21, TaPERK24, TaPERK25, TaPERK29, TaPERK30, TaPERK31, and TaPERK34 were up-regulated during the stripe rust infection. Furthermore, TaPERK3, TaPERK4, TaPERK18, and TaPERK32 were induced during the initial hours of heat stress. However, none of the genes were expressed in drought stress. It seems that the TaPERK family does not participate in drought stress, and only TaPERK4 and TaPERK18 genes were elevated during the combined drought and heat stress (Figure S9). Furthermore, TaPERK14, TaPERK19, TaPERK24, and TaPERK29 expression levels were significantly elevated in cold stress. Most TaPERKs respond similarly to biotic and abiotic stress; hence, all stress-responsive genes cluster together (Figure S8B). Several PERKS genes in A. thaliana, G. hirsutum, and B. rapa were responsive to diverse abiotic stresses, including cold, salt, heat, and PEG, indicating that PERK genes play a critical role in other plant species to adapt to different stress conditions [6,7,15,23]. PERK1 rapidly induced early perception and response to a wound stimulus in Chinese cabbage [12]. A PERK-like receptor kinase specifically interacts with the nuclear shuttle protein (NSP), led viral infection, and positively regulates the NSP function in cabbage leaf curl virus and geminivirus [32]. The expression profile of TaPERK genes under different stresses indicated that they might participate in the diverse biotic and abiotic stress tolerance in wheat. Therefore, these findings demonstrated that TaPERK genes respond to various stresses, and this might be used for breeding wheat lines to develop stress-tolerant varieties in wheat.
Many studies have shown that receptor kinases at the cell surface perceive a sudden changing environment that activates the various signaling pathways, regulating growth, reproduction, and response to diverse stress conditions [2,5,75]. Our protein–protein network analysis revealed that eighteen TaPERKs interacted with 10 different wheat proteins (Figure 10 and Table S9). TaPERK29 can interact with seven other wheat proteins (Traes_3AL_394642923.1, Traes_3AL_5E8DEE3E8.1, Traes_3B_514AAB5F3.1, Traes_3B_9EBD47B52.1, Traes_4AS_E28B34320.1, Traes_6BL_DFDCD5B11.1 and Traes_4DL_05BF7F181.1), which were cGMP-dependent protein kinase/PKG II, protein of unknown function (DUF1645) and BRASSINOSTEROID INSENSITIVE 1, playing a critical role in Brassinosteroid signaling. Brassinosteroids play an essential role in various cellular processes, such as cell division, seed germination, vascular differentiation, flowering, xylem cell differentiation, stomata formation, photomorphogenesis, and pollen tube growth [76,77,78,79]. In addition, TaPERK2 and TaPERK4 can interact with four other wheat proteins (Traes_7AL_5E0DD589E.1, Traes_4DL_E447FD9FD.1, Traes_6BL_DFDCD5B11.1 and Traes_7DL_909EA97B3.1) which were cGMP-dependent protein kinase/PKG II and non-specific serine/threonine-protein kinase/Threonine-specific protein kinase. cGMP-dependent protein kinase is a phosphorylated diverse biologically important pathway [45,46,47]. PKG activated by cGMP has been implicated in cell division regulation, nucleic acid synthesis response to biotic stress, stomata closure during osmotic stress, and the development of adventitious roots [45,46,47,48,49]. TaPERK29 was highly elevated during the stripe rust infection, and cold stress might interact with Traes_6BL_DFDCD5B11.1, that is, cGMP-dependent protein kinase/PKG II specifically activates the signal transduction pathways implicated in different stress tolerance in wheat.
Moreover, TaPERK2 was co-expressed with TaPERK4 in spike development (Figure 8 and Figure S9), indicating that TaPERK2 and TaPERK4 might have an essential function in spike development through interacting with each other. These results provided valuable insight and the complex biological functions of TaPERK genes. In summary, this study provides valuable information about the TaPERK gene family, functions in plant growth, and response to phytohormones and different stress. Therefore, the outcome of this work is significant to dissect and understand the precise functions of TaPERK in developmental processes and various biotic and abiotic stress in wheat.

4. Materials and Methods

4.1. Identification of PERK Genes in Wheat

To carry out the genome-wide survey in bread wheat (Triticum aestivum) cv. Chinese Spring, genome data including genomic, CDS, and protein sequences of TaPERK genes were downloaded from the Ensembl plants biomart (http://plants.ensembl.org/biomart/martview, accessed on 27 September 2021). Two approaches were used to identify the PERK genes family in wheat. In the first method, we prepared a local database of the wheat protein sequences in BioEdit ver. 7.2.6 [80]. The sixty-two PERK genes from A. thaliana, G. max, O. sativa, and Z. mays were used for BLASTp against the local database. To find PERK genes, the e-value of 10−5 and >100-bit score were kept cut-off, and eventually, the BLASTp result was tabulated. In the second approach, the protein sequences of PERK from the above plant species were retrieved from the Ensembl Plants (http://plants.ensembl.org/index.html, accessed on 27 September 2021) and the BLASTp search was performed against the T. aestivum proteome with an e-value 10−5 and bit-score > 100. Based on the above method, putative PERK candidates were selected. Further putative PERK candidates were confirmed for the presence of protein tyrosine kinase domain using other online databases: InterPro (https://www.ebi.ac.uk/interpro, accessed on 25 September 2021), Simple Modular 132 Architecture Research Tool tool (SMART, http://smart.emblheidelberg.de/, accessed on 25 September 2021), HMMscan (https://www.ebi.ac.uk/Tools/hmmer/search/hmmscan, accessed on 25 September 2021) and NCBI CDD (https://www.ncbi.nlm.nih.gov/Structure/cdd/cdd.shtml, accessed on 25 September 2021). Finally, the protein sequences with protein tyrosine kinase domains were taken and renamed according to their chromosomal positions.

4.2. Genomic Localization, Gene Duplication, and Synteny Analysis

To map the chromosomal locations of TaPERK genes, genomic positions of PERK genes were downloaded from Ensembl plants biomart (http://plants.ensembl.org/biomart/martview, accessed on 26 September 2021). The PERK genes were named with a ‘Ta’ prefix and numbered according to their chromosomal positions. PhenoGram was used to map the TaPERK genes on the chromosomes (http://visualization.ritchielab.org/phenograms/plot, accessed on 26 September 2021). MCScanX tool kit was used to examine gene duplication events and synteny analysis within species and other plant species [81,82]. We used default parameters in MCScanx for synteny analysis. The non-synonymous (Ka) and synonymous substitution (Ks) ratio was calculated to estimate the selection pressure of duplicated TaPERK genes using the TBtools [82].

4.3. Biophysical Characteristics, Subcellular Localization, and 3D Structure

The biophysical characteristics of TaPERK proteins were predicted using ExPASy [83] and an isoelectric point calculator [84]. Subcellular localization was evaluated using CELLO [85], softberry (www.softberry.com, accessed on 27 September 2021), and BUSCA [86]. Finally, the three-dimensional structure of TaPERKs was generated using the Phyre2 web server [87].

4.4. Exon/intron Structure, Protein Motif, and Gene Ontology Analysis

The coding sequence, genomic and protein sequences of TaPERK genes were downloaded from the Ensembl plants biomart (http://plants.ensembl.org/biomart/martview, accessed on 27 September 2021). Exon, intron positions, and untranslated regions were elucidated using the Gene Structure Display Server 2.0 (http://gsds.gao-lab.org/, accessed on 27 September 2021). The protein motifs in the TaPERK were visualized using MEME (Multiple Em for Motif Elicitation ver.5.3.3; http://meme-suite.org/tools/meme, accessed on 27 September 2021) with default settings. TaPERK protein sequences were explored to detect GO terms enrichment using EggNOG (http://eggnogdb.embl.de/#/app/emapper, accessed on 27 September 2021) and agriGO [88].

4.5. Promoter Cis-Acting Regulatory Elements (CAREs) and Protein Interaction Network Analysis

To identify cis-elements, 2 kb upstream sequences of PERK genes were retrieved from Ensembl plants and examined using a PlantCARE online webserver (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 28 September 2021). The number of occurrences for each cis-element motif was counted for TaPERK genes, and the most frequently occurring CAREs were used to generate Figure 7 using TBtools [82]. The TaPERK protein interaction network was predicted using the STRING webserver (https://string-db.org/cgi, accessed on 28 September 2021).

4.6. Expression Analysis of TaPERK Genes

Transcripts per million (TPM) values for five different tissues, including leaf, stem, root, spike, grain, and under various stress conditions, were downloaded from the Wheat Expression database (http://www.wheat-expression.com/, accessed on 30 September 2021). Heatmaps and principal component analysis (PCA) were performed using ClustVis [89] and TBtools software [82].

4.7. Plant Growth Conditions, Stress Treatment, and RT-qPCR Analysis

Wheat (Triticum aestivum L.) cv. HI 1612 was used for the experiments. Seeds of HI 1612 were sown on soil in plastic pots and reared in a greenhouse. Ten-day-old wheat seedlings were acclimatized for two days in growth chamber conditions. They were further subjected to drought and high-temperature stress (40 °C) for 1 h and 6 h [90], and cold stress for 3 days (4 °C). For the combined drought and high-temperature stress, first, wheat seedling was exposed for the drought stress, then given a heat shock for 1 h and 6 h at 40 °C in an incubator. Controls were kept at 25 °C. The cold, drought, and high-temperature stressed seedlings were collected for RNA extraction and stored at −80 °C. The RNA was isolated from control, drought, and heat-stressed seedlings, as described by [91,92]. cDNA was synthesized using the iScriptTM cDNA synthesis kit (Bio-Rad, Hercules, CA, USA). Quantitative real-time PCR (RT-qPCR) was performed using the Applied Biosystems 7500 Fast Real-Time PCR (Applied Biosystems) with the SYBR Premix (Toyobo, Osaka, Japan). Wheat actin (AB181991) was used as a control to normalize the gene expression data. Transcript abundance was analyzed using the RT-qPCR. Each qRT-PCR reaction was carried out with three biological samples with two technical replicates and repeated three times. The fold change was calculated based on mean 2−ΔΔCT values and, eventually, this fold value was used to plot the graph [93,94]. Furthermore, one-way ANOVA, followed by Tukey’s HSD for multiple pairwise comparisons were applied. Means, standard errors and statistical significances for each sample were represented in figures (* p < 0.05, ** p < 0.01). All primers used in this study are mentioned in Table S10.

5. Conclusions

Wheat is the most important cereal crop and widely consumed staple food worldwide. However, global warming is becoming a severe threat to food security due to the constant climate changes, largely influencing plant development and productivity. This has raised a major challenge for plant biologists to increase yield and improve wheat’s quality, biotic and abiotic stress tolerance. The PERKS gene family plays a critical role in plant development and responses to various stresses. We identified and characterized the PERK gene family in wheat in this work. Expression patterns also revealed the role of TaPERKs in different developmental stages and stress conditions. Thus, this study facilitates a detailed understanding of PERK genes’ biological functions in wheat under different developmental processes and stress conditions.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants11040496/s1. Figure S1: Molecular weight (kDa) vs. isoelectric point plots of TaPERK genes. The distinct round shape colors represent the TaPERK gene family members. Figure S2: Distribution of TaPERKs in a different group of the phylogenetic tree. The y-axis indicates the number of TaPERK genes, and the x-axis indicates the phylogenetic groups. Figure S3: Phylogenetic analysis of TaPERK genes. A phylogenetic tree was constructed using MEGAX with the neighbor-joining (NJ) method and 1000 bootstrap replications. A black asterisk indicates the duplicated genes. Figure S4: Chromosomal distribution and duplicated PERK gene pairs in wheat. Duplicated PERK gene pairs are connected with lines with distinct colors. The figure was generated using TB tools. Figure S5: Distribution of exon and introns in TaPERKs gene family. The y-axis indicates the number of exons and introns, and the x-axis indicates the TaPERK genes. Exons and introns are represented by purple and orange boxes, respectively. Figure S6: Alignment and 3-dimensional structure of the TaPERK protein sequences. A. The conserved protein tyrosine kinase domain is boxed with red color. Colored and shaded amino acids are chemically similar residues. Dashes indicate gaps introduced to maximize the alignment of the homologous region. B. Predicted 3D structures TaPERK proteins. Figure S7: Gene ontology term distribution TaPERK gene family predicted using AgriGO A. Biological Process. B. Cellular component. C. Molecular function. Figure S8: PCA plots displaying grouping of different (A) developmental stages (B) biotic and abiotic stress conditions based on the TaPERK expression pattern. DS: drought stress, HS: heat stress, Zt: Zymoseptoria tritici, PM: powdery mildew; SR: stripe rust, h: hour and d: days. Figure S9: Heatmaps representing the expression pattern of TaPERK genes in different stress conditions. TPM values were directly used to construct the heatmaps. DS: drought stress, HS: heat stress, Zt: Zymoseptoria tritici, PM: powdery mildew; SR: stripe rust, h: hour and d: days, Table S1: TaPERK genomic, CDS protein and promoter sequence. Table S2: PERK proteins from Arabidopsis, rice, soybean, and wheat used to generate a phylogenetic tree. Table S3: Ratio of Ka/Ks and distribution of duplicate wheat PERK genes. Table S4: Orthologous relationships of TaPERK genes with other PERK genes in B. distachyon, Ae. tauschii, T. dicoccoides, O. sativa and A. thaliana. Table S5: Domain organization of TaPERK genes predicted using Pfam with default parameters. Table S6: cis-regulatory elements present in the TaPERK gene promoter region. Table S7: Significant Go term predicted in TaPERK gene family by AgriGo analysis. Table S8: TaPERK gene annotation using eggNOGmapper. Table S9: The protein-protein interaction network between TaPERK and other proteins in wheat. Table S10: qRT-PCR primers of TaPERK genes.

Author Contributions

M.K. and M.S.K. designed and wrote the manuscript; M.K. acquired funding; M.K. and M.S.K. supervised the study; B.S.K., A.S., P.D., S.R., C.M., D.S., C.R., K.H.M.S., A.K., R.G. and S.-M.C. provided 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 Dongguk University Research Fund of 2021.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is available in the manuscript and in the Supplementary Materials.

Acknowledgments

M.K. would like to acknowledge the Department of Life Science, Dongguk University, for providing infrastructure to carry out experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Day, E.K.; Sosale, N.G.; Lazzara, M.J. Cell signaling regulation by protein phosphorylation: A multivariate, heterogeneous, and context-dependent process. Curr. Opin. Biotechnol. 2016, 40, 185–192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Hohmann, U.; Lau, K.; Hothorn, M. The structural basis of ligand perception and signal activation by receptor kinases. Annu. Rev. Plant Biol. 2017, 68, 109–137. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Sim, J.-S.; Kesawat, M.S.; Kumar, M.; Kim, S.-Y.; Mani, V.; Subramanian, P.; Park, S.; Lee, C.-M.; Kim, S.-R.; Hahn, B.-S. Lack of the α1, 3-fucosyltransferase gene (OsFucT) affects anther development and pollen viability in rice. Int. J. Mol. Sci. 2018, 19, 1225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Hanks, S.K.; Quinn, A.M.; Hunter, T. The protein kinase family: Conserved features and deduced phylogeny of the catalytic domains. Science 1988, 241, 42–52. [Google Scholar] [CrossRef]
  5. DeFalco, T.A.; Anne, P.; James, S.R.; Willoughby, A.; Johanndrees, O.; Genolet, Y.; Pullen, A.-M.; Zipfel, C.; Hardtke, C.S.; Nimchuk, Z.L. A conserved regulatory module regulates receptor kinase signaling in immunity and development. bioRxiv 2021. [Google Scholar] [CrossRef]
  6. Nakhamchik, A.; Zhao, Z.; Provart, N.J.; Shiu, S.-H.; Keatley, S.K.; Cameron, R.K.; Goring, D.R. A comprehensive expression analysis of the Arabidopsis proline-rich extensin-like receptor kinase gene family using bioinformatic and experimental approaches. Plant Cell Physiol. 2004, 45, 1875–1881. [Google Scholar] [CrossRef] [Green Version]
  7. Qanmber, G.; Liu, J.; Yu, D.; Liu, Z.; Lu, L.; Mo, H.; Ma, S.; Wang, Z.; Yang, Z. Genome-wide identification and characterization of the PERK gene family in Gossypium hirsutum reveals gene duplication and functional divergence. Int. J. Mol. Sci. 2019, 20, 1750. [Google Scholar] [CrossRef] [Green Version]
  8. Diévart, A.; Clark, S.E. LRR-containing receptors regulating plant development and defense. Development 2004, 131, 251–261. [Google Scholar] [CrossRef] [Green Version]
  9. Shiu, S.-H.; Bleecker, A.B. Receptor-like kinases from Arabidopsis form a monophyletic gene family related to animal receptor kinases. Proc. Natl. Acad. Sci. USA 2001, 98, 10763–10768. [Google Scholar] [CrossRef] [Green Version]
  10. Morris, E.R.; Walker, J.C. Receptor-like protein kinases: The keys to response. Curr. Opin. Plant Biol. 2003, 6, 339–342. [Google Scholar] [CrossRef]
  11. Shiu, S.-H.; Karlowski, W.M.; Pan, R.; Tzeng, Y.-H.; Mayer, K.F.; Li, W.-H. Comparative analysis of the receptor-like kinase family in Arabidopsis and rice. Plant Cell 2004, 16, 1220–1234. [Google Scholar] [CrossRef] [Green Version]
  12. Silva, N.F.; Goring, D.R. The proline-rich, extensin-like receptor kinase-1 (PERK1) gene is rapidly induced by wounding. Plant Mol. Biol. 2002, 50, 667–685. [Google Scholar] [CrossRef]
  13. Haffani, Y.Z.; Silva, N.F.; Goring, D.R. Receptor kinase signalling in plants. Can. J. Bot. 2004, 82, 1–15. [Google Scholar] [CrossRef]
  14. Bai, L.; Zhang, G.; Zhou, Y.; Zhang, Z.; Wang, W.; Du, Y.; Wu, Z.; Song, C.P. Plasma membrane-associated proline-rich extensin-like receptor kinase 4, a novel regulator of Ca2+ signalling, is required for abscisic acid responses in Arabidopsis thaliana. Plant J. 2009, 60, 314–327. [Google Scholar] [CrossRef]
  15. Chen, G.; Wang, J.; Wang, H.; Wang, C.; Tang, X.; Li, J.; Zhang, L.; Song, J.; Hou, J.; Yuan, L. Genome-wide analysis of proline-rich extension-like receptor protein kinase (PERK) in Brassica rapa and its association with the pollen development. BMC Genom. 2020, 21, 401. [Google Scholar] [CrossRef]
  16. Borassi, C.; Sede, A.; Mecchia, M.A.; Mangano, S.; Marzol, E.; Denita-Juarez, S.P.; Salter, J.D.S.; Velasquez, S.M.; Muschietti, J.P.; Estevez, J. Proline-rich Extensin-like Receptor Kinases PERK5 and PERK12 are involved in Pollen Tube Growth. bioRxiv 2021. [Google Scholar] [CrossRef]
  17. Shiu, S.-H.; Bleecker, A.B. Expansion of the receptor-like kinase/Pelle gene family and receptor-like proteins in Arabidopsis. Plant Physiol. 2003, 132, 530–543. [Google Scholar] [CrossRef] [Green Version]
  18. Li, J.; Chory, J. A putative leucine-rich repeat receptor kinase involved in brassinosteroid signal transduction. Cell 1997, 90, 929–938. [Google Scholar] [CrossRef] [Green Version]
  19. Li, J.; Wen, J.; Lease, K.A.; Doke, J.T.; Tax, F.E.; Walker, J.C. BAK1, an Arabidopsis LRR receptor-like protein kinase, interacts with BRI1 and modulates brassinosteroid signaling. Cell 2002, 110, 213–222. [Google Scholar] [CrossRef] [Green Version]
  20. Nam, K.H.; Li, J. BRI1/BAK1, a receptor kinase pair mediating brassinosteroid signaling. Cell 2002, 110, 203–212. [Google Scholar] [CrossRef] [Green Version]
  21. Haffani, Y.; Silva-Gagliardi, N.; Sewter, S.; Aldea, M.G.; Zhao, Z.; Nakhamchik, A.; Cameron, R.; Goring, D. Altered expression of PERK receptor kinases in Arabidopsis leads to changes in growth and floral organ formation. Plant Signal. Behav. 2006, 1, 251–260. [Google Scholar] [CrossRef] [Green Version]
  22. Won, S.-K.; Lee, Y.-J.; Lee, H.-Y.; Heo, Y.-K.; Cho, M.; Cho, H.-T. Cis-element-and transcriptome-based screening of root hair-specific genes and their functional characterization in Arabidopsis. Plant Physiol. 2009, 150, 1459–1473. [Google Scholar] [CrossRef] [Green Version]
  23. Humphrey, T.V.; Haasen, K.E.; Aldea-Brydges, M.G.; Sun, H.; Zayed, Y.; Indriolo, E.; Goring, D.R. PERK–KIPK–KCBP signalling negatively regulates root growth in Arabidopsis thaliana. J. Exp. Bot. 2015, 66, 71–83. [Google Scholar] [CrossRef] [Green Version]
  24. Hwang, I.; Kim, S.Y.; Kim, C.S.; Park, Y.; Tripathi, G.R.; Kim, S.-K.; Cheong, H. Over-expression of the IGI1 leading to altered shoot-branching development related to MAX pathway in Arabidopsis. Plant Mol. Biol. 2010, 73, 629–641. [Google Scholar] [CrossRef] [Green Version]
  25. Borassi, C.; Sede, A.R.; Mecchia, M.A.; Salgado Salter, J.D.; Marzol, E.; Muschietti, J.P.; Estevez, J.M. An update on cell surface proteins containing extensin-motifs. J. Exp. Bot. 2016, 67, 477–487. [Google Scholar] [CrossRef] [Green Version]
  26. Mansfield, T.; Hetherington, A.; Atkinson, C. Some current aspects of stomatal physiology. Annu. Rev. Plant Biol. 1990, 41, 55–75. [Google Scholar] [CrossRef]
  27. Webb, A.A.; McAinsh, M.R.; Taylor, J.E.; Hetherington, A.M. Calcium ions as intracellular second messengers in higher plants. Adv. Bot. Res. 1996, 22, 45–96. [Google Scholar] [CrossRef]
  28. Gong, M.; Chen, S.-N.; Song, Y.-Q.; Li, Z.-G. Effect of calcium and calmodulin on intrinsic heat tolerance in relation to antioxidant systems in maize seedlings. Funct. Plant Biol. 1997, 24, 371–379. [Google Scholar] [CrossRef]
  29. Hwang, Y.; Lee, H.; Lee, Y.S.; Cho, H.T. Cell wall-associated ROOT HAIR SPECIFIC 10, a proline-rich receptor-like kinase, is a negative modulator of Arabidopsis root hair growth. J. Exp. Bot. 2016, 67, 2007–2022. [Google Scholar] [CrossRef] [Green Version]
  30. Samuel, M.A.; Ellis, B.E. Double jeopardy: Both overexpression and suppression of a redox-activated plant mitogen-activated protein kinase render tobacco plants ozone sensitive. Plant Cell 2002, 14, 2059–2069. [Google Scholar] [CrossRef] [Green Version]
  31. Xing, Y.; Cao, Q.Q.; Zhang, Q.; Qin, L.; Jia, W.S.; Zhang, J.H. MKK5 Regulates High Light-Induced Gene Expression of Cu/Zn Superoxide Dismutase 1 and 2 in Arabidopsis. Plant Cell Physiol. 2013, 54, 1217–1227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Florentino, L.H.; Santos, A.A.; Fontenelle, M.R.; Pinheiro, G.L.; Zerbini, F.M.; Baracat-Pereira, M.C.; Fontes, E.P. A PERK-like receptor kinase interacts with the geminivirus nuclear shuttle protein and potentiates viral infection. J. Virol. 2006, 80, 6648–6656. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Singh, K.M.; Kumar, D.B.; Kumar, D.S. Manorama. Assessment of genetic diversity among Indian Sesame (Sesamum indicum L.) accessions using RAPD, ISSR and SSR markers. Res. J. Biotechnol. 2015, 10, 35–47. [Google Scholar]
  34. Kumar, M.; Kesawat, M.S.; Ali, A.; Lee, S.-C.; Gill, S.S.; Kim, H.U. Integration of abscisic acid signaling with other signaling pathways in plant stress responses and development. Plants 2019, 8, 592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kesawat, M.S.; Shivaraj, S.; Kim, D.K.; Kumar, M.; Hahn, B.S.; Deshmukh, R. Metalloids and Their Role in the Biological System. Met. Plants Adv. Future Prospect. 2020, 1–17. [Google Scholar] [CrossRef]
  36. The International Wheat Genome Sequencing Consortium. A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome. Science 2014, 345, 1251788. [Google Scholar] [CrossRef]
  37. Gill, B.S.; Appels, R.; Botha-Oberholster, A.-M.; Buell, C.R.; Bennetzen, J.L.; Chalhoub, B.; Chumley, F.; Dvorák, J.; Iwanaga, M.; Keller, B. A workshop report on wheat genome sequencing: International Genome Research on Wheat Consortium. Genetics 2004, 168, 1087–1096. [Google Scholar] [CrossRef] [Green Version]
  38. Afzal, F.; Chaudhari, S.K.; Gul, A.; Farooq, A.; Ali, H.; Nisar, S.; Sarfraz, B.; Shehzadi, K.J.; Mujeeb-Kazi, A. Bread wheat (Triticum aestivum L.) under biotic and abiotic stresses: An overview. Crop Prod. Glob. Environ. Issues 2015, 293–317. [Google Scholar] [CrossRef]
  39. Kumar, P.; Yadava, R.; Gollen, B.; Kumar, S.; Verma, R.K.; Yadav, S. Nutritional contents and medicinal properties of wheat: A review. Life Sci. Med. Res. 2011, 22, 1–10. [Google Scholar]
  40. Kumar, M.; Kherawat, B.S.; Dey, P.; Saha, D.; Singh, A.; Bhatia, S.K.; Ghodake, G.S.; Kadam, A.A.; Kim, H.-U.; Chung, S.-M. Genome-Wide Identification and Characterization of PIN-FORMED (PIN) Gene Family Reveals Role in Developmental and Various Stress Conditions in Triticum aestivum L. Int. J. Mol. Sci. 2021, 22, 7396. [Google Scholar] [CrossRef]
  41. Kesawat, M.S.; Das, B.K.; Bhaganagare, G.R. Genome-wide identification, evolutionary and expression analyses of putative Fe–S biogenesis genes in rice (Oryza sativa). Genome 2012, 55, 571–583. [Google Scholar] [CrossRef]
  42. Kesawat, M.S.; Kherawat, B.S.; Singh, A.; Dey, P.; Kabi, M.; Debnath, D.; Saha, D.; Khandual, A.; Rout, S.; Ali, A. Genome-wide identification and characterization of the brassinazole-resistant (BZR) gene family and its expression in the various developmental stage and stress conditions in wheat (Triticum aestivum L.). Int. J. Mol. Sci. 2021, 22, 8743. [Google Scholar] [CrossRef]
  43. Appels, R.; Eversole, K.; Stein, N.; Feuillet, C.; Keller, B.; Rogers, J.; Pozniak, C.J.; Choulet, F.; Distelfeld, A.; Poland, J. Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 2018, 361. [Google Scholar] [CrossRef] [Green Version]
  44. Lawton-Rauh, A. Evolutionary dynamics of duplicated genes in plants. Mol. Phylogenet. Evol. 2003, 29, 396–409. [Google Scholar] [CrossRef]
  45. Hou, Y.; Gupta, N.; Schoenlein, P.; Wong, E.; Martindale, R.; Ganapathy, V.; Browning, D. An anti-tumor role for cGMP-dependent protein kinase. Cancer Lett. 2006, 240, 60–68. [Google Scholar] [CrossRef]
  46. Gross, I.; Durner, J. In search of enzymes with a role in 3′, 5′-cyclic guanosine monophosphate metabolism in plants. Front. Plant Sci. 2016, 7, 576. [Google Scholar] [CrossRef] [Green Version]
  47. Shen, Q.; Zhan, X.; Yang, P.; Li, J.; Chen, J.; Tang, B.; Wang, X.; Hong, Y. Dual activities of plant cGMP-dependent protein kinase and its roles in gibberellin signaling and salt stress. Plant Cell 2019, 31, 3073–3091. [Google Scholar] [CrossRef] [Green Version]
  48. Wheeler, J.I.; Wong, A.; Marondedze, C.; Groen, A.J.; Kwezi, L.; Freihat, L.; Vyas, J.; Raji, M.A.; Irving, H.R.; Gehring, C. The brassinosteroid receptor BRI 1 can generate cGMP enabling cGMP-dependent downstream signaling. Plant J. 2017, 91, 590–600. [Google Scholar] [CrossRef] [Green Version]
  49. Isner, J.-C.; Maathuis, F.J. cGMP signalling in plants: From enigma to main stream. Funct. Plant Biol. 2016, 45, 93–101. [Google Scholar] [CrossRef]
  50. Lespinet, O.; Wolf, Y.I.; Koonin, E.V.; Aravind, L. The role of lineage-specific gene family expansion in the evolution of eukaryotes. Genome Res. 2002, 12, 1048–1059. [Google Scholar] [CrossRef] [Green Version]
  51. Jordan, I.K.; Makarova, K.S.; Spouge, J.L.; Wolf, Y.I.; Koonin, E.V. Lineage-specific gene expansions in bacterial and archaeal genomes. Genome Res. 2001, 11, 555–565. [Google Scholar] [CrossRef]
  52. Ramsey, J.; Schemske, D.W. Pathways, mechanisms, and rates of polyploid formation in flowering plants. Annu. Rev. Ecol. Syst. 1998, 29, 467–501. [Google Scholar] [CrossRef] [Green Version]
  53. Moore, R.C.; Purugganan, M.D. The early stages of duplicate gene evolution. Proc. Natl. Acad. Sci. USA 2003, 100, 15682–15687. [Google Scholar] [CrossRef] [Green Version]
  54. 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] [Green Version]
  55. Baumberger, N.; Doesseger, B.; Guyot, R.; Diet, A.; Parsons, R.L.; Clark, M.A.; Simmons, M.; Bedinger, P.; Goff, S.A.; Ringli, C. Whole-genome comparison of leucine-rich repeat extensins in Arabidopsis and rice. A conserved family of cell wall proteins form a vegetative and a reproductive clade. Plant Physiol. 2003, 131, 1313–1326. [Google Scholar] [CrossRef] [Green Version]
  56. Wang, D.; Guo, Y.; Wu, C.; Yang, G.; Li, Y.; Zheng, C. Genome-wide analysis of CCCH zinc finger family in Arabidopsis and rice. BMC Genom. 2008, 9, 44. [Google Scholar] [CrossRef] [Green Version]
  57. Yang, Z.; Gong, Q.; Qin, W.; Yang, Z.; Cheng, Y.; Lu, L.; Ge, X.; Zhang, C.; Wu, Z.; Li, F. Genome-wide analysis of WOX genes in upland cotton and their expression pattern under different stresses. BMC Plant Biol. 2017, 17, 113. [Google Scholar] [CrossRef]
  58. Yang, Z.; Gong, Q.; Wang, L.; Jin, Y.; Xi, J.; Li, Z.; Qin, W.; Yang, Z.; Lu, L.; Chen, Q. Genome-wide study of YABBY genes in upland cotton and their expression patterns under different stresses. Front. Genet. 2018, 9, 33. [Google Scholar] [CrossRef] [Green Version]
  59. Zhang, B.; Liu, J.; Yang, Z.E.; Chen, E.Y.; Zhang, C.J.; Zhang, X.Y.; Li, F.G. Genome-wide analysis of GRAS transcription factor gene family in Gossypium hirsutum L. BMC Genom. 2018, 19, 348. [Google Scholar] [CrossRef]
  60. Yin, G.; Xu, H.; Xiao, S.; Qin, Y.; Li, Y.; Yan, Y.; Hu, Y. The large soybean (Glycine max) WRKY TF family expanded by segmental duplication events and subsequent divergent selection among subgroups. BMC Plant Biol. 2013, 13, 148. [Google Scholar] [CrossRef] [Green Version]
  61. Ren, Z.; Yu, D.; Yang, Z.; Li, C.; Qanmber, G.; Li, Y.; Li, J.; Liu, Z.; Lu, L.; Wang, L. Genome-wide identification of the MIKC-type MADS-box gene family in Gossypium hirsutum L. unravels their roles in flowering. Front. Plant Sci. 2017, 8, 384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Dossa, K.; Diouf, D.; Cissé, N. Genome-wide investigation of Hsf genes in sesame reveals their segmental duplication expansion and their active role in drought stress response. Front. Plant Sci. 2016, 7, 1522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Roy, S.W.; Gilbert, W. The evolution of spliceosomal introns: Patterns, puzzles and progress. Nat. Rev. Genet. 2006, 7, 211–221. [Google Scholar] [PubMed]
  64. Roy, S.W.; Penny, D. Patterns of intron loss and gain in plants: Intron loss–dominated evolution and genome-wide comparison of O. sativa and A. thaliana. Mol. Biol. Evol. 2007, 24, 171–181. [Google Scholar] [CrossRef] [Green Version]
  65. Serrano, M.; Parra, S.; Alcaraz, L.D.; Guzmán, P. The ATL gene family from Arabidopsis thaliana and Oryza sativa comprises a large number of putative ubiquitin ligases of the RING-H2 type. J. Mol. Evol. 2006, 62, 434–445. [Google Scholar] [CrossRef]
  66. Dan, Y.; Niu, Y.; Wang, C.; Yan, M.; Liao, W. Genome-wide identification and expression analysis of the trehalose-6-phosphate synthase (TPS) gene family in cucumber (Cucumis sativus L.). PeerJ 2021, 9, e11398. [Google Scholar] [CrossRef]
  67. Lecharny, A.; Boudet, N.; Gy, I.; Aubourg, S.; Kreis, M. Introns in, introns out in plant gene families: A genomic approach of the dynamics of gene structure. J. Struct. Funct. Genom. 2003, 3, 111–116. [Google Scholar] [CrossRef]
  68. Hernandez-Garcia, C.M.; Finer, J.J. Identification and validation of promoters and cis-acting regulatory elements. Plant Sci. 2014, 217, 109–119. [Google Scholar] [CrossRef] [Green Version]
  69. Roy, A.L.; Sen, R.; Roeder, R.G. Enhancer–promoter communication and transcriptional regulation of Igh. Trends Immunol. 2011, 32, 532–539. [Google Scholar] [CrossRef] [Green Version]
  70. Roy, A.L.; Singer, D.S. Core promoters in transcription: Old problem, new insights. Trends Biochem. Sci. 2015, 40, 165–171. [Google Scholar] [CrossRef] [Green Version]
  71. Zhang, Y.; Wong, C.-H.; Birnbaum, R.Y.; Li, G.; Favaro, R.; Ngan, C.Y.; Lim, J.; Tai, E.; Poh, H.M.; Wong, E. Chromatin connectivity maps reveal dynamic promoter–enhancer long-range associations. Nature 2013, 504, 306–310. [Google Scholar] [CrossRef] [Green Version]
  72. Cheng, Y.; Tang, Q.; Li, Y.; Zhang, Y.; Zhao, C.; Yan, J.; You, H. Folding/unfolding kinetics of G-quadruplexes upstream of the P1 promoter of the human BCL-2 oncogene. J. Biol. Chem. 2019, 294, 5890–5895. [Google Scholar] [CrossRef]
  73. Fankhauser, C.; Chory, J. Light control of plant development. Annu. Rev. Cell Dev. Biol. 1997, 13, 203–229. [Google Scholar] [CrossRef] [Green Version]
  74. Dievart, A.; Gottin, C.; Périn, C.; Ranwez, V.; Chantret, N. Origin and diversity of plant receptor-like kinases. Annu. Rev. Plant Biol. 2020, 71, 131–156. [Google Scholar] [CrossRef] [Green Version]
  75. Li, S.; Liu, Z.; Chen, G.; Qanmber, G.; Lu, L.; Zhang, J.; Ma, S.; Yang, Z.; Li, F. Identification and analysis of GhEXO gene family indicated that GhEXO7_At promotes plant growth and development through brassinosteroid signaling in cotton (Gossypium hirsutum L.). Front. Plant Sci. 2021, 12, 719889. [Google Scholar] [CrossRef]
  76. Šimura, J.; Antoniadi, I.; Široká, J.; Tarkowská, D.e.; Strnad, M.; Ljung, K.; Novák, O. Plant hormonomics: Multiple phytohormone profiling by targeted metabolomics. Plant Physiol. 2018, 177, 476–489. [Google Scholar] [CrossRef] [Green Version]
  77. Clouse, S.D. Brassinosteroid signal transduction: From receptor kinase activation to transcriptional networks regulating plant development. Plant Cell 2011, 23, 1219–1230. [Google Scholar] [CrossRef] [Green Version]
  78. Clouse, S.D.; Langford, M.; McMorris, T.C. A brassinosteroid-insensitive mutant in Arabidopsis thaliana exhibits multiple defects in growth and development. Plant Physiol. 1996, 111, 671–678. [Google Scholar] [CrossRef] [Green Version]
  79. Kim, E.-J.; Russinova, E. Brassinosteroid signalling. Curr. Biol. 2020, 30, R294–R298. [Google Scholar] [CrossRef]
  80. Hall, T.A. BioEdit: A User-Friendly Biological Sequence Alignment Editor and Analysis Program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  81. 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] [Green Version]
  82. 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] [PubMed]
  83. Gasteiger, E.; Gattiker, A.; Hoogland, C.; Ivanyi, I.; Appel, R.D.; Bairoch, A. ExPASy: The proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003, 31, 3784–3788. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Kozlowski, L.P. IPC–isoelectric point calculator. Biol. Direct 2016, 11, 1–16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Yu, C.S.; Chen, Y.C.; Lu, C.H.; Hwang, J.K. Prediction of protein subcellular localization. Proteins Struct. Funct. Bioinform. 2006, 64, 643–651. [Google Scholar] [CrossRef] [PubMed]
  86. Savojardo, C.; Martelli, P.L.; Fariselli, P.; Profiti, G.; Casadio, R. BUSCA: An integrative web server to predict subcellular localization of proteins. Nucleic Acids Res. 2018, 46, W459–W466. [Google Scholar] [CrossRef] [PubMed]
  87. Kelley, L.A.; Mezulis, S.; Yates, C.M.; Wass, M.N.; Sternberg, M.J. The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc. 2015, 10, 845–858. [Google Scholar] [CrossRef] [Green Version]
  88. Tian, T.; Liu, Y.; Yan, H.; You, Q.; Yi, X.; Du, Z.; Xu, W.; Su, Z. agriGO v2. 0: A GO analysis toolkit for the agricultural community, 2017 update. Nucleic Acids Res. 2017, 45, W122–W129. [Google Scholar] [CrossRef]
  89. Metsalu, T.; Vilo, J. ClustVis: A web tool for visualizing clustering of multivariate data using Principal Component Analysis and heatmap. Nucleic Acids Res. 2015, 43, W566–W570. [Google Scholar] [CrossRef]
  90. Grigorova, B.; Vaseva, I.; Demirevska, K.; Feller, U. Combined Drought and Heat Stress in Wheat: Changes in Some Heat Shock Proteins. Biol. Plant. 2011, 55, 105–111. [Google Scholar] [CrossRef]
  91. Kim, D.K.; Kesawat, M.S.; Hong, C.B. One gene member of the ADP-ribosylation factor family is heat-inducible and enhances seed germination in Nicotiana tabacum. Genes Genom. 2017, 39, 1353–1365. [Google Scholar] [CrossRef]
  92. Kesawat, M.S.; Kim, D.K.; Zeba, N.; Suh, M.C.; Xia, X.; Hong, C.B. Ectopic RING zinc finger gene from hot pepper induces totally different genes in lettuce and tobacco. Mol. Breed. 2018, 38, 1–24. [Google Scholar] [CrossRef] [Green Version]
  93. 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]
  94. Narancio, R.; John, U.; Mason, J.; Spangenberg, G. Selection of optimal reference genes for quantitative RT-PCR transcript abundance analysis in white clover (Trifolium repens L.). Funct. Plant Biol. 2018, 45, 737–744. [Google Scholar] [CrossRef]
Plants 11 00496 g001 550
Figure 1. Phylogenetic analysis of TaPERK proteins with Arabidopsis (15), rice (8), and soybean (16). The phylogenetic analysis was executed using the ClustalW program as well as MEGAX software by the neighbor-joining method and bootstrap values of 1000 replicates. The numbers on the nodes indicate the bootstrap values. Distinct groups are represented by the different colors.
Figure 1. Phylogenetic analysis of TaPERK proteins with Arabidopsis (15), rice (8), and soybean (16). The phylogenetic analysis was executed using the ClustalW program as well as MEGAX software by the neighbor-joining method and bootstrap values of 1000 replicates. The numbers on the nodes indicate the bootstrap values. Distinct groups are represented by the different colors.
Plants 11 00496 g001
Plants 11 00496 g002a 550Plants 11 00496 g002b 550
Figure 2. Genomic distribution of identified PERK genes on the 21 chromosomes of wheat and within the three sub-genomes. (A) Schematic representations of the chromosomal distribution of PERK genes on the 21 chromosomes of wheat and the name of the gene on the right. The colored circles on the chromosomes indicate the position of the PERK genes. The chromosome numbers of the three sub-genomes are indicated at the top of each bar. (B) Distribution of PERK genes in the three sub-genomes. (C) Distribution of PERK genes across 21 chromosomes, Un: unaligned contig.
Figure 2. Genomic distribution of identified PERK genes on the 21 chromosomes of wheat and within the three sub-genomes. (A) Schematic representations of the chromosomal distribution of PERK genes on the 21 chromosomes of wheat and the name of the gene on the right. The colored circles on the chromosomes indicate the position of the PERK genes. The chromosome numbers of the three sub-genomes are indicated at the top of each bar. (B) Distribution of PERK genes in the three sub-genomes. (C) Distribution of PERK genes across 21 chromosomes, Un: unaligned contig.
Plants 11 00496 g002aPlants 11 00496 g002b
Plants 11 00496 g003 550
Figure 3. Syntenic relationships of TaPERK genes between Aegilops tauschii, Brachypodium distachyon, and Oryza sativa. The gray lines in the background represent the collinear blocks within Triticum aestivum and other plant genomes, while the red lines highlight the syntenic PERK gene pairs.
Figure 3. Syntenic relationships of TaPERK genes between Aegilops tauschii, Brachypodium distachyon, and Oryza sativa. The gray lines in the background represent the collinear blocks within Triticum aestivum and other plant genomes, while the red lines highlight the syntenic PERK gene pairs.
Plants 11 00496 g003
Plants 11 00496 g004 550
Figure 4. Diagrammatic representation of the exon–intron organization of the TaPERK genes. Yellow boxes represent exons, untranslated regions (UTRs) are indicated by blue boxes, and black lines represent introns. The lengths of the boxes and lines are scaled based on gene length. The exon and intron sizes can be estimated using the scale at the bottom.
Figure 4. Diagrammatic representation of the exon–intron organization of the TaPERK genes. Yellow boxes represent exons, untranslated regions (UTRs) are indicated by blue boxes, and black lines represent introns. The lengths of the boxes and lines are scaled based on gene length. The exon and intron sizes can be estimated using the scale at the bottom.
Plants 11 00496 g004
Plants 11 00496 g005a 550Plants 11 00496 g005b 550
Figure 5. Conserved motifs of TaPERK genes elucidated by MEME. Up to 10 motifs were shown in different colors. (A) Colored boxes representing different conserved motifs with different sequences and sizes. (B) Sequence logo conserved motif of the wheat PERK proteins. The overall height of each stack represents the degree of conservation at this position, while the height of individual letters within each stack indicates the relative frequency of the corresponding amino acids. The sequence of each motif, combined p-value, and length are shown on the left side of the figure. MEME Parameters: number of repetitions, any; maximum number of motifs, 10; optimum motif width, between 6 and 50.
Figure 5. Conserved motifs of TaPERK genes elucidated by MEME. Up to 10 motifs were shown in different colors. (A) Colored boxes representing different conserved motifs with different sequences and sizes. (B) Sequence logo conserved motif of the wheat PERK proteins. The overall height of each stack represents the degree of conservation at this position, while the height of individual letters within each stack indicates the relative frequency of the corresponding amino acids. The sequence of each motif, combined p-value, and length are shown on the left side of the figure. MEME Parameters: number of repetitions, any; maximum number of motifs, 10; optimum motif width, between 6 and 50.
Plants 11 00496 g005aPlants 11 00496 g005b
Plants 11 00496 g006a 550Plants 11 00496 g006b 550
Figure 6. Multiple sequence alignment of the TaPERK protein sequences. The conserved protein tyrosine kinase domain is boxed in red. Colored and shaded amino acids are chemically similar residues. Dashes indicate gaps introduced to maximize the alignment of the homologous region. * indicates positions which have a single, fully conserved residue.
Figure 6. Multiple sequence alignment of the TaPERK protein sequences. The conserved protein tyrosine kinase domain is boxed in red. Colored and shaded amino acids are chemically similar residues. Dashes indicate gaps introduced to maximize the alignment of the homologous region. * indicates positions which have a single, fully conserved residue.
Plants 11 00496 g006aPlants 11 00496 g006b
Plants 11 00496 g007 550
Figure 7. Cis-acting regulatory elements (CAREs) in the promoter region of the TaPERK genes family. The CAREs analysis was performed with a 2kb upstream region using PlantCARE online server. The different numbers of cis-regulatory elements represent different colors. (A) Hormone-responsive elements, stress-responsive elements, growth and development-related elements, light-responsive elements, and other elements with unknown functions are differentiated by color. (B) Most commonly occurring CAREs in TaPERKs.
Figure 7. Cis-acting regulatory elements (CAREs) in the promoter region of the TaPERK genes family. The CAREs analysis was performed with a 2kb upstream region using PlantCARE online server. The different numbers of cis-regulatory elements represent different colors. (A) Hormone-responsive elements, stress-responsive elements, growth and development-related elements, light-responsive elements, and other elements with unknown functions are differentiated by color. (B) Most commonly occurring CAREs in TaPERKs.
Plants 11 00496 g007
Plants 11 00496 g008 550
Figure 8. Heatmap representing expression profile of the TaPERK genes at various developmental stages. Columns represent genes, and rows represent different developmental stages. TPM values were used directly to create the heatmaps. The “z” nomenclature refers to Zadok’s growth stage.
Figure 8. Heatmap representing expression profile of the TaPERK genes at various developmental stages. Columns represent genes, and rows represent different developmental stages. TPM values were used directly to create the heatmaps. The “z” nomenclature refers to Zadok’s growth stage.
Plants 11 00496 g008
Plants 11 00496 g009 550
Figure 9. Quantitative real-time PCR analysis of selected TaPERK genes in response to drought stress (DS), heat stress (HS), and cold stress to verify RNA seq data. The wheat actin gene was used as the internal control to standardize the RNA samples for each reaction. Asterisks indicate significant differences compared with control. Bars represent results of Tukey’s HSD test at the <0.05 and <0.001 level (* p < 0.05, ** p lies in between the values of 0.05 and 0.001, and *** p < 0.001). Error bars show standard deviation. Data are mean ± SD (n = 3).
Figure 9. Quantitative real-time PCR analysis of selected TaPERK genes in response to drought stress (DS), heat stress (HS), and cold stress to verify RNA seq data. The wheat actin gene was used as the internal control to standardize the RNA samples for each reaction. Asterisks indicate significant differences compared with control. Bars represent results of Tukey’s HSD test at the <0.05 and <0.001 level (* p < 0.05, ** p lies in between the values of 0.05 and 0.001, and *** p < 0.001). Error bars show standard deviation. Data are mean ± SD (n = 3).
Plants 11 00496 g009
Plants 11 00496 g010 550
Figure 10. Protein–protein interaction analysis of TaPERKs proteins. Protein–protein interaction network produced by STRINGV9.1, each node represents a protein, and each edge represents an interaction, colored by evidence type. The figure highlights the connections between differentially represented proteins.
Figure 10. Protein–protein interaction analysis of TaPERKs proteins. Protein–protein interaction network produced by STRINGV9.1, each node represents a protein, and each edge represents an interaction, colored by evidence type. The figure highlights the connections between differentially represented proteins.
Plants 11 00496 g010
Table
Table 1. Nomenclature and characteristics of the putative proline-rich extensin-like receptor kinases (PERKs) proteins in wheat were predicted using various computational tools.
Table 1. Nomenclature and characteristics of the putative proline-rich extensin-like receptor kinases (PERKs) proteins in wheat were predicted using various computational tools.
Proposed Gene NameGene IDGenomic LocationOrientationCDS Length (bp)Intron NumberProtein Length (aa)Molecular Weight (KDa)Isoelectric Point (pI)GRAVYPredicted Subcellular Localization
TaPERK1TraesCS1A02G1279001A:155693812–155696618Forward1977765869.447.53−0.531Nucleus
TaPERK2TraesCS1B02G14700001B:209130189–209130266Reverse1431847652.146.17−0.5Nucleus
TaPERK3TraesCS1D02G004301D:2110107–2112027Forward1971765668.939.04−0.393Chloroplast outer
membrane
TaPERK4TraesCS1D02G1263001D:137437684–137440387Reverse1962765369.077.21−0.52Nucleus
TaPERK5TraesCS2A02G4182002A:674030843–674031911Forward3048231015110.336.33−0.193Plasma membrane
TaPERK6TraesCS2A02G4183002A:674050369–674051442Forward3042231013110.397.06−0.135Plasma membrane
TaPERK7TraesCS2A02G4184002A:674061248–674062244Forward3159231052113.735.96−0.127Plasma membrane
TaPERK8TraesCS2B02G4372002B:629023953–629025021Forward3045231014110.386.3−0.178Plasma membrane
TaPERK9TraesCS2B02G4373002B:629106216–629107285Forward3048231015110.626.61−0.147Plasma membrane
TaPERK10TraesCS2D02G4156002D:529537635–529538701Forward3048221015110.346.6−0.167Plasma membrane
TaPERK11TraesCS2D02G4157002D:529548057–529548998Forward277523924101.077.29−0.182Plasma membrane
TaPERK12TraesCS2D02G4158002D:529558487–529559547Forward3156231051113.456.11−0.102Plasma membrane
TaPERK13TraesCS3A02G0039003A:1925607–1927275Reverse2064768772.425.96−0.429Plasma membrane
TaPERK14TraesCS3A02G1522003A:142891955–142894634Forward1893763067.436.28−0.569Endomembrane
system
TaPERK15TraesCS3A02G2298003A:429615911–429617422Reverse2163672074.977.93−0.401Chloroplast
thylakoid lumen
TaPERK16TraesCS3A02G2781003A:507637093–507638935Reverse2028767572.427.31−0.481Plasma membrane
TaPERK17TraesCS3A02G2903003A:519244808–519246110Reverse2184772775.86.11−0.535Endomembrane
system
TaPERK18TraesCS3B02G0086003B:4324660–4326408Forward2061768671.885.97−0.437Plasma membrane
TaPERK19TraesCS3B02G1793003B:187347873–187350697Forward1896763167.466.35−0.569Endomembrane
system
TaPERK20TraesCS3B02G2591003B:416806224–416809608Reverse2097669872.987.63−0.448Plasma membrane
TaPERK21TraesCS3B02G3123003B:501498044–501499926Reverse2034767772.67.31−0.493Endomembrane
system
TaPERK22TraesCS3B02G3251003B:525990462–525991846Reverse2436781184.776.09−0.51Plasma membrane
TaPERK23TraesCS3D02G0054003D:2141185–2143272Forward1206640144.375.55−0.403Nucleus
TaPERK24TraesCS3D02G1600003D:130928685–130931461Forward1899763267.496.36−0.555Endomembrane
system
TaPERK25TraesCS3D02G2784003D:385473929–385474240Reverse2031867672.567.1−0.471Endomembrane
system
TaPERK26TraesCS3D02G2901003D:400311470–400312883Reverse1317643847.145.97−0.447Nucleus
TaPERK27TraesCS4A02G0775004A:76627667–76628358Forward1866562164.495.58−0.457Endomembrane
system
TaPERK28TraesCS4A02G4497004A:715718345–715719349Forward26041986794.757.03−0.145Plasma membrane
TaPERK29TraesCS4B02G2336004B:486206279–486206961Reverse1857561864.455.63−0.487Plasma membrane
TaPERK30TraesCS5A02G4113005A:599978835–599979642Reverse1722457360.337.96−0.387Chloroplast outer
membrane
TaPERK31TraesCS5B02G4150005B:589228532–589228944Reverse1842361364.687.86−0.368Chloroplast outer
membrane
TaPERK32TraesCS7A02G0386007A:17358644–17359648Reverse3030231009109.76.22−0.108Plasma membrane
TaPERK33TraesCS7A02G2319007A:202852283–202853761Reverse2187672876.245.32−0.484Nucleus
TaPERK34TraesCS7B02G1304007B:156752944–156754400Reverse2208673576.895.22−0.49Nucleus
TaPERK35TraesCS7D02G0348007D:17864178–17865182Reverse3021231006109.266.1−0.09Plasma membrane
TaPERK36TraesCS7D02G2327007D:194224547–194225929Forward2256675178.26.16−0.645Endomembrane
system
TaPERK37TraesCSU02G104700Un:92294980–92296477Reverse2205673476.555.32−0.483Nucleus
ID: identity; bp: base pair; aa: amino acids; pI: isoelectric point; MW: molecular weight; KDa: Kilo dalton.
Table
Table 2. Number of PERK proteins in different plant species.
Table 2. Number of PERK proteins in different plant species.
Plant SpeciesGenome Size (Approx.)Coding GenesPERK Genes
Triticum aestivum (6n)17 Gb107,89137
Arabidopsis thaliana (2n)135 Mb27,65515
Oryza sativa500 Mb37,9608
Zea mays (2n)2.4 Gb39,59123
Glycine max (2n)1.15 Gb55,89716
Sorghum bicolor (2n)730 Mb28,12015
Gossypium arboretum (2n)1746 Mb41,33015
Gossypium raimondii (2n)885 Mb40,97616
Gossypium hirsutum (4n)2.43 Gb75,37633
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kesawat, M.S.; Kherawat, B.S.; Singh, A.; Dey, P.; Routray, S.; Mohapatra, C.; Saha, D.; Ram, C.; Siddique, K.H.M.; Kumar, A.; et al. Genome-Wide Analysis and Characterization of the Proline-Rich Extensin-like Receptor Kinases (PERKs) Gene Family Reveals Their Role in Different Developmental Stages and Stress Conditions in Wheat (Triticum aestivum L.). Plants 2022, 11, 496. https://doi.org/10.3390/plants11040496

AMA Style

Kesawat MS, Kherawat BS, Singh A, Dey P, Routray S, Mohapatra C, Saha D, Ram C, Siddique KHM, Kumar A, et al. Genome-Wide Analysis and Characterization of the Proline-Rich Extensin-like Receptor Kinases (PERKs) Gene Family Reveals Their Role in Different Developmental Stages and Stress Conditions in Wheat (Triticum aestivum L.). Plants. 2022; 11(4):496. https://doi.org/10.3390/plants11040496

Chicago/Turabian Style

Kesawat, Mahipal Singh, Bhagwat Singh Kherawat, Anupama Singh, Prajjal Dey, Snehasish Routray, Chinmayee Mohapatra, Debanjana Saha, Chet Ram, Kadambot H. M. Siddique, Ajay Kumar, and et al. 2022. "Genome-Wide Analysis and Characterization of the Proline-Rich Extensin-like Receptor Kinases (PERKs) Gene Family Reveals Their Role in Different Developmental Stages and Stress Conditions in Wheat (Triticum aestivum L.)" Plants 11, no. 4: 496. https://doi.org/10.3390/plants11040496

APA Style

Kesawat, M. S., Kherawat, B. S., Singh, A., Dey, P., Routray, S., Mohapatra, C., Saha, D., Ram, C., Siddique, K. H. M., Kumar, A., Gupta, R., Chung, S.-M., & Kumar, M. (2022). Genome-Wide Analysis and Characterization of the Proline-Rich Extensin-like Receptor Kinases (PERKs) Gene Family Reveals Their Role in Different Developmental Stages and Stress Conditions in Wheat (Triticum aestivum L.). Plants, 11(4), 496. https://doi.org/10.3390/plants11040496

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