Skip to main content

Advertisement

A map of mobile DNA insertions in the NCI-60 human cancer cell panel

Mobile DNAvolume 7, Article number: 20 (2016) | Download Citation

Abstract

Background

The National Cancer Institute-60 (NCI-60) cell lines are among the most widely used models of human cancer. They provide a platform to integrate DNA sequence information, epigenetic data, RNA and protein expression, and pharmacologic susceptibilities in studies of cancer cell biology. Genome-wide studies of the complete panel have included exome sequencing, karyotyping, and copy number analyses but have not targeted repetitive sequences. Interspersed repeats derived from mobile DNAs are a significant source of heritable genetic variation, and insertions of active elements can occur somatically in malignancy.

Method

We used Transposon Insertion Profiling by microarray (TIP-chip) to map Long INterspersed Element-1 (LINE-1, L1) and Alu Short INterspersed Element (SINE) insertions in cancer genes in NCI-60 cells. We focused this discovery effort on annotated Cancer Gene Index loci.

Results

We catalogued a total of 749 and 2,100 loci corresponding to candidate LINE-1 and Alu insertion sites, respectively. As expected, these numbers encompass previously known insertions, polymorphisms shared in unrelated tumor cell lines, as well as unique, potentially tumor-specific insertions. We also conducted association analyses relating individual insertions to a variety of cellular phenotypes.

Conclusions

These data provide a resource for investigators with interests in specific cancer gene loci or mobile element insertion effects more broadly. Our data underscore that significant genetic variation in cancer genomes is owed to LINE-1 and Alu retrotransposons. Our findings also indicate that as large numbers of cancer genomes become available, it will be possible to associate individual transposable element insertion variants with molecular and phenotypic features of these malignancies.

Significance statement

Transposable elements are repetitive sequences that comprise much of our DNA. They create both inherited and somatically acquired structural variants. Here, we describe a first generation map of LINE-1 and Alu insertions in NCI-60 cancer cell lines. This provides a resource for discovering and testing functional consequences of these sequences.

Background

The National Cancer Institute-60 (NCI-60) cell panel was developed in the 1980s as a tool for pharmacologic screens and has become the most extensively studied collection of human cancers [1]. The panel comprises 59 cell lines encompassing nine tissue origins of malignancy, including blood, breast, colon, central nervous system, kidney, lung, ovary, prostate, and skin [2]. They have become a resource for high throughput characterizations and systems biology based approaches to cancer.

NCI-60 cell genomes have been described by targeted [3] and whole exome sequencing [4], karyotyping [5], and assays to detect copy number alteration [6], loss of heterozygosity [7], and DNA methylation [8]. Large scale mRNA [9] and microRNA [10] expression, protein abundance [11] and phosphorylation [12], and metabolomic [13] studies have also been conducted. Because assays are applied across the panel of cell lines in each case, datasets from orthogonal studies can be related to one another. For example, gene expression patterns have been found to be predictive of chemotherapeutic sensitivities [9].

Interspersed repeats have not been incorporated in these or many other genome-wide surveys. These repetitive sequences are dynamic constituents of human genomes and important sources of structural variation [14–20]. RNA transcribed from active elements can be reverse transcribed and integrated into the genome at new sites by proteins encoded by LINE-1 (Long INterspersed Element)-1 [21–23]. The result is that relatively recent insertions of LINE-1 (L1Hs) and Alu SINEs (AluYa5, AluYa8, AluYb8, AluYb9) are sources of genetic polymorphisms where both the pre-insertion allele and the insertion allele coexist in human populations. Moreover, LINE-1 sequences are hypomethylated [24–28] and express protein in a wide variety of human cancers [29], and somatic LINE-1 integrations have been reported in tumor genomes [15, 30–36].

It is well established that inherited and acquired mobile DNA insertions can affect gene expression; there is inherent potential for insertions to have effects on tumor biology. However, the large majority occur in intronic or intergenic regions. Strong biases in the distribution of insertion sites or recurrent ‘hotspots’ for insertions arising during tumor development are frequently not obvious, leading to the presumption that most are non-functional ‘passenger mutations’ [34, 36].

This is not such a tumor-normal comparison study, but rather, one aimed to identify potential functions of mobile DNAs in human cancer cells. Towards this end, we mapped LINE-1 and Alu insertions in the NCI-60 tumor cell panel. We used a method for interspersed repeat mapping, Transposon Insertion Profiling by microarray (TIP-chip), to identify insertion sites. We also use previous characterizations of the cell panel to associate specific insertions with cellular phenotypes.

Results

Transposon insertion profiling by microarray

To map mobile DNA insertions, we used a method we have termed transposon insertion profiling by microarray (TIP-chip), which uses vectorette PCR to amplify unknown sequence adjacent to a known primer-binding site (Fig. 1a). We surveyed three major currently active mobile DNAs in humans (L1Hs, AluYa5/8; and AluYb8/9) as previously described [14]. To focus on the potential functional impact of these sequences on cancer cell phenotypes, PCR amplicons were labeled and analyzed using a genomic tiling microarray designed to encompass 6,484 known Cancer Gene Index loci (+/- 10 kb) (Biomax™ Informatics), about 17 % of the genome. Peaks of signal intensity correspond to TE insertions (Fig. 1a, b); known LINE-1 and Alu elements incorporated in the reference genome assembly (hereafter, ‘reference insertions’) were used as a quality control metric and to set cut-offs for recognized peaks (Fig. 1c).

Fig. 1
figure1

