Search

The factors driving therapy resistance in diffuse glioma remain poorly understood. To identify treatment-associated cellular and genetic changes, we analyzed RNA and/or DNA sequencing data from the temporally separated tumor pairs of 304 adult patients with isocitrate dehydrogenase (IDH)-wild-type and IDH-mutant glioma. Tumors recurred in distinct manners that were dependent on IDH mutation status and attributable to changes in histological feature composition, somatic alterations, and microenvironment interactions. Hypermutation and acquired CDKN2A deletions were associated with an increase in proliferating neoplastic cells at recurrence in both glioma subtypes, reflecting active tumor growth. IDH-wild-type tumors were more invasive at recurrence, and their neoplastic cells exhibited increased expression of neuronal signaling programs that reflected a possible role for neuronal interactions in promoting glioma progression. Mesenchymal transition was associated with the presence of a myeloid cell state defined by specific ligand-receptor interactions with neoplastic cells. Collectively, these recurrence-associated phenotypes represent potential targets to alter disease progression.

Technological and computational advances in genomics and interactomics have made it possible to identify how disease mutations perturb protein–protein interaction (PPI) networks within human cells. Here, we show that disease-associated germline variants are significantly enriched in sequences encoding PPI interfaces compared to variants identified in healthy participants from the projects 1000 Genomes and ExAC. Somatic missense mutations are also significantly enriched in PPI interfaces compared to noninterfaces in 10,861 tumor exomes. We computationally identified 470 putative oncoPPIs in a pan-cancer analysis and demonstrate that oncoPPIs are highly correlated with patient survival and drug resistance/sensitivity. We experimentally validate the network effects of 13 oncoPPIs using a systematic binary interaction assay, and also demonstrate the functional consequences of two of these on tumor cell growth. In summary, this human interactome network framework provides a powerful tool for prioritization of alleles with PPI-perturbing mutations to inform pathobiological mechanism- and genotype-based therapeutic discovery. Human disease mutations affect protein–protein interfaces in a three-dimensional structurally resolved interaction network. Predicted oncoPPIs in cancer correlate with survival and drug sensitivity, and affect growth in vitro, supporting their relevance to disease pathogenesis.

Viral recombination is a major evolutionary mechanism driving adaptation processes, such as the ability of host-switching. Understanding global patterns of recombination could help to identify underlying mechanisms and to evaluate the potential risks of rapid adaptation. Conventional approaches (e.g., those based on linkage disequilibrium) are computationally demanding or even intractable when sequence alignments include hundreds of sequences, common in viral data sets. We present a comprehensive analysis of recombination across 30 genomic alignments from viruses infecting humans. In order to scale the analysis and avoid the computational limitations of conventional approaches, we apply newly developed topological data analysis methods able to infer recombination rates for large data sets. We show that viruses, such as ZEBOV and MARV, consistently displayed low levels of recombination, whereas high levels of recombination were observed in Sarbecoviruses, HBV, HEV, Rhinovirus A, and HIV. We observe that recombination is more common in positive single-stranded RNA viruses than in negatively single-stranded RNA ones. Interestingly, the comparison across multiple viruses suggests an inverse correlation between genome length and recombination rate. Positional analyses of recombination breakpoints along viral genomes, combined with our approach, detected at least 39 nonuniform patterns of recombination (i.e., cold or hotspots) in 18 viral groups. Among these, noteworthy hotspots are found in MERS-CoV and Sarbecoviruses (at spike, Nucleocapsid and ORF8). In summary, we have developed a fast pipeline to measure recombination that, combined with other approaches, has allowed us to find both common and lineage-specific patterns of recombination among viruses with potential relevance in viral adaptation.

The gene encoding the core spliceosomal protein SF3B1 is the most frequently mutated gene encoding a splicing factor in a variety of hematologic malignancies and solid tumors. SF3B1 mutations induce use of cryptic 3' splice sites (3'ss), and these splicing errors contribute to tumorigenesis. However, it is unclear how widespread this type of cryptic 3'ss usage is in cancers and what is the full spectrum of genetic mutations that cause such missplicing. To address this issue, we performed an unbiased pan-cancer analysis to identify genetic alterations that lead to the same aberrant splicing as observed with SF3B1 mutations. This analysis identified multiple mutations in another spliceosomal gene, SUGP1, that correlated with significant usage of cryptic 3'ss known to be utilized in mutant SF3B1 expressing cells. Remarkably, this is consistent with recent biochemical studies that identified a defective interaction between mutant SF3B1 and SUGP1 as the molecular defect responsible for cryptic 3'ss usage. Experimental validation revealed that five different SUGP1 mutations completely or partially recapitulated the 3'ss defects. Our analysis suggests that SUGP1 mutations in cancers can induce missplicing identical or similar to that observed in mutant SF3B1 cancers.

Background Macrophages are the most common infiltrating immune cells in gliomas and play a wide variety of pro-tumor and anti-tumor roles. However, the different subpopulations of macrophages and their effects on the tumor microenvironment remain poorly understood. Methods We combined new and previously published single-cell RNA-seq data from 98,015 single cells from a total of 66 gliomas to profile 19,331 individual macrophages. Results Unsupervised clustering revealed a pro-tumor subpopulation of bone marrow-derived macrophages characterized by the scavenger receptor MARCO, which is almost exclusively found in IDH1-wild-type glioblastomas. Previous studies have implicated MARCO as an unfavorable marker in melanoma and non-small cell lung cancer; here, we find that bulk MARCO expression is associated with worse prognosis and mesenchymal subtype. Furthermore, MARCO expression is significantly altered over the course of treatment with anti-PD1 checkpoint inhibitors in a response-dependent manner, which we validate with immunofluorescence imaging. Conclusions These findings illustrate a novel macrophage subpopulation that drives tumor progression in glioblastomas and suggest potential therapeutic targets to prevent their recruitment.

The evolutionary processes that drive universal therapeutic resistance in adult patients with diffuse glioma remain unclear1,2. Here we analysed temporally separated DNA-sequencing data and matched clinical annotation from 222 adult patients with glioma. By analysing mutations and copy numbers across the three major subtypes of diffuse glioma, we found that driver genes detected at the initial stage of disease were retained at recurrence, whereas there was little evidence of recurrence-specific gene alterations. Treatment with alkylating agents resulted in a hypermutator phenotype at different rates across the glioma subtypes, and hypermutation was not associated with differences in overall survival. Acquired aneuploidy was frequently detected in recurrent gliomas and was characterized by IDH mutation but without co-deletion of chromosome arms 1p/19q, and further converged with acquired alterations in the cell cycle and poor outcomes. The clonal architecture of each tumour remained similar over time, but the presence of subclonal selection was associated with decreased survival. Finally, there were no differences in the levels of immunoediting between initial and recurrent gliomas. Collectively, our results suggest that the strongest selective pressures occur during early glioma development and that current therapies shape this evolution in a largely stochastic manner.

