23249 results for "ega"

in 90.41 milliseconds.

• Systematic dissection of tumor-normal single-cell ecosystems across a thousand tumors of 30 cancer types

  The complexity of the tumor microenvironment poses significant challenges in cancer therapy. Here, to comprehensively investigate the tumor-normal ecosystems, we perform an integrative analysis of 4.9 million single-cell transcriptomes from 1,070 tumor and 493 normal samples in combination with pan-cancer 137 spatial transcriptomics, 8,887 TCGA, and 1,261 checkpoint inhibitor-treated bulk tumors. We define a myriad of cell states constituting the tumor-normal ecosystems and also identify hallmark gene signatures across different cell types and organs. Our atlas characterizes distinctions between inflammatory fibroblasts marked by AKR1C1 or WNT5A in terms of cellular interactions and spatial co-localization patterns. Co-occurrence analysis reveals interferon-enriched community states including tertiary lymphoid structure (TLS) components, which exhibit differential rewiring between tumor, adjacent normal, and healthy normal tissues. The favorable response of interferon-enriched community states to immunotherapy is validated using immunotherapy-treated cancers (n=1,261) including our lung cancer cohort (n=497). Deconvolution of spatial transcriptomes discriminates TLS-enriched from non-enriched cell types among immunotherapy-favorable components. Our systematic dissection of tumor-normal ecosystems provides a deeper understanding of inter- and intra-tumoral heterogeneity.

  Study EGAS50000000324

• Prediction of homologous recombination deficiency identifies colorectal tumors sensitive to PARP inhibition

  BRCAness is a well-established feature in breast, ovarian, prostate and pancreatic carcinomas, our recent findings indicate that up to 15% of colorectal cancers (CRC) also harbor defects in the HR pathway, presenting promising opportunities for innovative therapeutic strategies in CRC patients. We developed a new tool called HRDirect, which builds upon the HRDetect algorithm and is able to predict HRD from reference-free tumor samples. We validated HRDirect using matched breast cancer and CRC patient samples. Subsequently, we assessed its efficacy in predicting response to the PARP inhibitor olaparib by comparing it with two other commercial assays: AmoyDx HRD by Amoy Diagnostics and the TruSight Oncology (TSO) _500 HRD panel by Illumina NGS technology. While all three approaches successfully identified the most PARPi-sensitive CRC models, HRDirect demonstrated superior precision in distinguishing resistant models compared to AmoyDX and TSO500-HRD, which exhibited overlapping scores between sensitive and resistant cells. Furthermore, we propose integrating HRDirect scoring with ATM immunohistochemical analysis as part of our "composite biomarker approach" to enhance the identification of HRD tumors, with an immediate translational and clinical impact for CRC personalized treatment.

  Study EGAS50000000426

• Harnessing Epigenetic Regulators to improve HSC-based lentiviral gene therapy

  Curative benefits of autologous and allogeneic transplantation of hematopoietic stem cells (HSCs) have been proven for various diseases. However, the low number of true HSCs that can be collected from patients and subsequently in vitro maintenance and expansion of true HSCs for genetic correction remain challenging. Addressing this issue, we here focused on optimizing culture conditions to improve the ex vivo expansion of true HSCs for gene therapy purposes. In particular, we explore the use of epigenetic regulators to enhance the effectiveness of HSC-based lentiviral (LV) gene therapy. The HDAC inhibitor Quisinostat and the bromodomain inhibitor CPI203 each promote ex vivo expansion of functional HSCs, as validated by xenotransplantation assays and single cell RNA-sequencing analysis. We confirmed the stealth effect of LV transduction on the loss of HSC numbers in commonly used culture protocols, while addition of Quisinostat or CPI203 improved expansion of HSCs in transduction protocols. Of note, we demonstrated that addition of Quisinostat improved LV transduction efficiency of HSCs and early progenitors. Our suggested culture conditions highlight the potential therapeutic effect of epigenetic regulators in hematopoietic stem cell biology and their clinical applications to advance HSC-based gene correction.

  Study EGAS50000000175

• Identification of HER2-positive breast cancer molecular subtypes with potential clinical implications in the ALTTO clinical trial

  In HER2-positive breast cancer, clinical outcome and sensitivity to HER2-targeted therapies are influenced by both tumor and microenvironment features. However, we are currently unable to depict the molecular heterogeneity of this disease with sufficient granularity. Here, by performing gene expression profiling in HER2-positive breast cancers from patients receiving adjuvant trastuzumab in the ALTTO clinical trial, we identified and characterized five molecular subtypes associated with the risk of distant recurrence: immune-enriched (IM), proliferative/metabolic-enriched (P/Met), mesenchymal/stroma-enriched (Mes/S), luminal (LUM), and ERBB2-enrichedERBB2-dependent (ERBB2-EERBB2-D). We next validated the biological profiles of the subtypes and explored their prognostic/predictive value in external cohorts, namely the NeoALTTO trial, SCAN-B, I-SPY2, METABRIC and TCGA. IM tumors presented better survival outcomes, in contrast to Mes/S and P/Met tumors, while LUM tumors and ERBB2-EERBB2-D were characterized by low and high rates of pathological complete response, respectively. Of note, these molecular subtypes provide the rationale for treatment approaches leveraging the heterogeneous biology of HER2-positive breast cancer.

  Study EGAS50000000525

• The taxonomic composition of the human microbiome of CRC patients and healthy donors

  In this study, single homogenized human stool samples from five healthy donors and one CRC patient were analyzed by shotgun sequencing to generate a baseline for the establishment of humanized microbiome to be transferred into mice. Briefly, fecal DNA was isolated from these samples using the DNeasy PowerSoil DNA Isolation Kit (Qiagen), according to the kit manual, except that samples were homogenized using the TissueLyzerII for five minutes at 40Hz instead of intense vortexing. The extracted DNA was further purified and concentrated using the DNA Clean and Concentrator Kit (Zymo Research). This DNA was thawed and refrozen once at -80°C for 16s sequencing before being thawed again for WGS. 200 ng of DNA was used for metagenomic library preparation using the Next Generation Sequencing DNA Library Prep Kit (Westburg), according to the user manual using an average insert size of libraries of 450 base pairs. Following quantification and quality control (Qubit® dsDNA BR Assay Kit (Thermofisher Scientific); Bioanalyzer High Sensitivity DNA Analysis Kit (Agilent Technologies)), sequencing was performed on a NextSeq500 instrument using 2×150 base pair read lengths.

  Study EGAS50000000759

• Transcriptomic hallmarks and RNA isoform diversity in human neurodegenerative disease

  Single-cell RNA-sequencing (scRNA-seq) has enabled high-throughput RNA analysis at an unprecedented resolution. One staple analysis in scRNA-seq is identifying differentially expressed genes (DEGs) between two populations. However, the mRNA isoforms have often been overlooked in these transcriptomic analyses due to technical limitations. In order to investigate the mRNA isoforms at the single-cell level, we used a targeted single-nucleus long-read RNA-sequencing method combining single-nucleus short-read RNA-seq (snRNA-seq), target gene enrichment, and Iso-Seq. We performed snRNA-seq and targeted single-nucleus Iso-Seq (snIso-Seq) on 25 human frontal cortices from individuals with AD, DLB, or PD, and age-matched controls. Our short-read transcriptome analysis revealed that DLB samples had more overlapping DEGs with AD than PD despite sharing pathological hallmarks with PD. We observed vast RNA isoform diversity in the 50 genes targeted, including genes that were not differentially expressed in our short-read data. We also identified several novel isoforms that are not included in the reference annotation. Our data collectively characterize the transcriptome of various neurodegenerative diseases and illuminate the mRNA isoform diversity.

  Study EGAS50000000132

• Genetic coupling of enhancer activity and connectivity in gene expression control

  In this study, we investigated whether the genetic effects of eQTLs on enhancer activity and enhancer-promoter contacts (connectivity) are inherently coupled or distinct. To address this, we performed Capture Hi-C targeting ~1100 eQTL loci, along with ATAC-seq and RNA-seq, in monocytes from 34 healthy male donors. Using single-modality and multi-modality Bayesian statistical approaches, we demonstrated that eQTL effects are predominantly shared across enhancer activity, connectivity, and gene expression. However, we also identified a locus where an eQTL disrupts a CTCF binding site at a putative insulator, facilitating contacts between the PCK2 gene promoter and downstream enhancers. To validate these findings, we used 4C-seq to confirm the allelic effects of eQTLs on enhancer-promoter connectivity at the PCK2 locus and another eQTL region associated with THBS1 expression, both of which showed consistent eQTL effects across all three modalities. Additionally, we conducted Hi-C on monocytes from two donors to map reference TADs and performed ChIP-seq to confirm the allelic effects of the PCK2 eQTL on CTCF binding. Together, these data provide comprehensive insights into the genetic coupling of enhancer activity, connectivity, and gene expression.

  Study EGAS50000000756

• Atezolizumab Plus Chemotherapy With or Without Bevacizumab in Advanced Biliary Tract Cancer: Clinical and Biomarker Data From the Randomized Phase II IMbrave151 Trial

  PURPOSE Biliary tract cancers (BTCs) harbor an immunosuppressed tumor microenvironment (TME) and respond poorly to PD-1/PD-L1 inhibitors. Bevacizumab (anti-VEGF) plus chemotherapy can promote anticancer immunity, augmenting response to PD-L1 inhibition. PATIENTS AND METHODS This randomized, double-blind, proof-of-concept Phase II study enrolled patients (n=162) with previously untreated advanced BTC (IMbrave151; NCT04677504). Patients were randomized 1:1 to receive 3-weekly cycles of atezolizumab (1200 mg) plus bevacizumab (15 mg/kg) or atezolizumab plus placebo until disease progression or unacceptable toxicity. All patients received cisplatin (25 mg/m2) plus gemcitabine (1000 mg/m2; CisGem) on Days 1 and 8 every 3 weeks for up to 8 cycles. Stratification of patients was by disease status, geographic region, and primary tumor location. The primary endpoint was progression-free survival (PFS). No formal hypothesis testing was performed. Exploratory correlative biomarker analysis was undertaken using transcriptome analysis (n=95) and mutation profiling (n=102) on baseline tumor samples. RESULTS Between February and September 2021, 162 patients were enrolled. Median PFS was 8.3 months in the bevacizumab arm and 7.9 months in the placebo

  Study EGAS50000000387

• Distinct immune cell infiltration patterns in PDAC exhibit divergent immune cell selection and immunosuppressive mechanisms

  Pancreatic ductal adenocarcinoma has a dismal prognosis. A comprehensive analysis of single-cell multi-omic data from matched tumour-infiltrated CD45+ cells and peripheral blood in 12 patients, and two published datasets, reveals a complex immune infiltrate. Patients have either a myeloid-enriched or adaptive-enriched tumour microenvironment. Lymphocyte enrichment is intrinsically linked with highly distinct B and T cell clonal selection, diversification, and differentiation. Using TCR data, we see the largest clonal expansions in CD8 effector memory, senescent cells, and highly activated regulatory T cells which are induced within the tumour from naïve cells. We identify pathways that potentially lead to a suppressive microenvironement, including targets of immunotherapy and checkpoints including TIGIT/PVR and SIRPA/CD47. PDAC patient tumors from the APACT clinical trial revealed that patients with myeloid enrichment had the shorter overall survival compared to those adaptive cell enrichment. We identify multiple therapeutic targets that should form the basis for rational design including boosting of B cell responses, targeting immunosuppressive macrophages, and specific Treg depletion approaches. The dataset is scRNA-seq data from PDAC patients.

  Study EGAS50000000726

• Rapid brain tumor classification from sparse epigenomic data

  Although the intraoperative, molecular differential diagnosis of the approximately one hundred different brain tumor entities described to date has been a goal of neuropathology in the last decade, this has not yet been achieved in a clinically relevant time frame of less than one hour after biopsy collection. Recent advances in third-generation sequencing technologies have brought this once-elusive goal within reach. However, established machine learning techniques rely on concepts and methods, impractical for live diagnostic workflows in clinical applications. Here, we present MethyLYZR, a Naïve Bayesian framework enabling fully tractable live classification of cancer epigenomes. MethyLYZR can be run in parallel with an ongoing Nanopore experiment with negligible computational cost and provides clinically relevant and accurate cancer classification results within 15 minutes of sequencing. Therefore, only the time required for DNA extraction and the Nanopore sequencer's maximum parallel throughput remain limiting factors for even faster time-to-results. We demonstrate the potential utility of the MethyLYZR framework not only for the neurosurgical intraoperative use case but also for other oncologic indications and cell-free DNA from liquid biopsies.

  Study EGAS50000000559

• The genomic basis of childhood T-lineage acute lymphoblastic leukemia

  T-lineage acute lymphoblastic leukemia (T-ALL) is a high-risk tumor that has eluded comprehensive genomic characterization, in part due to the high frequency of non-coding genomic alterations resulting in oncogenic dysregulation. Here we report integrated genome and transcriptome sequencing of T-ALL tumor and remission samples obtained from over 1300 uniformly treated children with T-ALL, coupled with epigenomic and single cell analysis of malignant and normal T cell precursors. Integrated analysis identified 15 subtypes with distinct genomic drivers, gene expression, developmental state and outcome. Integration with chromatin topology analyses enabled elucidation of multiple mechanisms of enhancer deregulation that involve enhancers and genes in a subtype-specific fashion. We show that the immunophenotypically-described, high-risk entity of early T-cell precursor ALL is superseded by a broader category of leukemias of variable immunophenotype with diverse genomic alterations of a core set of genes encoding regulators of hematopoietic stem cell development. Numerous genetic alterations and disease subtypes emerged as independent predictors of survival and treatment failure in univariable and multivariable outcome models. These findings provide a roadmap for the classification, risk stratification and mechanistic understanding of this disease.

  Study EGAS50000000016

• Single cell atlas of large B-cell lymphoma

  Large B-cell lymphomas (LBCL) are clinically and biologically heterogeneous lymphoid malignancies, shaped by complex tumor and microenvironments that play a critical role in disease progression. There has been no direct and comprehensive single-cell characterization of tumor microenvironment in LBCL that simultaneously captures lymphoid, myeloid, and stromal populations across a sufficiently large cohort to reveal reproducible and biologically meaningful patterns. Here, we performed single-nucleus multiome profiling, coupled with bulk DNA and RNA sequencing, on 232 tumor and control biopsies to comprehensively map cell subpopulations of LBCL. Through this integrative approach, we identified Lymphoma Microenvironment Archetype Profiles (LymphoMAPs)—distinct, reproducible patterns defined by the co-occurrence of functionally specialized cell subsets. We further characterized the cell-cell communication pathways that shape these transcriptional phenotypes within archetypes. Importantly, we discovered associations between specific LymphoMAPs and clinical and molecular features of LBCL, including response to CD19 CAR T-cell therapy. Our study provides a valuable resource of lymphoid, myeloid, and non-hematopoietic cell states in LBCL and offers a framework for understanding how the abundance and functional state of cells within the tumor microenvironment may drive inter-patient variability in therapeutic outcome.

  Study EGAS50000001022

• Cambridge Neoadjuvant Cancer of the Prostate (CANCAP03): A Window Study into the Effects of Olaparib ± Degarelix in Primary Prostate Cancer.

  Study: Cambridge Neoadjuvant Cancer of the Prostate (CANCAP03): A Window Study into the Effects of Olaparib ± Degarelix in Primary Prostate Cancer. EudraCT Number: 2014-004417-86 Clinicaltrials.gov Number: NCT02324998 Purpose: Investigate combined PARP and androgen inhibition in primary prostate cancer. Aim: To understand the biological mechanism(s) underlying clinical efficacy, especially in the absence of mutations in homologous recombination (HR) repair pathways. Dataset: Single-nuclear RNA sequencing data from post-treatment radical prostatectomy tissue cores for six patients, 3 from each of olaparib and olaparib and degarelix cohorts. Samples were processed for nuclear isolation and extraction, prior to single-nuclear RNA sequencing. 10x libraries were constructed using the 10x supplied protocol, briefly: individual nuclei were trapped into droplets with a single bead, and poly A transcripts were barcoded during reverse transcription. Complementary DNA was amplified, and libraries were prepared for sequencing on NovaSeq 6000 using S4 flow cell (100k reads per cell). FASTQ files were processed using CellRanger pipeline (10X Genomics). Outputs from the CellRanger pipeline were analysed using the Seurat package (v4.9.9.9049) in R (v4.2.0). Processed data is provided as a Seurat object.

  Study EGAS50000000880

• Biallelic loss of function PAX4 variants are a Novel Cause of Transient Neonatal Diabetes

  The identification of genes in which loss of function (LoF) variants cause neonatal diabetes mellitus (NDM) can provide unique insights into the development and function of human pancreatic beta cells. We describe the identification of biallelic PAX4 LoF variants in 2 unrelated individuals with transient NDM: one due to a homozygous p.(Arg126*) PAX4 stop gain variant and a second due to a homozygous c.-352_104del 2.65kb deletion affecting the first 4 exons and promoter of the PAX4 Matched Annotation from NCBI and EMBL-EBI select transcript. We confirmed reduced transcript expression due to nonsense mediated decay for the e p.(Arg126*) variant in CRISPR-edited iPSC-derived pancreatic endoderm cells. Analysis of PAX4 binding through integration of CUT&RUN and RNA-sequencing in PAX4 depleted endocrine cells identified genes directly regulated by PAX4 involved in both pancreatic islet development and glucose-sensitive insulin secretion. The remitting and relapsing diabetes in both probands demonstrates that in contrast to mice PAX4 is not essential for the development of pancreatic beta-cells and supports differences in the requirements of functional beta-cell mass across the life-course.

  Study EGAS50000000857

• Ultra-Low-Pass Whole Genome Sequencing of Cell-Free DNA for Large B-Cell Lymphoma

  While deep targeted DNA sequencing of liquid biopsies has shown prognostic utility in large B-cell lymphoma (LBCL), the routine clinical adoption of these assays remains limited due to their high costs. Here, leveraging a well-annotated cohort of encompassing both frontline and relapsed/refractory (R/R) LBCL, we profiled patient plasma samples with two complementary modalities – ultra-low-pass whole genome sequencing (ulpWGS) and deep targeted DNA sequencing, the former being a cost-effective method to profile large scale chromosomal abnormalities and estimate tumor burden. Our findings revealed a strong association of high cell-free tumor burden by both genomic profiling modalities with established measures of tumor burden and patient survival. Notably, the associations with survival remained statistically significant after accounting for international prognostic index scoring. Furthermore, we showed that del(17p) in circulating tumor DNA as detected by ulpWGS was strongly associated with TP53 mutation status and predicted for significantly inferior patient outcome in frontline but not R/R LBCL patients. Our study demonstrates that ulpWGS can provide robust prognostic biomarkers for both frontline and R/R LBCL, highlighting its broad applicability for risk stratification.

  Study EGAS50000001027

• Final results and ctDNA analyses from the randomized phase 3 IMpassion031 trial evaluating peri-operative atezolizumab for early-stage triple-negative breast cancer

  In IMpassion031, adding atezolizumab to neoadjuvant chemotherapy significantly improved pathologic complete response (pCR) rate (primary endpoint) for stage II/III triple-negative breast cancer (TNBC). Patients with stage II/III TNBC were randomized 1:1 to neoadjuvant chemotherapy plus peri-operative atezolizumab or pre-operative placebo. Descriptive secondary endpoints included event-free, disease-free, and overall survival. Long-term outcomes favored the atezolizumab group (event-free survival HR 0.76, 95% CI: 0.47–1.21). Among patients without pCR, 14/70 (20%) atezolizumab-treated and 33/99 (33%) placebo-treated patients received additional adjuvant therapy, frequently capecitabine. In exploratory analyses, pCR was prognostic for long-term outcome. Patients with baseline ctDNA-negative status (6%) had excellent long-term outcomes. Most patients (87%) had cleared ctDNA at surgery. ctDNA-positive status at surgery identified a subset of non-pCR patients with poorest prognosis. Long-term safety was consistent with primary results. Adding atezolizumab to chemotherapy for stage II/III TNBC is associated with more favorable long-term outcomes. ctDNA dynamics provide prognostic value beyond pCR. ClinicalTrials.gov identifier: NCT03197935

  Study EGAS50000000974

• DynaTag for efficient profiling of transcription factors in small samples and single cells

  Systematic discovery of transcription factor (TF) landscapes in low-input samples and at single cell level is a major challenge in the fields of molecular biology, genetics, and epigenetics. To address this, we present here cleavage under Dynamic targets and Tagmentation (DynaTag), which allows robust mapping of specific TF-DNA interactions via the use of a physiological salt solution during sample preparation. Using this method, we revealed occupancy alterations for 15 different TFs in stem cell and cancer tissue models. Of note, we showcased for NANOG, MYC, and OCT4 differentially occupied genomic loci during stem cell differentiation, both in bulk and at single-cell resolution. DynaTag demonstrated superior signal-to-background ratio and occupancy resolution in comparison to CUT&RUN and ChIP-seq. In addition, a model for small cell lung cancer chemotherapy treatment in vivo showed increased chromatin occupancies of the gain-of-function mutant P53 p.R248Q, FOXA1, and MYC at altered gene pathways, such as epithelial-mesenchymal-transition. Altogether, these results highlight the technical advancement of DynaTag to uncover TF landscapes in cell models and in tissues of a complex disease model.

  Study EGAS50000001074

• TDP43 proteinopathy leads to divergent cryptic splicing in the cortex and spinal cord

  Mislocalization of the nuclear TAR DNA-binding protein 43 (TDP43) is a hallmark of ALS and FTD which leads to de-repression and inclusion of cryptic exons, promising biomarkers of TDP43 pathology in a spectrum of neurodegenerative diseases. However, most cryptic exons to date have been identified from in vitro models or a single cortical FTD dataset, and little is known about cryptic splicing in the spinal cord, or within different neuronal subtypes. We meta-analyzed published bulk RNAseq datasets representing 1,778 RNAseq profiles of ALS and FTD post-mortem tissue, and in vitro models with experimentally depleted TDP43. We identified 142 cryptic splices, including 76 novel events. We found a novel pattern of spinal cord cryptic splicing, validated in an independent cohort by qPCR, which differed significantly from cortical and in vitro splicing. Finally, leveraging multiple public single-nucleus RNAseq datasets of ALS and FTD motor and frontal cortex, we confirmed the elevation of cortical-enriched splices in disease and localized them to layer-specific neuronal populations. This catalog of cryptic splices could inform efforts to develop biomarkers for tissue-specific and cell type-specific TDP43 pathology.

  Study EGAS50000000575

• Vaccine-Expanded Plasmablast-like B Cells Are Associated with Response to Dendritic Cell Therapy in Metastatic Melanoma

  Dendritic cells (DCs), either monocyte-derived (moDCs) or naturally circulating, have been employed in cancer treatment across both metastatic and adjuvant settings. DC vaccines (DC Vax) are capable of inducing tumor-specific immune responses, yet clinical outcomes have thus far fallen short of expectations. Here, we performed a multi-omics analysis of the apheresis product, including CD14⁺ monocytes, and the final DCVax formulation. Among all the components analyzed, the most striking finding was the unexpected presence of CD19⁺ B cells within the DC Vax formulation, which distinguished Responders from Non-Responders to moDC vaccination in metastatic melanoma. The abundance of B cells in DC Vax correlated with their peripheral frequency; however, within the vaccine, these cells acquired a CD20lowCD27hiCD38hi plasmablast-like phenotype, distinct from their autologous circulating B-cell precursors. Notably, a parallel observation was made in Responders, whose pre-treatment lesions were enriched in B cells organized as Tertiary Lymphoid Structures. These results underscore the importance of B cells in both baseline immunity and DCVax formulation, hinting at unexplored mechanisms by which they may influence vaccine effectiveness.

  Study EGAS50000001177

• M116 Proteome Extracellular Vesicle Profiling

  For extracellular vesicles (ECV) extraction from plasma, 1M ammonium acetate was added to precipitate ECVs on ice for 45 min. Then, 100 mM ammonium acetate was added to the mixture, and ECVs were precipitated by centrifugation at 20,000g for 30 min. ECVs were washed with 50 mM ammonium bicarbonate (Sigma-Aldrich, Cat. 1066-33-7). Then, 600 µl of 1% ammonium deoxycholate (Sigma Aldrich, Cat. K2755-1MG) were added. The concentration of protein in each sample was measured using a bicinchoninic acid assay (BCA assay). Ammonium bicarbonate was used to dilute 500 µg of protein into a final volume of 500 µl. Next, dithiothreitol (Sigma-Aldrich, Cat. D9779) was added to obtain a final concentration of 20 mM, followed by iodoacetamide (Sigma-Aldrich, Cat. I6125) to a final concentration of 40 mM. Next, trypsin (Roche, Cat. RTRYP-RO) was added to the sample in a 1:25 protein ratio and incubated at 37°C overnight. Next day, formic acid (ThermoFisher Scientific, Cat. 28905) in a final concentration of 0.1% was added and extraction of proteins was done using Empore™ Solid Phase Extraction Cartridges (3M), following manufacturer’s instructions. The eluted samples were then centrifuged for 90 min. using a speed vacuum centrifuge (Thermo, RC1010), followed by snap freezing in liquid nitrogen. Then, the samples were kept in a freeze dryer (LyoDry Compact Benchtop, MechaTech) overnight. Next, the samples were reconstituted in 30 µl of 0.1% formic acid (FA) and an o-Phthaladehyde (Oparil) assay was performed to determine the concentration of each sample. After that, the sample was prepared in a concentration of 0.5 µg/µl using 0.1% FA and alcohol dehydrogenase (ADH). The samples were prepared in glass mass spectrometry vials for proteomic analysis using a Waters Synapt G2Si High-Definition Mass Spectrometry (Waters Corporation) operated by the MassLynx 4.1., 110 min. running time with 2 µl of an injection containing 1 µg of peptide. Quality controls were also run along with samples to guarantee consistency. Pooled quality controls were made from all samples, in which the samples were run at the beginning, middle and end of the mass spectrometry run. Samples were randomized before running the experiment. The proteomic data was then imported into Progenesis software 4.2 (Nonlinear Dynamic, UK) to identify and quantify peptides and proteins.

  Dataset EGAD50000001679

• Response to tumor-infiltrating lymphocyte adoptive therapy is associated with preexisting CD8+ T-myeloid cell networks in melanoma

  Adoptive cell therapy (ACT) using ex vivo–expanded tumor-infiltrating lymphocytes (TILs) can eliminate or shrink metastatic melanoma, but its long-term efficacy remains limited to a fraction of patients. Using longitudinal samples from 13 patients with metastatic melanoma treated with TIL-ACT in a phase 1 clinical study, we interrogated cellular states within the tumor microenvironment (TME) and their interactions. We performed bulk and single-cell RNA sequencing, whole-exome sequencing, and spatial proteomic analyses in pre- and post-ACT tumor tissues, finding that ACT responders exhibited higher basal tumor cell–intrinsic immunogenicity and mutational burden. Compared with nonresponders, CD8+ TILs exhibited increased cytotoxicity, exhaustion, and costimulation, whereas myeloid cells had increased type I interferon signaling in responders. Cell-cell interaction prediction analyses corroborated by spatial neighborhood analyses revealed that responders had rich baseline intratumoral and stromal tumor–reactive T cell networks with activated myeloid populations. Successful TIL-ACT therapy further reprogrammed the myeloid compartment and increased TIL-myeloid networks. Our systematic target discovery study identifies potential T-myeloid cell network–based biomarkers that could improve patient selection and guide the design of ACT clinical trials.

  Study EGAS50000001217

• Non-neuronal TGF-β–mediated extracellular matrix remodeling drives neurodegeneration in a PSP-Richardson syndrome model

  Progressive supranuclear palsy-Richardson syndrome (PSP-RS) is a devastating neurodegenerative disorder marked by rapid motor decline, cognitive impairment, and widespread accumulation of pathological tau. Yet, the molecular triggers of sporadic PSP-RS remain elusive. In this study, we combined patient-derived midbrain organoids with integrative multi-omics to uncover a previously unrecognized pathological cascade. At the core lies an early and aberrant activation of TGF-β signaling, which initiates profound extracellular matrix (ECM) remodeling and aberrant integrin-mediated focal adhesion. This cascade leads to persistent hyperactivation of PI3K/AKT and MAPK/ERK pathways, disrupting cytoskeletal integrity, impairing autophagy, promoting pathological tau phosphorylation and mislocalization, and ultimately driving synaptic degeneration. Single-cell RNA sequencing revealed an ectopic expansion of vascular leptomeningeal-like cells (VLMCs), a collagen-producing population, accompanied by a broad ECM gene dysregulation unique to PSP-RS. Remarkably, pharmacological blockade of TGF-β, AKT, ERK, or mTOR signaling reversed tau pathology and restored synaptic integrity. These findings position TGF-β-driven ECM remodeling as a master regulator of neurodegeneration in sporadic PSP-RS and unveil multiple actionable therapeutic targets for this currently untreatable disease.

  Study EGAS50000001236

• Pipeline Olympics: Continuous benchmarking of computational workflows for DNA methylation sequencing data

  DNA methylation is a widely studied epigenetic mark and a powerful biomarker of cell type, age, environmental exposure, and disease. Genome-wide DNA methylation profiling is usually performed using whole-genome sequencing following selective conversion of unmethylated cytosines into thymines via bisulfite treatment or enzymatic methods. Numerous software tools facilitate processing of raw sequencing data for DNA methylation profiling, but a comprehensive benchmarking has been lacking. In this study, we systematically compared methods and tools for processing DNA methylation sequencing data. We established a dedicated benchmarking dataset comprising genome-wide DNA methylation profiling data for four reference samples using five experimental protocols each. The results of data processing were compared to highly quantitative locus-specific DNA methylation measurements available for 46 genomic loci. Based on this gold-standard reference dataset, we identified workflows that consistently demonstrated high performance. We present our results in an interactive application that provides a “living benchmark” ready for continuous updating as new workflows are proposed and require benchmarking against existing methods. In summary, our study provides guidance for any laboratory studying DNA methylation on a genome-wide scale.

  Study EGAS50000000541

