Skip to main content

Genetic mutation profiles and immune microenvironment analysis of pulmonary enteric adenocarcinoma



Pulmonary enteric adenocarcinoma (PEAC) has distinctive clinical outcomes, radiographic, pathological and molecular characteristics. The prognosis of patients with PEAC was poor. However, molecular profiles and therapeutic biomarkers of PEAC remain elusive.


In the present study, the hospitalized patients with PEAC admitted to Tongji Hospital in Wuhan from January 1, 2014 to November 20, 2020 were retrospectively enrolled and followed until December 10, 2020. Comprehensive genomic profiling of tumor tissue from the PEAC patients were performed and compared with lung adenocarcinoma, colorectal cancer and metastatic colorectal carcinoma. Tumor immune microenvironment analysis were evaluated.


There were 10 patients with PEAC enrolled. 70% of patients were male and the median age of onset was 63 years (interquartile range, 55–72). There were six early-stage patients (Stage IA to IIB) and four stage IV patients. Molecular analysis revealed the most common gene mutations included TP53 (57%, 4/7) and KRAS (57%, 4/7) mutations. There were 40% mutations occurred in genes encoding receptor tyrosine kinases (RTKs). 100% of patients (8/8) were microsatellite stability (MSS). The median level of TMB was 6.0 (interquartile range, 4.5–7.0) mutations/Mb. Three of 10 patients showed low PD-L1 expression (tumor proportion score < 10%) and the others were PD-L1 negative. A small subset of CD8+, CD3+, CD68+ T cells were observed and were mainly distributed in the cancer stroma.


This study demonstrated that PEAC was characterized by low-frequency RTK gene mutation, high KRAS mutation, low PD-L1 expression, low TMB, and low CD8+ T cells infiltration.


Pulmonary enteric adenocarcinoma (PEAC) was classified as a rare variant of invasive adenocarcinoma in 2011 [1] and included in the World Health Organization classification of lung tumors until 2015 [2]. PEAC and pulmonary invasive mucinous adenocarcinoma showed common histomorphological characteristics featured by columnar cells with mucin production. Furthermore, PEAC are similar to metastatic colorectal carcinoma (MCC) in histomorphologic and immunohistochemical features with more than 50% enteric pattern. It is currently difficult to identify PEAC based on information from conventional morphology, immunohistochemistry and radiology. It is necessary to clinically exclude metastases from colorectal carcinoma before the final diagnosis. Investigations on molecular profiles of PEAC indicated that it could be used as potential optimal markers for differential diagnosis [3,4,5].

The treatment of lung adenocarcinoma has been improved in the past two decades with increased understanding in molecular genetics. Targeted therapies against driver mutations have been a major breakthrough in treatment of non-small cell lung cancer (NSCLC). The development of immune checkpoint inhibitor (ICI) represents another promising therapy for advanced NSCLC [6]. A number of biomarkers including programmed cell death ligand 1 (PD-L1) expression, tumor mutational burden (TMB), microsatellite instability (MSI), and CD8+ tumor stroma-infiltrating lymphocyte density have been used to predict response to immunotherapy [7,8,9,10,11].

Patients with PEAC were reported to have a higher rate of lymph node or distant metastases at diagnosis and worse prognosis than other types of lung adenocarcinomas [12, 13]. In view of poor survival after traditional therapy, the understanding of molecular genetics and inhibitory immune checkpoints might provide novel treatment options for PEAC. However, molecular profiles and therapeutic biomarkers of PEAC remain elusive. In the present study, we performed comprehensive genomic profiling and tumor immune microenvironment analysis of PEAC. The molecular profiles were compared with lung adenocarcinoma (LUAD), colorectal cancer (CRC) and MCC.


Patients and samples

This study was a retrospective cohort study of hospitalized patients with PEAC enrolled at Tongji Hospital, Huazhong University of Science and Technology in Wuhan from January 1, 2014 to November 20, 2020. This study was approved by Institutional Review Board of Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology. The diagnosis of PEAC was based on clinical, radiological and pathological information of each patient. Pathologic slides were evaluated by a pathologist to confirm the diagnosis of PEAC, which was defined as a pulmonary adenocarcinoma with an enteric differentiation component exceeding 50%. The possibility of intestinal cancer metastasis in each patient was excluded according to abdominal CT, FDG-PET, gastroscopy and colonoscopy. Patient data including demographic information, laboratory, radiological, and immunohistochemical examinations during hospitalization were extracted from electronic medical records. Patient outcomes were obtained by telephone interview. The last follow-up date was December 10, 2020.

