Lipoblastoma-Like Tumor and Fibrosarcoma-Like Lipomatous Neoplasm Represent the Same Entity: A Clinicopathologic and Molecular Genetic Study of 23 Cases Occurring in Both Men and Women at Diverse Locations.

Lipoblastoma-like tumor (LLT) is a benign soft tissue tumor demonstrating mixed morphologic features of lipoblastoma, myxoid liposarcoma, and spindle cell lipoma but lacking genetic alterations associated with those tumors. LLT was originally thought to be specific to the vulva but has since been reported in the paratesticular region. The morphologic features of LLT overlap with those of "fibrosarcoma-like lipomatous neoplasm" (FLLN), a rare, indolent adipocytic neoplasm considered by some to form part of the spectrum of atypical spindle cell and pleomorphic lipomatous tumor. We compared the morphologic, immunohistochemical, and genetic features of 23 tumors previously classified as LLT (n = 17) and FLLN (n = 6). The 23 tumors occurred in 13 women and 10 men (mean age, 42 years; range, 17 to 80 years). Eighteen (78%) cases arose in the inguinogenital region, whereas 5 tumors (22%) involved noninguinogenital soft tissue, including the flank (n = 1), shoulder (n = 1), foot (n = 1), forearm (n = 1), and chest wall (n = 1). Microscopically, the tumors were lobulated and septated, with variably collagenized fibromyxoid stroma, prominent thin-walled vessels, scattered univacuolated or bivacuolated lipoblasts, and a minor component of mature adipose tissue. Using immunohistochemistry, 5 tumors (42%) showed complete RB1 loss, with partial loss in 7 cases (58%). RNA sequencing, chromosomal microarray, and DNA next-generation sequencing study results were negative for significant alterations. There were no clinical, morphologic, immunohistochemical, or molecular genetic differences between cases previously classified as LLT or FLLN. Clinical follow-up (11 patients [48%]; range, 2-276 months; mean, 48.2 months) showed all patients were alive without disease, and only one patient had experienced a single local recurrence. We conclude that LLT and FLLN represent the same entity, for which "LLT" seems most appropriate. LLT may occur in either sex and any superficial soft tissue location. Careful morphologic study and appropriate ancillary testing should allow for the distinction of LLT from its potential mimics.

Multi-Omic Architecture of Obstructive Hypertrophic Cardiomyopathy.

BACKGROUND: Hypertrophic cardiomyopathy (HCM) is characterized by asymmetric left ventricular hypertrophy. Currently, hypertrophy pathways responsible for HCM have not been fully elucidated. Their identification could serve as a nidus for the generation of novel therapeutics aimed at halting disease development or progression. Herein, we performed a comprehensive multi-omic characterization of hypertrophy pathways in HCM. METHODS: Flash-frozen cardiac tissues were collected from genotyped HCM patients (n=97) undergoing surgical myectomy and tissue from 23 controls. RNA sequencing and mass spectrometry-enabled deep proteome and phosphoproteomic assessment were performed. Rigorous differential expression, gene set enrichment, and pathway analyses were performed to characterize HCM-mediated alterations with emphasis on hypertrophy pathways. RESULTS: We identified transcriptional dysregulation with 1246 (8%) differentially expressed genes and elucidated downregulation of 10 hypertrophy pathways. Deep proteomic analysis identified 411 proteins (9%) that differed between HCM and controls with strong dysregulation of metabolic pathways. Seven hypertrophy pathways were upregulated with antagonistic upregulation of 5 of 10 hypertrophy pathways shown to be downregulated in the transcriptome. Most upregulated hypertrophy pathways encompassed the rat sarcoma-mitogen-activated protein kinase signaling cascade. Phosphoproteomic analysis demonstrated hyperphosphorylation of the rat sarcoma-mitogen-activated protein kinase system suggesting activation of this signaling cascade. There was a common transcriptomic and proteomic profile regardless of genotype. CONCLUSIONS: At time of surgical myectomy, the ventricular proteome, independent of genotype, reveals widespread upregulation and activation of hypertrophy pathways, mainly involving the rat sarcoma-mitogen-activated protein kinase signaling cascade. In addition, there is a counterregulatory transcriptional downregulation of the same pathways. Rat sarcoma-mitogen-activated protein kinase activation may serve a crucial role in hypertrophy observed in HCM.

