23014 results for "ega"

in 92.25 milliseconds.

• Error-corrected flow-based sequencing at whole genome scale and its application to circulating cell-free DNA profiling

  Differentiating sequencing errors from true variants is a central genomics challenge calling for error suppression strategies that balance costs when sparse material or ultra-rare variants require high-accuracy methods. For example, circulating cell-free DNA (ccfDNA) sequencing for cancer monitoring is limited by sparsity of circulating tumor DNA (ctDNA), abundance of genomic material in samples, library preparation error rates and sequencing errors. Whole-genome sequencing (WGS) can overcome the low abundance of ccfDNA by integrating signal across the mutation landscape, but higher cost limits adoption compared with targeted panels, where ccfDNA abundance limits maximal coverage depths. Here we applied deep (~120x) lower-cost WGS (Ultima Genomics) with analytic error correction for tumor-informed ctDNA detection within the part-per-million range. We further developed high-coverage duplex error-corrected WGS of ccfDNA, achieving 7.7x10-8 error rates, to assess disease burden in melanoma and urothelial cancer patients without matched tumor sequencing. Together, we show that deeply sequenced error-corrected WGS accurately calls somatic variants, demonstrating the feasibility of WGS for tumor-informed and tumor-agnostic ctDNA detection.

  Study EGAS50000000844

• Whole exome sequencing of bladder tumors

  Bladder carcinomas are immunogenic, and patients with bladder cancer benefit from immune checkpoint therapy. This is correlated to a high tumor mutation burden, which provides a higher number of neoepitopes that can be recognized by tumor-specific CD8+ T cells. Intravesical Bacillus Calmette-Guérin (BCG) is used to treat non-muscle invasive bladder cancer (NMIBC), but its mechanism of action remains elusive. Most lymphocytes appearing in the urine of BCG-treated patients are CD4+ T cells though preclinical studies showed that CD8+ T cells are also necessary for BCG treatment efficacy. It is currently unknown which proportion of patients with non-metastatic bladder cancer develop a spontaneous antitumor CD8+ response, and if BCG treatment influences this response. were prepared using the xGen Dual Index UMI Adapters (IDT), and genomic DNA (gDNA) extracted (DNeasy kit, Qiagen) from formalin-fixed paraffin- embedded (FFPE) or snap-frozen tumor tissue, as well as autologous lymphocytes. Libraries were subjected to paired-end sequencing using Illumina HiSeq4000 with a coverage depth of 250X. Library preparation was carried out with the AllType FASTplex NGS 11 Loci Flex Kit (One Lambda), and sequencing was conducted on the Illumina MiSeq platform.

  Study EGAS50000001248

