BACKGROUND
once a plant recognizes a pathogen attack, the metabolism must balance demands for resources to support defense versus requirements for cellular maintenance, growth and reproduction  <cit> .
defense mechanisms involve a shift on metabolism, activating specific pathways such as synthesis of secondary metabolites, programmed cell death, ions translocation.
concomitantly, can occur a shutdown of other metabolic pathways not directly involved with defense response, such as those associated with growth and reproduction.
a resistance and/or tolerance trait is attributed whenever this defense response is successful in controlling pathogen or herbivore attack.
genetic control of metabolic re-programming is normally triggered by few resistance genes which are seek out to be transferred to other plants.
however, resistance response involves changes in other genes, not usually identified, with determinant roles in the overall response.
therefore for an effective transference of resistance traits, to know how these genes interact during re-programming of plant metabolism is essential.
among available methods for high-throughput analysis the microarray is a powerful tool for large-scale gene expression studies in many plant species with whole genome sequenced: potato  <cit> , tomato  <cit> , soybean  <cit> , wheat  <cit> , barley  <cit>  maize  <cit> , grape  <cit> , pine  <cit> , arabidopsis  <cit> .
the main advantage of microarray analyses is to evaluate the expression of large number of genes in different genotypes, organs, tissues, treatments, using the same set of genes.
these genes can be compared during different biological situations allowing both an association with metabolic pathways and establishement of their role on resistance response.
several studies have been carried out using microarray analysis to identify genes associated with plant defense  <cit> .
the leaf-miner, leucoptera coffeella   is a specialist parasite of coffea species.
upon oviposition on leaves, ecloded larvae feed directly from parenchyma tissues, leading to a reduction of foliar surface and eventual leaf drop  <cit> .
this damage results in reduction of photosynthetic area and plant survival.
in brazilian breeding programs resistance genes from c. racemosa have been transferred to the susceptible c. arabica by controlled crosses, and so far a large number of hybrid progenies are under selection for resistance to leaf-miner  <cit> .
although defense mechanism to leaf-miner is not understood yet, genetic analysis demonstrated that resistance to the insect is dominant and controlled by two complimentary genes  <cit> .
at the molecular level, there is little information regarding gene expression on coffee plants during defense response.
using subtractive hybridization methodology , mondego et al.
<cit>  found differentially expressed genes in coffee plants upon leaf-miner infestation, among which a miraculin-like encoding gene was significantly overexpressed in resistant coffee plants.
differential expression of defense-related genes such as lipoxygenase, glutathione transferase, protein-kinase receptor and glucanase was observed in response to leaf-miner infestation  <cit> .
however, the expression profiles indicate that differences results from gene expression timing along insect infection rather than with gene regulation.
despite the efforts of breeding programs to develop novel coffee cultivars bearing leaf-miner resistance, selection of progenies homozygotes for this trait is difficult  <cit> , as advanced generations are still producing susceptible plants.
therefore, information regarding molecular control of resistance response as well as identification of candidate-genes associated with these processes will contribute with assisted-selection.
in this context, the aim of our study was to explore transcriptomic differences throughout insect infestation, in susceptible and resistant c. arabica plants challenged by l. coffeella, using microarray technology.
the arrays were developed using coffee-specific oligoprobes designed based on gene sequences available at the brazilian coffee genome project  <cit> .
the database contains a collection of around  <dig>  gene sequences, covering most of the c. arabica genome  <cit> .
besides this, we selected a group of candidate-genes to be used as molecular markers for assisted-selection.
as far as we know, this is the first report of a large-scale transcriptional profile analysis used to study gene expression changes in coffee plants in the presence of an herbivore insect.
RESULTS
microarray analyses were performed to characterize large-scale gene expression profiles during leaf-miner development on coffee leaves.
the analyses included a hybridization of a 135 k array with  <dig> different samples, corresponding to time-course infestation stages in both resistant and susceptible plants.
the arrays contain sequences of around 33 k genes identified in est libraries prepared from different physiological and metabolic situations  <cit> .
a minimum of  <dig> - 24mer match probes for each selected gene were used for the array set up.
the arrays were hybridized with probes corresponding to  <dig> treatments of both susceptible  and resistant  leaves: non-infestated , after oviposition and egg-eclosion  and damaged by insect feeding .
initially, differential expression patterns were identified using statistical analysis, and specific transcriptional profiles were established for each evaluated interaction.
in a second approach, genes exhibiting differential expression among genotypes and treatments were submitted to in silico evaluations to classify and categorize those genes regarding their possible molecular functions and metabolic pathways.
finally, a group of  <dig> genes involved with defense-related mechanisms, exhibiting regulated expression, were further characterized using qrt-pcr.
microarray and statistical analysis
a total of  <dig>   <dig> and  <dig> differentially expressed genes were identified from in silico analyses of raw hybridization data considering  <dig> interactions: t0r x t0s, t1r x t1s and t2r x t2s .
comparing t0r x t0s we observed higher differential expression levels where foldchange values range from  <dig> and  <dig> times in up-regulated genes, and from  <dig> to  <dig> times in down-regulated genes in resistant leaves.
the other interactions exhibited foldchange values ranging from  <dig> and  <dig> times in down-regulated, and from  <dig> to  <dig> times in up-regulated genes .
different regulation profiles for defense response were observed among differentially expressed genes.
the first group includes chitinase and polygalacturonase genes, regulated differently in resistant and susceptible leaves from t <dig> to t <dig>  the second group includes genes that were up or down regulated in control resistant or susceptible leaves, but throughout insect development inverted their expression pattern.
example of this group is the gene encoding the enzyme polyphenol oxidase, up regulated at t <dig> in resistant plants and down regulated at t <dig> and t <dig>  and a third group includes genes that were up and/or down regulated in response to the leaf-miner infestation, either in susceptible or resistant leaves.
interaction between resistant and susceptible genotypes without insect infestation
since higher values for differential gene expression were observed when comparing t0r x t0s samples, we chose this interaction for further analyses and selection of candidate-genes for validation.
in this interaction,  <dig> genes exhibited differential expression,  <dig> were up regulated and  <dig> were down regulated.
regulated genes from the t <dig> samples were functionally characterized into three gene ontology categories – molecular function, component cellular and biological function - and grouped according their metabolic categories .
a description of the first one hundred most variable genes, both up and down regulated, is shown on tables  <dig> and  <dig> and illustrated in figures  <dig> and  <dig>  several contigs had no correspondence with defined categories.
a larger number of genes are up-regulated rather than down-regulated, and associated with primary and cellular metabolism, which included functions such as ion, protein and nucleic acid binding, hydrolase and transferase activities, among others.
as these differences were observed in non-infested leaves, possibly a different transcriptional programming takes place in resistant leaves.
this may result in a pre-defense status, enabling resistant plants to a faster defense response upon leaf-miner attack.
*the annotation of each sequence was established on the coffee genome database  <cit> .
**no match with any sequence on the solanaceae genomic database.
*the annotation of each sequence was established on the coffee genome database .
**no match with any sequence on the solanacea genomic database.
metabolic pathways
the categorization of annotated genes using kegg database was performed with the first  <dig> up and down regulated genes .
most of these genes are from starch and purine metabolism, and several others are involved in primary metabolism.
three main metabolic pathways are highlighted here: citrate metabolism, linoleic acid metabolism and phenylpropanoids metabolism.
we choose the citric acid cycle for further characterization as previous analyses using nmr indicated that lower levels of malate, a metabolite resulting from the citrate metabolism, are observed in resistant coffee leaves  <cit> .
several genes encoding citric acid cycle enzymes exhibit differential-expression .
expression of isocitrate lyase gene is repressed in resistant genotypes at t <dig> , suggesting that synthesis of malate may be deficient, and therefore low levels of malate may accumulate in those leaves.
however, this gene is up-regulated upon oviposition and egg ecclosion .
the linoleic acid is the first substrate of the jasmonic acid  pathway, a major signaling pathway during herbivore-defense responses.
control resistant plants  show up-regulation of jasmonate o-methyltransferase and lipoxygenase while differential expression for these genes was not observed at any stage in susceptible genotypes .
also,  <dig> genes from the alpha-linoleic acid metabolism and  <dig> genes from jasmonate biosynthesis were regulated in resistant plants .
for instance, transcripts of enoyl-coa hydratase and phospholipase a <dig> were four times more expressed in t <dig> than in t <dig> in resistant genotypes , but increased only at later stages in susceptible leaves.
phenylpropanoids are major plant phtytoalexins, part of the secondary metabolism .
twenty-seven genes from phenylpropanoids synthesis exhibited differential expression at t <dig>  with foldchange values ranging from  <dig> to - <dig>  transcript levels of phenylalanine ammonia lyase , the enzyme that catalyzes the first step of the pathway, is up-regulated only at t <dig> in resistant plants , and this level decreases along insect development.
in susceptible plants, pal transcript levels increase at final steps of insect infection, t <dig> .
however, genes fromlignin and isoflavones synthesis, downstream metabolites, such as cynnamyl alcoholdehydrogenase  and isoflavone reductase , are down-regulated in resistant plants.
on the other hand, genes from biosynthesis of anthocyanins and tannins, such as flavonoide 3’-hydroxylase  and leucoanthocyanidin dioxygenase  are up-regulated at all times in resistant plants.
this expression profile suggests that synthesis of anthocyanins and tannins is favored over synthesis of flavones.
also, activation of upstream genes such as phenylalanine ammonia lyase, chalcone synthase and flavonone dehydrogenase is observed at final stages in susceptible plants, indicating that phenylpropanoid biosynthesis is delayed.
several other defense-related genes are also positively regulated in resistant plants, including herbivore-response related genes glutathione-s-transferase and cysteine proteinase inhibitor.
apoptosis-related genes have a variable expression profile: catalase is up-regulated throughout insect development but citocrome c oxidase, superoxide dismutase and a senescence-associated protein have no differential expression, and polygalacturonase is up-regulated only at t <dig>
validation of expression profile for selected candidate-genes
we selected  <dig> genes for validation, listed with corresponding expression levels on table  <dig>  genes exhibited a consistent expression pattern when quantified by either microarray or qpcr, and the pearson coefficient for this comparison is  <dig>
fold-change
value
*
**
*microarray analysis value.
**qrt-pcr analysis value.
***pearson coefficient correlation value.
further qpcr analyses were performed to validate expression of selected leaf-miner resistance-associated candidate-genes .
these included genes from pathways described above and genes with either no significant hits or similarity to unknown proteins, which may represent coffee specific genes, not yet identified or characterized.
the putative caffeine synthase gene exhibited the greatest response to leaf-miner infection , as its expression was significantly higher  in control resistant leaves, and also in later stages .
the expression of gene sgn-e <dig>  encoding an unknown protein, is significantly higher  in control resistant leaves than in susceptible ones .
however, transcript levels dropped in resistant plants at first stages of infection, such as oviposition  and egg ecclosion .
this gene is a good candidate for differentiation of resistant and susceptible plants, although possibly is not related to resistance.
the isocitrate lyase gene is down-regulated  in resistant leaves at t <dig>  at initial steps of infection its transcript levels increased rapidly , but decreased during feeding stages.
other evaluated genes also exhibited a regulated expression.
the glycerol-3-phosphate acyltransferase  <dig> gene, a trans-membrane protein associated with synthesis of cutin, is up regulated .
a gene encoding metallothionein-like protein, a class of metal-chelators proteins with possible anti-oxidant role, is also up-regulated  at initial stages of insect development.
genes encoding unknown proteins with lipase protein domain, such as sgn-u <dig>  sl <dig> ch <dig> and sgn-u <dig> were activated in resistant plants, being up-regulated at oviposition and/or egg ecclosion with rq values of  <dig> ,  <dig>  and  <dig>  respectively.
as these genes are down-regulated on leaves at t <dig>  this profile during insect development suggests a role during defense response .
the expression profile of polygalacturonase gene was not confirmed by real-time pcr.
according to microarray in silico analysis this gene is activated at t <dig> , and afterwards transcript levels drop along insect development .
yet by real-time pcr analysis transcript levels are similar in susceptible and resistant plants at all times.
this result may reflect differences on genetic background of evaluated plants, once they are part of a population still segregating for some characteristics.
however, as this was the only observed discrepancy between all performed analyses, these genetic differences may not be associated with the resistance trait.
discussion
the use of resistant or tolerant cultivars represents an alternative for reducing the use of chemical defensives, the costs of production, and the negative impact over environment.
in coffee, breeding for herbivore resistance is important once coffea arabica, the main commercial coffea species, is susceptible to almost all known coffee pathogens.
then resistance must be transferred from other compatible species, which is not always possible due to limited efficiency of inter-specific crosses.
at the agronomic institute  leaf-miner resistance genes have been transferred from c. racemosa to c. arabica through traditional breeding strategies, and although the program is currently at advanced generations, no resistant commercial cultivar is available yet.
the lack of knowledge of molecular aspects controlling the resistance response, and the restricted genetic variability of breeding populations  <cit>  limit the development of genomic-based selection tools.
in this context, we aimed to provide information on molecular aspects of leaf-miner defense mechanisms and identification of reliable candidate-marker genes for assisted-selection.
those genomic tools associated with traditional breeding strategies guarantee that agronomical traits such as productivity and cup quality will be selected in advanced generations.
also, once novel genomic-markers associated with other desirable traits are developed for coffee genotypes, a genome-wide selection strategy will be possible to develop multiple-pathogen resistant cultivars.
we chose the microarray analysis approach to assess the co-expression of a large amount of genes, including those that are not looked at in common analyses.
nevertheless, the results described in this work indicated that pathways regularly activated in response to herbivory, such as linoleic acid cycle, phenylpropanoids synthesis and apoptosis, are also activated during coffee-response to leaf-miner.
genes associated with jasmonate  synthesis, such as lipoxygenase and enoyl-coa hydratase, and with flavonoids synthesis, such as chalcone synthase and flavanone 3-hydroxylase-like, are up-regulated in resistant plants.
also, pathways from the primary metabolism, such as the citric acid cycle are down-regulated during leaf-miner defense response in resistant plants, a profile observed also in conifers  <cit> .
results of this study provide evidence that most genes encoding enzymes from the citric acid cycle are down-regulated in resistant plants .
in a parallel analysis, a metabolite profile was established for resistant and susceptible genotypes using an nmr-based technique  <cit>  and indicated that malate levels on resistant leaves are lower than in susceptible ones.
malate results from conversion of either fumarate or glyoxylate.
expression of fumarase, that converts fumarate into malate, is down-regulated at t <dig> in resistant genotypes , what could explain the low malate levels.
production of malate from glyoxylate may also be deficient in resistant plants, once genes encoding for malate synthase, that converts glyoxylate into malate, and for isocitrate liase, the upstream enzyme that converts isocitrate into glyoxylate and succinate, are both down-regulated .
in contrast to this profile, susceptible plants exhibit a regular expression levels for these genes at t <dig>  therefore, both metabolic and transcriptional profiles support the affirmations that citrate cycle is down regulated in leaf-miner resistant coffee plants, and the model of down-regulation of primary metabolism in herbivore-resistant plants  <cit> .
biosynthesis of ja starts with alpha-linoleic acid release in non-injured tissues, triggered by systemin and phospholipase a <dig>  alpha-linoleic acid is then converted to ja after enzymatic steps performed by 13-lipoxygenase , allene oxide synthase , jasmonate o-methyltransferase and others  <cit> .
several genes from the ja biosynthesis pathway are up-regulated in resistant plants at t <dig>  including those from downstream steps such as jasmonate o-methyltransferase which expression is 10-fold higher than in susceptible plants.
all genes of the ja biosynthesis are either down-regulated or up-regulated at later stages in susceptible plants, as for instance expression of lox  increases only at t <dig>  these observations suggest that the ja signalling pathway, including intermediate signaling-molecules such as oxo-pentenyl-cyclopentane , may be impaired in susceptible plants.
down-regulation of genes from later steps of ja biosynthesis, such as allene oxide cyclase, allene oxide synthase, carboxyl methyltransferase, the enzyme that converts jasmonic acid into methyl jasmonate, is observed at t <dig> and t <dig> in resistant plants.
therefore, a feedback regulation may be activated, with a re-programming of transcriptional response upon leaf-miner infection.
genes associated with biosynthesis of secondary compounds are shown here to be regulated.
expression profile of genes from phenylpropanoids biosynthesis, both up-stream genes such as pal and chs, and downstream genes such as flavonoide 3’-hydroxylase, leucoanthocyanidin dioxygenase, reveals a preferential synthesis of tannins and anthocyanins instead of ligninin, flavones and isoflavones in resistant plants.
this profile indicates a direct defense strategy against leaf-miner, once among anthocyanins and tannins are found toxic compounds with antifeedants effects over insects  <cit> .
another gene linked to secondary metabolism is a putative caffeine synthase, which encodes one enzyme from caffeine biosynthetic pathway.
caffeine is an alkaloid distributed in coffee plant tissues and organs.
the fact that expression of a gene from its biosynthesis is significantly increased upon leaf-miner infection suggests that caffeine may have a role in defense response.
however, several studies regarding caffeine and leaf-miner development indicated that this compound has no effect on insect survival rates  <cit> .
once the ultimate goal of this study is to identify potential candidate for markers, several genes were selected for validation using real-time pcr.
potential candidates include: isocitrate lyase, which increased expression during initial steps of leaf-miner infection may be co-related with reduction of primary carbon metabolism; putative caffeine synthase, part of an important pathway of coffee plants; glycerol-3-phosphate acyltransferase  <dig>  a gene associated with lipid metabolism and part of cutin biosynthetic pathway, a secondary metabolite  <cit> ; and finally metallothionein-like protein gene, a metal-transporter protein family with an uncertain role in plant metabolism but previously associated with redox responses  <cit> .
future analyses include cloning and re-sequencing genomic regions of target genes from different genotypes in order to identify suitable polymorphisms.
among selected genes are those that have no similarity with any known reported gene or protein.
although they could not yet be associated with a biological process, their expression profile was very specific and related to defense response.
for instance, genes sgn-e <dig>  sgn-e <dig> sgn-e <dig> were activated upon leaf-miner infection in resistant plants, and therefore represent good candidates for further investigation.
another interesting unknown gene is sgn-e <dig>  which is highly activated in resistant plants at t <dig> but is repressed upon infection.
the expression profile indicates that this gene is useful for early differentiation between resistant and susceptible plants.
in summary, differential expression profiles between resistant and susceptible genotypes are observed even in the absence of leaf-miner, indicating that defense is already build up in resistant plants, as a priming mechanism.
then, a systemic defense response may be more rapidly activated in resistant plants, once basic compounds such as nitrogen and sugars are readily available as a result of repression of primary metabolism.
this shift in plant metabolism is common after a pathogen attack, where defense-related pathways are activated, resulting in reduction of growth and reproduction, and in changes on link-source relationship  <cit> .
during herbivore-defense transcript levels of genes involved in photosynthesis are also down-regulated  <cit> , probably as a strategy to liberate nitrogen compounds for the secondary metabolism.
maintenance of these physiological and metabolic states has a high energetic cost, and could represent a survival limitation if nutritional conditions on the field are depleted.
actually, field observations in cultivated areas demonstrated that leaf-miner resistant coffee plants, under a severe nutritional deficit, are attacked by the leaf-miner at the same intensity as susceptible plants  <cit> .
CONCLUSIONS
as a basal defense state is decisive for triggering a rapid resistance response, genes associated with priming validated here, represent key genes for assisted-selection.
future studies will focus on comparisons of selected genes genomic sequences, from both resistant and susceptible parental lines, to identify suitable marker polymorphisms.

 once a plant recognizes a pathogen attack, the metabolism must balance demands for resources to support defense versus requirements for cellular maintenance, growth and reproduction. defense mechanisms involve a shift on metabolism, activating specific pathways such as synthesis of secondary metabolites, programmed cell death, ions translocation. concomitantly, can occur a shutdown of other metabolic pathways not directly involved with defense response, such as those associated with growth and reproduction. a resistance and/or tolerance trait is attributed whenever this defense response is successful in controlling pathogen or herbivore attack. genetic control of metabolic re-programming is normally triggered by few resistance genes which are seek out to be transferred to other plants. however, resistance response involves changes in other genes, not usually identified, with determinant roles in the overall response. therefore for an effective transference of resistance traits, to know how these genes interact during re-programming of plant metabolism is essential. among available methods for high-throughput analysis the microarray is a powerful tool for large-scale gene expression studies in many plant species with whole genome sequenced: potato, tomato, soybean, wheat, barley maize, grape, pine, arabidopsis. the main advantage of microarray analyses is to evaluate the expression of large number of genes in different genotypes, organs, tissues, treatments, using the same set of genes. these genes can be compared during different biological situations allowing both an association with metabolic pathways and establishement of their role on resistance response. several studies have been carried out using microarray analysis to identify genes associated with plant defense. the leaf-miner, leucoptera coffeella is a specialist parasite of coffea species. upon oviposition on leaves, ecloded larvae feed directly from parenchyma tissues, leading to a reduction of foliar surface and eventual leaf drop. this damage results in reduction of photosynthetic area and plant survival. in brazilian breeding programs resistance genes from c. racemosa have been transferred to the susceptible c. arabica by controlled crosses, and so far a large number of hybrid progenies are under selection for resistance to leaf-miner. although defense mechanism to leaf-miner is not understood yet, genetic analysis demonstrated that resistance to the insect is dominant and controlled by two complimentary genes. at the molecular level, there is little information regarding gene expression on coffee plants during defense response. using subtractive hybridization methodology, mondego et al. found differentially expressed genes in coffee plants upon leaf-miner infestation, among which a miraculin-like encoding gene was significantly overexpressed in resistant coffee plants. differential expression of defense-related genes such as lipoxygenase, glutathione transferase, protein-kinase receptor and glucanase was observed in response to leaf-miner infestation. however, the expression profiles indicate that differences results from gene expression timing along insect infection rather than with gene regulation. despite the efforts of breeding programs to develop novel coffee cultivars bearing leaf-miner resistance, selection of progenies homozygotes for this trait is difficult, as advanced generations are still producing susceptible plants. therefore, information regarding molecular control of resistance response as well as identification of candidate-genes associated with these processes will contribute with assisted-selection. in this context, the aim of our study was to explore transcriptomic differences throughout insect infestation, in susceptible and resistant c. arabica plants challenged by l. coffeella, using microarray technology. the arrays were developed using coffee-specific oligoprobes designed based on gene sequences available at the brazilian coffee genome project. the database contains a collection of around gene sequences, covering most of the c. arabica genome. besides this, we selected a group of candidate-genes to be used as molecular markers for assisted-selection. as far as we know, this is the first report of a large-scale transcriptional profile analysis used to study gene expression changes in coffee plants in the presence of an herbivore insect.  microarray analyses were performed to characterize large-scale gene expression profiles during leaf-miner development on coffee leaves. the analyses included a hybridization of a 135 k array with different samples, corresponding to time-course infestation stages in both resistant and susceptible plants. the arrays contain sequences of around 33 k genes identified in est libraries prepared from different physiological and metabolic situations. a minimum of - 24mer match probes for each selected gene were used for the array set up. the arrays were hybridized with probes corresponding to treatments of both susceptible and resistant leaves: non-infestated, after oviposition and egg-eclosion and damaged by insect feeding. initially, differential expression patterns were identified using statistical analysis, and specific transcriptional profiles were established for each evaluated interaction. in a second approach, genes exhibiting differential expression among genotypes and treatments were submitted to in silico evaluations to classify and categorize those genes regarding their possible molecular functions and metabolic pathways. finally, a group of genes involved with defense-related mechanisms, exhibiting regulated expression, were further characterized using qrt-pcr. microarray and statistical analysis a total of and differentially expressed genes were identified from in silico analyses of raw hybridization data considering interactions: t0r x t0s, t1r x t1s and t2r x t2s. comparing t0r x t0s we observed higher differential expression levels where foldchange values range from and times in up-regulated genes, and from to times in down-regulated genes in resistant leaves. the other interactions exhibited foldchange values ranging from and times in down-regulated, and from to times in up-regulated genes. different regulation profiles for defense response were observed among differentially expressed genes. the first group includes chitinase and polygalacturonase genes, regulated differently in resistant and susceptible leaves from t to t the second group includes genes that were up or down regulated in control resistant or susceptible leaves, but throughout insect development inverted their expression pattern. example of this group is the gene encoding the enzyme polyphenol oxidase, up regulated at t in resistant plants and down regulated at t and t and a third group includes genes that were up and/or down regulated in response to the leaf-miner infestation, either in susceptible or resistant leaves. interaction between resistant and susceptible genotypes without insect infestation since higher values for differential gene expression were observed when comparing t0r x t0s samples, we chose this interaction for further analyses and selection of candidate-genes for validation. in this interaction, genes exhibited differential expression, were up regulated and were down regulated. regulated genes from the t samples were functionally characterized into three gene ontology categories – molecular function, component cellular and biological function - and grouped according their metabolic categories. a description of the first one hundred most variable genes, both up and down regulated, is shown on tables  and and illustrated in figures  and several contigs had no correspondence with defined categories. a larger number of genes are up-regulated rather than down-regulated, and associated with primary and cellular metabolism, which included functions such as ion, protein and nucleic acid binding, hydrolase and transferase activities, among others. as these differences were observed in non-infested leaves, possibly a different transcriptional programming takes place in resistant leaves. this may result in a pre-defense status, enabling resistant plants to a faster defense response upon leaf-miner attack. *the annotation of each sequence was established on the coffee genome database. **no match with any sequence on the solanaceae genomic database. *the annotation of each sequence was established on the coffee genome database. **no match with any sequence on the solanacea genomic database. metabolic pathways the categorization of annotated genes using kegg database was performed with the first up and down regulated genes. most of these genes are from starch and purine metabolism, and several others are involved in primary metabolism. three main metabolic pathways are highlighted here: citrate metabolism, linoleic acid metabolism and phenylpropanoids metabolism. we choose the citric acid cycle for further characterization as previous analyses using nmr indicated that lower levels of malate, a metabolite resulting from the citrate metabolism, are observed in resistant coffee leaves. several genes encoding citric acid cycle enzymes exhibit differential-expression. expression of isocitrate lyase gene is repressed in resistant genotypes at t, suggesting that synthesis of malate may be deficient, and therefore low levels of malate may accumulate in those leaves. however, this gene is up-regulated upon oviposition and egg ecclosion. the linoleic acid is the first substrate of the jasmonic acid pathway, a major signaling pathway during herbivore-defense responses. control resistant plants show up-regulation of jasmonate o-methyltransferase and lipoxygenase while differential expression for these genes was not observed at any stage in susceptible genotypes. also, genes from the alpha-linoleic acid metabolism and genes from jasmonate biosynthesis were regulated in resistant plants. for instance, transcripts of enoyl-coa hydratase and phospholipase a were four times more expressed in t than in t in resistant genotypes, but increased only at later stages in susceptible leaves. phenylpropanoids are major plant phtytoalexins, part of the secondary metabolism. twenty-seven genes from phenylpropanoids synthesis exhibited differential expression at t with foldchange values ranging from to - transcript levels of phenylalanine ammonia lyase, the enzyme that catalyzes the first step of the pathway, is up-regulated only at t in resistant plants, and this level decreases along insect development. in susceptible plants, pal transcript levels increase at final steps of insect infection, t. however, genes fromlignin and isoflavones synthesis, downstream metabolites, such as cynnamyl alcoholdehydrogenase and isoflavone reductase, are down-regulated in resistant plants. on the other hand, genes from biosynthesis of anthocyanins and tannins, such as flavonoide 3’-hydroxylase and leucoanthocyanidin dioxygenase are up-regulated at all times in resistant plants. this expression profile suggests that synthesis of anthocyanins and tannins is favored over synthesis of flavones. also, activation of upstream genes such as phenylalanine ammonia lyase, chalcone synthase and flavonone dehydrogenase is observed at final stages in susceptible plants, indicating that phenylpropanoid biosynthesis is delayed. several other defense-related genes are also positively regulated in resistant plants, including herbivore-response related genes glutathione-s-transferase and cysteine proteinase inhibitor. apoptosis-related genes have a variable expression profile: catalase is up-regulated throughout insect development but citocrome c oxidase, superoxide dismutase and a senescence-associated protein have no differential expression, and polygalacturonase is up-regulated only at t validation of expression profile for selected candidate-genes we selected genes for validation, listed with corresponding expression levels on table  genes exhibited a consistent expression pattern when quantified by either microarray or qpcr, and the pearson coefficient for this comparison is fold-change value * ** *microarray analysis value. **qrt-pcr analysis value. ***pearson coefficient correlation value. further qpcr analyses were performed to validate expression of selected leaf-miner resistance-associated candidate-genes. these included genes from pathways described above and genes with either no significant hits or similarity to unknown proteins, which may represent coffee specific genes, not yet identified or characterized. the putative caffeine synthase gene exhibited the greatest response to leaf-miner infection, as its expression was significantly higher in control resistant leaves, and also in later stages. the expression of gene sgn-e encoding an unknown protein, is significantly higher in control resistant leaves than in susceptible ones. however, transcript levels dropped in resistant plants at first stages of infection, such as oviposition and egg ecclosion. this gene is a good candidate for differentiation of resistant and susceptible plants, although possibly is not related to resistance. the isocitrate lyase gene is down-regulated in resistant leaves at t at initial steps of infection its transcript levels increased rapidly, but decreased during feeding stages. other evaluated genes also exhibited a regulated expression. the glycerol-3-phosphate acyltransferase gene, a trans-membrane protein associated with synthesis of cutin, is up regulated. a gene encoding metallothionein-like protein, a class of metal-chelators proteins with possible anti-oxidant role, is also up-regulated at initial stages of insect development. genes encoding unknown proteins with lipase protein domain, such as sgn-u sl ch and sgn-u were activated in resistant plants, being up-regulated at oviposition and/or egg ecclosion with rq values of, and respectively. as these genes are down-regulated on leaves at t this profile during insect development suggests a role during defense response. the expression profile of polygalacturonase gene was not confirmed by real-time pcr. according to microarray in silico analysis this gene is activated at t, and afterwards transcript levels drop along insect development. yet by real-time pcr analysis transcript levels are similar in susceptible and resistant plants at all times. this result may reflect differences on genetic background of evaluated plants, once they are part of a population still segregating for some characteristics. however, as this was the only observed discrepancy between all performed analyses, these genetic differences may not be associated with the resistance trait. discussion the use of resistant or tolerant cultivars represents an alternative for reducing the use of chemical defensives, the costs of production, and the negative impact over environment. in coffee, breeding for herbivore resistance is important once coffea arabica, the main commercial coffea species, is susceptible to almost all known coffee pathogens. then resistance must be transferred from other compatible species, which is not always possible due to limited efficiency of inter-specific crosses. at the agronomic institute leaf-miner resistance genes have been transferred from c. racemosa to c. arabica through traditional breeding strategies, and although the program is currently at advanced generations, no resistant commercial cultivar is available yet. the lack of knowledge of molecular aspects controlling the resistance response, and the restricted genetic variability of breeding populations limit the development of genomic-based selection tools. in this context, we aimed to provide information on molecular aspects of leaf-miner defense mechanisms and identification of reliable candidate-marker genes for assisted-selection. those genomic tools associated with traditional breeding strategies guarantee that agronomical traits such as productivity and cup quality will be selected in advanced generations. also, once novel genomic-markers associated with other desirable traits are developed for coffee genotypes, a genome-wide selection strategy will be possible to develop multiple-pathogen resistant cultivars. we chose the microarray analysis approach to assess the co-expression of a large amount of genes, including those that are not looked at in common analyses. nevertheless, the results described in this work indicated that pathways regularly activated in response to herbivory, such as linoleic acid cycle, phenylpropanoids synthesis and apoptosis, are also activated during coffee-response to leaf-miner. genes associated with jasmonate synthesis, such as lipoxygenase and enoyl-coa hydratase, and with flavonoids synthesis, such as chalcone synthase and flavanone 3-hydroxylase-like, are up-regulated in resistant plants. also, pathways from the primary metabolism, such as the citric acid cycle are down-regulated during leaf-miner defense response in resistant plants, a profile observed also in conifers. results of this study provide evidence that most genes encoding enzymes from the citric acid cycle are down-regulated in resistant plants. in a parallel analysis, a metabolite profile was established for resistant and susceptible genotypes using an nmr-based technique and indicated that malate levels on resistant leaves are lower than in susceptible ones. malate results from conversion of either fumarate or glyoxylate. expression of fumarase, that converts fumarate into malate, is down-regulated at t in resistant genotypes, what could explain the low malate levels. production of malate from glyoxylate may also be deficient in resistant plants, once genes encoding for malate synthase, that converts glyoxylate into malate, and for isocitrate liase, the upstream enzyme that converts isocitrate into glyoxylate and succinate, are both down-regulated. in contrast to this profile, susceptible plants exhibit a regular expression levels for these genes at t therefore, both metabolic and transcriptional profiles support the affirmations that citrate cycle is down regulated in leaf-miner resistant coffee plants, and the model of down-regulation of primary metabolism in herbivore-resistant plants. biosynthesis of ja starts with alpha-linoleic acid release in non-injured tissues, triggered by systemin and phospholipase a alpha-linoleic acid is then converted to ja after enzymatic steps performed by 13-lipoxygenase, allene oxide synthase, jasmonate o-methyltransferase and others. several genes from the ja biosynthesis pathway are up-regulated in resistant plants at t including those from downstream steps such as jasmonate o-methyltransferase which expression is 10-fold higher than in susceptible plants. all genes of the ja biosynthesis are either down-regulated or up-regulated at later stages in susceptible plants, as for instance expression of lox increases only at t these observations suggest that the ja signalling pathway, including intermediate signaling-molecules such as oxo-pentenyl-cyclopentane, may be impaired in susceptible plants. down-regulation of genes from later steps of ja biosynthesis, such as allene oxide cyclase, allene oxide synthase, carboxyl methyltransferase, the enzyme that converts jasmonic acid into methyl jasmonate, is observed at t and t in resistant plants. therefore, a feedback regulation may be activated, with a re-programming of transcriptional response upon leaf-miner infection. genes associated with biosynthesis of secondary compounds are shown here to be regulated. expression profile of genes from phenylpropanoids biosynthesis, both up-stream genes such as pal and chs, and downstream genes such as flavonoide 3’-hydroxylase, leucoanthocyanidin dioxygenase, reveals a preferential synthesis of tannins and anthocyanins instead of ligninin, flavones and isoflavones in resistant plants. this profile indicates a direct defense strategy against leaf-miner, once among anthocyanins and tannins are found toxic compounds with antifeedants effects over insects. another gene linked to secondary metabolism is a putative caffeine synthase, which encodes one enzyme from caffeine biosynthetic pathway. caffeine is an alkaloid distributed in coffee plant tissues and organs. the fact that expression of a gene from its biosynthesis is significantly increased upon leaf-miner infection suggests that caffeine may have a role in defense response. however, several studies regarding caffeine and leaf-miner development indicated that this compound has no effect on insect survival rates. once the ultimate goal of this study is to identify potential candidate for markers, several genes were selected for validation using real-time pcr. potential candidates include: isocitrate lyase, which increased expression during initial steps of leaf-miner infection may be co-related with reduction of primary carbon metabolism; putative caffeine synthase, part of an important pathway of coffee plants; glycerol-3-phosphate acyltransferase a gene associated with lipid metabolism and part of cutin biosynthetic pathway, a secondary metabolite; and finally metallothionein-like protein gene, a metal-transporter protein family with an uncertain role in plant metabolism but previously associated with redox responses. future analyses include cloning and re-sequencing genomic regions of target genes from different genotypes in order to identify suitable polymorphisms. among selected genes are those that have no similarity with any known reported gene or protein. although they could not yet be associated with a biological process, their expression profile was very specific and related to defense response. for instance, genes sgn-e sgn-e sgn-e were activated upon leaf-miner infection in resistant plants, and therefore represent good candidates for further investigation. another interesting unknown gene is sgn-e which is highly activated in resistant plants at t but is repressed upon infection. the expression profile indicates that this gene is useful for early differentiation between resistant and susceptible plants. in summary, differential expression profiles between resistant and susceptible genotypes are observed even in the absence of leaf-miner, indicating that defense is already build up in resistant plants, as a priming mechanism. then, a systemic defense response may be more rapidly activated in resistant plants, once basic compounds such as nitrogen and sugars are readily available as a result of repression of primary metabolism. this shift in plant metabolism is common after a pathogen attack, where defense-related pathways are activated, resulting in reduction of growth and reproduction, and in changes on link-source relationship. during herbivore-defense transcript levels of genes involved in photosynthesis are also down-regulated, probably as a strategy to liberate nitrogen compounds for the secondary metabolism. maintenance of these physiological and metabolic states has a high energetic cost, and could represent a survival limitation if nutritional conditions on the field are depleted. actually, field observations in cultivated areas demonstrated that leaf-miner resistant coffee plants, under a severe nutritional deficit, are attacked by the leaf-miner at the same intensity as susceptible plants.  as a basal defense state is decisive for triggering a rapid resistance response, genes associated with priming validated here, represent key genes for assisted-selection. future studies will focus on comparisons of selected genes genomic sequences, from both resistant and susceptible parental lines, to identify suitable marker polymorphisms. 
