Search Results - EGA European Genome-Phenome Archive

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

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

Study EGAS00001002270

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

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

Study EGAS00001002294

Botswana 15 autosomal unlinked microsatellites

Population structure was investigated in 990 Botswana individuals according to ethno-linguistics, Bantu and Khoisan, and geography (the nine administrative districts) using the AmpFLSTR Identifiler autosomal microsatellite markers. Genetic diversity and forensic parameters were calculated for the overall population, and according to ethno-linguistics and geography. The overall combined power of exclusion (CPE) was 0.9999965412 and the combined match probability 6,28 x 10-19. CPE was highest for the Khoisan Tuu ethnolinguistic group and the Northeast District at 0.9999582029 and 0.9999922652 respectively. CMP ranged from 6.28 x 10-19 (Khoisan Tuu) to 1,02 x 10-18 (Northwest district). Using pairwise genetic distances (FST), analysis of molecular variance (AMOVA), factorial correspondence analysis (FCA), and the unsupervised Bayesian clustering method found in STRUCTURE and landscape genetics software TESS, ethno-linguistics were found to have a greater influence on population structure than geography. FCA showed clustering between Bantu and Khoisan, and within the Bantu. This Bantu sub-structuring was not seen with STRUCTURE and TESS, which detected clustering only between Bantu and Khoisan. The patterns of population structure revealed highlight the need for regional reference databases that include ethno-linguistic and geographic location information. These markers have important potential for bio-anthropological studies as well as for forensic applications.

Study EGAS00001002380

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

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

Study EGAS00001002432

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

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

Study EGAS00001002454

Direct Detection of Early-stage Cancers Using Circulating Tumor DNA

Early detection and intervention are likely to be the most effective means for reducing morbidity and mortality of human cancer. However, development of methods for noninvasive detection of early-stage tumors has remained a challenge. We have developed an approach called targeted error correction sequencing (TEC-Seq) that allows ultrasensitive direct evaluation of sequence changes in circulating cell-free DNA using massively parallel sequencing. We have used this approach to examine 58 cancer-related genes encompassing 81 kb. Analysis of plasma from 44 healthy individuals identified genomic changes related to clonal hematopoiesis in 16% of asymptomatic individuals but no alterations in driver genes related to solid cancers. Evaluation of 200 patients with colorectal, breast, lung, or ovarian cancer detected somatic mutations in the plasma of 71, 59, 59, and 68%, respectively, of patients with stage I or II disease. Analyses of mutations in the circulation revealed high concordance with alterations in the tumors of these patients. In patients with resectable colorectal cancers, higher amounts of preoperative circulating tumor DNA were associated with disease recurrence and decreased overall survival. These analyses provide a broadly applicable approach for noninvasive detection of early-stage tumors that may be useful for screening and management of patients with cancer.

Study EGAS00001002577

Inhibiton of the GABPB1L-containing GABP tetramer is sufficient to reverse replicative immortality in TERT promoter mutant glioblastoma cells.

Over 50 types of cancer acquire TERT promoter mutations. These single point mutations reactivate telomerase, allowing for indefinite maintenance of telomere length and enabling cellular immortalization. The transcription factor binding site created by the point mutations specifically recruit the ETS factor GABP, a multimeric transcription factor composed of the GABPα and GABPβ subunits. GABP can form two functionally independent transcription factor species – a dimer or a tetramer – depending on which of the structurally distinct GABPβ isoforms is incorporated into the complex. In this study, we show that genetic disruption of GABPβ1L, a tetramer forming isoform of GABPβ that is dispensable in normal development, results in TERT silencing in a TERT promoter mutation dependent manner. Failure to activate TERT expression by GABPβ1L culminates in telomere dysfunction, DNA damage, and mitotic cell death exclusively in TERT promoter mutant cells. Furthermore, exogenous expression of TERT is sufficient to prevent telomere degradation and loss of cell viability in GABPβ1L-reduced lines bearing the mutant TERT promoter. Orthotopic injection of tetramer-deficient mutant TERT promoter GBM cells rendered lower tumor burden and prolonged the overall survival of the tumor-bearing mice. These results highlight the potentially widespread role of GABPB1L in enabling replicative immortality of TERT promoter cancers.

Study EGAS00001002582

WGS of 32 paired SRCC samples

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

Study EGAS00001002668

UK10K_RARE_NEUROMUSCULAR

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 neuromuscular disorder samples are part of the “rare disease” group, and will undergo exome sequencing. For further information with regard to this cohort please contact Francesco Muntoni (f.muntoni@ucl.ac.uk).

Study EGAS00001000101

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

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

