22929 results for "ega"

in 48.17 milliseconds.

• Evolution of an adenomacarcinoma in response to selection by targeted kinase inhibitors

  Background - Adenocarcinomas of the tongue are rare and represent the minority (20-25%) of salivary gland tumors affecting the tongue. We investigated the utility of massively parallel sequencing to characterize an adenocarcinoma of the tongue, before and after treatment. Results - In the pre-treatment tumor we identified 7,629 genes within regions of copy number gain. 1,078 genes exhibited increased expression relative to the blood and unrelated tumors and four genes contained somatic protein-coding mutations. Our analysis suggested the tumor cells were driven by the RET oncogene. Genes whose protein products are targeted by the RET inhibitors sunitinib and sorafenib correlated with being amplified and or highly expressed. Consistent with our observations administration of sunitinib was associated with stable disease lasting 4 months, after which the lung lesions began to grow. Administration of sorafenib and sulindac provided disease stabilization for an additional 3 months after which the cancer progressed and new lesions appeared. A recurring metastasis possessed 7,288 genes within copy number amplicons, 385 genes exhibiting increased expression relative to other tumours and 9 new somatic protein coding mutations. The observed mutations and amplif ications were consistent with therapeutic resistance arising through activation of the MAPK and AKT pathways. Conclusions - We conclude that complete genomic characterization of a rare tumour has the potential to aid in clinical decision making and identifying therapeutic approaches where no established treatment protocols exist. These results also provide direct in-vivo genomic evidence for mutational evolution within a tumor under drug selection and potential mechanisms of drug resistance accrual.

  Study EGAS00000000074

• Further details about data access are available at KNIH DAC.

  Dac EGAC00001000487

• Pharmacological improvement of CFTR function rescues airway epithelial homeostasis and host defense in cystic fibrosis children

  Rationale: Recent clinical trials have shown that elexacaftor/tezacaftor/ivacaftor (ETI) provides significant benefit to patients with cystic fibrosis (CF) by improving CFTR function. Despite these promising clinical findings, the specific effects on the airway inflammatory status leading to structural damage and impaired lung function, remain elusive. Objectives: We investigated the transcriptional changes at the single cell level of airway epithelial and immune cells that underlie the clinical improvement by ETI therapy in children with CF. Methods: Nasal swabs from thirteen children with CF and at least one F508del allele aged 6 to 12 years were collected at baseline and three months after initiation of ETI and subjected to scRNA-seq. In addition, we collected nasal swabs from 12 age- and sex-matched controls. Sweat chloride concentrations, spirometry and multiple breath washout were used for clinical evaluation of CF children. Measurements and main results: ETI significantly reduced sweat chloride concentrations (-52.2 ± 14.3 mmol/L, P < 0.001) and improved the lung clearance index (-1.3 ± 1.8, P < 0.05). Single cell transcriptomics revealed an impaired interferon signalling in epithelial cells of CF children at baseline, which was partially restored by ETI in conjunction with an ETI-induced increase in expression of MHC I and II encoding genes similar to control levels. Additionally, ETI markedly reduced the inflammatory phenotype of immune cells, particularly of neutrophils and macrophages. Conclusions: Improvement of CFTR function by ETI restores epithelial homeostasis and reduces immune cell inflammatory responses in the upper airways of children with CF, highlighting the potential of early initiation of ETI therapy.

  Study EGAS50000000128

• Systemic Lupus Erythematosus Serum Stimulation of Human Intestinal Organoids Induces Barrier Leakiness and Changes in Goblet Cell Differentiation

  Human intestinal epithelial cells are the interface between potentially harmful luminal content and basally residing immune cells. Their role is not only nutrient absorption but also the formation of a tight monolayer that constantly secretes mucus creating a multi-layered protective barrier. Alterations in this barrier can lead to increased gut permeability which is frequently seen in individuals with chronic extraintestinal autoimmune diseases, such as Systemic Lupus Erythematosus (SLE). Despite recent advances in identifying alterations in gut microbiota composition in SLE patients, not much attention has been given to the epithelial barrier itself. To date, it remains largely unexplored which role and function intestinal epithelial cells have in SLE pathology. Here, we present a unique near-physiologic in vitro model specifically designed to examine the effects of SLE on the epithelial cells. We utilize human colon organoids that are stimulated with serum obtained from SLE patients. Combining bulk and single-cell RNA transcriptomic analysis with functional assays revealed that SLE serum stimulation induced a unique expression profile marked by a reduction of goblet cell marker expression and mucus composition. Additionally, organoids exhibited imbalanced cellulat composition along with enhanced paracellular permeability,mitochondrial function, and a type I interferon gene signature. Similarly, transcriptomic analysis of SLE human colon biopsies revealed a downregulation of epithelial secretory markers. Our work uncovers a crucial connection between SLE and intestinal homeostasis that might be promoted in vivo through the blood, offering insights into the causal connection of barrier dysfunction and autoimmune diseases.

  Study EGAS50000000012

• Carboxylesterase 1 mediates a distinctive metabolic profile of dendritic cells to attain an inflammatory phenotype

  The metabolic profile of dendritic cells (DCs) shapes their phenotype and functions. Carboxylestrase 1 (CES1) enzyme is highly expressed in mononuclear myeloid cells however its exact role in DCs is elusive. We used a CES1 inhibitor (WWL113) and genetic overexpression to explore the role of CES1 in DCs differentiation in inflammatory models. CES1 expression was analyzed during CD14+ monocytes differentiation to DCs (MoDCs) using quantitative PCR. CES1 Inhibitor (WWL113) was applied during MoDCs differentiation. Surface markers, secreted cytokines, lactic acid production, phagocytic and T cell polarization capacity were analyzed. Transcriptomic and metabolic profile were assessed with RNA-sequencing and mass spectrometry. Cellular respiration was assessed with seahorse respirometry. Transgenic mice were used to assess CES1 overexpression in DCs in inflammatory models. CES1 expression peaks early during MoDCs differentiation. Pharmacological inhibition of CES1 led to higher expression of CD209, CD86 and MHCII. WWL113 treated MoDCs secreted higher quantities of IL6, IL8, TNF and IL10 and demonstrated stronger phagocytic ability and higher capacity to polarize Th17 differentiation in autologous DCs-T cells co-culture model. Transcriptomic profiling revealed enrichment of multiple inflammatory and metabolic pathways. Functional metabolic analysis shows impaired maximal mitochondrial respiration capacity, increased lactate production and decreased intracellular amino acids and TCA intermediates. Transgenic human CES1 overexpression in murine DCs generated less inflammatory phenotype and increased resistance to T cell mediated colitis. In conclusion, CES1 inhibition directs DCs differentiation towards more inflammatory phenotype, that shows stronger phagocytic capacity and supports Th17 skewing. This is associated with disrupted mitochondrial respiration and amino acids depletion.

  Study EGAS50000000230

• Cerebrospinal fluid cfDNA sequencing for classification of central nervous system glioma

  Primary central nervous system (CNS) gliomas can be classified by characteristic genetic alterations. In addition to solid tissue obtained by surgery or biopsy, cell-free DNA (cfDNA) derived from cerebrospinal fluid (CSF) is an alternative source of material for genomic analyses. Experimental design: We performed targeted next-generation sequencing (NGS) of CSF cfDNA in a representative cohort of 85 patients presenting with with suspicion of primary or recurrent glioma at two neurooncological centers. Copy-number variation (CNV) profiles, single nucleotide variants (SNVs), and small insertions/ deletions (indels) were combined into a molecular-guided tumor classification. Comparison with the solid tumor was performed for 38 cases with matching solid tissue available. Results: Cases were stratified into four groups: glioblastoma (n = 32), other glioma (n = 19), non-malignant (n = 17) and no confirmed diagnosis available (n = 17). We introduce a molecular-guided tumor classification, enabling identification of tumor entities and/ or cancer specific alterations in 75.0 % (n = 24) of glioblastoma and 52.6 % (n = 10) of other glioma cases. The overlap between CSF and matching solid tissue was highest for CNVs (26-48 %) and SNVs at pre-defined gene loci (44 %), followed by SNVs/indels identified via uninformed variant calling (8-14 %). A molecular-guided tumor classification was made 23.5 % (n = 4) of cases with no confirmed diagnosis available. Conclusions: We developed a workflow for targeted sequencing of CSF cfDNA as well as a strategy for interpretation and reporting of sequencing results based on a molecular-guided tumor classification in glioma.

  Study EGAS50000000060

