Next Article in Journal
Novel Inhibitor for Downstream Targeting of Transforming Growth Factor-β Signaling to Suppress Epithelial to Mesenchymal Transition and Cell Migration
Previous Article in Journal
Transcriptome Profiling Provides New Insights into the Molecular Mechanism Underlying the Sensitivity of Cotton Varieties to Mepiquat Chloride
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 9
10.3390/ijms23095048
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.)

by
Wanlu Zhang
1,2,3,
Qi Yuan
1,2,3,
Yiduo Wu
1,
Jing Zhang
1 and
Jingtao Nie
1,2,3,*
1
College of Horticulture Science, Zhejiang AF University, Hangzhou 311300, China
2
Key Laboratory of Quality and Safety Control for Subtropical Fruit and Vegetable, Ministry of Agriculture and Rural Affairs, Hangzhou 311300, China
3
Collaborative Innovation Center for Efficient and Green Production of Agriculture in Mountainous Areas of Zhejiang Province, College of Horticulture Science, Zhejiang AF University, Hangzhou 311300, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(9), 5048; https://doi.org/10.3390/ijms23095048
Submission received: 9 April 2022 / Revised: 26 April 2022 / Accepted: 29 April 2022 / Published: 2 May 2022
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Review Reports Versions Notes

Abstract

:
The NBS-LRR (NLR) gene family plays a pivotal role in regulating disease defense response in plants. Cucumber is one of the most important vegetable crops in the world, and various plant diseases, including powdery mildew (PM), cause severe losses in both cucumber productivity and quality annually. To characterize and understand the role of the CC-NBS-LRR(CNL) family of genes in disease defense response in cucumber plants, we performed bioinformatical analysis to characterize these genes systematically. We identified 33 members of the CNL gene family in cucumber plants, and they are distributed on each chromosome with chromosome 4 harboring the largest cluster of five different genes. The corresponding CNL family member varies in the number of amino acids and exons, molecular weight, theoretical isoelectric point (pI) and subcellular localization. Cis-acting element analysis of the CNL genes reveals the presence of multiple phytohormone, abiotic and biotic responsive elements in their promoters, suggesting that these genes might be responsive to plant hormones and stress. Phylogenetic and synteny analysis indicated that the CNL proteins are conserved evolutionarily in different plant species, and they can be divided into four subfamilies based on their conserved domains. MEME analysis and multiple sequence alignment showed that conserved motifs exist in the sequence of CNLs. Further DNA sequence analysis suggests that CsCNL genes might be subject to the regulation of different miRNAs upon PM infection. By mining available RNA-seq data followed by real-time quantitative PCR (qRT-PCR) analysis, we characterized expression patterns of the CNL genes, and found that those genes exhibit a temporospatial expression pattern, and their expression is also responsive to PM infection, ethylene, salicylic acid, and methyl jasmonate treatment in cucumber plants. Finally, the CNL genes targeted by miRNAs were predicted in cucumber plants. Our results in this study provided some basic information for further study of the functions of the CNL gene family in cucumber plants.

1. Introduction

Cucumber (C. sativus L.) is one of the most important vegetable crops in the world. During cucumber production, the cucumber is challenged with various serious plant diseases at different stages of development. For example, PM and downy mildew (DM) are the most common and serious diseases in cucumber plants, which lead to a significant reduction in the yield and quality of cucumbers. Therefore, one of the most desirable strategies to control these diseases is to breed resistant cultivars by cloning the resistance genes in cucumber plants. The NLR genes are the largest disease-resistant gene family in different crops by functioning to block pathogen invasion [1]. However, the function of these genes in cucumber plants has not been well documented. Therefore, it is of interest and importance to understand the function of the NLR gene family in disease response in cucumber plants.
Plants are constantly confronted by various pathogens and external adverse environment cues at different developmental stages. Therefore, plants have developed two layers of defensive mechanisms to resist the invasion of pathogens during this long-term evolutionary process [2]. PAMP-triggered immunity (PTI) is the first defensive mechanism which occurs on the plant cell membrane. It is an immunity reaction triggered by the recognition of pathogen-associated molecular patterns (PAMPs) by a variety of pattern recognition receptors (PRRs) on the plant cell membrane. However, pathogens secrete effectors in response to PTI immunity, leading some pathogens to disrupt the first line of defense of the cell membrane and enter the cell. As a stronger defense system in the cell, effector-triggered immunity (ETI) is immediately activated to prevent further invasion by pathogens. As receptors, the NLR genes directly or indirectly recognize the pathogen effectors [3,4]. At the same time, special responses are activated, including hypersensitivity responses (HR), cell death, and the accumulation of hydrogen peroxide.
The disease resistance gene (R gene) is a dominant resistant gene in plants and it can specifically detect pathogens to trigger resistance to disease [5]. The R gene is characterized by race specificity and higher efficiency. The NLR gene family is one of the largest R gene families in plants and contains nucleotide-binding sites (NBS) and leucine-rich repeat (LRR) domains [6]. They play an important role in protecting plants from various pathogens. Based on the N-terminal domain, this family can be divided into three categories: (coiled coil) CC-NBS-LRR (CNL), (toll/interleukin-1 receptor) TIR-NBS-LRR (TNL) and (resistance to powdery mildew 8, RPW8) RPW8-NBS-LRR (RNL) [7]. The sequences of this family of genes have a very important conserved domain—the NB-ARC domain (nucleotide-binding adaptor, shared by NOD-LRR proteins, APAF-1, R proteins, and the CED4 domain, also was named the NBS domain), which plays an important role in plant disease resistance and signal transduction [8,9,10,11]. The NBS domain consists of eight conserved motifs—P-LOOP (phosphate-binding loop), GLPL (Gly-Leu-Pro-Leu, also called kinase 3), RNBS-A (resistance nucleotide binding site-A), RNBS-B (resistance nucleotide binding site-B), RNBS-C (resistance nucleotide binding site-C), RNBS-D (resistance nucleotide binding site-D), kinase 2 and MHDV (Met-His-Asp-Val), and the conservation of these eight motifs is inconsistent across different plant species [12,13,14,15,16]. The C-terminal LRR structure is composed of various leucine or proline and aspartic acid amino acids that can specifically recognize proteins and participate in protein–protein interactions [17,18]. The LRR domains of NLR proteins often serve as detectors of pathogen invasion, either by directly interacting with pathogen-released effectors or by monitoring the status of effector-targeted host proteins such as RIN4 (RPM1-interacting protein 4) [19]. Upon recognition, the conformation of the NBS domain changes from an ADP-bound condensed state to an ATP-bound state exposed to the N-terminal domain, triggering a downstream hypersensitivity response that ultimately leads to apoptosis and spread of infected cells, and NBS proliferation inhibits pathogens [20].
In the NLR family, the CNL type is a class of important clades of plant disease resistant genes. Researches showed that CNL family genes can cause cell death in plants. For example, cell death was triggered in tobacco when a CNL protein (AT1G12290) was transiently expressed. However, in Arabidopsis, it was found that the Botrytis Susceptible1 Interactor (BOI) could regulate cell death by negatively regulating the level of AT1G12290 proteins [21,22]. The similar event of hybrid necrosis also occurs in cotton and wheat [23,24]. On the other hand, the CNL genes are also involved in the process of resisting invasion by different diseases in plants. For example, the protein encoded by the CNL genes was found to be resistant to Verticillium wilt in cotton and Arabidopsis [25,26]. In melon, CNL genes play a role in resistance to viral infection and aphid infection [27]. The role of the disease-resistant effects of CNL genes also exists in rice, wheat, and tomato [28,29,30,31,32]. In cucumber plants, two NLR genes, CsRSF1 and CsRSF2, were reported recently to positively regulate resistance to Sphaerotheca fuliginea [33]. Therefore, the CNL genes also play a role in disease resistance in cucumber plants. However, to what extent these genes contribute disease resistance and how they function in this process largely remain elusive.
In this study, 33 CNL genes were analyzed in cucumber plants through the bioinformatic method. The distribution, gene cluster, and characteristics of the CNL family members were analyzed. The Cis-elements of the promoters and conserved motifs of encoded proteins of these CNLs were predicted. Phylogenetic and synteny analysis of the CNL family members in different plant species indicated that this class of genes is evolutionarily conserved. Heatmap analysis using RNA-seq data indicated that the expression of CNL genes was induced by various biotic and abiotic stresses. Gene expression analysis further implied that these genes exhibit a temporospatial expression pattern, and their expression is also responsive to PM infection, ethylene, salicylic acid, and methyl jasmonate treatments in cucumber plants. Moreover, the CNL genes targeted by miRNAs were predicted. The results provide an understanding for the further study of the function of CNL family of genes in cucumber plants.

2. Results

2.1. Identification and Characterization of the CNL Family of Genes in Cucumber Plants

In total, 33 CNL genes were identified in cucumber plants using the genome of the Chinese Long cucumber, using Arabidopsis CNL genes as the query sequences [34]. Their protein sequences were verified using Pfam and NCBI databases for protein functional domains. We found that the physical chromosomal locations of the 33 CNL genes in cucumber plants are uneven (Figure 1). The majority of CNL genes in cucumber plants were located on chromosome 2 with 10 CNL genes, followed by chromosome 2 and 3 with six CNL genes. Five, two, and three CNL genes were located on chromosomes 5, 6, and 7, respectively. However, only one CNL gene was located on chromosome 1. Among the 33 CNL genes in cucumber plants, 16 genes existed in the form of five gene clusters, approximately accounting for 48.48% of the total CNL genes. These five gene clusters were located at the beginning or end of chromosomes 2, 3, 4, 5, and 7. However, no gene cluster existed on chromosomes 1 and 6 (Table S1). The largest gene cluster, consisting of five genes, was located on chromosome 4. Thus, it is likely that the probability of gene duplication on chromosome 4 was higher than that of other chromosomes in cucumber plants.
Additionally, the basic characteristics of the CNL genes and the encoded proteins were predicted (Table 1), including CDS length, the number of amino acids and exons, molecular weights, theoretical isoelectric point (pI), and subcellular localization. In all CsCNL genes, the CDS length was between 879 and 5403 bp and the number of exons ranged from one to seven. The length of the CNL proteins ranged from 293 to 1801 amino acids, and the isoelectric point was in the range from 5.57 to 9.15. In addition, the molecular weights varied significantly, ranging from 32.98491 kDa (Csa5G165310) to 207.11112 kDa (Csa2G014830). The result of subcellular localization showed that 19 (57.58%), 13 (39.39%) and 1 (3.03%) of these CNL proteins were located in the nucleus and cytoplasm, nucleus, or cytoplasm, respectively.

