14347 results for ""Pennsylvania+German"+ACS+study+'Socius'+open+access"

in 320.19 milliseconds.

• Renal_habitat_WXS

  This project aims at correlating genetic diversity (as well as transcriptomic, metabolic and immune diversity) observed in multiple targeted tumour biopsies and circulating tumour DNA with heterogeneity on multiparametric MRI. To this end, 20 patients with RCC (15 patients with localised/locally advanced disease and 5 with mRCC) will undergo multiparametric MRI including metabolic assessment by 13C-MRSI as well as measurements of tumour perfusion, oxygenation and water diffusion. Based on unique combinations of these imaging parameters, tumour habitats will be defined through cluster mapping. Intra-operative biopsy and immediate freezing of tumour cores from different imaging habitats will be performed. Multiregional ex-vivo tumour biopsies will be spatially registered to the MR data by using 3D printed tumour moulds generated from high-resolution volumetric sequences. These moulds allow slicing tumour specimens in identical planes to the MR acquisition and increase the accuracy of tissue sampling. The biopsies will then undergo molecular profiling such as metabolomic and RNA sequencing (including mRNA immune profiling), as well as DNA sequencing to assess driver and clonal mutations in different tumour habitats and guide design of appropriate cfDNA probes. The targeted cfDNA capture panel will be employed to determine the ability of cfDNA to capture the full genetic diversity observed in the primary tumour and to relate cfDNA to various imaging habitats in both renal and pancreatic cancer. Data integration will be undertaken between imaging habitats and molecular profiles for each region.

  Study EGAS00001003703

• Mutational landscape and patterns of clonal evolution in relapsed pediatric acute lymphoblastic leukemia

  Relapse of acute lymphoblastic leukemia (ALL) remains a leading cause of childhood death. Prior studies have shown clonal mutations at relapse often arise from relapse-fated subclones that exist at diagnosis. However, the genomic landscape, evolutionary trajectories and mutational mechanisms driving relapse are incompletely understood. In an analysis of 92 cases of relapsed childhood ALL, incorporating multimodal DNA and RNA sequencing, deep digital mutational tracking and xenografting to formally define clonal structure, we identify 50 significant targets of mutation with distinct patterns of mutational acquisition or enrichment. CREBBP, NOTCH1, and Ras signaling mutations rose from diagnosis subclones, whereas variants in NCOR2, USH2A and NT5C2 were exclusively observed at relapse. Evolutionary modeling and xenografting demonstrated that relapse-fated clones were minor (50%), major (27%) or multiclonal (18%) at diagnosis. Putative second leukemias, including those with lineage shift, were shown to most commonly represent relapse from an ancestral clone rather than a truly independent second primary leukemia. A subset of leukemias prone to repeated relapse exhibited hypermutation driven by at least three distinct mutational processes, resulting in heightened neoepitope burden and potential vulnerability to immunotherapy. Finally, relapse-driving sequence mutations were detected prior to relapse using deep digital PCR at levels comparable to orthogonal approaches to monitor levels of measurable residual disease. These results provide a genomic framework to anticipate and circumvent relapse by earlier detection and targeting of relapse-fated clones.

  Study EGAS00001003975

• The Southern African Human Genome Programme

  The Southern African Human Genome Programme (SAHGP) is a national and regional initiative that aims to unlock the unique genetic character of southern African populations. Its vision is to improve quality of life through understanding human genetic diversity. It is timely to advance the SAHGP initiative since global genomic initiatives are increasingly pointing to populations in Africa, more specifically the San and Khoe from southern Africa, as harbouring the most ancient genetic signatures among the living populations of the world. Studies on African populations have already generated a wealth of information upon which the SAHGP can build. This information will be used to promote and support genomic research programmes both nationally and regionally to address critical questions that would benefit the people of our region. The SAHGP aims to make a significant contribution to the understanding of DNA variation among southern Africans and how this impacts on the health of the people of this region. Aims of the SAHGP: To develop capacity for genomic research in southern Africa; To establish sustainable resources for genomic research in the region; and To translate knowledge and information into improvements in human health in the region. The objective is to derive the maximum benefit from the SAHGP, and that the data should be widely shared with the scientific community for the benefit of the people of the sub-continent and beyond.

  Study EGAS00001002639

• DNA hypermethylation and differential gene expression associated with Klinefelter syndrome

  The molecular basis for the phenotypic traits and morbidity in Klinefelter syndrome (KS) are not clarified. As DNA methylation affect gene expression and thereby play a role in disease susceptibility, we performed genome-wide DNA methylation profiling of leucocytes from peripheral blood samples from 67 KS patients, 67 male controls and 33 female controls, in addition to genome-wide RNA-sequencing profiling in a subset of 9 KS patients, 9 control males and 13 female controls. Characterization of the methylome as well as the transcriptome of both coding and non-coding genes identified a unique epigenetic and genetic landscape of both autosomal chromosomes as well as the X chromosome in KS. A subset of genes show significant correlation between methylation values and expression values. Gene set enrichment analysis of differentially methylated positions yielded terms associated with well-known comorbidities seen in KS. In addition, differentially expressed genes revealed enrichment for genes involved in the immune system, wnt-signaling pathway and neuron development. Based on our data we point towards many new candidate genes (AKAP17A, AMOT, APOB, DACT1, DDX58, DOCK7, EIF2S3, FIGNL1, G3BP1, HENMT1, HERC5, IFI44, IFI44L, IFIT1, IFIT3, ISG15, KANK1, LGALS1, NSD1, PEX10, PLSCR1, RSAD2, SHROOM2, SLC25A6, SPEG, SPON2, TXLNG), which may be implicated in the phenotype and further point towards non-coding genes (RP13-216E22.4, RP13-36G14.4, G087825, G088512), which may be involved in X chromosome inactivation in KS and in the regulation of escape genes.

  Study EGAS00001002797

• Atypical 3q26/MECOM rearrangements genocopy inv(3)/t(3;3) in acute myeloid leukemia

  Acute myeloid leukemia (AML) with inv(3)/t(3;3)(q21q26) is a distinct WHO recognized entity, characterized by its aggressive course and poor prognosis. In this subtype of AML, the translocation of a GATA2 enhancer (3q21) to MECOM (3q26) results in overexpression of the MECOM isoform EVI1 and monoallelic expression of GATA2 from the unaffected allele. The full-length MECOM transcript, MDS1-EVI1, is not expressed as the result of the 3q26 rearrangement. Besides the classical inv(3)/t(3;3), a number of other 3q26/MECOM rearrangements with poor treatment response have been reported in AML. Here we demonstrate, in a group of 33 AML patients with atypical 3q26 rearrangements, MECOM involvement with EVI1 overexpression, but no or low MDS1-EVI1 levels. Moreover, the 3q26 translocations in these AML patients often involve super-enhancers of genes active in myeloid development (e.g. CD164, PROM1, CDK6 or MYC). In more than 50% of these cases allele specific GATA2 expression was observed, either by copy number loss or by an unexplained allelic imbalance. Altogether, atypical 3q26 recapitulate the main leukemic mechanism of inv(3)/t(3;3) AML, namely EVI1 overexpression driven by enhancer hijacking, absent MDS1-EVI1 expression and potential GATA2 involvement. Therefore, we conclude that both atypical 3q26/MECOM and inv(3)/t(3;3) can be classified as a single entity of 3q26-rearranged AMLs. Routine analyses determining MECOM rearrangements, EVI1 and MDS1-EVI1 expression are required to recognize 3q-rearranged AML cases.

  Study EGAS00001004325

• Partner-independent fusion gene detection by multiplexed CRISPR/Cas9 enrichment and long-read Nanopore sequencing

  Fusion genes are hallmarks of various cancer types and important determinants for diagnosis, prognosis and treatment possibilities. The promiscuity of fusion genes with respect to partner choice and exact breakpoint-positions restricts their detection in the diagnostic setting, even for known and recurrent fusion gene configurations. To accurately identify these gene fusions in an unbiased manner, we developed FUDGE: a FUsion gene Detection assay from Gene Enrichment. FUDGE couples target-selected and strand-specific CRISPR/Cas9 activity for enrichment and detection of fusion gene drivers (e.g. BRAF, EWSR1, KMT2A/MLL) - without prior knowledge of fusion partner or breakpoint-location - to long-read Nanopore sequencing. FUDGE encompasses a dedicated bioinformatics approach (NanoFG) to detect fusion genes from Nanopore sequencing data. Our strategy is flexible with respect to target choice and enables multiplexed enrichment for simultaneous analysis of several genes in multiple samples in a single sequencing run. We observe on average a 508-fold on-target enrichment and identify fusion breakpoints at nucleotide resolution - all within two days. We demonstrate that FUDGE effectively identifies fusion genes in cancer cell lines, tumor samples and on whole genome amplified DNA irrespective of partner gene or breakpoint-position in 100% of cases. Furthermore, we show that FUDGE is superior to routine diagnostic methods for fusion gene detection. In summary, we have developed a rapid and versatile fusion gene detection assay, providing an unparalleled opportunity for pan-cancer detection of fusion genes in routine diagnostics.

  Study EGAS00001003964

• Memory-like HCV-specific CD8+ T cells retain a chronic scar after cure of chronic HCV infection

  Chronic hepatitis C virus (HCV) infection is associated with CD8+ T-cell exhaustion characterized by limited effector functions and thus compromised anti-viral activity. Exhausted HCV-specific CD8+ T cells are comprised of memory-like and terminally exhausted CD8+ T-cell subsets. So far, little is not known about the molecular profile and fate of these cells after elimination of chronic antigen stimulation by direct acting antiviral therapy (DAA). Here, we report an antigen-driven molecular core signature underlying exhausted CD8+ T-cell subset heterogeneity in chronic viral infection with a progenitor/progeny relationship of memory-like and terminally exhausted HCV-specific CD8+ T cells via an intermediate stage. Furthermore, transcriptional profiling reveals that the memory-like cells remain after DAA-mediated cure while terminally exhausted HCV-specific CD8+ T-cell subsets are lost. Thus, the memory polarization of the overall HCV-specific CD8+ T-cell response after cure does not result from re-differentiation of exhausted T cells. Consequently, antigen elimination has little impact on the exhausted core signature of memory-like CD8+ T cells that remains clearly different from bona fide T-cell memory. These results identify a molecular signature of T-cell exhaustion that is imprinted like a chronic scar in HCV-specific CD8+ T cells even after HCV cure, highlighting the requirement of re-programming to elicit full effector potential of exhausted T cells.

  Study EGAS00001004538

• Small variants in mtDNA Canary Islands - WES Illumina (ITER)

  The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

  Study EGAS00001005678

• Small variants in mtDNA Canary Islands - WGS Illumina (ITER)

  The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

  Study EGAS00001005679

• Non-coding mutations reveal cancer driver cistromes in luminal breast cancer

  Whole-genome sequencing of primary breast tumors enabled the identification of cancer driver genes and non-coding cancer driver plexuses from somatic mutations. However, differentiating between driver and passenger events among non-coding genetic variants remains a challenge to understand the etiology of cancer and inform the delivery of personalized cancer medicine. Herein, we reveal enrichment of non-coding mutations in cis-regulatory elements that cover a subset of transcription factors linked to tumor progression in luminal breast cancers. Using a cohort of 26 primary luminal ER+PR+ breast tumors, we compiled a catalogue of ~100,000 unique cis-regulatory elements from ATAC-seq data. Integrating this catalogue with somatic mutations from 350 publicly available breast tumor whole genomes, we identified four recurrently mutated individual cis-regulatory elements. By then partitioning the non-coding genome into cistromes, defined as the sum of binding sites for a transcription factor, we uncovered cancer driver cistromes for ten transcription factors, namely CTCF, ELF1, ESR1, FOSL2, FOXA1, FOXM1 GATA3, JUND, TFAP2A, and TFAP2C in luminal breast cancer. Nine of these ten transcription factors were shown to be essential for growth in breast cancer, with four exclusive to the luminal subtype. Collectively, we present a strategy to find cancer driver cistromes relying on quantifying the enrichment of non-coding mutations over cis-regulatory elements concatenated into a functional unit drawn from an accessible chromatin catalogue derived from primary cancer tissues.

  Study EGAS00001005235

• A novel orthotopic patient-derived xenograft model of radiation-induced glioma following medulloblastoma

  Radiation-induced glioma (RIG) is a highly aggressive brain cancer arising as a consequence of radiation therapy. We report a case of RIG that arose in the brain stem following treatment for paediatric medulloblastoma, and the development and characterisation of a matched orthotopic patient-derived xenograft (PDX) model (TK-RIG915). Patient and PDX tumours were analysed using DNA methylation profiling, whole genome sequencing (WGS) and RNA sequencing. While initially thought to be a diffuse intrinsic pontine glioma (DIPG) based on disease location, results from methylation profiling and WGS were not consistent with this diagnosis. Furthermore, clustering analyses based on RNA expression suggested the tumours were distinct from primary DIPG. Additional gene expression analysis demonstrated concordance with a published RIG expression profile. Multiple genetic alterations that enhance PI3K/AKT and Ras/Raf/MEK/ERK signalling were discovered in TK-RIG915 including an activating mutation in PIK3CA, upregulation of PDGFRA and AKT2, inactivating mutations in NF1, and a gain-of-function mutation in PTPN11. Additionally, deletion of CDKN2A/B, increased IDH1 expression, and decreased ARID1A expression were observed. Detection of phospho-S6, -4EBP1 and -ERK via immunohistochemistry confirmed PI3K pathway and ERK activation. Here, we report one of the first PDX models for RIG, which recapitulates the patient disease and is molecularly distinct from primary brain stem glioma. Genetic interrogation of this model has enabled the identification of potential therapeutic vulnerabilities in this currently incurable disease.

  Study EGAS00001004709

• Small variants in mtDNA Canary Islands - WGS Oxford Nanopore Technologies (ITER)

  The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. Besides, because of its relevance, we also assess the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short- read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

  Study EGAS00001005677

• Conserved features of TERT promoter duplications reveal an activation mechanism that mimics hotspot mutations in cancer

  Mutations in the TERT promoter are the single most common non-coding mutation in cancer and represent the genetic underpinnings of tumor cell immortality. Beyond the two most common point mutations, G228A and G250A, which selectively recruit the ETS factor GABP to activate TERT, the significance of other variants in the TERT promoter are unknown. We identified duplications of wildtype sequence within the core promoter region of TERT in 7 different cancer types that have strikingly similar features including size, insertion position, and inclusion of one of the native ETS motifs. Each duplication activates the TERT promoter to a similar level as G228A and G250A and is critically dependent on the insertion site. The GABP tetramer binds to the TERT duplicated promoter sequence by virtue of the native ETS motif and its duplicated version with precise spacing, and it is necessary for the transcriptional activation by all duplications tested. Spatiotemporal analysis in a multifocal glioblastoma shows the duplication is clonal and its activation of TERT is readily detectable at the single cell level and in bulk tumor tissue. We conclude that recurrent TERT promoter duplications of the native ETS sequence are functionally and mechanistically equivalent to the hotspot mutations that confer tumor cell immortality. The shared mechanism of these divergent somatic genetic alterations suggests a strong selective pressure for recruitment of the GABP tetramer to activate TERT.

  Study EGAS00001006118

• Prostate cancer ancestral genomic disparity

  Prostate cancer is characterised by significant global disparity; mortality rates in SubSaharan Africa are double to quadruple those in Eurasia. Hypothesising unknown interplay between genetic and non-genetic factors, tumour genome profiling envisages contributing mutational processes. Through whole-genome sequencing of treatment-naïve prostate cancer from 183 ethnically/globally distinct patients (African versus European), we generate the largest cancer genomics resource for Sub-Saharan Africa. Identifying ~2 million somatic variants, Africans carried the greatest burden. We describe a new molecular taxonomy using all mutational types and ethno39 geographic identifiers, including Asian. Defined as Global Mutational Subtypes (GMS) A–D, although Africans presented within all subtypes, we found GMS-B to be ‘African-specific’ and GMS-D ‘African-predominant’, including Admixed and European Africans. Conversely, Europeans from Australia, Africa and Brazil predominated within ‘mutationally-quiet’ and ethnically/globally ‘universal’ GMS-A, while European Australians shared a higher mutational burden with Africans in GMS45 C. GMS predicts clinical outcomes; reconstructing cancer timelines suggests four evolutionary trajectories with different mutation rates (GMS-A, low 0.968/year versus D, highest 1.315/year). Our data suggest both common genetic factors across extant populations and regional environmental factors contributing to carcinogenesis, analogous to gene-environment interaction defined here as a different effect of an environmental surrounding in persons with different ancestries or vice versa. We anticipate GMS acting as a proxy to intrinsic and extrinsic mutational processes in cancers, promoting global inclusion in landmark studies.

  Study EGAS00001006425

• Disease and phenotype relevant genetic variants identified from histone acetylomes in human hearts

  Identifying genetic markers for heterogeneous complex diseases such as heart failure has been challenging, and may require prohibitively large cohort sizes in genome-wide association studies (GWAS) in order to demonstrate statistical significance1. On the other hand, chromatin quantitative trait loci (QTL), elucidated by direct epigenetic profiling of specific human tissues, may contribute towards prioritising variants for disease-association. Here, we captured non-coding genetic variants by performing enhancer H3K27ac ChIP-seq in 70 human control and end-stage failing hearts, mapping out a comprehensive catalogue of 47,321 putative human heart enhancers. 3,897 differential acetylation peaks (FDR < 0.05) pointed to recognizable pathways altered in heart failure (HF). To identify cardiac histone acetylation QTLs (haQTLs), we regressed out confounding factors including HF disease status, and employed the G-SCI test2 to call out 1,680 haQTLs (FDR < 0.1). A subset of these showed significant association to gene expression, either in cis (180), or through long range interactions (81), identified by Hi-C and Hi-ChIP performed on a subset of hearts. Furthermore, a concordant relationship was found between the gain or disruption of specific transcription factor (TF) binding motifs, inferred from alternative alleles at the haQTLs, associated with altered H3K27ac peak heights. Finally, colocalisation of our haQTLs with heart-related GWAS datasets allowed us to identify 62 unique loci. Disease-association for these new loci may indeed be mediated through modification of H3K27-acetylation enrichment and their corresponding gene expression differences.

  Study EGAS00001003586

• Targeting the p53 pathway to treat Malignant Rhabdoid and High-Risk Atypical Teratoid Rhabdoid Tumors

  Rhabdoid tumors (RTs) predominantly affect young children and are among the deadliest pediatric solid tumors. Despite multimodal therapy consisting of surgery, radiation, and chemotherapy, children with these tumors have median survival of less than one year. RTs can arise throughout the body, including the central nervous system (CNS) where they are called atypical teratoid rhabdoid tumors (AT/RTs), and in extra-CNS locations such as the kidneys and other soft tissues where they are designated malignant RTs (MRTs). We previously identified MDM2 and MDM4 as therapeutic vulnerabilities in RTs and showed that treatment with the MDM2 inhibitor idasanutlin increased survival in mice bearing MRT xenografts. However, the therapeutic potential of idasanutlin in CNS RT tumors is unknown. Single agent therapies are prone to resistance and show limited clinical benefit on their own. Therefore, we sought to identify combination strategies incorporating idasanutlin that would be effective in both CNS and extra-CNS RTs. Here, we show that the XPO1 inhibitor selinexor increased nuclear retention of p53 and potentiated idasanutlin-induced p53 pathway activation and cytotoxicity in AT/RT and MRT cell lines in vitro. Importantly, combination therapy limited acquired resistance through TP53 mutation. In vivo, combination therapy was well-tolerated, reduced tumor burden, and increased survival in orthotopic models of both AT/RT and MRT. Our results demonstrate that combining idasanutlin with selinexor is a promising therapeutic strategy for children with rhabdoid tumors.

  Study EGAS00001007680