The emergence of SARS-CoV-2 underscores the need to better understand the evolutionary processes that drive the emergence and adaptation of zoonotic viruses in humans. In the betacoronavirus genus, which also includes SARS-CoV and MERS-CoV, recombination frequently encompasses the Receptor Binding Domain (RBD) of the Spike protein, which, in turn, is responsible for viral binding to host cell receptors. Here, we find evidence of a recombination event in the RBD involving ancestral linages to both SARS-CoV and SARS-CoV-2. Although we cannot specify the recombinant nor the parental strains, likely due to the ancestry of the event and potential undersampling, our statistical analyses in the space of phylogenetic trees support such an ancestral recombination. Consequently, SARS-CoV and SARS-CoV-2 share an RBD sequence that includes two insertions (positions 432-436 and 460-472), as well as the variants 427N and 436Y. Both 427N and 436Y belong to a helix that interacts directly with the human ACE2 (hACE2) receptor. Reconstruction of ancestral states, combined with protein-binding affinity analyses using the physics-based trRosetta algorithm, reveal that the recombination event involving ancestral strains of SARS-CoV and SARS-CoV-2 led to an increased affinity for hACE2 binding, and that alleles 427N and 436Y significantly enhanced affinity as well. Structural modeling indicates that ancestors of SARS-CoV-2 may have acquired the ability to infect humans decades ago. The binding affinity with the human receptor was subsequently boosted in SARS-CoV and SARS-CoV-2 through further mutations in RBD. In sum, we report an ancestral recombination event affecting the RBD of both SARS-CoV and SARS-CoV-2 that was associated with an increased binding affinity to hACE2.ImportanceThis paper addresses critical questions about the origin of the SARS-CoV-2 virus: what are the evolutionary mechanisms that led to the emergence of the virus, and how can we leverage such knowledge to assess the potential of SARS-like viruses to become pandemic strains? In this work, we demonstrate common mechanisms involved in the emergence of human-infecting SARS-like viruses: first, by acquiring a common haplotype in the RBD through recombination, and further, through increased specificity to the human ACE2 receptor through lineage specific mutations. We also show that the ancestors of SARS-CoV-2 already had the potential to infect humans at least a decade ago, suggesting that SARS-like viruses currently circulating in wild animal species constitute a source of potential pandemic re-emergence.

Glioblastoma is the most frequently occurring and invariably fatal primary brain tumor in adults. The vast majority of glioblastomas is characterized by chromosomal copy number alterations, including gain of whole chromosome 7 and loss of whole chromosome 10. Gain of whole chromosome 7 is an early event in gliomagenesis that occurs in proneural-like precursor cells, which give rise to all isocitrate dehydrogenase (IDH) wild-type glioblastoma transcriptional subtypes. Platelet-derived growth factor A (PDGFA) is one gene on chromosome 7 known to drive gliomagenesis, but, given its location near the end of 7p, there are likely several other genes located along chromosome 7 that select for its increased whole-chromosome copy number within glioblastoma cells. To identify other potential genes that could select for gain of whole chromosome 7, we developed an unbiased bioinformatics approach that identified homeobox A5 (HOXA5) as a gene whose expression correlated with gain of chromosome 7 and a more aggressive phenotype of the resulting glioma. High expression of HOXA5 in glioblastoma was associated with a proneural gene expression pattern and decreased overall survival in both human proneural and PDGF-driven mouse glioblastoma. Furthermore, HOXA5 overexpression promoted cellular proliferation and potentiated radioresistance. We also found enrichment of HOXA5 expression in recurrent human and mouse glioblastoma at first recurrence after radiotherapy. Overall, this study implicates HOXA5 as a chromosome 7-associated gene-level locus that promotes selection for gain of whole chromosome 7 and an aggressive phenotype in glioblastoma.

Relapsed acute lymphoblastic leukaemia (ALL) is associated with resistance to chemotherapy and poor prognosis. Gain-of-function mutations in the 5'-nucleotidase, cytosolic II (NT5C2) gene induce resistance to 6-mercaptopurine and are selectively present in relapsed ALL. Yet, the mechanisms involved in NT5C2 mutation-driven clonal evolution during the initiation of leukaemia, disease progression and relapse remain unknown. Here we use a conditional-and-inducible leukaemia model to demonstrate that expression of NT5C2(R367Q), a highly prevalent relapsed-ALL NT5C2 mutation, induces resistance to chemotherapy with 6-mercaptopurine at the cost of impaired leukaemia cell growth and leukaemia-initiating cell activity. The loss-of-fitness phenotype of NT5C2+/R367Q mutant cells is associated with excess export of purines to the extracellular space and depletion of the intracellular purine-nucleotide pool. Consequently, blocking guanosine synthesis by inhibition of inosine-5'-monophosphate dehydrogenase (IMPDH) induced increased cytotoxicity against NT5C2-mutant leukaemia lymphoblasts. These results identify the fitness cost of NT5C2 mutation and resistance to chemotherapy as key evolutionary drivers that shape clonal evolution in relapsed ALL and support a role for IMPDH inhibition in the treatment of ALL.

Precision medicine in cancer proposes that genomic characterization of tumors can inform personalized targeted therapies. However, this proposition is complicated by spatial and temporal heterogeneity. Here we study genomic and expression profiles across 127 multisector or longitudinal specimens from 52 individuals with glioblastoma (GBM). Using bulk and single-cell data, we find that samples from the same tumor mass share genomic and expression signatures, whereas geographically separated, multifocal tumors and/or long-term recurrent tumors are seeded from different clones. Chemical screening of patient-derived glioma cells (PDCs) shows that therapeutic response is associated with genetic similarity, and multifocal tumors that are enriched with PIK3CA mutations have a heterogeneous drug-response pattern. We show that targeting truncal events is more efficacious than targeting private events in reducing the tumor burden. In summary, this work demonstrates that evolutionary inference from integrated genomic analysis in multisector biopsies can inform targeted therapeutic interventions for patients with GBM.

Sato et al. demonstrate that MYC expression renders MCF10A cells sensitive to down-regulation of UVSSA, as it alleviates stress during MYC-dependent transcription. UVSSA is a promising novel target for attacking MYC-addicted cancer cells., Cancer cells develop strong genetic dependencies, enabling survival under oncogenic stress. MYC is a key oncogene activated across most cancers, and identifying associated synthetic lethality or sickness can provide important clues about its activity and potential therapeutic strategies. On the basis of previously conducted genome-wide screenings in MCF10A cells expressing MYC fused to an estrogen receptor fragment, we identified UVSSA, a gene involved in transcription-coupled repair, whose knockdown or knockout decreased cell viability when combined with MYC expression. Synthetic sick interactions between MYC expression and UVSSA down-regulation correlated with ATM/CHK2 activation, suggesting increased genome instability. We show that the synthetic sick interaction is diminished by attenuating RNA polymerase II (RNAPII) activity; yet, it is independent of UV-induced damage repair, suggesting that UVSSA has a critical function in regulating RNAPII in the absence of exogenous DNA damage. Supporting this hypothesis, RNAPII ChIP-seq revealed that MYC-dependent increases in RNAPII promoter occupancy are reduced or abrogated by UVSSA knockdown, suggesting that UVSSA influences RNAPII dynamics during MYC-dependent transcription. Taken together, our data show that the UVSSA complex has a significant function in supporting MYC-dependent RNAPII dynamics and maintaining cell survival during MYC addiction. While the role of UVSSA in regulating RNAPII has been documented thus far only in the context of UV-induced DNA damage repair, we propose that its activity is also required to cope with transcriptional changes induced by oncogene activation.

