Small cell carcinoma of the ovary of hypercalcemic type (SCCOHT) is an extremely rare, aggressive cancer affecting children and young women. We identified germline and somatic inactivating mutations in the SWI/SNF chromatin-remodeling gene SMARCA4 in 69% (9/13) of SCCOHT cases in addition to SMARCA4 protein loss in 82% (14/17) of SCCOHT tumors but in only 0.4% (2/485) of other primary ovarian tumors. These data implicate SMARCA4 in SCCOHT oncogenesis.
Diffuse large B-cell lymphoma (DLBCL) is an aggressive disease, with 30% to 40% of patients failing to be cured with available primary therapy. microRNAs (miRNAs) are RNA molecules that attenuate expression of their mRNA targets. To characterize the DLBCL miRNome, we sequenced miRNAs from 92 DLBCL and 15 benign centroblast fresh frozen samples and from 140 DLBCL formalin-fixed, paraffin-embedded tissue samples for validation.
We identify known and candidate novel miRNAs, 25 of which are associated with survival independently of cell-of-origin and International Prognostic Index scores, which are established indicators of outcome. Of these 25 miRNAs, six miRNAs are significantly associated with survival in our validation cohort. Abundant expression of miR-28-5p, miR-214-5p, miR-339-3p, and miR-5586-5p is associated with superior outcome, while abundant expression of miR-324-5p and NOVELM00203M is associated with inferior outcome. Comparison of DLBCL miRNA-seq expression profiles with those from other cancer types identifies miRNAs that were more abundant in B-cell contexts. Unsupervised clustering of miRNAs identifies two clusters of patients that have distinct differences in their outcomes. Our integrative miRNA and mRNA expression analyses reveal that miRNAs increased in abundance in DLBCL appear to regulate the expression of genes involved in metabolism, cell cycle, and protein modification. Additionally, these miRNAs, including one candidate novel miRNA, miR-10393-3p, appear to target chromatin modification genes that are frequent targets of somatic mutation in non-Hodgkin lymphomas.
Our comprehensive sequence analysis of the DLBCL miRNome identifies candidate novel miRNAs and miRNAs associated with survival, reinforces results from previous mutational analyses, and reveals regulatory networks of significance for lymphomagenesis.
Electronic supplementary material
The online version of this article (doi:10.1186/s13059-014-0568-y) contains supplementary material, which is available to authorized users.
Diffuse gliomas consist of both low- and high-grade varieties, each with distinct morphological and biological features. The often extended periods of relative indolence exhibited by low-grade gliomas (LGG; WHO grade II) differ sharply from the aggressive, rapidly fatal clinical course of primary glioblastoma (GBM; WHO grade IV). Nevertheless, until recently, the molecular foundations underlying this stark biological contrast between glioma variants remained largely unknown. The discoveries of distinctive and highly recurrent genomic and epigenomic abnormalities in LGG have both informed a more accurate classification scheme and pointed to viable avenues for therapeutic development. As such, the field of neuro-oncology now seems poised to capitalize on these gains to achieve significant benefit for LGG patients. This report will briefly recount the proceedings of a workshop held in January 2013 and hosted by Accelerate Brain Cancer Cure (ABC2) on the subject of LGG. While much of the meeting covered recent insights into LGG biology, its focus remained on how best to advance the clinical management, whether by improved preclinical modeling, more effective targeted therapeutics and clinical trial design, or innovative imaging technology.
clinical trials; genomics; low-grade glioma; personalized medicine
We describe the landscape of somatic genomic alterations based on multi-dimensional and comprehensive characterization of more than 500 glioblastoma tumors (GBMs). We identify several novel mutated genes as well as complex rearrangements of signature receptors including EGFR and PDGFRA. TERT promoter mutations are shown to correlate with elevated mRNA expression, supporting a role in telomerase reactivation. Correlative analyses confirm that the survival advantage of the proneural subtype is conferred by the G-CIMP phenotype, and MGMT DNA methylation may be a predictive biomarker for treatment response only in classical subtype GBM. Integrative analysis of genomic and proteomic profiles challenges the notion of therapeutic inhibition of a pathway as an alternative to inhibition of the target itself. These data will facilitate the discovery of therapeutic and diagnostic target candidates, the validation of research and clinical observations and the generation of unanticipated hypotheses that can advance our molecular understanding of this lethal cancer.
Gastric cancer is a leading cause of cancer deaths, but analysis of its molecular and clinical characteristics has been complicated by histological and aetiological heterogeneity. Here we describe a comprehensive molecular evaluation of 295 primary gastric adenocarcinomas as part of The Cancer Genome Atlas (TCGA) project. We propose a molecular classification dividing gastric cancer into four subtypes: tumours positive for Epstein–Barr virus, which display recurrent PIK3CA mutations, extreme DNA hypermethylation, and amplification of JAK2, CD274 (also known as PD-L1) and PDCD1LG2 (also knownasPD-L2); microsatellite unstable tumours, which show elevated mutation rates, including mutations of genes encoding targetable oncogenic signalling proteins; genomically stable tumours, which are enriched for the diffuse histological variant and mutations of RHOA or fusions involving RHO-family GTPase-activating proteins; and tumours with chromosomal instability, which show marked aneuploidy and focal amplification of receptor tyrosine kinases. Identification of these subtypes provides a roadmap for patient stratification and trials of targeted therapies.
The majority of oligodendrogliomas (ODGs) exhibit combined losses of chromosomes 1p and 19q and mutations of isocitrate dehydrogenase (IDH1-R132H or IDH2-R172K). Approximately 70% of ODGs with 1p19q co-deletions harbor somatic mutations in the Capicua Transcriptional Repressor (CIC) gene on chromosome 19q13.2. Here we show that endogenous long (CIC-L) and short (CIC-S) CIC proteins are predominantly localized to the nucleus or cytoplasm, respectively. Cytoplasmic CIC-S is found in close proximity to the mitochondria. To study wild type and mutant CIC function and motivated by the paucity of 1p19q co-deleted ODG lines, we created HEK293 and HOG stable cell lines ectopically co-expressing CIC and IDH1. Non-mutant lines displayed increased clonogenicity, but cells co-expressing the mutant IDH1-R132H with either CIC-S-R201W or -R1515H showed reduced clonogenicity in an additive manner, demonstrating cooperative effects in our assays. Expression of mutant CIC-R1515H increased cellular 2-Hydroxyglutarate (2HG) levels compared to wild type CIC in IDH1-R132H background. Levels of phosphorylated ATP-citrate Lyase (ACLY) were lower in cell lines expressing mutant CIC-S proteins compared to cells expressing wild type CIC-S, supporting a cytosolic citrate metabolism-related mechanism of reduced clonogenicity in our in vitro model systems. ACLY or phospho-ACLY were similarly reduced in CIC-mutant 1p19q co-deleted oligodendroglioma patient samples.
CIC; IDH1; ACLY; citrate; 2HG; Oligodendroglioma
Desmosterolosis is an autosomal recessive disorder of cholesterol biosynthesis caused by biallelic mutations of DHCR24 (homozygous or compound heterozygous), which encodes 3-β-hydroxysterol Δ-24-reductase. We report two sisters homozygous for the 571G>A (E191K) DHCR24 mutation. Comparison of the propositae to other reported individuals shows that psychomotor developmental delay, failure to thrive, dysgenesis of the corpus callosum, cerebral white matter atrophy and spasticity likely constitute the minimal desmosterolosis phenotype. The nonspecific features of desmosterolosis make it difficult to suspect clinically and therefore screening for it should be entertained early in the diagnostic evaluation.
DHCR24; Desmosterol; Intellectual disability; Cholesterol biosynthesis; Exome sequencing
While the methylation of DNA in 5′ promoters suppresses gene expression, the role of DNA methylation in gene bodies is unclear1–5. In mammals, tissue- and cell type-specific methylation is present in a small percentage of 5′ CpG island (CGI) promoters, while a far greater proportion occurs across gene bodies, coinciding with highly conserved sequences5–10. Tissue-specific intragenic methylation might reduce,3 or, paradoxically, enhance transcription elongation efficiency1,2,4,5. Capped analysis of gene expression (CAGE) experiments also indicate that transcription commonly initiates within and between genes11–15. To investigate the role of intragenic methylation, we generated a map of DNA methylation from human brain encompassing 24.7 million of the 28 million CpG sites. From the dense, high-resolution coverage of CpG islands, the majority of methylated CpG islands were revealed to be in intragenic and intergenic regions, while less than 3% of CpG islands in 5′ promoters were methylated. The CpG islands in all three locations overlapped with RNA markers of transcription initiation, and unmethylated CpG islands also overlapped significantly with trimethylation of H3K4, a histone modification enriched at promoters16. The general and CpG-island-specific patterns of methylation are conserved in mouse tissues. An in-depth investigation of the human SHANK3 locus17,18 and its mouse homologue demonstrated that this tissue-specific DNA methylation regulates intragenic promoter activity in vitro and in vivo. These methylation-regulated, alternative transcripts are expressed in a tissue and cell type-specific manner, and are expressed differentially within a single cell type from distinct brain regions. These results support a major role for intragenic methylation in regulating cell context-specific alternative promoters in gene bodies.
Intragenic DNA methylation; alternate promoters; comparative epigenomics; SHANK3
Tumor recurrence is a leading cause of cancer mortality. Therapies for recurrent disease may fail, at least in part, because the genomic alterations driving the growth of recurrences are distinct from those in the initial tumor. To explore this hypothesis, we sequenced the exomes of 23 initial low-grade gliomas and recurrent tumors resected from the same patients. In 43% of cases, at least half of the mutations in the initial tumor were undetected at recurrence, including driver mutations in TP53, ATRX, SMARCA4, and BRAF, suggesting recurrent tumors are often seeded by cells derived from the initial tumor at a very early stage of their evolution. Notably, tumors from 6 of 10 patients treated with the chemotherapeutic drug temozolomide (TMZ) followed an alternative evolutionary path to high-grade glioma. At recurrence, these tumors were hypermutated and harbored driver mutations in the RB and AKT-mTOR pathways that bore the signature of TMZ-induced mutagenesis.
Cryptococcus neoformans is a pathogenic basidiomycetous yeast responsible for more than 600,000 deaths each year. It occurs as two serotypes (A and D) representing two varieties (i.e. grubii and neoformans, respectively). Here, we sequenced the genome and performed an RNA-Seq-based analysis of the C. neoformans var. grubii transcriptome structure. We determined the chromosomal locations, analyzed the sequence/structural features of the centromeres, and identified origins of replication. The genome was annotated based on automated and manual curation. More than 40,000 introns populating more than 99% of the expressed genes were identified. Although most of these introns are located in the coding DNA sequences (CDS), over 2,000 introns in the untranslated regions (UTRs) were also identified. Poly(A)-containing reads were employed to locate the polyadenylation sites of more than 80% of the genes. Examination of the sequences around these sites revealed a new poly(A)-site-associated motif (AUGHAH). In addition, 1,197 miscRNAs were identified. These miscRNAs can be spliced and/or polyadenylated, but do not appear to have obvious coding capacities. Finally, this genome sequence enabled a comparative analysis of strain H99 variants obtained after laboratory passage. The spectrum of mutations identified provides insights into the genetics underlying the micro-evolution of a laboratory strain, and identifies mutations involved in stress responses, mating efficiency, and virulence.
Cryptococcus neoformans var. grubii is a major human pathogen responsible for deadly meningoencephalitis in immunocompromised patients. Here, we report the sequencing and annotation of its genome. Evidence for extensive intron splicing, antisense transcription, non-coding RNAs, and alternative polyadenylation indicates the potential for highly intricate regulation of gene expression in this opportunistic pathogen. In addition, detailed molecular, genetic, and genomic studies were performed to characterize structural features of the genome, including centromeres and origins of replication. Finally, the phenotypic and genome re-sequencing analysis of a collection of isolates of the reference H99 strain resulting from laboratory passage revealed that microevolutionary processes during in vitro culturing of pathogenic fungi can impact virulence.
Medulloblastoma is comprised of four distinct molecular variants: WNT, SHH, Group 3, and Group 4. We analyzed alternative splicing usage in 14 normal cerebellar samples and 103 medulloblastomas of known subgroup. Medulloblastoma samples have a statistically significant increase in alternative splicing as compared to normal fetal cerebella (2.3-times; P<6.47E-8). Splicing patterns are distinct and specific between molecular subgroups. Unsupervised hierarchical clustering of alternative splicing events accurately assigns medulloblastomas to their correct subgroup. Subgroup-specific splicing and alternative promoter usage was most prevalent in Group 3 (19.4%) and SHH (16.2%) medulloblastomas, while observed less frequently in WNT (3.2%), and Group 4 (9.3%) tumors. Functional annotation of alternatively spliced genes reveals over-representation of genes important for neuronal development. Alternative splicing events in medulloblastoma may be regulated in part by the correlative expression of antisense transcripts, suggesting a possible mechanism affecting subgroup specific alternative splicing. Our results identify additional candidate markers for medulloblastoma subgroup affiliation, further support the existence of distinct subgroups of the disease, and demonstrate an additional level of transcriptional heterogeneity between medulloblastoma subgroups.
medulloblastoma; alternative splicing; neuronal development; molecular subgroup; pediatric cancer
Leukemia stem cells (LSC) play a pivotal role in chronic myeloid leukemia (CML) tyrosine kinase inhibitor (TKI) resistance and progression to blast crisis (BC), in part, through alternative splicing of self-renewal and survival genes. To elucidate splice isoform regulators of human BC LSC maintenance, we performed whole transcriptome RNA sequencing; splice isoform-specific qRT-PCR, nanoproteomics, stromal co-culture and BC LSC xenotransplantation analyses. Cumulatively, these studies show that alternative splicing of multiple pro-survival BCL2 family genes promotes malignant transformation of myeloid progenitors into BC LSC that are quiescent in the marrow niche and contribute to therapeutic resistance. Notably, a novel pan-BCL2 inhibitor, sabutoclax, renders marrow niche-resident BC LSC sensitive to TKIs at doses that spare normal progenitors. These findings underscore the importance of alternative BCL2 family splice isoform expression in BC LSC maintenance and suggest that combinatorial inhibition of pro-survival BCL2 family proteins and BCR-ABL may eliminate dormant LSC and obviate resistance.
Pathways defining susceptibility of normal cells to oncogenic transformation may be valuable therapeutic targets. We characterized the cell of origin and its critical pathways in MN1-induced leukemias. Common myeloid (CMP), but not granulocyte-macrophage progenitors (GMP) could be transformed by MN1. Complementation studies of CMP-signature genes in GMPs demonstrated that MN1-leukemogenicity required the MEIS1/AbdB-like HOX-protein complex. ChIP-sequencing identified common target genes of MN1 and MEIS1, and demonstrated identical binding sites for a large proportion of their chromatin targets. Transcriptional repression of MEIS1 targets in established MN1 leukemias demonstrated antileukemic activity. As MN1 relies on but cannot activate expression of MEIS1/AbdB-like HOX proteins, transcriptional activity of these genes determines cellular susceptibility to MN1-induced transformation, and may represent a promising therapeutic target.
PMID: 21741595 CAMSID: cams3759
Recent sequencing efforts have described the mutational landscape of the pediatric brain tumor medulloblastoma. Although MLL2 is among the most frequent somatic single nucleotide variants (SNV), the clinical and biological significance of these mutations remains uncharacterized. Through targeted re-sequencing, we identified mutations of MLL2 in 8 % (14/175) of MBs, the majority of which were loss of function. Notably, we also report mutations affecting the MLL2-binding partner KDM6A, in 4 % (7/175) of tumors. While MLL2 mutations were independent of age, gender, histological subtype, M-stage or molecular subgroup, KDM6A mutations were most commonly identified in Group 4 MBs, and were mutually exclusive with MLL2 mutations. Immunohistochemical staining for H3K4me3 and H3K27me3, the chromatin effectors of MLL2 and KDM6A activity, respectively, demonstrated alterations of the histone code in 24 % (53/220) of MBs across all subgroups. Correlating these MLL2-and KDM6A-driven histone marks with prognosis, we identified populations of MB with improved (K4+/K27−) and dismal (K4−/K27−) outcomes, observed primarily within Group 3 and 4 MBs. Group 3 and 4 MBs demonstrate somatic copy number aberrations, and transcriptional profiles that converge on modifiers of H3K27-methylation (EZH2, KDM6A, KDM6B), leading to silencing of PRC2-target genes. As PRC2-mediated aberrant methylation of H3K27 has recently been targeted for therapy in other diseases, it represents an actionable target for a substantial percentage of medulloblastoma patients with aggressive forms of the disease.
MLL2; KDM6A; Histone lysine methylation; Medulloblastoma; PRC2
Chromosomal translocations are critically involved in the molecular pathogenesis of B-cell lymphomas, and highly recurrent and specific rearrangements have defined distinct molecular subtypes linked to unique clinicopathological features1,2. In contrast, several well-characterized lymphoma entities still lack disease-defining translocation events. To identify novel fusion transcripts resulting from translocations, we investigated two Hodgkin lymphoma cell lines by whole-transcriptome paired-end sequencing (RNA-seq). Here we show a highly expressed gene fusion involving the major histocompatibility complex (MHC) class II transactivator CIITA (MHC2TA) in KM-H2 cells. In a subsequent evaluation of 263 B-cell lymphomas, we also demonstrate that genomic CIITA breaks are highly recurrent in primary mediastinal B-cell lymphoma (38%) and classical Hodgkin lymphoma (cHL) (15%). Furthermore, we find that CIITA is a promiscuous partner of various in-frame gene fusions, and we report that CIITA gene alterations impact survival in primary mediastinal B-cell lymphoma (PMBCL). As functional consequences of CIITA gene fusions, we identify downregulation of surface HLA class II expression and overexpression of ligands of the receptor molecule programmed cell death 1 (CD274/PDL1 and CD273/PDL2). These receptor–ligand interactions have been shown to impact anti-tumour immune responses in several cancers3, whereas decreased MHC class II expression has been linked to reduced tumour cell immunogenicity4. Thus, our findings suggest that recurrent rearrangements of CIITA may represent a novel genetic mechanism underlying tumour–microenvironment interactions across a spectrum of lymphoid cancers.
Transposable element (TE) derived sequences comprise half of our genome and DNA methylome, and are presumed densely methylated and inactive. Examination of the genome-wide DNA methylation status within 928 TE subfamilies in human embryonic and adult tissues revealed unexpected tissue-specific and subfamily-specific hypomethylation signatures. Genes proximal to tissue-specific hypomethylated TE sequences were enriched for functions important for the tissue type and their expression correlated strongly with hypomethylation of the TEs. When hypomethylated, these TE sequences gained tissue-specific enhancer marks including H3K4me1 and occupancy by p300, and a majority exhibited enhancer activity in reporter gene assays. Many such TEs also harbored binding sites for transcription factors that are important for tissue-specific functions and exhibited evidence for evolutionary selection. These data suggest that sequences derived from TEs may be responsible for wiring tissue type-specific regulatory networks, and have acquired tissue-specific epigenetic regulation.
Primary triple negative breast cancers (TNBC) represent approximately 16% of all breast cancers1 and are a tumour type defined by exclusion, for which comprehensive landscapes of somatic mutation have not been determined. Here we show in 104 early TNBC cases, that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some exhibiting only a handful of somatic aberrations in a few pathways, whereas others contain hundreds of somatic events and multiple pathways implicated. Integration with matched whole transcriptome sequence data revealed that only ~36% of mutations are expressed. By examining single nucleotide variant (SNV) allelic abundance derived from deep re-sequencing (median >20,000 fold) measurements in 2414 somatic mutations, we determine for the first time in an epithelial tumour, the relative abundance of clonal genotypes among cases in the population. We show that TNBC vary widely and continuously in their clonal frequencies at the time of diagnosis, with basal subtype TNBC2,3 exhibiting more variation than non-basal TNBC. Although p53 and PIK3CA/PTEN somatic mutations appear clonally dominant compared with other pathways, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal and cell shape/motility proteins occurred at lower clonal frequencies, suggesting they occurred later during tumour progression. Taken together our results show that future attempts to dissect the biology and therapeutic responses of TNBC will require the determination of individual tumour clonal genotypes.
The extent to which a genomic test will be used in practice is affected by factors such as ability of the test to correctly predict response to treatment (i.e. sensitivity and specificity of the test), invasiveness of the testing procedure, test cost, and the probability and severity of side effects associated with treatment.
Using discrete choice experimentation (DCE), we elicited preferences of the public (Sample 1, N = 533 and Sample 2, N = 525) and cancer patients (Sample 3, N = 38) for different attributes of a hypothetical genomic test for guiding cancer treatment. Samples 1 and 3 considered the test/treatment in the context of an aggressive curable cancer (scenario A) while the scenario for sample 2 was based on a non-aggressive incurable cancer (scenario B).
In aggressive curable cancer (scenario A), everything else being equal, the odds ratio (OR) of choosing a test with 95% sensitivity was 1.41 (versus a test with 50% sensitivity) and willingness to pay (WTP) was $1331, on average, for this amount of improvement in test sensitivity. In this scenario, the OR of choosing a test with 95% specificity was 1.24 times that of a test with 50% specificity (WTP = $827). In non-aggressive incurable cancer (scenario B), the OR of choosing a test with 95% sensitivity was 1.65 (WTP = $1344), and the OR of choosing a test with 95% specificity was 1.50 (WTP = $1080). Reducing severity of treatment side effects from severe to mild was associated with large ORs in both scenarios (OR = 2.10 and 2.24 in scenario A and B, respectively). In contrast, patients had a very large preference for 95% sensitivity of the test (OR = 5.23).
The type and prognosis of cancer affected preferences for genomically-guided treatment. In aggressive curable cancer, individuals emphasized more on the sensitivity rather than the specificity of the test. In contrast, for a non-aggressive incurable cancer, individuals put similar emphasis on sensitivity and specificity of the test. While the public expressed strong preference toward lowering severity of side effects, improving sensitivity of the test had by far the largest influence on patients’ decision to use genomic testing.
Pharmacogenomics; Genomic medicine; Personalized medicine; Genetic testing; Discrete choice experiment; Conjoint analysis; Preference elicitation; Cancer treatment
Next-generation sequencing is making sequence-based molecular pathology and personalized oncology viable. We selected an individual initially diagnosed with conventional but aggressive prostate adenocarcinoma and sequenced the genome and transcriptome from primary and metastatic tissues collected prior to hormone therapy. The histology-pathology and copy number profiles were remarkably homogeneous, yet it was possible to propose the quadrant of the prostate tumour that likely seeded the metastatic diaspora. Despite a homogeneous cell type, our transcriptome analysis revealed signatures of both luminal and neuroendocrine cell types. Remarkably, the repertoire of expressed but apparently private gene fusions, including C15orf21:MYC, recapitulated this biology. We hypothesize that the amplification and over-expression of the stem cell gene MSI2 may have contributed to the stable hybrid cellular identity. This hybrid luminal-neuroendocrine tumour appears to represent a novel and highly aggressive case of prostate cancer with unique biological features and, conceivably, a propensity for rapid progression to castrate-resistance. Overall, this work highlights the importance of integrated analyses of genome, exome and transcriptome sequences for basic tumour biology, sequence-based molecular pathology and personalized oncology.
RNA sequencing; DNA sequencing; prostate cancer; fusion genes; neuroendocrine; personalized medicine; cancer genetics
Neuroblastoma is a malignancy of the developing sympathetic nervous system that often presents with widespread metastatic disease, resulting in survival rates of less than 50%1. To determine the spectrum of somatic mutation in high-risk neuroblastoma, we studied 240 cases using a combination of whole exome, genome and transcriptome sequencing as part of the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiative. Here we report a low median exonic mutation frequency of 0.60 per megabase (0.48 non-silent), and remarkably few recurrently mutated genes in these tumors. Genes with significant somatic mutation frequencies included ALK (9.2% of cases), PTPN11 (2.9%), ATRX (2.5%, an additional 7.1% had focal deletions), MYCN (1.7%, a recurrent p.Pro44Leu alteration), and NRAS (0.83%). Rare, potentially pathogenic germline variants were significantly enriched in ALK, CHEK2, PINK1, and BARD1. The relative paucity of recurrent somatic mutations in neuroblastoma challenges current therapeutic strategies reliant upon frequently altered oncogenic drivers.
Genomic profiling has identified a subtype of high-risk B-progenitor acute lymphoblastic leukemia (B-ALL) with alteration of IKZF1, a gene expression profile similar to BCR-ABL1-positive ALL and poor outcome (Ph-like ALL). The genetic alterations that activate kinase signaling in Ph-like ALL are poorly understood. We performed transcriptome and whole genome sequencing on 15 cases of Ph-like ALL, and identified rearrangements involving ABL1, JAK2, PDGFRB, CRLF2 and EPOR, activating mutations of IL7R and FLT3, and deletion of SH2B3, which encodes the JAK2 negative regulator LNK. Importantly, several of these alterations induce transformation that is attenuated with tyrosine kinase inhibitors, suggesting the treatment outcome of these patients may be improved with targeted therapy.
Summary: Despite recent progress, computational tools that identify gene fusions from next-generation whole transcriptome sequencing data are often limited in accuracy and scalability. Here, we present a software package, BreakFusion that combines the strength of reference alignment followed by read-pair analysis and de novo assembly to achieve a good balance in sensitivity, specificity and computational efficiency.
Supplementary data are available at Bioinformatics online
Castrate resistant prostate cancer (CRPC) and neuroendocrine carcinoma of the prostate are invariably fatal diseases for which only palliative therapies exist. As part of a prostate tumour sequencing program, a patient tumour was analyzed using Illumina genome sequencing and a matched renal capsule tumour xenograft was generated. Both tumour and xenograft had a homozygous 9p21 deletion spanning the MTAP, CDKN2 and ARF genes. It is rare for this deletion to occur in primary prostate tumours yet approximately 10% express decreased levels of MTAP mRNA. Decreased MTAP expression is a prognosticator for poor outcome. Moreover, it appears that this deletion is more common in CRPC than in primary prostate cancer. We show for the first time that treatment with methylthioadenosine and high dose 6-thioguanine causes marked inhibition of a patient derived neuroendocrine xenograft growth while protecting the host from 6-thioguanine toxicity. This therapeutic approach can be applied to other MTAP-deficient human cancers since deletion or hypermethylation of the MTAP gene occurs in a broad spectrum of tumours at high frequency. The combination of genome sequencing and patient-derived xenografts can identify candidate therapeutic agents and evaluate them for personalized oncology.
massively parallel sequencing; MTAP; patient-derived xenograft; genitourinary cancers: prostate; animal models of cancer; gene expression profiling; functional genomics; xenograft models
Medulloblastoma, the most common malignant pediatric brain tumour, is currently treated with non-specific cytotoxic therapies including surgery, whole brain radiation, and aggressive chemotherapy. As medulloblastoma exhibits marked intertumoural heterogeneity, with at least four distinct molecular variants, prior attempts to identify targets for therapy have been underpowered due to small samples sizes. Here we report somatic copy number aberrations (SCNAs) in 1087 unique medulloblastomas. SCNAs are common in medulloblastoma, and are predominantly subgroup enriched. The most common region of focal copy number gain is a tandem duplication of the Parkinson’s disease gene SNCAIP, which is exquisitely restricted to Group 4α. Recurrent translocations of PVT1, including PVT1-MYC and PVT1-NDRG1 that arise through chromothripsis are restricted to Group 3. Numerous targetable SCNAs, including recurrent events targeting TGFβ signaling in Group 3, and NF-κB signaling in Group 4 suggest future avenues for rational, targeted therapy.
Biallelic mutations of the DNA annealing helicase SMARCAL1 (SWI/SNF-related, matrix-associated, actin-dependent regulator of chromatin, subfamily a-like 1) cause Schimke immuno-osseous dysplasia (SIOD, MIM 242900), an incompletely penetrant autosomal recessive disorder. Using human, Drosophila and mouse models, we show that the proteins encoded by SMARCAL1 orthologs localize to transcriptionally active chromatin and modulate gene expression. We also show that, as found in SIOD patients, deficiency of the SMARCAL1 orthologs alone is insufficient to cause disease in fruit flies and mice, although such deficiency causes modest diffuse alterations in gene expression. Rather, disease manifests when SMARCAL1 deficiency interacts with genetic and environmental factors that further alter gene expression. We conclude that the SMARCAL1 annealing helicase buffers fluctuations in gene expression and that alterations in gene expression contribute to the penetrance of SIOD.