• Circulating tumor DNA dynamics reveal KRAS G12C mutation heterogeneity and response to treatment with the KRAS G12C inhibitor divarasib in solid tumors

  Purpose: To inform prognosis, treatment response, disease biology, and KRAS G12C mutation heterogeneity, we conducted exploratory circulating tumor DNA (ctDNA) profiling on 134 patients with solid tumors harboring a KRAS G12C mutation treated with single-agent divarasib (GDC-6036) in a phase 1 study. Experimental design: Plasma samples were collected for serial ctDNA profiling at baseline (Cycle 1 Day 1 prior to treatment) and multiple on-treatment time points (Cycle 1 Day 15 and Cycle 3 Day 1). Results: KRAS G12C ctDNA was detectable from plasma samples in 72.9% (43/59) and 92.6% (50/54) of patients with non-small cell lung cancer and colorectal cancer, respectively, the majority of whom were eligible for study participation based on a local test detecting the KRAS G12C mutation in tumor tissue. Baseline ctDNA tumor fraction was associated with tumor type, disease burden, and metastatic sites. A decline in ctDNA level was observed as early as Cycle 1 Day 15. Serial assessment showed a decline in ctDNA tumor fraction associated with response and progression-free survival. Except for a few cases of KRAS G12C sub-clonality, on-treatment changes in KRAS G12C variant allele frequency mirrored changes in the overall ctDNA tumor fraction. Conclusion: Across tumor types, the KRAS G12C mutation likely represents a truncal mutation in the majority of patients. Rapid and deep decline in ctDNA tumor fraction was observed in patients responding to divarasib treatment. Early on-treatment dynamics of ctDNA were associated with patient outcomes and tumor response to divarasib treatment.

  Study EGAS50000000315

• Lactate metabolism in cancer stem cell fate regulation

  Tumors arise as a result of uncontrolled cell proliferation driven by mutations in the signaling pathways that regulate the balance between cell division and differentiation. Despite of that, in the tumor, cancer stem cells seem to retain their capability to differentiate, which can lead to differential chemotherapy response. We previously showed that intestinal stem cells and differentiated cells display specific metabolic features and that they interact through metabolites. Moreover, we found that altering mitochondrial metabolism is sufficient to lead and define stem cell differentiation. In cancer, cellular metabolism is highly derailed. However, how metabolic changes in the cancer cell contribute to cell fate decisions remains unknown. Recently, the crosstalk between metabolism and epigenetics has emerged as one of the players in tumor development. Here we aim to address the impact of the altered metabolism of cancer cells on their epigenome as a regulatory axis of cell identity and fate. We use human-derived tumor organoids where we introduce genetic encoded fluorescent-based reporters to visualize cell identity and metabolism at single cell level. We combine 4D live imaging with a recently developed machine learning method to re-construct cell-lineages, trace stem cell state and metabolic changes during tumor organoid development. We find that, mimicking the tumor metabolic microenvironment (Warburg effect-driven low-glucose high-lactate) leads to changes cellular metabolism and histone acetylation. Strikingly, these alterations provoke profound changes in tumor development by increasing the population of cancer stem cells and their proliferation and differentiation dynamics.

  Study EGAS50000000063

• Tumor Evolution Analysis Uncovered Immune-Escape Related Mutations in Relapsed Diffuse Large B-Cell Lymphoma

  More than 60% of patients with primary diffuse large B-cell lymphoma (DLBCL) achieve durable remissions upon the chemoimmunotherapy with R-CHOP. However, the outcome of the relapse patients is extremely poor. To characterize the molecular features that underlie treatment resistance, we performed whole-exome sequencing (n=44) and deep targeted sequencing (n=27) of 44 diagnosis-relapse DLBCL pairs. We found that the overall landscape of genetic alterations was concordant between diagnosis and relapse DLBCL. By a deep comparison of the longitudinal sequencings, we identified relapse-enriched mutations that involved in evasion of immune surveillance (CD58, TNFRSF14, PKD1, FAS), resistance to DNA damaging agents (ATM, TP53), and critical tumor suppressor gene (ZFHX3). Knockdown of ATM and ZFHX3 resulted in increased drug resistance in both germinal center B cell like (GCB) and activated B cell like (ABC) DLBCL cell lines in vitro, while the knockdown of TNFRSF14 and FAS conferred resistance specifically in GCB cell lines. Furthermore, we summarized two evolutionary patterns of DLBCL relapse as the “durable retention” model with stable clone structure, and the “selective sweep” model which exhibited significant clonal switch. Of note, genetic alterations leading to attenuated immune responses usually presented as trunk mutations in the “durable retention” model (7/12) while were newly acquired or expanded upon relapse in the “selective sweep” model (9/16). Overall, our study demonstrated that relapse DLBCL derives from distinct clonal evolutionary patterns, and highlighted immune escape as an important contributor of tumor relapse.

  Study EGAS50000000032

• Surgery in combination with immune checkpoint therapy as an effective treatment for patients with metastatic cancer.

  Abstract Surgical removal of primary tumors was shown to reverse tumor-mediated immune suppression in pre-clinical models with metastatic disease. However, how cytoreductive surgery in the metastatic setting modulates the immune responses in patients, especially in the context of immune checkpoint therapy (ICT)-containing treatments is not understood. Here, we report the first prospective, non-comparative clinical trial (N=104) using three different ICT-containing strategies plus cytoreductive or “debulking” surgery to remove the primary tumor-bearing kidney or a metastatic lesion as a treatment for patients (N=43) with metastatic clear cell renal cell carcinoma (mccRCC). For those patients (N=61) who were not candidates for cytoreductive surgery, a biopsy was obtained instead for correlative biological studies. Our data demonstrated that the combination of ICT with cytoreductive surgery was safe and feasible in patients with mccRCC. The 2-year overall survival was 84% with a median OS of 54.7 months for patients who received ICT containing regimens plus surgery. Immune-monitoring studies with co-detection by indexing (CODEX) identified distinct tumor spatial conformation of cellular subsets as a novel and improved predictor of response to ICT. Importantly, single-cell RNA-sequencing data demonstrated that surgical removal of the tumor increased antigen-presenting dendritic cell population with a concurrent reduction in KDM6B expressing immune-suppressive myeloid cells in the peripheral blood. Together, this study highlighted the feasibility of combining ICT with cytoreductive surgery in a metastatic setting and demonstrated potential enhancement of immune responses following ICT plus cytoreductive surgery in patients with metastatic disease.

  Study EGAS00001005667

• Impact of Glycolysis Inhibition on the Epigenome of Synovial Fluid T Cells in Juvenile Idiopathic Arthritis

  Juvenile Idiopathic Arthritis (JIA) encompasses a group of childhood-onset autoimmune conditions characterized by joint inflammation and mononuclear cell infiltration, especially of activated CD4+ memory/effector T (Tmem/Teff) cells in synovial fluid (SF). JIA CD4+ T cells exhibit a distinct inflammation-associated epigenomic profile, though the underlying molecular mechanisms remain unclear. This study investigates the role of glycolysis in shaping the epigenomic landscape of JIA SF CD4+ T cells, with a focus on H3K27ac histone modifications. Using an ex vivo model of antigen stimulation, CD4+ T cells were isolated from peripheral blood (PB) of healthy controls (HC) and synovial fluid from JIA patients. Cells were activated in vitro with anti-CD3/CD28 for 24 hours, and changes in H3K27ac were assessed by chromatin immunoprecipitation followed by sequencing (ChIP-seq). To explore the specific role of glycolysis, cells were also treated with 2-deoxyglucose (2DG), an inhibitor of glucose-6-phosphate production. This inhibition impacted both H3K27ac acetylation patterns and gene expression associated with the activated epigenomic regions. Additionally, HC CD4+ T cells were activated in vitro in the presence of 30% SF or 30% plasma from JIA patients. RNA sequencing (RNA-seq) was performed to examine SF-induced transcriptional changes, providing insight into the metabolic dysregulation and epigenomic alterations observed in JIA SF CD4+ T cells. This dataset includes ChIP-seq and RNA-seq data that collectively advance our understanding of how glycolysis and the SF environment contribute to the epigenetic and transcriptional regulation in JIA CD4+ T cells.

  Study EGAS50000000808

• Circulating Tumor DNA Analysis in ERBB2-Amplified Colorectal Cancer: Biomarker Analysis of the MyPathway Trial

  Combination of two HER2-directed antibodies, pertuzumab and trastuzumab (P+T), has antitumor activity in HER2-positive colorectal cancer (CRC). Although liquid biopsies are increasingly being used in clinical oncology, the association between tumor and circulating tumor DNA (ctDNA) ERBB2 status and ctDNA monitoring for early response and resistance are unknown. Eighty-five patients with ERBB2-amplified and/or overexpressed CRC were treated with P+T in the MyPathway trial; 42 had ctDNA testing at Cycle 1 Day 1 (C1D1), and 38 had longitudinal plasma tested for ctDNA. We analyzed the ctDNA vs tissue ERBB2 concordance, genomic co-alterations, and ctDNA dynamics and association with response. Forty-one (98%) of 42 patients had genomic alterations detected in ctDNA at C1D1, and 29 (69%) had ERBB2 amplification in ctDNA. There was a strong correlation between ERBB2 copy number on next-generation sequencing in tissue and C1D1 ERBB2 ctDNA copy number. Thirty-seven percent achieved a molecular response by C3D1 on P+T, which was associated with prolonged progression-free survival and overall survival (OS). CDKN2A and KRAS mutations were associated with shorter OS, and a trend was seen with PIK3CA mutations. Several emerging co-alterations were identified in ctDNA at progression, including in the MAPK and PI3K pathways and other tyrosine receptor kinases. ctDNA can detect ERBB2 amplification in many, but not all, patients with ERBB2 amplification detected in tumor samples. ctDNA molecular response was associated with better survival, and ctDNA co-alterations may offer insights into mechanisms of intrinsic and acquired resistance.

  Study EGAS50000000916