• cfMethyl-seq data (cfSort study) from serial plasma samples of NSCLC patients

  The cfMethyl-Seq libraries of the serial plasma cfDNA samples from the four NSCLC patients were constructed following the standard protocol. 10 ng of cfDNA in the volume of 25 µl was used as input material. 5’-end dephosphorylation was done with 3 µl 10xCutSmart buffer and 2 µl quick CIP from NEB (Ipswich, MA) at 37 oC for 30 min then heat-inactivated at 80 oC for 5 min. The 3’-end blocking was done with 0.5 µl 10xCutSmart buffer, 3 µl 2.5 mM CoCl2, 1 µl terminal transferase (all from NEB), and 0.5 µl 1 mM ddGTP at 37 oC for 2 h followed by 75 oC for 20 min. The mixture was then purified with 2x AmpureXP beads (Beckman Coulter, Indianapolis, IN) and eluted in 21.5 µl RT-PCR grade water (Thermo-Fisher, Waltham, MA). Restriction digestion was done with 2.5 µl 10xCutSmart buffer and 1 µl MspI (NEB) for 18 h at 37  oC and 20 min at 65 oC . 0.5 µl 10xCutSmart buffer, 0.3 µl dACGTP mixture (100 mM dATP, 10 mM dCTP, 10 mM dGTP), 1 µl Klenow (exo-, 5U/µl, NEB) and 2.6 µl RT-PCR water, 0.6 µl 50 mM DTT (ThermoFisher) was added to the mixture for end repair and A-overhang addition with the program 30  oC for 20 min, 37 oC for 1 h and 75 oC for 20 min. Adapter ligation was then performed with 1 µl 10xThermoFisher HC T4 ligase buffer, 0.4 µl 100 mM ATP (ThermoFisher), 0.2 µl 50 mM DTT, 1 µl ThermoFisher HC T4 DNA ligase (30 Weiss Unit/µl), 5 ng home-made duplex UMI adapter with all the cytosines methylated (protocol adopted from Kennedy et al.) at 16 oC for 20 h and 65 oC for 20 min. Bisulfite conversion of the adapter-ligated product was carried out with QIAGEN EpiTect plus DNA bisulfite kit following their protocol for two rounds of conversion. The converted product was purified with Qiagen MinElute spin column and eluted with 20 µl RT-PCR water. PCR amplification was done using the NEBNext Multiplex Oligos for Illumina (2.5 µl of universal and index primer each) and 25 µl KAPA HiFi HotStart Uracil+ ReadyMix (Roche) with the following cycling conditions: 98  oC for 45 s, 15 cycles of 98  oC for 15 s, 60  oC for 30 s and 72  oC for 30 s, followed by a final extension at 72  oC for 5 min. The PCR product was purified with 1x AmpureXP beads and eluted with 30 µl EB buffer. DNA concentration was measured by Qubit 1xdsDNA HS assay. 5% TBE-UREA PAGE and bioanalyzer assay was performed as quality control on each library before sequencing.

  Dataset EGAD00001010881

• WGS and WES data of invasive micropapillary carcinoma of breast

  This study contains whole genome sequencing data and whole exon sequencing data of IMPC tumor and normal tissue sample.

  Dataset EGAD00001008450

