Search Results - EGA European Genome-Phenome Archive

Center for Common Disease Genomics [CCDG] - Cardiovascular: Genetics of Coronary Heart Disease - Characterizaton of Coronary Prone Pedigrees

The primary cohort includes samples from 499 men and 154 women who survived a myocardial infarction, percutaneous transluminal angioplasty, or coronary artery bypass grafting (CABG) by age 55 years for men or age 65 years for women. Patients were seen at least six months after their acute event. Samples from another 488 men and 541 women that showed a nearly identical age, gender, and measured blood pressure distribution were collected as population-based controls. These controls tended to live in urban Salt Lake City, where the vast majority of the population is Caucasian. Blood samples were collected in the morning after 12h to 16h of fasting and prepared according to guidelines of the Lipid Research Clinic's program Manual of Laboratory Operations. All samples have serum, DNA and/or white blood cells, and a host of measurements (e.g. standard clinical chemistries, LDL and HDL cholesterol, triglycerides, glucose, HbA1C, insulin, blood pressures, anthropometric data, information on medical conditions and drug use, etc.). Additional details about the cohort or sample preparation were published previously. Here we present in total 1,475 patients that were successfully sequenced. This number includes 228 women, 419 men, and 7 of unreported gender survivors as well as 339 women and 482 men that served as population-based controls.

Study phs001901

Whole exome and transcriptome analysis of UV-exposed epidermis and carcinoma in situ reveals early drivers of carcinogenesis

Squamous cell carcinoma in situ (SCCIS) is a prevalent precancerous lesion that can progress to cutaneous squamous cell carcinoma (cSCC). Although SCCIS is common, its pathogenesis remains poorly understood. To better understand SCCIS development, we performed laser-captured microdissection of human SCCIS and adjacent epidermis to isolate genomic DNA and RNA for next generation sequencing. Whole exome sequencing (WES) identified UV-signature mutations in multiple genes including Notch 1-3 (GeneID:4851, GeneID:4853, GeneID:4854) in the epidermis and SCCIS and oncogenic TP53 (GeneID:7157) mutations in SCCIS. Novel gene families, including SCHLAFEN genes, contained UV/oxidative-signature disruptive (UVD) epidermal mutations that manifested positive selection in SCCIS. The frequency and distribution of NOTCH and TP53 mutations indicate that NOTCH mutations may precede TP53 mutations. RNA sequencing identified 1166 differentially expressed genes; the top 5 enriched GO biological processes included: 1) immune response, 2) epidermal development, 3) protein phosphorylation, 4) regulation of catalytic activity, and 5) cytoskeletal regulation. The NEURL1 ubiquitin ligase, which targets Notch ligands for degradation, was upregulated in SCCIS. Neuralized 1 protein was found to be elevated in SCCIS suggesting that increased levels could represent a novel mechanism for downregulating Notch during UV-induced carcinogenesis. The data from DNA and RNA sequencing of epidermis and SCCIS provide insights regarding SCCIS formation.

Study phs002019

Human CD4 Memory T Cell Activation Time Course

We performed a high-resolution time course RNA-seq experiment on CD4+ memory T cells during activation with anti-CD3/CD28 beads. We performed this in 24 adult individuals of European ancestry with no autoimmune diseases. We used 1 million cells as starting material, per time point (up to 8 time points per individual from 0 to 72 hours). We identified widespread dynamic allele-specific expression (dynASE) across the genome, i.e. ASE that changes with time after activation, which reflects cell state-dependent genetic regulatory effects. DynASE events were enriched in risk loci for multiple autoimmune diseases. This highlights the importance of studying multiple cellular states in order to elucidate the mechanisms of autoimmune disease. It also provides a successful framework to do it in other cellular systems relevant for other traits. We also showed an unprecedented amount of regulatory variation for HLA genes. HLA-DQB1 has 3 distinct regulatory haplotypes with expression profiles that could come from 3 different genes. We performed ATAC-seq in replicate for CD4+ memory T cells stimulated for 72 hours from another healthy individual. These data, together with genetic fine-mapping, CRISPR editing, EMSA and luciferase assays, helped us find a SNP in the promoter of HLA-DQB1 that causes at least part of the large cell state-dependent change in expression, at the RNA and protein levels.

Study phs002259

Single-cell RNA-seq of rheumatoid arthritis synovial tissue using low-cost microfluidic instrumentation

Droplet-based single-cell RNA-seq has emerged as a powerful technique for massively parallel cellular profiling. While these approaches offer the exciting promise to deconvolute cellular heterogeneity in diseased tissues, the lack of cost-effective, reliable, and user-friendly instrumentation has hindered widespread adoption of droplet microfluidic techniques. To address this, we have developed a microfluidic control instrument that can be easily assembled from 3D printed parts and commercially available components costing approximately $575. We adapted this instrument for massively parallel scRNA-seq and deployed it in a clinical environment to perform single-cell transcriptome profiling of disaggregated synovial tissue from 5 rheumatoid arthritis patients. We sequenced 20,387 single cells from synovectomies, revealing 13 transcriptomically distinct clusters. These encompass a comprehensive and unbiased characterization of the autoimmune infiltrate, including inflammatory T and NK subsets that contribute to disease biology. Additionally, we identified fibroblast subpopulations that are demarcated via THY1 (CD90) and CD55 expression. Further experiments confirm that these represent synovial fibroblasts residing within the synovial intimal lining and subintimal lining, respectively, each under the influence of differing microenvironments. We envision that this instrument will have broad utility in basic and clinical settings, enabling low-cost and routine application of microfluidic techniques, and in particular single-cell transcriptome profiling. Reprinted from [Stephenson et al., Nature Communications, 2018], with permission from the Nature Publishing Group.

Study phs001529

OncoArray: Follow-up of Ovarian Cancer Genetic Association and Interaction Studies (FOCI)

The Follow-up of Ovarian Cancer Genetic Association and Interaction Studies (FOCI) was one of five projects funded in 2010 as part of the NCI's Genetic Associations and Mechanisms in Oncology (GAME-ON) initiative (http://epi.grants.cancer.gov/gameon/). FOCI represents a collective effort that builds upon the strengths and history of collaboration inherent in the Ovarian Cancer Association Consortium (OCAC), a multidisciplinary group comprised of epidemiologists, genetic epidemiologists, statistical geneticists, molecular and cell biologists and clinicians that was formed in 2005. The other four funded GAME-ON projects were: the ColoRectal TransdisciplinaryStudy (CORECT), Elucidating Loci Involved in Prostate Cancer Susceptibility (ELLIPSE), Discovery, Biology, and Risk of Inherited Variants in Breast Cancer (DRIVE), and Transdisciplinary Research in Cancer of the Lung (TRICL). As part of our aim to discover, expand, and replicate ovarian cancer susceptibility loci, the GAME-ON projects and other consortia formed the OncoArray network (http://epi.grants.cancer.gov/oncoarray/) to develop and genotype a new custom genotyping array in large numbers of cancer cases and controls (over 400,000 samples) across multiple cancer types. The FOCI data includes over 50,000 ovarian cancer cases and controls genotyped with the Oncoarray at the Center for Inherited Disease Research (CIDR). Genotype calling and quality control procedures were performed under a standardized protocol across the Oncoarray consortium, and over 490,000 SNPs passed QC and are included under this dbGaP submission.

Study phs001882

Identification of Cancer Predisposition Genes in Breast Cancer Families

Breast cancer has a significant genetic component. The strongest epidemiologic risk factor is family history of the disease. Familial studies have identified a number of genes contributing to increased risk. Of these, the highly penetrant BRCA1 and BRCA2 genes account for the majority of families in whom disease-causal mutations have been identified. However, mutations in other genes have been associated with predisposition to breast or breast and ovarian cancer including (CHEK2, ATM, PALB2, RAD51C, RAD51D). Together these genes account for less than 20% of excess familial breast cancer risk and do not explain ~50% of high-risk families referred for genetic testing. This significant gap in our understanding of the genetic basis of breast cancer poses a significant barrier to our ability to identify women at increased risk, who would benefit greatly from appropriate clinical management. To identify additional breast cancer predisposition genes whole exome sequencing of familial high-risk breast cancer patients was undertaken. DNA samples from 50 pairs of breast cancer patients from high-risk breast cancer families with no pathogenic mutations in the BRCA1 or BRCA2 breast cancer predisposition genes were subjected to whole exome sequencing. The goal was to select candidate predisposition genes for further validation studies based on sharing of putative protein inactivating mutations among both members of each family and/or the frequency of inactivating mutations among the 50 families.

Study phs000480

Uveal Melanoma Immunogenomics Predict Immunotherapy Resistance and Susceptibility

