14305 results for ""Pennsylvania+German"+ACS+study+'Socius'+open+access"

in 31.97 milliseconds.

• The Haplotype-Resolved Genome and Epigenome of the Aneuploid HeLa Cancer Cell Line

  The genomic architecture of HeLa remains largely unexplored beyond its karyotype, partly because like many cancers, its extensive aneuploidy renders such analyses challenging. We carried out haplotype-resolved whole-genome sequencing of the HeLa CCL-2 strain, examined point- and indel-mutation variations, mapped copy-number variations and loss of heterozygosity regions, and phased variants across full chromosome arms. We also investigated variation and copy-number profiles for HeLa S3 and eight additional strains. We find that HeLa is relatively stable in terms of point variation, with few new mutations accumulating after early passaging. Haplotype resolution facilitated reconstruction of an amplified, highly rearranged region of chromosome 8q24.21 at which integration of the human papilloma virus type 18 (HPV-18) genome occurred and that is likely to be the event that initiated tumorigenesis. We combined these maps with RNA-seq and ENCODE Project data sets to phase the HeLa epigenome. This revealed strong, haplotype specific activation of the proto-oncogene MYC by the integrated HPV-18 genome approximately 500 kilobases upstream, and enabled global analyses of the relationship between gene dosage and expression. These data provide an extensively phased, high-quality reference genome for past and future experiments relying on HeLa, and demonstrate the value of haplotype resolution for characterizing cancer genomes and epigenomes. Reprinted from Adey et. al. 2013, Nature, with permission from Nature Publishing Group.

  Study phs000642

