23014 results for "ega"

in 79.22 milliseconds.

• cfDNA dataset from the urine supernatant of ovarian cancer patients and healthy controls

  Background Ovarian cancer is the deadliest gynecological malignancy worldwide, due to frequent diagnosis at advanced stage. Simple and non-invasive methods for earlier and accurate detection are urgently needed. Here, the presence of tumor-derived DNA in home-collected urine, cervicovaginal self-samples and clinician-taken cervical scrapes of ovarian cancer patients was explored by DNA methylation and copy number analysis. Methods A total of 428 samples (urine, cervicovaginal self-samples, and clinician-taken cervical scrapes) of 110 healthy controls and 54 patients with benign (n=25) or malignant (n=29) ovarian masses were analyzed. Different urine fractions were examined (full void urine, urine supernatant, and urine sediment). All samples were tested for 12 methylation markers by quantitative methylation-specific PCR. Shallow whole-genome sequencing was performed to detect copy number aberrations and verify the presence of tumor-derived DNA. Results Full void urine was most discriminatory between healthy controls and ovarian cancer patients (C2CD4D, p=0.008; CDO1, p=0.022; MAL, p=0.008), followed by cervical scrapes (C2CD4D, p=0.001; CDO1, p=0.004). A significant difference between benign and malignant ovarian masses was only observed for GHSR in the urine sediment (p=0.024). Methylation levels in cervicovaginal self-samples did not discriminate between healthy controls, benign and malignant ovarian masses. Copy number changes were identified in 17% of urine supernatant samples and smaller fragment sizes were seen in urine supernatant samples with a high tumor fraction. Conclusions This study demonstrates the presence of ovarian cancer-derived DNA in home-collected urine. Additional studies are warranted to further explore the clinical applicability of this approach.

  Study EGAS00001007238

• Multi-omic analysis of the tumor microenvironment shows clinical correlations in Ph1 study of atezolizumab +/- SoC in MM

  Multiple myeloma (MM) remains incurable, and treatment of relapsed/refractory (R/R) disease is challenging. There is an unmet need for more targeted therapies in this setting; deep cellular and molecular phenotyping of the tumor and microenvironment in MM could help guide such therapies. This phase 1b study (NCT02431208) evaluated the safety and efficacy of the anti-programmed death-ligand 1 monoclonal antibody atezolizumab (Atezo) alone or in combination with the standard of care (SoC) treatments lenalidomide (Len) or pomalidomide (Pom) and/or daratumumab (Dara) in patients with R/R MM. Study endpoints included incidence of adverse events (AEs) and overall response rate (ORR). A novel unsupervised integrative multi-omic analysis was performed using RNA sequencing, mass cytometry immunophenotyping, and proteomic profiling of baseline and on-treatment bone marrow samples from patients receiving Atezo monotherapy or Atezo+Dara. A similarity network fusion (SNF) algorithm was applied to preprocessed data. Eighty-five patients were enrolled. Treatment-emergent deaths occurred in 2 patients; both deaths were considered unrelated to study treatment. ORRs ranged from 11.1% (Atezo+Len cohorts, n=18) to 83.3% (Atezo+Dara+Pom cohort, n=6). High-dimensional multi-omic profiling of the tumor microenvironment and integrative SNF analysis revealed novel correlations between cellular and molecular features of the tumor and immune microenvironment, patient selection criteria, and clinical outcome. Atezo monotherapy and SoC combinations were safe in this patient population and demonstrated some evidence of clinical efficacy. Integrative analysis of high dimensional genomics and immune data identified novel clinical correlations that may inform patient selection criteria and outcome assessment in future immunotherapy studies for myeloma.

  Study EGAS00001007286

• Whole genome sequencing of 108 epileptic patients from CENet cohort

  Epileptic patients have been sequenced as part of the CENet cohort which is composed of patients with Genetic Generalized Epilepsy (GGE) or Non-Acquired Focal Epilepsy (NAFE) collected in CHUM Research Center in Montreal. Patients were diagnosed by epileptologists. The clinical epilepsy phenotype was classified according to the current classification by the International League against Epilepsy (ILAE). Libraries preparation and whole-genome sequencing: gDNA was cleaned up using ZR-96 DNA Clean & ConcentratorTM-5 Kit (Zymo) prior to being quantified using the Quant-iTTM PicoGreen dsDNA Assay Kit (Life Technologies) and its integrity assessed on agarose gels. Libraries were generated using the TruSeq DNA PCR-Free Library Preparation Kit (Illumina) according to the manufacturer’s recommendations. Libraries were quantified using the Quant-iTTM PicoGreen dsDNA Assay Kit (Life Technologies) and the Kapa Illumina GA with Revised Primers-SYBR Fast Universal kit (Kapa Biosystems). Average size fragment was determined using a LabChip GX (PerkinElmer) instrument. The libraries were denatured in 0.05N NaOH and diluted to 8pM using HT1 buffer. The clustering was done on a Illumina cBot and the flowcell was ran on a HiSeq 2500 for 2x125 cycles (paired-end mode) using v4 chemistry and following the manufacturer's instructions. A phiX library was used as a control and mixed with libraries at 0.01 level. Bioinformatics: The Illumina control software was HCS 2.2.58, the real-time analysis program was RTA v. 1.18.64. Program bcl2fastq v1.8.4 was used to demultiplex samples and generate fastq reads. The filtered reads were aligned to reference Homo_sapiens assembly b37. Each readset was aligned to create a .cram file.

  Study EGAS00001007507

• Evaluation of somatic mutations in urine samples as a non-invasive method for the detection and molecular classification of endometrial cancer

  Purpose: Current diagnostic methods for endometrial cancer lack specificity, leading to many women undergoing invasive procedures. The aim of this study was to evaluate somatic mutations in urine to accurately discriminate endometrial cancer patients from controls. Experimental Design: Overall, 72 samples were analyzed using next-generation sequencing with molecular identifiers targeting 47 genes. We evaluated urine supernatant samples from women with endometrial cancer (n=19) and age-matched controls (n=20). Cell pellets from urine and plasma samples from seven cases were sequenced; further, we also evaluated paired tumor samples from all cases. Finally, immunohistochemical markers for molecular profiling were evaluated in all tumor samples. Results: Overall, we were able to identify mutations in DNA from urine supernatant samples in 100% of endometrial cancers. In contrast, only one control (5%) showed variants at a variant allele frequency (VAF)≥2% in the urine supernatant samples. The molecular classification obtained by using tumor samples and urine samples showed good agreement. Analyses in paired samples revealed a higher number of mutations and VAFs in urine supernatants than in urine cell pellets and blood samples. Conclusions: Evaluation of somatic mutations using urine samples may offer a user-friendly and reliable tool for endometrial cancer detection and molecular classification. The diagnostic performance for endometrial cancer detection was very high, and cases could be molecularly classified using these noninvasive and self-collected samples. Additional multi-center evaluations using larger sample sizes are needed to validate the results and understand the potential of urine samples for the early detection and prognosis of endometrial cancer.

  Study EGAS00001007396

• Regenerative effects of MSC treatment on busulfan-induced small intestine damaged organoids.

  Allogenic hematopoietic stem cell transplantation (HSCT) is a curative treatment for leukemia and a range of non-malignant disorders. The success of the therapy is hampered by occurrence of acute graft-versus-host disease (aGvHD); an inflammatory response damaging recipient organs, with gut, liver, and skin being the most susceptible. Intestinal GvHD injury is often a life-threatening complication in patients unresponsive to steroid treatment. Second-line available therapies are immunosuppressants or mesenchymal stromal/stem cell (MSC) infusions. Data from our institution and others demonstrate rescue of approximately 40-50% of patients suffering from aGvHD with mesenchymal stromal/stem cells (MSCs) and minor side effects. Although promising, the better understanding of MSC mode of action and patient response to MSC-based therapy is essential to improve this lifesaving treatment. Here, we developed a novel 3D co-culture model of human small intestinal organoids and MSCs, which allows to study the regenerative effects of MSCs on intestinal epithelium in a more physiologically relevant setting than existing in vitro systems. Using this model we mimicked chemotherapy-mediated damage of the intestinal epithelium. The treatment with busulfan, the chemotherapeutic commonly used as conditioning regiment before the HSCT, affected pathways regulating EMT, proliferation, and apoptosis in small intestinal organoids, as shown by transcriptomic and proteomic analysis. The co-culture of busulfan treated intestinal organoids with MSCs reversed the effects of busulfan on the transcriptome and proteome of intestinal epithelium, which we also confirmed by functional evaluation of proliferation and apoptosis. Collectively, we demonstrate that our novel in vitro co-culture system is a new valuable tool to facilitate the investigation of the molecular mechanisms behind the therapeutic effects of MSCs on damaged intestinal epithelium.

  Study EGAS00001007432

• Ultra-fast deep-learned pediatric CNS tumor classification.

  The primary treatment of CNS tumors starts with a neurosurgical resection in order to obtain tumor tissue for diagnosis and to reduce tumor load and mass effect. The neurosurgeon has to decide between radical resection versus a more conservative strategy to prevent surgical morbidity. The prognostic impact of a radical resection varies between tumor types. However due to a lack of pre-operative tissue-based diagnostics, limited knowledge of the precise tumor type is available at the time of surgery. Current standard practice includes preoperative imaging and intraoperative histological analysis, but these are not always conclusive. After surgery, histopathological and molecular tests are performed to diagnose the precise tumor type. The results may indicate that an additional surgery is needed or that the initial surgery could have been less radical. Using rapid Nanopore sequencing, a sparse methylation profile can be directly obtained during surgery, making it ideally suited to enable intraoperative diagnostics. We developed a state-of-the-art neural-network approach called Sturgeon, to deliver trained models that are lightweight and universally applicable across patients and sequencing depths. We demonstrate our method to be accurate and fast enough to provide a correct diagnosis with as little as 20 to 40 minutes of sequencing data in 45 out of 49 pediatric samples, and inconclusive results in the other four. In four intraoperative cases we achieved a turnaround time of 60-90 minutes from sample biopsy to result; well in time to impact surgical decision making. We conclude that machine-learned diagnosis based on intraoperative sequencing can assist neurosurgical decision making, allowing neurological comorbidity to be avoided or preventing additional surgeries.

  Study EGAS00001007475

• RNASeq Neoantigen Immunogenicity Landscape

  Immune checkpoint blockade (ICB) therapy is a cornerstone of oncologic treatment for patients with advanced stage non-small cell lung cancer (NSCLC) and other malignancies. Neoantigen immunoediting drives ICB efficacy, yet the fundamental physiochemical characteristics of neoantigens and how neoantigen immunogenicity shapes treatment response remains poorly understood. To help address these questions, a prospective clinical trial of NSCLC patients treated with nivolumab was conducted. We assessed genomic alterations in tumors from 58 patients and performed large-scale neoepitope immunogenicity analyses, before and during treatment (CheckMate153, CA209-153). Tumors were analyzed by whole-exome and transcriptome sequencing. In responding patients, loss of mutation and neoantigen burden early during therapy associated with clinical benefit. We evaluated the immunogenicity of 1,453 candidate neoantigens and identified 502 neopeptides that bound to MHC I and 196 neopeptides that were immunogenically recognized by T cells in the setting of nivolumab treatment. These T cell reactive neoantigens were differentially present in clonal populations that underwent distinctive evolutionary trajectories across responders and nonresponders. Mapping these neoantigens to tumor clonal dynamics and clinical response revealed strong selection against immunogenic neoantigen harboring clones compared to non-immunogenic clones. Using this large collection of neoantigens, we identified position specific amino acid features related to immunogenicity, which we used to develop and validate an immunogenicity score. Changes in the genomic and neoantigen immunogenicity landscapes were accompanied by temporal changes in the tumor microenvironment. Nivolumab-induced microenvironmental evolution in NSCLC shared some similarities with that in melanoma, yet critical differences in immunologic programs were apparent from comparative network analysis between tumor types. This study provides unprecedented molecular portraits of the neoantigen landscapes underlying nivolumab’s mechanism of action.

  Study EGAS00001007509

• MethylationEPIC BeadChip samples of pre- and post-5ZA-treated head and neck cancer patients refractory to anto-PD-1 therapy

  Background Immune checkpoint blockade (ICB) is an effective treatment in a subset of patients diagnosed with head and neck squamous cell carcinoma (HNSCC); however, most patients are refractory. Methods In a non-randomized, open-label phase 1b clinical trial, participants with recurrent and/or metastatic (R/M) HNSCC were treated with low-dose 5’azacytidine (5’aza) daily for either 5- or 10-days in combination with fixed doses of durvalumab and tremelimumab after progression on ICB therapy. The primary objective was to assess the biologically effective dose of 5’aza as determined by molecular changes in paired baseline and on-treatment tumor biopsies; the secondary objective was safety. Results Thirty-eight percent (3/8) of participants with evaluable paired tissue samples had a greater than 2-fold increase from baseline in interferon-gamma (IFN-ɣ) signature and CD274 (programmed cell death protein 1 ligand, PD-L1) expression within the tumor microenvironment (TME), which was associated with increased CD8+ T cell infiltration and decreased infiltration of CD4+ T regulatory cells. The mean neutrophil-to-lymphocyte ratio (NLR) decreased by greater than 50%, from 14.2 (SD 22.6) to 6.9 (SD 5.2). Median overall survival (OS) was 16.3 months (95% CI 1.9, NA) and 2-year OS rate was 24.7% (95% CI: 4.5%, 53.2%) with 58% (7/12) of treated participants demonstrated prolonged OS of greater than 12 months. Conclusion Our findings suggest that low-dose 5’aza can reprogram systemic host immune responses and the local TME to increase IFN-ɣ and PD-L1 expression. The increased expression of these established biomarkers correlated with prolonged OS upon re-challenge with ICB.

  Study EGAS00001007998

• Multi-Layered Molecular Profiling Informs the Diagnosis and Targeted Therapy of Desmoplastic Small Round Cell Tumor

  Desmoplastic small round cell tumor (DSRCT) is an ultra-rare soft tissue sarcoma with dismal prognosis and limited therapeutic options. Here, we identified individual therapeutic targets in a cohort of 30 DSRCT patients by use of multilayered molecular profiling (genomic, transcriptomic, and (phospho-)proteomic analyses) within the DKFZ NCT DKTK MASTER program. 28 of the 30 patients received at least one personalized treatment recommendation. In line with the observation that DSRCTs are an entity with a low mutational load, most of the treatment recommendations were based on target gene or -protein expression. 13 patients received a recommended therapy, resulting in four partial responses, four stable diseases, and five progressive diseases. Eight of 13 patients received the multi tyrosine kinase inhibitor pazopanib, with mixed response. Based on high ERBB2 expression levels, two patients received the antibody drug conjugate trastuzumab deruxtecan (T-DXd), showed exceptional, long-lasting responses and are currently (8/2024) still on treatment and progression-free. We detected high ERBB2 protein expression in all DSRCT samples with proteome data (n = 9), which suggests that ERBB2 may be a particularly useful therapeutic target for DSRCT patients. For one of the two exceptional responders we did not detect a substantial ERBB2/HER2 overexpression, which may indicate that T-DXd could also be efficacious in DSRCT patients with low levels of HER2, as also described for breast cancer patients. In conclusion, multi-omic molecular profiling revealed personalized therapeutic targets in patients with DSRCT with HER2 targeting, with T-DXd as a promising therapeutic option for this ultra-rare and difficult-to-treat tumor entity.

  Study EGAS00001007934

• Transcriptomics-driven classification of ALAL

  Acute leukemia of ambiguous lineage (ALAL) is a rare, poor-prognosis subtype of acute leukemia that cannot be assigned to a single hematopoietic lineage. Although ALAL patients are typically treated with acute myeloid leukemia (AML) or acute lymphoblastic leukemia (ALL) regimens, optimal treatment choice is hindered by their lineage ambiguity. Therefore, we investigated the added value of transcriptomics for improving lineage assignment, currently based mainly on surface markers. First, we used an in-house pipeline to comprehensively detect genetic lesions in RNA-sequencing data (n=23) with a sensitivity >90%, as confirmed by targeted DNA sequencing. Second, we compared ALAL gene expression profiles (GEPs) with representative AML (n=145), B-ALL (n=223) and T-ALL (n=85) cases. In a principal component analysis, ALALs did not form a separate group, but clustered with AML, B-ALL or T-ALL. Accordingly, a machine learning classifier trained with GEPs of acute leukemias segregated 20/23 ALALs into myeloid-, B- or T-lymphoid leukemia. These 20 cases harbored genetic abnormalities consistent with the classifier-assigned leukemia. Furthermore, ALAL GEPs were deconvoluted with single-cell transcriptional profiles of normal hematopoietic cells using CIBERSORTx, revealing enrichment for signatures of lineages corresponding to the leukemic type predicted by our algorithm. The classifier was validated on an external ALAL cohort (n=24), assigning 75% of the patients to a lineage matching their immunophenotypic and methylation profiles. In conclusion, integrative analysis of RNA-sequencing data can accurately classify most ALAL cases. The pipeline and classifier developed in our study are valuable tools to improve ALAL diagnosis and guide therapeutic decisions.

  Study EGAS00001007967

• Dynamics of circulating tumor DNA in acute myeloid leukemia (AML) patients who undergo allogeneic transplantation

  Allogeneic hematopoietic cell transplantation (allo-HCT) is the only curative option for many patients with acute myeloid leukemia (AML). With the increasing use of non-myeloablative conditioning regimens, it is widely assumed that the cytotoxic conditioning itself does not eradicate all leukemia cells. Instead, the donor-derived immune cells are believed to eliminate residual leukemia cells over time. This so-called graft-versus-leukemia (GvL) is postulated to be the primary mechanism by which allo-HCT achieves sustained remissions in patients with AML. However, the precise quantification of GvL effects and its contribution to leukemia eradication has not yet been demonstrated by sequential, highly sensitive measurements of residual disease at various time points following allo-HCT. In the current study, we employed a novel, highly sensitive and specific plasma-based assay to track mutations in the leukemia cells of 14 patients before and after allo-HCT. If the primary therapeutic benefit of allo-HCT is truly mediated through a GvL effect, we hypothesized that leukemia cells would persist following cytotoxic conditioning and would often be reduced after post-transplant immunosuppression was discontinued. In the majority of patients (11 of 14, 79%), the conditioning used for allo-HCT did not eradicate the leukemia cells. In seven of these 11 patients, the leukemia cells persisted until immunosuppression was discontinued, at which time the persistent leukemia cells were reduced or eradicated. These results demonstrate the feasibility of using circulating tumor DNA (ctDNA) as a dynamic marker of GvL as well as minimal residual disease for patients with AML undergoing alloHCT, and have substantial implications for their management.

  Study EGAS00001007969

• Whole genome sequencing of tumor samples from advanced pre-treated NSCLC patients recruited in a phase I/II single-arm trial utilising CVA21, an oncolytic coxsackie virus, in combination with pembrolizumab.

  Acquired or de novo resistance to immune checkpoint inhibitors occurs in the majority of advanced non-small cell lung cancers (NSCLC). There is an unmet need to improve outcomes for patients with this condition. Oncolytic viruses represent an attractive treatment approach due to their dual activity in inducing tumor cell lysis directly and potentially augmenting anti-tumor immunity. Here, we present the safety, efficacy and translational findings from a phase I/II single-arm trial utilising CVA21, an oncolytic coxsackie virus, in combination with pembrolizumab in patients with advanced pre-treated NSCLC. We performed paired biopsies pre- and post-treatment in 10 patients who received intravenous CVA21 and pembrolizumab, 8 of whom had prior treatment with immune checkpoint inhibitor (ICI) therapy. Whole genome sequencing and spatial proteomics were performed to comprehensively characterise the response to CVA21. Combination CVA21/pembrolizumab (anti PD-1) therapy was well tolerated with no serious treatment-related adverse events. Partial responses were seen in two patients with prior acquired anti-PD-1 resistance and disease stabilisation in 6 patients, giving a clinical benefit rate of 80%. High baseline tumor mutational burden and PD-L1 expression were observed in patients with better response to treatment. Interestingly, an increase in antigen presentation and CD8+ T cell infiltration was observed on treatment compared with baseline in patients with better progression-free survival. This study demonstrates the potential of CVA21 to modulate the immunogenicity of tumor cells and remodelling the tumor microenvironment, providing novel insights for patient selection for trials involving novel immunotherapeutic approaches.

  Study EGAS00001008258

• European BestAgeing Study on microRNA candidates for cardiovascular disease

  Background: Circulating miRNAs have emerged as promising biomarker candidates due to their stability and their role in regulating key pathological pathways in cardiovascular disease (CVD). Yet, large-scale, multicenter studies examining their diagnostic and prognostic potential are scarce. This study evaluates the potential of miRNA expression profiles to inform disease classification and risk stratification across major CVD phenotypes, including acute coronary syndrome (ACS), chronic coronary artery disease (CAD), dilated cardiomyopathy (DCM), and ischemic cardiomyopathy (ICM), in a large, multicenter European cohort. Methods: We assessed genome-wide miRNA expression profiles in a total of 1,209 cardiovascular patients and 848 controls in a uniform, standardized fashion, which renders this study one of the largest prospective miRNA studies. To focus on only the most biologically plausible miRNAs for clinical translation, we mined all original studies of miRNA candidates in CVD and performed differential miRNA expression and enrichment analysis. We then trained disease-specific binary classification models to evaluate the diagnostic potential of miRNA signatures. Finally, we evaluated prognosis and disease severity based on distinct miRNA levels. Results: 634 original abstracts were identified, detailing 166 ACS, 181 CAD, 56 DCM, and 182 ICM miRNAs. Without further optimization, the signatures of a priori miRNAs already yielded very good diagnostic performance with ROC AUC of 0.83 – 0.95. There was an improvement when considering additional miRNAs in a discovery setting. Interestingly, in ACS, CAD and DCM we observed a significantly worse prognosis in probands with higher scores in the miRNA signatures, indicating additional prognostic information. Conclusions: The European BestAgeing miRNA study reveals emerging associations of several miRNA signatures with cardiovascular disease discrimination and prognostication, providing a foundation for future external validation and potential clinical translation of this class of markers.

  Study EGAS00001008346

