3244 results for "*geneti*"

in 71.20 milliseconds.

• Genomic profiling of paediatric glioma cell lines

  Although paediatric gliomas resemble their adult counterparts in many ways, there appear to be distinct clinical and biological differences. One important factor hampering the development of new targeted therapies is the relative lack of cell lines derived from childhood glioma patients, as it is unclear whether the well-established adult lines commonly used are representative of the underlying molecular genetics of childhood tumours. We have carried out a detailed molecular profiling of paediatric glioma cell lines KNS42, SF188, UW479, RES259 and RES186.

  Study EGAS00001003006

• Molecular and Population Genetics (MPG) research at ICR DAC and POLA DAC

  Dac EGAC00001000328

• CTCF/cohesin-binding sites are frequently mutated in cancer

  Cohesin is present in almost all active enhancer regions, where it is associated with transcription factors. Cohesin frequently colocalizes with CTCF (CCCTC-binding factor), affecting genomic stability, expression and epigenetic homeostasis. Cohesin subunits are mutated in cancer, but CTCF/cohesin-binding sites (CBSs) in DNA have not been examined for mutations. Here we report frequent mutations at CBSs in cancers displaying a mutational signature where mutations in A•T base pairs predominate. Integration of whole-genome sequencing data from 213 colorectal cancer (CRC) samples and chromatin immunoprecipitation sequencing (ChIP-exo) data identified frequent point mutations at CBSs. In contrast, CRCs showing an ultramutator phenotype caused by defects in the exonuclease domain of DNA polymerase ɛ (POLE) displayed significantly fewer mutations at and adjacent to CBSs. Analysis of public data showed that multiple cancer types accumulate CBS mutations. CBSs are a major mutational hotspot in the noncoding cancer genome.

  Study EGAS00001003010

• Allele Balance Bias Identifies Systematic Genotyping Errors and False Disease Associations

  In recent years, Next Generation Sequencing (NGS) has become a cornerstone of clinical genetics and diagnostics. Many clinical applications require high precision, especially if rare events such as somatic mutations in cancer or genetic variants causing rare diseases need to be identified. Although random sequencing errors can be modeled statistically and deep sequencing minimizes their impact, systematic errors remain a problem even at high depth of coverage. Understanding their source is crucial to increase precision of clinical NGS applications. In this work, we studied the relation between recurrent biases in allele balance (AB), systematic errors and false positive variant calls across a large cohort of human samples analyzed by whole exome sequencing (WES). We have modeled the allele balance distribution for biallelic genotypes in 987 WES samples in order to identify positions recurrently deviating significantly from the expectation, a phenomenon we termed allele balance bias (ABB). Furthermore, we have developed a genotype callability score based on ABB for all positions of the human exome, which detects false positive variant calls that passed state-of-the-art filters. Finally, we demonstrate the use of ABB for detection of false associations proposed by rare variant association studies (RVAS).

  Study EGAS00001003027

• The evolutionary history of human colitis-associated colorectal cancer

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

  Study EGAS00001003028

• Patient-derived neuroblastoma model system OHC-NB1

  Accurate disease modeling presents a bottleneck against effective translation of novel neuroblastoma therapies into the clinic. Recently, considerable intra-tumor genetic heterogeneity has been observed in neuroblastoma, which calls for preclinical models that reflect this feature. Flexibility in a testing platform is also desirable to support assessment of different endpoints and tumor cell characteristics. We present the novel neuroblastoma model OHC-NB1, propagated from a bone marrow metastasis from a patient with INRG stage M, MYCN-amplified neuroblastoma at first diagnosis. To model different aspects of the disease, several preclinical models were developed from the bone marrow aspirate: a tissue plastic-adherent cell line, a 3-dimensional (3D) spheroid culture, and cell line-derived xenografts that were subcutaneously passaged in mice. Typical aspects of MYCN-amplified neuroblastoma, such as MYCN amplification in double minutes, 1p deletion, 17q gain, and a diploid karyotype, are present in OHC-NB1 and persist in all culture types. Whole exome sequencing revealed intra-tumor heterogeneity in the bone marrow metastasis reflected in the culture models. The set of OHC-NB1 cultures with a genetically heterogeneous MYCN-amplified neuroblastoma background provide a flexible preclinical testing system that recapitulates important features of metastatic neuroblastoma.

  Study EGAS00001003031

• SLC9A3R1 variant associated with age-related hearing loss

  Age-related hearing loss (ARHL) is the most common sensory impairment in the elderly affecting millions of people worldwide. To shed light on the genetics of ARHL, a large cohort of 464 Italian patients has been deeply analyzed at clinical and molecular level. A missense variant in SLC9A3R1 has been identified in two unrelated ARHL patients. A knock-in zebrafish model confirmed the pathogenic effect of the variant.

  Study EGAS00001003072

• Exome sequencing in bipolar disorder families

  Bipolar disorder (BD) is a major psychiatric disorder affecting around 1% of the global population. BD is characterized by recurrent manic and depressive episodes, and has an estimated heritability of around 70% Research has identified the first BD susceptibility genes. However, the underlying pathways and regulatory networks remain largely unknown. Since research has shown that the cumulative impact of common alleles with small effect appears to explain only around 25-38% of the phenotypic variance for BD, rare variants of high penetrance may also contribute to BD risk. The present study investigated the role of rare, nonsynonymous, and potentially functional variants via whole exome sequencing in 15 BD cases from two large, multiply affected families from Cuba. The high prevalence of BD in these pedigrees renders them promising in terms of the identification of genetic risk variants with large effect sizes. In addition, SNP array data were used to calculate polygenic risk scores for affected and unaffected family members. After correction for multiple testing, polygenic risk scores for common, BD-associated genetic variants were not significantly higher in BD cases than in healthy relatives. Exome sequencing identified a total of 17 rare and potentially damaging variants in 17 genes. The identified variants were shared by all investigated BD cases in the respective pedigree. The most promising variant was located in the gene SERPING1 (p.L349F), which has been reported as a genome-wide significant risk gene for schizophrenia. A further promising variant was located in CAPN2, which has been implicated in synaptic plasticity. The present data suggest novel candidate genes for BD susceptibility, and may facilitate the discovery of disease-relevant pathways and regulatory networks.

  Study EGAS00001003085

• Histone acetylome-wide association study on tuberculosis infection

  Epigenetic mechanisms are thought to play an important role in the pathogenesis of infectious diseases. Here, we describe the first histone acetylome-wide association study (HAWAS) of an infectious disease, based on genome-wide H3K27 acetylation profiling of peripheral granulocytes and monocytes from subjects with active Mycobacterium tuberculosis (Mtb) infection and healthy controls (135 ChIP-seq datasets).

  Study EGAS00001003118

• Heritable pulmonary arterial hypertension in a large Iberian family

  Data from a study of a large Iberian family (n=65 subjects, 5 generations) affected by pulmonary arterial hypertension (PAH) and segregating with the BMPR2 missense mutation p.Arg491Gln (rs137852749, c.1472G>A). PAH is a rare disease characterized by an abnormal rise in mean pulmonary arterial pressure (> or equal to 25 mmHg at rest), which, in turn, leads to a progressive increase in pulmonary vascular resistance and ultimately to death, due to right ventricular failure. Heritable PAH has an overall prevalence below 1 case per million adults and is defined by either the presence of a known genetic defect linked to the disease or a positive family history. Heritable PAH is inherited as an autosomal dominant disease. However, not all BMPR2 mutation carriers develop the disease, highlighting the presence of reduced penetrance. In this family, there are 22 mutation carriers from which 8 were diagnosed with heritable PAH and the other 14 were healthy at the time of examination.

  Study EGAS00001003123

• Whole-exome sequencing of extranodal NK/T cell lymphoma

  Extranodal NK/T cell lymphoma (ENKTCL) is an aggressive Epstein-Barr virus-associated T/NK neoplasm that predominantly affects Asians. To explore the causative somatic events, we conducted a comprehensive genetic analysis of 19 ENKTCL patients including whole-exome sequencing (N = 16). Commonly deregulated gene pathways in ENKTCL included epigenetic modifiers followed by human leukocyte antigen (HLA) and related genes and JAK-STAT pathway.

  Study EGAS00001004357

• Whole exome sequencing of small cell neuroendocrine cancer of the cervix

  In order to improve treatment selection for high grade neuroendocrine carcinomas of the cervix (NECC), we performed a comparative genomic analysis between this rare tumor type and other cervical cancer types, as well as extra-cervical neuroendocrine small cell carcinomas of the lung and bladder. We performed whole exome sequencing on fresh-frozen tissue from 15 NECCs and matched normal tissue. We then identified mutations and copy number variants using standard analysis pipelines. Published mutation tables from cervical cancers and extra-cervical small cell carcinomas were used for comparative analysis. Descriptive statistical methods were used and a two-sided threshold of P < .05 was used for significance. In the NECC cohort, we detected a median of 1.7 somatic mutations per megabase (range 1.0-20.9). PIK3CA p.E545K mutations were the most frequency observed oncogenic mutation (4/15 tumors, 27%). Activating MAPK pathway mutations in KRAS (p.G12D) and GNAS (p.R201C) co-occurred in two tumors (13%). In total we identified PI3-kinase or MAPK pathway activating mutations in 67% of NECC. When compared to NECC, lung and bladder small cell carcinomas exhibited a statistically significant higher rate of coding mutations (P < .001 for lung; P = .001 for bladder). Mutation of TP53 was uncommon in NECC (13%) and was more frequent in both lung (103 of 110 tumors [94%], P < .001) and bladder (18 of 19 tumors [95%], P < .001) small cell carcinoma. These comparative genomics data suggest that NECC may be genetically more similar to common cervical cancer subtypes than to extra-cervical small cell neuroendocrine carcinomas of the lung and bladder. These results may have implications for the selection of cytotoxic and targeted therapy regimens for this rare disease.

  Study EGAS00001003142

• High Altitude Pulmonary Hypertension

  The Central Asian Kyrgyz highland population provides a unique opportunity to address genetic diversity and understand the genetic mechanisms underlying hypoxia-induced high altitude pulmonary hypertension (HAPH). While a significant fraction of the population is unaffected, there are susceptible individuals who display HAPH in the absence of any lung, cardiac or hematologic disease. We report herein the analysis of the whole genome sequencing of healthy individuals compared with HAPH patients and other controls.

  Study EGAS00001003171

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

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

  Study EGAS00001003177

• Genome-wide association study of cervical cancer in East Asian populations

  The development of cervical cancer is initiated by human papillomavirus (HPV) infection, and involves both viral and host genetic factors. Genome-wide association studies (GWAS) of cervical cancer have identified associations in the HLA locus and two loci outside HLA, but the principal genes that control infection and pathogenesis have not been identified. In the present study, we performed GWAS of cervical cancer in East Asian populations, involving 2609 cases and 4712 controls in the discovery stage, 1461 cases and 3295 controls in the follow-up stage.

  Study EGAS00001003199

• FBSeq: RNA sequencing of human fetal brain.

  RNA sequencing of 120 human fetal brains (aged 12-19 post-conception weeks) was performed as part of a Medical Research Council (U.K.) funded project (MR/L010674/2) to investigate genetic effects on gene expression in the developing human brain and their role in neuropsychiatric disorders. Human fetal brain tissue from elective terminations of pregnancy was provided by the Human Developmental Biology Resource (HDBR) (http://www.hdbr.org), with ethical approval (#08/H0712/34), and with consent for use of fetal material for medical research provided by female donors. Total RNA was extracted from half of the available brain tissue from each fetus using Tri-Reagent and treated with TURBO DNase. RNA-Seq libraries were prepared using 1µg of purified total RNA and the Illumina TruSeq Stranded Total RNA Library Prep kit, following ribosomal RNA depletion. Libraries were sequenced on the Illumina HiSeq 2500 or HiSeq 4000 systems, generating at least 50 million read pairs per sample.

  Study EGAS00001003214

• An exome sequencing approach to defining the genetic risk factors for Achilles tendinopathy

  Musculoskeletal soft tissue injuries are complex phenotypes. Case control association studies using the candidate gene approach have predominantly been used to identify risk loci for these injuries. However, the ability to identify all risk conferring variants using this approach alone is unlikely. Therefore, this study aimed to define the genetic profile of these injuries using whole exome sequencing and a customized analysis pipeline. The exomes of ten exemplar asymptomatic controls and ten exemplar cases with Achilles tendinopathy were individually sequenced using a platform that included the coverage of the untranslated regions and miRBase miRNA genes. These results provide a proof of concept for the use of a customized pipeline for the exploration of a larger genomic dataset. This approach, using previous research to guide a targeted analysis of the data has generated new genetic signatures in the biology of musculoskeletal soft tissue injuries.

  Study EGAS00001003234

• Epigenetics/genetics of the patient-derived xenografts of pediatric T-cell leukemia

  We compared 12 pediatric T cell acute lymphoblastic leukemias collected at initial diagnosis and relapse with their corresponding PDX models. The analysis was performed on genomic level (whole genome sequencing (WGS), whole exome sequencing (WES), multiplex ligation probe amplification (MLPA), targeted sequencing) and the epigenetic level (DNA methylation, Assay for Transposase-Accessible Chromatin sequencing (ATAC-seq)). In sum, this study underlines the remarkable genomic stability, and for the first time documents the preservation of the epigenomic landscape in T-ALL-derived PDX models.

  Study EGAS00001003248

• 1 Intratumoral genetic heterogeneity and clonal evolution following neoadjuvant chemoradiotherapy (nCRT) in locally advanced rectal tumors.

  Purpose: Neoadjuvant chemoradiotherapy (nCRT) is the standard treatment forpatients with locally advanced rectal cancer. However, little is known about the effect of nCRT in the mutational landscape and intratumoral genetic heterogeneity (ITGH) of rectal tumors. Moreover, response to nCRT vary substantially among patients and the genetic determinants of treatment response are not fully understood. Here, we sought to determine the impact of nCRT in the mutational landscape and ITGH of rectal tumors and to identify putative drivers of therapeutic resistance. Patients and Methods: We analyzed whole-exome sequencing (WES) and clinical data from 79 primary non-treated rectal cancers obtained from The Cancer Genome Atlas (TCGA). We also compared the mutational landscape of a matched set of peripheral blood cells, irradiated tumor-adjacent colonic mucosa and pre and post-treatment tumor samples to determine the iatrogenic effects of nCRT and its impact in tissue genetic heterogeneity. Finally, we performed WES of 7 matched pre and post-treatment tumor samples to examine how nCRT affects the clonal architecture of rectal tumors. Results: We show that rectal tumors exhibit remarkable ITGH, which increases with disease progression. We also show that nCRT, per se, does not introduce novel somatic mutations, nor alters the clonal architecture of the irradiated colonic mucosa. Instead, we show that it acts as a potent selective pressure, favoring the expansion of tumor cell subpopulations more prone to resist therapy and increasing ITGH in the residual tumor. Conclusion: Our results reveal the highly heterogeneous and dynamic clonal architecture of rectal tumors undergoing nCRT. Our results also uncover putative genetic determinants of response to nCRT and a possible role for ITGH in predicting treatment outcome.

  Study EGAS00001003250

• 46 CLL Whole Genome Sequencing Study

  Chronic lymphocytic leukaemia (CLL) consists of two biologically and clinically distinct subtypes defined by the abundance of somatic hypermutation (SHM) affecting the Ig variable heavy-chain locus (IgHV). The molecular mechanisms underlying these subtypes are incompletely understood. Here, we present a comprehensive whole-genome sequencing analysis of somatically acquired genetic events from 46 CLL patients, including a systematic comparison of coding and non-coding single nucleotide variants, copy number variants and structural variants, regions of kataegis and mutation signatures between IgHVmut and IgHVunmut subtypes.

  Study EGAS00001003254

• Mutated H3 Histones Drive Human Pre-Leukemic Hematopoietic Stem Cell Expansion And Promote Leukemic Aggressiveness

  Leukemogenesis is a stepwise progression from mutated, pre-neoplastic hematopoietic stem cells (HSCs) to full-blown leukemia. Our ability to prevent or treat acute myeloid leukemia (AML) is limited by our incomplete understanding of the epigenetic disruption that is central to this process, including improper histone methylation. Comprehensive analysis of 16 histone H3 genes in 434 primary AML samples identified Q69H, A26P, R2Q, R8H and K27M/I mutations (1.6%), and a higher incidence in secondary AML (s-AML) (9%). We establish that these mutations are important early events in leukemogenesis. They occur in pre-leukemic HSCs, increase the frequency of functional human HSCs, and alter differentiation. The mutations are present in the major leukemic clones in primary samples, and the mutant histones amplify leukemic aggressiveness with increased proliferation, expansion of leukemic progenitor and blast cells, and superior competitiveness in vivo. These effects are dependent on the specific mutation. Genome-wide analysis of K27 mutants revealed increased expression of genes involved in erythrocyte and myeloid differentiation with a corresponding decrease in histone H3 K27 tri-methylation and increase in K27 acetylation. The functional impact of histone mutations is independent of RUNX1 mutations, although they can co-occur. These data establish the involvement of H3 mutations as initial drivers of pre-cancerous stem cell expansion and leukemogenesis.

  Study EGAS00001003288

• Personalised Mapping of Tumour Development in Synchronous Colorectal Cancer Patients

  We performed an in-depth characterisation of 12 tumours from 3 syCRC patients by analysing histopathological, whole genome sequencing and RNA-sequencing data. We assessed the extent of genetic overlap between synchronous tumours and examined associations between clinicopathological information and the molecular, microbial and immune features of each tumour genome.

  Study EGAS00001004413

• Genomic and transcriptomic determinants of therapy resistance and immune landscape evolution during anti-EGFR treatment in colorectal cancer

  Despite biomarker stratification, the anti-EGFR antibody cetuximab is only effective against a subgroup of colorectal cancers (CRC). This genomic and transcriptomic analysis of the cetuximab resistance landscape in 35 RAS wild-type CRCs identified associations of NF-1 and non-canonical RAS/RAF-aberrations with primary resistance and validated transcriptomic CRC-subtypes as non-genetic predictors of benefit. 64% of biopsies with acquired resistance harbored no genetic resistance drivers. Most of these had switched from a cetuximab-sensitive transcriptomic subtype pre-treatment to a fibroblast- and growth factor-rich subtype at progression. Fibroblast-supernatant conferred cetuximab resistance in vitro, confirming a major role for non-genetic resistance through stromal remodeling. Cetuximab treatment increased cytotoxic immune infiltrates and PDL1 and LAG3 immune-checkpoint expression, potentially providing opportunities to treat cetuximab-resistant CRCs with immunotherapy.

  Study EGAS00001003367

• 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

• Detection of mutational patterns in cell free DNA (cfDNA) of colorectal cancer by custom amplicon sequencing

  Monitoring the mutational patterns of solid tumors during cancer therapy is a major challenge in oncology. Analysis of mutations in cell free (cf) DNA offers a non-invasive approach to detect mutations that may be prognostic for disease survival or predictive for primary or secondary drug resistance. A main challenge for the application of cfDNA as a diagnostic tool is the diverse mutational landscape of cancer. Here, we developed a flexible end-to-end experimental and bioinformatics workflow to analyze mutations in cfDNA using custom amplicon sequencing. Our approach relies on open software tools to select primers suitable for multiplex PCR using minimal cfDNA as input. In addition, we developed a robust linear model to identify specific genetic alterations from sequencing data of cfDNA. We used our method to design a custom amplicon panel suitable for detection of hotspot mutations relevant for colorectal cancer and analyzed mutations in serial cfDNA samples from 34 patients with advanced colorectal cancer. Our data demonstrates that recurrent and patient-specific mutational patterns can be reliably detected for the majority of patients. Furthermore, we show that the allele frequency of mutations in cfDNA correlates well with disease progression. Finally, we demonstrate that monitoring of cfDNA can outperform the predictive power of currently used tumor markers and reveal mechanisms of resistance to anti-EGFR antibody treatment.

  Study EGAS00001003382

• Dissecting features of epigenetic variants underlying cardiometabolic risk using full-resolution epigenome profiling in regulatory elements

  Sparse profiling of CpG methylation in blood by microarrays have identified epigenetic links to common diseases. We apply methylC-capture sequencing (MCC-Seq) in a clinical population of ~200 adipose tissue and matched blood samples (Ntotal ~400), providing high-resolution methylation profiling (>1.3M CpGs) at regulatory elements. We link methylation to cardiometabolic risk through associations to circulating plasma lipid levels and identify lipid-associated CpGs with unique localization patterns in regulatory elements. We show distinct features of tissue-specific versus tissue-independent lipid-linked regulatory regions by contrasting with parallel assessments in ~800 independent adipose tissue and blood samples from the general population. We follow-up on adipose-specific regulatory regions under (1) genetic and (2) epigenetic (environmental) regulation via integrational studies. Overall, the comprehensive sequencing of regulatory element methylomes reveals a rich landscape of functional variants linked genetically as well as epigenetically to plasma lipid traits.

  Study EGAS00001003415

• 16S-based fecal microbiota composition of the Milieu Intérieur Cohort

  The Milieu Intérieur project is a population-based study coordinated by the Institut Pasteur, Paris. The objective is to dissect the interplay between genetics and environment and their impact on the immune system. One of the primary goals of the project is to ascertain the role of the microbiota in regulating immune programs.

  Study EGAS00001003419

• STREP GENE: Genetics and Severe Streptococcal Infections

  This study investigated host genetic susceptibility to invasive group A streptococcal disease in the UK. Cases were either survivors recruited retrospectively through a patient group or identified from the a tissue bank at Imperial College London. Those recruited through the patient group had survived an episode of invasive GAS disease at a UK hospital since 1980 with microbiological confirmation obtained either through Public Health England or from the treating hospital. Those identified from the tissue bank had been diagnosed with invasive GAS disease at the Hammersmith Hospital, London, UK, since 2006. All cases aged less than 65 years without long-term medical conditions genome-wide genotyped using either the Illumina HumanCore or the Global Screening Array platform.

  Study EGAS00001003421

• 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

• Functional and genomic heterogeneity of long-term self-renewing compartment as the origin of treatment resistance in pancreatic tumors

  Adaptive drug-resistance mechanisms allow human tumors to evade treatment through selection and expansion of treatment-resistant clones. Here, studying clonal evolution of tumor cells derived from human pancreatic tumors, we demonstrate that in vitro cultures and in vivo tumors are maintained by a common set of tumorigenic cells that can be used to establish clonal replica tumors (CRTs), large cohorts of animals bearing human tumors with identical clonal composition. Using CRTs to conduct quantitative assessments of adaptive responses to therapeutics, we uncovered a multitude of functionally heterogeneous subpopulations of cells with differential degrees of drug sensitivity. High-throughput isolation and deep characterization of unique clonal lineages showed genetic and transcriptomic diversity underlying functionally diverse subpopulations. Molecular annotation of gemcitabine-naïve clonal lineages with distinct responses to treatment in the context of CRTs generated signatures that can predict the response to chemotherapy, representing a potential biomarker to stratify patients with pancreatic cancer. RNA-Seq data on these clones are available on EBI ENA, under accession number: PRJEB27735.

  Study EGAS00001003442

• 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

• Whole genome sequencing of 25 South African individuals with myasthenia gravis

  This study investigated the genetic basis of an ophthalmoplegic subphenotype of MG (OP-MG) which occurs in a proportion of myasthenics with juvenile symptom onset and African genetic ancestry. This OP-MG subphenotype is characterised by treatment resistant weakness of the extraocular muscles (EOMs) and the pathogenetic mechanism(s) underlying its development is unknown. The aim of this study was to use a hypothesis-generating genome-wide case-control analysis to identify candidate OP-MG susceptibility genes and pathways. Whole genome sequencing (WGS) was performed on 25 AChR-antibody positive myasthenic individuals of African genetic ancestry sampled from the phenotypic extremes: 15 with OP-MG (cases) and 10 individuals with control MG (EOM treatment-responsive).

  Study EGAS00001003462

• Exome sequencing of stroke cases with good or bad recovery three months after stroke

  With the aim to identify rare genetic variation involved in stroke recovery, this is a pilot study on extreme phenotypes of good and bad recovery. Exome sequencing was performed on 41 and 49 samples in the bad and good recovery subset, respectively. Rare Variant association Analysis was performed on the exome sequencing results.

  Study EGAS00001003463

• Adenoma to Carcinoma transition in Colorectal Cancer

  Although colorectal carcinoma development is well accepted as a step-wise and sequential process, it is commonly investigated using non-sequential adenoma and carcinoma samples. Here we compared molecular pathologies of the earliest carcinoma stages with those in their adenomatous precursors to unravel common and divergent genetic traits between matched neoplasms with a shared evolutionary history but different contemporary biology.

  Study EGAS00001003468

• EXCEED Study

  EXCEED is a longitudinal population-based cohort which facilitates investigation of genetic, environmental and lifestyle-related determinants of a broad range of diseases and of multiple morbidity through data collected at baseline and via electronic healthcare record linkage. Recruitment has taken place in Leicester, Leicestershire and Rutland since 2013 and is ongoing, with approximately 10,000 participants aged 30-69 to date. The population of Leicester is diverse and additional recruitment from the local South Asian community is ongoing.

  Study EGAS00001003499

• Exome sequencing of synchronous colorectal cancers

  Approximately 4% of colorectal cancer (CRC) patients have at least two simultaneous cancers in the colon. Due to the shared environment, these synchronous CRCs (SCRCs) provide a unique setting to study colorectal carcinogenesis. Understanding whether these tumors are genetically similar or distinct is essential when designing therapeutic approaches. We performed exome sequencing of 47 primary cancers and corresponding normal samples from 23 patients. Additionally, we carried out a comprehensive mutational signature analysis to assess whether tumors had undergone similar mutational processes and the first immune cell score analysis (IS) of SCRC to analyze the interplay between immune cell invasion and mutation profile in both lesions of an individual. The tumor pairs shared only few mutations, favoring different mutations in known CRC genes and signaling pathways, and displayed variation in their signature content. Two tumor pairs had discordant mismatch repair statuses. In majority of the pairs, IS varied between primaries. Differences were not explained by any clinicopathological variable or mutation burden.

  Study EGAS00001003474

• Mutational patterns and regulatory networks in epigenetic subgroups of meningioma (H033)

  DNA methylation patterns delineate clinically relevant subgroups of meningioma. We previously established the six meningioma methylation classes (MC) benign-1, 2, 3, intermediate-A, B and malignant. Here, we set out to identify subgroup-specific mutational patterns and pathway regulation. Whole-genome-sequencing was performed on 62 samples across all MCs and WHO grades from 62 patients with matched blood control, including 40 sporadic and 22 radiation-induced (Mrad) meningiomas. RNA sequencing was added for 18 cases and chromatin-immunoprecipitation for the enhancer mark H3K27ac followed by sequencing (ChIP-seq) for 16 samples. Besides the known mutations in meningioma, structural alterations were found to contribute to the spectrum of mechanisms inactivating NF2 in sporadic meningioma similar to previous reports for Mrad. Aberrations of DMD were found to be enriched in MCs with NF2 mutations and DMD was among the most differentially upregulated genes in NF2 mutant compared to NF2 wild-type cases. The mutational signature AC3 was detected in both sporadic meningioma and Mrad, but distributed across the genome in sporadic cases and enriched near genomic breakpoints in Mrad. In general, the malignant MC presented a significantly higher exposure to AC3 and higher genomic instability beyond the mutational load than the other MCs. Pathway analysis of the ChIP-seq clusters revealed that the FOXM1 is most differentially activated in high-grade cases, along with super enhancer recruitment near HOX genes and respective upregulation of expression. This data further elucidate the biological mechanisms differentiating meningiomas of different WHO grade and epigenetic subgroups and suggest leveraging the genomic instability as novel therapeutic targets.

  Study EGAS00001003481

• Growth Hormone (GH) -secreting Pituitary Adenoma

  We profiled the somatic landscape of 21 growth hormone (GH) -secreting pituitary adenomas using somatic copy-number alteration (SCNA), whole-genome sequencing (WGS), bisulfate sequencing, and transcriptome approaches. See details in Valimaki et al. Genetic and epigenetic characterization of growth hormone (GH) - secreting pituitary tumors. Manuscript in preparation, 2019.

  Study EGAS00001003488

• 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

• The landscape of LAM disease

  Tumor biopsies from LAM disease were retrospectively analyzed by multiple techniques to characterize the alterations in patients ,to elucidate the landscape of genetic/genomic alterations.

  Study EGAS00001003534

• Clonal evolution in myeloma: the impact of maintenance lenalidomide and depth of response on the genetics and sub-clonal structure of relapsed disease in uniformly treated newly diagnosed patients

  The emergence of treatment resistant sub-clones is a key feature of relapse in multiple myeloma. Therapeutic attempts to extend remission and prevent relapse include the maximisation of response and use of maintenance therapy. We used whole exome sequencing to study the genetics of paired presentation and relapse samples from 56 newly diagnosed patients, following induction therapy, randomised to receive either lenalidomide maintenance or observation as part of the Myeloma XI trial. Patients included were considered high risk, relapsing within 30 months of maintenance randomisation. Patients achieving a complete response had predominantly branching evolutionary patterns leading to relapse, characterised by a greater mutational burden, an altered mutational profile, bi-allelic inactivation of tumour suppressor genes, and acquired structural aberrations. Conversely, in patients achieving a partial response the evolutionary features were predominantly stable with a similar mutational and structural profile. There were no significant differences between patients relapsing after maintenance lenalidomide vs observation. This study shows that the depth of response is a key determinant of the evolutionary patterns seen at relapse.

  Study EGAS00001003539

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

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

  Study EGAS00001003547

• Rare coding variants in lupus risk genes

  Systemic lupus erythematosus (SLE) is the prototypic systemic autoimmune disease. It is thought that many common variant gene loci of weak effect act additively to predispose to common autoimmune diseases, while the contribution of rare variants remains unclear. Here we describe that rare coding variants in lupus-risk genes are present in most SLE patients and healthy controls. We demonstrate the functional consequences of rare and low frequency missense variants in the interacting proteins BLK and BANK1, which are present alone, or in combination, in a substantial proportion of lupus patients. The rare variants found in patients, but not those found exclusively in controls, impair suppression of IRF5 and type-I IFN in human B cell lines and increase pathogenic lymphocytes in lupus-prone mice. Thus, rare gene variants are common in SLE and likely contribute to genetic risk.

  Study EGAS00001003548

• Genetic modification of primary human B cells to model high-grade lymphoma

  Study EGAS00001003560

• Genetic Alterations in Benign Breast Biopsies of Subsequent Breast Cancer Patients

  Fibrocystic changes are associated with an increased risk of breast cancer (approximately 1.5-2 times that of the general population). Although the relative risk seems low, because of the high frequency of such changes, the absolute risk is significant. Genetic alterations have been found in fibrocystic changes with or without epithelial changes, suggesting that critical oncogenic events are occurring at an early stage. We investigated a unique collective of 17 breast cancer patients who, prior to the diagnosis of invasive breast cancer, underwent open surgical biopsy showing fibrocystic changes of the breast. Massive parallel sequencing targeting genes frequently mutated in breast cancer was performed on the fibrocystic breast tissue as well as the ensuing cancer tissue.

  Study EGAS00001003563

• Probabilistic modeling of personalized drug combinations from integrated chemical screen and molecular data in sarcoma

  Background: Cancer patients with advanced disease exhaust available clinical regimens and lack actionable genomic medicine results, leaving a large patient population without effective treatments options when their disease inevitably progresses. To address the unmet clinical need for evidence-based therapy assignment when standard clinical approaches have failed, we have developed a probabilistic computational modeling approach which integrates sequencing data with functional assay data to develop patient-specific combination cancer treatments. Methods: Tissue taken from a murine model of alveolar rhabdomyosarcoma was used to perform single agent drug screening and DNA/RNA sequencing experiments; results integrated via our computational modeling approach identified a synergistic personalized two-drug combination. Cells derived from the tumor were allografted into mouse models and used to validate the personalized two-drug combination. Computational modeling of single agent drug screening and RNA sequencing of multiple heterogenous sites from a single patient’s epithelioid sarcoma identified a personalized two-drug combination effective across all tumor regions. The heterogeneity-consensus combination was validated in a xenograft model derived from the patient’s primary tumor. Cell cultures derived from human and canine undifferentiated pleomorphic sarcoma were assayed by drug screen; computational modeling identified a resistance-abrogating two-drug combination common to both cell cultures. This combination was validated in vitro via a cell regrowth assay. Results Our computational modeling approach addresses three major challenges in personalized cancer therapy: synergistic drug combination predictions (validated in vitro and in vivo in a genetically engineered murine cancer model), identification of unifying therapeutic targets to overcome intra-tumor heterogeneity (validated in vivo in a human cancer xenograft), and mitigation of cancer cell resistance and rewiring mechanisms (validated in vitro in a human and canine cancer model). Conclusions These proof-of-concept studies support the use of an integrative functional approach to personalized combination therapy prediction for the population of high-risk cancer patients lacking viable clinical options and without actionable DNA sequencing-based therapy.

  Study EGAS00001003564

• South Asia Rheumatic Heart Disease Genetics Network

  Rheumatic heart disease (RHD) is a major cause of cardiovascular death and disability amongst young adults. Building on a sample collection established for and genotyped in several earlier successful candidate gene studies, this case-control study investigated host genetic susceptibility to RHD using a genome-wide approach. The current dataset comprises over 850 individuals recruited in Uttar Pradesh, India, including cases of RHD based on echocardiographic diagnosis and controls recruited on the basis of normal echocardiograms. For this analysis all available samples were genotyped using the Illumina HumanCore-24 BeadChip platform. Efforts are now underway to establish new sample collections in other areas of India as well as Pakistan and later other neighbouring countries.

  Study EGAS00001003565

• Genetic dysregulation of gene expression and splicing during a ten-year period of human aging in the PIVUS study

  We performed RNA sequencing in whole-blood from the same individuals from the PIVUS study at ages 70 and 80 to quantify how gene expression, alternative splicing, and their genetic regulation are altered during this 10-year period of advanced aging.

  Study EGAS00001003583

• GDAP___Genome_Diversity_in_Africa_Project

  The Genomic Diversity in Africa Project (GDAP) started with the plan to develop a genomic resource from African populations, characterise genomic diversity and population history, and facilitate clinical studies in Africa. Currently, 25 individuals from 24 ethnolinguistic groups have been whole-genome sequenced at high depth totalling 585 individuals. An additional 41 individuals have been sequenced with 10X Genomics libraries. At this stage, the initial curation of this dataset has been finished and we are performing the analysis in coordination with our collaborators. The current state of the GDAP represents a very diverse panel of African populations that maximizes geographical and ethnic variation and represents a great starting point to achieve the aforementioned goals. However, southern sub-Saharan countries, Bantu speakers and hunter gatherer groups are currently underrepresented, despite being crucial to understand the evolutionary history of the continent. After extensive effort to collate studies documentation, we finally have the opportunity to sequence 600 new individuals from these groups, including countries as Gabon, Rwanda and Zambia, and address these deficiencies. We aim to proceed with the same strategy: to sequence at high depth 25 individuals with standard PCR free libraries, with 2 additional individuals with 10X Genomics Chromium libraries per ethnolinguistic group. The former allows a good representation of variants down to low frequency in any given population, and the latter allows accurate phasing and the analysis of structural variation. By including these new populations, we want to investigate three crucial questions in African history in addition to the initial objectives: the Bantu expansion, the evolutionary history of hunter gatherers and the transatlantic slave trade. Additionally, the expanded dataset will help us better discover the genetic variation present in Africa and characterize the African pangenome. 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 EGAS00001003602

• Discovery and capture of novel dynamic DNA methylation in human sperm with preferential links to altered folate metabolism

  The sperm DNA methylation landscape is unique and critical for offspring health. If gamete-derived DNA methylation escapes reprogramming in early embryos, epigenetic defects in sperm may be transmitted to the next generation. Current techniques to assess sperm DNA methylation show bias towards CpG dense regions and do not target areas of dynamic methylation, those predicted to be environmentally sensitive and tunable regulatory elements. Thus, this study aim to assess variation in human sperm DNA methylation and design a targeted capture panel to interrogate the human sperm methylome.

  Study EGAS00001003617

• Recapitulation of genetic predisposition to medulloblastoma in human neuroepithelial stem cells

  Human neural stem cell cultures provide progenitor cells that are potential cells of origin for brain cancers. However, the extent to which genetic predisposition to tumor formation can be faithfully captured in stem cell lines is uncertain. Here we evaluated neuroepithelial stem (NES) cells, representative of cerebellar progenitors. We transduced NES cells with MYCN, observing medulloblastoma upon orthotopic implantation in mice. Significantly, transcriptomes and patterns of DNA methylation from xenograft tumors were globally more representative of human medulloblastoma compared to a MYCN-driven genetically engineered mouse model. Orthotopic transplantation of NES cells generated from Gorlin syndrome, patients predisposed to medulloblastoma due to germline mutated PTCH1, also generated medulloblastoma. We engineered candidate cooperating mutations in Gorlin NES cells, with mutation of DDX3X or loss of GSE1 both accelerating tumorigenesis. These findings demonstrate that human NES cells provide a potent experimental resource for dissecting genetic causation in medulloblastoma.

  Study EGAS00001003620

• Genome-wide search to find the genetic elements underlying visual contour perception

  Contour integration is a key function of the visual system to deal with occlusion and discontinuity in natural scenes. The current study is the first one to investigate the genetic basis of this visual function in humans. A total of 2619 normal participants were tested on their ability to detect continuous contours embedded in a cluttered background with a standard psychophysical protocol. Multi-level genomic analysis was performed, involving heritability estimation based on single nucleotide polymorphism (SNP) and association testing at SNP and gene levels.

  Study EGAS00001003639

• Pangenomic classification of pituitary adenomas

  Pituitary neuroendocrine tumors (PitNETs) are common. They are classified in five main histological types: PitNETs of POU1F1 (Pit1) lineage -lactotroph, somatotroph and thyrotroph -, gonadotroph, and corticotroph PitNETs. Pituitary tumorigenesis is poorly understood. We report here a pangenomic characterization of 134 PitNETs from all types, based on exome, genome, transcriptome, miRnome and methylome. Pit1 lineage was associated with a specific epigenetic signature of diffuse DNA hypomethylation, and correlated with expression of Pit1 markers and chromosomal instability. 45% of somatotroph PitNETs showed GNAS mutations. GNAS-wild type somatotroph showed a distinct transcriptome profile, expressing gonadotroph transcription factor NR5A1 (SF1). Sparsely granulated somatotroph PitNETs was classified apart from somatotroph, owing to a thyrotroph transcriptome signature. 43% of corticotroph PitNETs showed USP8 mutations, sharing a low epithelial-mesenchimal transition (EMT) transcriptome signature. Silent corticotroph PitNETs showed a gonadotroph transcriptome signature evoking transdifferentiation, and lost somatostatin target SSTR5 expression.

  Study EGAS00001003642

• GWAS Results from Danjou et al, Nature Genetics 2015

  We report genome-wide association study results for the levels of A1, A2 and fetal hemoglobins, analyzed for the first time concurrently. Integrating high-density array genotyping and whole-genome sequencing in a large general population cohort from Sardinia, we detected 23 associations at 10 loci. Five signals are due to variants at previously undetected loci: MPHOSPH9, PLTP-PCIF1, ZFPM1 (FOG1), NFIX and CCND3. Among the signals at known loci, ten are new lead variants and four are new independent signals. Half of all variants also showed pleiotropic associations with different hemoglobins, which further corroborated some of the detected associations and identified features of coordinated hemoglobin species production.

  Study EGAS00001003645

• Latency and interval therapy affect the evolution in metastatic colorectal cancer

  While comparison of primary tumor and metastases has highlighted genomic heterogeneity in colorectal cancer (CRC), previous studies have focused on a single metastatic site or limited genomic testing. Combining data from whole exome and ultra-deep targeted sequencing, we explored possible evolutionary trajectories beyond the status of these mutations, particularly among patient-matched metastatic tumors. Our findings confirm the persistence of known clinically-relevant mutations (e.g. those of RAS family of oncogenes) in CRC primary and metastases, yet reveal that latency and interval systemic therapy affect the course of evolutionary events within metastatic lesions. Specifically, our analysis of patient-matched primary and multiple metastatic lesions, developed over time, showed a similar genetic composition for liver metastatic tumors, which were 21-months apart. This genetic makeup was different from those identified in lung metastases developed before manifestation of the second liver metastasis. These results underscore the role of dormancy in the evolutionary path of metastatic CRC.

  Study EGAS00001003646

• Whole genome sequencing of six ethnic groups from Burkina Faso, Cameroon, and Tanzania

  We conducted low-coverage (~10x) sequencing of six ethnic groups from three malaria-endemic countries in sub-Saharan Africa, with the aim of characterising genetic diversity and improving imputation performance for genotyped samples in these populations. This data is part of a pre-publication release and is available under managed access. Please see https://www.malariagen.net/data/terms-use/human-gwas-data for more information on MalariaGEN data access policies. The purpose of the project was to support the discovery and understanding of genetic variants that influence human disease. Specifically defined goals were: a. the discovery of single nucleotide variants at frequencies of 1% or higher in diverse populations, b. even more comprehensive discovery (variants down to frequencies of 0.1 - 0.5%) in functional gene regions, and c. discovery of structural variants, such as copy number variants, other insertions and deletions, and inversions, including sequence-level understanding of breakpoints.

  Study EGAS00001003648

• multi-OMICs study of a pair of infant monozygotic twins with concordant B-cell ALL

  B-cell acute lymphoblastic leukemia (B-cell ALL) is the most common cancer in childhood. Studying identical twins with B-cell ALL provides a unique and tractable model for deciphering the developmental timing of pre- and post-natal mutations contributing to clonal evolution. To date, this has mainly focused on major cytogenetic subgroups of childhood B-cell ALL, including MLL fusions, ETV6-RUNX1, hyperdiploidy, and BCR-ABL1. However, formal demonstration of the prenatal origin and “backtracking” the natural history of the leukemia remains understudied in “B-other”/Normal Karyotype (NK) B-cell ALL. To characterize the genetic and the epigenetic landscape of this particular leukemia subtype, we performed whole genome DNA-, B-cell receptor (BCR)-, and DNA bisulfite-sequencing on a pair of 8-month-old monozygotic twins diagnosed with concordant “B-other”/NK B-cell ALL.

  Study EGAS00001003651

• Genetic Skin Disease Group

  Dac EGAC00001001126

• RNA-seq as a tool for evaluating human embryo competence

  The majority of embryos that are created through IVF do not implant. One possibility for this inefficiency is an incomplete understanding of the molecular factors affecting embryo competence. Currently, the process of selecting an embryo for uterine transfer utilizes an ad-hoc combination of morphological criteria, the kinetics of development, and genetic testing for aneuploidy. However, no single criterion can ensure selection of a viable embryo. In contrast, RNA-sequencing of embryos could yield highly dimensional data, which may provide additional insight and illuminate the discrepancies among current selection criteria. Indeed, recent advances enabling the production of RNA-sequencing (RNA-seq) libraries from single cells have facilitated the application of this technique to the study of some transcriptional events in early human development. However, these studies have not assessed the quality of their constituent embryos relative to commonly used embryological criteria. Here, we perform proof-of-principle advancement to clinical selection procedures by generating high quality RNA-seq libraries from a trophectoderm biopsy as well as the remaining whole embryo. We combine state-of-the-art embryological methods with low-input RNA-seq to develop the first transcriptome-wide approach for use in future predictive embryology studies. Specifically, we demonstrate the capacity of RNA-seq as a promising tool in preimplantation screening by showing that biopsies of an embryo can capture valuable information content available in the whole embryo from which they are derived. Furthermore, we show that this technique can be used to generate a RNA-based digital karyotype, and to develop a foundational dataset for identifying candidate competence-associated genes. Together, these data establish the foundation for a future RNA-based diagnostic in IVF.

  Study EGAS00001003667

• Control iPSC lines with clinically annotated genetic variants for versatile multi-lineage differentiation

  Induced Pluripotent Stem Cells (iPSC) derived from healthy individuals are important controls for disease modeling studies. To create a resource of genetically annotated iPSC, we reprogrammed footprint-free lines from four volunteers in the Personal Genome Project Canada (PGPC). Multilineage directed differentiation efficiently produced functional cortical neurons, cardiomyocytes and hepatocytes. Pilot users further demonstrated line versatility by generating kidney organoids, T-cells and sensory neurons. A frameshift knockout was introduced into MYBPC3 and these cardiomyocytes exhibited a hypertrophic phenotype as expected. WGS annotation of PGPC lines revealed on average 20 coding variants. Importantly, nearly all annotated PGPC and HipSci lines harboured at least one pre-existing or acquired variant with cardiac, neurological or other disease associations. Overall, PGPC lines were efficiently differentiated by multiple users into cell types found in six tissues for disease modelling, and clinical annotation highlighted variant-preferred lines for use as unaffected controls in specific disease settings.

  Study EGAS00001003684

• UAMS Smoldering Myeloma Timeline Cohort

  To better understand the pattern of genetic changes over time, we performed whole exome sequencing of sequential bone marrow samples from 9 patients taken overtime including some paired SMM/newly diagnosed MM/Relapse MM samples. Samples from 9 patients (9 controls and 53 tumors) underwent whole exome sequencing with an additional capture for the IGH, IHK, IGL, and MYC loci. DNA was obtained from either CD138+ cells from the bone marrow of smoldering myeloma patients through time (tumor) or from stem cell harvests or peripheral blood cells from the same patient (control). 100 ng of DNA was fragmented, end-repaired, and adapters ligated using NimbleGen's MedExome. After PCR amplification hybridized libraries underwent further amplification before being sequenced on a NextSeq500 (Illumina) using 75 bp paired end reads.

  Study EGAS00001003687

• Infant Glioma Molecular Subtype

  Infant gliomas have completely different clinical behavior than gliomas in children and adults. Low-grade tumors are more aggressive with a higher mortality rate, while high-grade tumors are more benign. However, we have little understanding of their biology and so cannot explain this behavior and what constitutes optimal clinical management. To solve this enigma, we performed a comprehensive genetic analysis of the largest cohort to date of infantile gliomas. We uncovered several important findings which dramatically change the management of these cancers. Infant hemispheric gliomas harbor unique alterations in ROS, ALK, MET and other genes found in adult carcinomas. Unlike in adults, they are typically the only alteration present and some of these tumors regress from high- to low-grade partially explaining their favorable clinical outcome. Midline gliomas harbor RAS/MAPK pathway mutations, progress rapidly, and do not respond to current chemoradiation approaches, partially explaining their poorer outcome. These data suggest that infant gliomas require different therapeutic approaches including early use of targeted therapies.

  Study EGAS00001003714

• Human adipose tissue immune cells

  We describe a novel and conserved Trem2+ lipid-associated macrophage (LAM) subset and identify markers, spatial localization, origin, and functional pathways associated with these cells. Genetic ablation of Trem2 in mice globally inhibits the downstream molecular LAM program, leading to adipocyte hypertrophy as well as systemic hypercholesterolemia, body fat accumulation, and glucose intolerance. These findings identify Trem2 signaling as a major pathway by which macrophages respond to loss of tissue-level lipid homeostasis, highlighting Trem2 as a key sensor of metabolic pathologies across multiple tissues and a potential therapeutic target in metabolic diseases.

  Study EGAS00001003725

• From matrimonial practices to genetic diversity in Southeast Asian populations: the signature of the matrilineal puzzle

  This project aims at revisiting the "matrilineal puzzle" from a population genetics' point of view. The dataset consists of 719 individuals from 12 populations sampled in continental Southeast Asia and genotyped at 598764 loci. This dataset can only be accessed and used for non-commercial research purposes with a finality complying with the informed consent provided by donors. For accessing the data, please provide: 1. An explicit title for the project you wish to use the data for. 2. A detailed institutional affiliation of all recipients and users of the data. 3. A summary (500 words) of the research project, with emphasis on its scientific and societal finality, and predicted outputs. If access is granted by the DAC, data cannot be shared with a third party not initially involved in the project, without the third party officially applying for accessing the data.

  Study EGAS00001003727

• DNA Methylation loss coupled with mitotic cell division promotes immune evasion of tumours with high mutation load

  Mitotic cell division increases tumour mutation burden and copy number load with a positive and inverse correlation, respectively, with the clinical benefit of immunotherapy. Markers of cell division correlate also with genomic demethylation involving methylation loss in late-replicating partial methylation domains. Here we find that immunomodulatory pathway genes are concentrated in these domains and transcriptionally repressed in demethylated tumours with CpG island promoter hypermethylation. Global methylation loss correlated with immune evasion signatures independently of mutation burden and aneuploidy. Methylome data of our cohort (n = 60) and a published cohort (n = 81) in lung cancer and a melanoma cohort (n = 40) consistently demonstrated that genomic methylation alterations counteract the contribution of high mutation burden and increase immunotherapeutic resistance. Higher predictive power was observed for methylation loss than mutation burden. We also found that genomic hypomethylation correlates with the immune escape signatures of aneuploid tumours. Hence, DNA methylation alterations implicate epigenetic modulation as a combination regimens for precision immunotherapy.

  Study EGAS00001003731

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

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

  Study EGAS00001003737

• Glioblastoma initiating cells are sensitive to histone demethylase inhibition due to epigenetic deregulation

  Tumor-initiating cells are a subpopulation of cells that have self-renewal capacity to regenerate a tumor. Here, we identify stem cell-like chromatin features in human glioblastoma initiating cells (GICs) and link them to a loss of the repressive histone H3 lysine 9 trimethylation (H3K9me3) mark. Increasing H3K9me3 levels by histone demethylase inhibition led to cell death in GICs but not in their differentiated counterparts. The induction of apoptosis was accompanied by a loss of the activating H3 lysine 9 acetylation (H3K9ac) modification and accumulation of DNA damage and downregulation of DDR genes. Upon knockdown of histone demethylases KDM4C and KDM7A both differentiation and DNA damage was induced. Thus, the H3K9me3-H3K9ac equilibrium is crucial for GIC viability and represents a chromatin feature that can be exploited to specifically target this tumor subpopulation.

  Study EGAS00001003750

• Genetics and transcriptomes of pediatric B cell precursor leukemia with gain of chromosome 21

  Down syndrome (DS, constitutive trisomy 21) children have a 27-fold increased risk of developing B-ALL (DS-ALL), face a worse outcome due to treatment-related morbidities, and an increased rate of relapses compared to other children. This highlights the need to better understand the mechanisms of DS-associated leukemogenesis to develop more adapted treatment. In this study, we characterized the genetic and transcriptomic landscapes of DS-ALL and found a high incidence of somatic mutations leading to RAS/MAPK pathway activation in DS-ALL, as seen in other pediatric B-ALL presenting somatic gains of the chromosome 21 (B-ALL+21).

  Study EGAS00001003760

• 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

• Patient-Derived Lung Cancer Organoid

  Lung cancer shows substantial genetic and phenotypic heterogeneity across individuals, driving a need for personalized medicine. Here we report lung cancer organoids and normal bronchial organoids established from patients tissues comprising five histological subtypes of lung cancer and non-neoplastic bronchial mucosa as in vitro models representing individual patient. The lung cancer organoids recapitulate the tissue architecture of the primary lung tumours and maintain the genomic alterations of the original tumours during long-term expansion in vitro. The normal bronchial organoids maintain cellular components of normal bronchial mucosa. Lung cancer organoids respond to drugs based on their genomic alterations: a BRCA2 mutant organoid to olaparib, an EGFR mutant organoid to erlotinib, and an EGFR mutant/MET amplified organoid to crizotinib. Considering the short length of time from organoid establishment to drug testing, our newly developed model may prove useful for predicting patient-specific drug responses through in vitro patient-specific drug trials.

  Study EGAS00001003786

• Methylation CYLD cutaneous syndrome

  Patients with CYLD cutaneous syndrome (CCS; syn. Brooke-Spiegler syndrome) carry germline mutations in the tumor suppressor CYLD and develop multiple skin tumors with diverse histophenotypes. We comprehensively profiled the genomic landscape of 42 benign and malignant tumors across 13 individuals from four multigenerational families. Novel recurrent mutations were found in epigenetic modifiers DNMT3A and BCOR in 29% of benign tumors. Multi-level and microdissected sampling strikingly reveal that many clones with different DNMT3A mutations exist in these benign tumors, suggesting that intra-tumor heterogeneity is common. Integrated genomic, methylation and transcriptomic profiling in selected tumours suggest that isoform specific DNMT3A2 mutations are associated with dysregulated Wnt/ß-catenin pathway signaling. Phylogenetic and mutational signature analyses confirm cylindroma pulmonary metastases from primary skin tumours. These findings add novel dimensions to existing paradigms of cutaneous tumorigenesis and metastasis.

  Study EGAS00001003800

• UPF - Evolutionary Population Genetics lab Data Access Committee

  Dac EGAC00001000244

• Whole-exome sequencing of Fanconi anemia-like inherited bone marrow failure syndrome

  We identified three fibroblast cultures derived from aplastic anemia patients in the Japanese Cancer Research Resources Bank (JCRB) Cell Bank that showed high levels of sister chromatid exchanges (SCEs) in PHA-stimulated lymphoblasts. To elucidate the genetic basis of the high SCE, we carried out whole exome sequencing of these samples. Subsequently, germline variants were analysed for these patients.

  Study EGAS00001003809

• Epigenetic and metabolomic data from type 2 diabetes adolescents

  Objective: Rates of type 2 diabetes (T2D) among adolescents are on the rise. Epigenetic changes could be associated with the metabolic alterations in adolescents with T2D. Methods: We performed a cross-sectional integrated analysis of DNA methylation data from peripheral blood mononuclear cells with serum metabolomic data from First Nation adolescents with T2D and controls participating in the Improving Renal Complications in Adolescents with type 2 diabetes through REsearch (iCARE) cohort study, to explore the molecular changes in adolescents with T2D. Results: Our analysis showed that 43 serum metabolites and 36 differentially methylated regions (DMR) were associated with T2D. Several DMRs were located near the transcriptional start site of genes with established roles in metabolic disease and associated with altered serum metabolites (e.g. glucose, leucine, and gamma-glutamylisoleucine). These included the free fatty acid receptor-1 (FFAR1), upstream transcription factor-2 (USF2), and tumor necrosis factor-related protein-9 (C1QTNF9), among others. Conclusions: We identified DMRs and metabolites that merit further investigation to determine their significance in controlling gene expression and metabolism which could define T2D risk in adolescents.

  Study EGAS00001003816

• Whole exome sequencing CYLD cutaneous syndrome

  Patients with CYLD cutaneous syndrome (CCS; syn. Brooke-Spiegler syndrome) carry germline mutations in the tumor suppressor CYLD and develop multiple skin tumors with diverse histophenotypes. We comprehensively profiled the genomic landscape of 42 benign and malignant tumors across 13 individuals from four multigenerational families. Novel recurrent mutations were found in epigenetic modifiers DNMT3A and BCOR in 29% of benign tumors. Multi-level and microdissected sampling strikingly reveal that many clones with different DNMT3A mutations exist in these benign tumors, suggesting that intra-tumor heterogeneity is common. Integrated genomic, methylation and transcriptomic profiling in selected tumours suggest that isoform specific DNMT3A2 mutations are associated with dysregulated Wnt/ß-catenin pathway signaling. Phylogenetic and mutational signature analyses confirm cylindroma pulmonary metastases from primary skin tumours. These findings add novel dimensions to existing paradigms of cutaneous tumorigenesis and metastasis.

  Study EGAS00001003839

• Targeted sequencing CYLD cutaneous syndrome

  Patients with CYLD cutaneous syndrome (CCS; syn. Brooke-Spiegler syndrome) carry germline mutations in the tumor suppressor CYLD and develop multiple skin tumors with diverse histophenotypes. We comprehensively profiled the genomic landscape of 42 benign and malignant tumors across 13 individuals from four multigenerational families. Novel recurrent mutations were found in epigenetic modifiers DNMT3A and BCOR in 29% of benign tumors. Multi-level and microdissected sampling strikingly reveal that many clones with different DNMT3A mutations exist in these benign tumors, suggesting that intra-tumor heterogeneity is common. Integrated genomic, methylation and transcriptomic profiling in selected tumours suggest that isoform specific DNMT3A2 mutations are associated with dysregulated Wnt/ß-catenin pathway signaling. Phylogenetic and mutational signature analyses confirm cylindroma pulmonary metastases from primary skin tumours. These findings add novel dimensions to existing paradigms of cutaneous tumorigenesis and metastasis.

  Study EGAS00001003840

• RNA sequencing CYLD cutaneous syndrome

  Patients with CYLD cutaneous syndrome (CCS; syn. Brooke-Spiegler syndrome) carry germline mutations in the tumor suppressor CYLD and develop multiple skin tumors with diverse histophenotypes. We comprehensively profiled the genomic landscape of 42 benign and malignant tumors across 13 individuals from four multigenerational families. Novel recurrent mutations were found in epigenetic modifiers DNMT3A and BCOR in 29% of benign tumors. Multi-level and microdissected sampling strikingly reveal that many clones with different DNMT3A mutations exist in these benign tumors, suggesting that intra-tumor heterogeneity is common. Integrated genomic, methylation and transcriptomic profiling in selected tumours suggest that isoform specific DNMT3A2 mutations are associated with dysregulated Wnt/ß-catenin pathway signaling. Phylogenetic and mutational signature analyses confirm cylindroma pulmonary metastases from primary skin tumours. These findings add novel dimensions to existing paradigms of cutaneous tumorigenesis and metastasis.

  Study EGAS00001003841

• English Longitudinal Study for Ageing (ELSA) Genetic Data Access

  Dac EGAC00001000270

• Discovering genetic causes of optic atrophy syndromes through whole exome sequencing

  The genetic defects leading to optic atrophy range from mitochondrial DNA (mtDNA) point mutations in Leber’s hereditary optic neuropathy (LHON), to dominant and recessive mutations affecting a cluster of nuclear genes implicated in mitochondrial dynamics. We performed WES in patients with an optic atrophy spectrum disorder, including retinal macular dystrophy and kidney insufficiency leading to transplantation, associated with mitochondrial DNA (mtDNA) depletion without accumulation of multiple deletions, to identify the genetic causes of this syndrome.

  Study EGAS00001003850

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

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

  Study EGAS00001003861

• VIKING Health Study - Shetland

  The Viking Health Study – Shetland, together with ORCADES collectively termed VIKING, aims to identify genes influencing risk factors for common diseases such as heart disease, diabetes, glaucoma and stroke. The genetic research cohort study in Shetland was established to identify genetic variants, both common and rare, increasing risk of disease. This isolated population has a number of characteristics, the most important being the very large number of distant relatives, which are favourable for the identification of rare variants influencing disease risk. This study describes PCR-free paired end whole genome sequencing (TruSeq DNA PCR-Free, Illumina) run on a HiSeqX platform at 30X coverage on 500 individuals from VIKING. Unrelated individuals from the largest families were selected first, followed by those from smaller families, until eventually related individuals were selected to best represent the variation in the full cohort.

  Study EGAS00001003872

• Bioscientia Molecular Genetics DAC

  Dac EGAC00001000271

• Systematic comparative analysis of single-nucleotide variants detection methods from single-cell RNA sequencing data

  Systematic interrogation of single nucleotide variations (SNVs) is one of the most promising approaches to delineate the cellular heterogeneity and phylogenetic relationships at the single cell level. While SNV detection from abundant single cell RNA sequencing (scRNA-seq) data is applicable and cost-effective in identifying expressed variants, inferring sub-clones, and deciphering genotype-phenotype linkages, there is a lack of computational methods specifically developed for SNV calling in scRNA-seq. Although variant callers for bulk RNA-seq have been sporadically used in scRNA-seq, the performances of different tools have not been assessed. Here, we performed a systematic comparison of seven tools including SAMtools, the GATK pipeline, CTAT, FreeBayes, MuTect2, Strelka2 and VarScan2, using both simulation and scRNA-seq datasets, and identified multiple elements influencing their performance. Our study provided the first benchmarking to evaluate the performances of different SNV detection tools for scRNA-seq data.

  Study EGAS00001003883

• A Unifying Paradigm for Transcriptional Heterogeneity and Squamous Features in Pancreatic Ductal Adenocarcinoma

  Pancreatic cancer expression profiles largely reflect a classical or basal-like phenotype. The extent to which these profiles vary within a patient is unknown. We integrated evolutionary analysis and expression profiling in multiregion sampled metastatic pancreatic cancers, finding that squamous features are the histologic correlate of an RNA-seq defined basal-like subtype. In patients with coexisting basal/squamous and classical/glandular morphology, phylogenetic studies revealed that squamous morphology represented a subclonal population in an otherwise classical/glandular tumor. Cancers with squamous features were significantly more likely to have clonal mutations in chromatin modifiers, intercellular heterogeneity for MYC amplification, and entosis. These data provide a unifying paradigm for integrating basal-type expression profiles, squamous histology, and somatic mutations in chromatin modifier genes in the context of clonal evolution of pancreatic cancer.  

  Study EGAS00001003974

• Sequencing data of multiple sarcoma samples for personalized medicine and endotype identification

  Sequencing data derived from multiple patient samples was used to guide multiple patient-specific experiments and identify potential therapeutic options. These diseases are rare and high-risk (undifferentiated pleomorphic sarcoma, rhabdomyosarcoma) and thus have few available validated clinical options. This leaves the population of patients with few options if, and when, standard therapy fails to cure or control disease. An aspect of the personalize therapy design included the identification of broader disease endotypes (a combination of genetic and clinical data describing patient status) to identify subgroups of the individual disease that could better define treatments and could segment patient populations for deeper preclinical research.

  Study EGAS00001003981

• Lindsey E. Jones et al. Patient-derived cells from recurrent tumors that model the evolution of IDH-mutant glioma. Neuro-Oncology Advances, 2020. We established robust, versatile, and well-characterized Patient Derived Cells (PDC) of an IDH1-mutant astrocytoma and an IDH1-mutant oligodendroglioma that represents defined evolutionary stages of chemotherapy-induced hypermutation. The PDCs retain tumor subtype defining features over time as well as classic hallmarks of cancer, including anchorage independent growth and cell immortality. The integrated phylogenies composed of PDCs,single-cell-derived PDCs, patient-derived xenografts, and corresponding spatiotemporal tumor tissue samples also provide new insight into the clonality, evolutionary pattern, and immense mutational load of hypermutated IDH mutant gliomas.

  IDH-mutant lower-grade gliomas can undergo malignant progression via temozolomide-driven hypermutation. Patient-derived cells (PDC) that model the genetically distinct hypermutated (HM) tumor subgroup are generally lacking, and few if any human brain tumor cell models are from defined evolutionary time points. Here, we characterize multiple PDC derived from independent surgical specimens of IDH1-mutant recurrences, including an ATRX and TP53-mutant astrocytoma and a 1p/19q co-deleted and TERT promoter-mutant oligodendroglioma. We determined the evolutionary time points represented by each PDC using exome sequencing and phylogenetic reconstruction, comparing the PDC and single cell clones of the PDC (scPDC) to multiple spatiotemporal tumor tissue samples, PDC-derived xenografts (PDX) and patient-matched blood. The tumor samples exhibited TMZ-induced mutagenesis and a branching pattern of evolution. We found clear evidence of two fully independent founder HM clones in the tumor tissue that are faithfully represented by independent PDC. The PDC, scPDC and PDX also shared the mutagenesis signature and represent the mid and later evolutionary time points of their corresponding tumors. The PDC maintained the tumor subtype-defining features over many passages, including heterozygosity of the IDH1 R132H mutation, production of 2-hydroxyglutarate (2-HG), and subtype-specific telomere maintenance mechanisms. The PDC from both tumors exhibited anchorage-independent growth in soft agar. The oligodendroglioma PDC formed infiltrative intracranial tumors with characteristic oligodendroglioma histology, initially with a long period to tumor formation. We conclude that the PDC, scPDC and PDX faithfully model the heterogeneous clonal origins of the corresponding tumor tissue. The multilevel analysis also provides new insight into the intratumoral heterogeneity and vast mutational load of HM glioma. The PDC from multiple evolutionary time points presented in the context of full clinical timelines may be useful to model evolution and intratumoral heterogeneity, important sources of therapeutic failure.

  Study EGAS00001003992

• Mapping the breast cancer metastatic cascade onto circulating tumour DNA using genetic and epigenetic clonal tracking

  Circulating tumour DNA (ctDNA) allows tracking the evolution of human cancers at high resolution, overcoming many limitations of tissue biopsies. However, exploiting ctDNA to determine how a patient’s cancer is evolving in order to aid clinical decisions remains difficult. This is because ctDNA is a mix of fragmented alleles, and the contribution of different cancer cell subpopulations from distinct lesions to ctDNA remains largely unknown. Profiling of ctDNA almost invariably requires prior knowledge of what genomic alterations to track. Here we leverage on a rapid autopsy programme to demonstrate that unbiased genomic characterisation of several metastatic sites and concomitant ctDNA profiling at whole-genome resolution reveals the extent to which ctDNA is representative of widespread disease. We also present a ctDNA methylation profiling method that allows tracking evolutionary changes in the cancer cell population at single-molecule resolution without prior knowledge. These results have critical implications for the use of liquid biopsies to monitor cancer evolution in humans and potentially guide treatment.

  Study EGAS00001004014

• Uncovering tumor intrinsic and extrinsic factors that regulate hepatocellular carcinoma growth using patient derived xenograft assays

  Several major factors limit our understanding of hepatocellular carcinoma (HCC). First, human HCCs are infrequently biopsied for diagnosis and thus are not often biologically interrogated. Second, HCC initiation and progression are strongly influenced by the cirrhotic microenvironment, and the exact contributions of intrinsic and extrinsic tumor factors are unclear. A powerful approach to examine the personalized biology of liver cancers and the influence of host tissues is with patient derived xenograft (PDX) models. In Asia, HCCs from hepatitis B virus patients have been efficiently converted into PDXs, but few parallel efforts from the west have been reported. In a large-scale analysis, we implanted 93 HCCs and 8 cholangiocarcinomas (CCAs) to systematically analyze host factors and to define an optimized platform for PDX development from both surgical and biopsy samples. NSG mice that had undergone partial hepatectomy (PHx) represented the best combination of engraftability, growth, and passageability, but overall rates were low and indicative of a unique intrinsic biology for HCCs in the US. PDX models preserved the histology and genetic features of parental tumors, and ultimately, 9 new models were usable for pre-clinical studies. Intriguingly, HCC PDXs were differentially sensitive to regorafenib and sorafenib and CCA PDXs were also highly sensitive to regorafenib. PDX models functionalize early and advanced stage HCCs and revealed unique biological features of liver cancers from the US.

  Study EGAS00001004020

• Exome Sequencing of Spanish Patients with rare genetic diseases.

  Whole exome sequencing of Spanish patients suffering from a rare genetic disease. The study was carried out in 2013 as a part of a public call from CNAG (Centro Nacional de Análisis Genómico), where data of eight spanish families were sent to analyse. The study concluded with three solved cases: two of aniridia and one syndromic RP.

  Study EGAS00001004028

• Human Pancreatic Islet RNAseq - Lund

  Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, many key tissues and cell-types required for appropriate functional inference are absent from large-scale resources such as ENCODE and GTEx. We explored the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using RNA-Seq and genotyping data from 420 islet donors. We find: (a) eQTLs have a variable replication rate across the 44 GTEx tissues (<73%), indicating that our study captured islet-specific cis-eQTL signals; (b) islet eQTL signals show marked overlap with islet epigenome annotation, though eQTL effect size is reduced in the stretch enhancers most strongly implicated in GWAS signal location; (c) selective enrichment of islet eQTL overlap with the subset of T2D variants implicated in islet dysfunction; and (d) colocalization between islet eQTLs and variants influencing T2D or related glycemic traits, delivering candidate effector transcripts at 23 loci, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in tissues of greatest disease-relevance while expanding our mechanistic insights into complex traits association loci activity with an expanded list of putative transcripts implicated in T2D development.

  Study EGAS00001004042

• GWAS data on human pancreatic islets from 191 donors - Lund

  Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, many key tissues and cell-types required for appropriate functional inference are absent from large-scale resources such as ENCODE and GTEx. We explored the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using RNA-Seq and genotyping data from 420 islet donors. We find: (a) eQTLs have a variable replication rate across the 44 GTEx tissues (<73%), indicating that our study captured islet-specific cis-eQTL signals; (b) islet eQTL signals show marked overlap with islet epigenome annotation, though eQTL effect size is reduced in the stretch enhancers most strongly implicated in GWAS signal location; (c) selective enrichment of islet eQTL overlap with the subset of T2D variants implicated in islet dysfunction; and (d) colocalization between islet eQTLs and variants influencing T2D or related glycemic traits, delivering candidate effector transcripts at 23 loci, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in tissues of greatest disease-relevance while expanding our mechanistic insights into complex traits association loci activity with an expanded list of putative transcripts implicated in T2D development.

  Study EGAS00001004044

• Phenotyping data on human pancreatic islets from 191 donors - Lund

  Most signals detected by genome-wide association studies map to non-coding sequence and their tissue-specific effects influence transcriptional regulation. However, many key tissues and cell-types required for appropriate functional inference are absent from large-scale resources such as ENCODE and GTEx. We explored the relationship between genetic variants influencing predisposition to type 2 diabetes (T2D) and related glycemic traits, and human pancreatic islet transcription using RNA-Seq and genotyping data from 420 islet donors. We find: (a) eQTLs have a variable replication rate across the 44 GTEx tissues (<73%), indicating that our study captured islet-specific cis-eQTL signals; (b) islet eQTL signals show marked overlap with islet epigenome annotation, though eQTL effect size is reduced in the stretch enhancers most strongly implicated in GWAS signal location; (c) selective enrichment of islet eQTL overlap with the subset of T2D variants implicated in islet dysfunction; and (d) colocalization between islet eQTLs and variants influencing T2D or related glycemic traits, delivering candidate effector transcripts at 23 loci, including DGKB and TCF7L2. Our findings illustrate the advantages of performing functional and regulatory studies in tissues of greatest disease-relevance while expanding our mechanistic insights into complex traits association loci activity with an expanded list of putative transcripts implicated in T2D development.

  Study EGAS00001004056

• Expression profiles and genetic makeup of metastases of a cancer of unknown primary.

  Study EGAS00001004059

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

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

  Study EGAS00001004063

• RNA-seq study of longitudinal blood cell samples from children at risk of type 1 diabetes

  The appearance of type 1 diabetes (T1D)-associated autoantibodies is the first and only measurable parameter to predict progression toward T1D in genetically susceptible individuals. However, autoantibodies indicate an active autoimmune reaction, wherein the immune tolerance is already broken. Therefore, there is a clear and urgent need for new biomarkers that predict the onset of the autoimmune reaction preceding auto-antibody positivity or reflect progressive b-cell destruction. Here we report the mRNA sequencing–based analysis of 306 samples including fractionated samples of CD4+ and CD8+ T cells as well as CD4, CD8 cell fractions and unfractionated PBMC samples longitudinally collected from seven children who developed beta-cell autoimmunity (Case subjects) at a young age and matched control subjects.

  Study EGAS00001004071

• Evolutionary Origins of Recurrent Pancreatic Cancer

  Surgery is the only curative option for patients with Stage I/II pancreatic cancer, nonetheless most patients will recur after surgery and die of their disease. To identify novel opportunities for management of recurrent pancreatic cancer we performed whole exome or targeted sequencing of 10 resected primary cancers and their matched intrapancreatic recurrences or distant metastases. We identified that recurrent disease develops from a monophyletic or polyphyletic origin, and this distinction has potential clinical relevance. In all patients, treatment induced genetic bottlenecks led to a modified genetic landscape and subclonal heterogeneity for driver gene alterations in part due to intermetastatic seeding. In one patient what was believed to be recurrent disease was an independent (second) primary tumor. These findings strongly advocate for routine post-treatment sampling in the management of recurrent pancreatic cancer.

  Study EGAS00001004097

• Comprehensive pharmacogenomic characterization of gastric cancer

  Gastric cancer is among the most lethal human malignancies. Previous studies have identified molecular aberrations that constitute dynamic biological networks and genomic complexities of gastric tumors. However, the clinical translation of molecular-guided targeted therapy is hampered by challenges. Notably, solid tumors often harbor multiple genetic alterations, complicating the development of effective treatments. To address such challenges, we established a comprehensive dataset of molecularly annotated patient derivatives coupled with pharmacological profiles for 60 targeted agents to explore dynamic pharmacogenomic interactions in gastric cancers. We identified lineage-specific drug sensitivities based on histopathological and molecular subclassification, including substantial sensitivities toward VEGFR and EGFR inhibition therapies in diffuse- and signet ring-type gastric tumors, respectively. We identified potential therapeutic opportunities for WNT pathway inhibitors in ALK-mutant tumors, a significant association between PIK3CA-E542K mutation and AZD5363 response, and transcriptome expression of RNF11 as a potential predictor of response to gefitinib. Collectively, our results demonstrate the feasibility of drug screening combined with tumor molecular characterization to facilitate personalized therapeutic regimens for gastric tumors.

  Study EGAS00001004106

• Reference Exome Data for a Northern Brazilian population

  Whole exome sequencing (WES) is a popular technology widely used in the diagnosis of rare genetic diseases and provides useful variant data for analysis of complex diseases. However, there is not always adequate population-specific reference data available to assist clinical geneticists in assigning a diagnostic variant to a specific clinical condition. Here we provide a catalogue of variants called after sequencing the exomes of 45 babies from the State of Rio Grande do Nord in Brazil. Sequence data were processed using an ‘intersect-then-combine’ (ITC) approach, using GATK and SAMtools to call variants. A total of 612,761 variants were identified in at least one individual in this Brazilian Cohort, including 559,448 single nucleotide variants (SNVs) and 53,313 insertion/deletions. Of these, 58,111 overlapped with nonsynonymous (nsSNVs) or splice site (ssSNVs) SNVs in dbNSFP. As an aid to clinical diagnosis of rare diseases, we used the American College of Medicine Genetics and Genomics (ACMG) guidelines to assign pathogenic/likely pathogenic status to 439 (0.755%) of the 58,111 nsSNVs and ssSNVs. Our data set provides a useful reference point for diagnosis of rare diseases in Brazil.

  Study EGAS00001004112

• CD8+ T-cell exhaustion induced by leukemic cells drives progression in Chronic Lymphocytic Leukemia

  In order to shed light on mechanisms that lead to CLL progression, we investigated longitudinal changes in both the genetic and immunological landscapes. For that, we analyzed paired samples obtained from CLL patients at diagnosis and progression before front-line treatment or at diagnosis and long-term asymptomatic follow-up. By whole-exome sequencing (WES), we found that CLL cells displayed limited genetic evolution at progression. Conversely, T cells from progressing patients showed significant immunophenotypical and transcriptional changes over time not observed in patients that did not progress. Specifically, we found that CD8+ T cells at progression were characterized by an increasing terminally exhausted effector phenotype (T-betdim/- Eomeshi PD1hi) with high co-expression of inhibitory receptors (PD1, CD244 and CD160). This exhausted status in CD8+ T cells was induced by progressing leukemic cells mainly through a mechanism dependent on soluble factors, including IL-10. Thus, we demonstrate that immunological changes are of paramount importance for CLL clinical progression, thereby providing a rationale for the use of early immunotherapeutic intervention in this disease.

  Study EGAS00001004116

• Genomic landscape of inflammatory breast cancer by whole-genome sequencing

  Inflammatory breast cancer (IBC) is the most aggressive form of breast cancer, characterized by a highly invasive and metastatic phenotype, but little is known about its genetic drivers. Here we report whole-genome sequencing (WGS) of 20 newly diagnosed IBC biopsies to characterize for the first time the entire genomic landscape of this disease.

  Study EGAS00001004117