{'fraction_non_alphanumeric': 0.03880739542162922, 'fraction_numerical': 0.03673837502321017, 'mean_word_length': 5.048451788865715, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Transcriptome analysis of Hpa1Xoo transformed cotton revealed constitutive expression of genes in multiple signalling pathways related to disease resistance

The transcriptome profile in leaves and roots of the transgenic cotton line T-34 expressing hpa1Xoo from Xanthomonas oryzae pv. oryzae was analysed using a customized 12k cotton cDNA microarray. A total of 530 cDNA transcripts involved in 34 pathways were differentially expressed in the transgenic line T-34, in which 123 differentially expressed genes were related to the cotton defence responses including the hypersensitive reaction, defence responses associated with the recognition of pathogen-derived elicitors, and defence signalling pathways mediated by salicylic acid, jasmonic acid, ethylene, auxin, abscicic acid, and Ca2+. Furthermore, transcripts encoding various leucine-rich protein kinases and mitogen-activated protein kinases were up-regulated in the transgenic line T-34 and expression of transcripts related to the energy producing and consuming pathway was also increased, which suggested that the enhanced metabolism related to the host defence response in the transgenic line T-34 imposed an increased energy demand on the transgenic plant.


Introduction
Cotton (Gossypium spp.) is one of the most important textile fibre crops worldwide. A great many fundamental and applied studies have been conducted on this crop in the areas of weed and pest management due to the large and profitable market for growers. In the past two decades, genetically modified cotton has been developed in many countries, which has reduced reliance on pesticides for the production of this crop by ;80%, compared with the areas where conventional cotton varieties were grown (Phipps and Park, 2002). In China, genetically engineered cotton expressing d-endotoxins (Cry proteins) from Bacillus thuringiensis (Bt) was developed in 1997 and then adopted by 95% of cotton growers in northern China (Wu et al., 2008).
Although Bt cotton reduced the overall need for the insecticide spray initially, it was recently discovered that the wide application of Bt cotton in China had resulted in an increase in the mirid bug population on cotton due to the decrease of pesticide usage on this crop (Lu et al., 2010). Therefore, it is imperative to develop new sources of transgenes for genetically modified cotton that can provide a broader spectrum of resistance for pest management on cotton and increase the diversity of genetically engineered cotton in China.
Harpin is a heat-stable, glycine-rich and acidic protein secreted by Gram-negative plant pathogenic bacteria into the plant intercellular space through the type III secretion system (T3SS) (Wei et al., 1992;Perino et al., 1999). To date, a few harpins or harpin-like proteins, such as hrpNEa from Erwinia amylovora (Kim and Beer, 1998;Dong et al., 1999), hrpZ Pss from Pseudomonas syringae pv. syringae (He et al., 1993;Strobel et al., 1996), and hpa1 Xoo from Xanthomonas oryzae pv. oryzae  have been characterized and their roles in the non-host interaction have been investigated.
Harpins can induce hypersensitive response (HR) and systemic acquired resistance in non-host plants through exogenous application (Wei and Beer, 1993). In addition, they are capable of inducing many other host responses such as enhanced growth and drought tolerance in tobacco and Arabidopsis (Dong et al., 1999;Dong et al., 2004;Dong et al., 2005). Furthermore, harpins regulate the plasma membrane ion channels that are considered to be putative components of the signalling pathway related to the host defence response (El-Maarouf et al., 2001). Harpins can also induce the potassium efflux and rapid inhibition of ATP synthesis that result in plasma membrane depolarization and growth medium alkalinization in various plant species (Popham et al., 1995;Hoyos et al., 1996;Xie and Chen, 2000). More recently, it was found that HrpN of E. amylovora interacted with HrpN-interacting protein from Malus (HIPM) and its orthologue, AtHIPM, of Arabidopsis. Both genes contained functional signal peptides associated with the plasma membrane (Oh and Beer, 2007). It has been hypothesized that these plant phenotypes induced by harpins are related to the change in expression of genes involved in various signalling pathways, such as salicylic acid (SA)-, jasmonate (JA)-, and ethylene (ET)-mediated signalling pathways (Somssich and Hahlbrock, 1998;Dong et al., 1999).
In addition to exogenous application, it was reported that harpins expressed in transgenic plants also induced the host defence response . To date, several hrp genes have been transformed into plant species, including tobacco, rice, and Arabidopsis, and enhanced resistance to various bacterial and fungal pathogens was obtained in the transgenic plants Sohn et al., 2007). Several pathogenicity-related genes (pr), such as pr-1a, pr-1b, pr-2, pr-3, and Chia5, and genes related to the production of ET (e.g. NT-EFE26, NT-1A1C, DS321, NTACS1, and NT-ACS2) were up-regulated in transgenic plants expressing harpins Sohn et al., 2007;Shao et al., 2008). Nevertheless, the up-regulation of npr-1 (non-expressor of pr genes) in hrp-transformed plants was variable, which suggested that expression of npr-1 differed depending on the host and the origin of the hrp gene used in the transformation Sohn et al., 2007).
Although harpins induced the HR in non-host plants through exogenous application (He et al., 1993), no visible HR was found in the transgenic cotton line T-34 harbouring a constitutively expressed hpa1 Xoo in our previous study (Miao et al., 2010). The hpa1 Xoo -transformed cotton line T-34 showed improved resistance to Verticillum and Fusarium wilt caused by Verticillium dahliae and Fusarium oxysporum f.sp. vasinfectum, respectively. In addition, the oxidative burst and up-regulation of several key defencerelated genes were observed in transgenic T-34 in response to infection caused by V. dahliae, which suggested that transformation of cotton with hpa1 Xoo conferred enhanced defence response on pathogens through a priming mechanism. Since the mechanism of harpin-mediated plant responses is still largely unknown, the possibility of genome-wide modification of genetically modified plants or plants treated by harpins remains to be investigated. With the recent development of gene chip technology, it is possible to investigate the defence responses and signalling pathways involved in plants transformed with hpa1 Xoo on a genome-wide scale (Kim et al., 2006;Chibucos et al., 2009). In the present study, we investigated the pattern of global gene expression in the transgenic cotton line T-34 expressing hpa1 Xoo , and in the wild-type receptor (Z35) using a customized cotton 12k cDNA microarray. This study is of importance for understanding the function of hpa1 Xoo in genetically modified cotton and its effects on the cotton defence signalling pathway.