• MATISSE WES and bulk RNA-sequencing data

  High-risk cutaneous squamous cell carcinoma is frequently seen in frail, elderly patients and often requires mutilating surgical treatment and adjuvant radiotherapy (RT). The MATISSE trial (NCT04620200) was conducted to improve clinical prospects in CSCC patients eligible for extensive curative surgery, and evaluated early treatment response and immune dynamics following ultra-short neoadjuvant immunotherapy prior to standard of care surgery with or without radiotherapy. Fifty patients were randomized to receive two courses of nivolumab 3 mg/kg (NIVO, weeks 0&2) or nivolumab 3 mg/kg (weeks 0&2) and one course of ipilimumab 1 mg/kg (week 0 only, NIVO/IPI before surgery at week 4, with the primary objective being histopathological response based on residual viable tumor cells in the surgical specimen. Ten patients choose not to undergo surgery, of whom nine reached a clinical complete response during follow up, and were disease-free at median 34 months follow up. Serial tumor and blood samples at baseline, week 1, week 2 and week 4, enabled to investigate dynamical changes through bulk WES and RNA sequencing to potentially aid future de-escalation trials.

  Study EGAS50000001003

• Single-cell multidimensional profiling of tumor cell heterogeneity in supratentorial ependymomas

  Supratentorial ependymomas (ST-EPNs) are aggressive and treatment-resistant pediatric brain tumors that form along the human cortex, with suspected origin from radial glia cells. Although histologically similar, ST-EPNs are molecularly classified into six distinct subgroups, many of them remaining understudied. Here, we integrated molecular cell states and spatial architecture to resolve the cellular origins and aberrant differentiation trajectories of ST-EPNs at single-cell resolution. We show that, while all ST-EPNs display a unified origin from neuroepithelial cells, different subgroups harbor highly diverse cell states and differentiation potential associated with neural or ependymal differentiation. Spatial mapping uncovered a higher-order structural tumor organization driven by presence of the mesenchymal state, and refines classical histological features based on detailed molecular resolution. By benchmarking a morpho-molecular cell-state classification system in different disease models and human tumors, we identified four cellular subtypes defined by the architecture of neurite-like extensions that mediate either self-renewal or fast, saltatory invasive phenotypes reminiscent of neurons during neurodevelopment. Collectively, our study provides a comprehensive framework for the study of ST-EPN cellular states, and reveals biologically relevant tumor compartmentalization with potential clinical implication.

  Study EGAS50000001513

• Non-small cell lung cancer molecular subtypes and vulnerability to immunotherapy treatment combinations

  The phase 3 IMpower150 trial in treatment-naïve patients with metastatic non-small-cell lung cancer (NSCLC) demonstrated significantly longer progression-free (PFS) and overall survival (OS) with first-line atezolizumab (anti-PD-L1)-bevacizumab (anti-VEGF)-carboplatin-paclitaxel (ABCP) than with bevacizumab-carboplatin-paclitaxel (BCP). We characterised four molecular NSCLC subtypes identified by unsupervised clustering of transcriptomes of 564 pre-treatment primary tumour samples from IMpower150 using non-negative matrix factorization (NMF1-4). Each subtype had distinct tumour PD-L1 expression levels, epithelial characteristics, immune composition, and treatment outcomes. Both NMF2 (enriched in tumour proliferation signal, macrophages, and monocytes) and NMF4 (enriched in B-cells and T cells) had elevated tumour PD-L1 expression. Of these two, only NMF4 demonstrated PFS and OS benefits with ABCP versus either BCP or atezolizumab-carboplatin-paclitaxel (ACP). Patients with NMF1 (enriched in basal and squamous-like cells) had improved outcomes on ABCP compared with ACP or BCP; those with NMF3 (enriched in adenocarcinoma signatures) showed similar outcomes among treatments. These insights could help inform individualised first-line treatment for metastatic NSCLC.

  Study EGAS50000001272

• Dynamics of checkpoint receptors in γδ T cell subsets are associated with clinical responses during anti-PD-1 immunotherapies

  Gamma delta (γδ) T cells are innate-like lymphocytes with potent anti-tumor properties. Herein, we show that immune checkpoint receptors (ICRs) display differential expression and regulation by the JAK-STAT pathway in Vδ1 and Vδ2 cells and identify constitutive (e.g. TIGIT, PD-1) and inducible (e.g. TIM-3, LAG-3, CTLA-4) ICRs. In melanoma, all γδ T cell subsets downregulated AP-1 transcription factors, but Vδ1 cells specifically expressed high levels of ICR, TOX and inhibitory killer Ig-like receptor (KIR) transcripts, reminiscent of exhaustion. However, patient-derived cells were functionally competent, although induction of LAG-3 and CTLA-4 was impaired. During anti-PD-1 monotherapy, Vδ1 cells specifically bound high levels of therapeutic antibody but only in patients who responded to treatment, revealing a potential new prognostic marker for evaluating the efficacy of IC blockade (ICB) therapy. Finally, expression of KIR genes in Vδ1 cells was downregulated in response to successful ICB therapy. Collectively, our data indicate an intricate relationship between ICRs and γδ T cells and reveal novel approaches by which these cells can be harnessed in order to discern or improve cancer immunotherapy.

  Study EGAS50000001270

• FLT3-targeting and places FLT3 as an ideal therapeutic target to selectively eradicate LSCs, while sparing HSC.

  CRISPR/Cas9 FLT3 knock-out (FLT3-KO) in human HSCs and LSCs followed by functional xenograft assays to test their ability to regenerate human hematopoiesis and leukemia, respectively. FLT3-KO in HSCs from human fetal liver, cord blood and adult bone marrow showed no impairment in multilineage hematopoiesis in primary and secondary xenografts. By contrast, FLT3-KO LSCs from 6 of 7 FLT3-ITD mutated AMLs were able to generate short-term engraftment but were completely exhausted by 12 weeks. Thus, FLT3 is essential for LSC long-term propagation. This dependency was unique to FLT3-ITD AML samples as non-FLT3-ITD AML samples generated leukemic grafts upon FLT3-KO. Transcriptomic analysis revealed that FLT3-KO induced downregulation of DNA repair and cell cycle checkpoints, uniquely in FLT3-ITD AML, but not in healthy HSCs or other AMLs. Our research highlights a critical distinction between healthy HSCs and LSCs: whereas healthy hematopoiesis proceeds unperturbed upon FLT3-KO, FLT3-ITD leukemogenesis is impaired through elimination of LSCs. This evidence underscores the necessity for more potent FLT3-targeting and places FLT3 as an ideal therapeutic target to selectively eradicate LSCs, while sparing HSC.

  Study EGAS50000000437

• Proteomic predictors of individualized nutrient-specific insulin secretion in health and disease

  Population level variation and molecular mechanisms behind insulin secretion in response to carbohydrate, protein, and fat remain uncharacterized despite ramifications for personalized nutrition. We now define prototypical insulin secretion dynamics in response to the three macronutrients in islets from 140 cadaveric donors, including those diagnosed with type 2 diabetes. We leverage the insulin response heterogeneity and use transcriptomics and proteomics to identify molecular pathways of specific nutrient responsiveness. Surprisingly, we find robust insulin secretion to fatty acid stimulus in ~8% of donors, challenging the idea that fat has negligible effects on insulin release. Distinct islet proteomes with differences in metabolic signalling networks convey this hyper-responsiveness to fat relative to carbohydrate. By comparing human islets to human embryonic stem cell-derived islet clusters, we show that, unlike glucose-responsiveness, fat hyper-responsiveness is equivalent and may be a hallmark of functionally immature cells. Our study represents the first comparison of dynamic responses to nutrients and multi-omics analysis in human insulin secreting cells. Responses of different people’s islets to carbohydrate, protein, and fat lay the groundwork for personalized nutrition.

  Study EGAS00001007241

• Context-specific regulatory genetic variation in MTOR dampens neutrophil-T cell crosstalk in pneumonia-associated sepsis, modulating disease

  Sepsis is a heterogeneous clinical syndrome with a high mortality, requiring personalised stratification strategies. Here, we characterise genetic variation that modulates MTOR, a critical regulator of metabolism and immune responses in sepsis. The effects are highly context specific, involving a regulatory element that affects MTOR expression in activated T cells with opposite direction of effect in neutrophils. We show that the G-allele of the lead variant, rs4845987, associated with decreased risk of T2D, reduces MTOR expression in T cells and improves survival in sepsis due to pneumonia, with effects specific to sepsis endotype. Using ex vivo models, we demonstrate that activated T cells promote immunosuppressive sepsis neutrophils through released cytokines, a process dampened by hypoxia and the mTOR inhibitor rapamycin. Our work demonstrates a novel epigenetic mechanism that fine-tunes MTOR transcription and T cell activity via the variant-containing regulatory element, which further exhibits an allelic effect upon vitamin C treatment. These findings reveal how genetic variation interacts with disease state to modulate immune cell-cell communication, providing a patient stratification strategy to inform more effective sepsis treatment.

  Study EGAS50000000894

• A Single-Cell and Spatial Atlas of Early Human Olfactory Development

  The human nasal region is a complex anatomical structure comprising neural crest- and placode-derived lineages, yet its developmental dynamics in humans remain poorly understood due to limited access to fetal tissue and challenges in resolving its spatial and cellular complexity. We built a single-cell and spatial transcriptomic atlas of the human fetal nasal region from 10 fetuses spanning post-conceptional weeks 7 to 12. Using snRNA-seq, we identify 34 cell types in the fetal nasal region and reconstruct lineage trajectories in the olfactory epithelium (OE). Spatial analyses uncover regionalized OE organization, dorsoventral gradients of olfactory receptors (ORs) and high-resolution OR expression. Combining snRNAseq and MERFISH, we reveal early onset of OR gene expression, with most olfactory sensory neurons (OSNs) following the one neuron-one receptor rule already during the first trimester. This work provides the first integrated molecular and spatial framework of early human nasal region and it provides an unprecedented molecular blueprint of the human olfactory system development. In addition, this study offers a foundational resource for embryologists and developmental biologists investigating sensory neurogenesis, epithelial patterning, and congenital diseases.

  Study EGAS50000001203

• Whole-genome methylation profiling of menstrual stem cells identifies novel biomarkers for endometriosis

  Endometriosis is a highly prevalent condition that remains underdiagnosed and poorly managed due to the lack of clinically validated biomarkers and limited understanding of its underlying pathophysiology. Although menstrual blood–derived stem cells (MenSCs) have been implicated in disease pathogenesis, their potential as a diagnostic source has not been systematically explored. We conducted a clinical study including 42 women (19 with endometriosis and 23 controls) to evaluate whether DNA methylation profiles from freshly isolated MenSCs could identify disease-specific biomarkers. Whole-genome DNA methylation sequencing identified differentially methylated regions (DMRs) enriched in genes associated with key features of endometriosis, including inflammation, tissue remodelling, and developmental processes. These DMRs robustly discriminated cases from controls, independently of technical and clinical confounders. Machine learning models trained and validated on these regions demonstrated strong diagnostic performance (specificity 83%, sensitivity 79%). Integration with an independent single-cell RNA sequencing dataset suggested that the identified DMRs may influence gene expression, supporting their biological relevance. Collectively, these findings highlight MenSC DNA methylation profiling as a promising, non-invasive strategy for early diagnosis and personalised management of endometriosis. (PubMed IDs will be provided once the manuscript is published.)

  Study EGAS50000001640

• Non-canonical NF-κB signaling skews B cells away from germinal center to low-affinity effector fate

  The non-canonical (NC) NF-κB pathway is a central regulator of B-cell development and survival, yet its role during B-cell mediated adaptive responses remains uncharted. Here, we uncover an unanticipated regulatory function for NC NF-κB in shaping humoral immunity. High-resolution mapping of NC NF-kB status in B cells from human lymphoid tissues revealed low activity in germinal center B-cells (GCB) than effector subsets, suggesting a selective attenuation. Consistently, mice model harboring B-cells with enforced NC NF-κB activation revealed near- complete loss of GCB and high-affinity antibodies, following vaccination. GCB loss was not driven by intrinsic defects in formation or expansion, but by an immune-niche dysregulation marked by regulatory T-cell expansion. Mechanistically, NC NF-κB activation promoted B-cell- derived IL-10, and blocking this axis normalize the niche and rescue GCB. Concurrently, NC NF- κB signaling boosted low-affinity antibody responses by transcriptionally priming follicular B- cells for plasmacytic differentiation. Together, these findings position NC NF-κB as a tunable rheostat linking B-cell signaling to immune-niche control and humoral immunity.

  Study EGAS50000001646

• Submitters and requesters Statistics

  EGA Statistics Bibliography Growth Community Archive Distribution Catalogue Here, we show the amount and global distribution of submitters and requesters Requester accounts created by country The below figure represents the number of requester accounts by country. In most of the cases cases, country was inferred from the requester email domain. Those requesters whose email domain could not be assigned to a country such as @gmail.com, are not included. Note that a requester is a user that has had a granted permission at some point in time. Requester Accounts Created by Country Download as PNG Download as JPG Download as WEBP Download as CSV map pie mapChart('requester-country', 'https://stats.ega-archive.org/requester/country') pieChart('requester-country-pie-chart', 'https://stats.ega-archive.org/requester/country','Requester accounts created by country') Requester accounts created by year The below figure represents the number of requester accounts by year. Note that a requester is a user that has had a granted permission at some point in time. Requester Accounts Created by Year Download as PNG Download as JPG Download as WEBP Download as CSV barChart('requester-year', 'https://stats.ega-archive.org/requester/year', ['Requester Accounts', 'Years'], true) Industry requester accounts The below figure represents the top 15 industry requesting companies. Industry Requester Accounts Download as PNG Download as JPG Download as WEBP Download as CSV pieChart('industry-requester-accounts', 'https://stats.ega-archive.org/industry/requester/accounts', 'Industry Requester Accounts') Submitter accounts created by country The below figure represents the number of submitter accounts by country. Location was infered from the requester email domain. Submitter Accounts Created by Country Download as PNG Download as JPG Download as WEBP Download as CSV map pie mapChart('submitter-accounts', 'https://stats.ega-archive.org/submitter/accounts/country') pieChart('submitter-accounts-pie-chart', 'https://stats.ega-archive.org/submitter/accounts/country','Submitter accounts created by country') Submitter accounts created by year The below figure represents the number of submitter accounts by year. Submitter Accounts Created by Year Download as PNG Download as JPG Download as WEBP Download as CSV barChart('submitter-users-year', 'https://stats.ega-archive.org/submitter/users/year', ['Submitter Accounts', 'Years'], true) DAC Members by Country The below figure represents the number of DAC members by country. DAC Members by Country Download as PNG Download as JPG Download as WEBP Download as CSV map pie mapChart('dac-members-by-country', 'https://stats.ega-archive.org/dac_member/country') pieChart('dac-members-by-country-pie-chart', 'https://stats.ega-archive.org/dac_member/country','DAC members by country')

  Documentation about/statistics/community

• Adaptive Therapy Exploits Fitness Deficits in Chemotherapy-Resistant Ovarian Cancer to Achieve Long-Term Tumor Control Open Access

  Drug resistance results in poor outcomes for cancer patients. Adaptive therapy is a potential strategy to address drug resistance that exploits competitive interactions between sensitive and resistant subclones. Here, we showed that adapting carboplatin dose according to tumor response (adaptive therapy) significantly prolonged survival of murine ovarian cancer models compared to standard carboplatin dosing, without increasing mean daily drug dose or toxicity. Platinum resistant ovarian cancer cells exhibited diminished fitness when drug was absent in vitro and in vivo, which caused selective decline of resistant populations due to reduced proliferation and increased apoptosis. Conversely fitter, sensitive cells re-grew when drug was withdrawn. Using a bioinformatics pipeline that exploits copy number changes to quantify the emergence of treatment resistance, analysis of cell‐free DNA obtained longitudinally from ovarian cancer patients during treatment showed subclonal selection through therapy, and measurements of resistant population growth correlated strongly with disease burden. These preclinical findings pave the way for future clinical testing of personalized adaptive therapy regimens tailored to the evolution of carboplatin resistance in individual ovarian cancer patients.

  Study EGAS50000001142

• Whole Genome Sequencing of Neuroblastoma

  Neuroblastoma is a childhood tumour of the peripheral sympathetic nervous system. The pathogenesis has for a long time been quite enigmatic, as only very few gene defects were identified in this often lethal tumour1. Frequently detected gene alterations are limited to MYCN amplification (20%) and ALK activations (7%)2-5. Here we present a whole genome sequence analysis of 87 neuroblastoma of all stages. Only few recurrent amino-acid changing mutations were found. In contrast, analysis of structural defects identified in 18% of high stage neuroblastoma a local shredding of chromosomes, known as chromothripsis6. These tumours are associated with a poor outcome. Structural alterations recurrently affected ODZ3, PTPRD and CSMD1, which are involved in neuronal growth cone stabilization7-9. In addition, ATRX, TIAM1 and a series of regulators of the Rac Rho pathway were mutated, further implicating defects in neuritogenesis in neuroblastoma. Most tumors with defects in these genes were aggressive high stage neuroblastoma, but did not carry MYCN amplifications. The genomic landscape of neuroblastoma therefore reveals two novel molecular defects, chromotripsis and neuritogenesis gene alterations, which frequently occur in high risk tumors.

  Study EGAS00001000222

• Gut microbiome dynamics unravelled with metagenomics sequencing

  It has been already several years since we discovered that microbial therapeutics like Faecal Microbiota Transplantation can be used to treat Clostridium difficile infection (CDI) and Inflammatory bowel disease (IBD) with great success. The concept has challenged us to test the limits of our knowledge with the perspective to improve health in the most diverse conditions like Parkinson disease, depression, autism and obesity, just to mention a few. Note: the figure belongs to the original paper and is owned by the authors. Allyson L.Bird and colleagues from Genentech Inc. (San Francisco) recently published an interesting study building up on the key concept that “As microbial therapeutics are increasingly being tested in diverse patient populations, it is essential to understand the host and environmental factors influencing the microbiome.” The authors of this work analysed 1359 gut microbiome samples using shotgun metagenomics sequencing.  The results of the comparison done between microbiota composition across the different groups showed that biological sex is the strongest driver of diversity, while aging (age range studied: 20–69) appears to be the second factor. As smartly outlined in the authors’ own graphical abstract (figure to the left), Bacteroidota species consistently increased with donors' age while Actinobacteriota species, including Bifidobacterium, decreased. This study also integrated data from oncological patients to draw cancer-related evidences. After correcting the biases for technical and geographic differences, they found that cancer patients show significantly altered gut bacterial communities compared with their healthy counterparts. Eg. E.coli was found to be increased in cancer patients together with Bacteroidota/Firmicutes ratios. By comparing sequential samples taken from the same donors (413) within short-distance time points, the authors found that in the short term, intra-individual differences were less than the inter-individual ones, indicating a certain community stability being the norm. One of the pivotal points from this interesting study is that the data proceeds from extremely well-defined and selected donors. Specifically, the data for 949 selected healthy donors came from the Milieu Intérieur (MI) Consortium and extensive metadata, including demographic variables, serological measures, dietary information, and systemic immune profiles, are available. This precise stratification allowed for the robust detection of smaller changes, compared to former works. All sequencing data has been deposited in the EGA archive under the accession number EGAS00001004437. Enjoy this amazing open access paper published in the Experimental Journal of Medicine on January 2021.

  Blog gut-microbiome-dynamics-unravelled

• Whole-genome sequencing of two probands with hereditary spastic paraplegia reveals novel splice-donor region mutation and known pathogenic mutation in SPG11

  Hereditary spastic paraplegias (SPG) are a group of heterogeneous neurodegenerative disorders, which are often presented with overlapping phenotypes such as progressive paraparesis and spasticity. To assist the diagnosis of SPG subtypes, next-generation sequencing is often used to provide supporting evidence. In this study, we report the case of two probands from the same family with SPG symptoms, including bilateral lower limbs weakness, unsteady gait, cognitive decline, dysarthria and slurring of speech since age of 14. Subsequent whole-genome sequencing revealed that the patients are compound heterozygous for mutations in SPG11 gene, including the paternally inherited c.6856C>T (p.Arg2286*) mutation and the novel maternally inherited c.2316+5G>A splice donor region mutation. Mutations in SPG11 are the common cause of autosomal recessive Spastic Paraplegia type 11. According to the ClinVar database, there are already 101 reported pathogenic mutations in SPG11 that are associated with SPG. To our knowledge, this is the first report of SPG11 mutations in our local population. Novel splice mutation identified in this study would enrich the catalog of SPG11 mutations, potentially leading to better genetic diagnosis of SPG.

  Study EGAS00001001849

• Novel regional age-associated DNA methylation changes within human common disease-associated loci

  Background: Advancing age progressively impacts on risk and severity of chronic disease. It also modifies the epigenome, with changes in DNA methylation, due to both random drift and variation within specific functional loci. Results: In a discovery set of 2,238 peripheral-blood genome-wide DNA methylomes (MeDIP-seq) aged 19-82 years, we identified 71 age-associated Differentially Methylated Regions (a-DMRs, p < 1.85 x 10-8) within the Linkage Disequilibrium (LD) blocks of the NIH Catalogue of published GWAS SNPs. This included 52 novel regions, 29 within loci not covered by 450k or 27k Illumina array, and with marked enrichment for Poised Promoters and Enhancers across multiple cell types. In a replication set of 2,084 DNA methylomes, 95.7% of the a-DMRs showed the same direction of ageing effect, with 80.3% and 53.3% replicated to p < 0.05 and p < 1.85 x 10-8, respectively. Conclusion: By analysing the functionally enriched disease and trait-associated regions of the human genome, we identified novel epigenetic ageing changes, which could be useful biomarkers or provide mechanistic insights into age-related common diseases.

  Study EGAS00001001910

• Comprehensive molecular profiling identifies novel genetic drivers and subtypes underlying medulloblastoma

  Current therapies for medulloblastoma (MB), a highly malignant childhood brain tumor, impose debilitating effects on the developing child, warranting the development of molecularly targeted treatments with reduced toxicities. Prior studies have failed to disclose the full spectrum of driver genes and molecular processes operative in MB or adequately explain heterogeneity among molecular subgroups. Herein, we detail the somatic landscape across 491 sequenced MBs and molecular heterogeneity amongst 1,256 epigenetically analyzed cases, identifying subgroup-specific mutational signatures, driver genes, and pathway alterations including previously unappreciated actionable therapeutic targets. Integrative approaches for assigning driver genes to subgroups explained >65-70% of Group 3 and Group 4, doubling previous knowledge. Novel subtypes underlying Group 3 and Group 4 were differentially enriched for specific driver events, including hotspot in-frame indels targeting KBTBD4 and ‘enhancer hijacking’ driving activation of PRDM6. Thus, application of integrative genomics to an unprecedented cohort of clinical samples derived from a single childhood cancer entity has disclosed a series of new cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for treating MB.

  Study EGAS00001001953

• WTCCC2 Visceral Leishmaniasis (VL) samples

  A WTCCC2 project genome-wide association study for visceral leishmaniasis (VL) in individuals from India, Brazil and Sudan, genotyped on the custom Illumina 670k array. The WTCCC2 analysis of the Brazilian and Indian samples is described in Fakiola et al. [Nat Genet. 2013 Feb;45(2):208-13].It should be noted that due to expected family structure in the data, normal analyses of these data should include an estimation of the relatedness between the samples. For more details about sample collection for the project please refer to the Methods section of the paper above. The samples from India were all collected from Bihar state in northeastern India. The Brazilian samples were collected from a wide area of northern Brazil, and managed at two laboratories in Natal and Belem. For more details on the geographic distribution of the Brazil samples please see Jamieson et al. [Genes Immun. 2007 Jan;8(1):84-90]. We have not generated qc_passed or info files for the Sudanese data, a preliminary scan of the data indicated that there may be DNA quality issues with these samples.

  Study EGAS00001000773

• Global Anaplastic Thyroid Cancer Initiative

  Thyroid cancer is the most common endocrine malignancy. Most thyroid cancers are of the well-differentiated (non-aggressive) phenotype, and are almost always cured with standard treatments. In contrast, anaplastic thyroid cancer (ATC) is rare, accounting for only 1% of thyroid cancers, however it is perhaps the most lethal human malignancy with an average survival of 3 to 6 months. ATC presents with dramatic and rapid onset of airway and esophageal blockage with frequent spread to the lungs. Interestingly, 21–79% of ATCs have coexisting areas or a previous history of well-differentiated cancers. Although this suggests a progression from well-differentiated cancers to ATC, to date there is little molecular conformation that such a progression exists versus de novo generation of ATC. As ATC is extremely rare and most cases are inoperable, there are a paucity of tissue samples for study at any one center. The Global Anaplastic Thyroid Cancer Initiative (GATCI) aims to unite international institutions in order to pool samples (ATC with or without a paired well-differentiated component) for a comprehensive analysis of the genomic landscape of this disease.

  Study EGAS00001002234

• Field_effect_of_healthy_and_diseased_livers_WGS

  Recent work in the Campbell group has revealed somatic mutations present in normal, non-cancerous human skin. A subset of the mutations conferred selective advantages to the host cells, leading to clonal expansions and raising the risk for future cancer development. Capturing such somatic mutations in normal tissue is important to advance our understanding about carcinogenesis and could provide prospective medical insights. In this project, our goal is to detect somatic mutations in normal (pre-cancerous) liver tissue. Using Laser Microdissection technology, we will dissect individual liver lobules from patient samples and submit these to sequencing. For each patient sample, we aim to sequence multiple lobules to characterise the mutagenic burden. Samples will be taken from patients with different liver disease aetiologies, including alcoholism and obesity, with a view on distinguishing the prevalent mutation types occurring in each disease context. We will perform both whole genome and targeted sequencing, initially using the WTSI cancer panel. Later we aim to use a novel bait set that captures both cancer genes as well as genes relevant to the non-cancerous samples (ie. genes implicated in hereditary disorders, immune sequences).

  Study EGAS00001002413

• Cryptic Relatedness in the Singapore Living Biobank Project

  The Singapore Living Biobank is a collection of healthy population-based Chinese and Malay individuals, for the purpose of phenotype recall study of high-impact variant carriers. These individuals are sampled from two studies: Multi-Ethnic Cohort (MEC), and the Singapore Health 2012 (SH2012). The MEC is a population-based cohort initiated in 2007 to investigate the genetic and lifestyle factors that affect the risk of developing chronic diseases such as diabetes and cardiovascular outcomes in the three ethnic groups (Chinese, Malays, and Indians). The SH2012 study is a population-based cross-sectional survey conducted in Singapore between 2012 and 2013, with over-sampling of Malays and Indians. Participants in MEC and SH2012 completed a similar set of questionnaire components, health examination, and biochemisty panels. Description of the MEC and SH2012 studies can be found at http://blog.nus.edu.sg/sphs/. We generated whole-exome sequencing data and Illumina OmniExpress array genotyping data for 1,299 Chinese and 1,229 Malays from the Singapore Living Biobank. This study includes a subset of 762 individuals that were found to be closely related (≤3rd degree), including 263 Chinese and 499 Malays.

  Study EGAS00001002619

• Duplexseq_of_the_interstrand_crosslinks_WGS

  One of the most dangerous forms of DNA damage are interstrand crosslinks (ICLs), which covalently crosslink the two strands of the DNA double helix. The repair of these lesions is crucial for cellular survival due to their ability to block transcription and DNA replication. Initially, the major pathway that has been described in ICL repair involves a network of 22 genes that are mutated in a severe human genetic disease known as Fanconi Anemia (FA). Using synthetic lethality screens in the near-haploid human HAP1 cell line, we recently identified two potentially novel regulators of ICL repair, C1orf112 and THAP12. Loss of C1orf112 and THAP12 causes hypersensitivity to ICL-inducing DNA damaging agents, such as Mitomycin C (MMC). Additionally, C1orf112-depleted cells show elevated levels of micronuclei and accumulation of DNA damage in S-phase. To better understand how C1orf112 and THAP12 mediate the repair of ICLs, we want to perform mutational signature analysis, using the BotSeq method. Therefore, WT, C1orf112 and THAP12 knockout cells were cultured in vehicle or MMC treated conditions for 10 days and the genomic DNA was isolated. FANCA and FANCD2 knockout cells are taken along as controls in this experimental setting.

  Study EGAS00001006545

• SHANK2 mutations associated with Autism Spectrum Disorder cause hyperconnectivity of human neurons

  Heterozygous loss-of-function mutations in the synaptic scaffolding gene SHANK2 are strongly associated with autism spectrum disorder (ASD). To investigate their effect on synaptic connectivity, we generated cortical neurons from induced pluripotent stem cells (iPSC) derived from neurotypic and ASD-affected donors. We developed Sparse coculture for Connectivity (SparCon) assays where SHANK2 and control neurons were differentially labeled and sparsely seeded together on a lawn of unlabeled control neurons. We observed striking increases in dendrite length, dendrite complexity, total synapse number, and frequency of spontaneous excitatory postsynaptic currents. These findings were phenocopied in gene-edited homozygous SHANK2 knockout cells and rescued by gene correction of an ASD SHANK2 muation, supporting a role for SHANK2 as a regulator of connectivity in developing human neurons. Dendrite length increases were exacerbated by IGF1, TG003, or BDNF, and suppressed by DHPG treatment. The transcriptome in these isogenic SHANK2 neurons was deeply perturbed in synapse, plasticity, and neuronal morphogenesis gene sets and ASD gene modules, and activity-dependent dendrite extension was impaired. Our unexpected findings provide evidence for hyperconnectivity and profoundly altered transcriptome in SHANK2 neurons derived from ASD subjects.

  Study EGAS00001003436

• Pediatric Whole Genome Sequencing Diagnostic Utility

  The standard of care for first-tier clinical investigation of the etiology of congenital malformations and neurodevelopmental disorders is chromosome microarray analysis (CMA) for copy number variations (CNVs), often followed by gene(s)-specific sequencing searching for smaller insertion-deletions (indels) and single nucleotide variant (SNV) mutations. We compared diagnostic rate of WGS to CMA and targeted gene testing of 100 patients referred to The Hospital for Sick Children Genetics clinic in 2014. WGS identified genetic variants meeting clinical diagnostic criteria in 34% of cases, representing a 4-fold increase in diagnostic rate over CMA (8%) (p-value = 1.42e-05) alone and >2-fold increase in CMA plus targeted gene sequencing (13%) (p-value = 0.0009). WGS identified all rare clinically significant CNVs that were detected by CMA. In an additional 26 patients, WGS revealed indel and missense mutations presenting in a dominant (63%) or a recessive (37%) manner. We found four subjects with mutations in at least two genes associated with distinct genetic disorders, including two cases harboring a pathogenic CNV and SNV. When considering medically actionable secondary findings in addition to primary WGS findings, 38% of patients would benefit from genetic counselling.

  Study EGAS00001001623

• Spatial multi-omic map of human myocardial infarction

  Myocardial infarction is a leading cause of mortality worldwide 1 . While advances have been made in acute treatment, an incomplete understanding of remodelling processes has limited the effectiveness of therapies to reduce late-stage mortality 2 . Here, we generate an integrative high- resolution map of human cardiac remodelling after myocardial infarction using single-cell gene expression, chromatin accessibility, and spatial transcriptomic profiling of multiple physiological zones at distinct time points in myocardium from myocardial infarction and control patients. Multimodal data integration allowed us to evaluate cardiac cell-type compositions at increased resolution, yielding insights into changes of the cardiac transcriptome and epigenome through the identification of distinct tissue structures of injury, repair and remodelling. We identified and validated disease-specific cardiac cell-states of major cell-types and analysed them in their spatial context, evaluating their dependency on other cell-types. Our data elucidates molecular principles of human myocardial tissue organisation, recapitulating a gradual cardiomyocyte and myeloid continuum following ischemic injury. In total, our study provides an integrative molecular map of human myocardial infarction and represents an essential reference for the field and paves the way for advanced mechanistic and therapeutic studies of cardiac disease.

  Study EGAS00001006330

• scRNA-seq of HSPC treated with gemcitabine and carbplatin

  The study used cryopreserved hematopoietic stem and progenitor cells (HSPCs) harvested from one patient with chronic myelogenous leukemia (CML). The cells were thawed and acclimatized for 48 hours in StemMACS HSC Expansion Media XF (Miltenyi Biotec) supplemented with StemMACS HSC Expansion Cocktail (Miltenyi Biotec) at 37°C in a humidified atmosphere containing 5% CO2. 4 samples of 2000000 HSPCs in 2 ml media each were then prepared and treated for 24 hours: carboplatin-high: 150 µg/ml carboplatin carboplatin-low: 18.75 µg/ml carboplatin gemcitabine: 25 ng/ml gemcitabine control: no drug Single-cells were subsequently extracted using the ddSEQ™ Single-Cell Isolator (Bio-Rad) and later sequenced using the SureCell™ Whole Transcriptome Analysis 3' Library Prep Kit (Illumina) on the NextSeq 500 System (Illumina) using the NextSeq 500/550 High Output Kit v2.5 150 Cycles (Illumina) all following the manufacturer's instructions. FASTQ-files were then processed following the ddSeeker (Romagnoli, D., et al., BMC Genomics 2018, doi:10.1186/s12864-018-5249-x) and Drop-seq (Macosko, E.Z., et al., Cell 2015, doi:10.1016/j.cell.2015.05.002) protocols for processing scRNA-seq data to yield the final digital gene expression for each cell.

  Study EGAS00001004381

• Functional single-cell characterization of immune aplastic anemia shows convergence of NK and NK-like CD8+ T cells with disease-associated TCR signature

  Immune aplastic anemia (AA) is a life-threatening bone marrow (BM) failure disorder driven by an autoimmune T cell attack against the hematopoietic stem and progenitor cells (HSPCs). However, the autoantigen targets and the role of other immune cells are unknown. Here, by analyzing an international cohort of 156 AA patients with scRNA+TCRαβ-seq, TCRβ-seq, flow cytometry, and plasma cytokine profiling with comparisons to healthy controls and patients with other hematological disorders, we identify NK cells and CD8+ TEMRA cells expressing NK receptors with AA-associated TCRβ-seq motifs as the most dysregulated immune cell populations in AA BM. Functional co-culture experiments with scRNA-TCRαβ-seq readout utilizing primary HSPCs and immune cells provide evidence that NK cellscannot kill HSPCs alone, but sensitize them to CD8+ T cell-mediated killing. Furthermore, HSPCs induce strong activation of T cell clones with CD8+ TEMRA NK-like phenotype and AA-associated TCR motifs. Our results reveal the convergent evolution of innate and adaptive immune cells in AA, where NK cells support CD8+ T cell-mediated autoimmunity.

  Study EGAS00001007602

• Distinct portrayal of lesions in synchronous multifocal lung adenocarcinoma revealed by genome sequencing

  Distinguishing multiple primary lung cancers in the synchronous multifocal intrapulmonary lesions has important significance on clinical staging and therapeutic decision. To investigate genomic aberration profiles, we applied whole genome and whole exome sequencing, and microarray-based comparative genomic hybridization on 15 intrapulmonary tumors derived from six patients with synchronous multifocal lung cancers having similar histological diagnosis. Any pair of intrapulmonary tumors in a single patient, which shared the identical genetic background and environment, showed an extinctive heterogeneity between each other. Phylogenetic relationship analysis indicated an independently branched evolution among all the tumors, suggesting they were multiple primary lung cancers. EGFR or KRAS mutations were found in 7 or 3 out of the 15 tumors, from 3 or 2 patients, respectively. Somatic mutational heterogeneity of these two genes in a single patient was also observed. Our analysis indicates genomic aberration profiling is valuable for identification of multiple primary lung cancer, especially when high histopathological concordance was observed between lesions. We also suggest a thoroughly molecular diagnosis against therapeutic target genes should be taken for each accessible nodule before making a plan for adjuvant therapy.

  Study EGAS00001001572

• SINGLE-CELL RNA SEQUENCING Single-cell RNA sequencing was performed on 13 ‘mild-moderate’ and 10 ‘critical’ COVID19 PBMC samples

  Summary Epidemiological and clinical reports have indicated that SARS-CoV-2 virulence hinges upon the triggering of an aberrant host immune response, more so than on direct virus-induced cellular damage. To elucidate the immunopathology underlying COVID-19 severity, we performed cytokine and multiplex immune profiling in mild-moderate and critically ill COVID-19 patients. Hypercytokinemia in COVID-19 differed from the IFN-γ-driven cytokine storm in macrophage activation syndrome, and was more pronounced in critical versus mild-moderate COVID-19. Systems modelling of cytokine levels paired with deep-immune profiling showed that classical monocytes drive this hyper-inflammatory phenotype and that a reduction in T-lymphocytes correlates with disease severity, with CD8+ cells being disproportionately affected. Expression of antigen presenting machinery was also reduced in critical disease. Furthermore, we found that neutrophils contributed to disease severity and local tissue damage by amplification of hypercytokinemia and the formation of neutrophil extracellular traps. Together these findings suggest a myeloid-driven immunopathology, in which hyperactivated neutrophils and an ineffective adaptive immune system act as mediators of COVID-19 disease severity.

  Study EGAS00001005039

• SmMIP-tools:a computational toolset for processing and analysis of single-molecule molecular inversion probes derived data

  Single-molecule molecular inversion probes (smMIPs) provides a modular and cost-effective platform for high-multiplex targeted next-generation sequencing (NGS). Nevertheless, translating the raw smMIP-derived sequencing data into accurate and meaningful information currently requires proficient computational skills and a large amount of computational work, prohibiting wide-scale adoption of smMIP-based technologies. To enable easy, efficient, and accurate interrogation of smMIP-derived data, we developed SmMIP-tools, a computational toolset that combines the critical analytic steps for smMIP data interpretation into a single computational pipeline. Here, we describe in detail two major components of the software. The first is a read processing tool that performs quality control steps, generates read-smMIP linkages and retrieves molecular tags. The second is an error-aware variant caller capable of detecting single nucleotide variants (SNVs) and short insertions and deletions (indels). Using a cell-line DNA dilution series and a cohort of blood cancer patients, we benchmarked SmMIP-tools and evaluated its performance against clinical sequencing reports. We anticipate that SmMIP-tools will increase accessibility to smMIP-technology, enabling cost-effective genetic research to push personalized medicine forward.

  Study EGAS00001005359

• Transcriptomes of human CD4+ T lymphocytes - Metabolic project

  Human CD4+ T cells are essential mediators of immune responses. By altering the mitochondrial and metabolic states, we defined metabolic requirements of human CD4+ T cells for in vitro activation, expansion, and effector function. T cell activation and proliferation were reduced by inhibiting oxidative phosphorylation, while early cytokine production was maintained by either OXPHOS or glycolytic activity. Glucose deprivation in the presence of mild mitochondrial stress markedly reduced all three T cell functions, contrasting the exposure to resveratrol, an antioxidant and sirtuin-1 activator, which specifically inhibited cytokine production and T cell proliferation, but not T cell activation. Conditions that inhibited T cell activation were associated with the downregulation of 2′,5′-oligoadenylate synthetase genes via interferon response pathways. Our findings indicate that T cell function is grossly impaired by stressors combined with nutrient deprivation, suggesting that correcting nutrient availability, metabolic stress and/or the function of T cells in these conditions will improve the efficacy of T cell-based therapies.

  Study EGAS00001005565

• Single-nucleus transcriptomic profiling of aging Down Syndrome brains

  Down syndrome (DS) – trisomy of human chromosome 21 (HSA21) – is characterized by lifelong cognitive impairments and development of neuropathological hallmarks of Alzheimer’s disease (AD). The cellular and molecular modifications responsible for these effects are not understood. Here we performed single-nucleus RNA-sequencing (snRNA-seq) employing both short (Illumina) and long-read (Pacific Biosciences) sequencing technologies on a total of 29 DS and non-DS control prefrontal cortex samples. In DS, the ratio of inhibitory-to-excitatory neurons was significantly increased, which was not observed in previous reports examining sporadic AD. DS microglial transcriptomes displayed AD-related aging and activation signatures in advance of AD neuropathology, with increased microglial expression of C1q complement genes (associated with dendritic pruning) and the HSA21 transcription factor gene RUNX1. Long-read sequencing detected vast RNA isoform diversity within and amongst specific cell types including numerous novel sequences that differed between DS and normal brains. Notably, over 8,000 genes produced RNAs containing intra-exonic junctions (IEJs), including APP that had previously been associated with somatic gene recombination. These and related results illuminate large-scale cellular and transcriptomic alterations as novel features of the aging DS brain.

  Study EGAS00001005691

• Increased mutation accumulation during fetal development in Down syndrome

  Children with Down syndrome (DS) show a high incidence of leukemia, which is preceded by transient myeloproliferative disorder (TMD) in neonates. TMD is characterized by mutations in GATA1 that accumulate during fetal development, indicating that additional mutagenesis is required for leukemogenesis. Abnormalities in chromosome numbers can result in a mutator phenotype. However, it is unknown whether a trisomy of chromosome 21 (T21) results in enhanced mutation accumulation during human development. Here, we determined the mutational consequences of a constitutive T21 in DS fetuses compared to disomy 21 cells (D21) of D21 fetuses. We found that D21 fetal HPSCs had a higher somatic mutation rate compared to healthy adults HSPCs. Moreover, we observed an increase in somatic mutation load in T21 cells of DS fetuses, which was already apparent during the first cell divisions of embryogenesis. The mutation load could be explained by processes that are active during normal fetal development and were similar to those observed in TMD blast cells. Our study shows that T21 causes an increase in somatic load in fetal cells of DS fetuses, which may contribute to leukemia development early in life.

  Study EGAS00001003982

• The DNA methylation landscape of glioblastoma disease progression shows extensive heterogeneity in time and space

  Glioblastoma is characterized by widespread genetic and transcriptional heterogeneity, yet little is known about the role of the epigenome in glioblastoma disease progression. Here, we present genome-scale maps of DNA methylation in matched primary and recurring glioblastoma tumors, using data from a highly annotated clinical cohort that was selected through a national patient registry. We demonstrate the feasibility of DNA methylation mapping in a large set of routinely collected FFPE samples, and we validate bisulfite sequencing as a multipurpose assay that allowed us to infer a range of different genetic, epigenetic, and transcriptional characteristics of the profiled tumor samples. On the basis of these data, we identified subtle differences between primary and recurring tumors, links between DNA methylation and the tumor microenvironment, and an association of epigenetic tumor heterogeneity with patient survival. In summary, this study establishes an open resource for dissecting DNA methylation heterogeneity in a genetically diverse and heterogeneous cancer, and it demonstrates the feasibility of integrating epigenomics, radiology, and digital pathology for a national cohort, thereby leveraging existing samples and data collected as part of routine clinical practice.

  Study EGAS00001002538

• Cross-tissue transcriptomic analysis of human secondary lymphoid organ residing ILC3 reveals a default quiescent state in the absence of inflammation

  A cross tissue transcriptional comparison of human ILC3 isolated from non-inflamed lymph nodes (LN) and spleens, inflamed tonsils and the peripheral blood (PB). Using uniform gating and high-purity cell-sorting, ILC3 were sorted as Lineage(CD3,CD19,CD34,CD14,CD94)-CD45intCD127+CRTH2-CD117+. From LN, spleen, tonsils and PB, NKp44- ILC3 were isolated and, in addition, NKp44+ ILC3 were sorted from spleen and tonsils. RNA quality was measured on a RNA 6000 Pico Kit (Agilent Technologies) using a 2100 Bioanalyzer (Agilent Technologies). 100 pg high-quality RNA was isolated and cDNA was generated using the SMART-Seq (v3) Ultra Low Input RNA Kit (Clontech) with 15 cycles of amplification. Amplified and Covaris-fragmented cDNA was quality-checked on a High Sensitivity DNA chip (Agilent Technologies) and further processed according the TruSeq Nano DNA Library Preparation Kit (Illumina). Data was processed as follows: SMARTer adapters were trimmed with cutadapt and the resulting sequences were aligned to the human RefSeq transcriptome using 14 TopHat2 (Kim et al., 2013). Normalization and quantification was performed using Cufflinks (Trapnell et al., 2012). FPKM counts were determined per gene with HTSeq-count.

  Study EGAS00001002636

• The identification of genetic vulnerabilities in head and neck cancers for the development of novel treatments.

  The global cancer burden is expected to double by the year 2030 and this will disproportionately affect the Asian continent. Head and neck squamous cell carcinoma (HNSCC) is endemic in Asia. It is the sixth most common cancer worldwide and the second most common in Malaysia, with approximately 600,000 new cases per year globally. Longterm survival remains less than 50%, underscoring the need for better therapeutic strategies. In this proposal we aim to bring new technological advancements in the identification of critical genes in HNSCC and the identification of compounds for the development of new therapies in HNSCC. We propose to genomically characterise 16 unique Asian H&N cancer cell lines, to classify sub groups where identification of essential genes that can be targeted for tumor control can be applied in a more rational way. Further, we propose to use zebrafish, an established drug discovery platform to identify promising pathway-specific compounds that could be further studied using CRISPR/Cas9 lethality screens and biological models to better identify novel therapeutic strategies to combat this major health problem in Asia.

  Study EGAS00001002682

• Pan-neuroblastoma analysis reveals age- and signature-associated driver alterations

  Neuroblastoma is a pediatric malignancy with heterogeneous clinical outcomes depending on age and other factors. To better understand neuroblastoma pathogenesis, here we analyze whole-genome, whole-exome and/or transcriptome sequencing data from 702 neuroblastoma samples, half of which are new data. Over 95% of samples harbor at least one recurrent driver alteration and most aberrations, including MYCN, ATRX, and TERT alterations, differ in frequency by age. MYCN alterations occur at median 2.3 years of age, TERT at 3.8 years, and ATRX at 5.6 years, suggesting age-specific susceptibility to genetic alterations. Mutational signature analysis shows that reactive oxygen-species (ROS)-induced mutations are the most common cause of single nucleotide variant (SNV) drivers in neuroblastoma, including most ALK and Ras pathway SNVs. ROS-induced mutagenesis is both an early event, as it causes truncal variants, and continuous as it also induces relapse-specific mutations. ROS-induced mutations are more abundant in neuroblastomas with MYCN amplification, 17q gains, and increased expression of mitochondrial ribosome genes. We observe recurrent FGFR1 N546K or internal tandem duplication variants in 6 patients and ALK N-terminal structural alterations in 5 samples, identifying additional patients amenable to precision therapy.

  Study EGAS00001003931

• An alternative splicing modulator decreases mutant HTT and improves the molecular fingerprint in Huntington’s disease patient neurons

  Huntington's disease (HD) is a neurodegenerative disorder caused by poly-Q expansion in the Huntingtin (HTT) protein. Here, we delineate elevated mutant HTT (mHTT) levels in patient-derived cells including fibroblasts and iPSC derived cortical neurons using a GLP approved HTT assay. HD patients’ fibroblasts and cortical neurons recapitulate aberrant alternative splicing as a molecular fingerprint of HD. Branaplam is a splicing modulator currently tested in a phase II study in HD (NCT05111249). The drug lowers total HTT (tHTT) and mHTT levels in fibroblasts, iPSC, cortical progenitors, and neurons in a dose dependent manner at an IC50 consistently below 10nm without inducing cellular toxicity. Branaplam promotes inclusion of non-annotated novel exons. Amongst, a 115bp frameshift-inducing exon in the HTT transcript in Branaplam treated cells from Ctrl and HD patients leading to a profound reduction of HTT RNA and protein levels. Importantly, Branaplam ameliorates aberrant alternative splicing in HD patients’ fibroblasts and cortical neurons. These findings highlight the applicability of splicing modulators in the treatment of CAG repeat disorders and decipher their molecular effects associated with the pharmacokinetic and -dynamic properties in patient-derived cellular models.

  Study EGAS00001006289

• Study on the proliferation history of colorectal adenomas

  Mismatch repair (MMR) deficient colorectal adenomas are composed of transformed cells that descend from a common founder and progressively accumulate genomic alterations but their proliferation history is still largely unknown. Here we rebuilt the proliferation history of four MMR deficient colorectal adenomas from male individuals by deep sequencing the X chromosome and mapping the somatic mutations that were progressively acquired during tumor growth. We developed a novel method to call high and low frequency mutations, to cluster them according to their frequencies and rebuild the proliferation trees directly from these mutation clusters using a recursive algorithm. The trees of all four lesions were formed of a dominant subclone that coexisted with other genetically heterogeneous subpopulations of cells. Despite this similar hierarchical organization, however, the growth dynamics varied among and within tumors, likely depending on a combination of tumor-specific genetic and environmental factors. Our study provides insights into the biological properties of individual MMR deficient colorectal adenomas that may influence their growth and also the response to therapy. Extended to other solid tumors, our approach could inform on the mechanisms of cancer progression, thus providing insight into the best choice of cancer treatment methods.

  Study EGAS00001000883

• UK10K_OBESITY_GS

  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 GS samples are part of the Obesity group and will undergo exome sequencing.

  Study EGAS00001000242

• The mission of the BIOS Consortium is to create a large-scale data infrastructure and to bring together BBMRI researchers focusing on integrative omics studies in Dutch Biobanks.

  The advent of the genome-wide association study (GWAS) led to the successful identification of thousands of variants that are robustly associated with complex disease phenotypes. Dutch biobanks played a substantial role in these discoveries. For most of these variants, however, the mechanisms through which they contribute to these phenotypes remain unknown. The BIOS Consortium applies a functional genomics approach that integrates genome-wide genetic data with data on the epigenome and transcriptome to elucidate these mechanisms. Over 4000 samples from BBMRI-NL biobanks with in-depth information on disease phenotypes and GWAS data are being enriched with RNA-sequencing and genome-wide DNA methylation data. The same is true for samples with whole-genome sequencing data from GoNL. This unique data infrastructure provides a powerful platform to evaluate key questions in integrative omics from establishing comprehensive eQTL and meQTL catalogues to linking molecular pathways across omics levels to phenotypic outcomes.

  Study EGAS00001001077

• Comprehensive characterization of cell-free tumor DNA in plasma and urine of patients with renal tumors

  Cell-free tumour-derived DNA (ctDNA) allows non-invasive monitoring of various cancers but its utility in renal cell cancer (RCC) has not been well established. A combination of untargeted and targeted sequencing methods, applied to two independent cohorts of patients (n=90) with various renal tumour subtypes, were used to determine ctDNA content in plasma and urine. Our data revealed lower plasma ctDNA levels in RCC relative to other cancers, with untargeted detection of ~33% for both cohorts. A highly sensitive personalised approach, applied to plasma and urine from select patients (n=22), improved detection to ~50%, including in patients with early stage and even benign lesions. A machine-learning based model, applied to untargeted data, predicted this detection, potentially offering a means of triaging patient samples for personalised analysis. We observed that plasma, and for the first time, urine ctDNA may better represent tumour heterogeneity than a single tissue biopsy. Longitudinal sampling of >200 plasma samples revealed that ctDNA can track disease course. These data highlight low ctDNA levels in RCC but indicate potential clinical utility provided improvement in isolation and detection approaches.

  Study EGAS00001003530

• Genome-wide cell-free DNA fragmentation in patients with cancer

  Cell-free DNA (cfDNA) in the blood provides a noninvasive diagnostic avenue for patients with cancer. However, characteristics of the origins and molecular features of cfDNA are poorly understood. We developed an approach to evaluate fragmentation patterns of cfDNA across the genome and found that cfDNA profiles of healthy individuals reflected nucleosomal patterns of white blood cells, while patients with cancer had altered fragmentation profiles. We applied this method to analyze fragmentation profiles of 236 patients with breast, colorectal, lung, ovarian, pancreatic, gastric, or bile duct cancers and 245 healthy individuals. A machine learning model incorporating genome-wide fragmentation features had sensitivities of detection ranging from 57% to >99% among the seven cancer types at 98% specificity, with an overall AUC of 0.94. Fragmentation profiles could be used to identify the tissue of origin of the cancers to a limited number of sites in 75% of cases. Combining our approach with mutation based cfDNA analyses detected 91% of cancer patients. The results of these analyses highlight important properties of cfDNA and provide a proof of principle approach for screening, early detection, and monitoring of human cancer.

  Study EGAS00001003611

• Genome‐wide postnatal changes in immunity following fetal inflammatory response

  The fetal inflammatory response (FIR) increases the risk of perinatal brain injury, particularly in extremely low gestational age newborns (ELGANs, < 28 weeks of gestation). One of the mechanisms contributing to such a risk is a postnatal intermittent or sustained systemic inflammation (ISSI) following FIR. The link between prenatal and postnatal systemic inflammation is supported by the presence of well‐established inflammatory biomarkers in the umbilical cord and peripheral blood. However, the extent of molecular changes contributing to this association is unknown. Using RNA sequencing and mass spectrometry proteomics, we profiled the transcriptome and proteome of archived neonatal dried blood spot (DBS) specimens from 21 ELGANs. Comparing FIR‐affected and unaffected ELGANs, we identified 782 gene and 27 protein expression changes of 50% magnitude or more, and an experiment‐wide significance level below 5% false discovery rate. These expression changes confirm the robust postnatal activation of the innate immune system in FIR‐affected ELGANs and reveal for the first time an impairment of their adaptive immunity. In turn, the altered pathways provide clues about the molecular mechanisms triggering ISSI after FIR, and the onset of perinatal brain injury.

  Study EGAS00001003635

• Genomic analysis of Smoothened inhibitor resistance in basal cell carcinoma

  Smoothened inhibitors are currently being investigated for the treatment of several cancers. Vismodegib is approved for the treatment of locally advanced and metastatic basal cell carcinoma (BCC). The majority of BCC patients treated with vismodegib experience significant clinical benefit, however, a small number of patients develop resistance. Knowledge of resistance mechanisms can generate predictive biomarkers and is critical for the design of additional therapeutic strategies aimed at circumventing resistance. To investigate mechanisms of resistance to vismodegib in BCC we sequenced and profiled biopsies from patients who initially responded to treatment and subsequently progressed on drug.We find that resistance is associated with re-activation of the Hedgehog (Hh) pathway through diverse mechanisms, including concurrent copy number alterations in the downstream Hh pathway components SUFU and GLI2, as well as novel SMO mutations that reduce sensitivity to vismodegib. The latter fall into two classes: those that activate SMO and those that directly affect the vismodegib-binding pocket. Additionally, we observe evidence of clonal selection of SMO mutations and intra-tumor heterogeneity of resistance mechanisms, which implies diverse and combinatorial strategies are required to overcome resistance.

  Study EGAS00001000845

• Germline mutations in the transcription factor IKZF5 cause thrombocytopenia.

  To identify novel causes of hereditary thrombocytopenia, we performed a genetic association analysis of whole-genome sequencing data from 13,037 individuals enrolled in the NIHR BioResource, including 233 cases with isolated thrombocytopenia. We found an association between rare variants in the transcription factor encoding gene IKZF5 and thrombocytopenia. We report five causal missense variants in or near IKZF5 zinc fingers (Znfs), of which two occurred de novo and three co-segregated in three pedigrees. A canonical DNA-Znf binding model predicts that three of the variants alter DNA recognition. Expression studies showed that chromatin binding was disrupted in mutant compared to wild-type (WT) IKZF5 and electron microscopy revealed a reduced quantity of alpha granules in normally sized platelets. Proplatelet formation (PPF) was reduced in megakaryocytes (MKs) from seven cases relative to six controls. Comparison of RNA-seq data from platelets, monocytes, neutrophils and CD4 T-cells from three cases and 14 healthy controls showed 1,194 differentially expressed genes in platelets but only four DEGs in each of the other blood cell types. In conclusion, IKZF5 is a novel transcriptional regulator of megakaryopoiesis and the eighth transcription factor associated with dominant thrombocytopenia in humans.

  Study EGAS00001003736

• RNA sequencing and Illumina 2.5M SNP array data collected from 675 commonly used human cancer cell lines.

  Tumor-derived cell lines have served as vital models to advance our understanding of oncogene function and therapeutic response1. Although substantial effort has been directed to defining the genomic constitution of cancer cell line panels2–4, the transcriptome – which represents the active program of a cell – remains understudied. Here, we describe RNA sequencing and SNP array analysis of 675 commonly used human cancer cell lines. We explore numerous transcriptome features including coding and non-coding gene expression, transcribed mutations, gene fusion and expression of non-human sequences. Aside from many known aberrations we find new surprising characteristics, including more than 2200 unique fusion gene pairs representing a vast, testable repertoire of oncogenic fusions, many of which have analogs found in primary human tumors. We show that a combination of multiple genome and transcriptome features in a novel pathway-based approach enhances prediction of response to various targeted therapeutics. Our results provide valuable new insights into these critical pre-clinical models and provide added context for interpreting the numerous studies that employ these widely used cell lines.

  Study EGAS00001000610

• Reference Exome Data for a Northern Brazilian population

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

  Study EGAS00001004112

• Oncoprint GSCCs

  Glioblastoma (GBM) remains the most malignant primary brain tumor, with a median survival rarely exceeding 2 years. Tumor heterogeneity and an immunosuppressive microenvironment are key factors contributing to the poor response rates of current therapeutic approaches. GBM-associated macrophages (GAMs) often exhibit immunosuppressive features that promote tumor progression. However, their dynamic interactions with GBM tumor cells remain poorly understood. Here, we used patient-derived GBM stem cell cultures and combined single-cell RNA sequencing of GAM-GBM co-cultures and real-time in vivo monitoring of GAM-GBM interactions in orthotopic zebrafish xenograft models to provide insight into the cellular, molecular, and spatial heterogeneity. Our analyses revealed substantial heterogeneity across GBM patients in GBM-induced GAM polarization and the ability to attract and activate GAMs – features that correlated with patient survival. Differential gene expression analysis, immunohistochemistry on original tumor samples, and knock-out experiments in zebrafish subsequently identified LGALS1 as a primary regulator of immunosuppression. Overall, our work highlights that GAM-GBM interactions can be studied in a clinically relevant way using co-cultures and avatar models, while offering new opportunities to identify promising immune-modulating targets.

  Study EGAS00001007481

• Spatial_transcriptome_analysis_of_Paediatric_Thymus

  Cell atlas of Paediatric Thymus Cell atlas of Paediatric Thymus - The Human Cell Atlas initiative is an ambitious global research effort which aims to describe every cell in the human body across all developmental stages to generate a reference map to accelerate progress in science and medicine particularly relating to human development and disease. The initiative is led by investigators based at the Wellcome Sanger Institute and the Broad Institute of MIT and Harvard and funding support from the Wellcome and NIH. We aim to interrogate the cellular composition and molecular regulators underpinning a cells identity through; genetic and genomic analysis, expression analysis, non-identifiable video imaging and cell line and organoid production, as well as, computational algorithms to predict cellular developmental trajectories and cell-cell interactions. We will subsequently validate these findings in situ using spatial genomic and imaging technologies e.g. immunohistochemistry, RNAscope, in situ sequencing and spatial transcriptomic analyses. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001004281

• Widespread DNA hypomethylation and differential gene expression in Turner syndrome

  Adults with 45,X monosomy (Turner syndrome) reflect a surviving minority since more than 99% of fetuses with 45,X monosomy die in utero. In adulthood 45,X monosomy is associated with increased morbidity and mortality, although strikingly heterogeneous with some individuals left untouched while others suffer from cardiovascular disease, autoimmune disease and infertility. The present study investigates the leukocyte DNAmethylation profile by using the 450K-Illumina Infinium assay and the leukocyte RNA-expression profile in 45,X monosomy compared with karyotypically normal female and male controls. We present results illustrating that genome wide X-chromosome RNA-expression profile, autosomal DNA-methylation profile, and the X-chromosome methylation profile clearly distinguish Turner syndrome from controls. Our results reveal genome wide hypomethylation with most differentially methylated positions showing a medium level of methylation. Contrary to previous studies, applying a single loci specific analysis at well-defined DNA loci, our results indicate that the hypomethylation extend to repetitive elements. We describe novel candidate genes that could be involved in comorbidity in TS and explain congenital urinary malformations (PRKX), premature ovarian failure (KDM6A), and aortic aneurysm formation (ZFYVE9 and TIMP1)

  Study EGAS00001002190

• Determination_of_cell_specific_regulatory_enhancers_in_hematopoetic_models

  The BLUEPRINT project is a large-scale project investigating epigenetic mechanisms involved in blood formation, in health and disease. The human variation workpackage (WP10) of the project seeks to characterize the effect of common sequence variation on the epigenome status of a cell. To do this, the project will use highly purified blood cells to minimise “experimental noise” and therefore enhance the power to discover modest effects. Two peripheral blood cell types, the CD14+CD16- monocyte (an important central orchestrator of adaptive immunity and a bridge between innate and adaptive immunity) and the CD65+CD9- neutrophilic granulocyte (the frontline cell for innate immunity) have been selected for this purpose. The two types of cells will be obtained at high purity from adult blood (AB) of 200 healthy males and females, respectively. Cells will be purified by using already validated and fully operational protocols that are based on density gradient centrifugation of the buffy coat obtained from whole blood, followed by magnetic bead-based purification using monoclonal antibodies against Cluster of Differentiation (CD) lineage-specific cell surface markers. This data set contains functional genomics data for gene expression and chromatin state.

  Study EGAS00001000586

• We performed whole-exome sequencing of 20 samples (10 actinic keratosis and 10 cutaneous squamous cell carcinoma) to investigate a potential relationship between DNA methylation-based subtypes and genetic mutation patterns (Rodriguez-Paredes et al., Nat Commun 2017)

  Cutaneous squamous cell carcinoma (cSCC) is the second most common skin cancer type and arises from keratinocytes. Most cSCC progress from a UV-induced precancerous lesion termed actinic keratosis (AK). Despite various efforts to characterize these lesions molecularly, the etiology of AK and its progression to cSCC remain only partially understood. Here we have used Infinium MethylationEPIC BeadChips to interrogate the DNA methylation status of about 850.000 CpGs in epidermal preparations from healthy skin, AK and cSCC. Importantly, we found that the premalignant AK samples displayed classical features of cancer methylomes and were highly similar to cSCC methylomes. Further analysis identified typical features of stem cell methylomes, such as a reduced DNA methylation age, non-CpG methylation and stem cell-related keratin and enhancer methylation patterns. Interestingly, this signature was detected only in one half of the AK and cSCC samples, while the other half showed keratin and enhancer methylation patterns that were more closely related to the control epidermis. These findings suggest the existence of two distinct subclasses of AK and cSCC that originate from distinct keratinocyte differentiation stages.

  Study EGAS00001002670

• Gut microbiome modulates response to anti PD1 immunotherapy in metastatic melanoma patients

  There is a growing appreciation of the role of the microbiome in cancer, and evidence in pre-clinical models that the gut microbiome may modulate responses to immune checkpoint blockade though this has not been well-characterized in patients. We analyzed the oral (n=86)and gut (n=43) 16S microbiome in melanoma patients on PD-1 blockade. Significant differences were noted in the diversity and composition of the gut microbiome between responders and non-responders in patients with a fecal microbiome sample, with significantly higher alpha diversity and relative abundance of Ruminococcaceae bacteria) in R. Metagenomic studies (n=25) revealed functional differences in gut bacteria in R including enrichment of anabolic pathways. Immune profiling demonstrated enhanced systemic and anti-tumor immunity in patients with a favorable gut microbiome as well as in germ-free mice receiving fecal transplant from responding patients. Together, these data have important implications for treatment with immune checkpoint blockade in cancer. The fastq files associated with this dataset are stored at ENA under the following accession numbers: PRJEB22894 (Fecal 16S), PRJEB22874 (Oral 16S), PRJEB22895 (Murine 16S), PRJEB22893 (Fecal WGS).

  Study EGAS00001002698