Seven patient tumor tissue samples were sent to Burning Rock Biotech (Guangzhou, China) laboratory which is a College of American Pathologists-accredited and Clinical Laboratory Improvement Amendments-certified clinical laboratory for genetic profiling. The gene panel used in the present study consisting of 520 cancer-related genes (OncoScreen Plus). Two patient tissue samples were detected for common driver gene mutation. The comparison cohorts included LUAD and CRC patients from TCGA and MCC patients from MSKCC, which are available on the website [14, 15]. The LUAD cohort consists of 566 patients with lung adenocarcinoma. The CRC cohort consists of 594 patients with colorectal adenocarcinoma. The MCC cohort consists of 1099 patients with metastatic colorectal cancer. Patient dataset including clinical and genomic information were obtained from the cBioPoratal for cancer genomics (accessed on 14 August 2020).

DNA extraction

The QIAamp DNA FFPE Tissue Kit (Qiagen, UK) were used to extract tumor DNA from FFPE tumor samples according to the manufacturer’s instructions. The Qubit dsDNA HS Assay Kit (Life Technologies, Carlsbad, USA) were used to measure DNA concentration.

Library construction and sequencing

The M220 Focused-ultrasonicator (Covaris, Woburn, MA, USA) was used to shear DNA, followed by end repair, phosphorylation, and adaptor ligation. The Agencourt AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used to select DNA fragments with the range of 200–400 bp. Then, hybridization with capture probe baits, hybrid selection with magnetic beads, and PCR amplification were performed. Target capture was performed with a commercially-available panel of 520 cancer-related genes (OncoScreen Plus). DNA quality and fragment size were assessed by Bioanalyzer 2100 (Agilent, CA, USA). The indexed samples were sequenced on Illumina NextSeq 500 paired-end system (Illumina, Inc., Hayward, CA, USA).

Sequence data analysis

The paired-end reads were mapped to the human genome (hg19) by Burrows-Wheeler aligner v.0.7.10 [16]. Local alignment optimization, variant calling, and annotation were performed with the Genome Analysis Toolkit v.3.2 [17] and VarScan v.2.4.3 [18]. DNA translocation analysis was performed with Factera v.1.4.3 [19]. The variants were filtered with the VarScan filter pipeline, and loci with depths of less than 100 were filtered out. Germline mutations were also filtered out by sequencing matched white blood cells from the samples. Base-calling in tissue samples required at least eight supporting reads for single nucleotide variations and five supporting reads for insertion-deletion variations, respectively. Variants with population frequencies of over 0.1% on the Exome Aggregation Consortium, 1000 Genomes, dbSNP, and ESP6500SI-V2 databases were grouped as single-nucleotide polymorphisms and excluded from further analysis. The remaining variants were annotated with ANNOVAR (2016-02-01 release) [20] and SnpEff v.3.6 [21].

Tumor microenvironment analysis

Seven patients were evaluated for the PD-L1 expression and microenvironment in tumor slides with the PANO 7-plex IHC kit, cat 0004100100 (Panovue, Beijing, China). Expression of surface markers associated with the tumor microenvironment including PD-1, PD-L1, CD3, CD8, CD68, CD56 and CD163 were analyzed through multiplex staining and multispectral imaging [22]. Different primary antibodies were sequentially applied, followed by horseradish peroxidase-conjugated secondary antibody incubation and tyramide signal amplification. The slides were microwave heat-treated after each TSA operation. Nuclei were stained with 4′-6′-diamidino-2-phenylindole (DAPI, Sigma-Aldrich) after all the human antigens had been labelled. To obtain multispectral images, the stained slides were scanned using the Mantra System (PerkinElmer, Waltham, Massachusetts, US), which captures the fluorescent spectra at 20-nm wavelength intervals from 420 to 720 nm with identical exposure time; the scans were combined to build a single stack image. Images of unstained and single-stained sections were used to extract the spectrum of autofluorescence of tissues and each fluorescein, respectively. The extracted images were further used to establish a spectral library required for multispectral unmixing by inForm image analysis software (PerkinElmer, Waltham, Massachusetts, US). Using this spectral library, we obtained reconstructed images of sections with the autofluorescence removed. Three patients were evaluated for the PD-L1 expression of tumor with IHC 22C3 pharmaDx (Daka, Glostrup, Denmark). PD-L1 expression is determined using Tumor Proportion Score (TPS), which is the percentage of viable tumor cells showing partial or complete membrane staining at any intensity.

Microsatellite stability status

