Transcriptome response analysis of Arabidopsis thaliana to leafminer (Liriomyza huidobrensis)
© Zhang et al; licensee BioMed Central Ltd. 2012
Received: 26 May 2012
Accepted: 30 November 2012
Published: 11 December 2012
Plants have evolved a complicated resistance system and exhibit a variety of defense patterns in response to different attackers. Previous studies have shown that responses of plants to chewing insects and phloem-feeding insects are significantly different. Less is known, however, regarding molecular responses to leafminer insects. To investigate plant transcriptome response to leafminers, we selected the leafminer Liriomyza huidobrensis, which has a special feeding pattern more similar to pathogen damage than that of chewing insects, as a model insect, and Arabidopsis thaliana as a response plant.
We first investigated local and systemic responses of A. thaliana to leafminer feeding using an Affymetrix ATH1 genome array. Genes related to metabolic processes and stimulus responses were highly regulated. Most systemically-induced genes formed a subset of the local response genes. We then downloaded gene expression data from online databases and used hierarchical clustering to explore relationships among gene expression patterns in A. thaliana damaged by different attackers.
Our results demonstrate that plant response patterns are strongly coupled to damage patterns of attackers.
Plants have evolved a complicated resistance system to defend against damages from various types of attackers. Based on many studies devoted to plant defense signal transduction, three main plant defense hormones have been identified. They are salicylic acid (SA), jasmonic acid (JA), and ethylene (ET), which are the key signal molecules involved in defense against pathogens, insects, and fungi, respectively [1–3]. It has recently been shown, however, that these signals are correlated in a very complex fashion; sometimes they conflict, while at other times they cooperate [4, 5], indicating that plants express various defense patterns when damaged by different attackers. The factors determining these plant response patterns are still not clear, however.
Plant responses to chewing insects and phloem-feeding insects are significantly different [6, 7]. The two types of insects not only produce different elicitors, but also have different feeding guides. For example, wounding leads to leakage of plant cellular liquids, stimulating the mobilization of many defense pathways ; insect feeding causes similar damage to plants, but the elicitors in insect saliva can induce special plant defense proteins  or conversely suppress plant defense signals . Phloem-feeding insects cause little wounding but have long damage durations, and plant defenses to these insects are thus slight [11, 12]. In comparison to chewing insects and phloem-feeding insects, little is known about molecular responses to leafminers, which are insects with special feeding guides.
Pea leafminers (Liriomyza huidobrensis) feed on over 100 species in 22 plant families, including the model plant Arabidopsis thaliana. During the adult stage, the female fly uses her ovipositor to penetrate the epidermis of host plant leaves; she then either lays eggs inside the leaves or feeds at the wound site, which can greatly reduce photosynthesis and eventually kill young plants . Although male flies are unable to puncture leaves, they occasionally feed at the wounds and oviposition punctures made available by females . Plant cells around the oviposition holes usually die and form ayellow necrotic spot, which is similar to necrotic spots created by pathogens . The leafminer is thus a special insect with a damage pattern in plants somewhat similar to pathogens, and is a good model to test relationships between insect damage patterns and plant defense patterns. In addition to its importance as a model insect in plant defense studies, this pest can cause economic losses to host plant crops, as mining larvae consume foliage while dwelling inside leaves . Leafminer larvae consume mesophyll both in palisade and spongy tissues [13, 17]. Leafminer plant damage is very serious, but the insects are difficult to detect at early stages because they are well-hidden in the leaves; it is therefore important to explore the inherent defenses of plants to this insect. Consequently, our study of leafminer-plant interactions is not only important for exploring the mechanisms of plant-insect interactions, but also of value for leafminer pest management.
In this study, we used an Affymetrix ATH1 A. thaliana microarray, from an organism with a well-understood genomic background and thus capable of comprehensively representing transcriptome response, to study expression pattern changes in A. thaliana in response to both local (LI) and systemic (SI) pea leafminer damage. We found that more than 3000 genes were induced in the locally-damaged tissue, and that these genes could be divided into two categories: metabolic processes and stimulus response. Systemic defense of A. thaliana to pea leafminer was very similar to local defense, and the SI-induced genes were almost the same as LI-induced genes, but fewer in number and with lower fold changes. Our analysis of defense signal pathways to leafminer in A. thaliana revealed that signal responses to insects, bacteria, and fungi were all greatly induced. We then downloaded data from online databases and used hierarchical clustering to explore the relationships among A. thaliana expression patterns induced by different types of predators. Interestingly, two different types of data provided evidence that the response to pea leafminer in A. thaliana is more similar to that induced by pathogens than by insects, supporting our hypothesis that plant response patterns are closely related to the damage pattern of attackers.
Microarray expression patterns in leafminer-damaged A. thaliana
Further analysis indicated a high correlation between differentially-expressed genes from SI and LI samples. Although LI and SI expression patterns were very different (Figure 2A), 567 of the 625 genes regulated in SI were differentially-regulated in LI (Figure 2B). The up- and down-regulated genes of SI, with a few exceptions, comprised a subset of the genes differentially-regulated in LI (Figure 2C, D). These results indicate that systemic and local defense responses of A. thaliana to leafminers involve similar mechanisms, but with different gene expression fold change ranges. A GO analysis of genes exclusively regulated in SI revealed that the 47 genes exclusively up-regulated in SI fell into classes focused on transcriptional regulation and stimulus response to elicitors such as hormones and chitin (Additional file 2).
Functional classification of genes regulated by leafminer in A. thaliana
Differentially-expressed genes in both LI and SI were related to metabolism and stimulus response, but the question might be asked as to whether up-regulated genes regulate different subpathways than the down-regulated ones. To determine this, we dissected pathways regulated by genes up- and down-regulated by leafminer damage using EasyGO, a Gene Ontology-based annotation and functional enrichment analysis tool . Nearly all the up-regulated genes locally and systemically expressed in leafminer-damaged A. thaliana were clearly related to defense (Additional files 3 and 4). Although we used a very strict cutoff (P = 1 × 10-5), we found that many biological processes were enriched in up-regulated genes of LI and SI. Pathways directly related to defense were dramatically up-regulated, including responses to other organisms such as bacteria and fungi, wounding, abiotic stimuli such as osmotic stress, water, and cold, and chemical stimuli, including JA, ET, and chitin. Even metabolic processes that were enriched were related to stimulus response, such as cellular aromatic compound metabolic processes, typically indole and derivative metabolic processes. Response to abscisic acid stimulus and amino acid derivative metabolic processes were markedly up-regulated specifically in LI (Additional file 3). Down-regulated pathways were not as dramatically regulated as the up-regulated pathways (Additional files 5 and 6). Using less stringent cutoffs (P=0.00001 for LI, P=0.001 for SI) than those used in analysis of up-regulated genes, we found that the pathways of LI down-regulated genes were concentrated in the categories of cell surface receptor linked signal transduction, pigment biosynthetic processes, response to auxin stimulus, and response to some abiotic stimuli, such as heat and light (Additional file 5); SI down-regulated genes were enriched in response to temperature stimulus, secondary metabolic processes, and amino acid and derivative metabolic processes, the latter of which was similar to the locally down-regulated pathways (Additional file 6).
Comparison of expression patterns in A. thalianadamaged by leafminer and other typical plant damage organisms
Functional analysis of differentially-expressed genes in leafminer-damaged A. thaliana revealed, in addition to regulation of wounding and abiotic stimulus response genes, that many pathogen defense genes were also up-regulated, indicating that leafminers are different from other insects. Leafminers were already known to be unusual because their damage pattern is similar to that of pathogens. The question then arose: is the expression pattern from leafminer-damaged A. thaliana more similar to that of A. thaliana damaged by pathogens or damaged by other insects?
When we analyzed local and systemic expression patterns in leafminer-damaged A. thaliana, we found that there was a trade-off between the different signal pathways expressed. Many defense genes were up-regulated, whereas some stimulus response genes were down-regulated. After carefully analyzing these pathways, we found that the up-regulated genes were mainly involved in biotic stress response, while the down-regulated defense genes were related to abiotic stress. It can thus be seen that the plant defense system can deal with specific stimuli while simultaneously down-regulating other defense pathways to save energy.
Appropriate damage pattern recognition is very important for plants to express suitable defense genes. For example, the defense of A. thaliana to wounding is stronger than to Pieris rapae damage , because P. rapae can reduce tissue crushing and minimize cut leaf edges while removing maximum tissue mass. Indeed, two sets of variables play important roles in the host response: the exact nature of the physical injury and the extent to which elicitors are exposed to the host . Many insects and pathogens try to escape plant defenses by minimizing these two parameters [10, 11]. In this sense, the leafminer is not a good model of successful insect evolution, because the leaf-mining habit does not effectively diminish these two parameters [16, 22]. Finally, it should be noted that disease incidence in leaf-mining insects is remarkably lower than in external-feeding insects ; this suggests that the similarity of gene expression patterns between leafminer-damaged and bacterial pathogen-damaged A. thaliana is not caused by leafminer-derived disease.
Plant transcriptome responses to different attackers may be influenced not only by the attacker species, but also by plant physiological conditions. For example, inbred Solanum carolinense seedlings show markedly weaker defenses to herbivores than do outbred seedlings . When downloading data from TAIR databases and other published sources, we were therefore careful to select data from plants of the same genetic background and with similar physiological conditions and ages. Because responses of A. thaliana to leafminers and pathogens have similar patterns, just as their similar damage patterns we conclude that damage patterns of attackers play an important role in eliciting plant responses. Further investigation indicates secondary metabolite metabolism and some plant physiological responses, but not plant hormones, give rise to the defense pattern correlations among different attackers. In a previous study, A. thaliana response to silverleaf whitefly (Bemisia tabaci type B [SLWF]) involved regulation of several biotrophic pathogen defense pathways ; however, silverleaf whitefly damage strongly up-regulated SA signal pathways, whereas in our study leafminers primarily up-regulated JA signal pathways. Consequently, plant defense patterns to different attackers are shaped not only by interactions between different plant hormone signal pathways, but also different combinations of defense signal pathways, secondary metabolite metabolism, and important physiological responses.
We studied transcriptome response characteristics of A. thaliana to leafminer (L. huidobrensis) damage, and then analyzed its relationships with the response to other attackers. We found that leafminers induced many defense- and metabolism-related genes, with some associated with pathogen defense. Further analysis indicated that the response pattern of A. thaliana to leafminers was most similar to the response pattern of A. thaliana to bacterial pathogens, which was consistent with the similar damage patterns of leafminers and bacterial pathogens. The type of damage pattern thus appears to be an important determinant of plant response patterns.
Plants and insects
Seeds of A. thaliana (Col-0 ecotype) derived from selfed progenies of an Arabidopsis plant were surface-sterilized for 15 min in 10% bleach, washed four times with sterile water, and plated on half-strength Murashige and Skoog  medium. Plants were stratified at 4°C for 2 d in darkness and then transferred to a phytotron set at 22°C under a 16-h light/8-h dark photoperiod (light intensity 120 μmol m-2 s-1). After two to three weeks, seedlings were transplanted into a mixed peat-vermiculite (1:2) potting medium and placed in a growth chamber at 22°C under a 16-h light/8-h dark cycle (light intensity 120 μmol m-2 s-1). Four-week-old A. thaliana seedlings were used in subsequent experiments.
Seeds of Phaseolus vulgaris L. (cv. Naibai; Haizhong Vegetable Market, Beijing, China) were individually sown in 12-cm-diameter plastic pots containing a mixed peat-vermiculite (3:1) potting medium in an environmental chamber. Bean plants with two fully-developed true leaves were used as leafminer reproductive hosts.
Pea leafminers (L. huidobrensis) used in experiments had been under continuous culture for three years in our laboratory. During the course of the experiments, newly-emerged leafminers were reared on 10% honey for 2 d to permit mating before release onto A. thaliana plants.
Plant treatment and cDNA sample preparation
Approximately 50 mated L. huidobrensis adults were released for oviposition onto leaves of four-week-old A. thaliana plants in eight pots (four seedlings per pot). The A. thaliana seedlings were immature and had not yet begun to flower. Adult leafminers were removed within 4 h, by which time approximately half of the leaves had experienced damage. When leafminer larvae had reached the second instar stage, 96 h after oviposition, damaged leaves of A. thaliana containing leafminer larvae were collected for local damage analysis (LI); intact leaves adjacent to the damaged leaves on the same plant were collected for systemic damage analysis (SI). The samples were separately frozen in liquid nitrogen. Three replicate samples, each containing tissue from at least eight plants, were prepared for both LI and SI. Three control samples (H) were also prepared from leaves of same-aged healthy plants. Total RNAs were isolated from each replicate and control sample using an RNeasy Plant Mini Kit (Qiagen, Valencia, CA, USA). cDNA was synthesized from the isolated RNA following the manufacturer’s instructions (http://www.affymetrix.com/support/technical/manuals.affx). Expression Analysis Technical Manual For HT Array Plates Using the GeneChip® Array Station.
Microarray hybridization and data analysis
Affymetrix microarrays (A. thaliana ATH1 genome arrays) containing 22,810 probe sets were used in our experiments. Labeling and hybridization of the ATH1 microarrays (one sample per chip) was performed according to the manufacturer’s instructions (http://www.affymetrix.com/support/technical/manuals.affx). Expression Analysis Technical Manual For HT Array Plates Using the GeneChip® Array Station. Quality control was carried out according to Affymetrix microarray standards, which indicated that all samples met required criteria with respect to control signals, house-keeping gene signals, “present” percent numbers, and low background and noise values. The probe arrays were scanned and further analyzed using GENESPRING software (version 5.0; Silicon Genetics). Normalization per gene and per chip of the log2 values was performed to allow comparison of the three independent replicates performed for each set of experiments. In addition, normalization was performed separately for each experiment and plant tissue for all measurements using the flags “present”, “marginal”, or “absent” assigned by Affymetrix treatment of the arrays. However, only those transcripts that were declared “present” or “marginal” in at least two of three chips were taken into account. To identify differentially-expressed genes in each treatment, SAM analysis (Significance Analysis of Microarrays software package) was conducted on A. thaliana triplicate samples between treatments and controls using a q-value ≤ 0.05 and fold change ≥ 2 as cut-off criteria. We searched for GO information for the differently-expressed probe sets using EasyGO software (http://bioinformatics.cau.edu.cn/easygo/category_treeBrowse.html). For biological process searching, we applied χ2 tests with a false discovery rate (FDR)-adjusted cutoff of P < 0.0001. Cluster 3.0/Treeview software  was used to group and display genes with similar expression profiles (http://rana.stanford.edu/software/). Hierarchical clustering using default options and the uncentered correlation similarity metric was performed on the normalized data. The microarray data in MIAME-compliant format have been deposited in the GEO database (GEO record number GSE38281).
Downloading of microarray data from public databases
Transcriptome datasets used for expression cluster analysis of different damage patterns were downloaded from the Arabidopsis Information Resource (TAIR) public GeneChip database (ftp://ftp.arabidopsis.org/home/tair/Microarrays/Datasets/). Differential genes from each experiment were identified using the same method used for our own data. Detailed information regarding downloaded datasets from the TAIR microarray database is listed in Additional file 10. Other data was derived from the results of a published study .
Quantitative real-time PCR
Primers used for qRT-PCR to validate fold change of microarray hybridization
Sense primer 5′-3′
Anti-sense primer 5′-3′
We thank CapitalBio Corporation for assistance with Affymetrix microarray technology. This research was supported by the National Basic Research Program of China (973 Program) (No. 2012CB114105), and National Nature Science Foundation of China (30921063), Institute Special Fund for Basic Research, Institute of Forest Ecology, Environment, and Protection, Chinese Academy of Forestry (CAFRIFEEP201102-5).
- Ellis C, Turner JG: The Arabidopsis mutant cev1 has constitutively active jasmonate and ethylene signal pathways and enhanced resistance to pathogens. Plant Cell. 2001, 13 (5): 1025-1033.PubMedPubMed CentralView ArticleGoogle Scholar
- Delaney TP, Uknes S, Vernooij B, Friedrich L, Weymann K, Negrotto D, Gaffney T, Gut-Rella M, Kessmann H, Ward E, et al: A central role of salicylic acid in plant disease resistance. Science. 1994, 266 (5188): 1247-1250. 10.1126/science.266.5188.1247.PubMedView ArticleGoogle Scholar
- Turner JG, Ellis C, Devoto A: The jasmonate signal pathway. Plant Cell. 2002, 14: S153-S164.PubMedPubMed CentralGoogle Scholar
- Doares SH, Narvaezvasquez J, Conconi A, Ryan CA: Salicylic-acid inhibits synthesis of proteinase-inhibitors in tomato leaves induced by systemin and jasmonic acid. Plant Physiol. 1995, 108 (4): 1741-1746.PubMedPubMed CentralGoogle Scholar
- Reymond P, Farmer EE: Jasmonate and salicylate as global signals for defense gene expression. Curr Opin Plant Biol. 1998, 1 (5): 404-411. 10.1016/S1369-5266(98)80264-1.PubMedView ArticleGoogle Scholar
- Zheng S-J, Dicke M: Ecological genomics of plant-insect interactions: from gene to community. Plant Physiol. 2008, 146 (3): 812-817. 10.1104/pp.107.111542.PubMedPubMed CentralView ArticleGoogle Scholar
- Kempema LA, Cui X, Holzer FM, Walling LL: Arabidopsis transcriptome changes in response to phloem-feeding silverleaf whitefly nymphs. Similarities and distinctions in responses to aphids. Plant Physiol. 2007, 143 (2): 849-865.PubMedPubMed CentralView ArticleGoogle Scholar
- Cheong YH, Chang H-S, Gupta R, Wang X, Zhu T, Luan S: Transcriptional profiling reveals novel interactions between wounding, pathogen, abiotic stress, and hormonal responses in Arabidopsis. Plant Physiol. 2002, 129 (2): 661-677. 10.1104/pp.002857.PubMedPubMed CentralView ArticleGoogle Scholar
- Korth KL, Dixon RA: Evidence for chewing insect-specific molecular events distinct from a general wound response in leaves. Plant Physiol. 1997, 115 (4): 1299-1305.PubMedPubMed CentralGoogle Scholar
- Reymond P, Weber H, Damond M, Farmer EE: Differential gene expression in response to mechanical wounding and insect feeding in Arabidopsis. Plant Cell. 2000, 12 (5): 707-720.PubMedPubMed CentralView ArticleGoogle Scholar
- Moran PJ, Thompson GA: Molecular responses to aphid feeding in Arabidopsis in relation to plant defense pathways. Plant Physiol. 2001, 125 (2): 1074-1085. 10.1104/pp.125.2.1074.PubMedPubMed CentralView ArticleGoogle Scholar
- Thompson GA, Goggin FL: Transcriptomics and functional genomics of plant defence induction by phloem-feeding insects. J Exp Bot. 2006, 57 (4): 755-766. 10.1093/jxb/erj135.PubMedView ArticleGoogle Scholar
- Wei JN, Zhu JW, Kang L: Volatiles released from bean plants in response to agromyzid flies. Planta. 2006, 224 (2): 279-287. 10.1007/s00425-005-0212-x.PubMedView ArticleGoogle Scholar
- Parrella MP, Bethke JA: Biological Studies of Liriomyza-Huidobrensis (Diptera, Agromyzidae) on Chrysanthemum, Aster, and Pea. J Econ Entomol. 1984, 77 (2): 342-345.View ArticleGoogle Scholar
- Parrella MP: Biology of Liriomyza. Annu Rev Entomol. 1987, 32: 201-224. 10.1146/annurev.en.32.010187.001221.View ArticleGoogle Scholar
- Connor EF, Taverner MP: The evolution and adaptive significance of the leaf-mining habit. Oikos. 1997, 79 (1): 6-25. 10.2307/3546085.View ArticleGoogle Scholar
- Wei JN, Zou L, Kuang RP, He LP: Influence of leaf tissue structure on host feeding selection by pea leafminer Liriomyza huidobrensis (Diptera: Agromyzidae). Zool Stud. 2000, 39 (4): 295-300.Google Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34 (Web Server issue): W293-W297.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhou X, Su Z: EasyGO: gene ontology-based annotation and functional enrichment analysis tool for agronomical species. BMC Genomics. 2007, 8: 246-10.1186/1471-2164-8-246.PubMedPubMed CentralView ArticleGoogle Scholar
- De Vos M, Van Oosten VR, Van Poecke RM, Van Pelt JA, Pozo MJ, Mueller MJ, Buchala AJ, Metraux JP, Van Loon LC, Dicke M, et al: Signal signature and transcriptome changes of Arabidopsis during pathogen and insect attack. Mol Plant Microbe Interact. 2005, 18 (9): 923-937. 10.1094/MPMI-18-0923.PubMedView ArticleGoogle Scholar
- Farmer EE: Adding injury to insult: pathogen detection and responses. Genome Biol. 2000, 1 (2): 1012.1-1012.3.View ArticleGoogle Scholar
- Engelbrecht L, Orban U, Heese W: Leaf-miner caterpillars and cytokinins in the green Islands of Autumn leaves. Nature. 1969, 223 (5203): 319-321. 10.1038/223319a0.View ArticleGoogle Scholar
- Kariyat RR, Mena-Alí J, Forry B, Mescher MC, Moraes CMD, Stephenson AG: Inbreeding, herbivory, and the transcriptome of solanum carolinense. Entomol Exp Appl. 2012, 144: 134-144. 10.1111/j.1570-7458.2012.01269.x.View ArticleGoogle Scholar
- Murashige T, Skoog F: A Revised Medium for rapid growth and bio asssays with tobacco tissue cultures. Physiol Plant. 1962, 15 (3): 473-497. 10.1111/j.1399-3054.1962.tb08052.x.View ArticleGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to 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.