Background Recent large-scale cancer sequencing studies have discovered many novel cancer driver genes (CDGs) in human cancers. Some studies also suggest that CDG mutations contribute to cancer-associated epigenomic and transcriptomic alterations across many cancer types. Here we aim to improve our understanding of the connections between CDG mutations and altered cancer cell epigenomes and transcriptomes on pan-cancer level and how these connections contribute to the known association between epigenome and transcriptome. Method Using multi-omics data including somatic mutation, DNA methylation, and gene expression data of 20 cancer types from The Cancer Genome Atlas (TCGA) project, we conducted a pan-cancer analysis to identify CDGs, when mutated, have strong associations with genome-wide methylation or expression changes across cancer types, which we refer as methylation driver genes (MDGs) or expression driver genes (EDGs), respectively. Results We identified 32 MDGs, among which, eight are known chromatin modification or remodeling genes. Many of the remaining 24 MDGs are connected to chromatin regulators through either regulating their transcription or physically interacting with them as potential co-factors. We identified 29 EDGs, 26 of which are also MDGs. Further investigation on target genes’ promoters methylation and expression alteration patterns of these 26 overlapping driver genes shows that hyper-methylation of target genes’ promoters are significantly associated with down-regulation of the same target genes and hypo-methylation of target genes’ promoters are significantly associated with up-regulation of the same target genes. Conclusion This finding suggests a pivotal role for genetically driven changes in chromatin remodeling in shaping DNA methylation and gene expression patterns during tumor development. Electronic supplementary material The online version of this article (10.1186/s12920-018-0425-z) contains supplementary material, which is available to authorized users.

Topoisomerase IIB (TOP2B) is known to decoil and decatenate DNA, releasing chromatin torsional forces. Whereas TOP2B has been implicated in transcript splicing, direct evidence supporting this function, through localization of TOP2B to specific chromatin sites and in association with specific types of cancer, is modest at best. Here, we report our preliminary activities that are intended to address this knowledge gap. GSEA differential expression analysis using TCGA glioblastoma data showed that high TOP2B expression is associated with elevated expression of cell cycle-related genes (p-value= 1.67E(-05)). To investigate the possible association of TOP2B with specific gene sequences we used TOP2B ChIP-seq in analyzing glioma cell lines TS543 and BT142, and found TOP2B associations with the introns of especially long genes (gen length ≥ 100Kb, p-value= 0.001). Functional annotation analysis for genes with TOP2B intron binding showed a significant enrichment for genes with multiple known splice variants (p-value= 1.9E(-11)). To determine effects of TOP2B inhibition on splice variant expression, we treated BT142 cells with the TOP2 inhibitor ICRF-193 and examined treatment effects on splice variant expression using paired-end RNA-seq. We observed altered splicing of 319 genes following TOP2 inhibition, with increased expression of splice variant CDK4-009 being the most significant of all splice variant expression changes (p-value= 1E(-54)). This result corresponded well with TCGA data whose analysis showed a strong positive correlation between TOP2B and CDK4-009 expression (rho= -0.3 and p-value= 7.28E(-05)). CONCLUSIONS: Our results support TOP2B expression and activity as being influential in regulating the levels of CDK4 splice variant 009 in GBM.

Defining the transcriptomic identity of malignant cells is challenging in the absence of surface markers that distinguish cancer clones from one another, or from admixed non-neoplastic cells. To address this challenge, here we developed Genotyping of Transcriptomes (GoT), a method to integrate genotyping with high-throughput droplet-based single-cell RNA sequencing. We apply GoT to profile 38,290 CD34+ cells from patients with CALR-mutated myeloproliferative neoplasms to study how somatic mutations corrupt the complex process of human haematopoiesis. High-resolution mapping of malignant versus normal haematopoietic progenitors revealed an increasing fitness advantage with myeloid differentiation of cells with mutated CALR. We identified the unfolded protein response as a predominant outcome of CALR mutations, with a considerable dependency on cell identity, as well as upregulation of the NF-κB pathway specifically in uncommitted stem cells. We further extended the GoT toolkit to genotype multiple targets and loci that are distant from transcript ends. Together, these findings reveal that the transcriptional output of somatic mutations in myeloproliferative neoplasms is dependent on the native cell identity.

SummaryMeiotic recombination is a fundamental evolutionary process driving diversity in eukaryotes. In mammals, recombination is known to occur preferentially at specific genomic regions. Using topological data analysis (TDA), a branch of applied topology that extracts global features from large data sets, we developed an efficient method for mapping recombination at fine scales. When compared to standard linkage-based methods, TDA can deal with a larger number of SNPs and genomes without incurring prohibitive computational costs. We applied TDA to 1,000 Genomes Project data and constructed high-resolution whole-genome recombination maps of seven human populations. Our analysis shows that recombination is generally under-represented within transcription start sites. However, the binding sites of specific transcription factors are enriched for sites of recombination. These include transcription factors that regulate the expression of meiosis- and gametogenesis-specific genes, cell cycle progression, and differentiation blockage. Additionally, our analysis identifies an enrichment for sites of recombination at repeat-derived loci matched by piwi-interacting RNAs.

Protein synthesis in eukaryotes is regulated by diverse reprogramming mechanisms that expand the coding capacity of individual genes. Here, we exploit one such mechanism, termed -1 programmed ribosomal frameshifting (-1 PRF), to engineer ligand-responsive RNA switches that regulate protein expression. First, efficient -1 PRF stimulatory RNA elements were discovered by in vitro selection; then, ligand-responsive switches were constructed by coupling -1 PRF stimulatory elements to RNA aptamers using rational design and directed evolution in Saccharomyces cerevisiae. We demonstrate that -1 PRF switches tightly control the relative stoichiometry of two distinct protein outputs from a single mRNA, exhibiting consistent ligand response across whole populations of cells. Furthermore, -1 PRF switches were applied to build single-mRNA logic gates and an apoptosis module in yeast. Together, these results showcase the potential for harnessing translation-reprogramming mechanisms for synthetic biology, and they establish -1 PRF switches as powerful RNA tools for controlling protein synthesis in eukaryotes.

