14347 results for ""Pennsylvania+German"+ACS+study+'Socius'+open+access"

in 285.10 milliseconds.

• Clonal Lineage Tracing of Primary Human Cortical Progenitors, Cell Type Profiles in the Brain Vasculature and Genotyping Data for 22q11DS and Control iPS Lines

  We sought to characterize the cell types produced by individual human cortical progenitors and describe their lineage relationships. We labeled individual human cortical cells with a massively complex lentivirus library encoding a GFP reporter gene and a molecular barcode that could be captured by standard 3' scRNA-seq methodologies. We then cultured the labeled human cells for 6 weeks in either a 1) mouse cortical astrocyte co-culture system or 2) in the mouse cortex as a xenograft. In a separate experiment, we maintained human cortical tissue slices in organotypic slice culture for 12 days. GFP+ cells were then isolated by FACS and then captured on the 10x chromium scRNA-seq platform. An additional sample was labeled by lentivirus and briefly cultured for 3 days prior to being mixed with mouse 3T3 cells in order to perform a cell mixing experiment ("barnyard") to benchmark our methodologies. We found that human individual cortical progenitors gave rise to both excitatory and inhibitory neurons as well as glia.We also sought to characterize the adult brain vasculature in both vascular malformations and non-malformed tissue. We use the 10x chromium scRNA-seq platform across 10 individuals (5 control and 5 arteriovenous malformations). We find broad disregulation across multiple cell types including the immune cells and find a specific monocyte subtype to contribute to a rupture phenotype. We also find a novel population of fibromyocytes in both conditions.Finally, we confirmed the genotype status of iPS lines that were derived from either karyotypically normal control individuals or from DiGeorge Syndrome patients carrying the 22q11.2 deletion using the Illumina Infinium Global Screening Array v3.0.

  Study phs002624

• Enhancer profiling identifies epigenetic markers of endocrine resistance and reveals therapeutic options for metastatic castration-resistant prostate cancer patients (LuCaP cell line data)

  Androgen Receptor (AR) signaling inhibitors, including enzalutamide, are treatment options for patients with metastatic castration-resistant prostate cancer (mCRPC), but resistance inevitably develops. Using metastatic samples from a prospective phase II clinical trial, we epigenetically profiled enhancer/promoter activities with H3K27ac chromatin immunoprecipitation followed by sequencing, before and after AR-targeted therapy. We identified a distinct subset of H3K27ac-differentially marked regions that associated with treatment responsiveness. These data were successfully validated in mCRPC patient-derived xenograft models (PDX). In silico analyses revealed HDAC3 as a critical factor that can drive resistance to hormonal interventions, which we validated in vitro. Critically, we identified the pan-HDAC inhibitor, vorinostat to be effective in decreasing tumor cell proliferation both in vitro and in vivo. Moreover, we uncovered evidence for HDAC3 working together with Glucocorticoid Receptor (GR) as a potential mechanism for this therapeutic effect. These findings demonstrate rationale for new therapeutic strategies including HDAC inhibitors to improve patient outcome in advanced stages of mCRPC. Project Goal: We identified H3K27ac regions associated with resistance in patient samples. Extensive analysis of these regions suggested that HDAC3 may play a role in resistance. To test this, we analyzed ChIP-seq and RNA-seq from cell lines generated from PDX LuCaP samples (recapitulating resistant and non-resistant) to validate findings in patient samples. ChIP-seq samples were used to verify the signals identified from specific factors, such as H3K27ac, AR, FOXA1, GR and HDAC3, their association with resistance and possible mechanisms of action. RNA-seq samples were analyzed to determine the effect of AR signaling in these samples in specific genes in the H3K27ac resistance-associated regions.

  Study EGAS50000000917

• Cellular Profiling Identifies Targetable T Cell Phenotypes in Lymphocytic Variant Hypereosinophilic Syndrome

  Lymphocytic variant hypereosinophilic syndrome (LHES) is driven by aberrant Th2‑skewed lymphocyte clones that promote tissue eosinophilia. We performed multimodal single‑cell RNA‑sequencing (scRNA‑seq), cell‑surface protein (ADT) profiling, and paired αβ‑TCR sequencing on peripheral‑blood mononuclear cells (PBMCs) from two clinically characterized LHES patients and one idiopathic HES control. One LHES patient was sequenced at two timepoints, pre-JAK inhibitor treatment and post-JAK inhibitor treatment. Single‑cell regulatory‑network inference (SCENIC), CellChat ligand–receptor modelling, and scRNA variant calling were integrated with flow‑cytometric validation and 5,000× amplicon sequencing. Two clonally expanded CD3‑CD4⁺ central‑/effector‑memory T‑cell clusters were observed exclusively in LHES. These cells showed Th2‑associated regulons (GATA3, BATF, MAF), high expression of CCR4, CD52, PTGDR2 and IL13, and down‑regulation of CD247 (CD3ζ), a mechanism experimentally confirmed to abolish surface CD3/TCR. Autocrine IL‑7 signaling, predicted by CellChat, provides an antigen‑independent proliferative cue. Therapeutically actionable targets (CCR4, CD52, IL‑5, IL‑7, TNF) were prioritized; one refractory patient achieved complete clinical and cytologic remission after mogamulizumab (anti‑CCR4) therapy. Longitudinal sequencing revealed clonal evolution with acquisition of a STAT3 p.G618R mutation under JAK‑inhibition, highlighting a resistance mechanism. Raw scRNA‑seq (5′ cDNA, VDJ, ADT) data, and sample‑level metadata will be distributed through dbGaP. These resources provide the first high‑resolution immune atlas of LHES and a framework for rational, precision immunotherapy in eosinophilic disorders.

  Study phs004041

• Convergent evolution towards CD38 biallelic loss is a recurrent mechanism of resistance to anti-CD38 antibodies in multiple myeloma.

  Monoclonal antibodies directed against CD38 are a therapeutic mainstay in multiple myeloma (MM). While they have contributed to improved outcomes, most patients will still experience disease relapse. Despite their ubiquity, little is known about tumor-intrinsic mechanisms of resistance to these drugs. Antigen escape, whereby tumor cells alter expression of cell surface immunotherapy targets, has been proposed as a potential avenue, and has been seen with other immunotherapies in MM. Yet, it is unknown whether MM cells can develop permanent resistance to anti-CD38 antibodies by acquiring genomic events leading to biallelic disruption of the CD38 gene locus. Here, by using whole genome and whole exome sequencing data from 701 newly diagnosed patients, 72 patients at relapse with naivety to anti-CD38 antibodies, and 52 patients collected at relapse following anti-CD38 antibodies, we report biallelic inactivation of CD38 as a recurrent mechanism of resistance in 6% of patients, often via convergent evolution where distinct subclones are advantaged towards the same selective endpoint within the same patient. Monoallelic loss of the CD38 locus was seen in another 14% post-CD38-directed therapy. Three of the cases with biallelic events on CD38 presented with missense mutations and retained CD38 expression. Functional work demonstrated that two of these mutations decreased binding affinity and antibody-dependent cellular cytotoxicity of commercial antibodies Daratumumab and Isatuximab. However, the third mutation (R140G) provided relative resistance only towards one commercial product. Given that patients are often rechallenged with CD38 antibodies following disease progression, these data support the need for next generation sequencing to guide treatment selection.

  Study EGAS50000000709

• IDH-wildtype untreated human glioblastoma samples (GB-UK cohort), published in Noorani & Haughey et al 2025.

  Extrachromosomal DNA (ecDNA) oncogene amplification is associated with treatment resistance and shorter survival in cancer. Currently, the spatial dynamics of ecDNA, and their evolutionary impact, are poorly understood. Here, we investigate the spatiotemporal evolution of ecDNA by integrating computational modeling with samples from 94 treatment-naive human IDH wild-type glioblastoma patients. Random ecDNA segregation combined with ecDNA-conferred fitness advantages induce distinct spatial ecDNA copy number patterns which depend on ecDNA oncogenic makeup. EGFR-ecDNAs often reach high copy numbers and confer strong fitness advantages. In contrast, PDGFRA-ecDNAs reach lower copy numbers and confer weaker fitness advantages. EGFR structural variants occur exclusively on ecDNA, arise from and are intermixed with wild-type EGFR-ecDNAs. Computational modeling suggests wild-type and variant EGFR-ecDNAs often accumulate before clonal expansion. Using a genetically engineered mouse model, we confirm that the accumulation of oncogenic ecDNA prior to clonal expansion under strong positive selection is observable both in vivo and in vitro in mouse neural stem cells. Intermixed wild-type and variant ecDNAs remain prevalent at diagnosis, as less fit wild-type ecDNAs hitchhike on strongly selected variant EGFR-ecDNAs. A dN/dS analysis in the Genomics England cohort further suggests positive selection of variants on EGFR-ecDNAs but not PDGFRA-ecDNAs in glioblastoma. Overall, our results suggest ecDNA oncogenic makeup determines unique evolutionary trajectories in glioblastoma. New concepts like ecDNA tumor clonality, ecDNA heteroplasmy and multi-species ecDNAs require a refined evolutionary interpretation of genomic data in a large subset of glioblastoma patients and likely other solid cancers.

  Study EGAS00001008126

• RNAseq: Acquired non-permissive bone marrow microenvironment impairs hematopoietic stem cell proliferation and maintenance and B-cell development post-HSCT

  Defects in B-cell reconstitution upon hematopoietic stem-cell (HSC) transplantation (HSCT) are a common observation, yet the mechanism remains unexplained. The bone marrow (BM) stroma, including mesenchymal stromal cells (MSCs), guides HSC maintenance and B-lymphopoiesis by secreting crucial cytokines. We report acquired, permanent, selective and complete B-cell deficiency in the context of full donor-chimerism in a patient with X-linked lymphoproliferative disease and aimed to identify the contribution of the BM-microenvironment in disrupted B-cell reconstitution post-HSCT. We studied longitudinal BM samples from the patient and his identical twin, both of whom underwent HSCT with the same donor with opposite outcomes in B-cell reconstitution. In the index patient BM, we observed progressive loss of proliferation of HSCs and a selective block at the pre-BI cell stage. In vitro modeling studies showed limited survival of patient-HSCs and a relative accumulation of pre-B cells. Patient-derived MSCs failed to support survival and proliferation of HSCs and B-cell development of healthy-HSCs which was correlated with reduced CXCL12 levels. Using bulk RNA-sequencing of MSCs and in vitro functional studies, we showed global changes in the patients’ MSCs and a progressive loss of CXCL12 expression. Indeed, survival of patients HSCs improved supplementing in vitro development culture with CXCL12, suggesting a contribution of defective CXCL12 signaling to the phenotype. In summary, our data show that an acquired defect in the BM-stromal microenvironment and exhaustion of HSCs and committed progenitors may cause a permanent non-permissive state for normal B-cell development.

  Study EGAS50000001437

• ATAC-seq: Acquired non-permissive bone marrow microenvironment impairs hematopoietic stem cell proliferation and maintenance and B-cell development post-HSCT

  Defects in B-cell reconstitution upon hematopoietic stem-cell (HSC) transplantation (HSCT) are a common observation, yet the mechanism remains unexplained. The bone marrow (BM) stroma, including mesenchymal stromal cells (MSCs), guides HSC maintenance and B-lymphopoiesis by secreting crucial cytokines. We report acquired, permanent, selective and complete B-cell deficiency in the context of full donor-chimerism in a patient with X-linked lymphoproliferative disease and aimed to identify the contribution of the BM-microenvironment in disrupted B-cell reconstitution post-HSCT. We studied longitudinal BM samples from the patient and his identical twin, both of whom underwent HSCT with the same donor with opposite outcomes in B-cell reconstitution. In the index patient BM, we observed progressive loss of proliferation of HSCs and a selective block at the pre-BI cell stage. In vitro modeling studies showed limited survival of patient-HSCs and a relative accumulation of pre-B cells. Patient-derived MSCs failed to support survival and proliferation of HSCs and B-cell development of healthy-HSCs which was correlated with reduced CXCL12 levels. Using bulk RNA-sequencing of MSCs and in vitro functional studies, we showed global changes in the patients’ MSCs and a progressive loss of CXCL12 expression. Indeed, survival of patients HSCs improved supplementing in vitro development culture with CXCL12, suggesting a contribution of defective CXCL12 signaling to the phenotype. In summary, our data show that an acquired defect in the BM-stromal microenvironment and exhaustion of HSCs and committed progenitors may cause a permanent non-permissive state for normal B-cell development.

  Study EGAS50000001438

• Development and validation of multivariate predictors of primary endocrine resistance to tamoxifen and aromatase inhibitors in luminal breast cancer reveal drug-specific differences

  Background: Endocrine therapy is highly effective in blocking the estrogen receptor pathway in HR+/HER2- early breast cancer (EBC). However, up to 40% of patients experience relapse during or after adjuvant endocrine therapy. Here, we investigate molecular mechanisms associated with primary resistance to endocrine therapy and develop predictive models. Results: TP53 mutations were prominently associated with primary resistance to both tamoxifen (TAM) and aromatase inhibitors (AI), with AI non-responders exhibiting resistance in up to 32% of cases. Additionally, we identified distinct DNA methylation patterns in TAM and AI non-responders, with TAM non-responders showing global DNA methylation loss, associated with KRAS signaling, apical junctions and epithelial-mesenchymal transition (EMT). Conversely, we observed methylation gain in AI non-responders affecting developmental transcription factors, hypoxia and estrogen signaling. TAM or AI resistance was associated with increased methylation-inferred proportions of immune cells and decreased proportions of endothelial cells. Based on these findings and patient age, we developed the Predictive Endocrine ResistanCe Index (PERCI). PERCI stratified NR and R cases in both treatment groups and cohorts with high accuracy (ROC AUC TAM discovery 93.9%, validation 83%; AI discovery 98.6%, validation 76.9%). A simplified PERCI efficiently predicted progression-free survival in the TCGA-BRCA sub-cohort (Kaplan-Meier log-rank p-value = 0.03 between low and high PERCI groups). Conclusions: We identified genomic and epigenomic features associated with primary resistance to TMA and AI. By combining information on genomic alterations, patient age, differential methylation and tumor microenvironment (TME) composition, we developed PERCI TAM and PERCI AI as novel predictors of primary resistance, with potential additional prognostic value. Applying PERCI in a clinical setting may allow patient-specific drug selection to overcome resistance.

  Study EGAS50000001102

• CCR8-targeted specific depletion of clonally expanded Treg cells in tumor tissues evokes potent tumor immunity with long-lasting memory

  Foxp3-expressing CD25+CD4+ regulatory T cells (Tregs) are abundantly infiltrating into tumor tissues, hindering effective anti-tumor immunity. Here, assuming that such tumor Tregs were clonally expanding by recognizing tumor-associated antigens, we attempted to characterize them by their gene expression profiling and T-cell receptor (TCR) clonotyping at the single-cell level. We found that they predominantly expressed the chemokine receptor CCR8: ~40% of tumor Tregs, less than 10% of Tregs in other tissues, and few tumor-infiltrating conventional T cells expressed CCR8 in mice and humans. CCR8+ tumor Tregs were phenotypically highly differentiated and functionally stable. One injection of anti-CCR8 mAb with high antibody-dependent cellular cytotoxic activity was indeed able to specifically deplete clonally expanded tumor Tregs and thereby eradicate established tumors in mice. The depletion predominantly expanded tumor-infiltrating CD8+ effector T cells, which were different in activated or exhausted patterns from those seen after immune checkpoint blockade by anti-PD-1 mAb. The treatment also evoked potent secondary immune responses against the same tumor cell line inoculated several months after tumor eradication, indicating that the anti-CCR8-induced memory-type tumor-reactive effector T cells were, once induced, resistant to recovered Treg suppression. In contrast to such evocation of potent tumor immunity, anti-CCR8 treatment scarcely elicited discernable autoimmunity in mice, contrasting with a similar treatment for systemic whole Treg depletion, which induced histologically and serologically evident autoimmunity. Thus, temporary and specific depletion of clonally expanded Tregs in tumor tissues by cell-depleting anti-CCR8 mAb is sufficient to evoke potent tumor immunity with long-lasting memory and without serious autoimmunity.

  Study JGAS000312

• mutation analysys of Gorlin syndrome

  Gorlin syndrome is a genetic disorder of autosomal dominant inheritance that predisposes the affected individual to a variety of disorders that are attributed largely to heterozygous germline patched1 (PTCH1) mutations. PTCH1 is a hedgehog (Hh) receptor as well as a repressor, mutation of which leads to constitutive activation of Hh pathway. Hh pathway encompasses a wide variety of cellular signaling cascades, which involve several molecules; however, no associated genotype?phenotype correlations have been reported. Recently, mutations in Suppressor of fused homolog (SUFU) or PTCH2 were reported in patients with Gorlin syndrome. These facts suggest that multi-layered mutations in Hh pathway may contribute to the development of Gorlin syndrome. We demonstrated multiple mutations of Hh-related genes in addition to PTCH1, which possibly act in an additive or multiplicative manner and lead to Gorlin syndrome. High-throughput sequencing was performed to analyze exome sequences in four unrelated Gorlin syndrome patient genomes. Mutations in PTCH1 gene were detected in all four patients. Specific nucleotide variations or frameshift variations of PTCH1 were identified along with the inferred amino acid changes in all patients. We further filtered 84 different genes which are closely related to Hh signaling. Fifty three of these had enough coverage of over ?30. The sequencing results were filtered and compared to reduce the number of sequence variants identified in each of the affected individuals. We discovered three genes, PTCH2, BOC, and WNT9b, with mutations with a predicted functional impact assessed by MutationTaster2 or PolyPhen-2 (Polymorphism Phenotyping v2) analysis. It is noticeable that PTCH2 and BOC are Hh receptor molecules. No significant mutations were observed in SUFU. Multi-layered mutations in Hh pathway may change the activation level of the Hh signals, which may explain the wide phenotypic variability of Gorlin syndrome.

  Study JGAS000099

• Immunohistochemical and molecular pathological search in gastrointestinal tumors

  Purpose: RSPO fusions are recently identified recurrent genetic alterations in colorectal cancer (CRC), which lead to WNT pathway activation and are potential therapeutic targets. We aimed to investigate the clinicopathological and molecular features of RSPO fusion-positive CRCs. Experimental Design: We screened 1019 CRCs for RSPO fusions using multiplex reverse transcription polymerase chain reaction. RSPO fusion-positive tumors were characterized for their clinicopathological features and subjected to whole exome sequencing (WES). Results: Our analysis identified 29 CRCs with RSPO fusions (2.8%), consisting of 5 with an EIF3E-RSPO2 fusion and 24 with PTPRK-RSPO3 fusions. The patients were 17 females and 12 males with the median age of 64 years. Thirteen tumors were located in the right colon, 8 were in the left colon, and 8 were in the rectum. Histologically, most of the tumors (24/29, 83%) were classified as adenocarcinomas, not otherwise specified. However, about half of the tumors (13/29, 45%) had a focal or extensive mucinous component, which is significantly more frequent compared with RSPO fusion-negative tumors (13%; p = 8.1 �� 10-7). Four tumors (14%) were mismatch repair-deficient, three of which were associated with MLH1 promoter methylation. WES identified KRAS, BRAF, and NRAS mutations in a total of 27 cases (93%). In contrast, pathogenic mutations in major WNT pathway genes other than RSPO fusions were virtually absent. The presence of an RSPO fusion tended to show shorter overall and recurrence-free survivals in stage II and III patients although the differences were not significant. Conclusions: RSPO fusion-positive CRCs constitute a rare subgroup of CRCs with several characteristic clinicopathological and genetic features. The mutual exclusivity with other WNT pathway gene alterations supports RSPO fusions as valid therapeutic targets in CRC.

  Study JGAS000538

• Genome and gene analysis of gastrointestinal cancer and elucidation of its clinicopathological significance

  Cancer develops through the accumulation of genetic and epigenetic aberrations. To identify sequential molecular alterations that occur during the development of hepatocellular carcinoma (HCC), we molecularly compared 52 early and 108 overt HCC samples by genome sequencing. Gene mutations in the p53/RB1 pathway, WNT pathway, MLLs, SWI/SNF complexes, and AKT/PI3K pathway were common in HCC. In the early phase of all entities, TERT was the most frequently upregulated gene owing to diverse mechanisms. Despite frequent somatic mutations in driver genes, including CTNNB1 and TP53, early HCC was a separate molecular entity from overt HCC as each had a distinct expression profile. Notably, WNT target genes were not activated in early HCC, regardless of CTNNB1 mutation status, because ?-catenin did not translocate into the nucleus due to the E-Cadherin/?-catenin complex at the membrane. Conversely, WNT targets were definitively upregulated in overt HCC with CTNNB1 mutation associated with downregulation of CDH1 and hypomethylation of CpG islands in the target genes. Similarly, cell cycle genes downstream of the p53/RB pathway were upregulated only in overt HCC with TP53 and/or RB1 gene mutations associated with chromosomal deletion of 4q and/or 16q. Differing from the expression profile, HCC was epigenetically distinguished into four subclasses: normal-like methylation, global-hypomethylation (favorable prognosis), stem-like methylation (poor prognosis), and CpG island methylation. These methylation statuses were globally maintained through HCC progression. Taken together, additional molecular events besides driver gene mutations cooperatively contribute to transcriptional activation of downstream targets as HCC progresses, according to the methylation status of the associated genes.

  Study JGAS000233

