Gene expression profiling of alpha-radiation-induced rat osteosarcomas: Identification of dysregulated genes involved in radiation-induced tumorigenesis of bone

Effect of radiation exposure on survival after first solid cancer diagnosis in A-bomb survivors

BACKGROUND

Comparison of the estimated effect of atomic bomb radiation exposure on solid cancer incidence and solid cancer mortality in the RERF Life Span Study (LSS) reveals a difference in the magnitude and shape of the excess relative risk dose response. A possible contributing factor to this difference is pre-diagnosis radiation effect on post-diagnosis survival. Pre-diagnosis radiation exposure theoretically could influence post-diagnosis survival by affecting the genetic makeup and possibly aggressiveness of cancer, or by compromising tolerance for aggressive treatment for cancer.

METHODS

We analyze the radiation effect on post-diagnosis survival in 20,463 LSS subjects diagnosed with first-primary solid cancer between 1958 and 2009 with particular attention to whether death was caused by the first-primary cancer, other cancer, or non-cancer diseases.

RESULTS

From multivariable Cox regression analysis of cause-specific survival, the excess hazard at 1 Gy (EH_1Gy) for death from the first primary cancer was not significantly different from zero - p = 0.23, EH_1Gy = 0.038 (95 % CI: -0.023, 0.104). Death from other cancer and death from non-cancer diseases both were significantly associated with radiation dose: other cancer EH_1Gy = 0.38 (95 % CI: 0.24, 0.53); non-cancer EH_1Gy = 0.24 (95 % CI: 0.13, 0.36), both p < 0.001.

CONCLUSION

There is no detectable large effect of pre-diagnosis radiation exposure on post-diagnosis death from the first primary cancer in A-bomb survivors.

IMPACT

A direct effect of pre-diagnosis radiation exposure on cancer prognosis is ruled out as an explanation for the difference in incidence and mortality dose response in A-bomb survivors.

Pan-cancer analyses reveal IGSF10 as an immunological and prognostic biomarker

Background: IGSF10 is a member of the immunoglobulin superfamily. Over the previous decade, growing proof has validated definitive correlations between individuals of the immunoglobulin superfamily and human diseases. However, the function of IGSF10 in pan-cancer stays unclear. We aimed to analyze the immunological and prognostic value of IGSF10 in pan-cancer. Methods: We utilized a vary of bioinformatic ways to inspect the function of IGSF10 in pan-cancer, including its correlation with prognosis, immune cell infiltration, tumor mutational burden (TMB), microsatellite instability (MSI), mismatch repair (MMR), DNA methyltransferases, genetic alteration, drug sensitivity, etc. Results: We noticed low expression of IGSF10 in most cancer types. IGSF10 expression in tumor samples correlates with prognosis in most cancers. In most cancer types, IGSF10 expression was strongly related to immune cells infiltration, immune checkpoints, immune modulators, TMB, MSI, MMR, and DNA methyltransferases, among others. Functional enrichment analyses indicated that IGSF10 expression was involved in lymphocyte differentiation, cell molecules adhesion, etc. Furthermore, low IGSF10 expression could increase the drug sensitivity of many drugs. Conclusion: IGSF10 could serve as a novel prognostic marker and attainable immunotherapy target for several malignancies.

MicroRNA-106b-5p inhibits growth and progression of lung adenocarcinoma cells by downregulating IGSF10

In this study, we investigated the mechanistic role and prognostic significance of IGSF10 in lung adenocarcinoma. Oncomine database analysis showed that IGSF10 expression was significantly reduced in most cancer types, including lung adenocarcinoma (LUAD). In the TCGA-LUAD dataset, IGSF10 expression correlated positively with proportions of tumor-infiltrated B cells, CD4⁺ T cells, CD8⁺ T cells, neutrophils, macrophages, and dendritic cells. Kaplan-Meier survival analysis showed that overall survival of patients with low IGSF10 expression was significantly shorter than those with high IGSF10 expression. MiRWalk2.0 database analysis and dual luciferase reporter assays confirmed that miR-106b-5p suppressed IGSF10 expression by binding to its 3’UTR. MiR-106b-5p levels inversely correlated with IGSF10 expression in the TCGA-LUAD dataset. Moreover, inhibition of miR-106b-5p significantly decreased in vitro proliferation, migration, and invasion by LUAD cells, whereas miR-106b-5p overexpression reversed those effects. These results demonstrate that IGSF10 is an independent prognostic factor for LUAD. Furthermore, miR-106b-5p suppressed IGSF10 expression in LUAD tissues by binding to its 3’UTR, which makes IGSF10 and miR-106b-5p potential prognostic biomarkers and therapeutic targets in LUAD patients.

Comprehensive Analysis of Aerobic Exercise-Related Genes Identifies CDCA4 That Promotes the Progression of Osteosarcoma

Background

Exercise has a positive impact on patients with osteosarcoma, improving function, reducing disability, maintaining independence and quality of life. Exercise may also directly affect the effectiveness of cancer treatment. Cell division cycle-associated protein 4 (CDCA4) is reported to function importantly during numerous human cancers development. Nevertheless, the details toward CDCA4 function are still to be investigated.

Methods

This study comprehensively analyzed the GSE74194 database and obtained aerobic exercise-related genes. Protein-protein interaction network (PPI) and Gene Ontology (GO) analysis were performed on the differentially expressed genes (DEGs). Quantitative reverse transcription polymerase chain reaction (qRT-PCR) and tumor genome atlas (TCGA) data mining were applied to measure aerobic exercise-related gene CDCA4 level in osteosarcoma tissue. We conducted lots of functional experiments to uncover CDCA4 function and its corresponding mechanism in osteosarcoma.

