23249 results for "ega"

in 81.11 milliseconds.

• UK10K_RARE_SIR

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes.Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches.The SIR (Severe Insulin Resistance) samples are part of the “rare disease” group, and will undergo exome sequencing. For further information with regard to this cohort please contact Robert Semple (rks16@cam.ac.uk).

  Study EGAS00001000130

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

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

  Study EGAS00001004290

• Stratifying and Targeting Pediatric Medulloblastoma through Genomics

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

  Study EGAS00001000273

• Clonal_haematopoiesis_in_patients_with_AAA

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

  Study EGAS00001002873

• UK10K RARE CHD

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The CHD (congenital heart disease) samples will be part of the “rare disease” group, and will undergo exome sequencing.For further information with regard to this cohort please contact Shoumo Bhattacharya (shoumo@me.com).

  Study EGAS00001000125

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

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

  Study EGAS00001001028

• P647_Targeted_resequencing_project

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

  Study EGAS00001000305

• Mutational Patterns in Metastatic Cutaneous Squamous Cell Carcinoma

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

  Study EGAS00001003370

• Exome sequencing of a novel cervical cancer cell line

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

  Study EGAS00001003343

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

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

  Study EGAS00001003372

• Intratumoural Heterogeneity Underlies Distinct Therapy Responses and Treatment Resistance in Glioblastoma

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

  Study EGAS00001003438

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Colorectal_LCM

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

  Study EGAS00001003455

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

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

  Study EGAS00001003443

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

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

  Study EGAS00001004327

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

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

  Study EGAS00001004343

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

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

  Study EGAS00001001240

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

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

  Study EGAS00001001810

• A somatic reference standard for cancer genome sequencing with COLO829

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

  Study EGAS00001001385

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

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

  Study EGAS00001001515

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

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

  Study EGAS00001001690

• Ethiopia_Genome_Project__low_coverage_

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

  Study EGAS00001000238

• UK10K OBESITY TWINSUK

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. These TwinsUK samples are a subset of the TwinsUK cohort samples with a BMI>=40 that will be exome sequenced by the Obesity group.

  Study EGAS00001000306

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

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

  Study EGAS00001003689

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

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

  Study EGAS00001002654

• The Prediction and Prevention of Preeclampsia

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

  Study EGAS00001001898

• The_Genomic_Advances_in_Sepsis__GAinS__RNA_seq

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

  Study EGAS00001003772

• Extreme phenotypes define epigenetic and metabolic signatures in cardiovascular diseases.

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

  Study EGAS00001003780

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

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

  Study EGAS00001003860

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

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

  Study EGAS00001002406

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

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

  Study EGAS00001002700

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

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

  Study EGAS00001005053

• THE GENOMIC LANDSCAPE OF ACTINIC KERATOSIS

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

  Study EGAS00001004243

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

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

  Study EGAS00001004291

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

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

  Study EGAS00001002875

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

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

  Study EGAS00001004018

• Adult-type Granulosa Cell Tumour of the Ovary

  A recurrent mutation in FOXL2 (c.402C>G; p.C134W) is present in over 95% of adult-type granulosa cell tumours (AGCTs). In contrast, various loss-of-function mutations in FOXL2 lead to the development of blepharophimosis, ptosis and epicanthus inversus syndrome (BPES). BPES is characterized by an eyelid malformation often accompanied with primary ovarian insufficiency but not granulosa cell tumours. Two recent studies suggest that FOXL2 C402G is a gain or change-of-function mutation with altered DNA-binding specificity. Another study proposes that FOXL2 C402G is selectively targeted for degradation inducing somatic haploinsufficiency suggesting its role as a tumour suppressor. The latter study relies on data indicative of an FOXL2 allelic imbalance in AGCTs. Here we present RNA-seq data as genetic evidence that no real allelic imbalance is observed at the transcriptomic level in AGCTs. Additionally, there is no loss of protein expression in tumours harboring the mutated allele. These data and other features of this mutation as compared to other oncogenes and tumour suppressor genes argue strongly against FOXL2 being a tumour suppressor in this context. Given the likelihood that FOXL2 C402G is oncogenic, targeting the variant protein or its downstream consequences is the most viable path forward to identifying an effective treatment for this cancer.

  Study EGAS00001005414

• Hypoxia acts as an environmental cue for the human TRM differentiation program

  Tissue-resident memory T-cells (TRM) provide frontline defense against infectious diseases and contribute to anti-tumor immunity; however, aside from the necessity of TGF-β, knowledge regarding TRM-inductive cues remains incomplete, particularly for human cells. Oxygen tension is an environmental cue that distinguishes peripheral tissues from the circulation and here, we demonstrate that differentiation of human CD8+ T-cells in the presence of hypoxia and TGF-β1 led to the development of a TRM phenotype, characterized by a greater than five-fold increase in CD69+CD103+ cells expressing human TRM hallmarks and enrichment for endogenous human TRM gene signatures, including increased adhesion molecule expression and decreased expression of genes involved in recirculation. Hypoxia and TGF-β1 synergized to produce a significantly larger population of TRM phenotype cells than either condition alone, and comparison of these cells from the individual and combination conditions revealed distinct phenotypic and transcriptional profiles, indicating a programming response to milieu rather than a mere expansion. Our findings identify a previously unreported cue for the TRM differentiation program and can enable facile generation of human TRM phenotype cells in vitro for basic studies and translational applications such as adoptive cellular therapy.

  Study EGAS00001005286

• To determine the mutational impact of the in vitro culture, clonal human adult and pluripotent stem cell lines were subjected to a second clonal step after 3 months of culture. These subclones were whole genome sequenced to identify all the mutations that accumulated during the 3 month culture period.

  Genetic changes acquired during culture pose a potential risk for the successful application of stem cells. To assess the risk of in vitro expansion on mutation accumulation we have performed whole genome sequencing of clonally expanded human induced pluripotent stem cells (iPS cells) and adult stem cells (ASCs) to identify all mutations that accumulated over a fixed culture period. We find that ASCs acquire more single nucleotide variants (SNVs) and indels per population doubling than the iPS cells. When compared with ASCs, iPS cells are more vulnerable to mutations in genes and promoters. Mutational analysis revealed a clear in vitro induced mutational signature that is irrespective of stem cell type. This in vitro signature is characterized by C to G transversions that are probably caused by oxidative stress. Additionally, we observed stem cell specific mutational signatures and differences in transcriptional strand bias, indicating differential activity of DNA repair mechanisms between stem cell types in culture. In conclusion, in vitro culture induces mutation accumulation in iPS cells and ASCs. Culture under low-oxygen tension may help to reduce the number of culture-induced mutations.

  Study EGAS00001002955

• The Molecular Landscape of Asian Breast Cancers Reveals Clinically Relevant Population-Specific Differences

  Molecular profiles of breast cancer have contributed to an improved understanding of the disease, enabled development of molecular prognostic signatures to guide treatment decisions and unveiled new or more accurate therapeutic options for breast cancer patients. However, the extent to which differences in genetic, environmental and lifestyle factors influence molecular profiles in different populations remains poorly characterised, as relatively few large-scale molecular studies of breast tumours in non-Caucasian populations have hitherto been reported. Here, we present the molecular profiles of 560 Asian breast tumours and a comparative analysis of breast cancers arising in Asian and Caucasian women. Compared to the breast tumours in predominantly Caucasian women reported in TCGA and METABRIC, we show an increased prevalence of Her2-enriched molecular subtypes and higher prevalence of TP53 somatic mutations in ER+ Asian breast tumours. Using gene expression and immunohistochemistry, we observed elevated immune scores in Asian breast tumours, suggesting potential clinical response to immune checkpoint inhibitors. Whilst Her2-subtype and enriched immune score are associated with improved survival, presence of TP53 somatic mutations is associated with poorer survival in ER+ tumours. Taken together, these population differences unveil new opportunities to improve the understanding of this disease and lay the foundation for precision medicine in different populations.

  Study EGAS00001004518

• Host whole genome variations are associated with neurocognitive outcome in survivors of pediatric medulloblastoma

  Host whole genome analysis is a promising source of predictive information for long-term morbidity in cancer survivors. However, studies on genetic predictors of long-term outcome, particularly neurocognitive function following chemoradiation in pediatric oncology are limited. In the present study, we evaluated variations in host whole genome single nucleotide polymorphisms (SNPs) and their association with cognitive outcome. Whole-genome SNP analysis of host peripheral blood was conducted on 22 medulloblastoma long-term survivors, of whom 18 completed neuropsychological testing. First, unsupervised consensus clustering of the most variable SNPs within 409 genes involved in DNA repair was performed. Discrete variant groups were identified, although they were not associated with cognitive outcome, suggesting that variations in genes corresponding to a single functional group may be insufficient to predict long-term outcome alone. In support of this interpretation, unsupervised hierarchical clustering analysis using disease-associated gene variants by cognitive impairment status yielded two distinct variant clusters comprised of 36 variants, 34 of which were in noncoding regions. These findings illustrate for the first time that cognitively impaired survivors have a distinct variant profile compared to other medulloblastoma survivors. Future research in larger cohorts is needed to validate host genome predictors of cognitive impairment that may impact clinical management.

  Study EGAS00001002996

• Whole genome, exome and transcriptome analysis to guide individualized cancer interpretation

  Study Design: The Cancer Alliance study was a collaboration of nine academic medical institutions in the New York city area and led by the New York Genome Center. Tumor-normal paired samples were submitted by the institutions for whole genome, whole exome and transcriptome sequencing and bioinformatics analysis at New York Genome Center. Relevant clinical histories were also collected. Raw sequence files were available to all collaborators. Variant call files were analyzed both manually by interpreters at NYGC and by IBM Watson Genome Analytics. This study was approved by a central institutional review board (IRB), Biomedical Research Alliance of New York, and by local IRBs, including Stony Brook University and Northwell Health. The study was supported in part by a grant from the IBM corporation (IBM Watson Health) to the New York Genome Center, New York Genome Center philanthropic funds and Rockefeller University grant # UL1TR000043 from the National Center for Advancing Translational Sciences (NCATS), and the National Institutes of Health (NIH) Clinical and Translational Science Award (CTSA) program. The principal investigator Robert B. Darnell MD, PhD., is a Howard Hughes Medical Institute Investigator. Results: Pending Publication of Analysis of Study Conclusions: Pending Publication of Analysis of Study

  Study EGAS00001004013

