Next Article in Journal
An Endophytic Fungi-Based Biostimulant Modulates Volatile and Non-Volatile Secondary Metabolites and Yield of Greenhouse Basil (Ocimum basilicum L.) through Variable Mechanisms Dependent on Salinity Stress Level
Previous Article in Journal
Modern Development and Production of a New Live Attenuated Bacterial Vaccine, SCHU S4 ΔclpB, to Prevent Tularemia
Previous Article in Special Issue
Lack of Evidence for Transmission of Verticillium dahliae by the Olive Bark Beetle Phloeotribus scarabaeoides in Olive Trees
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Pathogens
Volume 10
Issue 7
10.3390/pathogens10070796
pathogens-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

The Transcriptomic Profile of Watermelon Is Affected by Zinc in the Presence of Fusarium oxysporum f. sp. niveum and Meloidogyne incognita

by
Kasmita Karki
1,
Tim Coolong
2,
Chandrasekar Kousik
3,
Aparna Petkar
1,
Brendon K. Myers
1,
Abolfazl Hajihassani
1,
Mihir Mandal
3,4,* and
Bhabesh Dutta
1,*
1
Department of Plant Pathology, University of Georgia, Tifton, GA 31793, USA
2
Department of Horticulture, University of Georgia, Athens, GA 30602, USA
3
Vegetable Laboratory, USDA, ARS, Charleston, SC 29414, USA
4
Department of Biology, Claflin University, Orangeburg, SC 29115, USA
*
Authors to whom correspondence should be addressed.
Pathogens 2021, 10(7), 796; https://doi.org/10.3390/pathogens10070796
Submission received: 31 December 2020 / Revised: 19 May 2021 / Accepted: 20 June 2021 / Published: 23 June 2021
(This article belongs to the Special Issue Soil-Borne Plant Pathogenic Fungi)
Browse Figures
Versions Notes

Abstract

:
Zinc (Zn) accumulation and deficiency affect plant response to pests and diseases differently in varying pathosystems. The concentrations of Zn in plants aid in priming defense signaling pathways and help in enhanced structural defenses against plant pathogens. Studies are lacking on how concentrations of Zn in watermelon plants influence defense against two important soil-borne pathogens: Fusarium oxysporum f. sp. niveum (FON) and southern root-knot nematode (RKN, Meloidogyne incognita). In this study a comparative transcriptomics evaluation of watermelon plants in response to high (1.2 ppm) and low (0.2 ppm) levels of Zn were determined. Differential transcript-level responses differed in watermelon plants when infected with FON or RKN or both under high- and low-Zn treatment regimes in a controlled hydroponics system. Higher numbers of differentially expressed genes (DEGs) were observed in high-Zn-treated than in low-Zn-treated non-inoculated plants, in plants inoculated with FON alone and in plants inoculated with RKN alone. However, in the co-inoculated system, low-Zn treatment had higher DEGs as compared to high-Zn treatment. In addition, most DEGs were significantly enriched in hormone signal transduction and MAPK signaling pathway, suggesting an induction of systemic resistance with high-Zn concentrations. Taken together, this study substantially expands transcriptome data resources and suggests a molecular potential framework for watermelon-Zn interaction in FON and RKN.

1. Introduction

The Fusarium wilt fungus (Fusarium oxysporum f. sp. niveum, FON) and the root-knot nematode (Meloidogyne incognita (Kofoid and White) Chitwood; RKN) are important soil-borne pathogens causing severe damage in watermelon production throughout the world [1,2,3,4]. FON can infect and induce symptoms in watermelon plants at any growth stage [5]. Symptoms of Fusarium wilt at the seedling stage include seedling dieback and a scorched appearance. In mature plants, tissue chlorosis and unilateral wilting can be observed after loss of turgor pressure. As the disease progresses, complete vine wilting is often observed. FON has four physiological races (0, 1, 2, and 3) based on variability in aggressiveness on differential cultivars [6]. RKN is an obligatory parasite that induces galls on the infected root system of a susceptible host. Galls disrupt the vascular system of plants making them grow poorly and can lead to plant death under heavy infestation [7,8]. Pre-plant fumigant application with methyl bromide (MeBr) was often used and had been effective in reducing both FON and RKN in the watermelon production system. Phasing out of MeBr in accordance with the Montreal Protocol (1993) resulted in limited options for managing both soil-borne pathogens.
Triploid seedless watermelon cultivars, that are currently cultivated widely, do not possess a high level of resistance against FON race 2 [9]. To further complicate the issue, FON race 3 has also been identified in several states of the U.S. [6]. Proline (prothioconazole; Bayer Crop Science) is the only effective fungicide labeled for use on watermelon to manage FON [10]. Increased selection pressure on FON pathogen populations might lead to increased resistance to the fungicide. Grafting susceptible watermelon scions to interspecific hybrid squash (Cucurbita maxima Duch. Ex Lam. × C. moschata Duch. Ex Poir) and bottle gourd (Lagenaria siceraria (Molina) Standl.) can be a viable option for FON management. Two drawbacks with bottle gourd and C. moschata × C. maxima hybrid rootstocks are their susceptibility to root-knot nematode and the higher expense of grafted plants [11]. Therefore, available disease control options should be combined with alternative strategies for the integrated management of FON.
Plant disease resistance is affected by the plants’ genetics and by the environment in which it is grown, including nutrient deficiencies and toxicity [12,13,14]. Induced resistance has been proposed as one of the explanations for the interaction between specific nutrients and degree of resistance to plant pathogens [15]. Induced resistance is a plant-based defense system that is elicited by specific environmental stimuli, which enables them to resist biotic challenges [16]. These defense responses include the oxidative burst, changes in cell composition and synthesis of antimicrobial compounds such as phytoalexins [17]. The involvement of mineral nutrients in inducing resistance has been studied in several plant-pathogen systems. For instance, several reports indicate that potassium salts can be used for the induction of host-resistance against powdery mildew on cucumber (Cucumis sativus) [18], pepper (Capsicum annuum) [19], tomato (Solanum lycopersicum) [20] and sugar beet (Beta vulgaris) [21]. Potassium application has also been shown to increase downy mildew in muskmelon (Cucumis melo), caused by Pseudoperonospora cubensis, and Microdochium patch caused by Microdochium nivale in creeping bentgrass (Agrotis stolonifera) [22,23]. Despite these reports, the impact of a particular nutrient cannot be generalized for all plant-pest/pathogen systems, either globally or individually [12,24]. Therefore, a more detailed understanding of the relationship between plant responses to nutritional status and biotic stresses is essential.
Zinc plays a crucial role in the plant response to pests and diseases. Plants naturally absorb high concentrations of metals such as Zn from the substrate as a self-defense mechanism against pathogens and herbivores [25,26,27]. Metal ions may activate defense reactions through a signaling pathway or by plant fortification. Several studies have demonstrated that Zn fertilization decreases disease severity of Phytophthora root rot (Phytophthora megasperma) and common leaf spot (Pseudopeziza medicaginis) in alfalfa (Medicago sativa L.), peach gummosis (Botryosphaeria dothidia) in peach (Prunus persica) and early blight (Alternaria solani) in potato (Solanum tuberosum) [28,29,30]. However, a protective concentration of Zn against certain pathogens can also induce increased susceptibility to another pathogen on the same host [31]. Despite these conflicting reports on Zn and disease resistance and susceptibility, information is lacking on how Zn would affect susceptibility and resistance of watermelon plants against FON and RKN. It is not understood if transcriptome level changes can occur especially in pathways that are related to induced resistance. In the present study, we utilized a hydroponic system to evaluate if varying levels of Zn with or without challenge inoculation with either FON or RKN or both can affect genes in pathways related to induced resistance.

2. Results

2.1. Overview of the RNA-Seq Results

In order to explore the transcript-level changes in watermelon amended with different levels of Zn and pathogen inoculation (FON or RKN or both), cDNA libraries (n = 45) were constructed. Three biological replicates were used for each treatment. A total of 800,231,485 raw reads were generated. After removing the low-quality reads and trimming the adapter sequences, we obtained 787,477,379 (98.41%) clean reads. Using HISAT2 software, 765,227,330 (97.30%) clean reads were mapped to the reference watermelon genome, and 724,923,705 (92.15%) unique reads were obtained. The Q20 percentage was over 98%, and the Q30 percentage was over 98.49% (Table 1).

2.2. Differentially Expressed Genes (DEGs)

A comparative gene expression analysis of high-Zn (HZ) and low-Zn (LZ) treated plants under inoculated (FON (F) or RKN (R) or both (F_R)) and non-inoculated treatments was conducted. The non-inoculated micronutrient control (Stnr) was used as a control nutrient solution. In non-inoculated plants, the number of DEGs in plants treated with high-Zn was higher than plants treated with low-Zn at both 7 and 11 days post-treatment (dpt) (Supplementary Figure S1 and Figure 1). At 7 dpt, under high-Zn treatment eight DEGs (5 up, 3 dn) were identified whereas in plants treated with low-Zn treatment six DEGs (3 up, 3dn) were obtained. At 11 dpt, there were 333 DEGs (113 up, 220 dn) in plants treated with high-Zn whereas 14 DEGs (13 up, 1dn) in plants treated with low-Zn were identified.
In plants inoculated with FON alone, a similar pattern was observed: the number of DEGs in plants treated with high-Zn (n = 524) was higher than Steiner (n = 100) and low-Zn (n = 104) treated plants. In high-Zn treated plants inoculated with FON, 447 genes were upregulated and 77 genes were downregulated. Steiner treated plants inoculated with FON had 77 upregulated genes and 23 downregulated genes. Low-Zn treated plants inoculated with RKN had 79 upregulated genes and 25 downregulated genes.
Results showed a higher number of DEGs in high-Zn (n = 2908) treated plants compared to Steiner (n = 0) and low-Zn (n = 2) also in plants inoculated solely with RKN. Among the 2908 DEGs in high-Zn treated plants inoculated with RKN, 1522 genes were upregulated, and 1386 genes were downregulated. The two genes in low-Zn treated plants inoculated with RKN were upregulated.
In plants inoculated with both FON and RKN, low-Zn (n = 489) treated plants had the highest number of DEGs, followed by Steiner (n = 19) and high-Zn (n = 1). The one DEG in high-Zn treated plants inoculated with FON and RKN was upregulated. There were 11 upregulated and 8 downregulated genes in Steiner treated plants inoculated with FON and RKN. In low-Zn treated plants inoculated with FON and RKN, 268 genes were upregulated and 221 genes were downregulated.

