**Background:** Lung cancer is the most prevalent cancer in China, with 816,563 new cases and 714,699 deaths in 2020. Emerging evidence links the lung microbiome to tumorigenesis through mechanisms including toxin secretion, immune barrier disruption, and oncogenic metabolite release. Previous studies using 16S rRNA sequencing have identified microbial differences in lung cancer, but metagenomic sequencing—which captures bacteria, viruses, and fungi—has been underutilized. This study aimed to screen microbial biomarkers in BALF (chosen over saliva to avoid oral contamination) using metagenomic sequencing to distinguish lung cancer from benign lung disease.
**Methods:** Between February and September 2021, 60 patients with suspicious lung nodules (0.8–3.0 cm) were enrolled at Tianjin Chest Hospital, including 29 with lung cancer (12 adenocarcinoma, 11 squamous cell carcinoma, 6 small cell lung cancer) and 31 with benign lung diseases (most commonly sarcoidosis). Exclusion criteria included acute lung infections, second primary tumors, COPD, pulmonary fibrosis, bronchiectasis, and antibiotic use within 1 month. BALF was collected via standardized bronchoscopy with a negative control (NC) wash. Genomic DNA was extracted using QIAamp DNA Microkit, libraries were constructed with QIAseq Ultralow Input Library Kit, and sequencing was performed on a Nova6000 PE150 platform. Raw reads were quality-filtered (Q<30, <35 bp removed), human reads were removed with Bowtie2, and taxonomic profiling used Kraken2 and Bracken. Contaminants were removed if RPM >50 in NC or ≥1/3 of total NC RPM. Alpha diversity (Chao1, Shannon, Gini–Simpson) and beta diversity (PCA, NMDS, MRPP) were assessed. LEfSe (LDA>2) and Metastat identified differentially abundant taxa. Metabolic pathway analysis used HUMAnN3 with MetaCyc and KEGG (GSVA, adjusted p<0.05, |log2FC|>0.2). Random forest classifiers with leave-one-out cross-validation were built using 40 differentially abundant genera and three tumor markers (CEA, NSE, CYF21-1).
**Key Results:** No significant differences in age, gender, BMI, or smoking status existed between groups. Tumor marker levels (CEA, NSE, CYF21-1) differed significantly. Alpha diversity was generally lower in malignant samples, with Chao1 index significantly different at genus level (p value not explicitly stated in text). MRPP confirmed intergroup differences exceeded intragroup differences. PCA showed malignant samples clustered more tightly than benign samples. MetaCyc analysis revealed amino acid-related metabolism (e.g., L-arginine, L-ornithine biosynthesis) was enriched in benign samples, while carbohydrate-related pathways (starch/sucrose metabolism, pentose/glucuronate interconversions, galactose metabolism) were enriched in malignant samples. LEfSe identified differentially abundant genera including Prevotella, Klebsiella, Mycobacterium, Gordonia, and Sphingobium. Metastat top genera included Achromobacter (benign: 0.000012±0.000008, malignant: 0.003713±0.002897, p=0.000999), Chryseobacterium, Herbaspirillum, Pedobacter, Thermomonas, Undibacterium, Caulobacter, Novosphingobium, Prevotella (benign: 0.189897±0.037106, malignant: 0.064073±0.018263, p=0.001998), and Dechloromonas. Random forest feature importance (MDA>5 or MDG>2) selected five key genera: Klebsiella, Mycobacterium, Pedobacter, Prevotella, and Xanthomonas. The five-genus classifier alone achieved AUC=0.898. Adding NSE yielded the best performance: AUC=0.959, specificity=85.7%, sensitivity=100%, cut-off=0.386. Comparators included all 40 genera (AUC=0.847), five genera+CEA (AUC=0.837), five genera+CYF21-1 (AUC=0.898), and five genera+NSE+CYF21-1 (AUC=0.878).
**Clinical Implications:** The five-genus+NSE model outperformed previously published diagnostic models (Jin et al. AUC=0.882; Cheng et al. AUC=0.845), suggesting strong potential as a clinical diagnostic tool. The differential enrichment of amino acid vs. carbohydrate metabolism pathways points to microbial involvement in lung cancer metabolic reprogramming. Specific genera (e.g., Prevotella linked to ERK/PI3K upregulation; Mycobacterium associated with prior TB and PD-1/PDL1 signaling) may represent therapeutic targets. Limitations include small sample size (n=60), single-center design, lack of validation cohort, no histological subtype stratification, and cross-sectional design preventing causal inference. BALF is less accessible than saliva, which may limit clinical translation. Further multi-center studies with larger cohorts are needed.