Index of Cancer-Associated Fibroblasts Is Superior to the Epithelial-Mesenchymal Transition Score in Prognosis Prediction.

In many solid tumors, tissue of the mesenchymal subtype is frequently associated with epithelial-mesenchymal transition (EMT), strong stromal infiltration, and poor prognosis. Emerging evidence from tumor ecosystem studies has revealed that the two main components of tumor stroma, namely, infiltrated immune cells and cancer-associated fibroblasts (CAFs), also express certain typical EMT genes and are not distinguishable from intrinsic tumor EMT, where bulk tissue is concerned. Transcriptomic analysis of xenograft tissues provides a unique advantage in dissecting genes of tumor (human) or stroma (murine) origins. By transcriptomic analysis of xenograft tissues, we found that oral squamous cell carcinoma (OSCC) tumor cells with a high EMT score, the computed mesenchymal likelihood based on the expression signature of canonical EMT markers, are associated with elevated stromal contents featured with fibronectin 1 (Fn1) and transforming growth factor-ß (Tgfß) axis gene expression. In conjugation with meta-analysis of these genes in clinical OSCC datasets, we further extracted a four-gene index, comprising FN1, TGFB2, TGFBR2, and TGFBI, as an indicator of CAF abundance. The CAF index is more powerful than the EMT score in predicting survival outcomes, not only for oral cancer but also for the cancer genome atlas (TCGA) pan-cancer cohort comprising 9356 patients from 32 cancer subtypes. Collectively, our results suggest that a further distinction and integration of the EMT score with the CAF index will enhance prognosis prediction, thus paving the way for curative medicine in clinical oncology.

Frequent overexpression of klotho in fusion-negative phosphaturic mesenchymal tumors with tumorigenic implications.

Phosphaturic mesenchymal tumors (PMT) are tumors that cause hypophosphatemia/osteomalacia chiefly by secreting FGF23. We have identified FN1-FGFR1/FGF1 fusion genes in nearly half of PMT, suggesting a central role of FGFR1 pathways in the pathogenesis of PMT. Tumorigenic drivers are unknown for tumors where previous study detected neither fusion, including many in bone, where FISH failed because of tissue decalcification. To identify alternative fusions in PMT without known fusions, as well as to validate the positive FISH results and characterize the fusion junctions, 34 PMT were studied, including 12 with known FN1-FGFR1 fusion by FISH (Group A), 2 with FN1-FGF1 (B), 12 with neither fusion (C), and 8 with previous acid-based decalcification and hence unknown fusion status (D). In total, 23 archival samples were subjected to anchored multiplex PCR-based RNA-sequencing (AMP-seq) with primers targeting FN1, genes encoding the FGF/FGFR families, and KL (α-Klotho); five Group C cases were also studied with whole-transcriptomic and exome-captured RNA sequencing, respectively. The AMP-seq results were consistent with previous FISH and/or transcriptomic sequencing data, except in one old Group A sample. One case had a novel FGFR1 exon 9 breakpoint, confirmed by genomic DNA sequencing. One Group D bone tumor was found to harbor FN1-FGF1. All 3 RNA-sequencing platforms failed to identify convincing fusion genes in Group C (N = 10), which instead expressed significantly higher levels of either KL or KLB. This result was further confirmed with KL and KLB RNA CISH semi-quantification (RNAscope). Our results demonstrated the utility of AMP-seq, which was compromised by decalcification and prolonged archiving. Of potential importance, fusion-negative PMT frequently overexpressed α-Klotho (or instead ß-Klotho less commonly), whose role as an obligatory co-receptor for FGF23-FGFR1 binding suggests its aberrant expression in osteocytes/osteoblasts might result in an FGF23-FGFR1 autocrine loop that in turn drives the overexpression of FGF23 and tumorigenesis through activated FGFR pathways.

A gene profiling deconvolution approach to estimating immune cell composition from complex tissues.

