23014 results for "ega"

in 43.82 milliseconds.

• Organoid BulkRNAseq

  Here, we aim to understand early determinants of brain organoid quality and generated a set of brain organoids from 12 different hPSC lines using an adaptation of the unguided Lancaster protocol (Lancaster & Knoblich, 2014). All cell lines were obtained from healthy donors. 12 donor hPSC lines were differentiated into brain organoids and a set of 72 individual organoids (6 per line) were randomly selected morphologically diverse organoids and subjected to bulk RNA sequencing for gene expression analysis. RNA was extracted using the RNeasy kit (Qiagen, Germany) according to the manufacturer’s instructions. A total of 1000 ng per sample was sent for RNA sequencing to Azenta Life Sciences (Genewiz Leipzig, Germany) for sequencing. Library preparation was performed using the NEBNext Ultra II RNA Library Prep Kit for Illumina following manufacturer’s instructions (NEB, Ipswich, MA, USA). mRNAs were first enriched with Oligo(dT) beads. The samples were sequenced using a 2x150 Pair-End (PE) configuration v1.5 at a depth of >20 million reads in each sample. Raw sequence data (.bcl files) generated from Illumina NovaSeq6000 was converted into fastq files and de-multiplexed using Illumina bcl2fastq program version 2.20. The raw sequence reads were aligned to the human reference genome GRCh38 using the STAR aligner (Dobin et al., 2013). Gene expression levels were quantified using the FeatureCounts tool (Liao et al., 2014).

  Study EGAS50000000659

• Nala TAS-LRS PGx Study

  Pharmacogenomics (PGx) testing improves medication safety and efficacy by identifying genetic variants that affect drug response. However, current technologies often fail to resolve complex loci, detect structural variants, or phase alleles accurately. Here, we present an end-to-end PGx workflow based on Targeted Adaptive Sampling-Long Read Sequencing (TAS-LRS), integrating a streamlined laboratory protocol with a bioinformatics pipeline that includes a novel CYP2D6 caller. Using 1,000 ng of DNA and three-sample multiplexing on a single PromethION flow cell, the assay achieves consistent on-target (25x) and off-target (3x) coverage, enabling accurate, haplotype-resolved testing of 35 pharmacogenes alongside genome-wide genotyping from off-target reads. We further developed the workflow into a clinically ready service and validated its performance across 17 reference and clinical samples. The assay demonstrated high concordance for small variants (99.9%) and structural variants (>95%), with phased diplotypes and metabolizer phenotypes reaching 97.7% and 98.0% concordance, respectively. Improved calls were observed in 12 genes due to enhanced genotyping, phasing, or novel allele detection. In addition, off-target reads supported accurate genome-wide imputation, comparable to short-read sequencing and superior to microarrays. These results establish the feasibility of long-read sequencing for clinical PGx testing and position TAS-LRS as a scalable solution combining both targeted and genome-wide utility.

  Study EGAS50000001116

• Somatic rearrangements causing oncogenic ectodomain deletions of FGFR1 in squamous cell lung cancer

  The discovery of frequent 8p11-p12 amplifications in squamous cell lung cancer (SQLC) has fueled hopes that FGFR1, located inside this amplicon, might be a therapeutic target. In a clinical trial, only 11% of patients with 8p11 amplification (detected by FISH) responded to FGFR kinase inhibitor treatment. To understand the mechanism of FGFR1 dependency, we performed deep genomic characterization of 52 SQLCs with 8p11-p12 amplification, including 10 tumors obtained from patients who had been treated with FGFR inhibitors. We discovered somatically altered variants of FGFR1 with deletion of exons 1-8 that resulted from intragenic tail-to-tail rearrangements. These ectodomain-deficient FGFR1 variants (ΔEC-FGFR1) were expressed in the affected tumors and were tumorigenic in both in vitro and in vivo models of lung cancer. Mechanistically, breakage-fusion-bridges were the source of 8p11-p12 amplification, resulting from frequent head-to-head and tail-to-tail rearrangements. Generally, tail-to-tail rearrangements within or in close proximity upstream of FGFR1 were associated with FGFR1 dependency. Thus, the genomic events shaping the architecture of the 8p11-p12 amplicon provide a mechanistic explanation for the emergence of FGFR1-driven SQLC. Specifically, we believe that FGFR1 ectodomain-deficient and FGFR1-centered amplifications caused by tail-to-tail rearrangements are a novel somatic genomic event that might be predictive of therapeutically relevant FGFR1 dependency.

  Study EGAS50000001368

• Transposable Elements Shape Stemness in Normal and Leukemic Hematopoiesis

  Despite most acute myeloid leukemia (AML) patients achieving complete remission after induction chemotherapy, two-thirds will relapse with fatal disease within five years. AML is organized as a cellular hierarchy sustained by leukemia stem cells (LSC) at the apex, with LSC properties directly linked to tumor progression, therapy failure, and disease relapse 1–5. Despite the central role of LSC in poor patient outcomes, little is known about the genetic determinants driving their stemness properties. As LSCs share many functional and molecular properties with normal hematopoietic stem cells (HSC) 6, we investigated accessible chromatin unique across normal hematopoietic and cancer cell states and identified transposable elements (TEs) as genetic determinants of both primitive populations in comparison with their downstream mature progeny. A clinically-relevant TE chromatin accessibility-based LSCTE121 signature was developed that enabled patient classification based on survival outcomes. Through functional assays, primitive cell specific-TE subfamilies were found to serve as docking sites for stem cell-associated regulators of genome topology or lineage-specific transcription factors, including LYL1 in LSCs. Finally, using chromatin editing tools, we establish that chromatin accessibility at LTR12C elements in LSCs are necessary to maintain stemness properties. Our work identifies TEs as genetic drivers of primitive versus mature cell states, where distinct TE subfamilies account for stemness properties in normal versus leukemic hematopoietic stem cells.

  Study EGAS50000001141

• Transposable Elements Shape Stemness in Normal and Leukemic Hematopoiesis

  Despite most acute myeloid leukemia (AML) patients achieving complete remission after induction chemotherapy, two-thirds will relapse with fatal disease within five years. AML is organized as a cellular hierarchy sustained by leukemia stem cells (LSC) at the apex, with LSC properties directly linked to tumor progression, therapy failure, and disease relapse 1–5. Despite the central role of LSC in poor patient outcomes, little is known about the genetic determinants driving their stemness properties. As LSCs share many functional and molecular properties with normal hematopoietic stem cells (HSC) 6, we investigated accessible chromatin unique across normal hematopoietic and cancer cell states and identified transposable elements (TEs) as genetic determinants of both primitive populations in comparison with their downstream mature progeny. A clinically-relevant TE chromatin accessibility-based LSCTE121 signature was developed that enabled patient classification based on survival outcomes. Through functional assays, primitive cell specific-TE subfamilies were found to serve as docking sites for stem cell-associated regulators of genome topology or lineage-specific transcription factors, including LYL1 in LSCs. Finally, using chromatin editing tools, we establish that chromatin accessibility at LTR12C elements in LSCs are necessary to maintain stemness properties. Our work identifies TEs as genetic drivers of primitive versus mature cell states, where distinct TE subfamilies account for stemness properties in normal versus leukemic hematopoietic stem cells.

  Study EGAS50000001147

• Neoadjuvant immune checkpoint blockade in mismatch-repair proficient colon cancers

  Immune checkpoint blockade has led to paradigm shifts in the treatment of various tumour types, yet limited efficacy has been observed in patients with metastatic mismatch-repair proficient (pMMR) colorectal cancer. Here we report an in-depth analysis of patients with early-stage pMMR colon cancer from the phase II NICHE study (ClinicalTrials.gov: NCT03026140). A total of 31 patients received neoadjuvant treatment of nivolumab plus ipilimumab followed by surgery. The response rate was 26% and included six patients with a major pathological response (≤10% residual viable tumour). One patient with an ongoing clinical complete response did not undergo surgery. Responses were observed despite a low tumour mutational burden in all tumours, while chromosomal genomic instability scores were significantly higher in responders compared to non-responders. Responding tumours were more likely to be TP53 mutant, KRASG12 wildtype and had significantly higher baseline expression of proliferation signatures and TCF1. Imaging mass cytometry also revealed a higher percentage of Ki-67+ cancer and Ki-67+ CD8+ T cells in responders compared to non-responders. Meanwhile, in non-responders, we observed an enrichment for an inflamed-fibrotic microenvironment associated with TGF-β expression, hinting towards possible resistance mechanisms. These results provide a comprehensive analysis of response to neoadjuvant ICB in early-stage pMMR colon cancers and identify potential biomarkers for patient selection.

  Study EGAS50000000856

