Skip Navigation


DNA Research Advance Access originally published online on July 8, 2008
DNA Research 2008 15(4):241-251; doi:10.1093/dnares/dsn015
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrowOA All Versions of this Article:
15/4/241    most recent
dsn015v2
dsn015v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Ohtsu, M.
Right arrow Articles by Murakami, Y.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ohtsu, M.
Right arrow Articles by Murakami, Y.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2008. Kazusa DNA Research Institute
The online version of this article has been published under an open access model. Users are entitled to use, reproduce, disseminate, or display the open access version of this article for non-commercial purposes provided that: the original authorship is properly and fully attributed; the Journal and Oxford University Press are attributed as the original place of publication with the correct citation details given; if an article is subsequently reproduced or disseminated not in its entirety but only in part or as a derivative work this must be clearly indicated. For commercial re-use, please contact journals.permissions@oxfordjournals.org

Novel DNA Microarray System for Analysis of Nascent mRNAs

Masaya Ohtsu1, Mika Kawate1, Masashi Fukuoka1, Wataru Gunji2, Fumio Hanaoka3, Takahiko Utsugi2, Fumitoshi Onoda1 and Yasufumi Murakami1,2,*

1 Faculty of Industrial Science and Technology, Department of Biological Science and Technology, Tokyo University of Science, 2641 Yamazaki, Noda-shi, Chiba 278-8510, Japan
2 Bio Matrix Research Inc., 105 Higashifukai, Nagareyama, Chiba 275-0101, Japan
3 Graduate School of Frontier Biosciences, Osaka University, 1-3 Yamada-oka, Suita, Osaka 565-0871, Japan

Received 3 May 2008; accepted 11 June 2008.


   Abstract
Top
 Abstract
1. Introduction
 2. Materials and Methods
 3. Results and Discussion
 Supplementary Data
 Funding
 References
 
Transcriptional activation and repression are a key step in the regulation of all cellular activities. The development of comprehensive analysis methods such as DNA microarray has advanced our understanding of the correlation between the regulation of transcription and that of cellular mechanisms. However, DNA microarray analysis based on steady-state mRNA (total mRNA) does not always correspond to transcriptional activation or repression. To comprehend these transcriptional regulations, the detection of nascent RNAs is more informative. Although the nuclear run-on assay can detect nascent RNAs, it has not been fully applied to DNA microarray analysis. In this study, we have developed a highly efficient method for isolating bromouridine-labeled nascent RNAs that can be successfully applied to DNA microarray analysis. This method can linearly amplify small amounts of mRNAs with little bias. Furthermore, we have applied this method to DNA microarray analysis from mouse G2-arrested cells and have identified several genes that exhibit novel expression profiles. This method will provide important information in the field of transcriptome analysis of various cellular processes.


   1. Introduction
Top
 Abstract
 1. Introduction
2. Materials and Methods
 3. Results and Discussion
 Supplementary Data
 Funding
 References
 
DNA microarray is an important tool for understanding regulatory networks. In fact, this technique was used to analyze changes in the amounts of mRNA in cellular phenomena such as cell differentiation, cellular senescence, and cell cycle progression in a comprehensive manner. Moreover, we have successfully improved the sensitivity and reproducibility of DNA microarray analysis.1Go Although many of the expression profile data obtained through DNA microarray analysis are available at various websites, such data are based on steady-state mRNA levels. It is obvious that both mRNA synthesis and degradation influence steady-state mRNA levels, but usually only the total mRNA is quantified. If we can quantify nascent mRNA in a real-time manner, it will become possible to estimate mRNA synthesis and degradation rates by comparing the nascent amount with the total amount of RNA measured by an ordinary DNA microarray system. One method of detecting nascent mRNA is the nuclear run-on assay.2Go Recently, the transcriptional profiling of radio-labeled RNAs using the nylon-membrane DNA microarray was reported.3Go–5Go The current standard platforms are GeneChip by Affymetrix, the Stanford-type DNA microarray, and the oligo-DNA microarray. In general, it is very difficult to compare data from such standard DNA microarray systems with data obtained by the nylon membrane DNA microarray. To comprehensively analyze nascent RNAs using these current standard platforms, improved methods for the isolation and labeling of nascent mRNAs were reported. Cleary et al. used uracil phosphoribosyltransferase gene-transformed human cells in Toxoplasma gondii. Nascent RNAs from these cells grown in 2,4-dithiouracil were labeled with thio-substrated UTP, purified, and analyzed by DNA microarray.6Go However, it is necessary to introduce this gene to the cell to be analyzed prior to DNA microarray analysis.

Bromouridine, which substitutes 5' of uridine to bromine, is inexpensive and becomes incorporated into cells with very high efficiency. Additional steps, such as electroporation or lipofection, are not necessary for labeling mRNA with this compound. After its incorporation into cells, bromouridine is converted to Br-UTP. During transcription, converted Br-UTP is recognized as the same substrate as UTP and is incorporated into nascent RNAs.7Go This labeling method was applied in the immunocytochemical analysis by using anti-bromodeoxyuridine (BrdU) antibody in several studies.8Go–11Go That antibody was used also to immunoprecipitate this BrU-labeled nascent RNA.12Go,13Go The same reports analyzed transcriptional changes of some genes in small amounts (on the order of the picogram) of nascent RNAs.12Go,13Go