• Paediatric Hepatic International Tumour Trial (JPLT2: PHITT)

  Hepatoblastoma (HBL) is the most common primary liver tumor in children. Despite their diversity in histopathological phenotype and clinical outcome, the underlying molecular aberrations have not been fully understood. Here we describe a comprehensive landscape of genomic, epigenomic, and transcriptomic heterogeneity among the 163 cases of HBLs and transitional liver cell tumors (TLCTs) collected in the multicenter clinical trial JPLT-2. The total numbers of mutational events in HBLs are very few (0.52 SNVs/Mb on exonic regions) but show a positive correlation with the age of onset. While dysregulation of Wnt-pathway related genes was observed (somatic deletion or mutations of CTNNB1 and germline or somatic mutations of APC) in 83.6% of the HBL cases, TERT promoter mutations are mainly in the elderly HBLs (> 8 years old), meeting with the proposed criteria of TLCT. These TLCT cases are clustered in the hypermethylated epigenotype in promoter and enhancer regions. Copy number analyses also revealed recurrent 11p15.5 UPD/LOHs including the IGF2/H19 region (31.6%), focal deletions in 4q35 (15.8%) and focal amplifications in 2q24.3 (7.0%) predominantly in the child and infant HBLs. Motif enrichment analysis of the specific hypomethylated enhancers in these HBL subtypes suggests specific binding of ASCL2, the regulatory transcription factor for definitive endoderm and digestive systems. Transcriptome and histochemical analysis revealed that ASCL2 is upregulated not in HCC and TLCT but in child and infant HBL, suggesting their ?cell of origin? derived from the premature hepatoblast with high proliferative potential such as intestinal epithelial cells. Systematic molecular profiling of HBL will provide a fundamental understanding of hepatoblast carcinogenesis and clues for risk stratification.

  Study JGAS000236

• Biallelic variants in the non-protein coding minor spliceosome components RNU6ATAC and RNU6ATAC cause syndromic monogenic autoimmune diabetes

  Non-protein coding genes are emerging as critical contributors to the aetiology of rare diseases, providing key insights into human biology and uncovering novel disease mechanisms. We identified 7 individuals from 4 families with early-onset diabetes (diagnosed <5 years) and immune dysregulatory features caused by biallelic variants in RNU6ATAC. RNU6ATAC encodes a small nuclear RNA (snRNA) acting as a catalytic component of the minor spliceosome, a protein-RNA complex mediating splicing of ~700 genes containing U12/minor-type introns. Variant screening of the other 64 minor spliceosome genes in 276 infants with diabetes identified 12 unrelated individuals with biallelic disease-causing variants in RNU4ATAC. Biallelic pathogenic RNU4ATAC variants are known to cause a variable spectrum of clinical features, which until now did not include diabetes. Clinically, 12/19 RNU6ATAC/RNU4ATAC patients had additional immune dysregulatory features, and 50% of patients tested were islet-autoantibody positive, strongly supporting an autoimmune aetiology for their diabetes. RNA-seq in 3 individuals with biallelic RNU6ATAC variants showed a pattern of intron retention in U12 intron-containing genes similar to that seen in RNU4ATAC-individuals (n=3), supporting a shared disease mechanism. Analysis of patients’ transcriptomic, methylation and immune data revealed impaired B cell development and maturation. We conclude that biallelic RNU6ATAC variants cause a novel syndrome of early-onset autoimmune diabetes and immune dysregulation. We further show that infancy-onset diabetes is a previously unrecognised feature of RNU4ATACopathy. Our work highlights the important role of two snRNAs critical to minor spliceosome function in immune system regulation, providing novel insights into the pathogenesis of autoimmune diabetes.

  Study EGAS50000001565

• Confronting historical legacies of biological anthropology in South Africa—Restitution, redress and community-centered science: The Sutherland Nine

  We describe a process of restitution of nine unethically acquired human skeletons to their families, together with attempts at redress. Between 1925–1927 C.E., the skeletonised remains of nine San or Khoekhoe people, eight of them known-in-life, were removed from their graves on the farm Kruisrivier, near Sutherland in the Northern Cape Province of South Africa. They were donated to the Anatomy Department at the University of Cape Town. This was done without the knowledge or permission of their families. The donor was a medical student who removed the remains from the labourers’ cemetery on his family farm. Nearly 100 years later, the remains are being returned to their community, accompanied by a range of community-driven interdisciplinary historical, archaeological and analytical (osteobiographic, craniofacial, ancient DNA, stable isotope) studies to document, as far as possible, their lives and deaths. The restitution process began by contacting families living in the same area with the same surnames as the deceased. The restitution and redress process prioritises the descendant families’ memories, wishes and desire to understand the situation, and learn more about their ancestors. The descendant families have described the process as helping them to reconnect with their ancestors. A richer appreciation of their ancestors’ lives, gained in part from scientific analyses, culminating with reburial, is hoped to aid the descendant families and wider community in [re-]connecting with their heritage and culture, and contribute to restorative justice, reconciliation and healing while confronting a traumatic historical moment. While these nine individuals were exhumed as specimens, they will be reburied as people.

  Study EGAS50000000971

• Repotrectinib in NTRK fusion–positive advanced solid tumors: a phase 1/2 trial

  Early-generation TRK tyrosine kinase inhibitors (TKIs) approved for treating NTRK fusion–positive (NTRK+) solid tumors provide clinical benefit; however, resistance emerges. Repotrectinib is a next-generation ROS1/TRK TKI with a compact macrocyclic structure designed to improve durability of response. TRIDENT-1 is a registrational phase 1/2 trial assessing repotrectinib, a next-generation ROS1/TRK TKI, in adults with advanced solid tumors, including NTRK+ disease. The primary endpoint was confirmed objective response; secondary endpoints included duration of response (DOR), progression-free survival (PFS), overall survival and safety. Median follow-up ranged between 21.3 months and 25.7 months. In the TKI-naive cohort (n = 51; 95% confidence interval (CI)), the response rate was 59% (44–72); the median DOR was not estimable (NE); and the median PFS was 30.3 months (9.0–NE). In the TKI-pretreated cohort (n = 69; 95% CI), the response rate was 48% (36–60); the median DOR was 9.8 months (7.4–13.0); and the median PFS was 7.4 months (3.9–9.7). Of 30 TKI-pretreated patients with NTRK solvent front mutations, 16 had a response (53%; 95% CI: 34–72). Intracranial responses were observed in two of three TKI-naive patients and in four of six TKI-pretreated patients with measurable intracranial disease at baseline. Among all treated patients (n = 565), the most common any-grade treatment-related adverse event (TRAE) was dizziness (57%); most TRAEs were low grade; and 4% discontinued repotrectinib due to a TRAE. Overall, repotrectinib demonstrated durable systemic and intracranial responses with generally low-grade adverse events in patients with NTRK+ solid tumors, including those with previous TRK TKI treatment and solvent front mutations. These results support the use repotrectinib to treat patients with NTRK+ solid tumors. ClinicalTrials.gov identifier: NCT03093116 .

  Study EGAS50000001572

• Recurrent somatic DICER1 mutations in non-epithelial ovarian tumors

  Background: Germline truncating mutations in DICER1, an RNase III-type endoribonuclease essential for processing of microRNAs, are seen in families with the pleuropulmonary blastoma-family tumor and dysplasia syndrome; these individuals have a propensity to develop non-epithelial ovarian tumors (i.e. sex-cord stromal and germ cell tumors). Methods: We sequenced the whole transcriptomes or exomes of 14 non-epithelial ovarian tumors and noted closely clustered mutations in the RNase IIIb domain of DICER1 in four cases. Using Sanger sequencing, we then analyzed additional ovarian tumors for mutations in this region. The impact of the mutations on the enzymatic activity of Dicer1 was determined using in vitro RNA cleavage assays. Results: RNase IIIb domain DICER1 mutations were found in 30/102 non-epithelial ovarian tumors (29%), predominantly Sertoli-Leydig cell tumors (26/43, 60%), including 4/4 from individuals with germline DICER1 mutations. The mutations were restricted to metal binding sites within the RNase IIIb catalytic centres critical for miRNA interaction and cleavage, and were somatic in all 14 cases where germline DNA was available for testing. Of 266 epithelial ovarian and endometrial cancers tested, one case, an ovarian carcinosarcoma, had a hotspot mutation. In vitro, the mutant DICER proteins showed reduced RNase IIIb activity but retained RNase IIIa activity. Conclusions: Somatic missense mutations of the RNase IIIb domain of DICER1 are common in non-epithelial ovarian tumors. Genetic, pathologic, and functional evidence suggest that these mutations do not obliterate DICER1 function, rather, they alter it in specific cell types, a novel mechanism through which perturbed microRNA expression is oncogenic.

  Study EGAS00001000135

• Driver mutations in histone H3.3 and chromatin remodelling genes in paediatric glioblastoma

  Glioblastoma multiforme (GBM) is a lethal brain tumour in adults and children. However, DNA copy number and gene expression signatures indicate differences between adult and paediatric cases. To explore the genetic events underlying this distinction, we sequenced the exomes of 48 paediatric GBM samples. Somatic mutations in the H3.3-ATRX-DAXX chromatin remodelling pathway were identified in 44% of tumours (21/48). Recurrent mutations in H3F3A, which encodes the replication-independent histone 3 variant H3.3, were observed in 31% of tumours, and led to amino acid substitutions at two critical positions within the histone tail (K27M, G34R/G34V) involved in key regulatory post-translational modifications. Mutations in ATRX (alpha-thalassaemia/mental retardation syndrome X-linked) and DAXX (death-domain associated protein), encoding two subunits of a chromatin remodelling complex required for H3.3 incorporation at pericentric heterochromatin and telomeres, were identified in 31% of samples overall, and in 100% of tumours harbouring a G34R or G34V H3.3 mutation. Somatic TP53 mutations were identified in 54% of all cases, and in 86% of samples with H3F3A and/or ATRX mutations. Screening of a large cohort of gliomas of various grades and histologies (n=784) showed H3F3A mutations to be specific to GBM and highly prevalent in children and young adults. Furthermore, the presence of H3F3A/ATRX-DAXX/TP53 mutations was strongly associated with alternative lengthening of telomeres and specific gene expression profiles. This is, to our knowledge, the first report to highlight recurrent mutations in a regulatory histone in humans, and our data suggest that defects of the chromatin architecture underlie paediatric and young adult GBM pathogenesis.

  Study EGAS00001000226

• The application of RNA sequencing for the diagnosis and genomic classification of pediatric acute lymphoblastic leukemia

  Acute lymphoblastic leukaemia (ALL) is the most common childhood malignancy, and implementation of risk-adapted therapy has been instrumental in the dramatic improvements in clinical outcomes. A key to risk-adapted therapies includes the identification of genomic features of individual tumors, including chromosome number (for hyper- and hypodiploidy) and gene fusions, notably ETV6-RUNX1, TCF3-PBX1 and BCR-ABL1 in B-cell ALL (B-ALL). RNA-sequencing (RNA-seq) of large ALL cohorts has expanded the number of recurrent gene fusions recognized as drivers in ALL, and identification of these new entities will contribute to refining ALL risk stratification. We used RNA-seq on 126 ALL patients from our clinical service to test the utility of including RNA-seq in standard of care diagnostic pipelines to detect gene rearrangements and IKZF1 deletions. RNA-seq identified 86% of rearrangements detected by standard of care diagnostics. KMT2A-rearrangements (MLL), whilst usually identified, were the most commonly missed by RNA-seq as a result of low expression. RNA-seq identified rearrangements that were not detected by standard of care testing in nine patients. These were found in patients who were not classifiable using standard molecular assessment. We developed an approach to detect the commonest IKZF1 deletion from RNA-seq data and validated this using a RQ-PCR assay. We applied an expression classifier to identify Ph-like B-ALL patients. T-ALL proved a rich source of novel gene fusions which have clinical implications or provide insights into disease biology. Our experience shows that RNA-seq can be implemented within an individual clinical service to enhance the current molecular diagnostic risk classification of ALL.

  Study EGAS00001004212

• A Multifactorial Tumor and Immune Cell Profile Determines Response to Immune Checkpoint blockade in Melanoma

  The complex crosstalk between tumor and immune cells during immune checkpoint blockade mandates the development of integrated models that interpret the anti-tumor immune response and predict clinical outcome. We have integrated genome-wide sequence and structural alterations with pre and on-therapy transcriptomic and T cell repertoire features for a cohort of 64 immunotherapy-naïve melanomas treated with anti-PD1 monotherapy or combined anti-PD1 and anti-CTLA4 therapy. While tumor mutation burden (TMB) was associated with improved response to therapy, expressed mutation burden was superior to TMB in predicting outcome. An increased pre-existing T cell density differentiated responding from non-responding tumors independent of therapy. Importantly, T cell repertoire reshaping determined by T cell receptor (TCR) clonotypic regressions and expansions reflected response to therapy such that a more dynamic repertoire predicted improved clinical outcome. Through whole-transcriptome analyses, we discovered differential abundance of B cell subsets in responding tumors, which highlights the importance of the interplay between pre-existing T and B cell immunity in shaping therapeutic response. Tumor temporal trajectories during therapy revealed distinct patterns of molecular response related to expressed mutation elimination or retention that accurately interpreted clinical response. High-dimensional genomic, transcriptomic and immune repertoire data were integrated by both a machine learning and a censored regression approach, resulting in a harmonized multi-modal predictor of response to immune checkpoint blockade. B cell abundance, expressed mutation load and tumor aneuploidy were combined to identify patients at high risk for recurrence, such that high risk patients had a significantly shorter progression-free survival (HR=9.18, 95% CI: 3.14-26.85, p=3.4e-06) especially in the anti-PD1/anti-CTLA4 group (HR=20.34, 95% CI: 4.08-101.32, p=1.68e-06).

  Study EGAS00001004548

• Comparison of capture-based method for transcriptome profiling of formalin-fixed paraffin embedded tumor samples

  Background: The need for fresh frozen (FF) tissue limits implementing RNA sequencing (RNA-seq) in the clinic. The majority of clinical samples are stored as formalin-fixed, paraffin-embedded (FFPE) tissues. Exome capture platforms have been developed for RNA-seq from FFPE samples. However, these methods have not been systematically compared. Methods: Transcriptomic analysis of 32 FFPE tumor samples from 11 patients was performed using three exome capture-based methods: Agilent SureSelect V6, TWIST NGS Exome, and IDT XGen Exome Research Panel. We compared these methods to TruSeq RNA-seq of fresh frozen (FF-TruSeq) tumor samples from the same patients. We assessed the recovery of clinically relevant biological features. Results: The Spearman's correlation coefficients between the global expression profiles of the three capture-based methods from FFPE and matched FF-TruSeq were high (rho = 0.72-0.9, p &lt; 0.05). A significant correlation between the expression of key immune genes between individual capture-based methods and FF-TruSeq (rho = 0.76-0.88, p &lt; 0.05) was observed. All exome capture-based methods reliably detected outlier expression of actionable genes, including ERBB2, MET, NTRK1, and PPARG. In urothelial cancer samples, the Agilent assay was associated with the highest molecular subtype concordance with FF-TruSeq (Cohen's k = 0.7, p &lt; 0.01). Both Agilent and IDT detected all the clinically relevant fusions that were initially identified in FF-TruSeq. Conclusion: All exome capture-based methods had comparable performance and concordance with FF-TruSeq. By enabling the interrogation of FFPE tumor samples, our findings will enable the implementation of RNA-seq in the clinic to guide precision oncology approaches.

  Study EGAS00001005255

• Are children born after medical assisted reproduction at greater risk of having an increased de novo mutation rate?

  Introduction: De novo mutations (DNMs) play a prominent role in sporadic disorders with reduced fitness such as infertility and intellectual disability. Advanced paternal age is known to increase disease risk in offspring by increasing the number of DNMs in their genome. Less is known about the effect of assisted reproduction techniques (ART) on the number of DNMs in offspring. With the on-going trend of delayed parenthood more children are now born both from older fathers and through ART. Materials and Methods: We investigated 49 trios (mother, father and child) and 2 quartets (mother, father and 2 siblings) divided into children born after spontaneous conception (n=18); born after in vitro fertilisation (IVF) (n=17) and born after intracytoplasmic sperm injection combined with testicular sperm extraction (ICSI-TESE) (n=18). Groups further divided by paternal age, young (&lt;35) or old (&gt;45 years of age at conception). Whole-genome sequencing was performed twice to independently detect and validate all DNMs in children. Results: A clear paternal age effect was observed, with 70 DNMs detected on average in children born to young fathers and 94 DNMs in those born to older fathers (p = 0.001). No significant differences were observed between different methods of conception (p = 1) with paternal age affecting all methods equally. Conclusions: Paternal age, not method of conception, had a major effect on the observed number of DNMs in offspring. Given the role DNMs in disease risk, this negative result is good news for IVF and ICSI-TESE born children, if replicated in larger cohorts.

  Study EGAS00001005569

• Stromal-induced epithelial-mesenchymal transition induces drug resistance in acute lymphoblastic leukemia

  Bone marrow (BM) microenvironment-induced drug resistance in acute lymphoblastic leukemia (ALL) is thought to be mediated by intercellular interaction of leukemic cells with BM niches, but the mechanistic basis of resistance is poorly understood. Here, we genomic and functional analyses, we show that ALL-BM mesenchymal stromal cells (MSC) interactions encompass direct cell-cell adhesion via integrin α4β1, uni-directional secretion of soluble factors from MSC to ALL, and bidirectional transfer of cellular vesicles through tunneling nanotubes. We show that this crosstalk drives an epithelial-mesenchymal transition (EMT)-like program in ALL cells adhering to MSC, thereby triggering stem cell-like features including drug resistance and survival advantage. Moreover, single-cell RNA sequencing of cell line and primary ALL samples identified a hybrid cell state derived from B-ALL adhered to MSC, that exhibits both B-ALL and MSC gene signatures that are orchestrated by an EMT-like program. Multiple EMT-activating pathways including WNT/β-catenin, cytokine/STAT3, TGFβ, NOTCH1, and NF-κB-dependent signaling pathways are activated in MSC-adherent ALL cells; however, only the inhibition of either WNT/β-catenin or cytokine/STAT3 pathway could attenuate the survival and drug resistance of ALL cells adherent to MSC. Notably, we identified that this hybrid cell state exhibited elevated WNT/β-catenin-mediated transcriptional activity. Finally, we show that blocking of the interaction of β-catenin and CBP interaction abrogates the survival and drug resistance of ALL cells adhered to MSC, indicating that targeting of an EMT-like program is potential therapeutic regimen in targeting cell extrinsically acquired drug resistance in ALL.

  Study EGAS00001006120

• The Genomic Landscape of Response to EGFR Blockade in Colorectal Cancer

  Colorectal cancer (CRC) is the third most common cancer world-wide with 1.2 million patients diagnosed yearly. In late stage CRC, the most commonly used targeted therapies are monoclonal antibodies cetuximab and panitumumab, which inactivate EGFR. Recent studies have identified alterations in KRAS and other genes as likely mechanisms of primary and secondary resistance to anti-EGFR antibody therapy. Despite these efforts, additional mechanisms of resistance to EGFR blockade are thought to be present in CRC and little is known about determinants of sensitivity to this therapy. To examine the effect of somatic genetic changes in CRC on response to anti-EGFR antibody therapy, we performed complete exome sequence and copy number analyses of 129 patient-derived tumorgrafts and targeted genomic analyses of 55 patient tumors, all of which were KRAS wild-type. We analyzed the response of tumors to anti-EGFR antibody blockade in tumorgraft models or in clinical settings. In addition to previously identified genes, we detected mutations in ERBB2, EGFR, FGFR1, PDGFRA, and MAP2K1 as potential mechanisms of primary resistance to this therapy. Novel alterations in the ectodomain of EGFR were identified in patients with acquired resistance to EGFR blockade. Amplifications and sequence changes in the tyrosine kinase receptor adaptor gene IRS2 were identified in tumors with increased sensitivity to anti-EGFR therapy. Therapeutic resistance to EGFR blockade could be overcome in tumorgraft models through combinatorial therapies targeting actionable genes. These analyses provide a systematic approach to evaluate response to targeted therapies in human cancer, highlight new mechanisms of responsiveness to anti-EGFR therapies, and provide new avenues for intervention in the management of CRC.

  Study EGAS00001001305

