23014 results for "ega"

in 253.78 milliseconds.

• APOBEC3 mutagenesis drives therapy resistance in breast cancer

  Acquired genetic alterations commonly drive resistance to endocrine and targeted therapies in metastatic breast cancer, however the underlying processes engendering these diverse alterations are largely uncharacterized. To identify the mutational processes operant in breast cancer and their impact on clinical outcomes, we utilized a well-annotated cohort of 3,880 patient samples with paired tumor-normal sequencing data. The mutational signatures associated with apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) enzymes were highly prevalent and enriched in post-treatment compared to treatment-naïve hormone receptor-positive (HR+) cancers. APOBEC3 mutational signatures were independently associated with shorter progression-free survival on antiestrogen plus CDK4/6 inhibitor combination therapy in patients with HR+ metastatic breast cancer. Whole genome sequencing (WGS) of breast cancer models and selected paired primary-metastatic samples demonstrated that active APOBEC3 mutagenesis promoted resistance to both endocrine and targeted therapies through characteristic alterations such as RB1 loss-of-function mutations. Evidence of APOBEC3 activity in pre-treatment samples illustrated a pervasive role for this mutational process in breast cancer evolution. The study reveals APOBEC3 mutagenesis to be a frequent mediator of therapy resistance in breast cancer and highlights its potential as a biomarker and target for overcoming resistance.

  Study EGAS50000000875

• Hematopoietic Tet2 inactivation enhances the response to checkpoint blockade immunotherapy

  Somatic mutations inactivating TET2 are among the most common drivers of clonal hematopoiesis (CH). While TET2 inactivation is associated with monocyte-derived inflammation and improved chimeric antigen-receptor-T cell function, its impact on immunotherapy response is unknown. In murine models, hematopoietic Tet2 mutation enhanced the immune checkpoint blockade (ICB) response via the combined presence of phagocytes, CD4 and CD8 T cells. This effect was lost with myeloid- or T-cell restricted Tet2 inactivation, or in mice with 20% Tet2-mutant hematopoiesis. Mechanistically, in Tet2-mutant tumor-infiltrating leukocytes (TILs), ICB preferentially restricted cell states linked to tumor progression while inducing anti-tumor states. Tet2-mutant monocytes activated costimulatory programs, while Tet2-mutant T cells showed enhanced T cell memory signatures, alongside decreased exhaustion and regulatory phenotypes. This murine data was clinically relevant, since tumors from colorectal cancer and melanoma patients with TET2 driver mutation-CH (TET2-CH) showed enhanced immune infiltration, inflammation, and T cell activation. In melanoma patients treated with ICB, TET2-CH was associated with 6-fold greater odds of clinical benefit. Collectively, we demonstrate that hematopoietic Tet2 inactivation primes leukocytes for anti-tumor states associated with immunotherapy response and provides a potential biomarker for personalized therapy.

  Study EGAS50000001191

• Drug Development against Tumor Microtube Networks in Glioblastoma

  Glioblastomas are aggressive brain tumors that resist current therapies by hijacking neurodevelopmental programs to invade the brain and form resilient multicellular networks. Tumor Microtubes (TMs), neurite-like processes extended by cancer cells, drive this pathobiology, but TM-inhibiting drugs are lacking. Here we developed a comprehensive, combined in vitro and in vivo drug screening approach, including novel AI-based analysis tools of TM parameters. Two Protein Kinase C (PKC) activators, TPPB and PMA, exhibited robust inhibition of TM formation in cell-based assays. Since unconnected glioblastoma cells exhibit increased sensitivity to cytotoxic therapy in vivo, brain tumor-bearing mice received TPPB, the most consistent TM-inhibiting compound, in combination with radiotherapy. Long-term intravital 2-photon microscopy confirmed anti-TM and anti-tumor effects of the co-treatment. Mechanistically, TPPB treatment inhibited several key molecular and functional components of TM-related glioblastoma cell network communication and decreased the expression of TTYH1, a key driver of invasive TMs. By establishing a novel screening pipeline for anti-TM drug development, our study identifies a druggable TM master regulator pathway and a small molecule with robust anti-TM effects. It provides a proof-of-principle that disrupting TMs and disconnecting tumor networks is a pharmacologically feasible strategy to make therapies against glioblastoma more effective.

  Study EGAS50000000477

• Recurrent intra-tumour heterogeneity is a hallmark of metastatic prostate cancer

  The evolution of cancers from low grade to metastatic disease is a major determinant of cancer mortality. Cancer evolution involves a complex interplay between cell intrinsic genetics and transcriptional alterations and the extrinsic microenvironment. To define the key mechanisms underpinning metastatic development, we focused on the most lethal form of prostate cancer, metastatic castration-resistant prostate cancer, and employed single-cell multi-omics and whole-genome sequencing to deeply profile 34 metastatic lesions from 9 patients obtained from rapid autopsy. We found that intra-tumour heterogeneity is an indicator of key evolutionary processes, characterised by recurrent tumour populations acting as critical functional components of the tumour ecosystem, irrespective of clonal and microenvironmental backgrounds. Unexpectedly, microenvironments across metastatic sites played a limited role while clonal evolution primarily promoted transcriptional noise. Intra-patient functional convergence of tumour ecosystems was observed across metastatic sites, pointing towards system-level selection pressures that drive the heterogeneity landscape of mCRPC. Our findings reveal functional evolutionary convergence of metastatic disease into units of intra-tumour heterogeneity identifying critical determinants of metastatic progression for therapeutic targeting. Note: matched WGS for some tumour single-cell experiments, as well as germline WGS for patients CA0027, CA0034, CA0035, CA0043, and CA0076, can be found in EGAS00001006598.

  Study EGAS50000001312

• Bruno et al.: Interferon gamma rebalances immunopathological signatures in Chronic Granulomatous Disease through metabolic rewiring

  Chronic granulomatous disease (CGD) is a primary immunodeficiency characterized by recurrent life-threatening infections and hyperinflammatory complications. It is caused by mutations in the NADPH oxidase complex and the consequent loss of reactive oxygen species (ROS) production. Recombinant human interferon gamma (rIFN-γ) prophylaxis reduces the risk of severe infections, but the mechanisms behind its efficacy in CGD are still an open question, as it does not restore NADPH oxidase-dependent ROS production. Here, we show that myeloid innate immune cells of CGD patients are transcriptionally and functionally reprogrammed to a hyperactive inflammatory status, displaying an impaired in vitro induction of trained immunity. CGD monocytes have reduced intracellular amino acids concentrations and profound functional metabolic defects, both at the level of glycolysis and mitochondrial respiration. Ex vivo and in vivo treatment with IFN-γ restored these metabolic defects and reduced excessive IL-1β and IL-6 production in response to fungal stimuli in CGD monocytes. These data suggest that prophylactic rIFNγ modulates the metabolic status of innate immune cells in CGD. These data shed light on the effects of NADPH-derived ROS deficiency to the metabolic programs of immune cells and pose the basis for targeting this immunometabolic axis, potentially beyond CGD, with IFN-γ immunotherapy.

  Study EGAS00001005463