• Error-corrected flow-based sequencing at whole genome scale and its application to circulating cell-free DNA profiling

  Differentiating sequencing errors from true variants is a central genomics challenge calling for error suppression strategies that balance costs when sparse material or ultra-rare variants require high-accuracy methods. For example, circulating cell-free DNA (ccfDNA) sequencing for cancer monitoring is limited by sparsity of circulating tumor DNA (ctDNA), abundance of genomic material in samples, library preparation error rates and sequencing errors. Whole-genome sequencing (WGS) can overcome the low abundance of ccfDNA by integrating signal across the mutation landscape, but higher cost limits adoption compared with targeted panels, where ccfDNA abundance limits maximal coverage depths. Here we applied deep (~120x) lower-cost WGS (Ultima Genomics) with analytic error correction for tumor-informed ctDNA detection within the part-per-million range. We further developed high-coverage duplex error-corrected WGS of ccfDNA, achieving 7.7x10-8 error rates, to assess disease burden in melanoma and urothelial cancer patients without matched tumor sequencing. Together, we show that deeply sequenced error-corrected WGS accurately calls somatic variants, demonstrating the feasibility of WGS for tumor-informed and tumor-agnostic ctDNA detection.

  Study EGAS50000000844

• Whole exome sequencing of bladder tumors

  Bladder carcinomas are immunogenic, and patients with bladder cancer benefit from immune checkpoint therapy. This is correlated to a high tumor mutation burden, which provides a higher number of neoepitopes that can be recognized by tumor-specific CD8+ T cells. Intravesical Bacillus Calmette-Guérin (BCG) is used to treat non-muscle invasive bladder cancer (NMIBC), but its mechanism of action remains elusive. Most lymphocytes appearing in the urine of BCG-treated patients are CD4+ T cells though preclinical studies showed that CD8+ T cells are also necessary for BCG treatment efficacy. It is currently unknown which proportion of patients with non-metastatic bladder cancer develop a spontaneous antitumor CD8+ response, and if BCG treatment influences this response. were prepared using the xGen Dual Index UMI Adapters (IDT), and genomic DNA (gDNA) extracted (DNeasy kit, Qiagen) from formalin-fixed paraffin- embedded (FFPE) or snap-frozen tumor tissue, as well as autologous lymphocytes. Libraries were subjected to paired-end sequencing using Illumina HiSeq4000 with a coverage depth of 250X. Library preparation was carried out with the AllType FASTplex NGS 11 Loci Flex Kit (One Lambda), and sequencing was conducted on the Illumina MiSeq platform.

  Study EGAS50000001248

