A Microarray-based Approach Identifies ADP Ribosylation Factor-like Protein 2 as a Target of microRNA-16*

microRNAs (miRNAs) are generally thought to negatively regulate the expression of their target genes by mRNA degradation or by translation repression. Here we show an efficient way to identify miRNA target genes by screening alterations in global mRNA levels following changes in miRNA levels. In this study, we used mRNA microarrays to measure global mRNA expression in three cell lines with increased or decreased levels of miR-16 and performed bioinformatics analysis based on multiple target prediction algorithms. For further investigation among the predicted miR-16 target genes, we selected genes that show an expression pattern opposite to that of miR-16. One of the candidate target genes that may interact with miR-16, ADP-ribosylation factor-like protein 2 (ARL2), was further investigated. First, ARL2 was deduced to be an ideal miR-16 target by computational predictions. Second, ARL2 mRNA and protein levels were significantly abolished by treatment with miR-16 precursors, whereas a miR-16 inhibitor increased ARL2 mRNA and protein levels. Third, a luciferase reporter assay confirmed that miR-16 directly recognizes the 3′-untranslated region (3′-UTR) of ARL2. Finally, we showed that miR-16 could regulate proliferation and induce a significant G0/G1 cell cycle arrest, which was due at least in part, to the down-regulation of ARL2. In summary, the present study suggests that integrating global mRNA profiling and bioinformatics tools may provide the basis for further investigation of the potential targets of a given miRNA. These results also illustrate a novel function of miR-16 targeting ARL2 in modulating proliferation and cell cycle progression.

Recently, a new class of RNA regulatory genes known as microRNAs (miRNAs) 5 has introduced a new layer of gene reg-ulation in eukaryotes (1,2). miRNAs are endogenous non-coding RNAs consisting of 19 -24 nucleotides. They play an important role in the negative regulation of gene expression by base pairing to complementary sites on target messenger RNAs (mRNAs), thus causing a block in translation or triggering the degradation of the target mRNAs (1,2). miRNAs are now reported to play fundamental roles in a wide variety of biological processes including developmental timing, apoptosis, proliferation, differentiation, organ development, carcinogenesis, immune response, and energy metabolism (1,2). miR-16 is a miRNA that has been extensively investigated. A study by Calin et al. (3) revealed that miR-15a and miR-16 -1 are down-regulated in the majority of chronic lymphocytic leukemia cases, implying that they may function as tumor suppressors in this disease. Cimmino et al. (4) further showed that miR-15a/miR-16 expression is inversely correlated to B cell lymphoma 2 (BCL2) expression, and BCL2 repression by these miRNAs induces spontaneous apoptosis in a leukemic cell line model. Moreover, it was shown that the miR-16 family negatively regulates cell cycle progression from G0/G1 to S by silencing multiple cell cycle genes including CCND1, CCND3, CCNE1, CDK6, CARD10, and CDC27 (5-7). Furthermore, Jing et al. (8) provided evidence that miR-16 is involved in controlling the turnover of AU-rich element-containing mRNA (ARE-RNA), suggesting that miRNAs may play an important role in ARE-mediated mRNA instability. Taken together, miR-16 is involved in the regulation of basal biological processes including cell cycle arrest, apoptosis, and tumorigenesis.
In the present study, we identified a large number of transcripts regulated by miR-16 directly or indirectly. We accomplished this by overexpressing or knocking down miR-16 in several cell lines and measuring global mRNA levels through microarray platforms. We identified ADP-ribosylation factorlike protein 2 (ARL2) among these candidate targets as a direct target of miR-16. The mechanism through which miR-16 executes its functions was also investigated in this study.

