Next Article in Journal
α-Lipoic Acid Reduces Iron-induced Toxicity and Oxidative Stress in a Model of Iron Overload
Next Article in Special Issue
Molecular Mechanisms Underlying Increase in Lysine Content of Waxy Maize through the Introgression of the opaque2 Allele
Previous Article in Journal
Cirsiliol Suppressed Epithelial to Mesenchymal Transition in B16F10 Malignant Melanoma Cells through Alteration of the PI3K/Akt/NF-κB Signaling Pathway
Previous Article in Special Issue
Comparative Transcriptome Analysis between a Resistant and a Susceptible Wild Tomato Accession in Response to Phytophthora parasitica
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 3
10.3390/ijms20030610
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

Transcriptome Analysis of Watermelon Leaves Reveals Candidate Genes Responsive to Cucumber green mottle mosaic virus Infection

by
Yuyan Sun
1,2,
Min Fan
1,* and
Yanjun He
1
1
Institute of Vegetables, Zhejiang Academy of Agricultural Sciences, Hangzhou 310021, China
2
Key Laboratory of Biology and Genetic Improvement of Horticultural Crops, Ministry of Agriculture, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(3), 610; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030610
Submission received: 11 January 2019 / Revised: 24 January 2019 / Accepted: 28 January 2019 / Published: 31 January 2019
(This article belongs to the Special Issue Transcriptional and Post-transcriptional Gene Regulation in Plants)
Browse Figures
Versions Notes

Abstract

:
Cucumber green mottle mosaic virus (CGMMV) is a member of the genus Tobamovirus, which cause diseases in cucurbits, especially watermelon. In watermelon, symptoms develop on the whole plant, including leaves, stems, peduncles, and fruit. To better understand the molecular mechanisms of watermelon early responses to CGMMV infection, a comparative transcriptome analysis of 24 h CGMMV-infected and mock-inoculated watermelon leaves was performed. A total of 1641 differently expressed genes (DEGs) were identified, with 886 DEGs upregulated and 755 DEGs downregulated after CGMMV infection. A functional analysis indicated that the DEGs were involved in photosynthesis, plant–pathogen interactions, secondary metabolism, and plant hormone signal transduction. In addition, a few transcription factor families, including WRKY, MYB, HLH, bZIP and NAC, were responsive to the CGMMV-induced stress. To confirm the high-throughput sequencing results, 15 DEGs were validated by qRT-PCR analysis. The results provide insights into the identification of candidate genes or pathways involved in the responses of watermelon leaves to CGMMV infection.

Graphical Abstract

1. Introduction

