Next Article in Journal
Early Prediction of Tumor Response to Neoadjuvant Chemotherapy and Clinical Outcome in Breast Cancer Using a Novel FDG-PET Parameter for Cancer Stem Cell Metabolism
Previous Article in Journal
Bioinformatics Data Mining Repurposes the JAK2 (Janus Kinase 2) Inhibitor Fedratinib for Treating Pancreatic Ductal Adenocarcinoma by Reversing the KRAS (Kirsten Rat Sarcoma 2 Viral Oncogene Homolog)-Driven Gene Signature
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JPM
Volume 10
Issue 3
10.3390/jpm10030131
jpm-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Liquid lncRNA Biopsy for the Evaluation of Locally Advanced and Metastatic Squamous Cell Carcinomas of the Head and Neck

by
Izabela Łasińska
1,2,*,†,
Tomasz Kolenda
3,4,*,†,
Kacper Guglas
4,5,
Magda Kopczyńska
3,4,
Joanna Sobocińska
3,
Anna Teresiak
4,
Norbert Oksza Strzelecki
3,
Katarzyna Lamperska
4,
Andrzej Mackiewicz
3,6 and
Jacek Mackiewicz
1,6,7,*
1
Department of Medical and Experimental Oncology, Heliodor Swiecicki Clinical Hospital, Poznan University of Medical Sciences, 16/18 Grunwaldzka Street, 60-786 Poznan, Poland
2
Specialist Nursing Laboratory, Faculty of Medicine and Health Science, University of Zielona Góra, Energetyków Street 2, 65-00 Zielona Gora, Poland
3
Department of Cancer Immunology, Chair of Medical Biotechnology, Poznan University of Medical Sciences, 8 Rokietnicka Street, 60-806 Poznan, Poland
4
Laboratory of Cancer Genetics, Greater Poland Cancer Centre, 15 Garbary Street, room 5025, 61-866 Poznan, Poland
5
Postgraduate School of Molecular Medicine, Medical University of Warsaw, 61 Zwirki i Wigury Street, 02-091 Warszawa, Poland
6
Department of Diagnostics and Cancer Immunology, Greater Poland Cancer Centre, 15 Garbary Street, 61-866 Poznan, Poland
7
Department of Oncology, Poznan University of Medical Sciences, 82-84 Szamarzewskiego, 60-569 Poznan, Poland
*
Authors to whom correspondence should be addressed.
†
These authors have contributed equally to this work.
J. Pers. Med. 2020, 10(3), 131; https://doi.org/10.3390/jpm10030131
Submission received: 22 July 2020 / Revised: 3 September 2020 / Accepted: 14 September 2020 / Published: 16 September 2020
Browse Figures
Versions Notes

Abstract

:
Background: Long non-coding RNA (lncRNA) are RNA molecules that are more than 200 nucleotides long and have the ability to modify the activity of genes. They can be found in both healthy and cancer tissues, as well as in plasma, saliva and other bodily fluids. They can also be used as biomarkers of early detection, prognosis and chemotherapy resistance in several cancer types. Treatment of head and neck squamous cell carcinoma (HNSCC) patients with locally advanced disease is still difficult, and choice of treatment should be based on more precise and available biomarkers, such as those obtained from a liquid biopsy. For improvement of treatment efficacy, identification and clinical implementation of new biomarkers are of the utmost importance. Methods: Plasma samples drawn before (p1) and three cycles post (p2) (TPF: docetaxel, cisplatin, 5-fluorouracil/PF: cisplatin, 5-fluorouracil) chemotherapy from 53 HNSCC patients (17 with locally advanced and 36 with metastatic disease) and 14 healthy volunteers were studied. Expression levels of 90 lncRNA expression were analyzed using the qRT-PCR method, and the obtained results were compared between proper groups. Statistical analyses were carried out using Jupyter Notebooks (5.7.2), Python (ver. 3.6) and GraphPad Prism 8. Results: The study demonstrated the differences between the expressions of several lncRNA in cancer patients’ and healthy volunteers’ plasma, as well as between locally advanced and metastatic patients’ groups. A correlation between the response to systemic therapy and lncRNA expression levels was observed. Patients with a (high/low) expression of Alpha 250 and Emx2os showed statistically significant differences in progression free survival (PFS), as well as for overall survival (OS) depending on the level of Alpha 250, snaR, SNHG1. The univariate and multivariate Cox regression model showed Alpha 250 as the best prognostic factor for HNSCC patients. Conclusions: Liquid biopsies based on lncRNAs are promising diagnostic tools that can be used to differentiate between those with cancer and healthy individuals. Additionally, they can also serve as biomarkers for chemotherapy resistance. An identified, circulating lncRNA Alpha 250 seems to prove the best prognostic biomarker, associated with extended PFS and OS, and should be validated in a larger cohort in the future.

1. Introduction

Head and neck cancers are mainly squamous cell carcinomas (HNSCCs, head and neck squamous cell carcinomas) with approximately 529,000 new cases diagnosed worldwide in 2012 [1]. The etiology of HNSCCs is complex and depends mostly on external risk factors, such as tobacco and/or alcohol use in older patient groups and infection with human papillomavirus (HPV) in younger patients. Furthermore, Epstein–Barr virus infection has been shown to play a role in the case of nasopharyngeal cancers [2]. According to The Cancer Genome Atlas (TCGA) the main genetic alterations present in HNSCC are associated with protein-coding genes involved in cell cycle (CDKN2A, TP53, CCND1), growth signaling (EGFR), survival (PIK3CA, PTEN), WNT signaling (FAT1, AJUBA, NOTCH1) and epigenetic regulation (KMT2D, NSD1). In human papillomavirus (HPV) infection, tumors showed changes in p16INK4A, TP53, RB or PI3CA [3].
Modern sequencing methods, as well as the creation and sharing of databases, enable faster and easier searching for molecular markers which can be used for HNSCC diagnostic, prognostic and prediction purposes [4]. The improvement of HNSCC treatment is still something to be achieved and current investigations are focusing on immunotherapy and specific biomarkers such as PD-L1/PD-L2 expression or the interferon-gamma gene signature [5]. However, a large part of the research on biomarkers in HNSCC concerns a transcriptome that does not encode a protein, which includes regulatory RNAs such as miRNA or lncRNA [6,7]. Moreover, when it comes to HNSCC serum, saliva and plasma are good sources of potential biomarkers [8].
It is accepted that only about 2% of the human genome encodes proteins, while the remaining 98% is classified as ‘junk’ DNA [9,10]. Further studies have shown that ‘junk’ DNA might be processed as RNA with non-protein coding potential (non-coding RNA, ncRNA) and play a role in cellular processes. It can also serve as a regulator of gene action as part of epigenetic mechanisms [11,12,13]. The changes in ncRNAs are the important elements of tumorigenesis and cancer progression [14]. One of the known types of ncRNAs is long non-coding RNAs (lncRNAs). lncRNAs are longer than 200 nucleotide transcripts. They can modify histone complexes and nucleus architecture and have an impact on RNA stability through epigenetic modifications, RNA transcription, translation and post-translational modification. Additionally, they interact with signaling and regulatory proteins and they act as “molecular sponges” [15,16]. There are about 50,000 types of lncRNAs in humans, which are transcribed mostly from introns. They are characterized by low conservatives across our species [17,18]. lncRNA transcripts are located in the nucleus, cytoplasm, or when circulating molecules form in various bodily fluids, such as plasma, saliva or gastric fluid [19,20,21]. It is estimated that about 40–90% of lncRNAs are translated, but these proteins are likely unstable [22].
Different lncRNAs take part in many important cellular processes, such as cell cycle, response to stress (including response to chemotherapeutics and irradiation), cellular senescence, cellular communication, hematopoiesis, embryogenesis, virus–host interactions and tissue microenvironment modification [23,24,25,26,27,28,29]. In cancerogenesis, changes in the expression levels of lncRNA are reported in tissues, as well as in bodily fluids [7]. lncRNA are responsible for the invasion, proliferation and metastatic potential in cancer cells as well as apoptosis in many cancers [15,30,31].
In HNSCC, more and more lncRNAs are linked with the biology of this cancer and are indicated as potential biomarkers [7,31]. For example, ZFAS1 regulates miR-150-5p and by its influence on cellular phenotype. Another lncRNA, EGOT, is connected with HPV infection. Both of these lncRNAs could serve as prognostic biomarkers in HNSCC [32,33]. However, the exact function of the most lncRNAs is still unknown and needs to be further investigated. It is acknowledged that patients’ response to treatment depends on their unique cancer genetic pattern, which is manifested by its specific phenotype and behavior [34]. lncRNAs seem to be an important regulator of cancer cells and are promising therapeutic targets [15].
Treatment of locally advanced HNSCC includes surgical resection, radiotherapy, or radiochemotherapy, which may be preceded with induction chemotherapy. In recurrent or metastatic disease, the treatment is based on chemotherapy or immunotherapy [35]. However, the capacity to make a molecular-based decision about treatment strategy is limited [36].
The monitoring of patients’ response to the chosen therapeutic strategies and prediction of the probability of treatment response or disease progression during cancer therapy is one of the most important challenges in contemporary oncology [36,37]. One promising strategy is liquid biopsy, used to detect circulating nucleic acids, DNA or RNA, which are linked with specific patients’ condition and their response to treatment. This method could potentially also specify the survival prognosis of individuals [38]. The different types of RNAs, including protein-coding and non-coding molecules, include: miRNA, tRNA-derived RNA fragments (tRFs), circRNA, lncRNAs. All of these, and others, can be found in serum or plasma. Furthermore, they are present as free acids, complexes with proteins or lipids, or incorporated into extracellular vesicles, as for example into exosomes [39]. Many studies indicate the potential role of circulating miRNAs as biomarkers, but the role of other non-coding RNAs, such as lncRNAs, seems to be unclear and should be explored more thoroughly to establish the best possible biomarker panel [39].
In this study, based on the commercially available lncRNA qRT-PCR plate, the expression of 90 transcripts was examined in plasma samples of healthy volunteers and HNSCC patients. Next, the clinical usefulness of expression levels of plasma lncRNAs as a diagnostic and prognostic tool in locally advanced and palliative HNSCC patients was evaluated.

