Discovery of novel DNA-damaging agents through phenotypic screening for DNA double-strand break.

DNA double-strand breaks (DSBs) seriously damage DNA and promote genomic instability that can lead to cell death. They are the source of conditions such as carcinogenesis and aging, but also have important applications in cancer therapy. Therefore, rapid detection and quantification of DSBs in cells are necessary for identifying carcinogenic and anticancer factors. In this study, we detected DSBs using a flow cytometry-based high-throughput method to analyze γH2AX intensity. We screened a chemical library containing 9600 compounds and detected multiple DNA-damaging compounds, although we could not identify mechanisms of action through this procedure. Thus, we also profiled a representative compound with the highest DSB potential, DNA-damaging agent-1 (DDA-1), using a bioinformatics-based method we termed "molecular profiling." Prediction and verification analysis revealed DDA-1 as a potential inhibitor of topoisomerase IIα, different from known inhibitors such as etoposide and doxorubicin. Additional investigation of DDA-1 analogs and xenograft models suggested that DDA-1 is a potential anticancer drug. In conclusion, our findings established that combining high-throughput DSB detection and molecular profiling to undertake phenotypic analysis is a viable method for efficient identification of novel DNA-damaging compounds for clinical applications.

Functional evaluation of BRCA1/2 variants of unknown significance with homologous recombination assay and integrative in silico prediction model.

Numerous variants of unknown significance (VUSs) exist in hereditary breast and ovarian cancers. Although multiple methods have been developed to assess the significance of BRCA1/2 variants, functional discrepancies among these approaches remain. Therefore, a comprehensive functional evaluation system for these variants should be established. We performed conventional homologous recombination (HR) assays for 50 BRCA1 and 108 BRCA2 VUSs and complementarily predicted VUSs using a statistical logistic regression prediction model that integrated six in silico functional prediction tools. BRCA1/2 VUSs were classified according to the results of the integrative in vitro and in silico analyses. Using HR assays, we identified 10 BRCA1 and 4 BRCA2 VUSs as low-functional pathogenic variants. For in silico prediction, the statistical prediction model showed high accuracy for both BRCA1 and BRCA2 compared with each in silico prediction tool individually and predicted nine BRCA1 and seven BRCA2 variants to be pathogenic. Integrative functional evaluation in this study and the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG/AMP) guidelines strongly suggested that seven BRCA1 variants (p.Glu272Gly, p.Lys1095Glu, p.Val1653Leu, p.Thr1681Pro, p.Phe1761Val, p.Thr1773Ile, and p.Gly1803Ser) and four BRCA2 variants (p.Trp31Gly, p.Ser2616Phe, p.Tyr2660Cys, and p.Leu2792Arg) were pathogenic. This study demonstrates that integrative evaluation using conventional HR assays and optimized in silico prediction comprehensively classified the significance of BRCA VUSs for future clinical applications.

BRCA1 transports the DNA damage signal for CDDP-induced centrosome amplification through the centrosomal Aurora A.

Breast cancer gene 1 (BRCA1) plays roles in DNA repair and centrosome regulation and is involved in DNA damage-induced centrosome amplification (DDICA). Here, the centrosomal localization of BRCA1 and the kinases involved in centrosome duplication were analyzed in each cell cycle phase after treatment with DNA crosslinker cisplatin (CDDP). CDDP treatment increased the centrosomal localization of BRCA1 in early S-G2 phase. BRCA1 contributed to the increased centrosomal localization of Aurora A in S phase and that of phosphorylated Polo-like kinase 1 (PLK1) in late S phase after CDDP treatment, resulting in centriole disengagement and overduplication. The increased centrosomal localization of BRCA1 and Aurora A induced by CDDP treatment involved the nuclear export of BRCA1 and BRCA1 phosphorylation by ataxia telangiectasia mutated (ATM). Patient-derived variants and mutations at phosphorylated residues of BRCA1 suppressed the interaction between BRCA1 and Aurora A, as well as the CDDP-induced increase in the centrosomal localization of BRCA1 and Aurora A. These results suggest that CDDP induces the phosphorylation of BRCA1 by ATM in the nucleus and its transport to the cytoplasm, thereby promoting the centrosomal localization Aurora A, which phosphorylates PLK1. The function of BRCA1 in the translocation of the DNA damage signal from the nucleus to the centrosome to induce centrosome amplification after CDDP treatment might support its role as a tumor suppressor.

The BRCA2 missense mutation K2497R suppressed self-degradation and increased ATP production and cell proliferation.