• Immune heterogeneity in small cell lung cancer and vulnerability to immune checkpoint blockade

  Following the pivotal Phase III IMpower133 study, atezolizumab (anti–PD-L1), combined with carboplatin and etoposide (CE), was the first immune checkpoint inhibitor approved for first-line treatment of extensive-stage small cell lung cancer (ES-SCLC) and is now a standard of care. A clearer understanding of therapeutically relevant SCLC subsets is needed to improve outcomes and to identify rational combination strategies. Transcriptomic analyses and nonnegative matrix factorization were conducted on 271 pre-treatment patient tumor samples from IMpower133 and four subsets with general concordance to previously reported SCLC subtypes were identified. Deeper investigation into the immune heterogeneity relevant to clinical outcomes uncover two subsets with differing neuroendocrine (NE) versus non-neuroendocrine (non-NE) phenotypes that demonstrated hallmarks of immune cell infiltration. The balance of tumor-associated macrophage (TAM) to T-effector signals distinguished these inflamed subsets. Tumors with low TAM but high T-effector signals had a NE phenotype and demonstrated longer overall survival with PD-L1 blockade and CE versus CE alone than did tumors with high TAM and high T-effector signal, which had a non-NE phenotype. The delineation of SCLC immune heterogeneity offers a clinically relevant approach to discriminate SCLC patients likely to benefit most from immunotherapeutic approaches and highlights the complex mechanisms underlying response and resistance to immune checkpoint blockade. 

  Study EGAS50000000138

• Phenotypic characterization and prognostic impact of CD103+ tissue-resident memory T cells in diffuse large B cell lymphoma

  Tissue-resident memory T (TRM) cells are a population of memory T cells that stably occupy tissues and play a key role in immunosurveillance, having been linked to favorable survival outcomes in various solid tumors. While TRM cells have been identified in lymph nodes, their phenotype and prognostic significance in B-cell non-Hodgkin lymphoma (B-NHL), including diffuse large B-cell lymphoma (DLBCL), remains poorly characterized. The frequency of CD103 expressing T cells in patient samples with DLCBL was quantified by immunofluorescence (IF) staining of tissue biopsies (n=306) and flow cytometry (n=252) of cell disaggregates, and linked with clinical outcome. The phenotype of TRM cells was characterized by spectral flow cytometry (n=10 DLBCL and 2 reactive lymph node (rLN) samples) and single-cell RNA sequencing (scRNAseq) (n=12 aggressive B cell lymphomas and 4 rLN samples). An additional 51 samples from an external dataset were included to verify the single cell transcriptome findings. The flow cytometry and scRNAseq findings in DLBCL were extrapolated to additional B-NHL entities. Through analysis of these datasets we were able to characterize CD103+ TRM-like cells in DLBCL and other B-NHL entities and identified them to represent a prognostically favorable population with an activated/cytotoxic T cell phenotype.

  Study EGAS50000000943

• Investigating Host Genetic Risk Factors for Tuberculosis in Highly Endemic South African Populations

  It is apparent from our recent population genetic and admixture mapping work 1–3 that a substantial part of the ancestral input into the South African Coloured (SAC) population is from the San or Khoe groups, which are not well represented in publicly available genetic databases. It is a reasonable assumption that the SAC population in the Western Cape may have derived more genetic input from the present KhoeSan population in the Northern Cape than the San residing in Namibia. It was concluded by our group that it is indeed the southern African KhoeSan group, the ≠Khomani, that best represent the KhoeSan contribution seen in the SAC 4. More recently, it has been shown that there are ancestry related increases in TB susceptibility, especially with increased Bantu-speaking African and KhoeSan ancestry 5. We aim to elucidate the epidemiological and human-host genetic risk factors for TB and the immunological pathways modulating TB infection. The proposed sub-study will enroll people evaluated for TB at Northern Cape community health clinics and their contacts from their households and/or community. We will conduct a demographic interview, collect saliva, which will be used for genetic analysis, and blood to determine latent TB infection status and capture immunological responses to mycobacterial infections at the cellular level and RNA sequencing.

  Study EGAS00001007850

• INSPIRE - Investigator-initiated Phase 2 Study of Pembrolizumab Immunological Response in Metastatic Solid Tumors

  Whole exome and RNA sequencing source data for: Clouthier, D.L., Lien, S.C., Yang, S.Y.C. et al. An interim report on the investigator-initiated phase 2 study of pembrolizumab immunological response evaluation (INSPIRE). j. immunotherapy cancer 7, 72 (2019). https://doi.org/10.1186/s40425-019-0541-0 Bratman S*, Yang SYC*, Iafolla MAJ, Liu Z, Hansen AR, Bedard P, Lheureux S, Spreafico A, Razak AA, Shchegrova S, Louie M, Billings P, Zimmerman B, Sethi H, Aleshi A, Torti D, Marsh K, Eagles J, Cirlan I, Hanna Y, Clouthier D, Lien SC, Ohashi PS, Xu W, Siu LL, Pugh TJ. Personalized circulating tumor DNA analysis as a predictive biomarker in solid tumor patients treated with pembrolizumab. Nature Cancer. 2020:1:873-881. https://doi.org/10.1038/s43018- 020-0096-5. Yang SYC, Lien SC, Wang BX, Clouthier DL, Hanna Y, Cirlan I, Zhu K, Bruce JP, El Ghamrasni S, Iafallo MAJ, Oliva M, Hansen AR, Spreafico A, Bedard PL, Lheureux S, Razak A, Speers V, Berman HK, Aleshin A, Haibe-Kains B, Brooks DG, McGaha TL, Butler MO, Bratman SV, Ohashi PS, Siu LL, and Pugh TJ. Pan-cancer analysis of longitudinal metastatic tumors reveals genomic alterations and immune landscape dynamics associated with pembrolizumab sensitivity. Nature Comm 2021:12, 5137. https://doi.org/10.1038/s41467-021-25432-7

  Study EGAS00001003280

• RBMX functional retrocopy safeguards brain development

  Retrotransposition has generated thousands of intronless gene copies in mammalian genomes, yet their contribution to brain development and evolution remains largely unexplored. We investigated the role for RBMX retrocopy in shaping neurodevelopment and modulating disease. First, we reported 8 distinct pathogenic variants in RBMX, an X-linked splicing regulator, in 11 males presenting with intellectual disability, microcephaly, and cortical malformations. Then, we combined transcriptomic profiling, protein–protein and protein–RNA interaction studies both in human cellular models and mouse embryonic cortices to assess the functional redundancy between RBMX and its retrocopy RBMXL1. Finally, we use mouse genetics to dissect RBMX function and its compensation by RBMXL1 in corticogenesis. We show that RBMX pathogenic variants disrupt cortical development through both loss- and gain-of-function mechanisms. Surprisingly, despite severe phenotypes in humans, Rbmx-deficient mice display only mild cortical abnormalities, a discrepancy likely due to compensation by Rbmxl1, a retrocopy that arose independently in mice and humans. We demonstrate that RBMX and RBMXL1 share protein and RNA partners and act redundantly in brain development, with RBMXL1 buffering the impact of RBMX deficiency. Together, these findings establish that RBMXL1 functions as a functional paralog of RBMX safeguarding neurodevelopment by buffering deleterious variation. More broadly, these results identify retrocopies as active contributors to neurodevelopmental robustness and suggest that functional retrocopies may have facilitated the evolutionary diversification of the mammalian brain.

  Study EGAS50000001650

