Qinhan Li
Zhenan Zhang
Yu Fan
Qian Zhang- Department of Urology, Peking University First Hospital, Institute of Urology, National Research Center for Genitourinary Oncology, Peking University, Beijing, China
The appearance of tyrosine kinase inhibitors (TKIs) has been a major breakthrough in renal cell carcinoma (RCC) therapy. Unfortunately, a portion of patients with TKIs resistance experience disease progression after TKIs therapy. Epigenetic alterations play an important role in the development of TKIs resistance. Current evidence suggests that epigenetic alterations occur frequently in RCC patients with poor response to TKIs therapy, and modulation of them could enhance the cytotoxic effect of antitumor therapy. In this review, we summarize the currently known epigenetic alterations relating to TKIs resistance in RCC, focusing on DNA methylation, non-coding RNAs (ncRNAs), histone modifications, and their interactions with TKIs treatment. In addition, we discuss application of epigenetic alteration analyses in the clinical setting to predict prognosis of patients with TKIs treatment, and the potential use of epigenetics-based therapies to surmount TKIs resistance.
Introduction
Renal cell carcinoma (RCC) is the most common type of renal cancer, causing more than 14,000 deaths yearly (Capitanio et al., 2019). For early stage of RCC, surgical excision is the recommended treatment. However, there are nearly 15% of patients with distant metastasis when diagnosed with RCC (Siegel et al., 2019).
Angiogenesis plays an important role in the biology and the pathogenesis of RCC. Loss of function of von Hippel–Lindau (VHL) tumor suppressor gene is a vital event in renal carcinogenesis and occurs in about 90% of all clear cell renal cell cancer (ccRCC; Nickerson et al., 2008). VHL encodes and forms a VHL protein complex, which acts as an essential factor in the oxygen-sensing pathway through ubiquitin-mediated degradation of hydroxylated hypoxia inducible factor 1 (HIF-1α) and HIF-2α (Maxwell et al., 1999; Kaelin, 2002). Loss of VHL function leads to the accumulation of HIF-1α and HIF-2α, which consequently facilitates transcription of the hypoxia response genes, such as genes in vascular endothelial growth factor (VEGF), platelet-derived growth factor (PDGF), and transforming growth factor alpha (TGF-α), eventually, resulting in angiogenesis and progression of tumor (Kourembanas et al., 1990; de Paulsen et al., 2001). The expression of VEGF and PDGF is significantly upregulated in RCC as a result of VHL inactivation, which, on the one hand, accelerates growth of tumor, on the other hand, is also its weakness. Tyrosine kinase inhibitors (TKIs), including sunitinib, pazopanib, axitinib, sorafenib, and cabozanitinib are thought to exert their major therapeutic effects in RCC by antagonism of VEGF receptor (VEGFR) and PDGF receptor (PDGFR), leading to a reduction of tumor angiogenesis.
For metastatic RCC (mRCC), sunitinib, pazopanib, and cabozantinib are approved for first-line treatment, while axitinib and sorafenib are chosen as second-line treatment. Sunitinib is the most commonly used TKIs which can delay tumor progression and improve patient survival. However, only 20–30% of patients respond to sunitinib treatment initially, and almost all initial responders develop resistance in 2 years (Morais, 2014). Subsequent antitumor therapies are followed by immune-checkpoint inhibitor and mammalian target of rapamycin (mTOR), such as nivolumab and everolimus. TKIs resistance poses a great challenge for the TKIs treatment. Therefore, understanding the distinct molecular mechanisms underlying TKIs resistance is vital to find efficient biomarkers to predict the effect of TKIs and facilitate the development of novel antitumor drugs which overcome this resistance.
An Overview of Epigenetic Modification
Epigenetics refers to the study of molecules and mechanisms that can control chromatin structure and influence gene expression or the propensity for genes to be transcribed within organisms in the context of the same DNA sequence. The ability of cells to retain and transmit their special gene expression patterns to the progeny cells, referred to as epigenetic memory, is governed by epigenetic marks, such as DNA methylations, histone modifications, and non-coding RNAs (ncRNAs; Thiagalingam, 2020). Epigenetic modification is heritable but reversible (Cavalli and Heard, 2019). The unique epigenome defining the genetic code associated with each individual gene regulates the expression status of that gene. Defects in epigenetic factors and epigenetic modifications could act as pushers for various diseases including cancer.
Epigenetic modification is associated with drug resistance in numerous types of cancer, including RCC (Chekhun et al., 2007; Knoechel et al., 2014; Adelaiye-Ogala et al., 2017; Leonetti et al., 2019), which regulates gene expression at the protein level (histone modification and nucleosome remodeling), DNA level (DNA methylation), and RNA level (ncRNA). Histones are the central component of nucleosomal subunit, including four types of histone proteins [histone 2A (H2A), H2B, H3, and H4], which are wrapped by a 147-base-pair segment of DNA (Santos-Rosa and Caldas, 2005; Audia and Campbell, 2016). Histone modifications mainly take place at histone tails, which are densely populated with basic lysine and arginine residues (Audia and Campbell, 2016). The acetylation and methylation of lysine residues are well-known. Acetylation can alter the charge on the lysine residues and weaken the interaction of these histones with DNA, making the chromatin structure more open and accessible (Dawson and Kouzarides, 2012). This process is regulated by two enzymatic families with competing activities: promoted by histone lysine acetyltransferases (HATs) and inhibited by the histone deacetylases (HDACs; Li et al., 2019). Methylation of lysine residues in histone tails contains three forms: monomethylation (me1), demethylation (me2), and trimethylation (me3), making activation or repression of transcription (Kouzarides, 2007). This process is also competitively regulated by histone lysine methyltransferases (KMTs) and histone lysine demethylases (KDMs). Generally, acetylation of lysine 14 of H3 (H3K14), monmethylation of H3K4, H3K9, and H3K79, and phosphorylation of serine 10 (H3S10) are all linked with transcriptional activity (Cheung et al., 2000; Lo et al., 2000; Barski et al., 2007), while trimethylation of H3K9, H3K79, and H3K27 marks transcriptional repression (Boyer et al., 2006; Barski et al., 2007).
At the DNA level, the methylation of the 5-carbon on cytosine CpG dinucleotides is considered as an important epigenetic marker. Catalyzed by DNA methyltransferases (DNMTs), 5-carbon of the cytosine ring on promoter CpG islands gets a methyl from S-adenosylmethionine, converting to 5-methycytosine (5mc). 5mc attracts HDACs and methy-CpG-binding domain proteins (MBDs) to the site, resulting in removal of acetyl groups from histone proteins, compact conformation of nucleosome, and downregulation of gene transcription (Robert et al., 2003; Feinberg and Tycko, 2004; Wang et al., 2009). This process can be reversed by ten-eleven translocation (TET) proteins, which oxidize 5mc into 5-hydroxymethylcytosine (5hmc) and subsequently into 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC) in an Fe(II)- and 2-oxoglutarate-dependent dioxygenases manner (Wu and Zhang, 2014; Rasmussen and Helin, 2016).
Non-coding RNAs regulate gene expression at RNA level, including mircoRNAs (miRNAs), small nucleolar RNAs (snoRNAs), piwiRNAs (piRNAs), and long ncRNA (lncRNA; Ma et al., 2013). In general, composed by about 19–25 nucleotides, miRNAs can lead to posttranscriptional gene silencing and translation stopping through binding to the 3'-untranslated region (3'-UTR) of the targeted messenger RNAs (mRNAs) and leading to its degradation or destabilization (Brennecke et al., 2005). lncRNAs are collectively defined as longer than 200 nucleotides in length, which modulate local or global gene expression in a neighboring (cis) or distal (trans) manner (Kopp and Mendell, 2018). For example, one classic cis-acting lncRNA is the X-inactive specific transcript (Xist) resulting in the X chromosome inaction (XCI) in mammals by recruiting various protein complexes to specific position (Lee and Bartolomei, 2013). Notably, lncRNAs can function as competing endogenous RNAs (ceRNAs) to compete with miRNAs by binding to their protein-coding transcripts, thereby antagonizing the repressive effects of miRNAs on mRNAs (Salmena et al., 2011; Du et al., 2016).
In this review, we summarize the currently known epigenetic alterations relating to TKIs resistance in RCC, focusing on DNA methylation, ncRNAs, histone modifications, and their interactions with TKIs treatment. In addition, we discuss the application of epigenetic alteration analyses in the clinical setting to predict prognosis of patients with TKIs treatment and develop new agents.
Mechanisms of Primary and Acquired Resistance to TKIs Treatment
There is no specific definition of TKIs resistance in RCC. Response to drug therapy is normally defined by the Response Evaluation Criteria In Solid Tumors (RECIST) criteria as evidence of tumor progression regardless of persistent treatment. Unfortunately, the current clinical studies depended on their own criteria to divide patients into responders and non-responders, which made their outcomes difficult to compare. Resistance to antiangiogenic therapy can be classified into intrinsic (primary) and acquired (secondary) resistance (Mollica et al., 2019). Intrinsic resistance is defined as an initial inefficacy of therapeutic agents, which may be attributed to the presence of resistant tumor clones prior to therapy due to inherited resistance or evolutionary clonal selection. Acquired resistance is classified as the progression of tumor after initial tumor regression during the therapy, which is often driven by the development of other pathways stimulating angiogenesis, such as AXL, MET, and PDGF/PDGFR, and thus the escape of cancer cells from VEGF/VEGFR blockade (Crawford et al., 2009; Zhou et al., 2016). While the explicit mechanisms of TKIs resistance are still being explored, several potential factors have been reported to be associated with TKIs resistance in RCC: lysosomal sequestration, mutations and modification of expression level, downstream signaling pathway activation, bypass or alternative pathway activation, ATP-binding cassette (ABC) efflux transporters, tumor microenvironment, epithelial-mesenchymal transition (EMT), and epigenetic modification (Housman et al., 2014; Makhov et al., 2018). Epigenetic regulation of the TKIs resistance is always linked to activation of downstream signaling pathways, promotion of EMT, and stimulation of alternative pathways.
Modulation of Downstream Signaling Pathways
Tyrosine kinase inhibitors exert their major antiangiogenic and antitumor effect in RCC by suppressing tyrosine kinase receptors on VEGFR and PDGFR and inhibiting their downstream signaling pathways. Therefore, RCC cells escape TKIs blockade through an important mechanism of activation of parallel downstream signaling pathways, among which PI3K/AKT and RAS/RAF/ERK are pivotal transduction cascades responsible for cell survival, proliferation, and invasion (Fresno Vara et al., 2004; Guo et al., 2015; Huang and Fu, 2015). The PI3K/AKT pathway is frequently activated in cancer and leads to the development and progression of numerous tumor types, including RCC (Samuels et al., 2004; Lawrence et al., 2014). PI3K, a family of lipid kinases, is normally activated by extracellular signals, such as growth factors, cytokines in physiologic conditions. Activated PI3K phosphorylates phosphatidylinositol 4,5-bisphosphate [Ptdlns(4,5)P2], propagating activation signals to downstream molecules (Hennessy et al., 2005). Phosphatase and tensin homolog deleted on chromosome 10 (PTEN) can turn off this pathway by inhibiting the phosphorylation of Ptdlns(4,5)P2 (Gewinner et al., 2009). AKT is the key mediator to respond to the PI3K signaling. The phosphorylated active AKT translocates from the cell membrane to other cell compartments to phosphorylate multiple downstream substrates, resulting in cell survival, growth, tumorigenesis, metastasis, and sunitinib resistance (Andjelković et al., 1997; Sakai et al., 2013; Fang et al., 2019). Activated by pAKT, mTOR complex 1 can lead to protein translation and lipid or nucleotide synthesis via phosphorylating numerous substrates, such as p70 ribosomal S6 kinase (p70S6K) and Eif4e-binding proteins (Manning and Cantley, 2003; Fruman and Rommel, 2014), eventually leading to the translation and accumulation of HIF-1α and HIF-2α. Acting as an inhibitory protein of the pathway, PTEN contributes to the downregulation of AKT activity, and loss of PTEN leads to sunitinib resistance due to lack of inhibitory input (Makhov et al., 2012). Sekino et al. (2019) identified miR-130 upregulation was associated with sunitinib resistance through suppression of PTEN.
Focal adhesion kinase (FAK) signaling plays an important role in activation of PI3K/AKT pathway by interacting with PI3K (Zhao and Guan, 2009; Poettler et al., 2013; Hung et al., 2017). Activation of FAK signaling contributes to the sorafenib and sunitinib resistance in a variety type of cancer, including RCC (Bai et al., 2012; Zhang et al., 2016; Zhou et al., 2017). The chromatin modifier enhancer of zeste homolog 2 (EZH2), a polycomb group protein homolog of Drosophila enhancer of zeste, is a histone methytransferase unit of polycomb repressive (PRC2), which can catalyze the trimethylation of H3K27, change chromatin configuration, and promote transcriptional silencing (Margueron and Reinberg, 2011; Di Croce and Helin, 2013). Adelaiye-Ogala et al. (2017) reported that increased EZH2 was associated with sunitinib resistance through redistribution in RCC cells, decreasingly binding to the PTK2 gene, which encodes the FAK, and increasingly binding to DAB2IP and PTPN2, which act as tumor suppressors to inhibit RAS/RAF/ERK and P13K/AKT signaling pathways.
Ras/Raf/ERK signaling pathway is other important transduction cascade transmitting EGFR signaling, responsible for cancer development, maintenance, progression and thus, poorer prognosis and TKIs resistance (Bridgeman et al., 2016; Mandal et al., 2016). The methylation of glutaminyl peptide cyclotransferase (QPCT) gene had been reported to associate with sunitinib resistance through Ras/Raf/ERK signaling pathway (Zhao et al., 2019). The QPCT gene encodes glutaminyl cyclase (QC), an enzyme that is involved in the posttranslational modification by converting the N-terminal glutaminyl and glutamyl into pyroglutamate through cyclization, making the protein more resistant to protease degradation, more hydrophobic, and more prone to aggregation and neurotoxicity (Khan et al., 2016; Vijayan and Zhang, 2019). Hypomethylated QPCT gene increased the expression of QC, the process promoted by the NF-κB signaling (p65; Kehlen et al., 2013), leading to upregulation of HRAS and activation of the Ras/Raf/ERK signaling pathway (Herrero et al., 2016; Michael et al., 2016; Zhao et al., 2019). Zhai et al. (2017) had observed that lncRNA-SARCC could regulate androgen receptor (AR) to increase miR-143-3p expression and inhibit its downstream signals, including AKT, MMP-13, K-RAS, and P-ERK. The expression of lncRNA-SARCC was upregulated in RCC cells treated with sunitinib, which was associated with decreased resistance to sunitinib.
Modulation of Epithelial-to-Mesenchymal Transition
Epithelial-to-mesenchymal transition is a biologic process that epithelial cells lose their cell–cell basement membrane contacts and their structural polarity to become spindle-shaped and morphologically similar to mesenchymal cell (He and Magi-Galluzzi, 2014). While potential mechanisms are not fully explicit, numerous studies indicate that EMT constitutes a relevant resistance mechanism to TKI treatment (Fang et al., 2019; Hwang et al., 2019; Zhu et al., 2019), and relates to the development of metastases in cancer (Bastid, 2012). Signal transduction affects EMT through the TGF-beta 1 (TGF-β1) in different mechanisms (Wendt et al., 2009; Feldkoren et al., 2017; Fardi et al., 2019). Schematically, TGF-β1 activates zinc finger E-box binding 1 (ZEB1) and ZEB2, which are responsible for a key transcriptional repressor of the cadherin 1 gene (CDH1). CDH1 encodes the cell-adhesion glycoprotein E-cadherin whose downregulation is a pivotal hallmark of EMT (Loh et al., 2019). As an activator of EMT, the expression of ZEB2 is regulated by miR-141 (Berkers et al., 2013). In detail, miR-141 downregulation induces EMT and hypoxia resistance through the upregulation of ZEB2 and suppression of E-cadherin, resulting in an unfavorable response to sunitinib resistance and poor prognosis (Berkers et al., 2013; Fang et al., 2013).
The overexpression of EZH2 is beneficial to EMT by repression of E-cadherin (Crea et al., 2012; Liu et al., 2016). Adelaiye-Ogala et al. (2017) reported that EZH2 expression was linked to sunitinib resistance in RCC through an adaptive kinome reprograming, such as increased global tyrosine and serine phosphorylation as well as increased phosphorylated FAK. SOX5, one of SOX family involving in the regulation of tumor progression, is thought to contribute to EMT in different types of cancer (Grimm et al., 2019). Liu et al. (2019) reported lncRNA-GAS5 was responsible to sorafenib resistance by functioning as ceRNA to repress miR-21, which controlled its downstream target SOX5.
The Wnt/β-catenin pathway acts as one of the signaling pathways controlling EMT through directly or indirectly targeting several key transcription factors regulating E-cadherin expression and/or the fate of other epithelial molecules (Valenta et al., 2012). SET and MYND domain-containing protein 2 (SMYD2), which acts as one of the SMYD-methyltransferase protein family and specifically methylates H3K4 through its SET domain (Abu-Farha et al., 2008), is deemed to regulate the expression of miR-125b (Yan et al., 2019). miR-125b bind directly to the 3'-UTR of DKK3, a key regulatory factor in the Wnt/β-catenin pathway which acts as a tumor suppressor in RCC (Lu et al., 2017). Thus, the activation of SMYD2/miRNA-EMT pathway weakens the effect sunitinib treatment and accelerates the tumor growth (Yan et al., 2019).
Activation of Bypass Pathways
Extra activation of bypass pathways driving angiogenesis is also one of the most important processes driving TKIs resistance. The activation of MET and AXL confers to the stimulation of their downstream signal cascades, including PI3K and RAS signaling pathway, resulting in sunitinib resistance (Huang and Fu, 2015; Zhou et al., 2016). lncRNA Activated in RCC with Sunitinib Resistance (lncARSR) functions as a sponge and competes for binding of miR-34 and miR-449 to their transcripts, leading to the upregulation of AXL/MET and the activation of STAT3, AKT, and ERK signaling (Qu et al., 2016). miR-32-5p can increase the efficacy of sunitinib by suppressing the testicular nuclear receptor 4 (TR4), which plays an important role in activation of HGF/MET signaling pathway (Wang et al., 2018).
Inactivation of VHL leads to increased HIF-1α and HIF-2α. In renal carcinogenesis, HIF-1α functions more as a tumor suppressor than a tumor promoter, whereas HIF-2α is deemed to predominantly promote tumor growth and angiogenesis (Raval et al., 2005). Specifically, HIF-1α inhibits interaction of MYC with its DNA-binding partners by displacing the SP1 transcription factor from MYC, while HIF-2α could enhance MYC activity by forming a complex with MAX, and thus stabilizing the MYC-MAX and MYC-MAX-SP1 complexes (Keith et al., 2011). HIF-2α/C-MYC axis relates to progression and TKIs resistance in RCC (Zhai et al., 2016; Maroto et al., 2017). Beuselinck et al. (2015) segregated specific groups of patients with ccRCC, who presented sunitinib resistance into four molecular tumor subtypes based on their mRNA expression data: ccrcc1 (c-myc-up), ccrcc2 (classical), ccrcc3 (normal-like), and ccrcc4 (c-myc-up and immune-up). ccrcc1/ccrcc4 subtypes posed a hypomethylation of MYC gene and a global hypermethylation level, with overexpression of MYC and down-expression of corresponding genes, such as PRC2 and SUZ12. Obviously, those two subgroups of patients experienced poor response to sunitinib treatment and shorter progression-free survival (PFS). Verbiest et al. (2018) reported the similar outcome in their study, which proved the resistance to pazopanib in ccrcc1/ccrcc4 subtypes.
In addition, lncRNA-SRLR overexpression is linked to sorafenib resistance through promotion of IL-6 transcription and activation of STAT3 (Xu et al., 2017). miR-942 is associated with sunitinib resistance by promoting the secretion of MMP9 and VEGF (Prior et al., 2014). miR-99a-3p, which targets ribonucleotide reductase regulatory subunit-M2 (RRM2), is downregulated in sunitinib-resistance RCC (Osako et al., 2019). Overexpression of breast cancer resistance protein BCRP/ABCG2, which is posttranscriptionally suppressed by miR-212-3p and miR-132-3p, is associated with superior response to sunitinib treatment in RCC patients (Reustle et al., 2018).
Some of the previously described epigenetic alterations associated with TKIs resistance are represented in Figure 1.
Figure 1. Mechanisms of tyrosine kinase inhibitor (TKI) resistance mediated by epigenetic alterations in renal cell carcinoma (RCC): vascular endothelial growth factor (VEGF) and platelet-derived growth factor (PDGF) bind to a tyrosine kinase receptor and activate the downstream focal adhesion kinase (FAK), PI3K, and RAS pathways. Activated FAK, PI3K, and RAS phosphorylate their downstream signaling cascade, eventually leading to the translation and accumulation of hydroxylated hypoxia inducible factor 1 (HIF-1α) and HIF-2α. In RCC, TKIs exert their influence on antiangiogenesis through inhibition of tyrosine kinase receptor. Epithelial-mesenchymal transition (EMT), activation of downstream signaling pathways and bypass pathways mediated by epigenetic alterations are responsible for the TKIs resistance. Long non-coding RNA (lncRNA)-SARCC increases miR-143-3p expression, thus inhibiting its downstream signals, including AKT, RAS, and ERK. miR-130 enhances HIF signaling by inhibition of PTEN. Hypermethylated QPCT reduces its protein level, leading to inhibition of RAS. EMT, a key transformation in TKIs resistance, is promoted by SOX5, zinc finger E-box binding 2 (ZEB2), and β-Catenin while inhibited by E-cadherin. lncRNA-GAS5 promotes EMT by competing with miR-21 which suppresses the expression of SOX5. miR-141 suppresses the expression of ZEB2 to inhibit its promotion of EMT. SET and MYND domain-containing protein 2 (SMYD2) leads to EMT by promoting the expression of miR-125b, which inhibits DDK3 and activates Wnt/β-catenin signaling pathway. The chromatin modifier enhancer of zeste homolog 2 (EZH2) can not only inhibit E-cadherin but also activate FAK signaling pathway to exert its influence on TKIs resistance. Except for VEGF receptor (VEGFR) and PDGF receptor (PDGFR), activation of MET, AXL, and IL-6 pathways can also lead to phosphorylation of downstream transduction cascades, such as PI3K, STAT3, and RAS. lncRNA Activated in RCC with Sunitinib Resistance (lncARSR) inhibits miR-34 and miR-449, and thus activates MET/AXL pathway. miR-32-5p inhibits MET pathway while lncRNA-SRLR activates interleukin-6 (IL-6)R pathway. RCC, renal cell carcinoma; VEGF, vascular endothelial growth factor; PDGF, platelet-derived growth factor; TKI, tyrosine kinase inhibitor; QPCT, the methylation of glutaminyl peptide cyclotransferase; AR, androgen receptor; EMT, epithelial-to-mesenchymal transition; EZH2, the chromatin modifier enhancer of zeste homolog 2; ZEB2, zinc finger E-box binding 2; SMYD2, SET and MYND domain-containing protein 2; GAS6, growth-arrest-specific protein 6; HGF, hepatocyte growth factor; IL-6, interleukin-6.
Clinical Implications of Epigenetics Analysis in Rcc
Tyrosine kinase inhibitors treatment has been established as first-line therapy for mRCC for a decade with 70–80% of disease control rate. However, approximately 20–30% of patients does not respond to TKIs treatment and experience disease progression within ≤3 months (Porta et al., 2012). Epigenetic alteration can act as a biomarker, which predicts the response of patient to antiangiogenic therapy, thus reducing unnecessary toxicities and costs and maximizing clinical benefit. Clinical investigations of a number of epigenetic alterations on FFPE/plasma samples and their correlation with response to TKIs therapies are listed in Table 1.
Table 1. Epigenetic biomarkers in RCC patients treated with TKIs.
On the histone modification level, tissue low EZH2 expression was associated with increased overall survival (OS) in RCC treated with sunitinib (p = 0.005; Adelaiye-Ogala et al., 2017). On the DNA methylation level, tissue hypomethylation level in the CpG sites of QPCT promoter region showed a poor response to sunitinib therapy (p < 0.05; Zhao et al., 2019). Hypermethylation of cystatin 6 (CST6), ladinin 1 (LAD1) and neurofilament heavy (NEFH) were all linked to shortened PFS (p = 0.009, p = 0.011, and p < 0.001, respectively) and OS (p = 0.011, p = 0.043, and p = 0.028, respectively) for antiangiogenic therapy, including sunitinib, sorafenib, axitinib, and bevacizumab, among which methylation of CST6 could predicted first-line therapy between response (0) and therapy failure (1) with an AUC of 0.88 and a sensitivity and specificity of 82 and 86%, respectively (Dubrowinskaja et al., 2014; Peters et al., 2014). Methylation level of VHL was found to be significantly upregulated after sunitinib therapy (p < 0.001; Stewart et al., 2016), while there was no correlation between VHL methylation and response to pazopanib (Choueiri et al., 2013). Beuselinck et al. (2015) reported the patients with tissue MYC overexpression and global CpG hypermethylation received a shorter PFS and OS after sunitinib treatment (p = 0.001 and p = 0.0003, respectively). In contrast, tissue unmethylation SYNPO2, the gene that encoded myopodin, discriminated progressing patients after TKIs treatment (sunitinib, sorafenib, and pazopanib) from those free of disease, and remained as an independent predictive factor for progression, disease-specific survival, and OS (p = 0.009, p = 0.006, and p = 0.01, respectively; Pompas-Veganzones et al., 2016).
On the ncRNA level, both miRNA and lncRNA showed their influences on the response to TKIs treatment. In their original study, Berkers et al. (2013) described the upregulation of miR-520 g, miR-155, and miR-526b and downregulation of miR-141, miR-376b in tissue were linked to the poor responders to sunitinib (p = 0.036, p = 0.04, p = 0.0067, p = 0.0098, and p = 0.032, respectively). In an observational prospective study, blood samples from 38 patients and 287 miRNAs were taken and evaluated before initiation of therapy and 14 days later in patients receiving sunitinib treatment for advanced RCC. Twenty eight miRNAs of the 287 were found to be significant differences of expressions between the poor response and response groups, among which, downexpression of miR-424 was linked with prolonged response (p = 0.016; Gámez-Pozo et al., 2012). Other researchers (Prior et al., 2014) explored a putative role of miRNAs in influencing sunitinib resistance to RCC in tissue, identifying that tissue overexpressed miR-942 was associated with sunitinib resistance, reduced time to progression (TTP) and OS (p = 0.0074, p = 0.003, and p = 0.0009, respectively), and predicted sunitinib efficacy with an AUC of 0.798 and a sensitivity and specificity of 92 and 50%, respectively. Lukamowicz-Rajska et al. (2016) reported that tissue decreased miR-99b-5p was associated with TKIs non-responders (sunitinib, sorafenib, and pazopanib) with a shorter PFS (<3 months, p < 0.0001). Similarly, Puente et al. (2017) identified that the expression of miR-23b, miR-27b, and miR-628-5p in tumor tissue was upregulated in long-term responders to sunitinib (p < 0.01, each), among which high level of miR-27b and miR-628-5p were associated with increased disease specific survival (p = 0.012 and p = 0.017, respectively). Nineteen miRNAs were explored to have different expressions in tissue, and lower level of miR-155 and miR-484 were associated with increased TTP in patients on sunitinib treatment (p < 0.01 and p < 0.05, respectively; Merhautova et al., 2015). Among 40 miRNAs of 232 found to be downregulated in sunitinib-treated RCC specimens compared with those in normal kidney tissues, miR-101 showed the most dramatic downregulation (p = 0.0013; Goto et al., 2016). Increased miR-9-5p and decreased miR-489-3p were found in non-responder patients of TKIs treatment, including sunitinib, sorafenib, and pazopanib compared to that in responder patients, and the AUC of miR-9-5p combined with clinicopathological variables to predict response(0)/non-response(1) to sunitinib treatment is 0.89 (Ralla et al., 2018). miR-212-3p and miR-132-30 were linked to shorter PFS of sunitinib therapy through interaction with BCRP/ABCG2 expression (Reustle et al., 2018). In a more recent study, high-throughput miRNA microarray performed on FFPE tumor specimens from 47 patients treated with sunitinib, 158 miRNAs were identified to have different expressions in patient with good and poor response (p < 0.05). Moreover, miR-376b was significantly upregulated in patients with a long-term response to sunitinib and could identify patients with long-term response with a sensitivity of 83% and specificity of 67% (p = 0.0002, AUC = 0.758; Kovacova et al., 2019).
Regarding lncRNAs, relevant studies had disclosed their influences on TKIs therapy and prognosis. lncARSR was exposed in plasma and tissue separately by Qu et al. (2016), and was deemed to act as a sponge to compete with miR-34 and miR-449. High level of lncARSR in pre-surgery plasma was an independent prognostic factor for patients with sunitinib treatment, and was correlated with decreased PFS (p = 0.02 and p = 0.014, respectively). Intriguingly, low level of lncARSR in tissue exhibited a superior PFS after receiving sunitinib therapy (p = 0.028). A microarray analysis conducted by Xu et al. (2017) revealed the similar outcome in lncRNA-SRLR. Briefly, high tissue lncRNA-SRLR was associated with poor response to sorafenib, and patients with low lncRNA-SRLR expression had a more significant improvement in PFS after receiving sorafenib treatment (p = 0.0198 and p = 0.0086, respectively). lncRNA-GAS5 was also found to be downregulated in RCC patients with sorafenib resistance (p < 0.01; Liu et al., 2019).
Overall, these studies demonstrate that epigenetic alterations could be promising predictive biomarkers for TKIs response, as they function as important roles involved in mediating resistance through regulating important mechanisms. However, most of these studies are involved in a small number of patients, which limits their reliability. Moreover, there is no accepted criterion on how long a PFS of good response should last, so each study divides patients into good responders and poor responders based on its own standard. The different criteria limit the application of those epigenetic biomarkers in clinical setting.
Epigenetic Alterations as Therapeutic Targets
Besides the predictive value, epigenetic alterations have potential to become targets themselves for drug development, in order to overcome the TKIs resistance in RCC. Preclinical studies on RCC cell lines demonstrate that reversions of epigenetic alterations are effective strategies to re-sensitize resistant clones to TKIs treatment, including demethylation, restoration of miRNA function, and inhibition of HDAC. Therefore, implementing epigenetics-based therapeutic strategies in patients is the next step, and relevant clinical trials are under way. Generally, there are two classes of epigenetics-based drugs in clinical trials: broad reprogrammers, which have a broad effect, and targeted therapies, which focus on specific miRNA expression or histone modifications (Jones et al., 2016). The formers are represented by the DNMT inhibitors (DNMTi) and the HDAC inhibitors (HDACi), and the latters are represented by EZH2 inhibitors (Jones et al., 2016; Joosten et al., 2018). The outcomes of current clinical trials concerning combination of epigenetics-based therapy with TKIs are listed in Table 2. So far, HDACi and EZH2 inhibitors are the most promising agents to reverse the TKI resistance with a vast of clinical studies completed or ongoing. Combination of HDACi and antiangiogenic agents is the most common trial to reverse the acquired resistance and re-sensitize tumors to antiangiogenic therapy. A phase I study evaluated the safety, tolerability, and preliminary efficacy of HDACi vorinostat plus sorafenib in patients with RCC and non-small cell lung cancer (NSCLC) and showed poor tolerance and no confirmed responses (Dasari et al., 2013). Other study focused on the combination vorinostat with pazopanib in advanced solid tumors including RCC, and identified that the treatment achieved stable disease for at least 6 months or partial response (PR; SD ≥ 6 months/PR) in 19% of all patients (n = 78), median PFS of 2.2 months, and median OS of 8.9 months (Fu et al., 2015). Furthermore, patients with detected hotspot TP53 mutations had a superior rate of SD ≥ 6 months/PR, median PFS, and OS compared with those with undetected hotspot TP53 mutations (45 vs. 16%, 3.5 vs. 2.0 months, and 12.7 vs. 7.4 months, respectively). In a phase I study, combination HDACi abexinostat with pazopanib in patient of RCC with tumor progression after received an average 2.5 lines of prior therapy and 1.6 lines of prior VEGF-targeting treatment received 27% of objective response rate and average 10.5 months of response duration (Aggarwal et al., 2017). Three patients with prior refractory disease to pazopanib monotherapy received durable minor or PR > 12 months treated with pazopanib plus abexinostat. Other clinical trial explored the effect of combination of HDACi with monoclonal antibody bevacizumab in advanced RCC. In a multicenter, single-arm phase I/II clinical trial, 33 patients with metastatic or unresectable ccRCC achieved 5.7 months of median PFS and 13.9 months of median OS, among which six patients achieved OR, including 1 CR and 5 PR (Pili et al., 2017). Regarding DNMTi, decitabine was the only agent tested by phase I trial and combination it with high-dose IL-2 achieved stable disease in three patients (Gollob et al., 2006). As the pharmacological defects of this DNMTi, such as short half-life and sensitivity to inactivation by cytidine deaminase, limit their clinical application, the second-generation DNMTi guadecitabine has been developed, which has shown promise in early preclinical models and clinical trial in patients with acute myeloid leukemia and myelodysplastic syndromes (Joosten et al., 2018).
Table 2. Epigenetic drugs plus TKIs in treatment of RCC.
The H3K27 histone N-methyltransferase EZH2 is a pusher of EMT leading to TKIs resistance, and its inhibitors may contribute to re-sensitize tumor to antiangiogenic treatment, which has been proven in preclinical test in RCC lines (Wagener et al., 2008; Adelaiye et al., 2015). The result of phase I trial that 64 patients including 21 with B-cell non-Hodgkin lymphoma and 43 with advanced solid tumors received EZH2 inhibitor tazemetostat showed the agent had a favorable safety profile and antitumor activity (Italiano et al., 2018).
Moreover, miRNAs also have the potential to become a target to reverse the TKIs resistance in RCC. Preclinical studies on RCC lines clearly demonstrated that both restoration of the tumor-suppressor miRNA function (by miRNA mimics) and inhibition of the oncogenic miRNAs (by antagomiRs) could re-sensitize resistance clones to TKIs. However, implementing miRNA-based therapies in clinic constitutes a significant challenge for clinicians and has not yet been realized. The main concerns fasten on the relative instability of miRNAs in body fluids and specific delivery of these miRNAs to tumor sites (Christopher et al., 2016; Leonetti et al., 2019). Recently, exosomes have been identified to function as carriers of miRNAs to deliver them from cell to extracellular milieu, which may become the sally port for miRNA-based therapy (Mathiyalagan and Sahoo, 2017; Rahbarghazi et al., 2019).
In addition, as epigenetic memory defines the ability of cells to retain and transmit their special gene expression status to the daughter cells, one differentiated somatic cell could become pluripotent and subsequently be reprogrammed into a different somatic cell through loss of its epigenetic memories responsible for its differentiated state. This process could serve as the basis for stem cell therapeutics by replacing one’s affected cells with his/her own cells, which may become the potential target of new agents (Thiagalingam, 2020).
Although agents targeting the epigenome could be a promising therapy strategy for TKIs resistance in mRCC because of the widespread epigenetic deregulation in this tumor type, there are several problems of those agents limiting their clinical application. For example, clinical activity of a drug is not only related to the original rationale but also attribute to the off-target effects. Studies about patients treated with epigenetic agents such as DNMTi revealed acute genome-wide demethylation (Yang et al., 2006), which may not only restore abnormally silenced expression but also activate normally silenced expression, leading to adverse off-targets effects. The individual responses of epigenetic agents are variable in different tumor types. So far, hypomethylating drugs are generally more effective in myeloid malignancy than in RCC. Furthermore, the majority of patients have been treated with DNMTi or HDACi for a shortened period of time, the long-term effects of those agents are not explicit for us. In addition, combination treatment might bring more severe and dose-limiting toxicities than monotherapy. As a result, additional trials are urged to future elaborate the interaction of those agents with mechanism of TKIs resistance and to assess their use in RCC patients.
Conclusion
Although the advent of TKIs therapy indeed provides concrete hope for patients with advenced RCC, a part of patients with intrinsic or acquired resistance to TKIs benefit a little from the therapy and experience tumor progression after treatment. Epigenetic alterations are involved in the mechanisms underlying this event and could act as excellent biomakers to predict the response of patients to TKIs treatment. However, no epigenetic biomarker is currently applied in clinical setting regardless of numerous epigenetic biomarkers reported. Low efficiency and high cost of them may be the cause of this event. Therefore, for the purpose of translation them into clinical practice, more high-quality epigenetic biomarker studies are needed. In view, the criteria of TKIs resistance are ambiguous, uniform defination of TKIs resistance is urgent affair. Assay, statistical methods, and study designs also need be standardized to optimize their practice. In addition, epigenetics-based therapies are in full swing, which hold great promise and may optimize the management of patients with advanced RCC.
Author Contributions
QL and ZZ searched the related published articles and wrote the manuscript. QZ and YF designed the work and instructed the progress of the review. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Abu-Farha, M., Lambert, J. P., Al-Madhoun, A. S., Elisma, F., Skerjanc, I. S., and Figeys, D. (2008). The tale of two domains: proteomics and genomics analysis of SMYD2, a new histone methyltransferase. Mol. Cell. Proteomics 7, 560–572. doi: 10.1074/mcp.M700271-MCP200
Adelaiye, R., Ciamporcero, E., Miles, K. M., Sotomayor, P., Bard, J., Tsompana, M., et al. (2015). Sunitinib dose escalation overcomes transient resistance in clear cell renal cell carcinoma and is associated with epigenetic modifications. Mol. Cancer Ther. 14, 513–522. doi: 10.1158/1535-7163.mct-14-0208
Adelaiye-Ogala, R., Budka, J., Damayanti, N. P., Arrington, J., Ferris, M., Hsu, C. C., et al. (2017). EZH2 modifies sunitinib resistance in renal cell carcinoma by kinome reprogramming. Cancer Res. 77, 6651–6666. doi: 10.1158/0008-5472.can-17-0899
Aggarwal, R., Thomas, S., Pawlowska, N., Bartelink, I., Grabowsky, J., Jahan, T., et al. (2017). Inhibiting histone deacetylase as a means to reverse resistance to angiogenesis inhibitors: phase I study of abexinostat plus pazopanib in advanced solid tumor malignancies. J. Clin. Oncol. 35, 1231–1239. doi: 10.1200/jco.2016.70.5350
Andjelković, M., Alessi, D. R., Meier, R., Fernandez, A., Lamb, N. J., Frech, M., et al. (1997). Role of translocation in the activation and function of protein kinase B. J. Biol. Chem. 272, 31515–31524. doi: 10.1074/jbc.272.50.31515
Audia, J. E., and Campbell, R. M. (2016). Histone modifications and cancer. Cold Spring Harb. Perspect. Biol. 8:a019521. doi: 10.1101/cshperspect.a019521
Bai, L., Yang, J. C., Ok, J. H., Mack, P. C., Kung, H. J., and Evans, C. P. (2012). Simultaneous targeting of Src kinase and receptor tyrosine kinase results in synergistic inhibition of renal cell carcinoma proliferation and migration. Int. J. Cancer 130, 2693–2702. doi: 10.1002/ijc.26303
Barski, A., Cuddapah, S., Cui, K., Roh, T. Y., Schones, D. E., Wang, Z., et al. (2007). High-resolution profiling of histone methylations in the human genome. Cell 129, 823–837. doi: 10.1016/j.cell.2007.05.009
Bastid, J. (2012). EMT in carcinoma progression and dissemination: facts, unanswered questions, and clinical considerations. Cancer Metastasis Rev. 31, 277–283. doi: 10.1007/s10555-011-9344-6
Berkers, J., Govaere, O., Wolter, P., Beuselinck, B., Schöffski, P., van Kempen, L. C., et al. (2013). A possible role for microRNA-141 down-regulation in sunitinib resistant metastatic clear cell renal cell carcinoma through induction of epithelial-to-mesenchymal transition and hypoxia resistance. J. Urol. 189, 1930–1938. doi: 10.1016/j.juro.2012.11.133
Beuselinck, B., Job, S., Becht, E., Karadimou, A., Verkarre, V., Couchy, G., et al. (2015). Molecular subtypes of clear cell renal cell carcinoma are associated with sunitinib response in the metastatic setting. Clin. Cancer Res. 21, 1329–1339. doi: 10.1158/1078-0432.ccr-14-1128
Boyer, L. A., Plath, K., Zeitlinger, J., Brambrink, T., Medeiros, L. A., Lee, T. I., et al. (2006). Polycomb complexes repress developmental regulators in murine embryonic stem cells. Nature 441, 349–353. doi: 10.1038/nature04733
Brennecke, J., Stark, A., Russell, R. B., and Cohen, S. M. (2005). Principles of microRNA-target recognition. PLoS Biol. 3:e85. doi: 10.1371/journal.pbio.0030085
Bridgeman, V. L., Wan, E., Foo, S., Nathan, M. R., Welti, J. C., Frentzas, S., et al. (2016). Preclinical evidence that trametinib enhances the response to antiangiogenic tyrosine kinase inhibitors in renal cell carcinoma. Mol. Cancer Ther. 15, 172–183. doi: 10.1158/1535-7163.mct-15-0170
Capitanio, U., Bensalah, K., Bex, A., Boorjian, S. A., Bray, F., Coleman, J., et al. (2019). Epidemiology of renal cell carcinoma. Eur. Urol. 75, 74–84. doi: 10.1016/j.eururo.2018.08.036
Cavalli, G., and Heard, E. (2019). Advances in epigenetics link genetics to the environment and disease. Nature 571, 489–499. doi: 10.1038/s41586-019-1411-0
Chekhun, V. F., Lukyanova, N. Y., Kovalchuk, O., Tryndyak, V. P., and Pogribny, I. P. (2007). Epigenetic profiling of multidrug-resistant human MCF-7 breast adenocarcinoma cells reveals novel hyper- and hypomethylated targets. Mol. Cancer Ther. 6, 1089–1098. doi: 10.1158/1535-7163.mct-06-0663
Cheung, P., Tanner, K. G., Cheung, W. L., Sassone-Corsi, P., Denu, J. M., and Allis, C. D. (2000). Synergistic coupling of histone H3 phosphorylation and acetylation in response to epidermal growth factor stimulation. Mol. Cell 5, 905–915. doi: 10.1016/s1097-2765(00)80256-7
Choueiri, T. K., Fay, A. P., Gagnon, R., Lin, Y., Bahamon, B., Brown, V., et al. (2013). The role of aberrant VHL/HIF pathway elements in predicting clinical outcome to pazopanib therapy in patients with metastatic clear-cell renal cell carcinoma. Clin. Cancer Res. 19, 5218–5226. doi: 10.1158/1078-0432.ccr-13-0491
Christopher, A. F., Kaur, R. P., Kaur, G., Kaur, A., Gupta, V., and Bansal, P. (2016). MicroRNA therapeutics: discovering novel targets and developing specific therapy. Perspect. Clin. Res. 7, 68–74. doi: 10.4103/2229-3485.179431
Crawford, Y., Kasman, I., Yu, L., Zhong, C., Wu, X., Modrusan, Z., et al. (2009). PDGF-C mediates the angiogenic and tumorigenic properties of fibroblasts associated with tumors refractory to anti-VEGF treatment. Cancer Cell 15, 21–34. doi: 10.1016/j.ccr.2008.12.004
Crea, F., Fornaro, L., Bocci, G., Sun, L., Farrar, W. L., Falcone, A., et al. (2012). EZH2 inhibition: targeting the crossroad of tumor invasion and angiogenesis. Cancer Metastasis Rev. 31, 753–761. doi: 10.1007/s10555-012-9387-3
Dasari, A., Gore, L., Messersmith, W. A., Diab, S., Jimeno, A., Weekes, C. D., et al. (2013). A phase I study of sorafenib and vorinostat in patients with advanced solid tumors with expanded cohorts in renal cell carcinoma and non-small cell lung cancer. Investig. New Drugs 31, 115–125. doi: 10.1007/s10637-012-9812-z
Dawson, M. A., and Kouzarides, T. (2012). Cancer epigenetics: from mechanism to therapy. Cell 150, 12–27. doi: 10.1016/j.cell.2012.06.013
de Paulsen, N., Brychzy, A., Fournier, M. C., Klausner, R. D., Gnarra, J. R., Pause, A., et al. (2001). Role of transforming growth factor-alpha in von Hippel--Lindau (VHL)(−/−) clear cell renal carcinoma cell proliferation: a possible mechanism coupling VHL tumor suppressor inactivation and tumorigenesis. Proc. Natl. Acad. Sci. U. S. A. 98, 1387–1392. doi: 10.1073/pnas.031587498
Di Croce, L., and Helin, K. (2013). Transcriptional regulation by Polycomb group proteins. Nat. Struct. Mol. Biol. 20, 1147–1155. doi: 10.1038/nsmb.2669
Du, Z., Sun, T., Hacisuleyman, E., Fei, T., Wang, X., Brown, M., et al. (2016). Integrative analyses reveal a long noncoding RNA-mediated sponge regulatory network in prostate cancer. Nat. Commun. 7:10982. doi: 10.1038/ncomms10982
Dubrowinskaja, N., Gebauer, K., Peters, I., Hennenlotter, J., Abbas, M., Scherer, R., et al. (2014). Neurofilament heavy polypeptide CpG island methylation associates with prognosis of renal cell carcinoma and prediction of antivascular endothelial growth factor therapy response. Cancer Med. 3, 300–309. doi: 10.1002/cam4.181
Fang, Y., Wei, J., Cao, J., Zhao, H., Liao, B., Qiu, S., et al. (2013). Protein expression of ZEB2 in renal cell carcinoma and its prognostic significance in patient survival. PLoS One 8:e62558. doi: 10.1371/journal.pone.0062558
Fang, L., Zhang, Y., Zang, Y., Chai, R., Zhong, G., Li, Z., et al. (2019). HP-1 inhibits the progression of ccRCC and enhances sunitinib therapeutic effects by suppressing EMT. Carbohydr. Polym. 223:115109. doi: 10.1016/j.carbpol.2019.115109
Fardi, M., Alivand, M., Baradaran, B., Farshdousti Hagh, M., and Solali, S. (2019). The crucial role of ZEB2: from development to epithelial-to-mesenchymal transition and cancer complexity. J. Cell. Physiol. doi: 10.1002/jcp.28277 [Epub ahead of print]
Feinberg, A. P., and Tycko, B. (2004). The history of cancer epigenetics. Nat. Rev. Cancer 4, 143–153. doi: 10.1038/nrc1279
Feldkoren, B., Hutchinson, R., Rapoport, Y., Mahajan, A., and Margulis, V. (2017). Integrin signaling potentiates transforming growth factor-beta 1 (TGF-β1) dependent down-regulation of E-cadherin expression—important implications for epithelial to mesenchymal transition (EMT) in renal cell carcinoma. Exp. Cell Res. 355, 57–66. doi: 10.1016/j.yexcr.2017.03.051
Fresno Vara, J. A., Casado, E., de Castro, J., Cejas, P., Belda-Iniesta, C., and González-Barón, M. (2004). PI3K/Akt signalling pathway and cancer. Cancer Treat. Rev. 30, 193–204. doi: 10.1016/j.ctrv.2003.07.007
Fruman, D. A., and Rommel, C. (2014). PI3K and cancer: lessons, challenges and opportunities. Nat. Rev. Drug Discov. 13, 140–156. doi: 10.1038/nrd4204
Fu, S., Hou, M. M., Naing, A., Janku, F., Hess, K., Zinner, R., et al. (2015). Phase I study of pazopanib and vorinostat: a therapeutic approach for inhibiting mutant p53-mediated angiogenesis and facilitating mutant p53 degradation. Ann. Oncol. 26, 1012–1018. doi: 10.1093/annonc/mdv066
Gámez-Pozo, A., Antón-Aparicio, L. M., Bayona, C., Borrega, P., Gallegos Sancho, M. I., García-Domínguez, R., et al. (2012). MicroRNA expression profiling of peripheral blood samples predicts resistance to first-line sunitinib in advanced renal cell carcinoma patients. Neoplasia 14, 1144–1152. doi: 10.1593/neo.12734
Gewinner, C., Wang, Z. C., Richardson, A., Teruya-Feldstein, J., Etemadmoghadam, D., Bowtell, D., et al. (2009). Evidence that inositol polyphosphate 4-phosphatase type II is a tumor suppressor that inhibits PI3K signaling. Cancer Cell 16, 115–125. doi: 10.1016/j.ccr.2009.06.006
Gollob, J. A., Sciambi, C. J., Peterson, B. L., Richmond, T., Thoreson, M., Moran, K., et al. (2006). Phase I trial of sequential low-dose 5-aza-2'-deoxycytidine plus high-dose intravenous bolus interleukin-2 in patients with melanoma or renal cell carcinoma. Clin. Cancer Res. 12, 4619–4627. doi: 10.1158/1078-0432.ccr-06-0883
Goto, Y., Kurozumi, A., Nohata, N., Kojima, S., Matsushita, R., Yoshino, H., et al. (2016). The microRNA signature of patients with sunitinib failure: regulation of UHRF1 pathways by microRNA-101 in renal cell carcinoma. Oncotarget 7, 59070–59086. doi: 10.18632/oncotarget.10887
Grimm, D., Bauer, J., Wise, P., Krüger, M., Simonsen, U., Wehland, M., et al. (2019). The role of SOX family members in solid tumours and metastasis. Semin. Cancer Biol. 67, 122–153. doi: 10.1016/j.semcancer.2019.03.004
Guo, H., German, P., Bai, S., Barnes, S., Guo, W., Qi, X., et al. (2015). The PI3K/AKT pathway and renal cell carcinoma. J. Genet. Genomics 42, 343–353. doi: 10.1016/j.jgg.2015.03.003
He, H., and Magi-Galluzzi, C. (2014). Epithelial-to-mesenchymal transition in renal neoplasms. Adv. Anat. Pathol. 21, 174–180. doi: 10.1097/pap.0000000000000018
Hennessy, B. T., Smith, D. L., Ram, P. T., Lu, Y., and Mills, G. B. (2005). Exploiting the PI3K/AKT pathway for cancer drug discovery. Nat. Rev. Drug Discov. 4, 988–1004. doi: 10.1038/nrd1902
Herrero, A., Casar, B., Colon-Bolea, P., Agudo-Ibanez, L., and Crespo, P. (2016). Defined spatiotemporal features of RAS-ERK signals dictate cell fate in MCF-7 mammary epithelial cells. Mol. Biol. Cell 27, 1958–1968. doi: 10.1091/mbc.E15-02-0118
Housman, G., Byler, S., Heerboth, S., Lapinska, K., Longacre, M., Snyder, N., et al. (2014). Drug resistance in cancer: an overview. Cancer 6, 1769–1792. doi: 10.3390/cancers6031769
Huang, L., and Fu, L. (2015). Mechanisms of resistance to EGFR tyrosine kinase inhibitors. Acta Pharm. Sin. B 5, 390–401. doi: 10.1016/j.apsb.2015.07.001
Hung, T. W., Chen, P. N., Wu, H. C., Wu, S. W., Tsai, P. Y., Hsieh, Y. S., et al. (2017). Kaempferol inhibits the invasion and migration of renal cancer cells through the downregulation of AKT and FAK pathways. Int. J. Med. Sci. 14, 984–993. doi: 10.7150/ijms.20336
Hwang, H. S., Go, H., Park, J. M., Yoon, S. Y., Lee, J. L., Jeong, S. U., et al. (2019). Epithelial-mesenchymal transition as a mechanism of resistance to tyrosine kinase inhibitors in clear cell renal cell carcinoma. Lab. Investig. 99, 659–670. doi: 10.1038/s41374-019-0188-y
Italiano, A., Soria, J. C., Toulmonde, M., Michot, J. M., Lucchesi, C., Varga, A., et al. (2018). Tazemetostat, an EZH2 inhibitor, in relapsed or refractory B-cell non-Hodgkin lymphoma and advanced solid tumours: a first-in-human, open-label, phase 1 study. Lancet Oncol. 19, 649–659. doi: 10.1016/s1470-2045(18)30145-1
Jones, P. A., Issa, J. P., and Baylin, S. (2016). Targeting the cancer epigenome for therapy. Nat. Rev. Genet. 17, 630–641. doi: 10.1038/nrg.2016.93
Joosten, S. C., Smits, K. M., Aarts, M. J., Melotte, V., Koch, A., Tjan-Heijnen, V. C., et al. (2018). Epigenetics in renal cell cancer: mechanisms and clinical applications. Nat. Rev. Urol. 15, 430–451. doi: 10.1038/s41585-018-0023-z
Kaelin, W. G. (2002). Molecular basis of the VHL hereditary cancer syndrome. Nat. Rev. Cancer 2, 673–682. doi: 10.1038/nrc885
Kehlen, A., Haegele, M., Menge, K., Gans, K., Immel, U. D., Hoang-Vu, C., et al. (2013). Role of glutaminyl cyclases in thyroid carcinomas. Endocr. Relat. Cancer 20, 79–90. doi: 10.1530/erc-12-0053
Keith, B., Johnson, R. S., and Simon, M. C. (2011). HIF1alpha and HIF2alpha: sibling rivalry in hypoxic tumour growth and progression. Nat. Rev. Cancer 12, 9–22. doi: 10.1038/nrc3183
Khan, R. A., Chen, J., Shen, J., Li, Z., Wang, M., Wen, Z., et al. (2016). Common variants in QPCT gene confer risk of schizophrenia in the Han Chinese population. Am. J. Med. Genet. B Neuropsychiatr. Genet. 171b, 237–242. doi: 10.1002/ajmg.b.32397
Knoechel, B., Roderick, J. E., Williamson, K. E., Zhu, J., Lohr, J. G., Cotton, M. J., et al. (2014). An epigenetic mechanism of resistance to targeted therapy in T cell acute lymphoblastic leukemia. Nat. Genet. 46, 364–370. doi: 10.1038/ng.2913
Kopp, F., and Mendell, J. T. (2018). Functional classification and experimental dissection of long noncoding RNAs. Cell 172, 393–407. doi: 10.1016/j.cell.2018.01.011
Kourembanas, S., Hannan, R. L., and Faller, D. V. (1990). Oxygen tension regulates the expression of the platelet-derived growth factor-B chain gene in human endothelial cells. J. Clin. Invest. 86, 670–674. doi: 10.1172/jci114759
Kouzarides, T. (2007). Chromatin modifications and their function. Cell 128, 693–705. doi: 10.1016/j.cell.2007.02.005
Kovacova, J., Juracek, J., Poprach, A., Kopecky, J., Fiala, O., Svoboda, M., et al. (2019). MiR-376b-3p is associated with long-term response to sunitinib in metastatic renal cell carcinoma patients. Cancer Genomics Proteomics 16, 353–359. doi: 10.21873/cgp.20140
Lawrence, M. S., Stojanov, P., Mermel, C. H., Robinson, J. T., Garraway, L. A., Golub, T. R., et al. (2014). Discovery and saturation analysis of cancer genes across 21 tumour types. Nature 505, 495–501. doi: 10.1038/nature12912
Lee, J. T., and Bartolomei, M. S. (2013). X-inactivation, imprinting, and long noncoding RNAs in health and disease. Cell 152, 1308–1323. doi: 10.1016/j.cell.2013.02.016
Leonetti, A., Assaraf, Y. G., Veltsista, P. D., El Hassouni, B., Tiseo, M., and Giovannetti, E. (2019). MicroRNAs as a drug resistance mechanism to targeted therapies in EGFR-mutated NSCLC: current implications and future directions. Drug Resist. Updat. 42, 1–11. doi: 10.1016/j.drup.2018.11.002
Li, P., Ge, J., and Li, H. (2019). Lysine acetyltransferases and lysine deacetylases as targets for cardiovascular disease. Nat. Rev. Cardiol. 17, 96–115. doi: 10.1038/s41569-019-0235-9
Liu, L., Pang, X., Shang, W., Xie, H., Feng, Y., and Feng, G. (2019). Long non-coding RNA GAS5 sensitizes renal cell carcinoma to sorafenib via miR-21/SOX5 pathway. Cell Cycle 18, 257–263. doi: 10.1080/15384101.2018.1475826
Liu, L., Xu, Z., Zhong, L., Wang, H., Jiang, S., Long, Q., et al. (2016). Enhancer of zeste homolog 2 (EZH2) promotes tumour cell migration and invasion via epigenetic repression of E-cadherin in renal cell carcinoma. BJU Int. 117, 351–362. doi: 10.1111/bju.12702
Lo, W. S., Trievel, R. C., Rojas, J. R., Duggan, L., Hsu, J. Y., Allis, C. D., et al. (2000). Phosphorylation of serine 10 in histone H3 is functionally linked in vitro and in vivo to Gcn5-mediated acetylation at lysine 14. Mol. Cell 5, 917–926. doi: 10.1016/s1097-2765(00)80257-9
Loh, C. Y., Chai, J. Y., Tang, T. F., Wong, W. F., Sethi, G., Shanmugam, M. K., et al. (2019). The E-cadherin and N-cadherin switch in epithelial-to-mesenchymal transition: signaling, therapeutic implications, and challenges. Cell 8:1118. doi: 10.3390/cells8101118
Lu, Y., Zhao, X., Liu, Q., Li, C., Graves-Deal, R., Cao, Z., et al. (2017). lncRNA MIR100HG-derived miR-100 and miR-125b mediate cetuximab resistance via Wnt/β-catenin signaling. Nat. Med. 23, 1331–1341. doi: 10.1038/nm.4424
Lukamowicz-Rajska, M., Mittmann, C., Prummer, M., Zhong, Q., Bedke, J., Hennenlotter, J., et al. (2016). MiR-99b-5p expression and response to tyrosine kinase inhibitor treatment in clear cell renal cell carcinoma patients. Oncotarget 7, 78433–78447. doi: 10.18632/oncotarget.12618
Ma, L., Bajic, V. B., and Zhang, Z. (2013). On the classification of long non-coding RNAs. RNA Biol. 10, 925–933. doi: 10.4161/rna.24604
Makhov, P. B., Golovine, K., Kutikov, A., Teper, E., Canter, D. J., Simhan, J., et al. (2012). Modulation of Akt/mTOR signaling overcomes sunitinib resistance in renal and prostate cancer cells. Mol. Cancer Ther. 11, 1510–1517. doi: 10.1158/1535-7163.mct-11-0907
Makhov, P., Joshi, S., Ghatalia, P., Kutikov, A., Uzzo, R. G., and Kolenko, V. M. (2018). Resistance to systemic therapies in clear cell renal cell carcinoma: mechanisms and management strategies. Mol. Cancer Ther. 17, 1355–1364. doi: 10.1158/1535-7163.mct-17-1299
Mandal, R., Becker, S., and Strebhardt, K. (2016). Stamping out RAF and MEK1/2 to inhibit the ERK1/2 pathway: an emerging threat to anticancer therapy. Oncogene 35, 2547–2561. doi: 10.1038/onc.2015.329
Manning, B. D., and Cantley, L. C. (2003). Rheb fills a GAP between TSC and TOR. Trends Biochem. Sci. 28, 573–576. doi: 10.1016/j.tibs.2003.09.003
Margueron, R., and Reinberg, D. (2011). The Polycomb complex PRC2 and its mark in life. Nature 469, 343–349. doi: 10.1038/nature09784
Maroto, P., Esteban, E., Parra, E. F., Mendez-Vidal, M. J., Domenech, M., Pérez-Valderrama, B., et al. (2017). HIF pathway and c-Myc as biomarkers for response to sunitinib in metastatic clear-cell renal cell carcinoma. Onco Targets Ther. 10, 4635–4643. doi: 10.2147/OTT.S137677
Mathiyalagan, P., and Sahoo, S. (2017). Exosomes-based gene therapy for microRNA delivery. Methods Mol. Biol. 1521, 139–152. doi: 10.1007/978-1-4939-6588-5_9
Maxwell, P. H., Wiesener, M. S., Chang, G. W., Clifford, S. C., Vaux, E. C., Cockman, M. E., et al. (1999). The tumour suppressor protein VHL targets hypoxia-inducible factors for oxygen-dependent proteolysis. Nature 399, 271–275. doi: 10.1038/20459
Merhautova, J., Hezova, R., Poprach, A., Kovarikova, A., Radova, L., Svoboda, M., et al. (2015). miR-155 and miR-484 are associated with time to progression in metastatic renal cell carcinoma treated with sunitinib. Biomed. Res. Int. 2015:941980. doi: 10.1155/2015/941980
Michael, J. V., Wurtzel, J. G., and Goldfinger, L. E. (2016). Inhibition of galectin-1 sensitizes HRAS-driven tumor growth to rapamycin treatment. Anticancer Res. 36, 5053–5061. doi: 10.21873/anticanres.11074
Mollica, V., Di Nunno, V., Gatto, L., Santoni, M., Scarpelli, M., Cimadamore, A., et al. (2019). Resistance to systemic agents in renal cell carcinoma predict and overcome genomic strategies adopted by tumor. Cancer 11:830. doi: 10.3390/cancers11060830
Morais, C. (2014). Sunitinib resistance in renal cell carcinoma. J. Kidney Cancer VHL 1, 1–11. doi: 10.15586/jkcvhl.2014.7
Nickerson, M. L., Jaeger, E., Shi, Y., Durocher, J. A., Mahurkar, S., Zaridze, D., et al. (2008). Improved identification of von Hippel-Lindau gene alterations in clear cell renal tumors. Clin. Cancer Res. 14, 4726–4734. doi: 10.1158/1078-0432.ccr-07-4921
Osako, Y., Yoshino, H., Sakaguchi, T., Sugita, S., Yonemori, M., Nakagawa, M., et al. (2019). Potential tumor-suppressive role of microRNA-99a-3p in sunitinib-resistant renal cell carcinoma cells through the regulation of RRM2. Int. J. Oncol. 54, 1759–1770. doi: 10.3892/ijo.2019.4736
Peters, I., Dubrowinskaja, N., Abbas, M., Seidel, C., Kogosov, M., Scherer, R., et al. (2014). DNA methylation biomarkers predict progression-free and overall survival of metastatic renal cell cancer (mRCC) treated with antiangiogenic therapies. PLoS One 9:e91440. doi: 10.1371/journal.pone.0091440
Pili, R., Liu, G., Chintala, S., Verheul, H., Rehman, S., Attwood, K., et al. (2017). Combination of the histone deacetylase inhibitor vorinostat with bevacizumab in patients with clear-cell renal cell carcinoma: a multicentre, single-arm phase I/II clinical trial. Br. J. Cancer 116, 874–883. doi: 10.1038/bjc.2017.33
Poettler, M., Unseld, M., Braemswig, K., Haitel, A., Zielinski, C. C., and Prager, G. W. (2013). CD98hc (SLC3A2) drives integrin-dependent renal cancer cell behavior. Mol. Cancer 12:169. doi: 10.1186/1476-4598-12-169
Pompas-Veganzones, N., Sandonis, V., Perez-Lanzac, A., Beltran, M., Beardo, P., Juárez, A., et al. (2016). Myopodin methylation is a prognostic biomarker and predicts antiangiogenic response in advanced kidney cancer. Tumour Biol. 37, 14301–14310. doi: 10.1007/s13277-016-5267-8
Porta, C., Sabbatini, R., Procopio, G., Paglino, C., Galligioni, E., and Ortega, C. (2012). Primary resistance to tyrosine kinase inhibitors in patients with advanced renal cell carcinoma: state-of-the-science. Expert. Rev. Anticancer. Ther. 12, 1571–1577. doi: 10.1586/era.12.81
Prior, C., Perez-Gracia, J. L., Garcia-Donas, J., Rodriguez-Antona, C., Guruceaga, E., Esteban, E., et al. (2014). Identification of tissue microRNAs predictive of sunitinib activity in patients with metastatic renal cell carcinoma. PLoS One 9:e86263. doi: 10.1371/journal.pone.0086263
Puente, J., Laínez, N., Dueñas, M., Méndez-Vidal, M. J., Esteban, E., Castellano, D., et al. (2017). Novel potential predictive markers of sunitinib outcomes in long-term responders versus primary refractory patients with metastatic clear-cell renal cell carcinoma. Oncotarget 8, 30410–30421. doi: 10.18632/oncotarget.16494
Qu, L., Ding, J., Chen, C., Wu, Z. J., Liu, B., Gao, Y., et al. (2016). Exosome-transmitted lncARSR promotes sunitinib resistance in renal cancer by acting as a competing endogenous RNA. Cancer Cell 29, 653–668. doi: 10.1016/j.ccell.2016.03.004
Rahbarghazi, R., Jabbari, N., Sani, N. A., Asghari, R., Salimi, L., Kalashani, S. A., et al. (2019). Tumor-derived extracellular vesicles: reliable tools for cancer diagnosis and clinical applications. Cell Commun. Signal 17:73. doi: 10.1186/s12964-019-0390-y
Ralla, B., Busch, J., Flörcken, A., Westermann, J., Zhao, Z., Kilic, E., et al. (2018). miR-9-5p in nephrectomy specimens is a potential predictor of primary resistance to first-line treatment with tyrosine kinase inhibitors in patients with metastatic renal cell carcinoma. Cancer 10:321. doi: 10.3390/cancers10090321
Rasmussen, K. D., and Helin, K. (2016). Role of TET enzymes in DNA methylation, development, and cancer. Genes Dev. 30, 733–750. doi: 10.1101/gad.276568.115
Raval, R. R., Lau, K. W., Tran, M. G., Sowter, H. M., Mandriota, S. J., Li, J. L., et al. (2005). Contrasting properties of hypoxia-inducible factor 1 (HIF-1) and HIF-2 in von Hippel-Lindau-associated renal cell carcinoma. Mol. Cell. Biol. 25, 5675–5686. doi: 10.1128/mcb.25.13.5675-5686.2005
Reustle, A., Fisel, P., Renner, O., Büttner, F., Winter, S., Rausch, S., et al. (2018). Characterization of the breast cancer resistance protein (BCRP/ABCG2) in clear cell renal cell carcinoma. Int. J. Cancer 143, 3181–3193. doi: 10.1002/ijc.31741
Robert, M. F., Morin, S., Beaulieu, N., Gauthier, F., Chute, I. C., Barsalou, A., et al. (2003). DNMT1 is required to maintain CpG methylation and aberrant gene silencing in human cancer cells. Nat. Genet. 33, 61–65. doi: 10.1038/ng1068
Sakai, I., Miyake, H., and Fujisawa, M. (2013). Acquired resistance to sunitinib in human renal cell carcinoma cells is mediated by constitutive activation of signal transduction pathways associated with tumour cell proliferation. BJU Int. 112, E211–E220. doi: 10.1111/j.1464-410X.2012.11655.x
Salmena, L., Poliseno, L., Tay, Y., Kats, L., and Pandolfi, P. P. (2011). A ceRNA hypothesis: the rosetta stone of a hidden RNA language? Cell 146, 353–358. doi: 10.1016/j.cell.2011.07.014
Samuels, Y., Wang, Z., Bardelli, A., Silliman, N., Ptak, J., Szabo, S., et al. (2004). High frequency of mutations of the PIK3CA gene in human cancers. Science 304:554. doi: 10.1126/science.1096502
Santos-Rosa, H., and Caldas, C. (2005). Chromatin modifier enzymes, the histone code and cancer. Eur. J. Cancer 41, 2381–2402. doi: 10.1016/j.ejca.2005.08.010
Sekino, Y., Sakamoto, N., Sentani, K., Oue, N., Teishima, J., Matsubara, A., et al. (2019). miR-130b promotes sunitinib resistance through regulation of PTEN in renal cell carcinoma. Oncology 97, 164–172. doi: 10.1159/000500605
Siegel, R. L., Miller, K. D., and Jemal, A. (2019). Cancer statistics, 2019. CA Cancer J. Clin. 69, 7–34. doi: 10.3322/caac.21551
Stewart, G. D., Powles, T., Van Neste, C., Meynert, A., O’Mahony, F., Laird, A., et al. (2016). Dynamic epigenetic changes to VHL occur with sunitinib in metastatic clear cell renal cancer. Oncotarget 7, 25241–25250. doi: 10.18632/oncotarget.8308
Thiagalingam, S. (2020). Epigenetic memory in development and disease: unraveling the mechanism. Biochim. Biophys. Acta Rev. Cancer 1873:188349. doi: 10.1016/j.bbcan.2020.188349
Valenta, T., Hausmann, G., and Basler, K. (2012). The many faces and functions of β-catenin. EMBO J. 31, 2714–2736. doi: 10.1038/emboj.2012.150
Verbiest, A., Couchy, G., Job, S., Zucman-Rossi, J., Caruana, L., Lerut, E., et al. (2018). Molecular subtypes of clear cell renal cell carcinoma are associated with outcome during pazopanib therapy in the metastatic setting. Clin. Genitourin. Cancer 16, e605–e612. doi: 10.1016/j.clgc.2017.10.017
Vijayan, D. K., and Zhang, K. Y. J. (2019). Human glutaminyl cyclase: structure, function, inhibitors and involvement in Alzheimer’s disease. Pharmacol. Res. 147:104342. doi: 10.1016/j.phrs.2019.104342
Wagener, N., Holland, D., Bulkescher, J., Crnkovic-Mertens, I., Hoppe-Seyler, K., Zentgraf, H., et al. (2008). The enhancer of zeste homolog 2 gene contributes to cell proliferation and apoptosis resistance in renal cell carcinoma cells. Int. J. Cancer 123, 1545–1550. doi: 10.1002/ijc.23683
Wang, M., Sun, Y., Xu, J., Lu, J., Wang, K., Yang, D. R., et al. (2018). Preclinical studies using miR-32-5p to suppress clear cell renal cell carcinoma metastasis via altering the miR-32-5p/TR4/HGF/met signaling. Int. J. Cancer 143, 100–112. doi: 10.1002/ijc.31289
Wang, Z., Zang, C., Cui, K., Schones, D. E., Barski, A., Peng, W., et al. (2009). Genome-wide mapping of HATs and HDACs reveals distinct functions in active and inactive genes. Cell 138, 1019–1031. doi: 10.1016/j.cell.2009.06.049
Wendt, M. K., Allington, T. M., and Schiemann, W. P. (2009). Mechanisms of the epithelial-mesenchymal transition by TGF-beta. Future Oncol. 5, 1145–1168. doi: 10.2217/fon.09.90
Wu, H., and Zhang, Y. (2014). Reversing DNA methylation: mechanisms, genomics, and biological functions. Cell 156, 45–68. doi: 10.1016/j.cell.2013.12.019
Xu, Z., Yang, F., Wei, D., Liu, B., Chen, C., Bao, Y., et al. (2017). Long noncoding RNA-SRLR elicits intrinsic sorafenib resistance via evoking IL-6/STAT3 axis in renal cell carcinoma. Oncogene 36, 1965–1977. doi: 10.1038/onc.2016.356
Yan, L., Ding, B., Liu, H., Zhang, Y., Zeng, J., Hu, J., et al. (2019). Inhibition of SMYD2 suppresses tumor progression by down-regulating microRNA-125b and attenuates multi-drug resistance in renal cell carcinoma. Theranostics 9, 8377–8391. doi: 10.7150/thno.37628
Yang, A. S., Doshi, K. D., Choi, S. W., Mason, J. B., Mannari, R. K., Gharybian, V., et al. (2006). DNA methylation changes after 5-aza-2'-deoxycytidine therapy in patients with leukemia. Cancer Res. 66, 5495–5503. doi: 10.1158/0008-5472.can-05-2385
Zhai, W., Sun, Y., Guo, C., Hu, G., Wang, M., Zheng, J., et al. (2017). LncRNA-SARCC suppresses renal cell carcinoma (RCC) progression via altering the androgen receptor(AR)/miRNA-143-3p signals. Cell Death Differ. 24, 1502–1517. doi: 10.1038/cdd.2017.74
Zhai, W., Sun, Y., Jiang, M., Wang, M., Gasiewicz, T. A., Zheng, J., et al. (2016). Differential regulation of LncRNA-SARCC suppresses VHL-mutant RCC cell proliferation yet promotes VHL-normal RCC cell proliferation via modulating androgen receptor/HIF-2α/C-MYC axis under hypoxia. Oncogene 35, 4866–4880. doi: 10.1038/onc.2016.19
Zhang, P. F., Li, K. S., Shen, Y. H., Gao, P. T., Dong, Z. R., Cai, J. B., et al. (2016). Galectin-1 induces hepatocellular carcinoma EMT and sorafenib resistance by activating FAK/PI3K/AKT signaling. Cell Death Dis. 7:e2201. doi: 10.1038/cddis.2015.324
Zhao, T., Bao, Y., Gan, X., Wang, J., Chen, Q., Dai, Z., et al. (2019). DNA methylation-regulated QPCT promotes sunitinib resistance by increasing HRAS stability in renal cell carcinoma. Theranostics 9, 6175–6190. doi: 10.7150/thno.35572
Zhao, J., and Guan, J. L. (2009). Signal transduction by focal adhesion kinase in cancer. Cancer Metastasis Rev. 28, 35–49. doi: 10.1007/s10555-008-9165-4
Zhou, Q., Guo, X., and Choksi, R. (2017). Activation of focal adhesion kinase and Src mediates acquired sorafenib resistance in A549 human lung adenocarcinoma xenografts. J. Pharmacol. Exp. Ther. 363, 428–443. doi: 10.1124/jpet.117.240507
Zhou, L., Liu, X. D., Sun, M., Zhang, X., German, P., Bai, S., et al. (2016). Targeting MET and AXL overcomes resistance to sunitinib therapy in renal cell carcinoma. Oncogene 35, 2687–2697. doi: 10.1038/onc.2015.343
Keywords: renal cell carcinoma, epigenetics, microRNA, long non-coding RNA, methylation, histone modification, target therapy, tyrosine kinase inhibitor
Citation: Li Q, Zhang Z, Fan Y and Zhang Q (2021) Epigenetic Alterations in Renal Cell Cancer With TKIs Resistance: From Mechanisms to Clinical Applications. Front. Genet. 11:562868. doi: 10.3389/fgene.2020.562868
Edited by:
Xiao Zhu, Guangdong Medical University, ChinaReviewed by:
David Jim Martino, University of Western Australia, AustraliaDuo Xu, Cornell University, United States
Copyright © 2021 Li, Zhang, Fan and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qian Zhang, emhhbmdxaWFuYmptdUAxMjYuY29t; Yu Fan, ZGFudGVmYW5ibXVAMTI2LmNvbQ==
†These authors have contributed equally to this work