22928 results for "ega"

in 84.77 milliseconds.

• Bacterial SNPs in the human gut microbiome associate with host BMI

  Genome-wide association studies (GWAS) have revolutionized the field of genetics, providing numerous associations between human SNPs and health traits. Likewise, metagenome-wide association studies (MWAS) between bacterial SNPs and human traits can suggest mechanistic links between the microbiota and its host. However, very few such studies have been done to date, in part due to the large number of samples required to address the variable coverage of bacteria across samples and the fact that most bacteria are present in only a subset of the population. Here, we devised an MWAS framework to detect SNPs and associate them with host phenotypes systematically. We recruited and obtained gut metagenomic samples from a cohort of 7,190 healthy individuals and discovered 1,358 statistically significant associations between a bacterial SNP and host body mass index (BMI). To address the possible effects of population structure and linkage disequilibrium, we applied a clumping procedure, and found 40 independent associations between SNPs and host BMI. We uncovered BMI-associated SNPs in the genomes of 27 bacterial species, even though the relative abundance of 12 of these species was not associated with the phenotype. We separated genome-wide sub-species variations from local associations in individual SNPs, and demonstrated how this framework can highlight specific bacterial functions of interest. Our results demonstrate the importance of considering nucleotide-level diversity in microbiome studies and pave the way toward a better understanding of interpersonal gut microbiome differences and their health implications.

  Study EGAS00001007204

• Atypical B cells and impaired SARS-CoV-2 neutralisation following heterologous vaccination in the elderly

  Age is a major risk factor for hospitalization and death after SARS-CoV-2 infection, even in vaccinees. Suboptimal responses to a primary vaccination course have been reported in the elderly, but there is little information regarding the impact of age on responses to booster third doses. Here we show that individuals 70 or older who received a primary two dose schedule with AZD1222 and booster third dose with mRNA vaccine achieved significantly lower neutralizing antibody responses against SARS-CoV-2 spike pseudotyped virus compared to those below 70 at one month post booster. Neither the concentration of serum binding anti spike IgG antibody, nor the frequency of spike- specific B cells showed differences by age grouping. Impaired neutralization potency and breadth post-third dose in the elderly was associated with enrichment of circulating “atypical” spike-specific B cells expressing CD11c and FCRL5. Presence of these atypical B cells was confirmed by scRNA seq. However, when considering individuals who received three doses of mRNA vaccine, we did not observe statistically significant differences in neutralization or enrichment in atypical B cells by age grouping following dose 3. Together, our data reveal impaired neutralising antibody responses in the elderly after heterologous vaccination with AdV followed by mRNA booster dose, but not in homologous three dose mRNA vaccination highlighting that different vaccine formats differentially instruct the memory B cell response.

  Study EGAS00001007385

• Oxidative phosphorylation is a key ontogenetic feature of monocyte immunometabolism promoting myeloid differentiation after birth

  Neonates primarily rely on innate immunity, yet their inflammatory responses to microbes, particularly those of monocytes, are usually restricted compared to adults. This is controversially interpreted as immaturity increasing the risk of sepsis or essential programming allowing immune adaptation after birth. Which changes the cellular immunometabolism undergoes after birth and its role in these concepts are poorly defined. Here, we applied transcriptomic, metabolic, and immunological approaches and found that monocytes exhibit an inverse ontogenetic balance of immunometabolism. Our data show that glycolysis in monocytes increases within the first year of life and fuels their inflammatory responsiveness. In contrast, oxidative phosphorylation (OXPHOS) is high in neonatal monocytes supporting myeloid differentiation but declines only gradually during the first five years of life. Treatment of neonatal monocytes with lipopolysaccharide induces an adult-like immunometabolic phenotype. Ketogenic diet restricts glycolysis in adult monocytes, but cannot reactivate OXPHOS and revive a neonatal-like immunometabolic phenotype, suggesting that neonatal metabolism is hardwired and cannot be restored by simple dietary changes. Transcriptional network and population-wide human variation analyses identified E2F1, MYB, STAT1 and FLI1 as important regulators of age-dependent cell functions and related energy demands in human monocytes. Collectively, our findings show that restricted glycolysis and increased OXPHOS are a physiological programming in neonatal monocytes. Premature switching to an adult-like metabolism could untimely enhance inflammatory responses and disrupt important myeloid differentiation, whereas microbial challenges accumulating during childhood seem to induce an essential gradual metabolic switch.

  Study EGAS00001007555

• Low-coverage whole genome sequencing for a highly selective cohort of severe COVID-19 patients

  Background Despite advances in identification of genetic markers associated to severe COVID-19, the full genetic characterisation of the disease remains elusive. Imputation of low-coverage whole genome sequencing (lcWGS) has emerged as a competitive method to study such disease-related genetic markers as they enable genotyping of most common genetic variants used for genome wide association studies. This study aims at exploring the potential use of imputation in lcWGS for a highly selected severe COVID-19 patient cohort. Findings We generated an imputed dataset of 79 variant call format (VCF) patient files using the GLIMPSE1 tool, each containing, on average, 9.5 million single nucleotide variants. The validation assessment of imputation accuracy yielded a squared Pearson correlation of approximately 0.97 across sequencing platforms, showing that GLIMPSE1 can be used to confidently impute variants with minor allele frequency up to approximately 2% in Spanish ancestry individuals. We conducted a comprehensive analysis on the patient cohort, examining hospitalisation and intensive care utilisation, sex and age-based differences, and clinical phenotypes using a standardised set of medical terms specifically developed to characterise severe COVID-19 symptoms for this cohort. Conclusion This dataset highlights the utility and accuracy of lcWGS imputation in the study of COVID-19 severity, setting a precedent for other applications in resource-constrained environments. The methods and findings presented here may be leveraged in future genomic projects, providing vital insights for health challenges like COVID-19.

  Study EGAS00001007573

• Epigenetic age deceleration reflects fitness improvements following a six-month endurance exercise intervention

  Epigenetic clocks are emerging as promising biomarkers of biological aging, yet their sensitivity to short-term interventions remains unclear. This pilot study investigates whether the GrimAge clock can capture the effects of a six-month cycling-based endurance exercise training intervention, with cardiorespiratory fitness (VO₂ max) and body composition as primary outcomes. We enrolled 42 adults aged 35-65, of whom 38 completed the study and 33 adhered to the protocol (>66% adherence). Participants demonstrated significant improvements in VO₂ max (+20%, P < 0.001) and body composition (P < 0.001). High-quality epigenetic data preprocessing yielded highly reproducible GrimAge estimates (< 2 months measurement error), which strongly correlated with chronological age (R² = 0.86, P < 0.001). On average, GrimAge decreased by 7.44 months relative to the expected trajectory (P = 0.012), reflecting improvements in VO₂ max (R² = 0.27, P = 0.002) but not body composition changes. Notably, GrimAge changes strongly correlated with fluctuations in leukocyte composition, particularly neutrophil fraction (R² = 0.74, P < 0.001). Adjusting for leukocyte composition improved consistency in GrimAge changes, aligning them with additional intervention outcomes and explaining up to 81% of variance. These findings demonstrate that GrimAge is responsive to short-term endurance training, serving as a meaningful biomarker of improved cardiorespiratory fitness, while also capturing immune system variability. This study supports the use of GrimAge in evaluating longevity interventions and highlights the importance of accounting for leukocyte composition in epigenetic aging research.

  Study EGAS00001008221

• The long term effects of chemotherapy on normal blood - WGS dataset

  In developed countries, ~10% of individuals are exposed to systemic chemotherapy for cancer and other diseases. Many chemotherapeutic agents act by increasing DNA damage in cancer cells, hence triggering cell death. However, there is limited understanding of the extent and consequences of collateral DNA damage to normal tissues. To investigate the impact of chemotherapy on mutation burdens and cell population structure of a normal tissue we sequenced blood cell genomes from 23 individuals, aged 3-80 years, treated with a range of chemotherapy regimens. Substantial additional mutation loads with characteristic mutational signatures were imposed by some chemotherapeutic agents, but there were differences in burden between different classes of agent, different agents of the same class and different blood cell types. Chemotherapy also induced premature changes in the cell population structure of normal blood, similar to those of normal ageing. The results constitute an initial survey of the long-term biological consequences of cytotoxic agents to which a substantial fraction of the population is exposed during the course of their disease management, raising mechanistic questions and highlighting opportunities for mitigation of adverse effects.

  Dataset EGAD00001015339

• cfMethyl-seq data (cfSort study) from serial plasma samples of NSCLC patients

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

  Dataset EGAD00001010881

