#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Simultaneous transcriptome analysis of oil palm clones and Phytophthora palmivora reveals oil palm defense strategies


Authors: Kelly Avila-Mendez aff001;  Ávila Rodrigo aff001;  Leonardo Araque aff001;  Hernán Mauricio Romero aff001
Authors place of work: Biology and Breeding Program, OiI Palm Research Center, Cenipalma, Bogotá, Colombia aff001;  Department of Biology, Universidad Nacional de Colombia, Bogotá, Colombia aff002
Published in the journal: PLoS ONE 14(9)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0222774

Summary

Phytophthora palmivora is an oomycete that causes oil palm bud rot disease. To understand the molecular mechanisms of this disease, palm clones with contrasting responses (Ortet 34, resistant and Ortet 57, susceptible) were inoculated with P. palmivora, and RNAseq gene expression analysis was performed. The transcriptome was obtained by sequencing using Illumina HiSeq2500 technology during the asymptomatic phase (24, 72 and 120 hours postinfection, hpi). A simultaneous analysis of differentially expressed gene (DEG) profiles in palm and P. palmivora was carried out. Additionally, Gene Ontology (GO) and gene network analysis revealed differences in the transcriptional profile of the two ortets, where a high specificity of the pathogen to colonize the susceptible ortet was found. The transcriptional analysis provided an overview of the genes involved in the recognition and signaling of this pathosystem, where different transcription factors, phytohormones, proteins associated with cell wall hardening and nitrogen metabolism contribute to the resistance of oil palm to P. palmivora. This research provides a description of the molecular response of oil palm to P. palmivora, thus becoming an important source of molecular markers for the study of genotypes resistant to bud rot disease.

Keywords:

Gene expression – Plant pathogens – Pathogens – Transcription factors – DNA transcription – Oil palm – Vegetable oils

Introduction

The vegetable oil obtained from the cultivation of oil palm (Elaeis guineensis Jacq.) is produced at the lowest cost due to the high yield of the palms. Oil palm can produce 5 to 10 times more oil than other oil crops [1]. The production of palm oil originates mainly from 5 countries: Malaysia, Indonesia, Thailand, Nigeria and Colombia; in 2016, Malaysia and Indonesia accounted for 77% of the world's share of palm oil production [2].

Oil palm crops are affected by various diseases. In Latin America, the main cause of crop losses is the oomycete Phytophthora palmivora [3, 4]. According to Cenipalma´s (Colombian Oil Palm Research Center) Plant Health Management Program, in 2017, 75,000 hectares were affected by the disease, which caused losses of more than USD 1,800 million in the four oil palm-growing areas of the country.

Fungi and oomycetes can cause the most devastating diseases in agriculture [5]. Among the oomycetes, the genus Phytophthora is very important because some of the worst crop diseases are caused by members of this genus [6, 7]. Phytophthora species have a hemibiotrophic life cycle, in which the first part of their development is biotrophic (the cells stay alive and the plant does not show symptoms) followed by a necrotrophic phase in which the pathogen has internally colonized host cells and the plant shows symptoms due to the death of the cells [8].

P. palmivora is a hemibiotrophic organism that causes diseases in various tropical crops, such as papaya, cocoa, coconut and durian, among others [9]. Since its identification in 2008 as the causal agent of oil palm bud rot, studies have focused on the development of inoculation methods in immature leaflets where the initial symptoms of the disease have been replicated [10, 11]. However, the underlying basis of the molecular response of oil palm to bud rot disease has not been studied. Understanding the molecular bases of the oil palm-P. palmivora interaction will not only allow a description of the mechanisms associated with the interaction but will also provide the possibility of developing mechanisms for the early detection of resistant materials in oil palm breeding programs.

The main objective of this research was to obtain information on the molecular mechanisms associated with oil palm resistance to P. palmivora. A transcriptomic strategy was developed to elucidate the mechanisms of palm resistance and pathogen infection. This information will refine the search for genotypes resistant to the disease and understanding of how the pathogen attacks the palm. This work provides the first description of the molecular mechanisms by which oil palm responds to infection by P. palmivora.

Materials and methods

Plant material

Oil palm clones were produced by a somatic embryogenesis technique in the Tissue Culture Laboratory of the Colombian Oil Palm Research Center (Cenipalma); the process to obtain clones takes three years. The produced shoots were then placed on rooting media (MS media) for four months. At this point, the clones were inoculated with P. palmivora. At the time of inoculation, the clones were 10 cm tall and had three to five leaves.

