14538 results for "PLOS+ONE+Pennsylvania+German+ACS+language+code+pdc+open+access+study"

in 38.57 milliseconds.

• Intratumoural heterogeneity and immune modulation in lung adenocarcinoma of female smokers and never smokers

  Lung cancer is still the leading cause of cancer death worldwide despite declining smoking prevalence in industrialised countries. Although lung cancer is highly associated with smoking status, a significant proportion of lung cancer cases develop in patients who never smoked, with an observable bias towards female never smokers. A better understanding of lung cancer heterogeneity and immune system involvement during tumour evolution and progression in never smokers is therefore highly warranted. We employed single nucleus transcriptomics of surgical lung adenocarcinoma (LADC) and normal lung tissue samples from patients with or without smoking history. Immune cells as well as fibroblasts and endothelial cells respond to tobacco smoke exposure by inducing a highly inflammatory state in normal lung tissue. In the presence of LADC, we identified differentially expressed transcriptional programmes in macrophages and cancer-associated fibroblasts, providing insight into how the niche favours tumour progression. Within tumours, we distinguished eight subpopulations of neoplastic cells in female smokers and never smokers. Through pseudotemporal ordering, we inferred a trajectory towards two differentiated tumour cell states implicated in cancer progression and invasiveness. A proliferating cell population sustaining tumour growth exhibits differential immune modulating signatures in both patient groups. Our results resolve cellular heterogeneity and immune interactions in LADC, with a special emphasis on female never smokers and implications for the design of therapeutic approaches.

  Study EGAS00001006331

• Bulk-tissue RNA-sequencing paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from post-mortem neuropathologically-confirmed control individuals.

  Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

  Study EGAS00001006380

• An Empirical Approach Leveraging Tumorgrafts to Dissect the Tumor Microenvironment in Renal Cell Carcinoma Identifies Missing Link to Prognostic Inflammatory Factors

  We undertook an empirical approach to dissect the tumor microenvironment by developing a bulk tumor dissection algorithm, DisHet, and leveraging RNA-Seq data of tumorgrafts (patient-derived tumors implanted in mice), in which only the tumor cell component expands. We found that approximately 65% of previously defined immune signature genes are not abundantly expressed in the Renal Cell Carcinoma (RCC) microenvironment, and we identified more than two times as many novel immune/stromal transcripts. By using refined immune/stroma-specific genes and genomics, electronic medical record data and imaging data of 1084 RCC patients, we discovered a highly-inflamed pan-RCC subtype enriched for Treg cells, NK cells, Th1 cells, neutrophils, macrophages, B cells, and CD8+ T cells. This inflamed subtype (IS) is enriched for aggressive RCCs, including BAP1-deficient clear-cell and type 2 papillary tumors. Interestingly, IS is correlated with systemic manifestations of inflammation in patients such as thrombocytosis and anemia, whose pathogenesis were poorly defined and have been predictors of poor prognosis in RCC. Indeed, we discovered that IS is a strong predictor of poor survival. Lastly, our analyses show that tumor cells may drive stromal immune response. Overall, these data provide a missing link between the tumor cells, the tumor microenvironment, and systemic factors.

  Study EGAS00001002786

• PAX5 biallelic genomic alterations define a novel subgroup of B cell precursor acute lymphoblastic leukemia

  Chromosomal rearrangements and specific aneuploidy patterns are initiating events and define subgroups in B cell precursor acute lymphoblastic leukemia (BCP-ALL). Here we analyzed 250 BCP-ALL cases and identified a novel subgroup (‘PAX5-plus’, n=19) by distinct DNA methylation and gene expression profiles. All patients in this subgroup harbored mutations in the B-lineage transcription factor PAX5, with p.P80R as hotspot. Mutations either affected two independent codons, consistent with compound heterozygosity, or suffered LOH predominantly through chromosome 9p aberrations. These biallelic events resulted in disruption of PAX5 transcriptional programs regulating B cell differentiation and tumor suppressor functions. Homozygous CDKN2A/B deletions and RAS activating hotspot mutations were highly enriched as cooperating events in the genomic profile of PAX5-plus ALL. Together this defined a specific pattern of triple alterations, exclusive to the novel subgroup. PAX5-plus ALL was observed in pediatric and adult patients. Although restricted by the limited sample size, a tendency for more favorable clinical outcome was observed with 10 of 12 adult PAX5-plus patients achieving long-term survival. PAX5-plus represents the first BCP-ALL subgroup defined by sequence alterations in contrast to gross chromosomal events and exemplifies how deregulated differentiation (PAX5), impaired cell cycle control (CDKN2A/B) together with sustained proliferative signaling (RAS) cooperatively drives leukemogenesis.

  Study EGAS00001003209

• 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

• 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

• 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

• 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

• Resistance to anti-EGFR therapy in colorectal cancer

  Monoclonal antibodies targeting the Epidermal Growth Factor Receptor (EGFR), such as cetuximab and panitumumab, have evolved to important therapeutic options in metastatic colorectal cancer (CRC). However, almost all patients with clinical response to anti-EGFR therapies show disease progression within a few months and little is known about mechanism and timing of resistance evolution. Here we performed whole genome sequencing of plasma DNA (plasma-Seq) from patients treated with anti-EGFR therapy. We demonstrate that development of resistance to anti-EGFR therapies is frequently associated with focal amplifications of KRAS, MET, and ERBB2. However, we also show that focal KRAS amplifications can be acquired in tumor genomes of patients under cytotoxic chemotherapy. Furthermore, we provide evidence that specific chromosomal polysomies, such as overrepresentations of 12p and 7p, harboring KRAS and EGFR, respectively, determine responsiveness to anti-EGFR therapy. In contrast, employing ultra-sensitive deep sequencing of genes associated with anti-EGFR resistance, such as KRAS, BRAF, PIK3CA, and EGFR, we did not observe the occurrence of novel, acquired mutations. Overall, whole-genome plasma DNA sequencing represents a non-invasive blood-based surrogate measure of changes in tumors. As such, plasma-Seq enables the identification of novel mutant clones and may therefore facilitate early adjustments of therapies that may delay or prevent disease progression.

  Study EGAS00001000582

• Whole-Genome and Epigenomic Landscapes of Etiologically Distinct Subtypes of Cholangiocarcinoma

  Cholangiocarcinoma (CCA) is a hepatobiliary malignancy exhibiting high incidence in countries with endemic liver-fluke infection. We analyzed 489 CCAs from 10 countries, combining whole-genome (71 cases), targeted/exome, copy-number, gene expression, and DNA methylation information. Integrative clustering defined 4 CCA clusters-fluke-positive CCAs (clusters 1/2) are enriched in ERBB2 amplifications and TP53 mutations; conversely, fluke-negative CCAs (clusters 3/4) exhibit high copy-number alterations and PD-1/PD-L2 expression, or epigenetic mutations (IDH1/2, BAP1) and FGFR/PRKA-related gene rearrangements. Whole-genome analysis highlighted FGFR2 3' untranslated region deletion as a mechanism of FGFR2 upregulation. Integration of noncoding promoter mutations with protein-DNA binding profiles demonstrates pervasive modulation of H3K27me3-associated sites in CCA. Clusters 1 and 4 exhibit distinct DNA hypermethylation patterns targeting either CpG islands or shores-mutation signature and subclonality analysis suggests that these reflect different mutational pathways. Our results exemplify how genetics, epigenetics, and environmental carcinogens can interplay across different geographies to generate distinct molecular subtypes of cancer.SIGNIFICANCE: Integrated whole-genome and epigenomic analysis of CCA on an international scale identifies new CCA driver genes, noncoding promoter mutations, and structural variants. CCA molecular landscapes differ radically by etiology, underscoring how distinct cancer subtypes in the same organ may arise through different extrinsic and intrinsic carcinogenic processes.

  Study EGAS00001001653