Breast cancer susceptibility gene 2 (BRCA2) mediates genome maintenance during the S phase of the cell cycle, with important roles in replication stress, centrosome replication, and cytokinesis. In this study, we showed that a small heat shock protein, HSP27, interacted with and participated in the degradation of BRCA2 in estrogen-treated MCF-7 cells. BRCA2 degradation reportedly requires ubiquitination of the C-terminal region; thus, fragments of amino acid (aa) residues 2241-2940 were produced and assayed for their degradation following cycloheximide (CHX) treatment. The results showed that aa 2491-2580 affected the degradation of BRCA2, especially lysine (Lys) 2497. Furthermore, the K2497 A/R mutation increased ATP production and the proliferation of DLD-1 (BRCA2 knockout) cells compared to the cells expressing wild-type BRCA2-FLAG. Notably, a single residue, Lys2497, affected BRCA2 degradation, and K2497R is reportedly a missense mutation in hereditary breast cancer.

BRCA2 represses the transcriptional activity of pS2 by E2-ERα.

Germline mutations to the breast cancer 2 (BRCA2) gene have been associated with hereditary breast cancer. In addition to estrogen uptake, BRCA2 expression increases in the S phase of the cell cycle and largely contributes to DNA damage repair associated with DNA replication. However, the role of BRCA2 in estrogen induction remains unclear. An expression plasmid was created to induce BRCA2 activation upon the addition of estradiol by introducing mutations to the binding sequences for the transcription factors USF1, E2F1, and NF-κB within the promoter region of BRCA2. Then, the estrogen receptor (ER) sites of the proteins that interact with BRCA2 upon the addition of estradiol were identified. Both proteins were bound by the helical domain of BRCA2 and activation function-2 of the ER, suggesting that this binding may regulate the transcriptional activity of pS2, a target gene of the estradiol-ER, by suppressing the binding of SRC-1, a coactivator required for activation of the transcription factor.

Ultradeep targeted sequencing of circulating tumor DNA in plasma of early and advanced breast cancer.

We present a study to evaluate the feasibility and clinical utility of amplicon-based Oncomine Pan-Cancer cell-free assay to detect circulating tumor DNA (ctDNA) in patients with early or advanced breast cancer. In this study, 109 early and metastatic breast cancer patients were recruited before the initiation of treatment. ctDNA mutation profiles were assessed through unique molecular tagging (UMT) and ultradeep next generation sequencing (NGS). For patients with mutations, DNA from corresponding white blood cells (WBC) was sequenced to exclude variants of clonal-hematopoietic (CH) origin. UMT targeted sequencing from plasma of 109 patients achieved a median total coverage of 55 498X and a median molecular coverage of 4187X. Among 53 ctDNA positive samples, 38% were mutation positive by WBC sequencing, indicating potentially false-positive results contributed by CH origin. Prevalence of CH-related mutations was associated with age (P = 7.51 × 10-4 ). After exclusion of CH mutations, ctDNA detection rates were 37% for local or locally advanced breast cancer (stage I-III) and 81% for metastatic or recurrent breast cancer. The ctDNA detection rate correlated with disease stage (P = 2.60 × 10-4 ), nodal spread (P = 6.49 × 10-3 ) and the status of distant metastases (P = 5.00 × 10-4 ). ctDNA variants were detected mostly in TP53, PIK3CA and AKT1 genes, with variants showing therapeutic relevance. This pilot study endorses the use of targeted NGS for non-invasive molecular profiling of breast cancer. Paired sequencing of plasma ctDNA and WBC should be implemented to improve accurate interpretation of liquid biopsy.

Pathogenicity assessment of variants for breast cancer susceptibility genes based on BRCAness of tumor sample.

Genes involved in the homologous recombination repair pathway-as exemplified by BRCA1, BRCA2, PALB2, ATM, and CHEK2-are frequently associated with hereditary breast and ovarian cancer syndrome. Germline mutations in the loci of these genes with loss of heterozygosity or additional somatic truncation at the WT allele lead to the development of breast cancers with characteristic clinicopathological features and prominent genomic features of homologous recombination deficiency, otherwise referred to as "BRCAness." Although clinical genetic testing for these and other genes has increased the chances of identifying pathogenic variants, there has also been an increase in the prevalence of variants of uncertain significance, which poses a challenge to patient care because of the difficulties associated with making further clinical decisions. To overcome this challenge, we sought to develop a methodology to reclassify the pathogenicity of these unknown variants using statistical modeling of BRCAness. The model was developed with Lasso logistic regression by comparing 116 genomic attributes derived from 37 BRCA1/2 biallelic mutant and 32 homologous recombination-quiescent breast cancer exomes. The model showed 95.8% and 86.7% accuracies in the training cohort and The Cancer Genome Atlas validation cohort, respectively. Through application of the model for variant reclassification of homologous recombination-associated hereditary breast and ovarian cancer causal genes and further assessment with clinicopathological features, we finally identified one likely pathogenic and five likely benign variants. As such, the BRCAness model developed from the tumor exome was robust and provided a reasonable basis for variant reclassification.

