23014 results for "ega"

in 74.80 milliseconds.

• DNA methylation and the adverse metabolic outcomes of adiposity

  Overweight and obesity affect ~1.5 billion people worldwide, and are major risk factors for type-2 diabetes (T2D), cardiovascular disease and related metabolic and inflammatory disturbances.1,2 Although the mechanisms linking adiposity to its clinical sequelae are poorly understood, recent studies suggest that adiposity may influence DNA methylation,3-6 a key regulator of gene expression and molecular phenotype.7 Here we use epigenome-wide association to show that body mass index (BMI, a key measure of adiposity) is associated with widespread changes in DNA methylation (187 genetic loci at P<1x10-7, range P=9.2x10-8 to 6.0x10-46; N=10,261 samples). Genetic association analyses demonstrate that the alterations in DNA methylation are predominantly the consequence of adiposity, rather than the cause. We find the methylation loci are enriched for functional genomic features in multiple tissues (P<0.05), and show that sentinel methylation markers identify gene expression signatures at 38 loci (P<9.0x10-6, range P=5.5x10-6 to 6.1x10-35, N=1,785 samples). The methylation loci identified highlight genes involved in lipid and lipoprotein metabolism, substrate transport, and inflammatory pathways. Finally, we show that the disturbances in DNA methylation predict future type-2 diabetes (relative risk per 1SD increase in Methylation Risk Score: 2.3 [2.07-2.56]; P=1.1x10-54). Our results provide new insights into the biologic pathways influenced by adiposity, and may enable development of new strategies for prediction and prevention of type-2 diabetes and other adverse clinical consequences of obesity.

  Study EGAS00001001922

• Immunodeficiency syndrome caused by LCP1 mutations

  Severe congenital neutropenia (SCN) is characterized by neutropenia, recurrent infections and an increased leukemia risk. Multiple genetic defects underlying SCN have been identified, but a genetic diagnosis is still lacking in a significant proportion of patients. Here, we report three independent pedigrees with autosomal dominant variants in LCP1. Variants c.740-1G>T and c.740-20_744del resulted in the same alternatively spliced RNA product, causing an in-frame deletion (p.A247_E254del). Variant c.509C>T in the third pedigree resulted in p.S170L. Variant carriers suffer from neutropenia, poor or complete lack of response to G-CSF treatment, and variable degrees of lymphopenia. Patients with A247_E254del presented with tetraploid cells in the bone marrow, indicative of disturbed cytokinesis. In one of these kindreds, two individuals developed acute leukemia. G-CSF non-responsiveness and defective cell cycling were repaired upon correction of the LCP1 A247_E254del variant in patient-derived induced pluripotent stem cells, supporting the monogenic origin of the disease. Indicative of their gain-of-function, both A247_E254del and S170L variants increased F-actin bundling and formation of abnormal protrusions. Single-cell transcriptome analysis of A247_E254del bone marrow-derived hematopoietic stem and progenitor cells (HSPCs) showed deregulation of signaling pathways controlling mitosis in multi-lineage and lymphoid-primed HSPC subsets. We conclude that activating LCP1 variants cause a new autosomal dominant hematopoietic disorder. Depending on the consequences of the LCP1 variants for protein structure, patients may suffer from G-CSF refractory severe neutropenia, lymphopenia, defective cytokinesis and predisposition to leukemia.

  Study EGAS00001008293

• The genomic landscape of pediatric acute lymphoblastic leukemia

  Acute lymphoblastic leukemia (ALL) is the most common pediatric cancer and arises from B- or T-lineage lymphocyte precursors. ALL comprises dozens of subtypes, and genomic analyses have largely been performed within these subtypes. Here we analyze pediatric ALL genomes across all subtypes, including 768 whole genomes, 1,729 exomes, and 1,889 transcriptomes in 2,754 patients. Most ALL subtypes harbor 4 or more driver alterations per sample, similar to adult cancers, despite low mutation burdens. Hyperdiploid B-ALL copy gains are likely acquired early and synchronously, with copy gains occurring before ultraviolet-induced mutations. By contrast, ultraviolet-induced mutations precede copy gains in iAMP21 B-ALL. Overall, we identified 378 putative ALL driver genes. Most driver alterations vary in prevalence across ALL subtypes, with B-ALL enriched for Ras and B-lineage-related alterations, and T-ALL enriched for PI3K, JAK, and cell cycle alterations. Most B-ALL (54.3%) and T-ALL (51.2%) samples bear at least one rare driver gene alteration (present in less than 2% of samples), including 30 putative novel cancer driver genes associated with ubiquitination, SUMOylation, non-coding, and other functions. Known or novel alterations associated with poor outcomes in specific subtypes include TBL1XR1 alterations in ETV6-RUNX1 patients, CREBBP in Ph-like-CRLF2, SETD2 in hyperdiploid, and PTEN in TAL1 patients. Intriguingly, DUX4 and KMT2A subtypes separate into CEBPA/FLT3- or NFATC4-expressing subgroups with potential clinical implications. Together, these results deepen understanding of ALL etiology and outcomes, and facilitate ALL model development.

  Study EGAS00001005250

• Methylation differences in trisomy 21 using monozygotic twins

  We have hypothesized that methylation differences induced by trisomy 21 (T21) contribute to the phenotypic characteristics and heterogeneity in T21. In order to determine the methylation differences in T21 without the interference of the interindividual genomic variation, we have used fetal skin fibroblasts from monozygotic (MZ) twins discordant for T21. We also used skin fibroblasts from MZ twins concordant for T21, normal MZ twins without T21, and unrelated normal and T21 individuals. We applied Reduced Representation Bisulfite Sequencing (RRBS) to generate genome wide nucleotide resolution of DNA methylation based on high throughput sequencing between each pair of samples. RRBS revealed differentially methylated promoter regions (DMRs) in MZ twins discordant for T21 that have also been observed in unrelated normal and T21 individuals. The identified DMRs are enriched for genes involved in embryonic organ morphogenesis. These DMRs are maintained in iPS cells generated from this twin pair and are correlated with the gene expression changes. We have also observed an increase in DNA methylation level in the T21 methylome compared to the normal euploid methylome. This observation is concordant with the up regulation of DNA methyltransferase enzymes (DNMT3B and DNMT3L) and down regulation of DNA demethylation enzymes (TET2 and TET3) in the iPSC of the T21 versus normal twin. Additionally, two sets of T21 MZ twins discordant for heart defect highlighted DMRs in genes that are associated with heart development. In conclusion the study of DNA methylation differences in MZ twins discordant for genomic abnormalities is a promising approach to understand the molecular pathophysiology of aneuploidies.

  Study EGAS00001001051

• A targeted gene panel that covers coding, noncoding, and short tandem repeat regions improves the diagnosis of patients with neurodegenerative diseases

  Genetic testing for diagnosis of neurodegenerative diseases (NDs) is highly challenging because of genotype heterogeneity and overlapping clinical manifestations. Targeted-gene panels (TGPs), coupled with next-generation sequencing (NGS), can facilitate the profiling of a large repertoire of ND-related genes and help with differential diagnosis. Due to the technical limitations inherent in NGS and TGPs, short tandem repeat (STR) variations are often ignored. However, STR expansions are known to cause such NDs as Huntington’s disease and several types of spinocerebellar ataxias, such as SCA type 3 (SCA3).Here, we studied the clinical utility of a custom-made TGP that targets 199 NDs and 311 ND-associated genes on 118 undiagnosed patients. At least one known pathogenic or likely pathogenic variation was found in 54 (45.8%) patients; 27 (22.9%) of the 118 patients demonstrated clinical profiles that matched the variants. Based on the findings from the TGP, we refined the original diagnosis of 16 patients. A high concordance of single-nucleotide polymorphisms (99.9%) and small insertions and deletions (93.2%) were observed when comparing the results from TGP and whole-exome sequencing of four patients. We tested an in-house STR detection algorithm, and it reached a specificity of 0.88 and a sensitivity of 0.82 in our SCA3 cohort. This study also uncovered a trove of novel and recurrent variants that may enrich the repertoire of ND-related genetic markers. We propose that a combination of comprehensive TGPs and the bioinformatics analysis pipeline can improve the clinical diagnosis and accelerate the diagnosis of NDs.

  Study EGAS00001003737

• Organoid cultures of early-onset colorectal cancers reveal distinct and rare genetic profiles

  Sporadic early-onset colorectal cancer (EOCRC) has bad prognosis, yet is poorly represented by cell line models. We examine the key mutational and transcriptomic alterations in an organoid biobank enriched in EOCRCs. We established paired cancer (n=32) and normal organoids (n=18) from 20 patients enriched in microsatellite-stable EOCRC. Exome and transcriptome analysis was performed. We observed a striking diversity of molecular phenotypes, including PTPRK-RSPO3 fusions. Transcriptionally, RSPO fusion organoids resembled normal colon organoids and were distinct from APC mutant organoids, with high BMP2 and low PTK7 expression. Single cell transcriptome analysis confirmed the similarity between RSPO fusion organoids and normal organoids, with a propensity for maturation on Wnt withdrawal, whereas the APC mutant organoids were locked in progenitor stages. CRISPR/Cas9 engineered mutation of APC in normal human colon organoids led to upregulation of PTK7 protein and suppression of BMP2, but less so with an engineered RNF43 mutation. The frequent co-occurrence of RSPO fusions with SMAD4 or BMPR1A mutation was confirmed in TCGA database searches. RNF43 mutation was found in organoid from a leukaemia survivor with a novel mutational signature; and organoids with POLE proofreading mutation displayed ultramutation. The cancer organoid genomes were stable over long culture periods, while normal human colon organoids tended to be subject to clonal dominance over time. These organoid models enriched in EOCRCs with linked genomic data fill a gap in existing CRC models and reveal distinct genetic profiles and novel pathway cooperativity.

  Study EGAS00001004063

• Single cell RNA-seq mapping of nasal and tracheobronchial airways in human healthy volunteers

  The respiratory tract constitutes an elaborated line of defense based on a unique cellular ecosystem. Single-cell profiling methods enable the investigation of cell population distributions and transcriptional changes along the airways. We have explored cellular heterogeneity of the human airway epithelium in 10 healthy living volunteers by single-cell RNA profiling. 77,969 cells were collected by bronchoscopy at 35 distinct locations, from the nose to the 12th division of the airway tree. The resulting atlas is composed of a high percentage of epithelial cells (89.1%), but also immune (6.2%) and stromal (4.7%) cells with peculiar cellular proportions in different sites of the airways. It reveals differential gene expression between identical cell types (suprabasal, secretory, and multiciliated cells) from the nose (MUC4, PI3, SIX3) and tracheobronchial (SCGB1A1, TFF3) airways. By contrast, cell-type specific gene expression was stable across all tracheobronchial samples. Our atlas improves the description of ionocytes, pulmonary neuro-endocrine (PNEC) and brush cells, which are likely derived from a common population of precursor cells. We also report a population of KRT13 positive cells with a high percentage of dividing cells which are reminiscent of “hillock” cells previously described in mouse. Robust characterization of this unprecedented large single-cell cohort establishes an important resource for future investigations. The precise description of the continuum existing from nasal epithelium to successive divisions of lung airways and the stable gene expression profile of these regions better defines conditions under which relevant tracheobronchial proxies of human respiratory diseases can be developed.

  Study EGAS00001004082

• Targeted and shallow whole genome sequencing identifies therapeutic opportunities in p53abn endometrial cancers

  Purpose: Shallow whole genome sequencing (sWGS) can detect copy number (CN) aberrations. In high-grade serous ovarian (HGSOC) sWGS identified CN signatures such as homologous recombination deficiency (HRD) to direct therapy. We applied sWGS with targeted sequencing to p53abn endometrial cancers (ECs) to identify additional prognostic stratification and therapeutic opportunities. Experimental Design: sWGS and targeted panel sequencing was performed on formalin-fixed paraffin-embedded p53abn ECs. CN alterations, mutational data and CN signatures were derived, and associations to clinicopathologic and outcomes data were assessed. Results: In 187 p53abn ECs, 5 distinct CN signatures were identified. Signature 5 was associated with BRCA1/2 CN loss with features similar to HGSOC HRD signature. 22% potential HRD cases were identified, 35 patients with signature 5, and 8 patients with BRCA1/2 somatic mutations. Signatures 3 and 4 were associated with whole genome duplication, and CCNE1, ERBB2 and MYC amplifications, with mutations in PIK3CA enriched in signature 3. We observed improved overall survival (OS) for patients with signature 2 and worse OS for signatures 1 and 3. 28% of patients had CCNE1 amplification and this subset was enriched with carcinosarcoma histotype. 34% of patients, across all histotypes, had ERBB2 amplification and/or HER2 overexpression on immunohistochemistry, which was associated with worse outcomes. Mutations in PPP2R1A (29%) and FBXW7 (16%) were among the top 5 most common mutations. Conclusions: sWGS and targeted sequencing identified therapeutic opportunities in 75% of p53abn EC patients. Further research is needed to determine the efficacy of treatments targeting these identified pathways within p53abn ECs.

  Study EGAS00001007661

• Pediatric Papillary Thyroid Carcinoma Whole Exome Sequencing

  BACKGROUND: Pediatric papillary thyroid carcinoma (PPTC) is distinct from adult-onset disease. Despite higher rates of regional and distal metastases and more frequent recurrence, there is, paradoxically, a highly favourable prognosis. Historically, only 40-50% of childhood PTC are driven by genomic variants common to adult PTC, while the oncogenic drivers in the remainder are unknown. METHODS: Fifty-two PPTCs were genotyped by candidate gene testing, followed, in a subset (n=32), by whole exome (WES) and transcriptome sequencing (RNASeq) to identify oncogenic variants in driver-unknown tumors and to explore gene expression patterns. FINDINGS: Thirty-one tumors (60%) carried variants identified by candidate gene testing, specifically: BRAFV600E(22%), NRASQ61R(2%), RET-CCDC6(19%), RET-NCOA4(15%) and PAX8-PPARG. Among the 21 remaining tumors with negative first-line genotyping, oncogenic variants were identified in 19 (90%). These were enriched with oncogenic fusions, with 11 non-recurrent fusion transcripts identified, including two that had not been previously described: STRN-RET and TG-PBF. Most fusions were associated with 3’ receptor tyrosine kinase (RTK) moieties: RET, MET, ALK and NTRK3. Gene expression analysis defined three distinct clusters, which demonstrated differing expression of genes involved in thyroid differentiation and MAPK signaling. INTERPRETATION: Exome and transcriptome-wide sequencing identify genomic drivers of PPTC in 98% of tumors, whereas conventional genotyping leaves many tumors genetically unexplained. PPTC is driven, in the majority of cases, by oncogenic fusion transcripts involving one of a limited number of RTKs. These observations have broad implications for diagnostic testing and for therapeutic options in children with advanced disease.

  Study EGAS00001005187

• WES analysis of a mixed cohort of pituitary tumors

  BACKGROUND: Sporadic pituitary adenomas are the second most common primary brain tumor, but their genetics remains incompletely understood. We performed deep DNA sequencing of 42 pituitary macroadenomas, including hormonally active and atypical adenomas. METHODS: Pituitary adenomas representative of a spectrum of histopathologic subtypes, primary and recurrent tumors, and MIB-1 proliferative indices were profiled. Whole-exome sequencing was performed on tumors and paired blood followed by analysis of somatic copy number alteration and mutation. RESULTS: On the whole, the tumors harbored a low rate of mutations, similar to other benign brain tumors. However, almost one-third (29%) exhibited somatic copy-number alterations across large fractions of the genome (up to 99%). Despite such widespread genomic disruption, these tumors had few focal events, which is atypical among highly-disrupted cancers. The other 71% of tumors formed a distinct molecular class, with somatic copy-number alterations involving less than 6% of the genome. Among the highly-disrupted group, 92% were functional or atypical null-cell adenomas, whereas 87% of the less-disrupted group were clinically non-functional adenomas. The only significantly recurrent events were arm-level losses of chromosomes 1, 2, 11, and 18. No significantly recurrent mutations or insertion/deletion events were identified, suggesting no genes are altered by exomic mutations in >15% of pituitary tumors. CONCLUSIONS: Genomic characterization of a large series of sporadic pituitary adenomas reveals two classes of tumor, based on genomic disruption from arm-level copy-number alterations, which associate with hormonal and histologic subtypes.

  Study EGAS00001001714

• High-Resolution Spatial Transcriptomics Uncover Epidermal-Dermal Divergences in Merkel Cell Carcinoma: Spatial Context Reshapes the Gene Expression Landscape

  Merkel cell carcinoma (MCC) is an aggressive malignancy with neuroendocrine differentiation marked by high plasticity, often manifesting as rapid therapy resistance. Although the cell-of-origin is presumed to be epithelial, epidermal localization of MCC is rarely observed, largely because in situ MCC is typically an incidental finding. Nevertheless, a subset of MCC tumors exhibits epidermotropism, wherein tumor cells are present in the epidermis. The behavior of cancer cells is profoundly influenced by the tumor microenvironment and interactions with neighboring cells. Notably, the normal counterparts of the cancer’s cell-of-origin have been shown to attenuate tumor aggressiveness. Thus, epidermotropic MCC presents a unique opportunity to explore the potential role of surrounding keratinocytes in modulating tumor cell behavior. While the epidermotropic tumor nests share histological resemblance with their dermal counterparts, their transcriptomic profiles remain unexplored. Here, we employed high-definition spatial and single cell transcriptomics to dissect the gene expression profiles of epidermotropic MCC cells, comparing them to MCC cells in the tumor core and those adjacent to blood vessels. Notably, epidermotropic MCC cells exhibit a transcriptomic signature reminiscent of cutaneous squamous cell carcinoma, characterized by upregulation of genes encoding keratins, S100A proteins, as well as calmodulin-like proteins 3 and 5. Mechanistically, this keratinocytic differentiation is associated with enhanced p63 activity, leading to the upregulation of PERP. Collectively, our study demonstrates that MCC cells can adopt a keratinocytic differentiation program in response to microenvironmental cues, underscoring the remarkable phenotypic plasticity of this malignancy and the importance of the tumor microenvironment for tumor characteristics.

  Study EGAS00001008157