2.2. Cis-Acting Element Analysis of the Promoters of the CNL Genes in Cucumber Plants

The promoters of the 33 genes in the cucumber CNL family were analyzed for cis-acting elements. The results showed that more than 20 different cis-elements were found (Figure S1 and Table S2). The elements contained five kinds of plant-hormone-related elements, such as gibberellin (P-box), salicylic acid (TCA-element), abscisic acid (ABRE), methyl jasmonate (TGACG-motif, CGTCA-motif), and auxin (TGA-element). It indicated that the cucumber CNL family of genes may be participating in the regulation of phytohormone response. The promoter of each sequence contained one to five kinds of plant-hormone-related elements. There were 23, 15, 24, 27, and 17 genes that contained abscisic-acid-, auxin-, gibberellin-, methyl-jasmonate-, and salicylic-acid-responsive elements, respectively. It showed that one gene may have more than two of the same hormone-related elements, with the most numerous being the element response to methyl jasmonate (Figure S1). On the other hand, other important cis-elements also existed in these promoters, responsive to biotic and abiotic stresses, such as low temperature stress (CCGAAA), responsive elements involved in light- and drought-induced MYB binding sites (MBS); endosperm-specific negative expression (AACA_motif), meristem expression (CAT-box), mesophyll cell differentiation (HD-Zip 1), and seed-specific regulation (RY-element).

2.3. Analysis of the Phylogenetic Relationship and the Conserved Motifs of CNL Proteins in Cucumber Plants

To analyze the evolutionary relationships of the CNL gene family members, phylogenetic tree was constructed using the neighbor-joining method. The CNL proteins of cucumber, cabbage, and Arabidopsis were divided into 4 clades (Figure 2). There were five, and nine CNLs in CNL-A and CNL-B clades, respectively. The remaining 19 cucumber CNLs were all clustered in the CNL-C clade, while the CNL-D clade did not contain any cucumber CNLs. The amino acid sequences of the 33 cucumber CNL proteins were extracted for motif conservation analysis, and 20 motifs were screened out (Figure 3A, Table S3). We found that most of the amino acid sequences were arranged by order of “motif13-motif9-motif7-motif16-motif14-motif20-motif11-motif2-motif1-motif10-motif3-motif5-motif8-motif6-motif19-motif4-motif17-motif15-motif18-motif12”. P-LOOP (motif7), RNBS-A (motif16), kinase-2 (motif14), RNBS-B (motif20), RNBS-C (motif11), GLPL (motif2), RNBS-D (motif1), and MHDV (motif5) were conserved regions of NBS in cucumber plants and other species (Figure 3B, Figure S2, Table S3). However, some motifs were lost during the evolution of the CNLs. Only six sequences of CNL protein contained all the 20 motifs, and the rest of them lacked several motifs, for example, the encoded proteins of Csa5G647580 and Csa5G165310 only contained four conserved motifs. In addition, motif6, motif9, motif12, and motif15 were obviously either increased or deleted and shifted. Combined with a phylogenetic tree analysis, it was found that the similarity of sequence was closer, the more related motifs were contained (Figure 3A). At the C-terminus of the sequence, the motifs of the LRR family are mainly arranged in a tandem crossover manner. There were seven motifs belonging to the LRR family in the sequences of CNL proteins (motif4, motif5, motif8, motif12, motif15, motif17, and motif19). The results indicated that there are multiple LRR family motifs in one CNL protein. Moreover, a cucumber CNL protein containing all 20 motifs was aligned with homologous protein sequences from other species. We found that the conserved motifs in the NBS domain existed in different species (Figure 3B, Figure S3). The results showed that the CNL proteins were conserved in dicotyledonous and monocotyledonous plants.

2.4. Synteny Analysis of the CNL Genes of Cucumber

To reveal the origin and evolution of the CNL gene family members, we carried out a comparative syntenic analysis on CsCNL genes with another four plants, including melon, watermelon, tomato, and soybean (Figure 4; Table S4). The results showed that a total of 14 CsCNL genes had syntenic relationship with those in melon, and 17 corresponding orthologs were identified in melon. Meanwhile, 11 CsCNL genes showed a syntenic relationship with those in watermelon, and 11 corresponding orthologs were identified in watermelon. Among these orthologous pairs, nine CsCNL genes (Csa2G012670, Csa2G014830, Csa2G075440, Csa2G433370, Csa3G172400, Csa3G814390, Csa4G015840, Csa4G638480, and Csa5G165310) had their corresponding orthologs both in melon and watermelon, suggesting that these genes might play an important role in the evolution of the CNL genes. Moreover, three and eight pairs of orthologous CNLs were identified between cucumber and tomato, and cucumber and soybean, respectively (Figure 4; Table S4). The results showed that Csa6G490170 was related to at least two pairs of homologs in tomato and soybean. The above results showed that the CNL genes showed stronger homology among cucumber, melon, and watermelon, than that among cucumber and two other species, which corresponded to the fact that they belonged to Cucurbitaceae crops.

2.5. Expression Analysis of the CNL Family Genes in Cucumber Plants

Using the published RNA-seq data on the CuGenDB website, the expression level of the 33 CNL genes in different cucumber organs or tissues (root, stem, leaf, male flower, female flower, ovary, and tendril) was summarized (Figure 5, Table S5). The result showed that Csa7G420890 was highly expressed in roots, leaves, and female flowers. However, some genes, such as Csa2G012670 and Csa2G014830, were only highly expressed in the roots. The expression levels of some genes in all seven tissues were very low, such as Csa5G266890, Csa5G647580, and Csa6G375730.
According to the analysis of the cis-acting elements of the promoter, the expression level of the CNL family of genes might be affected by biotic and abiotic stresses. To verify whether CNL genes are involved in responding to stress, the expression level of the CNL genes induced by various stresses in previous studies was summarized using heatmaps (Figure 6). The results showed that the expression level of some genes had no obvious difference under stress by DM and PM, such as Csa2G433370, Csa5G266890, Csa5G647580, Csa5G647550, and Csa6G375730. Some CNL genes were up-regulated after infection with DM in cucumber plants, such as Csa5G647590, and Csa3G172400 (Figure 6A, Table S6). We found that the expression level of all genes in resistant materials (SSL508-28) was higher than that in susceptible materials (D8) at the same time point after PM inoculation (Figure 6B, Table S7). For example, the expression level of the Csa4G016430 gene in the SSL508-28 inbred line was higher than that in the D8 inbred line at 0 and 24 hpi (hours post-inoculation). In addition, CNL genes were also involved in abiotic stress responses. The expression of cucumber CNL genes is up-regulated under salt treatments, especially the expression of Csa2G012670, which is significantly increased (Figure 6C, Table S8). Most CNL genes are generally up-regulated when induced by low temperatures, while the expression levels of minority genes were unchanged or down-regulated (Figure 6D, Table S9). Taken together, these results suggested that CNL genes might be involved in responding to biotic and abiotic stresses.

2.6. Analysis of the Expression of CNL Genes by qRT-PCR

Based on the above results, six typical CNL genes were selected to verify their tissue-specific expression patterns by qRT-PCR (Figure 7). The results showed that these genes were all highly expressed in roots, which was consistent with the above results. Csa3G684170 had relatively high expression in other tissues, except for the stem, cotyledons, and hypocotyl. Compared with the gene expression level in cotyledons, Csa4G016460 had higher expression levels in other tissues. The expression level of Csa4G638480 was lower in leaves and fruits, but higher in other tissues. The expression level of Csa7G420890 was relatively higher in roots, stems, hypocotyls, and female flowers, and the lowest in fruits. Csa2G014830 and Csa2G435460 had the highest expression in roots and low expression in other tissues. The results indicated that CNL genes might play different roles in the growth and development of cucumber.
On the other hand, the induced expression patterns of these genes responding to PM was performed using a pair of near-isogenic lines, S1003 (resistant inbred line) and NIL(Pm5.1) (susceptible inbred line) (Figure 8) [36]. The expression of Csa4G016460, Csa3G172400, and Csa4G638480 in S1003 were significantly lower than that in NIL(Pm5.1) at 0 hpi. Except for the above genes, the expression of the other three genes (Csa2G014830, Csa2G435460 and Csa3G684170) had no difference at 0 hpi between S1003 and NIL(Pm5.1). The expression of Csa4G016460 was continuously decreased after inoculation in two lines. The expression level of Csa4G016460 in S1003 was significantly higher than that in NIL(Pm5.1) at 12, and 24 hpi, but the expression of Csa4G016460 in NIL(Pm5.1) at 48 hpi was higher. Csa2G014830, Csa2G435460, Csa3G684170, Csa3G172400, and Csa4G638480 had similar expression patterns post-inoculation. Their expression level reached the highest at 12 hpi and then decreased in S1003 and NIL(Pm5.1). It also found that the expression of Csa2G014830 was significantly lower at 12 hpi, and 24 hpi in NIL(Pm5.1) than that in S1003, and there was no difference at 48 hpi between the two lines. Except for 0 and 24 hpi, there was a significant difference in the expression level of Csa2G435460 between S1003 and NIL(Pm5.1). The expression of Csa2G435460 in S1003 had a significantly higher level than that in NIL(Pm5.1) at 12 hpi, but had an opposite expression pattern at 24 hpi. Csa3G684170 had a similar expression pattern in two materials while the expression was significantly higher in NIL(Pm5.1) than S1003 at 24 hpi. There were no significant difference in the two lines at 48 hpi. The expression patterns of Csa3G172400 and Csa4G638480 were similar between S1003 and NIL(Pm5.1). The expression levels of Csa3G172400 and Csa4G638480 were significantly lower at 0 hpi at S1003. In addition, the expression levels of Csa3G172400 and Csa4G638480 were significantly lower at 12 hpi at NIL(Pm5.1), while they were significantly higher at 24 hpi and 48 hpi at NIL(Pm5.1). Five CNL genes reached highest expression level at 12 hpi, which showed that they might be corelated with the PM resistance conferred by Csmlo1. These results showed that the CNL genes might be involved in PM resistance.