• Nanopore sequencing of blood, saliva, buccal mucosa samples of patients with inherited retinal disorders

  X-linked genetic disorders typically affect females less severely than males due to the presence of a second X chromosome not carrying the deleterious variant. However, the phenotypic expression in females is highly variable, which may be explained by an allelic skew in X chromosome inactivation. Accurate measurement of X inactivation skew is crucial to understand and predict disease phenotype in carrier females, with prediction especially relevant for degenerative conditions. We propose a novel approach using nanopore sequencing to quantify skewed X inactivation accurately. By phasing sequence variants and methylation patterns, this single assay reveals the disease variant, X inactivation skew, its directionality, and is applicable to all patients and X-linked variants. Enrichment of X-chromosome reads through adaptive sampling enhances cost-efficiency. Our study includes a cohort of 16 X-linked variant carrier females affected by two X-linked inherited retinal diseases: choroideremia and RPGR-associated retinitis pigmentosa. As retinal DNA cannot be readily obtained, we instead determine the skew from peripheral samples (blood, saliva and buccal mucosa), and correlate it to phenotypic outcomes. This revealed a strong correlation between X inactivation skew and disease presentation, confirming the value in performing this assay and its potential as a way to prioritise patients for early intervention, such as gene therapy currently in clinical trials for these conditions. Our method of assessing skewed X inactivation is applicable to all long-read genomic datasets, providing insights into disease risk and severity and aiding in the development of individualised strategies for X-linked variant carrier females.

  Study EGAS50000000440

• PCR-free shallow whole genome sequencing for chromosomal copy number detection from plasma of cancer patients is an efficient alternative to the conventional PCR-based approach

  Somatic copy number alterations (SCNAs) can be detected in tumor cell-free DNA (cfDNA) by shallow whole genome sequencing (sWGS). Polymerase chain reaction (PCR) is typically included in library preparations, but a PCR-free method could be an effective alternative for research and diagnostics. Peripheral blood samples (n=22) were collected in EDTA-containing tubes from non-small cell lung cancer patients (n=10) and healthy donors (n=12). PCR and PCR-free library preparations underwent sWGS and the following metrics were compared: mapping quality, proportion of duplicate reads, genome coverage, sequencing read length and copy number profiles. The PCR-free method was further evaluated on archived plasma samples (n=33) of varying pre-analytical quality collected in EDTA or lithium heparin tubes. Archived samples were from acute myeloid leukaemia (n=15) and B-cell lymphoma (n=18) patients. The percentage of unique reads was significantly higher for the 22 samples analysed by PCR-free (95.7%) compared to PCR (84%). PCR-free libraries from 33 archived samples produced data of a comparable quality, with 95.4% unique reads. All other evaluated metrics were highly comparable for PCR and PCR-free library preparations. We conclude that SCNA detection can be performed by PCR-free sWGS on cfDNA from cancer patients and is feasible for diagnostic and archival plasma collected in EDTA- or lithium heparin-containing tubes. Our results pave the way for an automated cfDNA workflow for cancer plasma samples.

  Study EGAS00001004692

• Coding and regulatory somatic profiling of triple-negative breast cancer in Sub-Saharan African patients

  The burden of triple-negative breast cancer (TNBC) may be shaped by genetic factors, particularly inherited and somatic mutation profiles. However, data on this topic remain limited, especially for the African continent, where a higher TNBC incidence is observed. In the age of precision medicine, cataloguing TNBC diversity in African patients becomes imperative. We performed whole exome sequencing, including untranslated regions, on 30 samples from Angola and Cape Verde, which allowed to ascertain on potential regulatory mutations in TNBC for the first time. A high somatic burden was observed for the African cohort, with 86% of variants being so far unreported. Recurring to predictive functional algorithms, 17% of the somatic single nucleotide variants were predicted to be deleterious at the protein level, and 20% overlapped with candidate cis-regulatory elements controlling gene expression. Several of these somatic functionally-impactful mutations and copy number variation (mainly in 1q, 8q, 6 and 10p) occur in known BC- and all cancer-driver genes, enriched for several cancer mechanisms, including response to radiation and related DNA repair mechanisms. TP53 is the top of these known BC-driver genes, but our results identified possible novel TNBC driver genes that may play a main role in the African context, as TTN, CEACAM7, DEFB132, COPZ2 and GAS1. These findings emphasize the need to expand cancer omics screenings across the African continent, the region of the globe with highest genomic diversity, accelerating the discovery of new somatic mutations and cancer-related pathways.

  Study EGAS50000001013

• Whole-genome and transcriptome versus panel sequencing in precision oncology: A translational-clinical comparison

  Precision oncology offers new cancer treatment options. While the use of whole-genome (WGS) and transcriptome (TS) sequencing remains often confined to specific clinical programs, molecular analysis with gene panels is more common clinical practice. Methods: Twenty patients with rare or advanced tumors sequenced by WGS/TS within the DKFZ/NCT/DKTK-MASTER program from 2015 to 2020 (MASTER 1) were additionally sequenced with the TSO500 DNA and TST170 RNA Assay by Illumina. Original Molecular Tumor Board (MTB) therapy recommendations (TR) were juxtaposed with TR from NGS panel data analysis (comparison 1). To mitigate temporal differences between the two sequencings, WGS/TS was bioinformatically reanalyzed (MASTER 2), and new TR were again compared (comparison 2). Results: In comparison 1, 47.1% of the MASTER 1’s TR and 62.8% of the panel’s TR were identical, as were 45.9% of the TR of MASTER 2 and 54.9% of the panel in comparison 2. In both MASTER 1 and MASTER 2, 26.5% and 36.1% of TR, respectively, relied on Biomarkers (BM) which were not covered by the panel. Eight of the ten MTB-derived therapy implementations were also supported in the panel, the two remaining ones being based on BM absent from the panel. Conclusions: About half of the TR in WGS/TS and panel sequencing were identical. Around a third of TR issued by WGS/TS were based on BM which are not covered by the panel, leading to therapy implementations. The presented data underscores the clinical value of comprehensive molecular analyses in individuals with advanced cancers.

  Study EGAS50000000431

• Mutational signatures in normal tissues of patients treated with clinical stage G-quadruplex binder CX5461

  Drug induced photoreactions following exposure to UVA have been documented for several therapeutic compound classes, including psoralens and fluoroquinolones. CX5461, a quinolone-derived small molecule currently in phase 1 studies for cancer treatment, has been shown to be a clinically significant UV photosensitizer. Here we examined the mutagenic spectrum of CX5461 in normal tissues of patients undergoing treatment, by comparing deep whole genome sequencing (WGS) from skin biopsies taken from non-sun exposed areas pre and post drug exposure of patients in the Canadian Cancer Clinical Trials Group (CCTG) IND.231 trial (NCT02719977). We show that a low burden of SNVs (<100) and indels (<40) mutations was detected in post CX5461 skin biopsies, with a pattern resembling background or sequencing artefact signatures. To further characterize the mutational pattern of CX5461 photoreactivity, we compared mutational profiles of human retinal pigment epithelial cells exposed to UVA, CX5461, or a combination of both UV and CX5461. Treatment with CX5461 or UVA alone, resulted in a low SNV burden. In contrast, human cells exposed to UVA and CX5461 had a two to four orders higher mutational burden, with a distinct mutational pattern characterized by T>A, T>C mutations, dissimilar to that of UV alone or CX5461 alone. Our results suggest that human exposure to CX5461 alone is not associated with high levels of single base mutagenesis. However, CX5461 is a very potent UV photosensitizer, underscoring the need for adequate UV screening and controls when studying the effects of photosensitizing small molecules on single base DNA variation.

  Study EGAS50000001145

• Transcriptome profiling of megakaryocytes and platelets: application to GP9- and IKZF5-related thrombocytopenia

  Platelets are anucleate cells produced in the bone marrow and derived from large progenitor cells called megakaryocytes (MKs). Platelets receive RNA transcripts from their progenitorial MKs during thrombopoiesis. However, the correspondence between platelet and MK transcriptomes is poorly understood, particularly in the context of germline mutations that cause platelet formation defects or thrombocytopenia. We have studied the effects of two such mutations on MK and platelet transcriptomes. We generated immortalized MK cell line (imMKCL)-based models of Bernard-Soulier syndrome and IKZF5-related thrombocytopenia. MKs derived from imMKCLs with either a homozygous deletion of GP9 (GP9-/-) or a heterozygous Y121F variant in IKZF5 (IKZF5WT/Y121F) exhibited reduced proplatelet formation (reductions of 96% and 57%, respectively). Platelets from patients with either GP9-/- or IKZF5WT/Y121F genotypes had broad transcriptomic dysregulation, suggesting that (pro)platelet formation defects due to mutations in glycoprotein receptor and transcription factor genes such as GP9 and IKZF5 affect the trafficking of transcripts from MKs into platelets. RNA-seq data from MKs at four stages of differentiation revealed widespread but distinct changes in expression over time between the GP9-/- and the IKZF5WT/Y121F genotypes. Dysregulated genes in GP9-/- MKs were enriched for RNA metabolism and actin/tubulin folding pathways while those in IKZF5WT/Y121F MKs were enriched for cell cycle pathways. Most of these genes were also dysregulated in the platelets of patients with the corresponding diseases. Our results suggest that patients with inherited forms of thrombocytopenia present with a transcript trafficking defect during platelet formation.

  Study EGAS50000001276

