14658 results for "PLOS+ONE+Pennsylvania+German+ACS+language+code+pdc+open+access+study"

in 46.06 milliseconds.

• Mechanisms of Therapy Driven Clonal Evolution in Oncogenic RAS Mutant Relapsed Acute Lymphoblastic Leukemia

  Although multi-agent combination chemotherapy is curative in a significant fraction of childhood acute lymphoblastic leukemia (ALL) patients, 20% of cases relapse and most die due to chemo-refractory disease. Here we used whole-exome and whole-genome sequencing to analyze the mutational landscape and pattern of clonal evolution at relapse in pediatric ALL cases. These analyses showed that ALL relapses originate from a common ancestral precursor clone of the diagnosis and relapsed populations and frequently harbor mutations implicated in chemotherapy resistance. RAS-MAPK pathway activating mutations in NRAS, KRAS and PTPN11 were present in 24/55 (44%) cases in our series. Notably, while some cases showed emergence of RAS mutant clones at relapse, in others, RAS mutant clones present at diagnosis were replaced by RAS wild type populations. Mechanistically, functional dissection of mouse and human wild type Kras and mutant Kras (Kras G12D) isogenic leukemia cells demonstrated induction of methotrexate resistance, but also improved response to vincristine, in mutant Kras- expressing lymphoblasts. These results identify chemotherapy driven selection as a central mechanism of leukemia clonal evolution and pave the road for the development of tailored personalized therapies for the treatment of relapsed ALL.

  Study phs001072

• Transcriptome of Chronic Pain and Disease

  The Functional Genomics Section (FGS) is a part of the National Institute of Dental and Craniofacial Research (NIDCR) within the National Institutes of Health (NIH). Research at the FGS focuses on understanding the mechanics involved in nociception through both clinical samples from patients with chronic pain and mouse models of pain. Pain is relayed from the periphery to the central nervous system through specialized sensory neurons known as nociceptors. The cell body of nociceptors resides in neuronal structures that branch off of the spinal cord known as dorsal root ganglia (DRG). The DRG is a heterogenous tissue containing not only sensory neurons like nociceptors, but also glial cells, capillaries, fibroblasts, and resident macrophages. The interaction between the nociceptors and these DRG cell types can have a significant role in affecting neuronal hyperexcitability and downstream pain sensitivity. To better understand the cellular factors involved in producing conditions of chronic pain, DRGs from consenting organ donors were collected from individuals experiencing pain disorders such as painful diabetic neuropathy and rheumatoid arthritis. A multi-omics approach is being used to identify molecular mechanisms that lead to pain hypersensitivity in patients with chronic pain.  

  Study phs002548

• Detection of Genes Predisposing to Hematologic Malignancies

  We have been conducting genetic studies on families at high risk of different hematologic malignancies, in order to define the related tumors in the families, define precursor and other related conditions, and map and identify susceptibility genes. We have focused mainly on four types of lymphoid malignancies: chronic lymphocytic leukemia (CLL), Hodgkin lymphoma (HL), non-Hodgkin lymphoma (NHL), and Waldenström macroglobulinemia (WM). A few families with a rare lymphoma subtype, hairy cell leukemia (HCL) are included. In addition, single large pedigrees with acute myeloid leukemia (AML), and juvenile myelomocytic leukemia (JMML) are included. Families are ascertained for having at least two patients with the same hematologic malignancy and are classified by the type of malignancy that predominates in the family. Multiple types of lymphoid malignancies are often found in the same family. Other data has shown that these conditions aggregate together in families. Verification of cancer diagnoses is obtained through medical records, pathology reports, and flow cytometry. Family members with precursor traits are also included, monoclonal B-cell lymphocytosis (MBL) in CLL families and IgM monoclonal gammopathy of undetermined significance (MGUS) in WM families.

  Study phs001219

• Epigenomes and Transcriptomes of Brain Samples of Schizophrenics and Controls

  We examined the effects of antipsychotic medications on the cell-specific epigenomics and transcriptomics in the frontal cortex of schizophrenic subjects. We performed ChIP-seq (H3K4me3 and H3K27ac) and RNA-seq on neuronal and non-neuronal nuclei of postmortem frontal cortex samples from two cohorts of schizophrenics, either previously treated or not with antipsychotic medications, and control subjects individually matched by sex and age. In total, we tested 29 schizophrenics (15 untreated and 14 treated patients) and their matched control samples. Two replicates were produced on each brain sample for each assay. ChIP-seq data were generated using MOWChIP-seq (Zhu et al., PMID: 31666743) and RNA-seq data were produced using the combination of SMART-seq2 (Picelli et al., PMID: 24385147) and a Nextera XT DNA library preparation kit. The filename contains the information on the experimental group, the number of the sample within the group, the type of the brain nuclei, the epigenomic/transcriptomic data set (either H3K27ac or mRNA), and the number of replicate. For example, “SCZ-1-untreated-neuron-H3K4me3-R1” refers untreated schizophrenic sample #1, replicate 1 ChIP-seq data on H3K4me3.

  Study phs002487

• Clonal Evolution and Transcriptomic Analysis of Chronic Lymphocytic Leukemia Treated with Ibrutinib

  Disease progression has been increasingly appreciated with the use of ibrutinib to treat chronic lymphocytic leukemia (CLL). Understanding the pattern of clonal evolution and transcriptomic changes during the early treatment period with ibrutinib may provide insight into the mechanisms of resistance. We performed whole exome sequencing on sequential samples with matched germline controls from 61 patients with CLL during the first 12 months of ibrutinib monotherapy or ibrutinib in combination with rituximab. Additionally, RNA sequencing was done on a subset of 14 patients treated with ibrutinib monotherapy. The presence of clonal shifts was associated with an increased risk of disease progression, suggesting that resistant clones are more likely to emerge in CLL with greater evolutionary capacity.This cohort was later augmented with additional data (including data from 15 additional patients). In total, this second set of data contained whole exome sequencing (31 samples), RNA sequencing (6 samples) and reduced representation bisulfite sequencing (67 samples) for CLL samples and used them to identify expression subtypes of CLL. In an integrative analysis of genetic, transcriptomic, and epigenetic data, incorporating known and newly identified subtypes of CLL, we built new models to improve patient prognostication.

  Study phs001473

• E5103 Correlative Studies

  ECOG-ACRIN E5103 (E5103) was a phase III adjuvant breast cancer trial that randomized 4994 patients with node-positive or high-risk node-negative breast cancer to intravenous doxorubicin and cyclophosphamide every 2 or 3 weeks (at the discretion of treating physician), for four cycles of adriamycin and cyclophosphamide (AC) followed by 12 weeks of weekly paclitaxel (80 mg/m2) alone (Arm A), or to the same chemotherapy with either concurrent bevacizumab (Arm B), or concurrent plus sequential bevacizumab (Arm C). To develop predictive genetic biomarkers for chemotherapy efficacy and toxicities, we have successfully performed genome-wide genotyping and whole exome sequencing in the germline (blood) DNA of 3215 and 672 patients, respectively, from the E5103 trial to investigate the association of genetic variants with chemotherapy efficacy and toxicities. Case and control studies were performed with genome-wide analyses (GWAS) and whole exome analyses (WES) across all arms of E5103 to identify genotypes at single-nucleotide polymorphisms (SNPs) and deleterious gene mutations that were associated with paclitaxel-induced peripheral neuropathy (TIPN). Through GWAS, we have identified a SNP, rs3125923, and through WES, Cytochrome P450 oxidoreductase (POR), that was significantly associated with grade 3-4 TIPN.

  Study phs003201

• The Genomic Landscape of Endocrine Resistant Advanced Breast Cancers: Paired Pre- and Post-endocrine Therapy Samples.

  The genomic evolution of breast cancers exposed to systemic therapy and its effects on clinical outcome have not been broadly characterized. We integrated the genomic sequencing of 1918 breast cancers, including 1501 hormone receptor-positive tumors, with detailed clinical information and treatment outcomes. Functional mutations in ERBB2 and loss-of-function mutations in NF1 were more than twice as common in post-endocrine therapy tumors compared to treatment-naive tumors. Additional alterations in the MAPK pathway (EGFR, KRAS among others) and in estrogen receptor transcriptional regulators (MYC, CTCF, FOXA1 and TBX3) were also more common compared to hormonal therapy-naive tumors. To determine whether candidate genomic lesions in the MAPK pathway and estrogen receptor transcriptional regulators were present prior to therapy or whether such lesions arose under the selective pressure of therapy, we performed whole-exome sequencing in select patients who had adequate samples acquired prior to and after progression on endocrine therapy. In total, 95 specimens corresponding to 30 treatment-naive primary tumors, 35 post-treatment tumors, and 30 normal samples were sequenced from 30 patients with endocrine-resistant metastatic breast cancer.

  Study phs001674

• High response rate to anti PD-1 therapy in desmoplastic melanoma

  Desmoplastic melanoma (DM) is a rare subtype of melanoma characterized by dense fibrous stroma, resistance to chemotherapy and a lack of actionable driver mutations, but is highly associated with ultraviolet light DNA damage. We analysed 60 patients with advanced DM treated with programmed cell death 1 (PD-1) or PD-1 ligand (PD-L1) blocking antibody therapy. Objective tumor responses were observed in 42 of the 60 patients (70%, 95% confidence interval 57-81%), including 19 patients (32% overall) with a complete response. Whole-exome sequencing revealed a high mutational load and frequent NF-1 mutations (14 out of 17 cases). Immunohistochemistry (IHC) analysis from 19 DM and 13 non-DM revealed a higher percentage of PD-L1 positive cells in the tumor parenchyma in DM (p = 0.04), highly associated with increased CD8 density and PD-L1 expression in the tumor invasive margin. Therefore, patients with advanced DM derive significant clinical benefit from PD-1/PD-L1 immune checkpoint blockade therapy despite being a cancer defined by its dense desmoplastic fibrous stroma. The benefit is likely derived from the high mutational burden and a frequent pre-existing adaptive immune response limited by PD-L1 expression.

  Study phs001469

• Pharmacogenomics of HIV Therapy - Atazanavir Bilirubin-related Side Effects

  This is an evaluation of genetic associations with atazanavir (ATV) discontinuation for bilirubin-related causes within 12 months of treatment. Patients were treated at an HIV primary care clinic in Nashville TN from 1998 to 2012. A previously known SNP in UGT1A1 (rs887829) was used to define metabolizer genotypes (extensive, intermediate, slow metabolizer). Over 500,000 SNPs from genome-wide genotyping were used to define MDS (Multidimensional Scaling) coordinates to account for population stratification. Patients were defined as cases if they discontinued ATV for bilirubin-related causes within 12 months of treatment, otherwise they were defined as controls if they did not stop treatment. Among 321 evaluable patients, 15 (4.6%) had bilirubin-related atazanavir discontinuation within 12 months. Homozygosity for rs887829 T/T was present in 28.1% of black, 21.4% of Hispanic, and 8.6% of white patients. Among all patients the hazard ratio (HR) for bilirubin-related discontinuation with T/T versus C/C genotype was 7.3 (95% C.I.: 1.7 to 31.5; p = 0.007). Among 152 white patients the HR was 14.4 (95% CI: 2.6 to 78.7; p = 0.002), but among 153 black patients the HR was 0.8 (95% C.I. 0.05 to 12.7; p = 0.87).

  Study phs001484

• An APOBEC Cytidine Deaminase Mutagenesis in Human Cancers

  Reprinted from Roberts et al. "An APOBEC cytidine deaminase mutagenesis pattern is widespread in human cancers", Nature Genetics, 45:970-976, 2013, with permission of Nature Publishing Group: Recent studies indicate that a subclass of APOBEC cytidine deaminases, which convert cytosine to uracil during RNA editing and retrovirus or retrotransposon restriction, may induce mutation clusters in human tumors. We show here that throughout cancer genomes APOBEC-mediated mutagenesis is pervasive and correlates with APOBEC mRNA levels. Mutation clusters in whole-genome and exome data sets conformed to the stringent criteria indicative of an APOBEC mutation pattern. Applying these criteria to 954,247 mutations in 2,680 exomes from 14 cancer types, mostly from The Cancer Genome Atlas (TCGA), showed a significant presence of the APOBEC mutation pattern in bladder, cervical, breast, head and neck, and lung cancers, reaching 68% of all mutations in some samples. Within breast cancer, the HER2-enriched subtype was clearly enriched for tumors with the APOBEC mutation pattern, suggesting that this type of mutagenesis is functionally linked with cancer development. The APOBEC mutation pattern also extended to cancer-associated genes, implying that ubiquitous APOBEC-mediated mutagenesis is carcinogenic.

  Study phs000677

• NHLBI GO-ESP: Family Studies (Familial Interstitial Pneumonia)

  The NHLBI "Grand Opportunity" Exome Sequencing Project (GO-ESP), a signature project of the NHLBI Recovery Act investment, was designed to identify genetic variants in coding regions (exons) of the human genome (the "exome") that are associated with heart, lung and blood diseases. These and related diseases that are of high impact to public health and individuals from diverse racial and ethnic groups will be studied. These data may help researchers understand the causes of disease, contributing to better ways to prevent, diagnose, and treat diseases, as well as determine whether to tailor prevention and treatments to specific populations. This could lead to more effective treatments and reduce the likelihood of side effects. GO-ESP is comprised of five collaborative components: 3 cohort consortia - HeartGO, LungGO, and WHISP - and 2 sequencing centers - BroadGO and SeattleGO. The Familial Interstitial Pneumonia (FIP) project seeks to identify genetic variants in coding regions of the human genome that are linked to FIP by examining the coding regions among relatives with FIP. These data will be used in conjunction with our other genetic studies to help us better understand how and why some individuals develop pulmonary fibrosis.

  Study phs000582

• Genome-Wide Discovery of Novel Breast Cancer Predisposing Mutations

  Cancer risk shows clear heritability, but the inherited genetic factors remain largely unknown. This project will seek to identify new genes in which mutations confer hereditary risk for early onset breast cancer, with important implications for both biological understanding and clinical prediction and prevention. Identifying cancer predisposing mutations can provide new biological insights and significantly impact important clinical decision making regarding treatment, surveillance and preventive approaches. Examples include the BRCA1 breast cancer susceptibility gene, which affects treatment decisions (surgical management), surveillance (frequent breast MRI) and prevention (oophorectory for prevention of ovarian cancer). Hereditary early onset breast cancer patients commonly present to Cancer Genetics clinics, but the majority of these patients do not have identifiable mutations in known candidate genes. Here, we will perform whole exome sequencing to discover novel cancer susceptibility genes. The patients have been recruited from Memorial Sloan Kettering Cancer Center (PIs Kenneth Offit and Zsofia Stadler) and Massachusetts General Hospital (PI Daniel Haber). All exome sequencing was performed at the Broad Institute of Harvard and MIT; samples sequence capture was performed using Agilent SureSelect Human All Exon Kit v2 and sequencing was performed on an Illumina HiSeq 2000 or 2500.

  Study phs000822

• Altered chromosomal topology drives oncogenic programs in gastrointestinal stromal tumors

  Metabolic lesions with pleiotropic effects on epigenetic regulation and other cellular processes are widely implicated in cancer, yet their oncogenic mechanisms remain poorly understood. Succinate dehydrogenase (SDH) deficiency causes a subset of gastrointestinal stromal tumors (GISTs) with DNA hyper-methylation. Here we associate this hyper-methylation with changes in chromosome topology that activate oncogenic programs. To investigate epigenetic alterations in this disease, we systematically mapped DNA methylation, CTCF insulators, enhancers and chromosome topology in KIT-mutant, PDGFRA-mutant and SDH-deficient GISTs. Although these respective subtypes share similar enhancer landscapes, we identified hundreds of putative insulators where DNA methylation replaced CTCF binding in SDH-deficient GISTs. We focused on disrupted insulators that partitions super-enhancers from FGF3, FGF4 and the KIT oncogene. Recurrent loss of this insulator alters locus topology in SDH-deficient GISTs, allowing aberrant physical interaction between enhancers and oncogenes. CRISPR-mediated excision of the corresponding CTCF motif in an SDH-intact model disrupted the boundary and up-regulated FGFs and KIT expression. Our findings reveal how a metabolic lesion destabilizes chromatin structure to facilitate the initiation and selection of epigenetic alterations that drive oncogenic programs in the absence of canonical mutations.

  Study phs001906

• Clinical Resistance to Crenolanib in Acute Myeloid Leukemia Due to Diverse Molecular Mechanisms

  FLT3 mutations are commonly detected in Acute Myeloid Leukemia (AML) patients and are associated with poor prognosis. Crenolanib, a potent type I pan-FLT3 (GeneID:2322) inhibitor, is effective against both internal tandem duplications (ITD) and resistance-conferring tyrosine kinase domain (TKD) mutations. While crenolanib monotherapy has demonstrated significant clinical benefit in heavily pretreated relapsed/refractory AML patients, responses are transient and relapse eventually occurs. To investigate the mechanisms of crenolanib resistance, we performed whole exome sequencing of AML patient samples before and after crenolanib treatment (122 samples from 59 patients). Unlike other FLT3 inhibitors, crenolanib did not induce FLT3 activation loop mutations, and mutations of the FLT3 "gatekeeper" residue were infrequent. Instead, mutations of NRAS (GeneID:4893) and IDH2 (GeneID:3418) arose, mostly as FLT3-independent subclones, while TET2 (GeneID:54790) and IDH1 (GeneID:3417) predominantly co-occurred with the FLT3-mutant clone and were enriched in crenolanib poor-responders. The remaining patients exhibited post-crenolanib expansion of mutations associated with epigenetic regulators, transcription factors, and cohesion factors, suggesting diverse non-FLT3 genetic/epigenetic mechanisms of crenolanib resistance. Drug combinations in experimental models restored crenolanib sensitivity.

  Study phs001628

• Genetic Analysis of Desmoplastic Melanoma

  Desmoplastic melanoma is an infrequent variant of melanoma with sarcomatous histology, distinct clinical behavior, and unknown pathogenesis. We performed low-coverage genome and high-coverage exome sequencing of 20 desmoplastic melanomas, followed by targeted sequencing of 293 genes to validate candidate genes. A high mutation burden (median 62 mutations/Mb) ranked desmoplastic melanoma among the most highly mutated cancers. Mutation patterns strongly implicate UV-radiation as the dominant mutagen, indicating a superficially located cell of origin. Novel alterations included recurrent promoter mutations of NF-kappa B inhibitor epsilon, NFKBIE (IkBε) in 14.5% of samples. Commonly mutated oncogenes in melanomas, in particular BRAF(V600E) and NRAS(Q61K/R), were absent. Instead, other genetic alterations known to activate the MAPK and PI3K signaling cascades were identified in 73% of samples, affecting NF1, CBL, ERBB2, MAP2K1, MAP3K1, BRAF, EGFR, PTPN11, MET, RAC1, SOS2, NRAS, and PIK3CA, some of which being candidates for targeted therapies. Reprinted from: Shain, A. H. et al. Exome sequencing of desmoplastic melanoma identifies recurrent NFKBIE promoter mutations and diverse activating mutations in the MAPK pathway. Nat. Genet. With permission from Nature Genetics

  Study phs000977

• Characterization of Autosomal CNV Among the Negrito from Peninsular Malaysia

  Copy number variation (CNV) has been recognized as a major contributor to human genome diversity. It plays an important role in determining phenotypes and has been associated with a number of common and complex diseases. However the CNV data from diverse populations is rather limited. Here we report the first investigation of copy number variation (CNV) in the indigenous populations from Peninsular Malaysia. We genotyped 34 Negrito genomes from Peninsular Malaysia using the Affymetrix SNP 6.0 microarray and identified 62 putative novel CNVs, consisting of 25 gains and 37 losses. These CNVs appear unique to the Negrito population and were absent in the DGV, HapMap3 and Singapore Genome Variation Project (SGVP) datasets. Analysis of gene ontology revealed that genes within these CNVs were enriched in the immune system (GO:0002376), response to stimulus mechanisms (GO:0050896), as well as the metabolic pathways (GO:0001852). Copy number gains in CNVRs enriched with genes were significantly higher than the losses (P value <0.001). Therefore, in view of the small population size, relative isolation and semi-normadic lifestyles of this community, we speculate that these CNVs may be attributed to recent local adaptation of Negritos from Peninsular Malaysia.

  Study phs000664

• VA APOLLO Project - Research for Precision Oncology (RePOP)

  The Research for Precision Oncology Program (RePOP) is a research activity that establishes a cohort of Veterans diagnosed with cancer and who have had genomic analyses performed on their tumor tissue as part of standard of care. All data relevant to a patient's cancer and cancer care will be collected under RePOP, including patient demographics, co-morbidities, genomic analysis, treatments, medications, lab values, imaging studies, and outcomes. All RePOP participants will have signed/verbal informed consent and signed HIPAA authorization to have their data stored and shared from RePOP's Precision Oncology Program Data Repository (PODR). The Applied Proteogenomics OrganizationaL Learning and Outcomes (APOLLO) network is a collaboration between NCI, the Department of Defense (DoD), and the Department of Veterans Affairs (VA) to incorporate proteogenomics into patient care as a way of looking beyond the genome, to the activity and expression of the proteins that the genome encodes. The emerging field of proteogenomics aims to better predict how patients will respond to therapy by screening their tumors for both genetic abnormalities and protein information, an approach that has been made possible in recent years due to advances in proteomic technology.

  Study phs001374

• Integrated Analysis of Multimodal Single-Cell Data

  The simultaneous measurement of multiple modalities, known as multimodal analysis, represents an exciting frontier for single cell genomics and necessitates new computational methods that can define cellular states based on multiple data types. Here, we introduce 'weighted nearest neighbor' analysis, an unsupervised framework to learn the relative utility of each data type in each cell, enabling an integrative analysis of multiple modalities. We apply our procedure to a CITE-seq dataset of hundreds of thousands of human white blood cells alongside a panel of 228 antibodies to construct a multimodal reference atlas of the circulating immune system. We demonstrate that integrative analysis substantially improves our ability to resolve cell states and validate the presence of previously unreported lymphoid subpopulations. Moreover, we demonstrate how to leverage this reference to rapidly map new datasets, and to interpret immune responses to vaccination and COVID-19. Our approach represents a broadly applicable strategy to analyze single-cell multimodal datasets, including paired measurements of RNA and chromatin state, and to look beyond the transcriptome towards a unified and multimodal definition of cellular identity. All CITE-seq and ECCITE-seq raw matrices are available in GEO database under the accession number GEO: GSE164378.

  Study phs002315