2. Materials and Methods

2.1. Patients and Volunteers

In this prospective analysis, 53 plasma samples from HNSCC patients and 14 healthy volunteers (18–80 years old) were obtained between April 2016 and May 2017. The control group consisted of 9 men and 5 women with no comorbidities. Seventeen naive HNSCC patients were treated with chemotherapy, induction chemotherapy first then followed by chemoradiotherapy or radiotherapy. The remaining 36 HNSCC patients were treated before our study, and later received palliative chemotherapy only. All of the examined patients were HPV (human papillomavirus) negative. The patients’ characteristics are shown in Supplementary Table S1.

2.2. Treatment

All patients in the induction group (locally advanced disease) were treated with a TPF regimen, which contained the following: docetaxel (T—75 mg/m2 day 1), cisplatin (P—75 mg/m2 day 1), and 5-fluorouracil (F—750 mg/m2 days 1–5). It was repeated three times during a 21 day period. Following computed tomography (CT), they were treated at the Radiotherapy Department. All patients in the palliative group (recurrent and/or metastatic disease) were treated using a PF regimen, which contained a cisplatin (P—100 mg/m2 day 1) and 5-fluorouracil (F—1000 mg/m2 days 1–5) regimen repeated three times during a 21 day period (Cetuximab was not reimbursed in Poland during this trial). After the treatment course of 3 cycles of chemotherapy, CT was performed and appropriate treatment was administered, but this aspect is beyond the scope of this study. When a response or stable disease were observed, in line with RECIST (Response Evaluation Criteria in Solid Tumors) 1.1 [40], the patients received up to 3 more cycles of PF. For patients with progressive disease (PD) according to RECIST 1.1, after the 3 cycles of chemotherapy mentioned above, a new chemotherapy regimen (Methotrexate 40 mg/m2 day 1 weekly) was administered. Prior to each therapy cycle, the patients were tested for adequate marrow, kidney and liver function. Additionally, a performance status was evaluated. The treatment response according to RECIST 1.1 depending on type of chemotherapy, induction or palliative, is presented in Supplementary Table S2.

2.3. Blood Sampling Design

Whole blood samples (5 mL) were taken from all patients receiving induction and palliative chemotherapy at two points in time: (i) before the start of treatment (P1), and (ii) after 3 cycles of chemotherapy (P2). In the volunteer group, blood samples were taken only once, namely in April 2016. All participating patients signed consent forms pertaining to the study’s aim, as well as personal and genetic data protection. The study methodologies conformed to the standards of the Declaration of Helsinki. The study did not violate the rights of other persons or institutions. The Bioethical Committee of Medical University of Poznan approved this study in agreement No. 152/2016.

2.4. Sample Collection

Pure blood samples were collected into tubes containing ethylenediaminetetraacetic acid (EDTA) (SARSTED, Monovette EDTA K) and immediately centrifuged (10 min; 3500 RPM, RT). The upper plasma phase was transferred to a new tube without disturbing the intermediate buffy coat layer and frozen at a temperature of −80 °C.

2.5. RNA Isolation

Before RNA isolation was conducted, all samples were checked for hemolysis using conventional macroscopic classification methods for samples [41]. The plasma samples were later thawed and only one cycle of thawing was conducted before taking 200 μL of material for further processing. The total RNA was isolated from plasma samples using miRNeasy Serum/Plasma Kit (Qiagen) according to the protocol for total RNA. A NanoDrop spectrophotometer (Thermo Scientific) was used to estimate the quality and quantity of the RNA samples, and then they were taken for further analysis, while the remaining RNA samples were stored at a temperature of –80 °C.

2.6. cDNA Synthesis and qRT-PCR Reaction

In the study, the 90 lncRNAs potentially connected with cancer and well-annotated in the lncRNA database (www.lncrnadb.org) were analyzed using the commercially available LncProfiler qPCR Array Kit (System Biosciences). The reaction consisted of three steps: (i) poly-A tailing, (ii) annealing anchor dT adaptor, (iii) complementary cDNA synthesis [42]. For the poly-A tailing step, 5 μL of RNA (about 100 ng) was used and mixed with 2 μL 5x PolyA Buffer, 1 μL MnCl2, 1.5 μL ATP and 0.5 μL polyA polymerase and incubated for 30 min at a temperature of 37 °C. Next, 0.5 μL of Oligo(dT) adapter was added and the reaction was heated for 5 min at 60 °C, and later cooled to room temperature. Finally, 4 μL of reverse transcription (RT) Buffer, 2 μL deoxynucleotide mix, 1.5 μL 0.1 M dithiothreitol (DTT), 1.5 μL Random Primer Mix and 1 μL reverse transcriptase were added and incubated for one hour at 42 °C, followed by 10 min at 95 °C.
cDNA was used for the qRT-PCR reaction using LightCycler 480 SYBR Green I Master buffer (Roche) and lncRNA primers from Primer Plate (component of the LncProfiler qPCR Array Kit) according to the manufacturer’s protocols.
The quantitative PCR reaction was performed using a preincubation step (95 °C for 10 min), 60 cycles of two-step amplification (95 °C for 15 s and 60 °C for 60 s) and a melting step with LightCycler 96 (Roche). All qRT-PCR data of 90 lncRNAs were analyzed by calculating the ΔCt, normalized against the mean expression of reference genes from the LncProfiler qPCR Array Kit (SBI), Supplementary Table S3. The fold-change in lncRNA expression was determined using the equation 2−ΔCt and after that was compared to the appropriate group.

3. Statistical Analysis

Statistical analysis was performed with Jupyter Notebooks (5.7.2), Python (ver. 3.6), native libraries and ‘Lifelines’ library v.0.14.6 (Kaplan-Meier, COX) or GraphPad Prism 8. All data are presented as medians with 95% confidence intervals. Expression profiles of lncRNA were compared as follows: between HNSCC patients vs. healthy volunteers; HNSCC patients divided into two categories: (i) induction (locally advanced) vs. (ii) palliative treatment; localization of primary cancer lesion as well as treatment response in the palliative group were compared. The Shapiro–Wilk test was used to estimate data distribution, and, next, the T-test or Mann–Whitney tests were used to compare appropriate groups.
The correlation between clinical or laboratory parameters and lncRNA expression level was assessed using the Mann–Whitney and Kruskal–Wallis tests.
For diagnostic effectiveness of lncRNAs, ROC (receiver operating characteristic) calculation was performed and an optimal cut-off point was assessed according to the highest accuracy (minimal false negative and false positive rates). For prognostic purposes the value of AUC (area under the curve) of ROC for each lncRNA was calculated. For progression free survival (PFS) and overall survival (OS) in the palliative group, Kaplan–Meier curves were calculated.
The Cox proportional hazards regression model was used to determine which lncRNA is significantly associated with OS and PFS, and the results were presented as hazard ratios (HR) and 95% confidence intervals (CI). All of the statistical analyses were performed as two-tailed tests and considered significant at p < 0.05.

4. Results

4.1. Comparison of lncRNA Expression Levels between Healthy Volunteers and All HNSCC Patients

First, plasma lncRNA expression levels were compared in all HNSCC patients and all healthy volunteers. A significantly higher expression of 34 lncRNAs with a fold change of 0.49–2.1855 × 104 was found in cancer patients. lncRNAs included: 21A, Air, ANRIL, Alpha 280, BACE1AS, Dio3os, E2F4 antisense, Emx2os, EvF1 and EVF2, GAS5, H19, HOXA6as, Jpx, Kcnq1ot1, lincRNA-RoR, mascRNA, MER11C, ncR-uPAR, NRON, p53 mRNA, RNCR3, SCA8, SNHG1, SNHG3, SNHG4, SNHG6, SRA, ST7OT, TEA ncRNAs, TncRNA, UCA1, WT1-AS, YRNA1 and Zfas1. All data are presented in Figure 1 and Supplementary Table S4. Next, ROC analysis was performed with an assessment of AUC for a prognostic value of these 34 selected lncRNAs. The differentiation between cancer patients and healthy individuals was possible with a sensitivity of 38.5–100% and a specificity of 57.1–100%. Details are presented in Supplementary Table S4.

4.2. Comparison of lncRNA Expression Levels between Healthy Volunteers and Recurrent and/or Metastatic HNSCC Patients

