Genetic and epigenetic analysis of non-small cell lung cancer with NotI-microarrays

Study on Potential Differentially Expressed Genes in Idiopathic Pulmonary Fibrosis by Bioinformatics and Next-Generation Sequencing Data Analysis

Idiopathic pulmonary fibrosis (IPF) is a chronic progressive lung disease with reduced quality of life and earlier mortality, but its pathogenesis and key genes are still unclear. In this investigation, bioinformatics was used to deeply analyze the pathogenesis of IPF and related key genes, so as to investigate the potential molecular pathogenesis of IPF and provide guidance for clinical treatment. Next-generation sequencing dataset GSE213001 was obtained from Gene Expression Omnibus (GEO), and the differentially expressed genes (DEGs) were identified between IPF and normal control group. The DEGs between IPF and normal control group were screened with the DESeq2 package of R language. The Gene Ontology (GO) and REACTOME pathway enrichment analyses of the DEGs were performed. Using the g:Profiler, the function and pathway enrichment analyses of DEGs were performed. Then, a protein-protein interaction (PPI) network was constructed via the Integrated Interactions Database (IID) database. Cytoscape with Network Analyzer was used to identify the hub genes. miRNet and NetworkAnalyst databaseswereused to construct the targeted microRNAs (miRNAs), transcription factors (TFs), and small drug molecules. Finally, receiver operating characteristic (ROC) curve analysis was used to validate the hub genes. A total of 958 DEGs were screened out in this study, including 479 up regulated genes and 479 down regulated genes. Most of the DEGs were significantly enriched in response to stimulus, GPCR ligand binding, microtubule-based process, and defective GALNT3 causes HFTC. In combination with the results of the PPI network, miRNA-hub gene regulatory network and TF-hub gene regulatory network, hub genes including LRRK2, BMI1, EBP, MNDA, KBTBD7, KRT15, OTX1, TEKT4, SPAG8, and EFHC2 were selected. Cyclothiazide and rotigotinethe are predicted small drug molecules for IPF treatment. Our findings will contribute to identification of potential biomarkers and novel strategies for the treatment of IPF, and provide a novel strategy for clinical therapy.

Contribution of gene polymorphisms on 3p25 to salivary gland carcinoma, ameloblastoma, and odontogenic keratocyst in the Chinese Han population

OBJECTIVE

This study aimed to investigate the contribution of gene polymorphisms in 3p25 to salivary gland carcinoma (SGC), ameloblastoma (AM), and odontogenic keratocyst (OKC) in the Chinese Han population.

STUDY DESIGN

Sixteen tag-single nucleotide polymorphisms (SNPs) within 5 genes (SYN2, TIMP4, PPARG, RAF1, and IQSEC1) in 3p25 were genotyped in 411 individuals with or without SGC, AM, and OKC. Genotype, clinical phenotype, and bioinformatics analyses were performed to evaluate the function of candidate SNPs.

RESULTS

SYN2-rs3773364, TIMP4-rs3755724, PPARG-rs10865710, and PPARG-rs1175544 were related to decreased SGC susceptibility, whereas IQSEC1-rs2600322 and IQSEC1-rs2686742 decreased and increased AM risk, respectively. Stratification analysis revealed that the significance of the identified SNPs was stronger in females or individuals younger than 46 years in SGC. PPARG-rs10865710 and PPARG-rs1175544 were associated with lower lymph node metastasis. SYN2-rs3773364 and PPARG-rs1175544 were associated with favorable SGC patient survival. Functional assessments linked PPARG-rs1175544 to PPARG expression regulation. Linkage disequilibrium analysis revealed a haplotype (SYN2-rs3773364-A, TIMP4-rs3817004-A, and TIMP4-rs3755724-C) associated with decreased susceptibility to SGC. Generalized multifactor dimensionality reduction analysis indicated the gene-gene interactions among IQSEC1, TIMP4, and PPARG in SGC, AM, and OKC progression.

CONCLUSIONS

These variants play important roles in the progression of SGC, AM, and OKC in the Chinese Han population and may be considered biomarkers for early diagnosis and prognosis prediction.

Systematic Analysis of the Prognostic Significance and Roles of the Integrin Alpha Family in Non-Small Cell Lung Cancers

INTRODUCTION

Lung cancer is one of the most common cancer malignancies and the principal cause of cancer-associated deaths worldwide. Non-small cell lung cancers (NSCLCs) account for more than 80% of all lung cancer cases. Recent studies showed that the genes of the integrin alpha (α) (ITGA) subfamily play a fundamental role in various cancers. However, little is known about the expression and roles of distinct ITGA proteins in NSCLCs.

METHODS

Gene Expression Profiling Interactive Analysis and UALCAN (University of ALabama at Birmingham CANcer) web resources and The Cancer Genome Atlas (TCGA), ONCOMINE, cBioPortal, GeneMANIA, and Tumor Immune Estimation Resource databases were used to evaluate differential expression, correlations between the expression levels of individual genes, the prognostic value of overall survival (OS) and stage, genetic alterations, protein-protein interactions, and the immune cell infiltration of ITGAs in NSCLCs. We used R (v. 4.0.3) software to conduct gene correlation, gene enrichment, and clinical correlation of RNA sequencing data of 1016 NSCLCs from TCGA. To evaluate the expression of ITGA5/8/9/L at the expression and protein levels, qRT-PCR, immunohistochemistry (IHC), and hematoxylin and eosin (H&E) were performed, respectively.

RESULTS

Upregulated levels of ITGA11 messenger RNA and downregulated levels of ITGA1/3/5/7/8/9/L/M/X were observed in the NSCLC tissues. Lower expression of ITGA5/6/8/9/10/D/L was discovered to be expressively associated with advanced tumor stage or poor patient prognosis in patients with NSCLC. A high mutation rate (44%) of the ITGA family was observed in the NSCLCs. Gene Ontology functional enrichment analyses results revealed that the differentially expressed ITGAs could be involved in roles related to extracellular matrix (ECM) organization, collagen-containing ECM cellular components, and ECM structural constituent molecular functions. The results of the Kyoto Encyclopedia of Genes and Genomes analysis revealed that ITGAs may be involved in focal adhesion, ECM-receptor interaction, and amoebiasis; the expression of ITGAs was significantly correlated with the infiltration of diverse immune cells in NSCLCs. ITGA5/8/9/L was also highly correlated with PD-L1 expression. The validation results for marker gene expression in NSCLC tissues by qRT-PCR, IHC, and H&E staining indicated that the expression of ITGA5/8/9/L decreased compared with that in normal tissues.

CONCLUSION

As potential prognostic biomarkers in NSCLCs, ITGA5/8/9/L may fulfill important roles in regulating tumor progression and immune cell infiltration.

Exploratory Metabolomics Underscores the Folate Enzyme ALDH1L1 as a Regulator of Glycine and Methylation Reactions

Folate (vitamin B9) is involved in one-carbon transfer reactions and plays a significant role in nucleic acid synthesis and control of cellular proliferation, among other key cellular processes. It is now recognized that the role of folates in different stages of carcinogenesis is complex, and more research is needed to understand how folate reactions become dysregulated in cancers and the metabolic consequences that occur as a result. ALDH1L1 (cytosolic 10-formyltetrahydrofolate dehydrogenase), an enzyme of folate metabolism expressed in many tissues, is ubiquitously downregulated in cancers and is not expressed in cancer cell lines. The RT4 cell line (derived from papillary bladder cancer) which expresses high levels of ALDH1L1 represents an exception, providing an opportunity to explore the metabolic consequences of the loss of this enzyme. We have downregulated this protein in RT4 cells (shRNA driven knockdown or CRISPR driven knockout) and compared metabolomes of ALDH1L1-expressing and -deficient cells to determine if metabolic changes linked to the loss of this enzyme might provide proliferative and/or survival advantages for cancer cells. In this study, cell extracts were analyzed using Ultra High Performance Liquid Chromatography High Resolution Mass Spectrometry (UHPLC-HR-MS). A total of 13,339 signals were identified or annotated using an in-house library and public databases. Supervised and unsupervised multivariate analysis revealed metabolic differences between RT4 cells and ALDH1L1-deficient clones. Glycine (8-fold decrease) and metabolites derived from S-adenosylmethionine utilizing pathways were significantly decreased in the ALDH1L1-deficient clones, compared with RT4 cells. Other changes linked to ALDH1L1 downregulation include decreased levels of amino acids, Krebs cycle intermediates, and ribose-5-phosphate, and increased nicotinic acid. While the ALDH1L1-catalyzed reaction is directly linked to glycine biosynthesis and methyl group flux, its overall effect on cellular metabolism extends beyond immediate metabolic pathways controlled by this enzyme.