• mFAST-SeqS

  Recent progress in the analysis of cell-free DNA fragments (cell-free circulating tumor DNA, ctDNA) now allows monitoring of tumor genomes by non-invasive means. However, previous studies with plasma DNA from patients with cancer demonstrated highly variable allele frequencies of ctDNA. The comprehensive analysis of tumor genomes is greatly facilitated when plasma DNA has increased amounts of ctDNA. Therefore, a fast and cost-effective pre-screening method to identify such plasma samples without previous knowledge about alterations in the respective tumor genome could assist in the selection of samples suitable for further extensive qualitative analysis. To address this, we adapted the recently described FAST-SeqS method, which was originally established as a simple and effective, non-invasive screening method for fetal aneuploidy from maternal blood. We show that our modified FAST-SeqS method (mFAST-SeqS) can be used as a pre-screening tool for an estimation of the ctDNA percentage. Using a combined evaluation of genome-wide and chromosome-arm specific z-scores from dilution series with cell line DNA and by comparisons of plasma-Seq profiles with data from mFAST-SeqS, we established a detection limit of 10% or more of mutant alleles. Plasma samples with an mFAST-SeqS z-score above 5 showed highly concordant results compared to copy number profiles obtained from our previously described plasma-Seq approach.

  Study EGAS00001001133

• Homozygous loss-of-function variants in European cosmopolitan and isolate populations

  Homozygous loss of function (HLOF) variants provide a valuable window on gene function in humans, as well as an inventory of the human genes that are not essential for survival and reproduction. All humans carry at least a few HLOF variants, but the exact number of inactivated genes that can be tolerated is currently unknown-as are the phenotypic effects of losing function for most human genes. Here, we make use of 1432 whole exome sequences from five European populations to expand the catalogue of known human HLOF mutations; after stringent filtering of variants in our dataset, we identify a total of 173 HLOF mutations, 76 (44%) of which have not been observed previously. We find that population isolates are particularly well suited to surveys of novel HLOF genes because individuals in such populations carry extensive runs of homozygosity, which we show are enriched for novel, rare HLOF variants. Further, we make use of extensive phenotypic data to show that most HLOFs, ascertained in population-based samples, appear to have little detectable effect on the phenotype. On the contrary, we document several genes directly implicated in disease that seem to tolerate HLOF variants. Overall HLOF genes are enriched for olfactory receptor function and are expressed in testes more often than expected, consistent with reduced purifying selection and incipient pseudogenisation.

  Study EGAS00001001606

• Identifying genetic consequences of Epstein-Barr Virus transformation by comparing an individual’s genomic DNA with that of its lymphoblastoid cell line.

  A detailed analysis of whole genomes can be now achieved with next generation sequencing. Epstein Barr Virus (EBV) transformation is a widely used strategy in clinical research to obtain an unlimited source of a subject’s DNA. Although the mechanism of transformation and immortalization by EBV is relatively well known at the transcriptional and proteomic level, the genetic consequences of EBV transformation are less well understood. A detailed analysis of the genetic alterations introduced by EBV transformation is highly relevant, as it will inform on the usefulness and limitations of this approach. We used whole genome sequencing to assess the genomic signature of a low-passage lymphoblastoid cell line (LCL). Specifically, we sequenced the full genome (40X) of an individual using DNA purified from fresh whole blood as well as DNA from his LCL. A total of 217.33 Gb of sequence were generated from the cell line and 238.95 Gb from the normal genomic DNA. We determined with high confidence that 99.2% of the genomes were identical, with no reproducible changes in structural variation (chromosomal rearrangements and copy number variations) or insertion/deletion polymorphisms (indels). Our results suggest that, at this level of resolution, the LCL is genetically indistinguishable from its genomic counterpart and therefore their use in clinical research is not likely to introduce a significant bias.

  Study EGAS00001000323

• Mutated H3 Histones Drive Human Pre-Leukemic Hematopoietic Stem Cell Expansion And Promote Leukemic Aggressiveness

  Leukemogenesis is a stepwise progression from mutated, pre-neoplastic hematopoietic stem cells (HSCs) to full-blown leukemia. Our ability to prevent or treat acute myeloid leukemia (AML) is limited by our incomplete understanding of the epigenetic disruption that is central to this process, including improper histone methylation. Comprehensive analysis of 16 histone H3 genes in 434 primary AML samples identified Q69H, A26P, R2Q, R8H and K27M/I mutations (1.6%), and a higher incidence in secondary AML (s-AML) (9%). We establish that these mutations are important early events in leukemogenesis. They occur in pre-leukemic HSCs, increase the frequency of functional human HSCs, and alter differentiation. The mutations are present in the major leukemic clones in primary samples, and the mutant histones amplify leukemic aggressiveness with increased proliferation, expansion of leukemic progenitor and blast cells, and superior competitiveness in vivo. These effects are dependent on the specific mutation. Genome-wide analysis of K27 mutants revealed increased expression of genes involved in erythrocyte and myeloid differentiation with a corresponding decrease in histone H3 K27 tri-methylation and increase in K27 acetylation. The functional impact of histone mutations is independent of RUNX1 mutations, although they can co-occur. These data establish the involvement of H3 mutations as initial drivers of pre-cancerous stem cell expansion and leukemogenesis.

  Study EGAS00001003288

• Extracellular transcriptome in seminal plasma for non-invasive prostate cancer diagnosis

  A diagnostic non-invasive biomarker test for prostate cancer at an early stage, with high sensitivity and specificity, would improve diagnostic decision making. Extracellular RNAs present in seminal plasma might contain biomarker potential for the accurate detection of clinically significant prostate cancer. So far, the extracellular messenger RNA (mRNA) profile of seminal plasma has not been interrogated for its biomarker potential in the context of prostate cancer. Here, we investigate the mRNA transcriptome in seminal plasma samples obtained from prostate cancer patients (n=25), patients with benign prostate hyperplasia (n=26) and individuals without prostatic disease (n=6). Seminal plasma harbors a complex mRNA repertoire that reflects prostate as its tissue of origin. The endogenous RNA content is higher in the prostate cancer samples compared to the control samples. Prostate cancer antigen 3 (PCA3), a long non-coding RNA with prostate cancer-specific overexpression, and ATP-binding cassette transporter 1 (ABCA1), known to be involved in the prostate cancer pathogenesis, were more abundant in the prostate cancer group. In addition, twelve high confidence fusion transcripts could be detected in prostate cancer samples, including the bona-fide prostate cancer fusion transcript TMPRSS2-ERG. Our findings provide proof-of-principle that the extracellular transcriptome of seminal plasma can reveal information of an underlying prostate cancer.

  Study EGAS00001005277