BACKGROUND: A new emerged cancer treatment utilizes intrinsic immune surveillance mechanism that is silenced by those malicious cells. Hence, studies of tumor infiltrating lymphocyte populations (TILs) are key to the success of advanced treatments. In addition to laboratory methods such as immunohistochemistry and flow cytometry, in silico gene expression deconvolution methods are available for analyses of relative proportions of immune cell types. RESULTS: Herein, we used microarray data from the public domain to profile gene expression pattern of twenty-two immune cell types. Initially, outliers were detected based on the consistency of gene profiling clustering results and the original cell phenotype notation. Subsequently, we filtered out genes that are expressed in non-hematopoietic normal tissues and cancer cells. For every pair of immune cell types, we ran t-tests for each gene, and defined differentially expressed genes (DEGs) from this comparison. Equal numbers of DEGs were then collected as candidate lists and numbers of conditions and minimal values for building signature matrixes were calculated. Finally, we used v -Support Vector Regression to construct a deconvolution model. The performance of our system was finally evaluated using blood biopsies from 20 adults, in which 9 immune cell types were identified using flow cytometry. The present computations performed better than current state-of-the-art deconvolution methods. CONCLUSIONS: Finally, we implemented the proposed method into R and tested extensibility and usability on Windows, MacOS, and Linux operating systems. The method, MySort, is wrapped as the Galaxy platform pluggable tool and usage details are available at https://testtoolshed.g2.bx.psu.edu/view/moneycat/mysort/e3afe097e80a .

Characterization of FN1-FGFR1 and novel FN1-FGF1 fusion genes in a large series of phosphaturic mesenchymal tumors.

Phosphaturic mesenchymal tumors typically cause paraneoplastic osteomalacia, chiefly as a result of FGF23 secretion. In a prior study, we identified FN1-FGFR1 fusion in 9 of 15 phosphaturic mesenchymal tumors. In this study, a total of 66 phosphaturic mesenchymal tumors and 7 tumors resembling phosphaturic mesenchymal tumor but without known phosphaturia were studied. A novel FN1-FGF1 fusion gene was identified in two cases without FN1-FGFR1 fusion by RNA sequencing and cross-validated with direct sequencing and western blot. Fluorescence in situ hybridization analyses revealed FN1-FGFR1 fusion in 16 of 39 (41%) phosphaturic mesenchymal tumors and identified an additional case with FN1-FGF1 fusion. The two fusion genes were mutually exclusive. Combined with previous data, the overall prevalence of FN1-FGFR1 and FN1-FGF1 fusions was 42% (21/50) and 6% (3/50), respectively. FGFR1 immunohistochemistry was positive in 82% (45/55) of phosphaturic mesenchymal tumors regardless of fusion status. By contrast, 121 cases of potential morphologic mimics (belonging to 13 tumor types) rarely expressed FGFR1, the main exceptions being solitary fibrous tumors (positive in 40%), chondroblastomas (40%), and giant cell tumors of bone (38%), suggesting a possible role for FGFR1 immunohistochemistry in the diagnosis of phosphaturic mesenchymal tumor. With the exception of one case reported in our prior study, none of the remaining tumors resembling phosphaturic mesenchymal tumor had either fusion type or expressed significant FGFR1. Our findings provide insight into possible mechanisms underlying the pathogenesis of phosphaturic mesenchymal tumor and imply a central role of the FGF1-FGFR1 signaling pathway. The novel FN1-FGF1 protein is expected to be secreted and serves as a ligand that binds and activates FGFR1 to achieve an autocrine loop. Further study is required to determine the functions of these fusion proteins.

Over-expression of AURKA, SKA3 and DSN1 contributes to colorectal adenoma to carcinoma progression.

Development of colorectal cancer (CRC) involves sequential transformation of normal mucosal tissues into benign adenomas and then adenomas into malignant tumors. The identification of genes crucial for malignant transformation in colorectal adenomas (CRAs) has been based primarily on cross-sectional observations. In this study, we identified relevant genes using autologous samples. By performing genome-wide SNP genotyping and RNA sequencing analysis of adenocarcinomas, adenomatous polyps, and non-neoplastic colon tissues (referred as tri-part samples) from individual patients, we identified 68 genes with differential copy number alterations and progressively dysregulated expression. Aurora A, SKA3, and DSN1 protein levels were sequentially up-regulated in the samples, and this overexpression was associated with chromosome instability (CIN). Knockdown of SKA3 in CRC cells dramatically reduced cell growth rates and increased apoptosis. Depletion of SKA3 or DSN1 induced G2/M arrest and decreased migration, invasion, and anchorage-independent growth. AURKA and DSN1 are thus critical for chromosome 20q amplification-associated malignant transformation in CRA. Moreover, SKA3 at chromosome 13q was identified as a novel gene involved in promoting malignant transformation. Evaluating the expression of these genes may help identify patients with progressive adenomas, helping to improve treatment.