FFPE prepared sections were immunostained with automated immunostainer (Dako/Agilent Autostainer Link 48). Primary antibody specific for MLH1 (clone ES05, mouse), PMS2 (clone EP51, rabbit), MSH2 (clone FE11, mouse), and MSH6 (clone EP49, rabbit) (all Readyto-Use, from Dako, Glostrup, Denmark) was applied on the sections according to the manufacturer’s directions. Bound antibody was visualized using the EnVision Kit (Dako, Glostrup, Denmark). Nuclear staining in cancer cells with any intensity was recognized as positive. And nuclear staining in stromal or inflammatory cells was additionally recorded as an internal control. Suspected MMR deficiency was defined as complete loss of at least one of the 4 MMR proteins in all tumor cells.

Statistical analysis

The continuous variables were presented as mean or median. The categorical variables were presented as frequencies. Unpaired Wilcoxon signed-rank test was used to compare continuous variables, while two-sided Fisher’s exact tests were used to compare categorical variables, as appropriate. P < 0.05 was considered statistically significant. All bioinformatics analyses were performed with R (v.3.5.3, the R Foundation for Statistical Computing, Vienna, Austria).


Clinical and radiological characteristics

A total of 10 patients with PEAC were enrolled. Histopathological diagnoses of PEAC from tumor tissues (three biopsy samples from three patients P08-P10 and seven surgical resections from seven patients P01-P07) after hematoxylin and eosin staining were confirmed (Supplemental Fig. 1). 70% of patients (7/10) were male and the age at diagnosis ranged from 43 to 76 years with the median age at diagnosis was 63 years (interquartile range, 55–72). A half of patients (5/10) were current smokers with a median smoking history of 30 pack years. The serum levels of carcinoembryonic antigen were elevated in a half of patients. There were seven patients with tumor located in right lower lobes and three of them showed similar radiological features to pneumonic infiltrate or consolidation. There were six early-stage patients (Stage IA to IIB). Four patients showed lymph node and distant metastases (Stage IV), which affected the pleural, breast, pericardium and bone. There were seven patients underwent surgical treatment including one patient with pleural metastasis underwent lobectomy and pleurotomy. The median disease-free survival was 20.5 months (interquartile range, 16–28.3). Two patients with early-stage (IB and IIB) died 13–14 months after diagnosis and one of them died of pulmonary tuberculosis. The clinical and radiological data are summarized in Table 1.

Table 1 Demographic, clinical and radiological characteristics of the PEAC patients

Immunohistochemical analysis

Immunohistochemical analysis revealed that patients were 100% (10/10) positive for CK7, 60% (6/10) positive for TTF-1, 63% (5/8) positive for Napsin A; 50% (3/6) positive for CK20, 80% (8/10) positive for CDX-2, 70% (7/10) positive for villin, 78% (7/9) positive for P53, 17% (1/6) positive for SATB, whereas patients were 100% (0/8) negative for ALK fusion (Table 2). For Ki-67, there were eight patients had high (≥10%) expression, and two patients had low expression.

Table 2 Immunohistochemical features in the patients with PEAC

Genetic mutation spectrum analysis in the PEAC cohort

Targeted sequencing of 520 cancer related genes in seven patients (P01-P06 and P9) revealed 69 genomic alterations in 54 genes, including 34 missenses (49%), 9 frameshift (13%), 9 copy number amplifications (13%), 6 stop-gained mutations (9%), 2 splice acceptor variant (3%), 2 splice region variants (3%) and so on (Supplemental Table 1). Each patient carried at least one genetic mutation identified by sequencing. The most common gene mutations included TP53 (57%, 4/7), KRAS (57%, 4/7) mutations (Fig. 1). EGFR (19del) mutation was detected in a stage IA patient (P03, Fig. 1), an ERBB2 (20ins) mutation in a stage IV patient (P02), and a BRAF (V600M) mutation in a stage IIB patient (P04). An APC mutation was detected in a stage IIB patient (P05).

Fig. 1
figure 1

Mutation spectrum analysis heatmap in seven pulmonary enteric adenocarcinoma patients (P01-P06 and P09)

Common driver mutations (EGFR, ALK, BRAF, ERBB2, KRAS, MET, RET, and ROS1) analysis were performed in lung adenocarcinoma tissue samples from two PEAC patients (P08 and P10) using next-generation sequencing. Only KRAS mutation was detected in one patient (P08).

