LOXL2-mediated H3K4 oxidation reduces chromatin accessibility in triple-negative breast cancer cells

Oxidation of H3 at lysine 4 (H3K4ox) by lysyl oxidase-like 2 (LOXL2) generates an H3 modification with an unknown physiological function. We find that LOXL2 and H3K4ox are higher in triple-negative breast cancer (TNBC) cell lines and patient-derived xenografts (PDXs) than those from other breast cancer subtypes. ChIP-seq revealed that H3K4ox is located primarily in heterochromatin, where it is involved in chromatin compaction. Knocking down LOXL2 reduces H3K4ox levels and causes chromatin decompaction, resulting in a sustained activation of the DNA damage response (DDR) and increased susceptibility to anticancer agents. This critical role that LOXL2 and oxidized H3 play in chromatin compaction and DDR suggests that functionally targeting LOXL2 could be a way to sensitize TNBC cells to conventional therapy.


Introduction
Histone modifications contribute to gene regulation both by directly affecting chromatin structure and by recruiting effector proteins [1]. Deregulation of this enzymatic system can contribute to diseases, including cancer. The lysyl oxidase family of proteins are copper-and quinonedependent amine oxidases that oxidize the amino group located in the epsilon-position in lysines, thereby generating an aldehyde group [2]. One of the members of the LOX family, lysyl oxidase-like 2 (LOXL2), deaminates unmethylated and trimethylated lysine 4 in histone H3 (H3K4me3) through an amino-oxidase reaction that uses the Cu(II) ion and the internal cofactor lysine-tyrosylquinone, releasing the amino group and converting K4 into an allysine (H3K4ox) [3,4]. Generation of this peptidyl aldehyde likely alters the local macromolecular structure of chromatin and the nature of any protein-protein or protein-nucleic acid interactions. This is particularly relevant for gene regulation, as changes in the macromolecular status of histones can affect chromatin conformation [4][5][6].
LOXL2 is overexpressed in many tumors, and especially in breast cancers [7][8][9]. In this light, it is intriguing that some breast cancers are intrinsically resistant to chemotherapy; for these subtypes, chemotherapy induces a mesenchymal phenotype through the epithelial-to-mesenchymal transition Supplementary information The online version of this article (https:// doi.org/10.1038/s41388-019-0969-1) contains supplementary material, which is available to authorized users.
(EMT) [10]. EMT is likely to be a critical switch for tumor cell invasiveness and cell death resistance [11][12][13] and to involve chromatin reorganization, as it requires dramatic changes in cellular characteristics and gene expression [6,14]. Notably, the key transcription factor Snail1 interacts with LOXL2 [15], and LOXL2 H3K4 oxidase activity generates an H3K4ox that regulates the repression of the Ecadherin gene (CDH1) and heterochromatin transcription, which play roles in two essential steps of EMT [6,16].
Double-strand breaks (DSBs) are a major form of DNA damage and cause a specific signaling response, the DNA damage response (DDR), which can activate cell cycle checkpoint arrest and cell fate decisions, such as apoptosis or senescence. One of the first steps of DDR is the accumulation of DNA repair proteins at the damaged site [17]. Importantly, DDR activation can occur in the absence of DNA damage, by stable association of the repair factors with chromatin [18,19]. Moreover, increasing evidence suggests that higher-order chromatin structures affect DSB repair and signaling [20]. For example, DDR actively regulates decondensation of chromatin after DSBs [21], and it is amplified when chromatin is in an "open" state [20]. Similarly, DDR signaling is affected by chromatin compaction in a DNA damage-independent manner [18,19,[22][23][24].
We addressed the physiological functions of H3K4ox using an in-house generated antibody specific for this modification to analyze the H3K4ox levels in distinct breast cancer subtypes. Intriguingly, mesenchymal triple-negative breast cancer (TNBC) cell lines as well as breast cancer patient-derived xenografts (PDXs) had high H3K4ox levels that correlated with high LOXL2 expression, as compared with other subtypes. Using ChIP-seq to map its genome-wide localization, we found that H3K4ox was enriched in heterochromatin in TNBC cells, which are highly metastatic and resistant to chemotherapy. Decreasing LOXL2 levels reduced the amount of H3K4ox in chromatin, resulting in chromatin decondensation and a sustained activation of DDR. Further, both LOXL2 depletion and treatment of TNBC with chromatin-modifying drugs sensitized cancer cells to conventional treatments. Thus, targeting H3K4ox levels may open a new therapeutic window for this subtype of breast cancer.

Generating an H3K4ox-specific antibody
We initially generated a specific antibody for the recently discovered histone modification of H3K4ox, as a prerequisite for studying its physiological function. As the aldehyde group generated by LOXL2 reaction on trimethylated lysine 4 is highly reactive, and hence unfit for immunochemical studies, we hypothesized that a primary alcohol might provide a similar oxygen-bearing functionality that is less reactive, in order to generate a modificationspecific antibody that provides a readout of H3K4ox (Fig. 1a). We therefore synthesized a H3 peptide containing a 6-hydroxynorleucine residue as allysine at position 4 and used this for rabbit immunization (Fig. 1b). The resulting H3K4ox antibody was highly specific for the H3K4ox peptide, with very low cross-reactivity for unmodified H3 and no detected cross-reactivity for H3K9me3 or H3K4me3, in a wide range of experimental conditions (dot blots, western blots, and chromatin immunoprecipitation (ChIP) experiments) (Fig. 1c-f). Analysis of purified nucleosomes from 293T cells showed that H3K4ox levels increased, and H3K4me3 levels decreased, when nucleosomes were incubated with wild-type (wt) recombinant LOXL2 but not with a catalytically inactive LOXL2 (LOXL2mut) [3] (Fig. 1d, upper panel). Moreover, the levels of H3K4ox also increased in MCF-7 cells transfected with LOXL2 as compared with cells transfected with the empty vector (Fig. 1d, lower panel). Finally, MDA-MB-231 cells infected with an shRNA targeting the human LOXL2 (LOXL2 knockdown (KD)) showed a specific reduction in H3K4ox levels as compared with cells infected with an irrelevant shRNA (control), in both western blots and ChIP-PCR experiments using the E-cadherin gene promoter (CDH1), which is a well-known LOXL2-mediated H3K4 oxidation target promoter (Fig. 1e) [3,4]. Kinetics of the reaction using recombinant LOXL2 and nucleosomes revealed that levels of intermediate alcohol were maintained for 2 h, after which they were reduced (Fig. 1f). Biotinylated hydrazide was used to detect the generated aldehyde group, which appears after 4 h of reaction (Fig. 1f). Thus, as the intermediate alcohol is relatively stable, the antibody we generated can be used as a readout of the oxidized histone H3K4, although we cannot distinguish whether it recognizes the intermediate alcohol, the aldehyde group, or both.
To elucidate the function of H3K4ox in breast cancer cells, we first performed a ChIP-seq experiment with the Fig. 1 Quality control of the anti-H3K4ox antibody. a Schematic representation of the LOXL2 reaction. The red circle indicates the intermediate residue that is targeted by the in-house generated anti-H3K4ox antibody. b The artificial amino acid 6-hydroxynorleucine was used in the peptide to generate the anti-H3K4ox antibody. c The anti-H3K4ox antibody was found to be specific in western blot in two replicates of dot blots of dilution series of oxidized histone H3 peptide (H3K4ox) or unmodified H3 peptide (left panel), as well as in a representative dot blot of a dilution series of H3K9me3, H3K4ox, H3K4me3, or H3 peptides (right panel). d Nucleosomes were incubated with recombinant wild-type (wt) LOXL2 or a catalytically inactive LOXL2 (mut) purified from baculovirus to detect H3K4ox/ H3K4me3 levels (upper panel). Lysates of MCF-7 cells transfected with an empty vector (ø) or with LOXL2 were analyzed by western blotting, using the indicated antibodies (lower panel). e Western blot for LOXL2, H3K4ox, and total H3 from MDA-MB-231 cells infected with short hairpin RNA (shRNA) as a control, or a knockdown (KD) using a shRNA specific for LOXL2 (LOXL2 KD) (left panel). Anti-H3K4ox ChIP-PCR was used to analyze the E-cadherin gene (CDH1) promoter in MDA-MB-231 cells infected with shRNA for either control (green bar) or LOXL2 KD (blue bar). Data of qPCR amplifications were normalized to the input and to total H3 for each condition. Error bars indicate the SD from at least three independent experiments. **P < 0.01. f Western blot of H3K4ox and biotin incorporation (BTH-H3) in nucleosomes incubated with recombinant LOXL2 purified from baculovirus, after different incubation times anti-H3K4ox antibody to determine the genomic distribution of H3K4ox in the MDA-MB-231 cells. Peaks called using model-based analysis for ChIP-seq (MACS) [26] showed low differences in H3K4ox between two sequencing replicates, with a genome-wide Pearson correlation coefficient of the read count of the two replicates of 0.997 (Fig. 2d, upper panel). We observed that H3K4ox peaks were distributed throughout different genomic elements (Fig. 2d, lower panel). Using the ChromHMM tool [27], we assessed the statistical overrepresentation of the H3K4ox peaks through different chromatin states (promoter, weak or strong enhancer, insulator, repressed, or heterochromatin [nonrepetitive sequences]) and found that H3K4ox peaks were significantly overrepresented in heterochromatin (Fig. 2e, f). As generating an aldehyde in H3 removes a positive charge and creates a very reactive group, we hypothesized that this reaction affects chromatin structure. To test this, we used the assay for transposase-accessible chromatin [28,29] followed by deep sequencing (ATAC-seq), which exploits the ability of the prokaryotic transposase Tn5 to integrate preferentially into accessible (open) chromatin. ATAC-seq showed an increased ATAC signal in LOXL2 KD cells, but not with control cells, at H3K4ox-positive sites (Fig. 2g). These results were validated in selected regions by ChIP-qPCR and ATAC-qPCR in control and LOXL2 KD cells: H3K4ox enrichment decreased in LOXL2 KD cells, with a correlating increase of ATAC signal, in these regions (Fig. 2h). No changes were observed in a control promoter, POL2RA (Fig. 2h). These data demonstrated that, in the absence of LOXL2, H3K4ox levels decrease and chromatin adopts a more open conformation (Fig. 2g). Thus, our results showed that H3K4ox is enriched in heterochromatin and is directly linked with chromatin accessibility in those regions.
Chromatin structure alterations activate DDR in a LOXL2-dependent manner As the chromatin state can influence many aspects of DDR [30], we hypothesized that disruption of LOXL2 expression and impairment of H3K4ox generation might influence DDR by affecting chromatin accessibility. To test this, we analyzed by immunofluorescence MDA-MB-231 cells that had been infected with either LOXL2 KD or control lentiviruses, using two well-established markers of DDR: phosphorylated H2AX (γ-H2AX) and TP53-binding protein 1 (53BP1). Depletion of LOXL2 (using LOXL2 KD) led to more foci of both γ-H2AX and 53BP1 than in control cells, suggesting that LOXL2 KD cells may accumulate DNA breaks and/or activate DDR (Fig. 3a). To determine if the LOXL2 catalytic activity was involved in the observed phenotype, LOXL2 KD cells were complemented by reinfection with ectopic vector expressing either the wt LOXL2-IRES-GFP or a catalytically inactive LOXL2 (LOXL2mut-IRES-GFP), both of which were expressed at similar levels ( Fig. S1). Fewer γ-H2AX and 53BP1 foci were observed in LOXL2 KD cells after reintroduction of the wt (but not of the catalytically inactive) LOXL2 (Fig. 3b), establishing that suppressing DDR activation requires both the activity of LOXL2 and H3K4ox generation.

LOXL2 KD activates DDR independently of DNA damage
Increased DDR activation in LOXL2-depleted cells could be a consequence of more DSBs due to reduced H3K4ox levels and/or chromatin decondensation. Notably, aberrant silencing and conflicts between replication forks and transcription, as well as the presence of R-loops (a transcriptional intermediate), can result in DNA damage and are influenced by chromatin state, for example in cells lacking the linker histone H1 [31]. However, our analysis of RNAseq data revealed that LOXL2 KD cells did not have altered expression of repetitive elements (Table 1), and DDR activation in LOXL2 KD cells was not affected by cordycepin, an inhibitor of RNA synthesis that abolishes R-loop formation [32,33] (Fig. 3c). As R-loops can generate DNA damage during replication due to fork stalling and collapse [32,34], we next analyzed γ-H2AX and 53BP1 foci in noncycling LOXL2 KD cells and found DDR activation to also be increased (Fig. 3d). Overall, these data suggested that overexpression of repetitive elements, R-loop formation, and replication fork stalling were not responsible for activating DDR following LOXL2 depletion, when heterochromatin adopts a more open state. As chromatin structure alterations can trigger the DDR, even in the absence of DNA damage [22,35,36], we checked for the presence of DNA damage following LOXL2 depletion more directly, using the comet assay, in LOXL2 KD or control cells under alkaline conditions. No increases in DNA damage (due to either single-strand or double-strand DNA breaks) were observed in LOXL2 KD cells (Fig. 3e). Moreover, we did not observe any significant differences in chromosomal lesions between LOXL2 KD or control cells in metaphase spreads (Fig. 3f). Finally, analyzing for mitotic aberrations (anaphase bridges and micronuclei), which can be indicative of replication stress or DNA repair defects, we observed only a mild increase in anaphase bridges in the absence of LOXL2 (Fig. S2). Taken together, our data suggested that the combination of loss of LOXL2 and reduced H3K4ox levels in TNBC cells was sufficient to activate the DDR in the absence of detectable DNA lesions.

Alterations in chromatin compaction activate the DDR
To further address the origin of DDR signaling in LOXL2depleted cells, we analyzed the behavior of additional DDR signaling components. As both γ-H2AX and 53BP1 foci formation require the ATM kinase in some settings, we treated LOXL2 KD and control cells with the ATM inhibitor KU55933 and analyzed foci formation (Fig. 4a). Decreased foci of both markers were observed upon ATM inhibition, indicating that LOXL2-induced DDR was largely ATMdependent. Consistent with this, LOXL2 KD cells had increased phosphorylation of several ATM substrates, including KAP-1, CHK1, and CHK2 (Fig. 4b). However, we ruled out that increased DDR signaling was due to apoptosis in LOXL2 KD cells, as no cleaved caspase-3 signal was observed in either LOXL2 KD or control cells (Fig. S3a). As these data suggested that the LOXL2 KD cells activated a checkpoint response, we analyzed cell cycle progression following LOXL2 depletion. After synchronization with a double thymidine block, LOXL2 KD cells were not able to efficiently progress through the cell cycle (Fig. 4c), and western blotting for H3S10-P showed that this histone mark was undetectable in LOXL2 KD cells as compared with control cells. (Fig. 4d). These data strongly suggested that LOXL2 KD cells arrested primarily in G1. Consistent with this possibility, cell proliferation capacity of LOXL2 KD cells was blocked (Fig. 4e, upper panel), and their colony-formation capacity was strongly reduced after only a few passages (Fig. 4e, lower panel). This effect on proliferation was further confirm by gene ontology analysis of the gene expression pattern of LOXL2 KD cells ( Fig. S3b; Tables 2 and 3).
To test whether the effects on chromatin in LOXL2 KD cells directly activated ATM-dependent DDR signaling, we forced chromatin condensation in these cells by expressing the linker histone H1 or the H3K9 methyltransferase SUV-39H1 (as green fluorescent protein (GFP)-labeled proteins) and counted the number of γ-H2AX foci in GFP-positive cells. Notably, overexpression of either H1 GFP or SUV-39 GFP in LOXL2 KD cells reduced the number of γ-H2AX foci as compared with LOXL2 KD cells that expressed GFP alone (+MOCK GFP ) (Fig. 4f, g). This suggested that lack of LOXL2 and reduced H3K4ox levels affected the regulation of chromatin condensation (leading to decondensation) and activated DDR, even in the absence of DNA damage.

LOXL2 reduction enhances chemosensitivity of TNBC cells
We found that reducing H3K4ox levels via LOXL2 depletion led to chromatin decondensation, triggered DDR activation, and induced cell cycle arrest, suggesting that LOXL2 inhibition could be interesting as a breast cancer treatment. As no LOXL2-specific inhibitors are currently available, we therefore tested whether reducing the levels of functional LOXL2 would increase the sensitivity of breast cancer cells to chemotherapy. For this, we treated three different TNBC cell lines with doxorubicin (a topoisomerase inhibitor that generates DSBs) either alone or together with trichostatin A (TSA) (a general HDCA inhibitor [37] that increases chromatin accessibility [38]). Indeed, we observed that in all cases, treatment of doxorubicin together with TSA increased the percentage of cell death-indicative of an increased cell sensitivity to chemotherapy (Fig. 5a). Similar results were obtained with two TNBC_PDXs (PDX-549 and PDX-154) (Fig. 5a). To determine if these results were reproducible in vivo, PDX-549 cells were first grown ex vivo, and then 10 6 cells were subcutaneously implanted in nude mice. After tumor formation, mice were treated with TSA, doxorubicin, doxorubicin plus TSA, or (as a control) vehicle for 25 days. In agreement with our previous results, we observed that tumor growth was substantially impaired with the combined TSA/doxorubicin treatment (Fig. 5b).

Discussion
In this work, we found that several TNBC cell lines and PDXs express high levels of LOXL2 as compared with those of the luminal breast cancer subtype. With the use of a newly generated antibody, we were able to detect the H3K4ox modification produced by LOXL2 [3] and we observed how this histone H3 modification was higher in all TNBC cell lines and PDXs as compared with luminal subtypes. This H3 modification is enriched in heterochromatin and is required to maintain the condensed heterochromatin state. It is possible that during tumor progression, some cancer cells undergo the EMT program and start to express LOXL2. The transcription factor Snail1, together with LOXL2, would participate in downregulating the CDH1 gene and other heterochromatin transcripts, giving rise to transformation of cancer epithelial cells into mesenchymal cells (through EMT) [6,[39][40][41]. Recently, two different groups have suggested that EMT is dispensable for lung and pancreas metastasis but contributes significantly to the formation of recurrent metastasis after chemotherapy [42,43]. This finding is in agreement with our results, in which we observed that the most aggressive and resistant subtype of breast cancer (TNBC) has a  Contingency table of the  Fisher's exact test showed the statistical overrepresentation of the  H3K4ox peaks through different chromatin states. f Genome browser view of H3K4ox and H3K9me3-binding profiles at two representative heterochromatin regions. g Heatmaps show the ATAC signal in all peaks as well as in peaks that overlap with H3K4ox in LOXL2 KD or control cells. h H3K4ox ChIP-PCR validation of the selected genomic regions from the ChIP-seq from control or LOXL2 KD MDA-MB-231 cells. Data of qPCR amplification were normalized to the input and to total H3 (upper panel). ATAC-qPCR validation of the incorporation of the transposase Tn5 at the selected genomic regions from the ChIP-seq from control or LOXL2 KD MDA-MB-231 cells. Data of qPCR amplification were normalized to an unchanging genomic region (the HPRT promoter) and expressed as the fold-change relative to data obtained from control cells, which were set to 1 (lower panel). In both experiments, the RNA polymerase II (POL2RA) promoter was used as a negative control. Error bars indicate the SD from at least three independent experiments. *P < 0.05, **P < 0.01, ***P < 0.001 mesenchymal phenotype and is enriched for H3K4ox in its heterochromatin. We observed that reducing H3K4ox levels perturbed the balance between condensed and decondensed chromatin and activated the DDR. Similar phenotypes have been previously observed after depletion of SUV-39, PRDM2, or HP1, which are required for heterochromatin maintenance and cell survival after DNA damage [23,24,44,45]. H3K4ox levels have been correlated to the repression of CDH1 and major satellites (heterochromatin transcripts) during the EMT program [3,4,6]. CDH1 is upregulated in TNBC cells when LOXL2 is knocked down and the levels of H3K4ox are decreased; however, the H3K4ox modification does not seem to be involved in repression of heterochromatin transcripts, as we did not observe increased levels of heterochromatin transcripts in LOXL2 KD conditions. Notably, downregulation of heterochromatin transcripts during EMT occur at a specific timepoint during chromatin reorganization, after which the levels of transcription recovered despite the presence of H3K4ox. Intriguingly, although upstream DDR signaling is sustained when chromatin is maintained in a condensed state [22,46,47], we found that forcing chromatin decondensation was sufficient to activate DDR signaling in the absence of any detectable DNA lesions. It seems that an abnormal, unbalanced ratio between condensation and decondensation leads to persistent DDR signaling. As LOXL2 has multiple substrates besides histones [5], we cannot discard an additional putative role for LOXL2 in the direct regulation of the DDR machinery. However, we observed that the DDR induced in the absence of LOXL2 can be inhibited by forcing chromatin to condense using either of two different approaches, suggesting that the main molecular mechanism is the regulation per se of the chromatin compaction state.
While a correlation between chromatin compaction and DDR has been previously reported [22,46,47], we demonstrated here that H3 oxidation by LOXL2 is another molecular mechanism that maintains compaction. Notably, LOXL2 has been widely linked to cancer and the acquisition of cellular malignancy, as it is overexpressed in many tumors [15,[48][49][50][51] and has an important role in tumor formation [52].
Induction of chromatin compaction has been suggested as a potential therapeutic tool in gliomas, in which chromatin compaction limits DDR [46] and promotes damage resistance [53]. However, this does not appear to apply to TNBC cells. In support of this, we have now shown that several TNBC PDX tumors treated with TSA, a histone deacetylase inhibitor that generates an open chromatin state, enhanced tumor killing with doxorubicin in these cells. These observations suggest a rationale for studying combinations of "chromatin opening drugs," including inhibitors of histones deacetylases, EZH2, or LSD1, as a strategy for overcoming resistance to treatment in TNBC. ) as well as cells isolated from PDXs, were maintained in Dulbecco's modified Eagle's medium (Biowest; L0106-500); the TNBC line BT-549 (ATCC HTB 122) was maintained in Roswell Park Memorial Institute medium (Biowest; L0501-500). All cell cultures were supplemented with 10% fetal bovine serum (Gibco; 10270106), 2 mM Lglutamine (Biowest; X0550-100), and 1% penicillinstreptomycin (Gibco; 15140122) at 37°C in 5% CO 2 . All cell lines were regularly tested for the absence of mycoplasma using standard polymerase chain reaction (PCR). For LOXL2-Flag overexpression assays, MCF-7 cells were seeded for 24 h and transfected with 10 µg pcDNA3-hLOXL2-Flag vector using polyethylenimine polymer (Polysciences Inc; 23966-1). Lentiviral particles were produced in HEK293T cells (ATCC CRL-3216) for lentiviral infections to produce LOXL2 KD cells. Cells were grown to Fig. 3 Chromatin opening activates DDR via a mechanism that requires catalytically functional LOXL2 but can be independent of DNA damage. a γ-H2AX and 53BP1 staining and foci quantification are shown by immunofluorescence with a specific antibody for γ-H2AX (left image) or for 53BP1 (right image). Dot graphs indicate the number of foci for γ-H2AX (upper graph) and 53BP1 (lower graph) per cell in control and LOXL2 KD conditions. b γ-H2AX staining and foci quantification are shown by immunofluorescence with the indicated antibody after LOXL2 reinfection. MDA-MB-231 cells were first infected with control or LOXL2 KD lentivirus and then, after puromycin selection, again with GFP (MOCK GFP+ ), wild-type LOXL2-IRES-GFP (LOXL2 GFP+ ), or LOXL2mut-IRES-GFP (LOXL2mut GFP+ ). Cells were fixed after 24 h. Dot graphs indicate the number of γ-H2AX foci per GFP-positive cells containing MOCK GFP+ (upper graph), LOXL2 GFP+ (middle graph), or LOXL2mut GFP+ (lower graph). c Dot graphs indicate the number of γ-H2AX (upper graph) and 53BP1 (lower graph) foci per cell in control and LOXL2 KD cells after treatment with 200 μM cordycepin for the indicated timepoints. Error bars indicate the SD from at least three independent experiments. *P < 0.05, **P < 0.01, and ***P < 0.001, ****P < 0.0001 70% confluency and then transfected by dropwise addition of NaCl (150 mM), a DNA mixture of pLKO-short hairpin control (shCT)/shLOXL2 (7.5 µg), pCMV-VSVG (1.5 µg), pMDLg/pRRE (4.5 µg), and pRSV-Rev (1.5 µg), and polyethylenimine polymer (1 mg/ml) (Polysciences Inc), which had been preincubated for 15 min at room temperature. Transfection medium was replaced with fresh medium after 24 h (day 1). The supernatant from the transfection medium at days 2 and 3 was filtered with 0.45 μm filter unit (Merck Millipore) and stored at 4°C (both were mixed after day 2). The supernatant mixture was concentrated using Lenti-X Concentrator product (Clontech) and centrifuged at 1500 × g for 45 min at 4°C. The pellet was resuspended in 1 ml fresh medium, aliquoted (100 μl), and stored at -20°C. Aliquots were used to infect MDA-MB-231 cells [6]. For retroviral infections, HEK293 gag-pol cells were used to produce retroviral particles. Cells were transfected as for HEK293T cells with a mixture of NaCl (150 mM), DNA (2.5 µg of pCMV-VSV-G and 7.5 µg of pMSCV, pMSCV-LOXL2 wt-FLAG or pMSCV-LOXL2 mutFLAG IRES-GFP vectors) and polyethylenimine polymer (1 mg/ml) (Polysciences Inc; 23966-1) that had been preincubated for 15 min at room temperature; transfection  medium was replaced with fresh medium at 24 and 48 h and processed as explained above. Viral particles were concentrated using Retro-X Concentrator product (Clontech; 631456) and then used to infect MDA-MB-231 cells.

Antibodies and other reagents
The ). An anti-H3K4ox antibody was generated from a synthetic peptide in which fluorenylmethoxycarbonyl protecting group (Fmoc)-6-hydroxynorleucine (BAA1117, Iris Biotech) rather than Fmoc-lysine(Boc) was incorporated at position 4. The peptide was purified (>95%) by reversedphase high-performance liquid chromatography and its identity was confirmed by mass spectrometry, after which it was coupled to keyhole limpet hemocyanin (KLH) for antibody production. Antibody information can be found in the Supplementary Table 1.

Cell extracts, histone isolation, and PDX extraction
To obtain nuclear fractions of LOXL2-Flag-transfected MCF-7 cells and HEK293T cells, cells were first lysed in soft-lysis buffer (50 mM Tris, 2 mM EDTA, 0.1% NP-40, and 10% glycerol, supplemented with protease and phosphatase inhibitors) for 5 min on ice. Samples were centrifuged at 900 × g for 15 min, and the supernatant was discarded. The nuclear pellet was lysed in high-salt lysis buffer (20 mM HEPES pH 7.4, 350 mM NaCl, 1 mM MgCl 2 , 0.5% Triton X-100, and 10% glycerol, supplemented with protease and phosphatase inhibitors) for 30 min at 4°C, and samples were centrifuged at 16,000 × g for 10 min supernatant NaCl concentration was reduced to 300 mM NaCl by adding balance buffer (20 mM HEPES pH 7.4, 1 mM MgCl 2 , and 10 mM KCl). Cells extracts from the different breast cancer cell lines were obtained with SDS lysis buffer (2% SDS, 50 mM Tris-HCl, and 10% glycerol). Histones were isolated from the different breast cancer cell lines using acid precipitation. For this, cell monolayers were first scraped in 1 ml of lysis buffer (10 mM Tris pH 6.5,

ATAC-seq and ATAC-qPCR experiments
ATAC experiments were performed as described previously [28]. Cells were harvested and treated with transposase Tn5 (Nextera DNA Library Preparation Kit, Illumina; FC-121-1030). DNA was purified using MinElute PCR purification kit (Qiagen), samples were amplified by PCR using NEBNext High-Fidelity 2× PCR Master Mix (New England Biolabs; M5041), and DNA was again purified with the MinElute PCR purification kit. Reads produced by ATAC sequencing of two control (Control) replicates and two LOXL2 KD sequencing replicates (LOXL2 KD) were aligned to the hg19 build of the reference human genome using Bowtie 2 [56] with default parameters for pair-end sequencing. ATAC peaks were then called by combining aligned reads of both replicates of the control and the KD using MACS2. To allow for false discovery rate (FDR) threshold selection further downstream in the analysis, no FDR restrictions were imposed on the ATAC peak calling. For ATAC-qPCR experiments, the final elution product of the ATAC protocol was diluted 1:50. Incorporation of the transposase Tn5 to GRs was detected by quantitative staining with PCR SYBR Green (Quantabio; 95073). ATAC-qPCR results for the selected GRs were normalized to the same unaffected GR (HPRT promoter) for each condition. The sequences of the primers used can be found in Supplementary Table 2.

RNA-seq analysis
Reads produced by RNA-seq of the same replicates as described above of two controls (Control) and two LOXL2 KDs were aligned to the hg19 build of the reference human transcriptome using TopHat2 [57] with default parameters for pair-end sequencing. Aligned reads were then analyzed using a standard Cufflinks pipeline [58] to detect differentially expressed genes between the two conditions (LOXL2 KD and Control).

Differential expression analysis of transposable elements (TEs)
The trimmed RNA-seq reads from the two controls and the two LOXL2 samples were processed with the TEtools program [59] using the hg19 annotation of TEs. The obtained count table was then processed with NOIseq [60] to perform a differential expression analysis of LOXL2 against the control. Expression was normalized using the TMM method, whereby TEs with a probability higher than 0.95 of being differentially expressed were considered as statistically significant.

Replicate correlation
The read count (coverage) at each position of the hg19 human reference genome was computed for each replicate of the H3K4ox ChIP-seq and the control and LOXL2 KD ATAC sequencing, using the BEDtools genomecov capability [61]. Genomic positions with zero-read counts were filtered out. Replicate files of each experiment were merged to produce a single file aligned by genomic position, and the corresponding Pearson's correlation coefficient of read counts was computed. For the graphical representation of the correlation, 100,000 genomic positions were randomly selected.
Analysis of ATAC peaks that overlap with H3K4ox peaks All significant (P < 10 -5 ) ATAC peaks (LOXL2 KD versus control) and H3K4ox peaks were first intersected with the BEDtools intersect program [61]. Based on this intersection, ATAC peaks were classified as overlapping (if they intersected an H3K4ox peak) or orphan (if not). Only intersections involving more than 95% of the sequence of ATAC peaks were considered. Control and LOXL2 KD read counts over all genomic positions (see above) were intersected with both overlapping and orphan peaks. Read counts over genomic positions of control and experiment replicates were averaged. To carry out the heatmap representation, peak sequences (overlapping or orphan) were aligned by their summits. For linear representation, the average experimental read counts at each downstream and upstream position were summed for both the experimental and the control counts. Position-wise sums were then divided by the read count sum value obtained for the summit of control read counts, thus making all sums relative to the maximum control value.

Integrated analysis of H3K4ox and ATAC peaks, and differentially expressed genes
The differentially expressed genes detected through the RNA-seq analysis of control and LOXL2 KD cells were selected if they were in close vicinity (either up-or downstream) to overlapping ATAC peaks. Two different distance thresholds (0.5 and 1 Mb) were used to detect close differentially expressed genes.
Inhibition of RNA synthesis, ATM kinase, and cell cycle analysis MDA-MB-231-infected cells selected for 48 h with puromycin (2.5 μg/ml) were seeded and then maintained under puromycin selection.
For inhibition of RNA synthesis, cells were treated at 48 h with 200 µM cordycepin (Sigma; C3394) and fixed with 100% cold ethanol at designated timepoints.
For inhibition of the ATM kinase, cells were incubated at the 24-h timepoint with KU55933 (Sigma; SML1109) for a further 24 h at two different concentrations (either 5 or 10 µM), or for 2 h with doxorubicin (1 µM) as a positive control.
For cell cycle analysis, cells at 48 h after puromycin selection were first synchronized through the double thymidine block protocol. Specifically, cells were seeded to 50% confluency, incubated for 14 h with complete growth medium (supplemented with 2 mM thymidine), washed 2× with phosphate-buffered saline (PBS), and released by a 9-h incubation with complete medium growth. Cells were washed again 2× with PBS, incubated 14-h with 2-mM thymidine, and released with complete growth medium. Cells were then harvested at the designated timepoints and fixed with 100% cold ethanol. After two days, fixed cells were stained with propidium iodide (PI) and analyzed by flow cytometry using BD FACSCalibur (Becton Dickinson). Results were analyzed using BD CellQuest Pro software.

Non-replicative cell experiment
MDA-MB-231 cells were seeded in coverslips and maintained during all the experiment in Dulbecco's modified Eagle's medium (Biowest; L0106-500) with 0.5% fetal bovine serum (Gibco; 10270106) at 37°C in 5% CO 2 . After 24 h, cells were infected with lentiviral particles for LOXL2 KD. After 96 h under selection, cells were fixed with paraformaldehyde 4% with PBS.

Comet assay
MDA-MB-231-infected cells under puromycin selection (see above) were seeded at the 48-h timepoint and grown another 48 h more (still with puromycin at 2.5 µg/ml). DNA strand breaks were measured at the single-cell level for MDA-MB-231 control or LOXL2 KD cells using an alkaline comet assay with the CometAssay kit (Trevigen, Gaithersburg, MD) following the manufacturer's instructions. Briefly, 5000 cells in 50 μl PBS were combined with 500 μl molten LMAgarose (at 37°C), and 50 μl of this was immediately transferred to a CometSlide. After a 10-min incubation at 4°C, slides were immersed in 4°C lysis solution and incubated overnight at 4°C. Slides were then immersed in freshly prepared alkaline solution (200 mM NaOH, 1 mM EDTA pH > 13) for 1 h at 4°C in the dark. For electrophoresis, slides were placed in an electrophoresis slide tray containing alkaline solution and incubated with voltage (21 V) for 30 min. After that, slides were immersed twice in H 2 O for 5 min and then once in 70% ethanol for 5 min. Samples were dried at 37°C and stained with SYBR Green staining solution (SYBR Green I, Invitrogen) using a 1:10000 dilution in TE (10 mM Tris-HCl pH 7.5, 1 mM EDTA) for 30 min Finally, slides were dried at 37°C, and cells were imaged using an Olympus BX61 microscope.

Rescue experiments
MDA-MB-231-infected cells under puromycin selection (see above) were seeded at the 48-h timepoint, incubated under puromycin selection for an addition 24 h, and then transfected with the SUV-39H1-EGFP vector or reinfected with retroviral particles for LOXL2-FLAG or lentiviral particles for H1 expression. After 24 h, cells were fixed for immunofluorescence.

Immunofluorescence, image acquisition, and analysis
Cells were fixed with 4% paraformaldehid for 15 min at room temperature, blocked for 1 h with 1% PBS-bovine serum albumin, incubated at room temperature for 2 h with primary antibody, washed 3× with PBS, and then incubated for 1 h at room temperature with the secondary antibody. Cells were washed again 3× with PBS, incubated for 5 min with 4′,6diamidino-2-phenylindole (DAPI) (0.25 mg/ml) for cell nuclei staining, and then mounted with fluoromount. Fluorescence images corresponding to DAPI, GFP, γ-H2AX and 53BP1 were acquired in a Leica TCS SPE microscope using a Leica DFC300 FX camera and the Leica IM50 software.

Cellular viability experiment
Breast cancer cell lines and PDX cells were seeded in 96-well plates. Cellular viability was analyzed using Thiazolyl Blue Tetrazolium Bromide (MTT) (Sigma; M5655) at different timepoints. The concentrations used were: 0.1 µM doxorubicin for all samples; 250-nM TSA for the breast cancer cell lines; and 500-nM TSA for the PDXs cells. Absorbance was detected at 565 nm on Infinite® 200 PRO Series Multimode Reader (Tecan Group Ltd.) and analyzed with i-control™ Microplate Reader Software (Tecan Group Ltd.)

Cancer PDXs and treatments
Samples from patients with breast cancer were obtained from the operating room and transferred to the pathology department, where breast cancer samples were collected, transferred to the animal facility, and implanted into mice, to generate PDXs. All samples were implanted within 60 to 90 min after surgical removal. All patients willingly signed an informed consent, and the study was approved by the Ethics Committee of the Vall d'Hebron Hospital.
For in vivo experiments, a tumor from an established PDX was dissociated into single cells by enzymatic digestion (collagenase at 300 U/ml and hyaluronidase at 200 U/ ml) during 1 h at 37°C on a rotating wheel. The solution of digested tumor was treated with 0.025% trypsin and filtered sequentially using 100 and 40 µm strainers. Isolated cells were plated in culture dishes with DMEM-F12 supplemented with FBS, glutamine, and penicillin/streptomycin. Once the culture was established, 10 6 cells were injected into the number four fat pad of 6-week old NOD. CB17PrkdcSCID/J (NOD/SCID) female mice (Charles River) with Matrigel. For this, mice were anesthetized and shaved, the fourth and fifth sets of nipples were localized, and an inverted Y incision was made from the midline point between the fourth set of nipples, ending between the fourth and fifth sets to expose the fourth and fifth fat pads on one side. After the injection, animals were sutured, and analgesics injected. Animals were kept in a clean cage with drinking water supplemented with 1 µM 17-β-estradiol. Tumor xenografts were measured with callipers every 3 days, and tumor volume was determined using the formula: (length × width 2 ) × (pi/6). At the end of the experiment, animals were anesthetized with a 1.5% isoflorane-air mixture and were killed by cervical dislocation.
Treatments were administered intraperitoneally twice weekly. One week after cell injection, mice were randomized and treated with TSA (0.25 mg/kg), doxorubicin (2 mg/kg mouse weight), or doxorubicin plus TSA. The control group was injected with sterile PBS.
Mice were maintained and treated in accordance with institutional guidelines of Vall d'Hebron University Hospital Care and Use Committee.

GSE96064.
funding from the MINECO (a "Juan de la Cierva Incorporation" fellowship; IJCI-2014-20723). SP was a recipient of a Miguel Servet contract (ISCIII/FIS), and AI, JPC-C, LP-G, and GS-B are supported by contracts from Worldwide Cancer Research, Fundació La Marató TV3, Fundació FERO, and a FI Fellowship from the Generalitat de Catalunya, respectively.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons. org/licenses/by/4.0/.