Results

We screened a total of 547 DEGs related to aerobic exercise, of which 373 were up-regulated and 174 were down-regulated. PPI analysis revealed 90 genes that might play key roles. GO analysis showed that aerobic exercise-related DEGs were significantly enriched during the mitotic cell cycle, cell division, mitotic nuclear division and sister chromatid segregation, nuclear division, microtubule cytoskeleton organization involved protein, microtubule-based process, spindle organization, G2/M transition of mitotic cell cycle. Our results indicated that CDCA4 was increased in osteosarcoma tissues and cell lines, and its level had association with high mortality of osteosarcoma patients. Further studies revealed that absence of CDCA4 largely hindered osteosarcoma cancer cell proliferation, invasion, and migration.

Conclusion

Comprehensive bioinformatics analysis improves our understanding of the underlying molecular mechanisms of aerobic exercise on osteosarcoma. This provides evidence for the effect of aerobic exercise on CDCA4 expression. Our data suggested that CDCA4 could facilitate osteosarcoma development, and gave a hint that CDCA4 was a candidate target in the treatment of osteosarcoma, aerobic exercise might help the treatment and prognosis of patients with osteosarcoma.

Immunoglobulin superfamily member 10 is a novel prognostic biomarker for breast cancer

Background

Immunoglobulin superfamily member 10 (IGSF10) is a member of the immunoglobulin superfamily that is expressed at high levels in both the gallbladder and ovary. Currently, the role and possible mechanism of IGSF10 in breast cancer remain unclear.

Method

By applying real-time quantitative polymerase chain reaction (qRT-PCR) and immunohistochemistry (IHC), the expression of IGSF10 in breast cancer cells and tissues was detected. We collected the clinical information from 700 patients with breast cancer in The Cancer Genome Atlas (TCGA), and analyzed the relationship between IGSF10 expression and the clinicopathological features and survival outcomes of these patients. The potential mechanisms and pathways associated with IGSF10 in breast cancer were explored by performing a gene set enrichment analysis (GSEA).

Results

According to TCGA data, qRT-PCR and IHC experiments, levels of the IGSF10 mRNA and protein were significantly decreased in breast cancer tissues. IGSF10 expression was significantly correlated with age, tumor size, and tumor stage. Moreover, shorter overall survival (OS) and relapse-free survival (RFS) correlated with lower IGSF10 expression, according to the survival analysis. The multivariate analysis identified that IGSF10 as an independent prognostic factor for the OS (hazard ratio (HR) = 1.793, 95% confidence interval (CI) [1.141–2.815], P = 0.011) and RFS (HR = 2.298, 95% CI [1.317–4.010], P = 0.003) of patients with breast cancer. Based on the GSEA, IGSF10 was involved in DNA repair, cell cycle, and glycolysis. IGSF10 was also associated with the PI3K/Akt/mTOR and mTORC1 signaling pathways.

Conclusions

This study revealed a clear relationship between IGSF10 expression and the tumorigenesis of breast cancer for the first time. Therefore, further studies are needed to understand the mechanism of IGSF10 in breast cancer.

Identification of Metastasis-Associated Genes in Triple-Negative Breast Cancer Using Weighted Gene Co-expression Network Analysis

Triple-negative breast cancer (TNBC) is the most aggressive and fatal sub-type of breast cancer. This study aimed to identify metastasis-associated genes that could serve as biomarkers for TNBC diagnosis and prognosis. RNA-seq data and clinical information on TNBC from the Cancer Genome Atlas were used to conduct analyses. Expression data were used to establish co-expression modules using average linkage hierarchical clustering. We used weighted gene co-expression network analysis to explore the associations between gene sets and clinical features and to identify metastasis-associated candidate biomarkers. The K-M plotter website was used to explore the association between the expression of candidate biomarkers and patient survival. In addition, receiver operating characteristic curve analysis was used to illustrate the diagnostic performance of candidate genes. The pale turquoise module was significantly associated with the occurrence of metastasis. In this module, 64 genes were identified, and its functional enrichment analysis revealed that they were mainly associated with transcriptional misregulation in cancer, microRNAs in cancer, and negative regulation of angiogenesis. Further, 4 genes, IGSF10, RUNX1T1, XIST, and TSHZ2, which were negatively associated with relapse-free survival and have seldom been reported before in TNBC, were selected. In addition, the mRNA expression levels of the 4 candidate genes were significantly lower in TNBC tumor tissues compared with healthy tissues. Based on the K-M plotter, these 4 genes were correlated with poor prognosis of TNBC. The area under the curve of IGSF10, RUNX1T1, TSHZ2, and XIST was 0.918, 0.957, 0.977, and 0.749. These findings provide new insight into TNBC metastasis. IGSF10, RUNX1T1, TSHZ2, and XIST could be used as candidate biomarkers for the diagnosis and prognosis of TNBC metastasis.

Screening of Long Noncoding RNAs Induced by Radiation Using Microarray

