Authors: Yuanxiu Wang aff001; Guihong Xiong aff001; Zhe He aff001; Mingfeng Yan aff001; Manfei Zou aff001; Junxi Jiang aff001 Authors place of work: College of Agronomy, Jiangxi Agricultural University, Nanchang, China aff001; College of Bioscience and Bioengineering, Jiangxi Agricultural University, Nanchang, China aff002 Published in the journal: PLoS ONE 15(1) Category: Research Article doi: https://doi.org/10.1371/journal.pone.0227303
Ripe rot caused by Botryosphaeria dothidea causes extensive production losses in kiwifruit (Actinidia chinensis Planch.). Our previous study showed that kiwifruit variety “Jinyan” is resistant to B. dothidea while “Hongyang” is susceptible. For a comparative analysis of the response of these varieties to B. dothidea infection, we performed a transcriptome analysis by RNA sequencing. A total of 305.24 Gb of clean bases were generated from 36 libraries of which 175.76 Gb was from the resistant variety and 129.48 Gb from the susceptible variety. From the libraries generated, we identified 44,656 genes including 39,041 reference genes, 5,615 novel transcripts, and 13,898 differentially expressed genes (DEGs). Among these, 2,373 potentially defense-related genes linked to calcium signaling, mitogen-activated protein kinase (MAPK), cell wall modification, phytoalexin synthesis, transcription factors, pattern-recognition receptors, and pathogenesis-related proteins may regulate kiwifruit resistance to B. dothidea. DEGs involved in calcium signaling, MAPK, and cell wall modification in the resistant variety were induced at an earlier stage and at higher levels compared with the susceptible variety. Thirty DEGs involved in plant defense response were strongly induced in the resistant variety at all three time points. This study allowed the first comprehensive understanding of kiwifruit transcriptome in response to B. dothidea and may help identify key genes required for ripe rot resistance in kiwifruit.
Gene expression – Fruits – RNA sequencing – Arabidopsis thaliana – MAPK signaling cascades – Calcium signaling – Plant cell walls – Kiwifruit
Kiwifruit is an economically important fruit crop mainly grown in China, New Zealand, and Italy [1]. Ripe rot, caused by Botryosphaeria dothidea, is presently one of the most devastating diseases of kiwifruit in China and abroad, which restricts the sustainable development of kiwifruit industry [2–4]. In severe cases, the disease can cause up to 80% loss in production as occurred during 2011–2012 in Fengxin County, Jiangxi Province, China [5]. B. dothidea is a dominant species of the genus with worldwide distribution and a wide range of hosts. It causes dieback, branch cankers, and fruit rot in hosts including apple, pear, pistachio, and blueberry [6–9]. Fruit infection occurs mostly at the early fruiting stage. However, symptoms on fruit appear only from near maturity to storage, resulting in fruit drop and postharvest decay [2]. As B. dothidea is capable of infecting a large number of plant species and has latent infection features, developing kiwifruit varieties resistant to ripe rot through conventional breeding and biotechnology is considered one of the most effective management strategies.
Studies on the molecular mechanisms of kiwifruit resistance to ripe rot are limited. Furthermore, studies on the interaction between B. dothidea and other hosts are few and preliminary. An earlier study in Malus domestica reported the defensive role of PR4 (pathogenesis-related protein 4) against B. dothidea using RT-qPCR and SDS-PAGE [10]. In addition, Bai et al. reported an increased expression of XEGIP gene encoding xyloglucan-specific endo-(1–4)-beta-D-glucanase inhibitor protein in M. domestica in response to B. dothidea infection [11]. Zhang et al. reported a significant difference in CNGC among apple varieties with different resistance levels to B. dothidea [12]. However, these reports did not give a systematic description of the defense response mechanisms against fungal pathogens.
High-throughput RNA sequencing (RNA-seq) technology is a powerful and efficient method for transcriptome analysis with higher coverage and greater resolution. Researchers use RNA-seq to quantify, profile, and discover RNA transcripts. Studies have used transcriptomics technologies to study host-pathogen interactions, including those between banana and Fusarium oxysporum f. sp. cubense [13], maize and Sporisorium reilianum f. sp. zeae [14], pea and Aphanomyces euteiches [15], and cotton and Verticillium dahlia [16]. Therefore, we used RNA-seq to analyze the transcriptome profile of kiwifruit after B. dothidea inoculation to reveal the interaction mechanism between B. dothidea and kiwifruit.
In the present study, we explore the defense response of a susceptible variety (“Hongyang”, HY) and a resistant variety (“Jinyan”, JY) infected by B. dothidea using RNA-seq. Our findings will help understand the response of kiwifruit to B. dothidea infection and provide new theoretical basis for developing disease resistant variety by genetic engineering.
Two kiwifruit varieties, B. dothidea-susceptible “Hongyang” (HY) and -resistant “Jinyan” (JY), of Shankou kiwifruit orchard in Fengxin county of Jiangxi Province were used. HY fruits harvested at 140 days after flowering and JY fruits harvested at 180 days after flowering were selected for the experiments. A total of 210 B. dothidea strains were isolated from the lesions with the typical symptoms of ripe rot in the infected HY fruits. These strains were cultured at 27°C for 3 days, preserved on potato dextrose agar slants, and maintained in 20% glycerol (-80°C) at the College of Agronomy, Jiangxi Agricultural University (Jiangxi, China). After virulence assessment, “GF27”, the highly pathogenic strain of B. dothidea, was selected for inoculation in Shankou kiwifruit orchard.
B. dothidea strain GF27 was cultured on fresh potato dextrose agar at 27°C for 3 days and mycelial discs of 5 mm in diameter were punched out for inoculation. Healthy and ripe fruits on the trees were surface sterilized with 75% ethanol, peels were allowed to air-dry, and an epidermal tissue of 5 mm in diameter was removed from each fruit. Mycelial disc of B. dothidea was used to inoculate each wound. Control fruits received agar discs lacking mycelium. All treated and control fruits were covered with plastic bags to maintain humidity. We sampled control and treated fruits of the resistant and susceptible varieties for transcriptome analysis at 1, 3, and 6 days after inoculation. The flesh surrounding the discs were collected, frozen in liquid nitrogen, transported to the laboratory on dry ice, and stored at -80°C. Flesh surrounding the discs taken from five different fruits randomly selected from three different trees were polled as a biological replicate. Three independent biological replicates were prepared for each treatment. Samples collected at three different inoculation time points (1, 3, and 6 d) for “Hongyang” (HY, with B. dothidea inoculation; HC, without inoculation) and “Jinyan” (JY, with B. dothidea inoculation; JC, without inoculation) were used to construct thirty six fruit libraries for RNA-seq, which were named as HY11, HY12, HY13, HC11, HC12, HC13, HY31, HY32, HY33, HC31, HC32, HC33, HY61, HY62, HY63, HC61, HC62, HC63, JY11, JY12, JY13, JC11,JC12, JC13, JY31, JY32, JY33, JC31, JC32, JC33, JY61, JY62, JY63, JC61, JC62, and JC63.
Total RNA was isolated from frozen kiwifruits using TRIzol reagent (Invitrogen, CA, USA) following manufacturer’s instructions. We monitored RNA degradation and contamination using 1% agarose gels and assessed RNA integrity using RNA Nano 6000 Assay Kit with an Agilent Bioanalyzer 2100 (Agilent Technologies, CA, USA). The purity was assessed using a NanoPhotomere® UV-Vis spectrophotometer (Implen, CA, USA), and concentration was estimated using Qubit® RNA Assay Kit in Qubit® 2.0 Flurometer (Life Technologies, CA, USA).
Nuclear RNA (3 μg per sample) was used to generate cDNA libraries (n = 36) using NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, USA) according to manufacturer’s protocol. The mRNA was purified from total RNA using oligo (dT) magnetic beads and cleaved into small fragments. First strand cDNA was synthesized using random hexamer primers and M-MuLV reverse transcriptase (RNase H-), and second strand cDNA was synthesized using DNA polymerase I and RNase H. After adenylation of the 3’ ends of the cDNA fragments, NEBNext adapter oligonucleotides were ligated for hybridization. The library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA) to select cDNA fragments of preferentially 150–200 bp in length. The required fragments were enriched by PCR amplification using Phusion high-fidelity DNA polymerase, universal PCR primers, and index (X) primer. PCR products were purified (AMPure XP system), and the library quality was assessed on the Agilent Bioanalyzer 2100 system.
Clustering of the index-coded samples was performed on a cBot cluster generation system using TruSeq PE Cluster Kit v3-cBot-HS (Illumina) according to manufacturer’s instructions. The library preparations were sequenced on an Illumina HiSeq platform to generate 125 bp/150 bp paired-end reads. The raw sequencing data of this study have been deposited in the BIG Data Center GSA database (Accession No. CRA001649).
Raw reads were preprocessed to remove adapter sequences, reads containing ploy-N, and low-quality sequences (Q < 20). Clean reads were mapped to the kiwifruit reference genome (http://bioinfo.bti.cornell.edu/pub/kiwifruit) using TopHat (2.0.12) [17]. Novel transcripts were identified from TopHat alignment results by Cufflinks [18]. We counted the reads that mapped to each gene using HTSeqv0.6.1 (https://pypi.python.org/pypi/HTSeq) [19] and normalized to FPKM (fragments per kilobase of exon model per million mapped reads) using DESeqv1.10.1 software. Differential expression analysis of three biological replicates per treatment was done using DESeq R package (1.10.1) [20]. The resulting P-values were adjusted using Benjamini and Hochberg method [21] to control the false discovery rate. Genes with an adjusted P-value <0.05 found by DESeq were defined as DEGs.
Gene Ontology (GO) enrichment analysis of differentially expressed genes was implemented using GOseq R package [22] in which gene length bias was corrected. GO terms with corrected P-values <0.05 were considered significantly enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG; http://www.http://www.genome.jp/kegg/) pathway analysis of differentially expressed genes was performed by KOBAS software [23]. KEGG pathways with corrected P-values ≤0.05 were considered significantly enriched by DEGs.
To test the reliability of RNA-seq data, eight DEGs from the both varieties were selected for validation by RT-qPCR. Specific primers were designed using Primer 5.0 software (PREMIER Biosoft, Palo Alto, CA, USA) and synthesized by Sangon Biotech (Shanghai, China). Primer pairs are listed in S1 Table. Actin isoform B (ACTB) [24] served as the internal reference gene. The cDNA was synthesized from 4 μL of total RNA using GoScriptTM Reverse Transcription System (Promega, Madison, USA) in 20 μL of reaction mixture. RT-qPCR was performed using GoTaq® qPCR master mix (Promega, Madison, USA) on a Bio-Rad CFX 96 real time PCR system (Bio-Rad, Hercules, CA, USA). The total reaction volume was 20 μL including 0.5 μL F/R primer, 2.0 μL template, and 10 μL master mix. Reaction conditions were as follows: 3 min denaturation at 95°C, followed by 42 cycles of 95°C for 30 s and 60°C for 30 s. Following amplification, melting curve analysis was performed by increasing the temperature from 65°C to 95°C (0.5°C/5 s) to confirm the specificity of the PCR amplification. Three replicates with three technical repeats per experiment were maintained for each gene. A 2−ΔΔCT algorithm was applied for quantitative gene expression analysis [25].
A total of 1,219,912,954 raw reads were generated from 36 libraries with 305.24 Gb high-quality (Q > 20) clean bases that were selected for further analysis. More than 41 million reads were obtained for each library; 52.20%–86.26% of the total mapped reads were aligned onto the kiwifruit reference genome. Approximately, 1.50%–2.46% of the reads were mapped to multiple locations and 50.71%–83.98% were mapped to single locations (Table 1, S1 Fig). The number of mapped reads was not significantly different from the mapped chromosomes (S2 Fig).
All read counts were normalized to FPKM to obtain the relative level of expression. As shown in S2 Table, about 50% of the total number of genes (44,656) had FPKM ≥1 and approximately 5% of the genes had FPKM ≥60 in each library. The high correlation coefficient of the three biological replicates assured the reliability (S3 Fig).
A total of 13,898 DEGs (S3 Table) were detected in “Jinyan” (resistant, JY, R) and “Hongyang” (susceptible, HY, S) samples (Fig 1). In “Jinyan”, 579 (352 upregulated; 227 downregulated), 4,421 (2,680 upregulated; 1,741 downregulated), and 574 (538 upregulated; 36 downregulated) DEGs were detected in the pairwise comparisons JY1 vs. JC1, JY3 vs. JC3, and JY6 vs. JC6, respectively. Upregulated genes were more in number than downregulated genes at all three stages of infection (Fig 2A, S4 Fig). In“Hongyang”, 803 (639 upregulated; 164 downregulated), 1,109 (937 upregulated; 172 downregulated), and 11,998 (5,724 upregulated; 6,274 downregulated) DEGs were detected in the pairwise comparisons HY1 vs. HC1, HY3 vs. HC3, and HY6 vs. HC6, respectively. (Fig 2B, S4 Fig). Upregulated genes were more than downregulated genes at the first and the second time points (1 and 3 days after inoculation), while upregulated genes were less compared to downregulated genes at the third time point (6 days after inoculation). We detected 48 DEGs in R and 207 DEGs in S with sustained expression (Fig 2A and 2B). However, no DEGs were found in both the varieties at all three time points. Interestingly, 36 DEGs in R with sustained expression were detected in S at the third time point (S3 Table).
Hierarchical clustering of the DEGs was done based on log10(FPKM +1) of 36 samples (Fig 2C). The expression of DEGs was different in R and S before and after B. dothidea inoculation. These findings suggest the specific responses of R and S to B. dothidea. Similar expression patterns were found in the inoculated samples JY3 and HY6; distinct expression patterns were found in JY1, JY6, HY1, and HY3.
We performed GO term enrichment analysis (P < 0.05) to identify the functions of DEGs in each pairwise comparison at three time points. Through the analysis of the upregulated DEGs, we observed significantly enriched terms in three major GO categories such as biological process, molecular functions, and cellular component (S4 Table). Oxidation-reduction process (GO: 0055114), oxidoreductase activity (GO: 0016491), dioxygenase activity (GO: 0051213), and oxidoreductase activity (acting on paired donors) (GO: 0016706) associated with oxidation-reduction were found enriched at all three time points in S of which oxidation-reduction process (GO: 0055114) and oxidoreductase activity (GO: 0016491) were exclusively present on the sixth day in R. Polysaccharide metabolic process (GO: 0005976), xyloglucan (xyloglucosyl transferase activity) (GO: 0016762), and apoplast (GO: 0048046) correlated with the cell wall changes were found enriched in JY3 vs. JC3 and JY6 vs. JC6 pairwise comparisons, whereas no GO terms were enriched in JY1 vs. JC1. Based on these findings, we consider cell wall as a major player in kiwifruit in providing resistance against B. dothidea. Concurrently, we analyzed the GO terms of the downregulated DEGs in the inoculated samples of R and S varieties (S4 Table). There was no GO term associated with sustained downregulated genes at all three time points in R or S. Some GO terms (xyloglucan xyloglucosyl transferase activity; GO: 0016762) and apoplast; GO: 0048046) upregulated in R were downregulated in S. The differences in xyloglucan xyloglucosyl transferase activity (GO: 0016762) and apoplast (GO: 0048046) indicate their potential roles in regulating the resistance mechanism to B. dothidea in R.
To further identify the biological pathways in which the DEGs were involved, we performed KEGG analysis. In total, 6,557 DEGs were annotated in the KEGG database and assigned to 113 KEGG pathways. “Metabolic pathway” was the most enriched term and contained 223 DEGs (1.61% of the total DEGs, 13,898) followed by “protein processing” (136, 0.98%) and “carbon metabolism” (124, 0.89%) (S5 Table). A total of 24 pathways were significantly upregulated and 11 pathways were significantly downregulated (S5 Table) when treated samples were compared with the control samples of both the varieties. No pathway was enriched in the downregulated DEGs of HY1 vs. HC1 and JY1 vs. JC1 comparisons. The upregulated plant-pathogen interaction (ath04626) pathway was enriched in the R variety at all three time points (1, 3, and 6 days after inoculation); however, this pathway was enriched in S variety only at the third day after inoculation. Phenylpropanoid biosynthesis (ath00940) DEGs, which play a positive role in plant resistance response, were exclusively downregulated in HY3 vs. HC3. These results indicate that the defense response is more active in R than in S. Notably, photosynthesis (ath00195), photosynthesis-antenna proteins (ath00196), porphyrin and chlorophyll metabolism (ath00860), and carbon fixation in photosynthetic organisms (ath00710) were the critical downregulated pathways in JY3 vs. JC3. This observation confirms the role of photosynthesis-associated pathways in the resistance mechanism [26].
In order to investigate the defense mechanism of kiwifruits against B. dothidea infection, we identified 2,377 potential defense-related genes by searching the keywords in the gene annotation and referring to the literature on defense response. These defense-related genes included 519 PRRs, 32 MAPK, 583 TFs, 83 resistance proteins (R Proteins), 105 pathogenesis-related proteins (PRP), 217 calcium signaling genes, 312 cell wall modification-related genes, and 523 hormone metabolism genes (S6 Table). These DEGs involved in resistance mechanisms against B. dothidea showed distinct expression patterns when R and S varieties were compared. In the R variety, the upregulated DEGs were more than downregulated DEGs at the second time point and the downregulated DEGs were fewer at the third time point. In the S genotype, a similar number of DEGs were identified at the first and second time points. Although the number of upregulated genes increased by the third time point, the number of downregulated genes also increased (Fig 3). Collectively, 30 upregulated genes with sustained expression (PRRs, MAPK, calcium signaling, TFs, hormone metabolism, and cell wall modification) were identified to play roles in kiwifruit resistance against B. dothidea (Fig 4).
To confirm the accuracy of RNA-seq data, eight DEGs in both varieties (Two Protein phosphatase-2C (PP2C) genes (Achn251121 and Achn104901), Axi (auxin-independent growth promoter)-like protein gene (Achn040411), calcium-transporting ATPase gene (Achn012851), COBRA gene (Achn386421), GDSL esterase/lipase gene (Achn372801), mitogen activated protein kinase gene (Achn315051), calmodulin (CaM)-like protein gene (Achn327381)) were selected for validation by RT-qPCR. RT-qPCR results showed the same expression pattern as RNA-seq for these eight DEGs (Fig 4, Fig 5); however, the degree of expression varied between the two data sets because of the difference in sensitivity. These results suggested the reliability of RNA-seq to analyze the transcriptome of resistant and susceptible plants during pathogen infection.
This is the first study to use RNA-seq to identify genes in the resistant and susceptible kiwifruit varieties in response to B. dothidea invasion. A total of 13,898 DEGs were detected between “JY” and “HY” and 2,373 potential defense-related genes were identified by keyword search in gene annotations, which were determined by a literature search on the defense mechanisms against biotic stresses (S6 Table).
In the present study, we observed no significant difference in the expression levels of pattern recognition receptor (PRR) genes between R variety and S variety at the first time point (1 day after inoculation); however, the expression levels in R were more than in S at the second time point (3 days after inoculation) (S6 Table, Fig 4). Two genes encoding receptor-like kinase (Achn240671, Achn325051) sustained upregulated at all three time points in R genotype. In S genotype, Achn240671 was upregulated at the first and the third time points, and Achn325051 was upregulated only at the third time point (S7 Table; Fig 4). First layer of innate immune system in plants is based on a sensitive perception of pathogen or microbe-associated molecular patterns (PAMPs) through pattern recognition receptors (PRRs) at the cell surface, resulting in PAMP-triggered immunity (PTI) that halts further colonization [27]. These DEGs, which are potential candidate genes involved in the response to B. dothidea, have been reported to provoke a rapid immune response in plants [28]. It is speculated that PTI plays a role in kiwifruit resistance to B. dothidea. Achn278911 (PRR) sustained downregulated at all three time points in R variety; however, it was upregulated at the second time point (3 days after inoculation) in S variety. The contrasting expression pattern of Achn278911 (PRR) in R and S varieties emphasize the need to further study their role in kiwifruit resistance to B. dothidea.
In the present study, some DEGs encoding MAPK had different expression patterns in the two varieties. In S genotype, MAPK genes were almost not expressed except Achn135551 at the first and the second time points, and there was no significant difference in the number of upregulated and downregulated genes. Compared with the susceptible variety, MAPK cascades were induced early and almost upregulated at all three time points in R (Fig 3). MAPK cascades are universal signaling modules in eukaryotes that respond to various environmental stresses like cold, drought, and salinity [29]. They play crucial roles in signaling plant immune responses including phytoalexin biosynthesis, stomatal closure, and hypersensitive response (HR) [29, 30]. Three coexpressed genes (Achn091221 encoding mitogen-activated protein kinase kinase kinase (MAP3K), Achn315051 encoding MAPK, and Achn315061 encoding MAPK3) were upregulated in R. Achn315051 (MAPK), Achn315061 (MAPK3), and Achn343991 (encoding a kiwifruit anthranilate N-benzoyltransferase) (S7 Table, Fig 4) involved in the biosynthesis of phytoalexin sustained upregulated in R. In Arabidopsis, MAPK cascades were activated following the activation of WRKY-type transcription during Botrytis cinerea infection, and they regulated phytoalexin production [31]. Taken together, our analysis suggests that MAPK cascades may regulate the biosynthesis of phytoalexin in response to B. dothidea infection.
Calcium is a ubiquitous intracellular messenger, which regulates plant responses to abiotic and biotic stresses like heat, drought, salt, and pathogens [32]. Biotic stress induces increase in intracellular calcium ions (Ca2+) that combine with Ca2+-binding proteins leading to physiological and biochemical responses to pathogen invasion [33, 34]. Recent research has shown that Ca2+ signaling is necessary for stomatal closure, and it can activate cell wall modification following perception of oligogalacturonides (OGs), which are pectic fragments of plant cell wall [35, 36]. In the current study, calmodulin (CaM) gene Achn304251 and two calmodulin-binding family protein (CaMBP) encoding genes Achn199221 and Novel00884 were upregulated in R genotype; however, they were downregulated in S genotype (S6 Table). Seven CaMBP genes (Achn089411 (calmodulin-like protein), Achn201941 (calmodulin-binding protein), Achn327381 (calmodulin-like protein), Achn012851 (calcium-transporting ATPase), Achn124411 (guanylate kinase), Achn017121 (copine) and Achn024491 (src2-like protein)) sustained upregulated in R (S7 Table, Fig 4); however, did not sustained upregulated in S. These upregulated DEGs encoding Ca2+/CaM-binding proteins indicate the role of calcium signaling pathway in defense responses such as stomatal closure, phytoalexin biosynthesis, and accumulation of PR protein. Five genes associated with cell wall modification including Xyloglucan endotransglucosylase/hydrolase (XTH) (Achn367751 and Novel00171), 3-ketoacyl-CoA synthase (KCS) (Achn291611), extensin (Novel05351), and COBRA (Achn386421) (S7 Table, Fig 4) were upregulated at all three time points in R, which implies that Ca2+ might participate in cell wall modification in plant defense.
Transcription factors (TFs) play direct or indirect roles in different cell signaling pathways and regulate plant defense processes [37]. In this study, WRKY, GATA, and PLATZ transcription factors were proven to be involved in coping with pathogens. A total of 22 DEGs encoding WRKY (including WRKY 2, WRKY 3, WRKY 5, WRKY 8, WRKY 14 and WRKY 26) (S6 Table) TFs were upregulated only in R by second and third time points following inoculation with B. dothidea. This is in agreement with the previous observations in Arabidopsis on the role of WRKY 33 in providing resistance to necrotrophic fungal pathogens [38]. In our study, five DEGs encoding GATA TFs were upregulated in R; 14 DEGs encoding GATA TFs were induced in S, of which eight were downregulated at the third time point. Additionally, two GATA genes (Achn166151 and Achn294851) (S7 Table, Fig 4) sustained upregulated in R. The two DEGs encoding PLATZ TFs were induced in R, and one PLATZ gene (Achn020891) sustained downregulated in R.
Both PLATZ and GATA are plant-specific zinc-dependent DNA-binding proteins [39]. There are no reports to prove that GATA and PLATZ play roles in defense responses to pathogens. Previous studies have elucidated the roles of GATA 9, GNC, and CGA1/GNL in low temperature-induced stress responses [40], chlorophyll biosynthesis, and glutamate synthase, respectively, in Arabidopsis [41]. GmPLATZ1 gene was specifically induced by drought, high salinity, and abscisic acid (ABA) in soybean [42]. In the current study, two GATA genes (Achn166151 and Achn294851) sustained upregulated and Achn020891 (GATA) sustained downregulated in R. These observations suggest their roles in the resistance interaction between kiwifruit and B. dothidea.
In order to resist a wide range of microbial pathogens, plants produce various hormones including abscisic acid (ABA), salicylic acid (SA), jasmonic acid (JA), ethylene (ET), auxin, and brassinosteroids (BRs) that regulate physiological processes at low concentrations [43]. In the present study, SA and JA were not significantly different between R and S, which indicate that they may not play a role in defense.
ABA is considered as a negative regulator of plant defense responses against pathogens [44]. Protein phosphatase (PP2C) is often considered as a negative regulatory factor in the ABA signal transduction pathway [45]. However, studies also showed that core ABA signaling components like PP2C and SnRK2 genes were upregulated in Setaria viridis in response to abiotic stresses and proved PP2C as a positive regulator of callose deposition resistance limiting virus spread [46, 47]. In the current study (S6 Table), PYL and SnRK2 were not differentially expressed, while PP2C and bZIP were differentially expressed between the varieties after infection. Achn147981 (bZIP) was upregulated in R; however, it was downregulated in S. Novel04220 (bZIP) was induced to a very high level in R than in S. PP2C genes expressed earlier in R than in S. Five PP2C DEGs (Achn286741, Novel04947, Achn194931, Achn104901, and Achn251121) were detected that showed difference between R and S varieties. Two PP2C genes, Achn286741 and Novel04947, were upregulated in R and downregulated in S. Three PP2C genes, Achn194931, Achn104901, and Achn251121, sustained upregulated at all three time points in R and only at the third time point in S (S7 Table, Fig 4). Consistent with previous reports [47], our experiments revealed significant differences in the expression of ABA signaling-related genes in the two varieties. These findings indicate that ABA signaling pathway plays a significant role in the response to B. dothidea infection, PP2Cs act as positive regulators of ABA signaling pathway. Based on the upregulated Ca2+ and MAPK DEGs in the current study, we speculate that ABA may activate Ca2+ regulation and MAPK cascades to regulate stomatal closure or callose deposition. The three PP2C genes (Achn194931, Achn104901, and Achn251121) may be involved in stomatal closure that helps confine the infection and prevent the spread of B. dothidea.
Consistent with ABA, auxin signaling is also antagonistic to SA signaling and promotes plant susceptibility to the pathogen [48]. However, studies showed that auxin signaling can mediate resistance to Plectosphaerella cucumerina and P. cinnamomi in A. thaliana [49, 50]. Consistent with the previous reports, upregulation of auxin-related DEGs nearly doubled by third day in R, and the downregulation of auxin-related DEGs doubled by sixth day in S. Two auxin-related DEGs (auxin-responsive protein gene Achn030331 and auxin response factor gene Achn154151) were upregulated in R and downregulated in S during B. dothidea infection. The expression of Achn040411 encoding auxin-independent growth promoter increased continuously only in R (S7 Table, Fig 4).
Cyclin-dependent kinase inhibitor (ICK) gene SIAMESE-RELATED1 (SMR1) regulates cell cycle progression and effects growth in transgenic plants [51]. In this study, only four DEGs were induced in both the varieties. An ICK1 gene (Achn266081) sustained upregulated in R (S7 Table, Fig 4). Safae Hamdoun et al. demonstrated the positive role of ICK gene in regulating innate immunity against Pseudomonas syringae in A. thaliana [52]. Present study is consistent with the findings of Safae Hamdoun et al.; however, studies are needed to elucidate the exact role of Achn266081 (ICK1).
We found that five DEGs (S7 Table, Fig 4) involved in cell wall modification including XTH genes (Achn367751 and Novel00171), KCS gene (Achn291611), extensin gene (Novel05351), and COBRA gene (Achn386421) sustained upregulated in R compared with S. Plant cell wall is composed of cellulose, hemicellulose, pectin, and a few structural proteins [53]. It serves as a defense barrier that protects the plant from pathogen penetration. Liu et al. found that differentially accumulated proteins (DAPs) associated with cell wall reinforcement played a role in providing resistance to the necrotrophic fungal pathogen, B. cinerea in kiwifruit [54]. Miedes et al. reported that the reduced expression of xyloglucan endotransglucosylase/hydrolase (XTHs) facilitates Penicillium expansum to infect tomato [55]. KCS gene is involved in very long chain fatty acid (VLCFA) biosynthesis, with a proven role in pathogen resistance in Arabidopsis [56, 57]. Wei et al. demonstrated that Arabidopsis transgenic plants overexpressing extensin gene restrict Pseudomonas syringae invasiveness [58]. Similarly, the mutation in bc1 encoding a COBRA-like (COBL) protein caused cell wall thinness and reduced cellulose content in rice [59]. Consistent with these previous reports, our results indicate that the upregulation of five DEGs, XTHs (Achn367751 and Novel00171), KCS (Achn291611), extensin gene (Novel05351), and COBRA (Achn386421) might prevent B. dothidea infection by thickening cuticular wax, cross-linking extensin monomers, and strengthening plant cell walls.
Cellulose synthase (CesA) gene encodes cellulose enzymes involved in the biosynthesis of cell wall components; it was upregulated in watermelon during defense against Fusarium oxysporum f. sp. niveum [60]. Several structural cell wall proteins and extracellular remodeling enzymes were induced by transient downregulation of CesA genes associated with incompatible interaction of myb46 mutant Arabidopsis plants with B. cinerea [61]. CesA genes Achn299241 sustained upregulated and Achn343521 sustained downregulated in R, which indicate that the two CesA genes respond differently to the presence of B. dothidea. These results imply that cellulose synthase may regulate plant defense response by synthesizing cellulose to reinforce cell wall or by inhibiting cellulose synthesis to activate novel defense pathways in response to B. dothidea.
An earlier work showed that drastic silencing of GDSL1 featured with a conserved GDSL-motif at N-terminus and encoding a member of the esterase/lipase protein reduced cuticle thickness and cutin monomer content in tomato [62]. It was also reported that Arabidopsis GLIP1-elicited systemic resistance to Alternaria brassicicola is dependent on ethylene signaling [63]. In our study, Achn372801 encoding GDSL sustained upregulated at all three time points in R. However, no DEGs involved in ET signaling pathway exhibited sustained upregulation in R. The above analysis revealed that Achn372801 (GDSL) assists resistance to B. dothidea by thickening cuticular wax or disrupting fungal spore integrity in cooperation with signaling of other hormones.
Overall, the mechanism of plant defense response against B. dothidea is complex. Expression of resistance DEGs was detected in R earlier than in S, and the upregulated trend of genes was obvious in R than in S at the third time point. These findings may indicate the resistance specificity and prompt response of the resistant variety.
In this paper, we proposed a putative network underlying the sustained expression of defense-related DEGs in “Jinyan” (Fig 6). PRR proteins were activated by effector proteins, which in turn activated MAPK signaling or calcium signaling. Hormone metabolism and TFs pathways were also activated, which suggest their significant roles in kiwifruit resistance to B. dothidea. As a result, three major defense responses including cell wall modification, stomatal closure, and phytoalexin generation were triggered against B. dothidea.
In this study, we performed a transcriptome analysis to reveal the defense responses of both resistant and susceptible kiwifruit varieties during B. dothidea infection. A total of 305.24 Gb clean bases were generated and 13,898 DEGs were detected in 36 libraries. A total of 2,373 potential defense-related genes were identified; DEGs involved in PRRs, MAPK signaling, calcium signaling, hormone metabolism pathways, TFs pathways, and cell wall modification, which were reported previously as relevant to defense response, were explored and 30 candidate genes related to plant defense response were identified from these pathways. This study provides a better understanding of the molecular basis of defense against B. dothidea in kiwifruit, which may facilitate improvement in disease management via genetic engineering.
S1 Fig [rar] Classification of raw reads.
S2 Fig [rar] Read density on chromosomes.
S3 Fig [png] Pearson correlation between 36 sets of kiwifruit samples.
S4 Fig [rar] Volcano plot of differentially expressed genes during interaction with . .
S1 Table [docx] Primers used for RT-qPCR analysis of differentially expressed genes.
S2 Table [pdf] Number of genes at different FPKM interval after significance level correction.
S3 Table [xls] Annotation of differentially expressed genes.
S4 Table [xls] Significantly enriched GO terms.
S5 Table [xls] Significantly enriched KEGG terms.
S6 Table [xlsx] Expression analysis of kiwifruit defense-related genes.
S7 Table [txt] Transcript sequences mentioned in the paper.
1. Everett KR, Taylor RK, Romberg MK, Rees-George J, Fullerton RA, Vanneste JL, et al. First report of Pseudomonas syringae pv. actinidiae causing kiwifruit bacterial canker in New Zealand. Australas. Plant Pathol. 2011; 6(1): 67–71. https://doi.org/10.1007/s13314-011-0023-9.
2. Li C, Jiang JX, Leng JH, Li BM, Yu Q, Tu GQ. Isolation and identification of pathogenic fungi causing fruit rot of kiwifruit in fengxin county. Jiangxi Nongye Daxue Xuebao. 2012; 34(2): 259–263.
3. Koh YJ, Hur JS, Jung JS. Postharvest fruit rots of kiwifruit (Actinidia deliciosa) in Korea. N. Z. J. Crop Hortic. Sci. 2005; 33(3): 303–310. https://doi.org/10.1080/01140671. 2005.9514363.
4. Pennycook SR, Samuels GJ. Botryosphaeria and Fusicoccum species associated with ripe fruit rot of Actinidia deliciosa (Kiwifruit) in New Zealand. Mycotaxon. 1985; 24 : 445–458.
5. Li C, Wang L, Jiang JX, Wang YX, Zhao SG, Li BM,et al. Biological characteristics of the main pathogenic fungi causing fruit ripe rot of kiwifruit and disease resistance of kiwifruit cultivars. Jiangxi Nongye Daxue Xuebao. 2014; 36(5): 1061–1065. doi: 10.1093/nar/19.2.303.2014169
6. Jurick WM, Ii, Vico I, Gaskins VL, Janisiewicz WJ, Peter KA. First report of Botryosphaeria dothidea causing white rot on apple fruit in Maryland. Plant Dis. 2013; 97(7): 999. https://doi.org/10.1094/PDIS-01-13-0053-PDN 30722551
7. Garibaldi A, Bertetti D, Poli A, Gullino ML. First report of fruit rot in pear caused by Botryosphaeria dothidea in Italy. Plant Dis. 2012; 96(6): 910. doi: 10.1094/PDIS-02-12-0130-PDN 30727390
8. Wunderlich N, Costa SS, Tpoi RP, Ash GJ. First report of Botryosphaeria dothidea causing shoot blight and cankers of pistachio in Australia. Australas. Plant Dis. Notes. 2012; 7(1): 47–49. doi: 10.1007/s13314-012-0045-y
9. Choi IY. First report of bark dieback on blueberry caused by Botryosphaeria dothidea in Korea. Plant Dis. 2011; 95(2): 227. doi: 10.1094/PDIS-05-10-0371 30743439
10. Bai SH, Dong CH, Li BH, Dai HY. A PR-4 gene identified from Malus domestica is involved in the defense responses against Botryosphaeria dothidea. Plant Physiol Biochem. 2013; 62(1): 23–32. doi: 10.1016/j.plaphy.2012.10.016 23178481
11. Bai SH, Dong CH, Zhu J, Zhang YG, Dai HY. Identification of a xyloglucan-specific endo-(1–4)-beta - d-glucanase inhibitor protein from apple (Malus × domestica Borkh.) as a potential defense gene against Botryosphaeria dothidea. Plant Sci. 2015; 231 : 11–9. doi: 10.1016/j.plantsci.2014.11.003 25575987
12. Zhang WW, Dong CH, Zhang YG, Zhu J, Dai HY, Bai SH. An apple cyclic nucleotide-gated ion channel gene highly responsive to Botryosphaeria dothidea infection enhances the susceptibility of Nicotiana benthamiana to bacterial and fungal pathogens. Plant Sci. 2018; 269 : 94–105. doi: 10.1016/j.plantsci.2018.01.009 29606221
13. Li CY, Deng GM, Yang J, Viljoen A, Jin Y, Kuang RB, et al. Transcriptome profiling of resistant and susceptible Cavendish banana roots following inoculation with Fusarium oxysporum f. sp. cubense tropical race 4. BMC Genomics. 2012; 13(1): 374. doi: 10.1186/1471-2164-13-374 22863187
14. Zhang SP, Xiao YN, Zhao JR, Wang FG, Zheng YL. Digital gene expression analysis of early root infection resistance to Sporisorium reilianum f. sp. zeae in maize. Mol Genet Genomics. 2013; 288(1–2): 21–37. doi: 10.1007/s00438-012-0727-3 23196693
15. Hosseini S, Elfstrand M, Heyman F, Dan FJ, Karlsson M. Deciphering common and specific transcriptional immune responses in pea towards the oomycete pathogens Aphanomyces euteiches and Phytophthora pisi. BMC Genomics. 2015; 16(1): 1–18. doi: 10.1186/s12864-015-1829-1 26293353
16. Xu L, Zhu LF, Tu LL, Liu LL, Yuan DJ, Jin L, et al. Lignin metabolism has a central role in the resistance of cotton to the wilt fungus Verticillium dahliae as revealed by RNA-Seq-dependent transcriptional analysis and histochemistry. J Exp Bot. 2011; 62(15): 5607–5621. doi: 10.1093/jxb/err245 21862479
17. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009; 25(9):1105–1111. doi: 10.1093/bioinformatics/btp120 19289445
18. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012; 7(3): 562–578. doi: 10.1038/nprot.2012.016 22383036
19. Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015; 31(2): 166–169. doi: 10.1093/bioinformatics/btu638 25260700
20. Anders S, and Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010; 11(10):R106. doi: 10.1186/gb-2010-11-10-r106 20979621
21. Benjamini Y, Hochberg Y. Controlling the false discovery rate—a practical and powerful approach to multiple testing. Appl Stat. 1995; 57(1): 289–300. https://doi.org/10.2307/2346101
22. Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010, 11(2): R14. doi: 10.1186/gb-2010-11-2-r14 20132535
23. Mao X, Cai T, Olyarchuk JG, Wei L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005; 21(19): 3787–93. doi: 10.1093/bioinformatics/bti430 15817693
24. Jain M, Nijhawan A, Tyagi AK, Khurana JP. Validation of housekeeping genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochem Biophys Res Commun. 2006; 345(2): 646–651. doi: 10.1016/j.bbrc.2006.04.140 16690022
25. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001; 25(4): 402–408. doi: 10.1006/meth.2001.1262 11846609
26. Kangasjarvi S, Neukermans J, Li SC, Aro EM, Noctor G. Photosynthesis, photorespiration, and light signalling in defence responses. J Exp Bot. 2012; 63(4): 1619–1636. doi: 10.1093/jxb/err402 22282535
27. Jones JD, Dangl JL. The plant immune system. Nature. 2006; 444(7117): 323–329. doi: 10.1038/nature05286 17108957
28. Kishimoto K, Kouzai Y, Kaku H, Shibuya N, Minami E, Nishizawa Y. Perception of the chitin oligosaccharides contributes to disease resistance to blast fungus Magnaporthe oryzae in rice. The Plant Journal. 2010; 64(2): 343–354. doi: 10.1111/j.1365-313X.2010.04328.x 21070413.
29. Pathak RK, Taj G, Pandey D, Kumar A, Arora S. Modeling of the MAPK machinery activation in response to various abiotic and biotic stresses in plants by a system biology approach. Bioinformation. 2013; 9(9): 443–449. doi: 10.6026/97320630009443 23847397
30. Meng XZ, Zhang SQ. MAPK cascades in plant disease resistance signaling. Annu Rev Phytopathol. 2013; 51(1): 245–266. doi: 10.1146/annurev-phyto-082712-102314 23663002
31. Ren DT, Liu YD, Yang Kwang-Y, Han L, Mao GH, Glazebrook J, et al. A fungal-responsive MAPK cascade regulates phytoalexin biosynthesis in Arabidopsis. Proceedings of the National Academy of Sciences. 2008; 105(14): 5638–5643. doi: 10.1073/pnas.0711301105 18378893
32. Dodd AN, Kudla J, Sanders D. The language of calcium signaling. Annu Rev Plant Biol. 2010; 61(1): 593–620. doi: 10.1146/annurev-arplant-070109-104628 20192754
33. Lecourieux D, Ranjeva R, Pugin A. Calcium in plant defence-signalling pathways. New Phytol. 2006; 171(2): 249–69. doi: 10.1111/j.1469-8137.2006.01777.x 16866934
34. Reddy VS, Ali GS, Reddy ASN. Characterization of a pathogen-induced calmodulin-binding protein: mapping of four Ca2+-dependent calmodulin-binding domains. Plant Mol Biol. 2003; 52(1): 143–159. doi: 10.1023/a:1023993713849 12825696
35. Allen GJ, Chu SP, Schumacher K, Shimazaki CT, Vafeados D, Kemper A, et al. Alteration of stimulus-specific guard cell calcium oscillations and stomatal closing in Arabidopsis det3 mutant. Science. 2000; 289(5488): 2338–2342. doi: 10.1126/science.289.5488.2338 11009417
36. Moscatiello R, Mariani PD, Sanders D, Maathuis FJM. Transcriptional analysis of calcium-dependent and calcium-independent signalling pathways induced by oligogalacturonides. J Exp Bot. 2006; 57(11): 2847–2865. doi: 10.1093/jxb/erl043 16868046
37. Eulgen T. Regulation of the Arabidopsis defense transcriptome. Trends Plant Sci. 2005; 10(2): 71–78. doi: 10.1016/j.tplants.2004.12.006 15708344
38. Zheng ZY, Qamar SA, Chen ZX, Mengiste T. Arabidopsis WRKY33 transcription factor is required for resistance to necrotrophic fungal pathogens. Plant Journal. 2010; 48(4): 592–605. doi: 10.1111/j.1365-313X.2006.02901.x 17059405
39. Nagano Y, Furuhashi H, Inaba T, Sasaki Y. A novel class of plant-specific zinc-dependent DNA-binding protein that binds to A/T-rich DNA sequences. Nucleic Acids Res. 2001; 29(20): 4097–4105. doi: 10.1093/nar/29.20.4097 11600698
40. Bonthala VS, Mayes K, Moreton J, Blythe M, Wright V, May ST, et al. Identification of gene modules associated with low temperatures response in bambara groundnut by network-based analysis. PLoS One. 2016; 11(2): e0148771. doi: 10.1371/journal.pone.0148771 26859686
41. Hudson D, Guevara D, Yaish MW, Hannam C, Long N, Clarke JD, et al. GNC and CGA1 modulate chlorophyll biosynthesis and glutamate synthase (GLU1/Fd-GOGAT) expression in Arabidopsis. PLoS One. 2011; 6(11): e26765. doi: 10.1371/journal.pone.0026765 22102866
42. So HA, Choi SJ, Chung E, Lee JH. Molecular characterization of stress-inducible PLATZ. Plant Omics Journal. 2015; 8(6): 479–484.
43. Bari R, Jones JDG. Role of plant hormones in plant defence responses. Plant Mol Biol. 2009; 69(4): 473–488. doi: 10.1007/s11103-008-9435-0 19083153
44. Mauch-Mani B, Mauch F. The role of abscisic acid in plant-pathogen interactions. Curr Opin Plant Biol. 2005; 8(4): 409–414. doi: 10.1016/j.pbi.2005.05.015 15939661
45. Wang WM, Ma XF, Zhang Y, Luo MC, Wang GL, Bellizzi M, et al. PAPP2C interacts with the atypical disease resistance protein RPW8.2 and negatively regulates salicylic acid-dependent defense responses in Arabidopsis. Mol Plant. 2012; 5(5): 1125–1137. doi: 10.1093/mp/sss008 22334594
46. Duarte KE, de Souza WR, Santiaqo TR, Sampaio BL, Ribeiro AP, Cotta MG, et al. Identification and characterization of core abscisic acid (ABA) signaling components and their gene expression profile in response to abiotic stresses in Setaria viridis. Sci Rep. 2019; 9 : 4028. doi: 10.1038/s41598-019-40623-5 30858491
47. Seo JK, Kwon SJ, Cho WK, Choi HS, Kim KH. Type 2C protein phosphatase is a key regulator of antiviral extreme resistance limiting virus spread. Sci Rep. 2014; 4 : 5905. doi: 10.1038/srep05905 25082428
48. Wang D, Pajerowska-Mukhtar K, Culler AH, Dong XN. Salicylic acid inhibits pathogen growth in plants through repression of the auxin signaling pathway. Current Biology. 2007; 17(20): 1784–1790. doi: 10.1016/j.cub.2007.09.025 17919906
49. Llorente F, Muskett P, Sanchez-Vallet A, Lopez G, Ramos B, Sanchez-Rodriguez C, et al. Repression of the auxin response pathway increases Arabidopsis susceptibility to necrotrophic fungi. Mol Plant. 2008; 1(3): 496–509. doi: 10.1093/mp/ssn025 19825556
50. Eshraghi L, Anderson JP, Aryamanesh N, Mccomb JA, Shearer B, Hardy GS. Suppression of the auxin response pathway enhances susceptibility to Phytophthora cinnamomi while phosphite-mediated resistance stimulates the auxin signalling pathway. BMC Plant Biol. 2014; 14(1): 68. doi: 10.1186/1471-2229-14-68 24649892
51. Veylder LD, Beeckman T, Beemster GT, Krols L, Terras F, Landrieu I, et al. Functional analysis of cyclin-dependent kinase inhibitors of Arabidopsis. Plant Cell. 2001; 13(7): 1653–1667. doi: 10.1105/TPC.010087 11449057
52. Hamdoun S, Zhang C, Gill M, Kumar N, Churchman M, Larkin JC, et al. Differential roles of two homologous cyclin-dependent kinase inhibitor genes in regulating cell cycle and innate immunity in Arabidopsis. Plant Physiol. 2016; 170(1): 515–527. doi: 10.1104/pp.15.01466 26561564
53. Cosgrove DJ. Growth of the plant cell wall. Nat Rev Mol Cell Biol. 2005; 6(11): 850–861. doi: 10.1038/nrm1746 16261190
54. Liu J, Sui Y, Chen H, Liu Y, Liu Y. Proteomic Analysis of Kiwifruit in Response to the Postharvest Pathogen, Botrytis cinerea. Front Plant Sci. 2018; 9 : 158. doi: 10.3389/fpls.2018.00158 29497428
55. Eva M, Lorences EP. The implication of xyloglucan endotransglucosylase/hydrolase (XTHs) in tomato fruit infection by Penicillium expansum Link. A. J Agric Food Chem. 2007; 55(22): 9021–9026. doi: 10.1021/jf0718244 17960871
56. Todd J, Post-Beittenmiller D, Jaworski JG. KCS1 encodes a fatty acid elongase 3-ketoacyl-CoA synthase affecting wax biosynthesis in Arabidopsis thaliana. Plant J. 2010; 17(2): 119–130. doi: 10.1046/j.1365-313x.1999.00352.x 10074711
57. Raffaele S, Leger A, Roby D. Very long chain fatty acid and lipid signaling in the response of plants to pathogens. Plant Signal Behav. 2009; 4(2): 94–99. doi: 10.4161/psb.4.2.7580 19649180
58. Wei G, Shirsat AH. Extensin over-expression in Arabidopsis limits pathogen invasiveness. Mol Plant Pathol. 2010; 7(6): 579–592. doi: 10.1111/j.1364-3703.2006.00363.x 20507471
59. Li YH, Qian Q, Z YH, Yan MX, Sun L, Zhang M, et al. BRITTLE CULM1, which encodes a COBRA-like protein, affects the mechanical properties of rice plants. Plant Cell. 2003; 15(9): 2020–2031. doi: 10.1105/tpc.011775 12953108
60. Lü GY, Guo SG, Zhang HY, Geng LH, Song FM, Fei ZJ. Transcriptional profiling of watermelon during its incompatible interaction with Fusarium oxysporum f. sp. niveum. Eur J Plant Pathol. 2011; 131(4): 585–601. https://doi.org/10.1007/s10658-011-9833-z
61. Ramírez V, García-Andrade J, Vera P. Enhanced disease resistance to Botrytis cinerea in myb46 Arabidopsis plants is associated to an early downregulation of CesA genes. Plant Signal Behav. 2011; 6(6): 911–913. doi: 10.4161/psb.6.6.15354 21617373
62. Girard AL, Mounet F, Lemaire-Chamley M, Gaillard C, Elmorjani K, Vivancos J, et al. Tomato GDSL1 is required for cutin deposition in the fruit cuticle. Plant Cell. 2012; 24(7): 3106–3121. doi: 10.1105/tpc.112.099796 PMID: 22805434
63. Kwon SJ, Jin HC, Lee S, Nam MH, Chung JH, Kwon S, et al. GDSL lipase-like 1 regulates systemic resistance associated with ethylene signaling in Arabidopsis. Plant J. 2009; 58(2): 235–245. doi: 10.1111/j.1365-313X.2008.03772.x 19077166
Zvýšte si kvalifikáciu online z pohodlia domova
Nemáte účet? Registrujte sa
Zadajte e-mailovú adresu, s ktorou ste vytvárali účet. Budú Vám na ňu zasielané informácie k nastaveniu nového hesla.
