- Research article
- Open Access
Targeted nanopore sequencing for the identification of ABCB1 promoter translocations in cancer
BMC Cancer volume 20, Article number: 1075 (2020)
Resistance to chemotherapy is the most common cause of treatment failure in acute myeloid leukemia (AML) and the drug efflux pump ABCB1 is a critical mediator. Recent studies have identified promoter translocations as common drivers of high ABCB1 expression in recurrent, chemotherapy-treated high-grade serous ovarian cancer (HGSC) and breast cancer. These fusions place ABCB1 under the control of a strong promoter while leaving its open reading frame intact. The mechanisms controlling high ABCB1 expression in AML are largely unknown. We therefore established an experimental system and analysis pipeline to determine whether promoter translocations account for high ABCB1 expression in cases of relapsed human AML.
The human AML cell line THP-1 was used to create a model of chemotherapy resistance in which ABCB1 expression was driven by a promoter fusion. The THP-1 model was used to establish a targeted nanopore long-read sequencing approach that was then applied to cases of ABCB1high HGSC and AML. H3K27Ac ChIP sequencing was used to assess the activity of native promoters in cases of ABCB1high AML.
Prolonged in vitro daunorubicin exposure induced activating ABCB1 promoter translocations in human THP-1 AML cells, similar to those recently described in recurrent high-grade serous ovarian and breast cancers. Targeted nanopore sequencing proved an efficient method for identifying ABCB1 structural variants in THP-1 AML cells and HGSC; the promoter translocations identified in HGSC were both previously described and novel. In contrast, activating ABCB1 promoter translocations were not identified in ABCB1high AML; instead H3K27Ac ChIP sequencing demonstrated active native promoters in all cases studied.
Despite frequent high level expression of ABCB1 in relapsed primary AML we found no evidence of ABCB1 translocations and instead confirmed high-level activity of native ABCB1 promoters, consistent with endogenous regulation.
High level expression of the drug efflux pump ABCB1 leading to active efflux from blast cells of standard-of-care chemotherapy drugs such as daunorubicin and cytarabine is one of the strongest predictors for induction failure in acute myeloid leukemia (AML) . Interestingly, recent studies have shown promoter translocations to be common drivers of ABCB1 expression in recurrent, chemotherapy-treated high-grade serous ovarian (HGSC) and breast cancer [2, 3]. As a consequence ABCB1 is placed under the control of a strong promoter while leaving its open reading frame intact. The treatment-emergent selection of promoter-translocated ABCB1 variants provides compelling evidence for a central role of drug efflux in chemotherapy resistance, at least in cancers of ovary and breast. By contrast, the cellular mechanisms controlling ABCB1 expression in AML are largely unknown. Most publicly available genome sequencing data in AML were generated using whole-exome sequencing which does not detect structural variants (SVs) in non-coding sequences . Furthermore, few studies have included samples from relapsed patients previously treated with intensive chemotherapy, raising a question as to whether ABCB1 promoter SVs might be present in ABCB1high AML relapse cases. We therefore established an experimental system and analysis pipeline to determine whether promoter translocations account for high ABCB1 expression in cases of relapsed human AML.
THP-1 cells were from DSMZ (Braunschweig, Germany) and cultured in RPMI 1640 medium (Sigma Aldrich) supplemented with 2 mM L-Glutamine (Life Technologies, Carlsbad, CA) and 10% fetal bovine serum (Sigma Aldrich). Whilst under drug selection cells were counted and replated every third day. Cell lines were confirmed mycoplasma-free and authenticated by short tandem repeat DNA profiling.
Daunorubicin and verapamil were from Sigma Aldrich (St. Louis, MO); tariquidar was from Generon (Slough, UK). Compounds were resuspended in DMSO (tariquidar) or ddH20 (verapamil and daunorubicin), aliquoted and stored at − 20 °C. Final DMSO concentration was < 0.5% in all experiments.
Primary human samples
Primary human AML and HGSC samples were from the Manchester Cancer Research Centre Tissue Biobank (approved by the South Manchester Research Ethics Committee). Their use was authorized by the Tissue Biobank’s scientific sub-committee, with the informed consent of donors. For relapsed AML, the Biobank archive was searched for all samples. Where there were too few cryopreserved cells available for potential downstream analyses, or where the percentage of blasts in the cryopreserved sample was less than 80%, cases were excluded. This left 21 separate samples for analysis (Table S1). Three presentation cases, selected at random, were also included for comparison. For ChIP sequencing, selected samples were thawed and immediately crosslinked.
Processing of primary HGSC samples
Samples for ABCB1 quantitative PCR and nanopore sequencing were taken from ex vivo models established using primary ascitic fluid samples, as previously described . Briefly, ascites was centrifuged (500×g for 10 min at 4 °C) and cell pellets pooled in HBSS (Life Technologies). Red blood cells were removed using a red blood cell lysis buffer (Miltenyi Biotec) as per the manufacturer’s instructions. Tumour cells were plated into Primaria flasks containing OCMI. All cultures were incubated for 2–4 days at 37 °C in a humidified 5% CO2 and 5% O2 atmosphere. Media was replaced every 3–4 days. Upon cell attachment, stromal cells were separated from the mixed sample using 0.05% trypsin-EDTA. Once tumour cells reached 95% confluency, cells were passaged using 0.25% Trypsin-EDTA, centrifuged in DMEM containing 20% FBS and re-plated at a 1:2 ratio. For long-term storage, cells were frozen in Bambanker (Wako pure chemical).
Cell viability assays
5 × 103 cells were plated in each well of a 96-well plate with media containing a serial dilution of daunorubicin. Plates were incubated for 72 h at 37 °C. 20 μl of 140 μg/mL resazurin (Sigma Aldrich) was added to each well. Plates were then incubated for a further 4 h and read using a POLARstar Omega plate reader (BMG Labtech, Aylesbury, UK).
RNA sequencing and data analysis
Total RNA was extracted from 5 × 105 cells using QIAshredder spin columns and an RNeasy® Plus Micro kit (Qiagen, Manchester, UK). Prior to sequencing RNA integrity was checked using a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). PolyA libraries were prepared using a SureSelect Strand Specific RNA Library Prep kit (Agilent Technologies) and samples were then barcoded and pooled. Sequencing was performed using a NextSeq system (Illumina, San Diego, CA). Two replicates were sequenced for each cell line (THP-1_S and THP-1_R). A single run of 150 bp paired-end sequencing produced a mean of 55.1 M reads per sample. Reads were aligned to the human genome (hg38) using STAR v2.4.2a . DEseq2 was used to calculate FPKM (fragments per kilobase of transcript per million mapped reads) values for each transcript .
cDNA was generated using a High Capacity Reverse Transcription kit (Applied Biosystems, Foster City, CA). qPCR reactions were performed in MicroAmp® optical 384-well reaction plates and analysed using a QuantStudio® 5 PCR system (Applied Biosystems). Reactions were performed in quadruplicate and included primers for β-Actin (ACTB) as a housekeeping gene. Primers were designed using the Universal Probe Library (UPL) Assay Design Center (Roche, Basel, Switzerland) and purchased from Integrated DNA Technologies (Coralville, IA). Raw fluorescence data was converted to CT values using the Thermo Fisher Cloud facility (Waltham, MA) and normalised to ACTB. Primers were (i) ACTB (F) ATTGGCAATGAGCGGTTC, (R) GGATGCCACAGGACTCCAT, UPL probe #11 and (ii) ABCB1 (F) GGAAATTTAGAAGATCTGATGTCAAAC, (R) ACTGTAATAATAGGCATACCTGGTCA, UPL probe #65.
Chromatin immunoprecipitation (ChIP) and next generation sequencing
ChIP was performed using anti-H3K27Ac (ab4729, Abcam, Cambridge, UK). 108 cells (THP-1) or 107 cells (primary AML samples) were used for each precipitation using the method described by Lee et al. . Briefly, cells were cross-linked with 1% formaldehyde for 10 min at room temperature before the reaction was quenched with 0.125 M glycine. Cell pellets were washed twice with PBS and nuclear lysates sonicated for 6 cycles using a Bioruptor® Pico (Diagenode, Liege, Belgium). 10 μg of antibody bound to 100 μl of magnetic beads (Dynabeads Protein G, Invitrogen, Carlsbad, CA) was added to each sample and immunoprecipitation performed overnight on a rotator at 4 °C and 20 rpm. After five washes with RIPA buffer (50 mM HEPES pH 7.6, 1 mM EDTA, 0.7% Na deoxycholate, 1% NP-40, 0.5 M LiCl), chromatin IP-bound fractions were extracted by incubating for 15 min at 65 °C with elution buffer (50 mM TrisHCl pH 8, 10 mM EDTA, 1% SDS). Crosslinking was then reversed by incubation at 65 °C for 6 h. RNaseA (1 mg/ml) and proteinase K (20 mg/ml) were added to eliminate RNA and protein from the samples. DNA was extracted using phenol:chloroform:isoamyl alcohol and precipitated by adding 2 volumes of ice-cold 100% ethanol, glycogen (20 μg/μl), 200 mM NaCl and freezing at − 80 °C for at least 1 h. Pellets were washed with 70% ethanol and eluted in 50 μl 10 mM TrisHCl pH 8.0.
Libraries were prepared for sequencing using a Microplex Library Preparation Kit (Diagenode). 200-800 bp fragments were selected using AMPure XP beads (Beckman Coulter, Brea, CA). Sequencing was performed using a NextSeq desktop sequencing system (Illumina) with 75 bp, paired-end high output generating 65-80 M (THP-1_S and THP-1_R) and 34-45 M reads per sample (primary AML samples). Reads were aligned to the human genome (hg38) using BWA-MEM v0.7.15 . Read duplicates were removed using Picard v2.1.0. Reads were further filtered using Bedtools v2.25.0 to keep only paired reads that mapped to standard chromosomes and to remove reads with a mapping quality of less than 10. Reads mapped to blacklisted regions defined by ENCODE were then removed using Bedtools (http://mitra.stanford.edu/kundaje).
4C primer sequences and enzyme combinations were selected using the University of Chicago online tool (http://mnlab.uchicago.edu/4Cpd) with co-ordinates from the ABCB1 promoter active in THP-1_R cells (hg38, chr7:87,598,302-87,601,399). 4C sequencing was performed according to the protocol developed by Splinter et al. . Briefly, 107 cells were cross-linked with 2% formaldehyde for 10 min at room temperature before the reaction was quenched with 0.125 M glycine. Cells were lysed with buffer containing 50 mM Tris–HCl pH 7.5, 150 mM NaCl, 5 mM EDTA, 0.5% NP-40, 1% TX-100 and 1x complete protease inhibitors (Roche, #11245200). The cross-linked nuclear preparation was then incubated with DpnII. Digestion was confirmed by reversing crosslinking for an aliquot and running on a 0.6% agarose gel. Samples were then ligated overnight at 16 °C using T4 DNA ligase (Roche, #799009). Ligation efficiency was again confirmed with 0.6% agarose gel. Crosslinking was reversed and DNA extracted using phenol-chloroform. Samples were then subjected to a second digestion using Csp6I. Ligation was again performed overnight at 16 °C using T4 DNA ligase. DNA was then extracted using phenol-chloroform and purified with a QIAquick PCR purification kit (Qiagen, #28104). PCR primers were designed to incorporate 4C primers with a barcode and Illumina adapter sequences thus:
Reading primer: 5′ P5-Barcode-Primer 3′.
Non-reading primer: 5′ P7-Primer 3′.
Illumina adapter sequences:
P5: AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGA TCT.
Non-reading primer: CAAGCAGAAGACGGCATACGAGATAGGGTAGGTATTCCACTTTT.
Reading primer THP-1_S: AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTA CACGACGCTCTTCCGATCTGCCAATGAGATACCAGGTCTGATC (barcode GCCAAT).
Reading primer THP-1_R: AATGATACGGCGACCACCGAGATCTACACTCTTTCCCT.
ACACGACGCTCTTCCGATCTCTTGTAGAGATACCAGGTCTGATC (barcode CTTGTA).
PCR was performed with Expand Long Template Polymerase (Roche, #11759060001) using 3.2 μg of 4C template product and then purified using a High Pure PCR Product Purification Kit (Roche. #11732676001). Library quality was assessed using a 2100 Bioanalyzer (Agilent Technologies). Samples were sequenced with 10% phiX using a MiSeq desktop sequencing system (Illumina) with 75 bp single-end settings generating a mean of 1.4 M reads per sample. Sequencing data was deconvoluted using cutadapt v1.18. Reads were mapped and analysis performed using 4Cseqpipe .
Target enrichment and nanopore long-read sequencing
Genomic DNA was extracted from 5 × 106 cells using a QIAamp DNA mini kit (Qiagen, Manchester, UK). 1 μg of genomic DNA was sheared with a g-TUBE™ (Covaris, Inc., Woburn, MA) using 50 μl sample volume and a 30s spin at 11,000 rpm followed by inversion and a further 30s at 11,000 rpm. Barcoded libraries were prepared from 45 μl of fragmented DNA using a Ligation Sequencing Kit SQK-LSK108 / SQK-LSK109 (Oxford Nanopore Technologies, Oxford, UK) according to the 1D Sequence Capture protocol with the following modifications: the pre hybridisation PCR was replaced by a custom PCR with 6 cycles using Herculase II fusion DNA polymerase (Agilent Technologies). The subsequent AMPure XP purification was performed using a 0.6x bead to sample volume ratio. All purified product was used for the hybridization using an oligonucleotide designed using the online Agilent SureDesign service (https://earray.chem.agilent.com/suredesign/home.htm). The library targeted a 500 kb region that included ABCB1 and upstream neighbouring genes (RUNDC3B, SLC25A40, DBF4 and ADAM22; chr7:87,132,949-87,632,948 hg19). Following pull-down, the post-capture amplification was replaced by a custom PCR with 14 cycles as described above. Amplified libraries were quantified by Qubit (Invitrogen) and TapeStation 2200 (Agilent Technologies). An additional amplification was performed using 10 cycles of PCR and the product was quantified. Amplified libraries were again quantified by Qubit and 1 μg of each library used for the End-prep. AMPure XP purification was performed using a 0.47x bead to sample ratio. Four end-prepped libraries were pooled into one adapter ligation reaction using 125 ng of each, as measured by Tapestation. A 0.47x bead to sample ratio was used for the AMPure XP bead binding. Sequencing was performed on the MinION using the R9.4 / 9.4.1 Flow Cell FLO-MIN106 / FLO-MIN106D (Oxford Nanopore Technologies) according to the manufacturer’s instructions. Target enriched DNA from THP-1_R was also subjected to conventional sequencing using a MiSeq desktop sequencing system (Illumina) with 1% phiX and 251 bp, paired-end settings generating 1.92 M reads. Nanopore sequencing of THP-1_R generated 0.92 M reads with a mean read length of 1901 bp and a maximum read length of 224.8 kb. The coverage of the targeted region was 58.1x. Nanopore sequencing of AML/HGSC generated an average of 1.6 million reads per sample (range 0.8–3.0 million) with a mean read length of 1581 bp. FASTQ files were aligned against the GRCh38/hg38 reference human genome using LAST after training the aligner with a subsample of 10,000 reads . Structural variant calling was performed using NanoSV . Structural variants were confirmed through manual inspection and identification of chimeric reads using IGV (https://software.broadinstitute.org/software/igv/).
Flow cytometry & assessment of daunorubicin retention
Flow cytometry was performed using an LSR II flow cytometer (BD Biosciences, Franklin Lakes, NJ). FlowJo v10.1 (BD Biosciences) was used to analyze data. To assess daunorubicin retention 5 × 105 cells were resuspended in PBS containing 1 μM daunorubicin with 40 μM verapamil or 5 nM tariquidar or vehicle. Samples were incubated for 2 h at 37 °C then placed on ice to prevent further efflux. Daunorubicin accumulation was assessed by flow cytometry.
In the first instance we sought to establish whether activating promoter translocations could be experimentally induced in AML. ABCB1 promoter translocations have been described in several solid malignancy cancer cell lines but not in human myeloid leukemia cell lines . We exposed the human AML cell line THP-1 to escalating doses of daunorubicin over 142 days and generated a resistant line (THP-1_R) which exhibited a 28.3-fold increase in daunorubicin IC50 (Figs. S1A and S1B). RNA sequencing and confirmatory quantitative PCR indicated that ABCB1, which was not expressed in THP-1 cells, was one of the most highly expressed genes in THP-1_R cells (among the top 2.5% of genes by read count; Table S2, Fig. S1C). Parental drug-sensitive THP-1 (THP-1_S) cells did not efflux daunorubicin whereas THP-1_R cells exhibited robust drug efflux that was completely reversed by either verapamil (a non-specific ABC transporter substrate) or tariquidar (a highly specific inhibitor of ABCB1) , confirming that efflux was due to ABCB1 (Fig. S1D). Importantly, the resistance of THP-1_R cells to daunorubicin was a stable cellular phenotype and persisted over 3 months of continuous culture without daunorubicin selection, suggesting constitutive expression (Fig. S1E).
ABCB1 codes for four protein coding transcript variants and contains two promoters separated by a ~ 110 kb intron that is transcribed in variants 1–3. This large first intron, which overlaps part of the coding sequence for RUNDC3B (Fig. 1a), appears to be a common site for promoter translocations. In fact, all of the high-grade serous ovarian cancer samples with ABCB1 SVs characterised by Christie et al.  had at least one SV involving intron 1. We observed that THP-1_R cells expressed ABCB1 transcript variant 4 but inspection of RNA sequencing tracks revealed additional transcription of intron 1, starting ~ 48 kb upstream of the variant 4 promoter (Fig. 1a). ChIP sequencing for H3K27 acetylation, a histone modification that marks active promoters and enhancers, found no acetylation of either ABCB1 promoter confirming that neither was active despite high levels of transcription (Fig. 1a). To identify the co-opted promoter we performed 4C-sequencing of THP-1_S and THP-1_R cells with a viewpoint centered on the ABCB1 variant 4 promoter and observed in resistant cells both an abrupt loss of contacts within ABCB1 intron 1 and gain of a distant interaction with GTF2I (Fig. 1a and b). The GTF2I promoter is located ~ 12.94 million base pairs centromeric to the ABCB1 variant 4 promoter on chromosome 7q; it was strongly acetylated and GTF2I expression was among the top 25% of genes in both lines (Fig. 1a and Table S2). These data indicate that a daunorubicin-induced translocation fusing the GTF2I promoter with the coding sequence for ABCB1 was responsible for the acquisition of constitutive high level ABCB1 expression in THP-1_R cells.
Having demonstrated the presence of an experimentally initiated ABCB1-activating promoter translocation in THP-1_R cells, we developed a long-read sequencing technique to confirm this finding and to facilitate screening of primary patient AML samples for similar SVs. We performed targeted sequencing of a 500 kb region of chromosome 7 that contained both ABCB1 and SLC25A40, which is the most common translocation partner for ABCB1 in chemotherapy-treated breast and ovarian cancer patients . The region was enriched using a library of biotinylated oligonucleotides that provided 95.5% coverage of the targeted sequence (Fig. 1c). Uniquely mapping oligonucleotides cannot be produced for regions containing repetitive sequences. Given that breakpoints frequently contain repeats we first compared conventional sequencing to long-read nanopore sequencing, which generates reads of up to 15 kb that can bridge translocations involving repetitive sequences (Fig. 1c) . Both approaches identified the GTF2I promoter: ABCB1 translocation (Figs. 1a and 2a). However, despite lower coverage of the targeted region (58.1x versus 368.6x for MiSeq) nanopore long-read sequencing generated 52 chimeric reads that crossed the breakpoint and mapped unambiguously to GTF2I, compared with only six by conventional sequencing (Fig. 2b). These data demonstrate that long-read nanopore targeted sequencing efficiently identified the GTF2I promoter: ABCB1 translocation.
In chemotherapy-treated HGSC Christie et al.  identified ABCB1 fusions in 20/108 (18.5%) of cases, increasing to 31.3% if only ABCB1high cases were considered. To assess the sensitivity of our nanopore long-read sequencing approach we quantified ABCB1 expression in 33 samples of metastatic HGSC recovered from ascitic fluid (Fig. 2c), and sequenced the 15 cases with the highest expression (Table S3). In 5/15 cases we identified translocations involving splicing of non-coding 5′ regions of partner genes to exon 2 of ABCB1 as a result of an intron 1 breakpoint (Table S4 and Fig. 2c-d); in one sample three separate fusion genes were identified. In addition to expanding the repertoire of reported ABCB1 fusions in HGSC, our data also demonstrate that targeted nanopore sequencing is at least as sensitive as the multimodal approach taken by Christie et al. .
We next identified in our Biobank 21 cases of relapsed AML with sufficient cryopreserved cells for downstream analyses, and once more determined ABCB1 expression levels by quantitative PCR (Fig. 2e). All but two patients had previously been treated with anthracycline-containing chemotherapy and all had initially achieved a complete morphological remission. Patient characteristics and molecular genetic features are shown in Table S1. At original presentation 16/21 cases exhibited a normal karyotype, 2/21 an MLL gene rearrangement, 1/21 a t (8;21), 1/21 inv.  and 1/21 trisomy 11 (Table S1). Blood samples exhibited high blast cell percentages (median 93.5%, range 82–100%). Targeted nanopore sequencing of 12 cases with high level ABCB1 expression (Fig. 2e) achieved similar coverage to that achieved in THP-1_R (mean 45.3x, range 16.5–82.3), but failed to identify any functional translocations involving ABCB1 intron 1 (Figs. S2A and S2B). H3K27Ac ChIP sequencing of nine of 12 of the nanopore long read sequenced samples performed in parallel revealed acetylation of one or both ABCB1 promoters in every sample, demonstrating that high level ABCB1 expression in relapsed AML is driven by native promoters (Fig. 2f), at least in the samples evaluated. These results are consistent with our recent finding that ABCB1 expression in primary AML is highly dynamic and regulated by several active enhancers .
While we were able to experimentally induce in a human myeloid leukemia cell line the kind of treatment-emergent promoter-translocated ABCB1 variant found in chemotherapy-treated breast and ovarian cancers, we were unable to identify similar translocations in primary human relapsed AML samples from patients previously exposed to anthracycline-based chemotherapy. This is despite the strong association of ABCB1 with treatment failure and the high frequency of balanced translocations targeting a range of genes in AML [1, 18]. Instead, as we have recently demonstrated, ABCB1 expression appears to be regulated by its native promoters and a range of cooperating enhancer elements . One possible explanation for this difference is the relative genetic stability of AML which has one of the lowest mutation frequencies of any cancer . By contrast ovarian and breast cancers have much higher rates of mutation and more frequent microsatellite and chromosomal instability that predicts for poor outcome [20,21,22]. The mechanism of activation may also depend on the state of the ABCB1 promoter in different cancers and their cells of origin. Hematopoietic stem cells normally express ABCB1 and its regulatory mechanisms appear to be preserved in AML, allowing for dynamic expression and adaptation following treatment . In contrast, in tumours and tissues where the ABCB1 promoter is silenced, translocation may be the only mechanism by which expression can be induced, favouring the selection of structural variants that replace the silenced promoter. Despite the poor performance in clinical trials of inhibitors of ABCB1, recent findings continue to highlight the central clinical relevance of this drug efflux pump.
Our findings demonstrate the utility of long-read sequencing as an efficient technique for identifying SVs in human cancers and a promising approach for cost effective medical genetics . Whilst our approach expanded the repertoire of reported ABCB1 fusions in HGSC, we did not identify similar structural variants in ABCB1high AML, where expression appears to driven by native promoters. Further characterisation of disease-specific mechanisms of ABCB1 regulation are needed to inform novel approaches for overcoming chemotherapy resistance.
Availability of data and materials
Raw data files are available at the Gene Expression Omnibus with accession number GSE159512.
ABC sub-family B member 1
Acute myeloid leukaemia
Foetal bovine serum
Fragments per kilobase of transcript per million mapped reads
General transcription factor IIi
Histone H3 lysine 27 acetylation
High-grade serous ovarian cancer
Haematopoietic stem and progenitor cells
Quantitative polymerase chain reaction
Roswell Park Memorial Institute
Real-time polymerase chain reaction
RUN domain containing 3B
Solute carrier family 25 member 40
Universal Probe Library
Shaffer BC, Gillet JP, Patel C, Baer MR, Bates SE, Gottesman MM. Drug resistance: still a daunting challenge to the successful treatment of AML. Drug Resist Updat. 2012;15(1–2):62–9.
Patch AM, Christie EL, Etemadmoghadam D, Garsed DW, George J, Fereday S, et al. Whole-genome characterization of chemoresistant ovarian cancer. Nature. 2015;521(7553):489–94.
Christie EL, Pattnaik S, Beach J, Copeland A, Rashoo N, Fereday S, et al. Multiple ABCB1 transcriptional fusions in drug resistant high-grade serous ovarian and breast cancer. Nat Commun. 2019;10(1):1295.
TCGA. Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N Engl J Med. 2013;368(22):2059–74.
Nelson L, Tighe A, Golder A, et al. A living biobank of ovarian cancer ex vivo models reveals profound mitotic heterogeneity. Nat Commun. 2020;11(1):822.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Lee TI, Johnstone SE, Young RA. Chromatin immunoprecipitation and microarray-based analysis of protein location. Nat Protoc. 2006;1(2):729–48.
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60.
Splinter E, de Wit E, van de Werken HJ, Klous P, de Laat W. Determining long-range chromatin interactions for selected genomic sites using 4C-seq technology: from fixation to computation. Methods. 2012;58(3):221–30.
van de Werken HJ, Landan G, Holwerda SJ, Hoichman M, Klous P, Chachik R, et al. Robust 4C-seq data analysis to screen for regulatory DNA interactions. Nat Methods. 2012;9(10):969–72.
Kielbasa SM, Wan R, Sato K, Horton P, Frith MC. Adaptive seeds tame genomic sequence comparison. Genome Res. 2011;21(3):487–93.
Stancu MC, van Roosmalen MJ, Renkens I, Nieboer MM, Middelkamp S, de Ligt J, et al. Mapping and phasing of structural variation in patient genomes using nanopore sequencing. Nat Commun. 2017;8:1326–38.
Huff LM, Lee JS, Robey RW, Fojo T. Characterization of gene rearrangements leading to activation of MDR-1. J Biol Chem. 2006;281(48):36501–9.
Stewart A, Steiner J, Mellows G, Laguda B, Norris D, Bevan P. Phase I trial of XR9576 in healthy volunteers demonstrates modulation of P-glycoprotein in CD56+ lymphocytes after oral and intravenous administration. Clin Cancer Res. 2000;6(11):4186–91.
Bose P, Hermetz KE, Conneely KN, Rudd MK. Tandem Repeats and G-Rich Sequences Are Enriched at Human CNV Breakpoints. PLoS One. 2014;9(7):e101607.
Williams MS, Amaral FMR, Simeoni F, Somervaille TCP. Dynamic induction of drug resistance through a stress-responsive enhancer in acute myeloid leukemia. J Clin Investig. 2020;130(3):1217–32.
Mitelman F, Johansson B, Mertens F. The impact of translocations and gene fusions on cancer causation. Nat Rev Cancer. 2007;7(4):233–45.
Bailey MH, Tokheim C, Porta-Pardo E, Sengupta S, Bertrand D, Weerasinghe A, et al. Comprehensive characterization of Cancer driver genes and mutations. Cell. 2018;173(2):371–85 e18.
Bonneville R, Krook MA, Kautto EA, Miya J, Wing MR, Chen HZ, et al. Landscape of microsatellite instability across 39 Cancer types. JCO Precis Oncol. 2017. https://doi.org/10.1200/PO.17.00073.
McGranahan N, Burrell RA, Endesfelder D, Novelli MR, Swanton C. Cancer chromosomal instability: therapeutic and diagnostic challenges. EMBO Rep. 2012;13(6):528–38.
Wang ZC, Birkbak NJ, Culhane AC, Drapkin R, Fatima A, Tian R, et al. Profiles of genomic instability in high-grade serous ovarian cancer predict treatment outcome. Clin Cancer Res. 2012;18(20):5806–15.
Mantere T, Kersten S, Hoischen A. Long-read sequencing emerging in medical genetics. Front Genet. 2019;10:426. https://doi.org/10.3389/fgene.2019.00426.
This work was supported by grants from Cancer Research UK (C5759/A20971, A19842), Kay Kendall Leukaemia Fund (Clinical Training Fellowship to MW), an Educational grant from Imago Biosciences and The Christie Charity. The funding bodies had no role in: design of the study, collection of data, analysis of the data, interpretation of data, or in writing the manuscript.
Ethics approval and consent to participate
Primary human AML and HGSC samples were from the Manchester Cancer Research Centre Tissue Biobank (approved by the South Manchester Research Ethics Committee). Their use was authorized by the Tissue Biobank’s scientific sub-committee, with the written consent of donors.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Clinical details of AML samples used for targeted nanopore sequencing and ChIPseq.
Expression values (FPKM) for expressed protein coding genes.
Clinical details of HGSC samples used for targeted nanopore sequencing.
Translocations involving ABCB1 intron 1 detected by targeted nanopore sequencing of HGSC, AML and THP-1 samples.
Induction of ABCB1 in THP-1 AML cells with escalating daunorubicin exposure.
Nanopore sequencing of ABCB1high relapsed AML.
About this article
Cite this article
Williams, M.S., Basma, N.J., Amaral, F.M.R. et al. Targeted nanopore sequencing for the identification of ABCB1 promoter translocations in cancer. BMC Cancer 20, 1075 (2020). https://doi.org/10.1186/s12885-020-07571-0
- Acute myeloid leukemia
- Drug resistance
- Promoter translocation
- Ovarian cancer