DNA damage repair and G2/M arrest are the key factors regulating the survival of cancer cells exposed to radiation. Recent studies have shown that long noncoding RNAs (lncRNAs) play important roles in a variety of biological processes, including DNA repair, cell cycle regulation, differentiation, and epigenetic regulation. However, the knowledge about the genome scale of lncRNAs and their potential biological functions in tumor cells exposed to radiation are still unclear. In this study, we used LncRNA + mRNA Human Gene Expression Microarray V4.0 to profile lncRNA and messenger RNA (mRNA) from HeLa, MCF-7, and A549 cells after irradiation with 4 Gy of γ-radiation. We identified 230, 227, and 274 differentially expressed lncRNAs and 150, 214, and 274 differentially expressed mRNAs in HeLa, MCF-7, and A549 cells, respectively, among which there are 14 common differentially expressed lncRNAs and 22 common differentially expressed mRNAs in all of the 3 cell lines. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis indicated that these differentially expressed mRNAs were mainly associated with cell cycle. Further, we also predicted the target genes and functions of these differentially expressed lncRNAs. Our study on lncRNAs has greatly expanded the field of gene research in the relationship of radiation, cell cycle, and DNA damage.

Identification of prognostic markers of lung cancer through bioinformatics analysis and in vitro experiments

Lung cancer is one of the most common types of cancer worldwide. Understanding the molecular mechanisms underlying the development and progression of lung cancer may improve early diagnosis, treatment and prognosis. The aim of the present study was to examine the pathogenesis of lung cancer and to identify potentially novel biomarkers. Gene expression datasets of patients with lung cancer were obtained from the Gene Expression Omnibus. Genes which were most closely associated with lung cancer (core genes) were screened by weighted gene co‑expression network analysis. In vitro cell based experiments were further utilized to verify the effects of the core genes on the proliferation of lung cancer cells, adhesion between cells and the matrix, and the associated metabolic pathways. Based on WGCNA screening, two gene modules and five core genes closely associated with lung cancer, including immunoglobulin superfamily member 10 (IGSF10) from the turquoise module, and ribonucleotide reductase regulatory subunit M2, protein regulator of cytokinesis 1, kinesin family member (KIF)14 and KIF2C from the brown module were identified as relevant. Survival analysis and differential gene expression analysis showed that there were significant differences in IGSF10 expression levels between the healthy controls and patients with lung cancer. In patients with lung cancer, IGSF10 expression was decreased, and the overall survival time of patients with lung cancer was significantly shortened. An MTT and colony formation assay showed that IGSF10‑knockout significantly increased proliferation of lung cancer cells, and Transwell assays and adhesion experiments further suggested that the adhesion between cells and the matrix was significantly increased in IGSF10‑knockout cells. Gene Set Enrichment Analysis showed that the expression level of IGSF10 was significantly associated with the activation of the integrin‑β1/focal adhesion kinase (FAK) pathway. Western blotting revealed that knockout of IGSF10 resulted in the activation of the integrin‑β1/FAK pathway, as the protein expression levels of integrin‑β1, phosphorylated (p)‑FAK and p‑AKT were significantly upregulated. Activation of the integrin‑β1/FAK pathway, following knockout of IGSF10, affected the proliferation and adhesion of lung cancer cells. Therefore, IGSF10 my serve as a potential prognostic marker of lung cancer.

Leucine zipper and ICAT domain containing (LZIC) protein regulates cell cycle transitions in response to ionizing radiation

Common hallmarks of cancer include the dysregulation of cell cycle progression and the acquisition of genome instability. In tumors, G1 cell cycle checkpoint induction is often lost. This increases the reliance on a functional G2/M checkpoint to prevent progression through mitosis with damaged DNA, avoiding the introduction of potentially aberrant genetic alterations. Treatment of tumors with ionizing radiation (IR) utilizes this dependence on the G2/M checkpoint. Therefore, identification of factors which regulate this process could yield important biomarkers for refining this widely used cancer therapy. Leucine zipper and ICAT domain containing (LZIC) downregulation has been associated with the development of IR-induced tumors. However, despite LZIC being highly conserved, it has no known molecular function. We demonstrate that LZIC knockout (KO) cell lines show a dysregulated G2/M cell cycle checkpoint following IR treatment. In addition, we show that LZIC deficient cells competently activate the G1 and early G2/M checkpoint but fail to maintain the late G2/M checkpoint after IR exposure. Specifically, this defect was found to occur downstream of PIKK signaling. The LZIC KO cells demonstrated severe aneuploidy indicative of genomic instability. In addition, analysis of data from cancer patient databases uncovered a strong correlation between LZIC expression and poor prognosis in several cancers. Our findings suggest that LZIC is functionally involved in cellular response to IR, and its expression level could serve as a biomarker for patient stratification in clinical cancer practice.

Identification of BLCAP as a novel STAT3 interaction partner in bladder cancer

Bladder cancer associated protein (Blcap) expression is commonly down-regulated in invasive bladder cancer, and may have prognostic value given that its expression is negatively correlated with patient survival. We have previously investigated the expression patterns and cellular localization of Blcap in bladder cancer, where we found that about 20% of the lesions examined displayed strong nuclear expression of Blcap, and that this phenotype was associated with overall poor disease outcome. Here we report on the analysis of possible functional associations between nuclear expression of Blcap and canonical signaling pathways. We performed serial immunohistochemistry (IHC) analysis of bladder tissue samples, with serial sections stained with phospho-specific antibodies recognizing key signaling intermediates, such as P-Stat3, P-Akt, and P-Erk1/2, among others, in an immunophenotyping approach we have established and reported previously. Using this approach, we found that nuclear localization of Blcap was associated with expression of P-Stat3. A parallel analysis, cytokine profiling of bladder tumor interstitial fluids of samples expressing (or not) Blcap, showed interleukin (IL)-6, IL-8, and monocyte chemotactic protein 1 (MCP-1) to be correlated with nuclear expression of Blcap, independently supporting a role for Stat3 signaling in localization of Blcap. Multiple indirect immunofluorescence analysis of tissue biopsies confirmed that Blcap co-localized with Stat3. Furthermore, we could also demonstrate, using an in situ proximity ligation assay that Blcap and Stat3 are in close physical proximity of each other in bladder tissue, and that Blcap physically interacts with Stat3 as determined by co-immunoprecipitation of these proteins. Our data indicates that Blcap is a novel Stat3 interaction partner and suggests a role for Blcap in the Stat3-mediated progression of precancerous lesions to invasive tumors of the bladder.

