Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Transcriptional Analyses of Mandarins Seriously Infected by ‘Candidatus Liberibacter asiaticus’

  • Meirong Xu ,

    Contributed equally to this work with: Meirong Xu, Ya Li

    Affiliation Citrus Huanglongbing Research Laboratory, South China Agricultural University, Guangzhou, Guangdong, China

    ⨯
  • Ya Li ,

    Contributed equally to this work with: Meirong Xu, Ya Li

    Current Address: Agricultural College, Guangdong Ocean University, Zhanjiang, Guangdong, China

    Affiliation Citrus Huanglongbing Research Laboratory, South China Agricultural University, Guangzhou, Guangdong, China

    ⨯
  • Zheng Zheng,

    Affiliation Citrus Huanglongbing Research Laboratory, South China Agricultural University, Guangzhou, Guangdong, China

    ⨯
  • Zehan Dai,

    Affiliation Citrus Huanglongbing Research Laboratory, South China Agricultural University, Guangzhou, Guangdong, China

    ⨯
  • Yang Tao,

    Affiliation Citrus Huanglongbing Research Laboratory, South China Agricultural University, Guangzhou, Guangdong, China

    ⨯
  • Xiaoling Deng

    xldeng@scau.edu.cn

    Affiliation Citrus Huanglongbing Research Laboratory, South China Agricultural University, Guangzhou, Guangdong, China

    ⨯

Abstract

A range of leaf symptoms, including blotchy mottle, yellowing, and small, upright leaves with a variety of chlorotic patterns resembling those induced by zinc deficiencies, are associated with huanglongbing (HLB, yellow shoot disease), a worldwide destructive citrus disease. HLB is presumably caused by the phloem-limited fastidious prokaryotic α-proteobacterium ‘Candidatus Liberibacter spp.’ Previous studies focused on the proteome and transcriptome analyses of citrus 5 to 35 weeks after ‘Ca. L. spp.’ inoculation. In this study, gene expression profiles were analyzed from mandarin Citrus reticulate Blanco cv. jiaogan leaves after a 2 year infection with ‘Ca. L. asiaticus’. The Affymetrix microarray analysis explored 2,017 differentially expressed genes. Of the 1,364 genes had known functions, 938 (46.5%) were up-regulated. Genes related to photosynthesis, carbohydrate metabolic, and structure were mostly down-regulated, with rates of 92.7%, 61.0%, and 80.2%, respectively. Genes associated with oxidation-reduction and transport were mostly up-regulated with the rates of 75.0% and 64.6%, respectively. Our data analyses implied that the infection of ‘Ca. L. asiaticus’ could alter hormone crosstalk, inducing the jasmine acid pathway and depressing the ethylene and salicylic acid pathways in the citrus host. This study provides an enhanced insight into the host response of citrus to ‘Ca. L. asiaticus’ infection at a two-years infection stage.

Citation: Xu M, Li Y, Zheng Z, Dai Z, Tao Y, Deng X (2015) Transcriptional Analyses of Mandarins Seriously Infected by ‘Candidatus Liberibacter asiaticus’. PLoS ONE 10(7): e0133652. https://doi.org/10.1371/journal.pone.0133652

Editor: Ji-Hong Liu, Key Laboratory of Horticultural Plant Biology (MOE), CHINA

Received: April 9, 2015; Accepted: June 29, 2015; Published: July 21, 2015

Copyright: © 2015 Xu et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All CHP files or a matrix table are available from the GEO database (accession number GSE67376, http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE67376).

Funding: The experiment in this study was funded by National Natural Science Fund (Grant No.31201480) from National Natural Science Foundation of China and Modern Agricultural Industry Technology System (Grant No.C09058) from Chinese Ministry of Agriculture. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Citrus in tropical and sub-tropical areas is subjected to Huanglongbing (HLB, Yellow shoot disease), a highly descructive bacterial disease. Whole genome databases of two citrus varieties have been constructed (http://citrus.pw.usda.gov/, http://www.citrusgenomedb.org and http://citrus.hzau.edu.cn/orange/download/index.php) [1, 2]. Substantial research efforts have been made to study the citrus-’Candidatus Liberibacter spp.’ interaction using microarrays based on expression sequence tags, hi-seq based on transcriptome and comparative proteomic approaches [3–13] (Table 1).

thumbnail
Download:
Table 1. A list of citrus genomic, transcriptomic or proteomic studies.

https://doi.org/10.1371/journal.pone.0133652.t001

The use of genomic, transcriptomic, and proteomic approaches in the study of citrus responses to ‘Ca. L. spp.’ infections has improved our knowledge on HLB. The microarray method has been used to study the gene expression changes in different citrus leaves or fruits. Differentially expressed genes (DEGs) of ‘Ca. L. asiaticus infected ‘Valencia’ oranges involved in cell defense, transport, cellular organization, photosynthesis, and carbohydrate metabolism were most notable, especially the genes coding the phloem-specific lectin PP2-like protein and starch biosynthetic genes [3]. The isobaric tags for relative and absolute quantitation (iTRAQ) technique was used to characterize proteome changes in ‘Madam Vinous’ after ‘Ca. L. asiaticus’ inoculation, identifying potential targets of early infections, such as four miraculin-like proteins, chitinase, Cu/Zn superoxide dismutase and lipoxygenase [4]. The RNA-seq method has been used to characterize the gene regulatory network in ‘Valencia’ sweet orange’s response to ‘Ca. L. asiaticus’[7, 13]. The pathogen-induced accumulation of transcripts for phloem-specific lectin PP2-like protein, ADP-glucose pyrophosphorylase, starch synthase, granule-bound starch synthase, starch debranching enzyme, miraculin-like proteins, Cu/Zn superoxide dismutase, 2-oxoglutarate and Fe(II)-dependent oxygenases, WRKY transcription factors, sugar transporters, chitinases, lectin-related proteins, peroxiredoxins and a CAP 160 protein, were notable. The ‘Ca. L. asiaticus’ infection triggered a decrease in the production of photosynthesis-related proteins, defense-related pathogen-response proteins, oxygen-evolving enhancer proteins, a PSI 9 kDa protein and a Btf3-like protein were notable.

Most studies mentioned in Table 1 used leaf materials, including symptomatic (SY), asymptomatic (AS), apparently healthy (AH), healthy (H), immature leaves (IL) and mature leaves (ML). RNA samples of fruit, including flavedo, peel, vascular tissue, or juice vesicles, roots and stems were also extracted for the high-throughput analysis of citrus responding to ‘Ca. L. asiaticus’ infection (Table 1). The leaf, stem and root samples from the individual studies were collected within one year of ‘Ca. L. asiaticus’ infection. Moreover, fruit samples were collected from trees with an ambiguous infection time, probably more than 12 months after infection. In the citrus industry, the growers generally keep the diseased trees until almost no economic value is output, usually more than 2 years.

The objective of this study was to identify, the groups of citrus genes involved in the late stage of ‘Ca. L. asiaticus’ infections to provide information for HLB management.

Materials and Methods

Plant materials

These experiments were performed using Citrus reticulate Blanco cv. jiaogan in 2012. The Affymetrix GeneChip citrus genome were applied to study the molecular pathways mediated by ‘Ca. L. asiaticus’ inoculated 3-year-old jiaogan seedlings. Ten jiaogan trees scioned on Cleopatra mandarin (Citrus reticulata Blanco) rootstocks kept in screenhouse, under natural light conditions were selected. Each of them was graft-inoculated with one sweet orange scion with or without ‘Ca. L. asiaticus’ in Dectember, 2009. Two years later, PCR detection of ‘Ca. L. asiaticus’ using primer pairs OI1/OI2c [14] identified 5 positive trees with typical blotchy mottling symptoms, while the other 5 with healthy leaves were confirmed negative for Ca. L. asiaticus. Of them, three trees showing similar HLB symptoms and similar size were chosen as materials (Fig 1A). Three healthy trees without ‘Ca. L. asiaticus’ infection with similar growth conditions were selected as controls (Fig 1B). Fully expanded mature leaves with blotchy symptoms and leaves at the similar portions on healthy trees were collected for RNA extraction at 25 month after inoculation.

thumbnail
Download:
Fig 1. Phenotype symptoms of Citrus reticulata Blanco cv. Jiaogan 25 month-after-inoculation.

A, One of the three trees grafted with diseased scion; B, one of the three trees grafted with healthy scion.

https://doi.org/10.1371/journal.pone.0133652.g001

DNA extraction and conventional PCR detection

About 0.1 g midribs of leaves were cut into pieces and mixed with CXD buffer provided by HP plant DNA kit (D2485-02, OMEGA bio-tek). Samples were pulverized using MP FastPrep-24 machine, and DNA was extracted according to manufacturer’s instruction. The final DNA pellets were dissolved in 100 μL elution buffer and checked by electrophoresis in 1% (w/v) agarose gel. The concentration and quality of DNA were determined with an Agilent Model 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA), and all samples were adjusted to a final concentration of 50 ng/μL. Detection and identification of the ‘Ca. L. asiaticus’ was based on primer pairs OI1/OI2c [14].