• Bulk-tissue RNA-sequencing of anterior cingulate cortex samples derived from Lewy body disease patients

  Parkinson’s disease (PD), Parkinson’s disease with dementia (PDD) and dementia with Lewy bodies (DLB) are three clinically, genetically and neuropathologically overlapping neurodegenerative diseases collectively known as the Lewy body diseases (LBDs). A variety of molecular mechanisms have been implicated in PD pathogenesis, but the mechanisms underlying PDD and DLB remain largely unknown, a knowledge gap that presents an impediment to the discovery of disease-modifying therapies. Transcriptomic profiling can contribute to addressing this gap, but remains limited in the LBDs. Here, we applied paired bulk-tissue and single-nucleus RNA-sequencing to anterior cingulate cortex samples derived from 28 individuals, including healthy controls, PD, PDD and DLB cases (n = 7 per group), to transcriptomically profile the LBDs. Using this approach, we (i) found transcriptional alterations in multiple cell types across the LBDs; (ii) discovered evidence for widespread dysregulation of RNA splicing, particularly in PDD and DLB; (iii) identified potential splicing factors, with links to other dementia-related neurodegenerative diseases, coordinating this dysregulation; and (iv) identified transcriptomic commonalities and distinctions between the LBDs that inform understanding of the relationships between these three clinical disorders. Together, these findings have important implications for the design of RNA-targeted therapies for these diseases and highlight a potential molecular “window” of therapeutic opportunity between the initial onset of PD and subsequent development of Lewy body dementia.

  Study EGAS00001005305