Mapping transposable element (TE) insertion sites. a. A schematic illustrating the sequential steps of Transposon Insertion Profiling by microarray (TIP-chip). (1) An interval of double stranded genomic DNA with two TE insertions (boxes) oriented on opposing strands is shown; (2) the DNA is digested in parallel restriction enzyme reactions and ligated to vectorette oligonucleotides; (3) oligonucleotides complementary to the TE insertions prime first strand synthesis; (4) the elongating strands form reverse complements of the vectorette sequence; (5) there is exponential amplification of insertion site fragments; (6) these amplicons are labeled and hybridized to genomic tiling microarrays; and (7) ‘peaks’ of fluorescence intensity across several probes corresponding to contiguous genomic positions indicate a TE insertion. b. An example of a polymorphic Alu peak in two leukemia cell lines (SR and MOLT-4) in the third intron of the TCOF1 (Treacher Collins-Franceschetti syndrome 1) gene on chromosome 5. The upper panels show TIP-chip data for the insertion, which is present in the SR line and not the MOLT-4 cells. The Alu insertion is a minus (-) strand insertion to the right of the probe with the greatest intensity; an arrow is drawn to indicate its position and orientation, but the arrow is not drawn to scale. Alu insertions approximate 300 bp, and the width of the peak in this case is 5 kb. c. Peaks were recognized using a sliding window algorithm which identified adjacent probes above a threshold fluorescence intensity value. The threshold value was progressively lowered to identify peaks in a rank order. The graphs show the number of reference insertions identified verses peak rank for a representative LINE-1 and Alu TIP-chip. The cut-off for defining a candidate insertion was established using the inflection points (red arrows) of these plots

Full size image

A total of 749 and 2,100 peaks corresponding to candidate LINE-1 and Alu insertion sites respectively were recognized across the NCI-60 cell panel. These locations were cross-referenced to previously described insertions to define three categories: (i.) reference insertions, which include invariant insertions and insertion polymorphisms incorporated in the reference genome assembly; (ii.) inherited variants either previously described (known polymorphic) or newly discovered, but occurring in multiple, unrelated cell lines (novel polymorphic); and (iii.) novel, ‘singleton’ insertions seen uniquely in one cell line (Fig. 2a, b). The last category includes both insertions that were constitutive (germline) in the patient from whom the cell line was derived as well as somatic insertions acquired during tumor development or the propagation of these cell lines. A greater proportion of LINE-1 insertions were singletons (68 %) compared with Alu insertions (21 %). Density plots for both LINE-1 and Alu show most peaks fall into this last category, particularly for L1Hs, although a biphasic distribution was seen (Fig. 3a, b).

Fig. 2
figure2

Total TE insertions. The stacked bar plots show the relative numbers of novel variants, known variants, and reference insertions per cell line for LINE-1 (green, upper panel) and Alu (red, lower panel). The total number of insertions detected per cell line is similar across the tumor panel

Full size image
Fig. 3
figure3

Distribution of TE insertions across the NCI-60 panel. a. Individual insertions are arrayed in order of frequency horizontally, and cell lines are arrayed vertically. Yellow denotes presence of insertion; blue denotes absence. LINE-1 are on the upper plot, and Alu are on the lower. Cell types are listed for the lower panel, and the ordering is the same in the upper panel. b. The density plot shows proportions of insertions against the numbers of cell lines containing an insertion. For both Alu (red) and LINE-1 (green), there is a bimodal distribution. The leftmost density reflects a large number of polymorphic insertions with low allele frequencies and (for LINE-1 singletons) somatically acquired insertions. The rightmost increase in density shows common variants or fixed insertions present in most or all cell lines

Full size image

Our array encompassed 130 known reference LINE-1 and 1278 Alu insertions. A total of 112 LINE-1 and 1,160 Alu insertions detected were present in the reference genome assembly. A total of 697 LINE-1 and 1,147 Alu insertions were singleton or polymorphic (known and novel) segregating in human populations (Fig. 2a, b). Insertions incorporated in the reference genome that are known to be polymorphic are counted in both groups. A summary of insertion positions by tumor type and cell line can be found in Additional file 1: Table S1, Additional file 2: Table S2.

We found that each cell line had a unique transposable element (TE) insertion profile (Fig. 3a). After correcting for batch effects, a principal component analyses (PCA) did not show clustering by tumor type. As expected, however, pairs of cell lines derived from the same individual grouped together, and these pairs showed a high concordance of top-ranking peaks as compared to unrelated cell lines. We compared TE insertion profiles to described cytogenetic abnormalities. In some instances, insertions were informative of deletions; for example, a reference LINE-1 in the retinoblastoma 1 (RB1) locus was only absent in the MB468 breast cancer cell line, consistent with the homozygous deletion of RB1 reported for this cell line [37].

Insertions in genes involved in oncogenesis

In TIP-chip, probe spacing does not resolve insertions to the precise base, and insertion strandedness was not predicted for all peak intervals in this study. Despite these limitations, we identified peak intervals that partially or entirely overlapped exon intervals for further inspection. Partial overlaps were almost entirely attributable to insertions near an exon. We identified 9 insertions within exons, and all were located within gene 3’ untranslated regions (3’ UTRs); none affected protein open reading frames.

To begin to approach potential functional consequences of intronic insertions, we analyzed insertion sites in sets of genes with described roles in cancer. We considered collections of genes with TE insertions while grouping together malignant cell lines by tissue of origin. Interestingly, in breast cancer cell lines, we observed a significant enrichment of singleton and polymorphic LINE-1 and Alu insertions in “STOP genes”, defined in shRNA screens as suppressors of human mammary epithelial cell proliferation [38] (p = 1.23x10-9) (Fig. 4a). This result persisted when LINE-1 and Alu insertions were analyzed independently; LINE-1 singleton insertions but not Alu singleton insertions were also enriched in this gene set (Fig. 4b). Analysis of expression of these “STOP” genes shows that a preponderance of these genes are down-regulated; this result persists in those genes containing a TE insertion. The findings suggest that collectively, insertions may act to compromise expression of these genes.

Fig. 4
figure4