• Germline HAVCR2 mutations altering TIM-3 characterize subcutaneous panniculitis-like T-cell lymphomas with hemophagocytic lymphohistiocytic syndrome

  Sub-cutaneous panniculitis-like T-cell lymphomas (SPTCL), a non-Hodgkin lymphoma, can be associated with hemophagocytic lymphohistiocytosis (HLH), a life-threatening immune activation which adversely impacts survival1,2. T-cell-immunoglobulin mucin-3 (TIM-3) is a modulator of immune responses expressed on subgroups of T- and innate immune cells. We identify in ~60% of SPTCL cases germline, loss-of-function, missense variants in highly conserved residues of TIM-3, c.245A>G(p.Y82C) and c.291A>G(p.I97M), each with specific geographic distribution. Y82C-TIM-3 occurs on a potential founder chromosome in patients with East Asian and Polynesian ancestry, while I97M-TIM-3 occurs in Caucasians. Both variants induce protein misfolding and abrogate TIM-3’s plasma-membrane expression leading to persistent immune activation, increased production of inflammatory cytokines, including TNF-α and IL-1β, promoting HLH and SPTCL. Our findings highlight HLH/SPTCL as a new genetic entity and showcase TIM-3 mutations as the first causative genetic defect in SPTCL. While TIM-3-mutant HLH/SPTCL benefit from immunomodulation, therapeutic repression of the TIM-3 checkpoint may have adverse consequences.

  Study EGAS00001002765

