23014 results for "ega"

in 89.94 milliseconds.

• Coding and non-coding drivers of mantle cell lymphoma identified through exome and genome sequencing

  Mantle cell lymphoma (MCL) is an uncommon B-cell non-Hodgkin lymphoma (NHL) that is incurable with standard therapies. The genetic drivers of this cancer have not been firmly established and the features that contribute to differences in clinical course remain limited. To extend our understanding of the biological pathways involved in this malignancy, we performed a large-scale genomic analysis of MCL using data from 51 exomes and 34 genomes alongside previously published exome cohorts. To confirm our findings, we re-sequenced the genes identified in the exome cohort in 191 MCL tumors, each having clinical follow-up data. We confirmed the prognostic association of TP53 and NOTCH1 mutations. Our sequencing revealed novel recurrent non-coding mutations surrounding a single exon of the HNRNPH1 gene. In RNA-seq data from 103 of these cases, MCL tumors with these mutations had a distinct imbalance of HNRNPH1 isoforms. This altered splicing of HNRNPH1 was associated with inferior outcomes in MCL and showed a significant increase in protein expression by immunohistochemistry. We describe a functional role for these recurrent non-coding mutations in disrupting an auto-regulatory feedback mechanism, thereby deregulating HNRNPH1 protein expression. Taken together, these data strongly implicate a role for aberrant regulation of mRNA processing in MCL pathobiology.

  Study EGAS00001004289

• Patient-derived tumor organoids for personalized medicine in a rare case of hepatocellular carcinoma with neuroendocrine differentiation.

  We report the case of a 74-year-old man with a very rare subtype of hepatocellular carcinoma with neuroendocrine differentiation (HCC-NED). The patient presented with two independent tumors, a gastrointestinal stromal tumor in the stomach and a hepatocellular carcinoma in the liver. Both tumors were surgically removed in curative intent. Histopathological work-up of the liver tumor revealed poorly differentiated hepatocellular carcinoma (Edmondson-Steiner grade IV) with diffuse expression of neuroendocrine markers synaptophysin (SYP) and chromogranin (CHGA). Three months after resection, multifocal recurrence of the HCC with neuroendocrine differentiation (HCC-NED) was observed. In the meantime, tumor organoids have been generated from the resected HCC-NED and extensively characterized. Sensitivity to a number of drugs approved for the treatment of HCC or neuroendocrine carcinomas was tested in vitro. Based on their in vitro efficacy, etoposide and carboplatin were used as first line palliative combination treatment. Because genomic analysis revealed a NTRK1-mutation (kinase domain) and tumor organoids were sensitive to entrectinib, a pan-TRK inhibitor, the patient received entrectinib as second line therapy. After only two weeks, treatment had to be discontinued due to deterioration of the patient’s general condition. In conclusion, we demonstrate for the first time that preclinical drug testing using organoids is feasible in selected HCC cases.

  Study EGAS00001005887

• Transcriptomic analysis of Acute Myeloid Leukemia stem cells

  Acute Myeloid Leukemia (AML) is characterized by the accumulation of clonal myeloid blast cells unable to differentiate into mature leukocytes. Chemotherapy induces remission in the majority of patients, but relapse rates are high and lead to poor clinical outcomes. Since this is primarily caused by chemotherapy-resistant leukemic stem cells (LSCs), it is essential to eradicate LSCs to improve patient survival. LSCs have predominantly been studied at the transcript level, thus lacking information about post-transcriptionally regulated genes and associated networks. Here we extend our previous report on LSC proteomes to healthy age-matched hematopoietic stem and progenitor cells (HSPCs) and correlate the proteomes to the corresponding transcriptomes. By comparing LSCs to leukemic blasts and healthy HSPCs, we validate candidate LSC markers and highlight novel and potentially targetable proteins that are absent or only lowly expressed in HSPCs. In addition, our data provides strong evidence that LSCs harbor a characteristic energy metabolism, adhesion molecule composition, as well as RNA processing properties. Furthermore, correlating proteome and transcript data of the same individual samples highlights the strength of proteome analyses, which are particularly potent in detecting alterations in metabolic pathways. In summary, our study provides a comprehensive proteomic and transcriptomic characterization of functionally validated LSCs, blasts and healthy HSPCs, representing a valuable resource helping to design LSC-directed therapies.

  Study EGAS00001004402

• scRNA-seq of HGSC tumor and ascites

  Immune regulatory metabolites are key features of the tumor microenvironment (TME), yet with a few notable exceptions, their identities remain largely unknown. We uncovered the immune regulatory metabolic states and metabolomes of sorted tumor and stromal, CD4+, and CD8+ cells from the tumor and ascites of patients with high-grade serous ovarian cancer (HGSC) using high-dimensional flow cytometry and metabolomics supplemented with single cell RNA sequencing. Flow cytometry revealed that tumor cells show a consistently greater uptake of glucose than T cells, but similar mitochondrial activity. Cells within the ascites and tumor had pervasive metabolite differences, with a striking enrichment in 1-methylnicotinamide (MNA) in T cells infiltrating the tumor compared to ascites. Despite the elevated levels of MNA in T cells, the expression of nicotinamide N-methyltransferase, the gene encoding the enzyme that catalyses the transfer of a methyl group from S-adenosylmethionine to nicotinamide, was restricted to fibroblasts and tumor cells. Treatment of T cells with MNA resulted in an increase in T cell-mediated secretion of the tumor promoting cytokine tumor necrosis factor alpha. Thus, the TME-derived metabolite MNA contributes to an alternative and non-cell autonomous mechanism of immune modulation of T cells in HGSC. Collectively, uncovering the tumor-T cell metabolome may reveal metabolic vulnerabilities that can be exploited using T cell-based immunotherapies to treat human cancer.

  Study EGAS00001004829

• Genomics to select patients with metastatic breast cancer for targeted therapy (microarray_oncoscan)

  Cancer progression is driven in part by genomic alterations located in oncogenes or tumor suppressor genes. The genomic characterization of cancers has shown a large interpatient heterogeneity regarding the driver alterations, leading to the concept that generating a genomic profiling by multigene sequencing in patients with cancer could allow selecting effective targeted therapies. While this concept has been broadly implemented in daily practice, there is no evidence that such approach improves patient outcome, and how to optimally select the therapy to administer. A genomic profiling using next generation sequencing and copy number analyses was performed 1462 patients with Her2-non overexpressing metastatic breast cancer included in SAFIR02 Breast trial. 238 of these patients were randomized in two trials between a targeted therapy matched to genomic alteration and a maintenance chemotherapy. The trial shows that targeted therapies matched to genomics improves progression free survival (PFS) when genomic alterations are classified level I/II according to ESMO Scale for Clinical Actionability of Molecular Targets (ESCAT) (adjusted HR: 0.41, 90% CI: 0.27-0.61, p<0.001), but not when alterations are classified beyond level II (unadjusted HR: 1.15, 95% CI: 0.76-1.75). This trial provides evidence that the treatment decision led by genomics should be driven by a framework of target actionability in patients with mBC.

  Study EGAS00001005586

• Genomics to select patients with metastatic breast cancer for targeted therapy (microarray_agilent)

  Cancer progression is driven in part by genomic alterations located in oncogenes or tumor suppressor genes. The genomic characterization of cancers has shown a large interpatient heterogeneity regarding the driver alterations, leading to the concept that generating a genomic profiling by multigene sequencing in patients with cancer could allow selecting effective targeted therapies. While this concept has been broadly implemented in daily practice, there is no evidence that such approach improves patient outcome, and how to optimally select the therapy to administer. A genomic profiling using next generation sequencing and copy number analyses was performed 1462 patients with Her2-non overexpressing metastatic breast cancer included in SAFIR02 Breast trial. 238 of these patients were randomized in two trials between a targeted therapy matched to genomic alteration and a maintenance chemotherapy. The trial shows that targeted therapies matched to genomics improves progression free survival (PFS) when genomic alterations are classified level I/II according to ESMO Scale for Clinical Actionability of Molecular Targets (ESCAT) (adjusted HR: 0.41, 90% CI: 0.27-0.61, p<0.001), but not when alterations are classified beyond level II (unadjusted HR: 1.15, 95% CI: 0.76-1.75). This trial provides evidence that the treatment decision led by genomics should be driven by a framework of target actionability in patients with mBC.

  Study EGAS00001005587

• Targeted de-methylation of the FOXP3-TSDR

  CD4+ regulatory T cells (Tregs) are key mediators of immunological tolerance and promising effector cells for immuno-suppressive adoptive cellular therapy to fight autoimmunity and chronic inflammation. Their functional stability is critical for their clinical utility and has been correlated to the demethylated state of the TSDR/CNS2 enhancer element in the Treg lineage transcription factor FOXP3. However, proof for a causal contribution of the TSDR de-methylation to FOXP3 stability and Treg induction is so far lacking. We here established a powerful transient-transfection CRISPR-Cas9-based epigenetic-editing method for the selective de-methylation of the TSDR within the endogenous chromatin environment of a living cell. The induced de-methylated state was stable over weeks in clonal T cell proliferation cultures even after expression of the editing complex had ceased. Epigenetic editing of the TSDR resulted in FOXP3 expression, even in its physiological isoform distribution, proving a causal role for the de-methylated TSDR in FOXP3 regulation. However, successful FOXP3 induction was not associated with a switch towards a functional Treg phenotype, in contrast to what has been reported from FOXP3 overexpression approaches. Thus, TSDR de-methylation is required, but not sufficient for a stable Treg phenotype induction. Therefore, targeted demethylation of the TSDR may be a critical addition to published in vitro Treg induction protocols which so far lack FOXP3 stability.

  Study EGAS00001004867

• Single cell transcriptomics of hESC-derived midbrain dopaminergic neurons generated by a new human development-based protocol

  Recent advances in stem cell technologies provide new opportunities for regenerative medicine and Parkinson’s disease modelling. However, the precise composition and quality of human embryonic stem cell (hESC)-derived preparations at the single cell level, compared to endogenous brain standards, remains to be determined. In addition, multiple developmental pathways important for midbrain dopaminergic (mDA) neuron development remain to be implemented in hESCs. Here we use mouse developmental knowledge and human single cell RNA-sequencing (scRNA-seq) data to generate a new development-based protocol to differentiate hESCs into mDA neurons. We found that Laminin 511, dual WNT activation with CHIR99021 and WNT5A, together with FGF8b, improved different aspects of midbrain patterning. Moreover, mDA neurogenesis and differentiation were enhanced by activation of liver X receptors and inhibition of fibroblast growth factor signaling. Analysis of hESC-derived midbrain cultures by scRNA-seq revealed multiple hESC-derived clusters, which resembled those in the endogenous human ventral midbrain. Moreover, high quality mDA neurons were detected and found to progressively acquire functionality in vitro. Thus, we hereby provide a hESC differentiation protocol that recapitulates several aspects of endogenous mDA neuron development, and a strategy to determine cell composition and quality of hESC-derived preparations that can be useful for disease modeling and cell therapy.

  Study EGAS00001006313

• Intratumoural heterogeneity and immune modulation in lung adenocarcinoma of female smokers and never smokers

  Lung cancer is still the leading cause of cancer death worldwide despite declining smoking prevalence in industrialised countries. Although lung cancer is highly associated with smoking status, a significant proportion of lung cancer cases develop in patients who never smoked, with an observable bias towards female never smokers. A better understanding of lung cancer heterogeneity and immune system involvement during tumour evolution and progression in never smokers is therefore highly warranted. We employed single nucleus transcriptomics of surgical lung adenocarcinoma (LADC) and normal lung tissue samples from patients with or without smoking history. Immune cells as well as fibroblasts and endothelial cells respond to tobacco smoke exposure by inducing a highly inflammatory state in normal lung tissue. In the presence of LADC, we identified differentially expressed transcriptional programmes in macrophages and cancer-associated fibroblasts, providing insight into how the niche favours tumour progression. Within tumours, we distinguished eight subpopulations of neoplastic cells in female smokers and never smokers. Through pseudotemporal ordering, we inferred a trajectory towards two differentiated tumour cell states implicated in cancer progression and invasiveness. A proliferating cell population sustaining tumour growth exhibits differential immune modulating signatures in both patient groups. Our results resolve cellular heterogeneity and immune interactions in LADC, with a special emphasis on female never smokers and implications for the design of therapeutic approaches.

  Study EGAS00001006331

• Bulk-tissue RNA-sequencing paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from post-mortem neuropathologically-confirmed control individuals.

  Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

  Study EGAS00001006380

• 5- FU treated organoids

  5-Fluorouracil (5-FU) is a chemotherapeutic drug component that is commonly used for the treatment of solid cancers. The anticancer properties of 5-FU have been attributed to interference with nucleotide synthesis and through direct incorporation into the DNA as being a pyrimidine analog. As both mechanisms of action may have a mutational impact on surviving tumor cells, we performed three independent analyses to characterize the genomic consequences of 5-FU treatment: i) in vitro treatment of intestinal organoids with 5-FU followed by whole genome sequencing, ii) genome-wide mutation analyses in tumor samples from 5-FU treated patients with breast or colorectal cancer, and iii) analysis of paired biopsies from patients treated with 5-FU between biopsies. Our results demonstrate that both in vitro and in vivo 5-FU causes a mutational pattern that is dominated by T>G substitutions in a CTT context and with strong resemblance to COSMIC signature 17. Interestingly, this signature is also found in non 5-FU treated patients, most prominently in esophageal tumors, indicating that distinct endogenous and exogenous triggers can converge into highly similar mutational signatures. Furthermore, our results suggest that 5-FU may have adverse mutagenic effects on healthy cells and contributes to genetic variability in surviving cancer cells thereby contributing to tumor evolution.

  Study EGAS00001003592

• Real-time analysis of the cancer genome and fragmentome from plasma and urine short and long cell-free DNA using Nanopore sequencing

  Cell-free DNA (cfDNA) can be isolated from blood and/or urine of cancer patients and analyzed with sequencing. Unfortunately, most conventional short-read sequencing methods are technically challenging, labor intensive and time consuming, requiring several days but more typically weeks to obtain interpretable data which are limited by a bias for short cfDNA fragments. Here, we demonstrate that with Oxford Nanopore Technologies sequencing  we can achieve economical and ultra-fast delivery of clinical data from liquid biopsies.   Our ‘ITSFASTR’ approach is able to deliver copy number aberrations, and cfDNA fragmentation profiles in less than 24 hours from sample collection. The tumor-derived cfDNA fraction calculated from lung cancer patient plasma and urine from bladder cancer patients was highly correlated (R=0.98) to the tumor fraction calculated from conventional short-read sequencing of the same samples. cfDNA size profile and fragmentation patterns in plasma and urine exhibited the typical cfDNA features yet with a significantly higher proportion of informative fragments that exceed 300bp, exhibiting similar tumor fraction than shorter size fragments. Notably, comprehensive fragment-end composition and nucleosome profiling near transcription start sites can be determined from the same data.  We propose that ITSFASTR is the first point-of-care solution for obtaining genomic and fragmentomic results from liquid biopsies.

  Study EGAS00001006591

• Genomic insights into the pathogenesis of Epstein-Barr virus-associated diffuse large B-cell lymphoma by whole-genome and targeted amplicon-based sequencing

  Epstein–Barr virus (EBV)-associated diffuse large B-cell lymphoma not otherwise specified (DLBCL NOS) constitutes a distinct clinicopathological entity in the current World Health Organization (WHO) classification, yet its genomic features remain sparsely characterized. Combining whole-genome sequencing (WGS), targeted amplicon sequencing (tNGS) and fluorescence in-situ hybridization (FISH) of tumours from 47 cases of EBV+ DLBCL we delineate the genomic landscape of this rare disease. Interrogating WGS data by sequential gene-set enrichment and network-diffusion we describe an accumulation of mutations affecting the nuclear factor κB (NFκB) pathway (CSNK2A2, CARD10), IL-6/JAK/STAT (SOCS1/3, STAT3) and WNT signalling (FRAT1, SFRP5) alongside aberrations in immunological processes like interferon-response. Further, WGS uncovered structural aberrations including 6q deletions (5/8 patients), which were subsequently validated by FISH (14/32 cases). Integrated WGS and tNGS analysis emphasized a profile distinct from EBV-negative DLBCL enriched for mutations in ARID1A (45%), KMT2A/D (32/30%), ANKRD11 (32%), NOTCH2 (32%). Expanding on previous reports, we identified recurrent alterations in CCR6 (15%), DAPK1 (15%), TNFRSF21 (13%), CCR7 (11%) and YY1 (6%). Our comprehensive description of the mutational landscape in EBV+ DLBCL substantiates the evidence of its pathobiological independence. Further, these observations aid in the refinement of the molecular taxonomy for aggressive lymphomas and guides future therapeutic strategies.

  Study EGAS00001004941

