key: cord-0000058-9r62ffew authors: Pan, Yi-Shin; Lee, Yun-Shien; Lee, Yung-Lin; Lee, Wei-Chen; Hsieh, Sen-Yung title: Differentially profiling the low-expression transcriptomes of human hepatoma using a novel SSH/microarray approach date: 2006-05-31 journal: BMC Genomics DOI: 10.1186/1471-2164-7-131 sha: 4db3931c4b24ae432d68b7e24f593b1272445523 doc_id: 58 cord_uid: 9r62ffew BACKGROUND: The main limitation in performing genome-wide gene-expression profiling is the assay of low-expression genes. Approaches with high throughput and high sensitivity for assaying low-expression transcripts are urgently needed for functional genomic studies. Combination of the suppressive subtractive hybridization (SSH) and cDNA microarray techniques using the subtracted cDNA clones as probes printed on chips has greatly improved the efficiency for fishing out the differentially expressed clones and has been used before. However, it remains tedious and inefficient sequencing works for identifying genes including the great number of redundancy in the subtracted amplicons, and sacrifices the original advantages of high sensitivity of SSH in profiling low-expression transcriptomes. RESULTS: We modified the previous combination of SSH and microarray methods by directly using the subtracted amplicons as targets to hybridize the pre-made cDNA microarrays (named as "SSH/microarray"). mRNA prepared from three pairs of hepatoma and non-hepatoma liver tissues was subjected to the SSH/microarray assays, as well as directly to regular cDNA microarray assays for comparison. As compared to the original SSH and microarray combination assays, the modified SSH/microarray assays allowed for much easier inspection of the subtraction efficiency and identification of genes in the subtracted amplicons without tedious and inefficient sequencing work. On the other hand, 5015 of the 9376 genes originally filtered out by the regular cDNA microarray assays because of low expression became analyzable by the SSH/microarray assays. Moreover, the SSH/microarray assays detected about ten times more (701 vs. 69) HCC differentially expressed genes (at least a two-fold difference and P < 0.01), particularly for those with rare transcripts, than did the regular cDNA microarray assays. The differential expression was validated in 9 randomly selected genes in 18 pairs of hepatoma/non-hepatoma liver tissues using quantitative RT-PCR. The SSH/microarray approaches resulted in identifying many differentially expressed genes implicated in the regulation of cell cycle, cell death, signal transduction and cell morphogenesis, suggesting the involvement of multi-biological processes in hepato-carcinogenesis. CONCLUSION: The modified SSH/microarray approach is a simple but high-sensitive and high-efficient tool for differentially profiling the low-expression transcriptomes. It is most adequate for applying to functional genomic studies. Microarray is a powerful technique for simultaneously determining the expression of thousands of genes [1] [2] [3] . Such studies can quickly yield a genome-wide description of mRNA expression, called transcriptomes, in a given cell or tissue at a given physiologic or pathologic condition [4] . Even though, one of the main challenges in such genome-wide gene-expression profiling is the difficult inspection of genes with rare transcripts. On the other hand, PCR-based suppressive subtractive hybridization (SSH) techniques are highly sensitive for identifying differences in gene expression [5] [6] [7] [8] [9] . However, the potentials of SSH in assaying dynamic changes of gene expression in minute levels have never been addressed before. In addition, SSH techniques are also restricted in terms of limited specificity and difficulties in identifying enriched genes. Combining the SSH technique with highthroughput screening of the harvested clones through the use of cDNA microarrays could greatly reduce the tedious work for northern blot analysis, as well as the likelihood of false-positive clones enriched via SSH [10] . Such combined approaches by printing the clones obtained from the SSH amplicones on chips have been successfully used for profiling the differentiation of gene expression [10] [11] [12] [13] . Nevertheless, in such approaches, the genes of the subtracted clones remain to be sequenced for identification, and a large portion of redundancy in the enriched amplicons must also be identified. Moreover, since the targets used to hybridize the amplicon clones printed on chips were the un-enriched cDNA pools rather than the subtracted, enriched clones, such approaches would not increase the sensitivity for detection of the low-expression genes. Herein, we report our modifications using the subtracted amplicons as the targets to hybridize the pre-prepared microarray chips for SSH/microarray analysis. Since all of the probes on the microarray chips have been well characterized, the hundreds and thousands of genes in the subtracted amplicons can be determined by a single hybridization. Moreover, the relative expression status between the compared tissues for each gene can be augmented and easily determined by combining the use of the targets prepared both from forward and reverse subtractions of SSH. We named this modified approach "SSH/microarray" and used human hepatocellular carcinoma as a model to demonstrate the feasibility of this approach. The SSH/cDNA microarray versus regular cDNA microarray approaches To conduct the comparative transcriptomic studies on human hepatoma particularly for the genes at low-expres-sion levels, we set out an approach as shown in Figure 1 . We used the PCR-based suppressive subtractive hybridization to enrich the differentially expressed cDNA clones between human hepatoma and non-hepatoma liver tissues. The SSH was conducted both forwards using RNA prepared from hepatoma and non-hepatoma liver tissues as tester and driver, respectively, to yield the hepatoma (T-N) subtracted amplicons, and reversely using the mRNA prepared from non-hepatoma and hepatoma liver tissues as tester and driver respectively to yield the non-hepatoma (N-T) subtracted amplicons as well. Instead of the previously reported combination with cDNA microarray approaches, in which clones in the subtracted cDNA libraries were cloned and then printed on chips, we directly labeled the subtracted T-N and N-T amplicons and then used as targets to hybridize the pre-made microarray chips printed with the known 14,811 cDNA clones. We named this modified approach as "SSH/microarray". The same tissue RNA pairs and cDNA microarray chips were also used for the regular cDNA microarray assays for comparison. Since human hepatoma generally occurs in patients of chronic hepatitis with/without cirrhosis, to identify the genes implicated in the common biological processes leading to hepatocellular carcinogenesis, we selected three hepatoma patients with different underlying liver disease. The first was with chronic hepatitis C and cirrhosis, the second had chronic hepatitis B and cirrhosis, and the third had chronic hepatitis B but no clinical or histological evidence of cirrhosis. Figure 2 demonstrates the representative results obtained from the corresponding regular cDNA microarray and the SSH/microarray assays. Of note, in the SSH/microarray assays the genes in the subtracted T-N and N-T amplicons were readily identified without tremendous sequencing works, which were required for the original SSH and microarray combination assays. In addition, many of the clones, including those of the house-keeping genes, such as β-actin and GAPDH, with unremarkable difference in hybridization intensities by the regular cDNA microarray assays were efficiently excluded from the subtracted amplicons ( Figure 2 , indicated by white circles). This indicated the efficiency of subtraction hybridization. Most importantly, many of the clones with hybridization intensity below the evaluation threshold by the regular cDNA microarray approaches (hybridization intensities lower than or close to the noise) became detectable by the SSH/ microarray method ( Figure 2 , indicated with red circles). These observations suggested a more sensitive detection of the low abundance transcripts using the SSH/microarray assays. Figure 1 Experimental flowchart. Messenger RNA was prepared from hepatoma (T) and the corresponding non-hepatoma liver tissues (N), and then subjected to 1) PCR-based suppressive subtractive hybridization (SSH) followed by using the resulted subtracted cDNA libraries as targets for cDNA microarray analysis (the SSH/microarray), and 2) conventional cDNA microarray analysis. SSH was performed in both the forward (T as tester) and reverse (N as tester) direction to enrich up-regulated (T-N amplicon) as well as down-regulated transcriptomes (N-T amplicon) in human hepatoma, respectively. The two subtracted amplicons were labeled with fluorescent cy-dyes as targets for microarray analysis, in which dye-swapping approaches were used. The results thus obtained were then compared to those obtained from the conventional cDNA microarray assays. The differentially expressed genes were categorized into three groups: I) only detected by conventional cDNA microarray approaches, II) identified by both conventional cDNA microarray and SSH/microarray approaches, and III) only obtained from SSH/microarray assays. The results were further confirmed by qRT-PCR in 6 genes randomly selected from group III differentially expressed genes. Labeled and used to hybridize chips To further address whether the SSH/microarray approach was able to assay the expression of those genes with rare transcripts, we selectively inspected those genes with low expression. Of the total of 14,811 clones on the chips, 9376 clones that were filtered out for further analyses in the regular cDNA microarray studies due to the low hybridization intensity (intensities <500 unit after the noise subtracted) ( Figure 3A ) were selected for examination of their hybridization intensity in the SSH/microarray assays ( Figure 3B ). Of these, 5015 clones were enriched via the SSH/microarray assays and could be analyzed for differential expression (intensities >500 unit after the noise intensity subtracted) ( Figure 3B &3C ). That is, the SSH/microarray approach allowed about 100% more of the genes, most of which were low abundant, for the subsequent analyses. Since SSH specifically amplified the difference of gene expression, the SSH/microarray approach should also be able to identify more differentially expressed genes between hepatoma and non-hepatoma liver tissues, which were originally undetectable by the regular cDNA microarray methods. As shown in Fig 3C, of the 5015 lowexpression genes identified only by the SSH/microarray assays, we identified additional 512 and 1923 genes with at least a 2-fold decrease and increase in hepatoma tissues, respectively. We further examined the sensitivity for detection of differential expression by both approaches. We compared the ratio of Cy5/Cy3 intensity of each clone obtained from the regular cDNA microarray assays to that obtained from the SSH/microarray assays. As shown in Figure 4 , a total of 3028 genes with a ratio lower than 2 folds (-1 < log 2 < 1) using the cDNA microarray assays had a ratio greater than 2 folds (log 2 > 1 or < -1) using the SSH/microarray assays The representative results of the regular cDNA microarray vs Figure 2 The representative results of the regular cDNA microarray vs. SSH/microarray assays. RNA were prepared from the hepatoma and para-hepatoma liver tissues of a patient with HCC. Left panels are the representative results obtained via the regular cDNA microarray assays using the Cy 5 and Cy3 labeled aRNA generated from hepatoma and non-hepatoma liver tissues, respectively. Right panels are the representative results obtained via the SSH/microarray assays using Cy3 and Cy5 to label the forward and reverse subtracted amplicons, respectively. Lower panels are the close views of the upper panels. White circles indicate the clones, which had equal Cy3 and Cy5 hybridization intensities in the conventional cDNA microarray assays but been subtracted away from the subtracted amplicons. Red circles mark the clones, which had low Cy5 and Cy3 hybridization intensities in the conventional cDNA microarray but presented with differential expression in the subtracted amplicons. SSH/microarray ( Figure 4 , 1091 and 1937 genes indicated by the red box, respectively). That is, using the modified SSH/microarray approaches resulted in much more differentially expressed genes, particularly for those with rare transcripts that using the regular cDNA microarray approaches. A total of 69 and 701 genes differentially expressed in human hepatoma were finally identified (based on the data derived from the three tissue pairs in duplicate with log 2 > 1 or < -1, and P < 0.01) by the regular cDNA microarray assays and by the SSH/microarray assays, respec-tively ( Figure 5 & see additional file 1). That is, the SSH/ microarray assays detected 10-fold more differentially expressed genes, particularly for those low-expression genes, than did the conventional cDNA microarray. To confirm the differential expression of genes with low expression identified by the SSH/microarray assays (including 446 and 255 up-and down-regulated genes, respectively), we quantified the relative expression level of nine randomly selected genes, whose differential expression was detected by the SSH/microarray assays, in 18 pairs of HCC and the matched non-HCC liver tissues. We found the results were consistent with those obtained from SSH/microarray assays ( Figure 6 ). Detection of differential low-expression transcriptomes by the SSH/microarray Figure 3 Detection of differential low-expression transcriptomes by the SSH/microarray. A) Of the 14811 distinct genes on the chips, 9376 genes had hybridization intensities lower than 500 units (the defined low intensity threshold for exclusion from further microarray analysis) in a representative set of results of the conventional cDNA microarray assay. B)These low-expression genes were then subjected to re-calibration for their Cy3/Cy5 intensities in the SSH/microarray assay and 5015 genes became detectable in SSH assays (blue spots). C) The distribution of these 5015 low-expression genes identified only in the SSH/microarray assays, including 512 and 1923 genes down-and up-regulated in hepatoma tissues, respectively. The data presented in this figure were based on the results derived from patient 1, a case of chronic hepatitis C with cirrhosis and hepatoma. We compared the HCC up-and down-regulated transcriptomes identified by the SSH/microarray assays in accordance with their potential molecular functions, implicated biological processes and sub-cellular localization. As presented in Table 1 Microarray techniques are limited in the detection of genes with low expression [14] , while subtractive hybridization methods are restricted by their specificity and the tremendous work needed for validating the results, as well as for sequencing to identify genes. There have been many reports using microarray techniques for rapid and high throughput validation of the subtraction specificity of SSH by spotting the enriched clones on the chips or membranes [10] [11] [12] [13] [15] [16] [17] [18] [19] [20] [21] [22] [23] . However, such approaches would not only require tremendously sequencing works to identify genes including a great number of redundant clones in the subtracted amplicons, but also lose the sensitivity of SSH techniques to detect the low abundance transcripts, since the sensitivity of microarray analysis is determined by the targets used to hybridize chips but not by the probes printed on chips. Herein we report our modified approaches. Instead of using the subtracted clones as probes spotted on chips, we directly labeled the enriched amplicons and used them as targets to hybridize the pre-made microarray chips for microarray analysis (named as "SSH/microarray" in this report). This approach allowed us not only to readily identify genes without tremendous sequencing works in the subtracted amplicons regardless of many redundant clones, but also to easily evaluate the subtraction efficiency and specificity. Moreover, the SSH/microarray approach made it possible to conduct a transcriptomewide identification of differentially expressed genes particularly for those with low expression. It has been reported that the absolute expression level is not a crucial determinant for identifying genes, while the relative difference in expression levels does impact on whether or not a gene is recovered by subtractive hybridization [5, 24] . In this report using the SSH/microarray approach, we identified about ten times more of the differentially expressed genes, particularly for those with low expression, in human hepatoma. Our findings successfully demonstrated the transcriptome-wide assays of differentially expressed genes at low abundance. This simple but very powerful approach would greatly facilitate future researches on functional genomics [15, 19, 25] . One of the main concerns about SSH techniques is their specificity. However, as combined with microarray tech- Figure 4 Comparison of the relatively expression by regular cDNA microarray vs. the SSH/microarray. The distribution of the ratios of Cy5/Cy3 intensities obtained by the conventional cDNA microarray assays vs. the SSH/microarray assays is shown. Only the clones with significant Cy3 and Cy5 intensities in both assays were included in the analysis. The clones inside the red-boxes were those with their log 2 ratios < 1 and > -1 via the conventional cDNA microarray assays, while their log 2 ratios > 1 or < -1 by the SSH/microarray assays. The data presented were based on those derived from patient 1, a case of chronic hepatitis C with cirrhosis and hepatoma. Log 2 ( Cy5/Cy3) of SSH Log 2 ( Cy5/Cy3) of cDNA niques, the subtraction efficiency can be readily evaluated. As shown in Figure 2 , the subtraction efficiency was determined by the rare presence of clones with equal hybridization intensities between the forward and reverse amplicons, and by the consistent exclusion of the housekeeping genes in the subtracted amplicons. The specificity was further confirmed by quantification of the differential gene expression between hepatoma and para-hepatoma The differentially expressed genes Figure 5 The differentially expressed genes. The genes differentially expression in human hepatoma were identified by either only the conventional cDNA microarray assays (26 genes) or SSH/microarray assays (658 genes), or both (43 genes). The results were derived from three pairs of hepatoma liver tissues with at least two-fold difference in gene expression and P < 0.01. APEG1, TTN, HTN1, MEST, EVL3, VLDLR, MYOG, TNNI2, HOXC11, IGF2, AHSG, TTID, BIRC1, PAFAH1B1, SYK, LMNA, TNFRSF11B, FALZ, BMP1, ETS2, NRD1, CXCL1, JAG1, SEMA4G, HEY1, ZNF22, CRMP1, TEAD4, IGFBP1, IGFBP3, IGFBP7, PTGS1 , NEUGRIN, CUGBP1) related to the regulation of cell differentiation and embryonic development, and 59 genes (ASGR2, CAP2, RGS5, ITGB4BP, NCK1, EVL3, PLCG2, EPHA7, IGF2, IFNGR2, AHSG, CAP1, PTPRF, CXCL2, GNAI1, TGFBR1, NCSTN, DOK1, FRBB3, SYK, MC1R, CD79B, CSNKIE, CD69, AVPR1A, GNB2L1, ACVR1, BIRC2, FLT4, CXCL1, FAG1, PNOC, SLC9A3R1, LANCL1, GNB3, ADORA2B, FGR, NCK1, RAB11A, PLCG2, LOC91614, PRKAR1A, FLJ22595, ARF3, TYK2, MAP3K7IP1, MC1R, TNFSF10, ARF4, VAV2, AVPR1A, GNB2L1, MAPK10, STMN1, SNX17, ADORA2B, ECT2, RGS5, IGFBP3 ) associated with signal transduction in response to extra-cellular proliferation and growth stimuli were found. Our findings that differentially expressed genes were related to multi-biological processes suggest the complexity of the molecular mechanisms for hepatocellular carcinogenesis. In this study, we modified the previously method of the combination of the suppressive subtractive hybridization and microarray techniques to differentially profile the low-expression transcriptomes of human hepatocellular carcinoma by directly labeling both of the reciprocal subtracted amplicons as target for cDNA microarray assays. Compared to SSH or previous SSH in conjunction with microarray approaches, this modified approach provided us with three additional advantages: 1) easy inspection of the subtraction efficiency, 2) avoidance of tremendous sequencing work for gene identification, 3) high sensitiv-ity for identifying the low-expression, differentially expressed genes. This approach allowed for the detection of about ten times more of the differentially expressed genes than did the regular cDNA microarray approach in human hepatoma, particularly for those with low expression. Using this approach, we identified many genes potentially implicated in human hepatocarcinogenesis, which were not identified before. For its high efficiency and high sensitivity, this SSH/microarray approach is powerful for the rapid differentially profiling the lowexpression transcriptomes, and most adequate for applying to functional genomic studies. For SSH and cDNA microarray analysis, hepatoma and the corresponding non-cancerous liver tissues were obtained from 3 patients who had liver surgery at the Chang Gung Memorial Hospital. For reverse-transcription real-time PCR, hepatoma and the matched non-hepatoma liver tissues were obtained from additional 18 patients of hepatoma. Diagnoses of HCC and non-hepatoma liver tissues were based on histo-pathologic findings. The Internal Review Board for Medical Ethics of Chang Gung Memorial Hospital approved the specimen collection procedures and informed consent was obtained from each subject or subject's family. Total RNA from hepatoma and the para-hepatoma liver tissues and the poly(A) RNA was prepared as described before [26, 27] . SSH was performed with the Clontech PCR-Select cDNA Subtraction Kit (Clontech Laboratories Inc., Palo Alto, CA) as described by the manufacturer but with the following modifications. Starting material consisted of 2 μg hepatoma mRNA as tester and 2 μg nonhepatoma liver tissue mRNA as driver, and vice versa. Primary and secondary PCR conditions were altered to increase specificity of amplification according to either plan A or B. Both A and B reduced the extension time and the number of cycles of the primary PCR to 2 min and 26 cycles, respectively. The primary PCR products were diluted 1/50 prior to use in the secondary PCR. All other aspects of plan A were as per the instructions of the manufacturer. Plan B diverged from plan A in two ways. First, the initial cycle of primary PCR was performed using annealing and extension times that had been reduced to 15 s and 1.5 min, respectively. Second, for subsequent cycles, the denaturing time was increased to 10 s while the annealing and extension times were reduced to 15 s and 1.5 min, respectively. In this study, we used the GMRCL Human 15 K set, Version 2 chips as previously described [28] , which contained 14,811 sequence-verified, human cDNA clones mapped to 12,530 distinct genes. All of the samples for the regular cDNA microarray or SSH/microarray assays were performed with the dye-swapping microarray design for minimizing labeling bias and statistical variances of data. For the regular cDNA microarray experiment, we used 2 μg of the total RNA for labeling and hybridization using a 3DNA Array 350RP Detection kit (Genisphere, PA, USA). For the SSH experiment, 1 μl subtractive PCR products were labeled with Cy3 and Cy 5-dCTP (NEN, Boston, MA) using random primers. Unincorporated fluorescent nucle-otides were removed using a Qiaquick PCR purification kit (Qiagen). The fluorescent-labeled DNAs were mixed with 30 μg of human cot-1 DNA (Invitrogen) and 100 μg yeast tRNA, precipitated and then resuspended in 30 μl of Microarray Hybridization Buffer Version 2 (Amersham Pharmacia). The hybridization solution was heated to 80°C for 10 min to denature the DNA and was then incubated for 30 min at 37°C, allowing cot-1DNA and yeast tRNA to block the repetitive sequences in genome probes. The probes were hybridized to a human cDNA microarray (GMRCL Human 15 K). We scanned the slides with a confocal scanner ChipReader (Virtek, Canada) and acquired Quantification of relatively gene expression in human hepatoma Figure 6 Quantification of relatively gene expression in human hepatoma. A total of nine differentially expressed genes including six (GBA, PTPRF, GNL3, ERBB3, OGDHL, FMO3 belonging to Class III) identified only by the SSH/microarray assays and three (C9, MT1F, MT1X belonging to class II) identified by both of the SSH/microarray and regular cDNA microarray assays were assayed for their relative expression between HCC and the corresponding non-HCC liver tissues of 18 patients of hepatoma using real-time semi-quantitative RT-PCR. The results are presented with the log 2 values of the gene expression ratio's between HCC vs. the matched non-HCC liver tissues. Left panel is the relatively expression level of these genes initially determined by the SSH/microarray and regular cDNA microarray assays, respectively. Of note, the difference of the expression for the six genes belonging to Class III was originally only identified by the SSH/microarray assays, but not by the regular cDNA microarray assays. the spot and background intensities with the GenePix Pro 4.1 software (Axon Instruments, Inc., CA, USA). The within-slide normalization was done using programs written with MATLAB 6.5 software (The MathWorks, Inc., MA, USA). To validate the results obtained from the SSH/microarray assays, we randomly selected six differentially expressed genes identified only by the SSH/microarray assays and three differentially expressed genes identified by both the regular cDNA microarray and SSH/microarray assays for the comparison of gene expression between hepatoma and the corresponding non-hepatoma liver tissues in eighteen patients of HCC using reverse transcription realtime PCR (qRT-PCR). Total RNA was extracted from tissues with Trizol reagents and reverse transcribed using the SuperScript III first strand synthesis system (Invitrogen, Carlsbad, CA). qRT-PCR was conducted using the ABI PRISM 7000 sequence detection system (Applied Biosystems). Pre-designed Assays on Demand TaqMan probes and primer pairs for these 9 genes were obtained from Applied Biosystems Incoporated (ABI) (Foster City, CA). For each gene, two to four sets of Taq-Man probes and primers were tested. The probes contained a 6-carboxyfluorescein phosphoramidite (FAM dye) label at the 5' end of the gene and a minor groove binder and non-fluorescent quencher at the 3' end. These were designed to hybridize across exon junctions. As a result, no fluorescent signal was generated by these assays when genomic DNA was used as a substrate, which confirmed that the assays measured only mRNA. Equal amounts of RNA were used for all qRT-PCR reactions, which were performed in triplicate, and 18S ribosomal RNAs were used as internal controls. All data (derived from three pairs of HCC and non-HCC liver tissues in a dye-swapping approach) were filtered so that 446 up-regulated and 255 down-regulated genes in HCC with at least two-fold difference and p value less than 0.01 were included in the further studies. After removal of the un-annotated genes, a total of 360 and 202 HCC up-and down-regulated genes were subjected to the subsequent gene ontology analyses. The two lists of the differentially expressed genes were analyzed using the online software FatiGo for comparative gene ontology categories including molecular function, biological process and cellular component [29] . They were also imported into the on-line software KEGG2 for pathway mapping [30, 31] . The statistical significance was defined as P < 0.01 between the HCC up-and down-regulated gene groups using Fisher's exact test. Quantitative monitoring of gene expression patterns with a complementary DNA microarray Expression monitoring by hybridization to high-density oligonucleotide arrays Somogyi R: Large-scale temporal gene expression mapping of central nervous system development A DNA microarray survey of gene expression in normal human tissues Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries Suppression subtractive hybridization: a versatile method for identifying differentially expressed genes Differential screening and suppression subtractive hybridization identified genes differentiallyexpressed in an estrogen receptorpositive breast carcinoma cell line Characterization of minute differences between genomes of strains of Penicillium nalgiovense using subtractive suppression hybridization without cloning Suppression subtractive hybridization Combining SSH and cDNA microarrays for rapid identification of differentially expressed genes Detection of differentially expressed genes in primary tumor tissues using representational differences analysis coupled to microarray hybridization Tissue-wide expression profiling using cDNA subtraction and microarrays to identify tumor-specific genes Generation of subspecies level-specific microbial diagnostic microarrays using genes amplified from subtractive suppression hybridization as microarray probes Enhanced microarray performance using low complexity representations of the transcriptome Gene expression profiling of osteoclast differentiation by combined suppression subtractive hybridization (SSH) and cDNA microarray analysis Suppression subtractive hybridization and microarrayidentification of estrogen-regulated hypothalamic genes Suppression subtractive hybridization coupled with microarray analysis to examine differential expression of genes in virus infected cells Dys-regulation of clusterin in human hepatoma is not associated with tumorigenesis but is secondary to cell response to external tresses Differential Gene Expression Identified in Complete Hydatidiform Mole by Combining Suppression Subtractive Hybridization and cDNA Microarray Profile of differentially expressedgenes after transfer of chromosome 17 into the breast cancer cell line CAL51 Analysis of calciuminducible genes in keratinocytes using suppression subtractive hybridizationand cDNA microarray Identification of novel and known oocyte-specificgenes using complementary DNA subtraction and microarray analysis inthree different species Murphy D: Microarray screening of suppression subtractive hybridization-PCR cDNA libraries identifies novel RNAs regulated by dehydration in the rat supraoptic nucleus Identification of target genes for the Ewing's sarcoma EWS/FLI fusion protein by representational difference analysis Functional genomics approach to hypoxia signaling Polyadenylation of the mRNA of hepatitis delta virus is dependent on the structure of the nascent RNA and regulated by the small or large delta antigen Aberrant Caspase-activated DNase (CAD) gene in human hepatoma cells Molecular signature of clinical severity in recovering patients with severe acute respiratory syndrome coronavirus (SARS-CoV) FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes Reconstruction of amino acid biosynthesis pathways from the complete genome sequence KEGG 2: Kyoto Encyclopedia of Genes and Genomes All authors read and approved the final manuscript. This work was supported by a research grant from Chang Gung Memorial Hospital (CMRPG), Tao-Yuan, Taiwan. S-Y Hsieh conceived and designed this study and was responsible for the manuscript preparation. The list of the 69 and 701 differentially expressed genes. A total of 69 and 701 differentially expressed genes identified via conventional cDNA microarray and the modified SSH/microarray approaches, respectively, are provided (Pan YS et al., 2006) . Click here for file [http://www.biomedcentral.com/content/supplementary/1471-2164-7-131-S1.pdf]