• An integrative model of pathway convergence in genetically heterogeneous blast crisis chronic myeloid leukemia (CML)

  Targeted therapies against the BCR-ABL1 kinase have revolutionized treatment of chronic phase (CP) chronic myeloid leukemia (CML). In contrast, management of blast crisis (BC) CML remains challenging because BC cells acquire complex molecular alterations that confer stemness features to progenitor populations and resistance to BCR-ABL1 tyrosine kinase inhibitors. Comprehensive models of BC transformation have proved elusive because of the rarity and genetic heterogeneity of BC, but are important for developing biomarkers predicting BC progression and effective therapies. To better understand BC, we performed an integrated multiomics analysis of 74 CP and BC samples using whole-genome and exome sequencing, transcriptome and methylome profiling, and chromatin immunoprecipitation followed by high-throughput sequencing. Employing pathway-based analysis, we found the BC genome was significantly enriched for mutations affecting components of the polycomb repressive complex (PRC) pathway. While transcriptomically, BC progenitors were enriched and depleted for PRC1- and PRC2-related gene sets respectively. By integrating our data sets, we determined that BC progenitors undergo PRC-driven epigenetic reprogramming toward a convergent transcriptomic state. Specifically, PRC2 directs BC DNA hypermethylation, which in turn silences key genes involved in myeloid differentiation and tumor suppressor function via so-called epigenetic switching, whereas PRC1 represses an overlapping and distinct set of genes, including novel BC tumor suppressors. On the basis of these observations, we developed an integrated model of BC that facilitated the identification of combinatorial therapies capable of reversing BC reprogramming (decitabine+PRC1 inhibitors), novel PRC-silenced tumor suppressor genes (NR4A2), and gene expression signatures predictive of disease progression and drug resistance in CP.

  Study EGAS00001001751

• Genomic and epigenomic characterization of juvenile myelomonocytic leukemia (JMML)

  Juvenile myelomonocytic leukemia (JMML) is an aggressive myeloproliferative disorder of early childhood. While some cases show spontaneous remission, allogeneic hematopoietic stem cell transplantation (HSCT) remains the only curative treatment option for the majority of patients, however, the 5-year event-free survival reaches only about 50%. Hyperactive RAS signaling is assumed to be the main driving event in JMML. It is caused by genetic alterations in CBL, KRAS, NF1, NRAS, or PTPN11 in about 90% of patients. So far, there is no clear understanding of how RAS pathway mutations relate to the heterogeneous disease biology and variable clinical outcome seen in JMML patients. As a consequence, established clinical and genetic markers fail to fully represent the observed disease heterogeneity. We hypothesized that DNA methylation profiling, either alone or in combination with genetic alterations, might provide a molecular basis for disease classification. Genome wide DNA methylation analysis using the HumanMethylation450 Bead Chip array was performed in a discovery cohort of 20 patients (19 JMML, 1 Noonan syndrome). The clinical relevance of our findings was assessed in an unselected sample set consisting of 147 consecutive patients with JMML (n=129) or Noonan syndrome (n=18) associated myeloproliferative disorder (n=18) registered in the EWOG-MDS 1998 &amp; 2006 trials. Data integration was performed in a subset of patients with available exome sequencing (n=50) and expression profiling (n=15) data. This integrated approach identified three JMML subgroups characterized by distinct clinical and biological features. This analysis further suggested a molecular mechanism by which additional genetic events, presumably further activating the RAS-RAF-MEK-ERK pathway, mediate DNA hypermethylation via up-regulation of DNMTs in more aggressive JMML cases.

  Study EGAS00001002511

• MiR expression profiles of paired primary colorectal cancerand metastases by next-generation sequencing

  MicroRNAs (miRs) have been recognized as promising biomarkers. It is unknown to what extent tumor-derived miRs aredifferentially expressed between primary colorectal cancers (pCRCs) and metastatic lesions, and to what extent the expressionprofiles of tumor tissue differ from the surrounding normal tissue. Next-generation sequencing (NGS) of 220 fresh-frozen samples,including paired primary and metastatic tumor tissue and non-tumorous tissue from 38 patients, revealed expression of 2245known unique mature miRs and 515 novel candidate miRs. Unsupervised clustering of miR expression profiles of pCRC tissue withpaired metastases did not separate the two entities, whereas unsupervised clustering of miR expression profiles of pCRC withnormal colorectal mucosa demonstrated complete separation of the tumor samples from their paired normal mucosa. Twohundred and twenty-two miRs differentiated both pCRC and metastases from normal tissue samples (false discovery rate (FDR)o0.05). The highest expressed tumor-specific miRs were miR-21 and miR-92a, both previously described to be involved in CRC withpotential as circulating biomarker for early detection. Only eight miRs, 0.5% of the analysed miR transcriptome, were differentiallyexpressed between pCRC and the corresponding metastases (FDR o0.1), consisting of five known miRs (miR-320b, miR-320d,miR-3117, miR-1246 and miR-663b) and three novel candidate miRs (chr 1-2552-5p, chr 8-20656-5p and chr 10-25333-3p). Theseresults indicate that previously unrecognized candidate miRs expressed in advanced CRC were identified using NGS. In addition,miR expression profiles of pCRC and metastatic lesions are highly comparable and may be of similar predictive value for prognosisor response to treatment in patients with advanced CRC

  Study EGAS00001001127

• Molecular_diagnosis_of_albinism

  Albinism is genetically heterogeneous rare genetic condition affecting 1:17000 in the Western world (but more frequent in Africa) whose main feature is a profound visual impairment, characterised by foveal hypoplasia, abnormal chiasmatic connections, nystagmus and photofobia. All these features result in severly altered visual acuity (&lt;0,1), absent depth perception and poor night vision. People with albinism are primarily visually handicapped. In addition, for some types of albinism, the visual phenotype can be presented with partial or total hypopigmentation, hence resulting in a secondary phenotype which can lead to skin cancer if skin is not adequately protected. Recently a new syndrome has been described, FHONDA, with the same visual abnormalities of albinism but without pigment alteration. The traditional classification differentiates Oculoculatenous albinism (OCA), where hypopigmentation involves hair, skin and eyes versus Ocular Albinism (OA), where hypopigmentation only affects the eyes. These are non-sydrimic types of albinism. Some syndromic forms (Hermansky-Pudlak=HPS, Chediak-Higashi=CHS) affect cells beyond pigment cells, present in the lungs, immune system, platelets and intestines, resulting in more severe phenotypes that can be fatal. Mutations in at least 19 genes are assocaited with the corresponding types of albinism. Most hospitals will only diagnose the most frequent cases using traditional Sanger, MLPA approaches. Some will use CGH arrays. We aim to diagnose all cases of albinism through the Albinochip proposal, which combines a Sequenom first step of known mutations combined with subsequent NGS approaches. In some cases we fail to find a second mutation, these are good candidates for further full exome analyses. 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 EGAS00001002068

• PCGP Ph-like ALL

  Background: BCR-ABL1-like, or Ph-like acute lymphoblastic leukemia is characterized by a gene expression profile similar to BCR-ABL1 positive ALL, genetic alterations of lymphoid transcription factor genes, and poor outcome. Sequencing of small numbers of Ph-like ALL cases has identified genetic alterations activating kinase signaling suggesting Ph-like ALL may be amenable to treatment with tyrosine kinase inhibitors. However, the spectrum of genetic alterations in childhood and adult Ph-like ALL is incompletely understood. Methods: We performed genomic profiling of 1736 B-ALL cases and next-generation sequencing for 160 Ph-like cases. We examined the functional effects of chimeric fusion proteins in mouse cell lines. Results: The frequency of Ph-like ALL rose from 11% in children to 26% in young adults, and was associated with very poor outcome. Kinase-activating alterations were identified in 90% of Ph-like cases, most commonly fusions involving 10 kinase or cytokine receptor genes (ABL1, ABL2, CRLF2, CSF1R, EPOR, JAK2, NTRK3, PDGFRB, PTK2B and TYK2), and mutations involving FLT3, IL7R and SH2B3. Expression of ABL1, ABL2, CSF1R and JAK2 fusions resulted in cytokine-independent proliferation of cell lines and activation of pSTAT5. Cells expressing ABL1, ABL2, CSF1R and PDGFRB fusions were sensitive to dasatinib, and JAK2 fusions to ruxolitinib. Conclusions: Ph-like ALL is characterized by a diverse range of genomic alterations that converge on a limited number of kinase signaling pathways amenable to inhibition with currently available tyrosine kinase inhibitors. Trials identifying Ph-like ALL and testing the efficacy of tyrosine kinase inhibitor therapy are warranted

  Study EGAS00001000654

• Tracking the evolution of Therapy-Related Myeloid Neoplasms using chemotherapy signatures

  Patients treated with cytotoxic therapies, including autologous stem cell transplantation, are at risk for developing therapy-related myeloid neoplasms (tMN). Pre-leukemic clones (i.e., clonal hematopoiesis; CH) are detectable years before the development of these aggressive malignancies, though the genomic events leading to transformation and expansion are not well-defined. Here, leveraging distinctive chemotherapy-associated mutational signatures from whole-genome sequencing data and targeted sequencing of pre-chemotherapy samples, we reconstruct the evolutionary life-history of 39 therapy-related myeloid malignancies. A dichotomy is revealed, in which neoplasms with evidence of chemotherapy-induced mutagenesis from platinum and melphalan are hypermutated and enriched for complex structural variants (i.e., chromothripsis) while neoplasms with non-mutagenic chemotherapy exposures are genomically similar to de novo acute myeloid leukemia. Using chemotherapy-associated mutational signatures as temporal barcodes linked to a discrete clinical exposure in each patient’s life, we estimate that several complex events and genomic drivers are acquired after chemotherapy is administered. For patients with prior multiple myeloma who were treated with high-dose melphalan and autologous stem cell transplantation, we demonstrate that tMN can develop from either a reinfused CH clone that escapes chemotherapy exposure and is selected following reinfusion, or from TP53-mutant CH that survives direct myeloablative conditioning and acquires melphalan-induced DNA-damage. Overall, we reveal a novel mode of tMN progression that is not reliant on direct mutagenesis or even exposure to chemotherapy. Conversely, for tMN that evolve under the influence of chemotherapy-induced mutagenesis, distinct chemotherapies not only select pre-existing CH, but also promote the acquisition of recurrent genomic drivers.

  Study EGAS00001006903

• 'Targeted High Throughput Sequencing in Clinical Cancer Settings: Formaldehyde fixed-paraffin embedded (FFPE) tumor tissues, input amount and tumor heterogeneity'

  "Background: Massively parallel sequencing technologies have brought an enormous increase in sequencing throughput. However, these technologies need to be further improved with regard to reproducibility and applicability to clinical samples and settings. Methods: Using identification of genetic variations in prostate cancer as an example we address three crucial challenges in the field of targeted re-sequencing: Small nucleotide variation (SNV) detection in samples of formalin-fixed paraffin embedded (FFPE) tissue material, minimal amount of input sample and sampling in view of tissue heterogeneity.Results: We show that FFPE tissue material can supplement for fresh frozen tissues for the detection of SNVs and that solution-based enrichment experiments can be accomplished with small amounts of DNA with only minimal effects on enrichment uniformity and data variance. Finally, we address the question whether the heterogeneity of a tumor is reflected by different genetic alterations, e.g. different foci of a tumor display different genomic patterns. We show that the tumor heterogeneity plays an important role for the detection of copy number variations, but is of minor importance for the detection of somatic variations.Conclusions: The application of high throughput sequencing technologies in cancer genomics opens up a new dimension for the identification of disease mechanisms. In particular the ability to use small amounts of FFPE samples available from surgical tumor resections and histopathological examinations facilitates the collection of precious tissue materials. However, care needs to be taken in regard to the locations of the biopsies taken, which can have an influence on the prediction of copy number variations. Bearing these technological challenges in mind will significantly improve many large-scale sequencing studies and will - on a long-range - result in a more reliable prediction of individual cancer therapies."

  Study EGAS00001000136

• Deep single-cell RNA sequencing data for 12346 T cells from tumour, adjacent normal tissue and peripheral blood of treatment-naive NSCLC patients

  Cancer immunotherapies have shown sustained clinical responses in treating non-small cell lung cancer (NSCLC), but efficacy varies between patients and is believed to depend in part on the amount and properties of tumor infiltrating lymphocytes (TILs). To comprehensively depict and dissect the baseline landscape of the composition, lineage and functional states of TILs in lung cancer, here we generated deep single-cell RNA sequencing data for 12,346 T cells from the primary tumour, adjacent normal tissues and peripheral blood of 14 treatment-naive NSCLC patients. Combined expression and TCR-based lineage tracking revealed a significant proportion of effector T cells with common origins and similar functional states across peripheral blood and tumours pointing towards a highly migratory nature of these T cells. We also observed tumour-infiltrating CD8+ T cells undergoing extensive clonal expansion and exhaustion, with two clusters of cells exhibiting states preceding exhaustion. Survival analysis on independent datasets suggested that high ratio of pre-exhausted to exhausted T cells is associated with better prognosis of lung adenocarcinoma (LUAD). In addition, we observed further heterogeneity within the tumour regulatory T cells (Tregs), characterized by the bimodal distribution of TNFRSF9, an activation marker for antigen-specific Tregs. The gene signature of this group of activated tumour Tregs, which included IL1R2, correlated with poor prognosis in LUAD. The T cell clusters revealed by our single cell analyses provide a new approach for patient stratification, and the accompanying compendium of data will help the research community to gain further insight into the functional states and dynamics of T cell responses in lung cancer.

  Study EGAS00001002430

• Ultra-Fast Patient-Derived Xenografts Identify Functional and Spatial Tumour Heterogeneities that Drive Therapeutic Resistance

  Treatment of patients with advanced cancers increasingly relies on agents targeting specific molecular or cellular aberrations. These expensive therapies often fail because of pre-existing drug resistance. Intra-tumoural molecular and micro-environmental heterogeneity resistance can create resistance in one or more subclones of a tumour, leading to spatio-functional heterogeneity. We created a rapid approach for quantifying tumour functional phenotypes based on tumour-implantation into the chorioallantoic membrane of chick embryos. We apply this technique to 43 separate biopsies from 6 patients with renal cell carcinoma (RCC) with 36 replicates per biopsy, achieving a 93.6% engraftment rate (1449/1548). Challenging these models with sunitinib, an anti-angiogenic therapy frequently used in RCC allowed us to quantify spatio-functional heterogeneity in individual patients. Prediction of de novo drug resistance in advanced cancer patients would facilitate precision medicine by tailoring systemic therapies and improve clinical trial design outcome measures. We describe a patient-derived xenograft platform where fresh patient tumour samples are implanted into the chorioallantoic membrane (CAM) of chick embryos at high engraftment efficiencies, permitting large scale “tumour avatar” studies to guide the selection of drugs and to predict clinical drug resistance outcomes in renal cell carcinoma (RCC) patients. Applying this “tumour avatar” model in a co-clinical trial on RCC patients, we observe intra-tumoural functional heterogeneity with sunitinib treatment, revealing the presence of drug-resistant clones prior to the initiation of therapy. Multi-regional exome Matched DNA sequencing of the primary tumour and metastases revealed DNA mutational signatures associated with sunitinib resistance at baseline and after treatment. These These studies demonstrate a rapid and affordable method of findings confirm the existence of genetic and surveying functional tumour heterogeneities heterogeneity.

  Study EGAS00001002627

• Distinct evolutionary trajectories of primary high grade serous ovarian cancers revealed through spatial mutational profiling

  High-grade serous ovarian cancer (HGSC) is characterized by poor outcome, often attributed to emergence of treatment-resistant sub-clones. We sought to measure the degree of genomic diversity within primary, untreated HGSC to examine the natural state of tumor evolution prior to therapy. We performed exome sequencing, copy number analysis, targeted amplicon deep sequencing and gene expression profiling on thirty-one spatially and temporally separated HGSC tumor specimens (six patients) including ovarian masses, distant metastases, and fallopian tube lesions. We found widespread intra-tumoral variation in mutation, copy number, and gene expression profiles, with key driver alterations in genes present in only a subset of samples (e.g. PIK3CA, CTNNB1, NF1). On average, only 51.5% of mutations were present in every sample of a given case (range: 10.2% to 91.4%), with TP53 as the only somatic mutation consistently present in all samples. Complex segmental aneuploidies, such as whole genome doubling, were present in a subset of samples from the same individual, with divergent copy number changes segregating independently of point mutation acquisition. Reconstruction of evolutionary histories showed one patient with mixed HGSC and endometrioid histology with common etiologic origin in the fallopian tube and subsequent selection of different driver mutations in the histologically distinct samples. In this patient, we observed mixed cell populations in the early fallopian tube lesion, indicating diversity arises at early stages of tumorigenesis. Our results reveal that HGSC exhibit highly individual evolutionary trajectories and diverse genomic tapestries prior to therapy, exposing an essential biological characteristic to inform future design of personalized therapeutic solutions and investigation of drug resistance mechanisms.

  Study EGAS00001000547

• UK10K NEURO FSZ

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. These Finnish schizophrenia samples have been collected from a population cohort using national registers. The entire sample collection consists of 2756 individuals from 458 families of whom 931 are diagnosed with schizophrenia spectrum disorder, each family having at least two affected siblings. 170 families originate from an internal isolate (Kuusamo) with a three-fold lifetime risk for the trait. The genealogy of the internal isolate is well documented and the individuals form a “megapedigree” reaching to the 17th Century. All diagnoses are based on DSM-IV and for a large fraction of cases there is cognitive data. For further details/descriptions with regard to this data set please contact Tiina Paunio (tiina.paunio@thl.fi)

  Study EGAS00001000118

• A genome scan for genes underlying adult body size differences between Central African hunter-gatherers and farmers.

  The evolutionary and biological bases of the Central African pygmy phenotype, a characteristic of rainforest huntergatherers defined by reduced body size compared with neighboring farmers, remain largely unknown. Here, we perform a joint investigation in Central African hunter-gatherers and farmers of adult standing height, sitting height, leg length, and body mass index (BMI), considering 358 hunter-gatherers and 169 farmers with genotypes for 153,798 SNPs. In addition to reduced standing heights, hunter-gatherers have shorter sitting heights and leg lengths and higher sitting/standing height ratios than farmers and lower BMI for males. Standing height, sitting height, and leg length are strongly correlated with inferred levels of farmer genetic ancestry, whereas BMI is only weakly correlated, perhaps reflecting greater contributions of non-genetic factors to body weight than to height. Single- and multi-marker association tests identify one region and eight genes associated with hunter-gatherer/farmer status, and 24 genes associated with the height-related traits. Many of these genes have putative functions consistent with roles in determining their associated traits and the pygmy phenotype, and they include three associated with standing height in non-Africans (PRKG1, DSCAM, MAGI2). We find evidence that European height-associated SNPs or variants in linkage disequilibrium with them contribute to standing- and sitting-height determination in Central Africans, but not to the differential status of hunter-gatherers and farmers. These findings provide new insights into the biological basis of the pygmy phenotype, and they highlight the potential of cross-population studies for exploring the genetic basis of phenotypes that vary naturally across populations. Pemberton et al. Human Genetics 2018 https://doi.org/10.1007/s00439-018-1902-3

  Study EGAS00001002975

• Indonesian Genome Diversity Project

  Indonesia, the world’s fourth largest country by population, comprises an archipelago of about 900 permanently inhabited islands in tropical Asia, and hosts an astonishing array of human diversity that remains largely underrepresented in modern biological surveys. The region is unique for its key role in both the early and recent evolution of Homo outside of Africa. It has one of the first traces of anatomically modern humans in Eurasia, and was the focus of extensive interaction between archaic hominins and humans – both in the likely co-existence of H. floresiensis with modern humans, and in the introgression of Denisovans into the ancestors of its present inhabitants. More recently, Indonesia was a center of the spread of Neolithic culture by Austronesian speakers. Advancing maritime technologies allowed farming populations to treat this region as a springboard to reach oceanic islands as remote as Madagascar, Hawaii, Easter Island and New Zealand. Here we present complete human genomes from healthy individuals from multiple key islands in Indonesia and Papua. The genomes were sequenced to high depth (30x). The purpose of this dataset is to address different questions of human evolution in Island Southeast Asia and beyond, including early settlement of the area by anatomically modern humans, population structure, genetic admixture and adaptation, and contacts with archaic hominins. We employed a sampling strategy to capture a broad range of diversity across the region incorporating an average of 12 samples from each of 14 islands spanning Indonesia and Papua, from Sumatra in the west to New Britain in the east. Genomes from this dataset form part of an ongoing project to describe and understand human evolutionary history in Island Southeast Asia and Eurasia in general.

  Study EGAS00001003054