Immunotherapy using immune checkpoint inhibition (ICI) has shown success in treating metastatic cutaneous melanoma but has limited efficacy against metastatic uveal melanoma (UM), a rare variant arising from the immune privileged eye. To better understand this immune resistance, we performed comprehensive immunogenomic profiling of 100 human UM metastases using clinicogenomics, bulk and single cell transcriptomics, T cell receptor (TCR) repertoire analysis, and tumor infiltrating lymphocyte (TIL) potency assessment. We found that over half of these metastases harbored TIL with potent autologous tumor specificity, despite having low tumor mutational burden and resistance to prior immunotherapies, including ICI and the bispecific T cell engager tebentafusp. These T cell infiltrated metastases displayed activated antigen presenting cells, chronic interferon signaling, and diverse TCR repertoires. However, we observed strikingly low intratumoral TCR clonality and transcriptionally non-proliferative TIL within the tumor microenvironment even after ICI and tebentafusp therapy, demonstrating that these immunotherapies were insufficient to induce proliferation of the tumor reactive TIL. To harness the therapeutic potential of these quiescent TIL, we developed rapid tumor transcriptomic profiling to enable their selective in vivo identification and ex vivo liberation to counter their growth suppression. We demonstrate that adoptive transfer of these transcriptomic selected TIL can promote tumor immunity in patients with metastatic UM when other immunotherapies are incapable.Data that will be available through dbGaP includes raw bulk tumor transcriptomic sequencing data.

Study phs003330

NHLBI TOPMed: Best ADd-on Therapy Giving Effective Response (BADGER)

BADGER is a 56-week randomized, double-blind, three-treatment, three-period cross-over trial that will evaluate the differential improvement in control that is achieved following three separate treatment interventions in children whose asthma is not acceptably controlled on a low dose of ICS (per NAEPP guidelines). All participants will enter an 8-week run-in period during which time they will receive a dose of 1x ICS (fluticasone 200 μg/day). During this 8-week time period, running 2-week averages to establish the lack of acceptable asthma control will be calculated. Thus, a child could qualify for randomization at any time during this 8-week run-in period. This approach should maximize both patient safety and successful enrollment. Children will continue to receive 1x ICS during the entire treatment phase. During each period of the treatment phase, they also will receive one add-on therapy in the form of LABA, LTRA or additional 1x ICS. The order of the add-on therapy assignment will be determined by randomization into one of six treatment sequences (order determined randomly). Each treatment period will be 16 weeks in length; the initial 4 weeks of each period will be considered to be the washout period for the previous treatment. The primary outcome measures will be frequency of asthma exacerbations, asthma control days, and FEV1.

Study phs001728

Developing Allelic Imbalance Analysis from Single-Nucleus RNA-Seq Data

Single-cell RNA-seq (scRNA-seq) is emerging as a powerful tool for understanding gene function across diverse cells. Recently, this has included the use of allele-specific expression (ASE) analysis to better understand how variation in the human genome affects RNA expression at the single-cell level. We reasoned that, because intronic reads are more prevalent in single-nucleus RNA-Seq (snRNA-Seq) and introns are under lower purifying selection, and are thus enriched for genetic variants, that snRNA-seq should facilitate single-cell analysis of ASE. Here we demonstrate how experimental and computational choices can improve the results of allelic imbalance analysis. We explore how experimental choices, such as RNA source, read length, sequencing depth, genotyping, etc., impact the power of ASE-based methods. We developed a new suite of computational tools to process and analyze scRNA-seq and snRNA-seq for ASE. As hypothesized, we extracted more ASE information from reads in intronic regions than those in exonic regions and show how read length can be set to increase power. Additionally, hybrid selection improved our power to detect allelic imbalance in genes of interest. We also explored methods to recover allele-specific isoform expression levels from both long- and short-read snRNA-seq. Overall, we provide an end-to-end experimental and computational approach for future studies.

Study phs003749

National Institute on Aging Genetics of Alzheimer's Disease Data Storage Site (NIAGADS)

NIAGADS is an NIH controlled-access data repository that collaborates with researchers around the world to facilitate data access and sharing for genomic and associated omics data related to Alzheimer's Disease (AD), Alzheimer's Disease Related Dementias (ADRD), as well as healthy aging. Serving as the data coordinating center for the Alzheimer's Disease Sequencing Project (ADSP), NIAGADS hosts all ADSP-generated whole-exome, whole-genome and associated phenotypes. Other data types hosted by NIAGADS include genotype SNP arrays, functional omics, fluid biomarkers, aggregated imaging data (such as volumetric data for brain regions of interest), and summary statistics. In adherence with NIH security best practices and genome data sharing policies, NIAGADS supports data access requests through review by an NIH-appointed and registered Data Access Committee. Data stored within NIAGADS cannot be requested via dbGaP. Data requests must go through the NIAGADS Data Sharing Service (DSS) via the following NIAGADS url: https://dss.niagads.org/. Questions related to the request of data should be directed to the NIAGADS admin team at: niagads@pennmedicine.upenn.edu. DSS is a FISMA certified platform built by a team of data scientists at the University of Pennsylvania to support the request and access of AD/ADRD data by qualified researchers internationally. Questions related to the use of data hosted in the DSS should be directed to the Portal admin team at: help@niagads.org.

Study phs004180

The Haplotype-Resolved Genome and Epigenome of the Aneuploid HeLa Cancer Cell Line

The genomic architecture of HeLa remains largely unexplored beyond its karyotype, partly because like many cancers, its extensive aneuploidy renders such analyses challenging. We carried out haplotype-resolved whole-genome sequencing of the HeLa CCL-2 strain, examined point- and indel-mutation variations, mapped copy-number variations and loss of heterozygosity regions, and phased variants across full chromosome arms. We also investigated variation and copy-number profiles for HeLa S3 and eight additional strains. We find that HeLa is relatively stable in terms of point variation, with few new mutations accumulating after early passaging. Haplotype resolution facilitated reconstruction of an amplified, highly rearranged region of chromosome 8q24.21 at which integration of the human papilloma virus type 18 (HPV-18) genome occurred and that is likely to be the event that initiated tumorigenesis. We combined these maps with RNA-seq and ENCODE Project data sets to phase the HeLa epigenome. This revealed strong, haplotype specific activation of the proto-oncogene MYC by the integrated HPV-18 genome approximately 500 kilobases upstream, and enabled global analyses of the relationship between gene dosage and expression. These data provide an extensively phased, high-quality reference genome for past and future experiments relying on HeLa, and demonstrate the value of haplotype resolution for characterizing cancer genomes and epigenomes. Reprinted from Adey et. al. 2013, Nature, with permission from Nature Publishing Group.

Study phs000642

Chromatin remodeling enhances MAP3K8 expression in HAM: a key pathogenesis for therapeutic intervention

Human T-cell leukemia virus type 1 (HTLV-1)-associated myelopathy (HAM) is a debilitating neuroinflammatory disease with no available effective treatments. A hallmark of HAM is the transformation of HTLV-1-infected cells into T helper type 1 (Th1)-like cells, characterized by excessive interferon (IFN)-gamma production that drives chronic inflammation. However, the molecular mechanisms fuel this aberrant Th1-like transformation and sustained inflammation remain poorly understood. We hypothesized that HAM-characteristic chromatin remodeling plays a pivotal role in the overexpression of key genes driving inflammatory pathogenesis. Using transcriptomic analysis, chromatin accessibility profiling, and biomarker evaluations across HTLV-1-related diseases, we identify MAP3K8 as a key gene that defines the unique inflammatory profile of HAM. MAP3K8 overexpression promotes Th1-like differentiation and constitutively activates the MEK-ERK signaling pathway. Furthermore, we elucidate the mechanism by which HTLV-1 Tax, Fosl2, and c-Jun collaboratively induce HAM-characteristic chromatin remodeling at the enhancer region of the MAP3K8 locus. Crucially, we demonstrate that mitogen-activated protein kinase kinase (MEK) inhibitors effectively suppress the MAP3K8-MEK signaling cascade and significantly mitigated inflammatory pathogenesis in an ex vivo culture assay. Our findings provide critical insights into the virus-host interactions underpinning HAM and propose the MAP3K8-MEK-ERK axis as a promising therapeutic target for this challenging condition.

Study JGAS000835

Inherited chromosomally-integrated human herpesvirus 6A/B (HHV-6A/B) genome sequences in the Japanese population

The deposited sequences represent inherited chromosomally-integrated human herpesvirus 6 (iciHHV-6) genomes identified from 10 BioBank Japan (BBJ) subjects based on the analysis of the whole-genome sequencing (WGS) data. In brief, we re-aligned unmapped reads from the WGS to HHV-6A genome (Genbank accession number: KJ123690.1) using BWA-MEM algorithm. Ten out of 3,256 subjects were identified with high depth of coverage of HHV-6. We further distinguished subjects with integrated HHV-6A subtypes from those with integrated HHV-6B based on phylogenetic analysis of concatenated sequences of three viral genes (U27/U43/U83, which show high divergence between viral species). We reconstructed the HHV-6 viral genome of 10 subjects. First, reads that mapped to a reference HHV-6A genome (KJ123690.1) were further aligned against the integrated HHV-6A genome derived from a Japanese individual NA18999 (GenBank Accession number: KY316047.1) using BWA-MEM. Based on the alignment, variant calling was performed using freebayes (version: v1.2.0-2-g29c4002) with parameters ploidy = 1 and min-alternate-fraction = 0.8 {Garrison, 2012}. We generated subject-specific iciHHV-6 viral sequences by applying the resulting variants to the KY316047.1 reference genome using the FastaAlternateReferenceMaker function in the Genome Analysis Toolkit (GATK) v3.7. Please note that HHV-6B-specific variants are unlikely to be sensitively reconstructed using this method.