• NIPT samples for systematic evaluation of BinDel, a software tool for detecting clinically significant microdeletions in low-coverage WGS-based NIPT samples

  This dataset contains three sets of samples. The first sample set contains euploid fetus pregnancies reported by NIPTIFY screening test and postnatal evaluation. Dataset was processed similarly to previously published guidelines from KU Leuven, with modifications [1]. Briefly, peripheral blood samples were collected in cell-free DNA BCT tubes (Streck, USA), and plasma was separated with standard dual centrifugation. Cell-free DNA was extracted from 3 ml plasma using MagMAX Cell-Free DNA Isolation Kit (ThermoFisher Scientific). Whole-genome libraries were prepared using the FOCUS (Fragmented DNA Compact Sequencing Assay, Competence Centre on Health Technologies, Estonia) NIPT method protocol with 12 cycles for the final PCR enrichment step. In the following quantification, equal amounts of 36 samples were pooled, and the quality and quantity of the pool were assessed on Agilent 2200 TapeStation (Agilent Technologies, USA). Whole genome sequencing was performed on the NextSeq 550 instrument (Illumina Inc.) with an average coverage of 0.32× (minimum 0.08 and maximum 0.42) and producing 85 bp single-end reads. The second sample set contains a single NIPT sample postnatally diagnosed with Prader-Willi syndrome. The sample was sequenced with Illumina NextSeq 500 platform, producing 85 bp single-end reads with an average per-sample coverage of 0.32× at the University of Tartu, Institute of Genomics Core Facility, according to the manufacturer’s standard protocols, as described previously [2]. The third sample set contains samples SC005 (SeraCare Life Sciences Inc lot #10446565), SC0042 (#10571706), and SC016 (#10560229). These are SeraCare Life Sciences Inc circulating cell-free DNA (ccfDNA) like mixture of human genomic DNA that consists of matched maternal and fetus. SC005 and SC0042 consist of matched DNA of maternal and fetus with DiGeorge Syndrome. SC016 is a custom-ordered DNA Mix with fetus DNA having a pathogenic loss of the terminal region of 20p13 and a pathogenic 3q29 duplication. SC016 was processed as the first sample set was processed, and SC0042 was processed as the second sample set was processed. Sample SC005 was processed once as was sample set 1 and once as was sample set 2 processed. This study was performed with the approval of the Research Ethics Committee of the University of Tartu (#352/M-12). 1. Bayindir B, Dehaspe L, Brison N, Brady P, Ardui S, Kammoun M, et al. Noninvasive prenatal testing using a novel analysis pipeline to screen for all autosomal fetal aneuploidies improves pregnancy management. Eur J Hum Genet. 2015;23: 1286– 1293. doi:10.1038/ejhg.2014.282 2. Žilina O, Rekker K, Kaplinski L, Sauk M, Paluoja P, Teder H, et al. Creating basis for introducing noninvasive prenatal testing in the Estonian public health setting. Prenat Diagn. 2019;39: 1262–1268. doi:10.1002/pd.5578

  Dataset EGAD00001009512

• Bacterial SNPs in the human gut microbiome associate with host BMI

  Genome-wide association studies (GWAS) have revolutionized the field of genetics, providing numerous associations between human SNPs and health traits. Likewise, metagenome-wide association studies (MWAS) between bacterial SNPs and human traits can suggest mechanistic links between the microbiota and its host. However, very few such studies have been done to date, in part due to the large number of samples required to address the variable coverage of bacteria across samples and the fact that most bacteria are present in only a subset of the population. Here, we devised an MWAS framework to detect SNPs and associate them with host phenotypes systematically. We recruited and obtained gut metagenomic samples from a cohort of 7,190 healthy individuals and discovered 1,358 statistically significant associations between a bacterial SNP and host body mass index (BMI). To address the possible effects of population structure and linkage disequilibrium, we applied a clumping procedure, and found 40 independent associations between SNPs and host BMI. We uncovered BMI-associated SNPs in the genomes of 27 bacterial species, even though the relative abundance of 12 of these species was not associated with the phenotype. We separated genome-wide sub-species variations from local associations in individual SNPs, and demonstrated how this framework can highlight specific bacterial functions of interest. Our results demonstrate the importance of considering nucleotide-level diversity in microbiome studies and pave the way toward a better understanding of interpersonal gut microbiome differences and their health implications.

  Study EGAS00001007204

• Atypical B cells and impaired SARS-CoV-2 neutralisation following heterologous vaccination in the elderly

  Age is a major risk factor for hospitalization and death after SARS-CoV-2 infection, even in vaccinees. Suboptimal responses to a primary vaccination course have been reported in the elderly, but there is little information regarding the impact of age on responses to booster third doses. Here we show that individuals 70 or older who received a primary two dose schedule with AZD1222 and booster third dose with mRNA vaccine achieved significantly lower neutralizing antibody responses against SARS-CoV-2 spike pseudotyped virus compared to those below 70 at one month post booster. Neither the concentration of serum binding anti spike IgG antibody, nor the frequency of spike- specific B cells showed differences by age grouping. Impaired neutralization potency and breadth post-third dose in the elderly was associated with enrichment of circulating “atypical” spike-specific B cells expressing CD11c and FCRL5. Presence of these atypical B cells was confirmed by scRNA seq. However, when considering individuals who received three doses of mRNA vaccine, we did not observe statistically significant differences in neutralization or enrichment in atypical B cells by age grouping following dose 3. Together, our data reveal impaired neutralising antibody responses in the elderly after heterologous vaccination with AdV followed by mRNA booster dose, but not in homologous three dose mRNA vaccination highlighting that different vaccine formats differentially instruct the memory B cell response.

  Study EGAS00001007385

• Oxidative phosphorylation is a key ontogenetic feature of monocyte immunometabolism promoting myeloid differentiation after birth

  Neonates primarily rely on innate immunity, yet their inflammatory responses to microbes, particularly those of monocytes, are usually restricted compared to adults. This is controversially interpreted as immaturity increasing the risk of sepsis or essential programming allowing immune adaptation after birth. Which changes the cellular immunometabolism undergoes after birth and its role in these concepts are poorly defined. Here, we applied transcriptomic, metabolic, and immunological approaches and found that monocytes exhibit an inverse ontogenetic balance of immunometabolism. Our data show that glycolysis in monocytes increases within the first year of life and fuels their inflammatory responsiveness. In contrast, oxidative phosphorylation (OXPHOS) is high in neonatal monocytes supporting myeloid differentiation but declines only gradually during the first five years of life. Treatment of neonatal monocytes with lipopolysaccharide induces an adult-like immunometabolic phenotype. Ketogenic diet restricts glycolysis in adult monocytes, but cannot reactivate OXPHOS and revive a neonatal-like immunometabolic phenotype, suggesting that neonatal metabolism is hardwired and cannot be restored by simple dietary changes. Transcriptional network and population-wide human variation analyses identified E2F1, MYB, STAT1 and FLI1 as important regulators of age-dependent cell functions and related energy demands in human monocytes. Collectively, our findings show that restricted glycolysis and increased OXPHOS are a physiological programming in neonatal monocytes. Premature switching to an adult-like metabolism could untimely enhance inflammatory responses and disrupt important myeloid differentiation, whereas microbial challenges accumulating during childhood seem to induce an essential gradual metabolic switch.

  Study EGAS00001007555

• Ultra-long sequencing for contiguous haplotype resolution of the human immunoglobulin heavy chain locus

  Genetic diversity within the human immunoglobulin heavy chain (IGH) locus influences the expressed antibody repertoire and susceptibility to infectious and autoimmune diseases. However, repetitive sequences and complex structural variation pose significant challenges for large-scale characterization. Here, we introduce a method that combines Oxford Nanopore ultra-long sequencing and adaptive sampling with a bioinformatic pipeline to produce haplotype-resolved, annotated IGH assemblies. Notably, our strategy overcomes prior limitations in phasing resolution, enabling single-contig haplotype assemblies that span the entire IGH locus. We applied this method to four individuals and validated the accuracy of the IGH assemblies using Pacific Biosciences HiFi reads, demonstrating near-complete sequence congruence, detecting only indels with a high degree of confidence. Moreover, when applied to the reference material HG002, our pipeline revealed no base differences and a limited number of indels compared with the Telomere-to-Telomere genome benchmark across the IGH region. Importantly, in the four individuals, our approach uncovered 30 novel alleles and previously uncharacterized large structural variants, including a 120 kb duplication spanning IGHE to IGHA1 and an expanded seven-copy IGHV3-23 gene haplotype. These findings underscore the power of our method to resolve the full complexity of the IGH locus and uncover previously unrecognized variants that may affect immune function and disease susceptibility. Thus, our method provides a strong basis for future immunological research and translational applications.

  Study EGAS50000001042

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

  The structure, fragmentation pattern, length and terminal sequence of cell-free DNA (cfDNA) is under the influence of nucleases present in the blood. We hypothesized that differences in the diversity of bases at the end of cfDNA fragments can be leveraged on a genome-wide scale to enhance the sensitivity for detecting the presence of tumor signals in plasma. We surveyed the cfDNA termini in 72 plasma samples from 319 patients with 18 different cancer types using low-coverage whole genome sequencing. The fragment-end sequence and diversity were altered in all cancer types in comparison to 76 healthy controls. We converted the fragment end sequences into a quantitative metric and observed that this correlates with circulating tumor DNA tumor fraction (R = 0.58, p < 0.001, Spearman). Using these metrics, we were able to classify cancer samples from control at a low tumor content (AUROC of 91% at 1% tumor fraction) and shallow sequencing coverage (mean AUROC = 0.99 at >1M fragments). Combining fragment-end sequences and diversity using machine learning, we classified cancer from healthy controls (mean AUROC = 0.99, SD = 0.01). Using unsupervised clustering we showed that early-stage lung cancer can be classified from control or later stages based on fragment-end sequences. We observed that fragment-end sequences can be used for prognostication (hazard ratio: 0.49) and residual disease detection inresectable esophageal adenocarcinoma patients, moving fragmentomics toward a greater clinical implementation.

  Study EGAS00001005738

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

  The structure, fragmentation pattern, length and terminal sequence of cell-free DNA (cfDNA) is under the influence of nucleases present in the blood. We hypothesized that differences in the diversity of bases at the end of cfDNA fragments can be leveraged on a genome-wide scale to enhance the sensitivity for detecting the presence of tumor signals in plasma. We surveyed the cfDNA termini in 72 plasma samples from 319 patients with 18 different cancer types using low-coverage whole genome sequencing. The fragment-end sequence and diversity were altered in all cancer types in comparison to 76 healthy controls. We converted the fragment end sequences into a quantitative metric and observed that this correlates with circulating tumor DNA tumor fraction (R = 0.58, p < 0.001, Spearman). Using these metrics, we were able to classify cancer samples from control at a low tumor content (AUROC of 91% at 1% tumor fraction) and shallow sequencing coverage (mean AUROC = 0.99 at >1M fragments). Combining fragment-end sequences and diversity using machine learning, we classified cancer from healthy controls (mean AUROC = 0.99, SD = 0.01). Using unsupervised clustering we showed that early-stage lung cancer can be classified from control or later stages based on fragment-end sequences. We observed that fragment-end sequences can be used for prognostication (hazard ratio: 0.49) and residual disease detection inresectable esophageal adenocarcinoma patients, moving fragmentomics toward a greater clinical implementation.

  Study EGAS00001005747

• Non-muscle Invasive Bladder Cancer Molecular Subtypes Predict Differential Response to Intravesical Bacillus Calmette-Guérin

  The recommended treatment for patients with high-risk non-muscle invasive bladder cancer (HR-NMIBC) is tumor resection followed by adjuvant Bacillus Calmette-Guérin (BCG) bladder instillations. However, only 50% of patients benefit from this therapy. In case of progression to advanced disease, patients must undergo a radical cystectomy with significant morbidity and have a poor clinical outcome. Identifying tumors least likely to respond to BCG can translate into alternative treatments, such as early radical cystectomy or novel targeted or immunotherapies. Here we present molecular profiling of 132 BCG-naive, HR-NMIBC patients, and 44 post-BCG recurrences (34 matched), which uncovered three distinct BCG Response Subtypes (BRS1-3). Patients with BRS3 tumors have reduced recurrence and progression-free survival compared to BRS1-2. BRS3 tumors expressed high EMT-basal markers and had an immunosuppresive profile, which was confirmed with spatial proteomics. Tumors which recurred post-BCG were enriched for BRS3. BRS stratification was validated in a second cohort of 151 BCG-naive HR-NMIBC patients and the molecular subtypes outperformed guideline recommended risk stratification based on clinicopathological variables. For clinical application, we validated that a commercially approved assay was able to accurately predict BRS3 tumors (AUROC 0.86). Our findings provide a potential clinical tool for improved identification of HR-NMIBC patients at the highest risk of progression, which can be used to select patients for early radical cystectomy or novel subtype-directed therapies.

  Study EGAS00001006879

• H3Africa - Stroke Investigative Research and Education Networks

  Although stroke is a discrete phenotype, it is the clinical culmination of several complex and interacting biological processes, precipitated by various genetic and environmental factors, thereby making ready analyses of its underlying mechanisms a challenge. However, the formulation of successful tailor-made prevention strategies for limiting the immense personal and societal burden of stroke in Sub-Saharan Africa (SSA) is contingent on a better understanding of the predisposing risk factors for stroke in the region, which may not necessarily match the nature of stroke risk factors described among other people on other continents. Indeed, the underlying risk factors for stroke in SSA have hitherto not been comprehensively characterized, and currently there is a paucity of expertise, experience, and infrastructure in genomics and environmental research in SSA, to initiate and sustain such investigation. We seek to establish the Stroke Investigative Research and Education Network (SIREN) as a multidisciplinary collaborative research network of investigators in SSA and the United States collectively focused on exploring ways to promote a better understanding of the genetic and environmental risk factors for stroke among people of African ancestry. SIREN will evaluate unique contributors to stroke in SSA and compare these findings with those in a cohort of African American stroke survivors, all while simultaneously building sustainable capacities in phenomics, biobanking, genomics, biostatistics and bioinformatics for future high-level investigation of stroke and other vascular disease entities in SSA.

  Study EGAS00001007331

• Prognostic and therapeutic significance of leukemia subtypes and minimal residual disease measurements in pediatric acute lymphoblastic leukemia treated with contemporary risk-directed trial: a cohort study

  Background Genomic analyses to identify leukemia subtypes and refinement of minimal residual disease (MRD) determination are two recent advances that would have substantial impact on the management of children with acute lymphoblastic leukemia (ALL). However, the combined application of these two approaches to determine the prognostic and therapeutic significance of MRD measurement and full spectrum of leukemia subtypes in the context of contemporary risk-directed therapy has not been studied comprehensively.MethodsBetween October 29, 2007 and March 26, 2017, 598 consecutive children with newly diagnosed ALL were enrolled in a clinical trial at St. Jude Children’s Research Hospital. Patients were provisional classified according to the presenting clinical and laboratory features including immunophenotype and cytogenetics. Patients with ETV6-RUNX1 and high hyperdiploid B-progenitor ALL were provisionally included in the low-risk groups; those with TCF3-PBX1, hypodiploidy and T-cell ALL in the standard (intermediate)-risk group; and patients with BCR-ABL1, KMT2A and early T-cell precursor ALL in the high-risk group. Final risk assignment was based primarily on MRD levels measured by flow cytometry on Day 15 (≥1%) and Day 42 (<0.01% or ≥1%) of remission induction. Using genomic classification based on cytogenetics and genomic analysis (single nucleotide polymorphism array analysis and transcriptome, exome, and genome sequencing), patients were classified into 17 subtypes. The primary aim was to determine the event-free survival, overall survival and cumulative risk of relapse according to specific leukemia subtypes and MRD levels on Days 8, 15 and 42 of remission induction. This clinical trial was registered at ClinicalTrials.gov, number NCT00549848.FindingsPatients with ETV6-RUNX1, high hyperdiploidy and DUX4-rearranged B-ALL all received low- or standard-risk therapy, and had the highest 5-year event-free survival rates: 98.4% (95% CI, 95.9-100), 95.3% (91.2-99.4), and 95.0% (84.2-100), respectively. All 142 patients with Day 8 MRD<0.01%, except two with KMT2A-rearranged ALL and one with TCF3-PBX1, were in continuous complete remission. The TCF3-PBX1, PAX5alt, T-cell, early T-cell precursor, iAMP21 and hypodiploid ALL had intermediate outcome with 5-year event-free survival rates ranging from 80.0% (39.4-100) to 88.2% (71.7-100). Patients with KMT2A-rearranged, BCR-ABL1, BCR-ABL1-like and ETV6-RUNX1-like ALL had the worst outcomes, with 5-year event-free survival rates between 64.1% (43.9-84.3) and 76.2% (51.9-100), despite most patients receiving standard- or high-risk therapy. Treatment intensification based on MRD and genotype improved outcome. Relapse did not occur in any of the patients with ETV6-RUNX1, DUX4-rearranged, TCF3-PBX1, iAMP21, ZNF384-rearranged, hypodiploid, or BCR-ABL1-like B-ALL who had day 15 MRD≥1% and received standard-risk therapy. However, day 42 MRD<0.01% did not preclude relapse in patients with intermediate-risk or unfavorable subtypes including TCF3-PBX1, PAX5alt, T-cell, iAMP21, BCR-ABL1, BCR-ABL1-like, ETV6-RUNX1-like, KMT2A-rearranged, and MEF2D-rearranged ALL, despite intensified treatment in the standard- or high-risk arms. InterpretationContemporary risk-directed treatment based on genotypes and MRD can cure virtually all patients with ETV6-RUNX1, hyperdiploid>50 and DUX4-rearranged ALL. While intensified treatment improve outcome for patients with high MRD during induction, substantial proportions of patients with intermediate-risk or unfavorable subtypes would still relapse even achieving negative MRD status at the end of induction and require novel therapeutics to improve outcome.

  Study EGAS00001004739

• Evolution of an adenomacarcinoma in response to selection by targeted kinase inhibitors

  Background - Adenocarcinomas of the tongue are rare and represent the minority (20-25%) of salivary gland tumors affecting the tongue. We investigated the utility of massively parallel sequencing to characterize an adenocarcinoma of the tongue, before and after treatment. Results - In the pre-treatment tumor we identified 7,629 genes within regions of copy number gain. 1,078 genes exhibited increased expression relative to the blood and unrelated tumors and four genes contained somatic protein-coding mutations. Our analysis suggested the tumor cells were driven by the RET oncogene. Genes whose protein products are targeted by the RET inhibitors sunitinib and sorafenib correlated with being amplified and or highly expressed. Consistent with our observations administration of sunitinib was associated with stable disease lasting 4 months, after which the lung lesions began to grow. Administration of sorafenib and sulindac provided disease stabilization for an additional 3 months after which the cancer progressed and new lesions appeared. A recurring metastasis possessed 7,288 genes within copy number amplicons, 385 genes exhibiting increased expression relative to other tumours and 9 new somatic protein coding mutations. The observed mutations and amplif ications were consistent with therapeutic resistance arising through activation of the MAPK and AKT pathways. Conclusions - We conclude that complete genomic characterization of a rare tumour has the potential to aid in clinical decision making and identifying therapeutic approaches where no established treatment protocols exist. These results also provide direct in-vivo genomic evidence for mutational evolution within a tumor under drug selection and potential mechanisms of drug resistance accrual.

  Study EGAS00000000074

• Pharmacological improvement of CFTR function rescues airway epithelial homeostasis and host defense in cystic fibrosis children

  Rationale: Recent clinical trials have shown that elexacaftor/tezacaftor/ivacaftor (ETI) provides significant benefit to patients with cystic fibrosis (CF) by improving CFTR function. Despite these promising clinical findings, the specific effects on the airway inflammatory status leading to structural damage and impaired lung function, remain elusive. Objectives: We investigated the transcriptional changes at the single cell level of airway epithelial and immune cells that underlie the clinical improvement by ETI therapy in children with CF. Methods: Nasal swabs from thirteen children with CF and at least one F508del allele aged 6 to 12 years were collected at baseline and three months after initiation of ETI and subjected to scRNA-seq. In addition, we collected nasal swabs from 12 age- and sex-matched controls. Sweat chloride concentrations, spirometry and multiple breath washout were used for clinical evaluation of CF children. Measurements and main results: ETI significantly reduced sweat chloride concentrations (-52.2 ± 14.3 mmol/L, P < 0.001) and improved the lung clearance index (-1.3 ± 1.8, P < 0.05). Single cell transcriptomics revealed an impaired interferon signalling in epithelial cells of CF children at baseline, which was partially restored by ETI in conjunction with an ETI-induced increase in expression of MHC I and II encoding genes similar to control levels. Additionally, ETI markedly reduced the inflammatory phenotype of immune cells, particularly of neutrophils and macrophages. Conclusions: Improvement of CFTR function by ETI restores epithelial homeostasis and reduces immune cell inflammatory responses in the upper airways of children with CF, highlighting the potential of early initiation of ETI therapy.

  Study EGAS50000000128

• Systemic Lupus Erythematosus Serum Stimulation of Human Intestinal Organoids Induces Barrier Leakiness and Changes in Goblet Cell Differentiation

  Human intestinal epithelial cells are the interface between potentially harmful luminal content and basally residing immune cells. Their role is not only nutrient absorption but also the formation of a tight monolayer that constantly secretes mucus creating a multi-layered protective barrier. Alterations in this barrier can lead to increased gut permeability which is frequently seen in individuals with chronic extraintestinal autoimmune diseases, such as Systemic Lupus Erythematosus (SLE). Despite recent advances in identifying alterations in gut microbiota composition in SLE patients, not much attention has been given to the epithelial barrier itself. To date, it remains largely unexplored which role and function intestinal epithelial cells have in SLE pathology. Here, we present a unique near-physiologic in vitro model specifically designed to examine the effects of SLE on the epithelial cells. We utilize human colon organoids that are stimulated with serum obtained from SLE patients. Combining bulk and single-cell RNA transcriptomic analysis with functional assays revealed that SLE serum stimulation induced a unique expression profile marked by a reduction of goblet cell marker expression and mucus composition. Additionally, organoids exhibited imbalanced cellulat composition along with enhanced paracellular permeability,mitochondrial function, and a type I interferon gene signature. Similarly, transcriptomic analysis of SLE human colon biopsies revealed a downregulation of epithelial secretory markers. Our work uncovers a crucial connection between SLE and intestinal homeostasis that might be promoted in vivo through the blood, offering insights into the causal connection of barrier dysfunction and autoimmune diseases.

  Study EGAS50000000012

• Carboxylesterase 1 mediates a distinctive metabolic profile of dendritic cells to attain an inflammatory phenotype

  The metabolic profile of dendritic cells (DCs) shapes their phenotype and functions. Carboxylestrase 1 (CES1) enzyme is highly expressed in mononuclear myeloid cells however its exact role in DCs is elusive. We used a CES1 inhibitor (WWL113) and genetic overexpression to explore the role of CES1 in DCs differentiation in inflammatory models. CES1 expression was analyzed during CD14+ monocytes differentiation to DCs (MoDCs) using quantitative PCR. CES1 Inhibitor (WWL113) was applied during MoDCs differentiation. Surface markers, secreted cytokines, lactic acid production, phagocytic and T cell polarization capacity were analyzed. Transcriptomic and metabolic profile were assessed with RNA-sequencing and mass spectrometry. Cellular respiration was assessed with seahorse respirometry. Transgenic mice were used to assess CES1 overexpression in DCs in inflammatory models. CES1 expression peaks early during MoDCs differentiation. Pharmacological inhibition of CES1 led to higher expression of CD209, CD86 and MHCII. WWL113 treated MoDCs secreted higher quantities of IL6, IL8, TNF and IL10 and demonstrated stronger phagocytic ability and higher capacity to polarize Th17 differentiation in autologous DCs-T cells co-culture model. Transcriptomic profiling revealed enrichment of multiple inflammatory and metabolic pathways. Functional metabolic analysis shows impaired maximal mitochondrial respiration capacity, increased lactate production and decreased intracellular amino acids and TCA intermediates. Transgenic human CES1 overexpression in murine DCs generated less inflammatory phenotype and increased resistance to T cell mediated colitis. In conclusion, CES1 inhibition directs DCs differentiation towards more inflammatory phenotype, that shows stronger phagocytic capacity and supports Th17 skewing. This is associated with disrupted mitochondrial respiration and amino acids depletion.

  Study EGAS50000000230

• A Single Cell Transcriptomic Analysis of Human Neocortical Development

  Defining the number, proportion, or lineage of distinct cell types in the developing human brain is an important goal of modern brain research. We produced single cell transcriptomic profiles for 40,000 cells at mid-gestation to define deep expression profiles corresponding to all known major cell types at this developmental period and compare this with bulk tissue profiles. We identified multiple transcription factors (TFs) and co-factors expressed in specific cell types, including multiple new cell-type-specific relationships, providing an unprecedented resource for understanding human neocortical development and evolution. This includes the first single-cell characterization of human subplate neurons and subtypes of developing glutamatergic and GABAergic neurons. We also used these data to deconvolute single cell regulatory networks that connect regulatory elements and transcriptional drivers to single cell gene expression programs in the developing CNS. We characterized major developmental trajectories that tie cell cycle progression with early cell fate decisions during early neurogenesis. Remarkably, we found that differentiation occurs on a transcriptomic continuum, so that differentiating cells not only express the few key TFs that drive cell fates, but express broad, mixed cell-type transcriptomes prior to telophase. Finally, we mapped neuropsychiatric disease genes to specific cell types, implicating dysregulation of specific cell types in ASD, ID, and epilepsy, as the mechanistic underpinnings of several neurodevelopmental disorders. Together these results provide an extensive catalog of cell types in human neocortex and extend our understanding of early cortical development, human brain evolution and the cellular basis of neuropsychiatric disease.

  Study phs001836

• SLCO1B1 Variants and Methotrexate Clearance

  Methotrexate plasma concentration is related to its clinical effects. To identify the genetic basis of interindividual variability in methotrexate pharmacokinetics in children with newly diagnosed acute lymphoblastic leukemia (ALL), we performed a genome-wide analysis (GWAS) of 500,568 germline single-nucleotide polymorphisms (SNPs) in 434 children with ALL who received 3,014 courses of methotrexate at 2 to 5 g/m2. SNPs were validated in an independent cohort of 206 patients. Adjusting for age, race, sex, and methotrexate regimen, the most significant associations were with SNPs in the organic anion transporter polypeptide, SLCO1B1 (rs11045879 (P = 1.7 x 10-10) and rs4149081 (P = 1.7 x 10-9) (Trevino et al, PMID: 19901119). To test whether rare variants in SLCO1B1 could alter its function, we genotyped SLCO1B1 exons in a slightly larger group of 699 children with ALL who received methotrexate and identified 93 single nucleotide polymorphisms (SNPs). We found several common and rare non-synonymous (NS) SNPs associated with methotrexate clearance. NS SNPs predicted to be functionally damaging (common or rare) were more likely to be found among patients with the lowest adjusted methotrexate clearance (lowest 10%) than patients with high clearance (highest 10%). Four SLCO1B1 haplotypes were associated with reduced methotrexate clearance and we verified that these haplotypes have lower function with in vitro transport assays. We were able to quantitatively account for a third of the population variability in clearance of methotrexate with clinical and genetic covariates. This data set includes the dependent variable of methotrexate clearance and all of the SNP data available from arrays, sequencing, and genotyping.

  Study phs000426

• Investigation of the Causal Etiology in a Patient with T-B+NK+ Immunodeficiency

  We describe a male proband with normal B lymphocyte and natural killer (NK) cell numbers, but who had T lymphopenia beginning at 1 week of age. His peripheral blood contained very few T cell receptor excision circles (TREC), suggesting impaired thymic T cell development. He received an allogeneic hematopoietic stem cell marrow transplant, but despite engraftment, it did not fully correct the T lymphopenia. Sequencing excluded known immunodeficiency genes; however, whole exome sequencing revealed X-linked inheritance by the male proband of a missense mutation in MED14 (MED14V763A), a component of the mediator complex. Functional screening of the MED14V763A variant in zebrafish suggested that MED14 is essential for normal development, since morpholino (MO)-mediated loss of MED14 function attenuated T cell development, but it was rescued by ectopic expression of the wild type human MED14 but not by the MED14V763A variant. This defect in T cell repopulation is consistent with the incomplete T cell reconstitution observed upon transplantation of the proband with bone marrow from his healthy male sibling, who was found to have the same MED14V763A variant. However, the brother was normal, indicating that the altered MED14 could not by itself explain the disease. WES also revealed a mutation in the L1 cell adhesion molecule (L1CAMQ498H); however, introducing that mutation in vivo in mice did not disrupt T cell development. Consequently, the immunodeficiency in the patient is not fully elucidated. It is possible that MED14 contributes, but it is also possible that another gene is fully responsible for the disease.

  Study phs002990

• Common Variation in Candidate Genes in the Diabetes Prevention Program

  The Diabetes Prevention Program (DPP) was a multicenter controlled clinical trial examining the efficacy of an intensive lifestyle intervention or metformin to prevent or delay the development of diabetes in a population selected to be at high risk due to the presence of impaired glucose tolerance (IGT) and obesity. Development of diabetes, defined by 1997 American Diabetes Association (ADA) criteria, was the primary outcome while cardiovascular disease and its risk factors were important secondary outcomes. The pharmacological intervention was double blinded and placebo controlled. After randomization, participants had quarterly clinical evaluations and had, in addition, a fasting plasma glucose at semi-annual visits and a 75 gram oral glucose tolerance test at annual visits. Volunteers were recruited from populations known to be at particularly high risk for impaired glucose tolerance and type 2 diabetes including the following: persons with a family history of NIDDM, the elderly, overweight individuals, women with a history of diabetes during pregnancy ("gestational diabetes"), and minority populations including African Americans, Hispanic Americans, Asian and Pacific Island Americans, and Native Americans. The primary focus of the genetic investigations have been on candidate genes, including candidates derived from GWAS of diabetes, lipid and glycemia-related phenotypes. Please note more phenotype data for the DPPG cohort is available through the NIDDK Data Repository. Information on obtaining this phenotype data can be found by going to https://www.niddkrepository.org/studies/dppos/?query=DPP and https://www.niddkrepository.org/studies/dpp/?query=DPP. Data will need to be requested from both data sets. A linking file will be available to the NIDDK Data Repository.

  Study phs000681

• Role of Tobacco Smoke in Clear Cell Renal Cell Carcinoma

  Background. Clear cell renal cell carcinoma (ccRCC) is the most common histologically defined renal cancer. However, it is not a uniform disease and includes several genetic subtypes with different prognoses. ccRCC is also characterized by distinguished metabolic reprogramming. Tobacco smoking (TS) is an established risk factor for ccRCC with unknown effects on tumor pathobiology.Methods. We investigated the landscape of ccRCCs and paired normal kidney tissues (NKTs) using integrated transcriptomic, metabolomic and metallomic approaches in a cohort of never smokers (NS) and long-term current smokers (LTS) Caucasian males. Results. All three Omics domains consistently identified a distinct metabolic subtype of ccRCCs in LTS, characterized by activation of oxidative phosphorylation (OxPhos) coupled with reprogramming of the malate-aspartate shuttle and metabolism of aspartate, glutamate, glutamine and histidine. Cadmium, copper and inorganic arsenic accumulated in LTS tumors showing redistribution among intracellular pools, including relocation of copper into the cytochrome c oxidase complex. Gene expression signature based on the LTS metabolic subtype provided prognostic stratification of The Cancer Genome Atlas (TCGA) ccRCC tumors that was independent from genomic alterations.Conclusions. The work identified the TS related metabolic subtype of ccRCC with vulnerabilities that can be exploited for precision medicine approaches targeting metabolic pathways. The results provided rationale for the development of metabolic biomarkers with diagnostic and prognostic applications using evaluation of OxPhos status. The metallomic analysis revealed the role of disrupted metal homeostasis in ccRCC highlighting the importance of studying effects of metals from e-cigarettes and environmental exposures.

  Study phs002083

• Strabismus, CCDD and other anomalies

  Strabismus is one of the most common ophthalmological diseases affecting our population. It has many detrimental effects for affected individuals, including functional visual loss from amblyopia, secondary psychosocial difficulties, and limited employment opportunities. Strabismus also imposes a significant economical burden on society from the cost of screening programs, surgical correction, and lost productivity. Congenital incomitant strabismus is a broad term that encompasses congenital ocular motility disorders with restricted movement in one or more directions of gaze, and includes various forms of Duane&apos;s retraction syndrome, horizontal gaze palsy, congenital fibrosis of extraocular muscles, and Moebius syndrome. Incomitant strabismus together with a group of related disorders such as congenital ptosis and congenital facial palsy have been collectively redefined as congenital cranial dysinnervation disorders (CCDDs). This redefinition stemmed from previous studies in the Engle lab showing that these disorders can be caused by mutations in genes encoding transcription factors critical to ocular cranial motor neuron development and by genes that encode proteins essential to the proper growth and guidance of developing axons. Thus, while these disorders account for only a small percent of strabismus, understanding their genetic etiology has led to our ability to provide genetic testing and counseling, and has provided the first insights into the molecular etiologies of strabismus. While the Engle lab has identified all seven CCDD genes published to date, our cohort still includes &gt;100 familial CCDD pedigrees that are mutation-negative. We continue to utilize co-segregation analysis and whole exome or whole genome sequencing technologies to identify new disease genes in this mutation-negative cohort.

  Study phs000478

• Rb and p53-Deficient Myxofibrosarcoma and Undifferentiated Pleomorphic Sarcoma Require Skp2 for Survival

  Myxofibrosarcoma (MFS) and undifferentiated pleomorphic sarcoma (UPS) are highly genetically complex soft tissue sarcomas. Up to 50% of patients develop distant metastases, but current systemic therapies have limited efficacy. MFS and UPS have recently been shown to commonly harbor copy number alterations or mutations in the tumor suppressor genes RB1 and TP53. As these alterations have been shown to engender dependence on the oncogenic protein Skp2 for survival of transformed cells in mouse models, we sought to examine its function and potential as a therapeutic target in MFS/UPS. Comparative genomic hybridization (CGH) and next-generation sequencing (NGS) confirmed that a significant fraction of MFS and UPS patient samples (n=94) harbor chromosomal deletions and/or loss-of-function mutations in RB1 and TP53 (88% carry alterations in at least one gene; 60% carry alterations in both). Tissue microarray (TMA) analysis identified a correlation between absent Rb and p53 expression and positive expression of SKP2. Downregulation of SKP2 or treatment with Skp2 Inhibitor C1 revealed that Skp2 drives proliferation of patient-derived MFS/UPS cell lines deficient in both Rb and p53 by degrading p21 and p27. Inhibition of Skp2 using the neddylation-activating enzyme (NAE) inhibitor pevonedistat decreased growth of Rb/p53-negative patient-derived cell lines and mouse xenografts. These results demonstrate that loss of both Rb and p53 renders MFS and UPS dependent on Skp2, which can be therapeutically exploited and could provide the basis for promising novel systemic therapies for MFS and UPS.

  Study phs001982

• NHLBI GO-ESP: Family Studies (Familial Atrial Fibrillation)

  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. Large epidemiological studies have demonstrated a significant heritable component in atrial fibrillation (AF), especially the Lone forms, suggesting a monogenic syndrome. Although substantial genetic contribution has been made to the etiology of AF, the specific genes have not yet been identified. The familial form of this disease remains poorly characterized and largely undetermined. Here we seek to identify, characterize and determine the natural course of AF in our clinical practice. We identified four large multi-generation families (FAF 1-4). In FAF 1-2, most family members have symptomatic paroxysmal Atrial Fibrillation (AF) and were adequately treated with a combination of rate and rhythm therapies. By contrast, the AF substrate in FAF 3 and 4 was resistant to anti-arrhythmic drugs and ablation therapies.

  Study phs000362

• Liquid Biopsy Versus Tissue Biopsy to Assess Acquired Resistance and Tumor Heterogeneity in Gastrointestinal Cancers

  Tumor heterogeneity can drive the evolution of multiple tumor subclones harboring unique molecular resistance alterations in different metastatic lesions in an individual patient, under the selective pressure of therapy1-3. Studies have suggested that liquid biopsy may better capture the heterogeneity of acquired resistance, but systematic, direct comparisons of post-progression liquid vs. standard single-lesion tumor biopsies are lacking. In a prospective cohort of 44 patients with molecularly-defined gastrointestinal cancers and acquired resistance to targeted therapy, direct comparison of post-progression liquid vs. tumor biopsy revealed that liquid biopsy more frequently identified clinically-relevant resistance alterations and multiple resistance mechanisms, detecting resistance alterations not detected in matched tumor biopsy in 78% of cases. Whole-exome sequencing of serial cell-free DNA, tumor biopsies and rapid autopsy specimens elucidated geographic and evolutionary characteristics of heterogeneity. Our data suggest that acquired resistance is frequently characterized by profound tumor heterogeneity and highlight the potential clinical utility of liquid biopsy. Burrell, R. A. &#38; Swanton, C. Tumour heterogeneity and the evolution of polyclonal drug resistance. Molecular oncology 8, 1095-1111, doi:10.1016/j.molonc.2014.06.005 (2014). Misale, S. et al. Emergence of KRAS mutations and acquired resistance to anti-EGFR therapy in colorectal cancer. Nature 486, 532-536, doi:10.1038/nature11156 (2012). Ahronian, L. G. et al. Clinical Acquired Resistance to RAF Inhibitor Combinations in BRAF-Mutant Colorectal Cancer through MAPK Pathway Alterations. Cancer Disco 5, 358-367, doi:10.1158/2159-8290.Cd-14-1518 (2015).

  Study phs001853

• Whole Genome Bisulfite Sequencing of Circulating Cell-Free (CCF) DNA and its Cellular Contributors

  Background: Circulating cell free (ccf) fetal DNA has enabled non-invasive prenatal fetal aneuploidy testing without direct discrimination of the genetically distinct maternal and fetal DNA. Current testing may be improved by specifically enriching the sample material for fetal DNA. DNA methylation may allow for such a separation of DNA and thus support additional clinical opportunities; however, this depends on knowledge of the methylomes of ccf DNA and its cellular contributors. Results: Whole genome bisulfite sequencing was performed on a set of unmatched samples including ccf DNA from 8 non-pregnant (NP) and 7 pregnant female donors and genomic DNA from 7 buffy coat and 5 placenta samples. We found CpG cytosines within longer fragments were more likely to be methylated, linking DNA methylation and fragment size in ccf DNA. Comparison of the methylomes of placenta and NP ccf DNA revealed many of the 51,259 identified differentially methylated regions (DMRs) were located in domains exhibiting consistent placenta hypomethylation across millions of consecutive bases, regions we termed placenta hypomethylated domains. DMRs identified when comparing placenta to NP ccf DNA were recapitulated in pregnant ccf DNA, confirming the ability to detect differential methylation in ccf DNA mixtures. Conclusions: We generated methylome maps for four sample types at single base resolution, identified a link between DNA methylation and fragment length in ccf DNA, identified DMRs between sample groups, and uncovered the presence of megabase-size placenta hypomethylated domains. Furthermore, we anticipate these results to provide a foundation to which future studies using discriminatory DNA methylation may be compared.

  Study phs000846

• Stand Up To Cancer East Coast Prostate Cancer Research Group

  Most prostate cancer deaths are caused by metastatic, castration resistant disease (mCRPC). To develop a precision medicine framework for mCRPC, we established a multi-institutional, international clinical sequencing infrastructure to enroll and carry out prospective whole exome and transcriptome sequencing of tumors from a cohort of mCRPC patients. We obtained high quality DNA and RNA sequence data from 150 bone or soft tissue biopsies. Central pathology revealed high-grade adenocarcinoma with only four cases (3.6%) showing neuroendocrine differentiation. Aberrations of AR, ETS genes, TP53 and PTEN were frequent (40-60% of cases), with TP53 and AR alterations being the most enriched in mCRPC compared to primary prostate cancer. We identified novel genomic alterations in PIK3CA/B (fusions and mutations); R-spondin, BRAF and RAF1 (fusions); APC (inactivating mutations); delta-catenin (missense mutations); and ZBTB16/PLZF (homozygous deletions). Aberrations of BRCA2, BRCA1 and ATM were observed at substantially higher frequencies (19.3% overall) than seen in primary prostate cancers, with 56% of these being exclusively somatic. Putative driver gene alterations were identified in nearly all patients, and over half also harbored driver gene fusions, homozygous deletions and/or amplifications. Moreover, 89% of patients harbored a clinically actionable aberration including 62.7% with aberrations in AR, 65% in other cancer-related genes, and 8% with actionable pathogenic germline alterations. Overall, integrative clinical sequencing analysis can be safely and efficiently performed in mCRPC, yields findings that may be actionable for enrolling patients in clinical trials of targeted therapies, and may inform the basis of individual clinical responses.

  Study phs000915

• Molecular Evolution of Cancer

  Recently developed methods that utilize partitioning of long genomic DNA fragments, and barcoding of shorter fragments derived from them, have succeeded in retaining long-range information in short sequencing reads. These so-called read cloud approaches represent a powerful, accurate, and cost-effective alternative to single-molecule long-read sequencing. We developed software, GROC-SVs, that takes advantage of read clouds for structural variant detection and assembly. We apply the method to two 10x Genomics data sets, one chromothriptic sarcoma with several spatially separated samples, and one breast cancer cell line, all Illumina-sequenced to high coverage. Comparison to short-fragment data from the same samples, and validation by mate-pair data from a subset of the sarcoma samples, demonstrate substantial improvement in specificity of breakpoint detection compared to short-fragment sequencing, at comparable sensitivity, and vice versa. The embedded long-range information also facilitates sequence assembly of a large fraction of the breakpoints; importantly, consecutive breakpoints that are closer than the average length of the input DNA molecules can be assembled together and their order and arrangement reconstructed, with some events exhibiting remarkable complexity. These features facilitated an analysis of the structural evolution of the sarcoma. In the chromothripsis, rearrangements occurred before copy number amplifications, and using the phylogenetic tree built from point mutation data, we show that single nucleotide variants and structural variants are not correlated. We predict significant future advances in structural variant science using 10x data analyzed with GROC-SVs and other read cloud-specific methods.

  Study phs001255

• Breast Cancer Family Registry

  The Breast Cancer Family Registry (BCFR) is a multi-center prospective cohort, comprised of over 30,000 women and men from nearly 12,000 families from the United States, Canada, and Australia. Our BCFR resource has been used, and continues to be used, by breast cancer researchers around the world in order to find new ways to prevent, diagnose, and treat cancer. The BCFR provides an extensive and diverse range of resources, expertise, and specialized skills, and has several unique strengths: 1) the collection of a large number of individuals and families across a wide spectrum of breast cancer risk, including both affected and unaffected individuals; 2) the large collection of families with early-onset breast cancer; 3) the large collection of racial/ethnic minority families not replicated elsewhere; 4) the extensive molecular characterization performed to date; and 5) active follow-up of both probands and family members. Thus, the BCFR comprises a unique cohort of probands and family members at familial/genetic risk of breast cancer that will continue to facilitate a wide range of research studies, such as gene discovery, examination of cancer-related outcomes and risk factors in high-risk subjects, investigation of novel behavioral interventions, and cancer prevention trials among at-risk family members. Consequently, the BCFR Cohort, as one of the few cohorts available worldwide with biospecimens and extensive molecular and genetic characterization, combined with epidemiologic data and long-term follow up, will be an invaluable resource for translational research in the genetic epidemiology of breast cancer.

  Study phs002835