TE enrichment analyses. a. STOP and GO genes have been implicated in breast cancer as genes that appear to inhibit and promote tumor development, respectively. Using a hypergeometric distribution to assess enrichment, we found that TE insertions are enriched in STOP genes on the array (p = 1.23x10-9) but not in GO genes (p = 0.33). b. The bar graph shows enrichment by type of TE plotted as the negative log of the p-value. No GO gene enrichment is seen. STOP gene enrichment is seen considering all LINE-1 (p = 3.11x10-3); all Alu (p = 2.27x10-10); as well as LINE-1 singletons (p = 4.16x10-5). c. Insertions were also enriched in common insertion sites (CIS) (p = 1.46x10-4); COSMIC commonly mutated cancer genes (p = 7.74x10-10); and genes reported to acquire somatic LINE-1 insertions in cancer by Lee et al. (p = 5.34x10-14). d. Genes with TE insertions in ovarian cancer cell lines are more likely than other genes to be downregulated in ovarian cancer samples as compared to normal tissue controls. Randomly selected genes are shown for comparison (bottom panel)

Full size image

Consistent with this model, ovarian cancer cell lines showed a preponderance of insertions in genes that are down regulated in ovarian cancers as compared to normal tissue. A random set of genes from the array is shown as a histogram for comparison (Fig. 4d). This pattern was absent in other tumor types.

We saw an enrichment of singleton and polymorphic TEs in genes recurrently mutated in experimental cancer models and in human tumors. For the former, we considered common insertion sites (CIS) defined as gene loci recurrently interrupted by insertional mutagens in forward cancer gene screens in mice [39, 40] (p = 1.46x10-4). The latter was assessed using genes frequently mutated in human cancers taken from the Catalogue Of Somatic Mutations In Cancer (COSMIC) database [41] (p = 7.74x10-10) (Fig. 4c). We also compared our insertion profiles to sites of reported somatic TE insertions in human cancers. We analyzed novel (singleton and polymorphic) insertions and discovered that we had overlaps in 22 of the 64 genes noted by Lee et al. [32] and 23 of 76 from Solimini et al. [38](Fig. 4c). We anticipate the possibility that common insertion site loci will be identified as more insertion site mapping studies are conducted in human tumors.

Functional associations of individual insertions

An advantage of working with the NCI-60 cell lines is that these are well studied. To integrate our insertion site maps with other findings in these cells, we performed COMPARE analyses [42]. COMPARE is a pattern matching method developed specifically for NCI-60 cell lines that provides a p-value for each association (S5–25). Direct, local roles for TEs (in cis) were not observed for the majority of correlations. However, COMPARE did reveal three insertions associated with DNA hypermethylation within 30 kb of the insertion site. For example, a polymorphic Alu insertion in the SS18L1 (Synovial sarcoma translocation gene on chromosome 18-like 1) gene locus oriented anti-sense to the transcription of the gene, is associated with increased methylation of nearby CpG sites at the same gene locus (p = 6.67x10-6) (Fig. 5a).

Fig. 5
figure5

TE insertions associated with cellular phenotypes. a. Associations with DNA methylation. (Upper panel) Diagram of the SS18L1 (Synovial sarcoma translocation gene on chromosome 18-like 1) gene locus, which contains an antisense Alu associated with increased CpG methylation at that gene (i.e., in cis, p = 3.67x10-6) (Middle panel) Manhattan plot showing TE positions on the x-axis and strengths of association with gene methylation on the y-axis (Bonferroni-corrected p-values). Singleton insertions were excluded from association analyses. Gene abbreviations are given for both the gene in which the insertion is found (red) and the associated methylation site (black) in examples. The TE insertion at the CSRP2 (cysteine and glycine-rich protein 2) gene locus was associated with methylation at 22 distinct loci (Lower panel). The associated methylation sites are distributed throughout the genome. b. Manhattan plot showing associations with steady state mRNA levels. Gene abbreviations are given for both the gene in which the insertion is found (red) and the associated transcript level (black); in these two examples, the TE is associated with upregulation of the mRNA. c. Manhattan plot showing associations with drug sensitivity as measured by total cellular growth inhibition. The gene in which the insertion is found is given (red), as well as the associated pharmacologic agent (black)

Full size image

Manhattan plots illustrate highly significant correlations found in trans (Fig. 5a–c). A subset of insertions had multiple associations (vertical series of dots corresponding to one TE location), suggesting the possibility of pleomorphic effects of an insertion haplotype.

In addition, we encountered examples of single ‘driver’ mutations and cellular phenotypes that could be associated with multiple TE insertions. Five insertions correlated with a mutation in the ERBB2 gene (v-erb-b2 erythroblastic leukemia viral oncogene homolog 2, the HER2/neu locus), and more than 10 insertions were associated with thymidylate synthase activity (p values < 10-20). To probe relationships between multiple trans associated factors related to a single TE insertion, we performed pathway analyses on sets of genes, each encompassing the TE insertion locus and all RNAs and proteins with associated expression patterns. This yielded more than 250 curated pathways with enrichment p-values less than 10-4, supporting the concept that these are biologically relevant as opposed to spurious associations. All COMPARE results are provided in the (Additional file 3: Table S3).

Discussion

Our genomes are filled with highly repetitive DNA sequences derived from TEs. Tailored methods for their detection, including TIP-chip [14], targeted insertion site sequencing [15, 17, 18, 31, 36, 43], and algorithms for finding variants in whole genome sequencing [20, 34, 44] are revealing this previously masked dimension of genomic data. Collectively, these studies confirm that TEs are rich sources of genetic diversity in human populations, and provide evidence that they are somatically unstable in a variety of tumor types. Of the two most active germline elements, LINE-1 and Alu, (which is mobilized in trans by LINE-1-encoded proteins), LINE-1 has been more well documented to be active in cancer. Alu insertions account for more inherited polymorphisms. For both types of TEs, the vast majority of catalogued insertions are intronic and intergenic without clear function.

To begin a systematic survey for functionally consequential LINE-1 and Alu integrations in human neoplasias, we mapped these variants in the NCI-60 cell panel. NCI-60 is a unique resource for this, encompassing a variety of cancer cell lines that have the advantages of being well studied and readily available. We mapped LINE-1 and Alu insertion positions using a microarray-based approach over a large census of cancer genes. Even as TIP-chip is replaced by sequencing, we expect these data will provide a useful reference.