In our study of seven PEAC patients (P01-P06 and P09), genetic mutations in 414 of 520 targeted genes were comparable with LUAD, CRC, MCC patients from TCGA and the MSKCC databases. The mutation rates of PPP2R1A, KEAP1, KRAS, MED12, TP53 and GNAS were significantly different between PEAC patients and MCC patients (all P < 0.05) (Fig. 2). APC mutation was rare in LUAD but was common in CRC and MCC. Interestingly, mutations in KRAS were more frequently occurred in PEAC than other three tumor types.

Fig. 2
figure 2

Detection rate statistics in the PEAC cohort and comparison with TCGA and the MSKCC databases. CRC, colorectal cancer; LUAD, lung adenocarcinoma; MCC, metastatic colorectal carcinoma; PEAC, pulmonary enteric adenocarcinoma

Functional enrichment analysis based on the panel demonstrated that RAS family mutations account for half of gene mutations in PEAC. There were 40% mutations occurred in genes encoding receptor tyrosine kinases (RTKs) including EGFR, MET, BRAF, ERBB2, ALK, ROS1 and RET (Fig. 3). In general, frequencies of genetic mutation in PEAC patients were similar to those in LUAD patients but were different from those in CRC and MCC patients.

Fig. 3
figure 3

Functional enrichment analysis in the PEAC cohort and comparison with TCGA and the MSKCC databases. CRC, colorectal cancer; LUAD, lung adenocarcinoma; MCC, metastatic colorectal carcinoma; PEAC, pulmonary enteric adenocarcinoma

Single nucleotide variant analysis and mutation signature analysis

Single nucleotide variant analysis in PEAC showed that the C:G > A:T substitutions were common while the C:G > G:C substitutions were rare (Fig. 4). Comparisons in different tumor types indicated that C > A substitution subtype was more common in PEAC and LUAD, whereas C > T was more frequent in CRC and MCC.

Fig. 4
figure 4

Single nucleotide variant analysis in the PEAC cohort and comparison with TCGA and the MSKCC databases. CRC, colorectal cancer; LUAD, lung adenocarcinoma; MCC, metastatic colorectal carcinoma; PEAC, pulmonary enteric adenocarcinoma

The mutational signature was indicated by the six substitution subtypes (C > A, C > G, C > T, A > T, A > C, and A > G). In general, the profile of mutational signature in PEAC was similar to LUAD (Fig. 5). Signature 1 was presented in all cancer types as the result of an endogenous mutational process initiated by spontaneous deamination of 5-methylcytosine. Signature 4 and signature 29 were associated with tobacco use, however, their profiles in PEAC were distinctly different from LUAD. Signature 6 and signature 15 which were associated with defective DNA mismatch repair and microsatellite instability, were more common in colorectal cancer and rare in PEAC and LUAD.

Fig. 5
figure 5

Mutation signature analysis in the PEAC cohort and comparison with TCGA and the MSKCC databases. CRC, colorectal cancer; LUAD, lung adenocarcinoma; MCC, metastatic colorectal carcinoma; PEAC, pulmonary enteric adenocarcinoma

Tumor mutation burden and microsatellite stability status

As indicated in Table 3, the median level of TMB in seven PEAC patients was 5.98 (interquartile range, 4.49–6.98) mutations/Mb, ranging from 2.99 to 9.97 mutations/Mb. TMB was comparable between PEAC and lung adenocarcinoma across the TCGA cohort (median: 5.78; range: 0.5–48) [23].

Table 3 Tumor mutational burden, microsatellite stability status, and multiplex immunohistochemical predictive biomarkers

MSI is the molecular fingerprint of a deficient mismatch repair system. Immunohistochemistry analysis for DNA mismatch repair proteins (MLH1, PMS2, MSH2, MSH6) revealed that eight patients (P01-P07 and P09) were microsatellite stability (MSS) with positive for all microsatellite markers (Table 3).

Immunophenotypes and microenvironment analysis of PEAC

Tumor-infiltrating lymphocytes have been associated with response to immunotherapy. Multiplex immunohistochemistry assay was performed in seven patients (P01-P06 and P09) to simultaneously measure the expression of the multiple surface markers associated with the tumor microenvironment, such as PD-1, PD-L1, CD3, CD8, CD68, CD56 and CD163 (Table 3). 57% of the cases (4/7) showed positivity for PD-L1 at low expression level (1–49%). PD-1/PD-L1 expression was higher in cancer stroma compared with cancer nests in most patients. The CD8+, CD3+, CD68+ T cells were mainly distributed in the cancer stroma, whereas the CD56+ T cells were mainly located in cancer nests (Fig. 6).

Fig. 6
figure 6