• NIPT samples for systematic evaluation of BinDel, a software tool for detecting clinically significant microdeletions in low-coverage WGS-based NIPT samples

  This dataset contains three sets of samples. The first sample set contains euploid fetus pregnancies reported by NIPTIFY screening test and postnatal evaluation. Dataset was processed similarly to previously published guidelines from KU Leuven, with modifications [1]. Briefly, peripheral blood samples were collected in cell-free DNA BCT tubes (Streck, USA), and plasma was separated with standard dual centrifugation. Cell-free DNA was extracted from 3 ml plasma using MagMAX Cell-Free DNA Isolation Kit (ThermoFisher Scientific). Whole-genome libraries were prepared using the FOCUS (Fragmented DNA Compact Sequencing Assay, Competence Centre on Health Technologies, Estonia) NIPT method protocol with 12 cycles for the final PCR enrichment step. In the following quantification, equal amounts of 36 samples were pooled, and the quality and quantity of the pool were assessed on Agilent 2200 TapeStation (Agilent Technologies, USA). Whole genome sequencing was performed on the NextSeq 550 instrument (Illumina Inc.) with an average coverage of 0.32× (minimum 0.08 and maximum 0.42) and producing 85 bp single-end reads. The second sample set contains a single NIPT sample postnatally diagnosed with Prader-Willi syndrome. The sample was sequenced with Illumina NextSeq 500 platform, producing 85 bp single-end reads with an average per-sample coverage of 0.32× at the University of Tartu, Institute of Genomics Core Facility, according to the manufacturer’s standard protocols, as described previously [2]. The third sample set contains samples SC005 (SeraCare Life Sciences Inc lot #10446565), SC0042 (#10571706), and SC016 (#10560229). These are SeraCare Life Sciences Inc circulating cell-free DNA (ccfDNA) like mixture of human genomic DNA that consists of matched maternal and fetus. SC005 and SC0042 consist of matched DNA of maternal and fetus with DiGeorge Syndrome. SC016 is a custom-ordered DNA Mix with fetus DNA having a pathogenic loss of the terminal region of 20p13 and a pathogenic 3q29 duplication. SC016 was processed as the first sample set was processed, and SC0042 was processed as the second sample set was processed. Sample SC005 was processed once as was sample set 1 and once as was sample set 2 processed. This study was performed with the approval of the Research Ethics Committee of the University of Tartu (#352/M-12). 1. Bayindir B, Dehaspe L, Brison N, Brady P, Ardui S, Kammoun M, et al. Noninvasive prenatal testing using a novel analysis pipeline to screen for all autosomal fetal aneuploidies improves pregnancy management. Eur J Hum Genet. 2015;23: 1286– 1293. doi:10.1038/ejhg.2014.282 2. Žilina O, Rekker K, Kaplinski L, Sauk M, Paluoja P, Teder H, et al. Creating basis for introducing noninvasive prenatal testing in the Estonian public health setting. Prenat Diagn. 2019;39: 1262–1268. doi:10.1002/pd.5578

  Dataset EGAD00001009512

• Genome-wide cell-free DNA termini in patients with cancer

  The structure, fragmentation pattern, length and terminal sequence of cell-free DNA (cfDNA) is under the influence of nucleases present in the blood. We hypothesized that differences in the diversity of bases at the end of cfDNA fragments can be leveraged on a genome-wide scale to enhance the sensitivity for detecting the presence of tumor signals in plasma. We surveyed the cfDNA termini in 72 plasma samples from 319 patients with 18 different cancer types using low-coverage whole genome sequencing. The fragment-end sequence and diversity were altered in all cancer types in comparison to 76 healthy controls. We converted the fragment end sequences into a quantitative metric and observed that this correlates with circulating tumor DNA tumor fraction (R = 0.58, p < 0.001, Spearman). Using these metrics, we were able to classify cancer samples from control at a low tumor content (AUROC of 91% at 1% tumor fraction) and shallow sequencing coverage (mean AUROC = 0.99 at >1M fragments). Combining fragment-end sequences and diversity using machine learning, we classified cancer from healthy controls (mean AUROC = 0.99, SD = 0.01). Using unsupervised clustering we showed that early-stage lung cancer can be classified from control or later stages based on fragment-end sequences. We observed that fragment-end sequences can be used for prognostication (hazard ratio: 0.49) and residual disease detection inresectable esophageal adenocarcinoma patients, moving fragmentomics toward a greater clinical implementation.

  Study EGAS00001005738

• Genome-wide cell-free DNA termini in patients with cancer

  The structure, fragmentation pattern, length and terminal sequence of cell-free DNA (cfDNA) is under the influence of nucleases present in the blood. We hypothesized that differences in the diversity of bases at the end of cfDNA fragments can be leveraged on a genome-wide scale to enhance the sensitivity for detecting the presence of tumor signals in plasma. We surveyed the cfDNA termini in 72 plasma samples from 319 patients with 18 different cancer types using low-coverage whole genome sequencing. The fragment-end sequence and diversity were altered in all cancer types in comparison to 76 healthy controls. We converted the fragment end sequences into a quantitative metric and observed that this correlates with circulating tumor DNA tumor fraction (R = 0.58, p < 0.001, Spearman). Using these metrics, we were able to classify cancer samples from control at a low tumor content (AUROC of 91% at 1% tumor fraction) and shallow sequencing coverage (mean AUROC = 0.99 at >1M fragments). Combining fragment-end sequences and diversity using machine learning, we classified cancer from healthy controls (mean AUROC = 0.99, SD = 0.01). Using unsupervised clustering we showed that early-stage lung cancer can be classified from control or later stages based on fragment-end sequences. We observed that fragment-end sequences can be used for prognostication (hazard ratio: 0.49) and residual disease detection inresectable esophageal adenocarcinoma patients, moving fragmentomics toward a greater clinical implementation.

  Study EGAS00001005747

• Differentially methylated CpGs in response to growth hormone administration in children with idiopathic short stature

  Background: Recombinant human growth hormone (rhGH) has shown a great growth-promoting potential in children with idiopathic short stature (ISS). However, the response to rhGH differs across individuals, largely due to genetic and epigenetic heterogeneity. Since epigenetic marks on the methylome can be dynamically influenced by GH, we performed a comprehensive pharmaco-epigenomics analysis of DNA methylation changes associated with long-term rhGH administration in children with ISS. Results: We measured DNA methylation profiles before and after GH treatment (with a duration of ~18 months in average) on 47 healthy children using customized methylC-seq capture sequencing. Their changes were compared and associated with changes in plasma IGF1 by adjusting sex, age, treatment duration, and estimated blood proportions. We observed a considerable inter-individual heterogeneity of DNA methylation changes responding to GH treatment. We identified 267 response-associated differentially methylated cytosines (DMC) that were enriched in promoter regions, CpG islands and blood cell-type specific regulatory elements. Furthermore, the genes associated with these DMCs were enriched in the biology process of "cell development", "neuron differentiation" and "developmental growth", and in the TGF-beta signaling pathway, PPAR Alpha pathway, endoderm differentiation pathway, adipocytokine signaling pathway as well as PI3K-Akt signaling pathway, and cAMP signaling pathway. Conclusion: Our study provides a first insight in DNA methylation changes associated with rhGH administration, which may help understand mechanisms of epigenetic regulation on GH-responsive genes.

  Study EGAS00001006115

• Non-muscle Invasive Bladder Cancer Molecular Subtypes Predict Differential Response to Intravesical Bacillus Calmette-Guérin

  The recommended treatment for patients with high-risk non-muscle invasive bladder cancer (HR-NMIBC) is tumor resection followed by adjuvant Bacillus Calmette-Guérin (BCG) bladder instillations. However, only 50% of patients benefit from this therapy. In case of progression to advanced disease, patients must undergo a radical cystectomy with significant morbidity and have a poor clinical outcome. Identifying tumors least likely to respond to BCG can translate into alternative treatments, such as early radical cystectomy or novel targeted or immunotherapies. Here we present molecular profiling of 132 BCG-naive, HR-NMIBC patients, and 44 post-BCG recurrences (34 matched), which uncovered three distinct BCG Response Subtypes (BRS1-3). Patients with BRS3 tumors have reduced recurrence and progression-free survival compared to BRS1-2. BRS3 tumors expressed high EMT-basal markers and had an immunosuppresive profile, which was confirmed with spatial proteomics. Tumors which recurred post-BCG were enriched for BRS3. BRS stratification was validated in a second cohort of 151 BCG-naive HR-NMIBC patients and the molecular subtypes outperformed guideline recommended risk stratification based on clinicopathological variables. For clinical application, we validated that a commercially approved assay was able to accurately predict BRS3 tumors (AUROC 0.86). Our findings provide a potential clinical tool for improved identification of HR-NMIBC patients at the highest risk of progression, which can be used to select patients for early radical cystectomy or novel subtype-directed therapies.

  Study EGAS00001006879

• H3Africa - Stroke Investigative Research and Education Networks

  Although stroke is a discrete phenotype, it is the clinical culmination of several complex and interacting biological processes, precipitated by various genetic and environmental factors, thereby making ready analyses of its underlying mechanisms a challenge. However, the formulation of successful tailor-made prevention strategies for limiting the immense personal and societal burden of stroke in Sub-Saharan Africa (SSA) is contingent on a better understanding of the predisposing risk factors for stroke in the region, which may not necessarily match the nature of stroke risk factors described among other people on other continents. Indeed, the underlying risk factors for stroke in SSA have hitherto not been comprehensively characterized, and currently there is a paucity of expertise, experience, and infrastructure in genomics and environmental research in SSA, to initiate and sustain such investigation. We seek to establish the Stroke Investigative Research and Education Network (SIREN) as a multidisciplinary collaborative research network of investigators in SSA and the United States collectively focused on exploring ways to promote a better understanding of the genetic and environmental risk factors for stroke among people of African ancestry. SIREN will evaluate unique contributors to stroke in SSA and compare these findings with those in a cohort of African American stroke survivors, all while simultaneously building sustainable capacities in phenomics, biobanking, genomics, biostatistics and bioinformatics for future high-level investigation of stroke and other vascular disease entities in SSA.

  Study EGAS00001007331

• Genome-wide association study identifies common variants associated with breast cancer in South African Black women

  Genome-wide association studies (GWAS) have characterized the contribution of common variants to the development of breast cancer (BC) in populations of European descent, however no GWAS have been reported in resident African populations. This GWAS included 2,485 resident African BC cases and 1,101 population matched controls. Two risk loci were identified, one on chromosome 15, located between UNC13C and RAB27A (rs7181788, p = 1.01x10-08) and a second in an intron of USP22 on chromosome 17 (rs899342, p = 4.62x10-08). RAB27A is a member of the RAS oncogene family associated with invasive growth and metastasis, and UNC13C is involved in tumor progression. The rs899342 SNP is an eQTL for expression of USP22 in a wide range of tissues. Several genome-wide significant signals were detected in a hormone receptor subtype analysis. Comparison of a meta-analysis of other African ancestry BC GWAS datasets predominantly of West African ancestry with our South African data showed limited evidence of shared risk loci, suggesting potential heterogeneity in genetic risk factors for breast cancer in different African populations. A European ancestry derived polygenic risk model for BC was applied to our data, explaining only 0.79% of variance, with an area under the curve of 0.56. This study identified genetic risk factors for BC in a resident African population, but larger studies in this and other African populations are needed to further define the contribution of common variants to BC risk.

  Study EGAS00001008032

• 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

• Whole Exome Sequencing data of six chRCC tumors corresponding to three patients

  Mammalian target of rapamycin (mTOR) is a central regulator of mammalian metabolism and physiology. Aberrant hyperactivation of the mTOR pathway promotes tumor growth and metastasis, and can also promote tumor resistance to chemotherapy and cancer drugs; this makes mTOR an attractive cancer therapeutic target. mTOR inhibitors have been approved to treat cancer; however, the mechanisms underlying drug sensitivity remain poorly understood. Here, whole exome sequencing of three chromophobe renal cell carcinoma (chRCC) patients with exceptional mTOR inhibitor sensitivity revealed that all three patients shared somatic mutations in the deubiquitinase USP9X. The clonal characteristics of the mutations, which were amassed by studying multiple patients’ primary and metastatic samples from various years, together with the low USP9X mutation rate in unselected chRCC series, reinforced a causal link between USP9X and mTOR inhibitor sensitivity. Rapamycin treatment of USP9X-depleted HeLa and renal cancer 786-O cells, along with the pharmacological inhibition of USP9X, confirmed that this protein plays a role in patients’ sensitivity to mTOR inhibition. USP9X was not found to exert a direct effect on mTORC1, but subsequent ubiquitylome analyses identified p62 as a direct USP9X target. Increased p62 ubiquitination and the augmented rapamycin effect upon bortezomib treatment, together with the results of p62 and LC3 immunofluorescence assays, suggested that dysregulated autophagy in USP9X-depleted cells can have a synergistic effect with mTOR inhibitors. In summary, we show that USP9X constitutes a potential novel marker of sensitivity to mTOR inhibitors in chRCC patients, and represents a clinical strategy for increasing the sensitivity to these drugs.

  Study EGAS00001007104

• Cohesin Mutations in AML

  Cohesin shapes the nuclear chromatin architecture, including enhancer-promoter interactions, and its components, of which especially STAG2 and RAD21, are frequently mutated in myeloid malignancies. To elucidate mechanisms of leukemogenesis associated with cohesin mutations in humans, we comprehensively characterized genetic, epigenetic, transcriptional, and chromatin conformational changes in acute myeloid leukemia (AML). To corroborate our findings, we performed complementary siRNA-mediated depletion of STAG2, its paralogue STAG1 or RAD21 in cord blood-derived CD34+ primary human hematopoietic stem and progenitor cells (HSPCs). We show that STAG2 mutations consistently lead to the loss of STAG2 protein and are associated with a specific set of co-occurring mutations, while STAG1 was never mutated in AML. Loss of STAG2 was frequently compensated by STAG1. Still, specific loci displayed altered cohesin occupancy, gene expression and corresponding changes in local chromatin activation as measured by H3K27ac enrichment and chromatin accessibility. High-throughput chromosome conformation capture (in-situ Hi-C) revealed significantly altered chromatin looping in cohesin-mutated AMLs, including weakened enhancer-promoter contacts with reduced, cohesin-dependent promoter activity. In HSPCs, we detected transcriptomic and epigenetic effects overlapping STAG2-mutant AML-specific changes following STAG2 knockdown (KD), that were not invoked by the depletion of STAG1. We also found that STAG2 loss in cultured HSPCs impaired the differentiation capacity, especially erythroid colony formation which maintained HSPC-like gene expression. This work establishes STAG2 as a key regulator of cohesin-associated chromatin architecture, gene expression and differentiation in the human hematopoietic system and identifies candidate target genes that may be implicated in leukemogenesis.

  Study EGAS00001007405

• Evaluation of EBUS-TBNA Aspirates from Advanced NSCLC for Comprehensive Sequencing Platforms Including Whole Genome Sequencing

  Endobronchial ultrasound-guided transbronchial needle aspiration (EBUS-TBNA) is often the only source of tumor tissue from patients with advanced, inoperable lung cancer. EBUS-TBNA aspirates are used for the diagnosis, staging, and genomic testing to inform therapy options. Here we extracted DNA and RNA from 220 EBUS-TBNA aspirates to evaluate their suitability for whole genome (WGS), whole exome (WES), and comprehensive panel sequencing. For a subset of 40 cases, the same nucleic acid extraction was sequenced using WGS, WES, and the TruSight Oncology 500 assay. Genomic features were compared between sequencing platforms and compared with those reported by clinical testing. A total of 204 aspirates (92.7%) had sufficient DNA (100 ng) for comprehensive panel sequencing, and 109 aspirates (49.5%) had sufficient material for WGS. Comprehensive sequencing platforms detected all seven clinically reported tier 1 actionable mutations, an additional three (7%) tier 1 mutations, six (15%) tier 2–3 mutations, and biomarkers of potential immunotherapy benefit (tumor mutation burden and microsatellite instability). As expected, WGS was more suited for the detection and discovery of emerging novel biomarkers of treatment response. WGS could be performed in half of all EBUS-TBNA aspirates, which points to the enormous potential of EBUS-TBNA as source material for large, well-curated discovery-based studies for novel and more effective predictors of treatment response. Comprehensive panel sequencing is possible in the vast majority of fresh EBUS-TBNA aspirates and enhances the detection of actionable mutations over current clinical testing.

  Study EGAS00001007708

• RNA isoform repertoire of neuropsychiatric risk genes in human brain

  Neuropsychiatric disorders are highly complex conditions and the risk of developing a disorder has been tied to hundreds of genomic variants that alter the expression and/or products (isoforms) made by risk genes. However, how these genes contribute to disease risk and onset through altered expression and RNA splicing is not well understood. Here we show our current understanding of gene isoforms is far from complete and reveal the precise splicing profiles of neuropsychiatric disorder risk genes. Combining our new bioinformatic pipeline IsoLamp with nanopore long-read amplicon sequencing, we deeply profiled the RNA isoform repertoire of 31 high-confidence neuropsychiatric disorder risk genes in human brain. We show most risk genes are more complex than previously reported, identifying 443 novel isoforms and 28 novel exons, including isoforms which alter protein domains, and genes such as ATG13 and GATAD2A where most expression was from previously undiscovered isoforms. The greatest isoform diversity was present in the schizophrenia risk gene ITIH4. Mass spectrometry of brain protein isolates confirmed translation of a novel exon skipping event in ITIH4, suggesting a new regulatory mechanism for this gene in brain. Our results emphasize the widespread presence of previously undetected RNA and protein isoforms in brain and provide an effective approach to address this knowledge gap. Uncovering the isoform repertoire of neuropsychiatric risk genes will underpin future analyses of the functional impact these isoforms have on neuropsychiatric disorders, enabling the translation of genomic findings into a pathophysiological understanding of disease.

  Study EGAS00001007744

• The University of Hong Kong Intestinal Metaplasia Organoids Study

  Background: Gastric intestinal metaplasia (IM) is a pre-cancerous stage spanning a morphological spectrum that is poorly represented by human cell line models. Objective: We aim to establish and characterize human IM cell models to better understand IM progression along the cancer spectrum. Design: A large human gastric IM organoid (IMO) cohort (n=28), their clonal derivatives and normal gastric organoids (n=42) for comparison were established. Comprehensive multi-omics profiling and functional characterization were performed. Results: Single-cell transcriptomes revealed IMO cells spanning a spectrum from hybrid gastric/intestinal to advanced intestinal differentiation. Their lineage trajectories connected different cycling and quiescent stem and progenitors, highlighting differences in gastric to IM transition and the potential origin of IM from STMN1 cycling isthmus stem cells. Hybrid IMOs showed impaired differentiation potential, high lineage plasticity beyond gastric or intestinal fates, and reactivation of a fetal gene program. Cell populations in gastric IM and cancer tissues were highly similar to those derived from IMOs and exhibited a fetal signature. Genomically, IMOs showed elevated mutation burden, frequent chromosome 20 gain, and epigenetic de-regulation of many intestinal and gastric genes. Functionally, IMOs were FGF10 independent and showed downregulated FGFR2. Several IMOs exhibited a cell-matrix adhesion independent (CMi) subpopulation that displayed chromosome 20 gain but lacked key cancer driver mutations, potentially representing the earliest neoplastic precursor of IM-induced gastric cancer. Conclusions: Overall, our IM organoid biobank captured the heterogeneous nature of IM, revealing mechanistic insights on IM pathogenesis and progression, offering an ideal platform for studying early gastric neoplastic transformation and chemoprevention.

  Study EGAS00001007899

• Epigenome-wide methylation analysis of colorectal carcinoma, adenoma and normal tissue reveals novel biomarkers addressing unmet clinical needs

  Studies about CRC biomarker discovery have focused on methylation patterns in normal and colorectal tumor tissue, leading to a lack of knowledge of methylation in adenomas. Therefore, we performed the first epigenome-wide study to compare the epigenomes of all three tissue types combined, and to identify discriminatory biomarkers. Public Illumina MethylEPIC® data was collected from 552 normal, 118 carcinoma and 116 adenoma samples. Pairwise analyses were performed for discovery of differentially methylated probes (DMPs). We identified 693 813, 620 643 and 558 897 DMPs when comparing normal vs carcinoma, normal vs adenoma and adenoma vs carcinoma respectively. To double evidence (DE) these DMPs, analysis of Illumina 450K data from public datasets was performed. Through this analysis, 13 DE DMPs with │Δβ│≥0.3 were identified for adenoma vs carcinoma. A binary logistic regression model was generated to evaluate the discriminatory potential of these DE DMPS. The model was validated in an in-house experimental methylation dataset of 13 adenoma and 9 carcinoma samples. Our classifier reached a cross-validated AUC of 0.996 and 0.855 in the discovery and validation datasets respectively. A sensitivity and specificity of 96% and 95% were reached, with an overall accuracy of 96%. Our analysis show that methylation biomarkers have the potential to discriminate between normal, precursor and carcinoma tissues of the colorectum. More importantly, we highlight the power of the methylome, showing that DMPs can be used as biomarkers in the clinic to identify colorectal adenoma and carcinoma. 

  Study EGAS00001007017

• pancreatic ductal adenocarcinoma exomes in study: Pro-immunogenic Impact of MEK inhibition combined with an anti-CD40 immunostimulatory antibody

  Cancer types with lower mutational load and a non-permissive tumor microenvironment are intrinsically resistant to immune checkpoint blockade. While the combination of cytostatic drugs and immunostimulatory antibodies constitutes an attractive concept for overcoming this refractoriness, suppression of immune cell function by cytostatic drugs may limit therapeutic efficacy. Here we show that targeted inhibition of mitogen-activated protein kinase (MAPK) kinase (MEK) does not impair dendritic cell-mediated T-cell priming and activation. Accordingly, combining MEK inhibitors (MEKi) with agonist antibodies (Abs) targeting the immunostimulatory CD40 receptor resulted in potent synergistic anti-tumor efficacy. Detailed analysis of the mechanism of action of MEKi GDC-0623 by means of flow cytometric analysis of the tumor immune infiltrate and whole tumor transcriptomics showed that, in addition to its cytostatic impact on tumor cells, this drug exerts multiple pro-immunogenic effects, including the suppression of M2-type macrophages, myeloid derived suppressor cells and CD4+ T-regulatory cells. In addition, MEKi was found to induce tumor-cell intrinsic interferon signaling, which contributed to antigen presentation by tumor cells. Finally, the tumoridical impact of MEKi involves the activation of multiple pro-inflammatory pathways involved in immune cell effector function in the tumor microenvironment. Our data therefore indicate that the combination of MEK inhibition with agonist anti-CD40 Ab is a promising therapeutic concept, especially for the treatment of mutant Kras-driven tumors such as pancreatic ductal adenocarcinoma. The pancreatic ductal adenocarcinoma exomes were generated in the context of the Heidelberg Center for Personalized Oncology project HIPO-K28E.

  Study EGAS00001004196

• Targeted and shallow whole genome sequencing identifies therapeutic opportunities in p53abn endometrial cancers

  Purpose: Shallow whole genome sequencing (sWGS) can detect copy number (CN) aberrations. In high-grade serous ovarian (HGSOC) sWGS identified CN signatures such as homologous recombination deficiency (HRD) to direct therapy. We applied sWGS with targeted sequencing to p53abn endometrial cancers (ECs) to identify additional prognostic stratification and therapeutic opportunities. Experimental Design: sWGS and targeted panel sequencing was performed on formalin-fixed paraffin-embedded p53abn ECs. CN alterations, mutational data and CN signatures were derived, and associations to clinicopathologic and outcomes data were assessed. Results: In 187 p53abn ECs, 5 distinct CN signatures were identified. Signature 5 was associated with BRCA1/2 CN loss with features similar to HGSOC HRD signature. 22% potential HRD cases were identified, 35 patients with signature 5, and 8 patients with BRCA1/2 somatic mutations. Signatures 3 and 4 were associated with whole genome duplication, and CCNE1, ERBB2 and MYC amplifications, with mutations in PIK3CA enriched in signature 3. We observed improved overall survival (OS) for patients with signature 2 and worse OS for signatures 1 and 3. 28% of patients had CCNE1 amplification and this subset was enriched with carcinosarcoma histotype. 34% of patients, across all histotypes, had ERBB2 amplification and/or HER2 overexpression on immunohistochemistry, which was associated with worse outcomes. Mutations in PPP2R1A (29%) and FBXW7 (16%) were among the top 5 most common mutations. Conclusions: sWGS and targeted sequencing identified therapeutic opportunities in 75% of p53abn EC patients. Further research is needed to determine the efficacy of treatments targeting these identified pathways within p53abn ECs.

  Study EGAS00001007661

• Organoid cultures of early-onset colorectal cancers reveal distinct and rare genetic profiles

  Sporadic early-onset colorectal cancer (EOCRC) has bad prognosis, yet is poorly represented by cell line models. We examine the key mutational and transcriptomic alterations in an organoid biobank enriched in EOCRCs. We established paired cancer (n=32) and normal organoids (n=18) from 20 patients enriched in microsatellite-stable EOCRC. Exome and transcriptome analysis was performed. We observed a striking diversity of molecular phenotypes, including PTPRK-RSPO3 fusions. Transcriptionally, RSPO fusion organoids resembled normal colon organoids and were distinct from APC mutant organoids, with high BMP2 and low PTK7 expression. Single cell transcriptome analysis confirmed the similarity between RSPO fusion organoids and normal organoids, with a propensity for maturation on Wnt withdrawal, whereas the APC mutant organoids were locked in progenitor stages. CRISPR/Cas9 engineered mutation of APC in normal human colon organoids led to upregulation of PTK7 protein and suppression of BMP2, but less so with an engineered RNF43 mutation. The frequent co-occurrence of RSPO fusions with SMAD4 or BMPR1A mutation was confirmed in TCGA database searches. RNF43 mutation was found in organoid from a leukaemia survivor with a novel mutational signature; and organoids with POLE proofreading mutation displayed ultramutation. The cancer organoid genomes were stable over long culture periods, while normal human colon organoids tended to be subject to clonal dominance over time. These organoid models enriched in EOCRCs with linked genomic data fill a gap in existing CRC models and reveal distinct genetic profiles and novel pathway cooperativity.

  Study EGAS00001004063

• Pediatric Papillary Thyroid Carcinoma Whole Exome Sequencing

  BACKGROUND: Pediatric papillary thyroid carcinoma (PPTC) is distinct from adult-onset disease. Despite higher rates of regional and distal metastases and more frequent recurrence, there is, paradoxically, a highly favourable prognosis. Historically, only 40-50% of childhood PTC are driven by genomic variants common to adult PTC, while the oncogenic drivers in the remainder are unknown. METHODS: Fifty-two PPTCs were genotyped by candidate gene testing, followed, in a subset (n=32), by whole exome (WES) and transcriptome sequencing (RNASeq) to identify oncogenic variants in driver-unknown tumors and to explore gene expression patterns. FINDINGS: Thirty-one tumors (60%) carried variants identified by candidate gene testing, specifically: BRAFV600E(22%), NRASQ61R(2%), RET-CCDC6(19%), RET-NCOA4(15%) and PAX8-PPARG. Among the 21 remaining tumors with negative first-line genotyping, oncogenic variants were identified in 19 (90%). These were enriched with oncogenic fusions, with 11 non-recurrent fusion transcripts identified, including two that had not been previously described: STRN-RET and TG-PBF. Most fusions were associated with 3’ receptor tyrosine kinase (RTK) moieties: RET, MET, ALK and NTRK3. Gene expression analysis defined three distinct clusters, which demonstrated differing expression of genes involved in thyroid differentiation and MAPK signaling. INTERPRETATION: Exome and transcriptome-wide sequencing identify genomic drivers of PPTC in 98% of tumors, whereas conventional genotyping leaves many tumors genetically unexplained. PPTC is driven, in the majority of cases, by oncogenic fusion transcripts involving one of a limited number of RTKs. These observations have broad implications for diagnostic testing and for therapeutic options in children with advanced disease.

  Study EGAS00001005187

• Analysis of recent shared ancestry in a familial cohort identifies coding and noncoding autism spectrum disorder variants

  Autism spectrum disorder (ASD) is a collection of neurodevelopmental disorders characterized by deficits in social communication and restricted, repetitive patterns of behavior or interests. ASD is highly heritable, but genetically and phenotypically heterogeneous, reducing the power to identify causative genes. We performed whole genome sequencing (WGS) in an ASD cohort of 68 individuals from 22 families enriched for recent shared ancestry. We identified an average of 3.07 million variants per genome, of which an average of 112,512 were rare. We mapped runs of homozygosity (ROHs) in affected individuals and found an average genomic homozygosity of 9.65%, consistent with expectations for multiple generations of consanguineous unions. We identified potentially pathogenic rare exonic or splice site variants in 12 known (including KMT2C, SCN1A, SPTBN1, SYNE1, ZNF292) and 12 candidate (including CHD5, GRB10, PPP1R13B) ASD genes. Furthermore, we annotated noncoding variants in ROHs with brain-specific regulatory elements and identified putative disease-causing variants within brain-specific promoters and enhancers for 5 known ASD and neurodevelopmental disease genes (ACTG1, AUTS2, CTNND2, CNTNAP4, SPTBN4). We also identified copy number variants in two known ASD and neurodevelopmental disease loci in two affected individuals. In total we identified potentially etiological variants in known ASD or neurodevelopmental disease genes for ~61% (14/23) of affected individuals. We combined WGS with homozygosity mapping and regulatory element annotations to identify candidate ASD variants. Our analyses add to the growing number of ASD genes and variants and emphasize the importance of leveraging recent shared ancestry to map disease variants in complex neurodevelopmental disorders.

  Study EGAS00001006058

• Somatic mutations predict an aggressive phenotype in meningioma but do not drive grade progression

  A subset of meningiomas progress in histopathological grade and drivers of progression are poorly understood. We aimed to identify somatic mutations and copy number alterations (CNAs) associated with grade progression in a unique matched tumour dataset. Utilising a prospective database, we identified 10 patients with meningiomas that had undergone grade progression and for whom matched pre- and post-progression tissue (n=50 samples) was available for targeted next-generation sequencing. Mutations in NF2 were identified in 4/10 patients, 94% were non-skull base tumours. In one patient, three different NF2 mutations were identified in four tumours. NF2 mutated tumours showed large scale CNAs, with highly recurrent losses in 1p, 10, 22q, and frequent CNAs on chromosomes 2, 3 and 4. There was a correlation between grade and CNAs in a subset of patients. In two patients with sphenoid wing Grade III tumours without detected NF2 mutations, a STAT6-NAB2 inversion was identified. Two patients with tumours without detected NF2 mutations showed a combination of loss and high gain on chromosome 17q. Mutations in SETD2, TP53, TERT promoter and NF2 were not uniform across recurrent tumours and did not correspond with grade progression. Meningiomas that progress in grade generally have a mutational profile already detectable in the pre-progressed tumour, suggesting an aggressive phenotype . CNA profiling shows frequent alterations in NF2 mutated tumours compared to non NF2 mutated tumours and the pattern of CNAs may be associated with grade progression in a subset of cases.

  Study EGAS00001006585

• Transcriptional_reprogramming_from_innate_immune_functions_to_a_pro_thrombotic_signature_upon_SARS_CoV_2_sensing_by_monocytes_in_COVID_19

  Alterations in the myeloid immune compartment have been observed in COVID-19, but the specific mechanisms underlying these impairments are not completely understood. Here we examined the functionality of classical CD14+ monocytes as a main myeloid cell component in well-defined cohorts of patients with mild and moderate COVID-19 during the acute phase of infection and compared them to that of healthy individuals. We found that ex vivo isolated CD14+ monocytes from mild and moderate COVID-19 patients display specific expression patterns of costimulatory and inhibitory receptors that clearly distinguish them from healthy monocytes, as well as an altered metabolic profile. In addition, decreased NFB activation in COVID-19 monocytes ex vivo is accompanied by an intact type I IFN antiviral response. Secondary pathogen sensing ex vivo led to a state of functional unresponsiveness characterized by a defect in pro-inflammatory cytokine expression, NFB-driven cytokine responses and defective type I IFN response in moderate COVID-19 monocytes, together with defects in the metabolic reprogramming that innate immune cells usually undergo upon pathogen sensing. Transcriptionally and functionally, COVID-19 monocytes switched their gene expression signature from canonical innate immune functions to a pro-thrombotic phenotype characterized by enrichment of pathways involved in hemostasis, immunothrombosis, platelet aggregation and other accessory pathways to platelet activation, accumulation and clot formation, including extracellular matrix reorganization, integrin cell surface interactions and signaling by PDGF. These results provide a potential mechanism by which innate immune dysfunction in COVID-19 may contribute to disease pathology

  Study EGAS00001006788

• A Functional Network of Gastric-Cancer-Associated Splicing Events Controlled by Dysregulated Splicing Factors

  Except for a few genes known to have oncogenic spliceforms in gastric cancer, the full scope of aberrant splicing in gastric oncogenesis remains unclear. In this study, we elucidated the alternative splicing events associated with gastric carcinogenesis and identified upstream regulators governing these events. We performed RNA-seq on 19 matched tumor/normal pairs and a panel of gastric cancer cell lines, and identified 361 tumor associated (TA) alternative splicing events (ASEs), which included known oncogenic ASEs in genes such as INSR, FGFR2, CD44 and KRAS. We further identified 8 splicing factors dysregulated in gastric tumors, the expression of which are correlated to the splice ratio of the TA ASEs. We thereby constructed a dysregulated splicing network in gastric cancer, consisting of alternative splicing changes correlated to their potential regulators. This network featured three potential regulatory modules centered around the splicing factors ESRP2, MBNL1 and PTPB1. Knockdown of each splicing factor led to the expected changes in splicing of approximately half of the ASEs assayed. Thus, dysregulation of the splicing factors indeed was the cause of many of the TA ASEs. Pathway enrichment analysis showed that genes affected by tumor associated splicing were pre-dominantly involved in cytoskeletal organization. Concordant with this, knockdown of the splicing factors regulating the TA ASEs led to the expected changes in trans-well migration activity. Thus, we have established for the first time a comprehensive splicing network that likely drives gastric tumor progression, based on which novel therapeutic targets can be identified for further evaluation.

  Study EGAS00001002256

• Spatiotemporal evolution and inter-patient heterogeneity in primary and recurrent/metastatic head and neck squamous cell carcinoma

  Head and neck cancer squamous cell carcinomas (HNSCCs) are heterogeneous in terms of origin and aetiology. In addition, there is uncertainty about the spatiotemporal genetic evolution from initial diagnosis to recurrence/metastasis (R/M) after primary treatments and further disease progression following systemic treatment. Changes in the genetic profile have implications on the selection of appropriate treatments for patients, especially in the era of targeted therapies and immunotherapies. Here, we analysed a cohort of 11 HNSCC patients with metachronous R/M, among whom nine had paired primary and R/M samples suitable for next-generation sequencing (multiple R/M samples collected at different time points were available for some of the patients). At the genomic level, the R/M samples shared a fraction of the somatic single nucleotide variants (SNVs) with the index primary tumours, but they also acquired many additional mutations, while losing only a few others. A similar behaviour was also observed when examining the changes of mutational signatures between primary and R/M samples. Overall, R/M appeared thus more genetically diverse than the respective primary tumours. The transcriptomic analysis showed that R/M samples had lower immune cell infiltration and several genes related to immune response were significantly downregulated compared to the primary samples. Our results underline the importance of analysing multiple samples per patient to obtain a more complete picture of the patient’s tumour and advocates a re-biopsy in the event of R/M and treatment failure, in order to select the most appropriate therapeutic strategy.

  Study EGAS00001007464

• Benchmarking for alignment and variant calling

  Background: Lung cancer is a heterogeneous disease and the primary cause of cancer-related mortality worldwide. Somatic mutations, including large structural variants, are important biomarkers in lung cancer for selecting targeted therapy. Genomic studies in lung cancer have been conducted using short-read sequencing. Emerging long-read sequencing technologies are a promising alternative to study somatic structural variants, however there is no current consensus on how to process data and call somatic events. In this study, we preformed whole genome sequencing of lung cancer and matched non-tumour samples using long and short read sequencing to comprehensively benchmark three sequence aligners and six structural variant callers comprised of four generic callers (Sniffles2, cuteSV, SVIM and DELLY in generic mode) and two somatic callers (nanomonsv and DELLY in somatic modes). Results: Different combinations of aligners and variant callers influenced somatic structural variant detection. The choice of caller had a significant influence on somatic structural variant detection in terms of variant type, size, sensitivity, and accuracy. The performance of each variant caller was assessed by comparing to somatic structural variants identified by short read sequencing. When compared to somatic structural variants detected with short read sequencing, more events were detected with long read sequencing. The mean recall of somatic variant events identified by short read sequencing was higher for the somatic callers (64%) than generic callers (52%), with nanomonsv (69%) showing the best performance. Conclusion: Long read sequencing can identify somatic structural variants. The longer reads have the potential to improve our understanding of cancer development and inform personalized cancer treatment.

  Study EGAS00001007819

• High-Resolution Spatial Transcriptomics Uncover Epidermal-Dermal Divergences in Merkel Cell Carcinoma: Spatial Context Reshapes the Gene Expression Landscape

  Merkel cell carcinoma (MCC) is an aggressive malignancy with neuroendocrine differentiation marked by high plasticity, often manifesting as rapid therapy resistance. Although the cell-of-origin is presumed to be epithelial, epidermal localization of MCC is rarely observed, largely because in situ MCC is typically an incidental finding. Nevertheless, a subset of MCC tumors exhibits epidermotropism, wherein tumor cells are present in the epidermis. The behavior of cancer cells is profoundly influenced by the tumor microenvironment and interactions with neighboring cells. Notably, the normal counterparts of the cancer’s cell-of-origin have been shown to attenuate tumor aggressiveness. Thus, epidermotropic MCC presents a unique opportunity to explore the potential role of surrounding keratinocytes in modulating tumor cell behavior. While the epidermotropic tumor nests share histological resemblance with their dermal counterparts, their transcriptomic profiles remain unexplored. Here, we employed high-definition spatial and single cell transcriptomics to dissect the gene expression profiles of epidermotropic MCC cells, comparing them to MCC cells in the tumor core and those adjacent to blood vessels. Notably, epidermotropic MCC cells exhibit a transcriptomic signature reminiscent of cutaneous squamous cell carcinoma, characterized by upregulation of genes encoding keratins, S100A proteins, as well as calmodulin-like proteins 3 and 5. Mechanistically, this keratinocytic differentiation is associated with enhanced p63 activity, leading to the upregulation of PERP. Collectively, our study demonstrates that MCC cells can adopt a keratinocytic differentiation program in response to microenvironmental cues, underscoring the remarkable phenotypic plasticity of this malignancy and the importance of the tumor microenvironment for tumor characteristics.

  Study EGAS00001008157

• OMKar

  The whole genome karyotype refers to the sequence of large chromosomal segments that make up an individuals genotype. karyotype analysis, which includes descriptions of aneuploidies and other rearrangements is crucial for understanding genetic risk factors, for diagnosis, treatment decisions, and genetic counseling linked to constitutional disorders. The current karyotyping standard is based on microscopic examination of chromosomes, a complex process that requires high expertise and offers Mb scale resolution. Optical Genome Mapping (OGM) technology can identify large DNA lesions in a cost-effective manner. In this paper, we developed OMKar, a method that uses OGM data to create a virtual karyotype. OMKar processes Structural (SV) and Copy Number (CN) Variants as inputs and encodes them into a compact breakpoint graph. It recomputes copy numbers using Integer Linear Programming to maintain CN balance and then identifies constrained Eulerian paths representing entire donor chromosomes. In tests using 38 whole genome simulations of constitutional disorders, OMKar reconstructed the karyotype with 88% precision and 95% recall on SV concordance and 95% Jaccard score on CN concordance. We applied OMKar to 50 prenatal, 41 postnatal, and 63 parental samples from ten different sites. OMKar reconstructed the correct karyotype in 144 out of 154 samples, covering 25 of 25 aneuploidies, 32 of 32 balanced translocations, and 72 of 82 unbalanced variations. Detected constitutional disorders included Cri-du-chat, Wolf-Hirschhorn, Prader-Willi deletions, Down, and Turner syndromes. Importantly, it identified a plausible genetic mechanism for five cases of constitutional disorder that were not detected by other technologies. Together, these results demonstrate the robustness of OMKar for OGM-based karyotyping.

  Study EGAS00001008245

• Poly(A) RNA sequencing of hepatocellular carcinoma tumors and their matched noncancerous liver tissues

  Background &amp; Aims: The differentiation of distinct multifocal hepatocellular carcinoma (HCC): multicentric disease vs. intrahepatic metastases, in which the management and prognosis varies substantively, remains problematic. We aim to stratify multifocal HCC and identify novel diagnostic and prognostic biomarkers by performing whole genome and transcriptome sequencing, as part of a multi-omics strategy Methods: A complete collection of tumour and somatic specimens (intrahepatic HCC lesions, matched non-cancerous liver tissue and blood) were obtained from representative patientswith multifocal HCC exhibiting two distinct postsurgical courses. Whole-genome and transcriptome sequencing with genotyping were performed for each tissue specimen to con-trast genomic alterations, including hepatitis B virus integrations, somatic mutations, copy number variations, and structural variations. We then constructed a phylogenetic tree to visualise individual tumour evolution and performed functional enrichment analyses on select differentially expressed genes to elucidate biological processes involved in multifocal HCC development. Multi-omics data were integrated with detailed clinicopathological information to identify HCC biomarkers, which were further validated using a large cohort of HCC patients (n = 174).Results: The multi-omics profiling and tumour biomarkers could successfully distinguish the two multifocal HCC types, while accurately predicting clonality and aggressiveness. The dual specificity protein kinase TTK, which is a key mitotic checkpoint regulator with links to p53 signaling, was further shown to be a promising overall prognostic marker for HCC in the large patient cohort.Conclusions: Comprehensive multi-omics characterisation of multifocal tumour evolution may improve clinical decision-making, facilitate personalised medicine, and expedite identification of novel biomarkers and therapeutic targets in HCC.

  Study EGAS00001002337

• Genomic consequences of aberrant DNA repair mechanisms stratify ovarian cancer histotypes

  Aberrant DNA repair processes are a hallmark of human tumours. How these deficiencies variously impact the genomes of ovarian cancer between and within histotypes remains unknown. We studied the whole genome pointmutation and structural variation patterns of 133 tumours (59 high grade serous (HGSC), 35 clear cell (CCOC), 29 endometrioid (ENOC), and 10 adult granulosa cell tumours (GCT)) as a substrate for class discovery in ovarian cancer. Ab-initio clustering of integrated point mutation and structural variation signatures revealed seven novel subgroups, comprising between and within histotype stratification. Prevalence of fold-back inversions (FBI) co-localised with high level amplification events divided HGSC into two prognostically significant subgroups. This finding was recapitulated in two additional independent cohorts (total n=576 cases), and transcended gene based mutation status and gene expression as prognostically relevant features of HGSC. CCOC were divided into one group harbouring a signature reflective of active genome editing via APOBEC enzymes (26%) and a complementary group with an age-relatedmutational signature. ENOC were divided intomicrosatellite instable (MSI) cases (28%) with a distinctmismatch repairmutation signature and an outlying mutation rate, with the remainder of cases distributed amongst the remaining six groups. We suggest the FBI, MSI, and APOBEC groups represent biologic strata that could direct treatment strategies. Taken together, our work establishes the efficacy of the somatic genome as a biomarker to stratify ovarian cancers, simultaneously identifying patients that may benefit from emerging therapeutics and distinct subgroups within classic ovarian cancer histotypes.

  Study EGAS00001002390

• A_systems_biology_approach_to_understand_immunity_and_pathogenesis_of_malaria_in_children_exposed_to_endemic_Plasmodium_falciparum_transmission

  From long-term studies of children in Kenya, we know that some children have numerous malaria episodes through to older age, while other children become immune more rapidly. Our hypothesis is that these children with frequent malaria episodes are caught in a vicious circle, whereby malaria episodes lead to impaired immunity to malaria, which in turn leads to further episodes of malaria. We will therefore investigate immune responses in children with different life histories of malaria episodes. We will also include a group of children living in similar environmental conditions but who are completely unexposed to malaria. Aims: Using RNAseq we will characterise the transcriptome of a group of children with frequent and febrile malaria or infrequent/asymptomatic. These data will be combined with cytokine profiling data (NIMR) and used to build predictive models (Exeter). The predictions will be validated against a similar second group. A single snap shot of immune responses may represent only the endpoint of many immune processes and may not reflect those that are causally related to differences in malaria outcome. Therefore, we will collect and store samples from a larger group of children who will then be followed up for over 3-4 years. This extended follow up will allow us to identify the same groups as described above, and then examine samples taken before the development of those outcomes. 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 EGAS00001002978

• Study of the rare and low frequency variants in the Saguenay-Lac-Saint-Jean population

  The Saguenay–Lac-Saint-Jean (SLSJ) region is located in northeastern Quebec and is known for its unique demographic history and founder effect. Since founder populations are enriched with population-specific variants, we characterized the variants distribution in SLSJ and compared it with four European populations (Finnish, Sweden, United Kingdom and France), of which the Finnish population is another founder population. Targeted sequencing of the coding and non-coding immune regulatory regions of the SLSJ asthma familial cohort and the four European populations were performed. Rare and low-frequency coding and noncoding regulatory variants identified in the SLSJ population were then investigated for variant- and gene-level associations with asthma and allergy related-traits (eosinophil percentage, immunoglobulin (Ig)E levels and lung function). Our data showed that (1) rare or deleterious variants were not enriched in the two founder populations as compared to the three non-founder European populations; (2) a larger proportion of founder population-specific variants occurred with higher frequencies; and (3) low-frequency variants appeared to be more deleterious. Furthermore, a rare variant, rs1386931, located in the 3’UTR of CXCR6 and intron of FYCO1 was found to be associated with eosinophil percentage. Gene-based analyses identified NRP2, MRPL44 and SERPINE2 to be associated with various asthma and allergy related traits. Our study demonstrated the usefulness of using a founder population to identify new genes associated with asthma and allergy-related traits; thus better understand the genes and pathways implicated in pathophysiology.

  Study EGAS00001003103

• Breast cancer DNA repair

  Tumors of germline BRCA1/2 mutated carriers show homologous recombination (HR) deficiency (HRD), resulting in impaired DNA double strand break (DSB) repair and high sensitivity to Poly-(ADP-Ribose)-Polymerase (PARP) inhibitors. Although this therapy is expected to be effective beyond germline BRCA1/2 mutated carriers, a robust validated test to detect HRD tumors is lacking. In the present study we therefore evaluated a functional HR assay exploiting the formation of RAD51 foci in proliferating cells after ex vivo irradiation of fresh breast cancers (BrC) tissue: the RECAP test.Methods Fresh samples of 170 primary BrC were analyzed using the RECAP test. The molecular explanation for the HRD phenotype was investigated by exploring BRCA deficiencies, mutational signatures, tumor infiltrating lymphocytes (TILs) and microsatellite instability (MSI).Results RECAP was completed successfully in 148 out of 170 samples (87%). 24 tumors showed HRD (16%), while 6 tumors were HR intermediate (HRi) (4%). HRD was explained by BRCA deficiencies (mutations, promoter hypermethylation, deletions) in 16 cases. Several non-BRCA deficient HRD tumors showed BRCAness mutational signatures, suggesting that they are also bona fide HRD cases. HRD tumors showed an increased incidence of high TIL counts (p=0.023) compared to HR proficient (HRP) tumors and MSI was more frequently observed in the HRD group (2/20, 10%) than expected in BrC (1%) (p=0.017).Conclusion The RECAP test is a robust assay detecting both BRCA1/2 deficient and BRCA1/2 proficient HRD tumors, suggesting that this test identifies approximately 50% more patients that may benefit from PARPi treatment than BRCA gene testing only.

  Study EGAS00001002792

• A Single Complex Agpat2 Allele In A Patient With Partial Lipodystrophy

  Genetic lipodystrophies are a group of rare syndromes associated with major metabolic complications - including severe insulin resistance, type 2 diabetes mellitus, and hypertriglyceridemia - which are classified according to the distribution of adipose tissue. Lipodystrophies can be present at birth or develop during life and can range from local to partial and general. With at least 18 different genes implicated so far, definite diagnosis can be challenging due to clinical and genetic heterogeneity. In an adult female patient with clinical and metabolic features of partial lipodystrophy we identified via whole genome sequencing a single complex AGPAT2 allele [V67M;V167A], functionally equivalent to heterozygosity. AGPAT2 encodes for an acyltransferase implicated in the biosynthesis of triacylglycerol and glycerophospholipids. So far homozygous and compound heterozygous mutations in AGPAT2 have only been associated with generalized lipodystrophy. A SNP risk score analysis indicated that the index patient is not predisposed to lipodystrophy based on her genetic background. The partial phenotype in our patient is therefore more likely associated to the genetic variants in AGPAT2. To test whether the resulting double-mutant AGPAT2 protein is functional we analysed its in vitro enzymatic activity via mass spectrometry. The resulting AGPAT2 double mutant is enzymatically inactive. Our data support the view that the current classification of lipodystrophies as strictly local, partial or generalized may have to be re-evaluated and viewed more as a continuum, both in terms of clinical presentation and underlying genetic causes. Better molecular understanding of lipodystrophies may lead to new therapies to treat adipose tissue dysfunction in common and rare diseases.

  Study EGAS00001003177

• MYCN Amplification and ATRX Mutations are Incompatible in Neuroblastoma

  Age and stage of disease at diagnosis is associated with outcome in neuroblastoma. A previous study identified ATRX point mutations and in-frame deletions in older patients with stage 4 disease. To validate the frequency of ATRX mutations and assess the relation of mutations and known prognostic variables, we collaborated with the Children’s Oncology Group and sequenced the entire ATRX genomic locus in 475 neuroblastoma tumor samples to validate the frequency of ATRX mutations and assess the relation of mutations and known prognostic variables. Older age at diagnosis, stage 4 disease, and unfavorable histology were significantly associated with mutations in the ATRX gene; MYCN amplification was mutually exclusive from ATRX mutations. To directly test if ATRX mutations and MYCN amplification are incompatible in neuroblastoma, we mutated the ATRX gene in MYCN amplified neuroblastoma cell lines and ectopically expressed MYCN in ATRX mutant neuroblastoma cell lines. Both approaches showed that ATRX and MYCN amplification are incompatible in neuroblastoma in culture and in vivo. To gain a better understanding of the underlying molecular and cellular mechanisms of this incompatibility, we performed a series of experiments including telomere analysis, epigenetic profiling, electron microscopic analysis and metabolic profiling. We discovered that induction of MYCN expression in ATRX mutant neuroblastomas cells leads to metabolic reprogramming, mitochondrial dysfunction, increased reactive oxygen species and replicative stress. We propose that this leads to synthetic lethality in ATRX mutant neuroblastomas because of the underlying replicative stress in those cells as a result of dysfunction of this essential histone chaperone.

  Study EGAS00001003257

• Circumventing intratumoral heterogeneity to identify potential therapeutic targets in hepatocellular carcinoma. Details please contact lifuqiang@genomics and zhaoxin@genomics.cn.

  Background and aims: Intratumoral heterogeneity (ITH) challenges identifying mutations with target therapy potential whereas circulating cell-free DNAs (cfDNAs) could reflect nearly the entire mutation spectrum in given tumors. We investigated how to minimize the limit of ITH for profiling hepatocellular carcinoma (HCC).Methods: Thirty-two multi-regional HCC samples from five patients were subjected to whole exome sequencing (WES) and targeted deep sequencing (TDS). ITH extent was measured by the average percentage of non-ubiquitous mutations (present in parts of tumor regions). Matched cfDNAs were also analyzed by WES and TDS. Profiling efficiencies of single tumor specimen and cfDNA were compared and the one better depicted mutational landscape was selected to screen therapeutic targets.Results: We found variable extents of ITH in HCCs and observed branched and parallel evolution patterns. ITH level decreased at higher sequencing depth of TDS than that measured by WES (28.1% vs 34.9%, P &lt; 0.01) but it remained unchanged upon additional samples analyzed. TDS of single tumor specimen detected an average of 70% the total mutations in HCC. Although more mutations were detected in cfDNA under TDS than WES, an average of 47.2% total HCC mutations uncovered by cfDNA suggested tissue outperform cfDNA and the latter may serve as alternative in profiling HCC genome. Consequently, TDS of single tumor tissue in 66 patients and cfDNAs in four unresectable HCCs identified 38.6% (26/66 and 1/4) patients bearing therapeutic targets.Conclusions: TDS of single tumor specimen could largely circumvent ITH to uncover mutations indicative of target therapy in HCC.

  Study EGAS00001002207

• whole genome sequence data of multifocal hepatocellular carcinoma

  Background &amp; Aims: The differentiation of distinct multifocal hepatocellular carcinoma (HCC): multicentric disease vs. intrahepatic metastases, in which the management and prognosis varies substantively, remains problematic. We aim to stratify multifocal HCC and identify novel diagnostic and prognostic biomarkers by performing whole genome and transcriptome sequencing, as part of a multi-omics strategy Methods: A complete collection of tumour and somatic specimens (intrahepatic HCC lesions, matched non-cancerous liver tissue and blood) were obtained from representative patientswith multifocal HCC exhibiting two distinct postsurgical courses. Whole-genome and transcriptome sequencing with genotyping were performed for each tissue specimen to con-trast genomic alterations, including hepatitis B virus integrations, somatic mutations, copy number variations, and structural variations. We then constructed a phylogenetic tree to visualise individual tumour evolution and performed functional enrichment analyses on select differentially expressed genes to elucidate biological processes involved in multifocal HCC development. Multi-omics data were integrated with detailed clinicopathological information to identify HCC biomarkers, which were further validated using a large cohort of HCC patients (n = 174).Results: The multi-omics pro?ling and tumour biomarkers could successfully distinguish the two multifocal HCC types, while accurately predicting clonality and aggressiveness. The dual speci?city protein kinase TTK, which is a key mitotic checkpoint regulator with links to p53 signaling, was further shown to be a promising overall prognostic marker for HCC in the large patient cohort.Conclusions: Comprehensive multi-omics characterisation of multifocal tumour evolution may improve clinical decision-making, facilitate personalised medicine, and expedite identi?cation of novel biomarkers and therapeutic targets in HCC.

  Study EGAS00001002338

• Aging and genome-wide patterns of DNA methylation in an African rainforest hunter-gathering population

  Aging is associated with widespread changes in genome-wide patterns of DNA methylation. Thousands of CpG sites whose tissue-specific methylation levels are strongly correlated with chronological age have been previously identified. However, the majority of these studies have focused primarily on cosmopolitan populations living in the developed world; it is not known if age-related patterns of DNA methylation at these loci are similar across a broad range of human genetic and ecological diversity. We investigated genome-wide methylation patterns using saliva derived DNA from a traditionally hunting and gathering African populations: the Baka of the western Central African rainforest, together with the ≠Khomani San of the South African Kalahari Desert. We identify hundreds of CpG sites whose methylation levels are significantly associated with age, thousands that are significant in a meta-analysis, and replicate trends previously reported in populations of non-African descent. We confirm that an age associated site in the gene ELOVL2 shows a remarkably congruent relationship with aging in humans, despite extensive genetic and environmental variation across populations. We also demonstrate that genotype state at methylation quantitative trait loci (meQTLs) can affect methylation trends at some known age-associated CpG sites. Our study explores the relationship between CpG methylation and chronological age in populations of African hunter-gatherers, who rely on different diets across diverse ecologies. While many age-related CpG sites replicate across populations, we show that considering common genetic variation at meQTLs further improves our ability to detect previously identified age associations.

  Study EGAS00001002226

• A targeted gene panel that covers coding, noncoding, and short tandem repeat regions improves the diagnosis of patients with neurodegenerative diseases

  Genetic testing for diagnosis of neurodegenerative diseases (NDs) is highly challenging because of genotype heterogeneity and overlapping clinical manifestations. Targeted-gene panels (TGPs), coupled with next-generation sequencing (NGS), can facilitate the profiling of a large repertoire of ND-related genes and help with differential diagnosis. Due to the technical limitations inherent in NGS and TGPs, short tandem repeat (STR) variations are often ignored. However, STR expansions are known to cause such NDs as Huntington’s disease and several types of spinocerebellar ataxias, such as SCA type 3 (SCA3).Here, we studied the clinical utility of a custom-made TGP that targets 199 NDs and 311 ND-associated genes on 118 undiagnosed patients. At least one known pathogenic or likely pathogenic variation was found in 54 (45.8%) patients; 27 (22.9%) of the 118 patients demonstrated clinical profiles that matched the variants. Based on the findings from the TGP, we refined the original diagnosis of 16 patients. A high concordance of single-nucleotide polymorphisms (99.9%) and small insertions and deletions (93.2%) were observed when comparing the results from TGP and whole-exome sequencing of four patients. We tested an in-house STR detection algorithm, and it reached a specificity of 0.88 and a sensitivity of 0.82 in our SCA3 cohort. This study also uncovered a trove of novel and recurrent variants that may enrich the repertoire of ND-related genetic markers. We propose that a combination of comprehensive TGPs and the bioinformatics analysis pipeline can improve the clinical diagnosis and accelerate the diagnosis of NDs.

  Study EGAS00001003737

• Human genomic and phenotypic synthetic data for the study of rare diseases

  The purpose of this dataset is to facilitate development of technical implementations for rare disease data integration, analysis, discovery, and federated access. This synthetic dataset includes clinical and genomic data from 6 rare disease cases. It consists of 18 whole genomes (6 index cases with their parents) which have genetic background based on public human data sequenced in the context of the Illumina Platinum initiative (Eberle, MA et al. (2017)) and made available by the HapMap project (https://www.genome.gov/10001688/international-hapmap-project). In each of the cases, real causative variants correlating with the phenotypic data provided were spiked-in. The cases included in this synthetic dataset correspond to the following type of disorders: CASE 1- Congenital myasthenic syndrome (Autosomal Dominant -de novo variant) CASE 2- Macular dystrophy (Autosomal Dominant) CASE 3- Muscular dystrophy (Autosomal Recessive-compound heterozygous variants) CASE 4- Mitochondrial disorder (Autosomal Recessive-consanguineous case - homozygous variant) CASE 5- Breast cancer (Autosomal Dominant) CASE 6- Similar as case 1 for patient matchmaking tests: Congenital myasthenic syndrome (Autosomal Dominant-de novo variant) For each case you will be able to download the following data: clinical information (phenopackets per individual and pedigree per family), raw genomic data (FASTQ and BAMs) and processed genomic data (vcfs). When using the data, the following should be acknowledged: the RD-Connect GPAP (https://platform.rd-connect.eu/), EC H2020 project EJP-RD (grant # 825575), EC H2020 project B1MG (grant # 951724) and Generalitat de Catalunya VEIS project (grant # 001-P-001647).

  Study EGAS00001005702

• Single cell RNA-seq mapping of nasal and tracheobronchial airways in human healthy volunteers

  The respiratory tract constitutes an elaborated line of defense based on a unique cellular ecosystem. Single-cell profiling methods enable the investigation of cell population distributions and transcriptional changes along the airways. We have explored cellular heterogeneity of the human airway epithelium in 10 healthy living volunteers by single-cell RNA profiling. 77,969 cells were collected by bronchoscopy at 35 distinct locations, from the nose to the 12th division of the airway tree. The resulting atlas is composed of a high percentage of epithelial cells (89.1%), but also immune (6.2%) and stromal (4.7%) cells with peculiar cellular proportions in different sites of the airways. It reveals differential gene expression between identical cell types (suprabasal, secretory, and multiciliated cells) from the nose (MUC4, PI3, SIX3) and tracheobronchial (SCGB1A1, TFF3) airways. By contrast, cell-type specific gene expression was stable across all tracheobronchial samples. Our atlas improves the description of ionocytes, pulmonary neuro-endocrine (PNEC) and brush cells, which are likely derived from a common population of precursor cells. We also report a population of KRT13 positive cells with a high percentage of dividing cells which are reminiscent of “hillock” cells previously described in mouse. Robust characterization of this unprecedented large single-cell cohort establishes an important resource for future investigations. The precise description of the continuum existing from nasal epithelium to successive divisions of lung airways and the stable gene expression profile of these regions better defines conditions under which relevant tracheobronchial proxies of human respiratory diseases can be developed.

  Study EGAS00001004082

• Neoadjuvant immunotherapy leads to pathological responses in MMR proficient and MMR deficient early stage colon cancers

  PD-1 plus CTLA-4 blockade is highly effective in advanced-stage, mismatch repair (MMR)-deficient (dMMR) colorectal cancers, yet not in MMR-proficient (pMMR) tumors. We postulated a higher efficacy of neoadjuvant immunotherapy in early-stage colon cancers. In the exploratory NICHE study (ClinicalTrials.gov: NCT03026140), patients with dMMR or pMMR tumors received a single dose of ipilimumab and two doses of nivolumab before surgery, the pMMR group with or without celecoxib. The primary objective was safety and feasibility; 40 patients with 21 dMMR and 20 pMMR tumors were treated, and 3 patients received nivolumab monotherapy in the safety run-in. Treatment was well tolerated and all patients underwent radical resections without delays, meeting the primary endpoint. Of the patients who received ipilimumab + nivolumab (20 dMMR and 15 pMMR tumors), 35 were evaluable for efficacy and translational endpoints. Pathological response was observed in 20/20 (100%; 95% exact confidence interval (CI): 86–100%) dMMR tumors, with 19 major pathological responses (MPRs, ≤10% residual viable tumor) and 12 pathological complete responses. In pMMR tumors, 4/15 (27%; 95% exact CI: 8–55%) showed pathological responses, with 3 MPRs and 1 partial response. CD8+PD-1+ T cell infiltration was predictive of response in pMMR tumors. These data indicate that neoadjuvant immunotherapy may have the potential to become the standard of care for a defined group of colon cancer patients when validated in larger studies with at least 3 years of disease-free survival data.

  Study EGAS00001004160

• The genomic landscape of pediatric acute lymphoblastic leukemia

  Acute lymphoblastic leukemia (ALL) is the most common pediatric cancer and arises from B- or T-lineage lymphocyte precursors. ALL comprises dozens of subtypes, and genomic analyses have largely been performed within these subtypes. Here we analyze pediatric ALL genomes across all subtypes, including 768 whole genomes, 1,729 exomes, and 1,889 transcriptomes in 2,754 patients. Most ALL subtypes harbor 4 or more driver alterations per sample, similar to adult cancers, despite low mutation burdens. Hyperdiploid B-ALL copy gains are likely acquired early and synchronously, with copy gains occurring before ultraviolet-induced mutations. By contrast, ultraviolet-induced mutations precede copy gains in iAMP21 B-ALL. Overall, we identified 378 putative ALL driver genes. Most driver alterations vary in prevalence across ALL subtypes, with B-ALL enriched for Ras and B-lineage-related alterations, and T-ALL enriched for PI3K, JAK, and cell cycle alterations. Most B-ALL (54.3%) and T-ALL (51.2%) samples bear at least one rare driver gene alteration (present in less than 2% of samples), including 30 putative novel cancer driver genes associated with ubiquitination, SUMOylation, non-coding, and other functions. Known or novel alterations associated with poor outcomes in specific subtypes include TBL1XR1 alterations in ETV6-RUNX1 patients, CREBBP in Ph-like-CRLF2, SETD2 in hyperdiploid, and PTEN in TAL1 patients. Intriguingly, DUX4 and KMT2A subtypes separate into CEBPA/FLT3- or NFATC4-expressing subgroups with potential clinical implications. Together, these results deepen understanding of ALL etiology and outcomes, and facilitate ALL model development.

  Study EGAS00001005250

• Sequencing-based counting and size profiling of plasma Epstein-Barr virus DNA enhance population screening of nasopharyngeal carcinoma.

  Circulating tumor-derived DNA testing for cancer screening has recently been demonstrated in a prospective study on identification of nasopharyngeal carcinoma (NPC) among 20,174 asymptomatic individuals. Plasma EBV DNA, a marker for NPC, was detected using real-time PCR. While plasma EBV DNA was persistently detectable in 97.1% of the NPCs identified, ∼5% of the general population had transiently detectable plasma EBV DNA. We hypothesized that EBV DNA in plasma of subjects with or without NPC may have different molecular characteristics. We performed target-capture sequencing of plasma EBV DNA and identified differences in the abundance and size profiles of EBV DNA molecules within plasma of NPC and non-NPC subjects. NPC patients had significantly higher amounts of plasma EBV DNA, which showed longer fragment lengths. Cutoff values were established from an exploratory dataset and tested in a validation sample set. Adopting an algorithm that required a sample to concurrently pass cutoffs for EBV DNA counting and size measurements, NPCs were detected at a positive predictive value (PPV) of 19.6%. This represented superior performance compared with the PPV of 11.0% in the prospective screening study, which required participants with an initially detectable plasma EBV DNA result to be retested within 4 weeks. The observed differences in the molecular nature of EBV DNA molecules in plasma of subjects with or without NPC were successfully translated into a sequencing-based test that had a high PPV for NPC screening and achievable through single time-point testing.

  Study EGAS00001002707

• Characterization of a human iPSC-derived islet differentiation model

  To explore the contribution of islet development to T2D pathogenesis, we characterised the transcriptomes of 3 human iPSC lines (from independent donors) differentiated along the islet development lineage. Whole transcriptome RNA-seq was performed at 7 developmental stages: definitive endoderm, primitive gut tube, posterior foregut (PF), pancreatic endoderm, endocrine progenitors, endocrine-like cells, and beta-like cells (BLC). Differentially-expressed (ΔExp) genes were assigned to the stage in which they were most upregulated (versus baseline iPSC profile), and used to assess the enrichment of T2D GWAS genes, as well as to predict the upstream transcription factors at each developmental stage.We found 9409 ΔExp genes across all stages, including known markers of islet development (NEUROG3, INS). Genes ΔExp in the most-differentiated (BLC) stage were significantly enriched for genes within the 99% credible sets of T2D GWAS loci. Despite this enrichment in BLC only, over 70% of genes mapping within the credible sets were ΔExp before this stage, highlighting the relevance of key T2D GWAS genes in islet development: e.g. expression of TCF7L2 peaked during the PF stage (log2FC=1.2; q=8.5×10-10). REST and LMNA are also highlighted by the analyses as having a potential function in islet development, considering their pattern of expression and those of the genes they regulate. Transcriptomic signatures derived from this iPSC differentiation model highlight T2D-associated GWAS loci, monogenic diabetes genes and potential master regulators of gene networks with plausible function in islet development. All these provide clues for developmental mechanisms contributing to T2D diabetes pathology.

  Study EGAS00001002721

• Differential Presence of Exons in Cell-Free DNA Reveals Different Patterns in Colorectal Cancer Between Metastatic and Non-Metastatic Patients

  Tumor genotyping during disease follow-up is gaining an increasing interest for monitoring and clinical management in colorectal cancer (CRC). Growing evidences suggest that single biopsies are inefficient for molecular profiling, due to clonal evolution and intratumoral heterogeneity of primary tumors and/or metastasis. Circulating free DNA (cfDNA) has emerged as an alternative source of genetic material that seems to be representative of the continuously changing tumor molecular features. Thus, the so-called 'liquid biopsy' also offers the additional advantages of accessibility and non-invasiveness for the patient. Novel strategies based on next-generation sequencing (NGS) of cfDNA are being currently developed, with most studies focusing on targeted deep sequencing panels of potential clinically actionable genes. Here we describe an alternative approach, consisting of exome-sequencing in cfDNA at a relatively shallower depth to provide a more general overview of the circulating genome. In this study, two groups of colorectal cancer patients, with either disseminated (M) or localized disease (N), were compared using whole-exome sequencing and a RNA-seq bioinformatics analysis pipeline. Thus, a set of 379 exons present at different levels in cfDNA from both groups were identified, giving rise to a completely new concept in this field, termed as 'Differential Presence of Exons'. Differentially present exons (DPEs) were subsequently used for M and N patients clustering and classification, allowing for the design of a predictive algorithm with encouraging preliminary results in a small subset of patients that could not be initially classified attending to the selection criteria (U group, unassignable).

  Study EGAS00001002687

• Neoantigen responses to Immunotherapy in Prostate Cancer

  To determine whether antigen-specific T cell responses can be elicited following treatment with ICT in cancers that have a low TMB, a clinical trial with ipilimumab (anti-CTLA-4) was conducted in 30 patients with metastatic castration-resistant prostate cancer (CRPC). We identified 2 distinct groups of patients divided by clinical outcomes, which consisted of a cohort with Favorable and Unfavorable outcomes. Compared to the Unfavorable cohort, the patients in the Favorable cohort were more likely to have higher intratumoral CD8 T cell density and/or IFN-γ response gene signature, including five patients who had antigen-specific T cell responses that were only detectable after ipilimumab therapy. Three patients had T cell responses against tumor-associated antigens (PSMA and PAP). Importantly, two patients with a relatively low TMB had T cell responses against unique neoantigens. Collectively, our data indicate that a subset of patients with metastatic prostate cancer benefit from ICTs, particularly those who have immunological biomarkers associated with effector T cell responses. Publication: Sumit K. Subudhi, Luis Vence, Hao Zhao, Jorge Blando, Shalini S. Yadav, Qing Xiong, Alexandre Reuben, Ana Aparicio, Paul G. Corn, Brian F. Chapin, Louis L. Pisters, Patricia Troncoso, Rebecca Slack Tidwell, Peter Thall, Chang-Jiun Wu, Jianhua Zhang, Christopher L. Logothetis, Andrew Futreal, James P. Allison and Padmanee Sharma, 2020, Neoantigen responses, immune correlates, and favorable outcomes after ipilimumab treatment of patients with prostate cancer. Science Translational Medicine, 01 Apr 2020:Vol. 12, Issue 537, eaaz3577DOI: 10.1126/scitranslmed.aaz3577. https://pubmed.ncbi.nlm.nih.gov/32238575/

  Study EGAS00001004050

• TRACERx 100: whole exome data of the first 100 TRACERx tumours

  BACKGROUNDTRACERx (TRAcking Cancer Evolution through therapy (Rx)) is a prospective cohort study designed to investigate intratumor heterogeneity (ITH) in relation to clinical outcome, and to determine the clonal nature of driver events and evolutionary processes in early stage non-small cell lung cancer (NSCLC). METHODSMultiregion high-depth whole-exome sequencing (M-seq) was performed on 100 early stage NSCLC tumors resected prior to systemic therapy. A total of 327 tumor regions were sequenced and analyzed to define evolutionary histories, obtain a census of clonal and subclonal events, and assess the relationship between ITH and recurrence-free survival (RFS). RESULTSWidespread ITH was observed for both somatic copy number alterations (median 48% [0.03-88%]) and mutations (median 30% [0.5-93%]). Driver mutations in EGFR, MET, BRAF and TP53 were almost always clonal. However, heterogeneous driver alterations occurring later in evolution were found in over 75% of tumors and were common in PIK3CA, NF1 and genes involved in chromatin modification and DNA response and repair. Genome doubling and ongoing dynamic chromosomal instability (CIN), illustrated by mirrored subclonal allelic imbalance, were identified as causes of ITH resulting in parallel evolution of driver copy number events, including amplifications of CDK4, FOXA1, and BCL11A. Elevated copy number heterogeneity was associated with shorter RFS (HR=4.9, P=0.00044), which remained significant in a multivariate analysis.CONCLUSIONSITH mediated through CIN, rather than point mutational heterogeneity, was associated with increased risk of relapse, supporting its value as a prognostic predictor, and the need to target this high-risk phenotype.

  Study EGAS00001002247

• Biallelic tumor suppressor loss and DNA repair defects in de novo small cell prostate cancer

  Small cell prostatic carcinoma (SCPC) is an aggressive pathology that is managed similar to small cell lung cancer. SCPC can evolve from prostatic adenocarcinoma in response to androgen deprivation therapy but in rare cases is present at initial cancer diagnosis. The molecular etiology of de novo SCPC is incompletely understood due to the scarcity of tumor tissue and the brief life expectancy of patients. Through a retrospective search of our regional oncology pharmacy database, we identified 18 patients diagnosed with de novo SCPC between 2004 and 2017. Ten patients had pure SCPC pathology and the remainder had some admixed adenocarcinoma foci, but all were treated with first line platinum-based chemotherapy. Median overall survival was 28 months. We performed targeted DNA sequencing, whole exome sequencing and mRNA profiling on formalin-fixed paraffin embedded archival tumor tissue. We observed frequent biallelic deletion and/or mutation of tumor suppressors TP53, RB1 and PTEN, similar to treatment related SCPC. Indeed, at the RNA level pure de novo SCPC closely resembled treatment-related SCPC. However, five patients had biallelic loss of DNA repair genes including BRCA1, BRCA2, ATM and MSH2/6: potentially underlying the high genomic instability of this rare disease variant. Two patients with pure de novo SCPC harbored ETS gene rearrangements to androgen-driven promoters, consistent with evolution of de novo SCPC from an androgen-driven ancestor. Overall, our results reveal a highly aggressive molecular landscape underlies this unusual pathologic variant, and propose opportunities for targeted therapy strategies in a disease with few treatment options.

  Study EGAS00001003007

• Single-cell DNA amplicon sequencing reveals clonal heterogeneity and evolution in T-cell acute lymphoblastic leukemia

  T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive leukemia that is most frequent in children and is characterized by the presence of few chromosomal rearrangements and 10 to 20 somatic mutations in protein-coding regions at diagnosis. The majority of T-ALL cases harbor activating mutations in NOTCH1 together with mutations in genes implicated in kinase signaling, transcriptional regulation or protein translation. To obtain more insight in the level of clonal heterogeneity at diagnosis and during treatment, we used single-cell targeted DNA sequencing with the Tapestri platform. We designed a custom ALL panel and obtained accurate single-nucleotide variant and small insertion-deletion mutation calling for 305 amplicons covering 110 genes in about 4400 cells per sample and time point. A total of 108,188 cells were analyzed for 25 samples of 8 T-ALL patients. We typically observed a major clone at diagnosis (&gt;35% of the cells) accompanied by several minor clones of which some were less than 1% of the total number of cells. Four patients had &gt;2 NOTCH1 mutations some of which present in minor clones, indicating a strong pressure to acquire NOTCH1 mutations in developing T-ALL cells. By analyzing longitudinal samples, we detected the presence and clonal nature of residual leukemic cells as well as clones with a minor presence at diagnosis that evolved to clinically relevant major clones at later disease stages. Single-cell DNA amplicon sequencing is a sensitive assay to detect clonal architecture and evolution in T-ALL.

  Study EGAS00001004440

• Open_Targets_020_Epigenomes_of_Cell_Lines

  There is currently a drive to establish cell based assay systems of greater human biological and disease relevance through the use of well characterised transformed cell lines, primary cells and complex cellular models (e.g. co-culture, 3D models). However, although the field is gaining valuable experience in running more non-standard &amp; complex cell assays for target validation and compound pharmacology studies, there is the lack of a systematic approach to determine if this expansion in cell assay models is reflected in increased human biological and disease relevance. The increasing wealth of publically available transcriptomic, and epigenome (ENCODE and Epigenome Roadmap) data represents an ideal reference mechanism for determining the relationship between cell types used for target &amp; compound studies to primary human cells and tissues from both healthy volunteers &amp; patients. The CTTV020 epigenomes of cell line project aims to generate epigenetic and transcriptomic profiles of cell lines and compare these with existing and newly generated reference data sets from human tissue and cell types. The aim is to identify assay systems which will provide greater confidence in translating target biology and compound pharmacology to patients. Multiple cell types commonly used within research have been grouped according to biology. Examples include erythroid, lung epithelial, hepatocyte cell types and immortalised models of monocyte / macrophage biology. 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 EGAS00001003136

• Structure and evolution of double minutes in diagnosis and relapse brain tumors

  Double minute chromosomes are extrachromosomal circular DNA fragments frequently found in brain tumors. To understand their evolution, we characterized the double minutes in paired diagnosis and relapse tumors from a pediatric high-grade glioma and four adult glioblastoma patients. We determined the full structures of the major double minutes using a novel approach combining multiple types of supporting genomic evidence. Among the double minutes identified in the pediatric tumor, only one carrying EGFR was maintained at high abundance in both samples, whereas two others were resent in only trace amounts at diagnosis but abundant at relapse, and the rest were found either in the relapse sample only or in the diagnosis sample only. For the EGFR-carrying double minutes, we found a secondary somatic deletion in all copies at relapse, after erlotinib treatment, but this was rare at diagnosis, suggesting potential resistance to the EGFR inhibitor. This mutation caused an in-frame RNA transcript to skip exon 16, a novel transcript absent in EST database and in about 700 brain tumors that we reviewed. We observed similar patterns involving longitudinal copy number shift of double minutes in another four pairs (diagnosis / relapse) of adult glioblastoma. Overall, in three of five paired tumor samples, we found that although the same oncogenes were amplified at diagnosis and relapse, they were amplified on different double minutes. Our results suggest that double minutes readily evolve, increasing tumor heterogeneity rapidly. Understanding patterns of double minute evolution can shed light on future therapeutic solutions to brain tumors carrying such variants.

  Study EGAS00001003212

• Oncogenic fate conversion by PRDM16s causes acute myeloid leukemia

  Oncogenic mutations confer aberrant replicative capacity to cells with little or no replicative capacity, generating cancer stem cells that perpetuate the tumor through extensive self-replication and differentiation blockade. However, whether oncogenes disrupt the cellular identity ofcancer stem cell by altering the developmental potential is unknown. Fate conversion has been demonstrated by ectopic expression of master transcriptional regulators, such as PU.1 and C/EBP alpha that confers myeloid cell fate to other cell types, and PRDM16 that confers brown adipose fate to white adipocytes or myoblasts. Here we show that a transcriptional regulator overexpressed in human myeloid malignancies, PRDM16s, causes oncogenic fate conversion, by transforming cells fated to form platelets and erythrocytes into myeloid leukemia-initiating cells (LICs). Prdm16s expression in hematopoietic progenitor cells caused a myelodysplastic syndrome (MDS)-like disease that progressed to acute myelogenous leukemia (AML). The myeloid diseases caused by Prdm16s exhibited expansion of megakaryocyte-erythroid progenitors (MEPs) but not granulocyte-macrophage progenitors. MEPs from Prdm16s-induced leukemia possessed LIC potential, and expression of Prdm16s in normal MEPs was sufficient to convert them to myeloid LICs and blocked megakaryocytic/erythroid potential. Prdm16s induced the expression of myeloid master regulators, including PU.1 and C/EBP alpha, by interacting with their super enhancers. Ablation of myeloid master regulators attenuated the myeloid potential and reinstalled the megakaryocytic/erythroid potential of leukemic-MEPs in mouse models and in human AML with PRDM16 rearrangement. Our study demonstrates that oncogenic Prdm16 expression converts the fate of MEPs to a malignant myeloid fate by activating myeloid master transcription factors.

  Study EGAS00001003235

• C3 SNPs and outcome after lung transplantation

  Background: Development of chronic rejection is still a severe problem and causes high mortality rates after lung transplantation (LTx). Complement activation is important in the development of acute rejection (AR) and bronchiolitis obliterans syndrome, with C3 as a key complement factor. Methods: Since a single nucleotide polymorphism (SNP) in the C3 gene (rs2230199) is associated with long-term outcome after solid organ transplantation, we investigated this SNP in relation to LTx. In addition, we looked at local production of C3 by analyzing bronchoalveolar lavage fluid (BALF) of LTx patients using isoelectric focusing (IEF). Results: We demonstrated the presence of C3 in BALF and showed that this is produced by the donor lung based on the genotype of SNP rs2230199. We also analyzed donor and patient SNP configurations and observed a significant association between the SNP configuration in patients and episodes of AR during 4-years follow-up. Survival analysis showed a lower AR-free survival in homozygous C3 slow patients (p=0.005). Furthermore, we found a significant association between the SNP configuration in donors and BOS development. Patients receiving a graft from a donor with at least one C3 fast variant for rs2230199 had an inferior BOS-free survival (p=0.044). Conclusions: In conclusion, our data indicate local C3 production by donor lung cells. In addition, a single C3 SNP present in recipients affects short-term outcome after LTx, while this SNP in donors has an opposite effect on long-term outcome after LTx. These results could contribute to an improved risk stratification after transplantation.

  Study EGAS00001003843

• Longitudinal_profiling_of_the_immune_response_to_Plasmodium_vivax_in_naive_hosts_by_RNA_sequencing

  Plasmodium vivax offers unique challenges for control and elimination, and may prove a tougher hurdle to overcome than Plasmodium falciparum. And yet compared to P. falciparum we know very little about the innate and adaptive immune responses that need to be harnessed to reduce disease and transmission. We recently generated a blood bank of a new clonal field isolate of P. vivax (PvW1) for human challenge studies and used systems immunology tools to track the host response throughout infection and convalescence. As part of this study, RNA-sequencing was used to resolve changes in whole blood gene expression through time in 6 volunteers (7-9 time-points per volunteer). In summary, these data show that P. vivax induces two distinct transcriptional programmes in whole blood during and after infection. During infection, transcriptional profiling reveals the rapid mobilisation of an emergency myeloid response, which leads to systemic inflammation and the recruitment of all major T cell subsets into lymphoid tissues. Six days after infection, this innate response subsides and a transcriptional signature of proliferation is revealed. This most likely represents widespread activation of lymphocytes, which return to the circulation after parasite clearance - transcriptional profiling of T cells at this time-point could therefore reveal the outcomes of critical cell-cell interactions that take place within the spleen during infection. 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 EGAS00001003847

• Comparison of the diagnostic yield of aCGH and NGS across different neurodevelopmental disorders

  Current recommendations for the genetic diagnosis of neurodevelopmental disorders (NDDs) do not include the use of next generation sequencing (NGS) and are still based on chromosomal microarrays, such as comparative genomic hybridization array (aCGH). The present study compares the diagnostic yield obtained by aCGH and clinical exome sequencing in NDD globally and its spectrum of disorders. To that end, 1412 patients clinically diagnosed with NDDs and studied with aCGH were classified into phenotype categories: Global developmental delay/intellectual disability (GDD/ID); autism spectrum disorder (ASD); and Other NDDs. These categories were further subclassified into subcategories, based on the most frequent accompanying signs and symptoms: isolated forms, forms with epilepsy; forms with micro or macrocephaly and syndromic forms. A subset of 245 patients of the 1412 were subjected to clinical exome sequencing. Diagnostic yield of aCGH and clinical exome sequencing, expressed as the number of solved cases, was compared for each phenotype category and subcategory. Clinical exome sequencing was superior than aCGH for all cases except for isolated ASD, with no cases solved by NGS. Globally, clinical exome sequencing solved 20% of cases (versus 5.7% by aCGH) and the diagnostic yield was highest for all forms of GDD/ID and lowest for Other NDDs (7.1% versus 1.4% by aCGH) and ASD (6.1% versus 3% by aCGH). In the majority of cases, diagnostic yield was higher in the phenotype subcategories than in the mother category. These results support the use of NGS as the first-tier test in the diagnostic algorithm of all NDDs followed by aCGH when necessary.

  Study EGAS00001004949

• Expanding the neurodevelopmental phenotypes of individuals with de novo KMT2A variants

  De novo loss-of-function (LoF) variants in the KMT2A gene are associated with Wiedemann-Steiner Syndrome (WSS). Recently, de novo KMT2A variants have been identified in sequencing studies of cohorts of individuals with neurodevelopmental disorders (NDDs). However, most of these studies lack the detailed clinical information required to determine whether those individuals have isolated NDDs or WSS (i.e. syndromic NDDs). We performed thorough clinical and neurodevelopmental phenotyping on six individuals with de novo KMT2A variants. From these data, we found that all six patients met clinical criteria for WSS and we further define the neurodevelopmental phenotypes associated with KMT2A variants and WSS. In particular, we identified a subtype of Autism Spectrum Disorder (ASD) in five individuals, characterized by marked rigid, repetitive and inflexible behaviours, emotional dysregulation, externalizing behaviours, but relative social motivation. To further explore the clinical spectrum associated with KMT2A variants, we also conducted a meta-analysis of individuals with KMT2A variants reported in the published literature. We found that de novo LoF or missense variants in KMT2A were significantly more prevalent than predicted by a previously established statistical model of de novo mutation rate for KMT2A. Our genotype-phenotype findings better define the clinical spectrum associated with KMT2A variants and suggest that individuals with de novo LoF and missense variants likely have a clinically unrecognized diagnosis of WSS, rather than isolated NDD or ASD alone. This highlights the importance of a clinical genetic and neurodevelopmental assessment for individuals with such variants in KMT2A.

  Study EGAS00001003521

• Healthspan and lifespan extension by fecal microbiota transplantation in progeroid mice

  The gut microbiome is emerging as a key regulator of several metabolic, immune and neuroendocrine pathways. Gut microbiome deregulation has been implicated in major conditions such as obesity, type 2 diabetes, cardiovascular disease, non-alcoholic fatty acid liver disease or cancer, but its precise role in aging remains to be elucidated. Here, we characterize the gut microbiome profile of accelerated aging and show that its external modulation is sufficient to extend healthspan and lifespan. We found that two different mouse models of progeria present with intestinal dysbiosis, which is characterized, among other alterations, by an increment in proteobacteria, cyanobacteria and a loss in verrucomicrobia, whereas long-lived humans (i.e., centenarians) exhibit a remarkable increase in verrucomicrobia and a reduction in proteobacteria. Fecal microbiota transplantation proved to be effective to enhance healthspan and lifespan in both progeroid mouse models and, more importantly, the solely transplantation with the verrucomicrobia Akkermansia muciniphila was sufficient to exert beneficial effects. Our results demonstrate that intestinal dysbiosis is a phenotype associated with accelerated aging and its correction provides health benefits. Moreover, metabolomic analysis of ileal content points to the restoration of secondary bile acids as a possible mechanism for the beneficial outcome of reestablishing a healthy microbiome. Our results suggest the existence of a link between gut aging and the microbiota, and can help to gain insight into the rationale for microbiome-based interventions against age-related diseases. Mouse sequence data have been deposited in ENA (https://www.ebi.ac.uk/ena) under accession number PRJEB34214

  Study EGAS00001003656

• Detection of ctDNA in Plasma of Patients with Clinically Localised Prostate Cancer is Associated with Rapid Disease Progression

  Purpose: DNA originating from degenerate tumour cells can be detected in the circulation in many tumour types, where it can be used as a marker of disease burden as well as to monitor treatment response. Although circulating tumour DNA (ctDNA) measurement has prognostic/predictive value in metastatic prostate cancer, its utility in localised disease is unknown. Patients and Methods: We performed whole genome sequencing of tumour-normal pairs in eight patients with clinically localised disease undergoing prostatectomy, identifying high confidence genomic aberrations. A bespoke DNA capture and amplification panel against the highest prevalence, highest confidence aberrations for each individual was designed and used to interrogate ctDNA isolated from plasma prospectively obtained pre- and post- (24 hours and 6 weeks) surgery. In a separate cohort (n=189) we identified the presence of ctDNA TP53 mutations in pre-operative plasma in a retrospective cohort, and determined its association with biochemical- and metastasis-free survival. Results: Tumour variants in ctDNA were positively identified pre-treatment in two of eight patients, which in both cases remained detectable postoperatively. Patients with tumour variants in ctDNA had extremely rapid disease recurrence and progression compared to those where variants could not be detected. In terms of aberrations targeted, single nucleotide and structural variants outperformed indels and copy number aberrations. Detection of ctDNA TP53 mutations was associated with a significantly shorter metastasis-free survival (6.2 vs. 9.5 years (HR 2.4; 95% CIs 1.2-4.8, p = 0.014)). Conclusions: CtDNA is uncommonly detected in localised prostate cancer but its presence portends more rapidly progressive disease.

  Study EGAS00001004636

• Molecular classification improves risk assessment in adult B-lineage ALL: Patients on the international UKALLXII-ECOG2993 trial.

  Acute lymphoblastic leukemia (ALL) in adults is a rare disease with outcomes that remain significantly inferior to those of children with ALL. While genomic classification of pediatric B-lineage ALL (B-ALL) has led to improved risk assignment and therapeutic paradigms, similar advances in adults have been lacking. The international UKALLXII/E2993 (NCT00002514) trial accrued over 2000 adult (17-64 years old) patients over 13 years, including 1229 BCR/ABL negative B-ALLs. While 93% of B-ALL achieved a remission, 41% relapsed at a median time of 13 months (Range) and overall survival (OS) at 5 years was 42% (CI). Older age and higher white blood cell (WBC) count at presentation negatively affected OS. Patients with a sibling donor assigned to allogeneic transplantation had superior OS compared to patients without donor. Randomization of patients without a donor showed an advantage of chemotherapy over autologous transplantation in first remission. Transcriptome sequencing of 94% of E2993 patients enabled genomic subtyping. We found that one-third of patients belonged to favorable genotypes (DUX4-rearranged, ETV6-RUNX1, TCF3-PBX1, PAX5 P80R, high hyperdiploidy). Of 63.4% of patients classified as high-risk based age and WBC count at diagnosis, almost 40% harbored subtype-defining genetic associations associated with favorable- or intermediate outcomes. These patients did significantly better than expected from their protocol-defined risk assignment. Our data in a large adult B-ALL cohort treated on a single international trial show the important prognostic value of genomic analysis and the potential benefit for many of the novel genotypes.

  Study EGAS00001004638

• Sporadic and endemic Burkitt lymphoma have frequent FOXO1 mutations but distinct hotspots in the AKT recognition motif

  FOXO1 has an oncogenic role in adult germinal center derived lymphomas, in which mutations, predominately within the AKT recognition motif, cause nuclear retention of FOXO1 resulting in increased cell proliferation. To determine the prevalence and distribution of FOXO1 mutations in pediatric Burkitt lymphoma (BL), we sequenced a large number of sporadic and endemic BL patient samples. We report a high frequency of FOXO1 mutations in both sporadic and endemic BL at diagnosis, occurring in 23/78 (29%) and 48/89 (54%) samples respectively, as well as 8/16 (50%) cases at relapse. Mutations of T24 were the most common in sporadic BL but were rare in endemic cases, in which mutations of residue S22, also within the AKT recognition motif, were the most frequent. FOXO1 mutations were almost always present in the major tumor cell clone but were not associated with outcome. Analysis of other recurrent mutations reported in BL revealed that FOXO1 mutations were associated with mutations of DDX3X and ARID1A, but not MYC, TCF3/ID3 or members of the phosphatidylinositol 3’ OH kinase (PI3K) signaling pathway. We further show common nuclear retention of the FOXO1 protein, irrespective of mutation status, suggesting alternative unknown mechanisms for maintaining FOXO1 transcriptional activity in BL. CRISPR/Cas9 knockout of FOXO1 in an endemic cell line produced a significant decrease in cell proliferation, supporting an oncogenic role for FOXO1 in endemic BL. Thus, FOXO1 is frequently mutated in both sporadic and endemic BL and may offer a potential therapeutic target for pediatric BL patients worldwide.

  Study EGAS00001003719

• Next-Generation Sequencing of AV Nodal Reentry Tachycardia patients

  Atrioventricular nodal reentry tachycardia (AVNRT) is the most common form of regular paroxysmal supraventricular tachycardia. This arrhythmia affects women twice as frequently as men, and is often diagnosed in patients below 40 years of age. Familial clustering, early onset of symptoms, and lack of structural anomaly indicate involvement of genetic factors in AVNRT pathophysiology. We hypothesized that AVNRT patients have a high prevalence of variants in genes that are highly expressed in the atrioventricular conduction axis of the heart and potentially involved in arrhythmic diseases. Next-generation sequencing of 67 genes was applied to the DNA profile of 298 AVNRT patients and 10 AVNRT family members using HaloPlex Target Enrichment System.In total, we identified 229 variants in 60 genes; 215 missenses, four frame shifts, four codon deletions, three missense and splice sites, two stop-gain variants, and one start-lost variant. Sixty-five of these were not present in the Exome Aggregation Consortium (ExAC) database. Furthermore, we report two AVNRT families with co-segregating variants. Seventy-five of 284 AVNRT patients (26.4%) and three family members to different AVNRT probands had one or more variants in genes affecting the sodium handling. Fifty-four out of 284 AVNRT patients (19.0%) had variants in genes affecting the calcium handling of the heart. We furthermore find a large proportion of variants in the HCN1-4 genes. We did not detect a significant enrichment of rare variants in the tested genes.This could be an indication that AVNRT might be an electrical arrhythmic disease with abnormal sodium and calcium handling.

  Study EGAS00001002745

• Dynamics of genome architecture and chromatin function during human B cell differentiation and neoplastic transformation

  Despite recent advances, the dynamics of genome architecture and chromatin function during human cell differentiation and its potential reorganization upon neoplastic transformation remains poorly characterized. Here, we integrate in situ Hi-C and nine additional omic layers to define and biologically characterize the dynamic changes in three-dimensional (3D) genome architecture across normal B cell differentiation and in neoplastic cells from different subtypes of chronic lymphocytic leukemia (CLL) and mantle cell lymphoma (MCL) patients. Beyond conventional active (A) and inactive (B) compartments, an integrative analysis of Hi-C data uncovers a highly-dynamic intermediate compartment enriched in poised and polycomb-repressed chromatin. During B cell development, we detect that 28% of the compartments change at defined maturation stages and mostly involve the intermediate compartment. The transition from naive to germinal center B cells is associated with widespread chromatin activation, which mostly reverts into the naive state upon further maturation of germinal center cells into memory B cells. The analysis of CLL and MCL neoplastic cells points both to entity and subtype-specific alterations in chromosome organization. Remarkably, we observe that large chromatin blocks containing key disease-specific genes alter their 3D genome organization. These include the inactivation of a 2Mb region containing the EBF1 gene in CLL and the activation of a 6.1Mb region containing the SOX11 gene in clinically aggressive MCL. This study indicates that 3D genome interactions are extensively modulated during normal B cell differentiation and that the genome of B cell neoplasias acquires a tumor-specific 3D genome architecture.

  Study EGAS00001004763

• Transcriptional profiling of ovarian tumours and cell lines

  Purpose: Endogenous retroviruses (ERVs) play a role in a variety of biological processes, including embryogenesis and cancer. DNA methyltransferase inhibitors (DNMTi)-induced ERV expression triggers interferon responses in ovarian cancer cells via the viral sensing machinery. Baseline expression of ERVs also occurs in cancer cells, though this process is poorly understood and previously unexplored in epithelial ovarian cancer (EOC). Here, the prognostic and immunomodulatory consequences of baseline ERV expression was assessed in EOC. Experimental Design: ERV expression was assessed using EOC transcriptional data from TCGA and from an independent cohort (HH), as well as from untreated or DNMTi-treated EOC cell lines. LASSO logistic regression defined an ERV expression score to predict patient prognosis. Immunohistochemistry (IHC) was conducted on the HH cohort. Combination of DNMTi treatment with γδ T cells was tested in vitro, using EOC cell lines and patient-derived tumor cells. Results: ERV expression was found to define clinically relevant subsets of EOC patients. An ERV prognostic score was successfully generated in TCGA and validated in the independent cohort. In EOC patients from this cohort, a high ERV score was associated with better survival (log-rank p=0.0009) and correlated with infiltration of CD8+PD1+ T cells (r=0.46, p=0.0001). In vitro, baseline ERV expression dictates the level of ERV induction in response to DNMTi. Manipulation of ERV expression by DNMTi resulted in improved EOC cell killing by cytotoxic immune cells. Conclusions: These findings uncover the potential for baseline ERV expression to robustly inform EOC patient prognosis, influence tumor immune infiltration and affect anti-tumor immunity.

  Study EGAS00001004814

• Colorectal cancer organoids expressing BRAF (fusion) genes

  Fusion genes may be oncogenic drivers and potential targets for personalized therapies in a variety of cancer types. The BRAF gene is frequently involved in oncogenic fusions, with fusion frequencies of 0.2-3% throughout different cancers. BRAF can be fused to a wide variety of genes, however, BRAF fusions rarely occur in the exact same gene configuration, making assessment of the fusion relevance and decision for drug treatment challenging. Here, we devised a colorectal cancer organoid-based platform to functionally characterize diverse BRAF fusions, containing various partner genes (AGAP3, DLG1 and TRIM24). We compared these BRAF fusions with respect to cellular behaviour, downstream signaling activation and response to targeted therapies. We found that 5’ partner choice of BRAF affects cellular localization and intracellular signaling capacities of the fusion genes. The DLG1-BRAF fusion gene showed distinct localization to the plasma membrane and exhibited increased levels of MAPK pathway activation under unperturbed conditions. Furthermore, the different BRAF fusions showed varying sensitivities to the targeted inhibition of ERK and BRAF, with the DLG1-BRAF fusion being the most sensitive. Additionally, RNA-sequencing identified distinct subsets of genes affected by the DLG1-BRAF fusion gene. Importantly, all fusion genes conveyed resistance to the clinically relevant EGFR/HER2/HER4-inhibitor afatinib, suggesting that BRAF fusions should be screened alongside other MAPK pathway alterations to identify mCRC patients amenable to cetuximab treatment. In summary, we developed a platform to efficiently assess molecular and cellular effects of fusion genes, and revealed that differential drug responses and distinct gene expression profiles can be induced by different BRAF fusions .

  Study EGAS00001003558

• MYO5B mutations in Pheochromocytoma/Paraganglioma promote cancer progression

  Identification of additional cancer-associated genes and secondary mutations driving the metastatic progression in pheochromocytoma and paraganglioma (PPGL) is important. We recently reported novel recurrent nonsynonymous mutations in the MYO5B gene in metastatic PPGL. Here, we explored the functional impact of these MYO5B mutations, and analyzed MYO5B expression in primary PPGL tumor cases in relation to mutation status.Immunohistochemistry and mRNA expression analysis in 30 primary PPGL tumors revealed an increased MYO5B expression in metastatic compared to non-metastatic tumor cases. In addition, subcellular localization of MYO5B protein was altered from cytoplasmatic to membranous in some metastatic tumors, and the strongest and most abnormal expression pattern was observed in a PGL case harboring a MYO5B:p.G1611S mutation. Mutation analysis of additional PPGL tumors revealed additional, and hence recurrent, mutations in the paralog MYO5A. Three missense mutations in MYO5B were functionally studied using site directed mutagenesis and stable transfection into human neuroblastoma cells (SK-N-AS) and embryonic kidney cells (HEK293 cell line). Functional analysis of the mutants verified a significant increased proliferation rate in all MYO5B mutated clones (p.L587P, p.G1611S, and p.R1641C). The two somatically derived mutations p.L587P and p.G1611S were also found to increase the migration rate, and all three mutants affected intracellular transport by increased transferrin uptake. Downstream expression analysis of MYO5B mutants demonstrated a transcriptional activation of genes involved in glucose metabolism, proliferation, homeostasis, and migration. Our study validates the functional role of novel MYO5B mutations in proliferation and migration, and suggest the MYO5-pathway to be involved the malignant progression of PPGL tumors.

  Study EGAS00001003991

• Phylogenetic analysis of treatment-naive metastases using whole exome and genome sequencing data

  Metastases are responsible for the majority of cancer related deaths and are often difficult to treat successfully. Even for metastases that occur subsequent to treatment, we do not know whether such relapsing metastases were originally heterogeneous or if the observed heterogeneity is a consequence of therapy. Given that treatment can influence evolutionary dynamics by selecting resistant clones, imposing bottlenecks on cancer cell populations, and even inducing novel somatic mutations, characterizing the prior standing variation among metastases remains an important goal to predict initial therapeutic response. To quantify the heterogeneity at clinical presentation of advanced disease, we surveyed the literature for patients in which at least two treatment-naïve metastases underwent genome/exome-wide sequencing. Across all cancer types surveyed, only 18 subjects were found to fulfill this requirement. Including previously unpublished data from two subjects, we analyzed data of 74 untreated metastases and inferred cancer phylogenies. Putative driver gene mutations were acquired at a 2-fold higher rate on the trunk of all metastases than on branches. Mutations in driver genes along the trunk were strongly enriched for predicted functional consequences. Using a stochastic mathematical model, we find that driver gene heterogeneity among metastases mostly occurs in slowly growing cancers for highly advantageous driver gene mutations. Seeding efficiency only weakly affects the heterogeneity among metastases. These results explain why functional driver gene heterogeneity is uncommon in advanced disease prior to treatment, thus providing optimism for future therapies that depend on mutations across metastases.

  Study EGAS00001002777