• 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

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Oesophageal_adenocarcinoma

  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 (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 EGAS00001003702

• Small_molecule_inhibitors_in_melanoma___Kenski___Kong___WES

  Small-molecule inhibitors targeting the most commonly activated pathway in melanoma, MAPK pathway (either alone or in combination) are already given to melanoma patients for few years, and initially reduce tumour burden dramatically, eventually melanomas become resistant and tumours progress while on treatment. Resistance to this treatment occurs by acquisition of additional mutations or other alterations that affect the mitogen-activated protein kinase (MAPK) pathway by either direct or indirect signalling. Many resistance mechanisms somehow lead to reactivation of extracellular signal-regulated kinase (ERK), thereby restoring signalling of the oncogenic BRAF/MEK/ERK pathway. In addition, PI3K pathway activation contributes to resistance to BRAF inhibition. Less frequent but equally important to the phenomenon of targeted drug resistance is the observation that B15–20% of BRAF mutant melanoma patients fail to respond to BRAF inhibition already early on treatment, owing to intrinsic resistance. These patients have little therapeutic options, unless immunotherapy can be given. To better understand the resistance mechanisms in MAPK inhibitor-treated melanoma patients and melanoma biology, our lab generated a big panel of MAPK inhibitor resistant melanoma cell lines by continuous drug exposure. The understanding of the genetic landscape and gene expression as well as cross resistance to other treatment regimens, and other aspects of melanoma biology such as phenotype switch, will allow us to better exploit new therapeutic strategies for melanoma patients.

  Study EGAS00001002863

• Transcriptional and DNA binding profiling of CEBPE in REH acute lymphoblastic leukaemia cells using shRNA RNA-Seq and ChIP-Seq

  ChIP-Seq - CEBPE - REH. The ETV6/RUNX1 translocated acute lymphoblastic leukaemia cell line REH was used to perform ChIP-Seq using a CEBPE antibody. Cells were fixed in 1% formaldehyde for 10mins, prior to preparation of chromatin using Active Motif Express ChIP-IT. 2ug of antibody (anti CEBPE Atlas Antibodies HPA002928)was added to 25ug of chromatin O/N at 4C with rotation. Duplicate reactions were pooled and purified. 10ng of ChIPd and input DNA used for Illumina NGS preparation (NEBNext ChIP-Seq Library kit; New England Biolabs), CEBPE and Input DNA ChIP samples were sequenced on a MiSeq using 150bp Kit v3 paired end and a HiSeq 2500 using 2x101 version 4 paired end (Illumina) respectively. Reactions performed in duplicate. shCEBPE RNA-Seq - REH. REH cells were lentivirally transduced with a pTRIPZ shRNA vector for transcriptional profiling of CEBPE. Two controls (empty and non-targeting) and two CEBPE shRNAs (V3THS_150517(A13), V3THS_404312(G3) Dharmacon, GE) were transduced into REH cells. Cells were treated with 1ug/ml doxycyclin for 144hrs and total RNA purified. Libraries were prepared using NEBNext Ultra II Directional RNA Library Prep Kit and sequenced on an Illuimna HiSeq 2500 using 2x101 version 4 paired end chemistry. 3 biological replicates of each samples were prepared.

  Study EGAS00001002877

• scRNA-seq data of human nuclei collected from the temporal cortex of 17 individuals.

  Transcriptome analysis of single cells has tremendously increased our understanding of cellular heterogeneity in complex tissues. However, whole-cell transcriptomics cannot be applied to frozen tissue due to poor cell recovery after freezing. Human tissue is usually limited to tissue banks containing frozen material, and since single nuclei can be readily isolated from frozen tissues, human cell nuclei can substitute whole-cell analysis to extract single cell transcriptomes from human cells. Despite the importance of human tissue analysis, single-nucleus RNA-sequencing protocols are less developed than those for whole cells. SMARTSeq, either in form of custom-made SMARTSeq2 or as part of the Fluidigm C1 platform, provides a high resolution view of the transcriptome, has proven effective for whole-cell analysis and has recently been adapted to nuclei. Here, we show that the standard SMARTSeq2 protocol applied to single nuclei biases the amplification of single-nucleus transcriptomes due to unspecific capturing of transcripts, which can comprise the majority of reads in single-nucleus RNA-sequencing libraries. We demonstrate that modification of the SMARTSeq2 protocol with respect to template switching and cDNA amplification conditions generates sequencing libraries of higher complexity, allows for detection of more genes and reduces biases of the amplified transcriptome. In addition, modification of SMARTSeq2 chemistry allows for gene detection at shallower sequencing depth, thus increasing efficiency and decreasing cost.

  Study EGAS00001002882

• Glioblastoma epigenome profiling identifies SOX10 as a master regulator of molecular tumour subtype

  Glioblastoma in adult patients are classified into four subtypes, IDH, MES, RTK I, and RTK II, based on DNA methylation and RNA expression data. Tumour subtype transitions commonly occur during the treatment of glioblastoma patients, and transitions to the mesenchymal (MES) subtype have been associated with therapy resistance and adverse prognosis. Here, we generate DNA methylome and histone modification data of glioblastoma primary tumours and show that glioblastoma subtypes differ by their enhancer landscapes. Using Core Regulatory Circuitry analysis of H3K27ac data, and independent analysis of RNA-derived gene regulatory networks, we identified 38 subtype master regulators whose cell type-specific activities we mapped in single-cell RNA sequencing data. These analyses identified the chromatin modifier SRY-Box 10 (SOX10) as a master regulator in tumours of the RTK I subtype. In functional studies we demonstrate that SOX10 loss causes a transition of RTK I tumour cells to a mesenchymal cellular state by altering the accessibility of enhancer-rich chromatin regions and their occupancy by Bromodomain Containing 4 (BRD4), a dynamic constituent of super-enhancers. Treatment with the BRD4 inhibitor JQ1 blocked this transition. These data demonstrate the fundamental role of enhancer chromatin remodelling and master regulator activity in the establishment, maintenance and plasticity of glioblastoma cellular states. We also make our data publically browsable at https://dkfz-b060.github.io/gb_browse/

  Study EGAS00001003953

• Deep multi-region whole-genome sequencing reveals heterogeneity and gene-by-environment interactions in treatment-naive, metastatic lung cancer

  Our understanding of genomic heterogeneity in lung cancer is largely based on the analysis of early-stage surgical specimens. Here we used endoscopic sampling of paired primary and intrathoracic metastatic tumors from 11 lung cancer patients to map genomic heterogeneity inoperable lung cancer with deep whole-genome sequencing. Intra-patient heterogeneity in driver or targetable mutations was predominantly in the form of copy number gain. Private mutation signatures, including patterns consistent with defects in homologous recombination, were highly variable both within and between patients. Irrespective of histotype, we observed a smaller than expected number of private mutations, suggesting that ancestral clones accumulated large mutation burdens immediately prior to metastasis. Single-region whole-genome sequencing of from 20 patients showed that tumors in ever-smokers with the strongest tobacco signatures were associated with germline variants in genes implicated in the repair of cigarette-induced DNA damage. Our results suggest that lung cancer precursors in ever-smokers accumulate large numbers of mutations prior to the formation of frank malignancy followed by rapid metastatic spread. In advanced lung cancer, germline variants in DNA repair genes may interact with the airway environment to influence the pattern of founder mutations, whereas similar interactions with the tumor microenvironment may play a role in the acquisition of mutations following metastasis.

  Study EGAS00001003830

• Sequence Data from the phase 2 PrE0505 trial of Durvalumab with First Line Platinum-Pemetrexed for Unresectable Pleural Mesothelioma

  Mesothelioma is a rare and fatal cancer with limited therapeutic options until the recent approval of combination immune checkpoint blockade. We report the results of the phase 2 PrE0505 trial (NCT02899195) of the anti-PD-L1 antibody, durvalumab, plus platinum-based and pemetrexed chemotherapy for patients with previously untreated unresectable pleural mesothelioma. The primary endpoint was overall survival compared to historic control with cisplatin and pemetrexed chemotherapy. The combination of durvalumab with chemotherapy met the pre-specified primary endpoint reaching a median survival of 20.4 months vs. 12.1 months with historic control. Treatment emergent adverse events were consistent with known side effects of chemotherapy and all adverse events due to immunotherapy were grade ≤2. Integrated genomic and immune cell repertoire analyses revealed that a higher immunogenic mutation burden coupled with a more diverse T cell repertoire were linked with favorable clinical outcome. Structural genome-wide analyses demonstrated a higher degree of genomic instability in responding tumors of epithelioid histology. Patients with germline alterations in cancer predisposing genes, especially those involved in DNA repair, were more likely to attain long term survival. Our findings indicate that concurrent durvalumab with platinum-based chemotherapy has promising clinical activity and that responses are driven by the complex genomic background of malignant pleural mesothelioma.

  Study EGAS00001005426

• Comprehensive Genomic Characterization of Gene Therapy-Induced T-cell Acute Lymphoblastic Leukemia (H007)

  Six out of ten patients in a hematopoietic stem cell gene therapy (GT) trial for Wiskott-Aldrich syndrome (WAS) developed T-cell acute lymphoblastic leukemia (T-ALL) associated with gamma-retroviral vector integration at the LMO2 locus. We hypothesized that detailed molecular analysis of T-ALL driven by temporally defined first deregulation of LMO2 might offer valuable insights into the pathogenesis of GT-associated and sporadic T-ALL (HIPO project 007). In all six patients whole-genome and transcriptome sequencing demonstrated that in addition to oncogenic chromosomal translocations involving the T-cell receptor loci, retrovirally induced leukemias faithfully recapitulated secondary genomic events observed in sporadic T-ALL. Multiple genetic events culminated in the development of a rapidly proliferating leukemic clone and onset of clinical manifestation at 1.3 to 5 years after infusion of the gene corrected transplant. These genetic alterations affected NOTCH1 and PI3K signaling, cell cycle regulation as wells as T-cell transcription factors LEF1 and TCF1. Two of the six cases displayed a phenotype reminiscent of early immature T-ALLs, and lacked additional T-ALL typical alterations. Our analysis shows that insertional mutagenesis involving the LMO2 locus triggers T-ALL development due to diverse genetic alterations in large parts, but not exclusively, in known T-ALL driver genes, demonstrating the molecular diversity of GT-associated T-ALL.

  Study EGAS00001003870

• 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

• Coding and non-coding drivers of mantle cell lymphoma identified through exome and genome sequencing

  Mantle cell lymphoma (MCL) is an uncommon B-cell non-Hodgkin lymphoma (NHL) that is incurable with standard therapies. The genetic drivers of this cancer have not been firmly established and the features that contribute to differences in clinical course remain limited. To extend our understanding of the biological pathways involved in this malignancy, we performed a large-scale genomic analysis of MCL using data from 51 exomes and 34 genomes alongside previously published exome cohorts. To confirm our findings, we re-sequenced the genes identified in the exome cohort in 191 MCL tumors, each having clinical follow-up data. We confirmed the prognostic association of TP53 and NOTCH1 mutations. Our sequencing revealed novel recurrent non-coding mutations surrounding a single exon of the HNRNPH1 gene. In RNA-seq data from 103 of these cases, MCL tumors with these mutations had a distinct imbalance of HNRNPH1 isoforms. This altered splicing of HNRNPH1 was associated with inferior outcomes in MCL and showed a significant increase in protein expression by immunohistochemistry. We describe a functional role for these recurrent non-coding mutations in disrupting an auto-regulatory feedback mechanism, thereby deregulating HNRNPH1 protein expression. Taken together, these data strongly implicate a role for aberrant regulation of mRNA processing in MCL pathobiology.

  Study EGAS00001004289

• Single cell transcriptomics of hESC-derived midbrain dopaminergic neurons generated by a new human development-based protocol

  Recent advances in stem cell technologies provide new opportunities for regenerative medicine and Parkinson’s disease modelling. However, the precise composition and quality of human embryonic stem cell (hESC)-derived preparations at the single cell level, compared to endogenous brain standards, remains to be determined. In addition, multiple developmental pathways important for midbrain dopaminergic (mDA) neuron development remain to be implemented in hESCs. Here we use mouse developmental knowledge and human single cell RNA-sequencing (scRNA-seq) data to generate a new development-based protocol to differentiate hESCs into mDA neurons. We found that Laminin 511, dual WNT activation with CHIR99021 and WNT5A, together with FGF8b, improved different aspects of midbrain patterning. Moreover, mDA neurogenesis and differentiation were enhanced by activation of liver X receptors and inhibition of fibroblast growth factor signaling. Analysis of hESC-derived midbrain cultures by scRNA-seq revealed multiple hESC-derived clusters, which resembled those in the endogenous human ventral midbrain. Moreover, high quality mDA neurons were detected and found to progressively acquire functionality in vitro. Thus, we hereby provide a hESC differentiation protocol that recapitulates several aspects of endogenous mDA neuron development, and a strategy to determine cell composition and quality of hESC-derived preparations that can be useful for disease modeling and cell therapy.

  Study EGAS00001006313

• Hyperhaploid Multiple Myeloma

  Hyperhaploid multiple myeloma is a rare numerical aberration group defined by a range of 24-34 chromosomes which is associated with a poor prognosis. Hyperhaploid samples were sequenced by either exome or targeted panel sequencing. Samples had a median of 13 monosomies (range 12-14), which in general were those not associated with trisomies in hyperdiploid samples. The chromosomes traditionally trisomic in hyperdiploid myeloma were disomic in hyperhaploid myeloma with retention of heterodisomy. We examined the hyperhaploid samples for frequently mutated genes and found that 9/9 (100%) hyperhaploid samples had a mutation in TP53, exceeding the overall rate of mutation in newly diagnosed patients (5.5%), indicating an oncogenic dependency in this group. All samples with TP53 mutation also had monosomy of chromosome 17, indicating bi-allelic inactivation of TP53. As such, this high risk group is part of Double-Hit myeloma. DNA was obtained from either CD138+ cells from the bone marrow of multiple myeloma patients (tumor) or from stem cell harvests from the same patient (control). 100 ng of DNA was fragmented, end-repaired, and adapters ligated using the HyperPlus kit (KAPA Biosystems). After PCR amplification the libraries were hybridized with probes against either the entire exome (MedExome, Nimblegen) or a targeted panel of 140 genes using SeqCap reagents (Nimblegen). Hybridized libraries underwent further amplification before being sequenced on a NextSeq500 (Illumina) using 75 bp paired end reads.

  Study EGAS00001003203

• Patient-derived tumor organoids for personalized medicine in a rare case of hepatocellular carcinoma with neuroendocrine differentiation.

  We report the case of a 74-year-old man with a very rare subtype of hepatocellular carcinoma with neuroendocrine differentiation (HCC-NED). The patient presented with two independent tumors, a gastrointestinal stromal tumor in the stomach and a hepatocellular carcinoma in the liver. Both tumors were surgically removed in curative intent. Histopathological work-up of the liver tumor revealed poorly differentiated hepatocellular carcinoma (Edmondson-Steiner grade IV) with diffuse expression of neuroendocrine markers synaptophysin (SYP) and chromogranin (CHGA). Three months after resection, multifocal recurrence of the HCC with neuroendocrine differentiation (HCC-NED) was observed. In the meantime, tumor organoids have been generated from the resected HCC-NED and extensively characterized. Sensitivity to a number of drugs approved for the treatment of HCC or neuroendocrine carcinomas was tested in vitro. Based on their in vitro efficacy, etoposide and carboplatin were used as first line palliative combination treatment. Because genomic analysis revealed a NTRK1-mutation (kinase domain) and tumor organoids were sensitive to entrectinib, a pan-TRK inhibitor, the patient received entrectinib as second line therapy. After only two weeks, treatment had to be discontinued due to deterioration of the patient’s general condition. In conclusion, we demonstrate for the first time that preclinical drug testing using organoids is feasible in selected HCC cases.

  Study EGAS00001005887

• scRNA-seq of HGSC tumor and ascites

  Immune regulatory metabolites are key features of the tumor microenvironment (TME), yet with a few notable exceptions, their identities remain largely unknown. We uncovered the immune regulatory metabolic states and metabolomes of sorted tumor and stromal, CD4+, and CD8+ cells from the tumor and ascites of patients with high-grade serous ovarian cancer (HGSC) using high-dimensional flow cytometry and metabolomics supplemented with single cell RNA sequencing. Flow cytometry revealed that tumor cells show a consistently greater uptake of glucose than T cells, but similar mitochondrial activity. Cells within the ascites and tumor had pervasive metabolite differences, with a striking enrichment in 1-methylnicotinamide (MNA) in T cells infiltrating the tumor compared to ascites. Despite the elevated levels of MNA in T cells, the expression of nicotinamide N-methyltransferase, the gene encoding the enzyme that catalyses the transfer of a methyl group from S-adenosylmethionine to nicotinamide, was restricted to fibroblasts and tumor cells. Treatment of T cells with MNA resulted in an increase in T cell-mediated secretion of the tumor promoting cytokine tumor necrosis factor alpha. Thus, the TME-derived metabolite MNA contributes to an alternative and non-cell autonomous mechanism of immune modulation of T cells in HGSC. Collectively, uncovering the tumor-T cell metabolome may reveal metabolic vulnerabilities that can be exploited using T cell-based immunotherapies to treat human cancer.

  Study EGAS00001004829

• Genomics to select patients with metastatic breast cancer for targeted therapy (microarray_oncoscan)

  Cancer progression is driven in part by genomic alterations located in oncogenes or tumor suppressor genes. The genomic characterization of cancers has shown a large interpatient heterogeneity regarding the driver alterations, leading to the concept that generating a genomic profiling by multigene sequencing in patients with cancer could allow selecting effective targeted therapies. While this concept has been broadly implemented in daily practice, there is no evidence that such approach improves patient outcome, and how to optimally select the therapy to administer. A genomic profiling using next generation sequencing and copy number analyses was performed 1462 patients with Her2-non overexpressing metastatic breast cancer included in SAFIR02 Breast trial. 238 of these patients were randomized in two trials between a targeted therapy matched to genomic alteration and a maintenance chemotherapy. The trial shows that targeted therapies matched to genomics improves progression free survival (PFS) when genomic alterations are classified level I/II according to ESMO Scale for Clinical Actionability of Molecular Targets (ESCAT) (adjusted HR: 0.41, 90% CI: 0.27-0.61, p<0.001), but not when alterations are classified beyond level II (unadjusted HR: 1.15, 95% CI: 0.76-1.75). This trial provides evidence that the treatment decision led by genomics should be driven by a framework of target actionability in patients with mBC.

  Study EGAS00001005586

• Genomics to select patients with metastatic breast cancer for targeted therapy (microarray_agilent)

  Cancer progression is driven in part by genomic alterations located in oncogenes or tumor suppressor genes. The genomic characterization of cancers has shown a large interpatient heterogeneity regarding the driver alterations, leading to the concept that generating a genomic profiling by multigene sequencing in patients with cancer could allow selecting effective targeted therapies. While this concept has been broadly implemented in daily practice, there is no evidence that such approach improves patient outcome, and how to optimally select the therapy to administer. A genomic profiling using next generation sequencing and copy number analyses was performed 1462 patients with Her2-non overexpressing metastatic breast cancer included in SAFIR02 Breast trial. 238 of these patients were randomized in two trials between a targeted therapy matched to genomic alteration and a maintenance chemotherapy. The trial shows that targeted therapies matched to genomics improves progression free survival (PFS) when genomic alterations are classified level I/II according to ESMO Scale for Clinical Actionability of Molecular Targets (ESCAT) (adjusted HR: 0.41, 90% CI: 0.27-0.61, p<0.001), but not when alterations are classified beyond level II (unadjusted HR: 1.15, 95% CI: 0.76-1.75). This trial provides evidence that the treatment decision led by genomics should be driven by a framework of target actionability in patients with mBC.

  Study EGAS00001005587

• Targeted de-methylation of the FOXP3-TSDR

  CD4+ regulatory T cells (Tregs) are key mediators of immunological tolerance and promising effector cells for immuno-suppressive adoptive cellular therapy to fight autoimmunity and chronic inflammation. Their functional stability is critical for their clinical utility and has been correlated to the demethylated state of the TSDR/CNS2 enhancer element in the Treg lineage transcription factor FOXP3. However, proof for a causal contribution of the TSDR de-methylation to FOXP3 stability and Treg induction is so far lacking. We here established a powerful transient-transfection CRISPR-Cas9-based epigenetic-editing method for the selective de-methylation of the TSDR within the endogenous chromatin environment of a living cell. The induced de-methylated state was stable over weeks in clonal T cell proliferation cultures even after expression of the editing complex had ceased. Epigenetic editing of the TSDR resulted in FOXP3 expression, even in its physiological isoform distribution, proving a causal role for the de-methylated TSDR in FOXP3 regulation. However, successful FOXP3 induction was not associated with a switch towards a functional Treg phenotype, in contrast to what has been reported from FOXP3 overexpression approaches. Thus, TSDR de-methylation is required, but not sufficient for a stable Treg phenotype induction. Therefore, targeted demethylation of the TSDR may be a critical addition to published in vitro Treg induction protocols which so far lack FOXP3 stability.

  Study EGAS00001004867

• 5- FU treated organoids

  5-Fluorouracil (5-FU) is a chemotherapeutic drug component that is commonly used for the treatment of solid cancers. The anticancer properties of 5-FU have been attributed to interference with nucleotide synthesis and through direct incorporation into the DNA as being a pyrimidine analog. As both mechanisms of action may have a mutational impact on surviving tumor cells, we performed three independent analyses to characterize the genomic consequences of 5-FU treatment: i) in vitro treatment of intestinal organoids with 5-FU followed by whole genome sequencing, ii) genome-wide mutation analyses in tumor samples from 5-FU treated patients with breast or colorectal cancer, and iii) analysis of paired biopsies from patients treated with 5-FU between biopsies. Our results demonstrate that both in vitro and in vivo 5-FU causes a mutational pattern that is dominated by T>G substitutions in a CTT context and with strong resemblance to COSMIC signature 17. Interestingly, this signature is also found in non 5-FU treated patients, most prominently in esophageal tumors, indicating that distinct endogenous and exogenous triggers can converge into highly similar mutational signatures. Furthermore, our results suggest that 5-FU may have adverse mutagenic effects on healthy cells and contributes to genetic variability in surviving cancer cells thereby contributing to tumor evolution.

  Study EGAS00001003592