• MATISSE WES and bulk RNA-sequencing data

  High-risk cutaneous squamous cell carcinoma is frequently seen in frail, elderly patients and often requires mutilating surgical treatment and adjuvant radiotherapy (RT). The MATISSE trial (NCT04620200) was conducted to improve clinical prospects in CSCC patients eligible for extensive curative surgery, and evaluated early treatment response and immune dynamics following ultra-short neoadjuvant immunotherapy prior to standard of care surgery with or without radiotherapy. Fifty patients were randomized to receive two courses of nivolumab 3 mg/kg (NIVO, weeks 0&2) or nivolumab 3 mg/kg (weeks 0&2) and one course of ipilimumab 1 mg/kg (week 0 only, NIVO/IPI before surgery at week 4, with the primary objective being histopathological response based on residual viable tumor cells in the surgical specimen. Ten patients choose not to undergo surgery, of whom nine reached a clinical complete response during follow up, and were disease-free at median 34 months follow up. Serial tumor and blood samples at baseline, week 1, week 2 and week 4, enabled to investigate dynamical changes through bulk WES and RNA sequencing to potentially aid future de-escalation trials.

  Study EGAS50000001003

• Non-small cell lung cancer molecular subtypes and vulnerability to immunotherapy treatment combinations

  The phase 3 IMpower150 trial in treatment-naïve patients with metastatic non-small-cell lung cancer (NSCLC) demonstrated significantly longer progression-free (PFS) and overall survival (OS) with first-line atezolizumab (anti-PD-L1)-bevacizumab (anti-VEGF)-carboplatin-paclitaxel (ABCP) than with bevacizumab-carboplatin-paclitaxel (BCP). We characterised four molecular NSCLC subtypes identified by unsupervised clustering of transcriptomes of 564 pre-treatment primary tumour samples from IMpower150 using non-negative matrix factorization (NMF1-4). Each subtype had distinct tumour PD-L1 expression levels, epithelial characteristics, immune composition, and treatment outcomes. Both NMF2 (enriched in tumour proliferation signal, macrophages, and monocytes) and NMF4 (enriched in B-cells and T cells) had elevated tumour PD-L1 expression. Of these two, only NMF4 demonstrated PFS and OS benefits with ABCP versus either BCP or atezolizumab-carboplatin-paclitaxel (ACP). Patients with NMF1 (enriched in basal and squamous-like cells) had improved outcomes on ABCP compared with ACP or BCP; those with NMF3 (enriched in adenocarcinoma signatures) showed similar outcomes among treatments. These insights could help inform individualised first-line treatment for metastatic NSCLC.

  Study EGAS50000001272

• Dynamics of checkpoint receptors in γδ T cell subsets are associated with clinical responses during anti-PD-1 immunotherapies

  Gamma delta (γδ) T cells are innate-like lymphocytes with potent anti-tumor properties. Herein, we show that immune checkpoint receptors (ICRs) display differential expression and regulation by the JAK-STAT pathway in Vδ1 and Vδ2 cells and identify constitutive (e.g. TIGIT, PD-1) and inducible (e.g. TIM-3, LAG-3, CTLA-4) ICRs. In melanoma, all γδ T cell subsets downregulated AP-1 transcription factors, but Vδ1 cells specifically expressed high levels of ICR, TOX and inhibitory killer Ig-like receptor (KIR) transcripts, reminiscent of exhaustion. However, patient-derived cells were functionally competent, although induction of LAG-3 and CTLA-4 was impaired. During anti-PD-1 monotherapy, Vδ1 cells specifically bound high levels of therapeutic antibody but only in patients who responded to treatment, revealing a potential new prognostic marker for evaluating the efficacy of IC blockade (ICB) therapy. Finally, expression of KIR genes in Vδ1 cells was downregulated in response to successful ICB therapy. Collectively, our data indicate an intricate relationship between ICRs and γδ T cells and reveal novel approaches by which these cells can be harnessed in order to discern or improve cancer immunotherapy.

  Study EGAS50000001270

• FLT3-targeting and places FLT3 as an ideal therapeutic target to selectively eradicate LSCs, while sparing HSC.

  CRISPR/Cas9 FLT3 knock-out (FLT3-KO) in human HSCs and LSCs followed by functional xenograft assays to test their ability to regenerate human hematopoiesis and leukemia, respectively. FLT3-KO in HSCs from human fetal liver, cord blood and adult bone marrow showed no impairment in multilineage hematopoiesis in primary and secondary xenografts. By contrast, FLT3-KO LSCs from 6 of 7 FLT3-ITD mutated AMLs were able to generate short-term engraftment but were completely exhausted by 12 weeks. Thus, FLT3 is essential for LSC long-term propagation. This dependency was unique to FLT3-ITD AML samples as non-FLT3-ITD AML samples generated leukemic grafts upon FLT3-KO. Transcriptomic analysis revealed that FLT3-KO induced downregulation of DNA repair and cell cycle checkpoints, uniquely in FLT3-ITD AML, but not in healthy HSCs or other AMLs. Our research highlights a critical distinction between healthy HSCs and LSCs: whereas healthy hematopoiesis proceeds unperturbed upon FLT3-KO, FLT3-ITD leukemogenesis is impaired through elimination of LSCs. This evidence underscores the necessity for more potent FLT3-targeting and places FLT3 as an ideal therapeutic target to selectively eradicate LSCs, while sparing HSC.

  Study EGAS50000000437

• Proteomic predictors of individualized nutrient-specific insulin secretion in health and disease

  Population level variation and molecular mechanisms behind insulin secretion in response to carbohydrate, protein, and fat remain uncharacterized despite ramifications for personalized nutrition. We now define prototypical insulin secretion dynamics in response to the three macronutrients in islets from 140 cadaveric donors, including those diagnosed with type 2 diabetes. We leverage the insulin response heterogeneity and use transcriptomics and proteomics to identify molecular pathways of specific nutrient responsiveness. Surprisingly, we find robust insulin secretion to fatty acid stimulus in ~8% of donors, challenging the idea that fat has negligible effects on insulin release. Distinct islet proteomes with differences in metabolic signalling networks convey this hyper-responsiveness to fat relative to carbohydrate. By comparing human islets to human embryonic stem cell-derived islet clusters, we show that, unlike glucose-responsiveness, fat hyper-responsiveness is equivalent and may be a hallmark of functionally immature cells. Our study represents the first comparison of dynamic responses to nutrients and multi-omics analysis in human insulin secreting cells. Responses of different people’s islets to carbohydrate, protein, and fat lay the groundwork for personalized nutrition.

  Study EGAS00001007241

• Whole-genome methylation profiling of menstrual stem cells identifies novel biomarkers for endometriosis

  Endometriosis is a highly prevalent condition that remains underdiagnosed and poorly managed due to the lack of clinically validated biomarkers and limited understanding of its underlying pathophysiology. Although menstrual blood–derived stem cells (MenSCs) have been implicated in disease pathogenesis, their potential as a diagnostic source has not been systematically explored. We conducted a clinical study including 42 women (19 with endometriosis and 23 controls) to evaluate whether DNA methylation profiles from freshly isolated MenSCs could identify disease-specific biomarkers. Whole-genome DNA methylation sequencing identified differentially methylated regions (DMRs) enriched in genes associated with key features of endometriosis, including inflammation, tissue remodelling, and developmental processes. These DMRs robustly discriminated cases from controls, independently of technical and clinical confounders. Machine learning models trained and validated on these regions demonstrated strong diagnostic performance (specificity 83%, sensitivity 79%). Integration with an independent single-cell RNA sequencing dataset suggested that the identified DMRs may influence gene expression, supporting their biological relevance. Collectively, these findings highlight MenSC DNA methylation profiling as a promising, non-invasive strategy for early diagnosis and personalised management of endometriosis. (PubMed IDs will be provided once the manuscript is published.)

  Study EGAS50000001640

• Non-canonical NF-κB signaling skews B cells away from germinal center to low-affinity effector fate

  The non-canonical (NC) NF-κB pathway is a central regulator of B-cell development and survival, yet its role during B-cell mediated adaptive responses remains uncharted. Here, we uncover an unanticipated regulatory function for NC NF-κB in shaping humoral immunity. High-resolution mapping of NC NF-kB status in B cells from human lymphoid tissues revealed low activity in germinal center B-cells (GCB) than effector subsets, suggesting a selective attenuation. Consistently, mice model harboring B-cells with enforced NC NF-κB activation revealed near- complete loss of GCB and high-affinity antibodies, following vaccination. GCB loss was not driven by intrinsic defects in formation or expansion, but by an immune-niche dysregulation marked by regulatory T-cell expansion. Mechanistically, NC NF-κB activation promoted B-cell- derived IL-10, and blocking this axis normalize the niche and rescue GCB. Concurrently, NC NF- κB signaling boosted low-affinity antibody responses by transcriptionally priming follicular B- cells for plasmacytic differentiation. Together, these findings position NC NF-κB as a tunable rheostat linking B-cell signaling to immune-niche control and humoral immunity.

  Study EGAS50000001646

• Adaptive Therapy Exploits Fitness Deficits in Chemotherapy-Resistant Ovarian Cancer to Achieve Long-Term Tumor Control Open Access

  Drug resistance results in poor outcomes for cancer patients. Adaptive therapy is a potential strategy to address drug resistance that exploits competitive interactions between sensitive and resistant subclones. Here, we showed that adapting carboplatin dose according to tumor response (adaptive therapy) significantly prolonged survival of murine ovarian cancer models compared to standard carboplatin dosing, without increasing mean daily drug dose or toxicity. Platinum resistant ovarian cancer cells exhibited diminished fitness when drug was absent in vitro and in vivo, which caused selective decline of resistant populations due to reduced proliferation and increased apoptosis. Conversely fitter, sensitive cells re-grew when drug was withdrawn. Using a bioinformatics pipeline that exploits copy number changes to quantify the emergence of treatment resistance, analysis of cell‐free DNA obtained longitudinally from ovarian cancer patients during treatment showed subclonal selection through therapy, and measurements of resistant population growth correlated strongly with disease burden. These preclinical findings pave the way for future clinical testing of personalized adaptive therapy regimens tailored to the evolution of carboplatin resistance in individual ovarian cancer patients.

  Study EGAS50000001142

• Whole Genome Sequencing of Neuroblastoma

  Neuroblastoma is a childhood tumour of the peripheral sympathetic nervous system. The pathogenesis has for a long time been quite enigmatic, as only very few gene defects were identified in this often lethal tumour1. Frequently detected gene alterations are limited to MYCN amplification (20%) and ALK activations (7%)2-5. Here we present a whole genome sequence analysis of 87 neuroblastoma of all stages. Only few recurrent amino-acid changing mutations were found. In contrast, analysis of structural defects identified in 18% of high stage neuroblastoma a local shredding of chromosomes, known as chromothripsis6. These tumours are associated with a poor outcome. Structural alterations recurrently affected ODZ3, PTPRD and CSMD1, which are involved in neuronal growth cone stabilization7-9. In addition, ATRX, TIAM1 and a series of regulators of the Rac Rho pathway were mutated, further implicating defects in neuritogenesis in neuroblastoma. Most tumors with defects in these genes were aggressive high stage neuroblastoma, but did not carry MYCN amplifications. The genomic landscape of neuroblastoma therefore reveals two novel molecular defects, chromotripsis and neuritogenesis gene alterations, which frequently occur in high risk tumors.

  Study EGAS00001000222

• Gut microbiome dynamics unravelled with metagenomics sequencing

  It has been already several years since we discovered that microbial therapeutics like Faecal Microbiota Transplantation can be used to treat Clostridium difficile infection (CDI) and Inflammatory bowel disease (IBD) with great success. The concept has challenged us to test the limits of our knowledge with the perspective to improve health in the most diverse conditions like Parkinson disease, depression, autism and obesity, just to mention a few. Note: the figure belongs to the original paper and is owned by the authors. Allyson L.Bird and colleagues from Genentech Inc. (San Francisco) recently published an interesting study building up on the key concept that “As microbial therapeutics are increasingly being tested in diverse patient populations, it is essential to understand the host and environmental factors influencing the microbiome.” The authors of this work analysed 1359 gut microbiome samples using shotgun metagenomics sequencing.  The results of the comparison done between microbiota composition across the different groups showed that biological sex is the strongest driver of diversity, while aging (age range studied: 20–69) appears to be the second factor. As smartly outlined in the authors’ own graphical abstract (figure to the left), Bacteroidota species consistently increased with donors' age while Actinobacteriota species, including Bifidobacterium, decreased. This study also integrated data from oncological patients to draw cancer-related evidences. After correcting the biases for technical and geographic differences, they found that cancer patients show significantly altered gut bacterial communities compared with their healthy counterparts. Eg. E.coli was found to be increased in cancer patients together with Bacteroidota/Firmicutes ratios. By comparing sequential samples taken from the same donors (413) within short-distance time points, the authors found that in the short term, intra-individual differences were less than the inter-individual ones, indicating a certain community stability being the norm. One of the pivotal points from this interesting study is that the data proceeds from extremely well-defined and selected donors. Specifically, the data for 949 selected healthy donors came from the Milieu Intérieur (MI) Consortium and extensive metadata, including demographic variables, serological measures, dietary information, and systemic immune profiles, are available. This precise stratification allowed for the robust detection of smaller changes, compared to former works. All sequencing data has been deposited in the EGA archive under the accession number EGAS00001004437. Enjoy this amazing open access paper published in the Experimental Journal of Medicine on January 2021.

  Blog gut-microbiome-dynamics-unravelled

• Whole-genome sequencing of two probands with hereditary spastic paraplegia reveals novel splice-donor region mutation and known pathogenic mutation in SPG11

  Hereditary spastic paraplegias (SPG) are a group of heterogeneous neurodegenerative disorders, which are often presented with overlapping phenotypes such as progressive paraparesis and spasticity. To assist the diagnosis of SPG subtypes, next-generation sequencing is often used to provide supporting evidence. In this study, we report the case of two probands from the same family with SPG symptoms, including bilateral lower limbs weakness, unsteady gait, cognitive decline, dysarthria and slurring of speech since age of 14. Subsequent whole-genome sequencing revealed that the patients are compound heterozygous for mutations in SPG11 gene, including the paternally inherited c.6856C>T (p.Arg2286*) mutation and the novel maternally inherited c.2316+5G>A splice donor region mutation. Mutations in SPG11 are the common cause of autosomal recessive Spastic Paraplegia type 11. According to the ClinVar database, there are already 101 reported pathogenic mutations in SPG11 that are associated with SPG. To our knowledge, this is the first report of SPG11 mutations in our local population. Novel splice mutation identified in this study would enrich the catalog of SPG11 mutations, potentially leading to better genetic diagnosis of SPG.

  Study EGAS00001001849

• Novel regional age-associated DNA methylation changes within human common disease-associated loci

  Background: Advancing age progressively impacts on risk and severity of chronic disease. It also modifies the epigenome, with changes in DNA methylation, due to both random drift and variation within specific functional loci. Results: In a discovery set of 2,238 peripheral-blood genome-wide DNA methylomes (MeDIP-seq) aged 19-82 years, we identified 71 age-associated Differentially Methylated Regions (a-DMRs, p < 1.85 x 10-8) within the Linkage Disequilibrium (LD) blocks of the NIH Catalogue of published GWAS SNPs. This included 52 novel regions, 29 within loci not covered by 450k or 27k Illumina array, and with marked enrichment for Poised Promoters and Enhancers across multiple cell types. In a replication set of 2,084 DNA methylomes, 95.7% of the a-DMRs showed the same direction of ageing effect, with 80.3% and 53.3% replicated to p < 0.05 and p < 1.85 x 10-8, respectively. Conclusion: By analysing the functionally enriched disease and trait-associated regions of the human genome, we identified novel epigenetic ageing changes, which could be useful biomarkers or provide mechanistic insights into age-related common diseases.

  Study EGAS00001001910

• Comprehensive molecular profiling identifies novel genetic drivers and subtypes underlying medulloblastoma

  Current therapies for medulloblastoma (MB), a highly malignant childhood brain tumor, impose debilitating effects on the developing child, warranting the development of molecularly targeted treatments with reduced toxicities. Prior studies have failed to disclose the full spectrum of driver genes and molecular processes operative in MB or adequately explain heterogeneity among molecular subgroups. Herein, we detail the somatic landscape across 491 sequenced MBs and molecular heterogeneity amongst 1,256 epigenetically analyzed cases, identifying subgroup-specific mutational signatures, driver genes, and pathway alterations including previously unappreciated actionable therapeutic targets. Integrative approaches for assigning driver genes to subgroups explained >65-70% of Group 3 and Group 4, doubling previous knowledge. Novel subtypes underlying Group 3 and Group 4 were differentially enriched for specific driver events, including hotspot in-frame indels targeting KBTBD4 and ‘enhancer hijacking’ driving activation of PRDM6. Thus, application of integrative genomics to an unprecedented cohort of clinical samples derived from a single childhood cancer entity has disclosed a series of new cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for treating MB.

  Study EGAS00001001953

• WTCCC2 Visceral Leishmaniasis (VL) samples

  A WTCCC2 project genome-wide association study for visceral leishmaniasis (VL) in individuals from India, Brazil and Sudan, genotyped on the custom Illumina 670k array. The WTCCC2 analysis of the Brazilian and Indian samples is described in Fakiola et al. [Nat Genet. 2013 Feb;45(2):208-13].It should be noted that due to expected family structure in the data, normal analyses of these data should include an estimation of the relatedness between the samples. For more details about sample collection for the project please refer to the Methods section of the paper above. The samples from India were all collected from Bihar state in northeastern India. The Brazilian samples were collected from a wide area of northern Brazil, and managed at two laboratories in Natal and Belem. For more details on the geographic distribution of the Brazil samples please see Jamieson et al. [Genes Immun. 2007 Jan;8(1):84-90]. We have not generated qc_passed or info files for the Sudanese data, a preliminary scan of the data indicated that there may be DNA quality issues with these samples.

  Study EGAS00001000773

• Global Anaplastic Thyroid Cancer Initiative

  Thyroid cancer is the most common endocrine malignancy. Most thyroid cancers are of the well-differentiated (non-aggressive) phenotype, and are almost always cured with standard treatments. In contrast, anaplastic thyroid cancer (ATC) is rare, accounting for only 1% of thyroid cancers, however it is perhaps the most lethal human malignancy with an average survival of 3 to 6 months. ATC presents with dramatic and rapid onset of airway and esophageal blockage with frequent spread to the lungs. Interestingly, 21–79% of ATCs have coexisting areas or a previous history of well-differentiated cancers. Although this suggests a progression from well-differentiated cancers to ATC, to date there is little molecular conformation that such a progression exists versus de novo generation of ATC. As ATC is extremely rare and most cases are inoperable, there are a paucity of tissue samples for study at any one center. The Global Anaplastic Thyroid Cancer Initiative (GATCI) aims to unite international institutions in order to pool samples (ATC with or without a paired well-differentiated component) for a comprehensive analysis of the genomic landscape of this disease.

  Study EGAS00001002234

• Field_effect_of_healthy_and_diseased_livers_WGS

  Recent work in the Campbell group has revealed somatic mutations present in normal, non-cancerous human skin. A subset of the mutations conferred selective advantages to the host cells, leading to clonal expansions and raising the risk for future cancer development. Capturing such somatic mutations in normal tissue is important to advance our understanding about carcinogenesis and could provide prospective medical insights. In this project, our goal is to detect somatic mutations in normal (pre-cancerous) liver tissue. Using Laser Microdissection technology, we will dissect individual liver lobules from patient samples and submit these to sequencing. For each patient sample, we aim to sequence multiple lobules to characterise the mutagenic burden. Samples will be taken from patients with different liver disease aetiologies, including alcoholism and obesity, with a view on distinguishing the prevalent mutation types occurring in each disease context. We will perform both whole genome and targeted sequencing, initially using the WTSI cancer panel. Later we aim to use a novel bait set that captures both cancer genes as well as genes relevant to the non-cancerous samples (ie. genes implicated in hereditary disorders, immune sequences).

  Study EGAS00001002413

• Cryptic Relatedness in the Singapore Living Biobank Project

  The Singapore Living Biobank is a collection of healthy population-based Chinese and Malay individuals, for the purpose of phenotype recall study of high-impact variant carriers. These individuals are sampled from two studies: Multi-Ethnic Cohort (MEC), and the Singapore Health 2012 (SH2012). The MEC is a population-based cohort initiated in 2007 to investigate the genetic and lifestyle factors that affect the risk of developing chronic diseases such as diabetes and cardiovascular outcomes in the three ethnic groups (Chinese, Malays, and Indians). The SH2012 study is a population-based cross-sectional survey conducted in Singapore between 2012 and 2013, with over-sampling of Malays and Indians. Participants in MEC and SH2012 completed a similar set of questionnaire components, health examination, and biochemisty panels. Description of the MEC and SH2012 studies can be found at http://blog.nus.edu.sg/sphs/. We generated whole-exome sequencing data and Illumina OmniExpress array genotyping data for 1,299 Chinese and 1,229 Malays from the Singapore Living Biobank. This study includes a subset of 762 individuals that were found to be closely related (≤3rd degree), including 263 Chinese and 499 Malays.

  Study EGAS00001002619

• Duplexseq_of_the_interstrand_crosslinks_WGS

  One of the most dangerous forms of DNA damage are interstrand crosslinks (ICLs), which covalently crosslink the two strands of the DNA double helix. The repair of these lesions is crucial for cellular survival due to their ability to block transcription and DNA replication. Initially, the major pathway that has been described in ICL repair involves a network of 22 genes that are mutated in a severe human genetic disease known as Fanconi Anemia (FA). Using synthetic lethality screens in the near-haploid human HAP1 cell line, we recently identified two potentially novel regulators of ICL repair, C1orf112 and THAP12. Loss of C1orf112 and THAP12 causes hypersensitivity to ICL-inducing DNA damaging agents, such as Mitomycin C (MMC). Additionally, C1orf112-depleted cells show elevated levels of micronuclei and accumulation of DNA damage in S-phase. To better understand how C1orf112 and THAP12 mediate the repair of ICLs, we want to perform mutational signature analysis, using the BotSeq method. Therefore, WT, C1orf112 and THAP12 knockout cells were cultured in vehicle or MMC treated conditions for 10 days and the genomic DNA was isolated. FANCA and FANCD2 knockout cells are taken along as controls in this experimental setting.

  Study EGAS00001006545

• SHANK2 mutations associated with Autism Spectrum Disorder cause hyperconnectivity of human neurons

  Heterozygous loss-of-function mutations in the synaptic scaffolding gene SHANK2 are strongly associated with autism spectrum disorder (ASD). To investigate their effect on synaptic connectivity, we generated cortical neurons from induced pluripotent stem cells (iPSC) derived from neurotypic and ASD-affected donors. We developed Sparse coculture for Connectivity (SparCon) assays where SHANK2 and control neurons were differentially labeled and sparsely seeded together on a lawn of unlabeled control neurons. We observed striking increases in dendrite length, dendrite complexity, total synapse number, and frequency of spontaneous excitatory postsynaptic currents. These findings were phenocopied in gene-edited homozygous SHANK2 knockout cells and rescued by gene correction of an ASD SHANK2 muation, supporting a role for SHANK2 as a regulator of connectivity in developing human neurons. Dendrite length increases were exacerbated by IGF1, TG003, or BDNF, and suppressed by DHPG treatment. The transcriptome in these isogenic SHANK2 neurons was deeply perturbed in synapse, plasticity, and neuronal morphogenesis gene sets and ASD gene modules, and activity-dependent dendrite extension was impaired. Our unexpected findings provide evidence for hyperconnectivity and profoundly altered transcriptome in SHANK2 neurons derived from ASD subjects.

  Study EGAS00001003436

• Pediatric Whole Genome Sequencing Diagnostic Utility

  The standard of care for first-tier clinical investigation of the etiology of congenital malformations and neurodevelopmental disorders is chromosome microarray analysis (CMA) for copy number variations (CNVs), often followed by gene(s)-specific sequencing searching for smaller insertion-deletions (indels) and single nucleotide variant (SNV) mutations. We compared diagnostic rate of WGS to CMA and targeted gene testing of 100 patients referred to The Hospital for Sick Children Genetics clinic in 2014. WGS identified genetic variants meeting clinical diagnostic criteria in 34% of cases, representing a 4-fold increase in diagnostic rate over CMA (8%) (p-value = 1.42e-05) alone and >2-fold increase in CMA plus targeted gene sequencing (13%) (p-value = 0.0009). WGS identified all rare clinically significant CNVs that were detected by CMA. In an additional 26 patients, WGS revealed indel and missense mutations presenting in a dominant (63%) or a recessive (37%) manner. We found four subjects with mutations in at least two genes associated with distinct genetic disorders, including two cases harboring a pathogenic CNV and SNV. When considering medically actionable secondary findings in addition to primary WGS findings, 38% of patients would benefit from genetic counselling.

  Study EGAS00001001623

• Spatial multi-omic map of human myocardial infarction

  Myocardial infarction is a leading cause of mortality worldwide 1 . While advances have been made in acute treatment, an incomplete understanding of remodelling processes has limited the effectiveness of therapies to reduce late-stage mortality 2 . Here, we generate an integrative high- resolution map of human cardiac remodelling after myocardial infarction using single-cell gene expression, chromatin accessibility, and spatial transcriptomic profiling of multiple physiological zones at distinct time points in myocardium from myocardial infarction and control patients. Multimodal data integration allowed us to evaluate cardiac cell-type compositions at increased resolution, yielding insights into changes of the cardiac transcriptome and epigenome through the identification of distinct tissue structures of injury, repair and remodelling. We identified and validated disease-specific cardiac cell-states of major cell-types and analysed them in their spatial context, evaluating their dependency on other cell-types. Our data elucidates molecular principles of human myocardial tissue organisation, recapitulating a gradual cardiomyocyte and myeloid continuum following ischemic injury. In total, our study provides an integrative molecular map of human myocardial infarction and represents an essential reference for the field and paves the way for advanced mechanistic and therapeutic studies of cardiac disease.

  Study EGAS00001006330

• scRNA-seq of HSPC treated with gemcitabine and carbplatin

  The study used cryopreserved hematopoietic stem and progenitor cells (HSPCs) harvested from one patient with chronic myelogenous leukemia (CML). The cells were thawed and acclimatized for 48 hours in StemMACS HSC Expansion Media XF (Miltenyi Biotec) supplemented with StemMACS HSC Expansion Cocktail (Miltenyi Biotec) at 37°C in a humidified atmosphere containing 5% CO2. 4 samples of 2000000 HSPCs in 2 ml media each were then prepared and treated for 24 hours: carboplatin-high: 150 µg/ml carboplatin carboplatin-low: 18.75 µg/ml carboplatin gemcitabine: 25 ng/ml gemcitabine control: no drug Single-cells were subsequently extracted using the ddSEQ™ Single-Cell Isolator (Bio-Rad) and later sequenced using the SureCell™ Whole Transcriptome Analysis 3' Library Prep Kit (Illumina) on the NextSeq 500 System (Illumina) using the NextSeq 500/550 High Output Kit v2.5 150 Cycles (Illumina) all following the manufacturer's instructions. FASTQ-files were then processed following the ddSeeker (Romagnoli, D., et al., BMC Genomics 2018, doi:10.1186/s12864-018-5249-x) and Drop-seq (Macosko, E.Z., et al., Cell 2015, doi:10.1016/j.cell.2015.05.002) protocols for processing scRNA-seq data to yield the final digital gene expression for each cell.

  Study EGAS00001004381

• Functional single-cell characterization of immune aplastic anemia shows convergence of NK and NK-like CD8+ T cells with disease-associated TCR signature

  Immune aplastic anemia (AA) is a life-threatening bone marrow (BM) failure disorder driven by an autoimmune T cell attack against the hematopoietic stem and progenitor cells (HSPCs). However, the autoantigen targets and the role of other immune cells are unknown. Here, by analyzing an international cohort of 156 AA patients with scRNA+TCRαβ-seq, TCRβ-seq, flow cytometry, and plasma cytokine profiling with comparisons to healthy controls and patients with other hematological disorders, we identify NK cells and CD8+ TEMRA cells expressing NK receptors with AA-associated TCRβ-seq motifs as the most dysregulated immune cell populations in AA BM. Functional co-culture experiments with scRNA-TCRαβ-seq readout utilizing primary HSPCs and immune cells provide evidence that NK cellscannot kill HSPCs alone, but sensitize them to CD8+ T cell-mediated killing. Furthermore, HSPCs induce strong activation of T cell clones with CD8+ TEMRA NK-like phenotype and AA-associated TCR motifs. Our results reveal the convergent evolution of innate and adaptive immune cells in AA, where NK cells support CD8+ T cell-mediated autoimmunity.

  Study EGAS00001007602

• Distinct portrayal of lesions in synchronous multifocal lung adenocarcinoma revealed by genome sequencing

  Distinguishing multiple primary lung cancers in the synchronous multifocal intrapulmonary lesions has important significance on clinical staging and therapeutic decision. To investigate genomic aberration profiles, we applied whole genome and whole exome sequencing, and microarray-based comparative genomic hybridization on 15 intrapulmonary tumors derived from six patients with synchronous multifocal lung cancers having similar histological diagnosis. Any pair of intrapulmonary tumors in a single patient, which shared the identical genetic background and environment, showed an extinctive heterogeneity between each other. Phylogenetic relationship analysis indicated an independently branched evolution among all the tumors, suggesting they were multiple primary lung cancers. EGFR or KRAS mutations were found in 7 or 3 out of the 15 tumors, from 3 or 2 patients, respectively. Somatic mutational heterogeneity of these two genes in a single patient was also observed. Our analysis indicates genomic aberration profiling is valuable for identification of multiple primary lung cancer, especially when high histopathological concordance was observed between lesions. We also suggest a thoroughly molecular diagnosis against therapeutic target genes should be taken for each accessible nodule before making a plan for adjuvant therapy.

  Study EGAS00001001572

• SmMIP-tools:a computational toolset for processing and analysis of single-molecule molecular inversion probes derived data

  Single-molecule molecular inversion probes (smMIPs) provides a modular and cost-effective platform for high-multiplex targeted next-generation sequencing (NGS). Nevertheless, translating the raw smMIP-derived sequencing data into accurate and meaningful information currently requires proficient computational skills and a large amount of computational work, prohibiting wide-scale adoption of smMIP-based technologies. To enable easy, efficient, and accurate interrogation of smMIP-derived data, we developed SmMIP-tools, a computational toolset that combines the critical analytic steps for smMIP data interpretation into a single computational pipeline. Here, we describe in detail two major components of the software. The first is a read processing tool that performs quality control steps, generates read-smMIP linkages and retrieves molecular tags. The second is an error-aware variant caller capable of detecting single nucleotide variants (SNVs) and short insertions and deletions (indels). Using a cell-line DNA dilution series and a cohort of blood cancer patients, we benchmarked SmMIP-tools and evaluated its performance against clinical sequencing reports. We anticipate that SmMIP-tools will increase accessibility to smMIP-technology, enabling cost-effective genetic research to push personalized medicine forward.

  Study EGAS00001005359

• Transcriptomes of human CD4+ T lymphocytes - Metabolic project

  Human CD4+ T cells are essential mediators of immune responses. By altering the mitochondrial and metabolic states, we defined metabolic requirements of human CD4+ T cells for in vitro activation, expansion, and effector function. T cell activation and proliferation were reduced by inhibiting oxidative phosphorylation, while early cytokine production was maintained by either OXPHOS or glycolytic activity. Glucose deprivation in the presence of mild mitochondrial stress markedly reduced all three T cell functions, contrasting the exposure to resveratrol, an antioxidant and sirtuin-1 activator, which specifically inhibited cytokine production and T cell proliferation, but not T cell activation. Conditions that inhibited T cell activation were associated with the downregulation of 2′,5′-oligoadenylate synthetase genes via interferon response pathways. Our findings indicate that T cell function is grossly impaired by stressors combined with nutrient deprivation, suggesting that correcting nutrient availability, metabolic stress and/or the function of T cells in these conditions will improve the efficacy of T cell-based therapies.

  Study EGAS00001005565

• Single-nucleus transcriptomic profiling of aging Down Syndrome brains

  Down syndrome (DS) – trisomy of human chromosome 21 (HSA21) – is characterized by lifelong cognitive impairments and development of neuropathological hallmarks of Alzheimer’s disease (AD). The cellular and molecular modifications responsible for these effects are not understood. Here we performed single-nucleus RNA-sequencing (snRNA-seq) employing both short (Illumina) and long-read (Pacific Biosciences) sequencing technologies on a total of 29 DS and non-DS control prefrontal cortex samples. In DS, the ratio of inhibitory-to-excitatory neurons was significantly increased, which was not observed in previous reports examining sporadic AD. DS microglial transcriptomes displayed AD-related aging and activation signatures in advance of AD neuropathology, with increased microglial expression of C1q complement genes (associated with dendritic pruning) and the HSA21 transcription factor gene RUNX1. Long-read sequencing detected vast RNA isoform diversity within and amongst specific cell types including numerous novel sequences that differed between DS and normal brains. Notably, over 8,000 genes produced RNAs containing intra-exonic junctions (IEJs), including APP that had previously been associated with somatic gene recombination. These and related results illuminate large-scale cellular and transcriptomic alterations as novel features of the aging DS brain.

  Study EGAS00001005691

• Increased mutation accumulation during fetal development in Down syndrome

  Children with Down syndrome (DS) show a high incidence of leukemia, which is preceded by transient myeloproliferative disorder (TMD) in neonates. TMD is characterized by mutations in GATA1 that accumulate during fetal development, indicating that additional mutagenesis is required for leukemogenesis. Abnormalities in chromosome numbers can result in a mutator phenotype. However, it is unknown whether a trisomy of chromosome 21 (T21) results in enhanced mutation accumulation during human development. Here, we determined the mutational consequences of a constitutive T21 in DS fetuses compared to disomy 21 cells (D21) of D21 fetuses. We found that D21 fetal HPSCs had a higher somatic mutation rate compared to healthy adults HSPCs. Moreover, we observed an increase in somatic mutation load in T21 cells of DS fetuses, which was already apparent during the first cell divisions of embryogenesis. The mutation load could be explained by processes that are active during normal fetal development and were similar to those observed in TMD blast cells. Our study shows that T21 causes an increase in somatic load in fetal cells of DS fetuses, which may contribute to leukemia development early in life.

  Study EGAS00001003982

• The DNA methylation landscape of glioblastoma disease progression shows extensive heterogeneity in time and space

  Glioblastoma is characterized by widespread genetic and transcriptional heterogeneity, yet little is known about the role of the epigenome in glioblastoma disease progression. Here, we present genome-scale maps of DNA methylation in matched primary and recurring glioblastoma tumors, using data from a highly annotated clinical cohort that was selected through a national patient registry. We demonstrate the feasibility of DNA methylation mapping in a large set of routinely collected FFPE samples, and we validate bisulfite sequencing as a multipurpose assay that allowed us to infer a range of different genetic, epigenetic, and transcriptional characteristics of the profiled tumor samples. On the basis of these data, we identified subtle differences between primary and recurring tumors, links between DNA methylation and the tumor microenvironment, and an association of epigenetic tumor heterogeneity with patient survival. In summary, this study establishes an open resource for dissecting DNA methylation heterogeneity in a genetically diverse and heterogeneous cancer, and it demonstrates the feasibility of integrating epigenomics, radiology, and digital pathology for a national cohort, thereby leveraging existing samples and data collected as part of routine clinical practice.

  Study EGAS00001002538

• Cross-tissue transcriptomic analysis of human secondary lymphoid organ residing ILC3 reveals a default quiescent state in the absence of inflammation

  A cross tissue transcriptional comparison of human ILC3 isolated from non-inflamed lymph nodes (LN) and spleens, inflamed tonsils and the peripheral blood (PB). Using uniform gating and high-purity cell-sorting, ILC3 were sorted as Lineage(CD3,CD19,CD34,CD14,CD94)-CD45intCD127+CRTH2-CD117+. From LN, spleen, tonsils and PB, NKp44- ILC3 were isolated and, in addition, NKp44+ ILC3 were sorted from spleen and tonsils. RNA quality was measured on a RNA 6000 Pico Kit (Agilent Technologies) using a 2100 Bioanalyzer (Agilent Technologies). 100 pg high-quality RNA was isolated and cDNA was generated using the SMART-Seq (v3) Ultra Low Input RNA Kit (Clontech) with 15 cycles of amplification. Amplified and Covaris-fragmented cDNA was quality-checked on a High Sensitivity DNA chip (Agilent Technologies) and further processed according the TruSeq Nano DNA Library Preparation Kit (Illumina). Data was processed as follows: SMARTer adapters were trimmed with cutadapt and the resulting sequences were aligned to the human RefSeq transcriptome using 14 TopHat2 (Kim et al., 2013). Normalization and quantification was performed using Cufflinks (Trapnell et al., 2012). FPKM counts were determined per gene with HTSeq-count.

  Study EGAS00001002636

• The identification of genetic vulnerabilities in head and neck cancers for the development of novel treatments.

  The global cancer burden is expected to double by the year 2030 and this will disproportionately affect the Asian continent. Head and neck squamous cell carcinoma (HNSCC) is endemic in Asia. It is the sixth most common cancer worldwide and the second most common in Malaysia, with approximately 600,000 new cases per year globally. Longterm survival remains less than 50%, underscoring the need for better therapeutic strategies. In this proposal we aim to bring new technological advancements in the identification of critical genes in HNSCC and the identification of compounds for the development of new therapies in HNSCC. We propose to genomically characterise 16 unique Asian H&N cancer cell lines, to classify sub groups where identification of essential genes that can be targeted for tumor control can be applied in a more rational way. Further, we propose to use zebrafish, an established drug discovery platform to identify promising pathway-specific compounds that could be further studied using CRISPR/Cas9 lethality screens and biological models to better identify novel therapeutic strategies to combat this major health problem in Asia.

  Study EGAS00001002682

• Study on the proliferation history of colorectal adenomas

  Mismatch repair (MMR) deficient colorectal adenomas are composed of transformed cells that descend from a common founder and progressively accumulate genomic alterations but their proliferation history is still largely unknown. Here we rebuilt the proliferation history of four MMR deficient colorectal adenomas from male individuals by deep sequencing the X chromosome and mapping the somatic mutations that were progressively acquired during tumor growth. We developed a novel method to call high and low frequency mutations, to cluster them according to their frequencies and rebuild the proliferation trees directly from these mutation clusters using a recursive algorithm. The trees of all four lesions were formed of a dominant subclone that coexisted with other genetically heterogeneous subpopulations of cells. Despite this similar hierarchical organization, however, the growth dynamics varied among and within tumors, likely depending on a combination of tumor-specific genetic and environmental factors. Our study provides insights into the biological properties of individual MMR deficient colorectal adenomas that may influence their growth and also the response to therapy. Extended to other solid tumors, our approach could inform on the mechanisms of cancer progression, thus providing insight into the best choice of cancer treatment methods.

  Study EGAS00001000883

• UK10K_OBESITY_GS

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The GS samples are part of the Obesity group and will undergo exome sequencing.

  Study EGAS00001000242

• The mission of the BIOS Consortium is to create a large-scale data infrastructure and to bring together BBMRI researchers focusing on integrative omics studies in Dutch Biobanks.

  The advent of the genome-wide association study (GWAS) led to the successful identification of thousands of variants that are robustly associated with complex disease phenotypes. Dutch biobanks played a substantial role in these discoveries. For most of these variants, however, the mechanisms through which they contribute to these phenotypes remain unknown. The BIOS Consortium applies a functional genomics approach that integrates genome-wide genetic data with data on the epigenome and transcriptome to elucidate these mechanisms. Over 4000 samples from BBMRI-NL biobanks with in-depth information on disease phenotypes and GWAS data are being enriched with RNA-sequencing and genome-wide DNA methylation data. The same is true for samples with whole-genome sequencing data from GoNL. This unique data infrastructure provides a powerful platform to evaluate key questions in integrative omics from establishing comprehensive eQTL and meQTL catalogues to linking molecular pathways across omics levels to phenotypic outcomes.

  Study EGAS00001001077

• Comprehensive characterization of cell-free tumor DNA in plasma and urine of patients with renal tumors

  Cell-free tumour-derived DNA (ctDNA) allows non-invasive monitoring of various cancers but its utility in renal cell cancer (RCC) has not been well established. A combination of untargeted and targeted sequencing methods, applied to two independent cohorts of patients (n=90) with various renal tumour subtypes, were used to determine ctDNA content in plasma and urine. Our data revealed lower plasma ctDNA levels in RCC relative to other cancers, with untargeted detection of ~33% for both cohorts. A highly sensitive personalised approach, applied to plasma and urine from select patients (n=22), improved detection to ~50%, including in patients with early stage and even benign lesions. A machine-learning based model, applied to untargeted data, predicted this detection, potentially offering a means of triaging patient samples for personalised analysis. We observed that plasma, and for the first time, urine ctDNA may better represent tumour heterogeneity than a single tissue biopsy. Longitudinal sampling of >200 plasma samples revealed that ctDNA can track disease course. These data highlight low ctDNA levels in RCC but indicate potential clinical utility provided improvement in isolation and detection approaches.

  Study EGAS00001003530

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

  Cell-free DNA (cfDNA) in the blood provides a noninvasive diagnostic avenue for patients with cancer. However, characteristics of the origins and molecular features of cfDNA are poorly understood. We developed an approach to evaluate fragmentation patterns of cfDNA across the genome and found that cfDNA profiles of healthy individuals reflected nucleosomal patterns of white blood cells, while patients with cancer had altered fragmentation profiles. We applied this method to analyze fragmentation profiles of 236 patients with breast, colorectal, lung, ovarian, pancreatic, gastric, or bile duct cancers and 245 healthy individuals. A machine learning model incorporating genome-wide fragmentation features had sensitivities of detection ranging from 57% to >99% among the seven cancer types at 98% specificity, with an overall AUC of 0.94. Fragmentation profiles could be used to identify the tissue of origin of the cancers to a limited number of sites in 75% of cases. Combining our approach with mutation based cfDNA analyses detected 91% of cancer patients. The results of these analyses highlight important properties of cfDNA and provide a proof of principle approach for screening, early detection, and monitoring of human cancer.

  Study EGAS00001003611

• Genomic analysis of Smoothened inhibitor resistance in basal cell carcinoma

  Smoothened inhibitors are currently being investigated for the treatment of several cancers. Vismodegib is approved for the treatment of locally advanced and metastatic basal cell carcinoma (BCC). The majority of BCC patients treated with vismodegib experience significant clinical benefit, however, a small number of patients develop resistance. Knowledge of resistance mechanisms can generate predictive biomarkers and is critical for the design of additional therapeutic strategies aimed at circumventing resistance. To investigate mechanisms of resistance to vismodegib in BCC we sequenced and profiled biopsies from patients who initially responded to treatment and subsequently progressed on drug.We find that resistance is associated with re-activation of the Hedgehog (Hh) pathway through diverse mechanisms, including concurrent copy number alterations in the downstream Hh pathway components SUFU and GLI2, as well as novel SMO mutations that reduce sensitivity to vismodegib. The latter fall into two classes: those that activate SMO and those that directly affect the vismodegib-binding pocket. Additionally, we observe evidence of clonal selection of SMO mutations and intra-tumor heterogeneity of resistance mechanisms, which implies diverse and combinatorial strategies are required to overcome resistance.

  Study EGAS00001000845

• Germline mutations in the transcription factor IKZF5 cause thrombocytopenia.

  To identify novel causes of hereditary thrombocytopenia, we performed a genetic association analysis of whole-genome sequencing data from 13,037 individuals enrolled in the NIHR BioResource, including 233 cases with isolated thrombocytopenia. We found an association between rare variants in the transcription factor encoding gene IKZF5 and thrombocytopenia. We report five causal missense variants in or near IKZF5 zinc fingers (Znfs), of which two occurred de novo and three co-segregated in three pedigrees. A canonical DNA-Znf binding model predicts that three of the variants alter DNA recognition. Expression studies showed that chromatin binding was disrupted in mutant compared to wild-type (WT) IKZF5 and electron microscopy revealed a reduced quantity of alpha granules in normally sized platelets. Proplatelet formation (PPF) was reduced in megakaryocytes (MKs) from seven cases relative to six controls. Comparison of RNA-seq data from platelets, monocytes, neutrophils and CD4 T-cells from three cases and 14 healthy controls showed 1,194 differentially expressed genes in platelets but only four DEGs in each of the other blood cell types. In conclusion, IKZF5 is a novel transcriptional regulator of megakaryopoiesis and the eighth transcription factor associated with dominant thrombocytopenia in humans.

  Study EGAS00001003736

• RNA sequencing and Illumina 2.5M SNP array data collected from 675 commonly used human cancer cell lines.

  Tumor-derived cell lines have served as vital models to advance our understanding of oncogene function and therapeutic response1. Although substantial effort has been directed to defining the genomic constitution of cancer cell line panels2–4, the transcriptome – which represents the active program of a cell – remains understudied. Here, we describe RNA sequencing and SNP array analysis of 675 commonly used human cancer cell lines. We explore numerous transcriptome features including coding and non-coding gene expression, transcribed mutations, gene fusion and expression of non-human sequences. Aside from many known aberrations we find new surprising characteristics, including more than 2200 unique fusion gene pairs representing a vast, testable repertoire of oncogenic fusions, many of which have analogs found in primary human tumors. We show that a combination of multiple genome and transcriptome features in a novel pathway-based approach enhances prediction of response to various targeted therapeutics. Our results provide valuable new insights into these critical pre-clinical models and provide added context for interpreting the numerous studies that employ these widely used cell lines.

  Study EGAS00001000610

• Reference Exome Data for a Northern Brazilian population

  Whole exome sequencing (WES) is a popular technology widely used in the diagnosis of rare genetic diseases and provides useful variant data for analysis of complex diseases. However, there is not always adequate population-specific reference data available to assist clinical geneticists in assigning a diagnostic variant to a specific clinical condition. Here we provide a catalogue of variants called after sequencing the exomes of 45 babies from the State of Rio Grande do Nord in Brazil. Sequence data were processed using an ‘intersect-then-combine’ (ITC) approach, using GATK and SAMtools to call variants. A total of 612,761 variants were identified in at least one individual in this Brazilian Cohort, including 559,448 single nucleotide variants (SNVs) and 53,313 insertion/deletions. Of these, 58,111 overlapped with nonsynonymous (nsSNVs) or splice site (ssSNVs) SNVs in dbNSFP. As an aid to clinical diagnosis of rare diseases, we used the American College of Medicine Genetics and Genomics (ACMG) guidelines to assign pathogenic/likely pathogenic status to 439 (0.755%) of the 58,111 nsSNVs and ssSNVs. Our data set provides a useful reference point for diagnosis of rare diseases in Brazil.

  Study EGAS00001004112

• Oncoprint GSCCs

  Glioblastoma (GBM) remains the most malignant primary brain tumor, with a median survival rarely exceeding 2 years. Tumor heterogeneity and an immunosuppressive microenvironment are key factors contributing to the poor response rates of current therapeutic approaches. GBM-associated macrophages (GAMs) often exhibit immunosuppressive features that promote tumor progression. However, their dynamic interactions with GBM tumor cells remain poorly understood. Here, we used patient-derived GBM stem cell cultures and combined single-cell RNA sequencing of GAM-GBM co-cultures and real-time in vivo monitoring of GAM-GBM interactions in orthotopic zebrafish xenograft models to provide insight into the cellular, molecular, and spatial heterogeneity. Our analyses revealed substantial heterogeneity across GBM patients in GBM-induced GAM polarization and the ability to attract and activate GAMs – features that correlated with patient survival. Differential gene expression analysis, immunohistochemistry on original tumor samples, and knock-out experiments in zebrafish subsequently identified LGALS1 as a primary regulator of immunosuppression. Overall, our work highlights that GAM-GBM interactions can be studied in a clinically relevant way using co-cultures and avatar models, while offering new opportunities to identify promising immune-modulating targets.

  Study EGAS00001007481

• Spatial_transcriptome_analysis_of_Paediatric_Thymus

  Cell atlas of Paediatric Thymus Cell atlas of Paediatric Thymus - The Human Cell Atlas initiative is an ambitious global research effort which aims to describe every cell in the human body across all developmental stages to generate a reference map to accelerate progress in science and medicine particularly relating to human development and disease. The initiative is led by investigators based at the Wellcome Sanger Institute and the Broad Institute of MIT and Harvard and funding support from the Wellcome and NIH. We aim to interrogate the cellular composition and molecular regulators underpinning a cells identity through; genetic and genomic analysis, expression analysis, non-identifiable video imaging and cell line and organoid production, as well as, computational algorithms to predict cellular developmental trajectories and cell-cell interactions. We will subsequently validate these findings in situ using spatial genomic and imaging technologies e.g. immunohistochemistry, RNAscope, in situ sequencing and spatial transcriptomic analyses. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001004281

• Widespread DNA hypomethylation and differential gene expression in Turner syndrome

  Adults with 45,X monosomy (Turner syndrome) reflect a surviving minority since more than 99% of fetuses with 45,X monosomy die in utero. In adulthood 45,X monosomy is associated with increased morbidity and mortality, although strikingly heterogeneous with some individuals left untouched while others suffer from cardiovascular disease, autoimmune disease and infertility. The present study investigates the leukocyte DNAmethylation profile by using the 450K-Illumina Infinium assay and the leukocyte RNA-expression profile in 45,X monosomy compared with karyotypically normal female and male controls. We present results illustrating that genome wide X-chromosome RNA-expression profile, autosomal DNA-methylation profile, and the X-chromosome methylation profile clearly distinguish Turner syndrome from controls. Our results reveal genome wide hypomethylation with most differentially methylated positions showing a medium level of methylation. Contrary to previous studies, applying a single loci specific analysis at well-defined DNA loci, our results indicate that the hypomethylation extend to repetitive elements. We describe novel candidate genes that could be involved in comorbidity in TS and explain congenital urinary malformations (PRKX), premature ovarian failure (KDM6A), and aortic aneurysm formation (ZFYVE9 and TIMP1)

  Study EGAS00001002190

• Determination_of_cell_specific_regulatory_enhancers_in_hematopoetic_models

  The BLUEPRINT project is a large-scale project investigating epigenetic mechanisms involved in blood formation, in health and disease. The human variation workpackage (WP10) of the project seeks to characterize the effect of common sequence variation on the epigenome status of a cell. To do this, the project will use highly purified blood cells to minimise “experimental noise” and therefore enhance the power to discover modest effects. Two peripheral blood cell types, the CD14+CD16- monocyte (an important central orchestrator of adaptive immunity and a bridge between innate and adaptive immunity) and the CD65+CD9- neutrophilic granulocyte (the frontline cell for innate immunity) have been selected for this purpose. The two types of cells will be obtained at high purity from adult blood (AB) of 200 healthy males and females, respectively. Cells will be purified by using already validated and fully operational protocols that are based on density gradient centrifugation of the buffy coat obtained from whole blood, followed by magnetic bead-based purification using monoclonal antibodies against Cluster of Differentiation (CD) lineage-specific cell surface markers. This data set contains functional genomics data for gene expression and chromatin state.

  Study EGAS00001000586

• We performed whole-exome sequencing of 20 samples (10 actinic keratosis and 10 cutaneous squamous cell carcinoma) to investigate a potential relationship between DNA methylation-based subtypes and genetic mutation patterns (Rodriguez-Paredes et al., Nat Commun 2017)

  Cutaneous squamous cell carcinoma (cSCC) is the second most common skin cancer type and arises from keratinocytes. Most cSCC progress from a UV-induced precancerous lesion termed actinic keratosis (AK). Despite various efforts to characterize these lesions molecularly, the etiology of AK and its progression to cSCC remain only partially understood. Here we have used Infinium MethylationEPIC BeadChips to interrogate the DNA methylation status of about 850.000 CpGs in epidermal preparations from healthy skin, AK and cSCC. Importantly, we found that the premalignant AK samples displayed classical features of cancer methylomes and were highly similar to cSCC methylomes. Further analysis identified typical features of stem cell methylomes, such as a reduced DNA methylation age, non-CpG methylation and stem cell-related keratin and enhancer methylation patterns. Interestingly, this signature was detected only in one half of the AK and cSCC samples, while the other half showed keratin and enhancer methylation patterns that were more closely related to the control epidermis. These findings suggest the existence of two distinct subclasses of AK and cSCC that originate from distinct keratinocyte differentiation stages.

  Study EGAS00001002670

• Gut microbiome modulates response to anti PD1 immunotherapy in metastatic melanoma patients

  There is a growing appreciation of the role of the microbiome in cancer, and evidence in pre-clinical models that the gut microbiome may modulate responses to immune checkpoint blockade though this has not been well-characterized in patients. We analyzed the oral (n=86)and gut (n=43) 16S microbiome in melanoma patients on PD-1 blockade. Significant differences were noted in the diversity and composition of the gut microbiome between responders and non-responders in patients with a fecal microbiome sample, with significantly higher alpha diversity and relative abundance of Ruminococcaceae bacteria) in R. Metagenomic studies (n=25) revealed functional differences in gut bacteria in R including enrichment of anabolic pathways. Immune profiling demonstrated enhanced systemic and anti-tumor immunity in patients with a favorable gut microbiome as well as in germ-free mice receiving fecal transplant from responding patients. Together, these data have important implications for treatment with immune checkpoint blockade in cancer. The fastq files associated with this dataset are stored at ENA under the following accession numbers: PRJEB22894 (Fecal 16S), PRJEB22874 (Oral 16S), PRJEB22895 (Murine 16S), PRJEB22893 (Fecal WGS).

  Study EGAS00001002698