Hypoxia-Inducible Factor 2α Mutation-Related Paragangliomas Classify as Discrete Pseudohypoxic Subcluster

Recently, activating mutations of the hypoxia-inducible factor 2α gene (HIF2A/EPAS1) have been recognized to predispose to multiple paragangliomas (PGLs) and duodenal somatostatinomas associated with polycythemia, and ocular abnormalities. Previously, mutations in the SDHA/B/C/D, SDHAF2, VHL, FH, PHD1, and PHD2 genes have been associated with HIF activation and the development of pseudohypoxic (cluster-1) PGLs. These tumors overlap in terms of tumor location, syndromic presentation, and noradrenergic phenotype to a certain extent. However, they also differ especially by clinical outcome and by presence of other tumors or abnormalities. In the present study, we aimed to establish additional molecular differences between HIF2A and non-HIF2A pseudohypoxic PGLs. RNA expression patterns of HIF2A PGLs (n = 6) from 2 patients were compared with normal adrenal medullas (n = 8) and other hereditary pseudohypoxic PGLs (VHL: n = 13, SDHB: n = 15, and SDHD: n = 14). Unsupervised hierarchical clustering showed that HIF2A PGLs made up a separate cluster from other pseudohypoxic PGLs. Significance analysis of microarray yielded 875 differentially expressed genes between HIF2A and other pseudohypoxic PGLs after normalization to adrenal medulla (false discovery rate 0.01). Prediction analysis of microarray allowed correct classification of all HIF2A samples based on as little as three genes (TRHDE, LRRC63, IGSF10; error rate: 0.02). Genes with the highest expression difference between normal medulla and HIF2A PGLs were selected for confirmatory quantitative reverse transcriptase polymerase chain reaction. In conclusion, HIF2A PGLs show a characteristic expression signature that separates them from non-HIF2A pseudohypoxic PGLs. Unexpectedly, the most significantly differentially expressed genes have not been previously described as HIF target genes.

Soluble CD93 in Serum as a Marker of Allergic Inflammation

PURPOSE

CD93 is receiving renewed attention as a biomarker of inflammation. We aimed to evaluate the potential for serum sCD93 to serve as a novel biomarker for allergic inflammation.

MATERIALS AND METHODS

We enrolled 348 subjects with an allergic disease [allergic rhinitis (AR), chronic spontaneous urticaria (CSU), or bronchial asthma (BA)], including 14 steroid-naïve BA patients who were serially followed-up.

RESULTS

The serum sCD93 levels (ng/mL) in patients with exacerbated AR (mean±standard deviation, 153.1±58.4) were significantly higher than in patients without AR (132.2±49.0) or with stable AR (122.3±42.1). Serum sCD93 levels in exacerbated CSU (169.5±42.8) were also significantly higher than those in non-CSU (132.4±51.6) and stable CSU (122.8±36.2). This trend was also seen in BA. Serum levels in patients with ICS-naïve BA (161.4±53.1) were significantly higher than those in healthy controls without BA (112.2±30.8), low- and medium-dose ICS users. Serum sCD93 levels in high-dose ICS users (72.2±20.6) were significantly lower than those in low- and medium-dose users. The serum sCD93 levels in steroid-naïve patients with BA (195.1±72.7) decreased after ICS use for 4 weeks (134.4±42.8) and 8 weeks (100.7±13.4), serially.

CONCLUSION

Elevated serum sCD93 levels reflected exacerbated status of allergic diseases, including CSU, AR, and asthma. ICS use significantly diminished serum sCD93 levels in steroid-naïve patients with BA. This result may suggest sCD93 in serum as a therapeutic marker for allergic inflammation.

Combined treatment with X-ray irradiation and 5-aminolevulinic acid elicits better transcriptomic response of cell cycle-related factors than X-ray irradiation alone

PURPOSE

5-Aminolevulinic acid (ALA) is a precursor of the photosensitizer protoporphyrin (PpIX) used in photodynamic therapy. In our previous work, PpIX enhanced the generation of reactive oxygen species by X-ray irradiation. In this study, we evaluated the potential of ALA as an endogenous sensitizer to X-ray irradiation.

METHODOLOGY

Tumor-bearing C57BL/6J mice implanted with B16-BL6 melanoma cells were subsequently treated with irradiation (3 Gy/day for 10 days; total, 30 Gy) plus local administration of 50 mg/kg ALA 24 hours prior to each irradiation (ALA-XT). Tumor-bearing mice without treatment (NT), those treated with ALA only (ALAT), and those treated with X-ray irradiation only (XT) were used as controls.

RESULTS

ALA potentiated tumor suppression by X-ray irradiation. In microarray analyses using tumor tissue collected after 10 sessions of fractional irradiation, functional analysis revealed that the majority of dysregulated genes in the XT and ALA-XT groups were related to cell-cycle arrest. Finally, the XT and ALA-XT groups differed in the strength of expression, but not in the pattern of expression.

CONCLUSIONS

