22928 results for "ega"

in 87.16 milliseconds.

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Colorectal_LCM

  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. Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

  Study EGAS00001003455

• Exome sequencing of a novel cervical cancer cell line

  Human cancer cell lines are largely used in the searching for new antineoplastic agents. However, due to the artifacts of a long-term in culture, cell lines do not always represent the realistic tumor cell behavior. This has motivated the development of models that better mimetics the tumor tissue, among them, the establishment of primary cell cultures. In this work, we establish and characterized a low-passage cervix cancer cell line from a Brazilian patient with squamous cell carcinoma. The phenotype confirms the epithelial and tumor origin, through cytokeratins, EpCAM, and p16 staining. Whole exome sequencing showed relevant somatic mutations in several genes including BRCA2, TGFBR1, and IRX2 genes. CNV analysis by nanostring and WGS revealed amplification in genes related mainly with kinases proteins, involved in proliferation, migration and cell differentiation, such as EGFR, PIK3CA, and MAPK7. Overexpression of EGFR was confirmed by phospho RTK-array and western blot analysis. Furthermore, the cell was sensitivity to cisplatin, with IC50 13 times lower than SiHa cell line. In conclusion, this cervical cancer cell line presents molecular alterations that are an important tool for leading pre-clinical studies of new drugs that target one or more of the altered pathways.

  Study EGAS00001003343

• The evolutionary trajectory of a highly recurrent paediatric high grade neuroepithelial tumour with MN1:BEND2 fusion

  Astroblastomas are rare brain tumours which predominate in children and young adults, and have a controversial claim as a distinct entity, with no established WHO grade. Reports suggest a better outcome than high grade gliomas, though they frequently recur. Recently, they have been described to overlap with a newly-discovered group of tumours described as ’high grade neuroepithelial tumour with MN1 alteration’ (CNS HGNET-MN1), defined by global methylation patterns and strongly associated with gene fusions targeting MN1. We have studied a unique case of astroblastoma arising in a 6 year old girl, with multiple recurrences over a period of 10 years, with the pathognomonic MN1:BEND2 fusion. Exome sequencing allowed for a phylogenetic reconstruction of tumour evolution, which when integrated with clinical, pathological and radiological data provide for a thorough understanding of disease progression, with initial treatment driving tumour dissemination along four distinct trajectories. Infiltration of distant sites was associated with a later genome doubling, whilst there was evidence of convergent evolution of different lesions acquiring distinct alterations targeting NF-κB. These data represent an usual opportunity to understand the evolutionary history of a highly recurrent childhood brain tumour, and provide novel therapeutic targets for astroblastoma / CNS HGNET-MN1.

  Study EGAS00001002432

• Genomic profiling of Acute Lymphoblastic Leukemia in Ataxia Telangiectasia patients reveals tight link between ATM mutations and chromothripsis

  Recent developments in sequencing technologies led to the discovery of a novel form of genomic instability, termed chromothripsis. This catastrophic genomic event, involved in tumorigenesis, is characterized by tens to hundreds of simultaneously acquired locally clustered rearrangements on one chromosome. We hypothesized that leukemias developing in individuals with Ataxia Telangiectasia, who are born with two mutated copies of the ATM gene, an essential guardian of genome stability, would show a higher prevalence of chromothripsis due to the associated defect in DNA double-strand break repair. Using whole-genome sequencing, fluorescence in situ hybridization and RNA sequencing, we characterized the genomic landscape of Acute Lymphoblastic Leukemia (ALL) arising in patients with Ataxia Telangiectasia. We detected a high frequency of chromothriptic events in these tumors, specifically on acrocentric chromosomes, as compared to tumors from individuals with other types of DNA repair syndromes (27 cases total, 10 with Ataxia Telangiectasia). Our data suggest that the genomic landscape of Ataxia Telangiectasia ALL is clearly distinct from that of sporadic ALL. Mechanistically, short telomeres and compromised DNA damage response in cells of Ataxia Telangiectasia patients may be linked with frequent chromothripsis. Furthermore, we show that ATM loss is associated with increased chromothripsis prevalence in additional tumor entities.

  Study EGAS00001002270

• Whole-Genome Sequencing Suggests Schizophrenia Risk Mechanisms in Humans with 22q11.2 Deletion Syndrome

  Chromosome 22q11.2 microdeletions impart a high but incomplete risk for schizophrenia. Possible mechanisms include genome-wide effects of DGCR8 haploinsufficiency. In a proof-of-principle study to assess the power of this model, we used high-quality, whole-genome sequencing of nine individuals with 22q11.2 deletions and extreme phenotypes (schizophrenia, or no psychotic disorder at age >50 years). The schizophrenia group had a greater burden of rare, damaging variants impacting protein-coding neurofunctional genes, including genes involved in neuron projection (nominal P = 0.02, joint burden of three variant types). Variants in the intact 22q11.2 region were not major contributors. Restricting to genes affected by a DGCR8 mechanism tended to amplify between-group differences. Damaging variants in highly conserved long intergenic noncoding RNA genes also were enriched in the schizophrenia group (nominal P = 0.04). The findings support the 22q11.2 deletion model as a threshold-lowering first hit for schizophrenia risk. If applied to a larger and thus better-powered cohort, this appears to be a promising approach to identify genome-wide rare variants in coding and noncoding sequence that perturb gene networks relevant to idiopathic schizophrenia. Similarly designed studies exploiting genetic models may prove useful to help delineate the genetic architecture of other complex phenotypes.

  Study EGAS00001002344

• Natural genetic variation of the cardiac transcriptome in non-diseased donors and patients with dilated cardiomyopathy

  Background: Genetic variation is an important determinant of RNA transcription and splicing, which in turn contributes to variation in human traits including cardiovascular diseases.Results: Here we report the first in-depth survey of heart transcriptome variation using RNA-sequencing in 149 (97)* patients with dilated cardiomyopathy and 113 (108)* non-diseased controls. We reveal extensive differences of gene expression and splicing between dilated cardiomyopathy patients and controls, affecting known as well as novel dilated cardiomyopathy genes. Moreover, we show a widespread effect of genetic variation on the regulation of transcription, isoform usage and allele specific expression. Systematic annotation of genome wide association SNPs identifies 60 functional candidate genes for heart phenotypes, representing 20% of all published heart genome wide association loci. Focusing on the dilated cardiomyopathy phenotype we found that eQTL variants are also enriched for dilated cardiomyopathy genome wide association signals in two independent cohorts.Conclusions: RNA transcription, splicing and allele specific expression are each important determinants of the DCM phenotype and are controlled by genetic factors. Our results represent a powerful resource for the field of cardiovascular genetics. (*) numbers in parentheses represent samples that the publication is based on.

  Study EGAS00001002454

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Colorectal

  he 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. Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

  Study EGAS00001003774

• The evolutionary steps from primary to metastatic prostate cancer are largely uncharted, and the ability to use DNA present in body fluids as correlates of aggregate metastatic status is under-examined. We reconstructed phylogenies in ten prostate cancer patients with fatal disease using deep targeted sequencing of the prostate, adjacent and distant organs, as well as plasma, serum, and cerebrospinal fluid at various time points. A total of 163 samples are studied.

  The evolutionary steps from primary tumor to metastasis in prostate cancer are largely uncharted, and the ability to use serum, plasma, and cerebrospinal fluid as a correlate of aggregate metastatic tumor genomic status has not been tested. We used deep targeted sequencing to reconstruct tumor evolution in ten prostate cancer patients with fatal disease encompassing examination of the prostate and adjacent and distant organs, as well as plasma, serum, and cerebrospinal fluid at various time points. We show that there is substantial evolution from a common ancestor within the prostate that results in branching to multiple lineages which form an intermixed multi-clonal primary tumor mass. After the occurrence of key driver aberrations, one of these lineages will metastasize to multiple sites in a sequential fashion. These metastatic sites are then susceptible to being populated by cells from other intra-prostatic lineages or from other metastases. Genomic representation of metastases in body fluids is not uniform. Cerebrospinal fluid analysis can detect lineages not detected in circulating DNA, suggesting possible clinical utility.

  Study EGAS00001003848

• CancerDetector: Ultrasensitive and Non-Invasive Cancer Detection at the Resolution of Individual Reads using Cell-free DNA Methylation Sequencing Data

  The detection of tumor-derived cell-free DNA in plasma is one of the most promising directions in cancer diagnosis. The major challenge in such approach is how to identify the tiny amount of tumor DNAs out of total cell-free DNAs in blood. Here we propose an ultrasensitive cancer detection method, termed “CancerDetector”, using the DNA methylation profiles of cell-free DNAs. The key of our method is to probabilistically model the joint methylation patterns of multiple adjacent CpG sites on an individual sequencing read, in order to exploit the pervasive nature of DNA methylation for signal amplification. Therefore, CancerDetector can sensitively identify a trace amount of tumor cfDNAs in plasma, at the level of individual reads. We evaluated CancerDetector on the simulated data, and showed a high concordance of the predicted and true tumor burden. Testing CancerDetector on real plasma data demonstrated its high sensitivity and specificity in detecting tumor DNAs. In addition, the predicted tumor burden showed great consistency with tumor size and survival outcome. Note that all of those testing were performed on sequencing data at low to medium coverage (1X to 10X). Therefore, CancerDetector holds the great potential to detect cancer early and cost-effectively.

  Study EGAS00001002728

• WGS of 32 paired SRCC samples

  Signet-ring cell carcinoma (SRCC) has specific epidemiology and oncogenesis in gastric cancer, however with no systematical investigation for prognostic genomic features. We systematically investigated the clinical characteristics of 1,868 Chinese gastric cancer patients, and found that signet-ring cell content was significantly related to multiple clinical characteristics and treatment outcomes. Next, we performed whole genome sequencing for 32 SRCC patients with > 80% signet-ring cells in tumor, and identified frequent5 CLDN18-ARHGAP26/6 fusion (25%). With additional 797 patients for validation, we notice that prevalence of CLDN18-ARHGAP26/6 fusion is positively associated with signet-ring cell content (P = 4.1 x 10 -9 ), younger age at diagnosis (P = 4.2 x 10 -10 ), female/male ratio (P = 1.7 x 10 -9 ), and advanced TNM stage (P = 1.7 x 10 -5 ). Importantly, patients with CLDN18-ARHGAP26/6 fusion had a worse survival outcomes (P = 0.03), and got no benefit0 from platinum/fluoropyrimidines based chemotherapy (P = 0.92) compared to fusion-free patients (P = 0.001), which is consistent with the fact of oxaliplatin/fluoropyrimidine resistance acquired in CLDN18-ARHGAP26 introduced gastric cancer cell lines. Overall, this study provides new insights into the clinical and genomic features of SRCC, and highlights the importance of frequent CLDN18-ARHGAP26/6 fusions in chemotherapy response for SRCC.

  Study EGAS00001002668

• METABRIC: Data from Batra et al (2021); DNA methylation landscapes of 1538 breast cancers reveal a replication-linked clock, epigenomic instability and cis-regulation

  DNA methylation is aberrant in cancer, but the dynamics, regulatory role and clinical implications of such epigenetic changes are still poorly understood. Here, RRBS-based profiles of 1538 breast tumors and 244 normal breast tissues from the METABRIC cohort are reported, facilitating detailed analysis of DNA methylation within a rich context of genomic, transcriptional, and clinical data. Tumor methylation from immune and stromal signatures are deconvoluted leading to the discovery of a tumor replication-linked clock with genome-wide methylation loss in non-CpG island sites. Unexpectedly, methylation in most tumor CpG islands follows two replication-independent processes of gain (MG) or loss (ML) that we term epigenomic instability. Epigenomic instability is correlated with tumor grade and stage, TP53 mutations and poorer prognosis. After controlling for these global trans-acting trends, as well as for X-linked dosage compensation effects, cis-specific methylation and expression correlations are uncovered at hundreds of promoters and over a thousand distal elements. Some of these targeted known tumor suppressors and oncogenes. In conclusion, the study demonstrates that global epigenetic instability can erode cancer methylomes and expose them to localized methylation aberrations in-cis resulting in transcriptional changes seen in tumors.

  Study EGAS00001004327

• Characterizing microbiome-directed fibre snacks in gnotobiotic mice and humans

  Knowledge of the interrelationships between what we eat and the configurations of our gut microbial communities is providing important insights into how food components that are not directly metabolized by human enzymes are linked to our physiology and health status. Changing food preferences brought about by Westernization that have deleterious health effects1,2, plus rapid population expansion, ongoing challenges to sustainable agriculture, and other forces contributing to increased food insecurity, are catalyzing efforts to identify more nutritious and affordable foods3. The gut microbial community is complex, dynamic, and exhibits considerable intra- and interpersonal variation in its composition and functions. The massive number of potential interactions between its components makes it challenging to define the mechanisms by which food ingredients affect community properties. There is also a paucity of information about the ‘bioactive’ ingredients of foods that influence the fitness and expressed functions of community members. Here, plant fibres, from different sustainable sources and targeting distinct features of obese human gut microbiomes in gnotobiotic mice, were formulated into snack prototypes and used to supplement controlled diets consumed by overweight and obese adults; the results revealed fibre-specific changes in their microbiomes that were linked to changes in their plasma proteomes indicative of altered physiologic state.

  Study EGAS00001005268

• Molecular characteristics in Burkitt lymphoma over age groups

  Treatment optimization has improved the overall survival in Burkitt lymphoma (BL) patients but there are still relevant outcome differences between pediatric and adult patients. In both age groups survival remains poor for patients who suffer relapse. Therefore, new treatment approaches are warranted to improve survival. Thus, we conducted a comparative study in the largest pediatric cohort to date including 191 cases and 97 samples from adult patients and analyzed the mutational landscape. We detected overlapping molecular profiles but also significant differences in the mutational pattern. Samples from pediatric patients show a higher frequency of mutations in ID3, DDX3X and BAF-complex components ARID1A and SMARCA4, while YY1AP1 and BCL2 are more frequently mutated in cases from adult BL patients. Correlation of the mutational profile with clinical data was available for the pediatric cohort and identified TP53 mutations to be highly significantly associated with the incidence of relapse thus identifying the first molecular marker to be associated with the incidence of relapse in pediatric BL. As the hotspot analysis of TP53 mutations showed mostly gain of function mutations, it is tempting to speculate that drugs decreasing the amount of malignant TP53 or increase the abundance of wildtype TP53 could be a promising approach for future treatment regimens to further improve prognosis.

  Study EGAS00001005270

• Whole genome, whole exome, and targeted sequencing of high-grade meningioma tumor samples.

  High-grade meningiomas frequently recur and are associated with high rates of morbidity and mortality. To determine the factors that promote the development and evolution of these tumors, we analyzed the genomes of 274 high-grade meningiomas and compared this information with data from 456 low-grade meningiomas. High-grade meningiomas had a higher mutation burden than low-grade meningiomas but did not harbor any statistically significant mutated genes aside from NF2. High-grade meningiomas also possessed significantly elevated rates of chromosomal gains and losses, especially among tumors with monosomy 22. Meningiomas previously treated with adjuvant radiation had significantly more copy number alterations than radiation-induced or radiation-naïve meningiomas. Across serial recurrences, genomic disruption preceded the emergence of nearly all mutations, remained largely uniform across time, and when present in low-grade meningiomas, correlated with subsequent progression to a higher grade. In contrast to the largely stable copy number alterations, mutations were strikingly heterogeneous across tumor recurrences, likely due to extensive geographic heterogeneity in the primary tumor. While high-grade meningiomas harbored significantly fewer overtly targetable alterations than low-grade meningiomas, they contained numerous mutations that are predicted to be neoantigens, suggesting that immunologic targeting may be of therapeutic value.

  Study EGAS00001002294