• Germline HAVCR2 mutations altering TIM-3 characterize subcutaneous panniculitis-like T-cell lymphomas with hemophagocytic lymphohistiocytic syndrome

  Sub-cutaneous panniculitis-like T-cell lymphomas (SPTCL), a non-Hodgkin lymphoma, can be associated with hemophagocytic lymphohistiocytosis (HLH), a life-threatening immune activation which adversely impacts survival1,2. T-cell-immunoglobulin mucin-3 (TIM-3) is a modulator of immune responses expressed on subgroups of T- and innate immune cells. We identify in ~60% of SPTCL cases germline, loss-of-function, missense variants in highly conserved residues of TIM-3, c.245A>G(p.Y82C) and c.291A>G(p.I97M), each with specific geographic distribution. Y82C-TIM-3 occurs on a potential founder chromosome in patients with East Asian and Polynesian ancestry, while I97M-TIM-3 occurs in Caucasians. Both variants induce protein misfolding and abrogate TIM-3’s plasma-membrane expression leading to persistent immune activation, increased production of inflammatory cytokines, including TNF-α and IL-1β, promoting HLH and SPTCL. Our findings highlight HLH/SPTCL as a new genetic entity and showcase TIM-3 mutations as the first causative genetic defect in SPTCL. While TIM-3-mutant HLH/SPTCL benefit from immunomodulation, therapeutic repression of the TIM-3 checkpoint may have adverse consequences.

  Study EGAS00001002765