• Phylogenetic analysis of paired breast carcinomas identifies genetic events associated with clonal recurrence and invasive progression

  Development of ipsilateral breast carcinoma following diagnosis of breast ductal carcinoma in situ (DCIS) has been assumed to represent recurrence of the primary tumour. However, this may not be the case and it is important to know how often recurrences are new tumours. Ipsilateral primary-recurrence pairs (n=78) were sequenced to test their clonal relatedness. Shared genetic events were identified from whole exome sequencing (n=54 pairs) using haplotype-specific copy number and phylogenetic analysis. The remaining pairs were sequenced by a targeted panel or low-coverage whole genome sequencing. We included 32 non-recurrent DCIS to compare recurrent and non-recurrent disease. We found that 7% of DCIS recurrences were non-clonal by whole exome sequencing, indicative of a new breast carcinoma. Lower resolution methods detected a higher non-clonality rate (29%). Comparing primary DCIS with their recurrences found that evolution of DCIS to invasive disease was associated with increased ploidy and copy number events. TP53 mutations were enriched in DCIS with clonal recurrence compared with non-recurrent DCIS. Our results verify that de novo “recurrent tumours” of independent origin occur in patients who may be at high risk. (https://www.biorxiv.org/content/10.1101/2024.05.19.594731v2.full)

  Study EGAS50000001298

• Emergence of oncofetal plasticity is ubiquitous in early colorectal cancers

  Metastasis formation is classically considered a late-stage event in colorectal cancer (CRC) evolution. Yet, the time and patterning by which metastatic competence is acquired remain poorly understood. Here, we show that metastasis-associated oncofetal cell states already emerge at the earliest stages of CRC, concurrent with invasive front formation. However, while necessary for metastasis, we detect them ubiquitously among early non-metastatic cancers, highlighting additional bottlenecks like immune evasion. To understand how oncofetal cells first emerge, we generated multiregional organoid models that reflect successive tumor progression stages within individual early-stage CRCs. Whole-genome sequencing and growth factor-dependency assays exclude tumor cell-intrinsic acquired traits. In contrast, single-cell spatial atlases of the tumor-microenvironment before and after malignant transformation revealed stereotypic patterning of fibroblast subtypes resembling normal tissue architecture, resulting in distinct regional microenvironments. At the onset of malignant growth into the submucosa, the first cancer-associated fibroblasts (CAFs) to appear strongly resemble submucosal trophocytes and colocalize with oncofetal cell states at invasive fronts. Functionally, fibroblast-organoid co-cultures confirm that these trophocyte-like CAFs induce plastic transitioning to oncofetal states. Thus, interactions between tumor and submucosal fibroblasts directly following malignant transformation dictate the timing and location at which oncofetal plasticity first occurs during CRC progression.

  Study EGAS50000001532

• Assessing gene expression profiling from FFPE Patient Samples: A Comparison of Two Library Preparation Approaches and Recommendations

  Next-Generation Sequencing (NGS) has transformed cancer research and clinical practice, with Whole Exome Sequencing (WES) driving advances in mutational profiling and personalized oncology. Yet, transcriptomic signatures remain essential for understanding disease mechanisms, including therapy resistance pathways. RNA sequencing (RNA-seq), however, faces unique challenges when dealing with low-input or degraded RNA, as often found in archival formalin-fixed paraffin-embedded (FFPE) tissues. Although previous studies have compared library preparation protocols, rapidly evolving technologies call for updated evaluations. Here, we present a direct comparison of two FFPE-compatible stranded RNA-seq library preparation kits: TaKaRa™ SMARTer® Stranded Total RNA-Seq Kit v2 (Kit A) and Illumina Stranded Total RNA Prep Ligation with Ribo-Zero Plus (Kit B). Both kits generated high-quality RNA-seq data, yet important differences emerged. Notably, Kit A achieved comparable gene expression quantification to Kit B while requiring 20-fold less RNA input, a crucial advantage for limited samples, albeit with increased sequencing depth. We critically discuss these results in relation to RNA availability, technical performance, cost-effectiveness, processing time, and automation potential, offering practical guidance for selecting optimal RNA-seq strategies in clinical and translational research settings

  Study EGAS50000001066

• Myeloid cell networks govern re-establishment of original immune landscapes in recurrent ovarian cancer

  Immunotherapy has produced disappointing results in recurrent ovarian cancer (OC). However, the prognostic value of tumour-infiltrating lymphocytes (TILs) is largely based on the analysis of treatment-naive tumours. To understand the immunobiology of recurrent cancers, and their evolution, we profiled a total of 697 OC tumour samples (566 primary and 131 recurrent including 69 patient-matched samples) from 595 patients of five independent multi-institutional cohorts. By capturing heterogeneous TIL distributions, we identified four immune phenotypes associated with differential prognosis, TILs states and TILs:myeloid networks, which dictate malignant progression after chemotherapy and recurrence. Notably, recurrent tumours recapitulate the immunogenic landscapes of original cancers. Mirroring inflamed human OC, preclinical recurrent Brca1mut tumours maintained activated TILs:dendritic cells (DCs) niches but evaded immune control through COX/PGE2 signaling upregulation. Conversely, recurrent Brca1wt tumours displayed loss of TILs:DCs niches and accumulated immunosuppressive myeloid networks featuring Trem2/ApoEhigh TAMs and Nduf4l2high/Galectin3high malignant states. Our study highlights that persistent immunogenicity in recurrent OC is governed by the BRCA-dependent crosstalk between dissimilar myeloid cells and TILs. We propose new targets to modulate the malignant and TAM compartments which improve efficacy of standard of care chemotherapy in mouse OC models with HR deficiency.

  Study EGAS50000001069

• LySeqST: A targeted sequencing assay for robust genomic classification of diffuse large B-cell lymphoma

  Diffuse large B-cell lymphoma (DLBCL) is a biologically heterogeneous disease. Two genomic classification systems, LymphGen and DLBclass, are capable of reproducibly classifying single tumors into subtypes with prognostic relevance in the context of standard-of-care R-CHOP chemoimmunotherapy. Preliminary data from subgroup analyses of clinical trials suggest that the distinct biological features are targetable with drugs that exploit subtype-specific vulnerabilities, and large clinical trials are being designed to test this hypothesis. Although whole exome sequencing (WES) remains the gold standard for each of these classification systems, many laboratories may opt for smaller targeted sequencing panels that reduce costs associated with sequencing, data storage, and processing. Here, we describe the design and validation of a targeted sequencing panel (LySeqST) that captures the genomic features required specifically for LymphGen classification. We perform in silico and real-world validation vs. WES data and demonstrate that LySeqST can be used to accurately classify DLBCL tumors. In an unselected population-based cohort, we determine the real-world proportions and validate the prognostic relevance of LymphGen subtypes, confirming that only tumors expressing the dark zone gene expression signature are consistently associated with inferior outcomes. Our results support the use of LySeqST for accurate genomic classification of DLBCL.

  Study EGAS50000001601

• Genentech Colon Cancer Screen

  Identifying and understanding changes in cancer genomes is essential for the development of targeted therapeutics. Here we analyse systematically more than 70 pairs of primary human colon tumours by applying next-generation sequencing to characterize their exomes, transcriptomes and copy-number alterations. We have identified 36,303 protein-altering somatic changes that include several new recurrent mutations in the Wnt pathway gene TCF7L2, chromatin-remodelling genes such as TET2 and TET3 and receptor tyrosine kinases including ERBB3. Our analysis for significantly mutated cancer genes identified 23 candidates, including the cell cycle checkpoint kinase ATM. Copy-number and RNA-seq data analysis identified amplifications and corresponding overexpression of IGF2 in a subset of colon tumours. Furthermore, using RNA-seq data we identified multiple fusion transcripts including recurrent gene fusions involving R-spondin family members RSPO2 and RSPO3 that together occur in 10% of colon tumours. The RSPO fusions were mutually exclusive with APC mutations, indicating that they probably have a role in the activation of Wnt signalling and tumorigenesis. Consistent with this we show that the RSPO fusion proteins were capable of potentiating Wnt signalling. The R-spondin gene fusions and several other gene mutations identified in this study provide new potential opportunities for therapeutic intervention in colon cancer.

  Study EGAS00001000288

• Exome sequencing identifies mutation of the ribosome in T-cell acute lymphoblastic leukemia

  Acute lymphoblastic leukemia (ALL) is the most common pediatric malignancy and responds well to therapy in children, but not in adults. We used exome sequencing of 67 T-cell ALL (T- ALL) samples to gain insight in the mutational spectrum and age related differences of these leukemias. We detected protein-altering mutations in 508 genes, with an average of 8.2 mutations in pediatric and 21.0 in adult T-ALL. Based on stringent filtering, we predict 15 of these genes to be drivers, including 7 genes that were not previously implicated in T-ALL. We identify CNOT3 as a new tumor suppressor, as it is mutated in 12.2 % of adult T-ALLs and its knock-down causes de novo tumors in a sensitized drosophila eye tumor model. In addition, we identify mutations in RPL5 and RPL10, proteins of the 60S ribosomal complex2, in 18.2 % of pediatric T-ALLs. Mutation of arginine 98 in the X-linked RPL10 gene was most frequent. Yeast cells expressing the Rpl10 R98S mutant showed a clear ribosome biogenesis defect, excluding that this mutation is a passenger mutation. Our data provide insights in the mutational landscape of pediatric and adult T-ALL and identify the ribosome as a potential oncogenic factor in cancer.

  Study EGAS00001000296

• Single cell sequencing of breast cancer T cells reveals a tissue-resident memory subset associated with improved prognosis

  The quantity of tumor infiltrating lymphocytes (TILs) in breast cancer (BC) is a robust prognostic factor for improved patient survival particularly in triple negative and HER2-overexpressing BC subtypes. Whilst T cells are the predominant TIL population, the relationship between quantitative and qualitative differences in T cell subpopulations and patient prognosis remains unknown. We performed single-cell RNA sequencing (scRNA-seq) of 6311 T cells isolated from human BCs, and show for the first time that significant heterogeneity exists in the infiltrating T cell population. We demonstrated that high TIL BCs contained significant numbers of CD8+ T cells with features of tissue-resident memory T (TRM) cell differentiation and that these CD8+ TRM cells expressed high levels of immune checkpoint molecules and effector proteins. A CD8+ TRM gene signature developed from the scRNA-seq data was significantly associated with improved patient survival in early-stage triple negative breast cancer (TNBC) and provided better prognostication compared with CD8 expression alone. Our data suggests that CD8+ TRM contribute to BC immunosurveillance and are the key targets of modulation through immune checkpoint inhibition. Further understanding of the development, maintenance and regulation of TRM cells will be crucial for successful immunotherapeutic development in BC.

  Study EGAS00001002845

• Finding the way towards the eradication of therapy-related myeloid neoplasms

  Therapy-related neoplasms are directly caused by the drugs used to eradicate a previous cancer. It is well known that a subset of drugs designed to kill cancer cells through DNA-damage-induced apoptosis can intrinsically, as secondary effect, cause DNA-damage in other cells, thus inducing a neoplasm. Those neoplasms have especially bad prognosis. Why are these drugs still in use? Because for many diseases we do not have anything better, at the moment. Approximately 0.5-1% of the kids treated for various types of cancers develop a therapy-related myeloid neoplasm (tMN) as a result of the treatment. Despite having the same origin, their prognostic and molecular features are different from the adult ones, and a complete understanding of the pathogenesis of tMN in children is still missing. A total of 84 paediatrics neoplasms have been profiled by either Whole Genome or Whole Exome sequencing by Jason R. Schwartz and his colleagues, in a collaborative effort guided at the St. Jude Children's hospital, in Memphis. They identified a mean of 28 somatic mutations per patient, a number which is significantly bigger than the number of mutations found in paediatrics primary neoplasms, as depicted in the original paper's figure 1b, shown here. Their sequencing data identified that Ras/MAPK pathway mutations, alterations in RUNX1 or TP53, and KMT2A rearrangements as the most frequent drivers of the disease. Clonal evolution analysis (the comparative study of the behaviour of cells harbouring different mutations) highlighted that, unlike in adults’ tMN, there were not pre-existing leukemic clones, and virtually all mutated cells were found arising as a direct consequence of cytotoxic therapy. This represents yet another important difference between adult and paediatric tMN. This important study have been published in Nature communications in February 2021. The sequencing data have been deposited at EGA under the study number EGAS00001004850. Thanks to this, other researchers around the world could re-use it to uncover more information, integrate and validate their studies. Revealing the molecular features underlying the development of cancers is the only way towards the development of new generation targeted therapies, which will allow us to finally drop the cytotoxic drugs eliciting such unwanted and catastrophic secondary effects as tMN. It is not an overstatement to say that, while improving the treatment of primary tumours, we could entirely eradicate therapy-related neoplasms right from their origin. Looking forward.

  Blog eradication-of-therapy-related-myeloid-neoplasms

• Recurrent somatic JAK-STAT mutations within a novel RUNX1-mutated pedigree

  The acquisition of somatic mutations is an emerging field of investigation in familial leukemia. Currently, genetic profiles in familial MDS/AML are considered analogous to sporadic disease, although the patterns of clonal evolution within families are poorly defined. We performed whole exome profiling of tumour samples from a novel RUNX1 mutated family, to determine the stepwise evolution of MDS/AML across 4 young siblings. Three siblings developed monocytic AML/RAEB2 at 5 years of age, with hepatosplenomegaly and somatic mutations upregulating JAK-STAT signalling, the latter are typically detected in <5% of sporadic MDS/AML. Two siblings acquired the canonical JAK2 V617F mutation, while another acquired a unique missense mutation of SH2B3, a negative regulator of JAK2. Notably, 2/3 siblings demonstrated dosage amplification of these mutations due to acquired uniparental disomy of chromosomes 9p and 12q (encompassing JAK2 and SH2B3, respectively). All 4 siblings were heterozygous for the 46/1 JAK2 haplotype associated with predisposition to sporadic V617F myeloproliferative disorders, which likely influenced the acquisition of JAK2 mutations. Our findings provide further evidence that relatives with shared germline mutations may acquire somatic mutations in a non-random manner leading to convergent patterns of disease evolution.

  Study EGAS00001001862

• Comparison of EGF and PDGF driven glioblastomas.

  Glioblastoma multiforme (GBM) is the most common and aggressive type of primary brain tumor. Epidermal growth-factor (EGF) and platelet-derived growth-factor (PDGF) receptors are frequently amplified and/or possess gain-of-function mutations in these tumors. Despite years of research into tyrosine-kinase inhibitors, clinical trials have shown disappointing efficacy. One of the underlying reasons is intra-tumor heterogeneity. To assess the effect of intra-tumor heterogeneity on gene expression, we derived a general approach to map single-cell expression profiles to sequentially-acquired mutations identified from exome sequencing. Using 288 single cells, we constructed high-resolution phylogenies of EGF-driven and PDGF-driven GBMs, modeling transcriptional kinetics during tumor evolution. A deletion in a PDGF-receptor domain involved in dimerization correlated to an up-regulation of the PI3K/AKT pathway in a proneural GBM, and ectopic expression of this mutant PDGFRA in vitro enhanced proliferation and invasiveness. Descending the phylogenetic tree of the PDGFR-driven tumor corresponded to a progressive induction of an oligodendrocyte progenitor-like cell type that expresses growth and angiogenic factors. In contrast, phylogenetic analysis of an EGFR-amplified tumor showed a gradual up-regulation of pro-invasive genes.

  Study EGAS00001001900

• Molecular characteristics in Burkitt lymphoma over age groups

  Treatment optimization has improved the overall survival in Burkitt lymphoma (BL) patients but there are still relevant outcome differences between pediatric and adult patients. In both age groups survival remains poor for patients who suffer relapse. Therefore, new treatment approaches are warranted to improve survival. Thus, we conducted a comparative study in the largest pediatric cohort to date including 191 cases and 97 samples from adult patients and analyzed the mutational landscape. We detected overlapping molecular profiles but also significant differences in the mutational pattern. Samples from pediatric patients show a higher frequency of mutations in ID3, DDX3X and BAF-complex components ARID1A and SMARCA4, while YY1AP1 and BCL2 are more frequently mutated in cases from adult BL patients. Correlation of the mutational profile with clinical data was available for the pediatric cohort and identified TP53 mutations to be highly significantly associated with the incidence of relapse thus identifying the first molecular marker to be associated with the incidence of relapse in pediatric BL. As the hotspot analysis of TP53 mutations showed mostly gain of function mutations, it is tempting to speculate that drugs decreasing the amount of malignant TP53 or increase the abundance of wildtype TP53 could be a promising approach for future treatment regimens to further improve prognosis.

  Study EGAS00001005270

• The PEMDAC phase 2 study of pembrolizumab and entinostat in patients with metastatic uveal melanoma

  Patients with metastatic uveal melanoma (UM) have an abysmal prognosis. Preclinical studies have suggested that epigenetic therapy could enhance immunogenicity of cancer cells. Here we test if epigenetic therapy would enhance PD-1 immunotherapy in patients with metastatic UM. We report the results of the PEMDAC phase 2 clinical trial (n=29; NCT02697630) where the HDAC inhibitor entinostat was combined with the PD-1 inhibitor pembrolizumab in patients with metastatic UM. The primary endpoint was objective response rate (ORR), and was met with an ORR of 14%. The clinical benefit rate at 18 weeks was 28%, median progression free survival was 2.1 months and the median overall survival was 13.4 months. Toxicities were manageable, and there were no treatment-related deaths. Extensive genomics studies were performed using DNA/RNA and single cell sequencing and flow cytometry. Objective responses and prolonged survival were seen in three patients with BAP1 wildtype tumors, and in one patient with an iris melanoma that exhibited a UV signature. Longer survival also correlated with low baseline ctDNA levels or LDH. In conclusion, HDAC inhibition and anti-PD1 immunotherapy results in durable responses in a subset of patients with metastatic UM. Further exploration of combined immunotherapy and epigenetic therapy in metastatic UM is warranted.

  Study EGAS00001005478

• The Genomic Landscape of Childhood and Adolescent Melanoma

  Despite remarkable advances in the genomic characterization of adult melanoma, the molecular pathogenesis of pediatric melanoma remains largely unknown. We analyzed 15 conventional melanomas (CM), 3 melanomas arising in congenital nevi (CNM), and 5 spitzoid melanomas (SM), using various platforms, including whole genome or exome sequencing, molecular inversion probe assay, and/or targeted sequencing. CM demonstrated a high burden of somatic single nucleotide variations (SNV), with each case containing a TERT promoter (TERT-p) mutation, 13/15 containing an activating BRAF V600 mutation, and >80% of the identified SNV consistent with UV damage. By contrast, the 3 CNM contained an activating NRAS Q61 mutation and no TERT-p mutations. SM were characterized by chromosomal rearrangements resulting in activated kinase signaling in 40%, and an absence of TERT-p mutations except for the one SM that succumbed to hematogenous metastasis. We conclude that pediatric CM has a very similar UV-induced mutational spectrum to that found in the adult counterpart, emphasizing the need to promote sun protection practices in early life and to improve access to therapeutic agents being explored in adults in young patients. By contrast, the pathogenesis of CNM appears to be distinct. TERT-p mutations may identify the rare subset of spitzoid melanocytic lesions prone to disseminate.

  Study EGAS00001000901

• The evolutionary steps from primary to metastatic prostate cancer are largely uncharted, and the ability to use DNA present in body fluids as correlates of aggregate metastatic status is under-examined. We reconstructed phylogenies in ten prostate cancer patients with fatal disease using deep targeted sequencing of the prostate, adjacent and distant organs, as well as plasma, serum, and cerebrospinal fluid at various time points. A total of 163 samples are studied.

  The evolutionary steps from primary tumor to metastasis in prostate cancer are largely uncharted, and the ability to use serum, plasma, and cerebrospinal fluid as a correlate of aggregate metastatic tumor genomic status has not been tested. We used deep targeted sequencing to reconstruct tumor evolution in ten prostate cancer patients with fatal disease encompassing examination of the prostate and adjacent and distant organs, as well as plasma, serum, and cerebrospinal fluid at various time points. We show that there is substantial evolution from a common ancestor within the prostate that results in branching to multiple lineages which form an intermixed multi-clonal primary tumor mass. After the occurrence of key driver aberrations, one of these lineages will metastasize to multiple sites in a sequential fashion. These metastatic sites are then susceptible to being populated by cells from other intra-prostatic lineages or from other metastases. Genomic representation of metastases in body fluids is not uniform. Cerebrospinal fluid analysis can detect lineages not detected in circulating DNA, suggesting possible clinical utility.

  Study EGAS00001003848

• Pancreatic, Small-intestinal and Pulmonary Neuroendocrine Tumors

  Purpose: The primary origin of neuroendocrine tumor metastases can be difficult to determine by histopathology alone, but is critical for therapeutic decision making. DNA methylation-based profiling is now routinely used in the diagnostic workup of brain tumors. This has been enabled by the availability of cost-efficient array-based platforms. We have extended these efforts to augment histopathological diagnosis in neuroendocrine tumors. Experimental Design and Results: We compiled data of 69 small-intestinal, pulmonary, and pancreatic neuroendocrine tumors. These data were used to build a ridge regression calibrated random forest classification algorithm (NEN-ID) that predicts the origin of tumor samples with high accuracy (> 95%). The model was validated during 3x3 nested cross validation and tested in a local (n=26) and external (n=172) cohort. In addition, we show that our diagnostic approach is robust across a range of possible confounding experimental parameters such as tumor purity and array quality. A software infrastructure and online user interface was built to make the model available to the scientific community. Conclusions: This DNA methylation-based prediction model can be used in the workup for patients with neuroendocrine tumors of unknown primary. To facilitate validation and clinical implementation, we provide a user-friendly, publicly available web-based version of NEN-ID.

  Study EGAS00001004878

• RNA-Seq of vocal fold fibroblasts in Reinke’s edema and control subjects

  The voice disorder Reinke’s edema (RE) is a smoking- and voice-abuse associated benign lesion of the vocal folds, defined by an edema of the Reinke space, accompanied with pathological microvasculature changes and immune cell infiltration. Vocal fold fibroblasts (VFF) are the main cell type of the lamina propria and play a key role in the disease progression. Current therapy is restricted to symptomatic treatment. Hence, there is an urgent need for a better understanding of the molecular causes of the disease. In the present study, we investigated differential expression profiles of RE and control VFF by means of RNA sequencing. In addition, fast gene set enrichment analysis (FGSEA) was performed in order to obtain involved biological processes, mRNA and protein levels of targets of interest were further evaluated. We identified 74 differentially regulated genes in total, 19 of which were upregulated and 55 downregulated. Differential expression analysis and FGSEA revealed upregulated genes and pathways involved in extracellular matrix (ECM) remodeling, inflammation and fibrosis. Downregulated genes and pathways were involved in ECM degradation, cell cycle control and proliferation. The current study addressed for the first time a direct comparison of VFF from RE to control and evaluated immediate functional consequences.

  Study EGAS00001005130

• Sclerosing epithelioid fibrosarcoma case report

  Sclerosing epithelioid fibrosarcoma (SEF) is a rare and aggressive soft tissue sarcomathought to originate in fibroblasts of the tissues surrounding tendons, ligaments and muscles. Minimally responsive to conventional cytotoxic chemotherapies, greater than 50% of SEF patients experience relapse and/or metastatic disease. SEF is most commonly discovered in middle-aged and elderly adults, but also rarely in children. A common gene fusion occurring between the EWSR1 and CREB3L1 genes has been observed in 80-90% of SEF cases. We describe here the youngest SEF patient reported to date (a 3-year-old Caucasian male) who presented with numerous bony and lung metastases. Additionally, we perform a comprehensive literature review of all SEF-related articles published since the disease was first characterized. The patient described in this report experienced persistent disease progression despite being heavily treated with multiple surgeries, radiation, numerous chemotherapies and targeted therapeutics. Primary tumor, metastatic and relapse sites were sequenced by whole genome, whole exome, and deep transcriptome next generation sequencing with comparison to a patient-matched normal blood sample. Consistent across all sequencing analyses was the disease-defining EWSR1-CREB3L1 fusion as a single feature consensus. We provide an analysis of our genomic findings and a discussion as to potential therapeutic strategies for SEF.

  Study EGAS00001005214

• IDENTIFICATION AND TARGETED MANAGEMENT OF A PATIENT WITH A NEURODEGENERATIVE DISORDER CAUSED BY BIALLELIC MUTATIONS IN SLC5A6

  We describe a sibling pair displaying an early infantile-onset, progressive neurodegenerative phenotype, with symptoms of developmental delay and epileptic encephalopathy developing from 12 to 14 months of age. Using whole exome sequencing, compound heterozygous variants were identified in SLC5A6, which encodes the sodium-dependent multivitamin transporter (SMVT) protein. SMVT is an important transporter of the B-group vitamins biotin, pantothenate and lipoate. The protein is ubiquitously expressed and has major roles in vitamin uptake in the digestive system, as well as transport of these vitamins across the blood-brain barrier. Pathogenicity of the identified variants was demonstrated by impaired biotin uptake of mutant SMVT. Identification of this vitamin transporter as the genetic basis of this disorder guided targeted therapeutic intervention, resulting clinically in improvement of the patient’s neurocognitive and neuromotor function. This is the second report of biallelic mutations in SLC5A6 leading to a neurodegenerative disorder due to impaired biotin, pantothenate and lipoate uptake. The genetic and phenotypic overlap of these cases confirms mutations in SLC5A6 as the genetic cause of this disease phenotype. Recognition of the genetic disorder caused by SLC5A6 mutations is essential for early diagnosis and to facilitate timely intervention by triple vitamin (biotin, pantothenate and lipoate) replacement therapy.

  Study EGAS00001003861

• Patient-derived pediatric brain tumor organoids faithfully recapitulate primary tumors

  Extensive molecular analysis by many groups has revealed a heterogeneous landscape of embryonal brain tumors, including medulloblastoma (MB) and atypical teratoid/rhabdoid tumor (ATRT). MB, the most common malignant pediatric brain tumor, is comprised of four molecularly distinct subgroups: Wingless (WNT), Sonic Hedgehog (SHH), Group 3 (G3), and Group 4 (G4). Among these, MB-SHH harboring MYCN amplification and TP53 mutation and MB-G3 with MYC overexpression or amplification constitute the most aggressive tumors, have the poorest clinical outcomes, and lack effective therapeutic options for recurrent disease. Here we describe the establishment of tumor organoid (TO) models from patient-derived orthotopic xenografts (PDOX) of MYCN-amplified and TP53-mutant MB-SHH, MYC-amplified MB-G3 and two ATRT which all grew in mice when reimplanted into the cortices of naïve recipient animals (TOX). Methylation, bulk and single-cell RNA sequencing, and whole genome sequencing demonstrated that TOs and TOXs faithfully recapitulated the molecular landscape and heterogeneity of PDOXs and primary patient tumors compared to established cell lines. The development of TOs will facilitate in vitro high-throughput drug or CRISPR screens without the need for use tumor-bearing mice. This will also accelerate the identification and validation of new vulnerabilities and in vivo preclinical testing toward clinical trials.

  Study EGAS00001008305

• Characterizing microbiome-directed fibre snacks in gnotobiotic mice and humans

  Knowledge of the interrelationships between what we eat and the configurations of our gut microbial communities is providing important insights into how food components that are not directly metabolized by human enzymes are linked to our physiology and health status. Changing food preferences brought about by Westernization that have deleterious health effects1,2, plus rapid population expansion, ongoing challenges to sustainable agriculture, and other forces contributing to increased food insecurity, are catalyzing efforts to identify more nutritious and affordable foods3. The gut microbial community is complex, dynamic, and exhibits considerable intra- and interpersonal variation in its composition and functions. The massive number of potential interactions between its components makes it challenging to define the mechanisms by which food ingredients affect community properties. There is also a paucity of information about the ‘bioactive’ ingredients of foods that influence the fitness and expressed functions of community members. Here, plant fibres, from different sustainable sources and targeting distinct features of obese human gut microbiomes in gnotobiotic mice, were formulated into snack prototypes and used to supplement controlled diets consumed by overweight and obese adults; the results revealed fibre-specific changes in their microbiomes that were linked to changes in their plasma proteomes indicative of altered physiologic state.

  Study EGAS00001005268

• Diverse_outcomes_of_controlled_human_malaria_infection_originate_from_host_intrinsic_immune_variation_and_not_var_gene_switching

  Falciparum malaria is clinically heterogeneous and the relative contribution of parasite and host factors in shaping disease severity remains unclear. We set out to explore the interaction between host inflammation and parasite variant surface antigen (VSA) expression, asking whether this relationship underpins the variation observed in controlled human malaria infection. We uncovered marked heterogeneity in the response of naive hosts to blood challenge with some volunteers maintaining a state of immune quiescence, others triggering interferon-stimulated inflammation and a small group showing transcriptional evidence of myeloid cell suppression. Significantly, only inflammatory volunteers experienced hallmark symptoms of clinical malaria. When we then tracked temporal changes in parasite VSA expression to ask whether variants associated with severe disease preferentially expand in naive hosts (as predicted by current theory) we found that var gene profiles were essentially unchanged after 10-days of blood-stage infection. These data therefore show that the diverse clinical outcomes of CHMI largely originate from host-intrinsic variation and there is no evidence for switching or selection of var genes in naïve hosts. 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 EGAS00001003766

• Tumor-derived exosomes modulate PD-L1 expression in monocytes

  In chronic lymphocytic leukemia (CLL), monocytes and macrophages are skewed toward protumorigenic phenotypes, including the release of tumor-supportive cytokines and the expression of immunosuppressive molecules such as programmed cell death 1 ligand 1 (PD-L1). To understand the mechanism driving protumorigenic skewing in CLL, we evaluated the role of tumor cell–derived exosomes in the cross-talk with monocytes. We carried out RNA sequencing and proteome analyses of CLL-derived exosomes and identified noncoding Y RNA hY4 as a highly abundant RNA species that is enriched in exosomes from plasma of CLL patients compared with healthy donor samples. Transfer of CLL-derived exosomes or hY4 alone to monocytes resulted in key CLL-associated phenotypes, including the release of cytokines, such as C-C motif chemokine ligand 2 (CCL2), CCL4, and interleukin-6, and the expression of PD-L1. These responses were abolished in Toll-like receptor 7 (TLR7)–deficient monocytes, suggesting exosomal hY4 as a driver of TLR7 signaling. Pharmacologic inhibition of endosomal TLR signaling resulted in a substantially reduced activation of monocytes in vitro and attenuated CLL development in vivo. Our results indicate that exosome-mediated transfer of noncoding RNAs to monocytes contributes to cancer-related inflammation and concurrent immune escape via PD-L1 expression.

  Study EGAS00001002377

• Unifying recovery dynamics in heterogeneous diseases exemplified by COVID-19

  Studying the recovery dynamics in different diseases poses many difficulties, as patients often display high heterogeneity in their recovery courses. Moreover, most attempts to study disease dynamics focus on disease progression rather than disease recovery mechanisms. To model recovery dynamics, using severe COVID-19 as the example, we align heterogeneous recovery trajectories via a novel computational scheme applied to longitudinally sampled blood transcriptomes. We thus generate pseudotime trajectories, which we then link to cellular and molecular mechanisms based on cell deconvolution analysis and molecular pathway prediction, thus presenting a unique framework for studying recovery processes over time. Specifically, mature neutrophils displayed a gradual decrease during recovery, allowing superior useability for outcome prediction compared to currently used clinical markers. Further, we discovered a recovery-related regulatory change in gene programs resulting in immune rebalancing between interferon and NFkB activity and the restoration of cell homeostasis. We thus propose regulatory mechanisms governing COVID-19 recovery and suggest mature neutrophils and additional gene markers, as novel clinical biomarkers for disease outcome. Overall, we present a novel clinically relevant computational framework for modeling disease recovery, paving the way for future studies of the recovery dynamics in other diseases and tissues.

  Study EGAS00001005735

• The Spatial Heterogeneity in Multiple Myeloma - from the Subclonal Architecture to the Immune Microenvironment (partly hipo_K08K)

  Multiple myeloma (MM) cells show pronounced heterogeneity not only within a given patient but also between spatially separated tumor loci in the bone marrow. Understanding this spatial heterogeneity is emerging as a critical challenge for the successful treatment of the disease. Yet, our understanding of spatial differences in the subclonal architecture, molecular signatures and interactions with the tumor microenvironment remains very limited. To address this shortcoming, we performed bulk and single-cell multi-region sequencing, including random bone marrow samples from the iliac crest and paired imaging-guided focal lesion specimens from 15 newly diagnosed MM patients. We found a median of 5 subclones per patient and unique subclones in focal lesions. Central features of spatial heterogeneity included a consistent down-regulation of the chemoattractant cytokines CXCL7 and CXCL12 in focal lesion tumor cells as well as a significant depletion of macrophages in the focal lesion microenvironment. In contrast, a site-specific expansion of T-cell clones was not detected in focal lesions. In conclusion, our results demonstrate the relevance of considering spatial heterogeneity with potential implications for immunotherapies and molecular studies analyzing the role of gene signatures and MM-microenvironment interactions in disease progression.

  Study EGAS00001006090

• UK10K NEURO ABERDEEN

  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. This sample set comprises cases of schizophrenia with additional cognitive measurements, collected in Aberdeen, Scotland.For further information on the Aberdeen cohort please contact David St Clair (d.stclair@abdn.ac.uk).

  Study EGAS00001000109

• Ethiopia_Genome_Project__high_coverage_

  The Ethiopian area stands among the most ancient ones ever occupied by human populations and their ancestors. Particularly, according to archaeological evidences, it is possible to trace back the presence of Hominids up to at least 3 million years ago. Furthermore, the present day human populations show a great cultural, linguistic and historic diversity which makes them essential candidate to investigate a considerable part of the African variability. Following the typing of 300 Ethiopian samples on Illumina Omni 1M (see Human Variability in Ethiopia project, previously approved by the Genotyping committee) we now have a clearer idea on which populations living in the area include the most of the diversity.This project therefore aims to sequence the whole genome of 300 individuals at low (4-8x) depth belonging to the six most representative populations of the Ethiopian area to produce a unique catalogue of variants peculiar of the North East Africa. Furthermore 6 samples (one from each population) will also be sequenced at high (30x) depth to ensure full coverage of the diversity spectrum.The retrieved variants will be of great help in evaluating the demographic dynamics of those populations as well as shedding light on the migrations out of Africa.

  Study EGAS00001000237

• Exome_trios_in_patients_with_gastroschisis

  Gastroschisis (MIM 230750) is a herniation of the intestines through a defect of the abdominal wall lateral to the umbilicus (usually on the right side), and it is not covered by a membrane [Ledbetter, 2012]. Gastroschisis is a congenital anomaly with increasing incidence, easy prenatal diagnosis and extremely variable postnatal outcomes. On the basis of clinical manifestations, epidemiologic charateristics, and the presence and type of additional malformations, gastroschisis could be considered a heterogeneous condition with no gene/s discovered yet. This congenital anomaly affects approximately 1-3 infancts per 10,000 live births [Calzolari et al.1995;Parker et al.,2010] Current knowledge about causative mutations/variants. To date, no single gene has been linked to gastroschisis. Some publications have tried to link this malformation to variants in genes (such as AEBP1 (adipocyte enhancer binding protein) gene [Feldkamp et al,. 2012] or the VEGF-NOS3 pathway [Lammer et al., 2008]. Previously, a Scribble mutant mouse model (circletail) was reported to exhibit gastroschisis, however recent studies demonstrated that the Scribble knockout fetus exhibits exomphalos phenotype of gastroschisis [Carnagham et al., 2013]. 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 EGAS00001002664

• Multi-organ landscape of therapy-resistant melanoma

  Metastasis and failure of present-day therapies represent the most common causes of mortality in patients with cutaneous melanoma. To identify the underlying genetic and transcriptomic landscapes, here we analyzed multi-organ metastases and tumor-adjacent tissues from eleven rapid autopsies after treatment with MAPK inhibitor (MAPKi) and/or immune checkpoint blockade (ICB) and death due to acquired resistance. Either treatment elicits shared genetic alterations that suggest immune-evasive, cross-therapy resistance mechanisms. Large, non-clustered deletions, inversions, and inter-chromosomal translocations dominate rearrangements. Using 345 therapy-naïve TCGA melanoma and 131 patient-matched melanoma pre-and-post acquired resistance to only either treatment, we performed cross-cohort analyses to identify MAPKi and ICB as respective contributors to gene amplifications and deletions enriched in autopsy versus therapy-naïve tumors. Private/late mutations and structural variants display shifted mutational and rearrangement signatures, with MAPKi specifically selecting for signatures of defective homologous-recombination, mismatch, and base-excision repair. Transcriptomic signatures and cross-talks with tumor-adjacent macroenvironments nominate organ-specific adaptive pathways. An immune-desert, CD8+-macrophage-biased archetype, T-cell exhaustion, and type-2 immunity characterize the immune contexture. This multi-organ analysis of lethal, therapy-resistant melanoma presents preliminary insights with potential to improve therapeutic strategies.

  Study EGAS00001006644

• Allele Balance Bias Identifies Systematic Genotyping Errors and False Disease Associations

  In recent years, Next Generation Sequencing (NGS) has become a cornerstone of clinical genetics and diagnostics. Many clinical applications require high precision, especially if rare events such as somatic mutations in cancer or genetic variants causing rare diseases need to be identified. Although random sequencing errors can be modeled statistically and deep sequencing minimizes their impact, systematic errors remain a problem even at high depth of coverage. Understanding their source is crucial to increase precision of clinical NGS applications. In this work, we studied the relation between recurrent biases in allele balance (AB), systematic errors and false positive variant calls across a large cohort of human samples analyzed by whole exome sequencing (WES). We have modeled the allele balance distribution for biallelic genotypes in 987 WES samples in order to identify positions recurrently deviating significantly from the expectation, a phenomenon we termed allele balance bias (ABB). Furthermore, we have developed a genotype callability score based on ABB for all positions of the human exome, which detects false positive variant calls that passed state-of-the-art filters. Finally, we demonstrate the use of ABB for detection of false associations proposed by rare variant association studies (RVAS).

  Study EGAS00001003027

• Clinical outcomes and immune correlates of response to nivolumab plus chemoradiotherapy in women with locally-advanced cervical cancer – NiCOL study

  Concurrent chemoradiotherapy (CRT) with blockade of the PD-1 pathway may enhance immune-mediated tumor control through increased phagocytosis, cell death, and antigen presentation. The NiCOL phase 1 trial was designed to determine safety, tolerance, and immune correlates of response to the PD-1 blocking antibody, nivolumab, with and following concurrent CRT in 16 women with locally advanced cervical cancer (LACC). Three patients experienced grade 3 dose-limiting toxicities. The overall response rate was 93.8% [95%CI: 69.8%-99.8%] with a 2-year PFS of 75% [95% CI: 56.5%-99.5%]. Compared to patients with progressive disease (PD), progression-free (PF) subjects showed a brisker stromal immune infiltrate, higher proximity of tumor-infiltrating CD3+ T cells to PD-L1+ tumor cells and of Foxp3+ T cells to proliferating CD11c+ myeloid cells. PF had higher baseline levels of PD-1 on EMRA CD4+ T cells and ICOS-L on tumor-associated macrophages (TAMs) vs PD, which instead displayed enhanced PD-L1 expression on TAMs and higher, circulating frequencies of proliferating Tregs at baseline. Concomitant nivolumab plus definitive CRT is safe and associated with encouraging PFS rates. Further validation in the subset of LACC displaying pre-existing, adaptive immune activation is warranted.

  Study EGAS00001007297

• The Genomic Landscape of Core-Binding Factor Acute Myeloid Leukemias

  Pediatric de novo acute myeloid leukemia (AML) is a heterogeneous disease that can be divided into clinically distinct subtypes based on the presence of specific chromosomal abnormalities or gene alterations. One of the best characterized subtypes of AML involves leukemias with alterations of the core-binding factor (CBF)-complex, which comprises the FAB subtypes M2 and M4Eo and associates with a favorable outcome. Patients with the AML M2 subtype harbor a translocation between chromosomes 8 and 21 [t(8;21)] that yields the chimeric fusion gene RUNX1(AML1)-RUNX1T1(ETO), while patients with AML M4Eo express the chimeric fusion gene CBFβ-SMMHC(MYH11) as a result of an inversion/translocation event of chromosome 16 [inv(16)/t(16;16)]. In an effort to define the total complement of genetic changes in CBF-leukemia, we performed paired-end whole genome sequencing (WGS) on diagnostic leukemia blasts and matched germ line samples from 17 pediatric CBF-leukemia patients using the Illumina platform. Somatic alterations, including single nucleotide variations (SNVs) and structural variations (SVs), including insertions, deletions, inversions, and inter- and intra-chromosomal rearrangements, were detected using complementary analysis pipelines (Bambino, CREST and CONSERTING). Recurrent screening of identified mutations will be performed in a cohort of approximately 94 cases of CBF-leukemias.

  Study EGAS00001000349

• Long-term organoid culture of a small intestinal neuroendocrine tumor rna-seq

  Gastroenteropancreatic neuroendocrine tumors (GEP-NETs) are rare and highly heterogeneous neoplasms whose incidence has markedly increased over the last decades. A grading system based on the tumor cells' proliferation index predicts high-risk for G3 NETs. However, low-to-intermediate grade (G1/G2) NETs have an unpredictable clinical course that varies from indolent to highly malignant. Cultures of human cancer cells enable to perform functional perturbation analyses that are instrumental to enhance our understanding of cancer biology. To date, no tractable and reliable long-term culture of G1/G2 NET has been reported to permit disease modeling and pharmacological screens. Here, we report of the first long-term culture of a G2 metastatic small intestinal NET that preserves the main genetic drivers of the tumor and retains expression patterns of the endocrine cell lineage. Replicating the tissue, this long-term culture showed a low proliferation index, and yet it could be propagated continuously without dramatic changes in the karyotype. The model was readily available for pharmacological screens using targeted agents and as expected, showed low tumorigenic capacity in vivo. Overall, this is the first long-term culture of NETs to faithfully recapitulate many aspects of the original neuroendocrine tumor.

  Study EGAS00001007108

• ERG ALTERATIONS DEFINE A NOVEL SUBTYPE OF ACUTE LYMPHOBLASTIC LEUKEMIA

  Structural chromosomal alterations are a hallmark of acute lymphoblastic leukemia (ALL) yet many patients lack an identifiable abnormality on conventional cytogenetic analysis. Here, using integrated analysis of genome-wide DNA copy number and gene expression data of 1764 childhood and adult ALL cases, including whole-genome, exome and mRNA-sequencing of 72 cases, we report a novel subtype of B-progenitor ALL (B-ALL) representing 4-13% of B-ALL cases characterized by a distinct gene expression profile, and focal deletions of ERG (ETS related gene) in 56% of cases. These cases are characterized by expression of a ERG transcript that utilizes a novel exon in intron 6 of ERG spliced to the canonical reading frame of ERG exons 7-10, encoding a 27 kDa C-terminal ERG protein that retains the ETS and transactivating domains but lacks the N-terminal pointed and regulatory domains. This protein is a competitive inhibitor of wild-type ERG and is leukemogenic in vivo. ERG-altered ALL cases have a favourable outcome despite the presence of alterations associated with poor prognosis in non-ERG ALL, such as deletions of IKZF1. These findings identify a new subtype of leukemia characterized by a novel mechanism of ETS transcription factor dysregulation in cancer.

  Study EGAS00001000514

• Rare SNPs in receptor tyrosine kinases are negative outcome predictors in multiple myeloma

  Multiple myeloma (MM) is a plasma cell disorder that is characterized by a great genetic heterogeneity. Recent next generation sequencing studies revealed an accumulation of tumor-associated mutations in receptor tyrosine kinases (RTKs) which may also contribute to the activation of survival pathways in MM. To investigate the clinical role of RTK-mutations in MM, we deep-sequenced the coding DNA-sequence of EGFR, EPHA2, ERBB3, IGF1R, NTRK1 and NTRK2 which were previously found to be mutated in MM, in 75 uniformly treated MM patients of the "Deutsche Studiengruppe Multiples Myelom". Subsequently, we correlated the detected mutations with common cytogenetic alterations and clinical parameters. We identified 11 novel non-synonymous SNVs or rare patient-specific SNPs, not listed in the SNP databases 1000 genomes and dbSNP, in 10 primary MM cases. The mutations predominantly affected the tyrosine-kinase and ligand-binding domains and no correlation with cytogenetic parameters was found. Interestingly, however, patients with RTK-mutations, specifically those with rare patient-specific SNPs, showed a significant lower overall, event-free and progression-free survival. This indicates that RTK SNVs and rare patient-specific RTK SNPs are of prognostic relevance and suggests that MM patients with RTK-mutations could potentially profit from treatment with RTK-inhibitors.

  Study EGAS00001001665

• Cardiovascular disease biomarkers derived from circulating cell-free DNA methylation

  Acute coronary syndrome (ACS) remains a major cause of worldwide mortality. The syndrome occurs when blood flow to the heart muscle is decreased or blocked, causing muscle tissues to die or malfunction. There are three main types of ACS: Non-ST-elevation myocardial infarction, ST-elevation myocardial infarction, and unstable angina. The treatment depends on the type of ACS, and this is decided by a combination of clinical findings, such as electrocardiogram and plasma biomarkers. Circulating cell-free DNA (ccfDNA) is proposed as an additional marker for ACS since the damaged tissues can release DNA to the bloodstream. We used ccfDNA methylation profiles for differentiating between the ACS types and provided computational tools to repeat similar analysis for other diseases. We leveraged cell type specificity of DNA methylation to deconvolute the ccfDNA cell types of origin and to find methylation-based biomarkers that stratify patients. We identified hundreds of methylation markers associated with ACS types and validated them in an independent cohort. Many such markers were associated with genes involved in cardiovascular conditions and inflammation. ccfDNA methylation showed promise as a non-invasive diagnostic for acute coronary events. These methods are not limited to acute events, and may be used for chronic cardiovascular diseases as well.

  Study EGAS00001007263

• Developmental_Dysplasia_of_the_Hip__DDH_

  The National Joint Registry for England and Wales (NJR) has become established as one of the largest repositories of clinical information about patients undergoing hip and knee replacement. The primary aim of this proposal is to examine the feasibility of adding value to that clinical information by partnering it with a DNA dataset to allow its use to be extended to the study of the underlying pathogenesis of diseases of the hip and knee, and the complications of their treatment. We will do this 1) by establishing a mechanism for accessing and linking pseudo-anonymised demographic and phenotype data between the NJR database and the proposed DNA Biobank; 2) by validating the collection of postal DNA samples from consenting NJR patients; and 3) by genotyping genome-wide 900 patients with the condition developmental dysplasia of the hip (DDH) to find variants that play a role in this heritable condition. The successful demonstration that the NJR can be used to establish a DNA Biobank will provide a powerful platform to examine the genetic basis for a wide range of joint diseases, and the complications of their treatment. This proposal is to carry out the genotyping portion of this project on the current CoreExome array

  Study EGAS00001000916

