Molecular and Cellular Biology, November 2008, p. 6731-6745, Vol. 28, No. 22
0270-7306/08/$08.00+0 doi:10.1128/MCB.02103-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
INSERM U908, IFR 147 and Université Lille1, 59655 Villeneuve d'Ascq, France,1 Institut de Génétique Moléculaire, CNRS UMR 5535, Université Montpellier II, IFR 122, Montpellier, France,2 CNRS UMR 8161, IBL, Universités de Lille I et Lille II, Institut Pasteur, Lille, France3
Received 26 November 2007/ Returned for modification 14 January 2008/ Accepted 29 July 2008
|
|
|---|
|
|
|---|
Genomic imprinting plays a critical role in modulating gene expression during embryogenesis and normal development (36). It consists of an epigenetic modification that results in the silencing of a specific allele, depending on its parental origin. Most of the 80 mammalian imprinted genes identified thus far are organized into clusters (33). It has been suggested that the clustering of genes enables sharing of cis-acting elements located several kilobase pairs away (21). Differentially methylated regions associated with the imprinted clusters play a crucial role in the imprinted expression patterns, which are hence called imprinting control regions (ICRs) (22, 38). The H19/IGF2 locus is a well-characterized cluster. Both genes share a common set of enhancers located downstream from the H19 gene. The ICR, located 2 kb upstream of the H19 promoter, controls the monoallelic expression of the H19 and IGF2 genes by insulating communication between the 3' enhancers and the IGF2 promoter. The chromatin insulator property of the H19/IGF2 ICR is regulated by a chromatin insulator protein, CTCF (for CCTC-binding factor), which binds only the nonmethylated maternal allele. On the paternal allele, the ICR methylation does not allow CTCF binding and leads to IGF2 expression (for a review, see reference 24). The IGF2 gene itself also contains differentially methylated regions (DMRs), with DMR1 being a methylation-sensitive silencer and DMR2 being a methylation-sensitive activator (7, 30). Elsewhere, the role of CTCF in chromatin loops is now well demonstrated. Indeed, CTCF binding in the maternal ICR is able to regulate its interaction with matrix attachment region 3 (MAR3) and DMR1 at IGF2, thus forming a tight loop around the maternal IGF2 locus which may contribute to its silencing. These interactions may contribute to restrict the physical access of distal enhancers to IGF2 promoter (19, 30, 50). Additional mechanisms exist for an imprint mark, such as chromatin composition, organization, and histone acetylation or methylation state (12, 14), even if DNA methylation is by far the best candidate (4).
A more recently identified characteristic of imprinted genes is their association, in some cases, with noncoding antisense transcripts (ncRNAs), which have been suggested to constitute a new epigenetic regulatory system. Indeed, for the maternally expressed Igf2r and UBE3A genes, overlapping antisense transcripts have been found and are oppositely imprinted with respect to the protein coding gene. It has been proposed that antisense transcripts serve to regulate overlapping genes by promoter or transcript occlusion or by competing with these loci for regulatory elements such as transcription factors or enhancers (39, 40, 44, 54). Elsewhere, in the imprinted cluster at the distal end of the mouse chromosome 7 (orthologous to the human chromosome 11), the ICR located in the KCNQ1 gene contains the promoter of a long paternally expressed antisense transcript called KCNQ1OT1. This transcript takes part in gene silencing of the KCNQ1 domain, and its interruption by the insertion of a polyadenylation sequence caused a loss of methylation spreading in cultured cells (31, 47).
These ncRNAs are not yet clearly classified, but categories with known gene regulatory functions are emerging. It includes intergenic transcripts that regulate local chromatin activity, cis-acting long ncRNAs such as Xist involved in chromosome inactivation, and the imprinted Air and KCNQ1OT1 ncRNAs involved in domain silencing (15, 32).
In the present study we have identified and characterized a new transcript within the human and mouse H19/IGF2 locus. This long intergenic RNA that we named 91H is antisense to H19. We assessed the effect of the 91H knockdown by RNA interference on genomic imprinting and on H19 and IGF2 gene expression. The results reported here add further complexity to the locus regulation and are attractive within the context of imprinting maintenance and/or cancer development.
|
|
|---|
RNA extraction. Total RNA was extracted with Tri-Reagent (Euromedex) and treated for 1 h at 25°C with amplification-grade RNase-free DNase I (Invitrogen) or isolated by using a Nucleospin RNAII isolation kit (Macherey-Nagel).
Reverse transcription. Total reverse transcriptions were performed as follows: 1 µg of RNA, 4 µl of buffer containing random hexamers, and 1 µl of reverse transcriptase (RT; Bio-Rad) were incubated for 5 min at 70°C, for 30 min at 37°C, and for 5 min at 85°C in a final volume of 20 µl. For strand-specific reverse transcription, we used the Thermoscript RT-PCR System (Invitrogen). A total of 2 µg of total RNA was incubated 5 min at 65°C with a 1 µM concentration of the primer and a 10 µM concentration of deoxynucleoside triphosphate in a final volume of 12 µl. We then added 4 µl of buffer (5x), 1 µl of dithiothreitol (0.1 M), 1 µl of RNase Out (40 U/µl), and 2 µl of H2O. The sample was divided into two equal parts. Thermoscript (1 µl) was added in one of them (+RT), and 1 µl of H2O was added to the other to constitute the negative control (-RT). The reverse transcription was performed 20 min at 60°C to improve specificity. The reaction was stopped by 5 min at 85°C, and cDNA/RNA duplexes were removed by the addition of 1 µl of RNase H for 15 min at 37°C. The primer used for H19-specific RT is referred to as "H19-specific RT," and the primer used for 91H specific RT is referred to as U'sense (Table S1 in the supplemental material).
Real-time RT-PCR. Real-time PCR amplifications were performed by using a QuantiTect Sybr green PCR kit (Qiagen) with 2 µl of cDNA and 500 nM concentrations of the primers. The primers used for the H19 transcript, RPLP0 (for human acidic ribosomal phosphoprotein P0), GAPDH, U3snoRNA, MRPL23 mRNA, and 91H RNA (D and I) are described in Table S1 in the supplemental material. Other primer sequences are available upon request. The subsequent PCR conditions consisted of 40 cycles of 95°C for 15 s, 60°C for 20 s, and 72°C for 30 s. The data were analyzed by using an MX4000 PCR system software (Stratagene) with the Sybr green option (with dissociation curves). Standard curves were performed on serial dilutions of a PCR product for H19 and on serial dilutions of genomic human DNA for RPLP0. Values were obtained by the calculation methodology recommended by Pfaffl (34): ratio = (NC – b/a)target/(NC – b/a)reference, where a is the slope of the standard curve and b is the ordinate of origin.
5'RACE- and 3'RACE-PCR. Transcription start site and 3' end of the 91H transcript were determined by using RACE (rapid amplification of cDNA ends)-PCR experiments. A Smart RACE cDNA amplification kit (Clontech) was used according to the manufacturer's instructions with cDNA from T47D cells. Primers used for PCR and nested PCR were available upon request. PCR products were amplified with Titanium Taq DNA polymerase. RACE-PCR products were characterized by cloning and sequencing. For the 5'RACE, since the transcript is 120 kb in length, we performed the reverse transcription with a specific primer and not with the conventional oligo(dT). The genomic sequence upstream of the transcription start site was tested for promoter elements by using the TFsearch (http://molsun1.cbrc.aist.go.jp/research/db/TFSEARCH.html) (16) and AliBaba (http://www.gene-regulation.com/pub/programs.html) (13) programs.
Analysis of H19 and Igf2 restriction fragment length polymorphisms. PCRs encompassing AluI and ApaI sites for H19 and IGF2, respectively, were performed on genomic DNA (gDNA), cDNA of nontransfected cells, or cDNA of small interfering RNA (siRNA)-transfected cells. The primers used were HP1 and HP2 for H19 and Igf2ApaI sense and Igf2ApaI antisense for IGF2 (see Table S1 in the supplemental material). PCR conditions were carried out in the following manner: 38 cycles of 94°C for 30 s, 60°C for 30 s, and 72°C for 30 s for IGF2 and 38 cycles of 94°C for 30 s, 59°C for 30 s, and 72°C for 40 s for H19. Then, 15 µl of the PCR product was used for digestions with 2 µl of enzyme buffer 10x, 0.5 µg of bovine serum albumin/µl, and 10 U of the appropriate enzyme for 1 h at 37°C. For AluI polymorphism, digestion of the gDNA-derived PCR product yielded 228-, 128-, and 100-bp bands indicating heterozygosity. The cDNA restriction pattern corresponds to a 148-bp fragment in the allele that lacks the AluI site or two bands of 100 and 48 bp for the allele with the recognition site. For IGF2, ApaI digestion generates a 292-bp fragment in the allele without the recognition site or 218- and 74-bp fragments in the allele with the site.
Isolation of nuclei. Nuclei were isolated from human T47D cells and from proliferating and confluent murine C2C12 cells in growth medium and after 1 and 3 days in differentiation medium as follows. First, 106 cells were washed twice with phosphate-buffered saline (PBS). The cell layer was scraped into 500 µl of PBS and centrifuged 5 min at 1,200 x g. The supernatant was cleared, and this wash step was performed twice. The pellet was estimated to a volume (V) and resuspended in 5V of buffer A (TEA, 10 mM; KCl, 10 mM, MgCl2, 2.5 mM; dithiothreitol, 0.5 mM). After centrifugation under previously described conditions, the pellet was resuspended in 2V of buffer A and disrupted by using a Dounce homogenizer. The separation of nuclei from cytoplasmic compartments was done by centrifugation for 5 min at 3,000 x g. The supernatant containing cytoplasmic extracts was resuspended in 5V of Tri-Reagent. The pellet containing the nuclear extracts was washed in 5V buffer A and centrifuged as previously. The pellet was finally resuspended in 5V of Tri-Reagent. Both extracts were maintained at 4°C before RNA extractions. The efficiency of the fractionation was assessed by quantifying two cytoplasmic mRNAs (RPLP0 and GAPDH) and a nuclear RNA (snoRNA U3). We analyzed the RNA levels of two cytoplasmic control genes. About 80% of these RNA are detected in the cytoplasm. These data demonstrate the efficient separation of cellular compartments and the absence of contamination from one compartment to the other.
gDNA isolation. A total of 107 cells were washed twice with PBS. The cell layer was scraped into 500 µl of physiologic water and then centrifuged for 5 min at 1,000 x g. The supernatant was resuspended in 3 ml of STE buffer (Tris-HCl, 10 mM [pH 8.0]; NaCl, 100 mM; EDTA, 1 mM [pH 8.0]). Proteinase K was added at 200 µg/ml. The sample was incubated for 10 min at 37°C. Sodium dodecyl sulfate at a final concentration of 1% was added, and the sample was incubated again for 20 min. gDNA was extracted by addition of an equal volume of saturated phenol and centrifugation for 10 min at 10,000 x g. The aqueous phase was collected, and 0.2 M sodium acetate at and 2 volumes of ethanol were added. DNA was rolled up with a Pasteur pipette by mixing progressively the two phases. DNA was finally resuspended in 300 µl of H2O and stoked at 4°C.
RNA interference. RNA interference was carried out by using synthetic siRNA duplexes, as described by Elbashir et al. (10). Two synthetic siRNA duplexes (si91H 1 and si91H 2) corresponding to the 91H RNA sequences 5'-GGCGUCAUUCUGAUGGGACTT-3' and 5'-UUCAGGAGCUUAAGAUGCUTT-3', respectively were used to inhibit 91H RNA expression. A synthetic siRNA duplex (siGFP) corresponding to the green fluorescent protein mRNA sequence (5'-GCUGACCCUGAAGUUCAUCTT-3') and a sequence corresponding to an intergenic region located at the H19/IGF2 locus (+108603 bp of the H19 transcription start site; 5'-CGUGGGUGGAUGCAUGGAUTT-3') were used as a negative control. The siRNA duplexes were purchased from Eurogentec. Cells were grown on coverslips in six-well plates and transfected with 400 pmol of siRNAs using Jetsi (Eurogentec), as recommended by the manufacturer. To monitor the transfection efficacy, a tagged-siRNA duplex was transfected in parallel, and the transfection rate was evaluated by fluorescence-activated cell sorting and corresponds to 80 to 90% of transfected cells. After transfection, cells were lysed for total RNA isolation.
|
|
|---|
![]() ![]() View larger version (44K): [in a new window] |
FIG. 1. Detection of an intergenic transcriptional activity at the H19 locus in cancer cell lines. (A) Map of the H19 gene locus. The relative positions of the real-time PCR amplification fragments (A to R) are depicted in the map. CS represents conserved regions between mice and humans that possess enhancer activities (17). Black arrow represents a known EST that was identified after alignment of the human chromosome 11 sequence with EST databases. EST, BX 377296 cDNA source placenta. (B) Quantification of RNA levels were determined by real-time RT-PCR on total RNA samples in T47D and BT20 cell lines. The results were normalized to GAPDH expression levels. Agarose gel shows the amplification products obtained after the real-time PCR on T47D cells. Note that amplifications E to G may also potentially quantify H19 precursors RNAs. (C) Strand-specific RT in two regions showing a decrease of 91H transcription (region 1 and region 2). Antisense-specific RT analyses were performed with the primer C sense in region 1 and K sense in region 2. Transcription was detected by PCR using several contiguous primer pairs (D1, D, and D2 for region 1 and L, M1, M, and M2 for region 2). Results, expressed in terms of the cycle threshold (Ct) are indicated in the lower part of panel C. Cycle threshold values that are equivalent for all primer pairs indicate that there is no transcriptional arrest in both tested regions. Primer sequences are available in Table S1 in the supplemental material.
|
Transcriptional orientation of the intergenic transcriptional activity. We analyzed the transcriptional orientation at four positions flanking the H19 gene, as well as within the H19 gene. The principle of the experiments that we performed is described in Fig. 2A (see also Fig. S1A in the supplemental material). Briefly, primer AS allows the detection of transcripts in an antisense orientation to H19, while primer S is sense specific. These strand-specific primers were used for reverse transcription on total RNA samples from T47D cells. Real-time PCR amplifications were then performed on each side of the RT primers. At all positions analyzed around (see Fig. S2A in the supplemental material) or within (Fig. 2B and see Fig. S1B in the supplemental material) the H19 gene, only antisense-specific RT primers (AS primer) gave PCR amplifications. We conclude that, in humans, the detected transcripts are antisense to the H19 gene. This result is consistent with the data available from the EST described above. It should be noted that, within the H19 gene, no precursor of H19 (S primer) could be detected in the same experimental conditions where the antisense transcription is found (Fig. 2B).
![]() View larger version (33K): [in a new window] |
FIG. 2. Transcriptional orientation of the intergenic transcriptional activity at the H19 locus. (A) Map of the H19 locus. Black lines indicate the five regions (named S to W) where antisense transcription has been established using strand-specific reverse transcription and real-time PCR. Below the map of the H19 locus, the strategy used to analyze transcriptional orientation at the U region is depicted. The same strategy was used to analyze the S, T, V, and W regions. The locations of strand-specific primers used for reverse transcription (primer AS is sense to H19, and primer S is antisense to H19) and PCR amplifications (U and U' amplifications) are indicated. Strand-specific primers were used for reverse transcription on total RNA samples from T47D cells. Real-time PCRs were then performed on each side of the RT priming region to identify the orientation of 91H (U and U' PCR amplifications). (B) Agarose gel of PCR products obtained for the U region with the different RT-PCR combinations. Only RT analyses performed to detect antisense RNA (primer AS) gave significant PCR amplification and only when PCR amplified by the U' primer pair. To exclude any gDNA contaminations, we analyzed RT reactions performed without RT (-RT), and we ascertained the absence of self-priming analyzing reverse transcription performed without any strand-specific primers (H2O). The results obtained for the S, T, V, and W regions can be seen in Fig. S2A in the supplemental material. The primer sequences are available in Table S1 in the supplemental material.
|
The antisense H19 transcription defines a new transcript: the 91H RNA.
To further characterize the antisense H19 transcription, we then mapped its 5' and 3' limits. In humans, the 5' limit is located within intron 1 of the MRPL23 gene (Fig. 3A and B), a gene located
40 kb downstream of H19, which is transcribed in an opposite transcriptional orientation (Fig. 2A). Indeed, on total cDNA, we detected transcripts in the 3' flanking region of MRPL23 gene (MrpA amplification), as well as within the MRPL23 region (MrpC and MrpD amplifications), but not in the promoter region of MRPL23 gene (MrpG amplification) (Fig. 3B, total cDNA column). PCR amplifications performed on strand-specific RT with a MrpB primer (Fig. 3A and B) formally exclude any interference with MRPL23 pre-messengers and showed that the 5' limit of the antisense H19 transcription is located within MRPL23 intron 1, between MrpE and MrpF amplifications, about 3,000 to 1,900 bp downstream of the MRPL23 transcription start site. To define the precise beginning of the 91H transcript, we performed 5'RACE-PCR experiments. We performed PCR with specific primers, and we verified the specificity of the amplification by nested PCR. Cloning and sequencing of the PCR fragments obtained (Fig. 3C) confirmed that the transcription start site was located in intron 1 of the MRPL23 gene at bp 47348 of the transcription start site of H19. To identify transcriptional regulatory sequences, we analyzed the 5'-flanking region from 500 bp upstream of the transcription start site with a predicting computer program (TFsearch and AliBaba). The sequence analysis failed to identify a TATA-box (Fig. 3D), but we identified several putative E2F binding sites that are also known to regulate H19 expression (3). AML-1a, GATA-1/2, N-Myc, and other putative binding sites were also detected (Fig. 3D). Elsewhere, we found that the 3' limit was located about 70 kb upstream of H19 (Fig. 3E and F). Indeed, transcriptional activity has been detected on total cDNA with amplifications IgfA and IgfB but not with amplifications IgfC to IgfG. 3'RACE-PCR experiments mapped the 3' end at bp -72044 relative to the H19 transcription start site (Fig. 3G and I). We confirmed these data by PCR using primers located on each side of the 3' end, and we obtained amplification only with the sense and antisense 1 primer pair, but no fragment could be detected with the sense and antisense 2 primer pair (see Fig. 3H). Thus, antisense H19 transcription was detected throughout a large portion of the H19/IGF2 locus and has the potential to produce a single 119.392-kb transcript that we named the 91H RNA.
![]() ![]() View larger version (60K): [in a new window] |
FIG. 3. Determination of the 5' and 3' limits of the antisense H19 transcription. (A) Map of the MRPL23 gene region. The relative positions of the PCR amplifications and strand-specific primer used for reverse transcription (MrpB) are indicated. (B) Table indicates cycle thresholds obtained in real-time PCR for each PCR amplifications on total cDNA (random priming of the RT) or after strand-specific reverse transcription with MrpB primer. gDNA is used as a positive amplification control to assess for the amplification efficiency. (C) Determination of the transcription start site of the 91H gene by 5'RACE-PCR. RACE-PCR using T47D cell cDNA and 91H specific primers resulted in the amplification of an 1,300-bp DNA fragment. This fragment was cloned into TOPO plasmid for sequencing. (D) Nucleotide sequence of the 5' region of the 91H transcript. The transcription start site is indicated by a black arrow. Putative transcription factor binding sites identified by computer analysis are underlined. (E) Map of the H19/IGF2 locus. Boxes indicate PCR amplifications used to determine the 3' limits of the antisense transcription. Successful amplifications are depicted by full boxes, whereas abortive amplifications are shown as open boxes. (F) A table indicates the cycle thresholds (Ct) obtained for each PCR amplification. (G) Determination of the 3' end of the 91H RNA by 3'RACE-PCR. RACE-PCR using T47D cell cDNA resulted in the amplification of an 1,750 bp-DNA fragment. This fragment was cloned into the TOPO plasmid for sequencing. (H) Confirmation of the 3' end of 91H. RT-PCRs were performed with the primer pairs depicted in panel I: sense plus antisense 1 or sense plus antisense 2. (I) Nucleotide sequence of the 3' region of the 91H transcript. The H19 and IGF2 genes are indicated in boldface, and the 3' limit of 91H is pointed out by a black line. Positions relative to the H19 transcription start site (+1) are indicated. Primer sequences can be seen in Table S1 in the supplemental material.
|
|
View this table: [in a new window] |
TABLE 1. Cellular localization of human 91Ha
|
![]() View larger version (22K): [in a new window] |
FIG. 4. 91H expression is upregulated in breast cancer. (A to C) 91H RNA levels (A), H19 mRNA levels (B), and IGF2 mRNA levels (C) in human breast cells. RNA levels were determined by real-time RT-PCR on total RNA samples from three cancer cell lines (MCF-7, T47D, and BT20) and NBEC. (D to E) 91H RNA levels (D), H19 mRNA levels (E), and IGF2 mRNA levels (F) in human normal and cancer breast tissues. Expressions were normalized on RPLP0 gene expression. H19 RNA sense levels were quantified by using primers designed over exon/exon junctions, and 91H was detected with the H primer pair previously described.
|
91H stabilization in human cancer cells. One possibility to explain 91H RNA accumulation in cancer cells may be that it results from a stabilization of the 91H transcripts in these cells. To determine the stability of 91H RNA, inhibition of total cellular transcription by actinomycin D was used. MYC was chosen as a control for short half-life mRNA (Fig. 5B and D). Figure 5A and C shows the results from a quantitative RT-PCR analysis in normal and cancer cells treated from 30 min to 4 h with actinomycin D demonstrating that 91H RNA half-lives were dramatically increased in T47D cancer cells (>4 h) compared to normal cells (about 30 min). The same results were obtained in BT20 cells (data not shown). Therefore, 91H RNA accumulation in cancer cells likely results from transcript stabilization.
![]() View larger version (14K): [in a new window] |
FIG. 5. 91H RNA stability in cancer and normal breast cells. Cells were treated with 5 µg of actinomycin D/ml, and RNA levels were determined by real-time RT-PCR at the indicated times. cMYC was used as a positive control for actinomycin D treatment since this mRNA is known to be very unstable. 91H (A) and cMYC (B) RNA levels were determined in NBEC. 91H (C) and cMYC (D) RNA levels were determined in T47D cancer cells.
|
![]() View larger version (46K): [in a new window] |
FIG. 6. Allelic expression analysis of the human 91H RNA. Strand-specific reverse transcription was performed on total RNA from T47D cells with a reverse primer for H19-specific RT and a forward primer for 91H-specific RT (see the Materials and Methods). PCR were performed with primers located on each side of the AluI polymorphic restriction site as indicated in the figure (HP1 and HP2). The PCR products were then digested with AluI and separated into agarose gel. Note that, as expected, the size of the H19 amplification fragment is smaller than those obtained for 91H RNA amplification because of the removal of the 81-bp H19 intron 5 in the H19 RNA. ND, not digested; D, AluI digestion; -RT, amplifications on control reactions made without reverse transcription. gDNA, control amplification on T47D cell gDNA.
|
![]() View larger version (38K): [in a new window] |
FIG. 7. Allelic expression of the murine 91H RNA. (A) Map of H19 region indicating positions of PCR primer pairs used in real-time PCR (mC1 and mC2; mF1 and mF2). (B) RT using a H19-sense-intron 4 primer was performed on total RNA extracted from 5-day-old mouse liver issued from SDP711XM.m.dom. (left panel) or M.m.dom.XSDP711 hybrid mice (right panel). A 301-bp PCR fragment was then obtained by using the mC primer pair and digested (D) or undigested (ND) with BglI at an M. spretus-specific restriction site. (C) Reverse transcription using an mF1 sense primer was performed on the samples described above. A 344-bp PCR fragment was then obtained using the mF primer pair and cut by XbaI at a M. musculus domesticus-specific restriction site.
|
It was then of interest to determine whether, as observed in humans, the mouse 91H RNA was transcribed preferentially from one parental allele. We used a hybrid mouse (Mus musculus domesticus female x SDP711 or the reverse cross) to provide allele-specific restriction sites. SDP711 is a congenic mouse strain that has the distal part of chromosome 7 from Mus spretus origin). We first assessed 91H RNA allelic expression in H19 exon 5. In agreement with data obtained in humans, an M. spretus-specific BglI restriction site (Fig. 7A) allowed us to show that, in 5-day-old mouse liver, 91H RNA is expressed almost exclusively from the maternal allele (Fig. 7B). However, using an XbaI M. domesticus-specific restriction site located between the two endodermic enhancers, we show that the 91H RNA is expressed on both alleles in postnatal liver (Fig. 7C) and heart (data not shown). Elsewhere, as the regions located near the ICR (mA and mB) were recently found to be paternally transcribed (9, 42), we also determined the allelic expression in this region. While, within the ICR, some expression can be detected on the paternal allele, 91H RNA expression is only detected on the maternal allele in the intergenic region located between the ICR and the H19 gene (see Fig. S5 in the supplemental material).
In C2C12 myoblasts, 91H RNA is upregulated upon induction of myoblastic cell differentiation (see Fig. S4C in the supplemental material). 91H RNA levels are always much lower than those of the H19 RNAs. They are 4 to 5 orders of magnitude less abundant than H19 RNA levels and 10 to 100 times less abundant than H19 precursor RNAs. In mouse myoblastic C2C12 cells, the half-life of 91H was found to be very short since 91H RNA is largely depleted after 15 min of actinomycin D treatment, whereas H19 RNA shows no significant changes after 4 h of treatment (data not shown). Therefore, because of its short half-life and even if the 91H RNA levels are rather weak, 91H transcriptional activity may be important. Finally, in mouse proliferating and differentiated myoblastic C2C12 cells, we confirmed the nuclear localization of 91H, whereas the H19 RNA was mainly cytoplasmic (Table 2).
|
View this table: [in a new window] |
TABLE 2. Cellular localization of mouse H19 in proliferative and differentiated C2C12 myoblastic cellsa
|
50% reduction of IGF2 expression after si91H treatment (data not shown). Combining RT-PCR performed on cDNA from siRNA-treated cells and ApaI restriction, we showed that both IGF2 parental alleles are expressed (Fig. 8E). This result may be explained by a loss of IGF2 imprinting on the maternal allele. However, since IGF2 mRNA levels are weak in these cells and because they are decreased in siRNA-treated cells, we suspected that it may rather reflect a reduction of IGF2 expression on the nonimprinted paternal allele, thus uncovering a transcriptional leakage observed on the maternal allele. Using a previously described allele-specific quantitative RT-PCR technique (51), we demonstrated that downregulation of the 91H RNA by siRNA leads to a decrease of IGF2 expression exclusively on the paternal ApaI allele (Fig. 8F). Since this allele is the one that is preferentially expressed in PC3 cells, we can thus reasonably assume that it corresponds to the paternal allele.
![]() View larger version (37K): [in a new window] |
FIG. 8. 91H silencing by RNA interference reduces IGF2 expression on the paternal allele. (A) Map of H19 region indicating the positions of siRNA sequences and the PCR primer pairs used in real-time PCR. T47D cells were transfected with siRNA targeting either the green fluorescent protein (ctrl) or the 91H RNA (si91H 1 and 2 alone or in combination). 91H (B), H19 (C), and IGF2 (D) RNA levels were determined by real-time RT-PCR 24 h after transfection. The results were normalized to RPLP0 expression levels. (E) PC3 cells were transfected with si91H, and the imprinting status of the IGF2 alleles was examined using a RT-PCR amplification digested (D) or undigested (ND) with an ApaI polymorphic restriction site available in this cell line. (F) The IGF2 expression levels were determined on each parental allele using an allele-specific RT-PCR amplification method (51) from PC3 cells transfected with control siRNA or with si91H. (G) Model for 91H trans effect on IGF2 expression. Two sets of enhancers (HUC and ENH sequences) regulate IGF2 expression. Both would be required for full expression of the gene. In addition, the two IGF2 alleles would be competing for a common limited stock of regulatory elements (methylation/acetylation/transcription?). On the maternal allele, 91H would block the locus and prevent the HUC sequences from interaction with any regulatory factors. These would be then directed on the paternal allele, in the HUC region, and/or in the IGF2 promoter region and would cooperate with the cis endodermic enhancers, resulting in IGF2 enhanced expression. Upon siRNA treatment and 91H knockdown, the competition would be lost and, as a consequence, both IGF2 alleles would be now able to interact with the regulatory factors with similar efficiencies. Because these factors are in a limiting amount, a part would be depleted from the paternal IGF2 allele to interact with the maternal allele. This would lead to a paternal IGF2 expression decrease but unchanged maternal IGF2 expression because of the absence of functional ENH sequences. Arrows indicates positive regulations, whereas lines with bars correspond to inhibitions.
|
|
|
|---|
Interestingly, 91H RNA shows an increased stability in cancer cells, leading to its accumulation in cancer cell lines and breast cancer tissues. This transcript stabilization is reminiscent of the previously described stabilization of the H19 ncRNA, which is responsible for its accumulation during in vitro muscle cell differentiation (29).
We noticed a preferential 91H expression on the maternal allele in both mice and humans within the H19 gene region (Fig. 6 and 7) and between the ICR and the H19 cap site (Fig. S5B in the supplemental material). However, we detected a biallelic transcription in the 3' endodermic enhancers region (Fig. 7C) and in the ICR in mouse (Fig. S5A in the supplemental material) according to the results obtained by Drewell et al. (9) and Schoenfelder et al. (42), who demonstrated biallelic sense and antisense transcription in the mouse H19 ICR. These data show that, in the mouse, some regions, and particularly regulatory regions, probably generate transcription independently of 91H. However, in humans, using a HhaI polymorphic site, we only detected monoallelic antisense transcription in the ICR (Fig. S7 in the supplemental material).
The developmental expression patterns of 91H are very similar to those of the H19 RNA in murine liver, and both RNAs are upregulated upon induction of myoblastic cell differentiation. This suggests that both transcripts are subject to a common regulation of the locus which allows gene transcription in a time-restricted manner.
To address the functional role of the antisense transcript, we knocked down its expression by RNA interference assays in human cells. The results indicated that, whereas 91H silencing does not modify H19 imprinted status and slightly affects H19 expression, it strongly reduces the overall IGF2 mRNA levels. It was therefore somewhat surprising that semiquantitative RT-PCR detected transcriptional activity on both IGF2 alleles in siRNA-treated cells. To reconcile these data, we assumed that this IGF2 reexpression from the maternal allele may rather reflect a transcriptional leakage that was underestimated in untreated cells because of a high expression level on the paternal allele. Indeed, allele-specific quantitative RT-PCR demonstrated that downregulation of the 91H RNA leads to a decrease of IGF2 expression exclusively on the paternal allele. Therefore, we conclude that, in humans, 91H RNA produced from the maternal chromosome does not affect genomic imprinting but is involved in the maintenance of IGF2 gene expression in trans on the paternal allele. Importantly, two distinct siRNA probes, targeting different regions of the 91H transcript and used in separate experiments have identical effects on Igf2 mRNA levels. Therefore, the possibility that the siRNA probes interfere with known cis-acting regulatory DNA elements is very unlikely. IGF2 upregulation through a trans effect by ectopic expression of H19 antisense transcripts has been suggested (53). The fact that 91H knockdown only slightly affects H19 gene expression from the maternal allele seems to indicate that the 91H effect is not mediated through the ICR/CTCF complex. However, other candidates such as enhancer DNA regions may be considered for interacting with the antisense RNA. Indeed, 91H encompasses a region including a previously described nuclease hypersensitive site (approximately 30 kb upstream of the H19 gene) that can be found on both parental alleles (2) and that have enhancer functions in both imprinted and nonimprinted tissues (6). In addition, the biallelically transcribed H19 upstream conserved (HUC) sequences are regulatory elements located just upstream of the ICR that act as strong mesodermal enhancers and are supposed to activate IGF2 expression on the paternal allele (9).
The 91H trans effect may reflect a competition model that explains a communication between both parental chromosomes that we evoked in a previous work (11). Indeed, deletions of the active maternal allele of the H19 gene reduce the paternal IGF2 DMR2 methylation level and therefore probably have both cis and trans effects on IGF2 gene expression. These results are consistent with the present study and may be due to 91H transcript disruption. Elsewhere, it has been demonstrated that ncRNAs are able to guide the epigenetic regulation of genes. For the human HOX loci, their expression demarcates broad chromosomal domains of differential histone methylation and RNA polymerase accessibility (37). In particular, the HOTAIR ncRNA located in a regulatory boundary in the HOXC cluster represses transcription in trans across 40 kb of the HOXD locus. This ncRNA binds to an epigenetic regulatory complex and changes the methylation patterns of the target locus. We then proposed a model to explain the 91H trans effect on IGF2 expression in which cooperation between two sets of enhancer sequences would be required for full expression of the gene. In addition, both parental alleles would be competing for a common limiting pool of regulatory factors. In this model, 91H would be essentially involved in modifying the accessibility of HUC sequences to regulating proteins. Indeed, maternal expression of 91H would preferentially direct regulatory elements on the paternal allele within the HUC region. This region would cooperate with ENH sequences, resulting in paternal IGF2 enhanced expression (Fig. 8G). We cannot exclude the intervention of other enhancer sequences that may be responsible for the trans effect of 91H RNA on IGF2 expression (18).
Elsewhere, 91H may attract repressive chromatin modifications to the maternal allele by trapping factors responsible for DNA or histone modifications that become available for the paternal allele when 91H is disrupted, resulting in IGF2 expression decrease. Yu et al. (52) recently showed that many tumor suppressor genes have nearby antisense RNAs. These authors focused on the p15 antisense RNA, which is able to induce p15 gene silencing through heterochromatin formation (increase in dimethylation of H3K9 and decrease in dimethylation of H3K4). Thus, it will be useful to analyze histone methylation and/or acetylation patterns of the paternal IGF2 region in the presence or in the absence of the 91H RNA.
Other authors have already evoked this type of trans regulation within an imprinted locus. Takeda et al. (46) examined the callipyge phenotype (postnatal muscular hypertrophy of shape) caused by a mutation in the DLK1-GTL2 locus. The mutation and its "polar overdominance" mode of inheritance imply that only the heterozygous animals that have received the paternal mutation display the phenotype. The mutation, when paternally inherited, increases DLK1 and PEG11 protein levels, which results in the callipyge phenotype. However, in the case of a double maternal and paternal inherited mutation, it was suggested that the overexpression of ncRNAs from the maternal allele could inhibit in trans the translation of DLK1 and PEG11 genes, resulting in a normal phenotype. This trans inhibition have been formally demonstrated for PEG11 through a maternal miRNA. For DLK1, these authors propose that the posttranscriptional inhibition may be mediated by a long ncRNA which could, for instance, compete toward a trans factor responsible for DLK1 RNA localization or translation. Our model shares interesting similarities, although the trans regulation affecting DLK1 protein levels is different from those affecting IGF2 transcripts accumulation after 91H knockdown.
Interestingly, and in agreement with the role of 91H displayed in PC3 cells, we noticed an inverse correlation between 91H and IGF2 transcription in normal and cancer mammary cell lines. However, in breast cancer tissues, whereas 91H is always overexpressed, IGF2 and H19 expression profiles are more heterogeneous; this is probably due to various co- or posttranscriptional regulations. These data provide a link between antisense transcription and cancer.
Finally, 91H RNA discovery leads us to reconsider some previous results. Indeed, a set of studies in the mouse described deletions of sequences within the H19/IGF2 locus that did not affect H19 and IGF2 expression: deletion of the H19 promoter and structural gene (41), deletion of a G-rich repetitive element located immediately 3' to the ICR (48), deletion of a conserved direct tandem repeat 1 kb upstream of the ICR (23), and deletion of a 4.2-kb domain between the H19 transcription unit and the enhancers (49). We presume that these deletions did not impede 91H transcription and generated a transcript that could still assume its regulatory function.
In conclusion, the antisense transcript reported here adds further complexity to the cluster of imprinted genes in the human imprinted 11p15.5 region. 91H RNA is involved in the control of IGF2 expression in trans, suggesting that, beyond imprinting, other sophisticated mechanisms involved in the control of gene expression arise at this locus. Cooperation between individual regulatory elements (i.e., enhancers, silencers, antisense RNA, DNA, and histone modifications) seems to be required for fine-tuning H19 and IGF2 gene expression.
J.-J.C., E.A., and T.D. hold grants from the Fédération des Groupements des Entreprises Françaises dans la Lutte contre le Cancer, the Comités du Nord et de l'Aisne de la Ligue Nationale contre le Cancer, and the Association pour la Recherche sur le Cancer (ARC). N.B. is the recipient of fellowships from the Ministère de l'Education Nationale et de la Recherche, the ARC, and the Institut National du Cancer. This study was also supported by grants from the ARC (grant 4868), the GIS Longévité (contract GISLO401), the Fond National de la Science (ACI Jeune Chercheur), and the Agence Nationale de La Recherche (ANR-07-BLAN-0052-02 -T.F.-) and by funds from the Centre National de la Recherche Scientifique.
We dedicate this report to the memory of our friend Jean Coll.
Published ahead of print on 15 September 2008. ![]()
Supplemental material for this article may be found at http://mcb.asm.org/. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»