Seasonal influenza epidemics and occasional pandemics threaten public health worldwide. New alternative strategies for generating recombinant viruses with vaccine potential are needed. Interestingly, influenza viruses circulating in different hosts have been found to have distinct codon usage patterns, which may reflect host adaptation. We therefore hypothesized that it is possible to make a human seasonal influenza virus that is specifically attenuated in human cells but not in eggs by converting its codon usage so that it is similar to that observed from avian influenza viruses. This approach might help to generate human live attenuated viruses without affecting their yield in eggs. To test this hypothesis, over 300 silent mutations were introduced into the genome of a seasonal H1N1 influenza virus. The resultant mutant was significantly attenuated in mammalian cells and mice, yet it grew well in embryonated eggs. A single dose of intranasal vaccination induced potent innate, humoral, and cellular immune responses, and the mutant could protect mice against homologous and heterologous viral challenges. The attenuated mutant could also be used as a vaccine master donor strain by introducing hemagglutinin and neuraminidase genes derived from other strains. Thus, our approach is a successful strategy to generate attenuated viruses for future application as vaccines. IMPORTANCE Vaccination has been one of the best protective measures in combating influenza virus infection. Current licensed influenza vaccines and their production have various limitations. Our virus attenuation strategy makes use of the codon usage biases of human and avian influenza viruses to generate a human-derived influenza virus that is attenuated in mammalian hosts. This method, however, does not affect virus replication in eggs. This makes the resultant mutants highly compatible with existing egg-based vaccine production pipelines. The viral proteins generated from the codon bias mutants are identical to the wild-type viral proteins. In addition, our massive genome-wide mutational approach further minimizes the concern over reverse mutations. The potential use of this kind of codon bias mutant as a master donor strain to generate other live attenuated viruses is also demonstrated. These findings put forward a promising live attenuated influenza vaccine generation strategy to control influenza.

Outcomes of anticancer therapy vary dramatically among patients due to diverse genetic and molecular backgrounds, highlighting extensive inter-tumoral heterogeneity. The fundamental tenet of precision oncology defines molecular characterization of tumors to guide optimal patient-tailored therapy. Towards this goal, we have established a compilation of pharmacological landscape of 462 patient-derived tumor cells (PDCs) across 14 cancer types, together with genomic and transcriptomic profiling in 385 of these tumors. Compared to the traditional long-term cultured cancer cell-line models, PDCs recapitulate the molecular properties and biology of the diseases more precisely. In present study, we provide unprecedented insights into dynamic pharmacogenomic associations, including molecular determinants that elicit therapeutic resistance to EGFR inhibitors and potential repurposing of ibrutinib (currently being used in hematological malignancies) for EGFR-specific therapy in gliomas. Lastly, we present potential implementation of PDC-derived drug sensitivities for prediction of clinical response to targeted therapeutics using retrospective clinical studies.

