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.

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.

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.

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.

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.

The Role of the Histone Methyltransferase Enhancer of Zeste Homolog 2 (EZH2) in the Pathobiological Mechanisms Underlying Inflammatory Bowel Disease (IBD).

Regulatory T (Treg) cells expressing the transcription factor FOXP3 play a pivotal role in maintaining immunologic self-tolerance. We and others have shown previously that EZH2 is recruited to the FOXP3 promoter and its targets in Treg cells. To further address the role for EZH2 in Treg cellular function, we have now generated mice that lack EZH2 specifically in Treg cells (EZH2Δ/ΔFOXP3+). We find that EZH2 deficiency in FOXP3+ T cells results in lethal multiorgan autoimmunity. We further demonstrate that EZH2Δ/ΔFOXP3+ T cells lack a regulatory phenotype in vitro and secrete proinflammatory cytokines. Of special interest, EZH2Δ/ΔFOXP3+ mice develop spontaneous inflammatory bowel disease. Guided by these results, we assessed the FOXP3 and EZH2 gene networks by RNA sequencing in isolated intestinal CD4+ T cells from patients with Crohn's disease. Gene network analysis demonstrates that these CD4+ T cells display a Th1/Th17-like phenotype with an enrichment of gene targets shared by FOXP3 and EZH2. Combined, these results suggest that the inflammatory milieu found in Crohn's disease could lead to or result from deregulation of FOXP3/EZH2-enforced T cell gene networks contributing to the underlying intestinal inflammation.

Comprehensively evaluating cis-regulatory variation in the human prostate transcriptome by using gene-level allele-specific expression.

The identification of cis-acting regulatory variation in primary tissues has the potential to elucidate the genetic basis of complex traits and further our understanding of transcriptomic diversity across cell types. Expression quantitative trait locus (eQTL) association analysis using RNA sequencing (RNA-seq) data can improve upon the detection of cis-acting regulatory variation by leveraging allele-specific expression (ASE) patterns in association analysis. Here, we present a comprehensive evaluation of cis-acting eQTLs by analyzing RNA-seq gene-expression data and genome-wide high-density genotypes from 471 samples of normal primary prostate tissue. Using statistical models that integrate ASE information, we identified extensive cis-eQTLs across the prostate transcriptome and found that approximately 70% of expressed genes corresponded to a significant eQTL at a gene-level false-discovery rate of 0.05. Overall, cis-eQTLs were heavily concentrated near the transcription start and stop sites of affected genes, and effects were negatively correlated with distance. We identified multiple instances of cis-acting co-regulation by using phased genotype data and discovered 233 SNPs as the most strongly associated eQTLs for more than one gene. We also noted significant enrichment (25/50, p = 2E-5) of previously reported prostate cancer risk SNPs in prostate eQTLs. Our results illustrate the benefit of assessing ASE data in cis-eQTL analyses by showing better reproducibility of prior eQTL findings than of eQTL mapping based on total expression alone. Altogether, our analysis provides extensive functional context of thousands of SNPs in prostate tissue, and these results will be of critical value in guiding studies examining disease of the human prostate.

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.

Integrated genomic characterization reveals novel, therapeutically relevant drug targets in FGFR and EGFR pathways in sporadic intrahepatic cholangiocarcinoma.

Advanced cholangiocarcinoma continues to harbor a difficult prognosis and therapeutic options have been limited. During the course of a clinical trial of whole genomic sequencing seeking druggable targets, we examined six patients with advanced cholangiocarcinoma. Integrated genome-wide and whole transcriptome sequence analyses were performed on tumors from six patients with advanced, sporadic intrahepatic cholangiocarcinoma (SIC) to identify potential therapeutically actionable events. Among the somatic events captured in our analysis, we uncovered two novel therapeutically relevant genomic contexts that when acted upon, resulted in preliminary evidence of anti-tumor activity. Genome-wide structural analysis of sequence data revealed recurrent translocation events involving the FGFR2 locus in three of six assessed patients. These observations and supporting evidence triggered the use of FGFR inhibitors in these patients. In one example, preliminary anti-tumor activity of pazopanib (in vitro FGFR2 IC50≈350 nM) was noted in a patient with an FGFR2-TACC3 fusion. After progression on pazopanib, the same patient also had stable disease on ponatinib, a pan-FGFR inhibitor (in vitro, FGFR2 IC50≈8 nM). In an independent non-FGFR2 translocation patient, exome and transcriptome analysis revealed an allele specific somatic nonsense mutation (E384X) in ERRFI1, a direct negative regulator of EGFR activation. Rapid and robust disease regression was noted in this ERRFI1 inactivated tumor when treated with erlotinib, an EGFR kinase inhibitor. FGFR2 fusions and ERRFI mutations may represent novel targets in sporadic intrahepatic cholangiocarcinoma and trials should be characterized in larger cohorts of patients with these aberrations.