• Multimodal plasma and urinary cell-free DNA profiling improves risk stratification in newly diagnosed prostate cancer

  Prostate cancer (PCa) is a heterogeneous disease, impeding early detection and risk stratification. Liquid biopsies (LBx) enable minimally invasive tumor profiling, but circulating tumor-derived DNA (ctDNA) detection remains difficult, particularly in early-stage PCa. We developed a multimodal LBx approach combining genomic and epigenomic cell-free DNA (cfDNA) features in plasma and urine from newly diagnosed PCa patients to improve early characterization of PCa and risk stratification of aggressive disease. Plasma and urine samples from 55 localized PCa (lPCa) patients, 18 advanced PCa (aPCa) patients, and 36 cancer-free controls were subjected to low-coverage whole-genome sequencing and methylated DNA immunoprecipitation sequencing to assess fragmentation, chromosomal instability, and methylation in cfDNA. This complementary approach yielded a 45% ctDNA detection rate in newly diagnosed PCa. Major differences were observed between aPCa and controls, reflecting increasing signals with tumor progression. Epigenomic cfDNA features differentiated lPCa from aPCa, and ctDNA was detected in 46% of PCa patients with prostate-specific antigen <10 ng/ml, suggesting potential for risk stratification. This study highlights the value of multimodal LBx approaches for early characterization of primary PCa and identification of aggressive disease at initial diagnosis. Integration into clinical workflows could complement diagnostics and support personalized decision-making tailored to patients’ PCa risk profiles. 

  Study EGAS00001008195

• Whole genome sequencing of matched primary and metastatic acral melanomas

  Next generation sequencing has enabled systematic discovery of mutational spectra in cancer samples. Here we used whole genome sequencing to characterise somatic mutations and structural variation in a primary acral melanoma and its lymph node metastasis. Our data show that the somatic mutational rates in this acral melanoma sample pair was more comparable to the rate reported in cancer genomes not associated with mutagenic exposure than in the genome of a melanoma cell line or the transcriptome of melanoma short-term cultures. Despite the perception that acral skin is sun-protected, the dominant mutational signature in these samples is compatible with damage due to ultraviolet light exposure. A nonsense mutation in ERCC5 discovered in both the primary and metastatic tumours, could also have contributed to the mutational signature through accumulation of unrepaired dipyrimidine lesions. However, evidence of transcription-coupled repair was suggested by the lower mutational rate in the transcribed regions and expressed genes. The primary and the metastasis are highly similar at the level of global gene copy number alterations, loss of heterozygosity and single nucleotide variation (SNV). Furthermore, the majority of the SNVs in the primary tumour were propagated in the metastasis and one non-synonymous coding SNV and one splice site mutation appeared to arise de novo in the metastatic lesion.

  Study EGAS00001000169

• The Cardiogenics study

  Cardiogenics is an European collaborative project, it started in January 2007 and was funded by the European Commission through its Sixth Framework Program (reference LSHM -CT-2006-037593) until 2011. It involved several laboratories from France (INSERM UMRS937), Germany (Lübeck University; Regensburg University), United Kingdom (the Welcome Trust Sanger Institute; University of Leicester) and Sweden (Uppsala University Hospital) and was coordinated by a management group composed of scientists from these different institutions. Under the coordination of the management group the laboratories are still developing projects using Cardiogenics data sets. As part of the Cardiogenics project, RNA from monocytes and macrophages of 363 patients with coronary artery disease and 395 healthy individuals was prepared and genome-wide expression was assessed in both cell types using the Illumina HumanRef 8 v3 Beadchip . The DNA of all these individuals was genotyped using the Human 610 Quad custom arrays. Both data sets were used to conduct eQTL analyses. A subsample of subjects included in an allelic-imbalance substudy were also genotyped using the Illumina Human Custom 1.2M array. Detailed description of these datasets can be found in Shah S et al. Circ Cardiovs Genet 2011; Rotival M et al. Plos Genet 2011, Almlöf et al. Plos One 2012 and Garnier S et al. Plos Genet 2013

  Study EGAS00001000411

• A Cancer Cell-Line Titration Series for Evaluating Somatic Classification

  Accurate detection of somatic single nucleotide variants and small insertions and deletions from DNA sequencing experiments of tumour-normal pairs is a challenging task. Tumour samples are often contaminated with normal cells confounding the available evidence for the somatic variants. Furthermore, tumours are heterogeneous so sub-clonal variants are observed at reduced allele frequencies. We present here a cell-line titration series dataset that can be used to evaluate somatic variant calling pipelines with the goal of reliably calling true somatic mutations at low allele frequencies. Cell-line DNA was mixed with matched normal DNA at 8 different ratios to generate samples with known tumour cellularities, and exome sequenced on Illumina HiSeq to depths of >300x. The data was processed with several different variant calling pipelines and verification experiments were performed to assay >1,500 somatic variant candidates using Ion Torrent PGM as an orthogonal technology. By examining the variants called at varying cellularities and depths of coverage, we show that the best performing pipelines are able to maintain a high level of precision at any cellularity. In addition, we estimate the number of true somatic variants undetected as cellularity and coverage decrease. Our cell-line titration series dataset, along with the associated verification results, is effective for this evaluation and will serve as a valuable dataset for future somatic calling algorithm development.

  Study EGAS00001001016

• Genomic characterization of retinoblastoma (targeted sequencing)

  Retinoblastoma is the most frequent intraocular malignancy in children, originating from a maturing cone precursor in the developing retina. Little is known on the molecular basis underlying the biological and clinical behavior of this cancer. Here, using multi-omics data, we demonstrate the existence of two retinoblastoma subtypes. Subtype 1, of earlier onset, includes most of the heritable forms. It harbors few genetic alterations other than the initiating RB1 inactivation and corresponds to differentiated tumors expressing mature cone markers. By contrast, subtype 2 tumors harbor frequent recurrent genetic alterations including MYCN-amplification. They express markers of less differentiated cone together with neuronal/ganglion cell markers with marked inter- and intra-tumor heterogeneity. The cone dedifferentiation in subtype 2 is associated with stemness features including low immune and interferon response, E2F and MYC/MYCN activation and a higher propensity for metastasis. The recognition of these two subtypes, one maintaining a cone-differentiated state, and the other, more aggressive, associated with cone dedifferentiation and expression of neuronal markers, opens up important biological and clinical perspectives for retinoblastomas. In this study, we performed targeted sequencing of the exonic regions of RB1, BCOR and ARID1A on 37 samples from the series of 102 retinoblastomas. Among them, 23 samples were not subjected to whole-exome sequencing, and 14 samples have the whole-exome sequencing data (deposited in EGAS00001005248).

  Study EGAS00001005550

• UK10K NEURO ASD MGAS

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The MGAS (Molecular Genetics of Autism Study) samples are from a clinical sample seen by specialists at the Maudsley hospital and who have had detailed phenotypic assessments with ADI-R and ADOS.For further information on this cohort please contact Patrick Bolton (patrick.bolton@kcl.ac.uk).

  Study EGAS00001000113

• Papuan Y chromosome Diversity Panel

  Island Southeast Asia and Oceania host one of the world’s richest assemblages of human phenotypic, linguistic and cultural diversity. Despite this, the region’s male genetic lineages are globally among the last to remain unresolved. We compiled ~9.7 Mb of Y chromosome sequence from a diverse sample of over 380 men from this region, including 152 first reported here. The granularity of this dataset allows us to fully resolve and date the regional Y chromosome phylogeny. This new high-resolution tree confirms two main population bursts: multiple rapid diversifications following the region’s initial settlement ~50 kya, and extensive expansions <6 kya. Notably, ~40-25 kya the deep rooting local lineages of C-M130, M-P256, and S-B254 show almost no further branching events in Island Southeast Asia, New Guinea and Australia, matching a similar pause in diversification seen in maternal mitochondrial DNA lineages. The main local lineages start diversifying ~25 kya, at the time of the Last Glacial Maximum. This improved Y chromosome topology highlights localized events with important historical implications, including pre-Holocene contact between Mainland and Island Southeast Asia, potential interactions between Australia and the Papuan world, and a sustained period of diversification following the flooding of the ancient Sunda and Sahul continents as the insular landscape observed today formed. The high-resolution phylogeny of the Y chromosome presented here thus enables a detailed exploration of

  Study EGAS00001006025

• Enrichment of oral-derived bacteria in inflamed colorectal tumors and distinct associations of Fusobacterium in the mesenchymal subtype

  The association between colorectal cancer (CRC) clinical variables and Fusobacterium, but not other intra-tumoral bacteria, has been extensively studied. Here we leveraged whole-transcriptome sequencing from 807 CRC tumor samples from the AVANT phase III trial to dually characterize tumor gene expression and intra-tumoral bacteria. After stringent filtering, 74 high-confidence taxa were identified. 17 of these species, including 4 Fusobacterium spp., were classified as orally-derived and had a robust signal within right-sided, MSI-H, and BRAF mutant tumors. Across consensus molecular subtypes (CMS), integration of Fusobacterium animalis presence and tumor gene expression revealed that F. animalis had the greatest number of associations in mesenchymal CMS4 tumors, despite an overall lower prevalence than in immune CMS1 tumors. Pathway analysis within CMS4 revealed that F. animalis, but not other highly prevalent species, was uniquely associated with pathways for collagen degradation and formation as well as IL-6 and IL-1 cytokine signaling. These associations could explain in part why Fusobacterium has been previously associated with reduced survival in mesenchymal populations. Overall, this study identified oral-derived bacteria, including Fusobacterium, as broadly more prevalent in inflamed CRC tumors compared to the other subtypes, and the association of intra-tumoral bacteria on tumor gene expression is context- and species-specific.

  Study EGAS00001006757