Cancer is caused by germline and somatic mutations, which can share biological features such as amino acid change. However, integrated germline and somatic analysis remains uncommon. We present a framework that uses machine learning to learn features of recurrent somatic mutations to (1) predict somatic variants from tumor-only samples and (2) identify somatic-like germline variants for integrated analysis of tumor-normal DNA. Using data from 1769 patients from seven cancer types (bladder, glioblastoma, low-grade glioma, lung, melanoma, stomach, and pediatric glioma), we show that “somatic-like” germline variants are enriched for autosomal-dominant cancer-predisposition genes (p, Cancer: Analytic tool reveals inherited and non-inherited tumor-causing alterations Bladder cancer cells often harbor DNA mutations that occur after tumor development, including some mutations that affect DNA repair. Raul Rabadan, Jiguang Wang, and colleagues from Columbia University in New York, USA, developed an analytic framework for identifying genetic variants, both inherited and newly arisen, that contribute to tumor development. The machine-learning tool—known as Tumor-Only Boosting Identification, or TOBI—learns what’s a cancer-associated mutation from a small training set of tumor samples and matched healthy controls. The researchers then fed the algorithm data from the tumors of 1769 patients with cancers of the bladder, brain, lungs, stomach and skin. They found that TOBI pinpointed many inherited and non-inherited mutations known to contribute to cancer growth. In bladder cancer samples, the tool also revealed a previously unknown role for inherited mutations in BRCA2 and other DNA repair genes in the so-called Fanconi anemia pathway.

Bats harbor many viruses asymptomatically, including several notorious for causing extreme virulence in humans. To identify differences between antiviral mechanisms in humans and bats, we sequenced, assembled, and analyzed the genome of Rousettus aegyptiacus, a natural reservoir of Marburg virus and the only known reservoir for any filovirus. We found an expanded and diversified KLRC/KLRD family of natural killer cell receptors, MHC class I genes, and type I interferons, which dramatically differ from their functional counterparts in other mammals. Such concerted evolution of key components of bat immunity is strongly suggestive of novel modes of antiviral defense. An evaluation of the theoretical function of these genes suggests that an inhibitory immune state may exist in bats. Based on our findings, we hypothesize that tolerance of viral infection, rather than enhanced potency of antiviral defenses, may be a key mechanism by which bats asymptomatically host viruses that are pathogenic in humans.

Both conventional T (Tconv) cells and regulatory T (Treg) cells are activated through ligation of the T cell receptor (TCR) complex, leading to the induction of the transcription factor NF-κB. In Tconv cells, NF-κB regulates expression of genes essential for T cell activation, proliferation, and function. However the role of NF-κB in Treg function remains unclear. We conditionally deleted canonical NF-κB members p65 and c-Rel in developing and mature Treg cells and found they have unique but partially redundant roles. c-Rel was critical for thymic Treg development while p65 was essential for mature Treg identity and maintenance of immune tolerance. Transcriptome and NF-κB p65 binding analyses demonstrated a lineage specific, NF-κB-dependent transcriptional program, enabled by enhanced chromatin accessibility. These dual roles of canonical NF-κB in Tconv and Treg cells highlight the functional plasticity of the NF-κB signaling pathway and underscores the need for more selective strategies to therapeutically target NF-κB.

To ensure the completion of DNA replication and maintenance of genome integrity, DNA repair factors protect stalled replication forks upon replication stress. Previous studies have identified a critical role for the tumor suppressors BRCA1 and BRCA2 in preventing the degradation of nascent DNA by the MRE11 nuclease after replication stress. Here we show that depletion of SMARCAL1, a SNF2-family DNA translocase that remodels stalled forks, restores replication fork stability and reduces the formation of replication stress-induced DNA breaks and chromosomal aberrations in BRCA1/2-deficient cells. In addition to SMARCAL1, other SNF2-family fork remodelers, including ZRANB3 and HLTF, cause nascent DNA degradation and genomic instability in BRCA1/2-deficient cells upon replication stress. Our observations indicate that nascent DNA degradation in BRCA1/2-deficient cells occurs as a consequence of MRE11-dependent nucleolytic processing of reversed forks generated by fork remodelers. These studies provide mechanistic insights into the processes that cause genome instability in BRCA1/2- deficient cells.

The vast majority of the mammalian genome has the potential to expressnoncoding RNA (ncRNA). The 11-subunit RNA exosome complex is the main source of cellular 3′–5′ exoribonucleolytic activity and potentially regulates the mammalian noncoding transcriptome1. Here we generated a mouse model in which the essential subunit Exosc3 of the RNA exosome complex can be conditionally deleted. Exosc3-deficient B cells lack the ability to undergo normal levels of class switch recombination and somatic hypermutation, two mutagenic DNA processes used to generate antibody diversity via the B-cell mutator protein activation-induced cytidine deaminase (AID)2,3. The transcriptome of Exosc3-deficient B cells has revealed the presence of many novel RNA exosome substrate ncRNAs. RNA exosome substrate RNAs include xTSS-RNAs, transcription start site (TSS)-associated antisense transcripts that can exceed 500 base pairs in length and are transcribed divergently from cognate coding gene transcripts. xTSS-RNAs are most strongly expressed at genes that accumulate AID-mediated somatic mutations and/or are frequent translocation partners of DNA double-strand breaks generated at Igh in B cells4,5. Strikingly, translocations near TSSs or within gene bodies occur over regions of RNA exosome substrate ncRNA expression. These RNA exosome-regulated, antisense-transcribed regions of the B-cell genome recruit AID and accumulate single-strand DNA structures containing RNA–DNA hybrids. We propose that RNA exosome regulation of ncRNA recruits AID to single-strand DNA-forming sites of antisense and divergent transcription in the B-cell genome, thereby creating a link between ncRNA transcription and overall maintenance of B-cell genomic integrity.

SummaryFollicular lymphoma (FL) is an indolent disease, but 30%–40% of cases undergo histologic transformation to an aggressive malignancy, typically represented by diffuse large B cell lymphoma (DLBCL). The pathogenesis of this process remains largely unknown. Using whole-exome sequencing and copy-number analysis, we show here that the dominant clone of FL and transformed FL (tFL) arise by divergent evolution from a common mutated precursor through the acquisition of distinct genetic events. Mutations in epigenetic modifiers and antiapoptotic genes are introduced early in the common precursor, whereas tFL is specifically associated with alterations deregulating cell-cycle progression and DNA damage responses (CDKN2A/B, MYC, and TP53) as well as aberrant somatic hypermutation. The genomic profile of tFL shares similarities with that of germinal center B cell-type de novo DLBCL but also displays unique combinations of altered genes with diagnostic and therapeutic implications.

Amidst the growing literature on cancer genomics and intratumor heterogeneity, essential principles in evolutionary biology recur time and time again. Here we use these principles to guide the reader through major advances in cancer research, highlighting issues of “hit hard, hit early” treatment strategies, drug resistance, and metastasis. We distinguish between two frameworks for understanding heterogeneous tumors, both of which can inform treatment strategies: (1) The tumor as diverse ecosystem , a Darwinian population of sometimes-competing, sometimes-cooperating cells; (2) The tumor as tightly integrated, self-regulating organ , which may hijack developmental signals to restore functional heterogeneity after treatment. While the first framework dominates literature on cancer evolution, the second framework enjoys support as well. Throughout this review, we illustrate how mathematical models inform understanding of tumor progression and treatment outcomes. Connecting models to genomic data faces computational and technical hurdles, but high-throughput single-cell technologies show promise to clear these hurdles. This article is part of a Special Issue entitled: Evolutionary principles - heterogeneity in cancer?, edited by Dr. Robert A. Gatenby.

Patients with myelodysplasia or acute myeloid leukemia (AML) show decreased osteoblast numbers, and this is reproduced in mouse models of acute leukemia. Manipulating the osteoblast pool directly affects leukemia progression: osteoblast ablation increases leukemia burden, while maintenance of osteoblast numbers decreases it and prolongs survival. Here we show that the protective effect of osteoblasts against leukemia does not rely on a certain threshold of osteoblast numbers but rather requires a specific signaling pathway. Treatment of leukemic mice with parathyroid hormone, which increases osteoblast numbers, had no effect in leukemia progression. In contrast, maintaining the osteoblast pool in leukemic mice by decreasing the synthesis of gut-derived serotonin, reduced AML burden and prolonged survival. Using the human-AML MLL-AF9 model, we examined whether deletion of one of the main serotonin receptors expressed in osteoblasts affects leukemia progression. Global deletion of Htr1b in mice prevented leukemia. While removal of Htr1b in LepR+ positive mesenchymal stem cells had no effect in leukemia progression, its deletion in Col1a1-expressing osteoblast precursors significantly prevented lethality in the mutant mice. In order to dissect the crosstalk between leukemia-cells and osteoblasts, we co-cultured primary-human osteoblasts and AML cells. RNAseq analysis identified a signature of chemoattractants and pro-inflammatory cytokines secreted by osteoblasts upon exposure to AML cells, suggesting that the protective mechanism exerted by osteoblasts activates the anti-tumor immune response against leukemia. In summary, our data suggest that AML cells activate serotonin signaling in osteoblasts to promote their engraftment and expansion. Targeting this pathway may render the niche hostile to leukemia by enhancing the anti-tumor immune response elicited by osteoblasts.

Somatic mutations in hematopoietic precursors underlie the development of myeloid disorders, such as myeloproliferative neoplasms (MPN). However, our ability to interrogate the transcriptional impact of these mutations on human hematopoiesis is limited by the frequent admixing of mutant (MUT) with wildtype (WT) cells or with other subclones. Recently, digital single-cell RNA-sequencing has provided high-resolution maps of normal hematopoiesis. Nonetheless, due to their 3' bias, these methods do not capture the cell's mutational status. Efficient linking of single-cell genotype and transcriptomes would allow direct comparison of WT and MUT progenitors within the same sample, eliminating patient-specific and technical confounders. Thus, we developed single-cell Genotyping of Transcriptomes (GoT) to link genotypes of expressed genes to transcriptional profiling of thousands of cells by adapting the 10x Genomics platform. We capture the target locus from cDNA generated at an intermediate step, thus enabling linkage of genotype to whole transcriptomes via shared barcodes (Fig. 1A). We tested this approach via a species-mixing experiment, whereby mouse cells with MUT CALR were mixed with human cells with WT CALR. GoT of 1291 admixed cells provided genotyping for 97.5% of cells, and 96.7% matched the expected species (Fig. 1B). To demonstrate the ability of this technology to probe hematopoietic differentiation in MPN, we applied GoT to 20,908 CD34+ cells across five patients with CALR-mutated essential thrombocythemia (ET) or myelofibrosis (MF), resulting in genotyping of 82% of cells. We first performed clustering agnostic to genotype, based on transcriptome data alone, and found that cells clustered according to progenitor cell identity, rather than mutational status (Fig. 1C). Furthermore, projection of GoT data demonstrated that MUT cells were present across all progenitor clusters (Fig. 1D). However, the frequency of CALR-mutated cells was higher in committed progenitors, especially megakaryocytic progenitors (MkPs) compared to CD34+, CD38- hematopoietic stem progenitor cells (HSCPs, Fig. 1E). Thus, CALR mutation may confer a greater fitness impact in lineage-committed cells vs. HSPCs. Indeed, we found a significant increase in the number of MkPs in cell cycle in MUT cells compared to WT cells (Fig. 1F). Moreover, this increase in cell cycle activity correlated with the platelet count (Fig. 1G). This suggests that interrogation of MUT and WT progenitors may inform our understanding of patient phenotypic variability despite shared genotypes. GoT enables de novo differential expression discovery in MUT vs. WT cells within the same progenitor subset. MUT MkPs upregulated genes in the unfolded protein response, such as PDIA6, HSPA5 and XBP1 (Fig. 1H), consistent with the central role of CALR as a chaperone protein. On the other hand, MUT HSPCs showed upregulation of the NF-kB pathway (Fig. 1I), most significantly in the subcluster enriched with the earliest HSCs (Fig. 1J). Since the NF-kB pathway has been implicated in HSC self-renewal, our data provides a potential mechanism for clonal expansion and maintenance of CALR-mutated HSCs. Collectively, these findings demonstrate that the transcriptomic output of CALR mutations is closely dependent on cell identity. To further evaluate the potential of GoT to detect multiple genotypes in clonally complex neoplasms, we targeted three genes, clonal SF3B1 (VAF 47.5% by bulk exon sequencing), and subclonal CALR (43.5%) and NFE2 (33%), in CD34+ cells from a patient with MF (Fig. 1K). Through GoT, the subclonal transcriptional output was interrogated; for example, CALR mutation conferred proliferative advantage to megakaryocytic-erythroid progenitors even in the presence of SF3B1 mutation, while additional NFE2 mutation did not further increase cell cycle activity. In summary, GoT is a powerful tool for linking transcriptional changes to somatic genotypes at the single-cell level. Specifically, it uncovered the transcriptional impact of mutations in myeloid clonal growths in the context of distinct progenitor identities. Further application of GoT to additional MPN contexts as well as clonal hematopoiesis is thus anticipated to provide critical insights into the transcriptional programs that enable clonal expansion and evolution in human hematopoiesis. Disclosures Hoffman: Merus: Research Funding; Janssen: Research Funding; Summer Road: Research Funding; Incyte: Research Funding; Formation Biologics: Research Funding.

Programmed mutagenesis of the immunoglobulin locus of B lymphocytes during class switch recombination (CSR) and somatic hypermutation requires RNA polymerase II (polII) transcription complex-dependent targeting of the DNA mutator activation-induced cytidine deaminase (AID). AID deaminates cytidine residues on substrate sequences in the immunoglobulin (Ig) locus via a transcription-dependent mechanism, and this activity is stimulated by the RNA polII stalling cofactor Spt5 and the 11-subunit cellular noncoding RNA 3′–5′ exonucleolytic processing complex RNA exosome. The mechanism by which the RNA exosome recognizes immunoglobulin locus RNA substrates to stimulate AID DNA deamination activity on its in vivo substrate sequences is an important question. Here we report that E3-ubiquitin ligase Nedd4 destabilizes AID-associated RNA polII by a ubiquitination event, leading to generation of 3′ end free RNA exosome RNA substrates at the Ig locus and other AID target sequences genome-wide. We found that lack of Nedd4 activity in B cells leads to accumulation of RNA exosome substrates at AID target genes and defective CSR. Taken together, our study links noncoding RNA processing following RNA polII pausing with regulation of the mutator AID protein. Our study also identifies Nedd4 as a regulator of noncoding RNAs that are generated by stalled RNA polII genome-wide.

Glioblastoma is one of the most challenging forms of cancer to treat. Here we describe a computational platform that integrates the analysis of copy number variations and somatic mutations and unravels the landscape of in-frame gene fusions in glioblastoma. We found mutations with loss of heterozygosity in LZTR1, encoding an adaptor of CUL3-containing E3 ligase complexes. Mutations and deletions disrupt LZTR1 function, which restrains the self renewal and growth of glioma spheres that retain stem cell features. Loss-of-function mutations in CTNND2 target a neural-specific gene and are associated with the transformation of glioma cells along the very aggressive mesenchymal phenotype. We also report recurrent translocations that fuse the coding sequence of EGFR to several partners, with EGFR-SEPT14 being the most frequent functional gene fusion in human glioblastoma. EGFR-SEPT14 fusions activate STAT3 signaling and confer mitogen independence and sensitivity to EGFR inhibition. These results provide insights into the pathogenesis of glioblastoma and highlight new targets for therapeutic intervention. © 2013 Nature America, Inc. All rights reserved.

Background The extraordinary success of imatinib in the treatment of BCR-ABL1 associated cancers underscores the need to identify novel functional gene fusions in cancer. RNA sequencing offers a genome-wide view of expressed transcripts, uncovering biologically functional gene fusions. Although several bioinformatics tools are already available for the detection of putative fusion transcripts, candidate event lists are plagued with non-functional read-through events, reverse transcriptase template switching events, incorrect mapping, and other systematic errors. Such lists lack any indication of oncogenic relevance, and they are too large for exhaustive experimental validation. Results We have designed and implemented a pipeline, Pegasus, for the annotation and prediction of biologically functional gene fusion candidates. Pegasus provides a common interface for various gene fusion detection tools, reconstruction of novel fusion proteins, reading-frame-aware annotation of preserved/lost functional domains, and data-driven classification of oncogenic potential. Pegasus dramatically streamlines the search for oncogenic gene fusions, bridging the gap between raw RNA-Seq data and a final, tractable list of candidates for experimental validation. Conclusion We show the effectiveness of Pegasus in predicting new driver fusions in 176 RNA-Seq samples of glioblastoma multiforme (GBM) and 23 cases of anaplastic large cell lymphoma (ALCL). Contact: fa2306@columbia.edu.

Missense mutations in ATM kinase, a master regulator of DNA damage responses, are found in many cancers, but their impact on ATM function and implications for cancer therapy are largely unknown. Here we report that 72% of cancer-associated ATM mutations are missense mutations that are enriched around the kinase domain. Expression of kinase-dead ATM (AtmKD/-) is more oncogenic than loss of ATM (Atm-/-) in mouse models, leading to earlier and more frequent lymphomas with Pten deletions. Kinase-dead ATM protein (Atm-KD), but not loss of ATM (Atm-null), prevents replication-dependent removal of Topo-isomerase I-DNA adducts at the step of strand cleavage, leading to severe genomic instability and hypersensitivity to Topo-isomerase I inhibitors. Correspondingly, Topo-isomerase I inhibitors effectively and preferentially eliminate AtmKD/-, but not Atm-proficientor Atm-/- leukemia in animal models. These findings identify ATM kinase-domain missense mutations as a potent oncogenic event and a biomarker for Topo-isomerase I inhibitor based therapy. DOI: http://dx.doi.org/10.7554/eLife.14709.001, eLife digest Cancer is a genetic disease. To remain healthy, therefore, it is essential that cells do not accrue too many dangerous mutations in their DNA that allow cancers to grow and develop. An enzyme called ATM helps to do just that. DNA damage activates ATM, which, in turn, adds phosphate groups to other proteins. These newly tagged proteins then stop cells dividing until the DNA has been repaired. Human cancers often switch off ATM, either by completely deleting the enzyme or mutating it. This renders ATM unable to add phosphate groups to proteins, and so allows the cancer cells to continue proliferating even in the face of DNA damage. Yamamoto et al. wanted to know whether cancers that completely lack ATM behave differently from cancers that contain an inactive version of the enzyme. Studying mice that were engineered to have an inactive version of ATM in their blood cells showed that such mice developed blood cancers faster than mice have no ATM in their blood cells. In particular, cancer cells with the inactive form of the ATM enzyme accumulated more DNA damage than cells that lacked the enzyme completely. Using biochemical techniques, Yamamoto et al. then showed that the inactive form of ATM can prevent other enzymes from repairing DNA. Drugs that inhibit one of these repair enzymes – called topo-isomerase I – are already used in cancer treatments. These drugs were particularly effective on tumors with the inactive version of the ATM enzyme. As ATM is commonly mutated in human cancers, the next steps that follow on from this research are to develop methods to test which cancers contain the inactive form of the ATM enzyme. Clinical trials could then investigate how effectively topo-isomerase I inhibitors treat these specific types of cancer. DOI: http://dx.doi.org/10.7554/eLife.14709.002

Glioblastoma (GBM) is the most common and aggressive primary brain tumor. To better understand how GBM evolves, we analyzed longitudinal genomic and transcriptomic data from 114 patients. The analysis shows a highly branched evolutionary pattern in which 63% of patients experience expression-based subtype changes. The branching pattern, together with estimates of evolutionary rate, suggests that relapse-associated clones typically existed years before diagnosis. Fifteen percent of tumors present hypermutation at relapse in highly expressed genes, with a clear mutational signature. We find that 11% of recurrence tumors harbor mutations in LTBP4, which encodes a protein binding to TGF-β. Silencing LTBP4 in GBM cells leads to suppression of TGF-β activity and decreased cell proliferation. In recurrent GBM with wild-type IDH1, high LTBP4 expression is associated with worse prognosis, highlighting the TGF-β pathway as a potential therapeutic target in GBM.

Nodal marginal zone lymphoma (NMZL) is a rare, indolent B-cell tumor that is distinguished from splenic marginal zone lymphoma (SMZL) by the different pattern of dissemination. NMZL still lacks distinct markers and remains orphan of specific cancer gene lesions. By combining whole-exome sequencing, targeted sequencing of tumor-related genes, whole-transcriptome sequencing, and high-resolution single nucleotide polymorphism array analysis, we aimed at disclosing the pathways that are molecularly deregulated in NMZL and we compare the molecular profile of NMZL with that of SMZL. These analyses identified a distinctive pattern of nonsilent somatic lesions in NMZL. In 35 NMZL patients, 41 genes were found recurrently affected in ≥3 (9%) cases, including highly prevalent molecular lesions of MLL2 (also known as KMT2D; 34%), PTPRD (20%), NOTCH2 (20%), and KLF2 (17%). Mutations of PTPRD, a receptor-type protein tyrosine phosphatase regulating cell growth, were enriched in NMZL across mature B-cell tumors, functionally caused the loss of the phosphatase activity of PTPRD, and were associated with cell-cycle transcriptional program deregulation and increased proliferation index in NMZL. Although NMZL shared with SMZL a common mutation profile, NMZL harbored PTPRD lesions that were otherwise absent in SMZL. Collectively, these findings provide new insights into the genetics of NMZL, identify PTPRD lesions as a novel marker for this lymphoma across mature B-cell tumors, and support the distinction of NMZL as an independent clinicopathologic entity within the current lymphoma classification.

// Monica Messina 1,* , Sabina Chiaretti 1,* , Jiguang Wang 2,* , Anna Lucia Fedullo 1 , Nadia Peragine 1 , Valentina Gianfelici 1 , Alfonso Piciocchi 3 , Fulvia Brugnoletti 1 , Filomena Di Giacomo 4 , Simona Pauselli 1 , Antony B. Holmes 5 , Maria Cristina Puzzolo 1 , Giulia Ceglie 1 , Valerio Apicella 1 , Marco Mancini 1 , Geertruy te Kronnie 6 , Anna Maria Testi 1 , Antonella Vitale 1 , Marco Vignetti 3 , Anna Guarini 1 , Raul Rabadan 2 and Robin Foa 1 1 Hematology, Department of Cellular Biotechnologies and Hematology, Sapienza University, Rome, Italy 2 Department of Systems Biology, Biomedical Informatics and Center for Computational Biology and Bioinformatics, Columbia University, New York, NY, USA 3 GIMEMA Data Center, Rome, Italy 4 Department of Molecular Biotechnology and Health Science, and Center for Experimental Research and Medical Studies (CeRMS), University of Torino, Torino, Italy 5 Institute for Cancer Genetics and The Herbert Irving Comprehensive Cancer Center, Columbia University, New York, NY, USA 6 Department of Women’s and Children’s Health, University of Padova, Padova, Italy * These authors have contributed equally to this work Correspondence to: Robin Foa, email: // Keywords : acute lymphoblastic leukemia, next generation sequencing, copy number aberrations, novel prognostic markers, genetic-driven targeted therapy Received : October 23, 2015 Accepted : January 28, 2016 Published : February 12, 2016 Abstract To shed light into the molecular bases of B-lineage acute lymphoblastic leukemia lacking known fusion transcripts, i.e. BCR-ABL1 , ETV6-RUNX1 , E2A-PBX1 , and MLL rearrangements (B-NEG ALL) and the differences between children, adolescents/young adults (AYA) and adults, we analyzed 168 B-NEG ALLs by genome-wide technologies. This approach showed that B-NEG cases carry 10.5 mutations and 9.1 copy-number aberrations/sample. The most frequently mutated druggable pathways were those pertaining to RAS/RTK (26.8%) and JAK/STAT (12.5%) signaling. In particular, FLT3 and JAK/STAT mutations were detected mainly in AYA and adults, while KRAS and NRAS mutations were more frequent in children. RAS/RTK mutations negatively affected the outcome of AYA and adults, but not that of children. Furthermore, adult B-NEG ALL carrying JAK/STAT mutations had a shorter survival. In vitro experiments showed that FLT3 inhibitors reduced significantly the proliferation of FLT3 -mutated primary B-NEG ALL cells. Likewise, PI3K/mTOR inhibitors reduced the proliferation of primary cells harboring RAS and IL7R mutations. These results refine the genetic landscape of B-NEG ALL and suggest that the different distribution of lesions and their prognostic impact might sustain the diverse outcome between children, adults and partly AYA - whose genomic scenario is similar to adults - and open the way to targeted therapeutic strategies.

Background The Egyptian Rousette bat (Rousettus aegyptiacus), a common fruit bat species found throughout Africa and the Middle East, was recently identified as a natural reservoir host of Marburg virus. With Ebola virus, Marburg virus is a member of the family Filoviridae that causes severe hemorrhagic fever disease in humans and nonhuman primates, but results in little to no pathological consequences in bats. Understanding host-pathogen interactions within reservoir host species and how it differs from hosts that experience severe disease is an important aspect of evaluating viral pathogenesis and developing novel therapeutics and methods of prevention. Results Progress in studying bat reservoir host responses to virus infection is hampered by the lack of host-specific reagents required for immunological studies. In order to establish a basis for the design of reagents, we sequenced, assembled, and annotated the R. aegyptiacus transcriptome. We performed de novo transcriptome assembly using deep RNA sequencing data from 11 distinct tissues from one male and one female bat. We observed high similarity between this transcriptome and those available from other bat species. Gene expression analysis demonstrated clustering of expression profiles by tissue, where we also identified enrichment of tissue-specific gene ontology terms. In addition, we identified and experimentally validated the expression of novel coding transcripts that may be specific to this species. Conclusion We comprehensively characterized the R. aegyptiacus transcriptome de novo. This transcriptome will be an important resource for understanding bat immunology, physiology, disease pathogenesis, and virus transmission. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-2124-x) contains supplementary material, which is available to authorized users.