• Genomic insights into the pathogenesis of Epstein-Barr virus-associated diffuse large B-cell lymphoma by whole-genome and targeted amplicon-based sequencing

  Epstein–Barr virus (EBV)-associated diffuse large B-cell lymphoma not otherwise specified (DLBCL NOS) constitutes a distinct clinicopathological entity in the current World Health Organization (WHO) classification, yet its genomic features remain sparsely characterized. Combining whole-genome sequencing (WGS), targeted amplicon sequencing (tNGS) and fluorescence in-situ hybridization (FISH) of tumours from 47 cases of EBV+ DLBCL we delineate the genomic landscape of this rare disease. Interrogating WGS data by sequential gene-set enrichment and network-diffusion we describe an accumulation of mutations affecting the nuclear factor κB (NFκB) pathway (CSNK2A2, CARD10), IL-6/JAK/STAT (SOCS1/3, STAT3) and WNT signalling (FRAT1, SFRP5) alongside aberrations in immunological processes like interferon-response. Further, WGS uncovered structural aberrations including 6q deletions (5/8 patients), which were subsequently validated by FISH (14/32 cases). Integrated WGS and tNGS analysis emphasized a profile distinct from EBV-negative DLBCL enriched for mutations in ARID1A (45%), KMT2A/D (32/30%), ANKRD11 (32%), NOTCH2 (32%). Expanding on previous reports, we identified recurrent alterations in CCR6 (15%), DAPK1 (15%), TNFRSF21 (13%), CCR7 (11%) and YY1 (6%). Our comprehensive description of the mutational landscape in EBV+ DLBCL substantiates the evidence of its pathobiological independence. Further, these observations aid in the refinement of the molecular taxonomy for aggressive lymphomas and guides future therapeutic strategies.

  Study EGAS00001004941