• SINGLE_CELL_ANALYSIS_OF_IN_VITRO_ERYTHROPOIESIS

  We are developing a protocol to differentiate mouse and human induced pluripotent stem (IPS) and embryonic stem (ES) cells towards the haematopoietic pathway to generate erythrocytes in vitro. This system has many applications such as the study of the role of specific genes and human polymorphisms in infectious diseases such as malaria, as well as haematological diseases such as myelodysplastic syndrome. The nature of the in vitro differentiation process means that a heterogeneous population of cells is generated. In order to understand the types of cells produced with our protocol, we have performed a single cell analysis, which has the power to reveal the different populations of cells and their characteristics. For this, a cDNA library has been made that needs to be sequenced to obtain the gene expression profiles of the different cells. With this information we will be able to assess the quality of the differentiation protocol and improve it in order to produce better cells for the downstream applications.This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001000764

• H3Africa - Identification and characterization of novel hereditary neurological disease genes in Mali.

  Despite the vast diversity of its populations, genetic studies in Africa have been limited. African populations, Malians in particular, have a high rate of intra-ethnic and consanguineous marriage, resulting in increased prevalence of autosomal recessive diseases. Family-based genetic studies can be limited in developed countries due to small sibships. The average fertility rate in Mali is over 6 births per woman, offering a unique opportunity to find new disease genes or mutations that can then be studied in other populations. Hereditary neurological diseases are very debilitating diseases, and developing countries, particularly in Africa, pay a high price in terms of disability-adjusted life years. Although most are currently untreatable, increasing awareness about hereditary neurological disorders can reduce this burden. This study will help to identify and characterize novel hereditary neurological disease genes in Mali. These genes are likely to be important in the normal function of the nervous system and to have important pathophysiological implications for African and other populations. This study will also train physicians and students in the characterization of neurodegenerative diseases as well as in genetic technology and molecular biology, and build a suitable research environment that will retain them.

  Study EGAS00001003016