CDK1 inhibitor controls G2/M phase transition and reverses DNA damage sensitivity.

CDK1 plays key roles in cell cycle progression through the G2/M phase transition and activation of homologous recombination (HR) DNA repair pathway. Accordingly, various CDK1 inhibitors have been developed for cancer therapy that induce prolonged G2 arrest and/or sensitize cells to DNA damaging agents in tumor cells, resulting in cell death. However, CDK1 inhibition can induce resistance to DNA damage in certain conditions. The mechanism of different DNA damage sensitivity is not completely understood. We performed immunofluorescence and flow cytometry analysis to investigate DNA damage responses in human tumor cells during low and high dose treatments with RO-3306, a selective CDK1 inhibitor. This comparative investigation demonstrated that RO-3306-induced G2 arrest prevented cells with DNA double-strand breaks from transitioning into the M-phase and that the cells maintained their DNA repair capacity in G2-phase, even under RO-3306 dose-dependent DNA repair inhibition. These findings reveal that CDK1 inhibitor-induced DNA repair inhibition and cell cycle control, which regulate each other during the G2/M phase transition determine the cellular sensitivity to DNA damage, providing insight useful for developing clinical strategies targeting CDK1 inhibition in tumor cells.

Non-V600E BRAF mutations and EGFR signaling pathway in colorectal cancer.

The Raf murine sarcoma viral oncogene homolog B (BRAFV600E ) mutation (MT) in metastatic colorectal cancer (CRC) is a well-known prognostic indicator and a negative predictive biomarker for antiepidermal growth factor receptor (EGFR) treatment. However, the clinical characteristics and significance of BRAFnon-V600E MTs remain unclear. Here, we evaluated the clinical characteristics of BRAFnon-V600E MTs vs. those of other MTs in the EGFR signaling pathway, including BRAFV600E . Consecutive CRC patients in our institute from June 2012 to November 2013 were enrolled in our study. Multiplex genotyping of the EGFR pathway was performed with archival samples using a Luminex Assay for BRAFV600E /BRAFnon-V600E , KRAS/NRAS exons 2-4, and phosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic subunit alpha (PIK3CA). We analyzed correlations among the MT profiles, clinical data and primary tumor locations in CRC. All statistical analyses were performed using R software. CRC samples (824) from 374 (45.4%) male and 450 (54.6%) female patients were analyzed, of which 154 (18.7%), 202 (24.5%), 270 (32.8%) or 198 (24.0%) had Stages I, II, III or IV or recurrent CRC, respectively. The frequencies of BRAFV600E /BRAFnon-V600E , KRAS (including exons 2-4), NRAS and PIK3CA MTs were 5.3/1.7, 41.4, 3.3 and 9.6%, respectively. The characteristics of patients with the BRAFV600E MT were an age of ≥65 years old, a right-sided primary tumor location, poorly differentiated histology and an advanced disease stage. In contrast, the characteristics of patients with BRAFnon-V600E MTs were a left-sided primary tumor location and well-differentiated histology. BRAFnon-V600E MTs were relatively rare and showed different characteristics compared to the BRAFV600E MT. These results may contribute to future precision medicine.

Crosstalk of DNA double-strand break repair pathways in poly(ADP-ribose) polymerase inhibitor treatment of breast cancer susceptibility gene 1/2-mutated cancer.