• Genetic insights into the biological mechanisms governing human ovarian ageing

  Reproductive longevity is critical for fertility and impacts healthy ageing in women, yet insights into the underlying biological mechanisms and treatments to preserve it are limited. Here, we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in ~200,000 women of European ancestry. These common alleles influence clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations. Identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increase fertility and extend reproductive life in mice. Causal inference analyses using the identified genetic variants indicates that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases risks of hormone-sensitive cancers. These findings provide insight into the mechanisms governing ovarian ageing, when they act across the life-course, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

  Study EGAS00001004947

• Control of Focal Adhesion Kinase Activation by RUNX1-regulated miRNAs in high-risk AML

  We recently described a 16-gene expression signature for improved risk stratification of acute myeloid leukemia (AML) patients called the AML Prognostic Score (APS). A subset of APS high-risk AML patients showed increased levels of focal adhesion kinase (FAK), encoded by the Protein Tyrosine Kinase 2 (PTK2) gene, which was correlated with RUNX1 mutations. RUNX1 mutant cells are more sensitive to PTK2 inhibitors. As we were not able to detect RUNX1 binding sites in the PTK2 promoter, we hypothesized that RUNX1 might regulate micro(mi)RNAs that repress PTK2, such that loss-of-function RUNX1 mutations would result in reduced miRNA expression and derepression of PTK2. Examination of paired RNA-seq and miRNA-seq data from 301 AML cases revealed two miRNAs that positively correlated with RUNX1 expression, contained RUNX1 binding sites in their promoters and were predicted to target PTK2. We show that the hsa-let7a-2-3p and hsa-miR-135a-5p promoters are regulated by RUNX1, and that PTK2 is a direct target of both miRNAs. Even in the absence of RUNX1 mutations, hsa-let7a-2-3p and hsa-miR-135a-5p regulate PTK2 expression, and reduced expression of these two miRNAs sensitizes AML cells to PTK2 inhibition. These data explain how RUNX1 regulates PTK2, and identify potential miRNA biomarkers for targeting AML with PTK2 inhibitors.

  Study EGAS00001006491

• Uncovering tumor intrinsic and extrinsic factors that regulate hepatocellular carcinoma growth using patient derived xenograft assays

  Several major factors limit our understanding of hepatocellular carcinoma (HCC). First, human HCCs are infrequently biopsied for diagnosis and thus are not often biologically interrogated. Second, HCC initiation and progression are strongly influenced by the cirrhotic microenvironment, and the exact contributions of intrinsic and extrinsic tumor factors are unclear. A powerful approach to examine the personalized biology of liver cancers and the influence of host tissues is with patient derived xenograft (PDX) models. In Asia, HCCs from hepatitis B virus patients have been efficiently converted into PDXs, but few parallel efforts from the west have been reported. In a large-scale analysis, we implanted 93 HCCs and 8 cholangiocarcinomas (CCAs) to systematically analyze host factors and to define an optimized platform for PDX development from both surgical and biopsy samples. NSG mice that had undergone partial hepatectomy (PHx) represented the best combination of engraftability, growth, and passageability, but overall rates were low and indicative of a unique intrinsic biology for HCCs in the US. PDX models preserved the histology and genetic features of parental tumors, and ultimately, 9 new models were usable for pre-clinical studies. Intriguingly, HCC PDXs were differentially sensitive to regorafenib and sorafenib and CCA PDXs were also highly sensitive to regorafenib. PDX models functionalize early and advanced stage HCCs and revealed unique biological features of liver cancers from the US.

  Study EGAS00001004020

