14577 results for "PLOS+ONE+Pennsylvania+German+ACS+language+code+pdc+open+access+study"

in 64.51 milliseconds.

• University of Texas PDX Development and Trial Center Grant

  The goal for the University Texas PDX Development and Trial Center (UTPDTC) is to optimize personalized biomarker-based cancer therapy and identify effective targeted drugs based on the molecular characteristics of each tumor. Our short-term goals are to establish a biobank of clinically, and molecularly-annotated Patient-Derived Xenografts (PDXs) and to use PDXs as a platform for preclinical drug development and biomarker discovery. The primary goal for UTPDTC investigators will be to develop PDX trial strategies for preclinical testing of single agents and drug combinations. These models will allow the determination of the optimal treatments (single drugs or combinations) that should be tested in clinical trials in increasingly individualized, molecularly defined subsets of tumors. The goal of the Patient-Derived Xenograft Core is to provide high-quality clinically relevant and molecularly annotated PDX models for the research projects proposed in the University of Texas PDX Development and Trial Center (UTPDTC) grant application and to the research activity of the NCI PDX Development and Trial Centers Research Network (PDXNet) by leveraging PDX resources at our institutions and developing new PDX models from human cancer specimens using rigorous quality standards so that the models can be used to guide clinical trial development. The PDX models developed and/or characterized by the PDX Core will be available to other cancer researchers through the PDXNet and NCI Patient-Derived Models Repository (PDMR).

  Study phs001980

• Small Intestine Neuroendocrine Tumors (Carcinoid Tumors)

  The diagnosed incidence of small bowel neuroendocrine tumors (NETs) is increasing. While patients with localized disease can be treated surgically, those with metastatic disease currently have few treatment options. The success of biologically targeted therapies in other malignancies has led to interest in the molecular alterations underlying the pathogenesis of these NETs. To identify genetic aberrations in small intestine NETs, we generated copy number profiles from 31 primary and metastatic tumors and performed whole-exome sequencing on a subset of 29 primary small intestine NETs and 24 metastatic NETs in parallel with normal blood DNA. Whole-genome sequencing data was generated on 15 tumor/normal pairs and 5 primary/metastasis/normal trios. The global genetic landscape of small bowel NETs is relatively quiet. Consistent with previous studies, the overwhelming majority of tumors were characterized by loss of chromosome 18 and, to a lesser extent, other chromosome arm gains and losses. In stark contrast to arm-level alterations, recurrent high-level focal amplifications and deletions were much less prevalent in these tumors. High-throughput mutation screening and exome sequencing revealed similarly low rates of somatic mutation in NETs (median of 0.77 non-silent mutations per megabase (Mb) of coding DNA) compared to other recent cancer exome sequencing efforts. Our analysis of this cohort identified only a single, statistically significant recurrent somatic mutation targeting the cyclin-dependent kinase inhibitor gene, CDKN1B, encoding p27.

  Study phs000579

• Exome Sequencing of Esophageal Adenocarcinoma

  The incidence of esophageal adenocarcinoma (EAC) has risen 600% over the last 30 years. When coupled to the five-year survival rate of only 15% for this disease, identification of new therapeutic targets for EAC is of great importance. Here, we analyze the mutational spectra identified from the whole genome sequencing of 16 EACs and complete exome sequencing of 149 EAC tumors and matched germline samples. We identify a novel mutation signature marked by high prevalence of A to C transversions at AA*(N) trinucleotides. Notably, the development of EAC is preceded by pathologic intestinal metaplasia of the lower esophagus, Barrett's esophagus, itself a response to injury from gastric-esophageal reflux disease (GERD). Thus, we hypothesize to represent these AA*(N) mutations to a novel signature of GERD-induced mutagenesis. Analysis of the exome data identified 26 genes subject to statistically significant recurrent mutation. Of these 26 genes, only TP53, CDKN2A, SMAD4, and PIK3CA had been previously implicated in EAC. Novel mutations include those targeting chromatin modifying factors and novel candidate contributors to EAC: SPG20, TLR4, ELMO1 and DOCK2. Of these, ELMO1 and DOCK2 are notably direct dimerization partners that act to activate Rac1 to enhance cellular invasiveness, thus generating new hypotheses about the potential for the Rac1 pathway to be a novel target for therapy of EAC. "Reprinted from 'Unraveling the mutational complexity of esophageal adenocarcinoma', with permission from Nature Genetics."

  Study phs000598

• Region-Specific Neural Stem Cell Lineages Revealed by Single-Cell RNA-Sequences from Human Embryonic Stem Cells

  During development of the human brain, multiple cell types with diverse regional identities are generated. Here we report a system to generate early human brain forebrain and mid/hindbrain cell types from human embryonic stem cells (hESCs), and infer and experimentally confirm a lineage tree for the generation of these types based on single-cell RNA-Seq analysis. We engineered SOX2Cit/+ and DCXCit/Y hESC lines to target progenitors and neurons throughout neural differentiation for single-cell transcriptomic profiling, then identified discrete cell types consisting of both rostral (cortical) and caudal (mid/hindbrain) identities. Direct comparison of the cell types were made to primary tissues using gene expression atlases and fetal human brain single-cell gene expression data, and this established that the cell types resembled early human brain cell types, including preplate cells. From the single-cell transcriptomic data a Bayesian algorithm generated a unified lineage tree, and predicted novel regulatory transcription factors. The lineage tree highlighted a prominent bifurcation between cortical and mid/hindbrain cell types, confirmed by clonal analysis experiments. We demonstrated that cell types from either branch could preferentially be generated by manipulation of the canonical Wnt/beta-catenin pathway. In summary, we present an experimentally validated lineage tree that encompasses multiple brain regions, and our work sheds light on the molecular regulation of region-specific neural lineages during human brain development.

  Study phs001205

• Use of Whole Genome Sequencing for Diagnosis and Discovery in the Cancer Genetics Clinic

  Despite the potential of whole-genome sequencing (WGS) to improve patient diagnosis and care, the empirical value of WGS in the cancer genetics clinic is unknown. We performed WGS on members of two cohorts of cancer genetics patients: those with BRCA1/2 mutations (n = 176) and those without (n = 82). Initial analysis of potentially pathogenic variants (PPVs, defined as nonsynonymous variants with allele frequency < 1% in ESP6500) in 163 clinically-relevant genes suggested that WGS will provide useful clinical results. This is despite the fact that a majority of PPVs were novel missense variants likely to be classified as variants of unknown significance (VUS). Furthermore, previously reported pathogenic missense variants did not always associate with their predicted diseases in our patients. This suggests that the clinical use of WGS will require large-scale efforts to consolidate WGS and patient data to improve accuracy of interpretation of rare variants. While loss-of-function (LoF) variants represented only a small fraction of PPVs, WGS identified additional cancer risk LoF PPVs in patients with known BRCA1/2 mutations and led to cancer risk diagnoses in 21% of non-BRCA cancer genetics patients after expanding our analysis to 3209 ClinVar genes. These data illustrate how WGS can be used to improve our ability to discover patients' cancer genetic risks. "Reprinted from doi:10.1016/j.ebiom.2014.12.003, with permission from EBioMedicine."

  Study phs000942