RNA extraction, array analysis, and digital gene expression

Around 0.5 g of leaves grounded by mortar and pestle in liquid nitrogen was used for total RNA extraction. RNA manipulation was done according to SDS method with few modifications. Leaf samples were ground in 0.7 mL of extraction buffer mixed with 2% SDS, 0.02 mol/L Borax, 4% β- Mercaptoethanol, 3% Tris-phenol, and 2.2% Chloroform. Fine tissue was vortex mixed for 3 min at 6,000 rpm, centrifuged at 12,000 rpm for 5 min (4°C). The aqueous phase was transferred to a new tube, and then 0.6 mL of Tris-phenol and Chloroform were added separately, followed by 2 min of vortex and 5 min of spin respectively. The RNA was precipitated from the aqueous phase by mixing with 4M LiCl and ethanol. Diluted nucleic acid was treated with 3M NaAc and Ethanol on ice for 10 min. The RNA pellet was washed once with 75% Ethanol (in DEPC-treated water), inverted, mixed and centrifuged at 12,000 rpm for 30 minutes at 4°C. The RNA pellet was dried at room temperature and dissolved in 60 μL RNase-free water by passing the solution through the pipette tip for a few times. All RNA samples was qualified by Agilent Model 2100 Bioanalyzer and kept in -70°C for further use.

