Abstract
The epithelial-to-mesenchymal transition, the major process by which some cancer cells convert from an epithelial phenotype to a mesenchymal one, has been suggested to drive chemo-resistance and/or metastasis in patients with cancer. However, only a few studies have demonstrated the presence of CD45/CD326 doubly-positive cells (CD45/CD326 DPC) in cancer. We deployed a combination of cell surface markers to elucidate the phenotypic heterogeneity in non–small cell lung cancer (NSCLC) cells and identified a new subpopulation that is doubly-positive for epithelial and non–epithelial cell-surface markers in both NSCLC cells and patients' malignant pleural effusions.
We procured a total of 39 patients' samples, solid fresh lung cancer tissues from 21 patients and malignant pleural effusion samples from 18 others, and used FACS and fluorescence microscopy to check their surface markers. We also examined the EGFR mutations in patients with known acquired EGFR mutations.
Our data revealed that 0.4% to 17.9% of the solid tumor tissue cells and a higher percentage of malignant pleural effusion cells harbored CD45/CD326 DPC expressing both epithelial and nonepithelial surface markers. We selected 3 EGFR mutation patients and genetically confirmed that the newly identified cell population really originated from cancer cells. We also found that higher proportions of CD45/CD326 DPC are significantly associated with poor prognosis.
In conclusion, varying percentages of CD45/CD326 DPC exist in both solid cancer tissue and malignant pleural effusion in patients with NSCLC. This CD45/CD326 doubly-positive subpopulation can be an important key to clinical management of patients with NSCLC.
In this work, we report the identification of a novel cell population, CD45/CD326 doubly-positive cells (CD45/CD326 DPC), in non–small cell lung cancer (NSCLC), both in primary tumors and in metastatic pleural effusions. To elucidate the origin of these cells, we analyzed mutations of typical driver genes and found that these cells harbored the same mutations as those in corresponding primary tumors; these cells are thought to originate from the primary tumors. Moreover, patients with higher proportions of CD45/CD326 DPC are significantly associated with poor prognosis. Our results will shed light on the yet-to-be-discovered missing link in cancer development, and also give a valuable clue on clinical management of patients with NSCLC.
Introduction
Lung cancer is one of the most common types of malignant tumor in the world (1, 2), and non–small cell lung cancer (NSCLC) accounts for around 80% of lung cancer cases (3, 4). Many studies have shown that cancer cells are heterogeneous; such intratumoral diversities may help tumor cells to survive against immune system surveillance and a variety of treatment methods (5–7). A malignant effusion from patients with cancer is considered to be a poor prognostic sign, and tumor heterogeneity is also an important clinical indicator of patient outcomes. For example, higher tumor cell ratios in ascites from gastric cancer patients can contribute to poorer patient outcomes (8). The epithelial-to-mesenchymal transition (EMT) in cancer is known to be associated with therapeutic resistance, aggressiveness, and metastasis (4–6, 9, 10). However, the existence of doubly-positive for epithelial and nonepithelial surface markers, has not previously been demonstrated. Whether or not such CD45/CD326 doubly-positive cancer cells really exist is one of the most important questions in current oncology. If CD45/CD326 doubly-positive cells (CD45/CD326 DPC) can be identified in cancers, they would be promising targets for future research to elucidate the therapeutic resistance, aggressiveness, and metastatic ability of NSCLC.
In this study, we analyzed tissues from 39 patients with NSCLC. Twenty-one were samples of primary solid tumor tissues and 18 were samples of malignant pleural effusion. We attempted to identify CD45/CD326 doubly-positive cancer cells by FACS, fluorescence microscopy, and EGFR mutation analyses. The possible association between the frequency of CD45/CD326 DPC and poor prognosis was also examined. Our present results will contribute to understanding the effect of CD45/CD326 doubly-positive cancer cells on patient outcomes.
Materials and Methods
Sample collection
Samples from 39 patients with NSCLC were collected between 2010 and 2017 with written informed consent in accordance with the Declaration of Helsinki. The study was approved by the Ethics Committee of the Tohoku University School of Medicine under the accession number 2016-1-452. Tumors were procured within 2 to 3 hours after excision. Resected samples were kept in DMEM/F12 medium supplemented with 2% FBS (Thermo Fisher Scientific) for single cell dissociation. Solid tumors were minced into small pieces which could pass through a 5 mL pipette tip and digested with collagenase 1A solution (Stem Cell Technologies) in a 50 mL tube at 37°C on a shaker at 300 rpm agitation for 1 to 2 hours. The collagenase 1A was inactivated by addition of an equal amount of Iscove's modified Dulbecco's medium (IMDM; Thermo Fisher Scientific) supplemented with 5% FBS. Red blood cells were lysed in 1× ammonium–chloride–potassium (ACK) buffer (pH 7.2). As summarized in Supplementary Table S1, total of 39 patient samples, surgically resected cancer tissues from 21 patients and malignant pleural effusion samples from another 18 patients were used in this study.
Cell sorting for single-cell preparation
Single-cell suspensions were obtained by 2-step sequential filtration using a 100- and a 40-μm mesh and then counted. After the cell preparation, living cells were collected by FACSAria II (BD Bioscience). Debris was excluded by proper gating, and dead cells were excluded with propidium iodide (PI). The remaining living cells were further sorted into single cells and clumping cells by FACSAria II. The clumping cells were detected in groups of varying sizes by using forward scatter (FSC-W) and (FSC-H) and side scatter (SSC-W) and (SSC-H) sequentially. After these sorting processes, single cells without debris and clumping cells were separately purified for detection of genetic alterations and cell surface marker expression.
Epithelial and mesenchymal cell surface markers in FACS
In this study, CD326 (Miltenyi Biotec), an epithelial cell adhesion molecule, was used as the cell surface marker for epithelial cells (11). CD45 (Miltenyi Biotec), a leukocyte cell surface glycoprotein, was used as the surface marker for mesenchymal cells (12). The fluorescence-labeled antibodies against these surface markers (CD326 PE 130-113-264 1:50, CD45 FITC 130-113-17; 1:10) were incubated in the dark with samples at 4°C for 30 minutes, washed 3 times with FACS buffer, and filtered through a 35-μm mesh before the FACS analyses. Single cells from cancer samples were sorted according to the positivity of CD45 and CD326 on their stained cell surfaces and categorized into the following groups: (i) CD45neg/CD326neg cells, (ii) CD45neg/CD326pos epithelial cancer cells, (iii) CD45pos/CD326neg mesenchymal cells, and (iv) CD45pos/CD326pos cells. After sorting, the cells were collected and observed by fluorescence microscopy.
Immuno-fluorescence microscopy to detect CD45pos/CD326pos cells
Immuno-fluorescence microscopy used a N-SIM (Nikon Instech Co.), anti-pan Cytokeratin-antibody (KL1), which is commonly used for detection of epithelial cells (1:500; Cosmobio Ltd.), and anti-CD45-antibody (1:200) (CST 13917S, Cell Signaling Technology) to visually confirm the existence of CD45pos/CD326pos cancer cells in solid cancer samples. Two secondary antibodies were used: Alexa Fluor 488, a goat anti-mouse IgG antibody that emits a green color, was tagged on anti-KL1-antibody, (1:1,000, Fisher A-11001), and CY3 conjugated goat anti rabbit IgG antibody that emits a red color was tagged on anti-CD45-antibody (1:1,000, Abcam ab6939). Detection of cancer cells in pleural effusion employed 2 antibodies also used in the FACS analyses; PE-labeled anti-CD326-antibody (red) and FITC-labeled anti-CD45-antibody (green). Because merged immunofluorescent images of red and green colors result in yellow emission light, any CD45pos/CD326pos cells in the cancer samples will appear yellow in this protocol. 4′,6-Diamidino-2-phenylindole (DAPI) was used for nucleic acid staining to confirm that the substances bound by these labeled antibodies were actually living cells.
Detection of the EGFR gene mutations in CD45/CD326 DPC
To confirm that the CD45pos/CD326pos cells are truly cancer cells that harbor cancer-causing EGFR mutations, we performed the following analyses. Cells from the pleural effusion were categorized into the above-mentioned 4 groups, and their DNAs were purified using NucleoSpin Tissue XS (MACHEREY-NAGEL; GmbH & Co. KG), according to supplier's recommendation. Target EGFR regions were amplified by the GeneAmp PCR system 9700 [Applied Biosystems (13)]. Nucleotide sequences were determined by methods described previously (14) using an ABI PRISM310 Genetic Analyzer (Applied Biosystems). Because the EGFR mutations in our patients were all indel type, we could also simply have confirmed the mutations by polyacrylamide gel electrophoreses. Nucleotide sequences of the primers and PCR conditions are summarized in Supplementary Table S2.
Results
Confirmation of CD45pos/CD326pos cells by FACS
We detected a wide range of proportions of CD45pos/CD326pos single cells in our samples by FACS (Fig. 1A). These single cells were strong candidates for EMT, so we investigated them further. They were detected in both solid cancer tissue (Fig. 1B) and malignant pleural effusion (Fig. 1C). A serum specimen extracted from a healthy middle-aged male volunteer was used as the control; no CD45/CD326 DPC were detected in it (Fig. 1A). The average size of the CD45/CD326 DPC population was 5.5% (0.4%–17.9%) in solid tumors and nearly 7 times higher, 34.5% (1.1%–96.7%), in malignant pleural effusions.
Detection of CD45pos/CD326pos cells by FACS analysis in both solid tumor and pleural effusion. A, Typical examples of FACS analyses are shown. Results of 4 solid tumor tissue samples are displayed; various percentages of CD45pos/CD326poscells, indicated by red arrows, were detected. Such cells were not detected at all in the serum from a healthy male (blue arrow). In each quadrant, the number of cells and percentage are indicated at the top and bottom, respectively. B, Percentages of CD45pos (candidate mesenchymal) cells, CD326pos (candidate epithelial cancer) cells, and doubly-positive (CD45pos/CD326pos) cells (blue dots, 0.4%–17.9%) of the 21 analyzed solid tumors are shown. C, Percentages of CD45pos (candidate mesenchymal) cells, CD326pos (candidate epithelial cancer) cells, and doubly-positive (CD45pos/CD326pos) cells (blue dots, 1.1–96.7%) from 18 analyzed pleural effusions are shown.
Detection of CD45pos/CD326pos cells by FACS analysis in both solid tumor and pleural effusion. A, Typical examples of FACS analyses are shown. Results of 4 solid tumor tissue samples are displayed; various percentages of CD45pos/CD326poscells, indicated by red arrows, were detected. Such cells were not detected at all in the serum from a healthy male (blue arrow). In each quadrant, the number of cells and percentage are indicated at the top and bottom, respectively. B, Percentages of CD45pos (candidate mesenchymal) cells, CD326pos (candidate epithelial cancer) cells, and doubly-positive (CD45pos/CD326pos) cells (blue dots, 0.4%–17.9%) of the 21 analyzed solid tumors are shown. C, Percentages of CD45pos (candidate mesenchymal) cells, CD326pos (candidate epithelial cancer) cells, and doubly-positive (CD45pos/CD326pos) cells (blue dots, 1.1–96.7%) from 18 analyzed pleural effusions are shown.
Confirmation of CD45/CD326 doubly-positive cancer cells by fluorescence microscope and Papanicolaou staining
Using KL1 and CD326 antibodies for epithelial cell markers and an immuno-fluorescence microscope, we visually confirmed the existence of CD45/CD326 DPC (CD45pos/CD326pos cancer cells by FACS analyses) in samples with solid tumor tissue or malignant pleural effusion. Figure 2A illustrates immune-fluorescence microscopic images of an adenocarcinoma sample (T1aN0M0; stage IB). In the merged image, the KL-1pos/CD45pos doubly-positive cancer cells are clearly present (yellow arrows). The fluorescence of sorted pleural effusion cells from 4 quadrants (Q1–Q4) are shown in Fig. 2B; doubly-positive cells are present only in Q2, not in the other quadrants. Papanicolaou staining was done after sorting the malignant pleural effusion, and the sorted cells in each subpopulation were microscopically observed for positivity of CD45 and CD326 in Q1 through Q4. Q1 and Q2 contained malignant cells with large and irregular nuclei indicated by arrowheads. Based on these findings, CD45pos/CD326pos cells are considered to be CD45/CD326 doubly-positive malignant cells.
Immuno-fluorescence microscopy and Papanicolaou staining used to detect CD45pos/CD326pos doubly-positive cancer cells in solid tumors and pleural effusions. A, Immuno-fluorescence microscopic images of a solid cancer tissue sample from a patient with NSCLC are shown. DAPI staining indicates a nucleus (blue color, top left); KL-1 Alexa fluor-488 staining indicates an epithelial cell (green color, top right); CD45 Cy3 indicates a mesenchymal cell (red color, bottom left); and the merged image is bottom right. Arrows indicate KL-1pos/CD45pos doubly positive cells; candidate CD45/CD326 DPC. B, Results of immuno-fluorescence microscopic images of cells in pleural effusion from a patient with NSCLC. After gating by positivity of CD326 and/or CD45 expression, cells were divided into 4 groups; Q1 to Q4 (top center). Q1 column indicates CD45neg/CD326pos cells which express CD326 (red) but are negative for CD45 (green). Q2 column indicates CD45pos/CD326pos cells which express both CD326 (red) and CD45 (green). Yellow cells (doubly-positive for CD326 and CD45, i.e., CD45/CD326 doubly-positive cancer cells) in the merged images were only observed in Q2. Q3 and Q4 indicate non-epithelial cells that are negative for CD326 (red). C, Papanicolaou staining images of each group after cell sorting. As is seen with CD45neg/CD326pos epithelial cancer cells in Q1, CD45pos/CD326pos doubly-positive cancer cells also show malignant characteristics with large nucleus-to-cytoplasm ratios, as shown in Q2. On the other hand, CD326-negative groups (CD45pos/CD326neg in Q3 and CD45neg/CD326neg in Q4) do not show such malignant characteristics.
Immuno-fluorescence microscopy and Papanicolaou staining used to detect CD45pos/CD326pos doubly-positive cancer cells in solid tumors and pleural effusions. A, Immuno-fluorescence microscopic images of a solid cancer tissue sample from a patient with NSCLC are shown. DAPI staining indicates a nucleus (blue color, top left); KL-1 Alexa fluor-488 staining indicates an epithelial cell (green color, top right); CD45 Cy3 indicates a mesenchymal cell (red color, bottom left); and the merged image is bottom right. Arrows indicate KL-1pos/CD45pos doubly positive cells; candidate CD45/CD326 DPC. B, Results of immuno-fluorescence microscopic images of cells in pleural effusion from a patient with NSCLC. After gating by positivity of CD326 and/or CD45 expression, cells were divided into 4 groups; Q1 to Q4 (top center). Q1 column indicates CD45neg/CD326pos cells which express CD326 (red) but are negative for CD45 (green). Q2 column indicates CD45pos/CD326pos cells which express both CD326 (red) and CD45 (green). Yellow cells (doubly-positive for CD326 and CD45, i.e., CD45/CD326 doubly-positive cancer cells) in the merged images were only observed in Q2. Q3 and Q4 indicate non-epithelial cells that are negative for CD326 (red). C, Papanicolaou staining images of each group after cell sorting. As is seen with CD45neg/CD326pos epithelial cancer cells in Q1, CD45pos/CD326pos doubly-positive cancer cells also show malignant characteristics with large nucleus-to-cytoplasm ratios, as shown in Q2. On the other hand, CD326-negative groups (CD45pos/CD326neg in Q3 and CD45neg/CD326neg in Q4) do not show such malignant characteristics.
Confirmation of the EGFR gene mutations in CD45/CD326 DPC
Mutations of the EGFR gene were analyzed in all the 18 patients with pleural effusion. Five of them had disease-causing mutations, all of which were indel in either exon 19 or 20. We PCR-amplified target mutant regions in 3 of these 5 patients followed by 10% polyacrylamide gel electrophoresis. As shown in Fig. 3A, every DNA sample from CD45pos/CD326pos (doubly-positive) cells showed the same mutation observed in its corresponding patient. Notably, a strong upper mutant band was observed, but only a weak normal-sized band was detected in DNAs from Q1 and Q2 in Case #1; this strongly suggested a homozygous mutation in the cancer cells. Based on these results, the CD45/CD326 DPC were confirmed to be malignant from the genetic point of view. Nucleotide sequencing analyses were also performed; a typical example is shown in Fig. 3B. A 9-base insertion in exon 20 of EGFR, observed in Case #1, was detected only in Q1 and Q2, the CD326pos fraction. In contrast, the CD326neg group did not have such mutations. We summarized all the mutations detected in our study in Table S3. In addition, mutations of KRAS and TP53 as the typical driver genes for NSCLC were analyzed as shown in Supplementary Fig. S1.
Results of genetic analyses of EGFR. A, Polyacrylamide gel electrophoreses of PCR-amplified products show aberrations. The same 211-bp normal-sized band shown in WT is observed in 2 types of CD45neg/CD326pos and CD45neg/CD326neg cells but 9-bp larger bands are also obvious in lanes with CD45pos/CD326pos and CD45pos/CD326neg cells, as indicated by an arrow in samples from case #1. In the same manner, 9-bp shorter than normal-sized 146-bp bands in case #6 and 12-bp shorter bands in case #9 were observed only in lanes with CD45pos/CD326pos and CD45pos/CD326neg cells. The mutated bands are indicated by arrows. PC-9, which harbors a 15-base deletion in exon 19, was used as a mutant control. WT denotes the normal control, DNA from a healthy volunteer. B, Nucleotide sequencing results of case #1: a 9-base insertion (CCAGCGTGG) in EGFR exon 20 (yellow shaded) was clearly observed in the both CD45neg/CD326pos cells and CD45pos/CD326pos cells. The schema above indicates the insertion. No such mutant sequence was observed in CD45pos/CD326neg or CD45neg/CD326neg cells in the nonepithelial fraction, from the same patient.
Results of genetic analyses of EGFR. A, Polyacrylamide gel electrophoreses of PCR-amplified products show aberrations. The same 211-bp normal-sized band shown in WT is observed in 2 types of CD45neg/CD326pos and CD45neg/CD326neg cells but 9-bp larger bands are also obvious in lanes with CD45pos/CD326pos and CD45pos/CD326neg cells, as indicated by an arrow in samples from case #1. In the same manner, 9-bp shorter than normal-sized 146-bp bands in case #6 and 12-bp shorter bands in case #9 were observed only in lanes with CD45pos/CD326pos and CD45pos/CD326neg cells. The mutated bands are indicated by arrows. PC-9, which harbors a 15-base deletion in exon 19, was used as a mutant control. WT denotes the normal control, DNA from a healthy volunteer. B, Nucleotide sequencing results of case #1: a 9-base insertion (CCAGCGTGG) in EGFR exon 20 (yellow shaded) was clearly observed in the both CD45neg/CD326pos cells and CD45pos/CD326pos cells. The schema above indicates the insertion. No such mutant sequence was observed in CD45pos/CD326neg or CD45neg/CD326neg cells in the nonepithelial fraction, from the same patient.
CD45pos/CD326pos expression correlates with poor survival of the patients
Next, we examined the CD45pos/CD326pos cells impact on patient outcomes in the 18 patients with pleural effusion. The proportion of CD45pos/CD326pos cells in their entire CD326pos cell population were calculated and their relationship to patients' prognoses was plotted (Fig. 4). Two patients had over 90% of CD45pos/CD326pos cells, although number of cancer cells we could have analyzed was not many (data not shown). Four of the 18 (22%) patients remain alive during this cohort study. Three of these 4 survivors have the EGFR mutation and are being treated with tyrosine kinase inhibitor (TKI); the remaining survivor is being treated with PD-1. We further investigated whether there were any associations between CD45pos/CD326pos cell proportions and patients' prognoses. We divided the patients into 2 groups based upon their average CD45pos/CD326pos cell proportion and found that patients with lower proportions had better prognoses (P = 0.0213 by χ2 test, P = 0.0686 by Fisher exact test).
CD45pos/CD326pos expression may associate with poor survival of patients. Scatter plot of 18 patients with pleural effusion is shown. Single cells were isolated within 24 hours from the collection of pleural effusion. Y-axis indicates the percentage of the CD45pos/CD326pos cells in the whole CD326pos cell sample. X-axis indicates the days after the initial paracentesis. Patients labeled as alive indicate those who were alive at the last data collection; all other patients had died. Among the 4 surviving patients, 1 is being treated with PD-1, and 3 others with TKI. Ins indicates insertion and del indicates deletion.
CD45pos/CD326pos expression may associate with poor survival of patients. Scatter plot of 18 patients with pleural effusion is shown. Single cells were isolated within 24 hours from the collection of pleural effusion. Y-axis indicates the percentage of the CD45pos/CD326pos cells in the whole CD326pos cell sample. X-axis indicates the days after the initial paracentesis. Patients labeled as alive indicate those who were alive at the last data collection; all other patients had died. Among the 4 surviving patients, 1 is being treated with PD-1, and 3 others with TKI. Ins indicates insertion and del indicates deletion.
Discussion
In this study, we showed that varying percentages of CD45pos/CD326pos cells existed in all of the 21 solid tissue samples and 18 pleural effusion samples obtained from 39 patients with NSCLC. These doubly-positive cells are highly suspected to be EMT. They appeared to be malignant microscopically and had cancer-causing genetic mutations. Based on these results, we propose that CD45/CD326 doubly-positive cancer cells exist in both solid lung tumor tissues and pleural effusion in patients with NSCLC. Our present report is the first one demonstrating the evidence of such a cell type.
The phenomenon of EMT in cancer tissues is recognized as an important element in metastasis (15–19). However, EMT studied have searched at the tissue level or in cultured cells and/or tissues (20, 21). One previous report mentioned cells with detectable dual epithelial and mesenchymal characteristics in ascites from patients with gastric cancer (8). However, the existence and significance of CD45/CD326 DPC in primary tumors requires more scrutiny. This is the first report demonstrating the integrity of such cells. We used the EGFR mutation as the marker for cancer cell evaluation and monitored its presence in each analyzed sample and found that the CD45pos/CD326pos cells as well as the CD45neg/CD326pos (epithelial-only singly positive) cells have the same mutation as does each corresponding primary tumor. These indicate that the doubly-positive cells originated from cancer cells. Utilizing these doubly positive cells, we uncovered the intratumoral heterogeneity in the tissues from the patients in this study.
Elucidation of mechanisms of metastasis in solid tumors is a highly important area in cancer research. Epithelial malignant cells need to move via blood vessels or lymphatic ducts in spite of massive immune cell attacks; this is one of the most difficult of current topics (22, 23). EMTs have been suggested as enabling metastasis (24, 25). If some proportion of the epithelial cancer cells can transiently deform to mimic mesenchymal cells, they may be more able to move to other organs, to escape from attacks by immune systems, and to finally accomplish the process of metastasis. The existence of circulating tumor cells may be one example of an EMT process that can lead to worse outcomes (26–28), although precise mechanisms are yet to be determined. Such circulating tumor cells may include the CD45/CD326 DPC we observed in this study. Surprisingly, we observed that the percentages of these CD45/CD326 DPC are higher in malignant effusions than in the solid tumors (Fig. 1B and C). Furthermore, as shown in Figure 4, 2 cases with pleural effusion have over 90% CD45/CD326 DPC, and the total number of cancer cells was lower in these cases. It is possible that CD45/CD326 DPC are more able to infiltrate in pleural effusion than other cancer cells, although the precise mechanisms are not yet understood. It is even possible that such CD45/CD326 DPC may also be able to make the mesenchymal-to-epithelial transition (MET) in addition to the EMT, which is certainly suggested as a requirement in the process of metastasis (29).
In this study, we analyzed primary solid tumor tissue and pleural effusion, but metastatic tumor tissues and serum from these patients were not available. Further investigations using cells from metastatic lesions and/or patients' sera are necessary to verify our hypothesis of the contribution of CD45/CD326 DPC to metastasis. However, detection of CD45/CD326 DPC by liquid biopsy is not technically easy at present (30). Additional investigations using noncancerous samples such as pleural effusion from patients with nonmalignant diseases as controls will add valuable information to future investigations.
In conclusion, we have reported strong indications of the existence of CD45/CD326 DPC in a wide variety of concentrations in both the primary tumor and malignant pleural effusion from patients with NSCLC and shown that a high percentage of such CD45/CD326 DPC is an important prognostic factor. CD45/CD326 DPC will be promising targets for future investigations, particularly in research examining cancer metastasis.
Disclosure of Potential Conflicts of Interest
S. Kobayashi reports receiving speakers bureau honoraria from Bristol-Myers Squibb, Taiho Pharmaceutical Co., Ltd., and AstraZeneca. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: K. Ishizawa, T. Akaishi, A. Horii
Development of methodology: K. Ishizawa, M. Yamanaka, N.-A. Pham
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): K. Ishizawa, Y. Saiki, E. Miyauchi, S. Fukushige, T. Akaishi, A. Asao, R. Saito, Y. Tojo, A. Sakurada, Y. Okada, M.-S. Tsao
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): K. Ishizawa, M. Yamanaka, S. Fukushige, T. Akaishi, T. Mimori, N. Ishii, S. Kobayashi, M. Nagasaki
Writing, review, and/or revision of the manuscript: K. Ishizawa, M. Yamanaka, E. Miyauchi, S. Fukushige, T. Akaishi, T. Ishii, M.-S. Tsao, A. Horii
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): K. Ishizawa, M. Yamanaka, E. Miyauchi, A. Asao, M. Li, Y. Okada S. Kobayashi, M. Ichinose, A. Horii
Study supervision: K. Ishizawa, T. Mimori, R. Yamashita, M. Abe, T. Ishii, N. Ishii, M. Nagasaki, A. Horii
Acknowledgments
We are grateful to Dr. Barbara Lee Smith Pierce (Adjunct Professor, University of Maryland University College) for her editorial work in the preparation of the manuscript. We gratefully acknowledge Drs. Makoto Kobayashi, Osamu Usami, Taizo Hirano, Naoki Tode, Teruyuki Sato, Tomohiro Ichikawa, Yusaku Sasaki, and Tatsuma Okazaki from Department of Respiratory Medicine, Tohoku University, and Ms. Chiharu Kudo, Drs. Hikari Sato, and Yoshikazu Okubo from Japanese Red Cross Ishinomaki Hospital, for their collection of patients' samples. We are also grateful to Biomedical Research Core, Tohoku University School of Medicine for their technical supports. This work was supported in part by Grants-in-Aid for Scientific Research (JP17K09641, 26460468, and 17K17589) from the Japan Society for the Promotion of Science.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.