Comparison of biochemical, microbial and mucosal mRNA expression in bile acid diarrhoea and irritable bowel syndrome with diarrhoea.

OBJECTIVE: There are altered mucosal functions in irritable bowel syndrome with diarrhoea (IBS-D); ~30% of patients with IBS-D have abnormal bile acid (BA) metabolism (ABAM) and diarrhoea (summarised as BAD). AIM: To compare biochemical parameters, gastrointestinal and colonic transit, rectal sensation and pathobiological mechanisms in IBS-D without ABAM and in BAD (serum 7C4>52 ng/mL). DESIGN: In patients with Rome III criteria of IBS-D, we compared biochemical features, colonic transit, rectal sensation, deep genotype of five BA-related genes, ileal and colonic mucosal mRNA (differential expression (DE) analysis) and stool dysbiosis (including functional analysis of microbiome). Results in BAD were compared with IBS-D without ABAM. RESULTS: Compared with 161 patients with IBS-D without ABAM, 44 patients with BAD had significantly faster colonic transit, lower microbial alpha diversity, different compositional profile (beta diversity) and higher Firmicutes to Bacteroidetes ratio with evidence of decreased expression of bile acid thiol ligase (involved in transformation of primary to secondary BAs) and decreased sulfatases. In BAD (compared with IBS-D without ABAM), terminal ileal biopsies showed downregulation of SLC44A5 (a BA transporter), and ascending colon biopsies showed upregulation in barrier-weakening genes (CLDN2), serine protease inhibitors, immune activation, cellular differentiation and a cellular transporter (FABP6; BA binding). No DE of genes was documented in descending colon biopsies. The two groups had similar rectal sensation. CONCLUSION: Though sharing clinical symptoms with IBS-D, BAD is associated with biological differences and mechanisms that have potential to enhance diagnosis and treatment targeting barrier dysfunction, inflammatory and microbial changes.

Polymorphous Low-Grade Neuroepithelial Tumor of the Young (PLNTY): Molecular Profiling Confirms Frequent MAPK Pathway Activation.

Polymorphous low-grade neuroepithelial tumor of the young (PLNTY) is a recently described epileptogenic tumor characterized by oligodendroglioma-like components, aberrant CD34 expression, and frequent mitogen-activated protein kinase (MAPK) pathway activation. We molecularly profiled 13 cases with diagnostic histopathological features of PLNTY (10 female; median age, 16 years; range, 5-52). Patients frequently presented with seizures (9 of 12 with available history) and temporal lobe tumors (9 of 13). MAPK pathway activating alterations were identified in all 13 cases. Fusions were present in the 7 youngest patients: FGFR2-CTNNA3 (n = 2), FGFR2-KIAA1598 (FGFR2-SHTN1) (n = 1), FGFR2-INA (n = 1), FGFR2-MPRIP (n = 1), QKI-NTRK2 (n = 1), and KIAA1549-BRAF (n = 1). BRAF V600E mutation was present in 6 patients (17 years or older). Two fusion-positive cases additionally harbored TP53/RB1 abnormalities suggesting biallelic inactivation. Copy number changes predominantly involving whole chromosomes were observed in all 10 evaluated cases, with losses of chromosome 10q occurring with FGFR2-KIAA1598 (SHTN1)/CTNNA3 fusions. The KIAA1549-BRAF and QKI-NTRK2 fusions were associated respectively with a 7q34 deletion and 9q21 duplication. This study shows that despite its name, PLNTY also occurs in older adults, who frequently show BRAF V600E mutation. It also expands the spectrum of the MAPK pathway activating alterations associated with PLNTY and demonstrates recurrent chromosomal copy number changes consistent with chromosomal instability.

EWSR1-WT1 gene fusions in neoplasms other than desmoplastic small round cell tumor: a report of three unusual tumors involving the female genital tract and review of the literature.

Desmoplastic small round cell tumor (DSRCT) is a high-grade round cell sarcoma that typically arises in the abdominopelvic cavity of young males, co-expresses keratins and desmin, and carries a pathognomonic EWSR1-WT1 gene fusion. The EWSR1-WT1 gene fusion is generally considered specific for DSRCT, although there are two reports of this fusion in tumors otherwise lacking features of DSRCT. We report three female genital tract tumors with EWSR1-WT1 fusions but showing morphologic and immunohistochemical features incompatible with DSRCT. The tumors occurred in the uterine cervix, uterine corpus/ovaries, and vagina, respectively, of 46, 30, and 20-year-old women. Two tumors consisted of a sheet-like to fascicular proliferation of relatively uniform spindled to occasionally more epithelioid cells arrayed about thick-walled, hyalinized, and capillary-sized vessels, with distinctive areas of pseudovascular change, and absence of desmoplastic stroma. The third tumor resembled a monomorphic spindle cell sarcoma with necrosis. All had diffuse desmin and variable but more limited keratin expression, two of three expressed smooth muscle actin, and all were negative for h-caldesmon, CD10, estrogen receptor, myogenin, N-terminus WT-1, and S100 protein. One patient received neoadjuvant chemotherapy and radiation therapy followed by resection and is disease-free 42 months after diagnosis. Another patient was managed by resection only and is disease-free 9 months after initial diagnosis. The remaining patient recently underwent resection of multifocal pelvic disease. Comprehensive differential gene expression analysis on two tumors compared to two classic DSRCTs with known EWSR1-WT1 fusions resulted in 1726 genes that were differentially expressed (log2 fold change >2 or < -2) and statistically significant (FDR < 5%). In combination with previous reports, our findings suggest pleiotropy of the EWSR1-WT1 fusion is possible and not limited to DSRCT. Subsets of non-DSRCT EWSR1-WT1 positive tumors may represent discrete entities, but further study is necessary.

RNA-Seq Reveals Differences in Expressed Tumor Mutation Burden in Colorectal and Endometrial Cancers with and without Defective DNA-Mismatch Repair.

Tumor mutation burden (TMB) is an emerging biomarker of immunotherapy response. RNA sequencing in FFPE tissue samples was used for determining TMB in microsatellite-stable (MSS) and microsatellite instability-high (MSI-H) tumors in patients with colorectal or endometrial cancer. Tissue from tumors and paired normal tissue from 46 MSI-H and 12 MSS cases were included. Of the MSI-H tumors, 29 had defective DNA mismatch-repair mutations, and 17 had MLH1 promoter hypermethylation. TMB was measured using the expressed somatic nucleotide variants (eTMB). A method of accurate measurement of eTMB was developed that removes FFPE-derived artifacts by leveraging mutation signatures. There was a significant difference in the median eTMB values observed between MSI-H and MSS cases: 27.3 versus 6.7 mutations/megabase (mut/Mb) (P = 3.5 × 10-9). Among tumors with defective DNA-mismatch repair, those with mismatch-repair mutations had a significantly higher median eTMB than those with hypermethylation: 28.1 versus 17.5 mut/Mb (P = 0.037). Multivariate analysis showed that MSI status, tumor type (endometrial or colorectal), and age were significantly associated with eTMB. Additionally, using whole-exome sequencing in a subset of these patients, it was determined that DNA TMB correlated well with eTMB (Spearman correlation coefficient, 0.83). These results demonstrate that RNA sequencing can be used for measuring eTMB in FFPE tumor specimens.

Frequency of MicroRNA Response Elements Identifies Pathologically Relevant Signaling Pathways in Triple-Negative Breast Cancer.

Complex interactions between mRNAs and microRNAs influence cellular functions. The mRNA-microRNA interactions also determine the post-transcriptional availability of mRNAs and unbound microRNAs. MicroRNAs binds to one or more microRNA response elements (MREs) located on the 3'UTR of mRNAs. In this study, we leveraged MREs and their frequencies in cancer and matched normal tissues to obtain insights into disease-specific interactions between mRNAs and microRNAs. We developed a bioinformatics method "ReMIx" that utilizes RNA sequencing (RNA-Seq) data to quantify MRE frequencies across the transcriptome. We applied ReMIx to triple-negative (TN) breast cancer tumor-normal adjacent pairs and identified MREs specific to TN tumors. ReMIx identified candidate mRNAs and microRNAs in the MAPK signaling cascade. Further analysis of MAPK gene regulatory networks revealed microRNA partners that influence and modulate MAPK signaling. In conclusion, we demonstrate a novel method of using MREs in the identification of functionally relevant mRNA-microRNA interactions in TN breast cancer.