• Single_Cell_Sequencing_of_Sperm__scSperm_

  The aim of this project is to genotype and sequence single spermatozoa from two men, one in his twenties and the other in his seventies. The resulting data is used to quantify the mutations that have arisen in the gametes of both individuals in order to better understand the effect of aging on mutation rates and modes.Project Outline. In order to quantify mutations, semen from two individuals are sequenced. 48 single sperm cells are isolated from each individual, and their DNA is extracted. The resulting genomes are amplified using PicoPlex, GenomiPhi MDA, Repli-G MDA, and MALBAC. QC step is applied to check the quality of WGA DNA using standard Sequenom plex (26 SNPs). A subset of 32 amplification products which pass the intiall QC, are genotyped using Affymetrix SNP6 chips. 12 of the genotyped amplification products are also sequenced. In addition, one multi-cell sample per individual is sequenced as a reference and for validation purposes.Altogether, 12 single cell sperm genomes and two multi-cell genomes are sequenced, coming to a total of 14 genomes. Of the single cell sperm genomes, 2 are sequenced to 50x coverage, and the other 10 to 25x coverage. Both multi-cell genomes are sequenced to 25x coverage.

  Study EGAS00001000935

• An RCOR1 loss-associated gene expression signature identifies a prognostically significant DLBCL subgroup

  Effective treatment of diffuse large B cell lymphoma (DLBCL) is plagued by heterogeneous responses to standard therapy, and molecular mechanisms underlying unfavorable outcomes in lymphoma patients remain elusive. Here, we profiled 148 genomes with 91 matching transcriptomes in an R-CHOP-treated DLBCL cohort to uncover molecular subgroups linked to treatment failure. Systematic integration of high-resolution genotyping arrays and RNA-seq data revealed novel deletions in RCOR1 to be associated with unfavorable progression free survival (P = 0.001). Integration of expression data from the clinical samples with data from RCOR1 knockdowns in the lymphoma cell lines KM-H2 and Raji yielded an RCOR1 loss-associated gene signature comprising of 233 genes. This signature identified a subgroup of patients with unfavorable overall survival (OS) (P = 0.023). The prognostic significance of the 233-gene signature for OS was reproduced in an independent cohort comprising of 195 R- CHOP treated patients (P = 0.039). Additionally, we discovered that within the IPI low risk group, the gene signature provides additional prognostic value that was independent of the cell-of-origin phenotype. Taken together, we present a novel and reproducible molecular subgroup of DLBCL, impacting risk-stratification of R- CHOP treated DLBCL patients and revealing a new avenue for therapeutic intervention strategies.

  Study EGAS00001001000

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Colorectal

  he 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. 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 EGAS00001003774

