CAGE-TSSchip: promoter-based expression profiling using the 5'-leading label of capped transcripts
- Shintaro Katayama†1,
- Mutsumi Kanamori-Katayama†1,
- Kazumi Yamaguchi1, 2,
- Piero Carninci1, 3 and
- Yoshihide Hayashizaki1, 3Email author
© Katayama et al.; licensee BioMed Central Ltd. 2007
Received: 4 October 2006
Accepted: 26 March 2007
Published: 26 March 2007
Cap analysis gene expression (CAGE) technology has revealed numerous transcription start sites (TSSs) in mammals and has suggested complex promoter-based patterns of regulation. We developed the CAGE-TSSchip to detect promoter-based transcriptional activity. The CAGE-TSSchip is a customized oligonucleotide array that targets known TSSs identified by CAGE. A new labeling method, labeling capped transcripts from the 5'-end, had to be developed. The CAGE-TSSchip is accurate and sensitive, and represents the activity of each TSS.
Many genome sequencing projects of model species are finished and a large number of full-length cDNAs have been isolated. Trends in large-scale life science are changing from collection of essential elements to developing an understanding of global biologic mechanisms. Transcriptional regulatory pathways are among the basal functional mechanisms that remain largely unknown; estimation of promoter activity is an essential component of analysis of regulatory networks. Large-scale analysis of the human and mouse transcriptomes using cap analysis gene expression (CAGE) technology , revealed numerous transcription start sites (TSSs) [2, 3]. The TSSs are not randomly distributed; rather, they are concentrated at several short regions connected to each gene. On average there are five or more TSS clusters at one locus, and these are not only at the 5'-end of the gene but also within the open reading frame or 3'-untranslated region (UTR). Promoter-based expression clustering revealed that even TSS clusters in the same locus exhibit different expression patterns. This finding implies that the regulatory mechanism is defined by each TSS cluster. Measuring the transcriptional activity by using TSSs rather than genes would therefore lead to a better understanding of transcriptional regulatory mechanisms. Furthermore, promoter-based expression profiling is of benefit to the research community.
A tag-based approach for TSS analysis  such as CAGE requires deep sequencing when it is used to measure fluctuations in transcript expression, but deep sequencing is time consuming and expensive. Also, the various traditional expression profiling technologies did not represent the activity of each TSS but only the total activity of some TSSs. Searching among the microarray technologies for a technique that will permit large-scale promoter-by-promoter analysis, we modified our mature technology of purifying capped transcripts  and developed a new labeling method starting from the 5'-end of capped transcripts. This protocol made it possible for us to design an array for promoter-based expression profiling, which we named the CAGE-defined TSS chip (CAGE-TSSchip). We demonstrated its accuracy and sensitivity. Furthermore, by using CAGE-TSSchip we were able to predict principal regulatory factors.
Results and discussion
CAGE-TSSchip for mouse promoters
We isolated three total RNAs from mouse and conducted two comparisons using the CAGE-TSSchip; adult mouse liver versus mouse whole embryo in Theiler stage 17.5 (E17.5), and hepatocellular carcinoma cell line Hepa1-6 versus adult mouse healthy liver. We synthesized labeled RNAs using our 5'-leading method of capped transcripts and hybridized them to the CAGE-TSSchip. To estimate the reproducibility of our protocols, we designed dye swap experiments for these two comparisons. These experiments also helped us to reduce unavoidable technical variation . After elimination of control, non-uniform, non-significant, or saturated spots, we deleted the hybridization signal that did not exhibit similar values in each dye swap experiment. The scatter plots for each dye swap experiment revealed good correlation (r = 0.87-0.96; Additional data file 2). The variation caused by procedures (described in Materials and methods, below) including our 5'-leading label method is therefore small.
Accuracy and sensitivity: similar tendencies with qRT-PCR and CAGE
Cross-validation by qRT-PCR and CAGE in mouse liver versus E17.5
The frequency of CAGE tags reflects the activity of TSSs . We examined whether this TSS activity shown by CAGE was reflected in the CAGE-TSSchip. We counted CAGE tag numbers in liver and E17.5 at the region upstream from the CAGE-TSSchip probe position (see Materials and methods, below). We focused on 20 probes, which once again were compared with qRT-PCR. In this comparison CAGE tags corresponding to 17 probes exhibited similar positive log ratios, and two of the three remaining probes exhibited similar negative log ratios (Table 1). Therefore, the CAGE-TSSchip also shows an expression ratio similar to the frequency identified by CAGE tag.
CAGE or similar serial analysis technologies require deep sequencing if they are to recognize fluctuations in weak promoter activity. Therefore, sensitivity is an important issue for the CAGE-TSSchip. To estimate sensitivity, we evaluated whether results with the CAGE-TSSchip and the corresponding qRT-PCR were similar even when promoter activity is low. First, we selected some CAGE-TSSchip probes, without considering the log ratio values in the liver versus E17.5 comparison, and designed 88 primers (see Materials and methods, below) corresponding to these probes. We then measured expression levels using qRT-PCR and compared expression ratios (Additional data file 3). In this comparison we could identify a tendency toward large mathematical error (difference) in the log ratio between the CAGE-TSSchip and qRT-PCR at high maximum qRT-PCR Ct value in liver and E17.5 (Additional data file 4a). These findings mean that the log ratios of rare transcripts tend to differ between the two methods. This is intuitive because such large Ct values in qRT-PCR, especially 30 or greater, also exhibit technical variations in repetitive experiments. In our experiments, the Ct value of 30 is equal to one transcript per eight cells. However, the log ratios in the liver versus E17.5 comparison were well correlated between CAGE-TSSchip and qRT-PCR (r = 0.77) in the 42 probes with a maximum Ct value above 30 (Additional data file 4b). About two million tags are required to recognize promoter-level fluctuations in expression of such rare transcripts (Ct value >30) when using CAGE; this imposes considerable burdens in terms of time and money. In conclusion, the CAGE-TSSchip is fast, has a good cost/performance ratio, and exhibits acceptable sensitivity.
Observations: intensity of the CAGE TSSchip represents the activity of each TSS
The discrepancy in minor promoters in liver between CAGE-TSSchip and qRT-PCR was expected because our labeling method involves the extraction of active promoters. The CAGE-TSSchip results also suggest that the regulatory mechanisms between these two promoters are different, even though they belong to the same gene. Findings of hierarchical clustering in CAGE expression  support this suggestion, because the expression patterns in that study were clearly different; the former promoter belongs to expression cluster number 4 and the latter to number 1 (Additional data file 6). Therefore, the CAGE-TSSchip findings in terms of these two isoforms are reasonable.
We then examined the Aldh7a1 gene, which encodes an important enzyme (aldehyde dehydrogenase 7 family, member A1) that is highly expressed in liver. As for Bdh1, there are two isoforms in Aldh7a1; however, the first exon of the short isoform shares the third exon of the long isoform. The CAGE-TSSchip has five probes for Aldh7a1 (Figure 3b and Additional data file 5). The CAGE-TSSchip findings suggest that the major promoter of Aldh7a1 in liver is the first exon of the short isoform; validation by qRT-PCR supports this finding.
Based on our design of the CAGE-TSSchip, we expected the intensity of CAGE-TSSchip findings to represent the activity of each TSS, which would lead to a considerably greater difference for 5'-side probes than for 3'-side ones. This tendency could be seen in Bdh, in Aldh7a1, and in other genes (for example Ppp3ca, Scp2, Glo1, Fga, and Trf). Because of this, we believe that the CAGE-TSSchip, as a tool for measuring TSS activity, will perform as we expected it to.
E2F target genes were activated in Hepa1-6
We wished to demonstrate whether the CAGE-TSSchip makes it possible to analyze promoter-based regulatory mechanisms directly. It is noted that the functional regulatory elements that control transcription tend to be located close to TSSs . We designed the CAGE-TSSchip probes at the proximal downstream region of known TSSs, and our protocols including the 5'-leading label method can demonstrate which TSSs are controlled. Because of this, the CAGE-TSSchip can help to identify important promoters and control elements. Below, we describe a comparison of Hepa1-6 and 'normal' adult mouse liver, and demonstrate both regulated (target) gene screening and regulator prediction.
In this comparison, 117 nonredundant probes of 98 genes identified over-expression (log ratio >0.5) in Hepa1-6, and 47 nonredundant probes of 36 genes revealed under-expression (log ratio <-0.5; Additional data file 7). In the comparison of the cancer cell line with normal tissue, many promoters of genes related to cell proliferation are expected to be extracted. Actually, genes related to DNA metabolism, which form the superclass of DNA replication in Gene Ontology (GO), were the most significantly enriched among the former genes (87/98 genes had some GO annotation and 21/87 genes were annotated with GO:0006259; P = 3.22 × e-07 using GOstat ). Moreover, there were probes targeting the 5'-UTR or almost upstream of the genes encoding anaphase-promoting complex subunit 5 (Anapc5), minichromosome maintenance proteins (MCM2-7), and cyclin-dependent kinase 4 (Cdk4) in the former gene set. MCM genes have recently emerged as cancer biomarkers , and Cdk4 is important for cell cycle G1 phase progression. It is no surprise that abnormal proliferation occurs in the comparison between Hepa1-6 and 'normal' liver. Therefore, these target gene screen findings with CAGE-TSSchip agree well with current findings.
Over-represented TFBSs around the over-expressed promoter in Hepa1-6 compared with liver
|Log10(Hepa1-6/liver)| > 0.5
Hepa1-6 > liver
Liver > Hepa1-6
Number of probes
Number of probes
2.31 × e-06
6.04 × e-06
9.96 × e-05
2.50 × e-04
3.72 × e-04
We note that the first exon of E2F transcription factor 7 (E2F7) was over expressed in Hepa1-6. As de Bruin and coworkers  pointed out, this gene could block the E2F-dependent activation of a subset of E2F target genes. Zfp161 and Churc1 are novel candidate regulators of Hepa1-6 over-expressed genes because the CAGE-TSSchip analysis revealed that these TFBSs are also over-represented around the Hepa1-6 over-expressed promoters. These novel regulators might represent an alternative regulatory pathway for Hepa1-6 phenotype.
We believe CAGE-TSSchip to be a useful tool in promoter-by-promoter analysis of regulatory networks. When similar prediction is performed using a non-promoter-specific microarray-based gene expression technology, representative transcripts (for example, RefSeq sets) are used to identify the genomic region that regulates promoter activity. The 5'-end of these representative transcripts is assumed to be the candidate TSS. Furthermore, proximal regions of these TSSs are candidate regulatory regions of transcription when this type of technology is used. If this traditional technology yields a similar result, then the regulated TSSs identified by the CAGE-TSSchip should overlap with the 5'-end of the RefSeq transcripts. However, out of the 163 TSSs belonging to over-expressed or under-expressed genes, 74 did not overlap with the 5'-end of the RefSeq sets. All of the cDNAs can be used to capture all of the TSSs; in this case, many unregulated TSSs would be included. For example, the probe set of the Affymetrix MG-U74 v2 chip could not define the singular TSSs in 26 out of the 124 genes that exhibited over-expressed or under-expression. Such probes show the summation of activities in all alternative promoters, and the search space for regulatory elements expands. Therefore, although the prediction of important regulators using the traditional expression profiling technology might be able to achieve similar results as the CAGE-TSSchip, one could assume that the significance would be lower. The CAGE-TSSchip has been optimized for promoter-by-promoter analysis.
We developed the CAGE-TSSchip technology. This chip was able to identify the probes targeting the proximal region of the promoter defined by CAGE, and must be used with a new labeling method. This labeling method permitted labeling from the 5'-end of the capped transcripts. In a direct comparison between mouse liver and E17.5, CAGE-TSSchip identified expression ratios similar to those with qRT-PCR and CAGE, and had sufficient sensitivity to recognize the fluctuation in rare transcripts. Furthermore, the intensities of CAGE-TSSchip reflected the activity of each TSS, and so this technology may be useful in evaluating regulatory pathways.
CAGE-TSSchip permits promoter-based expression profiling with a favorable ratio of cost to performance and good accuracy by applying mature, two-color microarray technologies and equipment. Recently, several microarray platforms supporting one-color gene expression analysis for comparisons of many samples were developed. We were unfortunately unable to try these systems, but we will be able to change CAGE-TSSchip to a one-color analysis with minor modification.
In CAGE  and similar serial analysis technologies  for identification of novel TSSs, deep sequencing is necessary to identify promoters of rare transcripts or to compare expression levels in several samples. The current CAGE-TSSchip cannot identify novel promoters because we designed probes based on known transcripts and promoters, mainly defined by CAGE. However, several high-density microarray technologies will help us to identify novel promoters by combining them with our 5'-leading label method. A whole-genome tiling array is one approach to genome-wide promoter-based expression profiling.
An initial step in the analysis of transcriptional regulatory mechanisms is the identification of regulated elements and control elements. TSSs are just regulated elements, and major control elements are located around them. Therefore, promoter-based expression profiling is important in enhancing our understanding of regulatory mechanisms. CAGE-TSSchip and our 5'-leading label method is an alternative approach to promoter-based expression profiling, and it will help us to conduct promoter-by-promoter analysis of regulatory networks.
Materials and methods
Figure 2 is an overview of the CAGE-TSSchip probe design. First, we defined tag clusters (TCs) from transcripts and several tag-based resources. Furthermore, we chose the representative position of a TC as the most frequent TSS. We selected about 4,500 TCs from about 2,500 transcriptional units , which were mainly transcription factors  as defined by CAGE tags from E17.5. We then prepared 120-nucleotide long genomic sequence fragments located in the proximity downstream of the representative position of the TC, according to splicing patterns of known transcripts. Custom Microarray Design Services (Agilent Technologies, Santa Clara, California, United States)) proposed appropriate 60-mer probes from each fragment. We adopted two (redundant) probes from each fragment, and added several control probes and reference probes (reverse complement to the Agilent Catalog Array probes; the prefix of the probeID is 'A_'). All probe sequences and their annotations are available in Additional data file 1, and details of the probe design are available in Additional data file 8.
Tissues from adult male and embryos from C57BL/6J mice were extracted according to the RIKEN Institute's guidelines. The Hepa1-6 cell line was purchased from the RIKEN Cell Bank (Tsukuba, Ibaraki, Japan) and was cultured in Dulbecco's modified eagle medium supplemented with 10% heat-inactivated fetal bovine serum, 200 U/ml penicillin, and 200 μg/ml streptomycin. The total RNA was extracted using the acid phenol guanidinium thiocyanate-chloroform method. Details of the RNA preparation are available in Additional data file 8.
5'-Leading label and hybridization
Figure 1 shows the schematic procedure of the 5'-leading label and hybridization process. The cDNA synthesis was performed using 50 μg of total RNA and with first-strand cDNA primer (random sequence) using SuperScript II RT (Invitrogen, Carlsbad, California, United States). The full-length cDNAs were then selected with the biotinylated cap-trapper method. A specific linker was used that contained the T7 promoter sites 'upper oligonucleotide GN3' (sequence 5'-ACTAATACGACTCACTATAGGNNN-3') and 'upper oligonucleotide GGN2' (sequence 5'-ACTAATACGACTCACTATAGGGNN-3'), which were mixed at a ratio of 4:1. This mixture was in turn mixed at a ratio of 1:1 to the 'lower oligonucleotides' (sequence 5'-TGATTATGCTGAGTGATATCC-3'). The sequence was then ligated to the single-strand cDNA. The second strand of the cDNA was synthesized with the T7 promoter primer and the DNA polymerase (TaKaRa, Ohtsu, Shiga, Japan). Details of cDNA synthesis, cRNA amplification for the 5'-leading label, and the hybridization are available in Additional data file 8.
Quality check for the CAGE-TSSchip assay
Before analysis, control spots, saturated spots, non-uniform spots, and non-significant spots (according to Feature Extraction, the standard tool provided by Agilent for evaluating probe features) were removed. We also compared the Cy3 and Cy5 intensities of the same RNA samples in a dye swap experiment. We expected these signals to be correlated; however, low-intensity spots diverged somewhat from the regression line. We therefore excluded such probes that differed more than two times the standard residual deviation from the regression line. All intensity values and filtering results are available in Additional data file 1, and an overview can be found in Additional data file 2.
Validation with qRT-PCR
Primer pairs were designed using an optimal primer size of 20 bases and annealing temperature of 60°C, using Primer3 software . The uniqueness of the designed primers pairs was verified using the UCSC in silico PCR search in the UCSC Genome Browser Database . This method checks that homologous regions are not cross-amplified by the same primer pair. All primers were also verified by amplification with mouse genome DNA. First-strand cDNA synthesis (5 μg total RNA per 20 μl reaction) was carried out using a random primer and the ThermoScript RT-PCR System (Invitrogen), in accordance with the manufacturer's protocol. A qRT-PCR was carried out with first-strand cDNA corresponding to 12.5 ng total RNA per reaction well using the tailor-made reaction . The PCR reactions were performed with an ABI Prism (Applied Biosystems, Foster City, California, United States) using the following cycling protocols: 15 min hot start at 94°C, followed by 40 cycles of 15 s at 94°C, 30 s at 60°C, and 30 s at 72°C. The threshold cycle (Ct) value was calculated from amplification plots, in which the fluorescence signal detected was plotted against the PCR cycle. All primer sequences are available in Additional data files 3, 5 and 9.
Validation with CAGE
Transcripts overlapping with probes serve as guides for the assignment between probes and CAGE tags. The total number of CAGE tags located from the probe position to 100 nucleotides upstream of the 5'-end of the overlapping transcripts is the expression level as estimated by CAGE. If several transcripts overlap with the same probe, then the transcript transcribed from the most upstream position is chosen as a representative transcript. CAGE tags are classified by RNA samples. The target RNA library IDs in this study were CBR, CCM and IN, corresponding to liver, Hepa1-6 and E17.5, respectively. Finally, log ratio values were calculated according to CAGE-TSSchip assays. Dataset details are available in Additional data file 8.
Over-represented TFBS analysis
First, we chose probes exhibiting significant differences between Hepa1-6 and liver, with an absolute ratio above 0.5. After exclusion of redundant probes, we predicted the TFBSs around the probes in an area ranging from 1,000 nucleotides upstream to 200 nucleotides downstream using MATCH  from TRANSFAC  9.4, with minimum false-negative profiles (minFN94.prf). The over-representation of each binding matrix was evaluated by using Fisher's exact probability test . The matrices in Table 2 are the five most significantly over-represented ones in the regulatory regions of several genes, which exhibit higher expression in Hepa1-6 than in liver.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 provides TSSchip probe annotation and experimental results. Additional data file 2 shows the performance of 5'-leading label in dye swap experiments. Additional data file 3 provides details of sensitivity check with qRT-PCR. Additional data file 4 summarizes the sensitivity check with qRT-PCR. Additional data file 5 provides details of the alternative promoter check with qRT-PCR. Additional data file 6 provides CAGE expression clustering results of Bdh alternative promoters. Additional data file 7 summarizes over-expressed promoters in Hepa1-6 and liver. Additional data file 8 provides supplementary methods about the array probe design and whole protocols of wet experiments. Additional data file 9 gives details of cross-validation by qRT-PCR and CAGE in mouse liver versus E17.5.
We thank Yuki Tsujimura for her technical assistance, Yasumasa Kimura for the CAGE analysis, Noriko Ninomiya for the qRT-PCR analysis, Yayoi Kitazume for sample preparation, and Ann Karlsson and Hanna Daub for English editing. This study was supported by research grants from the Ministry of Education, Culture, Sports, Science and Technology of the Japanese Government to YH for the following: (1) the Genome Network Project and (2) the RIKEN Genome Exploration Research Project.
- Shiraki T, Kondo S, Katayama S, Waki K, Kasukawa T, Kawaji H, Kodzius R, Watahiki A, Nakamura M, Arakawa T, et al: Cap analysis gene expression for high-throughput analysis of transcriptional starting point and identification of promoter usage. Proc Natl Acad Sci USA. 2003, 100: 15776-15781. 10.1073/pnas.2136655100.PubMedPubMed CentralView ArticleGoogle Scholar
- Carninci P, Sandelin A, Lenhard B, Katayama S, Shimokawa K, Ponjavic J, Semple CA, Taylor MS, Engstrom PG, Frith MC, et al: Genome-wide analysis of mammalian promoter architecture and evolution. Nat Genet. 2006, 38: 626-635. 10.1038/ng1789.PubMedView ArticleGoogle Scholar
- Carninci P, Kasukawa T, Katayama S, Gough J, Frith MC, Maeda N, Oyama R, Ravasi T, Lenhard B, Wells C, et al: The transcriptional landscape of the mammalian genome. Science. 2005, 309: 1559-1563. 10.1126/science.1112014.PubMedView ArticleGoogle Scholar
- Harbers M, Carninci P: Tag-based approaches for transcriptome research and genome annotation. Nat Methods. 2005, 2: 495-502. 10.1038/nmeth768.PubMedView ArticleGoogle Scholar
- Katayama S, Hayashizaki Y: Complex transcription mechanisms in mammalian genomes: the transcriptome of FANTOM3. Curr Genomics. 2005, 6: 619-625. 10.2174/138920205775811443.View ArticleGoogle Scholar
- Carninci P, Kvam C, Kitamura A, Ohsumi T, Okazaki Y, Itoh M, Kamiya M, Shibata K, Sasaki N, Izawa M, et al: High-efficiency full-length cDNA cloning by biotinylated CAP trapper. Genomics. 1996, 37: 327-336. 10.1006/geno.1996.0567.PubMedView ArticleGoogle Scholar
- Carninci P, Westover A, Nishiyama Y, Ohsumi T, Itoh M, Nagaoka S, Sasaki N, Okazaki Y, Muramatsu M, Schneider C, et al: High efficiency selection of full-length cDNA by improved biotinylated cap trapper. DNA Res. 1997, 4: 61-66. 10.1093/dnares/4.1.61.PubMedView ArticleGoogle Scholar
- Matsumura H, Bin Nasir KH, Yoshida K, Ito A, Kahl G, Kruger DH, Terauchi R: SuperSAGE array: the direct use of 26-base-pair transcript tags in oligonucleotide arrays. Nat Methods. 2006, 3: 469-474. 10.1038/nmeth882.PubMedView ArticleGoogle Scholar
- Kanamori M, Konno H, Osato N, Kawai J, Hayashizaki Y, Suzuki H: A genome-wide and nonredundant mouse transcription factor database. Biochem Biophys Res Commun. 2004, 322: 787-793. 10.1016/j.bbrc.2004.07.179.PubMedView ArticleGoogle Scholar
- Kawaji H, Kasukawa T, Fukuda S, Katayama S, Kai C, Kawai J, Carninci P, Hayashizaki Y: CAGE Basic/Analysis Databases: the CAGE resource for comprehensive promoter analysis. Nucleic Acids Res. 2006, D632-D636. 10.1093/nar/gkj034. 34 DatabaseGoogle Scholar
- Churchill GA: Fundamentals of experimental design for cDNA microarrays. Nat Genet. 2002, 490-495. 10.1038/ng1031. 32 SupplGoogle Scholar
- Wasserman WW, Sandelin A: Applied bioinformatics for the identification of regulatory elements. Nat Rev Genet. 2004, 5: 276-287. 10.1038/nrg1315.PubMedView ArticleGoogle Scholar
- Beissbarth T, Speed TP: GOstat: find statistically overrepresented Gene Ontologies within a group of genes. Bioinformatics. 2004, 20: 1464-1465. 10.1093/bioinformatics/bth088.PubMedView ArticleGoogle Scholar
- Gonzalez MA, Tachibana KE, Laskey RA, Coleman N: Control of DNA replication and its potential clinical exploitation. Nat Rev Cancer. 2005, 5: 135-141. 10.1038/nrc1548.PubMedView ArticleGoogle Scholar
- Chu S, Ferro TJ: Sp1: regulation of gene expression by phosphorylation. Gene. 2005, 348: 1-11. 10.1016/j.gene.2005.01.013.PubMedView ArticleGoogle Scholar
- Kageyama R, Merlino GT, Pastan I: A transcription factor active on the epidermal growth factor receptor gene. Proc Natl Acad Sci USA. 1988, 85: 5016-5020. 10.1073/pnas.85.14.5016.PubMedPubMed CentralView ArticleGoogle Scholar
- de Bruin A, Maiti B, Jakoi L, Timmers C, Buerki R, Leone G: Identification and characterization of E2F7, a novel mammalian E2F family member capable of blocking cellular proliferation. J Biol Chem. 2003, 278: 42041-42049. 10.1074/jbc.M308105200.PubMedView ArticleGoogle Scholar
- Kasukawa T, Katayama S, Kawaji H, Suzuki H, Hume DA, Hayashizaki Y: Construction of representative transcript and protein sets of human, mouse, and rat as a platform for their transcriptome and proteome analysis. Genomics. 2004, 84: 913-921. 10.1016/j.ygeno.2004.08.011.PubMedView ArticleGoogle Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.PubMedGoogle Scholar
- Hinrichs AS, Karolchik D, Baertsch R, Barber GP, Bejerano G, Clawson H, Diekhans M, Furey TS, Harte RA, Hsu F, et al: The UCSC Genome Browser Database: update 2006. Nucleic Acids Res. 2006, D590-D598. 10.1093/nar/gkj144. 34 DatabaseGoogle Scholar
- Suzuki H, Okunishi R, Hashizume W, Katayama S, Ninomiya N, Osato N, Sato K, Nakamura M, Iida J, Kanamori M, et al: Identification of region-specific transcription factor genes in the adult mouse brain by medium-scale real-time RT-PCR. FEBS Lett. 2004, 573: 214-218. 10.1016/j.febslet.2004.07.068.PubMedView ArticleGoogle Scholar
- Kel AE, Gossling E, Reuter I, Cheremushkin E, Kel-Margoulis OV, Wingender E: MATCH: A tool for searching transcription factor binding sites in DNA sequences. Nucleic Acids Res. 2003, 31: 3576-3579. 10.1093/nar/gkg585.PubMedPubMed CentralView ArticleGoogle Scholar
- Matys V, Kel-Margoulis OV, Fricke E, Liebich I, Land S, Barre-Dirrie A, Reuter I, Chekmenev D, Krull M, Hornischer K, et al: TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 2006, D108-D110. 10.1093/nar/gkj143. 34 DatabaseGoogle Scholar
- Fisher LD, Belle G: Biostatistics: a Methodology for the Health Sciences. 1993, New York, NY: John Wiley & SonsGoogle Scholar
- Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucleic Acids Res. 2006, D16-D20. 10.1093/nar/gkj157. 34 DatabaseGoogle Scholar
- Pruitt KD, Tatusova T, Maglott DR: NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2005, D501-D504. 33 DatabaseGoogle Scholar
- Birney E, Andrews D, Caccamo M, Chen Y, Clarke L, Coates G, Cox T, Cunningham F, Curwen V, Cutts T, et al: Ensembl 2006. Nucleic Acids Res. 2006, D556-D561. 10.1093/nar/gkj133. 34 DatabaseGoogle Scholar
- Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, Agarwal P, Agarwala R, Ainscough R, Alexandersson M, An P, et al: Initial sequencing and comparative analysis of the mouse genome. Nature. 2002, 420: 520-562. 10.1038/nature01262.PubMedView ArticleGoogle Scholar
- Kent WJ: BLAT: the BLAST-like alignment tool. Genome Res. 2002, 12: 656-664. 10.1101/gr.229202. Article published online before March 2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Florea L, Hartzell G, Zhang Z, Rubin GM, Miller W: A computer program for aligning a cDNA sequence with a genomic DNA sequence. Genome Res. 1998, 8: 967-974.PubMedPubMed CentralGoogle Scholar
- McGinnis S, Madden TL: BLAST: at the core of a powerful and diverse set of sequence analysis tools. Nucleic Acids Res. 2004, W20-W25. 10.1093/nar/gkh435. 32 Web ServerGoogle Scholar
- Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987, 162: 156-159. 10.1016/0003-2697(87)90021-2.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.