• Pharmacogenomic Interactions in Glioblastoma Cell Line Models

  Glioblastoma (GBM) is a lethal disease without effective treatments. To advance toward effective therapeutic approaches that are biomarker-driven, there is a need for new targeted agents and to develop accurate preclinical models that encompass its cellular and molecular diversity. Historical cell line models poorly represent molecular diversity of GBMs, so we developed a living tissue bank encompassing 103 patient-derived GBM cell lines (PDGCLs) as a platform for preclinical pharmacogenomic investigation at scale. Thirty-four PDGCLs were submitted to a comprehensive molecular characterization of copy number alterations and mutations using whole exome sequencing (WES) and gene expression (microarrays), and tested against a panel of 390 small molecules at sixteen different duplicated concentrations. After integrating and systematically correlating molecular alterations with drug sensitivities, we yielded 8,083,140 pharmacogenomic interactions, and discovered novel drug and biomarker combinations which help to prioritize highly effective novel candidate treatments in GBM. Our findings demonstrate the utility of deep characterization of large numbers of models from a single disease to detect complex predictors of compound sensitivity, and indicate novel therapeutic approaches for a large fraction of patients with GBM.

  Study phs001793

• Methotrexate Clearance GWAS in Acute Lymphoblastic Leukemia (ALL) Patients

  We estimated the plasma clearance of methotrexate from 1279 acute lymphoblastic leukemia (ALL) patients treated on the Children Oncology Group 9904 and 9905 trials (http://clinicaltrials.gov/show/NCT00005585 and http://clinicaltrials.gov/show/NCT00005596). Patients received either a 24-hour infusion of a 1 g/m2 dose or 4-hour infusion of a 2 g/m2 dose). Methotrexate clearance was lower in older children (p=7 x 10-7), females (p=2.7 x 10-4), and patients who received a delayed intensification phase (p=0.0022). In a genome-wide analysis, methotrexate clearance was associated with polymorphisms in the SLCO1B1 gene (p=2.1 x 10-11), which encodes for an organic anion transporter. This replicates findings using different schedules of high-dose methotrexate in ALL patients treated on St. Jude protocols. A combined meta-analysis yields a p-value of 5.7 x 10-19 for the association of methotrexate clearance with SLCO1B1 SNP rs4149056 (Trevino et al, PMID 19901119 and Ramsey et al, PMID 23233662). This data set includes the dependent variable of methotrexate clearance and all of the SNP data available from arrays.

  Study phs000637

• Functional Enhancer Elements Drive Subclass-Selective Expression From Mouse to Human Neocortex

  A complete epigenomic map of human brain cell types will be necessary to fully understand their genetic regulation, and to create the next generation of non-species-restricted genetic tools for their functional characterization. Here we performed both bulk and single-nucleus ATAC-seq studies of fresh human neocortical samples excised during neurosurgery, in order to discover the functional genomic elements that discriminate cell types. This dataset contains 78 bulk layer-dissected ATAC-seq samples, and 3660 single-nucleus ATAC-seq samples from fresh acutely isolated human neurosurgical samples. Furthermore, from these data we identified and validated several enhancers for use in AAV context. We also provide validation single cell RNA-seq for 456 human acutely labeled neurosurgical neurons labeled by 3 unique enhancer-AAV vectors, and single cell RNA-seq for 1483 mouse neurons labeled in vivo by 11 unique enhancer-AAV vectors. Together these datasets provide high-resolution maps of single human neocortical cell types, as well as proof of the functional utility of the genomic elements discovered. As a whole these data will be important to furnish the next generation of cell type-specific functional tools and therapies in neuroscience.  

  Study phs002292

• Resistance Development in Basal Cell Nevus Syndrome through the Basal to Squamous Transition

  While basal cell carcinomas (BCCs) arise from ectopic hedgehog pathway activation and can be treated with pathway inhibitors, sporadic BCCs display high resistance rates, while tumors arising in Gorlin syndrome patients with germline Patched (PTCH1) mutations are uniformly suppressed by inhibitor therapy. In rare cases, Gorlin syndrome patients on long-term inhibitor therapy will develop individual resistant tumor clones that rapidly progress, but the basis of this resistance remains unstudied. Here we report a case of a Smoothened inhibitor (SMOi)-resistant tumor arising in a Gorlin patient on suppressive SMOi for nearly a decade. Using a combination of multi-omics and spatial transcriptomics, we define the tumor populations at the cellular and tissue level to conclude that Gorlin tumors can develop resistance to SMOi through the previously described basal to squamous cell carcinoma transition (BST). Intriguingly, through spatial whole exome genomic analysis, we nominate PCYT2, ETNK1, and the phosphatidylethanolamine biosynthetic pathway as novel genetic suppressors of BST resistance. These observations provide a general framework for studying tumor evolution and provide important clinical insight into mechanisms of resistance to SMOi for not only Gorlin syndrome but sporadic BCCs as well.

  Study phs003415

• Molecular Determinants of Esophageal Cancer in Tanzania

  We describe the genomic profiles of 61 esophageal squamous cell carcinoma (ESCC) cases from Tanzania and compare them to profiles from an existing cohort of ESCC cases from Malawi. We also provide a comparison to ESCC tumors in The Cancer Genome Atlas Program (TCGA). Despite a high burden of ESCC in Eastern Africa, investigations into the genomics in this region are nascent. This represents the largest comprehensive genomic analysis of ESCC from sub-Saharan Africa to date. We observed substantial transcriptional overlap with other squamous histologies via comparison with TCGA PanCan dataset. DNA analysis revealed known mutational patterns, both genome-wide as well as in genes known to be commonly mutated in ESCC. TP53 mutations were the most common somatic mutation in tumors from both Tanzania and Malawi but were detected at lower frequencies than previously reported in ESCC cases from other settings. In a combined analysis, two unique transcriptional clusters were identified: a proliferative/epithelial cluster and an invasive/migrative/mesenchymal cluster. Mutational signature analysis of the Tanzanian cohort revealed common signatures associated with aging and cytidine deaminase activity and an absence of signature 29, which was previously reported in the Malawi cohort.

  Study phs003217

• Single Duplex DNA Sequencing with CODEC Detects Mutations with High Sensitivity

  Detecting mutations from single DNA molecules is crucial in many fields but challenging. Next generation sequencing (NGS) affords tremendous throughput but cannot directly sequence double-stranded DNA molecules (single duplexes) to discern the true mutations on both strands. Here, we present Concatenating Original Duplex for Error Correction (CODEC) which confers single duplex resolution to NGS. CODEC affords 1,000-fold higher accuracy than NGS, using up to 100-fold fewer reads than Duplex Sequencing. CODEC revealed mutation frequencies of 2.72 x 10-8 in sperm of a 39-year-old individual, and somatic mutations acquired with age in blood cells. CODEC detected genome-wide clonal hematopoiesis mutations from single DNA molecules, single mutated duplexes from tumor genomes and liquid biopsies, microsatellite instability (MSI) with 10-fold greater sensitivity, and mutational signatures and specific tumor mutations with up to 100-fold fewer reads. CODEC enables more precise genetic testing and reveals biologically significant mutations which are commonly obscured by NGS errors. We applied CODEC to clinical colon and breast neoplasm samples and uncovered expected mutational signatures including homologous recombination deficiency (HRD) and MSI.

  Study phs003255

• A Microwell Platform for High-Throughput Longitudinal Phenotyping and Selective Retrieval of Organoids

  In traditional bulk culture, organoids tend to overlap, hindering the analysis of individual organoid characteristics in a high-throughput manner. Additionally, variations in the local spatial properties of the bulk matrix make it challenging to distinguish whether phenotypic differences between organoids stem from inherent cellular disparities or disparities in the microenvironment. To address these challenges, we have developed a microwell-based technique that facilitates the quantification of image-based parameters for organoids grown from single cells. Furthermore, these organoids can be easily retrieved from their respective microwells for molecular profiling. By employing a deep-learning image processing pipeline, we conducted a comprehensive assessment of various phenotypic traits in two CRISPR-engineered human gastric organoid models. This analysis encompassed growth rates, cellular movement, and apical-basal polarity. Through our investigation, we successfully identified genomic alterations associated with increased growth rate, as well as changes in accessibility and expression patterns correlated with apical-basal polarity. This novel methodology offers a valuable tool for studying organoid behavior at a single-cell level, enabling researchers to gain deeper insights into organoid biology and associated genetic modifications.

  Study phs003315

• Detecting and Subtyping Lung Cancer Through Analysis of Circulating Tumor DNA

  Small cell lung cancer (SCLC) is among the most aggressive tumors with poor clinical outcomes. Tumors are rarely resected, which is a major barrier to molecular profiling of tumor tissue; therefore, non-invasive approaches using circulating tumor DNA (ctDNA) from liquid biopsies are needed to help advance research and clinical care for patients with SCLC. We developed a targeted capture panel that profiles SCLC genes and transcriptional regulation sites, including transcription start sites (TSS) and transcription factor binding sites (TFBS). This is a single assay that detects genomic mutations and transcriptional regulation activity based on ctDNA fragment patterns reflecting epigenetic nucleosomal positions. We applied this assay to sequence patient-derived xenograft (PDX) and patient plasma samples. We also performed whole genome sequencing (WGS) for a subset of the PDX plasma samples to validate the capture data. We developed new algorithms to quantify activity of key SCLC marker genes of ASCL1, NEUROD1, POU2F3, and ATOH1 and to classify SCLC subtypes based on activity of these markers. Furthermore, we developed a model to classify SCLC from non-small cell lung cancer (NSCLC), which has potential important applications for monitoring transdifferentiation and tumor plasticity.

  Study phs003570

• Ischemia Reperfusion Responses in Human Lung Transplants at the Single Cell Resolution

  Ischemia reperfusion is an unavoidable step of organ transplantation. Development of therapeutics for lung injury during transplantation has proved challenging; understanding lung injury from human data at the single cell resolution is required to accelerate the development of therapeutics. Donor lung biopsies from six human lung transplant cases were collected at the end of cold preservation and 2-hour reperfusion and underwent single cell RNA sequencing. Donor and recipient origin of cells from the reperfusion timepoint were deconvolved. Gene expression profiles were (1) compared between each donor cell type between timepoints and (2) compared between donor and recipient cells. Inflammatory responses from donor lung macrophages were found after reperfusion with upregulation of multiple cytokines and chemokines, especially IL-1 and IL-1. Significant inflammatory responses were found in alveolar epithelial cells (featured by CXCL8) and lung endothelial cells (featured by IL-6 upregulation). Different inflammatory responses were noted between donor and recipient monocytes and CD8+ T cells. The inflammatory signals and differences between donor and recipient cells observed provide insight into the cellular and molecular mechanisms of ischemia reperfusion induced lung injury. Further investigations may lead to the development of novel targeted therapeutics.

  Study EGAS50000000490

• 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

• CONNECTS Master Protocol for Clinical Trials targeting Macro- and Micro-Immuno-Thrombosis, Vascular Hyperinflammation, and Hypercoagulability and Renin-Angiotensin-Aldosterone System (RAAS) in Hospitalized Patients with COVID-19 (ACTIV-4 Host Tissue)

  ACTIV 4 Host Tissue is a shared-placebo clinical trial platform evaluating therapies targeting the host response to COVID-19 in hospitalized patients. The overarching goal of the platform was to find effective strategies for inpatient management of patients with COVID-19. Therapeutic goals for patients hospitalized for COVID-19 include hastening recovery and preventing progression to critical illness, multiorgan failure, or death. Our objective is to determine whether modulating the host tissue response improves clinical outcomes among patients with COVID-19. Three trials were implemented to investigate three agents: TXA-127, TRV-027, and Fostamatinib. These agents all impact the host tissue response in COVID-19 via a number of unique mechanisms, including potential beneficial effects on the RAAS system and formation of neutrophil extracellular traps (NETs). Each trial evaluated the efficacy of these agents' ability to impact the host tissue response and improve outcomes in patients hospitalized with COVID-19. We found no evidence that any of the three agents significantly improved clinical outcomes as measured by the primary outcome, oxygen free days through day 28.

  Study phs003708

• 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

• Somatic Mutations in Individual Skin Cells

  Ultraviolet (UV) radiation is the primary cause of melanoma through the acquisition of driver mutations in melanocytes. It also plays a significant role in oncogenesis in other skin cells, such as keratinocytes. However, the mutational profiles vary widely among different skin cancers, and the early mutational processes in different cell types from the same skin, is not well understood. Additionally, our preliminary studies also suggest that the mutation burden in different skin cells vary significantly based on risk factors such as the anatomical site and the use of tanning beds. To explore this, we employed innovative Next Generation Sequencing (NGS) tools to perform a site-matched comparison of the genomic landscapes of melanocytes, keratinocytes, and fibroblasts. From 20 biopsies of physiologically normal skin, collected from various anatomical sites across 16 donors, we have performed DNA-sequencing on 106 melanocytes, 78 keratinocytes and 13 fibroblasts. Of these samples, we have also simultaneously performed RNA-sequencing on 78 melanocytes, 46 keratinocytes and 9 fibroblasts. For normal bulk DNA, we also performed targeted (UCSF500; n=9) or exome sequencing (n=8) of skin or saliva samples from the donors.

  Study phs003683

• Yale SPORE in Skin Cancer Project 2

  Understanding the relationship between tumor and peripheral immune environments could allow longitudinal immune monitoring in cancer. Here, we examined whether T cells that share the same TCRab and are found in both tumor and blood can be interrogated to gain insight into the ongoing tumor T cell response. Paired transcriptome and TCRab repertoire of circulating and tumor-infiltrating T cells were analyzed at the single cell level from matched tumor and blood from patients with metastatic melanoma. The data shared here consist of gene expression matrix and fastq files obtained from sequencing single-cell libraries generated with the Chromium (TM) 5' technique. For each of 18 immunotherapy naive and 11 previously treated melanoma patients, T cells (Live/CD45+/TcRab+) were FACS-sorted from PBMCs and freshly resected tumors. Metadata associated with each individual cell consist of Patient I.D., biological sex, immunotherapy status, amino acid and nucleotide CDR3 sequence of the TCR alpha and beta chain (obtained from matched VDJ region enriched libraries), and sequencing batch. Further details about previous treatment history and other clinical parameters can be found in the associated manuscripts (Lucca et al., 2021: PMID: 33651881; Lu et al., PMID: 39609626).

  Study phs002289

• Blood Group Genotypes and Phenotypes in Omani Blood Donors and Its Links With Susceptibility to Malaria

  The Omani population experienced endemic malaria prior to eradication in the late 1990's. Since malaria is known to be a strong selective pressure on the human genome, particularly at blood group associated loci, we aimed to understand how selective pressures from malaria may be affecting genetic variation at blood group loci. Blood samples were collected from 100 healthy, Omani blood donors. Serology testing for 10 clinically relevant blood groups and whole genome sequencing (WGS) to an average coverage of 16X were performed on all samples. Additionally, using the WGS data, genetic variants (SNPs and INDELs) were called for chromosomes 1 - 22 using GRCh38. We find signals of selection at the Plasmodium vivax associated locus, ACKR1 (the Duffy blood group), but not at any Plasmodium falciparum associated loci. We also characterized several previously undescribed blood group alleles in the Omani population affecting antigenic expression, mostly at the RHD, RHCE, and GYPB loci. The data available from this project include BAM files from 100 Omani samples aligned to GRCh38, per chromosome VCF files, and results from the serology tests for all 10 blood groups.

  Study phs003694

• 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

• 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

• Synovial Fibroblast Gene Expression in Response to Fibronectin Fragment

  Osteoarthritis (OA) is a prevalent joint disease characterized by the progressive degradation of articular cartilage, hypertrophy of neighboring bone, and varying degrees of synovial inflammation. Recent research has highlighted the potential role of synovial fibroblasts in the pathogenesis of OA, particularly their inflammatory properties that disrupt cartilage homeostasis. Notably, degraded cartilage matrix components, such as fibronectin fragments (FN-f), have been identified in the synovial fluid of OA patients, suggesting that FN-f signaling influences the synovium where synovial fibroblasts reside. Despite affecting over 500 million individuals globally, our limited understanding of the underlying mechanisms of OA has significantly impeded the development of effective treatment strategies. To investigate this, we obtained deidentified surgical waste material from individuals undergoing knee replacement for OA. Synovial tissue was used to isolate synovial fibroblasts for cell culture studies. We generated bulk RNA-seq data from primary synovial fibroblasts (N=5 individuals) treated at passage 3 with either Phosphate-Buffered Saline (PBS) as control or FN-f for 18h to model the OA synovial fibroblast phenotype. We identified differentially expressed genes at an adjusted p-value of

  Study phs003999

• 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

• Kids First: The Genomic Architecture of Hirschsprung Disease (HSCR)

  Hirschsprung disease (HSCR) is a male biased developmental disorder associated with a lack of innervation of the gastrointestinal tract. Genetic studies have been instrumental in understanding its multifactorial inheritance, high heritability, syndromic associations, and genetic heterogeneity with variable penetrance and expressivity. At least 24 genes and 9 loci with pathogenic alleles (PAs) are known to underlie HSCR pathogenesis and explain 62% of its population attributable risk (PAR). Despite this heterogeneity, there is functional unity in HSCR: ~53% of HSCR PAs disrupt RET and EDNRB signaling in the developing enteric nervous system (ENS) in which 11 genes comprising a gene regulatory network control RET and EDNRB gene expression. We propose to identify the remaining 30% PAR by studying 857 unrelated HSCR cases, their 125 affected and 1,446 unaffected firstdegree relatives by whole genome sequencing to increase statistical power of gene discovery through improved detection of all types of coding and regulatory PAs. HSCR arises from cell autonomous defects in enteric neural crest cell precursors (ENCCs) affecting their proliferation, differentiation and migration in the ENS, functional studies that will guide our detection of novel genes.

  Study phs003662

• 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

• A Therapeutic Vaccine for Fibrolamellar Hepatocellular Carcinoma

  Fibrolamellar hepatocellular carcinoma (FLC) is a rare form of liver cancer affecting 33 children and young adults that is driven by a chimeric protein, DNAJ-PKAc. We conducted a phase 1 clinical trial of a therapeutic vaccine targeting DNAJ-PKAc 38 (FLC-Vac), in combination with nivolumab and ipilimumab, in children and adults with advanced FLC. Among 16 patients enrolled, 12 completed the vaccine priming phase and were evaluable for both immunological and clinical endpoints. Grade 3 treatment-related adverse events were reported by six patients (37.5%). DNAJ43 PKAc-specific T cell responses were detected in 9/12 patients after treatment. In the subset of patients who completed the initial priming phase the disease control rate (DCR) was 75% (9/12), with three partial responses (25%). All patients with clinical responses also had DNAJ-PKAc specific T cell responses, from which we identified multiple class II-restricted T cell receptors (TCRs) with specificity for DNAJ-PKAc. Correlates of response included both functional neoantigen reactivity and changes in TCR repertoire features over time. Immune escape in two patients corresponded with immune exhaustion rather than neoantigen escape or HLA loss.

  Study phs003970

• Establishing stable brain tumor stem cell lines and translational research for new treatments

  Choroid plexus carcinomas (CPCs) have dismal prognoses with limited clinical options. A lack of characterized cell lines has hindered the development of new treatment strategies. Here, we perform whole genome sequencing analyses of two cell lines established from CPCs. Whole genome sequencing (WGS),was performed in CCHE-45 and NGT131 cell lines. Both cell lines harbored pathogenic TP53 point mutations; CCHE-45 additionally displayed TP53 loss. Furthermore, alterations of the NOTCH and WNT pathways were also detected in both cell lines. Two protein-coding gene fusions, BZW2-URGCP, and CTTNBP2-ERBB4, mutations of two oncodrivers, GBP-4 and KRTAP-12-2, and several copy number alterations were observed in CCHE-45, but not NGT131. The discovered difference's importance and implications highlight the possible diversity of choroid plexus carcinoma and call for additional research to fully understand disease pathogenesis. Whole genome sequencing revealed that the two choroid plexus carcinoma cell lines shared TP53 mutations and other common pathway alterations and activation of NOTCH and WNT pathways. Noticeable differences were also observed. These cell lines can serve as valuable models to propose new treatments in these rare but deadly brain tumors.

  Study JGAS000657

• High-coverage whole-genome sequencing reveals structural variations in triple-negative breast cancer

  Triple-negative breast cancer (TNBC) cells do not express estrogen receptors, progesterone receptors or human epidermal growth factor receptor 2 (HER2). Currently, apart from poly ADP-ribose polymerase (PARP) inhibitors, there are few effective therapeutic options for this cancer type. Here, we present a comprehensive characterization of genetic alterations in TNBC via high-coverage whole-genome sequencing along with transcriptomic and whole-exome sequencing. Silencing of BRCA1 gene impaired the homologous recombination (HR) pathway in a subset of TNBC, which exhibited similar phenotypes with tumors with BRCA1 mutations; they harbored many structural variations (SVs) with relative enrichment for tandem duplication. Clonal analysis suggested that TP53 mutations and methylation of CpG dinucleotides in the BRCA1 promoter were the early events of carcinogenesis. SVs were associated with driver oncogenic events such as amplification of MYC, NOTCH2 or NOTCH3 and impacted tumor suppressor genes including RB1, PTEN and KMT2C. We also identified a variety of oncogenes that could transform 3T3 mouse fibroblast cells, suggesting that individual TNBC tumors may undergo a unique driver event that can be targetable. Thus, we revealed several features of TNBC that had important clinical implications.

  Study JGAS000095

• 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

• Context-specific regulatory genetic variation in MTOR dampens neutrophil-T cell crosstalk in pneumonia-associated sepsis, modulating disease

  Sepsis is a heterogeneous clinical syndrome with a high mortality, requiring personalised stratification strategies. Here, we characterise genetic variation that modulates MTOR, a critical regulator of metabolism and immune responses in sepsis. The effects are highly context specific, involving a regulatory element that affects MTOR expression in activated T cells with opposite direction of effect in neutrophils. We show that the G-allele of the lead variant, rs4845987, associated with decreased risk of T2D, reduces MTOR expression in T cells and improves survival in sepsis due to pneumonia, with effects specific to sepsis endotype. Using ex vivo models, we demonstrate that activated T cells promote immunosuppressive sepsis neutrophils through released cytokines, a process dampened by hypoxia and the mTOR inhibitor rapamycin. Our work demonstrates a novel epigenetic mechanism that fine-tunes MTOR transcription and T cell activity via the variant-containing regulatory element, which further exhibits an allelic effect upon vitamin C treatment. These findings reveal how genetic variation interacts with disease state to modulate immune cell-cell communication, providing a patient stratification strategy to inform more effective sepsis treatment.

  Study EGAS50000000894

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

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

  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 EGAS00001001788

• De novo assembly of 150 Danish genomes reveals rich structural complexity

  Most known genetic variation in human genomes has been called from comparison of short reads to the reference genome, an approach biased against finding complex variation. We sequenced 150 individuals from 50 parent-offspring trios with multiple insert-size libraries to very high coverage. We show that each genome could be independently de novo assembled into a small number of high-quality scaffolds (median N50 > 21 Mb), each of quality comparable to long read assemblies while being very cost-effective. We show that our variant call set from comparing de novo assemblies is far more complete in terms of complex variation than previous studies. Importantly, even the complex 4-5 Mb extended MHC region was assembled and resolved into haplotypes, revealing >700kb novel sequence in this important region of the genome, and major parts of the Y chromosome including some palindromes were assembled with high accuracy. Finally, we show that our variant call-set allows for the genotyping of many more complex variants when used as a reference-panel for imputation into SNP-chip data or into previously resequenced genomes.

  Study EGAS00001002108

• Genome Diversity in Africa Project: Benin (2021-02-16)

  Globally, human populations show structured genetic diversity as a result of geographical dispersion, selection and drift. Understanding this genetic variation can provide insights into the evolutionary processes that shape both human adaptation and variation in disease. Populations from SSA have the highest levels of genetic diversity. This characteristic, in addition to historical genetic admixture, can lead to complexities in the design of studies assessing the genetic determinants of disease and human variation. However, such studies of African populations are also likely to provide new opportunities to discover novel disease susceptibility loci and variants and refine gene-disease association signals. A systematic assessment of genetic diversity within SSA would facilitate genomic epidemiological studies in the region. The Genome Diversity in Africa Project (GDAP) aims to produce a comprehensive catalogue of human genetic variation in SSA, including single nucleotide polymorphisms (SNPs), structural variants, and haplotypes. This resource will make a substantial contribution to understanding patterns of genetic diversity within and among populations in SSA, as well as providing a global resource to help design, implement and interpret genomic studies in SSA populations and studies comprising globally diverse populations, complementing existing genomic resources. Specifically, we plan to carry out high depth whole genome sequencing of up to 2000 individuals across Africa (25 individuals from each ethnolinguistic group). Our scientific objectives are to: 1) develop a resource that provides a comprehensive catalogue of genetic variation in populations from SSA accessible to the global scientific community; 2) characterise population genetic diversity, structure, gene flow and admixture across SSA; 3) develop a cost-efficient, next-generation genotype array for diverse populations across SSA; and 4) facilitate whole genome-sequencing association studies of complex traits and diseases by developing a reference panel for imputation and resource for enhancing fine-mapping disease susceptibility loci. These scientific objectives will be supported by cross-cutting operational activities, including network and management of the consortium, research ethics, and research capacity building in statistical genetics and bioinformatics 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/ . This dataset contains all the data available for this study on 2021-02-16.

  Dataset EGAD00001006970