• Rapid identification of somatic genome rearrangements as personalized biomarkers for blood-based cancer monitoring

  Increasing evidence shows the value of circulating tumour DNA (ctDNA) to detect cancer and monitor its progression. Somatic genomic structural variations (SVs) are promising personalized biomarkers for sensitive and specific detection of ctDNA in liquid biopsies. However, accurate, affordable, and fast identification of such SV biomarkers is challenging, which hinders routine use in the clinic. Here, we demonstrate a novel approach - termed SHARC - for rapid discovery of somatic SV breakpoints as personalized tumour biomarkers. SHARC combines low-coverage cancer genome sketching by using Oxford Nanopore portable sequencing with a random forest classification and a dedicated filtering pipeline to enrich for somatic SVs. Our method leverages the real-time and long-read capabilities of Nanopore sequencing to identify somatic SV breakpoints at nucleotide resolution from a tumour biopsy within two days. We applied SHARC to tumour samples of high-grade ovarian and prostate cancer and validated on average 10 somatic SVs per sample with the use of PCR mini-amplicons. Finally, we demonstrate that these somatic SV biomarkers can be used to detect tumour presence from liquid biopsies in a quantitative manner and we retrospectively monitored treatment response in patients with prostate cancer, demonstrating its potential benefit for clinical practice.

  Study EGAS00001003963

• Genome-wide DNA-methylation assessment by MethylCap-seq and Infinium HumanMethylation450 BeadChips: an independent large-scale comparison

  DNA-methylation is an important epigenetic feature in health and disease. Two cost-efficient genome-scale methodologies to assess DNA-methylation are MethylCap-seq and Illumina's Infinium HumanMethylation450 BeadChips (HM450). However, objective information regarding the best-suited methodology for a specific research question is scant. Therefore, we performed a large-scale evaluation on a set of 70 brain tissue samples obtained from 65 glioblastoma and 5 non-tumoral brain tissues, using a gold standard free Bayesian modeling procedure. While conditional specificity was adequate for both approaches, conditional sensitivity was systematically higher for HM450. Also the genome-wide characteristics were compared, revealing that the HM450 probes assess less than 10% of the regions identified as methylated by MethylCap-seq. Hence the latter method may detect more potentially relevant DNA-methylation, defined by either functional location or previously reported differentially methylated candidate regions. Our results therefore indicate that – at least for the tissue under study - both methodologies are complementary, with a higher sensitivity for HM450, but a far larger genome-wide coverage for MethylCap-seq. Note that here only the relevant MethylCap-seq data is deposited, for the HM450 data we refer to GEO (GSE60274).

  Study EGAS00001001191