Multiplex immunofluorescence staining for PD-1, PD-L1, CD3, CD8, CD68, CD56 and CD163 in the tumor tissue of the patients with PEAC

PD-L1 expressions were evaluated in other three patients (P07, P08 and P10) through immunohistochemical assay using anti-PD-L1 antibody clone 22C3. All three patients showed negative PD-L1 expression (TPS = 0%).


PEAC has distinctive clinical outcomes, radiographic, pathological and molecular characteristics. The prognosis of patients with PEAC was poor [12, 13]. 40% of patients in this study showed lymph node and distant metastases and were diagnosed with stage IV. Two patients with early-stage died 13–14 months after diagnosis.

The radiological study indicated that lesions in PEAC were larger and more solid compared to primary pulmonary invasive adenocarcinoma [24]. In this study, three patients were diagnosed as pneumonic-type of lung cancer. The tumors in seven patients were larger than 5 cm. Nine of the ten patients had unilateral lung carcinomas and seven of them were located in the right lower lobes.

The immunohistochemical markers of lung adenocarcinoma (CK7, TTF-1, Napsin A) and enteric differentiation (CDX-2, CK20, MUC2, villin) were both expressed in PEAC. The combination of CK7+/CDX-2+ was reported to have a high sensitivity (71.3%) and specificity (82%) in differential diagnosis of PEAC from colorectal adenocarcinoma, [25] while combining cadherin-17 (negative) and SATB homeobox 2 (negative) also showed high sensitivity (77.0%) and specificity (100%) [26]. Villin, β-catenin and SATB2 were served as useful immunohistochemical markers for differential diagnosis between PEAC and MCC [27, 28]. Our findings were in line with previous studies and showed positive rates of immunohistochemical markers in PEAC as following: CK7, 100% (10/10); TTF-1, 60% (6/10); Napsin A 63% (5/8); CK20, 50% (3/6); CDX-2, 80% (8/10); villin 70% (7/10).

The profile of driver mutations in PEAC was reported to be similar to NSCLC, but different from MCC [4]. Compared with pulmonary invasive adenocarcinoma, the PEAC indicated a higher incidence (14.3–63%) of KRAS mutation but fewer NSCLC RTK mutations [24, 25, 28]. [5, 29, 30] However, a study revealed that 92.31% (12/13) PEAC patients harbored mutations in well-established RTK genes (EGFR, ALK, ERBB2, BRAF) for NSCLC [4]. Molecular analysis in our study revealed KRAS and TP53 as the most frequently mutated genes (57%, 4/7). Only three cases harbored abnormalities affecting EGFR, ERBB2, and BRAF genes.

As for NSCLC patients without known driver mutation, immunotherapy represents a promising therapeutic strategy for patients with advanced or metastatic disease. Most of ICI targets the PD-1/PD-L1 axis to restore anti-tumor immunity. A number of clinical trials and real-word studies showed that ICI was an important treatment modality for patients with lung adenocarcinoma and significantly improved the prognosis of selected patients [31]. PD-1 targeting agent nivolumab showed treatment effective on NSCLC in the real-world study. Subgroup analyses showed patients with higher PD-L1 expression experienced a greater clinical benefit from nivolumab [32]. PD-L1 expression in tumor cells was identified as a valuable predictor of the efficacy of anti-PD-1/PD-L1 monotherapy in certain NSCLC patients [33]. A previous study in primary lung adenocarcinomas with intestinal differentiation showed that PD-L1 expression was limited to PEAC with a positive rate of 42.9% (3/7), and varied TPS, TPS 20, 50, and 80% respectively [12]. However, in our study four patients showed low PD-L1 expression and the remaining six patients were PD-L1 negative.

Furthermore, MSI, TMB, and CD8+ T-cell tumor-infiltrating lymphocytes have been identified as promising biomarkers to evaluate patients’ survival and response to PD-1/PD-L1 blockade [7,8,9,10,11]. The patients with high MSI were associated with benefit from immunotherapy [34]. Higher TMB also predicted favorable outcome to immunotherapy in NSCLC patients [9]. High MSI and high TMB both reflect instability in tumor cells, and usually occurred simultaneously in same. However, a recent study showed that MSS and TMB-high were more commonly occurred and might be benefit from immunotherapy [8]. The CD8+ T-cell tumor-infiltrating lymphocytes were associated with improved prognosis [10]. A recent study showed that PD-L1 expression on tumor cells in combination with CD8+ T-cell density were predictive biomarkers in patients with inoperable locally advanced NSCLC treated with concurrent chemoradiotherapy. The longest and shortest OS were observed in patients with PD-L1 negative/CD8+ low and PD-L1 positive/CD8+ low tumors respectively [35]. A previous study including 18 Chinese PEAC patients showed that compared with lung adenocarcinomas, PEAC had higher nonsynonymous TMB and MSI [25]. A study of seven Germany PEAC patients revealed a median TMB of 16.8 mutations/Mb which is much higher than that in our study (5.98 mutations/Mb), and none of them were MSI [12]. The mutation load in 13 Chinese PEAC patients in Zhang’s study was similar to our study [4]. The differences in TMB between the previous studies and our studies may be due to the limited sample size, different races, and genetic heterogeneity of PEAC patients. Above all, a majority of tumor in our study were MSS, low TMB, and low CD8+ T cells infiltration.

