Background Immunotherapy has become an important treatment option for patients with advanced non-small cell lung cancer (NSCLC). At present, none of these existing biomarkers can effectively stratify true responders and there is an urgent need for identifying novel biomarkers. Exosomes derived from the serum of patients with cancer have been proven to be reliable markers for cancer diagnosis. Here, we explored the possibility of using plasma-derived exosomal microRNAs as potential biomarkers for optimal selection of patients with advanced EGFR/ALK negative NSCLC to immunotherapy.
Methods From June 2017 to February 2019, 30 patients with advanced EGFR/ALK wild-type (WT) NSCLC who received PD-1/PD-L1 inhibitors were enrolled. The efficacy evaluation was conducted after every three cycles of treatment according to RECIST 1.1. Plasma samples of these patients were collected before the administration of PD-1/PD-L1 inhibitors as baseline, and after every three cycles if the patients achieved partial response (PR) or complete response. Plasma from seven healthy individuals was also collected as normal control. Exosomes were prepared by ultracentrifugation followed by total RNA extraction, and exosome-derived miRNAs were profiled using small RNA next-generation sequencing followed by differential expression analysis.
Results In order to identify biomarker for better response, all five patients who achieved PR and four patients with progressive disease (PD) at efficacy evaluation were included for differential expression analysis. Based on unsupervised hierarchical clustering, exosomal miRNA expression profile was significantly altered in patients with NSCLC compared with normal controls with a total of 155 differentially expressed exosomal miRNAs. Interestingly, hsa-miR-320d, hsa-miR-320c, and hsa-miR-320b were identified significantly upregulated in the PD groups compared with the PR group at baseline before the treatment. In addition, we identified that hsa-miR-125b-5p, a T-cell suppressor, showed a trend of increased expression in the PD group at baseline and was significantly downregulated in the post-treatment plasma exosomes compared with pre-treatment samples of the PR patients.
Conclusion Patients with NSCLC represent unique plasma exosomal miRNA profiles. Hsa-miR-320d, hsa-miR-320c, and hsa-miR-320b were identified as potential biomarkers for predicting the efficacy of immunotherapy in advanced NSCLCs. When T-cell suppressor hsa-miR-125b-5p was downregulated during the treatment, the patients may obtain increased T-cell function and respond well to immunotherapy.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
With the advent of immune checkpoint blockade therapy, a new therapeutic era has opened for patients with cancer. Immunotherapy has shown promising results in various types of cancer, and immune checkpoint inhibitor (ICI) drugs targeting the PD-1/PD-L1 (programmed cell death protein 1, CD279; programmed death-ligand 1, CD274)1 2 and CTLA-4 (cytotoxic T-lymphocyte-associated protein 4, CD152)3 have been approved for treating patients with advanced melanoma,4 5 metastatic renal cell carcinoma,6 and advanced non-small cell lung cancer (NSCLC).7 ICIs showed a survival advantage in comparison with chemotherapy in first-line and second-line setting of patients with NSCLC; however, overall response rate is around 20% with a large proportion of patients undergoing disease progression within the first week of treatment.8 9 To spare patients from ineffective therapy, the search for predictive biomarkers to select patients who would preferentially benefit from immunotherapy becomes important.
The expression of PD-L1 on tumor cell membranes has been most widely studied as the predictive biomarker for anti-PD-1/PD-L1 therapy.10 11 However, challenges remain mainly for several reasons. First, the use of different PD-L1 antibodies and corresponding cut-offs of its expression level makes it hard to compare and standardize.12 Second, the expression of PD-L1 can be dynamic and inducible during the recognition of T-cell antigen which may interfere with treatment response.13 Third, the presence of high heterogeneity of PD-L1 expression between primary tumor and metastatic tumor in the same patient limits the robustness of the testing.14 Finally, certain patients do not have tumor tissues available for PD-L1 test. High tumor mutational burden,15 16 high microsatellite instability, and mismatch repair deficiency17 18 are other emerging biomarkers that are significantly associated with response to PD-1/PD-L1 blockade immunotherapy although viable predictive values were observed within different cancer types and clinical trials. High degree of tumor-infiltrating lymphocytes has also been proven to be associated with improved survival in patients with NSCLC and is considered as a potential biomarker for checkpoint inhibitor immunotherapy.19 20 However, the establishment of a clinically feasible cut-off is very difficult due to the overlap of baseline T-cell density between patients with or without response. So far, the use of a single biomarker alone has not been sufficient for selection of patients with anti-PD-1/PD-L1 therapy. Further investigations for more robust biomarkers to better stratify the patients are needed.
Exosomes are small membrane vesicles with diameters of 30–100 nm that can be released by most cell types including cancer cells, which present in various body fluids including plasma, saliva, urine, and malignant effusions.21 Exosomes can transfer functional protein, mRNA, and microRNA (miRNA) to neighboring cells, fulfilling its role as the mediator of intercellular communication.22 miRNAs are enriched in exosomes and can regulate targeted mRNA translation or stability, and are implicated in the regulation of maturation, proliferation, differentiation, and activation of immune cells.23 24 The structure of exosomes protects miRNAs from degradation that is suitable for clinical detection.
In patients with cancer, exosomes derived from the serum of patients with cancer have been proven to be reliable markers for cancer diagnosis.25 26 Cancer-secreted exosomes can function as cell-to-cell mediators of metastasis, which prime pre-metastatic niche formation in distant organs27 28 and are responsible for engraftment and survival of the metastatic cells.29 In this study, we are aiming to explore potential plasma-derived exosomal miRNA signatures as potential biomarkers for the optimal selection of patients with advanced NSCLC to immunotherapy.
Patients and clinical sample collection
Thirty patients with EGFR/ALK negative lung cancer with available clinical information including age, gender, stage, treatment history, and baseline plasma samples were enrolled in this study in Guangdong Provincial People’s Hospital from June 2017 to February 2019. Peripheral blood was collected from each patient on a regular basis for routine clinical care. Plasma sample was prepared within 2 hours of blood drawn and then stored at −80°C. In this study, plasma samples of patients with advanced EGFR/ALK wild-type (WT) NSCLC were collected before the administration of PD-1/PD-L1 inhibitors as baseline. The efficacy evaluation was conducted after three cycles of treatment. For every three cycles, plasma samples were collected from patients until the disease progressed. Patients who achieved partial response (PR) or complete response (CR) were included in this study as responders, compared with patients with progressive disease (PD) on treatment as non-responders. Flow chart for patient selection and exclusion criteria is shown in online supplementary figure 1. No CR was observed in this patient cohort. In total, five patients who achieved PR and four patients with PD at efficacy evaluation were included in this study. Plasma samples from seven healthy individuals were also collected as normal controls.
Online supplementary table 1 lists the baseline clinicopathological characteristics of all the patients enrolled in this study. formalin-fixed paraffin-embedded(FFPE) slides of the tumor samples were prepared for PD-L1 expression evaluation using immunohistochemistry with the following PD-L1 antibody clones: SP142 or SP263 (Roche, USA), and 28-8 (Abcam, USA). For the nine patients (responders and non-responders) included in this study, seven patients were analyzed using SP263 antibody, while the other two were analyzed with SP142 or 28-8, respectively. Patients were treated with different immunotherapy drugs targeting PD-1/PD-L1. The median follow-up time was 8 months, ranging from 1 month to approximately 21 months. Except one patient with lymphoepithelioma-like carcinoma, five patients were diagnosed with adenocarcinoma, while the other three patients were diagnosed with squamous cell carcinoma. The median age at diagnosis is 55, ranging from 47 to 68.
Plasma exosomal RNA isolation and small RNA sequencing
One milliliter of plasma sample was centrifuged at 10,000×g for 30 min at 4°C to remove any cell debris. The collected supernatant was then subjected for ultra-high-speed centrifugation at 150,000×g for 70 min at 4°C. The pellet containing exosome was resuspended in 200 µL Phosphate buffered saline(PBS) for downstream applications. Total RNA including miRNA was extracted from plasma-derived exosome using miRNeasy Serum/Plasma Kit (QIAGEN) following the manufacturer’s instructions. The quantification and size distribution of the extraction were analyzed by Qubit V.4.0 and Agilent Bioanalyzer 2100 (Agilent), respectively. Quantified RNA was subjected for sequencing library preparation using NEBNext Small RNA Library Prep Set for Illumina (NEB Biolabs) following the manufacturer’s instructions. Briefly, isolated total RNA was subjected for 3′ and 5′ adaptor ligation, followed by 17 cycles of PCR amplification. PCR products from library preparation were subjected for gel electrophoresis on 6% Novex TBE PAGE gel (Thermo Fisher Scientific) and DNA fragments between 140 and 150 bp were recovered from the gel. Purified small RNA cDNA library was quantified by Qubit V.4.0 and the size distribution was analyzed on Agilent Bioanalyzer 2100, followed by sequencing on Illumina HiSeq4000 platform.
Characterization of purified exosomes
For verification of purified exosomes using electron microscopy, purified exosomes suspended in PBS were dropped on copper-coated grids. After staining with 2% uranyl acetate, grids were dried at room temperature and visualized using a Hitachi H-7650 transmission electron microscope. miRNA-seq data analysis miRNA identification and reads counting in each miRNA were performed using miRDeep2.30 After trimming the 3′ adaptor sequence, all sequences ranging in length from 18 to 26 nt were recorded in a non-redundant file along with reads count. To identify known miRNAs, the miRNA tags were aligned against miRNA precursor sequences reported in the miRNA database “miRBase” (release 21) using the “quantifier.pl” script within miRDeep2. Differential expression (DE) analysis of miRNA sequence data was performed with the Bioconductor package edgeR.31 miRNAs with read counts per million mapped reads (CPM) ≥5 in at least 20% of all samples were identified as expressed miRNAs. DE between different groups was evaluated by fitting a negative binomial generalized linear model and then adjusting the p value for multiple testing using the Benjamini-Hochberg correction with a false discovery rate (FDR) of 0.1 and a minimum log2(CPM) of 4.
Exosomal miRNA expression profiles
All five patients who achieved PR (responders) and four patients with PD (non-responders) were enrolled in the exosomal miRNA differential expression analysis. Baseline clinicopathological characteristics of patients enrolled are listed in online supplementary table 1 with detailed treatment history of the patients listed in online supplementary table 2. Nine baseline pre-treatment plasma samples and the post-treatment plasma samples from four responders as well as normal control samples from seven healthy individuals were collected for exosome purification and exosomal miRNA profiling. We purified exosomes from plasma by differential centrifugation and verified them as small double-leaflet membrane particles (30–200 nm) by transmission electron microscopy (TEM) (figure 1).
Raw reads of miRNA-seq from all plasma samples were normalized to CPM and 706 miRNAs were retained as expressed miRNAs in plasma exosomes (online supplementary table 3). Using hierarchical clustering, healthy individuals could be distinguished from patients with lung cancer based on expressed miRNAs profile in plasma-derived exosomes (figure 2). We further analyzed differentially expressed miRNAs (DE miRNAs) between healthy control samples and patient samples. As shown in online supplementary table 4, 155 DE miRNAs were identified and the top 10 upregulated and top 10 downregulated miRNAs are listed in online supplementary table 4a. Some of these top regulated miRNAs were highly related to immune pathways. Among them, hsa-miR-146a has been well studied in autoimmune disease pathogenesis and is related to NF-κB pathway involved in immune response.32 hsa-miR-221 and hsa-miR-224 have been reported to be related to anti-TNF treatment in inflammatory bowel disease.33 hsa-miR-224 has been proven to induce cell proliferation by target Smad4 through SMAD/TGF-β pathway.34 Meanwhile, three of the top regulated miRNAs which were downregulated in patients with lung cancer have been reported as tumor suppressors in different cancer types, including hsa-miR-144,35–37 hsa-miR-363-3p,38 and hsa-miR-150.39 40
Kyoto Encyclopedia of Genes and Genomes(KEGG) pathway analysis was performed with 155 differentially expressed exosomal miRNAs using DIANA tool (online supplementary table 4b). Cancer-related pathways were significantly enriched in lung cancer patient group including TGF-β, MAPK, Hippo, and Ras signaling pathways. Immune-related pathways including TNF signaling pathway, B-cell and T-cell receptor signaling pathways were also enriched in the patient group. In conclusion, plasma-derived exosomes’ miRNA expression profile was altered in patients with EGFR/ALK negative lung cancer and differently expressed miRNAs were enriched in cancer-related and immune-related pathways compared with healthy control.
Different exosomal miRNA expression profiles between PR and PD patients
To identify the potential predictors of anti-PD-1 therapy in patients with lung cancer, we divided patient samples into PR group and PD group. Then we compared the miRNA profile between PR group and PD group using exact test. With cut-off at logCPM>4, p<0.05, FDR≤0.1, hsa-miR-320d was considered as the most significant DE miRNA followed by hsa-miR-642a-3p, hsa-miR-320c, and hsa-miR-320b with FDR around 0.1 (table 1).
Next, we examined the miRNA expression level of hsa-miR-642a-3p, hsa-miR-320c, hsa-miR-320b, and hsa-miR-642a-3p in overall patient cohort and healthy population (H). We separated all the samples from patients into three groups: PR pre-treatment samples (PR-pre), PD pre-treatment samples (PD-pre), and PR post-treatment samples (PR-post). As shown in figure 2, hsa-miR-320d, hsa-miR-320c, and hsa-miR-320b belonged to hsa-miR-320 family and all showed significant difference in mRNA expression level between PD-pre group and PR-pre group. Moreover, compared with healthy control and PR group (PR-pre and PR-post), hsa-miR-320d, hsa-miR-320c, and hsa-miR-320b showed upregulation in PD-pre group (figure 3). hsa-miR-643a-3p was eliminated due to its low expression level and no significant difference in miRNA expression level among these three groups (figure 3D). In conclusion, high level of hsa-miR-320d, hsa-miR-320c, and hsa-miR-320b might correlate with unfavorable response to anti-PD-1 treatment and hsa-miR-320d, hsa-miR-320c, and hsa-miR-320b might be potential predictors of anti-PD-1 therapy.
Identification of plasma-derived exosomal miRNAs suppressed by anti-PD-1 therapy in patients with lung cancer
Paired pre-treatment and post-treatment plasma samples were collected from the four patients who achieved PR on anti-PD-1 therapy (online supplementary table 1). To identify potential exosomal miRNAs influenced by anti-PD-1 therapy during the treatment course, we compared the exosomal miRNA profiles before (PR-pre) and after (PR-post) anti-PD-1 treatment. With cut-off at logCPM>4, p<0.05, FDR≤0.1, only hsa-miR-125b-5p was differentially expressed between PR-pre and PR-post samples (table 2 and figure 4A), which was dramatically downregulated in the PR-post samples and suggested that the level of exosomal hsa-miR-125b-5p might be used as a potential indicator to monitor the efficacy of anti-PD-1 treatment. Interestingly, plasma samples from the PD group at baseline exhibited a trend of increased level of hsa-miR125b-5p compared with PR-pre group and healthy controls (figure 4A).
This study has demonstrated that hsa-miR-125b-5p could induce γδ T-cell apoptosis and downregulate γδ T-cell activation and cytotoxicity to tumor cells by decreasing degranulation and inhibiting the secretion of IFN-γ and TNF-ɑ41 (figure 4B). An increased exosomal hsa-miR-125b-5p level in PD group at baseline might indicate impaired γδ T-cell function and cytotoxicity to tumor cells that could confer resistance to immunotherapy. On the other hand, a significantly reduced exosomal hsa-miR-125b-5p level could indicate the activation of T cells and cytotoxicity to tumor cells in response to immunotherapy.
Here, we identified three miRNAs from hsa-miR-320 family as potential predictors and one potential target hsa-miR-125b-5p for anti-PD-1 treatment. The three miRNAs, hsa-miR-320d, hsa-miR-320c, and hsa-miR-320b, all exhibited upregulation in PD group compared with PR group and correlated with unfavorable response to anti-PD-1 treatment. hsa-miR-125b-5p was downregulated in patients with response to anti-PD-1 treatment which might be a target for anti-PD-1 therapy. Meanwhile, hsa-miR-125b-5p exhibited a trend of increase in PD group compared with PR group and healthy individual group. Based on the results, patients with low level of hsa-miR-320d, hsa-miR-320c, hsa-miR-320b, and hsa-miR-125b-5p might be better candidates for anti-PD-1 treatment. During the anti-PD-1 treatment, continuous decrease in the level of hsa-miR-125b-5p may be taken into consideration as a sign of better outcome and longer progression-free survival(PFS). Especially, in the patient with delayed clinical responses or pseudoprogression, the level of hsa-miR-125b-5p might also be used as a potential prognostic marker to indicate whether the patient would benefit from anti-PD-1 treatment.
The employment of plasma-derived exosomal miRNAs in anti-PD-1 therapy holds some advantages in clinical circumstances. First, the plasma-derived exosomes detected in patients are mainly released from tumor cells which can be a more accurate and dynamic reflection of the status and function of tumor cells, rather than tumor tissue. Second, miRNAs in exosomes are more stable compared with miRNAs in blood due to the double-membrane structure. Furthermore, the detection of exosomal miRNAs only requires patient blood sample which is able to overcome spatial heterogeneity of tissue samples and timely monitor tumor progression during the whole treatment course. In future application, we can go further to simplify this method by using qPCR to monitor miRNA biomarkers in plasma-derived exosomes.
The miR320 family which consists of five members (hsa-miR-320a, hsa-miR-320b, hsa-miR-320c, hsa-miR-320d, hsa-miR-320e) is conserved but only exists in vertebrates from Xenopus to humans. hsa-miR-320d and hsa-miR-320e only present in primates and humans.42 High expression of hsa-miR-320 has been reported to be associated with negative prognosis and high risk of metastasis in ovarian cancer.43 Overexpression of hsa-miR-320 promotes B-cell proliferation through suppressing Phosphatase and tensin homolog(PTEN) expression and promotes cyclin D1 expression.44 Our data also confirmed that high expression of hsa-miR-320d, hsa-miR-320c, and hsa-miR-320b has a negative prognosis and associated with bad response to anti-PD-1 treatment in patients with NSCLC. The observation of downregulated hsa-miR-125b-5p in PR-pre group also supported the role of hsa-miR-125b-5p as a T-cell suppressor in another study. We could further validate these observations by detecting the level of hsa-miR-320d, hsa-miR-320c, hsa-miR-320b, and hsa-miR-125b-5p in the PD group after anti-PD-1 treatment. However, plasma samples from PD group after anti-PD-1 treatment were not available for analysis in this study because the evaluation was stopped after the patient showed PD. Future studies need to be carried out for further investigation and validation.
Unfortunately, we did not compare these potential biomarkers with other markers to explore their prediction accuracy. The other limitation of our study is that the cohort we used is small. There are other clinicopathological factors such as smoking status45 and treatment history that may also contribute to the upregulation or downregulation of exosomal miRNAs. Therefore, these four miRNAs need to be tested in other cohorts to validate their roles as potential immunotherapy biomarkers.
In conclusion, we identified four plasma-derived exosomal miRNAs as potential biomarkers of anti-PD-1 therapy for patients with NSCLC which can be indicators for non-invasive longitudinal liquid biopsy to ruling in or ruling out patients undergo immunotherapy. Large cohort validation and functional study of these potential biomarkers are required before it can be applied to clinical use.
The exosomal miRNA profile of patients with NSCLC differs from healthy individuals. We identified three miRNAs from hsa-miR-320 family as potential predictors and one potential target hsa-miR-125b-5p for anti-PD-1 treatment in plasma exosomes of patients with NSCLC. Plasma-derived exosomal miRNAs may be potential biomarkers for anti-PD-1 therapy in EGFR/ALK wild-type advanced NSCLC.
X-XP and RY contributed equally.
Contributors X-XP collected and interpreted the patient data, and wrote the manuscript. R- YY analyzed data and wrote the manuscript. XW, S-YW, CP, Z-HC, X-CZ, YWS, LL, and Y-LW reviewed and edited the manuscript. C-YG was responsible for TEM data of the manuscript. QZ was responsible for conceptualization, funding, and methodology of the manuscript. All authors read and approved the final manuscript.
Funding The study was supported by National Key R&D Program of China (grant no. 2016YFC1303800 to QZ), National Natural Science Foundation of China (no. 81871891 to QZ), and High-level Hospital Construction Project (grant no. DFJH201810 to QZ).
Competing interests XW and R-YY are employees of Geneseeq Technology Inc. YWS, S-YW, and C-YG are shareholders or employees of Nanjing Geneseeq Technology Inc. Y-LW has received speaker fees from AstraZeneca, Roche, Eli Lilly, Pfizer, and Sanofi. QZ declares speaker fees from AstraZeneca and Roche.
Patient consent for publication Obtained.
Ethics approval The study was approved by the Ethics Committee at Guangdong Provincial People's Hospital.
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement Data are available on reasonable request.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.