• Genomics of enteropathy associated T cell lymphoma (EATL)

  Enteropathy-associated T cell lymphoma (EATL) is a lethal and the most common neoplastic complication of celiac disease. In this study, we applied whole exome sequencing in 69 EATL tumors to define the genetic landscape of EATLs. The most frequently mutated gene in EATL was SETD2 (32% of cases), with loss of function mutations comprising a large proportion (41%) of those mutations. The JAK-STAT pathway was the most frequently mutated signaling pathway, with mutations in STAT5B being the most common, followed by JAK1, JAK3, STAT3, and the negative regulator of the pathway, SOCS1. We also identified mutations in NRAS and KRAS, as well as TP53 and TERT. Interestingly, the two types of EATL had highly overlapping genetic alterations indicating shared mechanisms underlying their pathogenesis. Gene expression analysis of the two types of EATL identify robust differences in transcriptional profiles. We further modeled the effects of SETD2 loss in vivo by developing a T cell specific knockout mouse model. These mice manifested an expansion of γδ T cells that identify novel roles for SETD2 in T cell development and lymphomagenesis. Together, our data render the most comprehensive genetic portrait yet of this uncommon but lethal disease, and may inform future classification schemes.

  Study EGAS00001001954

• Comparison of structural variations from 10X Genomics linked-reads and conventional Illumina short-reads sequencing

  Structural variations (SVs) are large genomic rearrangements that can drive many diseases. Conventional short-reads whole genome sequencing (cWGS) allows their identification with base-pair resolution, but suffers from high false discovery rate. cWGS taps in short-range information from short-reads while linked-reads sequencing (10XWGS) utilizes long-range information. 10XWGS allows linkage of short-reads originating from the same large DNA molecule with a unique barcode captured in a gel bead in emulsion. This mitigates alignment-based artefacts from cWGS especially in repetitive regions. However, the false discovery rate of this technology is unclear. In this study, we performed a comprehensive analysis of different type and size of SVs predicted from these two technologies. The SVs common between both technologies were found to be highly specific by PCR and Sanger sequencing while validation rate dropped for uncommon events. Further, we propose a novel enrichment approach for filtering out false positive calls from both the technologies independently. To this end, we trained a machine learning model for respective technologies and used it to characterise SVs from MCF7 cell line and a primary breast cancer tumor with high precision. This approach would be valuable in understanding true mechanisms driven by SVs in various diseases.

  Study EGAS00001004093

• Native American gene flow into Polynesia predating Easter Island settlement

  The possibility of voyaging contact between prehistoric Polynesian and Native American populations has long intrigued researchers. Proponents have pointed to the existence of New World crops, such as the sweet potato and bottle gourd, in the Polynesian archaeological record, but nowhere else outside the pre-Columbian Americas while critics have argued that these botanical dispersals need not have been human mediated. The Norwegian explorer Thor Heyerdahl controversially suggested that prehistoric South American populations had an important role in the settlement of east Polynesia and particularly of Easter Island (Rapa Nui). Several limited molecular genetic studies have reached opposing conclusions, and the possibility continues to be as hotly contested today as it was when first suggested. Here we analyse genome-wide variation in individuals from islands across Polynesia for signs of Native American admixture, analysing 807 individuals from 17 island populations and 15 Pacific coast Native American groups. We find conclusive evidence for prehistoric contact of Polynesian individuals with Native American individuals (around AD 1200) contemporaneous with the settlement of remote Oceania. Our analyses suggest strongly that a single contact event occurred in eastern Polynesia, before the settlement of Rapa Nui, between Polynesian individuals and a Native American group most closely related to the indigenous inhabitants of present-day Colombia.

  Study EGAS00001004209

• Pancreatic tropism of metastatic renal cell carcinoma

  Renal cell carcinoma (RCC) is characterized by a particularly broad metastatic swath, and enigmatically, when the pancreas is a destination, the disease is associated with improved survival. Intrigued by this observation, we sought to characterize the clinical behavior, therapeutic implications, and underlying biology. While pancreatic metastases (PM) are infrequent, we identified 31 patients across two institutional cohorts and show that improved survival is independent of established prognostic variables, that these tumors are exquisitely sensitive to anti-angiogenic agents and resistant to immune checkpoint inhibitors (ICI), and that they are characterized by a distinctive biology. Primary tumors of patients with PM exhibited frequent PBRM1 mutations, 3p loss, and 5q amplification, along with a lower frequency of aggressive features such as BAP1 mutations and loss of 9p, 14q, and 4q. Gene expression analyses revealed constrained evolution with remarkable uniformity, reduced T-effector gene signatures and increased angiogenesis. Similar findings were observed histopathologically. Thus, RCC metastatic to the pancreas is characterized by indolent biology, heightened angiogenesis, and an uninflamed stroma, likely underlying its good prognosis, sensitivity to anti-angiogenic therapies, and refractoriness to ICI. These data suggest that metastatic organotropism may be an indicator of a particular biology with prognostic and treatment implications for patients.

  Study EGAS00001004208

• Pediatric Non-Down Syndrome Acute Megakaryoblastic Leukemia is Characterized by Distinct Genomic Subsets with Varying Outcomes

  Acute Megakaryoblastic Leukemia (AMKL) is a subtype of acute myeloid leukemia (AML) in which cells morphologically resemble abnormal megakaryoblasts.  While extremely rare in adults, AMKL accounts for 4-15% of newly diagnosed childhood AML1-3.  A significant proportion of pediatric AMKL cases occur in children with Down syndrome (DS).  These patients have excellent outcomes and are characterized at the genomic level by a founding GATA1 mutation4-6.  In contrast, AMKL in patients without DS (non-DS-AMKL) is frequently associated with poor outcomes.  Previous efforts have identified chimeric oncogenes in a significant number of cases, including RBM15-MKL1, CBFA2T3-GLIS2, KMT2A gene rearrangements and NUP98-KDM5A7-9.  The etiology of 30-40% of cases, however, remains unknown.  To better understand the genomic landscape of non-DS-AMKL, we performed RNA and exome sequencing on specimens from 99 patients (75 pediatric and 24 adult).  We demonstrate that pediatric non-DS-AMKL is a heterogeneous malignancy that can be divided into seven subgroups with varying outcomes.  These subgroups are characterized by chimeric oncogenes with cooperating mutations in epigenetic and kinase signaling genes.  Overall, these data shed light on the etiology of AMKL and provide useful information for treatment tailoring. 

  Study EGAS00001002183

• CancerLocator: Non-Invasive Cancer Diagnosis and Tissue-of-Origin Prediction Using Methylation Profiles of Cell-Free DNA

  Background: The detection and characterization of cell-free DNA in plasma is one of the most promising new areas in cancer diagnosis. Liquid biopsy, unlike traditional tissue biopsy, has the potential to diagnose a variety of different malignancies. Results: Here we propose a probabilistic method, CancerLocator, which exploits the diagnostic potential of cell-free DNA by determining not only the presence but also the location of tumors. CancerLocator simultaneously infers the proportions and the tissue-of-origin of tumor-derived cell-free DNA in a blood sample using genome-wide DNA methylation data. We comprehensively evaluate CancerLocator with simulations and real data, and compare its performance with that of two established multi-class classification methods. We show that the predicted tumor burdens are highly consistent with the true values. In addition, when the proportion of tumor-derived DNAs in the cell-free DNAs is low, the two popular machine learning methods completely fail for cancer diagnosis, while CancerLocator successfully overcomes the challenge. CancerLocator also achieves promising results on patient plasma samples, despite the fact that the DNA methylation data from these samples has very low sequencing coverage. Conclusions: CfDNA methylation may be developed as an important approach for non-invasive early cancer diagnosis.

  Study EGAS00001002211

• Whole genome sequencing of tumour and normal paired samples of diffuse intrinsic pontine gliomas

  Diffuse intrinsic pontine glioma (DIPG) have a universally dismal prognosis (median 9-12 months), with neither chemotherapeutic nor targeted agents showing any substantial survival benefit in clinical trials in children with these tumours. DIPG are highly infiltrative malignant glial neoplasms of the ventral pons, which due to their location within the brain, make them unsuitable for surgical resection and consequently have a universally dismal clinical outcome. Recent high-throughput sequencing approaches have revealed a striking prevalence of K27M mutations in the genes encoding the histone variants H3.3 (H3F3A) or H3.1 (HIST1H3B) in the childhood brain tumour DIPG. This K-to-M substitution confers a trans-dominant ablation of global H3K27 trimethylation, which likely profoundly alters gene expression through de-repression of polycomb repressive complex 2 (PRC2) target genes. Despite these advances in our understanding of the distinct biology of these tumours, approaches for specific novel therapeutic interventions are not clear, and little has been reported of the secondary mutations accompanying these changes. We carried out whole genome sequencing on a series of 20 tumour normal pairs of pre-treatment biopsy samples of DIPG, for which the patients underwent a safe stereotactic procedure, and whole exome sequencing on a further six tumour normal pairs obtained at autopsy.

  Study EGAS00001000572

• Accurate sample assignment in a multiplexed, ultra-sensitive, high-throughput sequencing assay for minimal residual disease

  High throughput sequencing (HTS) (next generation sequencing) of the rearranged immunoglobulin and T-cell receptor genes promises to be cheaper and more sensitive than current methods for monitoring minimal residual disease (MRD) in patients with acute lymphoblastic leukemia. However, adoption of new approaches by clinical laboratories requires careful evaluation of all potential sources of error and the development of strategies to ensure the highest accuracy. Timely and efficient clinical use of HTS platforms will depend on combining multiple samples (multiplexing) in each sequencing run. Here we examine immunoglobulin heavy chain gene HTS on the Illumina MiSeq platform for MRD (HTS-MRD). We identify errors associated with multiplexing that could potentially impact on the accuracy of MRD analysis. We optimise a strategy combining high purity, sequence-optimised oligonucleotides, dual-indexing and an error-aware demultiplexing approach to minimise errors and maximise sensitivity. We present a probability-based demultiplexing pipeline, Error-Aware Demultiplexer (EAD) - that is suitable for all MiSeq sequencing strategies and accurately assigns samples to the correct identifier without excessive loss of data. Finally using controls quantified by digital PCR, we show that HTS-MRD can accurately detect as few as 1 in 106 copies of specific leukemic MRD.

  Study EGAS00001001303

• Biallelic mutations in the ubiquitin ligase RFWD3 cause Fanconi anemia

  The WD40-containing E3 ubiquitin ligase RFWD3 has been recently linked to the repair of DNA damage by Homologous Recombination (HR). Here we show that an RFWD3 mutation within the WD40 domain is connected to the genetic disease Fanconi anemia (FA). An individual revealed congenital abnormalities characteristic for FA. Cells from the patient, carrying the compound heterozygous mutations c.205_206dupCC and c.1916T>A in RFWD3, show increased sensitivity to DNA interstrand cross-linking agents in terms of increased chromosomal breakage, reduced survival and cell cycle arrest in G2 phase. The cellular phenotype is mirrored in genetically engineered human and avian cells by inactivation of RFWD3 or by introduction of the patient-derived missense mutation, and is rescued by expression of wildtype RFWD3 protein. HR is disrupted in RFWD3 mutant cells, caused by impaired relocation of mutant RFWD3 to chromatin and defective physical interaction with RPA. Rfwd3 knockout mice exhibit increased embryonic lethality, are sub-fertile, show ovarian and testicular atrophy and have a reduced life span in which they resemble other FA mouse models. Although mutation was detected in a single child with FA yet, we propose RFWD3 as a novel FA gene, FANCW, supported by cellular paradigm systems and an animal model.

  Study EGAS00001002440

• High Grade Serous Ovarian Carcinomas Originate in the Fallopian Tube

  High-grade serous ovarian carcinoma (HGSOC) is the most frequent type of ovarian cancer and has a poor outcome. It has been proposed that fallopian tube cancers may be precursors of HGSOC but evolutionary evidence for this hypothesis has been limited. To provide insights into the origins of HGSOC, we We performed whole-exome sequence and copy number analyses of laser-capture microdissected fallopian tube lesions (p53 signatures, serous tubal intraepithelial carcinomas (STICs), and fallopian tube carcinomas), ovarian cancers, and metastases from nine patients with HGSOC. The majority of tumor-specific sequence and structural alterations in ovarian cancers were present in STICs, including those affecting TP53, BRCA1, BRCA2 or PTEN genes. An eEvolutionary analyseis revealed that p53 signatures and STICs were are the precursors of ovarian carcinoma which in turn gave rise to metastatic lesionsand identify a. In one patient we identified a second STIC as a metastasis in the fallopian tube opposite from the affected ovary. These analyses revealed a window of seven years between the development of a STIC and the initiation of ovarian carcinomas, with development of metastases following rapidly thereafter. Our results provide insights into the etiology of ovarian cancer and have implications for the prevention, early detection and therapeutic intervention of this disease.

  Study EGAS00001002589

• Whole-exome and targeted sequencing of pediatric hyperdiploid B-cell precursor acute lymphoblastic leukemia

  High hyperdiploidy defines the largest genetic entity of childhood acute lymphoblastic leukemia (ALL). Despite its relatively low recurrence risk, this subgroup generates a high proportion of relapses. The cause and origin of these relapses remains obscure. We therefore explored the mutational landscape in high hyperdiploid (HD) ALL with whole-exome (n=19) and subsequent targeted deep sequencing of 60 genes in 100 relapsing and 51 non-relapsing cases. We identified multiple clones at diagnosis that were primarily defined by a variety of mutations in receptor tyrosine kinase (RTK)/Ras pathway and chromatin-modifying genes. The relapse clones consisted of reappearing as well as new mutations, and overall contained more mutations. Although RTK/Ras pathway mutations were similarly frequent between diagnosis and relapse, both intergenic and intragenic heterogeneity was essentially lost at relapse. CREBBP mutations, however, increased from initially 18-30% at relapse, then commonly co-occurred with KRAS mutations (P<0.001) and these relapses appeared primarily early (P=0.012). Our results confirm the exceptional susceptibility of HD ALL to RTK/Ras pathway and CREBBP mutations, but, more importantly, suggest that mutant KRAS and CREBBP might cooperate and equip cells with the necessary capacity to evolve into a relapse-generating clone.

  Study EGAS00001001113