• Whole genome sequencing of Ewings Sarcoma

  Sarcomas are cancers of the bone and soft tissue often defined by their gene fusions. However, the timing, context, and processes by which these pathogenic fusions arise are unknown. We explored this in Ewing sarcoma, a cancer driven by EWSR1-ETS fusions, with very few cooperating mutations. Combining whole-genome sequencing with enhanced informatics, we found that the EWSR1-ETS fusion arose from striking rearrangement clusters in 42% of cases (52/124). Notably, these were organized in loops that universally contained the fusion at their center, while also weaving up to 18 genes together with it. We found the same pattern of rearrangements in three additional types of sarcoma. From these data, we define a new signature for sarcoma fusions that precedes other somatic changes, in the earliest replicating DNA of the genome. This dramatic, sudden process impinges on many genes – generating multiple coding changes that profoundly affect the transcriptome, with the disease-defining gene fusion at its core. These rearrangement loops emerge in an early ES clone from which both the primary tumor and the lethal relapse emerged, and then evolved in parallel until clinically detected.

  Study EGAS00001003385

• Onco-exaptation of an Endogenous Retroviral LTR Drives IRF5 Expression in Hodgkin Lymphoma

  The transcription factor interferon regulatory factor 5 (IRF5) is upregulated in Hodgkin lymphoma (HL) and a key regulator of the aberrant transcriptome characteristic of this disease. Here we show that IRF5 upregulation in HL is driven by transcriptional activation of a normally dormant endogenous retroviral LOR1a long terminal repeat (LTR) upstream of IRF5. Specifically, through screening of RNA-seq libraries, we detected LTR-IRF5 chimeric transcripts in multiple HL cell lines but not in normal B-cell controls. In HL, the LTR was in an open and hypomethylated epigenetic state and we further show the LTR is the site of transcriptional initiation. Among HL cell lines, usage of the LTR promoter strongly correlates with overall levels of IRF5 mRNA and protein, indicating that LTR transcriptional awakening is a major contributor to IRF5 upregulation in HL. Taken together, oncogenic IRF5 overexpression in HL is the result of a specific LTR transcriptional activation. We propose that such LTR derepression is a distinct mechanism of oncogene activation (“onco-exaptation”) and that such a mechanism warrants further investigation in molecular and cancer research.

  Study EGAS00001001205

• Dissection of the molecular complexity of colorectal cancer in pre-clinical models identifies predictive signatures of sensitivity to EGFR inhibitors

  Colorectal carcinoma (CRC) is a heterogeneous tumour entity, with only a fraction of the tumours responding to currently available therapies. A detailed understanding of the biology and genetics of CRC is a prerequisite for precision medicine. To address this challenge in clinical oncology, the OncoTrack consortium recruited 106 CRC patients (stages I-IV) and developed an integrated pre-clinical platform analysing drug sensitivity in the established matched patient-derived in vivo and in vitro models, 46 three-dimensional cell cultures and 59 xenografts, recapitulating many of the molecular features of the donor tumours. We generated extensive omics information on patient tumours and their derived models used to generate a compendium of drug sensitivity data. We could distinguish tumour molecular sub-groups responding differently to therapy and identified a novel independent predictor of sensitivity to the EGFR inhibitor cetuximab. The findings reported here will advance our understanding of CRC and its response to therapy.

  Study EGAS00001001752

• Whole genome and transcriptome analysis of medullary thyroid cancer

  Medullary thyroid cancer is the neoplasm of the calcitonin-producing parafollicular cells of the thyroid gland. The hormone calcitonin acts to lower the blood calcium level, counteracting the effects of the parathyroid hormone; it can also serve as a biomarker for the presence or recurrence of medullary thyroid cancer. Surgery is the only curative treatment for this cancer. Targeted therapeutics vandetanib and cabozantinib are approved for the treatment of aggressive and metastatic tumors that are not amenable to surgery. The use of these multikinase inhibitors are supported by the observed overactivation of the RET oncoprotein in a large sub-population of medullary thyroid cancers; however, not all patients carry oncogenic alterations of this kinase. Hence, there is a need for molecular characterization of medullary tumors on the genomic scale with the aim of identifying targetable mutations. Here, we describe the genomic profiles of two medullary thyroid cancers and report the presence of a putative oncogenic BRAF fusion in one. Such alterations previously observed in other malignancies and known targets of available drugs, can benefit patients who currently have no therapy options.

  Study EGAS00001001473

• Functional single-cell characterization of immune aplastic anemia shows convergence of NK and NK-like CD8+ T cells with disease-associated TCR signature

  Immune aplastic anemia (AA) is a life-threatening bone marrow (BM) failure disorder driven by an autoimmune T cell attack against the hematopoietic stem and progenitor cells (HSPCs). However, the autoantigen targets and the role of other immune cells are unknown. Here, by analyzing an international cohort of 156 AA patients with scRNA+TCRαβ-seq, TCRβ-seq, flow cytometry, and plasma cytokine profiling with comparisons to healthy controls and patients with other hematological disorders, we identify NK cells and CD8+ TEMRA cells expressing NK receptors with AA-associated TCRβ-seq motifs as the most dysregulated immune cell populations in AA BM. Functional co-culture experiments with scRNA-TCRαβ-seq readout utilizing primary HSPCs and immune cells provide evidence that NK cellscannot kill HSPCs alone, but sensitize them to CD8+ T cell-mediated killing. Furthermore, HSPCs induce strong activation of T cell clones with CD8+ TEMRA NK-like phenotype and AA-associated TCR motifs. Our results reveal the convergent evolution of innate and adaptive immune cells in AA, where NK cells support CD8+ T cell-mediated autoimmunity.

  Study EGAS00001007604

• noninvasive lung cancer subtyping

  Accurate diagnosis of lung cancer is important for treatment decision-making. Currently, the gold standard for diagnosing histological subtypes of lung cancer relies on tumor biopsies. Recently, liquid biopsy, particularly cell-free DNA (cfDNA), has shown promising results in cancer detection and classification. In this study, we investigated the potential of cfDNA methylome for the noninvasive classification of lung cancer histological subtypes. Specifically, we focused on the two most prevalent lung cancer subtypes, lung adenocarcinoma and lung squamous cell carcinoma. Using a fragment-based marker discovery approach, we identified robust subtype-specific methylation markers from tumor samples. These markers were successfully validated in independent cohorts and associated with subtype-specific transcription activity. Leveraging these markers, we constructed a subtype classification model using cfDNA methylation profiles, achieving an AUC of 0.808 in the cross-validation and an AUC of 0.747 in the independent validation. Additionally, tumor copy number variations inferred from cfDNA methylome analysis revealed potential implications for treatment selection. In summary, our study demonstrates the potential of cfDNA methylome analysis for noninvasively lung cancer subtyping, offering insights for cancer monitoring and early detection.

  Study EGAS00001007717

• An_exome_sequencing_study_of_the_HIV_elite_long_term_non_progressors_and_rapid_progressors__CASCADE_cohorts_

  Background: A rare subgroup of HIV infected individuals naturally controls infection withouttreatment. These ?elite controllers? constitute an important model for the natural control ofHIV infection. Indeed, the study of these individuals may provide insights into strategies forthe development of HIV vaccines. Although several HLA and chemokine alleles are knownto be over-represented in elite controllers, only a small portion of HIV phenotypic variation isexplained by known genetic variants. The elite controller phenotype is rare and distinct,representing the extreme of an infectious disease trait. As such, this phenotype may be partlyexplained by variation in host immune control, which may be characterized by differences inrare functional genetic variants. Genomic regions underlying elite control can be potentiallyidentified by comparing the presence or frequency of variants in this group to thatrepresenting the opposite extreme. In this context, ?rapid progressors? is a group defined byits rapid immunological and clinical disease progression.Aim: To extend an existing study, in order to identify DNA sequence variants involved in thecontrol of HIV infection with greater statistical resolution. Specifically, we aim to sequence upto 200 exomes from multiple cohort studies within the EuroCoord CASCADE collaboration (acollaboration of 25 HIV seroconversion cohort studies across Europe).

  Study EGAS00001000522

• The Ovarian Cancer Association Consortium OncoArray genome-wide association study

  A custom Illumina genotyping array, the OncoArray. The array includes a genome-wide backbone, comprising 230,000 SNPs tagging most common genetic variants, together with dense mapping of known susceptibility regions, rare variants from sequencing experiments, pharmacogenetic markers, and cancer-related traits. A total of 494,763 SNPs passed quality control steps with a sample success rate of 97% of the samples. Participating sites performed ancestry analysis using a common set of markers and a scoring algorithm based on principal components analysis as described in PMID:27697780.The array was used to genotype ovarian cancer cases and controls from multiple studies participating in the Ovarian Cancer Association Consortium, of which 54,747 passes QC.An initial principal components analysis was carried out using a method described in PMID:27697780 to assign subject to intercontinental ancestry groups: “European”, “African”, “Asian” and "Other". Within ancestry principal components analysis was carried using a set of ~33k unlinked markers. Subsequent analyses were adjusted for the ancestry specific PCs: 9 Eur, 10 Asian, 1 African and 6 Mixed. These PCs are included in the imputed dataset phenotypes file

  Study EGAS00001002305

• Whole Mitochondrial Sequencing of Gingivo-buccal Cancer: ICGC-India Project

  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. Detailed clinical characterization of the patients, collection of data on demographic and environmental exposures, and isolation of DNA samples from blood and tumor tissues collected from each patient are being done at the Advanced Centre for Research, Treatment and Education on Cancer, Mumbai. Whole genome and other sequencing experiments are being performed at the National Institute of Biomedical Genomics, Kalyani. Whole mitochondria sequencing of 89 patients was performed in this study. The mean depth of sequencing was at least 2500X, but variable across platforms, 2689X and 9894X for amplicon sequencing in Ion Torrent PGM (Tumor-blood pairs=16) and Illumina HiSeq 2500 (Tumor-blood pairs =48), respectively; and, 14,190X for whole genome sequencing in Illumina HiSeq 2500 (Tumor-blood pairs =25). Associations of mutation data with important clinical aspects have been found.

  Study EGAS00001002425

• The demographic history and mutational load of African hunter-gatherers and farmers

  The distribution of deleterious genetic variation across human populations is a key issue in evolutionary biology and medical genetics. However, the impact of different modes of subsistence on recent changes in population size, patterns of gene flow, and deleterious mutational load remains unclear. Here, we report high-coverage exomes from various populations of rainforest hunter-gatherers and farmers from central Africa. We find that the recent demographic histories of hunter-gatherers and farmers differed considerably, with population collapses for hunter-gatherers and expansions for farmers, accompanied by increased gene flow. We show that purifying selection against newly arising deleterious alleles is of similar efficiency across African populations, in contrast with Europeans where we detect weaker purifying selection. Furthermore, the per-individual mutation load of rainforest hunter-gatherers is similar to that of farmers, under both additive and recessive models. Our results indicate that differences in the cultural practices and demographic regimes of African populations have not resulted in large differences in mutational burden, and highlight the beneficial role of gene flow in reshaping the distribution of deleterious genetic variation across human populations.

  Study EGAS00001002457

• Evolutionary analysis of primary tumors and metastatic lesions from 20 breast cancer patients (99 samples in total) using exome sequencing data.

  Metastatic breast cancers are still incurable. Characterizing their evolutionary landscape including the role of metastatic axillary lymph nodes to seed distant organ metastasis can provide rational basis for effective treatments. Here, we describe the genomic analyses of the primary tumors and metastatic lesions from 20 breast cancer patients (99 samples). Our evolutionary analyses revealed diverse spreading and seeding patterns governing tumor progression. Although linear evolution to successive metastatic sites was common, parallel evolution from the primary tumor to multiple distant sites was also evident. Metastatic spreading was frequently coupled with polyclonal seeding, where multiple metastatic subclones originated from primary tumor and/or other distant metastases. Synchronous axillary lymph node metastasis, a well-established prognosticator of breast cancer, was not involved in seeding distant metastasis, suggesting haematogenous route for cancer dissemination. Clonal evolution coincided frequently with emerging driver alterations and evolving mutational processes, notably a significant increase in apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) associated mutagenesis. Our data provide the first genomic evidence regarding the role of axillary lymph node metastasis in seeding distant organ metastasis and elucidates the evolving mutational landscape during cancer progression.

  Study EGAS00001002737

• Defective T-cell expansion in RASGRP1 deficiency

  Inherited CTPS1- CD27- and CD70-deficiencies in humans have revealed key factors of T-lymphocyte expansion, that is a prerequisite for an efficient immunity to Epstein-Barr virus (EBV) infection. RASGRP1 is a T-lymphocyte specific nucleotide exchange factor known to activate the MAP kinases pathway. A deleterious homozygous mutation in RASGRP1 leading to the loss RASGRP1 expression was identified in two siblings who both developed a persistent EBV infection leading to Hodgkin lymphoma. RASGRP1-deficient T cells exhibited defective MAPK activation and impaired proliferation that was restored by expression of wild-type RASGRP1. Similar defects were observed in T cells from healthy individuals when RASGRP1 was downregulated. RASGRP1-deficient T cells also exhibited decreased CD27-dependent proliferation toward CD70-expressing EBV-transformed B cells, a crucial pathway required for expansion of antigen-specific T cells in anti-EBV immunity. Furthermore, RASGRP1-deficient T cells failed to upregulate CTPS1, an important enzyme required for DNA synthesis. These results show that RASGRP1 deficiency leads to susceptibility to EBV infection, and demonstrate the key role of RASGRP1 at the crossroad of pathways required for the expansion of activated T lymphocytes.

  Study EGAS00001002753

• Somatic IL4R Hotspot Mutations in Primary Mediastinal Large B-cell lymphoma lead to constitutive JAK-STAT activation

  Primary mediastinal large B cell lymphoma (PMBCL) is a distinct subtype of diffuse large B cell lymphoma thought to arise from thymic medullary B cells. Gene mutations underlying the molecular pathogenesis of the disease are incompletely characterized. Here, we describe novel somatic IL4R mutations in 15 out of 62 primary cases of PMBCL (24.2%) and in all PMBCL-derived cell lines tested. The majority of mutations (11/21; 52%) were hotspot single nucleotide variants in exon 8 leading to an I242N amino acid change in the transmembrane domain. Functional analyses establish this mutation as gain-of-function leading to constitutive activation of the JAK-STAT pathway and upregulation of downstream cytokine expression profiles and B cell specific antigens. Moreover, expression of I242N mutant IL4R in a mouse xenotransplantation model conferred growth advantage in vivo. The pattern of concurrent mutations within the JAK-STAT signaling pathway suggests additive/synergistic effects of these gene mutations contributing to lymphomagenesis. Our data establish IL4R mutations as novel driver alterations and provide a strong preclinical rationale for therapeutic targeting of JAK-STAT signaling in PMBCL.

  Study EGAS00001002796

• Hominin-specific NOTCH2 paralogs expand human cortical neurogenesis through regulation of Delta/Notch interactions.

  The human cerebral cortex has undergone rapid expansion and increased complexity during recent evolution. Hominid-specific gene duplications represent a major driving force of evolution, but their impact on human brain evolution remains unclear. Using tailored RNA sequencing (RNAseq), we profiled the spatial and temporal expression of Hominid-specific duplicated (HS) genes in the human fetal cortex, leading to the identification of a repertoire of 36 HS genes displaying robust and dynamic patterns during cortical neurogenesis. Among these we focused on NOTCH2NL, previously uncharacterized HS paralogs of NOTCH2. NOTCH2NL promote the clonal expansion of human cortical progenitors by increasing self-renewal, ultimately leading to higher neuronal output. NOTCH2NL function by activating the Notch pathway, through inhibition of Delta/Notch interactions. Our study uncovers a large repertoire of recently evolved genes linking genomic evolution to human brain development, and reveals how hominin-specific NOTCH paralogs may have contributed to the expansion of the human cortex. The dataset contains raw sequences (FASTQ files) from the Illumina 2x150bp paired-end RNA sequencing profiles of 11 fetal human brain samples at 7, 9, 12, 15 and 21 gestational weeks.

  Study EGAS00001002798

• PATL2 is a key actor of oocyte maturation whose invalidation causes infertility in women and mice

  The genetic causes of oocyte meiotic deficiency (OMD), a form of primary infertility characterised by the production of immature oocytes, remain largely unexplored. Using whole exome sequencing, we found that 26% of a cohort of 23 subjects with OMD harboured the same homozygous nonsense pathogenic mutation in PATL2, a gene encoding a putative RNA-binding protein. Using Patl2 knockout mice, we confirmed that PATL2 deficiency disturbs oocyte maturation, since oocytes and zygotes exhibit morphological and developmental defects respectively. PATL2's amphibian orthologue is involved in the regulation of oocyte mRNA as a partner of CPEB. However, Patl2's expression profile throughout oocyte development in mice, alongside colocalisation experiments with Cpeb1, Msy2 and Ddx6 (three oocyte RNA-regulators) suggest an original role for Patl2 in Mammals. Accordingly, transcriptomic analysis of oocytes from WT and Patl2-/- animals demonstrated that in the absence of Patl2, expression levels of a select number of highly relevant genes involved in oocyte maturation and early embryonic development are deregulated. In conclusion, PATL2 is a novel actor of mammalian oocyte maturation whose invalidation causes OMD in humans.

  Study EGAS00001002903

• Japan PBC-GWAS Haplotype Study

  Primary biliary cholangitis (PBC) susceptibility loci have largely been discovered through single SNP association testing. In this study, we report genic haplotype patterns associated with PBC risk genome-wide in two Japanese cohorts. Among the 74 genic PBC risk haplotype candidates we detected with a novel methodological approach in a discovery cohort of 1,937 Japanese, nearly two-thirds were replicated (49 haplotypes, Bonferroni-corrected P<6.8x10-4) in an independent Japanese cohort (N=949). Along with corroborating known PBC-associated loci (TNFSF15, HLA-DRA), risk haplotypes may potentially model cis-interactions that regulate gene expression. For example, one replicated haplotype association (9q32-9q33.1, OR=1.7, P=3.0x10-21) consists of intergenic SNPs outside of the human leukocyte antigen (HLA) region that overlap regulatory histone mark peaks in liver and blood cells, and are significantly associated with TNFSF8 expression in whole blood. We also replicated a novel haplotype association involving non-HLA SNPs mapped to UMAD1 (7p21.3; OR=15.2, P=3.9x10-9) that overlap enhancer peaks in liver and memory Th cells. Our analysis demonstrates the utility of haplotype association analyses in discovering and characterizing PBC susceptibility loci.

  Study EGAS00001002915

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

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

  Study EGAS00001003539

• Ewings Sarcoma RNA-Seq

  Sarcomas are cancers of the bone and soft tissue often defined by their gene fusions. However, the timing, context, and processes by which these pathogenic fusions arise are unknown. We explored this in Ewing sarcoma, a cancer driven by EWSR1-ETS fusions, with very few cooperating mutations. Combining whole-genome sequencing with enhanced informatics, we found that the EWSR1-ETS fusion arose from striking rearrangement clusters in 42% of cases (52/124). Notably, these were organized in loops that universally contained the fusion at their center, while also weaving up to 18 genes together with it. We found the same pattern of rearrangements in three additional types of sarcoma. From these data, we define a new signature for sarcoma fusions that precedes other somatic changes, in the earliest replicating DNA of the genome. This dramatic, sudden process impinges on many genes – generating multiple coding changes that profoundly affect the transcriptome, with the disease-defining gene fusion at its core. These rearrangement loops emerge in an early ES clone from which both the primary tumor and the lethal relapse emerged, and then evolved in parallel until clinically detected.

  Study EGAS00001003062

• Low and variable tumor reactivity of the intratumoral TCR repertoire in human cancers

  Infiltration of human cancers by T cells is generally interpreted as a sign of immune recognition, and there is a growing effort to reactivate dysfunctional T cells at such tumor sites. However, these efforts only have value if the intratumoral T cell receptor (TCR) repertoire of such cells is intrinsically tumor-reactive, and this has not been established in an unbiased manner for most human cancers. To address this issue, we analyzed the intrinsic tumor-reactivity of the intratumoral TCR repertoire of CD8+ T cells in ovarian and colorectal cancer – two tumor types for which T cell infiltrates form a positive prognostic marker. Data obtained demonstrate that a capacity to recognize autologous tumor is limited to approximately 10% of intratumoral CD8+ T cells. Furthermore, in two out of four patient samples tested, no tumor-reactive TCRs were identified, despite infiltration of their tumors by T cells. These data indicate that the intrinsic capacity of intratumoral T cells to recognize adjacent tumor tissue can be rare and variable, and suggest that clinical efforts to reactivate intratumoral T cells will benefit from approaches that simultaneously increase the quality of the intratumoral TCR repertoire.

  Study EGAS00001003119

• Genomic Profiling of Thyroid Cancer Reveals a Role for Thyroglobulin in Metastasis -Part II

  Papillary thyroid carcinoma (PTC) has a wide geographic variation in incidence, being highest in Saudi Arabia where it is only second to breast cancer as the most common cancer among females. Genomic profiling of PTC from Saudi Arabia has not been attempted previously. We performed whole-exome sequencing of 101 PTC samples and the corresponding genomic DNA to identify genes with recurrent somatic mutations, followed by sequencing of these genes using a next-generation gene panel approach in further 785 samples. In addition to BRAF, N-RAS and H-RAS, which have previously been shown to be recurrently mutated in PTC, our analysis highlights additional genes, including thyroglobulin (TG), which harbored somatic mutations in 3% of the entire cohort. Surprisingly, although TG mutations were not exclusive to mutations in the RAS-MAP kinase pathway, their presence was associated with a significantly worse clinical outcome, which suggests a pathogenic role beyond driving initial oncogenesis. Analysis of metastatic PTC tissue revealed significant enrichment for TG mutations (p<0.001), including events of apparent clonal expansion. Our results suggest a previously unknown role of TG somatic mutations in the pathogenesis of PTC and its malignant evolution.

  Study EGAS00001003357

• Integrated genomic characterization of adrenocortical carcinoma

  Adrenocortical carcinomas (ACC) are aggressive cancers originating in the cortex of the adrenal glands. Despite the overall poor prognosis, ACC outcome is heterogeneous. CTNNB1 and TP53 mutations are frequent in these tumors, but the complete spectrum of genetic changes remains undefined. Exome sequencing and SNP array analysis of 45 ACC revealed recurrent alterations in known drivers (CTNNB1, TP53, CDKN2A, RB1, MEN1) and genes not previously reported to be altered in ACC (ZNRF3, DAXX, TERT and MED12), which were validated in an independent cohort of 77 ACC. The cell-surface transmembrane E3 ubiquitin ligase ZNRF36 was the gene the most frequently altered (21%), and appears as a potential novel tumor suppressor gene related to the ß-catenin pathway. Our integrated genomic analyses led to the identification of two distinct molecular subgroups with opposite outcome. The C1A group of poor outcome ACC was characterized by numerous mutations and DNA methylation alterations, whereas the C1B group with good prognosis displayed a specific deregulation of two miRNA clusters. Thus, aggressive and indolent ACC correspond to two distinct molecular entities, driven by different oncogenic alterations.Array based dataset for this study has been made public at:http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE49280.

  Study EGAS00001000665

• Integrative_genome_profiling_in_AML

  AML emerges as a consequence of accumulating independent genetic aberrations that direct regulation and/or dysfunction of genes resulting in aberrant activation of signalling pathways, resistance to apoptosis and uncontrolled proliferation. Given the significant heterogeneity of AML genomes, AML patients demonstrate a highly variable response rate and poor median survival in response to current chemotherapy regimens. For the past 4 years we have conducted gene expression profiling on purified bone marrow populations equating to normal haematopoietic stem and progenitor cells from healthy subjects and patients with de novo AML in order to identify AML signatures of aberrantly expressed genes in cancer versus normal. We are now applying a series of bioinformatic methodologies combined with clinical and conventional diagnostic data to establish novel genomics strategies for improved prognostication of AML. Additionally, we use our AML signatures to unravel oncogenic signalling pathway activities in AML patients and test inhibitory drugs for these pathways inn preclinical therapeutic programmes. We consider that superimposing GEP and clinical data for our AML patient cohort with additional data on their mutational status will significantly improve the prognostic power of the study as well as unravel yet unknown mutations associated with aberrant signalling activities of oncogenic pathways.

  Study EGAS00001000858

• HCA_Thymus_Disease

  Single cell RNA sequencing on thymic cells from children with syndromic diseases such as Trisomy 21 and inflammatory/autoimmune diseases such as myasthenia gravis. Our group has previously performed atlasing of the healthy human thymus using single cell transcriptomics and obtained an unprecedented understanding of human thymus throughout development. In this study, we intend to survey diseased human thymus and compare this to healthy to obtain a better understanding of thymic dysfunction in human diseases. The diseases we will survey include congenital conditions such as Down's syndrome, CHARGE syndrome, DiGeorge syndrome as well as acquired diseases such as myasthenia gravis and thymic hyperplasia. We intend to perform genetic and genomic study on the samples, including expression analysis. This will not only help us to understand the cells and molecular pathways affected in dysfunctional thymus; but may also lead to generalisable understandings in human genetic diseases and inflammatory/autoimmune diseases. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001004310

• Biased allelic expression in human primary fibroblast single cells.

  The study of gene expression in mammalian single cells using genomic technologies now provides the possibility to investigate the patterns of allelic gene expression. We have used single-cell RNA sequencing to detect the allele-specific mRNA level in 203 single human primary fibroblast cells over 133,633 unique heterozygous single nucleotide variants (hetSNVs). We have observed that at the snapshot of analyses each cell contains mostly transcripts from one allele from the majority of genes; indeed 76.4% of the hetSNVs display stochastic monoallelic expression in single cells. Remarkably, adjacent hetSNVs exhibit haplotype consistent allelic ratio; in contrast distant sites located in two different genes are independent of the haplotype structure. Moreover, the allele-specific expression in single cells correlated with the abundance of the cellular transcript. We observed that genes expressing both alleles in majority of the single cells at a given time point are rare and enriched in highly expressed genes. The relative abundance of each allele in a cell is controlled by some regulatory mechanisms since we observed related single-cell allelic profiles according to genes. Overall, these results have direct implications in cellular phenotypic variability.

  Study EGAS00001001009

• Clinical significance of novel subtypes of acute lymphoblastic leukemia in the context of minimal residual disease-directed therapy

  We evaluate clinical significance of recently identified subtypes of acute lymphoblastic leukemia (ALL) in 598 children treated with minimal residual disease (MRD)-directed therapy. Among the 16 B-ALL and 8 T-ALL subtypes identified by next generation sequencing, ETV6-RUNX1, high-hyperdiploid and DUX4-rearranged B-ALL had the best five-year event-free survival rates (95% to 98.4%); TCF3-PBX1, PAX5alt, T-cell, ETP, iAMP21, and hypodiploid ALL intermediate rates (80.0% to 88.2%); and BCR-ABL1, BCR-ABL1-like and ETV6-RUNX1-like and KMT2A-rearranged ALL the worst rates (64.1% to 76.2%). All but three of the 142 patients with day-8 blood MRD <0.01% remained in remission. Among new subtypes, intensified therapy based on day-15 MRD≥1% improved outcome of DUX4-rearranged, BCR-ABL1-like, and ZNF384-rearranged ALL, and achievement of day-42 MRD<0.01% did not preclude relapse of PAX5alt, MEF2D-rearranged and ETV6-RUNX1-like ALL. Thus, new subtypes including DUX4-rearranged, PAX5alt, BCR-ABL1-like, ETV6-RUNX1-like, MEF2D-rearranged and ZNF384-rearranged ALL have important prognostic and therapeutic implications.

  Study EGAS00001005084

• An_evaluation_of_different_strategies_for_large_scale_pooled_sequencing_study_design_

  Second generation sequencing technology has enabled the design of large-scale sequencing experiments in targeted disease-related regions in thousands of cases and controls. Cost implications dictate a pooled DNA sequencing design, evaluation of allele frequency differences based on the resequenced samples and large-scale targeted follow-up in further samples. To evaluate the feasibility of this approach and to assess different study design strategies, we have carried out a pilot study which tests the feasibility of PCR-based and pull-down DNA pooling for (a) SNP discovery and (b) allele frequency comparison purposes. All samples are from the HapMap or 1958 BC. The sample composition of the pool is listed below: Pool 5: 31 HapMap Individuals (NA12249, NA12156, NA12004, NA11831, NA12716, NA11832, NA11993, NA12057, NA11995, NA12006, NA12144, NA12802, NA12146, NA12005, NA12003, NA07000, NA12043, NA12044, NA11992, NA11881, NA11994, NA07345, NA12154, NA06994, NA06985, NA12239, NA07022, NA07034, NA12155, NA07056, NA06993); 19 1958 BC Individuals (WTCCC88214, WTCCC88215, WTCCC88216, WTCCC88217, WTCCC88222, WTCCC88233, WTCCC88240, WTCCC88241, WTCCC88242, WTCCC88247, WTCCC88249, WTCCC88262, WTCCC88264, WTCCC88278, WTCCC88294, WTCCC88298, WTCCC88302, WTCCC88305, WTCCC88321) Some pools were carried out in duplicate to assess reproducibility. We focused on 10 (PD), 7(PCR) chromosomal regions, ~2(PD), ~1.6 (PCR) Mb in total.

  Study EGAS00001000134

• Rucaparib in patients presenting a metastatic breast cancer with Homologous Recombination Deficiency, without germline BRCA1/2 mutation, pronostic value of high LOH genomic score and predictive value of HRDetect (microarray)

  Breast cancer may present genomic alterations leading to homologous recombination deficiency. PARP inhibitors have proved their efficacy in patients with HER2-negative metastatic breast cancer (mBC) harboring germline (g) BRCA1/2 mutations. We conducted the phase 2 RUBY trial to assess the efficacy of rucaparib in HER2-negative mBC with high genomic loss of heterozygosity (LOH) score or somatic, without gBRCA1/2 mutation. 220 of 711 patients with mBC screened for LOH presented high LOH score which was associated with a higher likelihood of death (HR = 1.39, 95% CI: 1.11-1.75, p = 0.005). The primary objective was not reached with a clinical benefit rate (objective response or SD>16 weeks) of 13.5%. Two LOH-high patients, without somatic BRCA1/2 mutation, presented a complete and durable response (14 and 32 months). HRDetect tended to be associated with response to rucaparib, whithout reaching statistical significance (median HRDetect responders versus non responders: 0.465 versus, 0.040, p = 0.2135). Our data suggests that a small subset of patients with high LOH score could derive benefit from PARP inhibitors.

  Study EGAS00001004755

• UAMS Smoldering Myeloma Myeloma Sequencing

  Identifying high risk smoldering myeloma patients and progression mechanism is a prerequisite to implement effective inception strategies and curve myeloma related morbidity and mortality. We hypothesize that genomics may help identify determinants of progression that may help predict outcome and offer effective chemo preventive targets. Eighty-two patients underwent a custom targeted panel sequencing with an additional capture for the translocation loci. These results were compared to 223 newly diagnosed patients (EGAS00001003223 and EGAD00001004117) and 17 MGUS and 10 early myeloma patients. DNA was obtained from either CD138+ cells from the bone marrow of multiple myeloma patients (tumor) or from stem cell harvests or peripheral blood cells from the same patient (control). 100 ng of DNA was fragmented, end-repaired, and adapters ligated using the HyperPlus kit (KAPA Biosystems). After PCR amplification the libraries were hybridized with probes against either a targeted panel consisting of 140 genes and chromosomal regions (Nimblegen) using SeqCap reagents (Nimblegen). Hybridized libraries underwent further amplification before being sequenced on a NextSeq500 (Illumina) using 75 bp paired end reads. Overall, these data highlight the importance of dysregulation of the MAPK pathway in the progression to MM.

  Study EGAS00001003629

• BIOKEY: A single-cell catalogue of the dynamic changes underlying Checkpoint Immunotherapy response in Early Breast Cancer

  Checkpoint immunotherapy combined with neoadjuvant chemotherapy improves complete pathologic response in a subset of breast cancer patients. Here, we applied single-cell profiling to tumor biopsies collected before and during anti-PD1 therapy. One-third of tumors exhibited proliferative T-cells expanding along CD8+ or CD4+ lineages, which were either characterized by increased cytotoxicity and exhaustion or improved T-helper function, respectively. Lineage tracing in non-expanding tumors revealed at which point in the lineage T-cells were impaired, while gene expression modeling along these lineages revealed novel genes and underlying transcription factors involved in T-cell expansion. Interestingly, different dendritic and myeloid cell phenotypes could either stimulate or inhibit expanding T-cells, while cell-to-cell communication revealed an integrated immune context highly predictive of T-cell expansion, consisting of immune-stimulatory/-inhibitory interactions between cancer and various immune cell types. Our data yield unprecedented insights into the dynamic changes underlying checkpoint immunotherapy response in breast cancer.

  Study EGAS00001004809

• Androgen deprivation therapy promotes an inflammatory and obesity-like microenvironment in periprostatic fat

  Prostate cancer is a leading cause of cancer-related death and morbidity worldwide. Androgen deprivation therapy (ADT) is the cornerstone of management for advanced disease. The use of androgen deprivation therapies is associated with multiple side effects, including metabolic syndrome and truncal obesity. At the same time, obesity has been associated with both prostate cancer development and disease progression, linked to its effects on chronic inflammation at a tissue level. The connection between androgen deprivation therapy, obesity, inflammation, and prostate cancer progression is well-established in clinical settings; however, an understanding of the changes in adipose tissue at the molecular level induced by castrating therapies is missing. Here we investigated the transcriptional changes in periprostatic fat tissue induced by profound androgen deprivation therapy in a group of patients with high-risk tumours compared to a matching untreated cohort. We find that androgen deprivation therapy is associated with a pro-inflammatory and obesity-like adipose tissue microenvironment. This study suggests that the beneficial effect of androgen deprivation therapy may be partially counteracted by metabolic and inflammatory side effects in the adipose tissue surrounding the prostate.

  Study EGAS00001003286

• Genomic Profiling of Thyroid Cancer Reveals a Role for Thyroglobulin in Metastasis

  Papillary thyroid carcinoma (PTC) has a wide geographic variation in incidence, being highest in Saudi Arabia where it is only second to breast cancer as the most common cancer among females. Genomic profiling of PTC from Saudi Arabia has not been attempted previously. We performed whole-exome sequencing of 101 PTC samples and the corresponding genomic DNA to identify genes with recurrent somatic mutations, followed by sequencing of these genes using a next-generation gene panel approach in further 785 samples. In addition to BRAF, N-RAS and H-RAS, which have previously been shown to be recurrently mutated in PTC, our analysis highlights additional genes, including thyroglobulin (TG), which harbored somatic mutations in 3% of the entire cohort. Surprisingly, although TG mutations were not exclusive to mutations in the RAS-MAP kinase pathway, their presence was associated with a significantly worse clinical outcome, which suggests a pathogenic role beyond driving initial oncogenesis. Analysis of metastatic PTC tissue revealed significant enrichment for TG mutations (p<0.001), including events of apparent clonal expansion. Our results suggest a previously unknown role of TG somatic mutations in the pathogenesis of PTC and its malignant evolution.

  Study EGAS00001001788

• The proliferative history shapes the DNA methylome of B-cell tumor and predicts clinical outcome

  Here, we provide access to CLL and DLBCL DNA methylation and gene expression data in the context of a systematic analysis of the DNA methylation variability in 1,595 samples of normal cell subpopulations and 14 tumor subtypes spanning the entire human B-cell lineage. This approach showed that differential methylation among tumor entities relates to differences in cellular origin and to de novo epigenetic alterations, which allowed us to build an accurate machine learning-based diagnostic algorithm. We identify extensive patient-specific methylation variability in silenced chromatin associated with the proliferative history of normal and neoplastic B cells. Mitotic activity generally leaves both hyper- and hypomethylation imprints, but some B-cell neoplasms preferentially gain or lose DNA methylation. Subsequently, we construct a DNA methylation-based mitotic clock called epiCMIT, whose lapse magnitude represents a strong independent prognostic variable in B-cell tumors and is associated with particular driver genetic alterations. Our findings reveal DNA methylation as a holistic tracer of B-cell tumor developmental history, with implications in the differential diagnosis and prediction of clinical outcome.

  Study EGAS00001004640

• Identification of 19 novel loci reveals gene regulatory mechanisms determining susceptibility to testicular germ cell tumour

  Testicular germ cell tumour (TGCT) is the most common cancer in young men1,2, and is characterised by strong inherited genetic risk factors3. Here we have undertaken large-scale genome wide association study (GWAS) for TGCT, encompassing ~7,500 cases and ~23,000 controls, through the Oncoarray consortium. We identified 19 novel loci, approximately doubling the number of known TGCT risk loci to 44 (P<5x10-8)4-14 and conduct deep, high-throughput functional annotation of all risk loci. We establish a network of physical interactions for all risk SNPs to candidate casual genes in 3D space, using high-throughput chromosome conformation capture techniques (HiC) in TGCT cells. Firstly, functional evidence reveals widespread disruption of developmental transcriptional regulators, consistent with failed primordial germ cell differentiation as an initiating step in TGCT oncogenesis15. We secondly observe multiple risk loci associated with defective microtubule assembly, compatible with the high level of aneuploidy observed in TGCTs, suggesting gross chromosomal instability. Finally KIT-MAPK signalling features as a recurrently dysregulated pathway. In summary our findings substantially increase the number of known TGCT risk alleles, and provides a functional basis for disease susceptibility.

  Study EGAS00001001836

• Identification of rare germline variants in familial multiple myeloma

  The risk of developing multiple myeloma (MM) in first-degree relatives of affected individuals is 2 to 4 fold higher than in the general population, implying the existence of inherited susceptibility. While over 20 common, low-penetrance variants have been associated with MM risk, less is known about rare, high-penetrance germline predisposition. We aimed to identify and characterize high/moderated penetrance germline variants in MM families using whole genome sequencing. For that purpose we identified families with at least two individuals diagnosed with MM or its precursors MGUS and smoldering MM. Also, patients with solitary plasmacytoma and AL amyloidosis were enrolled. We also recruited unaffected family members, and they were screened to exclude undetected MM or its precursor stages. After whole genome sequencing of blood DNA and filtering for rare variants (minor allele frequency <0.1%) and family segregation, variants are prioritized by a range of published in silico tools and literature search. The most promising variants are subjected to experimental validation to evaluate their functional relevance. Our data will provide guidance to worldwide efforts to unravel germline predisposition in MM and validation of the proposed variants and genes.

  Study EGAS00001004734

• Familial psychosis associated with a missense mutation at MACF1 gene combined with the rare duplications DUP3p26.3 and DUP16p23.3, affecting the CNTN6 and CDH13 genes

  Psychosis is a highly heritable and heterogeneous psychiatric condition. Its genetic architecture is thought to be the result of the joint effect of common and rare variants. Families with high prevalence are an interesting approach to shed light on the rare variant's contribution without the need of collecting large cohorts. To unravel the genomic architecture of a family enriched for psychosis, with four affected individuals, we applied a system genomic approach based on karyotyping, genotyping by whole-exome sequencing to search for rare single nucleotide polymorphisms (SNPs) and SNP array to search for copy-number variants (CNVs). We identified a rare non-synonymous variant, g.39914279 C>G, in the MACF1 gene, segregating with psychosis. Rare variants in the MACF1 gene have been previously detected in SCZ patients. Besides, two rare CNVs, DUP3p26.3 and DUP16p23.3, were also identified in the family affecting relevant genes (CNTN6 and CDH13, respectively). We hypothesize that the co-segregation of these duplications with the rare variant g.39914279 C>G of MACF1 gene precipitated with schizophrenia and schizoaffective disorder.

  Study EGAS00001004791

• ARAF Mutations Confer Resistance to RAF Dimer Inhibitor Belvarafenib in NRAS- and BRAF- Mutant Melanoma

  Belvarafenib, a potent and selective RAF dimer (type II) inhibitor, exhibits clinical activity in BRAFV600E- and NRAS-mutant melanoma patients. Here, we report the first-in-human phase I study investigating maximum tolerated dose, assessing safety and preliminary efficacy of belvarafenib in BRAF- and RAS-mutated advanced solid tumors (NCT02405065, NCT03118817). Through generation of belvarafenib resistant NRAS- mutant melanoma cells and analysis of circulating tumor DNA from patients treated with belvarafenib, we identified novel recurrent mutations in ARAF within the kinase domain. ARAF mutants conferred resistance to belvarafenib in both a dimer- and kinase activity- dependent manner. Belvarafenib induced ARAF mutant dimers, and dimers containing mutant ARAF were active in the presence of inhibitor. ARAF mutations may serve as a general resistance mechanism for RAF dimer inhibitors as the mutants exhibit reduced sensitivity to a panel of type II RAF inhibitors. The combination of RAF plus MEK inhibition may be used to delay ARAF-driven resistance and suggests a rational combination for clinical use. Taken together, our findings reveal specific and compensatory functions for the ARAF isoform and implicate ARAF mutations as a driver of resistance to RAF dimer inhibitors.

  Study EGAS00001005086

• Risk and modifying factors in Frontotemporal Dementia

  Understanding the molecular mechanisms underlying frontotemporal dementia (FTD) is essential for the development of successful therapies. Here, we present Phase 1 of a multi-omics, multi-model data resource for FTD research which will allows in-depth molecular research into these mechanisms. We have integrated and analysed data from the frontal lobe of FTD patients with mutations in MAPT, GRN and C9orf72 and detected common and distinct dysregulated cellular pathways. Our results highlight that excitatory neurons are the most vulnerable neuronal cell type and that vascular aberrations are a common hallmark in FTD. Via integration of multi-omics data, we detected several transcription factors and pathways which regulate the strong neuroinflammation observed in FTD-GRN. Finally, using small RNA-seq data and verification experiments in cellular models, we identified several up-regulated miRNAs that inhibit cellular trafficking pathways in FTD and lead to microglial activation. In this work we shed light on novel mechanistic and pathophysiological hallmarks of FTD. In addition, we believe that this comprehensive, multi-omics data resource will further mechanistic FTD research by the community.

  Study EGAS00001004895

• Comparative transcriptome of CD34+ hematopoietic progenitors from myeloproliferative patients and control donors

  BCR-ABL-negative MPNs are related to mutations activating the JAK2/STAT pathway and occurring in hematopoietic stem cells (HSC). The most prevalent mutation is JAK2V617F. Calreticulin (CALR) and activating MPLW515 mutations are also observed in ET and PMF. Additional non-MPN mutations may be found, particularly in PMF, and can change disease phenotype or accelerate progression to myelofibrosis or leukemia. In this study, we wanted to compare the transcriptome of CD34+ progenitors from JAK2V617F postive MPN patients to the CD34+ progenitors from control donors. The patients were selected for their high JAK2V617F variant allele frequency (from 60 to 100% VAF). Peripheral blood from patients was collected in EDTA tubes. Hematopoietic progenitors (CD34+) were isolated from mononuclear cells (MNCs) by immunomagnetic enrichment (Miltenyi, Biotech) and were amplified for 5 days in serum-free medium in the presence of a cocktail of human recombinant cytokines containing EPO (1 U/mL) (Amgen Thousand Oaks, CA), TPO (20 ng/mL) (Kirin, Japan), SCF (25 ng/mL) (Biovitrum AB, Sweden), IL-3 (10 ng/mL), FLT3-L (10 ng/mL), G-CSF (20 ng/mL) and IL-6 (100 U/mL) (MiltenyiBiotec).

  Study EGAS00001005106

• Single cell RNA sequencing of synovial B cells in early Rheumatoid Arthritis

  Rheumatoid Arthritis (RA) is a prevalent autoimmune disease characterized by inflammation of the peripheral joints and occurence of autoantibodies, i.e. Rheumatoid Factor (RF) and anti-citrullinated protein antibodies (ACPA), in circulation. Previous studies have shown that there are memory B cells and autoantibody producing plasma cells present in the joint tissue of patients with long-standing RA. However, it has remained unclear, whether these B cell populations are already present in the joint tissue of patients at the disease onset and whether these are autoreactive. Here, we used a single cell RNA sequencing approach to dissect the B cell repertoire at this early timepoint. We found evidence for B and T cell interaction and presence of memory and plasma cell pools in ACPA- and ACPA+ RA. Our results demonstrated clonal relationships between the memory and plasma cell compartments and autoreactivity within the plasma cell pool. These findings challenge our understanding of the local adaptive immune response in the joints of ACPA- and ACPA+ RA patients with direct implications for B and T cell targeting therapy in both patient subgroups.

  Study EGAS00001005144

• Epigenome wide DNA methylation assay of gingivo-buccal oral squamous cell carcinoma using single base resolution high throughput array

  Gingivo-buccal oral squamous cell carcinoma (OSCC-GB) is the most common cancer among men in India and is associated with high mortality. Although OSCC-GB is known to be quite different from tongue cancer in its genomic presentation and its clinical behavior, it is treated identically as tongue cancer. Predictive markers of prognosis and therapy that are specific to OSCC-GB are, therefore, required. To achieve this, we have carried out epigenomic (n=87) profiling of paired tumor-normal tissues collected from OSCC-GB patients from India and genome-wide DNA methylation assays were performed. DNA was isolated using DNeasy Blood and Tissue Kit (Qiagen). The purity and concentration was estimated using NanoDrop 2000 (Thermo Fisher Scientific). Approximately 500 ng genomic DNA from each sample was used for sodium bisulfite conversion using the EZ DNA methylation Gold Kit (Zymo Research, USA). Genome-wide DNA methylation was assayed using iScan (Illumina), for paired tumour and adjacent normal samples of 25 patients using the Infinium Human Methylation450 BeadChip and of 62 patients using the Infinium MethylationEPIC BeadChip.

  Study EGAS00001003896

• Whole genome sequence and RNA-seq data from paired tumour and germline samples from mesothelioma patients.

  Malignant pleural mesothelioma (MPM) has a poor overall survival with few treatment options. Whole genome sequencing (WGS) combined with RNA sequencing (RNA-seq) analysis of the immune features of MPM offers the prospect of identifying changes that could inform future clinical trials. We analysed somatic mutation and RNA-seq data from 229 MPM samples, including 58 MPM samples that had undergone WGS from our own institutions, together with other published data. This combined analysis identified somatic driver genes, including newly identified candidate genes. Whole genome doubling was a frequent event that correlated with shorter survival. Mutational signature analysis revealed dominant signatures and showed that defects in homologous recombination repair were infrequent in our cohort. Within the tumour immune environment we identified high M2 macrophage infiltrate linked with MMP2, MMP14, TGFB1 and CCL2 expression, representing an immunosuppressive environment. A small subset of samples had a higher proportion of CD8 T cells and a high cytolytic score, suggesting a ‘hot’ immune environment which is independent of the somatic mutations. We propose that our findings on genomic changes and subtypes of immune microenvironments may influence therapeutic planning in the future.

  Study EGAS00001005196