Study JGAS000240

Treg cell subset-specific gene expression patterns in human head and neck cancer

Backgrand: The occurrence of head and neck cancer has been increased all over the world. Although the prognosis has been improved by using several treatment modalities including immunotherapy, it is still unsatisfying. It is still unclear why the immunotherapy is effective only for partial amount of patients. To resolve this problem and to make a new treatment strategy based on human immune system, further investigation of antitumor immunity is required. Regulatory T (Treg) cells play critical roles in human immune responses including antitumor immunity. Therefore, the research to elucidate Treg function is strongly urged to clarify the mechanisms of antitumor immunity. Methods: We performed RNA-seuencing analysis of regulatory T cells infiltrating head and neck cancer. We devided into four populations, DP, CD25+, CTLA4+, DN and analyzed gene expressions. Results: The gene expression patterns of Treg-related genes were different between CD25+ cells and CD25- cells. Furthermore, DP cells had unique gene expression pattern compared to other populations. Conclusions: We found Treg cells infiltrating head and neck cancer had high expression levels of Treg-related genes and DP cells had unique gene expression pattern that may become a target of new treatment strategy. Regulatory T (Treg) cells play critical roles in human immune responses including antitumor immunity. Therefore, the research to elucidate Treg function is strongly urged to clarify the mechanisms of antitumor immunity.

Study JGAS000135

RBMX functional retrocopy safeguards brain development

Retrotransposition has generated thousands of intronless gene copies in mammalian genomes, yet their contribution to brain development and evolution remains largely unexplored. We investigated the role for RBMX retrocopy in shaping neurodevelopment and modulating disease. First, we reported 8 distinct pathogenic variants in RBMX, an X-linked splicing regulator, in 11 males presenting with intellectual disability, microcephaly, and cortical malformations. Then, we combined transcriptomic profiling, protein–protein and protein–RNA interaction studies both in human cellular models and mouse embryonic cortices to assess the functional redundancy between RBMX and its retrocopy RBMXL1. Finally, we use mouse genetics to dissect RBMX function and its compensation by RBMXL1 in corticogenesis. We show that RBMX pathogenic variants disrupt cortical development through both loss- and gain-of-function mechanisms. Surprisingly, despite severe phenotypes in humans, Rbmx-deficient mice display only mild cortical abnormalities, a discrepancy likely due to compensation by Rbmxl1, a retrocopy that arose independently in mice and humans. We demonstrate that RBMX and RBMXL1 share protein and RNA partners and act redundantly in brain development, with RBMXL1 buffering the impact of RBMX deficiency. Together, these findings establish that RBMXL1 functions as a functional paralog of RBMX safeguarding neurodevelopment by buffering deleterious variation. More broadly, these results identify retrocopies as active contributors to neurodevelopmental robustness and suggest that functional retrocopies may have facilitated the evolutionary diversification of the mammalian brain.

Study EGAS50000001650

Whole genome sequencing of matched primary and metastatic acral melanomas

Next generation sequencing has enabled systematic discovery of mutational spectra in cancer samples. Here we used whole genome sequencing to characterise somatic mutations and structural variation in a primary acral melanoma and its lymph node metastasis. Our data show that the somatic mutational rates in this acral melanoma sample pair was more comparable to the rate reported in cancer genomes not associated with mutagenic exposure than in the genome of a melanoma cell line or the transcriptome of melanoma short-term cultures. Despite the perception that acral skin is sun-protected, the dominant mutational signature in these samples is compatible with damage due to ultraviolet light exposure. A nonsense mutation in ERCC5 discovered in both the primary and metastatic tumours, could also have contributed to the mutational signature through accumulation of unrepaired dipyrimidine lesions. However, evidence of transcription-coupled repair was suggested by the lower mutational rate in the transcribed regions and expressed genes. The primary and the metastasis are highly similar at the level of global gene copy number alterations, loss of heterozygosity and single nucleotide variation (SNV). Furthermore, the majority of the SNVs in the primary tumour were propagated in the metastasis and one non-synonymous coding SNV and one splice site mutation appeared to arise de novo in the metastatic lesion.

Study EGAS00001000169

A thymic ILC1-like progenitor with differentiation potential towards KIR+NKG2A- NK cells