EXPERIMENTAL PROCEDURES
Cells, Reagents, and Antibodies-Human lung adenocarcinoma A549 cells, human embryonic kidney 293A cells, and human breast cancer MCF7 cells were purchased from the China Cell Culture Center (Shanghai, China). 293A cells were cultured in RPMI 1640 medium (Invitrogen) supplemented with 10% fetal bovine serum (FBS) (Invitrogen). A549 and MCF7 cells were maintained in DMEM medium (Invitrogen) supplemented with 10% FBS. Cells were grown at 37°C in a  Table S1. 1 To whom correspondence may be addressed. E-mail: kzen@nju.edu.cn. 2 To whom correspondence may be addressed. E-mail: jfzhang@nju.edu.cn. 3 To whom correspondence may be addressed. E-mail: cyzhang@nju.edu.cn. 4 To whom correspondence may be addressed. E-mail: xichensunny@gmail.
RNA Isolation and Quantitative RT-PCR-Total RNA was extracted from the cultured cells using TRIzol Reagent (Invitrogen) according to the manufacturer's instructions. For quantitative RT-PCR analysis of ARL2 and ␤-actin, 1 g of total RNA was reverse transcribed to cDNA with oligdT and Thermoscript (TaKaRa, Dalian, China). Real-time PCR for ARL2 and ␤-actin was performed on an Applied Biosystems 7300 Sequence Detection System (Applied Biosystems, Foster City) using SYBR green dye (Invitrogen). A 20-l PCR reaction was used and included 1 l of RT product, 1ϫ QuantiTect SYBR green PCR Master Mix, and 0.5 M sense and antisense primers. The reactions were incubated in a 96-well plate at 95°C for 5 min followed by 40 cycles of 95°C for 30 s, 60°C for 30 s, and 72°C for 30 s. All reactions were run in triplicate. After the reactions were run, the threshold cycles (C T ) were determined using fixed threshold settings. The sequences of the sense and antisense primers used for amplification of ARL2 and ␤-actin were as follows: ARL2 (sense): 5Ј-GGGAGGACATCGACA-CCA-3Ј; ARL2 (antisense): 5Ј-AGGACCGCAGGGACTTCT-3Ј; ␤-actin (sense): 5Ј-AGGGAAATCGTGCGTGAC-3Ј; and ␤-actin (antisense): 5Ј-CGCTCATTGCCGATAGTG-3Ј.
Assays to quantify mature miR-16 were carried out using Taqman microRNA probes (Applied Biosystems) according to previous publications (9,10). Briefly, 5 l of total RNA was reverse-transcribed to cDNA using AMV reverse transcriptase (TaKaRa, Dalian, China) and a stem-loop RT primer (Applied Biosystems). Real-time PCR was performed using a TaqMan PCR kit on an Applied Biosystems 7300 Sequence Detection System (Applied Biosystems). All reactions, including notemplate controls, were run in triplicate. After the reactions, the C T values were determined using fixed threshold settings. In the experiments presented here, miRNA expression in cells was normalized to U6 snRNA, similar to many other reports (11). The relative amount of miR-16 to internal control U6 was calculated with the equation 2 Ϫ⌬CT , in which ⌬C T ϭ C T miR-16 Ϫ C T U6 .
Microarray Procedures-The commercially available 22K Human Genome Array was purchased from CapitalBio Corporation (Beijing, China). Labeling, hybridization, washing, and scanning were performed according to the standard operating procedure of CapitalBio Corporation. Briefly, total RNA was used to synthesize cDNA in an in vitro transcription reaction. cDNA was fluorescently labeled by Cy5 or Cy3-CPT with Klenow enzyme. After hybridization, the microarray was washed with two consecutive washing solutions (0.2% SDS and 2ϫ SSC at 42°C for 5 min, and 0.2% SSC for 5 min at room temperature). Subsequently, these arrays were scanned with a LuxScan 10KA confocal laser scanner (CapitalBio Corp.), and the obtained images were analyzed using LuxScan Version 3.0 (CapitalBio Corp.), which employed the LOWESS normalization method.
miR-16 Target Prediction-The analysis of miRNA predicted targets was determined with the algorithms from TargetScan (12), PicTar (13), and miRanda (14). Western Blotting-ARL2 protein level was quantified by Western blot analysis of whole cell extracts using antibodies against ARL2. These samples were normalized by blotting with an antibody against GAPDH.
Plasmid Construction and Luciferase Assay-The entire human ARL2 3Ј-untranslated region (3Ј-UTR) was amplified by PCR using human genomic DNA as a template. The PCR products were inserted into the p-MIR-report plasmid (Ambion). Efficient insertion was confirmed by sequencing. For luciferase reporter assays, cells were cultured in 6-well plates, and each well was transfected with 2 g of firefly luciferase reporter plasmid, 2 g of ␤-galactosidase expression vector (Ambion), and equal amounts of scrambled ncRNA, pre-miR-16, or anti-miR-16 using Lipofectamine 2000 (Invitrogen). The ␤-galactosidase vector was used as a transfection control. At 24 h post-transfection, cells were assayed using luciferase assay kits (Promega, Madison, WI). Data depicted are representative of three independent experiments performed on different days.
siRNA was transfected into A549 cells using Lipofectamine 2000 (Invitrogen) according to the manufacturer's instructions. Total RNA was isolated at 24 h post-transfection. The expression level of ARL2 mRNA was assessed by quantitative RT-PCR. The sequence with the best interfering effect (named siRNA-3) was selected and used in further studies.
Cell Viability Assay-A549 cells were plated at 2.5 ϫ 10 3 cells per well in 96-well plates and incubated overnight in DMEM medium supplemented with 10% FBS. After transfection, 20 l of 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyl tetrazolium bromide (MTT) (5 mg/ml) was added into a corresponding test well and incubated for 4 h. The supernatant was then discarded, and 200 l of DMSO was added to each well to dissolve the precipitate. Absorbance (A) was measured at a wavelength of 570 nm.
Cell Migration Assay-The migration ability of A549 cells was tested in a Transwell Boyden Chamber (6.5 mm, Costar, Cambridge, MA). The polycarbonate membranes (8-m pore size) on the bottom of the upper compartment of the Transwells were coated with a 0.1% gelatin matrix. Cells treated with ncRNA, pre-miR-16, or siRNA for 6 h were suspended in serum-free DMEM culture medium at a concentration of 4 ϫ 10 5 cells/ml and then added to the upper chamber (4 ϫ 10 4 cells/well). Simultaneously, 0.5 ml of DMEM with 10% FBS was added to the lower compartment, and the Transwell-containing plates were incubated for 8 h in a 5% CO 2 atmosphere saturated with H 2 O. At the end of the incubation, cells that had entered the lower surface of the filter membrane (migrant cells) were fixed with 90% ethanol for 15 min at room temperature, washed three times with distilled water, and stained with 0.1% crystal violet in 0.1 M borate and 2% ethanol for 15 min at room temperature. Cells remaining on the upper surface of the filter membrane (nonmigrant) were scraped off gently with a cotton swab. Images of migrant cells were captured by a photomicroscope (BX51, Olympus, Japan). Cell migration was quantified by blind counting of the migrated cells on the lower surface of the membrane with five fields per chamber.
Cell Cycle Assay-Cells were harvested, washed once with PBS, and fixed in 70% ethanol overnight. Staining for DNA content was performed with 50 mg/ml propidium iodide and 1 mg/ml RNase A for 30 min. Analysis was performed on a fluorescence-activated cell-sorting (FACS) flow cytometer (BD Biosciences, San Jose, CA) with Cell Quest Pro software. Cellcycle modeling was performed with Modfit 3.0 software (Verity Software House, Topsham, ME).
Statistical Analysis-All photo-images of Western blotting and semi-quantitative RT-PCR were representative of at least three independent experiments. Quantitative RT-PCR, luciferase reporter assay, and cell migration assays were performed in triplicate, and each experiment was repeated three to five times. Data shown are presented as mean Ϯ S.D. of at least three independent experiments. Statistical significance was considered at p Ͻ 0.05 using the Student's t test.