2.7. Expression Patterns of the CNL Genes Responding to Hormones in Cucumber Plants

The cis-elements of the predicted promoters of the CNL genes showed that the promoters of most CNLs contained 1–5 cis-elements responding to phytohormones. To further analyze the effect of phytohormones on the expression of CNL genes, the cucumber was treated with different exogenous hormones, which related to stress, and the expression patterns of seven CNL genes were selected to be monitored within a short period of time after hormone application.
The results showed that the expression level of seven CNL genes was changed, ranging from 0 h to 24 h after ethylene treatment (Figure 9). Specifically, it was found that expression of Csa3G684170, Csa2G014830, and Csa3G172400 was firstly decreased and then increased, and the expression reached the highest level at 12 h. The expression patterns of Csa2G420890 and Csa4G016460 were similar to the previous three genes, but the expression reached the highest level at 24 h after ethylene treatment. The gene expression levels of Csa4G638480 and Csa2G435460 were increased firstly and then decreased, and the expression reached the highest level at 12 h. It was speculated that, within 24 h of exogenous ethylene application, the CNL genes showed various expression patterns responding to ethylene, which might balance the relationship between development and the response to the environment.
Salicylic acid (SA) had been reported to respond to various disease resistance mechanisms in different plants [34,37,38]. Therefore, salicylic acid treatment was carried out on cucumber seedlings to further study the effect of salicylic acid on cucumber CNL gene expression by qRT-PCR (Figure 10). The expression pattern of the seven CNL genes showed a similar tendency responding to salicylic acid. The expression level of these genes increased continuously and then decreased back to normal. The difference in the pattern was that expression level of Csa3G684170, Csa2G014830, Csa3G172400, and Csa2G420890 reached the highest at 12 h, while the expression level of the other three genes reached the highest at 6 h after salicylic acid treatment. These results indicated that the expression of the CNL genes was induced by salicylic acid treatment in the early stages of treatment, especially at 6 h and 12 h.
Finally, aiming to analyze the effect of methyl jasmonate on the expression of CNL genes, the expression pattern of the CNL genes was monitored at 0–24 h after methyl jasmonate treatment, by qRT-PCR (Figure 11). We found that the expression levels of all genes kept relatively low at 0 h, 6 h, and 12 h after treatment, and the expression levels were largely increased at 24 h after treatment. The above results indicated that the expression of CNL genes was induced by different exogenous hormones.

2.8. Prediction Analysis of the Binding Site of CNL Genes Targeted by miRNA in Cucumber Plants

In order to determine whether the cucumber CNL genes were regulated by microRNA through a targeted binding site, the published miRNAs sequences were used to predict the binding site of CNL genes targeted by miRNA on the psRNATarget website [39,40]. The miRNA sequencing and prediction analysis has been performed in a previous study. The mature miRNA sequence was used for analysis of the binding site of the CNL genes. It turned out that 30 genes might be regulated by 33 kinds of microRNA (Table S10 and Figure 12). Among them, 17 CNL genes might be targeted by miRNA482, accounting for 17% of the CNL genes in cucumber plants. This is followed by the number of CNL genes targeted by miRNA2118, miRNA396, and miRNA156, accounting for 9%, 9%, and 8% of all CNL genes, respectively. Six, five, and five CNL genes were targeted by miRNA157, miRNA172, and miRNA395, respectively. In addition, it was shown that most of the interactions between the targeted CNLs and various miRNAs were at the transcriptional level, and only a few of them were at the translational level. We also found that several miRNAs could target one CNL gene, such as Csa4G638480, the homologous gene of AtADS1 in Arabidopsis, targeted by miRNA396, miRNA164, and miRNA390 (Table S10). Meanwhile, multiple CNL genes might be targeted by one miRNA, for example, the miRNA482 could target 17 CNL genes (Table S10).

3. Discussion

In this study, a total of 33 CNL genes were obtained from cucumber plants by the bioinformatics method (Figure 1 and Table 1), accounting for 0.135% of the whole gene numbers in cucumber plants. There were 455, 111, and 89 CNL genes in rice, tomato, and potato plants, respectively [41,42,43]. In addition, they accounted for 0.994%, 0.310%, and 0.153% of each gene numbers, which was larger than the proportion of CNL genes in cucumber plants. There were 69 CNL genes found in maize, with 0.089% proportion of the whole genes, which was lower than the proportion of CNL genes in cucumber plants [44]. In cucumber plants, 48.48% of the CNL genes were presented on the chromosome in the form of gene clusters (Table S1). This evolution phenomenon might have a positive effect on defending from various pathogens in nature. Evolutionary events, such as unequal crossover, insertion/deletion, and gene conversion, occurring in the NLR genes provided the possibility for increasing the rate of mutation and the formation of denser gene clusters on chromosomes, which has been found in lettuce and radish plants [45,46,47]. Based on the prediction of gene structure and protein physicochemical properties, we found that CsCNL genes had similar characteristics to the CNL genes in Chinese cabbage, cabbage and kiwifruit [1,13,48]. We speculated that the function of the CNL proteins may be affected by their characteristics.
Cis-elements play an important role in various life processes, such as plant growth and development, stress response, hormones response, and signal transduction. The promoters of CNL genes were analyzed for cis-acting elements (Figure S1, Table S2), and we found that many important cis-elements existed in the promoters, such as various plant-hormone-related elements, light-related elements. Among them, 16 gene members contained an element named “W-box”, which is the DNA binding site of the WRKY protein induced by salicylic acid [38]. Research suggested that it could positively regulate RPP8, a DM resistance gene, in Arabidopsis [49]. The results suggested that the CNL genes might be implicated in the defense against certain fungal diseases, and the CNL genes might also play a role in defense against various diseases in cucumber plants. To verify the roles of these hormonal elements in disease resistance, we monitored the expression patterns of the CNL genes under different hormone treatments by qRT-PCR (Figure 9, Figure 10 and Figure 11). We speculated that the expression level responding to methyl jasmonate at 24 h was particularly high due to the multiple elements related to methyl jasmonate response. These results indicated that CNL genes were induced by different hormones (ethylene, salicylic acid, and methyl jasmonate).
Motif conservation analysis was performed on the CNL gene family proteins, and 20 motifs were obtained (Figure 3, Figure S2). It was found that eight conserved motifs of the NBS domain existed in most of the proteins. In addition, motif6, motif 9, motif 12, and motif 15 were obviously either increased, or deleted and shifted, which might be related to the different genetic variation generated by natural selection in the face of the various adversities in cucumber plants (Figure 3A). However, some conserved motifs of individual CNL genes were lost during the evolution process, for example, Csa5G165310 lost 16 motifs in the sequence. Studies have reported that the loss rate of NLR genes in cucumber plants was higher than that of other plants. Because the copy number of cucumber NBS-LRR genes was affected by gene loss, this resulted in further diversity of the resistance conferred by them [50].
Synteny analysis showed that a higher synteny was maintained among cucurbit crops (Figure 4, Table S4). This result demonstrated the existence of more orthologous pairs between cucumber plants and other cucurbit crops, due to the closer relatedness of evolution. The results of tissue-specific expression analysis showed that the expression level of CNL genes was relatively high in the root, and the expression level in other tissues was varied (Figure 7). Some genes had different expression levels in different tissues. For example, Csa2G014830 has a high expression level in the roots, while its expression level is relatively low in tendrils and other tissues (Figure 4, Figure 7). Moreover, the expression levels of different genes in the same tissue were different, such as Csa4G016460 and Csa7G420890. Although the expression of some CNL genes in the leaves was relatively low, they were induced to higher levels when they were faced infection by PM and hormone treatment. This indicated the precise regulation of the CNL genes to balance the development and stress responses. Moreover, the results by qRT-PCR were slightly different from the previous studies, which might be caused by using a different cucumber inbred line. Heatmap analysis using published data found that the expression of cucumber CNL genes was also increased under abiotic stresses, such as salt stress and low temperature treatment (Figure 6). The results showed that these genes might not only confer resistance to pathogens, but also play a role in abiotic stresses. Inoculation experiments of PM pathogens were performed using one pair of NIL in cucumber plants, and gene expression at different time points after inoculation was detected by qRT-PCR (Figure 8). The expression was very low at 0 hpi, but the gene expression level of the CNL genes (such as Csa3G172400) increased rapidly after inoculation and reached the peak at 12 hpi. Previous studies have confirmed the existence of a recessive disease resistance gene in cucumber plants, Csmlo1, which could confer durable resistance to PM in cucumber plants [51]. The expression of the CsMLO1 gene reached the highest level at 12 hpi of PM, which indicated that the CsMLO1 gene played a key role in responding the attack of PM at this period post-inoculation of PM on cucumber plants. According to the induced expression results of the above CNL genes, we found that five genes have a similar expression pattern to CsMLO1. These results indicated that the functions of these genes might correlate with the Csmlo1 resistance pathway in the process of resisting the invasion of PM. The expression of these CNL genes was strictly controlled at a low level when there were no pathogens invading [51]. This was to prevent damage to plants caused by constitutive activation and spontaneous reaction, such as hypersensitivity reactions and reactive oxygen species (ROS) production [52]. On the other hand, the CNL family of genes were known as a class of dominant disease-resistance R genes, and they play a role in the plant’s defense mechanism (PTI and ETI) [4,53]. The above results indicated that the CNL genes might be involved in disease resistance pathways controlled by dominant genes, and might also be coupled with the Csmlo1 pathway that plays a role in PM resistance, or the other diseases resistance pathways. However, these assumptions still need further experiments to be carried out for confirmation. Therefore, the study of CNL genes is of great significance for cultivating durable and stable disease-resistant cucumber varieties.
The miRNAs are a class of small non-coding RNAs about 20-24nt in size, they mainly function to regulate gene expression at the post-transcription level in both plants and animals [54]. They play important roles in various life stages of plants, for example, the conserved miR156–miR172 were reported to regulate plant vegetative phase changes [55]. The miRNAs have been shown to be responsive to various biotic and abiotic stresses [56,57,58,59,60], including pathogen infection. Furthermore, miR396 negatively regulates rice blast resistance by inhibiting multiple OsGRFs genes [61], and miR482 is an ancient and conserved miRNA family of 22nt in length, and it plays an important role in the development and the process of resistance to diseases and stresses [62]. Studies have reported that miRNA482 could participate in the defense against pathogens by regulating their targeted genes, such as miR482 targets of the mRNA of CNL genes in cotton [39]. This study indicated that only 12% of NLR genes were targeted by the members of the miR482 family, and miRNA482 was inhibited when its targeted NLRs were induced to defend against pathogen invasion [39]. In the presence of pathogens, miR482 was silenced and the expression of the targeted genes NLR were increased, resulting in plant resistance. Specifically, miRNA482 cleaves the P-LOOP motif of the target NLR genes and produces a large amount of phasiRNA, which can enhance the silencing effect of miR482 on NLR genes. Thus, miR482 and its target genes NLR, as well as the generated phasiRNA, are involved in the process of plant immunity [37]. Moreover, the known miRNAs sequences were used to predict the possible target sites on cucumber CNL genes. It was found that one CNL gene could be targeted by different miRNAs, and the same miRNA could also target multiple CNL genes. Among them, miRNAs (such as miRNA156 and miRNA172) might regulate CNL genes and play an important role in growth and development. In addition, miRNA156 has been proven to function in the development of the age pathway and it is involved in responding to biotic stresses [55,63]. Importantly, miRNA482 was predicted to target 17 CNL genes in cucumber plants (Table S10). Previous studies have shown that miRNA482 could negatively regulate the CNL gene and play an important role in plant disease resistance [64]. These results indicated that miRNA482 and CNL genes might also be correlated in disease resistance in cucumber plants, and the CNL genes might be regulated by miRNA to balance the development and defense against pathogens. Therefore, the present study provided a direction for further research on the mechanism of action between miRNAs and CNL genes in cucumber plants. According to the previous studies, miRNA482 could cleave the target genes, for example, the CNL genes, and generate a large number of phasiRNAs, which were then involved in the response to pathogenic invasion in the plant’s ETI immune system. Therefore, we speculated that these CNL genes might be regulated by miRNAs, and they could cooperate with each other to play a role in disease resistance in cucumber plants.