Watermelon [Citrullus lanatus (Thunb.) Matsum. & Nakai] is an important cucurbit crop grown worldwide. Watermelon fruit provide glucose, fructose, malic acid, lycopene, vitamins, and amino acids, which are desired by consumers. In 2016, approximately 117 million tons of watermelon were produced, accounting for 10.9% of total global vegetable production (available online: http://www.fao.org/faostat). China was the leading country, producing approximately 79 million tons of watermelon, which constituted 67.7% of the worldwide production (available online: http://www.fao.org/faostat). Currently, watermelon production is threatened by the Cucumber green mottle mosaic virus (CGMMV).
CGMMV is a member of the genus Tobamovirus in the Virgaviridae family. It has a 6.4-kb single-stranded, positive-sense RNA genome, which produces a stiff, rod-shaped 300 × 18-nm particle [1]. CGMMV is readily spread by contact, and pollen and seed transmission [2]. Since it was first report in 1935 in England, where it was found to cause diseases in cucumber [3,4], CGMMV has achieved a global distribution, being present in more than 30 countries, and has become a major threat to cucurbit-based industries [5,6,7,8]. In CGMMV-infected watermelon, mottling and mosaic develop on the leaves, brown necrotic lesions develop on the stems and peduncles, and the flesh becomes spongy and rotting, with dirty red discolorations [9].
Several investigations of plant–virus interactions have focused on transcriptional or post-transcriptional responses to CGMMV infections [10,11,12,13,14,15]. For instance, miRNA sequencing analysis of cucumber and watermelon leaves infected by CGMMV revealed potential miRNAs and target genes involved in responses to CGMMV-induced stress [10,14]. Profiles of CGMMV-derived siRNAs in infected cucumber [11], bottle gourd [12] and watermelon [15] indicated that vsiRNAs have complicated functions during host–virus interactions. In addition, transcriptome analysis of watermelon fruit responding to CGMMV-induced stress revealed that some differentially expressed genes (DEGs) related to cell wall components and photosynthesis may be directly involved in the formation of diseased watermelon symptoms [13]. However, gene expression levels in leaves of CGMMV-infected watermelons have not been studied. Here, we analyzed the early responses of watermelon leaves to CGMMV infection using RNA sequencing (RNA-Seq) and investigated alterations in gene expression levels between healthy and CGMMV-infected plants. A subset of DEGs involved in photosynthesis, plant–pathogen interactions, secondary metabolism and plant hormone signal transduction were identified. The results provide insights into the identification of candidate genes or pathways associated with responses of watermelon leaves to CGMMV infection.

2. Results

2.1. Phenotypes and Confirmation of CGMMV in Watermelon Post Inoculation

At 24-h post-inoculation (24 hpi), there were no obvious disease symptoms exhibited in the watermelon leaves compared with the ck plants (Figure 1a,b). The accumulation of CGMMV in both ck and 24 hpi plants was confirmed by RT-PCR, the 24 hpi samples clearly amplified the specific 486-bp CGMMV CP fragment (Figure 1c). The results showed that even though no obvious symtoms developed in 24 hpi watermelon leaves, the virus has accumulated in the plants at the early infectious time.

2.2. Overview of the RNA-Seq Results

Here, triplicates of treatments for ck and 24 hpi samples were conducted, and six leaf samples were prepared for RNA extraction and Illumina sequencing. A total of 282,952,260 raw reads were generated, encompassing 42.44 Gb data, which were sufficient for the gene expression analysis. After removing the low-quality reads and adapter sequences, 276,176,444 (97.60%) clean reads were obtained. Using TopHat2 software, 266,058,279 (96.33%) clean reads were mapped to the watermelon genome, and 200,211,479 (72.49%) unique reads were mapped the reference genome. The Q20 percentage was over 99.60%, and the Q30 percentage was over 96.06% (Table 1).

2.3. Analysis of DEGs in Response to CGMMV-Induced Stress

Genes with a false discovery rate (FDR) < 0.05 and an estimated absolute log2 fold change (log2FC) ≥ 1 in sequence counts between libraries were considered significantly differentially expressed. Finally, 1641 DEGs were identified between the 24 hpi and ck leaf libraries (Figure 2; Table S1). Of these genes, 886 were up-regulated and 755 were down-regulated after CGMMV infection. To provide an expression overview of the DEGs, hierarchical clustering was generated and is shown in Figure 3.

2.4. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway Enrichment Analyses of DEGs

On the basis of the GO classifications, DEGs were separated into three categories: molecular function, cellular component, and biological process (Table S2). In the category of molecular function, pigment binding, chlorophyll binding, transcription factor (TF) activity, and sequence-specific DNA binding were significantly enriched (Figure 4; Table S2). In the category of cellular component, chloroplast thylakoid membrane, thylakoid, chloroplast stroma, chloroplast envelope, chloroplast, chloroplast thylakoid lumen, plastoglobule, apoplast, chloroplast thylakoid, thylakoid lumen, and cell wall were significantly enriched (Figure 4; Table S2). In the category of biological process, photosynthesis, light harvesting in photosystem (PS) I, response to abscisic acid (ABA), response to cold, response to heat, response to water deprivation, and photosynthesis were significantly enriched (Figure 4; Table S2).
In addition, 122 KEGG pathways for the DEGs were identified in watermelon leaves under CGMMV-induced stress (Table S3). Among them, the plant hormone signal transduction pathway (68 DEGs) was the largest group, followed by the plant–pathogen interaction pathway (56 DEGs) and the carbon metabolic pathway (39 DEGs) (Figure 5; Table S3). The top 20 enrichment pathways (p-value < 0.05) were photosynthesis p-antenna proteins (ko00196), glyoxylate and dicarboxylate metabolism (ko00630), alpha-linolenic acid metabolism (ko00592), fatty acid degradation (ko00071), glycerolipid metabolism (ko00561), photosynthesis (ko00195), galactose metabolism (ko00052), fatty acid metabolism (ko01212), alanine, aspartate and glutamate metabolism (ko00250), valine, leucine and isoleucine degradation (ko00280), cutin, suberine, and wax biosynthesis (ko00073), phenylpropanoid biosynthesis (ko00940), carotenoid biosynthesis (ko00906), isoquinoline alkaloid biosynthesis (ko00950), tyrosine metabolism (ko00350), ascorbate and aldarate metabolism (ko00053), beta-alanine metabolism (ko00410), circadian rhythm–plant (ko04712), biosynthesis of unsaturated fatty acids (ko01040), and plant hormone signal transduction (ko04075) (Figure 5; Table S3).

2.5. DEGs Involved in Photosynthesis

Photosynthesis, as the most fundamental and complex physiological process in plants, is disturbed by viral stress. Virus-infected plants are commonly characterized by decreased photosynthetic rates, which suggests that photosynthesis is a major activity suppressed by host antiviral defense responses [16,17]. In this study, 22 DEGs were associated with the pathway of photosynthesis (ko00195; ko00196) (Table 2). Among these, seven DEGs (Cla001715, Cla001790, Cla004698, Cla004699, Cla004704, Cla004703, and Cla007741), encoding the PsbQ, PsbO, PsbP, Psb27, and PsbW proteins, were involved in PSII. In addition, two DEGs (Cla007940 and Cla008429), encoding PSI reaction center subunits, were involved in PSI. Moreover, 10 DEGs (Cla004746, Cla011145, Cla011748, Cla012368, Cla013826, Cla018117, Cla019105, Cla019595, Cla022573, and Cla022963), encoding the chlorophyll a- and chlorophyll b-binding proteins, were involved in the light harvesting complex and serve as conduits of the excitation energy to the reaction center of PSII. All of these photosynthesis-related genes were down-regulated after CGMMV infection, suggesting a decreased photosynthetic efficiency upon CGMMV infection in the early stages.

2.6. DEGs Involved in Plant–Pathogen Interactions

Resistance to invading microorganisms is often governed by specific recognition between plant and pathogen proteins. In the present study, 55 DEGs (20 down- and 35 upregulated) were identified as participating in the pathway of plant–pathogen interaction (ko04626) (Table S4). These DEGs mainly encoded a receptor-like protein kinase, serine/threonine-protein kinase, disease resistance protein, calmodulin (CaM), CaM-like protein, and WRKY and MYB TFs. For example, DEGs encoding two CaMs (Cla003404 and Cla021803), one disease resistance protein (Cla020705), two WRKY TFs (Cla010918 and Cla014818), three MYB TFs (Cla005982, Cla015165, and Cla018610), 16 receptor-like protein kinases (Cla001490, Cla005698, Cla008301, Cla008723, Cla010837, Cla012567, Cla012568, Cla013903, Cla015191, Cla015419, Cla019852, Cla020279, Cla020764, Cla021452, Cla022195, and Cla022527) were upregulated after CGMMV infection. However, a few DEGs encoding receptor-like protein kinase, leucine-rich repeat-containing protein, and histone H2A were downregulated after CGMMV infection. These DEGs might play crucial roles in the responses to the CGMMV-induced stress.

2.7. DEGs Involved in Secondary Metabolism

Plant defense mechanisms against pathogens result in the differential regulation of various processes of primary and secondary metabolism. In the present study, a subset of DEGs involved in the pathways of secondary metabolism, such as phenylpropanoid biosynthesis (ko00940), carotenoid biosynthesis (ko00906), galactose metabolism (ko00052), and ascorbate and aldarate metabolism (ko00053), were significantly enriched (Table S5).
Pathogen infection often activates the phenylpropanoid pathway, leading to the production of both phytoalexins and lignin/suberin precursors for cell wall strengthening [18]. In this study, 37 DEGs were associated with the pathway of phenylpropanoid biosynthesis. Among which, 26 DEGs, encoding proteins such as cinnamoyl-CoA reductase (Cla000214, Cla017205, Cla017206, and Cla018764), cinnamyl alcohol dehydrogenase (Cla016048), caffeoyl-CoA O-methyltransferase (Cla022971), 4-coumarate CoA ligase-like (Cla022956), peroxidase (Cla001825, Cla003187, Cla003188, Cla003189, Cla003191, Cla003194, Cla013410, Cla014249, Cla017829, and Cla018276), were up-regulated after CGMMV infection. In addition, 11 DEGs, encoding proteins such as phenylalanine ammonia-lyase (Cla018297 and Cla018301), GDSL esterase/lipase (Cla018383 and Cla019629), and 4-coumarate-CoA ligase (Cla006818), were downregulated under the CGMMV-induced stress.
Carotenoids or carotenoid catabolites can act as stress-signaling molecules [19]. In total, 13 DEGs were associated with the carotenoid biosynthetic pathway in this study, and seven, which encoded an F-box protein (Cla005171), GDSL esterase/lipase (Cla010854), cytochrome P450 (Cla015020, Cla016011, Cla020673, and Cla022574), and UDP glycosyltransferase (Cla016026), showed up-regulated trends after CGMMV-induced stress. In contrast, six DEGs, encoding a short-chain dehydrogenase/reductase (Cla001908), phytoene synthase (Cla005425), lycopene cyclase (Cla016840), cytochrome P450 (Cla005457), GDSL esterase/lipase (Cla007748), and beta-hydroxysteroid dehydrogenase (Cla018423), were down-regulated under CGMMV-induced stress.
Ascorbic acid accumulates as a defense response to Turnip mosaic virus in resistant Brassica rapa cultivars [20]. In this study, 17 DEGs were involved in the ascorbate and aldarate metabolic pathway, and six of them, encoding laccase (Cla006606 and Cla014164), aldehyde dehydrogenase (Cla015539 and Cla019804), inositol oxygenase 1 (Cla016814), and mono dehydroascorbate reductase (Cla017711), were up-regulated under the CGMMV-induced stress conditions. Furthermore, 11 DEGs, encoding laccase (Cla001819, Cla008184, Cla011558, Cla015786, and Cla022385), GDSL esterase/lipase (Cla012640 and Cla012639), ascorbate peroxidase (Cla013254), UDP-D-glucose dehydrogenase (Cla019072), inositol oxygenase 1 (Cla011244), and inducer of CBF expression 1 (Cla002645), were down-regulated under CGMMV-induced stress conditions.
In total, 20 DEGs associated with the galactose metabolic pathway were identified. Among them, ten DEGs, encoding beta-galactosidase (Cla012613, Cla020797, Cla011211, and Cla011212), stachyose synthase (Cla017113 and Cla015152), galactinol synthase (Cla009222), reductase 1 (Cla018029), zinc finger (Cla014795), and alkaline alpha galactosidase I (Cla023372), were down-regulated. However, the other ten DEGs, aldose 1-epimerase (Cla018918), galactinol synthase (Cla010080, Cla010955, and Cla017336), NADP-dependent sorbitol 6-phosphate dehydrogenase (Cla015879 and Cla015878), stachyose synthase (Cla003446), beta-fructofuranosidase insoluble isoenzyme 2 (Cla015574), phosphofructokinase (Cla004847), and zinc finger (Cla018432), were up-regulated after the CGMMV infection.

2.8. DEGs Involved in Plant Hormone Signal Transduction

During viral infections, the disruption of the plant’s normal developmental physiology has been associated with alterations in phytohormone accumulation and signaling [21]. In this study, a few DEGs were involved in the plant hormone signal transduction pathway (K04075), particularly the accumulation and signaling of auxin and ABA (Table S6). For example, the expression levels of six auxin-responsive proteins (Cla002975, Cla006021, Cla014808, Cla019806, Cla021571, and Cla022332) increased under the CGMMV-induced stress conditions, whereas the expression of another five auxin-responsive proteins (Cla005087, Cla012819, Cla015002, Cla017664, and Cla020817) decreased. In addition, an auxin transporter-like protein (Cla004339) showed a decreased expression pattern. In general, ABA is involved in the negative regulation of plant defenses against various biotrophic and necrotrophic pathogens. In this study, ABA insensitive (Cla017696) and ABA receptor (Cla015009) were down-regulated after CGMMV infection. MYC and MYB proteins function in an ABA-dependent manner [22]. In the present study, genes encoding MYB (Cla012305) and MYC-like protein (Cla006058) were downregulated, whereas a MYB-like protein gene (Cla020369) was upregulated under the CGMMV-induced stress conditions, functioning in an ABA-dependent manner.

2.9. TFs Involved in CGMMV Stress Response

The invasion of plant viruses induces the expression of a variety of genes, which are usually regulated by TFs [23]. In this study, we observed 85 TFs that were differentially expressed after CGMMV infection (Figure 6; Table S7). These TFs included MYB (18 genes), NAC (15 genes), zinc finger (14 genes), HLH (11 genes), BZIP (7 genes), WRKY (5 genes), MADS box (5 genes), WD-40 (4 genes), ERF (3 genes), GRAS (2 genes), and SBP-box (1 gene). The WRKY family members were mainly associated with the plant–pathogen interaction pathway. The GRAS and SBP-box members were mainly involved in the plant hormone signal transduction pathway. The MYB family members were mainly involved in the plant hormone signal transduction, plant–pathogen interaction, oxidative phosphorylation, and RNA transport pathways. The HLH and bZIP family members were mainly involved in the plant hormone signal transduction and protein processing in the endoplasmic reticulum pathways. The NAC family members were involved in RNA degradation and glycerolipid metabolic pathways. The MADS box family members were involved in the RNA transport and protein export pathways. The WD-40 genes were involved in the ubiquitin-mediated proteolysis and ribosome biogenesis in eukaryotes pathways. The zinc finger genes were involved in the galactose metabolism, circadian rhythm, RNA transport, and ribosome pathways.

2.10. qRT-PCR Validation of DEGs

In total, 15 DEGs were selected for qRT-PCR validation (Figure 7). Cla010955 (galactinol synthase 3) was gradually up-regulated during the CGMMV-infection process. Cla001822 (MYB), Cla016019 (bZIP), Cla016840 (lycopene cyclase), and Cla017696 (ABA insensitive) were continually up-regulated at 0- (control, ck), 1-, and 6-h, slightly down-regulated at 18-h, and up-regulated again at 24- and 48-h post-CGMMV infection. The expression of Cla004746 (Chlorophyll a-b binding protein) decreased at ck, 1-, and 6-h, sharply increased at 18-h, and then decreased again at 24- and 48-h post-CGMMV infection. Cla003187 (peroxidase) was continually up-regulated at ck, 1-, 6-, 18-, and 24-h and then, was down-regulated at 48-h post-CGMMV infection. The expression patterns of Cla013731 (NAC) and Cla017392 (sucrose synthase) were similar, remaining relatively stably expressed at ck, 1-, and 6-h. They were then slightly down-regulated at 18-h and up-regulated at 24- and 48-h post-CGMMV infection. Cla008706 (XTH9) remained at a relatively stable expression level at ck, 1-, and 6-h, was sharply up-regulated at 18-h and then down-regulated at 24- and 48-h post-CGMMV infection. Cla009996 (HSP) remained relatively stably expressed at ck, 1-, and 6-h post-CGMMV infection. It was then sharply upregulated at 18-h, downregulated at 24-h, and upregulated again at 48-h post-CGMMV infection. Cla009996 (disease resistance-responsive protein) was upregulated at 1-h, downregulated at 6- and 18-h, sharply upregulated at 24-h and downregulated at 48-h post-CGMMV infection. Cla011558 (laccase 22) remained at a stable expression level at ck, 1-, and 6-h, decreased at 18- and 24-h, and was then upregulated at 48-h post-CGMMV infection. The expression of Cla012824 (cellulose synthase) was stable at ck, 1- and 6-h, increased at 18-h, and then decreased to a relatively stable level at 24- and 48-h post-CGMMV infection. Cla016011 (cytochrome P450) was downregulated at 1-h, upregulated at 6-h, then downregulated at 18-h, and upregulated at 24- and 48-h post-CGMMV infection. The qRT-PCR data were generally consistent with the high-throughput transcriptome results, indicating that the transcriptome sequencing data were reliable.

3. Discussion

First described in 1935 as infecting cucumber, CGMMV was one of the first plant viruses to be studied [4]. In CGMMV-infected watermelon, mottling and mosaicism develop on the leaves, brown necrotic lesions develop on the stems, and sponginess, rotting, and discolorations develop on the flesh. In addition to causing marketable yield losses owing to poor fruit quality, CGMMV also causes gross yield losses of > 50% in watermelon [24]. High-throughput sequencing technologies, such as RNA-Seq, are widely used to investigate the transcriptional profiles of different plants under virus-induced stress conditions [25,26]. RNA-Seq has also been utilized to identify the transcriptomic changes that occur during CGMMV-induced stress in bottle gourd leaves and fruit [27] and in watermelon fruit [13]. However, at present, there are no reports regarding the transcriptome early responses of watermelon leaves to CGMMV infection. Here, we performed transcriptome sequencing to assess the differences in gene expressions in the leaves of watermelon experiencing 24 h CGMMV-induced stress. The transcriptome data presented in this study provide insights into the early responses of watermelon leaves to CGMMV-induced stress and help us mine candidate genes or pathways to further illuminate the molecular mechanisms of watermelon responses to CGMMV-induced stress.

3.1. DEGs in Watermelon Leaves and Fruit

Under CGMMV-induced stress conditions, 1641 DEGs were identified in watermelon leaves, with 886 DEGs (53.99%) being upregulated and 755 DEGs (46.01%) being downregulated (Figure 3 and Table S1). These DEGs were mainly involved in photosynthesis, plant–pathogen interactions, secondary metabolism (phenylpropanoid biosynthesis, carotenoid biosynthesis, ascorbate and aldarate metabolism, and galactose metabolism), and plant hormone signal transduction pathways (Table S3). In a previous study, transcriptome profiles of CGMMV-infected and mock-inoculated watermelon fruit revealed 1621 DEGs, with 1052 up-regulated and 569 downregulated [13]. The DEG numbers in watermelon fruit was similar to the numbers in watermelon leaves. Consistent with our results, most of the DEGs in watermelon fruit were involved in metabolic processes, signal transduction, and plant–pathogen interactions [13]. Thus, common pathways were induced in the watermelon fruit and leaves by CGMMV infections.

3.2. Changes in Photosynthesis After CGMMV Infection

Changes in photosynthetic pathways after pathogen infection have been reported earlier [16,17]. Switching on the defense mechanisms and respiratory processes is a cost-intensive process, which might occur at the expense of photosynthesis [16,17]. In this study, 22 DEGs involved in the photosynthetic pathway were completely downregulated after CGMMV infection (Table 2). These DEGs encoded PsbO, PsbP, PsbQ, Psb27, PsbW, PSI reaction center subunits, and chlorophyll a- and b-binding proteins. The PSII of oxygenic organisms is a multi-subunit pigment–protein complex that catalyzes the photo-oxidation of water with a concomitant reduction of the catalysis of plastoquinone to plastoquinol [28]. PsbO, PsbP, and PsbQ proteins may perform alternative functions within PSII, and all three of the extrinsic proteins were required for the accumulation of PSII reaction centers in higher plants [29]. The PsbO component may exhibit carbonic anhydrase [30] and/or GTPase activity [31], the PsbP protein may be involved in binding the manganese required for photoactivation [32] and, along with PsbQ, participate in granal stack formation [33]. PsbW, a 6.1-kDa low-molecular-weight protein, is located close to the minor antenna of the PSII complex and is important for the contact and stability between several PSII–light harvesting complex II supercomplexes [34]. Psb27 plays a fundamental role in enabling plants to adapt to changes in light intensity independently of the formation of PSII supercomplexes and state transitions [35]. Chlorophyll a- and b-binding proteins are involved in the light harvesting complex and serve as conduits of excitation energy to the reaction center of PSII [36]. These reduced photosynthetic gene activity levels might be associated with chlorosis or decreased green-tissue surface areas after the systemic spread of a virus.

3.3. Changes in Plant–Pathogen Interactions After CGMMV Infection

Plants have evolved innate immune systems that recognize the presence of pathogens and initiate effective defense responses, whereas pathogens have evolved effector proteins that can suppress host immune responses. These form the so-called plant–pathogen interactions [37]. In the present study, 55 DEGs encoding receptor-like protein kinase, CaM-like protein, CaM, serine/threonine-protein kinase, and WRKY and MYB TFs were involved in the plant–pathogen interaction pathway (Table S4).
In response to several biotic and abiotic stresses, CaM plays an important role in sensing and transducing changes in cellular Ca2+ concentrations [38]. Several studies implicate CaMs in plant defense [39,40]. Transgenic tobacco (Nicotiana benthamiana) plants over-expressing SCaM4 and SCaM5 showed spontaneous lesions, increased PR genes expression and enhanced resistance to pathogens. The expression of tobacco NtCaM13, which is closely related to SCaM4 and SCaM5, was elevated at the RNA and protein levels in TMV-infected leaves [40]. In this study, two genes encoding CaMs were upregulated after the CGMMV infection, indicating these two CaMs were induced by the presence of CGMMV.
Receptor-like kinases/proteins (RLKs/RLPs) are on the front lines of the battle between plants and pathogens, because they are present at the plasma membrane and perceive signature molecules from either the invading pathogen or damaged plant tissue. These RLKs and RLPs perceive PAMPs from pathogens and initiate PAMP-triggered immunity, the first layer of plant innate immunity [41]. NbLRK1, a lectin RLK found in N. benthamiana, was identified as an interactor of the protein INF1, an elicitin from Phytophthora infestans [42]. Another RLK, was discovered in Nicotiana glutinosa and interacts directly with the fungal elicitin capsicein from Phytophthora capsici [43]. In this study, 26 RLKs were identified, with 16 upregulated and 10 downregulated after CGMMV infection. These RLKs might play crucial roles in watermelon’s innate immunity.

3.4. The Main TF Families Responding to CGMMV-Induced Stress

Briefly, TFs are composed of a DNA-binding domain that interacts with the cis-regulatory elements of its target genes and a protein–protein interaction domain that facilitates oligomerization between TFs and other regulators [44]. CGMMV-induced stress triggers extensive alterations in plant transcriptomes [13,27]. In a previous report, TFs, including those of the WRKY, MYB, bHLH, DREB, and ERF families, were significantly affected in watermelon fruit [13]. In our study, 85 differentially expressed TFs were identified in watermelon leaves after CGMMV infection (Figure 6; Table S7). The TF families identified in this study included MYB (18), NAC (15), zinc finger (14), HLH (11), BZIP (7), WRKY (5), MADS box (5), WD-40 (4), ERF (3), GRAS (2), and SBP-box (1). These TF families play vital roles in the transcriptional regulation of plants.
The plant-specific WRKY TFs form one of the largest TF families, and they are a class of DNA-binding proteins primarily found in plants. The WRKY domain contains approximately 60 amino acids, comprising a highly conserved short peptide, WRKYGQK, at the N-terminus and a C2H2 or C2HC zinc-binding motif at the C-terminus [45]. WRKYs play important roles in plant–pathogen interactions. WRKY27, a WRKY protein from pepper, positively regulates the stress-resistance response to a Ralstonia solanacearum infection by modulating SA-, JA- and ET-mediated signaling pathways in tobacco [46]. In tomato, six WRKY Group III TFs were identified, and these TFs respond to TYLCV infection [47]. In this study, we observed that two WRKY TFs (Cla017213 and Cla018197) were down-regulated and three WRKY TFs (Cla014433, Cla014818, and Cla010918) were up-regulated in the leaf tissues of watermelon after CGMMV infection (Figure 6; Table S7). Thus, WRKY TFs are also involved in CGMMV-induced stress responses.
NAC TFs belong to plant-specific TF family, and they have been identified in many plant species. The NAC protein contains a highly conserved DNA-binding domain at the N-terminus and diverse transcription activation or repression domains at the C-terminus [48,49]. In plants, the NAC TFs act as nodes of a regulatory network that responds to biotic stresses [50,51]. NAC TFs also take part in plant–pathogen interactions. Two wheat NAC TFs (TaNAC4 and TaNAC8) act as transcriptional activators in defense responses against strip rust pathogen infection [52,53]. The Arabidopsis NAC083 protein interacts with the Rep protein of Mungbean yellow mosaic India virus and plays important roles in viral DNA replication [51]. In tomato, six NAC TFs involved in the response to TYLCV infection were identified in resistant and susceptible cultivars [54]. In our study, 15 NAC TFs were affected by CGMMV infection in watermelon leaves (Figure 6; Table S7). Among them, 13 and 2 NAC TFs were up- and downregulated, respectively. The upregulated NAC TFs in leaf tissues might be involved in the adaptation to CGMMV-induced stress.
The MYB TFs are characterized by a DNA-binding MYB domain. The MYB domain is composed of approximately 52 amino acid residues that adopt a helix–turn–helix conformation that intercalates into the major groove of DNA [55,56]. After the first plant MYB gene, COLORED1, which is involved in anthocyanin biosynthesis, was identified in maize [57], a large number of MYB proteins were identified in different plant species. For example, in Arabidopsis, BOTRYTIS-SUSCEPTIBLE1, an R2R3-MYB gene (AtMYB108), is required for restricting the spread of two necrotrophic pathogens, Botrytis cinerea and Alternaria brassicicola, and is involved in the tolerance to osmotic and oxidative stresses [58]. In barley, the MYB TF MYB6 functions as positive regulator of basal and MLA-mediated immunity responses to Blumeria graminis [59]. In this study, we observed that 3 MYB TFs were down-regulated and 15 MYB TFs were up-regulated in the leaf tissues of watermelon after CGMMV infection (Figure 6; Table S7). In addition, we compared these MYB TFs in watermelon leaves with those in watermelon fruit. A MYB gene, Cla017179, was commonly up-regulated in both leaves and fruit after CGMMV infection, indicating the important role of Cla017179 in responses to CGMMV infection.
The bZIP TFs are vital players in plant innate immunity owing to their ability to regulate genes associated with PAMP-triggered immunity, effector-triggered immunity, and hormonal signaling networks [60,61]. One bZIP member, G/HBF-1 from Glycine max, binds to G-box and H-box motifs that have strong links to pathogen elicitors. These motifs have also been reported in the promoter of the chs15 gene, which is involved in flavonoid production in response to pathogen infection [62]. In this study, we observed that three bZIP TFs were down-regulated and four bZIP TFs were up-regulated in the leaf tissues of watermelon after CGMMV infection (Figure 6; Table S7). These bZIPs in watermelon act as positive or negative regulators and mediate the responses to CGMMV infection.

4. Materials and Methods

4.1. Plant Materials and the CGMMV Treatment

The seeds of the watermelon advanced inbred line ‘JJZ-M’, which is susceptive to CGMMV, were planted in a greenhouse and maintained at 25 °C. Seedlings at the two-true leaf stage were inoculated with the CGMMV virus. Mock plants were inoculated with sodium phosphate buffer (pH 7.2). Leaves were harvested at 24-h post-inoculation (24 hpi). The presence of CGMMV was confirmed by RT-PCR using CGMMV coat protein specific primers (F: 5′-ATGGCTTACAATCCGATCACAC-3′; R: 5′-CTAAGCTTTCGAGGTGGTAGCC-3′). Leaves from three plants were mixed to form a biological replicate, and three biological replicates for each treatment were used for the high-throughput RNA-Seq and analysis.

4.2. mRNA Library Construction and Sequencing

Total RNA was extracted using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s procedure. RNA quantity and purity were analyzed using a Bioanalyzer 2100 and RNA 6000 Nano Lab Chip Kit (Agilent, Santa Clara, CA, USA) with RIN number >7.0. Approximately 10 µg total RNA, representing a specific adipose type, was subjected to Poly(A) mRNA isolation using poly-T oligo-attached magnetic beads (Invitrogen, Carlsbad, CA, USA). Following purification, the mRNA was fragmented into small pieces using divalent cations under an elevated temperature. Then, the cleaved RNA fragments were reverse-transcribed to create the final cDNA library in accordance with the protocol for the mRNA-Seq sample preparation kit (Illumina, San Diego, CA, USA). The average insert size for the paired-end libraries was 300 bp (± 50 bp). Then, we performed the paired-end sequencing on an Illumina HiSeq 4000 at LC Sciences (Houston, TX, USA) following the vendor’s recommended protocol. Prior to assembly, the low-quality reads ((1) reads containing sequencing adaptors; (2) reads containing sequencing primers; and (3) nucleotides with q quality scores lower than 20) were removed. The raw sequence data was submitted to the NCBI Short Read Archive.

4.3. RNA-Seq Read Mapping

We aligned reads to the watermelon genome [63] using hierarchical indexing for spliced alignment of transcripts (HISAT) [64], which initially removed a portion of the reads based on the quality information accompanying each read and then mapped the reads to the reference genome. HISAT allowed multiple alignments per read (up to 20 by default) and a maximum of two mismatches when mapping the reads to the watermelon genome.

4.4. Transcript Abundance Estimation and Differential Expression Analysis

The mapped reads of each sample were assembled using StringTie [65]. Then, the transcriptomes from all the samples were merged to reconstruct a comprehensive transcriptome using perl scripts. After the final transcriptome was generated, StringTie and Ballgown [66] were used to estimate the expression levels of the transcripts. StringTie was used to determine expression levels for mRNAs by calculating fragments per kilobase of exon per million reads [67].
A differential expression analysis between the treatments (three biological replicates per treatment) was performed using the DESeq R package. DESeq provides statistical methods for determining differential expression in digital gene expression data using a model based on the negative binomial distribution. The differentially expressed mRNAs and genes were determined with log2 (fold change) > 1 or log2 (fold change) ≤ 1 and with statistical significance (p-value < 0.05) using the Ballgown R package.

4.5. GO and KEGG Enrichment Analyses of DEGs

The GO enrichment analysis of DEGs was performed using the GOseq R package in which gene length bias was corrected. GO terms with corrected p values < 0.05 were considered significantly enriched by DEGs. The KEGG database aids in understanding high-level functions and utilities of biological systems, such as the cell, organism and ecosystem, using molecular-level information, particularly large-scale datasets generated by genome sequencing and other high-throughput experimental technologies. We used KOBAS software to test the statistical enrichment of DEGs in KEGG pathways.

4.6. Validation of RNA-Seq Gene Expression Using qRT-PCR

To validate the gene expression levels revealed by the transcriptome data, leaves of watermelon plants were collected at 0-, 1-, 6-, 12-, 24- and 48-h post-CGMMV inoculation and used for RNA extraction. Total RNA (~2 μg) was reverse transcribed into cDNA using the TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix (TransGen Biotech, Beijing, China). cDNA was then used as the template for qRT-PCR, which was completed in a 20-μL sample containing 10 μL 2× TransStart Top Green qPCR SuperMix (TransGen Biotech, Beijing, China), 2.0 μL diluted cDNA, 0.4 μL 50× Passive Reference Dye I, and 0.4 μL each primer (10 mM). qRT-PCR was conducted using StepOne Real-Time PCR System (ABI, Foster City, CA, USA) with the following program: 95 °C for 30 s, 40 cycles of 95 °C for 5 s, 55 °C for 15 s, and 72 °C for 10 s. Each experiment was run in triplicate. The relative expression levels were calculated based on the expression of the housekeeping Tubulin gene. The fold changes in gene expression were estimated in terms of threshold cycles using the 2−∆∆CT method [68]. The primers designed for use in the qRT-PCR analysis are listed in Table S8.

5. Conclusions

In this study, we comparatively analyzed the gene expression profiles of 24 hpi and ck watermelon leaves. The CGMMV infection affected the gene expression levels of 1641 DEGs, with 886 upregulated and 755 downregulated. A functional analysis of DEGs based on a GO annotation and KEGG pathway analysis showed that these DEGs were mainly involved in photosynthesis, plant–pathogen interactions, secondary metabolism, and plant hormone signal transduction. The CGMMV infection also affected the expression levels of several TF families, including WRKY, NAC, MYB, bZIP, and zinc finger. Our genome-wide transcriptome analysis provides a basis for the further investigation of the molecular mechanisms underlying CGMMV infection in watermelon.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/20/3/610/s1. Table S1. Genes differentially expressed between mock- and CGMMV-infected watermelon leaf libraries. Table S2. GO term enrichment for differentially expressed genes. Table S3. KEGG pathway enrichment for differentially expressed genes. Table S4. Differentially expressed genes involved in the plant–pathogen interaction pathway. Table S5. Differentially expressed genes involved in the secondary biosynthetic pathway. Table S6. Differentially expressed genes involved in the plant hormone and signal transduction pathway. Table S7. Transcription factors involved in responses to CGMMV infection in watermelon leaves. Table S8. Primers used for qRT-PCR validation.

Availability of Sequencing Data

Sequencing data from this article were deposited at the sequence read archive (SRA) of the National Center for Biotechnology Information (NCBI) under the accession number of PRJNA514378.

Author Contributions

Y.S. and M.F. conceived and designed the experiments; Y.S. carried out the bioinformatics analysis and wrote the manuscript; M.F. and Y.H. gave insightful suggestions.

Acknowledgments

This work was funded by the Key Laboratory of Biology and Genetic Improvement of Horticultural Crops, Ministry of Agriculture, P. R. China (IVF201801), the Zhejiang Provincial Natural Science Foundation of China (LQ18C150002), the National Natural Science Foundation of China (31572145), and the Zhejiang Major Science and Technology Project for New Vegetable Variety Breeding (2016C02051-4-2). We thank Lesley Benyon, from Liwen Bianji, Edanz Group China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Compliance with Ethical Standards

This article does not contain any studies with human participants or animals performed by any of the authors.

References

  1. Tian, T.; Posis, K.; Maroon-Lango, C.J.; Mavrodieva, V.; Haymes, S.; Pitman, T.L.; Falk, W. First report of cucumber green mottle mosaic virus on melon in the United States. Plant Dis. 2014, 98, 1163. [Google Scholar] [CrossRef]
  2. Liu, H.W.; Luo, L.X.; Li, J.Q.; Liu, P.F.; Chen, X.Y.; Hao, J.J. Pollen and seed transmission of Cucumber green mottle mosaic virus in cucumber. Plant Pathol. 2014, 63, 72–77. [Google Scholar] [CrossRef]
  3. Mink, G.I. Pollen- and seed-transmitted viruses and viroids. Ann. Rev. Phytopathol. 1993, 31, 375–402. [Google Scholar] [CrossRef] [PubMed]
  4. Ainsworth, G.C. Mosaic disease of the cucumber. Ann. Appl. Biol. 1935, 22, 55–67. [Google Scholar] [CrossRef]
  5. Pop, I.; Jilaveanu, A. Identification of cucumber green mottle virus in Romania. Analele Institutului de Cercetari Pentru Protectia Plantelor 1985, 18, 43–47. [Google Scholar]
  6. Varveri, C.; Vassilakos, N.; Bem, F. Characterization and detection of Cucumber green mottle mosaic virus in Greece. Phytoparasitica 2002, 30, 493–501. [Google Scholar] [CrossRef]
  7. Li, R.; Zheng, Y.; Fei, Z.; Ling, K.S. First complete genome sequence of an emerging cucumber green mottle mosaic virus isolate in North America. Genome Announc. 2015, 3, e00452-15. [Google Scholar] [CrossRef]
  8. Tesoriero, L.A.; Chambers, G.; Srivastava, M.; Smith, S.; Conde, B.; Tran-Nguyen, L.T.T. First report of cucumber green mottle mosaic virus in Australia. Australas. Plant Dis. Notes 2016, 11, 1. [Google Scholar] [CrossRef]
  9. Dombrovsky, A.; Tran-Nguyen, L.T.T.; Jones, R.A.C. Cucumber green mottle mosaic virus: Rapidly increasing global distribution, etiology, epidemiology and management. Ann. Rev. Phytopathol. 2017, 55, 231–256. [Google Scholar] [CrossRef]
  10. Liu, H.W.; Luo, L.X.; Liang, C.Q.; Jiang, N.; Liu, P.F.; Li, J.Q. High-throughput sequencing identifies novel and conserved cucumber (Cucumis sativus L.) microRNAs in response to Cucumber green mottle mosaic virus infection. PLoS ONE 2015, 10, e0129002. [Google Scholar] [CrossRef]
  11. Li, Y.; Deng, C.; Shang, Q.; Zhao, X.; Liu, X.; Zhou, Q. Characterization of siRNAs derived from cucumber green mottle mosaic virus in infected cucumber plants. Arch. Virol. 2016, 161, 455–458. [Google Scholar] [CrossRef] [PubMed]
  12. Li, J.; Zheng, H.; Zhang, C.; Han, K.; Wang, S.; Peng, J.; Lu, Y.; Zhao, J.; Xu, P.; Wu, X.; et al. Different virus-derived siRNAs profiles between leaves and fruits in Cucumber green mottle mosaic virus-infected Lagenaria siceraria plants. Front. Microbiol. 2016, 7, 1797. [Google Scholar] [CrossRef] [PubMed]
  13. Li, X.; An, M.; Xia, Z.; Bai, X.; Wu, Y. Transcriptome analysis of watermelon (Citrullus lanatus) fruits in response to Cucumber green mottle mosaic virus (CGMMV) infection. Sci. Rep. 2017, 7, 16747. [Google Scholar] [CrossRef] [PubMed]
  14. Sun, Y.; Niu, X.; Fan, M. Genome-wide identification of cucumber green mottle mosaic virus-responsive microRNAs in watermelon. Arch. Virol. 2017, 162, 2591–2602. [Google Scholar] [CrossRef] [PubMed]
  15. Sun, Y.; Niu, X.; Cui, D.; Fan, M. High-throughput sequencing reveals vsiRNAs derived from Cucumber green mottle mosaic virus-infected watermelon. Hortic. Plant J. 2017, 3, 17–22. [Google Scholar] [CrossRef]
  16. Rahoutei, J.; García-Luque, I.; Barón, M. Inhibition of photosynthesis by viral infection: Effect on PSII structure and function. Physiol. Plant 2000, 110, 286–292. [Google Scholar] [CrossRef]
  17. Scharte, J.; Schon, H.; Weis, E. Photosynthesis and carbohydrate metabolism in tobacco leaves during an incompatible interaction with Phytophthora nicotianae. Plant Cell Environ. 2005, 28, 1421–1435. [Google Scholar] [CrossRef]
  18. Parker, D.; Beckmann, M.; Zubair, H.; Enot, D.P.; Caracuel-Rios, Z.; Overy, D.P.; Snowdon, S.; Talbot, N.J.; Draper, J. Metabolomic analysis reveals a common pattern of metabolic re-programming during invasion of three host plant species by Magnaporthe grisea. Plant J. 2009, 59, 723–737. [Google Scholar] [CrossRef]
  19. Fanciullino, A.L.; Bidel, L.P.R.; Urban, L. Carotenoid responses to environmental stimuli: Integrating redox and carbon controls into a fruit model. Plant Cell Environ. 2014, 37, 273–289. [Google Scholar] [CrossRef]
  20. Fujiwara, A.; Togawa, S.; Hikawa, T.; Matsuura, H.; Masuta, C.; Inukai, T. Ascorbic acid accumulates as a defense response to Turnip mosaic virus in resistant Brassica rapa cultivars. J. Exp. Bot. 2016, 67, 4391–4402. [Google Scholar] [CrossRef]
  21. Collum, T.D.; Culver, J.N. The impact of phytohormones on virus infection and disease. Curr. Opin. Virol. 2016, 17, 25–31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Abe, H.; Urao, T.; Ito, T.; Seki, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell 2003, 15, 63–78. [Google Scholar] [CrossRef] [PubMed]
  23. Nakashima, K.; Ito, Y.; Yamaguchi-Shinozaki, K. Transcriptional regulatory networks in response to abiotic stresses in Arabidopsis and grasses. Plant Physiol. 2009, 149, 88–95. [Google Scholar] [CrossRef] [PubMed]
  24. Reingold, V.; Lachman, O.; Koren, A.; Dombrovsky, A. First report of Cucumber green mottle mosaic virus (CGMMV) symptoms in watermelon used for the discrimination of non-marketable fruits in Israeli commercial fields. Plant Pathol. 2013, 28, 11. [Google Scholar]
  25. Fan, H.; Zhang, Y.; Sun, H.; Liu, J.; Wang, Y.; Wang, X.; Li, D.; Yu, J.; Han, C. Transcriptome analysis of beta macrocarpa and identification of differentially expressed transcripts in response to Beet necrotic yellow vein virus infection. PLoS ONE 2015, 10, e0132277. [Google Scholar] [CrossRef] [PubMed]
  26. Martin, K.; Singh, J.; Hill, J.H.; Whitham, S.A.; Cannon, S.B. Dynamic transcriptome profiling of Bean Common Mosaic Virus (BCMV) infection in Common Bean (Phaseolus vulgaris L.). BMC Genom. 2016, 17, 613. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, C.; Zheng, H.; Wu, X.; Xu, H.; Han, K.; Peng, J.; Lu, Y.; Lin, L.; Xu, P.; Wu, X.; et al. Genome-wide identification of new reference genes for RT-qPCR normalization in CGMMV-infected Lagenaria siceraria. Peer J. 2018, 6, e5642. [Google Scholar] [CrossRef] [PubMed]
  28. Nelson, N.; Yocum, C.F. Structure and function of Photosystems I and II. Ann. Rev. Plant Biol. 2006, 57, 521–565. [Google Scholar] [CrossRef]
  29. Bricker, T.M.; Frankel, L.K. Auxiliary functions of the PsbO, PsbP and PsbQ proteins of higher plant photosystem II: A critical analysis. J. Photochem. Photobiol. B Biol. 2011, 104, 165–178. [Google Scholar] [CrossRef]
  30. Lu, Y.K.; Stemler, A.J. Extrinsic Photosystem II carbonic anhydrase in maize mesophyll chloroplasts. Plant Physiol. 2002, 128, 643–649. [Google Scholar] [CrossRef]
  31. Spetea, C.; Hundal, T.; Lundin, B.; Heddad, M.; Adamska, I.; Andersson, B. Multiple evidence for nucleotide metabolism in the chloroplast thylakoid lumen. Proc. Natl. Acad. Sci. USA 2004, 101, 1409–1414. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Bondarava, N.; Krieger-Liszkay, A. Manganese binding to the 23 kDa extrinsic protein of Photosystem II. Biochim. Biophys. Acta 2007, 1767, 583–588. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. De Las Rivas, J.; Heredia, P.; Roman, A. Oxygen-evolving extrinsic proteins (PsbO, P, Q, R): Bioinformatic and functional analysis. Biochim. Biophys. Acta 2007, 1767, 575–582. [Google Scholar] [CrossRef] [PubMed]
  34. García-Cerdán, J.G.; Kovács, L.; Tóth, T.; Kereïche, S.; Aseeva, E.; Boekema, E.J.; Mamedov, F.; Funk, C.; Schröder, W.P. The PsbW protein stabilizes the supramolecular organization of photosystem II in higher plants. Plant J. 2011, 65, 368–381. [Google Scholar] [CrossRef] [PubMed]
  35. Hou, X.; Fu, A.; Garcia, V.J.; Buchanan, B.B.; Luan, S. PSB27: A thylakoid protein enabling Arabidopsis to adapt to changing light intensity. Proc. Natl. Acad. Sci. USA 2015, 112, 1613–1618. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Pietrzykowska, M.; Suorsa, M.; Semchonok, D.A.; Tikkanen, M.; Boekema, E.J.; Aro, E.M.; Jansson, S. The light-harvesting chlorophyll a/b binding proteins Lhcb1 and Lhcb2 play complementary roles during state transitions in Arabidopsis. Plant Cell 2014, 26, 3646–3660. [Google Scholar] [CrossRef] [PubMed]
  37. Dodds, P.N.; Rathjen, J.P. Plant immunity: Towards an integrated view of plant-pathogen interactions. Nat. Rev. Genet. 2010, 11, 539–548. [Google Scholar] [CrossRef]
  38. Yuenyong, W.; Chinpongpanich, A.; Comai, L.; Chadchawan, S.; Buaboocha, T. Downstream components of the calmodulin signaling pathway in the rice salt stress response revealed by transcriptome profiling and target identification. BMC Plant Biol. 2018, 18, 335. [Google Scholar] [CrossRef]
  39. Heo, W.D.; Lee, S.H.; Kim, M.C.; Kim, J.C.; Chung, W.S.; Chun, H.J.; Lee, K.J.; Park, C.Y.; Park, H.C.; Choi, J.Y.; et al. Involvement of specific calmodulin isoforms in salicylic acid-independent activation of plant disease resistance responses. Proc. Natl. Acad. Sci. USA 1999, 96, 766–771. [Google Scholar] [CrossRef] [Green Version]
  40. Yamakawa, H.; Mitsuhara, I.; Ito, N.; Seo, S.; Kamada, H.; Ohashi, Y. Transcriptionally and post-transcriptionally regulated response of 13 calmodulin genes to tobacco mosaic virus-induced cell death and wounding in tobacco plant. Eur. J. Biochem. 2001, 268, 3916–3929. [Google Scholar] [CrossRef] [Green Version]
  41. Yang, X.; Deng, F.; Ramonell, K.M. Receptor-like kinases and receptor-like proteins: Keys to pathogen recognition and defense signaling in plant innate immunity. Frontiers Biol. 2012, 7, 155–166. [Google Scholar] [CrossRef]
  42. Kanzaki, H.; Saitoh, H.; Takahashi, Y.; Berberich, T.; Ito, A.; Kamoun, S.; Terauchi, R. NbLRK1, a lectin-like receptor kinase protein of Nicotiana benthamiana, interacts with Phytophthora infestans INF1 elicitin and mediates INF1-induced cell death. Planta 2008, 228, 977–987. [Google Scholar] [CrossRef] [PubMed]
  43. Kim, Y.T.; Oh, J.; Kim, K.H.; Uhm, J.Y.; Lee, B.M. Isolation and characterization of NgRLK1, a receptor-like kinase of Nicotiana glutinosa that interacts with the elicitin of Phytophthora capsici. Mol. Biol. Rep. 2010, 37, 717–727. [Google Scholar] [CrossRef] [PubMed]
  44. Shiu, S.H.; Shih, M.C.; Li, W.H. Transcription factor families have much higher expansion rates in plants than in animals. Plant Physiol. 2005, 139, 18–26. [Google Scholar] [CrossRef] [PubMed]
  45. Eulgem, T.; Somssich, I.E. Networks of WRKY transcription factors in defense signaling. Curr. Opin. Plant Biol. 2007, 10, 366–371. [Google Scholar] [CrossRef] [Green Version]
  46. Dang, F.; Wang, Y.; She, J.; Lei, Y.; Liu, Z.; Eulgem, T.; Lai, Y.; Lin, J.; Yu, L.; Lei, D.; et al. Overexpression of CaWRKY27, a subgroup IIe WRKY transcription factor of Capsicum annuum, positively regulates tobacco resistance to Ralstonia solanacearum infection. Physiol. Plant 2014, 150, 397–411. [Google Scholar] [CrossRef]
  47. Huang, Y.; Li, M.Y.; Wu, P.; Xu, Z.S.; Que, F.; Wang, F.; Xiong, A.S. Members of WRKY Group III transcription factors are important in TYLCV defense signaling pathway in tomato (Solanum lycopersicum). BMC Genom. 2016, 17, 788. [Google Scholar] [CrossRef] [PubMed]
  48. Olsen, A.N.; Ernst, H.A.; Lo Leggio, L.; Skriver, K. NAC transcription factors: Structurally distinct, functionally diverse. Trends Plant Sci. 2005, 10, 1360–1385. [Google Scholar] [CrossRef]
  49. Hu, R.; Qi, G.; Kong, Y.; Kong, D.; Gao, Q.; Zhou, G. Comprehensive analysis of NAC domain transcription factor gene family in Populus trichocarpa. BMC Plant Biol. 2010, 10, 145. [Google Scholar] [CrossRef]
  50. Wang, X.; Basnayake, B.M.; Zhang, H.; Li, G.; Li, W.; Virk, N.; Mengiste, T.; Song, F. The Arabidopsis ATAF1, a NAC transcription factor, is a negative regulator of defense responses against necrotrophic fungal and bacterial pathogens. Mol. Plant-Microbe Interact. 2009, 22, 1227–1238. [Google Scholar] [CrossRef]
  51. Suyal, G.; Rana, V.S.; Mukherjee, S.K.; Wajid, S.; Choudhury, N.R. Arabidopsis thaliana NAC083 protein interacts with Mungbean yellow mosaic India virus (MYMIV) Rep protein. Virus Genes 2014, 48, 486–493. [Google Scholar] [CrossRef] [PubMed]
  52. Xia, N.; Zhang, G.; Liu, X.Y.; Deng, L.; Cai, G.L.; Zhang, Y.; Wang, X.J.; Zhao, J.; Huang, L.L.; Kang, Z.S. Characterization of a novel wheat NAC transcription factor gene involved in defense response against stripe rust pathogen infection and abiotic stresses. Mol. Biol. Rep. 2010, 37, e3703–e3712. [Google Scholar] [CrossRef] [PubMed]
  53. Xia, N.; Zhang, G.; Sun, Y.F.; Zhu, L.; Xu, L.S.; Chen, X.M.; Liu, B.; Yu, Y.T.; Wang, X.J.; Huang, L.L.; et al. TaNAC8, a novel NAC transcription factor gene in wheat, responds to stripe rust pathogen infection and abiotic stresses. Physiol. Mol. Plant 2010, 74, e394–e402. [Google Scholar] [CrossRef]
  54. Huang, Y.; Li, T.; Xu, Z.; Wang, F.; Xiong, A. Six NAC transcription factors involved in response to TYLCV infection in resistant and susceptible tomato cultivars. Plant Physiol. Biochem. 2017, 120, e61–e74. [Google Scholar] [CrossRef] [PubMed]
  55. Al-Attala, M.N.; Wang, X.; Abou-Attia, M.A.; Duan, X.; Kang, Z. A novel TaMYB4 transcription factor involved in the defence response against Puccinia striiformis f. sp. tritici and abiotic stresses. Plant Mol. Biol. 2014, 84, 589–603. [Google Scholar] [CrossRef] [PubMed]
  56. Dubos, C.; Stracke, R.; Grotewold, E.; Weisshaar, B.; Martin, C.; Lepiniec, L. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010, 15, 573–581. [Google Scholar] [CrossRef] [PubMed]
  57. Paz-Ares, J.; Ghosal, D.; Wienand, U.; Peterson, P.A.; Saedler, H. The regulatory c1 locus of Zea mays encodes a protein with homology to myb proto-oncogene products and with structural similarities to transcriptional activators. EMBO. J. 1987, 6, 3553–3558. [Google Scholar] [CrossRef]
  58. Mengiste, T.; Chen, X.; Salmeron, J.; Dietrich, R. The BOTRYTIS SUSCEPTIBLE1 gene encodes an R2R3MYB transcription factor protein that is required for biotic and abiotic stress responses in Arabidopsis. Plant Cell 2003, 15, 2551–2565. [Google Scholar] [CrossRef]
  59. Chang, C.; Yu, D.; Jiao, J.; Jing, S.; Schulze-Lefert, P.; Shen, Q.H.; Barley, M.L. A immune receptors directly interfere with antagonistically acting transcription factors to initiate disease resistance signaling. Plant Cell 2013, 25, 1158–1173. [Google Scholar] [CrossRef]
  60. Kuhlmann, M.; Horvay, K.; Stathmann, A.; Heinekamp, T.; Fischer, U.; Bottner, S.; Droge-Laser, W. The alpha-helical D1 domain of the bZIP transcription factor BZI-1 interacts with the ankyrin-repeat protein ANK1, and is essential for BZI-1 function, both in auxin signaling and pathogen response. J. Biol. Chem. 2003, 278, 8786–8794. [Google Scholar] [CrossRef]
  61. Murmu, J.; Bush, M.J.; DeLong, C.; Li, S.; Xu, M.; Khan, M.; Malcolmson, C.; Fobert, P.R.; Zachgo, S.; Hepworth, S.R. Arabidopsis basic leucine-zipper transcription factors TGA9 and TGA10 interact with floral glutaredoxins ROXY1 and ROXY2 and are redundantly required for anther development. Plant Physiol. 2010, 154, 1492–1504. [Google Scholar] [CrossRef] [PubMed]
  62. Noman, A.; Fahad, S.; Aqeel, M.; Ali, U.; Amanullah; Anwar, S.; Baloch, S.K.; Zainab, M. miRNAs: Major modulators for crop growth and development under abiotic stresses. Biotechnol. Lett. 2017, 39, 685–700. [Google Scholar] [CrossRef] [PubMed]
  63. Guo, S.; Zhang, J.; Sun, H.; Salse, J.; Lucas, W.J.; Zhang, H.; Zheng, Y.; Mao, L.; Ren, Y.; Wang, Z.; et al. The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions. Nat. Genet. 2013, 45, 51–58. [Google Scholar] [CrossRef]
  64. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed]
  65. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-Seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef] [PubMed]
  66. Frazee, A.C.; Pertea, G.; Jaffe, A.E.; Langmead, B.; Salzberg, S.L.; Leek, J.T. Ballgown bridges the gap between transcriptome assembly and expression analysis. Nat. Biotechnol. 2015, 33, 243–246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Trapnell, C.; Williams, B.A.; Pertea, G.; Mortazavi, A.; Kwan, G.; van Baren, M.J.; Salzberg, S.L.; Wold, B.J.; Pachter, L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 2010, 28, 511–515. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Ijms 20 00610 g001 550