mRNA analysis revealed that the combined use of ALA and X-ray irradiation sensitized tumors to X-ray treatment. Furthermore, the present results were consistent with ALA's tumor suppressive effects in vivo.

Comparative analysis of micro-RNAs in human papillomavirus-positive versus -negative oropharyngeal cancers

BACKGROUND

Oncogenic mechanisms of human papillomavirus (HPV)-positive oropharyngeal cancer are still poorly characterized. Analysis of their microRNA expression profile might provide valuable information.

METHODS

The microRNA expression profiles were analyzed by micro-arrays in 26 oropharyngeal cancers. A microRNA signature specific to HPV-status was identified by analyzing a learning/training set consisting of 16 oropharyngeal cancers. The robustness of this signature was further confirmed by blind case-by-case classification of a validation set composed of 10 independent tumors. Putative targeted molecular pathways were proposed using DIANA miRPath online software (http://microrna.gr/mirpath).

RESULTS

We have identified 25 miRNA signatures, which discriminates HPV16-positive oropharyngeal cancer from their HPV-negative counterparts. These 25 microRNAs play a potential role in Wnt and PI3K-pathways, cell-adhesion/cell-polarity, and the cytoskeleton regulation.

CONCLUSION

Our study contributes to a better understanding of pathogenic mechanisms involved in the development of HPV-positive oropharyngeal cancer and in the identification of potential therapeutic molecular targets. © 2016 Wiley Periodicals, Inc. Head Neck 38: 1708-1716, 2016.

Exome sequencing in one family with gastric- and rectal cancer

Background

Heritable factors are well known to increase the risk of cancer in families. Known susceptibility genes account for a small proportion of all colorectal cancer cases. The aim of this study was to identify the genetic background in a family suggested to segregate a dominant cancer syndrome with a high risk of rectal- and gastric cancer. We performed whole exome sequencing in three family members, 2 with rectal cancer and 1 with gastric cancer and followed it up in additional family members, other patients and controls.

Results

We identified 12 novel non-synonymous single nucleotide variants, which were shared among 5 affected members of this family. The mutations were found in 12 different genes; DZIP1L, PCOLCE2, IGSF10, SUCNR1, OR13C8, EPB41L4B, SEC16A, NOTCH1, TAS2R7, SF3A1, GAL3ST1, and TRIOBP. None of the mutations was suggested as a high penetrant mutation. It was not possible to completely rule out any of the mutations as contributing to disease, although seven were more unlikely than the others. Neither did we rule out the effect of all thousands of intronic, intergenic and synonymous variants shared between the three persons used for exome sequencing.

Conclusions

We propose this family, suggested to segregate dominant disease, could be an example of complex inheritance.

Downregulation of tumor suppressing STF cDNA 3 promotes epithelial-mesenchymal transition and tumor metastasis of osteosarcoma by the Wnt/GSK-3β/β-catenin/Snail signaling pathway

Epithelial to mesenchymal transition (EMT) has received considerable attention as a conceptual paradigm for explaining the invasive and metastatic behavior of cells during cancer progression. Our previous study showed that loss of expression of TSSC3 is positively associated with osteosarcoma malignancy and progression. However, whether TSSC3 mediates EMT in osteosarcoma is poorly understood. In the present study, we determined that TSSC3 downregulation induced cell migration and invasion ability and promoted mesenchymal transition of osteosarcoma cells by upregulating mesenchymal markers and inhibiting the epithelial markers. Furthermore, TSSC3 downregulation elicited a signaling cascade that included increased levels of Wnt3a and LRP5, inactivation of GSK-3β, accumulation of nuclear β-catenin and Snail, the augmented binding of β-catenin to TCF-4, and accordingly increased the expression of Wnt target genes (CD44, MMP7). The gene knockdown of these signaling proteins could inhibit TSSC3 downregulation-promoted EMT, migration, and invasion in osteosarcoma. Finally, TSSC3 overexpression obviously inhibited cell migration, invasion, and repressed mesenchymal phenotypes, reducing lung metastasis through GSK-3β activation. Collectively, TSSC3 downregulation promotes the EMT of osteosarcoma cells by regulating EMT markers via a signal transduction pathway that involves Snail, Wnt-β-catenin/TCF, and GSK-3β.

Effects of ionizing radiation on differentiation of murine bone marrow cells into mast cells

Mast cells, immune effector cells produced from bone marrow cells, play a major role in immunoglobulin E-mediated allergic responses. Ionizing radiation affects the functions of mast cells, which are involved in radiation-induced tissue damage. However, whether ionizing radiation affects the differential induction of mast cells is unknown. Here we investigated whether bone marrow cells of X-irradiated mice differentiated into mast cells. To induce mast cells, bone marrow cells from X-irradiated and unirradiated mice were cultured in the presence of cytokines required for mast cell induction. Although irradiation at 0.5 Gy and 2 Gy decreased the number of bone marrow cells 1 day post-irradiation, the cultured bone marrow cells of X-irradiated and unirradiated mice both expressed mast cell-related cell-surface antigens. However, the percentage of mast cells in the irradiated group was lower than in the unirradiated group. Similar decreases in the percentage of mast cells induced in the presence of X-irradiation were observed 10 days post irradiation, although the number of bone marrow cells in irradiated mice had recovered by this time. Analysis of mast cell function showed that degranulation of mast cells after immunoglobulin E-mediated allergen recognition was significantly higher in the X-irradiated group compared with in the unirradiated group. In conclusion, bone marrow cells of X-irradiated mice differentiated into mast cells, but ionizing radiation affected the differentiation efficiency and function of mast cells.

A predictive transcriptomic signature of oropharyngeal cancer according to HPV16 status exclusively

OBJECTIVE

Human-papillomaviruses (HPV) type 16 is a causative agent in an increasing subset of oropharyngeal squamous cell carcinomas (OPSCCs). These tumors have distinct oncogenic mechanisms and a more favorable prognosis than tobacco-induced OPSCCs. Although these differences emphasize the need for a specific therapeutic approach, HPV status is still not used to guide treatment. A better characterization of the molecular profile related to HPV16-induced OPSCC might help to develop personalized treatments.

PATIENTS AND METHODS

Using a human whole-genome DNA-microarray, we have examined the gene expression profiles in 15 HPV-negative and 15 transcriptionally-active HPV-positive OPSCCs. The study was conducted in two steps. Firstly, a learning/training-set consisting of 8 HPV16-positive and 8 HPV16-negative OPSCCs was analyzed to identify a specific signature. Potentially confounding factors (stage, sex and tobacco) were equally distributed in both groups. Subsequently the robustness of this signature was confirmed by blind case-by-case classification of a validation-set composed of the 14 remaining tumors.

RESULTS

We have identified a signature composed of 224 genes, which discriminates HPV16-induced OPSCC from their HPV-negative counterparts. After the blind classification of the 14 tumours, the viral status was revealed: 13 out of 14 tumors were correctly classified according to tumor etiology, 1/14 was not determined and none were misclassified. Several of the differentially expressed genes were involved in cell-cycle regulation, DNA replication and repair, transcription regulation, immune response and apoptosis.

CONCLUSION

Our study contributes to a better understanding of pathogenic mechanisms involved in the development of HPV-positive OPSCCs and in the identification of potential therapeutic targets.

The genetic basis for inactivation of Wnt pathway in human osteosarcoma

BACKGROUND

Osteosarcoma is a highly genetically unstable tumor with poor prognosis. We performed microarray-based comparative genomic hybridization (aCGH), transcriptome sequencing (RNA-seq), and pathway analysis to gain a systemic view of the pathway alterations of osteosarcoma.

METHODS

aCGH experiments were carried out on 10 fresh osteosarcoma samples. The output data (Gene Expression Omnibus Series accession number GSE19180) were pooled with published aCGH raw data (GSE9654) to determine recurrent copy number changes. These were analyzed using Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis to identify altered pathways in osteosarcoma. Transcriptome sequencing of six osteosarcomas was performed to detect the expression profile of Wnt signaling pathway genes. Protein expression of WNT1, β-catenin, c-myc, and cyclin D1 in the Wnt pathway was detected by immunohistochemistry (IHC) in an independent group of 46 osteosarcoma samples.

RESULTS

KEGG pathway analysis identified frequent deletions of Wnt and other Wnt signaling pathway genes. At the mRNA level, transcriptome sequencing found reduced levels of mRNA expression of Wnt signaling pathway transcripts. While WNT1 protein expression was detected by IHC in 69.6% (32/46) of the osteosarcomas, no β-catenin protein was detected in the nucleus. β-catenin protein expression was, however, detected in the membrane and cytoplasm of 69.6% (32/46) of the osteosarcomas. c-myc protein expression was detected in only 47.8% (22/46) and cyclin D1 protein expression in 52.2% (24/46) of osteosarcoma samples. Kaplan-Meier survival analysis showed that WNT1-negative patients had a trend towards longer disease free survival than WNT1-positive patients. Interestingly, in WNT1-negative patients, those who were also cyclin D1-negative had significantly longer disease free survival than cyclin D1-positive patients. However, there was no significant association between any of the investigated proteins and overall survival of human osteosarcoma patients.

CONCLUSIONS

Frequent deletions of Wnt and other Wnt signaling pathway genes suggest that the Wnt signaling pathway is genetically inactivated in human osteosarcoma.

Modeling distinct osteosarcoma subtypes in vivo using Cre:lox and lineage-restricted transgenic shRNA

Osteosarcoma is the most common primary cancer of bone and one that predominantly affects children and adolescents. Osteoblastic osteosarcoma represents the major subtype of this tumor, with approximately equal representation of fibroblastic and chondroblastic subtypes. We and others have previously described murine models of osteosarcoma based on osteoblast-restricted Cre:lox deletion of Trp53 (p53) and Rb1 (Rb), resulting in a phenotype most similar to fibroblastic osteosarcoma in humans. We now report a model of the most prevalent form of human osteosarcoma, the osteoblastic subtype. In contrast to other osteosarcoma models that have used Cre:lox mediated gene deletion, this model was generated through shRNA-based knockdown of p53. As is the case with the human disease the shRNA tumors most frequently present in the long bones and preferentially disseminate to the lungs; feature less consistently modeled using Cre:lox approaches. Our approach allowed direct comparison of the in vivo consequences of targeting the same genetic drivers using two different technologies, Cre:lox and shRNA. This demonstrated that the effects of Cre:lox and shRNA mediated knock-down are qualitatively different, at least in the context of osteosarcoma, and yielded distinct subtypes of osteosarcoma. Through the use of complementary genetic modification strategies we have established a model of the most common clinical subtype of osteosarcoma that was not previously represented and more fully recapitulated the clinical spectrum of this cancer.

Hck contributes to bone homeostasis by controlling the recruitment of osteoclast precursors

In osteoclasts, Src controls podosome organization and bone degradation, which leads to an osteopetrotic phenotype in src(-/-) mice. Since this phenotype was even more severe in src(-/-)hck(-/-) mice, we examined the individual contribution of Hck in bone homeostasis. Compared to wt mice, hck(-/-) mice exhibited an osteopetrotic phenotype characterized by an increased density of trabecular bone and decreased bone degradation, although osteoclastogenesis was not impaired. Podosome organization and matrix degradation were found to be defective in hck(-/-) osteoclast precursors (preosteoclast) but were normal in mature hck(-/-) osteoclasts, probably through compensation by Src, which was specifically overexpressed in mature osteoclasts. As a consequence of podosome defects, the 3-dimensional migration of hck(-/-) preosteoclasts was strongly affected in vitro. In vivo, this translated by altered bone homing of preosteoclasts in hck(-/-) mice: in metatarsals of 1-wk-old mice, when bone formation strongly depends on the recruitment of these cells, reduced numbers of osteoclasts and abnormal developing trabecular bone were observed. This phenotype was still detectable in adults. In summmary, Hck is one of the very few effectors of preosteoclast recruitment described to date and thereby plays a critical role in bone remodeling.

Immunoexpression analysis and prognostic value of BLCAP in breast cancer

Bladder Cancer Associated Protein (BLCAP, formerly Bc10), was identified by our laboratory as being down-regulated in bladder cancer with progression. BLCAP is ubiquitously expressed in different tissues, and several studies have found differential expression of BLCAP in various cancer types, such as cervical and renal cancer, as well as human tongue carcinoma and osteosarcoma. Here we report the first study of the expression patterns of BLCAP in breast tissue. We analyzed by immunohistochemistry tissue sections of normal and malignant specimens collected from 123 clinical high-risk breast cancer patients within the Danish Center for Translational Breast Cancer Research (DCTB) prospective study dataset. The staining pattern, the distribution of the immunostaining, and its intensity were studied in detail. We observed weak immunoreactivity for BLCAP in mammary epithelial cells, almost exclusively localizing to the cytoplasm and found that levels of expression of BLCAP were generally higher in malignant cells as compared to normal cells. Quantitative IHC analysis of BLCAP expression in breast tissues confirmed this differential BLCAP expression in tumor cells, and we could establish, in a 62-patient sample matched cohort, that immunostaining intensity for BLCAP was increased in tumors relative to normal tissue, in more than 45% of the cases examined, indicating that BLCAP may be of value as a marker for breast cancer. We also analyzed BLCAP expression and prognostic value using a set of tissue microarrays comprising an independent cohort of 2,197 breast cancer patients for which we had follow-up clinical information.

Emerging Roles for SSeCKS/Gravin/AKAP12 in the Control of Cell Proliferation, Cancer Malignancy, and Barriergenesis

Emerging data suggest that SSeCKS/Gravin/AKAP12 ("AKAP12"), originally identified as an autoantigen in cases of myasthenia gravis, controls multiple biological processes through its ability to scaffold key signaling proteins such as protein kinase (PK) C and A, calmodulin, cyclins, phosphoinositides, "long" β-1,4 galactosyltransferase (GalTase) isoform, Src, as well as the actin cytoskeleton in a spatiotemporal manner. Specialized functions attributed to AKAP12 include the suppression of cancer malignancy, especially aspects of metastatic progression, regulation of blood-brain and blood-retina barrier formation, and resensitization of β2-adrenergic pain receptors. Recent data identify a direct role for AKAP12 in cytokinesis completion, further suggesting a function as a negative regulator of cell senescence. The current review will discuss the emerging knowledge base of AKAP12-related biological roles and how the factors that affect AKAP12 expression or that interact with AKAP12 at the protein level control cancer progression and blood-tissue barrier formation.

Strategy to find molecular signatures in a small series of rare cancers: validation for radiation-induced breast and thyroid tumors

Methods of classification using transcriptome analysis for case-by-case tumor diagnosis could be limited by tumor heterogeneity and masked information in the gene expression profiles, especially as the number of tumors is small. We propose a new strategy, EMts_2PCA, based on: 1) The identification of a gene expression signature with a great potential for discriminating subgroups of tumors (EMts stage), which includes: a) a learning step, based on an expectation-maximization (EM) algorithm, to select sets of candidate genes whose expressions discriminate two subgroups, b) a training step to select from the sets of candidate genes those with the highest potential to classify training tumors, c) the compilation of genes selected during the training step, and standardization of their levels of expression to finalize the signature. 2) The predictive classification of independent prospective tumors, according to the two subgroups of interest, by the definition of a validation space based on a two-step principal component analysis (2PCA). The present method was evaluated by classifying three series of tumors and its robustness, in terms of tumor clustering and prediction, was further compared with that of three classification methods (Gene expression bar code, Top-scoring pair(s) and a PCA-based method). Results showed that EMts_2PCA was very efficient in tumor classification and prediction, with scores always better that those obtained by the most common methods of tumor clustering. Specifically, EMts_2PCA permitted identification of highly discriminating molecular signatures to differentiate post-Chernobyl thyroid or post-radiotherapy breast tumors from their sporadic counterparts that were previously unsuccessfully classified or classified with errors.