• MutWP1: CRUK Grand Challenge Mutographs of Cancer: Oesophageal adenocarcinoma

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

  Dataset EGAD00001006083

• Changes in alternative splicing and associated neo-antigens due to therapy

  Recent studies have identified alternative splicing (AS) as a novel source of neoantigens for immunotherapy. Surprisingly little is known about the AS milieu in recurrent glioblastoma (GBM), despite this being the venue for most clinical trials. We profiled 29 primary-recurrent paired human GBM specimens via RNA sequencing and re-analyzed RNA-sequencing data from non-malignant human brain tissues. From these data, we reconstructed the landscape of AS in GBM through recurrence and contrasted that to isoforms observed in non-malignant brain. The AS events we identified were cross-referenced with single-cell GBM atlases to determine cell-type specific splicing patterns. From this we identified novel splicing events in cell-surface proteins that are suitable targets for engineered T-cell therapies. We found recurrent-specific isoforms of mitogen-activated kinase pathway genes which are expressed exclusively by GBM stem-like cells that enhance invasiveness. These studies shed light on the effect of therapy on AS and identify novel targets for emerging immunotherapies.

  Study EGAS00001004524

• 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

• 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

• 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

• 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

• 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

• Germline biallelic mutation affecting the transcription factor Helios causes pleiotropic defects of immunity

  We studied a patient with recurrent respiratory infections and hypogammaglobulinemia and identified a germline homozygous missense mutation in IKZF2 encoding Helios (p.Ile325Val). We show that HeliosI325V retains DNA-binding and dimerization properties, but loses interaction with several partners, including epigenetic remodelers HDAC1, HDAC3 and the ATAC complex. Single-cell RNA-sequencing of peripheral blood mononuclear cells revealed gene expression signatures indicative of a shift towards pro-inflammatory, effector-like status in the patient’s T cells. We observed an upregulation of the anti-inflammatory gene, TSC22D3, encoding the glucocorticoid-induced leucine zipper (GILZ), whose expression was reported to increase upon IL-2 deprivation. We validated this finding in primary T cells, where we observed a pronounced reduction in IL-2 production upon their stimulation, together with a defect in their ability to proliferate, rendering them more anergic. Collectively, we identify a novel germline-encoded inborn error of immunity and define a role for Helios in conventional T cells, whereby interactions with specific binding partners is necessary to mediate the transcriptional programs that enable T-cell homeostasis in health and disease.

  Study EGAS00001005675