• Genome, epigenome and RNA sequences of monozygotic twins discordant for multiple sclerosis

  Monozygotic (MZ) twins have been widely employed for dissection of the relative contributions of genetics and environment in disease. In multiple sclerosis (MS), an autoimmune demyelinating disease that commonly causes neurodegeneration and disability in young adults, disease discordance in MZ twins has been interpreted to indicate environmental importance in pathogenesis. However, genetic and epigenetic differences between MZ twins have been described, challenging the accepted experimental paradigm in disambiguating effects of nature and nurture. Here, we report the genome sequences of one MS-discordant MZ twin pair and messenger RNA (mRNA) transcriptome and epigenome sequences of CD4+ lymphocytes from three MS-discordant, MZ twin pairs. No reproducible differences were detected between co-twins among ~3.6 million single nucleotide polymorphisms (SNPs) or ~0.2 million insertion-deletion polymorphisms (indels). Nor were any reproducible differences observed between siblings of the three twin pairs in HLA haplotypes, confirmed MS-susceptibility SNPs, copy number variations, mRNA and genomic SNP and indel genotypes, or expression of ~19,000 genes in CD4+ T cells. Only two to 176 differences in methylation of ~2 million CpG sites were detected between siblings of the three twin pairs, in contrast to ~800 differences in methylation between T cells of unrelated individuals and several thousand differences between tissues or normal and cancerous tissues. In the first systematic effort to estimate sequence variation among MZ co-twins, we did not find evidence for genetic, epigenetic or transcriptome differences that explained disease discordance. These are the first female, twin and autoimmune disease genome sequences reported.

  Study phs000239

• eMERGE: Genetics of Complex Pediatric Disorders from the Center for Applied Genomics

  The Center for Applied Genomics (CAG) is a specialized Center of Emphasis at the Children&#39;s Hospital of Philadelphia (CHOP) with the primary goal of translating basic research findings to medical innovations. The mission of CAG is to develop new and better ways to diagnose and treat children affected by rare and complex medical disorders. We aim to discover genetic causes for the most prevalent diseases of childhood including asthma, autism, diabetes, epilepsy, obesity, schizophrenia, pediatric cancer, and a range of rare diseases. Ultimately, our objective is to generate new diagnostic tests and to guide physicians to the most appropriate therapies. The CAG is one of the world&#39;s largest genetics research programs in pediatrics, and the only center at a pediatric hospital to have established a large-scale biobank of genotyped samples. The electronic Medical Records and Genomics (eMERGE) Network is a consortium of 9 clinical sites with EMR linked DNA biobanks, including Northwestern University and its NUgene biobank, funded by the NHGRI (National Human Genome Research Institute) to investigate the use of electronic medical record systems for genomic research. The goal of eMERGE is to conduct genome-wide association studies in approximately 100,000 individuals using EMR-derived phenotypes and DNA from linked biorepositories. Using electronic phenotyping methods, the consortium has been and is using DNA samples from all participating sites to explore the genetic determinants of approximately 80 phenotypes, including both diseases and traits, for which the electronic phenotyping algorithms have or are being published on PheKB.org.

  Study phs001165

• Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations

  It is well established that autism spectrum disorders (ASD) have a strong genetic component; however, for at least 70% of cases, the underlying genetic cause is unknown. Under the hypothesis that de novo mutations underlie a substantial fraction of the risk for developing ASD in families with no previous history of ASD or related phenotypes-so-called sporadic or simplex families-we sequenced all coding regions of the genome (the exome) for parent-child trios exhibiting sporadic ASD, including 189 new trios and 20 that were previously reported. Additionally, we also sequenced the exomes of 50 unaffected siblings corresponding to these new (n = 31) and previously reported trios (n = 19), for a total of 617 individual exomes from 209 families deposited in dbGaP. Here we show that de novo point mutations are overwhelmingly paternal in origin (4:1 bias) and positively correlated with paternal age, consistent with the modest increased risk for children of older fathers to develop ASD. Moreover, 39% (49 of 126) of the most severe or disruptive de novo mutations map to a highly interconnected beta-catenin/chromatin remodelling protein network ranked significantly for autism candidate genes. In proband exomes, recurrent protein-altering mutations were observed in two genes: CHD8 and NTNG1. Mutation screening of six candidate genes in 1,703 ASD probands identified additional de novo, protein-altering mutations in GRIN2B, LAMC3 and SCN1A. Combined with copy number variant (CNV) data, these results indicate extreme locus heterogeneity but also provide a target for future discovery, diagnostics and therapeutics.

  Study phs000482

• National Heart Lung and Blood Institute Exome sequencing in SCID

  Analysis of the molecular etiologies of severe combined immunodeficiency (SCID) has led to important insights into the control of immune cell development. Most cases of SCID result from either X-linked or autosomal recessive inheritance of mutations in a known causative gene. However, in some cases, the molecular etiology remains unclear. To identify the cause of SCID in a patient known to lack the protein tyrosine phosphatase CD45, we utilized single nucleotide polymorphisms (SNP) arrays and whole exome sequencing. The patient&#39;s mother was heterozygous for an inactivating mutation in CD45, while the paternal alleles lacked mutations. The patient exhibited a single CD45 mutation identical to the maternal allele. Patient SNP array analysis revealed no change in copy number but loss of heterozygosity for the entire length of chromosome 1 (Chr1), indicating that disease was caused by uniparental disomy (UPD) with isodisomy of the entire maternal Chr1 bearing the CD45 mutation. Non-lymphoid blood cells and other mesoderm and ectoderm-derived tissues retained UPD of the entire maternal Chr1 in this patient who had undergone successful bone marrow transplantation. Exome sequencing revealed mutations in 7 additional genes bearing nonsynonymous SNPs predicted to have deleterious effects. These findings represent the first reported case of SCID caused by UPD and suggest UPD should be considered in SCID and other recessive disorders, especially when the patient appears homozygous for an abnormal gene found in only one parent. Evaluation for alterations in other genes affected by UPD should also be considered in such cases.

  Study phs000479

• Transposable Elements in FTLD-TDP and ALS-TDP

  Amyotrophic Lateral Sclerosis (ALS) is a devastating neurodegenerative disease involving the rapid and progressive loss of motor neurons. While a number of inherited mutations have been identified as causing ALS, the vast majority of cases are sporadic with no family history of disease. This has led to the model that ALS may be a complex syndrome with an unknown number of causal molecular pathways. However, nearly all ALS patient tissues display aggregates of the RNA binding protein TDP-43 (TARDBP/TAR DNA-binding protein), and mutations in the gene encoding this protein are known to cause ALS. We show that one of the normal functions of TDP-43 is to bind and regulate retrotransposons, repetitive sequences in the genome that are normally silent, due in part to the function of TDP-43. Furthermore, we demonstrate that many ALS patients with TDP-43 pathology also display de-silencing of retrotransposon RNAs, consistent with a loss of normal TDP-43 function. Those ALS patients that exhibit high levels of retrotransposon expression form a distinct subset of ALS subjects that also show other correlated alterations to the transcriptome. Using unsupervised machine learning algorithms, we identified three distinct molecular subtypes from the whole genomes and transcriptomes of a large cohort of 148 ALS patient samples. The retrotransposon re-activated subset formed 20% of ALS samples; two additional groups displayed signatures of oxidative stress (61%), and glial dysfunction (19%), respectively. Together these results demonstrate that retrotransposon re-activation is common in a subset of ALS patients with TDP-43 pathology.

  Study phs001889

• Decoding Human Heart Morphogenesis through Single-cell Multi-modal Analyses

  The heart, which is the first organ to develop, is highly dependent on its form to function. However, how diverse cardiac cell types spatially coordinate to create complex morphological structures critical for heart function remains to be elucidated. Here, we show that integration of single cell RNA-sequencing with high-resolution multiplexed error-robust fluorescent in situ hybridization (MERFISH) not only resolves the identity of cardiac cell types developing the human heart but also provides a spatial mapping of individual cells that enables illumination of their organization into cellular communities forming distinct cardiac structures. We discovered that many of these cardiac cell types further specified into subpopulations exclusive to specific communities, supporting their specialization according to cellular ecosystem and anatomic region. In particular, ventricular cardiomyocyte subpopulations displayed an unexpected complex laminar organization across the ventricular wall and formed, with other cell subpopulations, several cellular communities. Interrogating cell-cell interactions within these communities revealed signaling pathways orchestrating the spatial organization of cardiac cell subpopulations during ventricular wall morphogenesis. In vivo conditional genetic mouse models and in vitro human pluripotent stem cell studies confirmed an intricate multicellular PLXN-SEMA crosstalk among specific ventricular cardiomyocyte, fibroblast and endothelial cell subpopulations that directs the compaction of the ventricular wall layers. Thus, these detailed findings into the cellular social interactions and specialization of cardiac cell types constructing and remodeling the human heart offer new insights into structural heart diseases as well as engineering complex multi-cellular tissues for human heart repair. Reprinted from Farah et al. (In press), with permission from Nature.

  Study phs002031