• SINGLE_CELL_ANALYSIS_OF_IN_VITRO_ERYTHROPOIESIS

  We are developing a protocol to differentiate mouse and human induced pluripotent stem (IPS) and embryonic stem (ES) cells towards the haematopoietic pathway to generate erythrocytes in vitro. This system has many applications such as the study of the role of specific genes and human polymorphisms in infectious diseases such as malaria, as well as haematological diseases such as myelodysplastic syndrome. The nature of the in vitro differentiation process means that a heterogeneous population of cells is generated. In order to understand the types of cells produced with our protocol, we have performed a single cell analysis, which has the power to reveal the different populations of cells and their characteristics. For this, a cDNA library has been made that needs to be sequenced to obtain the gene expression profiles of the different cells. With this information we will be able to assess the quality of the differentiation protocol and improve it in order to produce better cells for the downstream applications.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 EGAS00001000764

• Whole genome sequencing of Ewings Sarcoma

  Sarcomas are cancers of the bone and soft tissue often defined by their gene fusions. However, the timing, context, and processes by which these pathogenic fusions arise are unknown. We explored this in Ewing sarcoma, a cancer driven by EWSR1-ETS fusions, with very few cooperating mutations. Combining whole-genome sequencing with enhanced informatics, we found that the EWSR1-ETS fusion arose from striking rearrangement clusters in 42% of cases (52/124). Notably, these were organized in loops that universally contained the fusion at their center, while also weaving up to 18 genes together with it. We found the same pattern of rearrangements in three additional types of sarcoma. From these data, we define a new signature for sarcoma fusions that precedes other somatic changes, in the earliest replicating DNA of the genome. This dramatic, sudden process impinges on many genes – generating multiple coding changes that profoundly affect the transcriptome, with the disease-defining gene fusion at its core. These rearrangement loops emerge in an early ES clone from which both the primary tumor and the lethal relapse emerged, and then evolved in parallel until clinically detected.

  Study EGAS00001003385