• Novel CNV contribution to schizophrenia from a genome wide study of 41,321 subjects

  Genomic copy number variants (CNVs) have been strongly implicated in the etiology schizophrenia (SCZ). However, apart from a small number of risk variants, elucidation of the CNV contribution to risk has been difficult due to the very low frequencies of risk alleles, all occurring in less than 1% of patients. We sought to address this obstacle through a collaborative effort unprecedented in psychiatry in which we applied a centralized analysis pipeline to a large SCZ cohort of 21,094 cases and 20,227 controls. A global enrichment of CNV burden is observed in SCZ patients (OR=1.11, P=5.7e-15), and persists after excluding loci implicated in previous studies (OR=1.07, P=1.7e-6). CNV Burden is also enriched for genes associated with synaptic function (OR = 1.68, P = 2.8e-11) and neurobehavioral phenotypes in mouse (OR = 1.18, P= 7.3e-05). While only previously implicated loci surpass strict genome-wide testing correction, we identify “novel” CNVs (BH-FDR < 0.05) that confer both protective and risk effects on SCZ and are predominantly mediated by non-allelic homologous recombination (NAHR). Our combined data find genome-wide significant support for eight loci, including 1q21.1, 2p16.3 (NRXN1), 3q29, 7q11.2, 15q13.3, distal 16p11.2, proximal 16p11.2 and 22q11.2.

  Study EGAS00001001960

• WES of 2 human osteosarcoma and corresponding cell lines

  Progress in the systemic control of osteosarcoma has been limited over the past decades thus indicating the urgent clinical need for the development of novel treatment strategies. Therefore, we have recently developed new preclinical models to study promising novel agents for the treatment of pediatric osteosarcoma. The checkpoint kinase (chk) inhibitor prexasertib (LY2606368) and its salt form (LSN2940930) have recently been shown to be active in adult and pediatric malignancies, including sarcoma. We have now tested the potency of prexasertib in clonogenic survival assays in two new lines of primary patient-derived osteosarcoma cells and in two established osteosarcoma cell lines as a single agent and in combination with cisplatin and the poly ADP-ribose polymerase (PARP) inhibitor talazoparib. Prexasertib alone results in strongly reduced clonogenic survival at low nanomolar concentrations and acts by affecting cell cycle progression, induction of apoptosis and induction of double-stranded DNA breakage at concentrations that are well below clinically tolerable and safe plasma concentrations. In combination with cisplatin and talazoparib, prexasertib acts in a synergistic fashion. Chk1 inhibition by prexasertib and its combination with the DNA damaging agent cisplatin and the PARP-inhibitor talazoparib thus emerges as a potential new treatment option for pediatric osteosarcoma which will now have to be tested in preclinical primary patient derived in vivo models and clinical studies.

  Study EGAS00001003923

• Human Pancreatic Islet RNAseq - Lund

  Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, many key tissues and cell-types required for appropriate functional inference are absent from large-scale resources such as ENCODE and GTEx. We explored the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using RNA-Seq and genotyping data from 420 islet donors. We find: (a) eQTLs have a variable replication rate across the 44 GTEx tissues (<73%), indicating that our study captured islet-specific cis-eQTL signals; (b) islet eQTL signals show marked overlap with islet epigenome annotation, though eQTL effect size is reduced in the stretch enhancers most strongly implicated in GWAS signal location; (c) selective enrichment of islet eQTL overlap with the subset of T2D variants implicated in islet dysfunction; and (d) colocalization between islet eQTLs and variants influencing T2D or related glycemic traits, delivering candidate effector transcripts at 23 loci, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in tissues of greatest disease-relevance while expanding our mechanistic insights into complex traits association loci activity with an expanded list of putative transcripts implicated in T2D development.

  Study EGAS00001004042

• Phenotyping data on human pancreatic islets from 191 donors - Lund

  Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, many key tissues and cell-types required for appropriate functional inference are absent from large-scale resources such as ENCODE and GTEx. We explored the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using RNA-Seq and genotyping data from 420 islet donors. We find: (a) eQTLs have a variable replication rate across the 44 GTEx tissues (<73%), indicating that our study captured islet-specific cis-eQTL signals; (b) islet eQTL signals show marked overlap with islet epigenome annotation, though eQTL effect size is reduced in the stretch enhancers most strongly implicated in GWAS signal location; (c) selective enrichment of islet eQTL overlap with the subset of T2D variants implicated in islet dysfunction; and (d) colocalization between islet eQTLs and variants influencing T2D or related glycemic traits, delivering candidate effector transcripts at 23 loci, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in tissues of greatest disease-relevance while expanding our mechanistic insights into complex traits association loci activity with an expanded list of putative transcripts implicated in T2D development.

  Study EGAS00001004056

• Aggressive genomic features in clinicallyindolent primary HHV8-negative effusion-based lymphoma

  Human herpesvirus 8-negative effusion-based lymphoma (HHV8-negative EBL) is a distinct lymphoma entity related to fluid overload states caused by underlying medical conditions. Due to a phenotypical resemblance with secondary effusion of diffuse large B-cell lymphoma (DLBCL), HHV8-negative EBL is easily misdiagnosed. However, in contrast to DLBCL, HHV8-negative EBL follows a mild clinical course, underpinning the importance to differentiate between the two entities. Only limited molecular data of HHV8-negative EBL has been reported and characteristic genomic patterns are lacking. We performed extensive genomic profiling on 8 HHV8-negative EBL cases. A mean number of 33.6 copy number gains and losses per case were found, with recurrent focal deletions of chromosome 3p14.2 and 6q21 in 6/8 cases, encompassing tumor supressor genes FHIT and PRDM1, respectively. We observed a high rate of mutations across genes frequently involved in DLBCL, including patterns of somatic hypermutation. The most abundant hotspot mutation was found in MYD88L265P (3/8 cases). In addition, 5/8 cases had one or more translocations, involving MYC, BCL2, BCL6, TP63, EXOC2 and KMT2D, with Ig and non-Ig partners. We conclude that HHV8-negative EBL is characterized by genomic features that are commonly found in more aggressive B-cell lymphomas, but despite this follow an indolent clinical course.

  Study EGAS00001002743

• RNA-sequencing of six Pilocytic astrocytoma tumors

  Pilocytic astrocytoma (PA) is the most common pediatric brain tumor. A recurrent feature of PA is deregulation of the mitogen activated protein kinase (MAPK) pathway most often through KIAA1549-BRAF fusion, but also by other BRAF- or RAF1-gene fusions and point mutations (e.g. BRAFV600E). These features may serve as diagnostic and prognostic markers, and also facilitate development of targeted therapy. The aim of this study was to characterize the genetic alterations underlying the development of PA tumor in six cases, and evaluate methods for fusion oncogene detection. Using a combined analysis of RNA sequencing and copy number variation data we identified a new BRAF fusion involving the 5’ gene fusion partner GTF2I (7q11.23), not previously described in PA. The novel GTF2I-BRAF 19-10 fusion was found in one case, while the other five cases harbored the frequent KIAA1549-BRAF 16-9 fusion gene. Comparing fusion detection methods, Fluorescence in situ hybridization with BRAF break apart probe was the most sensitive method for detection of the two different BRAF rearrangements (GTF2I-BRAF and KIAA1549-BRAF). Our finding of a novel BRAF fusion in PA further emphasis the important role of B-Raf in tumorigenesis of these tumor types. Also, the growing list of BRAF/RAF gene fusions suggests these to be informative tumor markers in molecular diagnostics, which could guide future treatment strategies.

  Study EGAS00001002149

