Home Life Sciences Genome-wide identification of MAPK family genes and their response to abiotic stresses in tea plant (Camellia sinensis)
Article Open Access

Genome-wide identification of MAPK family genes and their response to abiotic stresses in tea plant (Camellia sinensis)

  • Xinhao Liu , Min Zhao , Caihua Gu , Haodong Jiang , Junyan Sun and Jie Li EMAIL logo
Published/Copyright: September 8, 2022

Abstract

Mitogen-activated protein kinase (MAPK) cascades are conserved and universal signal transduction modules that play important roles in regulating stress responses in plants. Although MAP3K, MP2K, and MPK family in tea plant (Camellia sinensis) have been investigated, little is known about MPK family genes responding to various abiotic stresses in tea plant. In this study, we performed a comprehensive genome-wide analysis of the tea plant MAPKs (CsMPKs) family gene based on the genomic data of tea plants by bioinformatics-based methods. Here, 21 putative CsMPK genes were identified in the tea plant and divided into 4 subfamilies according to the homologous to Arabidopsis and their phylogenetic relationships. The gene structure and conserved motifs of these CsMPKs in the same group showed high similarity, suggesting that they were highly conserved and might have a similar function. The expression profiles of the CsMPK genes were further investigated by quantitative real-time reverse transcription PCR, indicating that many CsMPK genes were involved in response to cold, drought, heat, or heat combined with drought treatment, suggesting their potential roles in abiotic stress responses in tea plant. These results would provide valuable information for further exploring the functional characterization of CsMPK genes in tea plants.

1 Introduction

As sessile organisms, plants often encounter various adverse conditions, such as high and low temperatures, drought, salt stresses, and pathogen invasion, during their growth and development. To coordinate the biotic and abiotic stresses, plants have evolved complex signaling networks to perceive and transmit environmental stimuli. Mitogen-activated protein kinase (MAPK) cascades play vital roles in controlling intracellular response to extracellular signals [1].

MAPK cascades classically are composed of three subsequently protein kinases, namely MAP kinase kinase kinase (MAPKKK/MAP3K/MEKK), MAP kinase kinase (MAPKK/MKK/MEK), and MAP kinase (MAPK/MPK), but sometimes contain an upstream MAP4K [2]. In general, MAPKs are phosphorylated by MAPK kinase at their conserved threonine and tyrosine residues in the activation loop (T-loop), and in turn, MAPKK is phosphorylated in the S/T-X3–5-S/T motif of their activation loop by MAPKKKs when development or environmental signals incurred [3]. Compared with MAPKKKs and MAPKKs, MAPKs act at the bottom of MAPK cascades showing more complexity and sequence diversity. Plant MAPKs can be divided into four groups, including A, B, C, and D subfamilies based on the phosphorylation motif and phylogenetic relationships of amino acid sequence. Members of the A, B, and C subfamilies have a TEY motif in their active sites, while members in the D subfamily possess a TDY motif in their phosphorylation site [1].

In plants, numerous researches demonstrated that MPKs play indispensable regulatory roles in response to abiotic stresses [4]. Such as MPK3 and MPK6 were reported to play a critical role in response to drought, cold, salt, and heat stress in Arabidopsis [4,5], Jammes et al. demonstrated that AtMPK9 and AtMPK12 were preferentially and highly expressed in guard cells and positively regulated ROS-mediated ABA signaling [6]. A role for the MAPK module consisting of MEKK1-MKK1/2-MPK4 has been confirmed in response to drought, cold, and salt stresses in Arabidopsis [7,8]. In rice, OsMPK5 was shown to be a double-faced player in the stress response network: a positive regulator of abiotic stress tolerance and ABA-mediated defense responses to brown spots but a negative regulator of ET-controlled resistance to Magnaporthe oryzae [9]. OsMPK33 was demonstrated to play a negative role in salt tolerance through unfavorable ion homeostasis [10]. The MKK1–MPK4 module was reported to mediate salt signaling in rice [11]. OsMKK6 and OsMPK3 constitute a moderately low-temperature signaling pathway and regulate the tolerance to cold stress in rice [12]. In addition, in maize, MPK4, MPK5, MPK6, and MPK7 were active by diverse stresses, such as drought, low temperature, and salt [13].

In light of the importance of MAPK genes, the identification and characterization of MAPK members have been conducted in various plant species. Such as Arabidopsis thaliana genome contains 20 MAPKs [3], there are 17 MAPKs in rice (Oryza sativa) genome [14], 54 MAPK genes were identified in wheat (Triticum aestivum) [15], 28 MAPKs in cotton (Gossypium raimondii) [16], and 28 MAPKs in sunflower (Helianthus annuus) [17]. Tea plant (Camellia sinensis), an evergreen woody plant, is an important economic crop that is widely distributed in subtropical to tropical climate regions [18]. Due to the local climate changes, tea plant frequently experiences various environmental stresses during its lifecycle, and heat, drought, and low-temperature stresses are the main factors that affect the yield and quality of tea products [19,20]. Recently, high-quality genome sequencing data for tea plants were presented, which provides convenience for further understanding the tea plant CsMPK gene family [21,22,23]. In the present study, 21 CsMPK genes were identified using the tea plant genome by exploring the genomic data. The analysis of identified CsMPK genes in the sequence features, phylogenetic relationships, cis-elements in promoters, and dynamic expression patterns in response to various abiotic stresses was conducted. The results of this study provided valuable information for further investigation of tea plant MPK gene family.

2 Materials and methods

2.1 Identification and characterization of CsMPK gene family