The analysis revealed a significantly higher expression of 41 lncRNAs in recurrent and/or metastatic HNSCC patients compared to healthy volunteers with a fold change of 3.2–2.7 × 103 for 21A, Air, Alpha 280, ANRIL, BACE1AS, Dios3os, E2F4 antisense, Emx2os, GAS5, H19, HAR1B, HOXA6as, HULC, Jpx, Kcnq1ot1, KRASP1, lincRNA-RoR, LUST, MALAT1, mascRNA, MEG9, MER11C, ncR-uPAR, NEAT1, NRON, p53 mRNA, PSF inhibiting RNA, PTENP1, RNCR3, SCA8, SNHG1, SNHG3, SNHG4, SNHG6, SRA, ST7OT, TncRNA, Tsix, UCA1, YRNA1 and Zfas1. All data are presented in Figure 2 and Supplementary Table S5. The ROC analysis with AUC estimation for prognostic value of 41 up-regulated lncRNAs was calculated. For the selected lncRNAs, the differentiation between studied groups was possible at a sensitivity and specificity between 38.9–100% and 42.9–100%, respectively (Supplementary Table S5).

4.3. Comparison of lncRNA Expression Levels between Healthy Volunteers and Locally Advanced HNSCC Patients

In cancer patients, a significantly higher expression of 47 lncRNAs with fold change 1.96–5.101 × 104 compared to healthy individuals was observed. However, a significantly lower expression of 2 lncRNAs (H19 upstream conserved 1 & 2, HOXA6as) with a fold change of 0.30–0.32 compared to healthy individuals was indicated, as shown in Figure 3 and Supplementary Table S6. The differentiation between study groups with high and low lncRNA expressions with ROC analysis was possible at a sensitivity and specificity of 37.5–100% and 71.4–100%, respectively. Details can be found in Supplementary Table S6.

4.4. Plasma lncRNA Expression Differs Depending on Primary Tumor Localization

The association between the expression level of plasma lncRNAs and the primary site of the tumor was indicated. Five lncRNAs, HAR1B, Jpx, LUST, NEAT1 and TncRNA, showed significantly changed expressions when pharynx and oral cavity localizations were compared. For cancers with pharyngeal and laryngeal localizations, lncRNA Air was only significantly altered. The most frequent changes in lncRNA expressions were detected between patients with salivary gland cancers in comparison to other tumor sites. All data are presented in Figure 4 and Supplementary Table S7.

4.5. Comparison of lncRNA Expression Levels According to Treatment Response in HNSCC Patients Receiving Palliative Chemotherapy

The different responses of the studied patients to treatment were assessed according to RECIST 1.1 criteria (Supplementary Table S2). Additionally, lncRNA expressions in specific subgroups were analyzed. In HNSCC, in the group receiving palliative chemotherapy, patients with PD (progressive disease) showed significantly high expressions of 14 lncRNAs compared to patients with DCR (disease control rate, CR+PR+SD), 7SK (p = 0.0142), CAR Intergenic 10 (p = 0.0022), HAR1A (p = 0.0279), NOR (p = 0.0171), IGF2AS (p = 0.017), Kcnq1ot1 (p = 0.0316), KRASP1 (p = 0.0171), L1PA16 (p = 0.0266), LOC285194 (p = 0.0374), MEG (p = 0.0012), NTT (p = 0.0192), Zeb2NAT (p = 0.0112), SAF (p = 0.0494), and SNHG6 (p = 0.0494), as shown in Figure 5. All patients responded to induction chemotherapy. Accordingly, it was not possible to conduct a comparative study.

4.6. Changes in lncRNA Expression Levels during Induction Chemotherapy and Palliative Chemotherapy

Subsequent analysis of changes in lncRNA expressions in a palliative group of eight patients with DCR at two time points—before treatment (P1) and after three cycles of chemotherapy (P2)—was performed. Significant downregulation of three lncRNAs, Jpx (p = 0.02), PRINS (p < 0.001) and PTENP1 (p < 0.0001), was observed, as shown in Figure 6A. In the remaining 28 patients included in the palliative group, a progression of disease (PD) was observed. In these patients, statistically significant downregulation of lncRNAs ANRIL (p = 0.045), MALAT1 (p = 0.0217) and NDM29 (p = 0.043) between two time points—(P1) and (P2)—was found, as shown in Figure 6B.
All patients (n = 17) with locally advanced HNSCC treated with induction chemotherapy presented a DCR response during this study. Analysis of lncRNA expressions in these patients revealed that only PRINS was significantly downregulated (p = 0.0131) after three cycles of chemotherapy (P2) compared to pretreatment (P1) levels (Figure 6C).

4.7. Assessment of lncRNA Expression in HNSCC Patients Receiving Palliative Chemotherapy Can Serve as a Prognostic Biomarker

Patients included in the palliative group were divided into two subgroups with high or low expressions of each lncRNA according to the mean expression level as a cut-off. Moreover, Kaplan–Meier curves were calculated for OS and PFS for each of the 90 examined lncRNAs.
Two of the lncRNAs, Alpha 250 and Emx2os, manifested a statistically significant difference in PFS in respect to high/low levels of expression. For patients with high expressions of Alpha 250 or Emx2os, a five month (95% CI (0.4–0.9) and 95% CI (0.1–0.9), respectively) median of PFS was noticed. Moreover, median PFS was not reached in the group with low expressions of these lncRNAs, and patients were characterized by a longer time without progression, as shown in Figure 7A.
Next, the OS in respect to high/low expression levels of all lncRNAs was examined, and three of them, snaR, SNHG1 and Alpha 250, showed statistically significant differences between groups of patients. In patients with low expressions of snaR or SNHG1, the survival of patients was shorter compared to the group with higher expressions of these lncRNAs: 30 month (95% CI (0.15–0.8), p = 0.03) and 22 month (95% CI (0.2–0.8), p = 0.043) mean of OS was respectively observed. In both cases, a median OS was not reached in patients with high expressions.
In the group of patients with high expressions of Alpha 250, a 17 month median OS was noticed (95% CI (0.2–0.9), p = 0.033). In the low expression group, the median was not reached, and patients revealed a longer OS compared to the group with high expression levels of this lncRNA, as shown in Figure 7B.
For further analysis, univariate and multivariate Cox regression models were performed to assess the independence of the prognostic value of Alpha 250, snaR and Emx2os, which were later indicated using Kaplan–Meier curves as prognostic factors. Univariate and multivariate models showed that Alpha 250 is only a statistically significant factor for probability of OS and PFS, and it is independent of age, neutrophil to lymphocyte ratio (NLR), body mass index (BMI) and primary localization of tumor. All data are presented in Table 1 and Table 2.

4.8. lncRNA Expression Levels and Clinical and Laboratory Parameters of HNSCC Patients

Expression levels of Alpha 250, Emx2os, snaR and SNHG1 were analyzed in terms of clinical and laboratory parameters that may be helpful in the assessment of patients’ prognosis. Significant differences between stage IVA and IVC for SNHG1 (0.0823 vs. 0.4601, p = 0.0401) and snaR in younger patients compared to older ones (2.514 vs. 12.1099, p = 0.0183) were observed. Regarding Alpha 250 and Emx2os, no significant differences in expression levels within the group of patients divided by age, sex, body mass index, nicotinism, alcoholism, localization of primary tumor, disease stage, neutrophil to lymphocyte ratio and platelets to lymphocyte ratio were noticed. All data are presented in Table 3.

5. Discussion