Proteomics analysis of chronic skin injuries caused by mustard gas

Sulfur mustard (SM) is an alkylating and forming chemical that was widely used by Iraqi forces during the Iran–Iraq wars. One of the target organs of SM is the skin. Understanding the mechanisms involved in the pathogenesis of SM may help better identify complications and find appropriate treatments. The current study collected ten SM-exposed patients with long-term skin complications and ten healthy individuals. Proteomics experiments were performed using the high-efficiency TMT10X method to evaluate the skin protein profile, and statistical bioinformatics methods were used to identify the differentially expressed proteins. One hundred twenty-nine proteins had different expressions between the two groups. Of these 129 proteins, 94 proteins had increased expression in veterans' skins, while the remaining 35 had decreased expression. The hub genes included RPS15, ACTN1, FLNA, HP, SDHC, and RPL29, and three modules were extracted from the PPI network analysis. Skin SM exposure can lead to oxidative stress, inflammation, apoptosis, and cell proliferation.

Sex-Specific Metabolic Effects of Dietary Folate Withdrawal in Wild-Type and Aldh1l1 Knockout Mice

ALDH1L1 (10-formyltetrahydrofolate dehydrogenase), an enzyme of folate metabolism, is highly expressed in the liver. It regulates the overall flux of folate-bound one-carbon groups by converting 10-formyltetrahydrofolate to tetrahydrofolate and CO2 in a NADP+-dependent reaction. Our previous study revealed that Aldh1l1 knockout (KO) mice have an altered liver metabotype with metabolic symptoms of folate deficiency when fed a standard chow diet containing 2 ppm folic acid. Here we performed untargeted metabolomic analysis of liver and plasma of KO and wild-type (WT) male and female mice fed for 16 weeks either standard or folate-deficient diet. OPLS-DA, a supervised multivariate technique that was applied to 6595 and 10,678 features for the liver and plasma datasets, respectively, indicated that genotype and diet, alone or in combination, gave distinct metabolic profiles in both types of biospecimens. A more detailed analysis of affected metabolic pathways based on most confidently identified metabolites in the liver and plasma (OL1 and OL2a ontology level) indicated that the dietary folate restriction itself does not fully recapitulate the metabolic effect of the KO. Of note, dietary folate withdrawal enhanced the metabolic perturbations linked to the ALDH1L1 loss only for a subset of metabolites. Importantly, both the ALDH1L1 loss and dietary folate deficiency produced sex-specific metabolic effects.

ITGA9: Potential Biomarkers and Therapeutic Targets in Different Tumors

Abstract:

Integrins are a class of a cell surface adhesion molecule which composed of α subunit (ITGA) and β subunit (ITGB). They belong to heterodimer transmembrane glycoproteins. Its main function in organisms is as the receptor of cell adhesion molecules (CAMs) and extracellular matrix (ECM). According to the current research integration analysis, integrin α9 (ITGA9) is one of the integrin subunits, and there are few studies on ITGA9 among integrins. ITGA9 can improve cell migration and regulate various cellular biological functions, such as tumor cell proliferation, adhesion, invasion, and angiogenesis. But its abnormal expression mechanism in cancer and its specific role in tumor growth and metastasis are still unknown to a great extent. This review reveals the role of ITGA9 in the complex pathogenesis of many tumors and cancers, providing a new direction for the treatment of tumors and cancers. Relevant studies were retrieved and collected through the PubMed system. After determining ITGA9 as the research object, we found the close relationship between ITGA9 and tumorigenesis through the analysis of the research articles on ITGA9 in the PubMed system in the last 15 years, and further determined the references mainly based on the influencing factors of the articles. Thus, the role of ITGA9 in tumor and cancer genesis, proliferation, and metastasis was reviewed and analyzed. ITGA9 is an integrin subunit, which has been proved to be abnormally expressed in many tumors. After sorting and analyzing the research data, it was found that the abnormal expression of ITGA9 in a variety of tumors, including glioblastoma, rhabdomyosarcoma, melanoma, hepatocellular carcinoma, nasopharyngeal carcinoma, multiple myeloma, non-small cell lung cancer, and prostate cancer, was closely related to the proliferation, metastasis, adhesion, and angiogenesis of tumor cells. These results suggest that ITGA9 plays an important role in the occurrence and development of tumors. The integrin subunit ITGA9 may serve as a biomarker for the diagnosis of tumors and a potential therapeutic target for anti-tumor therapies.

Structure of putative tumor suppressor ALDH1L1

Putative tumor suppressor ALDH1L1, the product of natural fusion of three unrelated genes, regulates folate metabolism by catalyzing NADP⁺-dependent conversion of 10-formyltetrahydrofolate to tetrahydrofolate and CO₂. Cryo-EM structures of tetrameric rat ALDH1L1 revealed the architecture and functional domain interactions of this complex enzyme. Highly mobile N-terminal domains, which remove formyl from 10-formyltetrahydrofolate, undergo multiple transient inter-domain interactions. The C-terminal aldehyde dehydrogenase domains, which convert formyl to CO₂, form unusually large interfaces with the intermediate domains, homologs of acyl/peptidyl carrier proteins (A/PCPs), which transfer the formyl group between the catalytic domains. The 4'-phosphopantetheine arm of the intermediate domain is fully extended and reaches deep into the catalytic pocket of the C-terminal domain. Remarkably, the tetrameric state of ALDH1L1 is indispensable for catalysis because the intermediate domain transfers formyl between the catalytic domains of different protomers. These findings emphasize the versatility of A/PCPs in complex, highly dynamic enzymatic systems.

Knockout of Putative Tumor Suppressor Aldh1l1 in Mice Reprograms Metabolism to Accelerate Growth of Tumors in a Diethylnitrosamine (DEN) Model of Liver Carcinogenesis

Simple Summary

Cancers often loose the enzyme of folate metabolism ALDH1L1. We proposed that such loss is advantageous for the malignant tumor growth and tested this hypothesis in mice proficient or deficient (gene knockout) in ALDH1L1 expression. Liver cancer in both groups was induced by injection of chemical carcinogen diethylnitrosamine. While the number of tumors observed in ALDH1L1 proficient and deficient mice was similar, tumors grew faster and to a larger size in the knockout mice. We conclude that the ALDH1L1 loss promotes liver tumor growth without affecting tumor initiation or multiplicity. Accelerated growth of tumors lacking the enzyme was linked to several metabolic pathways, which are beneficial for rapid proliferation.

Abstract

Cytosolic 10-formyltetrahydrofolate dehydrogenase (ALDH1L1) is commonly downregulated in human cancers through promoter methylation. We proposed that ALDH1L1 loss promotes malignant tumor growth. Here, we investigated the effect of the Aldh1l1 mouse knockout (Aldh1l1^−/−) on hepatocellular carcinoma using a chemical carcinogenesis model. Fifteen-day-old male Aldh1l1 knockout mice and their wild-type littermate controls (Aldh1l1^+/+) were injected intraperitoneally with 20 μg/g body weight of DEN (diethylnitrosamine). Mice were sacrificed 10, 20, 28, and 36 weeks post-DEN injection, and livers were examined for tumor multiplicity and size. We observed that while tumor multiplicity did not differ between Aldh1l1^−/− and Aldh1l1^+/+ animals, larger tumors grew in Aldh1l1^−/− compared to Aldh1l1^+/+ mice at 28 and 36 weeks. Profound differences between Aldh1l1^−/− and Aldh1l1^+/+ mice in the expression of inflammation-related genes were seen at 10 and 20 weeks. Of note, large tumors from wild-type mice showed a strong decrease of ALDH1L1 protein at 36 weeks. Metabolomic analysis of liver tissues at 20 weeks showed stronger differences in Aldh1l1^+/+ versus Aldh1l1^−/− metabotypes than at 10 weeks, which underscores metabolic pathways that respond to DEN in an ALDH1L1-dependent manner. Our study indicates that Aldh1l1 knockout promoted liver tumor growth without affecting tumor initiation or multiplicity.