Influence of β-catenin small interfering RNA on human osteosarcoma cells

This study examined the effect of small interfering RNA-mediated β-catenin knockdown on the survival, invasion and chemosensitivity of human osteosarcoma cells (U2-OS cells). The siRNA against β-catenin was constructed and transfected into U2-OS cells. The expression of β-catenin was detected by qRT-PCR and Western blotting. Cell growth and apoptosis was detected in the presence or absence of doxorubicin by MTT and flow cytometry, respectively. Cell invasion ability was measured by transwell assay. The results showed that the transfection of β-catenin siRNA resulted in decreased expression of β-catenin, suppression of invasion and motility of U2-OS cells, reduced chemosensitivity to doxorubicin in vitro, and little change in cell growth and apoptosis. Additionally, down-regulated MT1-MMP expression was found after transfection. It was concluded that knockdown of β-catenin gene may decrease the invasive ability of human osteosarcoma cells through down-regulated MT1-MMP expression, and the chemosensitivity of osteosarcoma cells against doxorubicin.

Gene expression signature discriminates sporadic from post-radiotherapy-induced thyroid tumors

Both external and internal exposure to ionizing radiation are strong risk factors for the development of thyroid tumors. Until now, the diagnosis of radiation-induced thyroid tumors has been deduced from a network of arguments taken together with the individual history of radiation exposure. Neither the histological features nor the genetic alterations observed in these tumors have been shown to be specific fingerprints of an exposure to radiation. The aim of our work is to define ionizing radiation-related molecular specificities in a series of secondary thyroid tumors developed in the radiation field of patients treated by radiotherapy. To identify molecular markers that could represent a radiation-induction signature, we compared 25K microarray transcriptome profiles of a learning set of 28 thyroid tumors, which comprised 14 follicular thyroid adenomas (FTA) and 14 papillary thyroid carcinomas (PTC), either sporadic or consecutive to external radiotherapy in childhood. We identified a signature composed of 322 genes which discriminates radiation-induced tumors (FTA and PTC) from their sporadic counterparts. The robustness of this signature was further confirmed by blind case-by-case classification of an independent set of 29 tumors (16 FTA and 13 PTC). After the histology code break by the clinicians, 26/29 tumors were well classified regarding tumor etiology, 1 was undetermined, and 2 were misclassified. Our results help shed light on radiation-induced thyroid carcinogenesis, since specific molecular pathways are deregulated in radiation-induced tumors.