• Single cell transcriptomic landscape of pediatric B-cell acute lymphoblastic leukemia: dissection of transcriptional heterogeneity and B-cell developmental state

  Introduction: Genomic analyses of bulk ALL samples have improved our understanding of the genetic basis and risk stratification of B-ALL, but do not directly examine intratumor heterogeneity or enable inference of leukemia developmental state and cell of origin. Methods: We profiled 89 B-ALL samples by single-cell RNA-seq (scRNA-seq) (10X Genomics 5’v2) and compared them to a scRNA-seq reference map of normal human B-cell development. Results: Analysis of heterogeneity of inferred DNA copy number alterations at single cell level showed that aneuploid ALL with near haploidy and hyperdiploid harbor chromosomal losses or gains in all blasts, consistent with early, synchronous origin rather than sequential accumulation. Intra-sample heterogeneity was driven by five gene expression signatures: cell cycle, metabolism, differentiation, inflammation, refining leukemic subtyping. An atlas of B cell development was established to map B-ALL samples to specific developmental states. Hyperdiploid samples resembled pro-B cells, while MEF2D-r and TCF3::PBX1 showed high pre-B abundance. Patients with high early lymphoid abundance had poorer outcomes. Conclusions: Understanding variation in transcriptional programs and developmental states of B-ALL blasts by sc transcriptome refines existing clinical and genomic classifications and provides novel prognostic markers.

  Study EGAS00001007512

• Whole-Genome Sequencing Suggests Schizophrenia Risk Mechanisms in Humans with 22q11.2 Deletion Syndrome

  Chromosome 22q11.2 microdeletions impart a high but incomplete risk for schizophrenia. Possible mechanisms include genome-wide effects of DGCR8 haploinsufficiency. In a proof-of-principle study to assess the power of this model, we used high-quality, whole-genome sequencing of nine individuals with 22q11.2 deletions and extreme phenotypes (schizophrenia, or no psychotic disorder at age >50 years). The schizophrenia group had a greater burden of rare, damaging variants impacting protein-coding neurofunctional genes, including genes involved in neuron projection (nominal P = 0.02, joint burden of three variant types). Variants in the intact 22q11.2 region were not major contributors. Restricting to genes affected by a DGCR8 mechanism tended to amplify between-group differences. Damaging variants in highly conserved long intergenic noncoding RNA genes also were enriched in the schizophrenia group (nominal P = 0.04). The findings support the 22q11.2 deletion model as a threshold-lowering first hit for schizophrenia risk. If applied to a larger and thus better-powered cohort, this appears to be a promising approach to identify genome-wide rare variants in coding and noncoding sequence that perturb gene networks relevant to idiopathic schizophrenia. Similarly designed studies exploiting genetic models may prove useful to help delineate the genetic architecture of other complex phenotypes.

  Study EGAS00001002344

• Novel genes for Intellectual Disability identified using whole genome sequence and pathway analysis

  Intellectual Disability (ID) is one of the most common global disorders. However for ~30% of patients with ID no cause is identified, despite a clinical diagnostic odyssey that includes genome wide clinical microarray (CMA). We carefully selected 8 trios, each composed of a child with ID and brain structural defect, and both normal parents for whole genome sequencing (WGS). We conducted WGS on the trio and filtered out de novo single nucleotide variants (SNVs) and then used a pathway-based refinement to select candidates. We confirmed a de novo pathogenic SNV in ARID1B, and likely pathogenic SNVs in CACNB3, SPRY4, PHF6, SQSTM1 and UPF1 in 5 of the 8 children. All genes except ARID1B and PHF6 are previously unreported for ID. We analysed our WGS data using 4 independent algorthims for copy number and structural variation including using de novo whole genome assembly. We confirmed one likely contributory 165kb de novo CNV (missed by CMA). We conducted a validation study of over 2000 exomes for our novel candidate genes and also present a strategy to analyze the non-coding sequence space. This study provides an extensive analysis of WGS in the context of ID and yielded causative variants in >60% of cases.

  Study EGAS00001001386

• Ewing and Ewing-like-sarcoma tumoroid biobank reveals distinct drug sensitivities: translocation makes the difference

  Ewing sarcoma (ES) and Ewing-like sarcomas (ELS) are highly aggressive tumors found in soft tissues and bone of mostly children and young adults. Despite being different in many aspects, including genetics, possible cell-of-origin, and pathology, patients with any of these entities all receive the same therapeutic regimen. Currently, there is a lack of preclinical cell models that can predict drug response of these entities in vitro. Here, we describe a pediatric ES and ELS tumor organoid (tumoroid) biobank containing tumoroids with different translocations, including EWSR1::FLI1, EWSR1::ERG, CIC::DUX4, and BCOR-rearrangements. Using histology, whole genome sequencing and RNA sequencing, we demonstrate that these tumoroids retain key properties of their matching patient tumors. In addition, we traced tumor subclones in the tumoroids across patient-matched longitudinal samples, which showed that cellular heterogeneity was maintained. We subsequently performed drug screening on the tumoroid models, which unveiled entity-specific drug sensitivity to various cytotoxic compounds and targeted compounds, including MCL-1 inhibitors for CIC::DUX4 sarcomas. Taken together, this newly established ES and ELS patient-derived tumoroid biobank generates a unique source of material for future basic cancer research and drug-screening.

  Study EGAS00001007941

• Pooled scRNA-seq of iPSC-derived neural stem cells from ADHD and control individuals

  This dataset comprises paired-end FASTQ files from single-cell RNA sequencing (scRNA-seq) of neural stem cells (NSCs) derived from induced pluripotent stem cells (iPSCs). The iPSCs were reprogrammed from somatic cells (keratinocytes, PBMCs, or lymphoblastoid cell lines) collected from individuals diagnosed with attention-deficit hyperactivity disorder (ADHD) and age-matched healthy controls (ages 6–18 years). Reprogramming was carried out using the non-integrative Sendai virus system (CytoTune-iPS 2.0, Thermo Fisher Scientific), introducing OCT3/4, SOX2, KLF4, and cMYC. iPSCs were cultured in Essential 8 medium on vitronectin and underwent rigorous quality control including immunocytochemistry, qRT-PCR, mycoplasma testing, and SNP genotyping (Illumina Global Screening Array). NSCs were generated using a standardized neural induction protocol and validated via qRT-PCR and immunostaining for markers including PAX6, SOX2, NESTIN, TUJ1, and FOXG1. Two independent scRNA-seq multiplexed experiments were performed: Experiment 1 (ADHD: n = 7; Control: n = 6) and Experiment 2 (ADHD: n = 9; Control: n = 7). Samples were pooled and demultiplexed using the Ensemblex algorithm, incorporating prior genotype data. After quality control filtering (>500 unique RNA molecules per cell), 30,433 single-cell transcriptomes were retained. Differential gene expression analysis was conducted to identify transcriptional signatures associated with ADHD compared to controls.

  Study EGAS00001008169

• Mapping_gene_environment_interactions_in_macrophages

  Immune cells sense and respond to external stimuli such as pathogens and initiate an inflammatory response. Genetic variants altering the behaviour and responses of immune cells may inform the results of disease association studies and provide an insight into the regulation of immune response. Studies mapping quantitative trait loci (QTL) are usually performed in quiescent or naïve cells and therefore any variants with a condition-specific manner (e.g. those with an effect only detectable in a specific time-point after cell stimulation) are likely to remain undiscovered. So far, only a handful of studies have attempted to study QTL in multiple states and even on those studies the number of states has remained relatively low (<5). Only very recently Horst and colleagues published the largest study to date investigating the effect of host variable (e.g. age, sex, contraceptive use etc.) on protein levels in macrophages stimulated with 19 different conditions (Horst et al, 2016). This project aims to develop a high-throughput experimental platform to detect QTL in a large number of states using iPSC derived macrophages. 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 EGAS00001002268

• Genomic profiling of Acute Lymphoblastic Leukemia in Ataxia Telangiectasia patients reveals tight link between ATM mutations and chromothripsis

  Recent developments in sequencing technologies led to the discovery of a novel form of genomic instability, termed chromothripsis. This catastrophic genomic event, involved in tumorigenesis, is characterized by tens to hundreds of simultaneously acquired locally clustered rearrangements on one chromosome. We hypothesized that leukemias developing in individuals with Ataxia Telangiectasia, who are born with two mutated copies of the ATM gene, an essential guardian of genome stability, would show a higher prevalence of chromothripsis due to the associated defect in DNA double-strand break repair. Using whole-genome sequencing, fluorescence in situ hybridization and RNA sequencing, we characterized the genomic landscape of Acute Lymphoblastic Leukemia (ALL) arising in patients with Ataxia Telangiectasia. We detected a high frequency of chromothriptic events in these tumors, specifically on acrocentric chromosomes, as compared to tumors from individuals with other types of DNA repair syndromes (27 cases total, 10 with Ataxia Telangiectasia). Our data suggest that the genomic landscape of Ataxia Telangiectasia ALL is clearly distinct from that of sporadic ALL. Mechanistically, short telomeres and compromised DNA damage response in cells of Ataxia Telangiectasia patients may be linked with frequent chromothripsis. Furthermore, we show that ATM loss is associated with increased chromothripsis prevalence in additional tumor entities.

  Study EGAS00001002270

• Whole genome, whole exome, and targeted sequencing of high-grade meningioma tumor samples.

  High-grade meningiomas frequently recur and are associated with high rates of morbidity and mortality. To determine the factors that promote the development and evolution of these tumors, we analyzed the genomes of 274 high-grade meningiomas and compared this information with data from 456 low-grade meningiomas. High-grade meningiomas had a higher mutation burden than low-grade meningiomas but did not harbor any statistically significant mutated genes aside from NF2. High-grade meningiomas also possessed significantly elevated rates of chromosomal gains and losses, especially among tumors with monosomy 22. Meningiomas previously treated with adjuvant radiation had significantly more copy number alterations than radiation-induced or radiation-naïve meningiomas. Across serial recurrences, genomic disruption preceded the emergence of nearly all mutations, remained largely uniform across time, and when present in low-grade meningiomas, correlated with subsequent progression to a higher grade. In contrast to the largely stable copy number alterations, mutations were strikingly heterogeneous across tumor recurrences, likely due to extensive geographic heterogeneity in the primary tumor. While high-grade meningiomas harbored significantly fewer overtly targetable alterations than low-grade meningiomas, they contained numerous mutations that are predicted to be neoantigens, suggesting that immunologic targeting may be of therapeutic value.

  Study EGAS00001002294

• Botswana 15 autosomal unlinked microsatellites

  Population structure was investigated in 990 Botswana individuals according to ethno-linguistics, Bantu and Khoisan, and geography (the nine administrative districts) using the AmpFLSTR Identifiler autosomal microsatellite markers. Genetic diversity and forensic parameters were calculated for the overall population, and according to ethno-linguistics and geography. The overall combined power of exclusion (CPE) was 0.9999965412 and the combined match probability 6,28 x 10-19. CPE was highest for the Khoisan Tuu ethnolinguistic group and the Northeast District at 0.9999582029 and 0.9999922652 respectively. CMP ranged from 6.28 x 10-19 (Khoisan Tuu) to 1,02 x 10-18 (Northwest district). Using pairwise genetic distances (FST), analysis of molecular variance (AMOVA), factorial correspondence analysis (FCA), and the unsupervised Bayesian clustering method found in STRUCTURE and landscape genetics software TESS, ethno-linguistics were found to have a greater influence on population structure than geography. FCA showed clustering between Bantu and Khoisan, and within the Bantu. This Bantu sub-structuring was not seen with STRUCTURE and TESS, which detected clustering only between Bantu and Khoisan. The patterns of population structure revealed highlight the need for regional reference databases that include ethno-linguistic and geographic location information. These markers have important potential for bio-anthropological studies as well as for forensic applications.

  Study EGAS00001002380

• The evolutionary trajectory of a highly recurrent paediatric high grade neuroepithelial tumour with MN1:BEND2 fusion

  Astroblastomas are rare brain tumours which predominate in children and young adults, and have a controversial claim as a distinct entity, with no established WHO grade. Reports suggest a better outcome than high grade gliomas, though they frequently recur. Recently, they have been described to overlap with a newly-discovered group of tumours described as ’high grade neuroepithelial tumour with MN1 alteration’ (CNS HGNET-MN1), defined by global methylation patterns and strongly associated with gene fusions targeting MN1. We have studied a unique case of astroblastoma arising in a 6 year old girl, with multiple recurrences over a period of 10 years, with the pathognomonic MN1:BEND2 fusion. Exome sequencing allowed for a phylogenetic reconstruction of tumour evolution, which when integrated with clinical, pathological and radiological data provide for a thorough understanding of disease progression, with initial treatment driving tumour dissemination along four distinct trajectories. Infiltration of distant sites was associated with a later genome doubling, whilst there was evidence of convergent evolution of different lesions acquiring distinct alterations targeting NF-κB. These data represent an usual opportunity to understand the evolutionary history of a highly recurrent childhood brain tumour, and provide novel therapeutic targets for astroblastoma / CNS HGNET-MN1.

  Study EGAS00001002432

• Natural genetic variation of the cardiac transcriptome in non-diseased donors and patients with dilated cardiomyopathy

  Background: Genetic variation is an important determinant of RNA transcription and splicing, which in turn contributes to variation in human traits including cardiovascular diseases.Results: Here we report the first in-depth survey of heart transcriptome variation using RNA-sequencing in 149 (97)* patients with dilated cardiomyopathy and 113 (108)* non-diseased controls. We reveal extensive differences of gene expression and splicing between dilated cardiomyopathy patients and controls, affecting known as well as novel dilated cardiomyopathy genes. Moreover, we show a widespread effect of genetic variation on the regulation of transcription, isoform usage and allele specific expression. Systematic annotation of genome wide association SNPs identifies 60 functional candidate genes for heart phenotypes, representing 20% of all published heart genome wide association loci. Focusing on the dilated cardiomyopathy phenotype we found that eQTL variants are also enriched for dilated cardiomyopathy genome wide association signals in two independent cohorts.Conclusions: RNA transcription, splicing and allele specific expression are each important determinants of the DCM phenotype and are controlled by genetic factors. Our results represent a powerful resource for the field of cardiovascular genetics. (*) numbers in parentheses represent samples that the publication is based on.

  Study EGAS00001002454

• Direct Detection of Early-stage Cancers Using Circulating Tumor DNA

  Early detection and intervention are likely to be the most effective means for reducing morbidity and mortality of human cancer. However, development of methods for noninvasive detection of early-stage tumors has remained a challenge. We have developed an approach called targeted error correction sequencing (TEC-Seq) that allows ultrasensitive direct evaluation of sequence changes in circulating cell-free DNA using massively parallel sequencing. We have used this approach to examine 58 cancer-related genes encompassing 81 kb. Analysis of plasma from 44 healthy individuals identified genomic changes related to clonal hematopoiesis in 16% of asymptomatic individuals but no alterations in driver genes related to solid cancers. Evaluation of 200 patients with colorectal, breast, lung, or ovarian cancer detected somatic mutations in the plasma of 71, 59, 59, and 68%, respectively, of patients with stage I or II disease. Analyses of mutations in the circulation revealed high concordance with alterations in the tumors of these patients. In patients with resectable colorectal cancers, higher amounts of preoperative circulating tumor DNA were associated with disease recurrence and decreased overall survival. These analyses provide a broadly applicable approach for noninvasive detection of early-stage tumors that may be useful for screening and management of patients with cancer.

  Study EGAS00001002577

