NgR1 is an NK cell inhibitory receptor that destabilizes the immunological synapse.

The formation of an immunological synapse (IS) is essential for natural killer (NK) cells to eliminate target cells. Despite an advanced understanding of the characteristics of the IS and its formation processes, the mechanisms that regulate its stability via the cytoskeleton are unclear. Here, we show that Nogo receptor 1 (NgR1) has an important function in modulating NK cell-mediated killing by destabilization of IS formation. NgR1 deficiency or blockade resulted in improved tumor control of NK cells by enhancing NK-to-target cell contact stability and regulating F-actin dynamics during IS formation. Patients with tumors expressing abundant NgR1 ligand had poor prognosis despite high levels of NK cell infiltration. Thus, our study identifies NgR1 as an immune checkpoint in IS formation and indicates a potential approach to improve the cytolytic function of NK cells in cancer immunotherapy.

Transcriptional Profiling of Advanced Urothelial Cancer Predicts Prognosis and Response to Immunotherapy.

Recent investigations reported that some subtypes from the Lund or The Cancer Genome Atlas (TCGA) classifications were most responsive to PD-L1 inhibitor treatment. However, the association between previously reported subtypes and immune checkpoint inhibitor (ICI) therapy responsiveness has been insufficiently explored. Despite these contributions, the ability to predict the clinical applicability of immune checkpoint inhibitor therapy in patients remains a major challenge. Here, we aimed to re-classify distinct subtypes focusing on ICI responsiveness using gene expression profiling in the IMvigor 210 cohort (n = 298). Based on the hierarchical clustering analysis, we divided advanced urothelial cancer patients into three subgroups. To confirm a prognostic impact, we performed survival analysis and estimated the prognostic value in the IMvigor 210 and TCGA cohort. The activation of CD8+ T effector cells was common for patients of classes 2 and 3 in the TCGA and IMvigor 210 cohort. Survival analysis showed that patients of class 3 in the TCGA cohort had a poor prognosis, while patients of class 3 showed considerably prolonged survival in the IMvigor 210 cohort. One of the distinct characteristics of patients in class 3 is the inactivation of the TGFß and YAP/TAZ pathways and activation of the cell cycle and DNA replication and DNA damage (DDR). Based on our identified transcriptional patterns and the clinical outcomes of advanced urothelial cancer patients, we constructed a schematic summary. When comparing clinical and transcriptome data, patients with downregulation of the TGFß and YAP/TAZ pathways and upregulation of the cell cycle and DDR may be more responsive to ICI therapy.

E2F1 Promotes Progression of Bladder Cancer by Modulating RAD54L Involved in Homologous Recombination Repair.

DNA repair defects are important factors in cancer development. High DNA repair activity can affect cancer progression and chemoresistance. DNA double-strand breaks in cancer cells caused by anticancer agents can be restored by non-homologous end joining (NHEJ) and homologous recombination repair (HRR). Our previous study has identified E2F1 as a key gene in bladder cancer progression. In this study, DNA repair genes related to E2F1 were analyzed, and RAD54L involved in HRR was identified. In gene expression analysis of bladder cancer patients, the survival of patients with high RAD54L expression was shorter with cancer progression than in patients with low RAD54L expression. This study also revealed that E2F1 directly binds to the promoter region of RAD54L and regulates the transcription of RAD54L related to the HRR pathway. This study also confirmed that DNA breaks are repaired by RAD54L induced by E2F1 in bladder cancer cells treated with MMC. In summary, RAD54L was identified as a new target directly regulated by E2F1. Our results suggest that, E2F1 and RAD54L could be used as diagnostic markers for bladder cancer progression and represent potential therapeutic targets.

Genomic predictors for recurrence patterns of hepatocellular carcinoma: model derivation and validation.

