Next Article in Journal
Activation of the miR-371/372/373 miRNA Cluster Enhances Oncogenicity and Drug Resistance in Oral Carcinoma Cells
Next Article in Special Issue
NLRP3 Inflammasome: A New Pharmacological Target for Reducing Testicular Damage Associated with Varicocele
Previous Article in Journal
Space-Dependent Glia–Neuron Interplay in the Hippocampus of Transgenic Models of β-Amyloid Deposition
Previous Article in Special Issue
Role of Nuclear Claudin-4 in Renal Cell Carcinoma
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 24
10.3390/ijms21249433
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Many Different LINE-1 Retroelements Are Activated in Bladder Cancer

by
Patcharawalai Whongsiri
1,
Wolfgang Goering
2,
Tobias Lautwein
3,
Christiane Hader
1,
Günter Niegisch
1,
Karl Köhrer
3,
Michèle J. Hoffmann
1 and
Wolfgang A. Schulz
1,*
1
Department of Urology, Medical Faculty, Heinrich Heine University, 40225 Dusseldorf, Germany
2
Institute of Pathology, Medical Faculty, Heinrich Heine University, 40225 Dusseldorf, Germany
3
Biological and Medical Research Center (BMFZ), Heinrich Heine University, 40225 Dusseldorf, Germany
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(24), 9433; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21249433
Submission received: 27 October 2020 / Revised: 4 December 2020 / Accepted: 7 December 2020 / Published: 11 December 2020
(This article belongs to the Special Issue Molecular Research on Urology 2.0)
Browse Figures
Versions Notes

Abstract

:
Human genomes contain about 100,000 LINE-1 (L1) retroelements, of which more than 100 are intact. L1s are normally tightly controlled by epigenetic mechanisms, which often fail in cancer. In bladder urothelial carcinoma (UC), particularly, L1s become DNA-hypomethylated, expressed and contribute to genomic instability and tumor growth. It is, however, unknown which individual L1s are activated. Following RNA-immunoprecipitation with a L1-specific antibody, third generation nanopore sequencing detected transcripts of 90 individual elements in the VM-Cub-1 UC line with high overall L1 expression. In total, 10 L1s accounted for >60% of the reads. Analysis of five specific L1s by RT-qPCR revealed generally increased expression in UC tissues and cell lines over normal controls, but variable expression among tumor cell lines from bladder, prostate and testicular cancer. Chromatin immunoprecipitation demonstrated active histone marks at L1 sequences with increased expression in VM-Cub-1, but not in a different UC cell line with low L1 expression. We conclude that many L1 elements are epigenetically activated in bladder cancer in a varied pattern. Our findings indicate that expression of individual L1s is highly heterogeneous between and among cancer types.

1. Introduction

Roughly one-sixth of the human genomes are derived from long interspersed element 1 (LINE-1, L1) retroelements. Most L1s are truncated or heavily mutated, but more than one hundred are intact and could in principle be capable of retrotransposition. Intact elements contain an internal promoter at their 5′-end and two open reading frames encoding the RNA-binding protein ORF1p and the reverse transcriptase/endonuclease ORF2p [1,2]. Their transcription is however restrained by epigenetic mechanisms including DNA methylation and heterochromatic histone modifications like H3K9me3 [2,3,4,5]. Together with additional mechanisms acting at the posttranscriptional level and at the genomic insertion step, L1 expression and retrotransposition events are limited to short periods during germ cell development and embryogenesis. However, these controls often fail during cancer development, evidenced by L1 DNA hypomethylation, increased expression of transcripts and encoded proteins and actual retrotransposition events [1,6,7,8]. Reactivation of L1 in cancer may contribute to tumor development and progression in various ways, e.g., by increasing genomic instability, disrupting tumor suppressor genes, altering transcription of adjacent genes and stabilizing telomeres to counter senescence [1,9].
An important issue in current research is to elucidate which specific intact L1 elements are active and expressed during normal development and in cancer. In particular, some studies have hinted that individual “master” elements predominate, perhaps in a cell type-specific manner [10,11]. Current knowledge on individual L1 expression in tumor cells and tissues is based largely on two kinds of approaches. Extensive whole genome sequencing data combined with refined bioinformatic analyses have allowed to delineate the frequency of retrotransposition events in various cancers and their consequences for the tumor genome; in many cases, moreover, the respective source elements could be identified, e.g., refs. [6,12,13]. These studies suggest that a limited number of different L1 elements contribute to retrotransposition in each tumor type, with some elements exhibiting particularly high activity, e.g., a L1 located at the TTC28 locus at chromosome 22q in colorectal cancers. This conclusion is further supported by studies based on comprehensive RNA analyses, with or without enrichment of L1 transcripts by various techniques, e.g., refs. [11,14,15]. Both types of studies thus suggest that the extent of L1 reactivation varies in between and among cancer types and that the individual elements being activated vary likewise. Specifically, Deininger et al. [15] reported that in commonly used cell lines like HeLa and HEK293, transcripts originating from a large number of elements can be detected, but that this repertoire is dominated by fewer than 10 elements with high expression. To date, however, no comprehensive delineation of the expressed L1 repertoire is available for any major cancer type.
Decreased DNA methylation at L1 promoters is particularly prevalent in urothelial carcinoma (UC), the most common cancer of the urinary bladder, and is associated with increased overall expression of full-length L1 elements [16,17,18]. Moreover, differences in the extent of active and repressive histone modifications at full-length L1 (flL1) sequences globally corresponded to differences in the levels of L1 mRNA and ORF1p expression in UC cell lines [19]. However, it is not known yet which individual L1 elements are expressed in this cancer type. Analysis of L1 hypomethylation at individual elements suggests a high variability [17]. Likewise, several different elements have been observed at retrotransposition sites in UC tissues [13].
Here, we defined the individual L1 elements that are activated and expressed in the UC cell line VM-Cub-1 with high overall L1 mRNA and ORF1p expression using a novel approach. We have previously shown that L1 expression in that cell line is functionally relevant contributing to cell proliferation and escape from senescence [18]. Since L1 ORF1p associates preferentially with full-length L1 transcripts [14], we precipitated ORF1p-associated RNA using a highly specific antibody against ORF1p that has become available recently [20]. The precipitated RNA was analyzed by nanopore long-read sequencing, which allows to directly identify individual L1 elements despite their close homology. This new approach also avoids contamination of sequencing data by L1 sequences from gene introns or fusion transcripts. In order to obtain a broader picture of the expression pattern of individual L1s, we measured the expression of five individual elements highlighted by nanopore sequencing using specific RT-qPCR assays across a large panel of cell lines as well as in UC and normal bladder tissues. Finally, we investigated the epigenetic status of three L1s by chromatin immunoprecipitation in two cell lines with high and low L1 expression, respectively.

2. Results

2.1. Delineation of L1 Expression by RNA Immunoprecipitation and Nanopore Sequencing

