As the largest gene family functioning in protein transport among human solute carriers, the SLC25 family (mitochondrial carrier family) can participate in development of cancer. However, a comprehensive exploration for the exactly roles of SLC family remains lacking. In the present study, a total of 15 functional SLC25 family genes were retrieved from all current publications. And multidimensional analyses were systematically performed based on the transcriptome and genome data of SLC25 family from a variety of online databases for their expression, immune cell infiltration, and cancer prognosis. Validation by qPCR and immunohistochemistry were further conducted for the expression of partial SLC25 family members in some tumor tissue. We found that the SLC25 family had strong correlation with immune cells, such as macrophages M2, CD8+ T cell, CD4+ T cell memory activated, and memory resting. Among them, SLC25A6 was most correlated with Macrophage M1 in uveal melanoma ( r = − 0.68 , P = 1.9 e − 0.5 ). Expression of mRNA level showed that SLC25A4 was downregulated in stomach adenocarcinoma and colon adenocarcinoma. SLC25A7 was highly expressed in stomach adenocarcinoma and colon adenocarcinoma. SLC25A23 was decreased in colon adenocarcinoma. qPCR and immunohistochemistry validation results were consistent with our bioinformatics prediction. SLC25A8 was associated with the prognosis of cancer. All these findings suggested that the SLC25 family might affects the immune microenvironment of the cancer and then had the potential to be predictive biomarkers for early diagnosis and prognosis as well as novel targets for individualized treatment of cancer.

BackgroundRecently, the incidence of cholangiocarcinoma (CCA) has gradually increased. As CCA has a poor prognosis, the ideal survival rate is scarce for patients. The abnormal expressed tsRNAs may regulate the progression of a variety of tumors, and tsRNAs is expected to become a new diagnostic biomarker of cancer.However, the expression of tsRNAs is obscure and should be elucidated in CCA. MethodsHigh-throughput RNA sequencing technology (RNA-seq) was utilized to determine the overall expression pro les of tsRNAs in 3 pairs CCA and adjacent normal tissues and to screen the tsRNAs that were differentially expressed. The target genes of dysregulated tsRNAs were predicted and the biological effects and potential signaling pathways of these target genes were explored by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to validate 11 differentially expressed tRFs with 12 pairs CCA and adjacent normal tissues. ResultsHigh-throughput RNA-seq totally demonstrated 535 dysregulated tsRNAs, of which 241 tsRNAs were upregulated and 294 tsRNAs were downregulated in CCA compared with adjacent normal tissues (|log2 (fold change) |≥1 and P value < 0.05). GO and KEGG enrichment analyses indicated that the target genes of dysregulated tRFs (tRF-34-JJ6RRNLIK898HR, tRF-38-0668K87SERM492V and tRF-39-0668K87SERM492E2) were mainly enriched in the Notch signaling pathway, Hippo signaling pathway, cAMP signaling pathway and in growth hormone synthesis, secretion and action, etc. qRT-PCR result showed that tRF-34-JJ6RRNLIK898HR/tRF-38-0668K87SERM492V/tRF-39-0668K87SERM492E2 was down-regulated (P = 0.021) and tRF-20-LE2WMK81 was up-regulated in CCA (P = 0.033). ConclusionDifferentially expressed tRFs in CCA are enriched in many pathways associated with neoplasms, which may impact the tumor progression and have potential to be diagnostic biomarkers and therapeutic targets of CCA.