Search Results - EGA European Genome-Phenome Archive

Mutated H3 Histones Drive Human Pre-Leukemic Hematopoietic Stem Cell Expansion And Promote Leukemic Aggressiveness

Leukemogenesis is a stepwise progression from mutated, pre-neoplastic hematopoietic stem cells (HSCs) to full-blown leukemia. Our ability to prevent or treat acute myeloid leukemia (AML) is limited by our incomplete understanding of the epigenetic disruption that is central to this process, including improper histone methylation. Comprehensive analysis of 16 histone H3 genes in 434 primary AML samples identified Q69H, A26P, R2Q, R8H and K27M/I mutations (1.6%), and a higher incidence in secondary AML (s-AML) (9%). We establish that these mutations are important early events in leukemogenesis. They occur in pre-leukemic HSCs, increase the frequency of functional human HSCs, and alter differentiation. The mutations are present in the major leukemic clones in primary samples, and the mutant histones amplify leukemic aggressiveness with increased proliferation, expansion of leukemic progenitor and blast cells, and superior competitiveness in vivo. These effects are dependent on the specific mutation. Genome-wide analysis of K27 mutants revealed increased expression of genes involved in erythrocyte and myeloid differentiation with a corresponding decrease in histone H3 K27 tri-methylation and increase in K27 acetylation. The functional impact of histone mutations is independent of RUNX1 mutations, although they can co-occur. These data establish the involvement of H3 mutations as initial drivers of pre-cancerous stem cell expansion and leukemogenesis.

Study EGAS00001003288

Landscape of somatic mutations and clonal evolution in mantle cell lymphoma

Mantle cell lymphoma (MCL) is an aggressive tumor, but a subset of patients may follow an indolent clinical course. To understand the mechanisms underlying this biological heterogeneity, we performed whole-genome and/or whole-exome sequencing on 29 MCL cases and their respective matched normal DNA, as well as 6 MCL cell lines. Recurrently mutated genes were investigated by targeted sequencing in an independent cohort of 172 MCL patients. We identified 25 significantly mutated genes, including known drivers such as ataxia-telangectasia mutated (ATM), cyclin D1 (CCND1), and the tumor suppressor TP53; mutated genes encoding the anti-apoptotic protein BIRC3 and Toll-like receptor 2 (TLR2); and the chromatin modifiers WHSC1, MLL2, and MEF2B. We also found NOTCH2 mutations as an alternative phenomenon to NOTCH1 mutations in aggressive tumors with a dismal prognosis. Analysis of two simultaneous or subsequent MCL samples by whole-genome/whole-exome (n = 8) or targeted (n = 19) sequencing revealed subclonal heterogeneity at diagnosis in samples from different topographic sites and modulation of the initial mutational profile at the progression of the disease. Some mutations were predominantly clonal or subclonal, indicating an early or late event in tumor evolution, respectively. Our study identifies molecular mechanisms contributing to MCL pathogenesis and offers potential targets for therapeutic intervention.

Study EGAS00001000510

Using_genetics_to_identify_cell_types_and_mechanisms_underlying_susceptibility_to_primary_sclerosing_cholangitis

Primary sclerosing cholangitis (PSC) is a T-cell mediated, chronic inflammatory condition of the biliary tree that is strongly associated with inflammatory bowel disease. Genome-wide association studies have identified 22 non-HLA genetic risk variants associated PSC. Identifying the genes impacted by these variants has proven difficult as the majority lie in non-coding regions of the genome. Knowledge of the genes and biological pathways these non-coding variants are perturbing is vital to understanding the disease biology. One means of assessing the impact of non-coding variants within disease associated loci upon genes is via colocalisation with eQTL. Many eQTL are cell-type specific, requiring the analysis of disease relevant cell types to detect colocalisation. We have collected PSC-relevant T-cell-subtypes from the peripheral blood of PSC patients via fluorescence activated cell sorting in preparation for RNA sequencing and mapping of eQTL. Samples were collected at the Norfolk and Norwich University Hopital, for which local ethical approval has been granted. Lysed cell samples will be transferred to WTSI and DNA/RNA will be extracted from lysed cell samples by T143 before genotyping (DNA) and custom library preparation and sequencing (RNA). This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

Study EGAS00001002642

Chromothripsis in human breast cancer (HIPO K26K/H017/A017)

Chromothripsis is a form of genome instability, by which a presumably single catastrophic event generates extensive genomic rearrangements of one or a few chromosome(s). Widely assumed to be an early event in tumor development, this phenomenon plays a prominent role in tumor onset. We analyzed chromothripsis in 252 human breast cancers from two patient cohorts (149 metastatic breast cancers, 63 untreated primary tumors, 29 local relapses, 11 longitudinal pairs) using whole-genome and whole-exome sequencing. We showed that chromothripsis affects a substantial proportion of human breast cancers, with a prevalence over 60% in a cohort of metastatic cases and 25% in a cohort comprising predominantly luminal breast cancers (cohorts from HIPO K26K and H017 and A017). In the vast majority of cases, multiple chromosomes per tumor are affected, with most chromothriptic events on chromosomes 11 and 17 including, among other significantly altered drivers, CCND1, ERBB2, CDK12 and BRCA1. Importantly, chromothripsis generates recurrent fusions that drive tumor development. Chromothripsis-related rearrangements are linked with univocal mutational signatures, with clusters of point mutations due to kataegis in close proximity to the genomic breakpoints, and with the activation of specific signaling pathways. Analysis of the temporal order of events in tumors with and without chromothripsis as well as longitudinal analysis of chromothriptic patterns in tumor pairs revealed important insights on the role of chromothriptic chromosomes in tumor evolution.

Study EGAS00001004662

Kibbutzim Family study