• C3 SNPs and outcome after lung transplantation

  Background: Development of chronic rejection is still a severe problem and causes high mortality rates after lung transplantation (LTx). Complement activation is important in the development of acute rejection (AR) and bronchiolitis obliterans syndrome, with C3 as a key complement factor. Methods: Since a single nucleotide polymorphism (SNP) in the C3 gene (rs2230199) is associated with long-term outcome after solid organ transplantation, we investigated this SNP in relation to LTx. In addition, we looked at local production of C3 by analyzing bronchoalveolar lavage fluid (BALF) of LTx patients using isoelectric focusing (IEF). Results: We demonstrated the presence of C3 in BALF and showed that this is produced by the donor lung based on the genotype of SNP rs2230199. We also analyzed donor and patient SNP configurations and observed a significant association between the SNP configuration in patients and episodes of AR during 4-years follow-up. Survival analysis showed a lower AR-free survival in homozygous C3 slow patients (p=0.005). Furthermore, we found a significant association between the SNP configuration in donors and BOS development. Patients receiving a graft from a donor with at least one C3 fast variant for rs2230199 had an inferior BOS-free survival (p=0.044). Conclusions: In conclusion, our data indicate local C3 production by donor lung cells. In addition, a single C3 SNP present in recipients affects short-term outcome after LTx, while this SNP in donors has an opposite effect on long-term outcome after LTx. These results could contribute to an improved risk stratification after transplantation.

  Study EGAS00001003843

• whole genome sequence data of multifocal hepatocellular carcinoma

  Background & Aims: The differentiation of distinct multifocal hepatocellular carcinoma (HCC): multicentric disease vs. intrahepatic metastases, in which the management and prognosis varies substantively, remains problematic. We aim to stratify multifocal HCC and identify novel diagnostic and prognostic biomarkers by performing whole genome and transcriptome sequencing, as part of a multi-omics strategy Methods: A complete collection of tumour and somatic specimens (intrahepatic HCC lesions, matched non-cancerous liver tissue and blood) were obtained from representative patientswith multifocal HCC exhibiting two distinct postsurgical courses. Whole-genome and transcriptome sequencing with genotyping were performed for each tissue specimen to con-trast genomic alterations, including hepatitis B virus integrations, somatic mutations, copy number variations, and structural variations. We then constructed a phylogenetic tree to visualise individual tumour evolution and performed functional enrichment analyses on select differentially expressed genes to elucidate biological processes involved in multifocal HCC development. Multi-omics data were integrated with detailed clinicopathological information to identify HCC biomarkers, which were further validated using a large cohort of HCC patients (n = 174).Results: The multi-omics pro?ling and tumour biomarkers could successfully distinguish the two multifocal HCC types, while accurately predicting clonality and aggressiveness. The dual speci?city protein kinase TTK, which is a key mitotic checkpoint regulator with links to p53 signaling, was further shown to be a promising overall prognostic marker for HCC in the large patient cohort.Conclusions: Comprehensive multi-omics characterisation of multifocal tumour evolution may improve clinical decision-making, facilitate personalised medicine, and expedite identi?cation of novel biomarkers and therapeutic targets in HCC.

  Study EGAS00001002338

• Analysis of recent shared ancestry in a familial cohort identifies coding and noncoding autism spectrum disorder variants

  Autism spectrum disorder (ASD) is a collection of neurodevelopmental disorders characterized by deficits in social communication and restricted, repetitive patterns of behavior or interests. ASD is highly heritable, but genetically and phenotypically heterogeneous, reducing the power to identify causative genes. We performed whole genome sequencing (WGS) in an ASD cohort of 68 individuals from 22 families enriched for recent shared ancestry. We identified an average of 3.07 million variants per genome, of which an average of 112,512 were rare. We mapped runs of homozygosity (ROHs) in affected individuals and found an average genomic homozygosity of 9.65%, consistent with expectations for multiple generations of consanguineous unions. We identified potentially pathogenic rare exonic or splice site variants in 12 known (including KMT2C, SCN1A, SPTBN1, SYNE1, ZNF292) and 12 candidate (including CHD5, GRB10, PPP1R13B) ASD genes. Furthermore, we annotated noncoding variants in ROHs with brain-specific regulatory elements and identified putative disease-causing variants within brain-specific promoters and enhancers for 5 known ASD and neurodevelopmental disease genes (ACTG1, AUTS2, CTNND2, CNTNAP4, SPTBN4). We also identified copy number variants in two known ASD and neurodevelopmental disease loci in two affected individuals. In total we identified potentially etiological variants in known ASD or neurodevelopmental disease genes for ~61% (14/23) of affected individuals. We combined WGS with homozygosity mapping and regulatory element annotations to identify candidate ASD variants. Our analyses add to the growing number of ASD genes and variants and emphasize the importance of leveraging recent shared ancestry to map disease variants in complex neurodevelopmental disorders.

  Study EGAS00001006058

• Somatic mutations predict an aggressive phenotype in meningioma but do not drive grade progression

  A subset of meningiomas progress in histopathological grade and drivers of progression are poorly understood. We aimed to identify somatic mutations and copy number alterations (CNAs) associated with grade progression in a unique matched tumour dataset. Utilising a prospective database, we identified 10 patients with meningiomas that had undergone grade progression and for whom matched pre- and post-progression tissue (n=50 samples) was available for targeted next-generation sequencing. Mutations in NF2 were identified in 4/10 patients, 94% were non-skull base tumours. In one patient, three different NF2 mutations were identified in four tumours. NF2 mutated tumours showed large scale CNAs, with highly recurrent losses in 1p, 10, 22q, and frequent CNAs on chromosomes 2, 3 and 4. There was a correlation between grade and CNAs in a subset of patients. In two patients with sphenoid wing Grade III tumours without detected NF2 mutations, a STAT6-NAB2 inversion was identified. Two patients with tumours without detected NF2 mutations showed a combination of loss and high gain on chromosome 17q. Mutations in SETD2, TP53, TERT promoter and NF2 were not uniform across recurrent tumours and did not correspond with grade progression. Meningiomas that progress in grade generally have a mutational profile already detectable in the pre-progressed tumour, suggesting an aggressive phenotype . CNA profiling shows frequent alterations in NF2 mutated tumours compared to non NF2 mutated tumours and the pattern of CNAs may be associated with grade progression in a subset of cases.

  Study EGAS00001006585

• Microsatellite data show recent demographic expansions in sedentary but not in nomadic human populations in Africa and Eurasia.

  The transition from hunting and gathering to plant and animal domestication was one of the most important cultural and technological revolutions in human history. According to archeologists and paleoanthropologists, this transition triggered major demographic expansions. However, few genetic studies have found traces of Neolithic expansions in the current repartition of genetic polymorphism, pointing rather toward Paleolithic expansions. Here, we used microsatellite autosomal data to investigate the past demographic history of 87 African and Eurasian human populations with contrasted lifestyles (nomadic hunter-gatherers, semi-nomadic herders and sedentary farmers). Likely due to the combination of a higher mutation rate and the possibility to analyze several loci as independent replicates of the coalescent process, the analysis of microsatellite data allowed us to infer more recent expansions than previous genetic studies, potentially resulting from the Neolithic transition. Despite the variability in their location and environment, we found consistent expansions for all sedentary farmers, while we inferred constant population sizes for all hunter-gatherers and most herders, which could result from constraints linked to a nomadic or semi-nomadic lifestyle and/or competition for land between herders and farmers. As an exception, we inferred expansions for Central Asian herders. This might be linked with the arid environment of this area, which may have been more favorable to nomadic herders than to sedentary farmers. Alternatively, current Central Asian herders may descent from populations who have first experienced a transition from hunter-gathering to sedentary agropastoralism, and then a second transition to nomadic herding.

  Study EGAS00001000652

• Clinical relevance of somatic copy-number alterations in plasma circulating tumor DNA from advanced EGFR-mutated NSCLC patients treated with osimertinib

  Background: Somatic copy-number alterations (SCNAs) are associated with drug resistance in patients with advanced epidermal growth factor receptor (EGFR)-mutated non-small cell lung cancer (NSCLC) treated with EGFR-tyrosine kinase inhibitors (TKIs). The purpose of our study was to assess the clinical relevance and development of genome-wide SCNAs in plasma circulating tumor DNA (ctDNA) under osimertinib therapy.Methods: We included 43 patients with advanced EGFR T790M-positive lung adenocarcinoma who were treated with osimertinib after progression under previous EGFR-TKI therapy. We performed genomic profiling of ctDNA in plasma samples from each patient obtained pre-osimertinib and after patients developed resistance to osimertinib. SCNAs were detected by shallow whole-genome plasma sequencing and EGFR mutations were assessed by droplet digital PCR.Results: SCNAs in resistance-related genes (rrSCNAs) were detected in 10 out of 39 (25.6%) evaluable patients before start of osimertinib. The presence of rrSCNAs in plasma before initiation of osimertinib therapy was associated with a lower response rate to osimertinib (50% versus 86%, p=0.02) and was an independent predictor for shorter progression-free survival (adjusted HR 3.26, 95% CI 1.37-7.75, p=0.008) and overall survival (adjusted HR 2.41, 95% CI 1.09-5.34, p=0.03).Conclusions: Genomic profiling of plasma ctDNA is clinically relevant and affects the efficacy and clinical outcome of osimertinib. Our approach allows comprehensive assessment of SCNAs in plasma samples of lung adenocarcinoma patients and may help to guide genotype-specific therapeutic strategies in the future.

  Study EGAS00001004539

• High-throughput Detection of Clinically Relevant Mutations in Archived Tumor Samples By Multiplexed PCR and Next Generation Sequencing

  PURPOSE: Tailoring cancer treatment to tumor molecular characteristics promises to make personalized medicine a reality. However, reliable genetic profiling of archived clinical specimens has been hindered by limited sensitivity and high false positive rates. Here, we describe a novel methodology, MMP-seq, which enables sensitive and specific high-throughput, high content genetic profiling in archived clinical samples. EXPERIMENTAL DESIGN: We first validated MMP-seq’s technical performance in 66 cancer cell lines and a Latin square cross-dilution of known somatic mutations. We next characterized the performance of MMP-seq in 17 formalin-fixed paraffin embedded (FFPE) clinical samples using matched fresh frozen (FF) tissue from the same tumors as benchmarks. To demonstrate the potential clinical utility of our methodology, we profiled FFPE tumor samples from 73 endometrial cancer patients. RESULTS: We demonstrated that MMP-seq enabled rapid and simultaneous profiling of a panel of 88 cancer genes in 48 samples, and detected variants at frequencies as low as 0.4%. We identified DNA degradation and deamination as the main error sources and developed practical and robust strategies for mitigating these issues, and dramatically reduced the false positive rate. Applying MMP-seq to a cohort of endometrial tumor samples identified extensive, potentially actionable alterations in the PI3K and RAS pathways, including novel PIK3R1 hotspot mutations that may disrupt negative regulation of PIK3CA. CONCLUSIONS: MMP-seq provides a robust solution for comprehensive, reliable and high- throughput genetic profiling of clinical tumor samples, paving the way for the incorporation of genomic-based testing into clinical investigation and practice.

  Study EGAS00001000674

• Combination of ribociclib and gemcitabine for the treatment of medulloblastoma

  Group 3 (G3) medulloblastoma (MB) is one of the deadliest forms of the disease for which novel treatment is desperately needed. Here we evaluate ribociclib, a highly selective CDK4/6 inhibitor, with gemcitabine in mouse and human G3 MBs. Ribociclib CNS penetration was assessed by in vivo microdialysis and by immunohistochemistry and gene expression studies. Survival studies to determine the efficacy of ribociclib and gemcitabine combination were performed on mice orthotopically implanted with luciferase labelled mouse and human G3 MB. Pharmacokinetic-pharmacodynamic outcomes and univariable survival models were analyzed to estimate survival. Gene activity inference using NetBID and tumor differentiation analysis investigated the effects of the combination after short and long-term treatments. Tumors from mice treated with oral ribociclib displayed inhibited RB phosphorylation, downregulated E2F target genes, and decreased proliferation. Treatment of mice with the combination of ribociclib and gemcitabine was well tolerated, slowed tumor progression and metastatic spread, and increased survival. Molecular analysis of treated versus untreated tumors showed a significant decrease in the activity and expression of genes involved in cell cycle progression and DNA damage response, and an increase in activity and expression of genes implicated in neuronal identity and neuronal differentiation. Ribociclib is CNS-penetrant. When administered/combined with gemcitabine in orthotopic G3 MB models resulted in improved survival. Our findings, with both mouse and human patient-derived-orthotopic xenograft models, suggest that this combination therapy has promise for children with G3 MB and may represent an effective treatment strategy for other CNS malignancies.

  Study EGAS00001006001

• Study of the rare and low frequency variants in the Saguenay-Lac-Saint-Jean population

  The Saguenay–Lac-Saint-Jean (SLSJ) region is located in northeastern Quebec and is known for its unique demographic history and founder effect. Since founder populations are enriched with population-specific variants, we characterized the variants distribution in SLSJ and compared it with four European populations (Finnish, Sweden, United Kingdom and France), of which the Finnish population is another founder population. Targeted sequencing of the coding and non-coding immune regulatory regions of the SLSJ asthma familial cohort and the four European populations were performed. Rare and low-frequency coding and noncoding regulatory variants identified in the SLSJ population were then investigated for variant- and gene-level associations with asthma and allergy related-traits (eosinophil percentage, immunoglobulin (Ig)E levels and lung function). Our data showed that (1) rare or deleterious variants were not enriched in the two founder populations as compared to the three non-founder European populations; (2) a larger proportion of founder population-specific variants occurred with higher frequencies; and (3) low-frequency variants appeared to be more deleterious. Furthermore, a rare variant, rs1386931, located in the 3’UTR of CXCR6 and intron of FYCO1 was found to be associated with eosinophil percentage. Gene-based analyses identified NRP2, MRPL44 and SERPINE2 to be associated with various asthma and allergy related traits. Our study demonstrated the usefulness of using a founder population to identify new genes associated with asthma and allergy-related traits; thus better understand the genes and pathways implicated in pathophysiology.

  Study EGAS00001003103

• Oncogenic fate conversion by PRDM16s causes acute myeloid leukemia

  Oncogenic mutations confer aberrant replicative capacity to cells with little or no replicative capacity, generating cancer stem cells that perpetuate the tumor through extensive self-replication and differentiation blockade. However, whether oncogenes disrupt the cellular identity ofcancer stem cell by altering the developmental potential is unknown. Fate conversion has been demonstrated by ectopic expression of master transcriptional regulators, such as PU.1 and C/EBP alpha that confers myeloid cell fate to other cell types, and PRDM16 that confers brown adipose fate to white adipocytes or myoblasts. Here we show that a transcriptional regulator overexpressed in human myeloid malignancies, PRDM16s, causes oncogenic fate conversion, by transforming cells fated to form platelets and erythrocytes into myeloid leukemia-initiating cells (LICs). Prdm16s expression in hematopoietic progenitor cells caused a myelodysplastic syndrome (MDS)-like disease that progressed to acute myelogenous leukemia (AML). The myeloid diseases caused by Prdm16s exhibited expansion of megakaryocyte-erythroid progenitors (MEPs) but not granulocyte-macrophage progenitors. MEPs from Prdm16s-induced leukemia possessed LIC potential, and expression of Prdm16s in normal MEPs was sufficient to convert them to myeloid LICs and blocked megakaryocytic/erythroid potential. Prdm16s induced the expression of myeloid master regulators, including PU.1 and C/EBP alpha, by interacting with their super enhancers. Ablation of myeloid master regulators attenuated the myeloid potential and reinstalled the megakaryocytic/erythroid potential of leukemic-MEPs in mouse models and in human AML with PRDM16 rearrangement. Our study demonstrates that oncogenic Prdm16 expression converts the fate of MEPs to a malignant myeloid fate by activating myeloid master transcription factors.

  Study EGAS00001003235

• Implementation of pediatric precision oncology into clinical practice: The individualized Therapies for Children with cancer program “iTHER”

  iTHER is a prospective national precision oncology program aiming to define tumor molecular profiles in children and adolescents with primary very high-risk, relapsed, or refractory pediatric tumors in order to identify relevant aberrations to inform treatment. From April 2017 to April 2021, 302 samples of 253 patients were included. Comprehensive molecular profiling utilizing lcWGS, WES, RNAseq, Affymetrix, and/or 850k methylation profiling was performed for 226 samples. Germline pathogenic variants were identified in 16% of patients. At least one somatic alteration was detected in 204 (90.3%), and 185 (81.9%) were considered druggable, with priority levels very high (6.1%), high (21.3%), moderate (26.0%), intermediate (36.1%), and borderline (10.5%). Diagnosis was revised or refined in 8 patients (3.5%). Spatial and temporal heterogeneity was observed in paired samples of 15 patients, indicating the value of sequential analysis. Of 137 patients with follow-up beyond 12 months, 21 molecularly matched treatments were applied in 19 patients (13.9%) with modest response. Most relevant barriers to not applying targeted therapies included poor performance status, as well as limited access to drugs within clinical trial. iTHER demonstrates feasibility of comprehensive molecular profiling across all ages, tumor types and stages in pediatric cancers, informing of diagnostic, prognostic, and targetable alterations as well as reportable germline variants within a clinically relevant timeframe. Therefore, WES and RNAseq is implemented into standard clinical care at the Princess Máxima Center for all children. Improved access to targeted treatments within biology-driven combination trials is required to eventually improve survival.

  Study EGAS00001007099

• Exome sequencing of uterine leiomyosarcomas

  Uterine leiomyosarcomas (ULMSs) are aggressive smooth muscle tumors associated with poor clinical outcome. Despite previous cytogenetic and molecular studies, their molecular background has remained elusive. To examine somatic variation in ULMS, we performed exome sequencing on 19 tumors. Altogether 43 genes were mutated in at least two ULMSs. Most frequently mutated genes included tumor protein P53 (TP53; 6/19; 33%), alpha thalassemia/mental retardation syndrome X-linked (ATRX; 5/19; 26%), and mediator complex subunit 12 (MED12; 4/19; 21%). Unlike ATRX mutations, both TP53 and MED12 alterations have repeatedly been associated with ULMSs. All the observed ATRX alterations were either nonsense or frameshift mutations. ATRX protein levels were reliably analyzed by immunohistochemistry in altogether 44 ULMSs, and the majority of tumors (23/44; 52%) showed clearly reduced expression. Loss of ATRX expression has been associated with alternative lengthening of telomeres (ALT), and thus the telomere length was analyzed with telomere-specific fluorescence in situ hybridization. The ALT phenotype was confirmed in all ULMSs showing diminished ATRX expression. Exome data also revealed one nonsense mutation in death-domain associated protein (DAXX), another gene previously associated with ALT, and the tumor showed ALT positivity. Aberrant expression of both TP53 and ATRX were associated with poor overall survival. In conclusion, exome sequencing revealed that TP53, ATRX, and MED12 are frequently mutated in ULMSs. ALT phenotype was commonly seen in tumors, indicating that ATR inhibitors, which were recently suggested as possible new drugs for ATRX-deficient tumors, could provide a potential novel therapeutic option for ULMS.

  Study EGAS00001001612