• Illumina HumanOmniExpress genotyping data from the TEENAGE (TEENs of Attica: Genes and Environment) study.

  The TEENAGE study is a cross-sectional study comprising healthy Greek adolescents. The study was conducted according to the guidelines specified in the Declaration of Helsinki and all procedures involving human subjects were approved by the Institutional Review Board of Harokopio University and the Greek Ministry of Education, Lifelong Learning and Religious Affairs. The TEENAGE study target population comprised adolescent students aged 13 to 15 years attending the first three classes of public secondary schools located in the wider Athens area of Attica. The regions across Attica from which students were recruited were physically similar and all urban in nature. Sampling of regions and schools selected to participate was random. Students were enrolled from the wider urban region to ensure that they were coming from families with varying educational background and socio-economic status. Moreover, in order to exclude any vast economic disparities, all participating schools were exclusively public schools, which are free of charge. A total of twenty-four schools participated in the study out of 285 junior secondary public schools in Attica. Anthropometric, medical and dietary information was obtained. DNA samples were genotyped using Illumina HumanOmniExpress BeadChips (Illumina, San Diego, CA, USA) at the Wellcome Trust Sanger Institute, Hinxton, UK. Genotype calling algorithm used was Illuminus (Teo et al., 2007).

  Study EGAS00001000996

• Finnish_population_cohort_genotyping_B

  The FINRISK cohorts comprise the respondents of representative, cross-sectional populationsurveys that are carried out every 5 years since 1972, to assess the risk factors of chronicdiseases (e.g. CVD, diabetes, obesity, cancer) and health behavior in the working agepopulation, in 3-5 large study areas of Finland. DNA samples were collected in the followingsurvey years: 1987, 1992, 1997, 2002, 2007, and 2012. The MONICA and EHES (EU)procedures were applied in phenotype collection (cf. MORGAM) and a wide spectrum oflaboratory tests was carried out from serum and plasma samples. Background information onsocioeconomic status, medical history, diet, exercise, measured anthropometric measures,etc. was collected by questionnaires and during a clinical visit. Plasma/serum samples arestill available for the 2002-2012 cohorts. The cohort sizes are 6000-8800 per survey. Thecohorts have been followed up by linking them to the national hospital discharge register,causes-of-death register and cancer register.This project is an extension to previous efforts to build a catalogue of Finnish genome widedata on population-based Finsrisk samples with rich phenotypic characterisations and healthregistry link-up. These samples will extend the current Sequencing Initiative Suomi (SISu)samples with a combination of genotyping using Illumina HumanCoreExome array and SISu-based imputation. This will lead to high confidence common and low frequency variantcatalogue. The project will be funded by Aarno Palotie’s remaining faculty fundscomplemented by Finnish funding from FIMM.

  Study EGAS00001001047

• Evolution of DNA repair defects during malignant progression of low-grade gliomas after temozolomide treatment.

  Temozolomide (TMZ) increases the overall survival of patients with glioblastoma (GBM), but its role in the clinical management of diffuse low-grade gliomas (LGG) is still being defined. DNA hypermethylation of the O 6 -methylguanine-DNA methyltransferase (MGMT) promoter is associated with an improved response to TMZ treatment, while inactivation of the DNA mismatch repair (MMR) pathway is associated with therapeutic resistance and TMZ-induced mutagenesis. We previously demonstrated that TMZ treatment of LGG induces driver mutations in the RB and AKT-mTOR pathways, which may drive malignant progression to secondary GBM. To better understand the mechanisms underlying TMZ-induced mutagenesis and malignant progression, we explored the evolution of MGMT methylation and genetic alterations affecting MMR genes in a cohort of 34 treatment-naïve LGGs and their recurrences. Recurrences with TMZ-associated hypermutation had increased MGMT methylation compared to their untreated initial tumors and higher overall MGMT methylation compared to TMZ-treated non-hypermutated recurrences. A TMZ-associated mutation in one or more MMR genes was observed in five out of six TMZ-treated hypermutated recurrences. In two cases, pre-existing heterozygous deletions encompassing MGMT, or an MMR gene, were followed by TMZ-associated mutations in one of the genes of interest. These results suggest that tumor cells with methylated MGMT may undergo positive selection during TMZ treatment in the context of MMR deficiency.

  Study EGAS00001001179

• JMML targeted sequencing (2013)

  Juvenile myelomonocytic leukemia (JMML) is a rare myelodysplastic myeloproliferative neoplasm of early childhood initiated by germline or somatic RAS-activating mutations (Chang, Dvorak et al. 2014). Clinical evolution is heterogeneous, with acute myeloid leukemia (AML) transformation in 1/3 of cases and frequent relapses after hematopoietic-stem-cell transplantation (HSCT), although some patients experience spontaneous remission (Locatelli, Nollke et al. 2005; Matsuda, Shimada et al. 2007; Takagi, Piao et al. 2013; Locatelli and Niemeyer 2015). The initiating lesion incompletely predicts outcome, suggesting the involvement of additional mutations. Genetic profiling of a large JMML cohort (n=118), including whole-exome sequencing (WES) in 30 cases, uncovered additional genetic abnormalities in 56 cases (47%). Somatic events were rare (0.38/case/Mb) and restricted to sporadic (47/78; 60%) or NF1-associated (8/8; 100%) JMML. Multiple concomitant genetic hits were identified targeting the RAS pathway in 13/78 (17%), RAC2D63V, which activates the PI3K/PDK1/AKT and the mTORC2 pathways, and a novel pathway for JMML, the polycomb repressive complex 2 (PRC2), in 26/78 (33%) of the sporadic JMML cases. Interestingly, PRC2 imbalance was recently shown to synergize with RAS activation (De Raedt, Beert et al. 2014). JMML outcome was associated to the mutational profile, suggesting a dose-dependent effect for RAS-pathway activation, distinguishing very aggressive JMML rapidly progressing to AML.

  Study EGAS00001001324

• GENETIC HISTORY OF ITALY

  Recent scientific literature has highlighted the relevance of population genetic studies both for disease association-mapping in admixed populations and for understanding the history of human migrations. Deeper insight into the history of Italian population is critical for understanding the peopling of Europe. Because of its crucial position at the centre of the Mediterranean basin, the Italian peninsula has experienced a complex history of colonization and migration, whose genetic signatures are still present in contemporary Italians. In this study, we investigated genomic variation in the Italian population using 2.5 million single nucleotide polymorphisms (SNPs) in a sample of more than 300 unrelated Italian subjects with well-defined geographical origins. We combined several analytical approaches to interpret genome-wide data on 1,272 individuals from European, Middle Eastern, and North African populations. We detected three major ancestral components contributing different proportions across the Italian peninsula, and signatures of continuous gene flow within Italy, which has produced remarkable genetic variability among contemporary Italians. In addition, we have extracted novel details about the Italian population’s ancestry: we identified genetic signatures of major historical events in Europe and the Mediterranean basin from the Neolithic (e.g., peopling of Sardinia) to recent times (e.g., ‘barbarian invasion’ of Northern and Central Italy). These results are valuable for further genetic, epidemiological and forensic studies in Italy and in Europe.

  Study EGAS00001001458

• Exome and RNA sequencing of relapsed TCF3-PBX1 t(1;19) acute lymphoblastic leukemia

  TCF3-PBX1 (E2A-PBX1) is a recurrent gene fusion in B-cell precursor lymphoblastic leukemia (BCP-ALL), which is caused by the translocation t(1;19)(q23;p13). TCF3-PBX1 BCP-ALL patients typically benefit from chemotherapy; however, many relapse and subsequently develop resistant disease with few effective treatment options. Mechanisms driving disease progression and therapy resistance have not been studied in TCF3-PBX1 BCP-ALL. Here, we aimed to identify novel treatment options for TCF3-PBX1 BCP-ALL by profiling leukemia cells from a relapsed patient, and determine molecular mechanisms underlying disease pathogenesis and progression. By drug sensitivity testing of leukemic blasts from the index patient, control samples and TCF3-PBX1 positive and negative BCP-ALL cell lines, we identified the phosphatidylinositide 3-kinase delta (p110δ) inhibitor idelalisib as an effective treatment for TCF3-PBX1 BCP-ALL. This was further supported by evidence showing TCF3-PBX1 directly regulates expression of PIK3CD, the gene encoding p110δ. Other somatic mutations to TP53 and MTOR, as well as aberrant expression of CXCR4, may influence additional drug sensitivities specific to the index patient and accompanied progression of the disease. Our results suggest that idelalisib is a promising treatment option for patients with TCF3-PBX1 BCP-ALL while other drugs could be useful depending on the genetic context of individual patients.

  Study EGAS00001001876

