Integrative analysis reveals methylenetetrahydrofolate dehydrogenase 1-like as an independent shared diagnostic and prognostic biomarker in five different human cancers

Background: Defects in methylenetetrahydrofolate dehydrogenase 1-like (MTHFD1L) expression have earlier been examined in only a few human cancers. Objectives: Multi-omics profiling of MTHFD1L as a shared biomarker in distinct subtypes of human cancers. Methods: In the current study, for the multi-omics analysis of MTHFD1L in 24 major subtypes of human cancers, a comprehensive in silico approach was adopted to mine different open access online databases including UALCAN, Kaplan–Meier (KM) plotter, LOGpc, GEPIA, Human Protein Atlas (HPA), Gene Expression across Normal and Tumor tissue (GENT2), MEXPRESS, cBioportal, STRING, DAVID, TIMER, and Comparative Toxicogenomics Database (CTD). Results: We noticed that the expression of MTHFD1L was significantly higher in all the analyzed 24 subtypes of human cancers as compared with the normal controls. Moreover, MTHDF1L overexpression was also found to be significantly associated with the reduced overall survival (OS) duration of Bladder urothelial cancer (BLCA), Head and neck cancer (HNSC), Kidney renal papillary cell carcinoma (KIRP), Lung adenocarcinoma (LUAD), and Uterine corpus endometrial carcinoma (UCEC). This implies that MTHFD1L plays a significant role in the development and progression of these cancers. We further noticed that MTHFD1L was also overexpressed in BLCA, HNSC, KIRP, LUAD, and UCEC patients of different clinicopathological features. Pathways enrichment analysis revealed the involvement of MTHFD1L-associated genes in five diverse pathways. We also explored few interesting correlations between MTHFD1L expression and its promoter methylation, genetic alterations, CNVs, and between CD8+ T immune cells level. Conclusion: In conclusion, our results elucidated that MTHFD1L can serve as a shared diagnostic and prognostic biomarker in BLCA, HNSC, KIRP, LUAD, and UCEC patients of different clinicopathological features.     

Prognostic and immunological role of Fam20C in pan-cancer

Background: The family with sequence similarity 20-member C (Fam20C) kinase plays important roles in physiopathological process and is responsible for majority of the secreted phosphoproteome, including substrates associated with tumor cell migration. However, it remains unclear whether Fam20C plays a role in cancers. Here, we aimed to analyze the expression and prognostic value of Fam20C in pan-cancer and to gain insights into the association between Fam20C and immune infiltration.Methods: We analyzed Fam20C expression patterns and the associations between Fam20C expression levels and prognosis in pan-cancer via the ONCOMINE, TIMER (Tumor Immune Estimation Resource), PrognoScan, GEPIA (Gene Expression Profiling Interactive Analysis), and Kaplan–Meier Plotter databases. After that, GEPIA and TIMER databases were applied to investigate the relations between Fam20C expression and immune infiltration across different cancer types, especially BLCA (bladder urothelial carcinoma), LGG (brain lower grade glioma), and STAD (stomach adenocarcinoma).Results: Compared with adjacent normal tissues, Fam20C was widely expressed across many cancers. In general, Fam20C showed a detrimental role in pan-cancer, it was positively associated with poor survival of BLCA, LGG, and STAD patients. Specifically, based on TCGA (The Cancer Genome Atlas) database, a high expression level of Fam20C was associated with worse prognostic value in stages T2–T4 and stages N0–N2 in the cohort of STAD patients. Moreover, Fam20C expression had positive associations with immune infiltration, including CD4⁺ T cells, macrophages, neutrophils, and dendritic cells, and other diverse immune cells in BLCA, LGG, and STAD.Conclusion: Fam20C may serve as a promising prognostic biomarker in pan-cancer and has positive associations with immune infiltrates.     

In silico assessment of EpCAM transcriptional expression and determination of the prognostic biomarker for human lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC)

