Association between miRNA signatures in serum samples from epidermal growth factor inhibitor treated patients and skin toxicity
Metrics: PDF 184 views | Full Text 718 views | ?
Sarah Kemski1,2, Vivien Molitor1, Michael Steffens1, Tim J. Nümm2, Nadine Herrmann2, Thorsten Hornung2, Thomas Bieber2, Christian Schumann3, Volker Kächele4, Thomas Seufferlein5, Volker Heinemann6, Catharina Scholl1 and Julia Carolin Stingl7
1 Research Division, Federal Institute for Drugs and Medical Devices (BfArM), Bonn, Germany
2 Department of Dermatology and Allergy, Christine Kühne Center for Allergy Research and Education (CK-CARE), University Hospital-Bonn, Bonn, Germany
3 Department of Pulmonology, Thorax Oncology, Sleep and Respiration Medicine, Hospital Group Allgäu, Kempten, Germany
4 Medical Centre for Haematology and Oncology, Ulm, Germany
5 Department of Internal Medicine I, University of Ulm, Ulm, Germany
6 Department of Internal Medicine III, Ludwig-Maximilians-University of Munich, Munich, Germany
7 Institute of Clinical Pharmacology, University Hospital of the RWTH Aachen, Aachen, Germany
|Julia Carolin Stingl,||email:||firstname.lastname@example.org|
Keywords: miRNA; EGFR; EGFRI; cancer; skin toxicity
Received: November 09, 2020 Accepted: April 19, 2021 Published: May 11, 2021
Objective: Epidermal growth factor receptor inhibitors (EGFRI) are used as targeted cancer therapy. On average 70% of patients treated with EGFRIs suffer from skin toxicity. Studies showed a correlation between overall survival and the appearance of a skin rash, which is used as a biomarker for therapy efficacy. Micro RNAs (miRNA) as tumor or resistance biomarkers for cancer therapy are also highly investigated. In our study, we searched for associations of miRNA expression profiles in serum, with the severity of skin rash, in order to identify tentative therapy predictive biomarkers.
Materials and Methods: Five candidate miRNAs were selected, based on an earlier in vitro next-generation-sequencing-experiment and after literature search. MiR-21, miR-31, miR-17, miR-106b and miR-520e were investigated in serum samples from patients (n = 254) treated with EGFRI. The quantitative expression of miRNA was tested for association with the occurrence/severity of the rash.
Results: In our cohort of patients treated with EGFR inhibiting monoclonal antibodies, miR-21 and miR-520e serum concentrations were negatively correlated with severity of skin rash (p-value 0.000582 and 1.53e-07 linear-trend-test) whereas for miR-31, a positive correlation was observed (p-value 9.01e-06 linear-trend-test).
Conclusions: This suggests that miR-21, miR-31 and miR-520e expression might be a treatment dependent marker for EGFRI induced skin rash.
Epidermal growth factor receptor inhibitors (EGFRI) belong to the so called targeted cancer therapy and are used for the treatment of different cancer types like non-small-lung-cancer (NSLC), head-and-neck-cancer (head-and-neck-ca), colon-rectal-cancer (colon-ca) and pancreas-cancer (pancreatic -ca) [1–4]. In most of these cancer types the EGF receptor (EGFR) is over expressed but a basal physiological expression of the receptor can also be found in epithelial cells such as keratinocytes . EGFRI can be divided into small molecule tyrosine kinase inhibitors (e.g., erlotinib and gefitinib), which bind to the intracellular tyrosine domain of the receptor and inhibit auto phosphorylation, and monoclonal antibodies (e.g., cetuximab and panitumumab), which block the ligand binding site of the receptor [1, 2, 6, 7]. About 70–80% of patients treated with an EGFRI develop a skin reaction within the second to fourth week of the treatment. The phenotype of EGFRI-induced skin reaction is characterized by inflammatory papulo-pustular follicular rash, skin xerosis, and pruritus occurring on the scalp, face, upper chest and back [8, 9]. Earlier studies showed a strong correlation between the appearance of the skin rash and improved tumor progression and overall survival (OS) of patients treated with EGFRI [10, 11]. Therefore, skin rash is currently used as a predictive biomarker for the efficacy of EGFRI therapy [11, 12]. However, severe skin reactions, especially grade 3 or 4, can cause a break or even the end of therapy. In order to reduce this burden different treatments to reduce skin rash are used in clinical practice. For example the Germane Cancer Association published a guideline for supportive care in 2016 including preemptive and reactive treatment. For a preemptive treatment, they recommend general behavior actions like sun protection and general skin care and an oral prophylaxis with tetracyclines. The recommended reactive treatment depends on the grade of skin rash. For grad 1 skin rash same measures as preemptive treatment plus topical antibiotics, for grade 2 same measures as grade 1 plus topic steroids; for grade 3 and 4 same measures as grad 2 plus systemic steroids or oral isotretionin (but not in combination with antibiotics) is recommended . The treatment recommendations of the guideline are based on different evince levels. The highest evidence can be found for a preemptive treatment of EGFRI induced skin rash with oral tetracycline. For this treatment a double-blind placebo-controlled randomized trial and different randomized trials showed a significant reduction of the severity of skin toxicity in patients treated with tetracycline’s [14–17]. The treatment of already existing skin rash with topic corticosteroids, which have an anti-inflammatory and anti-pruritic effects, shows lower evidence . Also with a low evidence, due to a missing control cohort in the study, orally administrated isotretinoin and clindamycin could reduce grade 2–3 skin rash to grade 0–1 . But with all these possible treatments of EGFRI induced skin rash, its treatment predictive value might get lost and new biomarkers will be needed.
EGFR ligands like epiregulin (EREG), amphiregulin (AREG), and hepatocyte growth factor (HGF) were investigated as biomarkers for skin rash and found to be inversely proportional to grades of skin toxicity . Also genetic variations were assessed as possible biomarkers but only a mutation in EGFR intron 1 showed significant associations toward skin toxicity .
The multi centric Dermatoxgen study was initiated by our research group with the aim to identify tentative biomarkers for prediction of the occurrence of EGFRI induced skin rash. Within the scope of this study, different aspects of the EGFRI-induced skin rash were investigated, all with the objective to better understand EGFRI induced skin rash and find predictive biomarkers. To avoid any bias in the appearance or grade of skin rash, patients in our study only received reactive treatment if necessary with topical corticosteroids, topical antibiotics, oral antibiotics and antihistamines. The study showed a correlation between the occurrence of the skin rash and the drug-metabolizing activity assessed by the erlotinib/O-desmethyl-erlotinib metabolic ratio . An investigation of genetic variations showed that patients carrying the HLA-A*02:01 or HLA-A*03:01 alleles are less likely to develop a skin rash . Also for the EGF receptor associations between skin toxicity of EGFRIs and mutations of the EGFR were found. For the SNP rs534124757 of the EGFR gene it was shown that in patients carrying at least one A allel skin toxicity was less frequently . Deep sequencing of the EGFR gene locus and its downstream signaling pathway revealed also an association between a SNP rs2293348 in the EGFR gene and the rash . But not only mutations in the EGFR itself seem to be associated with skin toxicity. A genotyping of cancer patients revealed an association between skin rash and a specific haplotype of phosphatidylinositol-3-kinase (PIK3CA) . An investigation of biomarkers in the plasma of the Dermatogen study group showed that patients with a lower IL-8 concentration are in risk to develop a more severe skin rash but had a longer survival. The same was found for patients with a low HGF concentration in the plasma [26, 27].
Gene regulatory biomarkers of EGFRI-induced skin rash
MicroRNAs (miRNAs) are short non-coding RNAs with a length of around 19–25 nucleotides, that are involved in the post transcriptional regulation of gene expression. One miRNA can regulate the expression of many different genes by inducing translational inhibition or transcript degradation. Therefore, miRNAs are important for cell processes like proliferation, differentiation, apoptosis, stress tolerance and immune response . The role and machinery of miRNAs in skin cells have been broadly investigated . MiR-21 for example is involved in migration processes of keratinocytes and fibroblasts . Some miRNAs are up- or down-regulated in pathophysiological processes. In inflammatory skin diseases like psoriasis and atopic dermatitis, miRNAs are suggested to be involved in inflammatory processes and immune dysfunction [30, 31]. However, miRNAs are not only highly involved in skin physiology and in pathophysiology; they are also in the focus of research in the field of tumor diagnostic and prognostic biomarkers and the resistance to cancer therapy [32, 33]. A set of 5 miRNAs (mir-17, mir-660, mir-92a, mir-106a, and mir-19b) is most commonly dysregulated in serum and tumor tissue samples from lung cancer patients and has therefore promising diagnostic potential . In a meta-analysis circulating miR-31 was determined as an effective biomarker for cancer detection and prognosis . MiRNAs may also be markers for therapy resistance. Different miRNAs mediate drug resistance in colon cancer by distinct mechanisms . 16 miRNAs were differentially expressed between EGFRI resistant patients and EGFRI sensitive patients and might therefor be dependent for the sensitivity to EGFRI treatment . Cetuximab resistance of colon cancer patients is connected to an upregulation of miR-199a-5p and miR-375 targeting PH domain and leucine-rich repeat protein phosphatase 1 (PHLPP1) a tumor suppressor .
All this shows the involvement of miRNAs in physiological and pathophysiological processes in skin cells, especially in inflammatory diseases like psoriasis or atopic dermatitis, and their already frequently discussed role as biomarkers in cancer. This leads to the assumption that miRNAs might also be involved in the development of EGFRI-induced skin rash and therefor can be used as possible biomarkers for therapeutic resistance. However, so far most of the studies in cancer patients focused on miRNAs expressed in cancer cells or secreted into body fluids most likely by cancer cells, while miRNAs associated with the skin rash have not been investigated. To elucidate the influence of miRNAs in the development of skin rash we previously investigated the miRNA expression in fibroblasts incubated with and without erlotinib . In this analysis now, we firstly identify the miRNA expression profiles in skin cells, namely keratinocytes and fibroblasts, depending on their capability to react to erlotinib incubation. From those we depicted a set of five miRNAs which concentrations were determined in serum samples, derived from EGFRI treated patients, and evaluated with respect to the appearance and severity of the skin rash.
Erlotinib induces a distinct miRNA profile in primary human keratinocytes and fibroblasts
In order to gain a first impression about the effect of EGFRI incubation on the miRNA expression in skin cells, we studied effects of in vitro EGFR inhibition with the (small molecule) EGFRI erlotinib on miRNA expression in human keratinocyte/and fibroblast cell samples from healthy donors. Cells were classified into cell samples more reactive on erlotinib incubation and less reactive ones. Detailed information on the classification is provided in the supplements (Supplementary Figures 1 and 2). Fifty-four miRNAs were exclusively up or down regulated in the cell samples that were more reactive on erlotinib incubation compared to less reactive ones (fold change ≥ 1.5 and ≤ 0.66; 32 miRNAs in keratinocytes with p-value ≤ 0.05; 22 miRNAs in fibroblasts with p-value < 0.01) (Figure 1; Supplementary Tables 1 and 2). In order to decide, which of those 54 differentially expressed miRNAs were the most promising ones to be further investigated, a literature search on miRNAs involved in regulation of EGFR pathways was conducted. The search included all publications in which the miRNAs were connected to EGFR inhibitor treatment, erlotinib treatment in particular and cancer cells or bio fluid samples from cancer patients. According to our cell culture findings and the literature reports, five miRNAs (miR-31, miR-17, miR-106b, miR-520e, miR-21) were chosen as candidates to be further investigated in serum samples from EGFRI treated patients for potential correlations with the development of a skin rash.
Figure 1: Erlotinib-sensitive keratinocytes and fibroblasts express a different miRNA pattern than those insensitive to erlotinib. Cells were treated with 5μM erlotinib or 0.05% DMSO and stimulated with 4 nM EGF. (A) Differentially expressed miRNAs between erlotinib sensitive and insensitive primary human keratinocyte cell samples (fold change ≥ 1.5 or fold change ≤ 0.66, p-value ≤ 0.05). (B) Differentially expressed miRNAs between erlotinib sensitive and insensitive human dermal fibroblast cell samples (fold change ≥ 1.5 or fold change ≤ 0.66, p-value < 0.01) Abbreviation: EGF: epidermal growth factor.
MiRNA serum levels of EGFRI treated patients are associated to EGFRI induced skin rash
Details about the patients who provided serum samples for the analysis of genetic and epigenetic biomarkers for EGFRI related skin rash from the Dermatoxgen study, are given in Tables 1 and 2. Patients were observed on week 1, 2, 3, 4, 26 and 53 after the start of the therapy. In week 4 serum, plasma and blood samples were taken. Serum samples were available from 254 patients and total RNA including miRNA was extracted.
Table 1: Patients characteristics (age, gender and smoking) and clinical characteristics (tumor type and skin rash)
|Patient characteristics||Patients||TKI treated patients||EGFR-mAbs treated patients|
|Sample size||N (254)||N (156)||N (98)|
|Age: Mean (SD)||66.4 (9.8)||67.3 (9.0)||65.0 (10.8)|
|Gender: N (%)|
|Male||161 (63.4)||92 (59.0)||69 (70.4)|
|Female||93 (36.6)||67 (41.0)||29 (29.6)|
|Smoking: N (%)|
|No||94 (38.1)||52 (34.3)||42 (43.8)|
|Yes (present)||30 (12.1)||15 (9.9)||15 (15.6)|
|Yes (former)||123 (49.8)||84 (55.6)||39 (40.6)|
|Lung-Cancer||136 (53.5)||107 (68.6)||29 (29.6)|
|Colon-Cancer||58 (22.8)||2 (1.3)||56 (57.1)|
|Head and Neck Cancer||11 (4.3)||0||11 (11.2)|
|Pancreatic Cancer||49 (19.3)||47 (30.1)||2 (2.0)|
|No||51 (20.1)||35 (22.4)||16 (16.3)|
|Grade 1||98 (38.6)||60 (38.5)||38 (38.8)|
|Grade 2||92 (36.2)||52 (33.3)||40 (40.8)|
|Grade 3||13 (5.1)||9 (5.8)||4 (4.1)|
Table 2: Distribution of therapeutic agents (EGFRIs) and cancer types in the patient cohort
|Head and Neck Cancer||11||0||0||0||11|
All five candidate miRNAs (miR-31, miR-17, miR-106b, miR-520e, miR-21), chosen from the in vitro NGS experiments and from the literature, were quantified in serum samples of EGFRI treated patients and correlated to the clinically observed severity and course of skin rash.
All four included EGFRIs have a skin rash as possible side effect and the skin rash appears tumor type and state independently. Therefor the association between the miRNA concentration and skin rash was analyzed independently from the tumor type and treatment, in all 254 patients.
Of the five miRNAs studied, miR-31, miR-21 and miR-520e showed correlations (p-value 0.00296, 0.0378 and 0.0199 [linear trend test], p-values not corrected for multiple testing) between serum concentration and the severity of the skin rash (Table 3).
Table 3: Associations between the serum concentration of the six tested miRNAs and the occurrence/severity of the skin rash
|miRNA||All Patients (N 254)||EGFR-mAbs treated patients (N 98)||TKI treaded patients (N 156)|
|miR-17||Appearance of skin rash||p = 0.376||p = 0.264||p = .429|
|Severity of skin rash||p = 0.594||p = 0.613||p = .62|
|miR-21||Appearance of skin rash||p = 0.098||p = 0.047||p = .439|
|Severity of skin rash||p = 0.002||p < 0.001||p = .115|
|miR-31||Appearance of skin rash||p = 0.564||p = 0.081||p = .472|
|Severity of skin rash||p = 0.037||p < 0.001||p = .43|
|miR-106b||Appearance of skin rash||p = 0.314||p = 0.938||p = .256|
|Severity of skin rash||p = 0.182||p = 0.699||p = 0.168|
|miR-520e||Appearance of skin rash||p = 0.126||p = 0.044||p = 0.651|
|Severity of skin rash||p = 0.019||p < 0.001||p = 0.473|
MiR-21 and miR-520 were negatively correlated with skin rash and showed the lowest expression in patients with severe grade 3 skin rash.
MiR-31 was positively correlated with skin rash with the highest expression in patients with grade 3 skin rash. For miR-17 and miR-106b, no significant association with severity of skin rash was found (Table 3).
Low miR-21, miR-520e and miR-31 serum concentrations are associated with longer survival
For those three miRNAs, which were correlated to the severity of the skin rash, a survival curve analysis was performed. The Kaplan-Meier-plots showed a significant correlation between a low miRNA serum level and a longer overall survival for miR-21 and miR-520e (p-value 0.00046 and 0.0088) (Figure 2A and 2C). For miR-31 a trend could be seen between low miR-31 serum level and longer survival (p-value 0.055) (Figure 2B).
Figure 2: A lower miRNA concentration is associated with a longer survival for miR-21, miR-31 and miR-520e. Patients were followed-up for 360 days after initiation of EGFRI therapy. The proportion of patients still alive is plotted over the observation period. Patients were split by delta CP values of each miRNA into four equally-sized groups. (A) Kaplan-Meier- Plot for miR-21; log rank test, p-value < 0.001; mean OS for miR-21 concentration 1: 221 days (SE: 15.3); mean OS for miR-21 concentration 2: 290 days (SE: 13.3); mean OS for miR-21 concentration 3: 271 days (SE: 14.1); mean OS for miR-21 concentration 4: 275 days (SE: 15.0). (B) Kaplan-Meier- Plot for miR-31; log rank test, p-value = 0.055; mean OS for miR-31 concentration 1: 233 days (SE: 14.8); mean OS for miR-31 concentration 2: 272 days (SE: 14.9) ; mean OS for miR-31 concentration 3: 275 days (SE: 15.0); mean OS for miR-31 concentration 4: 275 days (SE: 14.0). (C) Kaplan-Meier- Plot for miR-21; log rank test, p-value = 0.0088; mean OS for miR-520e concentration 1: 247 days (SE: 15.0); mean OS for miR-520e concentration 2: 261 days (SE: 14.8); mean OS for miR-520e concentration 3: 245 days (SE: 15.7); mean OS for miR-520e concentration 4: 302 days (SE: 12.7) Abbreviations: #: number of patients; dCP: delta crossing point; OS: overall survival.
Association between miR-21, miR-31 and miR-520e serum concentrations and skin rash is not treatment independent
Subgroup analyses of patients who had been treated with either small molecule or antibody tyrosine kinase inhibitors, revealed significant correlations between miR-21, miR-31, and miR-520e and the severity of the skin rash only in the antibody treatment group (p-values of 0.000582, 9.01e-06 and 1.53e-07 [linear trend test], not corrected for multiple testing) (Figure 3A–3C). As for the whole cohort a negative correlation of miR-21 and miR-520e serum concentration with the severity of the skin rash could be seen. For miR-31 a positive correlation with the severity of the skin rash was found. However, no significant correlation between miRNA serum concentration and skin rash in patients treated with small molecule tyrosine kinase inhibitors was observed.
Figure 3: Significant correlation between miRNA serum concentration and severity of the skin rash for patients treated with monoclonal antibody EGFRIs. MiRNA concentrations were determined by qPCR (n = 98). dCP values were calculated against miR-93. A high dCP value means the miRNA is down regulated. (A) miR-21 concentration plotted against the maximum severity of the skin rash during observation period; linear trend test, p-value < 0.001. (B) miR-31 concentration plotted against the maximum severity of the skin rash during observation period; linear trend test, p-value < 0.001. (C) miR-520e concentration plotted against the maximum severity of the skin rash during observation period; linear trend test, p-value < 0.001. (0 = no skin rash, 1 = light skin rash, 2 = mild skin rash, 3 = sever skin rash). Abbreviation: dCP: delta crossing point.
In our study, we analyzed the association between the typical skin rash developed under EGFRI therapy and the serum concentration of different miRNAs. Understanding this connection might lead to a better understanding of epigenetic mechanisms behind the skin rash and help to find predictive biomarkers for its development. Out of the five candidate miRNAs, chosen from a combination of literature search and prior in vitro experiments in a skin cell model, miR-21, miR-31 and miR-520e showed significant correlations between skin rash and serum concentration.
So far, there were many different studies searching for miRNAs as possible biomarkers in cancer development and disease outcome using tissue samples from the tumor, the surrounding tissue and/or blood samples [39, 40]. The objective of this study was to identify epigenetic miRNA associations with therapy-induced skin rash in patients receiving an EGFR targeting tumor treatment, since skin toxicity of EGFRI treatment has been shown to be a prognostic and therapy-predictive factor. In order to identify miRNAs that are involved in EGFR inhibition, we started by analyzing the effect of erlotinib incubation on the miRNA expression in primary human keratinocytes and fibroblasts from healthy donors. A literature search about which of the significantly regulated miRNAs from the in vitro experiments, were already connected to cancer and/or EGFRI therapy, gave us an idea, which miRNAs might be involved in the appearance of the skin rash and could potentially be found in serum samples from EGFRI treated patients.
Linear trend tests revealed negative correlations between the serum concentrations of miR-21 and miR-520e and the severity of the skin rash. Survival curve analysis also showed a negative correlation between miR-21 and miR-520e serum concentrations and overall survival time. For miR-31 a positive correlation between serum concentration and severity of the skin was observed. The survival curve for miR-31 showed a negative correlation between miR-31 serum concentration and survival time.
Because of the frequent mentioning of miR-21 in publications on miRNAs and cancer treatment resistance, this miRNA was the only one we investigated in our patient cohort even though it was not significantly differentially regulated in the in vitro cell experiments. Hua Shen et al. discovered that a high miR-21 and low phosphatase and tensin homolog (PTEN) concentration might indicate a poor gefitinib clinical response . Another study investigated the miR-21 expression in colon-cancer tissue compared to normal adjacent tumor tissues. The colon-cancer tissue had higher miR-21 content and tumor patients with a lower miR-21 expression had a longer overall survival (OS) and disease-free survival (DFS) [41, 42]. These results in tumor tissue fit well to our finding that lower miR-21 concentration correlated with higher grade of skin rash, since skin rash has a positive prognostic value. In our cohort, we also observed a correlation between miR-21 in serum and a longer survival of patients. However, in our study we looked at treatment related outcome of the patients, and therefor included different tumor and disease stages. The Kaplan-Meier plot demonstrates that a low miR-21 serum concentration was correlated to a longer survival and miR-21 serum concentration was inversely correlated with the grade of skin rash. Our results suggest a better response to the EGFRI treatment in patients with a low miR-21 serum concentration.
MiR-520 is mostly studied in patients with non-small-lung-cancer. In NSLC metastasis miR-520 seems to be upregulated and by this down regulates transforming growth factor-β receptor 2 (TGFBR2) which regulates carcinogenesis and metastasis via the transforming growth factor-β (TGF-β) pathway . In another study mmu-miR-291a-3p, a mouse homolog of miR-520e, inhibits cell senescence via the TGF-β receptor 2 signaling pathway by reducing protein expression of TGF-β receptor 2. Mmu-miR-291a-3p also seems to enhance the proliferative potential of senescent cells. . In our results, patients with EGFR inhibitor induced skin rash had lower expression of miR-520e indicating a better prognosis, and this was especially significant in patients treated with monoclonal antibodies. In our study, we used serum samples instead of NSCLC tissue for miRNA isolation, but in principle, the better prognosis was correlated with lower serum concentrations, which fits well to the findings in lung tissue. In addition, the fact that the correlation between lower miR-520e serum concentrations and skin rash gets even more significant looking at patients treated with monoclonal antibodies, shows that low miR-520e in our study might be an indicator for good response to the therapy because the miR-520e we measured might not originate from tumor cells. We could also see in the Kaplan- Meier plots that patients with low miR-520e expression have a longer survival time. A possible link between the grade of skin rash and lower miR-520e expression might be TGF-β. Patients with atopic dermatitis seem to have a higher expression of TGF-β and other inflammatory cytokines like IL-6, IL-8 and IL-10 compared to healthy individuals .
For miR-31 a study by Ning Xu et al. showed an upregulation in keratinocytes from psoriasis patients. They demonstrated that TGF-β1, which is also upregulated in psoriasis epidermis, upregulates miR-31. Furthermore, an inhibition of miR-31 suppressed NF-kB–driven promoter luciferase activity and the basal and tumor necrosis factor-α (TNF-α) induced production of IL-1b, CXCL1/growth-related oncogene-a, CXCL5/epithelial-derived neutrophil-activating peptide 78, and CXCL8/IL-8 in human primary keratinocytes . Regarding pathophysiology of the EGFRI induced skin rash, an inflammatory component is postulated by several publications. [8, 47, 48] This together with the effect of miR-31 on inflammatory cytokines is in line with our finding that a higher miR-31 concentration is associated with a higher grade of skin rash. A similar relationship can also be observed for miR-520e. On the other hand, we found a correlation between low miR-31 serum concentration and a longer survival. In a meta-analysis of 14 different studies including different tumor types, a significant correlation between low miR-31 concentration in blood and a longer overall survival was shown . Studies in tumor biopsies confirm that a low miR-31 expression correlates with a longer overall survival . Even for EGFR- inhibitor therapy a study comparing cetuximab treated patients with low or high miR-31 concentration in tumor tissue found that a low miR-31 expression leads to longer OS and PFS . The differences in correlation between miR-31 serum concentration and skin rash or OS might lead to the assumption, that miR-31 could have a different effect on skin cells than tumor cells under EGFRI treatment. Both, the study by Anandappa et al. and the one by Laurent-Puig et al. used tumor tissue instead of serum samples. The meta-analysis by Ma et al. used blood samples but they did not take into account the treatment of the cancer patients. Those two differences might support the suggestion that the miR-31 expression we found is treatment dependent and that miR-31 has different effects on skin cells like keratinocytes or fibroblasts than on tumor cells. MiR-31 expression is elevated in keratinocytes incubated with inflammatory cytokines like TNF-α or interferon-γ . Furthermore, an interferon-γ incubation of keratinocytes in which the EGFR was blocked resulted in an increased chemokine expression . This might be an explanation for the correlation between high miR-31 serum concentration and high grade of skin rash reflecting the inflammatory component of the skin rash. Blocking the EGF receptor with monoclonal antibodies might lead to an inflammatory reaction in the skin cells, which in turn leads to an elevated miR-31 expression in skin cells, which then might release miR-31 to the blood.
Looking at our results a correlation between miRNA serum levels and the skin rash under EGFRI treatment can be found. For three out of five candidate miRNAs, a significant association between the miRNA concentration and the severity of the skin rash was shown especially for patients treated with monoclonal antibodies. This suggests a specificity of miR-21, miR-31 and miR-520e for the therapy with monoclonal antibodies and expression of the analyzed miRNAs most likely by cells other than tumor cells, such as skin cells. As a possible mechanism, an inflammatory component could be identified by looking at the miRNAs we analyzed and their connection to other skin diseases. Inflammation is already discussed to be part of the mechanism of EGFRI induced skin toxicity . Taking all this together, it makes it quite interesting to further investigate the role of miR-21, miR-31 and miR-520e for the development of EGFR inhibitor induced skin rash in patients treated with monoclonal antibodies and to understand possible mechanisms behind it. A better understanding of the mechanism might than lead to possible predictive biomarkers for EGFRI induced skin rash, which are needed as an alternative to the skin rash is self as prognostic biomarker for therapy response. Because with a better treatment of the EGFRI induced skin rash it prognostic value might get lost.
The biggest limitation might be that there are no serum samples from patients before they were treated with an EGFR inhibitor. Hence, we cannot say if the effects we found are mostly due to the therapy with an EGFRI or if those effects can be seen in patients in general. However, because of the heterogeneity of the cohort concerning the tumor type and cancer state, a treatment specific effect is possible. Still these miRNAs give us a better understanding about the skin rash, which can be used for further investigations.
Materials and Methods
In vitro miRNA profiling in human primary skin cells
Primary human keratinocytes were isolated from skin samples derived from healthy controls of the Dermatox_Epigen study. This study aims to analyze epigenetic differences in skin samples of EGFR inhibitor treated cancer patients. As a control collective patients from the university women’s clinic Bonn were enrolled. These patients underwent plastic surgeries from which skin samples were derived to isolate keratinocytes. After obtaining the skin, a wash step was performed with a sterile gauze soaked in PBS to remove excess blood. Afterwards, the split thickness-skin (0.4 mm) was immediately prepared using a dermatome. The epidermal cell suspension was generated with freshly prepared trypsinization solution (PBS supplemented with 0.5% trypsin (Sigma-Aldrich) and 1% Antibiotic-Antimycotic (Gibco™, Waltham, MA, USA)). The floating split thickness-skin was incubated in the trypsin solution for 1 h at 37°C, whereby the dermal side was in contact with the trypsin solution. Afterwards, the dermis was removed and the epidermis was added to medium supplemented with 1% DNase to remove DNA and to avoid cell clumping. Next, the epidermal cell suspension was generated by several steps of pipetting up and down or by vortexing until the media became cloudy, indicating that the cells were in suspension. Written informed consent was obtained from all volunteers and the study was approved by the ethical boards of Bonn University. Donors were in average 38.8 (SD 9.6) years old and included 3 male and 9 female.
Primary human dermal fibroblasts were kindly provided be the Clinical Pharmacology Department of the University Medicine Göttingen (Dr. med. Markus Schirmer). Fibroblasts were isolated from healthy tissue removed by dermal excisions from patients at the clinic.
Keratinocytes were cultured in EpiLife medium with 1% human keratinocyte growth supplements (Gibco™, Waltham, MA, USA), 1% pen/strep and 0.1% amphotericin B, medium was changed every 2–3 days and cells were passaged when reaching a confluence of 75%. Fibroblasts were cultured in basal fibroblast growth medium 2 with supplement mix fibroblast growth medium 2 (2% FCS, 5 μg/ml Insulin, 0.001 μg/ml fibroblast growth factor (FGF)) (PromoCell, Heidelberg, Germany) and 1% pen/strep. The medium was changed every 2–3 days and cells were passaged when reaching a confluence of 75% every 5–7 days. Before miRNA profiling keratinocytes were incubated with 5μM erlotinib (SantaCruz Biotechnology, Dallas, TX, USA) or 0.05% DMSO for two hours and then stimulated with 4 nM EGF (PeproTech, Rocky Hill, NJ, USA) for 5 min. Fibroblasts were incubated with 5 μM erlotinib or 0.05% DMSO for 24 hours and then stimulated with 4nM EGF for 5 min.
RNA isolation and miRNA profiling with Next Generation Sequencing
Total RNA from keratinocytes and fibroblasts was extracted using Trizol/ Chloroform. RNA was precipitated from the aqueous phase with isopropanol, resuspended in RNase-free water and stored at –80°C.
Library preparation for miRNA profiling was performed with the NEB Next Multiplex Small RNA Library Prep Set for Illumina (New England Biolabs, Ipswich, MA, USA) following the manufacturer’s protocol. Resulting cDNA library was purified with QiaQuick PCR purification kit (Qiagen, Hilden, Germany). A size selection of the generated cDNA was done with gel electrophoresis. 60 μl cDNA probes were mixed with 12 μl loading dye and 15 μl were pipetted into 5 pockets of an 10% TBE-Gel (Invitrogen by Thermo Fisher, Carlsbad, CA, USA) and separated for 1 h 20 min at 150 V in TBE buffer. Ethidium bromide was used for the visualization of the cDNA bands and the one at ~ 140 bp was cut out. Gel extraction of the cDNA was done following the manufactures protocol of the NEB Next Multiplex Small RNA Library Prep Set for Illumina. For the quantification and qualification of the size-selected cDNA the Agilent high sensitivity DNA Kit (Agilent, Santa Clara, CA, USA)) was used and analyzed with the Agilent Bioanalyzer 2100. MiRNA sequencing was performed with the MiSeq from Illumina using the MiSeq Reagent Kit v3 (Illumina, San Diego, CA, USA) following the manufactures protocol. A final concentration 15 pMol of cDNA was used.
Analysis of sequencing data
Sequencing data was saved as fastq files which was then used for adapter trimming with the Cutadapt 1.9 software . Alignment of the trimmed sequencing data to the human genome (build37) was performed using bowtie1.0 . The same software was used for the aligment to published mature and precursor (hairpin) miRNAs collected in the miRBase data base (release 21). Read counts were normalised and miRNA expression was compared between keratinocytes/fibroblasts more reactive to erlotinib against less reactive cell samples. For further research miRNAs, which were significantly differentially expressed with a fold change ≥ 1.5 or a fold change ≤ 0.66 and a p value of 0.05 for keratinocytes or a p value of 0.01 for fibroblasts, were chosen. A cut of at a p-value of 0.01 was choosen for fibroblasts, because of the high amount of segnificant miRNAs at a cut of at a p-value of 0.05. For a better comparison of the three different treatments, Venn diagrams of differentially expressed miRNAS were generated. Data was uploaded on to the GEO database under the accession number GSE159602.
Literature search to choose miRNAs for further experiments
Fifty-four miRNAs were significantly differentially expressed exclusively after erlotinib incubation. A literature search was used to choose the most promising miRNAs for further experiments in serum samples from EGFRI treated patients. For the search the PubMed database was used. Search terms consisted of the name of one of the 54 miRNAs together with “cancer”, “EGFR inhibitors” or “erlotinib”. We did not discriminate between miRNA found in blood samples or in tumor cells/biopsies and we included all EGFR inhibitor therapies.
Analyzing miRNAs in patients of the dermatoxgen study
The Dermatoxgen study is a prospective, multicentric study with the aim to investigate pharmacogenetics factors of EGFR inhibitor induced skin toxicity. Patients with histologically confirmed solid tumors (pancreatic, colon, head and neck or non-small-cell lung cancer) who received an EGFR inhibitor (erlotinib, gefitinib, cetuximab or panitumumab) therapy for the first time were included in the Dermatoxgen study. Patient characteristics are shown in Tables 1 and 2. The study was approved by the ethical boards of Ulm University and the Ludwig-Maximilians-University of Munich and patients gave their written informed consent to participate. EGFRI admission was carried out according to approved indications, as previously described . The appearance and severity of the skin rash were documented once a week after treatment start for 4 weeks. The severity of the skin rash was graded following the Toxicity Criteria for Adverse Events of the American National Cancer Institute (NCICTCAE version 3.0, 2006). Blood, serum and plasma samples were collected after 4 weeks of treatment. The skin rash was only treated reactively if it was necessary during the treatment. As treatment topical corticosteroids, topical antibiotics, oral antibiotics and antihistamines were used depending on the grade of skin rash. A preemptive treatment was excluded to avoid any bias which may arise from the suppressing effect of the preemptive treatment. Patients were followed-up for 12 months, with visits after 6 months and 12 months. The survival status at 360 days after initiation of EGFRI treatment was used for Kaplan-Meier analyses. Serum samples from 254 patients were used for the miRNA analysis.
Collection of serum samples and miRNA extraction
Serum samples were collected during week 4 of the EGFRI treatment right before the admission of the next scheduled dose. About 7.5 ml blood were collected from each patient using blood-sampling tubes (Serum S-Monovette®, 7.5 ml, Sarstedt). After a rest of 20 min at room temperature blood samples were centrifuged at 2500× g for 10 min, aliquoted and stored immediately at –80°C.
miRNA was isolated with the miRNeasy serum/plasma kit (Qiagen, Hilden, Germany). Following the manufactures protocol 100 μl serum were used. 350 μl of the upper aquarious phase mixed with 525 μl 100% ethanol were pipetted onto the miRNeasy spin column. After several washing steps, RNA was eluted with 14 μl RNase-free water from the spin column.
Isolated miRNA was transcribed into cDNA using the miScript RT kit (Qiagen, Hilden, Germany). 5 μl miRNA sample were processed following the manufacturer’s protocol using the HiFlex Buffer. After reverse transcription the samples were diluted in 80μl RNase-free water. qPCR was performed using the miScript primer assays (miR-21, miR-31, miR-520e, miR-17, miR-106b) from Qiagen (Qiagen, Hilden, Germany). A master mix was prepared consisting of miScript Universal primer, miScript Primer assay, SyberGreen kit and RNase-free water (Qiagen, Hilden, Germany) following the manufacturer’s protocol. Samples were analyzed as triplicates in a 364 well plate with the Light Cycler 480 from Roche (Roche, Basel, Switzerland). As reference, miRNA-93 (Qiagen, Hilden, Germany) was used.
EGFRI: Epidermal growth factor receptor inhibitor; miRNA: micro RNA; NSLC: Non-small-lung-cancer; Head-and-neck-ca: head-and-neck-cancer; colon-ca: colon-rectal-cancer; pancreatic-ca: panctras-cancer; EGFR: epidermal growth factor receptor; ERFG: epiregulin; AREG: amphiregulin; SNP: Single Nucleotide Polymorphism; HGF: hepatocyte growth factor; PIK3CA: phosphatidylinositol-3-kinase; PHLPP1: PH domain and leucine-rich repeat protein phosphatase 1; EGFR-mAbs: epidermal growth factor receptor monoclonal antibodies; TKI: tyrosine kinase inhibitor; N: number of patients; NA: not known; dCP: delta crossing point; OS: overall survival; PTEN: phosphatase and tensin homolog; DFS: disease-free survival; TGFBR2: transforming growth factor-β receptor 2; TGF-β: transforming growth factor-β; TNF-α: tumor necrosis factor-α.
The authors thank the Wilhelm Sander Stiftung for founding the Dermatoxgen Study (DERMATOXGEN Nr. 2008.017). They also thank Dr. med. Markus Schirmer from the Clinical Pharmacology Department of the University Medicine Göttingen for providing the fibroblast cell lines. We thank Dr. Walgenbach of the Department of Plastic and Aesthetic Surgery, University of Bonn, Germany, for providing skin from healthy individuals for the isolation of keratinocytes.
CONFLICTS OF INTEREST
Authors have no conflicts of interest to declare.
The Dermatoxgen Study was founded by the Wilhelm Sander Stiftung (DERMATOXGEN Nr. 2008.017). Dermatox_Epigen study was supported by a collaborative grant from the Federal Institute for Drugs and Medical Devices (BfArM, Bonn, Germany, V-16700/68502/2016–2019). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
3. Moore MJ, Goldstein D, Hamm J, Figer A, Hecht JR, Gallinger S, Au HJ, Murawa P, Walde D, Wolff RA, Campos D, Lim R, Ding K, et al, and National Cancer Institute of Canada Clinical Trials Group. Erlotinib plus gemcitabine compared with gemcitabine alone in patients with advanced pancreatic cancer: a phase III trial of the National Cancer Institute of Canada Clinical Trials Group. J Clin Oncol. 2007; 25:1960–66. https://doi.org/10.1200/JCO.2006.07.9525. [PubMed].
8. Lacouture ME, Anadkat M, Jatoi A, Garawin T, Bohac C, Mitchell E. Dermatologic Toxicity Occurring During Anti-EGFR Monoclonal Inhibitor Therapy in Patients With Metastatic Colorectal Cancer: A Systematic Review. Clin Colorectal Cancer. 2018; 17:85–96. https://doi.org/10.1016/j.clcc.2017.12.004. [PubMed].
10. Perez-Soler R, Van Cutsem E. Clinical research of EGFR inhibitors and related dermatologic toxicities. Oncology (Williston Park). 2007; 21:10–6. [PubMed].
11. Wacker B, Nagrani T, Weinberg J, Witt K, Clark G, Cagnoni PJ. Correlation between development of rash and efficacy in patients treated with the epidermal growth factor receptor tyrosine kinase inhibitor erlotinib in two large phase III studies. Clin Cancer Res. 2007; 13:3913–21. https://doi.org/10.1158/1078-0432.CCR-06-2610. [PubMed].
12. Abdel-Rahman O, Fouad M. Correlation of cetuximab-induced skin rash and outcomes of solid tumor patients treated with cetuximab: a systematic review and meta-analysis. Crit Rev Oncol Hematol. 2015; 93:127–35. https://doi.org/10.1016/j.critrevonc.2014.07.005. [PubMed].
13. Franziska J. S3-Leitlinie Supportive Therapie bei onkologischen PatientInnen. Deutsche Krebsgesellschaft e.V. (DKG)). 2020.
14. Jatoi A, Rowland K, Sloan JA, Gross HM, Fishkin PA, Kahanic SP, Novotny PJ, Schaefer PL, Johnson DB, Tschetter LK, Loprinzi CL. Tetracycline to prevent epidermal growth factor receptor inhibitor-induced skin rashes: results of a placebo-controlled trial from the North Central Cancer Treatment Group (N03CB). Cancer. 2008; 113:847–53. https://doi.org/10.1002/cncr.23621. [PubMed].
15. Lacouture ME, Mitchell EP, Piperdi B, Pillai MV, Shearer H, Iannotti N, Xu F, Yassine M. Skin toxicity evaluation protocol with panitumumab (STEPP), a phase II, open-label, randomized trial evaluating the impact of a pre-Emptive Skin treatment regimen on skin toxicities and quality of life in patients with metastatic colorectal cancer. J Clin Oncol. 2010; 28:1351–57. https://doi.org/10.1200/JCO.2008.21.7828. [PubMed].
16. Melosky B, Anderson H, Burkes RL, Chu Q, Hao D, Ho V, Ho C, Lam W, Lee CW, Leighl NB, Murray N, Sun S, Winston R, Laskin JJ. Pan Canadian Rash Trial: A Randomized Phase III Trial Evaluating the Impact of a Prophylactic Skin Treatment Regimen on Epidermal Growth Factor Receptor-Tyrosine Kinase Inhibitor-Induced Skin Toxicities in Patients With Metastatic Lung Cancer. J Clin Oncol. 2016; 34:810–15. https://doi.org/10.1200/JCO.2015.62.3918. [PubMed].
17. Arrieta O, Vega-González MT, López-Macías D, Martínez-Hernández JN, Bacon-Fonseca L, Macedo-Pérez EO, Ramírez-Tirado LA, Flores-Estrada D, de la Garza-Salazar J. Randomized, open-label trial evaluating the preventive effect of tetracycline on afatinib induced-skin toxicities in non-small cell lung cancer patients. Lung Cancer. 2015; 88:282–88. https://doi.org/10.1016/j.lungcan.2015.03.019. [PubMed].
18. Annunziata MC, De Stefano A, Fabbrocini G, Leo S, Marchetti P, Romano MC, Romano I. Current Recommendations and Novel Strategies for the Management of Skin Toxicities Related to Anti-EGFR Therapies in Patients with Metastatic Colorectal Cancer. Clin Drug Investig. 2019; 39:825–34. https://doi.org/10.1007/s40261-019-00811-7. [PubMed].
19. Bidoli P, Cortinovis DL, Colombo I, Crippa A, Cicchiello F, Villa F, Cazzaniga ME, Altomare G. Isotretinoin plus clindamycin seem highly effective against severe erlotinib-induced skin rash in advanced non-small cell lung cancer. J Thorac Oncol. 2010; 5:1662–63. https://doi.org/10.1097/JTO.0b013e3181ec1729. [PubMed].
20. Takahashi N, Yamada Y, Furuta K, Nagashima K, Kubo A, Sasaki Y, Shoji H, Honma Y, Iwasa S, Okita N, Takashima A, Kato K, Hamaguchi T, Shimada Y. Association between serum ligands and the skin toxicity of anti-epidermal growth factor receptor antibody in metastatic colorectal cancer. Cancer Sci. 2015; 106:604–10. https://doi.org/10.1111/cas.12642. [PubMed].
21. Graziano F, Ruzzo A, Loupakis F, Canestrari E, Santini D, Catalano V, Bisonni R, Torresi U, Floriani I, Schiavon G, Andreoni F, Maltese P, Rulli E, et al. Pharmacogenetic profiling for cetuximab plus irinotecan therapy in patients with refractory advanced colorectal cancer. J Clin Oncol. 2008; 26:1427–34. https://doi.org/10.1200/JCO.2007.12.4602. [PubMed].
22. Steffens M, Paul T, Hichert V, Scholl C, von Mallek D, Stelzer C, Sörgel F, Reiser B, Schumann C, Rüdiger S, Boeck S, Heinemann V, Kächele V, et al. Dosing to rash?--The role of erlotinib metabolic ratio from patient serum in the search of predictive biomarkers for EGFR inhibitor-mediated skin rash. Eur J Cancer. 2016; 55:131–39. https://doi.org/10.1016/j.ejca.2015.11.022. [PubMed].
23. Fuerst D, Parmar S, Schumann C, Rüdiger S, Boeck S, Heinemann V, Kaechele V, Stiebel K, Paul T, Seufferlein T, Mytilineos J, Stingl JC. HLA polymorphisms influence the development of skin rash arising from treatment with EGF receptor inhibitors. Pharmacogenomics. 2012; 13:1469–76. https://doi.org/10.2217/pgs.12.122. [PubMed].
24. Parmar S, Schumann C, Rüdiger S, Boeck S, Heinemann V, Kächele V, Seeringer A, Paul T, Seufferlein T, Stingl JC. Pharmacogenetic predictors for EGFR-inhibitor-associated skin toxicity. Pharmacogenomics J. 2013; 13:181–88. https://doi.org/10.1038/tpj.2011.51. [PubMed].
25. Hasheminasab SM, Tzvetkov MV, Schumann C, Rüdiger S, Boeck S, Heinemann V, Kächele V, Steffens M, Scholl C, Hichert V, Seufferlein T, Brockmöller J, Stingl JC. High-throughput screening identified inherited genetic variations in the EGFR pathway contributing to skin toxicity of EGFR inhibitors. Pharmacogenomics. 2015; 16:1605–19. https://doi.org/10.2217/pgs.15.97. [PubMed].
26. Paul T, Schumann C, Rüdiger S, Boeck S, Heinemann V, Kächele V, Steffens M, Scholl C, Hichert V, Seufferlein T, Stingl JC. Cytokine regulation by epidermal growth factor receptor inhibitors and epidermal growth factor receptor inhibitor associated skin toxicity in cancer patients. Eur J Cancer. 2014; 50:1855–63. https://doi.org/10.1016/j.ejca.2014.04.026. [PubMed].
27. Hichert V, Scholl C, Steffens M, Paul T, Schumann C, Rüdiger S, Boeck S, Heinemann V, Kächele V, Seufferlein T, Stingl J. Predictive blood plasma biomarkers for EGFR inhibitor-induced skin rash. Oncotarget. 2017; 8:35193–204. https://doi.org/10.18632/oncotarget.17060. [PubMed].
30. Liu Q, Wu DH, Han L, Deng JW, Zhou L, He R, Lu CJ, Mi QS. Roles of microRNAs in psoriasis: Immunological functions and potential biomarkers. Exp Dermatol. 2017; 26:359–67. https://doi.org/10.1111/exd.13249. [PubMed].
32. Jiang C, Hu X, Alattar M, Zhao H. miRNA expression profiles associated with diagnosis and prognosis in lung cancer. Expert Rev Anticancer Ther. 2014; 14:453–61. https://doi.org/10.1586/14737140.2013.870037. [PubMed].
34. Boeri M, Verri C, Conte D, Roz L, Modena P, Facchinetti F, Calabrò E, Croce CM, Pastorino U, Sozzi G. MicroRNA signatures in tissues and plasma predict development and prognosis of computed tomography detected lung cancer. Proc Natl Acad Sci U S A. 2011; 108:3713–18. https://doi.org/10.1073/pnas.1100048108. [PubMed].
35. Ma Y, Chen Y, Lin J, Liu Y, Luo K, Cao Y, Wang T, Jin H, Su Z, Wu H, Chen X, Cheng J. Circulating miR-31 as an effective biomarker for detection and prognosis of human cancer: a meta-analysis. Oncotarget. 2017; 8:28660–71. https://doi.org/10.18632/oncotarget.15638. [PubMed].
36. Ma Y, Pan X, Xu P, Mi Y, Wang W, Wu X, He Q, Liu X, Tang W, An HX. Plasma microRNA alterations between EGFR-activating mutational NSCLC patients with and without primary resistance to TKI. Oncotarget. 2017; 8:88529–36. https://doi.org/10.18632/oncotarget.19874. [PubMed].
37. Mussnich P, Rosa R, Bianco R, Fusco A, D'Angelo D. MiR-199a-5p and miR-375 affect colon cancer cell sensitivity to cetuximab by targeting PHLPP1. Expert Opin Ther Targets. 2015; 19:1017–26. https://doi.org/10.1517/14728222.2015.1057569. [PubMed].
38. Hichert V, Steffens M, Paul T, Scholl C, Parmar S, Rüdiger S, Schumann C, Seufferlein T, Stingl JC. Gene regulatory biomarker identification for skin toxicities induced by EGFR inhibitor treatment. Int J Clin Pharmacol Ther. 2015; 53:1056–58. https://doi.org/10.5414/CPXCES14EA08. [PubMed].
39. Zhao F, Wei C, Cui MY, Xia QQ, Wang SB, Zhang Y. Prognostic value of microRNAs in pancreatic cancer: a meta-analysis. Aging (Albany NY). 2020; 12:9380–404. https://doi.org/10.18632/aging.103214. [PubMed].
40. Ferreira P, Roela RA, Lopez RVM, Del Pilar Estevez-Diz M. The prognostic role of microRNA in epithelial ovarian cancer: a systematic review of literature with an overall survival meta-analysis. Oncotarget. 2020; 11:1085–95. https://doi.org/10.18632/oncotarget.27246. [PubMed].
41. Shen H, Zhu F, Liu J, Xu T, Pei D, Wang R, Qian Y, Li Q, Wang L, Shi Z, Zheng J, Chen Q, Jiang B, Shu Y. Alteration in Mir-21/PTEN expression modulates gefitinib resistance in non-small cell lung cancer. PLoS One. 2014; 9:e103305. https://doi.org/10.1371/journal.pone.0103305. [PubMed].
42. Sabry D, El-Deek SEM, Maher M, El-Baz MAH, El-Bader HM, Amer E, Hassan EA, Fathy W, El-Deek HEM. Role of miRNA-210, miRNA-21 and miRNA-126 as diagnostic biomarkers in colorectal carcinoma: impact of HIF-1α-VEGF signaling pathway. Mol Cell Biochem. 2019; 454:177–89. https://doi.org/10.1007/s11010-018-3462-1. [PubMed].
43. Kucuksayan H, Akgun S, Ozes ON, Alikanoglu AS, Yildiz M, Dal E, Akca H. TGF-β-SMAD-miR-520e axis regulates NSCLC metastasis through a TGFBR2-mediated negative-feedback loop. Carcinogenesis. 2019; 40:695–705. https://doi.org/10.1093/carcin/bgy166. [PubMed].
44. Bae YU, Son Y, Kim CH, Kim KS, Hyun SH, Woo HG, Jee BA, Choi JH, Sung HK, Choi HC, Park SY, Bae JH, Doh KO, Kim JR. Embryonic Stem Cell-Derived mmu-miR-291a-3p Inhibits Cellular Senescence in Human Dermal Fibroblasts Through the TGF-β Receptor 2 Pathway. J Gerontol A Biol Sci Med Sci. 2019; 74:1359–67. https://doi.org/10.1093/gerona/gly208. [PubMed].
45. Fedenko ES, Elisyutina OG, Filimonova TM, Boldyreva MN, Burmenskaya OV, Rebrova OY, Yarilin AA, Khaitov RM. Cytokine gene expression in the skin and peripheral blood of atopic dermatitis patients and healthy individuals. Self Nonself. 2011; 2:120–24. https://doi.org/10.4161/self.2.2.16939. [PubMed].
46. Xu N, Meisgen F, Butler LM, Han G, Wang XJ, Söderberg-Nauclér C, Ståhle M, Pivarcsi A, Sonkoly E. MicroRNA-31 is overexpressed in psoriasis and modulates inflammatory cytokine and chemokine production in keratinocytes via targeting serine/threonine kinase 40. J Immunol. 2013; 190:678–88. https://doi.org/10.4049/jimmunol.1202695. [PubMed].
47. Guttman-Yassky E, Mita A, De Jonge M, Matthews L, McCarthy S, Iwata KK, Verweij J, Rowinsky EK, Krueger JG. Characterisation of the cutaneous pathology in non-small cell lung cancer (NSCLC) patients treated with the EGFR tyrosine kinase inhibitor erlotinib. Eur J Cancer. 2010; 46:2010–19. https://doi.org/10.1016/j.ejca.2010.04.028. [PubMed].
49. Anandappa G, Lampis A, Cunningham D, Khan KH, Kouvelakis K, Vlachogiannis G, Hedayat S, Tunariu N, Rao S, Watkins D, Starling N, Braconi C, Darvish-Damavandi M, et al. miR-31-3p Expression and Benefit from Anti-EGFR Inhibitors in Metastatic Colorectal Cancer Patients Enrolled in the Prospective Phase II PROSPECT-C Trial. Clin Cancer Res. 2019; 25:3830–38. https://doi.org/10.1158/1078-0432.CCR-18-3769. [PubMed].
50. Laurent-Puig P, Grisoni ML, Heinemann V, Liebaert F, Neureiter D, Jung A, Montestruc F, Gaston-Mathe Y, Thiébaut R, Stintzing S. Validation of miR-31-3p Expression to Predict Cetuximab Efficacy When Used as First-Line Treatment in RAS Wild-Type Metastatic Colorectal Cancer. Clin Cancer Res. 2019; 25:134–41. https://doi.org/10.1158/1078-0432.CCR-18-1324. [PubMed].
51. Yan S, Xu Z, Lou F, Zhang L, Ke F, Bai J, Liu Z, Liu J, Wang H, Zhu H, Sun Y, Cai W, Gao Y, et al. NF-κB-induced microRNA-31 promotes epidermal hyperplasia by repressing protein phosphatase 6 in psoriasis. Nat Commun. 2015; 6:7652. https://doi.org/10.1038/ncomms8652. [PubMed].
52. Mascia F, Mariani V, Girolomoni G, Pastore S. Blockade of the EGF receptor induces a deranged chemokine expression in keratinocytes leading to enhanced skin inflammation. Am J Pathol. 2003; 163:303–12. https://doi.org/10.1016/S0002-9440(10)63654-1. [PubMed].
53. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011; 17:10–12.
54. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009; 10:R25. https://doi.org/10.1186/gb-2009-10-3-r25. [PubMed].
All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 License.