• Targeted Linked-Read DNA-seq Analysis of Castration-Resistant Prostate Cancers

  Androgen receptor (AR) inhibition is standard of care for advanced prostate cancer (PC). However, its efficacy is limited by progression to castration-resistant PC (CRPC), usually due to AR re-activation via mechanisms that include AR amplification and structural rearrangement. These two classes of AR alterations often co-occur in CRPC tumors, but it is unclear whether this reflects intercellular or intracellular heterogeneity of AR. It is important to resolve the mechanisms of AR genomic rearrangements to develop new therapies and predictive biomarkers for PC. We developed a targeted linked-read DNA sequencing assay to interrogate alterations in AR and a panel of genes known to be amplified or mutated in CRPC. Linked-read DNA-sequencing libraries were generated using the 10X Genomics Chromium platform and subjected to hybridization enrichment using a custom-designed Agilent SureSelect bait panel. Post-capture DNA-seq libraries were pooled and sequenced on an Illumina HiSeq 2500 using 2x150 bp settings. This targeted linked-read DNA-seq assay was used to interrogate a set of 6 benchmarking samples harboring known AR gene rearrangements and other alterations in AR, as well a panel of 6 patient-derived xenograft CRPC tumors and 23 clinical CRPC tumors. In this cohort, we identified a subset of samples that developed complex, multiply-rearranged AR gene structures in conjunction with very high AR copy number. The data that will be made available through dbGaP are raw FASTQ files.

  Study phs003343

• Center for Sub-Cellular Genomics

  A cell is a highly complex system with distributed molecular physiologies in structured sub-cellular compartments whose interplay with the nuclear genome determine the functional characteristics of the cell. A classic example of distributed genomic processes is found in neurons. Learning and memory require modulation of individual synapses through RNA localization, localized translation, and localized metabolites, such as those from dendritic mitochondria. Dendrites of neurons integrate distributed synaptic signals into both electrical and nuclear transcriptional response. Dysfunction of these distributed genomic functions in neurons can result in a broad spectrum of neuropsychiatric diseases such as bipolar and depressive disorders, autism, among others. The Center for Sub-Cellular Genomics' (CSG) mission is to develop new technologies to enable the dissection of complex genomic interactions within a single cell. New technologies include methods to target DNA and RNA sampling at sub-cellular resolution in live cells; high-throughput isolation and microfluidics chemistry for sub-cellular fractions; sub-cellular resolution mass-spectrometry; and computational tools for analysis of sub-cellular genomics data. The new technologies developed by CSG will enable the broader genomics community to associate sub-cellular processes to cell function and dysfunction. These technologies will be utilized in the center to generate a comprehensive atlas of genomic signatures of sub-cellular compartments in the mammalian neuron system. The technologies and neuron-focused data obtained here will be relevant to neurodegenerative and neuropsychiatric diseases including schizophrenia, depression, autism spectrum disorder, traumatic brain injury, and Alzheimer's.

  Study phs002120

• Machine learning to detect the SINEs of cancer

  We previously described an approach called RealSeqS to evaluate aneuploidy in plasma cell-free DNA (cfDNA) through the amplification of ~350,000 repeated elements with a single primer. We hypothesized that an unbiased evaluation of the large amount of sequencing data obtained with RealSeqS might reveal other differences between plasma samples from patients with and without cancer. This hypothesis was tested through the development of a novel machine-learning approach called Alu Profile Learning Using Sequencing (A-PLUS) and its application to samples from 5108 individuals, 2037 with cancer and the remainder without cancer. Samples from cancer patients and controls were pre-specified into four cohorts used for: 1) model training, 2) analyte integration and threshold determination, 3) validation, and 4) reproducibility. A-PLUS alone provided a sensitivity of 40.5% across 11 different cancer types in the Validation Cohort, at a specificity of 98.5%. Combining A-PLUS with aneuploidy and 8 common protein biomarkers detected 51% of 1167 cancers at 98.9% specificity. We found that part of the power of A-PLUS could be ascribed to a single feature – the global reduction of AluS sub-family elements in the circulating DNA of cancer patients. We confirmed this reduction through the analysis of another independent dataset obtained with a very different approach (whole genome sequencing). The evaluation of Alu elements therefore has the potential to enhance the performance of several methods designed for the earlier detection of cancer.

  Study EGAS00001007169

• Uncovering Inversion Formation in the Human Genome and its Impact to Disease

  The relevance of inversions for disease causation, speciation and adaptation is broadly recognized, although the prevalence of inversions is unknown. In humans, de novo inversions are associated with congenital anomalies in approximately 9.6% of patients. Yet, despite the biological relevance of inversions, their molecular features, formation mechanism, impact to the genomic structure of carriers, as well as their contribution to clinical phenotypes, have not been further explored. Inversions are typically classified as balanced reciprocal events generated by ectopic recombination, although recent studies reveal a distinct picture whereby inversions originate from mechanisms that concomitantly generate copy number variants (CNVs). Surprisingly, those complex inversions underlie as much as 30% of neurodevelopmental-associated CNVs. The hypothesis of this project are: inversions are often generated de novo by mechanisms other than ectopic recombination a relevant fraction of inversions are associated with complex genomic rearrangements (CGRs), often overlooked in sporadic diseases, and inversions are a "hidden" type of structural variation for which contribution to a clinical phenotypes has been under assessed due to the lack of appropriate detection tools. A combined strategy using multiple genomic tools will be applied to characterize inversions and associated genomic alterations in individuals with neurodevelopmental disorders and controls: whole genome sequencing (WGS) with short-and long-reads; genome mapping; classical cytogenetics; and array CGH. This project will establish common ground to bridge studies of rare and common diseases, human evolution, and population genetics.

  Study phs002999

• Biopsy-Derived Organoids in Personalised Early Breast Cancer Care: Challenges of Tumour Purity and Normal Cell Overgrowth Cap Their Practical Utility

  The ability to establish organoids composed exclusively of tumour rather than healthy cells is essential for their implementation into clinical practice. Organoids have in recent years emerged as a powerful tool to expand patient material in culture and generate modifiable 3D models derived from humans or animal models. For translational research, they enable the creation of model systems for an ever-increasing number of cell types and diseases. And in personalised medicine, they potentially allow for functional drug testing with high predictive power in certain settings. We found that using biopsy material from untreated, early-stage primary breast cancer patients poses significant challenges for consistently culturing tumour cells as organoids. Specifically, we observed frequent outgrowth of genetically normal, non-cancerous epithelial cells. We analysed >100 biopsy samples from early-stage breast cancer and present our large collection of >70 organoid lines. We also show methods of assessing successful tumour cell culture in a time and cost-efficient manner, proving a high rate (>85%) of normal cell overgrowth in early-stage breast cancer organoids. Finally, we show a number of successful attempts to culture cancer organoids from mastectomy-derived tissue of advanced, metastatic breast cancer. We conclude that the usefulness of organoids from early breast cancer for translational research and personalised medicine, especially guidance of adjuvant or post-surgical maintenance therapy, is strongly limited by the low success rate of culturing cancerous cells under organoid conditions.

  Study EGAS50000000605

• Decoding Human Endogenous Retrovirus Expression in Liver Metastatic Colorectal Cancers: Implications for Diagnosis and Prognosis

  The association between Human Endogenous RetroVirus (HERV) expression and cancer can be considered to both reflect the epigenetic deregulation of cancer and facilitate oncogenesis through its promoter/enhancer action on genes. More than 50% of colorectal cancers develop liver metastases and HERV has never been studied in this setting. We collected 400 RNAseq samples from more than 200 patients with primary and liver metastases from published and unpublished databases. We corroborated our findings with tandem mass spectrometry from primary colorectal cancer and single-cell RNAseq from liver metastatic colorectal cancer (lmCRC). We identified global stability in HERV expression between liver metastases and primary colorectal cancers, suggesting an early oncogenic footprint. Therefore, we proposed a list of seventeen HERV loci for lmCRC characterization with ERV-derived peptides but no expression in normal tissues. Among them, we discovered six prognostic markers after surgical resection of lmCRC. Four of them, HERVH_Xp22.32a, HERVH_20p11.23b, HERVH_13q33.3, HERVH_13q31.3, showed poor prognostic value when expressed (log-rank p-value 0.028, 0.0083, 9e-4, 0.05, respectively). HERVH_Xp22.2c (log-rank p-value 0.032) and HERVH_8q21.3b (only in multivariate models) were associated with better prognosis when expressed. These markers showed a cumulative effect on survival when expressed. Some HERVs were associated with decreased cytotoxic immune cells and most of them correlated with cell cycle pathways. Our findings provide new insights into lmCRC classification and suggest novel prognostic biomarkers and potential therapeutic targets.

  Study EGAS50000000307