a total of est sequences from the coffee genome database were selected for a microarray analysis. up-regulated genes include those responsible for defense mechanisms, hypersensitive response and genes involved with cellular function and maintenance. the microarray analysis indicated that resistant plants differ in gene expression profile. expression changes were found in resistant plants for genes. also, our analyses indicated that differential expression profiles between resistant and susceptible genotypes are observed in the absence of leaf-miner, indicating that defense is already build up in resistant plants, as a priming mechanism. to understand the transcriptional basis of this response, we compare gene expression profiles of two coffee genotypes, susceptible and resistant to leaf miner. plants have evolved complex mechanisms to recognize such attacks, and to trigger a defense response. our results show evidences that coffee defense responses against leaf-miner attack are balanced with other cellular functions. fluorescence probes were synthesized using mrna from the infested and non-infested coffee plants.
 a successful development of herbivorous insects into plant tissues depends on coordination of metabolic processes. plants have evolved complex mechanisms to recognize such attacks, and to trigger a defense response. to understand the transcriptional basis of this response, we compare gene expression profiles of two coffee genotypes, susceptible and resistant to leaf miner. a total of est sequences from the coffee genome database were selected for a microarray analysis. fluorescence probes were synthesized using mrna from the infested and non-infested coffee plants. array hybridization, scanning and data normalization were performed using nimble scan® e arraystar® platforms. genes with foldchange values +/- were considered differentially expressed. a validation of differentially expressed genes was performed in infected plants using qrt-pcr approach.  the microarray analysis indicated that resistant plants differ in gene expression profile. we identified relevant transcriptional changes in defense strategies before insect attack. expression changes were found in resistant plants for genes. up-regulated genes include those responsible for defense mechanisms, hypersensitive response and genes involved with cellular function and maintenance. also, our analyses indicated that differential expression profiles between resistant and susceptible genotypes are observed in the absence of leaf-miner, indicating that defense is already build up in resistant plants, as a priming mechanism. validation of selected genes pointed to four selected genes as suitable candidates for markers in assisted-selection of novel cultivars.  our results show evidences that coffee defense responses against leaf-miner attack are balanced with other cellular functions. also analyses suggest a major metabolic reconfiguration that highlights the complexity of this response. coffea arabicaleaf minermicroarrayplant defense 