4. Materials and Methods

4.1. Identification of Cucumber CNL Family Genes

Firstly, blast homologous sequences of CNL in cucumber plants were searched on the CuGenDB website (http://cucurbitgenomics.org/, accessed on 14 March 2022) by using the Arabidopsis CNL family of genes as the query sequences [65]. Secondly, the NB-ARC domain was identified by PF00931 on the Pfam website (http://pfam.xfam.org/, accessed on 22 February 2022), then homologous protein sequences from Arabidopsis were scanned for “hmmsearch”. Then, the candidates and conserved domains of CNLs were confirmed by a NCBI CD search (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi, accessed on 8 March 2022). In addition, the coiled-coils structure was confirmed on the Paircoil2 website (http://cb.csail.mit.edu/cb/paircoil2/paircoil2.html, accessed on 25 February 2022), and the P-value parameter was set as 0.025 [66].

4.2. Analysis of Gene Characteristics, Genomic Distribution, and Cis-Acting Elements in Promoters

Information concerning the CNL genes, including exon numbers and CDS, were retrieved from the CuGenDB website and proved using FGENESH (http://linux1.softberry.com/berry.phtml?topic=fgenesh&group=programs&subgroup=gfind, accessed on 26 April 2022). The theoretical pI and molecular weight of the identified CNLs were calculated using the ProtParam website (https://web.expasy.org/protparam/, accessed on 22 April 2022). Moreover, we predicted the subcellular localization of CNL proteins by using the CELLO website (http://cello.life.nctu.edu.tw/, accessed on 8 April 2022). The cucumber genome annotation file was downloaded from the Ensembl Plants website (http://plants.ensembl.org/index.html, accessed on 16 April 2022) for analysis. According to the gene position, TBtools was used to map the physical location of the CNL genes. A gene cluster was regarded as the distance between two adjacent NLR genes being <200 kb, with ≤8 non-NBS-LRR genes between the two NLR genes [38]. Furthermore, the promoter sequence (2 kb upstream of the gene initiation codon) of CNL genes was submitted to the Plant CARE website for cis-element analysis (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 23 April 2022), and then visualized using TBtools.

4.3. Analysis of the Conserved Motif and the Synteny of the CNL

Using the MEME Suite website (https://meme-suite.org/meme/, accessed on 25 March 2022), the amino acid sequences of 33 CNL genes were added for motif prediction. The results were downloaded in a MEME.xml format, and then be input into TBtools software for visualization. The whole genome sequence files and gene annotation files of soybean, tomato, and maize from the Ensembl Plant website were downloaded, and the whole genome sequence files and gene annotation files of cucumber, melon, pumpkin, and watermelon from the CuGenDB website were downloaded. The collinear relationship between different species was plotted with the help of TBtools software.

4.4. Transcriptome Analysis of the Genes Expression of the CNLs

The cucumber transcriptome data (PRJNA80169) was found via the cucumber genome website, and the expression levels of 33 CNL family genes in seven tissues, including roots, stems, leaves, male flowers, female flowers, carpels, and tendrils, were analyzed. The expression levels of these genes were plotted into a heatmap using TBtools software. Similarly, heatmaps of CNL gene family expression were drawn for the cucumber transcriptome data of PM treatment (PRJNA321023) [67], DM treatment (PRJNA388584), salt stress treatment (PRJNA437579), and low temperature stress treatment (PRJNA438923).

4.5. Tissue-Expression Analysis

Cucumber 9930 materials and qRT-PCR analysis were used for tissue-specific expression experiments. The roots, stems, leaves, cotyledons, hypocotyls, male flowers, female flowers, and fruits of cucumber 9930 were sampled. The experiment had three biological replicates and each organ from 15 cucumber plants were harvested and used as a sample in each replicate.

4.6. Cucumber Materials and the Treatment of PM and Hormones

Cucumber seedlings were planted in a constant temperature incubator (16 h day and eight hours night; temperature 25 °C). Before the treatment, the seedlings were cultivated in an incubator free of pathogens. Cucumber 9930 materials were used for hormone-responding experiments, and S1003 and NIL(Pm5.1) materials were used for the inoculation experiments with PM [36,68]. At the two-leaf stage, different hormones, salicylic acid (2 mmol/L), methyl jasmonate (100 μmol/L), ethylene (200 ppmol/L), and abscisic acid (100 μmol/L) were sprayed on the cucumber 9930 leaves. Each treatment had three biological replicates and the leaves of 15 plants were harvested as samples at different time points after treatments (0 h, 6 h, 12 h, and 24 h) in each replicate. At the two-leaf stage, plants were inoculated with the PM pathogen by spraying a spore suspension (1 × 105 spores/mL) evenly onto the leaves. Samples were taken at 0 h, 12 h, 24 h, and 48 h after inoculation for expression analysis. The experiment had three biological replicates and the leaves of 10 plants were harvested and used as a sample at different time points in each replicate.

4.7. Total RNA Extraction and qRT-PCR Analysis

Total RNA from different cucumber tissues was extracted by the Triozol method, the DNA was removed, and then was reversely transcribed into cDNA using the PrimeScript first Strand cDNA Synthesis Kit (Takara, Japan). The primers for qRT-PCR were designed in Primer 5 and the CsActin was used as an internal control gene (Table S11). The qRT-PCR was carried out using the SYBR Premix Ex Taq II Kit (Takara, Japan), PCR was performed on a StepOne PlusTM (ABI, America) real-time PCR machine. The qRT-PCR parameters were set as 95 °C for 10 s, 60 °C for 15 s, and 72 °C for 25 s for 45 cycles, and melting curve analysis was performed with the default settings on the instrument. The 2−ΔΔCt calculation method was used for quantitative expression analysis. Data were analyzed using normalization algorithms.

4.8. Prediction of Binding Sites of the CNL Genes Targeted by miRNA

Downloaded miRNA sequences were from psRNATarget (https://www.zhaolab.org/psRNATarget/, accessed on 9 April 2022) [69]. The cDNA library was chosen, with the library entitled “Cucumis sativus (cucumber), cds, cucumber genome sequencing project, version 2”, and the rest of the parameters were set as default. The predicted binding sites of CNL genes were screened out in the Excel sheet.

5. Conclusions

In this study, a total of 33 CNL genes were identified in cucumber plants. The distribution of CNL genes on the chromosome was uneven, some genes existed on the chromosomes in the form of gene clusters. The corresponding CNL proteins were varied in the number of amino acids and exons, molecular weight, theoretical isoelectric point (pI) and subcellular localization. Analysis of cis-elements revealed that the promoters of the CNL genes included many important elements in the growth and development of plants, such as hormone-response elements, drought-response elements, cold-response elements and damage-inducing elements. Analysis of the phylogenetic tree and conversed motifs suggested that the CNL proteins are evolutionarily conserved. Syntenic analysis indicated that more orthologous pairs existed in the CNL gene family among the Cucurbitaceae crops (cucumber, melon, and watermelon), compared to tomato and soybean plants. Heatmap analysis and tissue-specific expression analysis of the CsCNL genes demonstrated their diverse spatiotemporal expression patterns. Moreover, heatmap analysis also indicated that the expression of CNL genes were induced by various biotic and abiotic stresses. The qRT-PCR analysis showed that those genes also responded to PM infection, and treatment with ethylene, salicylic acid, and methyl jasmonate in cucumber plants. Finally, we predicted that many CNL genes were targeted by miRNAs, especially the 17 CNL genes targeted by miRNA482. Our results provided the basis for further study of the function of CNL genes in cucumber plants.

Supplementary Materials

The following supporting information can be downloaded at: www.mdpi.com/article/10.3390/ijms1697963/s1.

Author Contributions

W.Z. and J.N. planned and designed the study and wrote the manuscript. W.Z. and Q.Y. performed the experiments. W.Z. implemented the software. Y.W. and J.Z. collected the sequences of genes in this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number No. 31701915, Zhejiang Province Public Welfare Technology Application Research Project, grant number No. LGN19C150007, and a Start-Up Found from Zhejiang Agriculture & Forestry University, grant number 2017FR006.

Institutional Review Board Statement

This study did not involve humans or animals.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (No. 31701915), Zhejiang Province Public Welfare Technology Application Research Project (No. LGN19C150007), and a Start-Up Fund from Zhejiang Agriculture & Forestry University (2017FR006). Special thanks to Professor Gang Wu for his help in editing this article.

Conflicts of Interest

The authors have declared that no competing interests exist.

References

  1. Liu, Y.; Li, D.; Yang, N.; Zhu, X.; Han, K.; Gu, R.; Bai, J.; Wang, A.; Zhang, Y. Genome-Wide Identification and Analysis of CC-NBS-LRR Family in Response to Downy Mildew and Black Rot in Chinese Cabbage. Int. J. Mol. Sci. 2021, 22, 4266. [Google Scholar] [CrossRef] [PubMed]
  2. Boller, T.; He, S.Y. Innate Immunity in Plants: An Arms Race Between Pattern Recognition Receptors in Plants and Effectors in Microbial Pathogens. Science 2009, 324, 742–744. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Jones, J.D.G.; Dangl, J.L. The plant immune system. Nature 2006, 444, 323–329. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Yuan, M.; Ngou, B.P.M.; Ding, P.; Xin, X.-F. PTI-ETI crosstalk: An integrative view of plant immunity. Curr. Opin. Plant Biol. 2021, 62, 102030. [Google Scholar] [CrossRef] [PubMed]
  5. Van der Hoorn, R.A.L.; Kamoun, S. From Guard to Decoy: A New Model for Perception of Plant Pathogen Effectors. The Plant Cell 2008, 20, 2009–2017. [Google Scholar] [CrossRef] [Green Version]
  6. Yue, J.; Meyers, B.C.; Chen, J.; Tian, D.; Yang, S. Tracing the origin and evolutionary history of plant nucleotide-binding site–leucine-rich repeat (NBS-LRR) genes. New Phytologist 2012, 193, 1049–1063. [Google Scholar] [CrossRef]
  7. Jones, J.D.G.; Vance, R.E.; Dangl, J.L. Intracellular innate immune surveillance devices in plants and animals. Science 2016, 354, aaf6395. [Google Scholar] [CrossRef] [Green Version]
  8. Liu, H.; Dong, S.; Gu, F.; Liu, W.; Yang, G.; Huang, M.; Xiao, W.; Liu, Y.; Guo, T.; Wang, H.; et al. NBS-LRR Protein Pik-H4 Interacts with OsBIHD1 to Balance Rice Blast Resistance and Growth by Coordinating Ethylene-Brassinosteroid Pathway. Front. Plant Sci. 2017, 8, 127. [Google Scholar] [CrossRef] [Green Version]
  9. Cai, H.; Wang, W.; Rui, L.; Han, L.; Luo, M.; Liu, N.; Tang, D. The TIR-NBS protein TN13 associates with the CC-NBS-LRR resistance protein RPS5 and contributes to RPS5-triggered immunity in Arabidopsis. Plant J 2021, 107, 775–786. [Google Scholar] [CrossRef]
  10. Chen, Y.; Zhong, G.; Cai, H.; Chen, R.; Liu, N.; Wang, W.; Tang, D. A Truncated TIR-NBS Protein TN10 Pairs with Two Clustered TIR-NBS-LRR Immune Receptors and Contributes to Plant Immunity in Arabidopsis. Int. J. Mol. Sci. 2021, 22, 4004. [Google Scholar] [CrossRef]
  11. Inoue, H.; Hayashi, N.; Matsushita, A.; Xinqiong, L.; Nakayama, A.; Sugano, S.; Jiang, C.-J.; Takatsuji, H. Blast resistance of CC-NB-LRR protein Pb1 is mediated by WRKY45 through protein–protein interaction. Proc. Natl. Acad. Sci. USA 2013, 110, 9577–9582. [Google Scholar] [CrossRef] [Green Version]
  12. Takken, F.L.; Albrecht, M.; Tameling, W.I. Resistance proteins: Molecular switches of plant defence. Curr. Opin. Plant Biol. 2006, 9, 383–390. [Google Scholar] [CrossRef] [PubMed]
  13. Liu, Z.; Xie, J.; Wang, H.; Zhong, X.; Li, H.; Yu, J.; Kang, J. Identification and expression profiling analysis of NBS–LRR genes involved in Fusarium oxysporum f.sp. conglutinans resistance in cabbage. 3 Biotech 2019, 9, 1–12. [Google Scholar] [CrossRef]
  14. Miller, R.N.G.; Bertioli, D.J.; Baurens, F.C.; Santos, C.M.R.; Alves, P.C.; Martins, N.F.; Togawa, R.C.; Souza Júnior, M.T.; Pappas Júnior, G.J. Analysis of non-TIR NBS-LRR resistance gene analogs in Musa acuminata Colla: Isolation, RFLP marker development, and physical mapping. BMC Plant Biology 2008, 8, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Saraste, M.; Sibbald, P.R.; Wittinghofer, A. The P-loop—A common motif in ATP- and GTP-binding proteins. Trends Biochem. Sci. 1990, 15, 430–434. [Google Scholar] [CrossRef]
  16. Meyers, B.C.; Dickerman, A.W.; Michelmore, R.W.; Sivaramakrishnan, S.; Sobral, B.W.; Young, N.D. Plant disease resistance genes encode members of an ancient and diverse protein family within the nucleotide-binding superfamily. Plant J. 1999, 20, 317–332. [Google Scholar] [CrossRef]
  17. Lozano, R.; Ponce, O.P.; Ramírez, M.; Mostajo, N.; Orjeda, G. Genome-Wide Identification and Mapping of NBS-Encoding Resistance Genes in Solanum tuberosum Group Phureja. PLoS ONE 2012, 7, e34775. [Google Scholar] [CrossRef]
  18. Wulff, B.B.; Thomas, C.M.; Smoker, M.; Grant, M.; Jones, J.D. Domain swapping and gene shuffling identify sequences required for induction of an Avr-dependent hypersensitive response by the tomato Cf-4 and Cf-9 proteins. Plant Cell 2001, 13, 255–272. [Google Scholar] [CrossRef] [Green Version]
  19. Ray, S.K.; Macoy, D.M.; Kim, W.-Y.; Lee, S.Y.; Kim, M.G. Role of RIN4 in Regulating PAMP-Triggered Immunity and Effector-Triggered Immunity: Current Status and Future Perspectives. Mol. Cells 2019, 42, 503–511. [Google Scholar] [CrossRef]
  20. Andersen, E.J.; Ali, S.; Byamukama, E.; Yen, Y.; Nepal, M.P. Disease Resistance Mechanisms in Plants. Genes 2018, 9, 339. [Google Scholar] [CrossRef] [Green Version]
  21. Huang, J.; Wu, X.; Sun, K.; Gao, Z. Structure and function analysis of a CC-NBS-LRR protein AT1G12290. Biochem. Biophys. Res. Commun. 2021, 534, 206–211. [Google Scholar] [CrossRef] [PubMed]
  22. Huang, J.; Wu, X.; Gao, Z. The RING-type protein BOI negatively regulates the protein level of a CC-NBS-LRR in Arabidopsis. Biochem. Biophys. Res. Commun. 2021, 578, 104–109. [Google Scholar] [CrossRef] [PubMed]
  23. Si, Y.; Zheng, S.; Niu, J.; Tian, S.; Gu, M.; Lu, Q.; He, Y.; Zhang, J.; Shi, X.; Li, Y.; et al. Ne2, a typical CC–NBS–LRR-type gene, is responsible for hybrid necrosis in wheat. New Phytol. 2021, 232, 279–289. [Google Scholar] [CrossRef] [PubMed]
  24. Deng, J.; Fang, L.; Zhu, X.; Zhou, B.; Zhang, T. A CC-NBS-LRR gene induces hybrid lethality in cotton. J. Exp. Bot. 2019, 70, 5145–5156. [Google Scholar] [CrossRef]
  25. Li, T.; Zhang, Q.; Jiang, X.; Li, R.; Dhar, N. Cotton CC-NBS-LRR Gene GbCNL130 Confers Resistance to Verticillium Wilt Across Different Species. Front. Plant Sci. 2021, 12, 695691. [Google Scholar] [CrossRef]
  26. Li, N.-Y.; Zhou, L.; Zhang, D.-D.; Klosterman, S.J.; Li, T.-G.; Gui, Y.-J.; Kong, Z.-Q.; Ma, X.-F.; Short, D.P.G.; Zhang, W.-Q.; et al. Heterologous Expression of the Cotton NBS-LRR Gene GbaNA1 Enhances Verticillium Wilt Resistance in Arabidopsis. Front. Plant Sci. 2018, 9, 119. [Google Scholar] [CrossRef] [Green Version]
  27. Dogimont, C.; Chovelon, V.; Pauquet, J.; Boualem, A.; Bendahmane, A. The Vat locus encodes for a CC-NBS-LRR protein that confers resistance to Aphis gossypii infestation and A. gossypii-mediated virus resistance. Plant J. 2014, 80, 993–1004. [Google Scholar] [CrossRef]
  28. Ma, J.; Lei, C.; Xu, X.; Hao, K.; Wang, J.; Cheng, Z.; Ma, X.; Zhou, K.; Zhang, X.; Guo, X.; et al. Pi64, Encoding a Novel CC-NBS-LRR Protein, Confers Resistance to Leaf and Neck Blast in Rice. Mol. Plant-Microbe Interact. 2015, 28, 558–568. [Google Scholar] [CrossRef] [Green Version]
  29. Yang, J.; Ji, L.; Zhu, B.; Yuan, X.; Jin, D.; Xie, G. OsCML16 interacts with a novel CC-NBS-LRR protein OsPi304 in the Ca(2+)/Mg(2+) dependent and independent manner in rice. Biochem Biophys Res Commun 2018, 504, 346–351. [Google Scholar] [CrossRef]
  30. Wang, Q.; Li, Y.; Kosami, K.; Liu, C.; Li, J.; Zhang, D.; Miki, D.; Kawano, Y. Three highly conserved hydrophobic residues in the predicted α2-helix of rice NLR protein Pit contribute to its localization and immune induction. Plant Cell Environ. 2022. [Google Scholar] [CrossRef]
  31. He, H.; Zhu, S.; Zhao, R.; Jiang, Z.; Ji, Y.; Ji, J.; Qiu, D.; Li, H.; Bie, T. Pm21, Encoding a Typical CC-NBS-LRR Protein, Confers Broad-Spectrum Resistance to Wheat Powdery Mildew Disease. Mol. Plant 2018, 11, 879–882. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Zhang, C.; Liu, L.; Wang, X.; Vossen, J.; Li, G.; Li, T.; Zheng, Z.; Gao, J.; Guo, Y.; Visser, R.G.F.; et al. The Ph-3 gene from Solanum pimpinellifolium encodes CC-NBS-LRR protein conferring resistance to Phytophthora infestans. Theor. Appl. Genet. 2014, 127, 1353–1364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Wang, X.; Chen, Q.; Huang, J.; Meng, X.; Cui, N.; Yu, Y.; Fan, H. Nucleotide-Binding Leucine-Rich Repeat Genes CsRSF1 and CsRSF2 Are Positive Modulators in the Cucumis sativus Defense Response to Sphaerotheca fuliginea. Int. J. Mol. Sci. 2021, 22, 3986. [Google Scholar] [CrossRef] [PubMed]
  34. Huang, S.; Li, R.; Zhang, Z.; Li, L.; Gu, X.; Fan, W.; Lucas, W.J.; Wang, X.; Xie, B.; Ni, P.; et al. The genome of the cucumber, Cucumis sativus L. Nat Genet 2009, 41, 1275–1281. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Li, Z.; Zhang, Z.; Yan, P.; Huang, S.; Fei, Z.; Lin, K. RNA-Seq improves annotation of protein-coding genes in the cucumber genome. BMC Genom. 2011, 12, 540. [Google Scholar] [CrossRef] [Green Version]
  36. Nie, J.; Wang, Y.; He, H.; Guo, C.; Zhu, W.; Pan, J.; Li, D.; Lian, H.; Pan, J.; Cai, R. Loss-of-Function Mutations in CsMLO1 Confer Durable Powdery Mildew Resistance in Cucumber (Cucumis sativus L.). Front Plant Sci 2015, 6, 1155. [Google Scholar] [CrossRef]
  37. Zhang, Y.; Waseem, M.; Zeng, Z.; Xu, J.; Chen, C.; Liu, Y.; Zhai, J.; Xia, R. MicroRNA482/2118, a miRNA superfamily essential for both disease resistance and plant development. New Phytol. 2021, 233, 2047–2057. [Google Scholar] [CrossRef]
  38. Richly, E.; Kurth, J.; Leister, D. thaliana Evolution. Mol. Biol. Evol. 2002, 19, 76–84. [Google Scholar] [CrossRef] [Green Version]
  39. Zhu, Q.-H.; Fan, L.; Liu, Y.; Xu, H.; Llewellyn, D.; Wilson, I. miR482 Regulation of NBS-LRR Defense Genes during Fungal Pathogen Infection in Cotton. PLoS ONE 2013, 8, e84390. [Google Scholar] [CrossRef] [Green Version]
  40. Nie, J.; Wang, H.; Zhang, W.; Teng, X.; Yu, C.; Cai, R.; Wu, G. Characterization of lncRNAs and mRNAs Involved in Powdery Mildew Resistance in Cucumber. Phytopathology 2021, 111, 1613–1624. [Google Scholar] [CrossRef]
  41. Zhou, T.; Wang, Y.; Chen, J.-Q.; Araki, H.; Jing, Z.; Jiang, K.; Shen, J.; Tian, D. Genome-wide identification of NBS genes in japonica rice reveals significant expansion of divergent non-TIR NBS-LRR genes. Mol. Genet. Genom. 2004, 271, 402–415. [Google Scholar] [CrossRef] [PubMed]
  42. Andolfo, G.; Sanseverino, W.; Rombauts, S.; Peer, Y.; Bradeen, J.M.; Carputo, D.; Frusciante, L.; Ercolano, M.R. Overview of tomato (Solanum lycopersicum) candidate pathogen recognition genes reveals importantSolanum R locus dynamics. New Phytol. 2012, 197, 223–237. [Google Scholar] [CrossRef] [PubMed]
  43. Jupe, F.; Pritchard, L.; Etherington, G.J.; MacKenzie, K.; Cock, P.J.; Wright, F.; Sharma, S.K.; Bolser, D.; Bryan, G.J.; Jones, J.D.; et al. Identification and localisation of the NB-LRR gene family within the potato genome. BMC Genom. 2012, 13, 75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Cheng, Y.; Li, X.; Jiang, H.; Ma, W.; Miao, W.; Yamada, T.; Zhang, M. Systematic analysis and comparison of nucleotide-binding site disease resistance genes in maize. FEBS J. 2012, 279, 2431–2443. [Google Scholar] [CrossRef] [PubMed]
  45. Kuang, H.; Woo, S.-S.; Meyers, B.; Nevo, E.; Michelmore, R.W. Multiple Genetic Processes Result in Heterogeneous Rates of Evolution within the Major Cluster Disease Resistance Genes in Lettuce. Plant Cell 2004, 16, 2870–2894. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Chin, D.B.; Arroyo-Garcia, R.; Ochoa, O.E.; Kesseli, R.V.; Lavelle, D.O.; Michelmore, R.W. Recombination and spontaneous mutation at the major cluster of resistance genes in lettuce (Lactuca sativa). Genetics 2001, 157, 831–849. [Google Scholar] [CrossRef]
  47. Ma, Y.; Chhapekar, S.S.; Lu, L.; Oh, S.; Singh, S.; Kim, C.S.; Kim, S.; Choi, G.J.; Lim, Y.P.; Choi, S.R. Genome-wide identification and characterization of NBS-encoding genes in Raphanus sativus L. and their roles related to Fusarium oxysporum resistance. BMC Plant Biol. 2021, 21, 1–17. [Google Scholar] [CrossRef]
  48. Wang, T.; Jia, Z.-H.; Zhang, J.-Y.; Liu, M.; Guo, Z.-R.; Wang, G. Identification and Analysis of NBS-LRR Genes in Actinidia chinensis Genome. Plants 2020, 9, 1350. [Google Scholar] [CrossRef]
  49. Mohr, T.J.; Mammarella, N.D.; Hoff, T.; Woffenden, B.J.; Jelesko, J.G.; McDowell, J.M. The Arabidopsis Downy Mildew Resistance Gene RPP8 Is Induced by Pathogens and Salicylic Acid and Is Regulated by W Box cis Elements. Mol. Plant-Microbe Interact. 2010, 23, 1303–1315. [Google Scholar] [CrossRef] [Green Version]
  50. Lin, X.; Zhang, Y.; Kuang, H.; Chen, J. Frequent loss of lineages and deficient duplications accounted for low copy number of disease resistance genes in Cucurbitaceae. BMC Genom. 2013, 14, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Korves, T.; Bergelson, J. A novel cost of R gene resistance in the presence of disease. Am. Nat. 2004, 163, 489–504. [Google Scholar] [CrossRef]
  52. Ariga, H.; Katori, T.; Tsuchimatsu, T.; Hirase, T.; Tajima, Y.; Parker, J.E.; Alcazar, R.; Koornneef, M.; Hoekenga, O.; Lipka, A.E.; et al. NLR locus-mediated trade-off between abiotic and biotic stress adaptation in Arabidopsis. Nat. Plants 2017, 3, 17072. [Google Scholar] [CrossRef] [Green Version]
  53. Yuan, M.; Jiang, Z.; Bi, G.; Nomura, K.; Liu, M.; Wang, Y.; Cai, B.; Zhou, J.-M.; He, S.Y.; Xin, X.-F. Pattern-recognition receptors are required for NLR-mediated plant immunity. Nature 2021, 592, 1–5. [Google Scholar] [CrossRef] [PubMed]
  54. Wu, G. Plant MicroRNAs and Development. J. Genet. Genom. 2013, 40, 217–230. [Google Scholar] [CrossRef]
  55. Wu, G.; Park, M.Y.; Conway, S.R.; Wang, J.-W.; Weigel, D.; Poethig, R.S. The Sequential Action of miR156 and miR172 Regulates Developmental Timing in Arabidopsis. Cell 2009, 138, 750–759. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Chen, S.; Wu, J.; Zhang, Y.; Zhao, Y.; Xu, W.; Li, Y.; Xie, J. Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164-NAC-mRNA Regulatory Pathway for Disease Defense in Populus. Front. Genet. 2021, 12, 668940. [Google Scholar] [CrossRef] [PubMed]
  57. Yao, S.; Yang, Z.; Yang, R.; Huang, Y.; Guo, G.; Kong, X.; Lan, Y.; Zhou, T.; Wang, H.; Wang, W.; et al. Transcriptional Regulation of miR528 by OsSPL9 Orchestrates Antiviral Response in Rice. Mol. Plant 2019, 12, 1114–1122. [Google Scholar] [CrossRef] [PubMed]
  58. Zhou, Z.; Cao, Y.; Li, T.; Wang, X.; Chen, J.; He, H.; Yao, W.; Wu, J.; Zhang, H. MicroRNAs Are Involved in Maize Immunity Against Fusarium verticillioides Ear Rot. Genom. Proteom. Bioinform. 2020, 18, 241–255. [Google Scholar] [CrossRef]
  59. Xu, T.; Zhang, L.; Yang, Z.; Wei, Y.; Dong, T. Identification and Functional Characterization of Plant MiRNA Under Salt Stress Shed Light on Salinity Resistance Improvement Through MiRNA Manipulation in Crops. Front. Plant Sci. 2021, 12, 665439. [Google Scholar] [CrossRef]
  60. Zhang, L.; Ding, H.; Jiang, H.; Wang, H.; Chen, K.; Duan, J.; Feng, S.; Wu, G. Regulation of cadmium tolerance and accumulation by miR156 in Arabidopsis. Chemosphere 2020, 242, 125168. [Google Scholar] [CrossRef] [PubMed]
  61. Chandran, V.; Wang, H.; Gao, F.; Cao, X.-L.; Chen, Y.-P.; Li, G.-B.; Zhu, Y.; Yang, X.-M.; Zhang, L.-L.; Zhao, Z.; et al. miR396-OsGRFs Module Balances Growth and Rice Blast Disease-Resistance. Front. Plant Sci. 2019, 9, 1999. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. de Vries, S.; Kukuk, A.; von Dahlen, J.K.; Schnake, A.; Kloesges, T.; Rose, L.E. Expression profiling across wild and cultivated tomatoes supports the relevance of early miR482/2118 suppression for Phytophthora resistance. Proc. Biol. Sci. 2018, 285, 1873. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Ge, Y.; Han, J.; Zhou, G.; Xu, Y.; Ding, Y.; Shi, M.; Guo, C.; Wu, G. Silencing of miR156 confers enhanced resistance to brown planthopper in rice. Planta 2018, 248, 813–826. [Google Scholar] [CrossRef]
  64. Shivaprasad, P.V.; Chen, H.-M.; Patel, K.; Bond, D.; Santos, B.; Baulcombe, D.C. A MicroRNA Superfamily Regulates Nucleotide Binding Site–Leucine-Rich Repeats and Other mRNAs. Plant Cell 2012, 24, 859–874. [Google Scholar] [CrossRef] [Green Version]
  65. Meyers, B.C.; Kozik, A.; Griego, A.; Kuang, H.; Michelmore, R.W. Genome-Wide Analysis of NBS-LRR–Encoding Genes in Arabidopsis. Plant Cell 2003, 15, 809–834. [Google Scholar] [CrossRef] [Green Version]
  66. Fiumara, F.; Fioriti, L.; Kandel, E.R.; Hendrickson, W.A. Essential Role of Coiled Coils for Aggregation and Activity of Q/N-Rich Prions and PolyQ Proteins. Cell 2010, 143, 1121–1135. [Google Scholar] [CrossRef] [Green Version]
  67. Xu, Q.; Xu, X.; Shi, Y.; Qi, X.; Chen, X. Elucidation of the molecular responses of a cucumber segment substitution line carrying Pm5.1 and its recurrent parent triggered by powdery mildew by comparative transcriptome profiling. BMC Genom. 2017, 18, 1–14. [Google Scholar] [CrossRef] [Green Version]
  68. Nie, J.; He, H.; Peng, J.; Yang, X.; Wang, Y.; Bie, B.; Zhao, J.; Wang, Y.; Si, L.; Pan, J.; et al. Identification and fine mapping of pm5.1: A recessive gene for powdery mildew resistance in cucumber (Cucumis sativus L.). Mol. Breed. 2015, 35, 7. [Google Scholar] [CrossRef]
  69. Dai, X.; Zhuang, Z.; Zhao, P.X. psRNATarget: A plant small RNA target analysis server (2017 release). Nucleic Acids Res. 2018, 46, W49–W54. [Google Scholar] [CrossRef] [Green Version]
Ijms 23 05048 g001 550
Figure 1. The distribution of the CNL genes located on the chromosomes in cucumber plants. The genetic distance of seven chromosomes were represented by the scale in megabases (Mb) on the left. The CNL genes are displayed using nomenclature for genome version 2 of Chinese Long cucumber. Blue lines represent the location of the gene on each chromosome.
Figure 1. The distribution of the CNL genes located on the chromosomes in cucumber plants. The genetic distance of seven chromosomes were represented by the scale in megabases (Mb) on the left. The CNL genes are displayed using nomenclature for genome version 2 of Chinese Long cucumber. Blue lines represent the location of the gene on each chromosome.
Ijms 23 05048 g001
Ijms 23 05048 g002 550
Figure 2. Phylogenetic relationships of CNL proteins among cucumber, cabbage and Arabidopsis. The unrooted phylogenetic tree was constructed by MEGA7.0 by neighbor-joining method with 1000 bootstrap replicates. The CNLs are divided into four major subfamilies. Different subfamily was indicated with different colors. The black triangles represented cabbage, black squares represented Arabidopsis, and black circles represented cucumber.
Figure 2. Phylogenetic relationships of CNL proteins among cucumber, cabbage and Arabidopsis. The unrooted phylogenetic tree was constructed by MEGA7.0 by neighbor-joining method with 1000 bootstrap replicates. The CNLs are divided into four major subfamilies. Different subfamily was indicated with different colors. The black triangles represented cabbage, black squares represented Arabidopsis, and black circles represented cucumber.
Ijms 23 05048 g002
Ijms 23 05048 g003 550
Figure 3. Analysis of the conserved motifs of CNL proteins. (A) MEME analysis of the CNL proteins in cucumber plants. (B) Alignment of the conserved NBS domain of CNL proteins between cucumber and other species. Black color: the homolog level was 100%; Pink color: the homolog level was greater or equal to 75%; Blue color: the homolog level was greater or equal to 50%. Maize, XP_008666515.1; bitter melon, XP_022145177.1; melon, XP_008453965.1; tomato, XP_004245923.1; soybean, XP_006599131.1; pumpkin, CmaCh14G014850.1; Arabidopsis, At1G10920; watermelon, Cla97C07G142810.1; and rice, XP_015651137.1.
Figure 3. Analysis of the conserved motifs of CNL proteins. (A) MEME analysis of the CNL proteins in cucumber plants. (B) Alignment of the conserved NBS domain of CNL proteins between cucumber and other species. Black color: the homolog level was 100%; Pink color: the homolog level was greater or equal to 75%; Blue color: the homolog level was greater or equal to 50%. Maize, XP_008666515.1; bitter melon, XP_022145177.1; melon, XP_008453965.1; tomato, XP_004245923.1; soybean, XP_006599131.1; pumpkin, CmaCh14G014850.1; Arabidopsis, At1G10920; watermelon, Cla97C07G142810.1; and rice, XP_015651137.1.
Ijms 23 05048 g003
Ijms 23 05048 g004 550
Figure 4. Synteny of the CNL genes among cucumber and other species. Gray lines in the background represent collinear blocks in cucumber plants and other genomes. The collinear gene pairs with CNL genes between different species were highlighted by the blue lines. Red inverted triangle indicated the locations of the CNL genes. Cucumber (C. sativus L.), melon (Cucumis melo L.), watermelon (Citrullus lanatus L.), tomato (Solanum lycopersicum L.), and soybean (Glycine max L.).
Figure 4. Synteny of the CNL genes among cucumber and other species. Gray lines in the background represent collinear blocks in cucumber plants and other genomes. The collinear gene pairs with CNL genes between different species were highlighted by the blue lines. Red inverted triangle indicated the locations of the CNL genes. Cucumber (C. sativus L.), melon (Cucumis melo L.), watermelon (Citrullus lanatus L.), tomato (Solanum lycopersicum L.), and soybean (Glycine max L.).
Ijms 23 05048 g004
Ijms 23 05048 g005 550
Figure 5. Tissue-specific expression of the CNL genes in cucumber plants. The transcriptional levels of CsCNL genes in seven tissues of cucumber 9930 were investigated based on public transcriptome data (PRJNA80169) [35]. The genome-wide expression of CsCNL genes were shown on a heatmap using a log2RPKM value, and −6.00 to 8.00 was artificially set with the color scale limits according to the normalized value. The color scale showed increasing expression levels from blue to red.
Figure 5. Tissue-specific expression of the CNL genes in cucumber plants. The transcriptional levels of CsCNL genes in seven tissues of cucumber 9930 were investigated based on public transcriptome data (PRJNA80169) [35]. The genome-wide expression of CsCNL genes were shown on a heatmap using a log2RPKM value, and −6.00 to 8.00 was artificially set with the color scale limits according to the normalized value. The color scale showed increasing expression levels from blue to red.
Ijms 23 05048 g005
Ijms 23 05048 g006 550
Figure 6. The expression patterns of the CNL genes under biotic and abiotic stresses in cucumber plants. (A) The expression of cucumber CNL genes in response to DM, post-inoculation at 0 h, 6 h, and 24 h. (B) The expression of cucumber CNL genes in response to PM, post-inoculation at 0 h, and 48 h. D8_CT, D8 line control; D8_PM, D8 line, 48 h post-inoculation with PM; SSL508-28_CT, SSL508-28 line as control; SSL508-28_PM, SSL508-28 line, 48 h post-inoculation with PM. (C) and (D) The expression of cucumber CNL genes in response to salt stress and chilling stresses. Hpi, hours post-inoculation. The genome-wide expression of CsCNL genes were shown on a heatmap using log2RPKM value, and −8.00 to 6.00, −9.00 to 9.00, −10.00 to 6.00 were artificially set with the color scale limits according to the normalized value. The color scale was shown increasing expression levels from blue to red.
Figure 6. The expression patterns of the CNL genes under biotic and abiotic stresses in cucumber plants. (A) The expression of cucumber CNL genes in response to DM, post-inoculation at 0 h, 6 h, and 24 h. (B) The expression of cucumber CNL genes in response to PM, post-inoculation at 0 h, and 48 h. D8_CT, D8 line control; D8_PM, D8 line, 48 h post-inoculation with PM; SSL508-28_CT, SSL508-28 line as control; SSL508-28_PM, SSL508-28 line, 48 h post-inoculation with PM. (C) and (D) The expression of cucumber CNL genes in response to salt stress and chilling stresses. Hpi, hours post-inoculation. The genome-wide expression of CsCNL genes were shown on a heatmap using log2RPKM value, and −8.00 to 6.00, −9.00 to 9.00, −10.00 to 6.00 were artificially set with the color scale limits according to the normalized value. The color scale was shown increasing expression levels from blue to red.
Ijms 23 05048 g006
Ijms 23 05048 g007 550
Figure 7. Tissue-specific expression analysis of the CNL genes in cucumber plants by qRT-PCR. R: root; S: stem; C: cotyledon; H: hypocotyl; L: leaf; M: male flower; Fe: female flower; and Fr: fruit. The vertical axis is relative to the expression level and x-axis represents different tissues. Values are mean ± SE (n = 3).
Figure 7. Tissue-specific expression analysis of the CNL genes in cucumber plants by qRT-PCR. R: root; S: stem; C: cotyledon; H: hypocotyl; L: leaf; M: male flower; Fe: female flower; and Fr: fruit. The vertical axis is relative to the expression level and x-axis represents different tissues. Values are mean ± SE (n = 3).
Ijms 23 05048 g007
Ijms 23 05048 g008 550
Figure 8. The expression patterns of the CNL genes responding to PM in cucumber plants. The expression levels of the CNL genes were detected in S1003 and NIL(Pm5.1) for 0 h, 12 h, 24 h, and 48 h after inoculation. S1003, resistance cucumber inbred line; NIL, and susceptible cucumber inbred line NIL(Pm5.1). The expression level of transcript at 0 h was set to a value of ‘1′. Values are mean ± SE (n = 3) (* and ** indicate significant differences between S1003 and NIL(Pm5.1) at p = 0.01 and 0.05, respectively).
Figure 8. The expression patterns of the CNL genes responding to PM in cucumber plants. The expression levels of the CNL genes were detected in S1003 and NIL(Pm5.1) for 0 h, 12 h, 24 h, and 48 h after inoculation. S1003, resistance cucumber inbred line; NIL, and susceptible cucumber inbred line NIL(Pm5.1). The expression level of transcript at 0 h was set to a value of ‘1′. Values are mean ± SE (n = 3) (* and ** indicate significant differences between S1003 and NIL(Pm5.1) at p = 0.01 and 0.05, respectively).
Ijms 23 05048 g008
Ijms 23 05048 g009 550
Figure 9. Expression analysis of the CNL genes responding to ethylene treatment in cucumber plants. The vertical axis is relative expression level and 0, 6, 12, 24 h on the x-axis indicate the treatment time. The expression level of transcript at 0 h was set to a value of ‘1′. Values are mean ± SE (n = 3).
Figure 9. Expression analysis of the CNL genes responding to ethylene treatment in cucumber plants. The vertical axis is relative expression level and 0, 6, 12, 24 h on the x-axis indicate the treatment time. The expression level of transcript at 0 h was set to a value of ‘1′. Values are mean ± SE (n = 3).
Ijms 23 05048 g009
Ijms 23 05048 g010 550
Figure 10. Expression analysis of the CNL genes responding to salicylic acid treatment in cucumber plants. The vertical axis is relative expression level and 0, 6, 12, 24 h on the x-axis indicate the treatment time. The expression level of transcript at 0 h was set to a value of ‘1′. Values are mean ± SE (n = 3).
Figure 10. Expression analysis of the CNL genes responding to salicylic acid treatment in cucumber plants. The vertical axis is relative expression level and 0, 6, 12, 24 h on the x-axis indicate the treatment time. The expression level of transcript at 0 h was set to a value of ‘1′. Values are mean ± SE (n = 3).
Ijms 23 05048 g010
Ijms 23 05048 g011 550
Figure 11. Expression analysis of the CNL genes responding to methyl jasmonate in cucumber plants. The vertical axis is relative expression level and 0, 6, 12, 24 h on the x-axis indicate the treatment time. The expression level of transcript at 0 h was set to a value of ‘1′. Values are mean ± SE (n = 3).
Figure 11. Expression analysis of the CNL genes responding to methyl jasmonate in cucumber plants. The vertical axis is relative expression level and 0, 6, 12, 24 h on the x-axis indicate the treatment time. The expression level of transcript at 0 h was set to a value of ‘1′. Values are mean ± SE (n = 3).
Ijms 23 05048 g011
Ijms 23 05048 g012 550
Figure 12. The predicted binding site of the CNL genes targeted by miRNAs in cucumber plants. Two dots indicated paired successfully between bases, and one dot indicated that there is also a pairing between U and G in the secondary structure. Blank space indicated that two bases failed to be paired.
Figure 12. The predicted binding site of the CNL genes targeted by miRNAs in cucumber plants. Two dots indicated paired successfully between bases, and one dot indicated that there is also a pairing between U and G in the secondary structure. Blank space indicated that two bases failed to be paired.
Ijms 23 05048 g012
Table
Table 1. Characteristics of the CNL family of genes and the corresponding proteins in cucumber plants.
Table 1. Characteristics of the CNL family of genes and the corresponding proteins in cucumber plants.
Gene IDLength of CDS (bp)Number of Amino AcidsNumber of ExonsMolecular Weight (kDa)Theoretical pISubcellular Localization
Csa1G201260
Csa2G008000		2409
2469		803
823		7
5		90.94932
95.32828		5.72
6.80		N/Cyto
Cyto/N
Csa2G0126702439813592.460417.93N
Csa2G014830540318017207.111125.57N
Csa2G07413027459152105.955858.29N
Csa2G075440405613524152.929066.35N/Cyto
Csa2G09693029139711110.913796.20N/Cyto
Csa2G403680396613224151.811256.54N
Csa2G43334028989662111.929628.38N/Cyto
Csa2G433370357911934136.252346.98N
Csa2G43546027219072104.056117.23N
Csa3G172400344411481131.344137.38N/Cyto
Csa3G684170326710891125.618996.10Cyto/N
Csa3G8143902409803193.331959.15N
Csa3G8144002466822194.426097.10N
Csa3G81540026768921103.116336.43N
Csa3G822360315910532120.249546.27N
Csa4G015840326110871124.356656.84Cyto/N
Csa4G015850324010801123.310247.36N/Cyto
Csa4G01636028749582109.571717.09N
Csa4G0164302520840396.126706.14Cyto/N
Csa4G01646028419472108.232337.33N
Csa4G6384802463821594.246556.58N
Csa5G165310879293132.984916.80Cyto
Csa5G266890314110471119.37916.04Cyto/N
Csa5G647550321310716121.786726.70N
Csa5G6475801143381344.469308.70N/Cyto
Csa5G6475901884628771.539995.97N
Csa6G37573027519172106.784558.21N
Csa6G4901702454818593.590076.05N
Csa7G239020362112072138.190077.59N/Cyto
Csa7G4208901665555264.716857.61N
Csa7G425940306010205116.400686.86N
N represents nucleus; Cyto represents cytoplasm.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, W.; Yuan, Q.; Wu, Y.; Zhang, J.; Nie, J. Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.). Int. J. Mol. Sci. 2022, 23, 5048. https://doi.org/10.3390/ijms23095048

AMA Style

Zhang W, Yuan Q, Wu Y, Zhang J, Nie J. Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.). International Journal of Molecular Sciences. 2022; 23(9):5048. https://doi.org/10.3390/ijms23095048

Chicago/Turabian Style

Zhang, Wanlu, Qi Yuan, Yiduo Wu, Jing Zhang, and Jingtao Nie. 2022. "Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.)" International Journal of Molecular Sciences 23, no. 9: 5048. https://doi.org/10.3390/ijms23095048

APA Style

Zhang, W., Yuan, Q., Wu, Y., Zhang, J., & Nie, J. (2022). Genome-Wide Identification and Characterization of the CC-NBS-LRR Gene Family in Cucumber (Cucumis sativus L.). International Journal of Molecular Sciences, 23(9), 5048. https://doi.org/10.3390/ijms23095048

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