Germline mutations in breast cancer susceptibility gene 1 or 2 (BRCA1 or BRCA2) significantly increase cancer risk in hereditary breast and ovarian cancer syndrome (HBOC). Both genes function in the homologous recombination (HR) pathway of the DNA double-strand break (DSB) repair process. Therefore, the DNA-repair defect characteristic of cancer cells brings about a therapeutic advantage for poly(ADP-ribose) polymerase (PARP) inhibitor-induced synthetic lethality. PARP inhibitor-based therapeutics initially cause cancer lethality but acquired resistance mechanisms have been found and need to be elucidated. In particular, it is essential to understand in detail the mechanism of DNA damage and repair to PARP inhibitor treatment. Further investigations have shown the roles of BRCA1/2 and its associations to other molecules in the DSB repair system. Notably, the repair pathway chosen in BRCA1-deficient cells could be entirely different from that in BRCA2-deficient cells after PARP inhibitor treatment. The present review describes synthetic lethality and acquired resistance mechanisms to PARP inhibitor through the DSB repair pathway and subsequent repair process. In addition, recent knowledge of resistance mechanisms is discussed. Our model should contribute to the development of novel therapeutic strategies.

Phase II trial of biweekly cetuximab and irinotecan as third-line therapy for pretreated KRAS exon 2 wild-type colorectal cancer.

Efficacy and safety of biweekly cetuximab plus irinotecan were evaluated to provide guidance for its use in Japan as third-line treatment for pretreated metastatic colorectal cancer (mCRC) patients harboring wild-type KRAS exon 2. Objective response rate (ORR) was used as primary endpoint based on an expected proportion of 0.23 with confidence width of 0.298 (95% CI, 0.105-0.403), which showed 35 to be the minimal participant number. Forty patients, refractory to first- and second-line chemotherapy containing irinotecan, oxaliplatin, and fluoropyrimidine, were enrolled. ORR and disease control rate were 25.0% (95% CI: 11.5-38.4) and 72.5% (95% CI: 56.8-86.4), respectively. Median progression-free survival (PFS), overall survival (OS), and number of courses were 5.70 months (95% CI: 2.7-7.9), 15.1 months (95% CI: 11.8-19.0), and 10.5 (range: 3.0-31.0), respectively. Grade 3 adverse events were skin toxicity (12.5%), diarrhea (10.0%), neutropenia (5.0%), febrile neutropenia (5.0%), nausea (5.0%), anorexia (5.0%), and fatigue (2.5%). Cmax mean was 723.2 µg/mL after first dose. High area under the curve (AUC)last variance was associated with t1/2 range of 131.2-1209.6 hours (median, 174.4 hours). Early tumor shrinkage (ETS) and median depth of response were 25.0% and 13.0%, respectively. Mutation frequencies in KRAS exon 3 or 4, NRAS, BRAF, and PIK3CA were 5.5%, 2.7%, 8.3%, and 5.5%, respectively. Multivariate Cox regression analysis assessed whether any gene mutations and ETS are predictors for PFS, and whether performance status, synchronous metastasis, and ETS are predictors for OS. Importantly, the data provide guidance for a biweekly cetuximab plus irinotecan regimen in mCRC patients.

Three DNA polymerases, recruited by different mechanisms, carry out NER repair synthesis in human cells.

Nucleotide excision repair (NER) is the most versatile DNA repair system that deals with the major UV photoproducts in DNA, as well as many other DNA adducts. The early steps of NER are well understood, whereas the later steps of repair synthesis and ligation are not. In particular, which polymerases are definitely involved in repair synthesis and how they are recruited to the damaged sites has not yet been established. We report that, in human fibroblasts, approximately half of the repair synthesis requires both pol kappa and pol delta, and both polymerases can be recovered in the same repair complexes. Pol kappa is recruited to repair sites by ubiquitinated PCNA and XRCC1 and pol delta by the classical replication factor complex RFC1-RFC, together with a polymerase accessory factor, p66, and unmodified PCNA. The remaining repair synthesis is dependent on pol epsilon, recruitment of which is dependent on the alternative clamp loader CTF18-RFC.

BRCA1 gene: function and deficiency.

The BRCA1 protein, a hereditary breast and ovarian cancer-causing gene product, is known as a multifunctional protein that performs various functions in cells. It is well known, along with BRCA 2, to cause hereditary breast and ovarian cancer, but here we will specifically focus on BRCA1. We introduce the mechanism and the latest report on homologous recombination repair, replication, involvement in checkpoint regulation, transcription, chromatin remodeling, and cytoplasmic function (centrosome regulation, apoptosis, selective autophagy), and consider the possibility of carcinogenesis from inhibition of the intracellular functions in each. We also consider the possibility of drug development based on each function. Finally, we will explain, from data obtained through basic research, that an appropriate regimen is important for raising the response rate for poly (ADP)-ribose polymerase inhibitors, in the case of low susceptibility, iatrogenic toxicity, tolerance, etc.

