Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 RESEARCH ARTICLE Open Access EST sequencing and gene expression profiling of defence-related genes from Persea americana infected with Phytophthora cinnamomi Waheed Mahomed1,2 and Noëlani van den Berg1,2* Abstract Background: Avocado (Persea americana) belongs to the Lauraceae family and is an important commercial fruit crop in over 50 countries The most serious pathogen affecting avocado production is Phytophthora cinnamomi which causes Phytophthora root rot (PRR) Root pathogens such as P cinnamomi and their interactions with hosts are poorly understood and despite the importance of both the avocado crop and the effect Phytophthora has on its cultivation, there is a lack of molecular knowledge underpinning our understanding of defence strategies against the pathogen In order to initiate a better understanding of host-specific defence we have generated EST data using 454 pyrosequencing and profiled nine defence-related genes from Pc-infected avocado roots Results: 2.0 Mb of data was generated consisting of ~10,000 reads on a single lane of the GS FLX platform Using the Newbler assembler 371 contigs were assembled, of which 367 are novel for Persea americana Genes were classified according to Gene Ontology terms In addition to identifying root-specific ESTs we were also able to identify and quantify the expression of nine defence-related genes that were differentially regulated in response to P cinnamomi Genes such as metallothionein, thaumatin and the pathogenesis related PsemI, mlo and profilin were found to be differentially regulated Conclusions: This is the first study in elucidating the avocado root transcriptome as well as identifying defence responses of avocado roots to the root pathogen P cinnamomi Our data is currently the only EST data that has been generated for avocado rootstocks, and the ESTs identified in this study have already been useful in identifying defence-related genes as well as providing gene information for other studies looking at processes such as ROS regulation as well as hypoxia in avocado roots Our EST data will aid in the elucidation of the avocado transcriptome and identification of markers for improved rootstock breeding and screening The characterization of the avocado transcriptome will furthermore form a basis for functional genomics of basal angiosperms Background Avocado (Persea americana Mill.) is an important agricultural crop in over 50 countries worldwide and is native to Mexico and Central America [1] It belongs to the genus-Persea, subgenus-Persea, family-Lauraceae and falls under the clade of magnoliids that are sister to eudicot and monodicot clades P americana is a diploid angiosperm consisting of 24 chromosomes with approximately 8.83 × 108 (883 Mb) base pairs (bp) To date, the avocado genome is not yet available and only a limited * Correspondence: noelani.vdberg@fabi.up.ac.za Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Pretoria, 0002, South Africa Full list of author information is available at the end of the article number (16558) of expressed sequence tags (ESTs) generated from only fruit and flowers have been sequenced, annotated and released on the NCBI database Phytophthora root rot (PRR), caused by Phytophthora cinnamomi Rands, is considered the most destructive pathogen-induced disease to the avocado industry [2-4] with production relying heavily on the use of phosphite trunk injections and tolerant rootstocks such as Dusa® [4,5] supported by planting in high organic matter soils and mulching to promote antagonistic microbial growth against P cinnamomi Metalaxyl has also showed promising results when used in conjunction with the tolerant rootstock Duke in California in the 1980s [6] However in South Africa it was reported that after © 2011 Mahomed and van den Berg; licensee BioMed Central Ltd This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 successful application for a period of two years, metalaxyl application became inefficient in controlling the disease [7] Recently P cinnamomi has shown a decrease in sensitivity to phosphite treatments after prolonged usage [8] The authors demonstrated that P cinnamomi isolates exposed to long periods of phosphite treatment in south west Australia showed reduced sensitivity to the fungicide when evaluated on avocado, lupins and eucalyptus The population of P cinnamomi isolated from phosphite treated sites colonized phosphite treated plant material easier than isolates not previously exposed to the fungicide This decreased sensitivity to phosphite could indicate the onset of resistance to the fungicide As early as 1926 avocado researchers identified that the success of the avocado industry lay in rootstock improvement [9] The world’s largest rootstock germplasm is maintained in California since 1957, with the hope of identifying more tolerant rootstocks for cultivation [6,10] To date a small number of rootstocks have been identified with partial resistance to P cinnamomi such as Thomas, Martin Grande, Barr Duke, Duke and D9 [11] In South Africa, the devastation caused by P cinnamomi in the 1970s prompted the importation of clonal rootstocks and the development of a large scale selection program in the 1980s For many years Duke remained the industry standard in South Africa, until 2002 with the release of the Dusa® rootstock by Westfalia Technological Services Dusa® gave avocado farmers a better alternative to Duke that showed improved tolerance to P cinnamomi as well as good fruit yields The avocado breeding program at Westfalia is a continuing process and uses previously identified tolerant rootstocks as parents to undergo open pollination Recently, field trials were conducted on a selection of rootstocks in Queensland, Australia with some selections such as ‘SHSR-02’, ‘SHSR-04’, un-grafted ‘Hass’ and Dusa ® demonstrating their tolerance to PRR [12] Despite the importance of avocado and a 60 year attempt to unravel the host pathogen interaction, our knowledge is based on; the analysis of root exudates[13], chemical analysis of roots [12], the application of chemicals to aid in suppression of the pathogen [14], and biochemical studies [15] Histological studies on roots infected with P cinnamomi have aimed to try and understand the plant pathogen interaction [16] It was observed that necrophylactic periderm and periclinal cell wall division occurred, which limited the pathogens progress but did not affect the viability of the pathogen or reduce its ability to infect the host plant P cinnamomi infect the plants roots via motile zoospores present in the soil The attraction of zoospores was investigated by Botha and Kotze in 1989 and it was found to be influenced by the composition of 14 amino Page of 14 acids in avocado root exudates [13] Sánchez-Pérez and colleagues tested crude root exudates for P cinnamomi mycelial inhibition and subsequently the compound known as stigmastan-3, 5-diene was identified as the inhibitory compound [12] García-Pineda et al (2010) investigated reactive oxygen species (ROS) formation and the role of nitric oxide (NO) against P cinnamomi [15] The authors observed an increase in ROS and NO levels and deduced that the increase in ROS observed may assist in weakening host tissue early in infection with the sharp increase in NO possibly resulting in salicylic acid (SA) accumulation This accumulation could cause an SA mediated H2O2 burst by the suppression of H2O2 degradation The authors hypothesize that (cytosolic tobacco catalase) CAT is bound by SA, which inhibits CATs H O degrading activity The effect of externally applied SA on root colonisation was also investigated and indicated that decreased root colonisation was associated with SA application SA has been implicated in regulating cell death, inducing resistance responses and activating various defence genes such as pathogenesis-related (PR) genes [17] but the mode of action has not been elucidated The production of NO and ROS have previously been demonstrated to activate cell death These early attempts on investigating the interaction between avocado and P cinnamomi have illustrated the complexity of the defence response, highlighting the need for the molecular elucidation of defence genes Molecular research on avocado has comprised of genetic relationship studies and the molecular characterization of the fruit and flowers There has been some gene characterization of avocado fruit ripening genes [18-22] The greater part of molecular detail exists due to a continuous effort in marker development to assist in either elucidating genetic relationships amongst scions [23-28], or scion improvement [29-32] There is currently a preliminary genetic map available based on microsatellites, random amplified polymorphic DNA (RAPD) markers and DNA fingerprint (DFP) markers [33] The most recent molecular development in the fight against PRR was the identification of 70 microsatellite markers that were developed from over 8000 ESTs in the hope of aiding in marker assisted breeding against PRR The ESTs were however from a floral gene database generated for comparative genomics research of basal angiosperms Their efficacy has yet to be tested for use in identifying tolerant rootstocks, but it is known that they amplify across all avocado varieties and can be used for investigation of genetic relations [34,35] The University of California Riverside (UCR) has recently employed 61 polymorphic AFLP markers to characterise PRR tolerance in 83 rootstocks from various locations including South Africa and Israel with the Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 majority of rootstocks from the UCR collection [36] The study concluded that resistance mechanisms vary between tolerant cultivars and no trend was observed in the cluster analysis The avocado/P cinnamomi interaction has not previously been elucidated on a molecular level Current knowledge is based on research of the non-host plant, Arabidopsis A study conducted on Arabidopsis infected with P cinnamomi revealed that ROS induction, HR activation, lignin synthesis and callose production was initiated upon infection The non-host showed activation of the ethylene and jasmonic acid pathways and only a minor involvement of the SA pathway [37] in contrast to the study conducted by García-Pineda et al (2010) on avocado which indicated that SA is a major inhibitor of pathogen colonisation Macroscopic changes such as callose production have also been observed during P cinnamomi infection in maize [38] and although model plants like Arabidopsis provide an insight to defence responses there are differences between non-host and host-specific defence responses In order to fully understand tolerance in avocado it is important to conduct molecular level studies on the host specific interaction between P americana and P cinnamomi Since there is no genome data available for avocado, the identification and characterisation of genes is difficult EST generation supplements the lack of genome data by providing transcript specific information and excluding non-coding regions of the genome Highthroughput sequencing is well suited for large scale EST discovery, providing a tool for gene discovery in nonmodel crops to evaluate the changes in gene expression to abiotic or biotic stresses [39-41] The cost of pyrosequencing is also lower than conventional EST sequencing, small transcripts are not lost, and the time in sequence generation from tissue isolation is greatly shortened [39] More specifically, it is advantageous for commercial crops that lack substantial molecular databases and will aid in their unconventional improvement [40] Avocado is one such commercial crop that is in need of development of molecular tools for the improvement of the crop Avocados’ importance as an agricultural crop has justified molecular investigation and the application of modern molecular tools for its improvement [26,27,31,33,42,43] The application of high-throughput sequencing to avocado is the next step in improving breeding of this economically important crop In this study ESTs of a tolerant avocado rootstock infected with Phytophthora cinnamomi were generated The 454 GS-FLX platform was used to generate sequence data for several time points including the uninfected, 6, 12, 24, 48 and 72 hours after infection, as well as to identify transcripts that were associated with Page of 14 the defence response We identified 371 transcripts from avocado and studied the gene expression of a selection of these ESTs, thereby providing the first molecular data for the avocado/P cinnamomi interaction Results 454 pyrosequencing and assembly Three cDNA libraries-uninfected (0 hr), library (6 & 12 hr) and library (24, 48 & 72 hr) were sequenced on a single lane on the GS FLX platform and generated a total of Mb of data (after trimming and quality control) consisting of 9953 reads and resulted in the assembly of 371 contigs (Table 1) These contigs comprised of 1407 reads from the uninfected library, 3584 reads from infection library and 4962 reads from infection library The average read lengths for the libraries were 216.4 bp for uninfected, 217.5 bp for library and 215.9 bp for the library (Figure 1) The pyrosequencing run was efficient based on the maximal amount of data obtainable being 2.5 Mb with a maximal read length of 250 bp EST identification and classifications After analysis using the dCAS software, 367 novel ESTs were identified for P americana The program used BLASTX amino acid comparisons to screen for homology of the contigs against the NCBI non-redundant (NR) database A large proportion of the sequences generated showed homology to hypothetical proteins and 45 of the 371 contigs had no similarity to previously annotated sequences (Table 2) Of the 371 contigs identified, only two sequences showed homology to previously identified fructose-bisphosphate aldolase and metallothionein type-II proteins from avocado, with the remaining 369 not having any avocado sequence homolog in the NR database Manual BLAST annotation did not influence transcript identification Contigs were grouped into functional classes according to the GO (Gene Ontology) and KOG (Eukaryotic Orthologous Groups) databases Nine percent of contigs were grouped into the unknown functional class in the KOG database while 44.5% of contigs from the GO classification were represented by unknown functions Table Excerpts of newblermetric reports from the uninfected, library and library libraries of Phytophthora cinnamomi infected avocado roots Uninfected response Library Library Total Reads 1407 3582 4961 Total Bases 288885 737254 1017064 Number of Contigs 43 139 189 The total number of reads and contigs are shown to illustrate the efficiency of the pyrosequencing run Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 Page of 14 Figure Read length distributions of uninfected, library and library infection Dusa® cDNA Pyrosequencing was performed on the 454 GS-FLX platform (a) Uninfected library contains reads with the highest frequency at around 245 bp (b) Library reads have the highest frequency at around 252 bp (c) Library reads have the highest frequency at around 240 bp (Figure & Table 3) The categories of post-translational modification; translation, ribosomal structure and biogenesis; signal transduction mechanisms and general function prediction contained a combined total of 34.8% of all contigs According to the GO database the Table Contig classification for cDNA libraries of Phytophthora cinnamomi infected avocado roots Uninfected Library Library Total Unidentified 16 24 45 Hypothetical protein 23 54 75 152 Genes identified 15 69 89 functional classifications of cell wall related; protein binding; stress response; ribosomal structural constituent; cytoplasmic biological processes; cellular component and other categories comprised 40.8% of all contigs with 3% (12/371) of the contigs linking directly with stress responses Over 20 putative defence related genes were identified ranging from general defence-related genes (metallothioneins, thaumatin and universal stress genes) to more specific oomycete defence-related genes (pathogenesis related protein PR10 and the oxysterolbinding gene) (Tables 4&5) 173 Contigs were classed as unidentified, identified or hypothetically identified if the sequence homology search revealed that there was no similarity, significant similarity or inferred structural function respectively The total number of genes identified was 173 with only 45 of the total 370 contigs generating no identification using the non-redundant (NR) database on the NCBI Species similarity between avocado and other plants We observed significant sequence homology between Vitis vinifera (grape) and avocado when the species origin of the sequence similarity was investigated The top three represented species according to amino acid Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 Page of 14 Figure KOG (euKaryotic Orthologous Groups) classifications of avocado transcripts identified in three cDNA libraries The contigs generated from the 454 data were compared against the KOG database to assign functional classifications homology on the NCBI were V vinifera, Arabidopsis thaliana and Oryza sativa, with V vinifera having the majority of the hits in all three libraries Twenty two percent of sequences showed homology to V vinifera sequences with 7.5% belonging to A thaliana and 7.8% of sequences to O sativa Homology to P americana was found in only 1% of sequences (4/371) (Figure 3) Only two genes were represented by the 1% in which Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 Page of 14 Table Contigs of avocado transcripts grouped into functional classes according to GO database Gene ontology Number of contigs Unknown 165 Other 48 ATP binding Biological process - cytoplasm 20 Cellular component 39 Response to stress 12 RNA binding Cell wall related Structural constituent of ribosome - translation - ribosome 10 12 Transcription factor activity - regulation of transcription Transferase activity/cell wall biogenesis Translation elongation factor activity - translation factor activity, nucleic acid binding Transporter activity - transport Water channel activity - transport - membrane Protein binding 10 Mitochondrion Kinase activity Protein folding - cellular component Membrane Lipid binding - lipid transport A large proportion of contigs (44.5%) fell into the category of unknown classification while contigs that link directly with stress responses constituted 3% of the total number of contigs the metallothionein transcript featured three times and fructose-bisphosphate aldolase featured once Grape vine featured among the top ten homologous hits of every contig that was annotated Thirty seven percent of the annotated contigs were represented by various plant species such as Prunus armeniaca, Solanum tuberosum, Hevea brasiliensis with the variety of plant species not biased to any particular family or order The majority of the species similarities relate to a large variety of plants that have been collectively categorised as other Quantitative gene expression analysis Expression analyses of nine genes were conducted at 0, 3, 6, 12, 24 and 48 hours post infection (hpi) to validate if the pyrosequencing data reflected their gene expression This was normalized against 18S and actin reference genes to give the relative gene expression The expression data was then compared against the pyrosequencing data which revealed that six of the nine genes showed similarities between the two methods, showing the highest expression at a time point belonging to the library from which the transcript emanated (Table 6) Thaumatin expression was significantly greater at 48 hpi (1.1) as oppose to the uninfected (0.4), as well as the & hpi (Figure 4a) The expression pattern indicated that thaumatin was only regulated in response to P cinnamomi by 12 hpi and increased by nearly threefold over a 36 hour period Thaumatin levels were significantly higher in the later infection time points when compared to the earlier time points-thus correlating with the pyrosequencing data The pathogenesis-related (PR-10) psemI gene showed significant increases in expression at & 24 hpi At 24 hpi psemI reached the highest expression level of 1.5 when compared to all time points (Figure 4b) By 48 hpi psemI expression had decreased significantly to 0.1, reaching levels comparable to hpi Cytochrome P450-like TBP (TATA box binding protein) showed a significant early response at hr after infection with P cinnamomi reacting with an increase of ten-fold At hpi the gene was significantly down-regulated followed by a substantial increase at 12 hpi- a similar level found at hpi The gene was then significantly down-regulated to 1.0 at 24 hpi and remained unchanged at 48 hpi (Figure 4c) Cytochrome P450-like TBP levels were constantly up- and down-regulated showing significant variation over time points The data was consistent with the pyrosequencing data for this transcript The gene encoding for a metallothionein-like protein was constitutively expressed at 0.6 and showed no significant changes in expression over the first hpi This was however followed by a significant increase in the expression at 12 hpi when compared to all other time points reaching levels of 3.2, the expression then decreased to 0.5 at 24 hpi and remained unchanged at Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 Page of 14 Table List of putative stress related genes isolated from Phytohpthora cinnamomi infected avocado roots Contig Putative identity E Value Species contig0020 leucine-rich repeat resistance protein-like protein 8e-15 Gossypium hirsutum contig0070 contig0088 4-coumarate-CoA ligase-like protein seven transmembrane protein Mlo 5e-26 1e-28 Arabidopsis thaliana Zea mays contig0106 pathogenesis-related protein PsemI 5e-14 Pseudotsuga menziesii contig0109 drought-induced protein 8e-13 Retama raetam contig0076 metallothionein-like protein type 6e-41 Persea americana contig00007 translationally controlled tumour protein like protein 2e-07 Nicotiana tabacum contig00011 thaumatin 2e-20 Vitis riparia contig00043 cinnamate-4-hydroxylase 3e-16 Gossypium arboreum contig00064 contig00065 metallothionein-like protein type AP2 domain containing protein 7e-41 2e-41 Persea americana Prunus armeniaca contig00073 oxysterol-binding protein 3e-28 Solanum tuberosum thaumatin-like protein, putative 3e-29 Arabidopsis thaliana contig00095 AP2 domain containing protein 7e-27 Prunus armeniaca contig00108 profilin-like protein 5e-17 Cinnamomum camphora contig00163 Translationally-controlled tumour protein homolog (TCTP) translationally controlled tumour protein 1e-29 Hevea brasiliensis contig00169 cysteine proteinase 8e-39 Elaeis guineensis contig00175 putative universal stress protein 2e-40 Cicer arietinum contig00187 contig00054 cytochrome P450 like TBP metallothionein-like protein class II 2e-60 4e-39 Nicotiana tabacum Nelumbo nucifera contig00057 dormancy/auxin associated family protein 6e-15 Arabidopsis thaliana contig00081 putative aquaporin PIP2-1 5e-76 Vitis berlandieri × Vitis rupestris library library library library Six genes were isolated from library cDNA while 15 genes were identified from the library cDNA These transcripts play a role in either biotic or abiotic stress responses showing the species to which the sequence showed homology 48 hpi (Figure 4d) The data showed similarity to the pyrosequencing data for this transcript The profilin-like gene was expressed constitutively at 1.9 prior to infection Three hours after infection the transcript was significantly down-regulated to 0.6 (a fold decrease) and remained unchanged at hpi There was a significant up-regulation from hpi to12 hpi with expression peaking at 2.6, followed by a significant decrease to 0.9 at 24 hpi and remained unchanged at 48 hpi as opposed to 12 hpi (Figure 4e) The MLO transmembrane protein encoding gene was constitutively expressed at 2.2 followed by a significant reduction at hpi compared to the uninfected time point At 12 hpi the tolerant rootstock responded with a significant increase to 5.7 Expression at 24 hpi was significantly down-regulated when compared to 12 hpi and Table Defence-related genes isolated from Phytophthora cinnamomi infected avocado roots Gene E-value -20 Response against thaumatin 2e Vitis riparia Armellaria mellea [65] metallothionein-like protein type 7e-41 Persea americana Agrobacterium rhizogenes [53] pathogenesis-related protein P sem I 5e-14 Pseudotsuga menziesii Phellinus weirii [49] putative universal stress protein profilin-like protein 2e-40 5e-17 Cicer arietinum Cinnamomum camphora Phytophthora infestans [60] oxysterol-binding protein 3e-28 Solanum tuberosum Phytophthora spp [66] LRR resistance protein-like protein 8e-15 Gossypium hirsutum Phytophthora infestans [64] seven transmembrane protein Mlo 1e-28 Zea mays Blumeria graminis f sp hordei [62] Certain transcripts could be related directly to oomycete infection in other plants Strong E-values indicate confidence in the identification of these defencerelated genes Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 Page of 14 Figure Number of contigs grouped according to sequence homology between avocado and other plant species The sequence similarities were analysed to establish which species was most represented by the 454 data There is an observable lack of avocado sequence data available on public databases reached a level of 1.5 at 48 hpi (Figure 4f) The mlo expression data was in agreement with the pyrosequencing data for this transcript The universal stress protein showed the maximum expression at 12 hpi but did not show a significant increase when compared to and hpi An overall increase in expression was viewed until 12 hpi, followed by a significant down-regulation to 0.4 at 48 hpi compared to 12 hpi (Figure 4g) Two genes encoding respectively for the thaumatinlike protein and the leucine rich repeat (LRR) resistance protein-like protein were constitutively expressed at 0.8 Both genes showed no statistically significant change in regulation over the 48 hour time course, however similar to the majority of genes in this experiment, both genes showed the highest increase in expression at 12 hpi The highest level of up-regulation achieved was 1.1 and 1.47 respectively (Figure 4h &4i) Table Similarities between pyrosequencing data and gene expression profiles of defence-related genes Sequence ID GenBank accession number cDNA Library Max qRT-PCR expression Similarities between 454 and qRTPCR Thaumatin LRR resistance PLP JO840464 JO840460 Library Library 48 hpi 12 hpi yes yes Metallothionein like protein JO840461 Uninfected, Library 12 hpi yes Thaumatin-like protein JO840465 Library 12 hpi yes Seven transmembrane protein MI0 JO840462 Library 12 hpi yes Pathogenesis-related protein PsemI JO840463 Library 24 hpi no Proflin-like protein JO840466 Library 12 hpi no Putative universal stress protein JO840467 Library 12 hpi no Cytochrome P450 like TBP JO840468 Library and Library 3/12 hpi yes All the genes chosen for expression profiling from the tolerant avocado rootstock infected with Phytophthora cinnamomi showed the highest expression from a time point related to the cDNA library of their identification except the pathogenesis-related protein, profilin-like protein and the universal stress protein Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 Page of 14 Figure Gene expression of Dusa® -a tolerant avocado rootstock, infected with Phytophthora cinnamomi Expression analysis was conducted at 0, 3, 6, 12, 24 & 48 hpi (hours post infection) with hr being the uninfected control The data was normalized using two reference genes-actin and 18S Expression analysis was performed in triplicate on three biological replicates (a) Thaumatin (b) Pathogenesisrelated protein psemI (PR10) (c) Cytochrome P450 (d) Metallothionein-like gene (e) Profilin-like gene (f) MLO transmembrane gene (g) The universal stress protein (h) The thaumatin-like gene (i) The LRR resistance protein-like protein Discussion We have sequenced the first set of avocado root transcriptomic data for the avocado-P cinnamomi interaction A single lane of pyrosequencing on the GS FLX platform generated 2.0 Mb (of a potential 2.5 Mb) of data, consisting of 9953 reads that assembled into 371 contigs of which 367 ESTs are novel for P americana and have not previously been identified In addition, we were also able to identify and quantify the expression of nine defence-related genes that were regulated in response to P cinnamomi The primary objective of this study was to generate EST data of a tolerant avocado rootstock infected with the root pathogen P cinnamomi This data identified the genes involved in cellular processes and defence mechanisms thereby providing the first platform for studying molecular mechanisms underlying tolerance in the roots of one of the important agricultural hosts of P cinnamomi The 371 contigs were grouped into 38 and 21 functional classes based on the KOG and GO databases respectively using the dCAS program As expected, the majority of sequences had unknown functions Due to the high sensitivity of sequencing, transcriptome studies identify many transcripts that have not yet been characterised and many that have unknown functions even when annotated using a database such as Gene ontology [44] The lack of EST and genome sequence data for avocado in general, specifically rootstocks, also accounts for the high frequency of unknown functions The top ten functional groupings according to the GO classification revealed that 44.5% of assembled contigs were represented by unknown functions followed by the functional groups of ‘other’, ‘cellular components’, ‘biological processes’, ‘stress responses’, ‘ribosome structure’, ‘cell wall related’, ‘protein binding’, ‘mitochondrion and ‘ATP-binding’ According to the KOG database much of sequence data matched categories of ‘general function prediction’, which means that these transcripts show homology to transcripts that are poorly characterised according to the NCBI The KOG database revealed that the top ten classes that the contigs grouped into were firstly; ‘general function prediction’ followed by ‘signal Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 transduction’, ‘unknown function’, ‘translation and ribosomal structure’, ‘chaperones’, ‘carbohydrate metabolism’, ‘intracellular trafficking’, ‘transcription’, ‘cytoskeleton’ and ‘inorganic ion transport and metabolism’ Furthermore, the presence of unidentified reads in this study is not unique to avocado, other studies have also produced sequence that did not align to any sequence present in NCBI datasets [39] We investigated the sequence homology between the avocado sequence data and the plant species that our data showed homology to Only 1% (4/371) of the sequenced contigs showed homology to P americana, while 20-30% of the contigs generated showed similarity to grapevine (V vinifera) The lack of similarities to any avocado sequence data observed in our study emphasizes the lack of genetic data on the NCBI The knowledge we have gained by sequencing avocado rootstock ESTs may provide some insight into other magnoliids or phylogenetically related plants The sequence and expression data generated in this study can form a basis for functional genomics of basal angiosperms - a group which has no other model [27] As expected we isolated a number of ESTs with homology to genes previously associated with defence responses in plants against pathogens and in some cases, against oomycetes Interesting defence ESTs included thaumatin, metallothionein, a PR10 pathogenesis-related protein, a mlo transmembrane protein and profilin Nine genes were quantified with qRT-PCR to elucidate the early gene response of a tolerant avocado rootstock infected with P cinnamomi as well as to validate the pyrosequencing data Thaumatin, a PR5 protein associated with the SA pathway [45,46], was significantly upregulated at 48 h in response to P cinnamomi infection The gene showed no changes in regulation during the first hours after inoculation with a mycelial suspension At 12 and 24 hpi expression showed an insignificant but steady increase in response to the infection PR5 is induced by biotic stress and further linked to increased pathogen resistance [47] García-Pineda et al (2010) showed decreased root colonization in the Arabidopsis-P cinnamomi system linked to SA The significant up-regulation of thaumatin in the P cinammomi tolerant avocado rootstock indicates the importance of the SA pathway in the early inhibition of the hemibiotroph P cinnamomi Hemibiotrophs have an initial biotrophic phase prior to becoming necrotrophs and PR5 gene activity in the SAdependant pathway has been previosuly shown to be effective against biotrophs [48] PsemI was highly expressed at 24 hpi in tolerant avocado roots infected with P cinnamomi The PR10 gene was identified in the response of Douglas -fir infected with Phellinus weirii [49] The authors showed that very Page 10 of 14 high concentrations of Pin m III (PsemI gene homologue) was responsible for resistance to the rust pathogen Cronartium ribicola This PR-10 gene has also been used as a marker in screening for P weirii resistance in Douglas-fir and could therefore be valuable in screening for PRR tolerance in avocado The gene encoding for the cytochrome P450-like TBP was the only transcript to be significantly induced by P cinnamomi as early as hpi This enzyme features in oxidative metabolism and the production of ROS This rapid response could be attributed to the universal nature of the protein in cell metabolism and growth Additionally it has been reported to be involved in biotic and abiotic environmental responses as well as in the HR response to infection [44,50-52] Our data revealed a noticeable host response at 12 hpi with the up-regulation of four transcripts, the metallothionein-like gene, the universal stress protein, profilin &mlo Metallothioneins inhibit programmed cell death (PCD) and Fumonisin B1-induced root death in tomato infected with Agrobacterium rhizogenes through interference of the ROS pathway ROS accumulation was significantly reduced under metallothionein over-expression, validating its function in ROS scavenging [53] The significant induction of metallothionein in the highly tolerant avocado rootstock at 12 hpi implies that this protein may play a role in conferring disease tolerance to P cinnamomi by scavenging ROS ROS generation is indicative of the activity of the hypersensitive response (HR), which leads to cell death and is effective against biotrophic organisms [54] The up-regulation of the universal stress protein at 12 hpi indicates the plant’s response to the stress of infection by P cinnamomi Universal stress proteins are mediated by ethylene [55], and our results may therefore implicate the involvement of the ethylene pathway in response to P cinnamomi, a pathway that has shown activation in the Arabidopsis/P cinnamomi interaction [37] Both profilin &mlo play a role in actin filament polarization [56,57] and actin rearrangement has been observed in plant-fungus interactions with successful pathogen infection resulting in the suppression of the rearrangement [58,59] Profilin is known to localize to the site beneath the cell wall, that is penetrated by the oomycetous appresorium [60] and either promotes or prevents actin polymerization in the actin cytoskeleton [56] Cell wall thickening during fungal attack also involves the re-orientation of actin filaments as a defense response in order to prevent pathogen ingress [61] The up-regulation of profilin in avocado roots suggests that profilin is being produced in response to P cinnamomi penetration MLO, a transmembrane protein, modulates actin cytoskeleton polarization in Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 resistant barley in response to a biotrophic fungus- Blumeria graminis f sp hordei [62] Successful defence against pathogens results in cell wall strengthening and is correlated with increased actin accumulation at sites of attempted penetration [57] The thaumatin-like gene expression showed no significant response to the oomycete infection This might be explained by the fact that it is usually induced by viral, bacterial and fungal infection [46], and is believed to destroy fungal cell walls using a variety of enzymatic activities [63] The LRR (leucine rich repeat) resistance protein-like gene also demonstrated no significant response to P cinnamomi Although resistance related LRR proteins have been found to interact specifically with other Phytophthora species, (Solanum tuberosumPhytophthora infestans interaction) [64], no such interaction has to our knowledge been described for P cinnamomi and its hosts, specifically avocado However, the identification of a LRR-like gene from P cinnamomi infected tolerant avocado roots warrants further investigation Conclusions This study identified an interesting set of genes regulated by the infection of the hemibiotroph P cinnamomi by using 454 pyrosequencing Defence genes included general defence-related transcripts (universal stress protein, metallothionein and thaumatin-like genes as well as cytochrome P450) and Phytophthora specific response genes such as (thaumatin, PR-10 and the LRR resistance protein-like gene) [44-47,49-52,55,64] We inoculated avocado roots with a mycelial suspension, based on the genetic response by the plant, we hypothesize that the plant response is delayed due to the slower infection rate of mycelia as opposed to zoospores that are able to germinate and encyst within hr of release from sporangia Despite this delay five of the nine defence-related transcripts showed a significant early response to the pathogen between and 12 hpi Based on this first set of transcriptome data we hypothesize that the tolerance of the rootstock in this study is most likely polygenic and based on the early detection of P cinnamomi followed by a response that included ROS and cell-wall strengthening Ongoing research in our laboratory has generated a second set of transcriptome data and has included a variety of rootstocks with different levels of PRR tolerance and susceptibility as well as additional time-points We have successfully produced the first molecular data for the avocado-Phytophthora cinnamomi interaction and believe that this data will contribute to the understanding of host defence against this devastating pathogen thereby aiding in the selection of tolerant avocado rootstocks Page 11 of 14 Methods Plant material inoculation Nine month old tolerant Dusa® clonal avocado plantlets were provided by Westfalia Technological Services (Tzaneen, South Africa) and inoculated with a Phytophthora cinnamomi mycelial suspension containing mycelia and zoospores A total of 33 g of mycelia was homogenised in 65 L of distilled water giving a final concentration of 0.5 g/L This was then mixed into 112 kg of vermiculite in a mistbed Plantlets were randomly grounded in vermiculite and constantly irrigated over a period of six weeks Scanning electron microscopy and confocal microscopy were used to confirm P cinnamomi infection, germination and root colonization which was then used to determine time points for root harvesting Based on the results obtained, zoospores germinated within hr & hyphae were visible at hr, we selected our time points accordingly Root material was harvested at hour (uninfected), 3, 6, 12, 24, 48 and 72 hours post infection (hpi), snap frozen in liquid nitrogen and stored on dry ice (-78°C) until the root material could be transported back to a -80°C freezer A subset of plants was left in the mistbed for six weeks as a positive control for disease Generation of cDNA libraries for pyrosequencing RNA isolations were done using the CTAB method [24] Roots were ground in liquid nitrogen and 2-3 g of root material was used per RNA extraction The Chloroform: isoamyl alcohol wash step was repeated times followed by washing with ethanol thrice Total RNA concentrations were quantified using the Nanodrop ND-100 Spectrophotometer (Nanodrop Technologies, Inc., Montchanin, USA) and verified on a 2% non-denaturing TAE agrose gel Total RNA from three biological replicates (per time point) was combined before mRNA purification Three technical replicates per biological replicate were performed for RNA isolation Prior to mRNA isolation, different time points were combined into three libraries for pyrosequencing and designated as uninfected, library and library The hr time point and was regarded as the uninfected library with library containing the & 12 hr infection time points while library contained the 24, 48 and 72 h time points Purification of mRNA was done according to manufacturer’s instructions using Oligotex (Oligotex™ mRNA kit, Qiagen, Valencia, California, USA) The mRNA purification was performed twice per sample to ensure the removal of any DNA contamination as well as to reduce the amount of rRNA available in the sample DNA contamination was assessed by using intron flanking primers F3H forward: (5’-TCTGATTTCGGAGATGACTCGC-3’) and F3H reverse: Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 (5’-TGTAGACTTGGGCCACCTCTTT-3’) to amplify a 300 bp fragment of the flavanone-3-hydroxylase gene from RNA as opposed to the 1200 bp fragment which is obtained from DNA cDNA libraries were synthesized with the Roche cDNA synthesis system (Roche Diagnostics, Mannheim, Germany) according to manufacturer’s instructions and purified using the Qiagen MinElute PCR Purification Kit (Qiagen, Valencia, California, USA) before being sequenced DNA contamination was assessed in cDNA as previously mentioned Pyrosequencing and Bioinformatics Libraries were sequenced by Inqaba Biotec (Sunnyside, South Africa) on the GS-FLX platform Approximately μg of cDNA was supplied for each library A third of each library was tagged with a different ten nucleotide tag (Uninfected tag- 5’CGTGTCTCTA’3, library infection tag- 5’CTCGCGTGTC’3, library infection tag5’TAGTATCAGC’3) and sequenced on a single lane A total of 9950 reads were assembled into 371 contigs using the Newbler assembler version 1.1.02.15 (Roche) Reads were trimmed before contig assembly Low quality reads were not included and the assembly was analyzed and annotated using dCAS (Desktop cDNA Annotation System) Version 1.4.1 Build 3791 and CLC Free Workbench software (CLC bio, Cambridge, MA) The BLASTX tool was used (using the PAM 30 matrix) in order to produce short and nearly exact matches The sequence data generated in this study is available on the NCBI Transcriptome Shotgun Assembly Sequence Database BioProjectID: PRJNA72155 Gene sequences for genes quantified are available through accession numbers TSA JO840460-JO840468 Quantitative gene expression analysis Nine genes were selected for gene expression analysis and included: Thaumatin, thaumatin-like, metallothionein- Page 12 of 14 like, leucine rich repeat resistance protein-like, pathogenesis-related protein PsemI, putative universal stress, profilin-like, transmembrane protein MLO and cytochrome P450-like TBP (TATA box binding protein) Reference genes chosen for the study were actin and 18s rRNA genes The time points used to analyze the gene expression were hours prior infection and, 3, 6, 12, 24 and 48 hours post infection cDNA synthesis was performed for qRT-PCR with starting material of 1-2 μg total RNA of each time point The ImProm-II™ single strand cDNA system from Promega (Promega Corporation, Madison, Wisconsin, USA) was used in conjunction with random hexamer primers from Invitrogen (Invitrogen Life Technologies, Mississauga, ON, Canada) Reverse transcription was carried out under the following conditions: 25°C for min, 42°C for 60 and 70°C for 10 Quantitative PCR Primers were designed using Primer Designer for Windows, version 4.2©, (scientific and educational software, Durham, NC) based on sequence homology Primers were designed to amplify products of no more than 150 base pairs, the selection was then auto adjusted using more aggressive criteria and primers were chosen within a 58-61°C range Primers were synthesised by Southern Cross Biotechnology (Cape Town, South Africa) (Table 1) Quantitative PCR was carried out using the Bio-Rad CFX96 real-time PCR detection system (Bio-Rad Laboratories, Hercules, CA) Reactions were performed in a 20 μl tube containing μl of each cDNA sample, 10 μl of iQ SYBR Green supermix (Bio-Rad), μl of each set of fragment specific forward and reverse primers and μl SABAX water Thermocycling was carried out at 95°C for 10 min, followed by 55 cycles of 95°C for 10 sec an fragment specific annealing temperatures (Table 7) for 15 sec and Table Primer sequences of selected putative avocado defence-related genes from Phytophthora cinnamomi infected avocado roots Sequence ID Fwd primer (5’-3’) Rev primer (5’-3’) Thaumatin CACCCTGTAGTTCACTCC CCAGATGCTTACAGTTACC Product length (bp) 75 Annealing °C 58.5 LRR resistance PLP GACATTCTTATAGCCATC ATAAACAATCTGATTTTG 135 56 Metallothionein like protein type AGTCTTCATCCCTAATACATATCCC GTTTGTGCGTGTCTGGTTTC 76 58.5 Thaumatin-like protein AAGCAGTCCTCAAGGTTC TTTCCGTTAGTGTCAAAGC 79 65 MLO protein TCGTGGATGGAAGGAGTG ATGGGCAAATCTAAATCTTGTTG 85 58.5 Pathogenesis-related protein PsemI Proflin-like protein GAAGATGGAGTACAAATAC TTCGGTATCTATGATGAG CACCTTGATGTGATAAAC ACGATATGACATTCAATAG 82 110 58.5 58.5 Putative universal stress protein GACATTCTTATAGCCATC ATAAACAATCTGATTTTG 135 56 Cytochrome P450 like TBP GTCAAAGTGAAGAAATTC AATCTCGTTAATCCATTC 119 58.5 114 58.5 18S GTCAAAGTGAAGAAATTC AATCTCGTTAATCCATTC Actin GAATCTGGACCATCTATTG TACCAACCAAACCAAATC 59 Sequences were chosen from the library and library based on differential expressions Two reference genes-actin and 18s were included Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 elongation of 72°C for 15 sec Three biological replicates were used with three technical replicates of each The specificity of each primer pair was thoroughly investigated by standard PCR before the quantitative PCR was conducted and then verified by the presence of a single melting temperature peak For the qPCR, the cycle threshold (CT) values were automatically determined using the accompanying Bio-Rad CFX manager (BioRad CFX Manager™ version 1.5) Page 13 of 14 13 14 15 16 Statistical analysis A Student T-test was carried out to determine significant differences between gene expression levels for quantitative gene expression analysis Statistical analysis was performed using the JMP ® program version 9.0.0 (SAS Institute, Inc., Cary, NC) with a 95% confidence interval 17 18 19 20 Acknowledgements We are grateful for the support of Westfalia Technological Services (WTS) for providing the plantlets This work was financially supported by the National Research Foundation (NRF) and Hans Merensky Foundation 21 22 Author details Forestry and Agricultural Biotechnology Institute (FABI), University of Pretoria, Pretoria, 0002, South Africa 2Department of Genetics, University of Pretoria, Pretoria, 0002, South Africa 23 24 Authors’ contributions WM and NVdB designed the experiments with the work being carried out by WM Both authors prepared and approved the final manuscript 25 Received: 22 August 2011 Accepted: 23 November 2011 Published: 23 November 2011 26 References FAOSTAT [http://faostat.fao.org/] Zentmyer GA: Diseases of the Avocado California Avocado Society Yearbook 1955, 39:44-58 Zentmyer G: Avocado diseases Tropical Pest Managment 1984, 30:388-400 Coffey DM: Phytophthora root rot of avocado-an integrated approach to control in California California Avocado Society Yearbook 1987, 71:121-137 Giblin F, Pegg K, Willingham S, Anderson J, Coates L, Cooke T, Dean J, Smith L: Phytophthora Revisited New Zealand and Australia Avocado Grower’s Conference Tauranga; 2005 Coffey DM, Guillemet FB: Avocado rootstocks California Avocado Society Yearbook 1987, 71:173-179 Darvas JM, Bezuidenhout JJ: Control of Phytophthora root rot of avocados by trunk injection South African Avocado Growers’ Association Yearbook 1987, 10:91-93 Dobrowolski MP, Shearer BL, Colquhoun IJ, O’Brien PA, Hardy GES: Selection for decreased sensitivity to phosphite in Phytophthora cinnamomi with prolonged use of fungicide Plant Pathology 2008, 57:928-936 Ben-Ya’acov A, Michelson E: Avocado rootstocks Horticultural Reviews 1995, 17:381-429 10 Zentmyer GA: The search for resistant rootstocks in Latin America California Avocado Society Yearbook 1957, 41:101-106 11 Gabor BK, Guillemet FB, Coffey DM: Comparison of field resistance to Phytophthora cinnamomi in twelve avocado rootstocks Hortscience 1990, 25:1655-1656 12 Sánchez-Pérez J, Jaimes-Lara M, Salgado-Garciglia R, López-Meza J: Root extracts from Mexican avocado Persea americana var drymifolia inhibit 27 28 29 30 31 32 33 34 35 36 the mycelial growth of the oomycete Phytophthora cinnamomi European Journal of Plant Pathology 2009, 124:595-601 Botha T, Kotze JM: Exudates of avocado rootstocks and their possible role in resistance to Phytophthora cinnamomi South African Avocado Growers’ Association Yearbook 1989, 12:64-65 Bekker TF, Kaiser C, Labuschagne N: Efficacy of water soluble silicon against Phytophthora cinnamomi root rot of avocado: A progress report South African Avocado Growers’ Association Yearbook 2006, 29:58-62 García-Pineda E, Benezer-Benezer M, Gutiérrez-Segundo A, RangelSánchez G, Arreola-Cortés A, Castro-Mercado E: Regulation of defence responses in avocado roots infected with Phytophthora cinnamomi (Rands) Plant and Soil 2010, 331:45-56 Phillips D, Grant BR, Weste G: Histological changes in the roots of an avocado cultivar, Duke 7, infected with Phytophthora cinnamomi Phytopathology 1987, 77:691-698 Dempsey DMA, Shah J, Klessig DF: Salicylic acid and disease resistance in plants Critical Reviews in Plant Sciences 1999, 18:547-575 Hammond-Kosack KE, Jones JDG: Plant disease resistance genes Annual Review of Plant Physiology and Plant Molecular Biology 1997, 48:575-607 Christoffersen RE, Tucker ML, Laties GG: Cellulase gene expression in ripening avocado fruit: The accumulation of cellulase mRNA and protein as demonstrated by cDNA hybridization and immunodetection Plant Molecular Biology 1984, 3:385-391 Cass LG, Kirven KA, Christoffersen RE: Isolation and characterization of a cellulase gene family member expressed during avocado fruit ripening Molecular and General Genetics MGG 1990, 223:76-86 Kupke T, Hernández-Acosta P, Culiáñez-Macià FA: 4’-phosphopantetheine and coenzyme A biosynthesis in plants Journal of Biological Chemistry 2003, 278:38229-38237 Chernys JT, Zeevaart JAD: Characterization of the 9-cis-epoxycarotenoid dioxygenase gene family and the regulation of abscisic acid biosynthesis in avocado Plant Physiology 2000, 124:343-354 Clegg MT, Davis JW: Molecular genetics of avocado California Avocado Society Yearbook 1989, 73:59-61 Chang S, Puryear J, Cairney J: A simple and efficient method for isolating RNA from pine trees Plant Molecular Biology Reporter 1993, 11:113-116 Davis J, Henderson D, Kobayashi M, Clegg MT: Genealogical relationships among cultivated avocado as revealed through RFLP analyses Journal of Heredity 1998, 89:319-323 Mhameed S, Sharon D, Kaufman D, Lahav E, Hillel J, Degani C, Lavi U: Genetic relationships within avocado (Persea americana Mill) cultivars and between Persea species Theoretical and Applied Genetics 1997, 94:279-286 Chanderbali AS, Albert VA, Ashworth VETM, Clegg MT, Litz RE, Soltis DE, Soltis PS: Persea americana (avocado): bringing ancient flowers to fruit in the genomics era BioEssays 2008, 30:386-396 Acheampong KA, Akromah R, Ofori FA, Takarama JF, Saada D, Bitton I, Lavi U: Genetic characterization of Ghanaian avocados Journal of the American Society for Horticultural Science 2008, 133:801-809 Clegg MT, Henderson D, Durbin M: The development of molecular markers in avocado California Avocado Society Yearbook 1992, 573 Clegg MT, Kobayashi M, Lin JZ: The use of molecular markers in the managment and improvement of avocado (Persea americana Mill.) Revista Chapingo Serie Horticultura 1999, 5:227-231 Lavi U, Hillel J, Vainstein A, Lahav E, Sharon D: Application of DNA fingerprints for identification and genetic analysis of avocado Journal of the American Society for Horticultural Science 1991, 116:1078-1081 Chen H, Ashworth VETM, Xu S, Clegg MT: Quantitative genetic analysis of growth rate in avocado Journal of the American Society for Horticultural Science 2007, 132:691-696 Sharon D, Cregan PB, Mhameed S, Kusharska M, Hillel J, Lahav E, Lavi U: An integrated genetic linkage map of avocado Theoretical and Applied Genetics 1997, 95:911-921 Borrone JW, Schnell RJ, Violi HA, Ploetz RC: Seventy microsatellite markers from Persea americana Miller (avocado) expressed sequence tags Molecular Ecology Notes 2007, 7:439-444 Albert V, Soltis D, Carlson J, Farmerie W, Wall PK, Ilut D, Solow T, Mueller L, Landherr L, Hu Y, et al: Floral gene resources from basal angiosperms for comparative genomics research BMC Plant Biology 2005, 5:5 Douhan G, Fuller E, McKee B, Pond E: Genetic diversity analysis of avocado (Persea americana Miller) rootstocks selected under greenhouse Mahomed and van den Berg BMC Plant Biology 2011, 11:167 http://www.biomedcentral.com/1471-2229/11/167 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 conditions for tolerance to phytophthora root rot caused by Phytophthora cinnamomi Euphytica 2011, 182:209-217 Rookes JE, Wright ML, Cahill DM: Elucidation of defence responses and signalling pathways induced in Arabidopsis thaliana following challenge with Phytophthora cinnamomi Physiological and Molecular Plant Pathology 2008, 72:151-161 Hinch JM, Clarke AE: Callose formation in Zea mays as a response to infection with Phytophthora cinnamomi Physiological Plant Pathology 1982, 21:113-124 Weber APM, Weber KL, Carr K, Wilkerson C, Ohlrogge JB: Sampling the Arabidopsis Transcriptome with Massively Parallel Pyrosequencing Plant Physiol 2007, 144:32-42 Priest HD, Fox SE, Filichkin SA, Mockler TC: Utility of next-generation sequencing for analysis of horticultural crop transcriptomes In International Symposium on Molecular Markers in Horticulture Volume 859 International Society for Horticultural Science Acta Horticulturae; 2010:283-288 Vera JC, Wheat CW, Fescemyer HW, Frilander MJ, Crawford DL, Hanski I, Marden JH: Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing Molecular Ecology 2008, 17:1636-1647 Mhameed S, Hillel J, Lahav E, Sharon D, Lavi U: Genetic association between DNA fingerprint fragments and loci controlling agriculturally important traits in avocado (Persea americana Mill.) Euphytica 1995, 84:81-87 Lavi U, Akkaya M, Bhagwat A, Lahav E, Cregan PB: Methodology of generation and characteristics of simple sequence repeat DNA markers in avocado (Persea americana M.) Euphytica 1994, 80:171-177 Coram TE, Wang M, Chen X: Transcriptome analysis of the wheat-Puccinia striiformis f sp tritici interaction Molecular Plant Pathology 2008, 9:157-169 Wang X, Tang C, Deng L, Cai G, Liu X, Liu B, Han Q, Buchenauer H, Wei G, Han D, et al: Characterization of a pathogenesis-related thaumatin-like protein gene TaPR5 from wheat induced by stripe rust fungus Physiologia Plantarum 2010, 139:27-38 Liu J-J, Sturrock R, Ekramoddoullah A: The superfamily of thaumatin-like proteins: its origin, evolution, and expression towards biological function Plant Cell Reports 2010, 29:419-436 Filippov A, Kuznetsova M, Rogozhin A, Spiglazova S, Smetanina T, Belousova M, Kamionskaya A, Skryabin K, Dolgov S: Increased resistance to late blight in transgenic potato expressing thaumatin II gene Ninth workshop of an European Network for development of an integrated control strategy of potato late blight 2005, 263-267 Thomma BPHJ, Eggermont K, Penninckx IAMA, Mauch-Mani B, Vogelsang R, Cammue BPA, Broekaert WF: Separate jasmonate-dependent and salicylate-dependent defense-response pathways in Arabidopsis are essential for resistance to distinct microbial pathogens Proceedings of the National Academy of Sciences 1998, 95:15107-15111 Ekramoddoullah AKM, Xueshu Y, Rona S, Arezoo Z, Doug T: Detection and seasonal expression pattern of a pathogenesis-related protein (PR-10) in Douglas-fir (Pseudotsuga menziesii) tissues Physiologia Plantarum 2000, 110:240-247 Thomas SW, Glaring MA, Rasmussen SW, Kinane JT, Oliver RP: Transcript profiling in the barley mildew pathogen Blumeria graminis by serial analysis of gene expression (SAGE) Molecular Plant-Microbe Interactions 2002, 15:847-856 Fan F, Li X-W, Wu Y-M, Xia Z-S, Li J-J, Zhu W, Liu J-X: Differential expression of expressed sequence tags in alfalfa roots under aluminium stress Acta Physiologiae Plantarum 2011, 33:539-546 Perrigault M, Tanguy A, Allam B: Identification and expression of differentially expressed genes in the hard clam, Mercenaria mercenaria, in response to quahog parasite unknown (QPX) BMC Genomics 2009, 10:377 Harvey J, Lincoln J, Gilchrist D: Programmed cell death suppression in transformed plant tissue by tomato cDNAs identified from an Agrobacterium rhizogenes-based functional screen Molecular Genetics and Genomics 2008, 279:509-521 Morel JB, Dangl JL: The hypersensitive response and the induction of cell death in plants Cell Death and Differentiation 1997, 4:671-683 Sauter M, Rzewuski G, Marwedel T, Lorbiecke R: The novel ethyleneregulated gene OsUsp1 from rice encodes a member of a plant protein Page 14 of 14 56 57 58 59 60 61 62 63 64 65 66 family related to prokaryotic universal stress proteins Journal of Experimental Botany 2002, 53:2325-2331 Sun H-Q, Kwiatkowska K, Yin HL: Actin monomer binding proteins Current Opinion in Cell Biology 1995, 7:102-110 Opalski KS, Schultheiss H, Kogel K-H, Hückelhoven R: The receptor-like MLO protein and the RAC/ROP family G-protein RACB modulate actin reorganization in barley attacked by the biotrophic powdery mildew fungus Blumeria graminis f.sp hordei The Plant Journal 2005, 41:291-303 Kobayashi I, Kobayashi Y, Hardham AR: Dynamic reorganization of microtubules and microfilaments in flax cells during the resistance response to flax rust infection Planta 1994, 195:237-247 Kobayashi Y, Kobayashi I, Funaki Y, Fujimoto S, Takemoto T, Kunoh H: Dynamic reorganization of microfilaments and microtubules is necessary for the expression of non-host resistance in barley coleoptile cells The Plant Journal 1997, 11:525-537 Schütz I, Gus-Mayer S, Schmelzer E: Profilin and Rop GTPases are localized at infection sites of plant cells Protoplasma 2006, 227:229-235 Schmelzer E: Cell polarization, a crucial process in fungal defence Trends in Plant Science 2002, 7:411-415 Büschges R, Hollricher K, Panstruga R, Simons G, Wolter M, Frijters A, van Daelen R, van der Lee T, Diergaarde P, Groenendijk J, et al: The barley Mlo gene: a novel control element of plant pathogen resistance Cell 1997, 88:695-705 Monteiro S, Barakat M, Picarra-Pereira MA, Teixeira AR, Ferreira RB: Osmotin and thaumatin from grape: a putative general defence mechanism against pathogenic fungi Phytopathology 2003, 93:1505-1512 Ballvora A, Ercolano MR, Weiß J, Meksem K, Bormann CA, Oberhagemann P, Salamini F, Gebhardt C: The R1 gene for potato resistance to late blight (Phytophthora infestans) belongs to the leucine zipper/NBS/LRR class of plant resistance genes The Plant Journal 2002, 30:361-371 Perazzolli M, Faccin S, Ciccotti AM, Schwarz F, Moser M, De Luca F, Velasco R, Gessler C, Pertot I, Moser C: Transcriptional analysis of the grape defence response against the root rot agent Armellaria mellea IXth International Conference on Grape Genetics and Breeding 2009, 827:619-622 Avrova AO, Taleb N, Rokka VM, Heilbronn J, Campbell E, Hein I, Gilroy EM, Cardle L, Bradshaw JE, Stewart HE, et al: Potato oxysterol binding protein and cathepsin B are rapidly up-regulated in independent defence pathways that distinguish R gene-mediated and field resistances to Phytophthora infestans Molecular Plant Pathology 2004, 5:45-56 doi:10.1186/1471-2229-11-167 Cite this article as: Mahomed and van den Berg: EST sequencing and gene expression profiling of defence-related genes from Persea americana infected with Phytophthora cinnamomi BMC Plant Biology 2011 11:167 Submit your next manuscript to BioMed Central and take full advantage of: • Convenient online submission • Thorough peer review • No space constraints or color figure charges • Immediate publication on acceptance • Inclusion in PubMed, CAS, Scopus and Google Scholar • Research which is freely available for redistribution Submit your manuscript at www.biomedcentral.com/submit ... doi:10.1186/1471-2229-11-167 Cite this article as: Mahomed and van den Berg: EST sequencing and gene expression profiling of defence-related genes from Persea americana infected with Phytophthora cinnamomi BMC Plant Biology... pathogenesis-related protein, a mlo transmembrane protein and profilin Nine genes were quantified with qRT-PCR to elucidate the early gene response of a tolerant avocado rootstock infected with. .. between pyrosequencing data and gene expression profiles of defence-related genes Sequence ID GenBank accession number cDNA Library Max qRT-PCR expression Similarities between 454 and qRTPCR Thaumatin