For a comprehensive analysis using DNA microarray, we have optimized the conditions for immunoprecipitation using BrdU antibody. We have improved the specificity of nascent mRNA isolation by orders of magnitude using adequate blocking agents in immunoprecipitation. As the amount of nascent RNA obtained was very small even under the optimized purification conditions, we sought to develop a method for amplifying small amounts of mRNA with little bias and have succeeded in this goal. Finally, we have performed transcriptional analysis of nascent RNAs using G2 phase-synchronized mouse mammary carcinoma FM3A cells, and we have obtained several interesting insights into cell cycle regulation of the G2 phase.

These results demonstrate that DNA microarray analysis with nascent RNAs obtained by our new method will provide valuable insights into the functions of genes in various types of cells.


   2. Materials and Methods
Top
 Abstract
 1. Introduction
 2. Materials and Methods
3. Results and Discussion
 Supplementary Data
 Funding
 References
 
2.1. Preparation of BrU-labeled eGFP cRNA and non-labeled luciferase cRNA
eGFP cRNA and luciferase cRNA were prepared by T7 in vitro transcription. DNA templates for each cRNA synthesis were constructed by PCR amplification from plasmid DNA containing eGFP or luciferase genes. These templates contain both T7 promoter sequence and polyA sequence. The DNA template of eGFP was amplified from pEGFP-c1 by PCR. The DNA template of luciferase was amplified from pTRE (Clontech, Mountain View, CA, USA) by PCR. Sequences of primers were described in ‘Supplementary data’. The PCR products were purified from agarose gel by using Wizard SV gel and PCR purification kit (Promega, Madison, WI, USA).

cRNAs were transcribed with Br-UTPs (Sigma-Aldrich, St Louis, MO, USA) and NTPs using the MAXIscript® T7 kit. cRNAs were purified according to the RNA-cleanup protocol of the RNeasy Mini kit (Qiagen, Hilden, Germany).

2.2. Preparation of Mouse anti-BrdU IgG binding dynabeads
Two micrograms of mouse anti-bromodeoxyuridine antibody (Roche Diagnostics, Indianapolis, IN, USA) were incubated with 25 µl Dynabeads® Goat anti-mouse IgG (Invitrogen, Carlsbad, CA, USA) in 2.0 ml collection tubes containing 100 µl DEPC-treated phosphate-buffered saline (PBS)/0.1% bovine serum albumin (BSA) solution. The tubes were rotated at room temperature for an hour. After the dynabeads were collected by a magnet rack, the dynabeads were washed three times with 1 ml DEPC-treated PBS/0.1% BSA, and 100 µl DEPC-treated PBS /0.1% BSA was added.

2.3. Immunoprecipitation of BrU-labeled RNA by antibody beads
Following steps were conducted in the dark. Two hundred nanograms of BrU-labeled eGFP cRNA and 200 ng of non-labeled luciferase cRNA were denatured at 80°C for 10 min. As the blocking agent, 20 µg FM3A total RNA or 200 µg Escherichia coli 16S and 23S ribosomal RNA (rRNA) or uridine (final concentration of 0.3 M) was added. The denatured RNAs were added to the beads containing 225 U/ml RNasin® Plus RNase inhibitor (Promega). PBS(-)/0.1% BSA was added to this solution to a volume of 250 µl. The beads were rotated for 1 h at room temperature. After the beads were washed with 0.8 ml DEPC-treated PBS/0.1% BSA containing 45 U/ml RNasin® Plus RNase inhibitor (Promega), they were also washed with 0.8 ml DEPC-treated PBS three times. After the addition of 10 µl of RNase-free water, RNAs which bound to beads were eluted by denaturation at 80°C for 10 min, and then were kept cold on ice for 2 min.

2.4. Reverse transcription of cRNA and Real-time RT–PCR
2.4.1. Reverse transcription of cRNA
Five microliters of the isolated RNAs were mixed with 1 µl of 50 µM oligo(dT)20 primer and denatured at 65°C for 5 min. First-strand cDNAs were synthesized from the isolated RNAs using the SuperScript® III First-Strand Synthesis System for RT–PCR (Invitrogen). After reverse transcription, 1 U of RNase H was reacted for 20 min at 37°C.

2.4.2. Real-time quantitative RT–PCR
Real-time quantitative PCR was performed in an ABI prism 7900HT Sequence Detection System (Applied Biosystems, Foster City, CA, USA) using SYBR® Premix Ex TaqTM (Perfect Real Time) (Takara Bio, Otsu, Japan). Sequences of each primer were described in ‘Supplementary data’. The PCR cycling program was set as follows: 95°C for 15 s followed by 40 cycles of 95°C for 5 s and 60°C for 34 s. One of cDNAs diluted 100 times was used in a 50 µl reaction. A standard graph for eGFP cDNA and luciferase cDNA was simultaneously generated using 104, 105, 106, 107, and 108 copies of the pEGFP plasmid and the pTRE plasmid. Each cDNA was analyzed in triplicate. Three experiments were conducted, and the averages of the cDNA copy numbers as well as standard deviation were calculated.

2.5. Cell culture of FM3A cells and labeling nascent RNA with BrU
Mouse FM3A cells and tsFT210 cells, Cdc2 temperature-sensitive mutant strain of FM3A cells,14Go,15Go were cultured in RPMI 1640 medium (Invitrogen) with 10% dialyzed bovine serum, 2 mM L-glutamine, 0.1 mM MEM non-essential amino acid (Invitrogen), and 50 U/ml penicillin–50 mg/ml streptomycin (Invitrogen) at 32°C in 5% CO2.