Identification of Transcripts Regulated by miR-16 by
Microarray Analysis-We first screened the expression levels of miR-16 in A549, 293A, and MCF7 cells. Our results showed that miR-16 can be readily detected in all cells (Fig. 1, A and B), consistent with the fact that miR-16 is ubiquitously expressed in many mammalian tissues and cells.
miRNAs are generally thought to negatively regulate the expression of their targets by mRNA degradation or by translational repression (1,2). Therefore, if a given miRNA can mediate the degradation of targeted mRNAs, this miRNA and its targets should have expression patterns in the opposite directions. In light of this mode of miRNA action, we postulated that candidate miRNAs affecting transcript levels could be deduced from mRNA microarray expression profiles. We therefore transfected A549, 293A, and MCF7 cells with equal doses of scrambled ncRNA, pre-miR-16, or anti-miR-16. We then surveyed potential genes that were inversely expressed relative to miR-16 using microarray analysis. As shown in Fig. 1C, the expression of miR-16 was abolished by the introduction of anti-miR-16, whereas pre-miR-16 significantly increased miR-16 levels in A549, 293A, and MCF7 cells. The mRNA microarray profiles clearly showed different mRNA expression patterns among pre-miR-16-and anti-miR-16-transfected cells compared with ncRNA-transfected cells. In total, 115, 389, and 188 genes were down-regulated (mean fold change Ͻ0.66) in miR-16-overexpressed A549, 293A, and MCF7 cells, respectively, while 45, 218, and 34 genes were up-regulated (mean fold change Ͼ 1.5) in miR-16-down-regulated A549, 293A, and MCF7 cells, respectively (Fig. 1D). To reduce the false positives and obtain a more accurate assessment of the genuine miR-16 targets, only the mRNAs that showed consistent expression patterns in all three cell lines were considered as candidate miR-16 targets. A set of 26 down-regulated and 6 up-regulated genes in these cells were identified as candidate miR-16 targets ( Fig. 1D and Tables 1 and 2).
Enrichment of miR-16 Target Genes among Differentially Expressed Genes-We next calculated whether the differentially regulated genes were predicted miR-16 targets using three widely used programs: TargetScan, PicTar, and MiRanda (12)(13)(14). Because individual computer-aided algorithms can generate a high number of false positives, the combination of these three approaches provided a more accurate assessment of the genuine miRNA targets than a single approach would. Therefore, these algorithms were used in combination to provide the highest probability of target identification, and only the genes predicted as miR-16 targets by at least two of the above-mentioned algorithms were considered positive. Among the list of the 26 down-regulated genes obtained from miR-16 overexpression, 13 were predicted as miR-16 targets (Table 3). From the six up-regulated genes obtained from miR-16 knockdown, 2 genes were identified as miR-16 targets (Table 3). In contrast, out of the 423 unchanged genes in the three cell lines (0.95 Ͻ mean fold change Ͻ 1.05), only 15 were identified as predicted miR-16 targets (Table 3). Moreover, seven genes showed the same direction of change with miR-16 (up-regulation upon overexpression of miR-16, mean fold change Ͼ 1.5; or downregulation upon knockdown of miR-16, mean fold change Ͻ 0.66); however, none were identified as predicted targets of miR-16 (Table 3). These results indicate that the predicted miR-16 targets are more frequently found among genes with expression patterns opposite to that of miR-16, further indicating that these relationships may be functional miRNA-target combinations.
The ARL2 3Ј-UTR Contains Four Putative Target Sites for miR-16-As noted, ARL2 was the most significantly down-regulated gene in all three cell lines when miR-16 was overexpressed (Table 1). Intriguingly, ARL2 was identified as a candidate target of miR-16 by all three programs (supplemental Table S1). The predicted interaction between miR-16 and its target sites in the ARL2 3Ј-UTR is illustrated in Fig. 2A. As shown in this figure, there are four potential target sites in the 3Ј-UTR of ARL2 for miR-16. The minimum free energy values of these hybrids are Ϫ14.2, Ϫ15.8, Ϫ15.5, and Ϫ14.8 kcal/mol, determined by RNA hybrid analysis (15), which are well within the range of authentic miRNA-target pairs ( Fig. 2A). Moreover, perfect base pairing between the "seeds" (the core sequence that encompasses the first 2-8 bases of the mature miRNA) and cognate targets ( Fig. 2A) was noted, and the miR-16 binding sequences at the ARL2 3Ј-UTR are highly conserved across species (Fig. 2B). The common criteria to decide whether a transcript is a target for a miRNA is base pairing between the "seed" and target, strict interspecies conservation of the miRNA binding sequence, and an inverse correlation between the miRNA expression levels and its target levels. Because miR-16 targeting ARL2 fulfills these criteria, we next focused on studying the role of miR-16 in regulating ARL2 expression.   Down-regulation of ARL2 Expression by miR-16-Microarray analysis showed that miR-16 was able to regulate the steadystate level of ARL2 mRNA (Table 1). We further validated this result by performing a quantitative RT-PCR assay with the same RNA preparations as discussed previously. A549, 293A, and MCF7 cells were transfected with ncRNA, pre-miR-16, or anti-miR-16 and analyzed for the expression of ARL2 mRNA by quantitative RT-PCR at 24 h post-transfection. All cells trans-fected with pre-miR-16 showed reduced levels of ARL2 mRNA relative to the cells transfected with ncRNA. In contrast, anti-miR-16 significantly increased the expression levels of ARL2 mRNA in these three cell lines (Fig. 3A).
To determine whether the overexpression or knockdown of miR-16 had an impact on ARL2 protein expression, we repeated the above-mentioned experiments and determined the expression of ARL2 protein by Western blotting at 24 h FIGURE 2. Identification of conserved miR-16 binding sites within the ARL2 3 UTR. A, schematic description of the hypothesized duplexes formed by interactions between the ARL2 3Ј-UTR binding sites and miR-16. Each binding site is numbered, and the predicted structure of each base-paired hybrid is diagrammed. Paired bases are indicated by a black oval, and G:U pairs are indicated by two dots. The predicted free energy of each hybrid is indicated. B, sequence alignment of the putative miR-16 binding sites across species. The seed complementary sites are marked in red, and all nucleotides in these regions are completely conserved in several species including human, mouse, and rat. post-transfection. As shown in Fig. 3B, the expression of ARL2 protein was significantly abolished by the introduction of pre-miR-16, whereas cells transfected with scrambled ncRNA maintained a considerable amount of ARL2 protein. In contrast, anti-miR-16 significantly increased ARL2 protein in A549 cells. Moreover, the ARL2 protein expression was inversely correlated to miR-16 in both 293A and MCF7 cells (data not shown). These results demonstrate that miR-16 regulates the expression of ARL2 at both the transcript and protein levels. ARL2 Is a Direct Target of miR-16-To determine whether the negative regulatory effects of miR-16 on ARL2 expression were mediated through binding to the presumed complementary sites at the 3Ј-UTR of ARL2, we fused the entire ARL2 3Ј-UTR into a downstream position of the firefly luciferase reporter plasmid. The resulting plasmid was introduced into A549 cells combined with a transfection control plasmid (␤-gal) and pre-miR-16, anti-miR-16, or scrambled ncRNA. As expected, overexpression of miR-16 resulted in a significant decrease in the luciferase reporter activity (normalized against ␤-gal activity) compared with the treatment with scrambled ncRNA; however, inhibition of miR-16 resulted in a significant increase in reporter activity (Fig. 3C). The observed alterations in luciferase activity were specific, as co-transfection of A549 cells with the parental luciferase plasmid (without the ARL2 3Ј-UTR) and pre-miR-16 or anti-miR-16 did not affect luciferase reporter activity (Fig. 3C). Furthermore, we introduced a point mutation into the corresponding seed complementary sites in the ARL2 3Ј-UTR to eliminate the predicted binding by miR-16 (Fig. 3D). As shown in Fig. 3C, mutations in the seed complementary sites almost fully rescued the repression of reporter activity by pre-miR-16. Also, the mutant luciferase reporter was unaffected by miR-16 inhibition in A549 cells. This suggests that these binding sites strongly contribute to the miRNA:mRNA interaction that mediates the post-transcriptional inhibition of ARL2 expression. Moreover, similar results were observed in experiments carried out in 293A and MCF7 cells (data not shown). In conclusion, our results demonstrate that miR-16 directly recognizes the 3Ј-UTR of the ARL2 transcript to down-regulate its expression.
The Role of miR-16 and ARL2 in Cell Cycle and Proliferation-To investigate the cellular phenotypes triggered by miR-16 targeting ARL2, A549 cells were transfected with miR-16 precursors or siRNA against ARL2 and analyzed for the changes of migration, proliferation, and cell cycle arrest. Cells transfected with ncRNA or control siRNA served as controls. Efficient interference of ARL2 expression is shown in Fig. 4, A and B.
We assessed the role of miR-16 on cell migration by Transwell assay. As shown in Fig. 4C, the migration rate of A549 cells transfected with pre-miR-16 was significantly decreased compared with control cells. However, no differences were observed in migration rates between the cells transfected with negative control siRNA or siRNA against ARL2. These results suggest that miR-16 could modulate cell migration by downregulating genes other than ARL2.
The proliferation rates of A549 cells with increased or decreased expression of miR-16 were determined via MTT Results are representative data from three independent experiments. Pictures of the Western blot assay were analyzed using Bandscan software, and a statistical analysis is present below. C, direct recognition of the ARL2 3Ј-UTR by miR-16. Firefly luciferase reporters containing either wt or mut ARL2 3Ј-UTR were co-transfected into A549 cells with scrambled ncRNA, pre-miR-16, or anti-miR-16. The parental luciferase plasmid was also transfected as a control. At 24 h post-transfection, cells were assayed using luciferase assay kits. For comparison, the luciferase activity in ncRNA-transfected cells was arbitrarily set at 1. The y axis shows arbitrary units representing relative luciferase activity. The results are presented as the mean Ϯ S.D. of three independent experiments (*, p Ͻ 0.05; ***, p Ͻ 0.001). D, diagram of the luciferase reporter plasmid carrying the firefly luciferase coding sequence attached to the human ARL2 3Ј-UTR with either wild-type (wt) or mutant (mut) miR-16 binding sites (the sequence that interacts with the 2-4 bases of miR-16 were mutated).
assay. Compared with ncRNA-transfected cells, cells transfected with anti-miR-16 proliferated at a significantly higher rate. In contrast, overexpression of miR-16 by transfection with pre-miR-16 resulted in significantly reduced proliferation (Fig.  4D). Furthermore, significant difference was observed in proliferation rates between the cells transfected with control siRNA and siRNA against ARL2. The results suggest that miR-16 might inhibit cell proliferation by silencing ARL2.
Finally, we investigated cell cycle distribution in cells with enhanced miR-16 or silenced ARL2 by flow cytometry analysis. Compared with cells transfected with ncRNA, cells transfected with pre-miR-16 triggered an accumulation of cells in the G0/G1 stage, whereas the numbers of cells in the S and G2/M phases decreased (Fig. 4E). The percentage of ncRNA-treated cells in the G0/G1 phase was 61.12 Ϯ 0.22%, while transfection with pre-miR-16 resulted in 75.33 Ϯ 0.85% of cells in the G0/G1 phase. The S phase fractions in control and pre-miR-16 groups were 28.11 Ϯ 1/19% and 17.90 Ϯ 1.27%, respectively, and the G2/M phase fractions in control and pre-miR-16 groups were 10.77 Ϯ 1.01% and 6.78 Ϯ 0.81%, respectively. Transfection with siRNA against ARL2 partially yielded the phenotype generated by overexpression of miR-16, but the effects were somewhat different, featured by a lower G0/G1 cell accumulation but a higher S and G2/M cell accumulation (Fig. 4E). In A549 cells, 61.26 Ϯ 0.30% of the cell population went into arrest in the G0/G1 phase when ARL2 was interfered with siRNA, whereas 67.60 Ϯ 0.17% of the control cells transfected with ncRNA were in this phase of the cell cycle. The S phase fractions in control and siRNA-treated groups were 27.34 Ϯ 0.92% and 24.18 Ϯ 0.67%, respectively, and the G2/M phase fractions in controland siRNA-treated groups were 11.40 Ϯ 0.73% and 8.22 Ϯ 0.54%, respectively. Comparable results were also observed in 293A and MCF7 cells (data not shown). Taken together, our data suggests that miR-16 might negatively regulate cell cycle progression from the G0/G1 phase to the S phase by silencing ARL2.