The analysis of potential genomic biomarkers of immunotherapy in 1000 Chinese patients with cancer showed that lung cancer patients with EGFR mutations had significantly lower TMB than those with wild-type EGFR [36]. The TMB of the patient (P03) with EGFR mutation was 3.99 mutations/Mb in our study. Lung adenocarcinoma with TP53 mutation/STK11-EGFR wild-type was reported to have higher CD8+ T-cell density and PD-L1 expression than other tumor subtypes, and the group of patients had a prolonged progression-free survival from immunotherapy [37]. There were three patients (P01, P05 and P06) harboring TP53 mutation/STK11-EGFR wild-type in our study, nevertheless, no differences were observed in CD8+ T-cell infiltration and PD-L1 expression in the three patients.

The main limitation of this study is the retrospective, observational nature of the study in a single-center with small sample size. It was difficult to draw a comprehensive conclusion regarding molecular genetics profile and predictive biomarkers for immunotherapy in PEAC. Large-scale prospective cohort studies are needed to validate our findings. Additional limit7ations include lacks of treatment response evaluation and analysis of prognostic factors due to the limited cases.


In summary, this study demonstrated that PEAC was characterized by low-frequency RTK gene mutation, high KRAS mutation, low PD-L1 expression, low TMB, and low CD8+ T cells infiltration, providing important information for the development of therapeutic strategies for patients with PEAC.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files.



Pulmonary enteric adenocarcinoma


Non-small cell lung cancer


Metastatic colorectal carcinoma


Programmed cell death ligand 1


Tumor mutational burden


Microsatellite instability


Lung adenocarcinoma


Colorectal cancer


Encoding receptor tyrosine kinases


Microsatellite stability


Tumor Proportion Score