Although histopathologic classification of adult gliomas is well established, accurate diagnosis suffers from high intra- and inter-observer variability. In order to investigate the molecular underpinnings we performed an unbiased pan-glioma analysis of molecular features. We assembled the largest possible primary glioma dataset to date consisting of 1,122 TCGA samples of grades II-IV glioma that have been analyzed through mRNA sequencing and expression arrays (n = 1045), DNA methylation arrays (n = 932), exome sequencing (n = 833), DNA copy number arrays (n = 1084) and whole genome sequencing (WGS, n = 141). In addition to 30 known glioma drivers like IDH1, TP53, and EGFR, we identified 70 novel significantly mutated genes which were mutated at 0.5-3% frequency, including KRAS, NOTCH2 and ARID2. Novel copy number drivers included SETD2, ARID2, and GIGYF2. Analysis of WGS data identified enrichment for mutations in the promoter regions of TERT, CACNG6 and TRIM28 that were associated with a significant difference in mRNA expression, suggesting functional consequences. We used WGS based telomere length estimates and found that ATRX mutated samples showed significantly longer telomeres while TERT mutant samples showed shortened telomeres, independent of differences in age and IDH status. Unsupervised clustering of DNA methylation and gene expression resulted in six methylation subtypes and four expression subtypes. These clusters overlapped significantly disregarding tumor histopathology and could be divided into three macro groups separated by IDH mutations, TERT promoter mutations and 1p/19q co-deletion. Functional enrichments within-cluster revealed that grade IV tumors are primarily marked by a hyperproliferation signature with respect to grade II-III tumors which are characterized by increased activation of genes related to glial-neuronal functioning. These findings suggest that tumor grade reflects a point in time in the temporal evolution within individual glioma subtypes. Our analyses identified new drivers; suggest a novel association between telomere length and maintenance mechanisms and shed light on glioma evolution.