After bromouridine was added to the medium to a final concentration of 100 µM, FM3A cells were cultured for an hour in the dark. FM3A cells were collected in 50 ml centrifuge tubes and washed twice with cold PBS, and total RNA was isolated by using an RNeasy midi kit (Qiagen). DNase I digestion was conducted during isolation of RNA.

2.6. Immunoprecipitation of nascent RNA from mouse FM3A cells
Twenty micrograms of total RNA from bromouridine-treated FM3A were immunoprecipitated with 200 µg of E. coli 16S and 23S rRNAs. The method of immunoprecipitation was described above. The amount of the eluted RNA was measured using the Quant-iTTM RiboGreen® RNA quantitation kit (Invitrogen). Five hundred nanometers of fluorescence of the RiboGreen reagent bound to RNA were measured using the fluorescence microplate reader ARVO (Perkin Elmer, Boston, MA, USA).

2.7. Spike-in experiment
Both BrU-labeled eGFP cRNA and non-labeled luciferase cRNA were added to 20 µg of total RNA from bromouridine-treated FM3A cells. We conducted three different kinds of experiments. Each experiment was performed with a different amount of cRNAs (2 ng, 200 pg, or 20 pg). These RNAs were immunoprecipitated with 200 µg E. coli 16S and 23S rRNAs as blocking agents.

2.8. Validation with RT–PCR of the DNA microarray data using each amplification method
Methods of experiences were described in ‘Supplementary data’.

2.9. Preparation of DNA Microarray
Oligonucleotide probe sets of Mus musculus AROS Version 3.0 (Operon Biotechnologies, Huntsville, AL, USA) were dissolved at 10 µM in 3x SSPE solution (Invitrogen), and spotted on {gamma}-amino propyl silane-coated UltraGAPS Coated Slides (Corning, New York, NY, USA) and UV cross-linked as described previously.1Go

2.10. cDNA synthesis and amplification
The amplification of RNA was performed using the Super SMARTTM PCR cDNA Synthesis Kit (Clontech) as follows. We have modified the manufacture's protocol and optimized it for our DNA microarray system. Total RNA (20 ng) was mixed with 7 µl of 12 µM 3' SMART CDS Primer II A and 7 µl of 12 µM SMART II A Oligonucleotide, incubated at 65°C for 2 min, and cooled to room temperature. To this mixture was added 20 µl of 5x First-Strand Buffer, 2 µl of 100 mM DTT, 10 µl of 10 mM of 50x dNTP, and 5 µl of PowerScriptTM Reverse Transcriptase, and the reaction was incubated at 42°C for 90 min. After the reaction, 1 µl of 0.5 M EDTA was added. First-strand cDNA was purified with the NucleoSpin Extraction Kit (Clontech). For the amplification by long-distance PCR using the Advantage® 2 PCR Kit (Clontech), first-strand cDNA was mixed with 10 µl of 10x Advantage® 2 PCR Buffer, 2 µl of 10 mM 50x dNTP, 2 µl of 12 µM of 5'PCR Primer II A, and 2 µl of 50x Advantage® 2 Polymerase Mix. The PCR conditions consisted of a 95°C for 1 min, followed by appropriate number of cycles of amplification at 95°C for 15 s, 65°C for 30 s, and 68°C for 6 min. The appropriate number of cycling on which the PCR product will remain in the exponential phase of amplification was determined by the cDNA amplification results of the range of PCR cycles obtained in pre-experiments. The PCR product was purified with Wizard SV Gel and PCR Clean-Up System (Promega). For aminoallyl labeling, 1 µg cDNA was mixed with 15 µl of Random Primers Buffer Mixture (Invitrogen), incubated at 95°C for 5 min, and cooled on ice. Then, 5 µl of an aminoallyl–dNTP mixture which contained 5 mM aminoallyl–dUTP (Sigma), 5 mM dATP (GE healthcare, Uppsala, Sweden), 5 mM dGTP (GE healthcare), and 5 mM dCTP (GE healthcare) was added. Next 9 U of Klenow Fragment (Invitrogen) was added and incubated at 37°C for 120 min. After the reaction, 5 µl of Stop Buffer (Invitrogen) was added. Aminoallyl-labeled cDNA was purified with ethanol precipitation following the first purification with Wizard SV Gel and PCR Clean-Up System. For the coupling reaction, aminoallyl-labeled cDNA was dissolved with 5 µl of 0.1 M NaHCO3 and mixed with 5 µl of Cy3 Mono-reactive Dye (GE healthcare) or Cy5 Mono-reactive Dye (GE healthcare) in 50 µl Atlas Glass Approved DMSO (BD Bioscience), and incubated for 60 min at room temperature. After coupling, 15 µl of 4 M hydroxylamine–HCl (Sigma-Aldrich) was added and incubated at room temperature for 15 min. Cy-Dye labeled cDNA was purified with the MinElute Reaction Cleanup Kit (Qiagen).

2.11. Hybridization
These procedures were based on TF oligo microarray hybridization as described previously.1Go Some of the conditions and compositions of solutions were changed: the hybridization solution consisted of a mixture of Cy3- and Cy5-labeled cDNAs, 40% formamide, 5x SSC, 0.4% SDS, and 1.9 µg/µl acetylated BSA (Invitrogen). Hybridization was performed with a hybridization chamber (Agilent Technologies, Palo Alto, CA, USA) and a hybridization gasket slide (Agilent Technologies). The temperature in both the pre-hybridization and hybridization steps was 45°C. Hybridization images were scanned using an Agilent DNA microarray scanner (Agilent Technologies).