• Non-Coding Autoimmune Risk Variant Defines Role for ICOS in T Peripheral Helper Cell Development

  Our fine-mapping of 76 non-MHC loci associated with risk for rheumatoid arthritis (RA, 11,475 cases, 15,870 controls) has recently identified rs117701653, a non-coding single nucleotide polymorphism (SNP) in the CTLA4/CD28/ICOS locus, as the variant most likely to modulate RA risk within this region, with the minor allele C showing disease protection compared to the major allele A. We found that this SNP exhibits allele-specific protein binding and transcriptional modulation, further supporting its regulatory nature. Here, we aimed to determine whether the functional non-coding variant rs117701653 modulates the expression of genes within the CD28/CTLA4/ICOS locus. To address this question, we employed low-input RNA sequencing using the Illumina SmartSeq2 platform on resting total CD4+ T cells from the blood of 24 healthy subjects (8 AA, 8 AC, 8 CC) recruited based on genotype from the Mass General Brigham biobank. We performed a targeted cis-eQTL mapping analysis for protein coding genes with transcription start sites (TSS) within a 1MB window of SNP rs117701653. These genes included WDR12, NBEAL1, CYP20A1, ABI2, RAPH1, CD28, CTLA4, ICOS, and PARD3B. Surprisingly, allelic variation at rs117701653 correlated strongly with expression of ICOS. This suggests that individuals carrying the protective allele C at rs117701653 express lower levels of ICOS, implicating ICOS as a plausible target gene influenced by this protective allele. The raw sequence data from the 24 healthy individuals are available through dbGAP. Gene expression matrices have been deposited to GEO accession GSE235868.

  Study phs003448

• ILyAD (Indolent Lymphoma And vitamin D)

  Indolent non-Hodgkin lymphomas are incurable diseases, and require intermittent and often morbid and expensive therapy during their prolonged natural history. Lower intensity and better-tolerated, cost-effective treatment strategies are needed for these patients. Low vitamin D levels at diagnosis of indolent lymphomas are strongly associated with inferior outcomes to treatment. Standard therapy for low tumor burden, indolent lymphoma includes a single agent monoclonal antibody (rituximab). In our double-blind, placebo-controlled, randomized trial we evaluated if vitamin D supplementation with 2000IU oral cholecalciferol daily for three years improved the progression-free survival of patients with indolent lymphoma treated with single agent rituximab. Stratification factors for randomization included histology and Follicular Lymphoma International Prognostic Index (FL-IPI) prognostic score. The primary endpoint of the trial was progression-free survival over three years; secondary endpoints included response at 13 weeks, and overall survival. Eligible participants were adults with biopsy proven follicular lymphoma (grades 1-3a), small lymphocytic lymphoma, marginal zone lymphoma and mucosal-associated lymphoid tissue with no prior anti-lymphoma systemic therapy. Participants had Ann Arbor stage II or greater, measurable disease defined by Lugano criteria. Participants with osteoporosis requiring prescription treatment, symptomatic primary hyperthyroidism, hypercalcemia, creatinine greater than 1.5 times the upper limit of normal or a history of calcium-related kidney stones were excluded. Germline DNA was isolated from saliva samples obtained from participants at baseline. Whole exome sequencing was performed and SNPs were analyzed using the Agilent SureSelectXT Human All Exon V6 +UTR array. SNP data from 185 participants are available through dbGaP.

  Study phs003503

• Cerebrospinal fluid cfDNA sequencing for classification of central nervous system glioma

  Primary central nervous system (CNS) gliomas can be classified by characteristic genetic alterations. In addition to solid tissue obtained by surgery or biopsy, cell-free DNA (cfDNA) derived from cerebrospinal fluid (CSF) is an alternative source of material for genomic analyses. Experimental design: We performed targeted next-generation sequencing (NGS) of CSF cfDNA in a representative cohort of 85 patients presenting with with suspicion of primary or recurrent glioma at two neurooncological centers. Copy-number variation (CNV) profiles, single nucleotide variants (SNVs), and small insertions/ deletions (indels) were combined into a molecular-guided tumor classification. Comparison with the solid tumor was performed for 38 cases with matching solid tissue available. Results: Cases were stratified into four groups: glioblastoma (n = 32), other glioma (n = 19), non-malignant (n = 17) and no confirmed diagnosis available (n = 17). We introduce a molecular-guided tumor classification, enabling identification of tumor entities and/ or cancer specific alterations in 75.0 % (n = 24) of glioblastoma and 52.6 % (n = 10) of other glioma cases. The overlap between CSF and matching solid tissue was highest for CNVs (26-48 %) and SNVs at pre-defined gene loci (44 %), followed by SNVs/indels identified via uninformed variant calling (8-14 %). A molecular-guided tumor classification was made 23.5 % (n = 4) of cases with no confirmed diagnosis available. Conclusions: We developed a workflow for targeted sequencing of CSF cfDNA as well as a strategy for interpretation and reporting of sequencing results based on a molecular-guided tumor classification in glioma.

  Study EGAS50000000060

• Lactate metabolism in cancer stem cell fate regulation

  Tumors arise as a result of uncontrolled cell proliferation driven by mutations in the signaling pathways that regulate the balance between cell division and differentiation. Despite of that, in the tumor, cancer stem cells seem to retain their capability to differentiate, which can lead to differential chemotherapy response. We previously showed that intestinal stem cells and differentiated cells display specific metabolic features and that they interact through metabolites. Moreover, we found that altering mitochondrial metabolism is sufficient to lead and define stem cell differentiation. In cancer, cellular metabolism is highly derailed. However, how metabolic changes in the cancer cell contribute to cell fate decisions remains unknown. Recently, the crosstalk between metabolism and epigenetics has emerged as one of the players in tumor development. Here we aim to address the impact of the altered metabolism of cancer cells on their epigenome as a regulatory axis of cell identity and fate. We use human-derived tumor organoids where we introduce genetic encoded fluorescent-based reporters to visualize cell identity and metabolism at single cell level. We combine 4D live imaging with a recently developed machine learning method to re-construct cell-lineages, trace stem cell state and metabolic changes during tumor organoid development. We find that, mimicking the tumor metabolic microenvironment (Warburg effect-driven low-glucose high-lactate) leads to changes cellular metabolism and histone acetylation. Strikingly, these alterations provoke profound changes in tumor development by increasing the population of cancer stem cells and their proliferation and differentiation dynamics.

  Study EGAS50000000063

• NCI Cancer Model Development for the Human Cancer Model Initiative (HCMI)

  The Cancer Model Development Center (CMDC) project is part of the Human Cancer Model Initiative (HCMI) and is managed by the Office of Cancer Genomics within the Center for Cancer Genomics at the National Cancer Institute (NCI). The goal of the HCMI is to generate approximately 1000 novel human "next generation" tumor-derived models as a community resource. The models aim to faithfully recapitulate the cellular complexity and heterogeneity of human tumors. The CMDCs are tasked by the NCI with producing a proportion of the cancer models. Wellcome Sanger Institute (WSI) is a contributing member of the HCMI. Models developed by HCMI are being derived from a number of tumor types and subtypes (see below), including rare adult, pediatric and understudied cancers. The number of models and the clinical sites from which patient samples are acquired is increasing. Normal tissue, originating tumor, and tumor-derived models are sequenced. The next generation cancer models are associated with rich clinical and genomic data through the NCI Genomic Data Commons for all researchers to advance cancer research. The CMDC and HCMI models will serve as excellent tools in research to provide insight into the pathways that influence tumor growth and progression, how tumors respond to therapy, etc. Ultimately, the data generated using these models will support the development of personalized oncology. There is a delay between registered IDs and model availability. The current model count is found in The HCMI Searchable Catalog which is updated when new models are available for distribution to the research community.

  Study phs001486

• Single Cell and Spatial Transcriptomic Profiling of Haemophilus ducreyi Infection

  We investigated the immune response in end-point pustules formed by Haemophilus ducreyi in volunteers using an experimental infection model. Six volunteers were infected in the skin of the upper arm at three sites with 41-85 CFU of H. ducreyi loaded on an allergy testing device; a sham-inoculated (e.g., wound) site was also performed. Five volunteers (4 men and 1 woman; 2 Asians and 3 Whites; aged 27-46 years) formed pustules 6-8 days post infection; 1 volunteer spontaneously resolved all sites. If a volunteer formed pustules, biopsies were performed to collect their pustule and wound sites. Biopsy pairs for all 5 pustule formers were subjected to single cell (sc) RNA-seq. Sequencing libraries were prepared using the 10X Genomics Single-cell 3' RNA-seq kit. Spatial transcriptomics was performed on biopsy pairs from 4 of the 5 pustule formers. Ten micron sections were placed on 10X Genomics Visium slides for fresh-frozen tissue, and libraries were prepared with the 10X Genomics Visium spatial 3' gene expression kit. All samples were sequenced on an Illumina NovaSeq 6000. Count tables were produced with CellRanger and SpaceRanger, respectively, and Seurat was used to analyze the data. We identified 13 major cell types by scRNA-seq; the proportion of some immune cell subsets was increased in pustules compared to wound sites. We also localized the major cell types within the tissue by using the scRNA-seq data to deconvolve the spatial transcriptomic data. The FASTQ files and count tables from CellRanger and SpaceRanger are provided.

  Study phs003754

• Single Chromatin Fiber Profiling in the Human Brain

  We developed a human postmortem brain-applicable &#8216;Fiber-seq' protocol for genome-scale, PCR-free single molecule sequencing of on average 8-10kb long chromatin fibers, to map nucleosome positioning and nucleosome-depleted regions which include, for example, active promoters and enhancers. We generated two Fiber-seq reference sets, one from sorted non-neuronal NeuN- prefrontal cortex (PFC) nuclei from a 32 year old female donor, with total of 11,143,795 reads/fibers, and median of 30-fold genomic coverage. The other reference set was generated on sorted prefrontal NeuN+ neurons, from two 24 year old brain donors (1 female/1 male), with a combined total read/fiber count of 4,749,556 and median of 20-fold genomic coverage. Each of the peaks in the Fiber-seq libraries seq signal correlated with the ATAC peak levels in the neuronal and non-neuronal ATAC-seq libraries, respectively (promoters, clusters 1, 3; NeuN+, r = 0.574; NeuN-, r=0.683, PHowever, Fiber-seq detected 19,978 highly significant peaks which, while showing some non-significant/near-background ATAC-Seq alignments, lacked a corresponding site-specific enrichment in conventional short-read ATAC- and ChIP-seq datasets. However, 60% of these peaks were in repetitive DNA elements. Therefore, long read single chromatin fiber sequencing, which exceed the base pair length of conventional short-read libraries by two orders of magnitude, could capture many actively regulated repeat elements at the site of specific gene loci, that otherwise would remain unmappable by conventional ATAC-, and histone- ChiP-seq.

  Study phs003771

• Multiple Myeloma Clinical Targeted Sequencing of Patient Samples

  Purpose: Molecular diagnostics for the risk stratification of patients with multiple myeloma (MM) have languished in the area of cytogenetics and fluorescence in situ hybidization (FISH). However, here we present the case for next generation sequencing of MM patient samples to not only identify high-risk markers, but also mutations and deletions relating to immunomodulatory drugs (IMiDs), and more importantly to guide the sequencing of immunotherapy regimens in response to intrinsic antigen escape as a means of treatment resistance and relapse. Experimental design: CD138+ cell selected samples (n=134) from bone marrow aspirates from patients with multiple myeloma were sequenced using a targeted panel in a clinical diagnostics laboratory. Data were analyzed for high-risk markers, treatment-related resistance mechanisms, and precision medicine targets to guide the future treatment of patients in the clinic. Results: High-risk markers including t(4;14), t(14;16), t(14;20), gain or amplification 1q, deletion of CDKN2C, and deletion or mutation of TP53 were identified in 15% of newly diagnosed myeloma patient samples. At relapse, alterations in the cereblon degradation pathway were found in 24.3% of IMiD treated patients. Deletions of 4p (CD38) were also enriched in patients who received anti-CD38 treatment (P=0.03) which were mostly monoallelic. Deletions and mutations were detected in TNFRSF17 encoding BCMA in patients treated with anti-BCMA regimens, and the information was used to change the treatment of the patients. Conclusions: Targeted sequencing of multiple myeloma patient diagnostic samples can be used for risk stratification, and also to monitor and adjust treatments as resistance mechanisms evolve.

  Study phs003908

• Circulating Tumor DNA Analysis in ERBB2-Amplified Colorectal Cancer: Biomarker Analysis of the MyPathway Trial

  Combination of two HER2-directed antibodies, pertuzumab and trastuzumab (P+T), has antitumor activity in HER2-positive colorectal cancer (CRC). Although liquid biopsies are increasingly being used in clinical oncology, the association between tumor and circulating tumor DNA (ctDNA) ERBB2 status and ctDNA monitoring for early response and resistance are unknown. Eighty-five patients with ERBB2-amplified and/or overexpressed CRC were treated with P+T in the MyPathway trial; 42 had ctDNA testing at Cycle 1 Day 1 (C1D1), and 38 had longitudinal plasma tested for ctDNA. We analyzed the ctDNA vs tissue ERBB2 concordance, genomic co-alterations, and ctDNA dynamics and association with response. Forty-one (98%) of 42 patients had genomic alterations detected in ctDNA at C1D1, and 29 (69%) had ERBB2 amplification in ctDNA. There was a strong correlation between ERBB2 copy number on next-generation sequencing in tissue and C1D1 ERBB2 ctDNA copy number. Thirty-seven percent achieved a molecular response by C3D1 on P+T, which was associated with prolonged progression-free survival and overall survival (OS). CDKN2A and KRAS mutations were associated with shorter OS, and a trend was seen with PIK3CA mutations. Several emerging co-alterations were identified in ctDNA at progression, including in the MAPK and PI3K pathways and other tyrosine receptor kinases. ctDNA can detect ERBB2 amplification in many, but not all, patients with ERBB2 amplification detected in tumor samples. ctDNA molecular response was associated with better survival, and ctDNA co-alterations may offer insights into mechanisms of intrinsic and acquired resistance.

  Study EGAS50000000916

• PCR-free shallow whole genome sequencing for chromosomal copy number detection from plasma of cancer patients is an efficient alternative to the conventional PCR-based approach

  Somatic copy number alterations (SCNAs) can be detected in tumor cell-free DNA (cfDNA) by shallow whole genome sequencing (sWGS). Polymerase chain reaction (PCR) is typically included in library preparations, but a PCR-free method could be an effective alternative for research and diagnostics. Peripheral blood samples (n=22) were collected in EDTA-containing tubes from non-small cell lung cancer patients (n=10) and healthy donors (n=12). PCR and PCR-free library preparations underwent sWGS and the following metrics were compared: mapping quality, proportion of duplicate reads, genome coverage, sequencing read length and copy number profiles. The PCR-free method was further evaluated on archived plasma samples (n=33) of varying pre-analytical quality collected in EDTA or lithium heparin tubes. Archived samples were from acute myeloid leukaemia (n=15) and B-cell lymphoma (n=18) patients. The percentage of unique reads was significantly higher for the 22 samples analysed by PCR-free (95.7%) compared to PCR (84%). PCR-free libraries from 33 archived samples produced data of a comparable quality, with 95.4% unique reads. All other evaluated metrics were highly comparable for PCR and PCR-free library preparations. We conclude that SCNA detection can be performed by PCR-free sWGS on cfDNA from cancer patients and is feasible for diagnostic and archival plasma collected in EDTA- or lithium heparin-containing tubes. Our results pave the way for an automated cfDNA workflow for cancer plasma samples.

  Study EGAS00001004692

• Coding and regulatory somatic profiling of triple-negative breast cancer in Sub-Saharan African patients

  The burden of triple-negative breast cancer (TNBC) may be shaped by genetic factors, particularly inherited and somatic mutation profiles. However, data on this topic remain limited, especially for the African continent, where a higher TNBC incidence is observed. In the age of precision medicine, cataloguing TNBC diversity in African patients becomes imperative. We performed whole exome sequencing, including untranslated regions, on 30 samples from Angola and Cape Verde, which allowed to ascertain on potential regulatory mutations in TNBC for the first time. A high somatic burden was observed for the African cohort, with 86% of variants being so far unreported. Recurring to predictive functional algorithms, 17% of the somatic single nucleotide variants were predicted to be deleterious at the protein level, and 20% overlapped with candidate cis-regulatory elements controlling gene expression. Several of these somatic functionally-impactful mutations and copy number variation (mainly in 1q, 8q, 6 and 10p) occur in known BC- and all cancer-driver genes, enriched for several cancer mechanisms, including response to radiation and related DNA repair mechanisms. TP53 is the top of these known BC-driver genes, but our results identified possible novel TNBC driver genes that may play a main role in the African context, as TTN, CEACAM7, DEFB132, COPZ2 and GAS1. These findings emphasize the need to expand cancer omics screenings across the African continent, the region of the globe with highest genomic diversity, accelerating the discovery of new somatic mutations and cancer-related pathways.

  Study EGAS50000001013

• Epigenetic Control of Topoisomerase 1 by MacroH2A1.1

  DNA transactions introduce torsional constraints that pose an inherent risk to genome integrity. While topoisomerase 1 (TOP1) activity is essential for removing DNA supercoiling, aberrant stabilization of TOP1:DNA cleavage complexes (TOP1ccs) can result in cytotoxic single-stranded DNA lesions (SSLs). What protects the genome from aberrant TOP1 activity remains unknown. Using a combination of CUT&#38;RUN and TOP1 Covalent Adduct Detection (CAD) Seq in MDA-MB-231 breast cancer cells (American Type Culture Collection, ATCC), we identified chromatin context as an essential means to ensure TOP1cc resolution at TOP1 hot spots. Through its ability to bind poly(ADP-ribose) (PAR), a protein modification required for TOP1cc repair, the histone variant macroH2A1.1 establishes a TOP1-permissive chromatin environment, while the alternatively spliced macroH2A1.2 isoform is unable to bind PAR or protect from TOP1ccs. MacroH2A1 isoform-specific analyses were based on reconstitution of macroH2A1.1 knockout cells with wild-type or PAR binding-deficient FLAG-tagged macroH2A1.1 in human breast cancer cells. Mechanistically, we find that macroH2A1.1 facilitates the recruitment of the TOP1cc repair factor XRCC1 in response to both endogenous and drug-induced topological stress. Impaired macroH2A1.1 splicing, a frequent cancer feature, was predictive of increased sensitivity to TOP1 poisons in a pharmaco-genomic screen in breast cancer cells, and macroH2A1.1 inactivation mirrored this effect in breast and ovarian cancer cells. Consistent with this, low macroH2A1.1 expression correlated with improved survival in cancer patients treated with TOP1 inhibitors (TOP1i). We propose that macroH2A1 alternative splicing serves as an epigenetic modulator of TOP1-associated genome maintenance and a potential cancer vulnerability.

  Study phs003729

• Genomic Characterization of Duke Melanoma Brain Metastases

  Previous work has shown that melanoma brain metastases (MBM) have a unique molecular profile compared to extracranial metastases (ECM). Description of the biology of MBM and associations between biological features and clinical outcomes will facilitate the design of rational therapies for patients with MBM. To better characterize the molecular profile of MBM, we conducted whole exome sequencing (WES) and total RNA sequencing on MBM tissues obtained from an independent cohort of 14 patients from Duke University Medical Center. This cohort adds valuable patient samples to the field and enables strengthened investigation of somatic mutations and gene expression profiles within and between cohorts. The mutational landscape and gene expression of MBM from the Duke cohort resembled those previously reported in a previously published dataset of WES and RNA sequencing in MBM and matched ECM tissues from UT MD Anderson Cancer Center (Fischer et al., 2019 30787016). Duke MBM demonstrated similar rates of mutations in genes associated with melanoma, including BRAF, NRAS, and PTEN, as the MD Anderson cohort. Immune populations in MBM were similar between Duke and MD Anderson MBM, though Duke samples demonstrated more resting CD4+ memory T cells. The ratio of M2 to M1 macrophages in the Duke cohort was similar to that of the MD Anderson cohort. In MBM, a higher M2:M1 ratio may contribute to an immune-suppressive tumor microenvironment. Sequencing data and sample attributes from MBM tissues obtained from 14 patients from Duke University Medical Center will be available through dbGaP.

  Study phs003009