• A thymic ILC1-like progenitor with differentiation potential towards KIR+NKG2A- NK cells

  KIR+NKG2A- Natural Killer (NK) cells have the unique ability to detect downregulation of single HLA-I allotypes, frequently occurring in malignantly-transformed and virus-infected cells in order to escape T cell recognition. We have recently shown that circulating Innate Lymphoid Cells 1 (cILC1) have the potential to generate such KIR+NKG2A- NK cells but their developmental origin was unknown. Here we demonstrate that the development of cILC1 is thymus-dependent and identify a putative progenitor of cILC1s in the thymus, called thyILC1. Single-cell RNAseq analysis revealed a close relationship of thyILC1s to CD34+ double-negative (DN) thymocytes, but in contrast to DNs, thyILC1s had lost T cell differentiation potential. However, like cILC1s they could be efficiently differentiated into KIR+NKG2A- NK cells. thyILC1s strongly expressed RAG1/2, PTCRA (preTCR, and NOTCH1/3, but lacked expression of TCR chains, suggesting that they failed to pass the -selection checkpoint of thymocyte development. Finally, analyses of patients with FOXN1het deficiency, showing thymic hypoplasia due to lack of mature thymic epithelial cells, exhibited a profound deficiency of cILC1s but not cILC2s and cILC3s, demonstrating thymus dependency of cILC1s. Together, the data suggest that thyILC1s are the source of a novel thymus-dependent NK cell differentiation pathway that promotes generation of KIR+NKG2A- NK cells.

  Study EGAS50000000760

• A thymic ILC1-like progenitor with differentiation potential towards KIR+NKG2A- NK cells - sc

  KIR+NKG2A- Natural Killer (NK) cells have the unique ability to detect downregulation of single HLA-I allotypes, frequently occurring in malignantly-transformed and virus-infected cells in order to escape T cell recognition. We have recently shown that circulating Innate Lymphoid Cells 1 (cILC1) have the potential to generate such KIR+NKG2A- NK cells but their developmental origin was unknown. Here we demonstrate that the development of cILC1 is thymus-dependent and identify a putative progenitor of cILC1s in the thymus, called thyILC1. Single-cell RNAseq analysis revealed a close relationship of thyILC1s to CD34+ double-negative (DN) thymocytes, but in contrast to DNs, thyILC1s had lost T cell differentiation potential. However, like cILC1s they could be efficiently differentiated into KIR+NKG2A- NK cells. thyILC1s strongly expressed RAG1/2, PTCRA (preTCR, and NOTCH1/3, but lacked expression of TCR chains, suggesting that they failed to pass the -selection checkpoint of thymocyte development. Finally, analyses of patients with FOXN1het deficiency, showing thymic hypoplasia due to lack of mature thymic epithelial cells, exhibited a profound deficiency of cILC1s but not cILC2s and cILC3s, demonstrating thymus dependency of cILC1s. Together, the data suggest that thyILC1s are the source of a novel thymus-dependent NK cell differentiation pathway that promotes generation of KIR+NKG2A- NK cells.

  Study EGAS50000000790

• Chromatin remodeling enhances MAP3K8 expression in HAM: a key pathogenesis for therapeutic intervention

  Human T-cell leukemia virus type 1 (HTLV-1)-associated myelopathy (HAM) is a debilitating neuroinflammatory disease with no available effective treatments. A hallmark of HAM is the transformation of HTLV-1-infected cells into T helper type 1 (Th1)-like cells, characterized by excessive interferon (IFN)-gamma production that drives chronic inflammation. However, the molecular mechanisms fuel this aberrant Th1-like transformation and sustained inflammation remain poorly understood. We hypothesized that HAM-characteristic chromatin remodeling plays a pivotal role in the overexpression of key genes driving inflammatory pathogenesis. Using transcriptomic analysis, chromatin accessibility profiling, and biomarker evaluations across HTLV-1-related diseases, we identify MAP3K8 as a key gene that defines the unique inflammatory profile of HAM. MAP3K8 overexpression promotes Th1-like differentiation and constitutively activates the MEK-ERK signaling pathway. Furthermore, we elucidate the mechanism by which HTLV-1 Tax, Fosl2, and c-Jun collaboratively induce HAM-characteristic chromatin remodeling at the enhancer region of the MAP3K8 locus. Crucially, we demonstrate that mitogen-activated protein kinase kinase (MEK) inhibitors effectively suppress the MAP3K8-MEK signaling cascade and significantly mitigated inflammatory pathogenesis in an ex vivo culture assay. Our findings provide critical insights into the virus-host interactions underpinning HAM and propose the MAP3K8-MEK-ERK axis as a promising therapeutic target for this challenging condition.

  Study JGAS000835

• MRCA and MRCE SNP genotypes

  Gene expression levels can be an important link DNA between variation and phenotypic manifestations. Our previous map of global gene expression based on approximately 400K SNPs and 50K transcripts in 400 sib pairs from the MRCA family panel has been widely used to interpret the results of GWAS. Here, we more than double the size of our initial dataset with expression data on 550 additional individuals from the MRCE family panel using the Illumina whole genome expression array. We have used new statistical methods for dimension reduction to account for non-genetic effects in estimates of expression levels and we have also included SNPs imputed from the 1000 Genomes Project. Our methods reduced false discovery rates and increased the number of eQTLs mapped either locally or at a distance (i.e. in cis or trans) from 1,534 in the MRCA dataset to 4,452 (with <5% FDR). Imputation of 1000 Genomes SNPs further increased the number of eQTLs to 7,302. Using the same methods and imputed SNPs in the newly acquired MRCE dataset we identified eQTLs for 9000 genes. The combined results identify strong local and distant effects for transcripts from 14,177 genes. Our eQTL database (http://www.hsph.harvard.edu/faculty/liming-liang/software/eqtl/) based on these results is freely available to help define the function of disease-associated variants.

  Study EGAS00000000137

• 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

• 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

• Inherited chromosomally-integrated human herpesvirus 6A/B (HHV-6A/B) genome sequences in the Japanese population

  The deposited sequences represent inherited chromosomally-integrated human herpesvirus 6 (iciHHV-6) genomes identified from 10 BioBank Japan (BBJ) subjects based on the analysis of the whole-genome sequencing (WGS) data. In brief, we re-aligned unmapped reads from the WGS to HHV-6A genome (Genbank accession number: KJ123690.1) using BWA-MEM algorithm. Ten out of 3,256 subjects were identified with high depth of coverage of HHV-6. We further distinguished subjects with integrated HHV-6A subtypes from those with integrated HHV-6B based on phylogenetic analysis of concatenated sequences of three viral genes (U27/U43/U83, which show high divergence between viral species). We reconstructed the HHV-6 viral genome of 10 subjects. First, reads that mapped to a reference HHV-6A genome (KJ123690.1) were further aligned against the integrated HHV-6A genome derived from a Japanese individual NA18999 (GenBank Accession number: KY316047.1) using BWA-MEM. Based on the alignment, variant calling was performed using freebayes (version: v1.2.0-2-g29c4002) with parameters ploidy = 1 and min-alternate-fraction = 0.8 {Garrison, 2012}. We generated subject-specific iciHHV-6 viral sequences by applying the resulting variants to the KY316047.1 reference genome using the FastaAlternateReferenceMaker function in the Genome Analysis Toolkit (GATK) v3.7. Please note that HHV-6B-specific variants are unlikely to be sensitively reconstructed using this method.

  Study JGAS000240

• Treg cell subset-specific gene expression patterns in human head and neck cancer

  Backgrand: The occurrence of head and neck cancer has been increased all over the world. Although the prognosis has been improved by using several treatment modalities including immunotherapy, it is still unsatisfying. It is still unclear why the immunotherapy is effective only for partial amount of patients. To resolve this problem and to make a new treatment strategy based on human immune system, further investigation of antitumor immunity is required. Regulatory T (Treg) cells play critical roles in human immune responses including antitumor immunity. Therefore, the research to elucidate Treg function is strongly urged to clarify the mechanisms of antitumor immunity. Methods: We performed RNA-seuencing analysis of regulatory T cells infiltrating head and neck cancer. We devided into four populations, DP, CD25+, CTLA4+, DN and analyzed gene expressions. Results: The gene expression patterns of Treg-related genes were different between CD25+ cells and CD25- cells. Furthermore, DP cells had unique gene expression pattern compared to other populations. Conclusions: We found Treg cells infiltrating head and neck cancer had high expression levels of Treg-related genes and DP cells had unique gene expression pattern that may become a target of new treatment strategy. Regulatory T (Treg) cells play critical roles in human immune responses including antitumor immunity. Therefore, the research to elucidate Treg function is strongly urged to clarify the mechanisms of antitumor immunity.

  Study JGAS000135

• 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

• 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

• 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

• 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

• 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&amp;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

• 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

• 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

• 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

• 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

• Analysis_of_genomic_integrity_of_disease_corrected_human_induced_pluripotent_stem_cells_by_exome_sequencing

  Induced pluripotent stem (iPS) cells hold great promise for autologous cell transplantation. In order to apply this therapy to monogenic disorder, disease-causing mutations must be corrected prior to transplantation. We generated iPS cells from patients with alpha-1 antitrypsin deficiency, which is caused by homozygous non-synonymous mutation in the SERPINA1 gene. In vitro differentiation of patient-derived iPS cells to hepatocytes showed recapitulation of the disease phenotype. We corrected the mutations by gene targeting using zinc finger nuclease (ZFN) and piggyBac transposon technologies. Our mutation correction is highly accurate. There is no exogenous sequence left behind at the targeted site. The corrected iPS cell-derived hepatocytes were indistinguishable from healthy hepatocyte, suggesting that the corrected gene was fully functional. The genomic integrity was also analyzed by CGH and SNP array, showing no genomic alterations. These are important points to ensure clinical safety as well as functional restoration of the corrected gene. One remaining question is if these corrected iPS cells are clinically safe at the single base-pair level. Recent publications suggested that cancer-related mutations are significantly accumulated during iPS cell propagation. We aim to assess the clinical safety of the corrected iPS cells by using exome sequencing. Outcome from this project would bring iPS cell technology one step closer to clinical applications.

  Study EGAS00001000055

• 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

• 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

• 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

• Real-time analysis of the cancer genome and fragmentome from plasma and urine short and long cell-free DNA using Nanopore sequencing

  Cell-free DNA (cfDNA) can be isolated from blood and/or urine of cancer patients and analyzed with sequencing. Unfortunately, most conventional short-read sequencing methods are technically challenging, labor intensive and time consuming, requiring several days but more typically weeks to obtain interpretable data which are limited by a bias for short cfDNA fragments. Here, we demonstrate that with Oxford Nanopore Technologies sequencing  we can achieve economical and ultra-fast delivery of clinical data from liquid biopsies.   Our ‘ITSFASTR’ approach is able to deliver copy number aberrations, and cfDNA fragmentation profiles in less than 24 hours from sample collection. The tumor-derived cfDNA fraction calculated from lung cancer patient plasma and urine from bladder cancer patients was highly correlated (R=0.98) to the tumor fraction calculated from conventional short-read sequencing of the same samples. cfDNA size profile and fragmentation patterns in plasma and urine exhibited the typical cfDNA features yet with a significantly higher proportion of informative fragments that exceed 300bp, exhibiting similar tumor fraction than shorter size fragments. Notably, comprehensive fragment-end composition and nucleosome profiling near transcription start sites can be determined from the same data.  We propose that ITSFASTR is the first point-of-care solution for obtaining genomic and fragmentomic results from liquid biopsies.

  Study EGAS00001006591

• Genome wide variation in the Angolan Namib desert reveals unique Pre-Bantu ancestry

  The populations of the Angolan Namib Desert have been largely neglected in previous surveys of the genomic landscape of southern Africa. Although at present, the Namib is culturally dominated by Southwest Bantu-speaking cattle-herders, the region exhibits an extraordinary ethnographic diversity which includes an array of semi-nomadic peoples whose subsistence strategies fall outside the traditional division between foraging and food production and can thus be referred to as “peripatetic”. Among these small-scale populations are the last speakers of the Kwadi branch of the Khoe-Kwadi language family associated with the introduction of pastoralism into southern Africa (Kwepe), as well as a range of groups whose origins remain enigmatic (Kwisi, Twa and Tjimba). Using genome-wide data from 208 individuals belonging to nine ethnically diverse groups from the Angolan Namib and adjacent areas (Kwepe, Kwisi, Twa, Tjimba, !Xun, Kuvale, Himba, Nyaneka, Ovimbundu) in combination with published data from other regions of Africa, we reconstruct in detail the histories of contact emerging from pre-historic migrations to southern Africa and show that peripatetic groups from southwestern Angola stand out for exhibiting elevated levels of an unique, regionally-specific and highly divergent Pre-Bantu ancestry. These findings highlight the importance of the Namib for understanding the deep genetic structure of Africa.

  Study EGAS00001007011

• 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

• Double mutant DNMT3A AML: a unique subtype

  Mutation of DNMT3A, encoding a de novo methyltransferase essential for cytosine methylation, is a common early event in clonal hematopoiesis (CH) and adult acute myeloid leukemia (AML). Spontaneous deamination of methylated cytosines damages DNA, which is repaired by the base excision repair (BER) enzymes MBD4 and TDG. Congenital MBD4-deficiency has been linked to early-onset CH and AML, and is marked by exceedingly high levels of DNA damage and mutation of DNMT3A. Strikingly, wildtype (WT) DNMT3A binds TDG, thereby potentiating its repair activity. Since TDG is the only remaining BER enzyme in MBD4-deficient AML patients capable of repairing methylation damage, we investigated whether mutant DNMT3A negatively affects the repair function of TDG. We found that, whereas WT DNMT3A stimulates TDG function, mutant DNMT3A impairs TDG-mediated repair of DNA damage in vitro. In light of this finding and to extrapolate our observations to the broader AML patient population, we investigate here the genetic profiles and survival outcomes of AML patients with single (SM) versus double mutant (DM) DNMT3A. DM DNMT3A AML patients show a characteristic driver mutation landscape and reduced overall survival when compared to SM DNMT3A AML patients. Importantly, whole-genome sequencing showed a trend for increased DNA damage in primary DM DNMT3A AML samples, especially when DNMT3A mutations are located at the DNMT3A-TDG interaction interface.

  Study EGAS00001007966

• deep-learning-powered tissue deconvolution for cfDNA

  Plasma cell-free DNA (cfDNA) is a noninvasive biomarker for cell death of all organs. Deciphering the tissue origin of cfDNA can reveal abnormal cell death because of diseases, which has great clinical potential in disease detection and monitoring. Despite the great promise, the sensitive and accurate quantification of tissue-derived cfDNA remains challenging to existing methods due to the limited characterization of tissue methylation and the reliance on unsupervised methods. To fully exploit the clinical potential of tissue-derived cfDNA, here we present one of the largest comprehensive and high-resolution methylation atlas based on 521 noncancer tissue samples spanning 29 major types of human tissues. We systematically identified fragment-level tissue-specific methylation patterns and extensively validated them in orthogonal datasets. Based on the rich tissue methylation atlas, we develop the first supervised tissue deconvolution approach, a deep-learning-powered model, cfSort, for sensitive and accurate tissue deconvolution in cfDNA. On the benchmarking data, cfSort showed superior sensitivity and accuracy compared to the existing methods. We further demonstrated the clinical utilities of cfSort with two potential applications: aiding disease diagnosis and monitoring treatment side effects. The tissue-derived cfDNA fraction estimated from cfSort reflected the clinical outcomes of the patients. In summary, the tissue methylation atlas and cfSort enhanced the performance of tissue deconvolution in cfDNA, thus facilitating cfDNA-based disease detection and longitudinal treatment monitoring.

  Study EGAS00001007213

• 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

• Insular Celtic population structure and genomic footprints of migration

  Previous studies of the genetic landscape of Ireland have suggested homogeneity, with population substructure undetectable using single-marker methods. Here we have harnessed the haplotype-based method fineSTRUCTURE in an Irish genome-wide SNP dataset, identifying 23 discrete genetic clusters which segregate with geographical provenance. Cluster diversity is pronounced in the west of Ireland but reduced in the east where older structure has been eroded by historical migrations. Accordingly, when populations from the neighbouring island of Britain are included, a west-east cline of Celtic-British ancestry is revealed along with a particularly striking correlation between haplotypes and geography across both islands. A strong relationship is revealed between subsets of Northern Irish and Scottish populations, where discordant genetic and geographic affinities reflect major migrations in recent centuries. Additionally, Irish genetic proximity of all Scottish samples likely reflects older strata of communication across the narrowest inter-island crossing. Using GLOBETROTTER we detected Irish admixture signals from Britain and Europe and estimated dates for events consistent with the historical migrations of the Norse-Vikings, the Anglo-Normans and the British Plantations. The influence of the former is greater than previously estimated from Y chromosome haplotypes. In all, we paint a new picture of the genetic landscape of Ireland, revealing structure which should be considered in the design of studies examining rare genetic variation and its association with traits.

  Study EGAS00001002769

• 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

• 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

• 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

• Single cell sequencing reveals the origin and the order of mutation acquisition in T-cell acute lymphoblastic leukemia.

  Next-generation sequencing has provided a detailed overview of the various genomic lesions implicated in the pathogenesis of T-cell acute lymphoblastic leukemia (T-ALL). Typically, 10 to 20 genomic lesions are found in T-ALL cells at diagnosis. However, it is currently unclear in which order these mutations are acquired and in which progenitor cells this is initiated. To address these questions, we used targeted single-cell sequencing of total bone marrow cells and CD34+CD38- multipotent progenitor cells for 4 T-ALL cases. Hierarchical clustering detected a dominant leukemia clone at diagnosis, accompanied by a few subclones harboring only a fraction of the mutations. We developed a graph-based algorithm to determine the order of mutation acquisition. Two of the four patients had an early event in a known oncogene (MED12, STAT5B) amongst various pre-leukemic events. Intermediate events included loss of 9p21 (CDKN2A/B) and acquisition of fusion genes, while NOTCH1 mutations were typically late events. Analysis of CD34+CD38- cells and myeloid progenitors revealed that in half of the cases somatic mutations were detectable in multipotent progenitor cells. We demonstrate that targeted single cell sequencing can elucidate the order of mutation acquisition in T-ALL and that T-ALL development can start in a multipotent progenitor cell.

  Study EGAS00001002830

• A genomic approach towards an understanding of clonal evolution and disease progression in multiple myeloma

  Multiple Myeloma (MM) is a largely incurable haematological malignancy defined by the clonal proliferation of malignant plasma cells within the bone marrow. Clonal heterogeneity has recently been established as a feature in MM, however, the subclonal evolution associated with disease progression has not been described. We used whole exome sequencing to analyse 10 paired patient samples, providing new insights into the progression from Monoclonal Gammopathy of Undetermined Significance (MGUS) and Smouldering MM (SMM), to symptomatic MM. We confirm that clonal heterogeneity is a common feature at diagnosis and that the driving events involved in disease progression are more complex than previously reported. While we observe some previously identified known “drivers” of MM, we find that the driving events involved in progression are complex and not limited to the known SNVs or CNVs. The RAS/MAPK pathway was found to be the most frequently deregulated pathway, with recurrent mutations in KRAS and NRAS observed in patients at both MGUS/SMM and MM stages. We reveal that MM evolution is mainly characterised by the phenomenon of clonal stability, where the subclonal plasma cell populations identified at MM progression are already present in the asymptomatic MGUS/SMM stages. These subclonal populations could be amenable to therapeutic intervention to arrest transformation to MM.

  Study EGAS00001002850

• 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

• 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

• 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 &lt;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

• 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 &gt;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

• 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

• 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

• 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

• 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 EGAS00001002642

• Chromothripsis in human breast cancer (HIPO K26K/H017/A017)

  Chromothripsis is a form of genome instability, by which a presumably single catastrophic event generates extensive genomic rearrangements of one or a few chromosome(s). Widely assumed to be an early event in tumor development, this phenomenon plays a prominent role in tumor onset. We analyzed chromothripsis in 252 human breast cancers from two patient cohorts (149 metastatic breast cancers, 63 untreated primary tumors, 29 local relapses, 11 longitudinal pairs) using whole-genome and whole-exome sequencing. We showed that chromothripsis affects a substantial proportion of human breast cancers, with a prevalence over 60% in a cohort of metastatic cases and 25% in a cohort comprising predominantly luminal breast cancers (cohorts from HIPO K26K and H017 and A017). In the vast majority of cases, multiple chromosomes per tumor are affected, with most chromothriptic events on chromosomes 11 and 17 including, among other significantly altered drivers, CCND1, ERBB2, CDK12 and BRCA1. Importantly, chromothripsis generates recurrent fusions that drive tumor development. Chromothripsis-related rearrangements are linked with univocal mutational signatures, with clusters of point mutations due to kataegis in close proximity to the genomic breakpoints, and with the activation of specific signaling pathways. Analysis of the temporal order of events in tumors with and without chromothripsis as well as longitudinal analysis of chromothriptic patterns in tumor pairs revealed important insights on the role of chromothriptic chromosomes in tumor evolution.

  Study EGAS00001004662