We hybridized four slides as self–self hybridization: Cy3-labeled cDNAs of nascent RNA from FM3A cells and Cy5-labeled cDNA of nascent RNA from FM3A cells.

In experiments using G2 phase-arrested cells, we hybridized four slides: Cy3-labeled cDNAs of asynchronous cells and Cy5-labeled cDNAs of G2 phase-arrested cells were hybridized in two slides, and dye-swapped cDNAs were hybridized in the other two.

2.12. DNA microarray analyses
Fluorescence intensities in spots were quantified using Feature Extraction software version GX (Agilent Technologies). Spots flagged as non-uniform outliners of either a feature (spot) or background were removed. In each array, genes whose expression levels were lower than those of the negative control spots under both conditions were removed. The ratios of signal intensities between the two cell types were calculated through Lowess normalization using GeneSpring GX (Agilent Technologies). We calculated the average and the standard deviation of the base 2 logarithm of the expression ratios of four spots and removed the expression ratios of genes whose standard deviation was >1. The expression ratios of genes whose average signal intensity was <50 were also removed.

2.13. Analyses of nascent RNA profiles in G2-arrested cells and asynchronous cells
tsFT210 cells were arrested in the G2 phase by incubation for 17 h at a non-permissive temperature (39°C). Bromouridine was added to the medium to a concentration of 100 µM at 16 h after the shift at 39°C. RNA was isolated using the RNeasy mini kit (Qiagen). DNase I digestion was conducted during isolation of RNA. The DNA microarray analyses are described above.

2.14. Cell cycle analysis
Cells were harvested, washed with PBS, and fixed in 70% ethanol at –20°C for at least 4 h. Fixed cells were washed twice with PBS and mixed with PBS containing 0.3% Triton X-100 (w/v) (Sigma-Aldrich). Cells were added with PBS (–) containing 50 µg/ml propidium iodide (Sigma-Aldrich) and 0.1 mg/ml RNase A (Sigma-Aldrich), and incubated for 20 min at 37°C in the dark before being filtered through a 41 µm nylon mesh. We analyzed the DNA contents by measuring the intensity of the fluorescence produced by propidium iodide using the FACSCalibur instrument (Beckon Dickinson, Mountain View, CA, USA) using instructions supplied by the manufacturer.


   3. Results and Discussion
Top
 Abstract
 1. Introduction
 2. Materials and Methods
 3. Results and Discussion
Supplementary Data
 Funding
 References
 
3.1. Experimental scheme of DNA microarray analysis for the nascent RNA
The new method for performing comprehensive analyses of nascent RNA is illustrated in Fig. 1. In the first step, bromouridine was added to the culture medium. This compound was incorporated into cells, which converted it to Br-UTP. As BrUTP is recognized as the same substrate as UTP, nascent mRNA is labeled by BrU7Go (Fig. 1). Total RNAs isolated from the BrU-treated cells were immunoprecipitated with anti-BrdU antibody beads.12Go,13Go The RNAs bound to the beads were eluted by heat denaturation after washing, and the eluted mRNAs were reverse-transcribed using oligo-dT primer. The cDNAs were linearly amplified using the template-switching (TS) PCR method. The amplified cDNAs labeled with Cy3/Cy5 were used for DNA microarray analysis (Fig. 1). Among these steps, the method for immunoprecipitating the BrU-labeled nascent RNAs was based on a method described previously12Go,13Go and optimized for DNA microarray analysis in the following manner. Approximately 50 pg of the BrU-labeled RNA from mouse pre-implantation embryos was obtained and subjected to the same analysis as used in a previous study.13Go The RNA amount should be >20 ng in total RNA. This is ~1000 times more than that in the previous study.


