- Research article
- Open Access
A noncoding RNA gene on chromosome 10p15.3 may function upstream of hTERT
BMC Molecular Biology volume 10, Article number: 5 (2009)
We attempted to clone candidate genes on 10p14–15 which may regulate hTERT expression, through exon trapping using 3 BAC clones covering the region. After obtaining 20 exons, we examined the function of RGM249 (RGM: RNA gene for miRNAs) we cloned from primary cultured human hepatocytes and hepatoma cell lines. We confirmed approximately 20 bp products digested by Dicer, and investigated the function of this cloned gene and its involvement in hTERT expression by transfecting the hepatoma cell lines with full-length dsRNA, gene-specific designed siRNA, and shRNA-generating plasmid.
RGM249 showed cancer-dominant intense expression similar to hTERT in cancer cell lines, whereas very weak expression was evident in human primary hepatocytes without telomerase activity. This gene was predicted to be a noncoding precursor RNA gene. Interestingly, RGM249 dsRNA, siRNA, and shRNA inhibited more than 80% of hTERT mRNA expression. In contrast, primary cultured cells overexpressing the gene showed no significant change in hTERT mRNA expression; the overexpression of the gene strongly suppressed hTERT mRNA in poorly differentiated cells.
These findings indicate that RGM249 might be a microRNA precursor gene involved in the differentiation and function upstream of hTERT.
Telomerase, which adds repeated telomere sequences to chromosome ends, is a ribonucleoprotein enzyme that includes an endogeneous RNA (hTR) which acts as a template and a reverse transcriptase (hTERT), and are crucial molecules [1–4]. Telomerase activity has been detected in immortalized and tumor cells and its strong expression represents an important difference between normal and cancer cells . Although it is thought that telomere maintenance is predominantly due to telomerase activity [6, 7], which is mainly regulated by hTERT expression in the presence of hTR [8, 9], some gaps exist in the steps between immortalization, carcinogenesis, and cellular senescence . Therefore, these unknown steps must first be fully identified before translational applications are possible [11–13]. As several reports have indicated that human chromosome 10p carries a set of genes involved in regulating telomerase activity [14–16], we performed microcell-mediated chromosome transfer, radiated microcell fusion, and bacterial artificial chromosome (BAC) transfer assuming loss of heterozygosity (LOH), and examined telomerase expression in transfectants to determine whether hTR- and hTERT-regulating genes are located on 10p15 . After we determined the telomerase regulatory regions on chromosome 10p15.1 using DNA markers, we confirmed that the chromosomal fragments containing the regions induce senescence in the transfectants [18, 19]. We then identified that radiated chromosomal fragments containing D10S1728 (a genomic region on chromosome 10p15) induced senescence in transfectants, and performed the sequencing of 3 BAC clones containing D10S1728 [20–22] (data not shown), exon trapping using the BAC clones [23–25], Northern blotting, and RT-PCR before cloning several full-length genes. One of the genes we cloned was named RGM249 (RNA gene for miRNAs which was 249 bp in length), and showed stronger gene expression in cancerous lesions than in normal liver.
We used the RNA interference (RNAi) method, which is a sequence-specific post-transcriptional gene silencing mechanism triggered by double-stranded RNA (dsRNA), to cause degradation of mRNAs homologous in sequence to dsRNA and to effectively inhibit gene-specific expression . Because dsRNA transfection is a common strategy for suppressing genes of interest, the small interfering RNA (siRNA) technique is an excellent gene-specific method for determining the function of the targeted gene by inhibiting its expression [27–30]. After focusing on RGM249 (RGM: RNA gene for miRNAs; 249 for 249 bp) among the genes cloned on 10p14–15 for further analyses, we examined the telomerase regulatory function of RGM249 using the overexpressing cassette and the gene-specific siRNA.
In this report, we present a summary of our cloning research and introduce functional analyses regarding the identified gene regulating hTERT mRNA.
Results and discussion
Screening of BAC library and transfection into HMc-Li7
We investigated the 3 BAC clones using a sequence-tagged-site marker (D10S1728) by PCR, based on previous reports [18, 19]. From these studies, we determined that the 3 BAC clones (M-6, H-11, and J-21) contained D10S1728.
Whole sequencing of H-11 and exon trapping
Approximately 90%~95% of the BAC clone H-11 sequence and supercontig were determined (studies were performed by the University of Colorado Cancer Center Sequencing Core). D10S1728 was regarded as the sequence on 10p15.3. We reconfirmed that the sequence and supercontig were derived from the original BAC clones by Southern blotting and PCR, and 20 exons were thereafter identified (fig 1-(a)). We performed a computer analysis for each contig, and the endosequences of the 3 BAC clones and the exon sequences we trapped showed that they closely overlapped each other (fig 1-(b)).
Northern blot analysis and one-step real-time RT-PCR
Using a probe against each exon sequence (47–376 bp), Northern blotting was performed, which demonstrated that several exons showed a very faint signal (data not shown). In particular, a very faint single hybridization signal in human Multiple Tissue Northern (MTN) blot (Clontech) migrating at about 250 bp was determined for RGM249. One-step real-time RT-PCR showed a very low expression rate, that is, less than 20,000 copies per 50 ng of total RNA extracted from human organs (fig 2-(b)). In normal liver, approximately 60 copies per 50 ng of total RNA of transcripts were expressed. In contrast, the study on human primary hepatocyte cultures, including nuclear and cytoplasmic compartments in 4 hepatoma cell lines (HepG2, Huh7, HMc-Li7 and Alexander), indicated that hepatoma cells expressed RGM249 at much higher levels than normal human hepatocytes (fig 2-(c)). In 39 hepatomas, the expression rate was higher (mean ± SEM: 609.2 ± 261.5 copy number per 100 ng/μl) than in the 35 adjacent noncancerous lesions (138.9 ± 46.5) (P = 0.007 by t-test and P = 0.040 by paired t-test) (fig 2-(d)).
5' & 3' RACE and cDNA library screening
5' and 3' RACE was performed using liver c DNA libraries (Clontech). Screening of the cDNA libraries yielded a clone with a full sequence length of 249 bp. A comparison of the cDNA sequence with the genomic sequence from H-11 revealed that the gene of interest was extended across 9 kb. The cDNA named RGM249 consisted of 2 exons and the RNA was not poly-adenylated (fig 1-(b)). The sequence with about 40% of the GC ratio is shown in fig 3-(a), indicating that this gene might be a noncoding RNA (ncRNA) gene, as shown in fig 3-(b). PCR-Single Strand Conformation Polymorphism (PCR-SSCP) and sequencing showed that the gene had no mutations other than a one-base addition at the last site of exon 1 in 10 normal liver tissue samples (10%; 1/10) and no mutations in 4 tumor cell lines (0%; 0/4). Although the cDNA showed no homology with either additional human loci except for the endogeneous locus, or mouse genome, the miRNA predicted to be digested from RGM249 by Dicer matched at least 3 miRNAs in the miRNA database (miRBase; mature miRNAs or stem-loop sequences) http://microrna.sanger.ac.uk/sequences/search.shtml. Although their functions remain unknown because we have not fully clarified whether RGM249 is actually processed by Dicer in vivo, given that it can be processed by Dicer in vitro. On the other hand, a motif search predicted no motifs in the RGM249 gene. Because this gene generated no protein, the secondary structure of RGM249 was drawn based on the structure provided at http://www.tbi.univie.ac.at/~ivo/RNA/ (fig 3-(b)). Moreover, because the full-length RNA was digested by Dicer enzymes, 20–23 bp sharp bands could be observed, indicating that nucleotides generated from RGM249 RNA may be miRNAs (fig 3-(c)).
GeneBank accession number
Sequence data have been deposited with the EMBL/GenBank Data Libraries under accession no. EF433558.
The effects of downregulating RGM249
In many hTERT- or hTR-regulating molecules reported previously, only two factors, hTERT and ribosomal protein L22 (RPL22), were transcriptionally suppressed in the hepatoma-derived cell lines (HLF) transfected with full-length RGM249 dsRNA (fig 4-(a)). To reconfirm the findings in dsRNA method, the siRNA method was performed. Results showed that their expressions were suppressed in HLF cells transfected with Stealth siRNA ① (fig 4-(b),(c)).
RGM249 siRNA suppressed more than 90% of hTERT mRNA, approximately 50% of RPL22 mRNA, and RGM249 transcripts (fig 4-(c)) unlike hTR expression. However, RPL22 expression in cells transfected with Stealth siRNA ① was weaker (50%) than in those transfected with full-length dsRNA (fig 4-(c)). siRNA ① + ② had less efficient knockdown due to the interaction between both siRNAs. These findings are consistent with our results obtained from dsRNA induction, an experiment performed under conditions of almost negligible cytokine storm. A similar result was obtained about 48 hours after transfection.
Subsequently, we examined the inhibitory effect of shRNA transfection. Compared with LacZ shRNA, we found that RGM249 shRNA and mt-1RGM249 shRNA significantly downregulated hTERT mRNA (P = 0.002 and P = 0.034, respectively, using the Mann-Whitney test; fig 4-(d)). Otherwise, no significant differences were observed, except for the effect on RGM249 mRNA. Similar to transfectants that received RGM249 siRNA, telomerase activity in transfectants with shRNA was suppressed in the same manner as mRNA expression, compared with that in transfectants with LacZ shRNA (fig 4-(e)). Telomere length was not measured due to the small number of cells available and the relatively short observation period. No immunoreactive effects on tumor cell proliferation resulting from cytokine storm were recognized in the GE profile by microarray analysis, which resulted in IFN-β being downregulated twofold.
The effects of overexpressing RGM249 mRNA
In the transfection of HLF or A172 cells showing RGM249 overexpression, as distinct from the transfection of other cell lines such as senescing or well-differentiated cells, hTERT mRNA was significantly downregulated. In particular, the growth of A172 transfectants was strongly suppressed and they were induced to undergo senescence. On the other hand, hTERT mRNA was upregulated in KMST cells, which was established by exposure to chemical oncogens (Additional file 1). However, the expression of hTR, β-actin, and RPL22 mRNA showed no significant differences between mock transfectants and RGM249 transfectants in KMST cells, similar to the transcriptional level of c-Myc, undifferentiated markers (Rex-1 and GATA4), and an initiation factor (eIF6). Additionally, the longevity and proliferation rate of transfectants were almost the same as those of the parental cells or mock transfectants. Interestingly, we examined the effect of RGM249 shRNA wrapped by cationized carriers on tumorigenesis in nude mice using local or venous injection, which resulted in more than 60% of the suppressive effect in vivo (manuscript in preparation). FACS analysis showed that RGM249 mRNA expression introduced more than half of RGM249 pRNAT/U6.1/neo/GFP-transfected populations into the preG1 phase (fig 5-(a)), senescence, and apoptosis (fig 5-(b)), compared with pRNAT/U6.1/neo/GFP-transfectants. These phenotypic changes were similar to the effects of RGM249 shRNA, compared with those of LacZ shRNA
LOH in human hepatocellular carcinoma (HCC) is rare, except at 10q26.11-qter, which is observed in 25% of HCC cases (32). However, relatively frequent alterations and deletions of chromosome 10p have been identified in human tumors, including glioblastoma multiforme (GBM) and malignant melanoma [33, 34]. In GBM, 3 common distinct regions of chromosome 10 deletion have been frequently identified; one of them is the distal half of the short arm [15, 35], indicating that several distinct tumor suppressor genes may be present on chromosome 10 and may contribute to the progression of malignant tumors. Structural alterations have only been observed in malignant melanoma (10p13-pter), gastric cancer (10p15), and GBM (10p14-pter). Unexpectedly, RGM249 with low GC ratio on 10p15.3 is potently miRNAs-generating gene, as well as cancer-related gene.
It is well known that a miRNA, generated from a miRNA precursor gene, can act as a new regulator which induces cancer cell lines to adopt a suppressive phenotype, except for the tumor suppressor gene. Because the phenotype we observed in a series of transfections with a chromosome (fragment) or BAC was suppressed, this miRNA might be a Suppressor-miR or it might have demonstrated a miRNA-directed inhibition. Actually, RGM249 repressed hTERT mRNA in A172 cells derived from GBM, moderately differentiated HCC (HMc-Li7), and poorly differentiated HLF. Beyond 60S RPL22, which we referred to in this study, several ribosomal protein mRNAs changed their expression level after the introduction of siRNA or shRNA into cell lines or tumors inoculated in nude mice (data not shown). In a miRNA-mediated mRNA repression study, RPL22, of which Greider et al. previously referred to as one of the binding protein to hTR that is essentially required for telomere maintenance and telomerase activity in normal tissues and disease, showed no significant alteration in expression level, whereas it showed an obvious change in siRNA-directed RNAi with reproducibility (2.08-fold downregulation in microarray analysis). We do not have any definite evidence that RPL22 functions in association with hTERT expression during carcinogenesis or in hypoxic tissues. However, while endogenous hTR expression is regulated at the transcriptional level, hTERT is regulated by alternative splicing under hypoxic conditions . It may be that RPL22 influences hTERT expression through hTR expression. RGM249 dsRNA repressed RPL22. However, RGM249 siRNA and shRNA did not significantly repress RPL22 because of the high homology shared by 2 parts of the RGM249 sequence which are different from the functional sequence, which was used in siRNA or shRNA and corresponds to one part of the RPL22 sequence. Eukaryotic initiation factor 6 (eIF6) has long been known to bind to free 60S ribosomal subunits, preventing the assembly of translationally competent 80S ribosome particles . Similar to eIF4E bound to Argonaute2 (Ago2), which regulates mRNA expression by cleaving its mRNA targets and competing with eIF6 for the mRNA cap , eIF6 may mediate siRNA-directed RNA interference rather than miRNA-directed translational repression. This study indicates that the sequence focused for siRNA or shRNA has a key function in this gene, and that RGM249 may functionally have a close link to gene-mediated processes at the transcriptional level, such as RNA biogenesis or gene regulation, and a link to ribonucleoproteins such as telomerase or RPL22. As a result, RNA molecules may induce suppressive or apoptotic events in cancer cells through an RGM249-related signaling pathway involved in a DNA-, RNA-, or protein-binding function (some of the significantly changed genes are shown in fig 6; there were 306 upregulated genes and 144 downregulated genes in the in vitro trial)
miRNA is physiologically generated by both Drosha and Dicer. We can not digest RGM249 mRNA by Drosha because Drosha protein is now unavailable commercially. RGM249 RNA processed by DICER, results in at least 3 small fragments of a product (17–23 bp) related to RNA species (fig 3-(c)). The sequencing of diced RGM249 mRNA proved that the sequence we used for siRNA or shRNA was included in fig 3-(b). Because the silencing of RGM249 induced the downregulation of hTERT expression, this gene may be involved in the oncogenic mechanism. In humans, ~800 miRNAs are predicted to exist ; each of these is thought to have hundreds of targets, suggesting that miRNAs play a critical regulatory role in cellular protein expression . In various human cancers as well as in developmental timing and patterning, differentiation, organogenesis, and stress response, alterations of miRNA profiles have been observed; these miRNAs may act as oncogenes, leading to the specific dysregulation of their target gene transcripts [41, 42]. Although alterations in miRNA expression have been described in cancer and may contribute to other disease states, it remains unclear what drives these alterations. Our findings might pave the way for defining unknown crucial roles that miRNAs located around this locus play in oncogenic regulation through hTERT regulation and, in particular, the role RGM249 may play during this process.
Although the way in which RGM249 regulates hTERT remains unidentified, the physiological expression of RGM249 synchronizes the upregulation of hTERT mRNA, and its excess expression using the RNA-expression vector downregulated hTERT in specific cancer cells. Moreover, its siRNA or shRNA suppressed hTERT. Taken together, RGM249 may function negatively or positively, reversibly or indirectly on unknown mediators upstream of hTERT. Using microarray analysis between transfectants with RGM249 siRNA and those with LacZ siRNA (fig 6), the comparative expression analysis demonstrated that ATM and TP53LP were upregulated, although p53 status (the presence of mutation) in each cell line did not correlate with hTERT mRNA expression, and Ras and several transcriptional factors were downregulated. As described in fig 6, RGM249 may have direct effects on several transcriptional factors. Rapid stress induced by siRNAs may induce upregulation of genes participating in both repairing action and apoptotic or growth inhibition induction. Comparative proteome analysis would be useful to explain the detailed mechanisms underlying these processes. The overexpression of RGM249 in normal cells did not alter hTERT expression level, suggesting that RNA-related molecules may be functionally different in normal cells and de-differentiated cells, and that a key molecule for uncovering the oncogenesis mechanism or a potential factor which could be applied as a therapeutic target may exist in the genome. The findings that both the overexpression and dsRNA knockdown of RGM249 suppressed the expression of hTERT are problematic. In addition, this gene may be indirectly associated with hTERT because its overexpression may not always reveal a normal function and RGM249 microRNAs are not predicted to bind to hTERT even with the use of available bioinformatics. Moreover, this gene is not mutated or homozygously deleted even in HCC cells and tissues. Based on the Celera assembly, it seems to have one base (T) addition at position 188 of the first exon, however, single nucleotide polymorphisms are not specific to cancer. Additional experiments such as promoter analysis are needed to elucidate the epigenetic network in which this gene works.
This novel noncoding gene is strongly expressed in poorly differentiated tumor cells. The genes (mRNAs) targeted commonly by both miRNAs and siRNAs in our study were MAPK6 and PRKCA, which can induce hTERT upregulation (fig 7-(a)). Representative genes (mRNAs) targeted by miRNAs with expression upregulated more than twofold by shRGM249RNA included DICER1, high mobility group A2 (HMGA2) and GAS7, as well as hsa-miR-let 7b, 7d, 7e, 7f, and 7g targets. In particular, RGM249 might be implicated in oncogenic transformation through the inhibition of HMGA2 interacting with let-7 miRNA (fig 7-(b), (c)) . RGM249 was not detected in any 'mirtrons or polycistron' detection program, genome, or miRNA database (MiRscan, PicTar, miRanda, etc.). We are now considering MED18 as a key molecule obtained from microarray analysis (fig 6 and 7) . The mediator head module stimulates basal RNA pol II transcription and enables transcriptional regulation. The head subunits MED8, MED18, and MED20 form a subcomplex (MED8/18/20) with 2 submodules. The highly conserved N-terminal domain of Med8 forms one submodule that binds the TATA box-binding protein (TBP) in vitro and is essential in vivo. Beyer et al. reported that the downregulation of another mediator complex (MED28) expression in NIH3T3 cells results in a significant induction of several genes associated with smooth muscle cell (SMC) differentiation and functions as a repressor of SMC differentiation . A mediator complex is essential for transcription induced by pol II in eukaryotes and is required for the recruitment of the transcription machinery after chromatin remodeling . MED18 may also be involved in the repression of differentiation, transcription, and telomere elongation. Microarray analysis revealed a strong upregulation of MED18 (×24.8), the RNA polymerase II mediator which functions as a coactivator involved in transcription regulation in nearly all RNA polymerase II-dependent genes, and other RNA-related or DNA-binding genes were upregulated (fig 6). In transfectants with MED18 siRNA, both RGM249 and hTERT mRNA were significantly upregulated (Supplement fig 2-(d)). MED18 siRNA transfected into A172 cells significantly upregulated hTERT mRNA and RGM249 expression (Additional file 1-(d)). Co-transfection by both MED18 siRNA and hTERT siRNA upregulated only RGM249 mRNA, suggesting that RGM249 may functionally locate upstream hTERT mRNA, and MED18 may closely interact with RGM249. Overexpressing RGM249 is carried out for the purpose of comprehending the possible suppressive effects of miRNAs and depleting RGM249 is performed for the purpose of understanding the physiological function of this gene. Our results suggest that RGM249 may not be a precursor gene such as Onco-miR but is a Suppressor-miR and might function to lower differentiation in its expressing cells by suppressing cell growth as an independent response to hTERT expression.
Taken together, the enhanced cell death in RGM249 transfectants likely resulted from an hTERT-independent mechanism or the presence of a mediator with crucial responsibility for oncogenesis, arguing that some other targets instead of hTERT might be the real players mediating these contradictory oncogenic and growth suppressive functions. It should be noted that mRNA studies do not always reflect processes at the protein level of gene regulation.
Further studies will be required to investigate these expression levels. Several ncRNAs in the exons we trapped might provide new insights into the coordinated and dynamic mechanisms of development, stress, transformation, and (de)differentiation, resulting in the possible prevention of oncogenesis and suppression of tumor progression.
DNA preparation and analysis
BAC clones (H-11, M-6, and J-21) were selected on the basis of results from a previous experiment (Research Genetics, Toronto, Canada) . BAC DNA was isolated using a modified alkaline lysis protocol . Pulsed-field electrophoresis analysis of the isolated BAC DNA was performed using the CHEF DRII system (Bio-Rad, CA, USA).
Sequence determination and analysis, and construction of a contig of BAC clones (M-6, H-11, and J-21)
Purified BAC DNAs were shared using the AERO-MIST nebulizer (CIS-US, Bedford, MA, USA) and after ligation of adapters, these DNAs were subcloned into the BstXI sites of pSHOT II, which was constructed by subcloning the BstXI sites from pcDNA II, prior to the disclosure of the endosequences of the 3 BAC clones (M-6, H-11, and J-21) (fig 1-(b)). A sequencing library was generated by sonicating BAC and cloning it into a phosphatased EcoRV site in Bluescript; initially, each clone was sequenced using Phil Greene's Phred and Phrap programs.
For all exon trapping experiments, we used the exon-trapping system kit (Gibco-BRL, Gaithersburg MD, USA) and the supplements recommended therein. We shotgunned M-6, H-11, and J-21 into pSPL3, employing BamHI and BglII or BamHI, BglII, and PstI for 2 ways of cutting BAC DNA using restriction enzymes. This method was performed as described previously . Thereafter, the BAC whole sequence or cDNA sequence was identified using a database http://www.ncbi.nlm.nih.gov/. In 1999 to 2005, the sequenced reference human genome draft was not definitely determined on 10p14–15 and altered several times during this interval.
Northern blotting, reverse-transcriptase polymerase chain reaction (RT-PCR), and onestep real-time RT-PCR for confirming the expression of exons in human organs
The probe, consisting of the respective exons, including two exons forming RGM249, radiolabeled by random priming using the Radprime kit (Gibco-BRL), was incubated with human multiple tissue northern blot (Human 12-Lane MTN Blot) (Clontech, Mountain View, CA, USA). One-step real-time RT-PCR experiments were performed using Human Total RNA Master Panels (Clontech) according to methods we previously reported using LightCycler (Roche, Basel, Switzerland) . The primer sequences used are shown in Table 1. Each quantification was confirmed by sequencing with high reproducibility (for RGM249; fig 2-(a), r = 0.998, for hTERT) .
cDNA cloning and exon/intron organization
Using the exons identified by Northern blotting or RT-PCR, we amplified cDNA fragments from human skeletal muscle and liver "Marathon-ready" cDNA, according to the manufacturer's protocol (Clontech). Amplicons were subcloned into the pcDNA3 vector (Invitrogen, Leek, The Netherlands). A comparison of cDNA and genomic sequences revealed a total of 20 exons and detected 3 alternative splice donor sites for 6 exons.
PCR-single-strand conformation polymorphism (SSCP) analysis for the detection of mutations in exons
To screen for point mutations in trapped exons, we performed PCR-SSCP assays using flanking intron primers, as previously described .
Measurement of telomerase activity
To confirm telomerase activity in transfectants after siRNA transfection, a telomerase detection kit, TRAP-eze, (Oncor Inc., Gaithersburg, MD, USA) used for RGM249 cloning and the TeloChaser kit (TOYOBO, Osaka, Japan) used for the expression analysis of shRNA transfectants, were employed for the telomeric repeat amplification protocol assay according to the manufacturer's protocol. Quantification of TRAP was estimated using an image analyzer (CS Analyzer 3) provided by ATTO Corp. (Tokyo, Japan).
Structural prediction of RGM249
To examine whether RGM249 can generate a protein, we performed an in vitro translation assay according to the manufacturer's instructions (Clontech). After examining the motif and homology http://www.ncbi.nlm.nih.gov/genome/seq/BlastGen/BlastGen.cgi?taxid=9606; http://hmmpfam.ddbj.nig.ac.jp/; http://blast.ncbi.nlm.nih.gov/Blast.cgi), we predicted the secondary structure of RGM249 using the Vienna RNA Secondary prediction and comparison program http://www.tbi.univie.ac.at/~ivo/RNA/; http://www.genebee.msu.su/services/rna2_reduced.html). Subsequently, we screened the micro RNA) (miRNA) database among species. To confirm whether the RGM249 gene functions as a miRNA precursor gene, we digested the full-length RNA with a Dicer enzyme using a Dicer siRNA generation kit according to the manufacturer's protocol (Gene Therapy Systems, Inc., CA, USA).
RGM249 mRNA expression in human organs, hepatocytes, hepatoma cell lines, and liver tissues
To confirm the transcriptional level (copy number) in human organs, liver-related normal cells, and immortalized cells, we performed one-step real-time RT-PCR using human hepatocytes (KAC, Kyoto, Japan), hepatoma cell lines (Tohoku University, Sendai, Japan), and 74 surgically resected liver tissue samples, including 33 pairs of hepatomas and corresponding adjacent tissues. We measured each copy number 3 times and calculated the average.
Effect of RGM249 overexpression on expression of other genes in RGM249 transfectants
After ligasing a synthesized RGM249 (Operon Biotechnologies, Tokyo, Japan) into the pRNAT/U6.1/Neo/GFP vector (GeneScript Corp., Piscataway, NJ, USA) digested at the KpnI/SacI site, it was transfected into several cell lines (A172, T98G, HLF, HT1080, HeLa, SW480, Huh7, HepG2, HMc-Li7, KMST-6, SUSM-1, and NIH3T3), with NH-12 as senescence-programmed cell lines, and keratinocytes and TIG 1–20 as primary cells. Immortalized cells were transfected with FuGene (Roche, Basel, Switzerland) to obtain stable cells overexpressing RGM249. After enhanced green fluorescent protein expression in the presence of a neomycin-resistant gene in all cell types was confirmed, β-gal staining (Senescence Detection kit; BioVision, CA, USA) and RT-PCR were performed. Primary cultured cells were electroporated using Nucleofector technology (Amaxa AG, Cologne, Germany) under optimal conditions. At least 5 transfectants were examined in mock or RGM249 transfectants. The respective gene expression (GE) ratio to β-actin mRNA in mocks or RGM249 tranfectants was analyzed using the Mann-Whitney test. Between both transfectants, FACS analysis of HLF cells was performed using Guava EasyCyte Mini according to the manufacturer's instructions (GE Healthcare UK Ltd., Buckinghamshire, England).
Transfection of total dsRNA and gene-specific siRNA derived from the RGM249 gene
dsRNA derived from the RGM249 sequence was generated from the expression vector with RGM249 in the MCS of Litmas28i (NEB, MA, USA) using an in vitro transcription method (T7 RiboMAX Express Large Scale RNA Production System, Promega, Madison, WI, USA) and was transfected into 4 hepatoma cell lines (HLF, Alexander, HepG2, and HMc-Li7) at a concentration range of 10–40 μg/ml medium. Subsequently, to avoid unspecific inhibition of dsRNA through interferon (INF) activation [50, 24], stealth siRNAs were designed on-line using RNAi designer provided by Invitrogen (Tokyo, Japan). At the optimal concentration (25 nM) of siRNA to the recipient (HLF), which expresses it strongly, we repeatedly performed transfections to examine the correlation of telomerase-related factors (hTR, TEP1, dyskerin, RPL22, etc.) with RGM249 expression . The expression of telomerase-related genes was estimated by RT-PCR after transfection of dsRNA (fig 4-(a)).
Two kinds of siRNAs (siRNA ① and siRNA ②) were designed using the Stealth RNAi designer and synthesized by Invitrogen (Invitrogen Corp, Carlsbad, CA, USA). The sequences are shown in Table 2. To suppress RGM249 expression with siRNA, we performed the following transfection procedures in HLF cells: transfection using Stealth siRNA ①, transfection employing Stealth siRNA ②, and transfection with both Stealth siRNA ① and Stealth siRNA ② at a final concentration of 25 nM as an optimal condition. We collected transfectants 24 hours after transfection for RNA extraction and confirmed the inhibitory effect on RGM249 expression by one-step RT-PCR. HLF cells treated with FuGene Reagent without siRNAs and with Stealth LacZ siRNA were used as controls. Transfectants were collected 48 hours, 1 week, and 2 weeks after transfection. Subsequently, we confirmed the suppressive effect of shRNA with the same functional sequence as that of siRNA on RGM249 expression, using the BLOCK-iT RNAi
Expressing vector (Invitrogen) which was designed to generate shRNA in transfectants. To examine the effect on hTERT expression 7 to 10 days later after selection with 200 μg/ml Zeocin, an shRNA-generating lentiviral vector with the functional sequence of RGM249 (fig 3-(b)) was transfected. Since the mock of the lentiviral vector was linearized according to the manufacturer's instructions, we used transfection reagents without RNA and the vector, LacZ shRNA as the controls, and wtRGM249 or mtRGM249 (mt-1 or mt-2) with one base (T or C) deletion in the wtRGM249 sequence. We evaluated GE in the transfectants semiquantitatively. Stealth siRNA and shRNA sequences for RGM249 are shown in Table 2.
Gene expression profile by microarray analysis
Microarray analysis was performed for total RNA purified from transfectants with LacZ shRNA and from those with RGM249 shRNA following confirmation of the downregulation of both RGM249 and hTERT mRNA expression by RT-PCR after strict quality control in accordance with the guidelines of the Hokkaido System Science (Sapporo, Japan). Data for genes that had changed more than twofold were additionally analyzed using Gene Ontology Tree Machine http://bioinfo.vanderbilt.edu/gotm/ and typical GE changes are shown in fig 6 and the left fig in Additional file 2. miRNA array analysis was performed using the same samples as the microarray analysis samples (Filgen, Inc., Nagoya, Japan). The results are shown in fig 7(b) and the right fig in Additional file 2. Genes targeted by miRNAs were investigated using DIANA-microT, Pictar, TargetScan, and miRanda http://diana.pcbi.upenn.edu/cgi-bin/TargetCombo.cgi.
telomerase RNA component
Bacterial artificial chromosome
human telomerase reverse transcriptase
short interfering RNA
short hairpin RNA
ribosomal protein L22
Greider CW, Blackburn EH: A telomeric sequence in the RNA of Tetrahymena telomerase required for telomere repeat synthesis. Nature. 1989, 337: 331-337. 10.1038/337331a0
Morin GB: Recognition of a chromosome truncation site associated with alpha-thalassaemia by human telomerase. Nature. 1991, 353: 454-456. 10.1038/353454a0
Fu D, Collins K: Purification of human telomerase complexes identifies factors involved in telomerase biogenesis and telomere length regulation. Mol Cell. 2007, 28: 775-785. 10.1016/j.molcel.2007.09.023.
Cohen SB, Graham ME, Lovrecz GO, Bache N, Robinson PJ, Reddel RR: Protein composition of catalytically active human telomerase from immortal cells. Science. 2007, 30: 1850-1853. 10.1126/science.1138596.
Kim NW, Piatyszek MA, Prowse KR, Harley CB, West MD, Ho PL, Coviello GM, Wright WE, Weinrich SL, Shay JW: Specific association of human telomerase activity with immortal cells and cancer. Science. 1994, 266: 2011-2015. 10.1126/science.7605428
Swiggers SJ, Nibbeling HA, Zeilemaker A, Kuijpers MA, Mattern KA, Zijlmans JM: Telomerase activity level, but not hTERT mRNA and hTR level, regulates telomere length in telomerase-reconstituted primary fibroblasts. Experimental Cell Res. 2004, 297: 434-443. 10.1016/j.yexcr.2004.03.028.
Chang JT, Chen YL, Yang HT, Chen CY, Cheng AJ: Differential regulation of telomerase activity by six telomerase subunits. Eur J Biochem/FEBS. 2002, 269: 3442-3450. 10.1046/j.1432-1033.2002.03025.x.
Chen CJ, Tsai NM, Liu YC, Ho LI, Hsieh HF, Yen CY, Harn HJ: Telomerase activity in human hepatocellular carcinoma: parallel correlation with human telomerase reverse transcriptase (hTERT) mRNA isoform expression but not with cell cycle modulators or c-Myc expression. Eur J Surg Oncol. 2002, 28: 225-234. 10.1053/ejso.2001.1237
Bachand F, Autexier C: Functional regions of human telomerase reverse transcriptase and human telomerase RNA required for telomerase activity and RNA-protein interactions. Mol Cell Biol. 2001, 21: 1888-1897. 10.1128/MCB.21.5.1888-1897.2001
de Souza Nascimento P, Alves G, Fiedler W: Telomerase inhibition by an siRNA directed against hTERT leads to telomere attrition in HT29 cells. Oncol Rep. 2006, 16 (2): 423-428.
Holt SE, Wright WE, Shay JW: Multiple pathways for the regulation of telomerase activity. Eur J Cancer. 1997, 33: 761-766. 10.1016/S0959-8049(97)00066-X
Oshimura M, Barrett JC: Multiple pathways to cellular senescence: role of telomerase repressors. Eur J Cancer. 1997, 33: 710-715. 10.1016/S0959-8049(97)00090-7
Horikawa I, Barrett JC: Transcriptional regulation of the telomerase hTERT gene as a target for cellular and viral oncogenic mechanisms. Carcinogenesis. 2003, 24: 1167-1176. 10.1093/carcin/bgg085
Stark M, Hayward N: Genome-wide loss of heterozygosity and copy number analysis in melanoma using high-density single-nucleotide polymorphism arrays. Cancer Res. 2007, 15: 2632-2642. 10.1158/0008-5472.CAN-06-4152.
Steck PA, Ligon AH, Cheong P, Yung WK, Pershouse MA: Two tumor suppressive loci on chromosome 10 involved in human glioblastomas. Genes Chromosomes Cancer. 1995, 12: 255-261. 10.1002/gcc.2870120404
Sasaki M, Honda T, Yamada H, Wake N, Barrett JC, Oshimura M: Evidence for multiple pathways to cellular senescence. Cancer Res. 1994, 54: 6090-6093.
Pirker C, Holzmann K, Spiegl-Kreinecker S, Elbling L, Thallinger C, Pehamberger H, Micksche M, Berger W: Chromosomal imbalances in primary and metastatic melanomas: over-representation of essential telomerase genes. Melanoma Res. 2003, 13: 483-492. 10.1097/00008390-200310000-00007
Miura N, Onuki N, Rathi A, Virmani A, Nakamoto S, Kishimoto Y, Murawaki Y, Kawasaki H, Hasegawa J, Oshimura M, Travis WD, Gazdar AF: hTR repressor-related gene on human chromosome 10p15.1. Br J Cancer. 2001, 85: 1510-1514. 10.1054/bjoc.2001.2121
Nishimoto A, Miura N, Horikawa I, Kugoh H, Murakami Y, Hirohashi S, Kawasaki H, Gazdar AF, Shay JW, Barrett JC, Oshimura M: Functional evidence for a telomerase repressor gene on human chromosome 10p15.1. Oncogene. 2001, 20: 828-835. 10.1038/sj.onc.1204165
Kim UJ, Birren BW, Slepak T, Mancino V, Boysen C, Kang HL, Simon MI, Shizuya H: Construction and characterization of a human bacterial artificial chromosome library. Genomics. 1996, 34: 213-218. 10.1006/geno.1996.0268
Mejia JE, Monaco AP: Retrofitting vectors for Escherichia coli-based artificial chromosomes (PACs and BACs) with markers for transfection studies. Genome Res. 1997, 7: 179-186. 10.1101/gr.7.2.179
Hejna JA, Johnstone PL, Kohler SL, Bruun DA, Reifsteck CA, Olson SB, Moses RE: Functional complementation by electroporation of human BACs into mammalian fibroblast cells. Nucleic Acids Res. 1998, 26: 1124-1125. 10.1093/nar/26.4.1124
Nisson PE, Ally A, Watkins PC: Protocols for trapping internal and 3'-terminal exons. PCR Methods Appl. 1994, 4 (1): S24-S39.
Church DM, Stotler CJ, Rutter JL, Murrell JR, Trofatter JA, Buckler AJ: Isolation of genes from complex sources of mammalian genomic DNA using exon amplification. Nature Genet. 1994, 6: 98-105. 10.1038/ng0194-98
Burn TC, Connors TD, Klinger KW, Landes GM: Increased exon-trapping efficiency through modifications to the pSPL3 splicing vector. Gene. 1995, 161: 183-187. 10.1016/0378-1119(95)00223-S
Tuschl T, Zamore PD, Lehmann R, Bartel DP, Sharp PA: Targeted mRNA degradation by double-stranded RNA in vitro. Genes Dev. 1999, 13: 3191-3197. 10.1101/gad.13.24.3191
Hamilton AJ, Baulcombe DC: A species of small antisense RNA in posttranscriptional gene silencing in plants. Science. 1999, 286: 950-952. 10.1126/science.286.5441.950
Martinez J, Patkaniowska A, Elbashir SM, Harborth J, Hossbach M, Urlaub H, Meyer J, Weber K, Vandenburgh K, Manninga H, Scaringe SA, Luehrmann R, Tuschl T: Analysis of mammalian gene function using small interfering RNAs. Nucleic Acids Res Suppl. 2003, 333-3
Harborth J, Elbashir SM, Vandenburgh K, Manninga H, Scaringe SA, Weber K, Tuschl T: Sequence, chemical, and structural variation of small interfering RNAs and short hairpin RNAs and the effect on mammalian gene silencing. Antisense Nucleic Acid Drug Dev. 2003, 13: 83-105. 10.1089/108729003321629638
Brummelkamp TR, Bernards R, Agami R: A system for stable expression of short interfering RNAs in mammalian cells. Science. 2002, 296: 550-553. 10.1126/science.1068999
Le S, Sternglanz R, Greider CW: Identification of two RNA-binding proteins associated with human telomerase RNA. Mol Biol Cell. 2000, 11: 999-1010.
Fujimori M, Tokino T, Hino O, Kitagawa T, Imamura T, Okamoto E, Mitsunobu M, Ishikawa T, Nakagama H, Harada H, Yagura M, Matsubara K, Nakamura Y: Allelotype study of primary hepatocellular carcinoma. Cancer Res. 1991, 51: 89-93.
Bigner SH, Mark J, Bigner DD: Cytogenetics of human brain tumors. Cancer Genet Cytogenet. 1990, 47: 141-154. 10.1016/0165-4608(90)90024-5
Isshiki K, Elder DE, Guerry D, Linnenbach AJ: Chromosome 10 allelic loss in malignant melanoma. Genes Chromosomes Cancer. 1993, 8: 178-184. 10.1002/gcc.2870080307
Karlbom AE, James CD, Boethius J, Cavenee WK, Collins VP, Nordenskjold M, Larsson C: Loss of heterozygosity in malignant gliomas involves at least three distinct regions on chromosome 10. Human Genet. 1993, 92: 169-174. 10.1007/BF00219686.
Anderson CJ, Hoare SF, Ashcroft M, Bilsland AE, Keith WN: Hypoxic regulation of telomerase gene expression by transcriptional and post-transcriptional mechanism. Oncogene. 2006, 5: 61-69.
Chendrimada TP, Finn KJ, Ji X, Baillat D, Gregory RI, Liebhaber SA, Pasquinelli AE, Shiekhattar R: MicroRNA silencing through RISC recruitment of eIF6. Nature. 2007, 447: 823-828. 10.1038/nature05841
Kiriakidou M, Tan GS, Lamprinaki S, De Planell-Saguer M, Nelson PT, Mourelatos Z: An mRNA m7G cap binding-like motif within human Ago2 represses translation. Cell. 2007, 129: 1141-1151. 10.1016/j.cell.2007.05.016
Bentwich I, Avniel A, Karov Y, Aharonov R, Gilad S, Barad O, Barzilai A, Einat P, Einav U, Meiri E, Sharon E, Spector Y, Bentwich Z: Identification of hundreds of conserved and nonconserved human microRNAs. Nature Genet. 2005, 37: 766-770. 10.1038/ng1590
Bartel DP, Chen CZ: Micromanagers of gene expression: the potentially widespread influence of metazoan microRNAs. Nature Reviews. 2004, 5: 396-400.
Michael MZ, SM OC, van Holst Pellekaan NG, Young GP, James RJ: Reduced accumulation of specific microRNAs in colorectal neoplasia. Mol Cancer Res. 2003, 1: 882-891.
Tong AW: Small RNAs and non-small cell lung cancer. Curr Mol Medicine. 2006, 6: 339-349. 10.2174/156652406776894554.
Mayr C, Hemann MT, Bartel DP: Disrupting the pairing between let-7 and Hmga2 enhances oncogenic transformation. Science. 2007, 315: 1576-1579. 10.1126/science.1137999
Wang X, Wang X: Systematic identification of microRNA functions by combining target prediction and expression profiling. Nucleic Acids Res. 2006, 34: 1646-1652. 10.1093/nar/gkl068
Lariviere L, Geiger S, Hoeppner S, Rother S, Strasser K, Cramer P: Structure and TBP binding of the Mediator head subcomplex Med8-Med18-Med20. Nature Structural Mol Biol. 2006, 13: 895-901. 10.1038/nsmb1143.
Beyer KS, Beauchamp RL, Lee MF, Gusella JF, Naar AM, Ramesh V: Mediator subunit MED28 (Magicin) is a repressor of smooth muscle cell differentiation. J Biol Chem. 2007, 282: 32152-32157. 10.1074/jbc.M706592200
He Q, Battistella L, Morse RH: Mediator requirement downstream of chromatin remodeling during transcriptional activation of CHA1 in yeast. J Biol Chem. 2008, 283 (9): 5276-5286. 10.1074/jbc.M708266200
Miura N, Nakamura H, Sato R, Tsukamoto T, Harada T, Takahashi S, Adachi Y, Shomori K, Sano A, Kishimoto Y, Ito H, Hasegawa J, Shiota G: Clinical usefulness of serum telomerase reverse transcriptase (hTERT) mRNA and epidermal growth factor receptor (EGFR) mRNA as a novel tumor marker for lung cancer. Cancer Sci. 2006, 97: 1366-1373. 10.1111/j.1349-7006.2006.00342.x
Takenoshita S, Tani M, Mogi A, Nagashima M, Nagamachi Y, Bennett WP, Hagiwara K, Harris CC, Yokota J: Mutation analysis of the Smad2 gene in human colon cancers using genomic DNA and intron primers. Carcinogenesis. 1998, 19: 803-807. 10.1093/carcin/19.5.803
Sharp PA, Zamore PD: Molecular biology. RNA interference. Science. 2000, 287: 2431-2433. 10.1126/science.287.5462.2431
This work was supported by a Grant-in-Aid from Takeda Science Foundation and by a Grant-in-Aid for Cancer Research from the Osaka Cancer Study Group. We thank Dr. Norimitsu Inoue and Dr. Zinne Andrew of the McDermott Center, University of Texas Southwestern Medical Center, USA, for screening the BAC library and Dr. Anthony P. Monaco for the kind gift of the exon-trapping vector. Cell lines used in this study were all provided by the Cell Resource Center for Biomedical Research, Institute of Development, Aging and Cancer, Tohoku University, Japan. All the PCR primers were designed by INTEC Web and Genome Informatics, Corporation (Tokyo, Japan).
The authors declare that they have no competing interests.
MN and SR carried out the molecular genetic studies. WJ and DH participated in the sequence alignment and drafted the manuscript. MJ kindly provided the exon trapping vector. SG kindly provided the clinical specimens. MN and VA participated in the study design and performed the statistical analysis. Other authors conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript. All contributors who do not meet the criteria for authorship are listed in Acknowledgements.
Norimasa Miura, Reina Sato contributed equally to this work.
Electronic supplementary material
Additional file 1: A comparison of hTERT mRNA expression between mock and RGM249 transfectants. In response to upregulated RGM249 mRNA, hTERT mRNA did not induce any significant change in senescence-programmed cells (keratinocytes, TIG1-20, and NH-12) except for KMST-6 cells. The dotted line indicates hTERT mRNA expression in parental cells. Six transfectants were examined in both mock and RGM249 transfectants. (b) Actin mRNA and RPL22 mRNA showed no significant difference among untransfected cells and mock or RGM249 transfectants in the 4 cell lines described in (a). (c) Comparison of hTERT mRNA expression between mock and RGM249 transfectants. In response to upregulated RGM249, hTERT mRNA induced a significant change in suppression in (left) HLF cells (poorly differentiated hepatoma) and (right) A172 cells (glioblastoma). The dotted line indicates mean mRNA expression in untransfected parental cells. hTR mRNA and b-actin mRNA showed no significant differences in b-actin, hTERT, RGM249, and hTR among untransfected cells and mock or RGM249 transfectants. (d) By transfection of MED18 siRNA or MED18/hTERT siRNA into A172 cells, an insight into the regulative network among 3 genes (RGM249, MED18, and hTERTmRNA) could be obtained. The number of respective inductions was more than 5. Following transfection of MED18 siRNA, both TERT mRNA and RGM249 mRNA were significantly upregulated (P < 0.01). Following co-transfection of MED18 siRNA and TERT siRNA, RGM249 mRNA was significantly upregulated (P < 0.05). (PDF 306 KB)
Additional file 2: A representative case of microarray analysis and microRNAarray analysis. Microarray analysis (a) and miRNA array analysis (b) were performed for total RNA purified from transfectants with LacZ shRNA and from those with RGM249 shRNA. The representative analyses are demonstrated, respectively. (PDF 533 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.