Approximately half of HNSCC cases are diagnosed in an advanced stage with poor outcome [43]. Accordingly, the identification of low invasive and simple biomarkers for these patients’ prognosis assessment is crucial. Recently, an intensive search, including liquid biopsy based on circulating nucleic acids, such as DNA and RNA [38,39,44], has been carried out. It was previously indicated that lncRNAs can be found in tumor tissue (including HNSCC) and in bodily fluids [44,45]. However, circulating lncRNA as potentially predictive and prognostic biomarkers in advanced or metastatic HNSCC were not reported.
The first important finding is that the level of plasma lncRNA expressions in HNSCC patients and healthy volunteers differs significantly. In addition, it directly correlates with the stage of the disease. Out of the 90 lncRNAs studied, 47 were high up- and 2 down-regulated in patients with locally advanced cancers. A total of 41 lncRNAs were up-regulated in patients with recurrent and/or metastatic cases, while 29 lncRNAs were down-regulated in locally advanced and recurrent and/or metastatic HNSCC. Altered lncRNA expressions enabled differentiation of HNSCC patients from a healthy population with high sensitivity and specificity, proving lncRNA to be a good biomarker. Yao et al. demonstrated the expression of lncRNAs in plasma and in cancer tissue using microarray and next-generation RNA-sequencing methods (NGS) validated by qRT-PCR. They reported changes in 432 lncRNA transcripts in plasma and 333 in tissue (fold changes > 4). However, the expression of only three lncRNA in the plasma of HNSCC patients was significantly higher than in healthy controls (HOXA11as, LINC00964 and MALAT1) and might be used as potential circulating biomarkers for early detection [45]. In our study, we observed changes in 34 out of 90 examined lncRNAs. In contrast to Yao et al., we did not observe expression changes in HOXA11as, LINC00964 and MALAT1 in all groups of HNSCC patients compared with healthy donors. However, up-regulation of HOXA11as and MALAT1 in patients with locally advanced and recurrent and/or metastatic HNSCC was indicated. Interesting, we observed multiple expression changes in lncRNAs associated with HOXA (homeobox A), a cluster referred to as HOXA3as, HOXA6as and HOXA11as, as well as with lncRNAs associated with SNHGs (small nucleolar RNA host genes) such as SNHG1, SNHG3, SNHG4, SNHG5 and SNHG6. However, until now there are no data to point to the role of HOXA3as and HOXA6as in cancer. Only HOXA11as is indicated as an important lncRNA in progression and metastasis of a number of different cancers by regulating proteins and miRNAs crucial in cellular processes such as cell cycle, apoptosis, epithelial-to-mesenchymal transition (EMT) process, invasion and metastasis [46,47]. SNHG1 was also described as an important modulator of cancer progression [48]. In laryngeal squamous cell carcinoma (LSCC) the expression of SNGH1 was up-regulated, connected with cell proliferation, EMT process, metastasis and poor patients’ prognosis [49]. SNHG1 inhibits LSCC growth and metastasis in vitro by promoting YAP1 expression and Hippo signaling activity through sponging miR-375 [50]. Wang et al. indicated that SNHG3 regulates the proliferation and migration of LSCC by way of modulation within the miR-384/WEE1 axis [51]. SNHG4 is not described for HNSCC, but in lung cancers it affects proliferation, migration and invasiveness, as well as EMT through miR-98-5p [52]. Similarly, SNHG5 in nasopharyngeal carcinoma regulates the expression of both the miR-1179/HMGB3 axis and cancer progression [53]. Finally, SNHG6 was indicated in the meta-analysis as an oncogene connected with cancer progression and could be used as a promising prognostic biomarker [54]. Only two studies described its role in HNSCC. It was shown that SNHG6 was up-regulated in tongue cancer, and its knockdown inhibited cell viability and proliferation, induced apoptosis and prevented EMT in vitro [55]. Guo et al., observed the potential role of SNHG6 in the oncogenesis of HPV-positive oropharyngeal squamous cell carcinoma (OPSCC) [56]. Moreover, we found up-regulation of Y RNA-1 (RNY1) transcripts in our HNSCC patients. Our observation confirmed a previous report by Victoria Martinez et al., who, using NGS, demonstrated changes in Y RNAs transcripts in serum from HNSCC patients as compared to healthy donors [57]. In our previous study, we demonstrated that in BRAF-mutant metastatic melanoma patients treated with a BRAF inhibitor, a higher expression of Y RNA-1 was associated with worse prognosis [58]. Furthermore, Guglas et al., based on HNSCC cell lines, patient samples and TCGA data, indicated that Y RNA-1 was significantly down-regulated, while its higher expression was associated with worse patient survival [59]. More and more data indicate that YRNAs as well as YRNA-derived fragments are important regulators in tumorigenesis, and they have a potential role in diagnostics [60].
In the group of locally advanced HNSCC patients, two lncRNAs were downregulated compared to healthy volunteers, described above HOXA6as, and H19 upstream conserved 1 and 2. Drewell et al. indicated that conserved elements upstream of the H19 gene are transcribed and act as mesodermal enhancers [61]. There is no information about H19 upstream-conserved 1 and 2, nor about its role in the context of circulating transcripts and immune cells.
Unfortunately, we have not elucidated the nature of the observed lncRNA expression alterations and the potential function of them. Still, the origin of circulating lncRNA is unclear, and some evidence suggests it could be tumor tissue, tumor cells, something related to the immune system or even blood cells [62,63]. Yao et al., indicated that only 12 lncRNA expressions are changed in the plasma and tissue of HNSCC patients [45]. This suggests the remaining changed lncRNAs could demonstrate different origins, but this requires more studies. We also observed, like others, that the expression of lncRNAs is changed depending on the patient’s status. Differences between lncRNA expression levels between cancer patients and healthy individuals, and in two patient groups, are probably due to changes in physiognomy during cancer progression, earlier treatment of metastatic disease and treatment resistance, as well as the stage of the disease. The relatively small group included in the study may also be a possible explanation of the results obtained.
Our comparison of primary tumor localization and lncRNA expression indicated some similarities, as well as differences in transcript expressions. In salivary glands, different lncRNA expressions were reported as in other tumors. This might be due to the different etiology of the salivary glands and larynx, pharynx and oral cavity cancers, which is mainly due to alcohol abuse and nicotinism [64].
An important finding of our study shows the lncRNA expression pattern in a metastatic group of patients with primary resistance to treatment. Reports about circulating lncRNAs and resistance to therapy are limited, and, to the best of our knowledge, this is the first identification of 14 up-regulated lncRNAs (7SK, CAR intergenic 10, HAR1A, IGF2AS, Kcnq1ot1, KRASP1, L1PA16, LOC285194, MEG9, NTT, NRON, Zeb2NAT, SAF, SNHG6). Moreover, these RNAs could probably be used as biomarkers of primary resistance in HNSCC. Furthermore, during palliative chemotherapy, reduction of Jpx, PRINS and PTENP1 expression levels after three treatment cycles correlated with treatment response. In contrast, the reduction of lncRNAs named ANRIL, MALAT1 and NDM29 was reported in these patients while developing progressive disease. All patients treated with induction chemotherapy responded to treatment, and these patients have up-regulated PRINS measured during the second blood collection. It is possible that the aforementioned effect was related to the treatment applied.
It is difficult to compare our results with other studies since there are no similar reports available. Only one study examined the potential role of plasma lincRNA-p21, GAS5 and HOTAIR as biomarkers to predict a response to chemoradiotherapy, which demonstrated that only GAS5 was useful when measured after chemotherapy [65]. In contrast, our study did not confirm these results. However, we observed significant differences in expression levels of GAS5 between HNSCC patients and healthy volunteers, as well as between HNSCC patients—locally advanced and recurrent and/or metastatic disease compared to non-cancer.
Most studies of chemical agents on lncRNA expressions were based on in vitro models using cell lines. Guglas et al. based their study on HNSCC cell lines (SCC-040, SCC-25, FaDu and Cal27) observed changes in lncRNA expressions after exposure to irradiation and chemotherapeutic drugs. These changes were dose (5, 10 or 20 Gy) and type (cisplatin or doxorubicin) dependent, respectively. Cisplatin down-regulated CAR Intergenic 10, PRINS, and PTENP1 expressions, while doxorubicin down-regulated MALAT1 and others in surviving cells, compared to those non-treated with toxic agents [25]. Another study indicated that ANRIL, HOTAIR and MALAT1 were significantly down-regulated in cancer tissue compared to adjacent normal tissues. Moreover, expressions of these lncRNAs were reduced in laryngeal cell lines (Hep-2 and AMC-HN8) after exposition to cisplatin and paclitaxel [66]. Ren et al., using NGS, compared two nasopharyngeal cell lines—parental CNE-2 and paclitaxel-resistant—and revealed 2670 known and 4820 novel lncRNAs between them. From the altered lncRNAs, they selected the n375709 transcript, and its knock-down in the NPC 5–8F and 6–10B cell lines increased cells’ sensitivity to paclitaxel in vitro [67]. The regulation of chemosensitivity in nasopharyngeal cancers also displays lncRNA-ROR. It is up-regulated in NPC tissue as well as cell lines and affects cell proliferation, cell cycle, apoptosis, migration, and invasion. In vitro studies indicated that lncRNA-ROR regulates the EMT process. Moreover, cis-dichlorodiammineplatinum (II) (DDP) is induced in the cell line expression of lncRNA-ROR. Artificial down-regulation of lncRNA-ROR is caused by DDP influence on the p53 pathway [68].
Clinical and laboratory parameters showed a statistically significant correlation between low expressions of SNHG1 and stage IVA of disease. The expression of snaR was associated with a younger age in patients as well. It is difficult to conclude these findings, as there is a lack of information on lncRNA expression levels based on the stage of disease and age of patients in the literature.
Our final important finding is the observation that patients with down-regulated Alpha 250 and Emx2os display longer PFS compared to those with upregulation. In addition, the down-regulation of Alpha 250 was favorable also for longer OS. It is the first report of Alpha 250 in the context of cancer and HNSCC. The origin of Alpha 250 and Alpha 280 transcripts is the first intron of the Rps14 gene, which encodes S14 proteins, a part of the 40S ribosome subunit responsible for the translation and synthesis of proteins. Up-regulated expressions of Alpha 250 are linked with cell proliferation [69]. Tang et al. identified that low expressions of Emx2os are closely associated with longer DSF and OS ratio in LSCC patients based on in silico analysis of databases. They also indicated that Emx2os regulates the expression of CALCA and GABRG2 through miR-124 [70]. We have also observed that high expressions of snaR and SNHG1 was correlated with a favorable OS in HNSCC patients. Liang et al. showed that an overexpression of plasma snaR is correlated with poor OS in LSCC patients and positively correlated with TGF-β1. Moreover, snaR overexpression coupled with TGF-β1 treatment promoted cell proliferation, migration, and invasion in vitro [71]. SNHG1 was indicated also in LSCC but, as mentioned above, not only in tissue and cell lines but also as circulating ncRNAs. Its higher expression was connected with poorer OS [49,50]. However, our study based on profiling of lncRNA expression in melanoma patients indicated that a low expression of SNHG1 seems to be a good marker of DSF [58].
Univariate and multivariate Cox regression models for selected lncRNAs indicated that only plasma Alpha 250 proved to be a prognostic biomarker, and it is the first time this biomarker for chemotherapy resistance in plasma was identified.

6. Conclusions, Limitation and Future Perspective