• Onco-exaptation of an Endogenous Retroviral LTR Drives IRF5 Expression in Hodgkin Lymphoma

  The transcription factor interferon regulatory factor 5 (IRF5) is upregulated in Hodgkin lymphoma (HL) and a key regulator of the aberrant transcriptome characteristic of this disease. Here we show that IRF5 upregulation in HL is driven by transcriptional activation of a normally dormant endogenous retroviral LOR1a long terminal repeat (LTR) upstream of IRF5. Specifically, through screening of RNA-seq libraries, we detected LTR-IRF5 chimeric transcripts in multiple HL cell lines but not in normal B-cell controls. In HL, the LTR was in an open and hypomethylated epigenetic state and we further show the LTR is the site of transcriptional initiation. Among HL cell lines, usage of the LTR promoter strongly correlates with overall levels of IRF5 mRNA and protein, indicating that LTR transcriptional awakening is a major contributor to IRF5 upregulation in HL. Taken together, oncogenic IRF5 overexpression in HL is the result of a specific LTR transcriptional activation. We propose that such LTR derepression is a distinct mechanism of oncogene activation (“onco-exaptation”) and that such a mechanism warrants further investigation in molecular and cancer research.

  Study EGAS00001001205

• Dissection of the molecular complexity of colorectal cancer in pre-clinical models identifies predictive signatures of sensitivity to EGFR inhibitors

  Colorectal carcinoma (CRC) is a heterogeneous tumour entity, with only a fraction of the tumours responding to currently available therapies. A detailed understanding of the biology and genetics of CRC is a prerequisite for precision medicine. To address this challenge in clinical oncology, the OncoTrack consortium recruited 106 CRC patients (stages I-IV) and developed an integrated pre-clinical platform analysing drug sensitivity in the established matched patient-derived in vivo and in vitro models, 46 three-dimensional cell cultures and 59 xenografts, recapitulating many of the molecular features of the donor tumours. We generated extensive omics information on patient tumours and their derived models used to generate a compendium of drug sensitivity data. We could distinguish tumour molecular sub-groups responding differently to therapy and identified a novel independent predictor of sensitivity to the EGFR inhibitor cetuximab. The findings reported here will advance our understanding of CRC and its response to therapy.

  Study EGAS00001001752

• Whole genome and transcriptome analysis of medullary thyroid cancer

  Medullary thyroid cancer is the neoplasm of the calcitonin-producing parafollicular cells of the thyroid gland. The hormone calcitonin acts to lower the blood calcium level, counteracting the effects of the parathyroid hormone; it can also serve as a biomarker for the presence or recurrence of medullary thyroid cancer. Surgery is the only curative treatment for this cancer. Targeted therapeutics vandetanib and cabozantinib are approved for the treatment of aggressive and metastatic tumors that are not amenable to surgery. The use of these multikinase inhibitors are supported by the observed overactivation of the RET oncoprotein in a large sub-population of medullary thyroid cancers; however, not all patients carry oncogenic alterations of this kinase. Hence, there is a need for molecular characterization of medullary tumors on the genomic scale with the aim of identifying targetable mutations. Here, we describe the genomic profiles of two medullary thyroid cancers and report the presence of a putative oncogenic BRAF fusion in one. Such alterations previously observed in other malignancies and known targets of available drugs, can benefit patients who currently have no therapy options.

  Study EGAS00001001473

• Functional single-cell characterization of immune aplastic anemia shows convergence of NK and NK-like CD8+ T cells with disease-associated TCR signature

  Immune aplastic anemia (AA) is a life-threatening bone marrow (BM) failure disorder driven by an autoimmune T cell attack against the hematopoietic stem and progenitor cells (HSPCs). However, the autoantigen targets and the role of other immune cells are unknown. Here, by analyzing an international cohort of 156 AA patients with scRNA+TCRαβ-seq, TCRβ-seq, flow cytometry, and plasma cytokine profiling with comparisons to healthy controls and patients with other hematological disorders, we identify NK cells and CD8+ TEMRA cells expressing NK receptors with AA-associated TCRβ-seq motifs as the most dysregulated immune cell populations in AA BM. Functional co-culture experiments with scRNA-TCRαβ-seq readout utilizing primary HSPCs and immune cells provide evidence that NK cellscannot kill HSPCs alone, but sensitize them to CD8+ T cell-mediated killing. Furthermore, HSPCs induce strong activation of T cell clones with CD8+ TEMRA NK-like phenotype and AA-associated TCR motifs. Our results reveal the convergent evolution of innate and adaptive immune cells in AA, where NK cells support CD8+ T cell-mediated autoimmunity.

  Study EGAS00001007604

• noninvasive lung cancer subtyping

  Accurate diagnosis of lung cancer is important for treatment decision-making. Currently, the gold standard for diagnosing histological subtypes of lung cancer relies on tumor biopsies. Recently, liquid biopsy, particularly cell-free DNA (cfDNA), has shown promising results in cancer detection and classification. In this study, we investigated the potential of cfDNA methylome for the noninvasive classification of lung cancer histological subtypes. Specifically, we focused on the two most prevalent lung cancer subtypes, lung adenocarcinoma and lung squamous cell carcinoma. Using a fragment-based marker discovery approach, we identified robust subtype-specific methylation markers from tumor samples. These markers were successfully validated in independent cohorts and associated with subtype-specific transcription activity. Leveraging these markers, we constructed a subtype classification model using cfDNA methylation profiles, achieving an AUC of 0.808 in the cross-validation and an AUC of 0.747 in the independent validation. Additionally, tumor copy number variations inferred from cfDNA methylome analysis revealed potential implications for treatment selection. In summary, our study demonstrates the potential of cfDNA methylome analysis for noninvasively lung cancer subtyping, offering insights for cancer monitoring and early detection.

  Study EGAS00001007717

• An_exome_sequencing_study_of_the_HIV_elite_long_term_non_progressors_and_rapid_progressors__CASCADE_cohorts_

  Background: A rare subgroup of HIV infected individuals naturally controls infection withouttreatment. These ?elite controllers? constitute an important model for the natural control ofHIV infection. Indeed, the study of these individuals may provide insights into strategies forthe development of HIV vaccines. Although several HLA and chemokine alleles are knownto be over-represented in elite controllers, only a small portion of HIV phenotypic variation isexplained by known genetic variants. The elite controller phenotype is rare and distinct,representing the extreme of an infectious disease trait. As such, this phenotype may be partlyexplained by variation in host immune control, which may be characterized by differences inrare functional genetic variants. Genomic regions underlying elite control can be potentiallyidentified by comparing the presence or frequency of variants in this group to thatrepresenting the opposite extreme. In this context, ?rapid progressors? is a group defined byits rapid immunological and clinical disease progression.Aim: To extend an existing study, in order to identify DNA sequence variants involved in thecontrol of HIV infection with greater statistical resolution. Specifically, we aim to sequence upto 200 exomes from multiple cohort studies within the EuroCoord CASCADE collaboration (acollaboration of 25 HIV seroconversion cohort studies across Europe).

  Study EGAS00001000522

• The Ovarian Cancer Association Consortium OncoArray genome-wide association study

  A custom Illumina genotyping array, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis as described in PMID:27697780.The array was used to genotype ovarian cancer cases and controls from multiple studies participating in the Ovarian Cancer Association Consortium, of which 54,747 passes QC.An initial principal components analysis was carried out using a method described in PMID:27697780 to assign subject to intercontinental ancestry groups: “European”, “African”, “Asian” and "Other". Within ancestry principal components analysis was carried using a set of ~33k unlinked markers. Subsequent analyses were adjusted for the ancestry specific PCs: 9 Eur, 10 Asian, 1 African and 6 Mixed. These PCs are included in the imputed dataset phenotypes file

  Study EGAS00001002305

• Whole Mitochondrial Sequencing of Gingivo-buccal Cancer: ICGC-India Project

  As a part of the ICGC, India has undertaken genomic studies on gingivobuccal cancer of the oral cavity, which is the most prevalent form of cancer among men in India. There are various known environmental (life-style) correlates of this cancer, the most important of which are tobacco chewing and HPV infection. Detailed clinical characterization of the patients, collection of data on demographic and environmental exposures, and isolation of DNA samples from blood and tumor tissues collected from each patient are being done at the Advanced Centre for Research, Treatment and Education on Cancer, Mumbai. Whole genome and other sequencing experiments are being performed at the National Institute of Biomedical Genomics, Kalyani. Whole mitochondria sequencing of 89 patients was performed in this study. The mean depth of sequencing was at least 2500X, but variable across platforms, 2689X and 9894X for amplicon sequencing in Ion Torrent PGM (Tumor-blood pairs=16) and Illumina HiSeq 2500 (Tumor-blood pairs =48), respectively; and, 14,190X for whole genome sequencing in Illumina HiSeq 2500 (Tumor-blood pairs =25). Associations of mutation data with important clinical aspects have been found.

  Study EGAS00001002425

• The demographic history and mutational load of African hunter-gatherers and farmers

  The distribution of deleterious genetic variation across human populations is a key issue in evolutionary biology and medical genetics. However, the impact of different modes of subsistence on recent changes in population size, patterns of gene flow, and deleterious mutational load remains unclear. Here, we report high-coverage exomes from various populations of rainforest hunter-gatherers and farmers from central Africa. We find that the recent demographic histories of hunter-gatherers and farmers differed considerably, with population collapses for hunter-gatherers and expansions for farmers, accompanied by increased gene flow. We show that purifying selection against newly arising deleterious alleles is of similar efficiency across African populations, in contrast with Europeans where we detect weaker purifying selection. Furthermore, the per-individual mutation load of rainforest hunter-gatherers is similar to that of farmers, under both additive and recessive models. Our results indicate that differences in the cultural practices and demographic regimes of African populations have not resulted in large differences in mutational burden, and highlight the beneficial role of gene flow in reshaping the distribution of deleterious genetic variation across human populations.

  Study EGAS00001002457

• Evolutionary analysis of primary tumors and metastatic lesions from 20 breast cancer patients (99 samples in total) using exome sequencing data.

  Metastatic breast cancers are still incurable. Characterizing their evolutionary landscape including the role of metastatic axillary lymph nodes to seed distant organ metastasis can provide rational basis for effective treatments. Here, we describe the genomic analyses of the primary tumors and metastatic lesions from 20 breast cancer patients (99 samples). Our evolutionary analyses revealed diverse spreading and seeding patterns governing tumor progression. Although linear evolution to successive metastatic sites was common, parallel evolution from the primary tumor to multiple distant sites was also evident. Metastatic spreading was frequently coupled with polyclonal seeding, where multiple metastatic subclones originated from primary tumor and/or other distant metastases. Synchronous axillary lymph node metastasis, a well-established prognosticator of breast cancer, was not involved in seeding distant metastasis, suggesting haematogenous route for cancer dissemination. Clonal evolution coincided frequently with emerging driver alterations and evolving mutational processes, notably a significant increase in apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) associated mutagenesis. Our data provide the first genomic evidence regarding the role of axillary lymph node metastasis in seeding distant organ metastasis and elucidates the evolving mutational landscape during cancer progression.

  Study EGAS00001002737

• Defective T-cell expansion in RASGRP1 deficiency

  Inherited CTPS1- CD27- and CD70-deficiencies in humans have revealed key factors of T-lymphocyte expansion, that is a prerequisite for an efficient immunity to Epstein-Barr virus (EBV) infection. RASGRP1 is a T-lymphocyte specific nucleotide exchange factor known to activate the MAP kinases pathway. A deleterious homozygous mutation in RASGRP1 leading to the loss RASGRP1 expression was identified in two siblings who both developed a persistent EBV infection leading to Hodgkin lymphoma. RASGRP1-deficient T cells exhibited defective MAPK activation and impaired proliferation that was restored by expression of wild-type RASGRP1. Similar defects were observed in T cells from healthy individuals when RASGRP1 was downregulated. RASGRP1-deficient T cells also exhibited decreased CD27-dependent proliferation toward CD70-expressing EBV-transformed B cells, a crucial pathway required for expansion of antigen-specific T cells in anti-EBV immunity. Furthermore, RASGRP1-deficient T cells failed to upregulate CTPS1, an important enzyme required for DNA synthesis. These results show that RASGRP1 deficiency leads to susceptibility to EBV infection, and demonstrate the key role of RASGRP1 at the crossroad of pathways required for the expansion of activated T lymphocytes.

  Study EGAS00001002753

• Somatic IL4R Hotspot Mutations in Primary Mediastinal Large B-cell lymphoma lead to constitutive JAK-STAT activation

  Primary mediastinal large B cell lymphoma (PMBCL) is a distinct subtype of diffuse large B cell lymphoma thought to arise from thymic medullary B cells. Gene mutations underlying the molecular pathogenesis of the disease are incompletely characterized. Here, we describe novel somatic IL4R mutations in 15 out of 62 primary cases of PMBCL (24.2%) and in all PMBCL-derived cell lines tested. The majority of mutations (11/21; 52%) were hotspot single nucleotide variants in exon 8 leading to an I242N amino acid change in the transmembrane domain. Functional analyses establish this mutation as gain-of-function leading to constitutive activation of the JAK-STAT pathway and upregulation of downstream cytokine expression profiles and B cell specific antigens. Moreover, expression of I242N mutant IL4R in a mouse xenotransplantation model conferred growth advantage in vivo. The pattern of concurrent mutations within the JAK-STAT signaling pathway suggests additive/synergistic effects of these gene mutations contributing to lymphomagenesis. Our data establish IL4R mutations as novel driver alterations and provide a strong preclinical rationale for therapeutic targeting of JAK-STAT signaling in PMBCL.

  Study EGAS00001002796

• Hominin-specific NOTCH2 paralogs expand human cortical neurogenesis through regulation of Delta/Notch interactions.

  The human cerebral cortex has undergone rapid expansion and increased complexity during recent evolution. Hominid-specific gene duplications represent a major driving force of evolution, but their impact on human brain evolution remains unclear. Using tailored RNA sequencing (RNAseq), we profiled the spatial and temporal expression of Hominid-specific duplicated (HS) genes in the human fetal cortex, leading to the identification of a repertoire of 36 HS genes displaying robust and dynamic patterns during cortical neurogenesis. Among these we focused on NOTCH2NL, previously uncharacterized HS paralogs of NOTCH2. NOTCH2NL promote the clonal expansion of human cortical progenitors by increasing self-renewal, ultimately leading to higher neuronal output. NOTCH2NL function by activating the Notch pathway, through inhibition of Delta/Notch interactions. Our study uncovers a large repertoire of recently evolved genes linking genomic evolution to human brain development, and reveals how hominin-specific NOTCH paralogs may have contributed to the expansion of the human cortex. The dataset contains raw sequences (FASTQ files) from the Illumina 2x150bp paired-end RNA sequencing profiles of 11 fetal human brain samples at 7, 9, 12, 15 and 21 gestational weeks.

  Study EGAS00001002798

• PATL2 is a key actor of oocyte maturation whose invalidation causes infertility in women and mice

  The genetic causes of oocyte meiotic deficiency (OMD), a form of primary infertility characterised by the production of immature oocytes, remain largely unexplored. Using whole exome sequencing, we found that 26% of a cohort of 23 subjects with OMD harboured the same homozygous nonsense pathogenic mutation in PATL2, a gene encoding a putative RNA-binding protein. Using Patl2 knockout mice, we confirmed that PATL2 deficiency disturbs oocyte maturation, since oocytes and zygotes exhibit morphological and developmental defects respectively. PATL2's amphibian orthologue is involved in the regulation of oocyte mRNA as a partner of CPEB. However, Patl2's expression profile throughout oocyte development in mice, alongside colocalisation experiments with Cpeb1, Msy2 and Ddx6 (three oocyte RNA-regulators) suggest an original role for Patl2 in Mammals. Accordingly, transcriptomic analysis of oocytes from WT and Patl2-/- animals demonstrated that in the absence of Patl2, expression levels of a select number of highly relevant genes involved in oocyte maturation and early embryonic development are deregulated. In conclusion, PATL2 is a novel actor of mammalian oocyte maturation whose invalidation causes OMD in humans.

  Study EGAS00001002903

