Transcriptomic Analysis of Avocado Hass (Persea americana Mill) in the Interaction System Fruit-Chitosan-Colletotrichum

Avocado (Persea americana) is one of the most important crops in Mexico as it is the main producer, consumer, and exporter of avocado fruit in the world. However, successful avocado commercialization is often reduced by large postharvest losses due to Colletotrichum sp., the causal agent of anthracnose. Chitosan is known to have a direct antifungal effect and acts also as an elicitor capable of stimulating a defense response in plants. However, there is little information regarding the genes that are either activated or repressed in fruits treated with chitosan. The aim of this study was to identify by RNA-seq the genes differentially regulated by the action of low molecular weight chitosan in the avocado-chitosan-Colletotrichum interaction system. The samples for RNA-seq were obtained from fruits treated with chitosan, fruits inoculated with Colletotrichum and fruits both treated with chitosan and inoculated with the fungus. Non-treated and non-inoculated fruits were also analyzed. Expression profiles showed that in short times, the fruit-chitosan system presented a greater number of differentially expressed genes, compared to the fruit-pathogen system. Gene Ontology analysis of differentially expressed genes showed a large number of metabolic processes regulated by chitosan, including those preventing the spread of Colletotrichum. It was also found that there is a high correlation between the expression of genes in silico and qPCR of several genes involved in different metabolic pathways.


INTRODUCTION
The avocado (Persea americana) cv Hass, is from the economic point of view, one of the most important fruits in Mexico because it is the main producer, exporter and consumer of this fruit in the world. However post-harvest diseases represent a serious problem for exports of avocado fruit. Anthracnose, caused by Colletotrichum sp., is one of the most economically important diseases of the avocado fruit during storage and marketing, with losses close to 20% of the total avocado production (Freeman et al., 1996;Rodríguez-López et al., 2009). To control this pathogen, chemical fungicides are normally used, even though they have the potential of affecting the environment and the health of consumers (Brent and Hollomon, 2007). For this reason it is important to seek biological alternatives to the use of fungicides, such as chitosan, which has no polluting characteristics and acts as an elicitor of plant defense responses. Chitosan, or poly [ß-(1-4) -2-amino-2-deoxy-D-glucopyranose], is a copolymer of N-acetyl-D-glucosamine units derivative from chitin by a deacetylation treatment in alkaline medium (Islam et al., 2011). The elicitor activity of chitosan is well documented, for instance, it has been reported that upon external application of chitosan the production of chitinases, protease inhibitors and phytoalexins can be induced, triggering a defense response in plants (Terry and Joyce, 2004;Amborabé et al., 2008). However, the mechanism by which chitosan elicits a defense response in plants is not fully elucidated, although some reports indicate that secondary metabolites could have an important role in this response. For example, chitosan application in pear (Pyrus pyrifolia L. cv. Xuehua) resulted in an increase of phenylalanine ammonia lyase (PAL), polyphenol oxidase (PPO) and peroxidase (POD) activities, as well as in the induction of the expression of β-1,3-glucanase and chitinase genes, which may also be involved in defense against Alternaria kikuchiana and Physalospora piricola (Meng et al., 2010). Similarly, using a chitosan film in tomato (Lycopersicon esculentum Mill) it was observed an increase in the PPO and POD activities, and in the production of phenolic compounds that correlated with an increase of resistance against Botrytis cinerea and Penicillium expansum (Liu et al., 2007). Moreover, in Arabidopsis seedlings chitosan regulated the expression of defense response genes against Botrytis cinerea and the genes involved in biosynthesis of camalexin (Povero et al., 2011). The activity of enzymes like PAL and tyrosine ammonialyase (TAL) and the expression of genes like those encoding for PPO, POD and Isoflavone Synthase I and II, have been reported to increase in the presence of chitosan in mango fruit (Khan et al., 2003;Berumen-Varela et al., 2015) and soybean (Chen et al., 2009). More recently, RNA-seq analysis has allowed the identification of defense genes induced in orange leaves upon treatment with chitosan, such as different transcription factors and genes involved in some hormone pathways (Coqueiro et al., 2015). In spite of these previous studies, the mechanism by which the chitosan is perceived by the plant and how it induces resistance remains unknown, therefore it is important to design studies that could help shed some light on these processes.
With the aim of having a better understanding of the mechanisms of action of chitosan in the induction of resistance to Colletotrichum sp. in avocado fruit, here we analyzed the global transcriptional profile of avocado fruits treated with chitosan and in the presence or absence of Colletotrichum. We also report on the differentially expressed genes in response to chitosan that could be involved in the resistance induced by this elicitor.