2.3. GO Enrichment Analyses

GO terms with corrected p-value less than 0.05 were considered significantly enriched by differentially expressed genes and are shown in Figure 2 (Supplementary Table S2). Significantly enriched GO terms were observed only at 11 dpt. In the biological process, movement of cell or subcellular component, microtubule-based process and microtubule-based movement were the significantly enriched terms common in high-Zn treated plants inoculated with FON alone, low-Zn treated plants inoculated with FON alone, high-Zn treated plants inoculated with RKN alone, and low-Zn treated plants inoculated with FON and RKN (Figure 2). In the cellular component, thylakoid part, thylakoid and photosystem II were significantly enriched terms in both high-Zn treated plants, with and without inoculation with RKN (alone). In high-Zn treated plants inoculated with FON alone, chromosome, chromosomal part and chromatin were the most enriched terms (Figure 3).
Moreover, in molecular function, exopeptidase activity was the only significantly enriched term in low-Zn treated plants (Figure 4). In high-Zn treated plants inoculated with FON alone, transcription factor activity sequence-specific DNA binding, nucleic acid binding transcription factor activity and protein dimerization activity were the most enriched terms. The most enriched term in Steiner treated plant inoculated with FON alone were transcription factor activity sequence-specific DNA binding, nucleic acid binding transcription factor activity and sequence-specific DNA binding. Similarly, tubulin binding, protein complex binding and motor activity were the most enriched terms in low-Zn treated plants inoculated with FON alone. High-Zn treated plants inoculated with RKN alone showed the most enrichment of cytoskeletal protein binding, macromolecular complex binding and protein complex binding. In low-Zn treated plants inoculated with FON and RKN, pyrophosphatase activity, nucleoside-triphosphatase activity and hydrolase activity acting on acid anhydrides in phosphorus-containing anhydrides were the most enriched terms.

2.4. KEGG Classification of DEGs

KEGG enrichment results showed that DEGs were classified into multiple signaling pathways (Table 2 and Supplementary Table S3). Significant KEGG (corrected p-value < 0.05) enrichment was observed in plants treated with high-Zn at 7 dpt. At 11 dpt plants in treatments: high-Zn inoculated with FON alone, low-Zn inoculated with FON alone, and high-Zn inoculated with RKN alone, showed significant enrichment.
Peroxisome and glycosaminoglycan degradation were significantly enriched with two and one DEGs, respectively, in plants treated with high-Zn at 7 dpt. Two pathways, plant hormone signal transduction and plant-pathogen interaction, were significantly enriched in high-Zn treated plants inoculated with FON alone. In plants treated with low-Zn that were inoculated with FON, alpha-linolenic acid metabolism, valine, leucine and isoleucine biosynthesis and linoleic acid metabolism were significantly enriched. Metabolic pathways, carbon fixation in photosynthetic organisms, photosynthesis, pentose phosphate pathway, starch and sucrose metabolism, ascorbate and aldarate metabolism, carbon metabolism, photosynthesis-antenna proteins and carotenoid biosynthesis were enriched in high-Zn treated plants inoculated with RKN alone.

2.5. Investigation of DEGs Associated with Plant Phytohormone Signaling Pathway in Plants Treated with High- and Low-Zn

Phytohormone signaling has been shown to have a crucial regulatory role in plants [32]. Therefore, we have studied the effect of different levels of Zn on changes in the phytohormone signaling pathway. As shown in Table 3, many genes involved in the synthesis and regulation of abscisic acid, auxin, brassinosteroid, cytokinin, ethylene, jasmonic acid and salicylic acid were affected. The data indicate a synergistic regulation of these hormonal related DEGS as shown in Figure 5 through protein–protein interactome map studies using STRING database https://string-db.org/ (accessed on 16 May 2021) with Arabidopsis thaliana as the reference organism.
Jasmonic acid (JA) is a major class of plant hormones involved in mediating plant response to biotic and abiotic stress. Four genes predicted to be jasmonate zim domain protein (JAZ: ClCG07G014490, ClCG07G005870, ClCG06G001800, ClCG03G014890) and one gene predicted to be transcription factor MYC2-like protein (MYC2: ClCG10G001590) were induced. Number as well as expression of these genes was higher in high-Zn treatments (high-Zn inoculated with FON and RKN alone) (Table 3).
Salicylic acid (SA) is a key phytohormone that induces plants to produce systemic acquired resistance (SAR) to defend against pathogens [33,34]. Nonexpressor of pathogenesis-related genes 1 (NPR1: ClCG04G004170) and pathogenesis-related protein1 (PR-1: ClCG02G007190) are key regulators in the SA-dependent pathway (Nobuta et al. 2007). In high-Zn treated plants inoculated with RKN alone, one NPR1 gene was significantly upregulated, while one PR-1 gene was significantly downregulated.
Ethylene plays a vital role in plant development and the initiation of defense mechanisms against pathogens [35]. Only one gene associated with the ethylene signaling pathway, ethylene-responsive transcription factor (ERF1/2: ClCG05G004370), was induced in high-Zn treated plants inoculated with RKN alone.
Genes related to abscisic acid, auxin, brassinosteroid and cytokinin signaling pathways were differentially expressed with the difference in Zn level. In the abscisic acid pathway, serine/threonine-protein kinase SRK2 (SnRK2: ClCG00G003890), abscisic acid receptor PYR/PYL family (PYR/PYL: ClCG05G017790, ClCG05G000910), and protein phosphatase 2C (PP2C: ClCG03G001230, ClCG03G016180, ClCG05G008820) were affected. Both the genes related to the abscisic acid pathway were downregulated in plants treated with high-Zn, while in high-Zn treated plants inoculated with RKN, 60% of these genes were upregulated.
In the auxin signaling pathway, SAUR family protein (SAUR: ClCG05G018470), auxin responsive GH3 gene family (GH3: ClCG07G012040, ClCG05G015030, ClCG05G010400, ClCG01G017190), auxin-responsive protein IAA (AUX/IAA: ClCG11G006270, ClCG11G003560, ClCG07G013340, ClCG07G000600, ClCG06G011420, ClCG06G008840, ClCG05G004670, ClCG01G018130), auxin transporter-like protein 3 (AUX1: ClCG02G011740, ClCG02G004210), and auxin response factor (ARF: ClCG09G009470, ClCG05G018490, ClCG07G011440) were affected. A maximum number of genes were differentially regulated in the high-Zn treatments inoculated with FON and RKN alone. In plants treated with high-Zn, one auxin response factor (ClCG09G009470) was upregulated while another one (ClCG05G004670) was downregulated. In high-Zn treated plants inoculated with FON, all six genes (ClCG02G011740, ClCG05G010400, ClCG05G015030, ClCG05G018490, ClCG06G008840, ClCG11G006270) were upregulated. In high-Zn treated plants inoculated with RKN, 12 out of 15 genes were upregulated. Both low-Zn treated plants inoculated with FON and Steiner treated plants inoculated with FON had one upregulated gene. In low-Zn treated plants inoculated with FON and RKN, two (ClCG05G010400, ClCG05G015030) out of four genes were upregulated. In the brassinosteroid pathway, cyclin D3 (CYCD3: ClCG08G016160, ClCG05G020230, ClCG02G006220, ClCG05G010220) and brassinosteroid signaling positive regulator family protein (BZR1/2: ClCG08G015300) were affected, and all of these genes were upregulated. Similarly, histidine phototransfer protein (AHP: ClCG09G000190), Arabidopsis histidine kinase 2/3/4 (cytokinin receptor) (CRE1: ClCG11G001010) and type-B Arabidopsis response regulator (B-ARR: ClCG05G002710) were affected in the cytokinin pathway.

2.6. Investigation of DEGs Associated with MAPK and Zinc Finger Proteins Signaling Pathway in Plants Treated with High- and Low-Zn