Real time PCR of 25 citrus genes (mentioned in [15]) reported to be differentially expressed in the previous microarray experiment [3–6] were carefully carried out using the 6 samples collected here to test the gene expression consistency within biological repeats. No significant difference was observed among the ‘Ca. L.asiaticus’ treated or non-treated samples for all 25 genes including 2 housekeeping candidates. Consequently, the same amounts (ng) of three infected RNA samples or healthy samples were mixed respectively. Mixed RNA samples were sent to a biological company (CapitalBio corperation, Beijing) for cDNA preparation, hybridization, scanning, and image analysis of the arrays according to the manufacturer’s recommendations (Affymetrix, Santa Clara, CA, USA). The microarray hybridization was based on 33,879 expressed sequence tag sequences from different citrus (http://www.affymetrix.com/products/arrays/specific/citrus.affx). Differentially expressed genes were screened with a cutoff threshold P value of 0.05 and log2 fold change (LFC) of≧2.00 or ≦-2.00.

Reverse transcription and quantitative RT-PCR

Three μg RNA of each sample was reverse transcribed using ReverTra Ace qPCR RT Kit (TOYOBO, Code No.: FSQ-101). The cDNA samples were quality-checked and quantified using Agilent Model 2100 Bioanalyzer for a second time to adjust the concentration of each sample to 1000 ng/μL. Recombinant DNase I(Takara, D2270A)was served for removal of residual genomic DNA fragments that would interrupt definition of semi-quantitative PCR. To screen the primers, cDNA samples of 2-year after inoculation, the 3-year after inoculation, and healthy controls were mixed separately as templates. Primers with good dissociation curve, amplification curve and without nonspecific amplification were selected for all samples. PCR products were recovered (BioTeke, DP1602) and qualified (Agilent Model 2100 Bioanalyzer) for the 1:10 dilution series to construct the standard curve of each gene.

The 18S rRNA gene was used as an internal control. Primers for 15 selected differentially genes designed using Primer Premier 5 were showed in S1 Table. Programs set for conversional PCR was an initial 5 minute denaturation at 95°C followed by 30 cycles of thermocycling (95°C, 30s; 55°C, 30s; 72°C, 30s), then extended at 72°C for 7 min. Amplifications were carried out in a BIO-RAD C1000 Thermal Cycler. Quantitative PCR (qPCR) was carried out in the Bio-Rad CFX Connect Real-Time System with SYBR Premix Ex Taq (Ti RNaseH, Plus) (DRR820B). Standard curve, melting curve, amplification curve, and Ct(CP), Ct(SDM), Qty(CP), and Qty(SDM) were collected for result analyses. Six replicates were used for each sample (three biological replicates × two technical replicates). The relative quantification of target gene transcripts was determined by the comparative CT method 2−ΔΔCT [16], where ΔΔCT = (CT, Target–CT, 18S rRNA)HLB−(CT, Target–CT, 18S rRNA)Control.

Data analysis

Expression signals were normalized by Robust Multichip Analysis (RMA) approach[17]. Comparison analysis was performed by the signal ration of treated sample and control sample to identify the differentially expressed genes (DEGs). Genes with ratio of ≤0.5 were considered down-regulated DEGs, while genes with ratio of ≥2 were up-regulated DEGs. The FDR-adjusted significance level cutoff was set at 0.05 for overrepresentation determination. The manual functional analysis was referred to the system of Adams, Kerlavage [18]. Functional enrichment/overrepresentation analysis was carried out using the agriGO database (http://bioinfo.cau.edu.cn/agriGO/) [19].

Lipoxygenase activity determination

‘Ca. L. asiaticus’ infected citrus leaves and healthy leaves collected every two months during one year (From May, 2013) were used for lipoxygenase (LOX) determination. The activity of LOX was detected according to Chen, Xu [20]. Briefly, every two gram of leaf tissue was grinded in 10 mL cold phosphate buffer (50 mmol/L, pH 7.0) with individual mortar and pestle on ice. Liquid supernatant was collected after centrifugation at 15000 g for 15 min for LOX activity detection. The hydroperoxy lipid product of the reaction contains a conjugated diene which absorbs strongly at 234 nm. The reaction was carried out in a solution containing 25μL Linoleic acid sodium salt (25 mmol/L), 2.775 mL acetic buffer (100 mmol/L, pH 5.5) in a final volume of 3 mL. By measuring product formation through the change in absorbance at this wavelength, the quantitative assay of LOX can thus be determined.

Light microscopy

Mature leaves, all in similar size, with typical mottling and Znic-deficiency-like symptoms were collected from mandarin (Citrus reticulate Blanco) trees. Leaves from healthy trees were collected as control. Quantitive real time PCR was conducted to confirm the present of ‘Ca. L. asiaticus’. Around 5 mm midribs of each leaf sample were transferred into FAA solution (Volume ratio of 70% ethanol, acetic acid, and 37% formaldehyde was 90: 5: 5) for 48 h. After fixative treatment, samples were gone through paraffin in sections workflow: dehydration, clearing, waxfilling, embedding, and sectioning. Sections were strained with 0.5% dye Fast Green and 1% dye Safranine. Light micrographs were scanned using a light micrograph (Carl Zeiss, West Germany) with an attached camera.

Results and Discussion

Host metabolisms were affected by ‘Ca. L. asiaticus’ infection

Of the 2,017 significantly regulated DEGs, 938 (46.5%) were up-regulated, and 1,080 (53.5%) were down-regulated in infected Jiaogan (All Affymatrix data was submitted to GenBank with GEO accession numbers of GSE67376, http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE67376). The 1,364 genes affected by ‘Ca. L. asiaticus’ in the late stage of infection having known functions were related to plant defense/stress response (23.6%), carbohydrate (sugar/starch/xylan) metabolism (7.9%), photosynthesis (4.4%), plant growth and root/steam/leaf/flower/seed development (7.9%), structure (cell wall metabolism/nucleus organization/ribosomal protein/fiber protein/membrane protein/vacuolar protein) (6.4%), signaling transduction/translation factors (14.5%), DNA/chromatin/nucleotide/ADP/ion/protein/carbohydrate binding (3.9%), sugar/protein/lipid/ion/oligopeptide/dipeptide/tripeptide/amino acid/ATP/ADP/bilirubin/carbohydrate/drug/electron transport (11.9%), oxidation-reduction process (3.2%) with oxidoreductase activity, other metabolic processes (12.8%), and others (10.6%). The expression pattern (up- or down-regulation) of these classified DEGs in the HLB infected trees are shown in Fig 2A. Genes related to photosynthesis and structure were mostly down-regulated by the ‘Ca. L. asiaticus’ infection, whereas genes associated with the oxidation-reduction process and transport were mostly up-regulated.

thumbnail
Download:
Fig 2. Functional categories and Gene ontology classification of all differentially expressed genes in response to infection with‘Candidatus Liberibacter asiaticus’.

A, Classification of specifically expressed genes into functional categories in response to ‘Ca. L. asiaticus’. B, Gene ontology classification of all differentially expressed genes. DEGs, differentially expressed genes; BP-1, cellular nitrogen compound metabolic process; BP-2, carbohydrate metabolic process; BP-3, polysaccharide metabolic process; BP-4, cellular carbohydrate metabolic process; BP-5, cellular polysaccharide metabolic process; BP-6, cellular glucan metabolic process; BP-7, glucan metabolic process; BP-8, starch metabolic process; BP-9, sucrose metabolic process; BP-10, oligosaccharide metabolic process; BP-11, disaccharide metabolic process; BP-12, electron transport; BP-13, glycoside metabolic process; BP-14, photosynthesis; BP-15, pigment biosynthetic process; BP-16, pigment metabolic process; BP-17, secondary metabolic process; BP-18, phenylpropanoid metabolic process; BP-19, heterocycle biosynthetic process; MF-1, hydrolase activity, hydrolyzing O-glycosyl compounds; MF-2, hydrolase activity, acting on glycosyl bonds; MF-3, transferase activity, transferring glycosyl groups; MF-4, oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen; MF-5, catalytic activity; MF-6, oxidoreductase activity; MF-7, transferase activity, transferring hexosyl groups; CC-1, membrane-bounded vesicle; CC-2, thylakoid; CC-3, cytoplasmic membrane-bounded vesicle; CC-4, vesicle; CC-5, cytoplasmic vesicle; CC-6, thylakoid part; CC-7, chloroplast thylakoid; CC-8, organelle subcompartment; CC-9, plastid thylakoid; CC-10, photosynthetic membrane; CC-11, chloroplast; CC-12, cell wall; CC-13, external encapsulating structure; CC-14, thylakoid membrane; CC-15, plastid part; CC-16, chloroplast part; CC-17, chloroplast thylakoid membrane; CC-18, plastid thylakoid membrane; CC-19, plastid; CC-20, membrane; CC-21, photosystem.

https://doi.org/10.1371/journal.pone.0133652.g002

Gene ontology (GO) assignments were used to classify the functions of the DEGs (S2, S3, and S4 Tables), which were assigned to 16 functional groups in the ‘biological processes’ main category, 21 groups in ‘cellular components’ and 7 groups in molecular function’ (Fig 2B). The majority of the GO annotations were positioned in the following categories: photosynthesis, electron transport and the metabolic processes of nitrogen compounds, carbohydrates, polysaccharides, glucans, starch, sucrose, oligosaccharides, disaccharides, glycosides, pigments, and phenylpropanoids and heterocycles.

The significant up-regulation (21.25-fold) of a gene coding the glutamic acid-rich protein cNBL1700, which may be related to pathogenesis [21], was observed. In the interaction with ‘Ca. L. asiaticus’, citrus expressed some defense-related genes in response to the stress. Two thi1 genes (Cit.465.1 and Cit.377.1), which may be involved in bacterial-induced DNA damage tolerance in plant cells [22], were up-regulated more than 5 times. Some heat shock protein genes, involved in the regulation of transcription were highly up-regulated, more than 10-fold. The stress related HSP genes were also mentioned in previous studies (induced in [3], but down-regulated in [9]), which functions in plant growth and development [23–27]. However, the expression levels of some DnaJ/Hsp genes were highly depressed.

‘Ca. L. asiaticus’ depressed the photosynthesis process of mandarin

The photosynthesis process was depressed as a result of HLB-induced incomplete coloring (S5 Table). Photosynthesis-related genes showed opposite expression patterns in leaves (mainly down-regulated) and fruit (mainly up-regulated) in the previous studies. Among the 55 genes related to the photosynthetic process, only 4 (Cit.3817.1, Cit.3817.1, Cit.7266.1 and Cit.17826.1) had increased expression levels after two years of ‘Ca. L. asiaticus’ infection. HLB interrupts the citrus metabolic system and causes complete (pale yellowing) or incomplete (blotchy mottle) coloring symptoms on leaves. HLB-infected samples collected in this study displayed typical blotching symptoms. Thus, it is not surprising that photosynthesis-related proteins showed the greatest decrease in response to ‘Ca. L. asiaticus’ infection. This is consistent with results from most previous studies that focused on HLB-affected leaves. Most DEGs in this category were photosystem II factors. Seven genes coding light-harvesting chlorophyll a/b binding proteins, including Lil3, cab-11, CP24 and CP29, were repressed in HLB samples compared to healthy samples. The chl a/b-binding protein plays an essential role in chlorophyll and tocopherol biosynthesis reactions in photosystem II [28]. Genes coding chlorophyll a oxygenase, lipoic acid synthase, oxygen-evolving enhancer protein, or Ribulose bisphosphate carboxylase/oxygenase (Rubisco) activase, were down-regulated. Chlorophyll a oxygenase is involved in chlorophyll b formation from chlorophyll a [29]. Additionally, we observed a significant decrease in the accumulation of some subunits of the photosynthetic system in the HLB leaves, including the cytochrome b6f complex subunit (Cit.9493.1 and Cit.40088.1), the photosystem I-N subunit (Cit.8714.1 and Cit.30551.1), oxygen-evolving enhancer protein 2 (Cit.24174.1) and magnesium-chelatase subunit chlI (Cit.21132.1). Violaxanthin de-epoxidase was reported to be involved in the xanthophyll cycle under stress [30]. It is consistent with results from a previous study by Nwugo [9] that Rubisco activase coding genes were repressed. Rubisco activase is a catalytic chaperone that affects photosynthetic induction, nucleotides and sugar phosphates [31, 32]. Photosynthesis is central to all aspects of plant biology. The strong reduction in photosynthesis may relate to the symptoms of the diseased leaves.

HLB induced the synthesis of starch and sucrose in leaves

Photosynthesis in plants directly influences the carbohydrate metabolic product. Numerous genes with GO terms from the carbohydrate metabolic process (starch anabolism, sucrose, cellulose catabolic process, allantoin catabolic process and sugar metabolism) were differentially expressed in response to HLB (S6 Table). The β-amylase and β-mannosidase enzymes are important for starch degradation [33], while starch branching enzyme, starch phosphorylase type H, granule-bound starch synthase, UDP-glucosyl transferase, UDP-glucose dehydrogenase and isoamylase [34]were related to starch synthesis. Most down-regulated DEGs in this study that were categorized in carbohydrate metabolic process functioned as degrading enzymes. The down-regulated DEGs that had been identified in previous studies were four β-amylase genes (coding an important enzyme in starch degradation [35]), alcohol acyl transferase genes and glycosyl transferase family 8 genes involved in sugar metabolism. The up-regulated genes, coding ADP-glucose pyrophosphorylase, α-amylase, β-glucosidase, glycosyl hydrolase family 38, granule-bound starch synthase, starch branching enzyme, starch phosphorylase type H, UDP-glucuronosyltransferase and starch synthase, had been identified in previous studies on HLB-infected leaves [3, 5, 6] but not on HLB-infected fruits [7]. Five of the differentially expressed β-glucosidase genes were over-expressed (3 by more than 18-fold) in HLB-affected leaf tissues. The genes of alcohol acyl transferase, alcohol dehydrogenase 2, glycerol-3-phosphate dehydrogenase, glycosyl hydrolase family protein, raffinose synthase, fructose-bisphosphate aldolase and trehalose-6-phosphate phosphatase were repressed, as previously reported [5, 9]. However, fructose-bisphosphate aldolase, β-xylosidase, glycosyl hydrolase family protein, β-xylosidase, and β-1,3-glucanase genes were up-regulated by HLB in the study of Albrecht and Bowman [3, 6] but showed no regular pattern at this extremely late infection stage.

Several DEGs that were classified in the carbohydrate metabolic process had not been previously reported in HLB infections. The depressed DEGs were mannosidase genes, xyloglucan-specific fungal endoglucanase inhibitor genes, ribose 5-phosphate isomerase, mannose-6-phosphate isomerase genes and glucan 1,3-β-glucosidase genes. The DEGs specially induced at this extremely late stage included PG1, β-fructofuranosidase, bifunctional lysine-ketoglutarate reductase, glucosyltransferase, exocellular acid invertase 1, limonoid UDP-glucosyltransferase, β-galactosidase and hexose transport protein HEX6. As a cellulose catabolic process gene, β-galactosidase may also play a role during abscission, and growth and development processes in flowers and fruitlets [36]. An exocellular acid invertase 1 with fructan β-fructosidase activity [33] and a limonoid UDP-glucosyltransferase with sinapate 1-glucosyltransferase activity had expression levels that were increased 6.5 and 5.4 times after HLB infection. Limonoid UDP-glucosyltransferase is supposed to be highly associated with the sweet or sour taste of citrus fruits [37]. The up-regulation of key starch synthases and sugar biosynthetic genes, together with the down-regulation of related degrading enzymes, likely leads to the accumulation of starch in HLB-affected leaves and the sour taste in the affected fruits.

Plant growth and development was influenced by ‘Ca. L. asiaticus’ infection

Photosynthesis is an important biological process in plants because it provides energy for plant growth. About 7.9% of the DEGs induced by ‘Ca. L. asiaticus’ at the very late stage were classified in the ‘plant growth and development’ category. All 4 differentially expressed flavonol synthase genes were induced. The flavonol synthase genes from the ‘Satsuma’ mandarin were differentially regulated in the fruit developmental stage in a tissue-specific manner [38]. At 25 months after infection, genes encoding bifunctional nuclease, CCAAT-box-binding transcription factor-related protein, GGPP synthase, GPI-anchored protein, NAC domain protein NAC2 and SRG1 protein were highly induced, while genes of a TCP family transcription factor, two rapid alkalinization factor 2s, one growth-regulating factor, two rapid specific-tissue protein 2s, four arabinogalactan proteins and two ent-kaurene oxidases were significantly suppressed. The CCT (for CONSTANS, CONSTANS-LIKE and TOC1) domain is involved in photoperiodic flowering, light signaling and circadian rhythm regulatory processes in Arabidopsis thaliana. It is similar to the yeast HEME ACTIVATOR PROTEIN2 (HAP2) that binds to CCAAT boxes in eukaryotic promoters [39]. Levels of bifunctional nuclease 1 mRNA were significantly different in different plant tissues and were induced during leaf and stem senescence [40]. Similarly, the up-regulated GGPP synthase functioned in the leaf senescence in Arabidopsis [41]. The GPI-anchored protein was involved in the pollen tube growth process [42]. SRG1 proteins involved in embryo development ending in seed dormancy or organ senescence were also induced in sweet orange infected by ‘‘Ca. L. asiaticus” [4]. Two GO results for the TCP family transcription factor are the negative regulation of leaf senescence (GO:1900056) and root development (GO:0048364). Additionally, the plant growth and development-related arabinogalactan proteins were reported to be repressed in ‘Valencia’ orange leaves responding to the HLB pathogen [3]. The nt-kaurene oxidase is a gibberellin biosynthesis enzyme [43] that functions in the gibberellin biosynthesis pathway, and subsequently influences plant growth. Rapid alkalinization factors were reported in the brassinosteroid-mediated signaling pathway that negatively regulates plant growth. No highly induced or suppressed genes were found, which may due to the long time interaction of citrus plants with the pathogen that may have allowed the host to develop a maintenance level towards this stress.

Plant cell organelles were damaged by the 2 year long HLB infection

Also noticeable among the highly down-regulated genes were those for ribosomal proteins, membrane proteins, thioredoxins and vacuolar proteins. In total, 81 DEGs were classified as structural proteins. Most (80.2%) of the structures were depressed in this very late stage of infection (S7 Table). Among the DEGs, the ribosomal proteins were especially prominent (24 of the 28 were down-regulated). However, the ribosomal proteins involved in protein metabolism were up-regulated in the study of Albrecht and Bowman [6]. Transcripts coding for hydroxyproline-rich glycoproteins were also noticeably down-regulated. Thioredoxins, components of the plant cell redox regulatory system that function as antioxidants in ROS scavenging, were up-regulated in HLB-affected sweet orange fruit [8] but down-regulated in lemon leaves [10].

The midribs of HLB-affected and uninfected mandarins were stained and observed with LM (Fig 3). Significant disorder of the cells and collapse of sieve tubes and companion cells were observed. The phloem layer and cambium layer of infected samples were thicker than uninfected midribs. Starch accumulation was more serious in the diseased midribs. These anatomical changes in ‘Ca. L. asiaticus’ affected samples were similar but stronger than previous report [5].

thumbnail
Download:
Fig 3. Anatomical analysis (Fast Green and Safranine O) of midrib phloem tissues of the huanglongbing (HLB) affected and healthy Citrus reticulata Blanco using light microscopy.

A, midrib of healthy leaf. Vascular bundle is in oval-shaped, all tissues are well developed, the interior of parenchyma cells are relatively clear. B and C, Huanglongbing infected midribs from leaves with blotchy mottle or Zinc-deficiency-like symptoms. Parenchyma contains more internal substances like starchs, half round of vascular bundle is poor developed and phloem cambium is thickened.

https://doi.org/10.1371/journal.pone.0133652.g003

Transcription/signal transduction-related genes were alerted

A total of 82 transcriptional genes or genes with transcription factor activity were differentially regulated after two years of ‘Ca. L. asiaticus’ infection (S8 Table). Of these, 46 were down-regulated. The MADS-box protein 4 coding gene (Cit.144.1) was induced by 13-fold. A TAZ zinc finger family protein (Cit.6280.1.S1), with a GO classification as a negative regulator of transcription from the RNA polymerase II promoter, was increased by11.8-fold. However, for the DEGs encoding certain transcription factors, 3 WRKY family, 5 myb family, 10 AP2 domain, 10 bHLH and 3 heat shock, as well as 6 genes annotated as unknown proteins with transcription factor activity and 3 zinc finger (B-box type) family proteins, the expression profiles were not consistent within the groups. Consequently, no conclusions about transcription factor-regulated signaling pathways were drawn, unlike in some studies focusing on earlier HLB infection stages. However, the transcripts of differentially regulated C2H2-type and GATA-type zinc finger proteins, auxin/indole-3-acetic acid proteins and squamosa promoter-binding proteins were depressed. Whereas, three genes encoding scarecrow transcription factor family proteins, three genes of the bZIP transcription factor family, and two genes encoding dehydration-responsive element binding proteins [44, 45], were all up-regulated.

In addition to the transcription factor genes, another 107 genes encoding proteins involved in the regulation of transcription processes or signal transduction were classified in the signaling category. For these genes, most (69 of 107) were down-regulated. The most prominent groups included 14 of 15 genes encoding unknown proteins, 6 genes for ethylene responsive element binding proteins or ethylene-induced esterases, and 4 dormancy-associated protein genes. The expression of other genes was also suppressed, including two copies of the genes encoding transducin-like, GTP cyclohydrolase, allergenic isoflavone reductase-like protein Bet v 6.0102 and gibberellin-regulated proteins. All three Avr9/Cf-9 rapidly elicited proteins (Cit.6781.1, Cit.15649.1 and Cit.16734.1) were decreased by more than 2.8-fold. Avr9/Cf-9 rapidly elicited proteins are thought to be involved in the ionotropic glutamate receptor signaling pathway and hypersensitive responses in plant cell death, a disease-resistance pattern [46].

The ‘Ca. L. asiaticus’ infection induced the transcription of 38 genes in signaling pathways, such as genes encoding 33 kDa secretory protein-related, α-expansin 3, TMKL1 precursor, ionotropic glutamate receptor homolog GLR4, response regulator 6, RESPONSE REGULATOR 9, cyclic nucleotides, calmodulin-regulated ion channel proteins and homogentisate geranylgeranyl transferase. Four no apical meristem (NAM)-like protein genes and 3 NAM family protein genes had increased transcript levels. No regular patterns for the responses of the 8 serine/threonine kinases, 4 isoflavone reductases or 11 protein kinases to ‘Ca. L. asiaticus’ were observed.

Other metabolic processes influenced by HLB

DEGs in this category were related to the biological processes of terpene synthase, lipid metabolism, cell/tissue development and cell wall metabolism, acyl-carrier-protein biosynthesis, asparagine biosynthesis, sterol biosynthesis, fatty acid biosynthesis, isoprenoid biosynthesis, chitin metabolism, Mo-molybdopterin cofactor biosynthesis, flavonol biosynthesis, UMP biosynthesis and other secondary metabolisms. Of these, 85 (48.9%) were depressed during the late stage of ‘Ca. L. asiaticus’ infection. Among the down-regulated genes, those encoding 3 acid phosphatases, 6 α- or β-tubulins, 3 early nodulinses, 4 pectate lyases, 7 pectin esterases, and 2 cinnamoyl-CoA reductases were apparent. The expression pattern of the pectin esterase-coding genes was different in fruit. Consistent with the study of Nwugo [9], cinnamoyl-CoA reductase genes and tubulin genes were significantly depressed in response to ‘Ca. L. asiaticus’. However, the energy/metabolisms-related genes of the cytosolic aldehyde dehydrogenase and dienelactone hydrolase family proteins/carboxymethylenebutenolidase were not increased in the study by Nwugo [9].

Significantly, levels of transcripts for the two CCR proteins [47, 48], two catechol O-methyltransferases (involved in synthesis of the flavor molecule guaiacol [49, 50]), and two caffeic acid O-methyltransferase IIs, which function in lignin, flavonoids and sinapoyl malate biosynthesis [51, 52], were increased more than 5-fold. Cit.21988.1.S1, which is involved in fatty acid/lipid and amino acid metabolism, was induced 3.4-fold in this study, but was decreased in fruit [8]. The four genes encoding cellulose synthases, related to cell organization and biogenesis, and the UDP-glucose metabolic process, were reported to be depressed in Albrecht and Bowman (3], but were induced here. A fatty acid/lipid and amino acid metabolism gene encoding family II extracellular lipase 3 was highly (7.13-fold) induced here, as well as in fruit [8]. The other 4.5-fold and 7.1-fold up-regulated genes FPh1, related to the phospholipase/carboxylesterase family leucoanthocyanidin dioxygenase-like protein [53] and O-methyltransferase family 2 protein, respectively, were also induced by ‘Ca. L. asiaticus’ in other studies. The up-regulation of oxidoreductase, 2-oxoglutarate and Fe(II) oxygenase family proteins and lipase class 3 family protein genes, as well as the down-regulation of monogalactosyldiacylglycerol synthase and H+-transporting ATP synthase genes were consistent with previous research on ‘Ca. L. asiaticus’[6, 8, 9].

HLB alteration of the hormone-mediated immune response

Hormone signaling pathways are reported to strongly relate to plant disease progress. A lot of phytohormone related genes were differentially expressed in mandarins plants infected with ‘Ca. L. asiaticus’ (S9 Table). The expression of genes encoding three ethylene-responsive element binding proteins, 1-aminocyclopropane-1-carboxylate oxidase, ethylene response factor ERF3b and ethylene-induced esterase, were all down-regulated by the ‘Ca. L. asiaticus’ infection. Additionally, all five genes encoding aquaporin proteins had decreased expression levels. Aquaporin PIP1 interacts with PIP2 in ethylene-regulated plant growth [54]. Only a gene encoding an ethylene-forming-enzyme-like dioxygenase-like protein was induced under the long term infection. Consequently, this study suggests that the ethylene pathway was depressed in this study, which is different from previous studies. Reduced transcription was observed for genes encoding a NPR1/NIM1-interacting protein 3, a BTB/POZ domain-containing protein, a S-adenosyl-L-methionine:salicylic acid carboxyl methyltransferase, a UDP-glucose:salicylic acid glucosyltransferase, a catalase, Nine glutathione peroxidase and a Cu/Zn superoxide dismutase that are involved in salicylic acid signaling. Moreover, the transcriptional patterns of salicylic acid-mediated defense response-related proteins, PR5-1, PR1 and PR10A, were different [55]. To conclude, a two year ‘Ca. L. asiaticus’ infection triggered a response via both salicylic acid and ethylene pathways.

Seventeen Of the 20 DEGs related to the jasmine acid signaling pathway (S9 Table)., including four LOX coding genes, three NAM family protein genes and seven glutathione S-transferase genes were predominately up-regulated compared with in non-infected trees. Three replicates of leaves collected every two months from HLB infected trees and healthy trees were used to determine the LOX activity. The one-year data showed that the concentration of LOX in the diseased trees was relatively higher than that in healthy trees among the year (Fig 4). This is consistent with the array data.

thumbnail
Download:
Fig 4. Evaluation of the time course lipoxygenase (LOX) activity in huanglongbing (HLB) affected and healthy Citrus reticulata Blanco leaves.

https://doi.org/10.1371/journal.pone.0133652.g004

Strikingly, eight genes corresponding to auxin-regulated proteins were down-regulated, while genes encoding auxin-repressing proteins were up-regulated (S9 Table). This indicates that the auxin-regulated pathways were suppressed after 2 years of a ‘Ca. L. asiaticus’ infection. Similarly, several genes involved in the signal transduction of gibberellins and abscisic acid were induced while some of them were depressed. A gibberellin-responsive GAST1 homolog was 11.1-fold up-regulated. This gene was predicted to be involved in brassinosteroid, gibberellin and abscisic acid responsive expression by GO analysis.

Other stress-related genes are differentially regulated by ‘Ca. L. asiaticus’

The two year ‘Ca. L. asiaticus’ infection up-regulated the expression of 149 plant defense or stress-response genes, including homologs of protein kinase Xa21, SAR related N-hydroxycinnamoyl/benzoyltransferase (4.35-fold), disease resistance protein (also up-regulated in previous studies of HLB-citrus interactions), glutathione S-transferase, EIX receptor 2 [56, 57], hypersensitive-related kinases, like HSR203J [58], abiotic stress-related calmodulin [59] and resistance-related RPP8-like protein. Consistent with some of the previous studies [4, 10], the following genes were up-regulated: a gene coding receptor-like kinase CHRK1 (6.3-fold), five ATPase genes (7.5-, 4.6-, 3-, 2.6- and 2.4-fold), and eight stress/defense response chitinase genes (10.4-, 6.8-, 5.2-, 5.1-, 4.5-, 4.1-, 2.6- and 2.5-fold), involved in cell organization and biogenesis.

Among the genes down-regulated more than 2-fold were enzymes associated with aquaporin (5), aspartyl protease family protein (2), subtilisin-like protease (3), invertase/pectin methylesterase inhibitor (4) [60, 61], leucine rich repeat protein (9 of 11), which was different from in Fan et al. [55, 62], and miraculin-like protein 2 (9 of 10) (induced in[10] and [4]). A gene putatively coding an alanine-glyoxylate aminotransferase, which functions in the rice-blast interaction was 7.3-fold induced [63].

qRT-PCR validation

To corroborate the microarray data, 15 genes classified in different functional groups were amplified using the quantitative real time PCR method. Samples used were from Jiaogan leaves infected by ‘Ca. L. asiaticus’ for 2 or 3 years, and their respective non-infected controls. The gene coding glutamic acid-rich protein cNBL1700, a pathogenesis related gene, which was induced 21.25 times when screened using a microarray, was 15.8±3.47 and 5.13±0.79 times up-regulated after 2 and 3 years, respectively. A chalcone synthase gene, related to plant growth and development, was detected to be highly depressed in the microarray result, and the trends were the same, 5.22±0.02 and 20.01±7.00 fold down-regulated after 2 or 3 years, respectively, for the real time PCR results. The expression profiles of all 15 genes were consistent with the microarray data when using the same samples tae at 2 years after infection (Fig 5). The expression patterns of most genes in the samples taken 3 years after infection were similar to, but stronger than, those taken after 2 years (Fig 5).

thumbnail
Download:
Fig 5. Real time PCR and microarray profiles of 15 selected differentially expressed genes (DEGs).

A, ethylene pathway related genes; B, jasmine acid pathway related genes; C, photosynthetic related genes; D, transport related genes; E, starch degradation or synthetic genes; F, plant defense and stress related genes. A-1, ethylene responsive element binding protein; A-2, ethylene response factor ERF3b; A-3, ethylene responsive element binding protein; B-1, no apical meristem (NAM) family protein; B-2, glutathione S-transferase GST 14; C-1, NADPH-protochlorophyllideoxidoreductase; C-2, light-harvesting chlorophyll a/b binding protein; C-3, chlorophyll synthase, putative; D-1, glucose-6-phosphate/phosphate translocator; D-2, glucose-6-phosphate/phosphate-translocator precursor; E-1, beta-amylase; E-2, starch branching enzyme; F-1, leucine rich repeat protein.

https://doi.org/10.1371/journal.pone.0133652.g005

Conclusions

The broad range of DEGs suggests that citrus trees were profoundly disturbed after 2 years of infection with the bacterial disease HLB. Firstly, the photosynthesis process was most dominantly depressed. Consequently, the carbohydrate metabolic process, and plant growth and development were influenced. The ‘Ca. L. asiaticus’ infection also down-regulated the transcription of genes encoding structural protein and disrupt the cell of midribs. Thirdly, the ethylene and salicylic acid signaling pathways were depressed while the jasmine acid signaling pathways were induced. Finally, an analysis of differentially-regulated signaling genes and enriched stress-related genes demonstrated that multiple regulatory pathways, including NAM, nectarin 5, tubulin, chitinase CHI1, and miraculin-like protein 2 regulons, are involved in the very late stage of ‘Ca. L. asiaticus’ infection. Compared to earlier infections, the metabolic processes of nitrogen compounds, polysaccharides, oligosaccharides, disaccharides, glycosides, pigments, phenylpropanoids and heterocycles were especially activated by ‘Ca. L. asiaticus.

Supporting Information

S1 Table. Quantitative RT-PCR primer sets used to determine differential expression between mandarins infected by ‘Candidatus Liberibacter asiaticus’ for 2 years and healthy cirrus.

https://doi.org/10.1371/journal.pone.0133652.s001

(XLS)

S2 Table. Gene ontology enrichment result of the differentially expressed genes.

P, Biological Process; F, Mollecular Function; C, Cellular Componen.

https://doi.org/10.1371/journal.pone.0133652.s002

(XLS)

S3 Table. Gene ontology enrichment result of the down regulated genes.

F, Mollecular Function.

https://doi.org/10.1371/journal.pone.0133652.s003

(XLS)

S4 Table. Gene ontology enrichment result of the up-regulated genes.

P, Biological Process; F, Mollecular Function; C, Cellular Componen.

https://doi.org/10.1371/journal.pone.0133652.s004

(XLS)

S5 Table. Functional classification of photosynthetic process related genes with differential expression (2-fold) in mandarins plants infected with ‘Ca. L. asiaticus’ compared with non-infected plants 25 months after inoculation.

https://doi.org/10.1371/journal.pone.0133652.s005

(XLS)

S6 Table. Huanglongbing regulation of genes related to carbohydrate metabolic.

https://doi.org/10.1371/journal.pone.0133652.s006

(XLS)

S7 Table. Huanglongbing regulation structure protein genes.

https://doi.org/10.1371/journal.pone.0133652.s007

(XLS)

S8 Table. Huanglongbing regulation of transcriptional factors and signaling genes.

https://doi.org/10.1371/journal.pone.0133652.s008

(XLS)

S9 Table. The plant hormone related genes regulated by 25 months’ Huanglongbing infection.

SA, Salicylic acid; JA, Jasmine acid; GI, Giberellin; ET, Ethylene; CT, Cytokinin; BR, Brassinosteroid; AU, Auxin; IAA, indole-3-acetic acid; ABA, Abscisic acid.

https://doi.org/10.1371/journal.pone.0133652.s009

(XLS)

Acknowledgments

The authors acknowledge Dr. Jianchi Chen from United Department of Agriculture, Agricultural Research Service, for editorial assistance.

Author Contributions

Conceived and designed the experiments: MX YL XD. Performed the experiments: MX YL ZD YT. Analyzed the data: MX ZZ. Contributed reagents/materials/analysis tools: MX YL XD. Wrote the paper: MX.

References

  1. 1. Talon M, Gmitter FG Jr. Citrus genomics. International journal of plant genomics. 2008;2008:528361. pmid:18509486; PubMed Central PMCID: PMC2396216.
  2. 2. Xu Q, Chen LL, Ruan X, Chen D, Zhu A, Chen C, et al. The draft genome of sweet orange (Citrus sinensis). Nature genetics. 2013;45(1):59–66. pmid:23179022.
  3. 3. Albrecht U, Bowman KD. Gene expression in Citrus sinensis (L.) Osbeck following infection with the bacterial pathogen Candidatus Liberibacter asiaticus causing Huanglongbing in Florida. Plant Science. 2008;175(3):291–306.
  4. 4. Fan J, Chen C, Yu Q, Brlansky RH, Li ZG, Gmitter FG Jr. Comparative iTRAQ proteome and transcriptome analyses of sweet orange infected by "Candidatus Liberibacter asiaticus". Physiologia plantarum. 2011;143(3):235–45. pmid:21838733.
  5. 5. Kim JS, Sagaram US, Burns JK, Li JL, Wang N. Response of sweet orange (Citrus sinensis) to 'Candidatus Liberibacter asiaticus' infection: microscopy and microarray analyses. Phytopathology. 2009;99(1):50–7. pmid:19055434.
  6. 6. Albrecht U, Bowman KD. Transcriptional response of susceptible and tolerant citrus to infection with Candidatus Liberibacter asiaticus. Plant science: an international journal of experimental plant biology. 2012;185–186:118–30. pmid:22325873.
  7. 7. Martinelli F, Uratsu SL, Albrecht U, Reagan RL, Phu ML, Britton M, et al. Transcriptome profiling of citrus fruit response to huanglongbing disease. PloS one. 2012;7(5):e38039. pmid:22675433; PubMed Central PMCID: PMC3364978.
  8. 8. Liao HL, Burns JK. Gene expression in Citrus sinensis fruit tissues harvested from huanglongbing-infected trees: comparison with girdled fruit. Journal of experimental botany. 2012;63(8):3307–19. pmid:22407645; PubMed Central PMCID: PMC3350938.
  9. 9. Nwugo CC, Lin H, Duan Y, Civerolo EL. The effect of 'Candidatus Liberibacter asiaticus' infection on the proteomic profiles and nutritional status of pre-symptomatic and symptomatic grapefruit (Citrus paradisi) plants. BMC plant biology. 2013;13:59. pmid:23578104; PubMed Central PMCID: PMC3668195.
  10. 10. Nwugo CC, Duan Y, Lin H. Study on citrus response to huanglongbing highlights a down-regulation of defense-related proteins in lemon plants upon 'Ca. Liberibacter asiaticus' infection. PloS one. 2013;8(6):e67442. pmid:23922636; PubMed Central PMCID: PMC3686881.
  11. 11. Mafra V, Martins PK, Francisco CS, Ribeiro-Alves M, Freitas-Astua J, Machado MA. Candidatus Liberibacter americanus induces significant reprogramming of the transcriptome of the susceptible citrus genotype. BMC genomics. 2013;14:247. pmid:23586643; PubMed Central PMCID: PMC3635983.
  12. 12. Aritua V, Achor D, Gmitter FG, Albrigo G, Wang N. Transcriptional and microscopic analyses of citrus stem and root responses to Candidatus Liberibacter asiaticus infection. PloS one. 2013;8(9):e73742. pmid:24058486; PubMed Central PMCID: PMC3772824.
  13. 13. Martinelli F, Reagan RL, Uratsu SL, Phu ML, Albrecht U, Zhao W, et al. Gene regulatory networks elucidating huanglongbing disease mechanisms. PloS one. 2013;8(9):e74256. pmid:24086326; PubMed Central PMCID: PMC3783430.
  14. 14. Jagoueix S, Bove J-M, Garnier M. The phloem-limited bacterium of greening disease of citrus is a member of the α subdivision of the Proteobacteria. Int J Syst Bacteriol. 1994;44(3):379–86. pmid:7520729
  15. 15. Xu MR, Liang MD, Chen JC, Xia YL, Zheng Z, Zhu Q, et al. Preliminary research on soil conditioner mediated citrus Huanglongbing mitigation in the field in Guangdong, China. Eur J Plant Pathol. 2013;137(2):283–93. pmid:WOS:000323747800010.
  16. 16. Livak KJ, Schmittgen TD. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2< sup>− ΔΔCT Method. methods. 2001;25(4):402–8. pmid:11846609
  17. 17. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4(2):249–64. pmid:WOS:000182894900007.
  18. 18. Adams MD, Kerlavage AR, Fleischmann RD, Fuldner RA, Bult CJ, Lee NH, et al. Initial assessment of human gene diversity and expression patterns based upon 83 million nucleotides of cDNA sequence. Nature. 1995;377(6547 Suppl):3–174. pmid:7566098
  19. 19. Du Z, Zhou X, Ling Y, Zhang ZH, Su Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;38:W64–W70. pmid:WOS:000284148900012.
  20. 20. Chen K, Xu C, Lou J, Zhang S, Ross G. Lipoxygenase in relation to the ripening and softening of Actinidia fruit. Acta phytophysiologica Sinica/Zhongguo zhi wu sheng li xue hui zhu bian. 1998;25(2):138–44.
  21. 21. Gao F, Zhou BJ, Li GY, Jia PS, Li H, Zhao YL, et al. A glutamic acid-rich protein identified in Verticillium dahliae from an insertional mutagenesis affects microsclerotial formation and pathogenicity. PloS one. 2010;5(12):e15319. pmid:21151869; PubMed Central PMCID: PMC2998422.
  22. 22. Machado CR, de Oliveira RL, Boiteux S, Praekelt UM, Meacock PA, Menck CF. Thi1, a thiamine biosynthetic gene in Arabidopsis thaliana, complements bacterial defects in DNA repair. Plant molecular biology. 1996;31(3):585–93. pmid:8790291.
  23. 23. Bekh-Ochir D, Shimada S, Yamagami A, Kanda S, Ogawa K, Nakazawa M, et al. A novel mitochondrial DnaJ/Hsp40 family protein BIL2 promotes plant growth and resistance against environmental stress in brassinosteroid signaling. Planta. 2013;237(6):1509–25. pmid:23494613; PubMed Central PMCID: PMC3664749.
  24. 24. Cho YS, Park SH, Kim CK, Oh KH. Induction of stress shock proteins DnaK and GroEL by phenoxyherbicide 2,4-D in Burkholderia sp. YK-2 isolated from rice field. Current microbiology. 2000;41(1):33–8. pmid:10919396.
  25. 25. Keith LM, Partridge JE, Bender CL. dnaK and the heat stress response of Pseudomonas syringae pv. glycinea. Molecular plant-microbe interactions: MPMI. 1999;12(7):563–74. pmid:10478477.
  26. 26. Kilstrup M, Jacobsen S, Hammer K, Vogensen FK. Induction of heat shock proteins DnaK, GroEL, and GroES by salt stress in Lactococcus lactis. Applied and environmental microbiology. 1997;63(5):1826–37. pmid:9143115; PubMed Central PMCID: PMC168475.
  27. 27. Liu JZ, Whitham SA. Overexpression of a soybean nuclear localized type-III DnaJ domain-containing HSP40 reveals its roles in cell death and disease resistance. The Plant journal: for cell and molecular biology. 2013;74(1):110–21. pmid:23289813.
  28. 28. Tanaka R, Rothbart M, Oka S, Takabayashi A, Takahashi K, Shibata M, et al. LIL3, a light-harvesting-like protein, plays an essential role in chlorophyll and tocopherol biosynthesis. P Natl Acad Sci USA. 2010;107(38):16721–5. pmid:WOS:000282003700061.
  29. 29. Tanaka A, Ito H, Tanaka R, Tanaka NK, Yoshida K, Okada K. Chlorophyll a oxygenase (CAO) is involved in chlorophyll b formation from chlorophyll a. Proceedings of the National Academy of Sciences. 1998;95(21):12719–23.
  30. 30. Han H, Gao S, Li B, Dong XC, Feng HL, Meng QW. Overexpression of violaxanthin de-epoxidase gene alleviates photoinhibition of PSII and PSI in tomato during high light and chilling stress. Journal of plant physiology. 2010;167(3):176–83. pmid:19767125.
  31. 31. Carmo-Silva AE, Salvucci ME. The Regulatory Properties of Rubisco Activase Differ among Species and Affect Photosynthetic Induction during Light Transitions. Plant Physiol. 2013;161(4):1645–55. pmid:WOS:000316987900006.
  32. 32. Fukayama H, Ueguchi C, Nishikawa K, Katoh N, Ishikawa C, Masumoto C, et al. Overexpression of rubisco activase decreases the photosynthetic CO2 assimilation rate by reducing rubisco content in rice leaves. Plant & cell physiology. 2012;53(6):976–86. pmid:22470057.
  33. 33. Mishra P, Dubey RS. Excess nickel modulates activities of carbohydrate metabolizing enzymes and induces accumulation of sugars by upregulating acid invertase and sucrose synthase in rice seedlings. Biometals: an international journal on the role of metal ions in biology, biochemistry, and medicine. 2013;26(1):97–111. pmid:23179408.
  34. 34. Wu AC, Ral J-P, Morell MK, Gilbert RG. New Perspectives on the Role of α-and β-Amylases in Transient Starch Synthesis. PloS one. 2014;9(6):e100498. pmid:24971464
  35. 35. Rejzek M, Stevenson CE, Southard AM, Stanley D, Denyer K, Smith AM, et al. Chemical genetics and cereal starch metabolism: structural basis of the non-covalent and covalent inhibition of barley β-amylase. Molecular BioSystems. 2011;7(3):718–30. pmid:21085740
  36. 36. Wu Z, Burns JK. A β-galactosidase gene is expressed during mature fruit abscission of ‘Valencia’orange (Citrus sinensis). Journal of experimental botany. 2004;55(402):1483–90. pmid:15208347
  37. 37. Zaare-Nahandi F, Hosseinkhani S, Zamani Z, Asadi-Abkenar A, Omidbaigi R. Delay expression of limonoid UDP-glucosyltransferase makes delayed bitterness in citrus. Biochemical and biophysical research communications. 2008;371(1):59–62. pmid:18407832
  38. 38. Moriguchi T, Kita M, Ogawa K, Tomono Y, Endo T, Omura M. Flavonol synthase gene expression during citrus fruit development. Physiologia plantarum. 2002;114(2):251–8. pmid:11903972
  39. 39. Wenkel S, Turck F, Singer K, Gissot L, Le Gourrierec J, Samach A, et al. CONSTANS and the CCAAT box binding complex share a functionally important domain and interact to regulate flowering of Arabidopsis. The Plant Cell Online. 2006;18(11):2971–84.
  40. 40. Pérez-Amador MA, Abler ML, De Rocher EJ, Thompson DM, van Hoof A, LeBrasseur ND, et al. Identification of BFN1, a bifunctional nuclease induced during leaf and stem senescence in Arabidopsis. Plant Physiology. 2000;122(1):169–80. pmid:10631260
  41. 41. Beck G, Coman D, Herren E, Ruiz-Sola MÁ, Rodríguez-Concepción M, Gruissem W, et al. Characterization of the GGPP synthase gene family in Arabidopsis thaliana. Plant molecular biology. 2013;82(4–5):393–416. pmid:23729351
  42. 42. Edstam MM, Blomqvist K, Eklof A, Wennergren U, Edqvist J. Coexpression patterns indicate that GPI-anchored non-specific lipid transfer proteins are involved in accumulation of cuticular wax, suberin and sporopollenin. Plant Mol Biol. 2013;83(6):625–49. pmid:WOS:000327093600008.
  43. 43. Itoh H, Tatsumi T, Sakamoto T, Otomo K, Toyomasu T, Kitano H, et al. A rice semi-dwarf gene, Tan-Ginbozu (D35), encodes the gibberellin biosynthesis enzyme, ent-kaurene oxidase. Plant molecular biology. 2004;54(4):533–47. pmid:15316288
  44. 44. Islam MS, Wang MH. Expression of dehydration responsive element-binding protein-3 (DREB3) under different abiotic stresses in tomato. BMB reports. 2009;42(9):611–6. pmid:19788864.
  45. 45. Zhang M, Liu W, Bi YP. [Dehydration-responsive element-binding (DREB) transcription factor in plants and its role during abiotic stresses]. Yi chuan = Hereditas / Zhongguo yi chuan xue hui bian ji. 2009;31(3):236–44. pmid:19273435.
  46. 46. Romeis T, Tang S, Hammond-Kosack K, Piedras P, Blatt M, Jones JD. Early signalling events in the Avr9/Cf-9-dependent plant defence response. Molecular plant pathology. 2000;1(1):3–8. pmid:20572944.
  47. 47. Lauvergeat V, Lacomme C, Lacombe E, Lasserre E, Roby D, Grima-Pettenati J. Two cinnamoyl-CoA reductase (CCR) genes from Arabidopsis thaliana are differentially expressed during development and in response to infection with pathogenic bacteria. Phytochemistry. 2001;57(7):1187–95. pmid:11430991.
  48. 48. Thévenin J, Pollet B, Letarnec B, Saulnier L, Gissot L, Maia-Grondard A, et al. The simultaneous repression of CCR and CAD, two enzymes of the lignin biosynthetic pathway, results in sterility and dwarfism in Arabidopsis thaliana. Molecular Plant. 2011;4(1):70–82. pmid:20829305
  49. 49. Mageroy MH, Tieman DM, Floystad A, Taylor MG, Klee HJ. A Solanum lycopersicum catechol‐O‐methyltransferase involved in synthesis of the flavor molecule guaiacol. The Plant Journal. 2012;69(6):1043–51. pmid:22103597
  50. 50. Zhu B. Catechol-O-Methyltransferase (COMT)-mediated methylation metabolism of endogenous bioactive catechols and modulation by endobiotics and xenobiotics: importance in pathophysiology and pathogenesis. Current drug metabolism. 2002;3(3):321–49. pmid:12083324
  51. 51. Do C-T, Pollet B, Thévenin J, Sibout R, Denoue D, Barrière Y, et al. Both caffeoyl Coenzyme A 3-O-methyltransferase 1 and caffeic acid O-methyltransferase 1 are involved in redundant functions for lignin, flavonoids and sinapoyl malate biosynthesis in Arabidopsis. Planta. 2007;226(5):1117–29. pmid:17594112
  52. 52. Trabucco GM, Matos DA, Lee SJ, Saathoff AJ, Priest HD, Mockler TC, et al. Functional characterization of cinnamyl alcohol dehydrogenase and caffeic acid O-methyltransferase in Brachypodium distachyon. BMC biotechnology. 2013;13:61. pmid:23902793; PubMed Central PMCID: PMC3734214.
  53. 53. Stracke R, De Vos RC, Bartelniewoehner L, Ishihara H, Sagasser M, Martens S, et al. Metabolomic and genetic analyses of flavonol synthesis in Arabidopsis thaliana support the in vivo involvement of leucoanthocyanidin dioxygenase. Planta. 2009;229(2):427–45. pmid:18998159
  54. 54. Chen W, Yin X, Wang L, Tian J, Yang R, Liu D, et al. Involvement of rose aquaporin RhPIP1; 1 in ethylene-regulated petal expansion through interaction with RhPIP2; 1. Plant molecular biology. 2013;83(3):219–33. pmid:23748738
  55. 55. Choi DS, Hwang IS, Hwang BK. Requirement of the cytosolic interaction between PATHOGENESIS-RELATED PROTEIN10 and LEUCINE-RICH REPEAT PROTEIN1 for cell death and defense signaling in pepper. The Plant Cell Online. 2012;24(4):1675–90.
  56. 56. Bar M, Sharfman M, Ron M, Avni A. BAK1 is required for the attenuation of ethylene-inducing xylanase (Eix)-induced defense responses by the decoy receptor LeEix1. The Plant journal: for cell and molecular biology. 2010;63(5):791–800. pmid:20561260.
  57. 57. Rotblat B, Enshell-Seijffers D, Gershoni JM, Schuster S, Avni A. Identification of an essential component of the elicitation active site of the EIX protein elicitor. The Plant journal: for cell and molecular biology. 2002;32(6):1049–55. pmid:12492845.
  58. 58. Pontier D, Balague C, Bezombes-Marion I, Tronchet M, Deslandes L, Roby D. Identification of a novel pathogen-responsive element in the promoter of the tobacco gene HSR203J, a molecular marker of the hypersensitive response. The Plant journal: for cell and molecular biology. 2001;26(5):495–507. pmid:11439136.
  59. 59. Magnan F, Ranty B, Charpenteau M, Sotta B, Galaud JP, Aldon D. Mutations in AtCML9, a calmodulin-like protein from Arabidopsis thaliana, alter plant responses to abiotic stress and abscisic acid. The Plant journal: for cell and molecular biology. 2008;56(4):575–89. pmid:18643966.
  60. 60. Hong MJ, Kim DY, Lee TG, Jeon WB, Seo YW. Functional characterization of pectin methylesterase inhibitor (PMEI) in wheat. Genes & genetic systems. 2010;85(2):97–106.
  61. 61. An SH, Sohn KH, Choi HW, Hwang IS, Lee SC, Hwang BK. Pepper pectin methylesterase inhibitor protein CaPMEI1 is required for antifungal activity, basal disease resistance and abiotic stress tolerance. Planta. 2008;228(1):61–78. pmid:18327607
  62. 62. Park YH, Choi C, Park EM, Kim HS, Park HJ, Bae SC, et al. Over-expression of rice leucine-rich repeat protein results in activation of defense response, thereby enhancing resistance to bacterial soft rot in Chinese cabbage. Plant cell reports. 2012;31(10):1845–50. pmid:22717673
  63. 63. Bhadauria V, Banniza S, Vandenberg A, Selvaraj G, Wei Y. Alanine: Glyoxylate aminotransferase 1 is required for mobilization and utilization of triglycerides during infection process of the rice blast pathogen, Magnaporthe oryzae. Plant signaling & behavior. 2012;7(9):1206–8.