Characterization by Suppression Subtractive Hybridization of Transcripts That Are Differentially Expressed in Leaves of Anthracnose-Resistant Ramie Cultivar.

For the purpose of screening putative anthracnose resistance-related genes of ramie (Boehmeria nivea L. Gaud), a cDNA library was constructed by suppression subtractive hybridization using anthracnose-resistant cultivar Huazhu no. 4. The cDNAs from Huazhu no. 4, which were infected with Colletotrichum gloeosporioides, were used as the tester and cDNAs from uninfected Huazhu no. 4 as the driver. Sequencing analysis and homology searching showed that these clones represented 132 single genes, which were assigned to functional categories, including 14 putative cellular functions, according to categories established for Arabidopsis. These 132 genes included 35 disease resistance and stress tolerance-related genes including putative heat-shock protein 90, metallothionein, PR-1.2 protein, catalase gene, WRKY family genes, and proteinase inhibitor-like protein. Partial disease-related genes were further analyzed by reverse transcription PCR and RNA gel blot. These expressed sequence tags are the first anthracnose resistance-related expressed sequence tags reported in ramie.

Introduction

Ramie (Boehmeria nivea L. Gaud), a perennial herbaceous plant grown for its fibers and a member of the Urticaceae family, is planted principally in China and other Asian countries including the Philippines, India, South Korea, and Thailand (Wang et al. 2008). Ramie fibers are long and strong and possess good durability and absorbency with excellent luster (Zheng et al. 2000). In China, about 130,000 ha are planted with ramie. Ramie crops are usually infected by diseases caused by fungi, viruses, or nematodes, which decrease the yield as well as the quality of the fibers. Anthracnose is a major disease of ramie caused by the fungus Colletotrichum, first described in 1831 by Corda. Its classification is phylum Deuteromycotina, class Coelomycetes, order Melanconiales, family Melanconiaceae. This important plant fungal disease is widespread in temperate zones, particularly in subtropical and tropical regions. It often damages fruit trees, vegetables, flowers, medicinal plants, and the fruit, stems, and leaves of field crops (Sutton 1992; Agostini et al. 1992), causing large losses in agricultural production. In China, Colletotrichum boehmeriae Sawada has been reported as a pathogen of ramie, but with no further investigations (Sawada 1914, 1919). In our former research, we have identified the disease caused by Colletotrichum gloeosporioides (Wang et al. 2010a, b) which can infect the stalk and foliage, thereby degrading not only the quantity but also the quality of the fiber. Generally, the disease causes yield losses of 20–40% and as high as 55% (Li and Ma 1993).

To study resistance genes in ramie, we constructed a suppression subtractive hybridization (SSH) (Diatchenko et al. 1996) library induced by C. gloeosporioides. The SSH technique enables specific cloning of expressed sequence tags (ESTs) representing genes that are differentially expressed in different mRNA populations and isolates genes without prior knowledge of their sequence or identity. It uses common molecular biological techniques that do not require specialized equipment or analyses (Moody 2001). This method was very useful in screening new genes and construction of cDNA library (Lin et al. 2010; Wang et al. 2011; Dipnarayan et al. 2011; Zhang et al. 2011a, b; Zhou et al. 2011; Prabu et al. 2011; Yang et al. 2011). Given the current understanding of ramie at the molecular level, this technique is a suitable method for the detection of low-abundance, differentially expressed transcripts that can lead to the isolation of putative disease resistance and defense-related genes.

Materials and Methods

Material Preparation and Isolation of Colletotrichum

Huazhu no. 4, an elite ramie cultivar in Hubei Province, China that shows high resistance to ramie anthracnose (C. gloeosporioides), was planted in the Ramie Germplasm Resources Garden of Huazhong Agricultural University. The plants used in this study were transplanted into the greenhouse under daylight with additional illumination provided by high-pressure sodium lamps to give a photoperiod of 12 h.