• Complex genotype-phenotype relationships shape the response to treatment of Down Syndrome Childhood Acute Lymphoblastic Leukaemia

  Extensive genetic and epigenetic variegation has been demonstrated in many malignancies. Importantly, their interplay has the potential to contribute to disease progression and treatment resistance. To shed light on the complex relationships between these different sources of intra-tumour heterogeneity, we explored their relative contributions to the evolutionary dynamics of Acute Lymphoblastic Leukaemia (ALL) in children with Down syndrome, which has particularly poor prognosis. We quantified the tumour propagating potential of genetically distinct sub-clones using serial transplantation assays and SNP-arrays. While most leukaemias were characterized by a single dominant subclone, others were highly heterogeneous. Importantly, we provide clear and direct evidence that genotypes and phenotypes with functional relevance to leukemic progression and treatment resistance can co-segregate within the disease. Hence, individual genetic lesions can be restricted to well-defined cell immunophenotypes, corresponding to different stages of the leukemic differentiation hierarchy and varied proliferation potentials. As a result of this difference in fitness, which can be accurately quantified via competitive transplantation assays, matching diagnostic, post-treatment, and relapse leukaemias can be dominated by different genotypes, including pre-leukemic clones persisting throughout the disease progression and treatment. Intriguingly, plasticity also appears to be a temporally defined property that can segregate with genotype. These results suggest that Down Syndrome ALL should be viewed as a complex matrix of cells exhibiting genetic and epigenetic heterogeneity that foster extensive clonal evolution and competition. Therapeutic intervention reshapes this ‘eco-system’ and may provide the right conditions for the preferential expansion of selected compartments and subsequently relapse.

  Study EGAS50000001287

• Heterogeneity of IKZF1 genomic alterations and risk of relapse in childhood B-cell precursor acute lymphoblastic leukemia

  Purpose Genomic alterations of IKZF1 are common and associated with adverse clinical features in B-ALL, but the relationship between the type of IKZF1 alteration, disease subtype and outcome are incompletely understood. Patients and Methods We used transcriptome and genomic sequencing, and SNP microarray data to determine leukemia subtype and genomic alterations for 688 pediatric patients with B-ALL enrolled into St. Jude Total Therapy 15 and 16 studies. Results IKZF1 alterations were identified in 115 (16.7%) patients, commonly in BCR::ABL1 (78%) and CRLF2-rearranged, BCR::ABL1-like B-ALL (70%), and were associated with five-year cumulative incidence of relapse (CIR) of 14.8 ± 3.3% compared to 5.0 ± 0.9% for patients without any IKZF1 alteration (P < 0.0001). IKZF1 deletions of exon 4-7 (P = 0.0002), genomic IKZF1plus with any IKZF1 deletion (P = 0.006) or with focal IKZF1 deletion (P = 0.0007), and unfavorable genomic subtypes (P < 0.005) were independently adversely prognostic. Patients with both IKZF1 exon 4-7 deletion and unfavorable genomic subtype had a significantly increased risk of relapse (HR = 58.3; 95% CI, 11.9 – 285.4, P < 0.0001) whereas IKZF1 sequence mutations were not independently predictive of outcome. Conclusions Genomic IKZF1plus with any IKZF1 deletion, IKZF1 deletion of exon 4-7, and unfavorable subtype confer increased risk of relapse. The combination of unfavorable genomic subtype and IKZF1 deletion of exon 4-7 identified patients at greatest risk of relapse despite MRD-directed therapy. The type of IKZF1 alteration together with subtype are informative for risk stratification and predict response in patients with B-ALL.

  Study EGAS50000000106

• WES analysis of a mixed cohort of pituitary tumors

  BACKGROUND: Sporadic pituitary adenomas are the second most common primary brain tumor, but their genetics remains incompletely understood. We performed deep DNA sequencing of 42 pituitary macroadenomas, including hormonally active and atypical adenomas. METHODS: Pituitary adenomas representative of a spectrum of histopathologic subtypes, primary and recurrent tumors, and MIB-1 proliferative indices were profiled. Whole-exome sequencing was performed on tumors and paired blood followed by analysis of somatic copy number alteration and mutation. RESULTS: On the whole, the tumors harbored a low rate of mutations, similar to other benign brain tumors. However, almost one-third (29%) exhibited somatic copy-number alterations across large fractions of the genome (up to 99%). Despite such widespread genomic disruption, these tumors had few focal events, which is atypical among highly-disrupted cancers. The other 71% of tumors formed a distinct molecular class, with somatic copy-number alterations involving less than 6% of the genome. Among the highly-disrupted group, 92% were functional or atypical null-cell adenomas, whereas 87% of the less-disrupted group were clinically non-functional adenomas. The only significantly recurrent events were arm-level losses of chromosomes 1, 2, 11, and 18. No significantly recurrent mutations or insertion/deletion events were identified, suggesting no genes are altered by exomic mutations in &gt;15% of pituitary tumors. CONCLUSIONS: Genomic characterization of a large series of sporadic pituitary adenomas reveals two classes of tumor, based on genomic disruption from arm-level copy-number alterations, which associate with hormonal and histologic subtypes.

  Study EGAS00001001714