The original explants were obtained from resistant (Ortet 34) and susceptible (Ortet 57) Elaeis guineensis donor palms. Ortet 34 was a Deli x Cameroon palm selected from the oil palm-growing area of Puerto Wilches (Santander, Colombia), an area that had suffered an outbreak of bud rot disease. Ortet 57 was an Angola pisifera palm selected from the Palmar de la Vizcaína Experimental Field.

The contrasting responses of the oil palm clones were determined using previously performed screening tests [12]. These tests showed statistically significant differences between the two ortets. Additionally, the behavior of the clonal materials was consistent with the behavior of the donor palms in the field.

Pathogen

Cultures of an isolate of P. palmivora from a collection of the Plant Pathology Laboratory of the Colombian Oil Palm Research Center–Cenipalma- (CPPhZC-05) were used. The isolate was inoculated on Petri dishes in clarified V8 juice medium (20% V8 juice, 5 g/l CaCO3, 50 mg/l B-sitosterol, 1.5% agar, rifampicin 1 μg/mL). The pathogen was kept under dark conditions for the first 3 days at 25ºC. After three days, the Petri dishes were placed under a 12-hour photoperiod until the cultures completed 10 days of growth.

Inoculation

The inoculations were made with 10-day-old cultures of P. palmivora. The in vitro clones were inoculated at the petiole base with 20 μL of a zoospore suspension at a concentration of 3x106 zoospores/mL. The release of zoospores was performed by thermal shock with sterile distilled water at 4ºC, and the quantification of the pathogen was carried out using a Neubauer Chamber.

The plants were kept in a Weiss® growth chamber at a temperature of 28ºC with a photoperiod of 12 hours of light (72% relative humidity) and 12 hours of dark (75% relative humidity). Eight individuals of the susceptible ortet (Ortet 57) and eight individuals of the resistant ortet (Ortet 34) were inoculated per replication (two biological replicates). Three individuals per ortet were used for the negative control (mock inoculation), which were inoculated with 20 μL of sterile distilled water.

Extraction of total RNA

Tissue from the petiole base of the palm clones was collected from 8 individuals per ortet per biological replicate for each time of infection (24, 72 and 120 hours post infection). RNA was extracted using the Ambion® (AM1912) kit (Ambion Inc., Foster City, Ca), following the manufacturer’s instructions.

For each ortet and time of infection, RNA pools were made for the construction of libraries. The libraries were prepared using the Truseq Stranded mRNA LT Sample Kit, and the samples were sequenced using the HiSeq2500 sequencing platform (126 bp paired-end sequences).

Bioinformatics analysis

The cleaning of the reads obtained from raw data was performed by removing i) low quality reads (Q≥30) and ii) adapter sequences using Trimmomatic software (Version 0.36).

After the quality parameters of the reads of each sample were confirmed, the reads were mapped against the reference genome of Elaeis guineensis using STAR software (Version 5.3a). The bam files for subsequent analyses were obtained using Samtools (Version 1.6–4).

