13342 results for ""60+mice"+"6+weeks"+diet+study+2014+2017+DOI"

in 56.94 milliseconds.

• Colorectal cancer organoids expressing BRAF (fusion) genes

  Fusion genes may be oncogenic drivers and potential targets for personalized therapies in a variety of cancer types. The BRAF gene is frequently involved in oncogenic fusions, with fusion frequencies of 0.2-3% throughout different cancers. BRAF can be fused to a wide variety of genes, however, BRAF fusions rarely occur in the exact same gene configuration, making assessment of the fusion relevance and decision for drug treatment challenging. Here, we devised a colorectal cancer organoid-based platform to functionally characterize diverse BRAF fusions, containing various partner genes (AGAP3, DLG1 and TRIM24). We compared these BRAF fusions with respect to cellular behaviour, downstream signaling activation and response to targeted therapies. We found that 5’ partner choice of BRAF affects cellular localization and intracellular signaling capacities of the fusion genes. The DLG1-BRAF fusion gene showed distinct localization to the plasma membrane and exhibited increased levels of MAPK pathway activation under unperturbed conditions. Furthermore, the different BRAF fusions showed varying sensitivities to the targeted inhibition of ERK and BRAF, with the DLG1-BRAF fusion being the most sensitive. Additionally, RNA-sequencing identified distinct subsets of genes affected by the DLG1-BRAF fusion gene. Importantly, all fusion genes conveyed resistance to the clinically relevant EGFR/HER2/HER4-inhibitor afatinib, suggesting that BRAF fusions should be screened alongside other MAPK pathway alterations to identify mCRC patients amenable to cetuximab treatment. In summary, we developed a platform to efficiently assess molecular and cellular effects of fusion genes, and revealed that differential drug responses and distinct gene expression profiles can be induced by different BRAF fusions .

  Study EGAS00001003558

• Phylogenetic analysis of treatment-naive metastases using whole exome and genome sequencing data

  Metastases are responsible for the majority of cancer related deaths and are often difficult to treat successfully. Even for metastases that occur subsequent to treatment, we do not know whether such relapsing metastases were originally heterogeneous or if the observed heterogeneity is a consequence of therapy. Given that treatment can influence evolutionary dynamics by selecting resistant clones, imposing bottlenecks on cancer cell populations, and even inducing novel somatic mutations, characterizing the prior standing variation among metastases remains an important goal to predict initial therapeutic response. To quantify the heterogeneity at clinical presentation of advanced disease, we surveyed the literature for patients in which at least two treatment-naïve metastases underwent genome/exome-wide sequencing. Across all cancer types surveyed, only 18 subjects were found to fulfill this requirement. Including previously unpublished data from two subjects, we analyzed data of 74 untreated metastases and inferred cancer phylogenies. Putative driver gene mutations were acquired at a 2-fold higher rate on the trunk of all metastases than on branches. Mutations in driver genes along the trunk were strongly enriched for predicted functional consequences. Using a stochastic mathematical model, we find that driver gene heterogeneity among metastases mostly occurs in slowly growing cancers for highly advantageous driver gene mutations. Seeding efficiency only weakly affects the heterogeneity among metastases. These results explain why functional driver gene heterogeneity is uncommon in advanced disease prior to treatment, thus providing optimism for future therapies that depend on mutations across metastases.

  Study EGAS00001002777

• Colorectal adenomas, carcinomas and adjacent normal NKI-AvL TGO series NGS-ProToCol

  Cancer is caused by changes in the genome of cells, resulting in aberrant transcripts and proteins. In order to understand why normal cells become cancer cells and how we can identify and destroy them, the CTMM-NGS-ProToCol project aimed to understand prostate and colorectal cancer through DNA and RNA molecular profiling, to further improve diagnosis, prognosis and treatment.The state-of-the-art within the field of genome sequencing makes it possible to identify all the modifications in our genome that encode for RNA and of the transcripts themselves. The latest development in RNA sequencing is that if one does not select for polyA RNA, but first removes rRNA from the total RNA and then performs ‘random’ RNA sequencing, 30-60% of the reads are intronic. Since about 40% of the genome consists of genes (and only 1.6% is exonic), deep RNA sequencing will cover not only the exons, but also a large portion of the introns and therefore much of the genome. The power of the combined sequencing of genome and transcriptome is the full view of copy number aberrations (CNA), exon mutations and the consequence and relevance of the mutations at the RNA level in the form of differential gene expression, alternative splicing, alternative promoter usage, novel transcripts, fusion transcripts, read-through transcripts and mutations. For the colorectal (tumor) tissue dataset, we applied molecular profiling to 30 colorectal adenomas, 30 colorectal carcinomas and 18 normal adjacent colon tissues.

  Study EGAS00001002854

• Integrated clinical, whole genome, and transcriptome analysis of multisampled lethal metastatic prostate cancer

  We report the first combined analysis of whole genome sequence, detailed clinical history, and transcriptome sequence of multiple prostate cancer metastases in a single patient (A21). Whole genome and transcriptome sequence was obtained from 9 anatomically separate metastases, and targeted DNA sequencing was performed in cancerous and noncancerous foci within the primary tumor specimen removed 5 years prior to death. Transcriptome analysis revealed increased expression of AR-regulated genes in liver metastases that harbored an AR p.L702H mutation, suggesting a dominant effect by the mutation despite being present in only 1 of an estimated 16 copies per cell. The metastases harbored several alterations to the PI3K/AKT pathway, including a clonal truncal mutation in PIK3CG and present in all metastatic sites studied. The list of truncal genomic alterations shared by all metastases included homozygous deletion of TP53, hemizygous deletion of RB1 and CHD1, and amplification of FGFR1. If the patient were treated today given this knowledge, use of second-generation androgen-directed therapies, cessation of glucocorticoid administration, and therapeutic inhibition of the PI3K/AKT pathway or FGFR1 receptor could provide personalized benefit. Three previously unreported truncal clonal missense mutations (ABCC4 p.R891L, ALDH9A1 p.W89R, and ASNA1 p.P75R) were expressed at the RNA level and assessed as druggable. The truncal status of mutations is critical for actionability, and can only be determined through analysis of multiple sites of metastasis. Our findings suggest that a large set of deeply analyzed cases could serve as powerful guide to more effective prostate cancer basic science and personalized cancer medicine clinical trials.

  Study EGAS00001001659

• Clonal fitness inferred from timeseries modeling of single cell cancer genomes

  Progress in defining genomic fitness landscapes in cancer by copy number alterations (CNA) has been impeded by lack of single cell and timeseries sampling. We generated 42,000 single cell whole genomes (scWGS) from breast epithelium and primary triple negative breast cancer (TNBC) patient-derived xenografts (PDX) collected during multi-year time series. Using a Wright-Fisher population genetics model, we inferred reproducible CNA-defined clonal fitness dynamics induced by TP53 mutation and cisplatin chemotherapy, with accurate forecasting of experimentally enforced clonal competition dynamics. Drug treatment in three long-term serially passaged TNBC PDX resulted in cisplatin resistant clones that had shown low fitness in the untreated setting. By contrast high fitness clones from treatment naive controls were eradicated reflecting an inversion of the fitness landscape. Upon drug release selective pressure dynamics were reversed indicating a fitness cost of treatment resistance. Taken together, our findings reveal clonal fitness dynamics linked to CNA and therapeutic resistance in polyclonal tumours.

  Study EGAS00001004448

• Versatile workflow for cell-type resolved transcriptional and epigenetic profiles from cryopreserved human lung

  The complexity of the lung microenvironment and changes in cellular composition during disease progression make it exceptionally hard to understand molecular mechanisms leading to the development of chronic lung diseases. Although recent advances in cell-type resolved and single-cell sequencing approaches hold great promise for studying complex diseases, their implementation relies on local access to fresh tissue, as traditional tissue storage methods do not allow viable cell isolation. To overcome these hurdles, we developed a novel, versatile workflow that allows long-term storage of human lung tissue with high cell viability, permits thorough sample quality check before cell isolation, and is compatible with sequencing-based profiling, including single-cell approaches. We demonstrate that cryopreservation is suitable for the isolation of multiple cell types from human lung and is applicable to both healthy and diseased tissue, including COPD and tumour samples. Basal cells isolated from cryopreserved airways retain the ability to differentiate, indicating that cellular identity is not altered by cryopreservation. Importantly, using RNA sequencing and EPIC Array, we show that genome-wide gene expression and DNA methylation signatures are preserved upon cryopreservation, emphasizing the suitability of our workflow for -omics profiling of human lung cells. In addition, we obtained high-quality single-cell RNA sequencing data of cells isolated from the cryopreserved human lung, demonstrating that cryopreservation empowers single-cell approaches. Overall, thanks to its simplicity, our cryopreservation workflow is well-suited for prospective tissue collection by academic collaborators and biobanks, opening worldwide access to human tissue.

  Study EGAS00001004477

• Open_Targets_020_Epigenomes_of_Cell_Lines

  There is currently a drive to establish cell based assay systems of greater human biological and disease relevance through the use of well characterised transformed cell lines, primary cells and complex cellular models (e.g. co-culture, 3D models). However, although the field is gaining valuable experience in running more non-standard & complex cell assays for target validation and compound pharmacology studies, there is the lack of a systematic approach to determine if this expansion in cell assay models is reflected in increased human biological and disease relevance. The increasing wealth of publically available transcriptomic, and epigenome (ENCODE and Epigenome Roadmap) data represents an ideal reference mechanism for determining the relationship between cell types used for target & compound studies to primary human cells and tissues from both healthy volunteers & patients. The CTTV020 epigenomes of cell line project aims to generate epigenetic and transcriptomic profiles of cell lines and compare these with existing and newly generated reference data sets from human tissue and cell types. The aim is to identify assay systems which will provide greater confidence in translating target biology and compound pharmacology to patients. Multiple cell types commonly used within research have been grouped according to biology. Examples include erythroid, lung epithelial, hepatocyte cell types and immortalised models of monocyte / macrophage biology. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001003137

• Open_Targets_020_Epigenomes_of_Cell_Lines

  There is currently a drive to establish cell based assay systems of greater human biological and disease relevance through the use of well characterised transformed cell lines, primary cells and complex cellular models (e.g. co-culture, 3D models). However, although the field is gaining valuable experience in running more non-standard & complex cell assays for target validation and compound pharmacology studies, there is the lack of a systematic approach to determine if this expansion in cell assay models is reflected in increased human biological and disease relevance. The increasing wealth of publically available transcriptomic, and epigenome (ENCODE and Epigenome Roadmap) data represents an ideal reference mechanism for determining the relationship between cell types used for target & compound studies to primary human cells and tissues from both healthy volunteers & patients. The CTTV020 epigenomes of cell line project aims to generate epigenetic and transcriptomic profiles of cell lines and compare these with existing and newly generated reference data sets from human tissue and cell types. The aim is to identify assay systems which will provide greater confidence in translating target biology and compound pharmacology to patients. Multiple cell types commonly used within research have been grouped according to biology. Examples include erythroid, lung epithelial, hepatocyte cell types and immortalised models of monocyte / macrophage biology. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001003139

• Novel optineurin frameshift insertion causing familial frontotemporal dementia and parkinsonism without amyotrophic lateral sclerosis

  Frontotemporal Dementia (FTD) is a common cause of Young Onset Dementia and has diverse clinical manifestations involving behaviour, executive function, language and motor function, including parkinsonism. Up to 50% of FTD patients report a positive family history, supporting a strong genetic basis, particularly in cases with both FTD and amyotrophic lateral sclerosis (FTD-ALS). Mutations in three genes are associated with the majority of familial FTD (fFTD) cases - microtubule associated protein tau gene (MAPT), progranulin gene (PGRN), and hexanucleotide repeat expansions in chromosome 9 open reading frame 72 (C9orf72), whilst mutations in other genes such as optineurin (OPTN) have rarely been reported. Mutations in OPTN have been reported mostly in familial and sporadic cases of ALS, or in rare cases of FTD-ALS, but not in association with pure or predominant FTD and/or parkinsonian phenotype. Here, we report for the first time, a family from the Philippines with 4 affected members harbouring the same OPTN frameshift insertion, presenting with FTD-related phenotypes, including one sibling with predominant parkinsonism resembling corticobasal syndrome. Notably, none of the affected members showed any evidence of motor neuron disease or ALS at the time of writing, both clinically and on electrophysiological testing, expanding the phenotypic spectrum of OPTN mutations. Close follow-up of mutation carriers for the development of new clinical features and wider investigation of additional family members with further genetic analyses will be conducted to investigate the possibility of other genetic modifiers in this family which could explain phenotypic heterogeneity.

  Study EGAS00001005220

• Exome-wide mutation analysis of cell-free DNA to simultaneously monitor the full spectrum of cancer treatment outcomes

  Purpose: Cell-free DNA (cfDNA) offers a non-invasive approach to monitor cancer. Here we develop a method using whole-exome sequencing (WES) of cfDNA for simultaneously monitoring the full spectrum of cancer treatment outcomes, including MRD, recurrence, evolution, and second primary cancers. Experimental Design: Three simulation datasets were generated from 26 cancer patients to benchmark the detection performance of MRD/recurrence and second primary cancers. For further validation, cfDNA samples (n=76) from cancer patients (n=35) with six different cancer types were used for performance validation during various treatments.Results: We present a cfDNA-based cancer monitoring method, named cfTrack. Taking advantage of the broad genome coverage of WES data, cfTrack can sensitively detect MRD and cancer recurrence by integrating signals across known clonal tumor mutations of a patient. In addition, cfTrack detects tumor evolution and second primary cancers by de novo identifying emerging tumor mutations. A series of machine learning and statistical denoising techniques are applied to enhance the detection power. On the simulation data, cfTrack achieved an average AUC of 99% on the validation dataset and 100% on the independent dataset in detecting recurrence in samples with tumor fraction ≥0.05%. In addition, cfTrack yielded an average AUC of 88% in detecting second primary cancers in samples with tumor fraction ≥0.2%. On real data, cfTrack accurately monitors tumor evolution during treatment, which cannot be accomplished by previous methods.Conclusion: Our results demonstrated that cfTrack can sensitively and specifically monitor the full spectrum of cancer treatment outcomes using exome-wide mutation analysis of cfDNA.

  Study EGAS00001005906

• Transcriptome Sequencing PPGL

  RNA library preparation of samples with RNA Integrity Number (RIN) > 5.5 (median 7.1; range 5.6-8.9) was performed as described in the TruSeq Stranded mRNA Library Prep Kit (Illumina, RS-122-2101). Briefly, 0.7-1 µg total RNA PolyA+ fraction was purified and randomly fragmented, converted to double stranded cDNA and processed through subsequent enzymatic treatments of end-repair, dA-tailing, and ligation to adapters. Adapter-ligated library was completed by PCR with Illumina PE primers. The resulting purified cDNA library was applied to an Illumina flow cell for cluster generation and sequenced on an Illumina HiSeq2500 on a 51bp single-read format following manufacturer's protocols. cDNA libraries from FFPE tumors and low integrity RNA’s (RIN<5.5) from FF tumors (median 2.4; range 1-5.5) but % of fragments > 200 nucleotides higher than 20% (median 57.5; range 20-96) were prepared using QuantSeq 3’ mRNA-Seq Library Prep Kit FWD for Illumina (Lexogen, 015) with a UMI Second Strand Synthesis Module for QuantSeq FWD (Lexogen, 081), vendor’s protocol for low input/ low quality/ FFPE RNA using 200-500 ng of total RNA, according to availability, was used as starting material. PCR Add-on Kit for Illumina (Lexogen, 020) was used to adjust the library amplification number of cycles. Libraries were applied to an Illumina flow cell for cluster generation and sequenced on NovaSeq6000. Image analysis, per-cycle basecalling and quality score assignment was performed with Illumina Real Time Analysis software. Demultiplexing of BCL files to FASTQ format was performed with the bcl2fastq Software (Illumina).

  Study EGAS00001006044

• Robust methylation-based classification of brain tumours using nanopore sequencing

  Background: DNA methylation-based classification of cancer provides a comprehensive molecular approach to diagnose tumours. In fact, DNA methylation profiling of human brain tumours already profoundly impacts clinical neuro-oncology. However, current implementations using hybridization microarrays are time-consuming and costly. We recently reported on shallow nanopore whole-genome sequencing for rapid and cost-effective generation of genome-wide 5-methylcytosine profiles as input to supervised classification. Here, we demonstrate that this approach allows to discriminate a wide spectrum of primary brain tumours. Results: Using public reference data of 82 distinct tumour entities, we performed nanopore genome sequencing on N=382 tissue samples covering 46 brain tumour (sub)types. Using bootstrap sampling in a cohort of N = 56 cases, we find that a minimum set of 1,000 random CpG features is sufficient for high-confidence classification by ad hoc random forests. We implemented score recalibration as confidence measure for interpretation in a clinical context and empirically determined a platform-specific threshold in a randomly sampled discovery cohort (N = 185). Applying this cut-off to an independent validation series (N = 184) yielded 148 classifiable cases (sensitivity 80.4%) and demonstrated 100 % specificity. Cross-lab validation demonstrated robustness with concordant results across four laboratories in 10/11 (90.9%) cases. In a prospective benchmarking (N = 15), median time to results was 21.1 hours. Conclusions: In conclusion, nanopore sequencing allows robust and rapid methylation-based classification across the full spectrum of brain tumours. Platform-specific confidence scores facilitate clinical implementation for which prospective evaluation is warranted and ongoing.

  Study EGAS00001006540

• Transcriptional_reprogramming_from_innate_immune_functions_to_a_pro_thrombotic_signature_upon_SARS_CoV_2_sensing_by_monocytes_in_COVID_19

  Alterations in the myeloid immune compartment have been observed in COVID-19, but the specific mechanisms underlying these impairments are not completely understood. Here we examined the functionality of classical CD14+ monocytes as a main myeloid cell component in well-defined cohorts of patients with mild and moderate COVID-19 during the acute phase of infection and compared them to that of healthy individuals. We found that ex vivo isolated CD14+ monocytes from mild and moderate COVID-19 patients display specific expression patterns of costimulatory and inhibitory receptors that clearly distinguish them from healthy monocytes, as well as an altered metabolic profile. In addition, decreased NFB activation in COVID-19 monocytes ex vivo is accompanied by an intact type I IFN antiviral response. Secondary pathogen sensing ex vivo led to a state of functional unresponsiveness characterized by a defect in pro-inflammatory cytokine expression, NFB-driven cytokine responses and defective type I IFN response in moderate COVID-19 monocytes, together with defects in the metabolic reprogramming that innate immune cells usually undergo upon pathogen sensing. Transcriptionally and functionally, COVID-19 monocytes switched their gene expression signature from canonical innate immune functions to a pro-thrombotic phenotype characterized by enrichment of pathways involved in hemostasis, immunothrombosis, platelet aggregation and other accessory pathways to platelet activation, accumulation and clot formation, including extracellular matrix reorganization, integrin cell surface interactions and signaling by PDGF. These results provide a potential mechanism by which innate immune dysfunction in COVID-19 may contribute to disease pathology

  Study EGAS00001006788

• Chronic lymphocytic leukemia IGHV somatic hypermutation detection by targeted capture next-generation sequencing

  Somatic hypermutation (SHM) status of the immunoglobulin heavy variable (IGHV) gene locus plays a crucial role in determining the prognosis and treatment of patients with chronic lymphocytic leukemia (CLL). A common approach for determining SHM status is multiplex PCR of the IGH gene followed by Sanger sequencing; however, this technique is low throughput, vulnerable to failure, and does not allow multiplexing with other molecular diagnostic assays. Here we designed and validated a DNA targeted capture approach that allows the assessment of IGHV SHM status as a sub-module of a larger next-generation sequencing (NGS) assay. We developed an end-to-end assay that takes as input raw FASTQ files and outputs a report containing IGHV SHM status and V-allele usage following the guidelines of the European Research Initiative on CLL (ERIC). We validated the approach on 35 CLL patient samples, 34 of which were previously characterized using Sanger sequencing. The targeted capture NGS panel identified the SHM status of 34 of 35 CLL patients. We showed 100% sensitivity and specificity among the 33 CLL samples with both NGS and Sanger calls. Furthermore, we demonstrated that this panel can be combined with additional targeted capture panels to detect common and prognostically important CLL single nucleotide variants (SNV) and insertion/deletion (indel) variants as well as copy number variants (del17p). Thus, a targeted capture approach to IGHV SHM detection can be easily integrated into broader sequencing panels, allowing broad CLL prognostication in a single molecular assay.

  Study EGAS00001006887