• RB Loss in Resistant EGFR Mutant Lung Adenocarcinomas that Transform to Small Cell Lung Cancer

  Tyrosine kinase inhibitors (TKIs) are effective treatments for non-small cell lung cancers (NSCLCs) with epidermal growth factor receptor (EGFR) mutations. However, they do not lead to cures, and, on average, relapse occurs typically occurs afterafter an average of one year of continuous treatment. In a subset of patients, aA fundamental histological transformation from NSCLC to small cell lung cancer (SCLC) is observed in a subset of the resistant cancers, but the molecular changes associated with this transformation remain unknown. Analysis of a cohort of tumor samples and cell lines derived from resistant EGFR mutant patients with SCLC transformation revealed that RB is lost in 100% of these SCLC transformed cases, but rarely in those that remain NSCLC. Global changes in gene expression, including increased neuroendocrine marker expression and absence of EGFR expression, are observed in resistant EGFR mutant cancers that transformed to SCLC. Consistent Further, with their similarities to classical SCLC at the genetic and gene expression levelincreased neuroendocrine marker and decreased EGFR expression as well as greater sensitivity to BCL2 family inhibition are observed in resistant SCLC transformed cancers compared to resistant NSCLCs. , cell lines derived from resistant EGFR mutant SCLC biopsies are substantially more sensitive to Bcl-2 family inhibition with ABT-263 compared to those derived from resistant EGFR mutant NSCLCs. Together, these findings suggest that despite developing initially as EGFR mutant adenocarcinomas, this subset of resistant cancers ultimately adopt many of the molecular and phenotypic characteristics of classical SCLC.

  Study EGAS00001001102

• Spatiotemporal evolution and inter-patient heterogeneity in primary and recurrent/metastatic head and neck squamous cell carcinoma

  Head and neck cancer squamous cell carcinomas (HNSCCs) are heterogeneous in terms of origin and aetiology. In addition, there is uncertainty about the spatiotemporal genetic evolution from initial diagnosis to recurrence/metastasis (R/M) after primary treatments and further disease progression following systemic treatment. Changes in the genetic profile have implications on the selection of appropriate treatments for patients, especially in the era of targeted therapies and immunotherapies. Here, we analysed a cohort of 11 HNSCC patients with metachronous R/M, among whom nine had paired primary and R/M samples suitable for next-generation sequencing (multiple R/M samples collected at different time points were available for some of the patients). At the genomic level, the R/M samples shared a fraction of the somatic single nucleotide variants (SNVs) with the index primary tumours, but they also acquired many additional mutations, while losing only a few others. A similar behaviour was also observed when examining the changes of mutational signatures between primary and R/M samples. Overall, R/M appeared thus more genetically diverse than the respective primary tumours. The transcriptomic analysis showed that R/M samples had lower immune cell infiltration and several genes related to immune response were significantly downregulated compared to the primary samples. Our results underline the importance of analysing multiple samples per patient to obtain a more complete picture of the patient’s tumour and advocates a re-biopsy in the event of R/M and treatment failure, in order to select the most appropriate therapeutic strategy.

  Study EGAS00001007464

• Neoadjuvant immunotherapy leads to pathological responses in MMR proficient and MMR deficient early stage colon cancers

  PD-1 plus CTLA-4 blockade is highly effective in advanced-stage, mismatch repair (MMR)-deficient (dMMR) colorectal cancers, yet not in MMR-proficient (pMMR) tumors. We postulated a higher efficacy of neoadjuvant immunotherapy in early-stage colon cancers. In the exploratory NICHE study (ClinicalTrials.gov: NCT03026140), patients with dMMR or pMMR tumors received a single dose of ipilimumab and two doses of nivolumab before surgery, the pMMR group with or without celecoxib. The primary objective was safety and feasibility; 40 patients with 21 dMMR and 20 pMMR tumors were treated, and 3 patients received nivolumab monotherapy in the safety run-in. Treatment was well tolerated and all patients underwent radical resections without delays, meeting the primary endpoint. Of the patients who received ipilimumab + nivolumab (20 dMMR and 15 pMMR tumors), 35 were evaluable for efficacy and translational endpoints. Pathological response was observed in 20/20 (100%; 95% exact confidence interval (CI): 86–100%) dMMR tumors, with 19 major pathological responses (MPRs, ≤10% residual viable tumor) and 12 pathological complete responses. In pMMR tumors, 4/15 (27%; 95% exact CI: 8–55%) showed pathological responses, with 3 MPRs and 1 partial response. CD8+PD-1+ T cell infiltration was predictive of response in pMMR tumors. These data indicate that neoadjuvant immunotherapy may have the potential to become the standard of care for a defined group of colon cancer patients when validated in larger studies with at least 3 years of disease-free survival data.

  Study EGAS00001004160

• pancreatic ductal adenocarcinoma exomes in study: Pro-immunogenic Impact of MEK inhibition combined with an anti-CD40 immunostimulatory antibody

  Cancer types with lower mutational load and a non-permissive tumor microenvironment are intrinsically resistant to immune checkpoint blockade. While the combination of cytostatic drugs and immunostimulatory antibodies constitutes an attractive concept for overcoming this refractoriness, suppression of immune cell function by cytostatic drugs may limit therapeutic efficacy. Here we show that targeted inhibition of mitogen-activated protein kinase (MAPK) kinase (MEK) does not impair dendritic cell-mediated T-cell priming and activation. Accordingly, combining MEK inhibitors (MEKi) with agonist antibodies (Abs) targeting the immunostimulatory CD40 receptor resulted in potent synergistic anti-tumor efficacy. Detailed analysis of the mechanism of action of MEKi GDC-0623 by means of flow cytometric analysis of the tumor immune infiltrate and whole tumor transcriptomics showed that, in addition to its cytostatic impact on tumor cells, this drug exerts multiple pro-immunogenic effects, including the suppression of M2-type macrophages, myeloid derived suppressor cells and CD4+ T-regulatory cells. In addition, MEKi was found to induce tumor-cell intrinsic interferon signaling, which contributed to antigen presentation by tumor cells. Finally, the tumoridical impact of MEKi involves the activation of multiple pro-inflammatory pathways involved in immune cell effector function in the tumor microenvironment. Our data therefore indicate that the combination of MEK inhibition with agonist anti-CD40 Ab is a promising therapeutic concept, especially for the treatment of mutant Kras-driven tumors such as pancreatic ductal adenocarcinoma. The pancreatic ductal adenocarcinoma exomes were generated in the context of the Heidelberg Center for Personalized Oncology project HIPO-K28E.

  Study EGAS00001004196

• Genome-Wide Association Study of aspirin-induced PUD in a UK cohort

  Aspirin is a widely prescribed non-steroidal anti-inflammatory drug used as an effective analgesic and at low doses as anti-platelet drug for the prevention of cardiovascular diseases. However, aspirin is associated with a significant incidence of upper gastrointestinal (GI) complications, notably peptic ulceration. This study aims to identify genetic risk factors associated with the pathogenesis of aspirin-induced upper GI ulceration.A cohort of 248 patients with aspirin-induced, endoscopiclly confirmed GI ulceration and 491 controls (ulcer, without aspirin or aspirin without ulcer) were genotyped for 1.524,956 SNPs using the Illumina Omni 2.5 SNP array. A replication cohort (206 NSAID-induced ulcer cases [of which 87 were aspirin] and 309 controls) was subsequently genotyped for identified SNP association signals using Sequenom iPLEX Massarry. Logistic regression analysis of identified two SNP signal with suggestive association with aspirin-induced GI ulceration. A SNP within the EYA1 locus (rs12678747) was found to be associated with aspirin GI ulceration in both the discovery cohort (p=3.79x10-7, OR=2.06 ([1.55-2.73]) and, to a lesser extent, in the replication cohort (p=.0005, OR=1.71[1.17-2.50]. Combined analysis of both cohorts yielded a genome wide significant association (p=1.46x10-8). The rs12678747 genotype was found to be associated with significantly higher EYA1 transcript in ulcer biopsy tissue (p<0.05).Data suggest that EYA1 may represent a risk locus for aspirin-induced GI ulceration. Though it is thought that EYA1 may play a key role in mediating cell death/ survival

  Study EGAS00001002052

• Search for new loci and low-frequency variants influencing glioma risk by exome-array analysis

  To identify protein altering variants (PAVs) for glioma we analysed Illumina HumanExome BeadChip exome array data on 1,882 glioma cases and 8,079 controls from three independent European populations. In addition to single variant tests we incorporated information on the predicted functional consequences of PAVs and analysed sets of genes with a higher likelihood of having a role in glioma on the basis of the profile of somatic mutations documented by large-scale sequencing initiatives. Globally tThere was a strong relationship between effect size and SNPs predicted to be damaging (P=2.29x10-49); however, these variants which are most likely to impact on risk, are rare (MAF<5%). While no single variant showed an association which was statistically significant at the genomewide threshold a number of represented promising Notable associations - were seen with the BRCA2:c.9976A>T, p.(Lys3326Ter)BRCA2 p.Lys332X variant, which has been shown documented to influence breast and lung cancer risks (OR=2.3, P=4.00x10-4 for glioblastoma [GBM]) and IDH2:c.782G>A, p.(Arg261His)IDH2 p.Arg261His (OR=3.21, P=7.67x10-3, for non-GBM). Additionally, gene burden tests revealed a statistically significant association for HARS2 and risk of GBM (P=2.20x10-6). Genome scans of low frequency PAVs represent a complementary strategy to identify disease-causing variants compared with scans based on tagSNPs. Strategies to lessen the multiple testing burden by restricting analysis to PAVs with higher priors affords an opportunity to maximise study power.

  Study EGAS00001001258

• Circumventing intratumoral heterogeneity to identify potential therapeutic targets in hepatocellular carcinoma. Details please contact lifuqiang@genomics and zhaoxin@genomics.cn.

  Background and aims: Intratumoral heterogeneity (ITH) challenges identifying mutations with target therapy potential whereas circulating cell-free DNAs (cfDNAs) could reflect nearly the entire mutation spectrum in given tumors. We investigated how to minimize the limit of ITH for profiling hepatocellular carcinoma (HCC).Methods: Thirty-two multi-regional HCC samples from five patients were subjected to whole exome sequencing (WES) and targeted deep sequencing (TDS). ITH extent was measured by the average percentage of non-ubiquitous mutations (present in parts of tumor regions). Matched cfDNAs were also analyzed by WES and TDS. Profiling efficiencies of single tumor specimen and cfDNA were compared and the one better depicted mutational landscape was selected to screen therapeutic targets.Results: We found variable extents of ITH in HCCs and observed branched and parallel evolution patterns. ITH level decreased at higher sequencing depth of TDS than that measured by WES (28.1% vs 34.9%, P < 0.01) but it remained unchanged upon additional samples analyzed. TDS of single tumor specimen detected an average of 70% the total mutations in HCC. Although more mutations were detected in cfDNA under TDS than WES, an average of 47.2% total HCC mutations uncovered by cfDNA suggested tissue outperform cfDNA and the latter may serve as alternative in profiling HCC genome. Consequently, TDS of single tumor tissue in 66 patients and cfDNAs in four unresectable HCCs identified 38.6% (26/66 and 1/4) patients bearing therapeutic targets.Conclusions: TDS of single tumor specimen could largely circumvent ITH to uncover mutations indicative of target therapy in HCC.

  Study EGAS00001002207

• Aging and genome-wide patterns of DNA methylation in an African rainforest hunter-gathering population

  Aging is associated with widespread changes in genome-wide patterns of DNA methylation. Thousands of CpG sites whose tissue-specific methylation levels are strongly correlated with chronological age have been previously identified. However, the majority of these studies have focused primarily on cosmopolitan populations living in the developed world; it is not known if age-related patterns of DNA methylation at these loci are similar across a broad range of human genetic and ecological diversity. We investigated genome-wide methylation patterns using saliva derived DNA from a traditionally hunting and gathering African populations: the Baka of the western Central African rainforest, together with the ≠Khomani San of the South African Kalahari Desert. We identify hundreds of CpG sites whose methylation levels are significantly associated with age, thousands that are significant in a meta-analysis, and replicate trends previously reported in populations of non-African descent. We confirm that an age associated site in the gene ELOVL2 shows a remarkably congruent relationship with aging in humans, despite extensive genetic and environmental variation across populations. We also demonstrate that genotype state at methylation quantitative trait loci (meQTLs) can affect methylation trends at some known age-associated CpG sites. Our study explores the relationship between CpG methylation and chronological age in populations of African hunter-gatherers, who rely on different diets across diverse ecologies. While many age-related CpG sites replicate across populations, we show that considering common genetic variation at meQTLs further improves our ability to detect previously identified age associations.

  Study EGAS00001002226

• Poly(A) RNA sequencing of hepatocellular carcinoma tumors and their matched noncancerous liver tissues

  Background & Aims: The differentiation of distinct multifocal hepatocellular carcinoma (HCC): multicentric disease vs. intrahepatic metastases, in which the management and prognosis varies substantively, remains problematic. We aim to stratify multifocal HCC and identify novel diagnostic and prognostic biomarkers by performing whole genome and transcriptome sequencing, as part of a multi-omics strategy Methods: A complete collection of tumour and somatic specimens (intrahepatic HCC lesions, matched non-cancerous liver tissue and blood) were obtained from representative patientswith multifocal HCC exhibiting two distinct postsurgical courses. Whole-genome and transcriptome sequencing with genotyping were performed for each tissue specimen to con-trast genomic alterations, including hepatitis B virus integrations, somatic mutations, copy number variations, and structural variations. We then constructed a phylogenetic tree to visualise individual tumour evolution and performed functional enrichment analyses on select differentially expressed genes to elucidate biological processes involved in multifocal HCC development. Multi-omics data were integrated with detailed clinicopathological information to identify HCC biomarkers, which were further validated using a large cohort of HCC patients (n = 174).Results: The multi-omics profiling and tumour biomarkers could successfully distinguish the two multifocal HCC types, while accurately predicting clonality and aggressiveness. The dual specificity protein kinase TTK, which is a key mitotic checkpoint regulator with links to p53 signaling, was further shown to be a promising overall prognostic marker for HCC in the large patient cohort.Conclusions: Comprehensive multi-omics characterisation of multifocal tumour evolution may improve clinical decision-making, facilitate personalised medicine, and expedite identification of novel biomarkers and therapeutic targets in HCC.

  Study EGAS00001002337

• Genomic consequences of aberrant DNA repair mechanisms stratify ovarian cancer histotypes

  Aberrant DNA repair processes are a hallmark of human tumours. How these deficiencies variously impact the genomes of ovarian cancer between and within histotypes remains unknown. We studied the whole genome pointmutation and structural variation patterns of 133 tumours (59 high grade serous (HGSC), 35 clear cell (CCOC), 29 endometrioid (ENOC), and 10 adult granulosa cell tumours (GCT)) as a substrate for class discovery in ovarian cancer. Ab-initio clustering of integrated point mutation and structural variation signatures revealed seven novel subgroups, comprising between and within histotype stratification. Prevalence of fold-back inversions (FBI) co-localised with high level amplification events divided HGSC into two prognostically significant subgroups. This finding was recapitulated in two additional independent cohorts (total n=576 cases), and transcended gene based mutation status and gene expression as prognostically relevant features of HGSC. CCOC were divided into one group harbouring a signature reflective of active genome editing via APOBEC enzymes (26%) and a complementary group with an age-relatedmutational signature. ENOC were divided intomicrosatellite instable (MSI) cases (28%) with a distinctmismatch repairmutation signature and an outlying mutation rate, with the remainder of cases distributed amongst the remaining six groups. We suggest the FBI, MSI, and APOBEC groups represent biologic strata that could direct treatment strategies. Taken together, our work establishes the efficacy of the somatic genome as a biomarker to stratify ovarian cancers, simultaneously identifying patients that may benefit from emerging therapeutics and distinct subgroups within classic ovarian cancer histotypes.

  Study EGAS00001002390

• Recurrent epimutations activate gene body promoters in primary glioblastoma

  Aberrant DNA hypomethylation may play an important role in the growth rate of glioblastoma (GBM), but the functional impact on transcription remains poorly understood. We assayed the GBM methylome with MeDIP-seq and MRE-seq, adjusting for copy number differences, in a small set of non-glioma CpG island methylator phenotype (non-G-CIMP) primary tumors. Recurrent hypomethylated loci were enriched within a region of chromosome 5p15 that is specified as a cancer amplicon and also encompasses TERT, encoding telomerase reverse transcriptase, which plays a critical role in tumorigenesis. Overall, 76 gene body promoters were recurrently hypomethylated, including TERT and the oncogenes GLI3 and TP73. Recurring hypomethylation also affected previously unannotated alternative promoters, and luciferase reporter assays for three of four of these promoters confirmed strong promoter activity in GBM cells. Histone H3 lysine 4 trimethylation (H3K4me3) ChIP-seq on tissue from the GBMs uncovered peaks that coincide precisely with tumor-specific decrease of DNA methylation at 200 loci, 133 of which are in gene bodies. Detailed investigation of TP73 and TERT gene body hypomethylation demonstrated increased expression of corresponding alternate transcripts, which in TP73 encodes a truncated p73 protein with oncogenic function and in TERT encodes a putative reverse transcriptase-null protein. Our findings suggest that recurring gene body promoter hypomethylation events, along with histone H3K4 trimethylation, alter the transcriptional landscape of GBM through the activation of a limited number of normally silenced promoters within gene bodies, in at least one case leading to expression of an oncogenic protein.

  Study EGAS00001000685

• Benchmarking for alignment and variant calling

  Background: Lung cancer is a heterogeneous disease and the primary cause of cancer-related mortality worldwide. Somatic mutations, including large structural variants, are important biomarkers in lung cancer for selecting targeted therapy. Genomic studies in lung cancer have been conducted using short-read sequencing. Emerging long-read sequencing technologies are a promising alternative to study somatic structural variants, however there is no current consensus on how to process data and call somatic events. In this study, we preformed whole genome sequencing of lung cancer and matched non-tumour samples using long and short read sequencing to comprehensively benchmark three sequence aligners and six structural variant callers comprised of four generic callers (Sniffles2, cuteSV, SVIM and DELLY in generic mode) and two somatic callers (nanomonsv and DELLY in somatic modes). Results: Different combinations of aligners and variant callers influenced somatic structural variant detection. The choice of caller had a significant influence on somatic structural variant detection in terms of variant type, size, sensitivity, and accuracy. The performance of each variant caller was assessed by comparing to somatic structural variants identified by short read sequencing. When compared to somatic structural variants detected with short read sequencing, more events were detected with long read sequencing. The mean recall of somatic variant events identified by short read sequencing was higher for the somatic callers (64%) than generic callers (52%), with nanomonsv (69%) showing the best performance. Conclusion: Long read sequencing can identify somatic structural variants. The longer reads have the potential to improve our understanding of cancer development and inform personalized cancer treatment.

  Study EGAS00001007819

• A Functional Network of Gastric-Cancer-Associated Splicing Events Controlled by Dysregulated Splicing Factors

  Except for a few genes known to have oncogenic spliceforms in gastric cancer, the full scope of aberrant splicing in gastric oncogenesis remains unclear. In this study, we elucidated the alternative splicing events associated with gastric carcinogenesis and identified upstream regulators governing these events. We performed RNA-seq on 19 matched tumor/normal pairs and a panel of gastric cancer cell lines, and identified 361 tumor associated (TA) alternative splicing events (ASEs), which included known oncogenic ASEs in genes such as INSR, FGFR2, CD44 and KRAS. We further identified 8 splicing factors dysregulated in gastric tumors, the expression of which are correlated to the splice ratio of the TA ASEs. We thereby constructed a dysregulated splicing network in gastric cancer, consisting of alternative splicing changes correlated to their potential regulators. This network featured three potential regulatory modules centered around the splicing factors ESRP2, MBNL1 and PTPB1. Knockdown of each splicing factor led to the expected changes in splicing of approximately half of the ASEs assayed. Thus, dysregulation of the splicing factors indeed was the cause of many of the TA ASEs. Pathway enrichment analysis showed that genes affected by tumor associated splicing were pre-dominantly involved in cytoskeletal organization. Concordant with this, knockdown of the splicing factors regulating the TA ASEs led to the expected changes in trans-well migration activity. Thus, we have established for the first time a comprehensive splicing network that likely drives gastric tumor progression, based on which novel therapeutic targets can be identified for further evaluation.

  Study EGAS00001002256

• Exome Sequencing Of 75 Individuals From Multiply Affected Coeliac Families

  Coeliac disease is a highly heritable common autoimmune disease involving chronic small intestinal inflammation in response to dietary wheat. The human leukocyte antigen (HLA) region, and 40 newer regions identified by genome wide association studies (GWAS) and dense fine mapping, account for ~40% of the disease heritability. We hypothesized that in pedigrees with multiple individuals with coeliac disease rare [minor allele frequency (MAF) <0.5%] mutations of larger effect size (odds ratios of ~ 2 – 5) might exist. We sequenced the exomes of 75 coeliac individuals of European ancestry from 55 multiply affected families. We selected interesting variants and genes for further follow up using a combination of: an assessment of shared variants between related subjects, a model-free linkage test, and gene burden tests for multiple, potentially causal, variants. We next performed highly multiplexed amplicon resequencing of all RefSeq exons from 24 candidate genes selected on the basis of the exome sequencing data in 2,248 unrelated coeliac cases and 2,230 controls. 1,335 variants with a 99.9% genotyping call rate were observed in 4,478 samples, of which 939 were present in coding regions of 24 genes (Ti/Tv 2.99). 91.7% of coding variants were rare (MAF <0.5%) and 60% were novel. Gene burden tests performed on rare functional variants identified no significant associations (P<1x10-3) in the resequenced candidate genes. Our strategy of sequencing multiply affected families with deep follow up of candidate genes has not identified any new coeliac disease risk mutations.

  Study EGAS00001001093

• Differential Presence of Exons in Cell-Free DNA Reveals Different Patterns in Colorectal Cancer Between Metastatic and Non-Metastatic Patients

  Tumor genotyping during disease follow-up is gaining an increasing interest for monitoring and clinical management in colorectal cancer (CRC). Growing evidences suggest that single biopsies are inefficient for molecular profiling, due to clonal evolution and intratumoral heterogeneity of primary tumors and/or metastasis. Circulating free DNA (cfDNA) has emerged as an alternative source of genetic material that seems to be representative of the continuously changing tumor molecular features. Thus, the so-called 'liquid biopsy' also offers the additional advantages of accessibility and non-invasiveness for the patient. Novel strategies based on next-generation sequencing (NGS) of cfDNA are being currently developed, with most studies focusing on targeted deep sequencing panels of potential clinically actionable genes. Here we describe an alternative approach, consisting of exome-sequencing in cfDNA at a relatively shallower depth to provide a more general overview of the circulating genome. In this study, two groups of colorectal cancer patients, with either disseminated (M) or localized disease (N), were compared using whole-exome sequencing and a RNA-seq bioinformatics analysis pipeline. Thus, a set of 379 exons present at different levels in cfDNA from both groups were identified, giving rise to a completely new concept in this field, termed as 'Differential Presence of Exons'. Differentially present exons (DPEs) were subsequently used for M and N patients clustering and classification, allowing for the design of a predictive algorithm with encouraging preliminary results in a small subset of patients that could not be initially classified attending to the selection criteria (U group, unassignable).

  Study EGAS00001002687

• UK10K NEURO FSZNK

  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 Finnish schizophrenia sample set has been collected from a population cohort using national registers. The entire sample collection consists of 2756 individuals from 458 families of whom 931 are diagnosed with schizophrenia spectrum disorder. Families outside Kuusamo (n=288) all had at least two affected siblings. All diagnoses are based on DSM-IV and for a large fraction of cases there is cognitive data.For further details/descriptions with regard to this data set please contact Tiina Paunio (tiina.paunio@thl.fi)

  Study EGAS00001000119

• UK10K NEURO ASD SKUSE

  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 sample set of UK origin consists of clinically identified subjects with Autism Spectrum Disorders, mostly without intellectual disability (ie. Verbal IQs >70). The subjects represent children and adults with Autism, Asperger syndrome or Atypical Autism, identified according to standardized research criteria (ADI-algorithm, ADOS). A minority has identified comorbid neurodevelopmental disorders (e.g. ADHD). Family histories are available, with measures of broader phenotype in first-degree relatives.For further information on this cohort please contact David Skuse (d.skuse@ucl.ac.uk).

  Study EGAS00001000114

• Sequencing-based counting and size profiling of plasma Epstein-Barr virus DNA enhance population screening of nasopharyngeal carcinoma.

  Circulating tumor-derived DNA testing for cancer screening has recently been demonstrated in a prospective study on identification of nasopharyngeal carcinoma (NPC) among 20,174 asymptomatic individuals. Plasma EBV DNA, a marker for NPC, was detected using real-time PCR. While plasma EBV DNA was persistently detectable in 97.1% of the NPCs identified, ∼5% of the general population had transiently detectable plasma EBV DNA. We hypothesized that EBV DNA in plasma of subjects with or without NPC may have different molecular characteristics. We performed target-capture sequencing of plasma EBV DNA and identified differences in the abundance and size profiles of EBV DNA molecules within plasma of NPC and non-NPC subjects. NPC patients had significantly higher amounts of plasma EBV DNA, which showed longer fragment lengths. Cutoff values were established from an exploratory dataset and tested in a validation sample set. Adopting an algorithm that required a sample to concurrently pass cutoffs for EBV DNA counting and size measurements, NPCs were detected at a positive predictive value (PPV) of 19.6%. This represented superior performance compared with the PPV of 11.0% in the prospective screening study, which required participants with an initially detectable plasma EBV DNA result to be retested within 4 weeks. The observed differences in the molecular nature of EBV DNA molecules in plasma of subjects with or without NPC were successfully translated into a sequencing-based test that had a high PPV for NPC screening and achievable through single time-point testing.

  Study EGAS00001002707

• Neoantigen responses to Immunotherapy in Prostate Cancer

  To determine whether antigen-specific T cell responses can be elicited following treatment with ICT in cancers that have a low TMB, a clinical trial with ipilimumab (anti-CTLA-4) was conducted in 30 patients with metastatic castration-resistant prostate cancer (CRPC). We identified 2 distinct groups of patients divided by clinical outcomes, which consisted of a cohort with Favorable and Unfavorable outcomes. Compared to the Unfavorable cohort, the patients in the Favorable cohort were more likely to have higher intratumoral CD8 T cell density and/or IFN-γ response gene signature, including five patients who had antigen-specific T cell responses that were only detectable after ipilimumab therapy. Three patients had T cell responses against tumor-associated antigens (PSMA and PAP). Importantly, two patients with a relatively low TMB had T cell responses against unique neoantigens. Collectively, our data indicate that a subset of patients with metastatic prostate cancer benefit from ICTs, particularly those who have immunological biomarkers associated with effector T cell responses. Publication: Sumit K. Subudhi, Luis Vence, Hao Zhao, Jorge Blando, Shalini S. Yadav, Qing Xiong, Alexandre Reuben, Ana Aparicio, Paul G. Corn, Brian F. Chapin, Louis L. Pisters, Patricia Troncoso, Rebecca Slack Tidwell, Peter Thall, Chang-Jiun Wu, Jianhua Zhang, Christopher L. Logothetis, Andrew Futreal, James P. Allison and Padmanee Sharma, 2020, Neoantigen responses, immune correlates, and favorable outcomes after ipilimumab treatment of patients with prostate cancer. Science Translational Medicine, 01 Apr 2020:Vol. 12, Issue 537, eaaz3577DOI: 10.1126/scitranslmed.aaz3577. https://pubmed.ncbi.nlm.nih.gov/32238575/

  Study EGAS00001004050

• UK10K_NEURO_MUIR

  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 sample selection consists of subjects with schizophrenia (SZ), autism, or other psychoses all with mental retardation (learning disability). The samples were initially collected under the leadership of Walter J Muir (deceased), now with Prof. Blackwood, Dr McKechanie and Prof McIntosh as custodians. These subjects represent the intersection of severe forms of neurodevelopmental disorders, appear to have a higher rate of familiality of SZ than typical, and are likely to have more serious and penetrant forms of mutations.For further information on this cohort please contact Andrew McIntosh (andrew.mcintosh@ed.ac.uk).

  Study EGAS00001000122

• Characterization of a human iPSC-derived islet differentiation model

  To explore the contribution of islet development to T2D pathogenesis, we characterised the transcriptomes of 3 human iPSC lines (from independent donors) differentiated along the islet development lineage. Whole transcriptome RNA-seq was performed at 7 developmental stages: definitive endoderm, primitive gut tube, posterior foregut (PF), pancreatic endoderm, endocrine progenitors, endocrine-like cells, and beta-like cells (BLC). Differentially-expressed (ΔExp) genes were assigned to the stage in which they were most upregulated (versus baseline iPSC profile), and used to assess the enrichment of T2D GWAS genes, as well as to predict the upstream transcription factors at each developmental stage.We found 9409 ΔExp genes across all stages, including known markers of islet development (NEUROG3, INS). Genes ΔExp in the most-differentiated (BLC) stage were significantly enriched for genes within the 99% credible sets of T2D GWAS loci. Despite this enrichment in BLC only, over 70% of genes mapping within the credible sets were ΔExp before this stage, highlighting the relevance of key T2D GWAS genes in islet development: e.g. expression of TCF7L2 peaked during the PF stage (log2FC=1.2; q=8.5×10-10). REST and LMNA are also highlighted by the analyses as having a potential function in islet development, considering their pattern of expression and those of the genes they regulate. Transcriptomic signatures derived from this iPSC differentiation model highlight T2D-associated GWAS loci, monogenic diabetes genes and potential master regulators of gene networks with plausible function in islet development. All these provide clues for developmental mechanisms contributing to T2D diabetes pathology.

  Study EGAS00001002721

• TRACERx 100: whole exome data of the first 100 TRACERx tumours

  BACKGROUNDTRACERx (TRAcking Cancer Evolution through therapy (Rx)) is a prospective cohort study designed to investigate intratumor heterogeneity (ITH) in relation to clinical outcome, and to determine the clonal nature of driver events and evolutionary processes in early stage non-small cell lung cancer (NSCLC). METHODSMultiregion high-depth whole-exome sequencing (M-seq) was performed on 100 early stage NSCLC tumors resected prior to systemic therapy. A total of 327 tumor regions were sequenced and analyzed to define evolutionary histories, obtain a census of clonal and subclonal events, and assess the relationship between ITH and recurrence-free survival (RFS). RESULTSWidespread ITH was observed for both somatic copy number alterations (median 48% [0.03-88%]) and mutations (median 30% [0.5-93%]). Driver mutations in EGFR, MET, BRAF and TP53 were almost always clonal. However, heterogeneous driver alterations occurring later in evolution were found in over 75% of tumors and were common in PIK3CA, NF1 and genes involved in chromatin modification and DNA response and repair. Genome doubling and ongoing dynamic chromosomal instability (CIN), illustrated by mirrored subclonal allelic imbalance, were identified as causes of ITH resulting in parallel evolution of driver copy number events, including amplifications of CDK4, FOXA1, and BCL11A. Elevated copy number heterogeneity was associated with shorter RFS (HR=4.9, P=0.00044), which remained significant in a multivariate analysis.CONCLUSIONSITH mediated through CIN, rather than point mutational heterogeneity, was associated with increased risk of relapse, supporting its value as a prognostic predictor, and the need to target this high-risk phenotype.

  Study EGAS00001002247

• Breast cancer DNA repair

  Tumors of germline BRCA1/2 mutated carriers show homologous recombination (HR) deficiency (HRD), resulting in impaired DNA double strand break (DSB) repair and high sensitivity to Poly-(ADP-Ribose)-Polymerase (PARP) inhibitors. Although this therapy is expected to be effective beyond germline BRCA1/2 mutated carriers, a robust validated test to detect HRD tumors is lacking. In the present study we therefore evaluated a functional HR assay exploiting the formation of RAD51 foci in proliferating cells after ex vivo irradiation of fresh breast cancers (BrC) tissue: the RECAP test.Methods Fresh samples of 170 primary BrC were analyzed using the RECAP test. The molecular explanation for the HRD phenotype was investigated by exploring BRCA deficiencies, mutational signatures, tumor infiltrating lymphocytes (TILs) and microsatellite instability (MSI).Results RECAP was completed successfully in 148 out of 170 samples (87%). 24 tumors showed HRD (16%), while 6 tumors were HR intermediate (HRi) (4%). HRD was explained by BRCA deficiencies (mutations, promoter hypermethylation, deletions) in 16 cases. Several non-BRCA deficient HRD tumors showed BRCAness mutational signatures, suggesting that they are also bona fide HRD cases. HRD tumors showed an increased incidence of high TIL counts (p=0.023) compared to HR proficient (HRP) tumors and MSI was more frequently observed in the HRD group (2/20, 10%) than expected in BrC (1%) (p=0.017).Conclusion The RECAP test is a robust assay detecting both BRCA1/2 deficient and BRCA1/2 proficient HRD tumors, suggesting that this test identifies approximately 50% more patients that may benefit from PARPi treatment than BRCA gene testing only.

  Study EGAS00001002792

• A_systems_biology_approach_to_understand_immunity_and_pathogenesis_of_malaria_in_children_exposed_to_endemic_Plasmodium_falciparum_transmission

  From long-term studies of children in Kenya, we know that some children have numerous malaria episodes through to older age, while other children become immune more rapidly. Our hypothesis is that these children with frequent malaria episodes are caught in a vicious circle, whereby malaria episodes lead to impaired immunity to malaria, which in turn leads to further episodes of malaria. We will therefore investigate immune responses in children with different life histories of malaria episodes. We will also include a group of children living in similar environmental conditions but who are completely unexposed to malaria. Aims: Using RNAseq we will characterise the transcriptome of a group of children with frequent and febrile malaria or infrequent/asymptomatic. These data will be combined with cytokine profiling data (NIMR) and used to build predictive models (Exeter). The predictions will be validated against a similar second group. A single snap shot of immune responses may represent only the endpoint of many immune processes and may not reflect those that are causally related to differences in malaria outcome. Therefore, we will collect and store samples from a larger group of children who will then be followed up for over 3-4 years. This extended follow up will allow us to identify the same groups as described above, and then examine samples taken before the development of those outcomes. 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 EGAS00001002978

• Biallelic tumor suppressor loss and DNA repair defects in de novo small cell prostate cancer

  Small cell prostatic carcinoma (SCPC) is an aggressive pathology that is managed similar to small cell lung cancer. SCPC can evolve from prostatic adenocarcinoma in response to androgen deprivation therapy but in rare cases is present at initial cancer diagnosis. The molecular etiology of de novo SCPC is incompletely understood due to the scarcity of tumor tissue and the brief life expectancy of patients. Through a retrospective search of our regional oncology pharmacy database, we identified 18 patients diagnosed with de novo SCPC between 2004 and 2017. Ten patients had pure SCPC pathology and the remainder had some admixed adenocarcinoma foci, but all were treated with first line platinum-based chemotherapy. Median overall survival was 28 months. We performed targeted DNA sequencing, whole exome sequencing and mRNA profiling on formalin-fixed paraffin embedded archival tumor tissue. We observed frequent biallelic deletion and/or mutation of tumor suppressors TP53, RB1 and PTEN, similar to treatment related SCPC. Indeed, at the RNA level pure de novo SCPC closely resembled treatment-related SCPC. However, five patients had biallelic loss of DNA repair genes including BRCA1, BRCA2, ATM and MSH2/6: potentially underlying the high genomic instability of this rare disease variant. Two patients with pure de novo SCPC harbored ETS gene rearrangements to androgen-driven promoters, consistent with evolution of de novo SCPC from an androgen-driven ancestor. Overall, our results reveal a highly aggressive molecular landscape underlies this unusual pathologic variant, and propose opportunities for targeted therapy strategies in a disease with few treatment options.

  Study EGAS00001003007

• Open_Targets_020_Epigenomes_of_Cell_Lines

  There is currently a drive to establish cell based assay systems of greater human biological and disease relevance through the use of well characterised transformed cell lines, primary cells and complex cellular models (e.g. co-culture, 3D models). However, although the field is gaining valuable experience in running more non-standard & complex cell assays for target validation and compound pharmacology studies, there is the lack of a systematic approach to determine if this expansion in cell assay models is reflected in increased human biological and disease relevance. The increasing wealth of publically available transcriptomic, and epigenome (ENCODE and Epigenome Roadmap) data represents an ideal reference mechanism for determining the relationship between cell types used for target & compound studies to primary human cells and tissues from both healthy volunteers & patients. The CTTV020 epigenomes of cell line project aims to generate epigenetic and transcriptomic profiles of cell lines and compare these with existing and newly generated reference data sets from human tissue and cell types. The aim is to identify assay systems which will provide greater confidence in translating target biology and compound pharmacology to patients. Multiple cell types commonly used within research have been grouped according to biology. Examples include erythroid, lung epithelial, hepatocyte cell types and immortalised models of monocyte / macrophage biology. 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 EGAS00001003136

• A Single Complex Agpat2 Allele In A Patient With Partial Lipodystrophy

  Genetic lipodystrophies are a group of rare syndromes associated with major metabolic complications - including severe insulin resistance, type 2 diabetes mellitus, and hypertriglyceridemia - which are classified according to the distribution of adipose tissue. Lipodystrophies can be present at birth or develop during life and can range from local to partial and general. With at least 18 different genes implicated so far, definite diagnosis can be challenging due to clinical and genetic heterogeneity. In an adult female patient with clinical and metabolic features of partial lipodystrophy we identified via whole genome sequencing a single complex AGPAT2 allele [V67M;V167A], functionally equivalent to heterozygosity. AGPAT2 encodes for an acyltransferase implicated in the biosynthesis of triacylglycerol and glycerophospholipids. So far homozygous and compound heterozygous mutations in AGPAT2 have only been associated with generalized lipodystrophy. A SNP risk score analysis indicated that the index patient is not predisposed to lipodystrophy based on her genetic background. The partial phenotype in our patient is therefore more likely associated to the genetic variants in AGPAT2. To test whether the resulting double-mutant AGPAT2 protein is functional we analysed its in vitro enzymatic activity via mass spectrometry. The resulting AGPAT2 double mutant is enzymatically inactive. Our data support the view that the current classification of lipodystrophies as strictly local, partial or generalized may have to be re-evaluated and viewed more as a continuum, both in terms of clinical presentation and underlying genetic causes. Better molecular understanding of lipodystrophies may lead to new therapies to treat adipose tissue dysfunction in common and rare diseases.

  Study EGAS00001003177

• Structure and evolution of double minutes in diagnosis and relapse brain tumors

  Double minute chromosomes are extrachromosomal circular DNA fragments frequently found in brain tumors. To understand their evolution, we characterized the double minutes in paired diagnosis and relapse tumors from a pediatric high-grade glioma and four adult glioblastoma patients. We determined the full structures of the major double minutes using a novel approach combining multiple types of supporting genomic evidence. Among the double minutes identified in the pediatric tumor, only one carrying EGFR was maintained at high abundance in both samples, whereas two others were resent in only trace amounts at diagnosis but abundant at relapse, and the rest were found either in the relapse sample only or in the diagnosis sample only. For the EGFR-carrying double minutes, we found a secondary somatic deletion in all copies at relapse, after erlotinib treatment, but this was rare at diagnosis, suggesting potential resistance to the EGFR inhibitor. This mutation caused an in-frame RNA transcript to skip exon 16, a novel transcript absent in EST database and in about 700 brain tumors that we reviewed. We observed similar patterns involving longitudinal copy number shift of double minutes in another four pairs (diagnosis / relapse) of adult glioblastoma. Overall, in three of five paired tumor samples, we found that although the same oncogenes were amplified at diagnosis and relapse, they were amplified on different double minutes. Our results suggest that double minutes readily evolve, increasing tumor heterogeneity rapidly. Understanding patterns of double minute evolution can shed light on future therapeutic solutions to brain tumors carrying such variants.

  Study EGAS00001003212

• Longitudinal_profiling_of_the_immune_response_to_Plasmodium_vivax_in_naive_hosts_by_RNA_sequencing

  Plasmodium vivax offers unique challenges for control and elimination, and may prove a tougher hurdle to overcome than Plasmodium falciparum. And yet compared to P. falciparum we know very little about the innate and adaptive immune responses that need to be harnessed to reduce disease and transmission. We recently generated a blood bank of a new clonal field isolate of P. vivax (PvW1) for human challenge studies and used systems immunology tools to track the host response throughout infection and convalescence. As part of this study, RNA-sequencing was used to resolve changes in whole blood gene expression through time in 6 volunteers (7-9 time-points per volunteer). In summary, these data show that P. vivax induces two distinct transcriptional programmes in whole blood during and after infection. During infection, transcriptional profiling reveals the rapid mobilisation of an emergency myeloid response, which leads to systemic inflammation and the recruitment of all major T cell subsets into lymphoid tissues. Six days after infection, this innate response subsides and a transcriptional signature of proliferation is revealed. This most likely represents widespread activation of lymphocytes, which return to the circulation after parasite clearance - transcriptional profiling of T cells at this time-point could therefore reveal the outcomes of critical cell-cell interactions that take place within the spleen during infection. 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 EGAS00001003847

• MYCN Amplification and ATRX Mutations are Incompatible in Neuroblastoma

  Age and stage of disease at diagnosis is associated with outcome in neuroblastoma. A previous study identified ATRX point mutations and in-frame deletions in older patients with stage 4 disease. To validate the frequency of ATRX mutations and assess the relation of mutations and known prognostic variables, we collaborated with the Children’s Oncology Group and sequenced the entire ATRX genomic locus in 475 neuroblastoma tumor samples to validate the frequency of ATRX mutations and assess the relation of mutations and known prognostic variables. Older age at diagnosis, stage 4 disease, and unfavorable histology were significantly associated with mutations in the ATRX gene; MYCN amplification was mutually exclusive from ATRX mutations. To directly test if ATRX mutations and MYCN amplification are incompatible in neuroblastoma, we mutated the ATRX gene in MYCN amplified neuroblastoma cell lines and ectopically expressed MYCN in ATRX mutant neuroblastoma cell lines. Both approaches showed that ATRX and MYCN amplification are incompatible in neuroblastoma in culture and in vivo. To gain a better understanding of the underlying molecular and cellular mechanisms of this incompatibility, we performed a series of experiments including telomere analysis, epigenetic profiling, electron microscopic analysis and metabolic profiling. We discovered that induction of MYCN expression in ATRX mutant neuroblastomas cells leads to metabolic reprogramming, mitochondrial dysfunction, increased reactive oxygen species and replicative stress. We propose that this leads to synthetic lethality in ATRX mutant neuroblastomas because of the underlying replicative stress in those cells as a result of dysfunction of this essential histone chaperone.

  Study EGAS00001003257

• These are from Korean HCC samples with exome sequencing

  Hepatic resection is the most curative treatment option for early-stage hepatocellular carcinoma, but is associated with a high recurrence rate, which exceeds 50% at 5 years after surgery. Understanding the genetic basis of hepatocellular carcinoma at surgically curable stages may enable the identification of new molecular biomarkers that accurately identify patients in need of additional early therapeutic interventions. Whole exome sequencing and copy number analysis was performed on 231 hepatocellular carcinomas (72% with hepatitis B viral infection) that were classified as early-stage hepatocellular carcinomas, candidates for surgical resection. Recurrent mutations were validated by Sanger sequencing. Unsupervised genomic analyses identified an association between specific genetic aberrations and postoperative clinical outcomes. Recurrent somatic mutations were identified in 9 genes, including TP53, CTNNB1, AXIN1, RPS6KA3, and RB1. Recurrent homozygous deletions in FAM123A, RB1, and CDKN2A, and high-copy amplifications in MYC, RSPO2, CCND1, and FGF19 were detected. Pathway analyses of these genes revealed aberrations in the p53, Wnt, PIK3/Ras, cell cycle, and chromatin remodelling pathways. RB1 mutations were significantly associated with cancer-specific and recurrence-free survival after resection (p = 0.016 and p = 0.001, respectively). FGF19 amplifications, known to activate Wnt signalling, were mutually exclusive with CTNNB1 and AXIN1 mutations, and significantly associated with cirrhosis (p = 0.017). RB1 mutations can be used as a prognostic molecular biomarker for resectable hepatocellular carcinoma. Further study is required to investigate the potential role of FGF19 amplification in driving hepatocarcinogenesis in patients with liver cirrhosis and to investigate the potential of anti-FGF19 treatment in these patients.

  Study EGAS00001000604

• Chromatin immunoprecipitation linked to next-generation whole genome sequencing (ChIP-Seq) for H3K36me3 in paediatric high grade glioma cell lines KKNS4 and SF188 with and without a G34V mutation in H3F3A

  Glioblastomas of children and young adults have a median survival of only 12-15months and are clinically and biologically distinct from histologically similar cancers in older adults1. They are defined by highly specific mutations in the gene encoding the histone H3.3 variant H3F3A2, occurring either at or close to key residues marked by methylation for regulation of transcription – K27 and G34. We performed chromatin immunoprecipitation linked to next-generation whole genome sequencing (ChIP-Seq) for H3K36me3 in order to test the hypothesis that, rather than total H3K36me3, the G34V mutation may instead result in differential binding of the trimethyl mark throughout the genome. Compared with H3F3A wild-type SF188 paediatric GBM cells, H3K36me3 was found to be significantly differentially bound in KNS42 cells at 5130 distinct regions of the genome corresponding to 156 genes. Concurrent whole genome DNA sequencing demonstrated that this was not confounded by copy number alterations. The cerebral hemispheric-specific G34 mutation drives a distinct expression signature through differential genomic binding of the K36 trimethylation mark (H3K36me3). The transcriptional program induced recapitulates that of the developing forebrain, and involves numerous markers of stem cell maintenance, cell fate decisions and self-renewal. Critically, H3F3A G34 mutations cause profound upregulation of MYCN, a potent oncogene which is causative of glioblastomas when expressed in the correct developmental context. This driving aberration is selectively targetable in this patient population by inhibiting kinases responsible for stabilisation of the protein.

  Study EGAS00001001437

• Spatial and temporal evolution of distal 10q deletion, a prognostically unfavorable event in diffuse low-grade gliomas

  Background: The disease course of patients with diffuse low-grade glioma is notoriously unpredictable. Temporal and spatially distinct samples may provide insight into the evolution of clinically relevant copy number aberrations (CNAs). The purpose of this study is to identify CNAs that are indicative of aggressive tumor behaviour and can thereby complement the prognostically favorable 1p/19q co-deletion. Results: Genome-wide, 50 base pair single-end, sequencing was performed to detect CNAs in a clinically well-characterized cohort of 98 formalin-fixed paraffin-embedded low-grade gliomas. CNAs are correlated with overall survival as an endpoint. Seventy-five additional samples from spatially distinct regions and paired recurrent tumors of the discovery cohort were analysed to interrogate the intratumoral heterogeneity and spatial evolution. Loss of 10q25.2-qter is a frequent subclonal event and significantly correlates with an unfavorable prognosis. A significant correlation is furthermore observed in a validation set of 126 and confirmation set of 184 patients. Loss of 10q25.2-qter arises in a longitudinal manner in paired recurrent tumor specimens, whereas the prognostically favorable 1p/ 19q co-deletion is the only CNA that is stable across spatial regions and recurrent tumors. Conclusions: CNAs in low-grade gliomas display extensive intratumoral heterogeneity. Distal loss of 10q is a late onset event and a marker for reduced overall survival in low-grade glioma patients. Intratumoral heterogeneity and higher frequencies of distal 10q loss in recurrences suggest this event is involved in outgrowth to the recurrent tumor.

  Study EGAS00001000643

• Study on the consequences of prenatal famine exposure on DNA methylation.

  The potential of the epigenome to undergo environmentally-induced changes during early development is central to its postulated role in human disease. Here, we report a genome-scale view of DNA methylation differences after early gestational malnutrition that was caused by the Dutch Hunger Winter, a severe famine at the end of World War II. We adopted an annotation guided analysis of reduced-representation bisulfite sequencing data generated in 24 middle-aged individuals who had been exposed to famine in utero and 24 unexposed sibling controls. We observed differential DNA methylation primarily at genomic regions with regulatory potential, including enhancers active during early development. In-depth annotation of 181 prenatally-induced differentially methylated regions (P-DMRs) showed that they were predominantly located in gene bodies and co-occurred with histone marks denoting an active state. In line with the early gestational timing of the exposure, the P-DMRs mapped to genes that were enriched for differential expression patterns during blastocyst development and early organogenesis. Validation and further explorative analyses of 6 P-DMRs mapping to SMAD7, CDH23, INSR, RFTN1, CPT1A and KLF13 highlighted the critical role of gestational timing, indicated that differential methylation extended along pathways involved in growth, development and metabolism, and suggested associations with birth weight and serum LDL cholesterol in adulthood and confirmed an enhancer function for the INSR and CPT1A P-DMRs. Our results point toward a link between prenatal malnutrition and epigenetic modulation of growth and lipid metabolism that may underlie the adverse metabolic phenotype of exposed individuals in later life.

  Study EGAS00001000668

• Dilgom_Exome

  The Finrisk sample sets are part of the National FINRISK Study. It is a large population survey on risk factors of chronic, noncommunicable diseases. The survey is carried out since 1972 every five years using independent, random and representative population samples from different parts of Finland. The main results from the previous FINRISK 2007 survey are published.The National FINRISK Study Survey was carried out in 5 areas in Finland and 2000 inhabitants aged 25-75 years were invited to participate in each year. Among findings were that Finns continue to gain weight.Data from FINRISK surveys are used for many different research projects and for national health monitoring needs. The recent research activities deal, in addition to cardiovascular diseases and the classical risk factors, also with e.g. asthma and allergy, alcohol, socioeconomic factors and genetic epidemiology.The FINRISK study is part of the MORGAM Project (MONICA Risk, Genetics, Archiving and Monigraph), sponsored by the EU and MDECODE (Molecular Diversity and Epidemiology of Common Disease) program coordinated by the University of Michigan.The exome sequencing study will be part of the Dilgom study, which is a part of the larger Finrisk population based health study performed in Finland. It consists of 5000 individuals with a prospective aspect of metabolic traits. The cohort has been extensively phenotyped for their cardiovascular and metabolic status. So far, we have performed a 660K Illumina GWAS and a full genome wide expression study from peripheral blood cells of 500 individuals. The cohort has also been in total genotyped by the cardiometabochip

  Study EGAS00001000086

• Quantification of chromosomal copy number aberrations by shallow whole-genome sequencing

  Detection of DNA copy number aberrations by shallow whole-genome sequencing (WGS) faces many challenges including lack of completion and errors in the human reference genome, repetitive sequences, polymorphisms, variable sample quality, and biases in the sequencing procedures.Formalin-fixed paraffin-embedded (FFPE) archival material, the analysis of which is important for studies of cancer, presents particular analytical difficulties due to degradation of the DNA and frequent lack of matched reference samples. We present a robust, cost-effective WGS method for DNA copy number analysis that addresses these challenges more successfully than currently available procedures. In practice very useful profiles can be obtained with 0.1 fold genome coverage. We improve on previous methods by; first, implementing a combined correction for sequence mappability and GC content, and second, applying this procedure to sequence data from the 1000 Genomes Project in order to develop a blacklist of problematic genome regions. A small subset of these blacklisted regions were previously identified by ENCODE, but the vast majority are novel unappreciated problematic regions. Our procedures are implemented in a pipeline called QDNAseq. We have analyzed over 1,000 samples, most of which were obtained from the fixed tissue archives of over 25 institutions.We demonstrate that for most samples our sequencing and analysis procedures yield genome profiles with noise levels near the statistical limit imposed by read counting. The described procedures also provide better correction of artifacts introduced by low DNA quality than prior approaches, and better copy number data than high-resolution microarrays at substantially lower cost.

  Study EGAS00001000642

• Prediction and quantification of splice events from RNA-seq data

  Analysis of splice variants from short read RNA-seq data remains a challenging problem. Here we present a novel method for the genome-guided prediction and quantification of splice events from RNA-seq data, which enables the analysis of unannotated and complex splice events. Splice junctions and exons are predicted from reads mapped to a reference genome and are assembled into a genome-wide splice graph. Splice events are identified recursively from the graph and are quantified locally based on reads extending across the start or end of each splice variant. We assess prediction accuracy based on simulated and real RNA-seq data, and illustrate how different read aligners (GSNAP, HISAT2, STAR, TopHat2) affect prediction results. We validate our approach for quantification based on simulated data, and compare local estimates of relative splice variant usage with those from other methods (MISO, Cufflinks) based on simulated and real RNA-seq data. In a proof-of-concept study of splice variants in 16 normal human tissues (Illumina Body Map 2.0) we identify 249 internal exons that belong to known genes but are not related to annotated exons. Using independent RNA samples from 14 matched normal human tissues, we validate 9/9 of these exons by RT-PCR and 216/249 by paired-end RNA-seq (2 x 250 bp). These results indicate that de novo prediction of splice variants remains beneficial even in well-studied systems. An implementation of our method is freely available as an R/Bioconductor package SGSeq.

  Study EGAS00001001026

• Comparison of the diagnostic yield of aCGH and NGS across different neurodevelopmental disorders

  Current recommendations for the genetic diagnosis of neurodevelopmental disorders (NDDs) do not include the use of next generation sequencing (NGS) and are still based on chromosomal microarrays, such as comparative genomic hybridization array (aCGH). The present study compares the diagnostic yield obtained by aCGH and clinical exome sequencing in NDD globally and its spectrum of disorders. To that end, 1412 patients clinically diagnosed with NDDs and studied with aCGH were classified into phenotype categories: Global developmental delay/intellectual disability (GDD/ID); autism spectrum disorder (ASD); and Other NDDs. These categories were further subclassified into subcategories, based on the most frequent accompanying signs and symptoms: isolated forms, forms with epilepsy; forms with micro or macrocephaly and syndromic forms. A subset of 245 patients of the 1412 were subjected to clinical exome sequencing. Diagnostic yield of aCGH and clinical exome sequencing, expressed as the number of solved cases, was compared for each phenotype category and subcategory. Clinical exome sequencing was superior than aCGH for all cases except for isolated ASD, with no cases solved by NGS. Globally, clinical exome sequencing solved 20% of cases (versus 5.7% by aCGH) and the diagnostic yield was highest for all forms of GDD/ID and lowest for Other NDDs (7.1% versus 1.4% by aCGH) and ASD (6.1% versus 3% by aCGH). In the majority of cases, diagnostic yield was higher in the phenotype subcategories than in the mother category. These results support the use of NGS as the first-tier test in the diagnostic algorithm of all NDDs followed by aCGH when necessary.

  Study EGAS00001004949

• Expanding the neurodevelopmental phenotypes of individuals with de novo KMT2A variants

  De novo loss-of-function (LoF) variants in the KMT2A gene are associated with Wiedemann-Steiner Syndrome (WSS). Recently, de novo KMT2A variants have been identified in sequencing studies of cohorts of individuals with neurodevelopmental disorders (NDDs). However, most of these studies lack the detailed clinical information required to determine whether those individuals have isolated NDDs or WSS (i.e. syndromic NDDs). We performed thorough clinical and neurodevelopmental phenotyping on six individuals with de novo KMT2A variants. From these data, we found that all six patients met clinical criteria for WSS and we further define the neurodevelopmental phenotypes associated with KMT2A variants and WSS. In particular, we identified a subtype of Autism Spectrum Disorder (ASD) in five individuals, characterized by marked rigid, repetitive and inflexible behaviours, emotional dysregulation, externalizing behaviours, but relative social motivation. To further explore the clinical spectrum associated with KMT2A variants, we also conducted a meta-analysis of individuals with KMT2A variants reported in the published literature. We found that de novo LoF or missense variants in KMT2A were significantly more prevalent than predicted by a previously established statistical model of de novo mutation rate for KMT2A. Our genotype-phenotype findings better define the clinical spectrum associated with KMT2A variants and suggest that individuals with de novo LoF and missense variants likely have a clinically unrecognized diagnosis of WSS, rather than isolated NDD or ASD alone. This highlights the importance of a clinical genetic and neurodevelopmental assessment for individuals with such variants in KMT2A.

  Study EGAS00001003521

• Prediction of pigmentation phenotypes by SNP typing in a Northern German population

  Human pigmentation traits are of great interest to many research areas, from ancient DNA analysis to forensic science. We aimed to develop a gene-based predictive model for pigmentation phenotypes in a realistic target population for forensic case work from Northern Germany. Our aim was to determine whether better prediction accuracy can be achieved, or fewer genetic markers may suffice, than in previously studied, genetically more heterogeneous populations. We investigated the association between eye, hair and skin colour, and 12 candidate single nucleotide polymorphisms (SNPs) from six genes. Our study comprised two samples of 300 and 100 individuals from Northern Germany who were carefully characterized with regard to pigmentation phenotypes. The first sample was used to select trait-associated SNPs whereas the second sample served to estimate odds ratios (ORs) and to quantify the predictive capability of the respective SNP genotypes. SNP rs12913832 in HERC2 was found to be strongly associated with blue eye colour (OR=15.6, p<1.2•10-4) and to yield reasonable predictive power (90% sensitivity, 63% specificity). SNP associations with hair and skin colour were weaker and genotypes less predictive. A comparison to two recently published sets of markers to predict eye and hair colour revealed that the consideration of additional SNPs with weak to moderate effect increases the predictive power in Northern Germans for eye colour, but not for hair colour. In addition, fine phenotyping and differentiation of hair colour (light / dark and red tint / no red tint) were found to increase the number of significant genotype-phenotype associations.

  Study EGAS00001001174

• Single-cell DNA amplicon sequencing reveals clonal heterogeneity and evolution in T-cell acute lymphoblastic leukemia

  T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive leukemia that is most frequent in children and is characterized by the presence of few chromosomal rearrangements and 10 to 20 somatic mutations in protein-coding regions at diagnosis. The majority of T-ALL cases harbor activating mutations in NOTCH1 together with mutations in genes implicated in kinase signaling, transcriptional regulation or protein translation. To obtain more insight in the level of clonal heterogeneity at diagnosis and during treatment, we used single-cell targeted DNA sequencing with the Tapestri platform. We designed a custom ALL panel and obtained accurate single-nucleotide variant and small insertion-deletion mutation calling for 305 amplicons covering 110 genes in about 4400 cells per sample and time point. A total of 108,188 cells were analyzed for 25 samples of 8 T-ALL patients. We typically observed a major clone at diagnosis (>35% of the cells) accompanied by several minor clones of which some were less than 1% of the total number of cells. Four patients had >2 NOTCH1 mutations some of which present in minor clones, indicating a strong pressure to acquire NOTCH1 mutations in developing T-ALL cells. By analyzing longitudinal samples, we detected the presence and clonal nature of residual leukemic cells as well as clones with a minor presence at diagnosis that evolved to clinically relevant major clones at later disease stages. Single-cell DNA amplicon sequencing is a sensitive assay to detect clonal architecture and evolution in T-ALL.

  Study EGAS00001004440

• Associations between APOE status and cognitive ability across the lifecourse

  The association between APOE genotype and cognitive function suggests a positive role for the e2 allele and a negative role for the e4 allele. Both alleles have relatively low frequencies in the general population, hence meta-analyses have been based on many small, heterogeneous studies. Here we report the APOE-cognition associations in the largest single analysis to date. APOE status and cognitive ability were measured in 18,337 participants from the Generation Scotland study between 2006 and 2011. The age range was 18-94 years with a mean of 47 (SD 15). Four cognitive domains were assessed: verbal declarative memory (paragraph recall), processing speed (digit symbol substitution), verbal fluency (phonemic verbal fluency), and vocabulary (Mill Hill synonyms). Linear regression was used to assess the associations between APOE genetic status and cognition. Possession of the e4 allele was associated with lower scores on the measures of memory and processing speed in subjects aged >60. Across all age ranges, the e4 allele was linked to better verbal fluency scores. In younger subjects (≤60 years) the e4 allele was linked to higher vocabulary scores. There were no associations between the e2 allele and cognitive ability. As seen in previous meta-analyses, the APOE e4 allele is linked to poorer cognitive performance in the domains of memory and processing speed. By contrast, positive associations were seen between the e4 allele and measures of verbal fluency and vocabulary. All associations were relatively small and, in many cases, nominally significant despite the very large sample size.

  Study EGAS00001001235

• Transcriptome_human_nasal_epithelium

  1. Odors are detected, firstly, by olfactory sensory neurons (OSNs) in the olfactory epithelium of the nose. This neurons then project directly to the olfactory bulb in the brain. Olfaction depends on cellular regeneration of the OE, olfactory bulb and hippocampus, and on their continual re-wiring. The olfactory neural pathway includes regions of the frontal, temporal and limbic brain, which in turn overlap with brain areas involved in brain disorders. OSNs are the only aspect of the human brain exposed to the external environment. This not only makes them vulnerable to environmental changes, but also accessible for biomedical studies.We have already sequenced and developed a protocol for analyzing the transcriptome of mouse main olfactory epithelium and single OSNs. We propose here to perform a similar study for samples from the human olfactory epithelium. We have developed a minimally invasive method for obtaining human OSNs, among other cells from the nasal epithelium. In this experiment, we have obtained cell samples from the olfactory epithelium, including OSN, from healthy volunteers. We would like to further characterize them by RNA sequencing. This will give us valuable insight into human olfaction. It will also provide a first step into a new avenue to study, and find biomarkers for, brain diseases though the analysis of these easily available neurons. 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 EGAS00001001294

• UK10K_RARE_THYROID

  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.Two cohorts of subjects are being analysed: Individuals with Congenital Hypothyroidism (CH) due either to dysgenesis or dyshormonogenesis; Patients with Resistance to Thyroid hormone (RTH), a disorder characterized by elevated thyroid hormones and variable tissue refractoriness to hormone action. The CH cohort has been enriched for genetic aetiologies by recruiting cases that are familial, on a consanguineous background or syndromic. The RTH cohort consists of cases in which candidate gene analyses have been negative.For further information with regard to this cohort please contact Krishna Chatterjee (kkc1@medschl.cam.ac.uk).

  Study EGAS00001000131

• A unidirectional histone code in bidirectional promoters across cell types

  Recent studies indicate that thousands of genes are expressed from bidirectional promoters (BDPs). Gene regulation at BDPs is poorly understood, in particular how the cell is able to regulate them differently. Here we investigate the effect of histone Modifications in BDP regulation. In this study, we model the gene activity using different gene expression assays, such as RNA-Seq, GRO-cap, and CAGE around the BDP transcription start sites (TSSs) using different histone modifications in various cell types. We develop a new statistical approach that links histone modifications to gene expression at BDPs. It improves over previous methods, because it is able to capture spatial dependencies of histone modifications along a promoter and leads to more interpretable results. We predict a general histone code that is independent of transcript orientation, cell type, and promoter configuration. The histone code at BDPs reveals that promoters are regulated unidirectionally, such that the majority of histone marks associated with gene expression occur downstream of the gene's TSS. By contrasting associations of histone marks with steady-state levels of capped RNAs in CAGE and nascent capped RNAs in GRO-cap, we show which histone marks have a preference for initiating polymerases and actively elongating polymerases. Using single-cell data we show that the bidirectional histone signal of activating marks is often due averaging of a heterogeneous cell population. Our results have implications for other experiments where the relationship between histone modification and gene initiation or gene elongation is investigated on a genome-wide scale.

  Study EGAS00001001656

• Healthspan and lifespan extension by fecal microbiota transplantation in progeroid mice

  The gut microbiome is emerging as a key regulator of several metabolic, immune and neuroendocrine pathways. Gut microbiome deregulation has been implicated in major conditions such as obesity, type 2 diabetes, cardiovascular disease, non-alcoholic fatty acid liver disease or cancer, but its precise role in aging remains to be elucidated. Here, we characterize the gut microbiome profile of accelerated aging and show that its external modulation is sufficient to extend healthspan and lifespan. We found that two different mouse models of progeria present with intestinal dysbiosis, which is characterized, among other alterations, by an increment in proteobacteria, cyanobacteria and a loss in verrucomicrobia, whereas long-lived humans (i.e., centenarians) exhibit a remarkable increase in verrucomicrobia and a reduction in proteobacteria. Fecal microbiota transplantation proved to be effective to enhance healthspan and lifespan in both progeroid mouse models and, more importantly, the solely transplantation with the verrucomicrobia Akkermansia muciniphila was sufficient to exert beneficial effects. Our results demonstrate that intestinal dysbiosis is a phenotype associated with accelerated aging and its correction provides health benefits. Moreover, metabolomic analysis of ileal content points to the restoration of secondary bile acids as a possible mechanism for the beneficial outcome of reestablishing a healthy microbiome. Our results suggest the existence of a link between gut aging and the microbiota, and can help to gain insight into the rationale for microbiome-based interventions against age-related diseases. Mouse sequence data have been deposited in ENA (https://www.ebi.ac.uk/ena) under accession number PRJEB34214

  Study EGAS00001003656

• Detection of ctDNA in Plasma of Patients with Clinically Localised Prostate Cancer is Associated with Rapid Disease Progression

  Purpose: DNA originating from degenerate tumour cells can be detected in the circulation in many tumour types, where it can be used as a marker of disease burden as well as to monitor treatment response. Although circulating tumour DNA (ctDNA) measurement has prognostic/predictive value in metastatic prostate cancer, its utility in localised disease is unknown. Patients and Methods: We performed whole genome sequencing of tumour-normal pairs in eight patients with clinically localised disease undergoing prostatectomy, identifying high confidence genomic aberrations. A bespoke DNA capture and amplification panel against the highest prevalence, highest confidence aberrations for each individual was designed and used to interrogate ctDNA isolated from plasma prospectively obtained pre- and post- (24 hours and 6 weeks) surgery. In a separate cohort (n=189) we identified the presence of ctDNA TP53 mutations in pre-operative plasma in a retrospective cohort, and determined its association with biochemical- and metastasis-free survival. Results: Tumour variants in ctDNA were positively identified pre-treatment in two of eight patients, which in both cases remained detectable postoperatively. Patients with tumour variants in ctDNA had extremely rapid disease recurrence and progression compared to those where variants could not be detected. In terms of aberrations targeted, single nucleotide and structural variants outperformed indels and copy number aberrations. Detection of ctDNA TP53 mutations was associated with a significantly shorter metastasis-free survival (6.2 vs. 9.5 years (HR 2.4; 95% CIs 1.2-4.8, p = 0.014)). Conclusions: CtDNA is uncommonly detected in localised prostate cancer but its presence portends more rapidly progressive disease.

  Study EGAS00001004636

• Molecular classification improves risk assessment in adult B-lineage ALL: Patients on the international UKALLXII-ECOG2993 trial.

  Acute lymphoblastic leukemia (ALL) in adults is a rare disease with outcomes that remain significantly inferior to those of children with ALL. While genomic classification of pediatric B-lineage ALL (B-ALL) has led to improved risk assignment and therapeutic paradigms, similar advances in adults have been lacking. The international UKALLXII/E2993 (NCT00002514) trial accrued over 2000 adult (17-64 years old) patients over 13 years, including 1229 BCR/ABL negative B-ALLs. While 93% of B-ALL achieved a remission, 41% relapsed at a median time of 13 months (Range) and overall survival (OS) at 5 years was 42% (CI). Older age and higher white blood cell (WBC) count at presentation negatively affected OS. Patients with a sibling donor assigned to allogeneic transplantation had superior OS compared to patients without donor. Randomization of patients without a donor showed an advantage of chemotherapy over autologous transplantation in first remission. Transcriptome sequencing of 94% of E2993 patients enabled genomic subtyping. We found that one-third of patients belonged to favorable genotypes (DUX4-rearranged, ETV6-RUNX1, TCF3-PBX1, PAX5 P80R, high hyperdiploidy). Of 63.4% of patients classified as high-risk based age and WBC count at diagnosis, almost 40% harbored subtype-defining genetic associations associated with favorable- or intermediate outcomes. These patients did significantly better than expected from their protocol-defined risk assignment. Our data in a large adult B-ALL cohort treated on a single international trial show the important prognostic value of genomic analysis and the potential benefit for many of the novel genotypes.

  Study EGAS00001004638

• Sporadic and endemic Burkitt lymphoma have frequent FOXO1 mutations but distinct hotspots in the AKT recognition motif

  FOXO1 has an oncogenic role in adult germinal center derived lymphomas, in which mutations, predominately within the AKT recognition motif, cause nuclear retention of FOXO1 resulting in increased cell proliferation. To determine the prevalence and distribution of FOXO1 mutations in pediatric Burkitt lymphoma (BL), we sequenced a large number of sporadic and endemic BL patient samples. We report a high frequency of FOXO1 mutations in both sporadic and endemic BL at diagnosis, occurring in 23/78 (29%) and 48/89 (54%) samples respectively, as well as 8/16 (50%) cases at relapse. Mutations of T24 were the most common in sporadic BL but were rare in endemic cases, in which mutations of residue S22, also within the AKT recognition motif, were the most frequent. FOXO1 mutations were almost always present in the major tumor cell clone but were not associated with outcome. Analysis of other recurrent mutations reported in BL revealed that FOXO1 mutations were associated with mutations of DDX3X and ARID1A, but not MYC, TCF3/ID3 or members of the phosphatidylinositol 3’ OH kinase (PI3K) signaling pathway. We further show common nuclear retention of the FOXO1 protein, irrespective of mutation status, suggesting alternative unknown mechanisms for maintaining FOXO1 transcriptional activity in BL. CRISPR/Cas9 knockout of FOXO1 in an endemic cell line produced a significant decrease in cell proliferation, supporting an oncogenic role for FOXO1 in endemic BL. Thus, FOXO1 is frequently mutated in both sporadic and endemic BL and may offer a potential therapeutic target for pediatric BL patients worldwide.

  Study EGAS00001003719

• Evolution of a FLT3-TKD mutated subclone at meningeal relapse in acute promyelocytic leukemia (H021)

  Here we report the case of an acute promyelocytic leukemia (APL) patient who, though negative for FLT3 mutations at diagnosis, developed isolated FLT3-TKD positive meningeal relapse, which, in retrospect, could be traced back to a minute bone marrow subclone present at first diagnosis. Initially, the 48-year old female diagnosed with high-risk APL had achieved complete molecular remission after standard treatment with all-trans retinoic acid (ATRA) and chemotherapy according to the AIDA protocol. 13 months after start of ATRA maintenance the patient suffered clinically overt meningeal relapse along with minute molecular traces of PML/RARA in the bone marrow. Following treatment with arsenic trioxide (ATO) and ATRA in combination with intrathecal cytarabine and methotrexate, the patient achieved a complete molecular remission in both cerebrospinal fluid (CSF) and bone marrow, which currently lasts for two years after completion of therapy. Whole exome sequencing and subsequent ultradeep targeted re-sequencing revealed a heterozygous FLT3-TKD mutation in CSF leukemic cells (p.D835Y, c.2503G>T, 1000/1961 reads (51%)), which was undetectable in the concurrent bone marrow sample. Interestingly, the FLT3-TKD mutated meningeal clone originated from a small bone marrow subclone present in a variant allele frequency of 0.4% (6/1553 reads) at initial diagnosis. This case highlights the concept of clonal evolution with a subclone harboring an additional mutation being selected as the “fittest” and leading to meningeal relapse. It also further supports earlier suggestions that FLT3 mutations may play a role for migration and clonal expansion in the CSF sanctuary site.H021

  Study EGAS00001001848

• Next-Generation Sequencing of AV Nodal Reentry Tachycardia patients

  Atrioventricular nodal reentry tachycardia (AVNRT) is the most common form of regular paroxysmal supraventricular tachycardia. This arrhythmia affects women twice as frequently as men, and is often diagnosed in patients below 40 years of age. Familial clustering, early onset of symptoms, and lack of structural anomaly indicate involvement of genetic factors in AVNRT pathophysiology. We hypothesized that AVNRT patients have a high prevalence of variants in genes that are highly expressed in the atrioventricular conduction axis of the heart and potentially involved in arrhythmic diseases. Next-generation sequencing of 67 genes was applied to the DNA profile of 298 AVNRT patients and 10 AVNRT family members using HaloPlex Target Enrichment System.In total, we identified 229 variants in 60 genes; 215 missenses, four frame shifts, four codon deletions, three missense and splice sites, two stop-gain variants, and one start-lost variant. Sixty-five of these were not present in the Exome Aggregation Consortium (ExAC) database. Furthermore, we report two AVNRT families with co-segregating variants. Seventy-five of 284 AVNRT patients (26.4%) and three family members to different AVNRT probands had one or more variants in genes affecting the sodium handling. Fifty-four out of 284 AVNRT patients (19.0%) had variants in genes affecting the calcium handling of the heart. We furthermore find a large proportion of variants in the HCN1-4 genes. We did not detect a significant enrichment of rare variants in the tested genes.This could be an indication that AVNRT might be an electrical arrhythmic disease with abnormal sodium and calcium handling.

  Study EGAS00001002745

• Dynamics of genome architecture and chromatin function during human B cell differentiation and neoplastic transformation

  Despite recent advances, the dynamics of genome architecture and chromatin function during human cell differentiation and its potential reorganization upon neoplastic transformation remains poorly characterized. Here, we integrate in situ Hi-C and nine additional omic layers to define and biologically characterize the dynamic changes in three-dimensional (3D) genome architecture across normal B cell differentiation and in neoplastic cells from different subtypes of chronic lymphocytic leukemia (CLL) and mantle cell lymphoma (MCL) patients. Beyond conventional active (A) and inactive (B) compartments, an integrative analysis of Hi-C data uncovers a highly-dynamic intermediate compartment enriched in poised and polycomb-repressed chromatin. During B cell development, we detect that 28% of the compartments change at defined maturation stages and mostly involve the intermediate compartment. The transition from naive to germinal center B cells is associated with widespread chromatin activation, which mostly reverts into the naive state upon further maturation of germinal center cells into memory B cells. The analysis of CLL and MCL neoplastic cells points both to entity and subtype-specific alterations in chromosome organization. Remarkably, we observe that large chromatin blocks containing key disease-specific genes alter their 3D genome organization. These include the inactivation of a 2Mb region containing the EBF1 gene in CLL and the activation of a 6.1Mb region containing the SOX11 gene in clinically aggressive MCL. This study indicates that 3D genome interactions are extensively modulated during normal B cell differentiation and that the genome of B cell neoplasias acquires a tumor-specific 3D genome architecture.

  Study EGAS00001004763

• Transcriptional profiling of ovarian tumours and cell lines

  Purpose: Endogenous retroviruses (ERVs) play a role in a variety of biological processes, including embryogenesis and cancer. DNA methyltransferase inhibitors (DNMTi)-induced ERV expression triggers interferon responses in ovarian cancer cells via the viral sensing machinery. Baseline expression of ERVs also occurs in cancer cells, though this process is poorly understood and previously unexplored in epithelial ovarian cancer (EOC). Here, the prognostic and immunomodulatory consequences of baseline ERV expression was assessed in EOC. Experimental Design: ERV expression was assessed using EOC transcriptional data from TCGA and from an independent cohort (HH), as well as from untreated or DNMTi-treated EOC cell lines. LASSO logistic regression defined an ERV expression score to predict patient prognosis. Immunohistochemistry (IHC) was conducted on the HH cohort. Combination of DNMTi treatment with γδ T cells was tested in vitro, using EOC cell lines and patient-derived tumor cells. Results: ERV expression was found to define clinically relevant subsets of EOC patients. An ERV prognostic score was successfully generated in TCGA and validated in the independent cohort. In EOC patients from this cohort, a high ERV score was associated with better survival (log-rank p=0.0009) and correlated with infiltration of CD8+PD1+ T cells (r=0.46, p=0.0001). In vitro, baseline ERV expression dictates the level of ERV induction in response to DNMTi. Manipulation of ERV expression by DNMTi resulted in improved EOC cell killing by cytotoxic immune cells. Conclusions: These findings uncover the potential for baseline ERV expression to robustly inform EOC patient prognosis, influence tumor immune infiltration and affect anti-tumor immunity.

  Study EGAS00001004814

• Colorectal cancer organoids expressing BRAF (fusion) genes

  Fusion genes may be oncogenic drivers and potential targets for personalized therapies in a variety of cancer types. The BRAF gene is frequently involved in oncogenic fusions, with fusion frequencies of 0.2-3% throughout different cancers. BRAF can be fused to a wide variety of genes, however, BRAF fusions rarely occur in the exact same gene configuration, making assessment of the fusion relevance and decision for drug treatment challenging. Here, we devised a colorectal cancer organoid-based platform to functionally characterize diverse BRAF fusions, containing various partner genes (AGAP3, DLG1 and TRIM24). We compared these BRAF fusions with respect to cellular behaviour, downstream signaling activation and response to targeted therapies. We found that 5’ partner choice of BRAF affects cellular localization and intracellular signaling capacities of the fusion genes. The DLG1-BRAF fusion gene showed distinct localization to the plasma membrane and exhibited increased levels of MAPK pathway activation under unperturbed conditions. Furthermore, the different BRAF fusions showed varying sensitivities to the targeted inhibition of ERK and BRAF, with the DLG1-BRAF fusion being the most sensitive. Additionally, RNA-sequencing identified distinct subsets of genes affected by the DLG1-BRAF fusion gene. Importantly, all fusion genes conveyed resistance to the clinically relevant EGFR/HER2/HER4-inhibitor afatinib, suggesting that BRAF fusions should be screened alongside other MAPK pathway alterations to identify mCRC patients amenable to cetuximab treatment. In summary, we developed a platform to efficiently assess molecular and cellular effects of fusion genes, and revealed that differential drug responses and distinct gene expression profiles can be induced by different BRAF fusions .

  Study EGAS00001003558

• MYO5B mutations in Pheochromocytoma/Paraganglioma promote cancer progression

  Identification of additional cancer-associated genes and secondary mutations driving the metastatic progression in pheochromocytoma and paraganglioma (PPGL) is important. We recently reported novel recurrent nonsynonymous mutations in the MYO5B gene in metastatic PPGL. Here, we explored the functional impact of these MYO5B mutations, and analyzed MYO5B expression in primary PPGL tumor cases in relation to mutation status.Immunohistochemistry and mRNA expression analysis in 30 primary PPGL tumors revealed an increased MYO5B expression in metastatic compared to non-metastatic tumor cases. In addition, subcellular localization of MYO5B protein was altered from cytoplasmatic to membranous in some metastatic tumors, and the strongest and most abnormal expression pattern was observed in a PGL case harboring a MYO5B:p.G1611S mutation. Mutation analysis of additional PPGL tumors revealed additional, and hence recurrent, mutations in the paralog MYO5A. Three missense mutations in MYO5B were functionally studied using site directed mutagenesis and stable transfection into human neuroblastoma cells (SK-N-AS) and embryonic kidney cells (HEK293 cell line). Functional analysis of the mutants verified a significant increased proliferation rate in all MYO5B mutated clones (p.L587P, p.G1611S, and p.R1641C). The two somatically derived mutations p.L587P and p.G1611S were also found to increase the migration rate, and all three mutants affected intracellular transport by increased transferrin uptake. Downstream expression analysis of MYO5B mutants demonstrated a transcriptional activation of genes involved in glucose metabolism, proliferation, homeostasis, and migration. Our study validates the functional role of novel MYO5B mutations in proliferation and migration, and suggest the MYO5-pathway to be involved the malignant progression of PPGL tumors.

  Study EGAS00001003991

• Phylogenetic analysis of treatment-naive metastases using whole exome and genome sequencing data

  Metastases are responsible for the majority of cancer related deaths and are often difficult to treat successfully. Even for metastases that occur subsequent to treatment, we do not know whether such relapsing metastases were originally heterogeneous or if the observed heterogeneity is a consequence of therapy. Given that treatment can influence evolutionary dynamics by selecting resistant clones, imposing bottlenecks on cancer cell populations, and even inducing novel somatic mutations, characterizing the prior standing variation among metastases remains an important goal to predict initial therapeutic response. To quantify the heterogeneity at clinical presentation of advanced disease, we surveyed the literature for patients in which at least two treatment-naïve metastases underwent genome/exome-wide sequencing. Across all cancer types surveyed, only 18 subjects were found to fulfill this requirement. Including previously unpublished data from two subjects, we analyzed data of 74 untreated metastases and inferred cancer phylogenies. Putative driver gene mutations were acquired at a 2-fold higher rate on the trunk of all metastases than on branches. Mutations in driver genes along the trunk were strongly enriched for predicted functional consequences. Using a stochastic mathematical model, we find that driver gene heterogeneity among metastases mostly occurs in slowly growing cancers for highly advantageous driver gene mutations. Seeding efficiency only weakly affects the heterogeneity among metastases. These results explain why functional driver gene heterogeneity is uncommon in advanced disease prior to treatment, thus providing optimism for future therapies that depend on mutations across metastases.

  Study EGAS00001002777

• A novel breast cancer cell line derived from a metastatic bone lesion of a breast cancer patient

  Purpose: We sought to generate and characterize a novel cell line from a breast cancer bone metastasis in order to better study the progression of the disease.Methods: The cell line, P7731, was derived from a metastatic bone lesion of a breast cancer patient and assessed for marker expression. P7731 was analyzed for DNA copy number variation, somatic mutations and gene expression and compared to the primary tumor. Results: P7731 cells are negative for estrogen receptor alpha (ERα), progesterone receptor (PR) and HER2 (triple-negative), strongly express vimentin (100% of cells positive) and also express cytokeratins 8/18 and 19 but at lower frequencies. Flow cytometry indicates P7731 cells are predominantly CD44+/CD49f+/EpCAM-, consistent with a primitive, mesenchymal-like phenotype. The cell line is tumorigenic in immunocompromised mice. Exome sequencing identified a total of 45 and 76 somatic mutations in the primary tumor and cell line respectively, of which 32 were identified in both samples and included mutations in known driver genes PIK3CA, TP53 and ARID1A. P7731 retains the DNA copy number alterations present in the matching primary tumor. Homozygous deletions detected in the cell line and in the primary tumor were found in regions containing three known (CDKN2A, CDKN2B and CDKN1B) and 23 putative tumor suppressor genes. Cell line specific gene amplification coupled with mRNA expression analysis revealed genes and pathways with potential pro-metastatic functions.Conclusion: This novel human breast cancer-bone metastasis cell line will be a useful model to study aspects of breast cancer biology, particularly metastasis-related changes from breast to bone.

  Study EGAS00001002840

• Colorectal adenomas, carcinomas and adjacent normal NKI-AvL TGO series NGS-ProToCol

  Cancer is caused by changes in the genome of cells, resulting in aberrant transcripts and proteins. In order to understand why normal cells become cancer cells and how we can identify and destroy them, the CTMM-NGS-ProToCol project aimed to understand prostate and colorectal cancer through DNA and RNA molecular profiling, to further improve diagnosis, prognosis and treatment.The state-of-the-art within the field of genome sequencing makes it possible to identify all the modifications in our genome that encode for RNA and of the transcripts themselves. The latest development in RNA sequencing is that if one does not select for polyA RNA, but first removes rRNA from the total RNA and then performs ‘random’ RNA sequencing, 30-60% of the reads are intronic. Since about 40% of the genome consists of genes (and only 1.6% is exonic), deep RNA sequencing will cover not only the exons, but also a large portion of the introns and therefore much of the genome. The power of the combined sequencing of genome and transcriptome is the full view of copy number aberrations (CNA), exon mutations and the consequence and relevance of the mutations at the RNA level in the form of differential gene expression, alternative splicing, alternative promoter usage, novel transcripts, fusion transcripts, read-through transcripts and mutations. For the colorectal (tumor) tissue dataset, we applied molecular profiling to 30 colorectal adenomas, 30 colorectal carcinomas and 18 normal adjacent colon tissues.

  Study EGAS00001002854

• Integrated clinical, whole genome, and transcriptome analysis of multisampled lethal metastatic prostate cancer

  We report the first combined analysis of whole genome sequence, detailed clinical history, and transcriptome sequence of multiple prostate cancer metastases in a single patient (A21). Whole genome and transcriptome sequence was obtained from 9 anatomically separate metastases, and targeted DNA sequencing was performed in cancerous and noncancerous foci within the primary tumor specimen removed 5 years prior to death. Transcriptome analysis revealed increased expression of AR-regulated genes in liver metastases that harbored an AR p.L702H mutation, suggesting a dominant effect by the mutation despite being present in only 1 of an estimated 16 copies per cell. The metastases harbored several alterations to the PI3K/AKT pathway, including a clonal truncal mutation in PIK3CG and present in all metastatic sites studied. The list of truncal genomic alterations shared by all metastases included homozygous deletion of TP53, hemizygous deletion of RB1 and CHD1, and amplification of FGFR1. If the patient were treated today given this knowledge, use of second-generation androgen-directed therapies, cessation of glucocorticoid administration, and therapeutic inhibition of the PI3K/AKT pathway or FGFR1 receptor could provide personalized benefit. Three previously unreported truncal clonal missense mutations (ABCC4 p.R891L, ALDH9A1 p.W89R, and ASNA1 p.P75R) were expressed at the RNA level and assessed as druggable. The truncal status of mutations is critical for actionability, and can only be determined through analysis of multiple sites of metastasis. Our findings suggest that a large set of deeply analyzed cases could serve as powerful guide to more effective prostate cancer basic science and personalized cancer medicine clinical trials.

  Study EGAS00001001659

• Clonal fitness inferred from timeseries modeling of single cell cancer genomes

  Progress in defining genomic fitness landscapes in cancer by copy number alterations (CNA) has been impeded by lack of single cell and timeseries sampling. We generated 42,000 single cell whole genomes (scWGS) from breast epithelium and primary triple negative breast cancer (TNBC) patient-derived xenografts (PDX) collected during multi-year time series. Using a Wright-Fisher population genetics model, we inferred reproducible CNA-defined clonal fitness dynamics induced by TP53 mutation and cisplatin chemotherapy, with accurate forecasting of experimentally enforced clonal competition dynamics. Drug treatment in three long-term serially passaged TNBC PDX resulted in cisplatin resistant clones that had shown low fitness in the untreated setting. By contrast high fitness clones from treatment naive controls were eradicated reflecting an inversion of the fitness landscape. Upon drug release selective pressure dynamics were reversed indicating a fitness cost of treatment resistance. Taken together, our findings reveal clonal fitness dynamics linked to CNA and therapeutic resistance in polyclonal tumours.

  Study EGAS00001004448

• Versatile workflow for cell-type resolved transcriptional and epigenetic profiles from cryopreserved human lung

  The complexity of the lung microenvironment and changes in cellular composition during disease progression make it exceptionally hard to understand molecular mechanisms leading to the development of chronic lung diseases. Although recent advances in cell-type resolved and single-cell sequencing approaches hold great promise for studying complex diseases, their implementation relies on local access to fresh tissue, as traditional tissue storage methods do not allow viable cell isolation. To overcome these hurdles, we developed a novel, versatile workflow that allows long-term storage of human lung tissue with high cell viability, permits thorough sample quality check before cell isolation, and is compatible with sequencing-based profiling, including single-cell approaches. We demonstrate that cryopreservation is suitable for the isolation of multiple cell types from human lung and is applicable to both healthy and diseased tissue, including COPD and tumour samples. Basal cells isolated from cryopreserved airways retain the ability to differentiate, indicating that cellular identity is not altered by cryopreservation. Importantly, using RNA sequencing and EPIC Array, we show that genome-wide gene expression and DNA methylation signatures are preserved upon cryopreservation, emphasizing the suitability of our workflow for -omics profiling of human lung cells. In addition, we obtained high-quality single-cell RNA sequencing data of cells isolated from the cryopreserved human lung, demonstrating that cryopreservation empowers single-cell approaches. Overall, thanks to its simplicity, our cryopreservation workflow is well-suited for prospective tissue collection by academic collaborators and biobanks, opening worldwide access to human tissue.

  Study EGAS00001004477

• Identification of Therapeutic Targets in Rhabdomyosarcoma Through Integrated Genomic, Epigenomic, and Proteomic Analyses

  Personalized cancer therapy based on the somatic mutations identified in patient tumors is becoming an increasingly emphasized approach to improve outcomes of patients with cancer. There are also examples of therapeutic vulnerabilities in cancer that result from changes in gene expression that are a direct or indirect result of tumor specific epigenetic perturbations. These genomic and epigenomic changes are ultimately manifest in the tumor proteome and phosphoproteome. In this study, we integrated genomic, epigenomic and proteomic data for rhabdomyosarcoma (RMS) to identify therapeutic vulnerabilities. RMS was selected for this analysis because RAS pathway mutations in rhabdomyosarcoma (RMS) are the most common potentially actionable lesions across pediatric solid tumors. The epigenomic data was useful for identifying deregulated developmental pathways in RMS including the WNT, HH, BMP, adenyl cyclase, p38/MAPK and PI3K pathways. Perturbations in those 6 myogenic signal transduction pathways were also evident in the proteome and phosphoteome data. In addition, the proteomic/phosphoproteomic data revealed that the cell cycle checkpoint, unfolded protein response and RB/E2F pathways were deregulated in RMS relative to normal muscle. Recent success targeting CDK4/6 and MEK in adult cancers with RAS mutations led us to test the value of these targets in RMS tumors. We also targeted the unfolded protein response pathway and cell cycle checkpoint pathway using molecular targeted therapeutics in orthotopic patient derived xenografts. Taken together, these data demonstrate the value of integrating epigenomic and proteomic data to identify tumor vulnerabilities that extend beyond somatic mutations identified in the genome.

  Study EGAS00001002967

• Detection and characterization of lung cancer using cell-free DNA fragmentomes

  Lung cancer remains the leading cause of cancer death world-wide, largely due to its late diagnosis. Non-invasive approaches for assessment of cell-free DNA (cfDNA) provide an opportunity for detection and intervention that may have broader accessibility than current imaging approaches. Using a machine learning model for detecting tumor-derived cfDNA through genome-wide analyses of cfDNA fragmentation, we examined a prospective study of 365 individuals at risk for lung cancer (Lung Cancer Diagnostic Study, LUCAS), including 129 individuals ultimately diagnosed with lung cancer and 236 individuals determined to not have lung cancer. We externally validated the cancer detection model using an independent cohort of 385 non-cancer individuals and 46 predominantly early stage lung cancer patients. Combining fragmentation features with clinical risk factors and CEA levels followed by CT imaging detected 94% of patients with cancer across stages and subtypes, including 91% of stage I/II and 96% of stage III/IV, at 80% specificity. Genome-wide fragmentation profiles across ~13,000 ASCL1 transcription factor binding sites could be used to distinguish individuals with small cell lung cancer (SCLC) from those with non-small cell lung cancer (NSCLC) with high accuracy (AUC=0.98). Among individuals with lung cancer, a higher cfDNA fragmentation score was associated with tumor size and invasion, and represented an independent prognostic indicator of survival. These studies provide a facile approach for non-invasive detection of lung cancer and clinical management of this disease.

  Study EGAS00001005340

• Human genomic and phenotypic synthetic data for the study of rare diseases

  The purpose of this dataset is to facilitate development of technical implementations for rare disease data integration, analysis, discovery, and federated access. This synthetic dataset includes clinical and genomic data from 6 rare disease cases. It consists of 18 whole genomes (6 index cases with their parents) which have genetic background based on public human data sequenced in the context of the Illumina Platinum initiative (Eberle, MA et al. (2017)) and made available by the HapMap project (https://www.genome.gov/10001688/international-hapmap-project). In each of the cases, real causative variants correlating with the phenotypic data provided were spiked-in. The cases included in this synthetic dataset correspond to the following type of disorders: CASE 1- Congenital myasthenic syndrome (Autosomal Dominant -de novo variant) CASE 2- Macular dystrophy (Autosomal Dominant) CASE 3- Muscular dystrophy (Autosomal Recessive-compound heterozygous variants) CASE 4- Mitochondrial disorder (Autosomal Recessive-consanguineous case - homozygous variant) CASE 5- Breast cancer (Autosomal Dominant) CASE 6- Similar as case 1 for patient matchmaking tests: Congenital myasthenic syndrome (Autosomal Dominant-de novo variant) For each case you will be able to download the following data: clinical information (phenopackets per individual and pedigree per family), raw genomic data (FASTQ and BAMs) and processed genomic data (vcfs). When using the data, the following should be acknowledged: the RD-Connect GPAP (https://platform.rd-connect.eu/), EC H2020 project EJP-RD (grant # 825575), EC H2020 project B1MG (grant # 951724) and Generalitat de Catalunya VEIS project (grant # 001-P-001647).

  Study EGAS00001005702

• Open_Targets_020_Epigenomes_of_Cell_Lines

  There is currently a drive to establish cell based assay systems of greater human biological and disease relevance through the use of well characterised transformed cell lines, primary cells and complex cellular models (e.g. co-culture, 3D models). However, although the field is gaining valuable experience in running more non-standard & complex cell assays for target validation and compound pharmacology studies, there is the lack of a systematic approach to determine if this expansion in cell assay models is reflected in increased human biological and disease relevance. The increasing wealth of publically available transcriptomic, and epigenome (ENCODE and Epigenome Roadmap) data represents an ideal reference mechanism for determining the relationship between cell types used for target & compound studies to primary human cells and tissues from both healthy volunteers & patients. The CTTV020 epigenomes of cell line project aims to generate epigenetic and transcriptomic profiles of cell lines and compare these with existing and newly generated reference data sets from human tissue and cell types. The aim is to identify assay systems which will provide greater confidence in translating target biology and compound pharmacology to patients. Multiple cell types commonly used within research have been grouped according to biology. Examples include erythroid, lung epithelial, hepatocyte cell types and immortalised models of monocyte / macrophage biology. 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 EGAS00001003137

• Open_Targets_020_Epigenomes_of_Cell_Lines

  There is currently a drive to establish cell based assay systems of greater human biological and disease relevance through the use of well characterised transformed cell lines, primary cells and complex cellular models (e.g. co-culture, 3D models). However, although the field is gaining valuable experience in running more non-standard & complex cell assays for target validation and compound pharmacology studies, there is the lack of a systematic approach to determine if this expansion in cell assay models is reflected in increased human biological and disease relevance. The increasing wealth of publically available transcriptomic, and epigenome (ENCODE and Epigenome Roadmap) data represents an ideal reference mechanism for determining the relationship between cell types used for target & compound studies to primary human cells and tissues from both healthy volunteers & patients. The CTTV020 epigenomes of cell line project aims to generate epigenetic and transcriptomic profiles of cell lines and compare these with existing and newly generated reference data sets from human tissue and cell types. The aim is to identify assay systems which will provide greater confidence in translating target biology and compound pharmacology to patients. Multiple cell types commonly used within research have been grouped according to biology. Examples include erythroid, lung epithelial, hepatocyte cell types and immortalised models of monocyte / macrophage biology. 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 EGAS00001003139

• Healthy_ageing_thymus

  Single cell RNA sequencing of immune and non-immune cells from healthy ageing thymus Thymus is a primary lymphoid organ that creates an environment for the T lymphocyte precursors differentiation into the naive T cells, which involves TCR gene rearrangement, negative and positive selection. Despite thymus critical function for the recognition of the pathogens it starts to atrophy very early in life. Thymic Epithelial Space (TES), where thymocyte education occurs, starts to decline rapidly just after birth, while with the onset of puberty thymus is also progressively replaced with adipose tissue. This leads to a reduced naive T cell output,which decreases an organism's ability to recognize pathogens. In our lab, we have already collected and characterized cell type composition and gene regulatory networks of the fetal tand adult human thymus (Park et al, Science 2020). In the current project we plan to expand the number of paediatric and adult samples in the study with the aim to understand how the process of thymic involution happens. We will use single-cell transcriptomics to zoom on changes that occur in immune and non-immune cell types in different phases of thymic involution with the hope to understand which gene expression changes are causal for the thymic involution process and how naïve T cell production can be increased in elderly. 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 EGAS00001004311

• Novel optineurin frameshift insertion causing familial frontotemporal dementia and parkinsonism without amyotrophic lateral sclerosis

  Frontotemporal Dementia (FTD) is a common cause of Young Onset Dementia and has diverse clinical manifestations involving behaviour, executive function, language and motor function, including parkinsonism. Up to 50% of FTD patients report a positive family history, supporting a strong genetic basis, particularly in cases with both FTD and amyotrophic lateral sclerosis (FTD-ALS). Mutations in three genes are associated with the majority of familial FTD (fFTD) cases - microtubule associated protein tau gene (MAPT), progranulin gene (PGRN), and hexanucleotide repeat expansions in chromosome 9 open reading frame 72 (C9orf72), whilst mutations in other genes such as optineurin (OPTN) have rarely been reported. Mutations in OPTN have been reported mostly in familial and sporadic cases of ALS, or in rare cases of FTD-ALS, but not in association with pure or predominant FTD and/or parkinsonian phenotype. Here, we report for the first time, a family from the Philippines with 4 affected members harbouring the same OPTN frameshift insertion, presenting with FTD-related phenotypes, including one sibling with predominant parkinsonism resembling corticobasal syndrome. Notably, none of the affected members showed any evidence of motor neuron disease or ALS at the time of writing, both clinically and on electrophysiological testing, expanding the phenotypic spectrum of OPTN mutations. Close follow-up of mutation carriers for the development of new clinical features and wider investigation of additional family members with further genetic analyses will be conducted to investigate the possibility of other genetic modifiers in this family which could explain phenotypic heterogeneity.

  Study EGAS00001005220

• Exome-wide mutation analysis of cell-free DNA to simultaneously monitor the full spectrum of cancer treatment outcomes

  Purpose: Cell-free DNA (cfDNA) offers a non-invasive approach to monitor cancer. Here we develop a method using whole-exome sequencing (WES) of cfDNA for simultaneously monitoring the full spectrum of cancer treatment outcomes, including MRD, recurrence, evolution, and second primary cancers. Experimental Design: Three simulation datasets were generated from 26 cancer patients to benchmark the detection performance of MRD/recurrence and second primary cancers. For further validation, cfDNA samples (n=76) from cancer patients (n=35) with six different cancer types were used for performance validation during various treatments.Results: We present a cfDNA-based cancer monitoring method, named cfTrack. Taking advantage of the broad genome coverage of WES data, cfTrack can sensitively detect MRD and cancer recurrence by integrating signals across known clonal tumor mutations of a patient. In addition, cfTrack detects tumor evolution and second primary cancers by de novo identifying emerging tumor mutations. A series of machine learning and statistical denoising techniques are applied to enhance the detection power. On the simulation data, cfTrack achieved an average AUC of 99% on the validation dataset and 100% on the independent dataset in detecting recurrence in samples with tumor fraction ≥0.05%. In addition, cfTrack yielded an average AUC of 88% in detecting second primary cancers in samples with tumor fraction ≥0.2%. On real data, cfTrack accurately monitors tumor evolution during treatment, which cannot be accomplished by previous methods.Conclusion: Our results demonstrated that cfTrack can sensitively and specifically monitor the full spectrum of cancer treatment outcomes using exome-wide mutation analysis of cfDNA.

  Study EGAS00001005906

• Transcriptome Sequencing PPGL

  RNA library preparation of samples with RNA Integrity Number (RIN) > 5.5 (median 7.1; range 5.6-8.9) was performed as described in the TruSeq Stranded mRNA Library Prep Kit (Illumina, RS-122-2101). Briefly, 0.7-1 µg total RNA PolyA+ fraction was purified and randomly fragmented, converted to double stranded cDNA and processed through subsequent enzymatic treatments of end-repair, dA-tailing, and ligation to adapters. Adapter-ligated library was completed by PCR with Illumina PE primers. The resulting purified cDNA library was applied to an Illumina flow cell for cluster generation and sequenced on an Illumina HiSeq2500 on a 51bp single-read format following manufacturer's protocols. cDNA libraries from FFPE tumors and low integrity RNA’s (RIN<5.5) from FF tumors (median 2.4; range 1-5.5) but % of fragments > 200 nucleotides higher than 20% (median 57.5; range 20-96) were prepared using QuantSeq 3’ mRNA-Seq Library Prep Kit FWD for Illumina (Lexogen, 015) with a UMI Second Strand Synthesis Module for QuantSeq FWD (Lexogen, 081), vendor’s protocol for low input/ low quality/ FFPE RNA using 200-500 ng of total RNA, according to availability, was used as starting material. PCR Add-on Kit for Illumina (Lexogen, 020) was used to adjust the library amplification number of cycles. Libraries were applied to an Illumina flow cell for cluster generation and sequenced on NovaSeq6000. Image analysis, per-cycle basecalling and quality score assignment was performed with Illumina Real Time Analysis software. Demultiplexing of BCL files to FASTQ format was performed with the bcl2fastq Software (Illumina).

  Study EGAS00001006044

• A higher ctDNA fraction decreases survival in regorafenib-treated metastatic colorectal cancer patients

  The predictive effect of circulating tumor DNA (ctDNA) in colorectal cancer (CRC) treatment is still highly discussed. The primary objective of this study was to investigate a possible prognostic/predictive value of ctDNA under regorafenib treatment. This prospective multicenter translational biomarker phase II pilot study enrolled 30 metastatic CRC patients (67% men, 33% women) treated with regorafenib. ctDNA was assessed in plasma before treatment start and at defined time points during administration. Measurement of tumor fraction as well as mutation and copy number analysis of CRC driver genes were performed by next generation sequencing approaches. Multivariate analyses for survival and treatment efficacy were adjusted to age, gender and ECOG. Disease control rate was 30%. Median tumor fraction at baseline was 18.5% (0-49.9). Mutations in CRC driver genes or genes involved in angiogenesis were identified in 25 patients (83.3%). KRAS mutations were detected in 13 out of 14 KRAS positive tumors, in three patients without KRAS mutation in the respective tumor acquired mutations as a consequence of prior anti-EGFR treatment were detected. In a subset of patients, novel occurring mutations or focal amplifications were detected. A tumor fraction of 5% and higher at baseline was significantly associated with a decreased OS (p=0.022; HHR 3.110 (95%CI: 1.2-8.2). ctDNA is detectable in a high proportion of mCRC patients. Higher ctDNA levels are associated with survival among regorafenib treatment. Moreover, our data highlight the benefit of a combined evaluation of mutations and somatic copy number alterations (SCNA) in advanced cancer patients.

  Study EGAS00001004491

• Butyrate producers as potential next-generation probiotics: safety assessment of the administration of Butyricicoccus pullicaecorum to healthy volunteers

  Advances in gut microbiota research have triggered interest in developing colon butyrate producers as niche-specific next-generation probiotics, targeted at increasing colon butyrate production and countering disease-associated microbiota alterations. Crucial steps in the development of next-generation probiotics are the design of formulations with a reasonable shelf life as well as the safety demonstration of an intervention in healthy volunteers. One such potential next-generation butyrate-producing probiotic is Butyricicoccus pullicaecorum 25-3T, with demonstrated safety in in vitro as well as in animal models. Here, we examined the strain's safety, tolerability, and impact on microbiota composition and metabolic activity in healthy volunteers in a randomized, double blind, placebo-controlled, crossover study in 30 healthy volunteers. Study design consisted of two four-week intervention periods (108 CFU B. pullicaecorum [treatment] or maltodextrin [placebo] per day) with a three-week washout in between. We assessed adverse events, blood parameters (primary endpoints), and fecal microbiota composition and metabolite profiles (secondary endpoints).The number of reported adverse events during the B. pullicaecorum treatment was similar to placebo intervention, as were observed changes in blood chemistry parameters, bowel habits, and fecal calprotectin concentrations. Administration of the strain did not induce any disruptive effect in microbiota composition or metabolic activity. In this first human intervention trial with a butyrate-producing Clostridium cluster IV isolate, we demonstrated B. pullicaecorum 25-3T administration to be both safe and well tolerated by healthy participants. This safety study paves the way for the further development of the strain as a next-generation probiotic.

  Study EGAS00001003276

• Transcriptomic intra-tumor heterogeneity of colorectal cancer varies depending on tumor location within the colorectum

  Background; Intra-tumor heterogeneity (ITH) of colorectal cancer (CRC) complicates molecular tumor classification, such as transcriptional subtyping. Differences in cellular states, biopsy cell composition, and tumor microenvironment may all lead to ITH. Here we analyze ITH at the transcriptomic and proteomic levels to ascertain whether subtype discordance between multiregional biopsies reflects relevant biological ITH or lack of classifier robustness. Further, we study the impact of tumor location on ITH.MethodsMultiregional biopsies from stage II and III CRC tumors were analyzed by RNA sequencing (41 biopsies, 14 tumors) and multiplex immune protein analysis (89 biopsies, 29 tumors). CRC subtyping were performed using consensus molecular subtypes (CMS), CRC intrinsic subtypes (CRIS), and TUMOR types. ITH-scores and network maps were defined to determine the origin of heterogeneity. A validation cohort was used with one biopsy per tumor (162 tumors).ResultsOverall, inter-tumor transcriptional variation exceeded ITH, and subtyping calls were frequently concordant between multiregional biopsies. Still, some tumors had high ITH and were classified discordantly. Proximal tumors showed ITH related to differences in the microenvironment, while ITH within distal tumors were cancer-cell related. This subtype discordancy reflected actual molecular ITH within the tumors. The relevance of the subtypes was reflected at protein level where several inflammation markers were significantly increased in immune related transcriptional subtypes, which was verified in an independent cohort (Wilcoxon rank sum test; p<0.05).ConclusionOur transcriptomic and proteomic analyses show that the tumor location along the colorectum influence the ITH of CRC, which again influence the concordance of subtyping.

  Study EGAS00001004668

• Exome sequencing

  Human telomere biology disorders (TBD), a heterogeneous group of disorders characterized by telomere attrition and premature aging phenotypes, are caused by inherited loss-of-function mutations in telomere-associated genes. Here, we identify three germline heterozygous missense variants in RPA1 gene in four unrelated probands presenting with short telomeres and varying clinical features of TBD including bone marrow failure, myelodysplastic syndrome, T- and B-cell lymphopenia, pulmonary fibrosis, and skin manifestations. All variants cluster to DNA binding domain A of RPA1 protein. RPA1 is an essential single-strand DNA-binding protein essential for DNA replication and repair. RPA1 has been implicated in telomere maintenance but the mechanism remains elusive. Thus far, no genetic disorder had been linked to RPA1. Using FRET assays, we showed that RPA1E240K and RPA1V227A exhibit increased binding to single-strand and telomeric DNA, implying a gain in DNA-binding function while RPA1T270A has binding properties similar to wild type. To study the mutational effect, RPA1E240K, with the highest binding capacity, was knocked-in iPSCs using CRISPR/Cas9, which showed severe telomere shortening and impaired hematopoietic differentiation. Furthermore, in the RPA1E240K patient, stable blood counts over two decades due coincided with clonal rescue hematopoiesis with an RNA-degrading truncating RPA1 mutation and a uniparental isodisomy 17p with loss of RPA1E240K allele. Using single cell DNA sequencing, the two somatic genetic rescue events were proven to be independently acquired in hematopoietic stem cells. In summary, we describe the first human disease caused by germline RPA1 variants in individuals with TBD.

  Study EGAS00001005761

• Mutation-specific non-canonical pathway of PTEN as a distinct therapeutic target for glioblastoma

  PTEN is one of most frequently altered tumor suppressor genes in malignant tumors. The dominant negative effect of PTEN alteration suggests that aberrant function of PTEN mutation might be more disastrous than deletion, the most frequent genomic event in glioblastoma (GBM). This study aimed to understand the functional properties of various PTEN missense mutations and to investigate their clinical relevance. The genomic landscape of PTEN alteration was analyzed using the Samsung Medical Center GBM cohort and validated via The Cancer Genome Atlas dataset. Several hot-spot mutations were identified, and their subcellular distributions and phenotypes were evaluated. We established a library of cancer cell lines that over-express these mutant proteins using the U87MG and patient-derived cell models lacking functional PTEN. PTEN mutations were categorized into two major subsets: missense mutations in the phosphatase domain and truncal mutations in the C2 domain. We determined the subcellular compartmentalization of four mutant proteins (H93Y, C124S, R130Q, and R173C) from the former group and found that they had distinct localizations; those associated with invasive phenotypes (‘edge mutations’) localized to the cell periphery, while the R173C mutant localized to the nucleus. Invasive phenotypes derived from edge substitutions were unaffected by an anti-PI3K/Akt agent but were disrupted by microtubule inhibitors. PTEN mutations exhibit distinct functional properties with regard to their subcellular localization. Further, some missense mutations (‘edge mutations’) in the phosphatase domain caused enhanced invasiveness associated with dysfunctional cytoskeletal assembly, thus suggesting it to be a potent therapeutic target.

  Study EGAS00001004753

• RaScALL: Rapid screening of RNA-seq in acute lymphoblastic leukaemia

  RNA-sequencing (RNA-seq) efforts in acute lymphoblastic leukaemia (ALL) have identified numerous prognostically significant genomic alterations which can guide diagnostic risk stratification and treatment choices when detected early. However, a full RNA-seq Bioinformatics workflow is time-consuming and costly in a clinical setting where rapid detection and accurate reporting of clinically relevant alterations are essential. To accelerate the identification of ALL-specific alterations (including gene fusions, single nucleotide variants and focal gene deletions), we developed the rapid screening tool RaScALL, capable of identifying more than 100 prognostically significant lesions directly from raw sequencing reads. RaScALL uses the k-mer based targeted detection tool km and known ALL variant information to achieve a high degree of accuracy for reporting subtype defining genomic alterations compared to standard alignment-based pipelines. Gene fusions, including difficult to detect fusions involving EPOR and DUX4, were accurately identified in 98% (164 samples) of reported cases in a 180-patient Australian study cohort and 95% (n=63) of samples in a North American validation cohort. Pathogenic sequence variants were correctly identified in 75% of tested samples, including all cases involving subtype defining variants PAX5 p.P80R (n=12) and IKZF1 p.N159Y (n=4). Accurate detection of intragenic IKZF1 deletions resulting in aberrant transcript isoforms was also detectable with 98% accuracy. Importantly, the median analysis time for detection of all targeted alterations averaged 22 minutes per sample, significantly shorter than standard alignment-based approaches, ensuring accerelated risk-stratification and therapeutic triage.

  Study EGAS00001006460

• Integrative genomic analysis identifies key pathogenic mechanisms in primary mediastinal large B-cell lymphoma

  Primary mediastinal large B-cell lymphoma (PMBL) represents a clinically and pathologically distinct subtype of large B-cell lymphomas. Furthermore, molecular studies, including global gene expression profiling, have provided evidence that PMBL is more closely related to classical Hodgkin lymphoma (cHL). Although targeted sequencing studies have revealed a number of mutations involved in PMBL pathogenesis, a comprehensive description of disease-associated genetic alterations and perturbed pathways is still lacking. Here, we performed whole-exome sequencing of 95 PMBL tumors to inform on oncogenic driver genes and recurrent copy number alterations. The integration of somatic gene mutations with gene expression signatures provides further insights into genotype-phenotype interrelation in PMBL. We identified highly recurrent oncogenic mutations in the JAK-STAT and NF-kB pathways, and provide additional evidence of the importance of immune evasion in PMBL (CIITA, CD58, B2M, CD274, PDCD1LG2). Our analyses highlight the IRF-pathway as a putative novel hallmark with frequent alterations in multiple pathway members (IRF2BP2, IRF4, IRF8). In addition, our integrative analysis illustrates the importance of JAK1, RELB and EP300 mutations driving oncogenic signaling. The identified driver genes were significantly more frequently mutated in PMBL as compared to diffuse large B-cell lymphoma, whereas only a limited number of genes were significantly different between PMBL and cHL, emphasizing the close relationship between these entities. Our study, performed on a large cohort of PMBL, highlights the importance of distinctive genetic alterations for disease taxonomy with relevance for diagnostic work-up and therapeutic decision-making.

  Study EGAS00001003746