BACKGROUND: Typically observed at 2 y after surgical resection, late recurrence is a major challenge in the management of hepatocellular carcinoma (HCC). We aimed to develop a genomic predictor that can identify patients at high risk for late recurrence and assess its clinical implications. METHODS AND FINDINGS: Systematic analysis of gene expression data from human liver undergoing hepatic injury and regeneration revealed a 233-gene signature that was significantly associated with late recurrence of HCC. Using this signature, we developed a prognostic predictor that can identify patients at high risk of late recurrence, and tested and validated the robustness of the predictor in patients (n = 396) who underwent surgery between 1990 and 2011 at four centers (210 recurrences during a median of 3.7 y of follow-up). In multivariate analysis, this signature was the strongest risk factor for late recurrence (hazard ratio, 2.2; 95% confidence interval, 1.3-3.7; p = 0.002). In contrast, our previously developed tumor-derived 65-gene risk score was significantly associated with early recurrence (p = 0.005) but not with late recurrence (p = 0.7). In multivariate analysis, the 65-gene risk score was the strongest risk factor for very early recurrence (<1 y after surgical resection) (hazard ratio, 1.7; 95% confidence interval, 1.1-2.6; p = 0.01). The potential significance of STAT3 activation in late recurrence was predicted by gene network analysis and validated later. We also developed and validated 4- and 20-gene predictors from the full 233-gene predictor. The main limitation of the study is that most of the patients in our study were hepatitis B virus-positive. Further investigations are needed to test our prediction models in patients with different etiologies of HCC, such as hepatitis C virus. CONCLUSIONS: Two independently developed predictors reflected well the differences between early and late recurrence of HCC at the molecular level and provided new biomarkers for risk stratification. Please see later in the article for the Editors' Summary.

Lipocalin-2 negatively modulates the epithelial-to-mesenchymal transition in hepatocellular carcinoma through the epidermal growth factor (TGF-beta1)/Lcn2/Twist1 pathway.

UNLABELLED: Lipocalin-2 (Lcn2) is preferentially expressed in hepatocellular carcinoma (HCC). However, the functional role of Lcn2 in HCC progression is still poorly understood, particularly with respect to its involvement in invasion and metastasis. The purpose of this study was to investigate whether Lcn2 is associated with the epithelial-mesenchymal transition (EMT) in HCC and to elucidate the underlying signaling pathway(s). Lcn2 was preferentially expressed in well-differentiated HCC versus liver cirrhosis tissues, and its expression was positively correlated with the stage of HCC. The characteristics of EMT were reversed by adenoviral transduction of Lcn2 into SH-J1 cells, including the down-regulation of N-cadherin, vimentin, alpha-smooth muscle actin, and fibronectin, and the concomitant up-regulation of CK8, CK18, and desmoplakin I/II. Knockdown of Lcn2 by short hairpin RNA (shRNA) in HKK-2 cells expressing high levels of Lcn2 was associated with EMT. Epidermal growth factor (EGF) or transforming growth factor beta1 (TGF-ß1) treatment resulted in down-regulation of Lcn2, accompanied by an increase in Twist1 expression and EMT in HCC cells. Stable Lcn2 expression in SH-J1 cells reduced Twist1 expression, inhibited cell proliferation and invasion in vitro, and suppressed tumor growth and metastasis in a mouse model. Furthermore, EGF or TGF-ß1 treatment barely changed EMT marker expression in SH-J1 cells ectopically expressing Lcn2. Ectopic expression of Twist1 induced EMT marker expression even in cells expressing Lcn2, indicating that Lcn2 functions downstream of growth factors and upstream of Twist1. CONCLUSION: Together, our findings indicate that Lcn2 can negatively modulate the EMT in HCC cells through an EGF (or TGF-ß1)/Lcn2/Twist1 pathway. Thus, Lcn2 may be a candidate metastasis suppressor and a potential therapeutic target in HCC.

Unraveling the role of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer by multi-omics analyses.