• Single cell RNA sequencing of synovial B cells in early Rheumatoid Arthritis

  This dataset includes 3432 paired single cell sequencing fastq files derived from synovial B cells of 5 early Rheumatoid Arthritis patients. Libraries were prepared using the SmartSeq2 protocol. All wells contained ERCC spike-ins. Libraries were sequenced on an Illumina NovaSeq instrument with 2x100bp paired-end reads yielding a median of 2.5M reads/well. Sample aliases ending on 368 and 384 represent empty control wells. Sample aliases starting with P11417_4 and P13157_1 belong to ACPA- patient A7, sample aliases starting with P11417_5 and P11417_6 belong to ACPA+ patient A1, sample aliases starting with P13157_2 belong to ACPA- patient A3, sample aliases starting with P13157_3 and P13157_6 belong to ACPA+ patient A2, sample aliases starting with P13157_4 and P13157_5 belong to ACPA- patient A4.

  Dataset EGAD00001009076

• WGS of cell-free DNA derived from plasma of patients with pediatric sarcoma and healthy controls, and lcWGS/RRBS of matched tumor tissue

  This dataset contains: i) 241 deep (median 12x) whole-genome sequencing profiles of 95 patients with Ewing sarcoma, 31 patients with other pediatric sarcomas, and 22 additional profiles from healthy controls. Sequencing was performed on a NovaSeq 6000 instrument using the NovaSeq S4 2x100 bp configuration. In addition, pilot experiments for 18 cfDNA samples were performed using Illumina HiSeq 2000/2500 machines (2x75 bp configuration). Data is provided as .fastq.gz files (2 files, .R1 and .R2, per sample). i) Low coverage whole-genome-sequencing on 43 tumor biopsy samples from patients with Ewing sarcoma with matched cfDNA samples. The samples were sequenced using a NovaSeq 6000 instrument with the NovaSeq S4 1x100 bp configuration. Data are provided as unmapped (raw) .bam files. iii) Reduced-representation bisulfite sequencing data for 38 tumor biopsy samples from patients with Ewing sarcoma with matched cfDNA samples, and 2 control samples. RRBS libraries were sequenced on Illumina HiSeq 2000/2500 machines. Data are provided as unmapped (raw) .bam files.

  Dataset EGAD00001007080

• Next Generation Children project

  Whole genome sequencing of sick children in neonatal and paediatric intensive care units, aligned to reference assembly GRCh38.

  Dataset EGAD00001007868

• panALL exome sequencing, data set2, 700 samples

  panALL exome sequencing, data set2, 700 samples

  Dataset EGAD00001008268

• EPICC: Evolutionary Predictions in Colorectal Cancer WGS-3

  WGS profiling bam files from colorectal carcinoma and adenoma.

  Dataset EGAD00001008642

• Clonal dynamics after allogeneic haematopoietic cell transplantation using genome-wide somatic mutations - WGS

  Allogeneic haematopoietic cell transplantation (HCT) replaces the stem cells responsible for blood production with those harvested from a donor, and is received by 40,000 patients worldwide each year. To quantify dynamics of long-term stem cell engraftment, we sequenced whole genomes of 2,824 single-cell-derived haematopoietic colonies from blood samples of 10 donor-recipient pairs taken 9-31 years after HLA-matched sibling HCT. With younger donors, 10,000-50,000 stem cells had engrafted and were still contributing to haematopoiesis at time of sampling, but estimates were 10-fold lower with older donors. Engrafted stem cells made multilineage contributions to myeloid, B-lymphoid and T-lymphoid populations, although individual clones often showed biases towards one or other mature cell type. Recipients had lower clonal diversity than matched donors, equivalent to ~10-15 years of additional ageing, arising from up to 25-fold greater expansion of stem cell clones. An HCT-related population bottleneck alone could not explain these differences: instead, phylogenetic trees evinced two distinct modes of HCT-specific selection. In 'pruning selection', cell divisions underpinning recipient-enriched clonal expansions had occurred in the donor, preceding transplant - their selective advantage derived from preferential mobilisation, harvest, survival ex vivo or initial homing. In 'growth selection', cell divisions underpinning clonal expansion occurred through proliferative advantage in the recipient's marrow after homing - clones with multiple driver mutations especially demonstrated this pattern. Uprooting stem cells from their native environment and transplanting them to foreign soil exaggerates selective pressures, distorting and accelerating the loss of clonal diversity compared to the unperturbed haematopoiesis of donors.

  Dataset EGAD00001010872

• A study of the genetic basis of evation by Acute Myeloid Leukaemia of Graft vs Leukaemia effects after allogeneic bone marrow transplantation

  Acute myeloid leukaemia (AML) is an aggressive and molecularly diverse disease with a poor overall survival of 20-25%. With an annual incidence of 2.9 per 100,000, AML is currently the commonest myeloid malignancy in Europe, yet the two main therapeutic options for this disease, anthracyclines and purine analogues, have remained unchanged for over 20 years. Currently patients are stratified at diagnosis according to a series of clinicopathological parameters (e.g. age, white cell count and presence/absence of previous clonal haematological disease) and molecular markers (e.g. chromosomal translocations/deletions, aneuploidy and mutations in genes such as FLT3 and NPM1). Patients with adverse prognostic features, whose prognosis is particularly poor (e.g. <15% long-term survival) are offered treatment with allogeneic bone marrow transplantation (allo-BMT) if a sibling or unrelated donor is available. This can significantly improve survival (e.g. up to 40% long-term survival in some contexts), albeit at the expense of significant toxicity and transplant-related mortality (TRM). Allo-BMT is thought to work in part by allowing the delivery of large doses of chemotherapy followed by haemopoietic "rescue" with donor haemopoietic stem cells (haemopoietic failure would otherwise ensue). However, potentially the most potent effect of allo-BMT is the cytotoxic effect of donor lymphocytes against AML blasts, a phenomenon known as graft-vs-leukaemia (GVL) effect. Increasingly, transplants using reduced chemotherapy intensity (mini-allografts) are being used that partially circumvent the toxicity from chemotherapy and rely on GVL to effect cure. Nevertheless, AML relapse after allo-BMT still occurs at a significant rate of up to 80% depending on the type of transplant. There is accumulating evidence that genetic events in residual leukaemic cells enable them to evade immunodetection and therefore survive the GVL effect and expand to cause relapse. The most striking example of this is the loss of HLA antigens after transplants in which donor and recipient are not fully HLA-matched. In these cases, the leukaemia "deletes" the genomic region containing the disparate HLA antigen which was preferentially targeted as "foreign" by the GVL effect. However, the genetic basis of immune evasion in the majority of transplants, which are fully HLA matched, is not known. One possibility is that loss of genes coding for antigens outside the HLA locus but which are also targets of GVL may operate, alternatively genetic events that affect processes downstream of immunological cytotoxicity may be responsible. The identification of genetic events that mediate immune evasion would not only facilitate the understanding of this process but can help plan therapeutic interventions that improve the outcomes of allogeneic transplantation for AML and other disorders. We intend to study this by conducting exome sequencing on 6 cases of AMLs from patients that attend my clinic at Addenbrooke's hospital and have relapsed after allogeneic transplantation. Samples from AML diagnosis, remission/normal and AML relapse (total n=18) will be studied to identify somatic mutations in the primary AML and those acquired by the relapsed clone. The 18 samples will also be studied by array CGH to detect regions of genomic amplification or deletion.

  Dataset EGAD00001000404

• Targeted sequencing of genes recurrently mutated in AML

  Background Massively parallel sequencing technology has transformed cancer genomics. It is now feasible, in a clinically relevant time-frame, for a clinically manageable cost, to screen DNA from patient tumours for mutations essentially genome-wide. The challenge for personalised medicine will be to increase the sample size to thousands or tens of thousands of well-characterised cases in order to attain sufficient statistical power to stratify patients accurately across the complexity and genomic heterogeneity expected for most of the common tumour types. Currently, whole genome sequencing on this scale is not feasible, and targeted sequencing of relevant portions of the genome will be required. Pilot data We have developed protocols for large-scale, multiplexed sequencing of 100-200 genes in thousands of samples. Essentially, using robotic technology, genomic DNA from the cancer specimen is processed into sequencing libraries with unique DNA barcodes, thereby allowing sequencing reads to be attributed to the sample they derive from. Currently, these sequencing libraries can be generated in a 96-well format using fully automated protocols, and we are exploring methods to expand this to a 384-well format. The sequencing libraries are pooled and hybridized to custom sets of RNA baits representing the genomic regions of interest. Sequencing of the pulled-down libraries is done in pools of 48-96 samples per lane of an Illumina Hi-Seq. This protocol is already implemented at the Sanger Institute. We have published proof that somatic mutations in novel cancer genes can be identified from exome-wide sequencing. In unpublished pilot data, we have established the feasibility of robotic library production, custom pull-down, and multiplexed sequencing of barcoded libraries for 100 known myeloid cancer genes across 760 myelodysplasia samples. Highlights of the data thus far analysed reveal that the coverage is remarkably even between samples; when 96 samples are run, average coverage per lane of sequencing is ~250, with 90-95% of targeted exons covered by >25 reads; known mutations can be discovered in the data set; and the protocol is amenable to whole genome amplified DNA. The bioinformatic algorithms for identification of substitutions and indels in pull-down data are well-established; we have pilot data proving that copy number changes, LOH and genomic rearrangements in specific regions of interest can also be identified by tiling of baits across the relevant loci. Proposal We propose to apply this methodology to 10000 samples from patients with AML enrolled in clinical trials over the last 10-20 years. Oncogenic point mutations and potentially genomic rearrangements will be identified, and linked to clinical outcome data, with a view to undertaking the following sorts of analyses: ? Identification of co-occurrence, mutual exclusivity and clusters of driver mutations. ? Correlation of prognosis with driver mutations and potentially gene-gene interactions ? Exploration of genomic markers of drug response Ultimately, we would like to be in a position to release the mutation data together with matched clinical outcome data to genuine medical researchers via a controlled access approach, possibly within the COSMIC framework (www.sanger.ac.uk/genetics/CGP/cosmic/). The vision here is to generate a portal whereby a clinician faced with an AML patient and his / her mutational profile can obtain a ?personalised? prediction of outcome, together with a fair assessment of the uncertainty of the estimate. With a sufficient sample size, there would also be the potential to develop decision support algorithms for therapeutic choices based on such data.

  Dataset EGAD00001000606

• Multi-omic analysis of Down Syndrome in thyroid

  Background: Down syndrome (DS) is a chromosomal disorder characterized by trisomy of chromosome 21 and is associated with various health problems. Congenital thyroid dysfunction is a common endocrine abnormality in DS, but its exact cause remains unknown. This work aimed at studying the etiology of congenital thyroid dysfunction in DS. Methods: Thyroid tissue of fetuses with DS (n=4) and fetuses without a genetic/developmental abnormality (n=5) were analyzed using histopathological evaluation, RNA sequencing (RNA-seq), and DNA methylation (DNAm) profiling. Results: Histological examination showed abnormally developed DS thyroid tissue compared to non-DS/healthy tissue. RNA-seq analysis showed numerous differences in gene expression between DS and healthy thyroid tissue, including significant downregulation of thyroid genes FOXE1, IYD and DIO2. Gene set enrichment analysis showed disruption of vital cellular processes and functions in DS fetal thyroid tissue. Analysis of DNAm and expression quantitative trait methylation (eQTM) analysis identified functional and biologically relevant associations, that can be linked to plausible disrupted pathways in DS. Conclusions: Based on these observations, we conclude that congenital non-autoimmune thyroid dysfunction in DS is probably a DS-specific form of thyroid dysfunction, characterized by abnormal thyroid development, altered expression of genes that regulate thyroid development and hormone production, and altered expression and DNAm of genes pertaining to numerous vital pathways leading to significant metabolic disturbances. The disturbances as investigated by multi-omics analyses (RNA-seq, DNAm and eQTM analyses) are widely spread over the genome, and presumably caused by the direct and downstream effects of overexpressed chromosome 21 genes (dosage effect), and (epi)genomic disturbances.

  Study EGAS00001007677

• Identification of intronic sequences promoting natural skipping of ABCA4 exon 15 using long-read transcriptomics and midigene assays

  Stargardt disease is an inherited retinal disease that leads to progressive loss of vision, and is caused by biallelic variants in the ABCA4 gene. A significant proportion of ABCA4 variants of variable penetrance affect pre-mRNA splicing. Previously, we observed the exclusion of certain ABCA4 exons into the final transcript in the absence of any disease-causing variants. In this study, we explored this phenomenon that we coined natural exon skipping (NES) for the entire ABCA4 transcript, using targeted RT-PCR analysis as well as long-read transcriptome sequencing of control retina samples. For the most abundant NES event, which is skipping of ABCA4 exon 15, overlapping ABCA4 deletion constructs were employed to identify responsible cis-regulatory sequence motifs. Subsequently, antisense oligonucleotides (AONs) were designed and tested to reverse NES and restore fully-spliced ABCA4 mRNA. We identified intronic regions on each side of exon 15 that, when deleted, increased the inclusion of exon 15 in the final transcripts. Interestingly, these motifs appear to work in a synergistic manner. Further fine-mapping of these regions revealed several cis-regulatory motifs that may be responsible for inducing NES, on each side of the exon. AONs targeting these motifs unfortunately did not decrease the skipping of exon 15. In conclusion, we discovered an intriguing phenomenon that, in the absence of disease-causing variants, results in the skipping of one or more exons in the ABCA4 transcript. These NES events appear to be regulated by cis-regulatory motifs within the flanking introns. This work captured the extend of NES throughout ABCA4 transcript and the role of cis-regulatory motifs in pre-mRNA splicing.. The use of AONs targeting cis-regulatory motifs is an innovative strategy for splicing modulation, suitable for complementing future variant-specific splicing modulation therapies in ABCA4 retinopathies.

  Study EGAS50000000071

• Blueprint Data Access Committee.

  Access to the sequence and alignment and genotype level data produced by the Blueprint Consortium is controlled by the Blueprint Data Access Committee (DAC). Access to data will be granted to qualified investigators for appropriate use. Please email blueprint-dac@ebi.ac.uk to request the Data Access Application and Data Access Agreement.

  Dac EGAC00001000135

• Integrating molecular imaging and transcriptomic profiling in advanced HER2-positive breast cancer receiving trastuzumab emtansine (T-DM1): an analysis of the ZEPHIR clinical trial

  Purpose: The ZEPHIR clinical trial evaluated the role of [89Zr]trastuzumab-PET/CT (HER2-PET/CT) and 2-[18F]fluoro-2-deoxy-D-glucose PET/CT ([18F]FDG-PET/CT) in predicting outcomes in patients with advanced HER2-positive breast cancer treated with trastuzumab emtansine (T-DM1). Here, we combined molecular/metabolic imaging and transcriptomic data to investigate the biological processes associated with [89Zr]trastuzumab and [18F]FDG uptake, and to dissect the mechanisms involved in T-DM1 resistance. Experimental design: RNA was extracted from metastasis biopsies obtained at study inclusion in the ZEPHIR trial. HER2-PET/CT and [18F]FDG-PET/CT imaging data of biopsied lesions were integrated with RNA sequencing data. Lesions were compared based on the level of [89Zr]trastuzumab uptake (dichotomized into high/low classes) as well as on the presence/absence of metabolic response, defined comparing baseline and on-treatment [18F]FDG-PET/CT. Results: We analyzed matched transcriptomic and molecular/metabolic imaging data for 24 metastases (one lesion/patient). Pathways involved in extracellular matrix (ECM) organization, glycosylphosphatidylinositol synthesis, and collagen formation were enriched in lesions presenting low [89Zr]trastuzumab uptake. [18F]FDG uptake at baseline correlated with proliferation and immune-related processes. Hypoxia and ECM-related processes were enriched in lesions that showed no metabolic response to T-DM1. Interestingly, immune-related processes (e.g., Fc gamma receptor activation, complement cascade) were associated with high [89Zr]trastuzumab uptake and metabolic response. Conclusions: To our knowledge, this study represents the first correlative analysis between [89Zr]trastuzumab tumor uptake and gene expression profiling in humans. Our findings point toward a role of ECM in impairing [89Zr]trastuzumab tumor uptake and T-DM1 metabolic response in advanced HER2-positive breast cancer, and highlight the potential of molecular imaging to depict tumor microenvironment features.

  Study EGAS50000000470

• Genomic Profiling of Pediatric Tumors by Cell Free DNA Sequencing

  BACKGROUND Analysis of cell-free DNA (cfDNA) from peripheral blood has been studied extensively for detection and monitoring therapy response in adult cancers, but utility in pediatric tumors is unclear. METHODS We studied 233 children with hematologic, solid and brain tumors. cfDNA was extracted from 1 mL of plasma collected at diagnosis. Circulating tumor DNA (ctDNA) in cfDNA was quantified and diverse classes of somatic genomic variants were detected using PeCan-Seq. PeCan-Seq findings were compared with matched tumor whole-genome and whole-exome sequencing data. Minimal residual disease (MRD) status was evaluated by cfDNA PeCan-Seq in 14 patients with B-cell acute lymphoblastic leukemia. RESULTS Plasma cfDNA yields at diagnosis were higher for hematologic malignancies than for solid and brain tumors. Patients with hematologic malignancy were 100% ctDNA-positive, the median ctDNA fraction in cfDNA was 0.77. We detected 97% variants that were present in tumors and covered by the capture panel, and 50 subclonal variants that were present in cfDNA only. Our approach enabled detection of sequence and structural variants in B-ALL ctDNA at varying levels of MRD. Nineteen of 38 patients with solid tumors were ctDNA-positive, with high prevalence in Ewing Sarcoma, liver tumors and osteosarcoma. One of 19 brain tumor patients was ct-DNA positive, indicating lack of utility for plasma-derived cfDNA in this disease. Notably, cfDNA analysis of a single plasma sample identified genomic features in ctDNA derived from B-ALL, concurrent primary neuroblastoma, and a later neuroblastoma relapse, indicating that genomic variations from multiple tumors could be detected simultaneously. CONCLUSIONS Pediatric patients with hematologic malignancy had an abundance of plasma ctDNA at diagnosis. Plasma ctDNA levels varied between different solid tumor types, more studies are needed to determine utility. Panel-based cfDNA sequencing could serve as a non-invasive approach for pediatric cancer diagnosis and disease monitoring.

  Study EGAS50000000072

• Integrative analysis reveals a macrophage-predominant, immunosuppressive immune microenvironment and subtype-specific therapeutic vulnerabilities in advanced salivary gland cancer

  Background Salivary gland cancers (SGC) are rare and heterogeneous malignant tumors. Advanced SGC lack established treatment options and show poor response to immunotherapy. Here, an integrative multi-omics analysis in a large cohort of advanced SGC revealed insights into the tumor immune microenvironment (TIM) and distinct mechanisms of immune evasion. Methods A total of 104 patients with recurrent/metastatic SGC from the DKTK MASTER program were included in this study. Whole-exome or whole-genome sequencing and RNA-sequencing was performed on fresh frozen tumor tissue. The tumor immune microenvironment was analyzed using CIBERSORT deconvolution analysis and immune gene expression scores in bulk RNA-sequencing data. Single-nuclei sequencing and immunohistochemistry analyses were performed in selected samples. Results were validated in bulk RNA-sequencing data of a previously published independent dataset. Results Bulk transcriptome analysis revealed an immune-deserted TIM in the majority of advanced SGC samples. Immune exclusion was most prominent in adenoid cystic carcinoma (ACC) subgroup 1 exhibiting a downregulation of the antigen processing machinery. Only a small subset of advanced SGC, including few adenoid cystic carcinoma, exhibited T-cell inflammation, which was correlated with tumor mutational burden in Non-ACC samples. Subtype specific expression of immune checkpoints as well as cancer testis antigens were identified with prominent expression of VTCN1 in luminal cells within ACC. Single-cell RNA-sequencing and bulk RNA-seq deconvolution analysis validated immune cell exclusion and revealed a TIM that was dominated by M2 macrophages across SGC subtypes. Among evaluable patients treated with immune checkpoint inhibitors, a high T-cell to macrophage ratio was associated with clinical benefit. Conclusions These data support biomarker-based development of immune-checkpoint inhibition and the development of novel immune-checkpoint inhibitors and cellular therapies in SGC.

  Study EGAS50000000809

• A prospective study identifies MisMatch Repair genes as candidate predisposing genes for uveal melanoma.

  Objective: The aim of this study is to identify new predisposing genes for uveal melanoma (UM), the most common intraocular malignancy in adults. Design: Prospective cohort study from Institut Curie. Participants: 381 unselected patients diagnosed with UM between July 2021 and February 2023 at Institut Curie. Methods: All patients received genetic counseling and consented to extended genetic testing. A panel of 122 genes was analyzed by targeted sequencing on germline DNA from these patients. Main outcomes measures: UM was diagnosed by ophthalmologists at Institut Curie, and a senior pathologist confirmed the diagnosis when tumor or biopsy was available. Classification of the germline genetic variants was performed according to the ACMG guidelines. Results: We identified at least one pathogenic variant (PV) in 79 patients. Among those variants, 20 were found in clinically relevant genes for the patient’s follow-up and genetic counseling. We found a significant enrichment of germline pathogenic variants in Mismatch Repair (MMR) MLH1, MSH2 and MSH6 genes in this UM cohort, suggesting that MMR germline alterations could predispose to UM. From this cohort, one tumor was available from a patient carrying a MLH1 germline mutation. The tumor was characterized by a monosomy 3, in line with a loss of heterozygosity (LOH) of MLH1, located on chromosome 3. Loss of expression of MLH1 and PMS2 was observed by immunohistochemistry, and MMR mutational signatures SBS6, ID1 and ID2 were identified from the whole genome sequencing of the tumor. These data suggest that MLH1 is implicated in the oncogenesis of this UM. Conclusions: This study is the largest prospective constitutional study performed on UM patients, and the first demonstrating that MMR germline alterations are significantly enriched among these patients. Further analyses will be needed to estimate the risk of developing UM for patients with Lynch syndrome

  Study EGAS50000000914

• Emirati Genome Project subset of 43,608 WGS samples for population-scale variant discovery and allele frequency mapping.

  Here, we present a comprehensive genomic characterization of a cohort of 43,608 Emirati genomes sequenced as part of the Emirati Genome Program (EGP). This study identified more than 421 million single nucleotide variants and indels and more than 600 million copy number- and structural variants. Small variants had 756 million molecular effects annotated. Of the identified variants, 69.04% were observed with more than one allele. Of 7.7 million polymorphic variants having an allele frequency (AF) of more than 5% in EGP, 1,348 have a predicted deleterious effect on a protein. Characterization of global variation reveals that EGP represents a genetic continuum encompassing the range of African, Asian, and European populations. It is best described by two Arabian, an Eurasian, and an African component, with the predominant Arabian component linked to mitochondrial haplogroups J and T, which are commonly attributed to the Middle East. Various aneuploidies of sex chromosomes were detected in 93 individuals overall and of chromosome 21 in 41 individuals. The median inbreeding coefficient and cumulative runs of homozygosity (ROHs) lengths were increased due to extensive consanguinity, being largest in the groups with Arabian main ancestry components and higher than reported for Qatar. Families were identified based on genetic relatedness and classified into 264 families with unrelated parents and 247 families with third- and fourth-degree consanguineous parents. Representative consanguineous pedigrees of families in EGP were outlined. Cumulative ROHs were affected by the main ancestry component and significantly increased in offspring of consanguineous parents, with a pronounced difference between 3rd and 4th-degree relatedness. Investigation of cumulative AFs of variants causing Mendelian diseases highlighted genes related to alpha- and beta-thalassemia (HBB, HBA2). It showed a high burden of variants causing severe recessive diseases, metabolic and retinal disorders, and hearing loss. In summary, EGP represents a landmark effort in characterizing the genetic diversity of the Emirati population, leveraging the largest Middle Eastern cohort reported to date.

  Study EGAS50000001071

• Single-Gene vs. Panel Sequencing in Advanced HR+/HER2− Breast Cancer

  Background: Hormone receptor-positive (HR+)/HER2-negative (HER2−) breast cancer is the most common subtype of breast cancer, with biomarker-driven therapies improving outcomes. Liquid biopsies, which analyze circulating tumor DNA (ctDNA), are increasingly essential in clinical management, offering a minimally invasive approach to assess somatic alterations and guide therapy decisions. However, the choice of assay directly affects clinical utility and global access to validated testing platforms remains inconsistent. To address these challenges, this study aimed to compare the performance of single-gene and panel-based sequencing for detecting PIK3CA mutations in HR+/HER2− breast cancer and evaluate the added value of broader genomic profiling in guiding treatment decisions. Patients and methods: We conducted a prospective, multicenter study comparing our previously established SiMSen-Seq (SSS) assay, a single-gene method targeting PIK3CA hotspot mutations, with the AVENIO ctDNA Expanded assay, enriching for 77 clinically relevant cancer genes. Additionally, tumor fractions were estimated with mFAST-SeqS by assessing chromosomal aneuploidy. We analyzed 161 plasma samples collected from 146 patients with HR+/HER2− advanced breast cancer before starting a new line of palliative therapy. Results: High concordance (92.6%) was observed between SSS and AVENIO for PIK3CA mutation detection, with accurate variant allele frequency (VAF) quantification confirmed by Bland-Altman analysis and linear mixed-effects modeling. AVENIO identified additional PIK3CA mutations in 19% of samples as well as other actionable alterations, including ESR1 mutations (17.5%) and PI3K pathway alterations (40.6%). Panel sequencing improved the interpretation of negative liquid biopsy results by integrating tumor fraction and VAF metrics. Conclusions: Panel sequencing assays such as AVENIO provide broader genomic coverage, enhance the detection of actionable mutations and refine the interpretation of negative liquid biopsy findings. These findings support the integration of comprehensive ctDNA profiling into routine clinical practice for HR+/HER2− advanced breast cancer while highlighting the need for improved assay accessibility and expanded genomic targets.

  Study EGAS00001008200

• CDK4/6 inhibition in advanced chordoma: final results of the NCT PMO-1601

  Background. Chordoma is a rare bone tumor with no approved systemic therapy. Inactivation of the tumor suppressor p16 (encoded by CDKN2A) leads to aberrant cyclin D-CDK4/6-RB pathway activity, which can be inhibited by palbociclib in chordoma cell lines. Patients and Methods. We conducted a phase II single-arm, open-labeled trial on palbociclib in adult patients with advanced chordomas with p16 loss (by immunohistochemistry) or CDKN2A loss (by genomic analysis) and CDK4/6 and RB1 presence (by immunohistochemistry or RNA-Sequencing). Patients received 125mg palbociclib once daily for 21 days in a 28-day cycle. The study used a Simon optimal 2-stage design. The primary endpoint was disease control rate after six cycles. Response was evaluated according to Response Evaluation Criteria in Solid Tumors (RECIST) version 1.1. Secondary endpoints included progression-free survival, overall survival, and biomarker analysis. The study was considered positive if 25% of patients reached the primary endpoint. Results. Between December 2017 and April 2023, 28 patients (median age 60 years, range 31-84) were enrolled. After a median follow-up of 28 months, the disease control rate was 39%, with 11 out of 28 patients achieving stable disease. No objective responses were observed. The median progression-free survival was 5.6 months, and the median overall survival was 24.6 months. Treatment was well tolerated without new safety signals. CDKN2A deletions were the most common genomic alterations underlying p16 loss and were accompanied by MTAP co-deletion. There was no correlation between responder phenotypes by immunohistochemistry and outcome. RNA-based analysis showed heterogeneous pathways deregulation at study entry suggesting underlying primary resistance mechanisms Conclusion. This trial met its primary endpoint with a 39% disease control rate, though antitumor activity was modest. Molecularly tailored combination therapies should be considered to improve efficacy. (ClinicalTrials.gov number, NCT03110744).

  Study EGAS00001007985

• A circulating biomarker for severity of facioscapulohumeral muscular dystrophy

  Facioscapulohumeral muscular dystrophy (FSHD) is a prevalent, incurable skeletal myopathy. Clinical trials for FSHD are hindered by heterogeneous biomarkers poorly associated with clinical severity, requiring invasive muscle biopsy. Macroscopically FSHD presents with slow fatty replacement of muscle, rapidly accelerated by inflammation. Mis-expression of the transcription factor DUX4 is currently accepted to underlie FSHD pathogenesis and mechanisms including PAX7 target gene repression have been proposed. Here we perform RNA-sequencing on MRI guided inflamed (n=24) and isogenic non-inflamed (n=24) muscle biopsies from the same clinically characterised FSHD patients, alongside isogenic peripheral blood mononucleated cells from a subset of patients (n=13), and unaffected controls. We employ multivariate models to evaluate the clinical associations of 5 published FSHD transcriptomic biomarkers. We demonstrate that PAX7 target gene repression can discriminate control, inflamed and non-inflamed FSHD muscle independently of age and sex (p<0.013), while the discriminatory power of DUX4 target genes is limited to distinguishing FSHD muscle from control. Importantly, the level of PAX7 target gene repression in non-inflamed muscle associates with clinical assessments of FSHD severity (p=0.04). DUX4 target gene biomarkers in FSHD muscle show associations with lower limb fat fraction and D4Z4 array length, but not clinical assessment. Lastly, PAX7 target gene repression in FSHD muscle correlates with the level in isogenic peripheral blood mononucleated cells (p=0.002). A refined PAX7 target gene biomarker comprising 143/601 PAX7 target genes computed in peripheral blood (the FSHD muscle-blood biomarker) associates with clinical severity in FSHD patients (p<0.036). Our new circulating biomarker validates as a classifier of clinical severity in an independent data set of 54 FSHD patient blood samples, with improved power in older patients (p=0.03). In summary, we present the minimally invasive FSHD muscle-blood biomarker of FSHD clinical severity valid in patient muscle and blood, with potential utility in routine disease monitoring and clinical trials.

  Study EGAS00001007350

• First genome-wide association study in an Australian Aboriginal population provides insights into genetic risk factors for Body Mass Index and Type 2 Diabetes

  A body mass index (BMI) >22kg/m2 is a risk factor for type 2 diabetes (T2D) in Aboriginal Australians. To identify loci associated with BMI and T2D we undertook a genome-wide association study using 1,355,133 quality-controlled single nucleotide polymorphisms (SNPs) genotyped (Illumina 2.5M Duo Beadchip) in 402 individuals in extended pedigrees from a Western Australian Aboriginal community. Imputation using the thousand genomes (1000G) reference panel extended the analysis to 6,724,284 post quality-control autosomal SNPs. No associations achieved genome-wide significance, commonly accepted as P<5x10-8. Nevertheless, genes/pathways in common with other ethnicities were identified despite the arrival of Aboriginal people in Australia >45,000 years ago. The top hit (rs10868204 Pgenotyped=1.50x10-6; rs11140653 Pimputed_1000G=2.90x10-7) for BMI lies 5’ of NTRK2, the type 2 neurotrophic tyrosine kinase receptor for brain-derived neurotrophic factor (BDNF) that regulates energy balance downstream of melanocortin-4 receptor (MC4R). PIK3C2G (rs12816270 Pgenotyped=8.06x10-6; rs10841048 Pimputed_1000G=6.28x10-7) was associated with BMI, but not with T2D as reported elsewhere. BMI also associated with CNTNAP2 (rs6960319 Pgenotyped=4.65x10-5; rs13225016 Pimputed_1000G=6.57x10-5), previously identified as the strongest gene-by-environment interaction for BMI in African-Americans. The top hit (rs11240074 Pgenotyped=5.59x10-6, Pimputed_1000G=5.73x10-6) for T2D lies 5’ of BCL9 that, along with TCFL2, promotes beta-catenin’s transcriptional activity in the WNT signaling pathway. Additional hits occurred in genes affecting pancreatic (KCNJ6, KCNA1) and/or GABA (GABRR1, KCNA1) functions. Notable associations observed for genes previously identified at genome-wide significance in other populations included MC4R (Pgenotyped=4.49x10-4) for BMI and IGF2BP2 Pimputed_1000G=2.55x10-6) for T2D. Our results may provide novel functional leads in understanding disease pathogenesis in this Australian Aboriginal population.

  Study EGAS00001001004

• Nala GSI GSAv3 PGx Study

  Background: Microarrays are a well-established and widely adopted technology capable of interrogating hundreds of thousands of loci across the human genome. Combined with imputation to cover common variants not included in the chip design, they offer a cost-effective solution for large-scale genetic studies. Beyond research applications, this technology can be applied for testing pharmacogenomics, nutrigenetics, and complex disease risk prediction. However, establishing clinical reporting workflows requires a thorough evaluation of the assay's performance, which is achieved through validation studies. In this study, we performed pre-clinical validation of a genetic testing workflow based on the Illumina Global Screening Array for 25 pharmacogenomic-related genes. Methods: To evaluate the accuracy of our workflow, we conducted multiple pre-clinical validation studies. Here, we present the results of accuracy and precision assessments, involving a total of 73 cell lines. These assessments encompass reference materials from the Genome-In-A-Bottle (GIAB), the Genetic Testing Reference Material Coordination Program (GeT-RM) projects, as well as additional samples from the 1000 Genomes project (1KGP). We conducted an accuracy assessment of genotype calls for target loci in each indication against established truth sets. Results: In our per-sample analysis, we observed a mean analytical sensitivity of 99.39% and specificity 99.98%. We further assessed the accuracy of star-allele calls by relying on established diplotypes in the GeT-RM catalogue or calls made based on 1KGP genotyping. On average, we detected a diplotype concordance rate of 96.47% across 14 pharmacogenomic-related genes with star allele-calls. Lastly, we evaluated the reproducibility of our findings across replicates and observed 99.48% diplotype and 100 % phenotype inter-run concordance. Conclusion: Our comprehensive validation study demonstrates the robustness and reliability of the developed workflow, supporting its readiness for further development for applied testing.

  Study EGAS00001007710

• mRNA capture sequencing and RT-qPCR for the detection of pathognomonic, novel and secondary fusion transcripts in formalin-fixed paraffin-embedded tissue: a sarcoma showcase

  RNA sequencing is widely used to study gene abundance patterns, but less frequently to assess structural alterations such as fusion transcripts. In this study, we assess the performance of messenger RNA (mRNA) capture sequencing to identify fusion transcripts in formalin-fixed paraffin-embedded (FFPE) tissue of different types of sarcoma, followed by orthogonal confirmation using RT-qPCR. To validate our workflow, we used the TruSight RNA Pan-Cancer Panel to analyze a first cohort of 6 positive control tumors for which the diagnostic workup demonstrated the presence of a specific chromosomal rearrangement. Fusion transcript calling by FusionCatcher confirmed all known aberrations and enabled the identification of both fusion gene partners and breakpoints. Next, whole transcriptome TruSeq RNA Exome sequencing was applied to a second cohort of 17 patients, diagnosed with fusion gene-negative alveolar rhabdomyosarcoma (ARMS) or undifferentiated round cell sarcoma (URCS), for whom fluorescence in situ hybridization (FISH) did not identify the classical pathognomonic rearrangements. Remarkably, for 6 out of 17 patients, mRNA capture sequencing readily detected a pathognomonic fusion transcript, i.e. PAX3-FOXO1 in 2 ARMS patients, and EWSR1-FLI1, EWSR1-ERG or EWSR1-NFATC2 in 4 URCS patients. These data show that RNA capture sequencing may enhance the detection rate of pathognomonic fusion genes in sarcoma. For the 11 patients with no known pathognomonic fusion transcript, 11 newly identified fusion transcripts were confirmed by RT-qPCR, including COPS3-TOM1L2, NCOA1-DTNB, WWTR1-LINC01986, PLAA-MOB3B, AP1B1-CHEK2 and BRD4-LEUTX fusion transcripts in ARMS patients. Additionally, recurrently detected secondary fusion transcripts in two patients diagnosed with EWSR1-NFATC2-positive sarcoma were confirmed (COPS4-TBC1D9, PICALM-SYTL2, SMG6-VPS53 and UBE2F-ALS2). Our results provide new insights into the underlying genetic causes of these malignancies.

  Study EGAS00001005202

• In this study, blood-brain barrier (BBB)-forming brain endothelial-like cells were generated from apolipoprotein E gene allele E4 (APOE4, high AD risk) and allele E3 (APOE3, lower AD risk) carrying patient-derived induced pluripotent stem cells (iPSCs). Cells were subsequently exposed to focused ultrasound and microbubbles (FUS+MB) to induce BBB opening and their transcriptome analysed. RNA sequencing (RNA-seq) results demonstrated minimal changes in the gene expression following FUS+MB suggesting safety of FUS+MB application in the clinical setting.

  The blood-brain barrier (BBB) is a major impediment to therapeutic intracranial drug delivery for the treatment of neurodegenerative diseases, including Alzheimer’s disease (AD). Focused ultrasound applied together with microbubbles (FUS+MB) is a novel technique to transiently open the BBB and increase drug delivery. Evidence suggests that FUS+MB is safe, however the effects of FUS+MB on human BBB cells, especially in relation to AD, remain sparsely investigated. Here we developed a human sporadic AD BBB cell platform to investigate the effects of FUS+MB on BBB cells and screen for the delivery of two potentially therapeutic AD antibodies. We generated BBB cells from apolipoprotein E gene allele E4 (APOE4, high AD risk) and allele E3 (APOE3, lower AD risk) carrying patient-derived induced pluripotent stem cells (iPSCs). We utilized this robust and reproducible platform to demonstrate increased delivery of potentially therapeutic antibodies across the BBB, including an analogue of Aducanumab (AduhelmTM) following FUS+MB treatment. Our results also demonstrate the safety of FUS+MB indicated by minimal changes in the cell transcriptome as well as little or no changes in cell viability and inflammatory responses within first 24 h post FUS+MB, with potentially important implications in a clinical setting. Finally, we report a more physiologically relevant hydrogel-based 2.5D BBB model as a key development for FUS+MB-mediated drug delivery screening, with potentially higher translational utility. Our results provide an important advancement for identifying FUS+MB-deliverable drugs and screening for cell- and patient-specific effects of FUS+MB, accelerating its use as a therapy in AD.

  Study EGAS00001005944

• Genomic profiling of paediatric high grade gliomas from the HERBY clinical trial

  The HERBY trial was a phase-II open-label, randomised, multicentre trial evaluating bevacizumab in patients with newly-diagnosed non-brainstem high grade glioma (HGG) between the ages of 3-18yrs. 121 patients were randomised with 1yr event-free survival as the primary end-point. Confirmation of HGG diagnosis by reference pathology was mandatory before randomisation, followed by review with five independent expert neuropathologists. We collected specimens from 89 patients consenting to translational research, and performed Sanger sequencing for H3F3A, Illumina 450K BeadArray methylation profiling, and whole exome sequencing and RNA sequencing where possible. 7/89 patients (8%) harboured H3F3A G34R/V mutations, whilst 24/89 (27%) harboured H3F3A K27M, the latter reflecting a high proportion of the novel entity recognised in the 2016 WHO classification of diffuse midline glioma with H3K27M mutation. Both histone mutations conferred a significantly shorter progression-free (p=0.0104) and overall survival (p=0.00159). 450K methylation subtyping additionally identified a number of subgroups in histone wild-type cases. These included infrequent (4/86, 4.6%) cases with IDH1 mutation in older children, further categorised into astrocytic ATRX-mutant (n=3) or 1p19q co-deleted with TERT promoter mutation (n=1). 9/74, (12%) of cases were classified as biologically resembling pleomorphic xanthoastrocytoma (PXA) by methylation profiling, with 5/9 harbouring BRAF V600E mutations and epithelioid histology, and 3/9 NF1 mutation and giant cell features. Three cases had methylation profiles more closely resembling low grade gliomas, though were morphologically high grade, and harboured MAPK dysregulation in two cases, with the third part of an additional infant (<3yrs) cohort. Finally, four cases had a mutational burden several orders of magnitude higher than the rest, with (2197-5332) somatic coding variants/sample reflecting a hypermutator phenotype. These data provide an important insight into the wide biological diversity of ‘HGG’ found in the paediatric age group and allow for a profound refinement of clinical trial interpretation.

  Study EGAS00001002328

• High-throughput telomere length measurement at nucleotide resolution using the PacBio high fidelity sequencing platform

  Telomeres are specialized nucleoprotein structures at the ends of linear chromosomes that are essential for the stability as well as the complete replication of the human genome. The progressive shortening of steady-state telomere length in normal human somatic cells is a promising biomarker for age-associated diseases. However, there remain substantial challenges in quantifying telomere length in clinics and research laboratories due to the lack of accurate and high-throughput telomere length measurement method. Here, we describe a novel workflow to capture and analyze telomeres accurately at the single-nucleotide resolution. By using newly designed telobaits that are complementary to the single-stranded G-rich 3’ telomeric overhangs, we can specifically and efficiently capture full-length telomeres and subject them to single-molecule real-time (SMRT) sequencing. We show that the telomere length can be measured at nucleotide resolution using the PacBio high fidelity (HiFi) sequencing platform, which provides the necessary accuracy and throughput to uncover the steady-state telomere length distribution in human culture cell lines as well as clinical patient samples. This new method is a valuable tool not only in population-wide epidemiology studies but also for the longitudinal assessment of telomere attrition in individuals. Our results also reveal the extreme heterogeneity of telomeric variant sequences (TVSs) that is dispersed throughout the telomere repeat region. Importantly, the unique distribution pattern and sequence of TVSs in individuals may be used as a fingerprint for sample identification. In addition, the presence of TVSs disrupts the continuity of the canonical (5’-TTAGGG-3’)n telomere repeats, which could affect the binding of shelterin complexes at the chromosomal ends. Therefore, TVSs may act as a form of genetic polymorphisms that impairs telomere protection. Together with the absolute telomere length in each individual, TVSs may have profound implications in human aging and diseases.

  Study EGAS00001006103

• Local In Time Statistics for processual research

  Background: Functional genomics in a processual analysis cover the time-dependent changes in transcriptomics and epigenetics before diagnosis of a disease, reflecting the changes in both life style and disease processes. The aim of this paper is to explore the dynamic, time-dependent mechanisms of the metastatic processes, using blood transcriptomics and including time in a continuous manner. For achieving this goal we develop new statistical methods based on statistics that are local in time. Methods: The new statistical method, Local In Time Statistics (LITS), is based on calculating statistics in moving windows and randomization. The method has been tested for the analysis of a dataset that collectively provides information on the blood transcriptome up to eight years before breast cancer diagnosis. The dataset from the NOWAC Post-genome Cohort consists of 467 case-control pairs matched on birth year and time of blood sampling. The data for a pair is the difference in log2 gene expression between the case and control. The stratified analyses are based on important biological differences like metastatic versus non-metastatic cancer, and the mode of cancer detection, i.e. screening detected versus clinically detected cancers. The dataset was used for examining whether the gene expression profile varies between cases and controls, with time, or between cases with and without metastases. Results: The null hypotheses of no differences between cases and controls, no time-dependent changes, and no differences between different strata were all rejected. For screening detected cancers the probability of correct prediction of metastasis status was best in year 1 before diagnosis compared to year 3 and 4 before diagnosis for clinically detected cancers. The predictor was not very sensitive to the number of genes included.Conclusions: Using a new statistical method, LITS, we have demonstrated time-dependent changes of the blood transcriptome up to eight years before breast cancer diagnosis.

  Study EGAS00001002520

• RNAseq data to study PRPF6 regulated splice forms in colon cancer cell lines

  Alternative splicing plays critical roles in differentiation, development, and cancer (Pettigrew et al., 2008; Chen and Manley, 2009). The recent identification of specific spliceosome inhibitors has generated interest in the therapeutic potential of targeting this cellular process (van Alphen et al., 2009). Using an integrated genomic approach, we have identified PRPF6, an RNA binding component of the pre-mRNA spliceosome, as an essential driver of oncogenesis in colon cancer. Importantly, PRPF6 is both amplified and overexpressed in colon cancer, and only colon cancer cells with high PRPF6 levels are sensitive to its loss. Our data clearly point to an important role for PRPF6 in colon cancer growth and suggest that a better understanding of its role in alternative splicing in colon cancer is warranted. To determine the specific alternative splice forms that PRPF6 regulates in colon cancer, we plan three experiments: 1. The first involves knocking down expression of PRPF6 in two different cancer cell lines with 3 different siRNAs, and then completing RNA-seq to determine the gene expression changes that occur relative to a non-targeting control siRNA. Because of the role for PRPF6 in pre-mRNA splicing, we especially want to quantify the changes in splice-specific forms of all genes genome-wide to identify genes whose splicing is altered upon PRPF6 knockdown. 2. The second involves immunoprecipitating PRPF6 from two different cancer cell lines and isolating any RNA that is bound to PRPF6, since PRPF6 is an RNA-binding protein. We then want to carry out RNA-seq to identify which RNA molecules co-immunoprecipitated with PRPF6. This will help us determine possible functions for PRPF6 in regulating colon cancer growth. 3. The third involves overexpressing PRPF6 in cell lines and then carrying out RNA-seq to identify any changes in splice-specific gene expression. This will allow us to determine whether increased PRPF6 expression is sufficient to drive alternative splicing changes.

  Study EGAS00001000761

• Genome-wide DNA methylation sequencing identifies epigenetic perturbations in the upper airways under long-term exposure to moderate levels of ambient air pollution

  While the link between exposure to high levels of ambient particulate matter (PM) and increased incidences of respiratory and cardiovascular diseases is widely recognized, recent epidemiological studies have shown that low PM concentrations are equally associated with adverse health effects. As DNA methylation is one of the main mechanisms by which cells regulate and stabilize gene expression, changes in the methylome could constitute early indicators of dysregulated signaling pathways. So far, little is known about PM-associated DNA methylation changes in the upper airways, the first point of contact between airborne pollutants and the human body. Here, we focused on cells of the upper respiratory tract and assessed their genome-wide DNA methylation pattern to explore exposure-associated early regulatory changes. Using a mobile epidemiological laboratory, nasal lavage samples were collected from a cohort of 60 adults that lived in districts with records of low (Simmerath) or moderate (Stuttgart) PM10 levels in Germany. PM10 concentrations were verified by particle measurements on the days of the sample collection and genome-wide DNA methylation was determined by enzymatic methyl sequencing at single-base resolution. We identified 231 differentially methylated regions (DMRs) between moderately and lowly PM10 exposed individuals. A high proportion of DMRs overlapped with regulatory elements, and DMR target genes were involved in pathways regulating cellular redox homeostasis and immune response. In addition, we found distinct changes in DNA methylation of the HOXA gene cluster whose methylation levels have previously been linked to air pollution exposure but also to carcinogenesis in several instances. The findings of this study suggest that regulatory changes in upper airway cells occur at PM10 levels below current European thresholds, some of which may be involved in the development of air pollution-related diseases.

  Study EGAS00001007374

• Assessing the impact of low frequency coding variants on disease risk using the Exomechip

  Following several rounds of Genome Wide Association (GWA) scans and subsequent meta-analyses of results multiple risk loci have been identified for many common diseases. However, in most instances the implicated common variants identified so far explain only a modest fraction of the genetic risk. In parallel to our efforts to identify the causative variants in the known loci it is necessary to assess the full spectrum of sequence variants for disease risk and in particular low frequency and rare variants that have not been tested in a comprehensive way so far. To this end an international effort by investigators who have performed whole exome sequencing in circa 12,000 individuals have assembled a set of ~250,000 exonic variants of low frequency – defined as seen in at least two studies and a minimum of three individuals (non-sysnonymous) or two individuals for SNPs altering splice sites / stop codons. The SNP content from the exome studies was complemented with additional interesting SNP sets totalling 25,000 markers (GWA tag SNPs, grid of common variants, HLA, Mitochondrial, Y-chromosome etc) and were used to generate a custom iSELECT array, the exome chip. Adequately powered association studies to test low frequency variants for association to disease risk are still very expensive if conducted by whole exome sequencing. The exome chip provides a cost efficient way to undertake such an experiment and the Wellcome Trust Cases Control Consortium will be applying this approach to eight diseases, type 1 and type 2 diabetes, coronary artery disease, hypertension, multiple sclerosis, rheumatoid arthritis, bipolar, and ankylosing spondylitis. The first objective will be to generate a large set of at least 10,000 common UK controls as analysis of lower frequency variants will require large sample sizes. This includes ~6000 samples from the 1958 Birth Cohort.

  Study EGAS00001000584

• Whole Genome Sequencing of Asian Lung Cancers: Second Hand Smoke is Not Responsible for Higher Incidence of Lung Cancer Among Asian Never-Smokers

  Background Asian nonsmoking populations have a higher incidence of lung cancer compared to their European counterparts. There is a long-standing hypothesis that the increase of lung cancer in Asian never-smokers is due to environmental factors such as second-hand smoke. Results We analyzed whole-genome sequencing of 30 Asian lung cancers. Unsupervised clustering of mutational signatures separated the patients into two categories: i) all of the never-smokers, and ii) only smokers or ex-smokers. Half of the ex-smokers / smokers were in the never-smoker-like cluster. The overall somatic variant profiles of Asian lung cancers were similar to that of European origin with G.C>T.A being predominant in smokers. We found EGFR and TP53 are the most frequently mutated genes with mutations in 50% and 27% of individuals, respectively. Among these Asian never-smokers, 71% had an EGFR mutation compared to 20% of the smokers in the smoking cluster. Other frequently mutated genes include RYR2, SATB2, C1orf88, FERMT1 and CTNNB1. Somatic alterations occurred in WNT signaling pathway genes, suggesting this pathway plays a role in both Western and Asian lung cancers. Conclusions Asian never-smokers have lung cancer signatures distinct from the smoker signature and their mutation profiles are similar to that of European never-smokers. The profiles of Asian and European smokers are also similar. This suggests that the same mutational mechanisms underlie the etiology for both ethnic groups, and the high incidence of lung cancer in Asian never-smokers might not be due to second hand smoke or other carcinogens that cause oxidative DNA damage. Half of the ex-smokers/smokers have a molecular phenotype similar to never-smokers; of which, 50% had EGFR mutations. This suggests that routine EGFR testing is warranted in the Asian population, regardless of smoking status.

  Study EGAS00001000621

• Probabilistic modeling of personalized drug combinations from integrated chemical screen and molecular data in sarcoma

  Background: Cancer patients with advanced disease exhaust available clinical regimens and lack actionable genomic medicine results, leaving a large patient population without effective treatments options when their disease inevitably progresses. To address the unmet clinical need for evidence-based therapy assignment when standard clinical approaches have failed, we have developed a probabilistic computational modeling approach which integrates sequencing data with functional assay data to develop patient-specific combination cancer treatments. Methods: Tissue taken from a murine model of alveolar rhabdomyosarcoma was used to perform single agent drug screening and DNA/RNA sequencing experiments; results integrated via our computational modeling approach identified a synergistic personalized two-drug combination. Cells derived from the tumor were allografted into mouse models and used to validate the personalized two-drug combination. Computational modeling of single agent drug screening and RNA sequencing of multiple heterogenous sites from a single patient’s epithelioid sarcoma identified a personalized two-drug combination effective across all tumor regions. The heterogeneity-consensus combination was validated in a xenograft model derived from the patient’s primary tumor. Cell cultures derived from human and canine undifferentiated pleomorphic sarcoma were assayed by drug screen; computational modeling identified a resistance-abrogating two-drug combination common to both cell cultures. This combination was validated in vitro via a cell regrowth assay. Results Our computational modeling approach addresses three major challenges in personalized cancer therapy: synergistic drug combination predictions (validated in vitro and in vivo in a genetically engineered murine cancer model), identification of unifying therapeutic targets to overcome intra-tumor heterogeneity (validated in vivo in a human cancer xenograft), and mitigation of cancer cell resistance and rewiring mechanisms (validated in vitro in a human and canine cancer model). Conclusions These proof-of-concept studies support the use of an integrative functional approach to personalized combination therapy prediction for the population of high-risk cancer patients lacking viable clinical options and without actionable DNA sequencing-based therapy.

  Study EGAS00001003564

• McGill Epigenomics Mapping Centre

  The McGill Epigenomics Mapping Centre (EMC) and Data Coordination Centre (EDCC) were established in 2012 at the McGill University and Genome Quebec Innovation Centre (Montreal, Canada) to support large-scale human epigenome mapping for a broad spectrum of cell types and diseases. Previous epigenetic studies have provided profiles of specific chromatin marks in relation to basic biological processes, and future studies of molecular mechanisms must build upon these proofs-of-concept by integrating sequence-based variation with multiple levels of epigenetic and transcriptional regulation across the genome of human tissues and animal disease models. This project leverages the high-throughput sequencing infrastructure and expertise in genomics and transcriptomics at the McGill Innovation Centre to carry out this research. Controlled access to the data by collaborators and the greater scientific community is gained via a portal, which takes advantage of Compute Canada high-performance computing cluster resources to manage the large volume of data associated from the generation of reference epigenome maps. McGill University is one of two Canadian mapping and data coordination centres, the other based at the Michael Smith Genome Sciences Centre in Vancouver, British Columbia. The generation of comprehensive epigenome maps at McGill University is part of a larger international effort that is coordinated by the International Human Epigenome Consortium (IHEC), whose overall long-term objective is to determine the extent to which the epigenome shapes human populations over generations and in response to the environment. This project is funded under the Canadian Epigenetics, Environment, and Health Research Consortium (CEEHRC) by the Canadian Institutes of Health Research and by Genome Quebec, with additional support from Genome Canada. The computing and networking infrastructure, and part of the software development, are provided by Compute Canada and CANARIE.

  Study EGAS00001000995

• Polymorphisms in the mitochondrial genome are associated with bullous pemphigoid in Germans

  Bullous pemphigoid (BP) is the most prevalent autoimmune skin blistering disease and is characterized by the generation of autoantibodies against the hemidesmosomal proteins BP180 (type XVII collagen) and BP230. Most intriguingly, BP is distinct from other autoimmune diseases because it predominantly affects elderly individuals above the age of 75 years, raising the question why autoantibodies and the clinical lesions of BP emerges mostly in this later stage of life, even in individuals harboring known putative BP-associated germline gene variants. The mitochondrial genome (mtDNA) is a potential candidate to provide additional insights into the BP etiology; however, the mtDNA has not been extensively explored to date. Therefore, we sequenced the whole mtDNA of German BP patients (n=180) and age- and sex-matched healthy controls (n=188) using next generation sequencing (NGS) technology, followed by the replication study using Sanger sequencing of an additional independent BP (n=89) and control cohort (n=104). While the BP and control groups showed comparable mitochondrial haplogroup distributions, the haplogroup T exhibited a tendency of higher frequency in BP patients suffering from neurodegenerative diseases (ND) compared to BP patients without ND (p= 0.1448, Fisher’s exact test). A total of four single nucleotide polymorphisms (SNPs) in the mtDNA, namely, m.16263T>C, m.11914G>A, m.16051A>G, and m.16162A>G, were found to be significantly associated with BP based on the meta-analysis of our NGS data and the Sanger sequencing data (p=0.0017, p=0.0132, p=0.0129, and p=0.0076, respectively, Peto’s test). In summary, our study is the first to interrogate the whole mtDNA in BP patients and controls and to implicate multiple novel mtDNA variants in disease susceptibility. Studies using larger cohorts and more diverse populations are warranted to explore the functional consequences of the mtDNA variants identified in this study on immune and skin cells to understand their contributions to BP pathology.

  Study EGAS00001003932

• UK10K COHORT ALSPAC

  The UK10K project proposes 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 directly analyse 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 Avon Longitudinal Study of Parents and Children (ALSPAC) is a two-generation prospective study. Pregnant women living in one of three health districts in the former county of Avon with an expected delivery date between April 1991 and December 1992 were eligible to be enrolled in the study, and this formed the initial point of contact for the development of a large, family based resource. Information has been collected on the children and the mothers through retrieval of biological materials (e.g. antenatal blood samples, placentas), biological sampling (e.g. collection of cord blood, umbilical cord, milk teeth, hair, toenails, blood and urine), self-administered questionnaires, data extraction from medical notes, linkage to routine information systems and at repeat research clinics.

  Study EGAS00001000090

• Genetic_Heterogeneity_of_the_familial_gastric_neuroendocrine_tumors

  Gastric neuroendocrine tumors (gNETs) occur with an estimated frequency of 2 per 100,000 in the general population. Type I gastric neuroendocrine tumors (NETs) represent the 75% of gNTEs and arise from gastric enterochromaffin-like (ECL) cells. They have late age of onset and usually benigh course. Classically, hypergastrinemia in patients who have autoimmune atrophic gastritis, causes hyperplasia of gastric ECL cells that progresses into type I gastric NETs and parietal cell (PC) destruction. The genetic bases in families with this disease are unknown. We performed an exome sequencing study of an atypical aggressive familial gNETs case (with early age onset, nodal infiltrations and gastric adenocarcinomas) that followed a recessive model. We identified a deleterious mutation in homozygosis in the ATP4A gene, which encodes the proton pump responsible for acid secretion by gastric parietal cells. This mutation lead to achlorhydria first, and hypergastrinemia and gNET developing as consequence (Calvete et al. 2014). Recently, two more families with gNETs, classical clinical traits and recessive model have been studies by WES but we didn't find any mutation in the ATP4a gene. However, putative mutations affecting genes that contribute to the development and the integrity of PC have been found suggesting that genetic alterations associated to this disorder target to a unique cell type (parietal cells). In order to cinfirm this hypothesis, it is necessary the search for new genes implicated in the gNETs, more familial cases are needed to be studied. We have identified four more new familial gNETs cases. Here, we propose their study by WES. The first family is formed by thress siblings with gNETs. The other families include two siblings with gNETs. 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 EGAS00001002273

• Genome_Diversity_in_Africa_Project__Uganda

  Globally, human populations show structured genetic diversity as a result of geographical dispersion, selection and drift. Understanding this genetic variation can provide insights into the evolutionary processes that shape both human adaptation and variation in disease. Populations from SSA have the highest levels of genetic diversity. This characteristic, in addition to historical genetic admixture, can lead to complexities in the design of studies assessing the genetic determinants of disease and human variation. However, such studies of African populations are also likely to provide new opportunities to discover novel disease susceptibility loci and variants and refine gene–disease association signals. A systematic assessment of genetic diversity within SSA would facilitate genomic epidemiological studies in the region. The Genome Diversity in Africa Project (GDAP) aims to produce a comprehensive catalogue of human genetic variation in SSA, including single nucleotide polymorphisms (SNPs), structural variants, and haplotypes. This resource will make a substantial contribution to understanding patterns of genetic diversity within and among populations in SSA, as well as providing a global resource to help design, implement and interpret genomic studies in SSA populations and studies comprising globally diverse populations, complementing existing genomic resources. Specifically, we plan to carry out high depth whole genome sequencing of up to 2000 individuals across Africa (25 individuals from each ethnolinguistic group). Our scientific objectives are to: 1) develop a resource that provides a comprehensive catalogue of genetic variation in populations from SSA accessible to the global scientific community; 2) characterise population genetic diversity, structure, gene flow and admixture across SSA; 3) develop a cost-efficient, next-generation genotype array for diverse populations across SSA; and 4) facilitate whole genome-sequencing association studies of complex traits and diseases by developing a reference panel for imputation and resource for enhancing fine-mapping disease susceptibility loci. These scientific objectives will be supported by cross-cutting operational activities, including network and management of the consortium, research ethics, and research capacity building in statistical genetics and bioinformatics

  Study EGAS00001002523

• Next-Generation Sequencing of RNA and DNA Isolated from Paired Fresh-Frozen and Formalin-Fixed Paraffin-Embedded Samples of Human Cancer and Normal Tissue

  Formalin-fixed, paraffin-embedded (FFPE) tissues are an invaluable resource for clinical research. However, nucleic acids extracted from FFPE tissues are fragmented and chemically modified making them challenging to use in molecular studies. We analysed 23 fresh-frozen (FF), 35 FFPE and 38 paired FF/FFPE specimens, representing six different human tissue types (bladder, prostate and colon carcinoma; liver and colon normal tissue; reactive tonsil) in order to examine the potential use of FFPE samples in next-generation sequencing (NGS) based retrospective and prospective clinical studies. Two methods for DNA and three methods for RNA extraction from FFPE tissues were compared and were found to affect nucleic acid quantity and quality. DNA and RNA from selected FFPE and paired FF/FFPE specimens were used for exome and transcriptome analysis. Preparations of DNA Exome-Seq libraries was more challenging (29.5 % success) than that of RNA-Seq libraries, presumably because of modifications to FFPE tissue-derived DNA. Libraries could still be prepared from RNA isolated from two-decade old FFPE tissues. Data were analysed using the CLC Bio Genomics Workbench and revealed systematic differences between FF and FFPE tissue-derived nucleic acid libraries. In spite of this, pairwise analysis of DNA Exome-Seq data showed concordance for 70-80 % of variants in FF and FFPE samples stored for fewer than three years. RNA-Seq data showed high correlation of expression profiles in FF/FFPE pairs (Pearson Correlations of 0.90 +/- 0.05), irrespective of storage time (up to 244 months) and tissue type. A common set of 1,494 genes was identified with expression profiles that were significantly different between paired FF and FFPE samples irrespective of tissue type. Our results are promising and suggest that NGS can be used to study FFPE specimens in both prospective and retrospective archive-based studies in which FF specimens are not available.

  Study EGAS00001000737

• Dynamics of multiple resistance mechanisms in plasma DNA and their clinical implications for NSCLC patients receiving EGFR-targeted therapies

  Background:The development of resistance mechanisms to targeted therapies is complicated by underlying tumour heterogeneity. When multiple resistance mechanisms co-exist, identifying the dominant driver(s) is critical for selection of treatment. Patients and methods:We studied the relative dynamics of multiple oncogenic drivers and resistance mechanisms in 50 patients with EGFR-mutant non-small cell lung cancer (NSCLC), during treatment with gefitinib and hydroxychloroquine (NCT00809237). We performed digital PCR and targeted deep sequencing of cell-free DNA from longitudinal plasma samples from all patients, and shallow whole genome sequencing of serial samples from three patients who underwent histological transformation to small-cell lung cancer (SCLC). ResultsEGFR activating mutations were assessed in tumour samples and accurately identified in plasma samples of 95% of patients (41/43). We identified additional mutations including EGFR T790M (31/50, 62%), TP53 (23/50, 46%), PIK3CA (7/50, 14%), and PTEN (4/50, 8%), and tracked their relative levels in plasma. Patients with both TP53 and EGFR mutations detected in plasma pre-treatment tended to have worse overall survival than those where only EGFR was detected. The resistance-conferring EGFR T790M mutation was identified in 62% of the patients who progressed. Changes in relative levels of T790M and activating mutations in EGFR corresponded to clinical responses in two patients who were given sequential treatments. Patients who progressed without T790M detected in plasma had worse PFS during TKI continuation, and developed alternative resistance mechanisms, that could be identified in plasma DNA, including SCLC-associated copy number changes and TP53 mutations, that tracked with progression and response to subsequent therapies. Conclusion:Longitudinal analysis of multiple oncogenic drivers in ctDNA may help identifying dominant mechanisms of resistance to EGFR-targeted therapies, and highlight the importance of monitoring genetic events that are not direct drug targets but provide non-invasive molecular information that may guide clinical management.

  Study EGAS00001002908

• Population Structure and Genetic Diversity in Argentinean populations

  Argentina was one of the last places peopled when humans entered the Americas. Its earliest archeological sites are dated from 13,000 years BP, and throughout its history the region has exhibited broad cultural diversity, involving both hunter-gatherer and agro-pastoral societies, different linguistic families, the Southern expansion of the Inca Empire, the arrival of the Spanish conquistadors and slave trade from Africa. More recently, Argentina has had the highest proportion of influx of transoceanic migration between 1850 and 1940 of any Latin American country. Furthermore, the geographical distribution of immigrants was oriented towards the Pampa region of central-eastern Argentina and the Northeast, so different patterns of mixed ancestry throughout the country are expected. In this present work we aim to describe the history and structure of the Argentinean population. Thus far, we have collected DNA samples from 2904 participants from 15 dispersed regions in Argentina, from the Andean populations in the Northwest to the river lands of the Northeast. Our preliminary data includes genome-scale Illumina Exome Array 250K genotypes, plus next-generation sequencing data from uniformly spaced genomic regions targeted with a genotype-by- sequencing strategy (which covers 1.5% of the genome), for 375 participants from 9 of the 15 populations. We will identify proportions and patterns of European, African and Indigenous American ancestry in the Argentinean cohort via local ancestry inference. The first series of analyses will focus on genomic tracts of Indigenous American ancestry for fine-scale examination of sub-continental structure by comparing to a panel of 493 samples from 52 Indigenous American populations from throughout the Americas. Furthermore, this work will serve as a test of different low-cost strategies (exome genotype array vs. genotype-by-sequencing) for accurate inference of population structure in extant populations.

  Study EGAS00001001663

• Copy number profiling of primary samples and cell lines of retinoblastoma

  While RB1 loss initiates retinoblastoma development, additional somatic copy number alterations (SCNAs) can drive tumor progression. Although SCNAs have been identified with good concordance between studies at a cytoband resolution, accurate identification of single genes for all recurrent SCNAs is still challenging. This study presents a comprehensive meta-analysis of genome-wide SCNAs integrated with gene expression profiling data, narrowing down the list of plausible retinoblastoma driver genes. We performed SCNA profiling of 45 primary retinoblastoma samples and 8 retinoblastoma cell lines by high-resolution microarrays. We combined our data with genomic, clinical and histopathological data of ten published genome-wide SCNA studies, which strongly enhanced the power of our analyses (N=310). Comprehensive recurrence analysis of SCNAs in all studies integrated with gene expression data allowed us to reduce candidate gene lists for 1q, 2p, 6p, 7q and 13q to a limited gene set. Besides the well-established driver genes RB1 (13q-loss) and MYCN (2p-gain) we identified CRB1 and NEK7 (1q-gain), SOX4 (6p-gain) and NUP205 (7q-gain) as novel retinoblastoma driver candidates. Depending on the sample subset and algorithms used, alternative candidates were identified including MIR181 (1q-gain) and DEK (6p gain). Remarkably, our study showed that copy number gains rarely exceeded change of one copy, even in pure tumor samples with 100% homozygosity at the RB1 locus (N=34), which is indicative for intra-tumor heterogeneity. In addition, profound between-tumor variability was observed that was associated with age at diagnosis and differentiation grades. Since focal alterations at commonly altered chromosome regions were rare except for 2p24.3 (MYCN), further functional validation of the oncogenic potential of the described candidate genes is now required. For further investigations, our study provides a refined and revised set of candidate retinoblastoma driver genes.

  Study EGAS00001001715

• NEOPREDICT-Lung: Neoadjuvant Immunotherapy for Resectable NSCLC

  Links via: DOI: https://doi.org/10.1038/s41591-024-02965-0 PubMed: https://pubmed.ncbi.nlm.nih.gov/38689060/ Journal: https://www.nature.com/articles/s41591-024-02965-0 Abstract Antibodies targeting the immune checkpoint molecules PD-1, PD-L1 and CTLA-4 administered alone or in combination with chemotherapy are standard of care in most patients with metastatic non-small cell lung cancers. When given prior to curative surgery, tumor responses and improved event-free survival are achieved. New antibody combinations may be more efficacious and tolerable. In an ongoing, open-label phase II study, 60 biomarker-unselected, treatment-naïve patients with resectable non-small cell lung cancer were randomized to receive two preoperative doses of nivolumab (anti-PD-1) with or without relatlimab (anti-LAG-3) antibody therapy. The primary study endpoint was feasibility of surgery within 43 days, which was met by all patients. Curative resection was achieved in 95% of patients. Secondary endpoints included pathological and radiographic response rates, pathologically complete resection rates, disease-free and overall survival rates, and safety. Major pathological (≤10% viable tumor cells) and objective radiographic responses were achieved in 27% and 10% (nivolumab), and 30% and 27% (nivolumab and relatlimab) of patients, respectively. In 100% (nivolumab) and 90% (nivolumab and relatlimab) of patients tumors and lymph nodes were pathologically completely resected. With 12 months median duration of follow-up, disease-free survival and overall survival rates at 12 months were 89% and 93% (nivolumab), and 93% and 100% (nivolumab and relatlimab). Both treatments were safe with grade ≥ 3 treatment-emergent adverse events reported in 10% and 13% of patients per study arm. Exploratory analyses provided insights into biological processes triggered by preoperative immunotherapy. This study establishes the feasibility and safety of dual targeting of PD-1 and LAG-3 prior to lung cancer surgery.

  Study EGAS00001007753

• The G2 gene expression signature and MYC overexpression are independent poor prognostic factors in childhood high-grade osteosarcoma

  Purpose In high-grade osteosarcoma, prognostic factors at diagnosis are insufficient for stratifying patients into relevant subgroups. Recently, a transcriptomic study developed the G1/G2 gene expression signature, in which the G2 signature was associated with unfavorable survival. An orthogonal study identified MYC amplification as an unfavorable prognostic factor using targeted next-generation sequencing. The purpose of this study was to validate the independent prognostic value and to investigate the combined prognostic value of the G1/G2 signature with, MYC amplification and/or MYC expression for survival prediction. Material and methods This study included pediatric and adolescent patients with high-grade osteosarcoma. RNA sequencing was performed in 48 patients. Whole-exome sequencing was performed in 40 patients. Gene expression signature scores, MYC amplification (defined as > 7 copies) and MYC expression levels were calculated. Multivariable Cox proportional hazards analysis was performed for event-free survival (EFS) (primary endpoint) and overall survival (OS) (secondary endpoint). Results In the full cohort, the 3-year event-free survival rate was 37%. In multivariable Cox regression analysis with metastatic disease stage (n=21, 44%) as covariate, the G2 signature and MYC expression were independently associated with worse outcomes in terms of EFS (Hazard ratio (HR) 3.32, 95% confidence interval (CI) 1.34-8.21 and HR 3.38, 95%-CI 1.71-6.66, respectively) and OS (HR 4.07, 95%-CI 1.19-13.9 and HR 2.88, 95%-CI 1.22-6.76, respectively). MYC amplification was not associated with EFS or OS in univariable analysis (HR 1.88, 95% -CI 0.74-4.77 and HR 0.79, 95%-CI 0.21-3.05, respectively). Conclusion The G2 gene expression signature and MYC expression were independently associated with unfavorable outcomes in a pediatric cohort of patients with high-grade osteosarcoma. The combined prognostic value warrants further prospective validation and could potentially serve as a stratification marker for future osteosarcoma treatment protocols.

  Study EGAS00001008073

• Whole exome and RNA sequencing of organoid samples derived from TRACERx patients

  Immune checkpoint blockade in non-small cell lung cancer (NSCLC) rarely leads to complete tumour remission, suggesting that tumours harbour subpopulations of cancer cells that can escape immune predation. However, while tumours consist of multiple genetically distinct clones, the extent to which separate tumour subclones differ intrinsically in their functional capacity for immune evasion remains largely unexplored. This has been challenging to address due to an inability to isolate and propagate individual tumour subclones from human cancers. Here, we leverage the multi-region TRACERx lung cancer evolution study to generate a patient-derived platform that allows the functional analysis of immune escape at single clone resolution. We generated libraries of >20 separate NSCLC organoid lines per patient by isolating individual (clonal) organoids established from multiple spatially separated tumour regions. Each organoid subline was co-cultured with autologous tumour infiltrating lymphocytes (TIL) to evaluate their capacity to elicit a T cell response. Organoid lines derived from separate tumour regions, or from single clones within individual regions, differ vastly in their capacity to elicit a T cell response, showing the coexistence of intrinsically immune evading (‘hot’) and non-evading (‘cold’) clones in the same primary tumour. Organoid lines escaping T cell recognition represent genetically and transcriptionally unique subclones with a distinct evolutionary history, showing that evolutionary sculpting of subclonal diversity can have functional consequences for anti-tumour immunity. Notably, the identification of immune-evasive subclones would not have been possible using computational approaches alone, highlighting the unique value of functional co-culture platforms. This patient-derived functional co-culture platform reveals that tumour evolution can give rise to distinct cancer clones with intrinsic differences in immune evasion capacity and provides an approach to prospectively identify, isolate and characterise immune evading subclones from patients with cancer.

  Study EGAS00001008092

• Raw bulk TCRseq data fromTIL and non-TIL expanded TRACERx samples in the manuscript Bulk TCRseq from TRACERx samples from 'Subclonal immune evasion in non-small cell lung cancer'

  Immune checkpoint blockade in non-small cell lung cancer (NSCLC) rarely leads to complete tumour remission, suggesting that tumours harbour subpopulations of cancer cells that can escape immune predation. However, while tumours consist of multiple genetically distinct clones, the extent to which separate tumour subclones differ intrinsically in their functional capacity for immune evasion remains largely unexplored. This has been challenging to address due to an inability to isolate and propagate individual tumour subclones from human cancers. Here, we leverage the multi-region TRACERx lung cancer evolution study to generate a patient-derived platform that allows the functional analysis of immune escape at single clone resolution. We generated libraries of >20 separate NSCLC organoid lines per patient by isolating individual (clonal) organoids established from multiple spatially separated tumour regions. Each organoid subline was co-cultured with autologous tumour infiltrating lymphocytes (TIL) to evaluate their capacity to elicit a T cell response. Organoid lines derived from separate tumour regions, or from single clones within individual regions, differ vastly in their capacity to elicit a T cell response, showing the coexistence of intrinsically immune evading (‘cold) and non-evading (‘hot) clones in the same primary tumour. Organoid lines escaping T cell recognition represent genetically and transcriptionally unique subclones with a distinct evolutionary history, showing that evolutionary sculpting of subclonal diversity can have functional consequences for anti-tumour immunity. Notably, the identification of immune-evasive subclones would not have been possible using computational approaches alone, highlighting the unique value of functional co-culture platforms. This patient-derived functional co-culture platform reveals that tumour evolution can give rise to distinct cancer clones with intrinsic differences in immune evasion capacity and provides an approach to prospectively identify, isolate and characterise immune evading subclones from patients with cancer.

  Study EGAS00001008118

• Successful immune checkpoint blockade in a patient with advanced stage microsatellite unstable biliary tract cancer (H021)

  Cancers acquire multiple somatic mutations that can lead to the generation of immunogenic mutation-induced neoantigens. These neoantigens can be recognized by the host’s immune system. However, continuous stimulation of immune cells against tumor antigens can lead to immune cell exhaustion, which allows uncontrolled outgrowth of tumor cells. Recently, immune checkpoint inhibitors have emerged as a novel approach to overcome immune cell exhaustion and re-activate anti-tumor immune responses. In particular, antibodies blocking the exhaustion-mediating programmed death receptor (PD-1)/PD-L1 pathway have shown clinical efficacy. The effects were particularly pronounced in tumors with DNA mismatch repair deficiency and a high mutational load, which typically occur in the colon and endometrium. Here, we report on a 24-year old woman diagnosed with extrahepatic cholangiocarcinoma who showed strong and durable response to the immune checkpoint inhibitor pembrolizumab, although treatment was initiated at an advanced stage of disease. The patient’s tumor displayed DNA mismatch repair deficiency and microsatellite instability (MSI), but lacked other features commonly discussed as predictors of response towards checkpoint blockade, such as PD-L1 expression or dense infiltration with cytotoxic T cells. Notably, high levels of HLA class I and II antigen expression were detected in the tumor, suggesting a potential causal relation between functionality of the tumor’s antigen presentation machinery and the success of immune checkpoint blockade. We suggest determining MSI status in combination with HLA class I and II antigen expression in tumors potentially eligible for immune checkpoint blockade even in the absence of conventional markers predictive for anti-PD-1/PD-L1 therapy and in entities not commonly linked to MSI phenotype. Further studies are required to determine the value of these markers for predicting the success of immune checkpoint blockade. (H021)

  Study EGAS00001002441

• Exome Sequencing to Define the Landscape of Plasma Cells in Systemic Light chain Amyloidosis

  Systemic light chain amyloidosis (AL) is characterized by the deposition of immunoglobulin light chains as amyloid fibrils in different organs, where they form toxic protein aggregates. The underlying disease is a plasma cell disorder, but limited whole exome data are available. We report the findings of an exome sequencing study in AL to define a plasma cell signature and compare this to monoclonal gammopathy of undefined significance (MGUS) and myeloma (MM). Twenty-four samples from unselected newly diagnosed untreated AL patients were analysed. CD138+ cells were isolated from bone marrow cells using MACSorting (Miltenyi Biotech, Bisley, UK). Cells were lysed in RLT+ buffer and DNA/RNA extracted using the AllPrep kit (Qiagen, Manchester, UK). Peripheral blood was isolated from patients, white blood cells purified by Ficoll-Pacque and DNA extracted using the QIAamp DNA mini kit (Qiagen). There were a median of 47 acquired nonsynonymous variants (range 4-253) per patient. Thirty percent of patients yielded a mutation in a myeloma driver gene. These mutations were either clonal or sub-clonal. There was evidence of MAPK activation with NRAS mutations, NF-κB activation and DNA repair pathway alterations. Exome data were used to determine the cytogenetic groups of AL samples and identified hyperdiploidy in 33% of cases and t(11;14) in 33% of cases. All translocations were generated via class switch recombination and breakpoints on chromosome 11 fell into the expected regions. A novel translocation, unlikely to be driving disease progression was identified involving RCC1. Finally, further analysis of the APCS gene encoding for the SAP protein did not show any biased genotype in the AL population in comparison to the general population. We conclude that exome sequencing identifies a genetic make-up of AL that is similar to other plasma cell disorders such as MM and MGUS. This not only includes copy number abnormalities and translocations but also nonsynonymous mutations.

  Study EGAS00001001418

• H3Africa - Genomic and Environmental Risk Factors for Cardiometabolic Disease in Africans

  The long-term vision of the AWI-Gen Collaborative Centre (CC) is to build sustainable capacity in Africa for research that leads to an understanding of the interplay between genetic, epigenetic and environmental risk factors for obesity and related cardiometabolic diseases (CMD) in sub-Saharan Africa. The CC will be consolidated under the auspices of the University of the Witwatersrand (Wits) and the International Network for the Demographic Evaluation of Populations and Their Health in Developing Countries (INDEPTH). It will capitalize on the unique strengths of existing longitudinal cohorts, including the urban Soweto and rural Agincourt studies in South Africa (Wits based), and the well established INDEPTH demographic health and surveillance centers in Kenya, Ghana, Burkina Faso and South Africa. The centers offer established infrastructure, trained fieldworkers, long-standing community engagement, and detailed longitudinal phenotypic data, focusing on obesity and cardiometabolic health. Key strengths are harmonized phenotyping across sites, building on strong existing cohorts, and representation of the geographic and social variability of African populations. We aim to: 1. Build sustainable infrastructure (biobanks and laboratories) and capabilities (well characterized population cohorts, genotyping and bioinformatics) for genomic research on the African continent; 2. Understand the genomic architecture of sub-Saharan populations from west, east and south Africa to guide genomic studies (genome sequencing and high throughput SNP and CNV arrays using unrelated individuals and family trios to improve the accuracy of haplotype analyses) and; 3. Investigate the independent and synergistic genomic contributions to body fat distribution (BMI, hip/waist circumference, subcutaneous and visceral fat) in these populations considering the relevant environmental and social contexts (rural/urban communities, quickly transitioning obesity prevalence, differential HIV, TB, and malaria infection histories). We will investigate the effect of obesity and fat distribution on the risk for CMD in the longitudinal cohorts.

  Study EGAS00001002482

• UK10K_NEURO_ASD_FI

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches.These samples are a subset of a nationwide collection of Finnish autism spectrum disorder (ASD) samples. The samples have been collected from Central Hospitals across Finland in collaboration with the University of Helsinki. The samples consist of 93 individuals with a diagnosis of autistic disorder or Asperger syndrome from 36 families with at least two affected individuals. Of these individuals, 16 can be genealogically connected to form two large pedigrees originating from Central Finland, suggesting possible genetic risk factors shared identical by descent within the pedigrees. All diagnoses are based on ICD-10 and DSM-IV diagnostic criteria for ASDs. Additional phenotypic data is available for a subset of the individuals.For further information with regard to this cohort please contact Aarno Palotie (Aarno.palotie@helsinki.fi).

  Study EGAS00001000110

• Liver biopsy derived induced pluripotent stem cells provide an unlimited supply for the generation of patient-specific hepatocyte-like cells

  Background: Hepatocyte-like cells (HLCs) differentiated from induced pluripotent stem cells (iPSCs) have emerged as a promising cell culture model to study metabolism, biotransformation, viral infections and genetic liver diseases. However the derivation of iPSCs is restricted by ethical and procedural constraints. Furthermore, incomplete HLC differentiation remains a major challenge. iPSC may carry-on a tissue of origin dependent expression memory influencing iPSC differentiation into different cell types. If liver derived iPSCs (Li-iPSCs) maintain a liver specific gene expression profile is not known. It is also not known, if Li-iPSCs would allow the generation of more fully differentiated HLCs.Methods: Primary liver cells (PLCs) were expanded from liver needle biopsies and reprogrammed into Li-iPSCs using a non-integrative Sendai virus-based system. Li-iPSCs were differentiated into HLCs using established differentiation protocols. The HLC phenotype was characterized at the protein, functional and transcriptional level. RNA sequencing data were generated from the originating liver biopsies, the Li-iPSCs, fibroblast derived iPSCs, and differentiated HLCs, and used to characterize and compare their transcriptome profiles.Results: PLCs were obtained from small pieces of liver biopsies and could be reprogrammed into Li-iPSCs. Li-iPSCs indeed retain a liver specific transcriptional footprint. Li-iPSCs can be propagated to provide an unlimited supply of cells for differentiation into Li-HLCs. Similar to HLCs derived from fibroblasts, Li-HLCs could not be fully differentiated into hepatocytes. Relative to the originating liver, Li-HLCs showed lower expression of liver specific transcription factors and increased expression of genes involved in the differentiation of other tissues.Conclusions: PLCs and Li-iPSCs obtained from small pieces of human needle liver biopsies constitute a novel unlimited source for the production of HLCs. Despite the preservation of a liver specific gene expression footprint in Li-iPSCs, the generation of fully differentiated hepatocytes cannot be achieved with the current differentiation protocols.

  Study EGAS00001002676

• Characterization of genetic intratumor heterogeneity in colorectal cancer and matching patient-derived spheroid cultures.

  Patient-derived in vitro cultures of colorectal cancer (CRC) may help guide treatment strategies prior to patient treatment. However, most previous studies have been performed on a single biopsy per tumor. The purpose of this study was to analyze multiple spatially distinct biopsies from CRCs and see how well intratumor heterogeneity (ITH) was recapitulated in matching patient-derived spheroids. Three to five biopsies were collected from six CRC tumors. Each biopsy was split in two; one half was used for spheroid culturing, the other half was used for direct DNA and RNA purification. For two of the patients, lymph node metastases were analyzed. Somatic mutations were called from whole exome sequencing data. Each tumor contained multiple mutations shared across all biopsies and spheroids, including major CRC drivers such as APC, KRAS, and TP53. At the same time, all tumors exhibited ITH on both mutation and copy number level. The concordance between biopsies and spheroids ranged between 40-70% for coding mutations. For three of the patients the biopsy and spheroid from matching areas clustered together, meaning that the spheroid resembled the area of origin more than the other areas. However, all biopsies and spheroids contained private mutations. Thereby, multiple cultures from spatially distinct sites of the tumor increase the insight into the genetic profile of the entire tumor. RNA sequencing based molecular subtyping was stable and seemed independent of genetic ITH. In conclusion, all examined CRC tumors contained genetic ITH on mutation and copy number level. Spheroid cultures partly reflected this ITH, meaning that multiple cultures from distinct sites of the tumor improved the representation of the different genetic subclones in the primary tumor. This should be taken in to account when establishing patient-derived models for drug screening.

  Study EGAS00001002684

• Intra-tumor heterogeneity of localized lung adenocarcinomas defined by multi-region sequencing

  Background: Intra-tumor heterogeneity (ITH) may have profound impacts on biopsy strategy, characterization of actionable targets, treatment planning, and drug resistance. ITH of lung cancers has not been well studied. Methods: We performed multi-region whole exome sequencing (average sequencing depth 277x) on 48 tumor regions from 11 resected lung adenocarcinomas (8 stage I, 2 stage II and one stage III tumors) followed by deep sequencing validation (average sequencing depth 863x). Results: All tumors showed clear evidence of spatial ITH. On average, 76% of all mutations and 20/21 known cancer gene mutations were present in all regions of individual tumors (trunk mutations). Copy number analysis did not show substantial difference in large-scale chromosome aberrations. A canonical KRAS p.G12C mutation was detected in only 1 of 4 tumor regions (branch mutations) at exome depth but detected in all four tumor regions (trunk mutations) at increased sequencing depth. In addition, many heterogeneous branch mutations defined by exome sequencing were detectable in all regions of individual tumors with increasing sequencing depth. Significant differences in mutational spectrum were observed in 6 tumors indicating that specific mutational processes were likely operative at different times during development of these tumors. With a median follow up of 21 months post surgery, 3 patients have relapsed and all 3 patients presented with significantly larger fractions of subclonal mutations in their primary tumors. Conclusions: Although multi-region sampling is still required to fully assess ITH complexity, single biopsy analysis at appropriate depth might be sufficient to identify the majority of known cancer gene mutations in localized lung adenocarcinomas. Larger subclonal mutation fractions may be associated with increased likelihood of postsurgical relapse in patients with localized lung adenocarcinomas. Studies on larger cohort, ideally with comprehensive clinical annotation and repeat biopsy at relapse are needed to fully understand the clinical impact of ITH.

  Study EGAS00001000930

• Exploration of CNV’s and SNV’s in cancers with well-known genetic rearrangements: Identification of additional genetic changes in rearrangements-driven cancer

  Specific cancers are well known to have specific gene rearrangements. Representative cancer is chronic myeloid leukemia (CML) which is famous for BCR/ABL gene rearrangements. By the way, also genetic rearrangements are not infrequently found in acute leukemia, and those include CBFB/MYH11 rearrangements, PML/RARA rearrangements, AML/ETO rearrangements, and BCR/ABL rearrangements. The situation is same is for lung adenocarcinoma of non-smokers and prostate cancers, and ALK and RET gene rearrangements are recently found. However, not all these rearrangements are sufficient to develop cancer, and it is suggested that some genetic aberrations are required for the development and progression of cancer in addition to these genetic rearrangements. In this study, we have a plan to explore the existence of additional genetic changes including single nucleotide variants (SNV’s) and copy number variations (CNV’s) in cancers with specific genetic rearrangements so as to define additional genetic changes that confer full oncogenic potential to cancer cells. So as to look at SNV’s and small indels, we are going to use whole exome sequencing (WES) data. And as for CNV’s, we are going to use array-CGH data. Also, to check the functionality of these genetic changes, we are going to use whole transcriptome data (WTS). In addition, in samples with whole genome sequencing (WGS ) data, we are going to explore genetic changes in intronic regions also. In fact, we sequenced 10 AML samples which are known to have specific genetic rearrangements. Analysis of these WGS, WES and WTS data are under the way. If we have a chance to access WES, array-CGH and WTS data of other cancer samples with specific genetic rearrangements, we would be able to compare additional genetic changes required for oncogenesis in cancer cells with various genetic rearrangements. We hope this result will enlighten our understanding of cancer genetics further.

  Study EGAS00001000673

• Whole exome sequencing of germline DNA was performed and subsequent polymorphisms in genes known and putatively involved in the innate immune response to fungi were identified

  Disease can be conceptualized as the result of interactions between infecting microbe and holobiont, the combination of a host and its microbial communities. It is likely that genomic variation in the host, infecting microbe, and commensal microbiota are key determinants of infectious disease clinical outcomes. However, until recently, simultaneous, multiomic investigation of infecting microbe and holobiont components has rarely been explored. Herein, we characterized the infecting microbe, host, micro- and mycobiomes leading up to infection onset in a leukemia patient that developed invasive mucormycosis. We discovered that the patient was infected with a strain of the recently described Mucor velutinosus species which we determined was hypervirulent in a Drosophila challenge model and has a predisposition for skin dissemination. After completing the infecting M. velutinosus genome and genomes from four other Mucor species, comparative pathogenomics was performed and assisted in identifying 66 M. velutinosus-specific putatively secreted proteins, including multiple novel secreted aspartyl proteinases which may contribute to the unique clinical presentation of skin dissemination. Whole exome sequencing of the patient revealed multiple non-synonymous polymorphisms in genes critical to control of fungal proliferation, such as TLR6 and PTX3. Moreover, the patient had a non-synonymous polymorphism in the NOD2 gene and a missense mutation in FUT2, which have been linked to microbial dysbiosis and microbiome diversity maintenance during physiologic stress, respectively. In concert with host genetic polymorphism data, the micro- and mycobiome analyses revealed that the infection developed amid a dysbiotic microbiome with low α-diversity, dominated by staphylococci. Additionally, longitudinal mycobiome data showed that M. velutinosus DNA was detectable in oral samples preceding disease onset. Our genome-level study of the host-infecting microbe-commensal triad extends the concept of personalized genomic medicine to the holobiont-infecting microbe interface thereby offering novel opportunities for using synergistic genetic methods to increase understanding of infectious diseases pathogenesis and clinical outcomes.

  Study EGAS00001001542

• Bulk-RNA Sequencing of high-grade pancreatic and non-pancreatic Neuroendocrine Neoplasms

  Therapeutic decisions in oncology depend on a precise pathological classification of individual neoplasms. Recent years have seen an intensification of research activities aimed at the extraction of clinically relevant information from patient-derived 'omics' data based on Machine-Learning models. However, a comprehensive training of Machine-Learning models requires sufficiently large numbers of training samples, which are usually not available for rare cancer types. The problem is worsened when individual tissues segregate into different cancer subtypes, as their discrimination would require even more training samples. Methods: Here, we report on a new data-augmentation technique to support the training of Machine-Learning models on ‘omics’ data from pancreatic neuroendocrine neoplasms (panNEN). PanNENs display all properties described above: Only about 2-3% of all pancreatic neoplasms are neuroendocrine and they fall into different subtypes with distinctly different prognosis, which makes the precise classification of such samples both difficult and important for therapy decisions. The approach reconstructs a given transcriptome based on healthy pancreatic cell type signatures and creates a Machine-Learning model that integrates the observed reconstruction error and predicted cell type proportions as features. Results: A benchmark of the deconvolution model predictions with the ground-truth found that the model could efficiently predict the sample grading, disease-related patient survival time, and differentiate between different subtypes in four panNEN and one mixed panNEN and non-pancreatic NENs dataset. We compared the predictive performance of the deconvolution-trained model to that of a model trained directly on the transcriptomic data, under inclusion of the Ki-67 classification gold standard biomarker, and found the performances to be comparable. Conclusion: Our approach serves as a data-augmentation technique to facilitate the training of Machine Learning models for rare cancer types via utilization of data from healthy origin. Additionally, the deconvolution results were clinically interpretable and further research can render the deconvolution approach an effective complementary asset for the clinical classification of neoplasms.

  Study EGAS00001004861

• Molecular Tumor Board to Inform on Personalized Medicine for a Man with Advanced Prostate Cancer

  Purpose: Molecular profiling of cancer is increasingly common as part of routine care in oncology, and germline and somatic profiling may provide insights and actionable targets for men with metastatic prostate cancer. However, all reported cases are of deidentified individuals without full medical and genomic data available in the public domain. Patient and Methods: We present a case of whole genome tumor and germline sequencing in a patient with advanced prostate cancer, who has agreed to make his genomic and clinical data publicly available in an accompanying personal testimonial. Results: We describe an 84-year-old Caucasian male with a Gleason 10 oligometastastic hormone-sensitive prostate cancer. Whole-genome sequencing provided insights into his tumor’s underlying mutational processes and the development of an SPOP mutation. It also revealed an androgen-receptor dependency of his cancer which was reflected in his durable response to radiation and hormonal therapy. Potentially actionable genomic lesions in the tumor were identified through a personalized medicine approach for potential future therapy, but at the moment, he remains in remission, illustrating the hormonal sensitivity of his SPOP-driven prostate cancer. We also placed this patient in the context of a large prostate-cancer cohort from the PCAWG (Pan-cancer Analysis of Whole Genomes) group. In this comparison, the patient's cancer appears typical in terms of the number and type of somatic mutations, but it has a somewhat larger contribution from the mutational process associated with aging. Conclusion: We combined the expertise of medical oncology and genomics approaches to develop a molecular tumor board to integrate the care and study of this patient, who continues to have an outstanding response to his combined modality treatment. This identifiable case potentially helps overcome barriers to clinical and genomic data sharing.

  Study EGAS00001004648

• Genetic analysis of HLA and immune escape genes in Diffuse Large B-cell Lymphoma

  Fifty percent of diffuse large B-cell lymphoma (DLBCL) cases lack cell-surface expression of the class-I major histocompatibility complex (MHC-I), thereby escaping immune recognition by cytotoxic T cells. In order to comprehensively identify the mechanisms involved in MHC-I loss, we first performed immunophenotypic analysis of both MHC class-I and -II in 657 cases across the spectrum of B-NHL, revealing that loss of MHC-I, but not MHC-II, is preferentially restricted to DLBCL. We then used whole exome and targeted deep-sequencing to examine genes involved in MHC-I expression in 74 DLBCL samples representative of MHC-I positive and negative cases. We show here that somatic biallelic or monoallelic inactivation of B2M and/or HLA-I is present in 80.9% (34/42) of MHC-I negative tumors. Furthermore, 68.8% (22/32) of MHC-I positive DLBCLs also harbored monoallelic HLA-I genetic alterations (MHC-I positivemono) that lead to allelic imbalance, suggesting allele-specific inactivation. Both MHC-I negative and MHC-I positivemono cases showed a significantly higher mutational burden as well as a higher number of inferred neo-antigens, suggesting co-selection of HLA-I loss and sustained neo-antigen production. Interestingly, the analysis of > 500.000 individuals in two databases revealed homozygosity of germline HLA-I genes in 26-38% of DLBCL patients, a frequency significantly higher than that observed in any other cancer type. In mice, germinal-center B cells lacking HLA-I expression did not progress to lymphoma and were counterselected in the context of oncogene-driven lymphomagenesis, suggesting that additional events are needed to license immune evasion. These results suggest a multi-step process of HLA-I loss including both germ-line and somatic events in DLBCL development, and have direct implications for the pathogenesis and immunotherapeutic targeting of this disease.

  Study EGAS00001005054

• Same-day genomic and epigenomic diagnosis of brain tumors using realtime nanopore sequencing

  Molecular classification of cancer has entered clinical routine to inform diagnosis, prognosis and treatment decisions. At the same time, new tumor entities have been identified that cannot be defined histologically. For central nervous systems tumors, the current World Health Organization classification explicitly demands molecular testing, e.g. for 1p/19q- codeletion or IDH mutations, to make an integrated histomolecular diagnosis. However, a plethora of sophisticated technologies is currently needed to assess different genomic and epigenomic alterations and turnaround times are in the range of weeks, which makes standardized and widespread implementation difficult and hinders timely decision making. Here, we explored the potential of a pocket-size nanopore sequencing device for multimodal and rapid molecular diagnostics of cancer. Low-pass whole genome sequencing was used to simultaneously generate copy number (CN) and methylation profiles from native tumor DNA in the same sequencing run. Single nucleotide variants in IDH1, IDH2, TP53, H3F3A and the TERT promoter region were identified using deep amplicon sequencing. Nanopore sequencing yielded ~0.1X genome coverage within six hours and resulting CN and epigenetic profiles correlated well with matched microarray data. Diagnostically relevant alterations, such as 1p/19q codeletion, and focal amplifications could be recapitulated. Using ad hoc random forests, we could perform supervised pan-cancer classification to distinguish gliomas, medulloblastomas and brain metastases of different primary sites. Single nucleotide variants in IDH1, IDH2 and H3F3A were identified using deep amplicon sequencing within minutes of sequencing. Detection of TP53 and TERT promoter mutations shows that sequencing of entire genes and GC-rich regions is feasible. Nanopore sequencing allows same-day detection of structural variants, point mutations and methylation profiling using a single device with negligible capital cost. It outperforms hybridization-based and current sequencing technologies with respect to time-to diagnosis and required laboratory equipment and expertise, aiming to make precision medicine possible for every cancer patient, even in resource restricted settings.

  Study EGAS00001002213

• Deep single-cell RNA sequencing data for 11138 T cells from tumour, adjacent normal tissue and peripheral blood of treatment-naive CRC patients. The DATA ACCESS AGREEMENT is provided at https://github.com/zhangyybio/single-T-cell-data-access. Applicants can request access to the data by directly downloading it or by sending an email to cancerpku@pku.edu.cn. The process that is used to approve an application includes verifying the institution, participants and research purposes of the application. In general this process will take about two weeks. In principal, any scientific research program complying with the laws and bioethic regulation policies of China will be approved.

  T cells are central players in cancer immunotherapy1, yet some of their fundamental properties such as development and migration within tumours remain elusive. The enormous T cell receptor (TCR) repertoire, required for recognising foreign and self-antigens2,3, could serve as lineage tags to track these T cells in tumours4. Here, we obtained transcriptomes of 11,138 single T cells from 12 colorectal cancer (CRC) patients and developed STARTRAC (Single T-cell Analysis by Rna-seq and Tcr TRACking) indices to quantitatively analyse dynamic relationships among 20 identified T cell subsets with distinct functions and clonalities. While both CD8+ effector and “exhausted” T cells exhibited high clonal expansion, they were independently connected with tumour-resident CD8+ effector memory cells, implicating a TCR-based fate decision. Of the CD4+ T cells, the majority of tumour-infiltrating Tregs showed clonal exclusivity, whereas certain Treg clones were developmentally linked to multiple TH clones. Notably, we identified two IFNG+ TH1-like clusters in tumours, the GZMK+ TEM and CXCL13+ TH1-like clusters, which were associated with distinct IFN-γ-regulating transcription factors, EOMES/RUNX3 and BHLHE40, respectively. Only CXCL13+BHLHE40+ TH1-like cells were preferentially enriched in tumours of microsatellite-instable (MSI) patients, which might explain their favourable response rates to immune-checkpoint blockade. Furthermore, we found IGFLR1 to be highly expressed in both CXCL13+BHLHE40+ TH1-like and CD8+ exhausted T cells and possessed co-stimulatory functions. Our integrated STARTRAC analyses provides a powerful avenue to comprehensively dissect the T cell properties in CRC, which could shed new insights into the dynamic relationships of T cells in other cancers

  Study EGAS00001002791

• Molecular characterization of Barrett’s esophagus at single cell resolution

  Backgrounds & Aims: Patient with Barrett’s esophagus (BE) have an increased risk to develop esophageal adenocarcinoma and are under periodical endoscopic surveillance for malignant transformation. Pathologists distinguish between different BE stages to identify high-risk patients, but this process remains difficult as stage-specific markers are still missing. We established BE organoid cultures, applied single cell sequencing approaches to identify differences in gene expression and genomic alterations, and validated candidate genes for their predictive value on histological sections. Methods: We collected 43 epithelial biopsies from 21 patients, established organoids for the clonal expansion of matching healthy and diseased tissue from two patients to identify genomic alterations by whole genome sequencing (WGS). We performed single cell DNA-sequencing (scDNAseq) to analyse DNA alterations within biopsies from 8 patients and single cell RNA-sequencing (scRNAseq) to identify gene expression differences between BE stages from 19 patients. Candidate marker genes were then validated by RNA in situ hybridization experiments on histological resection specimen. Results: Most BE biopsies were chromosomal stable (CS) and their single base substitution (SBS) signatures were indistinguishable from the healthy control tissues. Dysplastic BE cells contained areas of chromosomal instabilities (CIN), which differed by the presence of the SBS17 signatures from biopsy-matching CS cells or patient-matching healthy control cells. Such CIN areas affected efficient clustering of the scRNAseq data. We identified two sets of marker genes (SLC5A5, PSCA, LIPF and ANPEP, CEACAM6, REG4), which distinguish columnar BE epithelium from non-dysplastic/dysplastic stages. CLDN2 allowed the distinction between dysplastic and more advanced stages as its expression level increased and spread in dysplastic BE glands. Conclusions: Molecular characterization of BE samples identified a correlation between the presence of CIN and COSMIC signature 17, and we report seven marker genes that help in distinguishing different BE stages.

  Study EGAS00001005221

• WES_of_adult_intellectual_disabilities_with_co_morbid_psychiatric_disorders

  ost adults with intellectual disabilities (ID) do not undergo genetic diagnostic investigation as part of their clinical care and have 'missed the boat' with regard to the WES and WGS genetic testing that is now being provided for children with ID. There is a dramatically increased risk of psychistric disorders in adults with ID, e.g. the risk of psychoses is 10X higher than in the general population. It remains an open question as to how much of adult ID is genetic in origin and how similar the genetic forms of adult ID are to those being diagnosed in children, in part due to survivor bias. There is also the opportunity to identify adults with treatable forms of ID, of which over 80 have been described, thus improving their clinical management. Furthermore, analysis of medical records of adults with genetic forms of ID can help to characterise the 'natural history' of individual disorders, resulting in more accurate prognoses for diagnosed children and identifying opportunities for improved management and possibly therapeutic intervention (e.g. optimal anti-epileptic therapy). Here we propose to exome sequence (to ~50X coverage) 200 adults with ID and co-morbid psychiatric disorders. This cohort has previously been assayed with chromosomal microarrays (Wolfe et al 2017 EJHG, 25, 66–72) identifying a diagnostic yield of ~11% which is comparable to the CNV diagnostic yield in various child ID cohorts (10-15%). The authors observed no substantive biases in diagnostic yield between different psychiatric diagnostic classes. The WES data will be analysed using the diagnostic workflows developed in the DDD study to ensure comparability between child and adult ID datasets. This study is intended as a pilot study to demonstrate the value of WES in adults with ID. 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 EGAS00001002962

• Short and long-read genome sequencing methodologies for somatic variant detection; genomic analysis of a patient with diffuse large B-cell lymphoma

  Recent advances in throughput and accuracy mean that the Oxford Nanopore Technologies (ONT) PromethION platform is a now a viable solution for WGS. New bioinformatic methods have been developed to take advantage of this long read data, however much of the validation of these tools has focussed on calling germline variants (both SNVs and structural variants). Somatic variants are outnumbered many-fold by germline variants and their detection is further complicated due to their frequency varying depending on tumour purity/subclonality. Here, we evaluate the extent to which Nanopore WGS enables genome-wide detection and analysis of somatic variation. We do this through sequencing tumour and germline genomes for a patient with diffuse B-cell lymphoma. We examine the capability of currently available tools for calling somatic variants in ONT data by comparing the data with results from 150bp short-read sequencing of the same samples. We then conduct a detailed analysis of the performance of multiple long-read mappers and structural variant callers for calling large, somatic structural variants (SVs) in ONT data. Our protocol achieved yields of up to 96 mapped Gb per PromethION flow cell with average read lengths of ~5kb. Calling germline SNVs from these data achieved good specificity and sensitivity. However, results of somatic SNV calling highlight the need for the development of specialized joint calling algorithms. Our analysis of structural variants shows that the comparative performance of different tools varies significantly between SV types, and suggest long reads are especially advantageous for calling large somatic deletions and duplications. Finally, we highlight the utility of long reads for phasing clinically relevant variants by using the ONT data to confirm that a somatic 1.6Mb deletion and a p.(Arg249Met) mutation involving TP53 are oriented in trans.

  Study EGAS00001004266

• Reduced Representation Bisulfite Sequencing (RRBS) in various tissues to discover DNA methylation markers for the diagnosis of breast and ovarian cancer.

  Despite a myriad of attempts in the last three decades to diagnose ovarian cancer (OC) earlier, this clinical aim still remains a significant challenge. In addition, monitoring treatment and early detection of fatal breast cancer (BC) remains a major unmet need. Aberrant methylation patterns of linked CpGs analyzed in DNA fragments shed by cancers into the bloodstream (i.e. cell-free DNA) can provide highly specific signals indicating cancer presence.In order to discover the most relevant CpG regions, we used RRBS and analysed several tissue samples: 11 prospectively collected invasive epithelial ovarian cancer samples (high grade serous n=8, low grade serous n=1, endometrioid n=1, mucinous n=1, mean age = 54.7 years), 8 prospectively collected invasive ductal breast cancer samples (2/8 triple negative; mean age = 56.6 years), one benign tumor (papillary serous cystadenoma, age = 86 years), 18 non-neoplastic tissue samples (breast n=7 and adnexal n=11, mean age = 60.2 years), two non-neoplastic endometrial tissues (mean age = 68 years) and twenty three white blood cell samples (mean age = 57.8 years) were assessed by RRBS. Genome wide methylation analysis was performed by Reduced Representation Bisulfite Sequencing (RRBS) at GATC Biotech. DNA was digested with MspI followed by size selection of the library, providing enhanced coverage for the CpG-rich regions. The digested DNA was adapter ligated, bisulfite modified and PCR amplified. The libraries were sequenced on Illumina’s HiSeq 2500. Analysis of the first samples sequenced with 100bp paired-end mode showed that the library insert size was small. Therefore, the remaining samples were sequenced with 50 bp paired-end mode. Using Genedata Expressionist® for Genomic Profiling v9.1, we established a bioinformatics pipeline for the detection of cancer specific differentially methylated regions (DMRs). The most promising DMRs were taken forward for the development and validation of serum based clinical assays.

  Study EGAS00001002609

• Multiple Myeloma Total Therapy trial patient sequencing

  The study of multiple myeloma (MM) genomics has identified many abnormalities that are associated with poor progression free survival (PFS) and overall survival (OS). Copy number abnormalities have been extensively studied in many datasets with long follow-up, however, the prognostic impact of mutations have not been extensively studied and available datasets have generally had a relatively short follow-up. These analyses have identified a range of mutations that are associated with prognosis, making it important to extend these observations in larger studies with robust diagnostic technologies. Samples from newly diagnosed MM patients enrolled in Total Therapy trials (n=225) were sequenced on a targeted panel consisting of 140 genes and additional regions of interest for copy number, as well as tiling of the Ig and MYC loci for detection of translocations. Samples were sequenced to a median depth of 452x using 2x75 bp paired end reads. Reads were aligned to hg19 and mutations called using Strelka and filtered with fpfilter. Translocations were called by Manta, and copy number determined by read depth ratio and loss of heterozygosity comparison with a patient matched non-tumor sample. DNA was obtained from either CD138+ cells from the bone marrow of multiple myeloma patients (tumor) or from stem cell harvests or peripheral blood cells from the same patient (control). 100 ng of DNA was fragmented, end-repaired, and adapters ligated using the HyperPlus kit (KAPA Biosystems). After PCR amplification the libraries were hybridized with probes against either a targeted panel consisting of 140 genes and chromosomal regions (Nimblegen) using SeqCap reagents (Nimblegen). Hybridized libraries underwent further amplification before being sequenced on a NextSeq500 (Illumina) using 75 bp paired end reads There are 450 (225 tumor and 225 germline) samples in this study. 263 are available as part of this dataset. The remaining 187 are available with dataset accession id EGAD00001004117.

  Study EGAS00001003223