Figure 1
View larger version (30K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 1. Overview of the isolation of nascent RNAs for microarray analysis.

 
Therefore, it was necessary to develop a method for immunoprecipitating a large amount of BrU-labeled RNA with high specificity.

3.2. Determination of optimum conditions for enriching nascent RNAs using cRNAs synthesized in vitro
As the amount of BrU-labeled nascent RNAs is much smaller than that of non-labeled RNAs in the RNA fractions extracted from living cells, it was necessary to develop a method for obtaining BrU-labeled RNAs with high recovery as well as on a certain scale. In the first trial, a mixture of BrU-labeled eGFP cRNAs and non-labeled luciferase cRNAs prepared by T7 in vitro transcription was immunoprecipitated in the same tube. cRNAs eluted from the affinity beads were reverse-transcribed with oligo-dT primer, and the copy number of each cDNA was quantified by real-time quantitative PCR (Fig. 2A). We defined the ratio of the eGFP cDNA copy number to the luciferase cDNA copy number as the index of immunoprecipitation specificity. This index well reflects the enrichment efficiency of nascent mRNAs. We also examined the effects of adding blocking agents such as uridine, total RNA from FM3A cells, and E. coli 16S and 23S rRNAs. Although interferon-{tau} RNA was shown to work well as a blocker in a previous report,13Go we chose uridine, E. coli rRNA, and mouse total RNA as candidates for a blocking agent because we thought they would be more suitable for large-scale isolation of nascent RNAs. We evaluated the effects of these blocking agents with increasing amounts of cRNA. Under our conditions, the addition of uridine did not change the ratio. In contrast, the ratios observed in the experiments with E. coli rRNA or FM3A total RNA were approximately double those obtained in an experiment carried out without RNAs as the blocking agent (Fig. 2B). These data suggest that any kind of RNA can increase the specificity of the immunoaffinity purification of nascent RNAs. We thus chose E. coli 16S and 23S rRNAs as suitable blocking agents. These bacterial rRNAs are inexpensive and do not exist in mammalian cells. In addition, E. coli 16S and 23S rRNAs should not be reverse-transcribed when mRNA is reverse-transcribed by using oligo-dT primer.


Figure 2
View larger version (18K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 2. Binding specificity of BrU-labeled cRNA by in vitro transcription. (A) eGFP and luciferase gene were transcribed in vitro by T7 RNA polymerase. Br-UTP was used as the substrate instead of UTP during the transcription of eGFP. (B) Mixture of 200 ng of BrU-labeled eGFP and non-labeled luciferase cRNA was immunoprecipitated with addition to blocking agents. The eluted cRNAs were converted into cDNA, and then the cDNA copy number was measured using real-time quantitative PCR. The blocking agents used are as follows: (i) no addition, (ii) uridine, (iii) E. coli 16S and 23S RNAs, and (iv) total RNA from mouse FM3A cells. ‘Ratio (eGFP/Luci.)’ represents the ratio of the eGFP cDNA copy number to the luciferase copy number. Error bars represent standard errors of the means (n = 3).

 
3.3. Conditions for labeling nascent RNAs with BrU
Next, we treated living cells with bromouridine and immunoprecipitated BrU-labeled nascent RNAs from BrU-treated cells. We first tested and determined the appropriate concentration (100 µM) of bromouridine, i.e. that which did not affect proliferation and cell cycle progression in mouse FM3A cells (data not shown). Using DNA microarray analysis, we also confirmed that the gene expression profile of FM3A cells after treatment with 100 µM bromouridine for 4 h was extensively the same as that of untreated cells (data not shown). Approximately 20 µg total RNA was isolated from 2 x 106 FM3A cells that had been treated with 100 µM bromouridine for 1 h; this total RNA was then subjected to immunoaffinity purification. We used E. coli 16S and 23S rRNAs as blocking agents. The amount of eluted RNAs in the sample of BrU-treated cells was ~30 ng, so we had successfully obtained a sufficient amount of BrU-labeled RNA for DNA microarray analyses. We then evaluated the specificity of immunoaffinity purification.

3.4. Evaluation of purification efficiency with spike-in control
To more precisely estimate the efficiency with which nascent RNAs were purified, small amounts of T7 in vitro transcribed BrU-labeled eGFP cRNA and non-labeled luciferase cRNA were added to the RNAs obtained from BrU-treated FM3A cell samples as the ‘spike-in’ control. We immunoprecipitated these RNA mixtures and compared the cDNA copy numbers of eluted eGFP cDNA and luciferase cDNA. To determine whether or not the amount of the eluted BrU-labeled eGFP cRNA was proportional to that of the input eGFP cRNA, increasing amounts (20 pg, 200 pg, and 2 ng) of eGFP cRNA and luciferase cRNA were added to 20 µg of BrU-labeled total RNA obtained from FM3A cells (Fig. 3A). The ratio of eluted eGFP cDNA to eluted luciferase cDNA exceeded 83 under the conditions shown in Fig. 3B. In addition, the copy number of eluted BrU-labeled eGFP cRNA was proportional to that of the input eGFP cRNA. These results indicated that we were able to immunoaffinity-purify BrU-labeled RNA with high specificity and that the eluted RNA reflects the nascent RNA before immunoprecipitation. Taken together, these findings indicate that we have successfully optimized the conditions for immunoaffinity-purification with high specificity. It was already demonstrated that the enrichment process using anti-BrdU immunoaffinity beads does not have any bias.13Go As it is important to confirm that the cDNA amplification process does not have any significant bias, we performed the following experiments.


Figure 3
View larger version (27K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 3. Binding specificity in immunoprecipitation with spike-in controls. (A) Overview of the experiment. (B) Br-eGFP and luciferase cRNA were prepared as spike-in controls. ‘Ratio (eGFP/Luci.)’ represents the ratio of the eGFP cDNA copy number to the luciferase copy number. Error bars represent standard errors of the means (n = 3).

 
3.5. DNA microarray analyses of nascent RNAs
To carry out DNA microarray analysis, nascent RNAs of BrU-treated cells have to be amplified without bias. We previously developed a method for effectively amplifying small amounts of RNA using the TS-PCR reaction16Go (Fig. 4A). To examine whether or not this amplification method has any amplification bias, we compared the data obtained by DNA microarray analysis with those obtained by quantitative real-time PCR using RNA samples prepared from ES cells.1Go We considered that the efficiency of amplification would vary according to the differences in cellular RNA abundance. Genes with abundant RNA would be amplified without bias, but genes with low amounts of RNA might fail to be amplified accurately. As the expression levels of transcription-factor genes are generally much lower than those of genes in other categories, we focused on expression changes of the transcription factor genes. We used the TF oligo microarray, which can accurately determine expression changes of transcription factor genes.1Go As shown in Fig. 4B, the expression changes in TS-PCR amplification were detected accurately, without amplification bias, except for Prdx1, Rara. Especially, the fold changes of Dcx, Meis1, and Ttyh1 obtained by TS PCR were not only quite similar to the data obtained without amplification, but were also similar to the data from quantitative real-time PCR, whereas amplification by in vitro transcription with T7 promoter exhibited some discrepancies. Thus, amplification by TS PCR has very little bias (Fig. 4B).

However, there was a possibility that amplification efficiency and the reverse-transcription efficiency of BrU-labeled mRNA were different from that of non-labeled mRNA, so we compared the agarose-gel electrophoresis pattern of amplified cDNA from BrU-labeled mRNA (nascent mRNA) with that from non-labeled mRNA (steady-state mRNA). The result obtained by agarose-gel electrophoresis experiment of amplified cDNA from BrU-labeled mRNA shows a smear pattern identical to or very similar to that of amplified cDNA from non-labeled mRNA (data not shown). This result suggests that both the reverse-transcription efficiency and the amplification efficiency were almost identical between BrU-labeled RNA and non-labeled RNA. A self–self hybridization experiment was conducted with amplified cDNA labeled with Cy3/Cy5, and we detected gene expression levels with broad signal ranges. The ratios of Cy5 to Cy3 in all genes were in the range of 0.5–2.0 (Fig. 4C), indicating that our DNA microarray analysis is highly reproducible.


Figure 4
View larger version (35K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 4. (A) Eluted RNAs were converted to cDNA, and cDNA was amplified by TS-PCR. (B) Validation with quantitative RT–PCR of the DNA microarray data using each amplification method. Genes expressed at various levels in the non-amplification method, TS-PCR amplification method, and T7 in vitro transcription (IVT) amplification method were confirmed using quantitative RT–PCR. Quantitative RT–PCR was performed using the total RNA of these transcription factors. Gapd was used as the control. The same total RNA was used for quantitative RT–PCR and DNA microarray experiments. The ratio of the total target RNA levels was normalized by using those of the Tbp (Mus musculus TATA box binding protein) gene, and the fold change for each gene in the differentiated ES cells compared to the control ES cells was calculated. Error bars represent standard deviations of log (base 2) of ratio. (C) Amplified cDNA labeled with Cy3/Cy5 underwent self-self hybridization.

 
3.6. DNA microarray analysis of nascent RNA in mouse G2 phase-arrested FT210 cells
We then compared the expression profiles of nascent RNAs of 32 000 transcript products between asynchronous mouse FM3A cells and G2 phase-synchronized FM3A cells. To obtain a highly synchronized population of G2 phase cells, we chose tsFT210 cells, which are temperature-sensitive Cdc2 mutant cells of mouse mammary carcinoma FM3A cells.14Go,15Go tsFT210 cells are grown at 32°C (the permissive temperature), but in the G2 phase, they are arrested at 39°C (Fig. 5A). With this cell line, a highly synchronized population of cells is easily available by just shifting the culture temperature. It has already been shown that such synchronization is merely a result of the inactivation of just one protein (CDC2 protein).15Go We also analyzed steady-state RNA expression profiles of asynchronous cells and G2-synchronized cells. Scatter plots of the expression level are shown in Fig. 5B. Considering the variation of DNA microarray data, we regarded the genes that showed changes exceeding twofold as either up- or down-regulated genes. Figure 5C shows Venn diagram analyses of genes whose expression was up- or down-regulated more than twofold in nascent RNA and steady-state RNA. There was no gene whose nascent RNA was up-regulated, whereas its steady-state RNA was down-regulated, and vice versa.


Figure 5
View larger version (35K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		Figure 5. DNA microarray analyses of nascent RNAs. (A) FT210 cells were synchronized in the G2 phase by incubation at 39 °C for 17 h. Cell cycle profiles of asynchronous cells and G2-arrested cells were analyzed by propidium iodide staining and flow cytometry. (B) Scatter plots comparing the expression profiles of asynchronous cells with those of G2-arrested cells in a nascent RNA profile (left panel) and a steady-state RNA profile (right panel). (C) The Venn diagrams show the overlap of genes whose expression levels changed more than twofold in the nascent RNA profile and the steady-state RNA profile.

 
Among the genes whose nascent RNA levels changed more than twofold, cell cycle progression-related genes are listed in Table 1. Of the 18 genes listed, 14 genes exhibited fold changes in the analysis of nascent RNA but not in that of steady-state RNA.


View this table:
[in this window]
[in a new window]

  		Table 1. The partial lists of cell cycle related genes whose level of nascent RNAs changed more than twofold in G2 arrested cells

 
We found interesting expression profiles in Cdkn1a. Cdkn1a (p21) was up-regulated in both analyses (Table 1). Although this gene was previously reported to function as a signal transducer involved in the transition from the G1 phase to the S phase or from the G2 phase to the M phase,17Go the changes in the expression profiles of the gene at the G2 phase were not analyzed. As up-regulation of the steady-state RNA and nascent RNA of this gene is clearly observed, it follows that this gene may have an important function in the transition from the G2 phase to the M phase.

In addition to Cdkn1a, we have found that Pak1 and Pak2 also exhibited interesting expression profiles. The changes in expression levels of these genes were detected in the analysis of only the nascent RNAs. The PAK family is involved in microtubule and actin cytoskeleton reorganization.18Go PAK1 and PAK2 interact with centrosomal kinase Aurora-A, and activated PAK1 phosphorylates Aurora-A in mitosis.19Go On the other hand, the injection of 58 kDa active Pak1 into Xenopus oocytes inhibits cell cleavage by arresting G2/M progression,20Go suggesting that these products have different role at these phases. These results of DNA microarray analysis indicate that mRNAs of these genes were not only synthesized but also degraded between the G2 phase and the M phase. As it is likely that synthesis and degradation occur in the same period, it seems difficult to identify the changes in expression levels of these gene mRNAs just by analyzing the steady-state levels.

We also found that Rplp2, Rpl27, Rps25, Rpl9, and Rps29, which are categorized as ribosomal proteins, were down-regulated at the G2 phase (Table 2). Interestingly, the changes in the expression of these genes were not detected in steady-state RNA, although they were detected in the nascent RNA analysis (Table 2). Changes in gene expression were not detected in other ribosomal proteins in both analyses. Therefore, it is possible that these five ribosomal proteins have distinct features in their gene expression regulatory mechanisms. A previous study predicted binding sites in transcriptional regulatory regions or promoters of mammalian ribosomal protein genes for transcriptional factors.21Go Those authors observed that genes for these five ribosomal proteins lack TATA-box and A/T-rich sequences similar to the TATA box, whereas two-thirds of the genes for ribosomal proteins have the TATA box or its homolog in their promoter regions. Another study found that a unique sequence (a tandem repeat of the TCTCGCGAGA motif) does exist among many human genes that lack the TATA box.22Go This sequence motif was found at especially high frequencies in the genes for cell cycle regulators as well as in human ribosomal proteins. This sequence motif was also found in the upstream regions of three genes (Rpl9, Rps29, and Rps25) out of the five genes that exhibit changes in gene expression in only the nascent RNA analysis. As translational control is also regarded as an important mechanism for cell cycle regulation, it is possible to obtain new insight into cell cycle regulation in the G2 phase by performing functional analyses of these ribosomal proteins.


View this table:
[in this window]
[in a new window]

  		Table 2. The lists of ribosomal proteins whose level of nascent RNAs changed more than twofold in G2 arrested cells

 
3.7. Applications of comprehensive nascent RNA expression analyses
These examples indicate that our DNA microarray analysis of nascent RNAs has advantages over DNA microarray analysis of steady-state mRNAs. One advantage is that it can provide numerous insights into the regulatory mechanism underlying gene expression. For example, as the changes in expression levels of many transcription factors are relatively small in the process of cellular differentiation, our new method will give us novel insights into the structure of the gene expression regulatory network, since this method can detect slight expression changes. Collecting time-course data on nascent RNAs will also provide important information.

Recently, Kenzelmann et al.23Go reported that nascent RNAs could be analyzed by DNA microarray. In their report, the nascent RNAs were mRNAs labeled using 4-thiol uridine and then purified using affinity beads whose surfaces were coated with a toxic organomercury compound.23Go In our method, precipitation was performed using Dynabeads, which are commonly used in classical immunoprecipitation and do not contain a heavy metal. This is another advantage of our method. Although the DNA microarray used in this study is an oligo-DNA microarray, this method is applicable to other platforms such as cDNA microarray or Affymetrix's GeneChip.

In summary, our study provides very useful techniques for isolating nascent RNAs with high specificity. By comparing the data between steady-state RNAs and nascent RNAs, we can analyze changes in gene expression with much higher sensitivity. We can also determine the kinetics of transcription by combining the DNA microarray study of nascent mRNAs with that of existing mRNAs. The application of such an analysis will provide important information in the field of transcriptome analysis. Microarray data of transcripts whose nascent RNAs labels changed more than twofold in G2 arrested cells are showed in Supplementary data.


   Supplementary Data
Top
 Abstract
 1. Introduction
 2. Materials and Methods
 3. Results and Discussion
 Supplementary Data
Funding
 References
 
Supplementary data are available online at www.dnaresearch.oxfordjournals.org.


   Funding
Top
 Abstract
 1. Introduction
 2. Materials and Methods
 3. Results and Discussion
 Supplementary Data
 Funding
References
 
This work has been supported by the grants supplied by the Ministry of Education, Culture, Sports, Science and Technology (MEXT), the Japan Science and Technology Agency (JST), and Bio Matrix Research Inc.


   Footnotes
 
* To whom correspondence should be addressed. Tel. +81 4-7124-1501 ext. 4408. Fax. +81 4-7122-1360. E-mail: yasufumi{at}rs.noda.tus.ac.jp

Edited by Kazuo Shinozaki


   References
Top
 Abstract
 1. Introduction
 2. Materials and Methods
 3. Results and Discussion
 Supplementary Data
 Funding
 References
 

  1. Gunji W., Kai T., Sameshima E., et al. Global analysis of the expression patterns of transcriptional regulatory factors in formation of embryoid bodies using sensitive oligonucleotide microarray systems. Biochem. Biophys. Res. Commun. (2004) 325:265–275.[CrossRef][Web of Science][Medline]
  2. Brown A. M., Jeltsch J. M., Roberts M., Chambon P. Activation of pS2 gene transcription is a primary response to estrogen in the human breast cancer cell line MCF-7. Proc. Natl. Acad. Sci. USA (1984) 81:6344–6348.[Abstract/Free Full Text]
  3. Fan J., Yang X., Wang W., Wood W. H. III, Becker K. G., Gorospe M. Global analysis of stress-regulated mRNA turnover by using cDNA arrays. Proc. Natl. Acad. Sci. USA (2002) 99:10611–10616.[Abstract/Free Full Text]
  4. Garcia-Martinez J., Aranda A., Perez-Ortin J. E. Genomic run-on evaluates transcription rates for all yeast genes and identifies gene regulatory mechanisms. Mol. Cell (2004) 15:303–313.[CrossRef][Web of Science][Medline]
  5. Cheadle C., Fan J., Cho-Chung Y., et al. Control of gene expression during T cell activation: alternate regulation of mRNA transcription and mRNA stability. BMC Genom. (2005) 6:75.[CrossRef]
  6. Cleary M. D., Meiering C. D., Jan E., Guymon R., Boothroyd J. C. Biosynthetic labeling of RNA with uracil phosphoribosyltransferase allows cell-specific microarray analysis of mRNA synthesis and decay. Nat. Biotechnol. (2005) 23:232–237.[CrossRef][Web of Science][Medline]
  7. Wansink D. G., Schul W., van der Kraan I., van Steensel B., van Driel R., de Jong L. Fluorescent labeling of nascent RNA reveals transcription by RNA polymerase II in domains scattered throughout the nucleus. J. Cell Biol. (1993) 122:283–293.[Abstract/Free Full Text]
  8. Jackson D. A., Iborra F. J., Manders E. M., Cook P. R. Numbers and organization of RNA polymerases, nascent transcripts, and transcription units in HeLa nuclei. Mol. Biol. Cell (1998) 9:1523–1536.[Abstract/Free Full Text]
  9. Waksmundzka M., Debey P. Electric field-mediated BrUTP uptake by mouse oocytes, eggs, and embryos. Mol. Reprod. Dev. (2001) 58:173–179.[CrossRef][Web of Science][Medline]
  10. Hoshino M., Tagawa K., Okuda T., Okazawa H. General transcriptional repression by polyglutamine disease proteins is not directly linked to the presence of inclusion bodies. Biochem. Biophys. Res. Commun. (2004) 313:110–116.[CrossRef][Web of Science][Medline]
  11. Grande M. A., van der Kraan I., de Jong L., van Driel R. Nuclear distribution of transcription factors in relation to sites of transcription and RNA polymerase II. J. Cell Sci. (1997) 110:1781–1791.[Abstract]
  12. Haider S. R., Juan G., Traganos F., Darzynkiewicz Z. Immunoseparation and immunodetection of nucleic acids labeled with halogenated nucleotides. Exp. Cell Res. (1997) 234:498–506.[CrossRef][Web of Science][Medline]
  13. Kageyama S., Nagata M., Aoki F. Isolation of nascent messenger RNA from mouse preimplantation embryos. Biol. Reprod. (2004) 71:1948–1955.[Abstract/Free Full Text]
  14. Mineo C., Murakami Y., Ishimi Y., Hanaoka F., Yamada M. Isolation and analysis of a mammalian temperature-sensitive mutant defective in G2 functions. Exp. Cell Res. (1986) 167:53–62.[CrossRef][Web of Science][Medline]
  15. Th'ng J. P., Wright P. S., Hamaguchi J., et al. The FT210 cell line is a mouse G2 phase mutant with a temperature-sensitive CDC2 gene product. Cell (1990) 63:313–324.[CrossRef][Web of Science][Medline]
  16. Kageyama S., Gunji W., Nakasato M., Murakami Y., Nagata M., Aoki F. Analysis of transcription factor expression during oogenesis and preimplantation development in mice. Zygote (2007) 15:117–128.[CrossRef][Web of Science][Medline]
  17. Niculescu A. B. III, Chen X., Smeets M., Hengst L., Prives C., Reed S. I. Effects of p21Cip1/Waf1 at both the G1/S and the G2/M cell cycle transitions: pRb is a critical determinant in blocking DNA replication and in preventing endoreduplication. Mol. Cell. Biol. (1998) 18:629–643.[Abstract/Free Full Text]
  18. Wittmann T., Bokoch G. M., Waterman-Storer C. M. Regulation of leading edge microtubule and actin dynamics downstream of Rac1. J. Cell Biol. (2003) 161:845–851.[Abstract/Free Full Text]
  19. Zhao Z. S., Lim J. P., Ng Y. W., Lim L., Manser E. The GIT-associated kinase PAK targets to the centrosome and regulates Aurora-A. Mol. Cell (2005) 20:237–249.[CrossRef][Web of Science][Medline]
  20. Rooney R. D., Tuazon P. T., Meek W. E., et al. Cleavage arrest of early frog embryos by the G protein-activated protein kinase PAK I. J. Biol. Chem. (1996) 271:21498–21504.[Abstract/Free Full Text]
  21. Perry R. P. The architecture of mammalian ribosomal protein promoters. BMC Evol. Biol. (2005) 5:15.[CrossRef][Medline]
  22. Wyrwicz L. S., Gaj P., Hoffmann M., Rychlewski L., Ostrowski J. A common cis-element in promoters of protein synthesis and cell cycle genes. Acta Biochim. Pol. (2007) 54:89–98.[Web of Science][Medline]
  23. Kenzelmann M., Maertens S., Hergenhahn M., et al. Microarray analysis of newly synthesized RNA in cells and animals. Proc. Natl. Acad. Sci. USA (2007) 104:6164–6169.[Abstract/Free Full Text]

Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrowOA All Versions of this Article:
15/4/241    most recent
dsn015v2
dsn015v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Request Permissions
Google Scholar
Right arrow Articles by Ohtsu, M.
Right arrow Articles by Murakami, Y.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ohtsu, M.
Right arrow Articles by Murakami, Y.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?