Epithelial cell adhesion molecule (EpCAM) is a transmembrane glycoprotein which is involved in cell signaling, proliferation, maturation, and movement, all of which are crucial for the proper development of cells and tissues. Cleavage of the EpCAM protein leads to the up-regulation of c-myc, e-fabp, and cyclins A and E which promote tumorigenesis. EpCAM can act as potential diagnostic and prognostic biomarker for different types of cancers as it is also found to be expressed in epithelia and epithelial-derived neoplasms. Hence, we aimed to analyze the EpCAM gene expression and any associated feedback in the patients of two major types of lung cancer (LC) i.e., lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC), based on the publicly available online databases. In this study, server-based gene expression analysis represents the up-regulation of EpCAM in both LUAD and LUSC subtypes as compared to the corresponding normal tissues. Besides, the histological sections revealed the over-expression of EpCAM protein in cancerous tissues by depicting strong staining signals. Furthermore, mutation analysis suggested missense as the predominant type of mutation both in LUAD and LUSC in the EpCAM gene. A significant correlation (P-value < 0.05) between the higher EpCAM expression and lower patient survival was also found in this study. Finally, the co-expressed genes were identified with their ontological features and signaling pathways associated in LC development. The overall study suggests EpCAM to be a significant biomarker for human LC prognosis.     

Identification of a novel glycolysis-related gene signature that can predict the survival of patients with lung adenocarcinoma

Lung cancer is one of the most malignant cancers worldwide, and lung adenocarcinoma (LUAD) is the most common histologic subtype. Thousands of biomarkers related to the survival and prognosis of patients with this cancer type have been investigated through database mining; however, the prediction effect of a single gene biomarker is not satisfactorily specific or sensitive. Thus, the present study aimed to develop a novel gene signature of prognostic values for patients with LUAD. Using a data-mining method, we performed expression profiling of 1145 mRNAs in large cohorts with LUAD (n = 511) from The Cancer Genome Atlas database. Using the Gene Set Enrichment Analysis, we selected 198 genes related to GLYCOLYSIS, which is the most important enrichment gene set. Moreover, these genes were identified using Cox proportional regression modeling. We established a risk score staging system to predict the outcome of patients with LUAD and subsequently identified four genes (AGRN, AKR1A1, DDIT4, and HMMR) that were closely related to the prognosis of patients with LUAD. The identified genes allowed us to classify patients into the high-risk group (with poor outcome) and low-risk group (with better outcome). Compared with other clinical factors, the risk score has a better performance in predicting the outcome of patients with LUAD, particularly in the early stage of LUAD. In conclusion, we developed a four-gene signature related to glycolysis by utilizing the Cox regression model and a risk staging model for LUAD, which might prove valuable for the clinical management of patients with LUAD.     

Zinc Finger E-Box Binding Homeobox 2 as a Prognostic Biomarker in Various Cancers and Its Correlation with Infiltrating Immune Cells in Ovarian Cancer

This study investigated the expression of zinc finger E-box binding homeobox 2 (ZEB2), its prognostic significance in various cancers, and the correlation between ZEB2 and infiltrating immune cells and ZEB2-related proteins in ovarian cancer (OV). The Gene Expression Profiling Interactive Analysis tool was used to analyze RNA sequencing data and cancer survival rates, based on normal and tumor tissue data available in The Cancer Genome Atlas (TCGA) database. The Kaplan–Meier plotter and PrognoScan databases were used to analyze the prognostic value of ZEB2 in OV (n = 1144). The Tumor Immune Estimation Resource was used to investigate the correlation between ZEB2 and infiltrating immune cells in various cancers, including OV. High ZEB2 expression was associated with a poorer prognosis in OV. In OV, ZEB2 is positively correlated with CD8+T cells, neutrophils, macrophages, and dendritic cell invasion; and ZEB2 is negatively correlated with tumor-infiltrating B cells. The STRING database was used to investigate the correlations with ZEB2-related proteins. The results reveal that ZEB2 was positively correlated with SMAD1 and SMAD2 in OV. Our findings may serve as a potential prognostic biomarker, and provide novel insights into the tumor immunology in OV. Thus, ZEB2 may be a potential diagnostic and therapeutic target in OV.     

Expression and significance of SOX B1 genes in glioblastoma multiforme patients