In this study, we examined the expression profile of 90 different lncRNAs in plasma taken from 53 HNSCC patients (17 with locally advanced and 36 with metastatic disease) during chemotherapy and from 14 healthy volunteers. Based on these results, there are four major findings in our study: (i) plasma lncRNAs in HNSCC patients are markers for disease diagnosis; (ii) plasma lncRNA expression levels may serve as biomarkers of primary chemotherapy resistance in HNSCC patients; (iii) lncRNA expression levels bear prognostic biomarker potential for treatment efficacy; and (iv) lncRNA Alpha 250 is a prognostic biomarker for OS and PFS of HNSCC patients.
As observed, the plasma-based liquid biopsy and quantification of lncRNAs is a potentially clinically useful tool for the evaluation of locally advanced and metastatic HNSCC patients. However, our study is one of the first to describe plasma based lncRNAs in this specific group of patients, and it is difficult to discuss the results in a broader perspective. The lack of similar data for lncRNA Alpha 250, as well as for the most of indicated here lncRNAs in the context of HNSCC, makes our study limited. To verify our results, an independent study based on the large group of locally advanced and metastatic HNSCC patients should be undertaken.
The second important question is the source of the circulating transcripts. Unfortunately, we did not examine the tissue from analyzed patients, so it is difficult to answer if the observed lncRNAs are released from cancer tissue to the bloodstream or are a whole body response to the disease and treatment strategy. In the future, studies like these should be more comprehensive and profiling of transcripts in cancer tissue and circulating need to be obligatory. Moreover, future functional studies about the potential role of circulating lncRNAs could give a more comprehensive point of view in relation to whether these transcripts influence the potential regulatory role in other cells such as immune cells, or are only a consequence of still unknown processes. However, when it comes to biomarkers, the aforementioned questions are not so important, because their principal role is the description of the state of the patient and they provide clinical answers about a current or a potential treatment strategy.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4426/10/3/131/s1, Table S1: Baseline characteristic of patients included into the study. Table S2: Treatment response as defined by RECIST 1.1. Table S3: Examined lncRNAs and reference genes ncProfiler qPCR Array Kit (System Biosciences). Table S4: Median expression of lncRNA in healthy volunteers and HNSCC patients with Confidence Interval (CI), and prognostic value of lncRNAs with sensitivity and specificity features to distinguish the examined groups. Table S5: Median expression of lncRNA in volunteers and palliative treated HNSCC patients with Confidence Interval (CI), and prognostic value of lncRNAs with sensitivity and specificity features to distinguish the examined groups. Table S6: Median expression of lncRNA in volunteers and locally advanced HNSCC patients with Confidence Interval (CI), and prognostic value of lncRNAs with sensitivity and specificity features to distinguish the examined groups. Table S7: lncRNA expression according to the primary site of tumor.

Author Contributions

Authors’ individual contributions: conceptualization, I.Ł., T.K., and M.J.; methodology, T.K.; investigation, T.K., K.G., M.K., J.S., N.O.S. and I.Ł.; data curation, T.K., I.Ł. and J.M.; writing—original draft preparation, I.Ł. and T.K.; writing—review and editing, T.K., A.T., K.L., A.M. and J.M.; visualization, T.K., K.G. and M.K.; supervision, A.M., K.L. and J.M.; funding acquisition, K.L. and J.M.; I.Ł. and T.K. have contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Science Centre, Poland, allocated on the basis of decision no.: 2016/21/B/NZ7/01773.

Acknowledgments

Thanks to Greater Poland Cancer Center and Poznan University of Medical Sciences for supporting our work and providing a fully equipped laboratory to perform the necessary analyzes. This work was supported by the National Science Centre, Poland, allocated on the basis of decision no.: 2016/21/B/NZ7/01773. Language correction was conducted by Adam Slowinski from Universal Translations (www.biurotlumaczen.co).

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Availability of Data and Materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Ethics Approval

The study is based on human blood samples. All donors were informed about the aim of the study, personal data and genetic data protection, and they agreed to have their material used. The experiments were undertaken with the understanding and written consent of each subject. The study methodologies conformed to the standards set by the Declaration of Helsinki. The study did not violate the rights of other persons or institutions. The Bioethical Committee approved this study—agreement No. 152/2016.