The Kibbutzim Family Study (KFS) was established in 1992 to investigate the environmental and genetic determinants of cardiometabolic risk factors and their change over time. The participants belong to large families living in close-knit communities, called “Kibbutzim”, in Northern Israel. The Kibbutz has been a communal settlement, which has created a relatively homogeneous environment for its members. Kibbutz members are mostly of Ashkenazi Jewish ancestry, with the remaining members belonging to other Jewish subgroups. Participants were recruited in two phases from six Kibbutzim. In the first recruitment phase of the study (1992–1993) 500 individuals from 80 extended families (range 2 to 43) were examined. During the second phase (1999–2000), all participants from the first phase were invited for repeat examinations (80% response rate) and additional new participants were recruited, giving a total of 922 individuals from 150 extended families (range 2 to 55). Families were invited to participate if they consisted of at least four individuals who (i) lived in the Kibbutz, (ii) spanned at least two generations, and (iii) were at least 15 years old. Families were retained if at least two family members consented to participate in the study. Overall, 1033 participants were recruited; 111 were examined only in the first phase, 533 only in the second phase, and 389 were included in both. 901 individuals were successfully genotyped using Illumina HumanCoreExome BeadChip.

Study EGAS00001002782

Characterization_of_genomic_landscape_of_Peripheral_T_cell_Lymphomas__not_otherwise_specified__PTCL_NOS_

Peripheral T-cell lymphomas not otherwise specified (PTCL-NOS) represent a heterogeneous group of nodal and extra-nodal mature T-cell lymphomas, with a low prevalence in Western countries. PTCL-NOSs account for about 25% of all PTCLs and are currently diagnosed based on exclusion criteria, as this lymphomas lack unifying morphological, phenotypic and genomic features. Cytogenetic and FISH analysis of PTCL-NOS samples have not revealed recurrent pathogenetic abnormalities, while gene expression profiling has shown only partial ability to segregate cases representing homogeneous clinic-pathological entities. This underscores the need to look at PTCL-NOS with innovative and high-throughput approaches to identify recurrent genetic lesions that could further our understanding of the biology of this heterogeneous group of diseases, provide better diagnostic tools and perhaps new targets for innovative treatments. Our aim is to study ~15 patients affected by PTCL-NOS. Out study will be funded by a private, non-profit Italian cancer research fund (Associazione Italiana per la Ricerca sul Cancro, www.airc.it) based on a grant owned by Anna Dodero and Cristiana Carniti, hematologists at INT. Samples will be analysed by whole genome sequencing using Illumina X10 machines, on a 150bp-PE protocol. Data will be analysed using the pipeline available in Team 78, under the supervision of Peter Campbell, the WTSI faculty who will oversee the project, and by Francesco Maura, visiting scientist at the WTSI.

Study EGAS00001002057

Integrated analysis of relapsed B-cell precursor Acute Lymphoblastic Leukemia identifies subtype-specific cytokine and metabolic signatures

Recent efforts reclassified B-Cell Precursor Acute Lymphoblastic Leukemia (BCP-ALL) into more refined subtypes. Nevertheless, outcomes of relapsed BCP-ALL remain unsatisfactory, particularly in adult patients where the molecular basis of relapse is still poorly understood. To elucidate the evolution of relapse in BCP-ALL, we established a comprehensive multi-omics dataset including DNA-sequencing, RNA-sequencing, DNA methylation array and proteome MASS-spec data from matched diagnosis and relapse samples of BCP-ALL patients (n = 50) including the subtypes DUX4, Ph-like and two aneuploid subtypes. Relapse-specific alterations were enriched for chromatin modifiers, nucleotide and steroid metabolism including the novel candidates FPGS, AGBL and ZNF483. The proteome expression analysis unraveled deregulation of metabolic pathways at relapse including the key proteins G6PD, TKT, GPI and PGD. Moreover, we identified a novel relapse-specific gene signature specific for DUX4 BCP-ALL patients highlighting chemotaxis and cytokine environment as a possible driver event at relapse. This study presents novel insights at distinct molecular levels of relapsed BCP-ALL based on a comprehensive multi-omics integrated data set including a valuable proteomics data set. The relapse specific aberrations reveal metabolic signatures on genomic and proteomic levels in BCP-ALL relapse. Furthermore, the chemokine expression signature in DUX4 relapse underscores the distinct status of DUX4-fusion BCP-ALL.

Study EGAS00001002856

Mutational signatures of aflatoxin

Aflatoxin B1 (AFB1) is a mutagen and IARC (International Agency for Research on Cancer) Group 1 carcinogen that causes hepatocellular carcinoma (HCC). Here we present the first whole genome data on the mutational signatures of AFB1 exposure from a total of > 40,000 mutations in four experimental systems: two different human cell lines, and in liver tumors in wild-type mice and in mice that carried a hepatitis B surface antigen transgene – this to model the multiplicative effects of aflatoxin exposure and hepatitis B in causing HCC. AFB1 mutational signatures from all four experimental systems were remarkably similar. We integrated the experimental mutational signatures with data from newly-sequenced HCCs from Qidong County, China, a region of well-studied aflatoxin exposure. This indicated that COSMIC mutational signature 24, previously hypothesized to stem from aflatoxin exposure, indeed likely represents AFB1 exposure, possibly combined with other exposures. Among published somatic mutation data, we found evidence of AFB1 exposure in 0.7% of HCCs treated in North America, 1% of HCCs from Japan, but 16% of HCCs from Hong Kong. Thus, aflatoxin exposure apparently remains a substantial public health issue in some areas. This aspect of our study exemplifies the promise of future widespread resequencing of tumor genomes in providing new insights into the contribution of mutagenic exposures to cancer incidence.