Impact of RNA degradation on fusion detection by RNA-seq.

BACKGROUND: RNA-seq is a well-established method for studying the transcriptome. Popular methods for library preparation in RNA-seq such as Illumina TruSeq® RNA v2 kit use a poly-A pulldown strategy. Such methods can cause loss of coverage at the 5' end of genes, impacting the ability to detect fusions when used on degraded samples. The goal of this study was to quantify the effects RNA degradation has on fusion detection when using poly-A selected mRNA and to identify the variables involved in this process. RESULTS: Using both artificially and naturally degraded samples, we found that there is a reduced ability to detect fusions as the distance of the breakpoint from the 3' end of the gene increases. The median transcript coverage decreases exponentially as a function of the distance from the 3' end and there is a linear relationship between the coverage decay rate and the RNA integrity number (RIN). Based on these findings we developed plots that show the probability of detecting a gene fusion ("sensitivity") as a function of the distance of the fusion breakpoint from the 3' end. CONCLUSIONS: This study developed a strategy to assess the impact that RNA degradation has on the ability to detect gene fusions by RNA-seq.

Brain expression genome-wide association study (eGWAS) identifies human disease-associated variants.

Genetic variants that modify brain gene expression may also influence risk for human diseases. We measured expression levels of 24,526 transcripts in brain samples from the cerebellum and temporal cortex of autopsied subjects with Alzheimer's disease (AD, cerebellar n=197, temporal cortex n=202) and with other brain pathologies (non-AD, cerebellar n=177, temporal cortex n=197). We conducted an expression genome-wide association study (eGWAS) using 213,528 cisSNPs within ± 100 kb of the tested transcripts. We identified 2,980 cerebellar cisSNP/transcript level associations (2,596 unique cisSNPs) significant in both ADs and non-ADs (q<0.05, p=7.70 × 10(-5)-1.67 × 10(-82)). Of these, 2,089 were also significant in the temporal cortex (p=1.85 × 10(-5)-1.70 × 10(-141)). The top cerebellar cisSNPs had 2.4-fold enrichment for human disease-associated variants (p<10(-6)). We identified novel cisSNP/transcript associations for human disease-associated variants, including progressive supranuclear palsy SLCO1A2/rs11568563, Parkinson's disease (PD) MMRN1/rs6532197, Paget's disease OPTN/rs1561570; and we confirmed others, including PD MAPT/rs242557, systemic lupus erythematosus and ulcerative colitis IRF5/rs4728142, and type 1 diabetes mellitus RPS26/rs1701704. In our eGWAS, there was 2.9-3.3 fold enrichment (p<10(-6)) of significant cisSNPs with suggestive AD-risk association (p<10(-3)) in the Alzheimer's Disease Genetics Consortium GWAS. These results demonstrate the significant contributions of genetic factors to human brain gene expression, which are reliably detected across different brain regions and pathologies. The significant enrichment of brain cisSNPs among disease-associated variants advocates gene expression changes as a mechanism for many central nervous system (CNS) and non-CNS diseases. Combined assessment of expression and disease GWAS may provide complementary information in discovery of human disease variants with functional implications. Our findings have implications for the design and interpretation of eGWAS in general and the use of brain expression quantitative trait loci in the study of human disease genetics.

MAP-RSeq: Mayo Analysis Pipeline for RNA sequencing.

BACKGROUND: Although the costs of next generation sequencing technology have decreased over the past years, there is still a lack of simple-to-use applications, for a comprehensive analysis of RNA sequencing data. There is no one-stop shop for transcriptomic genomics. We have developed MAP-RSeq, a comprehensive computational workflow that can be used for obtaining genomic features from transcriptomic sequencing data, for any genome. RESULTS: For optimization of tools and parameters, MAP-RSeq was validated using both simulated and real datasets. MAP-RSeq workflow consists of six major modules such as alignment of reads, quality assessment of reads, gene expression assessment and exon read counting, identification of expressed single nucleotide variants (SNVs), detection of fusion transcripts, summarization of transcriptomics data and final report. This workflow is available for Human transcriptome analysis and can be easily adapted and used for other genomes. Several clinical and research projects at the Mayo Clinic have applied the MAP-RSeq workflow for RNA-Seq studies. The results from MAP-RSeq have thus far enabled clinicians and researchers to understand the transcriptomic landscape of diseases for better diagnosis and treatment of patients. CONCLUSIONS: Our software provides gene counts, exon counts, fusion candidates, expressed single nucleotide variants, mapping statistics, visualizations, and a detailed research data report for RNA-Seq. The workflow can be executed on a standalone virtual machine or on a parallel Sun Grid Engine cluster. The software can be downloaded from http://bioinformaticstools.mayo.edu/research/maprseq/.