References

  1. Shield, K.D.; Ferlay, F.; Jemal, A.; Sankaranarayanan, R.; Chaturvedi, A.K.; Bray, F.; Soerjomataram, I. The global incidence of lip, oral cavity, and pharyngeal cancers by subsite in 2012. CA Cancer J. Clin. 2017, 67, 51–56. [Google Scholar] [CrossRef]
  2. Tsao, S.W.; Tsang, C.M.; Lo, K.W. Epstein-Barr virus infection and nasopharyngeal carcinoma. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2017, 372, 20160270. [Google Scholar] [CrossRef] [PubMed]
  3. The Cancer Genome Atlas Network. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature 2015, 517, 576–582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Tonella, L.; Giannoccaro, M.; Alfieri, S.; Canevari, S.; De Cecco, L. Gene Expression Signatures for Head and Neck Cancer Patient Stratification: Are Results Ready for Clinical Application? Curr. Treat. Options Oncol. 2017, 18, 32. [Google Scholar] [CrossRef] [PubMed]
  5. Solomon, B.; Young, R.J.; Rischin, D. Head and neck squamous cell carcinoma: Genomics and emerging biomarkers for immunomodulatory cancer treatments. Semin. Cancer Biol. 2018, 52, 228–240. [Google Scholar] [CrossRef]
  6. Sannigrahi, M.K.; Sharma, R.; Panda, N.K.; Khullar, M. Role of non-coding RNAs in head and neck squamous cell carcinoma: A narrative review. Oral Dis. 2018, 24, 1417–1427. [Google Scholar] [CrossRef] [PubMed]
  7. Guglas, K.; Bogaczyńska, M.; Kolenda, T.; Ryś, M.; Teresiak, A.; Bliźniak, R.; Łasińska, I.; Mackiewicz, J.; Lamperska, K. lncRNA in HNSCC: Challenges and potential. Contemp. Oncol. 2017, 21, 259–266. [Google Scholar] [CrossRef] [Green Version]
  8. Arantes, L.M.R.B.; De Carvalho, A.C.; Melendez, M.E.; Carvalho, A.L. Serum, plasma and saliva biomarkers for head and neck cancer. Expert Rev. Mol. Diagn. 2018, 18, 85–112. [Google Scholar] [CrossRef]
  9. Orgel, L.E.; Crick, F.H.C. Selfish DNA: The ultimate parasite. Nature 1980, 284, 604–607. [Google Scholar] [CrossRef] [PubMed]
  10. Niu, D.K.; Jiang, L. Can ENCODE tell us how much junk DNA we carry in our genome? Biochem. Biophys. Res. Commun. 2013, 430, 1340–1343. [Google Scholar] [CrossRef] [PubMed]
  11. Kiss, T. Biogenesis of small nuclear RNPs. J. Cell Sci. 2004, 117, 5949–5951. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Prensner, J.R.; Chinnaiyan, A.M. The emergence of lncRNAs in cancer biology. Cancer Discov. 2011, 1, 391–407. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Gomes, A.Q.; Nolasco, S.; Soares, H. Non-Coding RNAs: Multi-Tasking Molecules in the Cell. Int. J. Mol. Sci. 2013, 14, 16010–16039. [Google Scholar] [CrossRef] [PubMed]
  14. Chi, Y.; Wang, D.; Yu, W.; Yang, J. Long Non-Coding RNA in the Pathogenesis of Cancers. Cells 2019, 8, 1015. [Google Scholar] [CrossRef] [Green Version]
  15. Schmitt, A.M.; Chang, H.Y. Long Noncoding RNAs in Cancer Pathways. Cancer Cell 2016, 29, 452–463. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Yao, R.W.; Wang, Y.; Chen, L.L. Cellular functions of long noncoding RNAs. Nat. Cell Biol. 2019, 21, 542–551. [Google Scholar] [CrossRef]
  17. Derrien, T.; Johnson, R.; Bussotti, G.; Tanzer, A.; Djebali, S.; Tilgner, H.; Guernec, G.; Martín, D.; Merkel, A.; Knowles, D.G.; et al. The GENCODE v7 catalog of human long noncoding RNAs: Analysis of their gene structure, evolution, and expression. Genome Res. 2012, 22, 1775–1789. [Google Scholar] [CrossRef] [Green Version]
  18. Johnsson, P.; Lipovich, L.; Grandér, D.; Morris, K.V. Evolutionary conservation of long non-coding RNAs; sequence, structure, function. Biochim. Biophys. Acta 2014, 1840, 1063–1071. [Google Scholar] [CrossRef] [Green Version]
  19. Hubbard, T.; Barker, D.; Birney, E.; Cameron, G.; Chen, Y.; Clark, L.; Cox, T.; Cuff, J.; Curwen, V.; Down, T.; et al. The Ensembl genome database project. Nucl. Acid Res. 2002, 30, 38–41. [Google Scholar] [CrossRef] [Green Version]
  20. Morris, K.V.; Mattick, J.S. The rise of regulatory RNA. Nat. Rev. Genet. 2014, 15, 423–437. [Google Scholar] [CrossRef] [Green Version]
  21. Iyer, M.K.; Niknafs, Y.S.; Malik, R. The landscape of long noncoding RNAs in the human transcriptome. Nat. Genet. 2015, 47, 199–208. [Google Scholar] [CrossRef] [PubMed]
  22. Ji, Z.; Song, R.; Regev, A.; Struhl, K. Many lncRNAs, 5′UTRs and pseudogenes are translated and some are likely to Express functional proteins. eLife 2015, 4, e08890. [Google Scholar] [CrossRef] [PubMed]
  23. Yan, L.; Yang, M.; Guo, H.; Yang, L.; Wu, J.; Li, R.; Liu, P.; Lian, Y.; Zheng, X.; Yan, J.; et al. Single-cell RNA-Seq profliling of human preimplantation embryos and embryonic stem cells. Nat. Struct. Mol. Biol. 2013, 20, 1131–1139. [Google Scholar] [CrossRef] [PubMed]
  24. Pirogov, S.A.; Gvozdev, V.A.; Klenov, M.S. Long Noncoding RNAs and Stress Response in the Nucleolus. Cells 2019, 8, 668. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Guglas, K.; Kolenda, T.; Teresiak, A.; Kopczyńska, M.; Łasińska, I.; Mackiewicz, J.; Mackiewicz, A.; Lamperska, K. lncRNA Expression after Irradiation and Chemoexposure of HNSCC Cell Lines. Noncoding RNA 2018, 4, 33. [Google Scholar] [CrossRef] [Green Version]
  26. Puvvula, P.K. LncRNAs Regulatory Networks in Cellular Senescence. Int. J. Mol. Sci. 2019, 20, 2615. [Google Scholar] [CrossRef] [Green Version]
  27. Dahariya, S.; Paddibhatla, I.; Kumar, S.; Raghuwanshi, S.; Pallepati, A.; Gutti, R.K. Long non-coding RNA: Classification, biogenesis and functions in blood cells. Mol. Immunol. 2019, 112, 82–92. [Google Scholar] [CrossRef]
  28. Liu, S.; Liu, X.; Li, J.; Zhou, H.; Carr, M.J.; Zhang, Z.; Shi, W. Long noncoding RNAs: Novel regulators of virus-host interactions. Rev. Med. Virol. 2019, 29, e2046. [Google Scholar] [CrossRef]
  29. Lin, Y.H.; Wu, M.H.; Yeh, C.T.; Lin, K.H. Long Non-Coding RNAs as Mediators of Tumor Microenvironment and Liver Cancer Cell Communication. Int. J. Mol. Sci. 2018, 19, 3742. [Google Scholar] [CrossRef] [Green Version]
  30. Kung, J.T.; Colognori, D.; Lee, J.T. Long noncoding RNAs: Past, present, and future. Genetics 2013, 193, 651–669. [Google Scholar] [CrossRef] [Green Version]
  31. Kolenda, T.; Guglas, T.; Ryś, M.; Bogaczyńska, M.; Teresiak, A.; Bliźniak, R.; Łasińska, I.; Mackiewicz, J.; Lamperska, K. Biological role of long non-coding RNA in head and neck cancers. Rep. Pract. Oncol. Radiother. 2017, 22, 378–388. [Google Scholar] [CrossRef]
  32. Kolenda, T.; Guglas, K.; Kopczyńska, M.; Teresiak, A.; Bliźniak, R.; Mackiewicz, A.; Lamperska, K.; Mackiewicz, J. Oncogenic Role of ZFAS1 lncRNA in Head and Neck Squamous Cell Carcinomas. Cells 2019, 8, 366. [Google Scholar] [CrossRef] [Green Version]
  33. Kolenda, T.; Kopczyńska, M.; Guglas, K.; Teresiak, A.; Bliźniak, R.; Łasińska, I.; Mackiewicz, J.; Lamperska, K. EGOT lncRNA in head and neck squamous cell carcinomas. Pol. J. Pathol. 2018, 69, 356–365. [Google Scholar] [CrossRef] [PubMed]
  34. Jin, J.; Wu, X.; Yin, J.; Li, M.; Shen, J.; Li, J.; Zhao, Y.; Zhao, Q.; Wu, J.; Wen, Q.; et al. Identification of Genetic Mutations in Cancer: Challenge and Opportunity in the New Era of Targeted Therapy. Front. Oncol. 2019, 9, 263. [Google Scholar] [CrossRef] [PubMed]
  35. National Comprehensive Cancer Network Website. Available online: https://www.nccn.org (accessed on 20 June 2020).
  36. Malone, E.; Siu, L.L. Precision Medicine in Head and Neck Cancer: Myth or Reality? Clin. Med. Insights Oncol. 2018, 12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Goossens, N.; Nakagawa, S.; Sun, X.; Hoshida, Y. Cancer biomarker discovery and validation. Transl. Cancer Res. 2015, 4, 256–269. [Google Scholar]
  38. Lianidou, E.; Pantel, K. Liquid biopsies. Genes Chromosomes Cancer 2019, 58, 219–232. [Google Scholar] [CrossRef]
  39. Sole, C.; Arnaiz, E.; Manterola, L.; Otaegui, D.; Lawrie, C.H. The circulating transcriptome as a source of cancer liquid biopsy biomarkers. Semin. Cancer Biol. 2019, 58, 100–108. [Google Scholar] [CrossRef]
  40. Eisenauer, E.A.; Therasse, P.; Bogaerts, J.; Schwartz, L.; Sargent, D.; Ford, R.; Dancey, J.; Arbuck, S.; Gwyther, S.; Mooney, M.; et al. New response evaluation criteria in solid tumors: Revised RECIST guideline (version 1.1). Eur. J. Cancer Suppl. 2009, 45, 228–247. [Google Scholar] [CrossRef]
  41. Lippi, G. Systematic Assessment of the Hemolysis Index: Pros and Cons. Adv. Clin. Chem. 2015, 71, 157–170. [Google Scholar]
  42. Kolenda, T.; Ryś, M.; Guglas, K.; Teresiak, A.; Bliźniak, R.; Mackiewicz, J.; Lamperska, K. Quantification of long non-coding RNAs using qRT-PCR: Comparison of different cDNA synthesis methods and RNA stability. Arch. Med. Sci. 2019, 15. [Google Scholar] [CrossRef] [Green Version]
  43. Thomas, G.R.; Shnayder, Y. Chapter 40—genomic evaluation of head and neck cancer A2—ginsburg, geoffrey S. In Essentials of Genomic and Personalized Medicine; Willard, H.F., Ed.; Academic Press: Cambridge, MA, USA, 2010; pp. 511–521. [Google Scholar]
  44. Kolenda, T.; Guglas, K.; Baranowski, D.; Sobocińska, J.; Kopczyńska, M.; Teresiak, A.; Bliźniak, R.; Lamperska, K. cfRNAs as biomarkers in oncology—Still experimental or applied tool for personalized medicine already? Rep. Pract. Oncol. Radiother. 2020, 25, 783–792. [Google Scholar] [CrossRef]
  45. Yao, Y.; Chen, X.; Lu, S.; Zhou, C.; Xu, G.; Yan, Z.; Yang, J.; Yu, T.; Chen, W.; Qian, Y.; et al. Circulating Long Noncoding RNAs as Biomarkers for Predicting Head and Neck Squamous Cell Carcinoma. Cell. Physiol. Biochem. 2018, 50, 1429–1440. [Google Scholar] [CrossRef]
  46. Lu, C.W.; Zhou, D.D.; Xie, T.; Hao, J.L.; Pant, O.P.; Lu, C.B.; Liu, X.F. HOXA11 antisense long noncoding RNA (HOXA11-AS): A promising lncRNA in human cancers. Cancer Med. 2018, 7, 3792–3799. [Google Scholar] [CrossRef] [PubMed]
  47. Xue, J.Y.; Huang, C.; Wang, W.; Li, H.B.; Sun, M.; Xie, M. HOXA11-AS: A novel regulator in human cancer proliferation and metastasis. Onco Targets Ther. 2018, 11, 4387–4393. [Google Scholar] [CrossRef] [Green Version]
  48. Thin, K.Z.; Tu, J.C.; Raveendran, S. Long non-coding SNHG1 in cancer. Clin. Chim. Acta 2019, 494, 38–47. [Google Scholar] [CrossRef] [PubMed]
  49. Lin, S.X.; Jiang, H.; Xiang, G.Z.; Zhang, W.R.; Weng, Y.H.; Qiu, F.D.; Wu, J.; Wang, H.G. Up-regulation of long non-coding RNA SNHG1 contributes to proliferation and metastasis in laryngeal squamous cell carcinoma. Eur. Rev. Med. Pharmacol. Sci. 2018, 22, 1333–1341. [Google Scholar]
  50. Gao, L.; Cao, H.; Cheng, X. A positive feedback regulation between long noncoding RNA SNHG1 and YAP1 modulates growth and metastasis in laryngeal squamous cell carcinoma. Am. J. Cancer Res. 2018, 8, 1712–1724. [Google Scholar]
  51. Wang, L.; Su, K.; Wu, H.; Li, J.; Song, D. LncRNA SNHG3 regulates laryngeal carcinoma proliferation and migration by modulating the miR-384/WEE1 axis. Life Sci. 2019, 232, 116597. [Google Scholar] [CrossRef]
  52. Tang, Y.; Wu, L.; Zhao, M.; Zhao, G.; Mao, S.; Wang, L.; Liu, S.; Wang, X. LncRNA SNHG4 promotes the proliferation, migration, invasiveness, and epithelial-mesenchymal transition of lung cancer cells by regulating miR-98-5p. Biochem. Cell Biol. 2019, 97, 767–776. [Google Scholar] [CrossRef] [PubMed]
  53. Liu, D.; Wang, Y.; Zhao, Y.; Gu, X. LncRNA SNHG5 promotes nasopharyngeal carcinoma progression by regulating miR-1179/HMGB3 axis. BMC Cancer 2020, 20, 178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Zhang, S.; Qiu, D.; Xie, X.; Shen, Y. Clinicopathological and prognostic value of SNHG6 in cancers: A systematic review and a meta-analysis. BMC Cancer 2020, 20, 343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Zhao, Y.; Wang, J.; Ma, K. Knockdown of lncRNA SNHG6 Inhibites the Proliferation and Epithelial Mesenchymal Transition in Tongue Cancer Cells. Xi Bao Yu Fen Zi Mian Yi Xue Za Zhi 2018, 34, 806–811. [Google Scholar] [PubMed]
  56. Guo, T.; Zambo, K.D.A.; Zamuner, F.T.; Ou, T.; Hopkins, C.; Kelley, D.Z.; Wulf, H.A.; Winkler, E.; Erbe, R.; Danilova, L.V.; et al. Chromatin structure regulates cancer-specific alternative splicing events in primary HPV-related oropharyngeal squamous cell carcinoma. Epigenetics 2020, 15, 959–971. [Google Scholar] [CrossRef]
  57. Martinez, B.V.; Dhahbi, J.M.; Lopez, Y.O.M.; Lamperska, K.; Golusinski, P.; Luczewski, L.; Kolenda, T.; Atamna, H.; Spindler, S.R.; Golusinski, W.; et al. Circulating small non-coding RNA signature in head and neck squamous cell carcinoma. Oncotarget 2015, 6, 19246–19263. [Google Scholar] [CrossRef]
  58. Kolenda, T.; Rutkowski, P.; Michalak, M.; Kozak, K.; Guglas, K.; Ryś, M.; Galus, Ł.; Woźniak, S.; Ługowska, I.; Gos, A.; et al. Plasma lncRNA expression profile as a prognostic tool in BRAF-mutant metastatic melanoma patients treated with BRAF inhibitor. Oncotarget 2019, 10, 3879–3893. [Google Scholar] [CrossRef] [Green Version]
  59. Guglas, K.; Kolenda, T.; Stasiak, M.; Kopczyńska, M.; Teresiak, A.; Ibbs, M.; Bliźniak, R.; Lamperska, K. YRNAs: New Insights and Potential Novel Approach in Head and Neck Squamous Cell Carcinoma. Cells 2020, 9, 1281. [Google Scholar] [CrossRef]
  60. Guglas, K.; Kołodziejczak, I.; Kolenda, T.; Kopczyńska, M.; Teresiak, A.; Sobocińska, J.; Bliźniak, R.; Lamperska, K. YRNAs and YRNA-Derived Fragments as New Players in Cancer Research and Their Potential Role in Diagnostics. Int. J. Mol. Sci. 2020, 21, 5682. [Google Scholar] [CrossRef]
  61. Drewell, R.A.; Arney, K.L.; Arima, T.; Barton, S.C.; Brenton, J.D.; Surani, M.A. Novel conserved elements upstream of the H19 gene are transcribed and act as mesodermal enhancers. Development 2002, 129, 1205–1213. [Google Scholar]
  62. Song, C.; Song, C.; Chen, K.; Zhang, X. Inhibition of long non-coding RNA IGF2AS protects apoptosis and neuronal loss in anesthetic-damaged mouse neural stem cell derived neurons. Biomed. Pharmacother. 2017, 85, 2018–2024. [Google Scholar] [CrossRef]
  63. Zhang, X.; Zhang, X.; Hu, R.; Hao, L. Retraction: Prognostic implication and functional role of long noncoding RNA IGF2AS in human non-small cell lung cancer. J. Cell. Biochem. 2019, 20, 12067. [Google Scholar]
  64. Cohen, N.; Fedewa, S.; Chen, A.Y. Epidemiology and demographics of the head and neck cancer population. Oral Maxillofac. Surg. Clin. N. Am. 2018, 30, 381–395. [Google Scholar] [CrossRef] [PubMed]
  65. Fayda, M.; Isin, M.; Tambas, M. Do circulating long non-coding RNAs (lncRNAs) (LincRNA-p21, GAS 5, HOTAIR) predict the treatment response in patients with head and neck cancer treated with chemoradiotherapy? Tumour Biol. 2016, 37, 3969–3978. [Google Scholar] [CrossRef] [PubMed]
  66. Chen, H.; Xin, Y.; Zhou, L. Cisplatin and paclitaxel target significant long noncoding RNAs in laryngeal squamous cell carcinoma. Med. Oncol. 2014, 31, 246. [Google Scholar] [CrossRef]
  67. Ren, S.; Li, G.; Liu, C. Next generation deep sequencing identified a novel lncRNA n375709 associated with paclitaxel resistance in nasopharyngeal carcinoma. Oncol. Rep. 2016, 36, 1861–1867. [Google Scholar] [CrossRef]
  68. Li, L.; Gu, M.; You, B. Long non-coding RNA ROR promotes proliferation, migration and chemoresistance of nasopharyngeal carcinoma. Cancer Sci. 2016, 107, 1215–1222. [Google Scholar] [CrossRef]
  69. Tasheva, E.S.; Roufa, D.S. Regulation of human RPS14 transcription by intronic antisense RNAs and ribosomal protein S14. Genes Dev. 1995, 9, 304–316. [Google Scholar] [CrossRef] [Green Version]
  70. Tang, Z.; Wei, G.; Zhang, L.; Xu, Z. Signature microRNAs and long noncoding RNAs in laryngeal cancer recurrence identified using a competing endogenous RNA network. Mol. Med. Rep. 2019, 19, 4806–4818. [Google Scholar] [CrossRef] [Green Version]
  71. Liang, K.; Yang, Y.; Zha, D.; Yue, B.; Qiu, J.; Zhang, C. Overexpression of lncRNA snaR is correlated with progression and predicts poor survival of laryngeal squamous cell carcinoma. J. Cell. Biochem. 2018. [Google Scholar] [CrossRef]