TIP-chip across the NCI-60 panel revealed numerous novel candidate TEs, totaling about 500 L1Hs and 1000 AluYa/Yb insertions distributed across the 60 cell lines. These include insertions that are unique to a cell line (‘singleton’) and novel polymorphic insertions (found in unrelated cell lines). Although ‘singletons’ may be enriched for tumor-specific, somatic insertion events, matched non-neoplastic cells for the corresponding patient cases are not available, and therefore we cannot definitively differentiate somatic from inherited variants. Similarly, these cell lines have undergone numerous passages since their creation, and somatic insertion events occurring in culture cannot be clearly recognized. We note a greater proportion of LINE-1 singletons (68 % of LINE-1 loci) than Alu singletons (21 % of Alu loci), consistent with ongoing LINE-1 retrotransposition in vivo or in vitro.

We approached the question of TE function by two avenues. We first tested for biases in the distribution of insertions with respect to known gene sets. We found a preferential accumulation of TE insertions in retained copies of ‘STOP genes’ in breast cancer cell lines; these gene loci function as inhibitors of mammary epithelial cell proliferation. Experimental models suggest that it is advantageous for tumor growth to compromise the function of these genes [38], and we speculate that TE insertions are enriched at these loci because they have a role in this process. These ‘STOP genes’ are downregulated in the breast cancer cell lines, as is the subset of ‘STOP genes’ containing TE insertions. We also found preferential TE accumulation in genes downregulated in ovarian cancers compared with normal ovarian tissue, which would be consistent with this model. Finally, genes with functional roles in cancer were also more commonly seen as insertion sites than expected. These included genes ‘hit’ recurrently by insertional mutagenesis in forward genetics screens in mice, the so-called common insertion sites (CIS), and in genes commonly mutated in human cancers (COSMIC catalog) [41].

We note that the exonizations of intronic LINE-1 [45] and Alu sequences [46] are being increasingly recognized using RNA-seq, and that many of the resulting transcripts have an altered protein coding capacity. It may be possible to identify aberrant mRNA species corresponding to these insertion loci and thus invoke a molecular mechanism to underlie this type of functional effect.

Our second approach relied on association studies. We used existing data in COMPARE analyses to test for relationships between TE insertion alleles and cellular phenotypes. In the case of DNA methylation only, cis effects could be seen relating individual TEs with local DNA hypermethylation. We identified three Alu integrations associated with DNA hypermethylation at the insertion site (+/- 30 kb). The most notable is a polymorphic Alu insertion in the first intron of the SS18L1 (synovial sarcoma translocation gene on chromosome 18-like 1) gene locus associated with CpG hypermethylation at the same locus (p = 3.67x10-6). SS18 and SS18L1 encode transcriptional regulators and are breakpoints in chromosomal translocations in synovial sarcoma [47]. These translocations are not seen in the NCI-60 panel tumors, and whether the epigenetic signature associated with the Alu insertion impacts expression of this gene is unknown. So, while it is not clear at this point that SS18L1 methylation is germane to the development of these malignancies, our ability to relate genotype and epigenetics at these sites demonstrate the value of this approach.

The large majority of statistically significant associations between insertions and cellular phenotypes appeared to involve indirect or trans effects that are difficult to test further. Pathway analyses suggest that many are not random, but reflect recognized, related gene sets. It may be that the indirect effects can be dissected for some insertion alleles; particularly promising may be those at loci of transcriptional regulators with definable target genes [29].

Conclusions

In summary, we profiled LINE-1 and Alu insertion sites in a panel of widely used cancer cell lines, the NCI-60. We expect maps such as these will be a useful resource for experimentalists with interests in how transposable element insertions interact with genes. Our analyses show that insertion sites can be integrated with other data to develop testable hypotheses about the function of mobile DNAs in cancer.

Methods

NCI-60 cell lines

The National Cancer Institute-60 (NCI-60) human cancer cells are a group of 60 cell lines representing nine different types of neoplasias(breast cancer, colon cancer, CNS tumor, leukemia, lung cancer, melanoma, ovarian cancer, prostate cancer, and renal cell carcinoma) composed of 54 individual cancer cases and three pairs of cell lines (ADR and OVCAR-8; MB-435 and M14; and SNB19 and U251) with each pair originating from the same patient [48, 49]. The NCI-60 panel has been extensively characterized in a breadth of molecular and pharmacologic assay [50]. Genomic DNA was obtained directly from the NCI.

Microarray design

A genomic tiling micorarrray was designed to cover the NCI Cancer Gene Index (disease list). A total of 6,484 RefSeq gene identifiers were extracted from the. XML file and converted to genomic coordinates corresponding to each transcript unit +/- 10 kb hg19 reference genome assembly (February 2009, GRCh37). UCSC Table Browser intervals were merged using GALAXY [51], and probes were chosen for the NimbleGen HD (2.1 M feature) array platform by the manufacturer (Roche NimbleGen, Madison, WI).

Transposon insertion profiling by microarray (TIP-chip)

Five micrograms of genomic DNA of each cell line was digested overnight in parallel reactions using four restriction enzymes (AseI, BspHI, HindIII, and Xbal). Sticky ends were ligated to annealed, partially complementary vectorette oligonucleotide adapters. Each template was aliquoted into 3 separate vectorette PCR reactions for L1Hs, AluYa5/8, and AluYb8/9 mobile DNA families. These were then labeled with Cy3-dUTP for LINE-1 and Cy5-dUTP for Alu and hybridized to Nimblegen genome tiling arrays according to the manufacturer’s instructions. Reference insertions are those incorporated in the Feb. 2009 assembly of the human genome (hg19, GRCh37 Genome Reference Consortium Human Reference 37, GCA_000001405.1).

Peak recognition