Endemic Burkitt lymphoma (eBL) is primarily found in children in equatorial regions and represents the first historical example of a virus-associated human malignancy. Although Epstein-Barr virus (EBV) infection and MYC translocations are hallmarks of the disease, it is unclear whether other factors may contribute to its development. We performed RNA-Seq on 20 eBL cases from Uganda and showed that the mutational and viral landscape of eBL is more complex than previously reported. First, we found the presence of other herpesviridae family members in 8 cases (40%), in particular human herpesvirus 5 and human herpesvirus 8 and confirmed their presence by immunohistochemistry in the adjacent non-neoplastic tissue. Second, we identified a distinct latency program in EBV involving lytic genes in association with TCF3 activity. Third, by comparing the eBL mutational landscape with published data on sporadic Burkitt lymphoma (sBL), we detected lower frequencies of mutations in MYC, ID3, TCF3 and TP53, and a higher frequency of mutation in ARID1A in eBL samples. Recurrent mutations in two genes not previously associated with eBL were identified in 20% of tumors: RHOA and cyclin F (CCNF). We also observed that polyviral samples showed lower numbers of somatic mutations in common altered genes in comparison to sBL specimens, suggesting dual mechanisms of transformation, mutation versus virus driven in sBL and eBL respectively., Author Summary Burkitt lymphoma is endemic in sub-Saharan Africa and affects primarily children of age 4–7 years. Historically, it was one of the first tumors associated with a virus (EBV) and bearing a translocation involving an oncogene, i.e. MYC. There are three distinct clinical variants of Burkitt lymphoma according to the World Health Organization: sporadic, endemic and immunodeficiency-related. Although there has been some recent work on the molecular characterization of sporadic Burkitt lymphomas, little is known about the pathogenesis of endemic cases. In this work, we analyzed 20 samples of RNASeq from Burkitt lymphoma collected in Lacor Hospital (Uganda, Africa) and validated in an extension panel of 73 samples from Uganda and Kenya. We identify the presence in the adjacent non-neoplastic tissue of other herpesviridae family members in 53% of the cases, namely cytomegalovirus (CMV) and Kaposi sarcoma herpesvirus (KSHV). We also demonstrate expression of EBV lytic genes in primary tumor samples and find an inverse association between EBV lytic expression and TCF3 activity. When studying the mutational profile of endemic Burkitt tumors, we find recurrent alterations in genes rarely mutated in sporadic Burkitt lymphomas, i.e. ARID1A, CCNF and RHOA, and lower numbers of mutations in genes previously reported to be commonly mutated in sporadic cases, i.e. MYC, ID3, TCF3, TP53. Together, these results illustrate a distinct genetic and viral profile of endemic Burkitt lymphoma, suggesting a dual mechanism of transformation (mutation versus virus driven in sBL and eBL respectively).