Immune checkpoint inhibitors


  1. Travis WD, Brambilla E, Noguchi M, Nicholson AG, Geisinger K, Yatabe Y, et al. International Association for the Study of Lung Cancer/American Thoracic Society/European Respiratory Society: international multidisciplinary classification of lung adenocarcinoma: executive summary. Proc Am Thorac Soc. 2011;8(5):381–5.

    Article  PubMed  Google Scholar 

  2. Travis WD, Brambilla E, Nicholson AG, Yatabe Y, Austin JHM, Beasley MB, et al. The 2015 World Health Organization classification of lung tumors: impact of genetic, clinical and radiologic advances since the 2004 classification. J Thorac Oncol. 2015;10(9):1243–60.

    Article  PubMed  Google Scholar 

  3. Jurmeister P, Scholer A, Arnold A, Klauschen F, Lenze D, Hummel M, et al. DNA methylation profiling reliably distinguishes pulmonary enteric adenocarcinoma from metastatic colorectal cancer. Mod Pathol. 2019;32(6):855–65.

    Article  CAS  PubMed  Google Scholar 

  4. Zhang J, Xiang C, Han Y, Teng H, Li X, Shao J, et al. Differential diagnosis of pulmonary enteric adenocarcinoma and metastatic colorectal carcinoma with the assistance of next-generation sequencing and immunohistochemistry. J Cancer Res Clin Oncol. 2019;145(1):269–79.

    Article  CAS  PubMed  Google Scholar 

  5. Garajova I, Funel N, Fiorentino M, Agostini V, Ferracin M, Negrini M, et al. MicroRNA profiling of primary pulmonary enteric adenocarcinoma in members from the same family reveals some similarities to pancreatic adenocarcinoma-a step towards personalized therapy. Clin Epigenetics. 2015;7(1):129.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Brody R, Zhang Y, Ballas M, Siddiqui MK, Gupta P, Barker C, et al. PD-L1 expression in advanced NSCLC: insights into risk stratification and treatment selection from a systematic literature review. Lung Cancer. 2017;112:200–15.

    Article  PubMed  Google Scholar 

  7. Dudley JC, Lin MT, Le DT, Eshleman JR. Microsatellite instability as a biomarker for PD-1 blockade. Clin Cancer Res. 2016;22(4):813–20.

    Article  CAS  PubMed  Google Scholar 

  8. Goodman AM, Sokol ES, Frampton GM, Lippman SM, Kurzrock R. Microsatellite-stable tumors with high mutational burden benefit from immunotherapy. Cancer Immunol Res. 2019;7(10):1570–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Chae YK, Davis AA, Raparia K, Agte S, Pan A, Mohindra N, et al. Association of Tumor Mutational Burden with DNA repair mutations and response to anti-PD-1/PD-L1 therapy in non-small-cell lung Cancer. Clin Lung Cancer. 2019;20(2):88–96.e6.

    Article  CAS  PubMed  Google Scholar 

  10. Yu Y, Zeng D, Ou Q, Liu S, Li A, Chen Y, et al. Association of Survival and Immune-Related Biomarkers with Immunotherapy in patients with non-small cell lung Cancer: a Meta-analysis and individual patient-level analysis. JAMA Netw Open. 2019;2(7):e196879.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Goodman AM, Kato S, Bazhenova L, Patel SP, Frampton GM, Miller V, et al. Tumor mutational burden as an independent predictor of response to immunotherapy in diverse cancers. Mol Cancer Ther. 2017;16(11):2598–608.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Jurmeister P, Vollbrecht C, Behnke A, Frost N, Arnold A, Treue D, et al. Next generation sequencing of lung adenocarcinoma subtypes with intestinal differentiation reveals distinct molecular signatures associated with histomorphology and therapeutic options. Lung Cancer. 2019;138:43–51.

    Article  PubMed  Google Scholar 

  13. Feng C, Feng M, Gao Y, Zhao X, Peng C, Yang X, et al. Clinicopathologic significance of intestinal-type molecules’ expression and different egfr gene status in pulmonary adenocarcinoma. Appl Immunohistochem Mol Morphol. 2019;27(5):364–72.

    Article  CAS  PubMed  Google Scholar 

  14. National Library of Medicine, National Center for Biotechnology Information. Available:

  15. Yaeger R, Chatila WK, Lipsyc MD, Hechtman JF, Cercek A, Sanchez-Vega F, et al. Clinical sequencing defines the genomic landscape of metastatic colorectal cancer. Cancer Cell. 2018;33(1):125–36

  16. Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Koboldt DC, Zhang Q, Larson DE, Shen D, McLellan MD, Lin L, et al. VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012;22(3):568–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Newman AM, Bratman SV, Stehr H, Lee LJ, Liu CL, Diehn M, et al. FACTERA: a practical method for the discovery of genomic rearrangements at breakpoint resolution. Bioinformatics. 2014;30(23):3390–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Cingolani P, Platts A, Wang le L, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012;6(2):80–92.

    Article  CAS  Google Scholar 

  22. Stack EC, Wang C, Roman KA, Hoyt CC. Multiplexed immunohistochemistry, imaging, and quantitation: a review, with an assessment of Tyramide signal amplification, multispectral imaging and multiplex analysis. Methods. 2014;70(1):46–58.

    Article  CAS  PubMed  Google Scholar 

  23. Cancer Genome Atlas Research N. Comprehensive molecular profiling of lung adenocarcinoma. Nature. 2014;511(7511):543–50.

    Article  CAS  Google Scholar 

  24. Zhao L, Huang S, Liu J, Zhao J, Li Q, Wang HQ. Clinicopathological, radiographic, and oncogenic features of primary pulmonary enteric adenocarcinoma in comparison with invasive adenocarcinoma in resection specimens. Medicine (Baltimore). 2017;96(39):e8153.

    Article  Google Scholar 

  25. Chen M, Liu P, Yan F, Xu S, Jiang Q, Pan J, et al. Distinctive features of immunostaining and mutational load in primary pulmonary enteric adenocarcinoma: implications for differential diagnosis and immunotherapy. J Transl Med. 2018;16(1):81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Bian T, Zhao J, Feng J, Zhang Q, Qian L, Liu J, et al. Combination of cadherin-17 and SATB homeobox 2 serves as potential optimal makers for the differential diagnosis of pulmonary enteric adenocarcinoma and metastatic colorectal adenocarcinoma. Oncotarget. 2017;8(38):63442–52.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Gu L, Wang XZ, Wen W, Lin J, Chen XF, Lai GX, et al. Clinical analysis of 23 patients pathologically diagnosed with primary and secondary pulmonary enteric adenocarcinoma. Chin Med J. 2019;132(11):1368–9.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Matsushima J, Yazawa T, Suzuki M, Takahashi Y, Ota S, Nakajima T, et al. Clinicopathological, immunohistochemical, and mutational analyses of pulmonary enteric adenocarcinoma: usefulness of SATB2 and beta-catenin immunostaining for differentiation from metastatic colorectal carcinoma. Hum Pathol. 2017;64:179–85.

    Article  CAS  PubMed  Google Scholar 

  29. Nottegar A, Tabbo F, Luchini C, Guerrera F, Gaudiano M, Bria E, et al. Pulmonary adenocarcinoma with enteric differentiation: dissecting oncogenic genes alterations with DNA sequencing and FISH analysis. Exp Mol Pathol. 2017;102(2):276–9.

    Article  CAS  PubMed  Google Scholar 

  30. Nottegar A, Tabbo F, Luchini C, Brunelli M, Bria E, Veronese N, et al. Pulmonary adenocarcinoma with enteric differentiation: immunohistochemistry and molecular morphology. Appl Immunohistochem Mol Morphol. 2018;26(6):383–7.

    Article  CAS  PubMed  Google Scholar 

  31. Niu Y, Lin A, Luo P, Zhu W, Wei T, Tang R, et al. Prognosis of lung adenocarcinoma patients with NTRK3 mutations to immune checkpoint inhibitors. Front Pharmacol. 2020;11:1213.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Morita R, Okishio K, Shimizu J, Saito H, Sakai H, Kim YH, et al. Real-world effectiveness and safety of nivolumab in patients with non-small cell lung cancer: a multicenter retrospective observational study in Japan. Lung Cancer. 2020;140:8–18.

    Article  PubMed  Google Scholar 

  33. Xu Y, Wan B, Chen X, Zhan P, Zhao Y, Zhang T, et al. The association of PD-L1 expression with the efficacy of anti-PD-1/PD-L1 immunotherapy and survival of non-small cell lung cancer patients: a meta-analysis of randomized controlled trials. Transl Lung Cancer Res. 2019;8(4):413–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Ren D, Hua Y, Yu B, Ye X, He Z, Li C, et al. Predictive biomarkers and mechanisms underlying resistance to PD1/PD-L1 blockade cancer immunotherapy. Mol Cancer. 2020;19(1):19.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Gennen K, Kasmann L, Taugner J, Eze C, Karin M, Roengvoraphoj O, et al. Prognostic value of PD-L1 expression on tumor cells combined with CD8+ TIL density in patients with locally advanced non-small cell lung cancer treated with concurrent chemoradiotherapy. Radiat Oncol. 2020;15(1):5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Zang YS, Dai C, Xu X, Cai X, Wang G, Wei J, et al. Comprehensive analysis of potential immunotherapy genomic biomarkers in 1000 Chinese patients with cancer. Cancer Med. 2019;8(10):4699–708.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Biton J, Mansuet-Lupo A, Pecuchet N, Alifano M, Ouakrim H, Arrondeau J, et al. TP53, STK11, and EGFR mutations predict tumor immune profile and the response to anti-PD-1 in lung adenocarcinoma. Clin Cancer Res. 2018;24(22):5710–23.

    Article  CAS  PubMed  Google Scholar 

Download references


We sincerely thank Yiran Cai, Zhiqiu Chen, and Yuxi Zhou who are working at Burning Rock Biotech company for their contribution to this study.


Not applicable.

Author information

Authors and Affiliations



Min Xie, Dong Kuang and Xiaochen Li performed study concept and design; Xiaochen Li and Min Xie performed development of methodology, acquisition, analysis, interpretation of data, statistical analysis, writing, review and revision of the paper; Dong Chen, Yong Li and Xiansheng Liu provided technical and material support. All authors read and approved the final paper.

Corresponding authors

Correspondence to Dong Kuang or Xiaochen Li.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the institutional review board of Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology.

Consent for publication

Not applicable.

Competing interests

The authors have no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplemental Fig. 1.

Hematoxylin and eosin staining of tumor tissues from all PEAC patients.

Additional file 2: Supplemental Table 1.

Genetic mutational profile of PEAC using a 520-gene OncoScreen Plus panel.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xie, M., Chen, D., Li, Y. et al. Genetic mutation profiles and immune microenvironment analysis of pulmonary enteric adenocarcinoma. Diagn Pathol 17, 30 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: