POLD2 and KSP37 (FGFBP2) correlate strongly with histology, stage and outcome in ovarian carcinomas

scRNA-seq revealed high stemness epithelial malignant cell clusters and prognostic models of lung adenocarcinoma

Lung adenocarcinoma (LUAD) is one of the sole causes of death in lung cancer patients. This study combined with single-cell RNA-seq analysis to identify tumor stem-related prognostic models to predict the prognosis of lung adenocarcinoma, chemotherapy agents, and immunotherapy efficacy. mRNA expression-based stemness index (mRNAsi) was determined by One Class Linear Regression (OCLR). Differentially expressed genes (DEGs) were detected by limma package. Single-cell RNA-seq analysis in GSE123902 dataset was performed using Seurat package. Weighted Co-Expression Network Analysis (WGCNA) was built by rms package. Cell differentiation ability was determined by CytoTRACE. Cell communication analysis was performed by CellCall and CellChat package. Prognosis model was constructed by 10 machine learning and 101 combinations. Drug predictive analysis was conducted by pRRophetic package. Immune microenvironment landscape was determined by ESTIMATE, MCP-Counter, ssGSEA analysis. Tumor samples have higher mRNAsi, and the high mRNAsi group presents a worse prognosis. Turquoise module was highly correlated with mRNAsi in TCGA-LUAD dataset. scRNA analysis showed that 22 epithelial cell clusters were obtained, and higher CSCs malignant epithelial cells have more complex cellular communication with other cells and presented dedifferentiation phenomenon. Cellular senescence and Hippo signaling pathway are the major difference pathways between high- and low CSCs malignant epithelial cells. The pseudo-temporal analysis shows that cluster1, 2, high CSC epithelial cells, are concentrated at the end of the differentiation trajectory. Finally, 13 genes were obtained by intersecting genes in turquoise module, Top200 genes in hdWGCNA, DEGs in high- and low- mRNAsi group as well as DEGs in tumor samples vs. normal group. Among 101 prognostic models, average c-index (0.71) was highest in CoxBoost + RSF model. The high-risk group samples had immunosuppressive status, higher tumor malignancy and low benefit from immunotherapy. This work found that malignant tumors and malignant epithelial cells have high CSC characteristics, and identified a model that could predict the prognosis, immune microenvironment, and immunotherapy of LUAD, based on CSC-related genes. These results provided reference value for the clinical diagnosis and treatment of LUAD.

Aging gene signature of memory CD8+ T cells is associated with neurocognitive functioning in Alzheimer's disease

BACKGROUND

Memory CD8+ T cells expand with age. We previously demonstrated an age-associated expansion of effector memory (EM) CD8+ T cells expressing low levels of IL-7 receptor alpha (IL-7Rαlow) and the presence of its gene signature (i.e., IL-7Rαlow aging genes) in peripheral blood of older adults without Alzheimer's disease (AD). Considering age as the strongest risk factor for AD and the recent finding of EM CD8+ T cell expansion, mostly IL-7Rαlow cells, in AD, we investigated whether subjects with AD have alterations in IL-7Rαlow aging gene signature, especially in relation to genes possibly associated with AD and disease severity.

RESULTS

We identified a set of 29 candidate genes (i.e., putative AD genes) which could be differentially expressed in peripheral blood of patients with AD through the systematic search of publicly available datasets. Of the 29 putative AD genes, 9 genes (31%) were IL-7Rαlow aging genes (P < 0.001), suggesting the possible implication of IL-7Rαlow aging genes in AD. These findings were validated by RT-qPCR analysis of 40 genes, including 29 putative AD genes, additional 9 top IL-7R⍺low aging but not the putative AD genes, and 2 inflammatory control genes in peripheral blood of cognitively normal persons (CN, 38 subjects) and patients with AD (40 mild cognitive impairment and 43 dementia subjects). The RT-qPCR results showed 8 differentially expressed genes between AD and CN groups; five (62.5%) of which were top IL-7Rαlow aging genes (FGFBP2, GZMH, NUAK1, PRSS23, TGFBR3) not previously reported to be altered in AD. Unbiased clustering analysis revealed 3 clusters of dementia patients with distinct expression levels of the 40 analyzed genes, including IL-7Rαlow aging genes, which were associated with neurocognitive function as determined by MoCA, CDRsob and neuropsychological testing.

CONCLUSIONS

We report differential expression of "normal" aging genes associated with IL-7Rαlow EM CD8+ T cells in peripheral blood of patients with AD, and the significance of such gene expression in clustering subjects with dementia due to AD into groups with different levels of cognitive functioning. These results provide a platform for studies investigating the possible implications of age-related immune changes, including those associated with CD8+ T cells, in AD.

Selinexor Synergistically Promotes the Antileukemia Activity of Venetoclax in Acute Myeloid Leukemia by Inhibiting Glycolytic Function and Downregulating the Expression of DNA Replication Genes

Introduction

The BCL-2 inhibitor venetoclax has been widely used in the treatment of acute myeloid leukemia (AML); however, AML patients treated with venetoclax gradually develop resistance. The exportin-1 (XPO1) inhibitor selinexor can synergistically promote the antileukemia activity of venetoclax, but the mechanism remains unclear.

Methods and Results

Annexin V/7-aminoactinomycin D assays were used to examine the effects of a combination of venetoclax and selinexor (VEN+SEL) on AML cell lines and primary AML cells. RNA sequencing and oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) determinations by a Seahorse XF analyzer were employed to investigate the molecular mechanism of the toxicity of the VEN+SEL combination to AML cells. The cytotoxicity of NK cell combined with VEN+SEL combination was assessed in vitro using flow cytometry. VEN+SEL enhanced the apoptosis of AML cells (KG-1A and THP-1) and primary AML samples in vitro. The ECAR and OCR results demonstrated that the VEN+SEL combination significantly inhibited glycolytic function. RNA sequencing of THP-1 cells demonstrated that DNA replication-related genes were downregulated after treatment with the VEN+SEL combination.

Conclusion

This study indicated that selinexor can synergistically enhance the antileukemia activity of venetoclax in AML cells in vitro by inhibiting glycolytic function and downregulating DNA replication-related genes. Based on our experimental data, combining selinexor with venetoclax is an appropriate advanced treatment option for AML patients.

Aging gene signature of IL-7 receptor alpha low effector memory CD8+ T cells is associated with neurocognitive functioning in Alzheimer's disease

CD45RA+ effector memory (EM) CD8+ T cell expansion was reported in Alzheimer's disease (AD). Such cells are IL-7 receptor alpha (IL-7Rα)low EM CD8+ T cells, which expand with age and have a unique aging gene signature (i.e., IL-7Rαlow aging genes). Here we investigated whether IL-7Rαlow aging genes and previously reported AD and memory (ADM) genes overlapped with clinical significance in AD patients. RT-qPCR analysis of 40 genes, including 29 ADM, 9 top IL-7Ralow aging and 2 control genes, showed 8 differentially expressed genes between AD and cognitively normal groups; five (62.5%) of which were top IL-7Rαlow aging genes. Over-representation analysis revealed that these genes were highly present in molecular and biological pathways associated with AD. Distinct expression levels of these genes were associated with neuropsychological testing performance in 3 subgroups of dementia participants. Our findings support the possible implication of the IL-7Rαlow aging gene signature with AD.

Early embryonic lethality of mice lacking POLD2

As a highly conserved DNA polymerase (Pol), Pol δ plays crucial roles in chromosomal DNA synthesis and various DNA repair pathways. However, the function of POLD2, the second small subunit of DNA Pol δ (p50 subunit), has not been characterized in vivo during mammalian development. Here, we report for the first time, the essential role of subunit POLD2 during early murine embryogenesis. Although Pold2 mutant mouse embryos exhibit normal morphology at E3.5 blastocyst stage, they cannot be recovered at gastrulation stages. Outgrowth assays reveal that mutant blastocysts cannot hatch from the zona pellucida, indicating impaired blastocyst function. Notably, these phenotypes can be recapitulated by small interfering RNA (siRNA)-mediated knockdown, which also exhibit slowed cellular proliferation together with skewed primitive endoderm and epiblast allocation during the second cell lineage specification. In summary, our study demonstrates that POLD2 is essential for the earliest steps of mammalian development, and the retarded proliferation and embryogenesis may also alter the following cell lineage specifications in the mouse blastocyst embryos.

Methylation-based reclassification and risk stratification of skull-base chordomas

BACKGROUND

Skull-base chordomas are rare malignant bone cancers originating from the remnant of the notochord. Survival is variable, and clinical or molecular factors cannot reliably predict their outcomes. This study therefore identified epigenetic subtypes that defined new chordoma epigenetic profiles and their corresponding characteristics.

METHODS

Methylation profiles of 46 chordoma-resected neoplasms between 2008 and 2014, along with clinical information, were collected. K-means consensus clustering and principal component analysis were used to identify and validate the clusters. Single-sample gene set enrichment analysis, methylCIBERSORT algorithm, and copy number analysis were used to identify the characteristics of the clusters.

RESULTS

Unsupervised clustering analysis confirmed two clusters with a progression-free survival difference. Gene set enrichment analysis indicated that the early and late estrogen response pathways and the hypoxia pathway were activated whereas the inflammatory and interferon gamma responses were suppressed. Forty-six potential therapeutic targets corresponding to differentially methylated sites were identified from chordoma patients. Subgroups with a worse outcome were characterized by low immune cell infiltration, higher tumor purity, and higher stemness indices. Moreover, copy number amplifications mostly occurred in cluster 1 tumors and the high-risk group. Additionally, the presence of a CCNE1 deletion was exclusively found in the group of chordoma patients with better outcome, whereas RB1 and CDKN2A/2B deletions were mainly found in the group of chordoma patients with worse outcome.

CONCLUSIONS

Chordoma prognostic epigenetic subtypes were identified, and their corresponding characteristics were found to be variable.

POLD2 is activated by E2F1 to promote triple-negative breast cancer proliferation

Triple-negative breast cancer (TNBC) is a highly malignant breast cancer subtype with a poor prognosis. Improved insight into the molecular biology basis of TNBC progression is urgently needed. Herein, we reported that POLD2 was highly expressed in TNBC and patients with high POLD2 expression in their tumors had poor clinical outcomes. In functional studies, knockdown of POLD2 inhibited the proliferation of TNBC. Mechanistically, we revealed that transcription factor E2F1 directly bound to the promoter of POLD2 and regulated its expression in TNBC cells, which in turn contributed to the proliferation of TNBC. Additionally, rescue experiments validated that E2F1-mediated cell proliferation in TNBC was dependent on POLD2. Taken together, our results elucidated a novel mechanism of the E2F1-POLD2 axis in TNBC proliferation, and POLD2 may be a potential therapeutic target for TNBC treatment.

An integrative analysis revealing POLD2 as a tumor suppressive immune protein and prognostic biomarker in pan-cancer

Introduction: POLD2 is an indispensable subunit of DNA polymerase δ, which is responsible for the synthesis of the backward accompanying strand in eukaryotic organisms. Current studies have found an association between POLD2 and the development of a variety of cancers. However, its value in cancer immunotherapy has not been fully established.

Methods: POLD2 expression was analyzed using RNA expression and clinical data from TCGA and GTEx databases. The prognostic impact of POLD2 on tumor patients was analyzed using clinical survival data from TCGA. Gene enrichment analysis was performed using the R package “cluster analyzer” to explore the role of POLD2. We used the TIMER2 database to analyze the relationship between immune cell infiltration and POLD2 expression in TCGA. We downloaded relevant data from TCGA and analyzed the relationship between POLD2 and immune checkpoints, immunosuppressive genes, immune activating genes, chemokines and chemokine receptors.

Results: POLD2 was significantly overexpressed in most tumors compared to normal tissue. High POLD2 expression was significantly associated with advanced tumor stage, significantly shorter overall survival and progression-free survival. Also, we found that POLD2 expression correlated strongly with immunomodulatory genes, and significantly negatively with most immune checkpoints (PD-L1, CTLA4, TIM3, and CD28). Pathway enrichment analysis suggests that low expression of POLD2 promotes immune regulation-related pathways and high expression promotes metabolic and DNA repair-related pathways. Furthermore, tumor microenvironment analysis suggests that high POLD2 expression inhibits infiltration of CD8⁺ T cells and CD4⁺ memory T cells.

Discussion: In conclusion, POLD2 may be a molecular biomarker for pan-cancer prognosis and immunotherapy. It may serve as a potential target for new insights in human tumor prognosis prediction and immunotherapy assessment.

The Mechanism and Prognostic Value of DNA Polymerase δ Subunits in Hepatocellular Carcinoma: Implications for Precision Therapy

PURPOSE

Polymerase δ (POLD) proteins is a pivotal B-family DNA polymerase in the process of genome replication and repair and are comprised of POLD1-4. The predictive value of POLDs in hepatocellular carcinoma (HCC) has not been evaluated until now.

PATIENTS AND METHODS

A total of 369 hepatocellular carcinoma samples and 50 adjacent normal samples were enrolled from the TCGA-LIHC database, and the GSE10186 database was also used. Transcription, methylation and genetic alteration status of HCC patients were evaluated by GEPIA, Kaplan-Meier plotter, cBioPortal, MethHC, MethSurv. SurvExpress was employed to generate the overall prognosis prediction signature of POLDs. POLDs coexpressed genes were explored and enriched in potential pathways. K-M curves were generated to compare the different survival results in different groups, while ROC curves were used to validate the efficiency of the POLD signature.

RESULTS

All four POLD subunits were highly expressed in HCC tumor tissues. POLD1-3 and increased mRNA levels were also positively associated with advanced tumor stage and OS prognosis. Methylation in the promoter of POLDs affects mRNA expression and OS, especially for some specific CpG sites. Meanwhile, POLDs could preferably predict the prognosis for patients who suffered from a high gene mutation burden. We evaluated the combined prognostic predictive value of four POLD subunits in both the TCGA-LIHC and GSE10186 databases and recognized the statistically significant HR of the high-risk group, along with the reliable predictive value. The coexpressed gene sets and annotation results showed that the POLD coexpressed genes were mostly associated with DNA repair and cell cycle regulation pathways.

CONCLUSION

POLD is an essential predictive factor for the prognosis of HCC. The united signature could precisely identify unfavorable clinical outcome of HCC.

IL12 integrated into the CAR exodomain converts CD8+ T cells to poly-functional NK-like cells with superior killing of antigen-loss tumors

Chimeric antigen receptor (CAR)-redirected T cell therapy often fails to control tumors in the long term due to selecting cancer cells that downregulated or lost CAR targeted antigen. To reprogram the functional capacities specifically of engineered CAR T cells, we inserted IL12 into the extracellular moiety of a CD28-ζ CAR; both the CAR endodomain and IL12 were functionally active, as indicated by antigen-redirected effector functions and STAT4 phosphorylation, respectively. The IL12-CAR reprogrammed CD8⁺ T cells toward a so far not recognized natural killer (NK) cell-like signature and a CD94⁺CD56⁺CD62L^high phenotype closely similar, but not identical, to NK and cytokine induced killer (CIK) cells. In contrast to conventional CAR T cells, IL12-CAR T cells acquired antigen-independent, human leukocyte antigen E (HLA-E) restricted cytotoxic capacities eliminating antigen-negative cancer cells in addition to eliminating cancer cells with CAR cognate antigen. Simultaneous signaling through both the CAR endodomain and IL12 were required for inducing maximal NK-like cytotoxicity; adding IL12 to conventional CAR T cells was not sufficient. Antigen-negative tumors were attacked by IL12-CAR T cells, but not by conventional CAR T cells. Overall, we present a prototype of a new family of CARs that augments tumor recognition and elimination through expanded functional capacities by an appropriate cytokine integrated into the CAR exodomain.

Identification of cuproptosis-related lncRNA prognostic signature for osteosarcoma

BACKGROUND

Copper is an indispensably mineral element involved in various metabolic processes and functions in the active sites of many metalloproteins. Copper dysregulation is associated with cancers such as osteosarcoma (OS), the most common primary bone malignancy with invasiveness and metastasis. However, the causality between cuproptosis and OS remains elusive. We aim to identify cuproptosis-related long non-coding RNAs (lncRNAs) for osteosarcomatous prognosis, immune microenvironment response, and immunotherapy.

METHODS

The Person correlation and differential expression analysis were used to identify differentially expressed cuproptosis-related lncRNAs (CRLs). The univariate, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analysis were performed to construct the CRL signature. The Kaplan-Meier (K-M) survival analysis, receiver operating characteristic (ROC) curve, internal validation, independent prognostic analysis, and nomograph were used to evaluate the prognostic value. The functional enrichment, tumor microenvironment, immunotherapy and chemotherapy response between the two distinct groups were further explored using a series of algorithms. The expression of signature CRLs was verified by real-time quantitative polymerase chain reaction (RT-qPCR) in OS cell lines.

RESULTS

A novel CRL signature consisting of four CRLs were successfully identified. The K-M survival analysis indicated that the OS patients in the low-risk groups had a better prognosis than that in the high-risk group. Then, the ROC curve and subgroup survival analysis confirmed the prognostic evaluation performance of the signature. Equally, the independent prognostic analysis demonstrated that the CRL signature was an independently predicted factor for OS. Friends analysis determined the hub genes that played a critical role in differentially expressed genes between two distinct risk groups. In addition, the risk score was related to immunity status, immunotherapy response, and chemotherapeutic drug sensitivity. Finally, the expression of these signature CRLs detected by RT-qPCR was consistent with the bioinformatic analysis results.

CONCLUSION

In summary, our study confirmed that the novel CRL signature could effectively evaluate prognosis, tumor immune microenvironment, and immunotherapy response in OS. It may benefit for clinical decision-making and provide new insights for personalized therapeutics.

GPR116 overexpression correlates with poor prognosis in gastric cancer

The G protein-coupled receptor 116 (GPR116) is an adhesion subtype of the G protein-coupled receptor family and has been reported to be involved in tumorigenesis and cancer progression. Moreover, it has been shown to predict poor prognosis in breast and colorectal cancers. However, little is known about the role of GPR116 in gastric cancer (GC). In this study, we aimed to investigate the expression and clinical prognostic significance of GPR116 in GC.The mRNA expression levels of GPR116 in GC were analyzed using Gene Expression Omnibus and UALCAN databases, and GPR116 protein expression in GC tissues was detected using immunohistochemistry. The relationship between GPR116 expression and prognosis in patients with GC was analyzed and further validated using the Kaplan-Meier Plotter database. The correlation between GPR116 and the differentially expressed genes identified was analyzed using the LinkedOmics database. Gene set enrichment analysis was performed using WebGestalt. The results revealed that GPR116 expression was significantly upregulated in GC tissues, which was positively correlated with tumor node metastasis (TNM) staging and tumor invasion. Prognostic analysis suggested that high GPR116 expression contributed to poor overall survival in GC patients. Multivariate Cox analysis indicated that GPR116 overexpression was an independent prognostic indicator in patients with GC (HR = 1.855, 95% CI 1.021-3.370, P = .043). Enrichment analysis showed that GPR116 co-expression genes were mainly involved in extracellular matrix-receptor interaction, focal adhesion, cell adhesion, PI3K-Akt signaling, DNA replication, and cell cycle pathways. In conclusion, GPR116 was highly expressed in GC tissues and associated with poor prognosis in patients with GC, Thus GPR116 may be a novel prognostic marker and a potential therapeutic target for GC treatment.

Development of prediction models for the sensitivity of oral squamous cell carcinomas to preoperative S-1 administration

S-1 is an anticancer agent that is comprised of tegafur, gimeracil, and oteracil potassium, and is widely used in various carcinomas including oral squamous cell carcinoma (OSCC). Although an established prediction tool is not available, we aimed to develop prediction models for the sensitivity of primary OSCC cases to the preoperative administration of S-1.

We performed DNA microarray analysis of 95 cases with OSCC. Using global gene expression data and the clinical data, we developed two different prediction models, namely, model 1 that comprised the complete response (CR) + the partial response (PR) versus stable disease (SD) + progressive disease (PD), and model 2 that comprised responders versus non-responders. Twelve and 18 genes were designated as feature genes (FGs) in models 1 and 2, respectively, and, of these, six genes were common to both models. The sensitivity was 96.3%, the specificity was 91.2%, and the accuracy was 92.6% for model 1, and the sensitivity was 95.6%, the specificity was 85.2%, and the accuracy was 92.6% for model 2. These models were validated using receiver operating characteristic analysis, and the areas under the curves were 0.967 and 0.949 in models 1 and 2, respectively. The data led to the development of models that can reliably predict the sensitivity of patients with OSCC to the preoperative administration of S-1. The mechanism that regulates S-1 sensitivity remains unclear; however, the prediction models developed provide hope that further functional investigations into the FGs will lead to a greater understanding of drug resistance.

ShRNA-based POLD2 expression knockdown sensitizes glioblastoma to DNA-Damaging therapeutics

Glioblastoma (GBM) has limited therapeutic options. DNA repair mechanisms contribute GBM cells to escape therapies and re-establish tumor growth. Multiple studies have shown that POLD2 plays a critical role in DNA replication, DNA repair and genomic stability. We demonstrate for the first time that POLD2 is highly expressed in human glioma specimens and that expression correlates with poor patient survival. siRNA or shRNA POLD2 inhibited GBM cell proliferation, cell cycle progression, invasiveness, sensitized GBM cells to chemo/radiation-induced cell death and reversed the cytoprotective effects of EGFR signaling. Conversely, forced POLD2 expression was found to induce GBM cell proliferation, colony formation, invasiveness and chemo/radiation resistance. POLD2 expression associated with stem-like cell subsets (CD133⁺ and SSEA-1⁺ cells) and positively correlated with Sox2 expression in clinical specimens. Its expression was induced by Sox2 and inhibited by the forced differentiation of GBM neurospheres. shRNA-POLD2 modestly inhibited GBM neurosphere-derived orthotopic xenografts growth, when combined with radiation, dramatically inhibited xenograft growth in a cooperative fashion. These novel findings identify POLD2 as a new potential therapeutic target for enhancing GBM response to current standard of care therapeutics.

Testing Mediation Effects in High-Dimensional Epigenetic Studies

Mediation analysis has been a powerful tool to identify factors mediating the association between exposure variables and outcomes. It has been applied to various genomic applications with the hope to gain novel insights into the underlying mechanism of various diseases. Given the high-dimensional nature of epigenetic data, recent effort on epigenetic mediation analysis is to first reduce the data dimension by applying high-dimensional variable selection techniques, then conducting testing in a low dimensional setup. In this paper, we propose to assess the mediation effect by adopting a high-dimensional testing procedure which can produce unbiased estimates of the regression coefficients and can properly handle correlations between variables. When the data dimension is ultra-high, we first reduce the data dimension from ultra-high to high by adopting a sure independence screening (SIS) method. We apply the method to two high-dimensional epigenetic studies: one is to assess how DNA methylations mediate the association between alcohol consumption and epithelial ovarian cancer (EOC) status; the other one is to assess how methylation signatures mediate the association between childhood maltreatment and post-traumatic stress disorder (PTSD) in adulthood. We compare the performance of the method with its counterpart via simulation studies. Our method can be applied to other high-dimensional mediation studies where high-dimensional mediation variables are collected.

Gene set enrichment analysis and meta-analysis identified 12 key genes regulating and controlling the prognosis of lung adenocarcinoma

The aim of the present study was to analyze lung adenocarcinoma-associated microarray data and identify potentially crucial genes. The gene expression profiles were downloaded from the Gene Expression Omnibus database and 6 datasets, of which 2 were discarded and 4 were retained, were preprocessed using packages in the R computing language. Subsequently, Gene Set Enrichment Analysis (GSEA) and meta-analysis was used to screen the common pathways and differentially expressed genes at the transcriptional level. The genes detected from GSEA through The Cancer Genome Atlas databases were subsequently examined, and the crucial genes by survival data were identified. Pathways of the crucial genes were obtained using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway of the online website Database for Annotation, Visualization and Integrated Discovery (DAVID) tool, and the pathways of crucial genes that were upregulated or downregulated were matched using the Venn method to identify the common crucial pathways. Furthermore, on the basis of the common crucial pathways, key genes that are closely associated with the development and progression of lung adenocarcinoma were identified with the KEGG pathway of DAVID. Additional information was obtained through Gene Ontology annotation. A total of two key pathways, including cell cycle and DNA replication, as well as 12 key genes [DNA polymerase δ subunit 2, DNA replication licensing factor MCM4, MCM6, mitotic checkpoint serine/threonine-protein kinase BUB1, BUB1β, mitotic spindle assembly checkpoint protein MAD2A, dual specificity protein kinase TTK, M-phase inducer phosphatase 1, cell division control protein 45 homolog, cyclin-dependent kinase inhibitor 1C, pituitary tumor-transforming gene 1 protein and polo-like kinase 1] were identified. These key pathways and genes may be studied in future studies involving gene transfection/knockdown, which may provide insights into the prognosis of lung adenocarcinoma. Additional studies are required to confirm their biological function.

CAD/POLD2 gene expression is associated with poor overall survival and chemoresistance in bladder urothelial carcinoma

Somatic mutations in DNA repair genes have been clinically associated with chemosensitivity, although few studies have interrogated the nucleotide synthesis pathways that supply DNA repair processes. Previous work suggests that bladder urothelial carcinoma is uniquely enriched for mutations in nucleotide excision repair genes, and that these mutations are associated with response to platinum-based therapy and favorable survival. Conversely, the de novo pyrimidine synthesis pathway has recently emerged as a putative clinical target. This anabolic process is thought to supply DNA repair processes such as nucleotide excision repair; that is, DNA repair enzymes may require a sufficient nucleotide supply available to reverse the intended genotoxic damage of systemic chemotherapy in rapidly proliferating cancer cells. Therefore, we explored the prognostic complementarity between de novo pyrimidine synthesis and nucleotide excision repair expression in a total of 570 bladder urothelial carcinoma patients. Ultimately, we show that the de novo pyrimidine synthesis gene CAD is associated with poor survival (P = 0.008) and is co-altered with the nucleotide excision repair gene POLD2. High expression of POLD2 was also associated with poor overall survival (P = 0.019) and was significantly correlated with CAD expression in pre-treatment patient tumor samples (P = 2.44e-4). Expression of each gene was associated with cisplatin-based therapy resistance, and accordingly, CADhighPOLD2high patients were associated with worse survival than CADhighPOLD2low and CADlowPOLD2high patients. Together, these biomarkers could help elucidate mechanisms of chemoresistance to further personalize therapeutic strategies in bladder urothelial carcinoma.

FGF binding proteins (FGFBPs): Modulators of FGF signaling in the developing, adult, and stressed nervous system

Members of the fibroblast growth factor (FGF) family are involved in a variety of cellular processes. In the nervous system, they affect the differentiation and migration of neurons, the formation and maturation of synapses, and the repair of neuronal circuits following insults. Because of the varied yet critical functions of FGF ligands, their availability and activity must be tightly regulated for the nervous system, as well as other tissues, to properly develop and function in adulthood. In this regard, FGF binding proteins (FGFBPs) have emerged as strong candidates for modulating the actions of secreted FGFs in neural and non-neural tissues. Here, we will review the roles of FGFBPs in the peripheral and central nervous systems.

Novel Mad2-targeting miR-493-3p controls mitotic fidelity and cancer cells' sensitivity to paclitaxel

The molecular pathways that contribute to the proliferation and drug response of cancer cells are highly complex and currently insufficiently characterized. We have identified a previously unknown microRNA-based mechanism that provides cancer cells means to stimulate tumorigenesis via increased genomic instability and, at the same time, evade the action of clinically utilized microtubule drugs. We demonstrate miR-493-3p to be a novel negative regulator of mitotic arrest deficient-2 (MAD2), an essential component of the spindle assembly checkpoint that monitors the fidelity of chromosome segregation. The microRNA targets the 3′ UTR of Mad2 mRNA thereby preventing translation of the Mad2 protein. In cancer cells, overexpression of miR-493-3p induced a premature mitotic exit that led to increased frequency of aneuploidy and cellular senescence in the progeny cells. Importantly, excess of the miR-493-3p conferred resistance of cancer cells to microtubule drugs. In human neoplasms, miR-493-3p and Mad2 expression alterations correlated with advanced ovarian cancer forms and high miR-493-3p levels were associated with reduced survival of ovarian and breast cancer patients with aggressive tumors, especially in the paclitaxel therapy arm. Our results suggest that intratumoral profiling of miR-493-3p and Mad2 levels can have diagnostic value in predicting the efficacy of taxane chemotherapy.

Global gene expression changes in the prefrontal cortex of rabbits with hypercholesterolemia and/or hypertension

Although many studies have identified a link between hypercholesterolemia or hypertension and cognitive deficits, till date, comprehensive gene expression analyses of the brain under these conditions is still lacking. The present study was carried out to elucidate differential gene expression changes in the prefrontal cortex (PFC) of New Zealand white rabbits exposed to hypercholesterolemia and/or hypertension with a view of identifying gene networks at risk. Microarray analyses of the PFC of hypercholesterolemic rabbits showed 850 differentially expressed genes (DEGs) in the cortex of hypercholesterolemic rabbits compared to controls, but only 5 DEGs in hypertensive rabbits compared to controls. Up-regulated genes in the PFC of hypercholesterolemic rabbits included CIDEC, ODF2, RNASEL, FSHR, CES3 and MAB21L3, and down-regulated genes included FAM184B, CUL3, LOC100351029, TMEM109, LOC100357097 and PFDN5. Comparison with our previous study on the middle cerebral artery (MCA) of the same rabbits showed many differentially expressed genes in common between the PFC and MCA, during hypercholesterolemia. Moreover, these genes tended to fall into the same functional networks, as revealed by IPA analyses, with many identical node molecules. These include: proteasome, insulin, Akt, ERK1/2, histone, IL12, interferon alpha and NFκB. Of these, PSMB4, PSMD4, PSMG1 were chosen as representatives of genes related to the proteasome for verification by quantitative RT-PCR. Results indicate significant downregulation of all three proteasome associated genes in the PFC. Immunostaining showed significantly increased number of Aβ labelled cells in layers III and V of the cortex after hypercholesterolemia and hypertension, which may be due to decreased proteasome activity and/or increased β- or γ-secretase activity. Knowledge of altered gene networks during hypercholesterolemia and/or hypertension could inform our understanding of the link between these conditions and cognitive deficits in vascular dementia or Alzheimer's disease.

A network-pathway based module identification for predicting the prognosis of ovarian cancer patients

Background

This study aimed to screen multiple genes biomarkers based on gene expression data for predicting the survival of ovarian cancer patients.

Methods

Two microarray data of ovarian cancer samples were collected from The Cancer Genome Atlas (TCGA) database. The data in the training set were used to construct Reactome functional interactions network, which then underwent Markov clustering, supervised principal components, Cox proportional hazard model to screen significantly prognosis related modules. The distinguishing ability of each module for survival was further evaluated by the testing set. Gene Ontology (GO) functional and pathway annotations were performed to identify the roles of genes in each module for ovarian cancer.

Results

The network based approach identified two 7-gene functional interaction modules (31: DCLRE1A, EXO1, KIAA0101, KIN, PCNA, POLD3, POLD2; 35: DKK3, FABP3, IRF1, AIM2, GBP1, GBP2, IRF2) that are associated with prognosis of ovarian cancer patients. These network modules are related to DNA repair, replication, immune and cytokine mediated signaling pathways.

Conclusions

The two 7-gene expression signatures may be accurate predictors of clinical outcome in patients with ovarian cancer and has the potential to develop new therapeutic strategies for ovarian cancer patients.

Gonadal transcriptomic analysis of yellow catfish (Pelteobagrus fulvidraco): identification of sex-related genes and genetic markers

Yellow catfish (Pelteobagrus fulvidraco) has been recognized as a vital freshwater aquaculture species in East and Southeast Asia. In addition to its commercial interest, it is also attracted much attention because of its value in studying sex-determination mechanisms. A comprehensive gonadal transcriptome analysis is believed to provide a resource for genome annotation, candidate gene identification, and molecular marker development. Herein, we performed a de novo assembly of yellow catfish gonad transcriptome by high-throughput Illumina sequencing. A total of 82,123 contigs were obtained, ranging from 351 to 21,268 bp, and N50 of 2,329 bp. Unigenes of 21,869 in total were identified. Of these, 229 and 1,188 genes were found to be specifically expressed in XY gonad tissue for 1 yr and 2 yr old yellow catfish, respectively; correspondingly, 51 and 40 genes were identified in XX gonad tissue at those two stages. Gene ontology and KEGG analysis were conducted and classified all contigs into different categories. A large number of unigenes involved in sex determination were identified, as well as microsatellites and SNP variants. The expression patterns of sex-related genes were then validated by quantitative real-time PCR (qRT-PCR) suggesting the high reliability of RNA-Seq results. In this study, the transcriptome of yellow catfish gonad was first sequenced, assembled, and characterized; it provides a valuable genomic resource for better understanding of yellow catfish sex determination as well as development of molecular markers, thereby assisting in the production of monosex yellow catfish for aquaculture.

BUB1 mRNA is significantly co-expressed with AURKA and AURKB mRNA in advanced-stage ovarian serous carcinoma

The objective of this study was to investigate the expression and clinical role of the spindle checkpoint kinase budding uninhibited by benzimidazole 1 (Bub1) in primary and metastatic advanced-stage ovarian serous carcinoma. BUB1 mRNA expression was analyzed in 178 tumors (88 effusions, 38 primary carcinomas, and 52 solid metastases) from 144 patients with advanced-stage disease using quantitative real-time polymerase chain reaction (PCR). Bub1 protein expression by Western blotting was studied in 63 carcinomas (30 effusions and 33 solid lesions). BUB1 mRNA expression at different anatomic sites was studied for association with clinicopathologic parameters, including chemotherapy resistance and survival. BUB1 mRNA was universally expressed in serous carcinomas, irrespective of anatomic site. BUB1 mRNA levels were uniformly low in six ovarian surface epithelium specimens analyzed for comparative purposes. Bub1 protein was expressed in 22/30 effusions and 28/33 solid lesions. BUB1 mRNA expression was significantly higher in chemo-naïve primary carcinomas and solid metastases compared to specimens obtained following neoadjuvant chemotherapy (p < 0.001) and was unrelated to chemotherapy exposure in effusions nor to chemoresponse or survival at any anatomic site. BUB1 mRNA levels in both effusions and solid lesions were strongly related to the mRNA levels of AURKA and AURKB previously studied in this cohort (p < 0.001 for both). Bub1 is widely expressed in primary and metastatic OC, suggesting a biological role in this cancer. BUB1 mRNA levels are lower following chemotherapy exposure in solid lesions, though its presence is unrelated to clinical behavior including response to chemotherapy and survival. BUB1 is co-expressed with AURKA and AURKB suggesting biological relationship between these spindle cell components.

Global miRNA expression analysis of serous and clear cell ovarian carcinomas identifies differentially expressed miRNAs including miR-200c-3p as a prognostic marker

BACKGROUND

Improved insight into the molecular characteristics of the different ovarian cancer subgroups is needed for developing a more individualized and optimized treatment regimen. The aim of this study was to a) identify differentially expressed miRNAs in high-grade serous ovarian carcinoma (HGSC), clear cell ovarian carcinoma (CCC) and ovarian surface epithelium (OSE), b) evaluate selected miRNAs for association with clinical parameters including survival and c) map miRNA-mRNA interactions.

METHODS

Differences in miRNA expression between HGSC, CCC and OSE were analyzed by global miRNA expression profiling (Affymetrix GeneChip miRNA 2.0 Arrays, n = 12, 9 and 9, respectively), validated by RT-qPCR (n = 35, 19 and 9, respectively), and evaluated for associations with clinical parameters. For HGSC, differentially expressed miRNAs were linked to differentially expressed mRNAs identified previously.

RESULTS

Differentially expressed miRNAs (n = 78) between HGSC, CCC and OSE were identified (FDR < 0.01%), of which 18 were validated (p < 0.01) using RT-qPCR in an extended cohort. Compared with OSE, miR-205-5p was the most overexpressed miRNA in HGSC. miR-200 family members and miR-182-5p were the most overexpressed in HGSC and CCC compared with OSE, whereas miR-383 was the most underexpressed. miR-205-5p and miR-200 members target epithelial-mesenchymal transition (EMT) regulators, apparently being important in tumor progression. miR-509-3-5p, miR-509-5p, miR-509-3p and miR-510 were among the strongest differentiators between HGSC and CCC, all being significantly overexpressed in CCC compared with HGSC. High miR-200c-3p expression was associated with poor progression-free (p = 0.031) and overall (p = 0.026) survival in HGSC patients. Interacting miRNA and mRNA targets, including those of a TP53-related pathway presented previously, were identified in HGSC.

CONCLUSIONS

Several miRNAs differentially expressed between HGSC, CCC and OSE have been identified, suggesting a carcinogenetic role for these miRNAs. miR-200 family members, targeting EMT drivers, were mostly overexpressed in both subgroups, among which miR-200c-3p was associated with survival in HGSC patients. A set of miRNAs differentiates CCC from HGSC, of which miR-509-3-5p and miR-509-5p are the strongest classifiers. Several interactions between miRNAs and mRNAs in HGSC were mapped.

Proteomic comparison of mcf-7 tumoursphere and monolayer cultures

Breast cancer is a heterogenous disease, composed of tumour cells with differing gene expressions and phenotypes. Very few antigens have been identified and a better understanding of tumour initiating-cells as targets for therapy is critically needed. Recently, a rare subpopulation of cells within tumours has been described with the ability to: (i) initiate and sustain tumour growth; (ii) resist traditional therapies and allow for secondary tumour dissemination; and (iii) display some of the characteristics of stem cells such as self-renewal. These cells are termed tumour-initiating cells or cancer stem cells, or alternatively, in the case of breast cancer, breast cancer stem cells. Previous studies have demonstrated that breast cancer stem cells can be enriched for in “tumoursphere” culture. Proteomics represents a novel way to investigate protein expression between cells. We hypothesise that characterisation of the proteome of the breast cancer line MCF-7 tumourspheres compared to adherent/differentiated cells identifies proteins of novel interest for further isolating or targeting breast cancer stem cells. We present evidence that: (i) the proteome of adherent cells is different to the proteome of cells grown in sphere medium from either early passage (passage 2) or late passage (passage 5) spheres; (ii) that spheres are enriched in expression of a variety of tumour-relevant proteins (including MUC1 and Galectin-3); and (iii) that targeting of one of these identified proteins (galectin-3) using an inhibitor (N-acetyllactosamine) decreases sphere formation/self-renewal of MCF-7 cancer stem cells in vitro and tumourigenicity in vivo. Hence, proteomic analysis of tumourspheres may find use in identifying novel targets for future therapy. The therapeutic targeting of breast cancer stem cells, a highly clinically relevant sub-population of tumour cells, has the potential to eliminate residual disease and may become an important component of a multi-modality treatment of cancer.

Novel therapeutic vaccines [(HSP65 + IL-12)DNA-, granulysin- and Ksp37-vaccine] against tuberculosis and synergistic effects in the combination with chemotherapy

PURPOSE

Multi-drug resistant tuberculosis (MDR-TB) and extremely drug resistant (XDR) TB are big problems in the world. We have developed novel TB therapeutic vaccines, HVJ-Envelope/HSP65 + IL-12 DNA vaccine (HSP65-vaccine), granulysin vaccine and killer specific secretory protein of 37kDa (Ksp37) vaccine.

METHODS AND RESULTS

HSP65 vaccine showed strong therapeutic effect against both MDR-TB and XDR-TB in mice. Intradermal immunization of HSP65-vaccine showed stronger therapeutic effect against TB than intramuscular or subcutaneous immunization. Furthermore, the synergistic therapeutic effect was observed when the vaccine was administrated in combination with Isoniazid (INH), which is a first line drug for chemotherapy. The combination of types of vaccines (HSP65- and granulysin- vaccines) also showed synergistic therapeutic effect. In the monkey model, granulysin-vaccine prolonged the survival period after the infection of TB and long-term survival was observed in vaccine-treated group. We examined the potential of two kinds of novel DNA vaccines (Ksp37-vaccine and granulysin-vaccine). Both vaccines augmented in vivo differentiation of CTL against TB. We measured the amount of Ksp37 protein in human serum and revealed that the level of Ksp37 protein of patients with tuberculosis was lower than that of healthy volunteers. Therefore, we established Ksp37 transgenic mice as well as granulysin transgenic mice to elucidate the function of those proteins. Both transgenic mice were resistant to TB infection.

CONCLUSION

These data indicate the potential of combinational therapy; the combination of two DNA vaccines or combination of DNA vaccine with antibiotic drug. Thus, it will provide a novel strategy for the treatment of MDR-TB.

The study of novel DNA vaccines against tuberculosis: induction of pathogen-specific CTL in the mouse and monkey models of tuberculosis

RESULTS

HSP65 + IL-12 DNA vaccine showed higher protective efficacy compared with BCG in both mouse and monkey models of TB. It induced the TB-specific CTL in the mouse model of TB, while little level of activity was observed after the injection of BCG. It also showed strong therapeutic efficacy against MDR-TB. In the monkey model, the vaccine augmented the production of IFN-γ and IL-2 from PBL and the therapeutic effect was correlated with the level of IL-2. We next evaluated the potential of DNA vaccine encoding a granulysin, which is an important defensive molecule expressed by human T cells. We found that granulysin-encoding vaccine induced the differentiation of the CTL in vitro and in vivo. It also showed therapeutic efficacy against TB in the monkey as well as the mouse model. The DNA vaccine encoding a Ksp37 also induced the TB-specific CTL in vitro and in vivo in the mouse model. It augmented the production of IL-2, IFN-γ and IL-6 from T cells and spleen cells. A synergistic effect on the activation of the TB-specific CTL was observed by the combination of Ksp37 DNA vaccine with granulysin DNA vaccine.

PURPOSE AND METHODS

Emergence of the multi-drug resistant (MDR) Mycobacterium tuberculosis (TB) is a big problem in the world. We have developed novel TB vaccines [DNA vaccines encoding HSP65 + IL-12, granulysin or killer-specific secretory protein of 37kDa (Ksp37)] using Hemagglutinating virus of Japan -envelope (HVJ-E). It is suggested that the activity of the TB-specific CTL is one of the most important factor for the resistance to TB and immunity for TB in chronic human TB disease. Therefore, we examined the level of activation of the TB-specific CTL after the administration of these vaccines.

CONCLUSION

These data indicate that our novel vaccines (HSP65 + IL-12 DNA, granulysin and Ksp37) have a capability to activate the TB-specific CTL and will be very strong protective and therapeutic vaccines against TB.

ZNF385B and VEGFA are strongly differentially expressed in serous ovarian carcinomas and correlate with survival

BACKGROUND

The oncogenesis of ovarian cancer is poorly understood. The aim of this study was to identify mRNAs differentially expressed between moderately and poorly differentiated (MD/PD) serous ovarian carcinomas (SC), serous ovarian borderline tumours (SBOT) and superficial scrapings from normal ovaries (SNO), and to correlate these mRNAs with clinical parameters including survival.

METHODS

Differences in mRNA expression between MD/PD SC, SBOT and SNO were analyzed by global gene expression profiling (n = 23), validated by RT-qPCR (n = 41) and correlated with clinical parameters.

RESULTS

Thirty mRNAs differentially expressed between MD/PD SC, SBOT and SNO were selected from the global gene expression analyses, and 21 were verified (p<0.01) by RT-qPCR. Of these, 13 mRNAs were differentially expressed in MD/PD SC compared with SNO (p<0.01) and were correlated with clinical parameters. ZNF385B was downregulated (FC = -130.5, p = 1.2×10(-7)) and correlated with overall survival (p = 0.03). VEGFA was upregulated (FC = 6.1, p = 6.0×10(-6)) and correlated with progression-free survival (p = 0.037). Increased levels of TPX2 and FOXM1 mRNAs (FC = 28.5, p = 2.7×10(-10) and FC = 46.2, p = 5.6×10(-4), respectively) correlated with normalization of CA125 (p = 0.03 and p = 0.044, respectively). Furthermore, we present a molecular pathway for MD/PD SC, including VEGFA, FOXM1, TPX2, BIRC5 and TOP2A, all significantly upregulated and directly interacting with TP53.

CONCLUSIONS

We have identified 21 mRNAs differentially expressed (p<0.01) between MD/PD SC, SBOT and SNO. Thirteen were differentially expressed in MD/PD SC, including ZNF385B and VEGFA correlating with survival, and FOXM1 and TPX2 with normalization of CA125. We also present a molecular pathway for MD/PD SC.