23014 results for "ega"

in 63.68 milliseconds.

• Genome of the Netherlands

  The Genome of the Netherlands (GoNL) Project characterizes DNA sequence variation, common and rare, for SNVs, short insertions and deletions (indels) and larger deletions in 769 individuals of Dutch ancestry selected from five biobanks under the auspices of the Dutch hub of the Biobanking and Biomolecular Research Infrastructure (BBMRI-NL). The samples come from a representative sample of 250 trio-families from all provinces in the Netherlands. The parent-offspring trios include adult individuals ranging in age from 19 to 87 years (mean=53 years; SD=16 years) from birth cohorts 1910-1994. Sequencing was done on blood-derived DNA from uncultured cells and accomplished coverage was 14-15x. Samples where contributed by LifeLines (http://lifelines.nl/lifelines-research/general), The Leiden Longevity Study (http://www.healthy-ageing.nl; http://www.langleven.net), The Netherlands Twin Registry (NTR: http://www.tweelingenregister. org), The Rotterdam studies, (http://www.erasmus-epidemiology.nl/rotterdamstudy) and the Genetic Research in Isolated Populations program (http://www.epib.nl/research/geneticepi/research.html#gip). The sequencing was carried out in collaboration with the Beijing Institute for Genomics (BGI). The analysis was done by a consortium lead by UMCG, LUMC, Erasmus MC, VU university and UMCU, see http://www.nlgenome.nl. Funding for the project was provided by the Netherlands Organization for Scientific Research under award number 184021007, dated July 9, 2009 and made available as a Rainbow Project of the Biobanking and Biomolecular Research Infrastructure Netherlands (BBMRI-NL).

  Study EGAS00001000644

• Cisplatin increases sensitivity to FGFR inhibition in patient-derived xenograft models of lung squamous cell carcinoma

  Lung squamous cell carcinoma (SqCC) is a molecularly complex and genomically unstable disease. No targeted therapy is currently approved for lung SqCC, although potential oncogenic drivers of SqCC have been identified, including amplification of the fibroblast growth factor receptor 1 (FGFR1). Reports from a recently completed clinical trial indicate low response rates in patients treated with FGFR tyrosine kinase inhibitors, suggesting inadequacy of FGFR1 amplification as a biomarker of response, or the need for combination treatment. We aimed to develop accurate models of lung SqCC and determine improved targeted therapies for these tumors. We show that detection of FGFR1 mRNA by RNA in situ hybridisation is a better predictor of response to FGFR inhibition than FGFR1 gene amplification using clinically relevant patient-derived xenograft (PDX) models of lung SqCC. FGFR1-overexpressing tumors were observed in all histological subtypes of non-small cell lung cancers (NSCLC) as assessed on a tissue microarray, indicating a broader range of tumors that may respond to FGFR inhibitors. In FGFR1-overexpressing PDX tumors, we observed increased differentiation and reduced proliferation following FGFR inhibition. Combination therapy with cisplatin was able to increase tumor cell death, and dramatically prolonged animal survival compared to single agent treatment. Our data suggest that FGFR tyrosine kinase inhibitors can benefit NSCLC patients with FGFR1-overexpressing tumors and provides a rationale for clinical trials combining cisplatin with FGFR inhibitors.

  Study EGAS00001002423

• Identification_of_synthetic_lethal_genes_by_CRISPR_Cas9_library

  Sequencing is performed on PCR products from the amplification of gRNA segments introduced into cell lines (some are commercial, some are derived by collaborators, as described in HMDMC 13/076). Cellular DNA or RNA are not sequenced. This data is defined as "Other Biological/Biochemical Assay Data" under the WTSI datasharing guidelines and consequently will be released at the time of publication and/or through faculty websites. SEQUENCING DETAILS: in October 2016 3 DNA samples will be submitted for PCR amplification and purification through CGP pipeline using Thompson's custom primer set and PCR settings and sequencing with Hiseq using a custom protocol (with 3 forward reads to read sgRNA1, sgRNA2 and the index, optimized with Illumina Bespoke Sanger team); read length 50bp. The following set of samples will be submitted in later 2016 and ealry 2017 as cell pellet to be processed through CGP Cas9 pipeline (Thompson's oligos and protocols) and sequences on Hiseq (same custom sequencing protocol described above) We are utilizing a custom CRISPR/Cas9 library to identify synthetic lethal gene pairs. Cells will be transduced with LV expressing Cas9, then infected with LV expressign the library of 42'000 sgRNA pairs. Cell pellet will be collected from the cell population at different timepoints. The samples will be processed by CGP pipeline customized for the specific protocol and then sequenced to identify the 19bp of the sgRNA and quantify their representation.

  Study EGAS00001002117

• The evolutionary history of human colitis-associated colorectal cancer

  Inflammatory bowel disease confers an increased lifetime risk of developing colorectal cancer (CRC), and colitis-associated CRC (CA-CRC) is molecularly distinct from sporadic CRC (S-CRC). Here we have dissected the evolutionary history of CA-CRC using multi-region sequencing. Exome sequencing was performed on fresh-frozen multiple regions of carcinoma, adjacent non-cancerous mucosa and blood from 12 CA-CRC patients (n=55 exomes), and key variants were validated with orthogonal methods. Genome-wide copy number profiling was performed using SNP arrays and low-pass whole genome sequencing on archival non-dysplastic mucosa (n=9), low-grade dysplasia (LGD, n=30), high-grade dysplasia (HGD, n=13), mixed LGD/HGD (n=7) and CA-CRC (n=19). Phylogenetic trees were reconstructed, and evolutionary analysis used to reveal the temporal sequence of events leading to CA-CRC. 10/12 tumors were microsatellite stable with a median mutation burden of 3.0 single nucleotide alterations (SNAs) per Mb, around 20% higher than S-CRC (2.5 SNAs/Mb), and consistent with elevated aging-associated mutational processes. Non-dysplastic mucosa had considerable mutation burden (median 47 SNAs), including mutations shared with the neighboring CA-CRC, indicating a pre-cancer mutational field. CA-CRCs were often near-triploid (42%) or near-tetraploid (21%) and phylogenetic analysis revealed that copy number alterations (CNAs) began to accrue in non-dysplastic bowel, but the LGD/HGD transition often involved a punctuated catastrophic CNA increase.

  Study EGAS00001003028

• Tissue-specific mutation accumulation in human adult stem cells during life

  Gradual accumulation of mutations in human adult stem cells during life is associated with various age-related diseases, including cancer. The number of stem cell divisions throughout life is believed to be a major determinant for mutation accumulation and could explain the extreme variation of cancer incidence across different organs. Yet, mutation patterns and rates of healthy adult stem cells remain unknown. Here, we determined genome-wide mutation patterns in primary adult stem cells of the small intestine, colon and liver of human donors with ages ranging from 3 to 87 years. We find that the number of mutations increases linearly with age up to several thousand mutations per cell at 87 years of age, while mutation spectra remain constant throughout life. Small intestine and colon stem cells have a 2-fold higher mutation rate per year compared with liver stem cells. These differences could be exclusively attributed to the mutagenic action of spontaneous deamination of cytosine residues and may reflect the high stem cell division rate in these tissues. The genomic distribution of somatic mutations is non-random and predominantly associated with DNA replication dynamics in the small intestine and colon, and with transcription in the liver. These results indicate that a stable balance between various mutagenic and DNA repair processes is maintained throughout life and that the activity of these processes in adult stem cells varies between tissues.

  Study EGAS00001001682

• Germline variants collaborate with somatic mutations and initiate and/or drive disease in primary myelodysplastic syndrome (MDS) and therapy-related myeloid neoplasms (t-MN)

  Therapy-related myeloid neoplasms (T-MN) are poorly characterized secondary hematological malignancies following chemotherapy/radiotherapy exposure. We compared the clinical and mutational characteristics of T-MN (n=129) and primary myelodysplastic syndrome (P-MDS, n=108) patients. Although the somatic mutation frequency was similar between T-MN and P-MDS patients (93% in both groups), the pattern was distinct. TP53 mutations were more frequent in T-MN (29.5% versus 7%), while spliceosomal complex mutations were more common in P-MDS (56.5% versus 25.6%).In contrast to P-MDS, the ring sideroblasts (RS) phenotype was not associated with better survival in T-MN, most probably due to genetic association with TP53 mutations. SF3B1 was mutated in 96% of P-MDS with ≥15% RS, but in only 32% T-MN. TP53 mutations were detected in 92% T-MN with ≥15% RS and SF3B1 wild-type cases. Interestingly, T-MN and P-MDS patients with “Very low” or “Low” Revised International Prognostic Scoring System (IPSS-R) showed similar biological and clinical characteristics. In a Cox-regression analysis, TP53 mutation was a poor prognostic factor in T-MN, independent of IPSS-R cytogenetic, disease-modifying therapy and NRAS mutation.Our data have direct implications for T-MN management and provide evidence that, in addition to conventional disease parameters, mutational analysis should be incorporated in T-MN risk stratification.

  Study EGAS00001003547

• Low-coverage whole genome sequencing for a highly selective cohort of severe COVID-19 patients

  Background Despite advances in identification of genetic markers associated to severe COVID-19, the full genetic characterisation of the disease remains elusive. Imputation of low-coverage whole genome sequencing (lcWGS) has emerged as a competitive method to study such disease-related genetic markers as they enable genotyping of most common genetic variants used for genome wide association studies. This study aims at exploring the potential use of imputation in lcWGS for a highly selected severe COVID-19 patient cohort. Findings We generated an imputed dataset of 79 variant call format (VCF) patient files using the GLIMPSE1 tool, each containing, on average, 9.5 million single nucleotide variants. The validation assessment of imputation accuracy yielded a squared Pearson correlation of approximately 0.97 across sequencing platforms, showing that GLIMPSE1 can be used to confidently impute variants with minor allele frequency up to approximately 2% in Spanish ancestry individuals. We conducted a comprehensive analysis on the patient cohort, examining hospitalisation and intensive care utilisation, sex and age-based differences, and clinical phenotypes using a standardised set of medical terms specifically developed to characterise severe COVID-19 symptoms for this cohort. Conclusion This dataset highlights the utility and accuracy of lcWGS imputation in the study of COVID-19 severity, setting a precedent for other applications in resource-constrained environments. The methods and findings presented here may be leveraged in future genomic projects, providing vital insights for health challenges like COVID-19.

  Study EGAS00001007573

• Preclinical Pediatric Molecular Analysis for Therapy Choice (MATCH)

  Purpose: Pediatric solid tumors arise from diverse tissues during development and exhibit a wide range of molecular, cellular and genetic features. This diversity, combined with the low incidence of pediatric cancer makes it increasingly difficult to personalize therapy for individual patients based on the unique features of their tumors. Therefore, well-credentialed preclinical models that capture the diversity and heterogeneity of pediatric solid tumors are essential for identifying molecular targeted therapeutics for precision medicine. Experimental Design: Here, we report 281 orthotopic patient derived xenografts (O-PDXs) from 224 patients representing 24 different types of pediatric solid tumors. We have performed genomic characterization of the O-PDXs and compared them to their corresponding patient tumors. To demonstrate the feasibility and utility of using such a diverse collection of O-PDXs in preclinical studies, we performed a preclinical pediatric precision medicine trial based on the NCI-COG Pediatric MATCH trial enrollment criteria. We also tested molecular targeted therapy for a novel oncogenic fusion recently reported in pediatric melanoma and precision drug delivery using nano-liposomal irinotecan. Results: In the preclinical match study, we identified one specific single agent response to LY3023414 in an osteosarcoma model with a TSC2 variant of unknown significance. Combinations of nano-liposomal irinotecan were effective in treating both Ewing sarcoma and rhabdomyosarcoma xenografts. Conclusion: Our studies demonstrate the value of large, well-credentialed preclinical models for future precision medicine in pediatric oncology using single agents, drug combinations and novel drug formulations.

  Study EGAS00001008011

• Genetic and epigenetic characterization of adenoid cystic carcinoma

  Translocation events are frequent in cancer and may create chimeric fusions or regulatory rearrangements that drive oncogene overexpression. Although regulatory rearrangements are increasingly recognized in hematopoietic and even solid tumors, the underlying mechanisms remain obscure. Here we identify super-enhancer translocations that drive overexpression of the oncogenic transcription factor MYB as a unifying theme in adenoid cystic carcinoma (ACC). Whole genome sequencing data and chromatin state maps for 13 primary tumors and xenografts reveal distinct chromosomal rearrangements that juxtapose super-enhancer clusters to the MYB locus. Chromosome conformation capture confirms that the translocated enhancers physically interact with the MYB promoter. Remarkably, the MYB product itself binds to the translocated enhancers, thus creating a positive feedback loop that sustains its own expression. MYB also binds a large number of active enhancers that drive opposing regulatory programs in alternate cell lineages in ACC. MYB cooperates with the transcription factor TP63 in the myoepithelial component of the tumors, but promotes a Notch program in the luminal epithelial component. Bromodomain inhibitors, which disrupt enhancer function, slow tumor growth in ACC primagraft models in vivo, but are ineffective against high grade tumors that harbor activating mutations in the Notch pathway. Thus, our study identifies super-enhancer translocations as a unifying feature of ACC, and provides insight into the mechanism by which sustained MYB overexpression drives alternate cell fates in this disease.Processed ChIP-seq data is available on GEO under accession number GSE76465.

  Study EGAS00001001457

• A no-stop mutation in MAGEB4 is a possible cause of rare X-linked azoospermia and oligozoospermia in a consanguineous Turkish family

  Purpose: The purpose of this study was to identify mutations that cause non-syndromic male infertility using whole exome sequencing of family cases.Methods: We recruited a consanguineous Turkish family comprising nine siblings with male triplets; two of the triplets were infertile as well as one younger infertile brother. Whole exome sequencing (WES) performed on two azoospermic brothers, identified a mutation in the MAGEB4 gene which was confirmed via Sanger sequencing and then screened for on control groups and unrelated infertile subjects. The effect of the mutation on mRNA and protein levels were tested after in vitro cell transfection. Structural features of MAGEB4 were predicted throughout the conserved MAGE domain.Results: The novel single base substitution (c.1041A>T) in the X-linked MAGEB4 gene was identified as a no-stop mutation. The mutation is predicted to add 24 amino acids to the C-terminus of MAGEB4. Our functional studies were unable to detect any effect either on mRNA stability, intracellular localization of the protein or the ability to homo / heterodimerize with other MAGE proteins. We thus hypothesize that these additional amino acids may affect the proper protein interactions with MAGEB4 partners. Conclusion: The whole exome analysis of a consanguineous Turkish family revealed MAGEB4 as a possible new X-linked cause of inherited male infertility. This study provides the first clue to the physiological function of a MAGE protein.

  Study EGAS00001002272

• Successful BRAF/MEK-inhibition in a young patient with BRAF V600E-mutated extrapancreatic acinar cell carcinoma (HIPO-021)

  Pancreatic acinar cell carcinoma (PAC) is a rare disease with a poor prognosis. Treatment options for metastatic PAC are limited and often follow chemotherapeutic regimens for pancreatic ductal adenocarcinoma. Although recurrent genomic alterations, such as BRAF fusions and defects in genes involved in homologous recombination DNA repair, have been described in PAC, data on the clinical efficacy of molecularly guided, targeted treatment are scarce. Here we describe the case of a 27-yr-old patient with BRAF V600E-mutated PAC who was successfully treated with a combination of BRAF and MEK inhibitors. The patient presented to our clinic with abdominal pain and weight loss. Imaging showed extensive retroperitoneal disease as well as mediastinal lymphadenopathy. Because of elevated alpha-fetoprotein (AFP) levels and inconclusive histologic findings, a germ cell tumor was suspected; however, PEI chemotherapy was unsuccessful. A repeat biopsy yielded the diagnosis of PAC and treatment with FOLFIRINOX was initiated. Comprehensive molecular profiling within the MASTER (Molecularly Aided Stratification for Tumor Eradication Research) precision oncology program revealed a somatic BRAF V600E mutation and a germline PALB2 stop-gain mutation. Therapy was therefore switched to BRAF/MEK inhibition, resulting in almost complete remission and disease control for 12 mo and a remarkable improvement in the patient's general condition. These results indicate that BRAF alterations are a valid therapeutic target in PAC that should be routinely assessed in this patient population.

  Study EGAS00001004282

• UK10K COHORT TWINSUK

  The UK10K project proposes 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 directly analyse 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 TwinsUK samples will be part of the cohort study and will undergo whole genome sequencing. For further information with regard to this cohort please contact Brent Richards (brent.richards@kcl.ac.uk) or Nicole Soranzo (ns6@sanger.ac.uk).

  Study EGAS00001000108

• UK10K NEURO ASD BIONED

  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 BioNED (Biomarkers for Childhood onset neuropsychiatric disorders) study has been carrying out detailed phenotypic assessments evaluating children with an autism spectrum disorder. These assessments included ADI-R, ADOS, neuropsychology, EEG etc. There are 56 DNA samples from this study (25 extracted from blood). For further information with regard to this cohort please contact Patrick Bolton (patrick.bolton@kcl.ac.uk).

  Study EGAS00001000111

• UK10K NEURO ASD GALLAGHER

  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. This is an Irish sample set of individuals with ASD (approximately 50% with comorbid intellectual disability). Individuals have been diagnosed with ADI/ ADOS, measures of cognition/ adaptive function. They represent a more severe, narrowly defined cohort of ASD subjects. Family histories are available for some with measures of broader phenotype. For further information on this cohort please contact Nadia Bolshakova (bolshakn@tcd.ie).

  Study EGAS00001000112