Previous studies have demonstrated that the MAPK signaling pathway widely exists in eukaryotic organisms, participates in plant growth, and responds to abiotic and biotic stress [36,37]. As shown in Table 4 and Figure 6, a total of 26 genes in MAPK signaling were significantly affected. Results showed that the effect in terms of number and induction level of genes was higher in high-Zn treatments (high-Zn, high-Zn inoculated with FON alone and high-Zn inoculated with RKN alone). Genes such as WRKY DNA-binding protein 33 (WRKY33: ClCG03G000130, ClCG10G022500), respiratory burst oxidase protein D (RbohD: ClCG10G013100), MYC2 (ClCG10G001590, ClCG07G007900) and 1-amicocyclopropane-1-carboxylate synthase 6 (ACS6: ClCG07G007900) were significantly upregulated in high-Zn treated plants inoculated with FON and RKN alone.
In addition, the current study also investigated the regulation of zinc finger proteins (ZFP) in plants treated with high- and low-Zn. A previous report has indicated that plant zinc finger transcription factors are positive regulators of plant immunity, and the association of zinc finger domains in NBS-LRR resistance proteins is well conserved in plants and regulates the plant defense mechanism against diverse plant pathogens [38,39]. As shown in Table S4, a significant change in the regulation of zinc finger proteins in response to high-zinc (80 DEGs) and low-zinc (5 DEGs) treatment (Table S4) was observed. Results showed that 40 genes were upregulated by >two-fold under high-Zn inoculated with RKN alone, 12 genes under high-Zn inoculated with FON alone, and three genes on high-Zn uninoculated treatments. A significant upregulation of ZFPs on high-Zn, as compared to low or normal zinc conditions, supports our hypothesis of priming immunity against RKN and FON by activating the SA mediated signaling with increased transcription of downstream pathogenesis related (PR) genes.

2.7. Confirmation of Infection in FON and RKN Inoculation Plants

FON was re-isolated from 100% of the plants that were inoculated with FON or co-inoculated with RKN but not from non-inoculated Steiner control or plants only inoculated with RKN. Isolates were tested for FON by conventional PCR as described earlier [40]. One hundred percent of the putatively isolated FON colonies from infected plants were confirmed as FON using the PCR assay. No difference in plants being colonized by FON for all treatments and no Fusarium wilt symptoms were observed. Also, no differences in root gall formation between treatments were observed, except for a decrease in total nematode counts in plants inoculated with both FON and RKN under high zinc treatment (Figure 7). Galls were present only on plants inoculated with RKN or RKN co-inoculated with FON but were absent on non-inoculated or plants only inoculated with FON.

3. Discussion

In the present study, a whole transcriptome analysis was performed using RNA-seq in inoculated and non-inoculated (FON, RKN or both) watermelon plants in response to high- and low-Zn levels. We found a higher number of DEGs in high-Zn than in low-Zn-treated plants that were not inoculated; also in plants that were inoculated with FON or RKN alone. However, the effect of high-Zn was not consistent in a co-inoculated system; low-Zn treated plants inoculated with FON and RKN had higher DEGs than high-Zn treated plants inoculated with FON and RKN. These findings indicate that high-Zn treatment can alter differential expression of genes in watermelon, but the level of expression is not consistent across individual (FON or RKN) or co-inoculated (FON and RKN) pathosystems. We identified that the expression level of many genes associated with the phytohormone signaling pathway and the MAPK signaling pathway was affected in high-Zn treated plants, particularly in plants inoculated with FON or RKN alone. As the goal of this study was to evaluate the systemic response of watermelon away from the point of inoculation/infection, leaf tissues instead of root tissues were used for assessments. Within the time frame of the greenhouse study, we did not observe any wilting in plants inoculated with FON, RKN or both. However, fungal isolation and confirmation from the hypocotyl region confirmed the roots were indeed colonized by FON.
Hormone signaling pathways play significant roles in regulating interactions between plants and microorganisms [41,42]. We observed an increased expression of JA pathway genes: JAZ and MYC2 in high-Zn treated plants inoculated with FON or RKN alone. The core signal transduction mechanism of JA signaling comprises JAZ and MYC. Specific JAZ/TFs are generated by JAZ and different TFs that explicitly regulate many downstream responses [43]. The JAZ-MYC module triggers the plant defense response against pathogen infection by increasing the concentration of defense compounds such as indole alkaloids, terpenoid phytoalexins and often through airborne signals (e.g., green leaf volatiles and volatile terpenes) [44]. Gallego et al. [45] also found that surplus Zn can potentiate plant defense responses, especially in the synthesis of JA and it’s signaling pathway, thus improving plant resistance in Arabidopsis thaliana against Alternaria brassicicola. In addition, we observed an increased expression of NPR1 in high-Zn treated plants inoculated with RKN. NPR1 plays an integral part in the efficacy of the plant defense response. It activates PR gene expression by recruiting TGA transcription factors [46]. Arabidopsis NPR1 mutants showed decreased PR gene expression and increased susceptibility to pathogens [47]. Also, the ERF1/2 gene involved in ethylene signaling was upregulated in high-Zn treated plants inoculated with RKN. Activation of ERF genes is known to enhance plant disease resistance [48]. This could mean that high-Zn concentration in watermelon plants modulated induced resistance against FON and RKN. It is further justified by the number of genes involved in synthesis and regulation of SA, JA, abscisic acid, auxin, brassinosteroid, and cytokinin that were significantly affected in high-Zn treated plants inoculated with FON or RKN alone.
Plant MAPKs participate in plant growth, development, and responses to endogenous and environmental cues. Studies have shown that MAPKs could be activated by external sensors for cellular reactions [49]. A study conducted by Bi and Zhou [50] demonstrated that several pathogen-secreted effectors inhibit the MAPK cascade, which confirms the involvement of MAPKs in plant-pathogenic interactions and their role in plant response to pathogen invasion. Our study demonstrated that many genes in the MAPK pathway were significantly affected in high-Zn treated plants inoculated with FON or RKN alone. The genes that were significantly upregulated include WRKY33, RbohD, MYC2, and ACS6. Genes in the WRKY family have been confirmed to perform important regulatory functions to modulate pathogen-triggered cellular responses in a variety of plants, and most WRKY factors participate in the salicylic acid signaling pathway. RbohD is responsible for ROS production, which is involved in regulation of immune function against various pathogens [51,52]. As mentioned earlier, MYC2 a triggers defense response via the JA pathway. Similarly, ACS regulates synthesis of ethylene, which plays a positive role in host resistance against fungal and bacterial pathogens [53,54].
In summary, we observed that a high-Zn level affects plant–pathogen interaction in watermelon by regulating many crucial plant defense pathways (Figure 8, Figure S1–S6). Multiple genes in several hormone signaling pathways, including SA, JA, ethylene, abscisic acid, auxin, brassinosteroid and cytokinin were significantly affected. However, this effect was limited to FON-only and RKN-only inoculated systems. In the FON and RKN co-inoculated system, the induction of resistance in gall formation was observed; however, changes in DEGs were not observed. It is well known that FON and RKN together can be more severe to watermelon plants and it is possible that such severe interactions might counteract induced host-resistance with altered response in presence of high zinc levels. Pathogens are also known to evolve strategies to counteract Zn-related plant defense [55]. However, plant infection was confirmed at the end of this experiment by the recovery of the fungus and morphological and molecular identification for FON [56] and by presence of galls for RKN in all inoculated plants. Significant differences in terms of disease severity were not observed among inoculated treatments except for the reduced nematode counts when FON and RKN are inoculated simultaneously in high-Zn treated plants compared with Steiner treated plants with both pathogens inoculated. Future research is required to understand the transcriptional and translational basis of this observation. Further study is also required to optimize the level of Zn that can successfully reduce disease severity at the field level based on the effect of Zn ions on natural rhizospheral microbial diversity in natural field conditions. To our knowledge, this is the first study to analyze the effect of Zn level on watermelon in relation to FON or RKN or both infection under controlled hydroponics conditions.

4. Materials and Methods

4.1. Experimental Set-Up

Seeds of watermelon cultivar “Sugar Baby”, which is susceptible to both FON and RKN, were sown into sheets of 2.54 cm2/cell Rockwool cubes (Grodan Inc., Hedehusene, Denmark) and covered with a thin layer of vermiculite. The sheet was kept moist with water as needed. After germinating, seedlings were watered as needed and fertilized once with a “Steiner solution”, which was modified from the Steiner universal nutrient solution. The composition of the Steiner solution was as follows (mg.L−1): N-256 mg.L−1 [NH4H2PO4, KNO3, Ca(NO3)2], P-48 mg.L−1 (NH4H2PO4), K-304 (KNO3), Ca-180 mg.L−1 [Ca(NO3)2], Mg-48 mg.L−1 (MgSO4), B-1 mg.L−1 (H3BO3), Cu-0.2 mg.L−1 (CuSO4), Mo-0.1 mg.L−1 (Na2MoO4.2H2O), Fe-3 mg.L−1 (Fe Chelate; Sequestrene 330), Mn-1 mg.L−1 (MnCl2) and Zn-0.4 mg.L−1 (ZnSO4•7H2O). Seedlings were maintained for three weeks on Rockwool and then were transferred to plastic containers (43.05 cm W × 30.48 cm D × 19.81 cm H) containing 8L solution of respective micronutrient treatments. The micronutrient treatments were prepared by modifying the Zn composition in the Steiner solution mentioned above. The micronutrient treatments included: high-Zn (HZ: 3X concentration of Zn in Steiner solution, (1.2 ppm)), low-Zn (LZ: 0.5X concentration of Zn in Steiner solution, (0.2 ppm)), and Steiner (Stnr: X concentration of Zn (0.4 ppm)).
Three seedlings along with Rockwool were placed equidistant in a nine-well Styrofoam tray with holes at the bottom so that roots could suspend down and be in contact with the nutrient solution. The Styrofoam tray was then placed over the plastic container filled with micronutrient treatments, as indicated above. Solutions were aerated through a 15.2 cm aquarium air-stone fitted with plastic tubes (0.3-mm diameter) and an air pump (Pentair aquatic eco-systems Inc., Apopka, FL, USA). Plastic containers were maintained to the original volume (8L) by adding dH2O every two days. At 8 days post-treatment (dpt), the micronutrient treatments were either inoculated with FON or RKN or co-inoculated with both FON and RKN. For FON inoculation, a 1 mL suspension containing 5 × 105 microconidia/mL (race 2) was applied at the base of the watermelon seedling. For nematode inoculation, a 1ml suspension containing 6000 active RKN J2s (race 3) was added at the base of the seedling. For treatments that were inoculated with both FON and RKN, a similar inoculation approach was adopted and the inocula were applied simultaneously. Plants were maintained at 28 °C mean greenhouse temperature. Three replicates (plastic container with three seedlings each) per treatment were used in the experiment and treatments were arranged in a completely randomized design.