The reads of inoculated samples not aligned to the palm genome were mapped against a P. palmivora reference genome (http://www.ncbi.nlm.nih.gov/bioproject/558734, obtained from a high quality genome assembly based on long and short read sequencing [13]. The platforms Illumina (Illumina, Inc., San Diego, CA) and PacBio (Pacific Biosciences of California, Inc.) were used for sequencing. 2x150 bp paired-end libraries were obtained for Illumina; and for PacBio, libraries were designed to obtain 50X sequencing depth. The PacBio long reads were assembled de novo using the software Canu v1.7. The alignment of the Illumina reads was carried out with bowtie2 v2.3.4.1. Following the assembly of the reference genome, its structural annotation was performed using the tool MAKER.

Differential expression analysis was performed using the statistical package DeSeq 1.18.0–1 (Debian Package). The statistical parameters were an FDR (false discovery rate) cutoff of 0.05 and a fold change (FC) Log2> | 2 |. The identification of unique and overlapping genes of the differentially expressed gene dataset was performed using R-Venn 3.4.3–1 and R-Bioconductor 2.38.0–1.

For the functional analysis, the lists of differentially expressed genes (DEGs) were analyzed with Blast2GO, which generated a GO list for later analysis of networks.

For the network analyses we used the GO list with its p-value for each ortert and infection time and these data were exported to Revigo (REVIGO is freely available at http://revigo.irb.hr/). The information was loaded to the platform and using the interactive graph interphase the .xgmml files were obtained (these files can be used offline). The subsequent visualization of each network was performed using Cytoscape software Version 3.6.0 (https://cytoscape.org/). Networks were based on a color scale (preloaded by default in Cytoscape), the most significant nodes were highlighted, where the red nodes had greater control and/or contribution to the network formed.

Gene validation by qRT-PCR

Twenty genes for the oil palm and four genes for P. palmivora were selected from the dataset to be validated by qRT-PCR. RNA was obtained from the inoculated and negative control clones of the contrasting ortets. First, the DNA was removed using the Ambion® DNAase-free Kit (Cat No AM1906), and the synthesis of the cDNA was carried out using the SuperScript IV First-Strand Synthesis System (Cat No 18901050), following the manufacturer’s instructions. The concentration of total RNA required for cDNA synthesis was 1 μg.

The qRT-PCR reactions were normalized with the GADPH (glyceraldehyde 3- phosphate dehydrogenase) gene for the oil palm qRT-PCR reactions, and a ribosomal protein gene was used for P. palmivora. For both genes, the expression value was verified against the obtained dataset. The primers used for the GADPH gene were those previously reported [14], and for the pathogen, the primer is listed in S1 Table. The qRT-PCR reactions were performed using a LightCycler 480 from Roche ®.

Fold change (FC) values of the qRT-PCR reactions were calculated using the Livak method [15]. RNA-seq FC, and qRT-PCR values were compared by linear regression analysis. The significance of the data was confirmed by Pearson`s rank correlation coefficient (after verifying the normality of the data) using Stata 14 software.

Results

Bioinformatics analysis—Mapping

After the cleaning of raw data, all the samples for each replicate showed a sufficient quantity of reads for the analysis (Table 1); an average of 120 million reads (clean data) was obtained for all samples, and more than 85% of the reads had a Q> 30, ensuring the quality of the sequencing process and the quality of the final reads obtained for the subsequent differential expression analysis.

Tab. 1. Reads, quality and mapping percentage of transcriptomic data.
Reads, quality and mapping percentage of transcriptomic data.

The percentage of mapping performed with the E. guineensis genome and mapping values over 85% were obtained for all samples (Table 1), where 27000 of the 37000 genes of the genome were mapped and used in the subsequent differential expression analyses.

Differential expression analysis

After the read quality verification process and the mapping of each sample were completed, the reads were used for differential expression analysis, which was performed independently for the plant reads and the pathogen reads.

The differential expression analysis (Fig 1) summarizes the expression profile found in the pathosystem evaluated. Principal component analysis (PCA) showed a differential distribution of the genes found per ortet and time of infection (Fig 1A), indicating that for each ortet and time of infection, the plant responded differently to the interaction with the pathogen.

Differential expression profile for pathosystem oil palm clon- <i>P</i>. <i>palmivora</i>.
Fig. 1. Differential expression profile for pathosystem oil palm clon- P. palmivora.
A: Component Principal Analysis (PCA). B: Heatmap for oil palm clones C: Heatmap for P. palmivora expressing during infection process. Ortet 57: susceptible. Ortet 34: resistant.

Heatmaps showed different sets of gene clusters for the plant (Fig 1B), with differential expression patterns for each ortet, indicating that once the interaction with the pathogen started, the signaling cascades between the two ortets were different for the same gene.

Regarding the pathogen gene expression profile (Fig 1C), a large cluster was detected, indicating that the genes used by the pathogen to infect a susceptible ortet and a resistant ortet exhibited a completely different expression pattern.

The Venn diagram analysis showed that out of the 37,000 genes in the palm genome, approximately 8,000 genes were expressed as unique genes for each ortet and time of infection for the conditions of this biological model (Fig 2).

Venn diagram for unigenes expressing differentially in contrasting ortets.
Fig. 2. Venn diagram for unigenes expressing differentially in contrasting ortets.
A. 24 hpi; B. 72 hpi; C. 120 hpi. Ortet 57: susceptible. Ortet 34: Resistant.

The unique genes were used to perform the functional and network analyses to identify the different molecular strategies used by resistant or susceptible ortets to respond to infection by the pathogen.

Functional annotation

The genes obtained from the differential expression analysis dataset were used to carry out functional annotation for each main GO group: Molecular Function (GO: 0003674); Biological Process (GO: 0008150) and Cell Component (GO: 0005575).

Regarding molecular function, in the three infection periods, a different GO distribution was found in the evaluated ortets; some of these were found as unique GOs for each material; for example, at 24 hours of infection, genes in the transferase category were found in the resistant ortet, while genes in the kinase category were found in the susceptible ortet (Fig 3).

GO distribution for molecular function of DEG in oil palm clones during interaction process with <i>P</i>. <i>palmivora</i>.
Fig. 3. GO distribution for molecular function of DEG in oil palm clones during interaction process with P. palmivora.
Top panel Ortet 34 (Resistant), Bottom panel Ortet 57 susceptible.

At 72 hours post infection, the susceptible ortet maintained an equal distribution of its genes for this GO category, while the resistant ortet showed new categories of genes such as genes associated with the expression of kinase-like proteins, production of enzymes from oxidoreductase activity, and ion-binding proteins. Finally, at 120 hpi, a similar distribution in the GO of the molecular function was found in the two ortets.

The GO associated with biological processes included the stress response category, which was found in both ortets. The main differences for this GO were related to other processes. Thus, at 24 hours post infection (hpi), the susceptible ortet showed a category associated with the metabolism of carbohydrates, while in the resistant material, the genes associated with the expression of nitrogen proteins were active at 72 hpi, a pathway that remained active until 120 hpi (Fig 4).

GO distribution for biological process of DEG in oil palm clones during interaction process with <i>P</i>. <i>palmivora</i>.
Fig. 4. GO distribution for biological process of DEG in oil palm clones during interaction process with P. palmivora.
Top panel Ortet 34 (Resistant), Bottom panel Ortet 57 susceptible.

Finally, regarding the GO of the cellular component, a higher proportion of genes in the cell membrane was found for the susceptible ortet, which was maintained at 24 and 72 hpi, while for the resistant ortet, there was a greater distribution of the location of the genes in the different organelles in the 3 infection periods (Fig 5). These results suggest that from the first hours of infection, the two ortets perceived the pathogen, but molecular communication channels were established only with the susceptible ortet, which explains the differences observed in each category analyzed.

GO distribution for cellular component of DEG in oil palm clones during interaction process with <i>P</i>. <i>palmivora</i>.
Fig. 5. GO distribution for cellular component of DEG in oil palm clones during interaction process with P. palmivora.
Top panel Ortet 34 (Resistant), Bottom panel Ortet 57 susceptible.

Network analysis

The distribution of the genes in the different GO categories made it possible to visualize the main differences between the susceptible and resistant ortets. To describe the palm-P. palmivora interaction model, we decided to generate networks associated with the GOs of the biological process (this category included the response to stress). In this way, it was possible to identify differences in the defense responses of the evaluated ortets.

At 24 hpi, both ortets showed differences in the networks. For the resistant ortet, two main networks were identified with three principal nodes (Fig 6A) associated with the macromolecule and cellular metabolism, while in the susceptible ortet, the networks did not have principal nodes established during this infection period (Fig 6B).

Interaction network of DEG in oil palm clones at 24 hpi of infection.
Fig. 6. Interaction network of DEG in oil palm clones at 24 hpi of infection.
A: Ortet 34 resistant B: Ortet 57 susceptible.

At 72 hpi, the network configuration changed significantly for both ortets, where principal nodes were detected. However, each ortet formed different molecular networks, which demonstrates the contrasting response of the ortets to P. palmivora.

The resistant ortet (Fig 7A) developed networks associated with the metabolism of the phenylpropanoid pathway, protein metabolism, protein phosphorylation and phosphate-containing compound metabolism. In contrast, the susceptible ortet (Fig 7B) showed a more complex network, with a principal node related to the response to stress, which was connected with a large number of nodes, one of them related to the negative regulation of the response to stimulus.

Interaction network of DEG in oil palm clones at 72 hpi of infection.
Fig. 7. Interaction network of DEG in oil palm clones at 72 hpi of infection.
A: Ortet 34 resistant B: Ortet 57 susceptible.

At the end of the biotrophic phase, each ortet continued to show large differences in the networks activated. The resistant ortet (Fig 8A) finished this phase with four networks, and a new one was formed related to homeostatic regulation. For the susceptible ortet (Fig 8B), only one network was present: the process related to the response to stress. However, principal nodes were not found at this time of infection.

Fig. 8.
Interaction network of DEG in oil palm clones at 120 hpi of infection A: Ortet 34 resistant B: Ortet 57 susceptible.

Gene validation

The selected genes for oil palm are reported in Table 2. These genes had a relationship with some of the different defense pathways in plants. For example, genes related to the reinforcement of the cell wall, such as chorismate synthase (CHO), wall-associated receptor kinase (WAK),

Tab. 2. List of validated genes for q-RTPCR.
List of validated genes for q-RTPCR.

The susceptible ortet showed overexpression of genes such as acidic endochitinase (AEL), pathogenesis-related protein (PR-1) and glucan endo-1,3-beta-glucosidase (GEB), genes that although involved in defense responses, failed to trigger an effective response against the pathogen advance. Additionally, at 24 hpi, this ortet showed overexpression of Mildew Resistance Locus O (MLO), a gene associated with susceptibility.

The genes selected for P. palmivora were: a gene related to sporulation process and, three genes related to pathogenesis, elicitin, protease and thioredoxin. In the four selected genes, the FC was over 2 when P. palmivora was inoculated on the susceptible clones. However, the FC of the same genes did not reach the FC threshold of 2 when the pathogen interacted with resistant clones (S1 Fig)

A Pearson correlation was conducted for the genes validated for oil palm and the pathogen. The obtained qRT-PCR values were correlated with the FC values obtained from the DEG data set. The Pearson correlation for the oil palm genes was 0.83, and for the pathogen genes, it was 0.87

Discussion

In this study, Illumina (HiSeq2500) technology was used to explore and describe the contrasting responses of oil palm clones to inoculation with P. palmivora. Currently, the molecular mechanisms of the defense response in the oil palm-P. palmivora interaction process are not known. The next two sections will include the analysis related to oil palm defense and P. palmivora pathogenicity strategies.

Oil palm molecular response

Previous studies conducted by our research team found that the biotrophic (asymptomatic) phase of the disease in oil palm clones lasts approximately 120 hours [12]. The transcriptomic analysis showed differences between the two evaluated ortets, where a coordinated series of molecular responses between the palm and P. palmivora were observed.

With the GO enrichment analyses, we found differences in the response to the pathogen between the evaluated ortets. Thus, the differential expression of genes encoding kinase-type proteins found in the susceptible ortet and in the pathogen for the same infection period (24 hpi) suggests the establishment of signaling cascades necessary for successful colonization, as in other pathosystems in which effectors that interact with kinase-type proteins have been reported to block the defense response of plants [16].

An interesting result was TIR1 (TRANSPORT INHIBITOR RESPONSE 1) overexpression in the susceptible ortet (Fig 9), a response that has been observed in other pathosystems in which TIR overexpression has been related to pathogen susceptibility [17, 18]. In fact, auxins have been related to susceptibility to biotrophic pathogens [19], with TIR1 playing a major role. Thus, the reduction of TIR1 transcript levels through posttranscriptional gene silencing or downregulation via the salicylic acid pathway [20], increases pathogen resistance [21, 22].

qRT-PCR validation of candidate genes.
Fig. 9. qRT-PCR validation of candidate genes.
Each graph represents the fold change obtained by qRT-PCR, including the ortet and its inoculation time according to expression profile of each gene. CHO: Chorismate synthase. WAL: Wall-associated receptor kinase. CAF: caffeoylshikimate esterase. WKRY: WRKY transcription factor 40. ACO: 1-aminocicloprppane-1-1carboxylate oxidase. GATA: GATA transcription factor. MYB: MYB transcription factor. RGA: resistance gene analog. AEL: acidic endochitinase-like. TIR: transport inhibitor response 1 like. CALLO: callose synthase. IER1: immediate early response. ETI: ethylene insensitive 3. CTP: CTP synthase. BHB Transcription factor bHLB. MLO. DIS: disease resistance protein. PHOS: 1-phospatidylinositol-3-3phosphate 5 kinase. GEB: glucan endo-1,3-beta glucosidase-like. PR-1: pathogen related.

In fact, gene expression profiles and the gene networks of the resistance ortet at every point after infection suggest the presence of a salicylic acid-mediated defense response through the phenylpropanoid pathway and aromatic compounds [23, 24]. The Chorismate synthase (CHO) gene was overexpressed at the three studied times after infection, indicating that this resistance mechanism was present at the times the pathogen should be in the biotrophic phase, and as a result, the pathogen could not advance or cause the disease.

In the resistant ortet, a whole pathway of nitrogen compounds was upregulated, which together with the phenylpropanoid pathway and aromatic compounds suggests a defense pathway associated with cell wall reinforcement as a first defense barrier, and in the late infection processes, it could be complemented by defense pathways controlled by transcription factors such as MYB and WRKY (Fig 9).

It is important to note that the cell wall reinforcement detected in the transcriptional profile of the resistant ortet may not be associated with the degradation of the membrane components of the pathogen, as suggested by the gene expression profile of glucan endo-1,3-beta-glucosidase (GEB). This gene was found to be overexpressed (Fig 9) in the susceptible ortet, which reinforces the hypothesis that in our pathosystem, P. palmivora uses a set of highly specialized pathogenic molecules to infect oil palm tissues.

At 120 hours post infection, the resistant ortet was not only able to control the entry of the pathogen but also ensured the continuity of the redox state of the cells [12], suggesting that the stability of the redox state is a defense mechanism consistent with the distribution of genes found in the cellular component, where redox-related genes are not only expressed at the membrane level but also in the different mitochondria and plastid organelles associated with the defense response [25].

P. palmivora molecular attack

The inoculation method used in this work did not include a physical open wound on the plant; therefore, the pathogen itself needed to find ways to colonize the tissue. This model shows that once the pathogen started its infection process, the relationship with the plant was different in the susceptible and resistant ortets. Somehow, the pathogen recognized that the ortets were different and released different mechanisms of pathogenicity.

The analyses carried out suggest that P. palmivora, in its infection process, manipulated all the transcriptional machinery of the susceptible ortet (S2 and S3 Figs) by developing fine molecular communication channels. According to reports in the literature, the successful infection process carried out by the pathogen is probably due to a greater control of the secreted effectors [2628]; categories such as transmembrane transport activity and kinases indicate that the pathogen colonized the susceptible ortet tissue from the first hours of infection (24 hpi). The expression of these genes is consistent with the expression of genes such as MLO genes (Fig 9).

It was interesting to note that when P. palmivora attempted to colonize the resistant ortet, it continuously used genes that encode hydrolase-type proteins, and the gene distribution found in the membrane and transmembrane transport activity (S3 Fig) suggests that the pathogen recognized the membrane components of the susceptible material [29, 30].

The presence of this type of protein in the susceptible ortet suggests that it may be a response to the stimulus created by the PAMP-type molecules of the pathogen [31, 32], with which it is able to exploit the molecular mechanisms and colonize the tissue to continue its infectious process to obtain primary and secondary metabolites according to the metabolic capacity reported for oomycetes [24, 33], which is consistent with the gene expression profile of the susceptible ortet in terms of the biological process component.

The biological processes found in the gene expression of the pathogen (S3 Fig) showed how, with small differences in their transcriptional profile, the pathogen is capable of colonizing the tissue of the susceptible ortet and blocking its defense response mechanism.

The analysis of gene expression performed for the biotrophic phase of the disease shows a fine molecular communication between the palm and the pathogen, and due to a coordinated series of events, the pathogen successfully colonizes the susceptible ortet, which itself fails to develop an effective defense mechanism, while the resistant ortet blocks the spread of the pathogen.

Supporting information

S1 Table [docx]
Gene ID, primers, annealing temperature and amplicon size of genes used for qRT-PCR validation.

S1 Fig [tiff]
qRT-PCR Validation of candidate genes for . .

S2 Fig [tiff]
GO distribution for molecular function of DEG in . .

S3 Fig [tiff]
GO distribution for biological process of DEG in . .


Zdroje

1. Vijay V, Pimm SL, Jenkins CN, Smith SJ. The Impacts of Oil Palm on Recent Deforestation and Biodiversity Loss. PloS one. 2016;11(7):e0159668. doi: 10.1371/journal.pone.0159668 Vijay2016. 27462984

2. Fedepalma. Statistical Yearbook 2018. Bogota, Colombia: Federación Colombiana de Cultivadores de Palma de Aceite; 2018. 208 p.

3. Ma Latifah. Cross-infectivity of oil palm by Phytophthora spp. isolated from perennial crops in Malaysia. For Pathol. 2017;47(6):1–6. doi: 10.1111/efp.12374 Latifah2017.

4. Sundram S, Intan-Nur AMA. South American Bud rot: A biosecurity threat to South East Asian oil palm. Crop Protect. 2017;101:58–67. doi: 10.1016/j.cropro.2017.07.010 Sundram2017.

5. Franceschetti M, Maqbool A, Jiménez-Dalmaroni MJ, Pennington HG, Kamoun S, Banfield MJ. Effectors of filamentous plant pathogens: commonalities amid diversity. Microbiol Mol Biol Rev. 2017;81(2):e00066–16. doi: 10.1128/MMBR.00066-16 28356329

6. Coaker G. Pathogen Specialization. science. 2014;343(6170):496–7. doi: 10.1126/science.1250171 24482473

7. Kong P, Hong C. Soil bacteria as sources of virulence signal providers promoting plant infection by Phytophthora pathogens. Scientific reports. 2016;6:33239. doi: 10.1038/srep33239 27616267

8. Chen X-R, Huang S-X, Zhang Y, Sheng G-L, Li Y-P, Zhu F. Identification and functional analysis of the NLP-encoding genes from the phytopathogenic oomycete Phytophthora capsici. Mol Genet Genomics. 2018:1–13.

9. Evangelisti E, Rey T, Schornack S. Cross-interference of plant development and plant-microbe interactions. Curr Opin Plant Biol. 2014;20:118–26. doi: 10.1016/j.pbi.2014.05.014 Evangelisti2014. 24922556

10. Sarria G, Martinez G, Varon F, Drenth A, Guest D. Histopathological studies of the process of Phytophthora palmivora infection in oil palm. Eur J Plant Pathol. 2016;145(1):39–51. doi: 10.1007/s10658-015-0810-9 Sarria2016.

11. Torres G, Sarria G, Martinez G, Varon F, Drenth A, Guest D. Bud Rot Caused by Phytophthora palmivora: A Destructive Emerging Disease of Oil Palm. Phytopathology. 2016;106(4):320–9. doi: 10.1094/PHYTO-09-15-0243-RVW Torres2016. 26714102

12. Ávila-Méndez K, Avila-Diazgranados R, Pardo A, Herrera M, Sarria G, Romero HM. Response of in vitro obtained oil palm and interspecific OxG hybrids to inoculation with Phytophthora palmivora. For Pathol. 2019;0(0):e12486. doi: 10.1111/efp.12486

13. Gil-Bedoya J. Genomic variability of Phytophthora palmivora isolates from different oil palm cultivation regions in Colombia [M.Sc.]. Bogota: Universidad de los Andes; 2019.

14. Chan PL, Rose RJa. Evaluation of reference genes for quantitative real-time PCR in oil palm elite planting materials propagated by tissue culture. PLoS ONE. 2014;9(6). doi: 10.1371/journal.pone.0099774 Chan2014. 24927412

15. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-$\Delta$ $\Delta$CT method. Methods. 2001;25(4):402–8. doi: 10.1006/meth.2001.1262 Livak2001. 11846609

16. Sharpee WC, Dean RA. Form and function of fungal and oomycete effectors. Fungal Biology Reviews. 2016;30(2):62–73. doi: 10.1016/j.fbr.2016.04.001 Sharpee2016.

17. Kidd BN, Kadoo NY, Dombrecht B, Tekeoglu M, Gardiner DM, Thatcher LF, et al. Auxin signaling and transport promote susceptibility to the root-infecting fungal pathogen Fusarium oxysporum in Arabidopsis. Mol Plant-Microbe Interact. 2011;24(6):733–48. doi: 10.1094/MPMI-08-10-0194 21281113

18. Llorente F, Muskett P, Sánchez-Vallet A, López G, Ramos B, Sánchez-Rodríguez C, et al. Repression of the auxin response pathway increases Arabidopsis susceptibility to necrotrophic fungi. Molecular plant. 2008;1(3):496–509. doi: 10.1093/mp/ssn025 19825556

19. Shigenaga AM, Argueso CT, editors. No hormone to rule them all: Interactions of plant hormones during the responses of plants to pathogens. Semin Cell Dev Biol; 2016: Elsevier.

20. Wang D, Pajerowska-Mukhtar K, Culler AH, Dong X. Salicylic acid inhibits pathogen growth in plants through repression of the auxin signaling pathway. Curr Biol. 2007;17(20):1784–90. doi: 10.1016/j.cub.2007.09.025 17919906

21. Fu J, Wang S. Insights into auxin signaling in plant-pathogen interactions. Front Plant Sci. 2011;2(November):74. doi: 10.3389/fpls.2011.00074 Fu2011. 22639609

22. Naseem M, Kaltdorf M, Dandekar T. The nexus between growth and defence signalling: Auxin and cytokinin modulate plant immune response pathways. J Exp Bot. 2015;66(16):4885–96. doi: 10.1093/jxb/erv297 Naseem2015. 26109575

23. Berens ML, Berry HM, Mine A, Argueso CT, Tsuda K. Evolution of Hormone Signaling Networks in Plant Defense. Annu Rev Phytopathol. 2017;55(1):annurev—phyto—080516—35544. doi: 10.1146/annurev-phyto-080516-035544 Berens2017. 28645231

24. Shine MB, Yang JW, El-Habbak M, Nagyabhyru P, Fu DQ, Navarre D, et al. Cooperative functioning between phenylalanine ammonia lyase and isochorismate synthase activities contributes to salicylic acid biosynthesis in soybean. New Phytol. 2016;212(3):627–36. doi: 10.1111/nph.14078 Shine2016. 27411159

25. Camejo D, Guzmn-Cedeo n, Moreno A. Reactive oxygen species, essential molecules, during plant-pathogen interactions. Plant Physiol Biochem. 2016;103:10–23. doi: 10.1016/j.plaphy.2016.02.035 Camejo2016. 26950921

26. de Vries S, von Dahlen JK, Uhlmann C, Schnake A, Kloesges T, Rose LE. Signatures of selection and host-adapted gene expression of the Phytophthora infestans RNA silencing suppressor PSR2. Mol Plant Pathol. 2017;18(1):110–24. doi: 10.1111/mpp.12465 DeVries2017. 27503598

27. Derevnina L, Dagdas YFa. Tansley insight Nine things to know about elicitins. New Phytol. 2016;212(4):888–95. doi: 10.1111/nph.14137 Derevnina2016. 27582271

28. Yang F, Li W, Jrgensen HJL. Transcriptional reprogramming of wheat and the hemibiotrophic pathogen Septoria tritici during two phases of the compatible interaction. PLoS ONE. 2013;8(11):1–15. doi: 10.1371/journal.pone.0081606 Yang2013. 24303057

29. Asai S, Shirasu K. Plant cells under siege: Plant immune system versus pathogen effectors. Curr Opin Plant Biol. 2015;28:1–8. doi: 10.1016/j.pbi.2015.08.008 Asai2015. 26343014

30. Li B, Meng X, Shan L, He P. Transcriptional Regulation of Pattern-Triggered Immunity in Plants. Cell Host and Microbe. 2016;19(5):641–50. doi: 10.1016/j.chom.2016.04.011 Li2016. 27173932

31. Bigeard J, Colcombet J, Hirt H. Signaling mechanisms in pattern-triggered immunity (PTI). Molecular Plant. 2015;8(4):521–39. doi: 10.1016/j.molp.2014.12.022 Bigeard2015. 25744358

32. Macho AP, Zipfel C. Plant PRRs and the activation of innate immune signaling. Mol Cell. 2014;54(2):263–72. doi: 10.1016/j.molcel.2014.03.028 Macho2014. 24766890

33. Birkenbihl RP, Liu S, Somssich IE. Transcriptional events defining plant immune responses. Curr Opin Plant Biol. 2017;38:1–9. doi: 10.1016/j.pbi.2017.04.004 Birkenbihl2017. 28458046


Článok vyšiel v časopise

PLOS One


2019 Číslo 9
Najčítanejšie tento týždeň
Najčítanejšie v tomto čísle
Kurzy

Zvýšte si kvalifikáciu online z pohodlia domova

Získaná hemofilie - Povědomí o nemoci a její diagnostika
nový kurz

Eozinofilní granulomatóza s polyangiitidou
Autori: doc. MUDr. Martina Doubková, Ph.D.

Všetky kurzy
Prihlásenie
Zabudnuté heslo

Zadajte e-mailovú adresu, s ktorou ste vytvárali účet. Budú Vám na ňu zasielané informácie k nastaveniu nového hesla.

Prihlásenie

Nemáte účet?  Registrujte sa

#ADS_BOTTOM_SCRIPTS#