• Performance comparison of three DNA extraction kits on human whole-exome data from formalin-fixed paraffin-embedded normal and tumor samples

  We generated 42 human whole-exome sequencing data sets from fresh-frozen (FF) and FFPE samples. These samples include normal and tumor tissues from two different organs (liver and colon), that we extracted with three different FFPE extraction kits (QIAamp DNA FFPE Tissue kit and GeneRead DNA FFPE kit from Qiagen, Maxwell\textsuperscript{TM} RSC DNA FFPE Kit from Promega). Variant calling analysis shows a very high rate of concordance between matched FF / FFPE pairs and equivalent performance for the three kits we analyzed. We find a significant variation in the difference of total number of variants called between FF and FFPE samples for the three different FFPE DNA extraction kits. Coverage analysis shows that FFPE samples have less good indicators than FF samples, yet the coverage quality remains above accepted thresholds. We detect limited but significant variations in coverage indicator values between the three FFPE extraction kits. Globally, the GeneRead and QIAamp kits have better variant calling and coverage indicators than the Maxwell kit on the samples used in this study, although this kit performs better on some indicators and has advantages in terms of practical usage. Taken together, our results confirm the potential of FFPE samples analysis for clinical genomic studies, but also indicate that the choice of a FFPE DNA extraction kit should be done with careful testing and analysis beforehand in order to maximize the accuracy of the results.

  Study EGAS00001002631

• Identification of a Fumarate-Hydratase Deficient-like RCC Subtype with Convergent Pathway Alterations

  The development of molecularly targeted therapies builds upon the identification of tumor clades with shared molecular abnormalities. Hereditary leiomyomatosis and renal cell carcinoma (HLRCC) is a familial syndrome resulting from germline mutation of the two-hit tumor suppressor gene, Fumarate Hydratase (FH). HLRCC renal tumors present often in young adults, exhibit unique histological features and are particularly aggressive. HLRCC is uncommon and poorly understood, and whether there are sporadic renal tumors with similar characteristics, is unclear. Notably, FH can also be inactivated somatically in RCC, and should be included among FH-deficient (FHD) kidney tumors. We identified a tumor clade without FH mutations but with similar morphological and molecular characteristics as well as aggressive clinical behavior, which we name FHDL (FHD-like). We show that FHD/FHDL tumors represent a unique subtype of CIMP (CpG Island Methylator Phenotype) with the highest levels of overall genome methylation and CDKN2A silencing. FHD/FHDL tumors are characterized by convergent activation of the Hippo pathway as determined by frequent and mutually exclusive mutations in NF2 and PTPN14. AKR1B10, an NRF2 target, was upregulated in FHD/FHDL tumors as well as GATA3, which suggests a shared origin in the renal medulla and is consistent with topological studies. There is loss of 5-hydroxymethylcytosine in both FHD and FHDL tumors. These data define a group of particularly aggressive renal cancers with shared molecular features and potential targets for therapeutic intervention.

  Study EGAS00001002646

• Whole genome study of Hurthle cell thyroid carcinoma

  Oncocytic thyroid carcinoma, also known as Hürthle cell thyroid carcinoma, accounts for only a small percentage of all thyroid cancers. However, this malignancy often presents at an advanced stage and poses unique challenges to patients and clinicians. Surgical resection of the tumor accompanied in some cases by radioactive iodine treatment, radiation and chemotherapy are the established modes of therapy. Knowledge of the perturbed oncogenic pathways can provide better understanding of the mechanism of disease and thus opportunities for more effective clinical management. Initially, two oncocytic thyroid carcinomas and their matched normal tissues were profiled using whole genome sequencing. Subsequently, 72 oncocytic thyroid carcinomas, one cell line and 5 Hürthle cell adenomas were examined by targeted sequencing for the presence of mutations in multiple endocrine neoplasia I (MEN1) gene. We report the identification of MEN1 loss-of-function mutations in approximately 7% of patients diagnosed with oncocytic thyroid carcinoma. Whole genome sequence data also revealed large regions of copy number variation encompassing nearly the entire genomes of these tumours. Menin, a ubiquitously expressed nuclear protein, is a well-characterized tumor suppressor whose loss is the cause of multiple endocrine neoplasia type 1 syndrome. Menin is involved in several major cellular pathways such as regulation of transcription, control of cell cycle, apoptosis and DNA damage repair pathways. Mutations of this gene in a subset of Hürthle cell tumors point to a potential role for this protein and its associated pathways in thyroid tumorigenesis.

  Study EGAS00001000940

• H3.3G34R/V-transformed interneuron progenitors co-opt PDGFRA for gliomagenesis

  Histone H3.3 glycine 34 to arginine/valine (H3.3G34R/V) mutations occur in deadly hemispheric high-grade gliomas. These tumors show exquisite regional and temporal specificity, suggesting a developmental context permissive to the effects of G34R/V mutations. Here we present the molecular landscape of G34R/V gliomas (n=85) and show that 50% bear activating mutations in PDGFRA, with strong selection pressure for PDGFRAMUT clones at recurrence. We show that G34R/V tumors arise in interneuron progenitors of the foetal ventral forebrain expressing GSX2 and the DLX family of homeobox transcription factors, where terminal neuronal differentiation is impaired through aberrant G34R/V-mediated H3K27me3. Frequent co-occurrence of G34R/V & PDGFRAMUT is facilitated in this interneuron lineage-of-origin as PDGFRA forms an aberrant chromatin loop with the adjacent GSX2, hijacking its active chromatin conformation. At the single-cell level, G34R/V tumours entirely lack oligodendroglial transcriptional programs prominent in other glioma entities, and instead harbour dual neuronal and astroglial compartments. CRISPR-removal of H3.3G34R/V does not impact tumorigenicity suggesting this mutation becomes dispensable, while PDGFRAMUT are potently oncogenic regardless of G34 mutation. Collectively, our results suggest that G34R/V gliomas arise in foetal interneuron progenitors unable to undergo terminal differentiation, enabling co-option of PDGFRA through inappropriate expression and activating mutations to promote gliogenesis and oncogenicity. Reliance on PDGFRA for oncogenesis may be of therapeutic opportunity in G34R/V glioma.

  Study EGAS00001004301

• The spatio-temporal evolution of lymph node spread in early breast cancer

  Purpose: The most significant prognostic factor in early breast cancer is lymph node involvement. This stage between localised and systemic disease is key to understanding breast cancer progression, however our knowledge of the evolution of lymph node malignant invasion remains limited, as most currently available data derive from primary tumours. Experimental design: In 11 treatment-naïve node positive early breast cancer patients without clinical evidence of distant metastasis, we investigated lymph node evolution using spatial multi-region sequencing (n=78 samples) of primary and lymph node deposits and genomic profiling of matched longitudinal circulating tumour DNA (ctDNA). Results: Linear evolution from primary to lymph node was rare (1/11) whereas the majority of cases displayed either early divergence between primary and nodes (4/11), or no detectable divergence (6/11) where both primary and nodal cells belonged to a single recent expansion of a metastatic clone. Divergence of metastatic subclones was driven in part by APOBEC. Longitudinal ctDNA samples from 2 subjects taken peri-operatively reflected the two major evolutionary patterns and demonstrate that private mutations can be detected even from early metastatic nodal deposits. Moreover, node removal resulted in disappearance of private lymph node mutations in ctDNA.Conclusions: This study sheds new light on a crucial evolutionary step in the natural history of breast cancer, demonstrating early establishment of axillary lymph node metastasis in a substantial proportion of patients.

  Study EGAS00001002947

• Genomic Alterations in Gingivo-buccal Cancer: ICGC-India Project_YR01

  As a part of the ICGC, India has undertaken genomic studies on gingivobuccal cancer of the oral cavity, which is the most prevalent form of cancer among men in India. There are various known environmental (life-style) correlates of this cancer, the most important of which are tobacco chewing and HPV infection. Paired DNA samples – isolated from the tumor tissue and from the blood of fifty patients – have been analyzed to catalog germline and somatic mutations. Association between each observed genomic alteration and exposure to environmental risk factors is being explored.Detailed clinical characterization of the patients, collection of data on demographic and environmental exposures, and isolation of DNA samples from blood and tumor tissues collected from each patient are being done at the Advanced Centre for Research, Treatment and Education on Cancer, Mumbai. Exome capture and deep resequencing are being performed at the National Institute of Biomedical Genomics, Kalyani. We have analyzed the exomes of paired blood and tumor DNA samples using two orthogonal platforms – exome capture by Roche-Nimblegen SeqCap Ez kit followed by sequencing on Roche-454 GS FLX, and exome capture by Illumina TrueSeq Exome Enrichment kit followed by sequencing on Illumina HiSeq-2000. Each exome is being sequenced at a mean depth of 30x. In addition, we have used Illumina Omni Quad SNP-chips to generate genotype data at 1.14 million genomic markers. The average concordance between SNP-chip and sequence data is 99.8%.

  Study EGAS00001000249

• GCAT | Genomes for life: cohort study of the genomes of Catalonia

  The GCAT Study have recruited 20 000 participants aged 40–65 years. Participants who agreed to take part in the study completed a self-administered computer-driven questionnaire, and underwent blood pressure, cardiac frequency and anthropometry measurements. For each participant, blood plasma, blood serum and white blood cells are collected at baseline. A total of 5459 genomic profiles have been characterised by comprehensive genotyping. Genome-wide genotypes have been generated using Illumina Infinium SNP-bead array technology. We chose the Multi-Ethnic Global (MEGAEX, V.2) consortium array, a multipurpose, multiethnic genotyping array with two million selected markers (including previously described germline mutations, insertions-deletions (InDels) and SNPs).We have strictly followed the standard manufacturer recommended automated protocol for the Infinium HTS Assay scanned with a HiScan confocal scanner (Illumina, San Diego, California, USA). Genome Studio V.2011.1 has been used for raw data analysis. Genotyping was performed at the Genomics and Bioinformatics Unit of the PMPPC Institute for Health Science Research Germans Trias i Pujol, in Badalona, Spain. Future plans The first follow-up study started in December 2017 and will end by March 2018. Residences of all subjects will be geocoded during the following year. Several genomic analyses are ongoing, and metabolomic and genomic integrations will be performed to identify underlying genetic variants, as well as environmental factors that influence metabolites. http://dx.doi.org/10.1136/bmjopen-2017-018324

  Study EGAS00001003018

• ICU_transcriptomics__Assessing_the_role_of_the_host_immune_response_in_patients_with_ventilator_associated_pneumonia

  The aim of the proposed project is to provide transcriptomic sequencing data analysis of cells isolated from bronchioalveolar lavage and blood samples from critically ill patients with pneumonia on an intensive care unit in order to investigate the host contribution to disease. The lung microbiome of these patients will be analysed in parallel by deep sequencing and 16S sequencing. These data will provide important information to the development of strategies aimed at improving patient outcomes, reducing the development of antimicrobial resistance and facilitating infection control measures in hospital. Specifically, the proposed work will be a proof of principle study of 100 patients. We will use transcriptomic data to define the host contribution to the symptoms of these critically ill patients diagnosed with pneumonia. The data will be correlated with cytokine/chemokine concentrations, routine microbiology analysis and pathogen genomic data from the same samples to determine the ecology of the lungs of these patients. We will look for associations between clinical features, host response and pathogen genomics in order to identify new diagnostic markers. This study should provide important information to improve the speed and accuracy of diagnosis - resulting in more targeted and successful treatment of critically ill patients in the ICU. 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 EGAS00001003074

• Whole Tumor spatial heterogeneity in metastatic melanoma

  Heterogeneous inter-tumoral responses, sustained periods of apparent clinical benefit despite lack of objective response to various therapies, and even spontaneous remission are well known within a subpopulation of advanced melanoma patients. The molecular and cellular dynamics facilitating long-term survival remain poorly defined, particularly in the current era ofexposure to multiple potentially active therapies. We studied an exceptional case of long-term survival in a patient with non-responding metastatic melanoma in order to characterize the clonal and microenvironmental factors.We performed immunogenomic analysis of 3 metachronous tumors, and of 67 intratumor sub-regions of one metastatic lesion, usingwhole exome sequencing, RNA-sequencing, immunohistochemistry, and T cell receptor sequencing. Longitudinal analysesidentified mutations in severalknownresistance genes affecting distinctmetastases. Intra-tumoral immune signatures revealed spatially-distinct pockets of immune activation and suppression throughout the metastasis. A specific T cell VBeta CDR3 rearrangement was identified as both dominant and recurrent not only across multiple spatial points within a single tumor mass, but also across metachronous tumors spanning the patient’s disease course. Immunophenotyping of this T cell population with single-cell sequencingsuggestedrepeatedT-cell priming events leading to the persistence of both activated and exhausted T cells bearing this same TCR-Beta at any given time.Our findings highlight genomic and immune heterogeneityin metastatic disease, extending to the intra-tumoral leveland provide a rationale for integrative spatiotemporal analysis of metastatic lesions in the context of immune checkpoint blockade.

  Study EGAS00001003292

• Exploiting evolutionary steering in cancer therapy

  Drug resistance, mediated by intra-tumour heterogeneity and clonal evolution, is arguably the biggest problem in cancer therapy today. However, evolving resistance to one drug may come at the cost of increased sensitivity to another due to so-called evolutionary trade-offs. This weakness can be exploited in the clinic using an approach called ‘evolutionary herding’ that aims at controlling the tumour cell population to delay or prevent resistance. However, model systems able to recapitulate cancer evolution experimentally are lacking and current in vitro techniques based on small populations, re-plating and escalating dose, are unsuitable to develop evolutionary herding strategies. We present a novel methodology for evolutionary herding in vitro and ex vivo using patient-derived organoids. Our approach is based on a combination of single-cell barcoding, very large populations of 108-109 cells grown without re-plating, realistic high drug doses, time-course monitoring of cancer clones, and mathematical modelling of tumour evolutionary dynamics. We demonstrate evolutionary herding in non-small cell lung cancer in vitro and in patient-derived colorectal cancer organoids (PDO). We show that herding causes controlled evolutionary bottlenecks that lead to collateral sensitivity. Through genomic analysis, we were also able to determine the mechanisms that drive such sensitivity. Our approach allows modelling evolutionary trade-offs experimentally to test patient-specific evolutionary herding strategies that can be translated into the clinic to control treatment resistance.

  Study EGAS00001003200

• Probabilistic cell type assignment of single-cell transcriptomic data reveals spatiotemporal microenvironment dynamics in human cancers

  Single-cell RNA sequencing (scRNA-seq) has transformed biomedical research, enabling decomposition of complex tissues into disaggregated, functionally distinct cell types. For many applications, investigators wish to identify cell types with known marker genes. Typically, such cell type assignments are performed through unsupervised clustering followed by manual annotation based on these marker genes, or via "mapping" procedures to existing data. However, the manual interpretation required in the former case scales poorly to large datasets, which are also often prone to batch effects, while existing data for purified cell types must be available for the latter. Furthermore, unsupervised clustering can be error-prone, leading to under- and over- clustering of the cell types of interest. To overcome these issues we present CellAssign, a probabilistic model that leverages prior knowledge of cell type marker genes to annotate scRNA-seq data into pre-defined and de novo cell types. CellAssign automates the process of assigning cells in a highly scalable manner across large datasets while simultaneously controlling for batch and patient effects. We demonstrate the analytical advantages of CellAssign through extensive simulations and exemplify real-world utility to profile the spatial dynamics of high-grade serous ovarian cancer and the temporal dynamics of follicular lymphoma. Our analysis reveals subclonal malignant phenotypes and points towards an evolutionary interplay between immune and cancer cell populations with cancer cells escaping immune recognition.

  Study EGAS00001003452

• Genome-wide study of the effect of blood collection tubes on the cell-free DNA methylome

  Background The methylation pattern of cfDNA, isolated from liquid biopsies, is gaining substantial interest for diagnosis and monitoring of diseases. We have evaluated the impact of type of blood collection tube and time delay between blood draw and plasma preparation on bisulfite-based cfDNA methylation profiling. Methods 15 tubes of blood were drawn from three healthy volunteer subjects (BD Vacutainer K2E EDTA spray tubes, Streck Cell-Free DNA BCT tubes, PAXgene Blood ccfDNA tubes, Roche Cell-Free DNA Collection tubes and Biomatrica LBgard blood tubes in triplicate). Samples were either immediately processed or stored at room temperature for 24 or 72 hours before plasma preparation. DNA fragment size was evaluated by capillary electrophoresis. Reduced representation bisulfite sequencing was performed on the cell-free DNA isolated from these plasma samples. We evaluated the impact of blood tube and time delay on several quality control metrics. Results All preservation tubes performed similar on the quality metrics that were evaluated. Furthermore, a considerable increase in cfDNA concentration and the fraction of it derived from NK cells was observed after a 72-hour time delay in EDTA tubes. Conclusion The methylation pattern of cfDNA is robust and reproducible in between the different preservation tubes. EDTA tubes processed as soon as possible, preferably within 24 hours, are the most cost effective. If immediate processing is not possible, preservation tubes are valid alternatives.

  Study EGAS00001004271

• Long-read trio sequencing of unsolved patients with intellectual disability

  Long-read sequencing (LRS) has the potential to comprehensively identify all medically-relevant genome variation, including variation commonly missed by short read sequencing (SRS) approaches. To determine this potential, we performed LRS up to 40x coverage using the Pacific Biosciences Sequel instrument for 5 patient-parent trios with intellectual disability, whose aetiology remained unresolved after SRS exomes and genomes. A first systematic assessment of LRS coverage showed that ~35Mb of the genome was only accessible by LRS and not SRS. Structural variant (SV) calling yielded on average 28,292 SV calls per individual, totalling 12.9 MB of sequence. Trio-based analyses showed concordance for up to 95% of these SV calls across the genome, and 80% of the LRS SV calls were not identified by SRS. De novo mutation analysis did not identify any de novo SVs, confirming these are rare events. Because of high sequence coverage, we were also able to call single nucleotide substitutions. On average, we identified 3 million substitutions per genome, with a Mendelian inheritance concordance of up to 97%. Of these, approximately 100,000 were located in the ~35Mb of the genome that was only captured by LRS. Moreover, these variants affected the coding sequence of 64 genes, including 32 known to cause Mendelian disorders. Although no clinical diagnoses were made in these 5 patients, our data show the potential added value of LRS compared to SRS for identifying medically relevant genome variation.

  Study EGAS00001004319