The role of the serine/glycine metabolic pathway (SGP) has recently been demonstrated in tumors; however, the pathological relevance of the SGP in thyroid cancer remains unexplored. Here, we perform metabolomic profiling of 17 tumor-normal pairs; bulk transcriptomics of 263 normal thyroid, 348 papillary, and 21 undifferentiated thyroid cancer samples; and single-cell transcriptomes from 15 cases, showing the impact of mitochondrial one-carbon metabolism in thyroid tumors. High expression of serine hydroxymethyltransferase-2 (SHMT2) and methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) is associated with low thyroid differentiation scores and poor clinical features. A subpopulation of tumor cells with high mitochondrial one-carbon pathway activity is observed in the single-cell dataset. SHMT2 inhibition significantly compromises mitochondrial respiration and decreases cell proliferation and tumor size in vitro and in vivo. Collectively, our results highlight the importance of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer and suggest that SHMT2 is a potent therapeutic target.

Histone deacetylase 6 functions as a tumor suppressor by activating c-Jun NH2-terminal kinase-mediated beclin 1-dependent autophagic cell death in liver cancer.

UNLABELLED: Ubiquitin-binding histone deacetylase 6 (HDAC6) is uniquely endowed with tubulin deacetylase activity and plays an important role in the clearance of misfolded protein by autophagy. In cancer, HDAC6 has become a target for drug development due to its major contribution to oncogenic cell transformation. In the present study we show that HDAC6 expression was down-regulated in a large cohort of human hepatocellular carcinoma (HCC) patients, and that low expression of HDAC6 was significantly associated with poor prognosis of HCC patients in 5-year overall, disease-free, and recurrence-free survival. Notably, we observed that ectopic overexpression of HDAC6 suppressed tumor cell growth and proliferation in various liver cancer cells, and elicited increased LC3B-II conversion and autophagic vacuole formation without causing apoptotic cell death or cell cycle inhibition. In addition, the sustained overexpression of HDAC6 reduced the in vivo tumor growth rate in a mouse xenograft model. It was also found that HDAC6 mediated autophagic cell death by way of Beclin 1 and activation of the LC3-II pathway in liver cancer cells, and that HDAC6 overexpression activated c-Jun NH2-terminal kinase (JNK) and increased the phosphorylation of c-Jun. In contrast, the induction of Beclin 1 expression was blocked by SP600125 (a specific inhibitor of JNK) or by small interfering RNA directed against HDAC6. CONCLUSION: Our findings suggest that loss of HDAC6 expression in human HCCs and tumor suppression by HDAC6 occur by way of activation of caspase-independent autophagic cell death through the JNK/Beclin 1 pathway in liver cancer and, thus, that a novel tumor suppressor function mechanism involving HDAC6 may be amenable to nonepigenetic regulation.

Sixty-five gene-based risk score classifier predicts overall survival in hepatocellular carcinoma.

UNLABELLED: Clinical application of the prognostic gene expression signature has been delayed due to the large number of genes and complexity of prediction algorithms. In the current study we aimed to develop an easy-to-use risk score with a limited number of genes that can robustly predict prognosis of patients with hepatocellular carcinoma (HCC). The risk score was developed using Cox coefficient values of 65 genes in the training set (n = 139) and its robustness was validated in test sets (n = 292). The risk score was a highly significant predictor of overall survival (OS) in the first test cohort (P = 5.6 × 10(-5), n = 100) and the second test cohort (P = 5.0 × 10(-5) , n = 192). In multivariate analysis, the risk score was a significant risk factor among clinical variables examined together (hazard ratio [HR], 1.36; 95% confidence interval [CI], 1.13-1.64; P = 0.001 for OS). CONCLUSION: The risk score classifier we have developed can identify two clinically distinct HCC subtypes at early and late stages of the disease in a simple and highly reproducible manner across multiple datasets.

A novel prognostic subtype of human hepatocellular carcinoma derived from hepatic progenitor cells.