Isolation and Identification of Colletotrichum sp. from Avocado Fruits
In order to isolate and identify the pathogen causing anthracnose symptoms in avocado, we used fruits that were harvested at physiological maturity stage from an orchard located in Tepic, Mexico. Fungus development was encouraged by placing fruits in chambers at a relative humidity of 90-95% and 25 • C for 5 days. Once disease symptoms were developed in the fruit, damaged tissue was washed by immersion in sodium hypochlorite solution (2%) for 2 min, rinsed with sterile water for 2 min and placed on filter paper to eliminate moisture. Moisture-free tissue sections were placed on potato dextrose agar (PDA) (DIBCO) and incubated at 25 • C for 24-72 h. Fungus genus was determined according to the taxonomic keys described by Barnett and Hunter (1998) and by molecular techniques using amplification of ITS1 and ITS4 primers according to White et al. (1990).

Fruit Maturity and Sampling Time
To establish the sampling times for the RNA-seq analysis, we first evaluated the development of anthracnose symptoms caused by Colletotrichum gloeosporioides in avocado fruits at different states of maturity: physiological, intermediate and consumer maturity. For this, 10 fruits (for each state of maturity) were inoculated with 4 µL spore suspension (1 × 10 6 mL −1 ) of C. gloeosporioides previously isolated. Development of the infection was monitored in a stereoscope to follow the fruit invasion by the fungus at 0, 1, 2, 6, 9, 24, 72, and 96 h postinoculation. For RNA-seq experiments we then selected the maturity state "intermediate" and sampling times at 0, 1, 6, 9, and 24 h post-inoculation. In order to study the elicitor effect of chitosan, another group of 10 fruits in a state of "intermediate" maturity were inoculated with 40 µL of the spore suspension and treated with a solution of low molecular weight chitosan (LMWC) 1.5% (w/v) (Sigma Aldrich; viscosity 35 cps in 1% chitosan solution; 96.1% deacetylation), the inoculated/treated fruits were kept at 25 • C for 7 days to record the disease development.

Chitosan Treatment
After determining the appropriate sampling conditions for RNAseq, four treatments were performed: untreated control fruits (C), fruits inoculated with C. gloeosporioides (P), fruits treated with LMWC without pathogen (Q) and fruits inoculated with C. gloeosporioides and treated with LMWC (QP). Pathogen inoculation (P and QP treatments) was performed using an insulin syringe as follows: 40 µL of a spore suspension (1 × 10 6 spores/mL) were inoculated by penetrating 3 mm into the fruit peel, C and Q treatments were inoculated with 40 µL of sterile water. After 30 min, fruits were immersed during 1 min into LMWC solution at 1.5% w/v, the treatments that did not contain chitosan (C and P) were immersed in sterile water also for 1 min. Sample collection was performed by cutting sections of 5 × 5 cm of peel and pulp of the avocado fruits, which were immediately frozen in liquid nitrogen and macerated to a fine powder in a sterile mortar. Samples of three biological replicates of each treatment were stored independently as powder at −80 • C for later use. For all these treatments samples were taken at: initial time (0 h), early response (1, 6, and 9 h) and late response (24 h). Samples 6 and 9 h were mixed, because no differences in the development of the disease were observed.

