• Open Access
    Original Article

    Liquid biopsy: an examination of platelet RNA obtained from head and neck squamous cell carcinoma patients for predictive molecular tumor markers

    Lisa T. Huber 1† , 
    Johann M. Kraus 2† , 
    Jasmin Ezić 1† , 
    Amin Wanli 1 , 
    Marco Groth 3 , 
    Simon Laban 1 , 
    Thomas K. Hoffmann 1 , 
    Barbara Wollenberg 4 , 
    Hans A. Kestler 2* , 
    Cornelia Brunner 1*

    Explor Target Antitumor Ther. 2023;4:422–446 DOI: https://doi.org/10.37349/etat.2023.00143

    Received: November 16, 2022 Accepted: March 01, 2023 Published: June 30, 2023

    Academic Editor: Giulia Martini, University of Campania “Luigi Vanvitelli”, Italy

    This article belongs to the special issue The Implementation of Liquid Biopsy in Clinical Practice for Different Solid Tumor

    Abstract

    Aim:

    Recently, a tumor cell-platelet interaction was identified in different tumor entities, resulting in a transfer of tumor-derived RNA into platelets, named further “tumor-educated platelets (TEP)”. The present pilot study aims to investigate whether such a tumor-platelet transfer of RNA occurs also in patients suffering from head and neck squamous cell carcinoma (HNSCC).

    Methods:

    Sequencing analysis of RNA derived from platelets of tumor patients (TPs) and healthy donors (HDs) were performed. Subsequently, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was used for verification of differentially expressed genes in platelets from TPs and HDs in a second cohort of patients and HDs. Data were analyzed by applying bioinformatic tools.

    Results:

    Sequencing of RNA derived from the tumor as well as from platelets of TPs and HDs revealed 426 significantly differentially existing RNA, at which 406 RNA were more and 20 RNA less abundant in platelets from TPs in comparison to that of HDs. In TPs’ platelets, abundantly existing RNA coding for 49 genes were detected, characteristically expressed in epithelial cells and RNA, the products of which are involved in tumor progression. Applying bioinformatic tools and verification on a second TP/HD cohort, collagen type I alpha 1 chain (COL1A1) and zinc finger protein 750 (ZNF750) were identified as the strongest potentially platelet-RNA-sequencing (RNA-seq)-based biomarkers for HNSCC.

    Conclusions:

    These results indicate a transfer of tumor-derived messenger RNA (mRNA) into platelets of HNSCC patients. Therefore, analyses of a patient’s platelet RNA could be an efficient option for liquid biopsy in order to diagnose HNSCC or to monitor tumorigenesis as well as therapeutic responses at any time and in real time.

    Keywords

    Tumor educated platelets, liquid biopsy, head and neck

    Introduction

    Head and neck squamous cell carcinoma (HNSCC) describes tumors that originate from the epithelium of the oral cavity, pharynx, and larynx. They represent the sixth most prevalent cancer type globally with 932,000 cases and a mortality of 467,000 deaths every year as of 2020 [1]. Risk factors primarily include smoking and excessive alcohol consumption, as well as aging, diet, oral hygiene, and family history. Additionally, infection with human papilloma virus (HPV), a circular DNA virus that targets the epithelial cells of the mucosa, plays an important role in the development of HNSCC. In general, treatment options include surgery, radio-, chemo- and immunotherapy [2]. Besides this, novel therapeutic approaches based on molecular tumor markers were developed like the targeting of epidermal growth factor receptor (EGFR), a highly upregulated receptor in HNSCC associated with poor prognosis. Additionally, as an immunotherapeutic approach, the checkpoint inhibitor pembrolizumab has been approved by the Food and Drug Administration (FDA) in 2019 as a first-line treatment for metastatic or recurrent HNSCC [3]. Despite recent advances in HNSCC therapeutic options, the mortality rate is still around 50%, which indicates the need for more effective therapies and diagnostic as well as disease monitoring tools. For this, reliable biomarkers should be defined that can predict the onset and progression of HNSCC or even the therapeutic outcome. Moreover, based on biomarker expression personalized therapeutic options could be developed and applied. However, presumably because of the high heterogeneity in HNSCC, in contrast to other tumor entities, a robust biomarker panel for HNSCC has not been defined so far, as applicable in clinical settings [4].

    There are different ways to obtain information about the genetic changes of the tumor. Although regarded as the gold standard for tumor diagnosis, tissue biopsy suffers from several limitations. Firstly, tissue biopsy is invasive in nature and can potentially lead to complications, causing limitations in the reproducibility and repeatability of the diagnosis. On top of that, tumor heterogeneity impedes the ability to obtain spatially and temporally accurate information about the tumors [5]. For this reason, recent efforts have been made to find new methods for cancer diagnosis and screening options.

    Liquid biopsy describes the analysis of biomarkers released by tumors into body fluids, such as blood, urine, or saliva. These biomarkers could be cells, nucleic acids, or proteins. Liquid biopsy offers a minimally invasive procedure, which is more easily applicable and could be repeated more frequently, offering a more real-time and dynamic snapshot of the tumor, which combats the spatial and temporal limitation of tissue biopsy [6]. This enables clinicians not only to utilize liquid biopsy as an effective diagnostic tool but also as an instrument for therapy and tumor progression monitoring. The advantage of analyzing genetic markers also offers a potential approach for personalized treatment [7].

    Information about potentially existing tumors could be obtained from several sources through liquid biopsy, like circulating tumor cells (CTCs), circulating cell-free tumor DNA (ctDNA), and exosomes [8–10]. Their analyses were suggested as sensitive and specific for cancer prediction, diagnosis, prognosis, and monitoring over time [11–13] including HNSCC [14–18]. Thus, a recent study on 21 patients with locally advanced oropharyngeal carcinoma treated with neoadjuvant chemotherapy suggested that the variation of CTC numbers might reflect the effectiveness of treatment [19]. Also, the potential of ctDNA as a biomarker for detecting minimal residual disease and recurrence in HNSCC was demonstrated recently [20]. Additionally, the analyses of circulating exosomes emerged as a potential tool for HNSCC diagnosis [21] and several studies correlated proteins identified in patients’ exosomes with the activity of tumor disease and therapy responses [21, 22].

    However, with regard to head and neck tumors, all these types of liquid biopsies are far from routine clinical use because, despite many advantages over standard diagnostics, they have some limitations. Over the last years, another form of liquid biopsy has received considerable attention—the analyses of so-called “tumor-educated platelets (TEP)”. Platelets interact with tumor cells leading to platelet degranulation, known as tumor cell-induced platelet aggregation. The release of their storage products supports tumor growth, angiogenesis, and tumor spread [23].

    Cancer-platelet interaction is a two-way communication path. Not only does platelet activation promote tumor progression, but also tumor cells affect platelets, specifically, their RNA profile. The idea of platelets having a distinct RNA profile in cancer patients compared with healthy controls was first introduced in 2010 [24]. In this study, Calverley et al. [24] demonstrated that platelets of metastatic lung cancer patients exhibit 200 differentially expressed genes compared to healthy controls. In 2011, Nilsson et al. [25] described the transfer of tumor cell RNA into platelets, leading to the formation of so-called TEP. Based on platelet RNA profiles of 283 patients, Best et al. [26] developed an algorithm that distinguished healthy and tumor patients (TPs) with an accuracy of 96%. In this study, six tumor entities (non-small cell lung cancer, colorectal carcinoma, breast carcinoma, hepatobiliary carcinoma, glioblastoma, and pancreatic carcinoma) have been differentiated.

    In recent years, great efforts have been made to define TEP-derived biomarkers for the diagnosis of a wide variety of tumors for different cancer entities. Thus, applying particle-swarm optimization-enhanced algorithms, Best et al. [27] defined a biomarker panel allowing the detection of early- and late-stage non-small cell lung cancer with an accuracy of 81% and 88%, respectively. Sheng et al. [28] identified a biomarker panel of 48 genes for the diagnostic of the same tumor entity with an accuracy of 89%, whereas Xing et al. [29] identified ITGA2B as a diagnostic and prognostic TEP-based marker for non-small cell lung cancer, just to mention some of the TEP-RNA-sequencing (RNA-seq) work on this cancer entity. Additionally, TEP-derived RNA panels, determined again by swarm intelligence, were applied to identify glioblastoma patients with an accuracy of 95% [30]. TEP-RNA-seq-derived biomarkers were also identified for sarcoma [31], colorectal [32], and liver cancer [33–36]. Moreover, various clinical studies have confirmed the analysis of TEP as a robust form of liquid biopsy and other clinical trials are currently recruiting. The study described here was initiated because of the advantages of the analysis of TEP compared to other forms of liquid biopsy and the lack of data on TEP in HNSCC.

    The present work represents a pilot study with the aim to investigate whether the analysis of platelet RNA is a diagnostic option also for HNSCC. Since differentially existing messenger RNA (mRNA) in platelets of HNSCC patients were observed, the present study provides the basis for designing further prospective studies considering cohort size, robust inclusion criteria, in- or exclusion of certain tumor sites as well as HPV status.

    Materials and methods

    Patients

    The study was approved by the local ethics committee of Ulm University [90/15], and written signed informed consent was obtained from all the participants prior to conducting the study. The patient characteristics of cohorts 1 and 2 are shown in Tables 1 and 2. The patients were selected based on their diagnosis. In order to achieve the most robust results possible, which are not influenced by age, gender, or side diagnoses, these data were collected, but they did not influence inclusion in this study. The only inclusion criterion as a TP was the diagnosis of HNSCC; the inclusion criterion as a healthy donor (HD) was the confirmed exclusion of a tumor disease on the day of the blood donation. For the RNA-seq experimental study design, the suggestions of Lin et al. [37] for cohort estimation were followed. In order to determine sample sizes for the pilot data, a sample size estimation based on the assumption of d = 2, alpha = 0.05, beta = 0.8 using a Wilcoxon-test (software: R 3.2.2) [38] was performed.

    TP’s and HD’s characteristics, cohort 1

    IDAge (years)SexSide diagnosesNotesTumor localisationTNMHPV status
    TP_158MaleNone40 py, 1 beer/dayTonguepT2pN0 cM0Negative
    TP_218MaleNoneNoneTonguepT2pN0 cM0Negative
    TP_361MaleNone40 py, 2 beers/dayGlossotonsillarpT2pN0 cM0Negative
    TP_475MaleType 2 diabetes mellitusNoneTonguepT1pN0 cM0Negative
    TP_567MaleLoss of sensorineural hearing2 Beers/dayTonguepT3pN0 cM0Negative
    HD_159MaleNone42 py---
    HD_222MaleNoneN.d.---
    HD_363MaleProstate adenoma, rectal adenoma, hypertension, obstructive sleepN.d.---
    HD_474MaleApnea, loss of sensorineural hearingNone---
    HD_569MalePolyarthritis, hypertensionN.d.---
    Display full size

    ID: identification; py: pack years; TNM: tumor, nodes, metastasis; p: pathological; c: clinical; N.d.: no data; -: not applicable

    TP’s and HD’s characteristics, cohort 2

    IDAge (years)SexSide diagnosesNotesTumor localisationTNMHPV status
    TP_668MaleNone15 py, 1 bottle of wine/dayOral cavitypT1pN0 cM0Negative
    TP_771MaleNoneNoneTonguepT2pN0 cM0Negative
    TP_861MaleNone4 py, 2 beers/dayOral cavity/oropharynxpT2pN0 cM0HPV 33
    TP_957MaleNone30 py, 5 beers/dayOral cavity/tonsilpT3pN0 cM0Negative
    TP_1051FemaleNone1 pyTonguepT3pN1 cM0Negative
    TP_1159MaleBenign prostate hyperplasia45 pyTonguepT3pN0 cM0Negative
    HD_628MaleNoneNone---
    HD_728MaleNoneNone---
    HD_871MaleChronic otitis media, presbyakusis, hypertensionNone---
    HD_952FemaleHypothyroidismNone---
    HD_1070MaleNoneNone---
    HD_1174FemaleLoss of sensorineural hearingNone---
    HD_1222MaleChronic tonsillitis1.5 py, occasional alcohol consumption---
    Display full size

    -: not applicable

    Blood collection and platelets isolation

    Blood (50 mL) was taken using 10 mL S-monovets (Sarstedt, Nürnbrecht, Germany) each with 3.2% citrate, from HDs and HNSCC TPs before the onset of any therapeutic measures. “Healthy” is defined here as “no diagnosed tumor” at the time of blood collection. Additionally, tissue material from the primary tumor was collected. Donors’ monovets were centrifuged at 120 × g for 20 min. The exact number of platelets that were obtained has not been determined routinely. During the examination of the platelets under the microscope, random samples showed regular platelet counts, from 150,000/µL to 450,000/µL blood. The platelet-rich plasma was centrifuged at 360 × g for 20 min. The resulting pellet was resuspended in 400 μL to 500 μL RNAlaterTM (Qiagen, Hilden, Sweden). To check the purity of the platelets, 10 μL of platelet suspension was diluted with 90 μL of Tuerk solution (Merck, Darmstadt, Germany) and checked microscopically. The preparation was defined as pure if at most 5 nuclear-containing cells per 10 million platelets could be identified.

    RNA isolation and RNA-seq

    For isolation of total RNA 200 μL of platelet suspension was used and centrifuged at 8,000 × g for 5 min. Total RNA from platelets and tumor tissue (80 mg) was isolated according to the protocol of the mirVana™ miRNA Isolation Kit (Thermo Fisher Scientific, Waltham, US). The pieces of tissue were lysed for 4 min at 50 Hz in the TissueLyser LT (Qiagen, Hilden, Sweden) using 5 mm beads. RNA-seq was performed using the next-generation sequencing technology of Illumina (US) [39]. Library preparation was carried out with the TruSeq Stranded mRNA Kit (illumina®, San Diego, US) according to the manufacturer’s instructions [40]. The quality of the sequencing libraries was controlled using the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, US) and the Agilent DNA 7500 Kit (Agilent Technologies, Santa Clara, US). Sequencing was performed on a Illumina Hiseq2500 sequencing system (illumina, San Diego, US) in 51 bp, single-end sequencing, and high-output mode. Fifteen samples were sequenced on 4 lanes and reached 50 ± 5 mio reads per sample. The data were extracted with the bcl2fastq conversion software v1.8.4 program (illumina®, San Diego, US).

    RNA-seq data analysis

    Obtained reads were aligned with the human genome version hg 19 using the Bowtie2 software and the standard parameters [41]. Read counts were determined with high-throughput sequencing (HTSeq)-count based on annotation GRCh37.75 [42]. Data analysis was performed in R (4.0.2) using the integrated development environment (IDE) Rstudio (1.3.1056). Normalization and analysis of differentially existing RNAs were performed using DESeq2 (1.16.0.), and the shrinkage estimator used was “apeglm” [43, 44]. Principal component analysis (PCA) was computed within DESeq2. Clustering in the heatmaps was done using hierarchical clustering within the heatmap (1.0.12) package. For gene ontology, gene set enrichment analysis (GSEA) was performed [45–47]. For the computation of distance matrices, Euclidean distance was used. Genes with a false discovery rate (FDR) < 0.05 were considered differentially expressed.

    Quantitative reverse transcription-polymerase chain reaction

    Fifty ng RNA of the platelets was used in a one-step-polymerase chain reaction (PCR) reaction using QuantiNova RT-PCR kit (ID: 208352; Qiagen, Hilden, Germany) according to the protocol of the producer. Reactions were running at the LightCycler® 96 System (Roche, Basel, Schweiz). Primers were obtained from biomers.net GmbH (Ulm, Germany) and are presented in Table 3. Data were analyzed by the formula of relative quantification according to Pfaffl [48]. For statistical analyses, the Welch’s t-test or Mann-Whitney-U test was used, depending on whether the data were normally distributed or not.

    Primer characteristics

    Gene name/NCBI accession NREfficiency test mRNAAmplicon size (bp)Primer forwardPrimer reverse
    GAPDH/NM_001357943.22.00/universal RNA115GCTCTCTGCTCCTCCTGTTCACGACCAAATCCGTTGACTC
    UBC/NM_021009.72.00/universal RNA105CTGATCAGCAGAGGTTGATCTTTGACGGAGTACCAGGTGCAAG
    HPRT1/NM_000194.31.98/universal RNA94TGACACTGGCAAAACAATGCAGGTCCTTTTCACCAGCAAGCT
    PGK1/NM_000291.42.00/universal RNA181GGACAATGGAGCCAAGTCGGTGGGTTGGCACAGGCTTTCT
    COL1A1/NM_000088.41.96/universal RNA185TGGTTCGTGACCGTGACCTCCAGGTTGCAGCCTTGGTTGG
    ZNF750/NM_024702.32.00/UDSCC583GCACAGAATGCCTACCTGCCAGGAGGGTCTCCGTTCACAAC
    MAGEA6/NM_005363.52.00/universal RNA182TCGGGGATCCCAAGAAGCTGAATGCGAGGTCCTCCACTGA
    Ly6D/NM_003695.31.89/UDSCC577GCCAGCTCTCGCTTCTGCAATCCGCACAGTCCTTCTTCACC
    WNT5a/XM_011534089.22.00/universal RNA66AGCGGGTTCCTGAGTGAATCCCCTGACCCCTCTCTAGTT
    FAT1/NM_005245.42.00/universal RNA153ATCCTGCAACCGGCTCTCTCCGTGGTCATTCGTGTCGCTG
    KRT5/NM_000424.42.00/universal RNA91GCAGATCAAGACCCTCAACAATCCACTTGGTGTCCAGAACCT
    EGFR/NM_005228.52.00/universal RNA200TCCTCATTGCCCTCAACACAGTGGCAGGGTTGTTGCTGAACC
    PDPN/NM_198389.21.95/universal RNA176CTCAAACGTGGCCACCAGTCGCCTTCCCGACATTTTTCGCA
    CDH11/NM_001330576.22.00/universal RNA195ACACGGCCAATGGACCAAGATGCTACTCATGGGCGGGATG
    CA12/NM_206925.32.00/universal RNA173AGTTTTCCGAAACCCCGTGCGGATGATGCCTTGGGAGAAGGA
    RBP1/NM_001130992.32.00/universal RNA161GCTCCAGTCACTCCCCGAAACACCGTCCTGCACGATCTCT
    PLA2G4E/XM_011521238.2, XM_011521239.2, NM_001206670.1, XM_011521238.2

    2.00/universal RNA

    UDSCC5

    		76, 143, 193, 197GCCAGCTACATCACCGGTCT, TGTCAACTCCAGCTACCCGC, CAGAGCCCACAAACGGATGA, TGCAGAAGCGGAAGGTCGTTGGAGGACCAGTCAGGGTCA, TGGGGAAGGGGATGTTCTGC, CAGCCTGCCGGACATTTTTCA, ACCGGTGATGTAGCTGGCAC
    Display full size

    NR: number; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; UBC: ubiquitin C; HPRT1: hypoxanthine phosphoribosyltransferase 1; PGK1: phosphoglycerate kinase 1; COL1A1: collagen type I alpha 1 chain; ZNF750: zinc finger protein 750; MAGEA6: MAGE family member A6; Ly6D: lymphocyte antigen 6 family member D; WNT5a: Wnt family member 5a; FAT1: FAT atypical cadherin 1; KRT5: keratin 5; PDPN: podoplanin; CDH11: cadherin 11; CA12: carbonic anhydrase 12; RBP1: retinol binding protein 1; PLA2G4E: phospholipase A2 group IVE

    Results

    RNA-seq data gaining

    In order to investigate whether HNSCC-derived platelets harbor RNA originally synthesized in tumor cells, RNA-seq was performed on RNA isolated from platelets of HDs or patients suffering from HNSCC. Patients and HD characteristics are depicted in Tables 1 and 2.

    Additionally, RNA derived from the corresponding tumor was sequenced. Each sequenced sample showed about 50 ± 5 million reads (read length 51 bp) and 94% of the bases reached a Phred-scaled quality score Q of 30 to 40. Per sample, around 30 million reads could be mapped to genes in the human genome. More than 90% of these reads were localized in exons.

    To verify the extraction of platelet RNA, platelet markers like beta 2-microglobulin (B2M), pro-platelet basic protein (PPBP), thymosin beta 4 X-linked (TMSB4X) and platelet factor 4 (PF4) could be detected in both groups—platelets from TPs as well as from HDs (Table 4).

    Normalized counts of platelet marker genes

    HD/TPPPBPTMSB4XB2MPF4
    HD_123,44424,6067,82114,659
    HD_242,12278,49727,72021,999
    HD_3144,119107,55678,81273,119
    HD_423,27961,19511,42017,009
    HD_532,41436,76227,74814,760
    TP_168,510172,49253,73861,160
    TP_240,87936,38623,58420,828
    TP_344,63799,44551,79134,089
    TP_458,88744,98145,31440,587
    TP_5197,738157,380126,643111,286
    Display full size

    Patients platelet transcriptome is different from that of HDs

    To get insights into the data set obtained by RNA-seq, PCA was carried out (Figure 1). The PCA revealed sample separation, which was in accordance with the biological groups (tumor patient-tumor tissue, tumor patient-platelets; HD-platelets). Since the PCA was computed based on the whole transcriptome profile, different RNA profiles in each group can be assumed. Moreover, the samples of the platelets derived from TPs scatter on both axes of the diagram (Figure 1) indicating that platelets of TPs contain a different cargo of RNA. Notably, the transcriptome of platelets of HD_3 is located within the group of platelets derived from TPs. In addition, the transcriptomes of platelets obtained from TP_1 and TP_3 separate from the samples of the other TPs and orient towards the tumor tissue samples.

    Display full size

    Two-dimensional visualization of the PCA. For each sample, the RNA expression profile, which was detected by RNA-seq, is displayed as one point in a two-dimensional space. The first two principal components capture 92% of the variability. Samples with similar expression profiles are grouped together. The three test groups (TP-tissue; TP-platelets; HD-platelets) separate from each other, suggesting differences in the expression profiles within the analyzed groups. PC: principal component

    In order to figure out the cause of the separation of the 2 groups of platelet samples, RNA-seq data were analyzed further. Comparing the RNA expression profiles of platelets derived from HDs and TPs resulted in 426 significantly differentially existing RNAs (FDR < 0.05; Table S1). Among these 426 mRNA, 406 were abundantly detected in the group of platelets derived from TPs, and interestingly, 20 mRNA were present in lower amounts in platelets derived from TPs in comparison to that of HDs.

    By the use of the “Gene Ontology” platform, an enrichment analysis was carried out in order to get insights into the biological functions of identified differentially existing RNAs, which were assigned to a list of predefined biological terms [45, 46]. Term analyses revealed an association of detected RNA with squamous epithelium and keratinization. Likewise, RNA was also associated with terms of cell adhesion, signal transduction, or tumor-specific properties, such as proliferation, angiogenesis, and migration. Based on the information extracted from “Uniprot” and “NCBI Gene”, the 426 highly abundantly detected RNAs were further manually categorized, which confirmed this analysis: 49 RNAs encode genes that are characteristically expressed in epithelial cells and could have been transmitted by epithelial tumor cells (Table 5).

    Overview of genes with epithelial association

    TermGenes
    Epithelial differentiationZNF750 [49], SFN [50], FAT2 [51, 52], SBSN [53], GRHL1 [54], GRHL2 [55], TGM3 [56]
    Epithelial cell adhesionEPCAM [57], DSP [58], CLDN4 [59], DSG3 [60], PKP1 [61], DSC3 [60], DSG2 [60], CLDN1 [59], CDH1 [62], LAD1 [63], CXADR [64], PKP3 [61], F11R [65], PERP [66]
    CytokeratinKRT5 [67], KRT6a [67], KRT6b [67], KRT14 [67], KRT17 [67], KRT15 [67], KRT75 [67], KRT80 [67]
    KeratinizationPPL [68], BNC1 [69], CNFN [70], KAZN [71], SPRR1A [72], SPRR2A [72], SPRR2B [72], SPRR2G [72], SPRR2F [72], SPRR2E [72], SERPINB13 [73], KRTDAP [74], DMKN [75], EVPL [76], IRF6 [77], FERMT1 [78]
    Other epithelial characteristicsESRP1 [62], ESRP2 [62], LCN2 [79], CerS3 [80]
    Display full size

    Among the 426 differentially existing RNAs, 49 RNAs were found encoding genes characteristically expressed by epithelial cells. SFN: stratifin; SBSN: suprabasin; GRHL1: grainyhead like transcription factor 1; TGM3: transglutaminase 3; EPCAM: epithelial cell adhesion molecule; DSP: desmoplakin; CLDN4: claudin 4; PKP1: plakophilin 1; DSC3: desmocollin 3; LAD1: ladinin 1; CXADR: CXADR Ig-like cell adhesion molecule; F11R: F11 receptor; PERP: p53 apoptosis effector related to PMP22; PPL: periplakin; BNC1: basonuclin 1; CNFN: cornifelin; KAZN: kazrin; SPRR1A: small proline rich protein 1A; SERPINB13: serpin family B member 13; KRTDAP: keratinocyte differentiation associated protein; DMKN: dermokine; EVPL: envoplakin; IRF6: interferon regulatory factor 6; FERMT1: fermitin family member 1; ESRP1: epithelial splicing regulatory protein 1; LCN2: lipocalin 2; CerS3: ceramide synthase 3

    In addition to epithelial markers, a great part of differentially existing RNAs code for genes associated with the hallmarks of cancer (Table 6) [81].

    Selection of significantly differentially expressed genes that correspond to the characteristics of tumors

    GenesFunction in tumor cellsLog2 fold changeP adjustedReference
    WNT5aOncogen and tumor suppressor, associated with invasion, proliferation, and EMT5.7110.001[82]
    MMP3Associated with metastasis5.3140.002[83]
    SOX2Associated with tumor growth and metastasis5.2320.005[84]
    TMPRSS4Associated with tumor proliferation and metastasis5.1530.002[85]
    FOXE1Associated with tumor growth4.9490.020[86]
    FOXQ1Associated with metastasis4.8820.005[87]
    TNCAssociated with metastasis and EMT4.6530.002[88]
    TBX2Associated with tumor proliferation4.5080.007[89]
    CRYABAntiapoptotic factor4.4740.007[90]
    MMP2Associated with metastasis4.2490.002[83]
    PTHLHAssociated with tumor growth4.2440.041[91]
    EMP2Associated with angiogenesis4.2110.001[92]
    CCN1Associated with tumor growth4.2030.006[93]
    SPP1Associated with tumor growth and metastasis4.1960.018[94]
    GPNMBAssociated with metastasis4.0640.003[95]
    ANO1Associated with metastasis4.0280.014[96]
    MMP12Associated with metastasis3.9520.083[83]
    YAP1Associated with tumor growth3.8700.001[97]
    MMP13Associated with metastasis3.8320.077[83]
    SNAI2Associated with EMT and metastasis3.8130.012[98]
    Display full size

    EMT: epithelial mesenchymal transition; MMP3: matrix metalloproteinase 3; SOX2: SRY-box transcription factor 2; TMPRSS4: transmembrane serine protease 4; FOXE1: forkhead box E1; TNC: tenascin C; TBX2: T-box 2; CRYAB: crystallin alpha B; PTHLH: parathyroid hormone like hormone; EMP2: epithelial membrane protein 2; CCN1: cellular communication network factor 1; SPP1: secreted phosphoprotein 1; GPNMB: glycoprotein nmb; ANO1: anoctamine 1; YAP1: Yes associated protein 1; SNAI2: snail family transcriptional repressor 2

    Among these genes, MMPs, transcription factors, anti-apoptotic factors, and signaling proteins were identified. Some of these genes are even associated with a characteristic role in the tumorigenesis of head and neck carcinoma (like SOX2 [84], CRYAB [90], and PTHLH [91]). Interestingly, 41 non-coding RNAs were significantly enriched. Two-thirds of these encode long non-coding RNAs.

    A comparison of the level of mRNA existing in platelets from an HD, in TPs’ platelets and tumor tissue (Figure 2) revealed in some cases that the increase in platelet mRNA coding for certain genes is unlikely caused by the transfer of tumor-derived mRNA into platelets since the level of tumor RNA expression is low. This concerns for example the mRNA KCNQ1 opposite strand/antisense transcript 1 (KCNQ1OT1), MT-CO1 pseudogene 12 (MTCO1P12), MTCO1P35, and MTCO1P40. Therefore, genes abundantly existing in platelets from TPs compared to healthy platelets and present in tumor tissue were chosen for our further analyses. In doing so, RNA coding for 254 genes were identified, which will be further called “tumor-related mRNA” or “tumor-related genes” (Table S1–3).

    Display full size

    Expression profile of the top 20 differentially existing mRNA comparing platelets of HD and platelets of TPs in all test groups (TP-tumor tissue; TP-platelets; HD-platelets). Presented are box-plots. Each dot represents the value of 1 sample. The minimum, the lower quartile, the upper quartile and the maximum are shown. Depicted lines within the boxes represent the median. VST: variance stabilizing transformation; ARK1B10: aldo-keto reductase family 1 member B10; ARHGAP29: Rho GTPase activating protein 29; LAMB3: laminin subunit beta 3; LINCO2224: long intergenic non-protein coding RNA 2224; MTRN2L8: MT-RNR2 like 8 (pseudogene); NAMPTP1: nicotinamide phosphoribosyltransferase pseudogene 1; PITX1: paired like homeodomain 1; SNORA77: small nucleolar RNA, H/ACA box 77

    The volcano (Figure 3) plot illustrates all significantly differentially existing RNAs including “tumor-related mRNA” detected in platelets of TPs vs. platelets of HDs. Differentially existing and “tumor-related mRNA” are marked in red. Interestingly, the mRNA abundantly detected in tumor platelets lies in the region with a high log2-fold change of expression.

    Display full size

    Volcano plot visualizing all significantly differentially existing RNAs detected in platelets of TPs vs. platelets of HDs. The x-axis represents the log2 fold change. The y-axis represents the –log10 of the adjusted P values (FDR). The horizontal line indicates the FDR threshold of 0.05. RNA shown in blue (20) appeared reduced in platelets from TPs compared to healthy. All other genes are upregulated in tumor platelets. Genes shown in red (254) are considered to be tumor-related (abundantly found in both platelets and tumor tissue from TPs). Genes shown in yellow (152) are abundantly found in platelets from TPs but not highly expressed in tumor tissue. NS: not significant

    In order to analyze the expression profile of “tumor-related mRNA” even more, a heatmap using hierarchical clustering has been computed, showing the expression of all 254 tumor-related genes (Figure 4). A dendrogram shows a cluster built by the gene profiles of HD_1, HD_2, HD_4, and HD_5. The samples HD_1 and HD_5 bear the greatest similarities concerning their expression profiles. The profiles of the three TPs TP_4, TP_5, and TP_2 are also of high similarity. Notably, as already revealed by PCA (Figure 1), the profile of the HD_3 is close to that of the TPs. Apart from that patient, a separation of TPs and HDs is visible. In the samples, TP_3 and TP_1, almost every existing RNA, are higher expressed than in all the other TP samples. Also, a clear gradient increase in expression can be observed across the 3 groups.

    Display full size

    Heatmap of RNA abundance showing the 254 tumor-related and platelet-derived, significantly enriched RNA between the three groups (TP-platelets, HD-platelets, and TP-tumor tissue). HD_‘number’_p: HD_‘number’_platelets; TP_‘number’_p: tumor patient_‘number’_platelets; TP_‘number’_t: tumor patient_‘number’_tumor tissue

    In search of platelet-associated predictive markers allowing the diagnosis of HNSCC

    In order to identify a biomarker set as a potential diagnostic tool, on the top 30 “tumor-related” differential existing mRNA in platelets of TPs vs. healthy controls was focused (Figure 5). In order to select a robust set of mRNAs allowing the prediction of a tumor in the head and neck region with a very low error rate, those mRNA with a very low abundance in platelets of HDs were selected (n = 14; indicated with red arrows in Figure 5). Additionally, two additional mRNA based on high log2 fold change and their specific function in cancer of epithelial origin, KRT5 [67, 99] and ZNF750 [100–102] (indicated with blue arrows in Figure 5) were selected.

    Display full size

    Boxplots revealing the expression profile of the top 30 differentially existing “tumor-related mRNA” comparing platelets of HDs and platelets of TPs in all test groups (TP-tumor tissue; TP-platelets; HD-platelets). The genes are sorted by absolute log2 fold change in descending order. Presented are box-plots. Each dot represents the value of 1 sample. The minimum, the lower quartile, the upper quartile and the maximum are shown. Depicted lines within the boxes represent the median. Selection of mRNA with a very low abundance in platelets of HDs (n = 14; indicated with red arrows) and of two additional mRNA based on high log2 fold change and their specific function in cancer (n = 2; indicated with blue arrows) for identification of platelet-based HNSCC-specific biomarkers. TPT1P8: TPT1 pseudogene 8; IGSF3: immunoglobulin superfamily member 3; ENSG000002289: RNA gene which is affiliated with the long non-coding RNA (lncRNA) class; SPARCL1: SPARC like 1; PLAT: plasminogen activator, tissue type; S100A2: S100 calcium binding protein A2; PTPRZ1: protein tyrosine phosphatase receptor type Z1; PLAZG4E: phospholipase A2 group IVE; DCN: decorin

    The abundance of these mRNA was determined in a second cohort of HDs and TPs. The characteristics of the analyzed cohort here are depicted in Tables 4–6. In addition, the level of existing mRNA in a set of housekeeping genes in order to exclude those housekeepers showing variation in the platelet-derived mRNA abundance between HDs and TPs were checked. Thus, for further analyses and data normalization, GAPDH, UBC, HPRT, and PGK1 have been chosen (Figure 6).

    Display full size

    Boxplots reveal the expression profile of eight selected housekeeping genes. Their abundance was estimated by comparing all test groups (TP-tumor tissue; TP-platelets; HD-platelets). Presented are box-plots. Each dot represents the value of 1 sample. The minimum, the lower quartile, the upper quartile and the maximum are shown. Depicted lines within the boxes represent the median. ACTB: actin beta; RPL13: ribosomal protein L13; TBP: TATA binding protein

    First, primer pairs, which were designed using NCBI primer blast (Table 3), were tested using human Universal Reference RNA (Thermo Fisher) as well as the HNSCC cell lines UDSCC1 or UDSCC5, for which in a previous RNA-seq project the expression of genes of interest was detected. For 12 out of 16 genes, primer pairs useful for quantitative reverse transcription (qRT)-PCR could be established. For PLA2G4E, TMPRSS4, PITX1, and TPT1P8, two to four primer pairs were tested. Unfortunately, none of these primer pairs amplified a fragment of the expected size. For 5 out of 12 primer pairs, amplifying EGFR, PDPN, CDH11, CA12, and RBP1, no amplicon could be detected analyzing platelet-derived mRNA. However, using the left 7 primer pairs amplifying COL1A1, ZNF750, MAGEA6, Ly6D, WNT5a, FAT1, and KRT5 in all cases, except FAT1, a gain of predicted mRNA in platelets of TPs vs. platelets of HDs was observed (Figure 7).

    Display full size

    The relative abundance of target mRNA coding for COL1A1, ZNF50, MAGEA6, Ly6D, WNT5a, FAT1, and KRT5 estimated relative to the expression of housekeeping genes GAPDH, UBC, HPRT1, and GGK1 in platelets of HDs vs. TPs. Depicted lines within the boxes represent the median, and “+” represents the arithmetical mean. When significant differences were observed, the P value is highlighted in bold. The P value was calculated according to Mann-Whitney-U tets or Welch’s t-test depending on the data’s normal distribution

    Thus, the RNA-seq-data could be verified. Despite the fact, that in the group of HDs, one sample (HD_12; Table 2) behaved like a platelet mRNA derived from a TP, significant differences in COL1A1 and ZNF750 mRNA content in platelets of TPs vs. healthy controls were revealed (Figure 7). Static tests were applied again excluding the data derived from person HD_12. In this case, in 4 out of 7 analyzed genes, significant differences between mRNA abundance in platelets of HDs vs. TPs were revealed using GAPDH for normalization (Figure 8). In some cases, the significance level could not reach caused by a low sample number rather than a high sample variation. Unfortunately, due to the limited amount of mRNA extracted from platelets, not all TPs, as well as healthy individuals could be included in all analyses. For the same reason, not all predefined housekeeping mRNA was used for normalization. Of note, there are high variations in the abundance of analyzed mRNA in the TP’s group, independent of the tumor stage (data not shown).

    Display full size

    Presented are the relative abundance of target mRNA coding for COL1A1, ZNF50, MAGEA6, Ly6D, WNT5a, FAT1, and KRT5 estimated relative to the expression of housekeeping genes GAPDH, UBC, HPRT1, and GGK1 in platelets of HDs vs. TPs by boxplots. Here, the young HD (HD_12) with chronic tonsillitis and a smoking history was excluded from the HD group. Depicted lines within the boxes represent the median, and “+” represents the arithmetical mean. When significant differences were observed, the P value is indicated in bold. The P value was calculated according to Mann-Whitney-U or Welch’s t-test depending on the data’s normal distribution

    Taken together, in this pilot study differentially existing mRNA profiles between platelets of HDs vs. HNSCC patients were identified. Moreover, the basis for a biomarker-panel development useful for diagnosis and possibly progress or therapy control of head and neck tumors using liquid biopsy-obtained platelets was created.

    Discussion

    The current knowledge of thrombocyte biology, especially in the context of cancer, still leaves many unanswered questions. In addition, the mechanism responsible for changes in the RNA profile in platelets of TPs is not completely understood. However, the quite big array of publications over recent years describing differences in mRNA expression between HDs and patients of different tumor entities is not to be neglected. Also, the comparative analysis of the RNA profiles of platelets of HDs vs. HNSCC TPs revealed 426 differentially existing RNA. Among them, 254 “tumor-related RNA”, coding for epithelial markers, and, importantly, for tumor-associated genes were identified. Particularly, cytokeratin 5 and 14 are characteristic elements of the squamous epithelium [67], the origin of HNSCC. Additionally, the presence of RNA, which encodes oncogenes or tumor suppressor genes that are characteristically expressed in HNSCC, like SOX2 [84], CRYAB [90], and PTHLH [91], further suggests transmission of tumor RNA into platelets [25, 81]. Therefore, the existence of epithelial and tumor marker RNA in TPs’ platelets may cause by the previous observations that tumor cells transfer their RNA into platelets [25].

    Furthermore, 41 non-coding RNA forms were significantly differentially enriched in platelets derived from TPs vs. HD platelets, which could be explained by the fact that some of these forms are also poly-adenylated and therefore addressed by this approach [103]. Also, in the study of Best et al. [26], 20 differentially expressed non-coding RNAs were detected. The significantly altered existence of long non-coding RNA in platelets of TPs was also confirmed by Luo et al. [104]. In the present work, as well as in Best et al. [26], among others, the long intergenic non-protein coding RNA 152 (LINC00152) was abundantly found in TPs’ platelets, a long non-coding RNA expressed in many types of cancer, including oral squamous cell carcinoma [105, 106], acting as a potent oncogene involved in cancer cell proliferation, invasiveness, and motility.

    Together, the presence of characteristics for cancer cells non-coding RNA as well as of RNA coding for epithelial and cancer-specific genes could indicate a transfer of RNA from nuclear tumor cells. However, these data should be interpreted with caution as the existence of platelet-coated CTC and an association of tumor-derived extracellular vesicles with platelets have been described [107, 108]. Therefore, contamination of platelet RNA with RNA from other tumor-derived sources cannot be excluded.

    Despite the small cohort size, qRT-PCR results could validate significantly differentially existing mRNAs for 2 out of the 7 examined genes, namely ZNF750 and COL1A1. These are very promising findings: whereas ZNF750 is a lineage-specific tumor suppressor in squamous cell carcinoma [49, 100–102], COL1A1 is associated with many tumor entities including oral and pancreatic cancer with prognostic relevance [109–111]. Additionally, in gastric cancer, COL1A1 is highly expressed and associated with cancer cell invasion, metastasis, and poor prognosis [112, 113]. COL1A is also a driver of oncogenesis in hepatocellular carcinoma [114]. Importantly, COL1A1 was recently identified as a candidate biomarker for HNSCC [99] and associated with cancer stem cell signature in HNSCC [115]. Thus, COL1A1 may serve as a possible diagnostic as well as a prognostic biomarker not only for HNSCC but also for other solid tumor entities. It, therefore, seems reasonable to implement the molecular COL1A1 characterization as a biomarker in the routine clinical/pathological molecular HNSCC diagnostic.

    Although non-significant, many genes demonstrate a clear trend, where platelets of TPs show elevated mRNA levels compared to those of healthy individuals, this could be seen in almost all cases in varying degrees, except for FAT1. Reasons for non-significant results could be variations in individual patients or the small sample sizes, which is a limitation of the present study. Although sample size calculations according to Lin et al. [37] were performed before starting the here described pilot study, not all hypotheses could be checked with the estimated sample size. However, it is remarkable that, despite this limitation, the abundance of tumor-related RNA could be verified in a second, equally small cohort.

    After taking a closer look at the individual results, two HDs appeared to behave like a TP (cohort 1: HD_3; cohort 2: HD_12), whose platelets were showing substantial “tumor-related mRNA” abundance. There was an attempt to acquire a more recent blood sample in order to investigate whether any malignancies are present now, but unfortunately, the donors have not agreed. While no further details on the health status of HD_3 is available, the HD HD_12 was a chronic tonsillitis patient, 22 years old of age at the time of the blood draw (5 years ago), with a tobacco and alcohol consumption history. Since the results obtained from platelet RNA of HD_12 were suspicious, it was hypothesized that this individual had a pathological condition influencing platelet physiology. This is likely since changes in gene expression profile, as well as platelet count in patients of chronic tonsillitis and other diseases, were observed [116–118].

    However, based on an HNSCC-causing history in the case of HD_12, the possibility cannot be excluded, that with the method applied here, a pre-cancerogenic processes in that person classified as a HD was identified. To prove this assumption, prospective studies should be conducted analyzing platelet mRNA content routinely and the patients should be observed for years.

    Additionally, one TP (TP_8) could be considered a HD. One could assume that with early tumor stage less tumor RNA can be taken up by platelets. However, this association turned out to be inconsistent in the here performed analyses. But certain drugs inhibit the function of platelets which could affect the interaction of tumor cells and platelets [119], and consequently the amount of transferred mRNA.

    Statistical procedures exist to detect outliers [120]. However, the data presented here stem from a small sample size and most do not follow a normal distribution, so it was not possible to reliably apply these methods. However, in order to determine sample sizes for this pilot study we followed the suggestions of Lin et al. [37].

    The first and fundamental question this work attempted to answer was whether a transfer of tumor RNA to platelets actually takes place since platelets from head and neck cancer patients have not been analyzed so far. And this study was able to answer this question positively using this small cohort. Additionally, applying this study design, a set of RNA that was specified as “tumor-related” in a second, again small cohort, was identified as at least enriched in patient platelets. Precisely because the data were collected using small cohorts, these findings show the great possible diagnostic potential of platelet analysis in head and neck cancer patients.

    A strength of the present study is the inclusion of tumor RNA into the sequencing analyses and not just comparing platelets of TPs with that of HDs. Doing so, it was revealed that differentially existing RNA in platelets is not exclusively caused by the transfer from tumor to platelets. This observation has to consider for further studies on TEP. There are many factors influencing the platelet transcriptome even under non-pathological conditions. According to Simon et al. [121], the mRNA profiles of platelets even differ depending on the age of the patients, their ethnicity, or their gender. As already mentioned, different diseases can cause transcriptomic changes in platelets [116–118]. Moreover, there is a certain platelet transcriptome variance observed even in HDs, probably depending on the platelet subpopulation [122, 123]. All these factors could contribute to individual variations and to the initial observation of the present study that not all differentially existing RNA observed in platelets of TPs are characterized by a gain of RNA. Thus, there is an urgent need for robust biomarkers allowing the identification of defined diseases, including HNSCC.

    Despite the limitations of this work, the presented results are very promising, because so far, there is a lack of effective screening options for the diagnostics and therapy monitoring of HNSCCs [124] and liquid biopsy poses a great potential for the future of cancer diagnostics. The results presented here provide a basis for a biomarker panel that shows a potential ability for platelets to distinguish between healthy and TPs in HNSCC tumors with a simple blood draw.

    With the help of future research, the analysis of platelet RNA in head and neck TPs could generate a powerful method for monitoring tumorigenesis at any time and in real-time. As was suggested also for other cancer entities, tumor-specific therapies could be applied in a targeted manner and modified with respect to response rate or resistance development [6, 34, 125].

    Abbreviations

    B2M:

    beta 2-microglobulin
    CDH11 :

    cadherin 11
    CLDN4 :

    claudin 4
    COL1A1:

    collagen type I alpha 1 chain
    CRYAB :

    crystallin alpha B
    CTCs:

    circulating tumor cells
    DSP :

    desmoplakin
    EGFR:

    epidermal growth factor receptor
    EMT:

    epithelial mesenchymal transition
    FAT1 :

    FAT atypical cadherin 1
    FDR:

    false discovery rate
    GAPDH :

    glyceraldehyde-3-phosphate dehydrogenase
    HDs:

    healthy donors
    HNSCC:

    head and neck squamous cell carcinoma
    HPRT1 :

    hypoxanthine phosphoribosyltransferase 1
    HPV:

    human papilloma virus
    KRT5 :

    keratin 5
    Ly6D :

    lymphocyte antigen 6 family member D
    MAGEA6 :

    MAGE family member A6
    MMP3 :

    matrix metalloproteinase 3
    mRNA:

    messenger RNA
    MTCO1P12 :

    MT-CO1 pseudogene 12
    NCBI:

    national center for biotechnology information
    PCA:

    principal component analysis
    PCR:

    polymerase chain reaction
    PDPN :

    podoplanin
    PGK1 :

    phosphoglycerate kinase 1
    PITX1 :

    paired like homeodomain 1
    PLA2G4E :

    phospholipase A2 group IVE
    PTHLH :

    parathyroid hormone like hormone
    py:

    pack years
    RBP1 :

    retinol binding protein 1
    RNA-seq:

    RNA-sequencing
    SOX2 :

    SRY-box transcription factor 2
    SPRR1A :

    small proline rich protein 1A
    TEP:

    tumor-educated platelets
    TMPRSS4 :

    transmembrane serine protease 4
    TPs:

    tumor patients
    UBC :

    ubiquitin C
    WNT5a :

    Wnt family member 5a
    ZNF750 :

    zinc finger protein 750

    Supplementary materials

    The supplementary material for this article is available at: https://www.explorationpub.com/uploads/Article/file/1002143_sup_1.xlsx.

    Declarations

    Acknowledgments

    Thanks are due to Ivonne Görlich, Katja Hasch, Gabriela Cudek and Monika Jerg for excellent technical assistance.

    Author contributions

    HAK and CB: Conceptualization, Methodology, Data Curation, Supervision, Project administration. HAK, MG, JMK and JE: Software. HAK, JMK, MG, LTH, AW, JE and CB: Validation, Investigation. JMK, MG, LTH, JE and AW: Formal analysis. HAK, TKH and CB: Resources. LTH and CB: Writing—original draft. HAK, JMK, MG, AW, JE, BW, SL and TKH: Writing—review and editing. JMK, JE and LTH: Visualization. All authors read and approved the submitted version.

    Conflicts of interest

    The authors declare that they have no conflicts of interest.

    Ethical approval

    The study was approved by the local ethics committee of the Ulm University [90/15] and complies with the Declaration of Helsinki.

    Consent to participate

    Informed consent to participate in the study was obtained from all participants.

    Consent to publication

    Not applicable.

    Availability of data and materials

    RNA-seq data are available from Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE156554) under Record ID GSE156554.

    Funding

    Not applicable.

    Copyright

    © The Author(s) 2023.

    References

    Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71:209–49. [DOI] [PubMed]
    Johnson DE, Burtness B, Leemans CR, Lui VWY, Bauman JE, Grandis JR. Head and neck squamous cell carcinoma. Nat Rev Dis Primers. 2020;6:92. [DOI] [PubMed] [PMC]
    Bhatia A, Burtness B. Treating head and neck cancer in the age of immunotherapy: a 2023 update. Drugs. 2023;83:217–48. [DOI] [PubMed]
    Ranganath K, Feng AL, Franco RA, Varvares MA, Faquin WC, Naunheim MR, et al. Molecular biomarkers of malignant transformation in head and neck dysplasia. Cancers (Basel). 2022;14:5581. [DOI] [PubMed] [PMC]
    Corcoran RB. Liquid biopsy versus tumor biopsy for clinical-trial recruitment. Nat Med. 2020;26:1815–6. [DOI] [PubMed]
    Li G, Sun Y. Liquid biopsy: advances, limitations and clinical applications. JSM Biotechnol Bioeng. 2017;4:1078.
    Di Meo A, Bartlett J, Cheng Y, Pasic MD, Yousef GM. Liquid biopsy: a step forward towards precision medicine in urologic malignancies. Mol Cancer. 2017;16:80. [DOI] [PubMed] [PMC]
    Stroun M, Anker P, Maurice P, Lyautey J, Lederrey C, Beljanski M. Neoplastic characteristics of the DNA found in the plasma of cancer patients. Oncology. 1989;46:318–22. [DOI] [PubMed]
    Skog J, Würdinger T, van Rijn S, Meijer DH, Gainche L, Curry WT Jr, et al. Glioblastoma microvesicles transport RNA and proteins that promote tumour growth and provide diagnostic biomarkers. Nat Cell Biol. 2008;10:1470–6. [DOI] [PubMed] [PMC]
    Racila E, Euhus D, Weiss AJ, Rao C, McConnell J, Terstappen LWMM, et al. Detection and characterization of carcinoma cells in the blood. Proc Natl Acad Sci U S A. 1998;95:4589–94. [DOI] [PubMed] [PMC]
    Cohen JD, Li L, Wang Y, Thoburn C, Danilova L, Douville C, et al. Detection and localization of surgically resectable cancers with a multi-analyte blood test. Science. 2018;359:926–30. [DOI] [PubMed] [PMC]
    Allard WJ, Matera J, Miller MC, Repollet M, Connelly MC, Rao C, et al. Tumor cells circulate in the peripheral blood of all major carcinomas but not in healthy subjects or patients with nonmalignant diseases. Clin Cancer Res. 2004;10:6897–904. [DOI] [PubMed]
    Hanjani NA, Esmaelizad N, Zanganeh S, Gharavi AT, Heidarizadeh P, Radfar M, et al. Emerging role of exosomes as biomarkers in cancer treatment and diagnosis. Crit Rev Oncol Hematol. 2022;169:103565. [DOI] [PubMed]
    Yang WY, Feng LF, Meng X, Chen R, Xu WH, Hou J, et al. Liquid biopsy in head and neck squamous cell carcinoma: circulating tumor cells, circulating tumor DNA, and exosomes. Expert Rev Mol Diagn. 2020;20:1213–27. [DOI] [PubMed]
    Payne K, Spruce R, Beggs A, Sharma N, Kong A, Martin T, et al. Circulating tumor DNA as a biomarker and liquid biopsy in head and neck squamous cell carcinoma. Head Neck. 2018;40:1598–604. [DOI] [PubMed]
    Liebs S, Eder T, Klauschen F, Schütte M, Yaspo ML, Keilholz U, et al. Applicability of liquid biopsies to represent the mutational profile of tumor tissue from different cancer entities. Oncogene. 2021;40:5204–12. [DOI] [PubMed] [PMC]
    Lousada-Fernandez F, Rapado-Gonzalez O, Lopez-Cedrun JL, Lopez-Lopez R, Muinelo-Romay L, Suarez-Cunqueiro MM. Liquid biopsy in oral cancer. Int J Mol Sci. 2018;19:1704. [DOI] [PubMed] [PMC]
    Han J, Lin Q, Dong C. Plasma cell-free circRNAs panel act as fingerprint predicts the occurrence of laryngeal squamous cell carcinoma. Aging (Albany NY). 2021;13:17328–36. [DOI] [PubMed] [PMC]
    Gauthier A, Philouze P, Lauret A, Alphonse G, Malesys C, Ardail D, et al. Circulating tumor cell detection during neoadjuvant chemotherapy to predict early response in locally advanced oropharyngeal cancers: a prospective pilot study. J Pers Med. 2022;12:445. [DOI] [PubMed] [PMC]
    Flach S, Howarth K, Hackinger S, Pipinikas C, Ellis P, McLay K, et al. Liquid BIOpsy for MiNimal RESidual DiSease Detection in Head and Neck Squamous Cell Carcinoma (LIONESS)—a personalised circulating tumour DNA analysis in head and neck squamous cell carcinoma. Br J Cancer. 2022;126:1186–95. [DOI] [PubMed] [PMC]
    Hofmann L, Ludwig S, Vahl JM, Brunner C, Hoffmann TK, Theodoraki MN. The emerging role of exosomes in diagnosis, prognosis, and therapy in head and neck cancer. Int J Mol Sci. 2020;21:4072. [DOI] [PubMed] [PMC]
    Teng Y, Gao L, Loveless R, Rodrigo JP, Strojan P, Willems SM, et al. The hidden link of exosomes to head and neck cancer. Cancers (Basel). 2021;13:5802. [DOI] [PubMed] [PMC]
    Battinelli EM, Markens BA, Italiano JE Jr. Release of angiogenesis regulatory proteins from platelet alpha granules: modulation of physiologic and pathologic angiogenesis. Blood. 2011;118:1359–69. [DOI] [PubMed] [PMC]
    Calverley DC, Phang TL, Choudhury QG, Gao B, Oton AB, Weyant MJ, et al. Significant downregulation of platelet gene expression in metastatic lung cancer. Clin Transl Sci. 2010;3:227–32. [DOI] [PubMed] [PMC]
    Nilsson RJA, Balaj L, Hulleman E, van Rijn S, Pegtel DM, Walraven M, et al. Blood platelets contain tumor-derived RNA biomarkers. Blood. 2011;118:3680–3. [DOI] [PubMed] [PMC]
    Best MG, Sol N, Kooi I, Tannous J, Westerman BA, Rustenburg F, et al. RNA-seq of tumor-educated platelets enables blood-based pan-cancer, multiclass, and molecular pathway cancer diagnostics. Cancer Cell. 2015;28:666–76. [DOI] [PubMed] [PMC]
    Best MG, Sol N, In ‘t Veld SGJG, Vancura A, Muller M, Niemeijer ALN, et al. Swarm intelligence-enhanced detection of non-small-cell lung cancer using tumor-educated platelets. Cancer Cell. 2017;32:238–52. [DOI] [PubMed] [PMC]
    Sheng M, Dong Z, Xie Y. Identification of tumor-educated platelet biomarkers of non-small-cell lung cancer. Onco Targets Ther. 2018;11:8143–51. [DOI] [PubMed] [PMC]
    Xing S, Zeng T, Xue N, He Y, Lai Yz, Li HL, et al. Development and validation of tumor-educated blood platelets integrin alpha 2b (ITGA2B) RNA for diagnosis and prognosis of non-small-cell lung cancer through RNA-seq. Int J Biol Sci. 2019;15:1977–92. [DOI] [PubMed] [PMC]
    Sol N, in ‘t Veld SGJG, Vancura A, Tjerkstra M, Leurs C, Rustenburg F, et al. Tumor-educated platelet RNA for the detection and (pseudo)progression monitoring of glioblastoma. Cell Rep Med. 2020;1:100101. [DOI] [PubMed] [PMC]
    Heinhuis KM, In ’t Veld SGJG, Dwarshuis G, van den Broek D, Sol N, Best MG, et al. RNA-sequencing of tumor-educated platelets, a novel biomarker for blood based sarcoma diagnostics. Cancer (Basel). 2020;12:1372. [DOI] [PubMed] [PMC]
    Yang L, Jiang Q, Li DZ, Zhou X, Yu DS, Zhong J. TIMP1 mRNA in tumor-educated platelets is diagnostic biomarker for colorectal cancer. Aging (Albany NY). 2019;11:8998–9012. [DOI] [PubMed] [PMC]
    Zhu B, Gu S, Wu X, He W, Zhou H. Bioinformatics analysis of tumor-educated platelet microRNAs in patients with hepatocellular carcinoma. Biosci Rep. 2021;41:BSR20211420. [DOI] [PubMed] [PMC]
    Waqar W, Asghar S, Manzoor S. Platelets’ RNA as biomarker trove for differentiation of early-stage hepatocellular carcinoma from underlying cirrhotic nodules. PLoS One. 2021;16:e0256739. [DOI] [PubMed] [PMC]
    Best MG, Wurdinger T. Tumor-educated platelets for the earlier detection of hepatocellular carcinoma. Clin Res Hepatol Gastroenterol. 2020;44:794–5. [DOI] [PubMed]
    Asghar S, Waqar W, Umar M, Manzoor S. Tumor educated platelets, a promising source for early detection of hepatocellular carcinoma: liquid biopsy an alternative approach to tissue biopsy. Clin Res Hepatol Gastroenterol. 2020;44:836–44. [DOI] [PubMed]
    Lin CW, Liao SG, Liu P, Lee MLT, Park YS, Tseng GC. RNASeqDesign: a framework for ribonucleic acid sequencing genomewide power calculation and study design issues. J R Stat Soc Ser C Appl Stat. 2019;68:683–704. [DOI] [PubMed] [PMC]
    Rosner B. Fundamentals of biostatistics. 8th ed. Boston: Cengage learning; 2015.
    Bentley DR, Balasubramanian S, Swerdlow HP, Smith GP, Milton J, Brown CG, et al. Accurate whole human genome sequencing using reversible terminator chemistry. Nature. 2008;456:53–9. [DOI] [PubMed] [PMC]
    TruSeq® stranded mRNA sample preparation guide [Internet]. San Diego: illunina®; c2013 [cited 2022 Nov 15]. Available from: https://ostr.ccr.cancer.gov/wp-content/uploads/2020/05/truseq-stranded-mrna-sample-prep-guide-15031047-e.pdf
    Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9. [DOI] [PubMed] [PMC]
    Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9. [DOI] [PubMed] [PMC]
    Zhu A, Ibrahim JG, Love MI. Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics. 2019;35:2084–92. [DOI] [PubMed] [PMC]
    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. [DOI] [PubMed] [PMC]
    Schmid F, Schmid M, Müssel C, Sträng JE, Buske C, Bullinger L, et al. GiANT: gene set uncertainty in enrichment analysis. Bioinformatics. 2016;32:1891–4. [DOI] [PubMed]
    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9. [DOI] [PubMed] [PMC]
    Lever J, Krzywinski M, Altman N. Points of significance: principal component analysis. Nat Methods. 2017;14:641–2. [DOI]
    Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45. [DOI] [PubMed] [PMC]
    Sen GL, Boxer LD, Webster DE, Bussat RT, Qu K, Zarnegar BJ, et al. ZNF750 is a p63 target gene that induces KLF4 to drive terminal epidermal differentiation. Dev Cell. 2012;22:669–77. [DOI] [PubMed] [PMC]
    Kim S, Wong P, Coulombe PA. A keratin cytoskeletal protein regulates protein synthesis and epithelial cell growth. Nature. 2006;441:362–5. [DOI] [PubMed]
    Matsui S, Utani A, Takahashi K, Mukoyama Y, Miyachi Y, Matsuyoshi N. Human Fat2 is localized at immature adherens junctions in epidermal keratinocytes. J Dermatol Sci. 2007;48:233–6. [DOI] [PubMed]
    Matsui S, Utani A, Takahashi K, Mukoyama Y, Miyachi Y, Matsuyoshi N. Knockdown of Fat2 by siRNA inhibits the migration of human squamous carcinoma cells. J Dermatol Sci. 2008;51:207–10. [DOI] [PubMed]
    Park GT, Lim SE, Jang SI, Morasso MI. Suprabasin, a novel epidermal differentiation marker and potential cornified envelope precursor. J Biol Chem. 2002;277:45195–202. [DOI] [PubMed] [PMC]
    Frisch SM, Farris JC, Pifer PM. Roles of Grainyhead-like transcription factors in cancer. Oncogene. 2017;36:6067–73. [DOI] [PubMed]
    Chen W, Liu ZX, Oh JE, Shin KH, Kim RH, Jiang M, et al. Grainyhead-like 2 (GRHL2) inhibits keratinocyte differentiation through epigenetic mechanism. Cell Death Dis. 2012;3:e450. [DOI] [PubMed] [PMC]
    Eckert RL, Sturniolo MT, Broome AM, Ruse M, Rorke EA. Transglutaminase function in epidermis. J Invest Dermatol. 2005;124:481–92. [DOI] [PubMed]
    Laimer K, Fong D, Gastl G, Obrist P, Kloss F, Tuli T, et al. EpCAM expression in squamous cell carcinoma of the oral cavity: frequency and relationship to clinicopathologic features. Oral Oncol. 2008;44:72–7. [DOI] [PubMed]
    Wanuske MT, Brantschen D, Schinner C, Stüdle C, Walter E, Hiermaier M, et al. Clustering of desmosomal cadherins by desmoplakin is essential for cell-cell adhesion. Acta Physiol (Oxf). 2022;231:e13609. [DOI] [PubMed]
    De Vicente JC, Fernández-Valle Á, Vivanco-Allende B, Rodríguez Santamarta T, Lequerica-Fernández P, Hernández-Vallejo G, et al. The prognostic role of claudins -1 and -4 in oral squamous cell carcinoma. Anticancer Res. 2015;35:2949–59. [PubMed]
    Harrison OJ, Brasch J, Lasso G, Katsamba PS, Ahlsen G, Honig B, et al. Structural basis of adhesive binding by desmocollins and desmogleins. Proc Natl Acad Sci U S A. 2016;113:7160–5. [DOI] [PubMed] [PMC]
    Hatzfeld M, Wolf A, Keil R. Plakophilins in desmosomal adhesion and signaling. Cell Commun Adhes. 2014;21:25–42. [DOI] [PubMed]
    Barriere G, Fici P, Gallerani G, Fabbri F, Zoli W, Rigaud M. Circulating tumor cells and epithelial, mesenchymal and stemness markers: characterization of cell subpopulations. Ann Transl Med. 2014;2:109. [DOI] [PubMed] [PMC]
    Marinkovich MP, Taylor TB, Keene DR, Burgeson RE, Zone JJ. LAD-1, the linear IgA bullous dermatosis autoantigen, is a novel 120-kDa anchoring filament protein synthesized by epidermal cells. J Invest Dermatol. 1996;106:734–8. [DOI] [PubMed]
    Cohen CJ, Shieh JTC, Pickles RJ, Okegawa T, Hsieh JT, Bergelson JM. The coxsackievirus and adenovirus receptor is a transmembrane component of the tight junction. Proc Natl Acad Sci U S A. 2001;98:15191–6. [DOI] [PubMed] [PMC]
    Itoh M, Sasaki H, Furuse M, Ozaki H, Kita T, Tsukita S. Junctional adhesion molecule (JAM) binds to PAR-3: a possible mechanism for the recruitment of PAR-3 to tight junctions. J Cell Biol. 2001;154:491–8. [DOI] [PubMed] [PMC]
    Ihrie RA, Marques MR, Nguyen BT, Horner JS, Papazoglu C, Bronson RT. Perp is a p63-regulated gene essential for epithelial integrity. Cell. 2005;120:843–56. [DOI] [PubMed]
    Moll R, Divo M, Langbein L. The human keratins: biology and pathology. Histochem Cell Biol. 2008;129:705–33. [DOI] [PubMed] [PMC]
    Ruhrberg C, Hajibagheri MAN, Parry DAD, Watt FM. Periplakin, a novel component of cornified envelopes and desmosomes that belongs to the plakin family and forms complexes with envoplakin. J Cell Biol. 1997;139:1835–49. [DOI] [PubMed] [PMC]
    Tseng H, Green H. Association of basonuclin with ability of keratinocytes to multiply and with absence of terminal differentiation. J Cell Biol. 1994;126:495–506. [DOI] [PubMed] [PMC]
    Michibata H, Chiba H, Wakimoto K, Seishima M, Kawasaki S, Okubo K, et al. Identification and characterization of a novel component of the cornified envelope, cornifelin. Biochem Biophys Res Commun. 2004;318:803–13. [DOI] [PubMed]
    Groot KR, Sevilla LM, Nishi K, DiColandrea T, Watt FM. Kazrin, a novel periplakin-interacting protein associated with desmosomes and the keratinocyte plasma membrane. J Cell Biol. 2004;166:653–9. [DOI] [PubMed] [PMC]
    Steinert PM, Marekov LN. The proteins elafin, filaggrin, keratin intermediate filaments, loricrin, and small proline-rich proteins 1 and 2 are isodipeptide cross-linked components of the human epidermal cornified cell envelope. J Biol Chem. 1995;270:17702–11. [DOI] [PubMed]
    Welss T, Sun J, Irving JA, Blum R, Smith AI, Whisstock JC, et al. Hurpin is a selective inhibitor of lysosomal cathepsin L and protects keratinocytes from ultraviolet-induced apoptosis. Biochemistry. 2003;42:7381–9. [DOI] [PubMed]
    Edqvist PHD, Fagerberg L, Hallström BM, Danielsson A, Edlund K, Uhlén M, et al. Expression of human skin-specific genes defined by transcriptomics and antibody-based profiling. J Histochem Cytochem. 2015;63:129–41. [DOI] [PubMed] [PMC]
    Utsunomiya A, Chino T, Utsunomiya N, Luong VH, Tokuriki A, Naganuma T, et al. Homeostatic function of dermokine in the skin barrier and inflammation. J Invest Dermatol. 2020;140:838–49.E9. [DOI] [PubMed]
    Ruhrberg C, Hajibagheri MA, Simon M, Dooley TP, Watt FM. Envoplakin, a novel precursor of the cornified envelope that has homology to desmoplakin. J Cell Biol. 1996;134:715–29. [DOI] [PubMed] [PMC]
    Richardson RJ, Dixon J, Malhotra S, Hardman MJ, Knowles L, Boot-Handford RP, et al. Irf6 is a key determinant of the keratinocyte proliferation-differentiation switch. Nat Genet. 2006;38:1329–34. [DOI] [PubMed]
    Herz C, Aumailley M, Schulte C, Schlötzer-Schrehardt U, Bruckner-Tuderman L, Has C. Kindlin-1 is a phosphoprotein involved in regulation of polarity, proliferation, and motility of epidermal keratinocytes. J Biol Chem. 2006;281:36082–90. [DOI] [PubMed]
    Mongre RK, Sodhi SS, Sharma N, Ghosh M, Kim JH, Kim N, et al. Epigenetic induction of epithelial to mesenchymal transition by LCN2 mediates metastasis and tumorigenesis, which is abrogated by NF-κB inhibitor BRM270 in a xenograft model of lung adenocarcinoma. Int J Oncol. 2016;48:84–98. [DOI] [PubMed] [PMC]
    Jennemann R, Rabionet M, Gorgas K, Epstein S, Dalpke A, Rothermel U, et al. Loss of ceramide synthase 3 causes lethal skin barrier disruption. Hum Mol Genet. 2012;21:586–608. [DOI] [PubMed]
    Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144:646–74. [DOI] [PubMed]
    Asem MS, Buechler S, Wates RB, Miller DL, Stack MS. Wnt5a signaling in cancer. Cancers (Basel). 2016;8:79. [DOI] [PubMed] [PMC]
    Stamenkovic I. Matrix metalloproteinases in tumor invasion and metastasis. Cancer Biol. 2000;10:415–33. [DOI] [PubMed]
    Tang X, Shen X, Li L, Zhang Y, Chen G. SOX2 overexpression correlates with poor prognosis in laryngeal squamous cell carcinoma. Auris Nasus Larynx. 2013;40:481–6. [DOI] [PubMed]
    Jung H, Lee KP, Park SJ, Park JH, Jang Ys, Choi SY, et al. TMPRSS4 promotes invasion, migration and metastasis of human tumor cells by facilitating an epithelial–mesenchymal transition. Oncogene. 2008;27:2635–47. [DOI] [PubMed]
    Fan Y, Ding Z, Yang Z, Deng X, Kang J, Wu B, et al. Expression and clinical significance of FOXE1 in papillary thyroid carcinoma. Mol Med Rep. 2013;8:123–7. [DOI] [PubMed]
    Liu JY, Wu XY, Wu GN, Liu FK, Yao XQ. FOXQ1 promotes cancer metastasis by PI3K/AKT signaling regulation in colorectal carcinoma. Am J Transl Res. 2017;9:2207–18. [PubMed] [PMC]
    Nagaharu K, Zhang X, Yoshida T, Katoh D, Hanamura N, Kozuka Y, et al. Tenascin C induces epithelial-mesenchymal transition-like change accompanied by SRC activation and focal adhesion kinase phosphorylation in human breast cancer cells. Am J Pathol. 2011;178:754–63. [DOI] [PubMed] [PMC]
    Abrahams A, Parker MI, Prince S. The T-box transcription factor Tbx2: its role in development and possible implication in cancer. IUBMB Life. 2010;62:92–102. [DOI] [PubMed]
    Chin D, Boyle GM, Williams RM, Ferguson K, Pandeya N, Pedley J, et al. Alpha B-crystallin, a new independent marker for poor prognosis in head and neck cancer. Laryngoscope. 2005;115:1239–42. [DOI] [PubMed]
    Chang WM, Lin YF, Su CY, Peng HY, Chang YC, Hsiao JR, et al. Parathyroid hormone-like hormone is a poor prognosis marker of head and neck cancer and promotes cell growth via RUNX2 regulation. Sci Rep. 2017;7:41131. [DOI] [PubMed] [PMC]
    Qin Y, Takahashi M, Sheets K, Soto H, Tsui J, Pelargos P, et al. Epithelial membrane protein-2 (EMP2) promotes angiogenesis in glioblastoma multiforme. J Neurooncol. 2017;134:29–40. [DOI] [PubMed] [PMC]
    Kok SH, Chang HH, Tsai JY, Hung HC, Lin CY, Chiang CP, et al. Expression of Cyr61 (CCN1) in human oral squamous cell carcinoma: an independent marker for poor prognosis. Head Neck. 2010;32:1665–73. [DOI] [PubMed]
    Shevde LA, Samant RS. Role of osteopontin in the pathophysiology of cancer. Matrix Biol. 2014;37:131–41. [DOI] [PubMed] [PMC]
    Arosarena OA, dela Cadena RA, Denny MF, Bryant E, Barr EW, Thorpe R, et al. Osteoactivin promotes migration of oral squamous cell carcinomas. J Cell Physiol. 2016;231:1761–70. [DOI] [PubMed] [PMC]
    Ruiz C, Martins JR, Rudin F, Schneider S, Dietsche T, Fischer CA, et al. Enhanced expression of ANO1 in head and neck squamous cell carcinoma causes cell migration and correlates with poor prognosis. PLoS One. 2012;7:e43265. [DOI] [PubMed] [PMC]
    Lee TF, Tseng YC, Chang WC, Chen YC, Kao YR, Chou TY, et al. YAP1 is essential for tumor growth and is a potential therapeutic target for EGFR-dependent lung adenocarcinomas. Oncotarget. 2017;8:89539–51. [DOI] [PubMed] [PMC]
    Peinado H, Marin F, Cubillo E, Stark HJ, Fusenig N, Nieto MA, et al. Snail and E47 repressors of E-cadherin induce distinct invasive and angiogenic properties in vivo. J Cell Sci. 2004;117:2827–39. [DOI] [PubMed]
    Kengkarn S, Petmitr S, Boonyuen U, Reamtong O, Poomsawat S, Sanguansin S. Identification of novel candidate biomarkers for oral squamous cell carcinoma based on whole gene expression profiling. Pathol Oncol Res. 2020;26:2315–25. [DOI] [PubMed]
    Hazawa M, Lin DC, Handral H, Xu L, Chen Y, Jiang YY, et al. ZNF750 is a lineage-specific tumour suppressor in squamous cell carcinoma. Oncogene. 2017;36:2243–54. [DOI] [PubMed] [PMC]
    Nambara S, Masuda T, Tobo T, Kidogami S, Komatsu H, Sugimachi K, et al. Clinical significance of ZNF750 gene expression, a novel tumor suppressor gene, in esophageal squamous cell carcinoma. Oncol Lett. 2017;14:1795–801. [DOI] [PubMed] [PMC]
    Otsuka R, Akutsu Y, Sakata H, Hanari N, Murakami K, Kano M, et al. ZNF750 expression as a novel candidate biomarker of chemoradiosensitivity in esophageal squamous cell carcinoma. Oncol. 2017;93:197–203. [DOI] [PubMed]
    Guttman M, Amit I, Garber M, French C, Lin MF, Feldser D, et al. Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature. 2009;458:223–7. [DOI] [PubMed] [PMC]
    Luo CL, Xu ZG, Chen H, Ji J, Wang YH, Hu W, et al. LncRNAs and EGFRvIII sequestered in TEPs enable blood-based NSCLC diagnosis. Cancer Manag Res. 2018;10:1449–59. [DOI] [PubMed] [PMC]
    Seo D, Kim D, Kim W. Long non-coding RNA linc00152 acting as a promising oncogene in cancer progression. Genomics Inform. 2019;17:e36. [DOI] [PubMed] [PMC]
    Zou X, Hu X, He F, Zhang M, Kong X, Rui S, et al. LncRNA LINC00152 promotes oral squamous cell carcinoma growth via enhancing upstream transcription factor 1 mediated mitochondrial ribosomal protein L52 transcription. J Oral Path Med. 2021;51:454–63. [DOI] [PubMed]
    Rikkert LG, van der Pol E, van Leeuwen TG, Nieuwland R, Coumans FAW. Centrifugation affects the purity of liquid biopsy-based tumor biomarkers. Cytometry A. 2018;93:1207–12. [DOI] [PubMed] [PMC]
    Placke T, Örgel M, Schaller M, Jung G, Rammensee HG, Kopp HG, et al. Platelet-derived MHC class I confers a pseudonormal phenotype to cancer cells that subverts the antitumor reactivity of natural killer immune cells. Cancer Res. 2012;72:440–8. [DOI] [PubMed]
    He B, Lin X, Tian F, Yu W, Qiao B. MiR-133a-3p inhibits oral squamous cell carcinoma (OSCC) proliferation and invasion by suppressing COL1A1. J Cell Biochem. 2018;119:338–46. [DOI] [PubMed]
    Wiechec E, Magan M, Matic N, Ansell-Schultz A, Kankainen M, Monni O, et al. Cancer-associated fibroblasts modulate transcriptional signatures involved in proliferation, differentiation and metastasis in head and neck squamous cell carcinoma. Cancers (Basel). 2021;13:3361. [DOI] [PubMed] [PMC]
    Li J, Ding Y, Li A. Identification of COL1A1 and COL1A2 as candidate prognostic factors in gastric cancer. World J Surg Oncol. 2016;14:297. [DOI] [PubMed] [PMC]
    Sun H. Identification of key genes associated with gastric cancer based on DNA microarray data. Oncol Lett. 2016;11:525–30. [DOI] [PubMed] [PMC]
    Li X, Liu L, Song X, Wang K, Niu L, Xie L, et al. TEP linc-GTF2H2-1, RP3-466P17.2, and lnc-ST8SIA4-12 as novel biomarkers for lung cancer diagnosis and progression prediction. J Cancer Res Clin Oncol. 2021;147:1609–22. [DOI] [PubMed]
    Ma HP, Chang HL, Bamodu OA, Yadav VK, Huang TY, Wu ATH, et al. Collagen 1A1 (COL1A1) is a reliable biomarker and putative therapeutic target for hepatocellular carcinogenesis and metastasis. Cancers (Basel). 2019;11:786. [DOI] [PubMed] [PMC]
    Wu ZH, Li C, Zhang YJ, Zhou W. Identification of a cancer stem cells signature of head and neck squamous cell carcinoma. Front Genet. 2022;13:81477. [DOI] [PubMed] [PMC]
    Lood C, Amisten S, Gullstrand B, Jönsen A, Allhorn M, Truedsson L, et al. Platelet transcriptional profile and protein expression in patients with systemic lupus erythematosus: up-regulation of the type I interferon system is strongly associated with vascular disease. Blood. 2010;116:1951–7. [DOI] [PubMed]
    Cengiz C, Erhan Y, Murat T, Ercan A, Ibrahim S, Ihsan G, et al. Values of mean platelet volume in patients with chronic tonsillitis and adenoid hypertrophy. Pakistan J Med Sci. 2013;29:569–72. [DOI] [PubMed] [PMC]
    Thompson RW, Gungor A. Immune thrombocytopenia of childhood responsive to tonsillectomy in the setting of chronic tonsillitis: a case report and literature review. Am J Otolaryngol. 2017;38:639–41. [DOI] [PubMed]
    Yan M, Lesyk G, Radziwon-Balicka A, Jurasz P. Pharmacological regulation of platelet factors that influence tumor angiogenesis. Semin Oncol. 2014;41:370–7. [DOI] [PubMed]
    Zimek A, Filzmoser P. There and back again: outlier detection between statistical reasoning and data mining algorithms. Wiley Interdiscip Rev Data Min Knowl Discov. 2018;8:e1280. [DOI]
    Simon LM, Edelstein LC, Nagalla S, Woodley AB, Chen ES, Kong X, et al. Human platelet microRNA-mRNA networks associated with age and gender revealed by integrated plateletomics. Blood. 2014;123:e37–45. [DOI] [PubMed] [PMC]
    Clancy L, Beaulieu LM, Tanriverdi K, Freedman JE. The role of RNA uptake in platelet heterogeneity. Thromb Haemost. 2017;117:948–61. [DOI] [PubMed]
    Supernat A, Popęda M, Pastuszak K, Best MG, Grešner P, Veld SI’, et al. Transcriptomic landscape of blood platelets in healthy donors. Sci Rep. 2021;11:15679. [DOI] [PubMed] [PMC]
    Kujan O, Glenny A, Duxbury A, Thakker N, Sloan P, Duxbury J. Screening programmes for the early detection and prevention of oral cancer. Cochrane Database Syst Rev. 2003;2003:CD004150. [DOI] [PubMed]
    Nilsson RJA, Karachaliou N, Berenguer J, Gimenez-Capitan A, Schellen P, Teixido C, et al. Rearranged EML4-ALK fusion transcripts sequester in circulating blood platelets and enable blood-based crizotinib response monitoring in non-small-cell lung cancer. Oncotarget. 2016;7:1066–75. [DOI] [PubMed] [PMC]
Download PDF
Share
Subscribe