• Integrative genomic profiling of hepatocellular adenomas reveals recurrent FRK activating mutations and mutational processes of malignant transformation

  Hepatocellular adenomas (HCA) are benign liver tumors predominantly developed in women using oral contraception. Here, exome sequencing identified novel recurrent somatic FRK mutations that induce constitutive kinase activity, STAT3 activation and cell proliferation sensitive to Src inhibitors. We also found new recurrent mutations activating JAK1, gp130 or ß-catenin. Chromosome copy-number and methylation profiling revealed patterns that correlated with specific gene mutations and tumor phenotypes. Finally, integrative analysis of HCA transformed to hepatocellular carcinoma revealed ß-catenin mutation as an early alteration and TERT promoter mutations as associated with the last step of the adenoma-carcinoma transition. In conclusion, we identified the genomic diversity in benign hepatocyte proliferation, new therapeutic targets and the key genomic determinants of the adenoma-carcinoma transformation sequence. Significance: Malignant transformation of benign adenoma into malignant carcinoma is frequently observed in several epithelial tumor types. Adenoma-carcinoma transition is crucial for the patient prognosis, however, little is known on the molecular mechanisms involved. Here, we performed an integrated genomic analysis of hepatocellular adenomas (HCA). Among several novel gene mutations, we identified recurrent somatic mutation activating FRK, a Src-like kinase. FRK activating mutations induce STAT3 activation and cell proliferation that can be targeted by Src inhibitors. Focusing on malignant transformation of HCA in hepatocellular carcinoma, we identified CTNNB1 and TERT promoter mutation as early and late genomic events involved in adenoma-carcinoma transition. In conclusion, we propose new therapeutic targets and an identification of adenomas with the highest risk of malignant transformation. Highlights: - FRK is activated by somatic recurrent mutations. - FRK mutants promote cell proliferation sensitive to src inhibitors. - CTNNB1 activation is an early event in HCA transformation. - TERT promoter mutations are associated with adenoma-carcinoma transition.

  Study EGAS00001000679

• Dynamics of genomic clones in breast cancer patient xenografts at single cell resolution

  Human cancers, including breast cancers, are comprised of clones differing in mutation content that evolve dynamically in space and time following principles of Darwinian evolution, underpinning important emergent features such as drug resistance and metastasis. Human breast cancer xenoengraftment is used as a means of capturing and studying tumour biology, and breast tumour xenografts are generally assumed to be reasonable models of the originating tumours. However the consequences and reproducibility of engraftment and propagation on the genomic clonal architecture of tumours has not been systematically examined at single cell resolution. Here we show by both deep genome and single cell sequencing methods, the clonal dynamics of initial engraftment and subsequent serial propagation of primary and metastatic human breast cancers in immunodeficient mice. In all cases examined, clonal selection on engraftment was observed in both primary and metastatic breast tumours, varying in degree from extreme selective engraftment of minor (&lt;5% of starting population) clones to moderate, polyclonal engraftment. Furthermore, ongoing clonal dynamics during serial passaging is a feature of tumours experiencing modest initial selection. Through single cell sequencing, we show that major mutation clusters estimated from tumour population sequencing relate predictably to the most abundant clonal genotypes, even in clonally complex and rapidly evolving cases. Finally, we show that similar clonal expansion patterns can emerge in independent grafts of the same starting tumour population, indicating that genomic aberrations can be reproducible determinants of evolutionary trajectories. Our results show that measurement of genomically defined clonal population dynamics will be highly informative for functional studies utilizing patient-derived breast cancer xenoengraftment.

  Study EGAS00001000952

• 1 Intratumoral genetic heterogeneity and clonal evolution following neoadjuvant chemoradiotherapy (nCRT) in locally advanced rectal tumors.

  Purpose: Neoadjuvant chemoradiotherapy (nCRT) is the standard treatment forpatients with locally advanced rectal cancer. However, little is known about the effect of nCRT in the mutational landscape and intratumoral genetic heterogeneity (ITGH) of rectal tumors. Moreover, response to nCRT vary substantially among patients and the genetic determinants of treatment response are not fully understood. Here, we sought to determine the impact of nCRT in the mutational landscape and ITGH of rectal tumors and to identify putative drivers of therapeutic resistance. Patients and Methods: We analyzed whole-exome sequencing (WES) and clinical data from 79 primary non-treated rectal cancers obtained from The Cancer Genome Atlas (TCGA). We also compared the mutational landscape of a matched set of peripheral blood cells, irradiated tumor-adjacent colonic mucosa and pre and post-treatment tumor samples to determine the iatrogenic effects of nCRT and its impact in tissue genetic heterogeneity. Finally, we performed WES of 7 matched pre and post-treatment tumor samples to examine how nCRT affects the clonal architecture of rectal tumors. Results: We show that rectal tumors exhibit remarkable ITGH, which increases with disease progression. We also show that nCRT, per se, does not introduce novel somatic mutations, nor alters the clonal architecture of the irradiated colonic mucosa. Instead, we show that it acts as a potent selective pressure, favoring the expansion of tumor cell subpopulations more prone to resist therapy and increasing ITGH in the residual tumor. Conclusion: Our results reveal the highly heterogeneous and dynamic clonal architecture of rectal tumors undergoing nCRT. Our results also uncover putative genetic determinants of response to nCRT and a possible role for ITGH in predicting treatment outcome.

  Study EGAS00001003250

• Genetic characterization of B-cell prolymphocytic leukemia: a hierarchical prognostic model involving MYC and TP53 abnormalities - RNA-seq

  B-cell prolymphocytic leukemia (B-PLL) is a rare disease, whose molecular pathology is largely unknown. We report here the cytogenetic and molecular findings in a large series of 34 B-PLL. Karyotype (K) was complex (≥3 abnormalities) in 73%, and highly complex (HCK≥5) in 45%. The most frequent chromosomal aberrations were: translocation targeting the MYC gene [t(MYC)] (62%), 17p deletion including TP53 gene (38%), trisomy 18/18q (30%), 13q14 deletion (29%), trisomy 3 (24%), trisomy 12 (24%) and 8p deletion (23%). Whole-Exome Sequencing performed in 16 patients revealed recurrent mutations in TP53 (6/16, 38%), MYD88 (n=4), BCOR (n=4), MYC (n=3), SF3B1 (n=3), FAT1 (n=3), SETD2 (n=2), CHD2 (n=2), CXCR4 (n=2) and BCLAF1 (n=2). The main group of patients (21/34, 62%) had a t(MYC) associated with a higher percentage of prolymphocytes (p=0.03), CD38 expression (p&lt;0.001), lower K complexity (p=0.0004), mutations in MYC and in genes involved in RNA metabolism and chromatin remodeling. Principal component analysis of gene expression data showed that patients with t(MYC) clustered together. A second group with MYC gain (5/34, 15%), was associated with HCK≥5 (p=0,01) and trisomy 3 (p=0,008). Altogether, 26/34 patients (76%) had a MYC activation, translocation or gain, that were mutually exclusive. We identified 3 distinct cytogenetic prognostic groups (p=0.0006): lower risk: absence of MYC activation (median not reached); intermediate risk: MYC activation without del17p (125 months); high risk: MYC activation + del17p (11 months). Our results show that cytogenetic analysis is a useful diagnostic tool in B-PLL that improves prognostic stratification.

  Study EGAS00001003274

• Systems biology of Colorectal Cancer

  Colorectal cancer (CRC) is one of the most common cancers in both males and females, and it is perhaps the best understood of all epithelial tumors in terms of its molecular origin. Yet, despite large amount of work that has concentrated on understanding of colon tumorigenesis, we still do not know the full complement of molecular lesions that are individually necessary – and together sufficient – to cause colorectal cancer. Neither do we understand why some specific mutations that are relatively rare in other tumors (e.g. loss of the APC tumor suppressor) are extremely common in colorectal cancer. We propose here to use the tools of systems biology to develop a quantitative and comprehensive model of colorectal tumorigenesis. The model will include identification of cell-type specific and oncogenic pathways that contribute to colon tumorigenesis, and explain in molecular detail how a genotype of an individual CRC leads to activation of downstream genes that drive uncontrolled cell growth. This model will subsequently be used to find novel therapeutic targets, to guide genetic screening to identify individuals with elevated risk for developing CRC, and to classify patients into molecular subgroups to select the treatment combination that is optimal for each patient (personalized medicine). The specific objectives of the SYSCOL project are: 1. Identify genetic markers for individual risk using genotyping and sequencing of constitutional DNA from sporadic and familial CRC cases and controls 2. Identify genes and regulatory elements that contribute to colorectal cancer cell growth 3. Use data from Aims 1-2 to develop a quantitative model for colorectal tumorigenesis 4. Apply the model for identification of high-risk individuals, for molecular classification of the disease, and for identification of novel molecular treatment targets

  Study EGAS00001000854

• Mutational patterns and regulatory networks in epigenetic subgroups of meningioma (H033)

  DNA methylation patterns delineate clinically relevant subgroups of meningioma. We previously established the six meningioma methylation classes (MC) benign-1, 2, 3, intermediate-A, B and malignant. Here, we set out to identify subgroup-specific mutational patterns and pathway regulation. Whole-genome-sequencing was performed on 62 samples across all MCs and WHO grades from 62 patients with matched blood control, including 40 sporadic and 22 radiation-induced (Mrad) meningiomas. RNA sequencing was added for 18 cases and chromatin-immunoprecipitation for the enhancer mark H3K27ac followed by sequencing (ChIP-seq) for 16 samples. Besides the known mutations in meningioma, structural alterations were found to contribute to the spectrum of mechanisms inactivating NF2 in sporadic meningioma similar to previous reports for Mrad. Aberrations of DMD were found to be enriched in MCs with NF2 mutations and DMD was among the most differentially upregulated genes in NF2 mutant compared to NF2 wild-type cases. The mutational signature AC3 was detected in both sporadic meningioma and Mrad, but distributed across the genome in sporadic cases and enriched near genomic breakpoints in Mrad. In general, the malignant MC presented a significantly higher exposure to AC3 and higher genomic instability beyond the mutational load than the other MCs. Pathway analysis of the ChIP-seq clusters revealed that the FOXM1 is most differentially activated in high-grade cases, along with super enhancer recruitment near HOX genes and respective upregulation of expression. This data further elucidate the biological mechanisms differentiating meningiomas of different WHO grade and epigenetic subgroups and suggest leveraging the genomic instability as novel therapeutic targets.

  Study EGAS00001003481

• Telomerase activation by genomic rearrangements in high-risk neuroblastoma

  Neuroblastoma is a malignant pediatric tumor of the sympathetic nervous system1. Roughly half of these tumors regress spontaneously or are cured by limited therapy. By contrast, high-risk neuroblastomas have an unfavorable clinical course despite intensive multimodal treatment, and their molecular basis has remained largely elusive2-4. We have performed whole-genome sequencing of 56 neuroblastomas (high-risk, n=39; low-risk, n=17) and discovered recurrent genomic rearrangements affecting a chromosomal region (5p15.22) proximal of the telomerase reverse transcriptase gene (TERT). These rearrangements occurred only in high-risk neuroblastomas (12/39, 31%) in mutually exclusive fashion with MYCN amplifications and ATRX mutations, which are known genetic events in this tumor type1,2,5. In an extended case series (n=217), TERT rearrangements defined a subgroup of high-risk tumors with particularly poor outcome. Despite the large diversity of these rearrangements, they all induced massive transcriptional upregulation of TERT. In the remaining high-risk tumors, TERT expression was also elevated in MYCN-amplified tumors, whereas alternative lengthening of telomeres was present in neuroblastomas without TERT or MYCN alterations, suggesting that telomere lengthening represents a central mechanism defining this subtype. The 5p15.22 rearrangements juxtapose the TERT coding sequence to strong enhancer elements, resulting in massive chromatin remodeling and DNA methylation of the affected region. Supporting a functional role of TERT, neuroblastoma cells bearing rearrangements or amplified MYCN exhibited both upregulated TERT expression and enzymatic telomerase activity. In summary, our findings show that remodeling of the genomic context abrogates transcriptional silencing of TERT in high-risk neuroblastoma and places telomerase activation in the center of transformation in a large fraction of these tumors.

  Study EGAS00001001308

• Evolutionary dynamics of residual disease in human glioblastoma

  Glioblastoma is the most common and aggressive adult brain malignancy against which conventional surgery and chemoradiation provide limited benefit. Even when a good treatment response is obtained, recurrence inevitably occurs either locally (c.80%) or distally (c.20%), driven by cancer clones that are often genomically distinct from those in the primary tumour. Glioblastoma cells display a characteristic infiltrative phenotype, invading the surrounding tissue and often spreading across the whole brain. We hypothesised that cancer cells responsible for relapse reside in two compartments of residual disease that are left behind after treatment: the infiltrated normal brain parenchyma, and the sub-ventricular zone (SVZ). In this model, residual disease subclones diverge early during glioblastoma evolution, may remain dormant in the normal parenchyma and are responsible for recurrence. To test our hypothesis, we performed whole-exome sequencing of 69 multi- region samples collected using fluorescence-guided resection from 11 patients, including the infiltrating tumour margin (M) and the SVZ for each patient, as well as matched blood. We used a phylogenomic approach to dissect the spatio-temporal evolution of each tumour and unveil the relation between residual disease and the main tumour mass. We also analysed two patients with paired primary-recurrence samples with matched residual disease. Our results suggest that the infiltrative phenotype is an early evolutionary trait of glioblastoma, giving rise to the tumour mass detected at surgery. After treatment, the same infiltrative clone may seed the growth of a recurrent tumour, thus representing the ‘missing link’ between the primary tumour and recurrent disease. These results are consistent with recognised clinical phenotypic behaviour and suggest that more specific therapeutic targeting of cells in the infiltrated brain parenchyma may improve patient’s outcome.

  Study EGAS00001003043

• The Genetic Landscape of BCL2 Break Negative Follicular Lymphoma

  Background: 15% of high and 50% of low stage FL lack the t(14;18). Nevertheless, most of these cases express BCL2 by so far unknown mechanisms. Principals/Methodology: To find molecular mechanisms, triggering FL-pathogenesis in t(14;18)-negative FL, the exonic SNV-profiles of 33 FL with and 16 FL without t(14;18) were assessed by whole exome-sequencing and external SNV-data integration. The SNV-profile of the whole exome-sequenced FL was correlated with copy number, LOH and gene-expression data. Finally, the N-glycosylation status was assessed (LympTrack assay) in a separate validation cohort. Results: SNVs in FL with and without t(14;18) affected different genes, but accumulated in similar pathways, e.g. apoptosis, epigenetic processes and cell cycle/proliferation. Mutations in genes associated with immune-response were more prominently and those associated with NFkB and N-glycosylation selectively enriched in t(14;18)-negative FL. ~50% of mutated genes showed BCL2-asociation at high confidence in FL with and without t(14;18). Moreover, among the genes mutated in t(14;18)-negative FL 555 (~1/3) were affected by copy number alterations and/or copy-neutral LOH, 125 were differently expressed between FL with and without t(14;18) (p&lt;0.01) and 62 were differentially expressed and affected by copy number alterations and/or copy-neutral LOH. N-glycosylation was strongly reduced in t(14;18)-negative FL as compared to results reported for historical FL cohorts. Conclusion: These results suggest a stronger crosstalk with the microenvironment that may compensate for the lack of N-glycosylation and a diverse portfolio of SNVs that may contribute to the upregulation of BCL2 and FL pathogenesis in t(14;18)-negative FL.

  Study EGAS00001002164

• THAP11 mutations in a patient with a cblX-like phenotype implicates THAP11 in the regulation of cobalamin metabolism and early vertebrate development

  Mutations in HCFC1 are associated with cblX (MIM309541), an X-linked recessive disorder, characterized by defects in cobalamin metabolism and other developmental defects. HCFC1 is a transcriptional co-regulator, which interacts with transcription factors to regulate the expression of a myriad of genes. Notably, HCFC1 regulates cobalamin metabolism via the regulation of MMACHC expression through its interaction with THAP11, a THAP domain-containing transcription factor. In addition, the HCFC1/THAP11 complex potentially regulates genes involved in diverse cellular functions including cell cycle, proliferation, and transcription. Thus, it is likely that mutations in THAP11 may also result in biochemical and other phenotypes similar to those observed in patients with cblX. We report a patient who presented with phenotypic features that overlap cblX, but did not have any mutations in either MMACHC or HCFC1. We sequenced THAP11 by Sanger sequencing and discovered a potentially pathogenic, homozygous variant in THAP11, c.240C&gt;G (p.Phe80Leu). Functional analysis in the developing zebrafish embryo, demonstrate that both THAP11 and HCFC1 regulate the proliferation and differentiation of neural precursors, suggesting important roles in normal brain development. Further, the loss of THAP11 in zebrafish embryos, results in craniofacial abnormalities including the complete loss of Meckel’s cartilage, the ceratohyal, and all of the ceratobranchial cartilages. These data are consistent with our previous work which demonstrated a role for HCFC1 in vertebrate craniofacial development. Furthermore, our high throughput RNA-sequencing analysis reveals several overlapping gene targets of HCFC1 and THAP11 that likely contribute to the phenotypes associated with cblX and related disorders. Thus, both HCFC1 and THAP11 play important roles in the regulation of cobalamin metabolism as well as other pathways involved in early vertebrate development.

  Study EGAS00001002201

• Suppressors and activators of JAK-STAT signaling at diagnosis and relapse of acute lymphoblastic leukemia in Down Syndrome

  Children with Down Syndrome (DS) are prone to development of high risk B-cell precursor ALL (DS-ALL) that differs genetically from most sporadic pediatric ALLs. Increased expression of CRLF2, the receptor to thymic stromal lymphopoietin (TSLP) characterizes about half of DS-ALLs and also a subgroup of sporadic “Philadelphia-like ALLs”. To understand the pathogenesis of relapsed DS-ALL we performed integrative genomic analysis of 25 matched diagnosis remission and relapse DS-ALLs. We found that the CRLF2 rearrangements are early events during DS-ALL evolution and generally stable between diagnoses and relapse. Secondary activating signaling events in JAK-STAT/RAS pathway were ubiquitous but highly redundant between diagnosis and relapse, suggesting that signaling is essential but that no specific mutations are “relapse driving”. We further found that activated JAK2 may be naturally suppressed in 25% of CRLF2 positive DS-ALL by loss-of-function aberrations in USP9X, a deubuitinase previously shown to stabilize the activated phosphorylated JAK2. Interrogation of large ALL genomic databases extended our findings up to 25% of CRLF2 positive, “Ph like” ALLs. Pharmacological or genetic inhibition of USP9X as well as treatment with low-dose Ruxolitinib enhanced the survival of pre-B ALL cells overexpressing mutated JAK2. Thus, somehow counterintuitive, we found that suppression of JAK-STAT “hyper-signaling” may be beneficial to leukemic B cell precursors. This and the reduction of JAK mutated clones at relapse suggest that the therapeutic effect of JAK specific inhibitors may be limited. Rather, combined signaling inhibitors or direct targeting of the TSLP receptor may be a useful therapeutic strategy for DS-ALL.

  Study EGAS00001002410

• Genomic landscape of oral cancers (Illumina RNA-Seq)

  Human papillomavirus (HPV) is a necessary but insufficient cause of a subset of oral squamous cell carcinomas (OSCC) that is increasing markedly in frequency. To identify contributory, secondary genetic alterations in these cancers, we used comprehensive genomics methods to compare 149 HPV-positive and 335 HPV-negative OSCC tumor/normal pairs. Different behavioral risk factors underlying the two OSCC types were reflected in distinctive genomic mutational signatures. In HPV-positive OSCC, the signatures of APOBEC cytosine deaminase editing, associated with anti-viral immunity, were strongly linked to overall mutational burden. By contrast, in HPV-negative OSCC, T&gt;C substitutions in the sequence context 5’-ATN-3’ correlated with tobacco exposure. Universal expression of HPV E6*1 and E7 oncogenes was a sine qua non of HPV-positive OSCC. Significant enrichment of somatic mutations was confirmed or newly identified in PIK3CA, KMT2D, FGFR3, FBXW7, DDX3X, PTEN, TRAF3, RB1, CYLD, RIPK4, ZNF750, EP300, CASZ1, TAF5, RBL1, IFNGR1 and NFKBIA. Of these, many affect host pathways already targeted by HPV oncoproteins, including the p53 and pRB pathways, or disrupt host defenses against viral infections, including interferon and nuclear factor-B signaling. Frequent copy number changes were associated with concordant changes in gene expression. Chr. 11q (including CCND1) and 14q (including DICER1 and AKT1) were recurrently lost in HPV-positive OSCC, in contrast to their gains in HPV-negative OSCC. High-ranking variant allele fractions implicated ZNF750, PIK3CA and EP300 mutations as candidate driver events in HPV-positive cancers. We conclude that virus-host interactions cooperatively shape the unique genetic features of these cancers, distinguishing them from their HPV-negative counterparts.

  Study EGAS00001003237

• Targeting the DNA Repair Pathway in Ewing Sarcoma

  Ewing’s Sarcoma is a bone and soft tissue tumor that primarily affects adolescents and young adults. With current therapies, 70% of patients with localized disease survive but survival for metastatic and recurrent disease is poor. Whole genome sequencing of 19 Ewing’s Sarcoma tumors showed that STAG2 was mutated in 10% (2/19) of the tumors and STAG2 protein was absent in 14% (13/106) tumors by immunohistochemical staining. Previous studies have shown that glioblastoma cells lacking STAG2 are more sensitive to poly-ADP ribose polymerase (PARP) inhibitors. We found that Ewing’s Sarcoma cell lines are sensitive to PARP inhibitors irrespective of STAG2 protein expression. Ewing Sarcoma cell lines are defective in double strand DNA break repair. PARP inhibitor cytotoxicity in Ewing’s Sarcoma cells was potentiated 10-1,000 fold by DNA damaging agents (irinotecan and temozolomide). To extend these studies in vivo, we developed an orthotopic Ewing’s Sarcoma mouse model and performed pharmacokinetic and pharmacodynamic studies with three different PARP inhibitors (BMN-673, olaparib and veliparib) in clinical development for pediatric cancer. Those data were used to design preclinical phase I studies to identify tolerable drug combinations for pilot efficacy testing (preclinical phase II). Based on the results of the preclinical phase I/II data, we performed a double blind, randomized, placebo controlled preclinical phase III trial with 274 mice in 15 treatment groups. Irinotecan administered in a low-dose protracted schedule optimized for pediatric patients was an effective DNA damaging agent to combine with olaparib and BMN-673 and was better tolerated than combinations with temozolomide. Combining olaparib or BMN-673 with irinotecan and temozolomide gave complete and durable responses in over 80% of the mice.

  Study EGAS00001000839