• The Medical Genome Reference Bank: a whole genome data resource of 4,000 healthy elderly individuals.

  Allele frequency data from human reference populations is of increasing value for filtering and assignment of pathogenicity to genetic variants. Aged and healthy populations are more likely to be selectively depleted of pathogenic alleles, and therefore particularly suitable as a reference populations for the major diseases of clinical and public health importance. However, reference studies of the healthy elderly have remained under represented in human genetics. We have developed the Medical Genome Reference Bank (MGRB), a large scale comprehensive whole genome dataset of confirmed healthy elderly individuals, to provide a publicly accessible resource for health related research, and for clinical genetics. It also represents a useful resource for studying the genetics of healthy aging. The MGRB comprises 4,000 healthy, older individuals with no reported history of cancer, cardiovascular disease or dementia, recruited from two Australian community based cohorts. DNA derived from blood samples will be subject to whole genome sequencing. The MGRB will measure genome wide genetic variation in 4,000 individuals, mostly of European decent, aged 60 to 95 years (mean age 75 years). The MGRB has committed to a policy of data sharing, employing a hierarchical data management system to maintain participant privacy and confidentiality, whilst maximizing research and clinical usage of the database. The MGRB will represent a dataset of international significance, broadly accessible to the clinical and genetic research community.

  Study EGAS00001003511

• Efficacy of JAK/STAT pathway inhibition in murine xenograft models of early T-cell precursor (ETP) acute lymphoblastic leukemia

  Early T-cell precursor (ETP) acute lymphoblastic leukemia (ALL) is a recently described subtype of T-ALL characterized by a unique immunophenotype and genomic profile as well as a high rate of induction failure. Frequent mutations in cytokine receptor and JAK/STAT signaling pathways led us to hypothesize that ETP-ALL is dependent on JAK/STAT signaling. Here we demonstrate aberrant activation of the JAK/STAT pathway in ETP-ALL blasts relative to non-ETP T-ALL. Moreover, ETP-ALL showed hyperactivation of STAT5 in response to IL7, an effect that was abrogated by the JAK1/2 inhibitor ruxolitinib. In vivo, ruxolitinib displayed activity in 6/6 patient-derived murine xenograft models of ETP-ALL, with profound single-agent efficacy in 5 models. Ruxolitinib treatment decreased peripheral blast counts relative to pre-treatment levels and compared to control (P<0.01) in 5/6 ETP-ALL xenografts, with marked reduction in mean splenic blast counts (P<0.01) in 6/6 samples. Surprisingly, both JAK/STAT pathway activation and ruxolitinib efficacy were independent of the presence of JAK/STAT pathway mutations, raising the possibility that the therapeutic potential of ruxolitinib in ETP-ALL extends beyond those cases with JAK mutations. These findings establish the preclinical in vivo efficacy of ruxolitinib in ETP-ALL, a biologically distinct subtype for which novel therapies are needed.

  Study EGAS00001001146

• Initial whole genome sequencing of plasma cell neoplasms in First Responders exposed to the World Trade Center attack of September 11, 2001

  The World Trade Center (WTC) attack of September 11, 2001 created an unprecedented environmental exposure to known and suspected carcinogens. High incidence of multiple myeloma (MM) and precursor conditions has been reported among first responders to the WTC disaster. To expand on our prior screening studies, and to characterize the genomic impact of the exposure to known and potential carcinogens in the WTC debris, we were motivated to perform whole genome sequencing (WGS) of WTC first responders and recovery workers who developed a plasma cell disorder after the attack. PATIENTS AND METHODS: We performed WGS of 9 CD138-positive bone marrow mononuclear samples from patients who were diagnosed with plasma cell disorders after the WTC disaster. RESULTS: No significant differences were observed in comparing the post-WTC driver and mutational signatures landscape with 110 previously published WGS from 56 patients with MM and the CoMMpass WGS cohort (n=752). Leveraging constant activity of the single base substitution mutational signatures 1 and 5 over time, we estimated that tumor-initiating chromosomal gains were windowed to both pre- and post-WTC exposure. CONCLUSIONS: Although limitations in sample size preclude any definitive conclusions, our findings suggest that the observed increased incidence of plasma cell neoplasms in this population is due to complex and heterogeneous effects of the WTC exposure that may have initiated or contributed to progression of malignancy.

  Study EGAS00001004467

• The clonal and mutational evolution spectrum of primary triple negative breast cancers

  Primary triple negative breast cancers (TNBC) represent approximately 16% of all breast cancers and are a tumour type defined by exclusion,for which comprehensive landscapes of somatic mutation have not been determined. Here we show in 104 early TNBC cases, that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some exhibiting only a handful of somatic aberrations in a few pathways, whereas others contain hundreds of somatic events and multiple pathways implicated. Integration with matched whole transcriptome sequence data revealed that only 36% of mutations are expressed. By examining single nucleotide variant (SNV) allelic abundance derived from deep re-sequencing (median > 20,000 fold) measurements in 2414 somatic mutations, we determine for the first time in an epithelial tumour, the relative abundance of clonal genotypes among cases in the population. We show that TNBC vary widely and continuously in their clonal frequencies at the time of diagnosis, with basal subtype TNBC exhibiting more variation than non-basal TNBC. Although p53 and somatic mutations appear clonally dominant compared with other pathways, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal and cell shape/motility proteins occurred at lower clonal frequencies, suggesting they occurred later during tumour progression. Taken together our results show that future attempts to dissect the biology and therapeutic responses of TNBC will require the determination of individual tumour clonal genotypes.

  Study EGAS00001000132

• 33 patients with Monoclonal Gammopathy of Undetermined Significance (MGUS) had whole exome sequencing performed. Cells were sorted by FACSAria using CD138, CD19 and CD56 to obtain a pure abnormal plasma cell population. Generated somatic variants were compared to a previous study of 463 multiple myeloma patients.

  Monoclonal gammopathy of undetermined significance (MGUS) is a premalignant precursor of multiple myeloma (MM) with a 1% risk of progression per year. Although targeted analyses have shown the presence of specific genetic abnormalities such as IGH translocations, RB1 deletion, 1q gain, hyperdiploidy or RAS genes mutations, little is known about molecular mechanism of malignant transformation. We have performed whole exome sequencing together with SNP array analysis in 33 flow-cytometry separated abnormal PC samples of MGUS patients to describe somatic gene mutations and chromosome changes at the genome-wide level. Non-synonymous mutations (NS-SNVs) and copy number alterations (CNAs) were present in 97.0% and in 63.6% of cases, respectively. Importantly, the number of somatic mutations was significantly lower in MGUS compared to MM (p<10-4) and we have identified 6 myeloma significantly mutated genes which are KRAS, NRAS, DIS3, HIST1H1E, EGR1 and LTB in the MGUS dataset. We also found a positive correlation with increasing chromosome changes and somatic mutations. IGH translocations were present in 27.3% of cases comprising t(4;14), t(11;14), t(14;16) or t(14;20) and were in a similar frequency to MM, which corresponded with primary lesion hypothesis. Data from this study showed MGUS is a genetically comprehensive disease, however overall genetic instability is significantly lower compared to MM.

  Study EGAS00001001658