Plant and fungal materials
The transgenic cotton line T-34 was developed previously through the genetic transformation of the susceptible cotton variety Zhong-Mian-35 (Gossypium hirsutum L., abbreviated Z35) with hpa1 Xoo derived from X. oryzae pv. oryzae using a modified Agrobacteriummediated method (Miao et al., 2010). Homozygous transgenic lines from T1 to T6 were screened using kanamycin resistance, PCR analysis, Southern and Northern analysis as previously described (Miao et al., 2010). Three T6 progeny of the transgenic T-34 used in this study were grown in a sterile incubator (MLR-351; Sanyo Electric Co., Ltd, Japan) at 28°C with a photoperiod of 12 h under incandescent light. Cotton leaves were collected at the true leaf stage and stored in liquid nitrogen prior to RNA extraction. Wildtype cotton Z35 was used as the negative control.
V. dahliae strain (V bps ) was provided by Dr Ling Lin (Jiangsu Academy of Agricultural Science, China). V. dahliae was maintained on potato dextrose agar (PDA) at 25°C. The method of Joost et al. (1995) was used to prepare the inoculum and the V. dahliae conidia suspension was adjusted to a concentration of 1 3 10 7 conidia/ml. Two millilitres of conidia suspension was used to inoculate leaves of T-34 and Z35 by freshly cutting at the four to five leaf stage or dipping petioles for 3 h.
Microarray procedure and data analysis Total RNA samples were extracted from 1 g of cotton leaves and roots using the Plant Total RNA Extract Kit (Autolab Biotech, Beijing, China) according to the manufacturer's instructions. RNA samples were collected from three biological replicates. The 12k cDNA microarray was conduced at CapitalBio Corp. (Beijing, China) using the method described by Shi et al. (2006). Briefly, PCR products from 11 236 cotton unigene expressed sequence tags (ESTs) were printed on to amino-silanized glass slides in triplicate for each PCR product. Total RNA (5 lg) was used to synthesize cDNA in an in vitro transcription reaction and fluorescently labelled using Klenow enzyme (Promega, Beijing, China). After hybridizations, arrays were scanned with a confocal LuxScanä scanner (CapitalBio, Beijing, China) and analysed using the software LuxScanTM 3.0 (CapitalBio, Beijing, China). For extraction of data from the individual channels, faint spots with intensity <400 units after subtraction of the background were removed. A space-and intensity-dependent normalization method based on a LOWESS program (Yang et al., 2002) was used. The gene annotation and pathway identification were performed according to the method described by Shi et al. (2006).