Methylation Regulation of TLR3 on Immune Parameters in Lung Adenocarcinoma

This study aims to analyze the methylation regulation of TLR3 in lung adenocarcinoma (LUAD) and to explore the association of TLR3 expression with immune microenvironment. TLR3 has a decreased expression in LUAD tissues and low expression of TLR3 is not only associated with poor prognosis in patients with LUAD, but also can be used as a diagnostic marker. Bisulfite sequencing PCR (BSP) results showed that the methylation level in the promoter of TLR3 was negatively correlated with the level of TLR3 mRNA in LUAD tissues. TIMER analysis showed that TLR3 was negatively correlated with the tumor purity of LUAD and positively with immune cell infiltration to some extent. ESTIMATE analysis also suggested that TLR3 expression and its methylation had significant correlation with immune score. The lower immune scores were associated with the late stage of LUAD and poor prognosis. The high expression of TLR3 might inhibit the development of LUAD by activating apoptosis pathway. The proteins interacted with TLR3 were mainly involved in the apoptosis pathway and positively correlated with the key genes (MYD88, Caspase 8, BIRC3, PIK3R1) in this pathway. Therefore, TLR3 as a key biomarker for prognosis and diagnosis in LUAD, might be considered as a potential epigenetic and immunotherapeutic target.

Long Non-Coding RNA HOXA11-AS Modulates Proliferation, Apoptosis, Metastasis and EMT in Cutaneous Melanoma Cells Partly via miR-152-3p/ITGA9 Axis

Background

Long non-coding RNA homeobox A11 antisense RNA (HOXA11-AS) was showed to participate in the progression of different kinds of tumors, but the specific role of HOXA11-AS in cutaneous melanoma is not entirely unambiguous.

Methods

The levels of HOXA11-AS, microRNA-152-3p (miR-152-3p) and integrin alpha9 (ITGA9) were measured by quantitative real-time polymerase chain reaction (qRT-PCR). Cell proliferation was detected via 3-(4, 5-dimethylthiazol-2-y1)-2, 5-diphenyl tetrazolium bromide (MTT), and apoptosis was measured by flow cytometry. The assessment of cell metastasis was performed by transwell migration and invasion assays. The protein levels were detected through Western blot. Dual-luciferase reporter assay was utilized to explore the target relationship among HOXA11-AS, miR-152-3p and ITGA9. The effect of HOXA11-AS on melanoma in vivo was investigated via xenograft experiment.

Results

HOXA11-AS and ITGA9 were up-regulated while miR-152-3p was down-regulated in melanoma. Knockdown of HOXA11-AS refrained cell proliferation, metastasis and epithelial-mesenchymal transition (EMT) but induced apoptosis in melanoma cells. HOXA11-AS targeted miR-152-3p and overexpression of HOXA11-AS mitigated the miR-152-3p-induced effects on melanoma cellular behaviors. ITGA9 was a target of miR-152-3p and miR-152-3p inhibitor relieved the repression on proliferation, metastasis and EMT while elevation on apoptosis caused by si-ITGA9 via elevating ITGA9. HOXA11-AS knockdown restrained ITGA9 expression via up-regulating miR-152-3p. Suppression of HOXA11-AS inhibited melanoma progression in part through increasing miR-152-3p and decreasing ITGA9 expression in vivo.

Conclusion

HOXA11-AS modulated proliferation, apoptosis, metastasis and EMT in melanoma cells by regulating miR-152-3p/ITGA9 axis in part. HOXA11-AS could promote melanoma development and be used as a promising biomarker in the diagnosis and treatment for cutaneous melanoma.

Methylation biomarkers with discriminating ability are potential therapeutic targets in lung adenocarcinoma

Aims: Given the reversibility of methylation, biomarkers with discriminating ability are of great interest for targeted therapeutic sites. Materials & methods: Methylation array data of 461 lung adenocarcinoma (LUAD) patients comprising of 458 tumor and 32 LUAD paracancerous samples were compared using partial least squares discrimination analysis and receiver operating characteristics analysis. Results: A six-DNA methylation signature (corresponding to five genes) was found to significantly discriminate normal and LUAD samples. Kyoto Encyclopedia of Genes and Genomes analysis indicated enrichment of methylation sites in the Wnt pathway in LUAD compared with controls. Conclusion: This six-DNA methylation signature demonstrated potential as a novel biomarker for diagnosis and therapeutic targets. Further, inhibition of Wnt signaling pathway may be an important step in LUAD progression.

Tumor suppressor properties of the small C-terminal domain phosphatases in non-small cell lung cancer