4.2. RNA Extraction, Quantification and Integrity Determination

As the goal of this study was to evaluate systemic response of watermelon away from the point of inoculation/infection, leaf tissues instead of root tissues were used for assessments. Leaf samples (n = 3 per replicate/treatment) were collected by cutting the 3rd or 4th leaf from the terminal with a pair of sterile scissors at 7 days post-treatment (dpt) and at 11 dpt (3 days post-inoculation (dpi)). Samples were immediately stored in liquid nitrogen and later transferred to a −80 °C freezer at the UGA Tifton Campus laboratory until needed for further analysis. Leaves were ground in liquid nitrogen and total RNA was extracted from 100 mg of ground leaf tissue using the manufacturer’s protocol (RNeasy Plant Mini Kit). Concentrations were determined using a NanoDropTM Lite (Thermo Scientific., Wilmington, DE, USA) and a Qubit® RNA Assay Kit and a Qubit® 2.0 Fluorometer (Life Technologies, Frederick, MD, USA). A total of 45 RNA samples were sent to Novogene Corporation Inc. Sacramento, CA, USA for library construction, sequencing and bioinformatic analysis. The RNA purity was checked using the NanoPhotometer® Spectrophotometer (IMPLEN, CA, USA). RNA integrity was confirmed (RIN > 7) using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, USA) with a minimum RNA integrated number of 8.

4.3. Library Preparation and Transcriptome Sequencing

A total amount of 1 µg RNA per sample was used to generate RNA-seq libraries using the NEBNext® UltraTM RNA library Prep Kit for Illimina® (NEB, Ipswich, MA, USA) following the recommendations of the manufacturer, and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA utilizing poly-T oligo-attached magnetic beads. Fragmentation was performed using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First-strand cDNA was synthesized using random hexamer primer and M- MuLV Reverse Transcriptase (RNAse H). Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends through exonuclease/polymerase activities. Following the adenylation of 3′ ends of DNA fragments, NEBNext Adaptor with hairpin loop structure was ligated to prepare for hybridization.
In order to select cDNA fragments of preferably 150–200 bp in length, the library fragments were purified using the AMPure XP system (Beckman Coulter, Beverly, MA, USA). Then, 3 µL USER Enzyme (NEB, Ipswich, MA, USA) was used with size-chosen, adapter ligated cDNA at 37 °C for 15 min followed by 5 min at 95 °C before PCR. PCR was then performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers, and Index (X) Primer. Finally, PCR products were purified (AMPure XP system) and the quality of the library was assessed on the Agilent Bioanalyzer 2100 system. The clustering of index-coded samples was performed on the cBOt Cluster Generation System using the P.E. Cluster Kit cBot-H.S. (Illumina) as instructed by the manufacturer. After cluster generation, the preparation of the library was sequenced on the Illumina platform and 150 bp paired-end reads were generated.

4.4. Sequence Read Cleanup and Mapping to Genome