• Immune induction strategies to enhance the sensitivity to PD-1 blockade in metastatic triple negative breast cancer: the TONIC-trial

  The response rate in metastatic triple negative breast cancer (TNBC) to PD-1 blockade is low, highlighting an urgent clinical need for strategies that render the TNBC tumor microenvironment (TME) more sensitive to PD-1 blockade. Immunomodulatory mechanisms have been proposed for both chemotherapy and irradiation, but it has not been established whether these therapies may improve efficacy of PD-1 blockade by favorably changing the TME. In the first stage of this adaptive, non-comparative phase II trial (TONIC-trial; NCT02499367), 67 patients with metastatic TNBC were randomized to nivolumab without induction or to one of four induction treatments, consisting of irradiation of one metastatic lesion (3x8 Gy) or a two- week low-dose regimen of cyclophosphamide, cisplatin or doxorubicin, all followed by nivolumab. In the overall cohort of 66 evaluable patients, the objective response rate (ORR; iRECIST16) was 20% (95% CI: 11%-31%), with two complete responses and 11 partial responses. The majority of clinical responses were observed on nivolumab in the cisplatin (ORR 23%; 3/13) and doxorubicin (ORR 35%; 6/17) cohorts. After doxorubicin and cisplatin induction, we detected an upregulation of immune-related genes, involved in PD-1/PD-L1, and T-cell cytotoxicity pathways. This was further supported by enrichment among upregulated genes related to inflammation, JAK-STAT and TNFα-signaling after doxorubicin. In addition, we observed a trend towards increased T-cell infiltration, measured using T-cell receptor (TCR) sequencing, after doxorubicin. Together, the clinical and translational data of this study testing different immune induction strategies prior to nivolumab suggest that short-term doxorubicin and cisplatin may induce a more favorable TME and increase the likelihood of response to PD-1 blockade in TNBC. These data warrant further confirmation in TNBC and exploration of induction treatments prior to PD-1 blockade in other cancer types.

  Study EGAS00001003535

• GDAP___Genome_Diversity_in_Africa_Project

  The Genomic Diversity in Africa Project (GDAP) started with the plan to develop a genomic resource from African populations, characterise genomic diversity and population history, and facilitate clinical studies in Africa. Currently, 25 individuals from 24 ethnolinguistic groups have been whole-genome sequenced at high depth totalling 585 individuals. An additional 41 individuals have been sequenced with 10X Genomics libraries. At this stage, the initial curation of this dataset has been finished and we are performing the analysis in coordination with our collaborators. The current state of the GDAP represents a very diverse panel of African populations that maximizes geographical and ethnic variation and represents a great starting point to achieve the aforementioned goals. However, southern sub-Saharan countries, Bantu speakers and hunter gatherer groups are currently underrepresented, despite being crucial to understand the evolutionary history of the continent. After extensive effort to collate studies documentation, we finally have the opportunity to sequence 600 new individuals from these groups, including countries as Gabon, Rwanda and Zambia, and address these deficiencies. We aim to proceed with the same strategy: to sequence at high depth 25 individuals with standard PCR free libraries, with 2 additional individuals with 10X Genomics Chromium libraries per ethnolinguistic group. The former allows a good representation of variants down to low frequency in any given population, and the latter allows accurate phasing and the analysis of structural variation. By including these new populations, we want to investigate three crucial questions in African history in addition to the initial objectives: the Bantu expansion, the evolutionary history of hunter gatherers and the transatlantic slave trade. Additionally, the expanded dataset will help us better discover the genetic variation present in Africa and characterize the African pangenome. 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 EGAS00001003602

• Repurposing azacitidine and carboplatin to prime immune checkpoint blockade-resistant melanoma for anti-PD-L1 re-challenge

  Purpose: Drug repurposing offers the opportunity for chemotherapy to be used to re-establish sensitivity to immune checkpoint blockade (ICB) therapy. Here we investigated the clinical and translational aspects of an early phase II study of azacitidine and carboplatin priming for anti-PDL1 immunotherapy (Avelumab) in patients with advanced ICB-resistant melanoma. Experimental Design: 20 participants with ICB resistant metastatic melanoma received 2 x 4-week cycles of azacitidine and carboplatin followed by ICB re-challenge with anti-PD-L1 avelumab. The primary objective was overall response rate after priming and ICB re-challenge. Secondary objectives were Clinical benefit rate (CBR), progression-free survival (PFS) and overall survival (OS). Translational correlation analysis of HLA-A and PD-L1 expression, RNA-seq and reduced representation bisulfite sequencing (RRBS) of biopsies at baseline, after priming and after 6 cycles of avelmuab was performed. Results: The ORR determined after azacitidine and carboplatin priming was 10% (2/20) with 2 partial responses (PR). The ORR determined after priming followed by 6 cycles of avelumab (week 22) was 10%, with 2/20 participants achieving iPR. The CBR for azacitidine and carboplatin priming was 65% (13/20) and after priming followed by 6 cycles of avelumab CBR was 35% (n = 7/20). The median PFS was 18.0 weeks (95% CI: 14.87 – 21.13 weeks) and the median OS was 47.86 weeks (95% CI: 9.67 – 86.06 weeks). Translational correlation analysis confirmed HLA-A generally increased after priming with azacitidine and carboplatin, particularly if it was absent at the start of treatment. Average methylation of CpGs across the HLA-A locus was decreased after priming and T-cells, in particular CD8+, showed the greatest increase in infiltration. Conclusions: Priming with azacitidine and carboplatin can induce disease stabilization and re-sensitisation to ICB for metastatic melanoma.

  Study EGAS00001006419

• Repurposing azacitidine and carboplatin to prime immune checkpoint blockade-resistant melanoma for anti-PD-L1 re-challenge

  Purpose: Drug repurposing offers the opportunity for chemotherapy to be used to re-establish sensitivity to immune checkpoint blockade (ICB) therapy. Here we investigated the clinical and translational aspects of an early phase II study of azacitidine and carboplatin priming for anti-PDL1 immunotherapy (Avelumab) in patients with advanced ICB-resistant melanoma. Experimental Design: 20 participants with ICB resistant metastatic melanoma received 2 x 4-week cycles of azacitidine and carboplatin followed by ICB re-challenge with anti-PD-L1 avelumab. The primary objective was overall response rate after priming and ICB re-challenge. Secondary objectives were Clinical benefit rate (CBR), progression-free survival (PFS) and overall survival (OS). Translational correlation analysis of HLA-A and PD-L1 expression, RNA-seq and reduced representation bisulfite sequencing (RRBS) of biopsies at baseline, after priming and after 6 cycles of avelmuab was performed. Results: The ORR determined after azacitidine and carboplatin priming was 10% (2/20) with 2 partial responses (PR). The ORR determined after priming followed by 6 cycles of avelumab (week 22) was 10%, with 2/20 participants achieving iPR. The CBR for azacitidine and carboplatin priming was 65% (13/20) and after priming followed by 6 cycles of avelumab CBR was 35% (n = 7/20). The median PFS was 18.0 weeks (95% CI: 14.87 – 21.13 weeks) and the median OS was 47.86 weeks (95% CI: 9.67 – 86.06 weeks). Translational correlation analysis confirmed HLA-A generally increased after priming with azacitidine and carboplatin, particularly if it was absent at the start of treatment. Average methylation of CpGs across the HLA-A locus was decreased after priming and T-cells, in particular CD8+, showed the greatest increase in infiltration. Conclusions: Priming with azacitidine and carboplatin can induce disease stabilization and re-sensitisation to ICB for metastatic melanoma.

  Study EGAS00001006420

• ctDNA quantification in Ewing sarcoma patients

  Background. Treatment stratification and response assessment in pediatric sarcomas has relied on imaging studies and surgical/histopathological evidence of vital tumor cells. Such studies and evidence collection processes often involve radiation and/or general anesthesia in children. Cell-free circulating tumor DNA (ctDNA) detection in blood plasma is one available method of so-called liquid biopsies that has been shown to correlate qualitatively and quantitatively to some extent with the existence of vital tumor cells in the body. Our clinical observational study focused on the utility and feasibility of ctDNA detection in pediatric Ewing sarcoma (EWS) as a marker of minimal residual disease (MRD). Patients and Methods. We performed whole genome sequencing to identify the exact breakpoints in tumors known to carry the EWS-FLI1 fusion gene. Patient-specific fusion breakpoints were tracked in peripheral blood plasma using digital droplet PCR (ddPCR) before, during, and after therapy in six consecutive children and young adults with EWS. Presence and levels of fusion breakpoints were correlated with clinical disease courses. The genetic changes in the tumors were defined, and inter- and intra-individual changes in the tumor DNA copy loads were observed over time. Results. We show that the detection of ctDNA in the peripheral blood of EWS patients i) is feasible in the clinical routine and ii) allows for the longitudinal real-time monitoring of MRD activity in children and young adults. Although changing ctDNA levels correlated well with clinical outcome within patients, between patients, a high variability was observed (inter-individually) for the six patients described here. Conclusions. ctDNA detection by ddPCR is a highly sensitive, specific, feasible, and highly accurate method that can be applied in EWS for follow-up assessments as an additional surrogate parameter for clinical MRD monitoring and, potentially, also for treatment stratification in the near future.

  Study EGAS00001006433

• Single Cell Dissection of the Tumour Microenvironment Reveals Dynamic Interplay Shaping the Tumour Immunity Continuum in Ovarian Cancer

  Diverse tumour immune phenotypes with distinct T cell infiltration patterns result in different responses to cancer immunotherapies. However, the key determinants and biology underpinning these phenotypes remain elusive. Here, we provide a high-resolution dissection of the tumour microenvironment representative of different tumour immune phenotypes. By performing single-cell RNA sequencing (scRNAseq) of the tumour, immune, and stromal compartments from 15 ovarian cancer patients, we characterized the diverse cellular and functional phenotypes, as well as dynamic interplay within and between these components. Infiltrated and excluded tumours markedly differed from desert tumours in several tumour cell intrinsic features, including high interferon response, antigen presentation and oxidative phosphorylation. Moreover, these tumours were characterized by tumour associated macrophage (TAM)-like myeloid subsets, while desert tumours were enriched in myeloid derived suppressive cell (MDSC)-like myeloid subsets. Infiltrated and excluded tumours differed in their T cell composition and fibroblast subsets. While dysfunctional CD8+ GZMB T cells and IL1-activated fibroblasts (IL1 CAFs) were associated with infiltrated tumours, pre-dysfunctional CD8+ GZMK T cells and TGFB-activated fibroblasts (TGFB CAFs) were enriched in excluded tumours. These findings validated in bulk transcriptomic profiles of 1071 ovarian tumours and using in situ hybridization. Furthermore, our study revealed chemokine receptor-ligand interactions within and across compartments as potential mechanisms mediating immune cell infiltration, exemplified by the tumour cell-T cell crosstalk via a CXC16-CXCR6 axis, stromal-immune cell crosstalk via CXCL12/14-CXCR4 signalling, and TAM-like macrophage-T cell crosstalk through CXCR3-CXCL9/10/11 signalling. Our in-depth characterization of the tumour environment and intercellular interactions provides novel insights into potential molecular mechanisms that shape the distinct biology of tumour immune phenotypes. Our findings may inform novel therapeutic strategies for potentially shifting tumours along the immunity continuum to improve clinical benefit from cancer immunotherapies.

  Study EGAS00001004935

• Somatic copy number alteration and fragmentation analysis in circulating tumor DNA for cancer screening and treatment monitoring in colorectal cancer patients

  Background: Analysis of circulating free DNA (cfDNA) is a promising tool for personalized management of colorectal cancer (CRC) patients. Untargeted cfDNA analysis using whole-genome sequencing (WGS) does not need a priori knowledge of the patient´s mutation profile. Methods: Here we established LIquid biopsy Fragmentation, Epigenetic signature and Copy Number Alteration analysis (LIFE-CNA) using WGS with ~6x coverage for detection of circulating tumor DNA (ctDNA) in CRC patients as a marker for CRC detection and monitoring. Results: We describe the analytical validity and a clinical proof-of-concept of LIFE-CNA using a total of 259 plasma samples collected from 50 patients with stage I-IV CRC and 61 healthy controls. To reliably distinguish CRC patients from healthy controls, we determined cutoffs for the detection of ctDNA based on global and regional cfDNA fragmentation patterns, transcriptionally active chromatin sites, and SCNAs. We further combined global and regional fragmentation pattern into a machine learning classifier to accurately predict ctDNA for cancer detection. By following individual patients throughout their course of disease, we show that LIFE-CNA enables the reliable prediction of response or resistance to treatment up to 3.5 months before commonly used CEA. Conclusion: In summary, we developed and validated a sensitive and cost-effective method for untargeted ctDNA detection at diagnosis as well as for treatment monitoring of all CRC patients based on genetic as well as non-genetic tumor-specific cfDNA features. Thus, once sensitivity and specificity have been externally validated, LIFE-CNA has the potential to be implemented into clinical practice. To the best of our knowledge, this is the first study to consider multiple genetic and non-genetic cfDNA features in combination with ML classifiers and to evaluate their potential in both cancer detection and treatment monitoring.

  Study EGAS00001006490

• The University of Hong Kong Colon Cancer GCV Study

  Background: Ganciclovir (GCV) is widely used in solid organ and hematopoietic stem cell transplant patients for prophylaxis and treatment of cytomegalovirus. It has long been considered a mutagen and carcinogen. However, the contribution of GCV to cancer incidence and other factors that influence its mutagenicity remains unknown. Methods: This retrospective cohort study analysed genomics data for 121,771 patients who had undergone targeted sequencing compiled by the Genomics Evidence Neoplasia Information Exchange (GENIE) or Foundation Medicine (FM). A statistical approach was developed to identify patients with GCV-associated mutational signature (GCVsig) from targeted sequenced data of tumour samples. Cell line exposure models were further used to quantify mutation burden and DNA damage caused by GCV and other antiviral and immunosuppressive drugs. Results: Mutational profiles from 22 of 121,771 patient samples in the GENIE and FM cohorts showed evidence of GCVsig. A diverse range of cancers was represented. All patients with detailed clinical history available had previously undergone solid organ transplantation and received GCV and mycophenolate treatment. RAS hotspot mutations associated with GCVsig¬¬ were present in 9 of the 22 samples, with all samples harbouring multiple GCV-associated protein-altering mutations in cancer driver genes. In vitro testing in cell lines showed that elevated DNA damage response and GCV¬sig are uniquely associated with GCV but not acyclovir, a structurally similar antiviral. Combination treatment of GCV with the immunosuppressant, mycophenolate mofetil (MMF), increased the misincorporation of GCV in genomic DNA and mutations attributed to GCVsig in cell lines and organoids. Conclusions: In summary, GCV can cause a diverse range of cancers. Its mutagenicity may be potentiated by other therapies, such as mycophenolate, commonly co-prescribed with GCV for post-transplant patients. Further investigation of the optimal use of these drugs could help reduce GCV-associated mutagenesis in post-transplant patients.

  Study EGAS00001006707

• High-throughput telomere length measurement at nucleotide resolution using the PacBio high fidelity sequencing platform, P49-P104 and run1 replicates

  Telomeres are specialized nucleoprotein structures at the ends of linear chromosomes that are essential for the stability as well as the complete replication of the human genome. The progressive shortening of steady-state telomere length in normal human somatic cells is a promising biomarker for age-associated diseases. However, there remain substantial challenges in quantifying telomere length in clinics and research laboratories due to the lack of accurate and high-throughput telomere length measurement method. Here, we describe a novel workflow to capture and analyze telomeres accurately at the single-nucleotide resolution. By using newly designed telobaits that are complementary to the single-stranded G-rich 3’ telomeric overhangs, we can specifically and efficiently capture full-length telomeres and subject them to single-molecule real-time (SMRT) sequencing. We show that the telomere length can be measured at nucleotide resolution using the PacBio high fidelity (HiFi) sequencing platform, which provides the necessary accuracy and throughput to uncover the steady-state telomere length distribution in human culture cell lines as well as clinical patient samples. This new method is a valuable tool not only in population-wide epidemiology studies but also for the longitudinal assessment of telomere attrition in individuals. Our results also reveal the extreme heterogeneity of telomeric variant sequences (TVSs) that is dispersed throughout the telomere repeat region. Importantly, the unique distribution pattern and sequence of TVSs in individuals may be used as a fingerprint for sample identification. In addition, the presence of TVSs disrupts the continuity of the canonical (5’-TTAGGG-3’)n telomere repeats, which could affect the binding of shelterin complexes at the chromosomal ends. Therefore, TVSs may act as a form of genetic polymorphisms that impairs telomere protection. Together with the absolute telomere length in each individual, TVSs may have profound implications in human aging and diseases.

  Study EGAS00001006595

• ctDNA to predict risk of progression and death after trifluridin/tipiracil therapy

  Background. Late-line treatment in metastatic colorectal cancer (mCRC) can improve the prognosis. However, not every patient has a benefit and may experience severe side effects. Thus, predictive/prognostic biomarkers are urgently needed. Despite many promising results, the role of circulating tumor DNA (ctDNA) analysis in trifluridin/tipiracil (FTD/TPI) treated mCRC patients is so far still elusive. Study design. This prospective non-interventional translational biomarker phase II study enrolled 30 metastatic CRC patients treated with FTD/TPI. The primary study aim was to identify a molecular biomarker to predict overall survival (OS) benefit upon FTD/TPI treatment. Methods. mCRC patients (60% male, 40% female) with FTD/TPI therapy were included. Molecular profiling in cell-free DNA from plasma was performed prior to treatment and on average 4, 8 and 12 weeks after FTD/TPI initiation. Tumor levels were assessed as the highest variant allele frequency (hVAF) and correlated with OS. Uni- and multivariate analyses for survival were performed with clinical variables. Multivariate analyses for survival and treatment efficacy were adjusted to age, gender and ECOG performance status. Results. Compared to previous data, our disease control rate of 30% was lower with a median OS of 8.1 months and a median PFS of 3.2 months. At baseline (BL) and the first two follow-ups (FU), a hVAF cut-off of 7.5% was most discriminative (BL HR 0.39 95% CI 0.16 - 0.98 P=0.0212, FU1 HR 0.30 95% CI 0.12 -0.77 P=0.0021, FU2 HR 0.23 95% CI 0.08-0.70 P=0.0004) whereas at the last follow-up the stratification was most pronounced at a 15% cut-off. An elevated hVAF trajectory over time was associated with a higher risk of death. Conclusion. ctDNA is detectable in a high proportion of mCRC patients. Changes in ctDNA levels are associated with survival among FTD/TPI therapy and may identify patients who benefit most from treatment in these late-stage disease.

  Study EGAS00001006883

• One-step generation of tumor models by base editor multiplexing in adult stem cell-derived organoids

  Recent optimization of CRISPR/Cas9-mediated genome engineering has resulted in the development of base editors that can efficiently mediate C>T and A>G transitions. Combining these genome engineering tools with human adult stem cell (ASC)-derived organoid technology holds promise for disease modeling. Here, we demonstrate the application of base editors for the generation of complex tumor models in human ASC-derived hepatocyte, endometrial and intestinal organoids. First, using conventional and evolved Cas9-variants, we show efficacy of both cytosine and adenine base editors and use them to model four hot-spot point mutations in CTNNB1 in hepatocyte organoids. Next, we apply C>T base editors in endometrial organoids to insert nonsense mutations in PTEN and demonstrate tumorigenicity even in the heterozygous state. Furthermore, we use cytosine base editors for simultaneous oncogene activation (PIK3CA) and tumor-suppressor inactivation (APC and TP53). To increase the flexibility of base editor multiplexing, we then combine SpCas9 and SaCas9 base editors for simultaneous C>T and A>G editing at individual target sites. Finally, we show the power of base editor multiplexing by modeling colorectal tumorigenesis in a single step by simultaneously transfecting sgRNA’s targeting four cancer genes. Seven clonal organoid lines and one bulk wild-type control sample were paired-end whole-genome sequenced using the Illumina Novaseq 6000 system. We sequenced four clonal intestinal organoid lines harbouring engineered TP53 and FBXW7 mutations as well as three lines targeted for oncogenic APC/TP53/PIK3CA/SMAD4 mutations. This WGS showed, as previously reported, a genome-wide increase in C>T mutations due to C>T base editor off-target activity, which is not enriched in predicted off-target regions based on the sgRNA sequences. Furthermore, we confirmed the absence of editing-induced driver mutations and lack of off-target mutational hotspots created by the genomic engineering.

  Study EGAS00001006886

• Separation of chronic myeloid leukemia stem cells from normal hematopoietic stem cells at single-cell resolution

  Chronic myeloid leukemia (CML) is a clonal disease, initiated and sustained by a subpopulation of CML stem cells (CML-SC) containing the BCR::ABL1 oncogene. Genetic markers that enrich for CML-SC have been identified, including CD26 being most specific. Genome sequencing analysis remain yet to be performed at single cell level in order to determine differences and heterogeneity in the genetic profile of CML-SC versus normal hematopoietic stem cells (HSC). Here we show that CD26 and CD33 separate CML-SC from normal HSC in freshly isolated bone marrow samples from newly diagnosed CML patients. DNA sequencing confirms that single CML-SC derived colonies are positive for BCR::ABL1 and harbor a recurrent set of patient-specific clonal somatic mutations across different CML-SC derived colonies. In contrast, marker-negative (CD26-CD33-) HSC from the same patients are negative for the BCR::ABL1 oncogene and do not carry recurrent, clonal mutation across different colonies. These results support the concept that CD26/CD33 expression separates genetically defined, clonal mutant CML-SC from normal HSC that carry individual sets of random mutations in CML patients. We provide whole exome DNA sequencing data for 23 clinical samples of CML-SC plus two buccal swipes derived normal samples. CML samples are comprised of 4 to 8 replicates from two patients, at diagnosis and after treatment. Single CML stem cells before treatment and single non-transformed HSC at remission were selected from bone marrow samples by FACS, according to genetic markers CD33+CD26+ at diagnosis and CD33+CD26-/CD33-CD26- at remission. WES libraries of colony forming assays derived CML-SC and HSC populations were prepared using Agilent SureSelect Human All Exon V6 kit and sequenced running 150 cycles (2x 75bp paired-end) on an Illumina NextSeq 500 platform. Following GATK best practices for somatic mutation calling, novel CML variants and random mutation events during remission phase were obtained.

  Study EGAS00001006904

• Molecular features of adenomas predicting metachronous colorectal cancer: a nested-case control study

  "Background & Aims: Current morphological features defining advanced adenomas (size ≥10 mm, high grade dysplasia or >25% villous component) cannot optimally distinguish individuals at high-risk or low-risk of metachronous colorectal cancer (me-CRC), which may result in suboptimal surveillance. Certain DNA copy number alterations are associated with adenoma-to-carcinoma progression (cancer-associated events; CAEs). We aimed to evaluate whether CAEs can better predict an individual’s risk of me-CRC than the morphological advanced adenoma features. Methods: In this nested case-control study, 529 individuals with a single adenoma at their first colonoscopy were selected from a Norwegian adenoma cohort. DNA copy number profiles were determined, by means of low-coverage whole genome sequencing. We assessed the prevalence of CAEs in advanced and non-advanced adenomas and its association with me-CRC. For the latter, cases (with CRC) were matched to controls (without CRC), with the same follow-up period as cases, as well as age and sex. Subsequently, odds ratios (OR) for the association of the presence of an advanced adenoma or a molecularly-defined high-risk adenoma with me-CRC were calculated. Results: Molecular high-risk features were observed in 85/267 (32%) of advanced adenomas and in 27/262 (10%) of non-advanced adenomas. Me-CRC was more strongly associated with advanced adenomas (OR 3.58; 95% CI 2.27-5.72) than with molecular high-risk adenomas (OR 1.90; 95% CI 1.14–3.16). After adjusting for other variables, only advanced adenoma remained significantly associated with me-CRC (OR 3.39; 95% CI 2.10–5.55). Conclusion: Molecularly-defined high-risk adenomas were associated with me-CRC, but the association of advanced adenoma with me-CRC was stronger. " Project goal: Evaluate the value of DNA copy numbers alterations detected in a baseline adenoma, in the prediction of risk of metachronous CRC

  Study EGAS00001007039

• Whole genome and whole exome sequencing of serial biopsies of relapsed/refractory diffuse large B-cell lymphoma.

  Purpose: Diffuse large B-cell lymphoma (DLBCL) is cured in more than 60% of patients, but outcomes remain poor for patients experiencing disease progression or relapse (refractory or relapsed DLBCL [rrDLBCL]), particularly if these events occur early. Although previous studies examining cohorts of rrDLBCL have identified features that are enriched at relapse, few have directly compared serial biopsies to uncover biological and evolutionary dynamics driving rrDLBCL. Here, we sought to confirm the relationship between relapse timing and outcomes after second-line (immuno)chemotherapy and determine the evolutionary dynamics that underpin that relationship. Patients and methods: Outcomes were examined in a population-based cohort of 221 patients with DLBCL who experienced progression/relapse after frontline treatment and were treated with second-line (immuno)chemotherapy with an intention-to-treat with autologous stem-cell transplantation (ASCT). Serial DLBCL biopsies from a partially overlapping cohort of 129 patients underwent molecular characterization, including whole-genome or whole-exome sequencing in 73 patients. Results: Outcomes to second-line therapy and ASCT are superior for late relapse (>2 years postdiagnosis) versus primary refractory (<9 months) or early relapse (9-24 months). Diagnostic and relapse biopsies were mostly concordant for cell-of-origin classification and genetics-based subgroup. Despite this concordance, the number of mutations exclusive to each biopsy increased with time since diagnosis, and late relapses shared few mutations with their diagnostic counterpart, demonstrating a branching evolution pattern. In patients with highly divergent tumors, many of the same genes acquired new mutations independently in each tumor, suggesting that the earliest mutations in a shared precursor cell constrain tumor evolution toward the same genetics-based subgroups at both diagnosis and relapse. Conclusion: These results suggest that late relapses commonly represent genetically distinct and chemotherapy-naïve disease and have implications for optimal patient management.

  Study EGAS00001007053

• Proteomic Analysis of Non-Muscle Invasive and Muscle Invasive Bladder Cancer Highlights Distinct Subgroups With Metabolic, Matrisomal, and Immune Hallmarks

  We present the proteome biology of bladder cancer based on a cohort comprising treatment-naïve non-muscle invasive (NMIBC, N=17) and treatment-naïve muscle-invasive bladder cancer (MIBC, N=51), as well as neoadjuvant-treated MIBC (N=11); using formalin-fixed, paraffin-embedded samples and data-independent acquisition yielding >8,000 protein identifications. The NMIBC proteome reveals two subgroups that represent cell adhesion or lipid metabolism biology. MIBC, compared to NMIBC, shows an altered metabolic proteome together with enrichment of extracellular matrix (ECM), immune response, and complement proteins; in line with three proteomic MIBC subgroups that represent ECM, metabolic, or immune signatures. The ECM and immune subgroups are enriched for markers of M2-like tumor-associated macrophages and dendritic cells, while markers for natural killer cells are exclusive for the ECM subgroup, and markers for cytotoxic T-cells are a hallmark of the immune subgroup. Infiltration of CD8+ PD1+ positive cytotoxic T cells correlated with PD-L1 expression of MIBC, as demonstrated by multiplex immunofluorescence analysis. The metabolic subgroup presents a comparably immune-depleted microenvironment. YAP1 signaling is particularly prominent in the ECM subgroup. The protein BSCL2 emerges as a prognostic marker for MIBC. We highlight increased endogenous proteolysis in MIBC alongside the upregulation of matrix metalloproteases exceeding the upregulation of their endogenous inhibitors. Genomic panel sequencing underlines prevalent FGFR3 mutations in NMIBC and prevalent TP53 mutations in MIBC, but a genomic underpinning of the MIBC subgroups remains elusive. Proteogenomic analysis yields protein-level evidence for more than 1,300 single amino acid variants. The proteomes of neoadjuvant-treated MIBC resembled treatment-naïve MIBC. Tumor-stroma interaction of MIBC was further investigated by proteomic analysis of patient-derived xenografts of MIBC. Our study adds to the proteome-level understanding of bladder cancer.

  Study EGAS00001007290

• Single cell TCRseq data fromTIL and non-TIL expanded TRACERx samples in the manuscript Bulk TCRseq from TRACERx samples from 'Subclonal immune evasion in non-small cell lung cancer'

  Immune checkpoint blockade in non-small cell lung cancer (NSCLC) rarely leads to complete tumour remission, suggesting that tumours harbour subpopulations of cancer cells that can escape immune predation. However, while tumours consist of multiple genetically distinct clones, the extent to which separate tumour subclones differ intrinsically in their functional capacity for immune evasion remains largely unexplored. This has been challenging to address due to an inability to isolate and propagate individual tumour subclones from human cancers. Here, we leverage the multi-region TRACERx lung cancer evolution study to generate a patient-derived platform that allows the functional analysis of immune escape at single clone resolution. We generated libraries of >20 separate NSCLC organoid lines per patient by isolating individual (clonal) organoids established from multiple spatially separated tumour regions. Each organoid subline was co-cultured with autologous tumour infiltrating lymphocytes (TIL) to evaluate their capacity to elicit a T cell response. Organoid lines derived from separate tumour regions, or from single clones within individual regions, differ vastly in their capacity to elicit a T cell response, showing the coexistence of intrinsically immune evading (‘cold) and non-evading (‘hot) clones in the same primary tumour. Organoid lines escaping T cell recognition represent genetically and transcriptionally unique subclones with a distinct evolutionary history, showing that evolutionary sculpting of subclonal diversity can have functional consequences for anti-tumour immunity. Notably, the identification of immune-evasive subclones would not have been possible using computational approaches alone, highlighting the unique value of functional co-culture platforms. This patient-derived functional co-culture platform reveals that tumour evolution can give rise to distinct cancer clones with intrinsic differences in immune evasion capacity and provides an approach to prospectively identify, isolate and characterise immune evading subclones from patients with cancer.

  Study EGAS00001008161

• Dataset for Sarcoma RNA sequencing linked from study EGAS00001004813

  Sequencing data of 410 sarcoma tumor runs, which were uploaded to EGAS00001004813. The sequencing was always paired.

  Dataset EGAD00001010276

• Substantial somatic genomic variation and selection for BCOR mutations in human induced pluripotent stem cells

  Human Induced Pluripotent Stem Cells (hiPSC) are an established patient-specific model system where opportunities are emerging for cell-based therapies. We compared and contrasted hiPSCs derived from different tissues, skin and blood, in the same individual. We show extensive single-nucleotide mutagenesis in all hiPSC lines, although fibroblast-derived hiPSCs (F-hiPSCs) are particularly heavily mutagenized by ultraviolet (UV)-related damage. We utilized genome sequencing data on 454 F-hiPSCs and 44 blood-derived hiPSCs (B-hiPSCs) to gain further insights. Across 324 whole genome sequenced (WGS) F-hiPSCs derived by the Human Induced Pluripotent Stem Cell Initiative (HipSci), UV-related damage is present in ~72% of cell lines, sometimes causing substantial mutagenesis (range 0.25-15 per Mb). Furthermore, we find remarkable genomic heterogeneity between independent F-hiPSC clones derived from the same reprogramming process in the same donor, due to oligoclonal populations within fibroblasts. Combining WGS and exome-sequencing data of 452 HipSci F-hiPSCs, we identify 272 predicted pathogenic mutations in cancer-related genes, of which 21 genes were hit recurrently three or more times, involving 77 (17%) lines. Notably, 151 of 272 mutations were present in starting fibroblast populations suggesting that more than half of putative driver events in F-hiPSCs were acquired in vivo. In contrast, B-hiPSCs reprogrammed from erythroblasts show lower levels of genome-wide mutations (range 0.28-1.4 per Mb), no UV damage, but a strikingly high prevalence of acquired BCOR mutations in ~57% of lines, indicative of strong selection pressure. All hiPSCs had otherwise stable, diploid genomes on karyotypic pre-screening, highlighting how copy-number-based approaches do not have the required resolution to detect widespread nucleotide mutagenesis. This work strongly suggests that models for cell-based therapies require detailed nucleotide-resolution characterization prior to clinical application.

  Dataset EGAD00001007029

• RNAseq dataset for MALT1 in MCL

  Manuscript Title: Co-targeting of BTK and MALT1 overcomes resistance to BTK inhibitors in mantle cell lymphoma Journal: Journal of Clinical Investigation Authors Vivian Changying Jiang1, Yang Liu1, Junwei Lian1, Shengjian Huang1, Alexa Jordan1, Qingsong Cai1, Fangfang Yan3, Joseph Mitchell McIntosh1, Yijing Li1, Yuxuan Che1, Zhihong Chen1, Jovanny Vargas1, Maria Badillo1, JohnNelson Bigcal1, Heng-Huan Lee1, Wei Wang1, Yixin Yao1, Lei Nie1, Christopher Flowers1, and Michael Wang1, 2* Abstract Bruton’s tyrosine kinase (BTK) is a proven target in mantle cell lymphoma (MCL), an aggressive subtype of non-Hodgkin lymphoma. However, resistance to BTK inhibitors is a major clinical challenge. We here report that MALT1 is one of the top overexpressed genes in ibrutinib-resistant MCL cells, while expression of CARD11, which is upstream of MALT1, is decreased. MALT1 genetic knockout or inhibition produced dramatic defects in MCL cell growth regardless of ibrutinib sensitivity. Conversely, CARD11 knockout cells showed anti-tumor effects only in ibrutinib-sensitive cells, suggesting that MALT1 overexpression could drive ibrutinib resistance via bypassing BTK-CARD11 signaling. Additionally, BTK knockdown and MALT1 knockout markedly impaired MCL tumor migration and dissemination, and MALT1 pharmacological inhibition decreased MCL cell viability, adhesion, and migration by suppressing NF-κB, PI3K-ATK-mTOR, and integrin signaling. Importantly, co-targeting MALT1 with safimaltib and BTK with pirtobrutinib induced potent anti-MCL activity in ibrutinib-resistant MCL cell lines and patient-derived xenografts. Therefore, we conclude that MALT1 overexpression associates with resistance to BTK inhibitors in MCL, targeting abnormal MALT1 activity could be a promising therapeutic strategy to overcome BTK inhibitor resistance, and co-targeting of MALT1 and BTK should improve MCL treatment efficacy and durability as well as patient outcomes. Dataset description: The bulk RNA-seq dataset was generated for the cell lines below and used for two major purposes: 1. DEG analysis and GSEA analysis comparing IBN-R and IBN-S cells 2. DEG analysis and GSEA analysis comparing MCL cells with/without MI-2 treatment. sample Cell MI-2 Ibrutinib (IBN) Venetoclax (VEN) Used for IBN-R vs IBN-S comparison Used for MI-2 vs untreated (DMSO) H9 Granta519 - R S yes H21 Granta519 - R S yes H33 Granta519 - R S yes H10 Granta519-VEN-R - R R yes H22 Granta519-VEN-R - R R yes H34 Granta519-VEN-R - R R yes H3 JeKo BTK KD_1 - R R yes yes H15 JeKo BTK KD_1 - R R yes yes H27 JeKo BTK KD_1 - R R yes yes H5 JeKo BTK KD_2 - R R yes yes H17 JeKo BTK KD_2 - R R yes yes H29 JeKo BTK KD_2 - R R yes yes H1 JeKo-1 - S R yes yes H13 JeKo-1 - S R yes yes H25 JeKo-1 - S R yes yes H7 Mino - S S yes H19 Mino - S S yes H31 Mino - S S yes H8 Mino-VEN-R - S R yes H20 Mino-VEN-R - S R yes H32 Mino-VEN-R - S R yes H11 Rec-1 - S S yes H23 Rec-1 - S S yes H12 Rec-VEN-R - S S yes H24 Rec-VEN-R - S R yes H36 Rec-VEN-R - S R yes H35 Rec-1 -- S R yes H4 JeKo BTK KD_1 + MI-2 + yes H16 JeKo BTK KD_1 + MI-2 + yes H28 JeKo BTK KD_1 + MI-2 + yes H6 JeKo BTK KD_2 + MI-2 + yes H18 JeKo BTK KD_2 + MI-2 + yes H30 JeKo BTK KD_2 + MI-2 + yes H2 JeKo-1 + MI-2 + yes H14 JeKo-1 + MI-2 + yes H26 JeKo-1 + MI-2 + yes

  Dataset EGAD00001009771

• Solve-RD_RND_cohort-1_DF1+2_V1

  Solve-RD data submitted to the ERN-RND cohort for re-analysis (Data freeze 1+2) v1

  Dataset EGAD00001009769

• Procardis study on novel susceptibility genes for coronary artery disease (CAD)

  An increased level of Lp(a) lipoprotein has been identified as a risk factor for coronary artery disease that is highly heritable. The genetic determinants of the Lp(a) lipoprotein level and their relevance for the risk of coronary disease are incompletely understood. METHODS: We used a novel gene chip containing 48,742 single-nucleotide polymorphisms (SNPs) in 2100 candidate genes to test for associations in 3145 case subjects with coronary disease and 3352 control subjects. Replication was tested in three independent populations involving 4846 additional case subjects with coronary disease and 4594 control subjects. RESULTS: Three chromosomal regions (6q26-27, 9p21, and 1p13) were strongly associated with the risk of coronary disease. The LPA locus on 6q26-27 encoding Lp(a) lipoprotein had the strongest association. We identified a common variant (rs10455872) at the LPA locus with an odds ratio for coronary disease of 1.70 (95% confidence interval [CI], 1.49 to 1.95) and another independent variant (rs3798220) with an odds ratio of 1.92 (95% CI, 1.48 to 2.49). Both variants were strongly associated with an increased level of Lp(a) lipoprotein, a reduced copy number in LPA (which determines the number of kringle IV-type 2 repeats), and a small Lp(a) lipoprotein size. Replication studies confirmed the effects of both variants on the Lp(a) lipoprotein level and the risk of coronary disease. A meta-analysis showed that with a genotype score involving both LPA SNPs, the odds ratios for coronary disease were 1.51 (95% CI, 1.38 to 1.66) for one variant and 2.57 (95% CI, 1.80 to 3.67) for two or more variants. After adjustment for the Lp(a) lipoprotein level, the association between the LPA genotype score and the risk of coronary disease was abolished. CONCLUSIONS: We identified two LPA variants that were strongly associated with both an increased level of Lp(a) lipoprotein and an increased risk of coronary disease. Our findings provide support for a causal role of Lp(a) lipoprotein in coronary disease.

  Study EGAS00000000055

• Wellcome Trust Sanger Institute

  Wellcome Trust Sanger Institute

  Dac EGAC00001000205

• Knoll et al.: The life-saving benefit of dexamethasone in severe COVID-19 is linked to a reversal of monocyte dysregulation

  Following the demonstration of its life-saving effect in clinical trials, dexamethasone quickly became standard-of-care in the treatment of severe COVID-19. Beneficial effects were reported for patients requiring supplemental oxygen or invasive mechanical ventilation. Yet, a substantial proportion of patients still progressed to a critical condition or succumbed to the disease, despite timely initiation of glucocorticoid treatment. No molecular or cellular correlate of this beneficial treatment response has been defined. Here, we identify distinct cellular and molecular changes in circulating immune cells in patients with COVID-19 elicited in response to dexamethasone treatment. The most profound transcriptional changes in response to dexamethasone treatment were noted in monocytes and in B cells. In monocytes, we observed a reversal of hallmark signatures previously associated with COVID-19 severity, and the induction of a specific monocyte substate enriched in glucocorticoid response gene expression which we here refer to as dexamethasone response cluster. Molecular responses to dexamethasone treatment were directly linked to clinical outcome, since the reversal of pathogenic signatures and induction of glucocorticoid response genes were enriched in patients who survived the disease in response to dexamethasone treatment. Cellular responsiveness of circulating monocytes was thus identified as a correlate of clinical response to dexamethasone treatment. Changes in monocyte transcriptomes could also be linked to epigenetic alterations highlighting early differences in monocytes of dexamethasone responders and non-responders. Further, monocyte single-cell transcriptome-derived signatures were enriched in whole blood transcriptomes from patients with fatal outcome in two independent cohorts, highlighting the potential clinical value for early identification of non-responders who are refractory to dexamethasone treatment. Overall, our findings indicate that the life-saving effects of dexamethasone in COVID-19 patients are linked to a specific immunomodulatory effect based on the reversion of monocyte dysregulation. Our study elucidates the potential mechanism of action of one of the most efficient and cost-effective drugs for the treatment of COVID-19. It further highlights the potential of single-cell omics for monitoring target engagement of immunomodulatory drugs and for the stratification of patients for precision medicine approaches.

  Study EGAS00001007461

• Neoadjuvant atezolizumab plus chemotherapy in gastric and gastroesophageal junction adenocarcinoma: the phase 2 PANDA trial

  Gastric and gastroesophageal junction (G/GEJ) cancers carry a poor prognosis, and despite recent advancements most patients die of their disease. While immune checkpoint blockade has become part of the standard of care for patients with metastatic G/GEJ cancers, its efficacy and impact on the tumor microenvironment in the early disease setting remain largely unknown. We hypothesized a higher efficacy of neoadjuvant immunotherapy plus chemotherapy in patients with non-metastatic G/GEJ cancer. In the phase 2 PANDA trial, patients with previously untreated resectable G/GEJ tumors (n=21) received neoadjuvant treatment with one cycle of atezolizumab monotherapy followed by 4 cycles of atezolizumab plus docetaxel, oxaliplatin and capecitabine. Treatment was well-tolerated, and all patients underwent resection without treatment-related delays, meeting the primary endpoint of safety and feasibility. Tissue was obtained at multiple time points, allowing analysis of the effects of single-agent PD-L1 blockade and the subsequent combination with chemotherapy on the tumor microenvironment. A total of 21 patients were included, of whom 20 patients underwent surgery and were evaluable for the secondary pathologic response and survival endpoints, while 19 were evaluable for exploratory translational analyses. A major pathologic response (MPR, ≤10% residual viable tumor) was observed in 14/20 (70%; 95% CI 46-88%) patients, including 9 (45%; 95% CI 23-68%) pathologic complete responses. At a median follow-up of 47 months, 13/14 responders were alive and disease-free, while 5/6 nonresponders had died due to disease recurrence. Notably, baseline PD-1+CD8+ T cell infiltration was significantly higher in responders than in nonresponders, and comparison of the tumor microenvironment alterations following anti-PD-L1 monotherapy versus the subsequent combination with chemotherapy revealed that increased immune activity was achieved upon single-agent PD-1/L1 axis blockade. Based on these data, anti-PD-L1 plus chemotherapy warrants further exploration in patients with non-metastatic G/GEJ cancer, and these results call for validation in a larger cohort.

  Study EGAS50000000168