• Cancer-associated fibroblasts promote drug resistance in ALK-driven lung adenocarcinoma cells by upregulating lipid biosynthesis

  Targeted therapy interventions using tyrosine kinase inhibitors (TKIs) provide encouraging treatment responses in ALK-rearranged lung adenocarcinomas, yet resistances occur almost inevitably. Apart from tumor cell-intrinsic resistance mechanisms, accumulating evidence supports a role of cancer-associated fibroblasts (CAFs) in affecting the therapeutic vulnerability of lung cancer cells. Here, we aimed to investigate underlying molecular networks shaping the therapeutic susceptibility of ALK-driven lung adenocarcinoma cells via tumor microenvironmental cues using three-dimensional (3D) spheroid co-culture settings. We show that CAFs promote therapy resistance of lung tumor cells against ALK inhibition by reducing apoptotic cell death and increasing cell proliferation. Using single-cell RNA-sequencing analysis, we show that genes involved in lipogenesis constitute the major transcriptional difference between TKI-treated homo- and heterotypic lung tumor spheroids. CAF-conditioned medium and CAF-secreted factors HGF and NRG1 were both able to promote resistance of 3D-cultured ALK-rearranged lung tumor cells via AKT signaling, which was accompanied by enhanced de novo lipogenesis and supression of lipid peroxidation. Notably, simultaneous targeting of ALK and SREBP-1 was able to overcome the established CAF-driven lipid metabolic-supportive niche of TKI-resistant lung tumor spheroids. Our findings highlight a crucial role of CAFs in mediating ALK-TKI resistance via lipid metabolic reprogramming and suggest new ways to overcome resistance towards molecular directed drugs by targeting vulnerabilities downstream of oncogenic signaling.

  Study EGAS50000000135

• Enhanced Adjuvanticity of a Personal Neoantigen Vaccine Generates Potent Neoantigen-Specific Immunity

  Personalized neoantigen-targeting vaccines can induce durable and diverse T cell responses and have been associated with improved outcomes. To augment immunogenicity, we tested a synthetic long peptide vaccine formulated with the adjuvants Montanide and poly-ICLC in 10 patients with melanoma. The vaccines, co-administered locally with ipilimumab and systemically with nivolumab, generated de novo ex vivo T cell responses against the majority of immunizing neoepitopes in all 9 fully vaccinated patients, and ex vivo CD8+ T cell responses in 6 of 9. Hundreds of circulating and intratumoral T cell receptor (TCR) clonotypes emerged following vaccination that were distinct from those arising after PD-1 inhibition. By linking vaccine-specificity and cellular phenotype of TCRs in post-treatment tumors at single cell resolution, we demonstrate remodeling of the intratumoral T cell repertoire following vaccination. These observations demonstrate the impact of multi-pronged immune adjuvanticity on the T cell repertoire in the context of neoantigen-targeting vaccines. We are submitting WES data from whole blood samples and tumor FFPE samples, and RNA-seq data from tumor FFPE samples. We are also submitting RNA-seq data from purified tumor samples. scRNA-seq data will also be submitted for skin samples (pre and post neoantigen vaccine dose #3), and tumor samples before starting treatment, before starting neoantigen vaccine Dose #1, and before neoantigen vaccine dose #4. Bulk TCR-seq data will be submitted for blood samples.

  Study phs003919

• 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

• HuBMAP: A Spatially Resolved Molecular Atlas of Human Endothelium

  The human circulatory system plays a vital role in the proper functioning of every organ in the human body. Because vascular related diseases, such as atherosclerosis, do not occur homogeneously throughout the vascular tree, differences in the endothelium of various organs must exist. There is even evidence of heterogeneity of adjacent endothelial cells. We are studying normal endothelium throughout the body and across individuals of different ages to characterize these differences. The first step will be to procure normal tissue samples from a wide range of organs. From disaggregated samples of these tissues, we will perform single cell combinatorial indexing (sci-) to generate chromatin accessibility (sci-ATAC-seq) and transcriptome (sci-RNA-seq) data from millions of single cells across hundreds of samples.  Molecular heterogeneity might correlate with multiple facets of the anatomical distribution of endothelium, including vascular vs. lymphatic, arteries vs. veins, capillaries vs. larger vessels, differences between organs, organ zonation, cell-to-cell differences, etc. Based on heterogeneous markers defined from deep profiling of disaggregated cells, we will employ seqFISH to spatially relate the distinct endothelial subpopulations arising from the single nuclei data to histologic preparations. For example, this will allow us to identify subpopulations with particular markers as arterial, capillary, venous, or lymphatic. The end product will include multiparameter images in which the complete gene expression and chromatin accessibility profile can be interactively visualized for every cell.

  Study phs002267

• Genomic landscapes of endometrioid and mucinous ovarian cancers and morphologically similar tumor types

  While endometrioid and mucinous ovarian cancers represent nearly a quarter of ovarian cancers, their molecular characteristics and pathologic origins are poorly understood. To characterize the genomic and epigenomic alterations of these ovarian cancer subtypes and evaluate links to morphologically similar tumors from other sites, we performed sequence, copy number, mutation signature, rearrangement, and methylation analyses from tumor samples and matched normal tissues of 133 patients. These analyses included patients with ovarian endometrioid (n=44), ovarian mucinous (n=43), uterine endometrial (n=15), and gastrointestinal mucinous cancers (n=31). In addition to identifying genes previously known to be involved in these tumors, we identified alterations in RAD51C, NOTCH4, SMARCA1/4 and JAK1 in ovarian endometrioid, ESR1 in uterine endometrioid, and SMARCA4 in ovarian mucinous cancers. Whole genome sequencing revealed rearrangements involving PTEN, NF1, and NF2 in ovarian endometrioid tumors, and NF1 and MED1 in ovarian mucinous cancers. The number of alterations and affected genes in ovarian and uterine endometrioid cancers were similar as were genome-wide methylation profiles, suggesting that these tumors may have a similar tissue of origin. By contrast, methylation patterns in ovarian mucinous cancers were more similar to other gastrointestinal mucinous cancers, including mucinous cancers of the stomach, colon, and pancreas. These analyses provide insights into the genomic landscapes and origins of mucinous and endometrioid ovarian carcinomas, providing new avenues for early clinical intervention and management of patients with these cancers.

  Study EGAS50000000523

• CRISPR screening in human trophoblast stem cells reveals both shared and distinct aspects of human and mouse placental development

  The placenta serves as the interface between the mother and fetus, facilitating the exchange of gases and nutrients between their separate blood circulation systems. Trophoblasts in the placenta play a central role in this process. Our current understanding of mammalian trophoblast development relies largely on mouse models. However, given the diversification of mammalian placentas, findings from the mouse placenta cannot be readily extrapolated to other mammalian species, including humans. To fill this knowledge gap, we performed CRISPR knockout (KO) screening in human trophoblast stem cells (hTSCs). We targeted genes essential for mouse placental development and identified more than 100 genes as critical regulators in both human hTSCs and mouse placentas. Among them, we further characterized in detail two transcription factors, DLX3 and GCM1, and revealed their essential roles in hTSC differentiation. Moreover, a gene function-based comparison between human and mouse trophoblast subtypes suggests that their relationship may differ significantly from previous assumptions based on tissue localization or cellular function. Notably, our data reveal that hTSCs may not be analogous to mouse TSCs or the extraembryonic ectoderm (ExE) in which in vivo TSCs reside. Instead, hTSCs may be analogous to progenitor cells in the mouse ectoplacental cone and chorion. This finding is consistent with the absence of ExE-like structures during human placental development. Our data not only deepen our understanding of human trophoblast development but also facilitate cross-species comparison of mammalian placentas.

  Study JGAS000659

• Genetic alterations that deregulate RB and PDGFRA signaling pathways drive tumor progression in IDH2-mutant astrocytoma

  In IDH-mutant astrocytomas, IDH2 mutation is quite rare and biological mechanisms underlying tumor progression in IDH2-mutant astrocytoma remain elusive. Here, we report a unique case of IDH2 mutant astrocytoma, CNS WHO grade 3 that developed tumor progression. We performed a comprehensive genomic and epigenomic analysis for primary and recurrent tumors and found that both tumors harbored recurrent IDH2R172K and TP53R248W mutation with CDKN2A/B hemizygous deletion. We also found amplifications of CDK4 and MDM2 with PDGFRA gain in the recurrent tumor and upregulated protein expressions of these genes. We further developed, for the first time, a xenograft mouse model of IDH2R172K mutant astrocytoma from the recurrent tumor, but not from the primary tumor. Consistent with parent recurrent tumor cells, amplifications of CDK4 and MDM2 with PDGFRA gain was found, while CDKN2A/B was identified as homozygous deletion in the xenografts, qualifying for integrated diagnosis of astrocytoma, IDH2-mutant, CNS WHO grade 4. Cell viability assay found that CDK4/6 inhibitor and PDGFR inhibitor potently decreased cell viability in recurrent tumor cells, as compared to primary tumor cells. These findings suggest that gene alterations that activate retinoblastoma (RB) signaling pathways and PDGFR may drive tumor progression and xenograft formation in IDH2-mutant astrocytoma, which is equivalent to progressive IDH1-mutant astrocytoma. Also, our findings suggest that these genomic alterations may represent therapeutic targets in IDH2-mutant astrocytoma.

  Study JGAS000646

• Homozygous inactivation of CHEK2 is linked to a familial case of multiple primary lung cancer with accompanying cancers in other organs

  In clinical practice, there are a number of cancer patients with clear family histories, but not belonging to existing familial malignancies. Recent advances in genomic technologies have enhanced the possibility of identifying causative genes in such cases. Two siblings, an elder sister and a younger brother, were found to have multiple primary lung cancers at the age of 60. The former subsequently developed breast cancer, and had a history of myoma uteri. The latter had initially developed prostate cancer at the age of 59, and had a history of colon cancer. SNP genotyping revealed that approximately 10% of the genomes were homozygous in both patients. Exome sequencing revealed non-synonymous mutations in four genes in the runs of homozygosity: CHEK2, FCGRT, INPP5J, and SFI1. Evolutionary conservation of primary protein structures suggested the functional importance of the CHEK2 mutation, R474C. This mutation altered the tertiary structure of CHK2 by disrupting the salt bridge between R474 and E394. No such structural changes were observed with the other three mutated genes. Subsequent cell-based transfection analysis revealed that CHK2 R474C was unstable and scarcely activated. We concluded that CHEK2 is a causative gene in this case of familial cancer. Although homozygous inactivation of CHEK2 in mice led to cancers in multiple organs, accumulation of additional human cases is needed to establish its pathogenic role in humans.

  Study JGAS000057

• DNA methylation at HBV integrants and flanking host genomes

  Integration of DNA viruses into the human genome plays an important role in various types of tumors, including hepatitis B virus (HBV)-related hepatocellular carcinoma. However, the molecular details and clinical impact of HBV integration on either the human or HBV epigenomes are unknown. Here, we show that methylation of the integrated HBV DNA is related to the methylation status of the flanking human genome. We developed a next-generation sequencing-based method for structural methylation analysis of integrated viral genomes (denoted G-NaVI). This method is a novel approach that enables enrichment of viral fragments for sequencing using unique baits based on the sequence of the HBV genome. We detected integrated HBV sequences in the genome of the PLC/PRF/5 cell line and found variable levels of methylation within the integrated HBV genomes. Allele-specific methylation analysis revealed that the HBV genome often became significantly methylated when integrated into highly methylated host sites. After integration into unmethylated human genome regions such as promoters, however, the HBV DNA remains unmethylated and may eventually play an important role in tumorigenesis. The observed dynamic changes in DNA methylation of the host and viral genomes may functionally affect the biological behavior of HBV. These findings may impact public health given that millions of people worldwide are carriers of HBV. We also believe our assay will be a powerful tool to increase our understanding of the various types of DNA virus-associated tumorigenesis.

  Study JGAS000015

• Identification_of_synthetic_lethal_genes_by_CRISPR_Cas9_library

  Sequencing is performed on PCR products from the amplification of gRNA segments introduced into cell lines (some are commercial, some are derived by collaborators, as described in HMDMC 13/076). Cellular DNA or RNA are not sequenced. This data is defined as "Other Biological/Biochemical Assay Data" under the WTSI datasharing guidelines and consequently will be released at the time of publication and/or through faculty websites. SEQUENCING DETAILS: in October 2016 3 DNA samples will be submitted for PCR amplification and purification through CGP pipeline using Thompson's custom primer set and PCR settings and sequencing with Hiseq using a custom protocol (with 3 forward reads to read sgRNA1, sgRNA2 and the index, optimized with Illumina Bespoke Sanger team); read length 50bp. The following set of samples will be submitted in later 2016 and ealry 2017 as cell pellet to be processed through CGP Cas9 pipeline (Thompson's oligos and protocols) and sequences on Hiseq (same custom sequencing protocol described above) We are utilizing a custom CRISPR/Cas9 library to identify synthetic lethal gene pairs. Cells will be transduced with LV expressing Cas9, then infected with LV expressign the library of 42'000 sgRNA pairs. Cell pellet will be collected from the cell population at different timepoints. The samples will be processed by CGP pipeline customized for the specific protocol and then sequenced to identify the 19bp of the sgRNA and quantify their representation.

  Study EGAS00001002117

