Comprehensive analysis of BTNL9 as a prognostic biomarker correlated with immune infiltrations in thyroid cancer

Background Thyroid cancer (THCA) is the most common type of endocrine cancers, and the disease recurrences were usually associated with the risks of metastasis and fatality. Butyrophilin-like protein 9 (BTNL9) is a member of the immunoglobulin families. This study investigated the prognostic role of BTNL9 in THCA. Methods Gene enhancers of BTNL9 were identified by interrogating H3K27ac ChIP-seq and RNA-seq data of papillary thyroid cancer (PTC) and benign thyroid nodule (BTN) tissues. Meanwhile, BTNL9 expression level was verified by qRT-PCR in 30 pairs of primary THCA and adjacent normal tissues. Clinicopathological and RNA sequencing data were obtained from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) to analyze the relations between BTNL9 expression and immune cell infiltration, chemokines/cytokines, immune checkpoint genes, clinical parameters and prognosis values. Besides, survival analysis combining BTNL9 expression and immune cell infiltration scores was conducted. Functional enrichment analysis was performed to investigate the potential biological mechanisms. Cox regression analyses were used to explore independent clinical indicators, and a nomogram model incorporating BTNL9 expression with clinical parameters was established. Results BTNL9 showed significantly stronger H3K27ac modifications in BTN than PTC tissues at the promoter region (chr5: 181,035,673–181,047,436) and gene body (chr5: 181,051,544–181,054,849). The expression levels of BTNL9 were significantly down-regulated in THCA samples compared to normal tissues, and were strongly associated with different tumor stages, immune cell infiltrations, chemokines/cytokines and immune checkpoint genes in THCA. Functional enrichment analyses indicated that BTNL9 was involved in immune-related and cancer-related pathways. The Kaplan–Meier analysis showed lower BTNL9 expression was associated with poorer progression-free interval (PFI). BTNL9 expression and pathologic stages were independent prognostic indicators of PFI in THCA. Conclusions The results implied an important role of BTNL9 in the tumor progression, with the possibility of serving as a novel prognostic biomarker and a potential therapeutic target for THCA.