• RRBS sequencing of 7 tumour regions and a normal sample from a single TRACERx patient.

  Background: Epigenetic heterogeneity within a tumour can play an important role in tumour evolution and the emergence of resistance to treatment. It is increasingly recognised that the study of DNA methylation (DNAm) patterns along the genome - so-called `epialleles' - offers greater insight into epigenetic dynamics than conventional analyses which examine DNAm marks individually. Results: We have developed a Bayesian model to infer which epialleles are present in multiple regions of the same tumour. We apply our method to reduced representation bisulfite sequencing (RRBS) data from multiple regions of one lung cancer tumour and a matched normal sample. The model borrows information from all tumour regions to leverage greater statistical power. The total number of epialleles, the epiallele DNAm patterns, and a noise hyperparameter are all automatically inferred from the data. Uncertainty as to which epiallele an observed sequencing read originated from is explicitly incorporated by marginalising over the appropriate posterior densities. The degree to which tumour samples are contaminated with normal tissue can be estimated and corrected for. By tracing the distribution of epialleles throughout the tumour we can infer the phylogenetic history of the tumour, identify epialleles that differ between normal and cancer tissue, and define a measure of global epigenetic disorder. Conclusions: Detection and comparison of epialleles within multiple tumour regions enables phylogenetic analyses, identification of differentially expressed epialleles, and provides a measure of epigenetic heterogeneity.

  Study EGAS00001002484

• Renal_habitat_WXS

  This project aims at correlating genetic diversity (as well as transcriptomic, metabolic and immune diversity) observed in multiple targeted tumour biopsies and circulating tumour DNA with heterogeneity on multiparametric MRI. To this end, 20 patients with RCC (15 patients with localised/locally advanced disease and 5 with mRCC) will undergo multiparametric MRI including metabolic assessment by 13C-MRSI as well as measurements of tumour perfusion, oxygenation and water diffusion. Based on unique combinations of these imaging parameters, tumour habitats will be defined through cluster mapping. Intra-operative biopsy and immediate freezing of tumour cores from different imaging habitats will be performed. Multiregional ex-vivo tumour biopsies will be spatially registered to the MR data by using 3D printed tumour moulds generated from high-resolution volumetric sequences. These moulds allow slicing tumour specimens in identical planes to the MR acquisition and increase the accuracy of tissue sampling. The biopsies will then undergo molecular profiling such as metabolomic and RNA sequencing (including mRNA immune profiling), as well as DNA sequencing to assess driver and clonal mutations in different tumour habitats and guide design of appropriate cfDNA probes. The targeted cfDNA capture panel will be employed to determine the ability of cfDNA to capture the full genetic diversity observed in the primary tumour and to relate cfDNA to various imaging habitats in both renal and pancreatic cancer. Data integration will be undertaken between imaging habitats and molecular profiles for each region.

  Study EGAS00001003703

• Combining transcription factor binding affinities with open chromatin data for accurate gene expression prediction

  The binding and contribution of transcription factors (TF) to cell specific gene expression is often deduced from open-chromatin measurements to avoid cost and labour intensive TF ChIP-seq assays.It is important to develop reliable and fast computational methods for accurate TF binding prediction in open-chromatin regions (OCRs). Here, we report a novel segmentation-based method, TEPIC, to predict TF binding by combining sets of OCRs with position weight matrices.TEPIC can be applied to various open-chromatin data, e.g. DNaseI-seq and NOMe-seq, using either peaks or footprints as input.In addition to open-chromatin data, also Histone-Marks (HMs) can be used in TEPIC to identify candidate TF binding sites.TEPIC computes TF affinities and uses open-chromatin/HM signal intensity as quantitative measures of TF binding strength.Using machine learning techniques, we show that incorporating low affinity binding sites improves our ability to explain gene expression variability compared to the standard presence/absence classification of binding sites.Further, we show that both footprints and peaks capture essential TF binding events and lead to a good prediction performance.In our application, gene-based scores computed by TEPIC with one open-chromatin assay nearly reach the quality of several TF ChIP-seq datasets.Finally, we show that these scores correctly predict known transcriptional regulators as illustrated by the application to novel DNaseI-seq and NOMe-seq data for primary human hepatocytes and CD4+ T-cells, respectively.

  Study EGAS00001002073

• ATRX mutant neuroblastoma is sensitive to EZH2 inhibition via modulation of neuronal differentiation.

  Mutations and structural alterations of the SWI/SNF-like chromatin remodeler ATRX have been reported at high frequency in a number of adult and pediatric tumors1. However, the consequences of ATRX (Alpha Thalassemia/Mental Retardation, X-linked) mutations in cancer and their underlying epigenetic sensitivities remain ill defined. Particularly intriguing are the large N-terminal deletions of ATRX in neuroblastoma that generate in-frame fusion (IFF) proteins1–3 devoid of key chromatin interaction domains. Here we demonstrate that neuroblastoma cells harbouring ATRX IFFs have distinct gene expression programs compared to neuroblastoma cells that are wild type for ATRX. This is due in part to H3K27me3-mediated silencing of REST (RE1 Silencing Transcription Factor) target genes involved in neuronal differentiation. In turn, we find that ATRX IFF cells display exquisite sensitivity to EZH2 inhibition in both adherent and tumorsphere conditions, due in part to derepression of neurogenesis genes, including REST targets. Examination of the epigenomic landscape of a pediatric neuroblastoma tumor harboring an ATRX IFF revealed that H3K27me3 occupies a subset of REST target genes that are transcriptionally silenced, and are sensitive to EZH2 inhibition in our cell-based assays. Thus, our study greatly advances our understanding of indolent neuroblastoma and identifies EZH2 inhibition as a potential therapy for patients with ATRX mutant disease. Further, these studies may be applicable to other ATRX mutant pediatric malignancies, particularly those carrying similar structural alterations of ATRX.

  Study EGAS00001002507

• Mutational landscape and patterns of clonal evolution in relapsed pediatric acute lymphoblastic leukemia

  Relapse of acute lymphoblastic leukemia (ALL) remains a leading cause of childhood death. Prior studies have shown clonal mutations at relapse often arise from relapse-fated subclones that exist at diagnosis. However, the genomic landscape, evolutionary trajectories and mutational mechanisms driving relapse are incompletely understood. In an analysis of 92 cases of relapsed childhood ALL, incorporating multimodal DNA and RNA sequencing, deep digital mutational tracking and xenografting to formally define clonal structure, we identify 50 significant targets of mutation with distinct patterns of mutational acquisition or enrichment. CREBBP, NOTCH1, and Ras signaling mutations rose from diagnosis subclones, whereas variants in NCOR2, USH2A and NT5C2 were exclusively observed at relapse. Evolutionary modeling and xenografting demonstrated that relapse-fated clones were minor (50%), major (27%) or multiclonal (18%) at diagnosis. Putative second leukemias, including those with lineage shift, were shown to most commonly represent relapse from an ancestral clone rather than a truly independent second primary leukemia. A subset of leukemias prone to repeated relapse exhibited hypermutation driven by at least three distinct mutational processes, resulting in heightened neoepitope burden and potential vulnerability to immunotherapy. Finally, relapse-driving sequence mutations were detected prior to relapse using deep digital PCR at levels comparable to orthogonal approaches to monitor levels of measurable residual disease. These results provide a genomic framework to anticipate and circumvent relapse by earlier detection and targeting of relapse-fated clones.

  Study EGAS00001003975

• The Southern African Human Genome Programme

  The Southern African Human Genome Programme (SAHGP) is a national and regional initiative that aims to unlock the unique genetic character of southern African populations. Its vision is to improve quality of life through understanding human genetic diversity. It is timely to advance the SAHGP initiative since global genomic initiatives are increasingly pointing to populations in Africa, more specifically the San and Khoe from southern Africa, as harbouring the most ancient genetic signatures among the living populations of the world. Studies on African populations have already generated a wealth of information upon which the SAHGP can build. This information will be used to promote and support genomic research programmes both nationally and regionally to address critical questions that would benefit the people of our region. The SAHGP aims to make a significant contribution to the understanding of DNA variation among southern Africans and how this impacts on the health of the people of this region. Aims of the SAHGP: To develop capacity for genomic research in southern Africa; To establish sustainable resources for genomic research in the region; and To translate knowledge and information into improvements in human health in the region. The objective is to derive the maximum benefit from the SAHGP, and that the data should be widely shared with the scientific community for the benefit of the people of the sub-continent and beyond.

  Study EGAS00001002639

• Defining the metastasome in colorectal cancer: Implications for hypotheses on metastasis evolution and personalized therapy (HIPO-032)

  Personalized cancer therapy aims at individual genetic changes characterizing primary cancers. Still, however, metastasis is the major clinically unmet problem, also due to an incomplete understanding of its molecular evolution.This study is the first to define the metastasis-specific whole genome landscape of human colorectal primary vs. matched metastatic lesions.Protein coding, ncRNA genes and 3’UTRs harbored about a third each of single nucleotide variations (SNVs)/indels, 19% of them being metastasis-specific. We found novel mutational hills among ncRNAs and 3’UTRs, copy number changes able to explain changes in microRNA expression in metastases, and novel metastasis-specific mutations in the 3’UTRs on important cancer signalling genes. Some metastatic lesions showed targetable mutations not detected in the primaries. Analysis of mutual exclusivity patterns supported a model that has progressive alterations in important colorectal cancer genes and revealed new partners. Metastatic lesions were enriched in SNVs, specific structural variations, and alterations in genes affecting cell adhesion and extracellular matrix interaction.Furthermore, a hepatic stellate activation cascade was enriched in metastases, suggesting genetic programs for site-specific colonization. Allele frequency distribution and mutational signature analysis suggests a common ancestor clone to both the primary tumor and the metastasis, with additional mutagenesis ongoing in both sites independently. Taken together, our results significantly add to hypothesis generation on metastasis evolution, and suggest novel metastasis-specific genomic changes which would be missed by current personalized therapy concepts. (HIPO-032)

  Study EGAS00001002717

• DNA hypermethylation and differential gene expression associated with Klinefelter syndrome

  The molecular basis for the phenotypic traits and morbidity in Klinefelter syndrome (KS) are not clarified. As DNA methylation affect gene expression and thereby play a role in disease susceptibility, we performed genome-wide DNA methylation profiling of leucocytes from peripheral blood samples from 67 KS patients, 67 male controls and 33 female controls, in addition to genome-wide RNA-sequencing profiling in a subset of 9 KS patients, 9 control males and 13 female controls. Characterization of the methylome as well as the transcriptome of both coding and non-coding genes identified a unique epigenetic and genetic landscape of both autosomal chromosomes as well as the X chromosome in KS. A subset of genes show significant correlation between methylation values and expression values. Gene set enrichment analysis of differentially methylated positions yielded terms associated with well-known comorbidities seen in KS. In addition, differentially expressed genes revealed enrichment for genes involved in the immune system, wnt-signaling pathway and neuron development. Based on our data we point towards many new candidate genes (AKAP17A, AMOT, APOB, DACT1, DDX58, DOCK7, EIF2S3, FIGNL1, G3BP1, HENMT1, HERC5, IFI44, IFI44L, IFIT1, IFIT3, ISG15, KANK1, LGALS1, NSD1, PEX10, PLSCR1, RSAD2, SHROOM2, SLC25A6, SPEG, SPON2, TXLNG), which may be implicated in the phenotype and further point towards non-coding genes (RP13-216E22.4, RP13-36G14.4, G087825, G088512), which may be involved in X chromosome inactivation in KS and in the regulation of escape genes.

  Study EGAS00001002797

• Atypical 3q26/MECOM rearrangements genocopy inv(3)/t(3;3) in acute myeloid leukemia

  Acute myeloid leukemia (AML) with inv(3)/t(3;3)(q21q26) is a distinct WHO recognized entity, characterized by its aggressive course and poor prognosis. In this subtype of AML, the translocation of a GATA2 enhancer (3q21) to MECOM (3q26) results in overexpression of the MECOM isoform EVI1 and monoallelic expression of GATA2 from the unaffected allele. The full-length MECOM transcript, MDS1-EVI1, is not expressed as the result of the 3q26 rearrangement. Besides the classical inv(3)/t(3;3), a number of other 3q26/MECOM rearrangements with poor treatment response have been reported in AML. Here we demonstrate, in a group of 33 AML patients with atypical 3q26 rearrangements, MECOM involvement with EVI1 overexpression, but no or low MDS1-EVI1 levels. Moreover, the 3q26 translocations in these AML patients often involve super-enhancers of genes active in myeloid development (e.g. CD164, PROM1, CDK6 or MYC). In more than 50% of these cases allele specific GATA2 expression was observed, either by copy number loss or by an unexplained allelic imbalance. Altogether, atypical 3q26 recapitulate the main leukemic mechanism of inv(3)/t(3;3) AML, namely EVI1 overexpression driven by enhancer hijacking, absent MDS1-EVI1 expression and potential GATA2 involvement. Therefore, we conclude that both atypical 3q26/MECOM and inv(3)/t(3;3) can be classified as a single entity of 3q26-rearranged AMLs. Routine analyses determining MECOM rearrangements, EVI1 and MDS1-EVI1 expression are required to recognize 3q-rearranged AML cases.

  Study EGAS00001004325

• Partner-independent fusion gene detection by multiplexed CRISPR/Cas9 enrichment and long-read Nanopore sequencing

  Fusion genes are hallmarks of various cancer types and important determinants for diagnosis, prognosis and treatment possibilities. The promiscuity of fusion genes with respect to partner choice and exact breakpoint-positions restricts their detection in the diagnostic setting, even for known and recurrent fusion gene configurations. To accurately identify these gene fusions in an unbiased manner, we developed FUDGE: a FUsion gene Detection assay from Gene Enrichment. FUDGE couples target-selected and strand-specific CRISPR/Cas9 activity for enrichment and detection of fusion gene drivers (e.g. BRAF, EWSR1, KMT2A/MLL) - without prior knowledge of fusion partner or breakpoint-location - to long-read Nanopore sequencing. FUDGE encompasses a dedicated bioinformatics approach (NanoFG) to detect fusion genes from Nanopore sequencing data. Our strategy is flexible with respect to target choice and enables multiplexed enrichment for simultaneous analysis of several genes in multiple samples in a single sequencing run. We observe on average a 508-fold on-target enrichment and identify fusion breakpoints at nucleotide resolution - all within two days. We demonstrate that FUDGE effectively identifies fusion genes in cancer cell lines, tumor samples and on whole genome amplified DNA irrespective of partner gene or breakpoint-position in 100% of cases. Furthermore, we show that FUDGE is superior to routine diagnostic methods for fusion gene detection. In summary, we have developed a rapid and versatile fusion gene detection assay, providing an unparalleled opportunity for pan-cancer detection of fusion genes in routine diagnostics.

  Study EGAS00001003964

• Memory-like HCV-specific CD8+ T cells retain a chronic scar after cure of chronic HCV infection

  Chronic hepatitis C virus (HCV) infection is associated with CD8+ T-cell exhaustion characterized by limited effector functions and thus compromised anti-viral activity. Exhausted HCV-specific CD8+ T cells are comprised of memory-like and terminally exhausted CD8+ T-cell subsets. So far, little is not known about the molecular profile and fate of these cells after elimination of chronic antigen stimulation by direct acting antiviral therapy (DAA). Here, we report an antigen-driven molecular core signature underlying exhausted CD8+ T-cell subset heterogeneity in chronic viral infection with a progenitor/progeny relationship of memory-like and terminally exhausted HCV-specific CD8+ T cells via an intermediate stage. Furthermore, transcriptional profiling reveals that the memory-like cells remain after DAA-mediated cure while terminally exhausted HCV-specific CD8+ T-cell subsets are lost. Thus, the memory polarization of the overall HCV-specific CD8+ T-cell response after cure does not result from re-differentiation of exhausted T cells. Consequently, antigen elimination has little impact on the exhausted core signature of memory-like CD8+ T cells that remains clearly different from bona fide T-cell memory. These results identify a molecular signature of T-cell exhaustion that is imprinted like a chronic scar in HCV-specific CD8+ T cells even after HCV cure, highlighting the requirement of re-programming to elicit full effector potential of exhausted T cells.

  Study EGAS00001004538

• Identification of the mutational consequences of precancerous liver disease (including alcohol abuse) on the genomes of human adult stem cells.

  Stem cells are believed to be the cells-of-origin cancer. Nevertheless, their genome remains remarkably stable during life. To increase insight into how these cells can drive tumorigenesis, we set out to determine the genomic integrity of human adult stem cells from patients with liver disease associated with increased cancer risk. To this end, we characterized genome-wide mutation profiles of single ASCs from livers of patients with alcoholic liver disease, non-alcoholic liver disease and primary sclerosing cholangitis. The patients included in our study received a liver transplant due to their liver disease. Strikingly, we found that ASCs in these livers maintain a stable genome and have comparable mutation rates as healthy human liver ASCs. We did observe more nonsense and nonsynonymous mutations in COSMIC cancer genes that are hit by dominant mutations in tumors in comparison to healthy ASCs (n.s.), including two ASCs in alcoholic liver disease patients with a nonsense mutation in PTPRK. Next, we sequenced 5 biopsies of an alcoholic tumour to identify the mutational patterns in the tumour-initiating ASC. We found that the mutational patterns are, again, similar to healthy liver, but the mutational load is slightly higher. This suggests that the increased cancer risk in the liver disease patients is not caused by an increased mutational load, but rather that the precancerous liver environment promotes outgrowth of cells that would normally be selected against.

  Study EGAS00001002983

• Mutational Signature and Transcriptomic Classification Analyses as the Decisive Diagnostic Tools for a Cancer of Unknown Primary with Neuroendocrine Differentiation

  PURPOSE: Cancer of unknown primary (CUP) is a group of metastatic tumors in which the standard diagnostic work-up fails to identify the site of origin of the tumor. The potential impact of precision oncology on this group of patients is large since their tumors might have actionable driver mutations that can provide treatment options otherwise not available for patients with these fatal cancers. This study investigated if comprehensive genomics analyses could inform on the origin of the tumor. PATIENT AND METHODS: Here we describe a patient whose tumor was misdiagnosed at least three times. Next-generation sequencing, a PDX mouse model and bioinformatics was used to identify an actionable mutation, predict resistance development to the targeted therapy, and to correctly diagnose the origin of the tumor. The Cancer Genome Atlas was used to benchmark the bioinformatics workflow. RESULTS: Despite the lack of a known primary tumor site and the absence of diagnostic immunohistochemical markers, the origin of the patient's tumor was established using the novel bioinformatics workflow. This included a mutational signature analysis of the sequenced metastases and comparison of their transcriptomic profiles to a pan-cancer panel of tumors from The Cancer Genome Atlas. We further discuss the strengths and limitations of the latter approaches in the context of three potentially incorrectly diagnosed TCGA lung tumors. CONCLUSION: Comprehensive genomics analyses could inform on the origin of tumors in patients suffering from CUP.

  Study EGAS00001003026

• Small variants in mtDNA Canary Islands - WES Illumina (ITER)

  The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

  Study EGAS00001005678

• Small variants in mtDNA Canary Islands - WGS Illumina (ITER)

  The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

  Study EGAS00001005679

• Non-coding mutations reveal cancer driver cistromes in luminal breast cancer

  Whole-genome sequencing of primary breast tumors enabled the identification of cancer driver genes and non-coding cancer driver plexuses from somatic mutations. However, differentiating between driver and passenger events among non-coding genetic variants remains a challenge to understand the etiology of cancer and inform the delivery of personalized cancer medicine. Herein, we reveal enrichment of non-coding mutations in cis-regulatory elements that cover a subset of transcription factors linked to tumor progression in luminal breast cancers. Using a cohort of 26 primary luminal ER+PR+ breast tumors, we compiled a catalogue of ~100,000 unique cis-regulatory elements from ATAC-seq data. Integrating this catalogue with somatic mutations from 350 publicly available breast tumor whole genomes, we identified four recurrently mutated individual cis-regulatory elements. By then partitioning the non-coding genome into cistromes, defined as the sum of binding sites for a transcription factor, we uncovered cancer driver cistromes for ten transcription factors, namely CTCF, ELF1, ESR1, FOSL2, FOXA1, FOXM1 GATA3, JUND, TFAP2A, and TFAP2C in luminal breast cancer. Nine of these ten transcription factors were shown to be essential for growth in breast cancer, with four exclusive to the luminal subtype. Collectively, we present a strategy to find cancer driver cistromes relying on quantifying the enrichment of non-coding mutations over cis-regulatory elements concatenated into a functional unit drawn from an accessible chromatin catalogue derived from primary cancer tissues.

  Study EGAS00001005235

• DNA methylation analysis on PBL obtained from male patients with UC-PSC, UC and HC

  Primary sclerosing cholangitis (PSC) is a chronic inflammatory liver disease affecting the intra- and extrahepatic bile ducts, and is strongly associated with ulcerative colitis (UC). In this study, we explored the peripheral blood DNA methylome and its immune cell composition in patients with PSC-UC, UC, and healthy controls (HC) with the aim to develop a predictive assay in distinguishing patients with PSC-UC from those with UC alone. The peripheral blood DNA methylome of male patients with PSC and concomitant UC, UC and HCs was profiled using the Illumina HumanMethylation Infinium EPIC BeadChip (850K) array. Differentially methylated CpG position (DMP) and region (DMR) analyses were performed alongside gradient boosting classification analyses to discern PSC-UC from UC patients. As observed differences in the DNA methylome could be the result of differences in cellular populations, we additionally employed mass cytometry to characterize the immune cell compositions. Genome wide methylation analysis revealed no differentially methylated positions between PSC-UC and UC patients. Nonetheless, using gradient boosting we were capable of discerning PSC-UC from UC with an area under the receiver operator curve (AUROC) of 0.80. Four CpG sites annotated to the NINJ2 gene were found to strongly contribute to the predictive performance. While mass cytometry analyses corroborated the largely similar blood cell composition among patients with PSC-UC, UC and HC, a higher abundance of myeloid cells was observed in UC compared to PSC-UC patients.

  Study EGAS00001005832

• Predictive value of chromosome 18q11.2-q12.1 loss for benefit from bevacizumab in metastatic colorectal cancer; a post-hoc analysis of the randomized controlled trial AGITG-MAX

  Patients and methods: We analyzed tumor samples of 256 patients enrolled in the AGITG-MAX 3-arm randomized controlled trial (RCT) who were randomly assigned to receive either capecitabine, or capecitabine and bevacizumab with- or without mitomycin. Chromosome 18q11.2-q12.1 copy number status was measured and correlated with PFS according to a predefined analysis plan with marker-treatment interaction as the primary end-point.Results: Chromosome 18q11.2-q12.1 losses were detected in 71% of patients (181/256) characteristic for mCRC. Consistent with the non-randomized study, significant PFS benefit of bevacizumab was observed in chromosome 18q11.2-q12.1 loss patients (P=0.009), and not in 18q no-loss patients (P=0.67). Although significance for marker-treatment interaction was not reached (Pinteraction = 0.28), hazard ratio and 95% confidence interval (HRinteraction 0.72, CI 0.39 - 1.32) show striking overlap with the non-randomized study cohorts (HRinteraction 0.41, CI 0.32 - 0.8) supported by a non-significant Cochrane chi-squared (P=0.11) for heterogeneity. Conclusions: PFS characteristics for marker positive- and negative patients randomized for treatment with- or without bevacizumab were consistent with those of the non-randomized discovery study, though the number of AGITG-MAX samples fell short to reach significance for marker-treatment interaction. Post-hoc analysis of the AGITG-MAX RCT provides supportive evidence for chromosome 18q11.2-q12.1 as a clinically relevant predictive marker for mCRC patients.

  Study EGAS00001005453

• A novel orthotopic patient-derived xenograft model of radiation-induced glioma following medulloblastoma

  Radiation-induced glioma (RIG) is a highly aggressive brain cancer arising as a consequence of radiation therapy. We report a case of RIG that arose in the brain stem following treatment for paediatric medulloblastoma, and the development and characterisation of a matched orthotopic patient-derived xenograft (PDX) model (TK-RIG915). Patient and PDX tumours were analysed using DNA methylation profiling, whole genome sequencing (WGS) and RNA sequencing. While initially thought to be a diffuse intrinsic pontine glioma (DIPG) based on disease location, results from methylation profiling and WGS were not consistent with this diagnosis. Furthermore, clustering analyses based on RNA expression suggested the tumours were distinct from primary DIPG. Additional gene expression analysis demonstrated concordance with a published RIG expression profile. Multiple genetic alterations that enhance PI3K/AKT and Ras/Raf/MEK/ERK signalling were discovered in TK-RIG915 including an activating mutation in PIK3CA, upregulation of PDGFRA and AKT2, inactivating mutations in NF1, and a gain-of-function mutation in PTPN11. Additionally, deletion of CDKN2A/B, increased IDH1 expression, and decreased ARID1A expression were observed. Detection of phospho-S6, -4EBP1 and -ERK via immunohistochemistry confirmed PI3K pathway and ERK activation. Here, we report one of the first PDX models for RIG, which recapitulates the patient disease and is molecularly distinct from primary brain stem glioma. Genetic interrogation of this model has enabled the identification of potential therapeutic vulnerabilities in this currently incurable disease.

  Study EGAS00001004709

• Small variants in mtDNA Canary Islands - WGS Oxford Nanopore Technologies (ITER)

  The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

  Study EGAS00001005677

• A Dose Escalation Study of Efmarodocokin Alfa (UTTR1147A) in Healthy Volunteers and Patients with Ulcerative Colitis

  Background: The interleukin-22 cytokine (IL-22) has demonstrated efficacy in nonclinical colitis models with a non-immunosuppressive mechanism-of-action. Efmarodocokin alfa (UTTR1147A) is a fusion protein agonist that links IL-22 to the crystallizable fragment (Fc) of human immunoglobulin (IgG4) for improved pharmacokinetic characteristics, but with a mutation to minimize Fc effector functions. Methods: This randomized, phase 1b study evaluated the safety, tolerability, pharmacokinetics, and pharmacodynamics of repeat intravenous dosing of efmarodocokin alfa in healthy volunteers (HVs; n=32) and patients with ulcerative colitis (UC; n=24) at 30–90 µg/kg doses given biweekly (Q2W) or monthly (Q4W) (6:2 active:placebo per cohort). Results: The most common adverse events (AEs) were on-target dermatological effects (dry skin, erythema, and pruritus) that were reversible. Dose-limiting non-serious dermatological AEs (severe dry skin, erythema, exfoliation, and discomfort) were seen in two HVs and one patient at 90 ug/kg Q2W. Pharmacokinetics were generally dose-proportional across the dose levels tested, but patients demonstrated lower drug exposures relative to HVs at the same dose. IL-22 serum biomarkers and IL-22 responsive genes in colon biopsies were induced with active treatment. Patients demonstrated changes in microbiota composition following active treatment, thereby reversing baseline dysbiosis. Clinical response was observed in 7/18 active- and 1/6 placebo-treated patients; clinical remission was observed in 5/18 active- and 0/6 placebo-treated patients.

  Study EGAS00001006172

• Single_cell_characterization_of_T_cell_lymphoma_

  Lymphomas are types of blood cancer derived from lymphocytes, and can be classified into two main categories, Hodgkin's lymphomas (HL) and non-Hodgkin lymphomas (NHL). T-cell lymphoma is a type of NHL that develops from T lymphocytes and accounts for 10-15% of total NHL cases. As a rare disease, T-cell lymphoma is functionally, pathologically, and clinically complex and heterogeneous. The 2016 revision of the WHO classification has listed 29 subtypes of mature T/NK-cell neoplasms, among which over 25 are T lymphocyte in origin. These subtypes of T-cell lymphomas have different biological behaviours and clinical prognosis. Clinically, T-cell lymphoma ranges from indolent (slow-growing) to aggressive (fast-growing and spreading) disease. According to the guideline booklet provided by the Lymphoma Research Foundation (https://lymphoma.org/aboutlymphoma/nhl/), indolent T-cell lymphoma subtypes include adult T-cell leukaemia/lymphoma (ATLL), cutaneous T-cell lymphoma (CTCL), and mycosis fungoides (MF), while the aggressive subtypes include anaplastic large cell lymphoma (ALCL), angioimmunoblastic T-cell lymphoma (AITL), peripheral T-cell lymphoma (PTCL), and T lymphoblastic lymphoma (TLBL). PTCL and ALCL represent the most common subtypes of aggressive T-cell lymphoma whereas TLBL is a rare and aggressive subtype mostly diagnosed in children. In this study, we will use single cell RNA sequencing and single cell multiomic methods three aggressive T-cell lymphoma subtypes (ALCL, PTCL and TLBL) and try to gain a deeper insight into the genotypic and phenotypic features of these T-cell lymphomas.

  Study EGAS00001005750

• Genomic stratification and liquid biopsy in a rare adrenocortical carcinoma (ACC) case, with dual lung metastases

  Adrenocortical carcinoma is a rare malignancy with a poor prognosis and few treatment options. Molecular characterization of this cancer remains limited. We present a case of an adrenocortical carcinoma (ACC) in a 37-yr-old female, with dual lung metastases identified 1 yr following commencement of adjuvant mitotane therapy. As standard therapeutic regimens are often unsuccessful in ACC, we undertook a comprehensive genomic study into this case to identify treatment options and monitor disease progress. We performed targeted and whole-genome sequencing of germline, primary tumor, and both metastatic tumors from this patient and monitored recurrence over 2 yr using liquid biopsy for ctDNA and steroid hormone measurements. Sequencing revealed the primary and metastatic tumors were hyperhaploid, with extensive loss of heterozygosity but few structural rearrangements. Loss-of-function mutations were identified in MSH2, TP53, RB1, and PTEN, resulting in tumors with mismatch repair signatures and microsatellite instability. At the cellular level, tumors were populated by mitochondria-rich oncocytes. Longitudinal ctDNA mutation and hormone profiles were unable to detect micrometastatic disease, consistent with clinical indicators of disease remission. The molecular signatures in our ACC case suggested immunotherapy in the event of disease progression; however, the patient remains free of cancer. The extensive molecular analysis presented here could be applied to other rare and/or poorly stratified cancers, to identify novel or repurpose existing therapeutic options, thereby broadly improving diagnoses, treatments, and prognoses.

  Study EGAS00001003497

• Conserved features of TERT promoter duplications reveal an activation mechanism that mimics hotspot mutations in cancer

  Mutations in the TERT promoter are the single most common non-coding mutation in cancer and represent the genetic underpinnings of tumor cell immortality. Beyond the two most common point mutations, G228A and G250A, which selectively recruit the ETS factor GABP to activate TERT, the significance of other variants in the TERT promoter are unknown. We identified duplications of wildtype sequence within the core promoter region of TERT in 7 different cancer types that have strikingly similar features including size, insertion position, and inclusion of one of the native ETS motifs. Each duplication activates the TERT promoter to a similar level as G228A and G250A and is critically dependent on the insertion site. The GABP tetramer binds to the TERT duplicated promoter sequence by virtue of the native ETS motif and its duplicated version with precise spacing, and it is necessary for the transcriptional activation by all duplications tested. Spatiotemporal analysis in a multifocal glioblastoma shows the duplication is clonal and its activation of TERT is readily detectable at the single cell level and in bulk tumor tissue. We conclude that recurrent TERT promoter duplications of the native ETS sequence are functionally and mechanistically equivalent to the hotspot mutations that confer tumor cell immortality. The shared mechanism of these divergent somatic genetic alterations suggests a strong selective pressure for recruitment of the GABP tetramer to activate TERT.

  Study EGAS00001006118

• Prostate cancer ancestral genomic disparity

  Prostate cancer is characterised by significant global disparity; mortality rates in SubSaharan Africa are double to quadruple those in Eurasia. Hypothesising unknown interplay between genetic and non-genetic factors, tumour genome profiling envisages contributing mutational processes. Through whole-genome sequencing of treatment-naïve prostate cancer from 183 ethnically/globally distinct patients (African versus European), we generate the largest cancer genomics resource for Sub-Saharan Africa. Identifying ~2 million somatic variants, Africans carried the greatest burden. We describe a new molecular taxonomy using all mutational types and ethno39 geographic identifiers, including Asian. Defined as Global Mutational Subtypes (GMS) A–D, although Africans presented within all subtypes, we found GMS-B to be ‘African-specific’ and GMS-D ‘African-predominant’, including Admixed and European Africans. Conversely, Europeans from Australia, Africa and Brazil predominated within ‘mutationally-quiet’ and ethnically/globally ‘universal’ GMS-A, while European Australians shared a higher mutational burden with Africans in GMS45 C. GMS predicts clinical outcomes; reconstructing cancer timelines suggests four evolutionary trajectories with different mutation rates (GMS-A, low 0.968/year versus D, highest 1.315/year). Our data suggest both common genetic factors across extant populations and regional environmental factors contributing to carcinogenesis, analogous to gene-environment interaction defined here as a different effect of an environmental surrounding in persons with different ancestries or vice versa. We anticipate GMS acting as a proxy to intrinsic and extrinsic mutational processes in cancers, promoting global inclusion in landmark studies.

  Study EGAS00001006425

• Disease and phenotype relevant genetic variants identified from histone acetylomes in human hearts

  Identifying genetic markers for heterogeneous complex diseases such as heart failure has been challenging, and may require prohibitively large cohort sizes in genome-wide association studies (GWAS) in order to demonstrate statistical significance1. On the other hand, chromatin quantitative trait loci (QTL), elucidated by direct epigenetic profiling of specific human tissues, may contribute towards prioritising variants for disease-association. Here, we captured non-coding genetic variants by performing enhancer H3K27ac ChIP-seq in 70 human control and end-stage failing hearts, mapping out a comprehensive catalogue of 47,321 putative human heart enhancers. 3,897 differential acetylation peaks (FDR < 0.05) pointed to recognizable pathways altered in heart failure (HF). To identify cardiac histone acetylation QTLs (haQTLs), we regressed out confounding factors including HF disease status, and employed the G-SCI test2 to call out 1,680 haQTLs (FDR < 0.1). A subset of these showed significant association to gene expression, either in cis (180), or through long range interactions (81), identified by Hi-C and Hi-ChIP performed on a subset of hearts. Furthermore, a concordant relationship was found between the gain or disruption of specific transcription factor (TF) binding motifs, inferred from alternative alleles at the haQTLs, associated with altered H3K27ac peak heights. Finally, colocalisation of our haQTLs with heart-related GWAS datasets allowed us to identify 62 unique loci. Disease-association for these new loci may indeed be mediated through modification of H3K27-acetylation enrichment and their corresponding gene expression differences.

  Study EGAS00001003586

• Single-cell ATAC-seq analysis for COVID19 patients

  While SARS-CoV-2 infection causes mild respiratory disease in most individuals, a small group of patients develops severe COVID-19. Dysfunctional innate immune responses have been identified to contribute to differences in COVID-19 severity, but the key regulators are still unknown. Here, we present an integrative single-cell epigenetics, transcriptomic, and genetics analysis of peripheral blood mononuclear cells from hospitalized and convalescent COVID-19 patients. In classical monocytes, we identified 41.3% of significantly up-regulated genes in hospitalized COVID-19 patients potentially induced by differential chromatin accessibility. Sub-clustering and motif-enrichment analyses of monocytes reveal disease condition-specific regulation by transcription factors, such as C/EBPs and SPI1, and their targets, including a long-noncoding RNA LUCAT1, which further regulates interferon responses and is associated with the need for oxygen supply of COVID-19 patients. The interaction between C/EBPs and LUCAT1 was validated through loss-of-function experiments. Finally, we investigated genetic risk variants that exhibit allele-specific open chromatin (ASoC) in promoters/enhancers of COVID-19 patients. Integrating our data with publicly available expression quantitative trait loci and chromosomal interactions indicates that ASoC SNP rs6800484-C is associated with lower expression of CCR2, which may contributeto higher viral loads in lungs and higher risk of COVID-19 hospitalization. Altogether, our study highlights the diverse genetic and epigenetic regulators that contribute to the innate immune responses of different COVID-19 patients.

  Study EGAS00001006559

• Single-cell RNA-seq analysis for COVID19 patients

  While SARS-CoV-2 infection causes mild respiratory disease in most individuals, a small group of patients develops severe COVID-19. Dysfunctional innate immune responses have been identified to contribute to differences in COVID-19 severity, but the key regulators are still unknown. Here, we present an integrative single-cell epigenetics, transcriptomic, and genetics analysis of peripheral blood mononuclear cells from hospitalized and convalescent COVID-19 patients. In classical monocytes, we identified 41.3% of significantly up-regulated genes in hospitalized COVID-19 patients potentially induced by differential chromatin accessibility. Sub-clustering and motif-enrichment analyses of monocytes reveal disease condition-specific regulation by transcription factors, such as C/EBPs and SPI1, and their targets, including a long-noncoding RNA LUCAT1, which further regulates interferon responses and is associated with the need for oxygen supply of COVID-19 patients. The interaction between C/EBPs and LUCAT1 was validated through loss-of-function experiments. Finally, we investigated genetic risk variants that exhibit allele-specific open chromatin (ASoC) in promoters/enhancers of COVID-19 patients. Integrating our data with publicly available expression quantitative trait loci and chromosomal interactions indicates that ASoC SNP rs6800484-C is associated with lower expression of CCR2, which may contributeto higher viral loads in lungs and higher risk of COVID-19 hospitalization. Altogether, our study highlights the diverse genetic and epigenetic regulators that contribute to the innate immune responses of different COVID-19 patients.

  Study EGAS00001006560

• Direct Comparative Analysis of 10X Genomics Chromium and Smart-seq2

  Single cell RNA sequencing (scRNA-seq) is widely used for profiling transcriptomes of individual cells. The droplet-based 10X Genomics Chromium (10X) approach and the plate-based Smart-seq2 full-length method are two frequently-used scRNA-seq platforms, yet there are only a few thorough and systematic comparisons of their advantages and limitations. Here, by directly comparing the scRNA-seq data by the two platforms from the same samples of CD45- cells, we systematically evaluated their features using a wide spectrum of analysis. Smart-seq2 detected more genes in a cell, especially low abundance transcripts as well as alternatively spliced transcripts, but captured higher proportion of mitochondrial genes. The composite of Smart-seq2 data also resembled bulk RNA-seq data better. For 10X-based data, we observed higher noise for mRNA in the low expression level. Despite the poly(A) enrichment, approximately 10-30% of all detected transcripts by both platforms were from non-coding genes, with lncRNA accounting for a higher proportion in 10X. 10X-based data displayed more severe dropout problem, especially for genes with lower expression levels. However, 10X-data can better detect rare cell types given its ability to cover a large number of cells. In addition, each platform detected different sets of differentially expressed genes between cell clusters, indicating the complementary nature of these technologies. Our comprehensive benchmark analysis offers the basis for selecting the optimal scRNA-seq strategy based on the objectives of each study.

  Study EGAS00001003973

• Whole-genome sequencing of normal Singaporean volunteers

  Sleep is associated with various health outcomes. Despite their growing adoption, the potential for consumer wearables to contribute sleep metrics to sleep-related biomedical research remains largely uncharacterized. Here we analyzed sleep tracking data, along with questionnaire responses and multi-modal phenotypic data generated from 482 normal volunteers. First, we compared wearable-derived and self-reported sleep metrics, particularly total sleep time (TST) and sleep efficiency (SE). We then identified demographic, socioeconomic and lifestyle factors associated with wearable-derived TST. Among others, male gender (β = -15.539, 95% confidence interval [CI] = -26.245 - -4.832, p = 0.005), older age (β = -0.493, CI = -0.941 - -0.044, p = 0.032) and manual labor (β = -26.856, CI = -49.715 - -3.997, p = 0.022) were associated with reduced TST, whereas alcohol consumption (β = 19.247, CI = 8.008 - 30.486, p = 8.54x10-04) was associated with increased TST. Multi-modal phenotypic data analysis showed that wearable-derived TST and SE were associated with various cardiovascular disease risk markers such as body mass index (β = -0.006, CI = -0.00 - -0.000, p = 0.040) and waist circumference (β = 0.001, CI = -0.015-0.017, p = 0.893), whereas self-reported measures were not. We then showed that insufficient sleep was associated with premature telomere attrition, using telomere length estimated using whole-genome sequencing (β = 1.275, CI = 0.187 - 2.363, p = 0.023) and quantitative PCR (β = 0.001, CI = 8.318x10-05 - 0.001, p = 0.028). Our study highlights the potential for sleep metrics from consumer wearables to provide novel insights into data generated from population cohort studies.

  Study EGAS00001004007

• Whole transcriptome profiling of liquid biopsies from tumour xenografted mouse models: validation cohort

  While cell-free DNA (cfDNA) in liquid biopsies is widely being used and investigated, free circulating RNA (extracellular RNA, exRNA) has the potential to improve therapy response monitoring and cancer detection due to its dynamic nature. A fundamental open question hampering the initiation of large-scale liquid biopsy collections for tumour exRNA analysis is that it remains unclear in which blood subcompartment tumour-derived exRNAs primarily reside. We set out to develop a host-xenograft deconvolution framework, exRNAxeno, with cDNA mapping strategies to either a combined human-mouse reference genome or both species genomes in parallel, that can be applied to exRNA sequencing data from liquid biopsies of human xenograft mouse models, enabling to distinguish (human) tumoural RNA from (murine) host RNA, and as such to specifically analyse tumour-derived exRNA. Subsequently, the preferred exRNAxeno combination pipeline was applied to total exRNA sequencing data from blood platelets and three plasma fractions from a breast cancer patient-derived and neuroblastoma cell line-derived xenograft mouse model. We show that tumoural exRNA concentrations are not determined by plasma platelet levels, while host exRNA concentrations increase with platelet content. Furthermore, a large variability in exRNA levels and gene content across individual mice is observed. In general, the tumoural gene detectability in plasma is correlated with the RNA expression levels in the tumour tissue or cancer cell line. Our results unravel new aspects of tumour-derived exRNA biology in xenograft mouse models.

  Study EGAS00001006582

• H3Africa - Stroke Investigative Research and Education Networks

  Although stroke is a discrete phenotype, it is the clinical culmination of several complex and interacting biological processes, precipitated by various genetic and environmental factors, thereby making ready analyses of its underlying mechanisms a challenge. However, the formulation of successful tailor-made prevention strategies for limiting the immense personal and societal burden of stroke in Sub-Saharan Africa (SSA) is contingent on a better understanding of the predisposing risk factors for stroke in the region, which may not necessarily match the nature of stroke risk factors described among other people on other continents. Indeed, the underlying risk factors for stroke in SSA have hitherto not been comprehensively characterized, and currently there is a paucity of expertise, experience, and infrastructure in genomics and environmental research in SSA, to initiate and sustain such investigation. We seek to establish the Stroke Investigative Research and Education Network (SIREN) as a multidisciplinary collaborative research network of investigators in SSA and the United States collectively focused on exploring ways to promote a better understanding of the genetic and environmental risk factors for stroke among people of African ancestry. SIREN will evaluate unique contributors to stroke in SSA and compare these findings with those in a cohort of African American stroke survivors, all while simultaneously building sustainable capacities in phenomics, biobanking, genomics, biostatistics and bioinformatics for future high-level investigation of stroke and other vascular disease entities in SSA.

  Study EGAS00001007331

• Atypical B cells and impaired SARS-CoV-2 neutralisation following heterologous vaccination in the elderly

  Age is a major risk factor for hospitalization and death after SARS-CoV-2 infection, even in vaccinees. Suboptimal responses to a primary vaccination course have been reported in the elderly, but there is little information regarding the impact of age on responses to booster third doses. Here we show that individuals 70 or older who received a primary two dose schedule with AZD1222 and booster third dose with mRNA vaccine achieved significantly lower neutralizing antibody responses against SARS-CoV-2 spike pseudotyped virus compared to those below 70 at one month post booster. Neither the concentration of serum binding anti spike IgG antibody, nor the frequency of spike- specific B cells showed differences by age grouping. Impaired neutralization potency and breadth post-third dose in the elderly was associated with enrichment of circulating “atypical” spike-specific B cells expressing CD11c and FCRL5. Presence of these atypical B cells was confirmed by scRNA seq. However, when considering individuals who received three doses of mRNA vaccine, we did not observe statistically significant differences in neutralization or enrichment in atypical B cells by age grouping following dose 3. Together, our data reveal impaired neutralising antibody responses in the elderly after heterologous vaccination with AdV followed by mRNA booster dose, but not in homologous three dose mRNA vaccination highlighting that different vaccine formats differentially instruct the memory B cell response.

  Study EGAS00001007385

• Oxidative phosphorylation is a key ontogenetic feature of monocyte immunometabolism promoting myeloid differentiation after birth

  Neonates primarily rely on innate immunity, yet their inflammatory responses to microbes, particularly those of monocytes, are usually restricted compared to adults. This is controversially interpreted as immaturity increasing the risk of sepsis or essential programming allowing immune adaptation after birth. Which changes the cellular immunometabolism undergoes after birth and its role in these concepts are poorly defined. Here, we applied transcriptomic, metabolic, and immunological approaches and found that monocytes exhibit an inverse ontogenetic balance of immunometabolism. Our data show that glycolysis in monocytes increases within the first year of life and fuels their inflammatory responsiveness. In contrast, oxidative phosphorylation (OXPHOS) is high in neonatal monocytes supporting myeloid differentiation but declines only gradually during the first five years of life. Treatment of neonatal monocytes with lipopolysaccharide induces an adult-like immunometabolic phenotype. Ketogenic diet restricts glycolysis in adult monocytes, but cannot reactivate OXPHOS and revive a neonatal-like immunometabolic phenotype, suggesting that neonatal metabolism is hardwired and cannot be restored by simple dietary changes. Transcriptional network and population-wide human variation analyses identified E2F1, MYB, STAT1 and FLI1 as important regulators of age-dependent cell functions and related energy demands in human monocytes. Collectively, our findings show that restricted glycolysis and increased OXPHOS are a physiological programming in neonatal monocytes. Premature switching to an adult-like metabolism could untimely enhance inflammatory responses and disrupt important myeloid differentiation, whereas microbial challenges accumulating during childhood seem to induce an essential gradual metabolic switch.

  Study EGAS00001007555

• Genome-wide association study identifies common variants associated with breast cancer in South African Black women

  Genome-wide association studies (GWAS) have characterized the contribution of common variants to the development of breast cancer (BC) in populations of European descent, however no GWAS have been reported in resident African populations. This GWAS included 2,485 resident African BC cases and 1,101 population matched controls. Two risk loci were identified, one on chromosome 15, located between UNC13C and RAB27A (rs7181788, p = 1.01x10-08) and a second in an intron of USP22 on chromosome 17 (rs899342, p = 4.62x10-08). RAB27A is a member of the RAS oncogene family associated with invasive growth and metastasis, and UNC13C is involved in tumor progression. The rs899342 SNP is an eQTL for expression of USP22 in a wide range of tissues. Several genome-wide significant signals were detected in a hormone receptor subtype analysis. Comparison of a meta-analysis of other African ancestry BC GWAS datasets predominantly of West African ancestry with our South African data showed limited evidence of shared risk loci, suggesting potential heterogeneity in genetic risk factors for breast cancer in different African populations. A European ancestry derived polygenic risk model for BC was applied to our data, explaining only 0.79% of variance, with an area under the curve of 0.56. This study identified genetic risk factors for BC in a resident African population, but larger studies in this and other African populations are needed to further define the contribution of common variants to BC risk.

  Study EGAS00001008032

• WGS and WES data of invasive micropapillary carcinoma of breast

  This study contains whole genome sequencing data and whole exon sequencing data of IMPC tumor and normal tissue sample.

  Dataset EGAD00001008450

• cfMethyl-seq data (cfSort study) from serial plasma samples of NSCLC patients

  The cfMethyl-Seq libraries of the serial plasma cfDNA samples from the four NSCLC patients were constructed following the standard protocol. 10 ng of cfDNA in the volume of 25 µl was used as input material. 5’-end dephosphorylation was done with 3 µl 10xCutSmart buffer and 2 µl quick CIP from NEB (Ipswich, MA) at 37 oC for 30 min then heat-inactivated at 80 oC for 5 min. The 3’-end blocking was done with 0.5 µl 10xCutSmart buffer, 3 µl 2.5 mM CoCl2, 1 µl terminal transferase (all from NEB), and 0.5 µl 1 mM ddGTP at 37 oC for 2 h followed by 75 oC for 20 min. The mixture was then purified with 2x AmpureXP beads (Beckman Coulter, Indianapolis, IN) and eluted in 21.5 µl RT-PCR grade water (Thermo-Fisher, Waltham, MA). Restriction digestion was done with 2.5 µl 10xCutSmart buffer and 1 µl MspI (NEB) for 18 h at 37  oC and 20 min at 65 oC . 0.5 µl 10xCutSmart buffer, 0.3 µl dACGTP mixture (100 mM dATP, 10 mM dCTP, 10 mM dGTP), 1 µl Klenow (exo-, 5U/µl, NEB) and 2.6 µl RT-PCR water, 0.6 µl 50 mM DTT (ThermoFisher) was added to the mixture for end repair and A-overhang addition with the program 30  oC for 20 min, 37 oC for 1 h and 75 oC for 20 min. Adapter ligation was then performed with 1 µl 10xThermoFisher HC T4 ligase buffer, 0.4 µl 100 mM ATP (ThermoFisher), 0.2 µl 50 mM DTT, 1 µl ThermoFisher HC T4 DNA ligase (30 Weiss Unit/µl), 5 ng home-made duplex UMI adapter with all the cytosines methylated (protocol adopted from Kennedy et al.) at 16 oC for 20 h and 65 oC for 20 min. Bisulfite conversion of the adapter-ligated product was carried out with QIAGEN EpiTect plus DNA bisulfite kit following their protocol for two rounds of conversion. The converted product was purified with Qiagen MinElute spin column and eluted with 20 µl RT-PCR water. PCR amplification was done using the NEBNext Multiplex Oligos for Illumina (2.5 µl of universal and index primer each) and 25 µl KAPA HiFi HotStart Uracil+ ReadyMix (Roche) with the following cycling conditions: 98  oC for 45 s, 15 cycles of 98  oC for 15 s, 60  oC for 30 s and 72  oC for 30 s, followed by a final extension at 72  oC for 5 min. The PCR product was purified with 1x AmpureXP beads and eluted with 30 µl EB buffer. DNA concentration was measured by Qubit 1xdsDNA HS assay. 5% TBE-UREA PAGE and bioanalyzer assay was performed as quality control on each library before sequencing.

  Dataset EGAD00001010881

• The long term effects of chemotherapy on normal blood - WGS dataset

  In developed countries, ~10% of individuals are exposed to systemic chemotherapy for cancer and other diseases. Many chemotherapeutic agents act by increasing DNA damage in cancer cells, hence triggering cell death. However, there is limited understanding of the extent and consequences of collateral DNA damage to normal tissues. To investigate the impact of chemotherapy on mutation burdens and cell population structure of a normal tissue we sequenced blood cell genomes from 23 individuals, aged 3-80 years, treated with a range of chemotherapy regimens. Substantial additional mutation loads with characteristic mutational signatures were imposed by some chemotherapeutic agents, but there were differences in burden between different classes of agent, different agents of the same class and different blood cell types. Chemotherapy also induced premature changes in the cell population structure of normal blood, similar to those of normal ageing. The results constitute an initial survey of the long-term biological consequences of cytotoxic agents to which a substantial fraction of the population is exposed during the course of their disease management, raising mechanistic questions and highlighting opportunities for mitigation of adverse effects.

  Dataset EGAD00001015339

• NIPT samples for systematic evaluation of BinDel, a software tool for detecting clinically significant microdeletions in low-coverage WGS-based NIPT samples

  This dataset contains three sets of samples. The first sample set contains euploid fetus pregnancies reported by NIPTIFY screening test and postnatal evaluation. Dataset was processed similarly to previously published guidelines from KU Leuven, with modifications [1]. Briefly, peripheral blood samples were collected in cell-free DNA BCT tubes (Streck, USA), and plasma was separated with standard dual centrifugation. Cell-free DNA was extracted from 3 ml plasma using MagMAX Cell-Free DNA Isolation Kit (ThermoFisher Scientific). Whole-genome libraries were prepared using the FOCUS (Fragmented DNA Compact Sequencing Assay, Competence Centre on Health Technologies, Estonia) NIPT method protocol with 12 cycles for the final PCR enrichment step. In the following quantification, equal amounts of 36 samples were pooled, and the quality and quantity of the pool were assessed on Agilent 2200 TapeStation (Agilent Technologies, USA). Whole genome sequencing was performed on the NextSeq 550 instrument (Illumina Inc.) with an average coverage of 0.32× (minimum 0.08 and maximum 0.42) and producing 85 bp single-end reads. The second sample set contains a single NIPT sample postnatally diagnosed with Prader-Willi syndrome. The sample was sequenced with Illumina NextSeq 500 platform, producing 85 bp single-end reads with an average per-sample coverage of 0.32× at the University of Tartu, Institute of Genomics Core Facility, according to the manufacturer’s standard protocols, as described previously [2]. The third sample set contains samples SC005 (SeraCare Life Sciences Inc lot #10446565), SC0042 (#10571706), and SC016 (#10560229). These are SeraCare Life Sciences Inc circulating cell-free DNA (ccfDNA) like mixture of human genomic DNA that consists of matched maternal and fetus. SC005 and SC0042 consist of matched DNA of maternal and fetus with DiGeorge Syndrome. SC016 is a custom-ordered DNA Mix with fetus DNA having a pathogenic loss of the terminal region of 20p13 and a pathogenic 3q29 duplication. SC016 was processed as the first sample set was processed, and SC0042 was processed as the second sample set was processed. Sample SC005 was processed once as was sample set 1 and once as was sample set 2 processed. This study was performed with the approval of the Research Ethics Committee of the University of Tartu (#352/M-12). 1. Bayindir B, Dehaspe L, Brison N, Brady P, Ardui S, Kammoun M, et al. Noninvasive prenatal testing using a novel analysis pipeline to screen for all autosomal fetal aneuploidies improves pregnancy management. Eur J Hum Genet. 2015;23: 1286– 1293. doi:10.1038/ejhg.2014.282 2. Žilina O, Rekker K, Kaplinski L, Sauk M, Paluoja P, Teder H, et al. Creating basis for introducing noninvasive prenatal testing in the Estonian public health setting. Prenat Diagn. 2019;39: 1262–1268. doi:10.1002/pd.5578

  Dataset EGAD00001009512

• 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

• Ultra-long sequencing for contiguous haplotype resolution of the human immunoglobulin heavy chain locus

  Genetic diversity within the human immunoglobulin heavy chain (IGH) locus influences the expressed antibody repertoire and susceptibility to infectious and autoimmune diseases. However, repetitive sequences and complex structural variation pose significant challenges for large-scale characterization. Here, we introduce a method that combines Oxford Nanopore ultra-long sequencing and adaptive sampling with a bioinformatic pipeline to produce haplotype-resolved, annotated IGH assemblies. Notably, our strategy overcomes prior limitations in phasing resolution, enabling single-contig haplotype assemblies that span the entire IGH locus. We applied this method to four individuals and validated the accuracy of the IGH assemblies using Pacific Biosciences HiFi reads, demonstrating near-complete sequence congruence, detecting only indels with a high degree of confidence. Moreover, when applied to the reference material HG002, our pipeline revealed no base differences and a limited number of indels compared with the Telomere-to-Telomere genome benchmark across the IGH region. Importantly, in the four individuals, our approach uncovered 30 novel alleles and previously uncharacterized large structural variants, including a 120 kb duplication spanning IGHE to IGHA1 and an expanded seven-copy IGHV3-23 gene haplotype. These findings underscore the power of our method to resolve the full complexity of the IGH locus and uncover previously unrecognized variants that may affect immune function and disease susceptibility. Thus, our method provides a strong basis for future immunological research and translational applications.

  Study EGAS50000001042

• Evolution of an adenomacarcinoma in response to selection by targeted kinase inhibitors

  Background - Adenocarcinomas of the tongue are rare and represent the minority (20-25%) of salivary gland tumors affecting the tongue. We investigated the utility of massively parallel sequencing to characterize an adenocarcinoma of the tongue, before and after treatment. Results - In the pre-treatment tumor we identified 7,629 genes within regions of copy number gain. 1,078 genes exhibited increased expression relative to the blood and unrelated tumors and four genes contained somatic protein-coding mutations. Our analysis suggested the tumor cells were driven by the RET oncogene. Genes whose protein products are targeted by the RET inhibitors sunitinib and sorafenib correlated with being amplified and or highly expressed. Consistent with our observations administration of sunitinib was associated with stable disease lasting 4 months, after which the lung lesions began to grow. Administration of sorafenib and sulindac provided disease stabilization for an additional 3 months after which the cancer progressed and new lesions appeared. A recurring metastasis possessed 7,288 genes within copy number amplicons, 385 genes exhibiting increased expression relative to other tumours and 9 new somatic protein coding mutations. The observed mutations and amplif ications were consistent with therapeutic resistance arising through activation of the MAPK and AKT pathways. Conclusions - We conclude that complete genomic characterization of a rare tumour has the potential to aid in clinical decision making and identifying therapeutic approaches where no established treatment protocols exist. These results also provide direct in-vivo genomic evidence for mutational evolution within a tumor under drug selection and potential mechanisms of drug resistance accrual.

  Study EGAS00000000074

• Further details about data access are available at KNIH DAC.

  Dac EGAC00001000487

• Pharmacological improvement of CFTR function rescues airway epithelial homeostasis and host defense in cystic fibrosis children

  Rationale: Recent clinical trials have shown that elexacaftor/tezacaftor/ivacaftor (ETI) provides significant benefit to patients with cystic fibrosis (CF) by improving CFTR function. Despite these promising clinical findings, the specific effects on the airway inflammatory status leading to structural damage and impaired lung function, remain elusive. Objectives: We investigated the transcriptional changes at the single cell level of airway epithelial and immune cells that underlie the clinical improvement by ETI therapy in children with CF. Methods: Nasal swabs from thirteen children with CF and at least one F508del allele aged 6 to 12 years were collected at baseline and three months after initiation of ETI and subjected to scRNA-seq. In addition, we collected nasal swabs from 12 age- and sex-matched controls. Sweat chloride concentrations, spirometry and multiple breath washout were used for clinical evaluation of CF children. Measurements and main results: ETI significantly reduced sweat chloride concentrations (-52.2 ± 14.3 mmol/L, P < 0.001) and improved the lung clearance index (-1.3 ± 1.8, P < 0.05). Single cell transcriptomics revealed an impaired interferon signalling in epithelial cells of CF children at baseline, which was partially restored by ETI in conjunction with an ETI-induced increase in expression of MHC I and II encoding genes similar to control levels. Additionally, ETI markedly reduced the inflammatory phenotype of immune cells, particularly of neutrophils and macrophages. Conclusions: Improvement of CFTR function by ETI restores epithelial homeostasis and reduces immune cell inflammatory responses in the upper airways of children with CF, highlighting the potential of early initiation of ETI therapy.

  Study EGAS50000000128

• Systemic Lupus Erythematosus Serum Stimulation of Human Intestinal Organoids Induces Barrier Leakiness and Changes in Goblet Cell Differentiation

  Human intestinal epithelial cells are the interface between potentially harmful luminal content and basally residing immune cells. Their role is not only nutrient absorption but also the formation of a tight monolayer that constantly secretes mucus creating a multi-layered protective barrier. Alterations in this barrier can lead to increased gut permeability which is frequently seen in individuals with chronic extraintestinal autoimmune diseases, such as Systemic Lupus Erythematosus (SLE). Despite recent advances in identifying alterations in gut microbiota composition in SLE patients, not much attention has been given to the epithelial barrier itself. To date, it remains largely unexplored which role and function intestinal epithelial cells have in SLE pathology. Here, we present a unique near-physiologic in vitro model specifically designed to examine the effects of SLE on the epithelial cells. We utilize human colon organoids that are stimulated with serum obtained from SLE patients. Combining bulk and single-cell RNA transcriptomic analysis with functional assays revealed that SLE serum stimulation induced a unique expression profile marked by a reduction of goblet cell marker expression and mucus composition. Additionally, organoids exhibited imbalanced cellulat composition along with enhanced paracellular permeability,mitochondrial function, and a type I interferon gene signature. Similarly, transcriptomic analysis of SLE human colon biopsies revealed a downregulation of epithelial secretory markers. Our work uncovers a crucial connection between SLE and intestinal homeostasis that might be promoted in vivo through the blood, offering insights into the causal connection of barrier dysfunction and autoimmune diseases.

  Study EGAS50000000012

• Carboxylesterase 1 mediates a distinctive metabolic profile of dendritic cells to attain an inflammatory phenotype

  The metabolic profile of dendritic cells (DCs) shapes their phenotype and functions. Carboxylestrase 1 (CES1) enzyme is highly expressed in mononuclear myeloid cells however its exact role in DCs is elusive. We used a CES1 inhibitor (WWL113) and genetic overexpression to explore the role of CES1 in DCs differentiation in inflammatory models. CES1 expression was analyzed during CD14+ monocytes differentiation to DCs (MoDCs) using quantitative PCR. CES1 Inhibitor (WWL113) was applied during MoDCs differentiation. Surface markers, secreted cytokines, lactic acid production, phagocytic and T cell polarization capacity were analyzed. Transcriptomic and metabolic profile were assessed with RNA-sequencing and mass spectrometry. Cellular respiration was assessed with seahorse respirometry. Transgenic mice were used to assess CES1 overexpression in DCs in inflammatory models. CES1 expression peaks early during MoDCs differentiation. Pharmacological inhibition of CES1 led to higher expression of CD209, CD86 and MHCII. WWL113 treated MoDCs secreted higher quantities of IL6, IL8, TNF and IL10 and demonstrated stronger phagocytic ability and higher capacity to polarize Th17 differentiation in autologous DCs-T cells co-culture model. Transcriptomic profiling revealed enrichment of multiple inflammatory and metabolic pathways. Functional metabolic analysis shows impaired maximal mitochondrial respiration capacity, increased lactate production and decreased intracellular amino acids and TCA intermediates. Transgenic human CES1 overexpression in murine DCs generated less inflammatory phenotype and increased resistance to T cell mediated colitis. In conclusion, CES1 inhibition directs DCs differentiation towards more inflammatory phenotype, that shows stronger phagocytic capacity and supports Th17 skewing. This is associated with disrupted mitochondrial respiration and amino acids depletion.

  Study EGAS50000000230

• Cerebrospinal fluid cfDNA sequencing for classification of central nervous system glioma

  Primary central nervous system (CNS) gliomas can be classified by characteristic genetic alterations. In addition to solid tissue obtained by surgery or biopsy, cell-free DNA (cfDNA) derived from cerebrospinal fluid (CSF) is an alternative source of material for genomic analyses. Experimental design: We performed targeted next-generation sequencing (NGS) of CSF cfDNA in a representative cohort of 85 patients presenting with with suspicion of primary or recurrent glioma at two neurooncological centers. Copy-number variation (CNV) profiles, single nucleotide variants (SNVs), and small insertions/ deletions (indels) were combined into a molecular-guided tumor classification. Comparison with the solid tumor was performed for 38 cases with matching solid tissue available. Results: Cases were stratified into four groups: glioblastoma (n = 32), other glioma (n = 19), non-malignant (n = 17) and no confirmed diagnosis available (n = 17). We introduce a molecular-guided tumor classification, enabling identification of tumor entities and/ or cancer specific alterations in 75.0 % (n = 24) of glioblastoma and 52.6 % (n = 10) of other glioma cases. The overlap between CSF and matching solid tissue was highest for CNVs (26-48 %) and SNVs at pre-defined gene loci (44 %), followed by SNVs/indels identified via uninformed variant calling (8-14 %). A molecular-guided tumor classification was made 23.5 % (n = 4) of cases with no confirmed diagnosis available. Conclusions: We developed a workflow for targeted sequencing of CSF cfDNA as well as a strategy for interpretation and reporting of sequencing results based on a molecular-guided tumor classification in glioma.

  Study EGAS50000000060

• Circulating tumor DNA dynamics reveal KRAS G12C mutation heterogeneity and response to treatment with the KRAS G12C inhibitor divarasib in solid tumors

  Purpose: To inform prognosis, treatment response, disease biology, and KRAS G12C mutation heterogeneity, we conducted exploratory circulating tumor DNA (ctDNA) profiling on 134 patients with solid tumors harboring a KRAS G12C mutation treated with single-agent divarasib (GDC-6036) in a phase 1 study. Experimental design: Plasma samples were collected for serial ctDNA profiling at baseline (Cycle 1 Day 1 prior to treatment) and multiple on-treatment time points (Cycle 1 Day 15 and Cycle 3 Day 1). Results: KRAS G12C ctDNA was detectable from plasma samples in 72.9% (43/59) and 92.6% (50/54) of patients with non-small cell lung cancer and colorectal cancer, respectively, the majority of whom were eligible for study participation based on a local test detecting the KRAS G12C mutation in tumor tissue. Baseline ctDNA tumor fraction was associated with tumor type, disease burden, and metastatic sites. A decline in ctDNA level was observed as early as Cycle 1 Day 15. Serial assessment showed a decline in ctDNA tumor fraction associated with response and progression-free survival. Except for a few cases of KRAS G12C sub-clonality, on-treatment changes in KRAS G12C variant allele frequency mirrored changes in the overall ctDNA tumor fraction. Conclusion: Across tumor types, the KRAS G12C mutation likely represents a truncal mutation in the majority of patients. Rapid and deep decline in ctDNA tumor fraction was observed in patients responding to divarasib treatment. Early on-treatment dynamics of ctDNA were associated with patient outcomes and tumor response to divarasib treatment.

  Study EGAS50000000315

• Lactate metabolism in cancer stem cell fate regulation

  Tumors arise as a result of uncontrolled cell proliferation driven by mutations in the signaling pathways that regulate the balance between cell division and differentiation. Despite of that, in the tumor, cancer stem cells seem to retain their capability to differentiate, which can lead to differential chemotherapy response. We previously showed that intestinal stem cells and differentiated cells display specific metabolic features and that they interact through metabolites. Moreover, we found that altering mitochondrial metabolism is sufficient to lead and define stem cell differentiation. In cancer, cellular metabolism is highly derailed. However, how metabolic changes in the cancer cell contribute to cell fate decisions remains unknown. Recently, the crosstalk between metabolism and epigenetics has emerged as one of the players in tumor development. Here we aim to address the impact of the altered metabolism of cancer cells on their epigenome as a regulatory axis of cell identity and fate. We use human-derived tumor organoids where we introduce genetic encoded fluorescent-based reporters to visualize cell identity and metabolism at single cell level. We combine 4D live imaging with a recently developed machine learning method to re-construct cell-lineages, trace stem cell state and metabolic changes during tumor organoid development. We find that, mimicking the tumor metabolic microenvironment (Warburg effect-driven low-glucose high-lactate) leads to changes cellular metabolism and histone acetylation. Strikingly, these alterations provoke profound changes in tumor development by increasing the population of cancer stem cells and their proliferation and differentiation dynamics.

  Study EGAS50000000063

• Tumor Evolution Analysis Uncovered Immune-Escape Related Mutations in Relapsed Diffuse Large B-Cell Lymphoma

  More than 60% of patients with primary diffuse large B-cell lymphoma (DLBCL) achieve durable remissions upon the chemoimmunotherapy with R-CHOP. However, the outcome of the relapse patients is extremely poor. To characterize the molecular features that underlie treatment resistance, we performed whole-exome sequencing (n=44) and deep targeted sequencing (n=27) of 44 diagnosis-relapse DLBCL pairs. We found that the overall landscape of genetic alterations was concordant between diagnosis and relapse DLBCL. By a deep comparison of the longitudinal sequencings, we identified relapse-enriched mutations that involved in evasion of immune surveillance (CD58, TNFRSF14, PKD1, FAS), resistance to DNA damaging agents (ATM, TP53), and critical tumor suppressor gene (ZFHX3). Knockdown of ATM and ZFHX3 resulted in increased drug resistance in both germinal center B cell like (GCB) and activated B cell like (ABC) DLBCL cell lines in vitro, while the knockdown of TNFRSF14 and FAS conferred resistance specifically in GCB cell lines. Furthermore, we summarized two evolutionary patterns of DLBCL relapse as the “durable retention” model with stable clone structure, and the “selective sweep” model which exhibited significant clonal switch. Of note, genetic alterations leading to attenuated immune responses usually presented as trunk mutations in the “durable retention” model (7/12) while were newly acquired or expanded upon relapse in the “selective sweep” model (9/16). Overall, our study demonstrated that relapse DLBCL derives from distinct clonal evolutionary patterns, and highlighted immune escape as an important contributor of tumor relapse.

  Study EGAS50000000032

• Surgery in combination with immune checkpoint therapy as an effective treatment for patients with metastatic cancer.

  Abstract Surgical removal of primary tumors was shown to reverse tumor-mediated immune suppression in pre-clinical models with metastatic disease. However, how cytoreductive surgery in the metastatic setting modulates the immune responses in patients, especially in the context of immune checkpoint therapy (ICT)-containing treatments is not understood. Here, we report the first prospective, non-comparative clinical trial (N=104) using three different ICT-containing strategies plus cytoreductive or “debulking” surgery to remove the primary tumor-bearing kidney or a metastatic lesion as a treatment for patients (N=43) with metastatic clear cell renal cell carcinoma (mccRCC). For those patients (N=61) who were not candidates for cytoreductive surgery, a biopsy was obtained instead for correlative biological studies. Our data demonstrated that the combination of ICT with cytoreductive surgery was safe and feasible in patients with mccRCC. The 2-year overall survival was 84% with a median OS of 54.7 months for patients who received ICT containing regimens plus surgery. Immune-monitoring studies with co-detection by indexing (CODEX) identified distinct tumor spatial conformation of cellular subsets as a novel and improved predictor of response to ICT. Importantly, single-cell RNA-sequencing data demonstrated that surgical removal of the tumor increased antigen-presenting dendritic cell population with a concurrent reduction in KDM6B expressing immune-suppressive myeloid cells in the peripheral blood. Together, this study highlighted the feasibility of combining ICT with cytoreductive surgery in a metastatic setting and demonstrated potential enhancement of immune responses following ICT plus cytoreductive surgery in patients with metastatic disease.

  Study EGAS00001005667

• Impact of Glycolysis Inhibition on the Epigenome of Synovial Fluid T Cells in Juvenile Idiopathic Arthritis

  Juvenile Idiopathic Arthritis (JIA) encompasses a group of childhood-onset autoimmune conditions characterized by joint inflammation and mononuclear cell infiltration, especially of activated CD4+ memory/effector T (Tmem/Teff) cells in synovial fluid (SF). JIA CD4+ T cells exhibit a distinct inflammation-associated epigenomic profile, though the underlying molecular mechanisms remain unclear. This study investigates the role of glycolysis in shaping the epigenomic landscape of JIA SF CD4+ T cells, with a focus on H3K27ac histone modifications. Using an ex vivo model of antigen stimulation, CD4+ T cells were isolated from peripheral blood (PB) of healthy controls (HC) and synovial fluid from JIA patients. Cells were activated in vitro with anti-CD3/CD28 for 24 hours, and changes in H3K27ac were assessed by chromatin immunoprecipitation followed by sequencing (ChIP-seq). To explore the specific role of glycolysis, cells were also treated with 2-deoxyglucose (2DG), an inhibitor of glucose-6-phosphate production. This inhibition impacted both H3K27ac acetylation patterns and gene expression associated with the activated epigenomic regions. Additionally, HC CD4+ T cells were activated in vitro in the presence of 30% SF or 30% plasma from JIA patients. RNA sequencing (RNA-seq) was performed to examine SF-induced transcriptional changes, providing insight into the metabolic dysregulation and epigenomic alterations observed in JIA SF CD4+ T cells. This dataset includes ChIP-seq and RNA-seq data that collectively advance our understanding of how glycolysis and the SF environment contribute to the epigenetic and transcriptional regulation in JIA CD4+ T cells.

  Study EGAS50000000808

• Circulating Tumor DNA Analysis in ERBB2-Amplified Colorectal Cancer: Biomarker Analysis of the MyPathway Trial

  Combination of two HER2-directed antibodies, pertuzumab and trastuzumab (P+T), has antitumor activity in HER2-positive colorectal cancer (CRC). Although liquid biopsies are increasingly being used in clinical oncology, the association between tumor and circulating tumor DNA (ctDNA) ERBB2 status and ctDNA monitoring for early response and resistance are unknown. Eighty-five patients with ERBB2-amplified and/or overexpressed CRC were treated with P+T in the MyPathway trial; 42 had ctDNA testing at Cycle 1 Day 1 (C1D1), and 38 had longitudinal plasma tested for ctDNA. We analyzed the ctDNA vs tissue ERBB2 concordance, genomic co-alterations, and ctDNA dynamics and association with response. Forty-one (98%) of 42 patients had genomic alterations detected in ctDNA at C1D1, and 29 (69%) had ERBB2 amplification in ctDNA. There was a strong correlation between ERBB2 copy number on next-generation sequencing in tissue and C1D1 ERBB2 ctDNA copy number. Thirty-seven percent achieved a molecular response by C3D1 on P+T, which was associated with prolonged progression-free survival and overall survival (OS). CDKN2A and KRAS mutations were associated with shorter OS, and a trend was seen with PIK3CA mutations. Several emerging co-alterations were identified in ctDNA at progression, including in the MAPK and PI3K pathways and other tyrosine receptor kinases. ctDNA can detect ERBB2 amplification in many, but not all, patients with ERBB2 amplification detected in tumor samples. ctDNA molecular response was associated with better survival, and ctDNA co-alterations may offer insights into mechanisms of intrinsic and acquired resistance.

  Study EGAS50000000916

• Nanopore sequencing of blood, saliva, buccal mucosa samples of patients with inherited retinal disorders

  X-linked genetic disorders typically affect females less severely than males due to the presence of a second X chromosome not carrying the deleterious variant. However, the phenotypic expression in females is highly variable, which may be explained by an allelic skew in X chromosome inactivation. Accurate measurement of X inactivation skew is crucial to understand and predict disease phenotype in carrier females, with prediction especially relevant for degenerative conditions. We propose a novel approach using nanopore sequencing to quantify skewed X inactivation accurately. By phasing sequence variants and methylation patterns, this single assay reveals the disease variant, X inactivation skew, its directionality, and is applicable to all patients and X-linked variants. Enrichment of X-chromosome reads through adaptive sampling enhances cost-efficiency. Our study includes a cohort of 16 X-linked variant carrier females affected by two X-linked inherited retinal diseases: choroideremia and RPGR-associated retinitis pigmentosa. As retinal DNA cannot be readily obtained, we instead determine the skew from peripheral samples (blood, saliva and buccal mucosa), and correlate it to phenotypic outcomes. This revealed a strong correlation between X inactivation skew and disease presentation, confirming the value in performing this assay and its potential as a way to prioritise patients for early intervention, such as gene therapy currently in clinical trials for these conditions. Our method of assessing skewed X inactivation is applicable to all long-read genomic datasets, providing insights into disease risk and severity and aiding in the development of individualised strategies for X-linked variant carrier females.

  Study EGAS50000000440

• PCR-free shallow whole genome sequencing for chromosomal copy number detection from plasma of cancer patients is an efficient alternative to the conventional PCR-based approach

  Somatic copy number alterations (SCNAs) can be detected in tumor cell-free DNA (cfDNA) by shallow whole genome sequencing (sWGS). Polymerase chain reaction (PCR) is typically included in library preparations, but a PCR-free method could be an effective alternative for research and diagnostics. Peripheral blood samples (n=22) were collected in EDTA-containing tubes from non-small cell lung cancer patients (n=10) and healthy donors (n=12). PCR and PCR-free library preparations underwent sWGS and the following metrics were compared: mapping quality, proportion of duplicate reads, genome coverage, sequencing read length and copy number profiles. The PCR-free method was further evaluated on archived plasma samples (n=33) of varying pre-analytical quality collected in EDTA or lithium heparin tubes. Archived samples were from acute myeloid leukaemia (n=15) and B-cell lymphoma (n=18) patients. The percentage of unique reads was significantly higher for the 22 samples analysed by PCR-free (95.7%) compared to PCR (84%). PCR-free libraries from 33 archived samples produced data of a comparable quality, with 95.4% unique reads. All other evaluated metrics were highly comparable for PCR and PCR-free library preparations. We conclude that SCNA detection can be performed by PCR-free sWGS on cfDNA from cancer patients and is feasible for diagnostic and archival plasma collected in EDTA- or lithium heparin-containing tubes. Our results pave the way for an automated cfDNA workflow for cancer plasma samples.

  Study EGAS00001004692

• Coding and regulatory somatic profiling of triple-negative breast cancer in Sub-Saharan African patients

  The burden of triple-negative breast cancer (TNBC) may be shaped by genetic factors, particularly inherited and somatic mutation profiles. However, data on this topic remain limited, especially for the African continent, where a higher TNBC incidence is observed. In the age of precision medicine, cataloguing TNBC diversity in African patients becomes imperative. We performed whole exome sequencing, including untranslated regions, on 30 samples from Angola and Cape Verde, which allowed to ascertain on potential regulatory mutations in TNBC for the first time. A high somatic burden was observed for the African cohort, with 86% of variants being so far unreported. Recurring to predictive functional algorithms, 17% of the somatic single nucleotide variants were predicted to be deleterious at the protein level, and 20% overlapped with candidate cis-regulatory elements controlling gene expression. Several of these somatic functionally-impactful mutations and copy number variation (mainly in 1q, 8q, 6 and 10p) occur in known BC- and all cancer-driver genes, enriched for several cancer mechanisms, including response to radiation and related DNA repair mechanisms. TP53 is the top of these known BC-driver genes, but our results identified possible novel TNBC driver genes that may play a main role in the African context, as TTN, CEACAM7, DEFB132, COPZ2 and GAS1. These findings emphasize the need to expand cancer omics screenings across the African continent, the region of the globe with highest genomic diversity, accelerating the discovery of new somatic mutations and cancer-related pathways.

  Study EGAS50000001013

• Whole-genome and transcriptome versus panel sequencing in precision oncology: A translational-clinical comparison

  Precision oncology offers new cancer treatment options. While the use of whole-genome (WGS) and transcriptome (TS) sequencing remains often confined to specific clinical programs, molecular analysis with gene panels is more common clinical practice. Methods: Twenty patients with rare or advanced tumors sequenced by WGS/TS within the DKFZ/NCT/DKTK-MASTER program from 2015 to 2020 (MASTER 1) were additionally sequenced with the TSO500 DNA and TST170 RNA Assay by Illumina. Original Molecular Tumor Board (MTB) therapy recommendations (TR) were juxtaposed with TR from NGS panel data analysis (comparison 1). To mitigate temporal differences between the two sequencings, WGS/TS was bioinformatically reanalyzed (MASTER 2), and new TR were again compared (comparison 2). Results: In comparison 1, 47.1% of the MASTER 1’s TR and 62.8% of the panel’s TR were identical, as were 45.9% of the TR of MASTER 2 and 54.9% of the panel in comparison 2. In both MASTER 1 and MASTER 2, 26.5% and 36.1% of TR, respectively, relied on Biomarkers (BM) which were not covered by the panel. Eight of the ten MTB-derived therapy implementations were also supported in the panel, the two remaining ones being based on BM absent from the panel. Conclusions: About half of the TR in WGS/TS and panel sequencing were identical. Around a third of TR issued by WGS/TS were based on BM which are not covered by the panel, leading to therapy implementations. The presented data underscores the clinical value of comprehensive molecular analyses in individuals with advanced cancers.

  Study EGAS50000000431

• Mutational signatures in normal tissues of patients treated with clinical stage G-quadruplex binder CX5461

  Drug induced photoreactions following exposure to UVA have been documented for several therapeutic compound classes, including psoralens and fluoroquinolones. CX5461, a quinolone-derived small molecule currently in phase 1 studies for cancer treatment, has been shown to be a clinically significant UV photosensitizer. Here we examined the mutagenic spectrum of CX5461 in normal tissues of patients undergoing treatment, by comparing deep whole genome sequencing (WGS) from skin biopsies taken from non-sun exposed areas pre and post drug exposure of patients in the Canadian Cancer Clinical Trials Group (CCTG) IND.231 trial (NCT02719977). We show that a low burden of SNVs (<100) and indels (<40) mutations was detected in post CX5461 skin biopsies, with a pattern resembling background or sequencing artefact signatures. To further characterize the mutational pattern of CX5461 photoreactivity, we compared mutational profiles of human retinal pigment epithelial cells exposed to UVA, CX5461, or a combination of both UV and CX5461. Treatment with CX5461 or UVA alone, resulted in a low SNV burden. In contrast, human cells exposed to UVA and CX5461 had a two to four orders higher mutational burden, with a distinct mutational pattern characterized by T>A, T>C mutations, dissimilar to that of UV alone or CX5461 alone. Our results suggest that human exposure to CX5461 alone is not associated with high levels of single base mutagenesis. However, CX5461 is a very potent UV photosensitizer, underscoring the need for adequate UV screening and controls when studying the effects of photosensitizing small molecules on single base DNA variation.

  Study EGAS50000001145

• Transcriptome profiling of megakaryocytes and platelets: application to GP9- and IKZF5-related thrombocytopenia

  Platelets are anucleate cells produced in the bone marrow and derived from large progenitor cells called megakaryocytes (MKs). Platelets receive RNA transcripts from their progenitorial MKs during thrombopoiesis. However, the correspondence between platelet and MK transcriptomes is poorly understood, particularly in the context of germline mutations that cause platelet formation defects or thrombocytopenia. We have studied the effects of two such mutations on MK and platelet transcriptomes. We generated immortalized MK cell line (imMKCL)-based models of Bernard-Soulier syndrome and IKZF5-related thrombocytopenia. MKs derived from imMKCLs with either a homozygous deletion of GP9 (GP9-/-) or a heterozygous Y121F variant in IKZF5 (IKZF5WT/Y121F) exhibited reduced proplatelet formation (reductions of 96% and 57%, respectively). Platelets from patients with either GP9-/- or IKZF5WT/Y121F genotypes had broad transcriptomic dysregulation, suggesting that (pro)platelet formation defects due to mutations in glycoprotein receptor and transcription factor genes such as GP9 and IKZF5 affect the trafficking of transcripts from MKs into platelets. RNA-seq data from MKs at four stages of differentiation revealed widespread but distinct changes in expression over time between the GP9-/- and the IKZF5WT/Y121F genotypes. Dysregulated genes in GP9-/- MKs were enriched for RNA metabolism and actin/tubulin folding pathways while those in IKZF5WT/Y121F MKs were enriched for cell cycle pathways. Most of these genes were also dysregulated in the platelets of patients with the corresponding diseases. Our results suggest that patients with inherited forms of thrombocytopenia present with a transcript trafficking defect during platelet formation.

  Study EGAS50000001276

• Complex genotype-phenotype relationships shape the response to treatment of Down Syndrome Childhood Acute Lymphoblastic Leukaemia

  Extensive genetic and epigenetic variegation has been demonstrated in many malignancies. Importantly, their interplay has the potential to contribute to disease progression and treatment resistance. To shed light on the complex relationships between these different sources of intra-tumour heterogeneity, we explored their relative contributions to the evolutionary dynamics of Acute Lymphoblastic Leukaemia (ALL) in children with Down syndrome, which has particularly poor prognosis. We quantified the tumour propagating potential of genetically distinct sub-clones using serial transplantation assays and SNP-arrays. While most leukaemias were characterized by a single dominant subclone, others were highly heterogeneous. Importantly, we provide clear and direct evidence that genotypes and phenotypes with functional relevance to leukemic progression and treatment resistance can co-segregate within the disease. Hence, individual genetic lesions can be restricted to well-defined cell immunophenotypes, corresponding to different stages of the leukemic differentiation hierarchy and varied proliferation potentials. As a result of this difference in fitness, which can be accurately quantified via competitive transplantation assays, matching diagnostic, post-treatment, and relapse leukaemias can be dominated by different genotypes, including pre-leukemic clones persisting throughout the disease progression and treatment. Intriguingly, plasticity also appears to be a temporally defined property that can segregate with genotype. These results suggest that Down Syndrome ALL should be viewed as a complex matrix of cells exhibiting genetic and epigenetic heterogeneity that foster extensive clonal evolution and competition. Therapeutic intervention reshapes this ‘eco-system’ and may provide the right conditions for the preferential expansion of selected compartments and subsequently relapse.

  Study EGAS50000001287

• Heterogeneity of IKZF1 genomic alterations and risk of relapse in childhood B-cell precursor acute lymphoblastic leukemia

  Purpose Genomic alterations of IKZF1 are common and associated with adverse clinical features in B-ALL, but the relationship between the type of IKZF1 alteration, disease subtype and outcome are incompletely understood. Patients and Methods We used transcriptome and genomic sequencing, and SNP microarray data to determine leukemia subtype and genomic alterations for 688 pediatric patients with B-ALL enrolled into St. Jude Total Therapy 15 and 16 studies. Results IKZF1 alterations were identified in 115 (16.7%) patients, commonly in BCR::ABL1 (78%) and CRLF2-rearranged, BCR::ABL1-like B-ALL (70%), and were associated with five-year cumulative incidence of relapse (CIR) of 14.8 ± 3.3% compared to 5.0 ± 0.9% for patients without any IKZF1 alteration (P < 0.0001). IKZF1 deletions of exon 4-7 (P = 0.0002), genomic IKZF1plus with any IKZF1 deletion (P = 0.006) or with focal IKZF1 deletion (P = 0.0007), and unfavorable genomic subtypes (P < 0.005) were independently adversely prognostic. Patients with both IKZF1 exon 4-7 deletion and unfavorable genomic subtype had a significantly increased risk of relapse (HR = 58.3; 95% CI, 11.9 – 285.4, P < 0.0001) whereas IKZF1 sequence mutations were not independently predictive of outcome. Conclusions Genomic IKZF1plus with any IKZF1 deletion, IKZF1 deletion of exon 4-7, and unfavorable subtype confer increased risk of relapse. The combination of unfavorable genomic subtype and IKZF1 deletion of exon 4-7 identified patients at greatest risk of relapse despite MRD-directed therapy. The type of IKZF1 alteration together with subtype are informative for risk stratification and predict response in patients with B-ALL.

  Study EGAS50000000106

• Case Report: Pre-Clinical Combination Targeting VEGF and PI3K in a Rare, Aggressive Mixed Endometrial Carcinoma

  We report a rare case of a young patient (VENUS 167) initially diagnosed with grade 1 endometrioid endometrial cancer, which, following endocrine treatment, presented with mixed aggressive carcinoma with three distinct histologic patterns: grade 1 endometrioid, large cell neuroendocrine, and undifferentiated carcinoma. The surgical specimen at the time of disease progression was used to establish OPTO.85, a patient-derived organoid (PDO), followed by a corresponding organoid-derived xenograft (ODX). Multi-omic analyses confirmed that OPTO.85 accurately reflected the patient’s tumor characteristics. WES analysis identified oncogenic alterations in PIK3CA, ARID1A, and CTNNB1. Further RNA-seq and ATAC-seq analyses revealed enrichment in VEGF and Wnt signaling pathways, suggesting therapeutic vulnerabilities. A high-throughput drug screen was conducted using ApexBio-Approved and Epigenetic Drug Libraries, along with Kinase Inhibitor and Tool Compound Libraries developed at the Ontario Institute of Cancer Research. The OPTO.85 PDO exhibited sensitivity to PI3K inhibitors and responsiveness to VEGF inhibition. Cediranib demonstrated synergy with BKM120, significantly reducing organoid growth. This combination also showed in vivo efficacy in the ODX model, where dual inhibitors significantly suppressed tumor growth compared to single compounds. This case exemplifies the impact of genomic profiling and patient-derived models in identifying actionable molecular changes in rare cancers with limited therapeutic options and poor prognosis. It highlights that high-throughput sequencing for individual patient tumors and generation of patient-derived models are feasible in endometrial cancer. This preclinical model may assist clinical decision and personalized therapy requiring validation in prospective studies.

  Study EGAS50000001665

• An autoinflammatory RIG-I variant causing Singleton-Merten Syndrome associates with small non-coding Y-RNAs

  The RNA sensor RIG-I performs a critical role in surveying the cytoplasm for the presence of viral nucleic acids and initiating the downstream anti-viral type I interferon pathway. Through recognition of specific features, such as the presence of a 5’ di/triphosphate motif or highly structured base-paired regions, RIG-I effectively differentiates between RNA of viral- and self-origin. In Singleton-Merten syndrome (SMS), gain-of-function variants in RIG-I lead to a breakdown in this surveillance system and results in aberrant sensing of self-RNAs and deleterious upregulation of type I interferons. The identity of the self-RNAs binding to gain-of-function RIG-I mutants has remained elusive and would provide a greater understanding of the aetiology of SMS. Here we used an infrared individual-nucleotide resolution UV-crosslinking and immunoprecipitation (irCLIP) approach to determine the RNA profile bound to a previously characterised ATPase-deficient SMS variant, RIG-IC268F. irCLIP identified a broad array of self-RNAs, primarily those transcribed by RNA polymerase III, that were bound to RIG-I. Subsequent native RNA immunoprecipitation confirmed a prominent and specific interaction between RIG-IC268F and Y-RNAs, a family of four structurally similar non-coding RNAs. Manipulation of Y-RNAs alone by targeting either Y-RNA transcripts or Y-RNA stabilising proteins was insufficient to negate RIG-I induce interferon, hinting at a broader profile of RNA polymerase III-derived RNAs being influential in driving the sterile activation of gain-of-function RIG-I variants.

  Study EGAS50000001660

• Sequencing and analysis of a South Asian-Indian personal genome

  India with over 1.3 billion people is estimated to carry three times more genetic diversity compared to Europe. Next-generation sequencing technologies have facilitated the understanding of diversity by enabling whole genome sequencing at greater speed and lower cost. While multiple genomes from people of European and a few of Asian descent have been sequenced, only recently a single male genome from the Indian subcontinent at sufficient depth and coverage was published. We have in this study sequenced and analyzed the genome of a South Asian Indian female (SAIF) from the Indian state of Kerala. We identified over 3.4 million SNPs in this genome including over 89,873 private variations. Comparison of the SAIF genome with several published personal genomes revealed that this individual shared ~50% of the SNPs with each of these genomes. Analysis of the SAIF mitochondrial genome showed that it is closely related to the U1a3 haplogroup which has been previously observed in Kerala. We assessed the SAIF genome for SNPs with health and disease consequences and found that the individual was at a higher risk for multiple sclerosis and a few other diseases. In assessing SNPs that modulate drug response we found a variation that predicts a favorable response to metformin used for treating diabetes. SNPs predictive of adverse reaction to warfarin indicated that the SAIF individual is not at risk for bleeding if treated with typical doses of warfarin. In addition to these, we report the presence of several other SNPs of medical relevance.

  Study EGAS00001000328

• Detection of Clinically Relevant Genetic Variantsin Autism Spectrum Disorderby Whole-Genome Sequencing

  Autism Spectrum Disorder (ASD) demonstrates high heritability and familial clustering, yet the genetic causes remain only partially understood as a result of extensive clinical and genomic heterogeneity. Whole-genome sequencing (WGS) shows promise as a tool for identifying ASD risk genes as well as unreported mutations in known loci, but an assessment of its full utility in an ASD group has not been performed. We used WGS to examine 32 families with ASD to detect de novo or rare inherited genetic variants predicted to be deleterious (loss-of-function and damaging missense mutations). Among ASD probands, we identified deleterious de novo muta- tions in six of 32 (19%) families and X-linked or autosomal inherited alterations in ten of 32 (31%) families (some had combinations of mutations). The proportion of families identified with such putative mutations was larger than has been previously reported; this yield was in part due to the comprehensive and uniform coverage afforded by WGS. Deleterious variants were found in four unrecognized, nine known, and eight candidate ASD risk genes. Examples include CAPRIN1 and AFF2 (both linked to FMR1, which is involved in fragile X syndrome), VIP (involved in social-cognitive deficits), and other genes such as SCN2A and KCNQ2 (linked to epilepsy), NRXN1, and CHD7, which causes ASD-associated CHARGE syndrome. Taken together, these results suggest that WGS and thorough bioinformatic analyses for de novo and rare inherited mutations will improve the detection of genetic variants likely to be associated with ASD or its accompanying clinical symptoms.

  Study EGAS00001000850

• Sequencing data from triple-negative breast cancer tumors

  Autoimmunity and anti-cancer immunity lie on the same biological continuum1,2, but their link remains obscure. The paraneoplastic neurological syndrome anti-NMDA receptor (NMDAR) encephalitis (ANRE) is a paradigm for their connectivity3 given that intratumoral NMDAR expression correlates with the generation of anti-NMDAR antibodies4,5. Here, we verify ectopic expression of GluN1 and GluN2B NMDAR sub-units in triple-negative breast cancer (TNBC)6 and model this using orthotopic TNBC tumors with inducible expression of GluN1-GluN2B NMDARs. We show that NMDAR expression is sufficient to induce B cell recruitment and their affinity maturation, consistent with an integrated adaptive immune response. Reconstruction of extended intratumoral B cell phylogenies and cryo-EM structural analyses demonstrated that affinity-matured hypermutated and class-switched antibodies emerged from pre-existing germline-configuration lower-affinity anti-NMDAR antibodies. Distinct matured antibodies targeted specific epitopes and induced conformational rearrangements within the NMDAR amino-terminal domain, predictive of their functional effects, ranging from inhibition to potentiation. Passive transfer of an NMDAR-potentiating antibody caused autonomic dysregulation and lowered the seizure threshold in healthy female mice, recapitulating key diagnostic criteria of ANRE4,5. We further identify a correlation between intratumoral NMDAR expression and anti-NMDAR antibody titers in TNBC patients. Taken together, our data establish a direct connection between intratumoral NMDAR expression, antibody maturation, and the onset of autoimmunity. These findings suggest that germline-encoded anti-NMDAR antibodies contribute to immune surveillance but can also trigger autoimmune disease upon maturation, revealing a mechanistic tradeoff between cancer immunity and neurotoxicity.

  Study EGAS50000001585

• Sequencing data from triple-negative breast cancer tumors

  Autoimmunity and anti-cancer immunity lie on the same biological continuum1,2, but their link remains obscure. The paraneoplastic neurological syndrome anti-NMDA receptor (NMDAR) encephalitis (ANRE) is a paradigm for their connectivity3 given that intratumoral NMDAR expression correlates with the generation of anti-NMDAR antibodies4,5. Here, we verify ectopic expression of GluN1 and GluN2B NMDAR sub-units in triple-negative breast cancer (TNBC)6 and model this using orthotopic TNBC tumors with inducible expression of GluN1-GluN2B NMDARs. We show that NMDAR expression is sufficient to induce B cell recruitment and their affinity maturation, consistent with an integrated adaptive immune response. Reconstruction of extended intratumoral B cell phylogenies and cryo-EM structural analyses demonstrated that affinity-matured hypermutated and class-switched antibodies emerged from pre-existing germline-configuration lower-affinity anti-NMDAR antibodies. Distinct matured antibodies targeted specific epitopes and induced conformational rearrangements within the NMDAR amino-terminal domain, predictive of their functional effects, ranging from inhibition to potentiation. Passive transfer of an NMDAR-potentiating antibody caused autonomic dysregulation and lowered the seizure threshold in healthy female mice, recapitulating key diagnostic criteria of ANRE4,5. We further identify a correlation between intratumoral NMDAR expression and anti-NMDAR antibody titers in TNBC patients. Taken together, our data establish a direct connection between intratumoral NMDAR expression, antibody maturation, and the onset of autoimmunity. These findings suggest that germline-encoded anti-NMDAR antibodies contribute to immune surveillance but can also trigger autoimmune disease upon maturation, revealing a mechanistic tradeoff between cancer immunity and neurotoxicity.

  Study EGAS50000001598

• Capturing disease severity in LIS1-lissencephaly reveals proteostasis dysregulation in patient-derived forebrain organoids

  Lissencephaly is a developmental cortical malformation characterized by reduced to absent gyri and a disorganized cortex, often leading to severe impairments in affected individuals and a reduced life expectancy. Heterozygous mutations in the LIS1 gene, encoding a regulator of the microtubule motor dynein, cause lissencephaly with different clinical severities. While the clinical disease spectrum correlates with the degree of lissencephaly, location and type of mutation may not. We leveraged forebrain-type organoids from LIS1-lissencephaly patients, diagnosed with mild, moderate or severe lissencephaly to investigate, in a cytoarchitecture and multi-omics approach, how the severity degree in patients might relate to specific mutations in the LIS1 gene. We questioned which processes during cortical development might be differentially affected by severity grade, and whether they could be pharmacologically targeted. We found alterations in neurodevelopment often with a severity-dependent gradient. Specifically, we identified alterations of the cytoarchitecture, progenitor cell homeostasis and neurogenesis. Particularly important disease-linked molecular mechanisms were microtubule destabilization, WNT-signaling, and perturbed cadherin- and unfolded protein-binding. Some mechanisms exhibited a severity-dependent gradient, or were specific to a severe grade. We present strategies to reverse phenotypic changes in LIS1-patient organoids, and an in silico approach with therapeutic potential. Thus, we show that different LIS1-severity grades can be recapitulated in vitro, that there is a direct link between the phenotype and genotype, that organoid-based disease modeling can identify molecular underpinnings of malformations of cortical development and that organoids provide a valid platform to develop and test therapeutic strategies.

  Study EGAS00001008227