Abstract
Botryosphaeria dieback on the grapevine is caused by Botryosphaeriaceae fungi, which threatens the yield and quality of grapes. At present, chemical control strategies are often observed to be ineffective in controlling the disease worldwide. Improving our understanding of the molecular mechanisms that confer resistance to pathogens would facilitate the development of more pathogen-tolerant grape varieties. Here, we used RNA sequencing analysis to profile the transcriptome of grapevine green shoots infected with Lasiodiplodia theobromae over a time course of 4, 8 and 12 hours post inoculation. A total of 5181 genes were identified as differentially expressed genes (DEGs), and DEGs were more abundant over time. Further analysis revealed that many of these DEGs are involved in plant-pathogen interactions, hormone signal transduction and phenylpropanoid biosynthesis pathways, suggesting that innate immunity, phytohormone signaling and many phenylpropanoid compounds, which constitute a complex defense network in plants, are involved in the response of grapevine against to L. theobromae infection. This study provides novel insights into the molecular mechanisms of plant–pathogen interactions that will be valuable for the genetic improvement of grapevines.
Similar content being viewed by others
Introduction
Grapes is one of the most cultivated fruits around the world, and fungal diseases are very serious threats to its production and quality. Grapevine Botryosphaeria dieback caused by species in the family Botryosphaeriaceae has become an economically important problem in almost all main grapevine growing countries in the last few decades1,2,3,4, with annual losses of more than $260 million USD in the state of California (USA) alone5. More than 20 Botryosphaeriaceae species are associated with grapevine Botryosphaeria dieback, and the most common symptoms are cankers and fruit rot1,6,7. However, the virulence and symptoms caused by a particular Botryosphaeriaceae species can differ greatly between grapevines based on geographical location3,8,9,10,11,12. Botryosphaeriaceae species are reported to infect their host through wounds or natural openings, such as lenticels and stomata13,14,15. These fungi can remain in a latent phase as endophytes without inducing apparent symptoms after entering their hosts and then change to the pathogenic phase as plant opportunistic fungal pathogens to cause disease symptoms when the hosts are stressed16,17,18. It is thus very difficult to determine the infection status of a host and conduct timely control strategies. Recent genomic and transcriptomic analyses have shown that these pathogens have evolved special adaptive mechanisms to infect woody plants15,18. These mechanisms include significant expansions of gene families associated with virulence and nutrient uptake, including cellular transporters, cell wall degrading enzymes (CWDEs), cytochrome P450s, putative effectors and biosynthesis of secondary metabolites. Furthermore, potential virulence functions of CAzymes enriched in Lasiodiplodia theobromae, which play crucial roles in early colonization and pathogenesis, were also revealed. Thus, understanding the molecular responses of grapevine to fungal infection could enable the development of disease-resistant cultivars as an effective and environmentally friendly strategy to control diseases.
To date, our knowledge of host defense responses to pathogens mainly comes from studies in the plant-pathogen model systems, such as Arabidopsis, tomato, tobacco and rice19,20,21,22. These studies have mentioned that the plant defense response is a complex network, including innate immunity, multiple signaling pathways, secondary metabolic products, and cross talks. Nonetheless, information on the host responses to the plant opportunistic fungal pathogens in woody tissues of perennial plants is very limited23. Most studies on woody tissues responding to pathogen infection have focused on the plant’s physical responses or metabolic changes, such as formation of tyloses and gels in xylem24, cell wall modifications with suberin25, production of resins26, and accumulation of peroxidases, superoxide dismutases, glutathione S-transferases, phenolic compounds, stilbenes, and phytoalexins20,27,28,29,30,31,32. Recently, defense responses, for example, cell death, expression of pathogenesis-related (PR) genes, and alkalinization of extracellular medium, were observed to be induced by purified secreted proteins isolated from two Botryosphaeriaceae species, Neofusicoccum parvum and Diplodia seriata, in Vitis suspension cells33. In recent decades, RNA sequencing (RNA-seq) has enabled the study of the transcriptional responses of woody plants to phytopathogens, revealing that genes involved in cell wall biosynthesis, hormone signal transduction and secondary metabolism processes play important roles34,35,36,37.
In the present study, we performed a transcriptome analysis of grapevine (Vitis vinifera) infected by L. theobromae, the most aggressive Botryosphaeriaceae species on grapevines in China3. Our objectives were to gain overall insight into the grapevine-L. theobromae interaction and identify the genes associated with key defense-related responses in grapevine during the early infection stages of L. theobromae. The comprehensive analysis of gene expression in this study will provide a better understanding of the resistance mechanism of grapevine against to Botryosphaeriaceae species and facilitate future research on the application of genomics-assisted breeding systems in grapevine.
Results
Sequencing data summary and prediction of novel transcripts
Lasiodiplodia theobromae successfully infected grapevine green shoots and produced visible lesions within 24 hours after inoculation. However, the defense mechanism of grapevine against fungal diseases is not very clear until now. Therefore, to understand the response of grapevine against L. theobromae during the early infection stages, we used Illumina next-generation sequencing technology (Illumina, San Diego, CA, USA) to analyze the transcriptomes of grapevine green shoots that had been inoculated with L. theobromae after 0 (control), 4, 8, and 12 hours post-inoculation (hpi), with two biological replicates in each of the four treatments (Supplementary Fig. S1; Supplementary Table S1). In total, we generated 135.68 MB reads with an average of 16.96 MB for each library. After strict filtration, ranging from 12.46 to 22.11 MB, the high-quality reads were mapped to the grape genome for further gene expression analysis (Supplementary Table S2). The high-quality reads were deposited in the NCBI SRA database under accession number SRP134225. We optimized the structure of 7,200 genes using the 29,971 annotated genes in the grape reference genome. In addition, we predicted 211 novel transcripts and quantified the expression levels of 13,354 genes (including 177 novel transcripts) across all eight plant samples (Supplementary Table S3).
Global gene expression profiling analysis and identification of differentially expressed genes (DEGs)
The correlation coefficient between two replicate samples was determined based on the FPKM values of genes. The R2 values of the four replicate groups (each representing one time point after inoculation) were 0.955, 0.962, 0.960 and 0.973. These results suggest that the replicate samples exhibited high repeatability and reliability. This result ensures our ability to obtain reliable differences in gene expression analyses between different L. theobromae-inoculated treatments and the 0 hpi sample (Supplementary Fig. S2).
To explore the differences among the transcriptomes induced at different time points after inoculation, we performed separate pairwise comparisons between each of the three L. theobromae-inoculated treatments and the 0 hpi sample. In total, 5181 genes were identified as DEGs with a |log2 fold change| > 2.0 and an adjusted p-value < 0.05 (Fig. 1A, Supplementary Tables S4–S7). Among these genes, 2,177 (1,311 up-regulated and 866 down-regulated), 3453 (1866 up-regulated and 1587 down-regulated) and 4661 (2197 up-regulated and 2464 down-regulated) DEGs were identified in 4 hpi vs 0 hpi, 8 hpi vs 0 hpi and 12 hpi vs 0 hpi, respectively (Fig. 1A). Groups of 245, 142, and 1,386 DEGs were identified as unique to 4 hpi vs 0 hpi, 8 hpi vs 0 hpi and 12 hpi vs 0 hpi, respectively, while 1,702 common genes were regulated at all three time points. A total of 1,706 genes were commonly expressed between at least two time points, with 133 genes between 4 hpi and 8 hpi, 1,476 between 8 hpi and 12 hpi, and 97 between 4 hpi and 12 hpi (Fig. 1B). The number of DEGs increased with time after inoculation. The 12 hpi treatment induced both the greatest number of DEGs and the maximum unique DEGs among the three time points. The unique DEGs expressed at different infection time points may play specific roles at the related infection time point.
To better understand the transcriptome changes in response to L. theobromae infection in grapevine, the expression patterns of all DEGs were analyzed. The expression patterns of all DEGs across the three inoculated time points indicated that the DEGs could be grouped into two main clusters (cluster I and cluster II) and further into six sub-clusters (Fig. 2). Approximately half (50.26%) of the DEGs belonged to cluster I and were down-regulated, whereas the remaining 49.74% of the DEGs, denoted cluster II, were up-regulated after infection (Fig. 2A). Notably, abundant DEGs were sharply induced at 4 hpi.
Gene Ontology classification and Kyoto Encyclopedia of Genes and Genomes analysis of DEGs
The functional categories of all DEGs were analyzed based on Gene Ontology (GO) annotation to gain functional information. In total, 1744 of 2177 DEGs in 4 hpi vs 0 hpi, 2832 of 3453 DEGs in 8 hpi vs 0 hpi, and 3766 of 4661 DEGs in 12 hpi vs 0 hpi were assigned to at least one term of the three major categories: biological process (BP), molecular function (MF) and cellular component (CC) (Supplementary Table S8). Of the 1,702 common DEGs, 1,382 were categorized into at least one GO term, among which 788 were up-regulated, whereas 594 were down-regulated. A large number of the up-regulated common DEGs were associated with catalytic and binding activity in the MF category and were found to participate in metabolic processes, defense responses, and secondary metabolic processes in the BP category.
A Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed to identify the biochemical pathways in which the DEGs were involved in grapevine (organism code “vvi”). The common DEGs were involved in both primary and secondary metabolisms, such as biosynthesis of amino acids (vvi 01230), biosynthesis of secondary metabolites (vvi 01110), carbon fixation in photosynthetic organisms (vvi 00710), carbon metabolism (vvi 01200), circadian rhythm-plant (vvi 04712), flavonoid biosynthesis (vvi 00941), glycolysis/gluconeogenesis (vvi 00010) and ribosome (vvi 03010), which were significantly enriched (adjusted p-value < 0.05). Moreover, several pathways related to disease defense response were also enriched early in infection. Three of these pathways are detailed below.
DEGs involved in the plant-pathogen interaction pathway
According to the KEGG pathway analysis, 83 DEGs identified in grapevine inoculated with L. theobromae were involved in the plant-pathogen interaction pathway, among which 32 were commonly differentially expressed across the time points (Fig. 3; Table S9). Cytoplasmic Ca2+ concentration rapidly accumulates when plants perceive pathogen-associated molecular patterns (PAMPs). The calcium signal is then transmitted by the calcium-dependent protein kinase (CDPK) and calmodulin (CaM)/calmodulin-like proteins (CML) to produce reactive oxygen species (ROS) and nitric oxide (NO) separately through a sequence of biochemical reactions to induce defense responses38. Here, we observed that cyclic nucleotide gated channels (CNGCs), CDPK and CaM/CML were up-regulated after L. theobromae infection, whereas the expression of genes encoding nitric oxide synthase was down-regulated. Typically, once plants recognize PAMPs, a tandem of mitogen-activated protein kinase (MAPK) cascades are activated followed by WRKY transcriptional factors, which prompt expression of defense-related genes21,39. We identified nine DEGs: flagellin-sensing 2 (FLS2), brassinosteroid insensitive 1-associated kinase 1 (BAK1), mitogen-activated protein kinase kinase 1 (MEKK1), mitogen-activated protein kinase 1 (MKK1/2), mitogen-activated protein kinase 4/5 (MKK4/5), mitogen-activated protein kinase 4 (MPK4), mitogen-activated protein kinase 3/6 (MPK3/6), WRKY transcription factor 22 (WRKY22/29), and WRKY transcription factor 33 (WRKY 25/33), respectively. All of these orthologs were reported to be involved in the pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) pathway. Among these DEGs, three MAPKs and two WRKYs were up-regulated upon L. theobromae infection (Fig. 3).
Plants can employ a second layer of defense response to avoid the invasion of pathogens, which is known as the effector-triggered immunity (ETI) pathway21,39. We identified eight DEGs in the ETI pathway upon L. theobromae infection. Among the ETI genes that were up-regulated were one negative regulator PRM1-interacting protein 4 (RIN4), two R genes (RPM1 and RPS2), and two disease-resistance associated genes (SGT1 and HSP90). Conversely, a disease resistance protein (PBS1) and an enhanced disease susceptibility protein (EDS1) were down-regulated.
DEGs involved in the plant hormone signal transduction pathway
Phytohormones act as critical signals for activating disease resistance in plants40. In this study, the signal transduction pathways of some main plant hormones such as auxin, cytokinin, jasmonic acid (JA), salicylic acid (SA), and ethylene (ET) were affected by L. theobromae infection. A total of 127 DEGs involved in the plant hormone signal transduction pathway were detected (Fig. 4, Supplementary Table S10). At 4 hpi, three DEGs in the JA transduction pathway were up-regulated (jasmonate ZIM-domain containing protein (JAZ), jasmonic acid-amino synthetase 1 (JAR1) and transcription factor MYC2). However, JAR1 was down-regulated at 8 hpi and 12 hpi, while JAZ and MYC2 continuously increased over time.
SA is an important plant hormone that regulates plant defense response40, and NPR1 is a key regulator of SA-mediated signal transduction. We found that NPR1 was up-regulated after L. theobromae infection. TGA transcription factors, which activate pathogenesis-related (PR) genes, such as PR-1, showed similar expression patterns with NPR1, whereas PR-1 was down-regulated at 4 hpi and then up-regulated at 8 hpi and 12 hpi. ET is another crucial phytohormone for inducing plant defense response41. We identified seven DEGs in the ET transduction pathway. Ethylene receptor (ETR), ethylene-insensitive protein 3 (EIN3), and ethylene-responsive factors were all up-regulated by L. theobromae. Nevertheless, MPK6 was down-regulated at all three time points. The genes encoding serine/threonine-protein kinase CTR1 were down-regulated at 4 hpi and 8 hpi but up-regulated at 12 hpi. Ethylene-insensitive protein 2 (EIN2) showed the opposite expression pattern of CTR1.
Auxin and cytokinin are associated with plant defense response42. In our study, we identified and analyzed 11 DEGs in the auxin and cytokinin signaling pathways. Ten out of 11 DEGs were up-regulated across all three treatment time points, and one, cytokinin response 1 (CRE1), was down-regulated at 4 hpi and subsequently up-regulated at 8 hpi and 12 hpi (Fig. 4).
DEGs involved in the phenylpropanoid biosynthesis pathway
The phenylpropanoid biosynthesis pathway, which is an important plant secondary metabolic pathway, plays significant roles in plant defense by enforcing or producing physical and chemical barriers against pathogen infection43. We identified 94 DEGs in the phenylpropanoid biosynthesis pathway in response to L. theobromae infection (Supplementary Table S11) and subsequently analyzed these DEGs using KEGG analysis (Fig. 5).
Phenylalanine ammonia-lyase (PAL), which is the first key enzyme in the phenylpropanoid pathway, is known to respond to pathogen infection. Out of a total of 18 transcripts encoding PAL, 16 were up-regulated upon L. theobromae infection. Furthermore, lignin is well characterized and plays an important role in defending against pathogens44. Here, in our data, a large number of enzymes associated with lignin biosynthesis were induced, including transcripts encoding cinnamyl-alcohol dehydrogenase (CAD), peroxidase (PO), 4-coumarate-CoA ligase (4CL), cinnamoyl-CoA reductase (CCR), caffeoyl-CoA O-methyltransferase, ferulate-5-hydroxylase (F5H) and trans-cinnamate 4-monooxygenase (CYP73A) (Fig. 5; Supplementary Table S8).
Genes uniquely expressed in grapevine during L. theobromae infection
The fragment per kilobase per million (FPKM) value of each transcript was estimated to evaluate the gene expression level. The genes of which FPKM was >4 in any of the infected samples and of which FPKM was <1 in the 0 hpi samples were identified as uniquely expressed genes (UEGs) upon L. theobromae infection. A total of 616 UEGs were identified, and we presume that these genes are those that are closely associated with systemic resistance (Fig. 6; Supplementary Table S12).
Validation of the gene expression of DEGs by RT-qPCR
To validate the reliability of the RNA sequencing (RNA-seq) data, the relative expression level of 16 genes was confirmed using RT-qPCR with specific primers after 0, 4, 8, and 12 hpi (Fig. 7; Supplementary Table S13). Ten genes involved in the three key defense pathways (plant-pathogen interaction pathway, plant hormone signal transduction pathway, and phenylpropanoid biosynthesis pathway) and six randomly selected genes (VIT_07s0191g00050, VIT_07s0005g03260, VIT_06s0004g07770, VIT_05s0102g00450, VIT_03s0038g00940 and VIT_01s0010g01970) were used for validation analysis. The relative expression levels of these genes that were obtained from the RT-qPCR assays were compared with our RNA-seq data. The correlation coefficients for the different time points were calculated to assess the correlation between RNA-seq and RT-qPCR analyses. The correlation coefficient values were 0.948, 0.606 and 0.959 for 4, 8 and 12 hpi, respectively, and the correlation was significant at the 0.05 level (Supplementary Fig. S3). The RT-qPCR results showed that the gene expression level of the six randomly selected genes was induced at 4 hpi and reached their highest levels of expression at 8 hpi, whereas RIN4 and JAZ increased in expression over time (Fig. 7). These results further indicated that the plant innate immune response, phytohormone signaling transduction and secondary metabolism products might function together against infection by L. theobromae in grapevine.
Discussion
Plants have evolved complicated signaling and defense pathways in response to pathogens over time. Identification and understanding of such processes will be a platform to develop genomic based breeding for economically important crops. Grapes have become one of the main fruits and beverage crops in China over the last few decades. Among several diseases associated with grapevines, the trunk disease caused by fungi belonging to the family Botryosphaeriaceae has a significant impact as it affects the perennial parts of the grapevine1. To date, the interaction mechanisms between grapevine and Botryosphaeria dieback are still poorly understood. Hence, the objectives of this study were to identify a set of genes responding to the early infection of L. theobromae and to understand signaling networks involved in grapevine. In this study, a total of 5181 DEGs were identified from the RNA-seq analysis performed, and KEGG analysis showed that pathways associated with primary metabolism, secondary metabolism, and defense response were enriched. Therefore, these results indicated that L. theobromae infection stimulated a wide range of responses in the host.
In plants, innate immunity, which consists of PTI and ETI, functions as the first line of defense responses21,39. FLS2 and CERK1, which are involved in the pathogen induced signal perception process, are two well-documented pattern recognition receptors (PRRs) located upstream of PTI. In the present study, CERK1 was up-regulated, whereas FLS2 was down-regulated after L. theobromae infection, indicating that CERK1 may mediate the plant immune pathway to restrain pathogen invasion; however, FLS2 may be recognized by an effector secreted by L. theobromae, leading to the suppression of fls2 expression as L. theobromae releases many kinds of effectors to facilitate its infection process18. BAK1 and FLS2 were previously reported to form a flagellin-induced complex, which is essential for the downstream flg22-signaling responses in Arabidopsis thaliana45. In our study, BAK1 and five members of two independent pathogen-responsive MAPK cascades responded to L. theobromae, suggesting that there are other proteins besides FLS2 that form complexes with BAK1 and mediate downstream MAPK cascades. MEKK1, which is a negative regulator of defense responses, functions together with MKK1/2 and MPK4 to negatively regulate plant innate immune responses46,47,48. MPK3/6, which are downstream of MKK4/5, can act as positive regulators of defense responses49,50. In our study, the MKK1/2 and MPK4 encoding genes were down-regulated, whereas genes encoding MKK4/5 and MPK3/6 were up-regulated, followed by the increased expression of WRKY25/33 and WRKY22/29, respectively, indicating that these two MAPK cascades participated in grapevine responses to L. theobromae infection. The up-regulated expression of VvWRKY33 in grapevine responding to L. theobromae is consistent with the induction of these genes by powdery mildew51 and downy mildew52 in grapevine.
Calcium plays essential roles in regulating plant responses to fungi53. Ca2+ can activate CDPK and Rboh protein, which subsequently induces the production of ROS54. In grapevine, Ca2+ influx was previously proven to be one of the defense reactions when inoculated with Botrytis cinerea55. Our data showed that CNGCs, CDPK and Rboh genes were up-regulated (Fig. 4), indicating that Ca2+ plays crucial roles in signal transduction during the early stages of grapevine defense response.
To infect plants, some pathogens have evolved secreted effectors that suppress PTI56. In response to these effectors, plants employ ETI. Some NB-LRR domain-containing proteins in plants have been identified as receptor proteins that recognize the effector proteins21,39. One receptor protein, RIN4, which is a negative regulator, can be sensed by two R proteins, RPM1 and RPS2, in the form of phosphorylation and degradation57. In grapevine, the RIN4, RPM1 and RPS2 encoding genes were up-regulated in response to L. theobromae infection (Fig. 3). HSP90, SGT1 and RAR1 are important cytosolic chaperones of many R proteins and are involved in the plant innate immunity response by forming a molecular chaperone complex58. We found that HSP90 and SGT1 were up-regulated and RAR1 was down-regulated in grapevine responding to L. theobromae. In addition, protein kinase PBS1, which activates R protein RPS5 via cleavage by an effector59, and EDS1, which is one of the essential components of basal resistance to pathogens, were both down-regulated in response to infection. In general, our data reveal that innate immunity consists of a large network of early responses to infection with L. theobromae.
Secondary metabolites are widespread in plants and contribute to plant disease resistance by functioning as defenses or as signaling molecules60. For example, plant hormones are essential signaling molecules that participate in plant growth, development and defense response processes. In our study, a large number of DEGs were enriched in some plant hormone signal transduction pathways, including the auxin, cytokinin, JA, SA and ET signaling pathways, indicating that these pathways may be involved in the responses of grapevines to L. theobromae infection. SA triggered in grapevine was degraded by enzymes, salicylate hydroxylase and intradiol ring cleavage dioxygenase, produced by L. theobromae61, which could explain why this pathogen successfully infects grapevine despite the genes of the SA signaling pathway continually increased expressed. SA is mainly involved in defending against biotrophic pathogens, whereas JA and ET usually participate in restraining necrotrophic pathogens62,63,64. The antagonism of SA and JA has been reported in rice and Arabidopsis65. For example, SA-response genes were activated by overexpressing OsNPR1, while JA marker genes were suppressed simultaneously in rice66. However, our results showed that most DEGs of the SA and JA signal transduction pathways were up-regulated in grapevine after L. theobromae infection, suggesting that the interaction between grapevine and L. theobromae may have different mechanisms compared with those observed in rice and Arabidopsis. Furthermore, this finding proved that L. theobromae has some unique infection mechanisms compared with other plant pathogens. Auxin signaling is usually a defense response in the plant against necrotrophic fungi, whereas it does not show an obvious relation to the hemibiotrophic pathogens67,68. In our study, all DEGs detected in the auxin transduction pathway were up-regulated, suggesting that grapevine may regulate auxin signaling to defend against L. theobromae. We proposed that L. theobromae secretes effectors into grapevine cells to increase auxin level and finally facilitate its infection, a speculation based on how Pseudomonas syringae invades A. thaliana67. Cytokinins acting as plant growth regulators are also involved in plant defense response. Cytokinins have been shown to induce resistance in A. thaliana against P. syringae by modulating SA signaling69. Recent research has shown that cytokinin signaling is activated by the P. syringae type III effector HopQ1, which leads to the suppression of the Arabidopsis defense response70. Here, in our study, DEGs of the cytokinin signaling pathway were up-regulated, indicating that L. theobromae may interfere with grapevine innate immunity through the secreted effectors.
Phenylpropanoid compounds can preform or induce physical and chemical barriers against pathogen infection, as well as secrete signal molecules that can ultimately induce defense-related gene expression in plants43. In the phenylpropanoid biosynthesis pathway analyzed in this study, most of the genes encoding PAL and enzymes involved in lignin biosynthesis were up-regulated after L. theobromae infection (Fig. 5). These outcomes were consistent with a previous study in which the PAL and 4CL genes of grapevine cells showed increased expression at 6 hpi when treated with Eutypa lata filtrates71. Similar results were also reported in rice, in which all of the transcripts encoding PAL were up-regulated during the early infection stage of Xanthomonas oryzae pv. oryzae72. Moreover, our results were consistent with previous studies of Arabidopsis73 and therefore indicated that grapevine may employ the PB pathway to form lignin or other small molecules to limit L. theobromae infection. Lasiodiplodia theobromae eliminated phenylpropanoid precursors to evade the plant defense response61; thus, it does make sense that grapevine continuously up-regulated genes related to phenylpropanoid compound biosynthesis to prevent infection of this pathogen.
In conclusion, we identified DEGs in the green shoots of grapevine inoculated with L. theobromae and obtained detailed expression profiles of the grapevine response over time. We identified a total of 5181 DEGs, of which 2177 (4 hpi), 3453 (8 hpi) and 4661 (12 hpi) were differentially expressed relative to those of the control (0 hpi). DEGs enriched in plant–pathogen interactions, hormone signal transduction, and the phenylpropanoid biosynthesis pathways were further analyzed. In conclusion, our results suggest that innate immunity, phytohormone signaling and many phenylpropanoid compounds (which constitute a complex defense network in plants) are involved in the response to infection. This study shows - sketchy information of grapevine against L. theobromae, deepens our knowledge of the plant molecular response to opportunistic fungal pathogens and will contribute to a better understanding of the molecular interactions between grapevine and L. theobromae. Additionally, the findings of this study will accelerate the research process on the resistance of grapevine against L. theobromae and will provide a better understanding of grapevine defense responses against plant pathogens.
Materials and Methods
Pathogen infection manipulation
Pathogenicity tests were conducted on 1-year-old green shoots of the grape cultivar ‘Summer Black’. Shoots were pruned to ~30 cm long. To inoculate with L. theobromae, 4-mm diameter wounds were made using a cork borer on each shoot. This procedure was followed by placing a 4 mm agar plug containing the fungal isolate on the wound. The wound with an agar plug was then wrapped with Parafilm (Bemis, Oshkosh, WI, USA). Ten shoots were used as biological replicates. Sterile potato dextrose agar plugs were used as control (designated 0 hpi). Inoculated shoots were planted in sterilized soil in small nursery pots and incubated under moist conditions in an incubation chamber at 25 °C with an 8/16 light/dark cycle.
Samples were collected from the shoots at 0, 4, 8, and 12 hpi. Two pieces of 0.5 cm2 tissues were cut above and below the inoculation sites from each shoot separately using flame-sterilized blades and forceps. The collected tissues were flash-frozen with liquid nitrogen and then immediately stored at −80 °C until RNA was extracted. Two biological replicates were collected at each time point. To obtain enough RNA for subsequent analyses, ten green shoots were pooled as one replicate.
Library construction and RNA-seq
An equal amount of collected tissues from each replicate per time point (0, 4, 8 and 12 hpi) was used for RNA extraction. Total RNA was extracted using TRIzol according to the manufacturer’s instructions (Invitrogen, Carlsbad, CA, USA). RNA integrity, purity, and concentrations were assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) and a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). Purification of messenger RNA (mRNA), construction of complementary DNA (cDNA) libraries, end repair of cDNA, adapter ligation, and cDNA amplification were performed by following methodologies for preparing Illumina RNA-seq libraries of Novogene Bioinformatics Technology Co., Ltd. (Beijing, China). The final quantified libraries were sequenced on an Illumina HiSeq. 2500 (Illumina) to generate paired-end reads of 150 bp.
Sequence alignment
Quality control of reads was assessed by running the FastQC program (v0.11.5), and Trimmomatic (v0.36)74 and low-quality sequences (phred ≤ 20, adapter and poly-N contamination) were filtered out. All subsequent analyses used only the high-quality reads. The reference genome and gene model annotation files were downloaded from the grape genome website (ftp://ftp.ensemblgenomes.org/pub/release-23/plants/fasta/vitis_vinifera/dna/). After building the index of the reference genome using Bowtie v2.2.3 (https://sourceforge.net/projects/bowtie-bio/files/bowtie2/2.3.3/, Broad Institute, Cambridge, MA, USA), paired-end reads were aligned to the reference genome using TopHat v2.0.12 (http://ccb.jhu.edu/software/tophat/manual.shtml, Broad Institute), and HTSeq v0.6.1 (https://files.pythonhosted.org/packages/3c/6e/f8dc3500933e036993645c3f854c4351c9028b180c6dcececde944022992/HTSeq-0.6.1p1.tar.gz, EMBL, Heidelberg, Germany) was used to count the read numbers mapped to each gene. The Cufflinks v2.1.1 (http://cole-trapnell-lab.github.io/cufflinks/) Reference Annotation Based Transcript (RABT) assembly method was used to construct and identify both known and novel transcripts from TopHat alignment results75,76. The reads per kilobase of transcript per million mapped reads (RPKM) of each gene was calculated based on the length of the gene and read count mapped to it. We carried out read alignment and expression quantification separately for each sample. Only genes with fragments per kilobase of transcript per million mapped reads (FPKM) that had a value >4 and that exhibited low variation across two biological replicates (coefficient of variation <30%) were considered reliable and were used in subsequent analyses77.
Identification and functional enrichment analyses of DEGs
The identification of DEGs among the six treatments was performed using the DESeq R package v1.18.0, which provides statistical routines for determining differential expression in digital gene expression data using a model based on the negative binomial distribution. The resulting p-values were adjusted using Benjamini and Hochberg’s approach to control the false discovery rate (FDR)78. Genes with an adjusted p-value < 0.05 found by DESeq. 2 were assigned as differentially expressed.
GO and KEGG pathway analysis
Gene Ontology (GO) enrichment analysis of differentially expressed genes was implemented by the clusterProfiler R package, in which gene length bias was corrected. GO terms with corrected p values less than 0.05 were considered as significantly enriched DEGs. To understand the molecular functions of DEGs, the clusterProfiler R package used to test the statistical enrichment of DEGs in KEGG pathways(http://www.genome.jp/kegg/).
Cluster analysis of gene expression differences
Cluster analysis is used to find genes with similar expression patterns under various experimental conditions. By clustering genes with similar expression patterns, it is possible to discern unknown functions of previously characterized genes or functions of unknown genes. In hierarchical clustering, areas of different colors denote different groups (clusters) of genes, and genes within each cluster may have similar functions or take part in the same biological process. In addition to the FPKM cluster, the H-cluster with complete linkage method was also used to cluster the log2 fold change. Genes within the same cluster exhibit the same trends in expression levels under different conditions.
Real-time qPCR analysis
To validate the RNA-seq results, we performed reverse-transcriptase quantitative PCR (RT-qPCR) using the same biological replicates. We treated 1 µg of RNA with DNase I (Thermo Fisher Scientific) to remove DNA before using SuperScript IV Reverse Transcriptase (Thermo Fisher Scientific) to perform a cDNA synthesis reaction according to the manufacturer’s protocols. A total of 16 genes were selected for qPCR validation of the RNA-seq results. The grapevine actin gene (Vvactin, accession D28123) was used as an endogenous control to normalize the gene expression data. Primers were designed using the Primer3 web tool (v0.4.0, http://bioinfo.ut.ee/primer3-0.4.0/). All primers are listed in Supplementary Table S14. RT-qPCR was performed on an ABI 7500 real-time PCR system (Applied Biosystems, Carlsbad, CA, USA) using SYBR® Premix Ex Taq™ (Tli RNaseH Plus) (Takara Bio Inc., Kusatsu, Shiga, Japan) according to the manufacturer’s instructions. Two biological replicates were used for RT-qPCR analysis for each treatment. At the same time, three technical replicates were used for each biological replicate. The ΔΔCT method was used for relative quantification of gene expression79. The standard deviations (SDs) and correlation coefficients were analyzed by Microsoft Excel 2016 (Microsoft Corporation, Redmond, WA, USA) and IBM SPSS Statistics 21.0 software (IBM China Company Ltd., Beijing, China), respectively.
Data Availability
Supplementary materials are uploaded with the manuscript, and they are available from the corresponding author upon reasonable request.
References
Urbez-Torres, J. R. The status of botryosphaeriaceae species infecting grapevines. Phytopathol. Mediterr. 50, S5––S45 (2011).
Baskarathevan, J., Jaspers, M. V., Jones, E. E. & Ridgway, H. J. Incidence and distribution of botryosphaeriaceous species in new zealand vineyards. Eur. J. Plant. Pathol. 132, 549–560, https://doi.org/10.1007/s10658-011-9900-5 (2012).
Yan, J. Y. et al. Species of botryosphaeriaceae involved in grapevine dieback in china. Fungal. Divers. 61, 221–236, https://doi.org/10.1007/s13225-013-0251-8 (2013).
Kaplana, J. T., Cooperc, M., Hillisd, V., Lubelld, M. & Baumgartner, K. Identifying economic hurdles to early adoption of preventativepractices: The case of trunk diseases in california winegrape vineyards. Wine Economics & Policy 5, 127–141, https://doi.org/10.1016/j.wep.2016.11.001 (2016).
Siebert, J. B. Eutypa: The economic toll on vineyards. Wines Vines. 50–56 (2001).
Phillips, A. J. et al. The botryosphaeriaceae: Genera and species known from culture. Stud. Mycol. 76, 51–167, https://doi.org/10.3114/sim0021 (2013).
Bertsch, C. et al. Grapevine trunk diseases: Complex and still poorly understood. Plant Pathol. 62, 243–265, https://doi.org/10.1111/j.1365-3059.2012.02674.x (2013).
Alves, A., Crous, P. W., Correia, A. & Phillips, A. J. L. Morphological and molecular data reveal cryptic speciation in lasiodiplodia theobromae. Fungal Divers. 28, 1–13 (2008).
Cai, L. et al. The evolution of species concepts and species recognition criteria in plant pathogenic fungi. Fungal Divers. 50, 121–133, https://doi.org/10.1007/s13225-011-0127-8 (2011).
Ko, T. W. K., Stephenson, S. L., Bahkali, A. H. & Hyde, K. D. From morphology to molecular biology: Can we use sequence data to identify fungal endophytes? Fungal Divers. 50, 113–120 (2011).
Ko, T. W. K., Stephenson, S. L., Hyde, K. D. & Lumyong, S. Influence of seasonality on the occurrence of myxomycetes. Chiang Mai J. Sci. 38, 71–84, https://doi.org/10.1007/s13225-011-0130-0 (2011).
Liu, J. K. et al. Towards a natural classification of botryosphaeriales. Fungal Divers. 57, 149–210, https://doi.org/10.1007/s13225-012-0207-4 (2012).
van Niekerk, J. M., Calitz, F. J., Halleen, F. & Fourie, P. H. Temporal spore dispersal patterns of grapevine trunk pathogens in south africa. Eur. J. Plant. Pathol. 127, 375–390, https://doi.org/10.1007/s10658-010-9604-2 (2010).
Urbez-Torres, J. R. & Gubler, W. D. Susceptibility of grapevine pruning wounds to infection by Lasiodiplodia theobromae and Neofusicoccum parvum. Plant Pathol. 60, 261–270, https://doi.org/10.1111/j.1365-3059.2010.02381.x (2011).
Massonnet, M. et al. Neofusicoccum parvum colonization of the grapevine woody stem triggers asynchronous host responses at the site of infection and in the leaves. Front. Plant. Sci. 8, https://doi.org/10.3389/fpls.2017.01117 (2017).
Czemmel, S. et al. Genes expressed in grapevine leaves reveal latent wood infection by the fungal pathogen Neofusicoccum parvum. Plos One 10, e0121828, https://doi.org/10.1371/journal.pone.0121828 (2015).
Chethana, K. W. T., Li, X. H., Zhang, W., Hyde, K. D. & Yan, J. Y. Trail of decryption of molecular research on Botryosphaeriaceae in woody plants. Phytopathol. Mediterr. 55, 147–171 (2016).
Yan, J. Y. et al. Comparative genome and transcriptome analyses reveal adaptations to opportunistic infections in woody plant degrading pathogens of botryosphaeriaceae. DNA Res. 25, 87–102, https://doi.org/10.1093/dnares/dsx040 (2018).
Fawke, S., Doumane, M. & Schornack, S. Oomycete interactions with plants: Infection strategies and resistance principles. Microbiol Mol. Biol. R. 79, 263–280, https://doi.org/10.1128/mmbr.00010-15 (2015).
Sela-Buurlage, M. B. et al. Genome-wide dissection of fusarium resistance in tomato reveals multiple complex loci. Mol. Genet. Genomics 265, 1104–1111, https://doi.org/10.1007/s004380100509 (2001).
Jones, J. D. G. & Dangl, J. L. The plant immune system. Nature 444, 323–329 (2006).
Onaga, G. W. K. Advances in plant tolerance to biotic stresses. InTech:, https://doi.org/10.5772/64351 (2016).
Tobias, P. A. & Guest, D. I. Tree immunity: Growing old without antibodies. Trends Plant Sci. 19, 367–370, https://doi.org/10.1016/j.tplants.2014.01.011 (2014).
Yadeta, K. & Thomma, B. P. H. J. The xylem as battleground for plant hosts and vascular wilt pathogens. Front. Plant. Sci. 4, https://doi.org/10.3389/fpls.2013.00097 (2013).
Pouzoulet, J., Jacques, A., Besson, X., Dayde, J. & Mailhac, N. Histopathological study of response of vitis vinifera cv. Cabernet sauvignon to bark and wood injury with and without inoculation by phaeomoniella chlamydospora. Phytopathol. Mediterr. 52, 313–323 (2013).
Martin, D., Tholl, D., Gershenzon, J. & Bohlmann, J. Methyl jasmonate induces traumatic resin ducts, terpenoid resin biosynthesis, and terpenoid accumulation in developing xylem of norway spruce stems. Plant Physiol. 129, 1003–1018, https://doi.org/10.1104/pp.011001 (2002).
Fossdal, C. G., Sharma, P. & Lonneborg, A. Isolation of the first putative peroxidase cdna from a conifer and the local and systemic accumulation of related proteins upon pathogen infection. Plant Mol. Biol. 47, 423–435 (2001).
Liu, J. J. & Ekramoddoullah, A. K. Characterization, expression and evolution of two novel subfamilies of pinus monticola cdnas encoding pathogenesis-related (pr)−10 proteins. Tree Physiol. 24, 1377–1385, https://doi.org/10.1093/treephys/24.12.1377 (2004).
del Rio, J. A. et al. Phenolic compounds have a role in the defence mechanism protecting grapevine against the fungi involved in petri disease. Phytopathol. Mediterr. 43, 87–94 (2004).
Zamani, A., Sturrock, R. N., Ekramoddoullah, A. K. M., Liu, J. J. & Yu, X. S. Gene cloning and tissue expression analysis of a pr-5 thaumatin-like protein in phellinus weirii-infected douglas-fir. Phytopathology 94, 1235–1243, https://doi.org/10.1094/phyto.2004.94.11.1235 (2004).
Wallis, C. et al. Systemic induction of phloem secondary metabolism and its relationship to resistance to a canker pathogen in austrian pine. New. Phytol. 177, 767–778, https://doi.org/10.1111/j.1469-8137.2007.02307.x (2008).
Spagnolo, A. et al. Differential responses of three grapevine cultivars to botryosphaeria dieback. Phytopathology 104, 1021–1035, https://doi.org/10.1094/phyto-01-14-0007-r (2014).
Stempien, E. et al. Secreted proteins produced by fungi associated with Botryosphaeria dieback trigger distinct defense responses in Vitis vinifera and Vitis rupestris cells. Protoplasma 255, 613–628, https://doi.org/10.1007/s00709-017-1175-z (2017).
Camps, C. et al. A transcriptomic study of grapevine (Vitis vinifera cv. Cabernet-sauvignon) interaction with the vascular ascomycete fungus eutypa lata. J. Exp. Bot. 61, 1719–1737, https://doi.org/10.1093/jxb/erq040 (2010).
Barakat, A. et al. Chestnut resistance to the blight disease: Insights from transcriptome analysis. Bmc Plant Biol. 12, https://doi.org/10.1186/1471-2229-12-38 (2012).
Xu, L. et al. Characterization of the Pinus massoniana transcriptional response to bursaphelenchus xylophilus infection using suppression subtractive hybridization. Int. J. Mol. Sci. 14, 11356–11375, https://doi.org/10.3390/ijms140611356 (2013).
Mangwanda, R. Myburg, A. A. & Naidoo, S. Transcriptome and hormone profiling reveals eucalyptus grandis defence responses against chrysoporthe austroafricana. Bmc Genomics 16, https://doi.org/10.1186/s12864-015-1529-x (2015)
Cheval, C., Aldon, D., Galaud, J. P. & Ranty, B. Calcium/calmodulin-mediated regulation of plant immunity. Bba-Mol. Cell. Res. 1833, 1766–1771 (2013).
Dodds, P. N. & Rathjen, J. P. Plant immunity: Towards an integrated view of plant-pathogen interactions. Nat. Rev. Genet. 11, 539–548, https://doi.org/10.1038/nrg2812 (2010).
Pieterse, C. M. J., V der Does, D., Zamioudis, C., Leon-Reyes, A. & Van Wees, S. C. M. Hormonal modulation of plant immunity. Annu. Rev. Cell Dev. Bi. 2012(28), 489–521 (2012).
van Loon, L. C., Geraats, B. P. J. & Linthorst, H. J. M. Ethylene as a modulator of disease resistance in plants. Trends Plant Sci. 11, 184–191, https://doi.org/10.1016/j.tplants.2006.02.005 (2006).
Fu, J. & Wang, S.P. Insights into auxin signaling in plant-pathogen interactions. Front. Plant. Sci. 2, https://doi.org/10.3389/fpls.2011.00074 (2011).
Dixon, R. A. et al. The phenylpropanoid pathway and plant defence - a genomics perspective. Mol Plant Pathol. 3, 371–390, https://doi.org/10.1046/j.1364-3703.2002.00131.x (2002).
Bhuiyan, N. H., Selvaraj, G., Wei, Y. & King, J. Role of lignification in plant defense. Plant Signal Behav. 4, 158–159, https://doi.org/10.4161/psb.4.2.7688 (2009).
Delphine, C. et al. A flagellin-induced complex of the receptor FLS2 and BAK1 initiates plant defence. Nature 448, 497–500, https://doi.org/10.1038/nature05999 (2007).
Nakagami, H., Soukupova, H., Schikora, A., Zarsky, V. & Hirt, H. A mitogen-activated protein kinase kinase kinase mediates reactive oxygen species homeostasis in arabidopsis. J. Biol. Chem. 281, 38697–38704, https://doi.org/10.1074/jbc.m605293200 (2006).
Suarez-Rodriguez, M. C. et al. Mekk1 is required for flg22-induced mpk4 activation in arabidopsis plants. Plant Physiology 143, 661–669, https://doi.org/10.1104/pp.106.091389 (2007).
Gao, M. H. et al. Mekk1, mkk1/mkk2 and mpk4 function together in a mitogen-activated protein kinase cascade to regulate innate immunity in plants. Cell Res. 18, 1190–1198, https://doi.org/10.1038/cr.2008.300 (2008).
Asai, T. et al. Map kinase signalling cascade in arabidopsis innate immunity. Nature 415, 977–983, https://doi.org/10.1038/415977a (2002).
Pedley, K. F. & Martin, G. B. Identification of mapks and their possible mapk kinase activators involved in the pto-mediated defense response of tomato. J. Biol. Chem. 279, 49229–49235, https://doi.org/10.1074/jbc.m410323200 (2004).
Welter, L. J. Genetic and molecular analysis of mildew disease resistance in grapevine. University Karlsruhe (KIT), Karlsruhe (2008).
Merz, P. R. et al. The transcription factor vvwrky33 is involved in the regulation of grapevine (Vitis vinifera) defense against the oomycete pathogen plasmopara viticola. Physiol. Plantarum 153, 365–380, https://doi.org/10.1111/ppl.12251 (2015).
Zipfel, C. Early molecular events in pamp-triggered immunity. Curr. Opin. Plant Biol. 12, 414–420, https://doi.org/10.1016/j.pbi.2009.06.003 (2009).
Sagi, M. & Fluhr, R. Superoxide production by plant homologues of the gp91(phox) nadph oxidase. Modulation of activity by calcium and by tobacco mosaic virus infection. Plant Physiol. 126, 1281–1290, https://doi.org/10.1104/pp.126.3.1281 (2001).
Poinssot, B. et al. The endopolygalacturonase 1 from Botrytis cinerea activates grapevine defense reactions unrelated to its enzymatic activity. Mol. Plant-Microbe Interact. 16, 553–564, https://doi.org/10.1094/MPMI.2003.16.6.553 (2003).
Boller, T. & He, S. Y. Innate immunity in plants: An arms race between pattern recognition receptors in plants and effectors in microbial pathogens. Science. 324, 742–744, https://doi.org/10.1126/science.1171647 (2009).
Kim, M. G. et al. Two pseudomonas syringae type iii effectors inhibit rin4-regulated basal defense in arabidopsis. Cell 121, 749–759, https://doi.org/10.1016/j.cell.2005.03.025 (2005).
Liu, Y. L., Burch-Smith, T., Schiff, M., Feng, S. H. & Dinesh-Kumar, S. P. Molecular chaperone hsp90 associates with resistance protein n and its signaling proteins sgt1 and rar1 to modulate an innate immune response in plants. J. Biol. Chem. 279, 2101–2108, https://doi.org/10.1074/jbc.m310029200 (2004).
Sun, J. H. et al. Comparative study of arabidopsis pbs1 and a wheat pbs1 homolog helps understand the mechanism of pbs1 functioning in innate immunity. Sci. Rep. 7
Wink, M. Evolution of secondary metabolites from an ecological and molecular phylogenetic perspective. Phytochemistry 64, 3–19, https://doi.org/10.1038/s41598-017-05904-x (2003).
Paolinelli-Alfonso, M. et al. Global transcriptional analysis suggests Lasiodiplodia theobromae pathogenicity factors involved in modulation of grapevine defensive response. BMC Genomics 17, 615, https://doi.org/10.1186/s12864-016-2952-3 (2016).
Glazebrook, J. Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Annu. Rev. Phytopathol. 43, 205–227, https://doi.org/10.1146/annurev.phyto.43.040204.135923 (2005).
Spoel, S. H., Johnson, J. S. & Dong, X. Regulation of tradeoffs between plant defenses against pathogens with different lifestyles. P. Natl. Acad. Sci. USA 104, 18842–18847, https://doi.org/10.1073/pnas.0708139104 (2007).
Bari, R. & Jones, J. Role of plant hormones in plant defence responses. Plant Mol. Biol. 69, 473–488, https://doi.org/10.1007/s11103-008-9435-0 (2009).
Robert-Seilaniantz, A., Grant, M. & Jones, J. D. G. Hormone crosstalk in plant disease and defense: More than just jasmonate-salicylate antagonism. Annu. Rev. Phytopatho. 49, 317–343, https://doi.org/10.1146/annurev-phyto-073009-114447 (2011).
Yuan, Y. X. et al. Functional analysis of rice NPR1-like genes reveals that OsNPR1/NH1 is the rice orthologue conferring disease resistance with enhanced herbivore susceptibility. Plant Biotechnol. J. 5, 313–324, https://doi.org/10.1111/j.1467-7652.2007.00243.x (2007).
Chen, Z. Y. et al. Pseudomonas syringae type iii effector avrrpt2 alters arabidopsis thaliana auxin physiology. P. Natl. Acad. Sci. USA 104, 20131–20136, https://doi.org/10.1073/pnas.0704901104 (2007).
Llorente, F. et al. A. Repression of the auxin response pathway increases arabidopsis susceptibility to necrotrophic fungi. Mol. Plant. 1, 496–509, https://doi.org/10.1093/mp/ssn025 (2008).
Choi, J. et al. The cytokinin-activated transcription factor arr2 promotes plant immunity via tga3/npr1-dependent salicylic acid signaling in arabidopsis. Dev. Cell 19, 284–295, https://doi.org/10.1016/j.devcel.2010.07.011 (2010).
Hann, D. R. et al. The pseudomonas type iii effector hopq1 activates cytokinin signaling and interferes with plant innate immunity. New Phytol. 201, 585–598, https://doi.org/10.1111/nph.12544 (2014).
Mutawila, C. et al. Response of Vitis vinifera cell cultures to Eutypa lata and Trichoderma atroviride culture filtrates: expression of defence-related genes and phenotypes. Protoplasma 254, 863–879, https://doi.org/10.1007/s00709-016-0997-4 (2017).
Cheng, X. J. et al. Transcriptome analysis confers a complex disease resistance network in wild rice oryza meyeriana against Xanthomonas oryzae pv. Oryzae. Sci. Rep. 6, https://doi.org/10.1038/srep38215 (2016).
Ali, M. B. & McNear, D. H. Induced transcriptional profiling of phenylpropanoid pathway genes increased flavonoid and lignin content in arabidopsis leaves in response to microbial products. Bmc Plant. Biol. 14–84, https://doi.org/10.1186/1471-2229-14-84 (2014).
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: A flexible trimmer for illumina sequence data. Bioinformatics 30, 2114–2120, https://doi.org/10.1093/bioinformatics/btu170 (2014).
Trapnell, C. et al. Differential gene and transcript expression analysis of rna-seq experiments with tophat and cufflinks. Nat. Protoc. 7, 562–578, https://doi.org/10.1038/nprot.2012.016 (2012).
Roberts, A., Pimentel, H., Trapnell, C. & Pachter, L. Identification of novel transcripts in annotated genomes using RNA-seq. Bioinformatics 27, 2325–2329, https://doi.org/10.1093/bioinformatics/btr355 (2011).
Trapnell, C. et al. Transcript assembly and quantification by rna-seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511–515, https://doi.org/10.1038/nbt.1621 (2010).
Benjamini, Y. & Yosef, H. Controlling the false discovery rate - a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society 57, 289–300 (1995).
Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative pcr and the 2 -δδ c t method. Methods 25, 402–408, https://doi.org/10.1006/meth.2001.1262 (2001).
Acknowledgements
This work was supported by the China Agriculture Research System (CARS-29), the technological innovation ability project of Beijing Academy of Agriculture and Forestry Sciences (KJCX20170412) and National Natural Science Foundation of China (31501597).
Author information
Authors and Affiliations
Contributions
W. Zhang and W. Zhao conceived and designed the experiments; W. Zhang and J. Yan performed the experiments; W. Zhao and J. Yan analyzed the data; X. Li, Q. Xing and T. Chethana contributed reagents/materials/analysis tools; W. Zhang wrote the paper.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Zhang, W., Yan, J., Li, X. et al. Transcriptional response of grapevine to infection with the fungal pathogen Lasiodiplodia theobromae. Sci Rep 9, 5387 (2019). https://doi.org/10.1038/s41598-019-41796-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-41796-9
This article is cited by
-
Transcriptome analysis of resistant and susceptible grapes reveals molecular mechanisms underlying resistance of white rot disease
Horticulture Advances (2023)
-
Single and combined effects of Pythium oligandrum Po37 and a consortium of three rhizobacterial strains on Sclerotinia stem rot severity and tomato growth promotion
Journal of Plant Pathology (2022)
-
Transcriptomic analysis of Dubas bug (Ommatissus lybicus Bergevin) infestation to Date Palm
Scientific Reports (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.