C. gloeosporioides samples were collected from anthracnose lesions on ramie, the Bast Fiber Crops Institute in Changsha, Hunan Province, China. After surface sterilization of lesions with 0.1% mercuric chloride for 40 s, followed by autoclaved water wash three times, small blocks (9 mm2) of diseased bark were aseptically transferred to potato dextrose agar (PDA) plates. Cultures were dark-cultured at 25°C incubator. C. gloeosporioides isolates were purified by single spore culturing (Choi et al. 1999). Spore masses were picked up with a sterilized wire loop and streaked on to the surface of water agar plates which were dark-cultured at 25°C incubator overnight. A single germinated spore was picked up with a sterilized needle and then transferred onto PDA plates. Pure cultures were stored at 4°C on PDA slants. Isolates were deposited in Huazhong Agricultural University, Hubei Province, China.

The ramie Huazhu no. 4 was sprayed with 20 ml of a solution of anthracnose spores (1 × 106/ml). The infected ramie plants and the control samples were cultured under conditions of 25°C and 80% humidity. Leaves were taken 6, 12, 24, 36, 48, and 72 h after inoculation. The samples were then quick-frozen in liquid nitrogen and stored at −80°C for RNA extraction.

Total RNA and mRNA Isolation

Total RNA was extracted from the samples by an optimized method (Wang et al. 2010a, b). The extracted total RNA was checked by electrophoresis on 1% agarose gels. mRNA was isolated using a PolyATract mRNA Isolation System kit (Promega, USA) according to the manufacturer’s instructions.

SSH Library Construction

An SSH library was constructed using a Clontech PCR-SelectTM cDNA Subtraction Kit (BD Biosciences, USA) based on the manufacturers’ instructions. The final PCR products were purified using an SBS Quick PCR Purification kit (SBS Genetech, Wuhan, China). The subtracted cDNA was inserted directly into pGEM®-T Easy Vector (Promega, USA) and then transformed into JM-109 Escherichia coli cells (Promega, USA) plated onto lysogeny broth (LB) agar containing 100 mg L−1 ampicillin, 1 mM isopropyl β-d-thiogalactopyranoside, and 80 mg L−1 5-bromo-4-chloro-3-indolyl β-d-galactopyranoside. Following incubation at 37°C overnight, positive transformants based on blue/white color selection were picked and arrayed into 96-well microplates and then cultured in LB medium containing 100 mg L−1 ampicillin. The resultant subtractive cDNA library was stored at −80°C with 15% glycerol.

Amplification of cDNA Inserts

The cDNA inserts were amplified by PCR in a 96-well plate (PTC-100, Germany) with nested PCR primer 1 and primer 2R, which were included in the PCR-SelectTM cDNA Subtraction Kit. PCR reactions (25 μL) contained 17.3 μL distilled water, 1 μL bacterial cultures, 1 μL each of nested PCR primer 1 and primer 2R (10 μM each), 2.5 μL 10× ExTaq buffer, 1.5 μL Mg2+ buffer (25 mM), 0.5 μL dNTPs (10 mM each), and 1 U ExTaq polymerase (SBS Genetech, Wuhan, China). Reaction samples were first denatured at 94°C for 5 min, followed by 30 cycles at 94°C for 20 s, 55°C for 45 s, and 72°C for 45 s, with a final extension at 72°C for 10 min. All PCR products were analyzed by agarose gel electrophoresis and clones of segments <200 bp or two or more segments eliminated.

Preparation of Probes and Differential Screening of the Subtracted Libraries

To exclude false-positive clones and to provide further data on the relative expression levels of the cloned cDNAs, initial screening of the library to remove false positives was performed by reverse northern analysis with total labeled cDNA from unsubtracted tester cDNA and unsubtracted driver cDNA as probes. Total RNAs were extracted from healthy ramie leaves and inoculated leaves. Then, they were reversed into cDNA separately. And then, they were labeled as tester cDNA and driver cDNA probes.

Five microliters of PCR products (about 100 ng) was mixed with 5 μL fresh 0.6 M NaOH for denaturation and 1.2 μL of the mixture was printed onto two Hybond-N+ nylon membranes (Boehringer, Mannheim, Germany). The membranes were neutralized in 0.5 M Tris–HCl (pH 7.5) for 5 min and rinsed in distilled water for 30 s. Samples were cross-linked to the membranes by baking for 2 h at 80°C and then stored at −4°C until use. Unsubtracted driver cDNA and tester cDNA probes were labeled using a DIG High Prime DNA Labeling and Detection Starter Kit II (Roche Diagnostics GmbH, Germany). Preparation of probes and hybridization were performed exactly according to the manufacturer’s instructions. The results of two hybridizations were recorded for each clone, and those showing the most marked differential expression were selected for sequencing (Sunny Biotechnology, Shanghai, China).