Identification of a novel FN1-FGFR1 genetic fusion as a frequent event in phosphaturic mesenchymal tumour.

Phosphaturic mesenchymal tumours (PMTs) are uncommon soft tissue and bone tumours that typically cause hypophosphataemia and tumour-induced osteomalacia (TIO) through secretion of phosphatonins including fibroblast growth factor 23 (FGF23). PMT has recently been accepted by the World Health Organization as a formal tumour entity. The genetic basis and oncogenic pathways underlying its tumourigenesis remain obscure. In this study, we identified a novel FN1-FGFR1 fusion gene in three out of four PMTs by next-generation RNA sequencing. The fusion transcripts and proteins were subsequently confirmed with RT-PCR and western blotting. Fluorescence in situ hybridization analysis showed six cases with FN1-FGFR1 fusion out of an additional 11 PMTs. Overall, nine out of 15 PMTs (60%) harboured this fusion. The FN1 gene possibly provides its constitutively active promoter and the encoded protein's oligomerization domains to overexpress and facilitate the activation of the FGFR1 kinase domain. Interestingly, unlike the prototypical leukaemia-inducing FGFR1 fusion genes, which are ligand-independent, the FN1-FGFR1 chimeric protein was predicted to preserve its ligand-binding domains, suggesting an advantage of the presence of its ligands (such as FGF23 secreted at high levels by the tumour) in the activation of the chimeric receptor tyrosine kinase, thus effecting an autocrine or a paracrine mechanism of tumourigenesis.

A comprehensive functional map of the hepatitis C virus genome provides a resource for probing viral proteins.

UNLABELLED: Pairing high-throughput sequencing technologies with high-throughput mutagenesis enables genome-wide investigations of pathogenic organisms. Knowledge of the specific functions of protein domains encoded by the genome of the hepatitis C virus (HCV), a major human pathogen that contributes to liver disease worldwide, remains limited to insight from small-scale studies. To enhance the capabilities of HCV researchers, we have obtained a high-resolution functional map of the entire viral genome by combining transposon-based insertional mutagenesis with next-generation sequencing. We generated a library of 8,398 mutagenized HCV clones, each containing one 15-nucleotide sequence inserted at a unique genomic position. We passaged this library in hepatic cells, recovered virus pools, and simultaneously assayed the abundance of mutant viruses in each pool by next-generation sequencing. To illustrate the validity of the functional profile, we compared the genetic footprints of viral proteins with previously solved protein structures. Moreover, we show the utility of these genetic footprints in the identification of candidate regions for epitope tag insertion. In a second application, we screened the genetic footprints for phenotypes that reflected defects in later steps of the viral life cycle. We confirmed that viruses with insertions in a region of the nonstructural protein NS4B had a defect in infectivity while maintaining genome replication. Overall, our genome-wide HCV mutant library and the genetic footprints obtained by high-resolution profiling represent valuable new resources for the research community that can direct the attention of investigators toward unidentified roles of individual protein domains. IMPORTANCE: Our insertional mutagenesis library provides a resource that illustrates the effects of relatively small insertions on local protein structure and HCV viability. We have also generated complementary resources, including a website (http://hangfei.bol.ucla.edu) and a panel of epitope-tagged mutant viruses that should enhance the research capabilities of investigators studying HCV. Researchers can now detect epitope-tagged viral proteins by established antibodies, which will allow biochemical studies of HCV proteins for which antibodies are not readily available. Furthermore, researchers can now quickly look up genotype-phenotype relationships and base further mechanistic studies on the residue-by-residue information from the functional profile. More broadly, this approach offers a general strategy for the systematic functional characterization of viruses on the genome scale.