• Inhibiton of the GABPB1L-containing GABP tetramer is sufficient to reverse replicative immortality in TERT promoter mutant glioblastoma cells.

  Over 50 types of cancer acquire TERT promoter mutations. These single point mutations reactivate telomerase, allowing for indefinite maintenance of telomere length and enabling cellular immortalization. The transcription factor binding site created by the point mutations specifically recruit the ETS factor GABP, a multimeric transcription factor composed of the GABPα and GABPβ subunits. GABP can form two functionally independent transcription factor species – a dimer or a tetramer – depending on which of the structurally distinct GABPβ isoforms is incorporated into the complex. In this study, we show that genetic disruption of GABPβ1L, a tetramer forming isoform of GABPβ that is dispensable in normal development, results in TERT silencing in a TERT promoter mutation dependent manner. Failure to activate TERT expression by GABPβ1L culminates in telomere dysfunction, DNA damage, and mitotic cell death exclusively in TERT promoter mutant cells. Furthermore, exogenous expression of TERT is sufficient to prevent telomere degradation and loss of cell viability in GABPβ1L-reduced lines bearing the mutant TERT promoter. Orthotopic injection of tetramer-deficient mutant TERT promoter GBM cells rendered lower tumor burden and prolonged the overall survival of the tumor-bearing mice. These results highlight the potentially widespread role of GABPB1L in enabling replicative immortality of TERT promoter cancers.

  Study EGAS00001002582

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__clear_cell_renal_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. 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 EGAS00001002608

• WGS of 32 paired SRCC samples

  Signet-ring cell carcinoma (SRCC) has specific epidemiology and oncogenesis in gastric cancer, however with no systematical investigation for prognostic genomic features. We systematically investigated the clinical characteristics of 1,868 Chinese gastric cancer patients, and found that signet-ring cell content was significantly related to multiple clinical characteristics and treatment outcomes. Next, we performed whole genome sequencing for 32 SRCC patients with > 80% signet-ring cells in tumor, and identified frequent5 CLDN18-ARHGAP26/6 fusion (25%). With additional 797 patients for validation, we notice that prevalence of CLDN18-ARHGAP26/6 fusion is positively associated with signet-ring cell content (P = 4.1 x 10 -9 ), younger age at diagnosis (P = 4.2 x 10 -10 ), female/male ratio (P = 1.7 x 10 -9 ), and advanced TNM stage (P = 1.7 x 10 -5 ). Importantly, patients with CLDN18-ARHGAP26/6 fusion had a worse survival outcomes (P = 0.03), and got no benefit0 from platinum/fluoropyrimidines based chemotherapy (P = 0.92) compared to fusion-free patients (P = 0.001), which is consistent with the fact of oxaliplatin/fluoropyrimidine resistance acquired in CLDN18-ARHGAP26 introduced gastric cancer cell lines. Overall, this study provides new insights into the clinical and genomic features of SRCC, and highlights the importance of frequent CLDN18-ARHGAP26/6 fusions in chemotherapy response for SRCC.

  Study EGAS00001002668

• UK10K_RARE_NEUROMUSCULAR

  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 neuromuscular disorder samples are part of the “rare disease” group, and will undergo exome sequencing. For further information with regard to this cohort please contact Francesco Muntoni (f.muntoni@ucl.ac.uk).

  Study EGAS00001000101

• CancerDetector: Ultrasensitive and Non-Invasive Cancer Detection at the Resolution of Individual Reads using Cell-free DNA Methylation Sequencing Data

  The detection of tumor-derived cell-free DNA in plasma is one of the most promising directions in cancer diagnosis. The major challenge in such approach is how to identify the tiny amount of tumor DNAs out of total cell-free DNAs in blood. Here we propose an ultrasensitive cancer detection method, termed “CancerDetector”, using the DNA methylation profiles of cell-free DNAs. The key of our method is to probabilistically model the joint methylation patterns of multiple adjacent CpG sites on an individual sequencing read, in order to exploit the pervasive nature of DNA methylation for signal amplification. Therefore, CancerDetector can sensitively identify a trace amount of tumor cfDNAs in plasma, at the level of individual reads. We evaluated CancerDetector on the simulated data, and showed a high concordance of the predicted and true tumor burden. Testing CancerDetector on real plasma data demonstrated its high sensitivity and specificity in detecting tumor DNAs. In addition, the predicted tumor burden showed great consistency with tumor size and survival outcome. Note that all of those testing were performed on sequencing data at low to medium coverage (1X to 10X). Therefore, CancerDetector holds the great potential to detect cancer early and cost-effectively.

  Study EGAS00001002728

• UK10K_RARE_FIND

  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 Raymond samples will be part of the “rare disease” group, and will undergo exome sequencing. For further information with regard to this cohort please contact Lucy Raymond (flr24@cam.ac.uk).

  Study EGAS00001000128

• UK10K_RARE_CILIOPATHIES

  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 Ciliopathies samples will be part of the rare disease group, and will undergo exome sequencing. For further information with regard to this cohort please contact Phil Beales (p.beales@ich.ucl.ac.uk).

  Study EGAS00001000126

• UK10K_RARE_SIR

  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 SIR (Severe Insulin Resistance) samples are part of the “rare disease” group, and will undergo exome sequencing. For further information with regard to this cohort please contact Robert Semple (rks16@cam.ac.uk).

  Study EGAS00001000130

• Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell RCC -- CA209-009

  Immune checkpoint inhibitors targeting the PD-1 pathway have transformed the management of many advanced malignancies, including clear cell renal cell carcinoma (ccRCC), but the drivers and resistors of PD-1 response remain incompletely elucidated. Here, we analyzed 592 tumors collected from advanced ccRCC patients enrolled in prospective clinical trials of treatment with PD-1 blockade (or mTOR inhibition as control arm) by whole-exome and RNA-sequencing, integrated with immunofluorescence analysis, to define the somatic alteration landscape of late-stage ccRCC and to uncover the immunogenomic determinants of therapeutic response. While conventional genomic markers (tumor mutation burden, neoantigen load) and degree of CD8+ T cell infiltration were not associated with clinical response, we discovered numerous chromosomal alterations in advanced ccRCC associated with response or resistance to PD-1 blockade. These advanced tumors were highly CD8+ T cell infiltrated, with only 22% and 5% with an immune desert and immune excluded phenotype, respectively. Our analysis revealed that CD8+ infiltrated tumors are depleted of favorable PBRM1 mutations and are enriched for unfavorable chromosomal losses of 9p21.3 when compared to non-infiltrated tumors. These data demonstrate how the interplay of somatic alterations and immunophenotypes impacts therapeutic efficacy.

  Study EGAS00001004290

• Stratifying and Targeting Pediatric Medulloblastoma through Genomics

  In this project, genomic analyses of pediatric medulloblastoma samples, obtained through the international medulloblastoma consortium, will be performed. RNA and miRNA expression profiles of 1000 samples, representing all four subgroups (Wnt, Shh, Group C, and D), will be studied to identify novel subtypes within each subgroup. The resulting subtype-specific expression profiles will support the development of reliable and robust biomarkers to more accurately and reliably classify medulloblastomas for treatment in clinical trials. For that purpose, two assays will be developed: an antibody-based immunohistochemical assay and an orthogonal nucleic acid-based hybridization assay. Additional genomic DNA analysis of the 300 high risk subgroup cases will support the discovery of subgroup specific somatic mutations in order to inform current clinical trials of targeted therapies, and to identify genes and pathways already targeted in other diseases. Such therapies could be rapidly transitioned to Phase II trials in medulloblastoma. Furthermore, the discovery of somatic mutations could be used for developing as well as validating specific biomarkers. The project team will also try to identify risk factors that predispose children to this type of cancer.

  Study EGAS00001000273

• Clonal_haematopoiesis_in_patients_with_AAA

  Recent advances in genomics have demonstrated that clonal haemopoiesis driven by leukaemia associated somatic mutations is a relatively common phenomenon that increases in frequency with advancing age. Whilst individuals with clonal haemopoiesis have an increased risk of developing haematological malignancies, they also have an increased mortality from other causes. Additionally, certain mutations are almost exclusively seen in individuals aged 70 years or older, whilst others are seen in individuals with non-haematological cancers including breast and ovarian. Recently, clonal haemopoiesis was found to be associated with a significantly increased risk of atherosclerotic cardiovascular disease. This association is thought to be causative with clonally-derived macrophages showing elevated expression of several chemokine and cytokine genes that contribute to atherosclerosis. Another vascular pathology, abdominal aortic aneurysm (AAA), increases with age and shares risk factors with atherosclerosis (including smoking, male sex, high cholesterol). However, the impact of these risk factors and the overlap between AAA and atherosclerosis is poorly understood. To investigate a possible link between clonal haemopoiesis and AAA, we will study DNA samples from 300 patients with AAA and up to 200 controls for evidence of clonal haemopoiesis. This will be done using target DNA enrichment with biotinylated RNA baits followed by high throughput sequencing.

  Study EGAS00001002873

• UK10K RARE CHD

  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 CHD (congenital heart disease) samples will be part of the “rare disease” group, and will undergo exome sequencing.For further information with regard to this cohort please contact Shoumo Bhattacharya (shoumo@me.com).

  Study EGAS00001000125

• Genomic Alterations in Gingivo-buccal Cancer: ICGC-India Project_YR02

  As a part of the ICGC, India has undertaken genomic studies on gingivobuccal cancer of the oral cavity, which is the most prevalent form of cancer among men in India. There are various known environmental (life-style) correlates of this cancer, the most important of which are tobacco chewing and HPV infection. Paired DNA samples – isolated from the tumour tissue and from the blood of fifty six patients – have been analyzed to catalog germline and somatic mutations. Association between each observed genomic alteration and exposure to environmental risk factors is being explored.Detailed clinical characterization of the patients, collection of data on demographic and environmental exposures, and isolation of DNA samples from blood and tumour tissues collected from each patient are being done at the Advanced Centre for Research, Treatment and Education on Cancer, Mumbai. Exome capture and deep resequencing are being performed at the National Institute of Biomedical Genomics, Kalyani. We have analyzed the exomes of paired blood and tumour DNA samples using exome capture by Illumina TrueSeq Exome Enrichment kit followed by sequencing on Illumina HiSeq-2000. Each exome is being sequenced at a mean depth of 70x for Tumour and 40x for Blood. In addition, we have used Illumina Omni 2.5 SNP-chips to generate genotype data.

  Study EGAS00001001028

• P647_Targeted_resequencing_project

  In collaboration with Dr Robert Semple we have identified a family harbouring an autosomal dominant variant, which leads to severe insulin resistance (SIR), short stature and facial dysmorphism. This family is unique within the SIR cohort in having normal lipid profiles, preserved adiponectin and normal INSR expression and phosphorylation. DNA is available for 7 affected and 7 unaffected family members across 3 generations. All 14 samples have been genotyped using microsatellites and the Affymetrix 6.0 SNP chip. Linkage analysis identified an 18.8Mb haplotype on chromosome 19 as a possible location of the causative variant. However, Exome sequencing of 3 affected and 1 unaffected family members has not identified the causative variant suggesting the possibility of an intronic or intergenic variant in this region or elsewhere in the genome. We propose to conduct whole genome sequencing of 5 members of the pedigree at a depth of 20X. The chosen samples are two sets of parents plus one member of an unaffected branch of the pedigree who shares the risk haplotype on chromosome 19. Sequencing of the two sets of parents will be used along with the genome-wide SNP data to impute 4 affected children giving an effect sample size of 6 affected individuals.

  Study EGAS00001000305

• Mutational Patterns in Metastatic Cutaneous Squamous Cell Carcinoma

  Cutaneous squamous cell carcinoma (cSCC) from the head and neck typically metastasizes to the lymph nodes of the neck and parotid glands. When a primary is not identified, they are difficult to distinguish from metastases of mucosal origin and primary salivary gland SCC. Ultraviolet radiation causes a mutation pattern that predominantly features cytosine to thymine transitions at dipyrimidine sites and has been associated with cSCC. In this study, we used whole genome sequencing data from 15 cSCC metastases and show that a UV signature mutation is pervasive across the cohort and distinct from mucosal SCC. The mutational burden was exceptionally high and concentrated in some regions of the genome, especially insulator elements (mean 162 mutations / Mb). We therefore evaluated the likely impact of UV-induced mutations on the dipyrimidine rich binding site of the main human insulator protein, CCCTC-binding factor (CTCF), and the possible implications on CTCF function and the spatial organization of the genome. Our findings suggest that mutation signature analysis may be useful in determining the origin of metastases in the neck and the parotid gland. Furthermore, UV-induced DNA damage to insulator binding sites may play a role in the carcinogenesis and progression of cSCC.

  Study EGAS00001003370

• Exome sequencing of a novel cervical cancer cell line

  Human cancer cell lines are largely used in the searching for new antineoplastic agents. However, due to the artifacts of a long-term in culture, cell lines do not always represent the realistic tumor cell behavior. This has motivated the development of models that better mimetics the tumor tissue, among them, the establishment of primary cell cultures. In this work, we establish and characterized a low-passage cervix cancer cell line from a Brazilian patient with squamous cell carcinoma. The phenotype confirms the epithelial and tumor origin, through cytokeratins, EpCAM, and p16 staining. Whole exome sequencing showed relevant somatic mutations in several genes including BRCA2, TGFBR1, and IRX2 genes. CNV analysis by nanostring and WGS revealed amplification in genes related mainly with kinases proteins, involved in proliferation, migration and cell differentiation, such as EGFR, PIK3CA, and MAPK7. Overexpression of EGFR was confirmed by phospho RTK-array and western blot analysis. Furthermore, the cell was sensitivity to cisplatin, with IC50 13 times lower than SiHa cell line. In conclusion, this cervical cancer cell line presents molecular alterations that are an important tool for leading pre-clinical studies of new drugs that target one or more of the altered pathways.

  Study EGAS00001003343