• Identification of a Fumarate-Hydratase Deficient-like RCC Subtype with Convergent Pathway Alterations

  The development of molecularly targeted therapies builds upon the identification of tumor clades with shared molecular abnormalities. Hereditary leiomyomatosis and renal cell carcinoma (HLRCC) is a familial syndrome resulting from germline mutation of the two-hit tumor suppressor gene, Fumarate Hydratase (FH). HLRCC renal tumors present often in young adults, exhibit unique histological features and are particularly aggressive. HLRCC is uncommon and poorly understood, and whether there are sporadic renal tumors with similar characteristics, is unclear. Notably, FH can also be inactivated somatically in RCC, and should be included among FH-deficient (FHD) kidney tumors. We identified a tumor clade without FH mutations but with similar morphological and molecular characteristics as well as aggressive clinical behavior, which we name FHDL (FHD-like). We show that FHD/FHDL tumors represent a unique subtype of CIMP (CpG Island Methylator Phenotype) with the highest levels of overall genome methylation and CDKN2A silencing. FHD/FHDL tumors are characterized by convergent activation of the Hippo pathway as determined by frequent and mutually exclusive mutations in NF2 and PTPN14. AKR1B10, an NRF2 target, was upregulated in FHD/FHDL tumors as well as GATA3, which suggests a shared origin in the renal medulla and is consistent with topological studies. There is loss of 5-hydroxymethylcytosine in both FHD and FHDL tumors. These data define a group of particularly aggressive renal cancers with shared molecular features and potential targets for therapeutic intervention.

  Study EGAS00001002646

• 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

• Negligible impact on missing heritability of autoimmune-locus rare coding-region variants

  Genome wide association studies (GWAS) have identified common variants of modest effect size at hundreds of loci for common autoimmune diseases - however a substantial fraction of heritability remains unexplained, and to which rare variants may contribute. To discover rare variants and test them for association with a phenotype, the majority of studies re-sequence a small initial sample size and then genotype the discovered variants in a larger sample set. This approach will fail to analyse a large fraction of the rare variants present in the entire sample set. Here we perform simultaneous amplicon sequencing-based variant discovery and genotyping for coding exons of 25 GWAS risk genes in 41,911 white-European origin UK subjects comprising 24,892 subjects with six autoimmune disease phenotypes and 17,019 controls, and show that rare coding-region variants at known loci play a negligible role in common autoimmune disease susceptibility. These results do not support the rare variant - synthetic genome-wide association hypothesis. Many known autoimmune disease risk loci contain multiple independent common and low frequency variant association signals, and so genes in these loci are a priori stronger candidates to harbor rare coding-region variants than other genes. Our data suggest that the missing heritability for common autoimmune diseases may not lie in the rare coding-region variant portion of the allelic spectrum, but perhaps as others have proposed in very many common variant signals of weak effect.

  Study EGAS00001000476

• Cisplatin increases sensitivity to FGFR inhibition in patient-derived xenograft models of lung squamous cell carcinoma

  Lung squamous cell carcinoma (SqCC) is a molecularly complex and genomically unstable disease. No targeted therapy is currently approved for lung SqCC, although potential oncogenic drivers of SqCC have been identified, including amplification of the fibroblast growth factor receptor 1 (FGFR1). Reports from a recently completed clinical trial indicate low response rates in patients treated with FGFR tyrosine kinase inhibitors, suggesting inadequacy of FGFR1 amplification as a biomarker of response, or the need for combination treatment. We aimed to develop accurate models of lung SqCC and determine improved targeted therapies for these tumors. We show that detection of FGFR1 mRNA by RNA in situ hybridisation is a better predictor of response to FGFR inhibition than FGFR1 gene amplification using clinically relevant patient-derived xenograft (PDX) models of lung SqCC. FGFR1-overexpressing tumors were observed in all histological subtypes of non-small cell lung cancers (NSCLC) as assessed on a tissue microarray, indicating a broader range of tumors that may respond to FGFR inhibitors. In FGFR1-overexpressing PDX tumors, we observed increased differentiation and reduced proliferation following FGFR inhibition. Combination therapy with cisplatin was able to increase tumor cell death, and dramatically prolonged animal survival compared to single agent treatment. Our data suggest that FGFR tyrosine kinase inhibitors can benefit NSCLC patients with FGFR1-overexpressing tumors and provides a rationale for clinical trials combining cisplatin with FGFR inhibitors.

  Study EGAS00001002423

• Evolution of the African pygmy phenotype

  The evolutionary history of the human ‘pygmy’ phenotype (small body size), a characteristic of African and Southeast Asian rainforest hunter-gatherers, is largely unknown. Here we use a genome-wide admixture mapping analysis to identify 16 genomic regions that are significantly associated with the pygmy phenotype in the Batwa, a rainforest hunter-gatherer population from Uganda (East Central Africa). The identified genomic regions have multiple attributes that provide supporting evidence of genuine association with the pygmy phenotype, including enrichments for SNPs previously associated with stature variation in Europeans and for genes with growth hormone receptor and regulation functions. To test adaptive evolutionary hypotheses we computed the haplotype-based iHS statistic and the level of population differentiation (FST) between the Batwa and their agricultural neighbors, the Bakiga, for each SNP. Both |iHS| and FST values were significantly higher for SNPs within the Batwa pygmy phenotype-associated regions than the remainder of the genome, a signature of polygenic adaptation. In contrast, when we expanded our analysis to include Baka rainforest hunter-gatherers from Cameroon and Gabon (West Central Africa) and Nzebi and Nzime neighboring agriculturalists, we did not observe elevated |iHS| or FST values in these genomic regions. Together, these results suggest adaptive and at least partially convergent origins of the pygmy phenotype even within Africa, supporting the hypothesis that small body size confers a selective advantage for tropical rainforest hunter-gatherers but raising questions about the antiquity of this behavior.

  Study EGAS00001000908

• Bioinformatic Methods and Bridging of Assay Results for Reliable Tumor Mutational Burden Assessment in Non-Small Cell Lung Cancer

  Tumor mutational burden (TMB) has emerged as a predictive biomarker of response to immune checkpoint inhibitors. Standardization of TMB measurement is essential for implementing diagnostic tools to guide treatment. Here we evaluate bioinformatic TMB analysis by whole exome sequencing (WES) in formalin-fixed, paraffin-embedded samples. In CheckMate 026, TMB was retrospectively assessed in 312 patients with non-small cell lung cancer (58% of the intent-to-treat population) who received first-line nivolumab treatment or chemotherapy. We examined the sensitivity of TMB assessment to bioinformatic filtering methods and assessed concordance between TMB data derived by WES and the FoundationOne CDx™ assay. TMB scores comprising synonymous, indel, frameshift, and nonsense mutations (all mutations) were 3.1-fold higher than data including missense mutations only, but values were highly correlated (Spearman’s r = 0.99). Scores from CheckMate 026 samples including missense mutations only were similar to those generated from data in The Cancer Genome Atlas, but those including all mutations were generally higher. Using databases for germline subtraction (instead of matched controls) showed a trend for race-dependent increases in TMB scores. Parameter variation can therefore impact TMB calculations, highlighting the need for standardization. Encouragingly, WES and FoundationOne CDx outputs were highly correlated (Spearman’s r = 0.90) and differences could be accounted for by empirical calibration, suggesting that reliable TMB assessment across assays, platforms and centers is achievable.

  Study EGAS00001003661

• Whole exome sequencing of an invasive diffuse intrinsic pontine glioma sampled from different sites

  The lateral ventricle subventricular zone (SVZ) is a frequent site of high-grade glioma (HGG) spread and tumor invasion of the SVZ is a predictor of a worse clinical prognosis. A range of HGG types invade the stem cell niche of the SVZ, including both adult glioblastoma and pediatric high-grade gliomas such as diffuse intrinsic pontine glioma (DIPG). The cellular and molecular mechanisms mediating this frequent invasion of the SVZ are poorly understood. Here we demonstrate that neural precursor cell (NPC):glioma cell communication underpins the propensity of glioma to colonize the SVZ. SVZ NPCs secrete chemoattractant signals toward which glioma cells home. Biochemical and proteomic analyses of factors secreted by SVZ NPCs revealed a discrete list of candidate proteins. Necessity and sufficiency testing implicated the neurite outgrowth-promoting factor pleiotrophin, along with three required binding partners (secreted protein acidic and rich in cysteine (SPARC), SPARC-like protein 1, and heat shock protein 90B), as key mediators of this chemoattractant effect. Pleiotrophin protein expression is strongly enriched in the SVZ, and knockdown of pleiotrophin expression starkly reduced glioma invasion of the SVZ in the adult murine brain. Pleiotrophin, in complex with the three binding partners, activated the Rho/ROCK pathway in DIPG cells, and inhibition of ROCK resulted in decreased DIPG invasion toward SVZ neural precursor cell-secreted factors. These findings demonstrate a pathogenic role for neural precursor cell – glioma interactions and potential therapeutic targets to limit glioma invasion.

  Study EGAS00001002326

• Successful BRAF/MEK-inhibition in a young patient with BRAF V600E-mutated extrapancreatic acinar cell carcinoma (HIPO-021)

  Pancreatic acinar cell carcinoma (PAC) is a rare disease with a poor prognosis. Treatment options for metastatic PAC are limited and often follow chemotherapeutic regimens for pancreatic ductal adenocarcinoma. Although recurrent genomic alterations, such as BRAF fusions and defects in genes involved in homologous recombination DNA repair, have been described in PAC, data on the clinical efficacy of molecularly guided, targeted treatment are scarce. Here we describe the case of a 27-yr-old patient with BRAF V600E-mutated PAC who was successfully treated with a combination of BRAF and MEK inhibitors. The patient presented to our clinic with abdominal pain and weight loss. Imaging showed extensive retroperitoneal disease as well as mediastinal lymphadenopathy. Because of elevated alpha-fetoprotein (AFP) levels and inconclusive histologic findings, a germ cell tumor was suspected; however, PEI chemotherapy was unsuccessful. A repeat biopsy yielded the diagnosis of PAC and treatment with FOLFIRINOX was initiated. Comprehensive molecular profiling within the MASTER (Molecularly Aided Stratification for Tumor Eradication Research) precision oncology program revealed a somatic BRAF V600E mutation and a germline PALB2 stop-gain mutation. Therapy was therefore switched to BRAF/MEK inhibition, resulting in almost complete remission and disease control for 12 mo and a remarkable improvement in the patient's general condition. These results indicate that BRAF alterations are a valid therapeutic target in PAC that should be routinely assessed in this patient population.

  Study EGAS00001004282

• Integration of genomics and histology reveals diagnosis and effective therapy of refractory cancer of unknown primary with PDL1 amplification (H021)

  Identification of the tissue of origin in cancer of unknown primary (CUP) poses a diagnostic challenge and is critical for directing site-specific therapy. Currently, clinical decision making in patients with CUP primarily relies on histopathology and clinical features. Comprehensive molecular profiling has the potential to contribute to diagnostic categorization and, most importantly, guide CUP therapy through identification of actionable lesions. We here report the case of an advanced-stage malignancy initially mimicking poorly differentiated soft-tissue sarcoma that did not respond to multi-agent chemotherapy. Molecular profiling within a clinical whole-exome and transcriptome sequencing program revealed a heterozygous, highly amplified KRAS G12S mutation, compound-heterozygous TP53 mutation/deletion, high mutational load, and focal amplification of chromosomes 9p (including PDL1 [CD274] and JAK2) and 10p (including GATA3). Integrated analysis of molecular data and conventional histopathology suggested a diagnosis of triple-negative breast cancer (TNBC) and provided a rationale for immune checkpoint inhibitor therapy with pembrolizumab, which resulted in rapid clinical improvement and a lasting partial remission. Analysis of 157 TNBC samples from The Cancer Genome Atlas revealed focal, high-level PDL1 amplification coinciding with excessive PDL1 mRNA expression in 24% of cases. Collectively, these results illustrate the diagnostic utility of multidimensional tumor profiling in cases with non-descript histology and immune phenotype, demonstrate the predictive power of genomic PDL1 amplification for immune checkpoint inhibition, and suggest a targeted therapeutic strategy in chromosome 9p24.1/PDL1-amplified cancers.H021

  Study EGAS00001001846

• Initial whole genome sequencing of plasma cell neoplasms in First Responders exposed to the World Trade Center attack of September 11, 2001

  The World Trade Center (WTC) attack of September 11, 2001 created an unprecedented environmental exposure to known and suspected carcinogens. High incidence of multiple myeloma (MM) and precursor conditions has been reported among first responders to the WTC disaster. To expand on our prior screening studies, and to characterize the genomic impact of the exposure to known and potential carcinogens in the WTC debris, we were motivated to perform whole genome sequencing (WGS) of WTC first responders and recovery workers who developed a plasma cell disorder after the attack. PATIENTS AND METHODS: We performed WGS of 9 CD138-positive bone marrow mononuclear samples from patients who were diagnosed with plasma cell disorders after the WTC disaster. RESULTS: No significant differences were observed in comparing the post-WTC driver and mutational signatures landscape with 110 previously published WGS from 56 patients with MM and the CoMMpass WGS cohort (n=752). Leveraging constant activity of the single base substitution mutational signatures 1 and 5 over time, we estimated that tumor-initiating chromosomal gains were windowed to both pre- and post-WTC exposure. CONCLUSIONS: Although limitations in sample size preclude any definitive conclusions, our findings suggest that the observed increased incidence of plasma cell neoplasms in this population is due to complex and heterogeneous effects of the WTC exposure that may have initiated or contributed to progression of malignancy.

  Study EGAS00001004467

• 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

• Massively parallel nanowell-based single-cell gene expression profiling

  Technological advances have enabled transcriptome characterization of cell types at the single-cell level providing new biological insights. New methods that enable simple yet high-throughput single-cell expression profiling are highly desirable. Here we report a novel nanowell-based single-cell RNA sequencing system, ICELL8, which enables processing of thousands of cells per sample. The system employs a 5184-nanowell-containing microchip to capture ~1300 single cells and process them. Each nanowell contains preprinted oligonucleotides encoding poly-d(T), a unique well barcode, and a unique molecular identifier. The ICELL8 system uses imaging software to identify nanowells containing viable single cells and only wells with single cells are processed into sequencing libraries. Here, we report the performance and utility of ICELL8 using samples of increasing complexity from cultured cells to mouse solid tissue samples. Our assessment of the system to discriminate between mixed human and mouse cells showed that ICELL8 has a low cell multiplet rate (< 3%) and low cross-cell contamination. We characterized single-cell transcriptomes of more than a thousand cultured human and mouse cells as well as 468 mouse pancreatic islets cells. We were able to identify distinct cell types in pancreatic islets, including alpha, beta, delta and gamma cells. Overall, ICELL8 provides efficient and cost-effective single-cell expression profiling of thousands of cells, allowing researchers to decipher single-cell transcriptomes within complex biological samples.

  Study EGAS00001002320

• Tissue-specific mutation accumulation in human adult stem cells during life

  Gradual accumulation of mutations in human adult stem cells during life is associated with various age-related diseases, including cancer. The number of stem cell divisions throughout life is believed to be a major determinant for mutation accumulation and could explain the extreme variation of cancer incidence across different organs. Yet, mutation patterns and rates of healthy adult stem cells remain unknown. Here, we determined genome-wide mutation patterns in primary adult stem cells of the small intestine, colon and liver of human donors with ages ranging from 3 to 87 years. We find that the number of mutations increases linearly with age up to several thousand mutations per cell at 87 years of age, while mutation spectra remain constant throughout life. Small intestine and colon stem cells have a 2-fold higher mutation rate per year compared with liver stem cells. These differences could be exclusively attributed to the mutagenic action of spontaneous deamination of cytosine residues and may reflect the high stem cell division rate in these tissues. The genomic distribution of somatic mutations is non-random and predominantly associated with DNA replication dynamics in the small intestine and colon, and with transcription in the liver. These results indicate that a stable balance between various mutagenic and DNA repair processes is maintained throughout life and that the activity of these processes in adult stem cells varies between tissues.

  Study EGAS00001001682

• The evolutionary history of human colitis-associated colorectal cancer

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

  Study EGAS00001003028

• Genome-wide analysis of genetic risk factors for rheumatic heart disease in Aboriginal Australians provides support for pathogenic molecular mimicry

  Background. Rheumatic heart disease (RHD) following Group A Streptococcus (GAS) infections is heritable and prevalent in Indigenous populations. Molecular mimicry between human and GAS proteins triggers pro-inflammatory cardiac valve-reactive T-cells. Methods. Genome-wide genetic analysis was undertaken in 1263 Aboriginal Australians (398 RHD cases; 865 controls). Single nucleotide polymorphisms (SNPs) were genotyped using Illumina HumanCoreExome BeadChips. Direct typing and imputation was used to fine-map the human leukocyte antigen (HLA) region. Epitope binding affinities were mapped for human cross-reactive GAS proteins, including M5 and M6. Results. The strongest genetic association was intronic to HLA-DQA1 (rs9272622; P=1.86x10-7). Conditional analyses showed rs9272622 and/or DQA1*AA16 account for the HLA signal. HLA-DQA1*0101_DQB1*0503 (OR 1.44, 95%CI 1.09-1.90, P=9.56x10-3) and HLA-DQA1*0103_DQB1*0601 (OR 1.27, 95%CI 1.07-1.52, P=7.15x10-3) were risk haplotypes; HLA_DQA1*0301-DQB1*0402 (OR 0.30, 95%CI 0.14-0.65, P=2.36x10-3) was protective. Human myosin cross-reactive N-terminal and B repeat epitopes of GAS M5/M6 bind with higher affinity to DQA1/DQB1 alpha/beta dimers for the two risk haplotypes than the protective haplotype. Conclusions. Variation at HLA_DQA1-DQB1 is the major genetic risk factor for RHD in Aboriginal Australians studied here. Cross-reactive epitopes bind with higher affinity to alpha/beta dimers formed by risk haplotypes, supporting molecular mimicry as the key mechanism of RHD pathogenesis.

  Study EGAS00001002678

• Whole genome sequencing with linked reads of pediatric glioblastoma samples

  Pediatric glioblastoma (pGBM) is a lethal cancer with no effective therapies. To understand mechanisms of tumor evolution in this cancer, we performed whole genome sequencing with linked reads on longitudinally resected pGBM samples. Our analyses showed that all diagnostic and recurrent samples were collections of genetically diverse subclones. Clonal composition rapidly evolved at recurrence, with less than 8% of non-synonymous single nucleotide variants being shared in diagnostic-recurrent pairs. In order to track the origins of the mutational events we observed in pGBM, we generated whole genome datasets for two patients and their parents. These trios showed that genetic variants could be (i) somatic, (ii) inherited from a healthy parent, or (iii) arose de novo in the germlines of pGBM patients. Analysis of variant allele frequencies supported a model of tumor growth involving slow-cycling cancer stem cells that give rise to fast-proliferating progenitor-like cells and to non-dividing cells. Interestingly, radiation and anti-mitotic chemotherapeutics did not increase overall tumor burden upon recurrence. These findings support an important role for slow-cycling stem cell populations in contributing to recurrences, since slow-cycling cell populations are expected to be less prone to genotoxic stress induced by these treatments and to accumulate few mutations. Our results highlight the need for new targeted treatments that account for the complex functional hierarchies and genomic heterogeneity of pGBM.

  Study EGAS00001003432

• UK10K NEURO ASD GALLAGHER

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. This is an Irish sample set of individuals with ASD (approximately 50% with comorbid intellectual disability). Individuals have been diagnosed with ADI/ ADOS, measures of cognition/ adaptive function. They represent a more severe, narrowly defined cohort of ASD subjects. Family histories are available for some with measures of broader phenotype. For further information on this cohort please contact Nadia Bolshakova (bolshakn@tcd.ie).

  Study EGAS00001000112

• Probabilistic cell type assignment of single-cell transcriptomic data reveals spatiotemporal microenvironment dynamics in human cancers

  Single-cell RNA sequencing (scRNA-seq) has transformed biomedical research, enabling decomposition of complex tissues into disaggregated, functionally distinct cell types. For many applications, investigators wish to identify cell types with known marker genes. Typically, such cell type assignments are performed through unsupervised clustering followed by manual annotation based on these marker genes, or via "mapping" procedures to existing data. However, the manual interpretation required in the former case scales poorly to large datasets, which are also often prone to batch effects, while existing data for purified cell types must be available for the latter. Furthermore, unsupervised clustering can be error-prone, leading to under- and over- clustering of the cell types of interest. To overcome these issues we present CellAssign, a probabilistic model that leverages prior knowledge of cell type marker genes to annotate scRNA-seq data into pre-defined and de novo cell types. CellAssign automates the process of assigning cells in a highly scalable manner across large datasets while simultaneously controlling for batch and patient effects. We demonstrate the analytical advantages of CellAssign through extensive simulations and exemplify real-world utility to profile the spatial dynamics of high-grade serous ovarian cancer and the temporal dynamics of follicular lymphoma. Our analysis reveals subclonal malignant phenotypes and points towards an evolutionary interplay between immune and cancer cell populations with cancer cells escaping immune recognition.

  Study EGAS00001003452

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

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

  Study EGAS00001003547

• Cistrome-partitioning reveals convergence of somatic mutations and risk-variants on master transcription regulators in primary prostate tumors

  Thousands of noncoding somatic Single Nucleotide Variants (SNVs) of unknown function are reported in tumors. Partitioning the genome according to cistromes, reveals the enrichment of somatic SNVs in prostate tumors as opposed to adjacent normal tissue cistromes of master transcription regulators, including AR, FOXA1 and HOXB13. This parallels enrichment of prostate cancer genetic predispositions over these transcription regulators’ tumor cistromes, exemplified at the 8q24 locus harboring both risk-variants and somatic SNVs in cis-regulatory elements, upregulating MYC expression and altering the binding of transcription regulators to DNA. However, Massively-Parallel Reporter Assays reveal that few SNVs can alter the transactivation potential of individual CREs. Instead, SNVs accumulate, similarly to inherited riskvariants, in cistromes of master transcription regulators required for prostate cancer development. Difficulties in inferring the biological significance of noncoding mutations have limited their inclusion in precision genomics medicine pipelines. Most attempts to delineate a role for noncoding mutations relied on detecting evidence for positive selection within individual CREs, such as reported for the TERT gene promoter. By considering the enrichment of noncoding mutations in cistromes as opposed to individual CREs, we reveal their specificity towards master transcription regulators that promote prostate cancer development, a feature shared with inherited risk-variants. Overall, our work provides a blueprint for the functional interpretation of noncoding mutations in genomic tests relying on defining cis-regulatory units according to cistrome-partitioning to identify cancer driver transcription regulators.

  Study EGAS00001003928

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

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

  Study EGAS00001000249

• Open_Targets_020_Epigenomes_of_Cell_Lines

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

  Study EGAS00001003138

• Long-read trio sequencing of unsolved patients with intellectual disability

  Long-read sequencing (LRS) has the potential to comprehensively identify all medically-relevant genome variation, including variation commonly missed by short read sequencing (SRS) approaches. To determine this potential, we performed LRS up to 40x coverage using the Pacific Biosciences Sequel instrument for 5 patient-parent trios with intellectual disability, whose aetiology remained unresolved after SRS exomes and genomes. A first systematic assessment of LRS coverage showed that ~35Mb of the genome was only accessible by LRS and not SRS. Structural variant (SV) calling yielded on average 28,292 SV calls per individual, totalling 12.9 MB of sequence. Trio-based analyses showed concordance for up to 95% of these SV calls across the genome, and 80% of the LRS SV calls were not identified by SRS. De novo mutation analysis did not identify any de novo SVs, confirming these are rare events. Because of high sequence coverage, we were also able to call single nucleotide substitutions. On average, we identified 3 million substitutions per genome, with a Mendelian inheritance concordance of up to 97%. Of these, approximately 100,000 were located in the ~35Mb of the genome that was only captured by LRS. Moreover, these variants affected the coding sequence of 64 genes, including 32 known to cause Mendelian disorders. Although no clinical diagnoses were made in these 5 patients, our data show the potential added value of LRS compared to SRS for identifying medically relevant genome variation.

  Study EGAS00001004319

• Targets of MEK inhibition in DIPG

  The survival of children with diffuse intrinsic pontine glioma (DIPG) remains dismal, with new treatments desperately needed. In the era of precision medicine, targeted therapies represent an exciting treatment opportunity, yet resistance can rapidly emerge, playing an important role in treatment failure. In a prospective biopsy-stratified clinical trial, we combined detailed molecular profiling (methylation BeadArray, exome, RNAseq, phospho-proteomics) linked to drug screening in newly-established patient-derived models of DIPG in vitro and in vivo. We identified a high degree of in vitro sensitivity to the MEK inhibitor trametinib in samples which harboured genetic alterations targeting the MAPK pathway, including the non-canonical BRAF_G469V mutation, and those affecting PIK3R1 and NF1. However, treatment of PDX models and of a patient with trametinib at relapse failed to elicit a significant response. We generated trametinib-resistant clones in the BRAF_G469V model through continuous drug exposure, and identified acquired mutations in MEK1/2 (MEK1_K57N, MEK1_I141S and MEK2_I115N) with sustained pathway up-regulation. These cells showed the hallmarks of mesenchymal transition, and expression signatures overlapping with inherently trametinib-insensitive primary patient-derived cells that predicted an observed sensitivity to dasatinib. Combinations of trametinib with dasatinib and the downstream ERK inhibitor ulixertinib showed highly synergistic effects in vitro. These data highlight the MAPK pathway as a therapeutic target in DIPG, and show the importance of parallel resistance modelling and rational combinatorial treatments likely to be required for meaningful clinical translation.

  Study EGAS00001004495

• The Medical Genome Reference Bank: a whole genome data resource of 4,000 healthy elderly individuals.

  Allele frequency data from human reference populations is of increasing value for filtering and assignment of pathogenicity to genetic variants. Aged and healthy populations are more likely to be selectively depleted of pathogenic alleles, and therefore particularly suitable as a reference populations for the major diseases of clinical and public health importance. However, reference studies of the healthy elderly have remained under represented in human genetics. We have developed the Medical Genome Reference Bank (MGRB), a large scale comprehensive whole genome dataset of confirmed healthy elderly individuals, to provide a publicly accessible resource for health related research, and for clinical genetics. It also represents a useful resource for studying the genetics of healthy aging. The MGRB comprises 4,000 healthy, older individuals with no reported history of cancer, cardiovascular disease or dementia, recruited from two Australian community based cohorts. DNA derived from blood samples will be subject to whole genome sequencing. The MGRB will measure genome wide genetic variation in 4,000 individuals, mostly of European decent, aged 60 to 95 years (mean age 75 years). The MGRB has committed to a policy of data sharing, employing a hierarchical data management system to maintain participant privacy and confidentiality, whilst maximizing research and clinical usage of the database. The MGRB will represent a dataset of international significance, broadly accessible to the clinical and genetic research community.

  Study EGAS00001003511

• Efficacy of JAK/STAT pathway inhibition in murine xenograft models of early T-cell precursor (ETP) acute lymphoblastic leukemia

  Early T-cell precursor (ETP) acute lymphoblastic leukemia (ALL) is a recently described subtype of T-ALL characterized by a unique immunophenotype and genomic profile as well as a high rate of induction failure. Frequent mutations in cytokine receptor and JAK/STAT signaling pathways led us to hypothesize that ETP-ALL is dependent on JAK/STAT signaling. Here we demonstrate aberrant activation of the JAK/STAT pathway in ETP-ALL blasts relative to non-ETP T-ALL. Moreover, ETP-ALL showed hyperactivation of STAT5 in response to IL7, an effect that was abrogated by the JAK1/2 inhibitor ruxolitinib. In vivo, ruxolitinib displayed activity in 6/6 patient-derived murine xenograft models of ETP-ALL, with profound single-agent efficacy in 5 models. Ruxolitinib treatment decreased peripheral blast counts relative to pre-treatment levels and compared to control (P<0.01) in 5/6 ETP-ALL xenografts, with marked reduction in mean splenic blast counts (P<0.01) in 6/6 samples. Surprisingly, both JAK/STAT pathway activation and ruxolitinib efficacy were independent of the presence of JAK/STAT pathway mutations, raising the possibility that the therapeutic potential of ruxolitinib in ETP-ALL extends beyond those cases with JAK mutations. These findings establish the preclinical in vivo efficacy of ruxolitinib in ETP-ALL, a biologically distinct subtype for which novel therapies are needed.

  Study EGAS00001001146

• The clonal and mutational evolution spectrum of primary triple negative breast cancers

  Primary triple negative breast cancers (TNBC) represent approximately 16% of all breast cancers and are a tumour type defined by exclusion,for which comprehensive landscapes of somatic mutation have not been determined. Here we show in 104 early TNBC cases, that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some exhibiting only a handful of somatic aberrations in a few pathways, whereas others contain hundreds of somatic events and multiple pathways implicated. Integration with matched whole transcriptome sequence data revealed that only 36% of mutations are expressed. By examining single nucleotide variant (SNV) allelic abundance derived from deep re-sequencing (median > 20,000 fold) measurements in 2414 somatic mutations, we determine for the first time in an epithelial tumour, the relative abundance of clonal genotypes among cases in the population. We show that TNBC vary widely and continuously in their clonal frequencies at the time of diagnosis, with basal subtype TNBC exhibiting more variation than non-basal TNBC. Although p53 and somatic mutations appear clonally dominant compared with other pathways, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal and cell shape/motility proteins occurred at lower clonal frequencies, suggesting they occurred later during tumour progression. Taken together our results show that future attempts to dissect the biology and therapeutic responses of TNBC will require the determination of individual tumour clonal genotypes.

  Study EGAS00001000132

• 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

• 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

• 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

• H3.3G34R/V-transformed interneuron progenitors co-opt PDGFRA for gliomagenesis

  Histone H3.3 glycine 34 to arginine/valine (H3.3G34R/V) mutations occur in deadly hemispheric high-grade gliomas. These tumors show exquisite regional and temporal specificity, suggesting a developmental context permissive to the effects of G34R/V mutations. Here we present the molecular landscape of G34R/V gliomas (n=85) and show that 50% bear activating mutations in PDGFRA, with strong selection pressure for PDGFRAMUT clones at recurrence. We show that G34R/V tumors arise in interneuron progenitors of the foetal ventral forebrain expressing GSX2 and the DLX family of homeobox transcription factors, where terminal neuronal differentiation is impaired through aberrant G34R/V-mediated H3K27me3. Frequent co-occurrence of G34R/V & PDGFRAMUT is facilitated in this interneuron lineage-of-origin as PDGFRA forms an aberrant chromatin loop with the adjacent GSX2, hijacking its active chromatin conformation. At the single-cell level, G34R/V tumours entirely lack oligodendroglial transcriptional programs prominent in other glioma entities, and instead harbour dual neuronal and astroglial compartments. CRISPR-removal of H3.3G34R/V does not impact tumorigenicity suggesting this mutation becomes dispensable, while PDGFRAMUT are potently oncogenic regardless of G34 mutation. Collectively, our results suggest that G34R/V gliomas arise in foetal interneuron progenitors unable to undergo terminal differentiation, enabling co-option of PDGFRA through inappropriate expression and activating mutations to promote gliogenesis and oncogenicity. Reliance on PDGFRA for oncogenesis may be of therapeutic opportunity in G34R/V glioma.

  Study EGAS00001004301

• 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

• 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

• 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

• Genetic variability in exon 1 of the glucocorticoid receptor gene NR3C1 is associated with postoperative complications

  Patients having major surgery experience inflammation post-operatively, which may contribute to post-operative morbidity. Endogenous glucocorticoids (GCs) are an essential part of the stress response, but this response varies between individuals, which may in turn affect clinical outcome in general, and post-operative inflammation, specifically. Exon 1 of the NR3C1 gene, encoding the glucocorticoid receptor, contains an established region of differential regulation. DNA methylation patterns in this region have been found to differ between individuals. We investigated the methylation status and genotype in the Cytosine-phosphate-Guanine (CpG) island in exon 1 of NR3C1 in 24 patients who underwent major abdominal surgery and correlated it with post-operative complications. DNA was extracted from peripheral blood leukocytes and underwent targeted bisulfite sequencing of the CpG island. Complications were graded according to the Clavien-Dindo classification and 14 out of 24 had post-operative complications. Multifactorial and partial least square methods were used to analyse the data. A homogenous demethylated pattern was observed in all patients and no single CpG methylation was associated with post-operative complications. Four single nucleotide polymorphisms (SNPs) were significantly correlated with higher Clavien-Dindo scores. Genetic variability in the chr5:143,402,505-143,405,805 region of exon 1 of the glucocorticoid receptor gene NR3C1, but not DNA methylation, was associated with more severe post-operative complications in patients having major abdominal surgery. This indicates that the response to GCs may be of clinical importance for inflammatory conditions.

  Study EGAS00001005737

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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

• 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. & 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

• 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

• 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

• 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

• 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

• 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&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

• 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