EST Sequencing and Analysis

Homology searches were conducted for all sequences using the GenBank database and the BLAST algorithms at the National Center for Biotechnology Information (NCBI) network service (http://www.ncbi.nlm.nih.gov/BLAST/). The cDNAs were classified according to the E-values generated in the BLAST searches. E-values <1e−10 were deemed to indicate significant homology and functional assignment. ESTs with E-values >1e−10 were deemed to have no significant homology to any known protein and assumed to be novel. EST sequences were assigned manually to functional categories based on a previous catalog system (Bevan et al. 1998).

RT-PCR Analysis

Total RNA was extracted from samples taken 6, 12, 24, 36, 48 and 72 h after inoculation. First-strand cDNAs were generated from 50 to 100 ng RNA samples using a reverse transcription PCR (RT-PCR) kit (Toyobo, Japan). To determine the expression of candidate genes, PCR was performed with one tenth of the first-strand cDNA template and the gene-specific primer pairs. Gene-specific RT-PCR primers were designed with Primer Premier 5.0 according to the cDNA sequences and were synthesized commercially (SBS Genetech, Wuhan, China). The following primer pairs were used: for biquitin-conjugating enzyme, fwd (5′-AAGCGACGATAATGGGTC-3′) and rev (5′-GCTATTGATGTTCGGGTGA-3′); for polypeptide for cytokinin-repressed mRNA, fwd (5′-AACACTATCGTCGGAGG-3′) and rev (5′-CAAGAAGGCACAGAGCAG-3′); for putative heat-shock protein 90, fwd (5′-AAGACCCTTAAACAACAAGA-3′) and rev (5′-AGTCGAGCCTAAACATCAG-3′); for actin factor, fwd (5′-GCCTGCTGCTTCCATTCC-3′) and rev (5′-TGGCTTACATTGCCCTTGA-3′); for pathogenesis-related protein, fwd (5′-CGGGCAACTACAATGGAGAA-3′) and rev (5′-AGCCATTTCAGAATCAAC-3′); for 3-hydroxyacyl-CoA dehydrogenase, fwd (5′-ATCCGTAACAGACCAAGA-3′) and rev (5′-GAGACAAACATCGCAAGT-3′); for Glycine max nodulin-26 mRNA, fwd (5′-GAGAAATAGCAAACCAACCC-3′) and rev (5′-CCCACCACGGACTACTGA-3′); for pathogenesis-related protein PR10, fwd (5′-GAAGGGCGGGCATGAGAT-3′) and rev (5′-TGTCGCTTTGTTTAGTTGTAGG-3′); and for stress and pathogenesis-related protein, fwd (5′-GGCATGAGATCAAGGAGG-3′) and rev (5′-TGTCGCTTTGTTTAGTTGTAGG-3′). The programs differed because of the different primers. The PCR products (8 μL) for each sample were electrophoresed in a 1.5% ethidium bromide agarose gel and viewed under ultraviolet light.

Northern Blot Analysis

Total RNA (20 μg) from samples taken at different times after inoculation were separated on a 1.0% agarose/formaldehyde gel and transferred to Hybond-N+ nylon membranes, by downward capillary transfer with 20× SSC. Specific probes were generated by purified PCR amplification of the relevant differentially expressed clones and labeled using a DIG High Prime DNA Labeling and Detection Starter Kit II (Roche Diagnostics GmbH, Germany). Labeling of probes and hybridization were performed exactly according to the manufacturer’s instructions.

Results and Discussion

A subtracted cDNA library was constructed using the healthy ramie Huazhu no. 4 cultivar and Huazhu no. 4 inoculated with C. gloeosporioides. The library comprised approximately 1,000 clones, with insertions ranging from 50 to 1,100 bp. One hundred seventy-two clones of greater signal intensity were found to be regulated and were sequenced. The EST sequences obtained were identified by homology searches in BLASTN and BLASTX databases. The partial results are summarized in Table 1. One hundred thirty-two single ESTs were assigned to functional categories and classified into 14 putative cellular functions, basing on the functional categories established for Arabidopsis (Bevan et al. 1998). The distribution of the genes is illustrated in Fig. 1. In the 132 single ESTs, 70 had high similarity to database entries. The largest set of genes (28%) was assigned to disease/defense. Unclear classification and unclassified were the second (17%) and the third (11%) largest groups, respectively. Based on the results of northern hybridization and EST sequencing, RT-PCR was performed to evaluate the differential expression of selected disease-related genes (Fig. 2). As indicated by the SSH results (Table 1), we selected nine clones: ubiquitin-conjugating enzyme 2 mRNA, polypeptide for cytokinin-repressed, putative heat-shock protein 90 (Hsp90), actin factor, pathogenesis-related protein, 3-hydroxyacyl-CoA dehydrogenase, nodulin-26, pathogenesis-related protein PR10, and pathogenesis-related protein, all of which exhibited high homology in the BLASTX comparison with the NCBI database. RT-PCR showed high levels of the genes encoding stress and pathogenesis-related protein at the time of induction of pathology and in the control group (CK) (Fig. 2). The clone related to actin factor showed peak expression after 24 h. Expression of some genes was induced after the induction of pathology: polypeptide for cytokinin-repressed mRNA showed differences after 12 h and peaked at 24 h, after which expression was maintained. However, another clone related to the stress and pathogenesis-related protein was expressed only at 12–48 h. We conclude that expression of this gene was brought about the induction of pathology. The clone with high homology to putative heat-shock protein 90 was highly expressed in CK but showed a downward trend after the induction of pathology. The Hsp90 proteins might play a role either in stabilization of proteins involved in disease resistance or in assembly of multi-subunit complexes required for activation of R proteins and intracellular disease resistance signaling pathways (Scofield et al. 2005). In the RNA gel blot analysis, we used probes synthesized from clone 645 (polypeptide for cytokinin-repressed) and clone 696 (pathogenesis-related protein PR10). From Fig. 3, it can be seen that it is likely that the corresponding genes were expressed at similar levels as observed by RT-PCR.

Table 1

Representative clones of SSH library from Boehmeria infected by C. gloeosporioides

Clone no.	Number of identical clonesa	Best E-valueb	Identity percentage (%)	Best match in database
Pathogenesis-related proteins/plant defense
15	emb|CT029927.1|	8e−10	92	Poplar cDNA sequences (SSH library infected with rust fungus)
39	gb|AAT02522.1|	0.006	62	Metallothionein 1a
42	gb|DQ907240.1|	0.059	90	Heat-shock transcription factor 1
46	emb|CT029469.1|	8e−10	92	Poplar cDNA sequences (SSH library infected with rust fungus)
72	ref|XM_002519612.1|	8e−95	83	(S)-2-hydroxy-acid oxidase
73	emb|CT029702.1|	0.063	100	Poplar cDNA sequences (SSH library infected with rust fungus)
88	gb|EU311423.1|	3e−21	77	Heat-shock protein 90 mRNA
141-1	ref|XM_002273749.1|	2e−09	75	Thaumatin-like protein
185	gb|AY919870.1|	4e−21	70	Betula platyphylla metallothionein-like protein mRNA
196	emb|CU231781.1|	6e−63	80	Populus EST from mildly drought-stressed leaves
205	gb|DQ186610.1|	3e−53	90	Catharanthus roseus calmodulin 2 mRNA
223	emb|AJ007349.1|TAE7349	0.002	61	Triticum aestivum mRNA for PR-1.2 protein
235	emb|CT028832.1|	0.019	100	Poplar cDNA sequences
248	emb|CAD42908.1|	3e−27	78	Catalase
367	gb|EU586200.1|	4e−135	86	Biquitin-conjugating enzyme 2 mRNA
372	gb|AAL36933.1|AF213481_1	2e−10	100	Gamma-tocopherol methyltransferase
425	emb|CT028832.1|	0.02	100	Poplar cDNA sequences (SSH library infected with rust fungus)
426	gb|AY847454.1|	1e−07	74	Metallothionein class II
459	gb|EU019559.1|	3e−52	83	Glycine max WRKY14
495	emb|CU231781.1|	1e−54	76	Populus EST from mildly drought-stressed leaves
578	gb|EU024476.1|	1e−07	97	Arachis hypogaea oxalate oxidase mRNA
585	gb|DQ497797.1|	8e−40	74	Lipoxygenase 2
589	emb|AJ007349.1|	4e−04	86	Triticum aestivum mRNA for PR-1.2 protein
599	dbj|AB091685.1|	3e−04	100	Daidzein 7-0-methyltransferase
661-2	emb|AM748435.1|	5e−07	84	Vigna unguiculata partial mRNA for putative plant disease resistance response protein family
640	ref|XM_002519612.1|	7e−95	84	Ricinus communis (S)-2-hydroxy-acid oxidase
667	gb|AY594295.1|	8e−17	68	Populus trichocarpa × Populus deltoides metallothionein 1a (MT1a)
696	gb|EU117126.1|	4e−24	81	Prunus domestica cultivar Violette pathogenesis-related protein PR10
702	emb|CU232440.1|	1e−07	91	Populus EST from mildly drought-stressed leaves
704	emb|AM748435.1|	5e−12	91	Putative plant disease resistance response protein family
719	emb|AJ130889.1|FSY130889	6e−14	81	Fagus sylvatica mRNA for stress and pathogenesis-related protein
761	gb|EU526202.1|	8e−13	73	PR10 protein
766	gb|AY158153.1|	4e−09	92	Citrus × paradisi type I proteinase inhibitor
792	emb|AJ496418.1|PPE496418	7e−42	79	Prunus persica mRNA for catalase
933	emb|CU226741.1|	2e−89	83	Populus EST from severely drought-stressed leaves
Metabolism
96	gb|EF568776.1|	7e−109	81	Prunus persica beta-galactosidase protein
288	gb|EU106889.1|	1e−14	79	Jatropha curcas chloroplast omega-6 fatty acid desaturase mRNA
600	ref|XM_002514233.1|	2e−51	77	Ricinus communis 3-hydroxyacyl-CoA dehydrogenase
664	dbj|AB221013.2|	5e−42	83	Beta vulgaris BvSDC1 mRNA for putative serine decarboxylase
707	ref|XM_002283715.1|	2e−25	85	Vitis vinifera malate dehydrogenase (NADP+)
740	ref|XM_002515921.1|	2e−53	85	Ricinus communis malic enzyme
Cell growth/division
178	gb|FJ226587.1|	4e−76	80	Glycine max late-embryogenesis abundant protein 2
233	gb|FJ156098.1|	4e−43	80	Translationally controlled tumor protein
587	ref|XM_002516884.1|	0.012	87	Translationally controlled tumor protein
645	dbj|D29684.1|CUSCR9	4e−21	83	Polypeptide for cytokinin-repressed mRNA
663	emb|AJ252088.1|	7e−05	77	mRNA for gibberellin 20-oxidase
757	emb|AJ238624.1|	2e−89	88	Translation initiation factor
Transporters
24	ref|XM_002336761.1|	2e−38	76	Populus trichocarpa oligopeptide transporter OPT family
153	emb|AJ534339.1|	0.008	86	Agaricus bisporus partial mRNA for putative inorganic phosphate transporter (ipt gene)
189	ref|XM_002336761.1|	8e−36	75	Populus trichocarpa oligopeptide transporter OPT family
457	emb|AJ534339.1|	3e−04	88	mRNA for putative inorganic phosphate transporter
660-2	emb|AJ272202.1|	2e−06	90	Arabidopsis thaliana mRNA for mitochondrial half-ABC transporter
665	gb|EF453694.1|	2e−70	78	Populus trichocarpa Mn-specific cation diffusion facilitator transporter (MTP11.2) mRNA
787	gb|DQ372573.1|	1e−05	84	Gossypium hirsutum transferring glycosyl protein (gt) mRNA
Transcription
218	emb|AM293617.1|	0.075	100	Platanus × acerifolia partial mRNA for putative transcription factor
227	emb|AM397241.1|	0.006	93	Ty1-copia reverse transcriptase (RT gene) and partial DfRedu pseudogene
230	dbj|AB235910.1|	0.20	91	Triticum aestivum gene for RNA ligase isoform 3
650	emb|AM293617.1|	0.002	100	Platanus × acerifolia partial mRNA for putative transcription factor
651	gb|AF190657.1|AF190657	3e−39	83	Nicotiana tabacum clone 7 poly(A)-binding protein (PABP) mRNA
926	gb|AY224135.1|	7e−11	77	Pelargonium × hortorum 4 coumarate CoA ligase mRNA
929	gb|AY224135.1|	5e−12	75	Pelargonium × hortorum 4 coumarate CoA ligase mRNA
936-1	gb|AY224135.1|	2e−16	81	Pelargonium × hortorum 4 coumarate CoA ligase mRNA
943-1	gb|AY224135.1|	2e−17	82	Pelargonium × hortorum 4 coumarate CoA ligase mRNA
971-1	gb|AY224135.1|	1e−06	71	Pelargonium × hortorum 4 coumarate CoA ligase mRNA
Signal transduction
205	gb|AF292108.1|AF292108	3e−97	90	Prunus avium calmodulin mRNA
367	gb|EU586200.1|	4e−135	86	Ubiquitin-conjugating enzyme 2
449	ref|XM_002275581.1|	3e−80	79	Vitis vinifera ubiquitin-activating enzyme
584	ref|XM_002512452.1|	2e−73	83	Ricinus communis GTP-binding protein sar1
Transposons
495-2	gb|EU558534.1|	8e−40	76	Arabidopsis lyrata clone SINE9 transposon-insertion display band
936-2	gb|EU558528.1|	1e−34	81	Arabidopsis lyrata clone CACTA2 transposon-insertion display band genomic sequence
944	gb|EU558534.1|	4e−25	72	Arabidopsis lyrata clone SINE9 transposon-insertion display band genomic sequence
Energy
118	ref|XM_002526592.1|	8e−23	70	Ricinus communis photosystem I reaction center subunit N
195	emb|AM748509.1|	0.003	100	Vigna unguiculata partial mRNA for putative single-stranded nucleic acid binding R3H
Protein synthesis
119	ref|XM_002519706.1|	8e−52	82	Ricinus communis 60S ribosomal protein L32
194	emb|AJ251333.1|	4e−71	81	Pisum sativum aair gene for acetohydroxy acid isomeroreductase
754	emb|AJ831469.1|	1e−09	76	Pisum sativum mRNA for putative glycine-rich protein precursor
756	emb|AJ697756.1|	0.026	88	Eucalyptus globulus subsp. globulus partial mRNA for putative ribosomal protein
Protein destination and storage
446	emb|AM748427.1|	0.005	91	Vigna unguiculata partial mRNA for putative cathepsin B-like cysteine protease
Intracellular traffic
609	gb|L12257.1|	6e−71	82	Glycine max nodulin-26 mRNA
654	gb|GQ354526.1|	0.002	91	Dekkera bruxellensis strain CBS 2499 mitochondrion
736	gb|L12257.1|	7e−71	82	Glycine max nodulin-26 mRNA
745	ref|XM_002522770.1|	2e−18	74	Ricinus communis mitochondrial ribosomal protein
Cell structure
204	gb|AY487315.1|	3e−97	89	Ficus carica actin 1 mRNA
943-2	emb|AJ420893.1|	4e−34	76	Entandrophragma cylindricum microsatellite DNA
Secondary metabolism
971-2	dbj|AB261987.1|	1e−04	79	Photinus pyralis mRNA for luciferase

The GenBank accession number of the clones in the Table 1: GW691509-GW691606; GO653137-GO653162

aBest match in the GenBank (gb), DNA Data Bank of Japan (dbj), European Molecular Biology Laboratory (emb), or NCBI database. The accession number of the match is in parentheses

bThe expected (E) value refers to the number of matches expected by chance alone. The lower the E value, the more strongly supported the match

Fig. 1

Functional categories of genes included in the anthracnose-resistant subtractive cDNA library, based on the classification described in Bevan et al. (1998)

Fig. 2

RT-PCR analysis of cDNAs representative of different types of gene expression patterns. CK, control leaves; 6h, 12h, 24h, 36h, 48h, 72h, 6, 12, 24, 36, and 72 h, respectively, after inoculation with C. gloeosporioides

Fig. 3

Northern hybridization of total RNA from anthracnose-resistant ramie cultivar Huazhu no. 4 using clones 645 and 696. A, B, ECL detection of the 645 and 696 probes, separately; CK, control leaves; 12h, 24h, 12 and 24 h after inoculation with C. gloeosporioides

Breeding for stable anthracnose resistance is difficult because of the variation among C. gloeosporioides for pathogenicity and virulence. In addition, the host–pathogen interaction is poorly understood. Few reports on anthracnose and the cloning and characterization of ramie resistance gene homologs have been published. In this study, the transcription profiles of uninfected and infected leaves of the anthracnose-resistant Huazhu no. 4 ramie cultivar were compared by SSH analysis. The SSH library was not very extensive, which might be due to the efficiency of the adapters ligation, which plays a key role in the subtractive library construction. Additionally, the efficiency depends upon several factors such as the quality of cDNA synthesis, the completeness of enzyme digestion, and the quality of purification. However, we obtained many disease resistance- and defense-related ESTs. There were only about 200 ESTs of ramie in the GenBank database, so our results are meaningful for future research into disease resistance in ramie. We found several transcripts encoding proteins belonging to groups including pathogenesis-related proteins, metallothioneins, severe drought-stressed proteins, Hsp90, thaumatin-like protein, catalase (CAT1), the WRKY family, and oxalate oxidase.

Although PCR-selective cDNA subtraction is a powerful tool for identifying differentially expressed genes, and subtractive products may contain cDNAs that are common to or have similar levels in the two tissues, false-positive clones may occur (Zeng et al. 2006). As a result, we performed RT-PCR for further identification. According to the RT-PCR results, the resistant ramie cultivar had many resistance genes before being infected, including ubiquitin-conjugating enzyme, PR-1.2 protein, nodulin-26 mRNA, and mRNAs for stress and pathogenesis-related protein. These proteins maintained their high levels before and after inoculation and played an important role in the resistance process. The pathogenesis-related proteins included several proteins that are involved in the detoxification of reactive oxygen species (Mishra et al. 2007). Another aspect of our findings was the decline of Hsp90, a molecular chaperone responsible for the folding and functions of important cellular proteins including steroid hormone receptors, protein kinases, and proteins controlling the cell cycle and apoptosis (Schulte et al. 1998). Wicklow et al. (2009) demonstrated that monorden from Colletotrichum graminicola is a potent selective inhibitor of Hsp90. Monorden might be a competitive ligand for the ATP-binding site of Hsp90, thus preventing the ATPase activity necessary for its role as a chaperone (Roe et al. 1999). It could be confirmed by further experiments on our pathogen, C. gloeosporioides. 3-Hydroxyacyl-CoA dehydrogenase (HADH) deficiency is an autosomal recessive metabolic disorder, resulting from mutations in the HADH gene in chromosome 4q22-q26 (http://www.ncbi.nlm.nih.gov/entrez/dispomim.cgi?id=231530; Vredendaal et al. 1998), which has so far only been described in a few patients (Martins et al. 2011). In plant of castor bean (Ricinus communis) endosperm and potato tubers, this enzyme had a high expression level (Courtois-verniquet and Douce 1993).

The conclusion of this study is that, in addition to the expression of several well-characterized pathogenesis-related proteins, high-level constitutive expression of genes encoding thaumatin-like protein, calmodulin, catalase, and proteinase inhibitor seems to be a major factor in rendering ramie resistant to infection by the fungus C. gloeosporioides. The recovery of many ESTs indicates that ramie quickly responds to C. gloeosporioides challenge and suggests that early transcriptional events can play an important role in determining whether ramie succumbs to or resists infection. This study is the first global analysis of genes in ramie with anthracnose. The ESTs identified following C. gloeosporioides challenge in this study should provide a useful genomic resource for biologists and plant breeders in developing new strategies for improving resistance to anthracnose in ramie.