• Automated system for scoring hematoxylin and eosin-stained ovarian cancer sections by identifying single cells uncovered that stromal cell ratio is a significant predictor for overall survival and progression-free survival.

  Concerted efforts in genomic studies have revealed profound insights in prognostic ovarian cancer subtypes. On the other hand, abundant histology slides have been generated to date, yet their uses remain very limited and largely qualitative. Our goal is to develop automated histology analysis as an alternative subtyping technology for ovarian cancer that is cost-efficient and do not rely on DNA quality. We develop an automated system for scoring hematoxylin and eosin-stained (H&E) primary tumour sections of 91 late-stage ovarian cancer to identify single cells including cancer and stromal cells. We demonstrated high accuracy of our system based on expert pathologists’ scores (cancer=97.1%, stromal=89.1%) as well as compared to immunohistochemistry scoring (correlation=0.87). Quantitative stromal cell ratio is significantly associated with poor overall survival after controlling for clinical parameters including debulking status and age (multivariate analysis p=0.0021, HR=2.54, CI=1.40-4.60) and progression-free survival (multivariate analysis p=0.022, HR=1.75, CI=1.09–2.82). We demonstrate how automated image analysis enables objective quantification of microenvironmental composition of ovarian tumours. Our analysis reveals a strong effect of the tumour microenvironment on ovarian cancer progression and highlights the potential of therapeutic interventions that target the stromal compartment or cancer-stroma signalling in the stroma-high, late-stage ovarian cancer subset.

  Study EGAS00001001694

• Transcriptome sequencing of gingivo-buccal oral squamous cell carcinoma for integrative analysis: alterations in expression of genes attributable to methylation changes

  Gingivo-buccal oral squamous cell carcinoma (OSCC-GB) is the most common cancer among men in India and is associated with high mortality. Although OSCC-GB is known to be quite different from tongue cancer in its genomic presentation and its clinical behavior, it is treated identically as tongue cancer. Predictive markers of prognosis and therapy that are specific to OSCC-GB are, therefore, required. To achieve this, we have carried out transcriptome (n=72) profiling of paired tumor-normal tissues collected from OSCC-GB patients from India. RNA sequencing was performed on tumor and paired adjacent normal samples from 72 OSCC-GB patients. From each patient, total RNA was extracted from tumor and normal tissue samples using AllPrep DNA/RNA Mini Kit (QIAGEN). The quality and concentration of isolated total RNA were checked using Agilent 2100 Bioanalyzer and NanoDrop 2000 (Thermo Fisher Scientific). The OD260/OD280 ratio was ≥ 2 and RNA Integrity Number (RIN) was ≥ 7.0 for all sequenced samples. Ribosomal RNA (rRNA) was removed from the RNA samples, using Ribo-Zero Magnetic Kit (Epicentre). Sequencing libraries were prepared from rRNA depleted samples using TrueSeq RNA Sample Preparation Kit (Illumina). Each triplex cDNA library pool was sequenced as 100bp paired-end on HiSeq-2000 or HiSeq-2500 (Illumina). Protocols suggested by the manufacturers were used for all assays.

  Study EGAS00001003893

• GWAS data on human pancreatic islets from 191 donors - Lund

  Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, many key tissues and cell-types required for appropriate functional inference are absent from large-scale resources such as ENCODE and GTEx. We explored the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using RNA-Seq and genotyping data from 420 islet donors. We find: (a) eQTLs have a variable replication rate across the 44 GTEx tissues (<73%), indicating that our study captured islet-specific cis-eQTL signals; (b) islet eQTL signals show marked overlap with islet epigenome annotation, though eQTL effect size is reduced in the stretch enhancers most strongly implicated in GWAS signal location; (c) selective enrichment of islet eQTL overlap with the subset of T2D variants implicated in islet dysfunction; and (d) colocalization between islet eQTLs and variants influencing T2D or related glycemic traits, delivering candidate effector transcripts at 23 loci, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in tissues of greatest disease-relevance while expanding our mechanistic insights into complex traits association loci activity with an expanded list of putative transcripts implicated in T2D development.

  Study EGAS00001004044

• Comparison of 133/144bp dominant phenotype vs 166bp dominant phenotype.

  Background: Previous studies found that cell-free DNA (cfDNA) generated from tumors was shorter than that from healthy cells, and selecting short cfDNA could enrich for tumor cfDNA and improve its usage in early cancer diagnosis and treatment monitoring; however, the underlying mechanism of shortened tumor cfDNA was still unknown, which potentially limits its further clinical application.Results: Using targeted sequencing of cfDNA in a large cohort of solid tumor patient, sequencing reads harboring tumor-specific somatic mutations were isolated to examine the exact size distribution of tumor cfDNA. For the majority of studied cases, 166bp remained as the peak size of tumor cfDNA, with tumor cfDNA showing an increased proportion of short fragments (100-150bp). Less than 1% of cfDNA samples were found to be peaked at 134/144bp and independent of tumor cfDNA purity. Using whole-genome sequencing of cfDNA, we discovered a positive correlation between cfDNA shortening and the magnitude of chromatin inaccessibility, as measured by transcription, DNase I hypersensitivity, and histone modifications. Tumor cfDNA shortening occurred simultaneously at both 5’ and 3’ ends of the DNA wrapped around nucleosomes. Conclusions: Tumor cfDNA shortening exhibited two distinctive modes. Tumor cfDNA purity and chromatin inaccessibility were contributing factors but insufficient to trigger a global transition from 166bp dominant to 134/144bp dominant phenotype.

  Study EGAS00001004135

• Mutational Landscape of Grey Zone Lymphoma

  The mutational landscape of grey zone lymphoma (GZL) has not yet been established and differences to related entities are mostly unknown. Here, we studied coding sequence mutations of 50 EBV-negative GZL and 20 polymorphic EBV-positive DLBCL NOS (poly-EBV-L) in comparison to classical Hodgkin lymphoma (cHL), primary mediastinal B cell lymphoma (PMBCL), and diffuse large B cell lymphoma (DLBCL). Paired exomes of 21 GZL and 7 poly-EBV-L cases were analyzed as a discovery cohort followed by targeted sequencing of 217 genes in an extension cohort of 29 GZL and 13 poly-EBV-L cases. GZL cases with thymic niche involvement displayed a mutation profile closely resembling cHL and PMBCL, with SOCS1 (45%), B2M (45%), TNFAIP3 (35%), GNA13 (35%) and NFKBIA (29%) being the most recurrently mutated genes. In contrast, GZL cases without thymic niche involvement (N=18) had a significantly distinct pattern, enriched in mutations related to apoptosis defects (TP53 (39%), BCL2 (28%), BIRC6 (22%)) and depleted in GNA13, XPO1 or NFKB signaling pathway mutations (TNFAIP3, NFKBIE, IKBKB, NFKBIA). Poly-EBV-L cases presented a distinct mutational profile including STAT3, PRKDC mutations and a significantly lower coding-mutation load in comparison to EBV-negative GZL. Our study highlights characteristic mutational patterns in GZL associated with presentation in the thymic niche suggesting a common cell of origin and disease evolution with related anterior mediastinal lymphomas.

  Study EGAS00001004482

• Similarity and diversity of the tumor microenvironment in multiple metastases: critical implications for overall and progression-free survival of high-grade serous ovarian cancer.

  The tumor microenvironment is pivotal in influencing cancer progression and metastasis. Different cells co-exist with high spatial diversity within a patient, yet their combinatorial effects are poorly understood. We investigate the similarity of the tumor microenvironment of 192 local metastatic lesions in 61 ovarian cancer patients. An ecologically inspired measure of microenvironmental diversity derived from multiple metastasis sites is correlated with clinicopathological characteristics and prognostic outcome. We demonstrate a high accuracy of our automated analysis across multiple sites. A low level of similarity in microenvironmental composition is observed between ovary tumor and corresponding local metastases (stromal ratio r = 0.30, lymphocyte ratio r = 0.37). We identify a new measure of microenvironmental diversity derived from Shannon entropy that is highly predictive of poor overall (p = 0.002, HR = 3.18, 95% CI = 1.51-6.68) and progression-free survival (p = 0.0036, HR = 2.83, 95% CI = 1.41-5.7), independent of and stronger than clinical variables, subtype stratifications based on single cell types alone and number of sites. Although stromal influence in ovary tumors is known to have significant clinical implications, our findings reveal an even stronger impact orchestrated by diverse cell types. Quantitative histology-based measures can further enable objective selection of patients who are in urgent need of new therapeutic strategies such as combinatorial treatments targeting heterogeneous tumor microenvironment.

  Study EGAS00001002065