Figure 1. Phenotypes and confirmation of Cucumber green mottle mosaic virus (CGMMV) in watermelon leaves. (a) ck watermelon leaf; (b) 24-h post-inoculation (24 hpi) watermelon leaf; (c) RT-PCR confirmation of CGMMV accumulation in watermelon leaves.
Figure 1. Phenotypes and confirmation of Cucumber green mottle mosaic virus (CGMMV) in watermelon leaves. (a) ck watermelon leaf; (b) 24-h post-inoculation (24 hpi) watermelon leaf; (c) RT-PCR confirmation of CGMMV accumulation in watermelon leaves.
Ijms 20 00610 g001
Ijms 20 00610 g002 550
Figure 2. Volcano map representing the number of differentially expressed genes of watermelon leaf between the 24 hpi and ck. Red dots represent up-regulated genes, and green dots represent down-regulated genes.
Figure 2. Volcano map representing the number of differentially expressed genes of watermelon leaf between the 24 hpi and ck. Red dots represent up-regulated genes, and green dots represent down-regulated genes.
Ijms 20 00610 g002
Ijms 20 00610 g003 550
Figure 3. Hierarchical clustering of differential expression profiles for 1641 genes between the 24 hpi and ck watermelon leaf libraries (false discovery rate (FDR) ≥ 0.05 and absolute value of the log2 ratio > 1).
Figure 3. Hierarchical clustering of differential expression profiles for 1641 genes between the 24 hpi and ck watermelon leaf libraries (false discovery rate (FDR) ≥ 0.05 and absolute value of the log2 ratio > 1).
Ijms 20 00610 g003
Ijms 20 00610 g004 550
Figure 4. Top 20 gene ontology (GO) terms enriched for different expressed genes. The size of each circle represents the number of significantly differentially expressed genes enriched in the corresponding term. The enrichment factor was calculated using the number of enriched genes divided by the total number of background genes in the corresponding pathway. The p-value was calculated using the Benjamini–Hochberg correction. A term with p < 0.05 is considered significantly overrepresented.
Figure 4. Top 20 gene ontology (GO) terms enriched for different expressed genes. The size of each circle represents the number of significantly differentially expressed genes enriched in the corresponding term. The enrichment factor was calculated using the number of enriched genes divided by the total number of background genes in the corresponding pathway. The p-value was calculated using the Benjamini–Hochberg correction. A term with p < 0.05 is considered significantly overrepresented.
Ijms 20 00610 g004
Ijms 20 00610 g005 550
Figure 5. Statistics of the top 20 pathways enriched for differentially expressed genes. The size of each circle represents the number of significantly differentially expressed genes enriched in the corresponding pathway. The enrichment factor was calculated using the number of enriched genes divided by the total number of background genes in the corresponding pathway. The p-value was calculated using the Benjamini–Hochberg correction. A pathway with p < 0.05 is considered significantly overrepresented.
Figure 5. Statistics of the top 20 pathways enriched for differentially expressed genes. The size of each circle represents the number of significantly differentially expressed genes enriched in the corresponding pathway. The enrichment factor was calculated using the number of enriched genes divided by the total number of background genes in the corresponding pathway. The p-value was calculated using the Benjamini–Hochberg correction. A pathway with p < 0.05 is considered significantly overrepresented.
Ijms 20 00610 g005
Ijms 20 00610 g006 550
Figure 6. Transcription factors involved in responses to CGMMV infection in watermelon leaves.
Figure 6. Transcription factors involved in responses to CGMMV infection in watermelon leaves.
Ijms 20 00610 g006
Ijms 20 00610 g007 550
Figure 7. qRT-PCR analysis of the expression of CGMMV response-related genes. The y-axis represents the log2 values (CGMMV/ck). Analyses were completed in triplicate, and the error bars represent the standard errors.
Figure 7. qRT-PCR analysis of the expression of CGMMV response-related genes. The y-axis represents the log2 values (CGMMV/ck). Analyses were completed in triplicate, and the error bars represent the standard errors.
Ijms 20 00610 g007
Table
Table 1. Overview of the watermelon leaf RNA sequencing (RNA-Seq) data.
Table 1. Overview of the watermelon leaf RNA sequencing (RNA-Seq) data.
SampleRaw ReadsRaw BaseClean ReadsMapped ReadsUnique Mapped ReadsQ20%Q30%
ck_152,659,470 7.90 G52,037,998 (98.82%)50,088,200 (96.25%)38,860,793 (74.68%)99.7596.64
ck_249,566,772 7.44 G49,101,272 (99.06%)46,572,421 (94.85%)32,009,543 (65.19%)99.6796.41
ck_343,812,198 6.57 G43,327,768 (98.89%)41,915,065 (96.74%)32,382,435 (74.74%)99.7196.34
24h_141,190,724 6.18 G36,924,504 (89.64%)35,371,699 (95.79%)27,425,778 (74.28%)99.6096.56
24h_246,758,1327.01 G46,328,460 (99.08%)44,965,010 (97.06%)33,083,469 (71.41%)99.8397.30
24h_348,964,964 7.34 G48,456,442 (98.96%)47,145,884 (97.30%)36,449,461 (75.22%)99.7196.06
Table
Table 2. Differentially expressed genes in watermelon leaves involved in photosynthesis.
Table 2. Differentially expressed genes in watermelon leaves involved in photosynthesis.
Gene NameDescriptionck_1ck_2ck_324h_124h_224h_3log2(fc)Regulation
Cla000152Sex-linked protein 9 (Fragment)275.6343.6292.0165.5117.5124.1−1.16down
Cla001715PsbQ 23.441.527.017.713.210.5−1.15down
Cla001790Oxygen-evolving enhancer protein 1 of photosystem II 223.666.0231.725.520.8144.1−1.45down
Cla004698PsbP 56.598.356.126.925.522.6−1.49down
Cla004704PsbP 47.660.144.426.317.621.4−1.22down
Cla004703PsbP 47.660.144.426.317.621.4−1.22down
Cla007741Photosystem II protein Psb27 33.540.937.720.212.013.7−1.29down
Cla007940Photosystem I reaction center subunit XI 13.917.215.06.46.57.3−1.19down
Cla008429Photosystem I reaction center subunit N 230.9231.3221.6128.392.9105.8−1.06down
Cla008898Ferredoxin--NADP reductase 607.9870.4595.8381.3298.7291.4−1.09down
Cla019799Photosystem II reaction center W protein25.111.212.54.77.28.2−1.27down
Cla019798Photosystem II reaction center W protein25.111.212.54.77.28.2−1.27down
Cla004746Chlorophyll a-b binding protein 6A308.0547.5409.9201.0146.4115.9−1.45down
Cla011145Chlorophyll a-b binding protein92.2154.387.847.634.851.2−1.32down
Cla011748Chlorophyll a-b binding protein 13192.7320.0218.0152.7105.874.4−1.13down
Cla012368Chlorophyll a-b binding protein 8733.71294.7799.6545.9441.2373.9−1.06down
Cla013826Chlorophyll a-b binding protein208.0246.4260.190.285.795.0−1.40down
Cla018117Chlorophyll a-b binding protein 634.131.536.913.911.815.4−1.32down
Cla019105Chlorophyll a-b binding protein P4903.51298.71142.1647.1447.7408.2−1.15down
Cla019595Chlorophyll a-b binding protein 214.720.26.34.21.62.5−1.89down
Cla022573Chlorophyll a-b binding protein 421.141.426.210.88.910.0−1.58down
Cla022963Chlorophyll a-b binding protein 71106.71615.01306.7786.1563.1628.6−1.03down

Share and Cite

MDPI and ACS Style

Sun, Y.; Fan, M.; He, Y. Transcriptome Analysis of Watermelon Leaves Reveals Candidate Genes Responsive to Cucumber green mottle mosaic virus Infection. Int. J. Mol. Sci. 2019, 20, 610. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030610

AMA Style

Sun Y, Fan M, He Y. Transcriptome Analysis of Watermelon Leaves Reveals Candidate Genes Responsive to Cucumber green mottle mosaic virus Infection. International Journal of Molecular Sciences. 2019; 20(3):610. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030610

Chicago/Turabian Style

Sun, Yuyan, Min Fan, and Yanjun He. 2019. "Transcriptome Analysis of Watermelon Leaves Reveals Candidate Genes Responsive to Cucumber green mottle mosaic virus Infection" International Journal of Molecular Sciences 20, no. 3: 610. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030610

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