• PD-L1 blockade in combination with carboplatin as immune induction in metastatic lobular breast cancer: the GELATO-trial

  Invasive lobular breast cancer (ILC) is the second most common histological breast cancer subtype but published data on trials specific for ILC are so far lacking. Translational research revealed that a subset of ILCs may be immune-related and more sensitive to DNA-damaging agents such as platinum. In murine ILC models, synergy between immune checkpoint blockade and platinum has been observed. Here, we tested this concept in the phase II, GELATO-trial (NCT03147040), in which patients with metastatic ILC were treated with weekly carboplatin (AUC 1.5) as immune induction treatment for 12 weeks and atezolizumab (PD-L1 blockade; every three weeks) from the third week onwards until disease progression. Four out of 23 evaluable patients had a partial response (17%, 95%CI 5-39%) and two patients had stable disease for at least 24 weeks, resulting in a clinical benefit rate of 26% (95%CI 10-48%). Out of these six patients, four patients had triple-negative ILC (TN-ILC). In serial biopsies of metastatic lesions, we observed higher CD8 T-cell infiltration, expression of immune checkpoints, and exhausted T cells upon carboplatin/PD-L1 blockade. This is the first report of a clinical trial specifically for ILC and we demonstrate promising anti-tumor activity of atezolizumab with carboplatin as immune induction, in particular for TN-ILC. While activity of carboplatin/PD-L1 blockade in classical ER+ ILC was limited, our translational data yield important insights for the design of highly needed clinical trials in ILC.

  Study EGAS00001006902

• A novel transcriptional signature identifies T-cell infiltration in high-risk paediatric cancer

  Background Molecular profiling of the tumour immune microenvironment (TIME) has enabled the rational choice of immunotherapies in some adult cancers. In contrast, the TIME of paediatric cancers is relatively unexplored. We speculated that a more refined appreciation of the TIME in childhood cancers, rather than a reliance on commonly used biomarkers such as tumour mutation burden (TMB), neoantigen load and PD-L1 expression, is an essential prerequisite for improved immunotherapies in childhood solid cancers. Methods We combined immunohistochemistry (IHC) with RNA-sequencing and whole genome sequencing across a diverse spectrum of high-risk paediatric cancers to develop an alternative, expression-based signature associated with CD8+ T-cell infiltration of the TIME. Further, we explored transcriptional features of immune archetypes, T-cell receptor sequencing diversity, assessed the relationship between CD8+ and CD4+ abundance by IHC and deconvolution predictions, and assessed common adult biomarkers such as neoantigen load and TMB. Results A novel 15-gene immune signature, Immune Paediatric Signature Score (IPASS), was identified. Using this signature, we estimate up to 31% of high-risk cancers harbour infiltrating T-cells. In addition, we showed that PD-L1 protein expression is poorly correlated with PD-L1 RNA expression and TMB and neoantigen load are not predictive of T-cell infiltration in paediatrics. Furthermore, deconvolution algorithms are only weakly correlated with IHC measurements of T-cells. Conclusions Our data provides new insights into the variable immune-suppressive mechanisms dampening responses in paediatric solid cancers. Effective immune-based interventions in high-risk paediatric cancer will require individualised analysis of the TIME.

  Study EGAS00001007029

• Treatment stratification and biomarker validation using patient-derived head and neck cancer organoids

  Background Organoids are three-dimensional structures that can be grown from patient-derived material including carcinomas. These in vitro models can be cryopreserved to generate “living” biobanks that reflect patient heterogeneity. Head and neck cancer (HNC) is a collective term used for tumors arising in the head and neck area, including squamous cell carcinomas (HNSCC) and adenocarcinoma of the salivary gland. Methods to culture patient-derived HNC organoids have previously been described. Methods Organoids were derived from resection or biopsy material of patients included between 2016 and 2022. Organoids were characterised with immunohistochemistry and whole exome sequencing (WES) and exposed to chemo-, radio- and targeted-therapies. CRISPR/Cas9-based gene editing was applied to validate proposed biomarkers. Results We have extended our previously published organoid biobank to a total of 110 organoid models, some of which were characterised by histology and WES. Organoids recapitulated both single nucleotide variants and copy number alterations as found in patient HNSCC samples. Organoid and patient response to radiotherapy (RT; primary (n=6) and adjuvant (n=15)) could be compared in 21 HNC cases. In the adjuvant RT group, organoids derived from patients that relapsed were more resistant to RT in vitro, showing potential for guiding treatment options in this subset of patients. Conclusion Taken together, these results indicate illustrate the potential of organoids as a diagnostic tool in personalised medicine for HNC and indicate the potential of this technology for biomarker discovery and validation within the context of relevant patient heterogeneity.

  Study EGAS00001007076

• Whole Exome Sequencing data of six chRCC tumors corresponding to three patients

  Mammalian target of rapamycin (mTOR) is a central regulator of mammalian metabolism and physiology. Aberrant hyperactivation of the mTOR pathway promotes tumor growth and metastasis, and can also promote tumor resistance to chemotherapy and cancer drugs; this makes mTOR an attractive cancer therapeutic target. mTOR inhibitors have been approved to treat cancer; however, the mechanisms underlying drug sensitivity remain poorly understood. Here, whole exome sequencing of three chromophobe renal cell carcinoma (chRCC) patients with exceptional mTOR inhibitor sensitivity revealed that all three patients shared somatic mutations in the deubiquitinase USP9X. The clonal characteristics of the mutations, which were amassed by studying multiple patients’ primary and metastatic samples from various years, together with the low USP9X mutation rate in unselected chRCC series, reinforced a causal link between USP9X and mTOR inhibitor sensitivity. Rapamycin treatment of USP9X-depleted HeLa and renal cancer 786-O cells, along with the pharmacological inhibition of USP9X, confirmed that this protein plays a role in patients’ sensitivity to mTOR inhibition. USP9X was not found to exert a direct effect on mTORC1, but subsequent ubiquitylome analyses identified p62 as a direct USP9X target. Increased p62 ubiquitination and the augmented rapamycin effect upon bortezomib treatment, together with the results of p62 and LC3 immunofluorescence assays, suggested that dysregulated autophagy in USP9X-depleted cells can have a synergistic effect with mTOR inhibitors. In summary, we show that USP9X constitutes a potential novel marker of sensitivity to mTOR inhibitors in chRCC patients, and represents a clinical strategy for increasing the sensitivity to these drugs.

  Study EGAS00001007104

• scRNAseq for patients with immunodeficiency and HCs

  The NFAT family of transcription factors plays central roles in adaptive immunity in murine models, however, their contribution to human immune homeostasis remains poorly defined. In a multigenerational pedigree, we identified three patients carrying germline biallelic missense variants in NFATC1, presenting with recurrent infections, hypogammaglobulinemia and decreased antibody responses. The compound heterozygous NFATC1 variants identified in the patients caused decreased stability and reduced binding of DNA and interacting proteins. We observed defects in early activation and proliferation of T and B cells from these patients, amenable to reconstitution upon genetic rescue. Following stimulation, T-cell activation and proliferation were impaired, reaching that of healthy controls with delay indicative of an adaptive capacity of the cells. Assessment of the metabolic capacity of patient T cells, revealed that NFATc1-dysfunction rendered T cells unable to engage in glycolysis following stimulation, although oxidative metabolic processes were intact. We hypothesized that NFATc1-mutant T cells could compensate for the energy deficit due to defective glycolysis by enhanced lipid metabolism as an adaptation, leading to a delayed, but not lost activation responses. Indeed, we observed increased 13C-labelled palmitate incorporation into citrate indicating higher fatty acid oxidation and we demonstrated that metformin and rosiglitazone improved patient T-cell effector functions. Collectively, enabled by our molecular dissection of NFATC1 mutations and extending the role of NFATc1 in human immunity beyond receptor signaling, and reveal evidence of metabolic plasticity in the context of impaired glycolysis observed in patient T cells to remedy delayed effector responses.

  Study EGAS00001007271

• Efficacy of nivolumab in pediatric cancers with high mutation burden and mismatch-repair deficiency

  Purpose: Checkpoint-inhibitors have limited efficacy for children with unselected solid and brain tumors. We report the first prospective pediatric trial (NCT02992964) using nivolumab exclusively for refractory non-hematological cancers harboring tumor mutation burden (TMB) ≥5 mutations/megabase (mut/Mb) and/or mismatch-repair deficiency (MMRD). Patients and methods: Twenty patients were screened, and ten ultimately included in the response cohort of whom nine had TMB >10mut/Mb (three initially eligible based on MMRD) and one patient had TMB between 5-10 mut/Mb. Results: Delayed immune responses contributed to best overall response of 50%, improving on initial objective responses (20%) and leading to 2-year overall survival (OS) of 50% (95% CI; 27, 93). Four children, including three with refractory malignant gliomas are in complete remission at a median follow-up of 37-months (range: 32.4-60), culminating in 2-year OS of 43% (95% CI; 18.2, 100). Biomarker analyses confirmed benefit in children with germline MMRD, microsatellite instability, higher activated and lower regulatory circulating T-cells. Stochastic mutation accumulation driven by underlying germline MMRD impacted the tumor microenvironment, contributing to delayed responses. No benefit was observed in the single patient with a MMR-proficient tumour and TMB 7.4 mut/Mb. Conclusions: Nivolumab resulted in durable responses and prolonged survival for the first time in a pediatric trial of refractory hypermutated cancers including malignant gliomas. Novel biomarkers identified here need to be translated rapidly to clinical care to identify children who can benefit from checkpoint-inhibitors, including for upfront management of their cancers.

  Study EGAS00001007393

• Cohesin Mutations in AML

  Cohesin shapes the nuclear chromatin architecture, including enhancer-promoter interactions, and its components, of which especially STAG2 and RAD21, are frequently mutated in myeloid malignancies. To elucidate mechanisms of leukemogenesis associated with cohesin mutations in humans, we comprehensively characterized genetic, epigenetic, transcriptional, and chromatin conformational changes in acute myeloid leukemia (AML). To corroborate our findings, we performed complementary siRNA-mediated depletion of STAG2, its paralogue STAG1 or RAD21 in cord blood-derived CD34+ primary human hematopoietic stem and progenitor cells (HSPCs). We show that STAG2 mutations consistently lead to the loss of STAG2 protein and are associated with a specific set of co-occurring mutations, while STAG1 was never mutated in AML. Loss of STAG2 was frequently compensated by STAG1. Still, specific loci displayed altered cohesin occupancy, gene expression and corresponding changes in local chromatin activation as measured by H3K27ac enrichment and chromatin accessibility. High-throughput chromosome conformation capture (in-situ Hi-C) revealed significantly altered chromatin looping in cohesin-mutated AMLs, including weakened enhancer-promoter contacts with reduced, cohesin-dependent promoter activity. In HSPCs, we detected transcriptomic and epigenetic effects overlapping STAG2-mutant AML-specific changes following STAG2 knockdown (KD), that were not invoked by the depletion of STAG1. We also found that STAG2 loss in cultured HSPCs impaired the differentiation capacity, especially erythroid colony formation which maintained HSPC-like gene expression. This work establishes STAG2 as a key regulator of cohesin-associated chromatin architecture, gene expression and differentiation in the human hematopoietic system and identifies candidate target genes that may be implicated in leukemogenesis.

  Study EGAS00001007405

• Peripheral blood DNA methylome in adalimumab-treated patients with rheumatoid arthritis

  Background and aims Rheumatoid arthritis (RA) patients are currently treated with biological agents mostly aimed at cytokine blockade such as anti-TNFα therapy. Currently, there are no biomarkers to predict therapy response to these agents. Here, we aimed to predict response to adalimumab (ADA) treatment in RA patients using DNA methylation in peripheral blood Methods DNA methylation profiling on whole peripheral blood from 92 RA patients before the start of adalimumab (ADA) was determined using Illumina HumanMethylationEPIC BeadChip Array. After 6 months, treatment response was assessed according to the Alliance of Associations for Rheumatology (EULAR) criteria for disease response. Patients were classified as responders (DAS28<3.2 or decrease of 1.2 points) or as non-responders (DAS28>5.1 or decrease of less than 0.6 points). Machine learning models were built with gradient boosting and stability selection models to predict response prior to ADA treatment with predictor DNA methylation markers. Results We demonstrated a 27-feature panel of DNA methylation markers or CpG classifiers that predict response in RA patients treated with ADA. Forty-nine patients were assigned as responder and 43 patients assigned as non-responders. We differentiate responders from non-responders with a high sensitivity (AUC 0.76). The predictor CpGs annotated to genes involved in immunological- and pathophysiological pathways related to RA such as T-cell signaling, B-cell pathology and angiogenesis. Conclusion Our findings indicate that DNA methylation signatures discriminate responders and non-responders to ADA treatment and can serve as a tool for therapy prediction.

  Study EGAS00001007578

• Evaluation of EBUS-TBNA Aspirates from Advanced NSCLC for Comprehensive Sequencing Platforms Including Whole Genome Sequencing

  Endobronchial ultrasound-guided transbronchial needle aspiration (EBUS-TBNA) is often the only source of tumor tissue from patients with advanced, inoperable lung cancer. EBUS-TBNA aspirates are used for the diagnosis, staging, and genomic testing to inform therapy options. Here we extracted DNA and RNA from 220 EBUS-TBNA aspirates to evaluate their suitability for whole genome (WGS), whole exome (WES), and comprehensive panel sequencing. For a subset of 40 cases, the same nucleic acid extraction was sequenced using WGS, WES, and the TruSight Oncology 500 assay. Genomic features were compared between sequencing platforms and compared with those reported by clinical testing. A total of 204 aspirates (92.7%) had sufficient DNA (100 ng) for comprehensive panel sequencing, and 109 aspirates (49.5%) had sufficient material for WGS. Comprehensive sequencing platforms detected all seven clinically reported tier 1 actionable mutations, an additional three (7%) tier 1 mutations, six (15%) tier 2–3 mutations, and biomarkers of potential immunotherapy benefit (tumor mutation burden and microsatellite instability). As expected, WGS was more suited for the detection and discovery of emerging novel biomarkers of treatment response. WGS could be performed in half of all EBUS-TBNA aspirates, which points to the enormous potential of EBUS-TBNA as source material for large, well-curated discovery-based studies for novel and more effective predictors of treatment response. Comprehensive panel sequencing is possible in the vast majority of fresh EBUS-TBNA aspirates and enhances the detection of actionable mutations over current clinical testing.

  Study EGAS00001007708

• RNA isoform repertoire of neuropsychiatric risk genes in human brain

  Neuropsychiatric disorders are highly complex conditions and the risk of developing a disorder has been tied to hundreds of genomic variants that alter the expression and/or products (isoforms) made by risk genes. However, how these genes contribute to disease risk and onset through altered expression and RNA splicing is not well understood. Here we show our current understanding of gene isoforms is far from complete and reveal the precise splicing profiles of neuropsychiatric disorder risk genes. Combining our new bioinformatic pipeline IsoLamp with nanopore long-read amplicon sequencing, we deeply profiled the RNA isoform repertoire of 31 high-confidence neuropsychiatric disorder risk genes in human brain. We show most risk genes are more complex than previously reported, identifying 443 novel isoforms and 28 novel exons, including isoforms which alter protein domains, and genes such as ATG13 and GATAD2A where most expression was from previously undiscovered isoforms. The greatest isoform diversity was present in the schizophrenia risk gene ITIH4. Mass spectrometry of brain protein isolates confirmed translation of a novel exon skipping event in ITIH4, suggesting a new regulatory mechanism for this gene in brain. Our results emphasize the widespread presence of previously undetected RNA and protein isoforms in brain and provide an effective approach to address this knowledge gap. Uncovering the isoform repertoire of neuropsychiatric risk genes will underpin future analyses of the functional impact these isoforms have on neuropsychiatric disorders, enabling the translation of genomic findings into a pathophysiological understanding of disease.

  Study EGAS00001007744

• OMKar

  The whole genome karyotype refers to the sequence of large chromosomal segments that make up an individuals genotype. karyotype analysis, which includes descriptions of aneuploidies and other rearrangements is crucial for understanding genetic risk factors, for diagnosis, treatment decisions, and genetic counseling linked to constitutional disorders. The current karyotyping standard is based on microscopic examination of chromosomes, a complex process that requires high expertise and offers Mb scale resolution. Optical Genome Mapping (OGM) technology can identify large DNA lesions in a cost-effective manner. In this paper, we developed OMKar, a method that uses OGM data to create a virtual karyotype. OMKar processes Structural (SV) and Copy Number (CN) Variants as inputs and encodes them into a compact breakpoint graph. It recomputes copy numbers using Integer Linear Programming to maintain CN balance and then identifies constrained Eulerian paths representing entire donor chromosomes. In tests using 38 whole genome simulations of constitutional disorders, OMKar reconstructed the karyotype with 88% precision and 95% recall on SV concordance and 95% Jaccard score on CN concordance. We applied OMKar to 50 prenatal, 41 postnatal, and 63 parental samples from ten different sites. OMKar reconstructed the correct karyotype in 144 out of 154 samples, covering 25 of 25 aneuploidies, 32 of 32 balanced translocations, and 72 of 82 unbalanced variations. Detected constitutional disorders included Cri-du-chat, Wolf-Hirschhorn, Prader-Willi deletions, Down, and Turner syndromes. Importantly, it identified a plausible genetic mechanism for five cases of constitutional disorder that were not detected by other technologies. Together, these results demonstrate the robustness of OMKar for OGM-based karyotyping.

  Study EGAS00001008245

• RNAseq data for Failure of Differentiation of the Rhombic Lip Constitutes Medulloblastoma

  We used single-cell transcriptomics to study cells from the developing human cerebellum, and show that different molecular subgroups of medulloblastoma resemble distinct glutamatergic progenitors.

  Dataset EGAD00001008458

• Combination of ribociclib and gemcitabine for the treatment of medulloblastoma

  Group 3 (G3) medulloblastoma (MB) is one of the deadliest forms of the disease for which novel treatment is desperately needed. Here we evaluate ribociclib, a highly selective CDK4/6 inhibitor, with gemcitabine in mouse and human G3 MBs. Ribociclib CNS penetration was assessed by in vivo microdialysis and by immunohistochemistry and gene expression studies. Survival studies to determine the efficacy of ribociclib and gemcitabine combination were performed on mice orthotopically implanted with luciferase labelled mouse and human G3 MB. Pharmacokinetic-pharmacodynamic outcomes and univariable survival models were analyzed to estimate survival. Gene activity inference using NetBID and tumor differentiation analysis investigated the effects of the combination after short and long-term treatments. Tumors from mice treated with oral ribociclib displayed inhibited RB phosphorylation, downregulated E2F target genes, and decreased proliferation. Treatment of mice with the combination of ribociclib and gemcitabine was well tolerated, slowed tumor progression and metastatic spread, and increased survival. Molecular analysis of treated versus untreated tumors showed a significant decrease in the activity and expression of genes involved in cell cycle progression and DNA damage response, and an increase in activity and expression of genes implicated in neuronal identity and neuronal differentiation. Ribociclib is CNS-penetrant. When administered/combined with gemcitabine in orthotopic G3 MB models resulted in improved survival. Our findings, with both mouse and human patient-derived-orthotopic xenograft models, suggest that this combination therapy has promise for children with G3 MB and may represent an effective treatment strategy for other CNS malignancies.

  Study EGAS00001006001