Wnts, bone and cancer

The Wnt pathway plays vital roles in bone and in cancer. In this issue of The Journal of Pathology, Cai and colleagues report results that suggest that the Wnt pathway is inactivated in bone cancers, a finding that could have significant implications for the development of Wnt agonists as bone anabolic agents. While these findings are at odds with the prevailing view that the Wnt pathway is oncogenic in all systems studied to date, they remind us how complex and still poorly understood this important signalling pathway remains. At the very least, these findings should provoke debate and stimulate further research into the role of Wnt signalling in osteosarcoma.

Cytogenetic and molecular characterization of plutonium-induced rat osteosarcomas

The association between ionizing radiation and the subsequent development of osteosarcoma has been well described, but little is known about the cytogenetic and molecular events, which could be involved in the formation of radiation-induced osteosarcomas. Here, we performed comparative genomic hybridization (CGH) to detect chromosomal copy number changes in a series of 16 rat osteosarcomas induced by injection of plutonium-238. Recurrent gains/amplifications were observed at chromosomal regions 3p12-q12, 3q41-qter, 4q41-qter, 6q12-q16, 7q22-q34, 8q11-q23, 9q11-q22, 10q32.1-qter, and 12q, whereas recurrent losses were observed at 1p, 1q, 3q23-q35, 5q21-q33, 8q24-q31, 10q22-q25, 15p, 15q, and 18q. The gained region at 7q22-q34 was homologous to human chromosome bands 12q13-q15/8q24/22q11-q13, including the loci of Mdm2, Cdk4, c-Myc and Pdgf-b genes. The lost regions at 5q21-q33, 10q22-q25 and 15q contained tumor suppressor genes such as p16INK4a/p19ARF, Tp53 and Rb1. To identify potential target gene(s) for the chromosomal aberrations, we compared the expression levels of several candidate genes, located within the regions of frequent chromosomal aberrations, between the tumors and normal osteoblasts by using quantitative RT-PCR analysis. The Cdk4, c-Myc, Pdgf-b and p57KIP2 genes were thought to be possible target genes for the frequent chromosomal gain at 7q22-34 and loss at 1q in the tumors, respectively. In addition, mutations of the Tp53 gene were found in 27% (4 of 15) osteosarcomas. Our data may contribute to further understanding of the molecular mechanisms underlying osteosarcomas induced by ionizing radiation in human.

Silencing of Cited2 and Akap12 genes in radiation-induced rat osteosarcomas

We have previously studied genomic copy number changes and global gene expression patterns in rat osteosarcomas (OS) induced by the bone-seeking alpha emitter (238)Pu by comparative genomic hybridization (CGH) and oligonucleotide microarray analyses, respectively. Among the previously identified genes that were down-regulated in radiation-induced rat OS tumors, Cited2 (Cbp/p300-interacting transactivator, with Glu/Asp-rich carboxy-terminal domain, 2) and Akap12 (a kinase anchoring protein, also known as src-suppressed C-kinase substrate, SSeCKS) genes mapped to the most frequently lost regions on chromosome 1p. In the present study, relative copy number losses of Cited2 and Akap12 genes were observed in 8 of 15 (53%) and 10 of 15 (67%) tumors by quantitative PCR analysis. Loss of Cited2 and Akap12 in the tumors was confirmed at the levels of mRNA and protein expression by quantitative RT-PCR and immunoblot analyses, respectively. These results indicate that Cited2 and Akap12 are silenced in radiation-induced OS, and therefore are novel candidate tumor-suppressor genes of this tumor.