The overall survival of glioblastoma multiforme (GBM) patients remains poor. To improve patient outcomes, effective diagnostic and prognostic biomarkers for GBM are needed. In this study, we first applied bioinformatic analyses to identify biomarkers for GBM, focusing on SOX (sex‐determining region on the Y chromosome (SRY)‐related high mobility group (HMG) box) B1 family members. The ONCOMINE, GEPIA, LinkedOmics and CCLE databases were used to assess mRNA expression levels of the SOX B1 family members in different cancers and normal tissue. Further bioinformatic analysis was performed using the ONCOMINE database in combination with the LinkedOmics data set to identify the prognostic value of SOX B1 family members for GBM. We found mRNA expression levels of all tested SOX B1 genes were significantly increased in GBM. In the LinkedOmics database, increased expression of SOX3 indicated a better overall survival. In GEPIA databases, increased expression of all SOX B1 family members suggested an improved overall survival, but none of them were statistically different. Then, Transwell assays and wound healing were employed to evaluate the motility and invasive captivity of U251 cells when silencing SOX2 and SOX3. We found exogenous inhibition of SOX2 appeared to reduce the migration and invasion of U251 cells in vitro. Collectively, our research suggested that SOX2 might serve as a cancer‐promoting gene to identify high‐risk GBM patients, and SOX3 had the potential to be a prognostic biomarker for GBM patients.     

通过生物信息学分析鉴定CYB561与浸润性乳腺癌不良预后相关

通过生物信息学研究细胞色素b561(Cytochromes b561,CYB561)与浸润性乳腺癌(Breast invasive carcinoma,BRCA)不良预后的相关性.CYB561氨基酸序列高度保守,在人体组织中广泛表达,可能预示着广泛的细胞功能.CYB561在多种癌症中的失调,提示其可能参与人类癌症发生发...     

INHA acts as a novel and potential biomarker in lung adenocarcinoma and shapes the immune-suppressive tumor microenvironment

BackgroundINHA expression has been correlated with the development, growth, and progression of multiple cancer types. However, the biological role of INHA has not been investigated in patients with lung adenocarcinoma (LUAD). Here, we performed a comprehensive bioinformatics analysis of the LUAD dataset to determine the mechanisms underlying the regulation of tumorigenesis by INHA.Materials and methodsINHA expression and clinical information of patients with LUAD were obtained from The Cancer Genome Atlas (TCGA) database. Protein levels in LUAD cell lines and human lung epithelial cells were examined by western blotting. Next, the prognostic value of INHA in LUAD was assessed using Cox regression analysis, while the potential biological functions and the impact on the immune microenvironment of INHA were investigated using gene set enrichment analysis (GSEA) and single sample GSEA (ssGSEA). Finally, the effect of INHA on LUAD cell proliferation and invasion was determined in vitro and in vivo.ResultsWe found significantly high mRNA and protein expression levels of INHA in LUAD tissues and cell lines. Additionally, a higher expression of INHA was linked to a shorter overall survival (OS) and a worse pathological stage, while INHA expression was associated with immune cell infiltration and immune-related markers in the LUAD tumor microenvironment. LUAD with high INHA expression tends to be a cold tumor. Furthermore, GO and KEGG enrichment analysis indicated that INHA-related genes were enriched in the cell adhesion and immune signaling pathways of LUAD. INHA promoted LUAD cell proliferation and invasion, in vitro and in vivo, by inducing the EGFR pathway.ConclusionOur findings revealed that INHA is overexpressed in LUAD and is linked to a poor prognosis. Our study demonstrates the potential of INHA as an immunotherapeutic and predictive biomarker in LUAD.     

A balance between activating and repressive histone modifications regulates cystic fibrosis transmembrane conductance regulator (CFTR) expression in vivo

The genetic mechanisms that regulate CFTR, the gene responsible for cystic fibrosis, have been widely investigated in cultured cells. However, mechanisms responsible for tissue-specific and time-specific expression are not completely elucidated in vivo. Through the survey of public databases, we found that the promoter of CFTR was associated with bivalent chromatin in human embryonic stem (ES) cells. In this work, we analyzed fetal (at different stages of pregnancy) and adult tissues and showed that, in digestive and lung tissues, which expressed CFTR, H3K4me3 was maintained in the promoter. Histone acetylation was high in the promoter and in two intronic enhancers, especially in fetal tissues. In contrast, in blood cells, which did not express CFTR, the bivalent chromatin was resolved (the promoter was labeled by the silencing mark H3K27me3). Cis-regulatory sequences were associated with lowly acetylated histones. We also provide evidence that the tissue-specific expression of CFTR is not regulated by dynamic changes of DNA methylation in the promoter. Overall, this work shows that a balance between activating and repressive histone modifications in the promoter and intronic enhancers results in the fine regulation of CFTR expression during development, thereby ensuring tissue specificity.     

Cell-type-specific long-range looping interactions identify distant regulatory elements of the CFTR gene

