Detection of early-stage lung cancer in sputum using automated flow cytometry and machine learning
Respiratory Research volume 24, Article number: 23 (2023)
Low-dose spiral computed tomography (LDCT) may not lead to a clear treatment path when small to intermediate-sized lung nodules are identified. We have combined flow cytometry and machine learning to develop a sputum-based test (CyPath Lung) that can assist physicians in decision-making in such cases.
Single cell suspensions prepared from induced sputum samples collected over three consecutive days were labeled with a viability dye to exclude dead cells, antibodies to distinguish cell types, and a porphyrin to label cancer-associated cells. The labeled cell suspension was run on a flow cytometer and the data collected. An analysis pipeline combining automated flow cytometry data processing with machine learning was developed to distinguish cancer from non-cancer samples from 150 patients at high risk of whom 28 had lung cancer. Flow data and patient features were evaluated to identify predictors of lung cancer. Random training and test sets were chosen to evaluate predictive variables iteratively until a robust model was identified. The final model was tested on a second, independent group of 32 samples, including six samples from patients diagnosed with lung cancer.
Automated analysis combined with machine learning resulted in a predictive model that achieved an area under the ROC curve (AUC) of 0.89 (95% CI 0.83–0.89). The sensitivity and specificity were 82% and 88%, respectively, and the negative and positive predictive values 96% and 61%, respectively. Importantly, the test was 92% sensitive and 87% specific in cases when nodules were < 20 mm (AUC of 0.94; 95% CI 0.89–0.99). Testing of the model on an independent second set of samples showed an AUC of 0.85 (95% CI 0.71–0.98) with an 83% sensitivity, 77% specificity, 95% negative predictive value and 45% positive predictive value. The model is robust to differences in sample processing and disease state.
CyPath Lung correctly classifies samples as cancer or non-cancer with high accuracy, including from participants at different disease stages and with nodules < 20 mm in diameter. This test is intended for use after lung cancer screening to improve early-stage lung cancer diagnosis.
Trial registration ClinicalTrials.gov ID: NCT03457415; March 7, 2018
Early detection of lung cancer through screening can increase survival and reduce morbidity [1, 2]. The US and regions of the UK recommend annual low-dose computed tomography (LDCT) screening for high-risk individuals . Although LDCT is very sensitive (93.8%) in detecting cancerous pulmonary nodules , its specificity is much lower (73.4%) because nodular images can be the result of various non-cancerous processes . A positive LDCT therefore requires follow-up tests to determine if the nodule is malignant . These medical procedures have inherent morbidity and mortality risks and can impose a serious burden on screening participants , while associated costs represent significant financial burdens to patients  and society [9, 10].
Efforts are underway to develop non-invasive tests that can be used after LDCT to improve screening’s predictive value [11, 12] or as stand-alone tests to identify people who should undergo screening [11, 13]. These tests aim to reduce unnecessary medical procedures while identifying those with lung cancer at an early stage. Sputum is easily accessible lung material that contains a variety of leukocytes and exfoliated bronchial epithelial cells , including premalignant and malignant cells in patients with lung cancer . We have previously reported on a slide-based microscopy assay that classified cancer and non-cancer patients using sputum stained with tetra(4-carboxyphenyl)porphyrin (TCPP) . Although 81% accurate, reading slides was time-consuming, subject to observer bias and could potentially miss key events by not evaluating the entire sample. We now report on a high-throughput approach using automated flow cytometry (FCM) instead of microscopy. FCM allows us to interrogate the entire sputum sample using TCPP and a panel of antibodies to capture cancer-predictive features in sputum. The field of automated FCM analysis has produced powerful software tools [17,18,19] that match or exceed human expertise in identifying cell populations of clinical importance . We adapted these tools to create an automated FCM sputum analysis platform, thereby eliminating potential operator bias .
Automated FCM techniques have been combined with machine learning approaches to distinguish leukemias from non-neoplastic cytopenias  and for biomarker discovery . We hypothesized that the same combination of approaches could be used to identify sputum features indicative for lung cancer in high-risk patients. Our objective was to develop an assay (referred to as CyPath Lung) that combines automated FCM data analysis of induced sputum samples with machine learning techniques to classify sputum samples as cancer or non-cancer. CyPath Lung is intended for use following screening to improve diagnosis of early-stage lung cancer.
Participants were identified and enrolled at a physician’s or study coordinator’s office at one of the following sites: Atlantic Health System, NJ; Mt. Sinai Hospital, NY; Radiology Associates of Albuquerque, NM; South Texas Veterans Healthcare System, TX; and Waterbury Pulmonary Associates, CT. Each site had received institutional approval to participate in the study. Samples were collected from April 2018 till November 2019 (LSRII sample set) and from July 2020 till November 2021 (Navios sample set).
Participants (males and females) were enrolled in one of two groups. The non-cancer group included participants (aged 52–79) who were either current smokers with a smoking history of at least 20 pack-years, or current non-smokers with a smoking history of at least 20 pack-years, who quit smoking within the past 15 years. The exceptions were two participants: one had quit smoking 26 years ago and one had smoked for 11.5 pack years. Most participants in the non-cancer group had received an LDCT result or other form of imaging that was not suspicious for cancer, and they were advised to return for LDCT screening in 12 months. In a few cases, participants initially placed in the non-cancer group underwent a follow-up LDCT, PET/CT or a biopsy. These participants were followed until their health status was confirmed. If they were diagnosed with lung cancer, they were switched to the cancer group.
Each participant in the cancer group had been evaluated by a physician as highly suspect of having lung cancer based on medical history and LDCT or other imaging results. The diagnosis was confirmed by biopsy after a sputum sample was provided. The exception was a patient who had developed a new nodule of 24 mm and who was too fragile to undergo biopsy. If biopsy showed no cancer, the participant was switched to the non-cancer group. There was no limitation of age or smoking history for enrollment in the cancer group.
For each participant we collected the following demographic data: gender (male or female); age (years); ethnicity (Hispanic/Latino or non-Hispanic/Latino); and race (American Indian/Alaska native; Asian; Black/African American; native Hawaiian/other Pacific islander; White; other). Data on smoking history was collected, as well as data on comorbidities (asthma, COPD, emphysema, chronic bronchitis) and previous cancer history. All participants needed to be willing to provide a primary care physician’s contact information and agree to have medical information released if requested. Exclusion criteria included the presence of severe obstructive lung disease and inability to cough with sufficient exertion to produce a sputum sample, angina with minimal exertion, and pregnancy.
Sample donors were trained on how to use the acapella assist device (Smiths Medical, St. Paul, MN), and expel their sample by coughing into a specimen cup, repeating this procedure at home for three consecutive days and storing their specimen cup in a refrigerator. Sample donors did not report experiencing any adverse events related to the sample collection procedure. Within one day after collection was completed, the sample was shipped overnight to the bioAffinity laboratory where further processing and FCM analysis took place by technicians blinded to the origin of the sample as well as the clinical information of the donor.
A set of 171 sputum samples was analyzed on the LSRII flow cytometer. One hundred and sixty-eight of the LSRII sample set were used for training and testing the model, as well as for developing the analysis pipeline. The final model was then validated on 150 LSRII samples that passed quality control. A second set of 45 samples was analyzed on the Navios EX. Thirty-two passed quality control and were used to independently measure the generality of the model/analysis pipeline by excluding a possible dependency on a particular flow cytometer or team of researchers. See Fig. 1 for more details.
This included four samples for which we did not have a definitive disease status because the addition of unlabeled samples had been shown to be helpful in model building .
Sample size considerations
Enrolment for the LSRII data set continued until sufficient cancer samples were collected to build a robust model for automated analysis. For test development and model training, we needed enough cancer samples in order to create subsets of samples through repeated randomization that would allow us to evaluate cancer predictor selection without repeatedly ending up with the same small number of cancer samples driving the model fitting. With 28 out of 150 samples being a cancer sample (~ 19%), we would be able to create > 3 million different training sets, each consisting of 100 samples (including 20 cancer samples, 20%). Three million different training sets are more than enough to test a wide range of potential parameters without a strong selection bias, while maintaining the cancer / non-cancer sample ratio of the entire sample set.
The Navios data set represents the number of samples recruited between establishing the analysis pipeline and the start of drafting this manuscript.
Sputum was dissociated and labeled using recently published protocols [25, 26]. Briefly, sputum samples were incubated with a mixture of 0.1% dithiothreitol and 0.5% N-acetyl-l-cysteine for 15 min at room temperature and neutralized with Hank’s Balanced Salt Solution (HBSS). Cells were then filtered through a 100-micron nylon strainer, washed and re-suspended in HBSS. Total cell yield was determined using trypan blue exclusion.
A small aliquot of cells was set aside for use as controls while the majority was divided into two tubes for the main analysis. Both tubes were labeled with Fixable Viability Stain 510 (FVS510) and CD45-PE. One tube, the “blood tube”, received CD66b-FITC (to identify granulocytes), CD3-Alexa-Fluor-488 (T-cells), CD19-Alexa-Fluor-488 (B-cells) and CD206-PE-CF594 (macrophages). In the other tube, the “epithelial tube”, cells were labeled with the epithelial cell markers pan-cytokeratin (Pan-CK)-Alexa-Fluor-488 and EpCAM-PE-CF594. Cells were incubated for 35 min on ice, protected from light. After washing with HBSS, cells were fixed and stored on ice until the next day, when a TCPP solution (20 µg/mL) was added (3.3 × 106 cells/ml; 1:1 v/v) for 1 h. Cells were washed twice and kept on ice and protected from light until analysis.
Sputum samples were acquired on a BD LSRII flow cytometer (BD Biosciences) equipped with four lasers (405 nm, 488 nm, 561 nm, and 633 nm) or on a Navios EX (Beckman Coulter Life Sciences) equipped with three lasers (405 nm, 488 nm and 638 nm). The settings used on each flow cytometer had been previously shown to generate similar sputum profiles .
Automated flow cytometric selection of viable single cells
The first stage of the CyPath Lung assay is the automated FCM identification of viable single cells (Fig. 2). The FCM component of the test consists of two assay tubes, one labelled with blood cell markers and one with epithelial cell markers . Cells in both tubes were also labeled with anti-CD45 antibodies which selectively bind leukocytes, a viability dye to eliminate dead cells, including squamous epithelial cells (SECs) , and TCPP to identify cancer or cancer-associated cells . Fluorescence intensities from antibody and TCPP staining were used exclusively for downstream numerical analysis.
Each sample run included polystyrene beads of known diameter (5–30 µm NIST beads), compensation tubes for each fluorochrome channel used, unstained sputum, and an antibody isotype sputum control. Each tube corresponded to a single Flow Cytometry Standard (FCS) file that contained sample metadata and per event values for each light and fluorescence channel acquired plus a Time parameter recorded as tubes are run.
As shown in Fig. 2A, the first step was to restrict events based on forward (FSC-A) and side (SSC-A) scatter area, both reasonable surrogates for cell size . A two-dimensional cluster gate identified the dominant peak of 5 µm NIST beads in FSC-A vs SSC-A (Fig. 2A(i)). The lower FSC-A limit of the bead cluster was set as the minimum sample FSC-A to exclude small particulates and debris. Upper limits of 2.5 × 105 were set on both FSC-A and SSC-A since events above those thresholds were found to be dead (FVS510+) cells (Fig. 2A and A(ii)).
Events within the bead size exclusion (BSE) gate were then restricted to exclude a population with unusual FSC and SSC height profiles (Fig. 2B) and a staining profile that might result in their inappropriate inclusion in downstream analyses (Fig. 2B(iii)). Exclusion of unusual-looking populations is warranted in general . In our case, it is important to exclude these spurious events to avoid including them in viable cell counts since the decision to include a sample for full analysis depends on total viable singlets.
Cells below the threshold for FVS510-positivity were retained (Fig. 2C). Heuristics based on subpopulations most likely to contain viable singlet cells (i.e., relatively small in area and height in light scatter channels) were used to guide its positioning (Fig. 3 and Additional file 1).
A “singlets” gate (FSC-Area vs FSC-Width) excluded cell doublets or small aggregates (Fig. 2D). In some samples, high SSC-A cells are included in the viable cell population and distort the results of the singlets gating algorithm. This can be corrected by fitting the gate to a temporary subpopulation excluding most of these high SSC-A events (Fig. 4A–D). In other samples, two populations can be seen within the viable cell gate, one with low FVS510 staining and the other just below the viability threshold and with a high side scatter profile (Fig. 4E–H). The correction involves resetting the viability gate and fitting the singlets gate to the more restricted viable population (see Additional file 1 for more details). Light scatter and fluorescence signal values were recorded for each single event and used for downstream model development and validation.
Development of the CyPath Lung cancer/non-cancer classifier
The second stage of the CyPath Lung assay is the analysis of light scatter and fluorescence signals from the viable single cells identified in the first stage by automated FCM. Logistic regression models describe a relationship between predictor variables and a categorical (in our case binary cancer/non-cancer) response variable. We measure our model’s performance by comparing its prediction of whether a sample is cancer or non-cancer (the response variable) to the known cancer/non-cancer status of the samples. Stepwise regression is a supervised machine learning process by which potentially predictive variables are added and removed and the resulting model examined for goodness of fit (see Additional file 1 for details). Demographic and clinical data (see Table 1) were included as potential predictors. Age was the only clinical parameter repeatedly rated as significant during forward and reverse stepwise regression.
Based on our earlier slide-based assay results, we anticipated that smoking history (or correlated factors like age) and TCPP signal density (as opposed to fluorescence intensity itself) would be important predictors . We therefore divided the fluorescence signals of all channels by log10 FSC-A or log10 SSC-A and partitioned the resulting density distribution into three regions (< 0.25, 0.25–0.6, > 0.6, Fig. 5A–D). Two such density signals proved informative for the classifier: TCPP/log10SSC-A (Fig. 5A, region 3 [R3]) and FVS510-A/log10FSC-A (Fig. 5C, region 2 [R2]). The predictive value of TCPP/log10SSC-A signal density was not imposed upon the stepwise regression but emerged spontaneously. The fact that FVS510-A signal density was also found to be informative is interesting and may be related to the fact that apoptotic cells can take up this dye at intermediate levels . No other single blood or epithelial fluorescence signal density was robustly and repeatedly identified as a predictor.
Combinations of lineage markers can identify subpopulations that single markers alone may not capture. Careful examination of patient samples revealed complex patterns of lineage marker expression in blood and epithelial tubes. Although there were differences between the non-cancer and cancer groups , we could not directly identify any subpopulation independently predictive of cancer by gating. Consequently, we decided to use a numerical approach to the analysis of pairwise markers by partitioning fluorescence based on signal distribution in blood (Fig. 5E, F) and epithelial tubes (data not shown). Signal intensity on the logicle scale was quantized into low (< 1.5), low-mid (1.5–2.5), mid (2.5–3), and high (> 3) windows. Events per 10,000 were tabulated for each area of the resulting 4 × 4 grid of CD206 vs CD3/CD19/CD66b (blood tube) or EpCAM vs Pan-CK (epithelial tube; not shown). One area in the blood signal intensity grid was found to be informative for the model (tan-shaded rectangle; low for CD206 and mid-level for CD3/CD19/CD66b, Fig. 5E). This population may indicate the presence of immune or inflammatory processes in the lung .
Once we had reduced the list of potential predictors to those that were most promising, we tested for pairwise interactions between them. One interaction proved informative: adding a negative value proportional to [age x number of events in FVS510-A/log10FSC-A R2] improved the classifier’s performance. Our interpretation of this interaction term is that it serves to moderate a possibly age-related accumulation of stressed cells in the non-cancer patient group as a consequence of smoking or health history .
Running the CyPath Lung assay pipeline
Having developed the two stages of the CyPath Lung assay, we could now assemble the full pipeline, including quality control steps, determination of predictive variable values, and classification of samples (Fig. 6).
Sample quality assessment begins by ensuring that the data file for each collection tube is readable and that its encoded data matrix is complete. Next, the Time signature is used to examine fluorescence channels in each tube and to remove anomalies in the flow rate arising from bubbles or clogs during sample acquisition . Fluorescence compensation tubes are used to derive the compensation matrix de novo (as opposed to using the compensation matrix encoded in the sample file metadata). Fluorescence signal is compensated and transformed to the logicle scale  to produce the sample data matrix used by the automated FCM gating to isolate viable singlet events. In order to have confidence in the downstream numerical analysis, given potentially small numbers of events in some analysis windows, we required that samples contain at least 10,000 viable singlets. We also required that at least 10 cells be present in the green-shaded area of Fig. 5E in which we find lung macrophages (CD206mid&highCD3/CD19/CD66blow−mid)  to confirm that the sputum sample originated in the lung.
The next step in the assay pipeline (Fig. 6) is to supply age in years and flow-based values from viable singlets to the classifier model. The likelihood of having cancer depends on four variables: age, number of events per 10,000 viable singlets (per 10K) with TCPP/log10SSC-A in Region 3 (Fig. 5A, R3), number of events per 10K with FVS510-A/log10FSC-A in Region 2 (Fig. 5C, R2), and number of events per 10K in the CD206lowCD3/CD19CD66bmid sector (Fig. 5E, tan shading). In addition, the model contains a negative term for the interaction between age and the FVS510 density variable and an “intercept” term (b0). The intercept term improves model fitting by not forcing the fitted line through 0 if all the variables are set to zero, but it is not directly interpretable as a biologically meaningful component of the classifier. The values of the coefficients (b1, b2, b3, b4, and b5) depend on the training set used for model fitting and provide weights for the variables. We did not normalize the data provided to the model in order to make interpretation of the model easier. For example, the model formula tells us that increasing the number of events with high TCPP density increases the likelihood of cancer, consistent with our previous results .
The final step of the pipeline (Fig. 6) is to make the cancer/non-cancer assignment. The model returns a value in the range of [0, 1]. Whether a given sample is classified as cancer depends on having the model return a value greater than a predetermined cutoff. If the value is less than or equal to the cutoff, the sample is classified as non-cancer. A rational cutoff can be selected by stepping through cutoff values between 0 and 1 and measuring the true positive and false positive calls as compared to the known group category at each step. Figure 7A shows the result of this process as a receiver operating characteristic (ROC) curve, with an area under the curve (AUC) of 0.89. The assay achieved its best performance at discriminating cancer from non-cancer at a threshold of 0.28 (Fig. 7B, solid vertical line).
Performance of CyPath Lung
We evaluated the performance of CyPath Lung for the 122 non-cancer and 28 cancer samples described in Table 1 and for an additional 32 samples (26 non-cancer and 6 cancer; Table 2) processed on a different FCM instrument (Navios EX). The same model with the same coefficients was used for both instruments, but the cutoff for the Navios samples was 0.5, based on the ROC curve for these samples. The results shown in Table 3 demonstrate that CyPath Lung performed very well with sensitivity, specificity, and accuracy all > 80% for the LSRII samples and very similar numbers for the smaller set of Navios EX samples. For both flow cytometry platforms we obtained a very robust negative predictive value (NPV) ≥ 95%.
The assay also performed remarkably well, with a sensitivity of 92% and specificity of 87% and an area under the ROC curve of 94%, if we restricted the analysis to cases where LDCT detected no nodules or only nodules < 20 mm in diameter (Table 3, “nodules all < 20 mm”). We do not consider the difference in the sensitivity and specificity between the full data set and the subset with nodules < 20 mm significant; however, it is evidence that the test performed equally well for difficult to diagnose individuals with smaller nodules. Furthermore, CyPath Lung performed well for all tumor types represented and at all disease stages, including I and II (Tables 4, 5).
Each of the retained predictors contributed significantly to the model (Wald Test P < 0.05) and removing them individually had a negative impact on the ability to correctly classify cancer and non-cancer samples (Table 6). Age is a well-established clinical correlate to lung cancer , as it is in our model; nevertheless, the correlation between age and the model value is not overwhelming in either LSRII or Navios EX samples (Fig. 8) with "cancer" called in some younger patients and "non-cancer" in many older ones. In fact, the exclusion of the CD206lowCD3/CD19CD66bmid signal resulted in as many misclassified samples as the exclusion of age and its interaction with FVS510-A/log10FSC-A R2 (Table 6).
To our knowledge, this study is the first that combines automated flow cytometric analysis with machine learning to predict the presence of lung cancer from sputum samples. Sputum as diagnostic material provides a snapshot of the tumor itself, of its microenvironment (ME), and of its field of cancerization (FoC). Expert cytological analysis of sputum can detect cancerous and pre-malignant cells , but it is a laborious approach that does not lend itself well to large-scale screening. Automated image processing has been used with some success to capture malignancy-associated changes in cells but is hampered by technical complexity and the low numbers of cells analyzed .
The case for moving to a high-throughput, automated flow-based approach combined with machine learning is thus compelling: (a) the assay can be put into routine lab use without requiring expert evaluation of samples or being subject to operator bias; (b) the entire sputum sample can be rapidly analyzed; and (c) numerical analysis can capture complex interactions between lung cancer, ME, and FoC cells which would be difficult for individuals to detect reliably. Our discovery during assay development of the predictive value of viability staining density, for example, suggests a link with apoptosis that merits further study. Our model also indicates that specific immune cell populations may be involved. Neither of these predictors was on our radar before we began model building, but in retrospect the importance of both processes early in cancer development is consistent with previous reports [42, 44].
Earlier versions of CyPath Lung relied solely on the porphyrin to distinguish cancer from non-cancer samples . The current automated flow cytometry-based test leverages viability staining and antibody profiling to capture additional aspects of tumorigenesis. One of the cancer predictors in CyPath Lung reflects an increase in immune cells in the cancer group. Since alterations in the immune system constitute an early response of the body to the presence of a tumor , it is possible that CyPath Lung can detect certain cancers before they are detectable by imaging. Others have shown that the performance of a sputum-based test for early lung cancer detection can significantly increase when different types of measurements are combined; for example, cytology with genetic mutations  or microRNAs and methylation biomarkers . Although we use one technology platform to measure different cancer-related processes, the additional parameters are likely contributing to the performance improvement from the slide-based assay to the flow cytometry-based assay (Fig. 7). Moreover, the flow cytometry-based assay reads the entire sample, which was also predicted to increase test performance .
Nearly 95% of participants in this study fulfilled the criteria for lung cancer screening most recently issued by the US Preventive Services Task Force . Although our study group can be considered a sample from those eligible for lung cancer screening (one of the target populations for CyPath Lung), the sampling was small with minorities being underrepresented, as were females in the cancer groups. Moreover, the cancer prevalence in our study was just below 19% for both data sets, which is considerably higher than in a lung cancer screening population  or in a patient group with lung nodules between 7 and 19 mm (the other target population for CyPath Lung) . Another limitation of our study is the lack of long-term follow-up of non-cancer participants to confirm they were indeed lung cancer-free. We intend to conduct a larger prospective clinical trial that addresses these limitations.
In its 2017 Official Policy Statement, the American Thoracic Society (ATS) stated that clinical usefulness of a novel biomarker should be evaluated by estimating the minimal accuracy required for that biomarker . A positive CyPath Lung test may help evaluate intermediate-sized lung nodules in LDCT-positive patients. The minimal accuracy ([sensitivity/(1-specificity] or positive diagnostic likelihood ratio (PDLR)) for CyPath Lung needs to be ≥ [(1-prevalence)/prevalence] × R/(1 − R)] according to the ATS statement. Assuming the threshold (R) – above which invasive follow up would be worthwhile - to be the frequency of cancer (4.8%) in the National Lung Screening Trial (NLST) population with intermediate nodules (7–19 mm in diameter) and using a cancer prevalence of 3.8% in the LDCT-positive population based on data of the NLST, we calculated the PDLR of CyPath Lung should be at least 1.28 [37, 50], which is a threshold met comfortably by our assay (Table 3).
The ATS statement also presents a use case where screening is expanded to include participants currently ineligible for LDCT screening . Using a hypothetical 1/500 prevalence of cancer and a harm threshold of 0.83%, a PDLR of 4.18 is estimated as the minimal accuracy for a useful test, a level met by the larger validation group of CyPath Lung (Table 3, LSRII). Using a hypothetical prevalence of 1/400 instead of 1/500 would yield a PDLR of 3.35, which both validation groups satisfy. When clinical utility is confirmed by future studies, CyPath Lung could serve to expand early lung cancer screening to relatively underserved populations such as younger females and male African American smokers [51, 52].
CyPath Lung is a non-invasive, sputum-based test for the early diagnosis of lung cancer. It uses a flow cytometric platform to analyze the cellular content of sputum with the analysis being fully automated and thus unbiased. The test is robust to differences in sample handling and processing and captures important predictive factors of early lung cancer carcinogenesis. The test performs well at 82% sensitivity and 88% specificity and achieves comparable performance when applied to an independent set of samples collected on a different flow cytometer. The test is also accurate in early stages (I and II) and in cases with nodules < 20 mm.
Availability of data and materials
Flow cytometry data files to reproduce the figures can be found in the FLOW repository (http://flowrepository.org) by accessing the following links: http://flowrepository.org/id/RvFrofFdPABAo3nGTsImVOuDSGrGDnt8PyAF8lc7PbLRmMTOWBbd7iDN0f9LalAG; http://flowrepository.org/id/RvFrTMVZGSin6hn2CH4oG2SZTPCTHni7Rs7vf5HxBzLH6s1sZdwhxDX4jXPxvFm7; http://flowrepository.org/id/RvFr9hmfgMXWnefKQbA9bkupqiEPy0DL6KAJxtpl7ZHnQzFWVU4LWDWojONbbkwD; http://flowrepository.org/id/RvFrF2ln2eH3W4IYYl9BOlaMhUlLd0zgsmxfwuhdzYhK67ve9iMwo1p7AYgUwZnJ. The entire data set used in this study is available from the corresponding author upon reasonable request.
Bead size exclusion
Epithelial cell adhesion molecule
Field of cancerization
Fixable Viability Stain 510
Hank’s Balanced Salt Solution
Squamous epithelial cells
National Lung Screening Trial Research Team, Aberle DR, Adams AM, Berg CD, Black WC, Clapp JD, et al. Reduced lung-cancer mortality with low-dose computed tomographic screening. N Engl J Med. 2011;365:395–409.
de Koning HJ, van der Aalst CM, de Jong PA, Scholten ET, Nackaerts K, Heuvelmans MA, et al. Reduced lung-cancer mortality with volume CT screening in a randomized trial. N Engl J Med. 2020;382:503–13. https://doi.org/10.1056/NEJMoa1911793.
Oudkerk M, Liu S, Heuvelmans MA, Walter JE, Field JK. Lung cancer LDCT screening and mortality reduction—evidence, pitfalls and future perspectives. Nat Rev Clin Oncol. 2021;18:135–51.
National Lung Screening Trial Research Team, Church TR, Black WC, Aberle DR, Berg CD, Clingan KL, et al. Results of initial low-dose computed tomographic screening for lung cancer. N Engl J Med. 2013;368:1980–91.
Seidelman JL, Myers JL, Quint LE. Incidental, subsolid pulmonary nodules at CT: etiology and management. Cancer Imaging. 2013;13:365–73.
Mazzone PJ, Lam L. Evaluating the patient with a pulmonary nodule: a review. JAMA. 2022;327:264–73.
Tanoue LT, Tanner NT, Gould MK, Silvestri GA. Lung cancer screening. Am J Respir Crit Care Med. 2015;191:19–33.
Tailor TD, Bell S, Fendrick AM, Carlos RC. Total and out-of-pocket costs of procedures after lung cancer screening in a national commercially insured population: estimating an episode of care. J Am Coll Radiol. 2022;19:35–46.
Black WC, Gareen IF, Soneji SS, Sicks JD, Keeler EB, Aberle DR, et al. Cost-effectiveness of CT screening in the National Lung Screening Trial. N Engl J Med. 2014;371:1793–802.
Huo J, Xu Y, Sheu T, Volk RJ, Shih Y-CT. Complication rates and downstream medical costs associated with invasive diagnostic procedures for lung abnormalities in the community setting. JAMA Intern Med. 2019;179:324–32.
Gaga M, Chorostowska-Wynimko J, Horváth I, Tammemagi MC, Shitrit D, Eisenberg VH, et al. Validation of lung EpiCheck, a novel methylation-based blood assay, for the detection of lung cancer in European and Chinese high-risk individuals. Eur Respir J. 2021;57:2002682.
Tanner NT, Springmeyer SC, Porter A, Jett JR, Mazzone P, Vachani A, et al. Assessment of integrated classifier’s ability to distinguish benign from malignant lung nodules: extended analyses and 2-year follow-up results of the PANOPTIC (Pulmonary Nodule Plasma Proteomic Classifier) trial. Chest. 2021;159:1283–7.
Sullivan FM, Mair FS, Anderson W, Armory P, Briggs A, Chew C, et al. Earlier diagnosis of lung cancer in a randomised trial of an autoantibody blood test followed by imaging. Eur Respir J. 2021;57:2000670.
Belda J, Leigh R, Parameswaran K, O’Byrne PM, Sears MR, Hargreave FE. Induced sputum cell counts in healthy adults. Am J Respir Crit Care Med. 2000;161:475–8.
Neumann T, Meyer M, Patten FW, Johnson FL, Erozan YS, Frable WJ, et al. Premalignant and malignant cells in sputum from lung cancer patients. Cancer. 2009;117:473–81.
Patriquin L, Merrick DT, Hill D, Holcomb RG, Lemieux ME, Bennett G, et al. Early detection of lung cancer with meso tetra (4-carboxyphenyl) porphyrin-labeled sputum. J Thorac Oncol. 2015;10:1311–8.
Montante S, Brinkman RR. Flow cytometry data analysis: recent tools and algorithms. Int J Lab Hematol. 2019;41(Suppl 1):56–62.
Pedreira CE, da Costa ES, Lecrevise Q, Grigore G, Fluxa R, Verde J, et al. From big flow cytometry datasets to smart diagnostic strategies: the EuroFlow approach. J Immunol Methods. 2019;475: 112631.
Wang S, Brinkman RR. Data-driven flow cytometry analysis. Methods Mol Biol. 2019;1989:245–65.
Conrad VK, Dubay CJ, Malek M, Brinkman RR, Koguchi Y, Redmond WL. Implementation and validation of an automated flow cytometry analysis pipeline for human immune profiling. Cytometry A. 2019;95:183–91. https://doi.org/10.1002/cyto.a.23664.
Brinkman RR. Improving the rigor and reproducibility of flow cytometry-based clinical research and trials through automated data analysis. Cytometry A. 2020;97:107–12.
Monaghan SA, Li J-L, Liu Y-C, Ko M-Y, Boyiadzis M, Chang T-Y, et al. A machine learning approach to the classification of acute leukemias and distinction from nonneoplastic cytopenias using flow cytometry data. Am J Clin Pathol. 2021;157: aqab148.
Botta C, Maia CDS, Garcés J-J, Termini R, Perez C, Manrique I, et al. FlowCT for the analysis of large immunophenotypic datasets and biomarker discovery in cancer immunology. Blood Adv. 2022;6(2):690–703.
Gibson BR, Rogers TT, Zhu X. Human semi-supervised learning. Top. Cogn Sci. 2013;5:132–72.
Grayson M, Lai S-C, Bederka LH, Araujo P, Sanchez J, Reveles XT, et al. Quality-controlled sputum analysis by flow cytometry. J Vis Exp. 2021. https://doi.org/10.3791/62785.
Bederka LH, Sanchez JR, Rebeles J, Araujo PR, Grayson MH, Lai S-C, et al. Sputum analysis by flow cytometry; an effective platform to analyze the lung environment. PLoS ONE. 2022;17: e0272069.
Kasai Y, Sugiyama H, Takagi R, Kondo M, Owaki T, Namiki H, et al. Brush biopsy of human oral mucosal epithelial cells as a quality control of the cell source for fabrication of transplantable epithelial cell sheets for regenerative medicine. Regen Ther. 2016;4:71–7.
El-Far M, Pimstone N. A comparative study of 28 porphyrins and their abilities to localize in mammary mouse carcinoma: uroporphyrin I superior to hematoporphyrin derivative. Prog Clin Biol Res. 1984;170:661–72.
Tzur A, Moore JK, Jorgensen P, Shapiro HM, Kirschner MW. Optimizing optical flow cytometry for cell volume-based sorting and analysis. PLoS ONE. 2011;6: e16053.
Cossarizza A, Chang H-D, Radbruch A, Acs A, Adam D, Adam-Klages S, et al. Guidelines for the use of flow cytometry and cell sorting in immunological studies (second edition). Eur J Immunol. 2019;49:1457–973.
Altorki NK, Markowitz GJ, Gao D, Port JL, Saxena A, Stiles B, et al. The lung microenvironment: an important regulator of tumour growth and metastasis. Nat Rev Cancer. 2019;19:9–31.
Bowdish DME. The aging lung: is lung health good health for older adults? Chest. 2019;155:391–400.
Watson JV. Time, a quality-control parameter in flow cytometry. Cytometry. 1987;8:646–9.
Parks DR, Roederer M, Moore WA. A new “Logicle” display method avoids deceptive effects of logarithmic scaling for low signals and compensated data. Cytometry A. 2006;69:541–51.
Liegeois M, Legrand C, Desmet CJ, Marichal T, Bureau F. The interstitial macrophage: a long-neglected piece in the puzzle of lung immunity. Cell Immunol. 2018;330:91–6.
Pepe MS, Janes H, Li CI, Bossuyt PM, Feng Z, Hilden J. Early-phase studies of biomarkers: what target sensitivity and specificity values might confer clinical utility? Clin Chem. 2016;62:737–42.
Signorell A, et mult. al. DescTools: tools for descriptive statistics. R package version 0.99.47. 2022. https://cran.r-project.org/package=DescTools. Accessed 26 Jan 2023
Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez J-C, et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.
Schaarschmidt F. bdpv: Inference and Design for Predictive Values in Diagnostic Tests. 2019. https://CRAN.R-project.org/package=bdpv. Accessed 18 Feb 2022
Mercaldo ND, Lau KF, Zhou XH. Confidence intervals for predictive values with an emphasis to case-control studies. Stat Med. 2007;26:2170–83.
White MC, Holman DM, Boehm JE, Peipins LA, Grossman M, Henley SJ. Age and cancer risk: a potentially modifiable relationship. Am J Prev Med. 2014;46:S7-15.
Li G, Guillaud M, LeRiche J, McWilliams A, Gazdar A, Lam S, et al. Automated sputum cytometry for detection of Intraepithelial neoplasias in the lung. Anal Cell Pathol. 2012;35:187–201.
Trinchieri G. Cancer and inflammation: an old intuition with rapidly evolving new concepts. Annu Rev Immunol. 2012;30:677–706.
Gonzalez H, Hagerling C, Werb Z. Roles of the immune system in cancer: from tumor initiation to metastatic progression. Genes Dev. 2018;32:1267–84.
Katz RL, Zaidi TM, Fernandez RL, Zhang J, He W, Acosta C, et al. Automated detection of genetic abnormalities combined with cytology in sputum is a sensitive predictor of lung cancer. Mod Pathol. 2008;21:950–60.
Su Y, Fang H, Jiang F. Integrating DNA methylation and microRNA biomarkers in sputum for lung cancer detection. Clin Epigenet. 2016;8:109.
US Preventive Services Task Force, Krist AH, Davidson KW, Mangione CM, Barry MJ, Cabana M, et al. Screening for lung cancer: US Preventive Services Task Force recommendation statement. JAMA. 2021;325:962–70.
Gierada DS, Pinsky P, Nath H, Chiles C, Duan F, Aberle DR. Projected outcomes using different nodule sizes to define a positive CT lung cancer screening examination. J Natl Cancer Inst. 2014;106: dju284.
Mazzone PJ, Sears CR, Arenberg DA, Gaga M, Gould MK, Massion PP, et al. Evaluating molecular biomarkers for the early detection of lung cancer: when is a biomarker ready for clinical use? An official American Thoracic Society policy statement. Am J Respir Crit Care Med. 2017;196:e15-29.
MacRosty CR, Rivera MP. Lung cancer in women: a modern epidemic. Clin Chest Med. 2020;41:53–65.
Prosper A, Brown K, Schussel B, Aberle D. Lung cancer screening in African Americans: the time to act is now. Radiol Imaging Cancer. 2020;2: e200107.
We would like to thank all sputum donors for their participation in this study and David Rodriguez for his assistance with figure preparations. Flow cytometric data acquisition was performed at the University of Texas Health San Antonio flow cytometry core facility (LSRII), which is supported by NIH-NCI P30 CA054174-20 and UL1 TR001120 grants, and at Precision Pathology Services (Navios EX).
This research was funded by bioAffinity Technologies.
Ethics approval and consent to participate
The study described herein is part of a minimal risk protocol registered with ClinicalTrials.gov (NCT03457415), reviewed and approved by the Sterling Institutional Review Board (Atlanta, GA), and conducted according to ethical principles of the Declaration of Helsinki (v 1996) and Good Clinical Practice guidelines. Participants were enrolled in the study upon signing a written consent after the study was explained to them by a study coordinator.
Consent for publication
MEL received payment from bioAffinity Technologies for services provided. XTR, JR, LHB, PRA, JS, MG, S-CL and VIR are current or past employees of bioAffinity. The respective affiliated hospitals and medical facilities of LRD, SAH, DGH, KL and LP received payment from bioAffinity for patient sample collection LRD, SAH, DGH, KL and LP performed. LP has consulted for bioAffinity on matters related to radiology. RPJ is the owner of Precision Pathology Services, a commercial collaborator of bioAffinity and holds stock in bioAffinity.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Lemieux, M.E., Reveles, X.T., Rebeles, J. et al. Detection of early-stage lung cancer in sputum using automated flow cytometry and machine learning. Respir Res 24, 23 (2023). https://doi.org/10.1186/s12931-023-02327-3