Real-time RT-PCR for gene expression
Three hours after inoculation, the complete laminae of cotton leaves at the V4 stage (four-leaf stage) were harvested from T-34 and Z35 and frozen in liquid nitrogen until the RNA extraction. RNA was extracted from cotton leaves using the RNAiso kit for polysaccharide-rich plant tissue (TaKaRa Ò Biotechnology, Dalian, China). The concentration of RNA was quantified using a biophotometer (Eppendorf AG, Hamburg, Germany). cDNA was synthesized using a PrimeScriptTM RT-PCR Kit (TaKaRa Ò Biotechnology, Dalian, China). Two-step real-time RT-PCR was performed on an ABI PRISM 7000 (ABI, Foster City, CA, USA) according to the procedure optimized for the SYBR Premix Ex TaqTM kit (TaKaRa Ò Biotechnology, Dalian, China). EF-1a, a conserved plant housekeeping gene , was used as the internal control to normalize the level of expression. All real-time RT-PCRs were performed in duplicate and three biological replicates were included. The expression data were normalized to EF-1a using the DDCT method described by Livak and Schmittgen (2001).

Results
The selection of plant materials for microarray analysis Three plants from T6 progeny of the transgenic T-34 were randomly selected. The presence and the constitutive expression of hpa1 Xoo in the selected T-34 were verified using PCR, Southern and Northern blot analysis. In all three plants, bands representing hpa1 Xoo insert, the 35s promoter, and the NOS terminator (420, 310, and 180 bp in length, respectively) were found in the PCR analysis with hpa1 Xoo -specific primers. In Southern blot analysis, three positive bands (4, 6, and 10 kb in length) were detected in all three T-34 progeny indicating the insertion of hpa1 Xoo at multiple chromosomal locations. The constitutive expression of hpa1 Xoo in selected T-34 plants was observed in the Western blot using the polyclonal antibody raised against purified hpa1 Xoo (data not shown).

The construction and verification of the cotton cDNA microarray analysis
Of 12 233 unigene ESTs initially selected, 11 236 ESTs were successfully amplified through PCR amplifications and printed onto aminosilane slides. These ESTs were used to construct a cDNA microarray for the identification of genes specifically or preferentially expressed in transgenic T-34. Differentially expressed genes (DEGs) that showed upregulation (2-fold increase in expression compared with expression in wild-type Z35) or down-regulation (2-fold decrease in expression compared with expression in wild-type Z35) in all three replicates of transgenic T-34 (corrected for the false discovery rate (FDR), P < 0.01) were considered as related to the transformation of hpa1 Xoo ( Figure. 1).
The quality of microarray data was verified based on the value of correlation coefficient (r values) and the swap-dye experiment. The correlation coefficient, which measured the biological reproducibility, was calculated from each replicate ( Table 1). The value from the swap-dye experiment measured the technical reproducibility. After the self-hybridization of Fig. 1. Hierarchical clustering of DEGs in leaves and roots of transgenic cotton T-34 compared with wild-type Z35 in the microarray analysis. Each column represents a single biological replicate and each row represents a differentially expressed probe set. L1, L2, and L3 represent biological replicates from leaves and R1, R2, and R3 represent biological replicates from roots. S/S represented self-to-self of Z35. The signal ratios were shown in a red-green colour scale, where red indicated up-regulation and green indicated down-regulation. Cy3-and Cy5-labelled probes prepared using the same RNA sample prepared from leaves of wild-type Z35 (V4 stage), all but 13 data points were scattered inside the 2-fold line, which indicated that the microarray was precisely executed (data not shown). In addition, we applied a spaceand intensity-dependent normalization based on a LOWESS program (Yang et al., 2002). Forty internal controls were deployed to obtain evenly distributed signal intensities (see Supplementary Table 1, at JXB Online).
To validate results from the microarray analysis, nine genes representing a wide range of relative fluorescencesignal intensity and different functions were selected from 530 DEGs found in leaves of T-34 and subjected to realtime RT-PCR analysis ( Table 2). The expression of eight genes analysed with real-time RT-PCR was consistent with the data obtained from the microarray analysis (Table 2).