• Genetic insights into the biological mechanisms governing human ovarian ageing

  Reproductive longevity is critical for fertility and impacts healthy ageing in women, yet insights into the underlying biological mechanisms and treatments to preserve it are limited. Here, we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in ~200,000 women of European ancestry. These common alleles influence clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations. Identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increase fertility and extend reproductive life in mice. Causal inference analyses using the identified genetic variants indicates that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases risks of hormone-sensitive cancers. These findings provide insight into the mechanisms governing ovarian ageing, when they act across the life-course, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

  Study EGAS00001004947

• Control of Focal Adhesion Kinase Activation by RUNX1-regulated miRNAs in high-risk AML

  We recently described a 16-gene expression signature for improved risk stratification of acute myeloid leukemia (AML) patients called the AML Prognostic Score (APS). A subset of APS high-risk AML patients showed increased levels of focal adhesion kinase (FAK), encoded by the Protein Tyrosine Kinase 2 (PTK2) gene, which was correlated with RUNX1 mutations. RUNX1 mutant cells are more sensitive to PTK2 inhibitors. As we were not able to detect RUNX1 binding sites in the PTK2 promoter, we hypothesized that RUNX1 might regulate micro(mi)RNAs that repress PTK2, such that loss-of-function RUNX1 mutations would result in reduced miRNA expression and derepression of PTK2. Examination of paired RNA-seq and miRNA-seq data from 301 AML cases revealed two miRNAs that positively correlated with RUNX1 expression, contained RUNX1 binding sites in their promoters and were predicted to target PTK2. We show that the hsa-let7a-2-3p and hsa-miR-135a-5p promoters are regulated by RUNX1, and that PTK2 is a direct target of both miRNAs. Even in the absence of RUNX1 mutations, hsa-let7a-2-3p and hsa-miR-135a-5p regulate PTK2 expression, and reduced expression of these two miRNAs sensitizes AML cells to PTK2 inhibition. These data explain how RUNX1 regulates PTK2, and identify potential miRNA biomarkers for targeting AML with PTK2 inhibitors.

  Study EGAS00001006491

• BinDel: software tool for detecting clinically significant microdeletions in low-coverage WGS-based NIPT samples

  Clinically pathogenic chromosomal microdeletions causing genetic disorders such as DiGeorge syndrome are rare genetic aberrations that can cause clinically relevant fetal and childhood developmental deficiencies. Clinical severity of such deficiencies depend on the exact genomic location and genes affected by the fetal chromosomal aberration. Here we present the BinDel, a novel region-aware microdeletion detection software package developed to infer clinically relevant microdeletion risk in low-coverage whole-genome sequencing NIPT data. To test BinDel, we quantified the impact of sequencing coverage, fetal DNA fraction, and region length on microdeletion risk detection accuracy. We also estimated BinDel accuracy on known microdeletion samples and clinically validated aneuploidy samples. BinDel identified each positive control sample as high risk. We also determined that it is critical to take into account that the sample with a detected high microdeletion risk does not have a full chromosome aneuploidy, as the latter can cause erroneous high microdeletion risk findings. We observed that lower sequencing coverage resulted in reduced microdeletion detection accuracy, and higher fetal fractions considerably increased the microdeletion detection accuracy, with coverage becoming increasingly relevant as fetal DNA fraction decreased. In conclusion, we developed an R package-based software tool BinDel for inferring fetal microdeletion risks, which accurately identified all positive control samples with microdeletion or -duplication aberrations as high-risk samples.

  Study EGAS00001006663

• Uncovering tumor intrinsic and extrinsic factors that regulate hepatocellular carcinoma growth using patient derived xenograft assays

  Several major factors limit our understanding of hepatocellular carcinoma (HCC). First, human HCCs are infrequently biopsied for diagnosis and thus are not often biologically interrogated. Second, HCC initiation and progression are strongly influenced by the cirrhotic microenvironment, and the exact contributions of intrinsic and extrinsic tumor factors are unclear. A powerful approach to examine the personalized biology of liver cancers and the influence of host tissues is with patient derived xenograft (PDX) models. In Asia, HCCs from hepatitis B virus patients have been efficiently converted into PDXs, but few parallel efforts from the west have been reported. In a large-scale analysis, we implanted 93 HCCs and 8 cholangiocarcinomas (CCAs) to systematically analyze host factors and to define an optimized platform for PDX development from both surgical and biopsy samples. NSG mice that had undergone partial hepatectomy (PHx) represented the best combination of engraftability, growth, and passageability, but overall rates were low and indicative of a unique intrinsic biology for HCCs in the US. PDX models preserved the histology and genetic features of parental tumors, and ultimately, 9 new models were usable for pre-clinical studies. Intriguingly, HCC PDXs were differentially sensitive to regorafenib and sorafenib and CCA PDXs were also highly sensitive to regorafenib. PDX models functionalize early and advanced stage HCCs and revealed unique biological features of liver cancers from the US.

  Study EGAS00001004020

• Cell-free DNA fragmentation patterns and personalized sequencing reveal circulating tumor DNA in urine and plasma of glioma patients

  Glioma-derived cell-free DNA (cfDNA) is challenging to detect using liquid biopsy as levels in body fluids are low. We determined the glioma-derived DNA fractions in tumor biopsies, in cerebrospinal fluid (CSF), plasma and urine samples, using deep sequencing of personalized capture panels. By sequencing cfDNA across thousands of mutations identified individually in each patient’s matched tumor we detected tumor-derived DNA in plasma (10/12) and urine samples (8/11). The median tumor fraction was 6.4x10-3 in CSF, 3.1x10-5 in plasma and 4.7x10-5 in urine. We identified a shift in the size distribution for mutant cfDNA fragments in these body fluids. Next, we analyzed cfDNA fragment sizes with paired-end shallow whole genome sequencing (WGS) in urine samples from 35 patients with gliomas, 8 individuals with non-malignant brain disorders, and 26 controls (n=69 individuals, 96 samples). cfDNA in urine of glioma patients was significantly more fragmented compared to urine from patients with non-malignant brain disorders (t-test, p=1.7x10-2) and compared to urine of controls (t-test, p=5.2x10-9). The proportion of DNA fragments <60 bp was higher in glioma patients urine and could be used for classification (AUC=0.93). Machine learning models integrating fragment lengths could identify urine samples from glioma patients (AUC=0.97 in cross-validation).

  Study EGAS00001004355

• Loss of LGR4/GPR48 causes severe neonatal salt-wasting due to disrupted WNT signaling altering adrenal zonation

  Disorders of isolated mineralocorticoid deficiency causing potentially life-threatening salt-wasting crisis early in life have been associated with gene variants of aldosterone biosynthesis or resistance, but in some patients no such variants are found. WNT/β-catenin signaling is crucial for differentiation and maintenance of the aldosterone producing adrenal zona glomerulosa (zG). We describe a highly consanguineous family with multiple perinatal deaths or infants presenting at birth with failure to thrive, severe salt-wasting crises associated with isolated hypoaldosteronism, nail anomalies, short stature, and deafness. Whole exome sequencing revealed a homozygous splice variant in the R-SPONDIN receptor LGR4 gene (c.618-1G>C) regulating WNT signaling. The resulting transcripts affected protein function and stability, and resulted in loss of Wnt/β-catenin signaling in vitro. The impact of LGR4 inactivation was analyzed by adrenal cortex specific ablation of Lgr4, using Lgr4Flox/Flox mated with Sf1:Cre mice. Inactivation of Lgr4 within the adrenal cortex in the mouse model caused decreased WNT signaling, aberrant zonation with deficient zG and reduced aldosterone production. Thus, human LGR4 mutations establish a direct link between LGR4 inactivation and decreased canonical WNT signaling with abnormal zG differentiation and endocrine function. Therefore, variants in WNT signaling and its regulators should systematically be considered in familial hyperreninemic hypoaldosteronism.

  Study EGAS00001006808

• Dynamics of sequence and structural cell-free DNA landscapes in small-cell lung cancer

  Patients with small-cell lung cancer (SCLC) have an exceptionally poor prognosis, calling for improved real-time non-invasive biomarkers of therapeutic response. We performed targeted error-correction sequencing on 171 serial plasmas and matched white blood cell (WBC) DNA from 33 patients with metastatic SCLC who received systemic treatment with chemotherapy or immunotherapy-containing regimens. Tumor-derived sequence alterations and plasma aneuploidy were evaluated serially and combined to assess changes in total cell-free tumor load (cfTL). Longitudinal dynamic changes in cfTL were monitored to determine circulating cell-free tumor DNA (ctDNA) molecular response during therapy. Combined tiered analyses of tumor-derived sequence alterations and plasma aneuploidy allowed for assessment of ctDNA molecular response in all patients. Molecular responses captured the therapeutic effect and long-term clinical outcomes in a more accurate and rapid manner compared to radiographic imaging. Patients with sustained molecular responses had longer overall (log-rank p=0.0006) and progression-free (log-rank p<0.0001) survival, with molecular responses detected on average 4 weeks earlier than imaging. ctDNA analyses provide a rapid approach for the assessment of early on-therapy molecular responses and have important implications for the management of patients with SCLC, including the development of improved strategies for real-time tumor burden monitoring.

  Study EGAS00001006831

• cfDNA dataset from extracelular vesicles and paired non-fractionated samples

  Liquid biopsies contain multiple analytes that can be mined to improve the detection and management of cancer. Beyond cell-free DNA (cfDNA), mutations have been detected in DNA associated with extracellular vesicles (EV-DNA). The genome-wide composition and structure of EV-DNA remains poorly characterized, and whether circulating EVs are enriched in tumor signal compared to unfractionated cfDNA is still unclear. Here, using whole genome sequencing from selected lung cancer patients with a relatively high cfDNA tumor content >5%, we determined that the tumor fraction and heterogeneity are comparable between DNA associated with small (<200 nm) EVs and matched plasma cfDNA. DNA in EV fractions, obtained with standardized size-exclusion chromatography, are comprised of short ~150-180 bp fragments and long >1000 bp fragments that are poor in tumor signal. Other fractions only exhibit short fragments with similar tumor DNA content. The composition in bases at the end of EV-DNA fragments, as well as their fragmentation patterns are similar to total plasma cfDNA. Mitochondrial DNA however is relatively enriched in EV fractions. Our results suggests that cfDNA in plasma may be of dual nature, either bound to proteins (including the nucleosome) but also associated to EVs. cfDNA associated to small EV (including exosomes) is however not preferentially enriched in tumor signal.

  Study EGAS00001006848

• CYP2C19 long-read sequencing

  Inter-individual differences in drug response based on genetic variations can lead to drug toxicity and treatment inefficacy. A large part of this variability is caused by genetic variants in pharmacogenes. Unfortunately, the Single Nucleotide Variant (SNV) panels currently used in clinical pharmacogenomic (PGx) testing are unable to detect all genetic variability in these genes. Long-read sequencing, on the other hand, has been shown to be able to resolve complex (pharmaco)genes. In this study we have assessed the added value of long-read sequencing for PGx focusing on the clinically important and highly polymorphic CYP2C19 gene. With a capture-based long-read sequencing panel we were able to characterize the entire region and assign variants to their allele of origin (phasing), resulting in the identification of 813 unique variants in 37 samples. To assess the clinical utility of this data we have compared the performance of three different *-allele tools (Aldy, PharmCat and PharmaKU) which are specifically designed to assign clinical haplotypes to pharmacogenes based on all input variants. We conclude that long-read sequencing can improve our ability to characterize the CYP2C19 locus, help to identify novel haplotypes and that *-allele tools are an useful asset in phenotype prediction. Ultimately, this approach could help to better predict an individual’s drug response and improve therapy outcomes.

  Study EGAS00001006929

• Genome wide variation in the Angolan Namib desert reveals unique Pre-Bantu ancestry

  The populations of the Angolan Namib Desert have been largely neglected in previous surveys of the genomic landscape of southern Africa. Although at present, the Namib is culturally dominated by Southwest Bantu-speaking cattle-herders, the region exhibits an extraordinary ethnographic diversity which includes an array of semi-nomadic peoples whose subsistence strategies fall outside the traditional division between foraging and food production and can thus be referred to as “peripatetic”. Among these small-scale populations are the last speakers of the Kwadi branch of the Khoe-Kwadi language family associated with the introduction of pastoralism into southern Africa (Kwepe), as well as a range of groups whose origins remain enigmatic (Kwisi, Twa and Tjimba). Using genome-wide data from 208 individuals belonging to nine ethnically diverse groups from the Angolan Namib and adjacent areas (Kwepe, Kwisi, Twa, Tjimba, !Xun, Kuvale, Himba, Nyaneka, Ovimbundu) in combination with published data from other regions of Africa, we reconstruct in detail the histories of contact emerging from pre-historic migrations to southern Africa and show that peripatetic groups from southwestern Angola stand out for exhibiting elevated levels of an unique, regionally-specific and highly divergent Pre-Bantu ancestry. These findings highlight the importance of the Namib for understanding the deep genetic structure of Africa.

  Study EGAS00001007011

• Germline whole genome sequencing of patients with Li-Fraumeni syndrome

  Five hundred ng to 1 ug of genomic DNA was submitted to The Centre for Applied Genomics (TCAG) at The Hospital for Sick Children for genomic library preparation and whole genome sequencing. DNA samples were quantified using the Qubit High Sensitivity Assay and purity was assessed using the Nanodrop OD 260/280 ratio. Approximately 500-700 ng of DNA was used as input material for library preparation using the Illumina TruSeq PCR-free DNA Library Prep Kit following the manufacturer’s recommended protocol. In brief, DNA was fragmented to 400 bp on average using sonication on a Covaris LE220 instrument. Fragmented DNA was end-repaired, A-tailed and indexed TruSeq Illumina adapters with overhang-T were added to the DNA. Libraries were then validated on a Bioanalyzer DNA High Sensitivity chip to check for size and absence of primer dimers and quantified by qPCR using Kapa Library Quantification Illumina/ABI Prism Kit protocol (KAPA Biosystems). Validated libraries were pooled in equimolar quantities and paired-end sequenced on an Illumina HiSeq X platform following Illumina’s recommended protocol to generate paired-end reads of 150-bases in length and an average depth of 40X. FASTQ files were aligned to the hg19 reference genome using BWA-mem v0.78. PCR duplicates were marked with Picard MarkDuplicates v1.1.08, and base recalibration and realignment was performed using GATK v2.8.1.

  Study EGAS00001007061

• BMP4 and temozolomide synergize in the majority of patient derived glioblastoma cultures

  The poor prognosis of patients with glioblastoma is partly due to resistance to standard treatment, which combines chemoradiation using temozolomide (TMZ). The concept of cancer stem cells provided new insight into tumor resistance and management, and isolated Glioblastoma stem-like cells (GSCs) which contribute to this therapy resistance. Bone morphogenetic protein 4 (BMP4) is a differentiation factor that stimulates astroglial differentiation of such GSCs, thereby reducing their self-renewal. Exposing GSCs to BMP4 could therefore sensitize them to chemotherapy with TMZ, but previous studies with various treatment paradigms have produced contradictory results in this matter. Therefore, we re-defined the optimal treatment paradigm for the combination of TMZ and BMP4 and subsequently re-assessed the inter-tumoral variability in response to co-treatment with TMZ+BMP4 in vitro. First, the simultaneous treatment was found more effective than sequential therapy. Second, after application of our optimized treatment protocol, 70% of a total of 20 patient-derived glioblastoma cultures showed synergy when the treatments were combined. The synergistic effect was a result of increased cellular apoptosis and not of decreased proliferation. Comparative RNA-sequencing suggested that combination treatment results in decreased MAPK signaling. We conclude that addition of BMP4 to standard therapy with TMZ could enlarge the group of patients with good therapeutic response.

  Study EGAS00001007095

• Structural variant analysis of homologous recombination-deficient genomes

  Homologous recombination (HR) deficiency causes DNA breaks and cytogenetic aberrations. Paradoxically, the types of DNA rearrangements specifically associated with HR-deficient cancers only minimally impact chromosomal structure. Addressing this, we combined a genome graph analysis of short-read whole genome sequencing (WGS) profiles across thousands of tumors with deep linked-read (LR) WGS of 46 BRCA1 or BRCA2 mutant breast cancers to discover a distinct class of HR deficiency-enriched rearrangements called reciprocal pairs. LR WGS showed that reciprocal pairs with identical rearrangement orientations gave rise to one of two distinct chromosomal outcomes, distinguishable only with long molecule data. While one (cis) outcome corresponded to the copy and pasting of a small segment to a distant site, a second (trans) outcome was a quasi-balanced translocation or multi-megabase inversion with substantial (10kb) duplications at each junction. The full spectrum of reciprocal pair outcomes could be explained by an HR-independent replication restart repair mechanism. LR WGS additionally identified single-strand annealing (SSA) as a BRCA2-deficiency specific repair pathway in human cancers. Replication restart- and SSA-associated SVs improved BRCA1- vs. BRCA2- deficiency classification and identified metastatic cancer cases with favorable chemotherapy responses. Our data reveal classes of BRCA1- and BRCA2-deficiency specific rearrangements as drivers of cytogenetic aberrations in HR deficient cells.

  Study EGAS00001007186