Non-Small Cell Lung Cancer (NSCLC) is responsible for the majority of deaths caused by cancer. Small C-terminal domain (CTD) phosphatases (SCP), CTDSP1, CTDSP2 and CTDSPL (CTDSPs) belong to SCP/CTDSP subfamily and are involved in many vital cellular processes and tumorigenesis. High similarity of their structures suggests similar functions. However their role in NSCLC remains insufficiently understood. For the first time we revealed the suppressor function of CTDSPs leading to a significant growth slowdown and senescence of A549 lung adenocarcinoma (ADC) cells in vitro. Their tumor-suppressive activity can be realized through increasing the proportion of the active form of Rb protein dephosphorylated at Ser807/811, Ser780, and Ser795 (P<0.05) thereby negatively regulating cancer cell proliferation. Moreover, we observed that a frequent (84%, 39/46) and highly concordant (Spearman's rank correlation coefficient (rs) = 0.53-0.62, P≤0.01) down-regulation of CTDSPs and RB1 is characteristic of primary NSCLC samples (n=46). A clear difference in their mRNA levels was found between lung ADCs with and without lymph node metastases, but not in squamous cell carcinomas (SCCs) (P≤0.05). Based on The Cancer Genome Atlas (TCGA) data and the results obtained using the CrossHub tool, we suggest that the well-known oncogenic cluster miR-96/182/183 could be a common expression regulator of CTDSPs. Indeed, according to our qPCR, the expression of CTDSPs negatively correlates with these miRs, but positively correlates with their intronic miR-26a/b. Our results reflect functional association of CTDSP1, CTDSP2, and CTDSPL, expand knowledge about their suppressor properties through Rb dephosphorylation and provide new insights into the regulation of NSCLC growth.

Biological effect of ribosomal protein L32 on human breast cancer cell behavior

Breast cancer (BC) is the most common malignancy among women worldwide. However, identifying effective biomarkers for the diagnosis and treatment of BC is challenging. Based on our previously developed 'humanized' mouse model of BC, microarray expression analysis was performed and multiple differentially expressed genes, including ribosomal protein (RP) L32, were screened. Recent reports have revealed that RPs are relevant to the development and progression of cancer. However, the expression and function of RPL32 in BC remains unknown. Therefore, in the present study, the role of RPL32 in the development of BC was explored. Immunohistochemical staining and reverse transcription‑quantitative PCR were used, and it was found that RPL32 was upregulated in human BC tissues and cells. Cell Counting Kit‑8, cell invasion and migration assays were performed, which demonstrated that RPL32 knockdown using lentivirus‑delivered small interfering RNA inhibited the migration and invasion of BC cells in vitro and in vivo (nude mouse model). Moreover, western blotting showed that RPL32 knockdown decreased the expression levels of matrix metalloproteinase (MMP)‑2 and MMP‑9. Thus, the present findings indicated a potential oncogenic role of RPL32, suggesting that it may be a novel target for molecular targeted therapy in patients with BC.

An Immune-Related Signature Predicts Survival in Patients With Lung Adenocarcinoma

We investigated the local immune status and its prognostic value in lung adenocarcinoma. In total, 513 lung adenocarcinoma samples from TCGA and ImmPort databases were collected and analyzed. The R package coxph was employed to mine immune-related genes that were significant prognostic indicators using both univariate and multivariate analyses. The R software package glmnet was then used for Lasso Cox regression analysis, and a prognosis prediction model was constructed for lung adenocarcinoma; clusterProfiler was selected for functional gene annotations and KEGG enrichment analysis. Finally, correlations between the RiskScore and clinical features or signaling pathways were established. Sixty-four immune-related genes remarkably correlated with patient prognosis and were further applied. Samples were hierarchically clustered into two subgroups. Accordingly, the LASSO regression algorithm was employed to screen the 14 most representative immune-related genes (PSMD11, PPIA, MIF, BMP5, DKK1, PDGFB, ANGPTL4, IL1R2, THRB, LTBR, TNFRSF1, TNFRSF17, IL20RB, and MC1R) with respect to patient prognosis. Then, the prognosis prediction model for lung adenocarcinoma patients (namely, the RiskScore equation) was constructed, and the training set samples were incorporated to evaluate the efficiency of this model to predict and classify patient prognosis. Subsequently, based on functional annotations and KEGG pathway analysis, the 14 immune-related genes were mainly enriched in pathways closely associated with lung adenocarcinoma and its immune microenvironment, such as cytokine–cytokine receptor interaction and human T-cell leukemia virus 1 infection. Furthermore, correlations between the RiskScore and clinical features of the training set samples and signaling pathways (such as p53, cell cycle, and DNA repair) were also demonstrated. Finally, the test set sample data were employed for independent testing and verifying the model. We established a prognostic prediction RiskScore model based on the expression profiles of 14 immune-related genes, which shows high prediction accuracy and stability in identifying immune features. This could provide clinical guidance for the diagnosis and prognosis of different immunophenotypes, and suggest multiple targets for precise advanced lung adenocarcinoma therapy based on subtype-specific immune molecules.

The Role of Single-Nucleotide Polymorphisms in the Function of Candidate Tumor Suppressor ALDH1L1

Folate (vitamin B9) is a common name for a group of coenzymes that function as carriers of chemical moieties called one-carbon groups in numerous biochemical reactions. The combination of these folate-dependent reactions constitutes one-carbon metabolism, the name synonymous to folate metabolism. Folate coenzymes and associated metabolic pathways are vital for cellular homeostasis due to their key roles in nucleic acid biosynthesis, DNA repair, methylation processes, amino acid biogenesis, and energy balance. Folate is an essential nutrient because humans are unable to synthesize this coenzyme and must obtain it from the diet. Insufficient folate intake can ultimately increase risk of certain diseases, most notably neural tube defects. More than 20 enzymes are known to participate in folate metabolism. Single-nucleotide polymorphisms (SNPs) in genes encoding for folate enzymes are associated with altered metabolism, changes in DNA methylation and modified risk for the development of human pathologies including cardiovascular diseases, birth defects, and cancer. ALDH1L1, one of the folate-metabolizing enzymes, serves a regulatory function in folate metabolism restricting the flux of one-carbon groups through biosynthetic processes. Numerous studies have established that ALDH1L1 is often silenced or strongly down-regulated in cancers. The loss of ALDH1L1 protein positively correlates with the occurrence of malignant tumors and tumor aggressiveness, hence the enzyme is viewed as a candidate tumor suppressor. ALDH1L1 has much higher frequency of non-synonymous exonic SNPs than most other genes for folate enzymes. Common SNPs at the polymorphic loci rs3796191, rs2886059, rs9282691, rs2276724, rs1127717, and rs4646750 in ALDH1L1 exons characterize more than 97% of Europeans while additional common variants are found in other ethnic populations. The effects of these SNPs on the enzyme is not clear but studies indicate that some coding and non-coding ALDH1L1 SNPs are associated with altered risk of certain cancer types and it is also likely that specific haplotypes define the metabolic response to dietary folate. This review discusses the role of ALDH1L1 in folate metabolism and etiology of diseases with the focus on non-synonymous coding ALDH1L1 SNPs and their effects on the enzyme structure/function, metabolic role and association with cancer.

Mutual suppression between BHLHE40/BHLHE41 and the MIR301B-MIR130B cluster is involved in epithelial-to-mesenchymal transition of endometrial cancer cells

BHLHE40 and BHLHE41 (BHLHE40/41) are basic helix-loop-helix type transcription factors involved in multiple cell activities including epithelial-to-mesenchymal transition (EMT). However, the expression mechanism of BHLHE40/41 in EMT remains unclear. In the present study, we showed that the expression levels of BHLHE40/41 were negatively correlated with those of the microRNA (MIR) 130 family in endometrial cancer (EC) specimens. Our in vitro assays indicated that the expression of BHLHE40/41 was suppressed directly by the MIR130 family in a 3'-untranslated region-mediated manner. In EC cells, the MIR130 family promoted EMT and tumor cell invasion by suppressing the expression of BHLHE40/41. We identified the critical promoter region of the MIR301B-MIR130B cluster for its basal transcription by the transcription factor, SP1. We also found that BHLHE40/41 suppressed the expression of MIR301B and MIR130B, and we identified a binding site in the promoter region for BHLHE40/41. This study is the first to report that BHLHE40/41 and the MIR301B-MIR130B cluster suppressed each other to regulate EMT and invasion of EC cells. We propose that BHLHE40/41 and the MIR130 family are excellent markers to predict the progression of EC cases, and that molecular therapy targeting the MIR130 family-BHLHE40/41 axis may effectively control EC extension.

Loss of ALDH1L1 folate enzyme confers a selective metabolic advantage for tumor progression

ALDH1L1 (cytosolic 10-formyltetrahydrofolate dehydrogenase) is the enzyme in folate metabolism commonly downregulated in human cancers. One of the mechanisms of the enzyme downregulation is methylation of the promoter of the ALDH1L1 gene. Recent studies underscored ALDH1L1 as a candidate tumor suppressor and potential marker of aggressive cancers. In agreement with the ALDH1L1 loss in cancer, its re-expression leads to inhibition of proliferation and to apoptosis, but also affects migration and invasion of cancer cells through a specific folate-dependent mechanism involved in invasive phenotype. A growing body of literature evaluated the prognostic value of ALDH1L1 expression for cancer disease, the regulatory role of the enzyme in cellular proliferation, and associated metabolic and signaling cellular responses. Overall, there is a strong indication that the ALDH1L1 silencing provides metabolic advantage for tumor progression at a later stage when unlimited proliferation and enhanced motility become critical processes for the tumor expansion. Whether the ALDH1L1 loss is involved in tumor initiation is still an open question.

CHIP E3 ligase mediates proteasomal degradation of the proliferation regulatory protein ALDH1L1 during the transition of NIH3T3 fibroblasts from G0/G1 to S-phase

ALDH1L1 is a folate-metabolizing enzyme abundant in liver and several other tissues. In human cancers and cell lines derived from malignant tumors, the ALDH1L1 gene is commonly silenced through the promoter methylation. It was suggested that ALDH1L1 limits proliferation capacity of the cell and thus functions as putative tumor suppressor. In contrast to cancer cells, mouse cell lines NIH3T3 and AML12 do express the ALDH1L1 protein. In the present study, we show that the levels of ALDH1L1 in these cell lines fluctuate throughout the cell cycle. During S-phase, ALDH1L1 is markedly down regulated at the protein level. As the cell cultures become confluent and cells experience increased contact inhibition, ALDH1L1 accumulates in the cells. In agreement with this finding, NIH3T3 cells arrested in G₁/S-phase by a thymidine block completely lose the ALDH1L1 protein. Treatment with the proteasome inhibitor MG-132 prevents such loss in proliferating NIH3T3 cells, suggesting the proteasomal degradation of the ALDH1L1 protein. The co-localization of ALDH1L1 with proteasomes, demonstrated by confocal microscopy, supports this mechanism. We further show that ALDH1L1 interacts with the chaperone-dependent E3 ligase CHIP, which plays a key role in the ALDH1L1 ubiquitination and degradation. In NIH3T3 cells, silencing of CHIP by siRNA halts, while transient expression of CHIP promotes, the ALDH1L1 loss. The downregulation of ALDH1L1 is associated with the accumulation of the ALDH1L1 substrate 10-formyltetrahydrofolate, which is required for de novo purine biosynthesis, a key pathway activated in S-phase. Overall, our data indicate that CHIP-mediated proteasomal degradation of ALDH1L1 facilitates cellular proliferation.

Deep Sequencing Revealed a CpG Methylation Pattern Associated With ALDH1L1 Suppression in Breast Cancer

Hypermethylation of promoter CpG islands is generally recognized epigenetic mechanism responsible for gene silencing in cancer. However, molecular details on how this epigenetic mark triggers the process of gene downregulation are still elusive. Here, we used deep bisulfite sequencing and qPCR analysis to investigate the pattern of CpG methylation of ALDH1L1 promoter region and its association with the gene expression level in 16 paired breast cancer (BC) samples of different clinical stages. Expression of ALDH1L1 gene was suppressed in all examined BC samples up to 200-fold, and average hypermethylation level of the promoter region correlated positively with ALDH1L1 downregulation. We determined the role of every individual CpG site within the ALDH1L1 promoter, including upstream untranscribed region, first untranslated exon, and the start of the first intron, in aberrant gene expression by correlation analysis. The search revealed CpG sites which methylation has the highest impact on intensity of gene transcription. The majority of such CpG sites are located in a compact region in the first intron of the ALDH1L1 gene. These results assist in unraveling of dynamic nature of CpG promoter hypermethylation as well as demonstrate the efficiency of deep bisulfite sequencing in search for novel epigenetic markers in cancer.

Efficacy of epidermal growth factor receptor-tyrosine kinase inhibitors for lung squamous carcinomas harboring EGFR mutation: A multicenter study and pooled analysis of published reports

Epidermal growth factor receptor (EGFR) mutations are common in lung adenocarcinoma (ADC) but rare in squamous cell carcinoma (SQC). The efficacy of EGFR-tyrosine kinase inhibitors (EGFR-TKIs) for SQC with EGFR mutations is unclear. The aim of this study was to evaluate the efficacy of EGFR-TKIs for these patients. We performed a retrospective matched-pair case-control study from 3 cancer centers, including 44 SQC and 44 ADC patients with EGFR mutation who were treated with EGFR-TKI. Subsequently, we performed a pooled analysis on the efficacy of EGFR-TKIs for EGFR-mutant SQC in 115 patients, including 71 patients selected from 25 published reports. In our multicenter study, EGFR-mutant SQC and ADC patients had similar objective response rate (ORR) (43.2% vs. 54.5%, p = 0.290), but SQC patients had lower disease control rate (DCR) (71.3% vs. 100%, p = 0.001), significant shorter median progression free survival (PFS) (5.1 vs. 13.0 months, p = 0.000) and median overall survival (OS) (17.2 vs. 23.6 months, p = 0.027). In pooled analysis, the ORR, DCR, PFS and OS of SQC patients were 39.1%, 71.3%, 5.6 months and 15.0 months, respectively. Performance status was the only independent predictor of PFS and erlotinib treatment was associated with a better survival. In conclusion, EGFR-TKI was less effective in EGFR-mutant SQC than in ADC but still has clinical benefit for SQC patients. Further study is need to evaluate the using of EGFR-TKIs in these SQC patients.

ALDH1L1 and ALDH1L2 Folate Regulatory Enzymes in Cancer

Epidemiological studies implicate excess ethanol ingestion as a risk factor for several cancers and support the concept of a synergistic effect of chronic alcohol consumption and folate deficiency on carcinogenesis. Alcohol consumption affects folate-related genes and enzymes including two major folate-metabolizing enzymes, ALDH1L1 and ALDH1L2. ALDH1L1 (cytosolic 10-formyltetrahydrofolate dehydrogenase) is a regulatory enzyme in folate metabolism that controls the overall flux of one-carbon groups in folate-dependent biosynthetic pathways. It is strongly and ubiquitously down-regulated in malignant tumors via promoter methylation, and recent studies underscored this enzyme as a candidate tumor suppressor and potential marker of aggressive cancers. A related enzyme, ALDH1L2, is the mitochondrial homolog of ALDH1L1 encoded by a separate gene. In contrast to its cytosolic counterpart, ALDH1L2 is expressed in malignant tumors and cancer cell lines and was implicated in metastasis regulation. This review discusses the link between folate and cancer, modifying effects of alcohol consumption on folate-associated carcinogenesis, and putative roles of ALDH1L1 and ALDH1L2 in this process.

Comprehensive epigenetic analysis of the signature genes in lung adenocarcinoma

Aim: This study aimed to explore the epigenetic modifications of signature genes in lung adenocarcinoma. Materials & methods: The data of miRNA expression, mRNA expression and DNA methylation were downloaded from The Cancer Genome Atlas. Differential analysis was performed, followed by correlation analysis of miRNA–mRNA and DNA methylation-mRNA. Results: A total of 14 significant inverse correlations between gene expression and DNA methylation were identified, the expressions of which were selected for further validation via GSE27262, displaying similar pattern with that of the integrated analysis. In addition, qRT-PCR results showed that the expression profiling results of six mRNAs and one miRNA were consistent with the findings of integrated analysis. Five genes showed higher diagnostic value, which was also associated with overall survival of patients. Conclusion: Taken together, the epigenetic alterations of signature genes may hold promise for becoming biomarkers for the early detection of lung adenocarcinoma.

ALDH1L1 variant rs2276724 and mRNA expression predict post-operative clinical outcomes and are associated with TP53 expression in HBV-related hepatocellular carcinoma

Aldehyde dehydrogenase 1 family member L1 (ALDH1L1) is downregulated in hepatocellular carcinoma (HCC) tumors, and its decreased expression is associated with the poor prognosis of HCC patients. We, therefore, evaluated the effect of single nucleotide polymorphisms (SNPs) of ALDH1L1, and its mRNA expression on the survival of hepatitis B virus (HBV)-related HCC patients and the association with tumor protein p53 (TP53) expression. ALDH1L1 SNPs in 415 HBV-related HCC patients were genotyped via direct sequencing. Expression profile chip datasets and survival information were obtained from GSE14520. The C allele (CT/CC) carriers of rs2276724 were significantly associated with a favorable prognosis [adjusted P=0.040; adjusted hazard ratio (HR)=0.725; 95% confidence interval (CI)=0.533–0.986]. Joint-effect analyses suggested that the CT/CC genotype of rs2276724 in TP53-negative patients was significantly associated with a decreased risk of death, compared to the TT genotype of rs2276724 in TP53-positive patients (adjusted P=0.037; adjusted HR=0.621; 95% CI=0.396–0.973). Furthermore, low expression of ALDH1L1 predicted a poor prognosis for the HBV-related HCC patients (adjusted P=0.04 for disease-free survival; adjusted P=0.001 for overall survival). Patients with high ALDH1L1 expression and low TP53 expression were significantly associated with a decreased risk of recurrence and death, and patients with a high TP53 expression were also significantly associated with a decreased risk of death in HBV-related HCC, compared with low ALDH1L1 and low TP53 expression. Our results suggest that ALDH1L1 may be a biomarker for predicting postoperative clinical outcomes. Moreover, ALDH1L1-rs2276724 and mRNA expression were associated with TP53 expression in HBV-related HCC patients.

Obesity and menopause modify the epigenomic profile of breast cancer

Obesity is a high risk factor for breast cancer. This relationship could be marked by a specific methylome. The current work was aimed to explore the impact of obesity and menopausal status on variation in breast cancer methylomes. Data from Infinium 450K array-based methylomes of 64 breast tumors were coupled with information on BMI and menopausal status. Additionally, DNA methylation results were validated in 18 non-tumor and 81 tumor breast samples. Breast tumors arising in either pre- or postmenopausal women stratified by BMI or menopausal status alone were not associated with a specific DNA methylation pattern. Intriguingly, the DNA methylation pattern identified in association with the high-risk group (postmenopausal women with high BMI (>25) and premenopausal women with normal or low BMI < 25) exclusively characterized by hypermethylation of 1287 CpG sites as compared with the low-risk group. These CpG sites included the promoter region of fourteen protein-coding genes of which CpG methylation over the ZNF577 promoter region represents the top scoring associated event. In an independent cohort, the ZNF577 promoter methylation remained statistically significant in association with the high-risk group. Additionally, the impact of ZNF577 promoter methylation on mRNA expression levels was demonstrated in breast cancer cell lines after treatment with a demethylating agent (5-azacytidine). In conclusion, the epigenome of breast tumors is affected by a complex interaction between BMI and menopausal status. The ZNF577 methylation quantification is clearly relevant for the development of novel biomarkers of precision therapy in breast cancer.

Downregulation of RPL15 may predict poor survival and associate with tumor progression in pancreatic ductal adenocarcinoma

Early diagnosis and treatment in pancreatic ductal adenocarcinoma (PDAC) is still a challenge worldwide. The poor survival of PDAC patients mainly due to early metastasis when first diagnosed and lack of prognostic biomarker. Ribosomal protein L15 (RPL15), an RNA-binding protein, is a component of ribosomal 60S subunit. It was reported that RPL15 is dysregulated in various type of cancers. However, little is known about the role of RPL15 in PDAC carcinogenesis and progression. Herein, we clarified RPL15 expression status may serve as an independent prognostic biomarker in three independent PDAC patient cohorts. We found that RPL15 was dramatically decreased in PDAC tissues and cell lines. The high expression of RPL15 was inversely correlated with TNM stage, histological differentiation, T classification and vascular invasion. Low expression of RPL15 was significantly associated with poor overall survival of PDAC patients. Furthermore, we demonstrated that the reduction of RPL15 may promote invasion ability of pancreatic cell by inducing EMT process. In conclusion, decreased RPL15 expression is associated with invasiveness of PDAC cells, and RPL15 expression status may serve as a reliable prognostic biomarker in PDAC patients.

Overexpression of CD44 is associated with the occurrence and migration of non-small cell lung cancer

Non-small cell lung cancer (NSCLC) is a potentially fatal disease and the incidence is increasing annually. In order to diagnose and treat NSCLC effectively, greater understanding of its molecular mechanism is required. In the present study, 36 NSCLC tissues and 10 normal tissues were selected. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was used to analyze the CD44 mRNA expression level in NSCLC tissue and DNA sequencing was performed to further verify the CD44 expression level. Differentially expressed genes between tumor tissues and controls were determined by DNA sequencing and the Gene_act_net between CD44 and its associated genes was constructed. Gene Ontology (GO) term enrichment analysis of the differentially expressed genes was performed by the Biological Networks Gene Ontology tool. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed based on the Expression Analysis Systematic Explorer test applied in the Database for Annotation, Visualization and Integrated Discovery. RT-qPCR results showed that CD34 was overexpressed in 21 of the 36 NSCLC tissues (58.3%). The Gene_act_net indicated that there were 20 differentially expressed genes with 17 upregulated and 3 downregulated. Among them, CD44, MET, ERBB2, EGFR, AKT1, IQGAP1 and STAT3 were associated with the occurrence and migration of NSCLC. In KEGG pathway analysis, extracellular matrix-receptor interaction and hematopoietic cell lineage pathways were the most affected by overexpressed CD44; and thus may be important in the development and migration of NSCLC. In conclusion, CD44 was overexpressed in NSCLC and the overexpression was associated with the occurrence of NSCLC and migration of NSCLC cells.

CrossHub: a tool for multi-way analysis of The Cancer Genome Atlas (TCGA) in the context of gene expression regulation mechanisms

The contribution of different mechanisms to the regulation of gene expression varies for different tissues and tumors. Complementation of predicted mRNA–miRNA and gene–transcription factor (TF) relationships with the results of expression correlation analyses derived for specific tumor types outlines the interactions with functional impact in the current biomaterial. We developed CrossHub software, which enables two-way identification of most possible TF–gene interactions: on the basis of ENCODE ChIP-Seq binding evidence or Jaspar prediction and co-expression according to the data of The Cancer Genome Atlas (TCGA) project, the largest cancer omics resource. Similarly, CrossHub identifies mRNA–miRNA pairs with predicted or validated binding sites (TargetScan, mirSVR, PicTar, DIANA microT, miRTarBase) and strong negative expression correlations. We observed partial consistency between ChIP-Seq or miRNA target predictions and gene–TF/miRNA co-expression, demonstrating a link between these indicators. Additionally, CrossHub expression-methylation correlation analysis can be used to identify hypermethylated CpG sites or regions with the greatest potential impact on gene expression. Thus, CrossHub is capable of outlining molecular portraits of a specific gene and determining the three most common sources of expression regulation: promoter/enhancer methylation, miRNA interference and TF-mediated activation or repression. CrossHub generates formatted Excel workbooks with the detailed results. CrossHub is freely available at https://sourceforge.net/projects/crosshub/.

Expression and DNA methylation alterations of seven cancer-associated 3p genes and their predicted regulator miRNAs (miR-129-2, miR-9-1) in breast and ovarian cancers

The methylation of promoter CpG islands and interactions between microRNAs (miRNAs) and messenger RNAs (mRNAs) of target genes are considered two crucial epigenetic mechanisms for inducing gene and pathway deregulation in tumors. Here, the expression levels of seven cancer-associated 3p genes (RASSF1(isoform A), RARB(isoform 2), SEMA3B, RHOA, GPX1, NKIRAS1, and CHL1) and their predicted regulator miRNAs (miR-129-2, miR-9-1) were analyzed in breast (BC, 40 samples) and ovarian (OC, 14 samples) cancers using RT-PCR and qPCR. We first revealed a negative correlation between the level of the miR-129-2 precursor and RASSF1(A) and GPX1 mRNA levels in BC (Spearman's correlation coefficient (rs) was − 0.26 in both cases). Similar results were observed for the miR-129-2 precursor and the RASSF1(A), GPX1, RARB(2), and CHL1 genes in OC (rs was in the range − 0.48 to − 0.54). Using methylation-specific PCR, a significant correlation was shown between promoter hypermethylation and the down-regulation of the RASSF1(A), GPX1, RARB(2), SEMA3B, MIR-129-2, and MIR-9-1 genes in BC (rs = 0.41 to 0.75) and of the RASSF1(A) gene in OC (rs = 0.67). We first demonstrated a high hypermethylation frequency of MIR-129-2 and SEMA3B (up to 45 to 48%) in both BC (69 samples) and OC (41 samples). Moreover, we observed a positive correlation between the hypermethylation of MIR-129-2 and the up-regulation of the RASSF1(A) and GPX1 genes in BC (rs = 0.38 and 0.42, respectively). QPCR analysis of the expression of RASSF1(A) and mature miR-129-2 in additional BC sample set (24 samples) revealed a negative correlation between them (rs = − 0.41) that strengthened the results obtained during the analysis of miR-129-2 precursor level. In summary, the obtained data indicate the involvement of methylation in the down-regulation of the studied coding and miRNA genes and suggest the involvement of miR-129-2 in the deregulation of RASSF1(A) via a direct interaction or/and mediators in common pathways (according to KEGG, Gene Ontology (FDR < 0.01), and GeneCards data) in the examined gynecological tumors.

Integrin α9 gene promoter is hypermethylated and downregulated in nasopharyngeal carcinoma

Epigenetic silencing of tumor suppressor genes (TSGs) by promoter methylation can be an early event in the multi-step process of carcinogenesis. Human chromosome 3 contains clusters of TSGs involved in many cancer types including nasopharyngeal carcinoma (NPC), the most common cancer in Southern China. Among ten candidate TSGs identified in chromosome 3 using NotI microarray, ITGA9 and WNT7A could be validated. 5'-aza-2' deoxycytidine treatment restored the expression of ITGA9 and WNT7A in two NPC cell lines. Immunostaining showed strong expression of these genes in the membrane and cytoplasm of adjacent control nasopharyngeal epithelium cells, while they were weakly expressed in NPC tumor cells. The ITGA9 promoter showed marked differentially methylation between tumor and control tissue, whereas no differentially methylation could be detected for the WNT7A promoter. The expression level of ITGA9 in NPC tumors was downregulated 4.9-fold, compared to the expression in control. ITGA9 methylation was detected by methylation specific PCR (MSP) in 56% of EBV positive NPC-cases with 100% specificity. Taken together, this suggests that ITGA9 might be a TSG in NPC that is involved in tumor cell biology. The possibility of using ITGA9 methylation as a marker for early detection of NPC should further be explored.

Identification of Novel Epigenetic Markers of Prostate Cancer by NotI-Microarray Analysis

A significant need for reliable and accurate cancer diagnostics and prognosis compels the search for novel biomarkers that would be able to discriminate between indolent and aggressive tumors at the early stages of disease. The aim of this work was identification of potential diagnostic biomarkers for characterization of different types of prostate tumors. NotI-microarrays with 180 clones associated with chromosome 3 genes/loci were applied to determine genetic and epigenetic alterations in 33 prostate tumors. For 88 clones, aberrations were detected in more than 10% of tumors. The major types of alterations were DNA methylation and/or deletions. Frequent methylation of the discovered loci was confirmed by bisulfite sequencing on selective sampling of genes: FGF12, GATA2, and LMCD1. Three genes (BHLHE40, BCL6, and ITGA9) were tested for expression level alterations using qPCR, and downregulation associated with hypermethylation was shown in the majority of tumors. Based on these data, we proposed the set of potential biomarkers for detection of prostate cancer and discrimination between prostate tumors with different malignancy and aggressiveness: BHLHE40, FOXP1, LOC285205, ITGA9, CTDSPL, FGF12, LOC440944/SETD5, VHL, CLCN2, OSBPL10/ZNF860, LMCD1, FAM19A4, CAND2, MAP4, KY, and LRRC58. Moreover, we probabilistically estimated putative functional relations between the genes within each set using the network enrichment analysis.

Membrane carbonic anhydrase IX expression and relapse risk in resected stage I-II non-small-cell lung cancer

BACKGROUND

Adjuvant chemotherapy reduces recurrences of non-small-cell lung cancer (NSCLC). To determine which patients need adjuvant chemotherapy, we assessed factors associated with time to relapse (TTR).

METHODS

In 230 resected stage I-II NSCLCs, we correlated immunohistochemistry scores for factors associated with cell growth rate, growth regulation, hypoxia, cell survival, and cell death with TTR.

RESULTS

With a median follow-up of 82 months (1-158) for those alive and relapse free at last follow-up, median time to recurrence was not reached. The 2- and 5-year probabilities of maintaining freedom from recurrence were 80.7% (95% confidence interval, 75.3%, 86.4%) and 74.6% (95% confidence interval, 68.6%, 81.2%), respectively. TTR curves flattened at an apparent cure rate of 70%. In multicovariate Cox models, factors correlating with shorter TTR were membranous carbonic anhydrase IX (mCAIX) staining (any versus none, hazard ratio = 2.083, p = 0.023) and node stage (N1 versus N0, hazard ratio = 2.591, p = 0.002). mCAIX scores correlated positively with tumor size, grade, squamous histology, necrosis, mitoses, Ki67, p53, nuclear DNA methyltransferase 1, and cytoplasmic enhancer-of-split-and-hairy-related protein, and they correlated inversely with papillary histology, epidermal growth factor receptor mutation (trend), copper transporter-1, and cytoplasmic hypoxia-inducible factor-1α, vascular endothelial growth factor, DNA methyltransferase 1, and excision repair cross-complementing rodent repair deficiency, complementation group 1.

CONCLUSION

Nodal stage and mCAIX immunohistochemistry were the strongest independent predictors of shorter TTR in resected NSCLCs. mCAIX correlated with tumor size, markers of tumor proliferation and necrosis, and tumor genetic characteristics, and it paradoxically correlated inversely with the hypoxia markers, hypoxia-inducible factor-1α and vascular endothelial growth factor. Presence of mCAIX could help determine patients with high risk of recurrence who might require adjuvant chemotherapy.

5-azacytidine enhances efficacy of multiple chemotherapy drugs in AML and lung cancer with modulation of CpG methylation

The DNA methyltransferase (DNMT) inhibitory drugs such as 5-azacytidine induce DNA hypomethylation by inhibiting DNA methyltransferases. While clinically effective, DNMT inhibitors are not curative. A combination with cytotoxic drugs might be beneficial, but this is largely unexplored. In the present study, we analyzed potential synergisms between cytotoxic drugs and 5-azacytidine in acute myeloid leukemia (AML) and non-small cell lung cancer (NSCLC) cells. Lung cancer and leukemia cell lines were exposed to low doses of 5-azacytidine with varying doses of cytarabine or etoposide for AML cells (U937 and HL60) as well as cisplatin or gemcitabine for NSCLC cells (A549 and HTB56) for 48 h. Drug interaction and potential synergism was analyzed according to the Chou-Talalay algorithm. Further analyses were based on soft agar colony formation assays, active caspase-3 staining and BrdU incorporation flow cytometry. To identify effects on DNA methylation patterns, we performed genome wide DNA methylation analysis using 450K bead arrays. Azacytidine at low doses was synergistic with cytotoxic drugs in NSCLC and in AML cell lines. Simultaneous exposure to 5-azacytidine with cytotoxic drugs showed strong synergistic activity. In colony formation assays these synergisms were repeatedly verified for 5-azacytidine (25 nM) with low doses of anticancer agents. 5-azacytidine neither affected the cell cycle nor increased apoptosis. 450K methylation bead arrays revealed 1,046 CpG sites in AML and 1,778 CpG sites in NSCLC cells with significant DNA hypomethylation (24-h exposure) to 5-azacytidine combined with the cytotoxic drugs. These CpG-sites were observed in the candidate tumor-suppressor genes MGMT and THRB. Additional incubation time after 24-h treatment led to a 4.1-fold increase of significant hypomethylated CpG-sites in NSCLC cells. These results suggest that the addition of DNA demethylating agents to cytotoxic anticancer drugs exhibits synergistic activity in AML and NSCLC. Dysregulation of an equilibrium of DNA methylation in cancer cells might increase the susceptibility for cytotoxic drugs.

A bayesian integrative model for genetical genomics with spatially informed variable selection

We consider a Bayesian hierarchical model for the integration of gene expression levels with comparative genomic hybridization (CGH) array measurements collected on the same subjects. The approach defines a measurement error model that relates the gene expression levels to latent copy number states. In turn, the latent states are related to the observed surrogate CGH measurements via a hidden Markov model. The model further incorporates variable selection with a spatial prior based on a probit link that exploits dependencies across adjacent DNA segments. Posterior inference is carried out via Markov chain Monte Carlo stochastic search techniques. We study the performance of the model in simulations and show better results than those achieved with recently proposed alternative priors. We also show an application to data from a genomic study on lung squamous cell carcinoma, where we identify potential candidates of associations between copy number variants and the transcriptional activity of target genes. Gene ontology (GO) analyses of our findings reveal enrichments in genes that code for proteins involved in cancer. Our model also identifies a number of potential candidate biomarkers for further experimental validation.

Gestational hypoxia and epigenetic programming of brain development disorders

Adverse environmental conditions faced by an individual early during its life, such as gestational hypoxia, can have a profound influence on the risk of diseases, such as neurological disorders, in later life. Clinical and preclinical studies suggest that epigenetic programming of gene expression patterns in response to maternal stress have a crucial role in the fetal origins of neurological diseases. Herein, we summarize recent studies regarding the role of epigenetic mechanisms in the developmental programming of neurological diseases in offspring, primarily focusing on DNA methylation/demethylation and miRNAs. Such information could increase our understanding of the fetal origins of adult diseases and help develop effective prevention and intervention against neurological diseases.

Development of a multiplex methylation specific PCR suitable for (early) detection of non-small cell lung cancer

Lung cancer is a worldwide health problem and a leading cause of cancer-related deaths. Silencing of potential tumor suppressor genes (TSGs) by aberrant promoter methylation is an early event in the initiation and development of cancer. Thus, methylated cancer type-specific TSGs in DNA can serve as useful biomarkers for early cancer detection. We have now developed a "Multiplex Methylation Specific PCR" (MMSP) assay for analysis of the methylation status of multiple potential TSGs by a single PCR reaction. This method will be useful for early diagnosis and treatment outcome studies of non-small cell lung cancer (NSCLC). Genome-wide CpG methylation and expression microarrays were performed on lung cancer tissues and matched distant non-cancerous tissues from three NSCLC patients from China. Thirty-eight potential TSGs were selected and analyzed by methylation PCR on bisulfite treated DNA. On the basis of sensitivity and specificity, six marker genes, HOXA9, TBX5, PITX2, CALCA, RASSF1A, and DLEC1, were selected to establish the MMSP assay. This assay was then used to analyze lung cancer tissues and matched distant non-cancerous tissues from 70 patients with NSCLC, as well as 24 patients with benign pulmonary lesion as controls. The sensitivity of the assay was 99% (69/70). HOXA9 and TBX5 were the 2 most sensitive marker genes: 87% (61/70) and 84% (59/70), respectively. RASSF1A and DLEC1 showed the highest specificity at 99% (69/70). Using the criterion of identifying at least any two methylated marker genes, 61/70 cancer samples were positive, corresponding to a sensitivity of 87% and a specificity of 94%. Early stage I or II NSCLC could even be detected with a 100% specificity and 86% sensitivity. In conclusion, MMSP has the potential to be developed into a population-based screening tool and can be useful for early diagnosis of NSCLC. It might also be suitable for monitoring treatment outcome and recurrence.

Epigenetic alterations of chromosome 3 revealed by NotI-microarrays in clear cell renal cell carcinoma

This study aimed to clarify epigenetic and genetic alterations that occur during renal carcinogenesis. The original method includes chromosome 3 specific NotI-microarrays containing 180 NotI-clones associated with 188 genes for hybridization with 23 paired normal/tumor DNA samples of primary clear cell renal cell carcinomas (ccRCC). Twenty-two genes showed methylation and/or deletion in 17–57% of tumors. These genes include tumor suppressors or candidates (VHL, CTDSPL, LRRC3B, ALDH1L1, and EPHB1) and genes that were not previously considered as cancer-associated (e.g., LRRN1, GORASP1, FGD5, and PLCL2). Bisulfite sequencing analysis confirmed methylation as a frequent event in ccRCC. A set of six markers (NKIRAS1/RPL15, LRRN1, LRRC3B, CTDSPL, GORASP1/TTC21A, and VHL) was suggested for ccRCC detection in renal biopsies. The mRNA level decrease was shown for 6 NotI-associated genes in ccRCC using quantitative PCR: LRRN1, GORASP1, FOXP1, FGD5, PLCL2, and ALDH1L1. The majority of examined genes showed distinct expression profiles in ccRCC and papillary RCC. The strongest extent and frequency of downregulation were shown for ALDH1L1 gene both in ccRCC and papillary RCC. Moreover, the extent of ALDH1L1 mRNA level decrease was more pronounced in both histological types of RCC stage III compared with stages I and II (P = 0.03). The same was observed for FGD5 gene in ccRCC (P < 0.06).

Dedicated to thememory of Eugene R. Zabarovsky

Identification and comparison of aberrant key regulatory networks in breast, colon, liver, lung, and stomach cancers through methylome database analysis

Aberrant methylation of specific CpG sites at the promoter is widely responsible for genesis and development of various cancer types. Even though the microarray-based methylome analyzing techniques have contributed to the elucidation of the methylation change at the genome-wide level, the identification of key methylation markers or top regulatory networks appearing common in highly incident cancers through comparison analysis is still limited. In this study, we in silico performed the genome-wide methylation analysis on each 10 sets of normal and cancer pairs of five tissues: breast, colon, liver, lung, and stomach. The methylation array covers 27,578 CpG sites, corresponding to 14,495 genes, and significantly hypermethylated or hypomethylated genes in the cancer were collected (FDR adjusted p-value <0.05; methylation difference >0.3). Analysis of the dataset confirmed the methylation of previously known methylation markers and further identified novel methylation markers, such as GPX2, CLDN15, and KL. Cluster analysis using the methylome dataset resulted in a diagram with a bipartite mode distinguishing cancer cells from normal cells regardless of tissue types. The analysis further revealed that breast cancer was closest with lung cancer, whereas it was farthest from colon cancer. Pathway analysis identified that either the “cancer” related network or the “cancer” related bio-function appeared as the highest confidence in all the five cancers, whereas each cancer type represents its tissue-specific gene sets. Our results contribute toward understanding the essential abnormal epigenetic pathways involved in carcinogenesis. Further, the novel methylation markers could be applied to establish markers for cancer prognosis.

BMP-4 genetic variants and protein expression are associated with platinum-based chemotherapy response and prognosis in NSCLC

To explore the role of genetic polymorphisms of bone morphogenic proteins 4 (BMP-4) in the response to platinum-based chemotherapy and the clinical outcome in patients with advanced nonsmall cell lung cancer (NSCLC), 938 patients with stage III (A+B) or IV NSCLC were enrolled in this study. We found that the variant genotypes of 6007C > T polymorphisms significantly associated with the chemotherapy response. The 6007CC genotype carriers had a higher chance to be responder to chemotherapy (adjusted odd ratio = 2.77; 95% CI: 1.83–4.18; adjusted < 0.001). The 6007C > T polymorphisms and BMP-4 expression also affect the prognosis of NSCLC. Patients with high BMP-4 expression had a significantly higher chance to be resistant to chemotherapy than those with low BMP-4 expression (OR = 2.81; 95% CI: 1.23–6.44; P = 0.01). The hazard ratio (HR) for 6007TT was 2.37 times higher than 6007CC (P = 0.003). In summary, the 6007C > T polymorphism of BMP-4 gene and BMP-4 tissue expression may be used as potential predictor for the chemotherapy response and prognosis of advanced NSCLC.

Classification of non-small cell lung cancer based on copy number alterations

Lung cancer is one of the leading causes of cancer mortality worldwide and non–small cell lung cancer (NSCLC) accounts for the most part. NSCLC can be further divided into adenocarcinoma (ACA) and squamous cell carcinoma (SCC). It is of great value to distinguish these two subgroups clinically. In this study, we compared the genome-wide copy number alterations (CNAs) patterns of 208 early stage ACA and 93 early stage SCC tumor samples. As a result, 266 CNA probes stood out for better discrimination of ACA and SCC. It was revealed that the genes corresponding to these 266 probes were enriched in lung cancer related pathways and enriched in the chromosome regions where CNA usually occur in lung cancer. This study sheds lights on the CNA study of NSCLC and provides some insights on the epigenetic of NSCLC.

Suppression of GNAI2 message in ovarian cancer

Background

Understanding the integration of hormone signaling and how it impacts oncogenesis is critical for improved cancer treatments. Here we elucidate GNAI2 message alterations in ovarian cancer (OvCa). GNAI2 is a heterotrimeric G protein which couples cell surface hormone receptors to intracellular enzymes, and is best characterized for its direct role in regulating cAMP response element-binding protein (CREB) function by decreasing intracellular cAMP through inhibiting adenylyl cyclase.

Methods

We probed the Origene human OvCa array for the presence of polymorphisms and gene expression alterations of GNAI2 using directing sequencing and qPCR. These data were supported by database mining of the [NCBI NIH GSE:6008, GSE:14764, GSE:29450, GDS:4066, GDS:3297, GSE:32474, and GSE:2003] datasets.

Results

No significant polymorphisms were found, including an absence of the gip2 oncogene. However, 85.9% of (506 of 589) OvCa patients had decreased GNAI2 message. Further characterization demonstrated that the GNAI2 message was on average decreased 54% and maximally decreased by 2.8 fold in clear cell carcinoma. GNAI2 message decreased in early stage cancer while message was increased compared to normal in advanced cancers. The changes in GNAI2 also correlated to deregulation of CREB, Fos, Myc, cyclins, Arf, the transition from estrogen dependence to independence, and metastatic potential.

Conclusion

These data strongly implicate GNAI2 as a critical regulator of oncogenesis and an upstream driver of cancer progression in OvCa.