Identification of regulatory elements and their target genes is complicated by the fact that regulatory elements can act over large genomic distances. Identification of long-range acting elements is particularly important in the case of disease genes as mutations in these elements can result in human disease. It is becoming increasingly clear that long-range control of gene expression is facilitated by chromatin looping interactions. These interactions can be detected by chromosome conformation capture (3C). Here, we employed 3C as a discovery tool for identification of long-range regulatory elements that control the cystic fibrosis transmembrane conductance regulator gene, CFTR. We identified four elements in a 460-kb region around the locus that loop specifically to the CFTR promoter exclusively in CFTR expressing cells. The elements are located 20 and 80 kb upstream; and 109 and 203 kb downstream of the CFTR promoter. These elements contain DNase I hypersensitive sites and histone modification patterns characteristic of enhancers. The elements also interact with each other and the latter two activate the CFTR promoter synergistically in reporter assays. Our results reveal novel long-range acting elements that control expression of CFTR and suggest that 3C-based approaches can be used for discovery of novel regulatory elements.     

ZNF76 predicts prognosis and response to platinum chemotherapy in human ovarian cancer

Ovarian cancer (OV) is the most lethal gynecologic malignancy. One major reason of the high mortality of the disease is due to platinum-based chemotherapy resistance. Increasing evidence reveal the important biological functions and clinical significance of zinc finger proteins (ZNFs) in OV. In the present study, the relationship between the zinc finger protein 76 (ZNF76) and clinical outcome and platinum resistance in patients with OV was explored. We further analyzed ZNF76 expression via multiple gene expression databases and identified its functional networks using cBioPortal. RT-qPCR and IHC assay shown that the ZNF76 mRNA and protein expression were significantly lower in OV tumor than that in normal ovary tissues. A strong relationship between ZNF76 expression and platinum resistance was determined in patients with OV. The low expression of ZNF76 was associated with worse survival in OV. Multivariable analysis showed that the low expression of ZNF76 was an independent factor predicting poor outcome in OV. The prognosis value of ZNF76 in pan-cancer was validated from multiple cohorts using the PrognoScan database and GEPIA 2. A gene-clinical nomogram was constructed by multivariate cox regression analysis, combined with clinical characterization and ZNF76 expression in TCGA. Functional network analysis suggested that ZNF76 was involved in several biology progressions which associated with OV. Ten hub genes (CDC5L, DHX16, SNRPC, LSM2, CUL7, PFDN6, VARS, HSD17B8, PPIL1, and RGL2) were identified as positively associated with the expression of ZNF76 in OV. In conclusion, ZNF76 may serve as a promising prognostic-related biomarker and predict the response to platinum in OV patients.     

ANKRD49 promotes the invasion and metastasis of lung adenocarcinoma via a P38/ATF‐2 signalling pathway

Lung adenocarcinoma (LUAD) is the most challenging neoplasm to treat in clinical practice. Ankyrin repeat domain 49 protein (ANKRD49) is highly expressed in several carcinomas; however, its pattern of expression and role in LUAD are not known. Tissue microarrays, immunohistochemistry, χ ² test, Spearman correlation analysis, Kaplan–Meier, log‐rank test, and Cox''s proportional hazard model were used to analyse the clinical cases. The effect of ANKRD49 on the LUAD was investigated using CCK‐8, clonal formation, would healing, transwell assays, and nude mice experiment. Expressions of ANKRD49 and its associated downstream protein molecules were verified by real‐time PCR, Western blot, immunohistochemistry, and/or immunofluorescence analyses. ANKRD49 expression was highly elevated in LUAD. The survival rate and Cox''s modelling analysis indicated that there may be an independent prognostic indicator for LUAD patients. We also found that ANKRD49 promoted the invasion and migration in both in in vitro and in vivo assays, through upregulating matrix metalloproteinase (MMP)‐2 and MMP‐9 activities via the P38/ATF‐2 signalling pathway Our findings suggest that ANKRD49 is a latent biomarker for evaluating LUAD prognosis and promotes the metastasis of A549 cells via upregulation of MMP‐2 and MMP‐9 in a P38/ATF‐2 pathway‐dependent manner.     

A Genotypic-Oriented View of CFTR Genetics Highlights Specific Mutational Patterns Underlying Clinical Macrocategories of Cystic Fibrosis