RNA Extraction and Sequencing
RNA extraction was performed for each of the biological replicates according to the methodology proposed by Djami-Tchatchou and Straker (2012). Then samples from the same treatment were mixed equimolarly. cDNA library preparation and sequencing was carried out at the National Laboratory of Genomics for Biodiversity (LANGEBIO-CINVESTAV-Irapuato, Mexico). Illumina HiSeq 2000 system was used and the sequencing run was 2 × 100 (paired-end). Raw data were obtained as FASTQ files, which include reads and their qualities defined based on PHRED, with an average quality >30. Raw sequencing data is available at NCBI under accession number SUB2518488.

Sequence Mapping
A reference transcriptome of avocado (Persea americana Mill) var drymifolia (TAD) (Ibarra-Laclette et al., 2015) and also a reference transcriptome of avocado fruit (Persea  americana Mill) cv Hass (TFH)  were used. The Bowtie2 program was used for mapping (Langmead and Salzberg, 2012), this tool also allows us to obtain the read count for any particular gene model.

Counting and determination of FPKM (Fragments Per
Kilobase of transcript per Million mapped reads) was performed using the program eXpress (Roberts and Pachter, 2013).

Statistical Analysis
Differentially expressed genes were identified by maximum likelihood analysis (P < 0.05) using the R statistical (Foundation for Statistical Computing, 2013) with the EdgeR library (McCarthy et al., 2012). The matrix was compared using information from the control samples (fruit not treated with chitosan and not inoculated at their respective sampling times).

Functional Annotation
To identify enriched GO terms (Process; Up-regulated and Down-regulated genes) Plant MetGenMAP (Ontology Process, multi-test correction FDR, p-value < 0.05) (Joung et al., 2009) as well as Classification SuperViewer (Provart and Tong, 2003) were used. For this analysis we used genes showing highest homology to A. thaliana genes, eliminating those locus that were repeated. Likewise, the MapMan program (Usadel et al., 2009) was used to identify genes that are specific of metabolic pathways of interest.