• Exome_trios_in_patients_with_gastroschisis

  Gastroschisis (MIM 230750) is a herniation of the intestines through a defect of the abdominal wall lateral to the umbilicus (usually on the right side), and it is not covered by a membrane [Ledbetter, 2012]. Gastroschisis is a congenital anomaly with increasing incidence, easy prenatal diagnosis and extremely variable postnatal outcomes. On the basis of clinical manifestations, epidemiologic charateristics, and the presence and type of additional malformations, gastroschisis could be considered a heterogeneous condition with no gene/s discovered yet. This congenital anomaly affects approximately 1-3 infancts per 10,000 live births [Calzolari et al.1995;Parker et al.,2010] Current knowledge about causative mutations/variants. To date, no single gene has been linked to gastroschisis. Some publications have tried to link this malformation to variants in genes (such as AEBP1 (adipocyte enhancer binding protein) gene [Feldkamp et al,. 2012] or the VEGF-NOS3 pathway [Lammer et al., 2008]. Previously, a Scribble mutant mouse model (circletail) was reported to exhibit gastroschisis, however recent studies demonstrated that the Scribble knockout fetus exhibits exomphalos phenotype of gastroschisis [Carnagham et al., 2013]. 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 EGAS00001002664

• Clonal_haematopoiesis_in_patients_with_AAA

  Recent advances in genomics have demonstrated that clonal haemopoiesis driven by leukaemia associated somatic mutations is a relatively common phenomenon that increases in frequency with advancing age. Whilst individuals with clonal haemopoiesis have an increased risk of developing haematological malignancies, they also have an increased mortality from other causes. Additionally, certain mutations are almost exclusively seen in individuals aged 70 years or older, whilst others are seen in individuals with non-haematological cancers including breast and ovarian. Recently, clonal haemopoiesis was found to be associated with a significantly increased risk of atherosclerotic cardiovascular disease. This association is thought to be causative with clonally-derived macrophages showing elevated expression of several chemokine and cytokine genes that contribute to atherosclerosis. Another vascular pathology, abdominal aortic aneurysm (AAA), increases with age and shares risk factors with atherosclerosis (including smoking, male sex, high cholesterol). However, the impact of these risk factors and the overlap between AAA and atherosclerosis is poorly understood. To investigate a possible link between clonal haemopoiesis and AAA, we will study DNA samples from 300 patients with AAA and up to 200 controls for evidence of clonal haemopoiesis. This will be done using target DNA enrichment with biotinylated RNA baits followed by high throughput sequencing.

  Study EGAS00001002873

• IDENTIFICATION AND TARGETED MANAGEMENT OF A PATIENT WITH A NEURODEGENERATIVE DISORDER CAUSED BY BIALLELIC MUTATIONS IN SLC5A6

  We describe a sibling pair displaying an early infantile-onset, progressive neurodegenerative phenotype, with symptoms of developmental delay and epileptic encephalopathy developing from 12 to 14 months of age. Using whole exome sequencing, compound heterozygous variants were identified in SLC5A6, which encodes the sodium-dependent multivitamin transporter (SMVT) protein. SMVT is an important transporter of the B-group vitamins biotin, pantothenate and lipoate. The protein is ubiquitously expressed and has major roles in vitamin uptake in the digestive system, as well as transport of these vitamins across the blood-brain barrier. Pathogenicity of the identified variants was demonstrated by impaired biotin uptake of mutant SMVT. Identification of this vitamin transporter as the genetic basis of this disorder guided targeted therapeutic intervention, resulting clinically in improvement of the patient’s neurocognitive and neuromotor function. This is the second report of biallelic mutations in SLC5A6 leading to a neurodegenerative disorder due to impaired biotin, pantothenate and lipoate uptake. The genetic and phenotypic overlap of these cases confirms mutations in SLC5A6 as the genetic cause of this disease phenotype. Recognition of the genetic disorder caused by SLC5A6 mutations is essential for early diagnosis and to facilitate timely intervention by triple vitamin (biotin, pantothenate and lipoate) replacement therapy.

  Study EGAS00001003861

• Rapid identification of somatic genome rearrangements as personalized biomarkers for blood-based cancer monitoring

  Increasing evidence shows the value of circulating tumour DNA (ctDNA) to detect cancer and monitor its progression. Somatic genomic structural variations (SVs) are promising personalized biomarkers for sensitive and specific detection of ctDNA in liquid biopsies. However, accurate, affordable, and fast identification of such SV biomarkers is challenging, which hinders routine use in the clinic. Here, we demonstrate a novel approach - termed SHARC - for rapid discovery of somatic SV breakpoints as personalized tumour biomarkers. SHARC combines low-coverage cancer genome sketching by using Oxford Nanopore portable sequencing with a random forest classification and a dedicated filtering pipeline to enrich for somatic SVs. Our method leverages the real-time and long-read capabilities of Nanopore sequencing to identify somatic SV breakpoints at nucleotide resolution from a tumour biopsy within two days. We applied SHARC to tumour samples of high-grade ovarian and prostate cancer and validated on average 10 somatic SVs per sample with the use of PCR mini-amplicons. Finally, we demonstrate that these somatic SV biomarkers can be used to detect tumour presence from liquid biopsies in a quantitative manner and we retrospectively monitored treatment response in patients with prostate cancer, demonstrating its potential benefit for clinical practice.

  Study EGAS00001003963

• Mutational Patterns in Metastatic Cutaneous Squamous Cell Carcinoma

  Cutaneous squamous cell carcinoma (cSCC) from the head and neck typically metastasizes to the lymph nodes of the neck and parotid glands. When a primary is not identified, they are difficult to distinguish from metastases of mucosal origin and primary salivary gland SCC. Ultraviolet radiation causes a mutation pattern that predominantly features cytosine to thymine transitions at dipyrimidine sites and has been associated with cSCC. In this study, we used whole genome sequencing data from 15 cSCC metastases and show that a UV signature mutation is pervasive across the cohort and distinct from mucosal SCC. The mutational burden was exceptionally high and concentrated in some regions of the genome, especially insulator elements (mean 162 mutations / Mb). We therefore evaluated the likely impact of UV-induced mutations on the dipyrimidine rich binding site of the main human insulator protein, CCCTC-binding factor (CTCF), and the possible implications on CTCF function and the spatial organization of the genome. Our findings suggest that mutation signature analysis may be useful in determining the origin of metastases in the neck and the parotid gland. Furthermore, UV-induced DNA damage to insulator binding sites may play a role in the carcinogenesis and progression of cSCC.

  Study EGAS00001003370

• Native American gene flow into Polynesia predating Easter Island settlement

  The possibility of voyaging contact between prehistoric Polynesian and Native American populations has long intrigued researchers. Proponents have pointed to the existence of New World crops, such as the sweet potato and bottle gourd, in the Polynesian archaeological record, but nowhere else outside the pre-Columbian Americas while critics have argued that these botanical dispersals need not have been human mediated. The Norwegian explorer Thor Heyerdahl controversially suggested that prehistoric South American populations had an important role in the settlement of east Polynesia and particularly of Easter Island (Rapa Nui). Several limited molecular genetic studies have reached opposing conclusions, and the possibility continues to be as hotly contested today as it was when first suggested. Here we analyse genome-wide variation in individuals from islands across Polynesia for signs of Native American admixture, analysing 807 individuals from 17 island populations and 15 Pacific coast Native American groups. We find conclusive evidence for prehistoric contact of Polynesian individuals with Native American individuals (around AD 1200) contemporaneous with the settlement of remote Oceania. Our analyses suggest strongly that a single contact event occurred in eastern Polynesia, before the settlement of Rapa Nui, between Polynesian individuals and a Native American group most closely related to the indigenous inhabitants of present-day Colombia.

  Study EGAS00001004209

• Reconstituting the transcriptome and DNA methylome landscapes of human implantation at single-cell resolution

  Implantation is a milestone event during mammalian embryogenesis. Due to the extreme difficulty of obtaining in vivo human early post-implantation embryos, the gene regulatory network and epigenetics controlling human embryo implantation remains elusive. Here, combining an in vitro culture system for human post-implantation development and single-cell omics sequencing technologies, over 10,000 single cells at five representative stages of pre/post-implantation development were systematically analyzed. Unsupervised dimensionality reduction and clustering algorithm of the transcriptome data show stepwise implantation routes for the epiblast, primitive endoderm, and trophectoderm lineages, suggesting preparation for the establishment of a mother-to-offspring connection after implantation. Female embryos showed asynchronous progress of dosage loss of X chromosomes during implantation. Furthermore, using the single cell trio-Seq (scTrio-Seq) strategy, re-methylation of the genomes of all the three lineages was unambiguously revealed. Surprisingly, the genome re-methylation of PE lineage were much slower than both EPI and TE lineages during the implantation process, indicating distinct methylome features between EPI and PE although both of which were derived from ICM. Collectively, our work paves the way for understanding the complex molecular mechanisms that regulate human embryo implantation, informing new insights and future efforts in early embryonic development and reproductive medicine.

  Study EGAS00001003443

• Molecular determinants of response to PD-L1 blockade across tumor types

  Immune checkpoint inhibitors targeting the PD-1/PD-L1 axis lead to durable clinical responses in subsets of cancer patients across multiple indications, including non-small cell lung cancer (NSCLC), urothelial carcinoma (UC) and renal cell carcinoma (RCC). Herein, we complement PD-L1 immunohistochemistry (IHC) and tumor mutation burden (TMB) with RNA-seq in 366 patients to identify unifying and indication-specific molecular profiles that can predict response to checkpoint blockade across these tumor types. Multiple machine learning approaches failed to identify a baseline transcriptional signature highly predictive of response across these indications. Signatures described previously for immune checkpoint inhibitors also failed to validate. At the pathway level, significant heterogeneity was observed between indications, in particular within the PD-L1+ tumors. mUC and NSCLC were molecularly aligned, with cell cycle and DNA damage repair genes associated with response in PD-L1- tumors. At the gene level, the CDK4/6 inhibitor CDKN2A was identified as a significant transcriptional correlate of response, highlighting the association of non-immune pathways to the outcome of checkpoint blockade. This cross-indication analysis revealed molecular heterogeneity between mUC, NSCLC and RCC tumors, suggesting that indication-specific molecular approaches should be prioritized to formulate treatment strategies.

  Study EGAS00001004343

• Large-scale viral genome analysis identifies novel clinical associations between hepatitis B virus and chronically infected patients

  Despite the high global prevalence of chronic hepatitis B (CHB) infection, datasets covering the whole hepatitis B viral genome from large patient cohorts are lacking, greatly limiting our understanding of the viral genetic factors involved in this deadly disease. We performed deep sequencing of viral samples from patients chronically infected with HBV to investigate the association between viral genome variation and patients’ clinical characteristics. We discovered novel viral variants strongly associated with viral load and HBeAg status. Patients with viral variants C1817T and A1838G had viral loads nearly three orders of magnitude lower than patients without those variants. These patients consequently experienced earlier viral suppression while on treatment. Furthermore, we identified novel variants that either independently or in combination with precore mutation G1896A were associated with the transition from HBeAg positive to negative phase of infection. These observations are consistent with the hypothesis that mutation of the HBeAg open reading frame is an important factor driving CHB patient’s HBeAg status. This analysis provides a detailed picture of HBV genetic variation in the largest patient cohort to date and highlights the diversity of plausible molecular mechanisms through which viral variation affects clinical phenotype.

  Study EGAS00001003689

• A system-wide approach to monitor responses to synergistic BRAF and EGFR inhibition in colorectal cancer cells

  Molecular heterogeneity of tumors, epigenetic changes and a diverse range of molecular mechanisms are main contributors to drug resistance, which represents one of the great challenges in cancer treatment. A deeper understanding of the molecular biology of cancer has resulted in better targeted therapies, where one or multiple drugs are adopted in novel therapies to tackle resistance. This beneficial effect of using combination treatments has to some extent also been observed in colorectal cancer (CRC) patients harboring the BRAF(V600E) mutation, whereby dual inhibition of BRAF(V600E) and EGFR increases antitumor activity. Notwithstanding this success, it is not clear whether this combination treatment is the only or most effective treatment to block resistance. Here, we investigate molecular responses upon single and multi-target treatments, over time, using BRAF(V600E) mutant CRC cells as a well-define model system. Through integration of transcriptomic, proteomic and phosphoproteomics data we obtain a comprehensive overview, revealing both well-known and novel responses. We primarily observe widespread upregulation of tyrosine kinases (RTKs) and and metabolic pathways. This points to by which the treated cells switch energy sources as a defensive response entering a quiescent like state, while activating signalling to re-activate the MAPK pathway.

  Study EGAS00001002654

• The_Genomic_Advances_in_Sepsis__GAinS__RNA_seq

  Sepsis is defined as life-threatening organ dysfunction caused by a dysregulated host response to infection. This cohort comprises a subset of patients enrolled in the Genomic Advances in Sepsis (GAinS) study, an established biobank of adult sepsis patients. Patients with sepsis due to community acquired pneumonia or faecal peritonitis were recruited from 34 hospitals across the UK from 2005-2018, with samples for functional genomics and detailed clinical information collected on the first, third and/or fifth day following ICU admission. RNA was extracted from leukocytes isolated at the bedside using LeukoLOCK kits. We have previously identified sepsis response signatures (SRSs), transcriptomic endotypes that are associated with differential early mortality (Davenport et al, Lancet Respir Med, 2016; Burnham et al, AJRCCM, 2017) and response to treatment in a clinical trial (Antcliffe et al, AJRCCM, 2018). We generated RNA sequencing data on 903 samples, including 134 samples repeated from our previously released microarray data. Libraries were prepared using NEB Ultra II Library Prep kits (Illumina) and sequenced on a NovaSeq 6000. Reads were aligned to the reference genome (GRCh38) using STAR and gene counts quantified using featureCounts (annotation Ensembl v99). Counts were TMM-normalised and log-transformed. Following QC, processed data were available on 864 samples from 667 unique patients.

  Study EGAS00001003772

• Extreme phenotypes define epigenetic and metabolic signatures in cardiovascular diseases.

  Improving the understanding of cardiometabolic syndrome pathophysiology and its relationship with thrombosis are ongoing healthcare challenges. Using plasma biomarkers analysis coupled with the transcriptional and epigenetic characterisation of cell types involved in thrombosis, obtained from two extreme phenotype groups (obese and lipodystrophy) and comparing these to lean individuals and blood donors, the present study identifies the molecular mechanisms at play, highlighting patterns of abnormal activation in innate immune phagocytic cells and shows that extreme phenotype groups could be distinguished from lean individuals, and from each other, across all data layers. The characterisation of the same obese group, six months after bariatric surgery shows the loss of the patterns of abnormal activation of innate immune cells previously observed. However, rather than reverting to the gene expression landscape of lean individuals, this occurs via the establishment of novel gene expression landscapes. Netosis and its control mechanisms emerge amongst the pathways that show an improvement after surgical intervention. Taken together, by integrating across data layers, the observed molecular and metabolic differences form a disease signature that is able to discriminate, amongst the blood donors, those individuals with a higher likelihood of having cardiometabolic syndrome, even when not presenting with the classic features.

  Study EGAS00001003780

• Single cell phenotypic profiling of 27 DLBCL cases reveals marked inter- and intra-tumoral heterogeneity

  Diffuse large B-cell lymphoma (DLBCL) is the most common histologic subtype of non-Hodgkin lymphoma and is notorious for its clinical heterogeneity. Patient outcomes can be predicted by cell-of-origin (COO) classification, demonstrating that the underlying transcriptional signature of malignant B-cells informs biological behavior in the context of standard combination chemotherapy regimens. In the current study, we used mass cytometry (CyTOF) to examine tumor phenotypes at the protein level with single cell resolution in a collection of 27 diagnostic DLBCL biopsy specimens from treatment naïve patients. We found that malignant B-cells from each patient occupied unique regions in 37-dimensional phenotypic space with no apparent clustering of samples into discrete subtypes. Interestingly, variable MHC class II expression was found to be the greatest contributor to phenotypic diversity. Within individual tumors, a subset of cases showed multiple phenotypic subpopulations, and in one case we were able to demonstrate direct correspondence between protein-level phenotypic subsets and DNA mutation-defined subclones. In summary, CyTOF analysis can resolve both inter- and intra-tumoral heterogeneity among primary samples, and reveals that each case of DLBCL is unique and may be comprised of multiple, genetically distinct subclones.

  Study EGAS00001003860

• Bisulfite-converted duplexes for the strand-specific detection and quantification of rare mutations

  The identification of mutations that are present at low frequencies in clinical samples is an essential component of precision medicine. The development of molecular barcoding for next generation sequencing has greatly enhanced the sensitivity of detecting such mutations by massively parallel sequencing. However, further improvements in specificity would be useful for a variety of applications. We herein describe a technology (BiSeqS) that can increase the specificity of sequencing by at least two orders of magnitude over and above that achieved with molecular barcoding and can be applied to any massively parallel sequencing instrument. BiSeqS employs bisulfite treatment to distinguish the two strands of molecularly barcoded DNA; its specificity arises from the requirement for the same mutation to be identified in both strands. Because no library preparation is required, the technology permits very efficient use of the template DNA as well as sequence reads, which are nearly all confined to the amplicons of interest. Such efficiency is critical for clinical samples, such as plasma, in which only tiny amounts of DNA are often available. We show here that BiSeqS can be applied to evaluate transversions, as well as small insertions or deletions, and can reliably detect one mutation among >10,000 wild type molecules.

  Study EGAS00001002406

• Biallelic mutations in the ubiquitin ligase RFWD3 cause Fanconi anemia

  The WD40-containing E3 ubiquitin ligase RFWD3 has been recently linked to the repair of DNA damage by Homologous Recombination (HR). Here we show that an RFWD3 mutation within the WD40 domain is connected to the genetic disease Fanconi anemia (FA). An individual revealed congenital abnormalities characteristic for FA. Cells from the patient, carrying the compound heterozygous mutations c.205_206dupCC and c.1916T>A in RFWD3, show increased sensitivity to DNA interstrand cross-linking agents in terms of increased chromosomal breakage, reduced survival and cell cycle arrest in G2 phase. The cellular phenotype is mirrored in genetically engineered human and avian cells by inactivation of RFWD3 or by introduction of the patient-derived missense mutation, and is rescued by expression of wildtype RFWD3 protein. HR is disrupted in RFWD3 mutant cells, caused by impaired relocation of mutant RFWD3 to chromatin and defective physical interaction with RPA. Rfwd3 knockout mice exhibit increased embryonic lethality, are sub-fertile, show ovarian and testicular atrophy and have a reduced life span in which they resemble other FA mouse models. Although mutation was detected in a single child with FA yet, we propose RFWD3 as a novel FA gene, FANCW, supported by cellular paradigm systems and an animal model.

  Study EGAS00001002440

• Genome-wide DNA Methylation is Predictive of Outcome in Juvenile Myelomonocytic Leukemia

  Juvenile myelomonocytic leukemia (JMML) is a myeloproliferative disorder of childhood caused by mutations in the Ras pathway. Outcomes in this disease vary dramatically from spontaneous resolution with little or no treatment to rapid relapse after hematopoietic stem cell transplantation. Given the high morbidity and late effects of transplant, it is critical to identify patients at diagnosis who can be observed rather than transplanted. We hypothesized that assessing DNA methylation status would help predict disease outcome. Genome-wide DNA methylation profiling using the Illumina 450k platform in a discovery cohort of 39 patients was performed. Unsupervised hierarchical clustering based on the most highly variable CpG sites identified three clusters of patients. Importantly, these clusters differed significantly in terms of 4-year event-free survival, with the lowest methylation cluster having the highest rates of survival. These findings were validated in an independent cohort of 40 patients. Of particular interest is that all but one of fourteen patients experiencing spontaneous resolution of their disease clustered together and closer to 22 healthy controls than the other JMML cases. This study demonstrates that DNA methylation patterns in JMML are predictive of outcome in this heterogeneously behaving disease and can identify patients who are most likely to experience spontaneous resolution.

  Study EGAS00001002700

• Hermansky-Pudlak syndrome type 1 causes impaired anti-microbial immunity through a pathogenic lipid metabolism-mTOR circuit

  Mendelian diseases that present with immune-mediated disorders can provide insights into the molecular mechanisms that drive inflammation. Hermansky-Pudlak syndrome (HPS) types 1 and 4 are caused by defective vesicle trafficking involving the BLOC-3 complex. The presence of inflammatory complications such as Crohn’s disease-like inflammation and lung fibrosis in these patients remains enigmatic. Using mass cytometry we observe an augmented inflammatory monocyte compartment in HPS1 patient peripheral blood that may be associated with a TNF - and IL-1α-dominated cytokine dysregulation. HPS1 patient monocyte-derived macrophages express an inflammatory TNF-OSM mRNA gene signature and changes in lipid metabolism. Using stimulation experiments and lysosomal proteomics we show that defective lipid metabolism drives RAB32-dependent mTOR signaling, facilitated by the accumulation of mTOR on lysosomes. This pathogenic circuit translates into aberrant bacterial clearance, which can be rescued with mTORC1 inhibition. We reveal that a pathogenic lipid-mTOR signaling circuit acts as a metabolic checkpoint for defective anti-microbial activity. This mechanism may be relevant to the complex pathology of HPS1 patients featuring macrophage lipid accumulation, granuloma formation, defective anti-microbial activity and tissue inflammation. Lastly, this circuit may be present in a wider group of disorders with defective lipid metabolism and cholesterol accumulation.

  Study EGAS00001005053

• THE GENOMIC LANDSCAPE OF ACTINIC KERATOSIS

  Actinic keratoses (AK) are lesions of epidermal keratinocyte dysplasia and are pre-cursor lesions to invasive cutaneous squamous cell carcinoma (cSCC). Identifying the specific genomic alterations driving the progression of normal skin to AK and invasive cSCC is challenging due to the massive, ultraviolet radiation-induced mutational burden characteristic to these lesions. Here, we present the largest whole exome sequencing study to date on AK with matched cSCC and demonstrate that AK and cSCC are almost indistinguishable at the genomic level, in terms of mutational burden, patterns of driver gene mutations and copy number alterations. We identified 44 significantly mutated AK driver genes through our established bioinformatics pipeline and demonstrate these genes are similarly mutated in cSCC. We also identified mutational signature-32 exclusively in AK from patients exposed to azathioprine, providing further, compelling evidence for this drug’s role in keratinocyte carcinogenesis, likely through its blocking of transcription coupled repair. We demonstrated that cSCC had higher levels of intra-sample heterogeneity than AK and that several signaling pathways, including immune-related signaling and TGF-beta signaling were significantly more mutated in cSCC. Integrating our findings with independent gene expression data confirms that dysregulated TGF-beta signaling may represent the critical pathway in the progression from AK to cSCC.

  Study EGAS00001004243

• Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell RCC -- CA209-010

  Immune checkpoint inhibitors targeting the PD-1 pathway have transformed the management of many advanced malignancies, including clear cell renal cell carcinoma (ccRCC), but the drivers and resistors of PD-1 response remain incompletely elucidated. Here, we analyzed 592 tumors collected from advanced ccRCC patients enrolled in prospective clinical trials of treatment with PD-1 blockade (or mTOR inhibition as control arm) by whole-exome and RNA-sequencing, integrated with immunofluorescence analysis, to define the somatic alteration landscape of late-stage ccRCC and to uncover the immunogenomic determinants of therapeutic response. While conventional genomic markers (tumor mutation burden, neoantigen load) and degree of CD8+ T cell infiltration were not associated with clinical response, we discovered numerous chromosomal alterations in advanced ccRCC associated with response or resistance to PD-1 blockade. These advanced tumors were highly CD8+ T cell infiltrated, with only 22% and 5% with an immune desert and immune excluded phenotype, respectively. Our analysis revealed that CD8+ infiltrated tumors are depleted of favorable PBRM1 mutations and are enriched for unfavorable chromosomal losses of 9p21.3 when compared to non-infiltrated tumors. These data demonstrate how the interplay of somatic alterations and immunophenotypes impacts therapeutic efficacy.

  Study EGAS00001004291

• This project aims to study human memory capacity, including short-term memory and long-term memory, systematically via genome-wide association studies

  Development of high-throughput genotyping platforms provides an opportunity to identify new genetic elements related to complex cognitive functions. Taking advantage of multi-level genomic analysis, here we studied the genetic basis of human short-term (STM, n=1,620) and long-term (LTM, n=1,526) memory functions. Heritability estimation based on single nucleotide polymorphism showed moderate heritability of short-term memory but very low heritability of long-term memory. In a two-step genome-wide association study, the markers rs13151012 and rs1558360 passed genome-wide significance (p < 5×10-8) in digit-span STM task and for the first principal component shared by two STM tasks; however, none of them survived the replication. In turn, we selected the ten most significant single nucleotide polymorphisms (SNPs) for replication tests. Among them, a SNP near ZFAT was significantly associated with STM performance in another independent population of 2,789; a polymorphism within BCAT2 was significantly associated with LTM in another independent population of 1,865. Furthermore, we performed a pathway analysis based on the current genomic data and found six pathways significantly associated with STM capacity and one pathway associated with LTM capacity.

  Study EGAS00001002875

• A Protein Neddylation Inhibitor MLN4924 Suppresses Patient-Derived Glioblastoma Cells via Inhibition of ERK and AKT Signaling

  Glioblastoma is a highly aggressive and lethal brain tumor, with limited treatment options. Abnormal activation of the neddylation pathway is observed in glioblastoma, and the NEDD8-activating enzyme (NAE) inhibitor, MLN4924 was previously shown to be effective in glioblastoma cell line models. However, its effect has not been tested in patient-derived glioblastoma stem cells. Immunohistochemistry and public data were analyzed to determine whether NEDD8 pathway proteins are important in glioblastoma development and patient survival. The effects of MLN4924 were evaluated in 4 glioblastoma cell lines and 15 patient-derived glioblastoma stem cells using high content analysis. Immunoblot analysis assessed the NEDD8 pathway proteins and biomarkers. Whole transcriptome analysis identified pathways associated with MLN4924 sensitivity. NAE1 and UBA3 levels increased in glioblastoma patients; high NEDD8 levels were associated with poor clinical outcomes. Glioblastoma cell lines and patient-derived stem cells were highly susceptible to MLN4924, while normal human astrocytes were resistant. Genomic analyses indicated that MLN4924 sensitive cells exhibited enrichment of ERK and AKT signaling. We verified that MLN4924 inhibits ERK and AKT phosphorylation in MLN4924 sensitive cells. Our findings suggest that patient-derived glioblastoma stem cells in the context of ERK and AKT activation are sensitive and highly regulated by neddylation inhibition.

  Study EGAS00001004018

• Hypoxia acts as an environmental cue for the human TRM differentiation program

  Tissue-resident memory T-cells (TRM) provide frontline defense against infectious diseases and contribute to anti-tumor immunity; however, aside from the necessity of TGF-β, knowledge regarding TRM-inductive cues remains incomplete, particularly for human cells. Oxygen tension is an environmental cue that distinguishes peripheral tissues from the circulation and here, we demonstrate that differentiation of human CD8+ T-cells in the presence of hypoxia and TGF-β1 led to the development of a TRM phenotype, characterized by a greater than five-fold increase in CD69+CD103+ cells expressing human TRM hallmarks and enrichment for endogenous human TRM gene signatures, including increased adhesion molecule expression and decreased expression of genes involved in recirculation. Hypoxia and TGF-β1 synergized to produce a significantly larger population of TRM phenotype cells than either condition alone, and comparison of these cells from the individual and combination conditions revealed distinct phenotypic and transcriptional profiles, indicating a programming response to milieu rather than a mere expansion. Our findings identify a previously unreported cue for the TRM differentiation program and can enable facile generation of human TRM phenotype cells in vitro for basic studies and translational applications such as adoptive cellular therapy.

  Study EGAS00001005286

• To determine the mutational impact of the in vitro culture, clonal human adult and pluripotent stem cell lines were subjected to a second clonal step after 3 months of culture. These subclones were whole genome sequenced to identify all the mutations that accumulated during the 3 month culture period.

  Genetic changes acquired during culture pose a potential risk for the successful application of stem cells. To assess the risk of in vitro expansion on mutation accumulation we have performed whole genome sequencing of clonally expanded human induced pluripotent stem cells (iPS cells) and adult stem cells (ASCs) to identify all mutations that accumulated over a fixed culture period. We find that ASCs acquire more single nucleotide variants (SNVs) and indels per population doubling than the iPS cells. When compared with ASCs, iPS cells are more vulnerable to mutations in genes and promoters. Mutational analysis revealed a clear in vitro induced mutational signature that is irrespective of stem cell type. This in vitro signature is characterized by C to G transversions that are probably caused by oxidative stress. Additionally, we observed stem cell specific mutational signatures and differences in transcriptional strand bias, indicating differential activity of DNA repair mechanisms between stem cell types in culture. In conclusion, in vitro culture induces mutation accumulation in iPS cells and ASCs. Culture under low-oxygen tension may help to reduce the number of culture-induced mutations.

  Study EGAS00001002955

• The Molecular Landscape of Asian Breast Cancers Reveals Clinically Relevant Population-Specific Differences

  Molecular profiles of breast cancer have contributed to an improved understanding of the disease, enabled development of molecular prognostic signatures to guide treatment decisions and unveiled new or more accurate therapeutic options for breast cancer patients. However, the extent to which differences in genetic, environmental and lifestyle factors influence molecular profiles in different populations remains poorly characterised, as relatively few large-scale molecular studies of breast tumours in non-Caucasian populations have hitherto been reported. Here, we present the molecular profiles of 560 Asian breast tumours and a comparative analysis of breast cancers arising in Asian and Caucasian women. Compared to the breast tumours in predominantly Caucasian women reported in TCGA and METABRIC, we show an increased prevalence of Her2-enriched molecular subtypes and higher prevalence of TP53 somatic mutations in ER+ Asian breast tumours. Using gene expression and immunohistochemistry, we observed elevated immune scores in Asian breast tumours, suggesting potential clinical response to immune checkpoint inhibitors. Whilst Her2-subtype and enriched immune score are associated with improved survival, presence of TP53 somatic mutations is associated with poorer survival in ER+ tumours. Taken together, these population differences unveil new opportunities to improve the understanding of this disease and lay the foundation for precision medicine in different populations.

  Study EGAS00001004518

• Host whole genome variations are associated with neurocognitive outcome in survivors of pediatric medulloblastoma

  Host whole genome analysis is a promising source of predictive information for long-term morbidity in cancer survivors. However, studies on genetic predictors of long-term outcome, particularly neurocognitive function following chemoradiation in pediatric oncology are limited. In the present study, we evaluated variations in host whole genome single nucleotide polymorphisms (SNPs) and their association with cognitive outcome. Whole-genome SNP analysis of host peripheral blood was conducted on 22 medulloblastoma long-term survivors, of whom 18 completed neuropsychological testing. First, unsupervised consensus clustering of the most variable SNPs within 409 genes involved in DNA repair was performed. Discrete variant groups were identified, although they were not associated with cognitive outcome, suggesting that variations in genes corresponding to a single functional group may be insufficient to predict long-term outcome alone. In support of this interpretation, unsupervised hierarchical clustering analysis using disease-associated gene variants by cognitive impairment status yielded two distinct variant clusters comprised of 36 variants, 34 of which were in noncoding regions. These findings illustrate for the first time that cognitively impaired survivors have a distinct variant profile compared to other medulloblastoma survivors. Future research in larger cohorts is needed to validate host genome predictors of cognitive impairment that may impact clinical management.

  Study EGAS00001002996

• Whole genome, exome and transcriptome analysis to guide individualized cancer interpretation

  Study Design: The Cancer Alliance study was a collaboration of nine academic medical institutions in the New York city area and led by the New York Genome Center. Tumor-normal paired samples were submitted by the institutions for whole genome, whole exome and transcriptome sequencing and bioinformatics analysis at New York Genome Center. Relevant clinical histories were also collected. Raw sequence files were available to all collaborators. Variant call files were analyzed both manually by interpreters at NYGC and by IBM Watson Genome Analytics. This study was approved by a central institutional review board (IRB), Biomedical Research Alliance of New York, and by local IRBs, including Stony Brook University and Northwell Health. The study was supported in part by a grant from the IBM corporation (IBM Watson Health) to the New York Genome Center, New York Genome Center philanthropic funds and Rockefeller University grant # UL1TR000043 from the National Center for Advancing Translational Sciences (NCATS), and the National Institutes of Health (NIH) Clinical and Translational Science Award (CTSA) program. The principal investigator Robert B. Darnell MD, PhD., is a Howard Hughes Medical Institute Investigator. Results: Pending Publication of Analysis of Study Conclusions: Pending Publication of Analysis of Study

  Study EGAS00001004013

• Hermansky-Pudlak syndrome type 1 causes impaired anti-microbial immunity through a pathogenic lipid metabolism-mTOR circuit - 10x Genomics scRNAseq

  Mendelian diseases that present with immune-mediated disorders can provide insights into the molecular mechanisms that drive inflammation. Hermansky-Pudlak syndrome (HPS) types 1 and 4 are caused by defective vesicle trafficking involving the BLOC-3 complex. The presence of inflammatory complications such as Crohn’s disease-like inflammation and lung fibrosis in these patients remains enigmatic. Using mass cytometry we observe an augmented inflammatory monocyte compartment in HPS1 patient peripheral blood that may be associated with a TNF - and IL-1α-dominated cytokine dysregulation. HPS1 patient monocyte-derived macrophages express an inflammatory TNF-OSM mRNA gene signature and changes in lipid metabolism. Using stimulation experiments and lysosomal proteomics we show that defective lipid metabolism drives RAB32-dependent mTOR signaling, facilitated by the accumulation of mTOR on lysosomes. This pathogenic circuit translates into aberrant bacterial clearance, which can be rescued with mTORC1 inhibition. We reveal that a pathogenic lipid-mTOR signaling circuit acts as a metabolic checkpoint for defective anti-microbial activity. This mechanism may be relevant to the complex pathology of HPS1 patients featuring macrophage lipid accumulation, granuloma formation, defective anti-microbial activity and tissue inflammation. Lastly, this circuit may be present in a wider group of disorders with defective lipid metabolism and cholesterol accumulation.

  Study EGAS00001005098

• Variation and transmission of the human gut microbiota across generations - shotgun data

  Although composition and functional potential of the human gut microbiota evolve over lifespan, kinship has been identified as a key covariate of microbial community diversification. To date, sharing of microbiota features within families has however mostly been assessed between parents and their direct offspring. Here, we investigate potential transmission and persistence of familial microbiome patterns and microbial genotypes in a family cohort (N=102) spanning three to five generations over the same female bloodline. We observe microbiome community composition to be associated with kinship, with seven (low-abundant) genera displaying familial distribution patterns. While kinship and current cohabitation emerged as closely entangled variables, our explorative analyses of microbial genotype distribution and transmission estimates point at the latter as a key covariate of strain dissemination. Highest potential transmission rates are estimated between sisters and mother-daughter pairs, decreasing with increasing daughter’s age, and being higher among cohabiting pairs than those living apart. Although rare, we do detect potential transmission events spanning three and four generations, primarily involving species of the genera Alistipes and Bacteroides. Overall, while our analyses confirm the existence of family-bound microbiome community profiles, transmission or co-acquisition of bacterial strains appears to be strongly linked to cohabitation.

  Study EGAS00001005649

• Variation and transmission of the human gut microbiota across generations - 16S data

  Although composition and functional potential of the human gut microbiota evolve over lifespan, kinship has been identified as a key covariate of microbial community diversification. To date, sharing of microbiota features within families has however mostly been assessed between parents and their direct offspring. Here, we investigate potential transmission and persistence of familial microbiome patterns and microbial genotypes in a family cohort (N=102) spanning three to five generations over the same female bloodline. We observe microbiome community composition to be associated with kinship, with seven (low-abundant) genera displaying familial distribution patterns. While kinship and current cohabitation emerged as closely entangled variables, our explorative analyses of microbial genotype distribution and transmission estimates point at the latter as a key covariate of strain dissemination. Highest potential transmission rates are estimated between sisters and mother-daughter pairs, decreasing with increasing daughter’s age, and being higher among cohabiting pairs than those living apart. Although rare, we do detect potential transmission events spanning three and four generations, primarily involving species of the genera Alistipes and Bacteroides. Overall, while our analyses confirm the existence of family-bound microbiome community profiles, transmission or co-acquisition of bacterial strains appears to be strongly linked to cohabitation.

  Study EGAS00001005651

• Nanopore_and_Illumina_sequencing_of_human_glioblastomas

  Liquid biopsy analysis refers to methods designed to detect tumour-specific material (e.g., circulating tumour cells or tumour DNA) in body fluids, such as urine or blood samples. A widely-used liquid biopsy approach consists of genotyping the circulating tumour DNA (ctDNA) through sequencing of plasma/serum DNA. Although useful in the context of metastatic disease (where the concentration of ctDNA is high), current liquid biopsy technologies show limited sensitivity of detection for the early detection of cancer, and low specificity, as cancer-related mutations accumulate in healthy tissues as part of the ageing process, thus making it difficult to distinguish these from tumour mutations, and because sequencing errors and true mutations show overlapping profiles. Therefore, developing liquid biopsy protocols with increased sensitivity and specificity represents an urgent clinical need. Here we harness extrachromosomal circular DNA (eccDNA) elements, which are circular DNA structures physically separated from the chromosomes of up to several Mbp long pervasive in human cancers, for liquid biopsy analysis. In this pilot study we will focus on the analysis of glioblastomas, because there is strong evidence for the presence of eccDNA in these tumour types, and because developing liquid biopsy approaches for brain tumours to reduce the invasiveness of brain tumour biopsies remains an unmet clinical need.

  Study EGAS00001005812

• A recessive PRDM13 mutation results in congenital hypogonadotropic hypogonadism and cerebellar hypoplasia

  The PRDM13 (PR Domain containing 13) putative chromatin modifier and transcriptional regulator functions downstream of the transcription factor PTF1A, which controls GABAergic fate in the spinal cord and neurogenesis in the hypothalamus. Here, we report a novel, recessive syndrome associated with PRDM13 mutation. Patients exhibited intellectual disability, ataxia with cerebellar hypoplasia, scoliosis and delayed puberty with congenital hypogonadotropic hypogonadism (CHH). Expression studies revealed Prdm13/PRDM13 transcripts in the developing hypothalamus and cerebellum in mouse and human. An analysis of hypothalamus and cerebellum development in mice homozygous for a Prdm13 mutant allele revealed a significant reduction in the number of Kisspeptin (Kiss1) neurons in the hypothalamus and PAX2+ progenitors emerging from the cerebellar ventricular zone. The latter was accompanied by ectopic expression of the glutamatergic lineage marker TLX3. Prdm13-deficient mice displayed cerebellar hypoplasia, normal gonadal structure, but delayed pubertal onset. Together, these findings identify PRDM13 as a critical regulator of GABAergic cell fate in the cerebellum and of hypothalamic kisspeptin neuron development, providing a mechanistic explanation for the co-occurrence of CHH and cerebellar hypoplasia in this syndrome. To our knowledge, this is the first evidence linking disrupted PRDM13-mediated regulation of Kiss1 neurons to CHH in humans.

  Study EGAS00001005878

• Sclerosing epithelioid fibrosarcoma case report

  Sclerosing epithelioid fibrosarcoma (SEF) is a rare and aggressive soft tissue sarcomathought to originate in fibroblasts of the tissues surrounding tendons, ligaments and muscles. Minimally responsive to conventional cytotoxic chemotherapies, greater than 50% of SEF patients experience relapse and/or metastatic disease. SEF is most commonly discovered in middle-aged and elderly adults, but also rarely in children. A common gene fusion occurring between the EWSR1 and CREB3L1 genes has been observed in 80-90% of SEF cases. We describe here the youngest SEF patient reported to date (a 3-year-old Caucasian male) who presented with numerous bony and lung metastases. Additionally, we perform a comprehensive literature review of all SEF-related articles published since the disease was first characterized. The patient described in this report experienced persistent disease progression despite being heavily treated with multiple surgeries, radiation, numerous chemotherapies and targeted therapeutics. Primary tumor, metastatic and relapse sites were sequenced by whole genome, whole exome, and deep transcriptome next generation sequencing with comparison to a patient-matched normal blood sample. Consistent across all sequencing analyses was the disease-defining EWSR1-CREB3L1 fusion as a single feature consensus. We provide an analysis of our genomic findings and a discussion as to potential therapeutic strategies for SEF.

  Study EGAS00001005215

• Comprehensive Genomic Characterization of Refractory Multiple Myeloma Reveals a Complex Mutational and Structural Landscape Associated with Drug Resistance (H067)

  Outcome of patients with multiple myeloma (MM) refractory to immunomodulatory agents (IMiDs) and proteasome inhibitors (PIs) is still very poor. To elucidate the genetic landscape of relapsed/refractory MM (RRMM) we performed whole genome and transcriptome sequencing of 39 heavily pretreated RRMM patients with at least double refractoriness revealing complex structural changes and a high mutational load (H067/K08K). We found a median of 67 structural variants and 116 nonsynonymous exonic mutations per patient. Assessment of oncogenic drivers showed NRAS, BRAF, TP53, KRAS, SLC4A7, FAM46C, DIS3, RB1, KMT2C, MLLT4, EWSR1, HCFC2, and COPS3 as significantly mutated. We found multiple genomic regions with bi-allelic events affecting known tumor suppressor genes, most frequently TP53, RB1, and TRAF3. Supervised analysis of mutational signatures showed a major contribution of Cosmic signature 6 associated with deficiency in DNA mismatch repair and the ‘BRCAness’ signature which coincided with mutations in genes known to confer PARP inhibitor sensitivity. Recurrently mutated gene networks included genes with relevance for PI and IMiD activity, the latter including in particular members of the Cereblon and the COP9 signalosome complex. In conclusion, this comprehensive genomic characterization revealed a complex mutational and structural landscape in RRMM and highlights potential implications for therapeutic strategies.

  Study EGAS00001004363

• tFL with a PMBL GE signature

  We investigated the clinicopathologic features of 5 follicular lymphomas (FL) that transformed to morphologic diffuse large B-cell lymphomas (DLBCL) and had a primary mediastinal large B-cell lymphoma (PMBL)-like gene expression profile. None of the transformed FL (tFL) arose in the mediastinum, all cases tested had a germinal center B-cell immunophenotype, 20% were CD30+, 60% CD23+, 80% MAL+, 20% CD200+, and 0% CD273/PDL2+. Whole exome sequencing detected alterations in genes associated with both FL/DLBCL (CREBBP, KMT2C, KMT2D, ARID1A, HIST1 members, and TNFRSF14) and PMBL (JAK-STAT pathway genes, B2M, and CD58). Copy number (CN) analysis detected gains/amplification of REL in 60% and STAT6 in 40%, gains of chromosome 16, including IL4R, in 40%, and both deletions and gains of 11q in 20%. This cohort, although limited in size, supports a subset of tFL that has blended features between FL/DLBCL and PMBL. Despite having some features that are less common in DLBCL (MAL and CD23 expression, JAK-STAT activation), these cases lack the most characteristic CN alteration seen in PMBL (9p24.1 gain/amplification). This cohort expands the biologic heterogeneity of tFL, illustrating a subset with blended FL/DLBCL and PMBL features, with potential treatment implications that include the use of novel PMBL-targeted therapies.

  Study EGAS00001005870

• The malaria-protective human glycophorin structural variant DUP4 shows somatic mosaicism and association with hemoglobin levels

  Glycophorin A and glycophorin B are red blood cell surface proteins that are both receptors for the parasite Plasmodium falciparum, which is the principal cause of malaria in sub-Saharan Africa. DUP4 is a complex structural genomic variant that carries extra copies of a glycophorin A - glycophorin B fusion gene, and has a dramatic effect on malaria risk by reducing the risk of severe malaria by up to 40% Using fiber-FISH and Illumina sequencing, we validate the structural arrangement of the glycophorin locus in the DUP4 variant, and reveal somatic variation in copy number of the glycophorin A- glycophorin B fusion gene. By developing a simple, specific, PCR-based assay for DUP4 we show the DUP4 variant reaches a frequency of 13% in a malaria-endemic village in south-eastern Tanzania. We genotype a substantial proportion of that village and demonstrate an association of DUP4 genotype with hemoglobin levels, a phenotype related to malaria, using a family-based association test. Taken together, we show that DUP4 is a complex structural variant that may be susceptible to somatic variation, and show that it is associated with a malarial-related phenotype in a non-hospitalizedlongitudinally-followed population.

  Study EGAS00001003239

• CRISPR/Cas9-mediated genome editing of Schistosoma mansoni acetylcholinesterase

  Here we report successful gene knock-in (KI) in the eggs of Schistosoma mansoni by combining CRISPR/Cas9 with single-stranded oligodeoxynucleotides (ssODNs). We targeted the acetylcholinesterase (AChE) gene of S. mansoni using two synthetic guide RNAs (gRNAs), X5 and X7, respectively. Liver eggs of S. mansoni were exposed to CRISPR-vector containing X5 or X7 by electroporation. Simultaneously, eggs were transfected with a ssODN donor encoding a stop codon in all six frames. Next generation sequencing analysis revealed that CRISPR/Cas9-mediated editing in S. mansoni eggs resulted in Homology-Directed Repair (HDR) when template DNA ssODN provided. Furthermore, soluble egg antigen (SEA) from AChE-modified eggs exhibited markedly reduced AChE activity compared with controls, indicative that programmed Cas9 cleavage mutated the AChE gene. Following injection of modified schistosome eggs into the tail veins of mice, a significant decrease in granuloma size in the lungs of these animals. Notably, an enhanced Th2 response induced by eggs in lung, and splenocytes small intestine-draining mesenteric lymph node cells was also generated in mice injected with X5-KI eggs in different methods. These findings further demonstrate the power and utility of CRISPR/Cas9-based genome editing for undertaking functional genomics studies in schistosomes.

  Study EGAS00001004455

• Sex chromosome aneuploidies give rise to changes in the circular RNA profile

  Purpose The landscape of circular RNAs (circRNAs), an important class of non-coding RNAs that regulate gene expression, has never been described in human disorders of sex chromosome aneuploidies. We profiled circRNAs in Turner syndrome females (45,X;TS) and Klinefelter syndrome males (47,XXY; KS) to investigate how circRNAs respond to a missing or an extra X chromosome. Methods Samples of blood, muscle and fat were collected from individuals with TS (n = 33) and KS (n = 22) and from male (n = 16) and female (n = 44) controls. CircRNAs were identified using a combination of circRNA identification pipelines (CIRI2, CIRCexplorer2 and circRNA_finder). Results Differential expression of circRNAs was observed throughout the genome in TS and KS, in all tissues. The host-genes from which several of these circRNAs were derived, were associated with known phenotypic traits. Furthermore, several differentially expressed circRNAs had the potential to capture micro RNAs that targeted protein-coding genes with altered expression in TS and KS. Conclusion Sex chromosome aneuploidies introduce pervasive changes in the circRNA transcriptome, demonstrating that the genomic changes in these syndromes are more complex than hitherto thought. CircRNAs may help explain some of the genomic and phenotypic traits observed in these syndromes.

  Study EGAS00001006404

• Multi-omics integration reveals only minor long-term molecular and functional sequelae in immune cells of individuals recovered from COVID-19

  The majority of COVID-19 patients experience mild to moderate disease course and recover within a few weeks. An increasing number of studies characterized the long-term changes in the specific anti-SARS-CoV-2 immune responses, but how COVID-19 shapes the innate and heterologous adaptive immune system after recovery is less well known. To comprehensively investigate the post-SARS-CoV-2 infection sequelae on the immune system, we performed a multi-omics study by integrating single-cell RNA-sequencing, single-cell ATAC-sequencing, genome-wide DNA methylation profiling, and functional validation experiments in convalescent COVID-19 and healthy individuals. We showed that immune responses generally recover without major sequelae after COVID-19. However, subtle differences persist at the transcriptomic level in monocytes, with downregulation of the interferon pathway, while DNA methylation also displays moderate changes in convalescent COVID-19 individuals. However, these differences did not affect the cytokine production capacity of PBMCs upon different bacterial, viral, and fungal stimuli, although baseline release of IL-1RA and IFN?? was higher in convalescent individuals. In conclusion, despite minor differences in epigenetic and transcriptional programs, the immune system of convalescent COVID-19 patients largely recovers to the homeostatic level of healthy individuals.

  Study EGAS00001005529

• Stereotyped B-cell responses are linked to IgG constant region polymorphisms in multiple sclerosis

  Clonally related B cells infiltrate the brain, meninges, and cerebrospinal fluid of multiple sclerosis (MS) patients, but the mechanisms driving the B-cell response and shaping the immunoglobulin repertoires remain unclear. Here, we used single-cell full-length RNA-seq and B-cell receptor reconstruction to simultaneously assess the phenotypes, isotypes, constant region polymorphisms, and the paired heavy- and light-chain repertoires in intrathecal B cells. We detected extensive clonal connections between the memory B cell and antibody-secreting cell compartments and observed clonally related cells of different isotypes, including IgM/IgG1, IgG1/IgA1, IgG1/IgG2, and IgM/IgA1. There was a strong dominance of the G1m1 allotype constant region polymorphisms in antibody-secreting cells, but not in memory B cells. Tightly linked to the G1m1 allotype, we found a preferential pairing of the immunoglobulin heavy-chain variable (IGHV)4 gene family with the κ variable (IGKV)1 gene family. The IGHV4-39 gene was most used and showed the highest frequency of pairing with IGKV1-5 and IGKV1(D)-33. These results link IgG constant region polymorphisms to stereotyped B-cell responses in MS and indicate that the intrathecal B-cell response in these patients could be directed against structurally similar epitopes.

  Study EGAS00001005745

• Regulatory T cell transcriptomic reprogramming characterizes adverse events by checkpoint inhibitors in solid tumors

  Immune checkpoint inhibitors (ICI), which target immune regulatory pathways to unleash antitumor responses, have revolutionized cancer immunotherapy. Despite the remarkable success of ICI immunotherapy, a significant proportion of patients whose tumors respond to these treatments develop immune-related adverse events (irAEs) resembling autoimmune diseases. Although the clinical spectrum of irAEs is well characterized, their successful management remains empiric. This is in part because the pathogenic mechanisms involved in the break-down of peripheral tolerance and induction of irAEs remain elusive. Herein, we focused on regulatory T cells (Tregs) in individuals with irAEs because these cells are vital for maintenance of peripheral tolerance, appear expanded in the peripheral blood of individuals with cancer and abundantly express checkpoint molecules, hence representing direct targets of ICI immunotherapy. Our data demonstrate an intense transcriptomic reprogramming of CD4+CD25+CD127– Tregs in the blood of individuals with advanced metastatic melanoma who develop irAEs following ICI immunotherapy, with a characteristic inflammatory, apoptotic and metabolic signature. This inflammatory signature was shared by Tregs from individuals with different types of cancer developing irAEs and individuals with autoimmune diseases. Our findings propose an inflammatory Treg reprogramming as a feature of immunotherapy-induced irAEs, that may facilitate translational approaches aiming to induce robust antitumor immunity without disturbing peripheral tolerance.

  Study EGAS00001004694

• Summary statistics from genome-wide association study in glioma of 12,488 cases and 18,169 controls.

  Genome-wide association studies (GWAS) have transformed our understanding of glioma susceptibility, but individual studies have had limited power to identify risk loci. We performed a meta-analysis of existing GWAS and two new GWAS, which totaled 12,488 cases and 18,169 controls. We identified five new loci for glioblastoma (GBM) at 1p31.3 (rs12752552; P = 2.04 × 10-9, odds ratio (OR) = 1.22), 11q14.1 (rs11233250; P = 9.95 × 10-10, OR = 1.24), 16p13.3 (rs2562152; P = 1.93 × 10-8, OR = 1.21), 16q12.1 (rs10852606; P = 1.29 × 10-11, OR = 1.18) and 22q13.1 (rs2235573; P = 1.76 × 10-10, OR = 1.15), as well as eight loci for non-GBM tumors at 1q32.1 (rs4252707; P = 3.34 × 10-9, OR = 1.19), 1q44 (rs12076373; P = 2.63 × 10-10, OR = 1.23), 2q33.3 (rs7572263; P = 2.18 × 10-10, OR = 1.20), 3p14.1 (rs11706832; P = 7.66 × 10-9, OR = 1.15), 10q24.33 (rs11598018; P = 3.39 × 10-8, OR = 1.14), 11q21 (rs7107785; P = 3.87 × 10-10, OR = 1.16), 14q12 (rs10131032; P = 5.07 × 10-11, OR = 1.33) and 16p13.3 (rs3751667; P = 2.61 × 10-9, OR = 1.18). These data substantiate that genetic susceptibility to GBM and non-GBM tumors are highly distinct, which likely reflects different etiology.

  Study EGAS00001003372

• RNA-seq study of human long-term and short-term hematopoietic stem cells from umblical cord blood with lentiviral overexpression of S1PR3

  To investigate the molecular and biological pathways altered by S1PR3OE in human hematopoietic stem cells (HSC), we performed RNA-sequencing (RNA-seq) of LT- and ST-HSC 3 days after transduction with control or S1PR3 overexpression (OE) lentiviral vectors. LT-HSC and ST-HSC from 3 pool of CB lin- were FACS-purified, cells were prestimulated for 4 hours and transduced with lentiviral vectors. At day 3, 2000-5300 BFP+ cells were FACS-purified for RNA isolation with a PicoPure kit. We were able to isolate only 1600-1800 BFP+ cells from LT-HSC control samples as opposed to 4000-5400 BFP+ cells from S1PR3OE samples. Thus, we pooled all control BFP+ LT-HSC cells into one sample for RNA-seq analysis. BFP- LT-HSC from control vector transduction were purified from CB1 as an additional LT-HSC control. Nextera libraries generated from 10 ng RNA from 5 LT-HSC samples (2 controls, 3 S1PR3OE) and 6 ST-HSC samples (3 controls, 3 S1PR3OE) were subjected to 125 bp, paired-end RNA-sequencing on the Illumina HiSeq 2500 with an average of 50 million reads/sample at the Center for Applied Genomics, Sick Kids Hospital.

  Study EGAS00001004798

• Proteogenomics of chronic lymphocytic leukemia

  Cancer heterogeneity at the proteome level may explain differences in therapy response and prognosis beyond the currently established genomic and transcriptomic based diagnostics. The relevance of proteomics for disease classifications remains to be established in clinically heterogeneous cancer entities such as chronic lymphocytic leukemia (CLL). Here, we characterized the proteome and transcriptome in-depth alongside genetic and ex-vivo drug response profiling in a clinically well annotated CLL discovery cohort (n= 68). Unsupervised clustering of the proteome data revealed six subgroups. Five of these proteomic groups were associated with genetic features, while one group was only detectable at the proteome level. This new group was characterized by accelerated disease progression, high spliceosomal protein abundances associated with aberrant splicing, and low B cell receptor signaling protein abundances (ASB-CLL). We developed classifiers to identify ASB-CLL based on its characteristic proteome or splicing signature in two independent cohorts (n= 165, n= 169) and confirmed that ASB-CLL comprises about 20 % of CLL patients. The inferior overall survival observed in ASB-CLL was independent of both TP53- and IGHV mutation status. Our multi-omics analysis refines the classification of CLL and highlights the potential of proteomics to improve cancer patient stratification beyond genetic and transcriptomic profiling. This is the transcriptomics data of the discovery cohort.

  Study EGAS00001005746

• Clinical genome sequencing uncovers potentially targetable truncations and fusions of MAP3K8 in spitzoid and other melanomas

  Spitzoid melanoma is a specific morphologic variant of melanoma that most commonly affects children and adolescents and ranges on the spectrum of malignancy from low-grade to overtly malignant. These tumors are generally driven by fusions of ALK, RET, NTRK1/3, MET, ROS1 and BRAF. However, in approximately 50% of cases, no genetic driver has been established. Clinical whole genome and transcriptome sequencing (RNA-Seq) of a spitzoid tumor from an adolescent revealed a novel gene fusion of MAP3K8, encoding a serine−threonine kinase that activates MEK. The patient who had exhausted all other therapeutic options was treated with a MEK inhibitor and underwent a transient clinical response. We subsequently analyzed spitzoid tumors from 49 patients by RNA-Seq and found in-frame fusions or C-terminal truncations in 33% of cases. The fusion transcripts and truncated genes all contained MAP3K8 exons 1-8 but lacked the autoinhibitory final exon. Data mining of RNA-Seq from the Cancer Genome Atlas (TCGA) uncovered analogous MAP3K8 rearrangements in 1.5% of adult melanomas. Thus, MAP3K8 rearrangements - uncovered by comprehensive clinical sequencing of a single case - are the most common genetic event in spitzoid melanoma, are present in adult melanomas, and could be amenable to MEK inhibition.

  Study EGAS00001003430

• Androgen receptor blockade promotes response to BRAF/MEK-targeted therapy

  We studied a group of melanoma patients treated with neoadjuvant BRAF/MEK-targeted therapy (NCT02231775, n=51), and observed significantly higher rates of major pathologic response (MPR= <10% viable tumor at resection) and improved recurrence-free survival (RFS) in females versus males (MPR-66% versus 14%, p=0.001; RFS-64% versus 32% at 2 years, p=0.021). Findings were validated in a several additional cohorts 2-4 patients with unresectable metastatic melanoma treated with BRAF and/or MEK-targeted therapy (n=664 patients in total), demonstrating improved progression-free survival (PFS) and overall survival (OS) in females versus males in several of these studies. Studies in pre-clinical models demonstrated significantly impaired anti-tumor activity in male versus female BRAF/MEK-treated mice (p=0.006), with significantly higher expression of androgen receptor (AR) in tumors of male and female BRAF/MEK-treated mice versus control (p=0.0006 and 0.0025). Pharmacologic inhibition of AR signaling improved responses to BRAF/MEK-targeted therapy in male and female mice (p=0.018 and p=0.003), whereas induction of AR signaling (via testosterone administration) was associated with significantly impaired response to BRAF/MEK-targeted therapy in males and females (p=0.021 and p<0.0001).

  Study EGAS00001006196

• WXS and RNA-seq for 22 patients treated with radiation + immunotherapy

  Over 500 active clinical trials are investigating combination radiation (RT) and immune checkpoint blockade (ICB) as a cancer treatment; however, the majority of trials have found no positive interaction. We performed a comprehensive molecular analysis of a randomized phase I clinical trial of patients with non-small cell lung cancer (NSCLC) treated with concurrent or sequential ablative RT and ICB. Concurrent treatment was superior to sequential treatment in augmenting local and distant tumor responses and improving overall survival in a subset of patients with highly aneuploid tumors, but not in those with less aneuploid tumors. Our analysis suggested that concurrent treatment eliminated immunologically cold aneuploid tumors better than sequential treatment and improved patient outcomes. In addition, we report that RT alone decreased intratumoral cytotoxic effector T cell and adaptive immune gene expression signatures, in contrast to upregulation of key immune pathways after concurrent administration of RT and ICB. Our findings distinguish the differential genomic and transcriptomic effects of RT versus RT and ICB and challenge the prevailing paradigm that local ablative RT beneficially stimulates the immune response. We propose the use of tumor aneuploidy as a novel biomarker and therapeutic target in personalizing treatment approaches for patients with NSCLC treated with RT and ICB.

  Study EGAS00001006212

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__clear_cell_renal_cell_carcinoma___2

  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. Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

  Study EGAS00001003542

• Whole genome sequencing delineates regulatory, structural, and cryptic splice variants in early onset cardiomyopathy

  Cardiomyopathy (CMP) is a heritable disorder. Over 50% cases are gene-elusive on clinical gene panel testing. The contribution of variants in non-coding DNA elements that result in cryptic splicing and regulate gene expression has not been explored. We analyzed whole genome sequencing (WGS) data in a discovery cohort of 209 pediatric CMP patients and 1,953 independent replication genomes and exomes. We searched for protein-coding variants, and non-coding variants predicted to affect the function or expression of genes. Thirty-nine % cases harbored pathogenic coding variants in known CMP genes, and 5% harbored high-risk loss-of-function (LoF) variants in additional candidate CMP genes. Fifteen % harbored high-risk regulatory variants in promoters and enhancers of CMP genes (Odds ratio 2.25, p=6.70×10-7 versus controls). Genes involved in α-dystroglycan glycosylation (FKTN, DTNA) and desmosomal signaling (DSC2, DSG2) were most highly enriched for regulatory variants (Odds ratio 6.7-58.1). Functional effects were confirmed in patient myocardium and reporter assays in human cardiomyocytes, and in zebrafish CRISPR knockouts. We provide strong evidence for the genomic contribution of functionally active variants in new genes and in regulatory elements of known CMP genes to early-onset CMP.

  Study EGAS00001004929

• Single-cell RNA sequencing reveals the mesangial identity and species diversity of glomerular cell transcriptomes

  Molecular characterization of the individual cell types in human kidney as well as model organisms are critical in defining organ function and understanding translational aspects of biomedical research. Previous studies have uncovered gene expression profiles of several kidney glomerular cell types, however, important cells, including mesangial (MCs) and glomerular parietal epithelial cells (PECs), were missing or incompletely described, and a systematic comparison between mouse and human kidney is lacking. To this end, we used Smart-seq2 to profile 4332 individual glomerulus-associated cells isolated from human living donor renal biopsies and mouse kidney. The analysis revealed genetic programs for all four glomerular cell types (podocytes, glomerular endothelial cells, MCs and PECs) as well as rare glomerulus-associated macula densa cells (MDCs). Importantly, we detected heterogeneity in glomerulus-associated Pdgfrb-expressing cells, including bona fide intraglomerular MCs with the functionally active phagocytic molecular machinery, as well as a unique mural cell type located in the central stalk region of the glomerulus tuft. Furthermore, we observed remarkable species differences in the individual gene expression profiles of defined glomerular cell types that highlight translational challenges in the field and provide a guide to design translational studies.

  Study EGAS00001004943

• Significant and pervasive effects of RNA degradation on Nanopore direct RNA sequencing

  Abstract Oxford Nanopore direct RNA sequencing (DRS) is capable of sequencing complete RNA molecules and accurately measuring gene and isoform expression. However, as DRS is designed to profile intact RNA, expression quantification may be more heavily dependent upon RNA integrity than alternate RNA-seq methodologies. It is currently unclear how RNA degradation impacts DRS or if it can be corrected for. To assess the impact of RNA integrity on DRS we performed a degradation time-series using SH-SY5Y neuroblastoma cells. Our results demonstrate that degradation is a significant and pervasive factor that can bias DRS measurements, including a reduction in library complexity resulting in an overrepresentation of short genes and isoforms. Degradation also biases differential expression analyses; however, we find explicit correction can almost fully recover meaningful biological signal. In addition, DRS provided less biased profiling of partially degraded samples than nanopore cDNA-PCR sequencing. Overall, we find samples with RIN > 9.5 can be treated as undegraded and samples with RIN > 7 can be utilised for DRS with appropriate correction. These results establish the suitability of DRS for a wide range of samples, including partially degraded in-vivo clinical and post-mortem samples, whilst limiting the confounding effect of degradation on expression quantification.

  Study EGAS00001006542

• Interplay of somatic alterations and immune infiltration modulates response to PD-1 blockade in advanced clear cell RCC-- CA209-025

  Immune checkpoint inhibitors targeting the PD-1 pathway have transformed the management of many advanced malignancies, including clear cell renal cell carcinoma (ccRCC), but the drivers and resistors of PD-1 response remain incompletely elucidated. Here, we analyzed 592 tumors collected from advanced ccRCC patients enrolled in prospective clinical trials of treatment with PD-1 blockade (or mTOR inhibition as control arm) by whole-exome and RNA-sequencing, integrated with immunofluorescence analysis, to define the somatic alteration landscape of late-stage ccRCC and to uncover the immunogenomic determinants of therapeutic response. While conventional genomic markers (tumor mutation burden, neoantigen load) and degree of CD8+ T cell infiltration were not associated with clinical response, we discovered numerous chromosomal alterations in advanced ccRCC associated with response or resistance to PD-1 blockade. These advanced tumors were highly CD8+ T cell infiltrated, with only 22% and 5% with an immune desert and immune excluded phenotype, respectively. Our analysis revealed that CD8+ infiltrated tumors are depleted of favorable PBRM1 mutations and are enriched for unfavorable chromosomal losses of 9p21.3 when compared to non-infiltrated tumors. These data demonstrate how the interplay of somatic alterations and immunophenotypes impacts therapeutic efficacy.

  Study EGAS00001004292

• T cell responses of ALS patients

  In this project, we aimed to understand how T cell responses contribute to the disease progression of amyotrophic lateral sclerosis (ALS). The present data is on single-cell sequencing isolated from human cerebrospinal fluid (CSF) cells from both ALS patients (n=5) and controls (n=4). This analysis was conducted as part of a bigger project which is summarized in the section below. We used flow cytometry to define T cell subsets and phenotypes in blood and CSF samples collected at the time of diagnosis on a cohort of 89 newly diagnosed ALS patients in Stockholm, Sweden. High frequency of CD4+FOXP3- effector T cells in blood and CSF was associated with poor survival whereas high frequency of activated regulatory T (Treg) cells and high ratio between activated and resting Treg cells in blood was associated with better survival. T cell profiles also predicted disease progression rate. On an independent cohort of cases and controls, we used single cell transcriptomics data to demonstrate that ALS patients had altered T cell gene expression patterns and clonally expanded CD4+ and CD8+ T cells in CSF. In summary, T cell responses contribute to disease progression of ALS, supporting modulation of adaptive immunity as a viable therapeutic option.

  Study EGAS00001006675

• Longitudinal evaluation of serum microRNAs as biomarkers for neuroblastoma burden and therapeutic p53 reactivation

  Accurate assessment of treatment response and residual disease is indispensable for the evaluation of cancer treatment efficacy. However, performing tissue biopsies for longitudinal follow-up poses a major challenge in the management of solid tumours like neuroblastoma. In the present study, we evaluated whether circulating miRNAs are suitable to monitor neuroblastoma tumour burden and whether treatment-induced changes of miRNA abundance in the tumour are detectable in serum. We performed small RNA sequencing on longitudinally collected serum samples from mice carrying orthotopic neuroblastoma xenografts that were exposed to treatment with idasanutlin or temsirolimus. We identified 57 serum miRNAs to be differentially expressed upon xenograft tumour manifestation, out of which 21 were also found specifically expressed in the serum of human high-risk neuroblastoma patients. The murine serum levels of these 57 miRNAs correlated with tumour tissue expression and tumour volume, suggesting potential utility for monitoring tumour burden. In addition, we describe serum miRNAs that dynamically respond to p53 activation following treatment of engrafted mice with idasanutlin. We identified idasanutlin-induced serum miRNA expression changes upon one day and 11 days of treatment. By limiting to miRNAs with a tumour-related induction, we put forward hsa-miR-34a-5p as a potential pharmacodynamic biomarker of p53 activation in serum.

  Study EGAS00001006678

• Hexanucleotide repeat expansions in C9orf72 alter microglial responses and prevent a coordinated glial reaction in ALS

  Neuroinflammation is an important hallmark in amyotrophic lateral sclerosis (ALS). Experimental evidence has highlighted a role of microglia in the modulation of motor neuron degeneration. However, the exact contribution of microglia to both sporadic and genetic forms of ALS is still unclear. We generated single nuclei profiles of spinal cord and motor cortex from sporadic and C9orf72 ALS patients, as well as controls. We particularly focused on the transcriptomic responses of both microglia and astrocytes. We confirmed that C9orf72 is highly expressed in microglia and shows a diminished expression in carriers of the hexanucleotide repeat expansion (HRE). This resulted in an impaired response to disease, with specific deficits in phagocytic and lysosomal transcriptional pathways. Astrocytes also displayed a dysregulated response in C9orf72 ALS patients, remaining in a homeostatic state. This suggests that C9orf72 HRE alters a coordinated glial response, which ultimately would increase the risk for developing ALS. Our results indicate that C9orf72 HRE results in a selective microglial loss-of-function, likely impairing microglial-astrocyte communication and preventing a global glial response. This is relevant as it indicates that sporadic and familial forms of ALS may present a different cellular substrate, which is of great importance for patient stratification and treatment.

  Study EGAS00001006711

• T cell landscape definition by multi-omics identifies galectin-9 as novel immunotherapy target in chronic lymphocytic leukemia (CLL)

  Failure of immunotherapy after applying checkpoint inhibitors or CAR-T cells is linked to T cell exhaustion. Here, we explored the T cell landscape in chronic lymphocytic leukemia (CLL) by single-cell omics analyses of blood, bone marrow and lymph node samples of patients and spleen samples of a CLL mouse model. By single-cell RNA- sequencing, mass cytometry (CyTOF), and multiplex image analysis of tissue microarrays, we defined the spectrum of phenotypes and transcriptional programs of T cells and and their differentiation state trajectories. We identified disease-specific accumulation of distinct regulatory T cell subsets and T cells harboring an exhausted phenotype exclusively in the CLL lymph node tissue. Integration of TCR data revealed a clonal expansion of CD8+ precursor exhausted T cells, suggesting their reactivity for CLL cells. Interactome analyses identified the TIM3 ligand Galectin-9 as novel immunoregulatory molecule in CLL. Blocking of Galectin-9 in CLL-bearing mice slowed down disease development and reduced the number of TIM3 expressing T cells. Galectin-9 expression correlated with shorter survival of CLL patients. Thus, Galectin-9 contributes to immune escape in CLL and represents a novel target for immunotherapy.

  Study EGAS00001006864

• Atypical teratoid/rhabdoid tumoroids reveal subgroup-specific drug vulnerabilities WGS

  Atypical teratoid/rhabdoid tumors (ATRTs) represent a rare, but aggressive pediatric brain tumor entity. They are genetically uniquely defined by alterations in the SWI/SNF chromatin remodeling complex members SMARCB1 or SMARCA4. ATRTs can be further classified in different molecular subgroups based on their epigenetic profiles. Although recent studies suggest that the different subgroups have distinct clinical features, subgroup-specific treatment regimens have not been developed thus far. This is hampered by the lack of pre-clinical in vitro models representative of the different molecular subgroups. Here, we describe the establishment of ATRT tumoroid models from the ATRT-MYC and ATRT-SHH subgroups. We demonstrate that ATRT tumoroids retain subgroup-specific epigenetic and gene expression profiles. High throughput drug screens on our ATRT tumoroids revealed distinct drug sensitivities between and within ATRT-MYC and ATRT-SHH subgroups. Whereas ATRT-MYC universally displayed high sensitivity to multi-targeted tyrosine kinase inhibitors, ATRT-SHH showed a more heterogeneous response with a subset showing high sensitivity to Notch inhibitors, which corresponded to high expression of Notch receptors. Our ATRT tumoroids represent the first pediatric brain tumor organoid model, providing a representative pre-clinical model which enables the development of subgroup-specific therapies.

  Study EGAS00001006865

• Atypical teratoid/rhabdoid tumoroids reveal subgroup-specific drug vulnerabilities RNA-Seq

  Atypical teratoid/rhabdoid tumors (ATRTs) represent a rare, but aggressive pediatric brain tumor entity. They are genetically uniquely defined by alterations in the SWI/SNF chromatin remodeling complex members SMARCB1 or SMARCA4. ATRTs can be further classified in different molecular subgroups based on their epigenetic profiles. Although recent studies suggest that the different subgroups have distinct clinical features, subgroup-specific treatment regimens have not been developed thus far. This is hampered by the lack of pre-clinical in vitro models representative of the different molecular subgroups. Here, we describe the establishment of ATRT tumoroid models from the ATRT-MYC and ATRT-SHH subgroups. We demonstrate that ATRT tumoroids retain subgroup-specific epigenetic and gene expression profiles. High throughput drug screens on our ATRT tumoroids revealed distinct drug sensitivities between and within ATRT-MYC and ATRT-SHH subgroups. Whereas ATRT-MYC universally displayed high sensitivity to multi-targeted tyrosine kinase inhibitors, ATRT-SHH showed a more heterogeneous response with a subset showing high sensitivity to Notch inhibitors, which corresponded to high expression of Notch receptors. Our ATRT tumoroids represent the first pediatric brain tumor organoid model, providing a representative pre-clinical model which enables the development of subgroup-specific therapies.

  Study EGAS00001006866

• Atypical teratoid/rhabdoid tumoroids reveal subgroup-specific drug vulnerabilities DNA-Methylation

  Atypical teratoid/rhabdoid tumors (ATRTs) represent a rare, but aggressive pediatric brain tumor entity. They are genetically uniquely defined by alterations in the SWI/SNF chromatin remodeling complex members SMARCB1 or SMARCA4. ATRTs can be further classified in different molecular subgroups based on their epigenetic profiles. Although recent studies suggest that the different subgroups have distinct clinical features, subgroup-specific treatment regimens have not been developed thus far. This is hampered by the lack of pre-clinical in vitro models representative of the different molecular subgroups. Here, we describe the establishment of ATRT tumoroid models from the ATRT-MYC and ATRT-SHH subgroups. We demonstrate that ATRT tumoroids retain subgroup-specific epigenetic and gene expression profiles. High throughput drug screens on our ATRT tumoroids revealed distinct drug sensitivities between and within ATRT-MYC and ATRT-SHH subgroups. Whereas ATRT-MYC universally displayed high sensitivity to multi-targeted tyrosine kinase inhibitors, ATRT-SHH showed a more heterogeneous response with a subset showing high sensitivity to Notch inhibitors, which corresponded to high expression of Notch receptors. Our ATRT tumoroids represent the first pediatric brain tumor organoid model, providing a representative pre-clinical model which enables the development of subgroup-specific therapies.

  Study EGAS00001006881

• Human leukocyte antigen alleles associate with COVID-19 vaccine immunogenicity and risk of breakthrough infection

  SARS-CoV-2 vaccine immunogenicity varies between individuals, and immune responses correlate with vaccine efficacy. Using data from 1,076 participants enrolled in ChAdOx1 nCov-19 vaccine efficacy trials in the United Kingdom, we find that inter-individual variation in normalised antibody responses against SARS-CoV-2 spike (S) and its receptor binding domain (RBD) at 28 days following first vaccination shows genome-wide significant association with major histocompatibility complex (MHC) class II alleles. The most statistically significant association with higher levels of anti-RBD antibody was HLA-DQB1*06 (P=3.2 x 10-9), which we replicate in 1,677 additional vaccinees. Individuals carrying HLA-DQB1*06 alleles were less likely to experience PCR-confirmed breakthrough infection during the Wuhan and Alpha-variant waves compared to non-carriers (HR 0.63, 0.42-0.93, P=0.02). We identify a distinct S-derived peptide that is predicted to bind effectively to HLA-DQB1*06 compared to other equivalent alleles, and find evidence of increased spike-specific memory B-cell responses in HLA-DQB1*06 carriers at 84 days post first vaccination. Our results demonstrate association of HLA type with COVID-19 vaccine antibody response and risk of breakthrough infection, with implications for future vaccine design and implementation.

  Study EGAS00001006909

• Quantitative analysis of a novel DNA hypermethylation panel for lung cancer diagnosis using bronchial specimens

  Non-diagnostic findings in Transbronchial Lung Biopsy (TBLB) and Endobronchial Ultrasound-guided Transbronchial Lung Biopsy (EBUS-TBLB) are not uncommon. A challenge is to improve lung cancer detection by TBLB/EBUS-TBLB. In order to improve the diagnostic yield of bronchoscopy, we used 850K methylation chip to identify methylation sites that distinguish malignant from benign lung nodules.We found that the combination of HOXA7, SHOX2,and SCT methylation analysis has the best diagnostic yield in 123 bronchial washing (sensitivity: 74·1%; AUC: 0·851) and 172 brushing samples (sensitivity: 86·1%; AUC: 0·915). We developed a kit comprising these three genes and validated the kit in 329 unique bronchial washing samples, 397 unique brushing samples, and 179 unique patients with both washing and brushing samples. The diagnostic accuracy of the panel alone in lung cancer diagnosis was 86·9%, 91·2% and 95% in bronchial washing, brushing, and washing+brushing samples, respectively. The panel's sensitivity in combination with cytology, rapid on-site evaluation (ROSE), and histology in lung cancer diagnosis was 90·8% and 95.8% in bronchial washing and brushing samples, respectively, and 100% in washing+brushing samples. Quantitative analysis of the three-gene panel improves lung cancer diagnosis by bronchoscopy.

  Study EGAS00001007089

• Long-term organoid culture of a small intestinal neuroendocrine tumor rna-seq

  Gastroenteropancreatic neuroendocrine tumors (GEP-NETs) are rare and highly heterogeneous neoplasms whose incidence has markedly increased over the last decades. A grading system based on the tumor cells' proliferation index predicts high-risk for G3 NETs. However, low-to-intermediate grade (G1/G2) NETs have an unpredictable clinical course that varies from indolent to highly malignant. Cultures of human cancer cells enable to perform functional perturbation analyses that are instrumental to enhance our understanding of cancer biology. To date, no tractable and reliable long-term culture of G1/G2 NET has been reported to permit disease modeling and pharmacological screens. Here, we report of the first long-term culture of a G2 metastatic small intestinal NET that preserves the main genetic drivers of the tumor and retains expression patterns of the endocrine cell lineage. Replicating the tissue, this long-term culture showed a low proliferation index, and yet it could be propagated continuously without dramatic changes in the karyotype. The model was readily available for pharmacological screens using targeted agents and as expected, showed low tumorigenic capacity in vivo. Overall, this is the first long-term culture of NETs to faithfully recapitulate many aspects of the original neuroendocrine tumor.

  Study EGAS00001007108

• Transcriptomics identifies blunted immunomodulatory effects of vitamin D in people with multiple sclerosis

  Vitamin D deficiency is a risk factor for developing multiple sclerosis (MS). However, the immune effects of vitamin D in people with MS are not well understood. We analyzed transcriptomic datasets generated by RNA sequencing of immune cell subsets (CD4+, CD8+ T cells, B cells, monocytes) from 33 healthy controls and 33 untreated MS cases. We utilized a traditional bioinformatic pipeline and weighted gene co-expression network analysis (WGCNA) to determine genes and pathways correlated with endogenous vitamin D. In controls, CD4+ and CD8+ T cells had 1079 and 1188 genes, respectively, whose expressions were correlated with plasma 25-hydroxyvitamin D level (P<0.05). Functional enrichment analysis identified association with TNF-alpha and MAPK signaling. In CD4+ T cells of controls, vitamin D level was associated with expression levels of several genes proximal to multiple sclerosis risk loci (P=0.01). Genes differentially associated with endogenous vitamin D by case-control status were enriched in TNF-alpha signaling via NF-κB. WGCNA suggested a blunted response to vitamin D in cases relative to controls. Collectively, our findings provide further evidence for the immune effects of vitamin D, and demonstrate a differential immune response to vitamin D in cases relative to controls, highlighting a possible mechanism contributing to MS pathophysiology.

  Study EGAS00001007254

• Targeting TRIP13 in Wilms Tumor with Nuclear Export Inhibitors

  Wilms tumor (WT) is the most common renal malignancy of childhood. Despite improvements in the overall survival, relapse occurs in ~15% of patients with favorable histology WT (FHWT). Half of these patients will succumb to their disease. Identifying novel targeted therapies in a systematic manner remains challenging in part due to the lack of faithful preclinical in vitro models. We established ten short-term patient-derived WT cell lines and characterized these models using low-coverage whole genome sequencing, whole exome sequencing and RNA-sequencing, which demonstrated that these ex-vivo models faithfully recapitulate WT biology. We then performed targeted RNAi and CRISPR-Cas9 loss-of-function screens and identified the nuclear export genes (XPO1 and KPNB1) as strong vulnerabilities. We observed that these models are sensitive to nuclear export inhibition using the FDA approved therapeutic agent, selinexor (KPT-330). Selinexor treatment of FHWT suppressed TRIP13 expression, which was required for survival. We further identified in vitro and in vivo synergy between selinexor and doxorubicin, a chemotherapy used in high risk FHWT. Taken together, we identified XPO1 inhibition with selinexor as a potential therapeutic option to treat FHWTs and in combination with doxorubicin, leads to durable remissions in vivo.

  Study EGAS00001007389

• Genomic evolution and transcriptional changes in the evolution of prostate cancer into neuroendocrine and ductal carcinoma types

  Prostate cancer is typically of acinar adenocarcinoma type but can occasionally present as neuro-endocrine and/or ductal type carcinoma. These are associated with clinically aggressive disease and the former often arises on a background of androgen deprivation therapy, although it can al-so arise de novo. Two prostate cancer cases were sequenced by exome capture from archival tis-sue. Case 1 was de novo small cell neuroendocrine carcinoma and ductal adenocarcinoma with 3 longitudinal samples over 5 years. Case 2 was a single time point after development of treatment related neuroendocrine prostate carcinoma. Case 1 showed whole genome doubling in all sam-ples and focal amplification of AR in all samples except the first time point. Phylogenetic analysis revealed a common ancestry for the ductal and small cell carcinoma. Case 2 showed 13q loss (involving RB1) in both adenocarcinoma and small cell carcinoma regions, and 3p gain, 4p loss, 17p loss (involving TP53) in the latter. By using highly curated samples, we demonstrate for the first time that small cell neuroendocrine and ductal prostatic carcinoma can have a common an-cestry and the process of evolution over time. We highlight whole genome doubling in a patient with prostate cancer relapse, reinforcing its poor prognostic nature.

  Study EGAS00001007412

• Genomic evolution and transcriptional changes in the evolution of prostate cancer into neuroendocrine and ductal carcinoma types (RNAseq)

  Prostate cancer is typically of acinar adenocarcinoma type but can occasionally present as neuro-endocrine and/or ductal type carcinoma. These are associated with clinically aggressive disease and the former often arises on a background of androgen deprivation therapy, although it can al-so arise de novo. Two prostate cancer cases were sequenced by exome capture from archival tis-sue. Case 1 was de novo small cell neuroendocrine carcinoma and ductal adenocarcinoma with 3 longitudinal samples over 5 years. Case 2 was a single time point after development of treatment related neuroendocrine prostate carcinoma. Case 1 showed whole genome doubling in all sam-ples and focal amplification of AR in all samples except the first time point. Phylogenetic analysis revealed a common ancestry for the ductal and small cell carcinoma. Case 2 showed 13q loss (involving RB1) in both adenocarcinoma and small cell carcinoma regions, and 3p gain, 4p loss, 17p loss (involving TP53) in the latter. By using highly curated samples, we demonstrate for the first time that small cell neuroendocrine and ductal prostatic carcinoma can have a common an-cestry and the process of evolution over time. We highlight whole genome doubling in a patient with prostate cancer relapse, reinforcing its poor prognostic nature.

  Study EGAS00001007428

• Integrative Multi-Omics and Drug Sensitivity Profiling Reveals Potential Predictive Biomarkers in Pediatric Solid Tumors from the INFORM Registry

  Cure rates for childhood malignancies using established therapy protocols have increased to an average of 80% but have reached a plateau. Moreover, survival rates are particularly low for some pediatric tumors (such as high-risk Group3 medulloblastomas, osteosarcomas, Ewing sarcomas, high-risk neuroblastomas and high-grade gliomas) and dismal for patients with relapsed malignancies. A functional drug response profiling platform for pediatric solid and brain tumors has been established within the INFORM program to identify patient-specific vulnerabilities, biomarkers and unravel molecular mechanisms associated with drug response profiles for clinical translation. Here, we perform a multi-omics analysis using drug sensitivity profiles, genomic, and transcriptomic data of 81 pediatric solid tumor samples. The integrative analysis revealed two multi-omics signatures with high drug sensitivity predictivity. One signature highlighted neuroblastoma samples with sensitivity to navitoclax (a BCL-2 family inhibitor). A second signature was significantly specific to a subset of Wilms tumors harboring the SIX1 (Q177R) hotspot mutation, which had high expression of MGAM, PTPN14, STAT4 and KDM2B and high sensitivity to MEK inhibitors. A patient-specific causal interaction network analysis elucidated the potential biological mechanism of action between MEK inhibitors and the SIX1 mutation in n = 6 Wilms tumor samples.

  Study EGAS00001008249

• Targeted sequencing of diffuse large B-cell lymphoma tumour samples from the UK's Haematological Malignancy Research Network

  Targeted sequencing was applied to an unselected population-based diffuse large B-cell lymphoma cohort (n=928) diagnosed in the UK's Haematological Malignancy Research Network catchment population of ~4 million (14 centres). DNA extracted from tumour samples was sequenced with a 293-gene panel using the Illumina HiSeq 2500. All data are provided in the CRAM format.

  Dataset EGAD00001008616

• Neutrophils and emergency granulopoiesis drive immune suppression and an extreme response endotype during sepsis

  This dataset contains raw sequencing files associated with the paper "Neutrophils and emergency granulopoiesis drive immune suppression and an extreme response endotype during sepsis" (https://doi.org/10.1038/s41590-023-01490-5). It is composed of two experiments, which are as follows: 1. Single-cell profiling of whole blood leukocytes: This experiment comprises sequencing files generated using the BD Rhapsody single-cell multi-omics profiling (RNA + protein) platform. This platform was used to profile the whole blood leukocyte population in a cohort of sepsis patients, cardiac surgery controls and healthy controls. There are 48 samples and four files per sample: two paired-end FASTQ files (R1 and R2) corresponding to the RNA-seq library, and two paired-end FASTQ files (R1 and R2) corresponding to the protein profiling (ADT-based AbSeq) library. 2. Single-cell profiling of circulating HSPCs: This experiment comprises sequencing files generated using the 10X single-cell multi-omics (RNA-seq + ATAC-seq) profiling platform. This platform was used to profile circulating HSPCs in a cohort of sepsis patients and healthy controls. There are 5 samples (or plexes), each of which consists of a pool of individuals for whom cells were multiplexed and sequenced as a single sample. There are four files per sample: two paired-end FASTQ files (R1 and R2) corresponding to the RNA-seq library, and two paired-end FASTQ files (R1 and R2) corresponding to the ATAC-seq library. Because these samples consists of multiplexed pools, index files are also provided. These can be used for sample deconvolution using the Cell Ranger pipeline. NOTE: Due to EGA's constrains in the number of files permitted per sample, index files (I1 and I2) for the HSPC data set are provided as a separate experiment named "FASTQ index files (I1 and I2) for deconvolution of 10X single-cell multiomics libraries".

  Dataset EGAD00001010927

• Dataset for soft_tissue_tumor-RNA

  Rare cancer sequencing data of 243 runs in tumor/control pairs, which were uploaded to umbrella studies. The sequencing was always paired

  Dataset EGAD00001008859

• EPICC: Evolutionary Predictions in Colorectal Cancer

  ATAC-seq profiling bam files from colorectal carcinoma and adenoma.

  Dataset EGAD00001008140

• Genetic drivers of epigenetic and transcriptional variation of human immune responses to infection (RNA-seq, ATAC-seq and ChIPmentation)

  RNA-seq, ATAC-seq and ChIPmentation data from monocyte-derived macrophages that were infected with Influenza A virus strain PR8WT, or a matching non-infected control.

  Dataset EGAD00001008422

• RNA-Seq in Patients with Primordial Dwarfism

  Transcriptome studies in patients with rare genetic diseases can potentially aid in the interpretation of likely causal genetic variation through identification of altered transcript abundance and/or structure. RNA-Seq is the most sensitive assay for both investigating transcript structure and abundance The primary aim of this pilot project is to investigate to what degree integrating exome-Seq and RNA-Seq data on the same individual can accelerate the identification of causal alleles for rare genetic diseases. There are two main strands to this: (i) identifying which variants discovered in exome-seq appear to be having a functional impact on transcripts, and (ii) identifying transcript outliers, especially among known causal genes, that may not necessarily have a causal variant identified from exome sequencing. The latter may identify the presence of causal variants that lie far from coding regions (e.g. the formation of cryptic splice sites deep within introns, or loss of long range regulatory elements), which can be confirmed with further targeted genetic assays. Just over 50% of all disease-causing variants recorded in the Human Gene Mutation Database (HGMD) affect transcript structure and abundance (e.g. nonsense SNVs, essential splice site SNVs, frameshifting indels, CNVs). This pilot project will study RNA from lymphoblastoid cell-lines from 12 patients with primordial dwarfism syndromes, for 10 of these samples we have previously generate exome data as part of our collaboration with the group of Prof Andrew Jackson. The two remaining samples are positive controls where the causal mutation is known, and is known to affect transcript structure and/or abundance. Primordial dwarfism is a prime candidate for these RNA-seq studies because all known causal mutations to date have key roles in DNA replication and thus, unsurprisingly, the products of the causal genes are typically ubiquitously expressed. Each RNA will be sequenced, with two technical replicates (independent RT-PCR and libraries) per sample, and each replicate run in 1/2 of a HiSeq lane using 100bp paired reads. Samples preparation was as follows :The cells were grown to confluency, then pellets frozen at -80. RNA samples were prepared using the Qiagen RNeasy kit, then nanodropped and analyzed using the bioanalyzer to determine concentration and purity. 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/

  Dataset EGAD00001000640

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__clear_cell_renal_cell_carcinoma

  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. Somatic mutational signatures are subsequently extracted by non-negative matrix factorisation methods and correlated with risk factors data. Through an enhanced understanding of cancer aetiology, Mutographs unprecedented effort is anticipated to outline modifiable risk factors, lead to new approaches to prevent cancer, and provide opportunities to empower early detection, refine high-risk groups and contribute to further therapeutic development.

  Study EGAS00001002608

• Diverse_outcomes_of_controlled_human_malaria_infection_originate_from_host_intrinsic_immune_variation_and_not_var_gene_switching

  Falciparum malaria is clinically heterogeneous and the relative contribution of parasite and host factors in shaping disease severity remains unclear. We set out to explore the interaction between host inflammation and parasite variant surface antigen (VSA) expression, asking whether this relationship underpins the variation observed in controlled human malaria infection. We uncovered marked heterogeneity in the response of naive hosts to blood challenge with some volunteers maintaining a state of immune quiescence, others triggering interferon-stimulated inflammation and a small group showing transcriptional evidence of myeloid cell suppression. Significantly, only inflammatory volunteers experienced hallmark symptoms of clinical malaria. When we then tracked temporal changes in parasite VSA expression to ask whether variants associated with severe disease preferentially expand in naive hosts (as predicted by current theory) we found that var gene profiles were essentially unchanged after 10-days of blood-stage infection. These data therefore show that the diverse clinical outcomes of CHMI largely originate from host-intrinsic variation and there is no evidence for switching or selection of var genes in naïve hosts. 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 EGAS00001003766

• RNA-Seq of vocal fold fibroblasts in Reinke’s edema and control subjects

  The voice disorder Reinke’s edema (RE) is a smoking- and voice-abuse associated benign lesion of the vocal folds, defined by an edema of the Reinke space, accompanied with pathological microvasculature changes and immune cell infiltration. Vocal fold fibroblasts (VFF) are the main cell type of the lamina propria and play a key role in the disease progression. Current therapy is restricted to symptomatic treatment. Hence, there is an urgent need for a better understanding of the molecular causes of the disease. In the present study, we investigated differential expression profiles of RE and control VFF by means of RNA sequencing. In addition, fast gene set enrichment analysis (FGSEA) was performed in order to obtain involved biological processes, mRNA and protein levels of targets of interest were further evaluated. We identified 74 differentially regulated genes in total, 19 of which were upregulated and 55 downregulated. Differential expression analysis and FGSEA revealed upregulated genes and pathways involved in extracellular matrix (ECM) remodeling, inflammation and fibrosis. Downregulated genes and pathways were involved in ECM degradation, cell cycle control and proliferation. The current study addressed for the first time a direct comparison of VFF from RE to control and evaluated immediate functional consequences.

  Study EGAS00001005130

• Identification of biomarkers of response to preoperative talazoparib monotherapy in treatment naïve gBRCA+ breast cancers

  Germline mutations in BRCA1 or BRCA2 exist in ~2-7% of breast cancer patients, which has led to the approval of PARP inhibitors in the advanced setting. We have previously reported a phase II neoadjuvant trial of single agent talazoparib for patients with germline BRCA pathogenic variants with a pathologic complete response (pCR) rate of 53%. As nearly half of the patients treated did not have pCR, better strategies are needed to overcome treatment resistance. To this end, we conducted multi-omic analysis of 13 treatment naïve breast cancer tumors from patients that went on to receive single-agent neoadjuvant talazoparib. We looked for biomarkers that were predictive of response (assessed by residual cancer burden) after 6 months of therapy. We found that all resistant tumors exhibited either the loss of SHLD2, expression of a hypoxia signature, or expression of a stem cell signature. These results indicate that the deep analysis of pre-treatment tumors can identify biomarkers that are predictive of response to talazoparib and potentially other PARP inhibitors, and provides a framework that will allow for better selection of patients for treatment, as well as a roadmap for the development of novel combination therapies to prevent emergence of resistance.

  Study EGAS00001005676

• Presence of rare potential pathogenic variants in subjects under 65 years old with very severe or fatal COVID‑19

  Rare variants affecting host defense against pathogens could be involved in COVID‑19 severity and may help explain fatal outcomes in young and middle‑aged patients. Our aim was to report the presence of rare genetic variants in certain genes, by using whole exome sequencing, in a selected group of COVID‑19 patients under 65 years who required intubation or resulting in death (n = 44). To this end, different etiopathogenic mechanisms were explored using gene prioritization‑based analysis in which genes involved in immune response, immunodeficiencies or blood coagulation were studied. We detected 44 different variants of interest, in 29 different patients (66%). Some of these variants were previously described as pathogenic and were located in genes mainly involved in immune response. A network analysis, including the 42 genes with candidate variants, showed three main components, consisting of 25 highly interconnected genes related to immune response and two additional networks composed by genes enriched in carbohydrate metabolism and in DNA metabolism and repair processes. In conclusion, we have detected candidate variants that may potentially influence COVID‑19 outcome in our cohort of patients. Further studies are needed to confirm the ultimate role of the genetic variants described in the present study on COVID‑19 severity.

  Study EGAS00001006372

• Multi-organ landscape of therapy-resistant melanoma

  Metastasis and failure of present-day therapies represent the most common causes of mortality in patients with cutaneous melanoma. To identify the underlying genetic and transcriptomic landscapes, here we analyzed multi-organ metastases and tumor-adjacent tissues from eleven rapid autopsies after treatment with MAPK inhibitor (MAPKi) and/or immune checkpoint blockade (ICB) and death due to acquired resistance. Either treatment elicits shared genetic alterations that suggest immune-evasive, cross-therapy resistance mechanisms. Large, non-clustered deletions, inversions, and inter-chromosomal translocations dominate rearrangements. Using 345 therapy-naïve TCGA melanoma and 131 patient-matched melanoma pre-and-post acquired resistance to only either treatment, we performed cross-cohort analyses to identify MAPKi and ICB as respective contributors to gene amplifications and deletions enriched in autopsy versus therapy-naïve tumors. Private/late mutations and structural variants display shifted mutational and rearrangement signatures, with MAPKi specifically selecting for signatures of defective homologous-recombination, mismatch, and base-excision repair. Transcriptomic signatures and cross-talks with tumor-adjacent macroenvironments nominate organ-specific adaptive pathways. An immune-desert, CD8+-macrophage-biased archetype, T-cell exhaustion, and type-2 immunity characterize the immune contexture. This multi-organ analysis of lethal, therapy-resistant melanoma presents preliminary insights with potential to improve therapeutic strategies.

  Study EGAS00001006644

• Clinical outcomes and immune correlates of response to nivolumab plus chemoradiotherapy in women with locally-advanced cervical cancer – NiCOL study

  Concurrent chemoradiotherapy (CRT) with blockade of the PD-1 pathway may enhance immune-mediated tumor control through increased phagocytosis, cell death, and antigen presentation. The NiCOL phase 1 trial was designed to determine safety, tolerance, and immune correlates of response to the PD-1 blocking antibody, nivolumab, with and following concurrent CRT in 16 women with locally advanced cervical cancer (LACC). Three patients experienced grade 3 dose-limiting toxicities. The overall response rate was 93.8% [95%CI: 69.8%-99.8%] with a 2-year PFS of 75% [95% CI: 56.5%-99.5%]. Compared to patients with progressive disease (PD), progression-free (PF) subjects showed a brisker stromal immune infiltrate, higher proximity of tumor-infiltrating CD3+ T cells to PD-L1+ tumor cells and of Foxp3+ T cells to proliferating CD11c+ myeloid cells. PF had higher baseline levels of PD-1 on EMRA CD4+ T cells and ICOS-L on tumor-associated macrophages (TAMs) vs PD, which instead displayed enhanced PD-L1 expression on TAMs and higher, circulating frequencies of proliferating Tregs at baseline. Concomitant nivolumab plus definitive CRT is safe and associated with encouraging PFS rates. Further validation in the subset of LACC displaying pre-existing, adaptive immune activation is warranted.

  Study EGAS00001007297

• Transcriptome of CD4+ T cells and CD8+ T cells in glioblastoma multiforme

  Patients included in this study were over 18 years of age and had a histology-confirmed diagnosis of glioblastoma multiforme (GBM). Exclusion criteria were the previous administration of any anti-tumor therapy including radiation therapy. All patients gave written informed consent. The study was approved by the local ethics committee (TUM Medical school) and conducted following the Declaration of Helsinki. During resection of the tumors, tumor tissue and tissue from normal appearing brain within the operative channel was collected. Blood was drawn during the surgical procedure. Single cell suspensions were prepared from the tumor tissue, the normal appearing brain, and the blood. CD4+ T cells and CD8+ T cells were sorted by flow cytometry. Only patients with a complete set of specimens (CD4+ tumor infiltrating lymphocytes (TIL), CD8+ TIL, CD4+ T cells from normal appearing brain, CD8+ T cells from normal appearing brain, blood-derived CD4+ and CD8+ T cells) containing a minimum of 1000 cells in each sorted sample were further analyzed (n=9). Total RNA was isolated from sorted cell populations using the RNAeasy Plus micro kit (Qiagen, 74034). Quality and integrity of total RNA was controlled on a Bioanalyzer 2100 (Agilent Technologies). Library preparation for bulk-sequencing of poly(A)-RNA was done as described previously (Parekh et al., 2016). Briefly, barcoded cDNA of each sample was generated with a Maxima RT polymerase (ThermoFisher Scientific, EP0742) using oligo-dT primer containing barcodes, unique molecular identifiers (UMIs) and an adaptor. Ends of the cDNAs were extended by a template switch oligo (TSO) and full-length cDNA was amplified with primers binding to the TSO-site and the adaptor. NEB UltraII FS kit was used to fragment cDNA. After end repair and A-tailing, a TruSeq adapter was ligated and 3'-end-fragments were finally amplified using primers with Illumina P5 and P7 overhangs. In comparison to previous descriptions (Parekh et al., 2016), the P5 and P7 sites were exchanged to allow sequencing of the cDNA in read 1 and barcodes and UMIs in read 2 to achieve a better cluster recognition. The library was sequenced on a NextSeq 500 (Illumina) with 59 cycles for the cDNA in read 1 and 16 cycles for the barcodes and UMIs in read 2. Data were processed using the published Drop-seq pipeline (v1.0) to generate sample- and gene-wise UMI tables (Macosko et al., 2015). Reference genome (GRCh38) was used for alignment. Transcript and gene definitions were used according to the Genecode Annotation Version 35.

  Dataset EGAD50000000226

• Liquid biopsy for molecular characterization of diffuse large B-cell lymphoma and early assessment of minimal residual disease

  Circulating tumor DNA (ctDNA) allows genotyping and minimal residual disease (MRD) detection in lymphomas. Using a NGS approach (Euroclonality-NDC), we evaluated clinical and prognostic value of ctDNA in a series of R-CHOP-treated DLBCL patients at baseline (n=68) and after 2-cycles (n=59), monitored by metabolic imaging (PET/CT). A molecular marker was identified in 61/68 (90%) ctDNA samples at diagnosis. Pre-treatment high ctDNA levels significantly correlated with elevated LDH, advanced stage, high risk IPI and a trend to shorter 2-year PFS. Valuable NGS data after 2-cycles of treatment were obtained in 44 cases, and 38 achieved major molecular response (MMR; 2.5-log drop in ctDNA). PFS curves displayed statistically significant differences among those achieving MMR vs. those not achieving MMR (2yr PFS of 76% vs. 0%, p<0.001). Similarly, more than 66% reduction in SUVmax by PET/CT identified two subgroups with different prognosis (2yr PFS of 83% vs. 38%; p<0.001). Combining both approaches MMR and SUVmax reduction, a better stratification was observed (2yr PFS of 84% vs. 17% vs. 0%, p<0.001). Euroclonality-NDC panel allows the detection of a molecular marker in the ctDNA in 90% of DLBCL. ctDNA reduction at 2 cycles and its combination with interim PET results improves patient prognosis stratification.

  Study EGAS50000000215

• Oncogenic Ectodomain deletion of FGFR1 caused by Breakage-Fusion-Bridges in Squamous Cell lung Cancer

  The identification of recurrent 8p11.23 amplifications including FGFR1 raised the hope of a treatable target in squamous cell lung cancer (SQLC). However, only a minority of patients with FGFR1-amplified tumors respond to single agent inhibitor therapy targeting FGFR. To understand the underlying mechanism of FGFR1 dependency, we performed whole genome and transcriptome sequencing of 25 FGFR1-amplified primary tumors with unknown response upon FGFR inhibition. In addition, we performed deep sequencing of 26 FGFR1-amplified samples whereof the response upon FGFR inhibition was known for 25 samples. In both cohorts we identified intra-chromosomal tail-to-tail breaks close to the FGFR1 transcription start site, being responsible for focal amplification of FGFR1. These specific breaks are caused by a Breakage-Fusion-Bridge-like (BFB-like) mechanism. Here, we associate these breaks with FGFR inhibitor sensitivity. Moreover, in some cases these breaks are located within the open reading frame of FGFR1, which leads to the expression of an ΔEC-FGFR1 transcript that lacks the ecto-domain. Overexpression of ΔEC-FGFR1 transforms Baf3 cells and lead to an FGFR1-dependent phenotype. Our results demonstrate that the truncation of the FGFR1 ectodomain is a frequent event in 8p11.23-amplified squamous cell lung cancer caused by tail-to-tail breaks. These breaks might be used as a predictive therapeutic marker to stratify patients for FGFR-inhibitor therapy.

  Study EGAS00001005059

• Detection of rare mutations, copy number variation, and DNA methylation in the same template DNA molecules

  The analysis of cell-free DNA (cfDNA) from plasma offers great promise for the earlier detection of cancer. At present, changes in DNA sequence, methylation, or copy number are the most sensitive ways to detect the presence of cancer. To further increase the sensitivity of such assays with limited amounts of sample, it would be useful to be able to evaluate the same template molecules for all these changes. Here we report an approach, called MethylSaferSeqS, that achieves this goal, and can be applied to any standard library preparation method suitable for massively parallel sequencing. The innovative step was to copy both strands of each DNA-barcoded molecule with a primer that allows the subsequent separation of the original strands (retaining their 5-methylcytosine residues) from the copied strands (in which the 5-methylcytosine residues are replaced with unmodified cytosine residues). The epigenetic and genetic alterations present in the DNA molecules can then be obtained from the original and copied strands, respectively. We applied this approach to plasma from 265 individuals, including 198 with cancers of the pancreas, ovary, lung and colon, and found the expected patterns of mutations, copy number alterations, and methylation. Furthermore, we could determine which original template DNA molecules were methylated and/or mutated. MethylSaferSeqS should be useful for addressing a variety of questions relating genetics and epigenetics in the future.

  Study EGAS00001006839

• Knoll et al Identification of drug candidates targeting monocyte reprogramming in people living with HIV

  Introduction People living with HIV (PLHIV) are characterized by functional reprogramming of innate immune cells even after long-term antiretroviral therapy (ART). In order to assess technical feasibility of omics technologies for application to larger cohorts, we compared multiple omics layers. Methods Bulk and single-cell transcriptomics, flow cytometry, proteomics, chromatin landscape analysis by ATAC-seq as well as ex vivo drug stimulation was performed in a small number of blood samples derived from PLHIV and healthy controls from the 200-HIV cohort study. Results Single-cell RNA-seq analysis revealed that most immune cells in peripheral blood of PLHIV are altered in their transcriptomes and that a specific functional monocyte state previously described in acute HIV infection is still existing in PLHIV while other monocyte cell states are only occurring acute infection. Moreover, a reverse transcriptome approach on a rather small number of PLHIV was sufficient to identify drug candidates for reversing the transcriptional phenotype of monocytes in PLHIV. Discussion These scientific findings and technological advancements for clinical application of single-cell transcriptomics form the basis for the larger 2000-HIV multicenter cohort study on PLHIV, for which a combination of bulk and single-cell transcriptomics will be included as the leading technology to determine disease endotypes in PLHIV and to predict disease trajectories and outcomes.

  Study EGAS00001007460

• The neo-open reading frame peptides that comprise the tumor framome are a rich source of neoantigens for cancer immunotherapy

  Identification of immunogenic cancer neoantigens as targets for therapy is challenging. Here, we integrate whole genome and long-read transcript sequencing of cancers to identify the collection of neo-open reading frame peptides (NOPs) expressed in tumors. We termed this collection of NOPs the tumor framome. NOPs represent tumor-specific peptides that are different from wild-type proteins and may be strongly immunogenic. We describe a class of hidden NOPs that derive from structural genomic variants involving an upstream protein coding gene driving expression and translation of non-coding regions of the genome downstream of a rearrangement breakpoint, i.e. where no gene annotation or evidence for transcription exists. The entire collection of NOPs represents a vast number of possible neoantigens particularly in tumors with many structural genomic variants and a low number of missense mutations. We show that NOPs are immunogenic and epitopes derived from NOPs can bind to MHC class I molecules. Finally, we provide evidence for the presence of memory T cells specific for hidden NOPs in peripheral blood from a patient with lung cancer. This work highlights NOPs as a major source of possible neoantigens for personalized cancer immunotherapy and provides a rationale for analyzing the complete cancer genome and transcriptome as a basis for detection of NOPs.

  Study EGAS00001006021

• Personalized peptide vaccination among 173 patients with IDH wildtype glioblastoma: a retrospective analysis

  Tumors from 173 GBM patients were analysed for somatic mutations to generate a personalized peptide vaccine targeting tumor-specific neoantigens. Exome libraries for 173 glioblastoma tumors and matched normal DNA were sequenced on Illumina platform, alongside whole transcriptomes from the tumor samples. Paper Abstract: Current treatment outcome of patients with glioblastoma (GBM) remains poor. Following standard therapy, recurrence is universal with limited survival. GBM tumors from 173 patients were analysed for somatic mutations to generate a personalized peptide vaccine targeting tumor-specific neoantigens. Among all patients, including 70 treated prior to progression (primary) and 103 treated after progression (recurrent), the median overall survival from first diagnosis was 31.9 months (95% CI: 25.0-36.5). Side effects were infrequent and were predominantly grade 1 or 2. A vaccine-induced immune response to at least one of the vaccinated peptides was detected in blood samples of 87 of 99 (88%) monitored patients. T-cell responses to vaccinated neoepitope peptides were durable in most patients. Significantly prolonged survival was observed for patients with multiple vaccine-induced immune responses (53 months) compared to those with no/low induced responses (27 months; P=0.03). Altogether, our results highlight that the application of personalized neoantigen-targeting peptide vaccine is feasible and represents a promising potential treatment option for GBM patients.

  Study EGAS50000000449

• Mid-pass Whole-genome Sequencing in a Malagasy Cohort Uncovers Body Composition Associations

  The majority of human genomic research studies have been conducted in European-ancestry cohorts, reducing the likelihood of detecting potentially novel and globally impactful findings. Here, we present mid-pass whole-genome sequencing data and a genome-wide association study in a cohort of 264 self-reported Malagasy individuals from three locations on the island of Madagascar. We describe genetic variation in this Malagasy cohort, providing insight into the shared and unique patterns of genetic variation across the island. We observe phenotypic variation by location, and find high rates of hypertension particularly in the Southern Highlands sampling site as well as elevated self-reported malaria prevalence in the West Coast site relative to other sites. After filtering to a subset of 214 minimally-related individuals, we find a number of genetic associations with body composition traits, including many variants that are only observed in African populations or populations with admixed African ancestry from the 1000 Genomes Project. This study highlights the importance of including diverse populations in genomic research for the potential to gain novel insights, even with small cohort sizes. This project was conducted in partnership and consultation with local stakeholders in Madagascar and serves as an example of genomic research that prioritizes community engagement and that has potential impacts on our understanding of human health and disease.

  Study EGAS50000000496

• Identification and characterisation of pathogenic and non-pathogenic FGF14 repeat expansions

  Repeat expansions in FGF14 cause autosomal dominant late-onset cerebellar ataxia (SCA27B) with estimated pathogenic thresholds of 250 (incomplete penetrance) and 300 AAG repeats (full penetrance), but the sequence of pathogenic and non-pathogenic expansions remains unexplored. Here, we demonstrate that STRling and ExpansionHunter accurately detect FGF14 expansions from short-read genomic data using outlier approaches. By combining long-range PCR and nanopore sequencing in 169 patients with cerebellar ataxia and 802 controls, we compare FGF14 expansion alleles, including interruptions and flanking regions. Uninterrupted AAG expansions are significantly enriched in patients with ataxia from a lower threshold (180-200 repeats) than previously reported based on expansion size alone. Conversely, AAGGAG hexameric expansions are equally frequent in patients and controls. Distinct 5' flanking regions, interruptions and pre-repeat sequences correlate with repeat size. Furthermore, pure AAG (pathogenic) and AAGGAG (non-pathogenic) repeats form different secondary structures. Regardless of expansion size, SCA27B is a recognizable clinical entity characterized by frequent episodic ataxia and downbeat nystagmus, similar to the presentation observed in a family with a previously unreported nonsense variant (SCA27A). Overall, this study suggests that SCA27B is a major overlooked cause of adult-onset ataxia, accounting for 23-31% of unsolved patients. We strongly recommend re-evaluating pathogenic thresholds and integrating expansion sequencing into the molecular diagnostic process.

  Study EGAS50000000481

• Raw sequencing files from WGS and RNA-Seq

  Acute leukemias are the most common malignancies in children with an onset of the disease among 0-6 years-old. The early development of theses malignancies has been proposed as an indicator of their prenatal origin, hypothesis that has been confirmed especially in the case of acute lymphoblastic leukemia. For acute myeloid leukemia (AML), comprising an heterogenous group of different cytogenetics subtypes, this prenatal origin has also been confirmed specially for MLL-rearranged and t(8;21) subgroups. AMLs carrying t(7;12) constitute up to 30% of cases under 2-years-old, frequently being associated with a trisomy of chromosome 19 and an overexpression of MNX1 gene. An in utero origin of this specific subgroup has been suggested but not clearly confirmed. Here, we present the case of a 5-months-old AML patient with t(7;12) in which we have fully characterized a novel breakpoint occurring among ETV6 and NOM1 genes. We demonstrate the presence of this specific ETV6-NOM1 fusion in CD34+ hematopoietic stem and progenitor cells isolated from the cord blood unit from the same patient and stored at birth. Moreover, trisomy +19 present at diagnosis was not detected at birth. Together, our results demonstrate the prenatal origin of the t(7;12) in infant AML, and establish trisomy +19 as a secondary event in this leukemia.

  Study EGAS50000000186

• A comprehensive DNA methylation landscape of human and mouse cell lines derived from hematological malignancies

  Human cancer cell lines constitute useful models to study the primary disease and many relevant findings in tumor biology have been originated from them since the first immortalized cell line (HeLa) was obtained. The easy experimental intervention, the purity of the transformed cells and their versatility to undergo high-throughput screenings represent advantageous features of the established cancer cell lines. In recent years, major efforts have characterized in detail the multiomics make-up of hundreds of cancer cell lines and studied their association with sensitivity to anticancer drugs. Most of these studies have been genetic-centric and have not characterized in detail the epigenetic profiles underlying the characteristics of these cells or their impact on the efficacy of antitumoral compounds. In this regard, a past version of a DNA methylation microarray or a less accessible and more time-consuming such as reduced representation bisulfite sequencing have been used in those attempts to interrogate the epigenetic setting. Herein, we have obtained the DNA methylation profiles of 210 cell lines derived from hematological malignancies utilizing comprehensive DNA methylation microarrays that interrogates more than 850,000 and 285,000 CpG sites from human and mouse genomes, respectively. Importantly, we also provide a pharmacoepigenetic example of the potential use of this resource by showing how DNA methylation can predict response to nucleoside analogues.

  Study EGAS50000000627

• DUX4 activates novel intergenic transcripts and isoforms in a tissue-specific context

  DUX4, a transcription factor driving zygotic genome activation during the embryonic cleavage stage in mammals, regulates the expression of genic and non-genic elements and modulates chromatin accessibility. Its misexpression in skeletal muscle causes Facioscapulohumeral muscular dystrophy (FSHD). The sporadic expression of DUX4 in FSHD muscle and the limitations of short-read RNA sequencing challenge, however, a comprehensive grasp of its consequences on the muscle transcriptome. Here, leveraging full-length isoform RNA sequencing with short-read RNA sequencing on DUX4-inducible myoblasts, we elucidated DUX4-induced transcriptomic changes within a myogenic context. DUX4 yields a complex transcriptome featuring numerous unannotated isoforms from known and novel loci. Moreover, we observed distinct isoform usage of DUX4 target genes between myogenic and embryonic cleavage states highlighting the distinct differences between DUX4 transcription factor activity in skeletal muscle and during the cleavage stage. We also identified ~1,400 unannotated intergenic isoforms predominantly originating from repetitive elements such as LTRs, LINEs, and SINEs. Integration with available datasets delineates the transcriptional and epigenetic landscape of these loci across myogenic and embryonic backgrounds. Like the cell type-specific isoform production, the usage of DUX4 binding sites at these intergenic loci is in part cell type-specific. Collectively, these findings demonstrate substantial and context-specific transcriptomic alterations induced by DUX4 further enriching our understanding of DUX4-induced muscle pathology.

  Study EGAS50000000503

• Localised colon cancer WES study contaning WBCs, tissue and plasma samples at different time points

  Decision-making based on circulating tumor DNA (ctDNA) status following curative-intent surgery in stage II-III colon cancer (CC) is where the strongest clinical evidence lies. Commercial or academic assays are predominantly used to detect ctDNA, aiming to reveal minimal residual disease (MRD) and identify patients at high risk of recurrence. However, while current tests demonstrate exceptional specificity in detecting MRD among patients who later experience recurrence, their sensitivity is notably limited, particularly immediately after surgery when decisions regarding ACT are crucial. Additionally, up to 70% of patients with stage II-III colon cancer fail to clear ctDNA after standard adjuvant treatment (ACT) and subsequently experience recurrence. In this study, we performed whole exome sequencing (WES) of ctDNA at different timepoints in relapsed CC patients using two independent cohorts, alongside transcriptomic and proteomic analyses of metastases to gain a deeper understanding of progression mechanisms. An assay based on a tumor-agnostic approach from WES could enhance sensitivity in detecting MRD disease compared to current assays. Immune evasion appears to be the primary driver of progression in the localized CC setting, facilitated by a functional mutational burden at relapse. This suggests that immunotherapy could extend its efficacy to patients with microsatellite stability, broadening treatment options. Moreover, organoid modeling affirms the promising potential of targeted therapy to surpass conventional treatments in eradicating MRD.

  Study EGAS50000000204

• A thymic ILC1-like progenitor with differentiation potential towards KIR+NKG2A- NK cells

  KIR+NKG2A- Natural Killer (NK) cells have the unique ability to detect downregulation of single HLA-I allotypes, frequently occurring in malignantly-transformed and virus-infected cells in order to escape T cell recognition. We have recently shown that circulating Innate Lymphoid Cells 1 (cILC1) have the potential to generate such KIR+NKG2A- NK cells but their developmental origin was unknown. Here we demonstrate that the development of cILC1 is thymus-dependent and identify a putative progenitor of cILC1s in the thymus, called thyILC1. Single-cell RNAseq analysis revealed a close relationship of thyILC1s to CD34+ double-negative (DN) thymocytes, but in contrast to DNs, thyILC1s had lost T cell differentiation potential. However, like cILC1s they could be efficiently differentiated into KIR+NKG2A- NK cells. thyILC1s strongly expressed RAG1/2, PTCRA (preTCR, and NOTCH1/3, but lacked expression of TCR chains, suggesting that they failed to pass the -selection checkpoint of thymocyte development. Finally, analyses of patients with FOXN1het deficiency, showing thymic hypoplasia due to lack of mature thymic epithelial cells, exhibited a profound deficiency of cILC1s but not cILC2s and cILC3s, demonstrating thymus dependency of cILC1s. Together, the data suggest that thyILC1s are the source of a novel thymus-dependent NK cell differentiation pathway that promotes generation of KIR+NKG2A- NK cells.

  Study EGAS50000000760

• A thymic ILC1-like progenitor with differentiation potential towards KIR+NKG2A- NK cells - sc

  KIR+NKG2A- Natural Killer (NK) cells have the unique ability to detect downregulation of single HLA-I allotypes, frequently occurring in malignantly-transformed and virus-infected cells in order to escape T cell recognition. We have recently shown that circulating Innate Lymphoid Cells 1 (cILC1) have the potential to generate such KIR+NKG2A- NK cells but their developmental origin was unknown. Here we demonstrate that the development of cILC1 is thymus-dependent and identify a putative progenitor of cILC1s in the thymus, called thyILC1. Single-cell RNAseq analysis revealed a close relationship of thyILC1s to CD34+ double-negative (DN) thymocytes, but in contrast to DNs, thyILC1s had lost T cell differentiation potential. However, like cILC1s they could be efficiently differentiated into KIR+NKG2A- NK cells. thyILC1s strongly expressed RAG1/2, PTCRA (preTCR, and NOTCH1/3, but lacked expression of TCR chains, suggesting that they failed to pass the -selection checkpoint of thymocyte development. Finally, analyses of patients with FOXN1het deficiency, showing thymic hypoplasia due to lack of mature thymic epithelial cells, exhibited a profound deficiency of cILC1s but not cILC2s and cILC3s, demonstrating thymus dependency of cILC1s. Together, the data suggest that thyILC1s are the source of a novel thymus-dependent NK cell differentiation pathway that promotes generation of KIR+NKG2A- NK cells.

  Study EGAS50000000790

• Rituximab-treated lymphoma patients show correlated deficiency in serological and T cell Spike-specific response after SARS-CoV-2 vaccination: insights from the CORSA Study.

  This study is part of the CORSA trial (NCT04345315), a prospective multicentric research project conducted by the Istituto Romagnolo per lo Studio dei Tumori “Dino Amadori” to investigate the immunologic response to mRNA SARS-CoV-2 vaccines in patients with hematologic malignancies (HM) and healthy individuals as controls (HC). The data in this repository refer to two subcohorts: 7 healthy controls and 10 lymphoma patients, randomly selected. Whole blood samples were collected to characterize T cell receptor (TCR) repertoire specificity against SARS-CoV-2 proteins before and after two BNT162b2 vaccination cycles, at baseline (before vaccination), 50 days after the first dose (V1), and 6 months after the first dose (V3). All participants were vaccinated under the Italian national campaign and provided informed consent following the Declaration of Helsinki. We investigated SARS-CoV-2-specific T cell responses through TCR repertoire profiling and examined correlations between serological and cellular responses in both patient and control groups. Using Next-Generation Sequencing (NGS) to characterize the CDR3 sequences of TCR β-chains from peripheral blood, we assessed clonal expansion (depth) and diversification (breadth) of SARS-CoV-2-specific T cell clonotypes over time, aiming to understand how cancer and its treatment influence vaccine-induced immunity in immunocompromised individuals.

  Study EGAS50000001205

• MRCA and MRCE SNP genotypes

  Gene expression levels can be an important link DNA between variation and phenotypic manifestations. Our previous map of global gene expression based on approximately 400K SNPs and 50K transcripts in 400 sib pairs from the MRCA family panel has been widely used to interpret the results of GWAS. Here, we more than double the size of our initial dataset with expression data on 550 additional individuals from the MRCE family panel using the Illumina whole genome expression array. We have used new statistical methods for dimension reduction to account for non-genetic effects in estimates of expression levels and we have also included SNPs imputed from the 1000 Genomes Project. Our methods reduced false discovery rates and increased the number of eQTLs mapped either locally or at a distance (i.e. in cis or trans) from 1,534 in the MRCA dataset to 4,452 (with <5% FDR). Imputation of 1000 Genomes SNPs further increased the number of eQTLs to 7,302. Using the same methods and imputed SNPs in the newly acquired MRCE dataset we identified eQTLs for 9000 genes. The combined results identify strong local and distant effects for transcripts from 14,177 genes. Our eQTL database (http://www.hsph.harvard.edu/faculty/liming-liang/software/eqtl/) based on these results is freely available to help define the function of disease-associated variants.

  Study EGAS00000000137

• Organoid BulkRNAseq

  Here, we aim to understand early determinants of brain organoid quality and generated a set of brain organoids from 12 different hPSC lines using an adaptation of the unguided Lancaster protocol (Lancaster & Knoblich, 2014). All cell lines were obtained from healthy donors. 12 donor hPSC lines were differentiated into brain organoids and a set of 72 individual organoids (6 per line) were randomly selected morphologically diverse organoids and subjected to bulk RNA sequencing for gene expression analysis. RNA was extracted using the RNeasy kit (Qiagen, Germany) according to the manufacturer’s instructions. A total of 1000 ng per sample was sent for RNA sequencing to Azenta Life Sciences (Genewiz Leipzig, Germany) for sequencing. Library preparation was performed using the NEBNext Ultra II RNA Library Prep Kit for Illumina following manufacturer’s instructions (NEB, Ipswich, MA, USA). mRNAs were first enriched with Oligo(dT) beads. The samples were sequenced using a 2x150 Pair-End (PE) configuration v1.5 at a depth of >20 million reads in each sample. Raw sequence data (.bcl files) generated from Illumina NovaSeq6000 was converted into fastq files and de-multiplexed using Illumina bcl2fastq program version 2.20. The raw sequence reads were aligned to the human reference genome GRCh38 using the STAR aligner (Dobin et al., 2013). Gene expression levels were quantified using the FeatureCounts tool (Liao et al., 2014).

  Study EGAS50000000659

• Nala TAS-LRS PGx Study

  Pharmacogenomics (PGx) testing improves medication safety and efficacy by identifying genetic variants that affect drug response. However, current technologies often fail to resolve complex loci, detect structural variants, or phase alleles accurately. Here, we present an end-to-end PGx workflow based on Targeted Adaptive Sampling-Long Read Sequencing (TAS-LRS), integrating a streamlined laboratory protocol with a bioinformatics pipeline that includes a novel CYP2D6 caller. Using 1,000 ng of DNA and three-sample multiplexing on a single PromethION flow cell, the assay achieves consistent on-target (25x) and off-target (3x) coverage, enabling accurate, haplotype-resolved testing of 35 pharmacogenes alongside genome-wide genotyping from off-target reads. We further developed the workflow into a clinically ready service and validated its performance across 17 reference and clinical samples. The assay demonstrated high concordance for small variants (99.9%) and structural variants (>95%), with phased diplotypes and metabolizer phenotypes reaching 97.7% and 98.0% concordance, respectively. Improved calls were observed in 12 genes due to enhanced genotyping, phasing, or novel allele detection. In addition, off-target reads supported accurate genome-wide imputation, comparable to short-read sequencing and superior to microarrays. These results establish the feasibility of long-read sequencing for clinical PGx testing and position TAS-LRS as a scalable solution combining both targeted and genome-wide utility.

  Study EGAS50000001116

• Somatic rearrangements causing oncogenic ectodomain deletions of FGFR1 in squamous cell lung cancer

  The discovery of frequent 8p11-p12 amplifications in squamous cell lung cancer (SQLC) has fueled hopes that FGFR1, located inside this amplicon, might be a therapeutic target. In a clinical trial, only 11% of patients with 8p11 amplification (detected by FISH) responded to FGFR kinase inhibitor treatment. To understand the mechanism of FGFR1 dependency, we performed deep genomic characterization of 52 SQLCs with 8p11-p12 amplification, including 10 tumors obtained from patients who had been treated with FGFR inhibitors. We discovered somatically altered variants of FGFR1 with deletion of exons 1-8 that resulted from intragenic tail-to-tail rearrangements. These ectodomain-deficient FGFR1 variants (ΔEC-FGFR1) were expressed in the affected tumors and were tumorigenic in both in vitro and in vivo models of lung cancer. Mechanistically, breakage-fusion-bridges were the source of 8p11-p12 amplification, resulting from frequent head-to-head and tail-to-tail rearrangements. Generally, tail-to-tail rearrangements within or in close proximity upstream of FGFR1 were associated with FGFR1 dependency. Thus, the genomic events shaping the architecture of the 8p11-p12 amplicon provide a mechanistic explanation for the emergence of FGFR1-driven SQLC. Specifically, we believe that FGFR1 ectodomain-deficient and FGFR1-centered amplifications caused by tail-to-tail rearrangements are a novel somatic genomic event that might be predictive of therapeutically relevant FGFR1 dependency.

  Study EGAS50000001368

• Transposable Elements Shape Stemness in Normal and Leukemic Hematopoiesis

  Despite most acute myeloid leukemia (AML) patients achieving complete remission after induction chemotherapy, two-thirds will relapse with fatal disease within five years. AML is organized as a cellular hierarchy sustained by leukemia stem cells (LSC) at the apex, with LSC properties directly linked to tumor progression, therapy failure, and disease relapse 1–5. Despite the central role of LSC in poor patient outcomes, little is known about the genetic determinants driving their stemness properties. As LSCs share many functional and molecular properties with normal hematopoietic stem cells (HSC) 6, we investigated accessible chromatin unique across normal hematopoietic and cancer cell states and identified transposable elements (TEs) as genetic determinants of both primitive populations in comparison with their downstream mature progeny. A clinically-relevant TE chromatin accessibility-based LSCTE121 signature was developed that enabled patient classification based on survival outcomes. Through functional assays, primitive cell specific-TE subfamilies were found to serve as docking sites for stem cell-associated regulators of genome topology or lineage-specific transcription factors, including LYL1 in LSCs. Finally, using chromatin editing tools, we establish that chromatin accessibility at LTR12C elements in LSCs are necessary to maintain stemness properties. Our work identifies TEs as genetic drivers of primitive versus mature cell states, where distinct TE subfamilies account for stemness properties in normal versus leukemic hematopoietic stem cells.

  Study EGAS50000001141

• Transposable Elements Shape Stemness in Normal and Leukemic Hematopoiesis

  Despite most acute myeloid leukemia (AML) patients achieving complete remission after induction chemotherapy, two-thirds will relapse with fatal disease within five years. AML is organized as a cellular hierarchy sustained by leukemia stem cells (LSC) at the apex, with LSC properties directly linked to tumor progression, therapy failure, and disease relapse 1–5. Despite the central role of LSC in poor patient outcomes, little is known about the genetic determinants driving their stemness properties. As LSCs share many functional and molecular properties with normal hematopoietic stem cells (HSC) 6, we investigated accessible chromatin unique across normal hematopoietic and cancer cell states and identified transposable elements (TEs) as genetic determinants of both primitive populations in comparison with their downstream mature progeny. A clinically-relevant TE chromatin accessibility-based LSCTE121 signature was developed that enabled patient classification based on survival outcomes. Through functional assays, primitive cell specific-TE subfamilies were found to serve as docking sites for stem cell-associated regulators of genome topology or lineage-specific transcription factors, including LYL1 in LSCs. Finally, using chromatin editing tools, we establish that chromatin accessibility at LTR12C elements in LSCs are necessary to maintain stemness properties. Our work identifies TEs as genetic drivers of primitive versus mature cell states, where distinct TE subfamilies account for stemness properties in normal versus leukemic hematopoietic stem cells.

  Study EGAS50000001147

• Neoadjuvant immune checkpoint blockade in mismatch-repair proficient colon cancers

  Immune checkpoint blockade has led to paradigm shifts in the treatment of various tumour types, yet limited efficacy has been observed in patients with metastatic mismatch-repair proficient (pMMR) colorectal cancer. Here we report an in-depth analysis of patients with early-stage pMMR colon cancer from the phase II NICHE study (ClinicalTrials.gov: NCT03026140). A total of 31 patients received neoadjuvant treatment of nivolumab plus ipilimumab followed by surgery. The response rate was 26% and included six patients with a major pathological response (≤10% residual viable tumour). One patient with an ongoing clinical complete response did not undergo surgery. Responses were observed despite a low tumour mutational burden in all tumours, while chromosomal genomic instability scores were significantly higher in responders compared to non-responders. Responding tumours were more likely to be TP53 mutant, KRASG12 wildtype and had significantly higher baseline expression of proliferation signatures and TCF1. Imaging mass cytometry also revealed a higher percentage of Ki-67+ cancer and Ki-67+ CD8+ T cells in responders compared to non-responders. Meanwhile, in non-responders, we observed an enrichment for an inflamed-fibrotic microenvironment associated with TGF-β expression, hinting towards possible resistance mechanisms. These results provide a comprehensive analysis of response to neoadjuvant ICB in early-stage pMMR colon cancers and identify potential biomarkers for patient selection.

  Study EGAS50000000856

• Immune heterogeneity in small cell lung cancer and vulnerability to immune checkpoint blockade

  Following the pivotal Phase III IMpower133 study, atezolizumab (anti–PD-L1), combined with carboplatin and etoposide (CE), was the first immune checkpoint inhibitor approved for first-line treatment of extensive-stage small cell lung cancer (ES-SCLC) and is now a standard of care. A clearer understanding of therapeutically relevant SCLC subsets is needed to improve outcomes and to identify rational combination strategies. Transcriptomic analyses and nonnegative matrix factorization were conducted on 271 pre-treatment patient tumor samples from IMpower133 and four subsets with general concordance to previously reported SCLC subtypes were identified. Deeper investigation into the immune heterogeneity relevant to clinical outcomes uncover two subsets with differing neuroendocrine (NE) versus non-neuroendocrine (non-NE) phenotypes that demonstrated hallmarks of immune cell infiltration. The balance of tumor-associated macrophage (TAM) to T-effector signals distinguished these inflamed subsets. Tumors with low TAM but high T-effector signals had a NE phenotype and demonstrated longer overall survival with PD-L1 blockade and CE versus CE alone than did tumors with high TAM and high T-effector signal, which had a non-NE phenotype. The delineation of SCLC immune heterogeneity offers a clinically relevant approach to discriminate SCLC patients likely to benefit most from immunotherapeutic approaches and highlights the complex mechanisms underlying response and resistance to immune checkpoint blockade. 

  Study EGAS50000000138

• Phenotypic characterization and prognostic impact of CD103+ tissue-resident memory T cells in diffuse large B cell lymphoma

  Tissue-resident memory T (TRM) cells are a population of memory T cells that stably occupy tissues and play a key role in immunosurveillance, having been linked to favorable survival outcomes in various solid tumors. While TRM cells have been identified in lymph nodes, their phenotype and prognostic significance in B-cell non-Hodgkin lymphoma (B-NHL), including diffuse large B-cell lymphoma (DLBCL), remains poorly characterized. The frequency of CD103 expressing T cells in patient samples with DLCBL was quantified by immunofluorescence (IF) staining of tissue biopsies (n=306) and flow cytometry (n=252) of cell disaggregates, and linked with clinical outcome. The phenotype of TRM cells was characterized by spectral flow cytometry (n=10 DLBCL and 2 reactive lymph node (rLN) samples) and single-cell RNA sequencing (scRNAseq) (n=12 aggressive B cell lymphomas and 4 rLN samples). An additional 51 samples from an external dataset were included to verify the single cell transcriptome findings. The flow cytometry and scRNAseq findings in DLBCL were extrapolated to additional B-NHL entities. Through analysis of these datasets we were able to characterize CD103+ TRM-like cells in DLBCL and other B-NHL entities and identified them to represent a prognostically favorable population with an activated/cytotoxic T cell phenotype.

  Study EGAS50000000943

• Investigating Host Genetic Risk Factors for Tuberculosis in Highly Endemic South African Populations

  It is apparent from our recent population genetic and admixture mapping work 1–3 that a substantial part of the ancestral input into the South African Coloured (SAC) population is from the San or Khoe groups, which are not well represented in publicly available genetic databases. It is a reasonable assumption that the SAC population in the Western Cape may have derived more genetic input from the present KhoeSan population in the Northern Cape than the San residing in Namibia. It was concluded by our group that it is indeed the southern African KhoeSan group, the ≠Khomani, that best represent the KhoeSan contribution seen in the SAC 4. More recently, it has been shown that there are ancestry related increases in TB susceptibility, especially with increased Bantu-speaking African and KhoeSan ancestry 5. We aim to elucidate the epidemiological and human-host genetic risk factors for TB and the immunological pathways modulating TB infection. The proposed sub-study will enroll people evaluated for TB at Northern Cape community health clinics and their contacts from their households and/or community. We will conduct a demographic interview, collect saliva, which will be used for genetic analysis, and blood to determine latent TB infection status and capture immunological responses to mycobacterial infections at the cellular level and RNA sequencing.

  Study EGAS00001007850

• Osteosarcoma 3D patient derived cultures to test genome-informed personalized treatment options; a feasibility study

  Background: Improving osteosarcoma treatment beyond conventional (neo)adjuvant chemotherapy and resection remains challenging. An urgent need for novel therapeutic options, particularly personalized and targeted approaches, has emerged due to high inter-patient molecular heterogeneity. A lack of representative in vitro and in vivo models impedes therapeutic development, therefore we aimed to create 3D in vitro long term culture models directly from patient material. Methods: Tumour cells from seven osteosarcoma patients were propagated in monolayer or collagen hydrogels, while whole-exome sequencing of corresponding primary tumour tissue was performed to identify potential drug targets. Established cultures were subsequently used to assess efficacy of the identified personalized treatment options. Results: Three out of seven hydrogel cultures harboured the same genetic alterations as the corresponding primary tumours. Among these, the L6565 patient line exhibited homozygous CDKN2A loss with retained Rb expression, rendering tumour cells sensitive to CDK4/CDK6 inhibition via palbociclib. This patient line was used to evaluate genome-informed therapy and to compare cell culture models of increasing complexity. Our findings demonstrate feasibility of establishing long-term patient-derived osteosarcoma cultures with a success rate of 14%. Tumour heterogeneity was reflected in advanced culture methods producing more variability in treatment response. Conclusion: These results highlight the potential of genome-informed therapies in osteosarcoma and the importance of refining culture techniques to enhance translational research and therapeutic outcomes.

  Study EGAS50000001583

• UK10K NEURO ASD MGAS

  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 MGAS (Molecular Genetics of Autism Study) samples are from a clinical sample seen by specialists at the Maudsley hospital and who have had detailed phenotypic assessments with ADI-R and ADOS.For further information on this cohort please contact Patrick Bolton (patrick.bolton@kcl.ac.uk).

  Study EGAS00001000113