• Oncogenome of Kaposi Sarcoma

  Kaposi Sarcoma (KS) is an aggressive cancer caused by the Kaposi Sarcoma Herpesvirus (KSHV/HHV-8). Individuals with immunodeficiencies, exemplified by HIV, have an elevated risk of developing KS. However, understanding of the genetic factors contributing to KS progression is still limited. To explore potential genetic alterations in KS that could offer biological or therapeutic insights, whole exome sequencing was performed on 78 KS tumors and matching normal skin samples from 59 adults with KS (46 with HIV-associated KS and 13 with HIV-negative KS) receiving treatment at the Uganda Cancer Institute in Kampala, Uganda. Results showed a very low mutational burden in all samples except one (median = 11 mutations), which is the smallest number of mutations found across all 33 cancer types in The Cancer Genome Atlas (TCGA). No recurrent mutations were identified. Mutational signatures included impaired DNA mismatch repair and smoking. There was no evidence suggesting that multiple tumors from the same patient originated from a single clone. The number of genome copy alterations per genome was higher in tumors from individuals without HIV and those with advanced-stage disease. This suggests that lesions that take longer to develop may accumulate more alterations, although the overall number of alterations remains low compared to other cancers. These findings suggest that KS pathogenesis differs from other cancers, with KSHV viral infection and the expression of viral oncogenes being the primary drivers of carcinogenesis, rather than clonal oncogenic transformation arising from genetic alterations of cancer-related cellular genes.

  Study phs003897

• Whole-genome and transcriptome versus panel sequencing in precision oncology: A translational-clinical comparison

  Precision oncology offers new cancer treatment options. While the use of whole-genome (WGS) and transcriptome (TS) sequencing remains often confined to specific clinical programs, molecular analysis with gene panels is more common clinical practice. Methods: Twenty patients with rare or advanced tumors sequenced by WGS/TS within the DKFZ/NCT/DKTK-MASTER program from 2015 to 2020 (MASTER 1) were additionally sequenced with the TSO500 DNA and TST170 RNA Assay by Illumina. Original Molecular Tumor Board (MTB) therapy recommendations (TR) were juxtaposed with TR from NGS panel data analysis (comparison 1). To mitigate temporal differences between the two sequencings, WGS/TS was bioinformatically reanalyzed (MASTER 2), and new TR were again compared (comparison 2). Results: In comparison 1, 47.1% of the MASTER 1’s TR and 62.8% of the panel’s TR were identical, as were 45.9% of the TR of MASTER 2 and 54.9% of the panel in comparison 2. In both MASTER 1 and MASTER 2, 26.5% and 36.1% of TR, respectively, relied on Biomarkers (BM) which were not covered by the panel. Eight of the ten MTB-derived therapy implementations were also supported in the panel, the two remaining ones being based on BM absent from the panel. Conclusions: About half of the TR in WGS/TS and panel sequencing were identical. Around a third of TR issued by WGS/TS were based on BM which are not covered by the panel, leading to therapy implementations. The presented data underscores the clinical value of comprehensive molecular analyses in individuals with advanced cancers.

  Study EGAS50000000431

• Singel-cell RNA sequencing and CUT&RUN sequencing of human RUNX2-deficient osteoblasts

  Runx2 is a master regulator of bone formation, and its dysfunction causes cleidocranial dysplasia (CCD) in humans. When iPS cells were generated from patients with CCD and Runx2-deficient iPS cells were generated using gene-editing techniques, abnormal laminopathy-like nuclei were observed. Runx2-deficient cells showed reduced Lamin A/C expression, but not protein levels. However, in Runx2-deficient cells, both the gene expression and protein levels of Nesprin-1 were reduced, perinuclear actin fibers were sparser, and nuclear stiffness was reduced. Forced expression of Lamin A/C increased nuclear stiffness but did not improve nuclear morphology. In contrast, the induction of Nesprin-1 expression alone normalized nuclear stiffness and restored nuclear morphology and perinuclear actin distribution. In Runx2-null cells, mechanical stress-induced phosphorylation of emerin was not observed. In contrast, forced expression of nesprin-1 in Runx2-null cells resulted in phosphorylation of emerin, indicating the restoration of intracellular tension. These observations were confirmed by atomic force microscopy. Therefore, the intracellular tension was inferred to pull the nuclear membrane into its normal shape. CUT&RUN assay and single RNA-seq analysis showed that an aberrant nuclear membrane caused loss of nuclear lamina gene regulation machinery, making the progression of normal osteogenic differentiation impossible; however, supplementation with Nesprin-1 restored gene regulation mechanisms and promoted pre-osteoblast formation with normal nuclear morphology. Nesprin-1 expression induced by Runx2 is essential for epigenetic regulation of the nuclear lamina. We propose CCD as a new type of laminopathy involving defective expression of Nesprin-1 regulated by Runx2.

  Study JGAS000663

• Mutational signatures in normal tissues of patients treated with clinical stage G-quadruplex binder CX5461

  Drug induced photoreactions following exposure to UVA have been documented for several therapeutic compound classes, including psoralens and fluoroquinolones. CX5461, a quinolone-derived small molecule currently in phase 1 studies for cancer treatment, has been shown to be a clinically significant UV photosensitizer. Here we examined the mutagenic spectrum of CX5461 in normal tissues of patients undergoing treatment, by comparing deep whole genome sequencing (WGS) from skin biopsies taken from non-sun exposed areas pre and post drug exposure of patients in the Canadian Cancer Clinical Trials Group (CCTG) IND.231 trial (NCT02719977). We show that a low burden of SNVs (<100) and indels (<40) mutations was detected in post CX5461 skin biopsies, with a pattern resembling background or sequencing artefact signatures. To further characterize the mutational pattern of CX5461 photoreactivity, we compared mutational profiles of human retinal pigment epithelial cells exposed to UVA, CX5461, or a combination of both UV and CX5461. Treatment with CX5461 or UVA alone, resulted in a low SNV burden. In contrast, human cells exposed to UVA and CX5461 had a two to four orders higher mutational burden, with a distinct mutational pattern characterized by T>A, T>C mutations, dissimilar to that of UV alone or CX5461 alone. Our results suggest that human exposure to CX5461 alone is not associated with high levels of single base mutagenesis. However, CX5461 is a very potent UV photosensitizer, underscoring the need for adequate UV screening and controls when studying the effects of photosensitizing small molecules on single base DNA variation.

  Study EGAS50000001145

• Transcriptome profiling of megakaryocytes and platelets: application to GP9- and IKZF5-related thrombocytopenia

  Platelets are anucleate cells produced in the bone marrow and derived from large progenitor cells called megakaryocytes (MKs). Platelets receive RNA transcripts from their progenitorial MKs during thrombopoiesis. However, the correspondence between platelet and MK transcriptomes is poorly understood, particularly in the context of germline mutations that cause platelet formation defects or thrombocytopenia. We have studied the effects of two such mutations on MK and platelet transcriptomes. We generated immortalized MK cell line (imMKCL)-based models of Bernard-Soulier syndrome and IKZF5-related thrombocytopenia. MKs derived from imMKCLs with either a homozygous deletion of GP9 (GP9-/-) or a heterozygous Y121F variant in IKZF5 (IKZF5WT/Y121F) exhibited reduced proplatelet formation (reductions of 96% and 57%, respectively). Platelets from patients with either GP9-/- or IKZF5WT/Y121F genotypes had broad transcriptomic dysregulation, suggesting that (pro)platelet formation defects due to mutations in glycoprotein receptor and transcription factor genes such as GP9 and IKZF5 affect the trafficking of transcripts from MKs into platelets. RNA-seq data from MKs at four stages of differentiation revealed widespread but distinct changes in expression over time between the GP9-/- and the IKZF5WT/Y121F genotypes. Dysregulated genes in GP9-/- MKs were enriched for RNA metabolism and actin/tubulin folding pathways while those in IKZF5WT/Y121F MKs were enriched for cell cycle pathways. Most of these genes were also dysregulated in the platelets of patients with the corresponding diseases. Our results suggest that patients with inherited forms of thrombocytopenia present with a transcript trafficking defect during platelet formation.

  Study EGAS50000001276

• Complex genotype-phenotype relationships shape the response to treatment of Down Syndrome Childhood Acute Lymphoblastic Leukaemia

  Extensive genetic and epigenetic variegation has been demonstrated in many malignancies. Importantly, their interplay has the potential to contribute to disease progression and treatment resistance. To shed light on the complex relationships between these different sources of intra-tumour heterogeneity, we explored their relative contributions to the evolutionary dynamics of Acute Lymphoblastic Leukaemia (ALL) in children with Down syndrome, which has particularly poor prognosis. We quantified the tumour propagating potential of genetically distinct sub-clones using serial transplantation assays and SNP-arrays. While most leukaemias were characterized by a single dominant subclone, others were highly heterogeneous. Importantly, we provide clear and direct evidence that genotypes and phenotypes with functional relevance to leukemic progression and treatment resistance can co-segregate within the disease. Hence, individual genetic lesions can be restricted to well-defined cell immunophenotypes, corresponding to different stages of the leukemic differentiation hierarchy and varied proliferation potentials. As a result of this difference in fitness, which can be accurately quantified via competitive transplantation assays, matching diagnostic, post-treatment, and relapse leukaemias can be dominated by different genotypes, including pre-leukemic clones persisting throughout the disease progression and treatment. Intriguingly, plasticity also appears to be a temporally defined property that can segregate with genotype. These results suggest that Down Syndrome ALL should be viewed as a complex matrix of cells exhibiting genetic and epigenetic heterogeneity that foster extensive clonal evolution and competition. Therapeutic intervention reshapes this ‘eco-system’ and may provide the right conditions for the preferential expansion of selected compartments and subsequently relapse.

  Study EGAS50000001287

• Comprehensive genomic analysis of colorectal cancer with microsatellite instability

  Microsatellite instability-high (MSI-H) colorectal cancers (CRCs) account for 10?15% of all CRCs. MSI-H CRC is characterized by a large number of somatic insertions/deletions (indels) resulting from either mutations within or the silencing of genes involved in the DNA mismatch repair (MMR) system. A subset of MSI-H CRCs is associated with Lynch syndrome (LS), which is caused by germline MMR gene mutations, leading to hereditary cancer predisposition. Other than frequent somatic mutations in BRAF, the transformation mechanisms underlying MSI-H CRC are largely unknown. Here, genomic DNA from 149 MSI-H CRC specimens was analyzed using whole-exome sequencing, and 93 of these samples were subjected to genome-wide DNA methylation analysis. Furthermore, transcriptome sequencing was conducted on 111 samples. Genomic/epigenomic analyses identified three subgroups within our cohort: (1) MSI-H CRCs with silenced MLH1 that share frequent indels, a specific mutation/copy number alteration profile and promoter DNA methylation; (2) LS-associated MSI-H CRCs with MMR genes containing germline mutations; and (3) the remaining MSI-H CRCs with frequent somatic disruptive mutations of MLH1 or MSH2. Unexpectedly, only the first group was found to frequently carry fusion-type protein kinases (15% of this group and 9.9% of all MSI-H CRCs) that are promising therapeutic targets. Thus, MSI-H CRCs can be classified into three subgroups with distinctive genomic as well as epigenomic statuses, probably reflecting the oncogenic processes of these cancers. Fusion-type kinases are enriched in MSI-H CRCs, shedding new light on potential treatment strategies.

  Study JGAS000113

• Frequent post-treatment monitoring of colorectal cancer using individualized ctDNA validated by multi-regional molecular profiling

  Purpose: Circulating tumor DNA (ctDNA) analysis has been proposed as an approach for prediction of post-treatment patient outcomes. However, whether a single platform will provide optimal information in all patients or alternatively a patient-specific monitoring approach based on assessment of mutations in the primary tumor from that patient remains an urgent question. Experimental Design: We conducted multiregional sequencing of 42 specimens of 14 colorectal tumors (Stage III and more) from 12 patients, including two double cancer cases, to identify the full spectrum of mutational heterogeneity and identify aberrations that could be used to develop personalized ctDNA assays. Results: ?Founder? mutations that occur in all regions of the sample were identified in 12/14 (85.7%) tumors. Subsequent phylogenetic analysis of each tumor showed that 12/14 tumors (85.7%) carried at least one ?truncal? mutation. Most founder and truncal mutations exhibited higher variant allele frequency (VAF) than ?non-founder? and ?branch? mutations. In addition, both founder and truncal mutations were more likely to be detected as ctDNA than non-founder and branch mutations. Synchronized ctDNA dynamics of multiple mutations suggested those mutations from the same clonal origin. For 10/12 patients (83.3%) with nearly 1,000 days of post-operative observation, the validity of frequent personalized ctDNA monitoring was confirmed in terms of early relapse prediction, treatment efficacy, and non-relapse corroboration. Conclusions: Personalized ctDNA monitoring based on aberrations with a high VAF in the primary tumor site should be explored in larger prospective clinical trials to determine the full clinical validity.

  Study JGAS000243

• Durable clinical impacts and mechanisms of action and resistance in histone K27 methylation-targeting epigenetic therapy

  Epigenomes allow the rectification of disordered cancer gene expressions, thereby providing new targets for pharmacological interventions in cancer therapy. The clinical utility of targeting histone H3 lysine tri-methylation (H3K27me3) as an epigenetic hallmark has been demonstrated. However, in actual therapeutic settings, the mechanism by which H3K27me3-targeting therapies exert their effects and the response of tumor cells remain unclear. Here we show the potency and mechanisms of action and resistance of a newly developed EZH1/2 dual inhibitor, valemetostat, in the clinical trials of patients with adult T-cell leukemia/lymphoma (ATL). Valemetostat administration reduced the tumor size and demonstrated durable clinical benefits. Integrative single-cell analyses of patients before and after the administration of valemetostat showed that valemetostat abolishes the highly condensed chromatin structure formed by the plastic H3K27me3 and neutralizes multiple gene loci, such as tumor suppressor genes. Nevertheless, subsequent long-term treatment encounters the emergence of resistant clones with reconstructed aggregate chromatin that closely resemble the pre-dose state. The acquired mutations at the EZH2???compound interface result in the propagation of the clones with re-increased H3K27me3. In patients free of this mutation, TET2 mutation or elevated DNMT3A expression causes similar chromatin recondensation through de novo DNA methylation at the H3K27me3-associated regions. We identified subpopulations with distinct metabolic and gene translation characteristics implicated in primary susceptibility until the acquisition of the heritable (epi)mutations. Targeting chromatin homeostasis might provide opportunities for further sustained epigenetic cancer therapies.

  Study JGAS000553

• CAGE-seq analysis of osteoblast derived from cleidocranial dysplasia human induced pluripotent stem cells

  Non-coding RNAs (ncRNAs) comprise a major portion of transcripts and serve an essential role in biological processes. Although the importance of major transcriptomes in osteogenesis has been extensively studied, the function of ncRNAs in human osteogenesis remains unclear. Previously, we developed hiPSCs from patients with cleidocranial dysplasia (CCD) caused by runt-related transcription factor 2 (RUNX2) haploinsufficiency. To gain insight into ncRNAs in osteogenesis, we surveyed differential ncRNA expression profiling and promoter differences of RUNX2 using patient-specific iPSCs and cap analysis gene expression (CAGE) technology to define the promoter landscape. Revertant iPSCs (Rev1 iPSCs) edited by CRISPR/Cas9 system to harbor mutation-corrected RUNX2 exhibited increased proximal promoter expression of RUNX2, while CCD iPSCs did not. We identified 2271 ncRNA genes with altered expression levels before and after differentiation, 31 of which showed at least 20-fold higher expression in Rev1 iPSCs. Bioinformatic analysis also categorized AC007392.3, LINC00379, RP11-122D10.1, and RP11-90J7.2 as enhancer regulatory regions, and HOXA-AS2, MIR219-2, and RP11-834C11.3 as dyadic regulatory regions of these ncRNAs. In addition, two miRNAs, termed MIR199A2 and MIR152, were found to have high enrichment of osteogenic-related terms. Upon further examination of the role of MIR152 on osteoblast differentiation, we found that MIR152 knockdown induced upregulation of ALP and COL1A1 in Saos-2 cells. Thus, ncRNAs were found to regulate the osteogenic differentiation potentials of hiPSCs that are used for bone regeneration and repair owing to their differentiation potentials. These data allow understanding ncRNA profiles of hiPSCs during osteogenesis.

  Study JGAS000248

• Defined Lifestyle and Germline Factors Predispose Asian Populations to Gastric Cancer

  Gastric cancer (GC) is one of the leading causes of cancer mortality throughout the world, with the highest incidence in east-Asia. Although the somatic genetics of GC has extensively been characterized by recent advances in cancer genome sequencing, the germline and environmental effects on GC and their ethnic differences have poorly been understood. Here, we performed genomic scale trans-ethnic analysis of 531 GC cases (319 Asian and 212 non-Asians), by integrating east-Asian and public GC datasets. There is one distinct GC subclass with clear alcohol-associated mutation signature and strong Asian specificity, and almost all the cases in the subclass are attributable to alcohol intake behavior, smoking habit, and Asian-specific ALDH2 defective allele. Alcohol-related GCs have low mutation burden and a characteristic immune profile of increased B-cell infiltration into the tumor area in relation to cancer cell-intrinsic CXCL13 chemokine expression. In addition to germline variants of genes involved in BRCA pathways, we found frequent (7.4%) germline CDH1 variants among Japanese GCs. Most of them were attributed to a few recurrent SNVs shared by both Japanese and Koreans, suggesting the existence of common ancestral events and widespread distributions of these pathogenic variants among east-Asian populations. Specifically, approximately one-fifth of diffuse-type GCs were attributable to the combination of alcohol intake and defective-ALDH2 allele or to germline CDH1 variants. These results revealed uncharacterized impacts of germline variants and their interplays with lifestyles in the high incidence areas.

  Study JGAS000229