First, raw data (raw reads) of the fastq format were processed. Clean reads were obtained by removing reads containing the adapter, reads containing poly-N and low-quality readings (<Q20) from raw data. At the same time, the percentages of reads with Q20, and Q30 were calculated. All downstream analyses were based on clean, high-quality ≥ Q20 data. Reference genome and gene model annotation files were downloaded directly from the genome website (http://cucurbitgenomics.org/pub/cucurbit/genome/watermelon/WCG/v2/) (Accessed on 15 April 2020). Index of the reference genome was developed using hisat2 2.1.0 and paired-end clean reads were aligned to the reference genome Charleston Gray [56,57] using HISAT2.

4.5. Quantification of Gene Expression Levels and Differential Expression of Analysis

FeatureCounts v1.5.0-p3 was used to count the read numbers mapped to each gene [55]. Fragments Per kilobase of transcript sequence per million or FPKM of each gene was then calculated based on the length of the gene and reads count mapped to this gene. Differential expression analysis was performed using the DESeq2 R package (1.14.1). The resulting p values were adjusted using the Benjamini and Hochberg method for controlling the false discovery rate (FDR). Corrected p-value of 0.05 and a log2 (foldchange) of 1 were set as the threshold for significantly differential expression.

4.6. GO and KEGG Enrichment Analysis of DEGs

Gene Ontology (GO) enrichment analysis of differentially expressed genes was performed by the clusterProfiler R package, in which gene length bias was corrected. GO terms with corrected p-value less than 0.05 were considered to be significantly enriched for a given set of genes. For pathway mapping, the Kyoto Encyclopedia of Genes and Genomics (KEGG) orthology database was adopted (http://www.genome.jp/kegg/) Accessed on 15 April 2020). We used clusterProfiler R package to test the statistical enrichment of differential expression genes in KEGG pathways.

4.7. Confirmation of Infection in FON and RKN Inoculation Plants

Visual symptoms of wilting were not observed in plants throughout the experiment. At the end of the experiment, one plant from each container (n = 2 containers/treatment), inoculated with FON alone, or co-inoculated with RKN and non-inoculated plants in the Steiner solution, were checked for the presence of FON. Stem pieces (0.5-cm-long) were cut from the base of the main stem of each plant using a pair of scissors sterilized by dipping in 70% ethanol. Stem pieces were surface-disinfested for 1.5 min with 0.6% sodium hypochlorite, rinsed in sterile water and placed on a semi-selective peptone pentachloronitrobenzene agar medium. Plates were incubated at 25 °C for 7 days. Fungal isolates were then microscopically identified based on morphological criteria [58] and further were confirmed as FON by PCR assay with FON specific primers Fon-1/Fon-2 [40]. Percentage of plants infested with FON as determined by morphological and PCR confirmatory assays were recorded. For plants inoculated with RKN, 10 cm of roots closest to the base of plants were sampled and evaluated, as galls were observed only in this portion. Rockwool plugs were removed from the root system, and roots were washed and rated visually for the presence of galls as described earlier [2,4]. Mean number of RKN adults for each treatment (t were calculated using Proc Glimmix with Beta distribution for percentage response variable in SAS 9.4 (SAS Institute). The effects of Zn application and pathogen exposure (FON or RKN or both) on RKN adult numbers in watermelon roots were also analyzed as described above.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/pathogens10070796/s1, Figure S1. Upregulated DEGs associated with diverse Biological Process BP under Hz_RKN condition. Figure S2. Upregulated DEGs associated with diverse Molecular Functions MF under Hz_RKN condition. Figure S3. Upregulated DEGs associated with diverse Cellular Component CC under Hz_RKN condition. Figure S4. Downregulated DEGs associated with diverse Biological Process BP under Hz_RKN condition. Figure S5. Downregulated DEGs associated with diverse Molecular Functions MF under Hz_RKN condition. Figure S6. Downregulated DEGs associated with diverse Cellular Component CC under Hz_RKN condition. Table S1. List of DEGs in different treatments. Table S2. Differentially expressed genes related to Gene Ontology enrichment analysis at 11 dpt. Table S3. KEGG analysis of differentially expressed genes at 11 dpt. Table S4. DEGs associated with zinc finger protein(s) and pathogenesis related (PR) proteins in response to high and low Zn treatment.

Author Contributions

Conceptualization, B.D. and K.K.; methodology, K.K., A.P.; validation, B.D., M.M., T.C., A.P. and C.K.; formal analysis, K.K; investigation, K.K.; resources, B.D., M.M., A.P. and T.C.; writing—original draft preparation, K.K.; writing—review and editing, K.K., B.D., M.M., T.C., C.K., A.P., B.K.M. and A.H.; visualization, K.K., M.M.; supervision, B.D., M.M.; project administration, B.D.; funding acquisition, B.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the USDA-NIFA-Methyl Bromide Transition Program; grant number 2018-51102-28743.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The transcriptomic sequencing data generated and analyzed in this study is available on request from [email protected].

Acknowledgments

We would like to thank Mike Foster and Michael Donahoo for their technical assistance in greenhouse.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Martyn, R.D. Fusarium wilt of watermelon: 120 years of research. Hortic. Rev. 2014, 42, 349–442. [Google Scholar]
  2. Davis, R.F. Effect of Meloidogyne incognita on watermelon yield. Nematropica 2007, 37, 287–294. [Google Scholar]
  3. Lynch, L.; Carpenter, J. The economic impacts of banning methyl bromide: Where do we need more research. In Proceedings of the 1999 Annual Research Conference on Methyl Bromide Alternatives and Emissions Reductions, San Diego, CA, USA, 1–4 November 1999; Available online: http://mbao.org/mbrpro99.html (accessed on 15 April 2021).
  4. Thies, J.A.; Ariss, J.J.; Hassell, R.L.; Olson, S.; Kousik, C.S.; Levi, A. Grafting for management of southern root-knot nematode, Meloidogyne incognita, in watermelon. Plant Dis. 2010, 94, 1195–1199. [Google Scholar] [CrossRef] [Green Version]
  5. Beattie, J.; Doolittle, S. Watermelons; US Government Printing Office: Washington, DC, USA, 1951. [Google Scholar]
  6. Zhou, X.; Everts, K.; Bruton, B. Race 3, a new and highly virulent race of Fusarium oxysporum f. sp. niveum causing Fusarium wilt in watermelon. Plant Dis. 2010, 94, 92–98. [Google Scholar] [CrossRef] [Green Version]
  7. Dale, P. Elimination of root-knot nematodes from roses by chemical bare-root dips. N. Z. J. Exp. Agric. 1973, 1, 121–122. [Google Scholar] [CrossRef]
  8. Newton, A.C.; Torrance, L.; Holden, N.; Toth, I.K.; Cooke, D.E.; Blok, V.; Gilroy, E.M. Climate change and defense against pathogens in plants. Adv. Appl. Microbiol. 2012, 81, 89–132. [Google Scholar]
  9. Wechter, W.P.; Kousik, C.; McMillan, M.; Levi, A.S. Identification of resistance to Fusarium oxysporum f. sp. niveum race 2 in Citrullus lanatus var. citroides plant introductions. HortScience 2012, 47, 334–338. [Google Scholar] [CrossRef] [Green Version]
  10. Miller, N.F.; Standish, J.R.; Quesada-Ocampo, L.M. Sensitivity of Fusarium oxysporum f. sp. niveum to prothioconazole and pydiflumetofen In vitro and efficacy for Fusarium wilt management in watermelon. Plant Health Prog. 2020, 21, 13–18. [Google Scholar] [CrossRef]
  11. Keinath, A.P.; Hassell, R.L. Suppression of Fusarium wilt caused by Fusarium oxysporum f. sp. niveum race 2 on grafted triploid watermelon. Plant Dis. 2014, 98, 1326–1332. [Google Scholar] [CrossRef]
  12. Agrios, G.N. Plant Pathology, 5th ed.; Elsevier Academic Press: San Diego, CA, USA, 2005. [Google Scholar]
  13. Krauss, A. Balanced nutrition and biotic stress. In Proceedings of the IFA Agricultural Conference on Managing Plant Nutrition, Barcelona, Spain, 29 June–2 July 1999. [Google Scholar]
  14. Marschner, H. Mineral Nutrition of Higher Plants; Academic Press: London, UK, 1995. [Google Scholar]
  15. Dordas, C. Role of nutrients in controlling plant diseases in sustainable agriculture. A review. Agron. Sustain. Dev. 2008, 28, 33–46. [Google Scholar] [CrossRef] [Green Version]
  16. Van Loon, L.; Bakker, P.; Pieterse, C. Systemic resistance induced by rhizosphere bacteria. Annu. Rev. Phytopathol. 1998, 36, 453–483. [Google Scholar] [CrossRef] [Green Version]
  17. Walters, D.; Walsh, D.; Newton, A.; Lyon, G. Induced resistance for plant disease control: Maximizing the efficacy of resistance elicitors. Phytopathology 2005, 95, 1368–1373. [Google Scholar] [CrossRef] [Green Version]
  18. Reuveni, R.; Dor, G.; Raviv, M.; Reuveni, M.; Tuzun, S. Systemic resistance against Sphaerotheca fuliginea in cucumber plants exposed to phosphate in hydroponics system, and its control by foliar spray of mono-potassium phosphate. Crop Prot. 2000, 19, 355–361. [Google Scholar] [CrossRef]
  19. Reuveni, M.; Oppenheim, D.; Reuveni, R. Integrated control of powdery mildew on apple trees by foliar sprays of mono-potassium phosphate fertilizer and sterol inhibiting fungicides. Crop Prot. 1998, 17, 563–568. [Google Scholar] [CrossRef]
  20. Ehret, D.; Utkhede, R.; Frey, B.; Menzies, J.; Bogdanoff, C. Foliar applications of fertilizer salts inhibit powdery mildew on tomato. Can. J. Plant Pathol. 2002, 24, 437–444. [Google Scholar] [CrossRef]
  21. Mosa, A. Management of sugar beet powdery mildew by foliar spraying of potassium phosphate salts. Arab Univ. J. Agric. Sci. 2002, 10, 1043–1057. [Google Scholar]
  22. Bains, S.; Jhooty, J. Effect of mineral nutrition of muskmelon on downy mildew caused by Pseudoperonospora cubensis. Plant Soil 1976, 29, 105. [Google Scholar]
  23. Bier, P.V.; Persche, M.; Koch, P.; Soldat, D.J. A long term evaluation of differential potassium fertilization of a creeping bentgrass putting green. Plant Soil 2018, 431, 303–316. [Google Scholar] [CrossRef]
  24. Huber, D. The role of mineral nutrition in defense. In Plant Disease: An Advanced Treatise: How Plants Defend Themselves; Horsfall, J.G., Cowling, E.B., Eds.; Academic Press: New York, NY, USA, 1980; pp. 381–406. [Google Scholar]
  25. Rolke, Y.; Liu, S.; Quidde, T.; Williamson, B.; Schouten, A.; Weltring, K.M.; Siewers, V.; Tenberge, K.B.; Tudzynski, B.; Tudzynski, P. Functional analysis of H2O2-generating systems in Botrytis cinerea: The major Cu-Zn-superoxide dismutase (BCSOD1) contributes to virulence on French bean, whereas a glucose oxidase (BCGOD1) is dispensable. Mol. Plant Pathol. 2004, 5, 17–27. [Google Scholar] [CrossRef]
  26. Babitha, M.; Bhat, S.; Prakash, H.; Shetty, H. Differential induction of superoxide dismutase in downy mildew-resistant and-susceptible genotypes of pearl millet. Plant Pathol. 2002, 51, 480–486. [Google Scholar] [CrossRef] [Green Version]
  27. Stolpe, C.; Giehren, F.; Krämer, U.; Müller, C. Both heavy metal-amendment of soil and aphid-infestation increase Cd and Zn concentrations in phloem exudates of a metal-hyperaccumulating plant. Phytochemistry 2017, 139, 109–117. [Google Scholar] [CrossRef]
  28. Grewal, H.S.; Graham, R.D.; Rengel, Z. Genotypic variation in zinc efficiency and resistance to crown rot disease (Fusarium graminearum Schw. Group 1) in wheat. Plant Soil. 1996, 186, 219–226. [Google Scholar] [CrossRef]
  29. Li, Z.; Fan, Y.; Gao, L.; Cao, X.; Ye, J.; Li, G. The dual roles of zinc sulfate in mitigating peach gummosis. Plant Dis. 2016, 100, 345–351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Machado, P.P.; Steiner, F.; Zuffo, A.M.; Machado, R.A. Could the supply of boron and zinc improve resistance of potato to early blight? Potato Res. 2018, 61, 169–182. [Google Scholar] [CrossRef]
  31. Helfenstein, J.; Pawlowski, M.L.; Hill, C.B.; Stewart, J.; Lagos-Kutz, D.; Bowen, C.R.; Frossard, E.; Hartman, G.L. Zinc deficiency alters soybean susceptibility to pathogens and pests. J. Plant Nutr. Soil Sci. 2015, 178, 896–903. [Google Scholar] [CrossRef]
  32. Flors, V.; Ton, J.; Van Doorn, R.; Jakab, G.; García-Agustín, P.; Mauch-Mani, B. Interplay between JA, SA and ABA signalling during basal and induced resistance against Pseudomonas syringae and Alternaria brassicicola. Plant J. 2008, 54, 81–92. [Google Scholar] [CrossRef]
  33. Kumar, D. Salicylic acid signaling in disease resistance. Plant Sci. 2014, 228, 127–134. [Google Scholar] [CrossRef]
  34. Maruri-López, I.; Aviles-Baltazar, N.Y.; Buchala, A.; Serrano, M. Intra and Extracellular Journey of the Phytohormone Salicylic Acid. Front. Plant Sci. 2019, 10, 423. [Google Scholar] [CrossRef]
  35. Blanco-Ulate, B.; Vincenti, E.; Powell, A.L.; Cantu, D. Tomato transcriptome and mutant analyses suggest a role for plant stress hormones in the interaction between fruit and Botrytis cinerea. Front. Plant Sci. 2013, 4, 142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Colcombet, J.; Hirt, H. Arabidopsis MAPKs: A complex signalling network involved in multiple biological processes. Biochem. J. 2008, 413, 217–226. [Google Scholar] [CrossRef]
  37. Rodriguez, M.C.; Petersen, M.; Mundy, J. Mitogen-activated protein kinase signaling in plants. Annu. Rev. Plant Biol. 2010, 61, 621–649. [Google Scholar] [PubMed]
  38. Cui, B.; Xu, S.; Li, Y.; Umbreen, S.; Frederickson, D.; Yuan, B.; Jiang, J.; Liu, F.; Pan, Q.; Loake, G.J. The Arabidopsis zinc finger proteins SRG2 and SRG3 are positive regulators of plant immunity and are differentially regulated by nitric oxide. New Phytol. 2021, 230, 259–274. [Google Scholar] [CrossRef] [PubMed]
  39. Gupta, S.K.; Rai, A.K.; Kanwar, S.S.; Sharma, T.R. Comparative analysis of zinc finger proteins involved in plant disease resistance. PLoS ONE 2012, 7, e42578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Lin, Y.H.; Chen, K.S.; Chang, J.Y.; Wan, Y.L.; Hsu, C.C.; Huang, J.W.; Chang, P.F.L. Development of the molecular methods for rapid detection and differentiation of Fusarium oxysporum and F. oxysporum f. sp. niveum in Taiwan. New Biotechnol. 2010, 27, 409–418. [Google Scholar] [CrossRef]
  41. Xie, S.; Jiang, H.; Ding, T.; Xu, Q.; Chai, W.; Cheng, B. Bacillus amyloliquefaciens FZB42 represses plant miR846 to induce systemic resistance via a jasmonic acid-dependent signalling pathway. Mol. Plant Pathol. 2018, 19, 1612–1623. [Google Scholar] [CrossRef] [Green Version]
  42. Vrabka, J.; Niehaus, E.-M.; Münsterkötter, M.; Proctor, R.H.; Brown, D.W.; Novák, O.; Pěnčik, A.; Tarkowská, D.; Hromadová, K.; Hradilová, M. Production and role of hormones during interaction of Fusarium species with maize (Zea mays L.) seedlings. Front. Plant Sci. 2019, 9, 1936. [Google Scholar] [CrossRef]
  43. Chini, A.; Gimenez-Ibanez, S.; Goossens, A.; Solano, R. Redundancy and specificity in jasmonate signalling. Curr. Opin. Plant Biol. 2016, 33, 147–156. [Google Scholar] [CrossRef]
  44. Havko, N.E.; Major, I.T.; Jewell, J.B.; Attaran, E.; Howe, G.A. Control of carbon assimilation and partitioning by jasmonate: An accounting of growth–defense tradeoffs. Plants 2016, 5, 7. [Google Scholar] [CrossRef] [Green Version]
  45. Gallego, B.; Martos, S.; Cabot, C.; Barceló, J.; Poschenrieder, C. Zinc hyperaccumulation substitutes for defense failures beyond salicylate and jasmonate signaling pathways of Alternaria brassicicola attack in Noccaea caerulescens. Physiol. Plant. 2017, 159, 401–415. [Google Scholar] [CrossRef]
  46. Després, C.; DeLong, C.; Glaze, S.; Liu, E.; Fobert, P.R. The Arabidopsis NPR1/NIM1 protein enhances the DNA binding activity of a subgroup of the TGA family of bZIP transcription factors. Plant Cell 2000, 12, 279–290. [Google Scholar] [CrossRef] [Green Version]
  47. Roetschi, A.; Si-Ammour, A.; Belbahri, L.; Mauch, F.; Mauch-Mani, B. Characterization of an Arabidopsis–Phytophthora pathosystem: Resistance requires a functional PAD2 gene and is independent of salicylic acid, ethylene and jasmonic acid signalling. Plant J. 2001, 28, 293–305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Pré, M.; Atallah, M.; Champion, A.; De Vos, M.; Pieterse, C.M.; Memelink, J. The AP2/ERF domain transcription factor ORA59 integrates jasmonic acid and ethylene signals in plant defense. Plant Physiol. 2008, 147, 1347–1357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Wu, Y.S.; Tang, K.X. MAP kinase cascades responding to environmental stress in plants. Acta Bot. Sin. 2004, 46, 127–136. [Google Scholar]
  50. Bi, G.; Zhou, J.-M. MAP kinase signaling pathways: A hub of plant-microbe interactions. Cell Host Microbe 2017, 21, 270–273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Torres, M.A.; Dangl, J.L.; Jones, J.D. Arabidopsis gp91phox homologues AtrbohD and AtrbohF are required for accumulation of reactive oxygen intermediates in the plant defense response. Proc. Natl. Acad. Sci. USA 2002, 99, 517–522. [Google Scholar] [CrossRef] [Green Version]
  52. Mittler, R.; Vanderauwera, S.; Gollery, M.; Van Breusegem, F. Reactive oxygen gene network of plants. Trends Plant Sci. 2004, 9, 490–498. [Google Scholar] [CrossRef]
  53. Wang, K.L.-C.; Li, H.; Ecker, J.R. Ethylene biosynthesis and signaling networks. Plant Cell 2002, 14, S131–S151. [Google Scholar] [CrossRef] [Green Version]
  54. Guan, R.; Su, J.; Meng, X.; Li, S.; Liu, Y.; Xu, J.; Zhang, S. Multilayered regulation of ethylene induction plays a positive role in Arabidopsis resistance against Pseudomonas syringae. Plant Physiol. 2015, 169, 299–312. [Google Scholar] [CrossRef] [Green Version]
  55. Cabot, C.; Martos, S.; Llugany, M.; Gallego, B.; Tolrà, R.; Poschenrieder, C. A role for zinc in plant defense against pathogens and herbivores. Front. Plant Sci. 2019, 10, 1171. [Google Scholar] [CrossRef]
  56. Wu, S.; Wang, X.; Reddy, U.; Sun, H.; Bao, K.; Gao, L.; Mao, L.; Patel, T.; Ortiz, C.; Abburi, V.L. Genome of ‘Charleston Gray’, the principal American watermelon cultivar, and genetic characterization of 1365 accessions in the US National Plant Germplasm System watermelon collection. Plant Biotechnol. J. 2019, 17, 2246–2258. [Google Scholar] [CrossRef] [Green Version]
  57. Liao, Y.; Smyth, G.K.; Shi, W. featureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef] [Green Version]
  58. Leslie, J.F.; Summerell, B.A. The Fusarium Laboratory Manual; Blackwell Publishing: Ames, IA, USA, 2008. [Google Scholar]
Pathogens 10 00796 g001 550
Figure 1. Scatter plot with differentially expressed genes (DEGs) at 7 and 11 days post-treatment (dpt), compared with non-inoculated micronutrient control (Stnr) in high-Zn (HZ) and low-Zn (LZ) treated plants inoculated with Fusarium oxysporum f. sp. niveum (F) or Meloidogyne incognita (R) or both (F_R). Red dot indicates upregulated genes while green dot indicates downregulated genes.
Figure 1. Scatter plot with differentially expressed genes (DEGs) at 7 and 11 days post-treatment (dpt), compared with non-inoculated micronutrient control (Stnr) in high-Zn (HZ) and low-Zn (LZ) treated plants inoculated with Fusarium oxysporum f. sp. niveum (F) or Meloidogyne incognita (R) or both (F_R). Red dot indicates upregulated genes while green dot indicates downregulated genes.
Pathogens 10 00796 g001
Pathogens 10 00796 g002 550
Figure 2. Gene Ontology enrichment analysis of biological process (BP) for upregulated and downregulated genes between high-Zn (HZ) or low-Zn (LZ) treated plants inoculated with Fusarium oxysporum f. sp. niveum (F) or Meloidogyne incognita (R) or both (F_R) and non-inoculated micronutrient control (Stnr) at 11 days post-treatment (dpt) and 3 days post-inoculation.
Figure 2. Gene Ontology enrichment analysis of biological process (BP) for upregulated and downregulated genes between high-Zn (HZ) or low-Zn (LZ) treated plants inoculated with Fusarium oxysporum f. sp. niveum (F) or Meloidogyne incognita (R) or both (F_R) and non-inoculated micronutrient control (Stnr) at 11 days post-treatment (dpt) and 3 days post-inoculation.
Pathogens 10 00796 g002
Pathogens 10 00796 g003 550
Figure 3. Gene Ontology enrichment analysis of cellular component (CC) for upregulated and downregulated genes between high-Zn (HZ) or low-Zn (LZ) treated plants inoculated with Fusarium oxysporum f. sp. niveum (F) or Meloidogyne incognita (R) or both (F_R) and non-inoculated micronutrient control (Stnr) at 11 days post-treatment (dpt) and 3 days post-inoculation.
Figure 3. Gene Ontology enrichment analysis of cellular component (CC) for upregulated and downregulated genes between high-Zn (HZ) or low-Zn (LZ) treated plants inoculated with Fusarium oxysporum f. sp. niveum (F) or Meloidogyne incognita (R) or both (F_R) and non-inoculated micronutrient control (Stnr) at 11 days post-treatment (dpt) and 3 days post-inoculation.
Pathogens 10 00796 g003
Pathogens 10 00796 g004 550
Figure 4. Gene Ontology enrichment analysis of molecular function (MF) for upregulated and downregulated genes between high-Zn (HZ) or low-Zn (LZ) treated plants inoculated with Fusarium oxysporum f. sp. niveum (F) or Meloidogyne incognita (R) or both (F_R) and non-inoculated micronutrient control (Stnr) at 11 days post-treatment (dpt) and 3 days post-inoculation.
Figure 4. Gene Ontology enrichment analysis of molecular function (MF) for upregulated and downregulated genes between high-Zn (HZ) or low-Zn (LZ) treated plants inoculated with Fusarium oxysporum f. sp. niveum (F) or Meloidogyne incognita (R) or both (F_R) and non-inoculated micronutrient control (Stnr) at 11 days post-treatment (dpt) and 3 days post-inoculation.
Pathogens 10 00796 g004
Pathogens 10 00796 g005 550
Figure 5. Protein–protein interactome map of hormone related genes affected by high-Zn (HZ) and low-Zn (LZ) in plants inoculated with Fusarium oxysporum f. sp. niveum (FON), Meloidogyne incognita (RKN) or both (FON_RKN), compared to non-inoculated micronutrient control (Stnr) at 11 days post-treatment (dpt) and 3 days post-inoculation.
Figure 5. Protein–protein interactome map of hormone related genes affected by high-Zn (HZ) and low-Zn (LZ) in plants inoculated with Fusarium oxysporum f. sp. niveum (FON), Meloidogyne incognita (RKN) or both (FON_RKN), compared to non-inoculated micronutrient control (Stnr) at 11 days post-treatment (dpt) and 3 days post-inoculation.
Pathogens 10 00796 g005
Pathogens 10 00796 g006 550
Figure 6. Protein–protein interactome map of genes involved in MAPK signaling affected by high-Zn (HZ) and low-Zn (LZ) in plants inoculated with Fusarium oxysporum f. sp. niveum (F), Meloidogyne incognita (R) or both (F_R), compared to non-inoculated micronutrient control (Stnr) at 11 days post-treatment (dpt) and 3 days post-inoculation.
Figure 6. Protein–protein interactome map of genes involved in MAPK signaling affected by high-Zn (HZ) and low-Zn (LZ) in plants inoculated with Fusarium oxysporum f. sp. niveum (F), Meloidogyne incognita (R) or both (F_R), compared to non-inoculated micronutrient control (Stnr) at 11 days post-treatment (dpt) and 3 days post-inoculation.
Pathogens 10 00796 g006
Pathogens 10 00796 g007 550
Figure 7. RKN gall formation symptoms and comparative mean nematode counts during RKN disease progression in watermelon roots under high-Zn (HZ) and low-Zn (LZ) conditions in plants inoculated along with Fusarium oxysporum f. sp. niveum (F) or Meloidogyne incognita (R) alone, as compared to non-inoculated control (micronutrient control (Stnr)) at 30dpt. Relative mean nematode counts (egg + J2) were counted in RKN infected samples (three biological replicates). Mean number of RKN adults for each treatment (t) were calculated using Proc Glimmix with Beta distribution for percentage response variable in SAS 9.4 (SAS Institute).
Figure 7. RKN gall formation symptoms and comparative mean nematode counts during RKN disease progression in watermelon roots under high-Zn (HZ) and low-Zn (LZ) conditions in plants inoculated along with Fusarium oxysporum f. sp. niveum (F) or Meloidogyne incognita (R) alone, as compared to non-inoculated control (micronutrient control (Stnr)) at 30dpt. Relative mean nematode counts (egg + J2) were counted in RKN infected samples (three biological replicates). Mean number of RKN adults for each treatment (t) were calculated using Proc Glimmix with Beta distribution for percentage response variable in SAS 9.4 (SAS Institute).
Pathogens 10 00796 g007
Pathogens 10 00796 g008 550
Figure 8. Proposed model of high zinc mediated resistance signaling in watermelon in response to FON and RKN. High level of zinc either in presence of RKN, FON or together induces EDS1 mediated salicylic acid signaling. The increased SA levels further potentiate downstream signaling components (WRKY, zinc finger proteins (ZFP) and map kinase (MAPK)) and trigger nuclear gene expressions of various pathogenesis related (PR1, PR2, PR5) and other defense related genes providing resistance to pathogen.
Figure 8. Proposed model of high zinc mediated resistance signaling in watermelon in response to FON and RKN. High level of zinc either in presence of RKN, FON or together induces EDS1 mediated salicylic acid signaling. The increased SA levels further potentiate downstream signaling components (WRKY, zinc finger proteins (ZFP) and map kinase (MAPK)) and trigger nuclear gene expressions of various pathogenesis related (PR1, PR2, PR5) and other defense related genes providing resistance to pathogen.
Pathogens 10 00796 g008
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.
Sample NameRaw ReadsRaw BasesClean ReadsTotal Mapped ReadsUniquely Mapped ReadsQ20
(%)		Q30
(%)
HighZn_7 dpt5.9 × 1078.935.8 × 107 (98.39%)5.5 × 107 (97.29%)5.2 × 107 (92.49%)98.3194.70
LowZn_7 dpt5.0 × 1077.504.9 × 107 (98.38%)4.7 × 107 (97.52%)4.5 × 107 (92.44%)98.3294.77
Stnr_7 dpt4.9 × 1077.504.9 × 107 (98.26%)4.9 × 107 (97.24%)4.6 × 107 (91.88%)98.2794.67
HighZn_11 dpt5.1 × 1077.765.0 × 107 (98.18%)4.9 × 107 (97.5%)4.7 × 107 (93.07%)98.1594.34
HighZn_F_11 dpt5.3 × 1078.105.3 × 107 (98.52%)5.1 × 107 (97.22%)4.9 × 107 (92.26%)98.2294.54
HighZn_F_R_11 dpt5.3 × 1078.055.2 × 107 (98.68%)5.1 × 107 (96.99%)4.8 × 107 (91.42%)98.3894.94
HighZn_R_11 dpt5.9 × 1079.005.9 × 107 (98.45%)5.6 × 107 (96.75%)5.2 × 107 (91.19%)98.2494.49
LowZn_11 dpt5.1 × 1077.705.0 × 107 (98.57%)4.9 × 107 (97.44%)4.6 × 107 (92.09%)98.3694.87
LowZn_F_11 dpt4.9 × 1077.334.8 × 107 (98.53%)4.6 × 107 (97.53%)4.4 × 107 (92.83%)98.2394.51
LowZn_F_R_11 dpt5.4 × 1078.205.3 × 107 (98.44%)5.0 × 107 (97.12%)4.7 × 107 (91.61%)98.3794.86
LowZn_R_11 dpt4.9 × 1077.434.8 × 107 (98.24%)4.9 × 107 (97.64%)4.6 × 107 (92.97%)98.3294.74
Stnr_11 dpt5.3 × 1077.965.2 × 107 (98.53%)5.0 × 107 (97.17%)4.7 × 107 (91.37%)98.0093.96
Stnr_F_11 dpt4.9 × 1077.474.8 × 107 (98.25%)4.5 × 107 (97.69%)4.3 × 107 (92.41%)98.2494.58
Stnr_F_R_11 dpt5.1 × 1077.675.0 × 107 (98.33%)6.8 × 107 (97.49%)6.5 × 107 (93.26%)98.3694.84
Stnr_R_11 dpt6.3 × 1079.506.2 × 107 (98.31%)4.4 × 107 (96.95%)4.1 × 107 (90.92%)98.3294.76
%Q20: The percentage of bases with phred quality score of 20 or higher; %Q30: The percentage of bases with phred quality score of 30 or higher.
Table
Table 2. Pathways and gene count of major changes in KEGG enrichment in DEGs.
Table 2. Pathways and gene count of major changes in KEGG enrichment in DEGs.
KEGGIDDescriptionp ValueCount
High-Zn_7 dpt
csv04146Peroxisome0.0021652
csv00531Glycosaminoglycan degradation0.016511
High-Zn_FON_11 dpt
csv04075Plant hormone signal transduction1.41 × 10−514
csv04626Plant-pathogen interaction0.00079
Low-Zn_FON_11 dpt
csv00592alpha-Linolenic acid metabolism0.002133
csv00290Valine, leucine and isoleucine biosynthesis0.003632
csv00591Linoleic acid metabolism0.004172
High-Zn_RKN_11 dpt
csv00710Carbon fixation in photosynthetic organisms1.73 × 10−725
csv00195Photosynthesis1.58 × 10−621
csv00030Pentose phosphate pathway0.0001517
csv00500Starch and sucrose metabolism0.0002131
csv00053Ascorbate and aldarate metabolism0.0005215
csv01200Carbon metabolism0.0013848
csv00196Photosynthesis—antenna proteins0.002468
csv00906Carotenoid biosynthesis0.0032711
Table
Table 3. Expression profile of hormone related DEGs.
Table 3. Expression profile of hormone related DEGs.
TreatmentGene IDDescriptionLog2Fold ChangeHormoneName
HighZn_11 dptClCG00G003890Kinase family protein−1.7089Abscissic acidSnRK2/OST1
ClCG01G001770Unknown protein0.86607
ClCG05G004670Auxin-responsive protein−1.02AuxinAUX/IAA/ATAUX2-11
ClCG05G010220Cyclin D3-11.9249BrassinosteroidCYCD3/CYCD3
ClCG05G017790Abscisic acid receptor PYR1−1.8846Abscissic acidPYR/PYL
ClCG09G000190Histidine-containing phosphotransfer protein, putative−1.9525CytokinineAHP/AHP4
ClCG09G009470Auxin response factor1.5998AuxinARF
HighZn_FON_11 dptClCG02G011740Auxin transporter-like protein 31.2816AuxinAUX1
ClCG03G014890Jasmonate-zim-domain protein 10.84193Jasmonic acidJAZ/JAZ1
ClCG05G010220Cyclin D3-11.8242BrassinosteroidCYCD3/CYCD3
ClCG05G010400GH31.7049AuxinGH3
ClCG05G015030GH31.4005AuxinGH3
ClCG05G018490Auxin response factor, putative1.1285AuxinARF
ClCG05G020230Cyclin D1.6576BrassinosteroidCYCD3
ClCG06G001800jasmonate-zim-domain protein 10 LENGTH = 1971.1178Jasmonic acidJAZ
ClCG06G008840Auxin-responsive protein1.1682AuxinAUX/IAA
ClCG07G014490Jasmonate ZIM-domain protein 3b2.033Jasmonic acidJAZ
ClCG08G015300BES1/BZR1 protein1.1164BrassinosteroidBZR1/2
ClCG10G001590Transcription factor MYC2-like protein2.4574Jasmonic acidMYC2
ClCG11G001010Histidine kinase 2−0.87669CytokinineCRE1/WOL
ClCG11G006270Auxin transporter-like protein1.0301AuxinAUX/IAA/AUX1
HighZn_RKN_11 dptClCG01G001770Unknown protein1.132
ClCG01G017190GH3-12.6134AuxinGH3/GH3.1
ClCG01G018130Auxin-responsive family protein1.2014AuxinAUX/IAA/IAA9
ClCG02G004210Auxin transporter-like protein 32.2592AuxinAUX1/LAX2
ClCG02G006220Cyclin D3.2 protein1.7357BrassinosteroidCYCD3
ClCG02G007190Cysteine-rich venom protein−1.3058Salicylic acidPR-1
ClCG02G011740Auxin transporter-like protein 32.6958AuxinAUX1/LAX2
ClCG03G001230Protein phosphatase 2C1.0321Abscissic acidPP2C/PP2CA
ClCG03G014890Jasmonate-zim-domain protein 11.8871Jasmonic acidJAZ/JAZ1
ClCG03G016180Protein phosphatase 2c, putative1.5072Abscissic acidPP2C/ABI1
ClCG04G004170BTB/POZ ankyrin repeat protein2.8402Salicylic acidNPR1/BOP2
ClCG05G000910Abscisic acid receptor PYR1−0.78244Abscissic acidPYR/PYL/RCAR1
ClCG05G002710Kinase family protein0.99557CytokinineB-ARR/BSK2
ClCG05G004370Ethylene-responsive transcription factor 1B, putative1.9806EthyleneERF1/2
ClCG05G008820Protein phosphatase 2c, putative2.2817Abscissic acidPP2C/PP2CA
ClCG05G010220Cyclin D3-13.8147BrassinosteroidCYCD3/CYCD3
ClCG05G010400GH32.8392AuxinGH3/GH3.1
ClCG05G015030GH32.6677AuxinGH3
ClCG05G017790Abscisic acid receptor PYR1−1.8103Abscissic acidPYR/PYL2
ClCG05G018470SAUR-like auxin-responsive protein family LENGTH = 1542.753AuxinSAUR
ClCG05G018490Auxin response factor, putative2.8171AuxinARF/MP
ClCG05G020230Cyclin D3.1068BrassinosteroidCYCD3
ClCG06G001800jasmonate-zim-domain protein 10 LENGTH = 1971.1222Jasmonic acidJAZ/JAZ10
ClCG06G008840Auxin-responsive protein1.7638AuxinAUX/IAA/IAA27
ClCG06G011420Indole-3-acetic acid inducible 29−1.2596AuxinAUX/IAA
ClCG07G000600Auxin-responsive protein−1.535AuxinAUX/IAA/IAA16
ClCG07G005870Jasmonate ZIM-domain protein 21.7448Jasmonic acidJAZ/JAZ3
ClCG07G011440auxin response factor 15 LENGTH = 5983.0858AuxinARF/ARF9
ClCG07G013340Auxin-responsive protein1.671AuxinAUX/IAA/IAA26
ClCG07G014490Jasmonate ZIM-domain protein 3b3.8539Jasmonic acidJAZ/JAZ6
ClCG08G015300BES1/BZR1 protein1.2391BrassinosteroidBZR1/2
ClCG08G016160Cyclin D3.1 protein0.77641BrassinosteroidCYCD3
ClCG09G009470Auxin response factor2.1743AuxinARF/ARF9
ClCG10G001590Transcription factor MYC2-like protein3.1504Jasmonic acidMYC2
ClCG11G003560response regulator 12 LENGTH = 596−1.018AuxinAUX/IAA/RR12
LowZn_FON_11 dptClCG05G015030GH31.389AuxinGH3
ClCG07G014490Jasmonate ZIM-domain protein 3b1.9189Jasmonic acidJAZ
ClCG07G005870Jasmonate ZIM-domain protein 21.5802Jasmonic acidJAZ
ClCG10G001590Transcription factor MYC2-like protein1.7156Jasmonic acidMYC2
LowZn_FON_RKN_11 dptClCG01G001770Unknown protein1.2498
ClCG05G004670Auxin-responsive protein−1.4502AuxinAUX/IAA/ATAUX2-11
ClCG05G010220Cyclin D3-12.3229BrassinosteroidCYCD3/CYCD3
ClCG05G010400GH31.4416AuxinGH3
ClCG05G015030GH31.9852AuxinGH3
ClCG05G020230Cyclin D2.6506BrassinosteroidCYCD3
ClCG07G005870Jasmonate ZIM-domain protein 21.7035Jasmonic acidJAZ
ClCG07G012040Auxin-responsive GH3 family protein−1.1949Auxin/GH3.17
ClCG10G001590Transcription factor MYC2-like protein1.8151Jasmonic acidMYC2
SteineRKN_FON_11 dptClCG05G010400GH31.2947AuxinGH3
ClCG07G014490Jasmonate ZIM-domain protein 3b2.2331Jasmonic acidJAZ
ClCG03G014890Jasmonate-zim-domain protein 11.285Jasmonic acidJAZ/JAZ1
ClCG10G001590Transcription factor MYC2-like protein1.8784Jasmonic acidMYC2
LZ_FON_11 dptClCG05G015030GH31.389AuxinGH3/DFL1
Table
Table 4. Expression profile of MAPK signaling DEGs.
Table 4. Expression profile of MAPK signaling DEGs.
TreatmentGene_idDescriptionLog2 Fold ChangeName
HighZn_11 dptClCG00G003890Kinase family protein−1.7089SnRK2
ClCG05G017790Abscisic acid receptor PYR1−1.8846PYR/PYL
ClCG07G011590Nucleoside diphosphate kinase−1.1574NDPK2
ClCG09G008560Kinase family protein−1.0076OXI1
ClCG11G001650Protein kinase1.2952ANP1
HighZn_FON_11 dptClCG03G000130WRKY transcription factor 071.5388WRKY33
ClCG07G0079001-aminocyclopropane-1-carboxylate synthase1.7681ACS6
ClCG10G001590Transcription factor MYC2-like protein2.4574MYC2
ClCG10G013100Respiratory burst oxidase-like protein1.5519RbohD
ClCG10G022500WRKY 7 transcription factor2.0214WRKY33
HighZn_RKN_11 dptClCG01G020770Chitinase−1.0037ChiB
ClCG02G007190Cysteine-rich venom protein−1.3058PR1
ClCG02G015820Protein kinase, putative0.98462ANP1
ClCG03G000130WRKY transcription factor 072.0512WRKY33
ClCG03G001230Protein phosphatase 2C1.0321PP2C
ClCG03G003540Nucleoside diphosphate kinase−1.1075NDPK2
ClCG03G016180Protein phosphatase 2c, putative1.5072PP2C
ClCG05G000910Abscisic acid receptor PYR1−0.78244PYR/PYL
ClCG05G004370Ethylene-responsive transcription factor 1B, putative1.9806ERF1
ClCG05G008820Protein phosphatase 2c, putative2.2817PP2C
ClCG05G015780Leucine-rich receptor-like protein kinase family protein LENGTH = 9762.7165ER/ERLs
ClCG05G017790Abscisic acid receptor PYR1−1.8103PYR/PYL
ClCG07G0079001-aminocyclopropane-1-carboxylate synthase2.5606ACS6
ClCG09G013300Receptor-like protein kinase2.4712ER/ERLs
ClCG09G015420Copper-transporting atpase paa1, putative−1.1722RAN1
ClCG10G001590Transcription factor MYC2-like protein3.1504MYC2
ClCG10G013100Respiratory burst oxidase-like protein2.6512RbohD
ClCG10G022500WRKY 7 transcription factor2.1625WRKY33
ClCG11G018720Catalase−1.4567CAT1
LZ_FON_11 dptClCG10G001590Transcription factor MYC2-like protein1.7156MYC2
LZ_FON_RKN_11 dptClCG07G011590Nucleoside diphosphate kinase−0.96588NDPK2
ClCG10G001590Transcription factor MYC2-like protein1.8151MYC2
ClCG11G001650Protein kinase1.5157ANP1
Stnr_FONClCG03G000130WRKY transcription factor 071.0872NDPK2
ClCG07G0079001-aminocyclopropane-1-carboxylate synthase1.2221MYC2
ClCG10G001590Transcription factor MYC2-like protein1.8784ANP1
ClCG10G013100Respiratory burst oxidase-like protein1.3053RbohD
ClCG10G022500WRKY 7 transcription factor2.0338WRKY33
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Karki, K.; Coolong, T.; Kousik, C.; Petkar, A.; Myers, B.K.; Hajihassani, A.; Mandal, M.; Dutta, B. The Transcriptomic Profile of Watermelon Is Affected by Zinc in the Presence of Fusarium oxysporum f. sp. niveum and Meloidogyne incognita. Pathogens 2021, 10, 796. https://doi.org/10.3390/pathogens10070796

AMA Style

Karki K, Coolong T, Kousik C, Petkar A, Myers BK, Hajihassani A, Mandal M, Dutta B. The Transcriptomic Profile of Watermelon Is Affected by Zinc in the Presence of Fusarium oxysporum f. sp. niveum and Meloidogyne incognita. Pathogens. 2021; 10(7):796. https://doi.org/10.3390/pathogens10070796

Chicago/Turabian Style

Karki, Kasmita, Tim Coolong, Chandrasekar Kousik, Aparna Petkar, Brendon K. Myers, Abolfazl Hajihassani, Mihir Mandal, and Bhabesh Dutta. 2021. "The Transcriptomic Profile of Watermelon Is Affected by Zinc in the Presence of Fusarium oxysporum f. sp. niveum and Meloidogyne incognita" Pathogens 10, no. 7: 796. https://doi.org/10.3390/pathogens10070796

APA Style

Karki, K., Coolong, T., Kousik, C., Petkar, A., Myers, B. K., Hajihassani, A., Mandal, M., & Dutta, B. (2021). The Transcriptomic Profile of Watermelon Is Affected by Zinc in the Presence of Fusarium oxysporum f. sp. niveum and Meloidogyne incognita. Pathogens, 10(7), 796. https://doi.org/10.3390/pathogens10070796

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