• Summary statistics from genome-wide association study in glioma of 12,488 cases and 18,169 controls.

  Genome-wide association studies (GWAS) have transformed our understanding of glioma susceptibility, but individual studies have had limited power to identify risk loci. We performed a meta-analysis of existing GWAS and two new GWAS, which totaled 12,488 cases and 18,169 controls. We identified five new loci for glioblastoma (GBM) at 1p31.3 (rs12752552; P = 2.04 × 10-9, odds ratio (OR) = 1.22), 11q14.1 (rs11233250; P = 9.95 × 10-10, OR = 1.24), 16p13.3 (rs2562152; P = 1.93 × 10-8, OR = 1.21), 16q12.1 (rs10852606; P = 1.29 × 10-11, OR = 1.18) and 22q13.1 (rs2235573; P = 1.76 × 10-10, OR = 1.15), as well as eight loci for non-GBM tumors at 1q32.1 (rs4252707; P = 3.34 × 10-9, OR = 1.19), 1q44 (rs12076373; P = 2.63 × 10-10, OR = 1.23), 2q33.3 (rs7572263; P = 2.18 × 10-10, OR = 1.20), 3p14.1 (rs11706832; P = 7.66 × 10-9, OR = 1.15), 10q24.33 (rs11598018; P = 3.39 × 10-8, OR = 1.14), 11q21 (rs7107785; P = 3.87 × 10-10, OR = 1.16), 14q12 (rs10131032; P = 5.07 × 10-11, OR = 1.33) and 16p13.3 (rs3751667; P = 2.61 × 10-9, OR = 1.18). These data substantiate that genetic susceptibility to GBM and non-GBM tumors are highly distinct, which likely reflects different etiology.

  Study EGAS00001003372

• Intratumoural Heterogeneity Underlies Distinct Therapy Responses and Treatment Resistance in Glioblastoma

  Glioblastomas are the most common and lethal neoplasms of the central nervous system. Neighbouring glioma cells maintain extreme degrees of genetic and phenotypic variation that forms intratumoural heterogeneity. This genetic diversity allows the most adaptive tumour clones to show treatment resistance and conserve the tumour bulk. We aimed to model this phenomenon and test the effectiveness of several targeted therapeutic interventions to overcome therapy resistance. Heterogenic tumour masses were first deconstructed into single tumour cells, which were expanded independently as single-cell clones. SNP arrays, whole-genome and RNA sequencing, and CpG methylation analysis validated the unique molecular profile of each tumour clone, which displayed distinct pathologic features, including cell morphology, growth rate, and resistance to temozolomide and ionizing radiation. We also identified variable sensitivities to AURK, CDK, and EGFR inhibitors, consistent with the heterogenic molecular alterations that each clone harboured. These targeted therapies effectively eliminated the temozolomide- and/or irradiation-resistant clones and also parental polyclonal cells. Our findings indicate that polyclonal tumours create a dynamic environment that consists of diverse tumour elements and treatment responses. Designing targeted therapies based on a range of molecular profiles can be a more effective strategy to eradicate treatment resistance, recurrence, and metastasis.

  Study EGAS00001003438

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Colorectal_LCM

  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. 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 EGAS00001003455

• Reconstituting the transcriptome and DNA methylome landscapes of human implantation at single-cell resolution

  Implantation is a milestone event during mammalian embryogenesis. Due to the extreme difficulty of obtaining in vivo human early post-implantation embryos, the gene regulatory network and epigenetics controlling human embryo implantation remains elusive. Here, combining an in vitro culture system for human post-implantation development and single-cell omics sequencing technologies, over 10,000 single cells at five representative stages of pre/post-implantation development were systematically analyzed. Unsupervised dimensionality reduction and clustering algorithm of the transcriptome data show stepwise implantation routes for the epiblast, primitive endoderm, and trophectoderm lineages, suggesting preparation for the establishment of a mother-to-offspring connection after implantation. Female embryos showed asynchronous progress of dosage loss of X chromosomes during implantation. Furthermore, using the single cell trio-Seq (scTrio-Seq) strategy, re-methylation of the genomes of all the three lineages was unambiguously revealed. Surprisingly, the genome re-methylation of PE lineage were much slower than both EPI and TE lineages during the implantation process, indicating distinct methylome features between EPI and PE although both of which were derived from ICM. Collectively, our work paves the way for understanding the complex molecular mechanisms that regulate human embryo implantation, informing new insights and future efforts in early embryonic development and reproductive medicine.

  Study EGAS00001003443

• METABRIC: Data from Batra et al (2021); DNA methylation landscapes of 1538 breast cancers reveal a replication-linked clock, epigenomic instability and cis-regulation

  DNA methylation is aberrant in cancer, but the dynamics, regulatory role and clinical implications of such epigenetic changes are still poorly understood. Here, RRBS-based profiles of 1538 breast tumors and 244 normal breast tissues from the METABRIC cohort are reported, facilitating detailed analysis of DNA methylation within a rich context of genomic, transcriptional, and clinical data. Tumor methylation from immune and stromal signatures are deconvoluted leading to the discovery of a tumor replication-linked clock with genome-wide methylation loss in non-CpG island sites. Unexpectedly, methylation in most tumor CpG islands follows two replication-independent processes of gain (MG) or loss (ML) that we term epigenomic instability. Epigenomic instability is correlated with tumor grade and stage, TP53 mutations and poorer prognosis. After controlling for these global trans-acting trends, as well as for X-linked dosage compensation effects, cis-specific methylation and expression correlations are uncovered at hundreds of promoters and over a thousand distal elements. Some of these targeted known tumor suppressors and oncogenes. In conclusion, the study demonstrates that global epigenetic instability can erode cancer methylomes and expose them to localized methylation aberrations in-cis resulting in transcriptional changes seen in tumors.

  Study EGAS00001004327

• Molecular determinants of response to PD-L1 blockade across tumor types

  Immune checkpoint inhibitors targeting the PD-1/PD-L1 axis lead to durable clinical responses in subsets of cancer patients across multiple indications, including non-small cell lung cancer (NSCLC), urothelial carcinoma (UC) and renal cell carcinoma (RCC). Herein, we complement PD-L1 immunohistochemistry (IHC) and tumor mutation burden (TMB) with RNA-seq in 366 patients to identify unifying and indication-specific molecular profiles that can predict response to checkpoint blockade across these tumor types. Multiple machine learning approaches failed to identify a baseline transcriptional signature highly predictive of response across these indications. Signatures described previously for immune checkpoint inhibitors also failed to validate. At the pathway level, significant heterogeneity was observed between indications, in particular within the PD-L1+ tumors. mUC and NSCLC were molecularly aligned, with cell cycle and DNA damage repair genes associated with response in PD-L1- tumors. At the gene level, the CDK4/6 inhibitor CDKN2A was identified as a significant transcriptional correlate of response, highlighting the association of non-immune pathways to the outcome of checkpoint blockade. This cross-indication analysis revealed molecular heterogeneity between mUC, NSCLC and RCC tumors, suggesting that indication-specific molecular approaches should be prioritized to formulate treatment strategies.

  Study EGAS00001004343

• Comprehensive genomic profiles of small cell lung cancer

  We have sequenced the genomes of 110 small cell lung cancers (SCLC), one of the deadliest human cancers. We found bi-allelic inactivation of TP53 and RB1 in nearly all the tumors analyzed, sometimes by complex genomic rearrangements. Two tumors with wild-type RB1 had evidence of chromothripsis leading to overexpression of Cyclin D1, revealing an alternative mechanism of Rb1 deregulation. Thus, loss of the tumor suppressors TP53 and RB1 is obligatory in SCLC. We further discovered somatic genomic rearrangements of TP73 that create an oncogenic version of this gene, TP73Δex2/3. In few cases, SCLC tumors exhibited kinase gene mutations, providing a possible therapeutic opportunity for individual patients. Finally, we observed inactivating mutations in NOTCH family genes in 25% of human SCLC. Accordingly, activation of Notch signaling in a pre-clinical SCLC mouse model dramatically reduced the number of tumors and extended the survival of the mutant mice. Furthermore, neuroendocrine gene expression was abrogated by Notch activity in SCLC cells. This first comprehensive study of somatic genome alterations in SCLC uncovers several key biological processes and identifies candidate therapeutic targets in this highly lethal form of cancer.

  Study EGAS00001000925

• Exploring the role of mtDNA variation in Multiple Sclerosis in a large cohort of discordant monozygotic twins

  Several lines of evidence indicate that mitochondrial DNA (mtDNA) variants might predispose to multiple sclerosis (MS). We examined this in 49 monozygotic (MZ) twin pairs clinically discordant for MS at study entry. Since the nuclear DNA of MZ twins is identical, our study provides a unique setting to functionally associate unique mtDNA variants and skewed heteroplasmy with MS development. In order to identify mtDNA variants in the twin cohort, we deeply sequenced the mitochondrial genome in blood of all 49 MZ twin pairs using next generation sequencing (Illumina, HiSeq or MiSeq) with an average coverage depth of ~25,000 sequences per base. Heteroplasmic variants were validated using targeted deep sequencing (TDS) (Illumina, MiSeq) in DNA isolated from blood and buccal swabs. In addition, whole blood was also available of 5 additional MS-affected or healthy siblings and 1 dizygotic twin pair discordant for MS, comprising in total 6 non-identical sib pairs. Moreover, all cases were screened for pathogenic mtDNA polymerase gamma (POLG) mutations. The data provides valuable insights in the possible involvement of mtDNA variants in the pathogenesis of MS, and gives information on the segregation of heteroplasmic variants within MZ twins and non-identical siblings and across different tissues.

  Study EGAS00001001240

• Clonal selection and double hit events involving tumor suppressor genes underlie relapse from total therapy

  In order to gain insights into the mechanisms of disease progression in multiple myeloma (MM) patients treated with multi-agent and high-dose chemotherapy we have performed a longitudinal study of 33 patients entered into Total Therapy protocols using gene expression profiling, high resolution copy number analyses and whole exome sequencing. The results of this study show the importance of acquired mutations in MM driver genes, bi-allelic events affecting tumor suppressor genes and increased proliferation rates in driving relapse, which is enhanced by intra-clonal heterogeneity and Darwinian type clonal evolution. Branching was the most frequent evolution pattern, characterized by clonal loss and appearance of new clones. We found a higher number of lost and acquired copy number aberrations and an excess of bi-allelic inactivation of tumor suppressor genes in GEP70 high risk (HiR) cases, consistent with genomic instability being a key feature of HiR MM. In conclusion, our study further stresses the impact of known prognostic markers and driver genes on the fitness level of MM cells to survive multi-agent chemotherapy. The frequent bi-allelic loss of tumor suppressor genes in HiR cases and the potentially resulting failure of treatment with DNA damaging agents highlight the need for new therapies for HiR MM.

  Study EGAS00001001810

• A somatic reference standard for cancer genome sequencing with COLO829

  Identification of somatic alterations in cancer has become feasible with the massive adoption of next generation sequencing. However, due to variability in sequencing and informatics pipelines, a common somatic reference is needed. We thus performed paired whole genome sequencing of a metastatic melanoma cell line (COLO829) and a matched lymphoblastoid line (COLO829BL) across three institutions (TGen, Illumina, Genome Sciences Centre at the British Columbia Cancer Agency). We performed a meta-analysis of all data, in combination with the originally reported analysis of these cell lines by Pleasance et al (PMID: 20016485), and report a somatic reference standard based on consensus events. DNA extractions, library preparation, sequencing, and analysis were separately performed at each site. Results were compiled with sequence data previously generated for the cell lines to identify true positive variants, collectively representing a somatic reference standard. Overall, common somatic events include point mutations, small insertion/deletions, and genes affected by concurrent copy number changes. We additionally show divergence of somatic mutations in COLO829 associated with differences in cell line lineage. We present this reference to the community as a standard for enabling interpretation and identification of somatic mutations across institutions and analytical pipelines.

  Study EGAS00001001385

• The HLA-DQβ1 insertion is a strong achalasia risk factor and displays a geospatial north-south gradient among Europeans

  Idiopathic achalasia is a severe motility disorder of the esophagus and is characterized by a failure of the lower esophageal sphincter to relax due to a loss of neurons in the myenteric plexus. Most recently, we identified an 8 amino acid insertion in the cytoplasmic tail of HLA-DQβ1 as strong achalasia risk factor in a sample set from Central Europe, Italy, and Spain. Here, we tested whether the HLA-DQβ1 insertion also confers achalasia risk in the Polish and Swedish population. We could replicate the initial findings and the insertion shows strong achalasia association in both samples (Poland P = 1.84×10-04, Sweden P = 7.44×10-05). Combining all five European data sets – Central Europe, Italy, Spain, Poland, and Sweden – the insertion is achalasia associated with Pcombined = 1.67×10-35. In addition, we observe that the frequency of the insertion shows a geospatial north-south gradient. The insertion is less common in northern (around 6-7% in patients and 2% in controls from Sweden and Poland) compared to southern Europeans (around 16% in patients and 8% in controls from Italy) and shows a stronger attributable risk (AR) in the southern European population. Our study provides evidence that the prevalence of achalasia may differ between populations.

  Study EGAS00001001515

• Whole exome sequencing on primary retinoblastoma tissues and matching lymphocyte DNA.

  Retinoblastoma is a rare childhood cancer initiated by RB1 mutation or MYCN amplification, while additional alterations may be required for tumor development. However, the view on single nucleotide variants is very limited. To better understand oncogenesis, we determined the genomic landscape of retinoblastoma. We performed exome sequencing of 71 retinoblastomas and matched blood DNA. Next, we determined the presence of single nucleotide variants, copy number alterations and viruses. Aside from RB1, recurrent gene mutations were very rare. Only a limited fraction of tumors showed BCOR (7/71, 10%) or CREBBP alterations (3/71, 4%). No evidence was found for the presence of viruses. Instead, specific somatic copy number alterations were more common, particularly in patients diagnosed at later age. Recurrent alterations of chromosomal arms often involved less than one copy, also in highly pure tumor samples, suggesting within-tumor heterogeneity. Our results show that retinoblastoma is among the least mutated cancers and signify the extreme sensitivity of the childhood retina for RB1 loss. We hypothesize that retinoblastomas arising later in retinal development benefit more from subclonal secondary alterations and therefore, these alterations are more selected for in these tumors. Targeted therapy based on these subclonal events might be insufficient for complete tumor control.

  Study EGAS00001001690

• Ethiopia_Genome_Project__low_coverage_

  The Ethiopian area stands among the most ancient ones ever occupied by human populations and their ancestors. Particularly, according to archaeological evidences, it is possible to trace back the presence of Hominids up to at least 3 million years ago. Furthermore, the present day human populations show a great cultural, linguistic and historic diversity which makes them essential candidate to investigate a considerable part of the African variability. Following the typing of 300 Ethiopian samples on Illumina Omni 1M (see Human Variability in Ethiopia project, previously approved by the Genotyping committee) we now have a clearer idea on which populations living in the area include the most of the diversity.This project therefore aims to sequence the whole genome of 300 individuals at low (4-8x) depth belonging to the six most representative populations of the Ethiopian area to produce a unique catalogue of variants peculiar of the North East Africa. Furthermore 6 samples (one from each population) will also be sequenced at high (30x) depth to ensure full coverage of the diversity spectrum.The retrieved variants will be of great help in evaluating the demographic dynamics of those populations as well as shedding light on the migrations out of Africa.

  Study EGAS00001000238

• UK10K OBESITY TWINSUK

  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. These TwinsUK samples are a subset of the TwinsUK cohort samples with a BMI>=40 that will be exome sequenced by the Obesity group.

  Study EGAS00001000306