• 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

• Case Report: Rare IKZF1 gene fusions identified in neonate with congenital KMT2A-rearranged Acute Lymphoblastic Leukemia

  Here we describe a rare case of congenital KMT2Ar ALL presenting with co-occurring IKZF1 gene fusions and a predictably aggressive disease trajectory. We report for the first time, the novel IKZF1::TUT1 and KDM2A::IKZF1 gene fusions. Rearrangements in-volving KMT2A are commonly retained in relapsed infant ALL, however, in this case the KMT2A::AFF1 gene fusion did not appear to be the lesion driving leukemic relapse. Instead, our data suggest that relapse was driven by IKZF1::TUT1. This gene fusion remained in all samples investigated, including the on-blinatumomab therapy sample taken immediately prior to relapse. Conversely, the KMT2A::AFF1 gene fusion was only detected in the diagnosis and refractory post-induction samples highlighting a key role for IKZF1::TUT1 in disease pathogenesis. Intriguingly, both IKZF1 gene fusions are predicted to be out-of-frame, however, our data demonstrate the IKZF1 gene is still expressed. This is not unprecedented and it has previously been observed that out-of-frame fusions can cause transcriptional activation/repression of genes involved in the fusions leading to increases or decreases of their expression and the as-sociated functional outcomes.

  Study EGAS00001006947

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• Clonal relationships of ductal carcinoma in situ and recurrent invasive breast cancers defined by genomic analysis

  Ductal Carcinoma In Situ (DCIS) is the most common form of pre-invasive breast cancer and despite treatment a small fraction (5-10%) of DCIS patients present with invasive disease many years later. A fundamental biologic question is whether the invasive disease recurring in the same breast is established by tumor cells in the initial DCIS or represents new unrelated disease. To address this question, we performed genomic analyses on the initial pure DCIS lesion and paired invasive recurrent tumors in 95 patients together with single cell DNA sequencing in a subset of cases. Our data shows that in 75% the invasive recurrence was clonally related to the initial DCIS, suggesting that the tumor cells were not eliminated during the initial treatment with surgery +/- radiotherapy. Surprisingly however, 18% were clonally unrelated to the DCIS, representing new independent lineages, and 7% of cases were ambiguous. Our findings show that although DCIS is often the precursor of invasive recurrence, a significant fraction of invasive recurrences are unrelated to the initial DCIS. This knowledge is essential for accurate risk evaluation of DCIS treatment de-escalation strategies and the identification of predictive biomarkers.

  Study EGAS00001005784

• 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

• Single cell RNA sequencing of synovial B cells in early Rheumatoid Arthritis

  Rheumatoid Arthritis (RA) is a prevalent autoimmune disease characterized by inflammation of the peripheral joints and occurence of autoantibodies, i.e. Rheumatoid Factor (RF) and anti-citrullinated protein antibodies (ACPA), in circulation. Previous studies have shown that there are memory B cells and autoantibody producing plasma cells present in the joint tissue of patients with long-standing RA. However, it has remained unclear, whether these B cell populations are already present in the joint tissue of patients at the disease onset and whether these are autoreactive. Here, we used a single cell RNA sequencing approach to dissect the B cell repertoire at this early timepoint. We found evidence for B and T cell interaction and presence of memory and plasma cell pools in ACPA- and ACPA+ RA. Our results demonstrated clonal relationships between the memory and plasma cell compartments and autoreactivity within the plasma cell pool. These findings challenge our understanding of the local adaptive immune response in the joints of ACPA- and ACPA+ RA patients with direct implications for B and T cell targeting therapy in both patient subgroups.

  Study EGAS00001005144

• UAMS Smoldering Myeloma Myeloma Sequencing

  Identifying high risk smoldering myeloma patients and progression mechanism is a prerequisite to implement effective inception strategies and curve myeloma related morbidity and mortality. We hypothesize that genomics may help identify determinants of progression that may help predict outcome and offer effective chemo preventive targets. Eighty-two patients underwent a custom targeted panel sequencing with an additional capture for the translocation loci. These results were compared to 223 newly diagnosed patients (EGAS00001003223 and EGAD00001004117) and 17 MGUS and 10 early myeloma patients. DNA was obtained from either CD138+ cells from the bone marrow of multiple myeloma patients (tumor) or from stem cell harvests or peripheral blood cells from the same patient (control). 100 ng of DNA was fragmented, end-repaired, and adapters ligated using the HyperPlus kit (KAPA Biosystems). After PCR amplification the libraries were hybridized with probes against either a targeted panel consisting of 140 genes and chromosomal regions (Nimblegen) using SeqCap reagents (Nimblegen). Hybridized libraries underwent further amplification before being sequenced on a NextSeq500 (Illumina) using 75 bp paired end reads. Overall, these data highlight the importance of dysregulation of the MAPK pathway in the progression to MM.

  Study EGAS00001003629

• Identification of rare germline variants in familial multiple myeloma

  The risk of developing multiple myeloma (MM) in first-degree relatives of affected individuals is 2 to 4 fold higher than in the general population, implying the existence of inherited susceptibility. While over 20 common, low-penetrance variants have been associated with MM risk, less is known about rare, high-penetrance germline predisposition. We aimed to identify and characterize high/moderated penetrance germline variants in MM families using whole genome sequencing. For that purpose we identified families with at least two individuals diagnosed with MM or its precursors MGUS and smoldering MM. Also, patients with solitary plasmacytoma and AL amyloidosis were enrolled. We also recruited unaffected family members, and they were screened to exclude undetected MM or its precursor stages. After whole genome sequencing of blood DNA and filtering for rare variants (minor allele frequency <0.1%) and family segregation, variants are prioritized by a range of published in silico tools and literature search. The most promising variants are subjected to experimental validation to evaluate their functional relevance. Our data will provide guidance to worldwide efforts to unravel germline predisposition in MM and validation of the proposed variants and genes.

  Study EGAS00001004734

• Familial psychosis associated with a missense mutation at MACF1 gene combined with the rare duplications DUP3p26.3 and DUP16p23.3, affecting the CNTN6 and CDH13 genes

  Psychosis is a highly heritable and heterogeneous psychiatric condition. Its genetic architecture is thought to be the result of the joint effect of common and rare variants. Families with high prevalence are an interesting approach to shed light on the rare variant's contribution without the need of collecting large cohorts. To unravel the genomic architecture of a family enriched for psychosis, with four affected individuals, we applied a system genomic approach based on karyotyping, genotyping by whole-exome sequencing to search for rare single nucleotide polymorphisms (SNPs) and SNP array to search for copy-number variants (CNVs). We identified a rare non-synonymous variant, g.39914279 C>G, in the MACF1 gene, segregating with psychosis. Rare variants in the MACF1 gene have been previously detected in SCZ patients. Besides, two rare CNVs, DUP3p26.3 and DUP16p23.3, were also identified in the family affecting relevant genes (CNTN6 and CDH13, respectively). We hypothesize that the co-segregation of these duplications with the rare variant g.39914279 C>G of MACF1 gene precipitated with schizophrenia and schizoaffective disorder.

  Study EGAS00001004791

• Risk and modifying factors in Frontotemporal Dementia

  Understanding the molecular mechanisms underlying frontotemporal dementia (FTD) is essential for the development of successful therapies. Here, we present Phase 1 of a multi-omics, multi-model data resource for FTD research which will allows in-depth molecular research into these mechanisms. We have integrated and analysed data from the frontal lobe of FTD patients with mutations in MAPT, GRN and C9orf72 and detected common and distinct dysregulated cellular pathways. Our results highlight that excitatory neurons are the most vulnerable neuronal cell type and that vascular aberrations are a common hallmark in FTD. Via integration of multi-omics data, we detected several transcription factors and pathways which regulate the strong neuroinflammation observed in FTD-GRN. Finally, using small RNA-seq data and verification experiments in cellular models, we identified several up-regulated miRNAs that inhibit cellular trafficking pathways in FTD and lead to microglial activation. In this work we shed light on novel mechanistic and pathophysiological hallmarks of FTD. In addition, we believe that this comprehensive, multi-omics data resource will further mechanistic FTD research by the community.

  Study EGAS00001004895

• Epigenome wide DNA methylation assay of gingivo-buccal oral squamous cell carcinoma using single base resolution high throughput array

  Gingivo-buccal oral squamous cell carcinoma (OSCC-GB) is the most common cancer among men in India and is associated with high mortality. Although OSCC-GB is known to be quite different from tongue cancer in its genomic presentation and its clinical behavior, it is treated identically as tongue cancer. Predictive markers of prognosis and therapy that are specific to OSCC-GB are, therefore, required. To achieve this, we have carried out epigenomic (n=87) profiling of paired tumor-normal tissues collected from OSCC-GB patients from India and genome-wide DNA methylation assays were performed. DNA was isolated using DNeasy Blood and Tissue Kit (Qiagen). The purity and concentration was estimated using NanoDrop 2000 (Thermo Fisher Scientific). Approximately 500 ng genomic DNA from each sample was used for sodium bisulfite conversion using the EZ DNA methylation Gold Kit (Zymo Research, USA). Genome-wide DNA methylation was assayed using iScan (Illumina), for paired tumour and adjacent normal samples of 25 patients using the Infinium Human Methylation450 BeadChip and of 62 patients using the Infinium MethylationEPIC BeadChip.

  Study EGAS00001003896

• The molecular basis of T-PLL is an actionable perturbation of TCL1/ATM- and epigenetically instructed damage responses

  T-cell prolymphocytic leukemia (T-PLL) is a rare and poor-prognostic mature T-cell malignancy. To address its incomplete molecular concept, we integrated large-scale profiling data of alterations in gene expression, allelic copy number (CN), and nucleotide sequences in 111 well-characterized patients. Besides prominent signatures of T-cell activation and prevalent clonal variants, we also identified novel hot-spots for CN variability, fusion molecules, alternative transcripts, and progression-associated dynamics. The overall lesional spectrum of T-PLL is mainly annotated to axes of DNA damage responses, T-cell receptor/cytokine signaling, and histone modulation. We formulate a multi-dimensional model of T-PLL pathogenesis centered around a unique combination of TCL1 overexpression with damaging ATM aberrations as initiating core lesions. The effects imposed by TCL1 cooperate with compromised ATM towards a leukemogenic phenotype of impaired DNA damage processing. Dysfunctional ATM appears inefficient in alleviating elevated redox burdens and telomere attrition and in evoking a p53-dependent apoptotic response to genotoxic insults. As non-genotoxic strategies, synergistic combinations of p53 reactivators and deacetylase inhibitors reinstate such cell death execution.

  Study EGAS00001002744