Desmoplastic Infantile Ganglioglioma: A MAPK Pathway-Driven and Microglia/Macrophage-Rich Neuroepithelial Tumor.

MAPK pathway activation has been recurrently observed in desmoplastic infantile ganglioglioma/astrocytoma (DIG/DIA) with reported disproportionally low mutation allele frequencies relative to the apparent high tumor content, suggesting that MAPK pathway alterations may be subclonal. We sought to expand the number of molecularly profiled cases and investigate if tumor cell composition could account for the observed low mutation allele frequencies. Molecular (targeted neuro-oncology next-generation sequencing/RNA sequencing and OncoScan microarray) and immunohistochemical (CD68-PGM1/CD163/CD14/CD11c/lysozyme/CD3/CD20/CD34/PD-L1) studies were performed in 7 DIG. Activating MAPK pathway alterations were identified in 4 (57%) cases: 3 had a BRAF mutation (V600E/V600D/V600_W604delinsDQTDG, at 8%-27% variant allele frequency) and 1 showed a TPM3-NRTK1 fusion. Copy number changes were infrequent and nonrecurrent. All tumors had at least 30% of cells morphologically and immunophenotypically consistent with microglial/macrophage lineage. Two subtotally resected tumors regrew; 1 was re-excised and received adjuvant treatment (chemotherapy/targeted therapy), with clinical response to targeted therapy only. Even with residual tumor, all patients are alive (median follow-up, 83 months; 19-139). This study further supports DIG as another MAPK pathway-driven neuroepithelial tumor, thus expanding potential treatment options for tumors not amenable to surgical cure, and suggests that DIG is a microglia/macrophage-rich neuroepithelial tumor with frequent low driver mutation allele frequencies.

An expanded variant list and assembly annotation identifies multiple novel coding and noncoding genes for prostate cancer risk using a normal prostate tissue eQTL data set.

Prostate cancer (PrCa) is highly heritable; 284 variants have been identified to date that are associated with increased prostate cancer risk, yet few genes contributing to its development are known. Expression quantitative trait loci (eQTL) studies link variants with affected genes, helping to determine how these variants might regulate gene expression and may influence prostate cancer risk. In the current study, we performed eQTL analysis on 471 normal prostate epithelium samples and 249 PrCa-risk variants in 196 risk loci, utilizing RNA sequencing transcriptome data based on ENSEMBL gene definition and genome-wide variant data. We identified a total of 213 genes associated with known PrCa-risk variants, including 141 protein-coding genes, 16 lncRNAs, and 56 other non-coding RNA species with differential expression. Compared to our previous analysis, where RefSeq was used for gene annotation, we identified an additional 130 expressed genes associated with known PrCa-risk variants. We detected an eQTL signal for more than half (n = 102, 52%) of the 196 loci tested; 52 (51%) of which were a Group 1 signal, indicating high linkage disequilibrium (LD) between the peak eQTL variant and the PrCa-risk variant (r2>0.5) and may help explain how risk variants influence the development of prostate cancer.

RNA sequencing identifies a novel USP9X-USP6 promoter swap gene fusion in a primary aneurysmal bone cyst.

Primary aneurysmal bone cyst (ABC) is a benign multiloculated cystic lesion of bone that is defined cytogenetically by USP6 gene rearrangements. Rearrangements involving USP6 are promoter swaps, usually generated by fusion of the noncoding upstream exons of different partner genes with exon 1 or 2 of USP6, thus leading to transcriptional upregulation of full-length USP6 coding sequence. Testing for USP6 rearrangements is used diagnostically to distinguish it from secondary ABC and other giant cell-rich primary bone tumors. In this report, we present a case of a 16-year-old male with a primary ABC of the left distal femur. USP6 break apart fluorescence in situ hybridization was positive for a rearrangement and conventional chromosome analysis identified a reciprocal X;17 translocation. In order to identify the putative USP6 fusion partner, we performed RNA sequencing and uncovered a novel USP9X-USP6 promoter swap fusion. This result was confirmed by reverse transcription-polymerase chain reaction (RT-PCR) and by mate pair sequencing thus showing the utility of these alternative methodologies in identifying novel fusion candidates. Ubiquitin-specific protease 9X (USP9X), like USP6, encodes a highly conserved substrate-specific deubiquitylating enzyme. USP9X is highly expressed in a number of tissue types and acts as both an oncogene and tumor suppressor in several human cancers. We conclude that oncogenic activation of USP6 via USP9X promoter exchange represents a novel driver of primary ABC formation.