• Cell-free DNA fragmentation patterns and personalized sequencing reveal circulating tumor DNA in urine and plasma of glioma patients

  Glioma-derived cell-free DNA (cfDNA) is challenging to detect using liquid biopsy as levels in body fluids are low. We determined the glioma-derived DNA fractions in tumor biopsies, in cerebrospinal fluid (CSF), plasma and urine samples, using deep sequencing of personalized capture panels. By sequencing cfDNA across thousands of mutations identified individually in each patient’s matched tumor we detected tumor-derived DNA in plasma (10/12) and urine samples (8/11). The median tumor fraction was 6.4x10-3 in CSF, 3.1x10-5 in plasma and 4.7x10-5 in urine. We identified a shift in the size distribution for mutant cfDNA fragments in these body fluids. Next, we analyzed cfDNA fragment sizes with paired-end shallow whole genome sequencing (WGS) in urine samples from 35 patients with gliomas, 8 individuals with non-malignant brain disorders, and 26 controls (n=69 individuals, 96 samples). cfDNA in urine of glioma patients was significantly more fragmented compared to urine from patients with non-malignant brain disorders (t-test, p=1.7x10-2) and compared to urine of controls (t-test, p=5.2x10-9). The proportion of DNA fragments <60 bp was higher in glioma patients urine and could be used for classification (AUC=0.93). Machine learning models integrating fragment lengths could identify urine samples from glioma patients (AUC=0.97 in cross-validation).

  Study EGAS00001004355

• Germline whole genome sequencing of patients with Li-Fraumeni syndrome

  Five hundred ng to 1 ug of genomic DNA was submitted to The Centre for Applied Genomics (TCAG) at The Hospital for Sick Children for genomic library preparation and whole genome sequencing. DNA samples were quantified using the Qubit High Sensitivity Assay and purity was assessed using the Nanodrop OD 260/280 ratio. Approximately 500-700 ng of DNA was used as input material for library preparation using the Illumina TruSeq PCR-free DNA Library Prep Kit following the manufacturer’s recommended protocol. In brief, DNA was fragmented to 400 bp on average using sonication on a Covaris LE220 instrument. Fragmented DNA was end-repaired, A-tailed and indexed TruSeq Illumina adapters with overhang-T were added to the DNA. Libraries were then validated on a Bioanalyzer DNA High Sensitivity chip to check for size and absence of primer dimers and quantified by qPCR using Kapa Library Quantification Illumina/ABI Prism Kit protocol (KAPA Biosystems). Validated libraries were pooled in equimolar quantities and paired-end sequenced on an Illumina HiSeq X platform following Illumina’s recommended protocol to generate paired-end reads of 150-bases in length and an average depth of 40X. FASTQ files were aligned to the hg19 reference genome using BWA-mem v0.78. PCR duplicates were marked with Picard MarkDuplicates v1.1.08, and base recalibration and realignment was performed using GATK v2.8.1.

  Study EGAS00001007061

• Machine learning guided signal enrichment for ultrasensitive plasma tumor burden monitoring

  In solid tumor oncology, circulating tumor DNA (ctDNA) is poised to transform care through accurate assessment of minimal residual disease (MRD) and therapeutic response monitoring. To overcome the sparsity of ctDNA fragments in low tumor fraction (TF) settings and increase MRD sensitivity, we previously leveraged genome-wide mutational integration through plasma whole genome sequencing (WGS). We now introduce MRD-EDGE, a composite machine learning-guided WGS ctDNA single nucleotide variant (SNV) and copy number variant (CNV) detection platform designed to increase signal enrichment. MRD-EDGESNV uses deep learning and a ctDNA-specific feature space to increase SNV signal-to-noise enrichment in WGS by 300X compared to previous WGS error suppression. MRD-EDGECNV also reduces the degree of aneuploidy needed for ultrasensitive CNV detection through WGS from 1 Gb to 200 Mb, thereby expanding its applicability to a wider range of solid tumors. We harness the improved performance to identify MRD following surgery in multiple cancer types, track changes in tumor burden in response to neoadjuvant immunotherapy in non-small cell lung cancer (NSCLC) and demonstrate ctDNA shedding in precancerous colorectal adenomas. Finally, the radical signal-to-noise enrichment in MRD-EDGESNV enables de novo mutation calling in melanoma and small cell lung cancer (SCLC) without matched tumor, yielding clinically informative TF monitoring for patients on immune checkpoint inhibition (ICI).

  Study EGAS00001007545

• RRBS data (cfSort study) from noncancer tissue DNA

  The RRBS libraries of the genomic DNA from the 521 tissue samples were constructed following the standard RRBS protocol. 100-200 ng of intact genomic DNA in the volume of 21.5 µl was used as input material. Restriction digestion was done with 2.5 µl 10xCutSmart buffer and 1 µl MspI (NEB) for 18 h at 37 oC and 20 min at 65 oC. 0.5 µl 10xCutSmart buffer, 0.3 µl dACGTP mixture (100 mM dATP, 10 mM dCTP, 10 mM dGTP), 1 µl Klenow (exo-, 5U/µl, NEB) and 2.6 µl RT-PCR water, 0.6 µl 50 mM DTT (ThermoFisher) was added to the mixture for end repair and A-overhang addition with the program 30 oC for 20 min, 37 oC for 1 h and 75 oC for 20 min. Adapter ligation was then performed with 1 µl 10xThermoFisher HC T4 ligase buffer, 0.4 µl 100 mM ATP (ThermoFisher), 0.2 µl 50 mM DTT, 1 µl ThermoFisher HC T4 DNA ligase (30 Weiss Unit/µl), 30 ng home-made duplex UMI adapter with all the cytosines methylated (protocol adopted from Kennedy et al.) at 16 oC for 20 h and 65 oC for 20 min. Bisulfite conversion of the adapter-ligated product was carried out with QIAGEN EpiTect plus DNA bisulfite kit following their protocol for two rounds of conversion. The converted product was purified with Qiagen MinElute spin column and eluted with 20 µl RT-PCR water. PCR amplification was done using the NEBNext Multiplex Oligos for Illumina (2.5 µl of universal and index primer each) and 25 µl KAPA HiFi HotStart Uracil+ ReadyMix (Roche) with the following cycling conditions: 98 oC for 45 s, 9 cycles of 98 oC for 15 s, 60 oC for 30 s and 72 oC for 30 s, followed by a final extension at 72 oC for 5 min. The PCR product was purified with 1x AmpureXP beads and eluted with 30 µl EB buffer. DNA concentration was measured by Qubit 1xdsDNA HS assay. 5% TBE-UREA PAGE and bioanalyzer assay was performed as quality control on each library before sequencing.

  Dataset EGAD00001010880

• The patterns and dynamics of genomic instability in metastatic pancreatic cancer

  Pancreatic cancer is an aggressive malignancy with a five-year mortality of 97–98%, usually due to widespread metastatic disease. Previous studies indicate that this disease has a complex genomic landscape, with frequent copy number changes and point mutations, but genomic rearrangements have not been characterized in detail. Despite the clinical importance of metastasis, there remain fundamental questions about the clonal structures of metastatic tumours, including phylogenetic relationships among metastases, the scale of ongoing parallel evolution in metastatic and primary sites, and how the tumour disseminates. Here we harness advances in DNA sequencing to annotate genomic rearrangements in 13 patients with pancreatic cancer and explore clonal relationships among metastases. We find that pancreatic cancer acquires rearrangements indicative of telomere dysfunction and abnormal cell-cycle control, namely dysregulated G1-to-S-phase transition with intact G2–M checkpoint. These initiate amplification of cancer genes and occur predominantly in early cancer development rather than the later stages of the disease. Genomic instability frequently persists after cancer dissemination, resulting in ongoing, parallel and even convergent evolution among different metastases. We find evidence that there is genetic heterogeneity among metastasis-initiating cells, that seeding metastasis may require driver mutations beyond those required for primary tumours, and that phylogenetic trees across metastases show organ-specific branches. These data attest to the richness of genetic variation in cancer, brought about by the tandem forces of genomic instability and evolutionary selection.

  Study EGAS00000000064