Gene Expression Validation
To validate the bioinformatics analysis on gene expression prediction based on the RNA-seq data, we used real-time RT-PCR for quantification of several avocado transcripts. For this, aliquots were taken from the same RNA that was used for sequencing, cDNAs were synthesized from 2 µg of total RNA using the Superscript II (Invitrogen), first strand synthesis procedure was followed by RNase H digestion (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. Primers for different genes encoding transcription factors, plant disease resistance proteins, genes involved in ethylene biosynthesis and secondary metabolites, as well as a gene as internal control (Table S1) were used. The primers were designed using the Primer Express 2.0 software (Applied Biosystems, Foster City, CA, USA) with melting temperatures (Tm) of 58-60 • C, primer lengths of 20-24 bp, and amplicon lengths of 61-150 bp. PCR reactions were performed in 48-well plates in a StepOne TM Real-Time PCR System (Applied Biosystems) using SYBR R Green. Three technical replicates were performed for each evaluation. Reactions were done in 25 µl containing 200 nM of each primer, 5 µl cDNA (corresponding to ∼3 ng), and 12.5 µl 2 × SYBR Green Master Mix Reagent (Applied Biosystems). PCR parameters were as recommended by the manufacturer: 95 • C for 10 min, 40 cycles of 95 • C for 15 s, and 60 • C for 1 min. Nontemplate controls were included for each primer pair and each PCR reaction was completed in triplicate. Dissociation curves for each amplicon were then analyzed to verify the specificity of each amplification reaction; the dissociation curve was obtained by heating the amplicon from 60 to 95 • C.

Isolation and Identification of the Pathogen Causing Anthracnose
Colletotrichum gloeosporioides is the causal agent of anthracnose in avocado fruits, whose characteristic is to remain dormant until the fruit ripens to the maturity consumer stage, when the defense system of the fruit is less active, allowing an easier and more effective infection (Rodríguez-López et al., 2009). Symptoms of anthracnose in avocado fruits were characterized as black sunken lesions with salmon-colored dots in the center and white mycelial growth on the surface. Pathogen identification at species level was performed by isolating the fungus from diseased avocado and grown in culture media. Morphological characterization was performed by evaluating mycelial growth at random, white mycelium, cottony little consistency, grayish colony, short hyphae on the edge of growth and yellowing in the bottom of the petri dish. Microscopic analysis showed that the isolate, according to taxonomic keys, belongs to the genus Colletotrichum (hyaline conidia, form straight, cylindrical and obtuse at the apex, measuring 9-24 um in length and 3-4.5 um width), while the amplification of ITS1 and ITS4 and subsequent fragment sequencing showed a 98% homology to C. gloeosporioides (File S2).

Determination of Ripeness of the Fruit for RNA-Seq
To determine the ripeness of the fruit for later transcriptomic analyses, a scale of maturity was established according to the fruit color change from the first postharvest day to maturity consumption stage (Figure 1). During the first 5 days the fruit is in physiological maturity, between 6 and 9 days we considered it as "intermediate maturity, " and after day 10 it is in consumption maturity (stored at 25 • C). In order to determine the sampling times for RNA extraction and sequencing, samples were collected from C. gloeosporioides inoculated avocado fruits (without chitosan treatment) at different stages of maturity (physiological, intermediate and consumption) and analyzed for disease symptoms (Figure 2). Inoculated fruits at the physiological maturity stage showed no clear symptoms of pathogen infection, possibly because some antifungal compounds present in the fruit during this stage of maturity prevent fungal infection (Prusky and Lichter, 2007). In fruits at intermediate maturity disease symptoms started to be visible after 9 h of inoculation and at 96 h postinoculation a limited area of the fruit showed clear disease symptoms. In contrast, in the fruits at the consumption maturity, disease spreads rapidly, becoming visible as early as 2 h postinoculation, which became widely spread after 96 h postinoculation. Since the development of infection in fruits at the intermediate maturity stage was slower and allowed to detect more concise changes during the infection process, we decided to use the intermediate fruit ripening stage to investigate the global transcriptional changes induced by chitosan treatment (Q condition), inoculation by C. gloeosporioides (P condition) and both, chitosan treatment and pathogen inoculation in the avocado fruits (QP condition).

Disease Reduction in Avocado Fruits Treated with Chitosan
Most control fruits developed anthracnose symptoms after 7 days, such as fleshiness and softening pulp (Figure 3A), indicating that C. gloeosporioides was present in its latent state. When fruits were inoculated with 40 mL of a suspension of 1 × 10 6 spores/mL, all of them developed anthracnose and it was observed the growth and invasion of the pathogen covering most of the fruit, showing more drastic symptoms than control fruits ( Figure 3B). In contrast, only 10% of fruits inoculated with C. gloeosporioides and treated with LMWC showed deterioration in quality and the pathogen failed to spread, maintaining a healthy fruit (Figure 3C). These results indicate that LMWC treatment readily increases the avocado fruit resistance against C. gloeosporioides.

RNA-Seq and Mapping to the Reference Transcriptomes
To obtain and analyze the global transcript profile of the avocado fruit upon treatment/inoculation with chitosan/C. gloeosporioides, RNA samples were obtained from fruits treated with chitosan (Q condition), inoculated with the fungus (P condition), chitosan-treated and inoculated with the fungus (QP condition) and a control (non-treated and non-inoculated). The RNA from each condition at different set times, was used for cDNA library preparation and sequenced using the ILLUMINA platform HiSeq 2000. Over 300 million sequences were obtained with an average length of 100 bp in paired-end format. Table 1 shows the number of paired-end reads obtained for each of the treatments. Read alignment to the TAD (Transcriptome Avocado var Drymifolia), resulted in an average of 93% mapped reads, while alignment to the TFH (Transcriptome Fruit var Hass) resulted in an average of 83% sequencing reads mapped. However, mapping to the TFH resulted in an increase of reads aligned concordantly exactly 1 time (ACE = 1) (an average of 62%) whereas mapping to the TAD only an average of 18% of reads ACE = 1 were obtained ( Table 1). Reads aligned concordantly exactly 1 time were used to avoid redundancy in further expression analysis (Langmead and Salzberg, 2012). These mapping differences are possibly due to the different characteristics of the two reference transcriptomes used here, TAD was derived from several organs of the drymifolia avocado plant (including fruit) while TFH was generated from the Hass avocado fruit only (Ibarra-Laclette et al., 2015;Kilaru et al., 2015).

Differential Gene Expression Profiles and Gene Ontology Categorization
The differential gene expression profiles were obtained to identify those genes in the avocado fruit that were up-or downregulated during treatment with chitosan or/and inoculation with the pathogen. Each of the sequenced libraries showed a similar pattern of the expression profiles using the two reference transcript datasets (TAD and TFH) ( Table 2). The results show that, when compared with control fruits (nontreated and non-inoculated), at 1 h in the condition P (inoculated vs. control) and Q (treated with LMWC vs. control) the highest amount of unigenes with differential expression is obtained, whereas in QP (inoculated and treated with LMWC vs. control) the highest amount is observed at 24 h. In contrast, the minimum amount of differentially expressed genes was obtained in conditions P_0 and P_6/9, whereas in conditions Q and QP the minimum amount was at time 6/9. Considering all differentially expressed genes, up-regulated and downregulated, using the two reference transcriptomes, we found that at the QP condition a higher differential expression was obtained compared to the control fruit, while in condition P a smaller amount of differentially expressed unigenes was obtained.
Gene expression analysis at different sampling times shows that at 24 h most genes of the QP condition were induced (FC > +2) (Figure 4), whereas in the P condition were mostly repressed and in the Q condition the number of genes induced and repressed is similar. More than 2400 genes were differentially expressed in the QP_24 condition, from this set of genes 2000 were induced (Figure 5). Interestingly, the highest number of differential gene expressed was observed in Q condition at initial times (0 and 1 h), in contrast we observed that in the same Q condition there was a minor number of differential expressed genes at long time (24 h) in comparison with P and QP conditions. This suggests that the fruit induced its defense system by regulating several genes in short times just after treatment with chitosan.
To identify and assign enriched Gene Ontology (GO) terms to our dataset (Plant MetGenMap, Ontology Process, multi-test correction FDR, p_value < 0.05), all of the differentially upand down-regulated unigenes showing sequence homology to A. thaliana (DNA sequences defined in terms by at least 60% similarity) were functionally categorized in biological process (BP), cellular component (CC), and molecular function (MF) ( Table S3). The gene expression profile from the different conditions, that were analyzed by GO categorization showed important differences. In the P_24 condition some downregulated unigenes were associated with "Systemic acquired resistance (SAR)" and "incompatible interaction with fungus, " suggesting that a decrease in the expression of these genes could be the cause of the observed disease development in the fruit. Other unigenes that are up-regulated in the same condition (P) are associated with "Defense response to fungus, " but such defense is not effective since the fruit is susceptible in this condition. In contrast, when only chitosan is present (condition Q_24) some of the unigenes that are categorized in biotic stimulus, specifically those involved in SAR response, are up-regulated. The QP_24 condition shows a possible combined response between the processes occurring in both Q and P separately, for example, the unigenes UN26087 (Log 2 FC = −4.37 and Log 2 FC = −3.17, P and QP conditions, respectively) and UN17145 (Log 2 FC = 3.09 and Log 2 FC = 5.30, Q and QP conditions, respectively), are both associated with SAR. Interestingly, in the QP condition a considerable increase of the expression was observed in the genes that are classified in this category (Table 3). GO analysis of the QP condition at 24 h showed processes as response to stimulus, response to stress, response to chemical stimulus and response to abiotic stimulus as the main processes that are affected due to treatment with chitosan and pathogen inoculation ( Table 4, and  Table S4).

Metabolic Pathways Related to Biotic Stress
In order to analyze changes in metabolic pathways related to biotic stress in the fruit for each of the treatments, we used the MapMan software (Usadel et al., 2009) based on the homology to A. thaliana genes of known function. Several subsets of genes encoding receptor-like kinases, signaling proteins, transcriptional factors, oxidative stress response elements, secondary metabolism factors, phytohormone-responsive genes and defense-related genes were identified. Notably, a greater number of genes were differentially expressed in the QP condition in all of the identified metabolic pathways. As described before, 24 h after inoculation fruits without chitosan showed characteristic symptoms of infection caused by C. gloeosporioides (Figure 2), whereas in chitosan treated fruits infection was not visible and a large number of genes related response to different metabolic processes were induced, mainly genes related to biotic stress responses (Figure 6).
An increase in the expression of unigenes actively involved in metabolic pathways related to biotic stress is observed from early times in fruits treated with LMW chitosan (Q_0) ( Table 5). GO classification shows the relationship of some induced unigenes related to SAR and those involved in the salicylic acid and jasmonic acid signaling pathways. Also, a slight induction of genes involved in response to both biotic and abiotic stress is observed, suggesting that a priming state could be acting where the fruit is conditioned for the super-activation of defenses activated by pathogen inoculation.

Gene Expression Validation by qRT-PCR
To validate and confirm the gene expression results obtained by bioinformatics analysis of our RNA-seq data, we performed real time quantitative RT-PCR of several avocado unigenes that showed different expression patterns ( Table 6). For this, we considered unigenes that mapped to both reference transcriptomes and that corresponded to various metabolic pathways such as ethylene, secondary metabolites, transcription factors and pathogenesis related. The results show a trend of expression as well as a high correlation between the expression results by RNA-seq and qRT-PCR.

RNA Sequencing and Differential Expression Profile
The use of two reference transcriptomes (TDH and TFH) allowed us to have a greater coverage in the analysis of our sequencing data. The comparative analysis of our data with both reference transcriptomes resulted in different percentage of mapped sequences (Table 1). On average 50% less of differentially expressed unigenes (induced and repressed) were obtained in the P condition when compared with those obtained in Q and QP conditions. On the other hand, the little difference of genes differentially expressed in the P condition in comparison with the control, could be explained by the fact that the control fruits are frequently naturally infected with the phytopathogen C. gloeosporioides, which is commonly found in the host from early stages of fruit development (Beno-Moualem and Prusky, 2000;Mendgen and Hahn, 2002). Thus, the development of the infection can be observed as the fruit ripens. However, in our experimental setting, we observed that at 1 h post inoculation, the fruit is highly resistant probably due to the presence of antifungal compounds ( Table 2) (Yakoby et al., 2002;Wharton and Diéguez-Uribeondo, 2004;Prusky and Lichter, 2007). By contrast, in the presence of chitosan and the pathogen (QP), there is a strong defense response even at 24 h post-inoculation ( Table 2). We hypothesize that the observed differences between the expression profile in QP condition in comparison with Q condition, are possibly due to potentiation of cellular signaling processes, as a consequence of the increased stress in the fruit induced by the chitosan-pathogen-fruit interaction.
The major GO categories for the P condition were the "response to stress to abiotic or biotic stimulus" and "metabolic processes, " whereas in the Q condition there were "cellular processes" and "response to stress to abiotic or biotic stimulus." Interestingly, when the fruit was inoculated with the pathogen and treated with chitosan (QP), a significant increase of differentially expressed genes that were located in several GO categories was observed, for example, "cellular processes, " "metabolic processes, " "response to stress abiotic or biotic stimulus, " "biological processes, " "transport, " "cell organization and biogenesis" and "signal transduction, " this indicate that the chitosan could be able to activate different metabolic pathways involved in the defense of the fruit.
In order to have specific annotation data from the profile expression obtained here, we used the information from A. thaliana based on sequence homology. Several avocado unigenes that were involved in SAR and that were induced in the P, Q, and QP conditions, have homologs in A. thaliana genes involved in the response to different pathogens. For example, AT3G02550 (leucine-rich repeat receptor-like kinase NILR1) and AT1G74360 (nematode-induced LRR-RLK 1) are required for induction of innate immunity to parasitic nematode (Mendy et al., 2017); AT5G17760 (P-loop containing nucleoside triphosphate hydrolases superfamily) is related to plant response during geminivirus infection (Ascencio-Ibáñez et al., 2008); AT1G15520 (ABCG40 atp-binding cassette g40) is related to Phytophthora resistance (Wang et al., 2015); and AT3G02260 (Calossin-like protein required for polar auxin transport) is involved in immune response to H. arabidopsidis (Meteignier et al., 2017) (Table 3). This suggests that the chitosan and/or phytopathogen perception could converge in some signal points activating common metabolic pathways in the fruit. Although different defense responses may be activated depending on the tissue or organ attacked by a pathogen, some common pathways  , 2008). Activation of these processes in short times possibly induces the "state of priming" related to sensitization of a cell or organism for enhanced defense; this condition causes faster and more robust activation of defense responses upon challenge with a pathogen (Conrath et al., 2015). Our results suggest that chitosan could induce a state of priming in short times after application, which promotes effective fruit resistance against C. gloeosporioides (Table 5).
On the other hand, several studies (Wang et al., 2010;Rodríguez-Carpena et al., 2011;Widsten et al., 2014) indicate that there are compounds in the epicarp of avocado fruits, mainly phenolic compounds such as catechins, epicatechin, proanthocyanidins, and quercetin, which have antibacterial and antifungal activity, and a high antioxidant capacity. It has also been identified the diene AFD (1-acetoxy-2-hidroxiy-4-oxoheneicosa 12, 16 diene) (Ardi et al., 1998) with high antifungal capacity. Our analysis showed that treatment of fruits with chitosan downregulates some genes involved in biosynthesis of phenylpropanoids (CHS and FLS) ( Table 6) whereas the gene encoding to 4-Coumarate coenzyme A ligase (4CL) is up-regulated. Likewise, we found transcription factors such as WRKY22 (Povero et al., 2011) andERF (Oñate-Sánchez andSingh, 2002;Vallejo-Reyna et al., 2015) that were induced in different treatments with chitosan compared to control treatments in avocado fruits. According with our results, the genes involved in the biosynthesis of diene AFD (Avfad1 and Avfae12-3) are up-regulated in the fruits treated with chitosan ( Table 6), suggesting a direct association between the diene AFD synthesis and chitosan.
The results presented in this study shows that chitosan acts as a molecule able to induce multiple metabolic responses in avocado fruit that collectively implements a defense system capable of counteract the infection by C. gloeosporioides. However, further studies are needed to experimentally determine the role and function of up-and down-regulated genes in fruits treated with chitosan and to dissect their participation in the resistance to pathogens as C. gloesporioides in avocado.

AUTHOR CONTRIBUTIONS
AC-L provided the idea of the work. AC-L and LX designed the experiments. PG and UL contributed with the in vivo assays in fruit. EC performed the identification of the phytopathogen. LX, SG, and JV performed the bioinformatics analysis. LH contributed in the sequencing and generation of libraries. LX and JV performed qRT-PCR assays. LX, AC-L, JV, and LH participated in the interpretation of results and critically reviewed the manuscript. LX wrote the paper. All authors read and approved the final manuscript.

ACKNOWLEDGMENTS
Authors express their gratitude to the Mexican National Council for Science and Technology (CONACYT) for financing this work through the Basic Science Research Project #252423. Partial funding was from DGAPA-UNAM: IA203313-IN214917, and CONCYTEG.