The present invention relates to biomarkers which comprise one or more genomic sequence(s) comprising epigenetic modification and their use in a method for predicting resistance to or assessing possible outcomes of cancer treatment, in particular for patient stratification.
The emergence of resistance to chemotherapy and targeted therapies is a major challenge for the treatment of cancer. Genetic heterogeneity within untreated tumors is now considered to be a key determinant of resistance; sub-population of cells bearing a mutation conveying resistance can survive and be selected in a Darwinian process (Schmitt, M. W. et al. 2016. Nat Rev Clin Oncol 13, 335-347). Deep sequencing and single-cell approaches have revealed the importance of genetic intra-tumor heterogeneity to tumor evolution (Roth, A. et al. 2016. Nat Methods 13, 573-576; Nik-Zainal, S. et al. 2012. Cell 149, 994-1007; McGranahan, N. & Swanton, C. 2017. Cell 168, 613-628) and shown that genetic heterogeneity within untreated tumors is a key factor in tumor resistance (Dagogo-Jack, I. & Shaw, A. 2018. Nat Rev Clin Oncol 15, 81-94).
However, in many cases genetic mechanisms driving resistance cannot be found, pointing to a role for non-genetic mechanisms (Salgia, R. & Kulkarni, P. 2018. Trends Cancer 4, 110-118). Non-genetic and particularly transcriptional and epigenetic mechanisms are anticipated to play a role in the adaptation of cancer cells confronted with environmental, metabolic or therapy-related stresses (Rathert, P. et al. 2015. Nature. 525, 543-547; Kim, C. et al. 2018. Cell 173, 879-893 e813).
Recent studies, using single-cell RNA sequencing (scRNA-seq), indicate that emergence of transcriptional subclones upon treatment may account for the adaptation of cancer cells to therapeutic pressure (Kim, C. et al. 2018. Cell 173, 879-893 e813; Horning, A. M. et al. 2018. Cancer Res 78, 853-864).
In contrast, only a few studies have tracked the clonal evolution of epigenetic alterations, exclusively analyzing DNA methylation at the population level (Mazor, T. et al. 2015. Cancer Cell 28, 307-317; Aryee, M. J. et al. 2013. Sci Transl Med 5, 169ra110) suggesting that DNA methylation alterations and genetic mutations shared a common evolutionary track.
Modulation of chromatin structure via histone modification is a major epigenetic mechanism and key regulator of gene expression. However, the contribution of chromatin heterogeneity to tumor evolution remains unknown.
Profiling histone modifications at single-cell resolution remains challenging, in part because the level of noise associated with non-specific binding during the immunoprecipitation tends to increase with low quantity of starting material. Immunoprecipitating chromatin from one single cell is not technically feasible and one therefore needs to tag the chromatin of each single-cell prior to immuno-precipitation, before pooling the chromatin fragments of several thousand cells and performing immuno-precipitation.
Until now, insufficient coverage has limited the applications of single-cell chromatin profiling to cell lines (Rotem, A. et al. 2015. Nat Biotechnol 33, 1165-1172; WO 2013/134261), preventing the study of the heterogeneity of chromatin states in complex biological systems such as tumors.
Using an improved single-cell ChIP method based on droplet microfluidics to profile chromatin landscapes of thousands of cells at single-cell resolution with a coverage of 10,000 loci/cell, the inventors detected the presence of relatively rare chromatin states within tumor samples. In patient-derived xenografts models of acquired resistance to chemotherapy and targeted therapy in breast cancer, the inventors found that respectively 3% and 16% of cells in the untreated tumors possessed chromatin features characteristic of resistant cancer cells. These cells, and cells from the resistant tumors, had lost chromatin marks (H3K27me3) on specific genomic sequences. Some of which are associated with stable transcriptional repression for genes known to promote resistance to treatment, potentially priming them for transcriptional activation.
In certain embodiments, the present invention relates to a biomarker for determining resistance to treatment of a cancer type with a cancer drug which comprises one or more genomic sequence(s) comprising a histone modification, wherein said one or more genomic sequence(s) is selected from the list of Tables 1 to 6. In particular, said histone is in a gene or in the proximity of said gene and said genomic sequences is selected from the list of Table 1 and 4. Preferably, the invention relates to a biomarker as described above for predicting resistance to treatment of a cancer type with a cancer drug prior administration of a treatment to a patient. In another aspect, the invention provides a method for determining the resistance to treatment of a cancer type with a cancer drug to a patient comprising: i) detecting a histone modification of at least one biomarker as described above which comprises one or more genomic sequence(s) selected from the list of Tables 1 to 6 in said patient tumor sample, and ii) determining from the presence or absence of histone modification of said biomarker, whether the patient is likely to be resistant or sensitive to said treatment. In certain embodiments, the present invention also provides a method for determining the resistance to treatment of a cancer type with a cancer drug to a patient comprising: i) determining in said patient tumor sample the expression level of a gene encoding by a biomarker as described above which comprises one or more genomic sequence(s) selected from the list of Tables 1 and 4, ii) determining from the expression level of said gene whether a patient is likely to be resistant or sensitive to said treatment.
Preferably, said method is realized prior to administration of any treatment or said treatment to the patient. Particularly, said cancer drug agent is a chemotherapy drug, preferably capecitabine and said genomic sequence is selected from the list of Tables 1 to 3. In another aspect, said cancer drug agent is an anti-hormonal drug, preferably tamoxifen and said genomic sequence is selected from the list of Tables 4 to 6.
Preferably, said histone modification is associated with transcriptional activation, more preferably said histone modification is a loss of transcriptional repressive chromatin marks, in particular H3K27me3, in said genomic sequence. Cancer according to the disclosure is preferably a breast cancer, preferably triple negative breast cancer.
In another aspect, the invention provides a combined preparation comprising a cancer drug and a compound that modulates the epigenetic status of the biomarker as described above, for use in cancer treatment, to reduce the development of resistance to said cancer treatment. Said combined preparation may comprise a cancer drug and a compound that modulates the epigenetic status of the biomarker as described above, for use in cancer treatment to reduce the development of resistance to said cancer treatment wherein said compound is administered before, after or concurrently with the therapeutic drug.
In the context of the invention, the term “treating” or “treatment”, as used herein, means reversing, alleviating, inhibiting the progress of, or preventing the disorder or condition to which such term applies, or reversing, alleviating, inhibiting the progress of, or preventing one or more symptoms of the disorder or condition to which such term applies.
“Treating cancer” includes, without limitation, reducing the number of cancer cells or the size of a tumor in the patient, reducing progression of a cancer to a more aggressive form (i.e. maintaining the cancer in a form that is susceptible to a therapeutic agent), reducing proliferation of cancer cells or reducing the speed of tumor growth, killing of cancer cells, reducing metastasis of cancer cells or reducing the likelihood of recurrence of a cancer in a subject. Treating a subject as used herein refers to any type of treatment that imparts a benefit to a subject afflicted with cancer or at risk of developing cancer or facing a cancer recurrence. Treatment includes improvement in the condition of the subject (e.g., in one or more symptoms), delay in the progression of the disease, delay in the onset of symptoms or slowing the progression of symptoms, etc.
As used herein, “drug” or “therapeutic agent” refers to a compound or agent that provides a desired biological or pharmacological effect when administered to a human or animal, particularly results in an intended therapeutic effect or response on the body to treat or prevent conditions or diseases. Therapeutic agents include any suitable biologically-active chemical compounds, biologically derived components such as for example small molecules, cells, proteins, peptides, antibodies, enzymes, polynucleotides, and radiochemical therapeutic agents, such as radioisotopes.
As used herein, a “therapeutic response” or “response to treatment with a drug” refers to a positive medical response characterized by objective parameters or criteria such as objective clinical signs of the disease, patient self-reported parameters and/or the increase of survival. The objective criteria for evaluating the response to drug-treatment will vary from one disease to another and can be determined easily by one skilled in the art by using clinical scores. A positive medical response to a drug can be readily verified in appropriate animal models of the disease which are well-known in the art and illustrated in the examples of the present application.
The term “determining resistance to a treatment with a drug”, as used herein, refers to an ability to assess whether the treatment of a patient with a drug will stop being effective in (e.g., stop providing a measurable benefit or positive medical response to) the subject after some time of administration of the treatment. In other terms, determining the resistance to a treatment refers to an ability to assess possible outcome of treatment of a cancer type. The resistance of a patient to a therapeutic agent may be determined by the lack of improvement in the disease state as measured by the absence of positive medical response, as compared to pre-treatment. In particular, such an ability to assess whether the treatment will stop being effective typically is exercised before treatment with the drug has begun in the subject.
However, it is also possible that such an ability to assess whether the treatment will stop being effective can be exercised after treatment has begun but before an indicator of ineffectiveness has been observed in the patient.
The term “predicting resistance to a treatment with a drug”, as used herein, refers to an ability to assess whether the treatment of a patient with a drug will stop being effective in (e.g., stop providing a measurable benefit or positive medical response to) the subject before treatment with the drug has begun in the subject.
The terms “subject” and “patient” are used interchangeably herein and refer to both human and non-human animals. As used herein, the term “patient” denotes a mammal, such as a rodent, a feline, a canine, and a primate. Preferably, a patient according to the invention is a human.
“a”, “an”, and “the” include plural referents, unless the context clearly indicates otherwise. As such, the term “a” (or “an”), “one or more” or “at least one” can be used interchangeably herein.
Biomarkers for Determining the Cancer Treatment Resistance
The inventors by using improved single-cell chromatin profiling identified epigenetic modifications in a specific genomic sequence characteristic of drug-resistant tumor cells.
These epigenetic modified genomic sequences can be used as biomarkers for determining drug-resistance in cancer patients, preferably before administration of cancer treatment and are listed in Tables 1 to 6. The genomic sequences listed in Tables 1 to 6 are identified by an identification number of the chromosomic region using the reference genome GRCh38 (Genome Reference Consortium Human Build 38 submitted in Dec. 17, 2013) (hg38) (GenBank assembly accession: GCA_000001405.15).
Certain embodiments of the present invention relate to a biomarker for determining resistance to treatment of a cancer type with a cancer drug which comprises one or more genomic sequence(s) selected from the list of Table 1 to 6, more preferably Tables 1 and 4. In a preferred embodiment, said genomic sequence comprises a histone modification, in particular a loss of H3K27me3, wherein said genomic sequence is selected from the list of Tables 1 to 6. In a preferred embodiment, the present invention relates to a predictive biomarker for predicting resistance to cancer treatment with a cancer drug prior administration of a treatment to a patient.
In another embodiment, the present invention relates to an isolated biomarker for determining resistance to or assessing possible outcome of treatment of a cancer type with a cancer drug in a patient which comprises purified nucleic acids having one or more genomic sequence(s) comprising a histone modification, wherein said one or more genomic sequence(s) is selected from the list of Tables 1 to 6.
The term “biomarker” refers to a distinctive biological or biologically derived indicator of a process, event or condition.
A “predictive biomarker” as used herein refers to a biomarker that can be used in advance of therapy to estimate the resistance of a patient suffering from a particular disease to a specific treatment of said disease. The biomarker for predicting the resistance of a patient to treatment with a drug prior administration to said treatment is herein referred to as pre-treatment predictive biomarker of drug-resistance.
By “epigenetic modification” or “epigenetic information” include with no limitations: histone modification, histone variant, DNA methylation, DNA modified bases and chromatin/DNA associated factors, preferably histone modification, histone variant and chromatin/DNA associated factors, more preferably histone modification.
By “characteristic epigenetic modification” or “specific epigenetic modification” it is meant an epigenetic modification that is present in a genomic sequence of tumor cells of a drug-resistant tumor and absent in the same genomic sequence of general population or from a selected population of subjects. The general population may comprise apparently healthy subjects, such as individuals who have not previously had any sign or symptoms indicating the presence of cancer. The term “healthy subjects” as used herein refers to a population of subjects who do not suffer from any known condition, and in particular, who are not affected with any cancer. In a preferred embodiment, the selected population may comprise subjects having an established cancer but who shows a clinically significant relief in a cancer type when treated with a cancer drug as described above.
In another embodiment, the epigenetic modification is present in a specific genomic sequence of tumor cells of a drug-resistant tumor and absent in the same genomic sequence of the majority of tumor cells of an untreated sensitive tumor from which the resistant-tumor is derived. As used herein, “the corresponding tumor cells” means tumor cells of the same tumor type. For example, the corresponding sensitive tumor cells may be the untreated tumor cells from which the resistant tumor cells are derived.
Preferably, epigenetic modifications are histone modifications. Histone modifications or histone post-translational modifications may be selected from the group comprising acetylation, amidation, deamidation, carboxylation, disulfide bond, formylation, glycosylation, hydroxylation, methylation, myristoylation, nitrosylation, phosphorylation, prenylation, ribosylation, sulphation, sumoylation, ubiquitination and derivatives thereof. Said histone modifications may be associated with transcriptional activation, such as for example histone H3 lysine 4 methylation (H3K4me3, H3K4me2 or H3K4me1) and histone acetylation. Alternatively, said histone modifications may be associated with transcriptional repression, such as for example histone H3 lysine 9 trimethylation (H3K9me3), H3K27me3 and H4K20me3.
In some embodiments, said specific epigenetic modification is loss of transcriptional repressive chromatin marks, in particular H3K27me3, in said genomic sequence.
Said specific epigenetic modification may be in a gene or in the proximity of said gene, i.e. less than 1 kb from the transcription start of said gene and said biomarker comprises one or more genomic sequence(s) selected from the list of Table 1 and 4. In a preferred embodiment, said biomarker comprising an epigenetic modification in a gene or in the proximity of said gene can be useful in determining resistance of a subject to a cancer treatment by determining the histone modifications of said biomarker and/or the expression level of said gene (e.g., mRNA or protein expression levels) in a patient sample.
In another embodiment, said specific epigenetic modification is not in a gene or in the proximity of said gene, i.e. more than 1 kb from the transcription start of said gene and said biomarker comprises one or more genomic sequence(s) selected from the list of Table 2 and 5. In a preferred embodiment, said biomarker comprising an epigenetic modification which is not in a gene or in the proximity of said gene can be useful in determining resistance of a subject to a cancer type by determining the epigenetic modifications of said biomarker, preferably by single-cell epigenetic profiling using a microfluidic system as described in examples of the present application.
The use of biomarker permits to determine the resistance to treatment of any cancer type, such as solid or liquid (or blood) cancer. In some embodiments, said biomarker determines the resistance to treatment of a cancer selected from the group consisting of: breast cancer, ovarian cancer, lung carcinoma, colorectal cancer, prostate cancer, pancreatic cancer and melanoma. In some preferred embodiments, said biomarker determines the resistance to treatment of breast cancer, preferably triple-negative breast cancer.
Targeted therapy includes the use of “targeted” drugs such as small molecule inhibitors or neutralizing monoclonal antibodies, that target proteins that are abnormally expressed in cancer cells and that are essential for their growth such as for example receptor and non-receptor tyrosine kinases, growth factors, hormone receptors, and others. Preferably, said targeted drugs are anti-hormonal drugs. Examples of anti-hormonal drugs include with no limitations: Tamoxifen, targeting the estrogen receptor; The cancer agent may be a drug for chemotherapy or targeted therapy. Chemotherapy includes the use of cytotoxic anti-neoplastic agents, such as alkylating agents, anti-metabolites, anti-microtubule agents, Topoisomerase inhibitors, cytotoxic antibiotics and others. Examples of chemotherapeutic drugs include with no limitations: Capecitabine, 5-FU, docetaxel, SN-38, CPT11, cisplatin, carboplatin, etc.
In some preferred embodiment, said cancer drug is a chemotherapy drug such as Capecitabine and said biomarker comprises one or more genomic sequence(s) selected from the list of Tables 1 to 3. Epigenetic modifications in the above listed genomic sequences, in particular loss of H3K27me3 in said genes, are found in Capecitabine resistant tumors, in particular Triple-negative breast cancer tumors.
In a preferred embodiment, said epigenetic modification is in a gene or in the proximity of said gene and the biomarker which determines the resistance to chemotherapy drug such as Capecitabine comprises one or more genomic sequence(s) selected from the list of Table 1.
In another preferred embodiments, said cancer drug is an anti-hormonal drug such as Tamoxifen, targeting the oestrogen receptor and said biomarker comprises one or more genomic sequence(s) selected from the list of Tables 4 to 6. Epigenetic modifications in the above-listed genomic sequences, in particular loss of H3K27me3 in said genes, are found in Tamoxifen resistant tumors, in particular luminal ER+ breast cancer tumors.
In a preferred embodiment, said epigenetic modification is in a gene or in the proximity of said gene and the biomarker which determines the resistance to anti-hormonal drug such as Tamoxifen comprises one or more genomic sequences selected from the list of Table 4.
In a preferred embodiment, the present disclosure provides a biomarker for predicting resistance to treatment of a cancer type prior treatment administration, wherein said biomarker comprises one or more epigenetic modified genomic sequence(s) as described above.
Method for Determining the Cancer Treatment Resistance
In another aspect, the present invention relates to a method for determining the resistance to treatment of a cancer type with a cancer drug to a patient comprising detecting a histone modification in said patient tumor sample of at least one biomarker as described above which comprises one or more genomic sequence selected from the list of Table 1 to 6.
In some preferred embodiments, said specific histone modification as disclosed above is a loss of transcriptional repressive chromatin marks, in particular H3K27me3.
According to the invention, the specific epigenetic modifications present in tumor cells can be identified by any methods well-known in the art including but not limiting to ChIP-qPCR, ChIPseq, ChIP on chip, EMSA, ATACseq, FISH, immunofluorescence, immuno-histochemistry CITEseq, Chem-Seq, DNAase-Seq, Hi-C, DAM-ID, TIRF microscopy (https://doi.org/10.1038/s4158), Split-seq. Preferably, specific epigenetic modifications, such as histone modification can be identified by single-cell epigenetic profiling using a microfluidic system as described in the examples of the present application.
In some embodiments, said at least one histone modification is present in at least 0.01%, preferably 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 1%, more preferably between 0.5% to 20% of tumor cells of the untreated sensitive tumor from which said drug-resistant tumor is derived.
In another embodiment, the method for determining the resistance to treatment of a cancer type with a cancer drug to a patient comprises determining in said patient tumor sample the expression level of at least one gene encoding by a biomarker as described above comprising one or more genomic sequence(s) selected from the list of Table 1 and 4.
Typically, the expression level of gene encoded by said biomarker in a patient sample is deemed to be higher or lower than the predetermined value obtained from the general population or from healthy subjects if the ratio of the expression level of said gene encoded by said biomarker in said patient to that of said predetermined value is higher or lower than 1.2, preferably 1.5, even more preferably 2, even more preferably 5, 10 or 20.
As used herein, the term “predetermined value of a biomarker” refers to the amount of the biomarker in biological samples obtained from the general population or from a selected population of subjects. For example, the general population may comprise apparently healthy subjects, such as individuals who have not previously had any sign or symptoms indicating the presence of cancer. The term “healthy subjects” as used herein refers to a population of subjects who do not suffer from any known condition, and in particular, who are not affected with any cancer. In another example, the predetermined value may be of the amount of biomarker obtained from selected population of subjects having an established cancer but who shows a clinically significant relief in a cancer type when treated with a cancer drug as described above. The predetermined value can be a threshold value, or a range. The predetermined value can be established based upon comparative measurements between apparently healthy subjects and subjects with established cancer.
The expression level of said gene encoded by a biomarker may be determined by any suitable methods known by skilled persons. Usually, these methods comprise measuring the quantity of mRNA or protein. Methods for determining the quantity of mRNA are well known in the art. For example the mRNA contained in the sample is first extracted according to standard methods, for example using lytic enzymes or chemical solutions or extracted by nucleic-acid-binding resins following the manufacturer's instructions. The extracted mRNA is then detected by hybridization (e.g., Northern blot analysis) and/or amplification (e.g., RT-PCR). Quantitative or semi-quantitative RT-PCR is preferred. In a preferred embodiment, the mRNA expression level is measured by RNA seq method, more preferably by single-cell RNA-seq. RNA seq can be used to analyse the cellular transcriptome. RNAseq, preferably single cell RNA seq can be performed for example in plate, micro or nano-wells, droplet-based microfluidics, microfluidics, tubes. The approach aims at deciphering tissue, sample, cell heterogeneity by interrogating genetic expression of whole transcriptome, or subset of gene, and classifying cells based on similar/closest transcriptomic expression pattern. The approach can also include as well interrogation of protein expression using nucleic acid tagged antibody against cell surface protein. One example of single cell RNA-seq method is illustrated in the
The level of the protein may be determined by any suitable methods known by skilled persons. Usually, these methods comprise contacting a cell sample, preferably a cell lysate, with a binding partner capable of selectively interacting with the protein present in the sample. The binding partner is generally a polyclonal or monoclonal antibodies, preferably monoclonal. The quantity of the protein may be measured, for example, by semi-quantitative Western blots, enzyme-labelled and mediated immunoassays, such as ELISAs, biotin/avidin type assays, radioimmunoassay, immunoelectrophoresis or immunoprecipitation or by protein or antibody arrays. The reactions generally include revealing labels such as fluorescent, chemiluminescent, radioactive, enzymatic labels or dye molecules, or other methods for detecting the formation of a complex between the antigen and the antibody or antibodies reacted therewith.
The term “patient sample” means any biological sample derived from a patient. Examples of such samples include fluids, tissues, cell samples, organs, biopsies, etc. Preferred biological samples are tumor sample.
The tumor sample may be from a patient tumor biopsy or a patient-derived xenograft (PDX) model of the cancer as disclosed in the examples of the present application. Drug-resistant tumor cells may be isolated directly from sample of patient drug-resistant tumor or generated from sample of patient untreated sensitive tumor or patient-derived xenograft by several rounds of drug treatment as disclosed in the examples of the present application. Various PDX models of cancers are available in the art. PDX models useful to perform the method of the present invention include with no limitations: luminal ER+ breast cancer (HBCx-22; Cottu et al., Breast Cancer Res. Treat., 2012, 133, 595-606) and derived Tamoxifen resistant cells (HBCx-22-TamR; Cottu et al., Clin. Cancer Res., 2014, 20, 4314-4325); Triple-negative breast cancer (HBCx-95; Cottu et al., Clin. Cancer Res., 2014, 20, 4314-4325; Marangoni et al., Clin. Cancer Res., 2018, 24, 2605-2615).
In some embodiments, said patient-derived tumor cells are from breast cancer, ovarian cancer, lung carcinoma, colorectal cancer, prostate cancer, pancreatic cancer and melanoma. Breast cancer include estrogen receptor positive (ER+), progesterone positive (PR+), HER2 positive (HER2+) and triple-negative (ER−, PR−, HER2−) breast cancer. In some preferred embodiments, said patient-derived tumor cells are from breast cancer, preferably triple-negative breast cancer.
“Triple-negative breast cancer” refers to any breast cancer that does not overexpress the genes for estrogen receptor (ER), progesterone receptor (PR) and HER2/Neu. This subtype of breast cancer is clinically characterized as more aggressive and less responsive to standard treatment and associated with poorer overall patient survival.
The presence of the biomarker(s) in the patient sample indicates that the patient is likely to be resistant to said cancer treatment, whereas the absence of the biomarker(s) indicates that the patient is likely to be responsive to said cancer treatment.
In some advantageous embodiments, the method comprises detecting at least one histone modification in at least one genomic sequence selected from the list of Table 1 to 6, and determining therefrom whether or not said patient is likely to be resistant to said cancer-drug treatment.
In another aspect, the present invention relates to a method for predicting the resistance to treatment of a cancer type with a cancer drug prior administration of said treatment to a patient, comprising:
In some advantageous embodiments, the method further comprises a step of classification of the patients into resistant and sensitive group based on the presence or absence of the epigenetic biomarker(s) according to the invention.
In some advantageous embodiments, when the patient is found likely to be resistant to said cancer-drug treatment, the method further includes administering a therapeutically effective amount of a compound that modulates the epigenetic status of the genomic region of interest comprising the epigenetic modification, to reduce the development of resistance to cancer treatment in said patient.
Therefore, the use of the prediction method of the invention increases the efficiency of cancer treatment by reducing the development of resistance to cancer treatment.
Combined Therapy to Reduce the Development of Drug-Resistance
In connection with the above method of prediction of resistance to cancer treatment, the present invention relates to a combined preparation comprising a cancer drug and a compound that modulates the epigenetic status of the biomarker as described above for use in cancer treatment to reduce the development of resistance to said cancer treatment. In some preferred embodiments, said combined preparation is used in a cancer patient previously classified as resistant to treatment with said cancer drug using the method for determining resistance to cancer treatment according to the invention.
The present invention relates also to a method of treating a cancer patient comprising administering to said patient a therapeutically effective amount of a cancer drug and a therapeutically effective amount of a compound that modulates the epigenetic status of the biomarker. In some preferred embodiments, said combined preparation is administered to a cancer patient previously classified as resistant to treatment with said cancer drug using the method for determining resistance to cancer treatment according to the invention.
As used herein, a “therapeutically effective amount” or an “effective amount” means the amount of a composition that, when administered to a subject for treating a state, disorder or condition is sufficient to effect a treatment. The therapeutically effective amount will vary depending on the compound, formulation or composition, the disease and its severity and the age, weight, physical condition and responsiveness of the subject to be treated.
Such compounds that modulate the epigenetic status of a genomic region of interest include histone deacetylase (HDAC) inhibitor, DNA methyltransferase (DNMT) inhibitors, and Histone Methyl Transferase (HMT) inhibitors. In a particular embodiment, said compound is a DNA demethylase inhibitor.
As used herein, a “demethylase inhibitor” is any agent capable of partially or fully inhibiting one or more of the biological activities of a histone demethylase protein including, without limitation, a polypeptide, a polynucleotide, or a small molecule. Histone demethylases are a family of enzymes that catalyze the removal of methyl groups from lysine and arginine residues on histone tails. A demethylase inhibitor, can be KDM1 inhibitors or JmjC KDM inhibitors.
In a more particular embodiment, Histone Lysine Demethylase (KDM) inhibitor which can be used includes but is not limited to tranylcypromine ((trans-2-phenylcyclopropyl-1-amine, trans-2-P CPA)) and analogs thereof, e.g., with substitutions at the benzene ring, e.g., tranylcypromine 7, trans-2-PCPA analogue 28 (trans-2-pentafluorophenylcyclopropylamine, 2-PFPA), and trans-2-PCPA analogues carrying 4-bromo, 4-methoxy, and 4-trifluoromethoxy substitutions at the benzene ring (see, e.g., Gooden et al, Bioorg Med Chem Lett. 2008; 18(10):3047-51; Binda et al, J Am Chem Soc. 2010; 132(19):6827-33; Mimasu et al, Biochemistry. 2010; 49(30):6494-503; Benelkebir et al, Bioorg Med Chem. 2011; 19(12):3709-16; and Mimasu et al. Biochem Biophys Res Commun. 2008; 366(1): 15-22) or other inhibitors, e.g., 2,4-pyridinedicarboxylic acid (2,4-PDCA) (see, e.g., Kristensen et al, FEBS J. 2012 June; 279(1 1): 1905-14), and inhibitors of jumonji C (jmjC)-containing KDMs, e.g., 5-Carboxy-8-hydroxyquinoline (IOX1) and n-octyl ester thereof, as described in Schiller et al, Chem Med Chem. 2014 March; 9(3):566-71, or other inhibitors, as described in Spannhoff et al, Chem Med Chem. 2009; 4(10):1568-82; Varier and Timmers, Biochim Biophys Acta. 2011; 1815(1):75-89; Luo et al, J Am Chem Soc. Jun. 22, 2011; 133(24): 9451-9456; and Rotili et al, J Med Chem. 2014 Jan. 9; 57(1):42-55.
A number of HDAC inhibitors are known in the art, including but not limited to: Sodium Butyrate, Trichostatin A, hydroxamic acids, cyclic tetrapeptides, trapoxin B, depsipeptides, benzamides, electrophilic ketones, aliphatic acid compounds, pyroxamide, valproic acid, phenylbutyrate, valproic acid, hydroxamic acids, romidepsin; CI-994 (N-acetyldinaline, also tacedinaline); vorinostat (SAHA), belinostat (PXD101), LAQ824, panobinostat (LBH589), Entinostat (SNDX-275; formerly MS-275), EVP-0334, SRT501, CUDC-101, JNJ-26481585, PC124781, Givinostat (ITF2357), and mocetinostat (MGCD0103).
A number of DNMT inhibitors are known in the art, including but not limited to azacytidine, decitabine, Zebularine (1-(β-D-ribofuranosyl)-1,2-dihydropyrimidin-2-one), procainamide, procaine, (−)-epigallocatechin-3-gallate, MG98, hydralazine, RG108, and Chlorogenic acid. See also Gros et al, Biochimie. 2012 November; 94(11):2280-96.
A number of EZH2/HMT inhibitors are known in the art, including but not limited to: EPZ005687; E7438; Ell (Qi et al, 2012, supra); EPZ-6438; GSK343; BLX-01294, U C0638, BRD4770, EPZ004777, AZ505 and PDB 4e47, and those described in Garapaty-Rao et al, Chem. Biol. 20(11): 1329-1339 (2013); Ceccaldi et al, ACS Chem Biol. 2013 Mar. 15; 8(3):543-8; US 20130303555; and WO2012/005805; see, e.g., Wagner and Jung, Nature Biotechnology 30:622-623(2012), and Yao et al., J Am Chem Soc. 2011 Oct. 26; 133(42): 16746-9. In some embodiments, inhibitors that act on the G9A H3K9 methyltransferase, are used, e.g., BIX-01294 or BRD4770.
The cancer drug is any drug for chemotherapy or targeted therapy as disclosed above. In some preferred embodiments, said cancer drug is Capecitabine, Tamoxifen, or others.
The cancer is any cancer type as disclosed above. In some embodiments said cancer is selected from the group consisting of: breast cancer, ovarian cancer, lung carcinoma, colorectal cancer, prostate cancer, pancreatic cancer and melanoma. In some preferred embodiments, said cancer is breast cancer, preferably triple-negative breast cancer.
The compound that modulates the epigenetic status of the genomic region of interest and the cancer drug may be used simultaneously, separately or sequentially. The compound that modulates the epigenetic status of the genomic region of interest may be administered before, after, or concurrently with the therapeutic drug. Preferably, to reduce the resistance to cancer treatment, the compound that modulates the epigenetic status of the genomic region of interest may be administered prior to the cancer drug by at least 6 hours, 12 hours, 1 days, 2 days, 3 days, 5 days, 1 week.
The compounds or cancer drugs described herein may be administered by any means known to those skilled in the art, including, without limitation, intravenously, orally, intra-tumoral, intra-lesional, intradermal, topical, intraperitoneal, intramuscular, parenteral, subcutaneous and topical administration. Thus the compositions may be formulated as an injectable, topical, ingestible, or suppository formulation. Administration of the compounds or therapeutic agents to a subject in accordance with the present invention may exhibit beneficial effects in a dose-dependent manner. Thus, within broad limits, administration of larger quantities of the compositions is expected to achieve increased beneficial biological effects than administration of a smaller amount. Moreover, efficacy is also contemplated at dosages below the level at which toxicity is seen.
It will be appreciated that the specific dosage of compounds or cancer drugs administered in any given case will be adjusted in accordance with the composition or compositions being administered, the volume of the composition that can be effectively delivered to the site of administration, the disease to be treated or inhibited, the condition of the subject, and other relevant medical factors that may modify the activity of the compositions or the response of the subject, as is well known by those skilled in the art.
For example, the specific dose of compounds or cancer drugs for a particular subject depends on age, body weight, general state of health, diet, the timing and mode of administration, the rate of excretion, medicaments used in combination and the severity of the particular disorder to which the therapy is applied. Dosages for a given patient can be determined using conventional considerations, e.g., by customary comparison of the differential activities of the compositions described herein and of a known agent, such as by means of an appropriate conventional pharmacological protocol. The compositions can be given in a single dose schedule, or in a multiple dose schedule.
Suitable dosage ranges for a compound that modulates the epigenetic status and/or cancer drug may be of the order of several hundred micrograms of the agent with a range from about 0.001 to 10 mg/kg/day, preferably in the range from about 0.01 to 1 mg/kg/day.
The invention will now be exemplified with the following examples, which are not limitative, with reference to the attached drawings in which:
Material and Methods
1. Cell Lines
Jurkat cells (ATCC, T18-125), an immortalized line of human T lymphocytes and Ramos cells (ATCC, CRL-1596), an immortalized line of human B lymphocytes, were grown in RPMI medium (ThermoFisher Scientific, #61870010) supplemented with 10% heat inactivated bovine serum (ThermoFisher Scientific, #16140071) and 1% Pen/Strep (ThermoFisher Scientific, #15140122). Mouse M300.19 cells (a gift from B. Moser), an immortalized line of mouse pre-B lymphocytes, were grown in RPMI 1640 medium (ThermoFisher Scientific, #61870010) supplemented with 10% fetal bovine serum (Fisher, #SH30070.03), 1% Pen/Strep (ThermoFisher Scientific, #15140122), 1% L-glutamine (ThermoFisher Scientific, #25030081) and 5×10−5 M P-mercapto-ethanol (ThermoFisher Scientific, #21985023).
2. Patient-Derived Xenografts (PDX)
Female Swiss nude mice were purchased from Charles River Laboratories and maintained under specific pathogen-free conditions. Their care and housing were in accordance with institutional guidelines and the rules of the French Ethics Committee (project authorization no. 02163.02). A PDX model of luminal breast cancer (HBCx-22) was previously established at Institut Curie from untreated early-stage luminal breast cancer with informed consent from the patient. Acquisition of a resistant phenotype for a derivative of HBCx-22, HBCx-22-TamR, was previously established and maintained. A PDX from a residual triple negative breast cancer post neo-adjuvant chemotherapy (HBCx-95) was previously established at Institut Curie with informed consent from the patient. HBCx-95 xenografts (n=6) were treated with Capecitabine (Xeloda, Roche Laboratories) orally at a dose of 540 mg/kg/day, 5 days a week for 6 weeks. Relative tumor volumes (RTV, mm3) were calculated. Mice with recurrent tumors were treated for a second round of Capecitabine when PDX reached a volume of over 200 mm3 (mice #35, #40 & #33). Mouse #40 did not respond to Capecitabine and PDX specimen was extracted at 1100 mm3 and tagged as HBCx-95-CapaR.
Prior to single-cell ChIP-seq, single-cell RNA-seq and bulk ChIP-seq, PDX were digested at 37° C. for 2 h with a cocktail of Collagenase I (Roche, #11088793001) and Hyaluronidase (Sigma, #H3506). Cells were further individualized at 37° C. using a cocktail of 0.25% trypsin/Versen (ThermoFisher Scientific, #15040-033), Dispase II (Sigma, #D4693) and Dnase I (Roche, #11284932001). Red Blood Cell lysis buffer (ThermoFisher Scientific, #00-4333-57) was then added to degrade red blood cells. To increase the viability of the cell suspension, dead cells were removed using the Dead Cell Removal kit (Miltenyi Biotec). Cells were re-suspended in PBS/0.04% BSA (ThermoFisher Scientific, #AM2616).
3. Single-Cell ChIP-Seq
3.1 Microfluidic Chips
Four microfluidic chips were used: i) to compartmentalize single cells with lysis reagents and MNase in droplets; ii) to produce hydrogel beads; iii) to compartmentalize single hydrogel beads in droplets, and iv) for one-to-one fusion of droplets containing digested nucleosomes (from single lysed cells) with droplets containing single hydrogel beads (
3.2 Microfluidic Operations
Droplet formation, fusion and fluorescence analysis was performed on a dedicated droplet microfluidic station, similar to Mazutis L. et al. 2013, Nat. Protoc. 8:873-891. The continuous oil phase for all droplet microfluidics experiments was Novec HFE7500 fluorinated oil (3M) containing 2% w/w 008-FluoroSurfactant (RAN Biotechnologies).
3.3 Cell Compartmentalization and Chromatin Digestion.
Cells were centrifuged (300 g, 5 min at 4° C.), labeled by 20 min incubation with 1 μM Calcein AM (ThermoFisher Scientific, #C3099). Then, cells are resuspended in cell suspension buffer, comprising DMEM/F12 (ThermoFisher Scientific) supplemented with 30% Percoll (Sigma, #P1644), 0.1% Pluronic F68 (ThermoFisher Scientific, #24040032), 25 mM Hepes pH 7.4 (ThermoFisher Scientific, #15630080) and 50 mM NaCl.
Cells were resuspended to give an average number of cells per droplets X of 0.1, resulting in 90.48% of empty droplets, 9.05% of droplets containing one cell and only 0.46% containing two or more cells due to Poisson distribution of the cells in droplets (Clausell-Tormos, J. et al. 2008. Chem Biol 15, 427-437). Overall, the inventors estimated that among non-empty droplets, 95.16% contained one cell and 4.84% contained two cells or more close to the expected values (94.92 and 5.08%, respectively).
The cells were co-flowed in a microfluidic chip (
The droplets were collected in a collection tube (1.5 ml Eppendorf tube filled with HFE-7500 fluorinated oil) and then incubated at 37° C. for 20 min.
3.4 Production of Hydrogel Beads Carrying Barcoded DNA Adaptors
3.4.1 DNA Barcoding Strategy
Hydrogel beads carrying barcoded DNA adaptors were produced by split-and-mix synthesis using a method similar to that previously described in Zilionis, R. et al. 2017. Nat Protoc 12:44-73; Klein, A. M. et al. 2015. Cell; 161:1187-1201. Briefly, polyethylene diacrylate (PEG-DA) hydrogel beads containing streptavidin acrylamide were produced and barcoded primers were added to the beads by split-and-pool synthesis using ligation (
The 9 μl droplets were produced at 4.5 kHz frequency and were exposed at 200 mW/cm2 with a 365 nm UV light source (OmniCure ac475-365) to trigger gel bead polymerization. Recovered gel beads were washed 10 times with washing buffer (100 mM Tris pH 7.4, 0.1% v/v Tween 20). Twelve million PEG-DA beads were incubated in 500 μl final volume for 1 h at room temperature with 50 μM of a photo-cleavable biotinylated dsDNA oligonucleotide (see SEQ ID NO: 1 and 2) and then distributed into a 96-well plate, each well containing 5 μl at 5 μM of a double-stranded DNA with a specific first index (index 1), for split-and-pool synthesis by ligation using T7 DNA ligase (NEB, #M0318) according to the manufacturer's instructions. At each round of split-and-pool, the hydrogel beads were pooled and washed (
3.4.2 Compartmentalization of Hydrogel Beads
The barcoded hydrogel beads were labeled by 30 min incubation with 10 μM Cy5-PEG3 biotin (Bioscience Interchim, #FP-1M1220) and washed with washing buffer (100 mM Tris pH 7.4, 0.1% v/v Tween 20), then suspended in bead mix (62.5 mM EGTA, 2 mM dNTPs, 1 mM ATP, 0.5 μM Sulforhodamine B). Barcoded hydrogel beads were co-flowed using the microfluidic device, with ligation mix (2× ligation buffer, 2 mM ATP, 1 μM Sulforhodamine B, 100 mM EGTA, 0.38 U/μl Fast-link DNA ligase [Lucigen, #LK0750H]) and EndRepair mix (4× ligation buffer, 4 mM dNTPs, 1 μM Sulforhodamine B, 0.08 U/μl Fast-link DNA ligase [Lucigen, #LK0750H], 0.15× ENDit repair mix [Lucigen, #ER0720]). The re-injection of close-packed barcoded hydrogel beads resulted in 65±5% of the droplets containing a single bead. The flow rates (150 μl/hr for the beads, 75 μl/hr for both ligation and EndRepair buffers, 150 μl/hr for the continuous oil phase) were calibrated to produce 100 pl droplets.
3.5 Fusion of Beads and Cell Droplets
Droplets containing fragmented chromatin and droplets containing barcoded hydrogel beads were re-injected into a microfluidic device with two aqueous inlets and one oil inlet for droplet fusion (
3.6 Nucleosomes Barcoding in Droplets
Fused droplets were collected and exposed for 90 seconds at 200 mW/cm2 with a 365 nm UV light source (OmniCure ac475-365). The ligation was performed at 16° C. overnight. The emulsion was then broken by addition of 1 volume of 80/20 v/v HFE-7500/1H,1H,2H,2H perfluoro-1-octanol (Sigma, #370533). The aqueous phase containing barcoded-nucleosomes was diluted by addition of 10 volumes of lysis dilution buffer (50 mM Tris-HCl pH 7.4, 1% Triton Tx-100, 0.1% DOC, 37.5 mM EDTA, 37.5 mM EGTA, 262.5 mM NaCl and 1.25 mM CaCl2)) and centrifuged 10 min at 10,000 g at 4° C. The soluble aqueous phase was used for the chromatin immunoprecipitation.
3.7 Immunoprecipitation of Barcoded-Nucleosomes
Protein-A magnetic particles (ThermoFisher Scientific, 10001D) were washed in blocking buffer comprising phosphate buffered saline (PBS) supplemented with 0.5% Tween 20, 0.5% BSA fraction V and incubated 4 hours at 4° C. in 1 ml blocking buffer with 2 μg of antibody (anti-H3K4me3 [Millipore, #07-473] and anti-H3K27me3 [Cell Signaling Technology, #9733]). After incubation, the particles were suspended with the barcoded-nucleosomes and incubated at 4° C. overnight. Magnetic particles were washed as described in Rotem, A. et al. 2015. Nat Biotechnol 33, 1165-1172 and immediately processed to prepare the sequencing library.
3.8 Sequencing Library Preparation & Sequencing
Concatemers of barcodes were digested by Pac restriction enzyme (NEB, #R0547), following the manufacturer's instructions. Immunoprecipitated chromatin was then treated with RNAse A (ThermoFisher Scientific, #EN0531) and with Proteinase K (ThermoFisher Scientific, #E00491). DNA was eluted from the magnetic particles with 1 volume of elution buffer (1% SDS, 10 mM Tris-HCl pH 8, 600 mM NaCl and 10 mM EDTA). Eluted DNA was purified with 1× AMPure XP beads (Beckman, #A63881) and eluted with RNAse/Dnase free water. Barcoded-nucleosomes were amplified by in vitro transcription using the T7 MegaScript kit (ThermoFisher Scientific, #AM1334). The resulting amplified RNA was purified using 1× RNAClean XP beads (Beckman, #A66514) and reverse transcribed using SEQ ID NO: 5 (
Single-cell ChIP-seq libraries were sequenced on an Illumina NextSeq 500 MidOutput 150 cycles. Cycles were distributed as follows: 50 bp (Read #1) were assigned for the genomic sequence and 100 bp (Read #2) were assigned to the barcode (
4. Single-Cell ChIP-Seq Data Analysis
Sequencing data were analyzed with Python (v2.7.12) and R (v3.3.3) using the reference genome hg38 (GenBank assembly accession: GCA_000001405.15 (Genome Reference Consortium Human Build 38 submitted in Dec. 17, 2013).
4.1 De-Multiplexing Cellular Barcodes
Barcodes were extracted from Reads #2 by first searching for the constant 4 bp linkers found between the 20-mer indices of the barcode allowing up to 1 mismatch in each linker (
4.2 Alignment, Filtering & Normalization
Reads #1 were aligned to mouse mm10 and human hg38 reference genomes using bowtie (v1.2.2) by keeping only reads having no more than one reportable alignments and 2 mismatches.
Raw reads are distributed according to a bimodal distribution, the lower peak most probably corresponding to droplets with barcoded beads but without cells (Rotem, A. et al. Nat Biotechnol; 2015, 33:1165-1172), and the right peak corresponding to droplets with cell with bead (
For each cell, reads were binned in non-overlapping 50 kb for H3K27me3, known to accumulate over broad genomic regions, and 5 kb genomic bins for H3K4me3, known to accumulate in narrow peaks around transcription start sites, spanning the genome to generate a n×m coverage matrix with n barcodes and m genomic bins. The inventors combined coverage matrices for each of the four analyses from the following samples: (i) Ramos and Jurkat, (ii) mouse cells from HBCx95 and HBCx-95-CapaR, (iii) human cells from HBCx95 and HBCx95-CapaR, and (iv) human cells from HBCx-22 and HBCx-22-TamR.
The inventors first removed cells with a total number of uniquely mapped reads within the upper percentile, considered as outliers, and filtered out genomic regions not represented in at least 1% of all cells. By PCA analysis, the inventors could group cells independently of coverage only if cells had at least 1,600 unique reads per cell. For all subsequent analyses, the inventors excluded cells with lower coverage. Coverage matrices were then normalized by dividing counts by the total number of reads per cell and multiplying by the average number of reads across all cells.
4.3 Unsupervised Clustering of Single-Cell ChIP-Seq Profiles
Normalized matrices were reduced by principal component analysis (n=50 first components selected for further analysis). To improve the stability of the clustering approaches, the inventors further limited the analysis to cells displaying a Pearson's pairwise correlation score above a threshold t with at least 1% of cells. Threshold t was defined as the upper percentile of Pearson's pairwise correlation scores for a randomized dataset.
The inventors used consensus clustering, Bioconductor ConsensusClusterPlus package (Wilkerson, M. D. & Hayes, D. N. 2010. Bioinformatics 26, 1572-1573), to examine the stability of the clusters and compute item consensus score for each cell. The inventors established consensus partitions of the data set in k clusters (for k=2, 3, . . . ), on the basis of 1,000 resampling iterations (80% of cells) of hierarchical clustering, with Pearson's dissimilarity as the distance metric and Ward's method for linkage analysis. The optimal number of clusters (k) was chosen to maximize intra-cluster correlation score. Clustering results were visualized with t-SNE plots (Van der Maaten, L. & Hinton, G. 2008. J Mach Learn Res 9, 2579-2605). To visualize chromatin profiles of subpopulations, the inventors aggregated reads of single-cells within each cluster and created enrichment profiles using the R package Sushi (Phanstiel, D. H. et al. 2014. Bioinformatics. 30, 2808-2810).
4.4 Differential Analysis of Single-Cell ChIP-Seg Profiles
To identify differentially enriched regions across single-cells for a given cluster, the inventors performed a non-parametric Wilcoxon rank sum test comparing normalized counts from individual cells from one cluster versus all other cells. The inventors testes for the null hypothesis that the distribution of normalized counts from the two compared groups have the same median, with a confidence interval 0.95. The inventors limited analysis to the windows selected for unsupervised analysis above.
P-values were corrected for multiple testing using Benjamini-Hochberg procedure (Benjamini, Y. & Hochberg, Y. 1995. J R Stat Soc 57, 289-300). Genomic regions were considered as ‘enriched’ or ‘depleted’ for H3K27me3 or H3K4me3 if adjusted p-values, ‘q-values’, were lower than 0.01 and the absolute log 2 fold change greater than 1.
4.5 scRNA-Seq Comparison.
For H3K27me3 scChIP-seq analysis, the inventors used peak annotation from bulk ChIP-seq datasets to further annotate 50 kbp windows and corresponding genes: for each window, the inventors kept for subsequent analyses (gene annotation and scRNA-seq comparison), genes with a transcription start site (TSS) overlapped by a peak in any condition, using bedtools (v2.17)50 and the reference annotation of the human transcriptome Gencode_hg38_v26, limited to protein coding, antisense and lncRNA genes.
5. Bulk ChIP-seq
ChIP experiments were performed as described previously in Vallot, C. et al. 2015. Cell Stem Cell 16, 533-546, on 106 cells from cell suspensions obtained above from HBCx-22, HBCx-22-TamR, HBCx-95 and HBCx-95-CapaR using anti-H3K27me3 antibody (Cell Signaling Technology, #9733). 2 ng of immune-precipitated and input DNA were used to prepare sequencing libraries using the Ovation Ultralow Library System V2 (Nugene) according to the manufacturer's instructions. Bulk ChIP-seq libraries were sequenced on an Illumina HiSeq 2500 in Rapid run mode SE50.
6. Bulk ChIP-Seq Data Analysis
Reads were aligned to mouse mm10 and human hg38 reference genomes using bowtie (v1.2.2) and the tool bamcmp was used to separate human from mouse sequences. Subsequent analysis was performed as explained previously. Only uniquely mapping reads were kept for the analysis, in addition PCR duplicates were removed using Picard Tools (https://broadinstitute.github.io/picard/). Data were binned in 5 kb (H3K4me3) or 50 kb (H3K27me3) consecutive genomic windows. For each window, log 2 RPM were computed, as the logged number of reads per millions of mapped reads.
7. Single-Cell RNA-Seq
Approximately 3,000 cells from each cell suspensions, HBCx-22, HBCx-22-TamR, HBCx-95 and HBCx-95-CapaR, were loaded on a Chromium Single Cell Controller Instrument (10× Genomics) according to the manufacturer's instructions. Samples and libraries were prepared according to the manufacterer's instructions. Libraries were sequenced on Illumina an HiSeq 2500 in Rapid run mode, using paired-end 26 bp-98 bp sequencing.
8. Single-Cell RNA-Seq Data Analysis
Single-cell sequencing files were processed using the Cell Ranger Single Cell Software Suite (v1.3.1) to perform quality control, sample de-multiplexing, barcode processing, and single-cell 3′ gene counting (http://software.10xgenomics.com/single-cell/overview/welcome) using the UCSC mouse (mm10) and human (hg19) transcriptome and genome with default parameters. 2,728 cells with an average coverage of 30,166 reads/cell (1,564 human and 1,191 mouse cells) for HBCx-22, 1,746 cells with an average coverage of 41,166 reads/cell (753 human and 1,013 mouse cells) for HBCx-22-TamR, 1,184 cells with an average coverage of 160,583 reads/cell (545 human and 647 mouse cells) for HBCx-95, 2,087 cells with an average coverage of 38,345 reads/cell (861 human and 1,242 mouse cells) for HBCx-95-CapaR were analyzed. Further analysis was performed in R (v3.3.3) using custom R scripts. Any cell with more than 10% of mitochondrial UMI counts was filtered out. The inventors only kept cells with a total count of UMI below 100,000 and of detected genes below 6,000 and over 1,000. The inventors then only kept genes with at least 1 transcript in at least 2 cells. Using the R package scater, scRNAseq count matrices were normalized for coverage and transformed by RLE, ‘Relative Log Expression’ method (McCarthy, D. J et al. Bioinformatics, 2017, 33:1179-1186). Using annotations from the R package ccRemover (Barron, M. et al. Sci Rep, 2016, 6:33892), the inventors removed genes related to cell cycle from subsequent clustering analyses to group cells according to cell identity and not cell-cycle related phenomena. Barnes-hut approximation to t-SNE was then performed on the n=50 first principal components (PCA) to visualize cells in a two-dimensional space. Clusters were identified using Consensus Clustering as for scChIP-seq analyses above. The inventors identified genes that were differentially expressed between clusters using edgeR GLM statistical models (Robinson, M. D. et al. Bioinformatics, 2010, 26:139-140). For stromal mouse cells, clusters were identified according to the differential expression of hallmark genes.
9. Copy Number Profiles of Bulk Tumor Cells
The R package HMMcopy was used to correct for copy number variation in non-treated versus resistant xenograft models. Reads from bulk input ChIP-seq samples were binned in 0.5 Mb non-overlapping regions spanning the genome. Regions with a deviation to the mean greater than n=2 standard deviations were removed for analysis.
Results
Characterization of Epigenetic and Transcriptional Makers
The inventors studied the heterogeneity of chromatin profiles among tumor cells from the same pair of triple-negative breast tumor samples (n=4,331 cells from HBCx_95 and HBCx_CapaR, with average coverage of 5,160 reads per cell).
The inventors removed from the analysis loci affected by copy-number variations, as identified from bulk DNA profiles, to focus on chromatin alterations (
While the chromatin profiles of sensitive cells were largely homogeneous, distinct chromatin states within the resistant population were apparent (
However, consensus clustering also showed that 3% of cells of from the untreated tumor (n=13 out of 484) robustly classify with resistant cells (resistant-like cells) (
Loci specifically devoid of H3K27me3 in cells from Chrom_c2 were enriched in genes targets of the Polycomb complex (
Interestingly, the inventors identified a genomic region including IGF2BP3, a gene known to promote resistance to chemotherapy ((
In addition, the inventors profiled a pair of luminal ER+ breast PDXs: HBCx-22, responsive to Tamoxifen and HBCx-22-TamR, a tumor derivative with previously characterized acquired resistance to Tamoxifen. To obtain a high average coverage of 10,228 reads per cell, the inventors limited the number of encapsulated cells (n=822 tumor cells,
Tumor cells displayed two major chromatin profiles related to their tumor of origin.
However, 16% (n=41 out of 255) of cells within the sensitive tumor shared chromatin features with all resistant cells (
Differential analysis of chromatin features revealed that resistant-like cells have predominantly lost H3K27me3 marks compared to sensitive-like cells (
Loci specifically devoid of H3K27me3 in cells from Chrom_c2 were enriched in genes targets of the Polycomb complex, and characteristic of basal-like signatures of the mammary epithelium (
With scRNA-seq, the inventors could only detect transcription in 2% of differentially enriched windows, and significant differential expression for 3 genes, all showing transcription activation in a fraction of resistant cells in mirror to their loss of H3K27me3 enrichment: EGFR, a gene implicated in resistance to Tamoxifen (Massarweh, S. et al. Cancer Res; 2008, 68:826-833; Ciupek, A. et al. Breast Cancer Res Treat; 2015, 154:225-237), IGFBP3 and ALCAM (
Parallel scRNA-seq analysis of the same samples revealed several clusters of cells within the resistant and sensitive tumor (
Conclusions
Profiling histone modifications at the single-cell level with high coverage, up to 10,000 loci in average per cell, was instrumental to reveal the presence of relatively rare chromatin states within tumor samples. This study suggests that rare cells with chromatin features characteristic of resistant cancer cells exist before treatment and could be selected for by cancer therapy.
Spontaneous heterogeneity of chromatin states in untreated cells may be a key molecular component in the acquisition of drug-resistance, regardless of the mechanism of action of the cancer treatment: here Tamoxifen, targeting the estrogen receptor, and Capecitabine, a classical chemotherapy inhibiting the synthesis of thymidine monophosphate. This study permits to discover new biomarkers for patient stratification and opens up perspectives for novel therapeutic strategies both in luminal and triple-negative breast cancer to counteract resistance. For example, preventing loss of repressive chromatin marks such as H3K27me3, as observed in resistant cells here, by combining treatment with Capecitabine and Tamoxifen with drugs such as demethylase inhibitors could be a strategy to consider to minimize resistance.
Number | Date | Country | Kind |
---|---|---|---|
19305527.4 | Apr 2019 | EP | regional |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2020/061459 | 4/24/2020 | WO | 00 |