• Japan PBC-GWAS Haplotype Study

  Primary biliary cholangitis (PBC) susceptibility loci have largely been discovered through single SNP association testing. In this study, we report genic haplotype patterns associated with PBC risk genome-wide in two Japanese cohorts. Among the 74 genic PBC risk haplotype candidates we detected with a novel methodological approach in a discovery cohort of 1,937 Japanese, nearly two-thirds were replicated (49 haplotypes, Bonferroni-corrected P<6.8x10-4) in an independent Japanese cohort (N=949). Along with corroborating known PBC-associated loci (TNFSF15, HLA-DRA), risk haplotypes may potentially model cis-interactions that regulate gene expression. For example, one replicated haplotype association (9q32-9q33.1, OR=1.7, P=3.0x10-21) consists of intergenic SNPs outside of the human leukocyte antigen (HLA) region that overlap regulatory histone mark peaks in liver and blood cells, and are significantly associated with TNFSF8 expression in whole blood. We also replicated a novel haplotype association involving non-HLA SNPs mapped to UMAD1 (7p21.3; OR=15.2, P=3.9x10-9) that overlap enhancer peaks in liver and memory Th cells. Our analysis demonstrates the utility of haplotype association analyses in discovering and characterizing PBC susceptibility loci.

  Study EGAS00001002915

• Clonal evolution in myeloma: the impact of maintenance lenalidomide and depth of response on the genetics and sub-clonal structure of relapsed disease in uniformly treated newly diagnosed patients

  The emergence of treatment resistant sub-clones is a key feature of relapse in multiple myeloma. Therapeutic attempts to extend remission and prevent relapse include the maximisation of response and use of maintenance therapy. We used whole exome sequencing to study the genetics of paired presentation and relapse samples from 56 newly diagnosed patients, following induction therapy, randomised to receive either lenalidomide maintenance or observation as part of the Myeloma XI trial. Patients included were considered high risk, relapsing within 30 months of maintenance randomisation. Patients achieving a complete response had predominantly branching evolutionary patterns leading to relapse, characterised by a greater mutational burden, an altered mutational profile, bi-allelic inactivation of tumour suppressor genes, and acquired structural aberrations. Conversely, in patients achieving a partial response the evolutionary features were predominantly stable with a similar mutational and structural profile. There were no significant differences between patients relapsing after maintenance lenalidomide vs observation. This study shows that the depth of response is a key determinant of the evolutionary patterns seen at relapse.

  Study EGAS00001003539

• Ewings Sarcoma RNA-Seq

  Sarcomas are cancers of the bone and soft tissue often defined by their gene fusions. However, the timing, context, and processes by which these pathogenic fusions arise are unknown. We explored this in Ewing sarcoma, a cancer driven by EWSR1-ETS fusions, with very few cooperating mutations. Combining whole-genome sequencing with enhanced informatics, we found that the EWSR1-ETS fusion arose from striking rearrangement clusters in 42% of cases (52/124). Notably, these were organized in loops that universally contained the fusion at their center, while also weaving up to 18 genes together with it. We found the same pattern of rearrangements in three additional types of sarcoma. From these data, we define a new signature for sarcoma fusions that precedes other somatic changes, in the earliest replicating DNA of the genome. This dramatic, sudden process impinges on many genes – generating multiple coding changes that profoundly affect the transcriptome, with the disease-defining gene fusion at its core. These rearrangement loops emerge in an early ES clone from which both the primary tumor and the lethal relapse emerged, and then evolved in parallel until clinically detected.

  Study EGAS00001003062

• Low and variable tumor reactivity of the intratumoral TCR repertoire in human cancers

  Infiltration of human cancers by T cells is generally interpreted as a sign of immune recognition, and there is a growing effort to reactivate dysfunctional T cells at such tumor sites. However, these efforts only have value if the intratumoral T cell receptor (TCR) repertoire of such cells is intrinsically tumor-reactive, and this has not been established in an unbiased manner for most human cancers. To address this issue, we analyzed the intrinsic tumor-reactivity of the intratumoral TCR repertoire of CD8+ T cells in ovarian and colorectal cancer – two tumor types for which T cell infiltrates form a positive prognostic marker. Data obtained demonstrate that a capacity to recognize autologous tumor is limited to approximately 10% of intratumoral CD8+ T cells. Furthermore, in two out of four patient samples tested, no tumor-reactive TCRs were identified, despite infiltration of their tumors by T cells. These data indicate that the intrinsic capacity of intratumoral T cells to recognize adjacent tumor tissue can be rare and variable, and suggest that clinical efforts to reactivate intratumoral T cells will benefit from approaches that simultaneously increase the quality of the intratumoral TCR repertoire.

  Study EGAS00001003119

• Genomic Profiling of Thyroid Cancer Reveals a Role for Thyroglobulin in Metastasis -Part II

  Papillary thyroid carcinoma (PTC) has a wide geographic variation in incidence, being highest in Saudi Arabia where it is only second to breast cancer as the most common cancer among females. Genomic profiling of PTC from Saudi Arabia has not been attempted previously. We performed whole-exome sequencing of 101 PTC samples and the corresponding genomic DNA to identify genes with recurrent somatic mutations, followed by sequencing of these genes using a next-generation gene panel approach in further 785 samples. In addition to BRAF, N-RAS and H-RAS, which have previously been shown to be recurrently mutated in PTC, our analysis highlights additional genes, including thyroglobulin (TG), which harbored somatic mutations in 3% of the entire cohort. Surprisingly, although TG mutations were not exclusive to mutations in the RAS-MAP kinase pathway, their presence was associated with a significantly worse clinical outcome, which suggests a pathogenic role beyond driving initial oncogenesis. Analysis of metastatic PTC tissue revealed significant enrichment for TG mutations (p<0.001), including events of apparent clonal expansion. Our results suggest a previously unknown role of TG somatic mutations in the pathogenesis of PTC and its malignant evolution.

  Study EGAS00001003357

• Integrated genomic characterization of adrenocortical carcinoma

  Adrenocortical carcinomas (ACC) are aggressive cancers originating in the cortex of the adrenal glands. Despite the overall poor prognosis, ACC outcome is heterogeneous. CTNNB1 and TP53 mutations are frequent in these tumors, but the complete spectrum of genetic changes remains undefined. Exome sequencing and SNP array analysis of 45 ACC revealed recurrent alterations in known drivers (CTNNB1, TP53, CDKN2A, RB1, MEN1) and genes not previously reported to be altered in ACC (ZNRF3, DAXX, TERT and MED12), which were validated in an independent cohort of 77 ACC. The cell-surface transmembrane E3 ubiquitin ligase ZNRF36 was the gene the most frequently altered (21%), and appears as a potential novel tumor suppressor gene related to the ß-catenin pathway. Our integrated genomic analyses led to the identification of two distinct molecular subgroups with opposite outcome. The C1A group of poor outcome ACC was characterized by numerous mutations and DNA methylation alterations, whereas the C1B group with good prognosis displayed a specific deregulation of two miRNA clusters. Thus, aggressive and indolent ACC correspond to two distinct molecular entities, driven by different oncogenic alterations.Array based dataset for this study has been made public at:http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE49280.

  Study EGAS00001000665

• Integrative_genome_profiling_in_AML

  AML emerges as a consequence of accumulating independent genetic aberrations that direct regulation and/or dysfunction of genes resulting in aberrant activation of signalling pathways, resistance to apoptosis and uncontrolled proliferation. Given the significant heterogeneity of AML genomes, AML patients demonstrate a highly variable response rate and poor median survival in response to current chemotherapy regimens. For the past 4 years we have conducted gene expression profiling on purified bone marrow populations equating to normal haematopoietic stem and progenitor cells from healthy subjects and patients with de novo AML in order to identify AML signatures of aberrantly expressed genes in cancer versus normal. We are now applying a series of bioinformatic methodologies combined with clinical and conventional diagnostic data to establish novel genomics strategies for improved prognostication of AML. Additionally, we use our AML signatures to unravel oncogenic signalling pathway activities in AML patients and test inhibitory drugs for these pathways inn preclinical therapeutic programmes. We consider that superimposing GEP and clinical data for our AML patient cohort with additional data on their mutational status will significantly improve the prognostic power of the study as well as unravel yet unknown mutations associated with aberrant signalling activities of oncogenic pathways.

  Study EGAS00001000858

• HCA_Thymus_Disease

  Single cell RNA sequencing on thymic cells from children with syndromic diseases such as Trisomy 21 and inflammatory/autoimmune diseases such as myasthenia gravis. Our group has previously performed atlasing of the healthy human thymus using single cell transcriptomics and obtained an unprecedented understanding of human thymus throughout development. In this study, we intend to survey diseased human thymus and compare this to healthy to obtain a better understanding of thymic dysfunction in human diseases. The diseases we will survey include congenital conditions such as Down's syndrome, CHARGE syndrome, DiGeorge syndrome as well as acquired diseases such as myasthenia gravis and thymic hyperplasia. We intend to perform genetic and genomic study on the samples, including expression analysis. This will not only help us to understand the cells and molecular pathways affected in dysfunctional thymus; but may also lead to generalisable understandings in human genetic diseases and inflammatory/autoimmune diseases. 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 EGAS00001004310

• Biased allelic expression in human primary fibroblast single cells.

  The study of gene expression in mammalian single cells using genomic technologies now provides the possibility to investigate the patterns of allelic gene expression. We have used single-cell RNA sequencing to detect the allele-specific mRNA level in 203 single human primary fibroblast cells over 133,633 unique heterozygous single nucleotide variants (hetSNVs). We have observed that at the snapshot of analyses each cell contains mostly transcripts from one allele from the majority of genes; indeed 76.4% of the hetSNVs display stochastic monoallelic expression in single cells. Remarkably, adjacent hetSNVs exhibit haplotype consistent allelic ratio; in contrast distant sites located in two different genes are independent of the haplotype structure. Moreover, the allele-specific expression in single cells correlated with the abundance of the cellular transcript. We observed that genes expressing both alleles in majority of the single cells at a given time point are rare and enriched in highly expressed genes. The relative abundance of each allele in a cell is controlled by some regulatory mechanisms since we observed related single-cell allelic profiles according to genes. Overall, these results have direct implications in cellular phenotypic variability.

  Study EGAS00001001009

• Clinical significance of novel subtypes of acute lymphoblastic leukemia in the context of minimal residual disease-directed therapy

  We evaluate clinical significance of recently identified subtypes of acute lymphoblastic leukemia (ALL) in 598 children treated with minimal residual disease (MRD)-directed therapy. Among the 16 B-ALL and 8 T-ALL subtypes identified by next generation sequencing, ETV6-RUNX1, high-hyperdiploid and DUX4-rearranged B-ALL had the best five-year event-free survival rates (95% to 98.4%); TCF3-PBX1, PAX5alt, T-cell, ETP, iAMP21, and hypodiploid ALL intermediate rates (80.0% to 88.2%); and BCR-ABL1, BCR-ABL1-like and ETV6-RUNX1-like and KMT2A-rearranged ALL the worst rates (64.1% to 76.2%). All but three of the 142 patients with day-8 blood MRD <0.01% remained in remission. Among new subtypes, intensified therapy based on day-15 MRD≥1% improved outcome of DUX4-rearranged, BCR-ABL1-like, and ZNF384-rearranged ALL, and achievement of day-42 MRD<0.01% did not preclude relapse of PAX5alt, MEF2D-rearranged and ETV6-RUNX1-like ALL. Thus, new subtypes including DUX4-rearranged, PAX5alt, BCR-ABL1-like, ETV6-RUNX1-like, MEF2D-rearranged and ZNF384-rearranged ALL have important prognostic and therapeutic implications.

  Study EGAS00001005084

• An_evaluation_of_different_strategies_for_large_scale_pooled_sequencing_study_design_

  Second generation sequencing technology has enabled the design of large-scale sequencing experiments in targeted disease-related regions in thousands of cases and controls. Cost implications dictate a pooled DNA sequencing design, evaluation of allele frequency differences based on the resequenced samples and large-scale targeted follow-up in further samples. To evaluate the feasibility of this approach and to assess different study design strategies, we have carried out a pilot study which tests the feasibility of PCR-based and pull-down DNA pooling for (a) SNP discovery and (b) allele frequency comparison purposes. All samples are from the HapMap or 1958 BC. The sample composition of the pool is listed below: Pool 5: 31 HapMap Individuals (NA12249, NA12156, NA12004, NA11831, NA12716, NA11832, NA11993, NA12057, NA11995, NA12006, NA12144, NA12802, NA12146, NA12005, NA12003, NA07000, NA12043, NA12044, NA11992, NA11881, NA11994, NA07345, NA12154, NA06994, NA06985, NA12239, NA07022, NA07034, NA12155, NA07056, NA06993); 19 1958 BC Individuals (WTCCC88214, WTCCC88215, WTCCC88216, WTCCC88217, WTCCC88222, WTCCC88233, WTCCC88240, WTCCC88241, WTCCC88242, WTCCC88247, WTCCC88249, WTCCC88262, WTCCC88264, WTCCC88278, WTCCC88294, WTCCC88298, WTCCC88302, WTCCC88305, WTCCC88321) Some pools were carried out in duplicate to assess reproducibility. We focused on 10 (PD), 7(PCR) chromosomal regions, ~2(PD), ~1.6 (PCR) Mb in total.

  Study EGAS00001000134

• Rucaparib in patients presenting a metastatic breast cancer with Homologous Recombination Deficiency, without germline BRCA1/2 mutation, pronostic value of high LOH genomic score and predictive value of HRDetect (microarray)

  Breast cancer may present genomic alterations leading to homologous recombination deficiency. PARP inhibitors have proved their efficacy in patients with HER2-negative metastatic breast cancer (mBC) harboring germline (g) BRCA1/2 mutations. We conducted the phase 2 RUBY trial to assess the efficacy of rucaparib in HER2-negative mBC with high genomic loss of heterozygosity (LOH) score or somatic, without gBRCA1/2 mutation. 220 of 711 patients with mBC screened for LOH presented high LOH score which was associated with a higher likelihood of death (HR = 1.39, 95% CI: 1.11-1.75, p = 0.005). The primary objective was not reached with a clinical benefit rate (objective response or SD>16 weeks) of 13.5%. Two LOH-high patients, without somatic BRCA1/2 mutation, presented a complete and durable response (14 and 32 months). HRDetect tended to be associated with response to rucaparib, whithout reaching statistical significance (median HRDetect responders versus non responders: 0.465 versus, 0.040, p = 0.2135). Our data suggests that a small subset of patients with high LOH score could derive benefit from PARP inhibitors.

  Study EGAS00001004755