Study EGAS00001002728

UK10K_RARE_FIND

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 Raymond samples will be part of the “rare disease” group, and will undergo exome sequencing. For further information with regard to this cohort please contact Lucy Raymond (flr24@cam.ac.uk).

Study EGAS00001000128

UK10K_RARE_CILIOPATHIES

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 Ciliopathies samples will be part of the rare disease group, and will undergo exome sequencing. For further information with regard to this cohort please contact Phil Beales (p.beales@ich.ucl.ac.uk).

Study EGAS00001000126

UK10K_RARE_SIR

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 SIR (Severe Insulin Resistance) samples are part of the “rare disease” group, and will undergo exome sequencing. For further information with regard to this cohort please contact Robert Semple (rks16@cam.ac.uk).

Study EGAS00001000130

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

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

Study EGAS00001004290

Stratifying and Targeting Pediatric Medulloblastoma through Genomics

In this project, genomic analyses of pediatric medulloblastoma samples, obtained through the international medulloblastoma consortium, will be performed. RNA and miRNA expression profiles of 1000 samples, representing all four subgroups (Wnt, Shh, Group C, and D), will be studied to identify novel subtypes within each subgroup. The resulting subtype-specific expression profiles will support the development of reliable and robust biomarkers to more accurately and reliably classify medulloblastomas for treatment in clinical trials. For that purpose, two assays will be developed: an antibody-based immunohistochemical assay and an orthogonal nucleic acid-based hybridization assay. Additional genomic DNA analysis of the 300 high risk subgroup cases will support the discovery of subgroup specific somatic mutations in order to inform current clinical trials of targeted therapies, and to identify genes and pathways already targeted in other diseases. Such therapies could be rapidly transitioned to Phase II trials in medulloblastoma. Furthermore, the discovery of somatic mutations could be used for developing as well as validating specific biomarkers. The project team will also try to identify risk factors that predispose children to this type of cancer.

Study EGAS00001000273

Clonal_haematopoiesis_in_patients_with_AAA

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

Study EGAS00001002873

UK10K RARE CHD

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 CHD (congenital heart disease) samples will be part of the “rare disease” group, and will undergo exome sequencing.For further information with regard to this cohort please contact Shoumo Bhattacharya (shoumo@me.com).

Study EGAS00001000125

Genomic Alterations in Gingivo-buccal Cancer: ICGC-India Project_YR02

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

Study EGAS00001001028

P647_Targeted_resequencing_project

In collaboration with Dr Robert Semple we have identified a family harbouring an autosomal dominant variant, which leads to severe insulin resistance (SIR), short stature and facial dysmorphism. This family is unique within the SIR cohort in having normal lipid profiles, preserved adiponectin and normal INSR expression and phosphorylation. DNA is available for 7 affected and 7 unaffected family members across 3 generations. All 14 samples have been genotyped using microsatellites and the Affymetrix 6.0 SNP chip. Linkage analysis identified an 18.8Mb haplotype on chromosome 19 as a possible location of the causative variant. However, Exome sequencing of 3 affected and 1 unaffected family members has not identified the causative variant suggesting the possibility of an intronic or intergenic variant in this region or elsewhere in the genome. We propose to conduct whole genome sequencing of 5 members of the pedigree at a depth of 20X. The chosen samples are two sets of parents plus one member of an unaffected branch of the pedigree who shares the risk haplotype on chromosome 19. Sequencing of the two sets of parents will be used along with the genome-wide SNP data to impute 4 affected children giving an effect sample size of 6 affected individuals.

Study EGAS00001000305

Exome sequencing of a novel cervical cancer cell line

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

Study EGAS00001003343

Mutational Patterns in Metastatic Cutaneous Squamous Cell Carcinoma

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

Study EGAS00001003370

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

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

Study EGAS00001003372

Intratumoural Heterogeneity Underlies Distinct Therapy Responses and Treatment Resistance in Glioblastoma

Glioblastomas are the most common and lethal neoplasms of the central nervous system. Neighbouring glioma cells maintain extreme degrees of genetic and phenotypic variation that forms intratumoural heterogeneity. This genetic diversity allows the most adaptive tumour clones to show treatment resistance and conserve the tumour bulk. We aimed to model this phenomenon and test the effectiveness of several targeted therapeutic interventions to overcome therapy resistance. Heterogenic tumour masses were first deconstructed into single tumour cells, which were expanded independently as single-cell clones. SNP arrays, whole-genome and RNA sequencing, and CpG methylation analysis validated the unique molecular profile of each tumour clone, which displayed distinct pathologic features, including cell morphology, growth rate, and resistance to temozolomide and ionizing radiation. We also identified variable sensitivities to AURK, CDK, and EGFR inhibitors, consistent with the heterogenic molecular alterations that each clone harboured. These targeted therapies effectively eliminated the temozolomide- and/or irradiation-resistant clones and also parental polyclonal cells. Our findings indicate that polyclonal tumours create a dynamic environment that consists of diverse tumour elements and treatment responses. Designing targeted therapies based on a range of molecular profiles can be a more effective strategy to eradicate treatment resistance, recurrence, and metastasis.

Study EGAS00001003438

MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Colorectal_LCM

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

Study EGAS00001003455

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

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

Study EGAS00001003443

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

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

Study EGAS00001004327

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

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

Study EGAS00001004343

Developmental_Dysplasia_of_the_Hip__DDH_

The National Joint Registry for England and Wales (NJR) has become established as one of the largest repositories of clinical information about patients undergoing hip and knee replacement. The primary aim of this proposal is to examine the feasibility of adding value to that clinical information by partnering it with a DNA dataset to allow its use to be extended to the study of the underlying pathogenesis of diseases of the hip and knee, and the complications of their treatment. We will do this 1) by establishing a mechanism for accessing and linking pseudo-anonymised demographic and phenotype data between the NJR database and the proposed DNA Biobank; 2) by validating the collection of postal DNA samples from consenting NJR patients; and 3) by genotyping genome-wide 900 patients with the condition developmental dysplasia of the hip (DDH) to find variants that play a role in this heritable condition. The successful demonstration that the NJR can be used to establish a DNA Biobank will provide a powerful platform to examine the genetic basis for a wide range of joint diseases, and the complications of their treatment. This proposal is to carry out the genotyping portion of this project on the current CoreExome array

Study EGAS00001000916

Comprehensive genomic profiles of small cell lung cancer

We have sequenced the genomes of 110 small cell lung cancers (SCLC), one of the deadliest human cancers. We found bi-allelic inactivation of TP53 and RB1 in nearly all the tumors analyzed, sometimes by complex genomic rearrangements. Two tumors with wild-type RB1 had evidence of chromothripsis leading to overexpression of Cyclin D1, revealing an alternative mechanism of Rb1 deregulation. Thus, loss of the tumor suppressors TP53 and RB1 is obligatory in SCLC. We further discovered somatic genomic rearrangements of TP73 that create an oncogenic version of this gene, TP73Δex2/3. In few cases, SCLC tumors exhibited kinase gene mutations, providing a possible therapeutic opportunity for individual patients. Finally, we observed inactivating mutations in NOTCH family genes in 25% of human SCLC. Accordingly, activation of Notch signaling in a pre-clinical SCLC mouse model dramatically reduced the number of tumors and extended the survival of the mutant mice. Furthermore, neuroendocrine gene expression was abrogated by Notch activity in SCLC cells. This first comprehensive study of somatic genome alterations in SCLC uncovers several key biological processes and identifies candidate therapeutic targets in this highly lethal form of cancer.

Study EGAS00001000925

Exploring the role of mtDNA variation in Multiple Sclerosis in a large cohort of discordant monozygotic twins

Several lines of evidence indicate that mitochondrial DNA (mtDNA) variants might predispose to multiple sclerosis (MS). We examined this in 49 monozygotic (MZ) twin pairs clinically discordant for MS at study entry. Since the nuclear DNA of MZ twins is identical, our study provides a unique setting to functionally associate unique mtDNA variants and skewed heteroplasmy with MS development. In order to identify mtDNA variants in the twin cohort, we deeply sequenced the mitochondrial genome in blood of all 49 MZ twin pairs using next generation sequencing (Illumina, HiSeq or MiSeq) with an average coverage depth of ~25,000 sequences per base. Heteroplasmic variants were validated using targeted deep sequencing (TDS) (Illumina, MiSeq) in DNA isolated from blood and buccal swabs. In addition, whole blood was also available of 5 additional MS-affected or healthy siblings and 1 dizygotic twin pair discordant for MS, comprising in total 6 non-identical sib pairs. Moreover, all cases were screened for pathogenic mtDNA polymerase gamma (POLG) mutations. The data provides valuable insights in the possible involvement of mtDNA variants in the pathogenesis of MS, and gives information on the segregation of heteroplasmic variants within MZ twins and non-identical siblings and across different tissues.

Study EGAS00001001240

Clonal selection and double hit events involving tumor suppressor genes underlie relapse from total therapy

In order to gain insights into the mechanisms of disease progression in multiple myeloma (MM) patients treated with multi-agent and high-dose chemotherapy we have performed a longitudinal study of 33 patients entered into Total Therapy protocols using gene expression profiling, high resolution copy number analyses and whole exome sequencing. The results of this study show the importance of acquired mutations in MM driver genes, bi-allelic events affecting tumor suppressor genes and increased proliferation rates in driving relapse, which is enhanced by intra-clonal heterogeneity and Darwinian type clonal evolution. Branching was the most frequent evolution pattern, characterized by clonal loss and appearance of new clones. We found a higher number of lost and acquired copy number aberrations and an excess of bi-allelic inactivation of tumor suppressor genes in GEP70 high risk (HiR) cases, consistent with genomic instability being a key feature of HiR MM. In conclusion, our study further stresses the impact of known prognostic markers and driver genes on the fitness level of MM cells to survive multi-agent chemotherapy. The frequent bi-allelic loss of tumor suppressor genes in HiR cases and the potentially resulting failure of treatment with DNA damaging agents highlight the need for new therapies for HiR MM.

Study EGAS00001001810

A somatic reference standard for cancer genome sequencing with COLO829

Identification of somatic alterations in cancer has become feasible with the massive adoption of next generation sequencing. However, due to variability in sequencing and informatics pipelines, a common somatic reference is needed. We thus performed paired whole genome sequencing of a metastatic melanoma cell line (COLO829) and a matched lymphoblastoid line (COLO829BL) across three institutions (TGen, Illumina, Genome Sciences Centre at the British Columbia Cancer Agency). We performed a meta-analysis of all data, in combination with the originally reported analysis of these cell lines by Pleasance et al (PMID: 20016485), and report a somatic reference standard based on consensus events. DNA extractions, library preparation, sequencing, and analysis were separately performed at each site. Results were compiled with sequence data previously generated for the cell lines to identify true positive variants, collectively representing a somatic reference standard. Overall, common somatic events include point mutations, small insertion/deletions, and genes affected by concurrent copy number changes. We additionally show divergence of somatic mutations in COLO829 associated with differences in cell line lineage. We present this reference to the community as a standard for enabling interpretation and identification of somatic mutations across institutions and analytical pipelines.

Study EGAS00001001385

The HLA-DQβ1 insertion is a strong achalasia risk factor and displays a geospatial north-south gradient among Europeans

Idiopathic achalasia is a severe motility disorder of the esophagus and is characterized by a failure of the lower esophageal sphincter to relax due to a loss of neurons in the myenteric plexus. Most recently, we identified an 8 amino acid insertion in the cytoplasmic tail of HLA-DQβ1 as strong achalasia risk factor in a sample set from Central Europe, Italy, and Spain. Here, we tested whether the HLA-DQβ1 insertion also confers achalasia risk in the Polish and Swedish population. We could replicate the initial findings and the insertion shows strong achalasia association in both samples (Poland P = 1.84×10-04, Sweden P = 7.44×10-05). Combining all five European data sets – Central Europe, Italy, Spain, Poland, and Sweden – the insertion is achalasia associated with Pcombined = 1.67×10-35. In addition, we observe that the frequency of the insertion shows a geospatial north-south gradient. The insertion is less common in northern (around 6-7% in patients and 2% in controls from Sweden and Poland) compared to southern Europeans (around 16% in patients and 8% in controls from Italy) and shows a stronger attributable risk (AR) in the southern European population. Our study provides evidence that the prevalence of achalasia may differ between populations.

Study EGAS00001001515

Whole exome sequencing on primary retinoblastoma tissues and matching lymphocyte DNA.

Retinoblastoma is a rare childhood cancer initiated by RB1 mutation or MYCN amplification, while additional alterations may be required for tumor development. However, the view on single nucleotide variants is very limited. To better understand oncogenesis, we determined the genomic landscape of retinoblastoma. We performed exome sequencing of 71 retinoblastomas and matched blood DNA. Next, we determined the presence of single nucleotide variants, copy number alterations and viruses. Aside from RB1, recurrent gene mutations were very rare. Only a limited fraction of tumors showed BCOR (7/71, 10%) or CREBBP alterations (3/71, 4%). No evidence was found for the presence of viruses. Instead, specific somatic copy number alterations were more common, particularly in patients diagnosed at later age. Recurrent alterations of chromosomal arms often involved less than one copy, also in highly pure tumor samples, suggesting within-tumor heterogeneity. Our results show that retinoblastoma is among the least mutated cancers and signify the extreme sensitivity of the childhood retina for RB1 loss. We hypothesize that retinoblastomas arising later in retinal development benefit more from subclonal secondary alterations and therefore, these alterations are more selected for in these tumors. Targeted therapy based on these subclonal events might be insufficient for complete tumor control.

Study EGAS00001001690

Ethiopia_Genome_Project__low_coverage_

The Ethiopian area stands among the most ancient ones ever occupied by human populations and their ancestors. Particularly, according to archaeological evidences, it is possible to trace back the presence of Hominids up to at least 3 million years ago. Furthermore, the present day human populations show a great cultural, linguistic and historic diversity which makes them essential candidate to investigate a considerable part of the African variability. Following the typing of 300 Ethiopian samples on Illumina Omni 1M (see Human Variability in Ethiopia project, previously approved by the Genotyping committee) we now have a clearer idea on which populations living in the area include the most of the diversity.This project therefore aims to sequence the whole genome of 300 individuals at low (4-8x) depth belonging to the six most representative populations of the Ethiopian area to produce a unique catalogue of variants peculiar of the North East Africa. Furthermore 6 samples (one from each population) will also be sequenced at high (30x) depth to ensure full coverage of the diversity spectrum.The retrieved variants will be of great help in evaluating the demographic dynamics of those populations as well as shedding light on the migrations out of Africa.

Study EGAS00001000238

UK10K OBESITY TWINSUK

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. These TwinsUK samples are a subset of the TwinsUK cohort samples with a BMI>=40 that will be exome sequenced by the Obesity group.

Study EGAS00001000306

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

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

Study EGAS00001003689

Pancreatic, Small-intestinal and Pulmonary Neuroendocrine Tumors

Purpose: The primary origin of neuroendocrine tumor metastases can be difficult to determine by histopathology alone, but is critical for therapeutic decision making. DNA methylation-based profiling is now routinely used in the diagnostic workup of brain tumors. This has been enabled by the availability of cost-efficient array-based platforms. We have extended these efforts to augment histopathological diagnosis in neuroendocrine tumors. Experimental Design and Results: We compiled data of 69 small-intestinal, pulmonary, and pancreatic neuroendocrine tumors. These data were used to build a ridge regression calibrated random forest classification algorithm (NEN-ID) that predicts the origin of tumor samples with high accuracy (&gt; 95%). The model was validated during 3x3 nested cross validation and tested in a local (n=26) and external (n=172) cohort. In addition, we show that our diagnostic approach is robust across a range of possible confounding experimental parameters such as tumor purity and array quality. A software infrastructure and online user interface was built to make the model available to the scientific community. Conclusions: This DNA methylation-based prediction model can be used in the workup for patients with neuroendocrine tumors of unknown primary. To facilitate validation and clinical implementation, we provide a user-friendly, publicly available web-based version of NEN-ID.

Study EGAS00001004878

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

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

Study EGAS00001002654

The Prediction and Prevention of Preeclampsia

The Prediction and Prevention of Preeclampsia (PREDO) is a prospective birth cohort study of Finnish women who were pregnant between 2005 and 2010 and their children. The PREDO study cohort was set up to identify novel risk factors and biomarkers in pregnant women associated with the development of preeclampsia and intrauterine growth restriction (IUGR). Two groups of pregnant women were enrolled during the first trimester of their pregnancy: first, pregnant women with a known clinical risk factor status for preeclampsia and IUGR (N=1079; 969 of these women had at least one and 110 had none of the known risk factors for preeclampsia and IUGR), and second, pregnant women who volunteered to participate regardless of their risk factor status for preeclampsia and IUGR (N=3698). Both samples were monitored during pregnancy with four clinical visits and bi-weekly self-reports. The post-delivery follow-up has taken place at approximately 2 weeks, 6 months, and 3.5 years after the delivery. The most recent follow-up started in 2016 and is ongoing. In the high-risk sample, cord blood samples were collected and genome-wide genotyping has been conducted with Illumina OmniExpressExome 1.2 and genome-wide methylation status has been analysed with Illumina 450k HM array.

Study EGAS00001001898

The_Genomic_Advances_in_SepsisGAinSRNA_seq

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

Study EGAS00001003772

Extreme phenotypes define epigenetic and metabolic signatures in cardiovascular diseases.

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

Study EGAS00001003780

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

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

Study EGAS00001003860

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

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

Study EGAS00001002406

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

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

Study EGAS00001002700

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

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

Study EGAS00001005053

THE GENOMIC LANDSCAPE OF ACTINIC KERATOSIS

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

Study EGAS00001004243

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

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

Study EGAS00001004291

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

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

Study EGAS00001002875

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

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

Study EGAS00001004018

23014 results for "ega"

in 58.53 milliseconds.