Cystic fibrosis (CF) is a monogenic disease caused by mutations of the cystic fibrosis transmembrane conductance regulator (CFTR) gene. The genotype–phenotype relationship in this disease is still unclear, and diagnostic, prognostic and therapeutic challenges persist. We enrolled 610 patients with different forms of CF and studied them from a clinical, biochemical, microbiological and genetic point of view. Overall, there were 125 different mutated alleles (11 with novel mutations and 10 with complex mutations) and 225 genotypes. A strong correlation between mutational patterns at the genotypic level and phenotypic macrocategories emerged. This specificity appears to largely depend on rare and individual mutations, as well as on the varying prevalence of common alleles in different clinical macrocategories. However, 19 genotypes appeared to underlie different clinical forms of the disease. The dissection of the pathway from the CFTR mutated genotype to the clinical phenotype allowed to identify at least two components of the variability usually found in the genotype–phenotype relationship. One component seems to depend on the genetic variation of CFTR, the other component on the cumulative effect of variations in other genes and cellular pathways independent from CFTR. The experimental dissection of the overall biological CFTR pathway appears to be a powerful approach for a better comprehension of the genotype–phenotype relationship. However, a change from an allele-oriented to a genotypic-oriented view of CFTR genetics is mandatory, as well as a better assessment of sources of variability within the CFTR pathway.     

RHOV promotes lung adenocarcinoma cell growth and metastasis through JNK/c-Jun pathway

Lung adenocarcinoma (LUAD) is a common type of lung cancer with high frequent metastasis and a high death rate. However, genes responsible for LUAD metastasis are still largely unknown. Here, we identify an important role of ras homolog family member V (RHOV) in LUAD metastasis using a combination of bioinformatic analysis and functional experiments. Bioinformatic analysis shows five hub LUAD metastasis driver genes (RHOV, ZIC5, CYP4B1, GPR18 and TCP10L2), among which RHOV is the most significant gene associated with LUAD metastasis. High RHOV expression predicted shorter overall survival in LUAD patients. RHOV overexpression promotes proliferation, migration, and invasion of LUAD cells, whereas RHOV knockdown inhibits these biological behaviors. Moreover, knockdown of RHOV suppresses LUAD tumor growth and metastasis in nude mice. Mechanistically, RHOV activates Jun N-terminal Kinase (JNK)/c-Jun signalling pathway, an important pathway in lung cancer development and progression, and regulates the expression of markers of epithelial-to-mesenchymal transition, a process involved in cancer cell migration, invasion and metastasis. RHOV-induced malignant biological behaviors are inhibited by pyrazolanthrone, a JNK inhibitor. Our findings indicate a critical role of RHOV in LUAD metastasis and may provide a biomarker for prognostic prediction and a target for LUAD therapy.     

Identification of KIF4A and its effect on the progression of lung adenocarcinoma based on the bioinformatics analysis

Background: Lung adenocarcinoma (LUAD) is the most frequent histological type of lung cancer, and its incidence has displayed an upward trend in recent years. Nevertheless, little is known regarding effective biomarkers for LUAD.Methods: The robust rank aggregation method was used to mine differentially expressed genes (DEGs) from the gene expression omnibus (GEO) datasets. The Search Tool for the Retrieval of Interacting Genes (STRING) database was used to extract hub genes from the protein–protein interaction (PPI) network. The expression of the hub genes was validated using expression profiles from TCGA and Oncomine databases and was verified by real-time quantitative PCR (qRT-PCR). The module and survival analyses of the hub genes were determined using Cytoscape and Kaplan–Meier curves. The function of KIF4A as a hub gene was investigated in LUAD cell lines.Results: The PPI analysis identified seven DEGs including BIRC5, DLGAP5, CENPF, KIF4A, TOP2A, AURKA, and CCNA2, which were significantly upregulated in Oncomine and TCGA LUAD datasets, and were verified by qRT-PCR in our clinical samples. We determined the overall and disease-free survival analysis of the seven hub genes using GEPIA. We further found that CENPF, DLGAP5, and KIF4A expressions were positively correlated with clinical stage. In LUAD cell lines, proliferation and migration were inhibited and apoptosis was promoted by knocking down KIF4A expression.Conclusion: We have identified new DEGs and functional pathways involved in LUAD. KIF4A, as a hub gene, promoted the progression of LUAD and might represent a potential therapeutic target for molecular cancer therapy.