Jpm 10 00131 g001 550
Figure 1. Significantly higher (p < 0.05) expression levels of 34 plasma lncRNAs between healthy volunteers (black) and head and neck squamous cell carcinoma (HNSCC) patients (red). Squares and circles represent means of expression and lines indicate 95% CI; * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 1. Significantly higher (p < 0.05) expression levels of 34 plasma lncRNAs between healthy volunteers (black) and head and neck squamous cell carcinoma (HNSCC) patients (red). Squares and circles represent means of expression and lines indicate 95% CI; * p < 0.05, ** p < 0.01, *** p < 0.001.
Jpm 10 00131 g001
Jpm 10 00131 g002 550
Figure 2. Significantly higher (p < 0.05) expression levels of 41 plasma lncRNAs between healthy volunteers (black) and recurrent and/or metastatic HNSCC patients (red). Squares and circles represent means of expression and lines indicate 95% CI; * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 2. Significantly higher (p < 0.05) expression levels of 41 plasma lncRNAs between healthy volunteers (black) and recurrent and/or metastatic HNSCC patients (red). Squares and circles represent means of expression and lines indicate 95% CI; * p < 0.05, ** p < 0.01, *** p < 0.001.
Jpm 10 00131 g002
Jpm 10 00131 g003 550
Figure 3. Significantly changed (p < 0.05) expression levels of 47 plasma lncRNAs between healthy volunteers (black) and locally advanced HNSCC patients (red). Squares and circles represent the mean of expression and lines indicate 95% CI; * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 3. Significantly changed (p < 0.05) expression levels of 47 plasma lncRNAs between healthy volunteers (black) and locally advanced HNSCC patients (red). Squares and circles represent the mean of expression and lines indicate 95% CI; * p < 0.05, ** p < 0.01, *** p < 0.001.
Jpm 10 00131 g003
Jpm 10 00131 g004 550
Figure 4. Changes in plasma lncRNA expression in HNSCC patients depending on the tumor localization in the pharynx, oral cavity, larynx and salivary gland sites. Only statistically significant changes (p < 0.05) in lncRNAs between groups are presented.
Figure 4. Changes in plasma lncRNA expression in HNSCC patients depending on the tumor localization in the pharynx, oral cavity, larynx and salivary gland sites. Only statistically significant changes (p < 0.05) in lncRNAs between groups are presented.
Jpm 10 00131 g004
Jpm 10 00131 g005 550
Figure 5. lncRNA expressions in HNSCC patients receiving palliative chemotherapy according to response to treatment (disease control rate vs. progression). Data present median expression of lncRNA with confidence interval (CI); p < 0.05.
Figure 5. lncRNA expressions in HNSCC patients receiving palliative chemotherapy according to response to treatment (disease control rate vs. progression). Data present median expression of lncRNA with confidence interval (CI); p < 0.05.
Jpm 10 00131 g005
Jpm 10 00131 g006 550
Figure 6. lncRNA expression in palliative group with (A) disease control rate (DCR) and (B) progressive disease (PD) response, and (C) in a locally advanced group, depending on two time points—before chemotherapy (P1) and after three cycles of chemotherapy (P2). Data present median expression of lncRNA with confidence interval (CI); p < 0.05.
Figure 6. lncRNA expression in palliative group with (A) disease control rate (DCR) and (B) progressive disease (PD) response, and (C) in a locally advanced group, depending on two time points—before chemotherapy (P1) and after three cycles of chemotherapy (P2). Data present median expression of lncRNA with confidence interval (CI); p < 0.05.
Jpm 10 00131 g006
Jpm 10 00131 g007 550
Figure 7. Progression free survival probability (A) and overall survival probability (B) of HNSCC patients receiving palliative chemotherapy depending on expression levels of Alpha 250, Emx2os, snaR and SNHG1 lncRNAs. Patients were divided into two subgroups according to high (red)/low (blue) expression levels of lncRNA, Kaplan–Meier analysis, and p < 0.05 considered as statistically significant.
Figure 7. Progression free survival probability (A) and overall survival probability (B) of HNSCC patients receiving palliative chemotherapy depending on expression levels of Alpha 250, Emx2os, snaR and SNHG1 lncRNAs. Patients were divided into two subgroups according to high (red)/low (blue) expression levels of lncRNA, Kaplan–Meier analysis, and p < 0.05 considered as statistically significant.
Jpm 10 00131 g007
Table
Table 1. Univariate Cox regression model for PFS and OS. PFS—progression free survival; OS—overall survival; HR—hazard ratio; CI—confidence interval; BMI—body mass index; NLR—neutrophil to lymphocyte ratio.
Table 1. Univariate Cox regression model for PFS and OS. PFS—progression free survival; OS—overall survival; HR—hazard ratio; CI—confidence interval; BMI—body mass index; NLR—neutrophil to lymphocyte ratio.
ParameterGroup StratificationPFSOS
p-ValHR95% CIp-ValHR95% CI
Age<61 vs. >610.820.800.16–4.20.300.410.06–2.61
NLR<5 vs. >50.540.600.12–3.060.900.510.19–6.25
BMI<18.49 vs. >250.891.100.22–5.750.462.150.28–16.8
Primary tumor localizationPharynx vs. oral cavity 0.400.670.01–30.10.440.230.0–10.48
Primary tumor localizationPharynx vs.
larynx		0.400.670.01–30.10.440.230.0–10.48
Primary tumor localizationPharynx vs. salivary gland0.400.670.01–30.10.440.230.0–10.48
Primary tumor localizationOral cavity vs. salivary gland0.461.800.36–9.270.440.220.14–8.60
Primary tumor localizationLarynx vs. salivary gland 0.901.100.21–5.760.901.090.14–8.60
lncRNA Alpha 250High vs. low expression0.0390.100.01–0.890.0420.200.01–0.88
lncRNA Emx2osHigh vs. low expression0.130.260.04–1.54---
lncRNA snaRHigh vs. low expression---0.203.00.45–20.50
lncRNA SNHG1High vs. low expression---0.106.90.66–74.12
Table
Table 2. Multivariate Cox regression model for PFS and OS. PFS—progression free survival; OS—overall survival; HR—hazard ratio; CI—confidence interval; BMI—body mass index; NLR—neutrophil to lymphocyte ratio.
Table 2. Multivariate Cox regression model for PFS and OS. PFS—progression free survival; OS—overall survival; HR—hazard ratio; CI—confidence interval; BMI—body mass index; NLR—neutrophil to lymphocyte ratio.
ParameterGroup StratificationPFSOS
p-ValHR95%CIp-ValHR95%CI
Age<61 vs. >610.900.840.02–30.310.600.480.03–9.05
NLR<5 vs. >50.900.860.02–33.560.900.800.06–12.33
BMI<18.49 vs. >250.901.090.22–47.520.901.010.02–43.42
Primary tumor localizationPharynx vs. oral cavity0.900.910–285.890.800.230–135.53
Primary tumor localizationPharynx vs.
larynx		0.900.910–285.890.800.690–135.53
Primary tumor localizationPharynx vs. salivary gland0.900.910–285.890.800.600–135.53
Primary tumor localizationOral cavity vs. salivary gland0.801.480.04–55.740.900.800.04–19.84
Primary tumor localizationLarynx vs. salivary gland 0.900.930.2–37.210.701.730.07–41.21
lncRNA Alpha 250High vs. low expression0.0390.200.01–0.870.0430.300.02–0.89
lncRNA Emx2osHigh vs. low expression0.800.670.02–24.62---
lncRNA snaRHigh vs. low expression---0.403.070.17–55.67
lncRNA SNHG1high vs. low expression---0.502.770.14–54.37
Table
Table 3. Expression levels of Alpha 250, Emx2os, snaR and SNHG1 lncRNAs depending on clinical and laboratory parameters in HNSCC patients; N—BMI between 18.5–24,8; PN—BMI under 18.5, BMI—body mass index, NLR—neutrophil to lymphocyte ratio, PLR—platelets to lymphocyte ratio. Mann–Whitney test and one-way ANOVA; p < 0.05 considered as statistically significant.
Table 3. Expression levels of Alpha 250, Emx2os, snaR and SNHG1 lncRNAs depending on clinical and laboratory parameters in HNSCC patients; N—BMI between 18.5–24,8; PN—BMI under 18.5, BMI—body mass index, NLR—neutrophil to lymphocyte ratio, PLR—platelets to lymphocyte ratio. Mann–Whitney test and one-way ANOVA; p < 0.05 considered as statistically significant.
ParameterAlpha 250Emx2ossnaRSNHG1
MedianLower QuartileUpper Quartilep-ValMedianLower QuartileUpper Quartilep-ValMedianLower QuartileUpper Quartilep-ValMedianLower QuartileUpper Quartilep-Val
Age
>610.00350.00140.08260.65250.24180.02570.59460.652512.10997.02426.74250.01830.50920.14191.39470.3046
<610.06790.0071.22190.08040.04620.38032.5140.73545.90420.28330.0371.1563
Sex
Female73.8260.0932155.0120.79740.15540.035240.3590.907112.2610.4948930.91140.69160.3480.022443.12570.7616
Male0.00950.00260.11510.11810.04580.58644.73721.282917.18130.41330.10421.1313
Nicotinism
Yes0.04460.00720.63220.28260.4840.04731.12990.729310.30911.842631.64330.47750.65520.37281.160.1216
No0.00650.00240.21380.10310.03110.26352.71320.388610.68720.16720.07150.6246
Alcoholism
Yes0.14860.00890.92240.81070.24180.05182.32020.67185.58671.815225.33110.64510.37280.06411.15510.5932
No0.00350.0020.03960.11060.03480.51875.54040.640612.10990.46010.08811.5281
BMI
N0.0240.00290.17240.56190.21320.09070.44040.92553.9342.081312.2610.7790.47560.1391.46590.779
PN0.00990.00240.5320.05950.02840.89537.09990.84422.50370.32530.07391.0851
NLR 1
>50.12540.00550.46690.28620.24180.05180.39170.27014.49060.660413.92620.87080.35240.0450.95960.7863
<50.00670.0010.07720.10220.02470.59265.54041.549129.78740.42980.10752.0887
PLR 1
>2150.03990.00320.25210.46020.26350.05550.75890.18862.41560.833811.36050.37750.21840.05650.75850.3219
<2150.00980.00180.20110.08040.02980.58099.11352.51435.87730.52570.1482.425
Stage
IVA0.00260.00050.00350.2150.10310.02570.58640.79160.9760.28239.11350.13780.08230.02130.48320.0401
IVC0.06790.00930.33680.16290.04580.57636.2682.117923.69760.46010.16721.1551
Localization
Pharynx0.20110.10210.48690.2450.24650.0760.5570.66523.9342.360217.0830.40770.39940.23910.77730.5987
Larynx0.02930.00370.21070.05180.03410.29383.81091.19669.43150.27780.07730.8421
Salivary gland0.00010.00010.00260.30610.017711.598169.190423.697681.00842.17470.427533.5545
Oral cavity0.0250.00820.24940.19080.07040.592.71321.306410.53610.29720.07390.7585