The variability in the prognosis of individuals with hepatocellular carcinoma (HCC) suggests that HCC may comprise several distinct biological phenotypes. These phenotypes may result from activation of different oncogenic pathways during tumorigenesis and/or from a different cell of origin. Here we address whether the transcriptional characteristics of HCC can provide insight into the cellular origin of the tumor. We integrated gene expression data from rat fetal hepatoblasts and adult hepatocytes with HCC from human and mouse models. Individuals with HCC who shared a gene expression pattern with fetal hepatoblasts had a poor prognosis. The gene expression program that distinguished this subtype from other types of HCC included markers of hepatic oval cells, suggesting that HCC of this subtype may arise from hepatic progenitor cells. Analyses of gene networks showed that activation of AP-1 transcription factors in this newly identified HCC subtype might have key roles in tumor development.

Application of comparative functional genomics to identify best-fit mouse models to study human cancer.

Genetically modified mice have been extensively used for analyzing the molecular events that occur during tumor development. In many, if not all, cases, however, it is uncertain to what extent the mouse models reproduce features observed in the corresponding human conditions. This is due largely to lack of precise methods for direct and comprehensive comparison at the molecular level of the mouse and human tumors. Here we use global gene expression patterns of 68 hepatocellular carcinomas (HCCs) from seven different mouse models and 91 human HCCs from predefined subclasses to obtain direct comparison of the molecular features of mouse and human HCCs. Gene expression patterns in HCCs from Myc, E2f1 and Myc E2f1 transgenic mice were most similar to those of the better survival group of human HCCs, whereas the expression patterns in HCCs from Myc Tgfa transgenic mice and in diethylnitrosamine-induced mouse HCCs were most similar to those of the poorer survival group of human HCCs. Gene expression patterns in HCCs from Acox1(-/-) mice and in ciprofibrate-induced HCCs were least similar to those observed in human HCCs. We conclude that our approach can effectively identify appropriate mouse models to study human cancers.

Prognostic discrimination using a 70-gene signature among patients with estrogen receptor-positive breast cancer and an intermediate 21-gene recurrence score.

The Oncotype DX® recurrence score (RS) predictor has been clinically utilized to appropriately select adjuvant chemotherapy for patients with estrogen receptor (ER)-positive early breast cancer. However, the selection of chemotherapy for patients with intermediate RSs remains controversial. We assessed the prognostic value of a 70-gene signature (70GS) among patients with ER-positive breast cancer and intermediate RSs. In addition, we sought to identify genes associated with poor 70GS scores based on gene expression profiling (GEP). GEP was performed using gene expression data from 186 patients with ER-positive breast cancer. The RS and 70GS score were calculated on the basis of GEP. Among 186 patients, 82 ER-positive patients with intermediate RSs were identified. These patients were stratified by 70GS, overall survival (OS) significantly differed according to 70GS (p=0.013). In a supervised hierarchical analysis according to 70GS, the expression of several representative genes for cell proliferation was significantly higher in the poor 70GS cluster than in the good 70GS cluster. Furthermore, among these patients, FOXM1, AURKA, AURKB, and BIRC5 displayed prognostic significance for OS. In conclusion, 70GS can help to discriminate survival differences among ER-positive patients with intermediate RSs. FOXM1, AURKA, AURKB, and BIRC5, are associated with poor 70GS scores.

Prognostic gene expression signature associated with two molecularly distinct subtypes of colorectal cancer.

AIMS: Despite continual efforts to develop prognostic and predictive models of colorectal cancer by using clinicopathological and genetic parameters, a clinical test that can discriminate between patients with good or poor outcome after treatment has not been established. Thus, the authors aim to uncover subtypes of colorectal cancer that have distinct biological characteristics associated with prognosis and identify potential biomarkers that best reflect the biological and clinical characteristics of subtypes. METHODS: Unsupervised hierarchical clustering analysis was applied to gene expression data from 177 patients with colorectal cancer to determine a prognostic gene expression signature. Validation of the signature was sought in two independent patient groups. The association between the signature and prognosis of patients was assessed by Kaplan-Meier plots, log-rank tests and the Cox model. RESULTS: The authors identified a gene signature that was associated with overall survival and disease-free survival in 177 patients and validated in two independent cohorts of 213 patients. In multivariate analysis, the signature was an independent risk factor (HR 3.08; 95% CI 1.33 to 7.14; p=0.008 for overall survival). Subset analysis of patients with AJCC (American Joint Committee on Cancer) stage III cancer revealed that the signature can also identify the patients who have better outcome with adjuvant chemotherapy (CTX). Adjuvant chemotherapy significantly affected disease-free survival in patients in subtype B (3-year rate, 71.2% (CTX) vs 41.9% (no CTX); p=0.004). However, such benefit of adjuvant chemotherapy was not significant for patients in subtype A. CONCLUSION: The gene signature is an independent predictor of response to chemotherapy and clinical outcome in patients with colorectal cancer.