• Genetic characterization of B-cell prolymphocytic leukemia: a hierarchical prognostic model involving MYC and TP53 abnormalities - WXS

  B-PLL is a rare disease, whose molecular pathology is largely unknown. We report here the cytogenetic and molecular findings in a large series of 34 B-PLL. Karyotype (K) was complex (≥3 abnormalities) in 73%, and highly complex (HCK≥5) in 45%. The most frequent chromosomal aberrations were: translocation targeting the MYC gene [t(MYC)] (62%), 17p deletion including TP53 gene (38%), trisomy 18/18q (30%), 13q14 deletion (29%), trisomy 3 (24%), trisomy 12 (24%) and 8p deletion (23%). Whole-Exome Sequencing performed in 16 patients revealed recurrent mutations in TP53 (6/16, 38%), MYD88 (n=4), BCOR (n=4), MYC (n=3), SF3B1 (n=3), FAT1 (n=3), SETD2 (n=2), CHD2 (n=2), CXCR4 (n=2) and BCLAF1 (n=2). The main group of patients (21/34, 62%) had a t(MYC) associated with a higher percentage of prolymphocytes (p=0.03), CD38 expression (p&lt;0.001), lower K complexity (p=0.0004), mutations in MYC and in genes involved in RNA metabolism and chromatin remodeling. Principal component analysis of gene expression data showed that patients with t(MYC) clustered together. A second group with MYC gain (5/34, 15%), was associated with HCK≥5 (p=0,01) and trisomy 3 (p=0,008). Altogether, 26/34 patients (76%) had a MYC activation, translocation or gain, that were mutually exclusive. We identified 3 distinct cytogenetic prognostic groups (p=0.0006): lower risk: absence of MYC activation (median not reached); intermediate risk: MYC activation without del17p (125 months); high risk: MYC activation + del17p (11 months). Our results show that cytogenetic analysis is a useful diagnostic tool in B-PLL that improves prognostic stratification.

  Study EGAS00001003275

• The Immune Microenvironment, Genome–Wide Copy Number Aberrations and Survival in Mesothelioma

  IntroductionResults of recent clinical studies of immune checkpoint inhibitors in malignant pleural mesothelioma (MPM) have dampened initial enthusiasm. However, the immune environment and targets of these treatments such as program cell death protein 1 (PD-1) and its ligand PD-L1 have not been well characterised in MPM. Using a large cohort of patients, we investigated PD-L1 expression, immune infiltrates and genome-wide copy number status and correlated them to clinicopathological features.Patients and methodsTissue microarrays were constructed and stained with PD-L1(clone E1L3N, CST), CD4, CD8 and FOXP3 antibodies. PD-L1 positivity was defined as ≥5% membranous staining regardless of intensity and PD-L1+hi as ≥50%. Genomic DNA from a representative subset of 113 patients was used for genome-wide copy number analysis. The percent genome alteration (PGA) was computed as a proxy for genomic instability, and statistical analyses were used to relate copy number aberrations (CNA) to other variables.ResultsAmongst 329 patients evaluated, PD-L1+ was detected in 130/311 (41.7%), but PD-L1+hi was only seen in 30/311 (9.6%). PD-L1+ correlated with non–epithelioid histology and increased infiltration with CD4+, CD8+ and FOXP3+ lymphocytes. PD-L1+hi expression correlated with worse prognosis (HR = 2.37; 95%CI= 1.57-3.56; P &lt; 0.001) on univariate analysis, but not in multivariate analysis. Higher PGA was associated with epithelioid histology and poorer survival (HR = 1.59; 95%CI= 1.01-2.5; P = 0.04), but not PD-L1 expression.ConclusionPD-L1 expression was associated with non-epithelioid MPM, poor clinical outcome and increased immunological infiltrates. Increased genomic instability did not correlate with PD-L1 expression, but was associated with poorer survival.

  Study EGAS00001002323

• Genomic landscape of oral cancers (Complete Genomics WGS)

  Human papillomavirus (HPV) is a necessary but insufficient cause of a subset of oral squamous cell carcinomas (OSCC) that is increasing markedly in frequency. To identify contributory, secondary genetic alterations in these cancers, we used comprehensive genomics methods to compare 149 HPV-positive and 335 HPV-negative OSCC tumor/normal pairs. Different behavioral risk factors underlying the two OSCC types were reflected in distinctive genomic mutational signatures. In HPV-positive OSCC, the signatures of APOBEC cytosine deaminase editing, associated with anti-viral immunity, were strongly linked to overall mutational burden. By contrast, in HPV-negative OSCC, T&gt;C substitutions in the sequence context 5’-ATN-3’ correlated with tobacco exposure. Universal expression of HPV E6*1 and E7 oncogenes was a sine qua non of HPV-positive OSCC. Significant enrichment of somatic mutations was confirmed or newly identified in PIK3CA, KMT2D, FGFR3, FBXW7, DDX3X, PTEN, TRAF3, RB1, CYLD, RIPK4, ZNF750, EP300, CASZ1, TAF5, RBL1, IFNGR1 and NFKBIA. Of these, many affect host pathways already targeted by HPV oncoproteins, including the p53 and pRB pathways, or disrupt host defenses against viral infections, including interferon and nuclear factor-B signaling. Frequent copy number changes were associated with concordant changes in gene expression. Chr. 11q (including CCND1) and 14q (including DICER1 and AKT1) were recurrently lost in HPV-positive OSCC, in contrast to their gains in HPV-negative OSCC. High-ranking variant allele fractions implicated ZNF750, PIK3CA and EP300 mutations as candidate driver events in HPV-positive cancers. We conclude that virus-host interactions cooperatively shape the unique genetic features of these cancers, distinguishing them from their HPV-negative counterparts.

  Study EGAS00001002393

• Intra-tumor heterogeneity and clonal evolution patterns towards platinum-resistant high-grade serous ovarian cancer

  Although many high-grade serous ovarian carcinoma (HGSOC) patients are initially sensitive to the standard-of-care platinum-based therapy, intra-tumor heterogeneity may fuel selection towards platinum-resistant clones over subsequent lines of therapy. We performed whole-exome coupled to ultra-deep re-sequencing and copy number alteration (CNA) profiling on 31 matched diagnosis and relapse tumors from HGSOC patients treated with platinum. Genetic heterogeneity within tumor pairs varied extremely, with some recurrent tumors being identical and others completely different from the primary tumor. All primary tumors showed evidence of homologous recombination (HR)-deficiency, but none of them carried reactivating mutations in the HR pathway at relapse, also not when becoming resistant to subsequent lines of platinum. We did also not identify other mutations, CNAs or pathways that were enriched in platinum-resistant clones. However, a platinum score consisting of 13 CNAs correlated significantly with platinum-free interval (PFI) after first-line therapy in multiple cohorts of &gt;450 primary HGSOCs. A relative increase of this score in the relapsed tumor correlated with a change in PFI during subsequent therapy. Clonal evolution under platinum therapy contributes to genetic heterogeneity within tumor pairs, implicating that genetic tests linked to targeted treatment options for recurrent HGSOC need to be performed on biopsies collected at relapse. Platinum resistance of an initially platinum-sensitive HR-deficient HGSOC at relapse is not caused by reactivation of the HR pathway, suggesting that these tumors are still responsive to PARP inhibitors, instead, several chromosomal aberrations contribute to a score predictive of disease outcome under platinum therapy, both in first and subsequent lines of platinum therapy.

  Study EGAS00001001244

• UK10K NEURO EDINBURGH

  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 subjects with schizophrenia recruited from psychiatric in-patient and out-patient facilities in Scotland. All diagnoses are based on standard research procedures and family histories are available. Patients have IQ&gt;70 and the cohort includes the following groups: 100 cases with detailed clinical, cognitive and structural and functional neuroimaging phenotypes; 138 familial cases who are the probands of families where DNA has been collected from other affected members; 162 unrelated individuals. In most cases patients and their families may be re contacted to take part in further studies.For further information on this cohort please contact Andrew McIntosh (andrew.mcintosh@ed.ac.uk).

  Study EGAS00001000117

• Targeted Gene Panel for 171 PTCLs

  Mature T-cell lymphomas, including peripheral T-cell lymphoma (PTCL) and extranodal NK/T-cell lymphoma (NKTL), represent a heterogeneous group of non-Hodgkin lymphomas with dismal outcomes and limited treatment options. To determine the extent of involvement of the JAK/STAT pathway in this malignancy, we performed targeted capture sequencing of 188 genes in this pathway in 171 PTCL and NKTL cases. A total of 272 nonsynonymous somatic mutations in 101 genes were identified in 73% of the samples, including 258 single-nucleotide variants and 14 insertions or deletions. Recurrent mutations were most frequently located in STAT3 and TP53 (15%), followed by JAK3 and JAK1 (6%) and SOCS1 (4%). A high prevalence of STAT3 mutation (21%) was observed specifically in NKTL. Novel STAT3 mutations (p.D427H, E616G, p.E616K, and p.E696K) were shown to increase STAT3 phosphorylation and transcriptional activity of STAT3 in the absence of cytokine, in which p.E616K induced programmed cell death-ligand 1 (PD-L1) expression by robust binding of activated STAT3 to the PD-L1 gene promoter. Consistent with these findings, PD-L1 was overexpressed in NKTL cell lines harboring hotspot STAT3 mutations, and similar findings were observed by the overexpression of p.E616K and p.E616G in the STAT3 wild-type NKTL cell line. Conversely, STAT3 silencing and inhibition decreased PD-L1 expression in STAT3 mutant NKTL cell lines. In NKTL tumors, STAT3 activation correlated significantly with PD-L1 expression. We demonstrated that STAT3 activation confers high PD-L1 expression, which may promote tumor immune evasion. The combination of PD-1/PD-L1 antibodies and STAT3 inhibitors might be a promising therapeutic approach for NKTL, and possibly PTCL.

  Study EGAS00001002740

• Exploiting immune cell receptor information to quantify index switching in single cell transcriptome sequencing experiment

  By offering high sequencing speed and ultra-high-throughput at a low price, Illumina next-generation sequencing platforms has been widely adopted in recent years. However, an experiment with multiplexed library could be at risk of molecular recombination, known as “index switching”, which causes certain amount of reads to be assigned to an incorrect library. It is reported that a new advance, exclusion amplification (ExAmp) in conjunction with the patterned flow cell technology introduced on HiSeq 3000/HiSeq 4000/HiSeq X sequencing systems, potentially suffers from a higher rate of index switching than conventional bridge amplification. We took advantage of the diverse but highly cell-unique expression of immune cell receptors to quantify index switching on single cell RNA-seq data that were sequenced on HiSeq 3500 and HiSeq 4000. By utilizing the unique T-cell receptor (TCR) expression, we could quantify the spread-of-signal from many different wells (n=51 from total three batches) due to index switching. We used TraCer to reconstruct full-length TCR from all samples, and then used Kallisto to quantify TCR gene expression. We found index switching in all three batches of experiments investigated. The median percentage of incorrectly detected markers was estimated to be 4.2% (interquartile range (IQR): 2.0%-8.7%). We did not detect any consistent pattern of certain indices to be more prone for switching than others, suggesting that index switching is a stochastic process. We confirm that index switching is a problem that affects all samples run in multiplexed libraries on Illumina HiSeq 3500 and HiSeq 4000 platforms.

  Study EGAS00001002911

• Integrative genomic analysis reveals cancer-associated mutations at diagnosis of CML in patients with high risk disease

  Genomic events associated with poor outcome in chronic myeloid leukemia (CML) are poorly understood. We performed whole exome sequencing, copy number variation and/or RNA-Seq for 65 patients to discover mutations at diagnosis and blast crisis (BC). Forty-six chronic phase patients with the extremes of outcome were studied at diagnosis. Cancer gene variants were detected in 15/27 patients (56%) with subsequent BC or poor outcome and in 3/19 optimal responders (16%), P=.007. Frequently mutated genes at diagnosis were ASXL1, IKZF1 and RUNX1. The methyltransferase SETD1B was a novel recurrently mutated gene. A novel class of variant associated with the Philadelphia translocation was detected at diagnosis in 11/46 patients (24%) comprising fusions and/or rearrangement of genes on the translocated chromosomes, with evidence of fragmentation, inversion and imperfect sequence reassembly. These were more frequent at diagnosis in patients with poor outcome: 9/27 (33%) versus 2/19 optimal responders (11%), P=.07. Thirty-nine patients were tested at BC and all had cancer gene variants, including ABL1 kinase domain mutations in 58%. However, ABL1 mutations co-occurred with other mutated cancer genes in 89% of cases, and these predated ABL1 mutations in 62% of evaluable patients. Gene fusions not associated with the Philadelphia translocation occurred in 42% of patients at BC and commonly involved fusion partners that were known cancer genes (78%). Genomic analysis revealed numerous relevant variants at diagnosis in patients with poor outcome and all patients at BC. Future refined biomarker testing of specific variants will likely provide prognostic information to facilitate a risk-adapted therapeutic approach.

  Study EGAS00001003071

• Whole exome sequencing of small cell neuroendocrine cancer of the cervix

  In order to improve treatment selection for high grade neuroendocrine carcinomas of the cervix (NECC), we performed a comparative genomic analysis between this rare tumor type and other cervical cancer types, as well as extra-cervical neuroendocrine small cell carcinomas of the lung and bladder. We performed whole exome sequencing on fresh-frozen tissue from 15 NECCs and matched normal tissue. We then identified mutations and copy number variants using standard analysis pipelines. Published mutation tables from cervical cancers and extra-cervical small cell carcinomas were used for comparative analysis. Descriptive statistical methods were used and a two-sided threshold of P &lt; .05 was used for significance. In the NECC cohort, we detected a median of 1.7 somatic mutations per megabase (range 1.0-20.9). PIK3CA p.E545K mutations were the most frequency observed oncogenic mutation (4/15 tumors, 27%). Activating MAPK pathway mutations in KRAS (p.G12D) and GNAS (p.R201C) co-occurred in two tumors (13%). In total we identified PI3-kinase or MAPK pathway activating mutations in 67% of NECC. When compared to NECC, lung and bladder small cell carcinomas exhibited a statistically significant higher rate of coding mutations (P &lt; .001 for lung; P = .001 for bladder). Mutation of TP53 was uncommon in NECC (13%) and was more frequent in both lung (103 of 110 tumors [94%], P &lt; .001) and bladder (18 of 19 tumors [95%], P &lt; .001) small cell carcinoma. These comparative genomics data suggest that NECC may be genetically more similar to common cervical cancer subtypes than to extra-cervical small cell neuroendocrine carcinomas of the lung and bladder. These results may have implications for the selection of cytotoxic and targeted therapy regimens for this rare disease.

  Study EGAS00001003142

• Combined MEK and BCL-2/XL inhibition is effective in high-grade serous ovarian cancer patient-derived xenograft models and BIM levels are predictive of responsiveness

  Most patients with late stage high-grade serous ovarian cancer (HGSOC) initially respond to chemotherapy but inevitably relapse and develop resistance, highlighting the need for novel therapies to improve patient outcomes. The MEK/ERK pathway is activated in a large subset of HGSOC, thus making it an attractive therapeutic target. Here, we systematically evaluated the extent of MEK/ERK pathway activation and efficacy of pathway inhibition in a large panel of well-annotated HGSOC patient-derived xenograft (PDX) models. The vast majority of models were nonresponsive to the MEK inhibitor cobimetinib (GDC-0973) despite effective pathway inhibition. Proteomic analyses of adaptive responses to GDC-0973 revealed that GDC-0973 upregulated the pro-apoptotic protein BIM, thus priming the cells for apoptosis regulated by BCL2-family proteins. Indeed, combination of both MEK inhibitor and dual BCL-2/XL inhibitor (ABT-263) significantly reduced cell number, increased cell death and displayed synergy in vitro in most models. In vivo, the GDC-0973 and ABT-263 combination was well tolerated and resulted in greater tumor growth inhibition than single agents. Detailed proteomic and correlation analyses identified two subsets of responsive models – those with high BIM at baseline that was increased with MEK inhibition and those with low basal Bim and high pERK levels. Models with low BIM and low pERK were non-responsive. Our findings demonstrate that combined MEK and BCL-2/XL inhibition has therapeutic activity in HGSOC models and provide a mechanistic rationale for clinical evaluation of this drug combination as well as the assessment of the extent to which BIM and/or pERK levels predict drug combination effectiveness in chemoresistant HGSOC.

  Study EGAS00001003427

• Plasma DNA aberrations in systemic lupus erythematosus revealed by genomic and methylomic sequencing

  We performed a high-resolution analysis of the biological characteristics of plasma DNA in systemic lupus erythematosus (SLE) patients using massively parallel genomic and methylomic sequencing. A number of plasma DNA abnormalities were found. First, aberrations in measured genomic representations (MGRs) were identified in the plasma DNA of SLE patients. The extent of the aberrations in MGRs correlated with anti-double–stranded DNA (anti-dsDNA) antibody level. Second, the plasma DNA of active SLE patients exhibited skewed molecular size-distribution profiles with a significantly increased proportion of short DNA fragments. The extent of plasma DNA shortening in SLE patients correlated with the SLE disease activity index (SLEDAI) and anti-dsDNA antibody level. Third, the plasma DNA of active SLE patients showed decreased methylation densities. The extent of hypomethylation correlated with SLEDAI and anti-dsDNA antibody level. To explore the impact of anti-dsDNA antibody on plasma DNA in SLE, a column-based protein G capture approach was used to fractionate the IgG-bound and non–IgG-bound DNA in plasma. Compared with healthy individuals, SLE patients had higher concentrations of IgG-bound DNA in plasma. More IgG binding occurs at genomic locations showing increased MGRs. Furthermore, the IgG-bound plasma DNA was shorter in size and more hypomethylated than the non–IgG-bound plasma DNA. These observations have enhanced our understanding of the spectrum of plasma DNA aberrations in SLE and may provide new molecular markers for SLE. Our results also suggest that caution should be exercised when interpreting plasma DNA-based noninvasive prenatal testing and cancer testing conducted for SLE patients.

  Study EGAS00001000962

• A molecular cell atlas of the human lung from single cell RNA sequencing

  Although single cell RNA sequencing studies have begun providing compendia of cell expression profiles, it has proven more difficult to systematically identify and localize all molecular cell types in individual organs to create a full molecular cell atlas. Here we describe droplet- and plate-based single cell RNA sequencing applied to ~75,000 human lung and blood cells, combined with a multi-pronged cell annotation approach, which have allowed us to define the gene expression profiles and anatomical locations of 58 cell populations in the human lung, including 41 of 45 previously known cell types or subtypes and 14 new ones. This comprehensive molecular atlas elucidates the biochemical functions of lung cell types and the cell-selective transcription factors and optimal markers for making and monitoring them; defines the cell targets of circulating hormones and predicts local signaling interactions including sources and targets of chemokines in immune cell trafficking and expression changes on lung homing; and identifies the cell types directly affected by lung disease genes and respiratory viruses. Comparison to mouse identified 17 molecular types that appear to have been gained or lost during lung evolution and others whose expression profiles have been substantially altered, revealing extensive plasticity of cell types and cell-type-specific gene expression during organ evolution including expression switches between cell types. This atlas provides the molecular foundation for investigating how lung cell identities, functions, and interactions are achieved in development and tissue engineering and altered in disease and evolution.

  Study EGAS00001004344

• We evaluate the potential for routine WGS using ONT by sequencing the well-characterised reference sample NA12878 and the genome of an individual with ataxia-pancytopenia syndrome accompanied by severe immune dysregulation.

  Whole-genome sequencing (WGS) is becoming widely used in clinical medicine in diagnostic contexts and to inform treatment choice. While current sequencing technologies have been extremely successful, their reliance on short read lengths necessarily involves some limitations: accurately assaying certain genomic regions and classes of variation can be problematic. Recent advances in throughput and cost have made WGS using the Oxford Nanopore Technologies (ONT) MinION long-read single-molecule sequencer a potential solution to these challenges. Here we evaluate the potential for routine WGS using ONT by sequencing the well-characterised reference sample NA12878 and the genome of an individual with ataxia-pancytopenia syndrome accompanied by severe immune dysregulation, to 82× and 30× respectively. For NA12878, we evaluated single-nucleotide variant (SNV) calls based on data from multiple base-calling algorithms. We demonstrate that phasing metrics from a novel, reference panel-free, long-read-based method can improve variant-calling performance from otherwise modest levels, resulting in a false discovery rate of 7.1% and false negative rate of 8.5%; remaining errors are concentrated near homopolymers and in regions of reduced sequencing coverage. In the clinical sample, we are able to identify and directly phase two non-synonymous de novo variants in SAMD9L, (OMIM #159550) inferring that they both lie on a common paternal haplotype by overlap with parental genotypes at nearby common variants. This work demonstrates that methodological innovation can substantially reduce variant-calling error rates in ONT data, and that with on-going improvements in throughput, base-calling and dedicated long-read-based SNV-calling methodology, ONT offers promise as an option for clinical WGS.

  Study EGAS00001003469

• Post_Mortem_Tissue_COVID19_RNA

  Single cell analysis of post mortem tissue samples from SARS-CoV2-infected patients We aim to directly identify the human cell types infected by SARS-CoV-2 and measure the cellular response to COVID19 infection across 20 different tissues from infected patient autopsies. We have examined the expression pattern of viral entry receptors across healthy human tissues to predict several candidate target cell types across the airway and heart. In addition, high prevalence of cardiac failure and abnormal renal function in COVID19 patients implicates heart and kidney involvement, but the pathogenesis of organ specific damage - whether via a direct cytopathic mechanism or an indirect inflammatory response - remains unknown . Currently, we lack confirmation of target cell types and cellular processes in infected tissues as autopsies are discouraged in most countries due to health and safety risks. Our collaborators Drs Michael Osborn and Brian Hanley (Imperial College) have outlined guidelines to perform post-mortem in COVID19 patients (Hanley et al., 2020) and have established a programme of autopsies for research to be performed in a high-risk facility at Westminster Public Mortuary. Here, we propose to identify infected cell types and aberrant molecular pathologies in this precious tissue resource using single cell and spatial genomics. We will prioritise three organ systems: the human airway, heart and the kidney. We will directly examine the cellular identities of SARS-CoV-2 infected cell types and identify the cellular responses to infection across these organs. This fundamental knowledge will help guide future treatment choices for COVID19. 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 EGAS00001004442

• ctDNA monitoring using patient-specific sequencing and integration of variant reads - Lung cohort

  Circulating tumor-derived DNA (ctDNA) can be used to monitor cancer dynamics noninvasively. Detection of ctDNA can be challenging in patients with low-volume or residual disease, where plasma contains very few tumor-derived DNA fragments. We show that sensitivity for ctDNA detection in plasma can be improved by analyzing hundreds to thousands of mutations that are first identified by tumor genotyping. We describe the INtegration of VAriant Reads (INVAR) pipeline, which combines custom error-suppression methods and signal-enrichment approaches based on biological features of ctDNA. With this approach, the detection limit in each sample can be estimated independently based on the number of informative reads sequenced across multiple patient-specific loci. We applied INVAR to custom hybrid-capture sequencing data from 176 plasma samples from 105 patients with melanoma, lung, renal, glioma, and breast cancer across both early and advanced disease. By integrating signal across a median of &gt;105 informative reads, ctDNA was routinely quantified to 1 mutant molecule per 100,000, and in some cases with high tumor mutation burden and/or plasma input material, to individual parts per million. This resulted in median Area Under the Curve (AUC) values of 0.98 in advanced cancers, and 0.80 in early stage and challenging settings for ctDNA detection. We generalized this method to whole-exome and whole-genome sequencing, showing that the INVAR may be applied without requiring personalized sequencing panels, so long as a tumor mutation list is available. As tumor sequencing becomes increasingly performed, such methods for personalized cancer monitoring may enhance the sensitivity of cancer liquid biopsies.

  Study EGAS00001004447

• UK10K NEURO UKSCZ

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. These samples have been collected from throughout the UK and Ireland. The samples fall into two main categories, approximately 500 have a full diagnostic work up. A proportion of these are cases with a positive family history of schizophrenia, either collected as sib-pairs or from multiplex kindred's. The second group consist mainly of &gt;300 samples that have been systematically collected within South Wales and in addition to full diagnostic work up have undergone detailed cognitive testing. All samples have obtained a DSM IV diagnosis of schizophrenia or schizoaffective disorder.For further information on this cohort please contact Michael Owen (OwenMJ@cardiff.ac.uk).

  Study EGAS00001000123

• Clinical Outcomes by Tumor Mutational Burden and Inflammatory Gene Expression With Combined Nivolumab and Ipilimumab or Monotherapy in Advanced Melanoma

  Purpose Exploratory analyses of CheckMate 066 and 067 trials were conducted to investigate associations of tumor mutational burden (TMB), a 4-gene inflammatory gene expression signature, and BRAF mutation status with tumor response, progression-free survival (PFS), and overall survival (OS) in patients with advanced melanoma. Patients and Methods Patients with known programmed death ligand 1 (PD-L1) expression and BRAF mutation status received nivolumab (NIVO) or dacarbazine in CheckMate 066 and either NIVO, ipilimumab (IPI), or NIVO+IPI in CheckMate 067. Whole exome sequencing and RNA sequencing were used to determine TMB and inflammatory gene expression signature scores, respectively. These biomarkers were evaluated in terms of their association with PFS and OS. Results In the NIVO, NIVO+IPI, and IPI arms of CheckMate 067, longer survival was associated with high (&gt; median) versus low (≤ median) TMB with hazard ratios (HRs) (95% confidence interval [CI]) for PFS of 0.45 (0.30–0.65), 0.55 (0.38–0.81), and 0.60 (0.43–0.82), and for OS of 0.46 (0.30–0.71), 0.53 (0.34–0.82), and 0.52 (0.36–0.74), respectively. For NIVO-treated patients, these results were confirmed in CheckMate 066. A survival benefit was observed with high TMB and absence of BRAF mutation. Survival was associated with high versus low inflammatory signature scores with HRs (95% CI) for PFS of 0.56 (0.34–0.94), 0.40 (0.23–0.72), and 0.43 (0.27–0.70), and for OS of 0.37 (0.20–0.66), 0.38 (0.19–0.74), and 0.46 (0.27–0.79), in the NIVO, NIVO+IPI, and IPI arms, respectively. Weak correlations were observed between PD-L1, TMB, and the inflammatory signature. Conclusions Combined assessment of TMB, inflammatory gene expression signature, and BRAF mutation status may be predictive for response to immunotherapy in advanced melanoma.

  Study EGAS00001004564

• Clinical Outcomes by Tumor Mutational Burden and Inflammatory Gene Expression With Combined Nivolumab and Ipilimumab or Monotherapy in Advanced Melanoma - CM066-WES

  Purpose Exploratory analyses of CheckMate 066 and 067 trials were conducted to investigate associations of tumor mutational burden (TMB), a 4-gene inflammatory gene expression signature, and BRAF mutation status with tumor response, progression-free survival (PFS), and overall survival (OS) in patients with advanced melanoma. Patients and Methods Patients with known programmed death ligand 1 (PD-L1) expression and BRAF mutation status received nivolumab (NIVO) or dacarbazine in CheckMate 066 and either NIVO, ipilimumab (IPI), or NIVO+IPI in CheckMate 067. Whole exome sequencing and RNA sequencing were used to determine TMB and inflammatory gene expression signature scores, respectively. These biomarkers were evaluated in terms of their association with PFS and OS. Results In the NIVO, NIVO+IPI, and IPI arms of CheckMate 067, longer survival was associated with high (&gt; median) versus low (≤ median) TMB with hazard ratios (HRs) (95% confidence interval [CI]) for PFS of 0.45 (0.30–0.65), 0.55 (0.38–0.81), and 0.60 (0.43–0.82), and for OS of 0.46 (0.30–0.71), 0.53 (0.34–0.82), and 0.52 (0.36–0.74), respectively. For NIVO-treated patients, these results were confirmed in CheckMate 066. A survival benefit was observed with high TMB and absence of BRAF mutation. Survival was associated with high versus low inflammatory signature scores with HRs (95% CI) for PFS of 0.56 (0.34–0.94), 0.40 (0.23–0.72), and 0.43 (0.27–0.70), and for OS of 0.37 (0.20–0.66), 0.38 (0.19–0.74), and 0.46 (0.27–0.79), in the NIVO, NIVO+IPI, and IPI arms, respectively. Weak correlations were observed between PD-L1, TMB, and the inflammatory signature. Conclusions Combined assessment of TMB, inflammatory gene expression signature, and BRAF mutation status may be predictive for response to immunotherapy in advanced melanoma.

  Study EGAS00001004567

• A reference map of potential determinants for the human serum metabolome

  The collection of metabolites circulating in the human blood, termed the serum metabolome, contains a plethora of biomarkers and causative agents. Although the origin of specific compounds is known, we have a poor understanding of the key determinants of most metabolites. Here, we measured the levels of 1251 circulating metabolites in serum samples from a healthy human cohort of 491 individuals, and devised machine learning algorithms to predict their levels in held-out subjects based on a comprehensive profile consisting of host genetics, gut microbiome, clinical parameters, diet, lifestyle, anthropometric measurements and medication data. Notably, we obtained statistically significant predictions for over 76% of the profiled metabolites. Despite using the strict out-of-sample prediction metric, which is a lower bound for the explained variance, diet and microbiome each explained hundreds of metabolites, with over 50% of the variance explained in some metabolites. We further validated the robustness of the microbiome related associations by showing a high replication rate in two geographically independent cohorts that were not available to us when developing the algorithms. We also demonstrate that some of these interactions are causal, as some metabolites we predicted to be positively associated with bread increased in level following a randomized clinical trial of bread intervention. Microbiome-explained metabolites were enriched with unnamed metabolites, and we devised an algorithm that accurately predicts their biological pathway, finding that they mainly include food components, aromatic amino acids and secondary bile acid derivatives. Overall, our results unravel potential determinants of over 800 metabolites, paving the way towards mechanistic understanding of alterations in metabolites under different conditions and to designing interventions for manipulating circulating metabolite levels.

  Study EGAS00001004512

• Tumor-derived cell lines as pharmacogenomic models to predict therapeutic vulnerabilities in hepatocellular carcinoma

  Hepatocellular carcinoma (HCC) is a heterogeneous aggressive malignancy with low efficacy of current therapies at advanced stages. We integrated molecular and pharmacological profiling of a large panel of liver cancer cell lines (LCCL) to assess their clinical relevance as HCC preclinical models and identify new effective therapies and biomarkers of response. We performed multi-omic analysis including whole-exome, RNA and miRNA sequencing and quantification of 126 proteins in a series 34 LCCL coupled with a screening of 29 anti-cancer agents. Molecular profiles of LCCL and primary HCC were compared and we searched for molecular features associated with drug response. Our panel of LCCL faithfully recapitulated the most aggressive molecular "proliferation class" of HCC. Genomic alterations of the RAS-MAPK pathway correlated with trametinib sensitivity that was more potent than FGFR4 inhibitors in tumor cells depending on the FGF19/FGFR4 pathway. Moreover, we showed that FGF19 amplification was only functional in tumor cells that retained a hepatocyte differentiation program with expression of FGFR4, KLB and NR1H4 highlighting the complex interplay between the genomic and transcriptomic context in drug response. We also identified drugs/combination that could target efficiently the transcriptomic "progenitor subclass of HCC", as well as specific genetic contexts such as inactivating mutations in TSC1/TSC2 and TP53 associated with higher sensitivity to the mTOR inhibitor rapamycin and the AURKA inhibitor alisertib, respectively. LCCL represent relevant preclinical models for drug-biomarker discovery in HCC and enabled to identify molecular contexts linked with particular therapeutic vulnerabilities that may be useful to stratify patients in future clinical trials

  Study EGAS00001003536

• Targeting FGFR1 for treatment of soft-tissue sarcoma (H021)

  Purpose: Altered FGFR1 signaling has emerged as therapeutic target in various epithelial malignancies. In contrast, the role of FGFR1 in soft-tissue sarcoma (STS) has not been established. Prompted by the detection and subsequent therapeutic inhibition of amplified FGFR1 in a patient with metastatic leiomyosarcoma, we investigated the oncogenic properties and potential as drug target of FGFR1 in STS.Experimental Design: The frequency of FGFR1 amplification and overexpression, as assessed by fluorescence in situ hybridization, microarray-based comparative genomic hybridization and mRNA expression profiling, SNP array profiling, and RNA sequencing, was determined in three independent patient cohorts. The sensitivity of STS cell lines with or without FGFR1 alterations to genetic and pharmacologic FGFR1 inhibition and the signaling pathways engaged by FGFR1 were investigated using short-term viability assays, long-term colony formation assays, and biochemical analysis.Results: Increased FGFR1 copy number was detected in 74 of 190 (38.9%; Cohort 1), 13 of 79 (16.5%; Cohort 2), and 80 of 254 (31.5%; Cohort 3) patients. FGFR1 overexpression occurred in 16 of 79 (20.2%, Cohort 2) and 39 of 254 (15.4%; Cohort 3) patients. Targeting of FGFR1 by short hairpin RNA-mediated knockdown and different small-molecule inhibitors (PD173074, AZD4547, BGJ398) revealed that the requirement for FGFR1 signaling in STS cells is primarily dictated by FGFR1 expression levels, and identified the MAPK-ERK1/2 axis as critical FGFR1 effector pathway.Conclusion: These data identify FGFR1 as driver gene in multiple STS subtypes and support FGFR1 inhibition, guided by patient selection according to FGFR1 expression and monitoring of MAPK-ERK1/2 signaling, as therapeutic option in this challenging group of diseases.H021

  Study EGAS00001001844

• Circulating Tumor DNA in Checkpoint Inhibitor treated Lung Cancer

  BackgroundImmune checkpoint inhibitors (ICI) are increasingly used to treat non-small cell lung cancer (NSCLC) and there is a critical need for novel predictive/prognostic biomarkers. We evaluated whether dynamics of circulating tumor (ctDNA) can identify NSCLC patients who benefit from single agent ICI.MethodsCirculating cell-free DNA (cfDNA) from 177 patients was analyzed prior ICI and at first treatment evaluation using a Commercial 77-gene panel. We performed de novo mutation calling and corrected for variants related clonal hematopoiesis. ctDNA changes were evaluated as predictors of tumor response estimated with PET/diagnostic CT and survival.ResultsMutations were detected in 159 of 174 (91.4%) analyzable patients. Most responders presented decreasing ctDNA levels, while non-responders had mixed molecular responses. Patients with a &gt;50% molecular ctDNA decrease or no detectable ctDNA demonstrated improved PFS (5.9 versus 1.8 months; HR 0.49, 95%CI, 0.35-0.67, P &lt; .001), and OS (18.9 versus 9.9 months; HR 0.41, 95%CI, 0.30-0.57, P &lt; .001) compared to patients not achieving these endpoints. Clearance of mutations (and no detected ctDNA) had a similar effect (PFS HR 0.45, 95%CI, 0.32-0.61, P &lt; .001; OS 0.37,95%CI, 0.27-0.52, P &lt; .001) compared to patients having stable or additional mutations. After adjusting for clinical factors both &gt;50% ctDNA decrease (HR 0.4, 95% CI, 0.3-0.6, P &lt; .001) and clearance of mutations (HR 0.42, 95% CI, 0.3-0.6, P &lt; .001) remained independent predictors for PFS and OS.ConclusionsctDNA testing reveals significant prognostic albeit not predictive information in ICI treated advanced NSCLC. This relevant information can be obtained with a commercial easy-to-implement ctDNA test that will facilitate extensive testing.

  Study EGAS00001004847

• Characterization of a human iPSC-derived endocrine pancreas model

  Directed differentiation of stem cells offers a scalable solution to the need for human cell models recapitulating islet biology and T2D pathogenesis. We profiled mRNA expression at six stages of an induced pluripotent stem cell (iPSC) model of endocrine pancreas development from two donors, and characterized the distinct transcriptomic profiles associated with each stage. Established regulators of endodermal lineage commitment, such as SOX17 (log2 fold change [FC] compared to iPSCs=14.2, p-value=4.9x10-5) and the pancreatic agenesis gene GATA6 (log2 FC=12.1, p-value=8.6x10-5), showed transcriptional variation consistent with their known developmental roles. However, these analyses highlighted many other genes with stage-specific expression patterns, some of which may be novel drivers or markers of islet development. For example, the leptin receptor gene, LEPR, was most highly expressed in published data from in vivo-matured cells compared to our endocrine pancreas-like cells (log2 FC=5.5, p-value=2.0x10-12), suggesting a role for the leptin pathway in the maturation process. Endocrine pancreas-like cells showed significant stage-selective expression of adult islet genes, including INS, ABCC8, and GLP1R, and enrichment of relevant GO-terms (e.g. "insulin secretion"; odds ratio=4.2, p-value=1.9x10-3): however, principal component analysis indicated that in vitro-differentiated cells were more immature than adult islets. Integration of the stage-specific expression information with genetic data from T2D genome-wide association studies revealed that 46 of 82 T2D-associated loci harbor genes present in at least one developmental stage, facilitating refinement of potential effector transcripts. Together, these data show that expression profiling in an iPSC islet development model can further understanding of islet biology and T2D pathogenesis.

  Study EGAS00001001803

• Breast cancer women lack normal lifelong immune response after full-term pregnancies

  Immunological tolerance during pregnancy is necessary for women’s bodies to accept the semi-allogenic fetus. Interestingly, “immune evasion”, which consists of some similar immunological changes, was recently proposed as a hallmark of cancer. Looking at these results, it seems paradoxical that increasing number of full-term pregnancies (FTP) reduces the risk of breast cancer in women. Our aim was to use a systems epidemiology approach to investigate this paradox by modelling breast cancer incidence as a function of parity in the Norwegian Women and Cancer (NOWAC) cohort, then testing the mathematical model using gene expression from blood in a subcohort of these women: the NOWAC postgenome cohort. Lastly, we undertook a gene set enrichment analysis for immunological gene sets. A linear trend nicely fitted the dataset showing an 8% decrease in risk for each FTP, independent of stratification on other risk factors and lasting for decades after a woman’s last FTP. However, when we looked at gene expression we found that hundreds of genes showed linear trends in cancer-free controls, but this was not the case in /any/ of the samples from breast cancer cases. In the gene set enrichment analysis, samples consisted of human blood or cord blood and mouse bone marrow or spleen. We found more adaptive gene sets than expected in humans and more innate gene sets than expected in mice. None of the gene sets provided information about the long-term effects of immunological tolerance during pregnancy, but the lack of relevant reductionist experiments hampers our interpretation. We suggest that, after resolution at birth, the immunological tolerance induced by pregnancy augments immune response to cancer development.

  Study EGAS00001002616

• DSRCT RNA genomic sequencing

  Purpose: To further understand the molecular pathogenesis of desmoplastic small round cell tumor (DSRCT), a rare but often fatal malignancy occurring primarily in young males, we used next-generation RNA sequencing to investigate the gene expression profiles intrinsic to this disease.Experimental Design: RNA from DSRCT tumor samples were obtained from the Cooperative Human Tissue Network biorepository. RNA sequencing was completed using the Illumina HiSeq 2000 system. Presence of the EWSR1-WT1 fusion was verified using both RNA-sequencing and PCR. Validation studies, including WT1 ChIP-seq, RNA sequencing, and EWS-WT1 knockdown experiments, were performed on an established DSRCT cell line, JN-DSRCT-1. Protein expression of clinically relevant genes determined by RNA sequencing was confirmed using immunohistochemistry. A panel of immune signature genes was also evaluated to identify possible immune therapeutic targets. Results: We analyzed the RNA from 14 tumor samples, of which two were confirmed negative for the diagnostic EWSR1-WT1 translocation. Principal component analysis of the sequencing data indicated that fusion-negative and fusion-positive tumors are transcriptionally distinct entities. Overlapping the significant ChIP-seq peaks with genes found to be highly expressed in the RNA sequencing data revealed that IGF2 and FGFR4 were both highly expressed, and targets of the EWS-WT1 fusion gene. This finding has potential clinical relevance for targeted therapeutic intervention. In addition, we identified the immune checkpoints CD200 and CD276 as potentially targetable genes whose expression is independent of the EWS-WT1 fusion gene in vitro.Conclusions: Despite the transcriptional complexity of this disease, we were able to identify IGF2, FGFR4, CD200, and CD276 as potential therapeutic targets for patients suffering from DSRCT.

  Study EGAS00001002770

• Combining a Universal Telomerase Based Cancer Vaccine with Ipilimumab in Patients with Metastatic Melanoma - Five-year Follow up of a Phase I/IIa Trial

  Background Ipilimumab improves survival for patients with metastatic malignant melanoma. Combining a therapeutic cancer vaccine with ipilimumab may increase efficacy by providing enhanced anti-tumor immune responses. UV1 consists of three synthetic long peptides from human telomerase reverse transcriptase (hTERT). These peptides comprise epitopes recognized by T cells from cancer patients experiencing long-term survival following treatment with a first-generation hTERT vaccine, and generate long-lasting immune responses in cancer patients when used as monotherapy. The objective of this trial was to investigate the safety and efficacy of combining UV1 with ipilimumab in metastatic melanoma. Patients and methods In this phase I/IIa, single center trial [NCT02275416], patients with metastatic melanoma received repeated UV1 vaccinations, with GM-CSF as an adjuvant, in combination with ipilimumab. Patients were evaluated for safety, efficacy and immune response. Immune responses against vaccine peptides were monitored in peripheral blood by measuring antigen-specific proliferation and IFN-γ production. Results Twelve patients were recruited. Adverse events were mainly diarrhea, injection site reaction, pruritus, rash, nausea and fatigue. Ten patients showed a Th1 immune response to UV1 peptides, occurring early and after few vaccinations. Three patients obtained a partial response and one patient a complete response. Overall survival was 50% at 5 years. Conclusion Treatment was well tolerated. The rapid expansion of UV1-specific Th1 cells in the majority of patients indicates synergy between UV1 vaccine and CTLA-4 blockade. This may have translated into clinical benefit, encouraging the combination of UV1 vaccination with standard of care treatment regimes containing ipilimumab/CTLA-4 blocking antibodies.

  Study EGAS00001005253

• Post_Mortem_Tissue_COVID19_spatial

  Single cell analysis of post mortem tissue samples from SARS-CoV2-infected patients We aim to directly identify the human cell types infected by SARS-CoV-2 and measure the cellular response to COVID19 infection across 20 different tissues from infected patient autopsies. We have examined the expression pattern of viral entry receptors across healthy human tissues to predict several candidate target cell types across the airway and heart. In addition, high prevalence of cardiac failure and abnormal renal function in COVID19 patients implicates heart and kidney involvement, but the pathogenesis of organ specific damage - whether via a direct cytopathic mechanism or an indirect inflammatory response - remains unknown . Currently, we lack confirmation of target cell types and cellular processes in infected tissues as autopsies are discouraged in most countries due to health and safety risks. Our collaborators Drs Michael Osborn and Brian Hanley (Imperial College) have outlined guidelines to perform post-mortem in COVID19 patients (Hanley et al., 2020) and have established a programme of autopsies for research to be performed in a high-risk facility at Westminster Public Mortuary. Here, we propose to identify infected cell types and aberrant molecular pathologies in this precious tissue resource using single cell and spatial genomics. We will prioritise three organ systems: the human airway, heart and the kidney. We will directly examine the cellular identities of SARS-CoV-2 infected cell types and identify the cellular responses to infection across these organs. This fundamental knowledge will help guide future treatment choices for COVID19. 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 EGAS00001004441

• Spiradenocarcinoma

  • Spiradenocarcinoma is a rare cutaneous sweat gland adnexal cancer with potential for aggressive behaviour. They are classified histologically into low- and high-grade tumours, with morphologically low-grade tumours thought to behave more favourably. However, limited information is available, with only 18 published cases. • We have collected morphologically low-grade spiroadenocarcinomas (one with a lung metastasis) and high-grade spiroadenocarcinomas, as well as some spiradenomas (benign lesions), cylindromas (another type of malignant cutaneous sweat gland adnexal tumour) and hybrid spiradenoma-cylindromas. • H&amp;E-stained sections were reviewed, follow-up was obtained, and immunohistochemistry for Ki-67, p53 and, MYB has been performed. The tumours were solitary, measuring 0.8-7?cm (median: 2.7?cm), with a predilection for the head and neck of elderly patients (median age: 72 years; range 53-92) without gender bias. Histologically, the tumours were multinodular and located in deep dermis and subcutis. A pre-existing spiradenoma was present in all cases. The malignant component was characterized by expansile growth with loss of the dual cell population, up to moderate cytological atypia and increased mitotic activity (median: 10/10 HPF; range 1-28). Additional findings included squamoid differentiation (n=9), necrosis (n=7), and ulceration (n=5). P53 expression was variable and no significant differences were noted in the benign compared with the malignant parts of the tumours. In contrast, in the malignant components the Ki-67 proliferative index was slightly increased, and MYB expression was lost. Follow-up (median: 67 months; range: 13-132) available for 16 patients (84%) revealed a local recurrence rate of 19% but no metastases or disease-related mortality. • Here we wish to exome sequence these cases to define the first genomic landscape for this malignancy.

  Study EGAS00001001799

• An analysis of humoral and cellular immune responses following COVID-19 vaccination.

  Two dose mRNA vaccination provides excellent protection against SARS-CoV-2. However, there are few data on vaccine efficacy in elderly individuals above the age of 801. Additionally, new variants of concern (VOC) with reduced sensitivity to neutralising antibodies have raised fears for vulnerable groups. Here we assessed humoral and cellular immune responses following vaccination with mRNA vaccine BNT162b22 in elderly participants prospectively recruited from the community and younger health care workers. Median age was 72 years and 51% were females amongst 140 participants. Neutralising antibody responses after the first vaccine dose diminished with increasing age, with a steep drop in participants over 80 years old. Sera from participants below and above 80 showed significantly lower neutralisation potency against B.1.1.7, B.1.351 and P.1. variants of concern as compared to wild type. Those over 80 were more likely to lack any neutralisation against VOC compared younger participants following first dose. Binding IgG and IgA antibodies were lower in the elderly, as was the frequency of SARS-CoV-2 Spike specific B- cells. We observed a trend towards lower somatic hypermutation in participants with suboptimal neutralisation, and elderly participants demonstrated clear reduction in class switched somatic hypermutation, driven by the IgA1/2 isotype. SARS-CoV-2 Spike specific T- cell IFN?? and IL-2 responses were secreted primarily by CD4 T cells, and impaired in the older age group. We conclude that the elderly are a high risk population that warrant specific measures in order to mitigate against vaccine failure, particularly where variants of concern are circulating. Specifically, the dosing interval should not be extended in this group.

  Study EGAS00001005380

• IMCISION DNAseq

  Standard-of-care (SOC) surgery for locoregionally advanced head and neck squamous cell carcinoma (HNSCC) is intensive and results in 30‒50% five-year overall survival. Anti-PD1 immune checkpoint blockade (ICB) durably improves survival rates in recurrent/metastatic HNSCC. We report on the non-randomized phase Ib/IIa IMCISION trial (NCT03003637) of 32 HNSCC patients treated with two doses neoadjuvant nivolumab (NIVO MONO, n=6, phase Ib arm A) or two doses of nivolumab plus a single dose of ipilimumab (COMBO, n=26, 6 in phase Ib arm B and 20 in phase IIa) prior to surgery and, if indicated, adjuvant radiotherapy. Neoadjuvant ICB was feasible in all phase Ib patients, meeting the phase Ib primary feasibility endpoint. One phase IIa patient had progressive disease precluding surgery. Primary tumor pathological response (phase IIa primary endpoint), defined as the % change in viable tumor cells from baseline biopsy to on-treatment resection, was evaluable in 29 patients. . We observed a major pathological response (MPR, 90‒100% response) in 8/23 (35%) evaluable COMBO and 1/6 (17%) NIVO MONO patients. None of the patients with MPR after NIVO MONO or COMBO developed a tumor relapse after 24.0 months median follow-up. FDG-PET-based total lesion glycolysis identified MPR patients prior to surgery. A baseline AID/APOBEC-associated mutational profile and an on-treatment decrease in hypoxia signature expression were observed in MPR patients. Our data indicate that neoadjuvant NIVO MONO and COMBO are safe in HNSCC, and that particularly COMBO ICB shows encouraging effectivity. IMCISION provides rationales for future HNSCC trials aiming to pre-operatively select patients with a likely MPR after neoadjuvant ICB for de-escalation of SOC.

  Study EGAS00001005466

• Characterization of DLBCL with a PMBL gene expression signature

  Primary mediastinal large B-cell lymphoma (PMBL) is a type of aggressive B-cell lymphoma that typicallyaffects young adults, characterized by presence of a bulky anterior mediastinal mass. Lymphomas withgene expression features of PMBL have been described in non-mediastinal sites, raising questions abouthow these tumors should be classified.Here, we investigated whether these "non-mediastinal PMBLs" are indeed PMBLs or instead represent adistinct group within DLBCL. From a cohort of 325 de novo DLBCL cases, we identified tumors frompatients without evidence of anterior mediastinal involvement that expressed a PMBL expression signature(nm-PMBLsig-pos, n=16, 5%). The majority of these tumors expressed MAL and CD23 – proteins typicallyobserved in bona fide PMBL (bf-PMBL). Evaluation of clinical features of nm-PMBLsig-pos cases revealedclose associations with DLBCL, and the majority displayed a germinal center B-cell-like cell-of-origin(GCB). In contrast to bf-PMBL, nm-PMBLsig-pos patients presented at an older age, did not show pleuraldisease, and bone/bone marrow involvement was observed in three cases. However, while clinicallydistinct from bf-PMBL, nm-PMBL-sig-pos tumors resembled bf-PMBL at the molecular level with upregulationof immune response, JAK-STAT, and NF-kB signatures. Mutational analysis revealed frequent somatic genemutations in SOCS1, IL4R, ITPKB and STAT6, as well as CD83 and BIRC3, with the latter genes beingsignificantly more frequently affected than in GCB-DLBCL and bf-PMBL.Our data establish nm-PMBLsig-pos lymphomas as a group of DLBCL with distinct phenotypic and geneticfeatures, and potential implications for gene expression- and mutation-based subtyping of aggressive Bcell lymphoma and related targeted therapies.

  Study EGAS00001005057

• Persistence of circulating tumor DNA in breast cancer patients during neoadjuvant treatment is a significant predictor of poor tumor response

  Purpose:Accurate assessment of response during neoadjuvant systemic treatment (NST) poses a major clinical challenge. Therefore, a minimally-invasive assessment of tumor response based on cell-free circulating tumor DNA (ctDNA) may be beneficial to guide treatment decisions.Experimental Design: We profiled 93 genes in tissue from 193 early breast cancer patients. Patient-specific assays were designed to track ctDNA during NST from 145 patients with available plasma. ctDNA presence and levels were correlated with complete pathological response (pCR) and residual cancer burden (RCB) as well as with clinicopathologic characteristics of the tumor to identify potential proxies for ctDNA release. Moreover, we tested the predictive value of driver events identified in the tumor.Results:At baseline, ctDNA could be detected in 63/145 (43.4%) patients and persisted in 25/63 (39.7%) patients at mid-therapy (MT) and 15/63 (23.8%) patients at the end of treatment. ctDNA detection at MT was significantly associated with higher RCB (OR 0.062, 95% CI 0.01-0.48, P=0.0077). Out of 31 patients with detectable ctDNA at MT, 30 patients (96.8%) were non-responders (RCB II, n=8; RCB III, n=22) and only one patient responded to the treatment (RCB I). Considering all 145 patients with baseline (BL) plasma, none of the patients with RCB 0 and only 6.7% of patients with RCB I had ctDNA detectable at MT, while 30.6% and 29.6% of patients with RCB II/III, respectively, had a positive ctDNA result. Conclusion:Overall, our results demonstrate that the detection and persistence of ctDNA at mid-therapy may have the potential to negatively predict response to neoadjuvant treatment and identify patients who will not achieve pCR or be classified with RCB II/III.

  Study EGAS00001005798

• Whole genome sequencing of a breast cancer cohort with known functional homologous recombination status

  Optimal patient selection for PARP inhibitors and double strand break (DSB) inducing chemotherapy is of great clinical importance. In principle, the use of these therapies can be extended beyond germline BRCA1/2 mutated cancers. Although several methods for detection of BRCA-like tumors have been developed, a golden standard is lacking. We have developed the RECAP test, a functional HR assay exploiting the formation of RAD51 foci in proliferating cells after ex vivo irradiation of fresh primary breast cancer (BrC) tissue. The RECAP test seems robust and has an advantage over genetic tests, due to its functional character for exploring the HR phenotype. Another type of HRD test is the BRCAness classifier, which is based on specific genomic patterns derived from copy number data of BRCA1/2 mutated BrCs that also occur in sporadic cancers. The BRCAness classifier is a robust test that is easily applicable in the clinic and has also been validated to predict in vivo response to high dose chemotherapy in clinical trials. More recently, a new HRD test based on genomic signatures which can be extracted from whole genome sequencing has been published, HRDetect. Studies comparing these different HRD tests within the same cohort of patients do not yet exist. Here, we have the unique opportunity to compare three HRD tests within the same cohort. We aim to compare three different HRD tests (RECAP, BRCAness classifier and HRDetect) within the same cohort of primary BrCs (n=41). As a secondary aim, we will analyze additional genetic alterations in HRD tumors (both BRCA-deficient and BRCA-proficient HRD tumors).

  Study EGAS00001005572

• Genetic Risk for Subsequent Neoplasms among Long-term Survivors of Childhood Cancer in the St. Jude Lifetime Cohort

  While childhood cancer survivors are known to be at increased risk of subsequent neoplasms (SNs), the contribution of germline mutations in cancer predisposition genes to the risk of SN is largely unknown. Whole-genome sequencing (30X) was performed on 3,006 5+ year survivors of childhood cancer (median age, 35.8 [range: 7.1-69.8] years) from the St. Jude Lifetime Cohort. Survivors underwent a comprehensive clinical assessment. Germline mutations in 156 cancer predisposition genes (SJCPG156), including a subset of 60 (SJCPG60) autosomal dominant genes with moderate to high penetrance, were classified for their pathogenicity. Piecewise exponential regression, stratified by radiation exposure, was used to evaluate the relative rate (RR) and 95% confidence interval (95% CI) of SN occurrence by mutation status. Pathogenic or likely pathogenic mutations in SJCPG156 and SJCPG60 were identified in 11.7% and 5.8% of survivors, respectively. The most frequently mutated SJCPG60 genes included RB1 (n=43), NF1 (n=22), BRCA2 (n=14), BRCA1 (n=12) and TP53 (n=10). Mutations in SJCPG60 were associated with the rate of subsequent sarcoma (RR, 10.3; 95% CI, 4.2 to 25.6) and breast cancer (RR, 13.3; 95% CI, 5.8 to 30.5) among irradiated survivors, and the rate of developing any SN (RR, 4.8; 95% CI, 2.4 to 9.4) and breast cancer (RR, 7.8; 95% CI, 2.5 to 24.0) among non-irradiated survivors. The rate of developing ≥2 histologically distinct SNs was increased 17.8-fold (95% CI, 3.5 to 89.8) among SJCPG60 mutation carriers treated without radiotherapy. Similar but attenuated associations were observed for mutations in SJCPG156. Survivors who develop a SN without prior radiotherapy or who are diagnosed with a subsequent sarcoma or breast cancer within a site of prior radiotherapy should be referred to genetic counseling.

  Study EGAS00001002499

• Monocyte and macrophage lipid accumulation results in downregulated type-I interferon responses

  Macrophages are critical components of atherosclerotic lesions and their pro- and anti-inflammatory responses influence atherogenesis. Type-I interferons (IFNs) are cytokines that play an essential role in antiviral responses and inflammatory activation and have been shown to promote atherosclerosis. Although the impact of type-I IFNs on macrophage foam cell formation is well-documented, the effect of lipid accumulation in monocytes and macrophages on type-I IFN responses remains unknown. Here we examined IFN stimulated (ISG) and non-ISG inflammatory gene expression in mouse and human macrophages that were loaded with acetylated LDL (acLDL), as a model for foam cell formation. We found that acLDL loading in mouse and human macrophages specifically suppressed expression of ISGs and IFN-β secretion, but not other pro-inflammatory genes. The downregulation of ISGs could be rescued by exogenous IFN-β supplementation. Activation of the cholesterol-sensing nuclear liver X receptor (LXR) recapitulated the cholesterol-initiated type-I IFN suppression. Additional analyses of murine in vitro and in vivo generated foam cells confirmed the suppressed IFN signalling pathways and suggest that this phenotype is mediated via downregulation of interferon regulatory factor binding at gene promoters. Finally, RNA-seq analysis of monocytes of familial hypercholesterolemia (FH) patients also showed type-I IFN suppression which was restored by lipid-lowering therapy and not present in monocytes of healthy donors. Taken together, we define type-I IFN suppression as an athero-protective characteristic of foamy macrophages. These data provide new insights into the mechanisms that control inflammatory responses in hyperlipidaemic settings and can support future therapeutic approaches focusing on reprogramming of macrophages to reduce atherosclerotic plaque progression and improve stability.

  Study EGAS00001005955

• Genomic landscape of oral cancers (Illumina WGS)

  Human papillomavirus (HPV) is a necessary but insufficient cause of a subset of oral squamous cell carcinomas (OSCC) that is increasing markedly in frequency. To identify contributory, secondary genetic alterations in these cancers, we used comprehensive genomics methods to compare 149 HPV-positive and 335 HPV-negative OSCC tumor/normal pairs. Different behavioral risk factors underlying the two OSCC types were reflected in distinctive genomic mutational signatures. In HPV-positive OSCC, the signatures of APOBEC cytosine deaminase editing, associated with anti-viral immunity, were strongly linked to overall mutational burden. By contrast, in HPV-negative OSCC, T&gt;C substitutions in the sequence context 5’-ATN-3’ correlated with tobacco exposure. Universal expression of HPV E6*1 and E7 oncogenes was a sine qua non of HPV-positive OSCC. Significant enrichment of somatic mutations was confirmed or newly identified in PIK3CA, KMT2D, FGFR3, FBXW7, DDX3X, PTEN, TRAF3, RB1, CYLD, RIPK4, ZNF750, EP300, CASZ1, TAF5, RBL1, IFNGR1 and NFKBIA. Of these, many affect host pathways already targeted by HPV oncoproteins, including the p53 and pRB pathways, or disrupt host defenses against viral infections, including interferon and nuclear factor-B signaling. Frequent copy number changes were associated with concordant changes in gene expression. Chr. 11q (including CCND1) and 14q (including DICER1 and AKT1) were recurrently lost in HPV-positive OSCC, in contrast to their gains in HPV-negative OSCC. High-ranking variant allele fractions implicated ZNF750, PIK3CA and EP300 mutations as candidate driver events in HPV-positive cancers. We conclude that virus-host interactions cooperatively shape the unique genetic features of these cancers, distinguishing them from their HPV-negative counterparts.

  Study EGAS00001003228

• Functional Mapping of AKT Signaling and Biomarkers of Response From the FAIRLANE Trial of Neoadjuvant Ipatasertib Plus Paclitaxel for Triple-Negative Breast Cancer

  Purpose: Despite extensive genomic and transcriptomic profiling, it remains unknown how signaling pathways are differentially activated and how tumors are differentially sensitized to certain perturbations. Here, we aim to characterize AKT signaling activity and its association with other genomic or immunohistochemistry-based PI3K/AKT pathway biomarkers as well as the clinical activity of ipatasertib (AKT inhibitor) in the FAIRLANE trial. Experimental Design: In FAIRLANE, 151 patients with early triple-negative breast cancer were randomized 1:1 to receive paclitaxel with ipatasertib or placebo for 12 weeks prior to surgery. Adding ipatasertib did not increase pathologic complete response rate and numerically improved overall response rate by magnetic resonance imaging (MRI). We used reverse-phase protein microarrays (RPPA) to examine the total level and/or phosphorylation states of over 100 proteins in various signaling or cell processes including PI3K/AKT and mTOR signaling. 125 baseline and 127 on-treatment samples were evaluable by RPPA, with 110 paired samples at both time points. Results: Tumors with genomic/protein alterations in PIK3CA/AKT1/PTEN were associated with higher levels of AKT phosphorylation. In addition, phosphorylated(p)AKT levels exhibited a significant association with enriched clinical benefit of ipatasertib, and identified patients who received benefit in the absence of PIK3CA/AKT1/PTEN alterations. Ipatasertib treatment led to a down-regulation of AKT/mTORC1 signaling, which was more pronounced among the tumors with PIK3CA/AKT1/PTEN alterations or among the responders to the treatment. Conclusions: We showed that the high baseline pAKT levels are associated with the alterations of PI3K/AKT pathway components and enriched benefit of ipatasertib in TNBC.

  Study EGAS00001005892

• The clinical utility of genomics in childhood cancer extends beyond targetable mutations

  Survival of children with relapsed cancer remains unacceptably low. Cancer genomic profiling may identify new therapeutic options, enhance diagnostic accuracy, and improve understanding of tumors’ origins and continued evolution. We deeply sequenced 864 cancer-associated genes and the complete genomes and transcriptomes for 300 pediatric and adolescent/young adult (AYA) patients with poor prognosis or rare tumors. Using integrative somatic-germline analyses, we assessed the clinical utility of cancer genomics in this population. We analyzed tumor mutational signatures and burden to interrogate the relevance of germline variants in DNA repair genes and the impact of therapeutic exposures on tumor evolution. We also evaluated the frequency of changes to tumor drivers and therapeutic targets at relapse in those with serial samples. Integrative analysis yielded clinically actionable variants in a majority of patients. Improved diagnostic accuracy led to modified management in a subset. Therapeutically targetable variants were of unanticipated timing and type, with many of these derived from the germline. Enrichment in corroborating mutational signature 3 (‘BRCAness’) in patients with germline homologous recombination defects suggests they are drivers of pediatric cancers and potential targets for PARP inhibition. We show that mutation burden, which may in part be driven by prior therapeutic exposures, was significantly elevated in a minority of patients. Over one-third of patients with sequential samples analyzed showed a change in therapeutically-targetable drivers. Comprehensive somatic-germline cancer genomic profiling is useful at multiple points in the care trajectory for pediatric/AYA patients, supporting its integration into early clinical management. Re-biopsy at the time of relapse should be considered, and specific attention to mutational burden is indicated. Defective DNA repair genes in the germline may represent important targetable drivers.

  Study EGAS00001006034

• ctDNA monitoring using patient-specific sequencing and integration of variant reads - Breast cohort

  Circulating tumor-derived DNA (ctDNA) can be used to monitor cancer dynamics noninvasively. Detection of ctDNA can be challenging in patients with low-volume or residual disease, where plasma contains very few tumor-derived DNA fragments. We show that sensitivity for ctDNA detection in plasma can be improved by analyzing hundreds to thousands of mutations that are first identified by tumor genotyping. We describe the INtegration of VAriant Reads (INVAR) pipeline, which combines custom error-suppression methods and signal-enrichment approaches based on biological features of ctDNA. With this approach, the detection limit in each sample can be estimated independently based on the number of informative reads sequenced across multiple patient-specific loci. We applied INVAR to custom hybrid-capture sequencing data from 176 plasma samples from 105 patients with melanoma, lung, renal, glioma, and breast cancer across both early and advanced disease. By integrating signal across a median of &gt;105 informative reads, ctDNA was routinely quantified to 1 mutant molecule per 100,000, and in some cases with high tumor mutation burden and/or plasma input material, to individual parts per million. This resulted in median Area Under the Curve (AUC) values of 0.98 in advanced cancers, and 0.80 in early stage and challenging settings for ctDNA detection. We generalized this method to whole-exome and whole-genome sequencing, showing that the INVAR may be applied without requiring personalized sequencing panels, so long as a tumor mutation list is available. As tumor sequencing becomes increasingly performed, such methods for personalized cancer monitoring may enhance the sensitivity of cancer liquid biopsies.

  Study EGAS00001004446

• Profiling the unique protective properties of intracranial arterial endothelial cells

  Cardiovascular disorders, like atherosclerosis and hypertension, are increasingly known to be associated with vascular cognitive impairment (VCI). In particular, intracranial atherosclerosis is one of the main causes of VCI, although plaque development occurs later in time and is structurally different compared to atherosclerosis in extracranial arteries. Recent data suggest that endothelial cells (ECs) that line the intracranial arteries may exert anti-atherosclerotic effects due to yet unidentified pathways. To gain insights into underlying mechanisms, we isolated post-mortem endothelial cells from both the intracranial basilar artery (BA) and the extracranial common carotid artery (CCA) from the same individual (total of 15 individuals) with laser capture microdissection. RNA sequencing revealed a distinct molecular signature of the two endothelial cell populations of which the most prominent ones were validated by means of qPCR. Our data reveal for the first time that intracranial artery ECs exert an immune quiescent phenotype. Secondly, genes known to be involved in the response of ECs to damage (inflammation, differentiation, adhesion, proliferation, permeability and oxidative stress) are differentially expressed in intracranial ECs compared to extracranial ECs. Finally, Desmoplakin (DSP) and Hop Homeobox (HOPX), two genes expressed at a higher level in intracranial ECs, and Sodium Voltage-Gated Channel Beta Subunit 3 (SCN3B), a gene expressed at a lower level in intracranial ECs compared to extracranial ECs, were shown to be responsive to shear stress and/or hypoxia. With our data we present a set of intracranial-specific endothelial genes that may contribute to its protective phenotype, thereby supporting proper perfusion and consequently may preserve cognitive function. Deciphering the molecular regulation of the vascular bed in the brain may lead to the identification of novel potential intervention strategies to halt vascular associated disorders, such as atherosclerosis and vascular cognitive dysfunction.

  Study EGAS00001004479

• Clinical Outcomes by Tumor Mutational Burden and Inflammatory Gene Expression With Combined Nivolumab and Ipilimumab or Monotherapy in Advanced Melanoma - CM67 WES

  Purpose Exploratory analyses of CheckMate 066 and 067 trials were conducted to investigate associations of tumor mutational burden (TMB), a 4-gene inflammatory gene expression signature, and BRAF mutation status with tumor response, progression-free survival (PFS), and overall survival (OS) in patients with advanced melanoma. Patients and Methods Patients with known programmed death ligand 1 (PD-L1) expression and BRAF mutation status received nivolumab (NIVO) or dacarbazine in CheckMate 066 and either NIVO, ipilimumab (IPI), or NIVO+IPI in CheckMate 067. Whole exome sequencing and RNA sequencing were used to determine TMB and inflammatory gene expression signature scores, respectively. These biomarkers were evaluated in terms of their association with PFS and OS. Results In the NIVO, NIVO+IPI, and IPI arms of CheckMate 067, longer survival was associated with high (&amp;amp;gt; median) versus low (≤ median) TMB with hazard ratios (HRs) (95% confidence interval [CI]) for PFS of 0.45 (0.30–0.65), 0.55 (0.38–0.81), and 0.60 (0.43–0.82), and for OS of 0.46 (0.30–0.71), 0.53 (0.34–0.82), and 0.52 (0.36–0.74), respectively. For NIVO-treated patients, these results were confirmed in CheckMate 066. A survival benefit was observed with high TMB and absence of BRAF mutation. Survival was associated with high versus low inflammatory signature scores with HRs (95% CI) for PFS of 0.56 (0.34–0.94), 0.40 (0.23–0.72), and 0.43 (0.27–0.70), and for OS of 0.37 (0.20–0.66), 0.38 (0.19–0.74), and 0.46 (0.27–0.79), in the NIVO, NIVO+IPI, and IPI arms, respectively. Weak correlations were observed between PD-L1, TMB, and the inflammatory signature. Conclusions Combined assessment of TMB, inflammatory gene expression signature, and BRAF mutation status may be predictive for response to immunotherapy in advanced melanoma.

  Study EGAS00001004555

• Exploration of coding and non-coding variants in cancer using GenomePaint.

  GenomePaint (https://proteinpaint.stjude.org/genomepaint) is a dynamic visualization platform for whole-genome, whole-exome, transcriptome, and epigenomic data, featuring a novel design that captures the inter-relatedness between DNA variations and RNA expression. Regulatory non-coding variants can be inspected and discovered along with coding variants, and their functional impact further explored by examining 3D genome and/or ChIP-seq data generated from cancer cell lines. Further, GenomePaint correlates mutation and expression patterns with patient outcomes, and can display external data such as adult cancer datasets and user-provided custom tracks. We used GenomePaint to analyze multi-omics data from 3,652 pediatric cancers representing 16 histotypes, and demonstrate the visualization features through examples, including two that led to new insights into oncogenic mechanisms in pediatric cancer. The first is the discovery of a new class of pathogenic recurrent variants that cause aberrant splicing, disrupting the RING domain of CREBBP, a driver gene frequently mutated in relapsed pediatric leukemia. The second is the cis-activation of the MYC oncogene in a subset of B-lineage acute lymphoblastic leukemia (B-ALL) via duplication of the NOTCH1-MYC enhancer (N-ME), previously discovered only in T-lineage ALL. The regulatory impact of N-ME enhancer amplification was initially confirmed by allelic imbalance in published gene expression and ChIP-seq data and verified by additional Capture-C and fluorescence in situ hybridization data generated by follow-up experiments. These examples demonstrate the power of GenomePaint in enabling not only data visualization but also integrative genomic analysis that can lead to novel biological insight for follow-up experimental validation.

  Study EGAS00001004669

• Identification of the dismal subtype of B-ALL with dysregulation of CDX2 and UBTF

  PURPOSE: Advance in the knowledge of genomic basis of B-cell acute lymphoblastic leukemia (B-ALL) has changed the treatment approaches and diagnostic assays and improved their outcome. Despite the recent development of next generation sequencing approaches, there are still some cases that their drivers are unknown. We aimed to define and describe the characteristics of an additional B-ALL subtype by integrating genomic, transcriptomic, and epigenomic approach. METHODS: More than 2,000 B-ALL cases of RNA-seq data from from Eastern Cooperative Oncology Group (ECOG) and the American College of Radiology Imaging Network (ACRIN) (n=764), the Children's Oncology Group (COG), Munich Leukemia Laboratory (MLL) (n=280), and other several collaborators were analyzed to identify the new subtype. Whole genome sequence was performed to detect mutations, stractural variants, and copy number alterations. The tridimentional analysis of chromatin was performed with HiChIP. RESULTS: We identified the new high-risk B-ALL subtype “CDX2/UBTF” that exhibited unique gene expression profiles, which is enriched in female (77.3%), adolescents and young adults (52.6%), and relapse cohort (3.3%). CD10 negativity and IgM positivity are hallmarks of this subtype with NTRK3 expression that can be an effective marker at diagnostic approach and a potential targeted therapy with TRK and mulkinase inbitors, larotrectinib and entrectinib. Genetically, two different alterations define CDX2/UBTF. One is UBTF-ATXN7L3 fusions caused by deletion of 17q21.31, and the other is extopic CDX2 expression through the enhancer hijacking mechanism induced by deletion of 13q12.2. Other genomic features include gain of 1q and PAX5 rearrangement (PAX5-ZCCHC7), might induce upregulation of histone cluster genes and PAX5. CONCLUSION: We described the novel B-ALL subtype “CDX2/UBTF” that has unique clinical and genomic characteristics.

  Study EGAS00001005863

• Comprehensive characterization of pre- and post-treatment samples of breast cancer reveal potential mechanisms of chemotherapy resistance

  When locally advanced breast cancer is treated with neoadjuvant chemotherapy, the recurrence risk is significantly higher if no complete pathologic response is achieved. Identification of the underlying resistance mechanisms is essential to select treatments with maximal efficacy and minimal toxicity. Here we employed gene expression profiles derived from 317 HER2-negative treatment-naïve breast cancer biopsies of patients who underwent neoadjuvant chemotherapy, deep whole exome and RNA-sequencing profiles of 22 matched pre- and post-treatment tumors, and treatment outcome data to identify biomarkers of response and resistance mechanisms. Molecular profiling of treatment-naïve breast cancer samples revealed that expression levels of proliferation, immune response and extracellular matrix (ECM) organization combined predict response to chemotherapy. Triple negative patients with high proliferation, high immune response and low ECM expression had a significantly better treatment response and survival benefit (HR 0.29, 95% CI 0.10-0.85; p=0.02), while in ER+ patients the opposite was seen (HR 4.73, 95% CI 1.51-14.8; 0=0.008). The characterization of paired pre-and post-treatment samples revealed that aberrations of known cancer genes were either only present in the pre-treatment sample (CDKN1B) or in the post-treatment sample(TP53, APC, CTNNB1). Proliferation-associated genes were frequently down-regulated in post-treatment ER+ tumors, but not in triple negative tumors. Genes involved in ECM were upregulated in the majority of post-chemotherapy samples. Genomic and transcriptomic differences between pre- and post-chemotherapy samples are common and may reveal potential mechanisms of therapy resistance. Our results show a wide range of distinct, but related mechanisms, with a prominent role for proliferation- and ECM-related genes.

  Study EGAS00001005876

• Combination of CDK4/6 with BET-bromodomain and PI3K/mTOR inhibitors in medulloblastoma in vitro and in vivo

  Despite improvement in the treatment of medulloblastoma (MB) over the last years, numerous patients with MYC- and MYCN-driven tumors still fail the current therapeutic regimen. The retinoblastoma (RB) pathway is often compromised in MB with RB remaining intact suggesting that CDK4/6 inhibition might represent a therapeutic strategy for which drug combination remains understudied. Using tumor lines grown in stem cell conditions referred as tumorspheres as well as tumor cells freshly dissociated from the brain of mice harboring intracranial patient-derived xenografts (PDX) of Group3 (G3) and Sonic Hedgehog (SHH) MB, we conducted high throughput drug combination screens using the CDK4/6 inhibitor (CDK4/6i) ribociclib as an anchor and 87 oncology drugs approved by FDA or currently in clinical trials. Bromodomain and extra terminal (BET) and PI3K/mTOR inhibitors potentiated ribociclib’s inhibition of proliferation in G3 and SHH MB. A reverse combination screen using BET inhibitor JQ1 as the anchor, revealed the three FDA approved CDK4/6i as the most potentiating drugs. Those synergistic combinations were tested in vivo in mice intracranially implanted with G3 and SHH MB tumor cells. Ribociclib showed single agent activity in several in vivo MB models. However, despite in vitro synergy, combinations of ribociclib with JQ1 and ribociclib with the PI3K/mTOR inhibitor paxalisib failed to improve the survival of MB bearing mice compared to ribociclib alone. Our data illustrates the difficulty in translating in vitro findings to in vivo for brain tumors partly due to insufficient free brain concentration and altered signaling pathways activation in vitro compared to in vivo.

  Study EGAS00001006286

• Long-term temporal stability of peripheral blood DNA methylation alterations in patients with inflammatory bowel disease.

  Introduction There is great current interest in the potential application of DNA methylation alterations in peripheral blood as biomarkers of susceptibility, progression, and treatment response in inflammatory bowel disease (IBD). To date, the intra-individual stability of peripheral blood leukocyte (PBL) methylation in IBD patients, a key consideration in biomarker development, has not been characterised. Here, we studied the long-term stability of all probes located on the Illumina HumanMethylation EPIC BeadChip array. Methods We followed a cohort of 46 adult IBD patients (36 Crohn’s disease (CD), 10 ulcerative colitis (UC), median age 44 (IQR: 27-56), 50% female) that received standard care without any intervention at the Amsterdam UMC. Paired PBL samples were collected at two time points with a median 7 (range: 2-9) years in between. Differential methylation and intra-class correlation (ICC) analyses were used to identify time-associated differences and temporally stable CpGs, respectively. Results Around 60% of all CpG loci on the EPIC array presented poor (ICC &lt;0.50); 119,388 (≈14%) and 41,274 (≈5%) showed good (ICC ≥0.75) or excellent (ICC ≥0.90) intra-individual stability. Functional over-representation analyses of the stable methylated CpGs implicated genes involved in cell-cell signaling, adhesion and neurogenesis. Focusing on previously identified consistently differentially methylated positions indicated that 22 CD-, 11 UC-, and 24 IBD-associated loci demonstrated good stability (ICC ≥0.75) over time, where we observed a marked stability of CpG loci associated to the HLA genes. Conclusion Our data provide insight into the long-term stability of the PBL DNA methylome within an IBD context, facilitating the selection of biologically relevant and robust IBD-associated epigenetic biomarkers with increased potential for independent validation.

  Study EGAS00001006501

• BCR Signaling in human BM PC

  The functions of bone marrow plasma cells (BMPC) beyond antibody production are not fully elucidated and distinct subsets of BMPC suggest potential different functions. Phenotypic differences were identified for human BMPC depending on CD19 expression. Since CD19 is a co-stimulatory molecule of the B-cell-receptor (BCR), and IgA+ and IgM+ BMPC express the BCR on their surface, we here studied whether CD19 expression affects cellular responses, such as BCR signaling and the expression of checkpoint molecules. We analyzed 132 BM samples from individuals undergoing routine total hip arthroplasty. We found that both CD19+ and CD19- BMPC expressed BCR signaling molecules. Notably, the BCR-associated kinase spleen tyrosine kinase (SYK) including pSYK was higher expressed in CD19+BMPC compared to CD19- BMPC. BCR stimulation also resulted in increased kinase phosphorylation downstream of the BCR while expression of CD19 remained stable afterwards. Interestingly, the BCR response was restricted to IgA+ BMPC independently of CD19 expression. With regard to the expression of checkpoint molecules, CD19- BMPC expressed higher levels of co-inhibitory molecule programmed cell death protein-1 (PD-1) than CD19+ BMPC. IgA+ BMPC characteristically upregulated PD-1 upon BCR stimulation in contrast to other PC subsets and inhibition of the kinase SYK abrogated PD-1 upregulation. In contrast, expression of PD-1 ligand, B and T lymphocyte attenuator (BTLA) and CD28 did not change upon BCR activation of IgA+ BMPC. Here, we identify a distinct functiotype of IgA+ BMPC that is independent of the phenotypic heterogeneity of the subsets according to their CD19 expression. The data suggest that IgA+ BMPC underlie different regulatory principles and/or exert distinct regulatory functions.

  Study EGAS00001004948

• Integrative analysis of exome-seq, RNA-seq, ATAC-seq (bulk and single-cell), and Hi-C data generated from 3-D spatially mapped samples acquired during surgical resection from 10 patients diagnosed with IDH-WT glioblastoma

  Glioblastoma (GBM), the most common primary brain cancer in adults, remains incurable with no targeted therapies approved despite decades of investigation into its molecular landscape. Treatment failure is attributed to intratumoral heterogeneity in the GBM genome and epigenome, which foster tumor evolution and selection of resistant clones. However, tumor evolution and intratumoral heterogeneity remain poorly understood on the level of the whole tumor as most studies are based on single samples per patient and lack spatial context. Here, we have used 3-D neuronavigation during surgical resection for 10 primary IDH-WT GBM patients to collect 102 samples representing maximal tumor diversity, each mapped by 3-D spatial coordinates. We have applied a strategic set of genomic and epigenomic assays spanning multiple levels of resolution to discover, orthogonally validate, and functionally assess drivers of tumor evolution and intratumoral heterogeneity. These include extrachromosomal DNA amplifications, chromothripsis events, inversions, and translocations that disrupt both the GBM genome and epigenome while generating fusion transcripts and opportunities for therapeutic intervention. We define epigenomic programs that contribute to GBM evolution and intratumoral heterogeneity, revealing their 3-D spatial patterning within whole tumors and their cell type(s) of origin in single-cell data from the same tumor samples. Notably, we distinguish neuronal, glial, and immune programs aberrantly active in tumor cells from their counterparts in normal cells and discover NEUROD1, JUN/FOS, and NF1 transcription factors as key drivers of GBM evolution and growth. Collectively, these data provide unprecedented insight into GBM evolution and intratumoral heterogeneity from single-cell to whole-tumor resolution, redefining current understanding and providing a rich resource of targets for therapeutic investigation.

  Study EGAS00001006785

• Single_cell_measurements_to_characterise_B_cell_repopulation_in_SLE_after_rituximab_therapy__a_pilot_study

  Systemic lupus erythematosus (SLE) is a systemic autoimmune disease characterised by autoantibodies against DNA and nuclear proteins. Kidney involvement (lupus nephritis) is one of the most serious manifestations and is the leading determinant of long-term prognosis. Therapies targeting B cells (e.g. rituximab) are used in the treatment of severe lupus, including lupus nephritis. Rituximab is a monoclonal antibody that targets the CD20 surface protein on the surface of B cells. Following rituximab treatment, there is near complete depletion of B cells in peripheral blood. Peripheral blood B cells slowly return over a period of months. The speed of repopulation varies but typically B cells have returned to baseline levels between 9 and 12 months post rituximab therapy. Return of B cells may be associated with increased disease activity, but this is not always the case. For example, some patients remain in clinical remission for several years following one cycle of rituximab, despite B cell repopulation, whilst others relapse as soon as the B cell number reaches the normal range. The reasons for the different response remain unclear. The B cell subpopulations that return following B cell depleting therapy have been previously studied using flow cytometry but have never been characterised using single cell -omic technologies and thus no data on their molecular signatures are available. Previous studies have shown that naïve B cells were the predominant repopulating B-cells, indicating de novo production of B cells. Obtaining an integrated understanding of the molecular programs of these cells at the single cell level will enhance our ability to predict disease outcome, to identify patients with poorer B-cell drug responses, understand disease heterogeneity and design more effective B-cell targeted therapies approach.

  Study EGAS00001006798

• Analysis of T-cell receptor clonotypes in tumor micro-environment identifies shared cancer type-specific signatures

  Despite the conventional view that a truly random V(D)J recombination process should generate a highly diverse immune repertoire, emerging reports suggest that there is a certain bias towards the generation of shared/public immune receptor chains. These studies were performed in viral diseases where public T-cell receptors (TCR) appear to confer better protective responses. Selective pressures generating common TCR clonotypes are currently not well understood but it is believed that they confer a growth advantage. As very little is known about public TCR clonotypes in cancer, here we set out to determine the extent of shared TCR clonotypes in the intra-tumor microenvironments of virus- and non-virus- driven head and neck cancers using TCR sequencing. We report that tumor-infiltrating T-cell clonotypes were indeed shared across individuals with the same cancer type, where the majority of shared sequences were specific to the cancer type (ie viral versus non-viral). These shared clonotypes were not particularly enriched in EBV-associated nasopharynx cancer, but in both cancers, exhibited distinct characteristics, namely shorter CDR3 lengths, restricted V- and J-gene usages and also demonstrated convergent V(D)J recombination. Many of these shared TCRs were expressed in patients with a shared HLA background. Pattern recognition of CDR3 amino acid sequences revealed strong convergence to specific motifs that were unique to each cancer type, suggesting they may be enriched for specificity to common antigens found in the tumor microenvironment. The identification of shared TCRs in infiltrating tumor T-cells not only adds to our understanding of the tumor adaptive immune recognition but could also serve as disease-specific biomarkers and guide the development of future immunotherapies.

  Study EGAS00001005480