• Sequence analysis of colorectal serrated lesions

  Recent studies have shown that colorectal serrated lesions, which include sessile serrated adenomas (SSAs) and traditional serrated adenomas (TSAs), are precursors of colorectal cancer. However, the molecular mechanisms underlying the carcinogenesis, particularly in TSAs, remain largely uncharacterized. To clarify their molecular and clinicopathological characteristics, we performed mutation and methylation analyses of cancer-associated genes in 78 serrated lesions, including TSAs, SSAs and microvesicular hyperplastic polyps. Target exon sequence analysis was performed with 39 genes, including genes known to be frequently mutated in colorectal cancers and/or serrated lesions. We also used bisulfite pyrosequencing to assess the methylation status of various cancer-associated genes and marker genes of the CpG island methylator phenotype (CIMP). The prevalence of mutations in genes associated with Wnt signaling was significantly higher in TSAs than SSAs (65% vs. 28%, p < 0.01). Among those, RNF43 mutations were observed in 38% of TSAs and 17% of SSAs. In immunohistochemical studies of 39 serrated lesions, the prevalence of abnormal nuclear ?-catenin accumulation was significantly higher in TSAs (57%) than SSAs (8%) (P = 0.01). SMOC1 methylation was detected in 54% of TSAs but in no SSAs (p < 0.01). Additionally, SMOC1 methylation was more prevalent among TSAs with KRAS mutation (82%) than with BRAF mutation (38%, p = 0.03). Lesions with CIMP-high or RNF43 mutations were detected only in TSAs with BRAF mutation, suggesting two distinct carcinogenic pathways in TSAs. Mutations in genes associated with Wnt signaling play a greater role in the carcinogenesis of TSAs than SSAs.

  Study JGAS000217

• Defined Lifestyle and Germline Factors Predispose Asian Populations to Gastric Cancer

  Gastric cancer (GC) is one of the leading causes of cancer mortality throughout the world, with the highest incidence in east-Asia. Although the somatic genetics of GC has extensively been characterized by recent advances in cancer genome sequencing, the germline and environmental effects on GC and their ethnic differences have poorly been understood. Here, we performed genomic scale trans-ethnic analysis of 531 GC cases (319 Asian and 212 non-Asians), by integrating east-Asian and public GC datasets. There is one distinct GC subclass with clear alcohol-associated mutation signature and strong Asian specificity, and almost all the cases in the subclass are attributable to alcohol intake behavior, smoking habit, and Asian-specific ALDH2 defective allele. Alcohol-related GCs have low mutation burden and a characteristic immune profile of increased B-cell infiltration into the tumor area in relation to cancer cell-intrinsic CXCL13 chemokine expression. In addition to germline variants of genes involved in BRCA pathways, we found frequent (7.4%) germline CDH1 variants among Japanese GCs. Most of them were attributed to a few recurrent SNVs shared by both Japanese and Koreans, suggesting the existence of common ancestral events and widespread distributions of these pathogenic variants among east-Asian populations. Specifically, approximately one-fifth of diffuse-type GCs were attributable to the combination of alcohol intake and defective-ALDH2 allele or to germline CDH1 variants. These results revealed uncharacterized impacts of germline variants and their interplays with lifestyles in the high incidence areas.

  Study JGAS000228

• Analysis of mutational and proteomic heterogeneity of gastric cancer to monitor post-treatment tumor burden using circulating tumor DNA

  It remains unclear whether complete spectrum of tumor heterogeneity is reflected in circulating tumor DNA (ctDNA). Here, we performed multi-regional target sequencing of primary tumors from 10 patients who underwent surgical resection for gastric cancer. A total of 126 non-synonymous mutations were identified from 30 samples from 10 tumors. Of these, 16 (12.7%) were likely mutations that occurred early in tumor development (i.e., present in all three regions) and were designated founder mutations. Among those validated with highly sensitive digital PCR (dPCR), 93% of non-founder mutations were found in all three regions. Variant allele frequencies (VAFs) of founder mutations were significantly higher than those of non-founder mutations. Phylogenetic analysis demonstrated that five of eight tumors (62.5%, excluding two hypermutator) had at least one identifiable truncal/founder mutation. These findings led us to prioritize founder mutations for ctDNA detection with individually-designed primer/probe sets for dPCR. Marked proteomic diversity was present across the tumor regions both within and between patients independent of mutational status. In preoperative plasma, the average ctDNA VAF of founder mutations was significantly higher than that of non-founder mutations (p = 0.039). Our present data demonstrated that high VAF mutations in the tumor likely represent founder or truncal mutations. However, mutation patterns alone appear insufficient to predict protein heterogeneity. These results suggest that, in practice, mutations having high VAF identified without multi-regional sequencing may be useful for ctDNA monitoring but are inadequate to fully predict proteomic heterogeneity of the tumor.

  Study JGAS000231

• Innate myeloid cell sbuset-specific gene expression patterns in the human colon are altered in Crohn's disease patients

  Background: There is a heterogeneous subset innate myeloid cells, such as macrophages and dendritic cells, in the human intestinal lamina propria. Several studies have demonstrated that these cells contribute to the maintenance of gut homeostasis through induction of inflammatory responses and tolerance via cell type-specific mechanisms; whereas, disrupted innate immune responses are implicated in the pathogenesis of Crohn?s disease. However, the detailed mechanisms by which each innate myeloid subset regulates gut homeostasis and inflammation largely remains unknown. Aim: We aimed to clarify the comprehensive gene expression profiles of innate myeloid cell subsets in the lamina propria from normal human colons (NC) and the inflamed colon sites from patients with Crohn?s disease (CDi). Methods: We performed RNA-sequencing analysis and precise bioinformatics analysis on three innate myeloid cell subsets, CD14-CD11c-, CD14-CD11c+, and CD14+CD11c+CD163low cells from NC and CDi. Results: Transcriptional analysis of the three subsets from the NC showed distinct gene expression patterns and gene ontology (GO) enrichment analysis revealed the associated innate myeloid subset-specific biological process (BP) terms. In addition, changes in gene expression patterns were observed in innate myeloid subsets from CDi. Furthermore, the core GO-BP terms for the genes upregulated in the innate myeloid cells from CDi were distinct from those found in NC. Conclusion: Our data identified the innate myeloid cell subset-specific transcriptomes and the associated enriched GO-BP terms in the NC and found these patterns were altered in CDi.

  Study JGAS000127

• DNA demethylation is associated with malignant progression of low-grade gliomas

  To elucidate the mechanisms of tumor progression is key to discovering strategies for treating patients with lower-grade glioma. Thus to characterize the changes in the molecular profile during malignant progression, we performed methylation array analysis (122 tumors), whole-exome sequencing (36 tumors and matched normal samples), and RNA sequencing (31 tumors). This cohort included 24 matched pairs of initial lower-grade gliomas and recurrent tumors, most of which showed malignant progression when they recurred. Nearly half of IDH-mutant glioblastomas that had progressed from lower-grade gliomas exhibited characteristic partial DNA demethylation in previously methylated genomic regions of their corresponding initial tumors, which had the glioma CpG island methylator phenotype (G-CIMP). In these glioblastomas, cell cycle-related genes were upregulated, RB and PI3K-AKT pathway genes were frequently altered. Notably, bioinformatics analysis revealed that late-replicating domain was significantly enriched in the demethylated regions that were mostly located in non-regulatory regions, suggesting that the loss of DNA methylation during malignant transformation may involve mainly passive demethylation due to a delay in maintenance of methylation during accelerated cell division. Nonetheless, a limited number of genes including IGF2BP3, which potentially drives cell proliferation, were presumed to be upregulated due to demethylation of their promoter. Thus, cell proliferation was further promoted by passive demethylation in the promoter region of these oncogenic genes. Our data indicated that demethylation of the G-CIMP profile found in a subset of recurrent gliomas reflects accelerated cell divisions accompanied by malignant transformation. Oncogenic genes activated by such epigenetic change represent potential therapeutic targets.

  Study JGAS000146

• Gene expression of human Th17 cells before and after activation

  Although gene expression of human Th17 clones has been studies, to date, molecular mechanisms for human Th17 development remain unclear. We have analyzed gene expression of human Th17 cells before and after activation. Th17-cell clones from peripheral blood CD4 T cells cultured with LPS-stimulated dendritic cells derived from monocytes stimulated by TNF were established by a limited dilution method. One of the Th17 clones was proliferated significantly, and high percentages (>70%) of Th17 cells in the cultured cells were maintained by addition of monocytes, IL-1b, and IL-15. CD161+CXCR3- cells (CD161+Th17) in the cultures were sorted and purified up to approximately 97% of IL-17-producing cells after treatment with PMA/ionomycin. Furthermore, these cells were treated with TCR-stimuli, IL-1b, and IL-23 for 2 days, and IL-17-producing Th17 cells (active-Th17) detected by IL-17-catch-up antibodies were sorted. Total RNA sampled from 1x 106 CD161+Th17, active-Th17, na?ve CD4 T cells, or CXCR3+ Th1 cells from combined CD4 T and dendritic cell cultures was used for SAGE. SAGE was performed by using SOLID TM SAGETM kit with barcoding adaptor molecules (Applied Biosystems). DNA sequence was analyzed by using AB5500 Genetic Analyzer (Applied Biosystems), and gene expression was analyzed by using SOLID? SAGETM software (Applied Biosystems). Our data shows high levels of IL-17A-mRNA expression by active-Th17; however, but only slight expression by CD161+Th17. Therefore, the data could provide alteration of gene expression by Th17 cells before and after stimulation.

  Study JGAS000005

• Identification of driver oncogenes in scirrhous-type gastric cancer cell lines

  Scirrhous type gastric cancer (SGC) is one of the most intractable cancer subtypes in human, and its therapeutic targets have been rarely identified to date. Exploration of somatic mutations in the SGC genome with the next generation sequencers has been hampered by low tumor contents by the markedly increased fibrous tissues. Thus, SGC cell lines may be useful resources to search for novel oncogenes in. Here we have conducted whole exome sequencing and RNA sequencing on two SGC cell lines, OCUM-8 and OCUM-9. Interestingly, most of the mutations thus identified have not been reported. In OCUM-8 cells, a novel CD44-IGF1R fusion gene is discovered, the protein product of which ligates the amino-terminus of CD44 to the transmembrane and the tyrosine-kinase domains of IGF1R. Further, both of CD44 and IGF1R are markedly amplified in the OCUM-8 genome and abundantly expressed. CD44-IGF1R has a transforming ability, and the suppression of its kinase activity leads to rapid cell death of OCUM-8. On the other hand, OCUM-9 seems to possess multiple oncogenic events in its genome. The MET gene is highly amplified, and the treatment of OCUM-9 with crizotinib, a MET inhibitor, results in cell death. Further, a novel BORCS5-ETV6 fusion gene is identified in the OCUM-9 genome. BORCS5-ETV6 possesses oncogenic activity, and suppression of its message partially inhibits cell growth. Prevalence of these novel fusion genes among SGC awaits further investigation, but we validate the significance of cell lines as appropriate reagents for detailed genomic analyses of SGC.

  Study JGAS000179

• Heterogeneity of IKZF1 genomic alterations and risk of relapse in childhood B-cell precursor acute lymphoblastic leukemia

  Purpose Genomic alterations of IKZF1 are common and associated with adverse clinical features in B-ALL, but the relationship between the type of IKZF1 alteration, disease subtype and outcome are incompletely understood. Patients and Methods We used transcriptome and genomic sequencing, and SNP microarray data to determine leukemia subtype and genomic alterations for 688 pediatric patients with B-ALL enrolled into St. Jude Total Therapy 15 and 16 studies. Results IKZF1 alterations were identified in 115 (16.7%) patients, commonly in BCR::ABL1 (78%) and CRLF2-rearranged, BCR::ABL1-like B-ALL (70%), and were associated with five-year cumulative incidence of relapse (CIR) of 14.8 ± 3.3% compared to 5.0 ± 0.9% for patients without any IKZF1 alteration (P < 0.0001). IKZF1 deletions of exon 4-7 (P = 0.0002), genomic IKZF1plus with any IKZF1 deletion (P = 0.006) or with focal IKZF1 deletion (P = 0.0007), and unfavorable genomic subtypes (P < 0.005) were independently adversely prognostic. Patients with both IKZF1 exon 4-7 deletion and unfavorable genomic subtype had a significantly increased risk of relapse (HR = 58.3; 95% CI, 11.9 – 285.4, P < 0.0001) whereas IKZF1 sequence mutations were not independently predictive of outcome. Conclusions Genomic IKZF1plus with any IKZF1 deletion, IKZF1 deletion of exon 4-7, and unfavorable subtype confer increased risk of relapse. The combination of unfavorable genomic subtype and IKZF1 deletion of exon 4-7 identified patients at greatest risk of relapse despite MRD-directed therapy. The type of IKZF1 alteration together with subtype are informative for risk stratification and predict response in patients with B-ALL.

  Study EGAS50000000106

• Case Report: Pre-Clinical Combination Targeting VEGF and PI3K in a Rare, Aggressive Mixed Endometrial Carcinoma

  We report a rare case of a young patient (VENUS 167) initially diagnosed with grade 1 endometrioid endometrial cancer, which, following endocrine treatment, presented with mixed aggressive carcinoma with three distinct histologic patterns: grade 1 endometrioid, large cell neuroendocrine, and undifferentiated carcinoma. The surgical specimen at the time of disease progression was used to establish OPTO.85, a patient-derived organoid (PDO), followed by a corresponding organoid-derived xenograft (ODX). Multi-omic analyses confirmed that OPTO.85 accurately reflected the patient’s tumor characteristics. WES analysis identified oncogenic alterations in PIK3CA, ARID1A, and CTNNB1. Further RNA-seq and ATAC-seq analyses revealed enrichment in VEGF and Wnt signaling pathways, suggesting therapeutic vulnerabilities. A high-throughput drug screen was conducted using ApexBio-Approved and Epigenetic Drug Libraries, along with Kinase Inhibitor and Tool Compound Libraries developed at the Ontario Institute of Cancer Research. The OPTO.85 PDO exhibited sensitivity to PI3K inhibitors and responsiveness to VEGF inhibition. Cediranib demonstrated synergy with BKM120, significantly reducing organoid growth. This combination also showed in vivo efficacy in the ODX model, where dual inhibitors significantly suppressed tumor growth compared to single compounds. This case exemplifies the impact of genomic profiling and patient-derived models in identifying actionable molecular changes in rare cancers with limited therapeutic options and poor prognosis. It highlights that high-throughput sequencing for individual patient tumors and generation of patient-derived models are feasible in endometrial cancer. This preclinical model may assist clinical decision and personalized therapy requiring validation in prospective studies.

  Study EGAS50000001665

• An autoinflammatory RIG-I variant causing Singleton-Merten Syndrome associates with small non-coding Y-RNAs

  The RNA sensor RIG-I performs a critical role in surveying the cytoplasm for the presence of viral nucleic acids and initiating the downstream anti-viral type I interferon pathway. Through recognition of specific features, such as the presence of a 5’ di/triphosphate motif or highly structured base-paired regions, RIG-I effectively differentiates between RNA of viral- and self-origin. In Singleton-Merten syndrome (SMS), gain-of-function variants in RIG-I lead to a breakdown in this surveillance system and results in aberrant sensing of self-RNAs and deleterious upregulation of type I interferons. The identity of the self-RNAs binding to gain-of-function RIG-I mutants has remained elusive and would provide a greater understanding of the aetiology of SMS. Here we used an infrared individual-nucleotide resolution UV-crosslinking and immunoprecipitation (irCLIP) approach to determine the RNA profile bound to a previously characterised ATPase-deficient SMS variant, RIG-IC268F. irCLIP identified a broad array of self-RNAs, primarily those transcribed by RNA polymerase III, that were bound to RIG-I. Subsequent native RNA immunoprecipitation confirmed a prominent and specific interaction between RIG-IC268F and Y-RNAs, a family of four structurally similar non-coding RNAs. Manipulation of Y-RNAs alone by targeting either Y-RNA transcripts or Y-RNA stabilising proteins was insufficient to negate RIG-I induce interferon, hinting at a broader profile of RNA polymerase III-derived RNAs being influential in driving the sterile activation of gain-of-function RIG-I variants.

  Study EGAS50000001660

• Detection of Clinically Relevant Genetic Variantsin Autism Spectrum Disorderby Whole-Genome Sequencing

  Autism Spectrum Disorder (ASD) demonstrates high heritability and familial clustering, yet the genetic causes remain only partially understood as a result of extensive clinical and genomic heterogeneity. Whole-genome sequencing (WGS) shows promise as a tool for identifying ASD risk genes as well as unreported mutations in known loci, but an assessment of its full utility in an ASD group has not been performed. We used WGS to examine 32 families with ASD to detect de novo or rare inherited genetic variants predicted to be deleterious (loss-of-function and damaging missense mutations). Among ASD probands, we identified deleterious de novo muta- tions in six of 32 (19%) families and X-linked or autosomal inherited alterations in ten of 32 (31%) families (some had combinations of mutations). The proportion of families identified with such putative mutations was larger than has been previously reported; this yield was in part due to the comprehensive and uniform coverage afforded by WGS. Deleterious variants were found in four unrecognized, nine known, and eight candidate ASD risk genes. Examples include CAPRIN1 and AFF2 (both linked to FMR1, which is involved in fragile X syndrome), VIP (involved in social-cognitive deficits), and other genes such as SCN2A and KCNQ2 (linked to epilepsy), NRXN1, and CHD7, which causes ASD-associated CHARGE syndrome. Taken together, these results suggest that WGS and thorough bioinformatic analyses for de novo and rare inherited mutations will improve the detection of genetic variants likely to be associated with ASD or its accompanying clinical symptoms.

  Study EGAS00001000850

• Sequencing data from triple-negative breast cancer tumors

  Autoimmunity and anti-cancer immunity lie on the same biological continuum1,2, but their link remains obscure. The paraneoplastic neurological syndrome anti-NMDA receptor (NMDAR) encephalitis (ANRE) is a paradigm for their connectivity3 given that intratumoral NMDAR expression correlates with the generation of anti-NMDAR antibodies4,5. Here, we verify ectopic expression of GluN1 and GluN2B NMDAR sub-units in triple-negative breast cancer (TNBC)6 and model this using orthotopic TNBC tumors with inducible expression of GluN1-GluN2B NMDARs. We show that NMDAR expression is sufficient to induce B cell recruitment and their affinity maturation, consistent with an integrated adaptive immune response. Reconstruction of extended intratumoral B cell phylogenies and cryo-EM structural analyses demonstrated that affinity-matured hypermutated and class-switched antibodies emerged from pre-existing germline-configuration lower-affinity anti-NMDAR antibodies. Distinct matured antibodies targeted specific epitopes and induced conformational rearrangements within the NMDAR amino-terminal domain, predictive of their functional effects, ranging from inhibition to potentiation. Passive transfer of an NMDAR-potentiating antibody caused autonomic dysregulation and lowered the seizure threshold in healthy female mice, recapitulating key diagnostic criteria of ANRE4,5. We further identify a correlation between intratumoral NMDAR expression and anti-NMDAR antibody titers in TNBC patients. Taken together, our data establish a direct connection between intratumoral NMDAR expression, antibody maturation, and the onset of autoimmunity. These findings suggest that germline-encoded anti-NMDAR antibodies contribute to immune surveillance but can also trigger autoimmune disease upon maturation, revealing a mechanistic tradeoff between cancer immunity and neurotoxicity.

  Study EGAS50000001585

• Sequencing data from triple-negative breast cancer tumors

  Autoimmunity and anti-cancer immunity lie on the same biological continuum1,2, but their link remains obscure. The paraneoplastic neurological syndrome anti-NMDA receptor (NMDAR) encephalitis (ANRE) is a paradigm for their connectivity3 given that intratumoral NMDAR expression correlates with the generation of anti-NMDAR antibodies4,5. Here, we verify ectopic expression of GluN1 and GluN2B NMDAR sub-units in triple-negative breast cancer (TNBC)6 and model this using orthotopic TNBC tumors with inducible expression of GluN1-GluN2B NMDARs. We show that NMDAR expression is sufficient to induce B cell recruitment and their affinity maturation, consistent with an integrated adaptive immune response. Reconstruction of extended intratumoral B cell phylogenies and cryo-EM structural analyses demonstrated that affinity-matured hypermutated and class-switched antibodies emerged from pre-existing germline-configuration lower-affinity anti-NMDAR antibodies. Distinct matured antibodies targeted specific epitopes and induced conformational rearrangements within the NMDAR amino-terminal domain, predictive of their functional effects, ranging from inhibition to potentiation. Passive transfer of an NMDAR-potentiating antibody caused autonomic dysregulation and lowered the seizure threshold in healthy female mice, recapitulating key diagnostic criteria of ANRE4,5. We further identify a correlation between intratumoral NMDAR expression and anti-NMDAR antibody titers in TNBC patients. Taken together, our data establish a direct connection between intratumoral NMDAR expression, antibody maturation, and the onset of autoimmunity. These findings suggest that germline-encoded anti-NMDAR antibodies contribute to immune surveillance but can also trigger autoimmune disease upon maturation, revealing a mechanistic tradeoff between cancer immunity and neurotoxicity.

  Study EGAS50000001598

• Single cell RNA-seq mapping of nasal and tracheobronchial airways in human healthy volunteers

  The respiratory tract constitutes an elaborated line of defense based on a unique cellular ecosystem. Single-cell profiling methods enable the investigation of cell population distributions and transcriptional changes along the airways. We have explored cellular heterogeneity of the human airway epithelium in 10 healthy living volunteers by single-cell RNA profiling. 77,969 cells were collected by bronchoscopy at 35 distinct locations, from the nose to the 12th division of the airway tree. The resulting atlas is composed of a high percentage of epithelial cells (89.1%), but also immune (6.2%) and stromal (4.7%) cells with peculiar cellular proportions in different sites of the airways. It reveals differential gene expression between identical cell types (suprabasal, secretory, and multiciliated cells) from the nose (MUC4, PI3, SIX3) and tracheobronchial (SCGB1A1, TFF3) airways. By contrast, cell-type specific gene expression was stable across all tracheobronchial samples. Our atlas improves the description of ionocytes, pulmonary neuro-endocrine (PNEC) and brush cells, which are likely derived from a common population of precursor cells. We also report a population of KRT13 positive cells with a high percentage of dividing cells which are reminiscent of “hillock” cells previously described in mouse. Robust characterization of this unprecedented large single-cell cohort establishes an important resource for future investigations. The precise description of the continuum existing from nasal epithelium to successive divisions of lung airways and the stable gene expression profile of these regions better defines conditions under which relevant tracheobronchial proxies of human respiratory diseases can be developed.

  Study EGAS00001004082

• The genomic landscape of pediatric acute lymphoblastic leukemia

  Acute lymphoblastic leukemia (ALL) is the most common pediatric cancer and arises from B- or T-lineage lymphocyte precursors. ALL comprises dozens of subtypes, and genomic analyses have largely been performed within these subtypes. Here we analyze pediatric ALL genomes across all subtypes, including 768 whole genomes, 1,729 exomes, and 1,889 transcriptomes in 2,754 patients. Most ALL subtypes harbor 4 or more driver alterations per sample, similar to adult cancers, despite low mutation burdens. Hyperdiploid B-ALL copy gains are likely acquired early and synchronously, with copy gains occurring before ultraviolet-induced mutations. By contrast, ultraviolet-induced mutations precede copy gains in iAMP21 B-ALL. Overall, we identified 378 putative ALL driver genes. Most driver alterations vary in prevalence across ALL subtypes, with B-ALL enriched for Ras and B-lineage-related alterations, and T-ALL enriched for PI3K, JAK, and cell cycle alterations. Most B-ALL (54.3%) and T-ALL (51.2%) samples bear at least one rare driver gene alteration (present in less than 2% of samples), including 30 putative novel cancer driver genes associated with ubiquitination, SUMOylation, non-coding, and other functions. Known or novel alterations associated with poor outcomes in specific subtypes include TBL1XR1 alterations in ETV6-RUNX1 patients, CREBBP in Ph-like-CRLF2, SETD2 in hyperdiploid, and PTEN in TAL1 patients. Intriguingly, DUX4 and KMT2A subtypes separate into CEBPA/FLT3- or NFATC4-expressing subgroups with potential clinical implications. Together, these results deepen understanding of ALL etiology and outcomes, and facilitate ALL model development.

  Study EGAS00001005250

• Analysis of recent shared ancestry in a familial cohort identifies coding and noncoding autism spectrum disorder variants

  Autism spectrum disorder (ASD) is a collection of neurodevelopmental disorders characterized by deficits in social communication and restricted, repetitive patterns of behavior or interests. ASD is highly heritable, but genetically and phenotypically heterogeneous, reducing the power to identify causative genes. We performed whole genome sequencing (WGS) in an ASD cohort of 68 individuals from 22 families enriched for recent shared ancestry. We identified an average of 3.07 million variants per genome, of which an average of 112,512 were rare. We mapped runs of homozygosity (ROHs) in affected individuals and found an average genomic homozygosity of 9.65%, consistent with expectations for multiple generations of consanguineous unions. We identified potentially pathogenic rare exonic or splice site variants in 12 known (including KMT2C, SCN1A, SPTBN1, SYNE1, ZNF292) and 12 candidate (including CHD5, GRB10, PPP1R13B) ASD genes. Furthermore, we annotated noncoding variants in ROHs with brain-specific regulatory elements and identified putative disease-causing variants within brain-specific promoters and enhancers for 5 known ASD and neurodevelopmental disease genes (ACTG1, AUTS2, CTNND2, CNTNAP4, SPTBN4). We also identified copy number variants in two known ASD and neurodevelopmental disease loci in two affected individuals. In total we identified potentially etiological variants in known ASD or neurodevelopmental disease genes for ~61% (14/23) of affected individuals. We combined WGS with homozygosity mapping and regulatory element annotations to identify candidate ASD variants. Our analyses add to the growing number of ASD genes and variants and emphasize the importance of leveraging recent shared ancestry to map disease variants in complex neurodevelopmental disorders.

  Study EGAS00001006058

• RB Loss in Resistant EGFR Mutant Lung Adenocarcinomas that Transform to Small Cell Lung Cancer

  Tyrosine kinase inhibitors (TKIs) are effective treatments for non-small cell lung cancers (NSCLCs) with epidermal growth factor receptor (EGFR) mutations. However, they do not lead to cures, and, on average, relapse occurs typically occurs afterafter an average of one year of continuous treatment. In a subset of patients, aA fundamental histological transformation from NSCLC to small cell lung cancer (SCLC) is observed in a subset of the resistant cancers, but the molecular changes associated with this transformation remain unknown. Analysis of a cohort of tumor samples and cell lines derived from resistant EGFR mutant patients with SCLC transformation revealed that RB is lost in 100% of these SCLC transformed cases, but rarely in those that remain NSCLC. Global changes in gene expression, including increased neuroendocrine marker expression and absence of EGFR expression, are observed in resistant EGFR mutant cancers that transformed to SCLC. Consistent Further, with their similarities to classical SCLC at the genetic and gene expression levelincreased neuroendocrine marker and decreased EGFR expression as well as greater sensitivity to BCL2 family inhibition are observed in resistant SCLC transformed cancers compared to resistant NSCLCs. , cell lines derived from resistant EGFR mutant SCLC biopsies are substantially more sensitive to Bcl-2 family inhibition with ABT-263 compared to those derived from resistant EGFR mutant NSCLCs. Together, these findings suggest that despite developing initially as EGFR mutant adenocarcinomas, this subset of resistant cancers ultimately adopt many of the molecular and phenotypic characteristics of classical SCLC.

  Study EGAS00001001102

• Pediatric Papillary Thyroid Carcinoma Whole Exome Sequencing

  BACKGROUND: Pediatric papillary thyroid carcinoma (PPTC) is distinct from adult-onset disease. Despite higher rates of regional and distal metastases and more frequent recurrence, there is, paradoxically, a highly favourable prognosis. Historically, only 40-50% of childhood PTC are driven by genomic variants common to adult PTC, while the oncogenic drivers in the remainder are unknown. METHODS: Fifty-two PPTCs were genotyped by candidate gene testing, followed, in a subset (n=32), by whole exome (WES) and transcriptome sequencing (RNASeq) to identify oncogenic variants in driver-unknown tumors and to explore gene expression patterns. FINDINGS: Thirty-one tumors (60%) carried variants identified by candidate gene testing, specifically: BRAFV600E(22%), NRASQ61R(2%), RET-CCDC6(19%), RET-NCOA4(15%) and PAX8-PPARG. Among the 21 remaining tumors with negative first-line genotyping, oncogenic variants were identified in 19 (90%). These were enriched with oncogenic fusions, with 11 non-recurrent fusion transcripts identified, including two that had not been previously described: STRN-RET and TG-PBF. Most fusions were associated with 3’ receptor tyrosine kinase (RTK) moieties: RET, MET, ALK and NTRK3. Gene expression analysis defined three distinct clusters, which demonstrated differing expression of genes involved in thyroid differentiation and MAPK signaling. INTERPRETATION: Exome and transcriptome-wide sequencing identify genomic drivers of PPTC in 98% of tumors, whereas conventional genotyping leaves many tumors genetically unexplained. PPTC is driven, in the majority of cases, by oncogenic fusion transcripts involving one of a limited number of RTKs. These observations have broad implications for diagnostic testing and for therapeutic options in children with advanced disease.

  Study EGAS00001005187

• WES analysis of a mixed cohort of pituitary tumors

  BACKGROUND: Sporadic pituitary adenomas are the second most common primary brain tumor, but their genetics remains incompletely understood. We performed deep DNA sequencing of 42 pituitary macroadenomas, including hormonally active and atypical adenomas. METHODS: Pituitary adenomas representative of a spectrum of histopathologic subtypes, primary and recurrent tumors, and MIB-1 proliferative indices were profiled. Whole-exome sequencing was performed on tumors and paired blood followed by analysis of somatic copy number alteration and mutation. RESULTS: On the whole, the tumors harbored a low rate of mutations, similar to other benign brain tumors. However, almost one-third (29%) exhibited somatic copy-number alterations across large fractions of the genome (up to 99%). Despite such widespread genomic disruption, these tumors had few focal events, which is atypical among highly-disrupted cancers. The other 71% of tumors formed a distinct molecular class, with somatic copy-number alterations involving less than 6% of the genome. Among the highly-disrupted group, 92% were functional or atypical null-cell adenomas, whereas 87% of the less-disrupted group were clinically non-functional adenomas. The only significantly recurrent events were arm-level losses of chromosomes 1, 2, 11, and 18. No significantly recurrent mutations or insertion/deletion events were identified, suggesting no genes are altered by exomic mutations in &gt;15% of pituitary tumors. CONCLUSIONS: Genomic characterization of a large series of sporadic pituitary adenomas reveals two classes of tumor, based on genomic disruption from arm-level copy-number alterations, which associate with hormonal and histologic subtypes.

  Study EGAS00001001714

• C3 SNPs and outcome after lung transplantation

  Background: Development of chronic rejection is still a severe problem and causes high mortality rates after lung transplantation (LTx). Complement activation is important in the development of acute rejection (AR) and bronchiolitis obliterans syndrome, with C3 as a key complement factor. Methods: Since a single nucleotide polymorphism (SNP) in the C3 gene (rs2230199) is associated with long-term outcome after solid organ transplantation, we investigated this SNP in relation to LTx. In addition, we looked at local production of C3 by analyzing bronchoalveolar lavage fluid (BALF) of LTx patients using isoelectric focusing (IEF). Results: We demonstrated the presence of C3 in BALF and showed that this is produced by the donor lung based on the genotype of SNP rs2230199. We also analyzed donor and patient SNP configurations and observed a significant association between the SNP configuration in patients and episodes of AR during 4-years follow-up. Survival analysis showed a lower AR-free survival in homozygous C3 slow patients (p=0.005). Furthermore, we found a significant association between the SNP configuration in donors and BOS development. Patients receiving a graft from a donor with at least one C3 fast variant for rs2230199 had an inferior BOS-free survival (p=0.044). Conclusions: In conclusion, our data indicate local C3 production by donor lung cells. In addition, a single C3 SNP present in recipients affects short-term outcome after LTx, while this SNP in donors has an opposite effect on long-term outcome after LTx. These results could contribute to an improved risk stratification after transplantation.

  Study EGAS00001003843

• whole genome sequence data of multifocal hepatocellular carcinoma

  Background &amp; Aims: The differentiation of distinct multifocal hepatocellular carcinoma (HCC): multicentric disease vs. intrahepatic metastases, in which the management and prognosis varies substantively, remains problematic. We aim to stratify multifocal HCC and identify novel diagnostic and prognostic biomarkers by performing whole genome and transcriptome sequencing, as part of a multi-omics strategy Methods: A complete collection of tumour and somatic specimens (intrahepatic HCC lesions, matched non-cancerous liver tissue and blood) were obtained from representative patientswith multifocal HCC exhibiting two distinct postsurgical courses. Whole-genome and transcriptome sequencing with genotyping were performed for each tissue specimen to con-trast genomic alterations, including hepatitis B virus integrations, somatic mutations, copy number variations, and structural variations. We then constructed a phylogenetic tree to visualise individual tumour evolution and performed functional enrichment analyses on select differentially expressed genes to elucidate biological processes involved in multifocal HCC development. Multi-omics data were integrated with detailed clinicopathological information to identify HCC biomarkers, which were further validated using a large cohort of HCC patients (n = 174).Results: The multi-omics pro?ling and tumour biomarkers could successfully distinguish the two multifocal HCC types, while accurately predicting clonality and aggressiveness. The dual speci?city protein kinase TTK, which is a key mitotic checkpoint regulator with links to p53 signaling, was further shown to be a promising overall prognostic marker for HCC in the large patient cohort.Conclusions: Comprehensive multi-omics characterisation of multifocal tumour evolution may improve clinical decision-making, facilitate personalised medicine, and expedite identi?cation of novel biomarkers and therapeutic targets in HCC.

  Study EGAS00001002338

• Microsatellite data show recent demographic expansions in sedentary but not in nomadic human populations in Africa and Eurasia.

  The transition from hunting and gathering to plant and animal domestication was one of the most important cultural and technological revolutions in human history. According to archeologists and paleoanthropologists, this transition triggered major demographic expansions. However, few genetic studies have found traces of Neolithic expansions in the current repartition of genetic polymorphism, pointing rather toward Paleolithic expansions. Here, we used microsatellite autosomal data to investigate the past demographic history of 87 African and Eurasian human populations with contrasted lifestyles (nomadic hunter-gatherers, semi-nomadic herders and sedentary farmers). Likely due to the combination of a higher mutation rate and the possibility to analyze several loci as independent replicates of the coalescent process, the analysis of microsatellite data allowed us to infer more recent expansions than previous genetic studies, potentially resulting from the Neolithic transition. Despite the variability in their location and environment, we found consistent expansions for all sedentary farmers, while we inferred constant population sizes for all hunter-gatherers and most herders, which could result from constraints linked to a nomadic or semi-nomadic lifestyle and/or competition for land between herders and farmers. As an exception, we inferred expansions for Central Asian herders. This might be linked with the arid environment of this area, which may have been more favorable to nomadic herders than to sedentary farmers. Alternatively, current Central Asian herders may descent from populations who have first experienced a transition from hunter-gathering to sedentary agropastoralism, and then a second transition to nomadic herding.

  Study EGAS00001000652

• High-throughput Detection of Clinically Relevant Mutations in Archived Tumor Samples By Multiplexed PCR and Next Generation Sequencing

  PURPOSE: Tailoring cancer treatment to tumor molecular characteristics promises to make personalized medicine a reality. However, reliable genetic profiling of archived clinical specimens has been hindered by limited sensitivity and high false positive rates. Here, we describe a novel methodology, MMP-seq, which enables sensitive and specific high-throughput, high content genetic profiling in archived clinical samples. EXPERIMENTAL DESIGN: We first validated MMP-seq’s technical performance in 66 cancer cell lines and a Latin square cross-dilution of known somatic mutations. We next characterized the performance of MMP-seq in 17 formalin-fixed paraffin embedded (FFPE) clinical samples using matched fresh frozen (FF) tissue from the same tumors as benchmarks. To demonstrate the potential clinical utility of our methodology, we profiled FFPE tumor samples from 73 endometrial cancer patients. RESULTS: We demonstrated that MMP-seq enabled rapid and simultaneous profiling of a panel of 88 cancer genes in 48 samples, and detected variants at frequencies as low as 0.4%. We identified DNA degradation and deamination as the main error sources and developed practical and robust strategies for mitigating these issues, and dramatically reduced the false positive rate. Applying MMP-seq to a cohort of endometrial tumor samples identified extensive, potentially actionable alterations in the PI3K and RAS pathways, including novel PIK3R1 hotspot mutations that may disrupt negative regulation of PIK3CA. CONCLUSIONS: MMP-seq provides a robust solution for comprehensive, reliable and high- throughput genetic profiling of clinical tumor samples, paving the way for the incorporation of genomic-based testing into clinical investigation and practice.

  Study EGAS00001000674

• Exome sequencing of uterine leiomyosarcomas

  Uterine leiomyosarcomas (ULMSs) are aggressive smooth muscle tumors associated with poor clinical outcome. Despite previous cytogenetic and molecular studies, their molecular background has remained elusive. To examine somatic variation in ULMS, we performed exome sequencing on 19 tumors. Altogether 43 genes were mutated in at least two ULMSs. Most frequently mutated genes included tumor protein P53 (TP53; 6/19; 33%), alpha thalassemia/mental retardation syndrome X-linked (ATRX; 5/19; 26%), and mediator complex subunit 12 (MED12; 4/19; 21%). Unlike ATRX mutations, both TP53 and MED12 alterations have repeatedly been associated with ULMSs. All the observed ATRX alterations were either nonsense or frameshift mutations. ATRX protein levels were reliably analyzed by immunohistochemistry in altogether 44 ULMSs, and the majority of tumors (23/44; 52%) showed clearly reduced expression. Loss of ATRX expression has been associated with alternative lengthening of telomeres (ALT), and thus the telomere length was analyzed with telomere-specific fluorescence in situ hybridization. The ALT phenotype was confirmed in all ULMSs showing diminished ATRX expression. Exome data also revealed one nonsense mutation in death-domain associated protein (DAXX), another gene previously associated with ALT, and the tumor showed ALT positivity. Aberrant expression of both TP53 and ATRX were associated with poor overall survival. In conclusion, exome sequencing revealed that TP53, ATRX, and MED12 are frequently mutated in ULMSs. ALT phenotype was commonly seen in tumors, indicating that ATR inhibitors, which were recently suggested as possible new drugs for ATRX-deficient tumors, could provide a potential novel therapeutic option for ULMS.

  Study EGAS00001001612

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

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

  Study EGAS00001004063