FKBP51 regulates cell motility and invasion via RhoA signaling.

FK506 binding protein 51 (FKBP51), a member of the immunophilin family, is involved in multiple signaling pathways, tumorigenesis, and chemoresistance. FKBP51 expression correlates with metastatic potential in melanoma and prostate cancer. However, the functions of FKBP51, particularly involving the regulation of cell motility and invasion, are not fully understood. We discovered two novel interacting partner proteins of FKBP51, i.e., deleted in liver cancer 1 (DLC1) and deleted in liver cancer 2 (DLC2), using immunoprecipitation and mass spectrometry. DLC1 and DLC2 are Rho GTPase-activating proteins that are frequently downregulated in various cancers. Next, we demonstrated that overexpression of FKBP51 enhances cell motility and invasion of U2OS cells via upregulation of RhoA activity and enhanced Rho-ROCK signaling. Moreover, FKBP51-depleted cells displayed a cortical distribution of actin filaments and decreased cell motility and invasion. Consistent with this phenotype, FKBP51 depletion caused a downregulation of RhoA activity. Considered together, our results demonstrate that FKBP51 positively controls cell motility by promoting RhoA and ROCK activation; thus, we have revealed a novel role for FKBP51 in cytoskeletal rearrangement and cell migration and invasion.

Mutation status of RAD51C, PALB2 and BRIP1 in 100 Japanese familial breast cancer cases without BRCA1 and BRCA2 mutations.

In addition to BRCA1 and BRCA2, RAD51C, PALB2 and BRIP1 are known as breast cancer susceptibility genes. However, the mutation status of these genes in Japanese familial breast cancer cases has not yet been evaluated. To this end, we analyzed the exon sequence and genomic rearrangement of RAD51C, PALB2 and BRIP1 in 100 Japanese patients diagnosed with familial breast and ovarian cancer and without BRCA1 and BRCA2 mutations. We detected a large deletion from exons 6 to 9 in RAD51C, 4 novel BRIP1 missense variants containing 3 novel non-synonymous variants, c.89A>C, c.736A>G and c.2131A>G, and a splice donor site variant c.918+2T>C. No deleterious variant of PALB2 was detected. The results of pedigree analysis showed that the proband with a large deletion on RAD51C had a family history of both breast and ovarian cancer, and the families of probands with novel BRIP1 missense variants included a male patient with breast cancer or many patients with breast cancer within the second-degree relatives. We showed that the mutation frequency of RAD51C in Japanese familial breast cancer cases was similar to that in Western countries and that the prevalence of deleterious mutation of PALB2 was possibly lower. Furthermore, our results suggested that BRIP1 mutation frequency in Japan might differ from that in Western countries.

LAMC2 is a predictive marker for the malignant progression of leukoplakia.

BACKGROUND: LAMC2 plays an important role in cancer invasion. The aim of this study was to (i) compare the immunoexpression of LAMC2 in different stages of oral squamous cell carcinoma (OSCC), early and advanced, and (ii) to evaluate LAMC2 as a marker of malignant transformation in leukoplakia. MATERIALS AND METHODS: The expression of LAMC2 was examined by immunohistochemistry in 50 surgical specimens of advanced OSCC assembled as tissue microarrays, and by cDNA microarray in 43 surgical specimens of advanced OSCC. LAMC2 expression was further examined in 39 surgical specimens of early OSCC and in 93 incisional biopsy specimens of leukoplakia of the tongue, which exhibited epithelial dysplasia. The relationship of LAMC2 expression score with clinico-pathological characteristics was analyzed. RESULTS: LAMC2 was remarkably upregulated in OSCC at the cancer-stroma interface. The grade of LAMC2 expression was significantly associated with the pattern and depth of invasion of OSCC. Foci of LAMC2-positive cells were observed in some cases of leukoplakia. The number and size of LAMC2-positive foci were significantly associated with the grade of dysplasia. The presence of LAMC2-positive foci was a significant predictive factor for the malignant progression of leukoplakia. LAMC2-positive leukoplakia had an approximately 11-fold increased risk of malignancy compared with LAMC2-negative leukoplakia. CONCLUSIONS: The results of this study highlight the value of LAMC2 as a marker of cancer invasion. LAMC2-positive foci in leukoplakia suggest an imminent risk of cancer. LAMC2 immunostaining is expected to contribute to a more precise assessment of the malignancy of leukoplakia.

[Molecular Diagnosis and Treatment of HBOC Syndrome].