DISCUSSION
Although the number of known miRNAs is continuously increasing, information regarding their precise cellular function remains limited. One of the main challenges in understanding the functions of miRNAs is to identify the genuine target genes of miRNAs. However, lack of reliable and specific methods for biological target validation hampers the full understanding of the mechanisms by which miRNAs execute their functions. Only a few miRNAs have thus far been assigned tar-get mRNAs, and the conventional methodologies are still labor intensive. Therefore, novel approaches for target identification are required to overcome this limitation.
Animal miRNAs were originally believed to block translational processes without affecting transcript levels (1, 2); however, recent studies have suggested that they can also induce mRNA degradation or destabilization even in the absence of extensive base pairing to their targets (16 -21). This was first clearly demonstrated by the observation that overexpression of a specific miRNA could reduce the expression of hundreds of mRNAs (16). Subsequently, the action of miRNAs on mRNA levels has been widely confirmed in other studies (17)(18)(19)(20). A recent publication reported the impact of miRNAs on global mRNA and protein expression and showed that the regulation of protein-coding genes by miRNAs is quite similar at both the transcript and protein levels (21). Inspired by this mode of miRNA action, we postulated that the high-throughput mRNA microarray assays capable of detecting such effects at the mRNA level would provide a promising way to assist in miRNA target identification. Indeed, by globally screening the expression patterns of transcripts in several cell lines using microarray analysis, we found that a large number of transcripts change after miRNA transfection ( Fig. 1 and Tables 1 and 2). To obtain a more accurate assessment of the genuine miR-16 targets, we combined bioinformatics programs to select candidate miR-16 targets from the differentially regulated genes. In contrast to unchanged gene expression and genes showing a similar direction of expression with miR-16, putative miRNA targets revealed significant enrichment in genes having expression patterns opposite to that of miR-16 (Table 3). These results suggest that many of these genes are genuine targets of miR-16 rather than genes secondarily regulated as part of the miRNA-mediated module controlled by the direct miRNA targets. Therefore, our strategies combining global gene expression analysis with bioinformatics prediction may provide additional confirmatory evidence to aid in miRNA target identification.
Among a panel of candidate targets that may interact with miR-16, ARL2, a 21-kDa GTPase belonging to the ADP-ribosylation factor (ARF) family (22,23), was further experimentally validated as a miR-16 target. The reasons we focused on miR-16 targeting ARL2 include the following: (a) ARL2 was down-regulated most significantly in all three cell lines when miR-16 was overexpressed and (b) ARL2 was deduced to be an ideal miR-16 target by computational predictions. Indeed, through several biological approaches, we demonstrated that miR-16 regulates the expression of ARL2 at both the transcript and protein levels. The results demonstrate that our approach is a powerful tool for miRNA target identification.
A number of genetic studies suggest that ARL2 is involved in tubulin folding and plays a role in microtubule dynamics (24 -30). In the present study, we showed that ARL2 had the potential to influence proliferation and cell cycle distribution in A549 cells. Moreover, we investigated whether the cellular phenotypes including migration, proliferation, and cell cycle arrest were regulated by miR-16 targeting ARL2. We showed that miR-16 could modulate cell migration by down-regulating genes other than ARL2, but that it negatively regulated proliferation and cell cycle progression from the G0/G1 phase to the S phase by silencing ARL2. However, siRNA against ARL2 only partially phenocopied the cellular phenotype of miR-16 on controlling cell cycle progression. One explanation for the more profound cell cycle arrest at the G0/G1 phase elicited by miR-16 rather than siRNA against ARL2 is that multiple cell cycle genes coordinately modulating cell cycle progression are targeted by miR-16. Indeed, several other cell cycle genes are regulated by miR-16, including CCND1, CCND3, CCNE1, CDK6, CARD10, and CDC27 (5-7). An emerging common theme is that multiple targets regulated by a single miRNA can act in concert, rather than individually, to regulate the same biological process. Coordinated regulation of many targets by a single miRNA may allow for a prompt cellular response to progress the cell cycle and also for rapid reversal of the miRNA-induced cell cycle regulation upon changes in miRNA synthesis, stability or localization. Taken together, our results combined with others demonstrate that the miR-16 family induces cell cycle arrest by silencing multiple downstream effectors simultaneously rather than an individual target.
In this study, we uncovered multiple targets of miR-16, including ARL2, by combining mRNA microarray profiles with bioinformatics analysis. We further experimentally validated ARL2 as a direct target of miR-16 and demonstrated that this targeting plays a critical role in regulating proliferation and cell cycle phenotypes.