• Hermansky-Pudlak syndrome type 1 causes impaired anti-microbial immunity through a pathogenic lipid metabolism-mTOR circuit - 10x Genomics scRNAseq

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

  Study EGAS00001005098

• Variation and transmission of the human gut microbiota across generations - shotgun data

  Although composition and functional potential of the human gut microbiota evolve over lifespan, kinship has been identified as a key covariate of microbial community diversification. To date, sharing of microbiota features within families has however mostly been assessed between parents and their direct offspring. Here, we investigate potential transmission and persistence of familial microbiome patterns and microbial genotypes in a family cohort (N=102) spanning three to five generations over the same female bloodline. We observe microbiome community composition to be associated with kinship, with seven (low-abundant) genera displaying familial distribution patterns. While kinship and current cohabitation emerged as closely entangled variables, our explorative analyses of microbial genotype distribution and transmission estimates point at the latter as a key covariate of strain dissemination. Highest potential transmission rates are estimated between sisters and mother-daughter pairs, decreasing with increasing daughter’s age, and being higher among cohabiting pairs than those living apart. Although rare, we do detect potential transmission events spanning three and four generations, primarily involving species of the genera Alistipes and Bacteroides. Overall, while our analyses confirm the existence of family-bound microbiome community profiles, transmission or co-acquisition of bacterial strains appears to be strongly linked to cohabitation.

  Study EGAS00001005649

• Variation and transmission of the human gut microbiota across generations - 16S data

  Although composition and functional potential of the human gut microbiota evolve over lifespan, kinship has been identified as a key covariate of microbial community diversification. To date, sharing of microbiota features within families has however mostly been assessed between parents and their direct offspring. Here, we investigate potential transmission and persistence of familial microbiome patterns and microbial genotypes in a family cohort (N=102) spanning three to five generations over the same female bloodline. We observe microbiome community composition to be associated with kinship, with seven (low-abundant) genera displaying familial distribution patterns. While kinship and current cohabitation emerged as closely entangled variables, our explorative analyses of microbial genotype distribution and transmission estimates point at the latter as a key covariate of strain dissemination. Highest potential transmission rates are estimated between sisters and mother-daughter pairs, decreasing with increasing daughter’s age, and being higher among cohabiting pairs than those living apart. Although rare, we do detect potential transmission events spanning three and four generations, primarily involving species of the genera Alistipes and Bacteroides. Overall, while our analyses confirm the existence of family-bound microbiome community profiles, transmission or co-acquisition of bacterial strains appears to be strongly linked to cohabitation.

  Study EGAS00001005651

• A recessive PRDM13 mutation results in congenital hypogonadotropic hypogonadism and cerebellar hypoplasia

  The PRDM13 (PR Domain containing 13) putative chromatin modifier and transcriptional regulator functions downstream of the transcription factor PTF1A, which controls GABAergic fate in the spinal cord and neurogenesis in the hypothalamus. Here, we report a novel, recessive syndrome associated with PRDM13 mutation. Patients exhibited intellectual disability, ataxia with cerebellar hypoplasia, scoliosis and delayed puberty with congenital hypogonadotropic hypogonadism (CHH). Expression studies revealed Prdm13/PRDM13 transcripts in the developing hypothalamus and cerebellum in mouse and human. An analysis of hypothalamus and cerebellum development in mice homozygous for a Prdm13 mutant allele revealed a significant reduction in the number of Kisspeptin (Kiss1) neurons in the hypothalamus and PAX2+ progenitors emerging from the cerebellar ventricular zone. The latter was accompanied by ectopic expression of the glutamatergic lineage marker TLX3. Prdm13-deficient mice displayed cerebellar hypoplasia, normal gonadal structure, but delayed pubertal onset. Together, these findings identify PRDM13 as a critical regulator of GABAergic cell fate in the cerebellum and of hypothalamic kisspeptin neuron development, providing a mechanistic explanation for the co-occurrence of CHH and cerebellar hypoplasia in this syndrome. To our knowledge, this is the first evidence linking disrupted PRDM13-mediated regulation of Kiss1 neurons to CHH in humans.

  Study EGAS00001005878

• Nanopore_and_Illumina_sequencing_of_human_glioblastomas

  Liquid biopsy analysis refers to methods designed to detect tumour-specific material (e.g., circulating tumour cells or tumour DNA) in body fluids, such as urine or blood samples. A widely-used liquid biopsy approach consists of genotyping the circulating tumour DNA (ctDNA) through sequencing of plasma/serum DNA. Although useful in the context of metastatic disease (where the concentration of ctDNA is high), current liquid biopsy technologies show limited sensitivity of detection for the early detection of cancer, and low specificity, as cancer-related mutations accumulate in healthy tissues as part of the ageing process, thus making it difficult to distinguish these from tumour mutations, and because sequencing errors and true mutations show overlapping profiles. Therefore, developing liquid biopsy protocols with increased sensitivity and specificity represents an urgent clinical need. Here we harness extrachromosomal circular DNA (eccDNA) elements, which are circular DNA structures physically separated from the chromosomes of up to several Mbp long pervasive in human cancers, for liquid biopsy analysis. In this pilot study we will focus on the analysis of glioblastomas, because there is strong evidence for the presence of eccDNA in these tumour types, and because developing liquid biopsy approaches for brain tumours to reduce the invasiveness of brain tumour biopsies remains an unmet clinical need.

  Study EGAS00001005812

• Sclerosing epithelioid fibrosarcoma case report

  Sclerosing epithelioid fibrosarcoma (SEF) is a rare and aggressive soft tissue sarcomathought to originate in fibroblasts of the tissues surrounding tendons, ligaments and muscles. Minimally responsive to conventional cytotoxic chemotherapies, greater than 50% of SEF patients experience relapse and/or metastatic disease. SEF is most commonly discovered in middle-aged and elderly adults, but also rarely in children. A common gene fusion occurring between the EWSR1 and CREB3L1 genes has been observed in 80-90% of SEF cases. We describe here the youngest SEF patient reported to date (a 3-year-old Caucasian male) who presented with numerous bony and lung metastases. Additionally, we perform a comprehensive literature review of all SEF-related articles published since the disease was first characterized. The patient described in this report experienced persistent disease progression despite being heavily treated with multiple surgeries, radiation, numerous chemotherapies and targeted therapeutics. Primary tumor, metastatic and relapse sites were sequenced by whole genome, whole exome, and deep transcriptome next generation sequencing with comparison to a patient-matched normal blood sample. Consistent across all sequencing analyses was the disease-defining EWSR1-CREB3L1 fusion as a single feature consensus. We provide an analysis of our genomic findings and a discussion as to potential therapeutic strategies for SEF.

  Study EGAS00001005215

• Comprehensive Genomic Characterization of Refractory Multiple Myeloma Reveals a Complex Mutational and Structural Landscape Associated with Drug Resistance (H067)

  Outcome of patients with multiple myeloma (MM) refractory to immunomodulatory agents (IMiDs) and proteasome inhibitors (PIs) is still very poor. To elucidate the genetic landscape of relapsed/refractory MM (RRMM) we performed whole genome and transcriptome sequencing of 39 heavily pretreated RRMM patients with at least double refractoriness revealing complex structural changes and a high mutational load (H067/K08K). We found a median of 67 structural variants and 116 nonsynonymous exonic mutations per patient. Assessment of oncogenic drivers showed NRAS, BRAF, TP53, KRAS, SLC4A7, FAM46C, DIS3, RB1, KMT2C, MLLT4, EWSR1, HCFC2, and COPS3 as significantly mutated. We found multiple genomic regions with bi-allelic events affecting known tumor suppressor genes, most frequently TP53, RB1, and TRAF3. Supervised analysis of mutational signatures showed a major contribution of Cosmic signature 6 associated with deficiency in DNA mismatch repair and the ‘BRCAness’ signature which coincided with mutations in genes known to confer PARP inhibitor sensitivity. Recurrently mutated gene networks included genes with relevance for PI and IMiD activity, the latter including in particular members of the Cereblon and the COP9 signalosome complex. In conclusion, this comprehensive genomic characterization revealed a complex mutational and structural landscape in RRMM and highlights potential implications for therapeutic strategies.

  Study EGAS00001004363

• tFL with a PMBL GE signature

  We investigated the clinicopathologic features of 5 follicular lymphomas (FL) that transformed to morphologic diffuse large B-cell lymphomas (DLBCL) and had a primary mediastinal large B-cell lymphoma (PMBL)-like gene expression profile. None of the transformed FL (tFL) arose in the mediastinum, all cases tested had a germinal center B-cell immunophenotype, 20% were CD30+, 60% CD23+, 80% MAL+, 20% CD200+, and 0% CD273/PDL2+. Whole exome sequencing detected alterations in genes associated with both FL/DLBCL (CREBBP, KMT2C, KMT2D, ARID1A, HIST1 members, and TNFRSF14) and PMBL (JAK-STAT pathway genes, B2M, and CD58). Copy number (CN) analysis detected gains/amplification of REL in 60% and STAT6 in 40%, gains of chromosome 16, including IL4R, in 40%, and both deletions and gains of 11q in 20%. This cohort, although limited in size, supports a subset of tFL that has blended features between FL/DLBCL and PMBL. Despite having some features that are less common in DLBCL (MAL and CD23 expression, JAK-STAT activation), these cases lack the most characteristic CN alteration seen in PMBL (9p24.1 gain/amplification). This cohort expands the biologic heterogeneity of tFL, illustrating a subset with blended FL/DLBCL and PMBL features, with potential treatment implications that include the use of novel PMBL-targeted therapies.

  Study EGAS00001005870