Characterization of three ciliopathy pedigrees expands the phenotype associated with biallelic C2CD3 variants.

Whole exome sequencing (WES) is utilized in diagnostic odyssey cases to identify the underlying genetic cause associated with complex phenotypes. Recent publications suggest that WES reveals the genetic cause in ~25% of these cases and is most successful when applied to children with neurological disease. The residual 75% of cases remain genetically elusive until more information becomes available in the literature or functional studies are pursued. WES performed on three families with presumed ciliopathy diagnoses, including orofaciodigital (OFD) syndrome, fetal encephalocele, or Joubert-related disorder, identified compound heterozygous variants in C2CD3. Biallelic variants in C2CD3 have previously been associated with ciliopathies, including OFD syndrome type 14 (OFD14; MIM: 615948). As three of the six identified variants were predicted to affect splicing, exon-skipping analysis using either RNA sequencing or PCR-based methods were completed to determine the pathogenicity of these variants, and showed that each of the splicing variants led to a frameshifted protein product. Using these studies in combination with the 2015 ACMG guidelines, each of the six identified variants were classified as either pathogenic or likely pathogenic, and are therefore likely responsible for our patients' phenotypes. Each of the families had a distinct clinical phenotype and severity of disease, extending from lethal to viable. These findings highlight that there is a broad phenotypic spectrum associated with C2CD3-mediated disease and not all patients present with the typical features of OFD14.

Transcriptomic signatures reveal immune dysregulation in human diabetic and idiopathic gastroparesis.

BACKGROUND: Cellular changes described in human gastroparesis have revealed a role for immune dysregulation, however, a mechanistic understanding of human gastroparesis and the signaling pathways involved are still unclear. METHODS: Diabetic gastroparetics, diabetic non-gastroparetic controls, idiopathic gastroparetics and non-diabetic non-gastroparetic controls underwent full-thickness gastric body biopsies. Deep RNA sequencing was performed and pathway analysis of differentially expressed transcripts was done using Ingenuity®. A subset of differentially expressed genes in diabetic gastroparesis was validated in a separate cohort using QT-PCR. RESULTS: 111 genes were differentially expressed in diabetic gastroparesis and 181 in idiopathic gastroparesis with a log2fold difference of | ≥ 2| and false detection rate (FDR) < 5%. Top canonical pathways in diabetic gastroparesis included genes involved with macrophages, fibroblasts and endothelial cells in rheumatoid arthritis, osteoarthritis pathway and differential regulation of cytokine production in macrophages and T helper cells by IL-17A and IL-17F. Top canonical pathways in idiopathic gastroparesis included genes involved in granulocyte adhesion and diapedesis, agranulocyte adhesion and diapedesis, and role of macrophages, fibroblasts and endothelial cells in rheumatoid arthritis. Sixty-five differentially expressed genes (log2fold difference | ≥ 2|, FDR < 5%) were common in both diabetic and idiopathic gastroparesis with genes in the top 5 canonical pathways associated with immune signaling. 4/5 highly differentially expressed genes (SGK1, APOLD1, CXCR4, CXCL2, and FOS) in diabetic gastroparesis were validated in a separate cohort of patients using RT-PCR. Immune profile analysis revealed that genes associated with M1 (pro inflammatory) macrophages were enriched in tissues from idiopathic gastroparesis tissues compared to controls (p < 0.05). CONCLUSIONS: Diabetic and idiopathic gastroparesis have both unique and overlapping transcriptomic signatures. Innate immune signaling likely plays a central role in pathogenesis of human gastroparesis.

