22928 results for "ega"

in 54.80 milliseconds.

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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