Malignant pleural effusion (MPE) is common in advanced lung cancer and often represents the only accessible source of tumor material. Beyond confirming malignancy, molecular characterization is increasingly important for guiding targeted therapies. The comparative performance of pleural fluid (PF) and plasma cell-free DNA (cfDNA) for detecting clinically relevant mutations remains uncertain.
MethodsWe conducted a prospective study of 70 patients with non-small cell lung cancer (NSCLC) and pleural effusion who underwent diagnostic thoracentesis. Paired PF and plasma samples were analyzed using next-generation sequencing (NGS) panels targeting clinically relevant driver and resistance mutations. Mutation detection rates, variant allele frequencies (VAFs), and concordance between PF and plasma were assessed. Pleural effusions were classified as definite or probable malignant according to pathological or predefined clinical criteria.
ResultsAmong patients with pathologically confirmed MPE (n=49), PF cfDNA detected clinically relevant mutations more frequently than plasma cfDNA (57% vs 41%), yielding a higher number of mutations (40 vs 27) and higher VAFs. Concordance between PF and plasma was low, with only 26% of mutations detected in both specimens. Among patients with probable MPE (n=21), clinically relevant mutations were identified in 11 cases (52%) when PF and plasma cfDNA were considered together. Despite frequent detection of actionable or resistance-associated mutations, only 16% of patients received targeted therapy.
ConclusionsPleural fluid cfDNA analysis detects clinically relevant mutations more frequently than plasma and provides complementary molecular information in lung cancer-associated pleural effusions, including cases with negative cytology, potentially expanding access to targeted therapies.
Malignant pleural effusion (MPE) is a frequent complication of advanced NSCLC and often provides the first or only accessible specimen for diagnosis. Although demonstration of malignant cells remains essential, the advent of targeted therapies and immunotherapy has made molecular profiling equally critical for treatment decision-making. Identification of oncogenic driver alterations, as well as assessment of programmed death ligand-1 (PD-L1) expression, may directly influence therapeutic selection [1].
Pleural fluid is particularly suitable for this purpose. Cell blocks permit immunohistochemical assessment, whereas analysis of cfDNA from the supernatant enables NGS to detect both therapeutically relevant driver mutations and resistance-associated mutations [2]. Molecular profiling of PF has shown high concordance with tumor tissue and may be performed when tissue sampling is limited, technically challenging, or contraindicated [3,4]. Moreover, emerging evidence suggests that cfDNA analysis may identify clinically relevant mutations even in cytology-negative effusions in selected cases, although current data remain limited and heterogeneous [4,5]. This raises the possibility of expanding molecular characterization without requiring additional invasive procedures; however, the diagnostic yield and clinical implications of this approach remain incompletely defined.
Plasma cfDNA is widely used for molecular profiling in advanced NSCLC. However, its performance relative to PF, particularly regarding mutation detection rates, VAF, and concordance, remains incompletely established.
This study aimed to compare PF and plasma cfDNA for detection of actionable mutations in lung cancer-associated MPE using a comprehensive NGS panel and to determine whether one specimen type is superior or whether both sources provide complementary molecular information.
MethodsStudy populationThis prospective study included consecutive patients with NSCLC and pleural effusion who underwent diagnostic thoracentesis at Hospital Universitari Arnau de Vilanova (Lleida, Spain) between 2018 and 2022. All patients had advanced disease and were classified as stage IV when MPE was pathologically confirmed. Patients with probable MPE were considered clinically consistent with advanced-stage disease. The study protocol was approved by the local ethics committee (CEIC No. 1947), and all patients provided written informed consent.
Diagnostic criteriaMalignant pleural effusion (MPE) was diagnosed by identification of malignant cells in PF smear or cell block specimens or by positive pleural biopsy when PF cytology was negative. Pleural effusion was classified as probably malignant when both of the following criteria were met: (1) histologically confirmed lung cancer obtained from lung tissue or an extrapleural metastatic site (e.g., liver, lymph nodes, or soft tissue); and (2) a lymphocytic exudate with negative PF cytology in which alternative causes of pleural effusion had been reasonably excluded through comprehensive PF analysis and imaging studies.
Data collectionClinical data collected included age, sex, radiographic size and laterality of pleural effusion, cytopathologic findings from PF and pleural biopsy specimens (when performed), histologic subtype of the primary lung tumor, PD-L1 status, and oncologic treatments administered. Information regarding prior systemic therapies and previous molecular testing was collected when available. Pleural effusion size was categorized as occupying<50% or ≥50% of the hemithorax on chest radiograph. NGS results were not systematically available to treating physicians at the time of clinical decision-making and were not used to guide treatment selection within the study protocol.
Pleural fluid supernatant and plasma cfDNA analysisPleural fluid and peripheral blood were collected at the time of diagnostic thoracentesis and processed within 2h. Pleural fluid and plasma samples were centrifuged using standard double-spin protocols, aliquoted, and stored at −80°C until analysis. Cell-free DNA was extracted from 5mL of PF supernatant and plasma using a commercially available kit (QIAamp Circulating Nucleic Acid Kit; Qiagen) according to the manufacturer's instructions. Cell-free DNA concentration was assessed by fluorometry, and fragment-size distribution and cfDNA quality were evaluated by capillary electrophoresis, including assessment of mono-nucleosomal fragment proportions. All samples met quality criteria for subsequent NGS analysis at the central sequencing facility.
Next-generation sequencing was performed at a centralized facility (Qiagen Genomic Services, Germany). Libraries were prepared using the QIAseq Targeted DNA Panel Library Prep Kit and sequenced on an Illumina NovaSeq 6000 platform using the Human Lung Cancer Panel (72 genes) (Table 1). Sequencing data were processed using standard bioinformatic pipelines for quality control, alignment, and somatic variant calling, with reads mapped to the human reference genome (hg19). Variant allele frequency was defined as the proportion of sequencing reads harboring a given variant relative to the total number of reads covering the genomic position of interest. Full technical details are provided in the Supplementary Material.
Genes included in the sequencing panel.
| AKT1 | CDKN2A | FBXO7 | KDR | MGA | NTRK2 | PTPRD | SETD2 |
| ALK | CDKN2B | FBXW7 | KEAP1 | MLH1 | NTRK3 | RARB | SMAD4 |
| AMER1 | CREBBP | FGFR1 | KIT | MUC16 | PDGFRA | RASSF1 | SMARCA4 |
| APC | CTNNB1 | FGFR2 | KMT2D | MYC | PIK3CA | RB1 | SOX2 |
| ARID1A | DDR2 | FGFR3 | KRAS | NF1 | PIK3CG | RBM10 | STK11 |
| ATM | EGFR | FHIT | LRP1B | NFE2L2 | PIK3R1 | RET | TNFAIP3 |
| BAI3 | EPHA5 | GRM8 | MAP2K1 | NOTCH1 | PIK3R2 | RIT1 | TP53 |
| BAP1 | ERBB2 | HRAS | MDM2 | NRAS | PKHD1 | ROS1 | TSC1 |
| BRAF | ERBB4 | JAK2 | MET | NTRK1 | PTEN | RUNX1T1 | U2AF1 |
Analyses focused on clinically relevant mutations with established therapeutic or resistance implications in NSCLC according to contemporary National Comprehensive Cancer Network (NCCN) guidelines [1]. These included actionable driver alterations (e.g., EGFR, ALK, ROS1, BRAF, KRAS, MET, RET, ERBB2, and NTRK fusions), as well as mutations associated with resistance to tyrosine kinase inhibitors (e.g., EGFR T790M and C797S, MET amplification, kinase domain mutations in ALK and ROS1, and KRAS variants). Alterations not covered by the sequencing panel (e.g., NRG1 gene fusion) were not evaluated.
PD-L1 assessmentProgrammed death ligand-1 expression was assessed by immunohistochemistry in available tumor tissue (primary or metastatic) and reported as tumor proportion score (TPS). In PF supernatants and plasma, PD-L1 levels were measured using a commercially available enzyme-linked immunosorbent assay (Human PD-L1 DuoSet ELISA, R&D Systems, Minneapolis, Minnesota, USA) and expressed as pg/mL.
Statistical analysisPaired comparisons between PF and plasma were performed using the Wilcoxon signed-rank test. Correlations were evaluated using Pearson or Spearman correlation coefficients, as appropriate. Agreement between sample types was assessed using contingency tables and Cohen κ coefficient. An oncoprint-style visualization was generated to descriptively illustrate patient-level mutation profiles across sample types. Statistical analyses were performed using SPSS version 29.0 and R software version 4.3.2.
ResultsBaseline characteristicsSeventy patients with stage IV NSCLC and pleural effusion were included. Adenocarcinoma was the predominant histologic subtype. The median age of the cohort was 70.5 years, and most patients were male. Malignant pleural effusion was pathologically confirmed in 49 patients (70%): 47 by positive PF cytology and 2 by pleural biopsy after negative cytology. In the remaining 21 patients (30%), pleural effusion was classified as probably malignant according to predefined clinical criteria.
During the study period, targeted therapy was administered to 11 patients (16%). Nine patients were treatment-naïve at the time of PF sampling and received EGFR-targeted therapy based on molecular findings from PF or tissue specimens. Two patients had received prior EGFR-targeted therapy before NGS analysis. In 1 previously treated patient, cfDNA analysis identified resistance-associated EGFR mutations (T790M and C797S). Median follow-up was 5.5 months, during which most patients died. Baseline characteristics are summarized in Table 2.
Baseline characteristics of the study population.
| Characteristic | Value (N=70) |
|---|---|
| Age, years | 70.5 (62–78) |
| Male sex | 51 (73%) |
| Lung cancer histology | |
| Adenocarcinoma | 55 (79%) |
| Squamous cell carcinoma | 15 (21%) |
| Pleural effusion side | |
| Right | 30 (43%) |
| Left | 33 (47%) |
| Bilateral | 7 (10%) |
| Pleural effusion≥½ hemithorax | 40 (57%) |
| Positive pleural fluid cytology | 47 (67%) |
| Malignant pleural effusion | |
| Definite | 49 (70%) |
| Probable | 21 (30%) |
| PD-L1 status | |
| Tissue, TPS, % (n=56) | 25 (0 to >50) |
| Pleural fluid, pg/mL (n=68) | 37 (0–56) |
| Plasma, pg/mL (n=61) | 0 (0–21) |
| Oncologic treatments | |
| Surgery | 2 (3%) |
| Chemotherapy | 42 (60%) |
| Immunotherapy | 30 (43%) |
| Targeted therapy | 11 (16%) |
| Follow-up, months | 5.5 (2–10) |
| Deaths | 67 (96%) |
Values are expressed as median (25–75th percentile) or No. (%) as appropriate. PD-L1, programmed death ligand-1; TPS, tumor proportion score.
Median cfDNA yield was significantly higher in PF than in plasma samples, with values of 8.82ng/mL (IQR, 3.38–28.03) and 1.64ng/mL (IQR, 1.15–3.04), respectively. Fragment-size analysis showed that both specimen types were enriched in mono-nucleosomal fragments (∼166bp), consistent with cfDNA profiles. The median proportion of mono-nucleosomal fragments was 47.11% (IQR, 36.87–55.15) in PF and 55.57% (IQR, 41.05–66.09) in plasma.
Detection of clinically relevant mutations in pleural fluid and plasmaClinically relevant mutations were detected in both PF and plasma cfDNA samples, with important differences according to specimen type and diagnostic category (Table 3), as illustrated in Fig. 1. Among patients with pathologically confirmed MPE, PF identified clinically relevant mutations in a higher proportion of patients than plasma (57% vs 41%). Overall, 67 distinct clinically relevant mutations were detected in this group, of which 40 were identified in PF and 27 in plasma, resulting in a significantly higher mutation yield for PF compared with plasma (P<.01).
Patients with clinically relevant mutations according to specimen type and diagnostic category.
| Definite MPE (n=49) | Probable MPE (n=21) | |||
|---|---|---|---|---|
| Pleural fluid | Plasma | Pleural fluid | Plasma | |
| Patients with ≥1 mutation, No. (%) | 28 (57%) | 20 (41%) | 6 (29%) | 5 (24%) |
| Patients with 1 mutation, No. (%) | 20 (41%) | 14 (29%) | 5 (24%) | 3 (14%) |
| Patients with ≥2 mutations, No. (%) | 8 (16%) | 6 (12%) | 1 (5%) | 2 (10%) |
| Total mutations detected, No. | 40 | 27 | 7 | 9 |
Clinically relevant mutations include targetable driver alterations and resistance-associated mutations. MPE, malignant pleural effusion.
Oncoprint representation of clinically relevant mutations detected in pleural fluid and plasma cfDNA. Each row represents an individual patient, and each column represents a gene. Mutations detected in pleural fluid and plasma cfDNA are indicated, allowing visualization of concordant and discordant alterations across sample types. Patients are stratified according to definite and probable malignant pleural effusion.
In patients with probable MPE, the prevalence of clinically relevant mutations was lower but not negligible. Mutations were detected in 29% of patients using PF and in 24% using plasma. In this subgroup, a total of 16 mutations were identified (7 in PF and 9 in plasma).
Across the entire cohort, PF overall detected a greater number of clinically relevant mutations than plasma (47 vs 36). EGFR alterations were the most frequently identified mutations, followed by KRAS variants. Mutations associated with resistance to targeted therapies, including EGFR T790M, EGFR C797S, and PIK3CA variants, were also detected. No ALK, MET, RET, ROS1, or NTRK fusions were identified in either specimen type (Table 4).
Spectrum of clinically relevant mutations detected in pleural fluid and plasma.
| Gene | Definite MPE (pleural fluid/plasma) | Probable MPE (pleural fluid/plasma) |
|---|---|---|
| EGFR | 27/21 | 3/4 |
| KRAS | 6/4 | 0/2 |
| ERBB2 | 1/1 | 0/1 |
| PIK3CA | 5/2 | 3/1 |
| BRAF | 1/0 | 0/0 |
| Total mutations | 40/27 | 7/9 |
Clinically relevant mutations include both targetable driver alterations and mutations associated with resistance to targeted therapies. Specific driver mutations detected included EGFR exon 7 (A289V), exon 18 (G719C), exon 19 deletions (E746_A750del and L747_A750delinsP) and point mutation (A750P), exon 21 mutations (L858R and L861Q); KRAS G12C; ERBB2 exon 20 insertion (Y772_A775dup); and BRAF V600E. Resistance-associated mutations identified in this study included EGFR exon 20 mutations (T790M and C797S) and PIK3CA exon 9 mutations (E542K and E545K). No ALK, MET, RET, ROS1, or NTRK gene fusions were identified.
MPE, malignant pleural effusion.
Variant allele frequency of clinically relevant mutations was significantly higher in PF cfDNA than in plasma cfDNA. Median VAF in PF was 0.136 (IQR, 0.021–0.465), compared with 0.038 (IQR, 0.017–0.121) in plasma (P<.001). This difference was consistently observed across mutation types and diagnostic categories.
Concordance between pleural fluid and plasma for mutation detectionOverall agreement between PF and plasma for the detection of clinically relevant mutations was low (Table 5). Among patients with pathologically confirmed MPE, concordant results between PF and plasma were observed in 26% of detected mutations. The corresponding Cohen κ coefficient was −0.49, indicating poor agreement beyond chance. In patients with probable MPE, no concordant mutations were identified between PF and plasma (0% agreement), with a Cohen κ coefficient of −0.97.
Concordance between pleural fluid and plasma for detection of clinically relevant mutations.
| Definite MPE (n=49) | Probable MPE (n=21) | |||
|---|---|---|---|---|
| Plasma: mutation detected | Plasma: no mutation detected | Plasma: mutation detected | Plasma: no mutation detected | |
| Pleural fluid: mutation detected | 14 | 26 | 0 | 7 |
| Pleural fluid: no mutation detected | 13 | 0 | 9 | 0 |
| Observed agreement, % | 26 | 0 | ||
| Cohen κ coefficient | −0.49 | −0.97 | ||
Clinically relevant mutations include both targetable driver alterations and mutations associated with resistance to targeted therapies. Concordance was assessed by comparing detection of each clinically relevant mutation within paired pleural fluid and plasma samples from the same patient.
MPE, malignant pleural effusion.
A subset of patients (n=14) had available formalin-fixed paraffin-embedded (FFPE) samples from PF cell blocks, pleural biopsies, or other tissue biopsies. In this subgroup, a total of 21 clinically relevant mutations were identified in FFPE samples. Overall concordance rates for detection of these mutations using liquid biopsies were higher for PF cfDNA than for plasma cfDNA (71% vs 43%). When stratified according to FFPE sample type, concordance rates with PF were 90% for cell blocks, 66% for pleural tissue, and 50% for other tissue biopsies. The corresponding concordance rates for plasma were 60%, 0%, and 38%, respectively (Supplementary Table S1).
Detection of clinically relevant mutations in probable malignant pleural effusionsAmong the 21 patients classified as having probable MPE, clinically relevant mutations were identified despite negative pleural cytology (Table 3), with patient-level data shown in Table 6. Using PF cfDNA analysis, at least 1 clinically relevant mutation was identified in 6 patients (29%), whereas plasma cfDNA analysis identified mutations in 5 patients (24%). Overall, a total of 7 mutations were detected in PF and 9 in plasma, involving genes associated with targetable alterations or resistance mechanisms, including EGFR, ERBB2, KRAS, and PIK3CA. When both minimally invasive specimens were considered together, clinically relevant mutations were identified in 11 of 21 patients (52%) with probable MPE.
Clinically relevant mutations detected in patients with probable malignant pleural effusion (n=21).
| Patient ID | Pleural fluid mutation(s) | Plasma mutation(s) |
|---|---|---|
| 1 | – | EGFR L861Q |
| 2 | – | EGFR T790M; EGFR C797S; EGFR L861Q; KRAS G12C |
| 3 | – | EGFR L861Q; ERBB2 Y772_A775dup |
| 4 | – | KRAS G12C |
| 5 | – | PIK3CA E542K |
| 6 | EGFR L861Q | – |
| 7 | EGFR L861Q | – |
| 8 | EGFR L861Q | – |
| 9 | EGFR G719C | – |
| 10 | PIK3CA E545K; PIK3CA E542K | – |
| 11 | PIK3CA E542K | – |
| 12–21 | – | – |
In this study, we demonstrated that PF cfDNA analysis detects clinically relevant mutations more frequently than plasma cfDNA in patients with lung cancer-associated MPE and yields higher VAFs, indicating a greater fraction of tumor-derived DNA within the pleural compartment. Concordance between PF and plasma was limited, indicating that both specimens are complementary rather than interchangeable. Importantly, clinically relevant mutations were identified in a substantial proportion of patients with probable MPE, underscoring the potential value of cfDNA analysis in this setting. Finally, although actionable alterations were frequently detected, only a limited proportion of patients received targeted therapy in routine clinical practice, revealing a gap between molecular detection and clinical implementation.
The proportion of patients in whom clinically relevant mutations were detected in PF cfDNA in our study (57% among confirmed MPEs) was lower than that reported in several previous series, in which detection rates approaching 80–95% have been described [3–7]. This discrepancy is likely explained by methodological rather than biological factors. Many prior studies reported all somatic alterations detected by NGS, including variants of uncertain significance, or applied sample-level selection criteria based on cfDNA yield, tumor fraction, or sequencing adequacy. In contrast, our analysis was performed on a per-patient basis and restricted to clinically relevant driver and resistance mutations, applying a more stringent and clinically oriented definition.
Despite differences in mutation yield, the degree of concordance between PF and plasma cfDNA observed in our cohort is consistent with previous reports [3–7]. Comparative studies have repeatedly shown incomplete overlap between specimens, with mutations detected exclusively in PF or plasma. In our study, only 26% of clinically relevant mutations were identified in both PF and plasma among patients with definite MPE, and no concordant mutations were observed in probable MPEs. These findings indicate that PF and plasma capture partially distinct molecular information. In a subset of patients with available tissue samples, concordance with FFPE-based molecular testing was higher for PF than for plasma, further supporting the greater sensitivity of PF cfDNA. In cases of discordant molecular findings between PF and other tumor sites, results should be interpreted within the context of spatial tumor heterogeneity and integrated into clinical decision-making rather than considered mutually exclusive. A minority of cases in our cohort reflected prior exposure to targeted therapies, as evidenced by the detection of resistance-associated mutations.
Evidence supporting the use of PF cfDNA in patients with cytology-negative pleural effusions remains limited and is derived mainly from small series. In most published studies, this subgroup comprised fewer than 10–15 patients and was included as an exploratory component within broader MPE cohorts [4,5,8,9]. In the study by Tong et al. [5,6], pleural effusions without detectable tumor cells on cytology were analyzed; PF cfDNA identified most tumor-derived mutations present in matched tissue, whereas pleural sediment DNA detected fewer than 10% of these alterations. Similarly, Mahmood et al. [4] reported oncogenic mutations in 7 cytology-negative effusions using PF cfDNA, highlighting that molecular alterations may be present despite the absence of cytologically identifiable malignant cells. More recently, Thakur et al. [10] demonstrated that cfDNA extracted from pleural effusion supernatant can reliably detect oncogenic mutations, including in cases with negative cytology. Against this background, our finding that clinically relevant mutations were identified in 11 of 21 patients (52%) with probable MPE provides additional quantitative support for the diagnostic value of PF cfDNA in this setting and represents one of the largest series focusing on clinically probable MPE.
Although clinically relevant mutations were detected in a substantial proportion of patients, only 16% received targeted therapy in our cohort, and all treated patients harbored sensitizing EGFR mutations. In contemporary series of advanced stage IV NSCLC, approximately 40–50% of patients are considered candidates for targeted therapies based on actionable oncogenic driver alterations [11]. This proportion is even higher in MPE secondary to NSCLC, with molecular analyses of pleural specimens consistently reporting actionable mutation rates ranging from 60% to 68%, largely driven by the high prevalence of EGFR mutations [12]. This range is consistent with our 57% detection rate among patients with pathologically confirmed MPE. The discrepancy between molecular eligibility and actual treatment in our cohort is explained mainly by the clinical context of the study period (2018–2022), during which comprehensive NGS-based profiling was not routinely implemented, analyses were frequently performed in research settings without real-time availability, and treatment decisions relied on limited EGFR testing.
Importantly, this finding should not be interpreted as a limitation of PF cfDNA performance. In the prospective study by Chang et al. [9], PF molecular analysis identified oncogenic driver mutations in 92.3% of patients, compared with 51.3% using standard diagnostic specimens, and cfDNA results were prospectively integrated into treatment selection, supporting the interpretation that the gap observed in our study reflects differences in clinical implementation rather than intrinsic limitations of PF-based cfDNA analysis. Identification of actionable mutations through liquid biopsy has clear therapeutic implications in advanced NSCLC. Current clinical guidelines support plasma-based genotyping when tissue is unavailable or insufficient, and our findings suggest that PF may represent an even more informative substrate in patients with MPE. In this context, liquid biopsy may facilitate timely initiation of targeted therapy and may be particularly valuable in patients who are not candidates for invasive procedures. However, negative results do not exclude targetable alterations, especially in settings of low tumor burden, and potential confounders such as clonal hematopoiesis and the need for appropriate variant filtering should be acknowledged when interpreting cfDNA results.
This study has several limitations. Although the cohort was prospectively enrolled, the analysis was conducted at a single center with a relatively limited sample size, which may affect generalizability. In addition, a proportion of pleural effusions were classified as probably malignant according to predefined clinical criteria rather than pathological confirmation, a strategy that reflects real-world clinical practice but may allow some degree of misclassification. Molecular analyses were performed in a research setting and were not systematically available to treating physicians at the time of clinical decision-making, precluding assessment of their direct impact on treatment selection. Finally, tissue-based molecular data were available only in a subset of patients and were analyzed in an exploratory manner; therefore, concordance results should be interpreted cautiously.
ConclusionsIn patients with lung cancer-associated pleural effusions, PF cfDNA analysis detects clinically relevant driver and resistance mutations more frequently than plasma and provides complementary molecular information. Clinically relevant mutations can be identified even in probable MPEs, underscoring the diagnostic value of PF-based molecular profiling. Systematic integration of PF cfDNA into routine clinical workflows has the potential to expand access to targeted therapies in advanced lung cancer [13].
CRediT authorship contribution statementJosé M. Porcel: Conceptualization, study design, data interpretation, and manuscript drafting. I. Hidalgo, M. Marqués, E. Parisi, A. Esquerda, and M.A. Sorolla: Sample processing and laboratory analysis. C. Farré and A. Salud: Clinical data collection. M.A. Sorolla: Data analysis. All authors reviewed and approved the final manuscript.
Declaration of generative AI and AI-assisted technologies in the writing processThe authors used Paperpal (Cactus Communications) as an AI-assisted tool for language editing and typographical corrections. In addition, generative artificial intelligence (ChatGPT, OpenAI) was used to assist in the design and visual development of the graphical abstract. All scientific content, study design, data analysis, interpretation of results, conclusions, and final approval of the manuscript and graphical abstract were performed exclusively by the authors.
FundingThis study was supported by grants from the Spanish Society of Pulmonology and Thoracic Surgery (SEPAR; Ref. 678/2018), Instituto de Salud Carlos III (PI20/01500 and PI23/00926), and the Agència de Gestió d’Ajuts Universitaris i de Recerca (AGAUR; 2021SGR00781).
Conflicts of interestNone declared.
This work was supported by the IRBLleida Biobank (B.0000682) and Plataforma Biobancos PT20/00021.