Each scanned array yielded a raw .tff file, which was processed using Nimblescan v2.5 (Roche Nimblegen, Madison, WI) to give genomic coordinates and probe intensities (.gff files). A PERL script removed probes overlapping repeats to reduce noise (RepeatMasking). Nimblescan called peaks using a sliding window threshold. Peaks were ranked by the threshold of the log2 transformed ratio of red (Alu) and green (L1) channels or the reciprocal (settings: percent (p) start = 90, p step = 1, #steps = 76, width of sliding window = 1500 bp, min probes > 4, all probes > 2). The top 5,000 L1 and Alu peaks were kept for evaluation.

Peak cut-off

Among these peaks, recovery of those corresponding to mobile DNA insertions in hg19 (reference insertions) was used as a proxy of assay performance. Reference insertion count was plotted against peaks recognized (Fig. 1c). A cut-off was imposed on the peak threshold value (p >70 for L1 and p > 60 for Alu) to include peaks up to the approximate inflection point of this curve in subsequent analyses. These threshold values were altered for outlier cell lines to reflect the curve inflection point. MYSQL was used to annotate peaks with respect to genes and known mobile DNA insertions (L1Hs, AluY, AluYa5, AluYa8, AluYb8, and AluYb9 using 1−2 kb margins). Lists of known insertions were obtained from previously published databases [14, 19, 52, 53].

Clustering and insertion profiles

Principle component analysis (PCA) (R-package) was used to remove batch effect. All insertions were sorted by density across the cell lines and plotted as a matrix. Cell lines lacking high-frequency insertions were assessed for karyotype abnormalities manually.

COMPARE analysis

Reference and non-reference insertions were analyzed using a COMPARE analysis [42] associating each with the CellMiner database of NCI-60 cell profiling studies. These have included DNA mutations and methylation; RNA and miRNA expression; protein expression, enzymatic activity; and drug inhibition studies. Associations for those insertions found in one cell line (singleton) were considered only for cis effects and were discarded from other associations due to their high false-positive rates. P-values for other insertions were corrected using Bonferroni multiple test correction and plotted using the start position of peak intervals to generate Manhattan plots (adaptation of Genetics Analysis Package, R-package).

Pathway analysis

Gene loci containing candidate non-reference (polymorphic and singleton) LINE-1 and Alu insertions and associated gene names from RNA and protein COMPARE analysis were uploaded in batch to the MSigDb ‘Investigate Gene Sets’ from the Broad Institute Gene Set Enrichment Analysis web interface [54] (using the C2 curated gene sets). Pathways were selected if the insertion locus was part of the pathway and the p-value of the pathway was less than 10-4. Interactome plots were used to visualize relationships between genes in pathways using Search Tool for the Retrieval of Interacting Genes/Protein (STRING) 9.0 [55]. Plots were adapted to show the gene locus containing the insertion (yellow) and the direction of related correlations (red for positive correlations with the insertion; purple for negative correlations).

Preferential integration sites

To investigate preferential transposable element insertion in genes implicated in oncogenesis and mouse common insertion sites, we used a hypergeometric distribution test (pHypr R-package) which controlled for genes tiled on the array. Results were plotted using the –log(p-value).

Tumor-normal gene expression studies

Tumor vs normal gene expression for genes containing candidate non-reference TE insertions was assessed for each tumor type using large tumor/normal gene expression databases. Tumor gene to normal gene expression ratios were obtained using NCBI GEO2R [56]. GEO2R was used to log2 transform expression data if datasets were not in log2 formats. Value distribution of all databases was assessed for median-centering prior to evaluation. Expression values for all insertion-containing genes was plotted as a horizontal bar plot. A random sample of 1000 genes from the array were evaluated in the same manner to serve as a control set. A histogram of random gene expression values was plotted. Databases (Breast = GSE5764, Ovarian = GSE26712, omitted samples with “no evidence of disease”, Colon = GSE6988, omitted non-primary tumors, Melanoma = GSE7553, CNS = GSE4290, non-tumor used as “normal” and non-glioblastomas omitted, Prostate = GSE3325, Renal = GSE11151, non-conventional tumors omitted, NSCL = GSE19188).

STOP gene expression in breast cancer cell lines

Expression of STOP genes containing candidate non-reference TE insertions was assessed using log2 transformed Agilent mRNA expression data [57] obtained from the CellMiner for the Breast cancer cell lines. The expression was averaged across all cell lines, sorted, and plotted as a horizontal bar plot. STOP genes tiled on the array, but without a TE insertion was plotted as well. Tumor-Normal expression for STOP genes was performed according to the methods used above in Tumor-Normal gene expression studies.

Abbreviations

LINE-1:

Long INterspersed Element-1

NCI:

National Cancer Institute

SINE:

Short INterspersed Element

TIP-chip:

Transposon insertion profiling by microarray

References

  1. 1.

    Shoemaker RH. The NCI60 human tumour cell line anticancer drug screen. Nat Rev Cancer. 2006;6:813–23.

  2. 2.

    Stinson SF, Alley MC, Kopp WC, Fiebig HH, Mullendore LA, Pittman AF, Kenney S, Keller J, Boyd MR. Morphological and immunocytochemical characteristics of human tumor cell lines for use in a disease-oriented anticancer drug screen. Anticancer Res. 1992;12:1035–53.

  3. 3.

    Ikediobi ON, Davies H, Bignell G, Edkins S, Stevens C, O'Meara S, Santarius T, Avis T, Barthorpe S, Brackenbury L, et al. Mutation analysis of 24 known cancer genes in the NCI-60 cell line set. Mol Cancer Ther. 2006;5:2606–12.

  4. 4.

    Abaan OD, Polley EC, Davis SR, Zhu YJ, Bilke S, Walker RL, Pineda M, Gindin Y, Jiang Y, Reinhold WC, et al. The exomes of the NCI-60 panel: a genomic resource for cancer biology and systems pharmacology. Cancer Res. 2013;73:4372–82.

  5. 5.

    Roschke AV, Tonon G, Gehlhaus KS, McTyre N, Bussey KJ, Lababidi S, Scudiero DA, Weinstein JN, Kirsch IR. Karyotypic complexity of the NCI-60 drug-screening panel. Cancer Res. 2003;63:8634–47.

  6. 6.

    Varma S, Pommier Y, Sunshine M, Weinstein JN, Reinhold WC. High resolution copy number variation data in the NCI-60 cancer cell lines from whole genome microarrays accessible through Cell Miner. PLoS ONE. 2014;9:e92047.

  7. 7.

    Ruan X, Kocher JP, Pommier Y, Liu H, Reinhold WC. Mass homozygotes accumulation in the NCI-60 cancer cell lines as compared to HapMap Trios, and relation to fragile site location. PLoS ONE. 2012;7:e31628.

  8. 8.

    Shen L, Kondo Y, Ahmed S, Boumber Y, Konishi K, Guo Y, Chen X, Vilaythong JN, Issa JP. Drug sensitivity prediction by CpG island methylation profile in the NCI-60 cancer cell line panel. Cancer Res. 2007;67:11335–43.

  9. 9.

    Staunton JE, Slonim DK, Coller HA, Tamayo P, Angelo MJ, Park J, Scherf U, Lee JK, Reinhold WO, Weinstein JN, et al. Chemosensitivity prediction by transcriptional profiling. Proc Natl Acad Sci U S A. 2001;98:10787–92.

  10. 10.

    Patnaik SK, Dahlgaard J, Mazin W, Kannisto E, Jensen T, Knudsen S, Yendamuri S. Expression of microRNAs in the NCI-60 cancer cell-lines. PLoS ONE. 2012;7:e49918.

  11. 11.

    Nishizuka S, Charboneau L, Young L, Major S, Reinhold WC, Waltham M, Kouros-Mehr H, Bussey KJ, Lee JK, Espina V, et al. Proteomic profiling of the NCI-60 cancer cell lines using new high-density reverse-phase lysate microarrays. Proc Natl Acad Sci U S A. 2003;100:14229–34.

  12. 12.

    Federici G, Gao X, Slawek J, Arodz T, Shitaye A, Wulfkuhle JD, De Maria R, Liotta LA, Petricoin 3rd EF. Systems analysis of the NCI-60 cancer cell lines by alignment of protein pathway activation modules with "-OMIC" data fields and therapeutic response signatures. Mol Cancer Res. 2013;11:676–85.

  13. 13.

    Jain M, Nilsson R, Sharma S, Madhusudhan N, Kitami T, Souza AL, Kafri R, Kirschner MW, Clish CB, Mootha VK. Metabolite profiling identifies a key role for glycine in rapid cancer cell proliferation. Science. 2012;336:1040–4.

  14. 14.

    Huang CR, Schneider AM, Lu Y, Niranjan T, Shen P, Robinson MA, Steranka JP, Valle D, Civin CI, Wang T, et al. Mobile interspersed repeats are major structural variants in the human genome. Cell. 2010;141:1171–82.

  15. 15.

    Iskow RC, McCabe MT, Mills RE, Torene S, Pittard WS, Neuwald AF, Van Meir EG, Vertino PM, Devine SE. Natural mutagenesis of human genomes by endogenous retrotransposons. Cell. 2010;141:1253–61.

  16. 16.

    Beck CR, Collier P, Macfarlane C, Malig M, Kidd JM, Eichler EE, Badge RM, Moran JV. LINE-1 retrotransposition activity in human genomes. Cell. 2010;141:1159–70.

  17. 17.

    Ewing AD, Kazazian Jr HH. High-throughput sequencing reveals extensive variation in human-specific L1 content in individual human genomes. Genome Res. 2010;20:1262–70.

  18. 18.

    Witherspoon DJ, Xing J, Zhang Y, Watkins WS, Batzer MA, Jorde LB. Mobile element scanning (ME-Scan) by targeted high-throughput sequencing. BMC Genomics. 2010;11:410.

  19. 19.

    Hormozdiari F, Alkan C, Ventura M, Hajirasouliha I, Malig M, Hach F, Yorukoglu D, Dao P, Bakhshi M, Sahinalp SC, Eichler EE. Alu repeat discovery and characterization within human genomes. Genome Res. 2011;21:840–9.

  20. 20.

    Stewart C, Kural D, Stromberg MP, Walker JA, Konkel MK, Stutz AM, Urban AE, Grubert F, Lam HY, Lee WP, et al. A comprehensive map of mobile element insertion polymorphisms in humans. PLoS Genet. 2011;7:e1002236.

  21. 21.

    Mathias SL, Scott AF, Kazazian Jr HH, Boeke JD, Gabriel A. Reverse transcriptase encoded by a human transposable element. Science. 1991;254:1808–10.

  22. 22.

    Feng Q, Moran JV, Kazazian Jr HH, Boeke JD. Human L1 retrotransposon encodes a conserved endonuclease required for retrotransposition. Cell. 1996;87:905–16.

  23. 23.

    Dewannieux M, Esnault C, Heidmann T. LINE-mediated retrotransposition of marked Alu sequences. Nat Genet. 2003;35:41–8.

  24. 24.

    Alves G, Tatro A, Fanning T. Differential methylation of human LINE-1 retrotransposons in malignant cells. Gene. 1996;176:39–44.

  25. 25.

    Jurgens B, Schmitz-Drager BJ, Schulz WA. Hypomethylation of L1 LINE sequences prevailing in human urothelial carcinoma. Cancer Res. 1996;56:5698–703.

  26. 26.

    Lin CH, Hsieh SY, Sheen IS, Lee WC, Chen TC, Shyu WC, Liaw YF. Genome-wide hypomethylation in hepatocellular carcinogenesis. Cancer Res. 2001;61:4238–43.

  27. 27.

    Chalitchagorn K, Shuangshoti S, Hourpai N, Kongruttanachok N, Tangkijvanich P, Thong-ngam D, Voravud N, Sriuranpong V, Mutirangura A. Distinctive pattern of LINE-1 methylation level in normal tissues and the association with carcinogenesis. Oncogene. 2004;23:8841–6.

  28. 28.

    Estecio MR, Gharibyan V, Shen L, Ibrahim AE, Doshi K, He R, Jelinek J, Yang AS, Yan PS, Huang TH, et al. LINE-1 hypomethylation in cancer is highly variable and inversely correlated with microsatellite instability. PLoS ONE. 2007;2:e399.

  29. 29.

    Rodic N, Sharma R, Sharma R, Zampella J, Dai L, Taylor MS, Hruban RH, Iacobuzio-Donahue CA, Maitra A, Torbenson MS, et al. Long interspersed element-1 protein expression is a hallmark of many human cancers. Am J Pathol. 2014;184:1280–6.

  30. 30.

    Miki Y, Nishisho I, Horii A, Miyoshi Y, Utsunomiya J, Kinzler KW, Vogelstein B, Nakamura Y. Disruption of the APC gene by a retrotransposal insertion of L1 sequence in a colon cancer. Cancer Res. 1992;52:643–5.

  31. 31.

    Solyom S, Ewing AD, Rahrmann EP, Doucet T, Nelson HH, Burns MB, Harris RS, Sigmon DF, Casella A, Erlanger B, et al. Extensive somatic L1 retrotransposition in colorectal tumors. Genome Res. 2012;22:2328–38.

  32. 32.

    Lee E, Iskow R, Yang L, Gokcumen O, Haseley P, Luquette 3rd LJ, Lohr JG, Harris CC, Ding L, Wilson RK, et al. Landscape of somatic retrotransposition in human cancers. Science. 2012;337:967–71.

  33. 33.

    Shukla R, Upton KR, Munoz-Lopez M, Gerhardt DJ, Fisher ME, Nguyen T, Brennan PM, Baillie JK, Collino A, Ghisletti S, et al. Endogenous retrotransposition activates oncogenic pathways in hepatocellular carcinoma. Cell. 2013;153:101–11.

  34. 34.

    Tubio JM, Li Y, Ju YS, Martincorena I, Cooke SL, Tojo M, Gundem G, Pipinikas CP, 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:1251343.

  35. 35.

    Doucet-O'Hare TT, Rodic N, Sharma R, Darbari I, Abril G, Choi JA, Young Ahn J, Cheng Y, Anders RA, Burns KH, et al. LINE-1 expression and retrotransposition in Barrett's esophagus and esophageal carcinoma. Proc Natl Acad Sci U S A. 2015;112:E4894–900.

  36. 36.

    Rodic N, Steranka JP, Makohon-Moore A, Moyer A, Shen P, Sharma R, Kohutek ZA, Huang CR, Ahn D, Mita P, et al. Retrotransposon insertions in the clonal evolution of pancreatic ductal adenocarcinoma. Nat Med. 2015;21:1060–4.

  37. 37.

    T'Ang A, Varley JM, Chakraborty S, Murphree AL, Fung YK. Structural rearrangement of the retinoblastoma gene in human breast carcinoma. Science. 1988;242:263–6.

  38. 38.

    Solimini NL, Xu Q, Mermel CH, Liang AC, Schlabach MR, Luo J, Burrows AE, Anselmo AN, Bredemeyer AL, Li MZ, et al. Recurrent hemizygous deletions in cancers may optimize proliferative potential. Science. 2012;337:104–9.

  39. 39.

    Akagi K, Suzuki T, Stephens RM, Jenkins NA, Copeland NG. RTCGD: retroviral tagged cancer gene database. Nucleic Acids Res. 2004;32:D523–7.

  40. 40.

    de Ridder J, Uren A, Kool J, Reinders M, Wessels L. Detecting statistically significant common insertion sites in retroviral insertional mutagenesis screens. PLoS Comput Biol. 2006;2:e166.

  41. 41.

    Bamford S, Dawson E, Forbes S, Clements J, Pettett R, Dogan A, Flanagan A, Teague J, Futreal PA, Stratton MR, Wooster R. The COSMIC (Catalogue of Somatic Mutations in Cancer) database and website. Br J Cancer. 2004;91:355–8.

  42. 42.

    Paull KD, Shoemaker RH, Hodes L, Monks A, Scudiero DA, Rubinstein L, Plowman J, Boyd MR. Display and analysis of patterns of differential activity of drugs against human tumor cell lines: development of mean graph and COMPARE algorithm. J Natl Cancer Inst. 1989;81:1088–92.

  43. 43.

    Baillie JK, Barnett MW, Upton KR, Gerhardt DJ, Richmond TA, De Sapio F, Brennan PM, Rizzu P, Smith S, Fell M, et al. Somatic retrotransposition alters the genetic landscape of the human brain. Nature. 2011;479:534–7.

  44. 44.

    Burns KH, Boeke JD. Human transposon tectonics. Cell. 2012;149:740–52.

  45. 45.

    Denli AM, Narvaiza I, Kerman BE, Pena M, Benner C, Marchetto MC, Diedrich JK, Aslanian A, Ma J, Moresco JJ, et al. Primate-specific ORF0 contributes to retrotransposon-mediated diversity. Cell. 2015;163:583–93.

  46. 46.

    Schwartz S, Gal-Mark N, Kfir N, Oren R, Kim E, Ast G. Alu exonization events reveal features required for precise recognition of exons by the splicing machinery. PLoS Comput Biol. 2009;5:e1000300.

  47. 47.

    Storlazzi CT, Mertens F, Mandahl N, Gisselsson D, Isaksson M, Gustafson P, Domanski HA, Panagopoulos I. A novel fusion gene, SS18L1/SSX1, in synovial sarcoma. Genes Chromosomes Cancer. 2003;37:195–200.

  48. 48.

    Ellison G, Klinowska T, Westwood RF, Docter E, French T, Fox JC. Further evidence to support the melanocytic origin of MDA-MB-435. Mol Pathol. 2002;55:294–9.

  49. 49.

    Garraway LA, Widlund HR, Rubin MA, Getz G, Berger AJ, Ramaswamy S, Beroukhim R, Milner DA, Granter SR, Du J, et al. Integrative genomic analyses identify MITF as a lineage survival oncogene amplified in malignant melanoma. Nature. 2005;436:117–22.

  50. 50.

    Ross DT, Scherf U, Eisen MB, Perou CM, Rees C, Spellman P, Iyer V, Jeffrey SS, Van de Rijn M, Waltham M, et al. Systematic variation in gene expression patterns in human cancer cell lines. Nat Genet. 2000;24:227–35.

  51. 51.

    Afgan E, Baker D, van den Beek M, Blankenberg D, Bouvier D, Cech M, Chilton J, Clements D, Coraor N, Eberhard C, et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2016 update. Nucleic Acids Res. 2016;44(W1):W3–W10.

  52. 52.

    Mir AA, Philippe C, Cristofari G. euL1db: the European database of L1HS retrotransposon insertions in humans. Nucleic Acids Res. 2015;43(Database issue):D43–7.

  53. 53.

    Wang J, Song L, Grover D, Azrak S, Batzer MA, Liang P. dbRIP: a highly integrated database of retrotransposon insertion polymorphisms in humans. Hum Mutat. 2006;27:323–9.

  54. 54.

    Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50.

  55. 55.

    Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, Doerks T, Stark M, Muller J, Bork P, et al. The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res. 2011;39:D561–8.

  56. 56.

    Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.

  57. 57.

    Liu H, D'Andrade P, Fulmer-Smentek S, Lorenzi P, Kohn KW, Weinstein JN, Pommier Y, Reinhold WC. mRNA and microRNA expression profiles of the NCI-60 integrated with drug activities. Mol Cancer Ther. 2010;9:1080–91.

Download references

Acknowledgements

We thank Peilin Shen, Jared Steranka, Youngran Park, Xuan Pham, and Ashley Castillo for technical assistance and thoughtful discussion of the project; Eitan Halper-Stromberg for computational assistance; and Beatriz Villaba-Martín for Fig. 1a. Susan L. Holbeck at the National Cancer Institute performed the COMPARE analysis, and Joel Bader provided approaches for pathway analyses.

Funding

This work was supported by R01CA163705, R01GM103999, and a Career Award for Medical Scientists from the Burroughs Welcome Foundation (to KHB); and P50GM107632 (to JDB). JGZ was supported by the Howard Hughes Medical Institute (HHMI) Medical Research Fellows Program.

Availability of data and material

Data are available as GSE83756.

Authors’ contributions

JGZ carried out the TIP-chip experiments; NR, WRY, CRLH, JW, VPG, and TCC optimized the protocols and/or conducted data analysis; JGZ, JDB, and KHB designed the study and secured funding for the project; JGZ and KHB drafted the manuscript. All authors contributed to the manuscript review.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Affiliations

  1. Department of Dermatology, Johns Hopkins University School of Medicine, 733 North Broadway, Miller Research Building Room 469, Baltimore, MD, 21205, USA

    • John G. Zampella

  2. Department of Pathology, Johns Hopkins University School of Medicine, 733 North Broadway, Miller Research Building Room 469, Baltimore, MD, 21205, USA

    • Nemanja Rodić
    • , Wan Rou Yang
    • , Toby C. Cornish
    •  & Kathleen H. Burns

  3. McKusick-Nathans Institute of Genetic Medicine, 733 North Broadway, Miller Research Building Room 469, Baltimore, MD, 21205, USA

    • Cheng Ran Lisa Huang
    • , Jane Welch
    • , Veena P. Gnanakkan
    • , Jef D. Boeke
    •  & Kathleen H. Burns

  4. High Throughput (HiT) Biology Center, 733 North Broadway, Miller Research Building Room 469, Baltimore, MD, 21205, USA

    • Jef D. Boeke
    •  & Kathleen H. Burns

  5. The Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, 733 North Broadway, Miller Research Building Room 469, Baltimore, MD, 21205, USA

    • Kathleen H. Burns

  6. Present address: Institute for Systems Genetics, NYU Langone University School of Medicine, New York, NY, 10016, USA

    • Jef D. Boeke

Authors

  1. Search for John G. Zampella in:

  2. Search for Nemanja Rodić in:

  3. Search for Wan Rou Yang in:

  4. Search for Cheng Ran Lisa Huang in:

  5. Search for Jane Welch in:

  6. Search for Veena P. Gnanakkan in:

  7. Search for Toby C. Cornish in:

  8. Search for Jef D. Boeke in:

  9. Search for Kathleen H. Burns in:

Corresponding author

Correspondence to Kathleen H. Burns.

Additional files

Additional file 1

A map of LINE-1 (L1) insertion site positions in the NCI-60 cell panel. Genomic coordinates of TIP-chip peaks are provided. Reference insertions are indicated in column D (hg19), and known polymorphic variants are indicated by a ‘Y’ in column E (Y/N, yes/no). For each cell line in columns G-BN, a ‘1’ indicates that the insertion is present, while ‘0’ indicates that the insertion is not found. (XLSX 85 kb)

Additional file 2

A map of Alu insertion site positions in the NCI-60 cell panel. Genomic coordinates of TIP-chip peaks are provided. Reference insertions are indicated in column D (hg19), and known polymorphic variants are indicated by a ‘Y’ in column E (Y/N, yes/no). For each cell line in columns G-BN, a ‘1’ indicates that the insertion is present, while ‘0’ indicates that the insertion is not found. (XLSX 380 kb)

Additional file 3

COMPARE analysis associating insertions with other cell characteristics. Different tabs are used for different datasets. Activity, enzyme activity measures; Decreased / Increased methylation, DNA methylation measures; Metabolome, metabolic intermediates; Drug Effect GI50, concentration for 50 % growth inhibition; Drug Effect TGI, concentration for total growth inhibition; miRNA, microRNA expression levels; RNA, mRNA expression levels; Mutations, somatically-acquired DNA mutations; Protein, protein expression. (XLSX 4049 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Breast Cancer Cell Line
  • Transposable Element Insertion
  • Reference Genome Assembly
  • Stop Gene
  • Mammary Epithelial Cell Proliferation

Mobile DNA

ISSN: 1759-8753