The identification of the CsMPK gene family of tea plant was conducted according to the method described by Wang et al. with some modifications [24]. The MPK protein sequences of A. thaliana, O. sativa, Vitis vinifera, and Populus trichocarpa were downloaded from TAIR (http://www.Arabidopsis.org), the Rice Genome Annotation Project (http://rice.plantbiology.msu.edu/index.shtml), the V. vinifera proteome 12 × database (http://www.genoscope.cns.fr/externe/GenomeBrowser/Vitis/), and Ensembl database (http://plants.ensembl.org), respectively. Whole genome sequences of tea plant were downloaded from the Tea Plant Genome Database (http://tpia.teaplant.org/index.html). The MPK sequences of A. thaliana, O. sativa, V. vinifera, and P. trichocarpa were used as queries to search against tea plant proteins using BLASTP program with an e-value of 1 × 10−5 as the threshold. All the acquired MPK protein sequences from the four plants mentioned above were used to build the local Hidden Markov Model-based searches (HMMER) that were used to identify the CsMPKs in tea plant with the threshold of E < 1 × 10−20. Then, the BLASTP results were integrated with the HMMER hits and parsed by manual editing to remove redundant. Furthermore, these identified CsMPK members were confirmed by CDD (https://www.ncbi.nlm.nih.gov/cdd/), PFAM (http://pfam.xfam.org/), InterProScan database (http://www.ebi.ac.uk/interpro/), and SMARAT (http://smart.embl-heidelberg.de/). The molecular weight (MW) and isoelectric point (pI) of the CsMPKs were predicted using ProtParam (http://web.expasy.org/protparam/).

2.2 Sequence alignments and phylogenetic tree construction of CsMPK genes

Multiple alignments of amino acid sequences of MPKs were performed using ClustalW program with the default parameters [25]. Phylogenetic tree of MPKs from C. sinensis, A. thaliana, O. sativa, V. vinifera, and P. trichocarpa was generated by MEGA 6.0 software using a neighbor-joining (NJ) method with 1,000 bootstrap replicates [26], which was visualized in Figtree (http://tree.bio.ed.ac.uk/software/figtree/).

2.3 Gene structure and conserved motif analysis

The exon/intron structure of the CsMPK genes was retrieved from a gene annotation file (http://www.plantkingdomgdb.com/tea_tree/data/gff3/), and the diagrams were displayed by comparing CDSs and their corresponding gene sequences from genome using the online program GSDS 2.0 (http://gsds.cbi.pku.edu.cn/) [27]. The conserved motifs of CsMAPKs were analyzed by the program of MEME with the following parameters: any number of repetitions, maximum of 10 motifs, and an optimum motif width of 6–50 amino acid residues [28].

2.4 cis-element analysis of putative promoter regions

The upstream sequences (1.5 kb) of CsMPK genes were submitted to the PlantCARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) to predict various cis-regulatory elements according to consensus sequences, positional matrices, and individual sites on particular promoter sequences [29].

2.5 Plant materials and treatments

Two-year-old cutting seedlings of the tea cultivar “Longjing43” were grown in a controlled environment chamber maintained for 12 h day (25 ± 3°C) and 12 h night (20 ± 3°C), and 75% relative humidity for 2 weeks. Then, cold, heat, drought, and heat combined with drought stress treatments were performed. Untreated tea plants were used as control. For cold treatment, the tea plants were kept in a growth chamber at 4°C. The first and second tender leaves were harvested 24 and 48 h after the treatment from control and cold-treated plants. The treatments of heat, drought, and heat combined with drought stress were conducted simultaneously. For heat treatment, the tea plants were transferred to chambers maintained at 38°C. For drought treatment, the tea plants were irrigated with 20% (w/v) polyethylene glycol (PEG) 6000. For heat combined with drought stress, the tea plants were irrigated PEG 6000 and then quickly put in chambers at 38°C. The first and second tender leaves treated with drought, heat, heat combined with drought stresses, and control leaves were collected 24 h post-treatment, respectively. The collected leaves were immediately frozen in liquid nitrogen and then stored at −80°C for gene expression analysis. For each treatment, three independent biological replicates were performed.

2.6 Total RNA isolation and quantitative real-time reverse transcription PCR (qRT-PCR) expression analysis

Total RNA was extracted from young leaves of tea plants using the “TaKaRa MiniBEST Plant RNA Extraction” Kit (TaKaRa, Dalian, China) in accordance with the manufacturer’s protocol. RNase-free DNaseI was added to each sample to digest the genomic DNA. Then, cDNA was synthesized following the kit instruction (Promega Corp., Madison, WI, USA) and stored at −20°C. The primer sequences used in this experiment were designed based on gene sequences and the Beacon designer software (NJ, USA). CsGAPDH was used as the reference gene [30]. The primer sequences are shown in Table S1. The qRT-PCR was performed on a CFX96 Real-Time PCR Detection System (Bio-Rad Laboratories, Inc.) with three independent biological replicates. The reaction conditions were: denaturation at 95°C for 30 s, 40 cycles of 95°C for 5 s, and 60°C extensions for 30 s. Relative gene expression levels were calculated using the 2−∆∆CT method [31].

3 Results

3.1 Identification and characterization of CsMPK family genes in tea plant

To identify MPK encoding genes from tea plant genome, both BLASTP and Hidden Markov Model (HMM) searches were conducted using A. thaliana, O. sativa, V. vinifera, and P. trichocarpa MPK proteins as query sequences and then used these search results to identify all the CsMPK family members by CDD, PFAM, and SMART programs. Finally, a total of 21 CsMPK members were identified from tea plant genome. To distinguish from these identified CsMPKs, we provisionally named them according to the homologous gene in A. thaliana (Table 1). The sequence data of all the above CsMPKs are shown in Table S2. The information analysis of CsMPKs showed that the open reading frame (ORF) length for CsMPK genes ranged from 850 to 6,765 bp, and they encoded proteins ranging from 284 to 1,482 amino acids in size. The deduced protein MWs ranged from 32.65 to 169.32 kDa, and the protein isoelectric points (pIs) were between 4.9 and 9.37 (Table 1).

Table 1

The information of MAPK gene family in tea plant

Gene Homologous gene Sequence ID Group ORF length (bp) Protein length (aa) MW (kDa) pI
CsMPK1-1 AtMPK1 TEA031435 C 1,262 421 47.94 5.97
CsMPK1-2 AtMPK1 TEA016315 C 1,152 384 44.23 6.70
CsMPK3-1 AtMPK3 TEA026040 A 1,113 372 42.75 5.25
CsMPK3-2 AtMPK3 TEA020852 A 1,092 365 41.68 5.89
CsMPK3-3 AtMPK3 TEA020851 A 850 284 32.65 5.62
CsMPK3-4 AtMPK3 TEA017807 A 2,461 822 91.59 8.64
CsMPK3-5 AtMPK3 TEA017811 A 1,110 371 41.43 6.27
CsMPK4-1 AtMPK4 TEA006273 B 6,765 1,482 169.32 4.90
CsMPK4-2 AtMPK4 TEA021759 B 988 330 38.21 5.41
CsMPK4-3 AtMPK4 TEA006436 B 1,098 367 42.24 6.28
CsMPK4-4 AtMPK4 TEA007103 B 2,379 798 90.21 6.18
CsMPK6 AtMPK6 TEA024415 A 1,516 507 57.24 6.06
CsMPK7 AtMPK7 TEA032724 C 1,105 368 42.22 7.60
CsMPK9-1 AtMPK9 TEA022268 D 1,554 517 59.34 6.75
CsMPK9-2 AtMPK9 TEA015676 D 1,698 569 64.24 6.43
CsMPK15 AtMPK15 TEA012905 D 1,753 587 66.87 6.99
CsMPK16-1 AtMPK16 TEA026883 D 2,916 975 110.3 9.22
CsMPK16-2 AtMPK16 TEA031053 D 1,752 587 66.62 9.09
CsMPK19-1 AtMPK19 TEA018880 D 1,791 599 67.28 9.22
CsMPK19-2 AtMPK19 TEA022253 D 2,363 792 90.12 9.25
CsMPK20 AtMPK20 TEA004137 D 855 286 33.27 9.37

3.2 Phylogenetic relationship analysis of CsMPK genes

To analyze the phylogenetic relationships of the MPKs, an NJ phylogenetic tree was constructed with the amino acid sequences of 20 AtMPKs from Arabidopsis, 17 OsMPKs from O. sativa, 12 VvMPKs from grapevine, 23 PtMPKs from poplar, and 21 putative CsMPKs from tea plant (Figure 1). Through phylogenetic analyses, these MPKs were classified into four groups, designated as A, B, C, and D. Group D was the first largest clade, including eight CsMPKs, followed by group A (six members). Group C (three members) had the fewest number of CsMPKs; the rest belonged to group B, containing four CsMPKs. CsMPKs belonging to groups A, B, and C possess a TEY motif, and the D group has a TDY motif at the activation site (Figure S1).

Figure 1 
                  Phylogenetic relationship of putative MPK genes from A. thaliana (AtMPK), O. sativa (OsMPK), V. vinifera (VvMPK), P. trichocarpa (PtMPK), and C. sinensis (CsMPK). The phylogenetic tree with 100 bootstrap replicates was constructed by MEGA6.0 program with the NJ method. Letters A–D named as different groups of MPKs.
Figure 1

Phylogenetic relationship of putative MPK genes from A. thaliana (AtMPK), O. sativa (OsMPK), V. vinifera (VvMPK), P. trichocarpa (PtMPK), and C. sinensis (CsMPK). The phylogenetic tree with 100 bootstrap replicates was constructed by MEGA6.0 program with the NJ method. Letters A–D named as different groups of MPKs.

3.3 Conserved motifs analysis of CsMPKs

The conserved protein domains in CsMPKs were identified and compared by the online MEME program. A total of 11 conserved motifs were detected (Figure 2). The lengths of the identified motifs varied from 11 to 50 amino acids (motif logos shown in Figure S2). Two CsMPKs (CsMPK4-1 and CsMPK19-2) contained the maximum number of motifs (16), while the minimum number of motifs (6) occurred in CsMPK20. Furthermore, all of CsMPKs contained motifs 2, 3, and 6 (contained the TXY signature motif), but motif 8 was found only in CsMPK4-1. The alignment analysis and sequence logo revealed that the sequences of CsMPKs in the motif regions were highly conserved (Figures S1 and S2). In addition, the members in the same subfamily were found to share similar conserved motifs. For example, the majority of CsMPKs in groups A, B, and C possessed 9 motifs, whereas most members in group D had at least 10 motifs. Interestingly, motif 6 in group D members occurred two times. These CsMPKs in the same subgroup showed similar motif distribution indicating that these genes might have relatively high conservation.

Figure 2 
                  The distribution of conserved motifs of CsMPKs from tea plant according to the phylogenetic relationship. The motifs were identified by MEME program with the complete amino acid sequences of the 21 CsMPKs. Each motif is represented by one color box.
Figure 2

The distribution of conserved motifs of CsMPKs from tea plant according to the phylogenetic relationship. The motifs were identified by MEME program with the complete amino acid sequences of the 21 CsMPKs. Each motif is represented by one color box.

3.4 Gene structure analysis of CsMPK genes

Gene structures of the identified 21 CsMPK genes were analyzed by aligning tea plant genomic DNA sequences. As shown in Figure 3, the exon/intron structures of CsMPK genes could be divided into four groups on the basis of their phylogenetic relationship, and those CsMPK genes belonged to the same clade and shared similar intron/exon organization structures. In group A, the number of introns in CsMPKs varied from 4 to 7. Intron numbers in group B varied widely, from 4 (CsMPK4-2) to 17 (CsMPK4-4). Group C had no more than 3 introns, while group D had at least 10 introns except CsMPK20 and CsMPK19-1, containing 5 and 8 introns, respectively.

Figure 3 
                  Phylogenetic relationship and gene structure of CsMPK genes. The phylogenetic tree was created using MEGA6.0 program with the NJ method. Exons and introns are indicated by blue bars and thin lines, respectively. Gene models were drawn to scale as indicated on bottom.
Figure 3

Phylogenetic relationship and gene structure of CsMPK genes. The phylogenetic tree was created using MEGA6.0 program with the NJ method. Exons and introns are indicated by blue bars and thin lines, respectively. Gene models were drawn to scale as indicated on bottom.

3.5 Stress-related cis-elements analysis of CsMPK genes

To further explore the regulatory mechanisms and potential functions of CsMPK genes under stress conditions, the 1.5 kb promoter sequences of the CsMPK genes were used to analyze the cis-regulatory elements. Two categories of cis-elements were found in the promoter region of CsMPK genes (Figure 4 and Table S3). The first class was constituted by hormone-responsive elements, such as ABRE, ERE, and P-box. TGACG-motif and CGTCA-motif are involved in methyl jasmonate responsiveness, which existed in 10 CsMPKs. ERE is related to an ethylene-responsive element found in 12 CsMPKs. The ABRE element existed in nine CsMPKs. The results suggested that some phytohormones (e.g., ABA and JA) are involved in regulating the expression of some CsMPK genes. The other category was constituted by abiotic stress response-related cis-elements, including LTR, MBS, MYB, MYC, and ARE. MYC responds to abiotic stress signals, which existed in all of CsMPK genes. MYB responds to dehydration and ABA signals and were also found in all CsMPK genes except CsMPK3-3. Eight CsMPKs possessed LTRE, which responds to cold, drought, and ABA signals. MBS is involved in drought response and detected in 8 CsMPKs.

Figure 4 
                  Predicted cis-elements in the 1.5 kb upstream promoter regions of CsMPK genes. Different cis-elements are represented by different colors.
Figure 4

Predicted cis-elements in the 1.5 kb upstream promoter regions of CsMPK genes. Different cis-elements are represented by different colors.

3.6 Expression patterns of CsMPK genes under abiotic stresses

To examine the potential role of CsMPK genes in various abiotic stresses, the expression patterns of all of these genes were investigated under heat, drought, or cold stress conditions (Figure 5a and b). For cold treatment (Figure 5a), the expressions of six CsMPK genes (CsMPK3-1, CsMPK3-2, CsMPK3-3, CsMPK3-4, CsMPK4-3, CsMPK15, and CsMPK19-1) were induced to present a significant upregulation, while CsMPK9-1, CsMPK9-2, and CsMPK16-2 were downregulated under cold stress condition. CsMPK4-1 was upregulated at 24 h of treatment and then downregulated at 48 h. Among the upregulated CsMPK genes, four genes belonged to group A, one in group B, and another two in group D. In response to drought treatment (Figure 5b), most group A and group C members were downregulated, whereas CsMPK genes in groups B and D showed a significantly upregulated expression. For heat treatment, only six CsMPK genes (CsMPK3-3, CsMPK3-4, CsMPK3-2, CsMPK3-1, CsMPK1-2, and CsMPK16-2) displayed downregulation, and no gene was found to show significant up-regulation (Figure 5). Under heat combined with drought treatments (Figure 5b), six CsMPK genes (CsMPK3-2, CsMPK4-2, CsMPK4-3, CsMPK19-1, CsMPK19-2, and CsMPK15) were significantly upregulated, and another six genes (CsMPK3-3, CsMPK3-4, CsMPK3-1, CsMPK1-1, CsMPK1-2, and CsMPK16-2) displayed a downregulated expression trend (Figure 5b). In addition, no changes in the expression of three CsMPK genes (CsMPK3-5, CsMPK4-4, and CsMPK20) were observed under cold, heat, drought, and heat combined with drought treatments. Notably, CsMPK3-2 showed differential expression under three stress conditions (cold, heat, and heat combined drought), indicating that it might play an important role in response to several stresses.

Figure 5 
                  Hierarchical clustering of the expression patterns of CsMPK genes in tea plant leaves under cold, drought, heat, and heat combined with drought stresses. The red and green colors indicate higher or lower expression levels, respectively. (a) Heatmap of the CsMPK genes expression under cold treatment. (b) Heatmap of the CsMPK genes expression under drought, heat, and heat combined with drought treatments. DT, HT, and HD were denoted as drought, heat, and heat combined with drought treatments, respectively.
Figure 5

Hierarchical clustering of the expression patterns of CsMPK genes in tea plant leaves under cold, drought, heat, and heat combined with drought stresses. The red and green colors indicate higher or lower expression levels, respectively. (a) Heatmap of the CsMPK genes expression under cold treatment. (b) Heatmap of the CsMPK genes expression under drought, heat, and heat combined with drought treatments. DT, HT, and HD were denoted as drought, heat, and heat combined with drought treatments, respectively.

4 Discussion

Due to the importance of plant MPKs in multifaceted biological processes such as growth, development, and stress responses [32,33], lots of MPK family members were identified in many plants through whole genome sequencing [34,35]. Recently, the tea plant genome database has been published [21,22,23], which provides the basis for the comprehensive analysis of MAPK cascade members. As the upstream components of the MAPK cascade, MAPKKK family members have been described in detail [36]. Chatterjee et al. analyzed the expression profiles and regulatory network of the MAPK cascade (MAPKK–MAPK) gene family of the tea plant [37]. Tea is the main industrial crop of the Xinyang district, which is often challenged by cold, heat, drought, and heat combined with drought stresses. To implore the tolerance mechanism of tea plant to above-mentioned stresses, we performed a detailed analysis of the CsMPK family genes and the expression characteristics of these genes under abiotic stresses. In this study, 21 CsMPK genes were identified in the tea plant genome. The number of CsMPKs is close to the MPK gene number in Arabidopsis (20 MPKs) and poplar (23 MPKs) (Nicole et al. identified 21 PtMPK genes from P. trichocarpa genome [38]. The variation in PtMPKs amount within the same species in the two studies might come as a result of different parameters employed during HMM profiling and motif analysis), but the genome size of tea plant (∼2.85 Gb) is approximately 23 times that of the Arabidopsis genome (∼125 Mb) and 6 times than poplar genome (∼485 Mb) [39]. The number of CsMPK members in tea plant was more than in grapevine (12 members), coffee (12 members), and kiwifruit (18 members) and slightly lower than that in apple (26 members). The genome sizes of these plants (grapevine is ∼490 Mb, coffee is ∼710 Mb, apple is ∼742 Mb, and kiwifruit is ∼616 Mb) are significantly smaller than the tea plant genome [40,41,42,43]. The phenomenon indicated that the number of MPK genes is irrelevant to the size of the plant genome.

The sequence alignment and phylogenetic tree analysis revealed four distinct groups, which were consistent with the MPKs previously identified in Arabidopsis [3], rice [14], poplar [38], and sunflower [17]. In group A, tea plants had four extra copies of MPK3 that might have been because of gene duplications. Such extra gene copies have also been observed in sunflowers [17] and apples [44]. Four copies of CsMPK4 were also found in group B. In plants, gene duplications play vital roles in the genomic expansion and diversification of gene function [45]. Two rounds of whole genome duplication events have occurred in the tea tree genome, which resulted in the expansion of several gene classes [22,23]. More copies of CsMPKs might suggest that their capabilities were broken down into more detailed functions.

The exon/intron structure plays an important role in organismal evolution [46]. In the study, we observed that the positions, sizes, and sequences of the introns differed between CsMPK genes (Figure 3). Members in group C had the fewest introns, while group D possessed the most introns. Similar exon/intron structures were found in apples [44]. Generally, genes with fewer introns could be rapidly activated in response to environmental stresses [47]. MPK genes are one of the rapidly induced genes under different environmental challenges [3]. The less of introns fit in with the needs for the rapid activation of MPK genes. In this study, we found that group A members with fewer introns were induced by various stress treatments, whereas genes with many introns, such as CsMPK4-4, CsMPK9-1, and CsMPK16-1, did not show any significant expression changes under cold, heat, and drought stresses (Figure 5). In addition, CsMPK members with similar conserved motifs and exon/intron structure showed a tendency to cluster into one subgroup (Figures 2 and 3), which implied functional consistency among these CsMPKs in the same subfamily. This correlation between motif arrangement and intron numbers supported the previous classification of the CsMPK genes. These results might provide an excellent reference to explore the functions of the CsMPKs.

Numerous reports demonstrated that plant MPK genes play important roles in determining response to abiotic stresses [48]. For example, AtMAPK3, AtMAPK4, ZmMAPK1, and SiMAPK7 were reported to regulate abiotic stress responses [49,50,51]. In this study, the expression levels of the CsMPK3-1, CsMPK3-3, and CsMPK3-4, which were the orthologue of the AtMAPK3 gene, were induced by cold (Figure 5). Furthermore, these genes possessed at least four types of cis-elements in the promoter region (Figure 4), suggesting that they were related to the transcriptional control of stress responses. Notably, CsMPK3-2 was significantly upregulated under cold and heat combined with drought, implying that CsMPK3-2 might be an important regulator in response to abiotic stresses. The expressions of CsMPK3-1, CsMPK3-2, CsMPK3-3, and CsMPK3-4 were downregulated under heat stress. These results are similar to the previous report in which the SlMPK3 was proposed to negatively regulate heat stress in tomato plants [52]. However, CsMPK3-5 had nine types of cis-elements and fewer introns, and no change was found in the expression level. The results suggested that the cis-elements of gene were not all of decisive factor in responding to stresses, or they could be induced by other stresses. CsMPK4-3, a member of group B, was also induced by drought stress, consisting of the AtMPK4 [53]. The CsMPK genes belonging to group D have not been as well studied as those of groups A and B. In the present study, the expression of CsMPK19-1 and CsMPK19-2 was increased in cold and heat combined with drought stress, and CsMPK15 and CsMPK9-2 were upregulated under drought stress, suggesting the possible roles of these CsMPK genes in abiotic stress responses. In addition, we found that most CsMPK genes possessed at least one hormone-responsive cis-element, including ABA, JA, and ethylene-responsive elements (Table S3). The relationship between MAPK signaling pathways and ABA, JA, SA, and ethylene in plant abiotic stress responses was characterized [54,55]. However, little is known about how the CsMPK genes respond to phytohormones in tea plant. These questions need further research. The regulatory mechanism of CsMPK genes under abiotic stress is complex. More research is needed to explore the specific functions of the CsMPK family genes in tea plants through additional experiments.

5 Conclusions

In this study, a genome-wide analysis of CsMPK family genes in tea plant was performed, and 21 putative CsMPK genes were identified. The identified CsMPK genes were classified into four subfamilies (groups A, B, C, and D). The gene structure and conserved motifs of the CsMPK genes supported the classification results. Numerous cis-elements associated with abiotic stresses were found in the promoter regions of CsMPK genes. Further expression profile analysis using qRT-PCR showed that most CsMPK genes have a positive or negative response to cold, drought, or heat stresses, indicating that CsMPK genes play important roles in abiotic stress response. This work provided useful information on the CsMPK family genes of tea plant and will contribute to investigating the function and regulatory mechanism of CsMPK genes.

  1. Funding information: This work was supported by the Youth Foundation of Xinyang Agriculture and Forestry University (2019LG007), the Foundation of Central Laboratory of Xinyang Agriculture and Forestry University (FCL202009 and FCL202015), Natural Science Foundation of Henan Province of China (202300410025), and Key Scientific Research Projects of Colleges and Universities of Henan Province of China (20B210019).

  2. Author contributions: Conceptualization, L.X.H. and L.J.; methodology, investigation, and formal analysis, L.X.H., Z.M., G.C.H., and J.H.D.; writing original draft, L.J.; writing review and editing, S.J.Y. and L.J.; visualization, supervision, and project administration, L.X.H., S.J.Y., and L.J.; funding acquisition, L.X.H. and L.J. All authors have read and agreed to the published version of the manuscript.

  3. Conflict of interest: Authors state no conflict of interest.

  4. Data availability statement: The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

References

[1] Cristina M, Petersen M, Mundy J. Mitogen-activated protein kinase signaling in plants. Annu Rev Plant Biol. 2010;61:621–49.10.1146/annurev-arplant-042809-112252Search in Google Scholar

[2] Champion A, Picaud A, Henry Y. Reassessing the MAP3K and MAP4K relationships. Trends Plant Sci. 2004;9:123–9.10.1016/j.tplants.2004.01.005Search in Google Scholar

[3] MAPK Group. Mitogen-activated protein kinase cascades in plants: A new nomenclature. Trends Plant Sci. 2002;7:301–8.10.1016/S1360-1385(02)02302-6Search in Google Scholar

[4] Raja V, Majeed U, Kang H, Andrabi KI, John R. Abiotic stress: Interplay between ROS, hormones and MAPKs. Env Exp Bot. 2017;137:142–57.10.1016/j.envexpbot.2017.02.010Search in Google Scholar

[5] de Zelicourt A, Colcombet J, Hirt H. The role of MAPK modules and ABA during abiotic stress signaling. Trends Plant Sci. 2016;21:677–85.10.1016/j.tplants.2016.04.004Search in Google Scholar PubMed

[6] Jammes F, Song C, Shin D, Munemasa S, Kwak JM. MAP kinases MPK9 and MPK12 are preferentially expressed in guard cells and positively regulate ROS-mediated ABA signaling. PNAS. 2009;106:20520–5.10.1073/pnas.0907205106Search in Google Scholar PubMed PubMed Central

[7] Furuya T, Matsuoka D, Nanmori T. Membrane rigidification functions upstream of the MEKK1-MKK2-MPK4 cascade during cold acclimation in Arabidopsis thaliana. FEBS Lett. 2014;588:2025–30.10.1016/j.febslet.2014.04.032Search in Google Scholar PubMed

[8] Qiu JL, Zhou L, Yun BW, Nielsen HB, Fiil BK, Petersen K, et al. Arabidopsis mitogen-activated protein kinase kinases MKK1 and MKK2 have overlapping functions in defense signaling mediated by MEKK1, MPK4, and MKS11. Plant Physiol. 2008;148:212–22.10.1104/pp.108.120006Search in Google Scholar PubMed PubMed Central

[9] Sharma R, Vleesschauwer D, Sharma MK, Ronald PC. Recent advances in dissecting stress-regulatory crosstalk in rice. Mol plant. 2013;6:250–60.10.1093/mp/sss147Search in Google Scholar PubMed

[10] Lee SK, Kim BG, Kwon TR, Jeong MJ, Park SR, Lee JW, et al. Overexpression of the mitogen-activated protein kinase gene OsMAPK33 enhances sensitivity to salt stress in rice (Oryza sativa L.). J Biosci. 2011;36:139–51.10.1007/s12038-011-9002-8Search in Google Scholar PubMed

[11] Wang F, Wen J, Zhang W. The mitogen-activated protein kinase cascade MKK1-MPK4 mediates salt signaling in rice. Plant Sci. 2014;227:181–9.10.1016/j.plantsci.2014.08.007Search in Google Scholar PubMed

[12] Xie G, Kato H, Imai R. Biochemical identification of the OsMKK6-OsMPK3 signalling pathway for chilling stress tolerance in rice. Biochem J. 2012;443:95–102.10.1042/BJ20111792Search in Google Scholar PubMed

[13] Wu T, Kong X, Zong X, Li D, Li D. Expression analysis of five maize MAP kinase genes in response to various abiotic stresses and signal molecules. Mol Biol Rep. 2011;38:3967–75.10.1007/s11033-010-0514-3Search in Google Scholar PubMed

[14] Reyna NS, Yang Y. Molecular analysis of the rice MAP kinase gene family in relation to Magnaporthe grisea infection. MPMI. 2006;19:530–40.10.1094/MPMI-19-0530Search in Google Scholar PubMed

[15] Zhan H, Yue H, Zhao X, Wang M, Song W, Nie X. Genome-wide identification and analysis of MAPK and MAPKK gene families in bread wheat (Triticum aestivum L.). Genes. 2017;8:284.10.3390/genes8100284Search in Google Scholar PubMed PubMed Central

[16] Zhang X, Xu X, Yu Y, Chen C, Wang J, Cai C, et al. Integration analysis of MKK and MAPK family members highlights potential MAPK signaling modules in cotton. Sci Rep. 2016;6:29781.10.1038/srep29781Search in Google Scholar PubMed PubMed Central

[17] Neupane S, Schweitzer SE, Neupane A. Identification and characterization of mitogen-activated protein kinase (MAPK) genes in sunflower (Helianthus annuus L.). Plants (Basel). 2019;8:28.10.3390/plants8020028Search in Google Scholar PubMed PubMed Central

[18] Liang C, Zhou Z, Yang Y. Genetic improvement and breeding of tea plant (Camellia sinensis) in China: from individual selection to hybridization and molecular breeding. Euphytica. 2007;154:239–48.10.1007/s10681-006-9292-3Search in Google Scholar

[19] Zheng C, Wang Y, Ding Z, Zhao L. Global transcriptional analysis reveals the complex relationship between tea quality, leaf senescence and the responses to cold-drought combined stress in Camellia sinensis. Front Plant Sci. 2016;7:1858.10.3389/fpls.2016.01858Search in Google Scholar PubMed PubMed Central

[20] Zhu X, Liao J, Xia X, Xiong F, Li Y, Shen J, et al. Physiological and iTRAQ-based proteomic analyses reveal the function of exogenous γ-aminobutyric acid (GABA) in improving tea plant (Camellia sinensis L.) tolerance at cold temperature. BMC plant biol. 2019;19:43.10.1186/s12870-019-1646-9Search in Google Scholar PubMed PubMed Central

[21] Wei C, Yang H, Wang S, Zhao J, Wan X. Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality. PNAS. 2018;115:4151–8.10.1073/pnas.1719622115Search in Google Scholar PubMed PubMed Central

[22] Xia E, Zhang H, Sheng J, Li K, Zhang Q, Kim C, et al. The tea tree genome provides insights into tea flavor and independent evolution of caffeine biosynthesis. Mol Plant. 2017;10:866–77.10.1016/j.molp.2017.04.002Search in Google Scholar PubMed

[23] Wang X, Feng H, Chang Y, Ma C, Wang L, Hao X, et al. Population sequencing enhances understanding of tea plant evolution. Nat Commun. 2020;11:4447.10.1038/s41467-020-18228-8Search in Google Scholar PubMed PubMed Central

[24] Wang G, Wang T, Jia ZH, Xuan JP, Pan DL, Guo ZR, et al. Genome-wide bioinformatics analysis of MAPK gene family in Kiwifruit (Actinidia Chinensis). Int J Mol Sci. 2018;19:2510–6.10.3390/ijms19092510Search in Google Scholar PubMed PubMed Central

[25] Hung JH, Weng Z. Sequence alignment and homology search with BLAST and ClustalW. Cold Spring Harb Protoc. 2016;11:1–6.10.1101/pdb.prot093088Search in Google Scholar PubMed

[26] Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis Version 6.0. Mol Biol Evol. 2013;30:2725–9.10.1093/molbev/mst197Search in Google Scholar PubMed PubMed Central

[27] Hu B, Jin J, Guo AY, Zhang H, Gao G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics. 2014;31:1296–7.10.1093/bioinformatics/btu817Search in Google Scholar PubMed PubMed Central

[28] Bailey TL, Mikael B, Buske FA, Martin F, Grant CE, Luca C, et al. MEME Suite: tools for motif discovery and searching. Nucleic Acids Res. 2009;37:202–8.10.1093/nar/gkp335Search in Google Scholar PubMed PubMed Central

[29] Lescot M. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002;30:325–7.10.1093/nar/30.1.325Search in Google Scholar PubMed PubMed Central

[30] Wu Z, Tian C, Jiang Q, Li X, Zhuang J. Selection of suitable reference genes for qRT-PCR normalization during leaf development and hormonal stimuli in tea plant (Camellia sinensis). Sci Rep. 2016;6:19748.10.1038/srep19748Search in Google Scholar PubMed PubMed Central

[31] Livak K, Schmittgen T. Analysis of relative gene expression data using real-time quantitative PCR and the 2−∆∆CT method. Methods. 2001;25:402–8.10.1006/meth.2001.1262Search in Google Scholar PubMed

[32] Andrasi N, Rigo G, Zsigmond L, Perz-Salamo I, Papdi C, Klement E, et al. The mitogen-activated protein kinase 4-phosphorylated heat shock factor A4A regulates responses to combined salt and heat stresses. J Exp Bot. 2019;70:4903–18.10.1093/jxb/erz217Search in Google Scholar PubMed PubMed Central

[33] Xu J, Zhang S. Mitogen-activated protein kinase cascades in signaling plant growth and development. Trends Plant Sci. 2015;20:56–64.10.1016/j.tplants.2014.10.001Search in Google Scholar PubMed

[34] McNeece BT, Sharma K, Lawrence GW, Lawrence KS, Klink VP. The mitogen activated protein kinase (MAPK) gene family functions as a cohort during the Glycine max defense response to Heterodera glycines. Plant Physiol Biochem. 2019;137:25–41.10.1016/j.plaphy.2019.01.018Search in Google Scholar PubMed

[35] Yan Y, Wang LZ, Ding Z, Tie W, Ding X, Zeng C, et al. Genome-wide identification and expression analysis of the mitogen-activated protein kinase gene family in Cassava. Front Plant Sci. 2016;7:1294.10.3389/fpls.2016.01294Search in Google Scholar PubMed PubMed Central

[36] Paul A, Srivastava AP, Subrahmanya S, Shen GX, Mishra N. In-silico genome wide analysis of Mitogen activated protein kinase kinase kinase gene family in C. sinensis. PLoS One. 2021;16(11):e0258657.10.1371/journal.pone.0258657Search in Google Scholar PubMed PubMed Central

[37] Chatterjee A, Paul A, Unnati GM, Rajput R, Biswas T, Kar T, et al. MAPK cascade gene family in Camellia sinensis: In-silico identification, expression profiles and regulatory network analysis. BMC Genomics. 2020;21:613.10.1186/s12864-020-07030-xSearch in Google Scholar PubMed PubMed Central

[38] Nicole MC, Hamel LP, Morency MJ, Beaudoin N, Ellis BE, Seguin A. MAP-ping genomic organization and organ-specific expression profiles of poplar MAP kinases and MAP kinase kinases. BMC Genomics. 2006;7:223.10.1186/1471-2164-7-223Search in Google Scholar PubMed PubMed Central

[39] Tuskan GA, DiFazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, et al. The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006;313:1596–604.10.1126/science.1128691Search in Google Scholar PubMed

[40] Daccord N, Celton JM, Linsmith G, Becker C, Choisne N, Schijlen E, et al. High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development. Nat Genet. 2017;49:1099–106.10.1038/ng.3886Search in Google Scholar PubMed

[41] Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007;449:463–7.10.1038/nature06148Search in Google Scholar PubMed

[42] Denoeud F, Carretero-Paulet L, Dereeper A, Droc G, Guyot R, Pietrella M, et al. The coffee genome provides insight into the convergent evolution of caffeine biosynthesis. Science. 2014;345:1181–4.10.1126/science.1255274Search in Google Scholar PubMed

[43] Huang SX, Ding J, Deng DJ, Tang W, Sun HH, Liu DY, et al. Draft genome of the kiwifruit Actinidia chinensis. Nat commun. 2013;4:2640.10.1038/ncomms3640Search in Google Scholar PubMed PubMed Central

[44] Zhang S, Xu R, Luo X, Jiang Z, Shu H. Genome-wide identification and expression analysis of MAPK and MAPKK gene family in Malus domestica. Gene. 2013;531:377–87.10.1016/j.gene.2013.07.107Search in Google Scholar PubMed

[45] Panchy N, Lehti-Shiu M, Shiu SH. Evolution of gene duplication in plants. Plant Physiol. 2016;171(4):2294–316.10.1104/pp.16.00523Search in Google Scholar PubMed PubMed Central

[46] Xu G, Guo C, Shan H, Kong H. Divergence of duplicate genes in exon-intron structure. PNAS. 2012;109:1187–92.10.1073/pnas.1109047109Search in Google Scholar PubMed PubMed Central

[47] Jeffares DC, Penkett CJ, Bahler J. Rapidly regulated genes are intron poor. Trends Genet. 2008;24:375–8.10.1016/j.tig.2008.05.006Search in Google Scholar PubMed

[48] Sinha AK, Jaggi M, Raghuram B, Tuteja N. Mitogen-activated protein kinase signaling in plants under abiotic stress. Plant Signa Behav. 2011;6:196–203.10.4161/psb.6.2.14701Search in Google Scholar PubMed PubMed Central

[49] Ortiz-Masia D, Perez-Amador MA, Carbonell J, Marcote MJ. Diverse stress signals activate the C1 subgroup MAP kinases of Arabidopsis. Febs Lett. 2007;581:1834–40.10.1016/j.febslet.2007.03.075Search in Google Scholar PubMed

[50] Yu L, Yan J, Yang Y, He L, Zhu W. Enhanced tolerance to chilling stress in tomato by overexpression of a mitogen-activated protein kinase, SlMPK7. Plant Mol Biol Rep. 2006;34:76–88.10.1007/s11105-015-0897-3Search in Google Scholar

[51] Wu L, Zu X, Zhang H, Wu L, Xi Z, Chen Y. Overexpression of ZmMAPK1 enhances drought and heat stress in transgenic Arabidopsis thaliana. Plant Mol Biol. 2015;88:429–43.10.1007/s11103-015-0333-ySearch in Google Scholar PubMed

[52] Yu W, Wang L, Zhao R, Sheng J, Zhang S, Li R, et al. Knockout of SlMAPK3 enhances tolerance to heat stress involving ROS homeostasis in tomato plants. BMC Plant Biol. 2019;19:354.10.1186/s12870-019-1939-zSearch in Google Scholar PubMed PubMed Central

[53] Ichimura K, Mizoguchi T, Yoshida R, Yuasa T, Shinozaki K. Various abiotic stresses rapidly activate Arabidopsis MAP kinases ATMPK4 and ATMPK6. Plant J. 2000;24:655–65.10.1046/j.1365-313x.2000.00913.xSearch in Google Scholar PubMed

[54] Danquah A, Zelicourt A, Colcombet J, Hirt H. The role of ABA and MAPK signaling pathways in plant abiotic stress responses. Biotechnol Adv. 2014;32:40–52.10.1016/j.biotechadv.2013.09.006Search in Google Scholar PubMed

[55] Meldau S, Ullman-Zeunert L, Govind G, Bartram S, Baldwin IT. MAPK-dependent JA and SA signalling in Nicotiana attenuata affects plant growth and fitness during competition with conspecifics. BMC Plant Biol. 2012;12:213.10.1186/1471-2229-12-213Search in Google Scholar PubMed PubMed Central

Received: 2022-01-03
Revised: 2022-05-18
Accepted: 2022-06-20
Published Online: 2022-09-08

© 2022 Xinhao Liu et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Biomedical Sciences
  2. Effects of direct oral anticoagulants dabigatran and rivaroxaban on the blood coagulation function in rabbits
  3. The mother of all battles: Viruses vs humans. Can humans avoid extinction in 50–100 years?
  4. Knockdown of G1P3 inhibits cell proliferation and enhances the cytotoxicity of dexamethasone in acute lymphoblastic leukemia
  5. LINC00665 regulates hepatocellular carcinoma by modulating mRNA via the m6A enzyme
  6. Association study of CLDN14 variations in patients with kidney stones
  7. Concanavalin A-induced autoimmune hepatitis model in mice: Mechanisms and future outlook
  8. Regulation of miR-30b in cancer development, apoptosis, and drug resistance
  9. Informatic analysis of the pulmonary microecology in non-cystic fibrosis bronchiectasis at three different stages
  10. Swimming attenuates tumor growth in CT-26 tumor-bearing mice and suppresses angiogenesis by mediating the HIF-1α/VEGFA pathway
  11. Characterization of intestinal microbiota and serum metabolites in patients with mild hepatic encephalopathy
  12. Functional conservation and divergence in plant-specific GRF gene family revealed by sequences and expression analysis
  13. Application of the FLP/LoxP-FRT recombination system to switch the eGFP expression in a model prokaryote
  14. Biomedical evaluation of antioxidant properties of lamb meat enriched with iodine and selenium
  15. Intravenous infusion of the exosomes derived from human umbilical cord mesenchymal stem cells enhance neurological recovery after traumatic brain injury via suppressing the NF-κB pathway
  16. Effect of dietary pattern on pregnant women with gestational diabetes mellitus and its clinical significance
  17. Potential regulatory mechanism of TNF-α/TNFR1/ANXA1 in glioma cells and its role in glioma cell proliferation
  18. Effect of the genetic mutant G71R in uridine diphosphate-glucuronosyltransferase 1A1 on the conjugation of bilirubin
  19. Quercetin inhibits cytotoxicity of PC12 cells induced by amyloid-beta 25–35 via stimulating estrogen receptor α, activating ERK1/2, and inhibiting apoptosis
  20. Nutrition intervention in the management of novel coronavirus pneumonia patients
  21. circ-CFH promotes the development of HCC by regulating cell proliferation, apoptosis, migration, invasion, and glycolysis through the miR-377-3p/RNF38 axis
  22. Bmi-1 directly upregulates glucose transporter 1 in human gastric adenocarcinoma
  23. Lacunar infarction aggravates the cognitive deficit in the elderly with white matter lesion
  24. Hydroxysafflor yellow A improved retinopathy via Nrf2/HO-1 pathway in rats
  25. Comparison of axon extension: PTFE versus PLA formed by a 3D printer
  26. Elevated IL-35 level and iTr35 subset increase the bacterial burden and lung lesions in Mycobacterium tuberculosis-infected mice
  27. A case report of CAT gene and HNF1β gene variations in a patient with early-onset diabetes
  28. Study on the mechanism of inhibiting patulin production by fengycin
  29. SOX4 promotes high-glucose-induced inflammation and angiogenesis of retinal endothelial cells by activating NF-κB signaling pathway
  30. Relationship between blood clots and COVID-19 vaccines: A literature review
  31. Analysis of genetic characteristics of 436 children with dysplasia and detailed analysis of rare karyotype
  32. Bioinformatics network analyses of growth differentiation factor 11
  33. NR4A1 inhibits the epithelial–mesenchymal transition of hepatic stellate cells: Involvement of TGF-β–Smad2/3/4–ZEB signaling
  34. Expression of Zeb1 in the differentiation of mouse embryonic stem cell
  35. Study on the genetic damage caused by cadmium sulfide quantum dots in human lymphocytes
  36. Association between single-nucleotide polymorphisms of NKX2.5 and congenital heart disease in Chinese population: A meta-analysis
  37. Assessment of the anesthetic effect of modified pentothal sodium solution on Sprague-Dawley rats
  38. Genetic susceptibility to high myopia in Han Chinese population
  39. Potential biomarkers and molecular mechanisms in preeclampsia progression
  40. Silencing circular RNA-friend leukemia virus integration 1 restrained malignancy of CC cells and oxaliplatin resistance by disturbing dyskeratosis congenita 1
  41. Endostar plus pembrolizumab combined with a platinum-based dual chemotherapy regime for advanced pulmonary large-cell neuroendocrine carcinoma as a first-line treatment: A case report
  42. The significance of PAK4 in signaling and clinicopathology: A review
  43. Sorafenib inhibits ovarian cancer cell proliferation and mobility and induces radiosensitivity by targeting the tumor cell epithelial–mesenchymal transition
  44. Characterization of rabbit polyclonal antibody against camel recombinant nanobodies
  45. Active legumain promotes invasion and migration of neuroblastoma by regulating epithelial-mesenchymal transition
  46. Effect of cell receptors in the pathogenesis of osteoarthritis: Current insights
  47. MT-12 inhibits the proliferation of bladder cells in vitro and in vivo by enhancing autophagy through mitochondrial dysfunction
  48. Study of hsa_circRNA_000121 and hsa_circRNA_004183 in papillary thyroid microcarcinoma
  49. BuyangHuanwu Decoction attenuates cerebral vasospasm caused by subarachnoid hemorrhage in rats via PI3K/AKT/eNOS axis
  50. Effects of the interaction of Notch and TLR4 pathways on inflammation and heart function in septic heart
  51. Monosodium iodoacetate-induced subchondral bone microstructure and inflammatory changes in an animal model of osteoarthritis
  52. A rare presentation of type II Abernethy malformation and nephrotic syndrome: Case report and review
  53. Rapid death due to pulmonary epithelioid haemangioendothelioma in several weeks: A case report
  54. Hepatoprotective role of peroxisome proliferator-activated receptor-α in non-cancerous hepatic tissues following transcatheter arterial embolization
  55. Correlation between peripheral blood lymphocyte subpopulations and primary systemic lupus erythematosus
  56. A novel SLC8A1-ALK fusion in lung adenocarcinoma confers sensitivity to alectinib: A case report
  57. β-Hydroxybutyrate upregulates FGF21 expression through inhibition of histone deacetylases in hepatocytes
  58. Identification of metabolic genes for the prediction of prognosis and tumor microenvironment infiltration in early-stage non-small cell lung cancer
  59. BTBD10 inhibits glioma tumorigenesis by downregulating cyclin D1 and p-Akt
  60. Mucormycosis co-infection in COVID-19 patients: An update
  61. Metagenomic next-generation sequencing in diagnosing Pneumocystis jirovecii pneumonia: A case report
  62. Long non-coding RNA HOXB-AS1 is a prognostic marker and promotes hepatocellular carcinoma cells’ proliferation and invasion
  63. Preparation and evaluation of LA-PEG-SPION, a targeted MRI contrast agent for liver cancer
  64. Proteomic analysis of the liver regulating lipid metabolism in Chaohu ducks using two-dimensional electrophoresis
  65. Nasopharyngeal tuberculosis: A case report
  66. Characterization and evaluation of anti-Salmonella enteritidis activity of indigenous probiotic lactobacilli in mice
  67. Aberrant pulmonary immune response of obese mice to periodontal infection
  68. Bacteriospermia – A formidable player in male subfertility
  69. In silico and in vivo analysis of TIPE1 expression in diffuse large B cell lymphoma
  70. Effects of KCa channels on biological behavior of trophoblasts
  71. Interleukin-17A influences the vulnerability rather than the size of established atherosclerotic plaques in apolipoprotein E-deficient mice
  72. Multiple organ failure and death caused by Staphylococcus aureus hip infection: A case report
  73. Prognostic signature related to the immune environment of oral squamous cell carcinoma
  74. Primary and metastatic squamous cell carcinoma of the thyroid gland: Two case reports
  75. Neuroprotective effects of crocin and crocin-loaded niosomes against the paraquat-induced oxidative brain damage in rats
  76. Role of MMP-2 and CD147 in kidney fibrosis
  77. Geometric basis of action potential of skeletal muscle cells and neurons
  78. Babesia microti-induced fulminant sepsis in an immunocompromised host: A case report and the case-specific literature review
  79. Role of cerebellar cortex in associative learning and memory in guinea pigs
  80. Application of metagenomic next-generation sequencing technique for diagnosing a specific case of necrotizing meningoencephalitis caused by human herpesvirus 2
  81. Case report: Quadruple primary malignant neoplasms including esophageal, ureteral, and lung in an elderly male
  82. Long non-coding RNA NEAT1 promotes angiogenesis in hepatoma carcinoma via the miR-125a-5p/VEGF pathway
  83. Osteogenic differentiation of periodontal membrane stem cells in inflammatory environments
  84. Knockdown of SHMT2 enhances the sensitivity of gastric cancer cells to radiotherapy through the Wnt/β-catenin pathway
  85. Continuous renal replacement therapy combined with double filtration plasmapheresis in the treatment of severe lupus complicated by serious bacterial infections in children: A case report
  86. Simultaneous triple primary malignancies, including bladder cancer, lymphoma, and lung cancer, in an elderly male: A case report
  87. Preclinical immunogenicity assessment of a cell-based inactivated whole-virion H5N1 influenza vaccine
  88. One case of iodine-125 therapy – A new minimally invasive treatment of intrahepatic cholangiocarcinoma
  89. S1P promotes corneal trigeminal neuron differentiation and corneal nerve repair via upregulating nerve growth factor expression in a mouse model
  90. Early cancer detection by a targeted methylation assay of circulating tumor DNA in plasma
  91. Calcifying nanoparticles initiate the calcification process of mesenchymal stem cells in vitro through the activation of the TGF-β1/Smad signaling pathway and promote the decay of echinococcosis
  92. Evaluation of prognostic markers in patients infected with SARS-CoV-2
  93. N6-Methyladenosine-related alternative splicing events play a role in bladder cancer
  94. Characterization of the structural, oxidative, and immunological features of testis tissue from Zucker diabetic fatty rats
  95. Effects of glucose and osmotic pressure on the proliferation and cell cycle of human chorionic trophoblast cells
  96. Investigation of genotype diversity of 7,804 norovirus sequences in humans and animals of China
  97. Characteristics and karyotype analysis of a patient with turner syndrome complicated with multiple-site tumors: A case report
  98. Aggravated renal fibrosis is positively associated with the activation of HMGB1-TLR2/4 signaling in STZ-induced diabetic mice
  99. Distribution characteristics of SARS-CoV-2 IgM/IgG in false-positive results detected by chemiluminescent immunoassay
  100. SRPX2 attenuated oxygen–glucose deprivation and reperfusion-induced injury in cardiomyocytes via alleviating endoplasmic reticulum stress-induced apoptosis through targeting PI3K/Akt/mTOR axis
  101. Aquaporin-8 overexpression is involved in vascular structure and function changes in placentas of gestational diabetes mellitus patients
  102. Relationship between CRP gene polymorphisms and ischemic stroke risk: A systematic review and meta-analysis
  103. Effects of growth hormone on lipid metabolism and sexual development in pubertal obese male rats
  104. Cloning and identification of the CTLA-4IgV gene and functional application of vaccine in Xinjiang sheep
  105. Antitumor activity of RUNX3: Upregulation of E-cadherin and downregulation of the epithelial–mesenchymal transition in clear-cell renal cell carcinoma
  106. PHF8 promotes osteogenic differentiation of BMSCs in old rat with osteoporosis by regulating Wnt/β-catenin pathway
  107. A review of the current state of the computer-aided diagnosis (CAD) systems for breast cancer diagnosis
  108. Bilateral dacryoadenitis in adult-onset Still’s disease: A case report
  109. A novel association between Bmi-1 protein expression and the SUVmax obtained by 18F-FDG PET/CT in patients with gastric adenocarcinoma
  110. The role of erythrocytes and erythroid progenitor cells in tumors
  111. Relationship between platelet activation markers and spontaneous abortion: A meta-analysis
  112. Abnormal methylation caused by folic acid deficiency in neural tube defects
  113. Silencing TLR4 using an ultrasound-targeted microbubble destruction-based shRNA system reduces ischemia-induced seizures in hyperglycemic rats
  114. Plant Sciences
  115. Seasonal succession of bacterial communities in cultured Caulerpa lentillifera detected by high-throughput sequencing
  116. Cloning and prokaryotic expression of WRKY48 from Caragana intermedia
  117. Novel Brassica hybrids with different resistance to Leptosphaeria maculans reveal unbalanced rDNA signal patterns
  118. Application of exogenous auxin and gibberellin regulates the bolting of lettuce (Lactuca sativa L.)
  119. Phytoremediation of pollutants from wastewater: A concise review
  120. Genome-wide identification and characterization of NBS-encoding genes in the sweet potato wild ancestor Ipomoea trifida (H.B.K.)
  121. Alleviative effects of magnetic Fe3O4 nanoparticles on the physiological toxicity of 3-nitrophenol to rice (Oryza sativa L.) seedlings
  122. Selection and functional identification of Dof genes expressed in response to nitrogen in Populus simonii × Populus nigra
  123. Study on pecan seed germination influenced by seed endocarp
  124. Identification of active compounds in Ophiopogonis Radix from different geographical origins by UPLC-Q/TOF-MS combined with GC-MS approaches
  125. The entire chloroplast genome sequence of Asparagus cochinchinensis and genetic comparison to Asparagus species
  126. Genome-wide identification of MAPK family genes and their response to abiotic stresses in tea plant (Camellia sinensis)
  127. Selection and validation of reference genes for RT-qPCR analysis of different organs at various development stages in Caragana intermedia
  128. Cloning and expression analysis of SERK1 gene in Diospyros lotus
  129. Integrated metabolomic and transcriptomic profiling revealed coping mechanisms of the edible and medicinal homologous plant Plantago asiatica L. cadmium resistance
  130. A missense variant in NCF1 is associated with susceptibility to unexplained recurrent spontaneous abortion
  131. Assessment of drought tolerance indices in faba bean genotypes under different irrigation regimes
  132. The entire chloroplast genome sequence of Asparagus setaceus (Kunth) Jessop: Genome structure, gene composition, and phylogenetic analysis in Asparagaceae
  133. Food Science
  134. Dietary food additive monosodium glutamate with or without high-lipid diet induces spleen anomaly: A mechanistic approach on rat model
  135. Binge eating disorder during COVID-19
  136. Potential of honey against the onset of autoimmune diabetes and its associated nephropathy, pancreatitis, and retinopathy in type 1 diabetic animal model
  137. FTO gene expression in diet-induced obesity is downregulated by Solanum fruit supplementation
  138. Physical activity enhances fecal lactobacilli in rats chronically drinking sweetened cola beverage
  139. Supercritical CO2 extraction, chemical composition, and antioxidant effects of Coreopsis tinctoria Nutt. oleoresin
  140. Functional constituents of plant-based foods boost immunity against acute and chronic disorders
  141. Effect of selenium and methods of protein extraction on the proteomic profile of Saccharomyces yeast
  142. Microbial diversity of milk ghee in southern Gansu and its effect on the formation of ghee flavor compounds
  143. Ecology and Environmental Sciences
  144. Effects of heavy metals on bacterial community surrounding Bijiashan mining area located in northwest China
  145. Microorganism community composition analysis coupling with 15N tracer experiments reveals the nitrification rate and N2O emissions in low pH soils in Southern China
  146. Genetic diversity and population structure of Cinnamomum balansae Lecomte inferred by microsatellites
  147. Preliminary screening of microplastic contamination in different marine fish species of Taif market, Saudi Arabia
  148. Plant volatile organic compounds attractive to Lygus pratensis
  149. Effects of organic materials on soil bacterial community structure in long-term continuous cropping of tomato in greenhouse
  150. Effects of soil treated fungicide fluopimomide on tomato (Solanum lycopersicum L.) disease control and plant growth
  151. Prevalence of Yersinia pestis among rodents captured in a semi-arid tropical ecosystem of south-western Zimbabwe
  152. Effects of irrigation and nitrogen fertilization on mitigating salt-induced Na+ toxicity and sustaining sea rice growth
  153. Bioengineering and Biotechnology
  154. Poly-l-lysine-caused cell adhesion induces pyroptosis in THP-1 monocytes
  155. Development of alkaline phosphatase-scFv and its use for one-step enzyme-linked immunosorbent assay for His-tagged protein detection
  156. Development and validation of a predictive model for immune-related genes in patients with tongue squamous cell carcinoma
  157. Agriculture
  158. Effects of chemical-based fertilizer replacement with biochar-based fertilizer on albic soil nutrient content and maize yield
  159. Genome-wide identification and expression analysis of CPP-like gene family in Triticum aestivum L. under different hormone and stress conditions
  160. Agronomic and economic performance of mung bean (Vigna radiata L.) varieties in response to rates of blended NPS fertilizer in Kindo Koysha district, Southern Ethiopia
  161. Influence of furrow irrigation regime on the yield and water consumption indicators of winter wheat based on a multi-level fuzzy comprehensive evaluation
  162. Discovery of exercise-related genes and pathway analysis based on comparative genomes of Mongolian originated Abaga and Wushen horse
  163. Lessons from integrated seasonal forecast-crop modelling in Africa: A systematic review
  164. Evolution trend of soil fertility in tobacco-planting area of Chenzhou, Hunan Province, China
  165. Animal Sciences
  166. Morphological and molecular characterization of Tatera indica Hardwicke 1807 (Rodentia: Muridae) from Pothwar, Pakistan
  167. Research on meat quality of Qianhua Mutton Merino sheep and Small-tail Han sheep
  168. SI: A Scientific Memoir
  169. Suggestions on leading an academic research laboratory group
  170. My scientific genealogy and the Toronto ACDC Laboratory, 1988–2022
  171. Erratum
  172. Erratum to “Changes of immune cells in patients with hepatocellular carcinoma treated by radiofrequency ablation and hepatectomy, a pilot study”
  173. Erratum to “A two-microRNA signature predicts the progression of male thyroid cancer”
  174. Retraction
  175. Retraction of “Lidocaine has antitumor effect on hepatocellular carcinoma via the circ_DYNC1H1/miR-520a-3p/USP14 axis”
Downloaded on 8.12.2025 from https://www.degruyterbrill.com/document/doi/10.1515/biol-2022-0466/html
Scroll to top button