Molecular characterization of early-stage multi-primary lung adenocarcinoma by transcriptome sequencing-a retrospective study.

[BACKGROUND] To investigate the molecular genetic features of multiple primary lung cancer (MPLC) to provide a basis and new methods for its identification, diagnosis, and treatment.

[METHODS] Transcriptome sequencing (RNA-seq) was performed on 16 tissue samples from eight patients with synchronous multiple primary lung adenocarcinoma (sMP-LUAD) and eight tissue samples from eight patients with single primary lung adenocarcinoma (SP-LUAD). Differentially expressed genes selected by bioinformatic methods were validated in 24 sets of sMP-LUAD and SP-LUAD samples using quantitative reverse transcription polymerase chain reaction (qRT-PCR). Based on The Cancer Genome Atlas (TCGA) database, the differentially expressed genes responsible for the biological behavior of lung adenocarcinoma and their clinical significance were analyzed.

[RESULTS] Overall, 194 differentially expressed genes were identified ( < 0.05), including 22 up-regulated and 172 down-regulated genes. Two up-regulated ( and ) and three down-regulated (, , and ) genes were selected for validation by qRT-PCR analysis. The qRT-PCR results showed that the expression of mRNA in the sMP-LUAD group was significantly higher ( < 0.05) than that in the SP-LUAD group; mRNA , , , and expression in the sMP-LUAD group was not statistically different from that in the SP-LUAD group ( > 0.05). Gene ontology (GO) enrichment analysis showed that mRNA was mainly enriched in the entries of positive regulation of wound healing and oxidation-reduction processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that DUOX1 can promote reactive oxygen species (ROS) production and be related to thyroid hormone production. Furthermore, based on the TCGA database, we analyzed the biological behavior and clinical significance of DUOX1 in lung adenocarcinoma using bioinformatics technology. mRNA expression was decreased in all stages and pathological subtypes of lung adenocarcinoma ( < 0.05). Immune infiltration analysis showed that DUOX1 with mast cells and eosinophils was positively correlated ( < 0.05), and Th2 cells were negatively correlated ( < 0.05). Logistic regression analysis showed that the expression of mRNA was significantly correlated with the patient's age, lymph node metastasis, and pathologic stage ( < 0.05). Kaplan-Meier survival plots showed that low DUOX1 expression was not significantly correlated with OS, DSS, and PFI ( > 0.05). Univariate and multivariate COX regression analysis revealed that mRNA could not be used as an independent factor for predicting prognosis ( > 0.05). Therefore, we developed a predictive nomogram model combining clinicopathological variables and mRNA to predict the survival of patients with lung adenocarcinoma.