• Prebiotic inulin-type fructans induce specific changes in the human gut microbiota

  ObjectiveContrary to the long-standing prerequisite of inducing selective (i.e. bifidogenic) effects, recent findings suggest that prebiotic interventions lead to ecosystem-wide microbiota shifts. Yet, a comprehensive characterization of this process is still lacking. Here, we apply 16S rDNA microbiota profiling and matching (GC-MS) metabolomics to assess the consequences of inulin fermentation both on the composition of the colon bacterial ecosystem and fecal metabolites profiles.DesignFecal samples collected during a double blind, randomized, cross-over intervention study (NCT02548247) set up to assess the effect of inulin consumption on stool frequency in healthy adults with mild constipation were analyzed. Fecal microbiota composition and metabolite profiles were linked to the study’s clinical outcome as well as to quality-of-life measurements recorded.ResultsWhile fecal metabolite profiles were not significantly altered by inulin consumption, our analyses did detect a modest effect on global microbiota composition. At the same time, specific inulin-induced changes in relative abundances of Anaerostipes, Bilophila, and Bifidobacterium were identified. The observed decrease in Bilophila abundances following inulin consumption was associated with both softer stools and a favorable change in constipation-specific quality of life measures.ConclusionsEcosystem-wide analysis of the effect of a dietary intervention with prebiotic inulin-type fructans on the colon microbiota revealed that this effect is specifically associated to three genera, one of which (Bilophila) representing a promising novel target for mechanistic research.

  Study EGAS00001002173

• Antisense long non-coding RNAs are deregulated in skin tissue of patients with systemic sclerosis

  Systemic sclerosis (SSc) is an autoimmune disease characterized by fibrosis of skin and multiple organs of which the pathogenesis is poorly understood. Here we studied differentially expressed coding and non-coding genes in relation to SSc pathogenesis with a specific focus on antisense non-coding RNAs. Skin biopsy-derived RNAs from fourteen early SSc patients and six healthy individuals were sequenced with ion-torrent and analysed using DEseq2. Overall, 4901 genes with a fold change >1.5 and a false discovery rate < 5% were detected in patients versus controls. Upregulated genes clustered in immunological, cell adhesion and keratin-related processes. Interestingly, 676 deregulated non-coding genes were detected, 257 of which were classified as antisense genes. Sense genes expressed opposite of these antisense genes were also deregulated in 42% of the observed sense-antisense gene pairs. The majority of the antisense genes had a similar effect sizes in an independent North American dataset with three genes (CTBP1-AS2, OTUD6B-AS1 and AGAP2-AS1) exceeding the study-wide Bonferroni-corrected ρ-value (PBonf<0.0023, Pcombined = 1.1x10-9, 1.4x10-8, 1.7x10-6, respectively). In this study, we highlight that together with coding genes, (antisense) long non-coding RNAs are deregulated in skin tissue of SSc patients suggesting a novel class of genes involved in pathogenesis of SSc.

  Study EGAS00001002751

• Genomic analysis reveals novel secondary drivers and progression pathways in skin basal cell carcinoma

  Basal cell carcinoma of the skin (BCC) is the most common malignant neoplasm in humans. BCC is primarily driven by the aberrant activation of the Sonic Hedgehog (Hh) pathway. However, its extensive phenotypical variation remains to be explained. The genetic profiling of 293 BCCs revealed the highest mutation rate observed in cancer (65 Mutations/Mb), with strong prevalence of UV-light signature mutations. 85% of BCCs harbored mutations in Hh pathway genes: mutually exclusive PTCH1 (73%) and SMO (20%) (P=6.6x10-8), SUFU (8%), and in TP53 (61%). 85% of BCCs also harbored secondary driver mutations in other genes implicated in BCC tumorigenesis. Recurrent secondary driver mutations were observed in MYCN (30%), PPP6C (15%), STK19 (10%), LATS1 (8%), ERBB2 (4%), PIK3CA (2%), RAC1 (1%) and N/K/H-RAS (2%). Loss of function (LoF) and deleterious missense mutations were observed in PTPN14 (23%), RB1 (8%) and FBXW7 (5%). In line with the mutational profiles detected by RNAseq, we observed activation of the Hh pathway as well as upregulation of target genes of the Hippo-YAP pathway and activation of MYCN target genes. The functional analysis of the novel tumorigenic driver mutations in MYCN, PTPN14 and LATS1 demonstrates their relevance in BCC tumorigenesis and provides an expanded molecular understanding of BCC.

  Study EGAS00001001540

• Identification_of_cardiovascular_biomarkers_through_an_integrative_omics_approach

  Recent GWAS studies have made extensive use of large eQTL data sets to functionallyannotate index SNPs. With a large number of association signals located outside codingregions there has been an intense search among sequence variants affecting geneexpression at the transcriptional level. However, little progress has been made in mappingregulatory variants that affect protein levels at the translational or post-translational level. It isnow possible to undertake a protein QTL scan for focused sets of e.g. oxidized proteins bymass spectrometry. We have established a collaboration with a longitudinal, family-basedstudy in France, the Stanislas cohort, which comprises circa 1000 nuclear families (4,295individuals) and has follow up data for 10 years (three visits). We have undertaken a pilotstudy in a focus set of 257 subjects from 79 families with the aim to integrate GWAS,transcriptomic and DNA methylation data with proteomic data on a set of 100 proteinsmeasured in PBMCs. We have already generated GWAS data using Illumina's core-exomechip as well as DNA methylation profiles with the 450K array. We propose to use RNA seq togenerate transcriptomic data of the corresponding PBMCs.This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001000711

• Comprehensive analysis of atypical teratoid rhabdoid tumour (ATRT) using genomic, epigenomic and transcriptomic techniques.

  ATRTs (Atypical teratoid rhabdoid tumours) represent one of the most aggressive pediatric brain cancers, but intriguingly exhibit few other recurrently mutated loci except for SMARCB1/hSNF5. We integrated whole genome (n=15), exome, copy number, gene expression and methylation analyses to comprehensively interrogate 64 ATRTs and observed that structural events were relatively frequent in the ATRT genome (~3 tumour). In addition to SMARCB1, which was targeted by structural events in a majority (49/64) of tumours, recurrent structural alterations targeting the LRP1B, CDH13, BCR and MKL1 loci were observed. We observed novel translocation events, including 2 targeting SMARCB1 which were detected only by NGS analyses. Significantly, integration of gene expression and methylation profiles with genomic analyses revealed ATRT comprise two sub-groups with distinct clinical and genetic features. Group 1 tumours were characterized by supra-tentorial brain location, focal intragenic alterations of SMARCB1 and a pro-neural gene expression signature with evidence of NOTCH pathway activation. While Group 2 ATRTs were characterized by infra-tentorial brain location, broad SMARCB1 alterations and an activation of the BMP signalling pathway. Our findings highlight epigenetic mechanisms as important determinants of ATRT tumour phenotypes, and demonstrate for the first time that ATRTs arising in different anatomical compartments comprise distinct molecular and therapeutic sub-groups.

  Study EGAS00001000506

• Cell fate mapping of human glioblastoma reveals an invariant stem cell hierarchy pre- and post-treatment

  Human glioblastomas (GBMs) are thought to harbour a subpopulation of glioblastoma stem cells (GSCs) that initiate tumour formation and regrowth following treatment. However, the origin of their proliferative heterogeneity during tumour growth pre- and post-treatment remains poorly understood. Here we study the clonal evolution of DNA barcoded GBMs following their serial propagation in xenograft mouse models to define the fate behaviour of individual GBM cells. Through quantitative analysis of clone sizes scored across multiple passages, we show that the growth of GBM clones in vivo are consistent with a remarkably neutral process resulting from a conserved proliferative hierarchy. In this model, slow-cycling stem-like cells give rise to a more rapidly cycling progenitor population with extensive self-maintenance capacity, that in turn generates short-lived cells that lack proliferative ability. We also identify rare "outlier" clones that deviate from these dynamics, and further show that chemotherapy reproducibly facilitates the selective expansion of pre-existing drug-resistant GBM lineages. Finally, we show that coexisting, functionally distinct GSCs may be targeted separately with epigenetic compounds. These findings show that, independent of an evolving mutational signature, the clonal dynamics of human GBM can be surprisingly robust, with intratumoural heterogeneity derived primarily from the defined stochastic fate behaviour of tumour cells within a conserved proliferative hierarchy.

  Study EGAS00001002424