Development and Verification of an RNA Sequencing (RNA-Seq) Assay for the Detection of Gene Fusions in Tumors.

We assessed the performance characteristics of an RNA sequencing (RNA-Seq) assay designed to detect gene fusions in 571 genes to help manage patients with cancer. Polyadenylated RNA was converted to cDNA, which was then used to prepare next-generation sequencing libraries that were sequenced on an Illumina HiSeq 2500 instrument and analyzed with an in-house developed bioinformatic pipeline. The assay identified 38 of 41 gene fusions detected by another method, such as fluorescence in situ hybridization or RT-PCR, for a sensitivity of 93%. No false-positive gene fusions were identified in 15 normal tissue specimens and 10 tumor specimens that were negative for fusions by RNA sequencing or Mate Pair NGS (100% specificity). The assay also identified 22 fusions in 17 tumor specimens that had not been detected by other methods. Eighteen of the 22 fusions had not previously been described. Good intra-assay and interassay reproducibility was observed with complete concordance for the presence or absence of gene fusions in replicates. The analytical sensitivity of the assay was tested by diluting RNA isolated from gene fusion-positive cases with fusion-negative RNA. Gene fusions were generally detectable down to 12.5% dilutions for most fusions and as little as 3% for some fusions. This assay can help identify fusions in patients with cancer; these patients may in turn benefit from both US Food and Drug Administration-approved and investigational targeted therapies.

Validation of Osteogenic Properties of Cytochalasin D by High-Resolution RNA-Sequencing in Mesenchymal Stem Cells Derived from Bone Marrow and Adipose Tissues.

Differentiation of mesenchymal stromal/stem cells (MSCs) involves a series of molecular signals and gene transcription events required for attaining cell lineage commitment. Modulation of the actin cytoskeleton using cytochalasin D (CytoD) drives osteogenesis at early timepoints in bone marrow-derived MSCs and also initiates a robust osteogenic differentiation program in adipose tissue-derived MSCs. To understand the molecular basis for these pronounced effects on osteogenic differentiation, we investigated global changes in gene expression in CytoD-treated murine and human MSCs by high-resolution RNA-sequencing (RNA-seq) analysis. A three-way bioinformatic comparison between human adipose tissue-derived MSCs (hAMSCs), human bone marrow-derived MSCs (hBMSCs), and mouse bone marrow-derived MSCs (mBMSCs) revealed significant upregulation of genes linked to extracellular matrix organization, cell adhesion and bone metabolism. As anticipated, the activation of these differentiation-related genes is accompanied by a downregulation of nuclear and cell cycle-related genes presumably reflecting cytostatic effects of CytoD. We also identified eight novel CytoD activated genes-VGLL4, ARHGAP24, KLHL24, RCBTB2, BDH2, SCARF2, ACAD10, HEPH-which are commonly upregulated across the two species and tissue sources of our MSC samples. We selected the Hippo pathway-related VGLL4 gene, which encodes the transcriptional co-factor Vestigial-like 4, for further study because this pathway is linked to osteogenesis. VGLL4 small interfering RNA depletion reduces mineralization of hAMSCs during CytoD-induced osteogenic differentiation. Together, our RNA-seq analyses suggest that while the stimulatory effects of CytoD on osteogenesis are pleiotropic and depend on the biological state of the cell type, a small group of genes including VGLL4 may contribute to MSC commitment toward the bone lineage.

Gene Expression Signatures Characterized by Longitudinal Stability and Interindividual Variability Delineate Baseline Phenotypic Groups with Distinct Responses to Immune Stimulation.

