23014 results for "ega"

in 66.62 milliseconds.

• 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

• 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

• 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

• 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

• 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