In previous work, we had observed that L1 mRNA and ORF1p expression vary strongly among UC cell lines [16,18]. As illustrated in Figure 1A, some cell lines, like VM-Cub-1 and BFTC-905, express high levels of ORF1p, albeit still less than some embryonal carcinoma cell lines like NCCIT. Other UC cell lines, like 5637 or UM-UC-3, express very little ORF1p, and the protein is undetectable in non-transformed urothelial cells like HBLAK.
To comprehensively identify all expressed individual L1s, we applied the strategy illustrated in Figure 1B to VM-Cub-1 cells. Initially, as ORF1p is known to bind co-translationally to the L1 transcripts that encode it [14], we optimized a protocol to immunoprecipitate all RNA species associated with the ORF1p protein using the highly specific antibody described by Rodic et al. [20]. Indeed, transcripts of full-length LINEs were several-fold enriched in VM-Cub-1 following immunoprecipitation with the ORF1p antibody, whereas no enrichment was achieved and only very low amounts of RNA were obtained if the same protocol was used on ORF1p low-expressing cells, indicating the specificity of the approach.
The immunoprecipitated RNA from VM-Cub-1 cells was then subjected to third generation nanopore sequencing. In two experiments using two different passages of VM-Cub-1 cells, 1843 and 2291 reads, respectively, could be assigned to full-length L1 elements of the evolutionarily young and active L1Hs family using stringent criteria (mapping quality ≥ 20). Exemplary alignments are illustrated in Figure 1C. In each experiment, around 100 of the 146 elements listed in L1Base (http://l1base.charite.de) were detected; following manual curation, 90 elements yielded unequivocal reads in both experiments (Table S1). The top 10 in each case comprised more than 60% (63.5% and 61.7%) of the reads, whereas the bottom 50 element made up only about 5% (Figure 1D). Notably, the top eight most frequently detected elements were the same between the experiments, although their order differed, likewise those ranking 9th–11th (Table 1).
Along with L1 transcripts, 3562 genes were at least 2-fold enriched by the RIP procedure, and 1024 at least 5-fold, compared to a total RNA-Seq analysis of VM-Cub-1 cells (Supporting data sets 1 and 2). According to a gene ontology analysis using the GOrilla database [21], genes related to positive regulation of defense response (GO:0031349, p = 0.000455), SCF-dependent proteasomal ubiquitin-dependent protein catabolic process (GO:0031146, p = 0.00077) and ubiquitin-dependent protein catabolic process (GO:0006511, p = 0.000879) were significantly overrepresented among the > 5-fold enriched genes.
L1 transcripts are known to accumulate in organelles associated with RNA turnover, like P-bodies and stress granules [22,23,24]. However, no significant overlap between the enriched gene set identified in our experiments and transcripts enriched in P-bodies [25] or stress granules [26], respectively, were observed. It should be considered, though, that these latter analyses were performed in different cell types.

2.2. Analysis of L1 Expression by qRT-PCR in Cancer Cell Lines and Tissues

For 5 of the top 10 L1s, namely UID-33, UID-59, UID-60, UID-66 and UID-108 (numbered according to L1Base), specific qRT-PCR assays could be designed. Using these assays, their expression was measured across a broader range of UC cell lines and normal controls (Figure 2) as well as in a series of urothelial bladder carcinoma tissues and morphologically normal tissues from cancerous tissues (Figure 3). The same series of samples was also investigated for expression of full-length L1Hs elements in toto (flL1) using the previously described 5′-LINE assay [27], which covers a few hundred full-length elements.
In keeping with our previous observations, flL1 expression was generally higher in bladder cancer cell lines than in benign controls (Figure 2a). As expected, all five individual elements were strongly expressed in VM-Cub-1 and were only weakly or not at all expressed in normal urothelial cells and in immortalized normal urothelial cell lines (HBLAK and TERT-NHUC), whereas their expression across other UC cell lines was very heterogeneous (Figure 2b–f). For instance, strong expression of UID-108 was restricted to a few UC cell lines (Figure 2f) and UID-33 was rather selectively expressed in a subset of the cancer cell lines (Figure 2b). The other elements were more widely expressed. It is instructive to compare the pattern of L1 expression among individual cell lines. Thus, as in VM-Cub-1; the five individual L1s were usually strongly expressed in other cell lines with high overall flL1 expression like BFTC-905 and SCaBER. In contrast, in 5637, with overall low flL1 and ORF1p expression (cf. Figure 1a), all individual L1s were expressed at low levels only, whereas in UM-UC-3, likewise low in flL1 and ORF1p, two elements were comparably strongly expressed.
In UC cancer tissues, flL1 expression was overall increased compared to benign bladder tissue samples, but with borderline statistical difference (p = 0.06 according to Mann–Whitney U-test) (Figure 3a). As in our previous study [16], the difference became statistically significant (p = 0.023) following adjustment for background transcription via an assay targeting the L1 3′-region (Figure S1). All individual elements were overexpressed in cancer tissues in a highly statistically significant fashion, except for UID-108 (Figure 3b–f). Notably, UID-33 and UID-66 expression were generally low in benign tissues, and high expression was restricted to cancer tissues, albeit not observed in every tumor (Figure 3b). The tumors with particularly high UID-33 expression were often advanced (i.e., stage pT3 or pT4); however, no statistically significant association of higher expression with tumor stage was observed for any other individual element (Figure S2). Expression of UID-33 and UID-66 was also higher in tumors with lymph node metastases, albeit with borderline statistical significance (Figure S3). No association with patient age or gender was observed.
Additionally, we investigated the expression of the individual L1 elements in testicular germ cell tumor (embryonal carcinoma) and prostate cancer cell lines in comparison to three bladder cancer cell lines with high, intermediate and low transcript levels. As among the UC cell lines, expression patterns were diverse, but all L1s could be detected in each cell line (Figure 4). Remarkably, several elements appeared most strongly expressed in the embryonal carcinoma cell lines, notably UID-108 (Figure 4e), whereas UID-33 expression was strongest in the BFTC-905 bladder cancer cell line (Figure 4a). Prostate cancer cell lines, which are all derived from metastatic tumors, also presented relatively high levels for most L1s (Figure 4).

2.3. Analysis of Histone Modifications at Individual L1s

To investigate whether histone modifications on individual L1 elements corresponded to expression ChIP-qPCR assay were carried out in the two bladder cancer cell lines VM-Cub-1 and 5637 with high and low L1 expression, respectively (Figure 5a). ChIP assay specificity was controlled by qRT-PCR on GAPDH as a constitutively expressed gene and CTCFL as a developmentally regulated gene weakly expressed in bladder cancer cells [28]. As expected, the active histone marks H3K18Ac and H3K4me3, but not the repressive histone marks H3K9me3 and H3K27me3 were enriched at the GAPDH promoter (Figure 5b). Further, as expected, H3K27me3 dominated on the CTCFL promoter in both cell lines, with lower levels of H3K9me3 and active marks (Figure 5c). The heterochromatic gene GRM6 was characterized by a strong H3K9me3 signal (Figure 5d). These controls validate the ChIP-qPCR assay in the present experiments. On L1-5′ promoters overall, we observed enrichment of active histone marks, especially H3K18Ac, in VM-Cub-1 compared to 5637, whereas H3K9me3 especially showed the converse pattern (Figure 5e). These findings are in accord with our previous report [19]. Specific assays could be designed for proximal upstream sequences of three individual L1s, UID-33, UID-66 and UID-108. As for LINE-1 promoters globally, the active histone marks, especially H3K18Ac, were clearly enriched at the promoters of these individual elements in VM-Cub-1 cells, but less so or not at all in 5637, fitting their higher expression in VM-Cub-1 (Figure 5f–h).

3. Discussion

To investigate the repertoire of individual L1 elements in bladder cancer, we used a new technical approach. The availability of a highly specific ORF1p antibody [20] allowed specific immunoprecipitation and third generation sequencing by nanopore technology, which yields long reads, permitted unambiguous alignment to individual L1 elements despite their extensive homology. Using this procedure, we detected transcripts from many loci, although some loci were more prominently expressed. This ‘long-tailed’ frequency distribution of individual elements from the RIP analysis of VM-Cub-1 UC cells is similar to that observed by Deininger et al. [15] in broadly used transformed cell lines like HeLa and HEK293, although their distribution was steeper.
Similarly, the recent pan-cancer analysis by Rodriguez-Martin identified 10 individual L1s as source elements for roughly two-thirds of all retrotransposition events, with the remaining third contributed by a broad range of others. However, retrotransposition events in the 23 bladder cancer samples could be assigned to only up to three source elements per sample and were only moderately frequent, despite the evidence for strong L1 activation in bladder cancer. This difference could be due to several factors. First, individual elements may differ in their actual competence for retrotransposition, as demonstrated by their different efficiency following expression from a vector in a model system [11]. Second, in addition to epigenetic repression at the transcriptional level, a variety of mechanisms act post-transcriptionally to restrain retroelement activity [2]. Such mechanisms may be active in UC cells since even ectopic overexpression of a L1 from a strong viral promoter in bladder cancer cell lines yielded only moderate increases in transcript levels [18]. Third, different experimental approaches will catch L1s at different phases of their expression. Thus, our ORF1p-based RIP would capture L1s after translation of their transcripts, when they could move on towards retrotransposition, but could also be sequestrated and degraded [22]. Notably, the co-precipitated transcripts in our experiments were not enriched for P-granule or stress granule associated transcripts, with the reservation that data from other cell types [25,26] had to be used for comparison. Fourth, retrotransposition is counteracted at the insertion step by various mechanisms, including checkpoint activation and cell death [1]. Many of these mechanisms involve the tumor suppressor p53, which is typically inactivated in advanced bladder cancers [29], but other mechanisms function independently of p53 [30,31]. Finally, there is evidence that the pattern of L1 expression and consequently retrotransposition may change during tumor progression, e.g., in prostate cancer [6]. Therefore, expression of individual elements in a tumor at any particular time may not necessarily indicate which elements will actually successfully retrotranspose into the genome.
Even with these considerations, it remains remarkable how little overlap we observed between the elements reported as frequent sources for retrotransposition and those strongly expressed in VM-Cub-1 and other bladder cancer cell lines. In fact, only one of the top elements reported by Deininger et al. [15] in HeLa and HEK-293 cells ranked also in our top 11 and likewise, only one of the frequent retrotransposition source elements across cancers identified by Rodriguez-Martin et al. [13]. The one exception, interestingly, was UID-135 from the TTC28 locus, which has been consistently reported as the most active element in colorectal and other cancers by several studies [6,11,12]. However, tellingly, in our analysis this element ranked only around 10th. We therefore feel that the results from our study and the literature are best explained by the assumption that a relatively large number of LINEs rather than one “master” element are expressed in many tumor cells, and that the L1 repertoire differs strongly between as well as within tumor types. The methodology presented in the present study, combining RIP and nanopore sequencing, should provide a straightforward approach to corroborate this conclusion, e.g., using embryonal carcinoma and prostate carcinoma cell lines, which appear to have a broad range of expressed L1s according to the results obtained by qRT-PCR.
This conclusion is underlined by the qRT-PCR analysis of individual L1 elements which demonstrated highly variable expression among UC cell lines, but also in cell lines from testicular germ cell cancers, which are known to express overall high levels of retroelements [32,33], and from prostate cancer, where expression appears to be stronger in advanced stage tumors and cell lines derived from them [27,34]. Of all L1s, UID-33 appeared overall most characteristic of bladder cancer, although it was not expressed strongly throughout all cell lines and tissues. Conversely, strong UID-108 expression was rather the exception in UC cell lines and tissues but was prominent in cell lines from testicular and prostate cancer.
Due to its high prevalence, L1 hypomethylation has been explored as a biomarker for bladder cancer diagnostics. For instance, one panel of methylation biomarkers for bladder cancer detection includes L1-MET, the L1 element next to the MET oncogene [35]. In a similar fashion, overexpression of L1s might be exploited as a cancer biomarker. However, measurements of overall L1 expression are confounded by background expression from a large number of non-intact L1 sequences from cellular gene promoters and do not distinguish well between tumor and benign tissues. In addition, they are highly sensitive to DNA contamination. These complications could be avoided by using assays for individual elements. Indeed, most qRT-PCR assays for individual L1s in our study significantly distinguished cancerous from normal tissues, promising good specificity for cancer detection. However, the high variability of the expression of each element on its own would limit assay sensitivity. We therefore assume that assays for individual L1s like UID-33 could be useful in bladder cancer diagnostics only as part of RNA panels in combination with other markers. Another issue to be more deeply explored is the association of individual L1 expression with clinical parameters, especially stage and prognosis. Our present study does not point at strong associations of expression of any element with clinical stage and the association with clinical outcome was not investigated. Nevertheless, such associations may be revealed by investigations on larger sample sets with long-term clinical follow-up. Moreover, as our tissue series comes from cystectomies, we could not include tissues from low grade and non-muscle-invasive tumors. L1 hypomethylation is an early event in urothelial carcinogenesis and some evidence hints at increased L1 expression [36], but a systematic investigation of total and individual L1 expression has not yet been performed and should be feasible with the techniques developed in the present study.
It is commonly thought that L1 expression is controlled by epigenetic mechanisms, especially DNA methylation at their moderately CpG-rich internal promoters and likely by histone modifications, especially H3K9me3 for repressed and H3 acetylation for active elements [2,4,37]. To date, however, too few studies are available to allow firm conclusions on how transcription of individual elements is enabled. In this regard, DNA methylation analyses of bladder cancer tissues suggest pronounced inter-tumor variability in the hypomethylation of individual elements [17]. Moreover, using a global assay for full-length L1Hs elements, we had previously observed that histone modification patterns in their 5′-region differed mostly by acetylation and, less clearly, by H3K4me3 methylation between strongly and weakly expressing bladder cancer cell lines [19]. Our present investigation of three individual elements confirms these findings. Thus, active elements appear to be characterized especially by histone acetylation, whereas H3K4me3 methylation is less prominent. Obviously, a larger number of individual elements and cell lines need to be investigated to elucidate the relationship between histone modifications and L1 transcriptional activity more precisely. On a technical note, ChIP analyses of individual L1s are complicated by the fact that they require primers for unique sequences within the first few hundred bp upstream of an element, but these upstream sequences contain repeats or other sequences unsuitable for qPCR for many L1 elements. Moreover, since primers have to be chosen upstream of the repetitive L1 sequence, H3K4me3 methylation at the actual L1 promoters may be underestimated by our analysis.

4. Materials and Methods

4.1. Tissue Samples

Urinary bladder cancer and normal tissues were obtained from cystectomies performed at the Dept. of Urology of the Heinrich Heine University. Overall, 51 cancer tissues, including 3 non muscle-invasive, 12 pT2, 22 pT3 and 14 pT4 cases as well as 31 morphologically normal bladder samples were analyzed. All tumors were high grade. In total 39 (76.5%) of the patients were male and 12 (23.5%) were female. Their ages ranged from 51 to 94 years, with a mean age of 72.4 ± 8.9 years. Of note, the sample set overlaps with that investigated in a previous study on total retroelement expression in bladder cancer [16]. Analysis of tissue samples was permitted by the ethics committee of the medical faculty of the university (study numbers #3836 and #4371, approval dates 12 April 2012 and 22 January 2013) and all patients agreed to the use of the tissues.

4.2. Cell Lines and Cell Culture

All bladder cancer cell lines (UCCs) (253J, 5637, 639-V, 647-V, BFTC-905, HT-1376, J82, MGHU4, RT4, RT-112, SCaBER, SD, SW1710, UM-UC-3, UM-UC-6, VM-Cub-1, T24) were cultured in DMEM GlutaMax (Gibco, Darmstadt, Germany), supplemented with 10% fetal calf serum [16]. BC61 was cultured as described [38]. The cell lines were obtained from the DSMZ (Braunschweig, Germany), except for UM-UC-3, kindly provided by Dr. Grossman, Houston. The telomerase-immortalized TERT-NHUC cell line was kindly provided by M. A. Knowles (Leeds, United Kingdom) and cultured as described. The spontaneously immortalized urothelial cell line HBLAK was cultured as described [39]. Primary urothelial cells cultures (UP) were established from ureters after nephrectomy and were routinely maintained in keratinocyte serum-free medium (KSFM, Gibco, Darmstadt, Germany) supplemented with 12.5 µg/mL bovine pituitary extract and 0.25 ng/mL epidermal growth factor as described [40]. Establishment and use of primary urothelial cells were approved by the ethics committee of the HHU Medical Faculty (study number #1788, approval date 23/4/2001). Prostate and testicular cancer cell lines were cultured as described [27].

4.3. RNA Immunoprecipitation via ORF1p

The RNA immunoprecipitation (RIP) procedure from a published protocol [41] with modifications as described [42] was adapted for L1 transcripts. In brief, 2 × 107 VM-Cub-1 or 5637 cells cultured under standard conditions were harvested by trypsinization, washed with PBS and lysed in RIPA buffer (50 mM Tris, pH 7.6, 150 mM NaCl, 1 mM EDTA, 0.1% sodium dodecylsulfate, 1% NP40, 0.5% DOC, supplemented with protease inhibitor (Sigma Aldrich, Munich, Germany). Following homogenization by douncing and centrifugation for 5 min at 12,000 g, aliquots were taken for protein quantification and RNA purification and the remaining supernatant was incubated with SureBeads Protein G Magnetic Beads (Biorad, Puchheim, Germany) preloaded with anti-ORF1p antibody (clone 4H1, MABC1152, Merck Millipore, Darmstadt, Germany) or an IgG control antibody (20 µg each) overnight at 4 °C. Following washing with RIPA buffer, the bound RNA was purified using acid phenol extraction, precipitated with ethanol and re-dissolved in RNA-free water.

4.4. Nanopore Sequencing and Data Evaluation

Library preparation for Nanopore Sequencing of p-ORF1 immunoprecipitated RNA was carried out using the PCR-cDNA Sequencing Kit (SQK-PCS108, Oxford, UK) from Oxford Nanopore Technologies following manufacturer’s guidelines. Sequencing was carried out on the GridION platform as per manufacturer’s guidelines using R9.4.1 flow cells. Each sample was sequenced on an individual flow cell and base calling was carried out by guppy v1.4.0. Sequencing yielded 6.3 million reads in the first experiment (sample 8VmCub1ORF1) and 7.0 million reads in the second experiment (sample 10VmCub1ORF1). The nanopore sequencing data are available at the SRA database, accession number PRJNA657989.
As a first step in data analysis, reads shorter than 100 bp were removed using the filtlong tool (v0.2.0, Melbourne, Australia). For further analyses, the filtered reads were mapped to a custom reference consisting of all 146 sequences of Human full-length, intact LINE1 Elements [fLL1] downloaded from L1Base (http://l1base.charite.de) [43] using minimap2 (v2.12). Only reads with a MAPQ score of ≥20 were used in subsequent analysis. Coverage analysis was completed using samtools v1.7 (Cambridge, UK). The individual L1 elements identified unambiguously in both runs are listed in Table S1.
To identify transcripts from protein-coding genes co-precipitated during RIP, htseq-count (v0.11, Heidelberg, Germany) was used to identify the highest expressed genes in the ORF1p precipitated RNA samples. Differential expression analysis was carried out via the R-package DESeq2 (v1.22.2, Boston, MA, USA), using a total RNA-Seq analysis of VM-Cub-1 cells as a reference.

4.5. RNA Extraction and Reverse Transcription

Total RNA was extracted from powdered tissues or cell cultures using acid phenol extraction following column purification to minimize DNA contamination as described [27]. Synthesis of complementary DNA was performed using the QuantiTec Reverse Transcription kit (Qiagen, Hilden, Germany) with another DNA removal step to remove any remaining DNA contamination.

4.6. Quantitative Reverse Transcription PCR (qRT-PCR)

qRT PCR was performed as described previously [27] on a Roche LightCycler 96 (Roche, Basel, Switzerland) using the QuantiTect SYBRGreen PCR Kit (Qiagen, Hilden, Germany). All qRT PCR data were adjusted to TATA-box-binding protein (TBP) mRNA (Table S2). For all other transcripts specifically designed primers (Table S2) were employed using the following PCR conditions: initial denaturation step at 95 °C for 15 min, followed by 40 amplification cycles consisting of denaturation at 95 °C for 15 s, annealing at the indicated temperature for 20 s and extension at 72 °C for 30 s. Assay specificity was controlled using UCSC in silico PCR [44] and by Sanger sequencing. Controls for DNA contamination were regularly conducted using mock cDNA samples prepared without addition of reverse transcriptase. All measurements were performed in at least duplicates; assay variance was <10%. Relative expression was calculated by the modified ΔΔCt method [45].

4.7. Chromatin Immunoprecipitation

ChIP was performed according to the kit instruction (ChIP-IT® Express Magnetic Chromatin Immunoprecipitation Kit, 53008, Active Motif, Waterloo, Belgium). In brief, 1.5 × 107 cells were treated with 1% paraformaldehyde to cross-link protein and DNA. Chromatin was extracted and sheared into 300–500 bp fragments (about 1–2 nucleosomes) by sonication. The sheared chromatin (25 µg) was incubated with primary antibody (3 µg) and magnetic Protein G-coupled beads at 4oC overnight. The primary antibodies employed in this study comprised Anti-Histone H3K4me3 (39915, Active Motif, Waterloo, Belgium), Anti-Histone H3K9me3 (ab8898, abcam, Cambridge, UK), Anti-Histone H3K18Ac (ab1191, abcam, Cambridge, UK), Anti-Histone H3 (39763, Active Motif, Waterloo, Belgium), Anti-Histone H3K27me3 (39535, Active Motif, Waterloo, Belgium), normal mouse IgG (sc-2025, Santa Cruz Biotechnology, Santa Cruz, USA) and normal rabbit IgG (cst2729, Cell Signaling Technology, Frankfurt, Germany). Antibody-bound chromatin complexes were precipitated by magnetic separation. Cross-linking reversal of chromatin was performed, and DNA was eluted from the magnetic beads.
Histone modifications were determined on the L1 promoter (5′-LINE-1), the indicated specific L1s, and a set of control genes, CTCFL (also known as BORIS), glutamate metabotropic receptor 6 (GRM6) and GAPDH by SYBRGreen-based qPCR using immunoprecipitated (IP) and input DNA samples as templates as previously described [19]. All samples were amplified by Light Cycler® 96 System (Roche, Mannheim, Germany) using QuantiTect SYBR Green PCR kit (Qiagen, Hilden, Germany) and specific primers (Table S3). The enrichment of chromatin marks at each gene was calculated and expressed as percentage input DNA.

4.8. Statistical Analysis

Two-sample t-test or Mann–Whitney test was used to determine the difference between two independent groups. p values < 0.05 were considered statistically significant.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/21/24/9433/s1.

Author Contributions

Conceptualization, W.A.S., W.G. and P.W.; methodology, T.L., C.H. and M.J.H.; formal analysis, P.W., T.L. and W.A.S.; investigation, P.W., W.G., T.L. and C.H.; resources, G.N., K.K. and M.J.H.; data curation, P.W. and T.L.; writing, W.A.S., P.W. and X.X.; supervision, K.K. and M.J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

Computational support of the Zentrum für Informations- und Medientechnologie, especially the HPC team (High Performance Computing) at the Heinrich-Heine University is gratefully acknowledged. P.W. was supported by a scholarship from the Development and Promotion of Science and Technology talented project (DPST), Thailand.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ChIPChromatin immunoprecipitation
flL1Full-le
Ngth LINE-1
L1LINE-1, long interspersed element 1
ORFOpen reading frame
PBSPhosphate-buffered saline
qRT-PCRQuantitative reverse transcription polymerase chain reaction
RIPRNA immunoprecipitation
UCUrothelial carcinoma

References

  1. Burns, K.H. Transposable elements in cancer. Nat. Rev. Cancer 2017, 17, 415–424. [Google Scholar] [CrossRef] [PubMed]
  2. Goodier, J.L.; Kazazian, H.H., Jr. Retrotransposons Revisited: The Restraint and Rehabilitation of Parasites. Cell 2008, 135, 23–35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Liu, N.; Lee, C.H.; Swigut, T.; Grow, E.; Gu, B.; Bassik, M.C.; Wysocka, J. Selective silencing of euchromatic L1s revealed by genome-wide screens for L1 regulators. Nature 2018, 553, 228–232. [Google Scholar] [CrossRef] [PubMed]
  4. Jönsson, M.E.; Brattås, P.L.; Gustafsson, C.; Petri, R.; Yudovich, D.; Pircs, K.; Verschuere, S.; Madsen, S.; Hansson, J.; Larsson, J.; et al. Activation of neuronal genes via LINE-1 elements upon global DNA demethylation in human neural progenitors. Nat. Commun. 2019, 10, 3182. [Google Scholar] [CrossRef] [Green Version]
  5. Sanchez-Luque, F.J.; Kempen, M.-J.H.; Gerdes, P.; Vargas-Landin, D.B.; Richardson, S.R.; Troskie, R.-L.; Jesuadian, J.S.; Cheetham, S.W.; Carreira, P.E.; Salvador-Palomeque, C.; et al. LINE-1 Evasion of Epigenetic Repression in Humans. Mol. Cell 2019, 75, 590–604.e12. [Google Scholar] [CrossRef]
  6. Tubio, J.M.C.; Li, Y.; Ju, Y.S.; Martincorena, I.; Cooke, S.L.; Tojo, M.; Gundem, G.; Pipinikas, C.P.; Zamora, J.; Raine, K.; et al. Mobile DNA in cancer. Extensive transduction of nonrepetitive DNA mediated by L1 retrotransposition in cancer genomes. Science 2014, 345, 125–134. [Google Scholar] [CrossRef] [Green Version]
  7. Chénais, B. Transposable Elements in Cancer and Other Human Diseases. Curr. Cancer Drug Targets 2015, 15, 227–242. [Google Scholar] [CrossRef]
  8. Ardeljan, D.; Taylor, M.S.; Ting, D.T.; Burns, K.H. The Human Long Interspersed Element-1 Retrotransposon: An Emerging Biomarker of Neoplasia. Clin. Chem. 2017, 63, 816–822. [Google Scholar] [CrossRef] [Green Version]
  9. Bergmann, M. A role of LINE-1 in telomere regulation. Front. Biosci. 2018, 23, 1310–1319. [Google Scholar] [CrossRef] [Green Version]
  10. Rangwala, S.H.; Zhang, L., Jr.; Many, H.H.K. LINE1 elements contribute to the transcriptome of human somatic cells. Genome Biol. 2009, 10, R100. [Google Scholar] [CrossRef] [Green Version]
  11. Philippe, C.; Vargas-Landin, D.B.; Doucet, A.J.; Van Essen, D.; Vera-Otarola, J.; Kuciak, M.; Corbin, A.; Nigumann, P.; Cristofari, G. Activation of individual L1 retrotransposon instances is restricted to cell-type dependent permissive loci. eLife 2016, 5, e13926. [Google Scholar] [CrossRef] [PubMed]
  12. Pitkänen, E.; Cajuso, T.; Katainen, R.; Kaasinen, E.; Välimäki, N.; Palin, K.; Taipale, J.; Aaltonen, L.A.; Kilpivaara, O. Frequent L1 retrotranspositions originating fromTTC28in colorectal cancer. Oncotarget 2014, 5, 853–859. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Rodriguez-Martin, B.; Alvarez, E.G.; Baez-Ortega, A.; Zamora, J.; Supek, F.; Demeulemeester, J.; Santamarina, M.; Ju, Y.S.; Temes, J.; Garcia-Souto, D.; et al. Pan-cancer analysis of whole genomes identifies driver rearrangements promoted by LINE-1 retrotransposition. Nat. Genet. 2020, 52, 306–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Taylor, M.S.; Lacava, J.; Mita, P.; Molloy, K.R.; Huang, C.R.L.; Li, N.; Adney, E.M.; Jiang, H.; Burns, K.H.; Chait, B.T.; et al. Affinity proteomics reveals human host factors implicated in discrete stages of LINE-1 retrotransposition. Cell 2013, 155, 1034–1048. [Google Scholar] [CrossRef] [Green Version]
  15. Deininger, P.; Morales, M.E.; White, T.B.; Baddoo, M.; Hedges, D.J.; Servant, G.; Srivastav, S.; Smither, M.E.; Concha, M.; Deharo, D.L.; et al. A comprehensive approach to expression of L1 loci. Nucleic Acids Res. 2017, 45, e31. [Google Scholar] [CrossRef]
  16. Ekreimer, U.; Schulz, W.A.; Ekoch, A.; Niegisch, G.; Egoering, W. HERV-K and LINE-1 DNA Methylation and Reexpression in Urothelial Carcinoma. Front. Oncol. 2013, 3. [Google Scholar] [CrossRef] [Green Version]
  17. Nüsgen, N.; Goering, W.; Dauksa, A.; Biswas, A.; Jamil, M.A.; Dimitriou, I.; Sharma, A.; Singer, H.; Fimmers, R.; Fröhlich, H.; et al. Inter-locus as well as intra-locus heterogeneity in LINE-1 promoter methylation in common human cancers suggests selective demethylation pressure at specific CpGs. Clin. Epigenetics 2015, 7, 17. [Google Scholar] [CrossRef] [Green Version]
  18. Vasudevan, A.A.J.; Kreimer, U.; Schulz, W.A.; Krikoni, A.; Schumann, G.G.; Häussinger, D.; Münk, C.; Egoering, W. APOBEC3B Activity Is Prevalent in Urothelial Carcinoma Cells and Only Slightly Affected by LINE-1 Expression. Front. Microbiol. 2018, 9. [Google Scholar] [CrossRef] [Green Version]
  19. Whongsiri, P.; Pimratana, C.; Wijitsettakul, U.; Sanpavat, A.; Jindatip, D.; Hoffmann, M.J.; Goering, W.; Schulz, W.A.; Boonla, C. Oxidative stress and LINE-1 reactivation in bladder cancer are epigenetically linked through active chromatin formation. Free. Radic. Biol. Med. 2019, 134, 419–428. [Google Scholar] [CrossRef]
  20. Rodić, N.; Sharma, R.; Sharma, R.; Zampella, J.; Dai, L.; Taylor, M.S.; Hruban, R.H.; Iacobuzio-Donahue, C.A.; Maitra, A.; Torbenson, M.S.; et al. Long Interspersed Element-1 Protein Expression Is a Hallmark of Many Human Cancers. Am. J. Pathol. 2014, 184, 1280–1286. [Google Scholar] [CrossRef] [Green Version]
  21. Eden, E.; Navon, R.; Steinfeld, I.; Lipson, D.; Yakhini, Z. GOrilla: A tool for discovery and visualization of enriched GO terms in ranked gene lists. BMC Bioinform. 2009, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Goodier, J.L.; Zhang, L.; Vetter, M.R.; Kazazian, H.H. LINE-1 ORF1 Protein Localizes in Stress Granules with Other RNA-Binding Proteins, Including Components of RNA Interference RNA-Induced Silencing Complex. Mol. Cell. Biol. 2007, 27, 6469–6483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Goodier, J.L.; Cheung, L.E.; Kazazian, J.H.H. Mapping the LINE1 ORF1 protein interactome reveals associated inhibitors of human retrotransposition. Nucleic Acids Res. 2013, 41, 7401–7419. [Google Scholar] [CrossRef] [PubMed]
  24. Hu, S.; Li, J.; Xu, F.; Mei, S.; Le Duff, Y.; Yin, L.; Pang, X.; Cen, S.; Jin, Q.; Liang, C.; et al. SAMHD1 Inhibits LINE-1 Retrotransposition by Promoting Stress Granule Formation. PLoS Genet. 2015, 11, e1005367. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Hubstenberger, A.; Courel, M.; Bénard, M.; Souquere, S.; Ernoult-Lange, M.; Chouaib, R.; Yi, Z.; Morlot, J.-B.; Munier, A.; Fradet, M.; et al. P-Body Purification Reveals the Condensation of Repressed mRNA Regulons. Mol. Cell 2017, 68, 144–157.e5. [Google Scholar] [CrossRef] [Green Version]
  26. Khong, A.; Matheny, T.; Jain, S.; Mitchell, S.F.; Wheeler, J.R.; Parker, R. The Stress Granule Transcriptome Reveals Principles of mRNA Accumulation in Stress Granules. Mol Cell 2017, 68, 808–820.e5. [Google Scholar] [CrossRef]
  27. Goering, W.; Ribarska, T.; Schulz, W.A. Selective changes of retroelement expression in human prostate cancer. Carcinogenesis 2011, 32, 1484–1492. [Google Scholar] [CrossRef]
  28. Hoffmann, M.J.; Müller, M.; Engers, R.; Schulz, W.A. Epigenetic control of CTCFL/BORIS and OCT4 expression in urogenital malignancies. Biochem. Pharmacol. 2006, 72, 1577–1588. [Google Scholar] [CrossRef]
  29. Knowles, M.A.; Hurst, C.D. Molecular biology of bladder cancer: New insights into pathogenesis and clinical diversity. Nat. Rev. Cancer 2015, 15, 25–41. [Google Scholar] [CrossRef]
  30. Ardeljan, D.; Steranka, J.P.; Liu, C.; Li, Z.; Taylor, M.S.; Payer, L.M.; Gorbounov, M.; Sarnecki, J.S.; Deshpande, V.; Hruban, R.H.; et al. Cell fitness screens reveal a conflict between LINE-1 retrotransposition and DNA replication. Nat. Struct. Mol. Biol. 2020, 27, 168–178. [Google Scholar] [CrossRef]
  31. Mita, P.; Sun, X.; Fenyö, D.; Kahler, D.J.; Li, D.; Agmon, N.; Wudzinska, A.; Keegan, S.; Bader, J.S.; Yun, C.; et al. BRCA1 and S phase DNA repair pathways restrict LINE-1 retrotransposition in human cells. Nat. Struct. Mol. Biol. 2020, 27, 179–191. [Google Scholar] [CrossRef] [PubMed]
  32. Woodcock, D.; Williamson, M.; Doherty, J. A Sensitive RNase Protection Assay to Detect Transcripts from Potentially Functional Human Endogenous L1 Retrotransposons. Biochem. Biophys. Res. Commun. 1996, 222, 460–465. [Google Scholar] [CrossRef] [PubMed]
  33. Belancio, V.P.; Roy-Engel, A.M.; Pochampally, R.R.; Deininger, P. Somatic expression of LINE-1 elements in human tissues. Nucleic Acids Res. 2010, 38, 3909–3922. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Briggs, E.M.; Ha, S.; Mita, P.; Brittingham, G.; Sciamanna, I.; Spadafora, C.; Logan, S.K. Long interspersed nuclear element-1 expression and retrotransposition in prostate cancer cells. Mob. DNA 2018, 9, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Erichsen, L.; Seifert, H.-H.; Schulz, W.A.; Hoffmann, M.J.; Niegisch, G.; Arauzo-Bravo, M.J.; Bendhack, M.L.; Poyet, C.; Hermanns, T.; Beermann, A.; et al. Basic Hallmarks of Urothelial Cancer Unleashed in Primary Uroepithelium by Interference with the Epigenetic Master Regulator ODC1. Sci. Rep. 2020, 10, 3808–3810. [Google Scholar] [CrossRef]
  36. Su, S.-F.; Abreu, A.L.D.C.; Chihara, Y.; Tsai, Y.; Andreu-Vieyra, C.; Daneshmand, S.; Skinner, E.C.; Jones, P.A.; Siegmund, K.; Liang, G. A Panel of Three Markers Hyper- and Hypomethylated in Urine Sediments Accurately Predicts Bladder Cancer Recurrence. Clin. Cancer Res. 2014, 20, 1978–1989. [Google Scholar] [CrossRef] [Green Version]
  37. Vafadar-Isfahani, N.; Parr, C.; McMillan, L.E.; Sanner, J.; Yeo, Z.; Saddington, S.; Peacock, O.; Cruickshanks, H.A.; Meehan, R.R.; Lund, J.N.; et al. Decoupling of DNA methylation and activity of intergenic LINE-1 promoters in colorectal cancer. Epigenetics 2017, 12, 465–475. [Google Scholar] [CrossRef]
  38. Seifert, H.-H.; Meyer, A.; Cronauer, M.V.; Hatina, J.; Müller, M.; Rieder, H.; Hoffmann, M.J.; Ackermann, R.; Schulz, W.A. A new and reliable culture system for superficial low-grade urothelial carcinoma of the bladder. World J. Urol. 2007, 25, 297–302. [Google Scholar] [CrossRef] [Green Version]
  39. Hoffmann, M.; Koutsogiannouli, E.; Skowron, M.A.; Pinkerneil, M.; Niegisch, G.; Brandt, A.; Stepanow, S.; Rieder, H.; Schulz, W.A. The New Immortalized Uroepithelial Cell Line HBLAK Contains Defined Genetic Aberrations Typical of Early Stage Urothelial Tumors. Bl. Cancer 2016, 2, 449–463. [Google Scholar] [CrossRef] [Green Version]
  40. Swiatkowski, S. Activities of MAP-Kinase Pathways in Normal Uroepithelial Cells and Urothelial Carcinoma Cell Lines. Exp. Cell Res. 2003, 282, 48–57. [Google Scholar] [CrossRef]
  41. Moran, V.A.; Niland, C.N.; Khalil, A.M. Co-Immunoprecipitation of Long Noncoding RNAs. Adv. Struct. Safe. Stud. 2012, 925, 219–228. [Google Scholar] [CrossRef]
  42. Hoffmann, M.J.; Dehn, J.; Droop, J.; Niegisch, G.; Niedworok, C.; Szarvas, T.; Schulz, W.A. Truncated Isoforms of lncRNA ANRIL Are Overexpressed in Bladder Cancer, But Do Not Contribute to Repression of INK4 Tumor Suppressors. Non-Coding RNA 2015, 1, 266–284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Penzkofer, T.; Dandekar, T.; Zemojtel, T. L1Base: From functional annotation to prediction of active LINE-1 elements. Nucleic Acids Res. 2004, 33, D498–D500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Rhead, B.L.; Karolchik, D.; Kuhn, R.M.; Hinrichs, A.S.; Zweig, A.S.; Fujita, P.A.; Diekhans, M.; Smith, K.E.; Rosenbloom, K.R.; Raney, B.J.; et al. The UCSC Genome Browser database: Update 2010. Nucleic Acids Res. 2010, 38, D613–D619. [Google Scholar] [CrossRef]
  45. Fleige, S.; Walf, V.; Huch, S.; Prgomet, C.; Sehm, J.; Pfaffl, M.W. Comparison of relative mRNA quantification models and the impact of RNA integrity in quantitative real-time RT-PCR. Biotechnol. Lett. 2006, 28, 1601–1613. [Google Scholar] [CrossRef]
Ijms 21 09433 g001 550
Figure 1. Comprehensive analysis of L1 expression in VM-Cub-1 UC cells. (A) Expression of ORF1p in representative UC cell lines analyzed by Western blotting. The non-transformed urothelial cell line HBLAK was used as a negative control and NCCIT embryonal carcinoma cells were used as a positive control. α-Tubulin was used as a loading control. (B) Experimental strategy to identify expression of individual L1s in bladder cancer cells. (C) Examples of mapping results from nanopore sequencing following RIP of VM-Cub-1 cells using the highly specific antibody against ORF1p. Read alignments are shown for the elements listed in Table 1. (D) Relative expression of the 90 L1 elements detected by RIP/nanopore sequencing in the two independent experiments (blue and orange) by rank.
Figure 1. Comprehensive analysis of L1 expression in VM-Cub-1 UC cells. (A) Expression of ORF1p in representative UC cell lines analyzed by Western blotting. The non-transformed urothelial cell line HBLAK was used as a negative control and NCCIT embryonal carcinoma cells were used as a positive control. α-Tubulin was used as a loading control. (B) Experimental strategy to identify expression of individual L1s in bladder cancer cells. (C) Examples of mapping results from nanopore sequencing following RIP of VM-Cub-1 cells using the highly specific antibody against ORF1p. Read alignments are shown for the elements listed in Table 1. (D) Relative expression of the 90 L1 elements detected by RIP/nanopore sequencing in the two independent experiments (blue and orange) by rank.
Ijms 21 09433 g001
Ijms 21 09433 g002 550
Figure 2. L1 expression across bladder cancer cells and normal urothelial cells. Measurements of RNA expression of overall LINE-1 (flL1) using the 5′-LINE1 assay (a) and the individual elements UID-33 (b), UID-59 (c), UID-60 (d), UID-66 (e) and UID-108 (f) were performed by quantitative real-time PCR. TBP was used as a reference gene.
Figure 2. L1 expression across bladder cancer cells and normal urothelial cells. Measurements of RNA expression of overall LINE-1 (flL1) using the 5′-LINE1 assay (a) and the individual elements UID-33 (b), UID-59 (c), UID-60 (d), UID-66 (e) and UID-108 (f) were performed by quantitative real-time PCR. TBP was used as a reference gene.
Ijms 21 09433 g002
Ijms 21 09433 g003 550
Figure 3. L1 expression in normal and cancerous bladder tissues. The expressions of L1s in bladder cancer (n = 51) and normal tissues (n = 36) were measured by qRT-PCR for overall long interspersed element 1 (LINE-1, L1) (a), and the individual elements UID-33 (b), UID-59 (c), UID-60 (d), UID-66 (e) and UID-108 (f). TBP was used as a reference gene. Statistical comparisons were performed by the Mann-Whitney U-test.
Figure 3. L1 expression in normal and cancerous bladder tissues. The expressions of L1s in bladder cancer (n = 51) and normal tissues (n = 36) were measured by qRT-PCR for overall long interspersed element 1 (LINE-1, L1) (a), and the individual elements UID-33 (b), UID-59 (c), UID-60 (d), UID-66 (e) and UID-108 (f). TBP was used as a reference gene. Statistical comparisons were performed by the Mann-Whitney U-test.
Ijms 21 09433 g003
Ijms 21 09433 g004 550
Figure 4. Expression of L1 elements in various cancer cell lines. The expression level of UID-33 (a), UID-59 (b), UID-60 (c), UID-66 (d) and UID-108 (e) were measured by qRT-PCR in different cancer cells including bladder cancer cell lines (BFTC-905, RT-112 and 5637), germ cell tumor cell lines (TERA1 and NCCIT) and prostate cancer cell lines (LNCaP and DU-145). TBP was used as a reference gene.
Figure 4. Expression of L1 elements in various cancer cell lines. The expression level of UID-33 (a), UID-59 (b), UID-60 (c), UID-66 (d) and UID-108 (e) were measured by qRT-PCR in different cancer cells including bladder cancer cell lines (BFTC-905, RT-112 and 5637), germ cell tumor cell lines (TERA1 and NCCIT) and prostate cancer cell lines (LNCaP and DU-145). TBP was used as a reference gene.
Ijms 21 09433 g004
Ijms 21 09433 g005 550
Figure 5. Histone modifications at intact L1s. ChIP-qPCR was performed in the cell lines VM-Cub-1 and 5637. (a) Expression of overall LINE-1, UID-33, UID-66 and UID-108 in VM-Cub-1 and 5637 cells as determined by qRT-PCR. ChIP was validated by qPCR for GAPDH (b), CTCFL (c) and GRM6 (d). (e)–(h) ChIP results for promoters of overall LINE-1 (e), UID-33 (f), UID-66 (g) and UID-108 (h). All ChIP results were calculated as percentages of input and the value of the respective normal IgG control (mouse or rabbit) was subtracted, which may yield negative values in cases of low enrichment. Mean values of one of two independent experiments are shown; each PCR was performed in triplicate.
Figure 5. Histone modifications at intact L1s. ChIP-qPCR was performed in the cell lines VM-Cub-1 and 5637. (a) Expression of overall LINE-1, UID-33, UID-66 and UID-108 in VM-Cub-1 and 5637 cells as determined by qRT-PCR. ChIP was validated by qPCR for GAPDH (b), CTCFL (c) and GRM6 (d). (e)–(h) ChIP results for promoters of overall LINE-1 (e), UID-33 (f), UID-66 (g) and UID-108 (h). All ChIP results were calculated as percentages of input and the value of the respective normal IgG control (mouse or rabbit) was subtracted, which may yield negative values in cases of low enrichment. Mean values of one of two independent experiments are shown; each PCR was performed in triplicate.
Ijms 21 09433 g005
Table
Table 1. Top L1 elements identified by RIP.
Table 1. Top L1 elements identified by RIP.
L1 IDRank Exp.1Expression Exp.1Rank Exp.2Expression Exp2.Localization ChromosomeLocalization BasesStrandClosest Gene
UID-108169.06423.90chr1p31.171,887,202–71,895,252+NEGR1
UID-66243.33130.76chr12q14.264,194,585–64,202,633+C12ORF66
UID-60339.94328.24chrXp21.136,464,177–36,472,219-Desert
UID-33436.08522.58ch17p13.19,614,984–9,623,031+CFAP52
UID-103524.15621.15chr1q25.1174,376,771–174,384,818-RABGAP1L
UID-140620.16812.48chr14q1226,628,254–26,636,301-NOVA1-AS1
UID-93719.25228.40chr3q23141,756,134–141,764,155-RNF7/GRK7
UID-59816.38718.01chrXp11.347,782,657–47,790,701-WASF4P
UID-135913.70115.94chr22q12.128,662,282–28,670,329+TTC28
UID-921010.74107.00chr3q25.1159,102,396–159,094,350-IQCJ-SCHIF1
UID-19119.74911.18chr16p12.321,049,706–21,041,661-DNAH3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Whongsiri, P.; Goering, W.; Lautwein, T.; Hader, C.; Niegisch, G.; Köhrer, K.; Hoffmann, M.J.; Schulz, W.A. Many Different LINE-1 Retroelements Are Activated in Bladder Cancer. Int. J. Mol. Sci. 2020, 21, 9433. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21249433

AMA Style

Whongsiri P, Goering W, Lautwein T, Hader C, Niegisch G, Köhrer K, Hoffmann MJ, Schulz WA. Many Different LINE-1 Retroelements Are Activated in Bladder Cancer. International Journal of Molecular Sciences. 2020; 21(24):9433. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21249433

Chicago/Turabian Style

Whongsiri, Patcharawalai, Wolfgang Goering, Tobias Lautwein, Christiane Hader, Günter Niegisch, Karl Köhrer, Michèle J. Hoffmann, and Wolfgang A. Schulz. 2020. "Many Different LINE-1 Retroelements Are Activated in Bladder Cancer" International Journal of Molecular Sciences 21, no. 24: 9433. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms21249433

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop