23014 results for "ega"

in 97.13 milliseconds.

• B cell receptor silencing reveals origin and dependencies of high-grade B cell lymphomas with MYC and BCL2 rearrangements

  High-grade B cell lymphomas with MYC and BCL2 rearrangements (HGBCL-DH-BCL2) represent highly aggressive Germinal center (GC) derived B cell malignancies with dismal prognosis and limited therapeutic options. In the present study, we investigated the expression of BCR complex components, including immunoglobulin (IG) heavy (H) and light (L) chains, and CD79A and CD79B signaling subunits in HGBCL-DH-BCL2, scoring 65% of cases with undetectable IGH immunoreactivity (IGHUND). Transcriptomics (bulk RNA-seq) and genetic (whole exome/genome sequencing) analyses indicated that IGHUND and IGH+ HGBCL-DH-BCL2 lymphomas represent distinct biological entities. Specifically, IGHUND HGBCL-DH-BCL2 exhibit a GC DZ-like molecular program, and select subset-specific gene mutations pointing to distinct evolutionary trajectories respect to the IGH+ counterparts. IG profiling revealed that IGHUND HGBCL-DH-BCL2 preserved potentially productive, hypermutated IGH variable (V) gene rearrangements, indicating a selective advantage conferred to the malignant B cells by residual IGH molecules. Accordingly, disruption of the productively rearranged IGHG1 gene in BCRnull HGBCL-DH-BCL2 cell line model confirmed the dependency on IGH expression for optimal fitness. Using transcriptomics and in situ RNA-scope analysis, we reported that IGHUND HGBCL-DH-BCL2 originated from IG-class switched GC B cells, while IGH-expressing tumors were selected for continuous IGM expression. Moreover, while IGH+ HGBCL-DH-BCL2 favor IGM/IG-Kappa expression, IGHUND counterparts complete IG-Lambda light chain rearrangements. Notably, BCR silencing, caused by accelerated IGH turnover and reduced IGH expression, precedes HGBCL-DH-BCL2 onset, inducing RAG1/2-dependent IG light chain editing and eventually facilitating t(8;22)/IGL::MYC translocations. Finally, IGH silencing observed in HGBCL-DH-BCL2 cell models, which leads to intracellular retention/silencing of CD79B, results in reduced sensitivity to the CD79B-targeting antibody-drug conjugate Polatuzumab-vedotin. Collectively, HGBCL-DH-BCL2 originates from isotype-switched t(14;18)+ GC B cells that undergo IG light chain editing fuelling intra-clonal diversification, BCR silencing, and the acquisition of t(8;22) translocation, while preserving IGH dependency.

  Study EGAS50000001047

• Expanding the SPTBN1 Spectrum: A Case of Neurodevelopmental Disorder with GI and Musculoskeletal Involvement

  Background: Mutations in the SPTBN1 gene, which encodes the βII-spectrin protein, are associated with a recently described neurodevelopmental syndrome characterized by developmental delay, intellectual disability, autistic traits, seizures, and dysmorphic features. While the phenotypic spectrum continues to evolve, gastrointestinal (GI) involvement and vertical transmission remain rarely reported features. Case Presentation: We describe a 13-year-old male from Georgia with a heterozygous missense mutation in SPTBN1 (c.3799C>G, p.Leu1267Val), inherited from his clinically unaffected mother. The patient presented with cognitive and motor delays, autistic traits, behavioral challenges, feeding difficulties, growth delay, musculoskeletal anomalies (pectus excavatum, pes planus), and persistent GI symptoms including abdominal pain and gallbladder sludge. Despite therapy with ursodeoxycholic acid and nutritional counseling, the gallbladder abnormalities persisted over a 12-month follow-up. Feeding therapy could not be initiated due to access barriers, although parental efforts led to modest dietary improvement. Upon recent visit pes planus was also confirmed. Conclusion: This case contributes to the growing body of literature on SPTBN1-related neurodevelopmental disorders, highlighting the potential for a broader phenotypic spectrum than currently recognized. In addition to the well-documented behavioral, cognitive, and musculoskeletal features, our patient presented with notable gastrointestinal involvement—specifically, gallbladder sludge and feeding disorder—which are not commonly reported in association with SPTBN1 mutations. Although causality cannot be established from a single case, the observed GI manifestations may be indirectly related to underlying neurobehavioral features such as ASD and associated eating disorders. The identification of a vertically transmitted missense variant (c.3799C>G, p.Leu1267Val) adds to the limited but growing evidence of inherited SPTBN1 mutations, expanding our understanding of transmission patterns. Given the novelty of the findings, we emphasize the need for multidisciplinary management, including feeding and behavioral therapy, as well as longitudinal follow-up to assess developmental, nutritional, and hepatobiliary outcomes. Future studies and case series are warranted to determine whether pes planus, pectus excavatum, gastrointestinal symptoms, including biliary tract involvement, represent an underrecognized component of the SPTBN1 phenotypic spectrum.

  Study EGAS50000001213

• German early-onset prostate cancer cohort of the Pan-Prostate Cancer Genome (PPCG) project

  Prostate cancer is the most frequent malignant tumor in males and the second most frequent cause of cancer-related death. Currently, in Germany, more than 60,000 prostate cancers are diagnosed every year. Although most of these patients are treated in a curative attempt, more than 10,000 German men die from prostate cancer annually. Owing to the demographic changes of our society, a further doubling of prostate cancer incidences during the next 20 years is expected. Prostate cancer is generally considered a tumor of elderly men. However, a fraction of prostate cancers are diagnosed at the age of 55 years or less. For several reasons, these “early onset prostate cancers” may represent a key entity for the understanding of prostate cancer biology. First, it is likely that early onset prostate cancers represent a distinct molecular subgroup of prostate cancer (PCa), potentially characterized by relatively small numbers of genetic changes, some of which may be particularly strong driver mutations for PCa development. Second, a fraction of prostate cancers in young individuals could represent classical prostate cancers that are detected at a very early stage and might therefore have accumulated molecular changes/mutations occurring that are most instrumental for prostate cancer early detection. Third, PCa with hereditary backgrounds are likely to accumulate in the age group below 55 years. A comparison with other sample sets (e.g. from other ICGC consortia), a systematic genomic analysis of young men with PCa could therefore lead to the detection of mechanisms for hereditary PCa. Fourth, a better understanding of these tumors is particularly relevant as finding optimal treatment regimens is most critical in young cancer patients. We are analyzing the entire genomic DNA sequences of at least 200 PCas (and matched non-tumorous DNA) of young men (≤ 55y), to at least 30 fold coverage, and integrate single-nucleotide variants and genomic structural variations with differential methylation, mRNA and miRNA expression data.Tissues were collected at the University Medical Center Hamburg-Eppendorf. Sequencing was performed at DKFZ and NCT (Heidelberg), EMBL (Heidelberg), and MPIMG (Berlin). Data management and bioinformatic data analyses is being conducted at DKFZ.

  Study EGAS00001003373

• Ancient Yakuts genomes and metagenomes

  Read data and bam alignments for 124 ancient yakuts

  Dataset EGAD50000001903

• METABRIC

  The elucidation of breast cancer subgroups and their molecular drivers requires integrated views of the genome and transcriptome from representative numbers of patients. We present an integrated analysis of copy number and gene expression in a discovery and validation set of 997 and 995 primary breast tumours, respectively, with long-term clinical follow-up. Inherited variants (CNVs, SNPs) and acquired somatic copy number aberrations (CNAs) were associated with expression in 40% of genes, although the landscape was dominated by cis and trans-acting CNAs. By delineating expression outlier genes driven in cis by CNAs, we identified putative cancer genes, including deletions in PPP2R2A, MTAP, and MAP2K4. Unsupervised analysis of paired DNA/RNA profiles revealed novel subgroups with distinct clinical outcomes, which reproduced in the validation cohort. These include a high-risk, ER-positive 11q13/14 cis-acting subgroup and a favourable prognosis subgroup devoid of CNAs. Trans-acting aberration 0152hotspots were found to modulate subgroup-specific gene networks, including a TCR deletion-mediated adaptive immune response in the 0152CNA-devoid sub-group and a Basal-specific chromosome 5 deletion-driven mitotic network. Our results provide a novel molecular stratification of the breast cancer population, derived from the impact of somatic copy number aberrations on the transcriptome.

  Study EGAS00000000083

• Lindsey E. Jones et al. Patient-derived cells from recurrent tumors that model the evolution of IDH-mutant glioma. Neuro-Oncology Advances, 2020. We established robust, versatile, and well-characterized Patient Derived Cells (PDC) of an IDH1-mutant astrocytoma and an IDH1-mutant oligodendroglioma that represents defined evolutionary stages of chemotherapy-induced hypermutation. The PDCs retain tumor subtype defining features over time as well as classic hallmarks of cancer, including anchorage independent growth and cell immortality. The integrated phylogenies composed of PDCs,single-cell-derived PDCs, patient-derived xenografts, and corresponding spatiotemporal tumor tissue samples also provide new insight into the clonality, evolutionary pattern, and immense mutational load of hypermutated IDH mutant gliomas.

  IDH-mutant lower-grade gliomas can undergo malignant progression via temozolomide-driven hypermutation. Patient-derived cells (PDC) that model the genetically distinct hypermutated (HM) tumor subgroup are generally lacking, and few if any human brain tumor cell models are from defined evolutionary time points. Here, we characterize multiple PDC derived from independent surgical specimens of IDH1-mutant recurrences, including an ATRX and TP53-mutant astrocytoma and a 1p/19q co-deleted and TERT promoter-mutant oligodendroglioma. We determined the evolutionary time points represented by each PDC using exome sequencing and phylogenetic reconstruction, comparing the PDC and single cell clones of the PDC (scPDC) to multiple spatiotemporal tumor tissue samples, PDC-derived xenografts (PDX) and patient-matched blood. The tumor samples exhibited TMZ-induced mutagenesis and a branching pattern of evolution. We found clear evidence of two fully independent founder HM clones in the tumor tissue that are faithfully represented by independent PDC. The PDC, scPDC and PDX also shared the mutagenesis signature and represent the mid and later evolutionary time points of their corresponding tumors. The PDC maintained the tumor subtype-defining features over many passages, including heterozygosity of the IDH1 R132H mutation, production of 2-hydroxyglutarate (2-HG), and subtype-specific telomere maintenance mechanisms. The PDC from both tumors exhibited anchorage-independent growth in soft agar. The oligodendroglioma PDC formed infiltrative intracranial tumors with characteristic oligodendroglioma histology, initially with a long period to tumor formation. We conclude that the PDC, scPDC and PDX faithfully model the heterogeneous clonal origins of the corresponding tumor tissue. The multilevel analysis also provides new insight into the intratumoral heterogeneity and vast mutational load of HM glioma. The PDC from multiple evolutionary time points presented in the context of full clinical timelines may be useful to model evolution and intratumoral heterogeneity, important sources of therapeutic failure.

  Study EGAS00001003992

• OpACIN-neo – A Multicenter Phase 2 Study to identify the Optimal neo-Adjuvant Combination scheme of Ipilimumab and Nivolumab – Whole exome sequencing

  The outcome of patients with macroscopic stage III melanoma is poor. Neoadjuvant treatment with ipilimumab plus nivolumab at the standard dosing schedule induced pathological responses in a high proportion of patients in two small independent early-phase trials. However, toxicity of the standard ipilimumab plus nivolumab dosing schedule was high, preventing its broader clinical use. The aim of the OpACIN-neo study (NCT02977052) to identify a dosing schedule of ipilimumab plus nivolumab that is less toxic but equally effective. 86 patients with macroscopic stage III melanoma were randomised to one of three dosing schedules: arm A, two cycles of ipilimumab 3 mg/kg plus nivolumab 1 mg/kg once every 3 weeks; arm B, two cycles of ipilimumab 1 mg/kg plus nivolumab 3 mg/kg once every 3 weeks; or arm C, two cycles of ipilimumab 3 mg/kg once every 3 weeks directly followed by two cycles of nivolumab 3mg/kg once every 2 weeks. Within the first 12 weeks, grade 3–4 immune-related adverse events were observed in 12 (40%) of 30 patients in group A, six (20%) of 30 in group B, and 13 (50%) of 26 in group C. Pathological responses occurred in 24 (80%) patients in group A, 23 (77%) in group B, and 17 (65%) in group C. The 2-year estimated relapse-free survival was 84% for all patients, 97% for patients achieving a pathologic response and 36% for nonresponders (p<0.001). High tumor mutational burden (TMB) and high interferon-gamma-related gene expression signature score (IFN-γ score) were associated with pathologic response and low risk of relapse; pathologic response rate (pRR) was 100% in patients with high IFN-γ score/high TMB, patients with high IFN-γ score/low TMB or low IFN-γ score/high TMB had pRRs of 91% and 88%; while patients with low IFN-γ score/low TMB had a pRR of only 39%. These data demonstrate long-term benefit in patients with a pathologic response and show the predictive potential of TMB and IFN-γ score.

  Study EGAS00001004832

• Automated machine-learning approach for next generation profiling of sequence alterations, mutation burden, microsatellite instability, and structural variants in human cancers

  Sequence and structural alterations together with tumor mutation burden (TMB) and microsatellite instability (MSI) have been identified as biomarkers for the determination of response to targeted and immune checkpoint inhibitor therapies. However, widespread clinical adoption of these biomarkers has historically been limited due to barriers such as evidence of clinical utility and reimbursement. We have developed 2.2 Mb targeted NGS system and an automated machine-learning analysis approach (PGDx elio™ tissue complete, ETC) that has been FDA cleared for examination of 500+ cancer-related genes and 68 mononucleotide repeats for identification of sequence and structural alterations, TMB, and MSI in solid cancers in a clinical setting. We designed and trained this approach using sequence data from 4,174 cancers and >124,000 in silico alterations and evaluated the methodology in >2,550 tumor or non-cancerous normal samples. Independent analyses of ETC sequence changes in 440 formalin fixed paraffin embedded (FFPE) tumor or cell line samples using MSK-IMPACT™, FoundationOne®, and ddPCR revealed a positive percent agreement (PPA) >97% with high sensitivity as low as 3% mutant allele fraction. We observed high concordance between panel-wide and whole-exome TMB for 307 pan-cancer FFPE tumors (Pearson r=0.95, p < 0.0001) using samples with ≥20% tumor cellularity. Comparison of the mutation context and repeat-based MSI approach in ETC with a multiplex MSI PCR assay in 223 samples revealed a PPA of 99% and negative predictive agreement (NPA) >99%. We confirmed the accuracy and precision of TMB and MSI measurements across three independent laboratories (CV of <5% and average PPA >99%, respectively). Finally, evaluation of amplifications and translocations against DNA and RNA-based approaches exhibited >98% NPA and PPA of 86% and 82% respectively. These results demonstrate high analytical performance for determination of sequence and structural changes, TMB, and MSI using a targeted NGS panel and provide a scalable and facile approach for evaluating these biomarkers in a clinical laboratory.

  Study EGAS00001005556

• OpACIN-neo – A Multicenter Phase 2 Study to identify the Optimal neo-Adjuvant Combination scheme of Ipilimumab and Nivolumab – RNA sequencing

  The outcome of patients with macroscopic stage III melanoma is poor. Neoadjuvant treatment with ipilimumab plus nivolumab at the standard dosing schedule induced pathological responses in a high proportion of patients in two small independent early-phase trials. However, toxicity of the standard ipilimumab plus nivolumab dosing schedule was high, preventing its broader clinical use. The aim of the OpACIN-neo study (NCT02977052) to identify a dosing schedule of ipilimumab plus nivolumab that is less toxic but equally effective. 86 patients with macroscopic stage III melanoma were randomised to one of three dosing schedules: arm A, two cycles of ipilimumab 3 mg/kg plus nivolumab 1 mg/kg once every 3 weeks; arm B, two cycles of ipilimumab 1 mg/kg plus nivolumab 3 mg/kg once every 3 weeks; or arm C, two cycles of ipilimumab 3 mg/kg once every 3 weeks directly followed by two cycles of nivolumab 3mg/kg once every 2 weeks. Within the first 12 weeks, grade 3–4 immune-related adverse events were observed in 12 (40%) of 30 patients in group A, six (20%) of 30 in group B, and 13 (50%) of 26 in group C. Pathological responses occurred in 24 (80%) patients in group A, 23 (77%) in group B, and 17 (65%) in group C. The 2-year estimated relapse-free survival was 84% for all patients, 97% for patients achieving a pathologic response and 36% for nonresponders (p<0.001). High tumor mutational burden (TMB) and high interferon-gamma-related gene expression signature score (IFN-γ score) were associated with pathologic response and low risk of relapse; pathologic response rate (pRR) was 100% in patients with high IFN-γ score/high TMB, patients with high IFN-γ score/low TMB or low IFN-γ score/high TMB had pRRs of 91% and 88%; while patients with low IFN-γ score/low TMB had a pRR of only 39%. These data demonstrate long-term benefit in patients with a pathologic response and show the predictive potential of TMB and IFN-γ score.

  Study EGAS00001004833

• ACT-Discover: identifying karyotype heterogeneity in pancreatic cancer evolution using ctDNA

  Background Liquid biopsies, and the dynamic tracking of somatic mutations within circulating tumour DNA (ctDNA) can provide insight into the dynamics of cancer evolution and the intra-tumour heterogeneity that fuels treatment resistance. However, identifying and tracking dynamic changes in somatic copy-number alterations (SCNAs), which have been associated with poor outcome and metastasis, using ctDNA is challenging. Pancreatic adenocarcinoma is a disease which has been considered to harbour early punctuated events in its evolution, leading to an early fitness peak, with minimal further subclonal evolution. Methods To interrogate the role of SCNAs in pancreatic adenocarcinoma cancer evolution, we applied whole-exome sequencing of 55 longitudinal cell-free DNA (cfDNA) samples taken from 24 patients (including 8 from whom a patient-derived xenograft (PDX) was derived) with metastatic disease prospectively recruited into a clinical trial. We developed a method, Aneuploidy in Circulating Tumour DNA (ACT-Discover), that leverages haplotype phasing of paired tumour biopsies or PDXs to identify SCNAs in cfDNA with greater sensitivity. Results SCNAs were observed within 28 of 47 evaluable cfDNA samples. Of these events, 30% could only be identified by harnessing the haplotype-aware approach leveraged in ACT-Discover. The exceptional purity of PDX tumours enabled near-complete phasing of genomic regions in allelic imbalance, highlighting an important auxiliary function of PDXs. Finally, although the classical model of pancreatic cancer evolution emphasises the importance of early, homogenous somatic events as a key requirement for cancer development, ACT-Discover identified substantial heterogeneity of SCNAs, including parallel focal and arm-level events, affecting different parental alleles within individual tumours. Indeed, ongoing acquisition of SCNAs was identified within tumours throughout the disease course, including within an untreated metastatic tumour. Conclusions This work demonstrates the power of haplotype phasing to study genomic variation in cfDNA samples and reveals undiscovered intra-tumour heterogeneity with important scientific and clinical implications. Implementation of ACT-Discover could lead to important insights from existing cohorts, or underpin future prospective studies seeking to characterise the landscape of tumour evolution through liquid biopsy.

  Study EGAS00001007077

• Effects of 2DG, Galactose, or Oligomycin on the epigenome remodeling induce by T cells activation.

  T cell activation is known to require a metabolic shift towards glycolysis, triggered by TCR-engagement . To substantiate the role of glycolysis in transcriptional reprogramming during antigen-driven T cell activation, we utilized a well-established ex-vivo model to recapitulate the process of antigen stimulation without the requirement for additional antigen-presenting cells . We initially explored and validated the conditions for both metabolic and epigenetic reprogramming. Antigen-mediated CD4+ T cells activation initiates context-specific gene-expression programs that drive effector functions and cell fates through changes in the epigenetic landscape . To evaluate global changes in histone modifications, we isolated CD4+ T cells from peripheral blood (PB) of healthy control (HC) human donors and either maintained them in resting conditions or stimulated with anti-CD3/CD28 for 24h. Rapid changes in the CD4+ epigenome and transcriptome were observed within 24h as measured by RNA-sequencing (RNA-seq) and H3K27ac chromatin immunoprecipitation followed by DNA-sequencing (ChIP-seq).To explore the requirement for metabolic reprogramming in regulating transcriptional responses after initial TCR engagement, CD4+ T cells from PB of healthy donors were isolated and activated in either the presence or absence of oligomycin (inhibitor of ATP synthase blocking the production of ATP by OXPHOS), 2-deoxyglucose (2DG; competitively inhibiting production of glucose-6-phosphate from glucose), or replacement of glucose with galactose. Galactose is metabolized through glycolysis, utilizing the Leloir pathway which is energy neutral, forcing cells to use all the pyruvate production to sustain the TCA cycle for OXPHOS and ATP synthesis at the expense of lactate production. Employing ChIP-seq and RNA-seq, inhibition of OXPHOS-mediated ATP production by oligomycin had only a limited effect on H3K27-acetylation or activation-induced transcription during the first 24h of activation. In contrast, inhibition of glycolysis by 2DG treatment compromised both H3K27-acetylation and transcriptional changes. However, under conditions of galactose-metabolism where lactate production is abrogated, activated CD4+ T cells were still able to globally remodel H3K27 acetylation and reprogram the transcriptome.

  Study EGAS00001007115

• Multi-omic and functional analysis for classification and treatment of sarcomas with FUS-TFCP2 or EWSR1-TFCP2 fusions(H021/INF)

  Many rare cancers are not well understood, and pathogenesis-directed therapies are often lacking, resulting in poor patient outcomes. Leveraging two (inter)national precision oncology trials that enroll large numbers of rare cancers, we investigated the clinical, histopathologic, molecular, and functional characteristics of rhabdomyosarcoma (RMS) with fusions of FUS or EWSR1 to the TFCP2 transcription factor, a recently discovered ultra-rare entity whose classification, pathogenesis, and optimal treatment are unclear. Unusually for fusion-driven sarcomas, most cases had highly rearranged genomes, including chromothripsis, and signs of defective homologous recombination DNA repair. All tumors showed extreme expression of a truncated TERT variant and the ALK receptor tyrosine kinase, which was additionally affected by intragenic deletions and aberrant splicing, resulting in the expression of shortened variants in 58% of cases. Three ALK variants were oncogenic in immortalized cells, and patient-derived tumor cells expressing two variants responded to certain ALK inhibitors. Other recurrent alterations included CDKN2A/MTAP co-deletions and mutations in PAPPA2, encoding an IGFBP5-specific proteinase, in 67% and 25% of cases, respectively. DNA methylation analysis, along with 19 other soft-tissue sarcoma classes, revealed a close relationship with undifferentiated sarcoma but not other RMS subtypes, suggesting that TFCP2-rearranged RMS is a separate entity, possibly arising from a distinct cell of origin. Expression of TFCP2 fusions in immortalized human cells blocked late myogenic differentiation and significantly induced genes that were also highly expressed in patient tumors, including ALK, TERT, and two known regulators of skeletal muscle cells, IGFBP5 and PTH1R. Genome-wide chromatin profiling demonstrated direct binding of FUS-TFCP2 to the ALK and TERT loci outside their regular promoters, which correlated with the expression of alternative transcript variants. Finally, FUS-TFCP2 inhibited the repair of DNA double-strand breaks in immortalized cells, rendering them sensitive to treatment with cisplatin. This study provides insight into the pathogenesis and therapeutic vulnerabilities of a new RMS subtype and illustrates the value of linking comprehensive molecular diagnostics and functional annotation within multicenter consortia to improve the understanding and clinical management of rare cancers.

  Study EGAS00001006939

• UK10K NEURO IOP COLLIER

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The Genetics and Psychosis (GAP) set consists of samples from subjects with schizophrenia, ascertained as a new-onset sample. This set is of UK origin, with data on cognition, brain imaging and other endophenotypes. The Maudsley twin series consists of probands ascertained from the Maudsley Twin Register, defined as patients of multiple birth who had suffered psychotic symptoms. This set is of UK origin, with data on cognition, brain imaging and other endophenotypes, with DNA available from an MZ or DZ affected or unaffected co-twin. The Maudsley family study (MFS) consists of over 250 families who have a history of schizophrenia or bipolar disorder Within the Maudsley Family Study, biological markers of psychosis include neuropsychological tests, Evoked Response Potentials Tests (ERPs), MRI scans, dermatoglyphics and eye tracking. Early risk factors for psychosis and clinical symptoms are also investigated. This set is of UK origin, with DNA available from both affected and unaffected relatives in many of the probands. For further information on this study please contact David Collier (david.collier@kcl.ac.uk).

  Study EGAS00001000121

• UK10K NEURO GURLING

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. This sample set consists of DNA from multiply affected schizophrenia families. The families have been diagnosed using the SADS-L clinical instrument which gives diagnoses at the probable level of the research diagnostic criteria (RDC). In addition all diagnoses are available using DSMIIIR criteria. These criteria are widely accepted as being valid and reliable for the diagnosis of schizophrenia. All families have been collected to ensure that they are uni-lineal for transmission of schizophrenia, i.e. they have only one affected parent with schizophrenia, or a relative of only one transmitting or obligate carrier parent with schizophrenia. Families with bi-lineal transmission of schizophrenia (i.e. with both parents being affected) were not sampled for this study. All families have multiple cases of schizophrenia and related disorders. The families have been selected to ensure there are no cases of bipolar disorder within them and that they do not contain bipolar disorder in any relatives on either side of the family.For further information on this cohort please contact Hugh Gurling (h.gurling@ucl.ac.uk).

  Study EGAS00001000225

• Therapeutic Resistance to PI3K-alpha Inhibitors

  We have studied the tumor genomic evolution in a patient with metastatic breast cancer bearing an activating PIK3CA mutation. The patient was treated with the PI3Kα inhibitor BYL719 and achieved a lasting clinical response on BYL719, but eventually progressed to treatment and died shortly thereafter. A rapid autopsy was performed and a total of 14 metastatic lesions were collected for further analysis. In order to identify possible genetic determinants of acquired resistance to PI3Kα inhibition, we took a three-step approach. First, we examined both the primary tumor (before BYL719 treatment) and the new lung metastasis by whole genome sequencing (DNA from the spleen was used as a normal control). Then, we analyzed the primary tumor, lung metastasis, and the peri-aortic lesion that remained stable (responding) at the time of progression to BYL719 therapy by whole exome sequencing. Finally, to confirm and expand our findings, we sequenced the primary tumor and all the metastatic lesions to >500-fold coverage using a custom targeted deep-sequencing assay, MSK-IMPACT. Using this 3-step approach, we have demonstrated that patient’s lesion underwent parallel genetic evolution at multiple metastatic sites with different PTEN genomic alterations leading to a convergent PTEN- null phenotype resistant to PI3Kα inhibition. In order to expand our observations, we analyzed paired samples (pre-treatment and at progression) from six additional patients enrolled in the BYL719 trial. Acquired bi-allelic loss of PTEN was found in one additional patient treated with BYL719 whereas in two patients PIK3CA mutations present in the primary tumor were no longer detected at the time of progression. To functionally characterize our findings, inducible PTEN knockdown in sensitive cells resulted in resistance to BYL719, while simultaneous PI3K-p110β blockade reverted this resistance phenotype, both in cell lines and in PTEN-null xenografts derived from our patient. We conclude that parallel genetic evolution of separate sites with different PTEN genomic alterations leads to a convergent PTEN-null phenotype resistant to PI3Kα inhibition.

  Study EGAS00001000991

• Genomic analysis of patient-derived xenograft models reveals intratumor-heterogeneity in endometrial cancer and can predict tumor growth inhibition with talazoparib

  Background: Endometrial cancer (EC) is a major gynecological cancer with increasing incidence. It comprised of four molecular subtypes with differing etiology, prognoses, and response to chemotherapy. In the future, clinical trials testing new single agents or combination therapies will be targeted to the molecular subtype most likely to respond. Pre-clinical models that faithfully represent the molecular subtypes of EC are urgently needed, we sought to develop and characterize a panel of novel EC patient-derived xenograft (PDX) models. Methods: Here, we report whole exome or whole genome sequencing of 11 PDX models and the matched primary tumor. Analysis of multiple PDX lineages and passages was performed to study tumor heterogeneity across lineages and/or passages. Based on recent reports of frequent defects in the homologous recombination (HR) pathway in EC, we assessed mutational signatures and HR deficiency scores and correlated these with in vivo responses to the PARP inhibitor (PARPi) talazoparib in six PDXs representing the copy number high/p53-mutant and mismatch-repair deficient molecular subtypes of EC. Results: PDX models were successfully generated with similar histomorphology from all four molecular subtypes of EC and uterine carcinosarcomas. We also observed a wide range of inter-tumor and intra-tumor heterogeneity across different PDX lineages, however one arm of most PDX models recapitulated the molecular landscape of the primary tumor without major genomic drift. An in vivo response to talazoparib was detected in four copy number high models. Two models (carcinosarcomas) showed a response consistent with stable disease and two models (one copy number high serous EC and another carcinosarcoma) showed significant tumor growth inhibition, albeit one consistent with progressive disease, however all lacked the HR deficiency genomic signature. Conclusions: EC PDX models represent the four molecular subtypes of disease and can capture intra-tumoral heterogeneity of the original primary tumor. PDXs of the copy number high molecular subtype showed sensitivity to PARPi, however, deeper and more durable responses will likely require combination of PARPi with other agents.

  Study EGAS00001004666

• Histological Transformation and Progression in Follicular Lymphoma: a Clonal Evolution Study

  Background: Follicular lymphoma (FL) is an indolent, yet incurable B-cell malignancy. A subset of patients experience increased mortality rate driven by two distinct clinical end points: histological transformation and early progression after immuno-chemotherapy. The nature of tumor clonal dynamics leading to these clinical endpoints is poorly understood and previously determined genetic alterations do not explain the majority of transformed cases or accurately predict early progressive disease. We contend that detailed knowledge of the expansion patterns of specific cell populations plus their associated mutations would provide insight into therapeutic strategies and disease biology over the timecourse of FL clinical histories. Methods and Findings: Using a combination of whole genome sequencing, targeted deep sequencing and digital droplet PCR on matched diagnostic and relapse specimens, we deciphered the constituent clonal populations in n=15 transformation patients and n=6 progressors, and measured their change in abundance over time. We observed widely divergent patterns of clonal dynamics in transformed cases relative to progressed cases. Transformation specimens were generally composed of clones that were rare or absent in diagnostic specimens, consistent with dramatic clonal expansions that came to dominate the transformation specimen. This pattern was independent of time to transformation and treatment modality. By contrast, progression specimens were composed of clones that were already present in the diagnostic specimens and exhibited only moderate clonal dynamics, even in the presence of immuno-chemotherapy. Analysis of somatic mutations impacting 94 genes was undertaken in a validation cohort consisting of 395 samples from 277 patients in order to decipher disrupted biology in the two clinical endpoints. We found 12 genes were more commonly mutated in transformed samples than in preceding FL, including TP53, B2M, CCND3, GNA13, S1PR2 and P2RY8. Moreover, 10 genes were more commonly mutated in diagnostic specimens of early progressors including TP53, BTG1, MKI67 and XBP1. Conclusions: Our results illuminate contrasting modes of evolution shaping the clinical histories of transformation and progression. They have implication for interpretation of evolutionary dynamics in the context of treatment-induced selective pressures, and indicate that transformation and progression will require different clinical management strategies.

  Study EGAS00001001709

• Foundation Medicine Genomic Data Used to Identify Prognostic Markers and Fusion Genes in Multiple Myeloma

  The purpose of this study is to identify prognostic markers and treatment targets using a clinically certified sequencing panel in multiple myeloma. Mutational burden was associated with maf and proliferation gene expression groups, and a high-mutational burden was associated with a poor prognosis. We identified homozygous deletions that were present in multiple myeloma within key genes, including CDKN2C, RB1, TRAF3, BIRC3 and TP53, and that bi-allelic inactivation was significantly enriched at relapse. Alterations in CDKN2C, TP53, RB1 and the t(4;14) were associated with poor prognosis. Alterations in RB1 were predominantly homozygous deletions and were associated with relapse and a poor prognosis which was independent of other genetic markers, including t(4;14), after multivariate analysis. Bi-allelic inactivation of key tumor suppressor genes in myeloma was enriched at relapse, especially in RB1, CDKN2C and TP53 where they have prognostic significance. In addition, chromosomal rearrangements that result in oncogenic kinase activation are present in many solid and hematological malignancies, but none have been reported in multiple myeloma (MM). Here we detected fusion genes in 1.5% of patients. These fusion genes were in-frame and the majority of them contained kinase domains from either receptor tyrosine kinases (ALK, ROS1, NTRK3, and FGFR1) or cytoplasmic kinases (BRAF, MAP3K14, and MAPK14) which would result in the activation of MEK/ERK, NF-κB or inflammatory signaling pathways. Fusion genes were present in smoldering MM, newly diagnosed MM and relapse patient samples indicating they are not solely late events. Most fusion genes were subclonal in nature, but one EML4-ALK fusion was clonal indicating it is a driver of disease pathogenesis. Samples with fusions of receptor tyrosine kinases were not found in conjunction with clonal Ras/Raf mutations indicating a parallel mechanism of MEK/ERK pathway activation. Fusion genes involving MAP3K14 (NIK), which regulates the NF-κB pathway, were detected as were t(14;17) rearrangements involving NIK in 2% of MM samples. Activation of kinases in myeloma through rearrangements presents an opportunity to use treatments existing in other cancers.

  Study EGAS00001002874

• Scottish High Grade Serous Ovarian Cancer

  Scottish HGSOC samples were collected via local Bioresource facilities at Edinburgh, Glasgow, Dundee and Aberdeen and stored in liquid Nitrogen until required. HGSOC patients were determined from pathology records and were included in the study where there was matched tumour and whole blood samples. On receipt of tumour material the tumour was processed as follows: firstly, the tumour sample was divided into two for DNA and RNA extraction. Slivers of tissue were cut from the front and rear faces of the DNA sample, then fixed in formalin and embedded in paraffin wax. Sections from the front and rear tissues from all samples were examined by H&E staining supplemented by WT1/p53 immunohistochemistry if required. Following pathology review, samples were only included if they met the following criteria: they were confirmed as HGSOC and there was greater than 40% tumour cellularity throughout the tumour, determined using the H&E sections. Somatic DNA was extracted using the Qiagen DNeasy Blood and tissue kit (cat no 69504). The tissue was initially homogenised using a Qiagen Bioruptor, followed by the manufacturers recommended protocol (including RNase digestion step). Germline DNA was extracted from 1-3ml whole blood using the Qiagen FlexiGene kit (cat no 51206) following the manufacturers recommended protocol. The resulting DNA underwent quality control as follows: firstly, A260 and A280nm were measured on a Denovix DS-11 Fx to qualitatively illustrate A260/280nm and A260/230nm ratios as surrogate measures of DNA purity. A260/280 had to be 1.8 or greater and A260/230 had to be 2.0 or greater. Then, DNA was quantified using LifeTechnologies Qubit dsDNA BR kit (cat no Q32850) and we required a minimum of 50ul at 25ng/ul for WGS. Thirdly, DNA was diluted to 25ng/ul and a representative sample was loaded onto a 0.8% TAE gel, ran at 100v for 60mins and then imaged using a BioRad ChemiDoc imaging system to visualise the DNA quality. Only when all 4 quality control requirements were satisfied was the DNA sequenced. The DNA was sequenced at the Glasgow Precision Oncology Laboratories.

  Study EGAS00001004410

• Germline elongator mutations in sonic hedgehog medulloblastoma

  Cancer genomics has illuminated a wide spectrum of genes and core molecular processes contributing to human malignancy. Still, the genetic and molecular basis of many cancers remains only partially explained. Genetic predisposition accounts for 5-10% of cancer diagnoses and genetic events cooperating with known somatic driver events are poorly understood. Analyzing established cancer predisposition genes in medulloblastoma (MB), a malignant childhood brain tumor, we recently identified pathogenic germline variants that account for 5% of all MB patients. Here, by extending our previous analysis to include all protein-coding genes, we discovered and replicated rare germline loss-of-function (LoF) variants across Elongator Complex Protein 1 (ELP1) on 9q31.3 in 15% of pediatric MBSHH cases, thus implicating ELP1 as the most common MB predisposition gene and increasing genetic predisposition to 40% for pediatric MBSHH. Inheritance was verified based on parent-offspring and pedigree analysis, which identified two families with a history of pediatric MB. ELP1-associated MBs were restricted to the molecular SHHa subtype and were characterized by universal biallelic inactivation of ELP1 due to somatic loss of chromosome 9q. The majority of ELP1-associated MBs exhibited co-occurring somatic PTCH1 (9q22.32) alterations, suggesting that ELP1-deficiency predisposes to tumor development in combination with constitutive activation of SHH signaling. ELP1 is an essential subunit of the evolutionary conserved Elongator complex, whose primary function is to enable efficient translational elongation through tRNAs modifications at the wobble (U34) position. Biochemical, transcriptional, and proteomic analyses revealed that ELP1-associated MBSHH are characterized by a destabilized core Elongator complex, loss of Elongator-dependent tRNA modifications, codon-dependent translational reprogramming, and induction of the unfolded protein response (UPR), consistent with deregulation of protein homeostasis due to Elongator-deficiency in model systems. Our findings suggest that genetic predisposition to proteome instability is a previously underappreciated determinant in the pathogenesis of pediatric brain cancer. These results provide strong rationale for further investigating the role of protein homeostasis in other pediatric and adult cancer types and potential opportunities for novel therapeutic interference.

  Study EGAS00001004126

• Gene expression-based prediction of pazopanib efficacy in sarcoma (HIPO, H021)

  The multi-receptor tyrosine kinase (RTK) inhibitor pazopanib is approved for the treatment of advanced non-adipocytic soft-tissue sarcoma and has also shown activity in other sarcoma subtypes. However, its clinical efficacy is highly variable, and no reliable predictors exist to date to select patients who are most likely to benefit from this drug. We analyzed the molecular profiles and clinical outcomes of pazopanib-treated sarcoma patients enrolled in a prospective observational study by the German Cancer Consortium, DKTK MASTER, that employs whole-genome/exome sequencing (WGS/WES) and transcriptome sequencing to inform the care of adults with advanced cancer across histologies who are younger than 51 and patients with rare tumors, including rare subtypes of more common entities, regardless of age.​ Among a total of 109 patients with available WGS/WES data, there was no correlation between clinical parameters, specific genetic alterations, or mutational signatures and clinical outcome. In contrast, analysis of a subcohort of 62 patients who underwent molecular analysis before pazopanib treatment and had transcriptome sequencing data available showed that mRNA levels of NTRK3 (hazard ratio [HR]=0.53, p=0.021), IGF1R (HR=1.82, p=0.027), and KDR (HR=0.50, p=0.011) were independently associated with progression-free survival (PFS). Based on the expression of these RTK genes, i.e., the features NTRK3-high, IGF1R-low, and KDR-high, we developed a pazopanib efficacy predictor (PEP) that stratified patients into three groups with significantly different PFS (p<0.0001). Application of the PEP to an independent cohort of pazopanib-treated sarcoma patients from DKTK MASTER (n=43) confirmed its potential to separate patient groups with significantly different PFS (p=0.02), whereas no such association was observed in sarcoma patients from DKTK MASTER (n=77) or The Cancer Genome Atlas sarcoma cohort (n=256) who were not treated with pazopanib. A score based on the combined expression of NTRK3, IGF1R, and KDR allows identification of sarcoma patients with good, intermediate, and poor outcome following pazopanib therapy and warrants investigation as a predictive tool in prospective studies to optimize the use of this drug in the clinic.

  Study EGAS00001005836

• Evolution and clinical impact of genetic epistasis within EGFR-mutant lung cancers

  The current understanding of tumorigenesis is largely centered on a monogenic driver oncogene model. This paradigm is incompatible with the prevailing clinical experience in most solid malignancies: monotherapy with a drug directed against an individual oncogenic driver typically results in incomplete clinical responses and eventual tumor progression1-7. By profiling the somatic genetic alterations present in over 2,000 cases of lung cancer, the leading cause of cancer mortality worldwide, we show that combinations of functional genetic alterations, i.e. genetic collectives dominate the landscape of advanced-stage disease. We highlight this polygenic landscape and evolution of advanced-stage non-small cell lung cancer (NSCLC) through the spatial-temporal genomic profiling of 7 distinct tumor biopsy specimens and 6 plasma specimens obtained from an EGFR-mutant NSCLC patient at (1) initial diagnosis of early-stage disease, (2) metastatic progression, (3) sequential treatment and resistance to 2 EGFR inhibitors, (4) death. The comprehensive genomic analysis of this case, coupled with circulating free (cf) tumor DNA profiling of additional advanced-stage EGFR-mutant NSCLC clinical cohorts with associated treatment responses uncovered features of evolutionary selection for multiple concurrent gene alterations: including the presence of EGFR inhibitor-sensitive (EGFRL858R;EGFRexon19del) or inhibitor-resistant (EGFRT790M;EGFRC797S) forms of oncogenic EGFR along with cell cycle gene alterations (e.g. in CDK4/6, CCNE1, RB1) and activating alterations in WNT/β-catenin and PI3K pathway genes, which our data suggest can cooperatively impart non-redundant functions to limit EGFR targeted therapy response and/or promote tumor progression. Moreover, evidence of an unanticipated parallel evolution of both EGFR T790M and two distinct forms of oncogenic PIK3CA was observed. Our study provides a large-scale clinical and genetic dataset of advanced-stage EGFR-mutant NSCLC, a rationale for specific polytherapy strategies such as EGFR and CDK4/6 inhibitor co-treatment to potentially enhance clinical outcomes, and prompts a re-evaluation of the prevailing paradigm of monogenic-based molecular stratification for targeted therapy. Instead, our findings highlight an alternative model of genetic collectives that operate through epistasis to drive lung cancer progression and therapy resistance.

  Study EGAS00001002604

• Identifying transcriptional programs underlying anti-EGFR small molecule response and resistance with TraCe-seq

  Genetic and non-genetic heterogeneity within cancer cell populations represents a major challenge to anti-cancer therapies. We currently lack robust methods to determine how pre-existing and adaptive features affect cellular responses to therapies. Here, by conducting clonal fitness mapping and transcriptional characterization using expressed barcodes and single-cell RNA-sequencing, we have developed TraCe-seq, a method that captures at clonal resolution the origin, fate, and differential early adaptive transcriptional programs of cells in a complex population in response to distinct treatments. We used TraCe-seq to benchmark how next-generation dual EGFR inhibitors-degraders compare to standard EGFR kinase inhibitors in EGFR-mutant lung cancer cells. To QC the TraCe-seq strategy, single-cell RNA-seq libraries were generated from a variety of human cancer cell lines transduced with the TraCe-seq library to validate the TraCe-seq strategy. Specifically, 5 different cell lines (PC9, MCF-10A, MDA-MB-231, NCI-H358, and NCI-H1373) were each transduced with a unique TraCe-seq barcode. The transduced cells were selected with puromycin only, dissociated to single cell suspensions, and then mixed together. The complex mixture of the 5 cell lines was profiled by 10X scRNA-seq. Furthermore, transduced NCI-H1373 cells were sorted by FACS to enrich for the top 50% of eGFP positive cells, and sorted cells were cultured briefly and used to construct scRNA-seq libraries and profiled by 10x scRNA-seq. To carry out the full TraCe-seq experiment, ~600 PC9 cells carrying unique TraCe-seq barcodes were expanded over 12 doublings to establish the barcoded population. A subset of the barcoded PC9 population was used to generate scRNA-seq libraries and profiled by 10x scRNA-seq prior to treatment to establish a baseline transcription profile for each barcoded clone. The rest of the cells were then treated for four days with 1 µM erlotinib, 1 µM GNE-069, or 1 µM GNE-104 respectively. scRNA-seq libraries were then generated form the treated cells and profiled by 10x scRNA-seq.

  Study EGAS00001005405

• Characterization of Cell Free Plasma Methyl-DNA From Xenografted Tumors to Guide the Selection of Diagnostic Markers for Early-Stage Cancers

  Circulating cell-free methyl-DNA (mcfDNA) contains promising cancer markers but its low abundance and possibly diverse origin pose challenges toward the accurate diagnosis of early stage cancers. By whole-genome bisulfite sequencing (WGBS) of cell-free DNA (cfDNA) from about 0.5 mL plasma of mice xenografted with human tumors, we obtained and aligned the reads to the human genome, filtered out the mouse and carrier bacterial sequences, and confirmed the tumor origin of methyl-cfDNA (mctDNA) by methylation-sensitive restriction enzyme digestion prior to species-specific PCR. We estimated that human tumor-specific reads (ctDNA) or mctDNA comprised about 0.29 or 0.01%, respectively of the xenograft mouse cfDNA, and about 0.029 or 0.001% of the cfDNA of human early stage cancer patients. Similar WGBS of early stage (0-II, node- and metastasis-free) breast, lung or colorectal cancer samples identified hundreds of specific DMRs (differentially methylated regions) compared to healthy controls. Their association with tumourigenesis was supported by stage-dependent methylation, tumor suppressor or oncogene clusters, and genes also identified in the xenograft samples. Using 20 three-cancer-common and 17 colorectal cancer-specific DMRs in combination (top 0.0018% of the WGBS methylation clusters) was sufficient to distinguish the stage I colorectal cancers from breast and lung cancers and healthy controls. Our data thus confirmed the tumor origin of mctDNA by sequence specificity, and provide a selection threshold for authentic tumor mctDNA markers toward precise diagnosis of early stage cancers solely by top DMRs in combination. Biological Materials were provided by the Ontario Tumour Bank, which is funded by the Ontario Institute for Cancer Research Reference for citation: Ling Liu, Jinghua Feng, Julian Polimeni, Manli Zhang, Hai Nguyen, Urmi Das, Xu Zhang, Harminder Singh, Xiao-Jian Yao, Etienne Leygue, Sam K.P. Kung, and Xie J. Characterization of cell free plasma methyl-DNA from xenografted tumours to guide the selection of diagnostic markers for early-stage cancers. Frontiers in Oncology, 2021 Feb. 5. DOI: 10.3389/fonc.2021.615821.

  Study EGAS00001006175

• Chromatin profiles classify castration-resistant prostate cancers suggesting therapeutic targets

  Untreated prostate cancers rely on androgen receptor (AR) signaling for growth and survival, forming the basis for the initial efficacy of androgen deprivation therapy (ADT). Yet, the disease can relapse and progress to a lethal stage termed castration-resistant prostate cancer (CRPC). Reactivation of AR signaling represents the most common driver of CRPC growth and next generation AR signaling inhibitors (ARSIs) are now used in combination with ADT as a first line therapy. However, ARSIs can result in selective pressure generating AR-independent tumors. The transition from AR-dependence frequently accompanies a change in phenotype resembling developmental trans-differentiation or ‘lineage plasticity’. Neuroendocrine prostate cancer, which lacks a defined pathologic classification, is the most studied type of lineage plasticity. However, most AR-null tumors do not exhibit neuroendocrine features and are classified as ‘double-negative prostate cancer’, the drivers of which are poorly defined. Lineage plasticity studies in CRPC are limited by the lack of genetically defined patient-derived models that recapitulate the disease spectrum. To address this, we developed a biobank of organoids generated from patient biopsies to study the landscape of metastatic CRPC and allow for functional validation assays. Proteins called transcription factors (TFs) are drivers of tumor lineage plasticity. To identify the key TFs that drive the growth of AR-independent tumors, we integrated epigenetic and transcriptomic data generated from CRPC models. By presenting a map of the chromatin accessibility and transcriptomic landscape of CRPC using a diverse biobank of organoids, PDXs and cell lines that recapitulate the heterogeneity of metastatic prostate cancer, we validate CRPC-AR and CRPC-NE subtypes and reveal two subtypes of AR-negative/low samples as well as their respective masterTFs. Additional analysis revealed a model in which YAP, TAZ, TEAD and AP-1 function together and drive oncogenic growth in CRPC-SCL samples. Overall, we show here how an approach to stratify CRPC patients into four subtypes using their transcriptomic signatures can potentially inform upon appropriate clinical decisions.

  Study EGAS00001006059

• HCA_Thymus_Disease_DiGeorge_Transplant_RNA

  Single cell transcriptomics study of thymic transplant biopsies Allogeneic thymus transplantation is the only curative therapy for complete DiGeorge Syndrome (cDGS), a rare severe primary immunodeficiency characterised by athymia. GOSH is one of only two centres worldwide to offer this treatment. Despite a lack of major histocompatibility complex (MHC)-matching between donor and host, transplanted thymus becomes repopulated by recipient bone marrow derived precursor cells and supports development of functional T-cells. The mechanisms underlying thymopoiesis in this context are poorly understood, but over time we observe reconstitution of T-cell immunity, with the ability to produce host naïve T-cells showing a broad T-cell receptor (TCR) repertoire and to generate MHC-restricted T-cell proliferative responses. Although lifesaving, the achieved immunological reconstitution is typically not complete with circulating T-cell numbers usually remaining below the age related normal ranges. Additionally, we observe persistence of donor-derived T-cells of unknown clinical significance. To gain more insight into the mechanisms by which MHC-mismatched transplanted thymus supports T-cell development with self-tolerance, as well as into the basis of suboptimal T-cell immunity, we now aim to investigate immune reconstitution after thymus transplantation in further detail by using single-cell transcriptomics, applied to thymic transplant biopsies and peripheral blood samples collected during standard post-transplant patient care. By identifying which lineages of host- derived cells repopulate the thymic tissue after transplantation, we will be able to address the role of MHC in positive and negative T-cell selection during T-cell differentiation. We will also be able to clarify the exact ontogeny of the persistent donor T-cells, as well as their possible role. Understanding the mechanisms of action of HLA-mismatched transplanted thymus will contribute to treatment optimisation. Additionally, our research provides a unique opportunity to further investigate key immunological concepts, such as tolerance and autoimmunity, challenging existing paradigms in thymus immunology. 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 EGAS00001004025

• High-resolution testing of ctDNA dynamics predicts survival in metastatic NSCLC

  One of the great challenges of therapeutic oncology is determining who might achieve survival benefit from a particular therapy. Circulating tumor DNA (ctDNA) provides real-time assessments of patient prognosis and response to treatment using a simple blood draw. While ctDNA positivity is established as a poor prognostic factor, studies on longitudinal ctDNA dynamics have been small and non-randomized, with ctDNA assessments done at disparate time points. To address this, we performed high-sensitivity longitudinal ctDNA testing in 466 patients across 5 time points (1,954 samples total) in a randomized phase III study comparing different chemotherapy-immunotherapy combinations. We leverage machine learning to jointly model multiple ctDNA metrics to predict overall survival in a training/testing framework. _ . Treatment initiation correlated with reductions in ctDNA levels, and training of our machine learning model suggests that assessment of ctDNA dynamics at C3D1 (cycle 3 day 1) of chemo-IO treatment may be optimal to predict OS. The model performs well in the hold-back test data, enabling stratification of patients with Stable Disease (SD) into high-risk vs low-intermediate-risk (HR = 3.2 [2.0-5.3], p <0.001; median 7.1 versus 22.3 months respectively); similarly, the model stratifies patients with a Partial Response (PR) (HR =3.3 [1.7-6.4], p <0.001; median 8.8 versus 28.6 months). Importantly, the model validates well in an external cohort of patients in a different treatment setting and assayed with a different ctDNA technology, in which model predictions similarly identified high-risk patients (OS HR=3.73 [1.83-7.60], logrank p=0.00012). Simulations of clinical trial scenarios employing our ctDNA model further suggest that early ctDNA testing outperforms early radiographic imaging for predicting trial outcomes (increasing the rate of ‘True Go’ decisions by 5.2 - 22.8% depending on the drug combination). Overall, we show that measuring ctDNA dynamics during the course of therapy dramatically improves patient risk stratification, and may provide a means to differentiate between competing therapies at an early time point during clinical trials.

  Study EGAS00001006703

• Influence of Genomic Landscape on Cancer Immunotherapy for Newly Diagnosed Ovarian Cancer: Biomarker Analyses from the IMagyn050 Randomized Clinical Trial

  The IMagyn050 trial (NCT03038100), which evaluated the efficacy of Atezo vs placebo (Pla) with carboplatin, paclitaxel and bevacizumab (CPB) in front line ovarian cancer patient (pts), did not meet its co-primary endpoints of PFS in ITT or PD-L1+ (Moore et al. JCO 2021). The aim of this exploratory study was to evaluate whether patients with BRCA1/2-mutated or homologous recombination-deficient (HRD) ovarian cancers benefitted from atezolizumab in the phase III IMagyn050 (NCT03038100) trial. Methods: Patients with newly diagnosed ovarian cancer were randomized to either atezolizumab or placebo with standard chemotherapy and bevacizumab. PD-L1 status of tumor-infiltrating immune cells was determined centrally (VENTANA SP142 assay). Genomic alterations, including deleterious BRCA1/2 alterations, genomic loss of heterozygosity (gLOH), tumor mutation burden (TMB), and microsatellite instability (MSI), were evaluated using the FoundationOne assay. HRD was defined as gLOH ≥16%, regardless of BRCA1/2 mutation status. Potential associations between progression-free survival (PFS) and genomic biomarkers were evaluated using standard correlation analyses and log-rank of Kaplan-Meier estimates. Results Among biomarker-evaluable samples, 22% (234/1050) harbored BRCA1/2 mutations and 46% (446/980) were HRD. Median TMB was low irrespective of BRCA1/2 or HRD. Only 3% (29/1024) had TMB ≥10 mut/Mb and 0.3% (3/1022) were MSI-high. PFS was better in BRCA2-mutated versus BRCA2-non-mutated tumors and in HRD versus proficient tumors. PD-L1 positivity (≥1% expression on immune cells) was associated with HRD but not BRCA1/2 mutations. PFS was not improved by adding atezolizumab in BRCA2-mutated or HRD tumors; there was a trend toward enhanced PFS with atezolizumab in BRCA1-mutated tumors. Conclusion Most ovarian tumors have low TMB despite BRCA1/2 mutations or HRD. Neither BRCA1/2 mutation nor HRD predicted enhanced benefit from atezolizumab. This is the first randomized double-blind trial in ovarian cancer demonstrating that genomic instability triggered by BRCA1/2 mutation or HRD is not associated with improved sensitivity to immune checkpoint inhibitors.

  Study EGAS00001006838

• TK-EPN862 - Patient-dervied xenograft model of Posterior Fossa A Ependymoma - RNAseq

  Ependymomas (EPN) are the third most common malignant brain cancer in children. Treatment strategies for pediatric EPN have remained unchanged over recent decades, with 10-year survival rates stagnating at just 67% for children aged 0-14 years. Moreover, patients who survive treatment for the most common subtype of pediatric EPN (posterior fossa A (PFA)), often suffer long-term neurological side effects as a result of therapy. It is evident that there is a need for safer, more effective treatments for pediatric EPN patients. There are ten distinct subgroups of EPN, each with their own molecular and prognostic features. To identify and facilitate the testing of new treatments for EPN, in vivo laboratory models representative of the diverse molecular subtypes are required. Here, we describe the establishment of a patient-derived orthotopic xenograft (PDOX) model of PFA EPN, derived from a metastatic cranial lesion. Patient and PDOX tumors were analyzed using immunohistochemistry, DNA methylation profiling, whole genome sequencing (WGS) and RNA sequencing. Both patient and PDOX tumors classified as PFA EPN by methylation profiling, and shared similar histological features consistent with this molecular subgroup. RNA sequencing revealed that gene expression patterns were maintained across the primary and metastatic tumors, as well as the PDOX. Copy number profiling revealed gains of chromosomes 7, 8 and 19, and loss of chromosomes 2q and 6q in the PDOX and matched patient tumor. No clinically significant single nucleotide variants were identified, consistent with the low mutation rates observed in PFA EPN. Overexpression of EZHIP RNA and protein, a common feature of PFA EPN, was also observed. Despite the aggressive nature of the tumor in the patient, this PDOX was unable to be maintained past two passages in vivo. Others who have successfully developed PDOX models report some of the lowest success rates for EPN compared to other pediatric brain cancer types attempted, with loss of tumorigenicity not uncommon, highlighting the challenge in propagating these tumors in the laboratory. Here, we discuss our collective experiences with PFA EPN PDOX model generation and propose potential approaches to improve future success in establishing preclinical EPN models

  Study EGAS00001006844

• Mapping the epigenomic landscape of human monocytes following innate immune activation reveals context-specific mechanisms driving endotoxin tolerance

  Background Monocytes are key mediators of innate immunity to infection, undergoing profound and dynamic changes in epigenetic state and immune function which are broadly protective but may be dysregulated in disease. Here, we aimed to advance understanding of epigenetic regulation following innate immune activation, acutely and in endotoxin tolerant states. Methods We exposed human primary monocytes from healthy donors (n=6) to interferon-γ or differing combinations of endotoxin (lipopolysaccharide), including acute response (2hr) and two models of endotoxin tolerance: repeated stimulations (6+6hr) and prolonged exposure to endotoxin (24hr). Another subset of monocytes was left untreated (naïve). We identified context-specific regulatory elements based on epigenetic signatures for chromatin accessibility (ATAC-seq) and regulatory non-coding RNAs from total RNA sequencing. Results We present an atlas of differential gene expression for endotoxin and interferon response, identifying widespread context specific changes. Across assayed states, only 24-29% of genes showing differential exon usage are also differential at the gene level. Overall, 19.9% (6884 of 34616) of repeatedly observed ATAC peaks were differential in at least one condition, the majority upregulated on stimulation and located in distal regions (64.1% vs 45.9% of non-differential peaks) within which sequences were less conserved than non-differential peaks. We identified enhancer-derived RNA signatures specific to different monocyte states that correlated with chromatin accessibility changes. The endotoxin tolerance models showed distinct chromatin accessibility and transcriptomic signatures, with integrated analysis identifying genes and pathways involved in the inflammatory response, detoxification, metabolism and wound healing. We leveraged eQTL mapping for the same monocyte activation states to link potential enhancers with specific genes, identifying 1,946 unique differential ATAC peaks with 1,340 expression associated genes. We further use this to inform understanding of reported GWAS, for example involving FCHO1 and coronary artery disease. Conclusion This study reports context-specific regulatory elements based on transcriptomic profiling and epigenetic signatures for enhancer-derived RNAs and chromatin accessibility in immune tolerant monocyte states, and demonstrates the informativeness of linking such elements and eQTL to inform future mechanistic studies aimed at defining therapeutic targets of immunosuppression and diseases.

  Study EGAS00001007362

• APOBEC mutagenesis is a common process in normal human small intestine

  The single base substitution mutational signatures SBS2 and SBS13, likely caused by APOBEC cytosine deaminases, are common in many human cancer types. However, the stimulus activating APOBEC mutagenesis is unknown and understanding of when it occurs in the progression from normal to cancer cell is limited. Here, as part of a wider survey of human tissues, we whole genome sequenced 342 microdissected normal epithelial crypts from the small intestines of 39 individuals. SBS2/13 mutations were present in 17% normal small intestine crypts and were likely due to APOBEC3A activity. Localised clusters of SBS2/13 mutations (kataegis) were also commonly found. APOBEC mutation burdens were variable between individuals and between crypts from the same individual. Crypts with SBS2/13 often had immediate crypt neighbours without SBS2/13, suggesting that the underlying cause of SBS2/13 is cell-intrinsic rather than a widely distributed microenvironmental exposure, or needs to be permitted by cell-intrinsic conditions. APOBEC mutagenesis occurred throughout the human lifespan, including in young children, and was episodic with a small number of episodes occurring during the life history of a single cell. The results indicate that APOBEC mutagenesis is more common in the small intestine epithelium than in many other cell types, and is an episodic process in vivo initiated or permitted by cell intrinsic factors.

  Dataset EGAD00001008764

• The Haemgen RBC study

  Anaemia is a major determinant of global ill-health. To refine our understanding of the genetic factors influencing red blood cell formation and function, we carried out a meta-analysis of genome-wide association studies (GWAS) for six red blood cell traits: haemoglobin (HB), mean cell haemoglobin (MCH), mean cell haemoglobin concentration (MCHC), mean cell volume (MCV), packed cell volume (PCV) and red blood cell count (RBC). We provide genome-wide association results for 62,553 people of European ancestry using up to 2,644,161 autosomal SNPs. Participants with extreme measurements (>+/-3SD from mean) were excluded on a per phenotype basis. Imputation was done using haplotypes from HapMap Phase 2. SNP associations with each phenotype were tested by linear regression using an additive genetic model. Associations were tested separately in each cohort, with principal components and other study specific factors as covariates to account of population substructure. We then carried out meta-analysis of results from the individual cohorts using z-scores weighted by square root of sample size. SNPs with MAF<1% (weighted average across cohorts) were removed, as were SNPs with weight <50% of phenotype sample size. Anaemia is a major determinant of global ill-health. To refine our understanding of the genetic factors influencing red blood cell formation and function, we carried out a meta-analysis of genome-wide association studies (GWAS) for six red blood cell traits: haemoglobin (HB), mean cell haemoglobin (MCH), mean cell haemoglobin concentration (MCHC), mean cell volume (MCV), packed cell volume (PCV) and red blood cell count (RBC). We provide genome-wide association results for 62,553 people of European ancestry using up to 2,644,161 autosomal SNPs. Participants with extreme measurements (>+/-3SD from mean) were excluded on a per phenotype basis. Imputation was done using haplotypes from HapMap Phase 2. SNP associations with each phenotype were tested by linear regression using an additive genetic model. Associations were tested separately in each cohort, with principal components and other study specific factors as covariates to account of population substructure. We then carried out meta-analysis of results from the individual cohorts using z-scores weighted by square root of sample size. SNPs with MAF<1% (weighted average across cohorts) were removed, as were SNPs with weight <50% of phenotype sample size.

  Study EGAS00000000132

• Warm_autopsy__mutational_signatures_and_clonal_units

  Aims and objectives 1. To produce a ‘first generation map’ of mutational burdens and mutational processes that operate in normal human tissues. 2. To examine clonal units and their distribution across normal human tissues with the emphasis on these particular areas: a) Tissues with an architecturally defined unit: Are architectural structures populated by a progeny derived from a single stem cell? Are clonal units confined to a given architectural structure? Examples of tissues where such approach will be used include prostate gland (acini and ducts), small bowel (crypts), thyroid gland (follicles) and others. b) Tissues without an architecturally defined unit: What is the average size of clonal populations? This group includes stratified/pseudostratified epithelia with a basal layer within which there are stem cells that give rise to a progeny of cells showing upwards migration and lateral expansion (e.g. urothelium, cervical squamous epithelium, respiratory epithelium). This approach can also be used in tissues with a sheet-like cell arrangement, for example brain and adrenal gland. Methods To answer these questions, we will use a range of normal human tissue samples which have already been collected from a single deceased individual. The samples are currently stored in a -80C freezer. Ethanol fixed frozen sections will be prepared from these biopsies. Putative clonal units will be identified and dissected out using laser capture microscopy. DNA extraction will be performed using the in-house (CGP) protease based method. In terms of sequencing methods, the project will be split into 2 phases (2 different approaches): - Phase 1: This part of the project will utilise a combination of low input library preparation (developed by Peter Ellis) and standard whole genome sequencing (HiSeqX). The phase has already started (December 2016). This approach will be particularly useful in investigation of mitotically active tissues. - Phase 2: This part of the project will utilise a combination of low input library method (Peter Ellis) and Bottle-neck sequencing (BotSeqS). BotSeqS is a type of duplex sequencing method that offers a more accurate variant calling and allows genome interrogation at a single cell level. This method is currently being developed and is expected to be implemented in 3-6 months (March-April 2017). This approach will be particularly useful for post-mitotic tissues and tissues without pre-defined architectural units.

  Study EGAS00001002216

• SABE_1171_WGS_GVCF_IDX

  Dataset EGAD00001008203

• CRCbiome: Gut metagenome of Norwegian screening participants using FIT sampling

  The study aims to characterize the gut metagenome of Norwegian screening participants and assess the ability of the gut metagenome to identify those at need for colonoscopy. Background: Colorectal cancer (CRC) screening reduces CRC incidence and mortality. However, current screening methods are either hampered by invasiveness or suboptimal performance, limiting their effectiveness as primary screening methods. To aid in the development of a non-invasive screening test with improved sensitivity and specificity, we have initiated a prospective biomarker study (CRCbiome), nested within a large randomized CRC screening trial in Norway. We aim to develop a microbiome-based classification algorithm to identify advanced colorectal lesions in screening participants testing positive for an immunochemical fecal occult blood test (FIT). We will also examine interactions with host factors, diet, lifestyle and prescription drugs. The prospective nature of the study also enables the analysis of changes in the gut microbiome following the removal of precancerous lesions. Methods: The CRCbiome study recruits participants enrolled in the Bowel Cancer Screening in Norway (BCSN) study, a randomized trial initiated in 2012 comparing once-only sigmoidoscopy to repeated biennial FIT, where women and men aged 50–74 years at study entry are invited to participate. Since 2017, participants randomized to FIT screening with a positive test result have been invited to join the CRCbiome study. Self-reported diet, lifestyle and demographic data are collected prior to colonoscopy after the positive FIT-test (baseline). Screening data, including colonoscopy findings are obtained from the BCSN database. Fecal samples for gut microbiome analyses are collected both before and 2 and 12 months after colonoscopy. Samples are analyzed using metagenome sequencing, with taxonomy profiles, and gene and pathway content as primary measures. CRCbiome data will also be linked to national registries to obtain information on prescription histories and cancer relevant outcomes occurring during the 10 year follow-up period. Discussion: The CRCbiome study will increase our understanding of how the gut microbiome, in combination with lifestyle and environmental factors, influences the early stages of colorectal carcinogenesis. This knowledge will be crucial to develop microbiome-based screening tools for CRC. By evaluating biomarker performance in a screening setting, using samples from the target population, the generalizability of the findings to future screening cohorts is likely to be high.

  Study EGAS50000000170

• Small RNA sequencing of endometrial cancer patients with samples collected prior to diagnosis and controls

  This dataset contains unpaired FASTQ files of 316 endometrial cancer cases and 316 matched controls representing circulating RNAs. RNA count files are provided as the output of sncRNA pipeline and represent 12 RNA types (isomiR, lncRNA, precursor miRNA, miRNA, miscRNA, mRNA, piRNA, scaRNA, snoRNA, snRNA, tRF, and tRNA). We also provide metadata of the samples.

  Dataset EGAD50000000391

• Identification of improved IL28B SNPs and haplotypes for prediction of drug response in treatment of hepatitis C using massively parallel sequencing in a cross-sectional European cohort

  The hepatitis C virus infects nearly 3% of the World’s population, causing severe liver disease in many. Standard of care therapy is currently pegylated interferon alpha and ribavirin (PegIFN/R), which is effective in less than half of those infected with the most common viral genotype. Two IL28B SNPs, rs8099917 and rs12979860, predict response to (PegIFN/R) therapy in treatment of hepatitis C virus infection. These SNPs were identified in genome wide analyses using Illumina genotyping chips. In people of European ancestry, there are 6 common (>1%) haplotypes for IL28B, one tagged by rs8099917 minor allele, four tagged by rs12979860. We used massive parallel sequencing of the IL28B and IL28A gene regions generated by PCR from pooled DNA samples from 100 responders and 99 non-responders to therapy, to identify common variants. Variants that had high odds ratios and were validated were then genotyped in a cohort of 905 responders and non-responders. Their predictive power was assessed, alone and in combination with HLA-C. Only SNPs in the IL28B linkage disequilibrium block predicted drug response. Eighteen SNPs were identified with evidence for association with drug response, and with a high degree of confidence in the sequence call. We found that two SNPs, rs4803221 (homozygote minor allele positive predictive value (PPV) of 77%) and rs7248668 (PPV 78%), predicted failure to respond better than the current best, rs8099917 (PPV 73%) and rs12979860 (PPV 68%) in this cross-sectional cohort. The best SNPs tagged a single common haplotype, haplotype 2. Genotypes predicted lack of response better than alleles. However, combination of IL28B haplotype 2 carrier status with the HLA-C C2C2 genotype, which has previousl y been reported to improve prediction in combination with IL28B, provides the highest PPV (80%). The haplotypes present alternative putative transcription factor binding and methylation sites. Massive parallel sequencing allowed identification and comparison of the best common SNPs for identifying treatment failure in therapy for HCV SNPs tagging a single haplotype have the highest PPV, especially in combination with HLA-C. The functional basis for the association may be due to altered regulation of the gene. These approaches have utility in improving diagnostic testing and identifying causal haplotypes or SNPs.

  Study EGAS00001000096

• Mutational dynamics of triple negative breast cancer over neoadjuvant chemotherapy treatment reveal frequent whole genome duplication events

  Background Triple-negative breast cancer (TNBC) is associated with poor survival rate and high genomic instability, generating complex tumour genomes. However, the processes generating this complexity are poorly studied over time. Here, we study the temporal dynamics of TNBC somatic mutations, revealing major transitions in tumour genome evolution, from diagnostic biopsies, through treatment, to cancer remission or recurrence. Methods Deep whole exome sequencing and CUTseq, a reduced representation whole genome sequencing approach, were performed in parallel, to comprehensively identify short nucleotide variants (SNVs), copy number alterations (CNAs) and aneuploidies genome-wide. We profiled (N=74) tumour samples from 22 patients across the course of disease, including spatially diverse samples from many tumours, allowing us to observe the gain and loss of candidate driver variants over time. Results As expected, tumour infiltrating lymphocyte (TIL) levels and residual cancer burden (RCB) score were negatively correlated, with good responses to neoadjuvant chemotherapy (NACT) seen in pre-treatment samples with high infiltration. In general, SNV mutational burden remained stable across disease progression and RCB classes. Recurrent SNVs across the cohort implicated several known TNBC driver genes in disease progression and response to treatment, with TP53, MICA, CYP2D6, BRCA1 and BRCA2 frequently altered. However, the presence of candidate driver SNVs in these genes often varied throughout a patient’s treatment, including the loss of those seen in a given driver gene in pre-treatment samples, but gain of novel variants in the same gene at later time points. Dramatic structural dynamics were seen in all tumours, including abundant CNA including chromosome arm aneuploidies seen in pre-treatment samples, followed by frequent loss of these alterations post-treatment, and their re-emergence at recurrence. Whole genome duplication (WGD) events appear to drive these dynamics, with a higher frequency of pre-treatment WGD seen in patients with the best response to NACT. Conclusions Comprehensive longitudinal profiling demonstrates the complex interplay of SNVs and copy number alterations during TNBC progression, leading to diverse evolutionary trajectories and patient outcomes. Complex mutational patterns emerge as important features during progression, with WGD events emerging as potential candidate biomarkers of response at both tumour establishment and recurrence.

  Study EGAS00001008261

• Sensitive and robust liquid biopsy-based detection of PIK3CA mutations

  Over the last years, the treatment of hormone receptor-positive metastatic breast cancer patients has dramatically changed. Combination strategies attempting to overcome resistance to endocrine therapy such as targeting the PI3 kinase pathway are gaining importance. The phase III SOLAR-1 trial demonstrated the benefit of the addition of alpelisib to fulvestrant, thereby providing further clinical evidence of the increasing importance of PIK3CA testing. Here, we performed a comparison of liquid biopsy and tissue-based detection of PIK3CA mutations in HR-positive metastatic breast cancer patients. Plasma samples and the latest available tumor tissue from 69 patients with metastatic hormone-receptor positive breast cancer were analyzed for PIK3CA hotspot mutations. A high-resolution NGS assay (SiMSen-Seq) was used for plasma samples covering 11 recurrent PIK3CA mutations, used for stratification in SOLAR-1. Additionally, mFAST-SeqS was used to estimate the tumor fractions in the plasma samples. For tumor tissue samples, targeted Ion Torrent NGS was employed. Matched tissue and plasma samples were available from 63 patients. PIK3CA mutations were detected in 46.1% of the tissue samples and in 50.7% of the plasma samples. The variant allele frequency (VAF) ranged from 4.36% to 72.93% with a median of 28.60% in tissue samples, and from 0.28% to 49.85% with a median of 1.17% in plasma samples. The most frequent variants in tissue and plasma were H1047R, E545K, and E542K. We detected at least one of 11 PIK3CA mutations from the SOLAR-1 trial in 51/63 patients, both in tissue and plasma. In 7 patients, such mutations were found only in plasma. In 5 patients, PIK3CA mutations found in tissue were not detectable in ctDNA, 2 of which had a low tumor fraction.Together, in 33/67 plasma samples without detectable PIK3CA mutations, 18 samples had elevated tumor fraction, implicating true negative results. In the remaining 15 patients with low tumor fraction, additional tissue analysis was needed. SiMSen-Seq-based detection of PIK3CA mutations in plasma shows advantageous concordance with the tissue analyses. A combination with an untargeted, mutation-independent approach for detecting ctDNA fractions may confirm a negative PIK3CA result and thereby enhance the performance of the SiMSen-Seq test. This combinatory approach for selecting suitable patients for alpelisib treatment should be validated in larger patient cohorts.

  Study EGAS00001004940

• Dissociation of solid tumour tissues with cold active protease for single-cell RNA-seq minimizes conserved collagenase-associated stress responses

  Background Single-cell RNA sequencing (scRNAseq) is a powerful tool for studying complex biological systems, such as tumour heterogeneity and tissue microenvironments. However, the sources of technical and biological variation in primary solid tumour tissues and patient-derived mouse xenografts for scRNAseq, are not well understood. Here, we used low temperature (6°C) protease and collagenase (37°C) to identify the transcriptional signatures associated with tissue dissociation across a diverse scRNAseq dataset comprising 128,481 cells from patient cancer tissues, patient-derived breast cancer xenografts and cancer cell lines.Results We observe substantial variation in standard quality control (QC) metrics of cell viability across conditions and tissues. From FACS sorted populations gated for cell viability, we identify a sub-population of dead cells that would pass standard data filtering practices, and quantify the extent to which their transcriptomes differ from live cells. We identify a further subpopulation of transcriptomically “dying” cells that exhibit up-regulation of MHC class I transcripts, in contrast with live and fully dead cells. From the contrast between tissue protease dissociation at 37°C or 6°C, we observe that collagenase digestion results in a stress response. We derive a core gene set of 512 heat shock and stress response genes, including FOS and JUN, induced by collagenase (37°C), which are minimized by dissociation with a cold active protease (6°C). While induction of these genes was highly conserved across all cell types, cell type-specific responses to collagenase digestion were observed in patient tissues. We observe that the yield of cancer and non-cancer cell types varies between tissues and dissociation methods.Conclusions The method and conditions of tumour dissociation influence cell yield and transcriptome state and are both tissue and cell type dependent. Interpretation of stress pathway expression differences in cancer single cell studies, including components of surface immune recognition such as MHC class I, may be especially confounded. We define a core set of 512 genes that can assist with identification of such effects in dissociated scRNA-seq experiments.

  Study EGAS00001003753

• Molecular_characterization_of_invasive_lobular_carcinoma

  Invasive lobular carcinoma (ILC) is the second most common histological subtype of breast cancer accounting for 10-15% of cases. ILC differs from invasive ductal carcinoma (IDC)with respect to epidemiology, histology, and clinical presentation. Moreover, ILC is lesssensitive to chemotherapy, more frequently bilateral, and more prone to form gastrointestinal, peritoneal, and ovarian metastases than IDCs. In contrast to IDC, the prognostic value ofhistological grade (HG) in ILC is controversial. One of the three major components of histological grading (tubule formation) is missing in ILC which hinders the process of gradingin this histological subtype and results in the classification of approximately two thirds of ILC as HG 2.Over the last decade, a number of gene expression signatures have shed light onto breast cancer classification, allowing breast cancer care to become more personalized. Withrespect to the management of estrogen receptor (ER)-positive breast cancer, several gene expression signatures provide prognostic and/or predictive information beyond what is possible with current classical clinico-pathological parameters alone. Nevertheless, most studies using gene expression signature have not considered different histologic subtypesseparately. Recently, a comprehensive research program has elucidated some of the biological underpinnings of invasive lobular carcinoma. Genetic material extracted from 200 ILC tumor samples were studied using gene expression profiling and identified ILCmolecular subtypes. These proliferation-driven gene signatures of ILC appear to have prognostic significance. In particular, the Genomic Grade (GG) gene signature improved upon HG in ILC and added prognostic value to classic clinico-pathologic factors. In addition this study demonstrated that most ILC are molecularly characterized as luminal-A (~75%)followed by luminal-B (~20%) and HER2-positve tumors (~5%). Moreover, we investigated the prognostic value of known gene signatures/ gene modules in the same cohort of ILC. As a second step within the scope of this project, we aim to investigate the interactionsbetween somatic ILC tumor mutations to observed transcriptome findings. To this end, we aim to perform somatic mutation analysis for the ILC tumors for which Affymetrix gene expression profiling is available. To this end, we will use a gene screen assay, which specifically interrogates the mutational status of a few hundreds of cancer genes. We believe that this pioneering effort will be fundamental for a tailored treatment of ILC withimprovement in patients' outcome.

  Study EGAS00001000292

• Genome_Diversity_in_Africa_Project___GemCode_libraries

  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 Africa 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 Africa would facilitate genomic epidemiological studies in the region.The Genome Diversity in Africa Project (GDAP) will extend and expand the African Genome Variation (AGV) project, an international collaboration aimed to produce a comprehensive catalogue of human genetic variation in Sub-Saharan Africa (SSA). Using a sequencing-based approach, GDA project aims to capture human genetic variation in Africa, 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 Africa, as well as providing a global resource to help design, implement and interpret genomic studies in African populations and studies comprising globally diverse populations, complementing existing genomic resources. Specifically, we plan to carry out whole genome sequencing at high depth (30x) of up to 1000 individuals across Africa (25 individuals from each ethnolinguistic group). We have already completed sequencing from 5 ethno-linguistic groups and we are now adding others.Our scientific objectives are to: 1) develop a resource that provides a comprehensive catalogue of genetic variation in populations from Africa accessible to the global scientific community; 2) characterise population genetic diversity, structure, gene flow and admixture across Africa; 3) develop a cost-efficient, next-generation genotype array for diverse populations across Africa; 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 EGAS00001001589

• Clinical and Molecular Investigation of Familial CEBPA-mutated Acute Myeloid Leukaemia

  Myelodysplastic syndromes (MDS) and acute myeloid leukaemia (AML) are clonal diseases arising from the sequential acquisition of genetic aberrations in haematopoietic stem cells. Whilst the aetiology is predominantly sporadic, rare reports describe the autosomal dominant inheritance of MDS and AML (familial MDS/AML), often presenting in younger patients (&lt;40 years of age) with a germline mutation in one of three haematopoietic transcription factors: RUNX1, CEBPA and GATA2. The clinical recognition of familial disease can be challenging due to wide variations in phenotype, latency and transmission, complicated further by the patients’ awareness and reporting of affected relatives. Distinct clinical manifestations have however been reported for each genetic subtype and whilst RUNX1 and GATA2 mutations can harbour prolonged dysplastic or cytopenic phases, germline CEBPA mutations universally present with AML de novo. Our group were the first to describe germline CEBPA mutations in 3 related individuals, all of whom developed AML between 10-30 years of age. The initial pedigree together with those described in the intervening period, were predominantly isolated case reports with limited follow up, hence long term outcomes for this group of patients remain poorly defined. This contrasts with CEBPA mutations occurring in 10-15% of sporadic AML with normal karytoype (NK-AML), where double mutations (CEBPAdm, found in approximately 50% of all CEBPA-AML) are considered a unique molecular entity, associated with favourable clinical outcomes. In sporadic and familial CEBPAdm AML, N’ terminal mutations enforce translation of the shorter P30-kDa isoform and are typically accompanied by C’ terminal insertions or deletions (disrupting the DNA-binding and/or Leucine zipper domains). Definition of the full molecular landscape awaits extensive whole exome sequencing, although GATA2 mutations (targeting the Zinc Finger 1, ZF1, domain) have been reported in 20-40% of sporadic CEBPAdm AML and recent targeted resequencing analysis revealed a further association with WT1 mutations. With minimal clinical and molecular information available for familial CEBPA-AML, we are often consulted regarding the management of this unique patient population. To address this need, we have investigated a multi-centre cohort of families with germline CEBPA mutations, providing extended clinical follow up, whole exome profiling of secondary mutations and molecular analysis of relapsed disease.

  Study EGAS00001000949

• Conjunctival fibrosis and the innate barriers to Chlamydia trachomatis intracellular infection: a genome wide association study

  This data comes from a genome-wide association study in 1090 cases of scarring trachoma and 1531 controls from a non-selected multi-ethnic sampling of The Gambia. The data constitutes hard base calls at 1457295 features from the HumanOmni2.5-8 (Omni2.5) BeadChip. Three genotype-calling algorithms were used : Illuminus (Bioinformatics 2007, 23:2741–6), GenCall (Patent US7035740, 2006) and GenoSNP (Bioinformatics 2008, 24:2209–14). Data from each calling algorithm was filtered to retain only SNPs with a call rate ≥ 0.98. To obtain a merged set of SNPs, all genotypes that matched across call-sets were retained, whilst those that mismatched between call-sets were set to missing. Genotypes that were present in one call set and missing in others were also retained. The merged SNP set contained 1467876 SNPs, which were filtered for call rate ≥ 0.99 and Hardy-Weinberg equilibrium P value &lt; 5 x 10-8. The case and control groups were approximately equivalent with respect to gender and ethnicity. 70% of cases and 63% of controls were female. Self-described ethnic composition of the cases was 29% Jola, 27% Mandinka, 21% Wolof, 10% Fula and 13% other/no data. In controls the composition was 25% Jola, 24% Mandinka, 21% Wolof, 8% Fula and 22% other/no data. Median age was 49 (range 32-60) in cases and 37 (range 12-52) in controls (t = -3.6019, P = 0.0003). This data set has substantial levels of population structure, which reflects extensive networks of pairwise kinship between participants. First pass tests of association in EMMAX (Nat Genet 2010, 42:348–54) led to genome-wide deflation of the test statistics (λ = 0.982 SE = 1.74e-05) that could be directly attributed to the differing age and gender distributions between cases and controls. Modelling the phenotype to adjust for age and gender successfully controlled for this and no genome-wide deviation from the null was subsequently observed (λ = 1.001, SE = 8.7 x 10-7). Neither principal components analysis (PCA) nor tests of proportional identity by state (IBS) variance in PLINK identified significant levels of within or between group genetic variance. The phenotypes submitted in the phenotypes file include both uncorrected (case/control status) and corrected (age and gender corrected case/control status) phenotypes, as well as age and self-described ethnicity.

  Study EGAS00001001516

• Genome_Diversity_in_Africa_Project__Benin

  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 bioinformaticsThis 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 EGAS00001001688

• Multi-omics data of 1000 Inflammatory Bowel Disease patients

  Portal available at https://1000ibd.org Inflammatory bowel disease (IBD) is a chronic complex disease of the gastrointestinal (GI) tract. Patients with IBD can experience a wide range of symptoms, but the pathophysiological mechanisms that cause these individual differences in clinical presentation remain largely unknown. In consequence, IBD is currently classified into subtypes using clinical characteristics. If we are to develop a more targeted treatment approach, molecular subtypes of IBD need to be discovered that can be used as new drug targets. To achieve this, we need multiple layers of molecular data are generated from the same IBD patients.We initiated the 1000IBD project to prospectively follow more than 1000 IBD patients from the Northern provinces of the Netherlands. For these patients, we have collected a uniquely large number of phenotypes and generated multi-omics profiles. To date, 1,215 participants have been enrolled in the project and enrolment is on-going, with 609 patients being present in the first data release. Phenotype data collected for these participants includes information on dietary and environmental factors, drug responses, and adverse drug events. Genome information has been generated using genotyping (ImmunoChip, Global Screening Array and HumanExomeChip) and sequencing (whole exome sequencing and targeted resequencing of IBD susceptibility loci), transcriptome information generated using RNA-sequencing of intestinal biopsies and microbiome information generated using both sequencing of the 16S rRNA gene and whole genome shotgun metagenomic sequencing.All molecular data generated within the 1000IBD project will be shared in multiple data on the European Genome-Phenome Archive (www.ega-archive.org). The first data release, will contain basic phenotypes for 609 participants, genotypes of 314 participants, gut microbiome data generated by tag sequencing the 16S gene of feces from 315 participants and intestinal biopsies from 107 participants, and gut microbiome metagenomic sequencing of feces of 355 participants . Future releases will comprise many more additional phenotypes and -omics data layers. 1000IBD data can be used by other researchers as a replication cohort, a dataset to test new software tools, or a dataset for applying new statistical models.

  Study EGAS00001002702

• UK10K RARE COLOBOMA

  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. Ocular coloboma is the most common significant developmental eye defect with an incidence of ~1 in 5,000 live births. It results from failure of optic fissure closure during embryogenesis. The position and extent of the fusion failure dictates the clinical appearance and functional effect. ~30% of coloboma cases are associated with other systemic malformations. These UK10K samples will mostly comprise isolated coloboma cases without systemic involvement (aka non-syndromal coloboma). There is strong evidence from family studies that coloboma has a major genetic component with autosomal dominance being the most common pattern of inheritance. However, many cases are isolated or show complex patterns of familial clustering. The genes responsible for isolated coloboma are largely unknown, but in a small number of families mutations in SHH, CHX10, and PAX6 have been identified indicating marked genetic heterogeneity. Thus in addition to the clinical benefits of achieving a molecular diagnosis there are also major scientific advantages to identifying coloboma genes, as these are likely to provide insights into the complex process of optic fissure closure, that is critical to normal eye development. In the longer term, understanding the molecular basis of the disease may provide clues to therapeutic strategies.For further information with regard to this cohort please contact David Fitzpatrick (david.fitzpatrick@ed.ac.uk).

  Study EGAS00001000127

• UK10K_OBESITY_SCOOP

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches.The SCOOP samples are part of the Obesity group and will undergo exome sequencing. Severe Childhood Onset Obesity Project (SCOOP) is a sub-cohort of the Genetics Of Obesity Study (GOOS) cohort established by Sadaf Farooqi and Steve O’Rahilly at the University of Cambridge over the last 12 years. The GOOS cohort contains &gt;4,000 patients of diverse geographic origin, many of whom have monogenic and syndromic forms of obesity, and includes patients that are offspring of consanguineous union. SCOOP is a subset of &gt;1500 UK Caucasian patients with severe, early onset obesity (all patients have a BMI Standard Deviation Score (SDS) &gt; 3 and obesity onset before the age of 10 years), in whom all known monogenic causes of obesity have been excluded. GWAS data on the SCOOP cohort will be available (WTCCC2 independent study) at the time of the start of this study. Data from this cohort has demonstrated that the prevalence of the common obesity risk alleles in FTO, MC4R and NEGR1 are amongst the highest within SCOOP, demonstrating its value also in the study of genetic variants with an impact on more common obesity. For further information with regard to this cohort please contact Sadaf Farooqi (isf20@medschl.cam.ac.uk).

  Study EGAS00001000124

• Clinical outcomes in ctDNA-positive urothelial carcinoma patients treated with adjuvant immunotherapy

  Minimally invasive approaches to detect residual disease after surgery are urgently needed to select patients at highest risk for metastatic relapse for additional therapies. Circulating tumour DNA (ctDNA) holds promise as a biomarker for molecular residual disease (MRD) and relapse,1-3 but its clinical value has yet to be demonstrated in a randomised clinical trial. We evaluated outcomes in post-surgical ctDNA-positive (+) patients in a randomised phase III trial of adjuvant atezolizumab versus observation. IMvigor010 enrolled 809 patients with muscle-invasive urothelial carcinoma and did not meet its primary endpoint of disease-free survival (DFS) in the intent-to-treat population. Within the study, an exploratory planned analysis of prospectively collected plasma was performed, which tested the utility of ctDNA to identify patients who may benefit from adjuvant atezolizumab treatment. ctDNA was measured at the start of therapy (cycle 1 day 1; C1D1) and at week 6 (cycle 3 day 1; C3D1), and 581 patients were evaluable for ctDNA. The prevalence of ctDNA positivity at C1D1 was 37% (n=214), and ctDNA positivity identified patients with poor prognosis (observation arm DFS HR= 6.19 (4.29, 8.91), p&amp;lt;0.0001). Here we show that ctDNA(+) patients had improved DFS and overall survival (OS) with atezolizumab versus observation (DFS HR= 0.56 (0.41-0.77); p=0.0003 and OS HR= 0.58 (0.4-0.86); p=0.0063). No difference in DFS or OS between arms was noted for ctDNA-negative patients. The rate of ctDNA clearance was higher with atezolizumab (18%) versus observation (4%) (p=0.0041). Transcriptomic analysis revealed that tumours from ctDNA(+) patients had higher expression of cell cycle and keratin genes. Within the ctDNA(+) patient population in the atezolizumab arm, non-relapsing patients were further enriched in prominent immune response signatures including PD-L1, IFNG, CXCL9, and high tumour mutational burden, whereas relapse was associated with angiogenesis and fibroblast-transforming growth factor-beta signatures (F-TBRS). TCGA molecular subset analysis revealed increased efficacy of atezolizumab in patients with basal-squamous tumours, consistent with underlying tumour-immune contexture. Together these findings suggest that adjuvant atezolizumab may be associated with improved outcomes compared with observation in this high-risk ctDNA(+) population. These findings, if validated in other settings, would shift approaches to post-operative cancer care.

  Study EGAS00001004997