• Transcriptional regulatory networks of tumor-associated macrophages that drive malignancy in mesenchymal glioblastoma

  Background: Glioblastoma (GBM) is a complex disease with extensive molecular and transcriptional heterogeneity. GBM can be subcategorized into four distinct subtypes; tumors that shift towards the mesenchymal phenotype upon recurrence are generally associated with treatment resistance, unfavorable prognosis, and the infiltration of pro-tumorigenic macrophages. Results: We explore the transcriptional regulatory networks of mesenchymal-associated tumor-associated macrophages (MA-TAMs), which drive the malignant phenotypic state of GBM, and identify macrophage receptor with collagenous structure (MARCO) as the most highly differentially expressed gene. MARCOhigh TAMs induce a phenotypic shift towards mesenchymal cellular state of glioma stem cells, promoting both invasive and proliferative activities, as well as therapeutic resistance to irradiation. MARCOhigh TAMs also significantly accelerate tumor engraftment and growth in vivo. Moreover, both MA-TAM master regulators and their target genes are significantly correlated with poor clinical outcomes, and are often associated with genomic aberrations in neurofibromin 1 (NF1) and phosphoinositide 3-kinases/mammalian target of rapamycin/Akt pathway (PI3K-mTOR-AKT)-related genes. We further demonstrate origination of MA-TAMs from peripheral blood, as well as their potential association with tumor-induced polarization states and immunosuppressive environments. Conclusions: Collectively, our study characterizes the global transcriptional profile of TAMs driving mesenchymal GBM pathogenesis, providing potential therapeutic targets for improving the effectiveness of GBM immunotherapy.

  Study EGAS00001002443

• Coding and small non-coding transcriptional landscape of tuberous sclerosis complex cortical tubers: implications for pathophysiology and treatment

  Tuberous Sclerosis Complex (TSC) is a rare genetic disorder that results from a mutation in the TSC1 or TSC2 genes leading to constitutive activation of the mechanistic target of rapamycin complex 1 (mTORC1). TSC is associated with autism, intellectual disability and severe epilepsy. Cortical tubers are believed to represent the neuropathological substrates of these disabling manifestations in TSC. In the presented study we used high-throughput RNA sequencing in combination with systems-based computational approaches to investigate the complexity of the TSC molecular network. Overall we detected 438 differentially expressed genes and 991 differentially expressed small non-coding RNAs in cortical tubers compared to autopsy control brain tissue. We observed increased expression of genes associated with inflammatory, innate and adaptive immune responses. In contrast, we observed a down-regulation of genes associated with neurogenesis and glutamate receptor signaling. MicroRNAs represented the largest class of over-expressed small non-coding RNA species in tubers. In particular, our analysis revealed that the miR-34 family (including miR-34a, miR-34b and miR-34c) was significantly over-expressed. Functional studies demonstrated the ability of miR-34b to modulate neurite outgrowth in mouse primary hippocampal neuronal cultures. This study provides new insights into the TSC transcriptomic network along with the identification of potential new treatment targets.

  Study EGAS00001002485

• UK10K NEURO ASD TAMPERE

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The Tampere Autism sample set consists of samples from Finnish subjects with ASD (autism spectrum disorders) with IQs over 70 recruited from a clinical centre for the diagnosis and treatment of children with ASD. For further information on this cohort please contact either Terho Lehtimaki (terho.lehtimaki@uta.fi) or Kaija Puura (kaija.puura@pshp.fi).

  Study EGAS00001000115

• Sensitive Monogenic Noninvasive Prenatal Diagnosis by Targeted Haplotyping

  During pregnancy, cell-free DNA (cfDNA) in maternal blood encompasses a small percentage of cell-free fetal DNA (cffDNA), an easilyaccessible source for determination of fetal disease status in risk families through non-invasive procedures. In case of monogenic heritabledisease, background maternal cfDNA prohibits direct observation of the maternally inherited allele. Non-invasive prenatal diagnostics(NIPD) of monogenic diseases therefore relies on parental haplotyping and statistical assessment of inherited alleles from cffDNA,techniques currently unavailable for routine clinical practice. Here, we present monogenic NIPD (MG-NIPD), which requires a bloodsample from both parents, for targeted locus amplification (TLA)-based phasing of heterozygous variants selectively at a gene of interest.Capture probes-based targeted sequencing of cfDNA from the pregnant mother and a tailored statistical analysis enables predicting fetalgene inheritance. MG-NIPD was validated for 18 pregnancies, focusing on CFTR, CYP21A2, and HBB. In all cases we could predict theinherited alleles with >98% confidence, even at relatively early stages (8 weeks) of pregnancy. This prediction and the accuracy ofparental haplotyping was confirmed by sequencing of fetal material obtained by parallel invasive procedures. MG-NIPD is a robustmethod that requires standard instrumentation and can be implemented in any clinic to provide families carrying a severe monogenicdisease with a prenatal diagnostic test based on a simple blood draw.

  Study EGAS00001002622

• Using_genetics_to_identify_cell_types_and_mechanisms_underlying_susceptibility_to_primary_sclerosing_cholangitis

  Primary sclerosing cholangitis (PSC) is a T-cell mediated, chronic inflammatory condition of the biliary tree that is strongly associated with inflammatory bowel disease. Genome-wide association studies have identified 22 non-HLA genetic risk variants associated PSC. Identifying the genes impacted by these variants has proven difficult as the majority lie in non-coding regions of the genome. Knowledge of the genes and biological pathways these non-coding variants are perturbing is vital to understanding the disease biology. One means of assessing the impact of non-coding variants within disease associated loci upon genes is via colocalisation with eQTL. Many eQTL are cell-type specific, requiring the analysis of disease relevant cell types to detect colocalisation. We have collected PSC-relevant T-cell-subtypes from the peripheral blood of PSC patients via fluorescence activated cell sorting in preparation for RNA sequencing and mapping of eQTL. Samples were collected at the Norfolk and Norwich University Hopital, for which local ethical approval has been granted. Lysed cell samples will be transferred to WTSI and DNA/RNA will be extracted from lysed cell samples by T143 before genotyping (DNA) and custom library preparation and sequencing (RNA). This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001002643

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__pancreas

  The Mutographs project aims to advance our understanding of the causes of cancer through studies of mutational signatures. Led by Mike Stratton, together with Paul Brennan, Ludmil Alexandrov, Allan Balmain, David Phillips and Peter Campbell, this large-scale international research endeavour was awarded a Cancer Research UK Grand Challenge. Different patterns of somatic mutation are generated by the different environmental, lifestyle and genetic factors that cause cancer, many of them are still unknown. Within Mutographs, the International Agency for Research on Cancer is coordinating the recruitment of 5000 individuals with cancer (colorectal, renal, pancreatic, oesophageal adenocarcinoma or oesophageal squamous cancers) across 5 continents to explore whether different mutational signatures explain marked variation in incidence. In brief, through an international network of collaborators around the world, biological materials are collected, along with demographic, histological, clinical and questionnaire data. Whole genome sequences of tumour-germline DNA pairs are generated at the Wellcome Trust Sanger Institute (Illumina HiSeqX, 40X and 20X depth respectively). Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

  Study EGAS00001002726

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Oesophageal_squamous_cell_carcinoma

  The Mutographs project aims to advance our understanding of the causes of cancer through studies of mutational signatures. Led by Mike Stratton, together with Paul Brennan, Ludmil Alexandrov, Allan Balmain, David Phillips and Peter Campbell, this large-scale international research endeavour was awarded a Cancer Research UK Grand Challenge. Different patterns of somatic mutation are generated by the different environmental, lifestyle and genetic factors that cause cancer, many of them are still unknown. Within Mutographs, the International Agency for Research on Cancer is coordinating the recruitment of 5000 individuals with cancer (colorectal, renal, pancreatic, oesophageal adenocarcinoma or oesophageal squamous cancers) across 5 continents to explore whether different mutational signatures explain marked variation in incidence. In brief, through an international network of collaborators around the world, biological materials are collected, along with demographic, histological, clinical and questionnaire data. Whole genome sequences of tumour-germline DNA pairs are generated at the Wellcome Trust Sanger Institute (Illumina HiSeqX, 40X and 20X depth respectively). Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

  Study EGAS00001002725