Global cell expression patterns in the transgenic line T-34 and wild-type receptor Z35
The microarray data were deposited in the GEO database of the National Center for Biotechnology (accession number GSE20446). In leaves and roots of transgenic T-34, a total of 1552 and 672 DEGs, respectively, were related to the hpa1 Xoo transformation with the exhibited ratio of >2.0 (2-fold increase in expression) or <0.5 (2-fold decrease in expression) in the microarray (Fig. 2). In leaves of T-34, 530 DEGs related to the hpa1 Xoo transformation were found in all three replicates, in which 229 transcripts (43.2%) were up-regulated and 301 transcripts (56.7%) were downregulated. In comparison, 65 DEGs related to the hpa1 Xoo transformation were identified in three root replicates, in which 49 (74.5%) transcripts were up-regulated and 16 (24.5%) transcripts were down-regulated. The comparison of DEGs in leaves and roots from transgenic T-34 and wildtype Z35 showed that 36 DEGs were up/down-regulated in both leaves and roots of transgenic T-34. Of these DEGs, 31 DEGs were up-regulated whereas only one DEG was down-regulated. Four DEGs were up-/down-regulated conversely in leaves and roots of transgenic T-34 ( Fig. 2 and Table 3).
Five hundred and thirty DEGs found in leaves of transgenic T-34 were annotated through the blast against Kegg (www.genome.ad.jp/kegg), Biocyc (www.biocyc.org), and Kobas (kobas.cbi.pku.edu.cn:8080) databases (Fig. 3). The largest class of annotated DEGs was transcription factors (10.47%) followed by DEGs encoding various enzymes (9.69, and proteins involved in the protein degradation process (6.39%). In addition, 6.17% of DEGs were related to responses induced by biotic or abiotic stress and 2.42% of DEGs were involved in hormone biosynthesis. Among 36 DEGs that were differentially expressed in leaves and roots of transgenic T-34, five DEGs (CM105C02, CM105B08, CM105F10, CM101F11, and CM049G04) were Housekeeping gene used as internal reference related to responses induced by biotic or abiotic stress and two DEGs (CM104G02 and CM105G10) were transcription factors (Table 3). Twenty-one DEGs were involved in cytological and biochemical processes related to the defences response, such as programmed cell death (PCD), phenylpropanoid and flavonoid biosynthesis, lignin biosynthesis, pathogenesisrelated proteins, and NAD/NADH phosphorylation and dephosphorylation (Table 4). In addition, a large group of DEGs involved in the signalling transduction cascade mediated by protein kinases, G-proteins, and calmodulinlike proteins were also identified (Table 5).

Differentially expressed transcripts related to the defence response in leaves of transgenic T-34
Five up-regulated DEGs were related to PCD. The upregulation of CM105A03 was the most significant followed by CM103E04 and CM025G05, which were increased by 6.74-, 2.45-and 2.09-fold, respectively. CM105A03 encoded an E3 ligase, an inhibitor of PCD in plants, and CM103E04 had sequence similarity to an aspartic protease gene involved in PCD by hydrolysing important cell components and/or activating other proteinases. CM025G05 was similar to the hypersensitivity-related gene hsr201 found in Nicotiana tabacum (Table 4).
Among five DEGs related to phenylpropanoid biosynthesis, the up-regulation of CM107B04, which encoded a phenylalanine ammonia lyase (PAL) catalysing the initial reaction of the phenylpropanoid biosynthetic pathway, was the most significant (2.63-fold) ( Table 4). In addition, expression of CM048C05 was increased by 2.48-fold in transgenic T-34. CM048C05 was nearly identical to a quinate hydroxycinnamoyl transferase (HCT), which synthesized hydroxycinnamoylesters in the phenylpropanoid pathway. In comparison, two DEGs (CM092H07 and CM099F12) were down-regulated in leaves of transgenic T-34. CM092H07 was similar to a 3-hydroxy-3-methylglutaryl-CoA reductase, which was the rate-limiting enzyme in the mevalonate pathway related to isoprenoid biosynthesis. CM099F12 encoded a prephenate dehydratase participating in phenylalanine, tyrosine, and tryptophan biosynthesis.
Four DEGs related to flavonoid biosynthesis were also upregulated in leaves of transgenic T-34 (2.09-to 3.63-fold increase in expression). Three DEGs (CM087F05, CM100D01, and CM113E06) encoded products similar to chalcone synthase. Chalcone synthase is important in the production of chalcones, a class of organic compounds commonly found in plants related to natural defence mechanisms. CM092E03 encoded a flavanone 3-dioxygenase. Furthermore, expression of three DEGs (CM072F03, CM089D11, and CM078C09) related to the biosynthesis of lignin was also increased in leaves of transgenic T-34 (5.21-, 2.39-, and 2.9-fold, respectivley). CM072F03 and CM089D11 were both similar to the gene encoding cinnamoyl-CoA reductase whereas CM078C09 had similarity to a gene encoding dihydroflavonol 4-reductase. Both of these two genes are very important for lignin biosynthesis in plants (Table 4).
Surprisingly, three DEGs with similarities to previously described pathogenesis-related (pr) genes were down-regulated in transgenic T-34. CM011H11 and CM113A02, which encoded products similar to a basic endochitinase (PR-3), were down-regulated by 0.44-and 0.29-fold, respectively. A decrease in expression was also found for CM072H01, which encodes a putative peroxidase (PR-9) ( Table 4).

Differentially expressed transcripts related to signal transduction in leaves of transgenic T-34
Receptor protein kinases play a fundamental signalling role in plant growth, plant development and host responses to biotic and abiotic stress. In the present study, 11 DEGs with similarities to various receptor protein kinases were identified in transgenic T-34 in which the leucine-rich repeat receptor kinases represented the largest group (CM113G02, CM030A02, CM057G05, CM023F07, CM019F07, CM089A03, CM076F10, CM093F12, CM081B01, and CM105D02). Additionally, expression of two DEGs (CM079E09 and CM057G05) similar to lectin protein kinases was also up-regulated 3.09-and 2.3-fold, respectively. In comparison, CM081B01, which encoded a mitogen-activated protein kinase (MAPK), was down-regulated (0.41-fold) ( Table 5). It has been reported that MAPKs relay signals from receptors and activate downstream defence responses in the host defence response (Desikan et al., 1999;Dievart and Clark, 2004).
Ethylene (ET), auxin, and abscisic acid (ABA) are three major plant hormones with important roles in signalling related to plant growth and the defence response. Three DEGs (CM105E12, CM038C12, and CM058G10), which encoded products similar to ET-responsive ERF/AP2 transcription factors, were up-regulated in transgenic T-34 (Table 5). The increase in expression of CM105E12 was the highest (8.83-fold) followed by CM038C12 and CM058G10, which were increased by 2.16-and 2.12-fold, respectively. CM111F04, CM024E02, and CM024A04 were all annotated as auxin-responsive transcription factors and their expression was increased 2.0-to 3.5-fold in leaves of transgenic T-34. In addition, one DEG (CM105G04) with similarity to a ROP GTPase related to an ABA-mediated signalling pathway and two DEGs similar to myb transcription factors were also up-regulated in transgenic T-34 (Table 5). Myb family transcription factors play regulatory roles in plant developmental processes and host responses to one or multiple types of hormone and stress treatment.
Small GTP-binding proteins are another class of signalling molecules that play a critical role in the regulation of a range of cellular processes including growth, differentiation, and intracellular transportation. In our study, three DEGs (CM019F05, CM105C07, and CM028G08) were identified as Ras-related GTP-binding proteins (Table 5). Although expression of CM019F05 was increased in transgenic T-34 (2.69fold), CM105C07 and CM028G08 were down-regulated. Five DEGs involved in the intracellular signalling cascade were also differentially regulated in transgenic plants. CM025C08, CM006A06, CM052F01, and CM092E09 encoded calciumbinding proteins and CM033E11 encoded a phosphoinositidespecific phospholipase.

Differentially expressed transcripts related to energy production in leaves of transgenic T-34
Since plant defence responses are often energy consuming, up-or down-regulation of DEGs related to the defence response, the biosynthesis of secondary metabolites, transcription factors and cellular signalling molecules may alter the state of the energy production in the transgenic plant. In our study, several DEGs related to various energy-producing pathways were up-regulated (>2.0-fold increase in expression) ( Table 6). CM103E07 (phosphopyruvate hydratase) and CM019F08 (peroxisomal biogenesis factor) were related to glycolysis. CM110G09 (aconitate hydratase)  and CM024A06 (isocitrate dehydrogenase) were involved in the tricarboxylic acid (TCA) cycle. CM102G09, CM006G10, and CM032E11, which were annotated as acetyl-CoA carboxylase, biotin carboxyl carrier protein, and the malic enzyme, were important for fatty acid biosynthesis. Furthermore, three DEGs (CM031A04, CM101H09, and CM097C03), related to NAD/NADH phosphorylation and dephosphorylation, were also up-regulated in transgenic T-34.

Enhanced expression of defence-related genes from T-34 after challenge by V. dahliae
Based on the result of microarray, six DEGs including CM011H11 (basic chitinase), CM048G (harpin-induced family protein), CM021D04 (S-adenosylmethionine synthetase), CM031A04 (NADH-ubiquinone oxidoreductase), CM107B04 (PAL2), and CM025C08 (calmodulin-like protein) were selected from different pathways and their expression in transgenic T-34 in response to V. dahliae infection was analysed using real-time RT-PCR analysis ( Fig. 4 and Table 3). All six DEGs were significantly upregulated in transgenic T-34, a >100 to 1000-fold increase in expression, after the transgenic T-34 plants were challenged with V. dahliae.

Discussion
Microarray is a useful technology to measure the transcriptome profile, which has been used successfully in various plant species to examine the host response to abiotic and biotic stress (Clarke and Zhu, 2006;Miyama and Tada, 2008;Van Hoewyk et al., 2008). In the present study, we used the microarray to investigate the transcriptome profile in leaves and roots of transgenic T-34 and wild-type Z35. Among 11 236 unigene ESTs included in the 12k cDNA microarray, 4.7% and 0.57% of ESTs were differentially expressed in leaves and roots of transgenic T-34, respectively. The ratios of DEGs to total ESTs included in the microarray chip are similar to that found in transgenic Arabidopsis overexpressing a stress response brassinosteroid receptor gene (BR1) (Kim et al., 2010) but significantly higher than that found in transgenic rice transformed with choline oxidase (0.032%) (Su et al., 2006), and with alanine aminotransferase (0.11%) (Shrawat et al., 2008). Although it is possible that this difference is merely due to the relatively larger number of transcripts included in gene chips used in previous studies, the higher ratios of DEGs found in transgenic T-34 could indicate that a larger portion of genes was differentially regulated in transgenic T-34. Five hundred and thirty and 65 DEGs related to hpa1 Xoo transformation were found in leaves and roots, respectively, of transgenic T-34, which indicateded that the transcriptome profile was altered not only in leaves but also in roots of transgenic T-34. Since no tissue specific promoter was added in front of the hpa1 Xoo transgene in the transformation and our expression analysis indicated that hpa1 Xoo was expressed in both leaves and roots of T-34 (data not shown), it is unlikely that the change in the transcriptome profile in different tissues of T-34 was due to the systemic effect of hpa1 Xoo . Although 36 common DEGs were found in leaves and roots of transgenic T-34, which indicated that cells from different tissues of T-34 had a similar reaction to the transformation of hpa1 Xoo , some differences were found between DEGs from leaves and roots of T-34. For examples, transcripts CM049G04, CM110G08, CM111E02, and CM022B06 were transversely regulated in leaves and roots of transgenic T-34. Some of DEGs found in leaves were not present in roots of T-34 (data not shown). These differences suggested the presence of certain tissue specificity in cotton in response to the hpa1 Xoo transformation.
In our study, several DEGs (e.g. CM025G05, CM009E10, and CM048G07) with similarities to genes encoding HRrelated proteins (HSR203 and proteins with signature TIR-NBS-LRR domains) were up-regulated in transgenic T-34. The up-regulation of these HR-related genes in transgenic T-34 is in agreement with our previous finding that micro-HRs were observed in leaves of transgenic T-34 in the absence of the pathogen (Miao et al., 2010). Up to date, three types of hrp-induced HR have been described including the visible HRs in response to foliar infiltration (He et al., 1993), microscopic HRs induced by foliar spray , and up-regulated HR-related marker genes in the absence of visible HR phenotypes in transgenic plants expressing harpins . This difference in HR phenotypes suggests that hrp-mediated HRs in hrptransformed transgenic plants differ from plants treated exogenously with harpins. It is interesting that CM103E04, which encoded an aspartic protease, was up-regulated in transgenic T-34 2.45-fold, compared with its expression in the wild-type Z35. Ge et al. (2005) previously reported that the aspartic protease functioned as an anti-cell death component in the reproduction and embryogenesis of Arabidopsis. Similarly, a rice nucellin gene (OsAsp1), which encoded an aspartic protease, was found to play an important role in PCD of the rice nucleolus (Bi et al., 2005). It is possible that the up-regulation of DEGs encoding anti-cell death inhibitor-like products in transgenic T-34 indicates the counter-response of the hpa1 Xoo -transformed plants to reduce the excess rate of the cell death since hpa1 Xoo is an effector protein inducing hypertensive reaction in non-host plants. Nevertheless we cannot rule out the possibility that the absence of visible HR phenotypes in transgenic T-34 is due merely to the dosage effect of hpa1 Xoo expressed in the transgenic plant (Miao et al., 2010).
In addition to DEGs related to HRs, DEGs involved in the biosynthesis of secondary metabolites were also upregulated in transgenic T-34. For example, CM107B04, which encoded a phenylalanine ammonia-lyase (pal), was up-regulated in transgenic T-34. The expression of three transcripts (CM087F05, CM100D01, and CM113E06) similar to a chalcone synthase (Chs) was also increased in transgenic T-34. Both Pal and Chs are considered as the classical marker genes for the defence response since end products of phenylpropanoid and flavonoid biosynthesis, such as flavonoid phytoalexins and lignin, often play diverse roles in host responses to biotic or abiotic stimuli (Sewalt et al., 1997). Accordingly, DEGs related to flavonoid and lignin biosynthesis were up-regulated in transgenic T-34. The up-regulation of DEGs related to the biosynthesis of defencerelated secondary metabolites suggested the constitutive activation of the host defence response in transgenic T-34. This result is in line with our previous finding that transgenic T-34 has enhanced resistance to a range of pathogens such as F. oxysporium and V. dahliae (Miao et al., 2010). It is surprising that several DEGs with similarities to basic endochitinase (PR-3) and peroxidase (PR-9) were down-regulated in transgenic T-34. PR proteins can either promote the degradation of fungal cell wall or release endogenous elicitors from the plant cell wall, which further stimulate the defence response (Baron and Zambryski, 1995). It was unlikely that the down-regulation of DEGs similar to pr proteins in transgenic T-34 was due to the suppression of the host defence response, since DEGs related to HR, flavonoid, and lignin biosynthesis were upregulated in transgenic T-34. In addition, our previous study showed that several key defence-related genes were also up-regulated in transgenic T-34 in response to inoculation with V. dahliae (Miao et al., 2010). Since only specific isoforms of basic chitinases and peroxidases exhibited antifungal activity (Sela-Buurlage et al., 1993), these pr-like DEGs may be not be related to the host defence response. For example, basic endochitinases are also expressed in tobacco and tomato during flower formation. Furthermore, several DEGs with similarity to endo-b-1,4-glucanases and glucan 1,3-b-glucosidases were also down-regulated in transgenic T-34 (data not shown). It is possible that the down-regulation of these DEGs is the result of a shift in carbohydrate metabolism in the transgenic plant due to the increase in secondary metabolite biosynthesis and energy production.
In addition to the up-regulation of DEGs related to the defence response, a number of DEGs involved in different signalling pathways were also up-regulated in transgenic T-34. For example, CM105A03, which encodes an E3 ubiquitin ligase, was up-regulated by 6.74-fold in transgenic T-34. There is emerging evidence that E3 ligase is important in the signalling of the plant immune system. For example, Yang et al. (2006) reported that E3 ligases were positive regulators of cell death and the host defence response across the Solanaceae and Brassicaceae families. Similarly, Abramovitch et al. (2006) found that the ubiquitin liagse activity of bacterial TypIII effector AvrPtoB suppressed plant cell death and immunity. To date, the regulation of defence gene expression by ET and JA has been linked to ubiquitination, which supports the hypothesis that multiple pathogen perception systems converge on common ubiquitination-based signalling pathways (Devoto et al., 2003). The up-regulation of CM105A03 in transgenic T-34 suggests that the ubiquitination-based signalling pathway is important for hpa1 Xoo -mediated signalling in transgenic cotton T-34.
It has been reported that various plant hormones are involved in hrp-regulated plant PCD (Dong, 1998). For example, the ET-mediated signalling pathway was associated with pathogen and insect resistance in Arabidopsis treated with harpin Ea from E. amylovora . In our study, three DEGs with similarities to the ET-responsive ERF/AP2 transcription factor were up-regulated in transgenic T-34. The ERF/AP2 transcription factor binds specifically to the GCC box present in many pr genes, to modulate their expression and participate in disease-resistance signalling pathways (Brown et al., 2003). The up-regulation of ERF/AP2 transcription factors has been previously reported in the defence response of Arabidopsis to infection caused by P. syringae pv. tomato DC3000 (avrRpt2) (Onate-Sanchez et al., 2007). Sasaki et al. (2007) reported that overexpression of ERF genes resulted in broad-spectrum resistance in transgenic plants. Dong et al. (2004) found that harpin activated ET signalling in Arabidopsis through EIN2 and EIN5, which conferred insect resistance and plant root growth in Arabidopsis. Nevertheless, two ET-responsive proteins, EIN2 and EIN5, were not differentially expressed in transgenic T-34 (data not shown). This discrepancy suggests that the EIN gene family may be not involved in the signalling pathway in transgenic T-34 in response to hpa1 Xoo transformation.
The plant hormone auxin regulates diverse aspects of plant reactions, including plant growth and development, and HRs (Dharmasiri and Estelle, 2004;Badescu and Napier, 2006). Gopalan (2008) reported that PCD initiated by harpin Ea was reversed by auxins without affecting the expression of marker genes related to local and systemic immunity. So, an interesting question is whether enhanced expression of auxin can be observed in genetically modified cotton. In our study, several auxin-related DEGs were upregulated in transgenic T-34 which supported the hypothesis that the auxin-mediated signalling cascade played an important role in the hpa1 Xoo -mediated response in transgenic T-34.
ABA, as an essential signal molecule modulating the plant response to abiotic and biotic stress, has divergent effects on defence responses (Asselbergh et al., 2004;Adie et al., 2007). The results of microarray analysis herein provided further evidence of the complexity of hormonally controlled signalling networks. For example, DEGs encoding products similar to myb transcription factors and GTPases responsive to ABA stimulus were up-regulated in transgenic T-34. These genes have been reported to be also responsive to other plant hormones such as SA and JA (Dong et al., 2005). Additionally, DEGs with similarities to GTP-and calucium-binding proteins related to the intercellular signalling cascade were also up-regulated in transgenic T-34. Although the mechanism for the upregulation of these signalling-related products in transgenic T-34 is still not clear, our previous study showed that the transformation of hpa1 Xoo conferred enhanced resistance in transgenic cotton through a priming mechanism (Miao et al., 2010). It is possible that the up-regulation of these signalling-related transcripts in hpa1 Xoo -transformed cotton increased the sensitivity of the transgenic plant to pathogen infection. This hypothesis is supported by the observation that multiple DEGs with similarities to leucine-rich repeat (LRR) protein were up-regulated in transgenic T-34. It is well known that LRR proteins are not only essential for the activation of defence genes after recognition between the plant and the product of an avirulence gene from the pathogen but also play important roles in various actions after ligand recognition. Nevertheless, further study is required to better understand the role and interplay of multiple hormones and other signalling molecules involved in pathways regulated by hpa1 Xoo in the transgenic plant.
It is interesting that several DEGs related to energy usage were also up-regulated in transgenic T-34. Although the effects of the up-regulation of DEGs in relation to the energy production and consumption in transgenic cotton is still not clear, our previous study (Miao et al., 2010) showed that the hpa1 Xoo -transformed transgenic line T-34 shared similar phenotypic characteristics, such as leaf morphology and fibre quality. The height of the T-34 line was lower than wild-type Z35 before the flowering stage. Nevertheless, no significant difference was found between the height of T-34 and Z35 after the flowering stage. Further study will be required to investigate whether the gain of cotton disease resistance resulting from the transformation by hpa1 Xoo will outweigh the possible decreased yield due to the change in the defence energy cost.
In summary, comprehensive information has been provided on transcriptome analysis of the hpa1 Xoo -transformed transgenic cotton T-34. Transcripts related to the defence response, secondary metabolite biosynthesis, and various signalling pathways were differentially expressed in the transgenic plant in response to hpa1 Xoo transformation. In addition, transcripts related to the energy producing pathways were also up-regulated in transgenic plants, which suggested that a high energy demand was imposed on transgenic cotton to support the cellular energy consumption resulting from the induction of multiple metabolic and host defence responses. This is somewhat similar to race non-specific resistance mediated by Lr34 as reported by Bolton et al. (2008).