Human immunity exhibits remarkable heterogeneity among individuals, which engenders variable responses to immune perturbations in human populations. Population studies reveal that, in addition to interindividual heterogeneity, systemic immune signatures display longitudinal stability within individuals, and these signatures may reliably dictate how given individuals respond to immune perturbations. We hypothesize that analyzing relationships among these signatures at the population level may uncover baseline immune phenotypes that correspond with response outcomes to immune stimuli. To test this, we quantified global gene expression in peripheral blood CD4+ cells from healthy individuals at baseline and following CD3/CD28 stimulation at two time points 1 mo apart. Systemic CD4+ cell baseline and poststimulation molecular immune response signatures (MIRS) were defined by identifying genes expressed at levels that were stable between time points within individuals and differential among individuals in each state. Iterative differential gene expression analyses between all possible phenotypic groupings of at least three individuals using the baseline and stimulated MIRS gene sets revealed shared baseline and response phenotypic groupings, indicating the baseline MIRS contained determinants of immune responsiveness. Furthermore, significant numbers of shared phenotype-defining sets of determinants were identified in baseline data across independent healthy cohorts. Combining the cohorts and repeating the analyses resulted in identification of over 6000 baseline immune phenotypic groups, implying that the MIRS concept may be useful in many immune perturbation contexts. These findings demonstrate that patterns in complex gene expression variability can be used to define immune phenotypes and discover determinants of immune responsiveness.

Network-directed cis-mediator analysis of normal prostate tissue expression profiles reveals downstream regulatory associations of prostate cancer susceptibility loci.

Large-scale genome-wide association studies have identified multiple single-nucleotide polymorphisms associated with risk of prostate cancer. Many of these genetic variants are presumed to be regulatory in nature; however, follow-up expression quantitative trait loci (eQTL) association studies have to-date been restricted largely to cis-acting associations due to study limitations. While trans-eQTL scans suffer from high testing dimensionality, recent evidence indicates most trans-eQTL associations are mediated by cis-regulated genes, such as transcription factors. Leveraging a data-driven gene co-expression network, we conducted a comprehensive cis-mediator analysis using RNA-Seq data from 471 normal prostate tissue samples to identify downstream regulatory associations of previously identified prostate cancer risk variants. We discovered multiple trans-eQTL associations that were significantly mediated by cis-regulated transcripts, four of which involved risk locus 17q12, proximal transcription factor HNF1B, and target trans-genes with known HNF response elements (MIA2, SRC, SEMA6A, KIF12). We additionally identified evidence of cis-acting down-regulation of MSMB via rs10993994 corresponding to reduced co-expression of NDRG1. The majority of these cis-mediator relationships demonstrated trans-eQTL replicability in 87 prostate tissue samples from the Gene-Tissue Expression Project. These findings provide further biological context to known risk loci and outline new hypotheses for investigation into the etiology of prostate cancer.

Gastroblastoma harbors a recurrent somatic MALAT1-GLI1 fusion gene.

Gastroblastoma is a rare distinctive biphasic tumor of the stomach. The molecular biology of gastroblastoma has not been studied, and no affirmative diagnostic markers have been developed. We retrieved two gastroblastomas from the consultation practices of the authors and performed transcriptome sequencing on formalin-fixed paraffin-embedded tissue. Recurrent predicted fusion genes were validated at genomic and RNA levels. The presence of the fusion gene was confirmed on two additional paraffin-embedded cases of gastroblastoma. Control cases of histologic mimics (biphasic synovial sarcoma, leiomyoma, leiomyosarcoma, desmoid-type fibromatosis, EWSR1-FLI1-positive Ewing sarcoma, Wilms' tumor, gastrointestinal stromal tumor, plexiform fibromyxoma, Sonic hedgehog-type medulloblastomas, and normal gastric mucosa and muscularis propria were also analyzed. The gastroblastomas affected two males and two females aged 9-56 years. Transcriptome sequencing identified recurrent somatic MALAT1-GLI1 fusion genes, which were predicted to retain the key domains of GLI1. The MALAT1-GLI1 fusion gene was validated by break-apart and dual-fusion FISH and RT-PCR. The additional two gastroblastomas were also positive for the MALAT1-GLI1 fusion gene. None of the other control cases harbored MALAT1-GLI1. Overexpression of GLI1 in the cases of gastroblastomas was confirmed at RNA and protein levels. Pathway analysis revealed activation of the Sonic hedgehog pathway in gastroblastoma and gene expression profiling showed that gastroblastomas grouped together and were most similar to Sonic hedgehog-type medulloblastomas. In summary, we have identified an oncogenic MALAT1-GLI1 fusion gene in all cases of gastroblastoma that may serve as a diagnostic biomarker. The fusion gene is predicted to encode a protein that includes the zinc finger domains of GLI1 and results in overexpression of GLI1 protein and activation of the Sonic hedgehog pathway.