• The malaria-protective human glycophorin structural variant DUP4 shows somatic mosaicism and association with hemoglobin levels

  Glycophorin A and glycophorin B are red blood cell surface proteins that are both receptors for the parasite Plasmodium falciparum, which is the principal cause of malaria in sub-Saharan Africa. DUP4 is a complex structural genomic variant that carries extra copies of a glycophorin A - glycophorin B fusion gene, and has a dramatic effect on malaria risk by reducing the risk of severe malaria by up to 40% Using fiber-FISH and Illumina sequencing, we validate the structural arrangement of the glycophorin locus in the DUP4 variant, and reveal somatic variation in copy number of the glycophorin A- glycophorin B fusion gene. By developing a simple, specific, PCR-based assay for DUP4 we show the DUP4 variant reaches a frequency of 13% in a malaria-endemic village in south-eastern Tanzania. We genotype a substantial proportion of that village and demonstrate an association of DUP4 genotype with hemoglobin levels, a phenotype related to malaria, using a family-based association test. Taken together, we show that DUP4 is a complex structural variant that may be susceptible to somatic variation, and show that it is associated with a malarial-related phenotype in a non-hospitalizedlongitudinally-followed population.

  Study EGAS00001003239

• CRISPR/Cas9-mediated genome editing of Schistosoma mansoni acetylcholinesterase

  Here we report successful gene knock-in (KI) in the eggs of Schistosoma mansoni by combining CRISPR/Cas9 with single-stranded oligodeoxynucleotides (ssODNs). We targeted the acetylcholinesterase (AChE) gene of S. mansoni using two synthetic guide RNAs (gRNAs), X5 and X7, respectively. Liver eggs of S. mansoni were exposed to CRISPR-vector containing X5 or X7 by electroporation. Simultaneously, eggs were transfected with a ssODN donor encoding a stop codon in all six frames. Next generation sequencing analysis revealed that CRISPR/Cas9-mediated editing in S. mansoni eggs resulted in Homology-Directed Repair (HDR) when template DNA ssODN provided. Furthermore, soluble egg antigen (SEA) from AChE-modified eggs exhibited markedly reduced AChE activity compared with controls, indicative that programmed Cas9 cleavage mutated the AChE gene. Following injection of modified schistosome eggs into the tail veins of mice, a significant decrease in granuloma size in the lungs of these animals. Notably, an enhanced Th2 response induced by eggs in lung, and splenocytes small intestine-draining mesenteric lymph node cells was also generated in mice injected with X5-KI eggs in different methods. These findings further demonstrate the power and utility of CRISPR/Cas9-based genome editing for undertaking functional genomics studies in schistosomes.

  Study EGAS00001004455

• Whole Exome Sequencing PPGL

  Whole exome sequencing (WES) libraries were prepared with the SureSelectXT Human All Exon V6+COSMIC target enrichment system (Agilent, 5190-9307) by following kit manufacturer`s instructions. 200 ng of genomic DNA from FFPE tumoral samples and 1 µg from FF or germline samples were used as starting material. SureSelectXT HS (Agilent, G9704A) and SureSelectXT2 (Agilent, G9621A) reagent kits were alternatively used for FFPE derived and good quality samples, respectively. DNAs were fragmented on a Covaris S2 sharing instrument based on initial sample’s integrity. Fragmented DNA samples were processed through subsequent enzymatic treatments of end-repair, dA-tailing, and ligation to Illumina's adapters. Adapter-ligated libraries were PCR amplified with Illumina PE primers and subsequently hybridized to SureSelect Oligo Capture Library Mix. Biotinylated hybrids were captured and the enriched fraction eluted. Libraries were completed by limited-cycle PCR with Illumina TruSeq primers and KAPA HiFi HotStart DNA pol (Roche KK2501). Sample-specific libraries were combined into equimolarly balanced pools and applied to Illumina flow cells for cluster generation. Sequencing was performed on either the Illumina HiSeq2500 or NovaSeq6000 in a paired-end 100bp reads mode to a median target coverage of ~200x and ~100x for tumors and germline, respectively.

  Study EGAS00001006043

• Multi-omics integration reveals only minor long-term molecular and functional sequelae in immune cells of individuals recovered from COVID-19

  The majority of COVID-19 patients experience mild to moderate disease course and recover within a few weeks. An increasing number of studies characterized the long-term changes in the specific anti-SARS-CoV-2 immune responses, but how COVID-19 shapes the innate and heterologous adaptive immune system after recovery is less well known. To comprehensively investigate the post-SARS-CoV-2 infection sequelae on the immune system, we performed a multi-omics study by integrating single-cell RNA-sequencing, single-cell ATAC-sequencing, genome-wide DNA methylation profiling, and functional validation experiments in convalescent COVID-19 and healthy individuals. We showed that immune responses generally recover without major sequelae after COVID-19. However, subtle differences persist at the transcriptomic level in monocytes, with downregulation of the interferon pathway, while DNA methylation also displays moderate changes in convalescent COVID-19 individuals. However, these differences did not affect the cytokine production capacity of PBMCs upon different bacterial, viral, and fungal stimuli, although baseline release of IL-1RA and IFN?? was higher in convalescent individuals. In conclusion, despite minor differences in epigenetic and transcriptional programs, the immune system of convalescent COVID-19 patients largely recovers to the homeostatic level of healthy individuals.

  Study EGAS00001005529

• Stereotyped B-cell responses are linked to IgG constant region polymorphisms in multiple sclerosis

  Clonally related B cells infiltrate the brain, meninges, and cerebrospinal fluid of multiple sclerosis (MS) patients, but the mechanisms driving the B-cell response and shaping the immunoglobulin repertoires remain unclear. Here, we used single-cell full-length RNA-seq and B-cell receptor reconstruction to simultaneously assess the phenotypes, isotypes, constant region polymorphisms, and the paired heavy- and light-chain repertoires in intrathecal B cells. We detected extensive clonal connections between the memory B cell and antibody-secreting cell compartments and observed clonally related cells of different isotypes, including IgM/IgG1, IgG1/IgA1, IgG1/IgG2, and IgM/IgA1. There was a strong dominance of the G1m1 allotype constant region polymorphisms in antibody-secreting cells, but not in memory B cells. Tightly linked to the G1m1 allotype, we found a preferential pairing of the immunoglobulin heavy-chain variable (IGHV)4 gene family with the κ variable (IGKV)1 gene family. The IGHV4-39 gene was most used and showed the highest frequency of pairing with IGKV1-5 and IGKV1(D)-33. These results link IgG constant region polymorphisms to stereotyped B-cell responses in MS and indicate that the intrathecal B-cell response in these patients could be directed against structurally similar epitopes.

  Study EGAS00001005745

• Regulatory T cell transcriptomic reprogramming characterizes adverse events by checkpoint inhibitors in solid tumors

  Immune checkpoint inhibitors (ICI), which target immune regulatory pathways to unleash antitumor responses, have revolutionized cancer immunotherapy. Despite the remarkable success of ICI immunotherapy, a significant proportion of patients whose tumors respond to these treatments develop immune-related adverse events (irAEs) resembling autoimmune diseases. Although the clinical spectrum of irAEs is well characterized, their successful management remains empiric. This is in part because the pathogenic mechanisms involved in the break-down of peripheral tolerance and induction of irAEs remain elusive. Herein, we focused on regulatory T cells (Tregs) in individuals with irAEs because these cells are vital for maintenance of peripheral tolerance, appear expanded in the peripheral blood of individuals with cancer and abundantly express checkpoint molecules, hence representing direct targets of ICI immunotherapy. Our data demonstrate an intense transcriptomic reprogramming of CD4+CD25+CD127– Tregs in the blood of individuals with advanced metastatic melanoma who develop irAEs following ICI immunotherapy, with a characteristic inflammatory, apoptotic and metabolic signature. This inflammatory signature was shared by Tregs from individuals with different types of cancer developing irAEs and individuals with autoimmune diseases. Our findings propose an inflammatory Treg reprogramming as a feature of immunotherapy-induced irAEs, that may facilitate translational approaches aiming to induce robust antitumor immunity without disturbing peripheral tolerance.

  Study EGAS00001004694

• RNA-seq study of human long-term and short-term hematopoietic stem cells from umblical cord blood with lentiviral overexpression of S1PR3

  To investigate the molecular and biological pathways altered by S1PR3OE in human hematopoietic stem cells (HSC), we performed RNA-sequencing (RNA-seq) of LT- and ST-HSC 3 days after transduction with control or S1PR3 overexpression (OE) lentiviral vectors. LT-HSC and ST-HSC from 3 pool of CB lin- were FACS-purified, cells were prestimulated for 4 hours and transduced with lentiviral vectors. At day 3, 2000-5300 BFP+ cells were FACS-purified for RNA isolation with a PicoPure kit. We were able to isolate only 1600-1800 BFP+ cells from LT-HSC control samples as opposed to 4000-5400 BFP+ cells from S1PR3OE samples. Thus, we pooled all control BFP+ LT-HSC cells into one sample for RNA-seq analysis. BFP- LT-HSC from control vector transduction were purified from CB1 as an additional LT-HSC control. Nextera libraries generated from 10 ng RNA from 5 LT-HSC samples (2 controls, 3 S1PR3OE) and 6 ST-HSC samples (3 controls, 3 S1PR3OE) were subjected to 125 bp, paired-end RNA-sequencing on the Illumina HiSeq 2500 with an average of 50 million reads/sample at the Center for Applied Genomics, Sick Kids Hospital.

  Study EGAS00001004798

• Proteogenomics of chronic lymphocytic leukemia

  Cancer heterogeneity at the proteome level may explain differences in therapy response and prognosis beyond the currently established genomic and transcriptomic based diagnostics. The relevance of proteomics for disease classifications remains to be established in clinically heterogeneous cancer entities such as chronic lymphocytic leukemia (CLL). Here, we characterized the proteome and transcriptome in-depth alongside genetic and ex-vivo drug response profiling in a clinically well annotated CLL discovery cohort (n= 68). Unsupervised clustering of the proteome data revealed six subgroups. Five of these proteomic groups were associated with genetic features, while one group was only detectable at the proteome level. This new group was characterized by accelerated disease progression, high spliceosomal protein abundances associated with aberrant splicing, and low B cell receptor signaling protein abundances (ASB-CLL). We developed classifiers to identify ASB-CLL based on its characteristic proteome or splicing signature in two independent cohorts (n= 165, n= 169) and confirmed that ASB-CLL comprises about 20 % of CLL patients. The inferior overall survival observed in ASB-CLL was independent of both TP53- and IGHV mutation status. Our multi-omics analysis refines the classification of CLL and highlights the potential of proteomics to improve cancer patient stratification beyond genetic and transcriptomic profiling. This is the transcriptomics data of the discovery cohort.

  Study EGAS00001005746

• Clinical genome sequencing uncovers potentially targetable truncations and fusions of MAP3K8 in spitzoid and other melanomas

  Spitzoid melanoma is a specific morphologic variant of melanoma that most commonly affects children and adolescents and ranges on the spectrum of malignancy from low-grade to overtly malignant. These tumors are generally driven by fusions of ALK, RET, NTRK1/3, MET, ROS1 and BRAF. However, in approximately 50% of cases, no genetic driver has been established. Clinical whole genome and transcriptome sequencing (RNA-Seq) of a spitzoid tumor from an adolescent revealed a novel gene fusion of MAP3K8, encoding a serine−threonine kinase that activates MEK. The patient who had exhausted all other therapeutic options was treated with a MEK inhibitor and underwent a transient clinical response. We subsequently analyzed spitzoid tumors from 49 patients by RNA-Seq and found in-frame fusions or C-terminal truncations in 33% of cases. The fusion transcripts and truncated genes all contained MAP3K8 exons 1-8 but lacked the autoinhibitory final exon. Data mining of RNA-Seq from the Cancer Genome Atlas (TCGA) uncovered analogous MAP3K8 rearrangements in 1.5% of adult melanomas. Thus, MAP3K8 rearrangements - uncovered by comprehensive clinical sequencing of a single case - are the most common genetic event in spitzoid melanoma, are present in adult melanomas, and could be amenable to MEK inhibition.

  Study EGAS00001003430

• Androgen receptor blockade promotes response to BRAF/MEK-targeted therapy

  We studied a group of melanoma patients treated with neoadjuvant BRAF/MEK-targeted therapy (NCT02231775, n=51), and observed significantly higher rates of major pathologic response (MPR= <10% viable tumor at resection) and improved recurrence-free survival (RFS) in females versus males (MPR-66% versus 14%, p=0.001; RFS-64% versus 32% at 2 years, p=0.021). Findings were validated in a several additional cohorts 2-4 patients with unresectable metastatic melanoma treated with BRAF and/or MEK-targeted therapy (n=664 patients in total), demonstrating improved progression-free survival (PFS) and overall survival (OS) in females versus males in several of these studies. Studies in pre-clinical models demonstrated significantly impaired anti-tumor activity in male versus female BRAF/MEK-treated mice (p=0.006), with significantly higher expression of androgen receptor (AR) in tumors of male and female BRAF/MEK-treated mice versus control (p=0.0006 and 0.0025). Pharmacologic inhibition of AR signaling improved responses to BRAF/MEK-targeted therapy in male and female mice (p=0.018 and p=0.003), whereas induction of AR signaling (via testosterone administration) was associated with significantly impaired response to BRAF/MEK-targeted therapy in males and females (p=0.021 and p<0.0001).

  Study EGAS00001006196

• WXS and RNA-seq for 22 patients treated with radiation + immunotherapy

  Over 500 active clinical trials are investigating combination radiation (RT) and immune checkpoint blockade (ICB) as a cancer treatment; however, the majority of trials have found no positive interaction. We performed a comprehensive molecular analysis of a randomized phase I clinical trial of patients with non-small cell lung cancer (NSCLC) treated with concurrent or sequential ablative RT and ICB. Concurrent treatment was superior to sequential treatment in augmenting local and distant tumor responses and improving overall survival in a subset of patients with highly aneuploid tumors, but not in those with less aneuploid tumors. Our analysis suggested that concurrent treatment eliminated immunologically cold aneuploid tumors better than sequential treatment and improved patient outcomes. In addition, we report that RT alone decreased intratumoral cytotoxic effector T cell and adaptive immune gene expression signatures, in contrast to upregulation of key immune pathways after concurrent administration of RT and ICB. Our findings distinguish the differential genomic and transcriptomic effects of RT versus RT and ICB and challenge the prevailing paradigm that local ablative RT beneficially stimulates the immune response. We propose the use of tumor aneuploidy as a novel biomarker and therapeutic target in personalizing treatment approaches for patients with NSCLC treated with RT and ICB.

  Study EGAS00001006212

• Sex chromosome aneuploidies give rise to changes in the circular RNA profile

  Purpose The landscape of circular RNAs (circRNAs), an important class of non-coding RNAs that regulate gene expression, has never been described in human disorders of sex chromosome aneuploidies. We profiled circRNAs in Turner syndrome females (45,X;TS) and Klinefelter syndrome males (47,XXY; KS) to investigate how circRNAs respond to a missing or an extra X chromosome. Methods Samples of blood, muscle and fat were collected from individuals with TS (n = 33) and KS (n = 22) and from male (n = 16) and female (n = 44) controls. CircRNAs were identified using a combination of circRNA identification pipelines (CIRI2, CIRCexplorer2 and circRNA_finder). Results Differential expression of circRNAs was observed throughout the genome in TS and KS, in all tissues. The host-genes from which several of these circRNAs were derived, were associated with known phenotypic traits. Furthermore, several differentially expressed circRNAs had the potential to capture micro RNAs that targeted protein-coding genes with altered expression in TS and KS. Conclusion Sex chromosome aneuploidies introduce pervasive changes in the circRNA transcriptome, demonstrating that the genomic changes in these syndromes are more complex than hitherto thought. CircRNAs may help explain some of the genomic and phenotypic traits observed in these syndromes.

  Study EGAS00001006404

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__clear_cell_renal_cell_carcinoma___2

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

  Study EGAS00001003542

• Whole genome sequencing delineates regulatory, structural, and cryptic splice variants in early onset cardiomyopathy

  Cardiomyopathy (CMP) is a heritable disorder. Over 50% cases are gene-elusive on clinical gene panel testing. The contribution of variants in non-coding DNA elements that result in cryptic splicing and regulate gene expression has not been explored. We analyzed whole genome sequencing (WGS) data in a discovery cohort of 209 pediatric CMP patients and 1,953 independent replication genomes and exomes. We searched for protein-coding variants, and non-coding variants predicted to affect the function or expression of genes. Thirty-nine % cases harbored pathogenic coding variants in known CMP genes, and 5% harbored high-risk loss-of-function (LoF) variants in additional candidate CMP genes. Fifteen % harbored high-risk regulatory variants in promoters and enhancers of CMP genes (Odds ratio 2.25, p=6.70×10-7 versus controls). Genes involved in α-dystroglycan glycosylation (FKTN, DTNA) and desmosomal signaling (DSC2, DSG2) were most highly enriched for regulatory variants (Odds ratio 6.7-58.1). Functional effects were confirmed in patient myocardium and reporter assays in human cardiomyocytes, and in zebrafish CRISPR knockouts. We provide strong evidence for the genomic contribution of functionally active variants in new genes and in regulatory elements of known CMP genes to early-onset CMP.

  Study EGAS00001004929

• Single-cell RNA sequencing reveals the mesangial identity and species diversity of glomerular cell transcriptomes

  Molecular characterization of the individual cell types in human kidney as well as model organisms are critical in defining organ function and understanding translational aspects of biomedical research. Previous studies have uncovered gene expression profiles of several kidney glomerular cell types, however, important cells, including mesangial (MCs) and glomerular parietal epithelial cells (PECs), were missing or incompletely described, and a systematic comparison between mouse and human kidney is lacking. To this end, we used Smart-seq2 to profile 4332 individual glomerulus-associated cells isolated from human living donor renal biopsies and mouse kidney. The analysis revealed genetic programs for all four glomerular cell types (podocytes, glomerular endothelial cells, MCs and PECs) as well as rare glomerulus-associated macula densa cells (MDCs). Importantly, we detected heterogeneity in glomerulus-associated Pdgfrb-expressing cells, including bona fide intraglomerular MCs with the functionally active phagocytic molecular machinery, as well as a unique mural cell type located in the central stalk region of the glomerulus tuft. Furthermore, we observed remarkable species differences in the individual gene expression profiles of defined glomerular cell types that highlight translational challenges in the field and provide a guide to design translational studies.

  Study EGAS00001004943

• Significant and pervasive effects of RNA degradation on Nanopore direct RNA sequencing

  Abstract Oxford Nanopore direct RNA sequencing (DRS) is capable of sequencing complete RNA molecules and accurately measuring gene and isoform expression. However, as DRS is designed to profile intact RNA, expression quantification may be more heavily dependent upon RNA integrity than alternate RNA-seq methodologies. It is currently unclear how RNA degradation impacts DRS or if it can be corrected for. To assess the impact of RNA integrity on DRS we performed a degradation time-series using SH-SY5Y neuroblastoma cells. Our results demonstrate that degradation is a significant and pervasive factor that can bias DRS measurements, including a reduction in library complexity resulting in an overrepresentation of short genes and isoforms. Degradation also biases differential expression analyses; however, we find explicit correction can almost fully recover meaningful biological signal. In addition, DRS provided less biased profiling of partially degraded samples than nanopore cDNA-PCR sequencing. Overall, we find samples with RIN > 9.5 can be treated as undegraded and samples with RIN > 7 can be utilised for DRS with appropriate correction. These results establish the suitability of DRS for a wide range of samples, including partially degraded in-vivo clinical and post-mortem samples, whilst limiting the confounding effect of degradation on expression quantification.

  Study EGAS00001006542

• Host pathogen interaction long read transcriptome

  Our incomplete knowledge of the human transcriptome impairs the detection of disease-causing variants, in particular in transcripts only expressed under certain conditions. These transcripts are often lacking from reference transcript sets, such as Ensembl/GENCODE and RefSeq, and could be relevant for establishing genetic diagnoses. We present SUsPECT (Solving Unsolved Patient Exomes/gEnomes using Custom Transcriptomes), a pipeline based on the Ensembl Variant Effect Predictor (VEP) to predict variant impact on custom transcript sets, such as those generated by long-read RNA-sequencing, for downstream prioritization. Our pipeline predicts the functional consequence and likely deleteriousness scores for missense variants in the context of novel open reading frames predicted from any transcriptome. We demonstrate the utility of SUsPECT by uncovering potential mutational mechanisms of pathogenic variants in ClinVar that are predicted to be benign using the reference transcript annotation. In further support of SUsPECT’s utility, we identified an enrichment of immune-related variants predicted to have a more severe molecular consequence when annotating with a newly generated transcriptome from stimulated immune cells instead of the reference transcriptome. Our pipeline outputs crucial information for further prioritization of potentially disease-causing variants for any disease and will become increasingly useful as more long-read RNA sequencing datasets become available.

  Study EGAS00001006779

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

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

  Study EGAS00001004292

• T cell responses of ALS patients

  In this project, we aimed to understand how T cell responses contribute to the disease progression of amyotrophic lateral sclerosis (ALS). The present data is on single-cell sequencing isolated from human cerebrospinal fluid (CSF) cells from both ALS patients (n=5) and controls (n=4). This analysis was conducted as part of a bigger project which is summarized in the section below. We used flow cytometry to define T cell subsets and phenotypes in blood and CSF samples collected at the time of diagnosis on a cohort of 89 newly diagnosed ALS patients in Stockholm, Sweden. High frequency of CD4+FOXP3- effector T cells in blood and CSF was associated with poor survival whereas high frequency of activated regulatory T (Treg) cells and high ratio between activated and resting Treg cells in blood was associated with better survival. T cell profiles also predicted disease progression rate. On an independent cohort of cases and controls, we used single cell transcriptomics data to demonstrate that ALS patients had altered T cell gene expression patterns and clonally expanded CD4+ and CD8+ T cells in CSF. In summary, T cell responses contribute to disease progression of ALS, supporting modulation of adaptive immunity as a viable therapeutic option.

  Study EGAS00001006675

• Longitudinal evaluation of serum microRNAs as biomarkers for neuroblastoma burden and therapeutic p53 reactivation

  Accurate assessment of treatment response and residual disease is indispensable for the evaluation of cancer treatment efficacy. However, performing tissue biopsies for longitudinal follow-up poses a major challenge in the management of solid tumours like neuroblastoma. In the present study, we evaluated whether circulating miRNAs are suitable to monitor neuroblastoma tumour burden and whether treatment-induced changes of miRNA abundance in the tumour are detectable in serum. We performed small RNA sequencing on longitudinally collected serum samples from mice carrying orthotopic neuroblastoma xenografts that were exposed to treatment with idasanutlin or temsirolimus. We identified 57 serum miRNAs to be differentially expressed upon xenograft tumour manifestation, out of which 21 were also found specifically expressed in the serum of human high-risk neuroblastoma patients. The murine serum levels of these 57 miRNAs correlated with tumour tissue expression and tumour volume, suggesting potential utility for monitoring tumour burden. In addition, we describe serum miRNAs that dynamically respond to p53 activation following treatment of engrafted mice with idasanutlin. We identified idasanutlin-induced serum miRNA expression changes upon one day and 11 days of treatment. By limiting to miRNAs with a tumour-related induction, we put forward hsa-miR-34a-5p as a potential pharmacodynamic biomarker of p53 activation in serum.

  Study EGAS00001006678

• Hexanucleotide repeat expansions in C9orf72 alter microglial responses and prevent a coordinated glial reaction in ALS

  Neuroinflammation is an important hallmark in amyotrophic lateral sclerosis (ALS). Experimental evidence has highlighted a role of microglia in the modulation of motor neuron degeneration. However, the exact contribution of microglia to both sporadic and genetic forms of ALS is still unclear. We generated single nuclei profiles of spinal cord and motor cortex from sporadic and C9orf72 ALS patients, as well as controls. We particularly focused on the transcriptomic responses of both microglia and astrocytes. We confirmed that C9orf72 is highly expressed in microglia and shows a diminished expression in carriers of the hexanucleotide repeat expansion (HRE). This resulted in an impaired response to disease, with specific deficits in phagocytic and lysosomal transcriptional pathways. Astrocytes also displayed a dysregulated response in C9orf72 ALS patients, remaining in a homeostatic state. This suggests that C9orf72 HRE alters a coordinated glial response, which ultimately would increase the risk for developing ALS. Our results indicate that C9orf72 HRE results in a selective microglial loss-of-function, likely impairing microglial-astrocyte communication and preventing a global glial response. This is relevant as it indicates that sporadic and familial forms of ALS may present a different cellular substrate, which is of great importance for patient stratification and treatment.

  Study EGAS00001006711

• T cell landscape definition by multi-omics identifies galectin-9 as novel immunotherapy target in chronic lymphocytic leukemia (CLL)

  Failure of immunotherapy after applying checkpoint inhibitors or CAR-T cells is linked to T cell exhaustion. Here, we explored the T cell landscape in chronic lymphocytic leukemia (CLL) by single-cell omics analyses of blood, bone marrow and lymph node samples of patients and spleen samples of a CLL mouse model. By single-cell RNA- sequencing, mass cytometry (CyTOF), and multiplex image analysis of tissue microarrays, we defined the spectrum of phenotypes and transcriptional programs of T cells and and their differentiation state trajectories. We identified disease-specific accumulation of distinct regulatory T cell subsets and T cells harboring an exhausted phenotype exclusively in the CLL lymph node tissue. Integration of TCR data revealed a clonal expansion of CD8+ precursor exhausted T cells, suggesting their reactivity for CLL cells. Interactome analyses identified the TIM3 ligand Galectin-9 as novel immunoregulatory molecule in CLL. Blocking of Galectin-9 in CLL-bearing mice slowed down disease development and reduced the number of TIM3 expressing T cells. Galectin-9 expression correlated with shorter survival of CLL patients. Thus, Galectin-9 contributes to immune escape in CLL and represents a novel target for immunotherapy.

  Study EGAS00001006864

• Atypical teratoid/rhabdoid tumoroids reveal subgroup-specific drug vulnerabilities WGS

  Atypical teratoid/rhabdoid tumors (ATRTs) represent a rare, but aggressive pediatric brain tumor entity. They are genetically uniquely defined by alterations in the SWI/SNF chromatin remodeling complex members SMARCB1 or SMARCA4. ATRTs can be further classified in different molecular subgroups based on their epigenetic profiles. Although recent studies suggest that the different subgroups have distinct clinical features, subgroup-specific treatment regimens have not been developed thus far. This is hampered by the lack of pre-clinical in vitro models representative of the different molecular subgroups. Here, we describe the establishment of ATRT tumoroid models from the ATRT-MYC and ATRT-SHH subgroups. We demonstrate that ATRT tumoroids retain subgroup-specific epigenetic and gene expression profiles. High throughput drug screens on our ATRT tumoroids revealed distinct drug sensitivities between and within ATRT-MYC and ATRT-SHH subgroups. Whereas ATRT-MYC universally displayed high sensitivity to multi-targeted tyrosine kinase inhibitors, ATRT-SHH showed a more heterogeneous response with a subset showing high sensitivity to Notch inhibitors, which corresponded to high expression of Notch receptors. Our ATRT tumoroids represent the first pediatric brain tumor organoid model, providing a representative pre-clinical model which enables the development of subgroup-specific therapies.

  Study EGAS00001006865

• Atypical teratoid/rhabdoid tumoroids reveal subgroup-specific drug vulnerabilities RNA-Seq

  Atypical teratoid/rhabdoid tumors (ATRTs) represent a rare, but aggressive pediatric brain tumor entity. They are genetically uniquely defined by alterations in the SWI/SNF chromatin remodeling complex members SMARCB1 or SMARCA4. ATRTs can be further classified in different molecular subgroups based on their epigenetic profiles. Although recent studies suggest that the different subgroups have distinct clinical features, subgroup-specific treatment regimens have not been developed thus far. This is hampered by the lack of pre-clinical in vitro models representative of the different molecular subgroups. Here, we describe the establishment of ATRT tumoroid models from the ATRT-MYC and ATRT-SHH subgroups. We demonstrate that ATRT tumoroids retain subgroup-specific epigenetic and gene expression profiles. High throughput drug screens on our ATRT tumoroids revealed distinct drug sensitivities between and within ATRT-MYC and ATRT-SHH subgroups. Whereas ATRT-MYC universally displayed high sensitivity to multi-targeted tyrosine kinase inhibitors, ATRT-SHH showed a more heterogeneous response with a subset showing high sensitivity to Notch inhibitors, which corresponded to high expression of Notch receptors. Our ATRT tumoroids represent the first pediatric brain tumor organoid model, providing a representative pre-clinical model which enables the development of subgroup-specific therapies.

  Study EGAS00001006866

• Atypical teratoid/rhabdoid tumoroids reveal subgroup-specific drug vulnerabilities DNA-Methylation

  Atypical teratoid/rhabdoid tumors (ATRTs) represent a rare, but aggressive pediatric brain tumor entity. They are genetically uniquely defined by alterations in the SWI/SNF chromatin remodeling complex members SMARCB1 or SMARCA4. ATRTs can be further classified in different molecular subgroups based on their epigenetic profiles. Although recent studies suggest that the different subgroups have distinct clinical features, subgroup-specific treatment regimens have not been developed thus far. This is hampered by the lack of pre-clinical in vitro models representative of the different molecular subgroups. Here, we describe the establishment of ATRT tumoroid models from the ATRT-MYC and ATRT-SHH subgroups. We demonstrate that ATRT tumoroids retain subgroup-specific epigenetic and gene expression profiles. High throughput drug screens on our ATRT tumoroids revealed distinct drug sensitivities between and within ATRT-MYC and ATRT-SHH subgroups. Whereas ATRT-MYC universally displayed high sensitivity to multi-targeted tyrosine kinase inhibitors, ATRT-SHH showed a more heterogeneous response with a subset showing high sensitivity to Notch inhibitors, which corresponded to high expression of Notch receptors. Our ATRT tumoroids represent the first pediatric brain tumor organoid model, providing a representative pre-clinical model which enables the development of subgroup-specific therapies.

  Study EGAS00001006881

• Human leukocyte antigen alleles associate with COVID-19 vaccine immunogenicity and risk of breakthrough infection

  SARS-CoV-2 vaccine immunogenicity varies between individuals, and immune responses correlate with vaccine efficacy. Using data from 1,076 participants enrolled in ChAdOx1 nCov-19 vaccine efficacy trials in the United Kingdom, we find that inter-individual variation in normalised antibody responses against SARS-CoV-2 spike (S) and its receptor binding domain (RBD) at 28 days following first vaccination shows genome-wide significant association with major histocompatibility complex (MHC) class II alleles. The most statistically significant association with higher levels of anti-RBD antibody was HLA-DQB1*06 (P=3.2 x 10-9), which we replicate in 1,677 additional vaccinees. Individuals carrying HLA-DQB1*06 alleles were less likely to experience PCR-confirmed breakthrough infection during the Wuhan and Alpha-variant waves compared to non-carriers (HR 0.63, 0.42-0.93, P=0.02). We identify a distinct S-derived peptide that is predicted to bind effectively to HLA-DQB1*06 compared to other equivalent alleles, and find evidence of increased spike-specific memory B-cell responses in HLA-DQB1*06 carriers at 84 days post first vaccination. Our results demonstrate association of HLA type with COVID-19 vaccine antibody response and risk of breakthrough infection, with implications for future vaccine design and implementation.

  Study EGAS00001006909

• Quantitative analysis of a novel DNA hypermethylation panel for lung cancer diagnosis using bronchial specimens

  Non-diagnostic findings in Transbronchial Lung Biopsy (TBLB) and Endobronchial Ultrasound-guided Transbronchial Lung Biopsy (EBUS-TBLB) are not uncommon. A challenge is to improve lung cancer detection by TBLB/EBUS-TBLB. In order to improve the diagnostic yield of bronchoscopy, we used 850K methylation chip to identify methylation sites that distinguish malignant from benign lung nodules.We found that the combination of HOXA7, SHOX2,and SCT methylation analysis has the best diagnostic yield in 123 bronchial washing (sensitivity: 74·1%; AUC: 0·851) and 172 brushing samples (sensitivity: 86·1%; AUC: 0·915). We developed a kit comprising these three genes and validated the kit in 329 unique bronchial washing samples, 397 unique brushing samples, and 179 unique patients with both washing and brushing samples. The diagnostic accuracy of the panel alone in lung cancer diagnosis was 86·9%, 91·2% and 95% in bronchial washing, brushing, and washing+brushing samples, respectively. The panel's sensitivity in combination with cytology, rapid on-site evaluation (ROSE), and histology in lung cancer diagnosis was 90·8% and 95.8% in bronchial washing and brushing samples, respectively, and 100% in washing+brushing samples. Quantitative analysis of the three-gene panel improves lung cancer diagnosis by bronchoscopy.

  Study EGAS00001007089

• Long-term organoid culture of a small intestinal neuroendocrine tumor

  Gastroenteropancreatic neuroendocrine tumors (GEP-NETs) are rare and highly heterogeneous neoplasms whose incidence has markedly increased over the last decades. A grading system based on the tumor cells' proliferation index predicts high-risk for G3 NETs. However, low-to-intermediate grade (G1/G2) NETs have an unpredictable clinical course that varies from indolent to highly malignant. Cultures of human cancer cells enable to perform functional perturbation analyses that are instrumental to enhance our understanding of cancer biology. To date, no tractable and reliable long-term culture of G1/G2 NET has been reported to permit disease modeling and pharmacological screens. Here, we report of the first long-term culture of a G2 metastatic small intestinal NET that preserves the main genetic drivers of the tumor and retains expression patterns of the endocrine cell lineage. Replicating the tissue, this long-term culture showed a low proliferation index, and yet it could be propagated continuously without dramatic changes in the karyotype. The model was readily available for pharmacological screens using targeted agents and as expected, showed low tumorigenic capacity in vivo. Overall, this is the first long-term culture of NETs to faithfully recapitulate many aspects of the original neuroendocrine tumor.

  Study EGAS00001007093

• Human liver NPCs single cell project

  Independent of their inflammatory phenotype, macrophages are key orchestrators of hepatic metabolism. Non-alcoholic fatty liver disease (NAFLD) often occurs in obese individuals and is among the most common causes of cirrhosis, the terminal chronic liver disease that may necessitate liver transplantation. While multiple populations of macrophages have been described in the human liver, their function and turnover in obese patients at high risk of developing NAFLD and cirrhosis is currently unknown. Herein we identified a specific human population of resident liver myeloid cells that protects against the metabolic impairment associated with obesity. By studying the turnover of liver myeloid cells in individuals undergoing liver transplantation using markers of donor-recipient mismatch, we made the novel discovery that liver myeloid cell turnover differs between humans and mice. Using single cell techniques and flow cytometry we determined that the proportion of the protective resident liver myeloid cells, denoted liver myeloid cells 2 (LM2), decreases during obesity. Functional validation approaches using human 2D and 3D cultures revealed that the presence of LM2 ameliorates the oxidative stress associated with obese conditions. Our study indicates that resident myeloid cells could be a therapeutic target to decrease the oxidative stress associated with NAFLD.

  Study EGAS00001007194

• Transcriptomics identifies blunted immunomodulatory effects of vitamin D in people with multiple sclerosis

  Vitamin D deficiency is a risk factor for developing multiple sclerosis (MS). However, the immune effects of vitamin D in people with MS are not well understood. We analyzed transcriptomic datasets generated by RNA sequencing of immune cell subsets (CD4+, CD8+ T cells, B cells, monocytes) from 33 healthy controls and 33 untreated MS cases. We utilized a traditional bioinformatic pipeline and weighted gene co-expression network analysis (WGCNA) to determine genes and pathways correlated with endogenous vitamin D. In controls, CD4+ and CD8+ T cells had 1079 and 1188 genes, respectively, whose expressions were correlated with plasma 25-hydroxyvitamin D level (P<0.05). Functional enrichment analysis identified association with TNF-alpha and MAPK signaling. In CD4+ T cells of controls, vitamin D level was associated with expression levels of several genes proximal to multiple sclerosis risk loci (P=0.01). Genes differentially associated with endogenous vitamin D by case-control status were enriched in TNF-alpha signaling via NF-κB. WGCNA suggested a blunted response to vitamin D in cases relative to controls. Collectively, our findings provide further evidence for the immune effects of vitamin D, and demonstrate a differential immune response to vitamin D in cases relative to controls, highlighting a possible mechanism contributing to MS pathophysiology.

  Study EGAS00001007254

• Targeting TRIP13 in Wilms Tumor with Nuclear Export Inhibitors

  Wilms tumor (WT) is the most common renal malignancy of childhood. Despite improvements in the overall survival, relapse occurs in ~15% of patients with favorable histology WT (FHWT). Half of these patients will succumb to their disease. Identifying novel targeted therapies in a systematic manner remains challenging in part due to the lack of faithful preclinical in vitro models. We established ten short-term patient-derived WT cell lines and characterized these models using low-coverage whole genome sequencing, whole exome sequencing and RNA-sequencing, which demonstrated that these ex-vivo models faithfully recapitulate WT biology. We then performed targeted RNAi and CRISPR-Cas9 loss-of-function screens and identified the nuclear export genes (XPO1 and KPNB1) as strong vulnerabilities. We observed that these models are sensitive to nuclear export inhibition using the FDA approved therapeutic agent, selinexor (KPT-330). Selinexor treatment of FHWT suppressed TRIP13 expression, which was required for survival. We further identified in vitro and in vivo synergy between selinexor and doxorubicin, a chemotherapy used in high risk FHWT. Taken together, we identified XPO1 inhibition with selinexor as a potential therapeutic option to treat FHWTs and in combination with doxorubicin, leads to durable remissions in vivo.

  Study EGAS00001007389

• Genomic evolution and transcriptional changes in the evolution of prostate cancer into neuroendocrine and ductal carcinoma types

  Prostate cancer is typically of acinar adenocarcinoma type but can occasionally present as neuro-endocrine and/or ductal type carcinoma. These are associated with clinically aggressive disease and the former often arises on a background of androgen deprivation therapy, although it can al-so arise de novo. Two prostate cancer cases were sequenced by exome capture from archival tis-sue. Case 1 was de novo small cell neuroendocrine carcinoma and ductal adenocarcinoma with 3 longitudinal samples over 5 years. Case 2 was a single time point after development of treatment related neuroendocrine prostate carcinoma. Case 1 showed whole genome doubling in all sam-ples and focal amplification of AR in all samples except the first time point. Phylogenetic analysis revealed a common ancestry for the ductal and small cell carcinoma. Case 2 showed 13q loss (involving RB1) in both adenocarcinoma and small cell carcinoma regions, and 3p gain, 4p loss, 17p loss (involving TP53) in the latter. By using highly curated samples, we demonstrate for the first time that small cell neuroendocrine and ductal prostatic carcinoma can have a common an-cestry and the process of evolution over time. We highlight whole genome doubling in a patient with prostate cancer relapse, reinforcing its poor prognostic nature.

  Study EGAS00001007412

• Genomic evolution and transcriptional changes in the evolution of prostate cancer into neuroendocrine and ductal carcinoma types (RNAseq)

  Prostate cancer is typically of acinar adenocarcinoma type but can occasionally present as neuro-endocrine and/or ductal type carcinoma. These are associated with clinically aggressive disease and the former often arises on a background of androgen deprivation therapy, although it can al-so arise de novo. Two prostate cancer cases were sequenced by exome capture from archival tis-sue. Case 1 was de novo small cell neuroendocrine carcinoma and ductal adenocarcinoma with 3 longitudinal samples over 5 years. Case 2 was a single time point after development of treatment related neuroendocrine prostate carcinoma. Case 1 showed whole genome doubling in all sam-ples and focal amplification of AR in all samples except the first time point. Phylogenetic analysis revealed a common ancestry for the ductal and small cell carcinoma. Case 2 showed 13q loss (involving RB1) in both adenocarcinoma and small cell carcinoma regions, and 3p gain, 4p loss, 17p loss (involving TP53) in the latter. By using highly curated samples, we demonstrate for the first time that small cell neuroendocrine and ductal prostatic carcinoma can have a common an-cestry and the process of evolution over time. We highlight whole genome doubling in a patient with prostate cancer relapse, reinforcing its poor prognostic nature.

  Study EGAS00001007428

• Integrated genomic analysis for HCC

  Liver cancer is a common malignancy and is the third cause of cancer-related death in the world and hepatocellular carcinoma (HCC) is the major type of liver cancer, accounting for ~90% of cases. Only 20 ~ 30% of patients can benefit from combined immunotherapy and targeted therapy due to the high inter-tumor heterogeneity of HCC. This work aims to identify homogeneous and robust molecular subtypes in HCC based on a large, homogenous and well-annotated cohort. We performed whole genome/exome sequencing (WGS/WES) and/or RNA sequencing (RNA-seq) on 529 HCC from 461 patients collected mainly in France. Based on the great consistence between genomic and transcriptomic data, we identified 9 robust HCC subtypes mainly based on driver mutations. And we further characterized HCC subtypes using genomic, transcriptomic and clinicopathological features. Among 9 subtypes, 5 subtypes belonged to chromosome instable tumors, while 3 subtypes belonged to chromosome stable tumors. In addition, our subtypes were associated with prognosis and we also found different distributions of subtypes across various features, such as etiology and gender. Our study provided a robust molecular classification based on a large HCC dataset from Western country, it improves our understanding of the mechanisms of carcinogenesis and will facilitate development of genome-based precision medicine in HCC.

  Study EGAS00001007957

• Detection of cancers three years prior to diagnosis using plasma cell-free DNA

  A major goal of early cancer detection is to identify subclinical disease when the tumor burden is low, so that treatments are more effective. But how early can cancers be detected prior to clinical signs or symptoms? This question can be answered only through the evaluation of participants whose clinical course has not been altered by the study itself. We here describe such an evaluation, performed on prospectively collected plasma samples from the Atherosclerosis Risk in Communities (ARIC) study, including 26 participants diagnosed with cancer and 26 matched controls. At the index time point, eight of these 52 participants scored positively with a multicancer early detection (MCED) blood test. All eight of these participants were diagnosed with cancer within 4 months after blood collection. In six of these 8 participants, we were able to assess an earlier plasma sample collected 3.1 to 3.5 years prior to clinical diagnosis. In four of these six participants, the same mutations detected by the MCED test could be identified, but at 8.6 to 79-fold lower levels. These results demonstrate that it is possible to detect circulating tumor DNA (ctDNA) more than three years prior to clinical diagnosis, and provide benchmark sensitivities required for the success of ctDNA-based tests for this purpose.

  Study EGAS00001008068

• Integrative Multi-Omics and Drug Sensitivity Profiling Reveals Potential Predictive Biomarkers in Pediatric Solid Tumors from the INFORM Registry

  Cure rates for childhood malignancies using established therapy protocols have increased to an average of 80% but have reached a plateau. Moreover, survival rates are particularly low for some pediatric tumors (such as high-risk Group3 medulloblastomas, osteosarcomas, Ewing sarcomas, high-risk neuroblastomas and high-grade gliomas) and dismal for patients with relapsed malignancies. A functional drug response profiling platform for pediatric solid and brain tumors has been established within the INFORM program to identify patient-specific vulnerabilities, biomarkers and unravel molecular mechanisms associated with drug response profiles for clinical translation. Here, we perform a multi-omics analysis using drug sensitivity profiles, genomic, and transcriptomic data of 81 pediatric solid tumor samples. The integrative analysis revealed two multi-omics signatures with high drug sensitivity predictivity. One signature highlighted neuroblastoma samples with sensitivity to navitoclax (a BCL-2 family inhibitor). A second signature was significantly specific to a subset of Wilms tumors harboring the SIX1 (Q177R) hotspot mutation, which had high expression of MGAM, PTPN14, STAT4 and KDM2B and high sensitivity to MEK inhibitors. A patient-specific causal interaction network analysis elucidated the potential biological mechanism of action between MEK inhibitors and the SIX1 mutation in n = 6 Wilms tumor samples.

  Study EGAS00001008249

• Spatiotemporal Genomic Profiling of Intestinal Metaplasia Reveals Clonal Dynamics of Gastric Cancer Progression

  Intestinal metaplasia (IM) is a pre-malignant condition of the gastric mucosa associated with increased gastric cancer (GC) risk. We analyzed 1256 gastric samples (1152 IMs) from 692 subjects from a 10-year prospective study. We identified 26 IM driver genes in diverse pathways including chromatin regulation (ARID1A) and intestinal homeostasis (SOX9), largely occurring as subclonal events. Analysis of clonal dynamics between and within subjects, and also longitudinally across time, revealed that IM clones are likely transient but increase in size upon progression to dysplasia, with eventual transmission of genetic events to paired GCs. Single-cell and spatial profiling highlighted changes in tissue ecology and lineage heterogeneity in IM, including an intestinal stem-cell dominant cellular compartment linked to early malignancy. Expanded transcriptome profiling revealed expression-based molecular subtypes of IM, including a body-resident “pseudoantralized” subtype associated with incomplete histology, antral/intestinal cell types, ARID1A mutations, inflammation, and microbial communities normally associated with the healthy oral tract. We demonstrate that combined clinical-genomic models outperform clinical-only models in predicting IMs likely to progress. Our results raise opportunities for GC precision prevention and interception by highlighting strategies for accurately identifying IM patients at high GC risk and a role for microbial dysbiosis in IM progression.

  Study EGAS00001007067

• Swarm Learning to identify COVID-19, tuberculosis and leukemia patients based on blood transcriptomes

  Identification of patients with life-threatening diseases including leukemias or infections such as tuberculosis or COVID-19 is an important goal of modern precision medicine. However, there is an increasing divide between what is technically possible and what is allowed because of privacy legislation. We have recently illustrated that classical machine learning can identify leukemia patients based on their blood transcriptomes. To facilitate integration of any omics data from any data owner world-wide without violating privacy laws, we here introduce Swarm Learning (SL), a decentralized machine learning approach uniting edge computing, artificial intelligence (AI), blockchain and privacy protection without the need for a central coordinator thereby going beyond federated learning. To illustrate its feasibility, using more than 12,000 transcriptomes from peripheral blood mononuclear cells and more than 2,000 peripheral blood transcriptomes we demonstrate that SL of omics data distributed across different individual sites leads to disease classifiers that outperform those developed at individual sites. Yet, SL completely protects local privacy regulations by design. We propose this approach to noticeably accelerate the introduction of precision medicine.

  Dataset EGAD00001006231

• NIPT samples for systematic evaluation of NIPT aneuploidy detection software tools

  This dataset contains two sets of samples. The reference sample set consists of a total of 669 samples that had been reported previously to be euploid by the NIPTIFY screening test. The validation sample set is based on a previously published validation study by Zilina et al. (1), consisting of 423 samples, of which 259 were high-risk pregnancies that had undergone diagnostic invasive prenatal analysis (1). All samples were sequenced with Illumina NextSeq 500 platform, producing 85 bp single-end reads with an average per-sample coverage of 0.32× at the University of Tartu, Institute of Genomics Core Facility, according to the manufacturer’s standard protocols, as described previously (1). This study was performed with the approval of the Research Ethics Committee of the University of Tartu (#315/T-13). 1. Zilina O, Rekker K, Kaplinski L, Sauk M, Paluoja P, Teder H, et al. Creating basis for introducing non‐invasive prenatal testing in the Estonian public health setting. Prenat Diagn [Internet]. 2019 Dec 6;39(13):1262-8. Available from: https://onlinelibrary.wiley.com/doi/abs/10.1002/pd.5578

  Dataset EGAD00001007712

• Convergent somatic mutations in effectors of insulin signalling in chronic liver disease

  Chronic liver disease is associated with metabolic dysregulation, liver failure and hepatocellular carcinoma. We analysed somatic mutations from 1202 genomes across 32 liver samples, including normal controls, alcohol-related and non-alcoholic fatty liver disease. Five of 27 patients with liver disease carried hotspot driver mutations in FOXO1, the major transcription factor downstream of insulin signalling. FOXO1 mutations were independently acquired by up to 5 distinct clones within the same patient’s sample, and impaired insulin-mediated nuclear export of FOXO1. GPAM, which produces storage triacylglycerol from dietary calories, also had significant excess of mutations, similarly exhibiting convergent evolution within biopsies. Telomeres were shorter in diseased than normal liver, with attrition more pronounced in larger clones. Multiple independent acquisitions of drivers within one small liver sample imply that such mutations could affect hundreds of grams of tissue across the whole organ, potentially contributing to systemic metabolic dysfunction.

  Dataset EGAD00001006255

• PSCP_bisulphite analysis in hESCs (2018-08-13)

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progeny. 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/ . This dataset contains all the data available for this study on 2018-08-13.

  Dataset EGAD00001004295

• FASTQ files of total RNA-Seq data from the POPS PET (pre-eclamptic) samples

  All the samples were obtained from the Pregnancy Outcome Prediction–a prospective cohort study of nulliparous women attending the Rosie Hospital, Cambridge (UK) for their dating ultrasound scan between January 14, 2008, and July 31, 2012. Ethical approval for the study was given by the Cambridgeshire 2 Research Ethics Committee (reference number 07/H0308/163) and all participants provided written informed consent. Cases of preeclampsia (PET) were defined on the basis of the 2013 ACOG criteria and cases of small for gestational age (SGA)infants were confined to severe SGA, i.e. a customized birth weight <5th percentile. Chorionic villi from the corresponding placentas (free from decidua, visible infarction, calcification, hematoma, or damage) were collected and processed within 30 minutes of separation from the uterus. After repeated washes in chilled phosphate buffered saline, the samples were placed in RNA later (Applied Biosystems) and stored at -80°C. Total placental RNA was extracted using mirVana Isolation Kit (Ambion). For each placenta, approximately 5 mg of tissue were homogenized in the Lysis/Binding solution for 20 sec at 6 m/s using a bead beater (FastPrep24) and Lysing Matrix D Tubes (MP Biomedicals). The samples were then spun at 13,000 rpm for 5 min at 4°C and the supernatants recovered. Afterwards, the manufacturer's instructions were followed. Immediately after the RNA extraction, placental RNA samples were DNase-treated using DNA-free DNA Removal Kit (Ambion), aliquoted, and stored in -80°C. Quantity and quality of the RNA samples were assessed using the Agilent 2100 Bioanalyzer, the Agilent RNA 6000 Nano Kit (Agilent Technologies), and Qubit fluorometer. Libraries were prepared starting with 300-500 ng of good quality total RNA (RIN ≥7.5) using the TruSeq Stranded Total RNA Library Prep Kit with Ribo-Zero Human/Mouse/Rat (Illumina), according to the manufacturer's instructions. The kit contains 96 uniquely indexed adapter combinations in order to allow pooling of multiple samples prior to sequencing. After determining their size (with the Agilent 2100 Bioanalyzer and the Agilent High Sensitivity DNA Kit by Agilent Technologies) and concentration (by qPCR with the KAPA Illumina ABI Prism Library Quantification Kit, Kapa Biosystems), libraries have been pooled and sequenced (single-end, 125 bp) using a Single End V4 Cluster Kit and an Illumina HiSeq2500 or HiSeq4000 instrument.

  Dataset EGAD00001003508

• Genetics causes of male infertility in 185 patient-parent trios from Netherlands and UK

  The dataset is composed by the raw and processed sequencing data generated from 185 Patients affected by azoospermia or severe oligozoospermia recruited from the Netherlands and the UK.

  Dataset EGAD00001007862

• Targeted sequencing of diffuse large B-cell lymphoma tumour samples from the UK's Haematological Malignancy Research Network

  Targeted sequencing was applied to an unselected population-based diffuse large B-cell lymphoma cohort (n=928) diagnosed in the UK's Haematological Malignancy Research Network catchment population of ~4 million (14 centres). DNA extracted from tumour samples was sequenced with a 293-gene panel using the Illumina HiSeq 2500. All data are provided in the CRAM format.

  Dataset EGAD00001008616

• Mutational landscape of normal epithelial cells in Lynch Syndrome patients

  Lynch Syndrome (LS) is an autosomal dominant disease conferring a high risk of colorectal cancer due to germline heterozygous mutations in a DNA mismatch repair (MMR) gene. Although cancers in LS patients show elevated somatic mutation burdens, information on mutation rates in normal tissues and understanding of the trajectory from normal to cancer cell is limited. Here we whole-genome sequenced 152 crypts from normal and neoplastic epithelial tissues from LS patients. In normal tissues the repertoire of mutational processes and mutation rates were similar to those found in wild type individuals. A morphologically normal colonic crypt with an increased mutation burden and mutational signatures consistent with MMR deficiency was identified, which may represent a very early stage of LS pathogenesis. Phylogenetic tress of tumour crypts indicated that the most recent ancestor cell of each tumour was already MMR deficient and had experienced multiple clonal evolution cycles. This study demonstrates the genomic stability of epithelial cells with heterozygous germline MMR gene mutations and highlights important differences in the pathogenesis of LS from other colorectal cancer predisposition syndromes.

  Dataset EGAD00001008092

• RNA-seq analysis of transcriptome variation with human ESC subclones

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progeny 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/ .

  Dataset EGAD00001004266

• Dataset for soft_tissue_tumor-RNA

  Rare cancer sequencing data of 243 runs in tumor/control pairs, which were uploaded to umbrella studies. The sequencing was always paired

  Dataset EGAD00001008859

• Neutrophils and emergency granulopoiesis drive immune suppression and an extreme response endotype during sepsis

  This dataset contains raw sequencing files associated with the paper "Neutrophils and emergency granulopoiesis drive immune suppression and an extreme response endotype during sepsis" (https://doi.org/10.1038/s41590-023-01490-5). It is composed of two experiments, which are as follows: 1. Single-cell profiling of whole blood leukocytes: This experiment comprises sequencing files generated using the BD Rhapsody single-cell multi-omics profiling (RNA + protein) platform. This platform was used to profile the whole blood leukocyte population in a cohort of sepsis patients, cardiac surgery controls and healthy controls. There are 48 samples and four files per sample: two paired-end FASTQ files (R1 and R2) corresponding to the RNA-seq library, and two paired-end FASTQ files (R1 and R2) corresponding to the protein profiling (ADT-based AbSeq) library. 2. Single-cell profiling of circulating HSPCs: This experiment comprises sequencing files generated using the 10X single-cell multi-omics (RNA-seq + ATAC-seq) profiling platform. This platform was used to profile circulating HSPCs in a cohort of sepsis patients and healthy controls. There are 5 samples (or plexes), each of which consists of a pool of individuals for whom cells were multiplexed and sequenced as a single sample. There are four files per sample: two paired-end FASTQ files (R1 and R2) corresponding to the RNA-seq library, and two paired-end FASTQ files (R1 and R2) corresponding to the ATAC-seq library. Because these samples consists of multiplexed pools, index files are also provided. These can be used for sample deconvolution using the Cell Ranger pipeline. NOTE: Due to EGA's constrains in the number of files permitted per sample, index files (I1 and I2) for the HSPC data set are provided as a separate experiment named "FASTQ index files (I1 and I2) for deconvolution of 10X single-cell multiomics libraries".

  Dataset EGAD00001010927

• EPICC: Evolutionary Predictions in Colorectal Cancer

  ATAC-seq profiling bam files from colorectal carcinoma and adenoma.

  Dataset EGAD00001008140

• FASTQ files of total RNA-Seq data from the POPS SGA (Small for Gestational Age) samples

  All the samples were obtained from the Pregnancy Outcome Prediction–a prospective cohort study of nulliparous women attending the Rosie Hospital, Cambridge (UK) for their dating ultrasound scan between January 14, 2008, and July 31, 2012. Ethical approval for the study was given by the Cambridgeshire 2 Research Ethics Committee (reference number 07/H0308/163) and all participants provided written informed consent. Cases of preeclampsia (PET) were defined on the basis of the 2013 ACOG criteria and cases of small for gestational age (SGA)infants were confined to severe SGA, i.e. a customized birth weight <5th percentile. Chorionic villi from the corresponding placentas (free from decidua, visible infarction, calcification, hematoma, or damage) were collected and processed within 30 minutes of separation from the uterus. After repeated washes in chilled phosphate buffered saline, the samples were placed in RNA later (Applied Biosystems) and stored at -80°C. Total placental RNA was extracted using mirVana Isolation Kit (Ambion). For each placenta, approximately 5 mg of tissue were homogenized in the Lysis/Binding solution for 20 sec at 6 m/s using a bead beater (FastPrep24) and Lysing Matrix D Tubes (MP Biomedicals). The samples were then spun at 13,000 rpm for 5 min at 4°C and the supernatants recovered. Afterwards, the manufacturer's instructions were followed. Immediately after the RNA extraction, placental RNA samples were DNase-treated using DNA-free DNA Removal Kit (Ambion), aliquoted, and stored in -80°C. Quantity and quality of the RNA samples were assessed using the Agilent 2100 Bioanalyzer, the Agilent RNA 6000 Nano Kit (Agilent Technologies), and Qubit fluorometer. Libraries were prepared starting with 300-500 ng of good quality total RNA (RIN ≥7.5) using the TruSeq Stranded Total RNA Library Prep Kit with Ribo-Zero Human/Mouse/Rat (Illumina), according to the manufacturer's instructions. The kit contains 96 uniquely indexed adapter combinations in order to allow pooling of multiple samples prior to sequencing. After determining their size (with the Agilent 2100 Bioanalyzer and the Agilent High Sensitivity DNA Kit by Agilent Technologies) and concentration (by qPCR with the KAPA Illumina ABI Prism Library Quantification Kit, Kapa Biosystems), libraries have been pooled and sequenced (single-end, 125 bp) using a Single End V4 Cluster Kit and an Illumina HiSeq2500 or HiSeq4000 instrument.

  Dataset EGAD00001003507

• Genetic drivers of epigenetic and transcriptional variation of human immune responses to infection (RNA-seq, ATAC-seq and ChIPmentation)

  RNA-seq, ATAC-seq and ChIPmentation data from monocyte-derived macrophages that were infected with Influenza A virus strain PR8WT, or a matching non-infected control.

  Dataset EGAD00001008422