• Characterizing the Role of the Immune Microenvironment in Multiple Myeloma Progression at a Single Cell Level

  The plasma cell disease multiple myeloma (MM) is highly dependent on its bone marrow microenvironment (BMME). Early alterations within the bone marrow microenvironment that contribute to the progression of multiple myeloma (MM) from its precursor stages could hold the key to identifying novel therapeutic approaches, yet the intrinsic variability in cellular populations between patients together with differences in sample processing and analysis methods have made it difficult to identify consistent changes between datasets. Here, we used single-cell RNA sequencing of BM cells from precursor stages, monoclonal gammopathy of unknown significance (MGUS), smoldering MM (SMM) and newly diagnosed MM and analyzed our data in combination with a previously published dataset that used a similar patient population and sample processing. Despite vast interpatient heterogeneity, some alterations were seen consistently in both datasets. We identified changes in immune cell populations as disease progressed that were characterized by a substantial decrease in memory and naïve CD4 T cells and an increase in CD8+ effector T cells and T-regulatory cells. These alterations were further accompanied by an enrichment of non-clonal memory B cells and an increase of CD14 and CD16 monocytes in MM compared to its precursor stages. These results provide crucial information of the immune changes associated with the progression to clinical MM and can help to develop immune based strategies for patient stratification and early therapeutic intervention. 

  Study phs002756

• Integrated Epigenetic Maps of Human Embryonic, Extraembryonic and Adult Cells

  Working cooperatively with other Mapping Centers and the Data Coordination Center (EDACC) funded by this Roadmap mechanism we have comprehensively mapped the epigenomes of select human cells with significant relevance to complex human disease. Our effort has been focused on cells relevant to human health and complex disease including cells from the blood, brain, breast, skin, male germ cells, extraembryonic tissues and trophoblast, and U.S. Government-approved lines of human embryonic stem cells. We have incorporated high quality, homogeneous cells from males and females, and two predominant racial groups, and biological replicates of each cell type. The resulting comprehensive maps include up to 6 histone modifications selected for their opposing roles in regulating active and inactive chromatin (H3K4me1, H3K4me3, H3K9me3, H3K9ac, H3K27ac, H3K27me3, H3K36me3), DNA methylation, miRNA and gene expression, and for select cases whole genome shotgun sequence. Our cell and data production pipeline incorporates verification and data validation with independent methods, and operates under a model motivated by increased data production and decrease cost. Our group in conjunction with the other REMC teams, the EDACC, ENCODE, Epigenetics of Human Health and Disease centers and the NIH Roadmap program are developing methods, tools and reference epigenome maps for the research community that will make the promise of epigenetics in understand and treating human complex disease a reality. Our reference epigenomes will enable new disciplines including human population epigenetics, comparative epigenomics, neuroepigenetics, and therapeutic epigenetics for tissue regeneration and reversal of disease.

  Study phs000791

• NIMH (National Institute of Mental Health) De Novo Mutation Identification in Taiwanese Schizophrenia Trios

  The substantial reproductive impact of schizophrenia, for which affected individuals have fewer than half as many offspring as unaffected individuals do, implies that mutations of largest effect will frequently be de novo mutations. Ascertaining exome sequence variation in father-mother-offspring trios allows such mutations to be identified and distinguished from the far-larger amount of rare variation that is inherited by each individual. The pursuit of this approach in a large, well-powered cohort of trios can also provide lessons that inform the development of such gene discovery strategies more generally in human genetics. Schizophrenia trios from the Taiwanese population are being collected by Dr. Ming Tsuang (PI, UC San Diego, California) and investigators in Taiwan (PI, Dr. Hai Gwo Hwu; both funded by NIMH grant 1R01MH085560; Expanding Rapid Ascertainment Networks of Schizophrenia Families in Taiwan). A total of 3800 trios are anticipated to be collected by May 2013. This represents a highly homogenous national sample from the same ancestral population. DNA samples will be obtained from the NIMH Repository, Rutgers University Cell and DNA Repository (described below) and stored at the Broad Institute. Genetic and data analyses will be performed at the Broad Institute. We propose to sequence the whole exome of trios by hybrid capture and Illumina next generation sequencing and perform targeted genotyping and validation of variants (SNPs, indels and CNVs) using several molecular methods, to include emulsion-based PCR and Sanger sequencing.

  Study phs001196

• Genetics of Disorders Affecting Tooth Structure, Number, Morphology and Eruption

  The long-term objective of this research project is to discover the complete spectrum of genes that contribute to the etiologies of genetic disorders of dentition. Inherited tooth defects can be isolated (the only symptom) or syndromic (having additional symptoms), and can be caused by mutations in many different genes. Clinically, it can be difficult to distinguish among the diverse forms. Making an accurate diagnosis would be greatly improved if reliable tests could identify the specific genetic defect that causes the disease in each case. Current technology can identify the numerous sequence variations throughout a patient's exome, but the disease-causing variation among them cannot be identified unless it is already on the list of genes that potentially cause the disease. A diagnosis based upon the genetic defect that causes the disorder distinguishes between isolated and syndromic tooth defects. A precise diagnosis permits formulation of a more accurate prognosis, which informs treatment decisions. The benefits of accomplishing our aims are to 1) enable genetic testing to achieve a precise and accurate diagnosis, and thereby optimize treatment, and 2) improve our understanding of normal and pathological tooth formation, leading to future improvements in treatment and eventually cures. Our short-term objectives are to 1) recruit kindreds with inherited disorders affecting their dentition, characterize their dental phenotype (and non-dental phenotypes when applicable) and the distribution of the disorder throughout the family, and 2) determine the genetic etiology of their disorders.

  Study phs001491

• RNA-Sequencing of B-Lymphoblastic Leukemia with Glucocorticoids and PI3K Delta Inhibition

  Mononuclear cells were isolated from bone marrow or peripheral blood of 7 pediatric subjects with B-lymphoblastic leukemia at the time of diagnosis. To determine the effect of PI3K delta inhibition on glucocorticoid-induced gene regulation, B-lymphoblastic leukemia cells were treated with glucocorticoids +/- the PI3K delta inhibitor, idelalisib, for 24 hours prior to RNA extraction and sequencing. Biological replicates were performed for 4 of the 7 subjects; 2 of these subjects had specimens demonstrating an additive response to the combination of glucocorticoids and idelalisib in viability testing in vitro, and 2 subjects' specimens demonstrated a synergistic response. All 7 subjects were analyzed together and the 4 specimens with biological replicates were also evaluated in an additional analysis. An average of 40 million reads per specimen and treatment condition were obtained. We compared gene expression in each treatment condition to vehicle-treated cells and expression in glucocorticoid or idelalisib only treated cells to combination treated cells. Using a high dose of dexamethasone (25 or 50 nM) and an equivalent dose of prednisolone, we found that the addition of idelalisib (500 nM) enhanced prednisolone-induced gene regulation but not dexamethasone-induced gene regulation. Surprisingly, the specimens with additive response to combination therapy demonstrated enhancement of prednisolone-induced gene regulation with the addition of idelalisib, while the specimens with synergistic responses did not. Raw sequencing data (FASTQ) and count tables will be available through dbGaP.

  Study phs003085