• Surgical Treatment for Ischemic Heart Failure (STICH-BioLINCC)

  Data Access NOTE: Please refer to the “Authorized Access” section below for information about how access to the data from this accession differs from many other dbGaP accessions. ObjectiveTo compare medical therapy with coronary bypass surgery and/or surgical ventricular reconstruction for patients with congestive heart failure and coronary artery disease.BackgroundCoronary artery disease (CAD) is the most common cause of heart failure, which in turn is a major cause of death and disability globally. Evidence from previous clinical trials supports the use of coronary-artery bypass grafting (CABG) to relieve disabling symptoms of angina, particularly among high-risk subgroups with extensive CAD. However, the studies did not include patients with severe left ventricular dysfunction, and developments in medical therapy have since led to updated guidelines. In addition, the benefits of CABG in patients with ischemic cardiomyopathy had still not been clearly established at the time of this study. STICH sought to evaluate the role of CABG in the treatment of patients with CAD and left ventricular systolic dysfunction.Reduced left ventricular function may occur after myocardial infarction, often in conjunction with left ventricular remodeling, including left ventricular enlargement and changes in chamber geometry. Left ventricular remodeling is correlated with progression of heart failure and a poor prognosis. Therefore, a surgical approach to remodeling through left ventricular volume reduction could improve outcomes for patients with CAD and heart failure. Surgical ventricular reconstruction (SVR) has been shown to reduce the left ventricular volume, increase the ejection fraction, and improve ventricular function. There is also evidence that SVR performed with CABG may reduce the rate of hospitalization and improve ventricular function, as compared to CABG alone. As part of a second hypothesis, STICH additionally investigated whether SVR when added to CABG would improve outcomes in patients with heart failure and CAD.ParticipantsA total of 2,136 participants were enrolled in STICH. 1,212 were enrolled in the hypothesis 1 component of the trial, with 602 participants assigned to receive medical therapy alone, and 610 participants assigned to receive medical therapy plus CABG. 1,000 were enrolled in the hypothesis 2 component, with 499 participants assigned to receive medical therapy plus CABG, and 501 participants assigned to receive medical therapy plus CABG and SVR. 76 participants that were assigned to the CABG with medical therapy treatment were enrolled in both hypothesis components.DesignAfter initial determination of overall eligibility, participants were evaluated to determine which component of the STICH program was appropriate for them on the basis of suitable therapeutic options. All participants underwent cardiac imaging for assessment of left ventricular function and wall motion. Participants in the hypothesis 1 component were randomly assigned to receive either medical therapy alone or medical therapy plus CABG. Participants in the hypothesis 2 component were randomly assigned to receive either medical therapy plus CABG or medical therapy plus CABG and SVR.At baseline, demographic factors and clinical characteristics were assessed, including current medications and prior diagnostic and other cardiovascular procedures, and a physical examination was performed. Guideline-based recommendations for drug and device use were emphasized for all participants. All participants underwent follow-up evaluations at the time of discharge or at 30 days for participants still hospitalized, every 4 months for the first year, and every 6 months thereafter. The median length of follow-up was 56 months for hypothesis 1 and 48 months for hypothesis 2.For participants receiving CABG, arterial grafting for stenosis of the left anterior descending coronary artery was required for all participants without specific contraindications. The use of additional arterial conduits supplemented by vein grafts was recommended for revascularization of all major vessels with clinically significant stenoses. Concurrent mitral-valve surgery for regurgitation was performed at the discretion of the surgeon. For participants receiving SVR, the operation was most commonly performed during a single period of cardioplegic arrest after construction of bypass grafts. However, the procedure could also be performed with the heart beating in order to facilitate identification of the noncontractile zone of scarring. After an anterior left ventriculotomy was centered in the zone of anterior asynergy, a suture was placed in the interior of the ventricle to encircle the scar at the boundary between the akinetic and viable tissue. Visual inspection and palpation facilitated the judgment of whether a patch was needed to optimize the chamber size without deforming the left ventricle during closure of the ventriculotomy.The primary outcome for hypothesis 1 was the rate of death from any cause. The primary outcome for hypothesis two was a composite of death from any cause and hospitalization for cardiac causes.ConclusionsThere was no significant difference between medical therapy alone and medical therapy plus CABG for death from any cause, though participants that underwent CABG had lower rates of death from cardiovascular causes, death from any cause, or hospitalization for cardiovascular causes. Adding SVR to CABG reduced the left ventricular volume, as compared with CABG alone. However, this anatomical change was not associated with a greater improvement in symptoms or exercise tolerance, or with a reduction in the rate of death or hospitalization for cardiac causes.

  Study phs003493

• Whole_genome_sequencing_of_100_DDD_trios_with_suspected_noncoding_causal_mutations

  We would like to whole-genome sequence 100 trios (300 individuals) from the DDD project in which we suspect the proband has a pathogenic noncoding mutation, to 30X coverage. These have been carefully selected to maximise the chances of finding pathogenic noncoding mutations. Specifically, we have selected probands who have no likely causal mutation within the coding region but who have a heterozygous loss-of-function mutation inherited from one parent in a known recessive developmental disorder gene, and whose phenotype is a good match for the typical phenotype caused by recessive variants in that gene. We suspect that these probands may have a heterozygous noncoding mutation on the other allele which reduces expression of the gene so that, together, these two mutations are effectively equivalent to a homozygous loss-of-function mutation. About a quarter of the probands have affected siblings and, if those siblings have been exome sequenced, we have checked that they are identical-by-descent in the region of the gene. We propose to search in the promoter region, intronic and annotated regulatory regions around these genes of interest for rare variants inherited from the other parent (i.e. not the one who transmitted the loss-of-function variant), or de novo variants. If we find such variants, and demonstrate an enrichment of this compound heterozygous LoF noncoding mechanism above that expected by chance, we intend to follow up specific candidates by RNAseq and potentially CRISPR/Cas9 experiments. 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 please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001002452

• NHLBI GO-ESP Family Studies: Pulmonary Arterial Hypertension

  The NHLBI "Grand Opportunity" Exome Sequencing Project (GO-ESP), a signature project of the NHLBI Recovery Act investment, was designed to identify genetic variants in coding regions (exons) of the human genome (the "exome") that are associated with heart, lung and blood diseases. These and related diseases that are of high impact to public health and individuals from diverse racial and ethnic groups will be studied. These data may help researchers understand the causes of disease, contributing to better ways to prevent, diagnose, and treat diseases, as well as determine whether to tailor prevention and treatments to specific populations. This could lead to more effective treatments and reduce the likelihood of side effects. GO-ESP is comprised of five collaborative components: 3 cohort consortia - HeartGO, LungGO, and WHISP - and 2 sequencing centers - BroadGO and SeattleGO. Pulmonary arterial hypertension (PAH) is a progressive disease characterized by widespread occlusion of the smallest arteries of the lungs. Pulmonary vascular obstruction leads to increased pulmonary vascular resistance, which subsequently causes heart failure with mean survival of 3 years. PAH occurs at all ages and affects women more than twice as frequently as men. Sporadic Idiopathic pulmonary arterial hypertension (IPAH), comprises 94% of what was formerly known as primary pulmonary hypertension, and is clinically and pathologically indistinguishable from familial PAH (FPAH). Most FPAH is due to mutation in BMPR2, including more than 120 families in the US. Our goal here is to find other genes that are a basis for FPAH, so we selected for exome sequencing 5 families among 40 who do not have mutation in BMPR2, or other known genes (ACVRL1, SMAD8, ENG) that rarely are the basis for FPAH.

  Study phs000354

• Whole Transcriptome Sequencing of Resectable Stage III/IV Melanoma Evaluated After Starting Hu14.18-IL2 Predicts Outcome

  Background: We analyzed whole transcriptome sequencing in tumors from 23 patients with stage III or IV melanoma from a pilot trial of the anti-GD2 immunocytokine, hu14.18-IL2, to identify predictive immune and/or tumor biomarkers in melanoma patients at high risk for recurrence. Methods: Patients were randomized to receive the first of 3 monthly-courses of hu14.18-IL2 immunotherapy either before (Group A) or after (Group B) complete surgical resection of all known disease. Tumors were evaluated by histology and whole transcriptome sequencing. Results: We report here that tumor infiltrating lymphocyte (TIL) levels directly associate with relapse-free survival (RFS) and overall survival (OS) in resected tumors from Group A, where early responses to the immunotherapy agent could be assessed. TIL levels directly associated with a previously reported immune signature which associated with RFS and OS, particularly in Group A tumors. In Group A tumors there were decreased cell cycling gene RNA transcripts, but increased RNA transcripts for repair and growth genes. We found that outcome (RFS and OS) was directly associated with several immune signatures and immune-related RNA transcripts and inversely associated with several tumor growth-associated transcripts, particularly in Group A tumors. Most of these associations were not seen in Group B tumors. Conclusion: We interpret these data to signify that both immunologic and tumoral cell processes, as measured by RNAseq analyses detected shortly after initiation of hu14.18-IL2 therapy are associated with long term survival and could potentially be used as prognostic biomarkers in tumor resection specimens obtained after initiating neoadjuvant immunotherapy. Note: RNA sequences and expression data can be found at NCBI GEO GSE133713

  Study phs001947

• Genomic Analysis of Mycosis Fungoides and Sézary Syndrome Identifies Recurrent Alterations in TNFR2

  Mycosis Fungoides (MF) and Sézary Syndrome (Sz) comprise the majority of Cutaneous T-Cell Lymphoma (CTCL) cases and are characterized by clinical heterogeneity. This array of symptoms includes skin patches, plaques and tumors as well as blood involvement and erythroderma. Because the genetic basis of CTCL is still poorly understood, we performed whole-exome sequencing on 11 MF/Sz samples and their matched, normal control DNA. Upon analyzing this data, we distilled a list of an additional 494 genes to be sequenced at depth. The majority of these 494 genes were sequenced in the exons alone, however, for a small subset of these genes we sequenced the entire genomic locus in search of structural variation events. This 494 gene targeted resequencing effort was performed upon 72 patients with MF/Sz and their matched, normal control DNA. An additional 5 MF/Sz and matched, normal DNA samples underwent whole genome amplification (WGA) prior to library preparation, which appears to have introduced mutations into these samples, however we also include them here. We additionally include 10 MF/Sz samples without a matched, normal control samples, 1 of which underwent WGA. In addition to primary disease samples, we also performed targeted resequencing on 4 CTCL cell lines (HH, MyLa, SeAx and Hut78) as well as on H9 cells (a subclone of Hut78 cells) and Jurkat cells, a T-ALL derived cell line that was used as an experimental system for our work. In sum, this data deposition includes exome and/or targeted resequencing from 92 individuals with 192 samples submitted.

  Study phs000913

• Genetics of Neuropsychiatric and Neurodevelopmental Disorders

  Although schizophrenia is a highly heritable disease, relatively little progress had been made in securely identifying the genetic causes of this disorder, and most instances of schizophrenia in the general population remain unexplained. One avenue of explanation for the genetic basis of schizophrenia, however, has been effectively closed by recent research. Genome-wide association studies (GWAS) have now shown that common variation makes at most a modest contribution to the risk of schizophrenia. At the same time that the role of common variation has been circumscribed by GWAS, however, the analysis of copy number variants that are detectable on a genome-wide scale has revealed and replicated a number of very rare variants that associate with schizophrenia. These rare copy number variants that have been implicated in schizophrenia, however, have one striking feature in common: they are all risk factors for other brain related disorders beyond schizophrenia such as mental retardation, autism and epilepsy. These findings argue that genetic risk factors may confer a highly penetrant vulnerability to neuropsychiatric disorder, which is then further modified by interacting genetic or environmental factors to determine the ultimate manifestation. Most schizophrenia collections that are being studied today, however, have been selected precisely for their homogeneity: including only schizophrenia patients with no comorbidities, or schizophrenia patients with relatives who have schizophrenia but no other neuropsychiatric conditions. These selection criteria are inconsistent with what we now know about the bulk of the genetic differences that have been associated with disease. The central hypothesis of this project is that there are rare genetic variants that strongly elevate the risk of various neuropsychiatric diseases, and that these risk factors can be identified most readily in families segregating multiple neuropsychiatric conditions.

  Study phs000682

• CSER: South-Seq: DNA Sequencing for Newborn Nurseries in the South

  SouthSeq, part of the Clinical Sequencing Evidence-Generating Research (CSER2) effort, aims to perform genome sequencing for infants in nurseries at hospitals in Alabama, Mississippi, Louisiana, and Kentucky. SouthSeq aims to sequence 600 newborn babies presenting with congenital anomalies, or other signs of a rare genetic disorder. The primary goals of SouthSeq are to provide early diagnoses to children with rare genetic conditions, to demonstrate the utility of genome sequencing as a frontline test with potential to improve intermediate and long-term clinical outcomes, and to equip non-genetics providers to deliver test results in order to facilitate the conduction of genome sequencing outside of major medical centers where patients may not have access to a genetics provider. SouthSeq also aims to ensure that genomic testing is applied across diverse communities; recruitment efforts focus on patients from rural and minority populations that have been historically under-represented in both medicine and genomics research. Sequencing and analysis for SouthSeq are being conducted at the HudsonAlpha Institute for Biotechnology in Huntsville, Alabama (https://www.hudsonalpha.org/). Enrollment of infants with signs of genetic disease is approved at five clinical sites, including nurseries at the University of Alabama at Birmingham, the University of Mississippi Medical Center, Woman's Hospital in Baton Rouge, Children's Hospital in New Orleans, and the University of Louisville. Both biological parents are enrolled when available so that follow-up testing can be performed to determine inheritance. Further, participant families may choose to opt-into return of secondary findings identified in the proband, which focus on pathogenic and likely pathogenic variation identified within ACMG SFv3.0 genes.

  Study phs002307

• Genomic Evolution of Low- and High-Grade Glioma

  Malignant transformation (MT) of IDH-mutant low-grade glioma (LGG) to aggressive high-grade tumors is an event of major clinical significance, eventually leading to death in the majority of LGG patients. While patients with primary glioblastoma (GBM) have benefited from increased overall survival due to treatment with the alkylating chemotherapy temozolomide (TMZ), the benefit of TMZ for patients with diffuse low-grade gliomas (LGG) remains unknown. To further explore the relationships among TMZ treatment, hypermutation and malignant progression, we studied tumor evolution in an expanded cohort of untreated and TMZ-treated patients by exome-sequencing of paired IDH1/2-mutant LGGs and their post-TMZ recurrences. Together, these findings further our understanding of the molecular features and clinical behavior of hypermutated clones. The heterogeneity present in individual tumors, i.e. intratumoral heterogeneity (ITH), is thought to be a major reason why targeted treatments eventually fail. However, we know little about ITH in most human tumors because studies typically analyze one biopsy per tumor without knowledge of where in the tumor the sample was taken. Understanding how malignant and non-malignant cellular populations are distributed in 3D tumor space is foundational knowledge for resolving how human tumors grow, with translational implications for identifying representative biopsy specimens and predicting the impact of targeted therapy. Diffuse gliomas are incurable brain tumors with variably aggressive molecular subtypes. To determine the extent of ITH and its 3D organization in diffuse glioma, we prospectively collected a median of 10 spatially mapped samples from newly diagnosed and recurrent patient tumors, including aggressive canonical IDH1/IDH2 wild-type (IDH-wtc) glioblastoma, slower-growing IDH1/IDH2 mutant (IDHmut) glioma and a noncanonical IDH-wt glioblastoma (IDH-wtnc) lacking the TERT promoter mutation typical of canonical Glioblastoma.

  Study phs002034

• Germ Cell Tumor Molecular Heterogeneity Revealed through Analysis of Primary and Metastasis Pairs

  Purpose: While primary germ cell tumors (GCTs) have been extensively characterized, molecular analysis of metastatic sites has been limited. We performed whole-exome sequencing and targeted next-generation sequencing on paired primary and metastatic GCT samples in a patient cohort enriched for cisplatin-resistant disease. Experimental Design: Tissue sequencing was performed on 100 tumor specimens from 50 patients with metastatic GCT, and sequencing of plasma cell-free DNA (cfDNA) was performed for a subset of patients. Data is being provided for a subset of 10 patients whose samples were analyzed using whole-exome sequencing. Results: The mutational landscape of primary and metastatic pairs from GCT patients was highly discordant (68% of all somatic mutations were discordant). Whereas genome duplication was common and highly concordant between primary and metastatic samples, only 25% of primary-metastasis pairs had =50% concordance at the level of DNA copy number alterations (CNAs). Evolutionary-based analyses revealed that most mutations arose after CNAs at the respective loci in both primary and metastatic samples, with oncogenic mutations enriched in the set of early occurring mutations versus variants of unknown significance (VUS). TP53 pathway alterations were identified in nine cisplatin-resistant cases and had the highest degree of concordance in primary and metastatic specimens, consistent with their association with this treatment-resistant phenotype. Conclusions: Analysis of paired primary and metastatic GCT specimens revealed significant molecular heterogeneity for both CNAs and somatic mutations. Among loci demonstrating serial genetic evolution, most somatic mutations arose after CNAs but oncogenic mutations were enriched in the set of early occurring mutations as compared to VUS. Alterations in TP53 were clonal when present and shared among primary-metastasis pairs. Reprinted from pending publication, with permission from JCO Precision Oncology.

  Study phs002229

• Genetic Progression of Head and Neck Squamous Cell Carcinoma

  These data are tumor/normal-matched exome sequencing results on 45 participants (median age of 58.6 years) treated for squamous cell carcinoma of the oropharynx at Mass Eye and Ear (MEE; Boston, MA, USA). Some tumors had two portions sequenced separately. All research complied with all relevant ethical regulations. Human subject research on patients at MEE was approved by the MEE Human Studies Committee under protocol HSC 11-024H, with informed consent obtained from participants. Participants were not compensated. Clinical data, including HPV status, were abstracted from patient records. Sequencing and clinical data were combined with TCGA data on 486 cases of head and neck squamous cell carcinoma (HNSCC). The PhylogicNDT suite of genomic analysis tools was used to infer the relative timing of genetic events during tumor progression for major molecularly defined classes of HNSCC: HNSCC negative for human papillomavirus (HPV), HPV-associated HNSCC, and HPV-negative HNSCC involving mutations in CASP8, NSD1, and NOTCH1. The timing of the development of whole-genome somatic events leading to tetraploidy or triploidy and to increased measures of intratumor genetic heterogeneity was determined. Associations of heterogeneity measures with clinical outcomes were assessed. The analysis identified genes with early mutations that likely play roles in HPV-tumor initiation (FBXW7, NOTCH1, CASP8, FAT1, NSD1, HRAS, EP300, CREBBP, KMT2D). Genes with mutations at later median relative timing (PK3CA, NFE2L2, HLA-A, SPEN, and KMT2C) were more likely to be important for progressing to later stages of the disease. Gain of chromosome arm 3q and loss of arm 11q were found to occur more frequent in HPV+ HNSCC cases. A high prevalence of 3p and 17p loss were also observed in HPV+ HNSCC cases.

  Study phs003139

• Characterization of Human Transcriptome by Computational and HTS Approaches

  We have deeply sequenced total RNA and whole exome of six different post-mortem human snap-frozen tissues (brain, liver, lung, striated muscle, kidney, heart) from three unrelated healthy Caucasian individuals (males, aged 47-54) with the aim to characterize post-transcriptional events due to alternative splicing and RNA editing. In addition, to facilitate the detection of RNA editing sites we have also resequenced the entire genome of each individual. Data are also used to investigate tissue-specific abundance of mitochondrial DNA from exomes and its correlation with mitochondrial transcription, mass and respiratory activity. Tissues were obtained from Cureline (South San Francisco, CA, USA). DNA and RNA were purified using the DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany) and the RNeasy Plus Mini Kit (Qiagen, Hilden, Germany), respectively. Exome capture was performed using the TruSeq Exome Enrichment Kit (Illumina, San Diego, CA) and the sequencing was carried out at IGA Technology Services in Udine (Italy) on Illumina HiSeq 2000 (generating for each tissue approximately 40 millions of 100bp paired-end reads). Strand-oriented RNA reads (2x100bp) were generated using the TruSeq Stranded Total RNA Sample Prep Kit (Illumina, San Diego, CA, USA) and sequenced at IGA Technology Services in Udine (Italy) on Illumina HiSeq 2000 platform. Whole genome resequencing was performed at Personal Genomics in Italy on an Illumina HiSeq 2000 platform (2x100bp) with a mean coverage of 30X. MiRNA-Seq was instead performed at IBBE Institute (Italy) using the TruSeq Small RNA Sample Prep Kit on Illumina MiSeq platform. Preliminary bioinformatics analyses have been performed to check quality using FASTQC software and inconsistent read regions have been trimmed using trim_galore. All reads were mapped onto the human reference genome (version hg19) using GSNAP.

  Study phs000870

• Exome Sequencing to Identify Medically Relevant Associations in Finnish Sub-Isolate Samples from the FINRISK Cohort

  We propose to undertake, at the McDonnell Genome Institute of Washington University, a collaborative exome sequencing project of ~10,000 samples from a large Finnish population cohort, FINRISK, a participating component of the Sequencing Initiative Suomi (SISu) consortium. Together these sequenced samples will provide a unique resource for identification of infrequent or rare variants with a large impact on cardiovascular and metabolic disorders as well as on a wide range of heritable, disease-related quantitative phenotypes. Our collaborative group has already shown, in exome sequencing studies of smaller Finnish samples, that numerous loss of function variants are dramatically more frequent in Finland compared to other European populations, and that such variants demonstrate strong associations to several disease-related traits. Finland is the largest population isolate in Europe and the enrichment of these variants reflects its rapid growth from a severe population bottleneck about 100 generations ago. Northern and Eastern parts of Finland experienced additional, more recent bottlenecks, and therefore the enrichment of particular high-impact variants in these regions is even more extreme (Stoll et al., 2013, PMID: 23912948). We therefore propose to focus exome sequencing efforts on about 10,000 members of FINRISK who are specifically drawn from these sub-isolate regions. As a very wide range of phenotypes are available from this cohort (including cardiovascular and metabolic disease outcomes obtained from national registries) we hypothesize that the proposed exome sequencing will identify numerous new disease-related associations. Within the National Biobanks of Finland (www.nationalbiobanks.fi) we have DNA and tissue samples from over 200,000 individuals of whom 50,000 have GWAS data; these samples provide an exceptional resource for both imputation from the sequencing studies proposed here, as well as for replication of associations.

  Study phs000756