• deep-learning-powered tissue deconvolution for cfDNA

  Plasma cell-free DNA (cfDNA) is a noninvasive biomarker for cell death of all organs. Deciphering the tissue origin of cfDNA can reveal abnormal cell death because of diseases, which has great clinical potential in disease detection and monitoring. Despite the great promise, the sensitive and accurate quantification of tissue-derived cfDNA remains challenging to existing methods due to the limited characterization of tissue methylation and the reliance on unsupervised methods. To fully exploit the clinical potential of tissue-derived cfDNA, here we present one of the largest comprehensive and high-resolution methylation atlas based on 521 noncancer tissue samples spanning 29 major types of human tissues. We systematically identified fragment-level tissue-specific methylation patterns and extensively validated them in orthogonal datasets. Based on the rich tissue methylation atlas, we develop the first supervised tissue deconvolution approach, a deep-learning-powered model, cfSort, for sensitive and accurate tissue deconvolution in cfDNA. On the benchmarking data, cfSort showed superior sensitivity and accuracy compared to the existing methods. We further demonstrated the clinical utilities of cfSort with two potential applications: aiding disease diagnosis and monitoring treatment side effects. The tissue-derived cfDNA fraction estimated from cfSort reflected the clinical outcomes of the patients. In summary, the tissue methylation atlas and cfSort enhanced the performance of tissue deconvolution in cfDNA, thus facilitating cfDNA-based disease detection and longitudinal treatment monitoring.

  Study EGAS00001007213

• ATAC-Seq of CD4 T cell subsets

  Blood samples were taken from 4 age-matched healthy volunteers (two males and two females), and naïve CD4 T cells were isolated using the naïve CD4+ T Cell Isolation Kit II, human (Miltenyi Biotec). Two hundred thousand cells were plated in a 96-well plate, activated with Dynabeads Human T-activator CD3/CD28 (bead-to-cell ratio of 1:1, as recommended by the manufacturer) (ThermoFisher Scientific) and polarized towards CD4 T cell subtypes (Teff, Th1, Th2, rTh17, pTh17, and Treg cells) with specific cocktails of recombinant cytokines and neutralizing antibodies (see details in the material and methods) for each subset. At day 4 post-stimulation, forty thousand cells were used to perform ATAC-Seq on CD4 T cell subtypes. Again, crude nuclei of live CD4+ T cell subsets were treated with Tagment DNA buffer and Tagment DNA Enzyme (Nextera DNA Library Prep Kit, Illumina), and then the DNA was purified by MinElute PCR Purification Kit (Qiagen). Transposed DNA fragments were amplified using specific adapters followed by purification with MinElute PCR Purification Kit (Qiagen). Fragments from 240-360pb were selected in the PippinHT system (Sage Science). The quality of the library and its DNA concentration were assessed by Bioanalyzer instruments (Agilent Technologies) and ultimately submitted for sequencing using Illumina HiSeq 2500 sequencer, V4 chemistry.

  Study EGAS00001007345

• Towards standardized whole exome sequencing (WES) for cancer patients: lessons from a multicentric pilot study

  A growing number of druggable targets and national initiatives for precision oncology necessitate broader genomic profiling for many cancer patients. Whole exome sequencing (WES) offers unbiased analysis of the entire coding sequence, segmentation-based detection of copy number alterations (CNAs), and accurate determination of complex biomarkers including tumor mutational burden (TMB), homologous recombination repair deficiency (HRD), and microsatellite instability (MSI). To assess the inter-institution variability of clinical WES, we performed a comparative pilot study between German Centers of Personalized Medicine (ZPMs) from five participating institutions. Tumor and matched normal DNA from 30 patients were analyzed using custom sequencing protocols and bioinformatic pipelines. Calling of somatic variants was highly concordant with a sensitivity between 91-95% and a positive predictive value (PPV) of 82-95% compared with a three-institution consensus and full agreement for 16 of 17 druggable targets. Deviations were observed for cases with low VAF or coverage, differing annotations, and different filter protocols. CNAs showed overall agreement in 76% for the genomic sequence with high wet-lab variability. Complex biomarkers correlated strongly between institutions (HRD: 0.79-1, TMB: 0.97-0.99) and all institutions agreed on microsatellite instability. This study will contribute to the development of quality control frameworks for comprehensive genomic profiling and sheds light onto parameters that require stringent standardization.

  Study EGAS00001007363

• Machine Learning Signal Enrichment for Ultrasensitive Plasma Tumor Burden Monitoring

  In solid tumor oncology, circulating tumor DNA (ctDNA) is poised to transform care through accurate assessment of minimal residual disease (MRD) and therapeutic response monitoring. To overcome the sparsity of ctDNA fragments in low tumor fraction (TF) settings and increase MRD sensitivity, we previously leveraged genome-wide mutational integration through plasma whole genome sequencing (WGS). We now introduce MRD-EDGE, a composite machine learning-guided WGS ctDNA single nucleotide variant (SNV) and copy number variant (CNV) detection platform designed to increase signal enrichment. MRD-EDGESNV uses deep learning and a ctDNA-specific feature space to increase SNV signal-to-noise enrichment in WGS by 300X compared to previous WGS error suppression. MRD-EDGECNV also reduces the degree of aneuploidy needed for ultrasensitive CNV detection through WGS from 1 Gb to 200 Mb, thereby expanding its applicability to a wider range of solid tumors. We harness the improved performance to identify MRD following surgery in multiple cancer types, track changes in tumor burden in response to neoadjuvant immunotherapy in non-small cell lung cancer (NSCLC) and demonstrate ctDNA shedding in precancerous colorectal adenomas. Finally, the radical signal-to-noise enrichment in MRD-EDGESNV enables de novo mutation calling in melanoma and small cell lung cancer (SCLC) without matched tumor, yielding clinically informative TF monitoring for patients on immune checkpoint inhibition (ICI).

  Study EGAS00001007451

• Machine learning guided signal enrichment for ultrasensitive plasma tumor burden monitoring

  In solid tumor oncology, circulating tumor DNA (ctDNA) is poised to transform care through accurate assessment of minimal residual disease (MRD) and therapeutic response monitoring. To overcome the sparsity of ctDNA fragments in low tumor fraction (TF) settings and increase MRD sensitivity, we previously leveraged genome-wide mutational integration through plasma whole genome sequencing (WGS). We now introduce MRD-EDGE, a composite machine learning-guided WGS ctDNA single nucleotide variant (SNV) and copy number variant (CNV) detection platform designed to increase signal enrichment. MRD-EDGESNV uses deep learning and a ctDNA-specific feature space to increase SNV signal-to-noise enrichment in WGS by 300X compared to previous WGS error suppression. MRD-EDGECNV also reduces the degree of aneuploidy needed for ultrasensitive CNV detection through WGS from 1 Gb to 200 Mb, thereby expanding its applicability to a wider range of solid tumors. We harness the improved performance to identify MRD following surgery in multiple cancer types, track changes in tumor burden in response to neoadjuvant immunotherapy in non-small cell lung cancer (NSCLC) and demonstrate ctDNA shedding in precancerous colorectal adenomas. Finally, the radical signal-to-noise enrichment in MRD-EDGESNV enables de novo mutation calling in melanoma and small cell lung cancer (SCLC) without matched tumor, yielding clinically informative TF monitoring for patients on immune checkpoint inhibition (ICI).

  Study EGAS00001007545

• Cell culture differentiation and proliferation conditions influence the in vitro regeneration of the human airway epithelium

  The human airway mucociliary epithelium can be recapitulated in vitro using primary cells cultured in an Air-Liquid Interface (ALI), a reliable surrogate to perform pathophysiological studies. As tremendous variations exist between media used for ALI-cultured human airway epithelial cells, our study aimed to evaluate at a single-cell level the impact of several media (e.g. BEGMTM, PneumaCultTM, “Half&Half” and “Clancy”) on cell type distribution. Our work revealed the impact of these media on cell composition, gene expression profile, cell signaling or epithelial shape. We found higher proportions of multiciliated cells in PneumaCultTM-ALI and Half&Half, stronger EGF signaling from basal cells in BEGMTM-ALI, differential expression of the SARS-CoV-2 entry factor ACE2, and distinct secretomes transcripts depending on media used. We also established that proliferation in PneumaCultTM-Ex Plus favored secretory cells fate, showing the key influence of proliferation media on late differentiation epithelial characteristics. Altogether, our data offer a comprehensive repertoire for evaluating the effects of culture conditions on airway epithelial differentiation and will help to choose the most relevant medium according to the phenomenon to be investigated such as cilia or mucus biology or viral infection. We detail useful parameters that should be explored to document airway epithelial cell fate and shaping.

  Study EGAS00001007572

• Patient-derived tumoroids from CIC::DUX4 rearranged sarcoma identify MCL1 as a therapeutic target

  High-risk sarcomas, such as metastasized and relapsed Ewing and CIC-rearranged sarcomas, currently lack effective curative systemic treatment options. Precision medicine approaches offer hope and among them ex vivo drug response profiling of patient-derived tumor cells emerges as a promising novel tool to identify effective drugs for individual patients. Here, we established ex vivo culture conditions to propagate Ewing and CIC::DUX4-sarcomas as tumoroids in a co-clinical manner. These models retained their original molecular and functional characteristics, including recurrent ARID1A mutations, and serve as tumor avatars suitable for large-scale drug testing. Indeed, screening a large drug library on a small living biobank of such tumors not only revealed distinct response differences between the two entities, but also identified an exceptional dependency of CIC::DUX4-sarcoma cells on MCL1. Mechanistically. MCL1 was identified as a direct target gene of the CIC::DUX4 fusion protein. Genetic and pharmacological inhibition of MCL1 induced rapid apoptosis in CIC::DUX4-sarcoma cells and inhibited tumor growth in a xenograft model. Thus, MCL1 represents a potential novel therapeutic target for CIC::DUX4-sarcoma. Overall, our study highlights the feasibility of co-clinical drug response profiling for individual sarcoma cases and suggests that prospective tests, to evaluate its clinical benefit are highly warranted.

  Study EGAS00001008039

• This study generated Oxford Nanopore long read sequencing data of cancer cell line mixtures for validating long-read variant calls in cancer genomics

  Long-read sequencing (LRS) improves genome alignment and facilitates resolving variants in complex genomic regions, making it a promising approach for cancer variant detection and biomarker discovery. Here, we evaluate the performance of LRS to detect somatic variants across different tumour purities and sequencing depths by comparing to somatic variants from short-read sequencing. We generated experimental mixtures of cancer cell lines and matched normal cell lines to simulate 10 tumour purities (ranging from 0 to 100% tumour content). This resulted in 22 samples which were sequenced using whole genome LRS to a targeted depth of 60x and somatic variants were identified. We downsampled sequencing data to explore the optimal read depth for somatic variant detection. Our results show that long-read variant calling tools achieve recall rates comparable to short-read gold standards. While a tumour sequencing depth of 30x to 60× is generally sufficient for detecting common variants, particularly structural variants (SVs), sequencing the matched normal sample at adequate depth is crucial for accuracy. Notably, we found variants unique to LRS that may represent real and previously undetected events. This study highlights key factors for optimising cancer genome sequencing with LRS and provides a comprehensive dataset of cell line mixtures for future research and tool development.

  Study EGAS00001008107

• WES data (cfTrack study) from NSCLC patients and OC patients

  We also collected samples from 8 NSCLC patients and 4 ovarian cancer patients and. For all 8 NSCLC patients, a tumor biopsy sample, a WBC sample, and three plasma samples were collected. For all 4 ovarian cancer patients, a WBC sample and two serum samples were collected. We collected tumor tissue sample from one ovarian cancer patient (OV4). The cfDNA was extracted from their plasma samples using the QIAamp circulating nucleic acid kit from QIAGEN (Germantown, MD). For serum cfDNA, ampure XP beads size selection was further performed to eliminate gDNA contamination. In brief, 0.5 volume of beads were first added to the cfDNA samples. After incubation, the supernatant was transferred to a new tube and an additional 2.0 volume of beads were added. After 80% ethanol wash, cfDNA was eluted from the beads. FA assays (Agilent Technologies) were performed to rule out the contamination of gDNA in the size selected samples. The cfDNA WES library of all patients and the genomic DNA WES library of the 4 ovarian cancer patients were constructed with the SureSelect XT HS kit from Agilent Technologies (Santa Clara, CA) according to the manufacturer’s protocol. In brief, 10ng of cfDNA was used as input material. After end repair/dA-tailing of cfDNA, the adaptor was ligated. The ligation product was purified with Ampure XP beads (Beckman-Coulter, Atlanta, GA) and the adaptor-ligated library was amplified with index primer in 10-cycle PCR. The amplified library was purified again with Ampure XP beads, and the amount of amplified DNA was measured using the Qubit 1xdsDNA HS assay kit (ThermoFisher, Waltham, MA). 700-1000 ng of DNA sample was hybridized to the capture library and pulled down by streptavidin-coated beads. After washing the beads, the DNA library captured on the beads was re-amplified with 10-cycle PCR. The final libraries were purified by Ampure XP beads. The library concentration was measured by Qubit, and the quality was further examined with Agilent Bioanalyzer before the final step of 2x150bp paired-end sequencing at an average coverage of 200. Whole-exome capture libraries of genomic DNA from the 8 NSCLC patients were constructed via Roche SeqCap EZ Exome V6 (Roche). Enriched exome libraries were sequenced on the Illumina HiSeq 3000 platform (Illumina) to generate 2x100bp paired-end reads at an average coverage of 200.

  Dataset EGAD00001008454

• Single cell transcriptional characterization of human megakaryocyte lineage commitment and maturation

  Single cell RNA-sequencing of sternal bone marrow reciding Hematopoietic Stem Cells (HSCs) and Megakaryocytes (MKs) from individuals undergoing elective open heart valve replacement. HSCs were defined as Lineage-, CD34+, CD38-, CD45RA-, CD90+, CD49f+ cells. MKs where CD41a+, CD42b+ and ploidy was determined with Hoechst. A sternal bone marrow scraping was taken directly following median sternotomy using a Volkmann’s spoon. The sample was collected into an EDTA Vacutainer tube containing 1.8mg/ml EDTA. 4mL of Dulbecco’s phosphate buffered saline (PBS, Sigma) containing 10% human serum albumin (HSA, Gemini Bio Products) was added and the whole volume was resuspended by pipetting 2-3 times. The sample was then put on metallic thermal beads (ThermoFisher Scientific) at a temperature between 0-4°C and transported to the University of Cambridge for further processing. For HSC isolation the cells were stained with the following antibody cocktail: PECy5 conjugated anti-lineage specific antibodies: CD2 (BD), CD3 (BD), CD10 (BD), CD11b (BD), CD11c (BD), CD19 (BD), CD20 (BD), CD56 (BD), biotinylated CD42b (Pab5, NHS Blood and Transplant, International Blood Group Reference Laboratory [IBGRL]), biotinylated GP6 (Pab5, NHS Blood and Transplant, International Blood Group Reference Laboratory [IBGRL]) used in combination with PECy5 conjugated streptavidin (Biolegend). Alexa Fluor 700 conjugated anti-CD34 (BD), PerCP-Cy5.5 conjugated anti-CD38 (BD), Pacific Blue conjugated anti-CD45RA (Invitrogen), PECy7 conjugated anti-CD90 (BD),PE conjugated anti-CD49f (BD). After staining cells were kept at 4°C before sorting using a FACS Aria Fusion flow sorter (BD). Single HSCs defined as Lineage-, CD34+, CD38-, CD45RA-, CD90+, CD49f+ cells were sorted by FACS directly into individual wells of a 96-well plate. Index sort data was collected for each single cell. For MK isolation the cells were stained for surface MK markers with mouse anti-human CD41a APC conjugated antibody (BD) and mouse anti-human CD42b PE conjugated antibody (BD) and for ploidy analysis with 1ug/ml Hoechst 33342 (Invitrogen). After incubation at 37°C for 30 minutes, the cells were kept at 4°C before sorting using a FACS Aria Fusion flow sorter (BD). Single cells and MK pools of 20 cells were sorted by FACS according to ploidy level using a 100uM nozzle directly into individual wells of a 96-well plate. cDNA synthesis and poly(A) enrichment was performed following the G&T-seq protocol (Macaulay et al. 2015), a variation of the Smart-seq2 protocol1. ERCC spike-in RNA (Ambion) was added to the lysis buffer in a dilution of 1:4,000,000.

  Dataset EGAD00001006677

• ATAC-seq of a selected group of AML cases

  ATAC-seq data were generated for a number of selected patients to investigate changes in enhancer and promoter regions. ATAC-seq was essentially carried out as described in31. Briefly, prior to transposition the viability of the cells was assessed and 1 × 106 cells were treated in culture medium with DNase I (Sigma) at a final concentration of 200 U ml−1 for 30 minutes at 37 °C. After Dnase I treatment, cells were washed twice with ice-cold PBS, and cell viability and the corresponding cell count were assessed. 5 × 104 cells were aliquoted into a new tube and spun down at 500 × g for 5 minutes at 4 °C, before the supernatant was discarded completely. The cell pellet was resuspended in 50 µl of ATAC-RSB buffer (10 mM Tris-HCl pH 7.4, 10 mM NaCl, 3 mM MgCl2) containing 0.1% NP-40, 0.1% Tween-20, and 1% Digitonin (Promega), and was incubated on ice for 3 minutes to lyse the cells. Lysis was washed out with 1 ml of ATAC-RSB buffer containing 0.1% Tween-20. Nuclei were pelleted at 500 × g for 10 minutes at 4 °C. The supernatant was discarded carefully and the cell pellet was resuspended in 50 µl of transposition mixture (25 µl 2× tagment DNA buffer, 2.5 µl transposase (100 nM final; Illumina), 16.5 µl PBS, 0.5 µl 1% digitonin, 0.5 µl 10% Tween-20, 5 µl H2O) by pipetting up and down six times. The reaction was incubated at 37 °C for 30 minutes with mixing before the DNA was purified using the Monarch PCR & DNA Cleanup Kit (NEB) according to the manufacturer’s instructions. Purified DNA was eluted in 20 µl elution buffer (EB) and 10 µl purified sample was objected to a ten-cycle PCR amplification using Nextera i7- and i5-index primers (Illumina). Purification and size selection of the amplified DNA were carried out with Agencourt AMPure XP beads. For purification the ratio of sample to beads was set to 1:1.8, whereas for size selection the ratio was set to 1:0.55. Purified samples were eluted in 15 µl of EB. Quality and concentration of the generated ATAC libraries were analyzed using the High Sensitivity D1000 ScreenTape Kit (Agilent) and libraries were sequenced paired-end on a NovaSeq (Illumina). ATAC-seq reads were aligned to the human reference genome build hg19 with bowtie2

  Dataset EGAD00001007583

• RRBS data (cfSort study) from noncancer tissue DNA

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

  Dataset EGAD00001010880

• RNA-sequencing of human skin samples obtained from SSc patients and healthy controls.

  Patient information SSc patients were recruited at the Department of Rheumatology of the Leiden University Medical Center (Leiden, The Netherlands). All patients met the American Rheumatism Association classification criteria for SSc (Subcommittee for scleroderma criteria 1980), and were classified according to LeRoy and Medsger criteria as either limited or diffuse cutaneous disease (LeRoy EC, Black C, Fleischmajer R, Jablonska S, Krieg T, Medsger TA Jr, Rowell N 1988). Institutional review board approval and written informed consent was obtained before patients entered this study. Two 4 mm skin biopsies were taken from a standardized location on the most proximal part of the lower arm, distal from the elbow. In 10 patients the skin biopsy came from a clinically affected area and in 4 patients the skin was locally unaffected. One sample was used for RNA sequencing and one sample was used for immunohistochemistry. Skin biopsies from healthy individuals were commercially sourced (Tissue Solutions, UK) and collected from donors undergoing skin resection surgery and after informed consent. To match the healthy skin with patients as much as possible, skin biopsies from healthy controls were also taken from a similar position (the under-arm (for 4 controls) and leg (for 2 controls)). Healthy skin donors were selected to match the age and sex of the SSc patient cohort. Biopsies from patients and controls were equally treated and were both stored at -80°C until RNA isolation was performed. RNA from frozen skin biopsies was isolated using RNeasy kit from fibrous tissue (Qiagen, the Netherlands). RNA quantity was determined by using SimplyNano 2000 and quality was assessed on Tapestation (Agilent, the Netherlands). All samples included in the study had a RIN score above 7.0. Transcriptome characterisation and analysis RNA sequencing was performed using polyA selection and a stranded protocol using Ion Torrent next generation sequencing technology (Service XS, The Netherlands). The Ion PI Template OT2 200 Kit v3 and Ion PI Sequencing 200 Kit v3 were used according to the manufacturer’s instructions. 20 samples were run on 11 PI chips. PI chip analyses, base calling and quality checks were performed using the Torrent Server Suite. An average of 42 million 100 bp reads was generated per sample. Following quality control, reads were aligned to the human genome (Homo sapiens GRh38.78) using Bowtie2 and STAR (Dobin et al. 2013; Langmead and Salzberg 2012). Reads were first aligned with STAR. For the unmapped reads from STAR, a second alignment step was performed using bowtie2 (local very sensitive options)

  Dataset EGAD00001003832

• Whole Genome Sequencing of clonal expansions of single healthy somatic cells (human, female)

  HSC73_clone: Bone marrow mononuclear cells from the healthy 73 years old female were thawed and labeled with Alexa-Fluor 488-conjugated anti-CD34 (581, Biolegend), Alexa-Fluor 700-conjugated anti-CD38 (HIT2, eBioscience), a cocktail of APC-conjugated lineage antibodies consisting of anti-CD4 (RPA-T4), anti-CD8 (RPA-T8), anti-CD11b (ICRF44), anti-CD20 (2H7), anti-CD56 (B159, all BD Biosciences), anti-CD14 (61D3), anti-CD19 (HIB19) and anti-CD235a (HIR2, all eBiocience) and 1 micro-gram/ml propidium iodide (Sigma). Using a BD FACSAria cell sorter, single Lin-CD34+CD38-PI- cells were individually sorted into low-adhesion 96-well tissue culture plates (Corning) containing 100micro-litre of StemSpan Serum-Free Expansion Medium (Stemcell technologies) supplemented with 100ng/ml of human SCF and FLT-3L, 50ng/ml of human TPO, 20ng/ml of human IL-3, IL-6 and G-CSF (all cytokines from Peprotech) and 50U/ml of penicillin and 50μg/ml of streptomycin (Sigma). Cells were incubated at 37 degrees C in a humidified atmosphere with 5% CO2 in air. After 5 days in culture, another 100micro litres of cytokine-containing medium were added. 13 days after seeding, clones B6 and G2 had expanded to approx. 105 cells and were selected for whole genome sequencing (2x101bp, paired-end, Illumina HiSeq2500) after tagmentation-based library preparation (see Extended Experimental Procedures) for clone B6 and standard library preparation for clone G2. For germline-control ~106 unsorted bone marrow mononuclear cells from the same donor were used for sequencing. An average of 30-fold sequence coverage for each the clones and the matching control were obtained.L4clone: A progenitor cell clone was raised from a peripheral blood sample of the 39 year old healthy female. Frozen peripheral blood mononuclear cells (PBMCs) were isolated from 2 ml heparinised peripheral blood via Ficoll Paque density centrifugation. A methylcellulose assay was performed as described earlier (Weisse et al., 2012). In brief, non-adherent mononuclear cells were incubated in the presence of the recombinant human cytokines IL-3, IL-5 and GM-CSF (R&D systems) over 14 days to induce colony formation. Colonies were detected under an inverted light microscope, and plucked by a pipette when colonies had approximately 10,000 cells/CFU. Each colony was washed three times in PBS and finally frozen as a cell pellet in -80 degrees C. Genomic DNA was isolated using the QIAamp DNA micro kit according to the instructions of the manufacturer (Qiagen, Hilden, Germany). Whole genome sequencing (2x101bp, paired-end, Illumina HiSeq2500) was performed for colony 4 after tagmentation-based library preparation and resulted in 15-fold sequence coverage for each the colony and the matching whole blood.

  Dataset EGAD00001000666

• The Dutch Microbiome Project Batch 1

  The Dutch Microbiome Project (DMP) data includes shotgun metagenomic sequencing of faecal samples 8,208 Dutch individuals. Paired-end sequencing was performed using Illumina HiSeq 2000 platform. Data is archived in two batches to facilitate easier data access and upload to EGA. Batch 1 of DMP includes 4396 samples.

  Dataset EGAD00001007027

• The patterns and dynamics of genomic instability in metastatic pancreatic cancer

  Pancreatic cancer is an aggressive malignancy with a five-year mortality of 97–98%, usually due to widespread metastatic disease. Previous studies indicate that this disease has a complex genomic landscape, with frequent copy number changes and point mutations, but genomic rearrangements have not been characterized in detail. Despite the clinical importance of metastasis, there remain fundamental questions about the clonal structures of metastatic tumours, including phylogenetic relationships among metastases, the scale of ongoing parallel evolution in metastatic and primary sites, and how the tumour disseminates. Here we harness advances in DNA sequencing to annotate genomic rearrangements in 13 patients with pancreatic cancer and explore clonal relationships among metastases. We find that pancreatic cancer acquires rearrangements indicative of telomere dysfunction and abnormal cell-cycle control, namely dysregulated G1-to-S-phase transition with intact G2–M checkpoint. These initiate amplification of cancer genes and occur predominantly in early cancer development rather than the later stages of the disease. Genomic instability frequently persists after cancer dissemination, resulting in ongoing, parallel and even convergent evolution among different metastases. We find evidence that there is genetic heterogeneity among metastasis-initiating cells, that seeding metastasis may require driver mutations beyond those required for primary tumours, and that phylogenetic trees across metastases show organ-specific branches. These data attest to the richness of genetic variation in cancer, brought about by the tandem forces of genomic instability and evolutionary selection.

  Study EGAS00000000064

• Nanopore sequencing enables allelic phasing of FLG loss-of-function variants, intragenic copy number variation and methylation status in atopic dermatitis and ichthyosis vulgaris

  Loss-of-function (LoF) variants in the FLG gene are associated with ichthyosis vulgaris (IV) and atopic dermatitis (AD). IV and AD patients with two FLG LoF variants are reported to have a more severe disease course. However, subsets of compound heterozygous patients display phenotypic heterogeneity despite a prediction of a severe phenotype from their genotype. The underlying genetic mechanism of this phenotypic variability is unknown. One possibility is that two LoF FLG variants could be located in cis (on the same allele), thus altering the predicted gene dosage. To investigate the allelic contribution of compound heterozygous FLG LoF mutations to IV phenotypic variability, we employ long reads from Oxford Nanopore Technologies, utilising two methodologies to phase the entire coding region of FLG into parental alleles: an amplicon-sequencing approach and an adaptive-sampling strategy. We investigated the allelic phasing of FLG LoF variants and intragenic-repeat copy-number variations (CNV) of 21 subjects with known compound heterozygous genotypes and mixed IV severity. Our results suggest that IV severity is not linked to the presence of cis variants in our IV patient cohort. Nanopore reads enabled detection and accurate phasing of CNVs and LoF variants into parental alleles, providing a framework for future in-depth genetic analysis. These results highlight the usefulness of utilising long reads for interrogating the FLG gene and other highly repetitive genes for genetic stratification related to clinical phenotypes.

  Study EGAS50000000166

• The tissue origin of highly elevated cell-free DNA in patients with and without cancer

  Forty years ago, the first measurements of cell free DNA (cfDNA) from cancer patients revealed that its total concentration was often elevated compared to healthy controls. The basis for this increase, and in particular the source of the extra cfDNA in such patients, has never been determined. To address this issue, we assessed methylation patterns of cfDNA from 57 individuals, 34 of whom had unequivocally elevated cfDNA concentrations (i.e., >15 ng/mL of plasma). We found that the excess cfDNA in cancer patients did not come from either neoplastic cells or surrounding cells from the tumor's tissue of origin. Instead, 80% of the cfDNA was derived from neutrophils, B cells or T cells, with neutrophils accounting for ~2/3. B cells contributed more than T cells despite the great excess of T cells over B cells in the blood. To see if this result was generalizable to other clinical situations in which elevated cfDNA concentrations have been described, we evaluated nine patients before and one day following surgery. After surgical resection, cfDNA concentration increased by an average of 16.1-fold, as expected. The majority of this increase (55%) was contributed by leukocytes, though damaged hepatocytes contributed a substantial minority (31% of increase). These data suggest that the elevated concentrations of cfDNA in cancer patients, as with patients after surgery, arises largely from dying leukocytes—not from the neoplastic cells or the tissue of origin of the cancer.

  Study EGAS00001005400

• METABRIC miRNA landscape

  MicroRNAs (miRNAs) are key players in cancer pathogenesis1. They display differential expression across breast cancer subtypes, and play oncogenic and tumor-suppressive roles2-6. A systematic study of the global miRNA landscape in breast cancer and the factors shaping it has not been undertaken. Here we report the results of profiling miRNA expression in 1,302 breast tumors with detailed clinical annotation and long-term follow-up, and for which matching genomic and mRNA expression data were available7. This provides a comprehensive view of the quantity, distribution and variation of the miRNA population and dissects how much genomic, transcriptional and post-transcriptional events contribute to miRNA expression architecture. The key clinical parameters and cellular pathways related to the miRNA landscape are identified, exemplified by context-dependent interaction with cell adhesion and Wnt signaling. Strikingly, in breast cancer miRNAs appear to act as modulators of mRNA-mRNA interactions rather than molecular switches. We demonstrate an important modulatory role for miRNAs in the biology of breast tumors devoid of somatic copy number aberrations, a common subtype where immune response is prominent. Remarkably, prognostic miRNA signatures derived in this subtype are consistently prognostic across several other subtypes and validate in external cohorts. These findings represent a new framework to study the biology of miRNA in human cancer.

  Study EGAS00000000122

• Machine learning to detect the SINEs of cancer

  We previously described an approach called RealSeqS to evaluate aneuploidy in plasma cell-free DNA (cfDNA) through the amplification of ~350,000 repeated elements with a single primer. We hypothesized that an unbiased evaluation of the large amount of sequencing data obtained with RealSeqS might reveal other differences between plasma samples from patients with and without cancer. This hypothesis was tested through the development of a novel machine-learning approach called Alu Profile Learning Using Sequencing (A-PLUS) and its application to samples from 5108 individuals, 2037 with cancer and the remainder without cancer. Samples from cancer patients and controls were pre-specified into four cohorts used for: 1) model training, 2) analyte integration and threshold determination, 3) validation, and 4) reproducibility. A-PLUS alone provided a sensitivity of 40.5% across 11 different cancer types in the Validation Cohort, at a specificity of 98.5%. Combining A-PLUS with aneuploidy and 8 common protein biomarkers detected 51% of 1167 cancers at 98.9% specificity. We found that part of the power of A-PLUS could be ascribed to a single feature – the global reduction of AluS sub-family elements in the circulating DNA of cancer patients. We confirmed this reduction through the analysis of another independent dataset obtained with a very different approach (whole genome sequencing). The evaluation of Alu elements therefore has the potential to enhance the performance of several methods designed for the earlier detection of cancer.

  Study EGAS00001007169

• Neoadjuvant immune checkpoint blockade in women with mismatch repair deficient endometrial cancer

  Neoadjuvant immune checkpoint blockade (ICB) has shown unprecedented activity in mismatch repair deficient (MMRd) colorectal cancers, but its effectiveness in MMRd endometrial cancer (EC) remains unknown. In this investigator-driven, phase I, feasibility study (NCT04262089), 10 women with MMRd EC of any grade, planned for primary surgery, received two cycles of neoadjuvant pembrolizumab (200 mg IV) every three weeks. A pathologic response (primary objective) was observed in 5/10 patients, with 2 patients showing a major pathologic response. No patient achieved a complete pathologic response. A partial radiologic response (secondary objective) was observed in 3/10 patients, 5/10 patients had stable disease and 2/10 patients were non-evaluable on Magnetic Resonance Imaging (MRI). All patients completed treatment without severe toxicity (exploratory objective). At median duration of follow-up of 22.5 months, two non-responders experienced disease recurrence. In-depth analysis of the loco-regional and systemic immune response (predefined exploratory objective) showed that monoclonal T cell expansion significantly correlated with treatment response. Tumour-draining lymph nodes displayed clonal overlap with intra-tumoural T cell expansion. All pre-specified endpoints, efficacy in terms of pathologic response as primary endpoint, radiologic response as secondary outcome and safety and tolerability as exploratory endpoint, were reached. Neoadjuvant ICB with pembrolizumab proved safe and induced pathologic, radiologic, and immunologic responses in MMRd EC, warranting further exploration of extended neoadjuvant treatment.

  Study EGAS50000000483

• Biopsy-Derived Organoids in Personalised Early Breast Cancer Care: Challenges of Tumour Purity and Normal Cell Overgrowth Cap Their Practical Utility

  The ability to establish organoids composed exclusively of tumour rather than healthy cells is essential for their implementation into clinical practice. Organoids have in recent years emerged as a powerful tool to expand patient material in culture and generate modifiable 3D models derived from humans or animal models. For translational research, they enable the creation of model systems for an ever-increasing number of cell types and diseases. And in personalised medicine, they potentially allow for functional drug testing with high predictive power in certain settings. We found that using biopsy material from untreated, early-stage primary breast cancer patients poses significant challenges for consistently culturing tumour cells as organoids. Specifically, we observed frequent outgrowth of genetically normal, non-cancerous epithelial cells. We analysed >100 biopsy samples from early-stage breast cancer and present our large collection of >70 organoid lines. We also show methods of assessing successful tumour cell culture in a time and cost-efficient manner, proving a high rate (>85%) of normal cell overgrowth in early-stage breast cancer organoids. Finally, we show a number of successful attempts to culture cancer organoids from mastectomy-derived tissue of advanced, metastatic breast cancer. We conclude that the usefulness of organoids from early breast cancer for translational research and personalised medicine, especially guidance of adjuvant or post-surgical maintenance therapy, is strongly limited by the low success rate of culturing cancerous cells under organoid conditions.

  Study EGAS50000000605

• Single cell transcriptomics in expanded Tregs of APS-1 patients

  Aim of study: Exploring single cell transcriptomics in expanded Tregs of with autoimmune polyendocrine syndrome type 1 (APS-1) patients (global gene expression, immune panel, TCR). N= 9 APS-1 patients and 9 healthy controls Protocol 1. Regulatory T cells (Tregs) were isolated from EDTA-blood, activated with anti-CD3/CD28 and IL2, expanded in cell culture for 14 days and stored at -150℃ in AB serum or FBS + 5% DMSO until use. 2. Cells were thawed, run through a Live/Dead column and carefully counted. 10000 Tregs were used as input for the 10x protocol. 3. cDNA was generated and amplified using Single Cell VDJ 5’ Gel Beads, Chromium Next GEM Chip K Single cell Kit, Next GEM Single Cell 5’ GEM Kit v2 and 5’ Feature Barcode Kit (all from 10x genomics). 4. Gene expression (GEX) libraries were generated by using the Library Construction Kit from 10x genomics. 5. T-cell receptor (TCR) libraries were generated using the Single Cell Human TCR Amplification Kit and the Library Construction Kit from 10x genomics. 6. 10x Genomics Target Hybridization Kit and Human Immunology Panel was used with GEX libraries. Only 8 patient samples and 8 control samples were used for the immune panel. 7. All libraries have a unique sample index (Dual Index TT Set A). 8. All samples were sequenced on Illumina NovaSeq 6000 on a NovaSeq S2 flow cell.

  Study EGAS50000000181

• Multi-omics bulk and single-cell profiling of epithelioid sarcoma

  Epithelioid sarcoma (EpS) is an aggressive sarcoma, characterized by the loss of SMARCB1 expression. EpS is traditionally classified as distal or proximal according to clinicopathological features, but its molecular characteristics remain largely unknown. To establish an EpS molecular classification and uncover determinants of inter- and intra-patient heterogeneity in EpS, we used multi-omics profiling and integrated the genomic, transcriptional and methylome landscapes with single-cell RNA sequencing on fresh samples as well as spatial transcriptomics. We identified two molecular subtypes of EpS: “distal-like” and “proximal-like”, which were distinct from the histological subtypes. Distal-like tumors include all morphological distal subtypes and a subset of proximal and hybrid EpS. They express a specific molecular single-cell derived epithelial-to-mesenchymal transition signature, which associates with improved patient survival and includes Desmoglein 2 (DSG2), which we identify as a potential routine immunohistochemical biomarker to improve diagnosis. Distal-like EpS also display increased peri-tumoral CD8+ T cell infiltrates and specific tumor-immune cell interactions. Conversely, proximal-like EpS harbor a higher inter-tumoral molecular heterogeneity with some cases resembling other SMARCB1-deficient tumors by DNA methylation profiling, and increased intra-tumoral pro-tumoral macrophage infiltrates. Our study introduces a novel molecular classification of EpS with prognostic value, associated diagnostic biomarkers and distinct tumor immune microenvironments, which extends beyond traditional clinicopathological classification. This paves the way for precision medicine-based therapeutic strategies for patients with EpS.

  Study EGAS50000000973

• Decoding Human Endogenous Retrovirus Expression in Liver Metastatic Colorectal Cancers: Implications for Diagnosis and Prognosis

  The association between Human Endogenous RetroVirus (HERV) expression and cancer can be considered to both reflect the epigenetic deregulation of cancer and facilitate oncogenesis through its promoter/enhancer action on genes. More than 50% of colorectal cancers develop liver metastases and HERV has never been studied in this setting. We collected 400 RNAseq samples from more than 200 patients with primary and liver metastases from published and unpublished databases. We corroborated our findings with tandem mass spectrometry from primary colorectal cancer and single-cell RNAseq from liver metastatic colorectal cancer (lmCRC). We identified global stability in HERV expression between liver metastases and primary colorectal cancers, suggesting an early oncogenic footprint. Therefore, we proposed a list of seventeen HERV loci for lmCRC characterization with ERV-derived peptides but no expression in normal tissues. Among them, we discovered six prognostic markers after surgical resection of lmCRC. Four of them, HERVH_Xp22.32a, HERVH_20p11.23b, HERVH_13q33.3, HERVH_13q31.3, showed poor prognostic value when expressed (log-rank p-value 0.028, 0.0083, 9e-4, 0.05, respectively). HERVH_Xp22.2c (log-rank p-value 0.032) and HERVH_8q21.3b (only in multivariate models) were associated with better prognosis when expressed. These markers showed a cumulative effect on survival when expressed. Some HERVs were associated with decreased cytotoxic immune cells and most of them correlated with cell cycle pathways. Our findings provide new insights into lmCRC classification and suggest novel prognostic biomarkers and potential therapeutic targets.

  Study EGAS50000000307

• The Dynamic Immune Behavior of Primary and Metastatic Ovarian Carcinoma

  Patients with high-grade serous ovarian carcinoma are usually diagnosed at an advanced stage and often develop resistance to conventional therapy. Epithelial ovarian tumors have been shown to be rich in CD8+ lymphocytes, which infiltrate both primary and omental sites. In recent years, several prominent studies have highlighted a positive association between the presence of these immune cells and a good prognosis; however, clinical trials of immune checkpoint inhibitors have so far shown unsatisfactory results in ovarian cancer. In this study, we used an innovative combination of single-cell sequencing and spatial transcriptomics to identify the molecular mechanisms that lead to immunosuppression in primary and metastatic HGSC. Primary tumors consistently showed a more active immune microenvironment than did omental tumors. In addition, while several immune cell subtypes were present in all samples, their relative abundance differed between primary and omental tumors. Finally, we found that primary tumors were mostly populated by dysfunctional CD4 and CD8 T cells in later stages of differentiation, while omental tumors were mostly populated by naïve immune cells; this, in turn, was correlated with expression changes in the interferon α and γ pathways in epithelial cells, showing that cross-communication between the epithelial and immune compartments is important for immune suppression in ovarian cancer. These findings can impact the timing and choice of immune-modulating drugs and should be considered when planning treatment for patients with high-grade serous ovarian carcinoma.

  Study EGAS50000000038

• Cancer-associated fibroblasts promote drug resistance in ALK-driven lung adenocarcinoma cells by upregulating lipid biosynthesis

  Targeted therapy interventions using tyrosine kinase inhibitors (TKIs) provide encouraging treatment responses in ALK-rearranged lung adenocarcinomas, yet resistances occur almost inevitably. Apart from tumor cell-intrinsic resistance mechanisms, accumulating evidence supports a role of cancer-associated fibroblasts (CAFs) in affecting the therapeutic vulnerability of lung cancer cells. Here, we aimed to investigate underlying molecular networks shaping the therapeutic susceptibility of ALK-driven lung adenocarcinoma cells via tumor microenvironmental cues using three-dimensional (3D) spheroid co-culture settings. We show that CAFs promote therapy resistance of lung tumor cells against ALK inhibition by reducing apoptotic cell death and increasing cell proliferation. Using single-cell RNA-sequencing analysis, we show that genes involved in lipogenesis constitute the major transcriptional difference between TKI-treated homo- and heterotypic lung tumor spheroids. CAF-conditioned medium and CAF-secreted factors HGF and NRG1 were both able to promote resistance of 3D-cultured ALK-rearranged lung tumor cells via AKT signaling, which was accompanied by enhanced de novo lipogenesis and supression of lipid peroxidation. Notably, simultaneous targeting of ALK and SREBP-1 was able to overcome the established CAF-driven lipid metabolic-supportive niche of TKI-resistant lung tumor spheroids. Our findings highlight a crucial role of CAFs in mediating ALK-TKI resistance via lipid metabolic reprogramming and suggest new ways to overcome resistance towards molecular directed drugs by targeting vulnerabilities downstream of oncogenic signaling.

  Study EGAS50000000135

• Diagnostic Genomic Analysis is Prognostic in AYA ALL Patients Treated on a MRD-Stratified Paediatric Protocol

  Background: Over the last decade, next generation sequencing has enabled classification of multiple new recurrent genomic drivers of acute lymphoblastic leukemia (ALL). We provide a comprehensive genomic survey of adolescent and young adult (AYA) patients treated on the Australasian Leukaemia and Lymphoma Group (ALLG) ALL06 study. ALL06 clinical trial described the efficacy and tolerability of AYA ALL patients aged 15-39 years treated on a single minimal residual disease (MRD)-stratified platform. In this study we have taken the ALL06 cohort and retrospectively performed comprehensive genomic analyses assessing response of AYA patients with various genomic subtypes treated with a paediatric approach. Methods: Genomic risk classification of B- and T-ALL patients (n=57) enrolled to the study was based on mRNA-Sequencing, MLPA, immunophenotyping and cytogenetics. Results: Using this approach, 36/40 (90%) B- and 13/17 (76.5%) T-ALL patients were classified according to genomic risk. We have linked outcome data with each patient’s genomic risk and found a strong correlation between adverse genomic risk and MRD at the end of consolidation, translating to inferior overall and relapse free survival. Patients with adverse risk genomics who achieved negative MRD status had improved responses compared to those with persistent MRD. Patients with standard risk genomics had excellent responses regardless of their MRD status. Conclusion: Our data argue strongly for incorporation of a genomic risk classification into future ALL treatment paradigms at the time of diagnosis.

  Study EGAS50000000752

• Ultra-sensitive tumor-informed ctDNA monitoring of treatment response in advanced esophagogastric cancer patients

  To explore whether ultra-sensitive circulating tumor DNA (ctDNA) profiling enables earlier prediction of treatment response and detection of disease progression, we applied NeXT Personal, an ultra-sensitive bespoke tumor-informed liquid biopsy platform, to profile tumor samples from the KeyLargo study, a phase II trial in which metastatic esophagogastric cancer (mEGC) patients received capecitabine, oxaliplatin, and pembrolizumab. A total of 24 patients were evaluated, and all were ctDNA-positive at baseline. ctDNA levels varied from 406,067 down to 1.5 parts per million (PPM) with a median limit of detection of 2.03 PPM. ctDNA dynamics were highly correlated with changes in tumor size (ρ = 0.59, p = 7.3×10-9). Lack of early molecular response (50% or greater decrease in ctDNA levels at first available time point after 30 days, C2D1 or C4D1) was associated with worse overall survival (OS) (HR 6.6, 95% CI 1.8-24.1, p = 0.005) and progression-free survival (PFS) (HR 15.4, 95% CI 2.7-87.0, p = 0.002). Lack of molecular clearance of ctDNA was associated with worse OS (HR 6.9, 95% CI 1.5-30.8, p = 0.012) and PFS (HR 19.2, 95% CI 2.4-152.8, p = 0.005). Molecular progression (ctDNA increase) preceded imaging-derived progression by a median lead time of 65 days. These results suggest that ultra-sensitive liquid biopsy approaches could improve treatment decision-making for mEGC patients receiving chemotherapy and immunotherapy.

  Study EGAS50000001224

• Delineating pediatric brain tumor progression using single-nuclei sequencing

  Pediatric brain cancers are the most common solid tumors in children with varying survival rates depending on molecular features of cancer cells and associated microenvironment. Several studies have used single-cell RNA sequencing to define transcriptional programs active in cancer cells and highlighted supportive roles of non-cancerous cells embedded in these tumors. Although application of single-cell RNA sequencing technology has expanded our understanding of intratumoral heterogeneity and underlying transcriptional mechanisms in brain tumors at diagnosis, mechanisms of disease remains elusive. Therefore, the purpose of this proposal is to characterize cellular and transcriptional states of malignant and non-malignant compartments at primary diagnosis and at relapse. In this study, we performed single-nuclei sequencing of 42 patient samples taken at various disease stages (initial diagnosis, progressive disease, and relapsed disease) from 20 patients with low-grade glioma or Atypical Teratoid Rhabdoid Tumor (ATRT). Bulk whole genome and RNA-seq data for these samples are available from the Children’s Brain Tumor Network and the Gabriella Miller Kids First Data Resource Center. Findings from our study will help define the mechanistic shifts of tumor and non-cancerous cells over course of disease progression. Collectively with the single-cell Pediatric Cancer Atlas funded, this project lays the foundation of an important data set focused on lethal relapsed cancers that can continue to be expanded to identify cancer-intrinsic and extrinsic microenvironmental factors supporting recurrence of pediatric brain tumors.

  Study EGAS50000001288

• Gut metagenome associations with extensive digital health data in a volunteer-based EstMB cohort second timepoint data

  As part of the Estonian Biobank, we established the Estonian Microbiome Cohort which includes stool, oral and plasma samples from 2,509 participants and is supplemented with multi-omic measurements, questionnaires, and regular linkages to national electronic health records. Out of the 2,509 Estonian Microbiome Project participants, a sub-cohort of over 300 individuals provided an additional stool sample after a median follow-up period of 4.4 years. All participants of the EstMB cohort gave informed consent for the data and samples to be used for scientific purposes. In the second time point, the participants were instructed to send the samples by post immediately after taking the sample and time their sample collection to avoid shipping on the weekends. In the study, by utilising retrospective medication usage data from the electronic health records and the Estonian microbiome cohort shotgun metagenomics data set (n = 2,509), we systematically evaluate the long-term effects of antibiotics and human-targeted medications on the gut microbiome. We show that past usage of medications is associated with the gut microbiome. For example, the effects of antibiotics, psycholeptics, antidepressants, proton pump inhibitors, and beta-blockers are detectable several years after use. Furthermore, by analysing a subcohort (n >300 individuals) from the second time point, we show that similar changes in the gut microbiome occur after treatment initiation or discontinuation, possibly indicating causal effects.

  Study EGAS50000001181

• Genomic landscapes of endometrioid and mucinous ovarian cancers and morphologically similar tumor types

  While endometrioid and mucinous ovarian cancers represent nearly a quarter of ovarian cancers, their molecular characteristics and pathologic origins are poorly understood. To characterize the genomic and epigenomic alterations of these ovarian cancer subtypes and evaluate links to morphologically similar tumors from other sites, we performed sequence, copy number, mutation signature, rearrangement, and methylation analyses from tumor samples and matched normal tissues of 133 patients. These analyses included patients with ovarian endometrioid (n=44), ovarian mucinous (n=43), uterine endometrial (n=15), and gastrointestinal mucinous cancers (n=31). In addition to identifying genes previously known to be involved in these tumors, we identified alterations in RAD51C, NOTCH4, SMARCA1/4 and JAK1 in ovarian endometrioid, ESR1 in uterine endometrioid, and SMARCA4 in ovarian mucinous cancers. Whole genome sequencing revealed rearrangements involving PTEN, NF1, and NF2 in ovarian endometrioid tumors, and NF1 and MED1 in ovarian mucinous cancers. The number of alterations and affected genes in ovarian and uterine endometrioid cancers were similar as were genome-wide methylation profiles, suggesting that these tumors may have a similar tissue of origin. By contrast, methylation patterns in ovarian mucinous cancers were more similar to other gastrointestinal mucinous cancers, including mucinous cancers of the stomach, colon, and pancreas. These analyses provide insights into the genomic landscapes and origins of mucinous and endometrioid ovarian carcinomas, providing new avenues for early clinical intervention and management of patients with these cancers.

  Study EGAS50000000523

• Systematic Growth Factor Profiling Platform for 3D Tumor Models Reveals EstradiolResponsive Cellular Mechanisms of Immunotherapy Resistance

  Epithelial ovarian cancer (EOC) exhibits significant molecular heterogeneity, which is thought to be influenced by the tumor microenvironment (TME). However, the specific factors driving this heterogeneity remain poorly understood. In this study, we integrated growth factoromics and genomic data to identify TME factors contributing to EOC heterogeneity and immunosuppression. Our analysis revealed distinct subgroups of EOC samples with differential responsiveness to specific growth factors, particularly estradiol and R-spondin. Notably, samples with high responsiveness to estradiol and Wnt (C1 cluster) were enriched in the Mesenchymal (MES) molecular subtype, while samples with high responsiveness to R-spondin (C2 cluster) were associated with the Differentiated (DIF) and Immunoreactive (IMR) subtypes. We discovered that the G protein-coupled estrogen receptor (GPER), rather than the classical estrogen receptor alpha (ERα), was the key mediator of estrogen responsiveness and a determinant of molecular subtype classification. Single-cell sequencing analysis identified subsets of estrogen-responsive malignant cell (MAL.PDCD5) and cancer-associated fibroblast (FB.TNFSF10), which were associated with immunosuppressive patterns in the TCGA cohort. Importantly, the interaction between these cell populations significantly increased upon estrogen treatment, suggesting a potential mechanism by which estrogen signaling may contribute to the development of an immunosuppressive TME. Our findings highlight the critical role of estrogen signaling in shaping the molecular heterogeneity and immunosuppressive TME of EOC and provide novel insights into potential therapeutic targets for overcoming immunosuppression in this malignancy.

  Study EGAS50000000422

• Tumor-reactive heterotypic CD8 T cell clusters from clinical samples

  Emerging evidence suggests a correlation between CD8+ T cell-tumor cell proximity and anti-tumor immune response. However, it remains unknown whether these cells exist as functional clusters that can be isolated from clinical samples. By conventional and imaging flow cytometry, we show here that from 21/21 human melanoma metastases, we could isolate heterotypic clusters, comprising CD8+ T cells interacting with one or more tumor cells and/or antigen-presenting cells (APCs). Single-cell RNA-sequencing revealed that T cells from clusters were enriched for gene signatures associated with tumor reactivity and exhaustion. Clustered T cells exhibited increased TCR clonality indicative of expansion, while TCR-matched T cells showed more exhaustion and co-modulation when conjugated to APCs than to tumor cells. T cells that were expanded from clusters ex vivo exerted on average 9-fold increased killing activity towards autologous melanomas, which was accompanied by enhanced cytokine production. Upon adoptive cell transfer into mice, T cells from clusters showed improved patient-derived melanoma control, which was associated with increased T cell infiltration and activation. Together, these results demonstrate that tumor-reactive CD8+ T cells are enriched in functional clusters with tumor cells and/or APCs and that they can be isolated and expanded from clinical samples. Typically excluded by single cell gating in flow cytometry, these distinct heterotypic T cell clusters serve as a valuable source to decipher functional tumor-immune cell interactions, while they may also be therapeutically explored.

  Study EGAS50000000785

• The circulating cell-free DNA landscape in sepsis is dominated by impaired liver clearance

  Circulating cell-free DNA (cfDNA) is a promising molecular biomarker. However, its utility in severe infection remains poorly understood. Here, we isolated cfDNA from sepsis patients and controls, demonstrating a 41-fold increase in the amount of cfDNA in circulation during disease. We used sequencing to reconstruct cfDNA methylomes, fragmentation profiles, and nucleosome footprints across 58 samples. We observed no difference in cfDNA composition between patients and controls, challenging the idea that cfDNA increases due to higher immune cell death during sepsis. Instead, we suggest that liver dysfunction prevents efficient clearance of cfDNA during disease. This was supported by fragmentation and end-motif patterns, both of which showed evidence of cfDNA being exposed to circulating nucleases for a prolonged period, proportionally to the extent of liver dysfunction. Variation in cfDNA of megakaryocyte-erythroid progenitor origin was also a significant contributor in sepsis, increasing over time. Moreover, we showed that cfDNA retains nucleosome footprints with cell type-specific gene activity information. We developed a novel approach to study nucleosome phasing that successfully recovers tissue-specific signatures. By combining this with single-cell data, we demonstrated that sepsis patients with liver dysfunction have higher amounts of cfDNA derived from Kupffer cells and the liver parenchyma. In conclusion, we present the first high-throughput multi-modal study of cfDNA during sepsis, which will serve as a reference point for future studies on the role of this biomarker in critical illness.

  Study EGAS50000001033

• Relapse CHL study

  Classical Hodgkin Lymphoma (CHL) is characterized by a complex tumor microenvironment (TME) that supports disease progression. While immune cell recruitment by Hodgkin and Reed–Sternberg (HRS) cells is well-documented, the role of non-malignant B cells in relapse remains unclear. Using single-cell RNA sequencing (scRNA-seq) on paired diagnostic and relapsed CHL samples, we identified distinct shifts in B cell populations, particularly an enrichment of naïve B cells and a reduction of memory B cells in early-relapse compared to late-relapse and newly diagnosed CHL. The enrichment of naïve B cells in early relapse biopsies was confirmed in independent validation cohorts using scRNA-seq and immunohistochemistry. Notably, naïve B cells in early-relapse samples exhibited high expression of LGALS9, an immunosuppressive molecule known to interact with T cell immunoglobulin and mucin domain 3 (TIM-3) on regulatory T cells (Tregs). Cell-cell interaction analysis revealed the importance of interactions between LGALS9+ naïve B cells and TIM-3+ Tregs in the early-relapse setting. Spatial analysis by imaging mass cytometry confirmed close proximity of LGALS9+ naïve B cells with TIM-3⁺ CD4⁺ T cells and HRS cells, pointing to their role in shaping an immunosuppressive niche. Our findings highlight a previously unrecognized population of LGALS9⁺ naïve B cells with immunoregulatory potential in early-relapse CHL and provide new insights into the spatial and transcriptional architecture of the relapsed TME in CHL.

  Study EGAS00001008222

• Noninvasive prenatal molecular karyotyping from maternal plasma

  Fetal DNA is present in the plasma of pregnant women. Massively parallel sequencing of maternal plasma DNA has been used to detect fetal trisomies 21, 18, 13 and selected sex chromosomal aneuploidies noninvasively. Case reports describing the detection of fetal microdeletions from maternal plasma using massively parallel sequencing have been reported. However, these previous reports were either polymorphism-dependent or used statistical analyses which were confined to one or a small number of selected parts of the genome. In this report, we reported a procedure for performing noninvasive prenatal karyotyping at 3 Mb resolution across the whole genome through the massively parallel sequencing of maternal plasma DNA. This method has been used to analyze the plasma obtained from 6 cases. In 5 cases, fetal microduplications or microdeletions have been detected successfully from maternal plasma. The two cases with fetal microduplications represented the first noninvasive prenatal detection of such changes from maternal plasma. In the remaining case, the plasma DNA sequencing result was consistent with the pregnant mother being a carrier of a microduplication. Simulation analyses were performed for determining the number of plasma DNA molecules that would need to be sequenced and aligned for enhancing the diagnostic resolution of noninvasive prenatal karyotyping to 2 Mb and 1 Mb. In conclusion, noninvasive prenatal molecular karyotyping from maternal plasma by massively parallel sequencing is feasible and would enhance the diagnostic spectrum of noninvasive prenatal testing.

  Study EGAS00001000439

• Targeting the p53 pathway to treat Malignant Rhabdoid and High-Risk Atypical Teratoid Rhabdoid Tumors

  Rhabdoid tumors (RTs) predominantly affect young children and are among the deadliest pediatric solid tumors. Despite multimodal therapy consisting of surgery, radiation, and chemotherapy, children with these tumors have median survival of less than one year. RTs can arise throughout the body, including the central nervous system (CNS) where they are called atypical teratoid rhabdoid tumors (AT/RTs), and in extra-CNS locations such as the kidneys and other soft tissues where they are designated malignant RTs (MRTs). We previously identified MDM2 and MDM4 as therapeutic vulnerabilities in RTs and showed that treatment with the MDM2 inhibitor idasanutlin increased survival in mice bearing MRT xenografts. However, the therapeutic potential of idasanutlin in CNS RT tumors is unknown. Single agent therapies are prone to resistance and show limited clinical benefit on their own. Therefore, we sought to identify combination strategies incorporating idasanutlin that would be effective in both CNS and extra-CNS RTs. Here, we show that the XPO1 inhibitor selinexor increased nuclear retention of p53 and potentiated idasanutlin-induced p53 pathway activation and cytotoxicity in AT/RT and MRT cell lines in vitro. Importantly, combination therapy limited acquired resistance through TP53 mutation. In vivo, combination therapy was well-tolerated, reduced tumor burden, and increased survival in orthotopic models of both AT/RT and MRT. Our results demonstrate that combining idasanutlin with selinexor is a promising therapeutic strategy for children with rhabdoid tumors.

  Study EGAS00001007680

• Whole exome sequencing data for paired non-serous endometrial and ovarian carcinomas from 27 patients with concurrent tumours.

  Introduction: Concurrent non-serous endometrial and ovarian tumours are often managed clinically as two separate primary tumours, but almost all exhibit evidence of a genomic relationship. Methodology: This study investigates the extent of relatedness using whole exome sequencing, which was performed on paired non-serous endometrial and ovarian carcinomas from 27 patients with concurrent tumours in a cohort with detailed clinicopathological annotation. Four whole exome sequencing-derived parameters (mutation, mutational burden, mutational signatures and mutant allele tumour heterogeneity scores) were used to develop a novel unsupervised model for the assessment of genomic similarity to infer genomic relatedness of paired tumours. Results: This novel model demonstrated genomic relatedness across all four parameters in all tumour pairs. Mutations in PTEN, ARID1A, CTNNB1, KMT2D and PIK3CA occurred most frequently and 24 of 27 (89 %) tumour pairs shared identical mutations in at least one of these genes, with all pairs sharing mutations in a number of other genes. Ovarian endometriosis, CTNNB1 exon 3 mutation, and progression and death from disease were present across the similarity ranking. Mismatch repair deficiency was associated with less genomically similar pairs. Discussion: Although there was diversity across the cohort, the presence of genomic similarity in all paired tumours supports the hypothesis that concurrent non-serous endometrial and ovarian carcinomas are genomically related and may have arisen from a common origin. Keywords: Co-existing; Concurrent; Endometrial carcinoma; Genomics; Mutation; Ovarian carcinoma; Paired tumours; Tumour mutational burden.

  Study EGAS00001008259

• Pseudodiastrophic dysplasia expands the known phenotypic spectrum of defects in proteoglycan biosynthesis

  Pseudodiastrophic dysplasia (PDD) is a severe skeletal dysplasia associated with prenatal manifestation and early lethality. Clinically, PDD is classified as a ‘dysplasia with multiple joint dislocations’ however, the molecular aetiology of the disorder is currently unknown. In order to identify the genetic defects underlying PDD, we performed whole exome sequencing (WES) on 3 patients from 2 unrelated families, clinically diagnosed with PDD. WES resulted in the identification of bi-allelic variants in the established skeletal dysplasia genes, B3GAT3 (Family 1) and CANT1 (Family 2). Mutations in these genes have previously been reported to cause ‘multiple joint dislocations, short stature, and craniofacial dysmorphism with or without congenital heart defects’ (B3GAT3) and Desbuquois dysplasia 1 (CANT1); disorders in the same nosological group as PDD. Follow-up of the B3GAT3 variants demonstrated significantly reduced B3GAT3/GlcAT-I expression. Downstream in vitro functional analysis revealed abolished biosynthesis of glycosaminoglycan side chains on proteoglycans. Functional evaluation of the CANT1 variant showed impaired nucleotidase activity, which results in inhibition of glycosaminoglycan synthesis through accumulation of uridine diphosphate. For the families described in this study, the PDD phenotype was caused by mutations in the known skeletal dysplasia genes B3GAT3 and CANT1, demonstrating the advantage of genomic analyses in delineating the molecular diagnosis of skeletal dysplasias. This finding expands the phenotypic spectrum of B3GAT3- and CANT1- related skeletal dysplasias to include PDD, and highlights the significant phenotypic overlap of conditions within the proteoglycan biosynthesis pathway.

  Study EGAS00001004145

• Copy number profiling of putative cancer stem from pleural effusion aspirates from breast cancer patients

  A common symptom during late stage breast cancer disease is pleural effusion, which is related to poor prognosis. In pleural effusions malignant cells can be detected, which indicates a metastatic spread from the primary tumor site. Beside the usefulness for investigation of the metastatic process, it has been a source to isolate cells with putative cancer stem cell (CSC) properties. For this study, pleural effusion aspirate from 25 metastatic breast cancer patients were processed to establish in vitro culture. Patient-derived aspirates were cultured under sphere forming conditions and isolated primary cultures were sorted for cancer stem cell populations ALDH1+ and CD44+CD24-/low and sphere forming efficiency of CSC and non-CSC subpopulation was determined. Unsorted cells had a higher potential to form spheres compared to CSC subpopulations. Subsequently, DNA was isolated from pleural effusions before and after cell sorting, to compare their copy number profiles generated by low-coverage whole genome sequencing to DNA from FFPE primary tumors or bone metastasis. Unbalanced profiles could be detected in four pleural samples before sorting, in two ALDH1+ sorted samples, in one CD44+CD24-/low and in all control samples. Additionally, spheroid culture from one pleural sample could be passaged on and was characterized. Therefore, pleural effusions are an appropriate media for enrichment of putative CSC. However, low-coverage sequencing is difficult due to low quantities of abnormal cells compared to normal cells.

  Study EGAS00001002343

• Genome-wide analysis of genetic risk factors for rheumatic heart disease in Aboriginal Australians provides support for pathogenic molecular mimicry

  Background. Rheumatic heart disease (RHD) following Group A Streptococcus (GAS) infections is heritable and prevalent in Indigenous populations. Molecular mimicry between human and GAS proteins triggers pro-inflammatory cardiac valve-reactive T-cells. Methods. Genome-wide genetic analysis was undertaken in 1263 Aboriginal Australians (398 RHD cases; 865 controls). Single nucleotide polymorphisms (SNPs) were genotyped using Illumina HumanCoreExome BeadChips. Direct typing and imputation was used to fine-map the human leukocyte antigen (HLA) region. Epitope binding affinities were mapped for human cross-reactive GAS proteins, including M5 and M6. Results. The strongest genetic association was intronic to HLA-DQA1 (rs9272622; P=1.86x10-7). Conditional analyses showed rs9272622 and/or DQA1*AA16 account for the HLA signal. HLA-DQA1*0101_DQB1*0503 (OR 1.44, 95%CI 1.09-1.90, P=9.56x10-3) and HLA-DQA1*0103_DQB1*0601 (OR 1.27, 95%CI 1.07-1.52, P=7.15x10-3) were risk haplotypes; HLA_DQA1*0301-DQB1*0402 (OR 0.30, 95%CI 0.14-0.65, P=2.36x10-3) was protective. Human myosin cross-reactive N-terminal and B repeat epitopes of GAS M5/M6 bind with higher affinity to DQA1/DQB1 alpha/beta dimers for the two risk haplotypes than the protective haplotype. Conclusions. Variation at HLA_DQA1-DQB1 is the major genetic risk factor for RHD in Aboriginal Australians studied here. Cross-reactive epitopes bind with higher affinity to alpha/beta dimers formed by risk haplotypes, supporting molecular mimicry as the key mechanism of RHD pathogenesis.

  Study EGAS00001002678

• Open_Targets_020_Epigenomes_of_Cell_Lines

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

  Study EGAS00001003138

• Evolution of the African pygmy phenotype

  The evolutionary history of the human ‘pygmy’ phenotype (small body size), a characteristic of African and Southeast Asian rainforest hunter-gatherers, is largely unknown. Here we use a genome-wide admixture mapping analysis to identify 16 genomic regions that are significantly associated with the pygmy phenotype in the Batwa, a rainforest hunter-gatherer population from Uganda (East Central Africa). The identified genomic regions have multiple attributes that provide supporting evidence of genuine association with the pygmy phenotype, including enrichments for SNPs previously associated with stature variation in Europeans and for genes with growth hormone receptor and regulation functions. To test adaptive evolutionary hypotheses we computed the haplotype-based iHS statistic and the level of population differentiation (FST) between the Batwa and their agricultural neighbors, the Bakiga, for each SNP. Both |iHS| and FST values were significantly higher for SNPs within the Batwa pygmy phenotype-associated regions than the remainder of the genome, a signature of polygenic adaptation. In contrast, when we expanded our analysis to include Baka rainforest hunter-gatherers from Cameroon and Gabon (West Central Africa) and Nzebi and Nzime neighboring agriculturalists, we did not observe elevated |iHS| or FST values in these genomic regions. Together, these results suggest adaptive and at least partially convergent origins of the pygmy phenotype even within Africa, supporting the hypothesis that small body size confers a selective advantage for tropical rainforest hunter-gatherers but raising questions about the antiquity of this behavior.

  Study EGAS00001000908

• Illumina HumanCoreExome genotyping data from the TEENAGE (TEENs of Attica: Genes and Environment) study.

  The TEENAGE study is a cross-sectional study comprising healthy Greek adolescents. The study was conducted according to the guidelines specified in the Declaration of Helsinki and all procedures involving human subjects were approved by the Institutional Review Board of Harokopio University and the Greek Ministry of Education, Lifelong Learning and Religious Affairs. The TEENAGE study target population comprised adolescent students aged 13 to 15 years attending the first three classes of public secondary schools located in the wider Athens area of Attica. The regions across Attica from which students were recruited were physically similar and all urban in nature. Sampling of regions and schools selected to participate was random. Students were enrolled from the wider urban region to ensure that they were coming from families with varying educational background and socio-economic status. Moreover, in order to exclude any vast economic disparities, all participating schools were exclusively public schools, which are free of charge. A total of twenty-four schools participated in the study out of 285 junior secondary public schools in Attica. Anthropometric, medical and dietary information was obtained. DNA samples were genotyped using Illumina HumanCoreExome BeadChips (Illumina, San Diego, CA, USA) at the Wellcome Trust Sanger Institute, Hinxton, UK. Genotypes were called using GenCall (Illumina Genome Studio) followed by zCall and quality control (QC) was performed in two stages (pre- and post-zCall).

  Study EGAS00001001001

• Whole genome sequencing with linked reads of pediatric glioblastoma samples

  Pediatric glioblastoma (pGBM) is a lethal cancer with no effective therapies. To understand mechanisms of tumor evolution in this cancer, we performed whole genome sequencing with linked reads on longitudinally resected pGBM samples. Our analyses showed that all diagnostic and recurrent samples were collections of genetically diverse subclones. Clonal composition rapidly evolved at recurrence, with less than 8% of non-synonymous single nucleotide variants being shared in diagnostic-recurrent pairs. In order to track the origins of the mutational events we observed in pGBM, we generated whole genome datasets for two patients and their parents. These trios showed that genetic variants could be (i) somatic, (ii) inherited from a healthy parent, or (iii) arose de novo in the germlines of pGBM patients. Analysis of variant allele frequencies supported a model of tumor growth involving slow-cycling cancer stem cells that give rise to fast-proliferating progenitor-like cells and to non-dividing cells. Interestingly, radiation and anti-mitotic chemotherapeutics did not increase overall tumor burden upon recurrence. These findings support an important role for slow-cycling stem cell populations in contributing to recurrences, since slow-cycling cell populations are expected to be less prone to genotoxic stress induced by these treatments and to accumulate few mutations. Our results highlight the need for new targeted treatments that account for the complex functional hierarchies and genomic heterogeneity of pGBM.

  Study EGAS00001003432

• Bioinformatic Methods and Bridging of Assay Results for Reliable Tumor Mutational Burden Assessment in Non-Small Cell Lung Cancer

  Tumor mutational burden (TMB) has emerged as a predictive biomarker of response to immune checkpoint inhibitors. Standardization of TMB measurement is essential for implementing diagnostic tools to guide treatment. Here we evaluate bioinformatic TMB analysis by whole exome sequencing (WES) in formalin-fixed, paraffin-embedded samples. In CheckMate 026, TMB was retrospectively assessed in 312 patients with non-small cell lung cancer (58% of the intent-to-treat population) who received first-line nivolumab treatment or chemotherapy. We examined the sensitivity of TMB assessment to bioinformatic filtering methods and assessed concordance between TMB data derived by WES and the FoundationOne CDx™ assay. TMB scores comprising synonymous, indel, frameshift, and nonsense mutations (all mutations) were 3.1-fold higher than data including missense mutations only, but values were highly correlated (Spearman’s r = 0.99). Scores from CheckMate 026 samples including missense mutations only were similar to those generated from data in The Cancer Genome Atlas, but those including all mutations were generally higher. Using databases for germline subtraction (instead of matched controls) showed a trend for race-dependent increases in TMB scores. Parameter variation can therefore impact TMB calculations, highlighting the need for standardization. Encouragingly, WES and FoundationOne CDx outputs were highly correlated (Spearman’s r = 0.90) and differences could be accounted for by empirical calibration, suggesting that reliable TMB assessment across assays, platforms and centers is achievable.

  Study EGAS00001003661

• Neoadjuvant combination PD-L1 plus CTLA-4 blockade in patients with cisplatin-ineligible operable urothelial carcinoma

  Immune checkpoint therapy (ICT) is being tested in the neoadjuvant setting for patients with localized urothelial carcinoma (UC), with one study reporting data in cisplatin-ineligible patients who received anti PD-L1 monotherapy. The study reported that patients with bulky tumors, a known high-risk feature defined as greater than clinical T2 disease, had fewer responses, with pathologic complete response (pCR) rate of 17%. Here, we report on the first pilot combination neoadjuvant trial (NCT02812420) with anti-PD-L1 (durvalumab) plus anti-CTLA-4 (tremelimumab) in cisplatin-ineligible patients, with all tumors identified as having high risk features (N=28). Primary endpoint was safety and we observed 6 of 28 patients (21%) with grade ≥3 immune-related adverse events, consisting of asymptomatic laboratory abnormalities (N=4), hepatitis and colitis (N=2). We also observed pCR of 37.5% and downstaging to pT1 or less in 58% of patients who completed surgery (N=24). In summary, we provide initial safety, efficacy and biomarker data with neoadjuvant combination anti-PD-L1 plus anti-CTLA-4, which warrants further development for patients with localized UC, especially cisplatin-ineligible patients with high-risk features who do not currently have an established standard-of-care neoadjuvant treatment. Publication: Gao et al. Neoadjuvant PD-L1 plus CTLA-4 blockade in patients with cisplatin-ineligible operable high-risk urothelial carcinoma. Nature Medicine volume 26, pages1845–1851(2020)

  Study EGAS00001004074

• Targets of MEK inhibition in DIPG

  The survival of children with diffuse intrinsic pontine glioma (DIPG) remains dismal, with new treatments desperately needed. In the era of precision medicine, targeted therapies represent an exciting treatment opportunity, yet resistance can rapidly emerge, playing an important role in treatment failure. In a prospective biopsy-stratified clinical trial, we combined detailed molecular profiling (methylation BeadArray, exome, RNAseq, phospho-proteomics) linked to drug screening in newly-established patient-derived models of DIPG in vitro and in vivo. We identified a high degree of in vitro sensitivity to the MEK inhibitor trametinib in samples which harboured genetic alterations targeting the MAPK pathway, including the non-canonical BRAF_G469V mutation, and those affecting PIK3R1 and NF1. However, treatment of PDX models and of a patient with trametinib at relapse failed to elicit a significant response. We generated trametinib-resistant clones in the BRAF_G469V model through continuous drug exposure, and identified acquired mutations in MEK1/2 (MEK1_K57N, MEK1_I141S and MEK2_I115N) with sustained pathway up-regulation. These cells showed the hallmarks of mesenchymal transition, and expression signatures overlapping with inherently trametinib-insensitive primary patient-derived cells that predicted an observed sensitivity to dasatinib. Combinations of trametinib with dasatinib and the downstream ERK inhibitor ulixertinib showed highly synergistic effects in vitro. These data highlight the MAPK pathway as a therapeutic target in DIPG, and show the importance of parallel resistance modelling and rational combinatorial treatments likely to be required for meaningful clinical translation.

  Study EGAS00001004495

• Assessing mitochondrial bioenergetics in coronary artery disease: A translational left ventricular tissue study in humans (The AMBITION study).

  Severe or recurrent myocardial ischemia can lead to long-term left ventricular (LV) dysfunction and heart failure in patients with coronary artery disease (CAD). In this study, we assessed for chronic metabolomic and transcriptomic adaptations within LV myocardium of patients with evidence of preoperative inducible ischemia undergoing coronary artery bypass grafting (CABG). During surgery, paired transmural LV biopsies were acquired on the beating heart from a region of inducible ischemia (determined by stress perfusion cardiovascular magnetic resonance), and a remote LV segment. From 33 patients, 63 LV biopsies were acquired on the beating heart during CABG, comparing to analysis of LV samples from 11 donor hearts. The global myocardial ATP/ADP ratio was reduced in CAD patients as compared to donor LV tissue with increased expression of oxidative phosphorylation (OXPHOS) genes encoding the electron transport chain complexes across multiple cell types. Paired analyses of biopsies obtained during CABG from LV segments with or without inducible ischemia revealed no significant difference in the ATP/ADP ratio, broader metabolic profile or expression of ventricular cardiomyocyte genes implicated in OXPHOS. Finally, differential metabolite analysis suggested dysregulation of several intermediates in patients with impaired LVEF; this notably included the tricyclic acid cycle metabolite, succinate. Overall, our results suggest that viable human myocardium in patients with stable CAD has global alterations in bioenergetic and transcriptional profile without large regional differences between areas with or without inducible ischemia.

  Study EGAS00001007351

• Spatio-temporal Profiling of a Rhabdoid Tumor Case Study

  Background Immune checkpoint inhibition (ICI) has revolutionized oncology, offering extended survival and long-term remission in previously incurable cancers. While highly effective in tumors with high mutational burden, lowly mutated cancers, including pediatric malignancies, present low response rate and limited predictive biomarkers. Patients and methods We present a framework for the identification and validation of tumor-reactive T cells as a biomarker to quantify ICI efficacy and as candidates for a personalized TCR-T cell therapy. Therefore, we profiled a pediatric malignant rhabdoid tumor patient with complete remission after ICI therapy using deep single-cell T cell receptor (TCR) repertoire sequencing of the tumor microenvironment (TME) and the peripheral blood. Results Tracking T cell dynamics longitudinally from the tumor to cells in circulation over a time course of 12 months revealed a systemic response and durable clonal expansion of tumor-resident and ICI-induced TCR clonotypes. We functionally validated tumor reactivity of TCRs identified from the TME and the blood by co-culturing patient-derived tumor cells with TCR-engineered autologous T cells. Here, we observed unexpectedly high frequencies of tumor-reactive TCR clonotypes in the TME and confirmed T cell dynamics in the blood post-ICI to predict tumor-reactivity. Conclusion These findings strongly support spatio-temporal tracking of T cell activity in response to ICI to inform therapy efficacy and to serve as a source of tumor-reactive TCRs for personalized TCR-T designs.

  Study EGAS00001008174

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

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

  Study EGAS00001008221

• Whole exome sequencing of an invasive diffuse intrinsic pontine glioma sampled from different sites

  The lateral ventricle subventricular zone (SVZ) is a frequent site of high-grade glioma (HGG) spread and tumor invasion of the SVZ is a predictor of a worse clinical prognosis. A range of HGG types invade the stem cell niche of the SVZ, including both adult glioblastoma and pediatric high-grade gliomas such as diffuse intrinsic pontine glioma (DIPG). The cellular and molecular mechanisms mediating this frequent invasion of the SVZ are poorly understood. Here we demonstrate that neural precursor cell (NPC):glioma cell communication underpins the propensity of glioma to colonize the SVZ. SVZ NPCs secrete chemoattractant signals toward which glioma cells home. Biochemical and proteomic analyses of factors secreted by SVZ NPCs revealed a discrete list of candidate proteins. Necessity and sufficiency testing implicated the neurite outgrowth-promoting factor pleiotrophin, along with three required binding partners (secreted protein acidic and rich in cysteine (SPARC), SPARC-like protein 1, and heat shock protein 90B), as key mediators of this chemoattractant effect. Pleiotrophin protein expression is strongly enriched in the SVZ, and knockdown of pleiotrophin expression starkly reduced glioma invasion of the SVZ in the adult murine brain. Pleiotrophin, in complex with the three binding partners, activated the Rho/ROCK pathway in DIPG cells, and inhibition of ROCK resulted in decreased DIPG invasion toward SVZ neural precursor cell-secreted factors. These findings demonstrate a pathogenic role for neural precursor cell – glioma interactions and potential therapeutic targets to limit glioma invasion.

  Study EGAS00001002326

• Targeted sequencing of candidate regions on chromosome 22q predisposing to multiple schwannomas

  Constitutional LZTR1 or SMARCB1 pathogenic variants have been found in ~86% of familial and ~40% of sporadic schwannomatosis cases. Hence, we performed massively parallel sequencing of the entire LZTR1, SMARCB1 and NF2 genomic loci in 35 individuals with schwannomas negative for constitutional first-hit pathogenic variants in the LZTR1/SMARCB1/NF2 coding sequences, however with 22q deletion and/or a different NF2 pathogenic variant in each tumor, including six cases with only one tumor available. Furthermore, we verified whether any other LZTR1/SMARCB1/NF2 (likely) pathogenic variant(s) could be found in 16 cases carrying a SMARCB1 constitutional variant in the 3’ untranslated region (3’-UTR) c.*17C>T, c.*70C>T or c.*82C>T. As no additional variants were found, functional studies were performed to clarify the effect of these 3’-UTR variants on the transcript. The 3’-UTR variants c.*17C>T and c.*82C>T showed pathogenicity by negatively affecting SMARCB1 transcript level. Two novel deep intronic SMARCB1 variants, c.500+883T>G and c.500+887G>A, resulting in out-of-frame missplicing of intron 4, were identified in two unrelated individuals. Further resequencing of the entire repeat - masked genomics sequences of chromosome 22q in individuals negative for pathogenic variants in the SMARCB1/LZTR1/NF2 coding- and non-coding regions revealed five potential schwannomatosis-predisposing candidate genes, i.e. MYO18B, NEFH, SGSM1, SGSM3 and SBF1, pending further verification.

  Study EGAS00001005680

• Genetic variability in exon 1 of the glucocorticoid receptor gene NR3C1 is associated with postoperative complications

  Patients having major surgery experience inflammation post-operatively, which may contribute to post-operative morbidity. Endogenous glucocorticoids (GCs) are an essential part of the stress response, but this response varies between individuals, which may in turn affect clinical outcome in general, and post-operative inflammation, specifically. Exon 1 of the NR3C1 gene, encoding the glucocorticoid receptor, contains an established region of differential regulation. DNA methylation patterns in this region have been found to differ between individuals. We investigated the methylation status and genotype in the Cytosine-phosphate-Guanine (CpG) island in exon 1 of NR3C1 in 24 patients who underwent major abdominal surgery and correlated it with post-operative complications. DNA was extracted from peripheral blood leukocytes and underwent targeted bisulfite sequencing of the CpG island. Complications were graded according to the Clavien-Dindo classification and 14 out of 24 had post-operative complications. Multifactorial and partial least square methods were used to analyse the data. A homogenous demethylated pattern was observed in all patients and no single CpG methylation was associated with post-operative complications. Four single nucleotide polymorphisms (SNPs) were significantly correlated with higher Clavien-Dindo scores. Genetic variability in the chr5:143,402,505-143,405,805 region of exon 1 of the glucocorticoid receptor gene NR3C1, but not DNA methylation, was associated with more severe post-operative complications in patients having major abdominal surgery. This indicates that the response to GCs may be of clinical importance for inflammatory conditions.

  Study EGAS00001005737

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

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

  Study EGAS00001005738

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

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

  Study EGAS00001005747

• Platelet RNAseq data for SLFN14 K219N patients

  Pathogenic missense variants in SLFN14, which encodes an RNA endoribonuclease protein that regulates rRNA degradation, are known to cause inherited thrombocytopenia with impaired platelet aggregation and ATP secretion. Despite rather mild laboratory defects, these patients display an obvious bleeding phenotype. The function of SLFN14 in megakaryocyte (MK) and platelet biology is unknown. This study aims to characterize the platelet transcriptome in patients with a SLFN14 K219N variant and model the disease in the immortalized megakaryocyte cell line imMKCL. Total platelet RNA was sequenced for two patients and 19 healthy controls. Differential gene expression analysis yielded a total of 2999 and 2888 significantly (|log2FC|>1, FDR<0.05) up- and downregulated genes, respectively. Remarkably, these downregulated genes were not enriched for any biological pathway while upregulated genes were enriched for pathways involved in (mitochondrial) translation and transcription with a significant upregulation of 134 ribosomal protein genes (RPG). Heterozygous and homozygous SLFN14 K219N imMKCL showed a defect in MK differentiation that aggravated during cell passaging of undifferentiated cells and in proplatelet formation. SLFN14 defective platelets and MK showed signs of rRNA degradation while this was absent in undifferentiated imMKCL cells. Upregulation of RPG through increased mTOR signaling in SLFN14 K219N MK seems to be a compensatory response for rRNA degradation. Indeed, mTOR inhibition with rapamycin resulted in further enhanced rRNA degradation in SLFN14 K219N MK. Taken together, our study indicates dysregulation of mTOR-regulated ribosomal biogenesis as the disease mechanism for SLFN14-related thrombocytopenia.

  Study EGAS00001006339

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

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

  Study EGAS00001006115

• Plasma-Seq of patients with metastatic prostate cancer

  Study 1 In this study we analysed patients with metastatic prostate cancer to scan their tumor genomes noninvasively in plasma DNA. We wanted to make whole-genome analysis from plasma DNA amenable to clinical routine applications and developed an approach based on a benchtop high-throughput platform, i.e. Illuminas MiSeq instrument. We performed whole-genome sequencing from plasma at a shallow sequencing depth to establish a genome-wide copy number profile of the tumor at low costs within 2 days. The genome-wide profiling in the plasma of our patients revealed multiple copy number aberrations including those previously reported in prostate tumors, such as losses in 8p and gains in 8q. High-level copy number gains in the AR locus were observed in patients with CRPC but not with CSPC disease. We identified the TMPRSS2-ERG rearrangement associated 3-Mbp deletion on chromosome 21 and found corresponding fusion plasma fragments in these cases. In an index case multiregional sequencing of the primary tumor identified different copy number changes in each sector, suggesting multifocal disease. Our plasma analyses of this index case, performed 13 years after resection of the primary tumor, revealed novel chromosomal rearrangements, which were stable in serial plasma analyses over a 9 months period, which is consistent with the presence of one metastatic clone. Our approach provides specific genomic signatures within 2 days which may therefore serve as “liquid biopsy”.

  Study EGAS00001000451

• Transcriptomic profiles of chronic lymphocytic leukemia before and after frontline therapy: 5-year results from the randomized CLL14 study

  CLL14 (NCT02242942) is a multinational, open-label Phase III study designed to compare venetoclax-obinutuzumab (Ven-Obi) versus chlorambucil-obinutuzumab (Clb-Obi) in previously untreated patients with CLL and co-existing conditions. Between Aug 7, 2015, and Aug 4, 2016, 432 patients were enrolled and randomly assigned to receive either venetoclax plus obinutuzumab (Ven-Obi, n=216) or chlorambucil plus obinutuzumab (Clb-Obi, n=216). In both study arms, patients received fixed-duration treatment of overall 12 cycles, each cycle 28 days long, and entered subsequent post-treatment surveillance. Peripheral blood samples were collected before the start of therapy, as well as every 3-6 months during and after therapy, including at the time of disease progression. The primary endpoint of the study was progression-free survival, secondary endpoints included overall response rate, minimal residual disease response, and overall survival. The primary read-out was conducted after a median follow-up of 28.1 months, when the primary endpoint was met with 30 primary end-point events (disease progression or death) in the Ven-Obi arm and 77 in the Clb-Obi arm (hazard ratio, 0.35; 95% confidence interval [CI], 0.23 to 0.53;P<0.001). Follow-up is ongoing until 9 years after the last patient enrollment. In an exploratory analysis, bulk RNAseq was conducted with CD19-enriched blood from patients before start of therapy (418 patients) and at first relapse (44 patients).

  Study EGAS00001006596

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

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

  Study EGAS00001006879

• RB1 Gene Inactivation by Chromothripsis in Human Retinoblastoma

  Retinoblastoma is a rare childhood cancer of the developing retina. Most retinoblastomas initiate with biallelic inactivation of the RB1 gene through diverse mechanisms including point mutations, nucleotide insertions, deletions, loss of heterozygosity and promoter hypermethylation. Recently, a novel mechanism of retinoblastoma initiation was proposed. Gallie and colleagues discovered that a small proportion of retinoblastomas lack RB1 mutations and had MYCN amplification [1]. In this study, we identified recurrent chromosomal, regional and focal genomic lesions in 94 primary retinoblastomas with their matched normal DNA using SNP 6.0 chips. We also analyzed the RB1 gene mutations and compared the mechanism of RB1 inactivation to the recurrent copy number variations in the retinoblastoma genome. In addition to the previously described focal amplification of MYCN and deletions in RB1 and BCOR, we also identified recurrent focal amplification of OTX2, a transcription factor required for retinal photoreceptor development. We identified 10 retinoblastomas in our cohort that lacked RB1 point mutations or indels. We performed whole genome sequencing on those 10 tumors and their corresponding germline DNA. In one of the tumors, the RB1 gene was unaltered, the MYCN gene was amplified and RB1 protein was expressed in the nuclei of the tumor cells. In addition, several tumors had complex patterns of structural variations and we identified 3 tumors with chromothripsis at the RB1 locus. This is the first report of chromothripsis as a mechanism for RB1 gene inactivation in cancer.

  Study EGAS00001000598

• T cell receptor repertoire sequencing reveals chemotherapy-driven clonal expansion in colorectal liver metastases

  Colorectal liver metastasis (CLM) is a leading cause of colorectal cancer mortality, and the response to immune checkpoint inhibition (ICI) in microsatellite stable CRC has been disappointing. Administration of cytotoxic chemotherapy may cause increased density of tumour infiltrating T cells, which has been associated with improved response to ICI. This study aimed to quantify and characterize T cell infiltration in CLM using T cell receptor (TCR) repertoire sequencing. Eighty-five resected CLM from patients included in the Oslo CoMet study were subjected to TCR repertoire sequencing. Thirty-five and 15 patients had received neoadjuvant chemotherapy (NACT) within a short or long interval, respectively, prior to resection, while 35 patients had not been exposed to NACT. T cell fractions were calculated, repertoire clonality was analysed based on Hill evenness curves, and TCR sequential convergence was assessed using network analysis. Increased T cell fractions (10.6% vs 6.3%) were detected in CLM exposed to NACT within a short interval prior to resection, while modestly increased clonality was observed in NACT exposed tumours independently of the timing of NACT administration and surgery. While private clones made up >90% of detected clones, network connectivity analysis revealed that public clones contributed the majority of TCR sequence convergence. TCR repertoire sequencing can be used to quantify T cell infiltration and clonality in clinical samples. This study provides evidence to support chemotherapy-driven T cell clonal expansion in CLM in a clinical context.

  Study EGAS00001007136

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

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

  Study EGAS00001007204

• Evolutionary origin and methylation status of human intronic CpG islands that are not present in mouse

  Imprinting of the human RB1 gene is due to the presence of a differentially methylated CpG island (CGI) in intron 2, which is part of a retrogene derived from the PPP1R26 gene on chromosome 9. The murine Rb1 gene does not have this retrogene and is not imprinted. We have investigated whether the RB1/Rb1 locus is unique with respect to these differences. For this we have compared the CGIs from human and mouse by in silico analyses and found out that there are at least 2033 human intronic CGIs that are not present in the mouse. Among these CGIs, 104 show sequences similarities elsewhere in the human genome, which suggests that they arose by retrotransposition. We could narrow down the time points when most of these CGIs appeared during evolution. Their methylation status was analysed in two monocyte methylome datasets from whole genome bisulfite sequencing and in 18 published methylomes. Five CGIs appear to be differentially methylated, which suggests that they are subject of genomic imprinting. They are located in introns of the following genes: RB1, ASRGL1, PARP11, PDXDC1 and MYO1D. Four of these CGIs are parts of retrogenes. Our study supports the notion that imprinting builds on host defence mechanisms by which the genome protects itself against retrotransposed and foreign DNA and that the epigenetic fate of the retrotransposed DNA depends on selective forces at the integration site.

  Study EGAS00001000719

• Cistrome-partitioning reveals convergence of somatic mutations and risk-variants on master transcription regulators in primary prostate tumors

  Thousands of noncoding somatic Single Nucleotide Variants (SNVs) of unknown function are reported in tumors. Partitioning the genome according to cistromes, reveals the enrichment of somatic SNVs in prostate tumors as opposed to adjacent normal tissue cistromes of master transcription regulators, including AR, FOXA1 and HOXB13. This parallels enrichment of prostate cancer genetic predispositions over these transcription regulators’ tumor cistromes, exemplified at the 8q24 locus harboring both risk-variants and somatic SNVs in cis-regulatory elements, upregulating MYC expression and altering the binding of transcription regulators to DNA. However, Massively-Parallel Reporter Assays reveal that few SNVs can alter the transactivation potential of individual CREs. Instead, SNVs accumulate, similarly to inherited riskvariants, in cistromes of master transcription regulators required for prostate cancer development. Difficulties in inferring the biological significance of noncoding mutations have limited their inclusion in precision genomics medicine pipelines. Most attempts to delineate a role for noncoding mutations relied on detecting evidence for positive selection within individual CREs, such as reported for the TERT gene promoter. By considering the enrichment of noncoding mutations in cistromes as opposed to individual CREs, we reveal their specificity towards master transcription regulators that promote prostate cancer development, a feature shared with inherited risk-variants. Overall, our work provides a blueprint for the functional interpretation of noncoding mutations in genomic tests relying on defining cis-regulatory units according to cistrome-partitioning to identify cancer driver transcription regulators.

  Study EGAS00001003928

• Characterization of HCV-specific CD4 T cells during DAA-Therapy

  Background: Chronic HCV-infection is characterized by a severe impairment of HCV-specific CD4 T cell help that is driven by chronic antigen stimulation. We aimed to study the fate of HCV-specific CD4 T cells after viral elimination. Methods: HCV-specific CD4 T cell responses were longitudinally analyzed using MHC class II tetramer-technology, multicolor flow cytometry and RNA sequencing in a cohort of chronically HCV-infected patients undergoing therapy with direct-acting antivirals. In addition, HCV-specific neutralizing antibodies and CXCL13 levels were analyzed. Results: We observed that the frequency of HCV-specific CD4 T cells increased within two weeks after initiation of DAA therapy. Multicolor flow cytometry revealed a downregulation of exhaustion and activation markers and an upregulation of memory-associated markers. While cells with a Th1 phenotype were the predominant subset at baseline, cells with phenotypic and transcriptional characteristics of follicular T helper cells increasingly shaped the circulating HCV-specific CD4 T cell repertoire, suggesting antigen-independent survival of this subset. These changes were accompanied by a decline of HCV-specific neutralizing antibodies and the germinal center activity. Conclusion: We identified a population of HCV-specific CD4 T cells with a follicular T helper cell signature that is maintained after therapy-induced elimination of persistent infection and may constitute an important target population for vaccination efforts to prevent re-infection and immunotherapeutic approaches for persistent viral infections.

  Study EGAS00001003950

• Molecular Attributes Underlying Central Nervous System and Systemic Relapse in Diffuse Large B-cell Lymphoma

  Central nervous system (CNS) relapse in diffuse large B-cell lymphoma remains a significant conundrum. While clinical risk scores and biomarkers allow the identification of patients with higher risk of CNS relapse, their predictive ability is low. Herein, we report findings from a comprehensive study examining gene expression changes, mutations and clonal evolution. Firstly, considering diagnostic samples from patients with CNS relapse, systemic relapse or no relapse, we found a striking enrichment of gene expression changes reflecting protein biosynthesis in patients who experienced CNS or systemic relapse. Furthermore, diagnostic specimens of relapsing cases were characterized by low immune infiltrates. The most discriminative gene feature in our analysis was SLAMF1, high expression of which was associated with favourable outcomes. Next, we assessed whether specific gene mutations were associated with subsequent CNS or systemic relapse. Activated B-cell-like genes such as PIM1, CD79B and PRDM1 were more frequently mutated in cases with CNS relapse, whereas TP53, BCL2 and FOXO1 mutations were associated with systemic relapse. Lastly, we performed tumor evolution analysis in five patients, identifying significant clonal divergence that accompanied CNS relapse. In summary, we identified both common biology underlying both CNS and systemic relapse, independent of cell-of-origin (COO), and specific COO-enriched gene mutations that conferred a higher risk of either CNS or systemic relapse. Moreover, we identified a pattern of clonal evolution that can explain difficulties in both predicting and treating CNS relapse.

  Study EGAS00001004057

• UK10K_RARE_HYPERCHOL

  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.Familial Hypercholesterolemia is a condition where the affected person has a consistently high level of LDL which can lead to early clogging of the coronary arteries. All patients selected for this study will have been found not to carry the common APOB and PCSK9 mutations, and to have no detectable LDLR mutations by testing for 18 common mutations. For further information with regard to this cohort please contact Steve Humphries(steve.humphries@ucl.ac.uk).

  Study EGAS00001000129

• Integration of genomics and histology reveals diagnosis and effective therapy of refractory cancer of unknown primary with PDL1 amplification (H021)

  Identification of the tissue of origin in cancer of unknown primary (CUP) poses a diagnostic challenge and is critical for directing site-specific therapy. Currently, clinical decision making in patients with CUP primarily relies on histopathology and clinical features. Comprehensive molecular profiling has the potential to contribute to diagnostic categorization and, most importantly, guide CUP therapy through identification of actionable lesions. We here report the case of an advanced-stage malignancy initially mimicking poorly differentiated soft-tissue sarcoma that did not respond to multi-agent chemotherapy. Molecular profiling within a clinical whole-exome and transcriptome sequencing program revealed a heterozygous, highly amplified KRAS G12S mutation, compound-heterozygous TP53 mutation/deletion, high mutational load, and focal amplification of chromosomes 9p (including PDL1 [CD274] and JAK2) and 10p (including GATA3). Integrated analysis of molecular data and conventional histopathology suggested a diagnosis of triple-negative breast cancer (TNBC) and provided a rationale for immune checkpoint inhibitor therapy with pembrolizumab, which resulted in rapid clinical improvement and a lasting partial remission. Analysis of 157 TNBC samples from The Cancer Genome Atlas revealed focal, high-level PDL1 amplification coinciding with excessive PDL1 mRNA expression in 24% of cases. Collectively, these results illustrate the diagnostic utility of multidimensional tumor profiling in cases with non-descript histology and immune phenotype, demonstrate the predictive power of genomic PDL1 amplification for immune checkpoint inhibition, and suggest a targeted therapeutic strategy in chromosome 9p24.1/PDL1-amplified cancers.H021

  Study EGAS00001001846

• Persistent STAG2 mutation in recurrent pediatric glioblastoma

  Similar to their adult counterparts, the prognosis for pediatric patients with high-grade gliomas remains poor. At time of recurrence, treatment options are limited and remain without consensus. This report describes the genetic findings, obtained from whole exome sequencing of a pediatric patient with glioblastoma who underwent multiple surgical resections and treatment with standard chemoradiation, as well as a novel recombinant poliovirus vaccine therapy. Strikingly, despite the variety of treatments, there was persistence of a tumor clone, characterized by a deleterious STAG2 mutation, whose deficiency in pre-clinical studies can cause aneuploidy and aberrant mitotic progression, but remains understudied in the clinical setting. There was near elimination of an EGFR mutated and amplified tumor clone after gross total resection, standard chemoradiation, and poliovirus therapy, followed by the emergence of a persistently STAG2 mutated clone, with rare mutations in PTPN11 and BRAF, the latter composed of a novel deleterious mutation previously not reported in pediatric glioblastoma (p.D594G). This was accompanied by a mutation signature shift towards one characterized by increased DNA damage repair defects, consistent with the known underlying STAG2 deficiency. As such, this case represents a novel report following the clinical and genetic progression of a STAG2 mutated glioblastoma, including treatment with a novel and emerging immunotherapy. Although STAG2 deficiency comprises only a small subset of gliomas, this case adds clinical evidence to existing pre-clinical data supporting a role for STAG2 mutations in tumor formation and resistance to standard therapies.

  Study EGAS00001004340

• Negligible impact on missing heritability of autoimmune-locus rare coding-region variants

  Genome wide association studies (GWAS) have identified common variants of modest effect size at hundreds of loci for common autoimmune diseases - however a substantial fraction of heritability remains unexplained, and to which rare variants may contribute. To discover rare variants and test them for association with a phenotype, the majority of studies re-sequence a small initial sample size and then genotype the discovered variants in a larger sample set. This approach will fail to analyse a large fraction of the rare variants present in the entire sample set. Here we perform simultaneous amplicon sequencing-based variant discovery and genotyping for coding exons of 25 GWAS risk genes in 41,911 white-European origin UK subjects comprising 24,892 subjects with six autoimmune disease phenotypes and 17,019 controls, and show that rare coding-region variants at known loci play a negligible role in common autoimmune disease susceptibility. These results do not support the rare variant - synthetic genome-wide association hypothesis. Many known autoimmune disease risk loci contain multiple independent common and low frequency variant association signals, and so genes in these loci are a priori stronger candidates to harbor rare coding-region variants than other genes. Our data suggest that the missing heritability for common autoimmune diseases may not lie in the rare coding-region variant portion of the allelic spectrum, but perhaps as others have proposed in very many common variant signals of weak effect.

  Study EGAS00001000476

• Massively parallel nanowell-based single-cell gene expression profiling

  Technological advances have enabled transcriptome characterization of cell types at the single-cell level providing new biological insights. New methods that enable simple yet high-throughput single-cell expression profiling are highly desirable. Here we report a novel nanowell-based single-cell RNA sequencing system, ICELL8, which enables processing of thousands of cells per sample. The system employs a 5184-nanowell-containing microchip to capture ~1300 single cells and process them. Each nanowell contains preprinted oligonucleotides encoding poly-d(T), a unique well barcode, and a unique molecular identifier. The ICELL8 system uses imaging software to identify nanowells containing viable single cells and only wells with single cells are processed into sequencing libraries. Here, we report the performance and utility of ICELL8 using samples of increasing complexity from cultured cells to mouse solid tissue samples. Our assessment of the system to discriminate between mixed human and mouse cells showed that ICELL8 has a low cell multiplet rate (< 3%) and low cross-cell contamination. We characterized single-cell transcriptomes of more than a thousand cultured human and mouse cells as well as 468 mouse pancreatic islets cells. We were able to identify distinct cell types in pancreatic islets, including alpha, beta, delta and gamma cells. Overall, ICELL8 provides efficient and cost-effective single-cell expression profiling of thousands of cells, allowing researchers to decipher single-cell transcriptomes within complex biological samples.

  Study EGAS00001002320