Study EGAS00001002490

Exon4-mutations in KRAS Affect MEK-ERK and Pi3K-AKT Signaling in Multiple Myeloma

Approximately 40% of MM cases harbor a point mutation in either KRAS or NRAS. However, it is still unclear whether RAS-mutations are associated with disease outcome. Specifically, no information exists whether KRAS-mutations have an impact on survival in MM patients at diagnosis that were treated with novel agents (e.g. bortezomib, lenalidomide). Targeting KRAS directly for therapeutic purposes is still impossible, but recently the first covalent mutant-specific inhibitor targeting KRASG12C made it into clinical trials. However, other KRAS-hotspot-mutations exist in MM patients, including the less common exon4-mutations. Within the current study, the coding regions of KRAS in 80 MM patients at diagnosis, uniformly treated with bortezomib and high-dose chemotherapy, were deep-sequenced. Moreover, the functional impact of KRASG12A as well as the two exon4-mutations A146T (AMO1) and A146V (MM patient) on different survival pathways was investigated. Specifically, an siRNA-mediated knockdown in AMO1 was performed and KRASWT, KRASG12A, KRAS146T and KRAS146V were overexpressed in the KRASWT cell lines JJN3 and OPM2 using lentiviral transduction and the transposase based Sleeping Beauty system. KRAS-mutations were not associated with survival. Yet, the expression of all KRAS-mutants had an impact on MEK/ERK-signaling in two independent MM cell lines and on AKT/PI3K-signaling in one MM cell line.

Study EGAS00001004110

Extracellular transcriptome in seminal plasma for non-invasive prostate cancer diagnosis

A diagnostic non-invasive biomarker test for prostate cancer at an early stage, with high sensitivity and specificity, would improve diagnostic decision making. Extracellular RNAs present in seminal plasma might contain biomarker potential for the accurate detection of clinically significant prostate cancer. So far, the extracellular messenger RNA (mRNA) profile of seminal plasma has not been interrogated for its biomarker potential in the context of prostate cancer. Here, we investigate the mRNA transcriptome in seminal plasma samples obtained from prostate cancer patients (n=25), patients with benign prostate hyperplasia (n=26) and individuals without prostatic disease (n=6). Seminal plasma harbors a complex mRNA repertoire that reflects prostate as its tissue of origin. The endogenous RNA content is higher in the prostate cancer samples compared to the control samples. Prostate cancer antigen 3 (PCA3), a long non-coding RNA with prostate cancer-specific overexpression, and ATP-binding cassette transporter 1 (ABCA1), known to be involved in the prostate cancer pathogenesis, were more abundant in the prostate cancer group. In addition, twelve high confidence fusion transcripts could be detected in prostate cancer samples, including the bona-fide prostate cancer fusion transcript TMPRSS2-ERG. Our findings provide proof-of-principle that the extracellular transcriptome of seminal plasma can reveal information of an underlying prostate cancer.

Study EGAS00001005277

MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Oesophageal_adenocarcinoma

The Mutographs project aims to advance our understanding of the causes of cancer through studies of mutational signatures. Led by Mike Stratton, together with Paul Brennan, Ludmil Alexandrov, Allan Balmain, David Phillips and Peter Campbell, this large-scale international research endeavour was awarded a Cancer Research UK Grand Challenge. Different patterns of somatic mutation are generated by the different environmental, lifestyle and genetic factors that cause cancer, many of them are still unknown. Within Mutographs, the International Agency for Research on Cancer is coordinating the recruitment of 5000 individuals with cancer (colorectal, renal, pancreatic, oesophageal adenocarcinoma or oesophageal squamous cancers) across 5 continents to explore whether different mutational signatures explain marked variation in incidence. In brief, through an international network of collaborators around the world, biological materials are collected, along with demographic, histological, clinical and questionnaire data. Whole genome sequences of tumour-germline DNA pairs are generated at the Wellcome Trust Sanger Institute (40X and 20X depth respectively). Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

Study EGAS00001003702

Small_molecule_inhibitors_in_melanoma_Kenski_Kong___WES

Small-molecule inhibitors targeting the most commonly activated pathway in melanoma, MAPK pathway (either alone or in combination) are already given to melanoma patients for few years, and initially reduce tumour burden dramatically, eventually melanomas become resistant and tumours progress while on treatment. Resistance to this treatment occurs by acquisition of additional mutations or other alterations that affect the mitogen-activated protein kinase (MAPK) pathway by either direct or indirect signalling. Many resistance mechanisms somehow lead to reactivation of extracellular signal-regulated kinase (ERK), thereby restoring signalling of the oncogenic BRAF/MEK/ERK pathway. In addition, PI3K pathway activation contributes to resistance to BRAF inhibition. Less frequent but equally important to the phenomenon of targeted drug resistance is the observation that B15–20% of BRAF mutant melanoma patients fail to respond to BRAF inhibition already early on treatment, owing to intrinsic resistance. These patients have little therapeutic options, unless immunotherapy can be given. To better understand the resistance mechanisms in MAPK inhibitor-treated melanoma patients and melanoma biology, our lab generated a big panel of MAPK inhibitor resistant melanoma cell lines by continuous drug exposure. The understanding of the genetic landscape and gene expression as well as cross resistance to other treatment regimens, and other aspects of melanoma biology such as phenotype switch, will allow us to better exploit new therapeutic strategies for melanoma patients.

Study EGAS00001002863

Transcriptional and DNA binding profiling of CEBPE in REH acute lymphoblastic leukaemia cells using shRNA RNA-Seq and ChIP-Seq

ChIP-Seq - CEBPE - REH. The ETV6/RUNX1 translocated acute lymphoblastic leukaemia cell line REH was used to perform ChIP-Seq using a CEBPE antibody. Cells were fixed in 1% formaldehyde for 10mins, prior to preparation of chromatin using Active Motif Express ChIP-IT. 2ug of antibody (anti CEBPE Atlas Antibodies HPA002928)was added to 25ug of chromatin O/N at 4C with rotation. Duplicate reactions were pooled and purified. 10ng of ChIPd and input DNA used for Illumina NGS preparation (NEBNext ChIP-Seq Library kit; New England Biolabs), CEBPE and Input DNA ChIP samples were sequenced on a MiSeq using 150bp Kit v3 paired end and a HiSeq 2500 using 2x101 version 4 paired end (Illumina) respectively. Reactions performed in duplicate. shCEBPE RNA-Seq - REH. REH cells were lentivirally transduced with a pTRIPZ shRNA vector for transcriptional profiling of CEBPE. Two controls (empty and non-targeting) and two CEBPE shRNAs (V3THS_150517(A13), V3THS_404312(G3) Dharmacon, GE) were transduced into REH cells. Cells were treated with 1ug/ml doxycyclin for 144hrs and total RNA purified. Libraries were prepared using NEBNext Ultra II Directional RNA Library Prep Kit and sequenced on an Illuimna HiSeq 2500 using 2x101 version 4 paired end chemistry. 3 biological replicates of each samples were prepared.

Study EGAS00001002877

EXOME-WIDE ASSOCIATION ANALYSIS OF CORONARY ARTERY DISEASE IN THE KINGDOM OF SAUDI ARABIA POPULATION

Coronary Artery Disease (CAD) remains the leading cause of mortality worldwide. Mortality rates associated with CAD have shown an exceptional increase particularly in fast developing economies like the Kingdom of Saudi Arabia (KSA). Over the past twenty years, CAD has become the leading cause of death in KSA and has reached epidemic proportions. This rise is undoubtedly caused by fast urbanization that is associated with a life-style that promotes CAD. However, the question remains whether genetics play a significant role and whether genetic susceptibility is increased in KSA compared to the well-studied Western European populations. Therefore, we performed an Exome-wide association study (EWAS) in 832 patients and 1,076 controls of Saudi Arabian origin to test whether population specific, strong genetic risk factors for CAD exist, or whether the polygenic risk score for known genetic risk factors for CAD, lipids, and Type 2 Diabetes show evidence for an enriched genetic burden. Our results do not show significant associations for a single genetic locus. However, the heritability estimate for CAD for this population was high (h2 = 0.53, S.E. = 0.1, p = 4e-12) and we observed a significant association of the polygenic risk score for CAD that demonstrates that the population of KSA, at least in part, shares the genetic risk associated to CAD in Western populations.

Study EGAS00001001645

scRNA-seq data of human nuclei collected from the temporal cortex of 17 individuals.

Transcriptome analysis of single cells has tremendously increased our understanding of cellular heterogeneity in complex tissues. However, whole-cell transcriptomics cannot be applied to frozen tissue due to poor cell recovery after freezing. Human tissue is usually limited to tissue banks containing frozen material, and since single nuclei can be readily isolated from frozen tissues, human cell nuclei can substitute whole-cell analysis to extract single cell transcriptomes from human cells. Despite the importance of human tissue analysis, single-nucleus RNA-sequencing protocols are less developed than those for whole cells. SMARTSeq, either in form of custom-made SMARTSeq2 or as part of the Fluidigm C1 platform, provides a high resolution view of the transcriptome, has proven effective for whole-cell analysis and has recently been adapted to nuclei. Here, we show that the standard SMARTSeq2 protocol applied to single nuclei biases the amplification of single-nucleus transcriptomes due to unspecific capturing of transcripts, which can comprise the majority of reads in single-nucleus RNA-sequencing libraries. We demonstrate that modification of the SMARTSeq2 protocol with respect to template switching and cDNA amplification conditions generates sequencing libraries of higher complexity, allows for detection of more genes and reduces biases of the amplified transcriptome. In addition, modification of SMARTSeq2 chemistry allows for gene detection at shallower sequencing depth, thus increasing efficiency and decreasing cost.

Study EGAS00001002882

Glioblastoma epigenome profiling identifies SOX10 as a master regulator of molecular tumour subtype

Glioblastoma in adult patients are classified into four subtypes, IDH, MES, RTK I, and RTK II, based on DNA methylation and RNA expression data. Tumour subtype transitions commonly occur during the treatment of glioblastoma patients, and transitions to the mesenchymal (MES) subtype have been associated with therapy resistance and adverse prognosis. Here, we generate DNA methylome and histone modification data of glioblastoma primary tumours and show that glioblastoma subtypes differ by their enhancer landscapes. Using Core Regulatory Circuitry analysis of H3K27ac data, and independent analysis of RNA-derived gene regulatory networks, we identified 38 subtype master regulators whose cell type-specific activities we mapped in single-cell RNA sequencing data. These analyses identified the chromatin modifier SRY-Box 10 (SOX10) as a master regulator in tumours of the RTK I subtype. In functional studies we demonstrate that SOX10 loss causes a transition of RTK I tumour cells to a mesenchymal cellular state by altering the accessibility of enhancer-rich chromatin regions and their occupancy by Bromodomain Containing 4 (BRD4), a dynamic constituent of super-enhancers. Treatment with the BRD4 inhibitor JQ1 blocked this transition. These data demonstrate the fundamental role of enhancer chromatin remodelling and master regulator activity in the establishment, maintenance and plasticity of glioblastoma cellular states. We also make our data publically browsable at https://dkfz-b060.github.io/gb_browse/

Study EGAS00001003953

Deep multi-region whole-genome sequencing reveals heterogeneity and gene-by-environment interactions in treatment-naive, metastatic lung cancer

Our understanding of genomic heterogeneity in lung cancer is largely based on the analysis of early-stage surgical specimens. Here we used endoscopic sampling of paired primary and intrathoracic metastatic tumors from 11 lung cancer patients to map genomic heterogeneity inoperable lung cancer with deep whole-genome sequencing. Intra-patient heterogeneity in driver or targetable mutations was predominantly in the form of copy number gain. Private mutation signatures, including patterns consistent with defects in homologous recombination, were highly variable both within and between patients. Irrespective of histotype, we observed a smaller than expected number of private mutations, suggesting that ancestral clones accumulated large mutation burdens immediately prior to metastasis. Single-region whole-genome sequencing of from 20 patients showed that tumors in ever-smokers with the strongest tobacco signatures were associated with germline variants in genes implicated in the repair of cigarette-induced DNA damage. Our results suggest that lung cancer precursors in ever-smokers accumulate large numbers of mutations prior to the formation of frank malignancy followed by rapid metastatic spread. In advanced lung cancer, germline variants in DNA repair genes may interact with the airway environment to influence the pattern of founder mutations, whereas similar interactions with the tumor microenvironment may play a role in the acquisition of mutations following metastasis.

Study EGAS00001003830

Spatial concordance of DNA methylation classification in diffuse glioma

Background: Intratumoral heterogeneity is a hallmark of diffuse gliomas. DNA methylation profiling is an emerging approach in the clinical classification of brain tumors. The goal of this study is to investigate the effects of intratumoral heterogeneity on classification confidence. Methods: We used neuronavigation to acquire 133 image-guided and spatially-separated stereotactic biopsy samples from 16 adult patients with a diffuse glioma (7 IDH-wildtype and 2 IDH-mutant glioblastoma, 6 diffuse astrocytoma, IDH-mutant and 1 oligodendroglioma, IDH-mutant and 1p19q codeleted), which we characterized using DNA methylation arrays. Samples were obtained from regions with and without abnormalities on contrast enhanced T1 weighted and fluid-attenuated inversion recovery MRI. Methylation profiles were analyzed to devise a three-dimensional reconstruction of (epi)genetic heterogeneity. Tumor purity was assessed from clonal methylation sites. Results: Molecular aberrations indicated that tumor was found outside imaging abnormalities, underlining the infiltrative nature of this tumor and the limitations of current routine imaging modalities. We demonstrate that tumor purity is highly variable between samples and explains a substantial part of apparent epigenetic spatial heterogeneity. We observed that DNA methylation subtypes are often, but not always, conserved in space taking tumor purity and prediction accuracy into account. Conclusion: Our results underscore the infiltrative nature of diffuse gliomas and suggest that DNA methylation subtypes are relatively concordant in this tumor type, although some heterogeneity exists.

Study EGAS00001005434

Bulk-tissue RNA-sequencing of anterior cingulate cortex samples derived from Lewy body disease patients

Parkinson’s disease (PD), Parkinson’s disease with dementia (PDD) and dementia with Lewy bodies (DLB) are three clinically, genetically and neuropathologically overlapping neurodegenerative diseases collectively known as the Lewy body diseases (LBDs). A variety of molecular mechanisms have been implicated in PD pathogenesis, but the mechanisms underlying PDD and DLB remain largely unknown, a knowledge gap that presents an impediment to the discovery of disease-modifying therapies. Transcriptomic profiling can contribute to addressing this gap, but remains limited in the LBDs. Here, we applied paired bulk-tissue and single-nucleus RNA-sequencing to anterior cingulate cortex samples derived from 28 individuals, including healthy controls, PD, PDD and DLB cases (n = 7 per group), to transcriptomically profile the LBDs. Using this approach, we (i) found transcriptional alterations in multiple cell types across the LBDs; (ii) discovered evidence for widespread dysregulation of RNA splicing, particularly in PDD and DLB; (iii) identified potential splicing factors, with links to other dementia-related neurodegenerative diseases, coordinating this dysregulation; and (iv) identified transcriptomic commonalities and distinctions between the LBDs that inform understanding of the relationships between these three clinical disorders. Together, these findings have important implications for the design of RNA-targeted therapies for these diseases and highlight a potential molecular “window” of therapeutic opportunity between the initial onset of PD and subsequent development of Lewy body dementia.

Study EGAS00001005305

Unraveling the transcriptomic profile of utero-tubal lavage fluid of ovarian cancer patients

The majority of ovarian cancers are diagnosed at an advanced stage and have a high mortality rate. Current screening strategies fail to improve prognosis because markers of early stage disease are lacking. This medical need justifies the search for biomarkers enabling early detection of ovarian cancer. In this study, we explore the extracellular transcriptome of utero-tubal lavage fluid obtained from 26 ovarian cancer patients and 48 controls using messenger RNA (mRNA) capture and small RNA sequencing. Because of its contact with the fallopian tube and the ovaries, collection of utero-tubal lavage might be the preferred non-invasive method to identify localized lesions compared with analyses based on peripheral blood samples. We observed an enrichment of ovarian and fallopian tube specific messenger RNAs in utero-tubal lavage fluid compared to other human biofluids. Over 300 mRNAs and 41 miRNAs were upregulated in ovarian cancer samples compared with controls. Upregulated genes were enriched for genes involved in cell cycle activation and proliferation, hinting at a tumor derived signal. In conclusion, we provide proof-of-principle that mRNA capture sequencing of utero-tubal lavage fluid is technically feasible, and that the extracellular transcriptome of utero-tubal lavage should be further explored in larger cohorts to assess the diagnostic value of the biomarkers identified in this study.

Study EGAS00001005498

SG Peranakan Project: Genetic admixture in the culturally unique Peranakan Chinese population in Southeast Asia

The Peranakan Chinese are culturally unique descendants of immigrants from China who settled in the Malay Archipelago ~300-500 years ago. Today, among large communities in Southeast Asia, the Peranakans have preserved Chinese traditions with strong influence from the local indigenous Malays. Yet, whether or to what extent genetic admixture co-occurred with the cultural mixture has been a topic of ongoing debate. We performed whole-genome sequencing (WGS) on 177 Singapore (SG) Peranakans and analyzed the data jointly with WGS data of Asian and European populations. We estimated that Peranakan Chinese inherited ~5.62% (95% confidence interval [CI]: 4.75-6.46%) Malay ancestry, much higher than that in SG Chinese (1.08%, 0.69-1.53%), southern Chinese (0.86%, 0.57-1.31%), and northern Chinese (0.25%, 0.18-0.33%). A sex-biased admixture history, in which the Malay ancestry was contributed primarily by females, was supported by X chromosomal variants, and mitochondrial (MT) and Y haplogroups. Finally, we identified an ancient admixture event shared by Peranakan Chinese and SG Chinese ~1,612 (95% CI: 1,345-1,923) years ago, coinciding with the settlement history of Han Chinese in southern China, apart from the recent admixture event with Malays unique to Peranakan Chinese ~190 (159-213) years ago. These findings greatly advance our understanding of the dispersal history of Chinese and their interaction with indigenous populations in Southeast Asia.

Study EGAS00001005379

Sequence Data from the phase 2 PrE0505 trial of Durvalumab with First Line Platinum-Pemetrexed for Unresectable Pleural Mesothelioma

Mesothelioma is a rare and fatal cancer with limited therapeutic options until the recent approval of combination immune checkpoint blockade. We report the results of the phase 2 PrE0505 trial (NCT02899195) of the anti-PD-L1 antibody, durvalumab, plus platinum-based and pemetrexed chemotherapy for patients with previously untreated unresectable pleural mesothelioma. The primary endpoint was overall survival compared to historic control with cisplatin and pemetrexed chemotherapy. The combination of durvalumab with chemotherapy met the pre-specified primary endpoint reaching a median survival of 20.4 months vs. 12.1 months with historic control. Treatment emergent adverse events were consistent with known side effects of chemotherapy and all adverse events due to immunotherapy were grade ≤2. Integrated genomic and immune cell repertoire analyses revealed that a higher immunogenic mutation burden coupled with a more diverse T cell repertoire were linked with favorable clinical outcome. Structural genome-wide analyses demonstrated a higher degree of genomic instability in responding tumors of epithelioid histology. Patients with germline alterations in cancer predisposing genes, especially those involved in DNA repair, were more likely to attain long term survival. Our findings indicate that concurrent durvalumab with platinum-based chemotherapy has promising clinical activity and that responses are driven by the complex genomic background of malignant pleural mesothelioma.

Study EGAS00001005426

MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancerhead_neck

The Mutographs project aims to advance our understanding of the causes of cancer through studies of mutational signatures. Led by Mike Stratton, together with Paul Brennan, Ludmil Alexandrov, Allan Balmain, David Phillips and Peter Campbell, this large-scale international research endeavour was awarded a Cancer Research UK Grand Challenge. Different patterns of somatic mutation are generated by the different environmental, lifestyle and genetic factors that cause cancer, many of them are still unknown. Within Mutographs, the International Agency for Research on Cancer is coordinating the recruitment of 5000 individuals with cancer (colorectal, renal, pancreatic, oesophageal adenocarcinoma or oesophageal squamous cancers) across 5 continents to explore whether different mutational signatures explain marked variation in incidence. In brief, through an international network of collaborators around the world, biological materials are collected, along with demographic, histological, clinical and questionnaire data. Whole genome sequences of tumour-germline DNA pairs are generated at the Wellcome Trust Sanger Institute. Non-tumour DNA is also studied for some cancer types. Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

Study EGAS00001005450

Comprehensive Genomic Characterization of Gene Therapy-Induced T-cell Acute Lymphoblastic Leukemia (H007)

Six out of ten patients in a hematopoietic stem cell gene therapy (GT) trial for Wiskott-Aldrich syndrome (WAS) developed T-cell acute lymphoblastic leukemia (T-ALL) associated with gamma-retroviral vector integration at the LMO2 locus. We hypothesized that detailed molecular analysis of T-ALL driven by temporally defined first deregulation of LMO2 might offer valuable insights into the pathogenesis of GT-associated and sporadic T-ALL (HIPO project 007). In all six patients whole-genome and transcriptome sequencing demonstrated that in addition to oncogenic chromosomal translocations involving the T-cell receptor loci, retrovirally induced leukemias faithfully recapitulated secondary genomic events observed in sporadic T-ALL. Multiple genetic events culminated in the development of a rapidly proliferating leukemic clone and onset of clinical manifestation at 1.3 to 5 years after infusion of the gene corrected transplant. These genetic alterations affected NOTCH1 and PI3K signaling, cell cycle regulation as wells as T-cell transcription factors LEF1 and TCF1. Two of the six cases displayed a phenotype reminiscent of early immature T-ALLs, and lacked additional T-ALL typical alterations. Our analysis shows that insertional mutagenesis involving the LMO2 locus triggers T-ALL development due to diverse genetic alterations in large parts, but not exclusively, in known T-ALL driver genes, demonstrating the molecular diversity of GT-associated T-ALL.

Study EGAS00001003870

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

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 EGAS00001005584

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

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

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

MutWP1__CRUK_Grand_Challenge__normal_kidney_Nanoseq

The Mutographs project aims to advance our understanding of the causes of cancer through studies of mutational signatures. Led by Mike Stratton, together with Paul Brennan, Ludmil Alexandrov, Allan Balmain, David Phillips and Peter Campbell, this large-scale international research endeavour was awarded a Cancer Research UK Grand Challenge. Different patterns of somatic mutation are generated by the different environmental, lifestyle and genetic factors that cause cancer, many of them are still unknown. Within Mutographs, the International Agency for Research on Cancer is coordinating the recruitment of 5000 individuals with cancer (colorectal, renal, pancreatic, oesophageal adenocarcinoma or oesophageal squamous cancers) across 5 continents to explore whether different mutational signatures explain marked variation in incidence. In brief, through an international network of collaborators around the world, biological materials are collected, along with demographic, histological, clinical and questionnaire data. Whole genome sequences of tumour-germline DNA pairs are generated at the Wellcome Trust Sanger Institute. Non-tumour DNA is also studied for some cancer types. Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

Study EGAS00001005451

MutWP1__CRUK_Grand_Challenge__matched_blood_Nanoseq

The Mutographs project aims to advance our understanding of the causes of cancer through studies of mutational signatures. Led by Mike Stratton, together with Paul Brennan, Ludmil Alexandrov, Allan Balmain, David Phillips and Peter Campbell, this large-scale international research endeavour was awarded a Cancer Research UK Grand Challenge. Different patterns of somatic mutation are generated by the different environmental, lifestyle and genetic factors that cause cancer, many of them are still unknown. Within Mutographs, the International Agency for Research on Cancer is coordinating the recruitment of 5000 individuals with cancer (colorectal, renal, pancreatic, oesophageal adenocarcinoma or oesophageal squamous cancers) across 5 continents to explore whether different mutational signatures explain marked variation in incidence. In brief, through an international network of collaborators around the world, biological materials are collected, along with demographic, histological, clinical and questionnaire data. Whole genome sequences of tumour-germline DNA pairs are generated at the Wellcome Trust Sanger Institute. Non-tumour DNA is also studied for some cancer types. Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

Study EGAS00001006774

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

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

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

Target gene sequencing for human normal endometrial glands

The raw fastq files target sequencing of 112 genes for 1,298 endometrial glands and matched blood samples. The paired-end sequencing data sets (R1 and R2) are deposited. ABCC1, ACRC, ANK3, ARHGAP35, ARID1A, ARID5B, ATCAY, ATM, ATR, BARD1, BCOR, BRCA1, BRCA2, BRD4, BRIP1, CAMTA1, CDC23, CDYL, CFAP54, CHD4, CHEK1, CHEK2, CTCF, CTNNB1, CUX1, DGKA, DISP2, DYNC2H1, EMSY, FAAP24, FAM135B, FAM175A, FAM65C, FANCA, FANCB, FANCC, FANCD2, FANCE, FANCF, FANCG, FANCI, FANCL, FANCM, FAT1, FAT3, FBN2, FBXW7, FGFR2, FRG1, GPR50, HEATR1, HIST1H4B, HNRNPCL1, HOOK3, KIAA1109, KIF26A, KMT2B, KMT2C, KRAS, LAMA2, LRP1B, MLH1, MON2, MRE11A, MSH2, MSH6, MTOR, NBN, PALB2, PHEX, PIK3CA, PIK3R1, PLXNB2, PLXND1, PMS2, POLE, POLR3B, PPP2R1A, PTEN, PTPN13, RAD50, RAD51, RAD51B, RAD51C, RAD51D, RAD52, RAD54B, RAD54L, RICTOR, SACS, SIGLEC9, SLC19A1, SLX4, SPEG, STT3A, TAF1, TAF2, TAS2R31, TFAP2C, TNC, TONSL, TP53, TTC6, UBA7, VNN1, WT1, XIRP2, ZBED6, ZC3H13, ZFHX3, ZFHX4, ZMYM4.

Dataset EGAD00001008465

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

A prognostic human brain network for diffuse midline glioma

Diffuse midline glioma (DMG) are near-universally lethal tumours of the paediatric central nervous system. In animal models, DMG form brain-wide, integrated networks through neuron-to-glioma synapses and glioma-to-glioma gap junctional coupling. This extensive connectivity robustly promotes the growth and invasion of DMG and other glial malignancies through paracrine mechanisms and direct neuron-to-glioma synapses. The organisation and clinical implications of these connections in the living human brain, however, remain to be elucidated. Here, we develop tumour network mapping to compute the brain-wide connectivity profile of DMG: defining a conserved brain network across pontine and thalamic DMG associated with patient short-term survival (DMG network). Tumour functional connectivity with the DMG network was independently predictive of patient overall survival across two external validation cohorts. Tumour growth mapped to DMG network-specific trajectories and peak, in-network neurometabolic changes across development spatiotemporally aligned with the peak age incidence of DMG. Single-nucleus RNA-sequencing confirmed diverse synaptic gene enrichment in high-connectivity DMG. Strikingly, incidental surgical resection of high-connectivity thalamic DMG tissue conferred a significant survival advantage. Collectively, these data define a conserved and prognostically important brain network in children with DMG, consistent with the hypothesis that DMG exploit otherwise healthy brain circuits to promote tumour growth.

Study EGAS50000001752

PAX5 biallelic genomic alterations define a novel subgroup of B cell precursor acute lymphoblastic leukemia

Chromosomal rearrangements and specific aneuploidy patterns are initiating events and define subgroups in B cell precursor acute lymphoblastic leukemia (BCP-ALL). Here we analyzed 250 BCP-ALL cases and identified a novel subgroup (‘PAX5-plus’, n=19) by distinct DNA methylation and gene expression profiles. All patients in this subgroup harbored mutations in the B-lineage transcription factor PAX5, with p.P80R as hotspot. Mutations either affected two independent codons, consistent with compound heterozygosity, or suffered LOH predominantly through chromosome 9p aberrations. These biallelic events resulted in disruption of PAX5 transcriptional programs regulating B cell differentiation and tumor suppressor functions. Homozygous CDKN2A/B deletions and RAS activating hotspot mutations were highly enriched as cooperating events in the genomic profile of PAX5-plus ALL. Together this defined a specific pattern of triple alterations, exclusive to the novel subgroup. PAX5-plus ALL was observed in pediatric and adult patients. Although restricted by the limited sample size, a tendency for more favorable clinical outcome was observed with 10 of 12 adult PAX5-plus patients achieving long-term survival. PAX5-plus represents the first BCP-ALL subgroup defined by sequence alterations in contrast to gross chromosomal events and exemplifies how deregulated differentiation (PAX5), impaired cell cycle control (CDKN2A/B) together with sustained proliferative signaling (RAS) cooperatively drives leukemogenesis.

Study EGAS00001003209

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

Prognostic relevance of microenvironmental factors CD163 and CD8 combined with EZH2 and chromosome 18 gain in a validation cohort of follicular lymphoma patients of the Lunenburg Lymphoma Biomarker Consortium

Background: In follicular lymphoma (FL), studies addressing the prognostic value of microenvironment-related immunohistochemical (IHC) markers and tumor cell-related genetic markers have yielded conflicting results, precluding implementation in practice. Therefore, the Lunenburg Lymphoma Biomarker Consortium (LLBC) performed a validation study for published markers. Methods: To maximize sensitivity, an end-of-spectrum design was applied for 122 uniformly immunochemotherapy-treated FL patients retrieved from international trials and registries; early failure (EF): progression or lymphoma-related death <2 years versus long remission: response duration of >5 years. IHC staining for T-cells and macrophages was performed on tissue microarrays from initial biopsy and scored with a validated computer-assisted protocol. Shallow whole-genome and deep targeted sequencing was performed on the same samples. Results: 96/122 cases with complete molecular and immunohistochemical data were included in the analysis. EZH2 wild-type (p=0.006), gain of chromosome 18 (p=0.002), low percentages of CD8+ cells (p=0.011) and CD163+ areas (p=0.038) were associated with EF. No significant differences in other markers were observed, thereby refuting previous claims on their prognostic significance. Conclusion: Using an optimized study design, this LLBC study validates wild-type EZH2 status, gain of chromosome 18, low percentages of CD8+ cells and CD163+ area as predictors of EF to immunochemotherapy in FL.

Study EGAS00001002049

The double-hit signature identifies double-hit diffuse large B-cell lymphoma with genetic events cryptic to FISH

High-grade B-cell lymphomas with MYC and BCL2 and/or BCL6 rearrangements (HGBL-DH/THs) include a group of diffuse large B-cell lymphomas (DLBCLs) with inferior outcomes after standard chemoimmunotherapy. We recently described a gene expression signature that identifies 27% of germinal center B-cell DLBCLs (GCB-DLBCLs) as having a double-hit-like expression pattern (DHITsig) and inferior outcomes; however, only half of these cases have both MYC and BCL2 translocations identifiable using standard breakapart fluorescence in situ hybridization (FISH). Here, 20 DHITsig+ GCB-DLBCLs apparently lacking MYC and/or BCL2 rearrangements underwent whole-genome sequencing. This revealed 6 tumors with MYC or BCL2 rearrangements that were cryptic to breakapart FISH. Copy-number analysis identified 3 tumors with MYC and 6 tumors with MIR17HG gains or amplifications, both of which may contribute to dysregulation of MYC and its downstream pathways. Focal deletions of the PVT1 promoter were observed exclusively among DHITsig+ tumors lacking MYC translocations; this may also contribute to MYC overexpression. These results highlight that FISH fails to identify all HGBL-DH/THs, while revealing a range of other genetic mechanisms potentially underlying MYC dysregulation in DHITsig+ DLBCL, suggesting that gene expression profiling is more sensitive for identifying the biology underlying poor outcomes in GCB-DLBCL.

Study EGAS00001004285

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

Inquiring the potential of donor derived CD19-CAR TSCM to treat B-cell malignancies.

Relapse after allogeneic hematopoietic stem cell transplantation (alloHSCT) is associated with poor survival in advanced B-cell malignancies. Donor-derived CD19-CAR T cells offer the potential for durable remissions but are limited by poor engraftment, expansion, and persistence. We report results from a sequential Phase 1 study (NCT01087294) in which patients with B-cell malignancies relapsed after alloHSCT received either donor-derived standard or stem cell memory T (TSCM)-enriched CD19-CAR T cells without lymphodepletion preconditioning. CAR TSCM cells exhibited delayed but superior expansion and persistence, leading to higher complete remission rates (45% vs. 10%) for doses < 3×108 CAR T cells. Despite their potency, CAR TSCM cells induced milder cytokine-mediated toxicities, driven by IFN-γ rather than IL-6, as observed in standard CAR T cell recipients. While both T cell products differentiated into effectors, only CAR TSCM cells robustly reconstituted the stem-like compartment. Longitudinal clonal-tracking revealed CAR TSCM cells were maintained by clonal succession, whereas standard CAR T cells persisted through maintenance or attrition of early-expanded clones. Unlike standard CAR T cells, CAR TSCM treatment failure was primarily due to tumor- and host-related factors rather than intrinsic T cell defects. These findings establish TSCM cells as a safe and highly effective CAR T cell platform.

Study EGAS00001008119

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

23249 results for "ega"

in 81.74 milliseconds.