High-resolution molecular validation of self-renewal and spontaneous differentiation in clinical-grade adipose-tissue derived human mesenchymal stem cells.

Improving the effectiveness of adipose-tissue derived human mesenchymal stromal/stem cells (AMSCs) for skeletal therapies requires a detailed characterization of mechanisms supporting cell proliferation and multi-potency. We investigated the molecular phenotype of AMSCs that were either actively proliferating in platelet lysate or in a basal non-proliferative state. Flow cytometry combined with high-throughput RNA sequencing (RNASeq) and RT-qPCR analyses validate that AMSCs express classic mesenchymal cell surface markers (e.g., CD44, CD73/NT5E, CD90/THY1, and CD105/ENG). Expression of CD90 is selectively elevated at confluence. Self-renewing AMSCs express a standard cell cycle program that successively mediates DNA replication, chromatin packaging, cyto-architectural enlargement, and mitotic division. Confluent AMSCs preferentially express genes involved in extracellular matrix (ECM) formation and cellular communication. For example, cell cycle-related biomarkers (e.g., cyclins E2 and B2, transcription factor E2F1) and histone-related genes (e.g., H4, HINFP, NPAT) are elevated in proliferating AMSCs, while ECM genes are strongly upregulated (>10-fold) in quiescent AMSCs. AMSCs also express pluripotency genes (e.g., POU5F1, NANOG, KLF4) and early mesenchymal markers (e.g., NES, ACTA2) consistent with their multipotent phenotype. Strikingly, AMSCs modulate expression of WNT signaling components and switch production of WNT ligands (from WNT5A/WNT5B/WNT7B to WNT2/WNT2B), while upregulating WNT-related genes (WISP2, SFRP2, and SFRP4). Furthermore, post-proliferative AMSCs spontaneously express fibroblastic, osteogenic, chondrogenic, and adipogenic biomarkers when maintained in confluent cultures. Our findings validate the biological properties of self-renewing and multi-potent AMSCs by providing high-resolution quality control data that support their clinical versatility.

RNA sequencing shows transcriptomic changes in rectosigmoid mucosa in patients with irritable bowel syndrome-diarrhea: a pilot case-control study.

Our aim was to conduct a pilot case-control study of RNA expression profile using RNA sequencing of rectosigmoid mucosa of nine females with -diarrhea-predominant irritable bowel syndrome (IBS-D) with accelerated colonic transit and nine female healthy controls. Mucosal total RNA was isolated and purified, and next-generation pair-end sequencing was performed using Illumina TruSeq. Analysis was carried out using a targeted approach toward 12 genes previously associated with IBS and a hypothesis-generating approach. Of the 12 targeted genes tested, patients with IBS-D had decreased mRNA expression of TNFSF15 (fold change controls to IBS-D: 1.53, P = 0.01). Overall, up- and downregulated mRNA expressions of 21 genes (P = 10(-5) to 10(-8); P values with false detection rates are shown) were potentially relevant to IBS-D including the following: neurotransmitters [P2RY4 (P = 0.001), vasoactive intestinal peptide (VIP, P = 0.02)]; cytokines [CCL20 (P = 0.019)]; immune function [C4BPA complement cascade (P = 0.0187)]; interferon-related [IFIT3 (P = 0.016)]; mucosal repair and cell adhesion [trefoil protein (TFF1, P = 0.012)], retinol binding protein [RBP2 (P = 0.017)]; fibronectin (FN1, P = 0.009); and ion channel functions [guanylate cyclase (GUCA2B, P = 0.017), PDZ domain-containing protein 3 (PDZD3, P = 0.029)]. Ten genes associated with functions related to pathobiology of IBS-D were validated by RT-PCR. There was significant correlation in fold changes of the selected genes (Rs = 0.73, P = 0.013). Up- or downregulation of P2RY4, GUC2AB, RBP2, FNI, and C4BPA genes were confirmed on RT-PCR, which also revealed upregulation of farnesoid X receptor (FXR) and apical sodium-coupled bile acid transporter (IBAT/ASBT). RNA-Seq and RT-PCR analysis of rectosigmoid mucosa in IBS-D show transcriptome changes that provide the rationale for validation studies to explore the role of mucosal factors in the pathobiology of IBS-D.

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.

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.

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.

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.