Share and Cite

MDPI and ACS Style

Łasińska, I.; Kolenda, T.; Guglas, K.; Kopczyńska, M.; Sobocińska, J.; Teresiak, A.; Strzelecki, N.O.; Lamperska, K.; Mackiewicz, A.; Mackiewicz, J. Liquid lncRNA Biopsy for the Evaluation of Locally Advanced and Metastatic Squamous Cell Carcinomas of the Head and Neck. J. Pers. Med. 2020, 10, 131. https://doi.org/10.3390/jpm10030131

AMA Style

Łasińska I, Kolenda T, Guglas K, Kopczyńska M, Sobocińska J, Teresiak A, Strzelecki NO, Lamperska K, Mackiewicz A, Mackiewicz J. Liquid lncRNA Biopsy for the Evaluation of Locally Advanced and Metastatic Squamous Cell Carcinomas of the Head and Neck. Journal of Personalized Medicine. 2020; 10(3):131. https://doi.org/10.3390/jpm10030131

Chicago/Turabian Style

Łasińska, Izabela, Tomasz Kolenda, Kacper Guglas, Magda Kopczyńska, Joanna Sobocińska, Anna Teresiak, Norbert Oksza Strzelecki, Katarzyna Lamperska, Andrzej Mackiewicz, and Jacek Mackiewicz. 2020. "Liquid lncRNA Biopsy for the Evaluation of Locally Advanced and Metastatic Squamous Cell Carcinomas of the Head and Neck" Journal of Personalized Medicine 10, no. 3: 131. https://doi.org/10.3390/jpm10030131

APA Style

Łasińska, I., Kolenda, T., Guglas, K., Kopczyńska, M., Sobocińska, J., Teresiak, A., Strzelecki, N. O., Lamperska, K., Mackiewicz, A., & Mackiewicz, J. (2020). Liquid lncRNA Biopsy for the Evaluation of Locally Advanced and Metastatic Squamous Cell Carcinomas of the Head and Neck. Journal of Personalized Medicine, 10(3), 131. https://doi.org/10.3390/jpm10030131

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop