Search Results - EGA European Genome-Phenome Archive

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

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

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

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

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

23014 results for "ega"

in 59.33 milliseconds.