• RNA_Seq_in_Patients_with_Primordial_Dwarfism

  Transcriptome studies in patients with rare genetic diseases can potentially aid in theinterpretation of likely causal genetic variation through identification of altered transcriptabundance and/or structure. RNA-Seq is the most sensitive assay for both investigatingtranscript structure and abundanceThe primary aim of this pilot project is to investigate to what degree integrating exome-Seqand RNA-Seq data on the same individual can accelerate the identification of causal allelesfor rare genetic diseases. There are two main strands to this: (i) identifying which variantsdiscovered in exome-seq appear to be having a functional impact on transcripts, and (ii)identifying transcript outliers, especially among known causal genes, that may not necessarilyhave a causal variant identified from exome sequencing. The latter may identify the presenceof causal variants that lie far from coding regions (e.g. the formation of cryptic splice sitesdeep within introns, or loss of long range regulatory elements), which can be confirmed withfurther targeted genetic assays. Just over 50% of all disease-causing variants recorded in theHuman Gene Mutation Database (HGMD) affect transcript structure and abundance (e.g.nonsense SNVs, essential splice site SNVs, frameshifting indels, CNVs).This pilot project will study RNA from lymphoblastoid cell-lines from 12 patients withprimordial dwarfism syndromes, for 10 of these samples we have previously generate exomedata as part of our collaboration with the group of Prof Andrew Jackson. The two remainingsamples are positive controls where the causal mutation is known, and is known to affecttranscript structure and/or abundance.Primordial dwarfism is a prime candidate for these RNA-seq studies because all knowncausal mutations to date have key roles in DNA replication and thus, unsurprisingly, theproducts of the causal genes are typically ubiquitously expressed.Each RNA will be sequenced, with two technical replicates (independent RT-PCR and libraries) persample, and each replicate run in 1/2 of a HiSeq lane using 100bp paired reads. Samples preparation was as follows :The cells were grown to confluency, then pellets frozen at -80. RNA samples were prepared using the Qiagen RNeasy kit, then nanodropped and analyzed using the bioanalyzer to determine concentration and purity.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 EGAS00001000283

• Defective Homologous Recombination DNA Repair as Therapeutic Target in Advanced-Stage Chordoma (HIPO_021)

  HIPO project: HIPO_021 Importance: Chordomas are rare tumors of the axial skeleton and skull base with few therapeutic options and no clinically validated molecular drug targets. The value of comprehensive genomic analyses for guiding medical therapy of patients with advanced-stage chordoma is unknown.Objective: To identify therapeutically tractable genetic lesions in a cohort of chordoma patients within a genomics-guided precision oncology program and to document the outcome of individualized, molecularly targeted chordoma therapy.Design, Setting, and Participants: We performed whole-exome sequencing of tumor and matched germline control samples from seven patients with locally advanced or metastatic chordoma who were enrolled in a cross-institutional molecular stratification registry trial for younger adults with advanced-stage cancer across all histologies and patients with rare tumors. All patients were heavily pretreated and had progressive disease prior to molecular analysis.Interventions: Individualized medical therapy was administered according to the patients’ molecular profiles.Main Outcomes and Measures: Candidate therapeutic targets identified by whole-exome sequencing and response to genotype-directed therapy.Results: All patients harbored alterations of two or more genes known to be involved in DNA repair via homologous recombination (HR), including heterozygous deletions of ERCC6, FANC family members, RAD51L (n = 6), BRCA2 (n = 5), ATR, CHEK2, RAD18, RAD51B, and XRCC3 (n = 4); inactivating PTEN mutations coupled with loss of heterozygosity (n = 2); and pathogenic germline variants in BRCA2 (n = 1), NBN (n = 1), and CHEK2 (n = 1) that were accompanied by somatic deletion of the corresponding wild-type alleles. Consistently, a mutational signature associated with defective HR was enriched in all samples and co-occurred with extensive genomic instability, as evidenced by HR deficiency scores and high numbers of large-scale state transitions. These results prompted off-label treatment with the poly(ADP-ribose) polymerase (PARP) inhibitor olaparib in a patient whose tumor was refractory to irradiation and systemic treatment with imatinib, which led to a prolonged response and substantial clinical improvement.Conclusions and Relevance: Advanced-stage chordomas are frequently characterized by genomic imprints of defective HR DNA repair. HR deficiency represents a new therapeutic opportunity in this intractable disease through repositioning of PARP inhibitors that warrants further exploration in clinical trials.

  Study EGAS00001002720

• RNA-seq of DF149 cells – a patient-derived xenograft model of ascites-derived, homologous recombination repair-proficient, high-grade serous ovarian carcinoma – cultured in vitro and isolated after 8 hours treatment with DMSO control (3 x biological replicates) and 2.5 µM CBL0137 (3 x biological replicates)

  Background: High-grade serous ovarian carcinomas (HGSCs) are a heterogeneous subtype of epithelial ovarian cancers and include serous cancers arising in the fallopian tube and peritoneum. These cancers are now subdivided into homologous recombination repair (HR)-deficient and proficient subgroups as this classification impacts on management and prognosis. PARP inhibitors (PARPi) have shown significant clinical efficacy, particularly as maintenance therapy following response to platinum based chemotherapy in BRCA-mutant or homologous recombination (HR)-deficient HGSCs in both the 1st and 2nd line settings. However, PARPi have limited clinical benefit in HR-proficient HGSCs which make up almost 50% of HGSC and improving outcomes in these patients is now a high priority due to the poor prognosis with current standard of care. There are a number of potential lines of investigation including efforts in sensitizing HR-proficient tumors to PARPi. Herein, we aimed to develop a novel combination therapy by targeting SSRP1 using a small molecule inhibitor CBL0137 with PARPi in HR-proficient HGSOCs. Experimental design: We tested anti-cancer activity of CBL0137 monotherapy using a panel of HGSOC cell lines and patient-derived tumor cells in vitro. RNA sequencing was used to map global transcriptomic changes in CBL0137-treated patient-derived HR-proficient HGSOC cells. We tested efficacy of CBL0137 in combination with PARPi using HGSOC cell lines and patient-derived tumor cells in vitro and in vivo. Results: We show that SSRP1 inhibition using a small molecule, CBL0137, that traps SSRP1 onto chromatin, exerts a significant anti-growth activity in vitro against HGSC cell lines and patient-derived tumor cells, and also reduces tumor burden in vivo. CBL0137 induced DNA repair deficiency via inhibition of the HR repair pathway and sensitized SSRP1-high HR-proficient HGSC cell lines and patient-derived tumor cells/xenografts to the PARPi, Olaparib in vitro and in vivo. CBL0137 also enhanced the efficacy of DNA damaging platinum-based chemotherapy in HGSC patient-derived xenografts. Conclusion: Our findings strongly suggest that combination of CBL0137 and PARP inhibition represents a novel therapeutic strategy for HR-proficient HGSCs that express high levels of SSRP1 and should be investigated in the clinic.

  Study EGAS00001006662

• Exploring high-throughput drug sensitivity testing in neuroblastoma cell lines and patient-derived tumor organoids in the era of precision medicine

  Despite druggable events in 80% of neuroblastoma patients within the Princess Máxima Center precision medicine program 'iTHER', clinical uptake of treatment recommendations has been low, and the clinical impact for individual patients remains hard to predict. This stresses the need for a method integrating genomics and transcriptomics with functional approaches into therapeutic decision making. We aimed to launch an online repository integrating genomics and transcriptomics with high-throughput drug screening (HTS) of nineteen commonly used neuroblastoma cell lines and fourteen generated neuroblastoma patient-derived organoids (NBL-PDOs) to improve identification of molecularly matched therapies and support clinical uptake. Cell lines, NBL-PDOs and their parental tumors were characterized utilizing (lc)WGS, WES and RNAseq. Cells were exposed to a ~200 compound library, and viability was assessed using CellTiterGlo or MTT. Molecular and HTS results were transferred to the R2 platform (http://r2platform.com/pmc_nb_drugs/), in which a comprehensive suite of visualizations and analysis options were implemented. As a result, a powerful reference set of cell lines is available, reflecting distinct known pharmacologic vulnerabilities, including resistance of TP53-inactivated lines to MDM2 inhibitor idasanutlin; ALK-inhibitor sensitivity in ALK-activated lines, and sensitivity of an ATM-deleted cell line to PARP- and CHEK1-inhibitors. Novel established NBL-PDOs retained molecular features of their parental tumor and again reflected established drug sensitivities. However, HTS also identified additional therapeutic vulnerabilities in vitro, such as a striking correlation between a positive mesenchymal signature and sensitivity to BCL2-inhibitor venetoclax. Finally, we explored personalized drug sensitivities within iTHER, demonstrating HTS can support genomic and transcriptomic results, thereby strengthening the rationale for clinical uptake. Furthermore, in vitro insensitivity may avoid ineffective treatments. Lastly, novel treatment options may be identified, as indicated by sensitivity to ponatinib correlating with overexpression of RET or KIT. In conclusion, we established a dynamic publicly available dataset with detailed genomic, transcriptomic, and pharmacological annotation of classical neuroblastoma cell lines as well as NBL-PDOs, representing the heterogeneous landscape of neuroblastoma. We anticipate that in vitro drug screening will be complementary to genomic-guided precision medicine by supporting clinical decision making, thereby improving prognosis for all neuroblastoma patients in the future.

  Study EGAS00001007160

• Discrepancies in Tumour Mutation Burden (TMB) reporting from sequential Endobronchial ultrasound trans bronchial needle aspiration (EBUS TBNA) samples within single lymph node stations

  Introduction: Tumour Mutation Burden (TMB) is a potential biomarker for immune therapies. Identifying factors that may affect TMB prediction is key to its utility as biomarker for treatment options. Here we investigated parameters that might affect TMB using cytology smears obtained from endobronchial ultrasound transbronchial needle aspiration (EBUS TBNA)-sampled malignant lymph nodes, which are often the only diagnostic samples collected in patients with advanced lung cancer. Methods: Individual Diff-Quik cytology smears were prepared for each needle pass from the same enlarged malignant mediastinal lymph nodes. DNA extracted from smears of each pass underwent sequencing using a large gene panel sequencing (TSO500, Illumina). TMB was estimated for each individual needle pass using the TSO500 Local App v. 2.0 (Illumina). Results: Twenty patients had two or more Diff-Quik smears (total 45 smears) which passed sequencing quality control. Average TMB for all smears was 8.7 ± 5.0 mutations per megabase (Mb). Sixteen of the 20 patients had paired samples with minimal differences in TMB score (average difference 1.3 ± 0.85). Paired samples from 13 patients had concordant TMB in terms of scores being below or above a threshold of 10 mutations/Mb. Samples from seven patients had discrepant TMB, with four cases in particular exhibiting an average difference of 11.3 ± 2.7 mutations/Mb. Factors affecting TMB calling involved tumour content of the samples, the amount of DNA used in sequencing as well as bone fide intra tumour heterogeneity between paired samples. Conclusions: TMB assessment is feasible from EBUS-TBNA derived Diff-Quik slides from a single needle pass. Repeated samples of one or more lymph nodes have minimal variation in TMB in most cases (80%), suggesting accuracy of the test overall. To enhance accuracy of the result, care should be taken in relation to the tumour content of the cytology smear. Further, clinicians should be aware that even slight differences in lymph node site selection could give rise to variable genomics and discrepant TMB due to tumour heterogeneity. Therefore, it would be reasonable for multi-lymph node site analysis of TMB to be performed using EBUS TBNA aspirates from different easily accessible lymph nodes. Further investigations need to determine whether these samples should be combined as one sample or tested separately.

  Study EGAS00001007484

• Plasma MicroRNA Signatures of Aging

  Background. MicroRNAs are small non-coding RNAs that regulate gene expression post-transcriptionally and show differential expression in various tissues with aging phenotypes. Detectable in circulation, extracellular microRNAs reflect (patho)physiological processes and hold promise as biomarkers for healthy aging and age-related diseases. This study aimed to explore plasma extracellular microRNAs as biological aging indicator and their associations with health outcomes using population-level data. Methods. We quantified plasma expression levels of 2,083 extracellular microRNAs using targeted RNA-sequencing in 2,684 participants from the population-based Rotterdam Study cohort. The training and test sets included 1,930 participants from the advanced-aged initial and second subcohort (RS-I/RS-II; median age: 70.6), while the validation set comprised 754 participants from the middle-aged fourth subcohort (RS-IV; median age: 53.5). Based on 591 microRNAs well-expressed in plasma, we examined differential expression of microRNAs with chronological age, PhenoAge –a composite score of age and nine multi-system blood biomarkers–, the frailty index, and mortality. Next, elastic net models were employed to construct composite microRNA-based aging biomarkers predicting chronological age (mirAge), PhenoAge (mirPA), frailty index (mirFI), and mortality (mirMort). The association of these aging biomarkers with different age-related health outcomes was assessed using Cox Proportional Hazard, linear regression, and logistic regression models in the test and validation sets. Results. We identified 188 microRNAs differentially expressed with chronological age within the RS-I/RS-II advanced-aged population (ntraining=1158, ntest=772), of which 177 microRNAs (94.1%) were replicated in the middle-aged RS-IV subcohort (nvalidation=754). Moreover, 227 miRNAs showed robust associations with PhenoAge, 61 with FI, and 16 with ten-year mortality independent of chronological age. Subsequently, we constructed four plasma microRNA-based aging biomarkers: mirAge with 108, mirPA with 153, mirFI with 81, and mirMort with 50 miRNAs. Elevated scores on these microRNA-based aging biomarkers were associated with unfavorable health outcomes, including lower subjective physical functioning and self-reported health and increased mortality and frailty risk, but not with first- or multi-morbidity. Overall, larger effect estimates were observed for mirPA, mirFI, and mirMort compared to mirAge. Conclusions. This study describes distinct plasma microRNA-aging signatures and introduces four microRNA-based aging biomarkers with potential to identify accelerated aging and age-related decline, providing insights into the intricate process of human aging.

  Study EGAS00001008117

• Whole Genome Sequencing and Variant Calling for Autism Families

  We extracted DNA from whole blood or lymphoblast-derived cell lines and assessed the DNA quality with PicoGreenTM and gel electrophoresis. Whole genome sequencing was performed (Illumina HiSeq2000 and Illumina HiSeq X). WGS reads were mapped to the human reference genome assembly hg19 (GRCh37) using Burrows-Wheeler Aligner v.0.7.12 (TCAG) or Isaac v.2.0.13 (Macrogen). For each genome, we performed local realignment and quality recalibration and detected SNVs and small indels using GATK Haplotype Caller v.3.4.6 without genotype refinement. We detected CNVs using ERDS (estimation by read depth with single nucleotide variants) and CNVnator. We detected structural variants using Manta v.0.29.6. When available by the variant caller (i.e. GATK and Manta), trio-based joint variant calling was conducted for each family.

  Dataset EGAD00001008452

• A body map of somatic mutagenesis in morphologically normal human tissues

  This dataset contained raw sequencing fastq data of the article "A body map of somatic mutagenesis in morphologically normal human tissues". We sampled morphologically normal tissue biopsies from 5 donors. We performed low-depth WGS on 1,764 samples, high-depth WGS on 48 samples and WES on 1,772 samples.

  Dataset EGAD00001007859

• Transcriptome of CD4+ T cells and CD8+ T cells in glioblastoma multiforme

  Patients included in this study were over 18 years of age and had a histology-confirmed diagnosis of glioblastoma multiforme (GBM). Exclusion criteria were the previous administration of any anti-tumor therapy including radiation therapy. All patients gave written informed consent. The study was approved by the local ethics committee (TUM Medical school) and conducted following the Declaration of Helsinki. During resection of the tumors, tumor tissue and tissue from normal appearing brain within the operative channel was collected. Blood was drawn during the surgical procedure. Single cell suspensions were prepared from the tumor tissue, the normal appearing brain, and the blood. CD4+ T cells and CD8+ T cells were sorted by flow cytometry. Only patients with a complete set of specimens (CD4+ tumor infiltrating lymphocytes (TIL), CD8+ TIL, CD4+ T cells from normal appearing brain, CD8+ T cells from normal appearing brain, blood-derived CD4+ and CD8+ T cells) containing a minimum of 1000 cells in each sorted sample were further analyzed (n=9). Total RNA was isolated from sorted cell populations using the RNAeasy Plus micro kit (Qiagen, 74034). Quality and integrity of total RNA was controlled on a Bioanalyzer 2100 (Agilent Technologies). Library preparation for bulk-sequencing of poly(A)-RNA was done as described previously (Parekh et al., 2016). Briefly, barcoded cDNA of each sample was generated with a Maxima RT polymerase (ThermoFisher Scientific, EP0742) using oligo-dT primer containing barcodes, unique molecular identifiers (UMIs) and an adaptor. Ends of the cDNAs were extended by a template switch oligo (TSO) and full-length cDNA was amplified with primers binding to the TSO-site and the adaptor. NEB UltraII FS kit was used to fragment cDNA. After end repair and A-tailing, a TruSeq adapter was ligated and 3'-end-fragments were finally amplified using primers with Illumina P5 and P7 overhangs. In comparison to previous descriptions (Parekh et al., 2016), the P5 and P7 sites were exchanged to allow sequencing of the cDNA in read 1 and barcodes and UMIs in read 2 to achieve a better cluster recognition. The library was sequenced on a NextSeq 500 (Illumina) with 59 cycles for the cDNA in read 1 and 16 cycles for the barcodes and UMIs in read 2. Data were processed using the published Drop-seq pipeline (v1.0) to generate sample- and gene-wise UMI tables (Macosko et al., 2015). Reference genome (GRCh38) was used for alignment. Transcript and gene definitions were used according to the Genecode Annotation Version 35.

  Study EGAS50000000156

• MediMer: A versatile do-it-yourself peptide-receptive MHC class I multimer platform for tumor neoantigen-specific T cell detection

  Peptide-loaded MHC class I (pMHC-I) multimers have revolutionized our capabilities to monitor disease-associated T cell responses with high sensitivity and specificity. To improve the discovery of T cell receptors (TCR) targeting neoantigens of individual tumor patients with recombinant MHC molecules, we developed a peptide-loadable MHC class I platform termed MediMer. MediMers are based on soluble disulfide-stabilized β2-microglobulin/heavy chain ectodomain single-chain dimers (dsSCD) that can be easily produced in large quantities in eukaryotic cells and tailored to individual patients’ HLA allotypes with only little hands-on time. Upon transient expression in CHO-S cells together with ER-targeted BirA biotin ligase, biotinylated dsSCD are purified from the cell supernatant and are ready to use. We show that CHO-produced dsSCD are free of endogenous peptide ligands. Empty dsSCD from more than 30 different HLA-A,B,C allotypes, that were produced and validated so far, can be loaded with synthetic peptides matching the known binding criteria of the respective allotypes, and stored at low temperature without loss of binding activity. We demonstrate the usability of peptide-loaded dsSCD multimers for the detection of human antigen-specific T cells with comparable sensitivities as multimers generated with peptide-tethered β2m-HLA heavy chain single-chain trimers (SCT) and wild-type peptide-MHC-I complexes prior formed in small-scale refolding reactions. Using allotype-specific, fluorophore-labeled competitor peptides, we present a novel dsSCD-based peptide binding assay capable of interrogating large libraries of in silico predicted neoepitope peptides by flow cytometry in a high-throughput and rapid format. We discovered rare T cell populations with specificity for tumor neoepitopes and epitopes from shared tumor-associated antigens in peripheral blood of a melanoma patient including a so far unreported HLA-C*08:02-restricted NY-ESO-1-specific CD8+ T cell population. Two representative TCR of this T cell population, which could be of potential value for a broader spectrum of patients, were identified by dsSCD-guided single-cell sequencing and were validated by cognate pMHC-I multimer staining and functional responses to autologous peptide-pulsed antigen presenting cells. By deploying the novel CHO producer cell-based dsSCD MHC-I MediMer platform, we hope to significantly improve success rates for the discovery of personalized neoepitope- specific TCR in the future by being able to also cover rare HLA allotypes.

  Study EGAS50000000065

• All you need to know about our new DAC Portal v2

  An updated version to boost the management of DACs (Data Access Committees), policies, and data access requests. In September of 2023, we launched new services for all EGA users, including a DAC Portal. Since then, our team has been working to refine it and improve the user experience. The version 2 of our DAC Portal is now ready, with a range of new features aimed at streamlining your workflow and boosting productivity. If you are interested in learning more about the latest version of our DAC Portal, then you are in the right place! In this post, you will find information about the implemented features and more. The full documentation is available on our website. Moreover, we highly recommend exploring the DAC Portal through the Take the Tour. This said, let's get cracking! Customise your data access requests table This table offers the possibility to tailor information about the data access requests by making visible or hiding different columns regarding the Dataset Persistent Identifier, the Organisation or the requestor's name s, among others. You can choose which columns to see, depending on your needs and interests. Additionally, we have implemented two new features columns with two new functionalities, DAC Comment and Expiration Date. You can find more details in the following sections! In this table, all data access requests are displayed by default. With the latest version of the Portal, you can now apply filters, combine them, and save them for future use. Leave internal comments on pending requests and visible to all the DAC Members We understand that many Data Access Committees (DACs) cannot instantly approve or deny data access requests. Often, there's an intermediary phase where discussions and the signing of legal documents are necessary. To facilitate this process, we've introduced the DAC comment feature! When you navigate to the DAC column, you'll find two tabs: "User" and "DAC". Under "User", you'll see the message provided by the requester at the time of requesting access to the data. As a DAC member/admin, you can then add an internal comment in the "DAC" tab. These comments are visible only to DAC members/admins, ensuring confidentiality. Requesters won't have access to DAC comments. By adding a DAC comment to a pending request, we'll know that you've begun reviewing it. This allows you to filter pending requests without DAC comments to focus on genuinely new requests, and it prevents you from receiving email notifications about ongoing data requests. Set expiration dates for permissions and revoke them automatically The latest update to the DAC Portal introduces a new feature allowing you to set expiration dates for new access requests. For instance, if you anticipate approving a data access request but only wish to grant access for a limited time, you can now specify an expiration date for the granted permissions. Once this date is reached, the permissions will automatically be revoked. By default, permissions are granted indefinitely without an expiration date. However, if you ever need to revoke permissions, you can easily do so from the History page. Audit your metadata objects Managing and auditing your metadata objects will be easier from now on as the new DAC Portal enhances and facilitates your role as a Data Controller. This is possible thanks to the tables available in each of the primary tabs visible on the homepage: DACs, Policies and Datasets. Starting with DACs, this section allows to check the policies linked to a specific DAC. In Policies, you can obtain a list of the datasets connected to a particular policy. In the Dataset tab, it is possible to group datasets by either DAC or Policy, as well as the release status. Now you have the option of deprecating metadata objects that you don’t want to use anymore, or that were created by mistake. You will find the deprecate button on the EDIT page of the DAC, as well as for registered policies. You can find more information here. What's neat about this feature is that once a DAC or policy is deprecated, it becomes hidden within the DAC Portal. This ensures that you only see information that is consistently up to date! Define your preferences regarding email notifications For pending requests, you will be able to decide how often you want to get the notifications: daily, weekly or fortnightly. Remember, you will not receive email notifications for pending requests with a DAC comment! Furthermore, you have the flexibility to decide which notifications you wish to receive, such as invitations to join a new DAC or notifications regarding DAC acceptance by the Helpdesk, among other options. Check the User preference to see the complete list of available notification settings! What’s next: incoming DAC Portal API In addition to the new DAC Portal, we are excited to announce the release of the DAC API! This enables users to programmatically manage permissions. If you are interested in learning more about the technical specifications, you can consult the DAC API specification. This post was written by Aina Jené and Ana T. Alonso.

  Blog new-dac-portal-v2

• RNA_seq_analysis_of_transcriptome_variation_with_human_ESC_subclones

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progenyThis 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 EGAS00001001655

• The subclonal architecture of metastatic breast cancer: Results from a prospective community-based rapid autopsy program 'CASCADE'

  Background: Understanding the cancer genome is seen as a key step in improving outcomes for cancer patients. Genomic assays are emerging as a possible avenue to personalised medicine in breast cancer. The majority of work in this area has targeted primary tumours however, and very few studies have performed comprehensive profiling of advanced disease. Evolution of the cancer genome during the natural history of breast cancer is largely unknown, as is the profile of disease at death. We sought to study in detail these aspects of advanced breast cancers that have resulted in lethal disease. Methods and Findings: Three patients with ER-positive, HER2-negative breast cancer and one patient with triple negative breast cancer underwent rapid autopsy as part of an institutional prospective community-based rapid autopsy program. Cases represented a range of management problems in breast cancer, including late relapse after early stage disease; de novo metastatic disease; discordant disease response and disease refractory to treatment. Between 5 and 12 metastatic sites were collected at autopsy together with available primary tumours and longitudinal metastatic biopsies taken during life. Samples underwent paired tumour-normal whole exome sequencing and single nucleotide polymorphism arrays. Subclonal architectures were inferred by jointly analysing all samples from each patient. Mutations were validated using high depth amplicon sequencing.Between cases, there were significant differences in mutational burden, driver mutations, mutational processes and copy number variation. Within each case, we found dramatic heterogeneity in subclonal structure from primary to metastatic disease and between metastatic sites, such that no single lesion captured the breadth of disease. Metastatic cross seeding was found in each case and treatment drove subclonal diversification. Subclones displayed parallel evolution of treatment resistance in some cases, and apparent augmentation of key oncogenic drivers as an alternative resistance mechanism. We also observed the key role of mutational processes in subclonal evolution.Limitations of this study include the potential for bias introduced by joint analysis of formalin fixed archival specimens with fresh specimens, and the difficulties in resolving subclones with whole exome sequencing. Other alterations that could define subclones such as structural variants or epigenetic modifications were not assessed. Conclusions: This study highlights the variety of mechanisms that shape the genome of metastatic breast cancer, and the value of studying advanced disease in detail. Treatment drives significant genomic heterogeneity in breast cancers which has implications for disease monitoring and treatment selection in the personalised medicine paradigm.

  Study EGAS00001002153

• DNA methylomes of monozygotic twins clinically discordant for multiple sclerosis

  Multiple sclerosis (MS) is an inflammatory demyelinating disease of the central nervous system. Although genetic susceptibility is important, a modest concordance rate for MS in monozygotic (MZ) twins suggests that interaction with other risk factors is required to develop clinical symptoms. In this study, we examined whether DNA methylation differences contribute to the discordant clinical manifestation of MS in MZ twins, and studied the impact of MS treatments on the DNA methylome. Genome-wide DNA methylation profiles of peripheral blood mononuclear cells (PBMCs) of 45 MZ twins clinically discordant for MS were generated using the Illumina HumanMethylationEPIC array. Repetitive element methylation and selected differentially methylated positions (DMPs) were validated using targeted deep bisulfite sequencing (TDBS). In addition, we performed whole genome bisulfite sequencing (WGBS) to profile CD4+ memory T-cells of a subset of four MS discordant MZ twins. Our results show that interferon treatment causes robust DNA methylation changes and several epigenetic biomarkers for interferon treatment response were identified. However, overall the PBMC-based methylomes of the MS discordant MZ twins were highly similar, since large systematic methylation differences (&gt;5%) were absent in the data. This suggests that previously reported large methylation changes are probably caused by genetic rather than epigenetic differences. In addition, our data does not support the hypothesis that the observed maternal parent-of-origin effect in MS is due to genomic imprinting errors, and no evidence was found that copy number variations explain the discordant phenotype in these MZ twins. Although not genome-wide significant, a couple of DMPs associated with the MS phenotype were identified and successfully technically replicated using TDBS, including a differentially methylated region (DMR) in the promoter of the transmembrane protein encoding gene TMEM232. Another MS-associated DMP, located in an enhancer in the gene body of the transcription factor ZBTB16, was also associated with medium-term glucocorticoid treatment history. WGBS analysis confirmed this DMP as a promising epigenetic biomarker for glucocorticoid treatment. In conclusion, this epigenome-wide association study (EWAS) in genetically identical twins identified a DMR in the TMEM232 promoter as a candidate loci associated with the clinical manifestation of MS. In addition, epigenetic biomarkers for MS treatments were identified, revealing that not only short-term, but also medium-term treatment effects are detectable in immune cells, which should be taken into account in future EWAS designs.

  Study EGAS00001003147

• PSCP_bisulphite_analysis_in_hESCs

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progeny. 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 EGAS00001001625

• PSCP_mutation_analysis_in_hESCs

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progeny. This data is part of a pre-publication release. For information on the proper use of pre-publication data shred by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001001561

• Pediatric Papillary Thyroid Carcinoma RNA-Seq

  Papillary thyroid carcinoma (PTC) is the most common malignancy amongst adolescent and young adult women. Clinical and genomic features of PTC in children differ from those in adults. Historically, a substantial proportion (~50%) of pediatric PTC (PPTC) lack any of the common driver mutations typical of adult PTC. Thus, molecular diagnostics developed for adults may not be valid in children. We applied novel bioinformatic pipelines to RNASeq analysis of pediatric thyroid tumors to identify known and novel fusion oncogenes as drivers of tumorigenesis and to define gene expression patterns among tumors. Analysis of 38 samples resulted in identification of 16 previously reported oncogenic fusions and three novel candidate fusions. In conducting these analyses, we reduced the proportion of cancers with unknown drivers from approximately 50% to 6%. Gene expression analysis identified three distinct categories of tumors. These studies represent the first unbiased genomic approach to childhood PTC. Novel bioinformatic pipeline. GEC Identified 3 distinct classes. DICER1, PTEN segregated with benign tumors. &gt;90% of PPTC carry identifiable ongognic driver variants, however, the landscape differs from that in adults and existing tools for tumor genotyping in adults may not be applicable to pediatric tumors. BACKGROUND: Papillary thyroid carcinoma demonstrates clinically distinct behaviour in pediatric patients when compared to adults, with increased regional and distal metastases, more frequent recurrence, yet overall highly favourable prognosis. Recent tumor genotyping has demonstrated that in pediatric tumors, ***-*** do not carry the most common genomic variants found in adult PTC, leaving a substantial proportion of driver-unknown or “dark matter” tumors. METHODS: We established a cohort of 52 pediatric PTC and performed quantitative-PCR genotyping followed in a subset by whole exome and whole transcriptome (RNASeq) sequencing using a novel, robust bioinformatic pipeline to identify oncogenic fusion transcripts. Gene expression analysis was used to derive distinct clusters. RESULTS: In 52 pediatric PTC, genotyping identified BRAFV600E(21%), NRASQ61R(2%), RET/PTC1 (21%), RET/PTC3 (17%) or PAX8-PPAR in 33 (63%) of tumors. Among the remaining 19 tumors, oncogenic variants were identified in /19 (%). The variants included SNVs in DICER1, ETV6-NTRK3, PTEN and 9 additional oncogenic fusions including two previously-undescribed fusion transcripts, *** and ***. CONCLUSIONS: The genomic landscape of PTC differs between children and adults. In children, there is a substantial fraction of tumors that are unexplained by variants typical of adult PTC. This discrepancy is explained largely by uncommon oncogenic fusions which can be identified by whole-transcriptome analysis.

  Study EGAS00001005182

• Single cell genomic variation induced by mutational processes in cancer

  How cell-to-cell copy number alterations that underpin genomic instability in human cancers drive genomic and phenotypic variation, and consequently cancer evolution, remains understudied. Here, applying scaled single cell whole genome sequencing to wildtype, TP53-deficient and TP53/BRCA1/2-deficient mammary epithelial cells (n=13,818 genomes) and primary triple negative breast (TNBC) and high grade serous ovarian cancers (HGSC) (n=22,057 genomes), we identified three distinct ‘foreground’ mutational patterns defined by cell-to-cell copy number and structural variation properties occurring in the ‘background’ of cancer-associated genomic instability. Cell and clone-specific high level amplifications (HLAMPs), parallel allele-specific copy number alterations and copy number segment length variation (serrate structural variations) exhibited measurable phenotypic and evolutionary consequences. Clone-specific HLAMPs in known oncogenes were highly prevalent in fold back inversion (FBI)-bearing tumors and associated with increased clone-to-clone phenotypic variation in gene expression inferred from matched single cell RNASeq. Parallel allele-specific alterations were linked with phylogenetically determined evolutionary diversity and clone-specific mono-allelic expression. Serrate patterns present across all clones were increased in FBI tumors, most prominently on a tetraploid background and were highly correlated with increased genomic diversity of cellular populations. Taken together, we show that cell-to-cell structural-copy number variation significantly impacts the phenotypes and evolutionary diversity of TNBC and HGSC, revealing previously hidden genomic states of cancer cells.

  Study EGAS00001006343

• TNBC ctDNA Targeted Panel

  Circulating tumour DNA (ctDNA) detection in liquid biopsy is an emerging alternative to tissue biopsy, but its utility in treatment response monitoring and prognosis in triple negative breast cancer (TNBC) remains to be well understood. In this study, we determined the presence of early stage TNBC ctDNA detectable actionable mutations with a clinically validated hotspot treatment indication panel. Sequencing of plasma DNA from 130 TNBC patients collected within 7 months of primary treatment completion revealed that 7.7% had detectable residual disease with a hotspot panel. Among neoadjuvant treated patients, we observed a trend where patients with incomplete pathologic response and positive ctDNA within 7 months of treatment completion were at risk of reduced progression free survival. We propose that a high risk subset of early TNBC patients treated in NAT protocols may be identifiable by combining tissue response and sensitive ctDNA detection.

  Study EGAS00001006937

• TRACERx NSCLC - Whole exome multiregion sequencing data from the 421 cohort

  TRACERx NSCLC - Whole exome multiregion sequencing data from the 421 cohort

  Dataset EGAD00001009825

• PanCuRx Translational Research Initiative

  The goal of the PanCuRx TRI is to seek solutions to the high fatality rate of pancreatic adenocarcinoma (PDAC), the most common type of pancreatic cancer, by generating new knowledge about the genetics and biology of the disease, mechanisms of how tumours grow and tailored treatment options. PDAC is the fourth leading cause of cancer death in Canada and is expected to become the second-leading cause of cancer death within the next decade. The current five-year survival rate of eight per cent is the lowest of all epithelial cancers. There are major clinical challenges to treating the disease, including the fact that PDAC spreads to other parts of the body early so surgical removal of the tumour benefits few patients, that PDAC has proven relatively resistant to systemic therapies such as chemotherapy and that there has been limited benefit from the use of newer molecular targeted agents. The PanCuRx team has created the world’s first resource of tumour-purified PDAC whole genomes; more than 500 PDAC samples, with combined clinical annotation, have been successfully sequenced and analyzed at the Ontario Institute for Cancer research in Toronto. The resulting clinical genomic resource created from this data has led to multiple new findings that have established a framework to better understand PDAC biology and lays the foundation for new genomic approaches to better understanding the disease. The team also launched a clinical trial called COMPASS (Comprehensive Molecular Characterization of Advanced Ductal Pancreas Adenocarcinoma for Better Treatment Selection: A Prospective Study), which is recruiting patients with advanced PDAC (both locally advanced and metastatic) at The Princess Margaret Cancer Centre (PM) in Toronto. The study uses the molecular results from individual patient tumours to guide selection of better second line treatment on a case-by-case basis. The primary objective of the initial study was to assess the feasibility of collecting tumour materials from patients with advanced PDAC and delivering results within 56 days of biopsy. Having achieved this objective, the trial has expanded to other cancer centres across Canada. Pancreatic cancer is one of the most deadly forms of cancer and outcomes have not improved in four decades. The PanCuRx initiative is providing a better understanding of the disease and actively applying new knowledge to cancer patients through clinical trials. This knowledge will help to improve quality of life and prolong survival for pancreatic cancer patients worldwide.

  Study EGAS00001002543

• Whole transcriptome RNA sequencing as comprehensive diagnostic tool for acute myeloid leukemia.

  BackgroundAcute myeloid leukemia (AML) is characterized by uncontrolled proliferation of malignant hematopoietic cells in the bone marrow that are arrested in differentiation. In AML pathogenesis, hematopoietic progenitor cells acquire multiple genetic aberrations often occurring in the same set of genes that ultimately lead to malignant transformation. In the WHO classification 2016, six AML classes with different prognosis are identified by chromosomal translocations measured by standard cytogenetics. All balanced translocations produce fusion genes, except for t(3;3)/inv(3;3), which lead to overexpression of MECOM/EVI1 associated with poor prognosis. For accurate risk assessment of cytogenetically normal AML, four genes need to be screened to distinguish AML with mutated NPM1 in the absence of FLT3 mutations and AML with bi-allelic CEBPA mutations, which have a favorable prognosis, and AML with RUNX1 mutations which have a poor prognosis. Recently, a full genomic classification system has been proposed by Papaemmanuil et al. (NEJM 2016). In this system, the same six classes with chromosomal translocations are identified by standard cytogenetics as well as five additional classes characterized by genetic mutations in 14 different genes. AimThe aim of the study was to investigate whether whole transcriptome RNA-sequencing (RNAseq) can be used as single technology for classification of AML. MethodA panel of 100 AML were analyzed and a bio-informatics pipeline called HAMLET (Human AML Expedited Transcriptomics) was developed in which 4 modules are integrated to detect (1) fusion genes, (2) small variants in 13 recurrently mutated genes, (3) internal tandem duplications in FLT3 (FLT3-ITD) and partial tandem duplications in KMT2A (KMT2A-PTD) and (4) overexpression of MECOM/EVI1. All mutations that were called by HAMLET were validated by diagnostic data as available for all 100 AML or targeted PCR followed by Sanger or next generation sequencing on all positive AML and at least an equal number of negative cases. Results &amp; discussionThe data showed that HAMLET accurately identified all genetic aberrations with high sensitivity and specificity. In addition, in 7 AML, fusion transcripts were detected that are not measured by standard cytogenetics including three cases that lack any class-defining lesion according to Papaemmanuil et al. Moreover, overexpression of MECOM/EVI1 was detected in two AML with inv(3) as well as in five cases without inv(3)/t(3;3), and a gene signature was measured to distinguish AML with CEBPA mutations with favorable prognosis. In conclusion, HAMLET provides a comprehensive and flexible pipeline for RNAseq analysis to retrieve all relevant information for current classification of AML as well as additional information that may improve classification in the future.

  Study EGAS00001003096

• Molecularly matched targeted therapies plus radiotherapy in patients with newly diagnosed glioblastoma without MGMT promoter hypermethylation (N2M2/NOA-20 phase I/IIa umbrella trial)

  Background: Patients with glioblastoma without MGMT promoter hypermethylation are unlikely to benefit from temozolomide (TMZ). Trials aiming at replacing TMZ with targeted agents in not molecularly selected patient populations have failed. Methods: This phase I/IIa umbrella trial aimed at showing safety, feasibility, and preliminary efficacy of targeted compounds in addition to standard radiotherapy initiated within 42 days postoperatively. Molecular diagnostics and bioinformatic evaluation are performed within 28 days after surgery. Stratification for treatment takes place in five subtrials, including alectinib, idasanutlin, palbociclib, vismogedib and temsirolimus, according to the best matching molecular alteration. Patients without matching alterations are randomized between subtrials without strong biomarkers using atezolizumab and asunercept and TMZ as standard of care. Primary objective of the phase I parts of the trial was dose finding or dose validation. In the phase IIa trials, centrally determined progression-free survival at six months (PFS-6) is used as endpoint for efficacy with interim analyses for futility (H0: p=0.231). Results: From May 2018 through July 2022, 301 patients were enrolled and 228 treated in 13 German NOA sites. The alectinib and vismodegib subtrials were closed since no molecularly matching patients were accrued; the idasanutlin subtrial was closed prior to the optimal dose at nine patients at discretion of the company providing the drug. The TMZ subtrial showed a PFS-6 of 18.52% (10/54 patients) (p=0.831) and a median overall survival (OS) of 12.1 months. Asunercept: PFS-6 of 15.4% (4/26) (p=0.8825) and OS of 12.8 months. Atezolizumab: PFS-6 of 21.4% (9/42) (p=0.660) and OS of 11.7 months. Palbociclib with patients demonstrating CDK4 amplification or CDKN2A/B codeletion: PFS-6 of 24.4% (10/41) (p=0.4823) and OS of 12.6 months. Temsirolimus with patients demonstrating mTOR activation: PFS-6 of 39.1% (18/46) (p=0.0109) and OS of 15.4 months. The regimen-limiting toxicity (RLT)-rate is 34.8%, which is insignificantly above the predefined unacceptable rate for RLTs of 30%. Most RLTs had severity grade 3, one RLT had severity grade 4. No RLTs resulted in death. Conclusions: N2M2 allows for elaborate molecular testing being integrated into the treatment decision and efficient determination of treatment activity for patients with newly diagnosed glioblastoma. There is clinical activity of temsirolimus in patients with tumors harboring an activated mTOR pathway although this is not positively prognostic without mTOR inhibition; there is no clinical activity for asunercept and atezolizumab in not molecularly selected patients and also palbociclib in molecularly selected patients.

  Study EGAS00001008033

• sn-RNAseq profiling of the impact of a cytokine storm model in human cardiac organoids

  We developed conditions to model conditions of diastolic dysfunction caused by inflammatory factors in human cardiac organoids. Our organoid model is based on conditions promoting maturation (Mills et al., PNAS, 2017 and Voges et al., In Revision) and incorporating human pluripotent stem cell derived cardiomyocytes, epicardial cells, fibroblasts, endothelial cells and pericytes. Following 15 days of differentiation (Voges et al., Development, 2017) we formed human cardiac organoids incorporating both our multicellular differentiation (Mills et al., PNAS, 2017) and additional endothelial cells (Voges et al., In Revision). After another 12 days of maturation in 5 days maturation medium and 5 days weaning medium (Voges et al., In Revision) human cardiac organoids were either treated with weaning medium or “cardiac cytokine storm” comprised of 100 ng/ml IFNg, 10 ng/ml IL-1 β and 10 µg/ml poly(I:C). After 48 hours human cardiac organoids were harvested and snap frozen in -80C before nuclei isolation and single nuclei RNA sequencing.Pooled hCO (~40) were homogenized in 4 mL lysis buffer (300 mM sucrose, 10 mM Tris-HCl (pH = 8), 5 mM CaCl2, 5 mM magnesium acetate, 2 mM EDTA, 0.5 mM EGTA, 1 mM DTT)(all Sigma-Aldrich) with 30 strokes of a dounce tissue grinder (Wheaton). Large pieces of hCO were allowed to settle and homogenate was passed through pre-wetted 40 µm cell strainers (Becton Dickinson). Remaining hCO material in the douncer was resuspended in 4 mL and the douncing and filtering steps were repeated twice. All steps of the homogenization were performed on ice. The filtered homogenate was centrifuged at 1500 x g for 5 min at 4oC. Nuclei pellets were then re-suspended in PBS. A fraction of resuspended nuclei were then stained with hoechst nuclear stain (1:500 dilution) and counted on a haemocytometer under fluorescent microscope.The nuclei were then re-centrifuged (1500 x g for 5 min at 4°C) and resuspended at a density to load ~5,000 nuclei per sample. Cell were loaded into the Chromium Controller (10X Genomics) for gel bead emulsion (GEM) formation. Library preparation was conducted according to the manufacturer’s recommended protocol using the Chromium Next GEM Single Cell 3’ GEM, Library &amp; Gel Bead Kit v3.1. Libraries were sequenced on the NextSeq 500/550 v2 (Illumina) with 150 bp reads.In the snRNA-seq we note the lower than expected percentage of non-myocytes and the loss of endothelial cells (Mills et al., PNAS, 2017 and Voges et al., In Revision), indicating the protocol requires further optimization for hCO samples.

  Study EGAS00001005174

• Beyond BRCA deficiency: Clinical and molecular predictors of survival in patients with BRCA-deficient tubo-ovarian high-grade serous carcinoma

  Homologous recombination DNA repair deficient (HRD) tumors account for about 50% of tubo-ovarian high-grade serous carcinomas (HGSC) and are associated with a high response to platinum-based chemotherapy. However, a subset of patients with HRD tumors have unexpected resistance to chemotherapy and experience short survival. Our study compared BRCA-deficient and proficient HGSC among patients with short (&lt;3 years, the lowest quartile of overall survival (OS)) versus longer-term (&gt;3 year OS) to identify prognostic markers and potential therapeutic opportunities. Primary tumors from a cohort of patients with advanced-stage HGSC (n=154) enriched for those with HRD tumors and short survival (n=42) were analyzed using whole-genome sequencing, bulk RNA-sequencing, and DNA methylation profiling. A subset of HRD tumors (n=143) were characterized for immune cell densities by multiplex immunofluorescence. Clinicopathological features were assessed in 1,389 patients with HGSC, including 282 with pathogenic germline BRCA variants (gBRCApv). Findings were validated in two international HGSC cohorts: the Ovarian Tumor Tissue Analysis consortium (n=5,875) and Multisciplinary Ovarian Cancer Outcome Group (n=895). In gBRCApv-carriers, residual disease following surgical cytoreduction was not associated with shorter OS unlike in non-carriers where it was strongly associated. While most patients with gBRCApv had a more favorable OS compared to non-carriers, a survival advantage was not apparent in those with variants located outside of BRCA1 exon 10/23 and BRCA2 exon 11/27. Patient survival correlated with HRD score, such that those with a score ≥63 had significantly longer OS compared to those with HRD scores of 42-63 and &lt;42. While higher HRD scores were associated with better survival outcomes, most gBRCApv carriers with a short survival had tumor HRD scores that exceeded the threshold typically considered predictive of a good response to carboplatin or PARP inhibitors. BRCA2-deficient HGSC with loss of the tumor suppressor NF1 was associated with exceptionally favorable outcome (median OS 17 years), whereas amplification of both PIK3CA and RAD21, was associated with poor outcome in BRCA2-deficient HGSC (median OS 2.9 years). The density of tumor-infiltrating CD8+ PD-1+ T cells was prognostic for OS in gBRCApv-carriers. In summary, rather than simply the presence or absence of a germline BRCA pathological variant influencing patient survival, the co-occurrence of specific somatic gene alterations, the extent of immune cell infiltration, the level of homologous recombination deficiency, and BRCA-variant location may modify clinical outcomes in patients with BRCA-deficient HGSC. Analysis of exceptional outcome patients, in this case gBRCApv-carriers with an atypically short survival, may identify determinants of outcome applicable to the wider HGSC patient population.

  Study EGAS00001008059

• Covid19 Gene expression

  The dataset contains the gene expression profile of each individual along with gene fusion events.

  Dataset EGAD00001010185

• DDD DATAFREEZE 2013-12-18: 1133 trios - VCF files (Ref: DDD Nature 2015)

  DDD DATAFREEZE 2013-12-18: 1133 trios - VCF files (Ref: DDD Nature 2015)

  Dataset EGAD00001001355

• DDD DATAFREEZE 2013-12-18: 1133 trios - exome sequence BAM files (Ref: DDD Nature 2015)

  DDD DATAFREEZE 2013-12-18: 1133 trios - exome sequence BAM files (Ref: DDD Nature 2015)

  Dataset EGAD00001001114

• Robust detection of translocations in lymphoma FFPE samples using Targeted Locus Capture-based sequencing

  Chromosomal translocations with immunoglobin (IG) loci are the classic drivers in a large subset of B-cell lymphomas. Detection of these translocations is important for confirmation of diagnosis and for prognosis and therapy decisions. Currently, molecular diagnosis of translocations in lymphomas is not addressed well by Next Generation Sequencing (NGS). The standard method for detection of translocations is Fluorescence In Situ Hybridization (FISH), which is labor-intensive, and can be difficult to interpret. There is a need for a robust technology that can be standardized. Targeted Locus Capture (TLC) selectively enriches and sequences entire genes based on the crosslinking of physically proximal sequences, and thereby enables complete sequencing of genes of interest, including detection of large structural variants. Because the technology is based on the crosslinking and fragmenting of DNA, it has particular advantages in the analysis of Formalin-Fixed, Paraffin-Embedded (FFPE) samples, in which DNA is inherently crosslinked and fragmented. In order to validate the FFPE-TLC technology as a novel approach for translocation detection in lymphoma samples, we have developed a panel assay containing genes with frequent translocations (MYC, BCL2, BCL6, IG loci). With this assay we have analyzed &gt;140 lymphoma and control FFPE samples of variable input amounts and qualities that had previously been analyzed with FISH, and a subset also with standard targeted NGS. Good concordance with FISH results was observed for both translocation positive and negative samples. In 10 cases for which FFPE-TLC detected a relevant fusion and FISH had been called negative, discordance could be explained by higher sensitivity of FFPE-TLC or by inconclusive FISH results. In a specific case, FFPE-TLC detected a small-distance rearrangement on chromosome 3 that caused a BCL6 fusion but led to insufficient and therefore undetectable break-apart with FISH. Secondly, the FFPE-TLC approach was tested on a set of 19 B-cell lymphoma FFPE samples that had previously been analyzed using standard targeted NGS and FISH and was enriched for discordant results between these methods. FFPE-TLC-based NGS enables more robust translocation calling as the detection relies on broad sequencing coverage across the translocation partner rather than on breakpoint sequences only. In 3 cases, FFPE-TLC could proof false negative calls in standard targeted NGS due to breakpoints located in regions difficult to capture or to sequence. In 1 case, standard targeted NGS had made a false positive call on a breakpoint sequence that was shown to be caused by a small insertion rather than a genuine translocation. This study shows that FFPE-TLC promises to be a robust alternative for FISH analysis and standard targeted NGS procedures in lymphoma diagnostics and possibly in other cancers with frequent structural variants. The FFPE-TLC approach enables a single, DNA-based NGS test detecting both small mutations and translocations.

  Study EGAS00001004760

• Capturing sex-specific and infertility-linked effects of assisted reproductive technologies on the cord blood DNA methylome

  Background: Children conceived through assisted reproduction are at an increased risk for growth and genomic imprinting disorders, often linked to DNA methylation defects. It has been suggested that assisted reproductive technologies (ART) and the underlying parental infertility can induce epigenetic instability, specifically interfering with DNA methylation reprogramming events during germ cell and preimplantation development. To date, human studies exploring the association between ART and DNA methylation defects report inconsistent or inconclusive results, likely due to population heterogeneity and the use of technologies with limited coverage of the epigenome. In our study, we explore the epigenetic risk of ART by comprehensively profiling the DNA methylome of 73 human cord blood samples from singleton pregnancies (n=36 control group, n=37 ART/infertility group) from a human prospective longitudinal birth cohort, the Design, Develop, Discover (3D) Study using a high resolution sequencing based custom capture panel that interrogates over 2.4 million autosomal CpGs in the genome. Results: the main findings We identified striking sex-specific effects of ART on cord blood DNA methylation patterning both genome-wide and for imprinted genes. When compared to less invasive procedures such as in utero insemination, more invasive ARTs (IVF, ICSI, embryo culture) resulted in more marked and distinct effects on the cord blood DNA methylome, including targeting of loci involved in development. Background: Children conceived through assisted reproduction are at an increased risk for growth and genomic imprinting disorders, often linked to DNA methylation defects. It has been suggested that assisted reproductive technologies (ART) and the underlying parental infertility can induce epigenetic instability, specifically interfering with DNA methylation reprogramming events during germ cell and preimplantation development. To date, human studies exploring the association between ART and DNA methylation defects report inconsistent or inconclusive results, likely due to population heterogeneity and the use of technologies with limited coverage of the epigenome. In our study, we explore the epigenetic risk of ART by comprehensively profiling the DNA methylome of 73 human cord blood samples from singleton pregnancies (n=36 control group, n=37 ART/infertility group) from a human prospective longitudinal birth cohort, the Design, Develop, Discover (3D) Study, using a high resolution sequencing based custom capture panel that interrogates over 2.4 million autosomal CpGs in the genome. Results: We identified striking sex-specific effects of ART on cord blood DNA methylation patterning both genome-wide and for imprinted genes. When compared to less invasive procedures such as in utero insemination, more invasive ARTs (in vitro fertilization (IVF), intracytoplasmic sperm injection (ICSI), embryo culture) resulted in more marked and distinct effects on the cord blood DNA methylome, including targeting of loci involved in development. Conclusions: Our study highlights the ability of a sensitive, targeted, sequencing based approach to uncover DNA methylation perturbations in cord blood associated with infertility and ART and influenced by offspring sex and ART technique invasiveness.

  Study EGAS00001006643

• Solve-RD_ITHACA_cohort-1_DF1+2_V1

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

  Dataset EGAD00001009770

• Koean HCC exome sequencing

  231 HCC exome sequencing with Sureselect 50Mb

  Dataset EGAD00001005361

• covid19 RNASeq raw fastq read

  446 samples of covid19 patients. Raw Reads in fastq format.

  Dataset EGAD00001010079

• Irish Covid19 RNASeq Alignment Data files

  The Dataset contains 446 Covid19 patient's RNASeq Alignment files in BAM format for both genomic and transcriptomic alignment.

  Dataset EGAD00001010099

• DDD DATAFREEZE 2013-12-18: 1133 trios - README, family trios, phenotypes, validated DNMs (Ref: DDD Nature 2015)

  DDD DATAFREEZE 2013-12-18: 1133 trios - README, family trios, phenotypes, validated DNMs (Ref: DDD Nature 2015)

  Dataset EGAD00001001413

• Institut Curie Neuroblastoma Whole Genome Sequencing Diagnosis Relapse

  Neuroblastoma, a clinically heterogeneous pediatric cancer, is characterized by distinct genomic profiles but few recurrent mutations. As neuroblastoma is expected to have high degree of genetic heterogeneity, study of neuroblastoma's clonal evolution with deep coverage whole-genome sequencing of diagnosis and relapse samples will lead to a better understanding of the molecular events associated with relapse. Samples were included in this study if sufficient DNA from constitutional, diagnosis and relapse tumors was available for WGS. Whole genome sequencing was performed on trios (constitutional, diagnose and relapse DNA) from eight patients using Illumina Hi-seq2500 leading to paired-ends (PE) 90x90 for 6 of them and 100x100 for two. Expected coverage for sample NB0175 100x100bp was 30X for tumor and constitutional samples. For the seven other patients expected coverage was 80X for tumor samples with PE 100x100, 100X in the other tumor samples and 50X for all constitutional samples (see table 1). Following alignment with BWA (Li et al., Oxford J, 2009 Jul) allowing up to 4% of mismatches, bam files were cleaned up according to the Genome Analysis Toolkit (GATK) recommendations (Van der Auwera et al., Current Protocols in Bioinformatics, 2013, picard-1.45, GenomeAnalysisTK-2.2-16). Variant calling was performed in parallel using 3 variant callers: GenomeAnalysisTK-2.2-16, Samtools-0.1.18 and MuTect-1.1.4 (McKenna et al., Genome Res, 2010; Li et al., Oxford J, 2009 Aug; Cibulskis et al., Nature, 2013). Annovar-v2012-10-23 with cosmic-v64 and dbsnp-v137 were used for the annotation and RefSeq for the structural annotation. For GATK and Samtools, single nucleotide variants (SNVs) with a quality under 30, a depth of coverage under 6 or with less than 2 reads supporting the variant were filter out. MuTect with parameters following GATK and Samtools thresholds have been used to filter our irrelevant variants. .SNVs within and around exons of coding genes overlapping splice sites.. Then,variants reported in more than 1% of the population in the 1000 genomes (1000gAprl_2012) or Exome Sequencing Project (ESP6500) have been discarded in order to filter polymorphisms. Finally, synonymous variants were filtered out. MuTect focuses on somatic by filtering with constitutional sample. Mpileup comparison between constitutional and somatic DNAs allowed us to focus also on tumor specific SNVs with GATK and Samtools. Finally, every SNV called by our pipeline and also supported in any constitutional samples were filtered our in order to prevent putative constitutional DNA coverage deficiency. Then we analyzed CNVs (copy number variants) with HMMcopy-v0.1.1 (Gavin et al., Genome Res, 2012) and control-FREEC-v6.7 (Boeva et al., Bioinformatics 2011) with a respective window of 2000bp and 1000 bp, and auto-correction of normal contamination of tumor samples for Control-FREEC. Finally we explored Structural variants (SVs) including deletions, inversions, tandem duplications and translocations using DELLY-v0.5.5 with standard parameters (Rausch et al., Oxford J, 2012). In tumors, at least 10 supporting reads were required to make a call and 5 supporting reads for the sample NB0175 with a coverage of only 40X (see table 2). To predict SVs in constitutional samples for subsequent somatic filtering, only 2 supporting reads were required in order not to miss one. To identify somatic events, all the SVs in each normal sample were first flanked by 500 bp in both directions and any SVs called in a tumor sample which was in the combined flanked regions of respective normal sample was removed (see graph 1). Deletions with more than 5 genes impacted or larger than 1Mb and inversions or tandem duplications covering more than 4 genes, were removed. We focused on exonic and splicing events for deletions, inversions, and tandem duplications. For translocation, we keep all SVs that occurred in intronic, exonic, 5'UTR, upstream or splicing regions. Bioinformatics detection of variations with Deep sequencing approach Once PE reads merged and adaptors trimmed by SeqPrep with default parameters, merged reads were aligned via the BWA (Li H. and Durbin R. 2009 PMID 19451168) allowing up to 1 differences in the 22-base-long seeds and reporting only unique alignments. Only reads having a mapping quality 20 or more have been further analysed. Variant calling software was not used, since we aimed to predict variations at low frequencies, observed in less than 1% of reads. Such variants require a custom approach. Using DepthOfCoverage functions of the Genome Analysis Toolkit (GATK) v2.13.2 (McKenna A, et al., 2010 Genome Research PMID: 20644199), we focused on high quality coverage of bases A, C, G and T at the targeted variant position. Depth of coverage of each base following a mapping quality higher than 20 and a base quality higher than 10 have been taken into account in order to focus only on high quality data. Aiming to determine the background level of variability at the studied regions, 10 control samples were included in the analysis. The same approach and filtering criteria have been applied as introduced above over the entire amplicons. In order to highlight variants, for each sample the frequencies of each bases at each amplicon position were then compared to those observed in the set of controls. Statistical analyses were performed with the R statistical software (http://www.R-project.org). Fisher’s exact two-sided tests with a Bonferroni correction were performed to compare percentages of bases between the data sets, i.e. for a given base between a case and the controls. Finally, significant variations were filtered-in once (i) a significant increase in the percentage of avariant base and (ii) a significant decrease in the percentage of it's reference base following our p.values criteria was observed (p.val &lt; 0.05).

  Dataset EGAD00001001356

• A_study_of_the_genetic_basis_of_evation_by_Acute_Myeloid_Leukaemia_of_Graft_vs_Leukaemia_effects_after_allogeneic_bone_marrow_transplantation

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

  Study EGAS00001000145

• The BC Cancer Agency's Personalized Onco-Genomics Project

  The BC Cancer Agency's Personalized Onco-Genomics (POG) Project involves molecular and clinical characterization of individual cancer patients to guide their treatment. The project aims to use whole genome and transcriptome data together with patient information to match treatments to cancer driver alterations and support clinical decision-making. Incorporating genomic data into decision-making should enable more rational development of treatment plans, facilitate enrollment into appropriate clinical trials and potentially identify less toxic and more effective therapeutic options.

  Study EGAS00001001159

• Genomic Profiling Reveals Spatial Intra-tumour Heterogeneity in Follicular Lymphoma

  Follicular lymphoma (FL) is an incurable indolent non-Hodgkin’s lymphoma (NHL) that behaves heterogeneously and is synonymous with the reciprocal translocation t(14;18)(q32;q21), leading to aberrant constitutive over-expression of BCL2 [ref]. High risk groups with poor overall survival include patients that transform to a high grade lymphoma, most commonly diffuse large B cell lymphoma (DLBCL) or progress early (within 2 years) of receiving treatment. Recently several studies have shed light on the diverse genomic landscape of FL, with frequently mutated genes including those involved in chromatin remodelling, (CREBBP, KMT2D, EZH2, EP300), immune modulation (TNFRSF14), JAK-STAT (STAT6, SOCS1), NF-ĸB (CARD11, MYD88, TNFAIP3) and mTORC1 signalling (RRAGC) [ref]. To date, the majority of these studies have either focussed on inter-patient genetic heterogeneity or established the extent of clonal heterogeneity as a patient’s tumour evolves over time. Furthermore, these temporal studies have alluded to the existence of a tumour-repopulating population (referred to as the common progenitor cell (CPC)) that evades treatment and acts as a reservoir, seeding each subsequent relapse and transformation event. Although these previous studies have illuminated the longitudinal clonal dynamics that occur in FL, our understanding of the degree of spatial or intra-tumour heterogeneity (ITH) that exists within an individual patient thus far remains limited. Earlier studies have predominantly focused on differences in cytological grade and immunophenotype between spatially separated lymph nodes and bone marrow, the most common site of extra-nodal involvement in FL, with genomic profiling restricted to examining somatic hypermutation patterns in IgVH genes [ref]. In contrast next generation sequencing has been increasingly utilised in solid organ malignancies to derive comprehensive genomic profiles from spatially separated sites. The seminal study by Gerlinger et al revealed profound differences in the genetic make-up between individual primary and metastatic sites of renal cell carcinoma [ref]. Significant heterogeneity has since been demonstrated both between and within lesions in the same patient in lung, pancreas and breast cancer with presumed driver mutations distributed within the “branches” of the evolutionary phylogenetic tree, suggesting that a single biopsy is incapable of capturing the full genomic heterogeneity of an individual’s malignancy [ref]. This geographical heterogeneity not only adds to the diverse pool of tumour subclones that may contribute to drug resistance mechanisms but importantly, presents a major obstacle for precision-based approaches focussed on targeting specific lesions within a single biopsy. In FL, the exponential increase in clinical trials using novel agents such as EZH2, PI3K and BTK inhibitors reflects this shift in cancer care and along with the development of molecular prognostic scores such as the m7-FLIPI highlights the clinical need to accurately define genomic alterations with clinical relevance. As the majority of FL patients manifest disseminated tumour involvement, we sought to uncover the extent and clinical importance of spatial heterogeneity in FL by comprehensively genetically profiling 22 synchronously removed spatially separated biopsies from 9 patients. Using a combination of whole exome and targeted deep sequencing, our results inferred the complex subclonal architecture within these tumours and the mutational differences between anatomical sites, demonstrating cases with significant intrinsic genetic diversity.

  Study EGAS00001002492

• Transcriptomic analysis of liver CD8+ T cells

  Sorted single CD8+T cells expressing CD14 from human liver for SMARTSeq2. Livers processed: Kucykowicz et al STAR Prot 2022:pubmed.ncbi.nlm.nih.gov/35516846/ Published: Pallett et al Nature 2022 Tissue CD14+CD8+T-cells are reprogrammed by myeloid cells and modulated by LPS The liver is bathed in bacterial products including LPS transported from the intestinal portal vasculature, but maintains a state of tolerance exploited by persistent pathogens and tumours1,2,3,4. The cellular basis mediating this tolerance, yet allowing a switch to immunity or immunopathology, needs to be better understood for successful immunotherapy of liver diseases. Here we show that a variable proportion of CD8+T-cells compartmentalised in the human liver co-stain for CD14 and other prototypic myeloid membrane proteins; CD14+CD8+T-cells preferentially accumulate within the donor pool in liver allografts, amongst hepatic virus-specific and tumour-infiltrating responses, and in cirrhotic ascites. CD14+CD8+T-cells exhibit increased turnover, activation and constitutive immunomodulatory features with high homeostatic IL-10/IL-2 production directly ex vivo, and an enhanced capacity to rapidly produce antiviral/anti-tumour cytokines upon TCR engagement. This CD14+CD8+T-cell profile can be recapitulated by acquisition of membrane proteins, including the LPS receptor complex, from mononuclear phagocytes, resulting in augmented tumour killing by TCR-redirected T-cells in vitro. CD14+CD8+T-cells express integrins and chemokine receptors favouring interactions with the local stroma, which can promote their induction through CXCL12. LPS can also increase the frequency of CD14+CD8+T-cells in vitro and in vivo, and skew their function towards the production of chemotactic and regenerative cytokines. Thus, bacterial products in the gut-liver axis and tissue stromal factors can tune liver immunity by driving myeloid instruction of CD8+T-cells with immunomodulatory capacity. For the analysis of the single cell transcriptome two intrahepatic perfusate samples were isolated and frozen at -80 'C for storage prior to transportation to Newcastle, UK. Once thawed in pre-warmed HBSS containing 0.001% DNaseI, intrahepatic leukocytes were sorted by FACS to index classical mononuclear phagocytes (MNP) defined as Singlets, CD3-CD56-CD19-CD14hiHLA-DR+ or CD14± CD8+T-cells as follows: Singlets, CD3+CD56-CD19-CD4-CD8+V7.2- and split in to CD14± based on an isotype control. Single cells were sorted into 96-well skirted LoBind plates (Eppendorf) containing 5 µL lysis buffer (TCL buffer [Qiagen], supplemented with 1% b-mercaptoethanol [Sigma Aldrich; Merck]) on a BD Bioscience FACSFusion. The sorted plates were subsequently immediately sealed, centrifuged, and placed on dry ice to flash-freeze the cell lysate and stored at -80 'C. A modified SMART-seq2 protocol57 was performed on the single flow cytometry sorted-cells as previously described58. After cDNA generation, libraries were prepared (384 cells per library) using the Illumina Nextera XT kit (Illumina). Each library was sequenced to achieve a minimum depth of 1-2 million raw reads per cell using an Illumina HiSeq 4000 using v. 4 SBS chemistry to generate 75-bp paired-end reads.

  Study EGAS00001006885

• Covid19 WGS variants analysis

  483 samples were collected for WGS and aligned with GRCh38 human genome. The variants were called using GATK (Sentieon v. 201808.03) in GVCF and VCF format.

  Dataset EGAD00001010138

• GCAT| WGS BAM V1

  This dataset include BAM files of 808 samples from GCAT cohort. Technology used HiSeq 4000, read length 150 bp, inner mate distance 300 bp. For each sample the paired -ends are generated in separated files. Each FASTQ is splitted in multiple LANEs and grouped by the Multiplex index.

  Dataset EGAD00001008202

• Irish Covid19 WGS Raw Reads

  Dataset contains 483 Irish origin Individuals with Covid19. paired end sequencing has been performed For WGS and later processed by using sention v 201808.03.

  Dataset EGAD00001010186

• Characterization of individual foci of multicentric/multifocal breast cancer using targeted next generation sequencing

  Multifocality or multicentricity in breast cancer may be defined as the presence of two or more tumor foci within a single quadrant of the breast or within different quadrants of the same breast, respectively. This original classification of the breast cancer as multicentric or multifocal was based on the assumption that cancers arising in the same quadrant were more likely to arise from the same ductal structures than those occurring in separate areas of the breast. The problem with these definitions is that the ?quadrants? of the breast are arbitrary external designations, as no internal boundaries do exist. This project will therefore focus both on synchronous multifocal and multicentric tumors. The incidence of multifocal and multicentric breast cancers was reported to be between 13 and 75% depending on the definition used, the extent of the pathologic sampling of the breast and whether in situ disease is considered evidence of multicentricity (1). Although this incidence is variable, those figures show that it is a frequent phenomenon. Multiple (multifocal/multicentric) breast carcinomas, especially when occurring in the same breast, represent a real challenge for both pathologists and clinicians in terms of identifying the cellular origin and the best therapeutic management of the cancer. Multifocality or multicentricity has been associated with a number of more aggressive features including an increased rate of regional lymph node metastases and adverse patient outcome when compared with unifocal tumors (2-3), and a possible increased risk of local recurrence following breast conserving surgery (4). For the moment, the literature is divided on whether there is a corresponding impact on survival outcomes. Today, the current convention to stage and to treat multifocal and multicentric tumors is the classical tumor-node-metastasis (TNM) staging guidelines with which tumor size is assessed by the largest tumor focus without taking other foci of disease into consideration. If some papers, as the recent one from Lynch and colleagues, support the current staging convention (3), others, however, as Boyages et al. suggested that aggregate size and not the size of the largest lesion should be considered in order to refine the prognostic assessment of those tumors (5). On the top of that, the question whether multifocal/multicentric carcinomas are due to the spread of a single carcinoma throughout the breast or is due to multiple carcinomas arising simultaneously has been a matter of debate. Some studies suggested that multifocal breast cancer may result from either intramammary spread from a single primary tumor or multiple synchronous primary tumors; whereas others suggest that multiple breast carcinomas always arise from the same clone (6-8). Recently, Pietri and colleagues analyzed the biological characterization of a series of 113 multifocal/multicentric breast cancers (8) which were diagnosed over a 5-year period. The expression of estrogen (ER) and progesterone (PgR) receptors, Ki-67 proliferative index, expression of HER2 and tumor grading were prospectively determined in each tumor focus, and mismatches among foci were recorded. Mismatches in ER status were present in 5 (4.4%) cases and PgR in 18 (15.9%) cases. Mismatches in tumor grading were present in 21 cases (18.6%), proliferative index (Ki-67) in 17 (15%) cases and HER2 status in 11 (9.7%) cases. Interestingly, this heterogeneity among foci has led to 14 (12.4%) patients receiving different adjuvant treatments compared with what would have been indicated if we had only taken into account the biologic status of the primary tumor. This study therefore showed that differences in biological characteristics of multifocal/multicentric lesions play a crucial role in the adjuvant treatment decision making process. In this study, we will concentrate on a larger series of patients with multifocal invasive ductal breast cancer lesions. We aim at: 1. Evaluating the incidence of multifocality according to the different breast cancer molecular subtypes (ER-/HER2-, HER2+, ER+/HER2-). 2. Evaluating the incidence of multifocality in patients with hereditary breast cancer disease (presence of germline BRCA1 or BRCA2 mutations). Moreover, we would like to investigate if multifocal lesions with BRCA1 or BRCA2 mutations exhibit a characteristic combination of substitution mutation signatures and a distinctive profile of deletions as demonstrated recently by Nik-Zainal and colleagues (9). 3. Correlating multifocality with clinical information in order to define its influence on patients? survival (DFS and OS). 4. Carrying high coverage targeted gene sequencing of driver cancer genes and genes whose mutation is of therapeutic importance in order to compare clinically-relevant genetic differences between several multifocal breast cancer lesions. 5. Evaluating the impact of the distance between the different lesions on the clinical outcome but also on the genetic differences. 6. Comparing gene expression patterns between several multifocal breast cancer lesions and correlate them with the results of the targeted genes screen. 7. Characterizing the genomic and transcriptomic status of cancer related genes in metastatic lesions (local recurrence, positive lymph node or distant metastatic sites) from the same multifocal invasive ductal breast cancer patients in order to evaluate the consequence of genomic and transcriptomic heterogeneity of multifocal lesions on metastatic lesions. Multiple (multifocal/multicentric) breast carcinomas, especially when occurring in the same breast, represent a real challenge for both pathologists and clinicians in terms of identifying the cellular origin and the best therapeutic choice. This project has the potential to identify genetic/transcriptomic differences existing between several lesions constituting multifocal breast cancers, which in the routine clinical practice are usually considered to be homogeneous among them. We foresee validating significant results in a larger series of patients and this, in turn, could have a remarkable impact on the treatment and clinical management of multifocal breast cancers. Indeed, we hope to provide some evidence whether or not each focus matters in multifocal and multicentric breast cancer to define the adequate therapeutic approach, especially in the context of targeted therapies. The work to be done at Sanger will be target gene screen pooling of 1400 samples.

  Dataset EGAD00001000624

• Effect of ETS2 modulation on chromatin accessibility and enhancer activity in inflammatory (TPP) macrophages

  This dataset consists of H3K27ac ChIP-seq and ATAC-seq in primary human macrophages following modulation of ETS2 expression. There are three datasets: H3K27ac ChIP-seq in ETS2-edited and unedited inflammatory (TPP) macrophages, H3K27ac ChIP-seq in ETS2-overexpressing and control M0 macrophages, and ATAC-seq in ETS2-edited and unedited inflammatory macrophages. For editing experiments, monocytes were transfected with CRISPR-Cas9 ribonucleoproteins containing a gRNA that was validated to edit ETS2, or a non-targeting control gRNA. Inflammatory macrophage differentiation was performed using conditions that model chronic inflammation (TPP): 3 days GM-CSF (50ng/mL) followed by 3 days GM-CSF, TNFa (50ng/mL), PGE2 (1mg/mL), and Pam3CSK4 (1mg/mL). For overexpression experiments, controlled overexpression of ETS2 mRNA or control mRNA (an equivalent amount of mRNA encoding the reverse complement of ETS2 – thereby controlling for the quantity, length and purine/pyrimidine composition of the transfected RNA but with a transcript that would not be translated) was induced in resting, non-activated (M0) macrophages by transfecting 500ng mRNA. To minimise non-specific activation due to the transfected RNA, in vitro transcription was performed using co-transcriptional capping (to minimise uncapped products), and incorporating modified, minimally immunogenic nucleotides (replacing uridine with N1-methyl-pseudouridine and cytidine with methylcytidine). After 18 hours, transfected cells were were activated with low dose LPS and harvested 6 hours later. For H3K27ac ChIP-seq, cells were harvested, cross-linked, quenched, lysed, and sheared. Immunoprecipitation of histone-DNA complexes was performed overnight at 4C with rotation using an anti-H3K27ac antibody and the SimpleChIP Plus Sonication ChIP kit (Cell Signaling Technology). Following reverse cross-linking, 50ng of immunoprecipitated DNA or input DNA were used to prepare sequencing libraries using the iDeal Library Preparation kit (Diagenode), according to manufacturer instructions. 10 PCR cycles were used for the amplification step and size selection was not performed. The quality and molarity of all libraries was assessed using a BioAnalyzer 2100 (Agilent) and the libraries were sequenced on an Illumina NovaSeq 6000. Raw data are provided as raw and aligned sequencing reads from H3K27ac-bound DNA and the input chromatin. ATAC-seq in TPP macrophages was performed using the Omni-ATAC protocol (Corces et al 2017) with the following modifications: cell number was increased to 75,000 cells, cell lysis time was increased to 5 minutes; volume of Tn5 transposase in the transposition mixture was doubled; duration of the transposition step was extended to 40 minutes. Sequencing libraries were prepared by PCR with the NEBNext Ultra II Q5Master Mix (NEB) as described previously (Corces et al 2017). Amplified libraries were cleaned using AMPure XP beads (Beckman Coulter). Sequencing reads were trimmed using TrimGalore! (Phred score 24), filtered to remove reads < 20 bp and aligned to the human genome (hg19) using Burrows-Wheeler Aligner (BWA). Aligned reads were converted to BAM files, sorted and indexed using SAMtools (Li et al, 2009). PCR duplicates and unmapped reads were removed, and the resulting BAM files were re-sorted and indexed. Raw data are provided as 100 bp paired-end Illumina reads from n = 3 donors for each experiment. BAM files are also provided.

  Study EGAS50000000109

• Small molecule inhibitors of LOXL synergize with 5-AZA to restore erythropoiesis in myeloid neoplasms

  Objectives: Bone marrow (BM) stroma and microenvironmental factors in the BM milieu of myeloid neoplasms are increasingly being recognized as novel perspective therapeutic targets. Lysyl oxidases (LOX/LOXL) are crucial enzymes that cross-link collagen and contribute to the deposition of aberrant stiff extracellular matrix (ECM) in the BM of patients with myelodysplastic syndromes (MDS) and myeloproliferative neoplasms (MPN). Here, we studied whether inhibition of LOX/LOXL enzymes produced by BM-derived MSCs favors the differentiation of MDS/MPN hematopoietic stem cell and progenitor cells (HPSCs) and enhances the effects of 5-azacytidine (5-AZA). Methods: Primary co-cultures of MSCs and HPSCs were established from n=31 MDS/MPN patients and treated with the pan-LOX/LOXL inhibitor PXS-5505, 5-AZA or the combination of both for 4 days. At the end of co-culture, HSPC differentiation was assessed by CFU assays and flow cytometry. In some of the experiments, MSCs monolayers were pre-treated or not with PXS-5505 for 6 days and decellularized using 20 mM ammonia. Subsequently, autologous HSPCs were cultured on ECM matrices in the presence of 5-AZA or PXS-5505 + 5-AZA combination with or without inhibitors of integrin signaling such as blebbistatin, Y-27632 and αvβ3-integrin blocking antibodies. Patient-derived xenografts (PDX) were established in NSG mice based on intrafemoral co-injection of HSPCs and autologous MSCs. PDX were treated with 3 cycles of 5-AZA (75 mg/m2), PXS-5505 (30 mg/kg on alternate days for the duration of 5-AZA treatment) or combination of both. Results: MSCs of MDS and MPN patients displayed significant overexpression of LOX/LOXL transcripts as compared to healthy age-matched controls. Remarkably, LOX/LOXL inhibition in MSCs/HSPCs co-cultures increased erythroid differentiation of patient HSPCs that had suboptimal or no response to 5-AZA alone. From a total of n=31 advanced stage MDS/MPN patients, n=11 (35.5%) responded to the combination treatment, which was evidenced by a significant increase of CD235a+CD45- erythroid cell production compared to 5-AZA/PXS/untreated arms (p&lt;0.05). Of note, combination treatment with 5-AZA + PXS-5505 facilitated differentiation of CD235+ erythroid progenitors from HSPCs sub-clones with low mutational VAFs. The presence of MSCs was necessary for the effects of the combination treatment and required direct cell contacts between MSCs and HSPCs. Moreover, PXS-5505-mediated alterations of the ECM composition alone were sufficient to induce erythroid differentiation. Inhibition of integrin-mediated signaling using blebbistatin, Y-27632 and αvβ3-integrin antibodies completely abrogated the effects of PXS-5505 induced erythroid differentiation. In vivo, the combination of PXS-5505 + 5-AZA favored erythroid differentiation of HSPCs and confirmed the in-vitro results. The combination treatment of 5-AZA + PXS-5505 induced the greatest reduction of disease burden as assessed by decrease of spleen indexes, total engraftment rates, BM fibrosis and mutational VAF profiles. Conclusions: The combination of 5-AZA + LOX/LOXL inhibitor PXS-5505 in MDS and MPN synergistically restores erythropoietic differentiation of primary HSPCs. The effect is dependent of contact interaction between HSPCs and niche components, such as MSCs and MSCs-derived ECM. Our data propose a strong rationale for a BM niche targeted combination treatment with 5-AZA + PXS-5505 in transfusion dependent MDS and MPN, which will be assessed in forthcoming clinical trials.

  Study EGAS00001006174

• Targeted_sequencing_of_genes_recurrently_mutated_in_AML

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

  Study EGAS00001000408

• Prognostic and therapeutic significance of leukemia subtypes and minimal residual disease measurements in pediatric acute lymphoblastic leukemia treated with contemporary risk-directed trial: a cohort study

  Background Genomic analyses to identify leukemia subtypes and refinement of minimal residual disease (MRD) determination are two recent advances that would have substantial impact on the management of children with acute lymphoblastic leukemia (ALL). However, the combined application of these two approaches to determine the prognostic and therapeutic significance of MRD measurement and full spectrum of leukemia subtypes in the context of contemporary risk-directed therapy has not been studied comprehensively.MethodsBetween October 29, 2007 and March 26, 2017, 598 consecutive children with newly diagnosed ALL were enrolled in a clinical trial at St. Jude Children’s Research Hospital. Patients were provisional classified according to the presenting clinical and laboratory features including immunophenotype and cytogenetics. Patients with ETV6-RUNX1 and high hyperdiploid B-progenitor ALL were provisionally included in the low-risk groups; those with TCF3-PBX1, hypodiploidy and T-cell ALL in the standard (intermediate)-risk group; and patients with BCR-ABL1, KMT2A and early T-cell precursor ALL in the high-risk group. Final risk assignment was based primarily on MRD levels measured by flow cytometry on Day 15 (≥1%) and Day 42 (&lt;0.01% or ≥1%) of remission induction. Using genomic classification based on cytogenetics and genomic analysis (single nucleotide polymorphism array analysis and transcriptome, exome, and genome sequencing), patients were classified into 17 subtypes. The primary aim was to determine the event-free survival, overall survival and cumulative risk of relapse according to specific leukemia subtypes and MRD levels on Days 8, 15 and 42 of remission induction. This clinical trial was registered at ClinicalTrials.gov, number NCT00549848.FindingsPatients with ETV6-RUNX1, high hyperdiploidy and DUX4-rearranged B-ALL all received low- or standard-risk therapy, and had the highest 5-year event-free survival rates: 98.4% (95% CI, 95.9-100), 95.3% (91.2-99.4), and 95.0% (84.2-100), respectively. All 142 patients with Day 8 MRD&lt;0.01%, except two with KMT2A-rearranged ALL and one with TCF3-PBX1, were in continuous complete remission. The TCF3-PBX1, PAX5alt, T-cell, early T-cell precursor, iAMP21 and hypodiploid ALL had intermediate outcome with 5-year event-free survival rates ranging from 80.0% (39.4-100) to 88.2% (71.7-100). Patients with KMT2A-rearranged, BCR-ABL1, BCR-ABL1-like and ETV6-RUNX1-like ALL had the worst outcomes, with 5-year event-free survival rates between 64.1% (43.9-84.3) and 76.2% (51.9-100), despite most patients receiving standard- or high-risk therapy. Treatment intensification based on MRD and genotype improved outcome. Relapse did not occur in any of the patients with ETV6-RUNX1, DUX4-rearranged, TCF3-PBX1, iAMP21, ZNF384-rearranged, hypodiploid, or BCR-ABL1-like B-ALL who had day 15 MRD≥1% and received standard-risk therapy. However, day 42 MRD&lt;0.01% did not preclude relapse in patients with intermediate-risk or unfavorable subtypes including TCF3-PBX1, PAX5alt, T-cell, iAMP21, BCR-ABL1, BCR-ABL1-like, ETV6-RUNX1-like, KMT2A-rearranged, and MEF2D-rearranged ALL, despite intensified treatment in the standard- or high-risk arms. InterpretationContemporary risk-directed treatment based on genotypes and MRD can cure virtually all patients with ETV6-RUNX1, hyperdiploid&gt;50 and DUX4-rearranged ALL. While intensified treatment improve outcome for patients with high MRD during induction, substantial proportions of patients with intermediate-risk or unfavorable subtypes would still relapse even achieving negative MRD status at the end of induction and require novel therapeutics to improve outcome.

  Study EGAS00001004739

• Exome sequencing reveals pathogenic variants in known and novel candidate genes for severe sperm motility disorders

  STUDY QUESTION: What are the causative genetic variants in patients with male infertility due to severe sperm motility disorders? SUMMARY ANSWER: We identified high confidence disease-causing variants in multiple genes previously associated with severe sperm motility disorders in 9 out of 21 patients (43%). We also identified variants in novel candidate genes in 8 additional patients (38%). WHAT IS KNOWN ALREADY: Severe sperm motility disorders are a form of male infertility characterised by vital, but immotile sperm often in combination with a spectrum of structural abnormalities of the sperm flagellum. Examples of this spectrum of disorders include dysplasia of the fibrous sheath / multiple morphological abnormalities of the sperm flagella (DFS-MMAF) and severe forms of asthenozoospermia associated to non-specific structural defects. These disorders frequently recur in families and may include chronic respiratory disease. As such, genetic causes are suspected and many genes have recently been reported. Currently, depending on the clinical sub-categorisation, up to 50% of causality in patients with severe sperm motility disorders can be explained by pathogenic variants in at least 21 known genes. STUDY DESIGN, SIZE, DURATION: We performed exome sequencing in patients with severe sperm motility disorders from two different clinics, aiming to identify the underlying genetic defects and establishing a diagnostic yield in both. PARTICIPANTS/MATERIALS, SETTING, METHOD: Two groups of infertile men, one from Argentina (n= 9) and one from Australia (n=12), with clinically defined sperm motility disorders were included. All patients in the Argentine cohort were diagnosed with DFS-MMAF, based on transmission electron microscopy. Exome sequencing was performed in all 21 patients and variants with an allele frequency of &lt;1% in the gnomAD population were prioritised and interpreted. MAIN RESULTS AND ROLE OF CHANCE: In 9/21 patients (43%), we identified pathogenic variants in known sperm motility disorder genes: CFAP43 (3 patients); CFAP44 (2 patients), QRICH2 (2 patients), DNAH1 (1 patient) and DNAH6 (1 patient). The diagnostic rate did not differ markedly between the Argentinian and the Australian cohort (44% and 42%, respectively). Furthermore, we identified patients with variants in the novel candidate sperm motility genes CFAP58, DNAH12, DRC1, MDC1, PACRG, SSPL2C and TPTE2. One patient presented with two possibly damaging variants in four candidate genes and it remains unclear which variants were responsible for the severe sperm motility defect in this patient. LIMITATIONS, REASONS FOR CAUTION: In this study, we described patients with potential bi-allelic variants in known and novel candidate genes for severe sperm motility disorders. Due to unavailability of parental DNA, we have not assessed de novo or maternally inherited dominant variants and could not phase the mutations to establish in all cases that the mutations occur on both alleles. WIDER IMPLICATIONS OF THE FINDINGS: Our results confirm an important role for five known genes for sperm motility and we demonstrate that exome sequencing is an effective method to diagnose patients with severe sperm motility defects (9/21 diagnosed; 43%). Furthermore, our analysis revealed seven novel candidate genes for severe sperm motility disorders. Genome-wide sequencing of additional patient cohorts and re-analysis of exome data of currently unsolved cases may reveal additional variants in these novel candidate genes.

  Study EGAS00001005018

• Recurrent DNMT3B gene rearrangements are associated with unfavorable outcome in dicentric (9;20)-positive pediatric BCP-ALL

  The dicentric chromosome dic(9;20)(p11∼13;q11), occurs in 2% of pediatric acute lymphoblastic leukemia (ALL) cases. However, the prognostic significance of the dic(9;20) alterations and the mechanism in which these alterations drive leukemogenesis remains elusive. Our study describes the demographic, clinical, prognostic, and molecular characteristics of the dic(9;20)-positive ALL in a cohort of 57 pediatric BCP-ALL patients. Targeted RNA-seq and Array CGH analysis unraveled heterogeneous breakpoints on chromosomes 9 and 20 and frequent deletions of the IKZF1 gene. Together with the intrinsic susceptibility for deletions in CDKN2A, CDKN2B, and PAX5 genes in the dic(9;20)-positive ALL, these patients are eligible for treatment escalation in the AIEOP-BFM ALL 2017 clinical study. Survival analysis revealed poor 5-year event-free survival (pEFS) of 69% in dic(9;20)-positive ALL cases, compared to BCP-ALL patients enrolled in modern treatment studies. Strikingly, we identified six dic(9;20)-positive cases with rearrangements involving the DNMT3B gene, and poor 5-year pEFS of 25% (SE = 20%), compared to 79% (SE = 9%) for the remaining cohort (P = 0.011). Despite the loss of the methyltransferase domain of the DNMT3B gene in most cases with DNMT3B rearrangements, we did not observe major changes in the global methylation patterns. Differential gene expression analysis of DNMT3B-rearranged ALL unraveled a dysregulation of genes involved in hematopoiesis, lymphocyte maturation, and the PI3K/AKT signaling pathway. After confirmation in an independent cohort, DNMT3B-rearranged ALL cases should be considered high-risk for relapse and treated accordingly. We provide clinical data sets of Array CGH, targeted RNA-seq, total RNA-seq, whole genome bisulfite (WGBS) and whole genome DNA sequencing (WGS) obtained from bone marrow (BM) or peripheral blood (PB) mononuclear cells of 57 pediatric patients with dicentric chromosome dic(9;20) positive Acute lymphocytic leukemia (ALL), from which in 6 cases DNMT3B gene rearrangement was identified. This data is complemented by total RNA-seq and WGBS of samples from 4 additional ALL patients with a t(12;21) translocation and ETV6-RUNX1 gene fusion. DNA was isolated from BM or PB lymphocytes using Qiagen QIAamp DNA Blood Midi Kit to perform i) Array CGH of 58 dic(9;20) positive samples by hybridizing 500ng DNA using a Agilent 400K SurePrint G3 Custom CGH Human Genome Microarray (e-Array design 84704) ii) WGBS of 6 DNMT3B rearrangement positive samples and 4 ETV6-RUNX1 positive samples using Tecan TrueMethyl oxBS-Seq module for library preparation and Illumina NovaSeq 6000 platform to run 2x151 cycles iii) WGS of DNMT3B rearrangement positive samples using Illumina Lotus DNA Library Prep Kit followed by sequencing running 2x160 cycles on an Illumina NovaSeq 6000 platform. RNA was isolated from PB lymphocytes using the PerkinElmer Chemagic 360 instrument, followed by i) targeted RNA-seq of 56 dic(9;20) positive samples prepared using Illumina TruSight RNA Pan-Cancer Panel and sequenced on an Illumina MiSeq platform running 2x75 cycles ii) total RNA-seq of 6 DNMT3B rearrangement positive samples and 4 ETV6-RUNX1 positive samples utilizing TruSeq Stranded Total RNA Library Prep Gold kit and running 2x100 cycles on an Illumina NovaSeq 6000 platform.

  Study EGAS00001007383

• The British Autozygosity Populations BioResource (2022-04-26)

  This data set contains whole exome sequences of individuals from 8086 (mostly British Pakistani/British Bangladeshi, mostly self-reported parentally related) individuals from the following studies: 1. 5236 British Pakistani/British Bangladeshi adults from East London Genes &amp; Health, now known as Genes &amp; Health 2. 2624 British South Asian mothers from Born in Bradford (mostly Pakistani) 3. 1061 British South Asian adults from Birmingham (mostly Pakistani) This dataset contains all the exome sequence data available for this study on 2022-04-26

  Dataset EGAD00001008736