KIR+NKG2A- Natural Killer (NK) cells have the unique ability to detect downregulation of single HLA-I allotypes, frequently occurring in malignantly-transformed and virus-infected cells in order to escape T cell recognition. We have recently shown that circulating Innate Lymphoid Cells 1 (cILC1) have the potential to generate such KIR+NKG2A- NK cells but their developmental origin was unknown. Here we demonstrate that the development of cILC1 is thymus-dependent and identify a putative progenitor of cILC1s in the thymus, called thyILC1. Single-cell RNAseq analysis revealed a close relationship of thyILC1s to CD34+ double-negative (DN) thymocytes, but in contrast to DNs, thyILC1s had lost T cell differentiation potential. However, like cILC1s they could be efficiently differentiated into KIR+NKG2A- NK cells. thyILC1s strongly expressed RAG1/2, PTCRA (preTCR, and NOTCH1/3, but lacked expression of TCR chains, suggesting that they failed to pass the -selection checkpoint of thymocyte development. Finally, analyses of patients with FOXN1het deficiency, showing thymic hypoplasia due to lack of mature thymic epithelial cells, exhibited a profound deficiency of cILC1s but not cILC2s and cILC3s, demonstrating thymus dependency of cILC1s. Together, the data suggest that thyILC1s are the source of a novel thymus-dependent NK cell differentiation pathway that promotes generation of KIR+NKG2A- NK cells.

Study EGAS50000000760

Somatic rearrangements causing oncogenic ectodomain deletions of FGFR1 in squamous cell lung cancer

The discovery of frequent 8p11-p12 amplifications in squamous cell lung cancer (SQLC) has fueled hopes that FGFR1, located inside this amplicon, might be a therapeutic target. In a clinical trial, only 11% of patients with 8p11 amplification (detected by FISH) responded to FGFR kinase inhibitor treatment. To understand the mechanism of FGFR1 dependency, we performed deep genomic characterization of 52 SQLCs with 8p11-p12 amplification, including 10 tumors obtained from patients who had been treated with FGFR inhibitors. We discovered somatically altered variants of FGFR1 with deletion of exons 1-8 that resulted from intragenic tail-to-tail rearrangements. These ectodomain-deficient FGFR1 variants (ΔEC-FGFR1) were expressed in the affected tumors and were tumorigenic in both in vitro and in vivo models of lung cancer. Mechanistically, breakage-fusion-bridges were the source of 8p11-p12 amplification, resulting from frequent head-to-head and tail-to-tail rearrangements. Generally, tail-to-tail rearrangements within or in close proximity upstream of FGFR1 were associated with FGFR1 dependency. Thus, the genomic events shaping the architecture of the 8p11-p12 amplicon provide a mechanistic explanation for the emergence of FGFR1-driven SQLC. Specifically, we believe that FGFR1 ectodomain-deficient and FGFR1-centered amplifications caused by tail-to-tail rearrangements are a novel somatic genomic event that might be predictive of therapeutically relevant FGFR1 dependency.

Study EGAS50000001368

MRCA and MRCE SNP genotypes

Gene expression levels can be an important link DNA between variation and phenotypic manifestations. Our previous map of global gene expression based on approximately 400K SNPs and 50K transcripts in 400 sib pairs from the MRCA family panel has been widely used to interpret the results of GWAS. Here, we more than double the size of our initial dataset with expression data on 550 additional individuals from the MRCE family panel using the Illumina whole genome expression array. We have used new statistical methods for dimension reduction to account for non-genetic effects in estimates of expression levels and we have also included SNPs imputed from the 1000 Genomes Project. Our methods reduced false discovery rates and increased the number of eQTLs mapped either locally or at a distance (i.e. in cis or trans) from 1,534 in the MRCA dataset to 4,452 (with <5% FDR). Imputation of 1000 Genomes SNPs further increased the number of eQTLs to 7,302. Using the same methods and imputed SNPs in the newly acquired MRCE dataset we identified eQTLs for 9000 genes. The combined results identify strong local and distant effects for transcripts from 14,177 genes. Our eQTL database (http://www.hsph.harvard.edu/faculty/liming-liang/software/eqtl/) based on these results is freely available to help define the function of disease-associated variants.

Study EGAS00000000137

A prognostic human brain network for diffuse midline glioma

Diffuse midline glioma (DMG) are near-universally lethal tumours of the paediatric central nervous system. In animal models, DMG form brain-wide, integrated networks through neuron-to-glioma synapses and glioma-to-glioma gap junctional coupling. This extensive connectivity robustly promotes the growth and invasion of DMG and other glial malignancies through paracrine mechanisms and direct neuron-to-glioma synapses. The organisation and clinical implications of these connections in the living human brain, however, remain to be elucidated. Here, we develop tumour network mapping to compute the brain-wide connectivity profile of DMG: defining a conserved brain network across pontine and thalamic DMG associated with patient short-term survival (DMG network). Tumour functional connectivity with the DMG network was independently predictive of patient overall survival across two external validation cohorts. Tumour growth mapped to DMG network-specific trajectories and peak, in-network neurometabolic changes across development spatiotemporally aligned with the peak age incidence of DMG. Single-nucleus RNA-sequencing confirmed diverse synaptic gene enrichment in high-connectivity DMG. Strikingly, incidental surgical resection of high-connectivity thalamic DMG tissue conferred a significant survival advantage. Collectively, these data define a conserved and prognostically important brain network in children with DMG, consistent with the hypothesis that DMG exploit otherwise healthy brain circuits to promote tumour growth.

Study EGAS50000001752

DUX4 activates novel intergenic transcripts and isoforms in a tissue-specific context

DUX4, a transcription factor driving zygotic genome activation during the embryonic cleavage stage in mammals, regulates the expression of genic and non-genic elements and modulates chromatin accessibility. Its misexpression in skeletal muscle causes Facioscapulohumeral muscular dystrophy (FSHD). The sporadic expression of DUX4 in FSHD muscle and the limitations of short-read RNA sequencing challenge, however, a comprehensive grasp of its consequences on the muscle transcriptome. Here, leveraging full-length isoform RNA sequencing with short-read RNA sequencing on DUX4-inducible myoblasts, we elucidated DUX4-induced transcriptomic changes within a myogenic context. DUX4 yields a complex transcriptome featuring numerous unannotated isoforms from known and novel loci. Moreover, we observed distinct isoform usage of DUX4 target genes between myogenic and embryonic cleavage states highlighting the distinct differences between DUX4 transcription factor activity in skeletal muscle and during the cleavage stage. We also identified ~1,400 unannotated intergenic isoforms predominantly originating from repetitive elements such as LTRs, LINEs, and SINEs. Integration with available datasets delineates the transcriptional and epigenetic landscape of these loci across myogenic and embryonic backgrounds. Like the cell type-specific isoform production, the usage of DUX4 binding sites at these intergenic loci is in part cell type-specific. Collectively, these findings demonstrate substantial and context-specific transcriptomic alterations induced by DUX4 further enriching our understanding of DUX4-induced muscle pathology.

Study EGAS50000000503

Phenotypic characterization and prognostic impact of CD103+ tissue-resident memory T cells in diffuse large B cell lymphoma

Tissue-resident memory T (TRM) cells are a population of memory T cells that stably occupy tissues and play a key role in immunosurveillance, having been linked to favorable survival outcomes in various solid tumors. While TRM cells have been identified in lymph nodes, their phenotype and prognostic significance in B-cell non-Hodgkin lymphoma (B-NHL), including diffuse large B-cell lymphoma (DLBCL), remains poorly characterized. The frequency of CD103 expressing T cells in patient samples with DLCBL was quantified by immunofluorescence (IF) staining of tissue biopsies (n=306) and flow cytometry (n=252) of cell disaggregates, and linked with clinical outcome. The phenotype of TRM cells was characterized by spectral flow cytometry (n=10 DLBCL and 2 reactive lymph node (rLN) samples) and single-cell RNA sequencing (scRNAseq) (n=12 aggressive B cell lymphomas and 4 rLN samples). An additional 51 samples from an external dataset were included to verify the single cell transcriptome findings. The flow cytometry and scRNAseq findings in DLBCL were extrapolated to additional B-NHL entities. Through analysis of these datasets we were able to characterize CD103+ TRM-like cells in DLBCL and other B-NHL entities and identified them to represent a prognostically favorable population with an activated/cytotoxic T cell phenotype.

Study EGAS50000000943

Liquid biopsy for molecular characterization of diffuse large B-cell lymphoma and early assessment of minimal residual disease

Circulating tumor DNA (ctDNA) allows genotyping and minimal residual disease (MRD) detection in lymphomas. Using a NGS approach (Euroclonality-NDC), we evaluated clinical and prognostic value of ctDNA in a series of R-CHOP-treated DLBCL patients at baseline (n=68) and after 2-cycles (n=59), monitored by metabolic imaging (PET/CT). A molecular marker was identified in 61/68 (90%) ctDNA samples at diagnosis. Pre-treatment high ctDNA levels significantly correlated with elevated LDH, advanced stage, high risk IPI and a trend to shorter 2-year PFS. Valuable NGS data after 2-cycles of treatment were obtained in 44 cases, and 38 achieved major molecular response (MMR; 2.5-log drop in ctDNA). PFS curves displayed statistically significant differences among those achieving MMR vs. those not achieving MMR (2yr PFS of 76% vs. 0%, p<0.001). Similarly, more than 66% reduction in SUVmax by PET/CT identified two subgroups with different prognosis (2yr PFS of 83% vs. 38%; p<0.001). Combining both approaches MMR and SUVmax reduction, a better stratification was observed (2yr PFS of 84% vs. 17% vs. 0%, p<0.001). Euroclonality-NDC panel allows the detection of a molecular marker in the ctDNA in 90% of DLBCL. ctDNA reduction at 2 cycles and its combination with interim PET results improves patient prognosis stratification.

Study EGAS50000000215

The neo-open reading frame peptides that comprise the tumor framome are a rich source of neoantigens for cancer immunotherapy

Identification of immunogenic cancer neoantigens as targets for therapy is challenging. Here, we integrate whole genome and long-read transcript sequencing of cancers to identify the collection of neo-open reading frame peptides (NOPs) expressed in tumors. We termed this collection of NOPs the tumor framome. NOPs represent tumor-specific peptides that are different from wild-type proteins and may be strongly immunogenic. We describe a class of hidden NOPs that derive from structural genomic variants involving an upstream protein coding gene driving expression and translation of non-coding regions of the genome downstream of a rearrangement breakpoint, i.e. where no gene annotation or evidence for transcription exists. The entire collection of NOPs represents a vast number of possible neoantigens particularly in tumors with many structural genomic variants and a low number of missense mutations. We show that NOPs are immunogenic and epitopes derived from NOPs can bind to MHC class I molecules. Finally, we provide evidence for the presence of memory T cells specific for hidden NOPs in peripheral blood from a patient with lung cancer. This work highlights NOPs as a major source of possible neoantigens for personalized cancer immunotherapy and provides a rationale for analyzing the complete cancer genome and transcriptome as a basis for detection of NOPs.

Study EGAS00001006021

Personalized peptide vaccination among 173 patients with IDH wildtype glioblastoma: a retrospective analysis

Tumors from 173 GBM patients were analysed for somatic mutations to generate a personalized peptide vaccine targeting tumor-specific neoantigens. Exome libraries for 173 glioblastoma tumors and matched normal DNA were sequenced on Illumina platform, alongside whole transcriptomes from the tumor samples. Paper Abstract: Current treatment outcome of patients with glioblastoma (GBM) remains poor. Following standard therapy, recurrence is universal with limited survival. GBM tumors from 173 patients were analysed for somatic mutations to generate a personalized peptide vaccine targeting tumor-specific neoantigens. Among all patients, including 70 treated prior to progression (primary) and 103 treated after progression (recurrent), the median overall survival from first diagnosis was 31.9 months (95% CI: 25.0-36.5). Side effects were infrequent and were predominantly grade 1 or 2. A vaccine-induced immune response to at least one of the vaccinated peptides was detected in blood samples of 87 of 99 (88%) monitored patients. T-cell responses to vaccinated neoepitope peptides were durable in most patients. Significantly prolonged survival was observed for patients with multiple vaccine-induced immune responses (53 months) compared to those with no/low induced responses (27 months; P=0.03). Altogether, our results highlight that the application of personalized neoantigen-targeting peptide vaccine is feasible and represents a promising potential treatment option for GBM patients.

Study EGAS50000000449

Organoid BulkRNAseq

Here, we aim to understand early determinants of brain organoid quality and generated a set of brain organoids from 12 different hPSC lines using an adaptation of the unguided Lancaster protocol (Lancaster & Knoblich, 2014). All cell lines were obtained from healthy donors. 12 donor hPSC lines were differentiated into brain organoids and a set of 72 individual organoids (6 per line) were randomly selected morphologically diverse organoids and subjected to bulk RNA sequencing for gene expression analysis. RNA was extracted using the RNeasy kit (Qiagen, Germany) according to the manufacturer’s instructions. A total of 1000 ng per sample was sent for RNA sequencing to Azenta Life Sciences (Genewiz Leipzig, Germany) for sequencing. Library preparation was performed using the NEBNext Ultra II RNA Library Prep Kit for Illumina following manufacturer’s instructions (NEB, Ipswich, MA, USA). mRNAs were first enriched with Oligo(dT) beads. The samples were sequenced using a 2x150 Pair-End (PE) configuration v1.5 at a depth of >20 million reads in each sample. Raw sequence data (.bcl files) generated from Illumina NovaSeq6000 was converted into fastq files and de-multiplexed using Illumina bcl2fastq program version 2.20. The raw sequence reads were aligned to the human reference genome GRCh38 using the STAR aligner (Dobin et al., 2013). Gene expression levels were quantified using the FeatureCounts tool (Liao et al., 2014).

Study EGAS50000000659

A thymic ILC1-like progenitor with differentiation potential towards KIR+NKG2A- NK cells - sc

KIR+NKG2A- Natural Killer (NK) cells have the unique ability to detect downregulation of single HLA-I allotypes, frequently occurring in malignantly-transformed and virus-infected cells in order to escape T cell recognition. We have recently shown that circulating Innate Lymphoid Cells 1 (cILC1) have the potential to generate such KIR+NKG2A- NK cells but their developmental origin was unknown. Here we demonstrate that the development of cILC1 is thymus-dependent and identify a putative progenitor of cILC1s in the thymus, called thyILC1. Single-cell RNAseq analysis revealed a close relationship of thyILC1s to CD34+ double-negative (DN) thymocytes, but in contrast to DNs, thyILC1s had lost T cell differentiation potential. However, like cILC1s they could be efficiently differentiated into KIR+NKG2A- NK cells. thyILC1s strongly expressed RAG1/2, PTCRA (preTCR, and NOTCH1/3, but lacked expression of TCR chains, suggesting that they failed to pass the -selection checkpoint of thymocyte development. Finally, analyses of patients with FOXN1het deficiency, showing thymic hypoplasia due to lack of mature thymic epithelial cells, exhibited a profound deficiency of cILC1s but not cILC2s and cILC3s, demonstrating thymus dependency of cILC1s. Together, the data suggest that thyILC1s are the source of a novel thymus-dependent NK cell differentiation pathway that promotes generation of KIR+NKG2A- NK cells.

Study EGAS50000000790

Raw sequencing files from WGS and RNA-Seq

Acute leukemias are the most common malignancies in children with an onset of the disease among 0-6 years-old. The early development of theses malignancies has been proposed as an indicator of their prenatal origin, hypothesis that has been confirmed especially in the case of acute lymphoblastic leukemia. For acute myeloid leukemia (AML), comprising an heterogenous group of different cytogenetics subtypes, this prenatal origin has also been confirmed specially for MLL-rearranged and t(8;21) subgroups. AMLs carrying t(7;12) constitute up to 30% of cases under 2-years-old, frequently being associated with a trisomy of chromosome 19 and an overexpression of MNX1 gene. An in utero origin of this specific subgroup has been suggested but not clearly confirmed. Here, we present the case of a 5-months-old AML patient with t(7;12) in which we have fully characterized a novel breakpoint occurring among ETV6 and NOM1 genes. We demonstrate the presence of this specific ETV6-NOM1 fusion in CD34+ hematopoietic stem and progenitor cells isolated from the cord blood unit from the same patient and stored at birth. Moreover, trisomy +19 present at diagnosis was not detected at birth. Together, our results demonstrate the prenatal origin of the t(7;12) in infant AML, and establish trisomy +19 as a secondary event in this leukemia.

Study EGAS50000000186

A Cancer Cell-Line Titration Series for Evaluating Somatic Classification

Accurate detection of somatic single nucleotide variants and small insertions and deletions from DNA sequencing experiments of tumour-normal pairs is a challenging task. Tumour samples are often contaminated with normal cells confounding the available evidence for the somatic variants. Furthermore, tumours are heterogeneous so sub-clonal variants are observed at reduced allele frequencies. We present here a cell-line titration series dataset that can be used to evaluate somatic variant calling pipelines with the goal of reliably calling true somatic mutations at low allele frequencies. Cell-line DNA was mixed with matched normal DNA at 8 different ratios to generate samples with known tumour cellularities, and exome sequenced on Illumina HiSeq to depths of >300x. The data was processed with several different variant calling pipelines and verification experiments were performed to assay >1,500 somatic variant candidates using Ion Torrent PGM as an orthogonal technology. By examining the variants called at varying cellularities and depths of coverage, we show that the best performing pipelines are able to maintain a high level of precision at any cellularity. In addition, we estimate the number of true somatic variants undetected as cellularity and coverage decrease. Our cell-line titration series dataset, along with the associated verification results, is effective for this evaluation and will serve as a valuable dataset for future somatic calling algorithm development.

Study EGAS00001001016

Transposable Elements Shape Stemness in Normal and Leukemic Hematopoiesis

Despite most acute myeloid leukemia (AML) patients achieving complete remission after induction chemotherapy, two-thirds will relapse with fatal disease within five years. AML is organized as a cellular hierarchy sustained by leukemia stem cells (LSC) at the apex, with LSC properties directly linked to tumor progression, therapy failure, and disease relapse 1–5. Despite the central role of LSC in poor patient outcomes, little is known about the genetic determinants driving their stemness properties. As LSCs share many functional and molecular properties with normal hematopoietic stem cells (HSC) 6, we investigated accessible chromatin unique across normal hematopoietic and cancer cell states and identified transposable elements (TEs) as genetic determinants of both primitive populations in comparison with their downstream mature progeny. A clinically-relevant TE chromatin accessibility-based LSCTE121 signature was developed that enabled patient classification based on survival outcomes. Through functional assays, primitive cell specific-TE subfamilies were found to serve as docking sites for stem cell-associated regulators of genome topology or lineage-specific transcription factors, including LYL1 in LSCs. Finally, using chromatin editing tools, we establish that chromatin accessibility at LTR12C elements in LSCs are necessary to maintain stemness properties. Our work identifies TEs as genetic drivers of primitive versus mature cell states, where distinct TE subfamilies account for stemness properties in normal versus leukemic hematopoietic stem cells.

Study EGAS50000001141

Transposable Elements Shape Stemness in Normal and Leukemic Hematopoiesis

Despite most acute myeloid leukemia (AML) patients achieving complete remission after induction chemotherapy, two-thirds will relapse with fatal disease within five years. AML is organized as a cellular hierarchy sustained by leukemia stem cells (LSC) at the apex, with LSC properties directly linked to tumor progression, therapy failure, and disease relapse 1–5. Despite the central role of LSC in poor patient outcomes, little is known about the genetic determinants driving their stemness properties. As LSCs share many functional and molecular properties with normal hematopoietic stem cells (HSC) 6, we investigated accessible chromatin unique across normal hematopoietic and cancer cell states and identified transposable elements (TEs) as genetic determinants of both primitive populations in comparison with their downstream mature progeny. A clinically-relevant TE chromatin accessibility-based LSCTE121 signature was developed that enabled patient classification based on survival outcomes. Through functional assays, primitive cell specific-TE subfamilies were found to serve as docking sites for stem cell-associated regulators of genome topology or lineage-specific transcription factors, including LYL1 in LSCs. Finally, using chromatin editing tools, we establish that chromatin accessibility at LTR12C elements in LSCs are necessary to maintain stemness properties. Our work identifies TEs as genetic drivers of primitive versus mature cell states, where distinct TE subfamilies account for stemness properties in normal versus leukemic hematopoietic stem cells.

Study EGAS50000001147

Aggressive genomic features in clinicallyindolent primary HHV8-negative effusion-based lymphoma

Human herpesvirus 8-negative effusion-based lymphoma (HHV8-negative EBL) is a distinct lymphoma entity related to fluid overload states caused by underlying medical conditions. Due to a phenotypical resemblance with secondary effusion of diffuse large B-cell lymphoma (DLBCL), HHV8-negative EBL is easily misdiagnosed. However, in contrast to DLBCL, HHV8-negative EBL follows a mild clinical course, underpinning the importance to differentiate between the two entities. Only limited molecular data of HHV8-negative EBL has been reported and characteristic genomic patterns are lacking. We performed extensive genomic profiling on 8 HHV8-negative EBL cases. A mean number of 33.6 copy number gains and losses per case were found, with recurrent focal deletions of chromosome 3p14.2 and 6q21 in 6/8 cases, encompassing tumor supressor genes FHIT and PRDM1, respectively. We observed a high rate of mutations across genes frequently involved in DLBCL, including patterns of somatic hypermutation. The most abundant hotspot mutation was found in MYD88L265P (3/8 cases). In addition, 5/8 cases had one or more translocations, involving MYC, BCL2, BCL6, TP63, EXOC2 and KMT2D, with Ig and non-Ig partners. We conclude that HHV8-negative EBL is characterized by genomic features that are commonly found in more aggressive B-cell lymphomas, but despite this follow an indolent clinical course.

Study EGAS00001002743

Analysis_of_genomic_integrity_of_disease_corrected_human_induced_pluripotent_stem_cells_by_exome_sequencing

Induced pluripotent stem (iPS) cells hold great promise for autologous cell transplantation. In order to apply this therapy to monogenic disorder, disease-causing mutations must be corrected prior to transplantation. We generated iPS cells from patients with alpha-1 antitrypsin deficiency, which is caused by homozygous non-synonymous mutation in the SERPINA1 gene. In vitro differentiation of patient-derived iPS cells to hepatocytes showed recapitulation of the disease phenotype. We corrected the mutations by gene targeting using zinc finger nuclease (ZFN) and piggyBac transposon technologies. Our mutation correction is highly accurate. There is no exogenous sequence left behind at the targeted site. The corrected iPS cell-derived hepatocytes were indistinguishable from healthy hepatocyte, suggesting that the corrected gene was fully functional. The genomic integrity was also analyzed by CGH and SNP array, showing no genomic alterations. These are important points to ensure clinical safety as well as functional restoration of the corrected gene. One remaining question is if these corrected iPS cells are clinically safe at the single base-pair level. Recent publications suggested that cancer-related mutations are significantly accumulated during iPS cell propagation. We aim to assess the clinical safety of the corrected iPS cells by using exome sequencing. Outcome from this project would bring iPS cell technology one step closer to clinical applications.

Study EGAS00001000055

Evolution of DNA repair defects during malignant progression of low-grade gliomas after temozolomide treatment.

Temozolomide (TMZ) increases the overall survival of patients with glioblastoma (GBM), but its role in the clinical management of diffuse low-grade gliomas (LGG) is still being defined. DNA hypermethylation of the O 6 -methylguanine-DNA methyltransferase (MGMT) promoter is associated with an improved response to TMZ treatment, while inactivation of the DNA mismatch repair (MMR) pathway is associated with therapeutic resistance and TMZ-induced mutagenesis. We previously demonstrated that TMZ treatment of LGG induces driver mutations in the RB and AKT-mTOR pathways, which may drive malignant progression to secondary GBM. To better understand the mechanisms underlying TMZ-induced mutagenesis and malignant progression, we explored the evolution of MGMT methylation and genetic alterations affecting MMR genes in a cohort of 34 treatment-naïve LGGs and their recurrences. Recurrences with TMZ-associated hypermutation had increased MGMT methylation compared to their untreated initial tumors and higher overall MGMT methylation compared to TMZ-treated non-hypermutated recurrences. A TMZ-associated mutation in one or more MMR genes was observed in five out of six TMZ-treated hypermutated recurrences. In two cases, pre-existing heterozygous deletions encompassing MGMT, or an MMR gene, were followed by TMZ-associated mutations in one of the genes of interest. These results suggest that tumor cells with methylated MGMT may undergo positive selection during TMZ treatment in the context of MMR deficiency.

Study EGAS00001001179

Exome and RNA sequencing of relapsed TCF3-PBX1 t(1;19) acute lymphoblastic leukemia

TCF3-PBX1 (E2A-PBX1) is a recurrent gene fusion in B-cell precursor lymphoblastic leukemia (BCP-ALL), which is caused by the translocation t(1;19)(q23;p13). TCF3-PBX1 BCP-ALL patients typically benefit from chemotherapy; however, many relapse and subsequently develop resistant disease with few effective treatment options. Mechanisms driving disease progression and therapy resistance have not been studied in TCF3-PBX1 BCP-ALL. Here, we aimed to identify novel treatment options for TCF3-PBX1 BCP-ALL by profiling leukemia cells from a relapsed patient, and determine molecular mechanisms underlying disease pathogenesis and progression. By drug sensitivity testing of leukemic blasts from the index patient, control samples and TCF3-PBX1 positive and negative BCP-ALL cell lines, we identified the phosphatidylinositide 3-kinase delta (p110δ) inhibitor idelalisib as an effective treatment for TCF3-PBX1 BCP-ALL. This was further supported by evidence showing TCF3-PBX1 directly regulates expression of PIK3CD, the gene encoding p110δ. Other somatic mutations to TP53 and MTOR, as well as aberrant expression of CXCR4, may influence additional drug sensitivities specific to the index patient and accompanied progression of the disease. Our results suggest that idelalisib is a promising treatment option for patients with TCF3-PBX1 BCP-ALL while other drugs could be useful depending on the genetic context of individual patients.

Study EGAS00001001876

Comparison of 133/144bp dominant phenotype vs 166bp dominant phenotype.

Background: Previous studies found that cell-free DNA (cfDNA) generated from tumors was shorter than that from healthy cells, and selecting short cfDNA could enrich for tumor cfDNA and improve its usage in early cancer diagnosis and treatment monitoring; however, the underlying mechanism of shortened tumor cfDNA was still unknown, which potentially limits its further clinical application.Results: Using targeted sequencing of cfDNA in a large cohort of solid tumor patient, sequencing reads harboring tumor-specific somatic mutations were isolated to examine the exact size distribution of tumor cfDNA. For the majority of studied cases, 166bp remained as the peak size of tumor cfDNA, with tumor cfDNA showing an increased proportion of short fragments (100-150bp). Less than 1% of cfDNA samples were found to be peaked at 134/144bp and independent of tumor cfDNA purity. Using whole-genome sequencing of cfDNA, we discovered a positive correlation between cfDNA shortening and the magnitude of chromatin inaccessibility, as measured by transcription, DNase I hypersensitivity, and histone modifications. Tumor cfDNA shortening occurred simultaneously at both 5’ and 3’ ends of the DNA wrapped around nucleosomes. Conclusions: Tumor cfDNA shortening exhibited two distinctive modes. Tumor cfDNA purity and chromatin inaccessibility were contributing factors but insufficient to trigger a global transition from 166bp dominant to 134/144bp dominant phenotype.

Study EGAS00001004135

Genomic analysis reveals novel secondary drivers and progression pathways in skin basal cell carcinoma

Basal cell carcinoma of the skin (BCC) is the most common malignant neoplasm in humans. BCC is primarily driven by the aberrant activation of the Sonic Hedgehog (Hh) pathway. However, its extensive phenotypical variation remains to be explained. The genetic profiling of 293 BCCs revealed the highest mutation rate observed in cancer (65 Mutations/Mb), with strong prevalence of UV-light signature mutations. 85% of BCCs harbored mutations in Hh pathway genes: mutually exclusive PTCH1 (73%) and SMO (20%) (P=6.6x10-8), SUFU (8%), and in TP53 (61%). 85% of BCCs also harbored secondary driver mutations in other genes implicated in BCC tumorigenesis. Recurrent secondary driver mutations were observed in MYCN (30%), PPP6C (15%), STK19 (10%), LATS1 (8%), ERBB2 (4%), PIK3CA (2%), RAC1 (1%) and N/K/H-RAS (2%). Loss of function (LoF) and deleterious missense mutations were observed in PTPN14 (23%), RB1 (8%) and FBXW7 (5%). In line with the mutational profiles detected by RNAseq, we observed activation of the Hh pathway as well as upregulation of target genes of the Hippo-YAP pathway and activation of MYCN target genes. The functional analysis of the novel tumorigenic driver mutations in MYCN, PTPN14 and LATS1 demonstrates their relevance in BCC tumorigenesis and provides an expanded molecular understanding of BCC.

Study EGAS00001001540

Identification_of_cardiovascular_biomarkers_through_an_integrative_omics_approach

Recent GWAS studies have made extensive use of large eQTL data sets to functionallyannotate index SNPs. With a large number of association signals located outside codingregions there has been an intense search among sequence variants affecting geneexpression at the transcriptional level. However, little progress has been made in mappingregulatory variants that affect protein levels at the translational or post-translational level. It isnow possible to undertake a protein QTL scan for focused sets of e.g. oxidized proteins bymass spectrometry. We have established a collaboration with a longitudinal, family-basedstudy in France, the Stanislas cohort, which comprises circa 1000 nuclear families (4,295individuals) and has follow up data for 10 years (three visits). We have undertaken a pilotstudy in a focus set of 257 subjects from 79 families with the aim to integrate GWAS,transcriptomic and DNA methylation data with proteomic data on a set of 100 proteinsmeasured in PBMCs. We have already generated GWAS data using Illumina's core-exomechip as well as DNA methylation profiles with the 450K array. We propose to use RNA seq togenerate transcriptomic data of the corresponding PBMCs.This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

Study EGAS00001000711

Comprehensive analysis of atypical teratoid rhabdoid tumour (ATRT) using genomic, epigenomic and transcriptomic techniques.

ATRTs (Atypical teratoid rhabdoid tumours) represent one of the most aggressive pediatric brain cancers, but intriguingly exhibit few other recurrently mutated loci except for SMARCB1/hSNF5. We integrated whole genome (n=15), exome, copy number, gene expression and methylation analyses to comprehensively interrogate 64 ATRTs and observed that structural events were relatively frequent in the ATRT genome (~3 tumour). In addition to SMARCB1, which was targeted by structural events in a majority (49/64) of tumours, recurrent structural alterations targeting the LRP1B, CDH13, BCR and MKL1 loci were observed. We observed novel translocation events, including 2 targeting SMARCB1 which were detected only by NGS analyses. Significantly, integration of gene expression and methylation profiles with genomic analyses revealed ATRT comprise two sub-groups with distinct clinical and genetic features. Group 1 tumours were characterized by supra-tentorial brain location, focal intragenic alterations of SMARCB1 and a pro-neural gene expression signature with evidence of NOTCH pathway activation. While Group 2 ATRTs were characterized by infra-tentorial brain location, broad SMARCB1 alterations and an activation of the BMP signalling pathway. Our findings highlight epigenetic mechanisms as important determinants of ATRT tumour phenotypes, and demonstrate for the first time that ATRTs arising in different anatomical compartments comprise distinct molecular and therapeutic sub-groups.

Study EGAS00001000506

JMML targeted sequencing (2013)

Juvenile myelomonocytic leukemia (JMML) is a rare myelodysplastic myeloproliferative neoplasm of early childhood initiated by germline or somatic RAS-activating mutations (Chang, Dvorak et al. 2014). Clinical evolution is heterogeneous, with acute myeloid leukemia (AML) transformation in 1/3 of cases and frequent relapses after hematopoietic-stem-cell transplantation (HSCT), although some patients experience spontaneous remission (Locatelli, Nollke et al. 2005; Matsuda, Shimada et al. 2007; Takagi, Piao et al. 2013; Locatelli and Niemeyer 2015). The initiating lesion incompletely predicts outcome, suggesting the involvement of additional mutations. Genetic profiling of a large JMML cohort (n=118), including whole-exome sequencing (WES) in 30 cases, uncovered additional genetic abnormalities in 56 cases (47%). Somatic events were rare (0.38/case/Mb) and restricted to sporadic (47/78; 60%) or NF1-associated (8/8; 100%) JMML. Multiple concomitant genetic hits were identified targeting the RAS pathway in 13/78 (17%), RAC2D63V, which activates the PI3K/PDK1/AKT and the mTORC2 pathways, and a novel pathway for JMML, the polycomb repressive complex 2 (PRC2), in 26/78 (33%) of the sporadic JMML cases. Interestingly, PRC2 imbalance was recently shown to synergize with RAS activation (De Raedt, Beert et al. 2014). JMML outcome was associated to the mutational profile, suggesting a dose-dependent effect for RAS-pathway activation, distinguishing very aggressive JMML rapidly progressing to AML.

Study EGAS00001001324

Automated system for scoring hematoxylin and eosin-stained ovarian cancer sections by identifying single cells uncovered that stromal cell ratio is a significant predictor for overall survival and progression-free survival.

Concerted efforts in genomic studies have revealed profound insights in prognostic ovarian cancer subtypes. On the other hand, abundant histology slides have been generated to date, yet their uses remain very limited and largely qualitative. Our goal is to develop automated histology analysis as an alternative subtyping technology for ovarian cancer that is cost-efficient and do not rely on DNA quality. We develop an automated system for scoring hematoxylin and eosin-stained (H&E) primary tumour sections of 91 late-stage ovarian cancer to identify single cells including cancer and stromal cells. We demonstrated high accuracy of our system based on expert pathologists’ scores (cancer=97.1%, stromal=89.1%) as well as compared to immunohistochemistry scoring (correlation=0.87). Quantitative stromal cell ratio is significantly associated with poor overall survival after controlling for clinical parameters including debulking status and age (multivariate analysis p=0.0021, HR=2.54, CI=1.40-4.60) and progression-free survival (multivariate analysis p=0.022, HR=1.75, CI=1.09–2.82). We demonstrate how automated image analysis enables objective quantification of microenvironmental composition of ovarian tumours. Our analysis reveals a strong effect of the tumour microenvironment on ovarian cancer progression and highlights the potential of therapeutic interventions that target the stromal compartment or cancer-stroma signalling in the stroma-high, late-stage ovarian cancer subset.

Study EGAS00001001694

Similarity and diversity of the tumor microenvironment in multiple metastases: critical implications for overall and progression-free survival of high-grade serous ovarian cancer.

The tumor microenvironment is pivotal in influencing cancer progression and metastasis. Different cells co-exist with high spatial diversity within a patient, yet their combinatorial effects are poorly understood. We investigate the similarity of the tumor microenvironment of 192 local metastatic lesions in 61 ovarian cancer patients. An ecologically inspired measure of microenvironmental diversity derived from multiple metastasis sites is correlated with clinicopathological characteristics and prognostic outcome. We demonstrate a high accuracy of our automated analysis across multiple sites. A low level of similarity in microenvironmental composition is observed between ovary tumor and corresponding local metastases (stromal ratio r = 0.30, lymphocyte ratio r = 0.37). We identify a new measure of microenvironmental diversity derived from Shannon entropy that is highly predictive of poor overall (p = 0.002, HR = 3.18, 95% CI = 1.51-6.68) and progression-free survival (p = 0.0036, HR = 2.83, 95% CI = 1.41-5.7), independent of and stronger than clinical variables, subtype stratifications based on single cell types alone and number of sites. Although stromal influence in ovary tumors is known to have significant clinical implications, our findings reveal an even stronger impact orchestrated by diverse cell types. Quantitative histology-based measures can further enable objective selection of patients who are in urgent need of new therapeutic strategies such as combinatorial treatments targeting heterogeneous tumor microenvironment.

Study EGAS00001002065

UK10K NEURO ASD TAMPERE

In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. The Tampere Autism sample set consists of samples from Finnish subjects with ASD (autism spectrum disorders) with IQs over 70 recruited from a clinical centre for the diagnosis and treatment of children with ASD. For further information on this cohort please contact either Terho Lehtimaki (terho.lehtimaki@uta.fi) or Kaija Puura (kaija.puura@pshp.fi).

Study EGAS00001000115

MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__pancreas

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

Study EGAS00001002726

Insular Celtic population structure and genomic footprints of migration

Previous studies of the genetic landscape of Ireland have suggested homogeneity, with population substructure undetectable using single-marker methods. Here we have harnessed the haplotype-based method fineSTRUCTURE in an Irish genome-wide SNP dataset, identifying 23 discrete genetic clusters which segregate with geographical provenance. Cluster diversity is pronounced in the west of Ireland but reduced in the east where older structure has been eroded by historical migrations. Accordingly, when populations from the neighbouring island of Britain are included, a west-east cline of Celtic-British ancestry is revealed along with a particularly striking correlation between haplotypes and geography across both islands. A strong relationship is revealed between subsets of Northern Irish and Scottish populations, where discordant genetic and geographic affinities reflect major migrations in recent centuries. Additionally, Irish genetic proximity of all Scottish samples likely reflects older strata of communication across the narrowest inter-island crossing. Using GLOBETROTTER we detected Irish admixture signals from Britain and Europe and estimated dates for events consistent with the historical migrations of the Norse-Vikings, the Anglo-Normans and the British Plantations. The influence of the former is greater than previously estimated from Y chromosome haplotypes. In all, we paint a new picture of the genetic landscape of Ireland, revealing structure which should be considered in the design of studies examining rare genetic variation and its association with traits.

Study EGAS00001002769

Single cell sequencing reveals the origin and the order of mutation acquisition in T-cell acute lymphoblastic leukemia.

Next-generation sequencing has provided a detailed overview of the various genomic lesions implicated in the pathogenesis of T-cell acute lymphoblastic leukemia (T-ALL). Typically, 10 to 20 genomic lesions are found in T-ALL cells at diagnosis. However, it is currently unclear in which order these mutations are acquired and in which progenitor cells this is initiated. To address these questions, we used targeted single-cell sequencing of total bone marrow cells and CD34+CD38- multipotent progenitor cells for 4 T-ALL cases. Hierarchical clustering detected a dominant leukemia clone at diagnosis, accompanied by a few subclones harboring only a fraction of the mutations. We developed a graph-based algorithm to determine the order of mutation acquisition. Two of the four patients had an early event in a known oncogene (MED12, STAT5B) amongst various pre-leukemic events. Intermediate events included loss of 9p21 (CDKN2A/B) and acquisition of fusion genes, while NOTCH1 mutations were typically late events. Analysis of CD34+CD38- cells and myeloid progenitors revealed that in half of the cases somatic mutations were detectable in multipotent progenitor cells. We demonstrate that targeted single cell sequencing can elucidate the order of mutation acquisition in T-ALL and that T-ALL development can start in a multipotent progenitor cell.

Study EGAS00001002830

Sensitive Monogenic Noninvasive Prenatal Diagnosis by Targeted Haplotyping

During pregnancy, cell-free DNA (cfDNA) in maternal blood encompasses a small percentage of cell-free fetal DNA (cffDNA), an easilyaccessible source for determination of fetal disease status in risk families through non-invasive procedures. In case of monogenic heritabledisease, background maternal cfDNA prohibits direct observation of the maternally inherited allele. Non-invasive prenatal diagnostics(NIPD) of monogenic diseases therefore relies on parental haplotyping and statistical assessment of inherited alleles from cffDNA,techniques currently unavailable for routine clinical practice. Here, we present monogenic NIPD (MG-NIPD), which requires a bloodsample from both parents, for targeted locus amplification (TLA)-based phasing of heterozygous variants selectively at a gene of interest.Capture probes-based targeted sequencing of cfDNA from the pregnant mother and a tailored statistical analysis enables predicting fetalgene inheritance. MG-NIPD was validated for 18 pregnancies, focusing on CFTR, CYP21A2, and HBB. In all cases we could predict theinherited alleles with >98% confidence, even at relatively early stages (8 weeks) of pregnancy. This prediction and the accuracy ofparental haplotyping was confirmed by sequencing of fetal material obtained by parallel invasive procedures. MG-NIPD is a robustmethod that requires standard instrumentation and can be implemented in any clinic to provide families carrying a severe monogenicdisease with a prenatal diagnostic test based on a simple blood draw.

Study EGAS00001002622

Using_genetics_to_identify_cell_types_and_mechanisms_underlying_susceptibility_to_primary_sclerosing_cholangitis

Primary sclerosing cholangitis (PSC) is a T-cell mediated, chronic inflammatory condition of the biliary tree that is strongly associated with inflammatory bowel disease. Genome-wide association studies have identified 22 non-HLA genetic risk variants associated PSC. Identifying the genes impacted by these variants has proven difficult as the majority lie in non-coding regions of the genome. Knowledge of the genes and biological pathways these non-coding variants are perturbing is vital to understanding the disease biology. One means of assessing the impact of non-coding variants within disease associated loci upon genes is via colocalisation with eQTL. Many eQTL are cell-type specific, requiring the analysis of disease relevant cell types to detect colocalisation. We have collected PSC-relevant T-cell-subtypes from the peripheral blood of PSC patients via fluorescence activated cell sorting in preparation for RNA sequencing and mapping of eQTL. Samples were collected at the Norfolk and Norwich University Hopital, for which local ethical approval has been granted. Lysed cell samples will be transferred to WTSI and DNA/RNA will be extracted from lysed cell samples by T143 before genotyping (DNA) and custom library preparation and sequencing (RNA). This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/

Study EGAS00001002643

MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Oesophageal_squamous_cell_carcinoma

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

Study EGAS00001002725

An Empirical Approach Leveraging Tumorgrafts to Dissect the Tumor Microenvironment in Renal Cell Carcinoma Identifies Missing Link to Prognostic Inflammatory Factors

We undertook an empirical approach to dissect the tumor microenvironment by developing a bulk tumor dissection algorithm, DisHet, and leveraging RNA-Seq data of tumorgrafts (patient-derived tumors implanted in mice), in which only the tumor cell component expands. We found that approximately 65% of previously defined immune signature genes are not abundantly expressed in the Renal Cell Carcinoma (RCC) microenvironment, and we identified more than two times as many novel immune/stromal transcripts. By using refined immune/stroma-specific genes and genomics, electronic medical record data and imaging data of 1084 RCC patients, we discovered a highly-inflamed pan-RCC subtype enriched for Treg cells, NK cells, Th1 cells, neutrophils, macrophages, B cells, and CD8+ T cells. This inflamed subtype (IS) is enriched for aggressive RCCs, including BAP1-deficient clear-cell and type 2 papillary tumors. Interestingly, IS is correlated with systemic manifestations of inflammation in patients such as thrombocytosis and anemia, whose pathogenesis were poorly defined and have been predictors of poor prognosis in RCC. Indeed, we discovered that IS is a strong predictor of poor survival. Lastly, our analyses show that tumor cells may drive stromal immune response. Overall, these data provide a missing link between the tumor cells, the tumor microenvironment, and systemic factors.

Study EGAS00001002786

13380 results for ""pdc"+"heritage+language"+census+study"

in 56.27 milliseconds.

Center for Common Disease Genomics [CCDG] - Cardiovascular: Genetics of Coronary Heart Disease - Characterizaton of Coronary Prone Pedigrees

Whole exome and transcriptome analysis of UV-exposed epidermis and carcinoma in situ reveals early drivers of carcinogenesis

Human CD4 Memory T Cell Activation Time Course

Single-cell RNA-seq of rheumatoid arthritis synovial tissue using low-cost microfluidic instrumentation

OncoArray: Follow-up of Ovarian Cancer Genetic Association and Interaction Studies (FOCI)

Identification of Cancer Predisposition Genes in Breast Cancer Families

Uveal Melanoma Immunogenomics Predict Immunotherapy Resistance and Susceptibility

NHLBI TOPMed: Best ADd-on Therapy Giving Effective Response (BADGER)

Developing Allelic Imbalance Analysis from Single-Nucleus RNA-Seq Data

National Institute on Aging Genetics of Alzheimer's Disease Data Storage Site (NIAGADS)

The Haplotype-Resolved Genome and Epigenome of the Aneuploid HeLa Cancer Cell Line

Chromatin remodeling enhances MAP3K8 expression in HAM: a key pathogenesis for therapeutic intervention

Inherited chromosomally-integrated human herpesvirus 6A/B (HHV-6A/B) genome sequences in the Japanese population

Treg cell subset-specific gene expression patterns in human head and neck cancer

RBMX functional retrocopy safeguards brain development

Whole genome sequencing of matched primary and metastatic acral melanomas

A thymic ILC1-like progenitor with differentiation potential towards KIR+NKG2A- NK cells

Somatic rearrangements causing oncogenic ectodomain deletions of FGFR1 in squamous cell lung cancer

MRCA and MRCE SNP genotypes

A prognostic human brain network for diffuse midline glioma

DUX4 activates novel intergenic transcripts and isoforms in a tissue-specific context

Phenotypic characterization and prognostic impact of CD103+ tissue-resident memory T cells in diffuse large B cell lymphoma

Liquid biopsy for molecular characterization of diffuse large B-cell lymphoma and early assessment of minimal residual disease

The neo-open reading frame peptides that comprise the tumor framome are a rich source of neoantigens for cancer immunotherapy

Personalized peptide vaccination among 173 patients with IDH wildtype glioblastoma: a retrospective analysis

Organoid BulkRNAseq

A thymic ILC1-like progenitor with differentiation potential towards KIR+NKG2A- NK cells - sc

Raw sequencing files from WGS and RNA-Seq

A Cancer Cell-Line Titration Series for Evaluating Somatic Classification

Transposable Elements Shape Stemness in Normal and Leukemic Hematopoiesis

Transposable Elements Shape Stemness in Normal and Leukemic Hematopoiesis

Aggressive genomic features in clinicallyindolent primary HHV8-negative effusion-based lymphoma

Analysis_of_genomic_integrity_of_disease_corrected_human_induced_pluripotent_stem_cells_by_exome_sequencing

Evolution of DNA repair defects during malignant progression of low-grade gliomas after temozolomide treatment.

Exome and RNA sequencing of relapsed TCF3-PBX1 t(1;19) acute lymphoblastic leukemia

Comparison of 133/144bp dominant phenotype vs 166bp dominant phenotype.

Genomic analysis reveals novel secondary drivers and progression pathways in skin basal cell carcinoma

Identification_of_cardiovascular_biomarkers_through_an_integrative_omics_approach

Comprehensive analysis of atypical teratoid rhabdoid tumour (ATRT) using genomic, epigenomic and transcriptomic techniques.

JMML targeted sequencing (2013)

Automated system for scoring hematoxylin and eosin-stained ovarian cancer sections by identifying single cells uncovered that stromal cell ratio is a significant predictor for overall survival and progression-free survival.

Similarity and diversity of the tumor microenvironment in multiple metastases: critical implications for overall and progression-free survival of high-grade serous ovarian cancer.

UK10K NEURO ASD TAMPERE

MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__pancreas

Insular Celtic population structure and genomic footprints of migration

Single cell sequencing reveals the origin and the order of mutation acquisition in T-cell acute lymphoblastic leukemia.

Sensitive Monogenic Noninvasive Prenatal Diagnosis by Targeted Haplotyping

Using_genetics_to_identify_cell_types_and_mechanisms_underlying_susceptibility_to_primary_sclerosing_cholangitis

MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__Oesophageal_squamous_cell_carcinoma

An Empirical Approach Leveraging Tumorgrafts to Dissect the Tumor Microenvironment in Renal Cell Carcinoma Identifies Missing Link to Prognostic Inflammatory Factors