• Genome-Wide Association Study of aspirin-induced PUD in a UK cohort

  Aspirin is a widely prescribed non-steroidal anti-inflammatory drug used as an effective analgesic and at low doses as anti-platelet drug for the prevention of cardiovascular diseases. However, aspirin is associated with a significant incidence of upper gastrointestinal (GI) complications, notably peptic ulceration. This study aims to identify genetic risk factors associated with the pathogenesis of aspirin-induced upper GI ulceration.A cohort of 248 patients with aspirin-induced, endoscopiclly confirmed GI ulceration and 491 controls (ulcer, without aspirin or aspirin without ulcer) were genotyped for 1.524,956 SNPs using the Illumina Omni 2.5 SNP array. A replication cohort (206 NSAID-induced ulcer cases [of which 87 were aspirin] and 309 controls) was subsequently genotyped for identified SNP association signals using Sequenom iPLEX Massarry. Logistic regression analysis of identified two SNP signal with suggestive association with aspirin-induced GI ulceration. A SNP within the EYA1 locus (rs12678747) was found to be associated with aspirin GI ulceration in both the discovery cohort (p=3.79x10-7, OR=2.06 ([1.55-2.73]) and, to a lesser extent, in the replication cohort (p=.0005, OR=1.71[1.17-2.50]. Combined analysis of both cohorts yielded a genome wide significant association (p=1.46x10-8). The rs12678747 genotype was found to be associated with significantly higher EYA1 transcript in ulcer biopsy tissue (p&lt;0.05).Data suggest that EYA1 may represent a risk locus for aspirin-induced GI ulceration. Though it is thought that EYA1 may play a key role in mediating cell death/ survival

  Study EGAS00001002052

• Microsatellite data show recent demographic expansions in sedentary but not in nomadic human populations in Africa and Eurasia.

  The transition from hunting and gathering to plant and animal domestication was one of the most important cultural and technological revolutions in human history. According to archeologists and paleoanthropologists, this transition triggered major demographic expansions. However, few genetic studies have found traces of Neolithic expansions in the current repartition of genetic polymorphism, pointing rather toward Paleolithic expansions. Here, we used microsatellite autosomal data to investigate the past demographic history of 87 African and Eurasian human populations with contrasted lifestyles (nomadic hunter-gatherers, semi-nomadic herders and sedentary farmers). Likely due to the combination of a higher mutation rate and the possibility to analyze several loci as independent replicates of the coalescent process, the analysis of microsatellite data allowed us to infer more recent expansions than previous genetic studies, potentially resulting from the Neolithic transition. Despite the variability in their location and environment, we found consistent expansions for all sedentary farmers, while we inferred constant population sizes for all hunter-gatherers and most herders, which could result from constraints linked to a nomadic or semi-nomadic lifestyle and/or competition for land between herders and farmers. As an exception, we inferred expansions for Central Asian herders. This might be linked with the arid environment of this area, which may have been more favorable to nomadic herders than to sedentary farmers. Alternatively, current Central Asian herders may descent from populations who have first experienced a transition from hunter-gathering to sedentary agropastoralism, and then a second transition to nomadic herding.

  Study EGAS00001000652

• Exome sequencing of uterine leiomyosarcomas

  Uterine leiomyosarcomas (ULMSs) are aggressive smooth muscle tumors associated with poor clinical outcome. Despite previous cytogenetic and molecular studies, their molecular background has remained elusive. To examine somatic variation in ULMS, we performed exome sequencing on 19 tumors. Altogether 43 genes were mutated in at least two ULMSs. Most frequently mutated genes included tumor protein P53 (TP53; 6/19; 33%), alpha thalassemia/mental retardation syndrome X-linked (ATRX; 5/19; 26%), and mediator complex subunit 12 (MED12; 4/19; 21%). Unlike ATRX mutations, both TP53 and MED12 alterations have repeatedly been associated with ULMSs. All the observed ATRX alterations were either nonsense or frameshift mutations. ATRX protein levels were reliably analyzed by immunohistochemistry in altogether 44 ULMSs, and the majority of tumors (23/44; 52%) showed clearly reduced expression. Loss of ATRX expression has been associated with alternative lengthening of telomeres (ALT), and thus the telomere length was analyzed with telomere-specific fluorescence in situ hybridization. The ALT phenotype was confirmed in all ULMSs showing diminished ATRX expression. Exome data also revealed one nonsense mutation in death-domain associated protein (DAXX), another gene previously associated with ALT, and the tumor showed ALT positivity. Aberrant expression of both TP53 and ATRX were associated with poor overall survival. In conclusion, exome sequencing revealed that TP53, ATRX, and MED12 are frequently mutated in ULMSs. ALT phenotype was commonly seen in tumors, indicating that ATR inhibitors, which were recently suggested as possible new drugs for ATRX-deficient tumors, could provide a potential novel therapeutic option for ULMS.

  Study EGAS00001001612

• Recurrent epimutations activate gene body promoters in primary glioblastoma

  Aberrant DNA hypomethylation may play an important role in the growth rate of glioblastoma (GBM), but the functional impact on transcription remains poorly understood. We assayed the GBM methylome with MeDIP-seq and MRE-seq, adjusting for copy number differences, in a small set of non-glioma CpG island methylator phenotype (non-G-CIMP) primary tumors. Recurrent hypomethylated loci were enriched within a region of chromosome 5p15 that is specified as a cancer amplicon and also encompasses TERT, encoding telomerase reverse transcriptase, which plays a critical role in tumorigenesis. Overall, 76 gene body promoters were recurrently hypomethylated, including TERT and the oncogenes GLI3 and TP73. Recurring hypomethylation also affected previously unannotated alternative promoters, and luciferase reporter assays for three of four of these promoters confirmed strong promoter activity in GBM cells. Histone H3 lysine 4 trimethylation (H3K4me3) ChIP-seq on tissue from the GBMs uncovered peaks that coincide precisely with tumor-specific decrease of DNA methylation at 200 loci, 133 of which are in gene bodies. Detailed investigation of TP73 and TERT gene body hypomethylation demonstrated increased expression of corresponding alternate transcripts, which in TP73 encodes a truncated p73 protein with oncogenic function and in TERT encodes a putative reverse transcriptase-null protein. Our findings suggest that recurring gene body promoter hypomethylation events, along with histone H3K4 trimethylation, alter the transcriptional landscape of GBM through the activation of a limited number of normally silenced promoters within gene bodies, in at least one case leading to expression of an oncogenic protein.

  Study EGAS00001000685

• DNA methylation and the adverse metabolic outcomes of adiposity

  Overweight and obesity affect ~1.5 billion people worldwide, and are major risk factors for type-2 diabetes (T2D), cardiovascular disease and related metabolic and inflammatory disturbances.1,2 Although the mechanisms linking adiposity to its clinical sequelae are poorly understood, recent studies suggest that adiposity may influence DNA methylation,3-6 a key regulator of gene expression and molecular phenotype.7 Here we use epigenome-wide association to show that body mass index (BMI, a key measure of adiposity) is associated with widespread changes in DNA methylation (187 genetic loci at P&lt;1x10-7, range P=9.2x10-8 to 6.0x10-46; N=10,261 samples). Genetic association analyses demonstrate that the alterations in DNA methylation are predominantly the consequence of adiposity, rather than the cause. We find the methylation loci are enriched for functional genomic features in multiple tissues (P&lt;0.05), and show that sentinel methylation markers identify gene expression signatures at 38 loci (P&lt;9.0x10-6, range P=5.5x10-6 to 6.1x10-35, N=1,785 samples). The methylation loci identified highlight genes involved in lipid and lipoprotein metabolism, substrate transport, and inflammatory pathways. Finally, we show that the disturbances in DNA methylation predict future type-2 diabetes (relative risk per 1SD increase in Methylation Risk Score: 2.3 [2.07-2.56]; P=1.1x10-54). Our results provide new insights into the biologic pathways influenced by adiposity, and may enable development of new strategies for prediction and prevention of type-2 diabetes and other adverse clinical consequences of obesity.

  Study EGAS00001001922

• Methylation differences in trisomy 21 using monozygotic twins

  We have hypothesized that methylation differences induced by trisomy 21 (T21) contribute to the phenotypic characteristics and heterogeneity in T21. In order to determine the methylation differences in T21 without the interference of the interindividual genomic variation, we have used fetal skin fibroblasts from monozygotic (MZ) twins discordant for T21. We also used skin fibroblasts from MZ twins concordant for T21, normal MZ twins without T21, and unrelated normal and T21 individuals. We applied Reduced Representation Bisulfite Sequencing (RRBS) to generate genome wide nucleotide resolution of DNA methylation based on high throughput sequencing between each pair of samples. RRBS revealed differentially methylated promoter regions (DMRs) in MZ twins discordant for T21 that have also been observed in unrelated normal and T21 individuals. The identified DMRs are enriched for genes involved in embryonic organ morphogenesis. These DMRs are maintained in iPS cells generated from this twin pair and are correlated with the gene expression changes. We have also observed an increase in DNA methylation level in the T21 methylome compared to the normal euploid methylome. This observation is concordant with the up regulation of DNA methyltransferase enzymes (DNMT3B and DNMT3L) and down regulation of DNA demethylation enzymes (TET2 and TET3) in the iPSC of the T21 versus normal twin. Additionally, two sets of T21 MZ twins discordant for heart defect highlighted DMRs in genes that are associated with heart development. In conclusion the study of DNA methylation differences in MZ twins discordant for genomic abnormalities is a promising approach to understand the molecular pathophysiology of aneuploidies.

  Study EGAS00001001051

• RB Loss in Resistant EGFR Mutant Lung Adenocarcinomas that Transform to Small Cell Lung Cancer

  Tyrosine kinase inhibitors (TKIs) are effective treatments for non-small cell lung cancers (NSCLCs) with epidermal growth factor receptor (EGFR) mutations. However, they do not lead to cures, and, on average, relapse occurs typically occurs afterafter an average of one year of continuous treatment. In a subset of patients, aA fundamental histological transformation from NSCLC to small cell lung cancer (SCLC) is observed in a subset of the resistant cancers, but the molecular changes associated with this transformation remain unknown. Analysis of a cohort of tumor samples and cell lines derived from resistant EGFR mutant patients with SCLC transformation revealed that RB is lost in 100% of these SCLC transformed cases, but rarely in those that remain NSCLC. Global changes in gene expression, including increased neuroendocrine marker expression and absence of EGFR expression, are observed in resistant EGFR mutant cancers that transformed to SCLC. Consistent Further, with their similarities to classical SCLC at the genetic and gene expression levelincreased neuroendocrine marker and decreased EGFR expression as well as greater sensitivity to BCL2 family inhibition are observed in resistant SCLC transformed cancers compared to resistant NSCLCs. , cell lines derived from resistant EGFR mutant SCLC biopsies are substantially more sensitive to Bcl-2 family inhibition with ABT-263 compared to those derived from resistant EGFR mutant NSCLCs. Together, these findings suggest that despite developing initially as EGFR mutant adenocarcinomas, this subset of resistant cancers ultimately adopt many of the molecular and phenotypic characteristics of classical SCLC.

  Study EGAS00001001102

• High-throughput Detection of Clinically Relevant Mutations in Archived Tumor Samples By Multiplexed PCR and Next Generation Sequencing

  PURPOSE: Tailoring cancer treatment to tumor molecular characteristics promises to make personalized medicine a reality. However, reliable genetic profiling of archived clinical specimens has been hindered by limited sensitivity and high false positive rates. Here, we describe a novel methodology, MMP-seq, which enables sensitive and specific high-throughput, high content genetic profiling in archived clinical samples. EXPERIMENTAL DESIGN: We first validated MMP-seq’s technical performance in 66 cancer cell lines and a Latin square cross-dilution of known somatic mutations. We next characterized the performance of MMP-seq in 17 formalin-fixed paraffin embedded (FFPE) clinical samples using matched fresh frozen (FF) tissue from the same tumors as benchmarks. To demonstrate the potential clinical utility of our methodology, we profiled FFPE tumor samples from 73 endometrial cancer patients. RESULTS: We demonstrated that MMP-seq enabled rapid and simultaneous profiling of a panel of 88 cancer genes in 48 samples, and detected variants at frequencies as low as 0.4%. We identified DNA degradation and deamination as the main error sources and developed practical and robust strategies for mitigating these issues, and dramatically reduced the false positive rate. Applying MMP-seq to a cohort of endometrial tumor samples identified extensive, potentially actionable alterations in the PI3K and RAS pathways, including novel PIK3R1 hotspot mutations that may disrupt negative regulation of PIK3CA. CONCLUSIONS: MMP-seq provides a robust solution for comprehensive, reliable and high- throughput genetic profiling of clinical tumor samples, paving the way for the incorporation of genomic-based testing into clinical investigation and practice.

  Study EGAS00001000674

• Exome Sequencing Of 75 Individuals From Multiply Affected Coeliac Families

  Coeliac disease is a highly heritable common autoimmune disease involving chronic small intestinal inflammation in response to dietary wheat. The human leukocyte antigen (HLA) region, and 40 newer regions identified by genome wide association studies (GWAS) and dense fine mapping, account for ~40% of the disease heritability. We hypothesized that in pedigrees with multiple individuals with coeliac disease rare [minor allele frequency (MAF) &lt;0.5%] mutations of larger effect size (odds ratios of ~ 2 – 5) might exist. We sequenced the exomes of 75 coeliac individuals of European ancestry from 55 multiply affected families. We selected interesting variants and genes for further follow up using a combination of: an assessment of shared variants between related subjects, a model-free linkage test, and gene burden tests for multiple, potentially causal, variants. We next performed highly multiplexed amplicon resequencing of all RefSeq exons from 24 candidate genes selected on the basis of the exome sequencing data in 2,248 unrelated coeliac cases and 2,230 controls. 1,335 variants with a 99.9% genotyping call rate were observed in 4,478 samples, of which 939 were present in coding regions of 24 genes (Ti/Tv 2.99). 91.7% of coding variants were rare (MAF &lt;0.5%) and 60% were novel. Gene burden tests performed on rare functional variants identified no significant associations (P&lt;1x10-3) in the resequenced candidate genes. Our strategy of sequencing multiply affected families with deep follow up of candidate genes has not identified any new coeliac disease risk mutations.

  Study EGAS00001001093

• UK10K NEURO FSZNK

  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 Finnish schizophrenia sample set has been collected from a population cohort using national registers. The entire sample collection consists of 2756 individuals from 458 families of whom 931 are diagnosed with schizophrenia spectrum disorder. Families outside Kuusamo (n=288) all had at least two affected siblings. All diagnoses are based on DSM-IV and for a large fraction of cases there is cognitive data.For further details/descriptions with regard to this data set please contact Tiina Paunio (tiina.paunio@thl.fi)

  Study EGAS00001000119

• UK10K NEURO ASD SKUSE

  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 of UK origin consists of clinically identified subjects with Autism Spectrum Disorders, mostly without intellectual disability (ie. Verbal IQs &gt;70). The subjects represent children and adults with Autism, Asperger syndrome or Atypical Autism, identified according to standardized research criteria (ADI-algorithm, ADOS). A minority has identified comorbid neurodevelopmental disorders (e.g. ADHD). Family histories are available, with measures of broader phenotype in first-degree relatives.For further information on this cohort please contact David Skuse (d.skuse@ucl.ac.uk).

  Study EGAS00001000114

• UK10K_NEURO_MUIR

  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 sample selection consists of subjects with schizophrenia (SZ), autism, or other psychoses all with mental retardation (learning disability). The samples were initially collected under the leadership of Walter J Muir (deceased), now with Prof. Blackwood, Dr McKechanie and Prof McIntosh as custodians. These subjects represent the intersection of severe forms of neurodevelopmental disorders, appear to have a higher rate of familiality of SZ than typical, and are likely to have more serious and penetrant forms of mutations.For further information on this cohort please contact Andrew McIntosh (andrew.mcintosh@ed.ac.uk).

  Study EGAS00001000122

• These are from Korean HCC samples with exome sequencing

  Hepatic resection is the most curative treatment option for early-stage hepatocellular carcinoma, but is associated with a high recurrence rate, which exceeds 50% at 5 years after surgery. Understanding the genetic basis of hepatocellular carcinoma at surgically curable stages may enable the identification of new molecular biomarkers that accurately identify patients in need of additional early therapeutic interventions. Whole exome sequencing and copy number analysis was performed on 231 hepatocellular carcinomas (72% with hepatitis B viral infection) that were classified as early-stage hepatocellular carcinomas, candidates for surgical resection. Recurrent mutations were validated by Sanger sequencing. Unsupervised genomic analyses identified an association between specific genetic aberrations and postoperative clinical outcomes. Recurrent somatic mutations were identified in 9 genes, including TP53, CTNNB1, AXIN1, RPS6KA3, and RB1. Recurrent homozygous deletions in FAM123A, RB1, and CDKN2A, and high-copy amplifications in MYC, RSPO2, CCND1, and FGF19 were detected. Pathway analyses of these genes revealed aberrations in the p53, Wnt, PIK3/Ras, cell cycle, and chromatin remodelling pathways. RB1 mutations were significantly associated with cancer-specific and recurrence-free survival after resection (p = 0.016 and p = 0.001, respectively). FGF19 amplifications, known to activate Wnt signalling, were mutually exclusive with CTNNB1 and AXIN1 mutations, and significantly associated with cirrhosis (p = 0.017). RB1 mutations can be used as a prognostic molecular biomarker for resectable hepatocellular carcinoma. Further study is required to investigate the potential role of FGF19 amplification in driving hepatocarcinogenesis in patients with liver cirrhosis and to investigate the potential of anti-FGF19 treatment in these patients.

  Study EGAS00001000604

• Chromatin immunoprecipitation linked to next-generation whole genome sequencing (ChIP-Seq) for H3K36me3 in paediatric high grade glioma cell lines KKNS4 and SF188 with and without a G34V mutation in H3F3A

  Glioblastomas of children and young adults have a median survival of only 12-15months and are clinically and biologically distinct from histologically similar cancers in older adults1. They are defined by highly specific mutations in the gene encoding the histone H3.3 variant H3F3A2, occurring either at or close to key residues marked by methylation for regulation of transcription – K27 and G34. We performed chromatin immunoprecipitation linked to next-generation whole genome sequencing (ChIP-Seq) for H3K36me3 in order to test the hypothesis that, rather than total H3K36me3, the G34V mutation may instead result in differential binding of the trimethyl mark throughout the genome. Compared with H3F3A wild-type SF188 paediatric GBM cells, H3K36me3 was found to be significantly differentially bound in KNS42 cells at 5130 distinct regions of the genome corresponding to 156 genes. Concurrent whole genome DNA sequencing demonstrated that this was not confounded by copy number alterations. The cerebral hemispheric-specific G34 mutation drives a distinct expression signature through differential genomic binding of the K36 trimethylation mark (H3K36me3). The transcriptional program induced recapitulates that of the developing forebrain, and involves numerous markers of stem cell maintenance, cell fate decisions and self-renewal. Critically, H3F3A G34 mutations cause profound upregulation of MYCN, a potent oncogene which is causative of glioblastomas when expressed in the correct developmental context. This driving aberration is selectively targetable in this patient population by inhibiting kinases responsible for stabilisation of the protein.

  Study EGAS00001001437

• Spatial and temporal evolution of distal 10q deletion, a prognostically unfavorable event in diffuse low-grade gliomas

  Background: The disease course of patients with diffuse low-grade glioma is notoriously unpredictable. Temporal and spatially distinct samples may provide insight into the evolution of clinically relevant copy number aberrations (CNAs). The purpose of this study is to identify CNAs that are indicative of aggressive tumor behaviour and can thereby complement the prognostically favorable 1p/19q co-deletion. Results: Genome-wide, 50 base pair single-end, sequencing was performed to detect CNAs in a clinically well-characterized cohort of 98 formalin-fixed paraffin-embedded low-grade gliomas. CNAs are correlated with overall survival as an endpoint. Seventy-five additional samples from spatially distinct regions and paired recurrent tumors of the discovery cohort were analysed to interrogate the intratumoral heterogeneity and spatial evolution. Loss of 10q25.2-qter is a frequent subclonal event and significantly correlates with an unfavorable prognosis. A significant correlation is furthermore observed in a validation set of 126 and confirmation set of 184 patients. Loss of 10q25.2-qter arises in a longitudinal manner in paired recurrent tumor specimens, whereas the prognostically favorable 1p/ 19q co-deletion is the only CNA that is stable across spatial regions and recurrent tumors. Conclusions: CNAs in low-grade gliomas display extensive intratumoral heterogeneity. Distal loss of 10q is a late onset event and a marker for reduced overall survival in low-grade glioma patients. Intratumoral heterogeneity and higher frequencies of distal 10q loss in recurrences suggest this event is involved in outgrowth to the recurrent tumor.

  Study EGAS00001000643

• Study on the consequences of prenatal famine exposure on DNA methylation.

  The potential of the epigenome to undergo environmentally-induced changes during early development is central to its postulated role in human disease. Here, we report a genome-scale view of DNA methylation differences after early gestational malnutrition that was caused by the Dutch Hunger Winter, a severe famine at the end of World War II. We adopted an annotation guided analysis of reduced-representation bisulfite sequencing data generated in 24 middle-aged individuals who had been exposed to famine in utero and 24 unexposed sibling controls. We observed differential DNA methylation primarily at genomic regions with regulatory potential, including enhancers active during early development. In-depth annotation of 181 prenatally-induced differentially methylated regions (P-DMRs) showed that they were predominantly located in gene bodies and co-occurred with histone marks denoting an active state. In line with the early gestational timing of the exposure, the P-DMRs mapped to genes that were enriched for differential expression patterns during blastocyst development and early organogenesis. Validation and further explorative analyses of 6 P-DMRs mapping to SMAD7, CDH23, INSR, RFTN1, CPT1A and KLF13 highlighted the critical role of gestational timing, indicated that differential methylation extended along pathways involved in growth, development and metabolism, and suggested associations with birth weight and serum LDL cholesterol in adulthood and confirmed an enhancer function for the INSR and CPT1A P-DMRs. Our results point toward a link between prenatal malnutrition and epigenetic modulation of growth and lipid metabolism that may underlie the adverse metabolic phenotype of exposed individuals in later life.

  Study EGAS00001000668

• Dilgom_Exome

  The Finrisk sample sets are part of the National FINRISK Study. It is a large population survey on risk factors of chronic, noncommunicable diseases. The survey is carried out since 1972 every five years using independent, random and representative population samples from different parts of Finland. The main results from the previous FINRISK 2007 survey are published.The National FINRISK Study Survey was carried out in 5 areas in Finland and 2000 inhabitants aged 25-75 years were invited to participate in each year. Among findings were that Finns continue to gain weight.Data from FINRISK surveys are used for many different research projects and for national health monitoring needs. The recent research activities deal, in addition to cardiovascular diseases and the classical risk factors, also with e.g. asthma and allergy, alcohol, socioeconomic factors and genetic epidemiology.The FINRISK study is part of the MORGAM Project (MONICA Risk, Genetics, Archiving and Monigraph), sponsored by the EU and MDECODE (Molecular Diversity and Epidemiology of Common Disease) program coordinated by the University of Michigan.The exome sequencing study will be part of the Dilgom study, which is a part of the larger Finrisk population based health study performed in Finland. It consists of 5000 individuals with a prospective aspect of metabolic traits. The cohort has been extensively phenotyped for their cardiovascular and metabolic status. So far, we have performed a 660K Illumina GWAS and a full genome wide expression study from peripheral blood cells of 500 individuals. The cohort has also been in total genotyped by the cardiometabochip

  Study EGAS00001000086

• Quantification of chromosomal copy number aberrations by shallow whole-genome sequencing

  Detection of DNA copy number aberrations by shallow whole-genome sequencing (WGS) faces many challenges including lack of completion and errors in the human reference genome, repetitive sequences, polymorphisms, variable sample quality, and biases in the sequencing procedures.Formalin-fixed paraffin-embedded (FFPE) archival material, the analysis of which is important for studies of cancer, presents particular analytical difficulties due to degradation of the DNA and frequent lack of matched reference samples. We present a robust, cost-effective WGS method for DNA copy number analysis that addresses these challenges more successfully than currently available procedures. In practice very useful profiles can be obtained with 0.1 fold genome coverage. We improve on previous methods by; first, implementing a combined correction for sequence mappability and GC content, and second, applying this procedure to sequence data from the 1000 Genomes Project in order to develop a blacklist of problematic genome regions. A small subset of these blacklisted regions were previously identified by ENCODE, but the vast majority are novel unappreciated problematic regions. Our procedures are implemented in a pipeline called QDNAseq. We have analyzed over 1,000 samples, most of which were obtained from the fixed tissue archives of over 25 institutions.We demonstrate that for most samples our sequencing and analysis procedures yield genome profiles with noise levels near the statistical limit imposed by read counting. The described procedures also provide better correction of artifacts introduced by low DNA quality than prior approaches, and better copy number data than high-resolution microarrays at substantially lower cost.

  Study EGAS00001000642

• Prediction and quantification of splice events from RNA-seq data

  Analysis of splice variants from short read RNA-seq data remains a challenging problem. Here we present a novel method for the genome-guided prediction and quantification of splice events from RNA-seq data, which enables the analysis of unannotated and complex splice events. Splice junctions and exons are predicted from reads mapped to a reference genome and are assembled into a genome-wide splice graph. Splice events are identified recursively from the graph and are quantified locally based on reads extending across the start or end of each splice variant. We assess prediction accuracy based on simulated and real RNA-seq data, and illustrate how different read aligners (GSNAP, HISAT2, STAR, TopHat2) affect prediction results. We validate our approach for quantification based on simulated data, and compare local estimates of relative splice variant usage with those from other methods (MISO, Cufflinks) based on simulated and real RNA-seq data. In a proof-of-concept study of splice variants in 16 normal human tissues (Illumina Body Map 2.0) we identify 249 internal exons that belong to known genes but are not related to annotated exons. Using independent RNA samples from 14 matched normal human tissues, we validate 9/9 of these exons by RT-PCR and 216/249 by paired-end RNA-seq (2 x 250 bp). These results indicate that de novo prediction of splice variants remains beneficial even in well-studied systems. An implementation of our method is freely available as an R/Bioconductor package SGSeq.

  Study EGAS00001001026

• Prediction of pigmentation phenotypes by SNP typing in a Northern German population

  Human pigmentation traits are of great interest to many research areas, from ancient DNA analysis to forensic science. We aimed to develop a gene-based predictive model for pigmentation phenotypes in a realistic target population for forensic case work from Northern Germany. Our aim was to determine whether better prediction accuracy can be achieved, or fewer genetic markers may suffice, than in previously studied, genetically more heterogeneous populations. We investigated the association between eye, hair and skin colour, and 12 candidate single nucleotide polymorphisms (SNPs) from six genes. Our study comprised two samples of 300 and 100 individuals from Northern Germany who were carefully characterized with regard to pigmentation phenotypes. The first sample was used to select trait-associated SNPs whereas the second sample served to estimate odds ratios (ORs) and to quantify the predictive capability of the respective SNP genotypes. SNP rs12913832 in HERC2 was found to be strongly associated with blue eye colour (OR=15.6, p&lt;1.2•10-4) and to yield reasonable predictive power (90% sensitivity, 63% specificity). SNP associations with hair and skin colour were weaker and genotypes less predictive. A comparison to two recently published sets of markers to predict eye and hair colour revealed that the consideration of additional SNPs with weak to moderate effect increases the predictive power in Northern Germans for eye colour, but not for hair colour. In addition, fine phenotyping and differentiation of hair colour (light / dark and red tint / no red tint) were found to increase the number of significant genotype-phenotype associations.

  Study EGAS00001001174

• Associations between APOE status and cognitive ability across the lifecourse

  The association between APOE genotype and cognitive function suggests a positive role for the e2 allele and a negative role for the e4 allele. Both alleles have relatively low frequencies in the general population, hence meta-analyses have been based on many small, heterogeneous studies. Here we report the APOE-cognition associations in the largest single analysis to date. APOE status and cognitive ability were measured in 18,337 participants from the Generation Scotland study between 2006 and 2011. The age range was 18-94 years with a mean of 47 (SD 15). Four cognitive domains were assessed: verbal declarative memory (paragraph recall), processing speed (digit symbol substitution), verbal fluency (phonemic verbal fluency), and vocabulary (Mill Hill synonyms). Linear regression was used to assess the associations between APOE genetic status and cognition. Possession of the e4 allele was associated with lower scores on the measures of memory and processing speed in subjects aged &gt;60. Across all age ranges, the e4 allele was linked to better verbal fluency scores. In younger subjects (≤60 years) the e4 allele was linked to higher vocabulary scores. There were no associations between the e2 allele and cognitive ability. As seen in previous meta-analyses, the APOE e4 allele is linked to poorer cognitive performance in the domains of memory and processing speed. By contrast, positive associations were seen between the e4 allele and measures of verbal fluency and vocabulary. All associations were relatively small and, in many cases, nominally significant despite the very large sample size.

  Study EGAS00001001235

• Search for new loci and low-frequency variants influencing glioma risk by exome-array analysis

  To identify protein altering variants (PAVs) for glioma we analysed Illumina HumanExome BeadChip exome array data on 1,882 glioma cases and 8,079 controls from three independent European populations. In addition to single variant tests we incorporated information on the predicted functional consequences of PAVs and analysed sets of genes with a higher likelihood of having a role in glioma on the basis of the profile of somatic mutations documented by large-scale sequencing initiatives. Globally tThere was a strong relationship between effect size and SNPs predicted to be damaging (P=2.29x10-49); however, these variants which are most likely to impact on risk, are rare (MAF&lt;5%). While no single variant showed an association which was statistically significant at the genomewide threshold a number of represented promising Notable associations - were seen with the BRCA2:c.9976A&gt;T, p.(Lys3326Ter)BRCA2 p.Lys332X variant, which has been shown documented to influence breast and lung cancer risks (OR=2.3, P=4.00x10-4 for glioblastoma [GBM]) and IDH2:c.782G&gt;A, p.(Arg261His)IDH2 p.Arg261His (OR=3.21, P=7.67x10-3, for non-GBM). Additionally, gene burden tests revealed a statistically significant association for HARS2 and risk of GBM (P=2.20x10-6). Genome scans of low frequency PAVs represent a complementary strategy to identify disease-causing variants compared with scans based on tagSNPs. Strategies to lessen the multiple testing burden by restricting analysis to PAVs with higher priors affords an opportunity to maximise study power.

  Study EGAS00001001258

• Transcriptome_human_nasal_epithelium

  1. Odors are detected, firstly, by olfactory sensory neurons (OSNs) in the olfactory epithelium of the nose. This neurons then project directly to the olfactory bulb in the brain. Olfaction depends on cellular regeneration of the OE, olfactory bulb and hippocampus, and on their continual re-wiring. The olfactory neural pathway includes regions of the frontal, temporal and limbic brain, which in turn overlap with brain areas involved in brain disorders. OSNs are the only aspect of the human brain exposed to the external environment. This not only makes them vulnerable to environmental changes, but also accessible for biomedical studies.We have already sequenced and developed a protocol for analyzing the transcriptome of mouse main olfactory epithelium and single OSNs. We propose here to perform a similar study for samples from the human olfactory epithelium. We have developed a minimally invasive method for obtaining human OSNs, among other cells from the nasal epithelium. In this experiment, we have obtained cell samples from the olfactory epithelium, including OSN, from healthy volunteers. We would like to further characterize them by RNA sequencing. This will give us valuable insight into human olfaction. It will also provide a first step into a new avenue to study, and find biomarkers for, brain diseases though the analysis of these easily available neurons. 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 EGAS00001001294

• UK10K_RARE_THYROID

  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.Two cohorts of subjects are being analysed: Individuals with Congenital Hypothyroidism (CH) due either to dysgenesis or dyshormonogenesis; Patients with Resistance to Thyroid hormone (RTH), a disorder characterized by elevated thyroid hormones and variable tissue refractoriness to hormone action. The CH cohort has been enriched for genetic aetiologies by recruiting cases that are familial, on a consanguineous background or syndromic. The RTH cohort consists of cases in which candidate gene analyses have been negative.For further information with regard to this cohort please contact Krishna Chatterjee (kkc1@medschl.cam.ac.uk).

  Study EGAS00001000131

• A unidirectional histone code in bidirectional promoters across cell types

  Recent studies indicate that thousands of genes are expressed from bidirectional promoters (BDPs). Gene regulation at BDPs is poorly understood, in particular how the cell is able to regulate them differently. Here we investigate the effect of histone Modifications in BDP regulation. In this study, we model the gene activity using different gene expression assays, such as RNA-Seq, GRO-cap, and CAGE around the BDP transcription start sites (TSSs) using different histone modifications in various cell types. We develop a new statistical approach that links histone modifications to gene expression at BDPs. It improves over previous methods, because it is able to capture spatial dependencies of histone modifications along a promoter and leads to more interpretable results. We predict a general histone code that is independent of transcript orientation, cell type, and promoter configuration. The histone code at BDPs reveals that promoters are regulated unidirectionally, such that the majority of histone marks associated with gene expression occur downstream of the gene's TSS. By contrasting associations of histone marks with steady-state levels of capped RNAs in CAGE and nascent capped RNAs in GRO-cap, we show which histone marks have a preference for initiating polymerases and actively elongating polymerases. Using single-cell data we show that the bidirectional histone signal of activating marks is often due averaging of a heterogeneous cell population. Our results have implications for other experiments where the relationship between histone modification and gene initiation or gene elongation is investigated on a genome-wide scale.

  Study EGAS00001001656

• Evolution of a FLT3-TKD mutated subclone at meningeal relapse in acute promyelocytic leukemia (H021)

  Here we report the case of an acute promyelocytic leukemia (APL) patient who, though negative for FLT3 mutations at diagnosis, developed isolated FLT3-TKD positive meningeal relapse, which, in retrospect, could be traced back to a minute bone marrow subclone present at first diagnosis. Initially, the 48-year old female diagnosed with high-risk APL had achieved complete molecular remission after standard treatment with all-trans retinoic acid (ATRA) and chemotherapy according to the AIDA protocol. 13 months after start of ATRA maintenance the patient suffered clinically overt meningeal relapse along with minute molecular traces of PML/RARA in the bone marrow. Following treatment with arsenic trioxide (ATO) and ATRA in combination with intrathecal cytarabine and methotrexate, the patient achieved a complete molecular remission in both cerebrospinal fluid (CSF) and bone marrow, which currently lasts for two years after completion of therapy. Whole exome sequencing and subsequent ultradeep targeted re-sequencing revealed a heterozygous FLT3-TKD mutation in CSF leukemic cells (p.D835Y, c.2503G&gt;T, 1000/1961 reads (51%)), which was undetectable in the concurrent bone marrow sample. Interestingly, the FLT3-TKD mutated meningeal clone originated from a small bone marrow subclone present in a variant allele frequency of 0.4% (6/1553 reads) at initial diagnosis. This case highlights the concept of clonal evolution with a subclone harboring an additional mutation being selected as the “fittest” and leading to meningeal relapse. It also further supports earlier suggestions that FLT3 mutations may play a role for migration and clonal expansion in the CSF sanctuary site.H021

  Study EGAS00001001848

• Integrated clinical, whole genome, and transcriptome analysis of multisampled lethal metastatic prostate cancer

  We report the first combined analysis of whole genome sequence, detailed clinical history, and transcriptome sequence of multiple prostate cancer metastases in a single patient (A21). Whole genome and transcriptome sequence was obtained from 9 anatomically separate metastases, and targeted DNA sequencing was performed in cancerous and noncancerous foci within the primary tumor specimen removed 5 years prior to death. Transcriptome analysis revealed increased expression of AR-regulated genes in liver metastases that harbored an AR p.L702H mutation, suggesting a dominant effect by the mutation despite being present in only 1 of an estimated 16 copies per cell. The metastases harbored several alterations to the PI3K/AKT pathway, including a clonal truncal mutation in PIK3CG and present in all metastatic sites studied. The list of truncal genomic alterations shared by all metastases included homozygous deletion of TP53, hemizygous deletion of RB1 and CHD1, and amplification of FGFR1. If the patient were treated today given this knowledge, use of second-generation androgen-directed therapies, cessation of glucocorticoid administration, and therapeutic inhibition of the PI3K/AKT pathway or FGFR1 receptor could provide personalized benefit. Three previously unreported truncal clonal missense mutations (ABCC4 p.R891L, ALDH9A1 p.W89R, and ASNA1 p.P75R) were expressed at the RNA level and assessed as druggable. The truncal status of mutations is critical for actionability, and can only be determined through analysis of multiple sites of metastasis. Our findings suggest that a large set of deeply analyzed cases could serve as powerful guide to more effective prostate cancer basic science and personalized cancer medicine clinical trials.

  Study EGAS00001001659

• Transcriptome Sequencing PPGL

  RNA library preparation of samples with RNA Integrity Number (RIN) &gt; 5.5 (median 7.1; range 5.6-8.9) was performed as described in the TruSeq Stranded mRNA Library Prep Kit (Illumina, RS-122-2101). Briefly, 0.7-1 µg total RNA PolyA+ fraction was purified and randomly fragmented, converted to double stranded cDNA and processed through subsequent enzymatic treatments of end-repair, dA-tailing, and ligation to adapters. Adapter-ligated library was completed by PCR with Illumina PE primers. The resulting purified cDNA library was applied to an Illumina flow cell for cluster generation and sequenced on an Illumina HiSeq2500 on a 51bp single-read format following manufacturer's protocols. cDNA libraries from FFPE tumors and low integrity RNA’s (RIN&lt;5.5) from FF tumors (median 2.4; range 1-5.5) but % of fragments &gt; 200 nucleotides higher than 20% (median 57.5; range 20-96) were prepared using QuantSeq 3’ mRNA-Seq Library Prep Kit FWD for Illumina (Lexogen, 015) with a UMI Second Strand Synthesis Module for QuantSeq FWD (Lexogen, 081), vendor’s protocol for low input/ low quality/ FFPE RNA using 200-500 ng of total RNA, according to availability, was used as starting material. PCR Add-on Kit for Illumina (Lexogen, 020) was used to adjust the library amplification number of cycles. Libraries were applied to an Illumina flow cell for cluster generation and sequenced on NovaSeq6000. Image analysis, per-cycle basecalling and quality score assignment was performed with Illumina Real Time Analysis software. Demultiplexing of BCL files to FASTQ format was performed with the bcl2fastq Software (Illumina).

  Study EGAS00001006044

• Non-viral precision T cell receptor replacement for personalized cell therapy

  The T cell receptor (TCR) provides the fine specificity of T cells to recognize mutations in cancer cells 1-3. We developed a clinical-grade approach based on CRISPR/Cas9 non-viral precision genome editing to simultaneously knock-out the two endogenous TCR genes, TCRα (TRAC) and TCRβ (TRBC), and insert in the TRAC locus the two chains of a neoantigen-specific TCR (neoTCR), isolated from the patient’s own circulating T cells using a personalized library of soluble predicted neoantigen-HLA capture reagents. Sixteen patients with refractory solid cancers received up to three distinct neoTCR-transgenic cell products, each expressing a patient-specific neoTCR, in a cell dose-escalation, first-in-human phase 1 clinical trial (NCT03970382). One patient had grade 1 cytokine release syndrome, and one grade 3 encephalitis. All had the expected side effects from the lymphodepleting chemotherapy. Five patients had stable disease, and the other 11 had disease progression as best response on therapy. NeoTCR-transgenic T cells were detected in tumour biopsies post-infusion at frequencies higher than the native TCRs pre-infusion. This study demonstrates the feasibility of isolating and cloning multiple TCRs recognizing mutational neoantigens, the simultaneous knock-out of the endogenous TCR and knock-in of the neoTCRs using single-step, non-viral precision genome editing, the manufacturing of neoTCR engineered T cells at clinical grade, the safety of infusing up to three gene edited neoTCR T cell products, and the ability of the transgenic T cells to traffic to the patients’ tumours.

  Study EGAS00001006898