Hereditary breast and ovarian cancer(HBOC)is an inherited cancer caused by mutations of the BRCA1 or BRCA2 genes. BRCA genetic testing is used for HBOC diagnosis and continues to progress such as the annotation of VUS. In HBOC clinical practice, surveillance methods have been established through collaboration between genetic medicine and cancer medicine, and treatment, including options based on genetic diagnosis, has advanced significantly. Furthermore, the analysis of BRCA1 and BRCA2 function has progressed, and a novel therapeutic method based on synthetic lethality, such as a PARP inhibitor use, has been developed. Furthermore, BRCA genetic testing is going to be used as a PGx test for the selection of sensitive cases. Meanwhile, familial breast cancer and ovarian cancer, in which cases of breast and ovarian cancer accumulate in the family, vary from some patients carrying a single mutated gene, such as BRCA1, to families that have multifactorial predisposing causes. Responsible genes of each group have been identified as high, moderate, and low susceptibility genes, and there are a number of cases where the responsible genes are unknown. Such genes need to be identified and a new diagnostic system needs to be established.

RAS mutation is a prognostic biomarker in colorectal cancer patients with metastasectomy.

Studies have demonstrated a relationship between clinical outcomes after curative resection for colorectal cancer (CRC) and gene mutations of the EGFR pathway; however, no studies have examined metastatic CRC (mCRC) patients with metastasectomy. The aim of this study was to evaluate the relationship between gene mutations of EGFR pathway and clinical outcomes after metastasectomy in mCRC patients. A total of 1,053 patients histopathologically confirmed CRC received a genotyping test for the EGFR pathway from February 2012 to October 2013. Detailed information was obtained through review of medical records. Gene mutations of EGFR pathway were analyzed by Luminex assay. Overall survival (OS) and recurrence free survival were estimated by the Kaplan-Meier method and the log-rank test was used to compare the survival outcomes by gene mutation status. A total of 132 patients received metastasectomy. The frequencies of KRAS exon 2, KRAS exon 3.4, NRAS, BRAF, and PIK3CA mutations were 38.6% (51/132), 3.6% (5/132), 5.1% (7/132), 5.1% (7/132), and 8.7% (12/132), respectively. With a median follow-up of 84.1 months (57.2-NA) for a survivor, the 4-year OS rate was 65.6% for mCRC with RAS mutation, and 81.3% for mCRC with wild-type RAS (p < 0.05). We observed a statistically significant correlation for only the RAS mutation and OS. In multivariate analysis, RAS mutation and liver metastasis were independent factors for shorter OS. There were no significant differences between gene mutations of EGFR pathway and recurrence free survival. RAS mutation in mCRC metastasectomy patients was associated with shorter overall survival.

THBS1 is induced by TGFB1 in the cancer stroma and promotes invasion of oral squamous cell carcinoma.

BACKGROUND: THBS1 (thrombospondin-1) is the extracellular matrix (ECM) protein that affects diverse cellular activities. It constitutes the tumor stroma, but the role of THBS1 in oral squamous cell carcinoma (OSCC) development is unclear. The aim of this study was to clarify the relevance of THBS1 in the pathogenesis of OSCC. MATERIALS AND METHODS: The expression of THBS1 was examined in 44 OSCC by immunohistochemical analysis and in 43 OSCC by cDNA microarray analysis. Cell culture experiments were conducted using human OSCC cell lines HSC3 and HO1N1 and mouse fibroblast ST2 cells to examine the effect of TGFB1 on THBS1 expression, and the effect of THBS1 on cellular behaviors. RESULTS: THBS1 was specifically induced in the tumor microenvironment of OSCC. THBS1 appeared to be produced mainly by the stromal cells, but also by OSCC cells. TGFB1 stimulated THBS1 expression in ST2, primary fibroblasts, and the OSCC cells. THBS1 promoted migration and invasion of HSC3 and HO1N1 in transwell migration assays. THBS1 stimulated the expression of MMP3 (matrix metalloprotease 3), MMP9, MMP11, and MMP13 in ST2 cells and MMP3, MMP11, and MMP13 in HO1N1 cells. The RGD peptide suppressed the THBS1-stimulated migration and upregulation of MMP11 and MMP13. CONCLUSIONS: THBS1 is a tumor-specific ECM protein that is induced by TGFB1 and promotes migration of cancer cells and stimulates the expression of MMPs partly through the integrin signaling, thereby favoring OSCC invasion.