Next generation sequencing and copy number analysis provide insights into the complexity of the CLL coding genome, and reveal an association between NOTCH1 mutational activation and poor prognosis., The pathogenesis of chronic lymphocytic leukemia (CLL), the most common leukemia in adults, is still largely unknown. The full spectrum of genetic lesions that are present in the CLL genome, and therefore the number and identity of dysregulated cellular pathways, have not been identified. By combining next-generation sequencing and copy number analysis, we show here that the typical CLL coding genome contains

Human populations contain a functional coding polymorphism (codon 72) in the p53 gene. To explore whether this polymorphism alters the transcriptional pattern of p53-regulated genes, the human isogenic cell lines harboring p53 with either the proline or arginine at codon 72 were employed to activate p53-mediated transcription. Thirty-four p53-regulated genes were assayed for their increased levels of mRNA using quantitative real-time PCR. The largest difference between p53-arginine and p53-proline was found with the PERP gene involved in cell-cell adhesion and apoptosis. The most common set of genes that are transcribed better by the p53-arginine protein than the p53-proline protein was found in the apoptotic function (DR-4, NOXA, PUMA, and PIG-3). LIF, a cytokine that is required for optimal reproductive function, was produced at 2x higher levels by the p53-arginine than the p53-proline allele. The genes that induced their mRNAs at the highest levels compared to the baseline tended to be synthesized better by the p53-arginine protein than the p53-proline protein. These molecular studies may help to explain the complicated associations observed between this polymorphism and the incidence of some cancers, the longevity of some populations, and the fecundity of different groups.

Segments of the 2009 human H1N1 strains have coexisted in swine influenza virus strains for more than 10 years. Vladimir Trifonov, Hossein Khiabanian, and Raul Rabadan describe the evolution of the H1N1 virus.