Single Cell Analysis of Human Thyroid Reveals the Transcriptional Signatures of Aging.

The thyroid gland plays a critical role in the maintenance of whole-body metabolism. However, aging frequently impairs homeostatic maintenance by thyroid hormones due to increased prevalence of subclinical hypothyroidism associated with mitochondrial dysfunction, inflammation, and fibrosis. To understand the specific aging-related changes of endocrine function in thyroid epithelial cells, we performed single-cell RNA sequencing (RNA-seq) of 54 726 cells derived from pathologically normal thyroid tissues from 7 patients who underwent thyroidectomy. Thyroid endocrine epithelial cells were clustered into 5 distinct subpopulations, and a subset of cells was found to be particularly vulnerable with aging, showing functional deterioration associated with the expression of metallothionein (MT) and major histocompatibility complex class II genes. We further validated that increased expression of MT family genes are highly correlated with thyroid gland aging in bulk RNAseq datasets. This study provides evidence that aging induces specific transcriptomic changes across multiple cell populations in the human thyroid gland.

Induction of pluripotent stem cells from adult somatic cells by protein-based reprogramming without genetic manipulation.

The concept of reprogramming of somatic cells has opened a new era in regenerative medicine. Transduction of defined factors has successfully achieved pluripotency. However, during the generation process of induced pluripotent stem (iPS) cells, genetic manipulation of certain factors may cause tumorigenicity, which limits further application. We report that that a single transfer of embryonic stem (ES) cell-derived proteins into primarily cultured adult mouse fibroblasts, rather than repeated transfer or prolonged exposure to materials, can achieve full reprogramming up to the pluripotent state without the forced expression of ectopic transgenes. During the process, gene expression and epigenetic status were converted from somatic to ES-equivalent status. We verified that protein-based reprogramming was neither by the contamination of protein donor ES cell nor by DNA/RNA from donor ES cell. Protein-iPS cells were biologically and functionally very similar to ES cells and differentiated into 3 germ layers in vitro. Furthermore, protein-iPS cells possessed in vivo differentiation (well-differentiated teratoma formation) and development (chimeric mice generation and a tetraploid blastocyst complementation) potentials. Our results provide an alternative and safe strategy for the reprogramming of somatic cells that can be used to facilitate pluripotent stem cell-based cell therapy.

Overexpression of Snail is associated with lymph node metastasis and poor prognosis in patients with gastric cancer.

BACKGROUND: Epithelial-mesenchymal transition (EMT) plays a significant role in tumor progression and invasion. Snail is a known regulator of EMT in various malignant tumors. This study investigated the role of Snail in gastric cancer. METHODS: We examined the effects of silenced or overexpressed Snail using lenti-viral constructs in gastric cancer cells. Immunohistochemical analysis of tissue microarrays from 314 patients with gastric adenocarcinoma (GC) was used to determine Snail's clinicopathological and prognostic significance. Differential gene expression in 45 GC specimens with Snail overexpression was investigated using cDNA microarray analysis. RESULTS: Silencing of Snail by shRNA decreased invasion and migration in GC cell lines. Conversely, Snail overexpression increased invasion and migration of gastric cancer cells, in line with increased VEGF and MMP11. Snail overexpression (≥75% positive nuclear staining) was also significantly associated with tumor progression (P < 0.001), lymph node metastases (P = 0.002), lymphovascular invasion (P = 0.002), and perineural invasion (P = 0.002) in the 314 GC patients, and with shorter survival (P = 0.023). cDNA microarray analysis revealed 213 differentially expressed genes in GC tissues with Snail overexpression, including genes related to metastasis and invasion. CONCLUSION: Snail significantly affects invasiveness/migratory ability of GCs, and may also be used as a predictive biomarker for prognosis or aggressiveness of GCs.

Sulfatase 1 and sulfatase 2 in hepatocellular carcinoma: associated signaling pathways, tumor phenotypes, and survival.

The heparin-degrading endosulfatases sulfatase 1 (SULF1) and sulfatase 2 (SULF2) have opposing effects in hepatocarcinogenesis despite structural similarity. Using mRNA expression arrays, we analyzed the correlations of SULF expression with signaling networks in human hepatocellular carcinomas (HCCs) and the associations of SULF expression with tumor phenotype and patient survival. Data from two mRNA microarray analyses of 139 and 36 HCCs and adjacent tissues were used as training and validation sets. Partek and Metacore software were used to identify SULF correlated genes and their associated signaling pathways. Associations between SULF expression, the hepatoblast subtype of HCC, and survival were examined. Both SULF1 and 2 had strong positive correlations with periostin, IQGAP1, TGFB1, and vimentin and inverse correlations with HNF4A and IQGAP2. Genes correlated with both SULFs were highly associated with the cell adhesion, cytoskeletal remodeling, blood coagulation, TGFB, and Wnt/ß-catenin and epithelial mesenchymal transition signaling pathways. Genes uniquely correlated with SULF2 were more associated with neoplastic processes than genes uniquely correlated with SULF1. High SULF expression was associated with the hepatoblast subtype of HCC. There was a bimodal effect of SULF1 expression on prognosis, with patients in the lowest or highest tertile having a worse prognosis than those in the middle tertile. SULFs have complex effects on HCC signaling and patient survival. There are functionally similar associations with cell adhesion, ECM remodeling, TGFB, and WNT pathways, but also unique associations of SULF1 and SULF2. The roles and targeting of the SULFs in cancer require further investigation.

Stepwise molecular mechanisms responsible for chemoresistance in bladder cancer cells.

Chemotherapy resistance is an obstacle to cancer therapy and is considered a major cause of recurrence. Thus, understanding the mechanisms of chemoresistance is critical to improving the prognosis of patients. Here, we have established a stepwise gemcitabine-resistant T24 bladder cancer cell line to understand the molecular mechanisms of chemoresistance within cancer cells. The characteristics of the stepwise chemoresistance cell line were divided into 4 phases (parental, early, intermediate, and late phases). These four phase cells showed increasingly aggressive phenotypes in vitro and in vivo experiments with increasing phases and revealed the molecular properties of the biological process from parent cells to phased gemcitabine-resistant cell line (GRC). Taken together, through the analysis of gene expression profile data, we have characterized gene set of each phase indicating the response to anticancer drug treatment. Specifically, we identified a multigene signature (23 genes including GATA3, APOBEC3G, NT5E, MYC, STC1, FOXD1, SMAD9) and developed a chemoresistance score consisting of that could predict eventual responsiveness to gemcitabine treatment. Our data will contribute to predicting chemoresistance and improving the prognosis of bladder cancer patients.

YAP1 activation is associated with the progression and response to immunotherapy of non-muscle invasive bladder cancer.

BACKGROUND: Despite the availability of several treatments for non-muscle-invasive bladder cancer (NMIBC), many patients are still not responsive to treatments, and the disease progresses. A new prognostic classifier can differentiate between treatment response and progression, and it could be used as a very important tool in patient decision-making regarding treatment options. In this study, we focused on the activation of Yes-associated protein 1 (YAP1), which is known to play a pivotal role in tumour progression and serves as a factor contributing to the mechanism of resistance to various relevant therapeutic agents. We further evaluated its potential as a novel prognostic agent. METHODS: We identified YAP1-associated gene signatures based on UC3-siYAP1 cells (n=8) and NMIBC cohort (n=460). Cross-validation was performed using 5 independent bladder cancer patient cohorts (n=1006). We also experimentally validated the changes of gene expression levels representing each subgroup. FINDINGS: The 976-gene signature based on YAP1-activation redefined three subgroups and had the benefits of Bacillus Calmette-Guérin (BCG) treatment in patients with NMIBC (hazard ratio 3.32, 95% CI 1.29-8.56, p = 0.01). The integrated analysis revealed that YAP1 activation was associated with the characterization of patients with high-risk NMIBC and the response to immunotherapy. INTERPRETATION: This study suggests that YAP1 activation has an important prognostic effect on bladder cancer progression and might be useful in the selection of immunotherapy. FUNDING: A funding list that contributed to this research can be found in the Acknowledgements section.

Transcriptomic Analysis of Papillary Thyroid Cancer: A Focus on Immune-Subtyping, Oncogenic Fusion, and Recurrence.

OBJECTIVES: Thyroid cancer is the most common endocrine tumor, with rapidly increasing incidence worldwide. However, its transcriptomic characteristics associated with immunological signatures, driver fusions, and recurrence markers remain unclear. We aimed to investigate the transcriptomic characteristics of advanced papillary thyroid cancer. METHODS: This study included 282 papillary thyroid cancer tumor samples and 155 normal samples from Chungnam National University Hospital and Seoul National University Hospital. Transcriptomic quantification was determined by high-throughput RNA sequencing. We investigated the associations of clinical parameters and molecular signatures using RNA sequencing. We validated predictive biomarkers using the Cancer Genome Atlas database. RESULTS: Through a comparison of differentially expressed genes, gene sets, and pathways in papillary thyroid cancer compared to normal tumor-adjacent tissue, we found increased immune signaling associated with cytokines or T cells and decreased thyroid hormone synthetic pathways. In addition, patients with recurrence presented increased CD8+ T-cell and Th1-cell signatures. Interestingly, we found differentially overexpressed genes related to immune-escape signaling such as CTLA4, IDO1, LAG3, and PDCD1 in advanced papillary thyroid cancer with a low thyroid differentiation score. Fusion analysis showed that the PI3K and mitogen-activated protein kinase (MAPK) signaling pathways were regulated differently according to the RET fusion partner genes (CCDC6 or NCOA4). Finally, we identified HOXD9 as a novel molecular biomarker that predicts the recurrence of thyroid cancer in addition to known risk factors (tumor size, lymph node metastasis, and extrathyroidal extension). CONCLUSION: We identified a high association with immune-escape signaling in the immune-hot group with aggressive clinical characteristics among Korean thyroid cancer patients. Moreover, RET fusion differentially regulated PI3K and MAPK signaling depending on the partner gene of RET, and HOXD9 was found to be a recurrence marker for advanced papillary thyroid cancer.

Analysis of differential proteomes of induced pluripotent stem cells by protein-based reprogramming of fibroblasts.

The recent generation of induced pluripotent stem (iPS) cells represents a novel opportunity to complement embryonic stem (ES) cell-based approaches. iPS cells can be generated by viral transduction of specific transcription factors, but there is a potential risk of tumorigenicity by random retroviral integration. We have generated novel iPS (sFB-protein-iPS) cells from murine dermal fibroblasts (FVB-sFB) that have ES cell characteristics, using ES cell-derived cell extracts instead of performing viral transduction. Notably, only cell extracts from an ES cell line (C57-mES) on the C57/BL6 background generated iPS cells in our protocol-not an ES cell line (E14-mES) on the 129 background. Hypothesizing that determining the differences in these 2 mES cell lines will provide vital insight into the reprogramming machinery, we performed proteomic and global gene expression analysis by iTRAQ and mRNA microarray, respectively. We observed that pluripotent ES cells and ES cell extract-derived iPS cells had differential proteomes and global gene expression patterns. Notably, reprogramming-competent C57-mES cells highly expressed proteins that regulate protein synthesis and metabolism, compared with reprogramming-incompetent 129-mES cells, suggesting that there is a threshold that protein synthetic machinery must exceed to initiate reprogramming.