• Chromothripsis in Patient WHIM-09

  We are studying the natural history, pathogenesis and treatment of patients with WHIM syndrome, an immunodeficiency disorder characterized by warts, hypogammaglobulinemia, recurrent infections and neutropenia usually due to autosomal dominant gain-of-function mutations in chemokine receptor CXCR4. We have identified a patient born with WHIM syndrome and the WHIM mutation CXCR4R334X who has been disease-free for 20 years and who lacks CXCR4R334X in myeloid cells, the cells that drive disease manifestations. She is a genetic and hematopoietic mosaic, since she still has the mutation in lymphoid cells and non-hematopoietic cells. Cytogenetics and microarray analysis revealed that the mechanism of loss of the mutation was deletion of the mutant allele from one copy of chromosome 2. Whole genome sequencing of patient neutrophil and skin fibroblast genomic DNA revealed that the mechanism of deletion was chromothripsis, a process of chromosome shattering resulting in deletions and rearrangements of the non-deleted chromosomal segments. In the patient, this process evidently occurred in a single hematopoietic stem cell (HSC), resulting in deletion of the disease allele CXCR4R334X and one copy of 163 other genes on chromosome 2. This HSC evidently acquired a growth advantage and repopulated the HSC population and the myeloid lineage. Consistent with this, studies using gene targeted mice in competitive bone marrow transplantation experiments revealed that selective Cxcr4 haploinsufficiency (inactivation of one copy of Cxcr4 and not of any other genes) was sufficient to confer a strong engraftment advantage over bone marrow cells from wild type mice as well as bone marrow cells from a mouse model of WHIM syndrome. These results suggest that CXCR4 knockdown may be a useful strategy to enhance bone marrow engraftment in the absence of toxic bone marrow conditioning regimens.

  Study phs000856

• CIP: Differential Response to Hydroxyurea and Incidence of Stroke in Sickle Cell Disease

  Sickle cell disease (SCD) is a severe debilitating hematological disorder associated with a high degree of morbidity and mortality. There are approximately 200,000 babies born with sickle cell disease each year, with the disease predominately affecting individuals in Africa. The overall global burden of the disease is tremendous, with more than 100,000 patients currently in the US and further millions worldwide. The governing bodies of the World Health Organization have recently adopted a resolution to strengthen the response to sickle disease in all affected countries and there is a definite need for high quality sickle cell disease research that has the potential to improve the treatment and prognosis of patients with this devastating disease. The clinical manifestations of SCD arise from a complex pathophysiology that includes hemolysis, acute vaso-occlusion, endothelial dysfunction, inflammation, and chronic organ damage. While the individual clinical course of this disease is highly variable, many of the associated complications demonstrate some degree of heritability. Intensive research into identifying genetic modifiers that can affect the pathophysiology of SCD has been limited to date and there is an urgent need to improve of our knowledge the molecular mechanisms underlying the clinical complications of SCD. The Sickle cell CIP project is investigating complication of stroke and pharmacogenomics of hydroxyurea response in patients with sickle cell anemia. The major benefit of hydroxyurea comes from its ability to induce fetal hemoglobin (HbF) and higher HbF levels are associated with reduced morbidity and mortality in SCA patients. We will perform whole exome and whole genome sequencing of SCA patients in order to identify genome variants associated with incidences of stroke and HbF response to hydroxyurea.

  Study phs000691

• Non-coding RNAs Activated by the Wnt/Beta-catenin Signaling Pathway in Hepatoblastoma

  Hepatoblastoma (HB) is the most common pediatric liver tumor, affecting mostly children under 3 years of age. This rare tumor represents 1% of all pediatric cancers. Genetic studies have shown that HB is characterized by high frequency mutations of the CTNNB1 gene encoding beta-catenin (around 75%) and relative genomic stability. Here we have analyzed the transcriptional profile of 21 HBs compared to matched non-tumor livers by Cap Analysis of Gene Expression (CAGE), which provides accurate and quantitative profiling of all transcripts. CAGE analysis revealed strong upregulation of known Wnt target coding genes in most tumors analyzed, consistent with previous transcriptomic studies. To better define the Wnt-dependent transcriptional landscape of HB, we integrated CAGE data with TCF4 ChIP-seq data from a CTNNB1-mutated cancer cell line and with the FANTOM5 genomic coordinates of TCF/LEF binding motifs. Both TCF/LEF binding motifs and ChIP-seq peaks were strongly enriched in the immediate upstream region, not only for protein-coding genes, but also for non-coding transcripts. Among the selected 112 top Wnt target genes at FDR<1.0E-6 and fold change>8, we found clear over-representation (66%) of distant transcription start sites (TSSs) representing lncRNAs and enhancer RNAs, which raises the question of their role in HB pathogenesis. Analysis of the 112 promoters using CAGEd-oPOSSUM confirmed the predominant involvement of Tcf/Lef transcription factors, together with HNF4G, GATA2, Sox3 and Ets-related genes. Finally, the 112 Wnt target signature defined 3 tumor subclasses, T1, T2 and T3, characterized by progressive alteration of the non-coding part of the transcriptome and significant differences in clinical behavior.

  Study phs001433

• Mitochondrial DNA deletion detection in POLG patients

  Accumulation of deletions in the mitochondrial genome is a pathological hallmark of both aging and a broad range of neurometabolic disorders. Prior mtDNA deletion studies have been limited by reliance upon either long range PCR or NGS strategies that target specific regions of the mitochondrial genome. We developed an ultrasensitive high throughput sequencing process (LostArc) to identify and quantify all deletions in circular DNA molecules while minimizing interference from linear DNA and PCR artefacts. We applied this process to skeletal muscle biopsies from 19 unaffected individuals (aged 17 to 93 years) and from 22 patients with heritable mitochondrial disease (17 to 80 years) bearing pathogenic variants in POLG, the nuclear gene encoding the catalytic subunit of the mitochondrial replicative DNA polymerase ɣ (Pol ɣ). To our astonishment, LostArc revealed 35 million mtDNA deletions with ~470,000 unique deletion spans. This unprecedented volume of data enabled detailed and highly informative bioinformatics analyses. We show that ablation of mtDNA is sufficient to explain skeletal muscle phenotypes associated with aging and POLG-derived disease, with a clear threshold separating symptomatic and asymptomatic individuals. Unsupervised hierarchical clustering and principle component analysis reveal distinct age- and disease-correlated deletion patterns that implicate DNA replication by Pol ɣ in the formation of mtDNA deletions. Deletion patterns also point to active mtDNA replication but little or no purifying selection against deleted mtDNA by mitophagy in postmitotic muscle fibers. Monte Carlo simulations of varied hypotheses for deletion formation, based in part on competing models of mtDNA replication, suggest the observed patterns are most consistent with deletions initiated by replication fork stalling during strand displacement mtDNA synthesis. Note that sequence data for HEK293 can be found under BioProject PRJNA642884.

  Study phs002052

• Peripheral Blood Transcriptome Analysis of ALS Patients

  Amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) are neurodegenerative disorders which overlap in their clinical presentation, pathology, and genetic origin. An over-representation of autoimmune disorders exists in both ALS and FTD, but this remains an unexplained epidemiologic observation. Expansions of a hexanucleotide repeat (GGGGCC) in the C9orf72 gene (GeneID:203228) are the most common cause of familial ALS and FTD (C9-ALS/FTD), leading to both repeat containing RNA and dipeptide accumulation coupled with decreased C9orf72 protein expression in brain and peripheral blood cells. Here we show that loss of C9orf72 in myeloid cells is sufficient to recapitulate the age dependent lymphoid hypertrophy and autoinflammation seen in complete C9orf72-/- mice. Dendritic cells isolated from C9orf72-/- mice showed marked early activation of the type I interferon response, and C9orf72-/- myeloid cells were primarily hyperresponsive to activators of stimulator of interferon gene (STING1 GeneID:340061), a key regulator of the innate immune response to cytosolic DNA. STING degradation through the autolysosomal pathway was diminished in C9orf72-/- mice myeloid cells, and blocking STING suppressed hyperactive type I interferon responses in C9orf72-/- myeloid cells and splenomegaly and inflammation in C9orf72-/- mice. Additionally, mice lacking one or both copies of C9orf72 were more susceptible to experimental autoimmune encephalitis, mirroring the susceptibility to autoimmune diseases seen in C9-ALS/FTD patients. Finally, we found that blood derived macrophages, peripheral blood monocytes, and brain tissue from C9-ALS/FTD patients all showed a hyperactive type I interferon signature compared to sporadic ALS/FTD subjects. Collectively, our results indicate that C9-ALS/FTD patients have an altered immunophenotype due to loss of C9orf72 regulation of STING mediated type I interferon activation.

  Study phs002055

• Intraductal transplantation models of human pancreatic ductal adenocarcinoma reveal progressive transition of molecular subtypes

  Pancreatic ductal adenocarcinoma (PDAC) is the most lethal common malignancy, with little improvement in patient outcomes over the past few decades. Recently, subtypes of pancreatic cancer with different prognoses have been elaborated. However, the inability to model these subtypes has precluded mechanistic investigation of their origins. Here, we present an orthotopic xenotransplantation model of human patient-derived PDAC organoids in which organoids are grafted via direct injection into the pancreatic duct of immunocompromised (NSG) mice. Tumors that arise from these intraductally grafted organoids (IGOs) enable the distinction of the two main PDAC subtypes: intraepithelial neoplasms from this model progress in an indolent or invasive manner representing the classical or basal-like subtypes of PDAC, respectively, depending on which organoid they originate from. RNA-sequencing was used to characterize the tumors that resulted from IGO xenografts, and to compare the fast- and slow-progressing IGO tumors. Tumors generated using a traditional orthotopic xenograft model, in which organoids are injected into the parenchyma of the pancreas (i.e. orthotopically grafted organoids, or OGOs), were also sequenced for comparison with the IGO tumors. To further characterize the model, laser-capture microdissection followed by RNA-sequencing was used to compare intraductal to invasive lesions resulting from IGOs. To better understand the connection between KRAS signaling and the fast-progressing phenotype, RNA-sequencing was used to compare cultured organoids in which oncogenic KRAS (GeneID: 3845) was hyperactivated to control organoids, and to compare laser-capture microdissected lesions derived from IGOs of KRAS-hyperactivated organoids to lesions derived from control organoids. To ascertain copy number, low-coverage, whole genome DNA-sequencing was performed on organoids cultured in vitro, and on laser-capture microdissected lesions derived from IGO tumors.

  Study phs002045

• A novel recurrent mutation in MITF predisposes to familial and sporadic melanoma

  We conducted whole-genome sequencing of probands from several melanoma families, identifying one individual carrying a novel germline variant (c.G1075A, NM_000248.3; p.E318K, NP_000239.1; rs149617956) in the melanoma lineage-specific oncogene MITF. While the variant cosegregated with melanoma in some, but not all cases in the family, linkage analysis of 31 families subsequently identified to carry the variant generated a LOD score of 2.7 under a dominant model, suggesting E318K as a possible intermediate risk variant. Consistent with this, E318K was significantly associated with melanoma in a large Australian case-control sample, giving an odds ratio (OR) of 2.33, 95% CI 1.21-4.70 (case and control carrier frequency, 0.0165 and 0.0072, respectively; P=0.008). Likewise, it was similarly associated in an independent case-control sample from the United Kingdom (UK P=0.012; combined P=0.0003, OR 2.19, 95% CI 1.41-3.45). In the Australian sample, the variant allele was significantly over-represented in cases with a family history of melanoma (OR 2.95, 95% CI 1.23-6.92), multiple primary melanomas (OR 4.22, 95% CI 1.52-10.91), or both (OR 8.37, 95% CI 2.58-23.80). The variant allele was also associated with increased nevus count (combined P=0.002, OR 2.54, 95% CI 1.42-4.55) and non-blue eye color (combined P=0.008, OR 2.01, 95% CI 1.11-3.81). Functional analysis of E318K showed that MITF encoded by the variant allele had impaired sumoylation and differentially regulated several MITF targets. These data indicate that MITF is a melanoma predisposition gene and highlights the utility of whole-genome sequencing to identify novel rare variants associated with disease susceptibility.

  Study phs000419

• Resistance to Checkpoint Blockade Therapy Through Inactivation of Antigen Presentation

  Although immune checkpoint blockade (CPB) leads to prolonged responses in 15-40% of patients with metastatic melanoma, treatment refractory disease and progression after initial response remain major causes of mortality. While predictors of response have been reported, the common mechanisms of both primary and acquired resistance are poorly understood. To identify mechanisms of resistance and examine the evolving landscape in response to CPB, we performed whole exome sequencing (WES), immunohistochemistry (IHC), and RNA-sequencing (RNAseq) of longitudinal tumor biopsies from 17 metastatic melanoma patients treated with various CPB therapies. We found no significant changes in both mutational and neoantigen loads over time between responders and nonresponders. However, we identified abnormalities in one gene, beta-2-microglobulin (B2M), an essential component of MHC Class I antigen presentation, that were present in samples during disease progression but not regression. In total, we identified B2M aberrations in 29.4% of patients, including multiple early frameshift mutations, loss of heterozygosity (LOH) overlapping B2M, and absence of tumor-specific B2M protein expression. Additional defects in the antigen presentation and IFNγ pathways were identified but were not restricted to progressing lesions in our cohort. In two independent cohorts of 105 and 38 melanoma patients treated with ipilimumab (anti-CTLA4) and pembrolizumab (anti-PD1) respectively, we found that B2M LOH was enriched 3-fold in nonresponders (~30%) vs. responders (~10%) and associated with poorer overall survival (log-rank p=0.01, p=0.006). Loss of both copies of B2M was found only in nonresponders. We also found evidence for association of LOH overlapping IFNGR1 with poorer overall survival exclusively in the anti-PD1 cohort. Thus, B2M loss is likely a common mechanism of primary and acquired resistance to therapies targeting CTLA4 or PD-1.

  Study phs001427

• Single-cell RNA-Seq of human lymphoma reveals malignant B cell diversity and patterns of T cell immune checkpoint co-expression

  Follicular lymphoma (FL) is a generally incurable B-cell malignancy which has the potential to transform into highly aggressive lymphomas. Genomic studies indicate it is often a small subpopulation rather than the dominant population in the FL that gives rise to the more aggressive subtype. To resolve the underlying transcriptional networks of follicular B-cell lymphomas at single molecule and cell resolution, we leveraged droplet-based barcoding technology for highly parallel single cell RNA-Seq. We analyzed the transcriptomes from tens of thousands of cells derived from five primary FL tumors. Simultaneously, we conducted multi-dimensional flow cell sorting to validate our characterizing of cellular lineages and critical expressed proteins. For each tumor, we identified multiple cellular subpopulations, matching known hematopoietic lineages. Comparison of gene expression by matched malignant and normal B cells from the same patient revealed tumor-specific features. Malignant B cells exhibited restricted immunoglobulin light chain expression (either Ig Kappa or Ig Lambda), as well the expected upregulation of the BCL2 gene, but also down-regulation of the FCER2, CD52 and MHC class II genes. By leveraging the single-cell resolution on large numbers of cells per patient, we were able to examine tumor-resident T cells. We identified pairs of immune checkpoint molecules that were co-expressed, providing a potentially useful strategy for selection of patient-tailored combination immunotherapies. In summary, massively parallel measurement of single-cell expression in thousands of tumor cells and tumor-resident lymphocytes can be used to obtain a systems-level view of the tumor microenvironment and identify new avenues for therapeutic development.

  Study phs001378

• Investigating the Role of Neddylation in the Repair of Topoisomerase I-Mediated DNA Damage in Colorectal Cancer

  Colorectal cancer is one of the most lethal cancers worldwide. First-line chemotherapy for metastatic CRC (mCRC) entails a combination of irinotecan with 5-fluorouracil and leucovorin (FOLFIRI). Irinotecan is a prodrug that is converted within the cell to its active metabolite SN38, a potent TOP1 inhibitor. TOP1 relieves DNA torsional strain arising from DNA metabolism by cleaving one strand of the DNA duplex and forming an enzyme-DNA covalent intermediate that can be trapped by inhibitors. The resulting TOP1 DNA-protein crosslinks (TOP1-DPCs) interfere with replication and transcription and induce cell death if left unrepaired. I have established the role of the ubiquitin-proteasome pathway in TOP1-DPC repair, but it remains unknown how ubiquitylation is activated. One possibility is that collision between advancing replication forks and the DPCs triggers a specific ubiquitin ligase to target DPCs for degradation. The DDB1-CUL4 (Cullin 4)-RBX1 (CRL4) ubiquitin ligase complex performs a salient role in both replication and DNA repair upon its activation by auto-NEDD8 modification (neddylation), but it is unknown whether it plays a role for TOP1-DPC repair. We plan to define the role of neddylation for TOP1-DPC repair during replication and examine irinotecan plus pevonedistat (PEV), a NEDD8-activating enzyme inhibitor in multiple clinical trials, in CRC patient-derived preclinical models including 3D organoids. The three in-house organoids, from which samples have been submitted to dbGaP, were sequenced as a means of verifying their identity as colon cancer organoids, to allow us to investigate the combination of TOP1 inhibitors and PEV and the molecular mechanisms underlying this drug combination in these organoids.

  Study phs003257

• A Polygenic Score for Acute Vaso-Occlusive Pain in Pediatric Sickle Cell Disease

  Recurrent acute pain, or vaso-occlusive pain crisis (VOP), is the most common complication of sickle cell disease and correlates strongly with increased hospital visits and early mortality (PMID: 1710777). While the genetics of VOP in sickle cell patients has been studied (PMID: 29205277, 27883292, 25102390, 30079801, 22925497, 29531649, 29620434, 19468207, 20172753, 22576309, 30031848, 24136375, 27603703, 29559808), it is not fully understood.Using whole genome sequencing, we interrogated the α-thalassemia deletion -α3.7 and 133 candidate single-nucleotide polymorphisms (SNPs) across 66 genes for associations with VOP in 327 SCCRIP participants followed longitudinally over 6 years. Twenty-one SNPs in 9 loci were associated with VOP, including 3 (BCL11A, MYB, and the β-like globin gene cluster) that regulate erythrocyte fetal hemoglobin (HbF) levels and 6 (COMT, TBC1D1, KCNJ6, FAAH, NR3C1, and IL1A) that were associated previously with various pain syndromes. An unweighted PGS integrating all 21 SNPs was associated with the VOP event rate (estimate, 0.35; standard error, 0.04; P = 5.9 × 10−14) and VOP event occurrence (estimate, 0.42; standard error, 0.06; P = 4.1 × 10−13). These associations were stronger than those of any single locus. Our findings provide insights into the genetic modulation of VOP in children with SCD. More generally, we demonstrate the utility of WGS for investigating genetic contributions to the variable expression of SCD-associated morbidities.For publications, written, and/or oral presentations resulting from use of this data, please cite the following information: "The original data description for A Polygenic Score for Acute Vaso-Occlusive Pain in Pediatric Sickle Cell Disease can be found here: (Journal Citation Blood Adv. 2021 Jul 27;5(14):2839-2851) (PMID: 34283174)."

  Study phs002470

• NHLBI GO-ESP: Lung Cohorts Exome Sequencing Project (Asthma): Genetic variants affecting susceptibility and severity

  The NHLBI "Grand Opportunity" Exome Sequencing Project (GO-ESP), a signature project of the NHLBI Recovery Act investment, was designed to identify genetic variants in coding regions (exons) of the human genome (the "exome") that are associated with heart, lung and blood diseases. These and related diseases that are of high impact to public health and individuals from diverse racial and ethnic groups will be studied. These data may help researchers understand the causes of disease, contributing to better ways to prevent, diagnose, and treat diseases, as well as determine whether to tailor prevention and treatments to specific populations. This could lead to more effective treatments and reduce the likelihood of side effects. GO-ESP is comprised of five collaborative components: 3 cohort consortia - HeartGO, LungGO, and WHISP - and 2 sequencing centers - BroadGO and SeattleGO. The exome sequencing asthma project includes 200 African-Americans with asthma from the NHLBI multicenter Severe Asthma Research Program (SARP). SARP participants were recruited at the NHLBI SARP sites with an emphasis on recruiting severe asthmatics (Moore et al., Am J Respir Crit Care Med, 2010. PMID: 19892860). Asthma status was based on both a physician's diagnosis and either bronchodilator reversibility or hyper-responsiveness to methacholine as well as less than 5 pack years of smoking. All subjects were carefully characterized using the standardized SARP protocol which included spirometry (medication withheld), maximum bronchodilator reversibility, hyper-responsiveness to methacholine (not performed in subjects with low baseline FEV1), skin-tests to common allergens, questionnaires on health care utilization and medication use and sputum, lung imaging and bronchoscopy in a subset. In addition GWAS data are available (phs000355, Illumina platform).

  Study phs000422

• The Human Pancreas Analysis Program (HPAP)

  The past decade has seen a dramatic improvement in our ability to phenotype and molecularly profile human tissues with unprecedented resolution at the genomic, epigenomic, protein, and functional levels. The Human Pancreas Analysis Program (HPAP), part of the Human Islet Research Network and supported by the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) through multiple NIH grants, is performing deep phenotyping of the human endocrine pancreas to better understand the cellular and molecular events that precede and lead to beta-cell loss and/or dysfunction in type 1 diabetes (T1D) and type 2 diabetes (T2D) as well as to accumulate, analyze, and distribute high-value data sets to the diabetes research community through the HPAP PANC-DB database (and additional details provided at that site). To this end, HPAP employs state-of-the-art technologies to perform comprehensive analyses of pancreas biology as it pertains to organ donors with T1D, autoantibody-positive donors without diabetes, donors with T2D, and control donors. Pancreas procurement and analyses take advantage of the expertise and extensive network of Organ Procurement Organizations (OPOs) and autoantibody screening centers established through the JDRF/The Lenona M. and Harry B. Helmsley Charitable Trust (HCT) – funded nPOD program. In contrast to nPOD, the major product provided by HPAP is not archived biomaterial subject to broad distribution, but rather, the delivery of extensive and high-quality molecular data sets to the diabetes research community in order to facilitate further discovery. Together, nPOD and HPAP are complementary programs that assist the diabetes community and afford the maximal opportunity for advancing knowledge about the pathogenesis of T1D and T2D.

  Study phs002465

• AC-ICAM: An Atlas and Compass of Immune-CAncer-Microbiome Interactions in Colon Cancer

  Fresh frozen tumor samples and matched normal colon specimens from chemotherapy naïve, consecutive patients with a histological diagnosis of colon adenocarcinoma at Leiden University Medical Center (LUMC) were subjected to multi-omics profiling at Sidra Integrated Genomics Services, part of Sidra Medicine, Doha, Qatar. This resulted in a high-quality dataset of samples collected between 2001 and 2015. Multi-step QC was performed, including assessment of the quality of the nucleic acid extraction, sequencing, and genetic concordance between DNA/RNA sequences from the same patient. After QC, this resulted in RNA sequencing on 348 tumor samples, whole exome sequencing (WES) on 281 tumor-normal sample pairs (n=562 samples), TCR sequencing of the TCR β locus (TRB) for 114 tumor and 10 matched normal samples, and 16S rDNA sequencing on 246 tumor-normal sample pairs (n=492 samples) and 42 tumor only samples, and 167 WGS non human-aligned reads of the tumor samples. Extensive clinicopathological and survival data of the cohort are available. The median follow-up time for this cohort is 4.6 years, and the distribution of gender and pathological stage are within expected ranges for this patient population with a mean age of 68. We refer to this resource as the Sidra-LUMC (SILU) colon adenocarcinoma (COAD) Atlas and Compass of Immune-Colon cancer-Microbiome interactions (AC-ICAM)."Reprinted from Nature Medicine,Roelands, J., Kuppen, P.J.K., Ahmed, E.I. et al. An integrated tumor, immune and microbiome atlas of colon cancer. Nat Med (2023). https://doi.org/10.1038/s41591-023-02324-5, with permission from Nature."

  Study phs002978

• Genetic Aberrations and Subclonal Structure Impact Chronic Lymphocytic Leukemia

  Large-scale whole-exome sequencing (WES) of primary tumor samples enables the unbiased discovery of recurrent putative driver events and patterns of clonal evolution. We report the identification of 44 recurrently mutated genes and 11 recurrent CNVs through the WES of 538 chronic lymphocytic leukemia (CLL) and matched germline DNAs. These include previously unrecognized cancer drivers (e.g., RPS15, IKZF3), and collectively identify nuclear export, MYC activity and MAPK signaling as central pathways affected by somatic mutation in CLL. A clonality analysis of this large dataset further enabled the reconstruction of temporal relationships between these driver events. Several drivers were associated with shorter progression-free survival (PFS) in 280 samples that were collected prior to uniform treatment with front line chemo-immunotherapy, with mature follow up of greater than 10 years. Direct comparison between matched pretreatment and relapse CLL from 59 samples demonstrated marked clonal evolution occurring in more than 95% of these patients. Distinct patterns of clonal evolution in relationship to specific gene alteration were observed, suggesting a hierarchy of fitness amongst mutations. Thus, large WES datasets of clinically informative samples enable the discovery of novel driver genes as well as the network of relationships between the drivers and their impact on disease relapse and clinical outcome.Additionally, we performed RNA-seq for 268 CLL samples (including 26 follow-up samples) and used them to identify expression subtypes of CLL. RRBS for 30 of these samples was also generated. In an integrative analysis of genetic, transcriptomic, and epigenetic data, incorporating known and newly identified subtypes of CLL, we built new models to improve patient prognostication.

  Study phs000922

• Center for Inherited Disease Research (CIDR)-National Institute on Aging (NIA) Whole Exome Analysis of Ehlers-Danlos Syndrome

  The Ehlers-Danlos syndromes (EDS) are a family of hereditary disorders of connective tissue characterized by joint hypermobility and skin involvement with soft, translucent, stretchy skin, and atrophic scarring. To date, over 12 different types of EDS are known, and the genes causing almost all of them have been found. However, the gene(s) underlying the most common of Ehlers-Danlos syndrome, known as the hypermobile type, remain elusive. The hypermobile type of Ehlers-Danlos syndrome is associated with multiple co-morbidities, including dysautonomia, chronic musculoskeletal pain, gastrointestinal dysmotility, mast cell activation and multiple neurologic complications. Affected persons often suffer from chronic fatigue and may be chronically disabled. While the precise frequency of the condition is not known, it may be as common as between 1/500 and 1/1000 persons. There is a significant overlap with the phenotype of fibromyalgia, and multiple persons diagnosed with fibromyalgia may well have the hypermobile type of Ehlers-Danlos syndrome. The phenotype of hypermobile EDS segregates as a single dominant trait with complete penetrance, variable expressivity, and a markedly evolving phenotype over time. Identification of the genetic cause of this hereditary disorder has the potential to inform rational therapy for this disabling condition. It is hoped that the insight into the pathogenesis of the hypermobile type of EDS will lead to a similar breakthrough in therapy for this disorder. Toward this end, whole exome sequencing (WES) of DNA from a total of 167 patients will be performed at the Center for Inherited Disease Research (CIDR) to identify genes underlying the hypermobile type of Ehlers-Danlos syndrome and those with EDS of an uncharacterized type.

  Study phs001779

• Genotyping NIGMS CEPH Samples from the United States, Venezuela, and France

  The Human Genetic Cell Repository (HGCR) is sponsored by the National Institute of General Medical Sciences (NIGMS) with the mission of supplying scientists with the materials for accelerating disease gene discovery and functional studies. The resources available include highly-characterized, contaminant-free fibroblast cell lines, transformed lymphoblastoid cell lines (LCLs), and DNA samples derived from these cultures. While the Repository has a major emphasis on inherited diseases, it also contains large collections of cell lines gathered from populations around the world that are intended for use in facilitating the understanding of human variation. Prominent among these resources is the CEPH Reference Family collection contributed by the Centre d'Etude du Polymorphisme Humain (CEPH), Foundation Jean Dausset, Paris, France. The CEPH collection includes families collected by R. White (Utah), J. Dausset (French), J. Gusella (Venezuelan), and J. Egeland (Amish). There are a total of 809 individuals accounting for 832 pedigree positions in the reference families. Family relationships for the 61 reference families were verified at Coriell and approved by CEPH. In an effort to enhance the value of this cell culture resource available from the Repository, the Coriell Genotyping and Microarray Center used the Affymetrix Genome-Wide Human SNP Array 6.0 platform to genotype 181 samples from the NIGMS HGCR CEPH collection. Included are thirteen families from the United States (Amish Pedigree 884 and Utah Pedigrees 1331, 1356, 1400, 1416, 1424, 1427, 1477, and 1582), France (Pedigrees 35 and 66), and Venezuela (Pedigrees 102 and 104). Twelve of the families consist of 3 generations and one family consists of 2 generations. Pedigree charts and sample descriptions are available on the Repository catalog website (see below).

  Study phs000268

• The ClinSeq Project: Piloting Large-Scale Genome Sequencing for Research in Genomic Medicine

  ClinSeq is a pilot project to investigate the use of whole-genome sequencing as a tool for clinical research. By piloting the acquisition of large amounts of DNA sequence data from individual human subjects, we are fostering the development of hypothesis-generating approaches for performing research in genomic medicine, including the exploration of issues related to the genetic architecture of disease, implementation of genomic technology, informed consent, disclosure of genetic information, and archiving, analyzing, and displaying sequence data. In the initial phase of ClinSeq, we are enrolling roughly 1000 participants; the evaluation of each includes obtaining a detailed family and medical history, as well as a clinical evaluation. The participants are being consented broadly for research on many traits and for whole-genome sequencing. Initially, Sanger-based sequencing of 300-400 genes thought to be relevant to atherosclerosis is being performed, with the resulting data analyzed for rare, high-penetrance variants associated with specific clinical traits. The participants are also being consented to allow the contact of family members for additional studies of sequence variants to explore their potential association with specific phenotypes. Here, we present the general considerations in designing ClinSeq, preliminary results based on the generation of an initial 826 Mb of sequence data, the findings for several genes that serve as positive controls for the project, and our views about the potential implications of ClinSeq. The early experiences with ClinSeq illustrate how large-scale medical sequencing can be a practical, productive, and critical component of research in genomic medicine. Reprinted from Genome Res. 2009 Sep; 19(9): 1665-1674, with permission from Genome Research, PMID: 19602640.

  Study phs000971

• Epigenetic Profiling of Human Colorectal Cancer

  DNA methylation, together with chromatin modifications, constitute the epigenome that functions to regulate gene expression and genome integrity. DNA methylation alterations are ubiquitous in human cancers, as many genes acquire DNA methylation in a cancer-specific manner. DNA methylation at these sites in the genome of cancer cells not only serves as a marker for tumor identification, but together with gene mutation and gene expression data, can also be used to describe subsets of tumors of the same organ source. We have previously shown distinct human colorectal cancer subtypes based on DNA methylation differences. Correlating these DNA methylation differences with clinical co-variates will serve to further understand how these distinct subtypes are generated. We have collected 100 colorectal tumor tissues (for which clinical information is known) and have obtained (unprotected) genome-wide DNA methylation and chromatin modification information for each sample for the purposes of identifying and classifying unique tumor subtypes of colorectal cancers. In addition, we have determined mutations of key genes relevant to colorectal cancer as well as gene expression profiles. We will use the clinical data for each de-identified sample to correlate with the DNA methylation, mutation and gene expression information so as to understand the driving forces behind these distinct colorectal subtypes. We have selected the most promising tumors for whole-genome bisulfite sequencing using next-generation sequencing technology to obtain complete maps of colon cancer methylomes. Researchers will be unable to identify the subjects because the samples and associated information have been de-identified and anonymized by the tissue source site. In addition, upon receipt by the USC Epigenome Center, we have assigned new random identifiers for each sample. The data generated using these new codes are not traceable to the patient identity.

  Study phs000385

• Sudden Death in the Young Case Registry

  Knowledge gaps persist about the incidence of and risk factors for sudden death in the young (SDY). The SDY Case Registry is a collaborative effort between the National Institutes of Health, the Centers for Disease Control and Prevention, and the Michigan Public Health Institute. Its goals are to: (1) describe the incidence of SDY in the United States by using population-based surveillance: (2) compile data from SDY cases to create a resource of information and DNA samples for research; (3) encourage standardized approaches to investigation, autopsy, and categorization of SDY cases; (4) develop partnerships between local, state, and federal stakeholders toward a common goal of understanding and preventing SDY; and (5) support families who have lost loved ones to SDY by providing resources on bereavement and medical evaluation of surviving family members. Built on existing Child Death Review programs and as an expansion of the Sudden Unexpected Infant Death Case Registry, the SDY Case Registry achieves its goals by identifying SDY cases, providing guidance to medical examiners/coroners in conducting comprehensive autopsies, evaluating cases through child death review and an advanced review by clinical specialists, and classifying cases according to a standardized algorithm. The SDY Case Registry also includes a process to obtain informed consent from next-of-kin to save DNA for research, banking, and, in some cases, diagnostic genetic testing. The SDY Case Registry will provide valuable incidence data and will enhance understanding of the characteristics of SDY cases to inform the development of targeted prevention efforts. For cases where informed consent has been obtained for research genomic sequencing has been performed. On these cases genomic data, phenotype information and on a limited number of cases redacted autopsy reports will be available through dbGaP.

  Study phs003221

• Clinical and Molecular Features of Acquired Resistance to Immunotherapy in Non-Small Cell Lung Cancer

  PD-(L)1 blockade has been shown to promote durable responses in non-small cell lung cancer (NSCLC) patients and has been rapidly incorporated into the treatment regimens for many with advanced disease. Despite this, mechanisms of acquired resistance (AR) to PD-(L)1 have not been well understood. To address the clinical and molecular landscape of AR to PD-(L)1 blockade in NSCLC patients, a large clinical cohort of AR to PD-(L)1 blockade in lung cancer was examined. In conjunction, a systematic genomic and transcriptomic analysis in a subset of patients (n = 29) with available tissue samples was conducted. To investigate the molecular mechanisms of AR to PD-1 blockade from pre- and/or post-treatment tumors, microarray-based whole transcriptome expression data for 29 tumor samples, and whole exome sequencing (WES) data from 34 tumor samples was generated from this subset. Microarray data was uploaded to the Gene Expression Omnibus (GEO) database. Most AR tumors showed retained or increased inflammatory characteristics with significant upregulation of IFNγ which was suggestive of persistent albeit insufficient anti-tumor immune response. This finding expounds upon the prior pre-clinical and translational data supporting the intricate role of IFNγ in sensitivity and resistance to immunotherapy. Whereas initial IFNγ exposure may be fundamental to T cell activation (a hallmark of immune response), persistent IFNγ related effects and upregulation may potentially signal immune dysfunction and IFNγ insensitivity in contrast. The inflammatory phenotypes identified have implications for future rational development of new immunotherapy strategies for patients with AR.

  Study phs002834

• Genomic Answers for Kids (GA4K)

  Identification of complex genetic variants including defects in gene regulatory circuits and uncharacterized genes present challenges for rare disease diagnosis. In Genomic Answers for Kids program we apply the joint interpretation of patient genome sequence with genomic endophenotypes to expand the clinically actionable genome among pediatric cases of suspected genetic disease. Our objective is to examine the role of novel and under-interpreted sequence variation by combining complete DNA sequences, personal epigenetic variations and massively parallel functional screens. To achieve this we perform augmented whole genome sequence (WGS) interpretation including reanalysis of clinical exomes as well as generating WGS for patients with negative exome results all be subject to family-based semi-automated recall pipeline to discover missed diagnostic variation. We use linked-read and long-read sequencing technologies to focus on putative structural variants missed in short-read genome and exome analysis by the optimized integration of linked and long read technologies that also improve identification of transmission patterns of all variants, as well as resolving genomic regions resistant to standard alignment. We also include tissue DNA sequencing to investigate somatic mosaicism. To further assist with WGS interpretation for uncommon variation we capture snapshots of patient transcriptomes and epigenomes in individual cells using single-cell RNA (scRNA) and sc open chromatin (scATAC) as well as bulk whole genome bisulphite genome sequencing for methylome interpretation. Alternative splicing is functionally assessed in available patient tissues using RNA-seq, including full length cDNA sequences by IsoSeq (PacBio) methodology. Our overarching goal is an increase rate of diagnostic genomic findings up to two-fold among rare disease, resulting in the majority (>50%) of patient cases resolved by the integrated system developed in our program.

  Study phs002206

• Blood Transcriptome Profiling Following Seizures

  The retrospective diagnosis or confirmation of a seizure, after the event, is very challenging. Many times clinicians can only use a patient's self-reported history. To further complicate the diagnosis, it is known that 20-30% of patients reporting a seizure actually suffer from psychogenic non-epileptic seizure (PNES). The patient has the behavior of a seizure, but does not suffer from brain seizure activity. Considerable time and resources are used to monitor PNES, and many patients undergo unnecessary therapy with anti-epileptic drugs. These resources could be more effectively used if we could confirm a seizure after the event (retrospective diagnosis). Here we propose a novel approach to distinguish seizures from PNES after the seizure event has happened.We have developed an accurate blood biomarker technology to diagnose acute brain injury, and assess chronic progression following brain injury. Previous studies support the premise that seizures affect blood RNA expression; here we assess temporal RNA expression pattern changes following seizure. We present preliminary preclinical data showing the discrimination of electric-stimulation evoked animal seizures. Here, we validate this maturing technology to assist with the discrimination of epileptic seizures from psychogenic non-epileptic seizures. We will analyze blood samples from patients undergoing video EEG monitoring in an epilepsy monitoring unit (Emory University Hospital, Atlanta GA). We expect 20% of these patients will have PNES. Once a patient has an EEG verified seizure, we will collect blood at various time points and identify biomarker RNA molecules in the blood. Blood samples will be collected from patients with non-verified seizures as control. These RNA expression values will be used to develop a predictive test for the occurrence of seizures. A biomarker gene panel will be developed from this exploratory project for validation in a larger patient population.

  Study phs003460

• Genetic Investigations of Attention-Deficit/Hyperactivity Disorder

  Design and Aims:We conducted whole-exome DNA sequencing on 152 parent-child trios (456 individuals total) that included a child meeting DSM-IV or DSM-5 criteria for ADHD along with both biological parents to identify rare de novo damaging DNA variants associated with ADHD risk. Mutation rates were compared between 147 ADHD trios and 780 unaffected control trios. We aimed to identify specific ADHD risk genes by combining the de novo trio data with rare variant data from an independent case-control cohort.Population Information:ADHD trios were collected from four sites. The 780 control trios were unaffected siblings from the Simons Simplex Collection who scored in the normal range on the Child Behavior Checklist (CBCL) or Adult Behavior Checklist (ABCL) ADHD subscale.Molecular Technologies:Whole-exome DNA sequencing was performed on the ADHD trios using the IDT xGen V1 or Agilent SureSelect All Exon V7 capture and the Illumina NovaSeq6000 sequencer. Control trios were previously sequenced with the NimbleGen SeqCap EzExomeV2 capture and Illumina HiSeq 2000 sequencer. Joint variant calling was performed for all trios. Rare de novo variants with allele frequency Principal Findings:We found a 1.67-fold enrichment of rare de novo damaging variants and a 1.93-fold enrichment of ultra-rare de novo damaging variants in ADHD cases vs controls. By integrating the de novo data with an independent case-control cohort, KDM5B was identified as a high-confidence ADHD risk gene (FDR=0.04). Three potential risk genes were also identified. Twenty-three genes with ultra-rare damaging de novo variants in ADHD overlapped with previously reported risk genes for other neuropsychiatric conditions. Exploratory analyses showed enrichment of these genes in early neurodevelopmental pathways.

  Study phs003647

• Chromatin Landscape of BET Inhibitor-Treated CD8+ T-cells from Chronic Lymphocytic Leukemia Patients

  The pan-BET (bromodomain and extra-terminal) inhibitor, OPN-51107 (OPN5) was found to alleviate features of T-cell exhaustion in models of chronic lymphocytic leukemia (CLL). To evaluate if OPN5 treatment could relieve fixed chromatin states associated with terminal T-cell dysfunction, ATAC-sequencing was performed on CLL patient T-cells. Peripheral blood mononuclear cells (PBMCs) from 4 treatment-naïve CLL patients were treated ex vivo with 0.5 μM OPN5 or equivalent vehicle (DMSO) in the presence of T-cell activating stimuli (10 μg/mL anti-CD3 and 5 μg/mL anti-CD28) for 16 hours. Live CD8+ T-cells were then isolated via fluorescence-activated cell sorting, and approximately 100,000 viable CD8+ T-cells were subjected to ATAC-sequencing. Libraries were prepared using the Active Motif ATAC-seq kit and underwent paired-end sequencing (2 x 50 bp, 100 cycles) on an Illumina NovaSeq 6000 SP sequencer. Analysis of this data revealed that OPN5 treatment induces a unique chromatin landscape in CLL patient T-cells that is less like that of terminally exhausted T-cells and more like that of effector and memory T-cells. In particular, OPN5 treatment reduced chromatin accessibility at promoter regions of terminal exhaustion-associated genes (e.g., CD101, CCR5, LAG3, BATF3) and increased chromatin accessibility at promoter regions of genes associated with progenitor T-cell function (e.g., TCF7, IL2, CXCR5, IFNG). Individual raw sequencing (fastq) files and processed bigwig files (aligned to the hg38 human reference genome), as well as bigwig files for grouped OPN5- and vehicle-treated samples (aggregate of all 4 patient samples per treatment group) are available via authorized access request through dbGaP.

  Study phs003613

• Toxigenic Clostridium perfringens isolated from at-risk pediatric inflammatory bowel disease patients

  With the goal of studying mucosal-associated bacteria and products detrimental to inflammatory bowel disease (IBD) patients, we directly cultured bacterial communities from mucosal biopsies and examined for detrimental traits. Focusing on red blood cell lysis, a virulence trait, communities from 5 of 9 patients exhibited hemolytic activity. Clostridium perfringens was isolated from each hemolytic community, and whole-genome sequencing confirmed presence of hemolytic toxins – alpha toxin in all 5 strains, perfringolysin O (PFO) in 4 strains, alveolysin in 2 strains – and absence of enterotoxins. In vitro cytotoxicity experiments demonstrated strain supernatants were toxic to various cell types located beneath the intestinal epithelial barrier, including endothelial, neuroblasts, and neutrophils, while impact on epithelial cells was less pronounced. These findings suggest C. perfringens may be damaging particularly when intestinal barrier integrity is compromised, such as in IBD patients. Further characterization using purified toxins and genetic insertion mutants confirmed PFO toxin was sufficient for cell toxicity. C. perfringens supernatants also induced activation of neuroblast and dorsal root ganglion (DRG) neuron cells in vitro, suggesting C. perfringens in inflamed mucosal tissue may directly contribute to abdominal pain, a frequent IBD symptom. Strikingly, C. perfringens was detected in biopsies of 8 of 9 pediatric patients while not in all corresponding stool samples, and in other datasets, C. perfringens was in stool samples of IBD adults (18.7-27.1%) much more frequently versus healthy (5.1%). Our findings showcase C. perfringens may have an important pathogenic impact on IBD patients. Similar to Clostridioides difficile, routine monitoring of C. perfringens and PFO toxins in patients may be beneficial.

  Study EGAS50000000334

• METABRIC

  Solid tumors are complex tissues composed of a mixture of cancer and normal cells, which complicates the interpretation of their molecular profiles. Normal cell contamination can dilute cancer cell information and tissue architecture is generally not reflected in molecular assays. To address these challenges, we developed a computational approach based on standard Haematoxylin and Eosin-stained sections and demonstrated its power in a discovery cohort of 323 breast tumors and an independent validation cohort of 241 tumors. First, to deconvolute cellular heterogeneity and detect subtle genomic aberrations, we introduced an algorithm based on tumor cellularity to increase the comparability of copy-number profiles between samples. Second, we demonstrated that a predictor for survival integrating image-based and molecular features significantly outperforms classifiers based on single data types. Third, we described and validated a novel, independent prognostic factor based on quantitative analysis of spatial patterns between stromal cells, which are not detectable by molecular assays. Our quantitative methods refine and complement molecular assays of tumor samples and could benefit all large-scale cancer studies.First, to deconvolute cellular heterogeneity and detect subtle genomic aberrations, we introduced an algorithm based on tumor cellularity to increase the comparability of copy-number profiles between samples. Second, we demonstrated that a predictor for survival integrating image-based and molecular features significantly outperforms classifiers based on single data types. Third, we described and validated a novel, independent prognostic factor based on quantitative analysis of spatial patterns between stromal cells, which are not detectable by molecular assays. Our quantitative methods refine and complement molecular assays of tumor samples and could benefit all large-scale cancer studies.

  Study EGAS00000000098

• Cell-Free DNA Genomic and Fragmentomic Features for Early Outcome Prediction in Diffuse Large B-Cell Lymphoma

  Purpose: Diffuse large B-cell lymphoma (DLBCL) patients need an accurate and early risk stratification strategy as prompt therapy escalation may improve outcome. Patients and methods: We evaluated cell free DNA (cfDNA) genomic and fragmentomic features in 190 patients with large B-cell lymphoma from more than 40 hospitals in the Benelux using Whole Genome Sequencing (WGS). We defined a new metric called ACT score (Aberrations, Contribution of short fragments, Terminal motif analyses) by combining four cfDNA signals and evaluated its performance to predict end-of-treatment (EOT) response and survival. Results: Individual cfDNA features and ACT score were altered per EOT response after one cycle of treatment (T1). The ACT score at T1 outperformed all individual cfDNA features in predicting EOT response (AUC = 0.74). Patients with a positive ACT score had an inferior outcome compared to ACT score negative patients [progression-free survival (PFS): HR 5.0 (95% CI 3.1-8.2), log-rank test, p < 0.0001) and overall survival (OS): HR 6.5 (95% CI 3.5-12.2), log-rank test, p < 0.0001]. The 2-year PFS in ACT score positive and negative patients were 29.8% and 78.1%, respectively. The prognostic value of the ACT score is independent of the International Prognostic Index and interim imaging. Conclusions: ACT score computed from a single plasma sample collected after one cycle of treatment can predict clinical outcome. This low-cost and easy-to-interpret test does not require tissue biopsies or a priori knowledge of mutations and has the potential to aid patient selection in interventional clinical trials and risk-adapted treatment strategies.

  Study EGAS50000000412

• Luminal breast epithelial cells from wildtype and BRCA mutation carriers harbor copy number alterations commonly associated with breast cancer

  Cancer-associated mutations have been documented in normal tissues, but the prevalence and nature of somatic copy number alterations and their role in tumor initiation and evolution is not well understood. Here, using single cell DNA sequencing, we describe the landscape of CNA variants in >48,000 breast epithelial cells from women with normal or high risk of developing breast cancer. Individual cells with one or two of a specific subset of CNAs (e.g. 1q gain and 16q, 22q, 7q, and 10q loss) are detectable in almost all breast tissues and, in those from BRCA1 or BRCA2 mutations carriers, occurs prior to loss of heterozygosity (LOH) of the wildtype alleles. These CNAs, which are among the most common associated with ductal carcinoma in situ (DCIS) and malignant breast tumors, are enriched almost exclusively in luminal cells and are generally absent from basal myoepithelial cells. Allele-specific analysis of the enriched CNAs reveals that each allele was independently altered, demonstrating convergent evolution of these CNAs in an individual breast. Tissues from BRCA1 or BRCA2 mutation carriers contained a small percentage of cells with extreme aneuploidy, featuring loss of TP53, LOH of BRCA1 or BRCA2, and multiple breast cancer-associated CNAs in addition to one or more of the common CNAs in 1q, 10q or 16q. Notably, cells with intermediate levels of CNAs were not detected, arguing against a stepwise gradual accumulation of CNAs. Overall, our findings demonstrate that chromosomal alterations in normal breast epithelium including from high risk women, link early stage CNAs to clonally expanded tumor genomes .

  Study EGAS00001007716

• Bone marrow breakout lesions act as key sites for tumor-immune cell diversification and exhaustion in multiple myeloma

  Multiple myeloma is a disease characterized by the expansion of cancerous plasma cells in the bone marrow. The bone marrow microenvironment plays a pivotal role in supporting myeloma growth and modulating tumor immunity. As the disease progresses, myeloma cells may become independent of the bone marrow environment, leading to extramedullary disease associated with poor prognosis. However, the early processes associated with bone marrow independence and its implications for disease and immune control remain poorly understood. Here, we employed comprehensive single-cell and spatial mapping to identify the disruption of the cortical bone and subsequent expansion of myeloma cells in breakout lesions as a key event in multiple myeloma pathogenesis and tumor immunity. Breakout lesions harbor an adapted niche and a unique immune microenvironment conferred by sustained immune-tumor interactions. In particular, tumor-reactive T cells with a highly exhausted phenotype, alongside unique NK cell and macrophage subtypes, specifically expanded in breakout lesions. Spatially-resolved genomic, transcriptomic and cellular analyses uncovered extensive intra-lesion heterogeneity, suggesting a divergent co-evolution of genomic variation and locally confined T cell responses in distinct subregions. Ultra-high plex imaging revealed that tumor-reactive T cells, together with other immune-regulatory cells and dendritic cells, co-expand and exhaust upon sustained tumor interactions in locally confined immune islands. In contrast, immune-regulatory macrophage and NK cell subsets operate spatially separated within myeloma cell deserts. Jointly, our analyses uncover bone marrow breakout lesions as a hotspot for tumor-immune cell interactions, diversification and exhaustion, representing a key event in myeloma pathogenesis with significant therapeutic implications.

  Study EGAS50000000304

• Transcriptional Response to Hypoxia in iPSC-Derived Endothelial Cells from a High Altitude Adapted Population

  We established lymphoblastoid cell lines from 10 adult individuals of Tibetan ancestry living in the Chicago area. These lines were reprogrammed to induced pluripotent stem cells (iPSC). We validated genetic ancestry by genome-wide SNP array genotyping and population genetic analyses. After being subjected to standard quality control testing iPSCs were differentiated into endothelial cells. Following differentiation, endothelial cells were isolated via a pull down protocol using beads coated with an antibody for a canonical surface cell marker of the vascular endothelium (CD144). Purified vascular endothelial cells were cultured in parallel in normoxia (20% O2) and hypoxia (1% O2) for 48 hours prior to harvesting and processing for bulk RNA-sequencing. We provide imputed genotype data and bulk RNA-sequencing data in normoxia and in hypoxia for each of the 10 individuals. We also provide multiplexed scRNA-sequencing data for the 10 Tibetan individuals pooled with 10 CHB individuals from the 1000 Genome Project. Additionally, we provide scRNAseq data for the cardiomyocytes derived from the same iPSC lines from the Tibetan individuals pooled together with cardiomyocytes derived from iPSC lines of 10 CHB individuals from the 1000 Genome Project. Lymphoblastoid cell lines from the CHB were reprogrammed to iPSCs and differentiated to cardiomyocytes using the same protocols as used for the Tibetans. The batches of 4 lines were balanced by sex and population. These iPSC-derived cardiomyocytes were cultured in hypoxia (1% O2) for 48 hours and subjected to scRNAseq in pooled batches of 3 or 4 lines. The CHB individuals included in this analysis were: NA18528, NA18531, NA18557, NA18596, NA18606, NA18608, NA18614, NA18619, NA18633, and NA18748.

  Study phs003758

• A Personalized Neoantigen Vaccine Generates Anti-Tumor Immunity in High-Risk Renal Cell Carcinoma

  Personalized cancer vaccines (PCVs) can generate circulating immune responses against predicted neoantigens. However, whether such responses can target cancer driver mutations, lead to immune recognition of a patient's tumor, and result in clinical activity are largely unknown, particularly in patients with lower mutation burden tumors. We conducted a phase I trial testing a neoantigen-targeting PCV in patients with high-risk, fully resected clear cell renal cell carcinoma (RCC; stage III or IV), with or without ipilimumab administered adjacent to the vaccine. At a median follow-up of 40.2 months from surgery, none of the 9 initial participants experienced a recurrence of RCC. No dose-limiting toxicities were observed. All patients generated T cell immune responses against PCV antigens, including RCC driver mutations in VHL, PBRM1, BAP1, KDM5C, and PIK3CA. Following vaccination, there was a durable expansion of peripheral T cell clones, and T cell reactivity against autologous tumor was detected in 7 of the 9 patients. Our results demonstrate that neoantigen-targeting PCVs in high-risk RCC are highly immunogenic, capable of targeting key driver mutations, and induce anti-tumor immunity. These observations, in conjunction with the absence of recurrence in all 9 vaccinated patients, highlights the promise of PCV as effective adjuvant therapy in RCC. We are submitting WES data from whole blood samples and tumor FFPE samples, and RNA-Seq data from tumor FFPE samples. We are also submitting RNA-Seq data from purified tumor samples. scRNA-seq data will be submitted for skin samples (vaccine site biopsies at Weeks 0, 4, and 13), tumor samples (including resected metastatic tumor samples for two patients), and normal kidney samples. Bulk TCR-seq data will be submitted for blood samples.

  Study phs003710

• Discovery and validation of an ancillary genomic test of malignancy for primary melanocytic tumors

  Background:Pathological diagnosis of melanocytic tumors can be difficult and prone to error. More accurate and reliable ancillary investigations than those currently available are required to improve diagnostic accuracy. Goals: (1) to generate a genomic model for distinguishing benign (nevi) from malignant (melanoma) melanocytic skin tumors, developed from a customized gene panel evaluated on a cohort of 250 nevi and melanomas. (2) To evaluate the performance of the genomic model in melanocytic tumors of uncertain diagnosis/behavior. Details: The cohort of 250 melanomas and naevi was split into a discovery set and a validation set. Our genomic model obtained an area under the curve of 95% in the discovery cohort and 96% for the validation cohort. Based on a Youden index threshold, the resulting classifier shows high specificity (>95% in the discovery and validation cohorts), with sensitivities >80% in the discovery and validation cohorts. Evaluation of the genomic model on a cohort of 110 melanocytic tumors identifies key driver mutations in agreement with the likely pathway of origin of the tumors. The probability of melanoma (p=0.02) from the genomic model was significantly associated with pathology-based estimates of malignancy for melanocytic tumors lacking pathway-defining genomic aberrations (conventional tumors). The genomic model identified 8/36 (22%) conventional tumors with unresolved malignancy to have a high probability of melanoma, with 7/8 of these cases showing negative/non-aberrant results with either fluorescent in-situ hybridization, PRAME and/or p16 immunostains. Furthermore, no evidence of malignant behavior based on the genomic model was seen in 4 conventional borderline tumors with more than 1 year follow-up. Conclusion: This genomic model can become a clinically useful ancillary tool that is highly specific for differentiating melanomas from nevi.

  Study EGAS50000000887

• Enhancer profiling identifies epigenetic markers of endocrine resistance and reveals therapeutic options for metastatic castration-resistant prostate cancer patients (LuCaP cell line data)

  Androgen Receptor (AR) signaling inhibitors, including enzalutamide, are treatment options for patients with metastatic castration-resistant prostate cancer (mCRPC), but resistance inevitably develops. Using metastatic samples from a prospective phase II clinical trial, we epigenetically profiled enhancer/promoter activities with H3K27ac chromatin immunoprecipitation followed by sequencing, before and after AR-targeted therapy. We identified a distinct subset of H3K27ac-differentially marked regions that associated with treatment responsiveness. These data were successfully validated in mCRPC patient-derived xenograft models (PDX). In silico analyses revealed HDAC3 as a critical factor that can drive resistance to hormonal interventions, which we validated in vitro. Critically, we identified the pan-HDAC inhibitor, vorinostat to be effective in decreasing tumor cell proliferation both in vitro and in vivo. Moreover, we uncovered evidence for HDAC3 working together with Glucocorticoid Receptor (GR) as a potential mechanism for this therapeutic effect. These findings demonstrate rationale for new therapeutic strategies including HDAC inhibitors to improve patient outcome in advanced stages of mCRPC. Project Goal: We identified H3K27ac regions associated with resistance in patient samples. Extensive analysis of these regions suggested that HDAC3 may play a role in resistance. To test this, we analyzed ChIP-seq and RNA-seq from cell lines generated from PDX LuCaP samples (recapitulating resistant and non-resistant) to validate findings in patient samples. ChIP-seq samples were used to verify the signals identified from specific factors, such as H3K27ac, AR, FOXA1, GR and HDAC3, their association with resistance and possible mechanisms of action. RNA-seq samples were analyzed to determine the effect of AR signaling in these samples in specific genes in the H3K27ac resistance-associated regions.

  Study EGAS50000000917

• Cellular Profiling Identifies Targetable T Cell Phenotypes in Lymphocytic Variant Hypereosinophilic Syndrome

  Lymphocytic variant hypereosinophilic syndrome (LHES) is driven by aberrant Th2&#8209;skewed lymphocyte clones that promote tissue eosinophilia. We performed multimodal single&#8209;cell RNA&#8209;sequencing (scRNA&#8209;seq), cell&#8209;surface protein (ADT) profiling, and paired &#945;&#946;&#8209;TCR sequencing on peripheral&#8209;blood mononuclear cells (PBMCs) from two clinically characterized LHES patients and one idiopathic HES control. One LHES patient was sequenced at two timepoints, pre-JAK inhibitor treatment and post-JAK inhibitor treatment. Single&#8209;cell regulatory&#8209;network inference (SCENIC), CellChat ligand&#8211;receptor modelling, and scRNA variant calling were integrated with flow&#8209;cytometric validation and 5,000&#215; amplicon sequencing. Two clonally expanded CD3&#8209;CD4&#8314; central&#8209;/effector&#8209;memory T&#8209;cell clusters were observed exclusively in LHES. These cells showed Th2&#8209;associated regulons (GATA3, BATF, MAF), high expression of CCR4, CD52, PTGDR2 and IL13, and down&#8209;regulation of CD247 (CD3&#950;), a mechanism experimentally confirmed to abolish surface CD3/TCR. Autocrine IL&#8209;7 signaling, predicted by CellChat, provides an antigen&#8209;independent proliferative cue. Therapeutically actionable targets (CCR4, CD52, IL&#8209;5, IL&#8209;7, TNF) were prioritized; one refractory patient achieved complete clinical and cytologic remission after mogamulizumab (anti&#8209;CCR4) therapy. Longitudinal sequencing revealed clonal evolution with acquisition of a STAT3 p.G618R mutation under JAK&#8209;inhibition, highlighting a resistance mechanism. Raw scRNA&#8209;seq (5&#8242; cDNA, VDJ, ADT) data, and sample&#8209;level metadata will be distributed through dbGaP. These resources provide the first high&#8209;resolution immune atlas of LHES and a framework for rational, precision immunotherapy in eosinophilic disorders.

  Study phs004041

• Convergent evolution towards CD38 biallelic loss is a recurrent mechanism of resistance to anti-CD38 antibodies in multiple myeloma.

  Monoclonal antibodies directed against CD38 are a therapeutic mainstay in multiple myeloma (MM). While they have contributed to improved outcomes, most patients will still experience disease relapse. Despite their ubiquity, little is known about tumor-intrinsic mechanisms of resistance to these drugs. Antigen escape, whereby tumor cells alter expression of cell surface immunotherapy targets, has been proposed as a potential avenue, and has been seen with other immunotherapies in MM. Yet, it is unknown whether MM cells can develop permanent resistance to anti-CD38 antibodies by acquiring genomic events leading to biallelic disruption of the CD38 gene locus. Here, by using whole genome and whole exome sequencing data from 701 newly diagnosed patients, 72 patients at relapse with naivety to anti-CD38 antibodies, and 52 patients collected at relapse following anti-CD38 antibodies, we report biallelic inactivation of CD38 as a recurrent mechanism of resistance in 6% of patients, often via convergent evolution where distinct subclones are advantaged towards the same selective endpoint within the same patient. Monoallelic loss of the CD38 locus was seen in another 14% post-CD38-directed therapy. Three of the cases with biallelic events on CD38 presented with missense mutations and retained CD38 expression. Functional work demonstrated that two of these mutations decreased binding affinity and antibody-dependent cellular cytotoxicity of commercial antibodies Daratumumab and Isatuximab. However, the third mutation (R140G) provided relative resistance only towards one commercial product. Given that patients are often rechallenged with CD38 antibodies following disease progression, these data support the need for next generation sequencing to guide treatment selection.

  Study EGAS50000000709

• IDH-wildtype untreated human glioblastoma samples (GB-UK cohort), published in Noorani & Haughey et al 2025.

  Extrachromosomal DNA (ecDNA) oncogene amplification is associated with treatment resistance and shorter survival in cancer. Currently, the spatial dynamics of ecDNA, and their evolutionary impact, are poorly understood. Here, we investigate the spatiotemporal evolution of ecDNA by integrating computational modeling with samples from 94 treatment-naive human IDH wild-type glioblastoma patients. Random ecDNA segregation combined with ecDNA-conferred fitness advantages induce distinct spatial ecDNA copy number patterns which depend on ecDNA oncogenic makeup. EGFR-ecDNAs often reach high copy numbers and confer strong fitness advantages. In contrast, PDGFRA-ecDNAs reach lower copy numbers and confer weaker fitness advantages. EGFR structural variants occur exclusively on ecDNA, arise from and are intermixed with wild-type EGFR-ecDNAs. Computational modeling suggests wild-type and variant EGFR-ecDNAs often accumulate before clonal expansion. Using a genetically engineered mouse model, we confirm that the accumulation of oncogenic ecDNA prior to clonal expansion under strong positive selection is observable both in vivo and in vitro in mouse neural stem cells. Intermixed wild-type and variant ecDNAs remain prevalent at diagnosis, as less fit wild-type ecDNAs hitchhike on strongly selected variant EGFR-ecDNAs. A dN/dS analysis in the Genomics England cohort further suggests positive selection of variants on EGFR-ecDNAs but not PDGFRA-ecDNAs in glioblastoma. Overall, our results suggest ecDNA oncogenic makeup determines unique evolutionary trajectories in glioblastoma. New concepts like ecDNA tumor clonality, ecDNA heteroplasmy and multi-species ecDNAs require a refined evolutionary interpretation of genomic data in a large subset of glioblastoma patients and likely other solid cancers.

  Study EGAS00001008126

• RNAseq: Acquired non-permissive bone marrow microenvironment impairs hematopoietic stem cell proliferation and maintenance and B-cell development post-HSCT

  Defects in B-cell reconstitution upon hematopoietic stem-cell (HSC) transplantation (HSCT) are a common observation, yet the mechanism remains unexplained. The bone marrow (BM) stroma, including mesenchymal stromal cells (MSCs), guides HSC maintenance and B-lymphopoiesis by secreting crucial cytokines. We report acquired, permanent, selective and complete B-cell deficiency in the context of full donor-chimerism in a patient with X-linked lymphoproliferative disease and aimed to identify the contribution of the BM-microenvironment in disrupted B-cell reconstitution post-HSCT. We studied longitudinal BM samples from the patient and his identical twin, both of whom underwent HSCT with the same donor with opposite outcomes in B-cell reconstitution. In the index patient BM, we observed progressive loss of proliferation of HSCs and a selective block at the pre-BI cell stage. In vitro modeling studies showed limited survival of patient-HSCs and a relative accumulation of pre-B cells. Patient-derived MSCs failed to support survival and proliferation of HSCs and B-cell development of healthy-HSCs which was correlated with reduced CXCL12 levels. Using bulk RNA-sequencing of MSCs and in vitro functional studies, we showed global changes in the patients’ MSCs and a progressive loss of CXCL12 expression. Indeed, survival of patients HSCs improved supplementing in vitro development culture with CXCL12, suggesting a contribution of defective CXCL12 signaling to the phenotype. In summary, our data show that an acquired defect in the BM-stromal microenvironment and exhaustion of HSCs and committed progenitors may cause a permanent non-permissive state for normal B-cell development.

  Study EGAS50000001437

• ATAC-seq: Acquired non-permissive bone marrow microenvironment impairs hematopoietic stem cell proliferation and maintenance and B-cell development post-HSCT

  Defects in B-cell reconstitution upon hematopoietic stem-cell (HSC) transplantation (HSCT) are a common observation, yet the mechanism remains unexplained. The bone marrow (BM) stroma, including mesenchymal stromal cells (MSCs), guides HSC maintenance and B-lymphopoiesis by secreting crucial cytokines. We report acquired, permanent, selective and complete B-cell deficiency in the context of full donor-chimerism in a patient with X-linked lymphoproliferative disease and aimed to identify the contribution of the BM-microenvironment in disrupted B-cell reconstitution post-HSCT. We studied longitudinal BM samples from the patient and his identical twin, both of whom underwent HSCT with the same donor with opposite outcomes in B-cell reconstitution. In the index patient BM, we observed progressive loss of proliferation of HSCs and a selective block at the pre-BI cell stage. In vitro modeling studies showed limited survival of patient-HSCs and a relative accumulation of pre-B cells. Patient-derived MSCs failed to support survival and proliferation of HSCs and B-cell development of healthy-HSCs which was correlated with reduced CXCL12 levels. Using bulk RNA-sequencing of MSCs and in vitro functional studies, we showed global changes in the patients’ MSCs and a progressive loss of CXCL12 expression. Indeed, survival of patients HSCs improved supplementing in vitro development culture with CXCL12, suggesting a contribution of defective CXCL12 signaling to the phenotype. In summary, our data show that an acquired defect in the BM-stromal microenvironment and exhaustion of HSCs and committed progenitors may cause a permanent non-permissive state for normal B-cell development.

  Study EGAS50000001438

• Development and validation of multivariate predictors of primary endocrine resistance to tamoxifen and aromatase inhibitors in luminal breast cancer reveal drug-specific differences

  Background: Endocrine therapy is highly effective in blocking the estrogen receptor pathway in HR+/HER2- early breast cancer (EBC). However, up to 40% of patients experience relapse during or after adjuvant endocrine therapy. Here, we investigate molecular mechanisms associated with primary resistance to endocrine therapy and develop predictive models. Results: TP53 mutations were prominently associated with primary resistance to both tamoxifen (TAM) and aromatase inhibitors (AI), with AI non-responders exhibiting resistance in up to 32% of cases. Additionally, we identified distinct DNA methylation patterns in TAM and AI non-responders, with TAM non-responders showing global DNA methylation loss, associated with KRAS signaling, apical junctions and epithelial-mesenchymal transition (EMT). Conversely, we observed methylation gain in AI non-responders affecting developmental transcription factors, hypoxia and estrogen signaling. TAM or AI resistance was associated with increased methylation-inferred proportions of immune cells and decreased proportions of endothelial cells. Based on these findings and patient age, we developed the Predictive Endocrine ResistanCe Index (PERCI). PERCI stratified NR and R cases in both treatment groups and cohorts with high accuracy (ROC AUC TAM discovery 93.9%, validation 83%; AI discovery 98.6%, validation 76.9%). A simplified PERCI efficiently predicted progression-free survival in the TCGA-BRCA sub-cohort (Kaplan-Meier log-rank p-value = 0.03 between low and high PERCI groups). Conclusions: We identified genomic and epigenomic features associated with primary resistance to TMA and AI. By combining information on genomic alterations, patient age, differential methylation and tumor microenvironment (TME) composition, we developed PERCI TAM and PERCI AI as novel predictors of primary resistance, with potential additional prognostic value. Applying PERCI in a clinical setting may allow patient-specific drug selection to overcome resistance.

  Study EGAS50000001102

• CCR8-targeted specific depletion of clonally expanded Treg cells in tumor tissues evokes potent tumor immunity with long-lasting memory

  Foxp3-expressing CD25+CD4+ regulatory T cells (Tregs) are abundantly infiltrating into tumor tissues, hindering effective anti-tumor immunity. Here, assuming that such tumor Tregs were clonally expanding by recognizing tumor-associated antigens, we attempted to characterize them by their gene expression profiling and T-cell receptor (TCR) clonotyping at the single-cell level. We found that they predominantly expressed the chemokine receptor CCR8: ~40% of tumor Tregs, less than 10% of Tregs in other tissues, and few tumor-infiltrating conventional T cells expressed CCR8 in mice and humans. CCR8+ tumor Tregs were phenotypically highly differentiated and functionally stable. One injection of anti-CCR8 mAb with high antibody-dependent cellular cytotoxic activity was indeed able to specifically deplete clonally expanded tumor Tregs and thereby eradicate established tumors in mice. The depletion predominantly expanded tumor-infiltrating CD8+ effector T cells, which were different in activated or exhausted patterns from those seen after immune checkpoint blockade by anti-PD-1 mAb. The treatment also evoked potent secondary immune responses against the same tumor cell line inoculated several months after tumor eradication, indicating that the anti-CCR8-induced memory-type tumor-reactive effector T cells were, once induced, resistant to recovered Treg suppression. In contrast to such evocation of potent tumor immunity, anti-CCR8 treatment scarcely elicited discernable autoimmunity in mice, contrasting with a similar treatment for systemic whole Treg depletion, which induced histologically and serologically evident autoimmunity. Thus, temporary and specific depletion of clonally expanded Tregs in tumor tissues by cell-depleting anti-CCR8 mAb is sufficient to evoke potent tumor immunity with long-lasting memory and without serious autoimmunity.

  Study JGAS000312

• mutation analysys of Gorlin syndrome

  Gorlin syndrome is a genetic disorder of autosomal dominant inheritance that predisposes the affected individual to a variety of disorders that are attributed largely to heterozygous germline patched1 (PTCH1) mutations. PTCH1 is a hedgehog (Hh) receptor as well as a repressor, mutation of which leads to constitutive activation of Hh pathway. Hh pathway encompasses a wide variety of cellular signaling cascades, which involve several molecules; however, no associated genotype?phenotype correlations have been reported. Recently, mutations in Suppressor of fused homolog (SUFU) or PTCH2 were reported in patients with Gorlin syndrome. These facts suggest that multi-layered mutations in Hh pathway may contribute to the development of Gorlin syndrome. We demonstrated multiple mutations of Hh-related genes in addition to PTCH1, which possibly act in an additive or multiplicative manner and lead to Gorlin syndrome. High-throughput sequencing was performed to analyze exome sequences in four unrelated Gorlin syndrome patient genomes. Mutations in PTCH1 gene were detected in all four patients. Specific nucleotide variations or frameshift variations of PTCH1 were identified along with the inferred amino acid changes in all patients. We further filtered 84 different genes which are closely related to Hh signaling. Fifty three of these had enough coverage of over ?30. The sequencing results were filtered and compared to reduce the number of sequence variants identified in each of the affected individuals. We discovered three genes, PTCH2, BOC, and WNT9b, with mutations with a predicted functional impact assessed by MutationTaster2 or PolyPhen-2 (Polymorphism Phenotyping v2) analysis. It is noticeable that PTCH2 and BOC are Hh receptor molecules. No significant mutations were observed in SUFU. Multi-layered mutations in Hh pathway may change the activation level of the Hh signals, which may explain the wide phenotypic variability of Gorlin syndrome.

  Study JGAS000099

• Immunohistochemical and molecular pathological search in gastrointestinal tumors

  Purpose: RSPO fusions are recently identified recurrent genetic alterations in colorectal cancer (CRC), which lead to WNT pathway activation and are potential therapeutic targets. We aimed to investigate the clinicopathological and molecular features of RSPO fusion-positive CRCs. Experimental Design: We screened 1019 CRCs for RSPO fusions using multiplex reverse transcription polymerase chain reaction. RSPO fusion-positive tumors were characterized for their clinicopathological features and subjected to whole exome sequencing (WES). Results: Our analysis identified 29 CRCs with RSPO fusions (2.8%), consisting of 5 with an EIF3E-RSPO2 fusion and 24 with PTPRK-RSPO3 fusions. The patients were 17 females and 12 males with the median age of 64 years. Thirteen tumors were located in the right colon, 8 were in the left colon, and 8 were in the rectum. Histologically, most of the tumors (24/29, 83%) were classified as adenocarcinomas, not otherwise specified. However, about half of the tumors (13/29, 45%) had a focal or extensive mucinous component, which is significantly more frequent compared with RSPO fusion-negative tumors (13%; p = 8.1 �� 10-7). Four tumors (14%) were mismatch repair-deficient, three of which were associated with MLH1 promoter methylation. WES identified KRAS, BRAF, and NRAS mutations in a total of 27 cases (93%). In contrast, pathogenic mutations in major WNT pathway genes other than RSPO fusions were virtually absent. The presence of an RSPO fusion tended to show shorter overall and recurrence-free survivals in stage II and III patients although the differences were not significant. Conclusions: RSPO fusion-positive CRCs constitute a rare subgroup of CRCs with several characteristic clinicopathological and genetic features. The mutual exclusivity with other WNT pathway gene alterations supports RSPO fusions as valid therapeutic targets in CRC.

  Study JGAS000538

• Genome and gene analysis of gastrointestinal cancer and elucidation of its clinicopathological significance

  Cancer develops through the accumulation of genetic and epigenetic aberrations. To identify sequential molecular alterations that occur during the development of hepatocellular carcinoma (HCC), we molecularly compared 52 early and 108 overt HCC samples by genome sequencing. Gene mutations in the p53/RB1 pathway, WNT pathway, MLLs, SWI/SNF complexes, and AKT/PI3K pathway were common in HCC. In the early phase of all entities, TERT was the most frequently upregulated gene owing to diverse mechanisms. Despite frequent somatic mutations in driver genes, including CTNNB1 and TP53, early HCC was a separate molecular entity from overt HCC as each had a distinct expression profile. Notably, WNT target genes were not activated in early HCC, regardless of CTNNB1 mutation status, because ?-catenin did not translocate into the nucleus due to the E-Cadherin/?-catenin complex at the membrane. Conversely, WNT targets were definitively upregulated in overt HCC with CTNNB1 mutation associated with downregulation of CDH1 and hypomethylation of CpG islands in the target genes. Similarly, cell cycle genes downstream of the p53/RB pathway were upregulated only in overt HCC with TP53 and/or RB1 gene mutations associated with chromosomal deletion of 4q and/or 16q. Differing from the expression profile, HCC was epigenetically distinguished into four subclasses: normal-like methylation, global-hypomethylation (favorable prognosis), stem-like methylation (poor prognosis), and CpG island methylation. These methylation statuses were globally maintained through HCC progression. Taken together, additional molecular events besides driver gene mutations cooperatively contribute to transcriptional activation of downstream targets as HCC progresses, according to the methylation status of the associated genes.

  Study JGAS000233

• Paediatric Hepatic International Tumour Trial (JPLT2: PHITT)

  Hepatoblastoma (HBL) is the most common primary liver tumor in children. Despite their diversity in histopathological phenotype and clinical outcome, the underlying molecular aberrations have not been fully understood. Here we describe a comprehensive landscape of genomic, epigenomic, and transcriptomic heterogeneity among the 163 cases of HBLs and transitional liver cell tumors (TLCTs) collected in the multicenter clinical trial JPLT-2. The total numbers of mutational events in HBLs are very few (0.52 SNVs/Mb on exonic regions) but show a positive correlation with the age of onset. While dysregulation of Wnt-pathway related genes was observed (somatic deletion or mutations of CTNNB1 and germline or somatic mutations of APC) in 83.6% of the HBL cases, TERT promoter mutations are mainly in the elderly HBLs (> 8 years old), meeting with the proposed criteria of TLCT. These TLCT cases are clustered in the hypermethylated epigenotype in promoter and enhancer regions. Copy number analyses also revealed recurrent 11p15.5 UPD/LOHs including the IGF2/H19 region (31.6%), focal deletions in 4q35 (15.8%) and focal amplifications in 2q24.3 (7.0%) predominantly in the child and infant HBLs. Motif enrichment analysis of the specific hypomethylated enhancers in these HBL subtypes suggests specific binding of ASCL2, the regulatory transcription factor for definitive endoderm and digestive systems. Transcriptome and histochemical analysis revealed that ASCL2 is upregulated not in HCC and TLCT but in child and infant HBL, suggesting their ?cell of origin? derived from the premature hepatoblast with high proliferative potential such as intestinal epithelial cells. Systematic molecular profiling of HBL will provide a fundamental understanding of hepatoblast carcinogenesis and clues for risk stratification.

  Study JGAS000236

• Biallelic variants in the non-protein coding minor spliceosome components RNU6ATAC and RNU6ATAC cause syndromic monogenic autoimmune diabetes

  Non-protein coding genes are emerging as critical contributors to the aetiology of rare diseases, providing key insights into human biology and uncovering novel disease mechanisms. We identified 7 individuals from 4 families with early-onset diabetes (diagnosed <5 years) and immune dysregulatory features caused by biallelic variants in RNU6ATAC. RNU6ATAC encodes a small nuclear RNA (snRNA) acting as a catalytic component of the minor spliceosome, a protein-RNA complex mediating splicing of ~700 genes containing U12/minor-type introns. Variant screening of the other 64 minor spliceosome genes in 276 infants with diabetes identified 12 unrelated individuals with biallelic disease-causing variants in RNU4ATAC. Biallelic pathogenic RNU4ATAC variants are known to cause a variable spectrum of clinical features, which until now did not include diabetes. Clinically, 12/19 RNU6ATAC/RNU4ATAC patients had additional immune dysregulatory features, and 50% of patients tested were islet-autoantibody positive, strongly supporting an autoimmune aetiology for their diabetes. RNA-seq in 3 individuals with biallelic RNU6ATAC variants showed a pattern of intron retention in U12 intron-containing genes similar to that seen in RNU4ATAC-individuals (n=3), supporting a shared disease mechanism. Analysis of patients’ transcriptomic, methylation and immune data revealed impaired B cell development and maturation. We conclude that biallelic RNU6ATAC variants cause a novel syndrome of early-onset autoimmune diabetes and immune dysregulation. We further show that infancy-onset diabetes is a previously unrecognised feature of RNU4ATACopathy. Our work highlights the important role of two snRNAs critical to minor spliceosome function in immune system regulation, providing novel insights into the pathogenesis of autoimmune diabetes.

  Study EGAS50000001565

• Confronting historical legacies of biological anthropology in South Africa—Restitution, redress and community-centered science: The Sutherland Nine

  We describe a process of restitution of nine unethically acquired human skeletons to their families, together with attempts at redress. Between 1925–1927 C.E., the skeletonised remains of nine San or Khoekhoe people, eight of them known-in-life, were removed from their graves on the farm Kruisrivier, near Sutherland in the Northern Cape Province of South Africa. They were donated to the Anatomy Department at the University of Cape Town. This was done without the knowledge or permission of their families. The donor was a medical student who removed the remains from the labourers’ cemetery on his family farm. Nearly 100 years later, the remains are being returned to their community, accompanied by a range of community-driven interdisciplinary historical, archaeological and analytical (osteobiographic, craniofacial, ancient DNA, stable isotope) studies to document, as far as possible, their lives and deaths. The restitution process began by contacting families living in the same area with the same surnames as the deceased. The restitution and redress process prioritises the descendant families’ memories, wishes and desire to understand the situation, and learn more about their ancestors. The descendant families have described the process as helping them to reconnect with their ancestors. A richer appreciation of their ancestors’ lives, gained in part from scientific analyses, culminating with reburial, is hoped to aid the descendant families and wider community in [re-]connecting with their heritage and culture, and contribute to restorative justice, reconciliation and healing while confronting a traumatic historical moment. While these nine individuals were exhumed as specimens, they will be reburied as people.

  Study EGAS50000000971

• Repotrectinib in NTRK fusion–positive advanced solid tumors: a phase 1/2 trial

  Early-generation TRK tyrosine kinase inhibitors (TKIs) approved for treating NTRK fusion–positive (NTRK+) solid tumors provide clinical benefit; however, resistance emerges. Repotrectinib is a next-generation ROS1/TRK TKI with a compact macrocyclic structure designed to improve durability of response. TRIDENT-1 is a registrational phase 1/2 trial assessing repotrectinib, a next-generation ROS1/TRK TKI, in adults with advanced solid tumors, including NTRK+ disease. The primary endpoint was confirmed objective response; secondary endpoints included duration of response (DOR), progression-free survival (PFS), overall survival and safety. Median follow-up ranged between 21.3 months and 25.7 months. In the TKI-naive cohort (n = 51; 95% confidence interval (CI)), the response rate was 59% (44–72); the median DOR was not estimable (NE); and the median PFS was 30.3 months (9.0–NE). In the TKI-pretreated cohort (n = 69; 95% CI), the response rate was 48% (36–60); the median DOR was 9.8 months (7.4–13.0); and the median PFS was 7.4 months (3.9–9.7). Of 30 TKI-pretreated patients with NTRK solvent front mutations, 16 had a response (53%; 95% CI: 34–72). Intracranial responses were observed in two of three TKI-naive patients and in four of six TKI-pretreated patients with measurable intracranial disease at baseline. Among all treated patients (n = 565), the most common any-grade treatment-related adverse event (TRAE) was dizziness (57%); most TRAEs were low grade; and 4% discontinued repotrectinib due to a TRAE. Overall, repotrectinib demonstrated durable systemic and intracranial responses with generally low-grade adverse events in patients with NTRK+ solid tumors, including those with previous TRK TKI treatment and solvent front mutations. These results support the use repotrectinib to treat patients with NTRK+ solid tumors. ClinicalTrials.gov identifier: NCT03093116 .

  Study EGAS50000001572

• Recurrent somatic DICER1 mutations in non-epithelial ovarian tumors

  Background: Germline truncating mutations in DICER1, an RNase III-type endoribonuclease essential for processing of microRNAs, are seen in families with the pleuropulmonary blastoma-family tumor and dysplasia syndrome; these individuals have a propensity to develop non-epithelial ovarian tumors (i.e. sex-cord stromal and germ cell tumors). Methods: We sequenced the whole transcriptomes or exomes of 14 non-epithelial ovarian tumors and noted closely clustered mutations in the RNase IIIb domain of DICER1 in four cases. Using Sanger sequencing, we then analyzed additional ovarian tumors for mutations in this region. The impact of the mutations on the enzymatic activity of Dicer1 was determined using in vitro RNA cleavage assays. Results: RNase IIIb domain DICER1 mutations were found in 30/102 non-epithelial ovarian tumors (29%), predominantly Sertoli-Leydig cell tumors (26/43, 60%), including 4/4 from individuals with germline DICER1 mutations. The mutations were restricted to metal binding sites within the RNase IIIb catalytic centres critical for miRNA interaction and cleavage, and were somatic in all 14 cases where germline DNA was available for testing. Of 266 epithelial ovarian and endometrial cancers tested, one case, an ovarian carcinosarcoma, had a hotspot mutation. In vitro, the mutant DICER proteins showed reduced RNase IIIb activity but retained RNase IIIa activity. Conclusions: Somatic missense mutations of the RNase IIIb domain of DICER1 are common in non-epithelial ovarian tumors. Genetic, pathologic, and functional evidence suggest that these mutations do not obliterate DICER1 function, rather, they alter it in specific cell types, a novel mechanism through which perturbed microRNA expression is oncogenic.

  Study EGAS00001000135

• Driver mutations in histone H3.3 and chromatin remodelling genes in paediatric glioblastoma

  Glioblastoma multiforme (GBM) is a lethal brain tumour in adults and children. However, DNA copy number and gene expression signatures indicate differences between adult and paediatric cases. To explore the genetic events underlying this distinction, we sequenced the exomes of 48 paediatric GBM samples. Somatic mutations in the H3.3-ATRX-DAXX chromatin remodelling pathway were identified in 44% of tumours (21/48). Recurrent mutations in H3F3A, which encodes the replication-independent histone 3 variant H3.3, were observed in 31% of tumours, and led to amino acid substitutions at two critical positions within the histone tail (K27M, G34R/G34V) involved in key regulatory post-translational modifications. Mutations in ATRX (alpha-thalassaemia/mental retardation syndrome X-linked) and DAXX (death-domain associated protein), encoding two subunits of a chromatin remodelling complex required for H3.3 incorporation at pericentric heterochromatin and telomeres, were identified in 31% of samples overall, and in 100% of tumours harbouring a G34R or G34V H3.3 mutation. Somatic TP53 mutations were identified in 54% of all cases, and in 86% of samples with H3F3A and/or ATRX mutations. Screening of a large cohort of gliomas of various grades and histologies (n=784) showed H3F3A mutations to be specific to GBM and highly prevalent in children and young adults. Furthermore, the presence of H3F3A/ATRX-DAXX/TP53 mutations was strongly associated with alternative lengthening of telomeres and specific gene expression profiles. This is, to our knowledge, the first report to highlight recurrent mutations in a regulatory histone in humans, and our data suggest that defects of the chromatin architecture underlie paediatric and young adult GBM pathogenesis.

  Study EGAS00001000226

• Molecular_diagnosis_of_albinism

  Albinism is genetically heterogeneous rare genetic condition affecting 1:17000 in the Western world (but more frequent in Africa) whose main feature is a profound visual impairment, characterised by foveal hypoplasia, abnormal chiasmatic connections, nystagmus and photofobia. All these features result in severly altered visual acuity (&lt;0,1), absent depth perception and poor night vision. People with albinism are primarily visually handicapped. In addition, for some types of albinism, the visual phenotype can be presented with partial or total hypopigmentation, hence resulting in a secondary phenotype which can lead to skin cancer if skin is not adequately protected. Recently a new syndrome has been described, FHONDA, with the same visual abnormalities of albinism but without pigment alteration. The traditional classification differentiates Oculoculatenous albinism (OCA), where hypopigmentation involves hair, skin and eyes versus Ocular Albinism (OA), where hypopigmentation only affects the eyes. These are non-sydrimic types of albinism. Some syndromic forms (Hermansky-Pudlak=HPS, Chediak-Higashi=CHS) affect cells beyond pigment cells, present in the lungs, immune system, platelets and intestines, resulting in more severe phenotypes that can be fatal. Mutations in at least 19 genes are assocaited with the corresponding types of albinism. Most hospitals will only diagnose the most frequent cases using traditional Sanger, MLPA approaches. Some will use CGH arrays. We aim to diagnose all cases of albinism through the Albinochip proposal, which combines a Sequenom first step of known mutations combined with subsequent NGS approaches. In some cases we fail to find a second mutation, these are good candidates for further full exome 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 EGAS00001002068

• Resolving the Full Spectrum of Human Genome Variation using Linked-Reads

  Large-scale population based analyses coupled with advances in technology have demonstrated that the human genome is more diverse than originally thought. Standard short-read approaches, used primarily due to accuracy, throughput and costs, fail to give a complete picture of a genome. They struggle to identify large, balanced structural events, cannot access repetitive regions of the genome and fail to resolve the human genome into its two haplotypes. Here we describe an approach that retains long range information while harnessing the power of short reads. Starting from only 1ng of DNA, we produce barcoded short read libraries. The use of novel informatic approaches allows for the barcoded short reads to be associated with the long molecules of origin producing a novel datatype known as 'Linked-Reads'. This approach allows for simultaneous detection of small and large variants from a single Linked-Read library. We have previously demonstrated the utility of whole genome Linked-Reads (lrWGS) for performing diploid, de novo assembly of individual genomes. In this manuscript, we show the utility of reference based analysis using a single Linked-Read library for full spectrum genome analysis. We demonstrate the ability of Linked-Reads to reconstruct megabase scale haplotypes and to recover parts of the genome that are typically inaccessible to short reads, including phenotypically important genes such as STRC, SMN1 and SMN2. We demonstrate the ability of both lrWGS and Linked-Read Whole Exome Sequencing (lrWES) to identify complex structural variations, including balanced events, single exon deletions, and single exon duplications. The data presented here show that Linked-Reads provide a scalable approach for comprehensive genome analysis that is not possible using short reads alone.

  Study EGAS00001003121

• A genome scan for genes underlying adult body size differences between Central African hunter-gatherers and farmers.

  The evolutionary and biological bases of the Central African pygmy phenotype, a characteristic of rainforest huntergatherers defined by reduced body size compared with neighboring farmers, remain largely unknown. Here, we perform a joint investigation in Central African hunter-gatherers and farmers of adult standing height, sitting height, leg length, and body mass index (BMI), considering 358 hunter-gatherers and 169 farmers with genotypes for 153,798 SNPs. In addition to reduced standing heights, hunter-gatherers have shorter sitting heights and leg lengths and higher sitting/standing height ratios than farmers and lower BMI for males. Standing height, sitting height, and leg length are strongly correlated with inferred levels of farmer genetic ancestry, whereas BMI is only weakly correlated, perhaps reflecting greater contributions of non-genetic factors to body weight than to height. Single- and multi-marker association tests identify one region and eight genes associated with hunter-gatherer/farmer status, and 24 genes associated with the height-related traits. Many of these genes have putative functions consistent with roles in determining their associated traits and the pygmy phenotype, and they include three associated with standing height in non-Africans (PRKG1, DSCAM, MAGI2). We find evidence that European height-associated SNPs or variants in linkage disequilibrium with them contribute to standing- and sitting-height determination in Central Africans, but not to the differential status of hunter-gatherers and farmers. These findings provide new insights into the biological basis of the pygmy phenotype, and they highlight the potential of cross-population studies for exploring the genetic basis of phenotypes that vary naturally across populations. Pemberton et al. Human Genetics 2018 https://doi.org/10.1007/s00439-018-1902-3

  Study EGAS00001002975

• Indonesian Genome Diversity Project

  Indonesia, the world’s fourth largest country by population, comprises an archipelago of about 900 permanently inhabited islands in tropical Asia, and hosts an astonishing array of human diversity that remains largely underrepresented in modern biological surveys. The region is unique for its key role in both the early and recent evolution of Homo outside of Africa. It has one of the first traces of anatomically modern humans in Eurasia, and was the focus of extensive interaction between archaic hominins and humans – both in the likely co-existence of H. floresiensis with modern humans, and in the introgression of Denisovans into the ancestors of its present inhabitants. More recently, Indonesia was a center of the spread of Neolithic culture by Austronesian speakers. Advancing maritime technologies allowed farming populations to treat this region as a springboard to reach oceanic islands as remote as Madagascar, Hawaii, Easter Island and New Zealand. Here we present complete human genomes from healthy individuals from multiple key islands in Indonesia and Papua. The genomes were sequenced to high depth (30x). The purpose of this dataset is to address different questions of human evolution in Island Southeast Asia and beyond, including early settlement of the area by anatomically modern humans, population structure, genetic admixture and adaptation, and contacts with archaic hominins. We employed a sampling strategy to capture a broad range of diversity across the region incorporating an average of 12 samples from each of 14 islands spanning Indonesia and Papua, from Sumatra in the west to New Britain in the east. Genomes from this dataset form part of an ongoing project to describe and understand human evolutionary history in Island Southeast Asia and Eurasia in general.

  Study EGAS00001003054

• Genetic characterization of B-cell prolymphocytic leukemia: a hierarchical prognostic model involving MYC and TP53 abnormalities - WXS

  B-PLL is a rare disease, whose molecular pathology is largely unknown. We report here the cytogenetic and molecular findings in a large series of 34 B-PLL. Karyotype (K) was complex (≥3 abnormalities) in 73%, and highly complex (HCK≥5) in 45%. The most frequent chromosomal aberrations were: translocation targeting the MYC gene [t(MYC)] (62%), 17p deletion including TP53 gene (38%), trisomy 18/18q (30%), 13q14 deletion (29%), trisomy 3 (24%), trisomy 12 (24%) and 8p deletion (23%). Whole-Exome Sequencing performed in 16 patients revealed recurrent mutations in TP53 (6/16, 38%), MYD88 (n=4), BCOR (n=4), MYC (n=3), SF3B1 (n=3), FAT1 (n=3), SETD2 (n=2), CHD2 (n=2), CXCR4 (n=2) and BCLAF1 (n=2). The main group of patients (21/34, 62%) had a t(MYC) associated with a higher percentage of prolymphocytes (p=0.03), CD38 expression (p&lt;0.001), lower K complexity (p=0.0004), mutations in MYC and in genes involved in RNA metabolism and chromatin remodeling. Principal component analysis of gene expression data showed that patients with t(MYC) clustered together. A second group with MYC gain (5/34, 15%), was associated with HCK≥5 (p=0,01) and trisomy 3 (p=0,008). Altogether, 26/34 patients (76%) had a MYC activation, translocation or gain, that were mutually exclusive. We identified 3 distinct cytogenetic prognostic groups (p=0.0006): lower risk: absence of MYC activation (median not reached); intermediate risk: MYC activation without del17p (125 months); high risk: MYC activation + del17p (11 months). Our results show that cytogenetic analysis is a useful diagnostic tool in B-PLL that improves prognostic stratification.

  Study EGAS00001003275

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

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

  Study EGAS00001003481

• An integrative model of pathway convergence in genetically heterogeneous blast crisis chronic myeloid leukemia (CML)

  Targeted therapies against the BCR-ABL1 kinase have revolutionized treatment of chronic phase (CP) chronic myeloid leukemia (CML). In contrast, management of blast crisis (BC) CML remains challenging because BC cells acquire complex molecular alterations that confer stemness features to progenitor populations and resistance to BCR-ABL1 tyrosine kinase inhibitors. Comprehensive models of BC transformation have proved elusive because of the rarity and genetic heterogeneity of BC, but are important for developing biomarkers predicting BC progression and effective therapies. To better understand BC, we performed an integrated multiomics analysis of 74 CP and BC samples using whole-genome and exome sequencing, transcriptome and methylome profiling, and chromatin immunoprecipitation followed by high-throughput sequencing. Employing pathway-based analysis, we found the BC genome was significantly enriched for mutations affecting components of the polycomb repressive complex (PRC) pathway. While transcriptomically, BC progenitors were enriched and depleted for PRC1- and PRC2-related gene sets respectively. By integrating our data sets, we determined that BC progenitors undergo PRC-driven epigenetic reprogramming toward a convergent transcriptomic state. Specifically, PRC2 directs BC DNA hypermethylation, which in turn silences key genes involved in myeloid differentiation and tumor suppressor function via so-called epigenetic switching, whereas PRC1 represses an overlapping and distinct set of genes, including novel BC tumor suppressors. On the basis of these observations, we developed an integrated model of BC that facilitated the identification of combinatorial therapies capable of reversing BC reprogramming (decitabine+PRC1 inhibitors), novel PRC-silenced tumor suppressor genes (NR4A2), and gene expression signatures predictive of disease progression and drug resistance in CP.

  Study EGAS00001001751

• Genomic and epigenomic characterization of juvenile myelomonocytic leukemia (JMML)

  Juvenile myelomonocytic leukemia (JMML) is an aggressive myeloproliferative disorder of early childhood. While some cases show spontaneous remission, allogeneic hematopoietic stem cell transplantation (HSCT) remains the only curative treatment option for the majority of patients, however, the 5-year event-free survival reaches only about 50%. Hyperactive RAS signaling is assumed to be the main driving event in JMML. It is caused by genetic alterations in CBL, KRAS, NF1, NRAS, or PTPN11 in about 90% of patients. So far, there is no clear understanding of how RAS pathway mutations relate to the heterogeneous disease biology and variable clinical outcome seen in JMML patients. As a consequence, established clinical and genetic markers fail to fully represent the observed disease heterogeneity. We hypothesized that DNA methylation profiling, either alone or in combination with genetic alterations, might provide a molecular basis for disease classification. Genome wide DNA methylation analysis using the HumanMethylation450 Bead Chip array was performed in a discovery cohort of 20 patients (19 JMML, 1 Noonan syndrome). The clinical relevance of our findings was assessed in an unselected sample set consisting of 147 consecutive patients with JMML (n=129) or Noonan syndrome (n=18) associated myeloproliferative disorder (n=18) registered in the EWOG-MDS 1998 &amp; 2006 trials. Data integration was performed in a subset of patients with available exome sequencing (n=50) and expression profiling (n=15) data. This integrated approach identified three JMML subgroups characterized by distinct clinical and biological features. This analysis further suggested a molecular mechanism by which additional genetic events, presumably further activating the RAS-RAF-MEK-ERK pathway, mediate DNA hypermethylation via up-regulation of DNMTs in more aggressive JMML cases.

  Study EGAS00001002511

• The application of RNA sequencing for the diagnosis and genomic classification of pediatric acute lymphoblastic leukemia

  Acute lymphoblastic leukaemia (ALL) is the most common childhood malignancy, and implementation of risk-adapted therapy has been instrumental in the dramatic improvements in clinical outcomes. A key to risk-adapted therapies includes the identification of genomic features of individual tumors, including chromosome number (for hyper- and hypodiploidy) and gene fusions, notably ETV6-RUNX1, TCF3-PBX1 and BCR-ABL1 in B-cell ALL (B-ALL). RNA-sequencing (RNA-seq) of large ALL cohorts has expanded the number of recurrent gene fusions recognized as drivers in ALL, and identification of these new entities will contribute to refining ALL risk stratification. We used RNA-seq on 126 ALL patients from our clinical service to test the utility of including RNA-seq in standard of care diagnostic pipelines to detect gene rearrangements and IKZF1 deletions. RNA-seq identified 86% of rearrangements detected by standard of care diagnostics. KMT2A-rearrangements (MLL), whilst usually identified, were the most commonly missed by RNA-seq as a result of low expression. RNA-seq identified rearrangements that were not detected by standard of care testing in nine patients. These were found in patients who were not classifiable using standard molecular assessment. We developed an approach to detect the commonest IKZF1 deletion from RNA-seq data and validated this using a RQ-PCR assay. We applied an expression classifier to identify Ph-like B-ALL patients. T-ALL proved a rich source of novel gene fusions which have clinical implications or provide insights into disease biology. Our experience shows that RNA-seq can be implemented within an individual clinical service to enhance the current molecular diagnostic risk classification of ALL.

  Study EGAS00001004212

• A Multifactorial Tumor and Immune Cell Profile Determines Response to Immune Checkpoint blockade in Melanoma

  The complex crosstalk between tumor and immune cells during immune checkpoint blockade mandates the development of integrated models that interpret the anti-tumor immune response and predict clinical outcome. We have integrated genome-wide sequence and structural alterations with pre and on-therapy transcriptomic and T cell repertoire features for a cohort of 64 immunotherapy-naïve melanomas treated with anti-PD1 monotherapy or combined anti-PD1 and anti-CTLA4 therapy. While tumor mutation burden (TMB) was associated with improved response to therapy, expressed mutation burden was superior to TMB in predicting outcome. An increased pre-existing T cell density differentiated responding from non-responding tumors independent of therapy. Importantly, T cell repertoire reshaping determined by T cell receptor (TCR) clonotypic regressions and expansions reflected response to therapy such that a more dynamic repertoire predicted improved clinical outcome. Through whole-transcriptome analyses, we discovered differential abundance of B cell subsets in responding tumors, which highlights the importance of the interplay between pre-existing T and B cell immunity in shaping therapeutic response. Tumor temporal trajectories during therapy revealed distinct patterns of molecular response related to expressed mutation elimination or retention that accurately interpreted clinical response. High-dimensional genomic, transcriptomic and immune repertoire data were integrated by both a machine learning and a censored regression approach, resulting in a harmonized multi-modal predictor of response to immune checkpoint blockade. B cell abundance, expressed mutation load and tumor aneuploidy were combined to identify patients at high risk for recurrence, such that high risk patients had a significantly shorter progression-free survival (HR=9.18, 95% CI: 3.14-26.85, p=3.4e-06) especially in the anti-PD1/anti-CTLA4 group (HR=20.34, 95% CI: 4.08-101.32, p=1.68e-06).

  Study EGAS00001004548

• Circulating Tumor DNA in Checkpoint Inhibitor treated Lung Cancer

  BackgroundImmune checkpoint inhibitors (ICI) are increasingly used to treat non-small cell lung cancer (NSCLC) and there is a critical need for novel predictive/prognostic biomarkers. We evaluated whether dynamics of circulating tumor (ctDNA) can identify NSCLC patients who benefit from single agent ICI.MethodsCirculating cell-free DNA (cfDNA) from 177 patients was analyzed prior ICI and at first treatment evaluation using a Commercial 77-gene panel. We performed de novo mutation calling and corrected for variants related clonal hematopoiesis. ctDNA changes were evaluated as predictors of tumor response estimated with PET/diagnostic CT and survival.ResultsMutations were detected in 159 of 174 (91.4%) analyzable patients. Most responders presented decreasing ctDNA levels, while non-responders had mixed molecular responses. Patients with a &gt;50% molecular ctDNA decrease or no detectable ctDNA demonstrated improved PFS (5.9 versus 1.8 months; HR 0.49, 95%CI, 0.35-0.67, P &lt; .001), and OS (18.9 versus 9.9 months; HR 0.41, 95%CI, 0.30-0.57, P &lt; .001) compared to patients not achieving these endpoints. Clearance of mutations (and no detected ctDNA) had a similar effect (PFS HR 0.45, 95%CI, 0.32-0.61, P &lt; .001; OS 0.37,95%CI, 0.27-0.52, P &lt; .001) compared to patients having stable or additional mutations. After adjusting for clinical factors both &gt;50% ctDNA decrease (HR 0.4, 95% CI, 0.3-0.6, P &lt; .001) and clearance of mutations (HR 0.42, 95% CI, 0.3-0.6, P &lt; .001) remained independent predictors for PFS and OS.ConclusionsctDNA testing reveals significant prognostic albeit not predictive information in ICI treated advanced NSCLC. This relevant information can be obtained with a commercial easy-to-implement ctDNA test that will facilitate extensive testing.

  Study EGAS00001004847