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

in 35.21 milliseconds.

• Single cell RNAseq of stenotic, inflamed and non-inflamed transmural lesions from patients with Crohn's disease

  This dataset contains a gene-cell matrix derived from single-cell RNA sequencing (scRNA-seq) data of ileal tissue from Crohn's disease (CD) patients and colorectal cancer (CRC) patients. It includes: Crohn's Disease Patients: A trio of transmural lesions (stenotic, inflamed, and non-inflamed) from each patient. Colorectal Cancer Patients: Unaffected ileal tissue used as external non-inflamed control. Cell Level Metadata: The dataset includes relevant cell-level metadata such as cell type annotations used in the study. Experimental Details: Platform: 10x Genomics Chromium Single Cell 3' GEX Sequencing: Illumina NovaSeq Processing: Data processed with Cell Ranger software. Resulting count matrices were merged for downstream analysis, including integration and dimensionality reduction. Dataset Composition: Crohn's Disease Patients: 10 patients with 3 samples each (non-inflamed, inflamed, stenotic), totaling 30 samples. Colorectal Cancer Patients: 5 patients with 1 sample each of unaffected tissue, totaling 5 samples. Data Provided: Merged Raw Count Matrix: The final merged raw count matrix used for downstream analysis. Cell Metadata File: Contains details of sample, tissue, and patient for each cell in the count matrix. Barcodes File: Indicate each cell barcode which also encodes the sample, tissue, and patient details for each cell. CD.S_Inf: Stenotic Corhn's disease inflamed samples CD.S_Sten: Stenotic CD patient stenosis sample CD.S_Prox: Stenotic CD Patient - proximal non-inflamed sample CC.C_Prox: CRC Patient proximal unaffected sample eg: A barcode 'CC.C_1_Prox_AAGTCGTAGACCCTTA' indicates CRC Patient unaffected proximal sampe from CRC Patient no.1 and the nucleic acid sequence indicate a unique cell from this sample. Total Samples: Crohn's Disease (CD) Patients: 30 samples Colorectal Cancer (CRC) Patients: 5 samples Patient_no Sample Sample_type 1 CC.C_1 CC.C_1_Prox CC.C_Prox 2 CD.S_1 CD.S_1_Prox CD.S_Prox 3 CD.S_1 CD.S_1_Infl CD.S_Infl 4 CD.S_1 CD.S_1_Sten CD.S_Sten 5 CC.C_2 CC.C_2_Prox CC.C_Prox 6 CD.S_2 CD.S_2_Prox CD.S_Prox 7 CD.S_2 CD.S_2_Infl CD.S_Infl 8 CD.S_2 CD.S_2_Sten CD.S_Sten 9 CC.C_3 CC.C_3_Prox CC.C_Prox 10 CC.C_4 CC.C_4_Prox CC.C_Prox 11 CD.S_3 CD.S_3_Prox CD.S_Prox 12 CD.S_3 CD.S_3_Infl CD.S_Infl 13 CD.S_3 CD.S_3_Sten CD.S_Sten 14 CD.S_4 CD.S_4_Prox CD.S_Prox 15 CD.S_4 CD.S_4_Infl CD.S_Infl 16 CD.S_4 CD.S_4_Sten CD.S_Sten 17 CC.C_5 CC.C_5_Prox CC.C_Prox 18 CD.S_5 CD.S_5_Prox CD.S_Prox 19 CD.S_5 CD.S_5_Infl CD.S_Infl 20 CD.S_5 CD.S_5_Sten CD.S_Sten 21 CD.S_6 CD.S_6_Prox CD.S_Prox 22 CD.S_6 CD.S_6_Infl CD.S_Infl 23 CD.S_6 CD.S_6_Sten CD.S_Sten 24 CD.S_7 CD.S_7_Prox CD.S_Prox 25 CD.S_7 CD.S_7_Infl CD.S_Infl 26 CD.S_7 CD.S_7_Sten CD.S_Sten 27 CD.S_8 CD.S_8_Prox CD.S_Prox 28 CD.S_8 CD.S_8_Infl CD.S_Infl 29 CD.S_8 CD.S_8_Sten CD.S_Sten 30 CD.S_9 CD.S_9_Prox CD.S_Prox 31 CD.S_9 CD.S_9_Infl CD.S_Infl 32 CD.S_9 CD.S_9_Sten CD.S_Sten 33 CD.S_10 CD.S_10_Prox CD.S_Prox 34 CD.S_10 CD.S_10_Infl CD.S_Infl 35 CD.S_10 CD.S_10_Sten CD.S_Sten

  Dataset EGAD50000000559

• Molecular characterization of invasive lobular carcinoma

  Invasive lobular carcinoma (ILC) is the second most common histological subtype of breast cancer accounting for 10-15% of cases. ILC differs from invasive ductal carcinoma (IDC)with respect to epidemiology, histology, and clinical presentation. Moreover, ILC is lesssensitive to chemotherapy, more frequently bilateral, and more prone to form gastrointestinal, peritoneal, and ovarian metastases than IDCs. In contrast to IDC, the prognostic value ofhistological grade (HG) in ILC is controversial. One of the three major components of histological grading (tubule formation) is missing in ILC which hinders the process of gradingin this histological subtype and results in the classification of approximately two thirds of ILC as HG 2.Over the last decade, a number of gene expression signatures have shed light onto breast cancer classification, allowing breast cancer care to become more personalized. Withrespect to the management of estrogen receptor (ER)-positive breast cancer, several gene expression signatures provide prognostic and/or predictive information beyond what is possible with current classical clinico-pathological parameters alone. Nevertheless, most studies using gene expression signature have not considered different histologic subtypesseparately. Recently, a comprehensive research program has elucidated some of the biological underpinnings of invasive lobular carcinoma. Genetic material extracted from 200 ILC tumor samples were studied using gene expression profiling and identified ILCmolecular subtypes. These proliferation-driven gene signatures of ILC appear to have prognostic significance. In particular, the Genomic Grade (GG) gene signature improved upon HG in ILC and added prognostic value to classic clinico-pathologic factors. In addition this study demonstrated that most ILC are molecularly characterized as luminal-A (~75%)followed by luminal-B (~20%) and HER2-positve tumors (~5%). Moreover, we investigated the prognostic value of known gene signatures/ gene modules in the same cohort of ILC. As a second step within the scope of this project, we aim to investigate the interactionsbetween somatic ILC tumor mutations to observed transcriptome findings. To this end, we aim to perform somatic mutation analysis for the ILC tumors for which Affymetrix gene expression profiling is available. To this end, we will use a gene screen assay, which specifically interrogates the mutational status of a few hundreds of cancer genes. We believe that this pioneering effort will be fundamental for a tailored treatment of ILC withimprovement in patients' outcome.

  Dataset EGAD00001000288

• PSCP_bisulphite analysis in hESCs (2018-08-13)

  Many studies over the past 10 years, culminating in the recent report of the International Stem Cell Initiative (ISCI, 2011) have shown that hPSC acquire genetic and epigenetic changes during their time in culture. Many of the genetic changes are non-random and recurrent, probably because they provide a selective growth advantage to the undifferentiated cells. Some are shared by embryonal carcinoma cells, the malignant counterparts of ES cells. The origins of these growth advantages are poorly understood, but may come from altered cell cycle dynamics, resistance to apoptosis or altered patterns of differentiation. Less is known about the nature and consequences of epigenetic changes, but it is likely that these similarly affect hPSC behaviour; e.g., enhanced expression of DLK1, an imprinted gene, is associated with altered hPSC growth (Enver et al 2005). Inevitably, these genetic and epigenetic changes will impact on our ability to use hPSC for regenerative medicine, either because malignant transformation of the undifferentiated cells or their differentiated derivatives to be used for transplantation compromises safety, or because they impede the function of those differentiated derivatives, or because they affect the efficiency with which the undifferentiated cells can be expanded and differentiated into desired cell types. Focusing initially upon the existing clinical grade hESC lines, later moving to iPSC, we will Consolidate and extend knowledge of the rate, type and functional impact of the genetic variations that occur during hPSC culture. We will use whole genome and exome sequencing as well as SNP arrays, together with clonal analysis and other cytogenetics techniques. Common changes will be compared with those found in the normal human population, at low frequency in the original cell population or observed during iPSC generation in the HIPSCI project currently based at the WTSI. These studies will provide a better understanding of the range of genetic changes that occur in hPSC beyond the CNVs already identified. In conjunction with cancer genome resources and expertise at WTSI, bioinformatic analyses of these hPSC data will allow us to assess potential impact on hPSC behaviour pertinent to applications in regenerative medicine, notably the likelihood that specific changes arising in undifferentiated PSC cultures may be associated with potential malignant transformation of differentiated progeny. 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/ . This dataset contains all the data available for this study on 2018-08-13.

  Dataset EGAD00001004295

• MutWP1: CRUK Grand Challenge Mutographs of Cancer: Oesophageal adenocarcinoma

  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 (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. This dataset contains all the data available for this study on 2021-09-27.

  Dataset EGAD00001006083

• Mutational landscape of normal epithelial cells in Lynch Syndrome patients

  Lynch Syndrome (LS) is an autosomal dominant disease conferring a high risk of colorectal cancer due to germline heterozygous mutations in a DNA mismatch repair (MMR) gene. Although cancers in LS patients show elevated somatic mutation burdens, information on mutation rates in normal tissues and understanding of the trajectory from normal to cancer cell is limited. Here we whole-genome sequenced 152 crypts from normal and neoplastic epithelial tissues from LS patients. In normal tissues the repertoire of mutational processes and mutation rates were similar to those found in wild type individuals. A morphologically normal colonic crypt with an increased mutation burden and mutational signatures consistent with MMR deficiency was identified, which may represent a very early stage of LS pathogenesis. Phylogenetic tress of tumour crypts indicated that the most recent ancestor cell of each tumour was already MMR deficient and had experienced multiple clonal evolution cycles. This study demonstrates the genomic stability of epithelial cells with heterozygous germline MMR gene mutations and highlights important differences in the pathogenesis of LS from other colorectal cancer predisposition syndromes.

  Dataset EGAD00001008092

• Comprehensive Genomic Characterization of Refractory Multiple Myeloma Reveals a Complex Mutational and Structural Landscape Associated with Drug Resistance (H067)

  Outcome of patients with multiple myeloma (MM) refractory to immunomodulatory agents (IMiDs) and proteasome inhibitors (PIs) is still very poor. To elucidate the genetic landscape of relapsed/refractory MM (RRMM) we performed whole genome and transcriptome sequencing of 39 heavily pretreated RRMM patients with at least double refractoriness revealing complex structural changes and a high mutational load (H067/K08K). We found a median of 67 structural variants and 116 nonsynonymous exonic mutations per patient. Assessment of oncogenic drivers showed NRAS, BRAF, TP53, KRAS, SLC4A7, FAM46C, DIS3, RB1, KMT2C, MLLT4, EWSR1, HCFC2, and COPS3 as significantly mutated. We found multiple genomic regions with bi-allelic events affecting known tumor suppressor genes, most frequently TP53, RB1, and TRAF3. Supervised analysis of mutational signatures showed a major contribution of Cosmic signature 6 associated with deficiency in DNA mismatch repair and the ‘BRCAness’ signature which coincided with mutations in genes known to confer PARP inhibitor sensitivity. Recurrently mutated gene networks included genes with relevance for PI and IMiD activity, the latter including in particular members of the Cereblon and the COP9 signalosome complex. In conclusion, this comprehensive genomic characterization revealed a complex mutational and structural landscape in RRMM and highlights potential implications for therapeutic strategies.

  Study EGAS00001004363

• Stratifying and Targeting Pediatric Medulloblastoma through Genomics

  In this project, genomic analyses of pediatric medulloblastoma samples, obtained through the international medulloblastoma consortium, will be performed. RNA and miRNA expression profiles of 1000 samples, representing all four subgroups (Wnt, Shh, Group C, and D), will be studied to identify novel subtypes within each subgroup. The resulting subtype-specific expression profiles will support the development of reliable and robust biomarkers to more accurately and reliably classify medulloblastomas for treatment in clinical trials. For that purpose, two assays will be developed: an antibody-based immunohistochemical assay and an orthogonal nucleic acid-based hybridization assay. Additional genomic DNA analysis of the 300 high risk subgroup cases will support the discovery of subgroup specific somatic mutations in order to inform current clinical trials of targeted therapies, and to identify genes and pathways already targeted in other diseases. Such therapies could be rapidly transitioned to Phase II trials in medulloblastoma. Furthermore, the discovery of somatic mutations could be used for developing as well as validating specific biomarkers. The project team will also try to identify risk factors that predispose children to this type of cancer.

  Study EGAS00001000273

• MutWP1__CRUK_Grand_Challenge_Mutographs_of_Cancer__clear_cell_renal_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. 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 EGAS00001002608

• Diverse_outcomes_of_controlled_human_malaria_infection_originate_from_host_intrinsic_immune_variation_and_not_var_gene_switching

  Falciparum malaria is clinically heterogeneous and the relative contribution of parasite and host factors in shaping disease severity remains unclear. We set out to explore the interaction between host inflammation and parasite variant surface antigen (VSA) expression, asking whether this relationship underpins the variation observed in controlled human malaria infection. We uncovered marked heterogeneity in the response of naive hosts to blood challenge with some volunteers maintaining a state of immune quiescence, others triggering interferon-stimulated inflammation and a small group showing transcriptional evidence of myeloid cell suppression. Significantly, only inflammatory volunteers experienced hallmark symptoms of clinical malaria. When we then tracked temporal changes in parasite VSA expression to ask whether variants associated with severe disease preferentially expand in naive hosts (as predicted by current theory) we found that var gene profiles were essentially unchanged after 10-days of blood-stage infection. These data therefore show that the diverse clinical outcomes of CHMI largely originate from host-intrinsic variation and there is no evidence for switching or selection of var genes in naïve hosts. 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 EGAS00001003766

• Identification of biomarkers of response to preoperative talazoparib monotherapy in treatment naïve gBRCA+ breast cancers

  Germline mutations in BRCA1 or BRCA2 exist in ~2-7% of breast cancer patients, which has led to the approval of PARP inhibitors in the advanced setting. We have previously reported a phase II neoadjuvant trial of single agent talazoparib for patients with germline BRCA pathogenic variants with a pathologic complete response (pCR) rate of 53%. As nearly half of the patients treated did not have pCR, better strategies are needed to overcome treatment resistance. To this end, we conducted multi-omic analysis of 13 treatment naïve breast cancer tumors from patients that went on to receive single-agent neoadjuvant talazoparib. We looked for biomarkers that were predictive of response (assessed by residual cancer burden) after 6 months of therapy. We found that all resistant tumors exhibited either the loss of SHLD2, expression of a hypoxia signature, or expression of a stem cell signature. These results indicate that the deep analysis of pre-treatment tumors can identify biomarkers that are predictive of response to talazoparib and potentially other PARP inhibitors, and provides a framework that will allow for better selection of patients for treatment, as well as a roadmap for the development of novel combination therapies to prevent emergence of resistance.

  Study EGAS00001005676

• Multi-organ landscape of therapy-resistant melanoma

  Metastasis and failure of present-day therapies represent the most common causes of mortality in patients with cutaneous melanoma. To identify the underlying genetic and transcriptomic landscapes, here we analyzed multi-organ metastases and tumor-adjacent tissues from eleven rapid autopsies after treatment with MAPK inhibitor (MAPKi) and/or immune checkpoint blockade (ICB) and death due to acquired resistance. Either treatment elicits shared genetic alterations that suggest immune-evasive, cross-therapy resistance mechanisms. Large, non-clustered deletions, inversions, and inter-chromosomal translocations dominate rearrangements. Using 345 therapy-naïve TCGA melanoma and 131 patient-matched melanoma pre-and-post acquired resistance to only either treatment, we performed cross-cohort analyses to identify MAPKi and ICB as respective contributors to gene amplifications and deletions enriched in autopsy versus therapy-naïve tumors. Private/late mutations and structural variants display shifted mutational and rearrangement signatures, with MAPKi specifically selecting for signatures of defective homologous-recombination, mismatch, and base-excision repair. Transcriptomic signatures and cross-talks with tumor-adjacent macroenvironments nominate organ-specific adaptive pathways. An immune-desert, CD8+-macrophage-biased archetype, T-cell exhaustion, and type-2 immunity characterize the immune contexture. This multi-organ analysis of lethal, therapy-resistant melanoma presents preliminary insights with potential to improve therapeutic strategies.

  Study EGAS00001006644

• Institut Curie Neuroblastoma Whole Genome Sequencing Diagnosis Relapse

  Neuroblastoma, a clinically heterogeneous pediatric cancer, is characterized by distinct genomic profiles but few recurrent mutations. As neuroblastoma is expected to have high degree of genetic heterogeneity, study of neuroblastoma's clonal evolution with deep coverage whole-genome sequencing of diagnosis and relapse samples will lead to a better understanding of the molecular events associated with relapse. Samples were included in this study if sufficient DNA from constitutional, diagnosis and relapse tumors was available for WGS. Whole genome sequencing was performed on trios (constitutional, diagnose and relapse DNA) from eight patients using Illumina Hi-seq2500 leading to paired-ends (PE) 90x90 for 6 of them and 100x100 for two. Expected coverage for sample NB0175 100x100bp was 30X for tumor and constitutional samples. For the seven other patients expected coverage was 80X for tumor samples with PE 100x100, 100X in the other tumor samples and 50X for all constitutional samples (see table 1). Following alignment with BWA (Li et al., Oxford J, 2009 Jul) allowing up to 4% of mismatches, bam files were cleaned up according to the Genome Analysis Toolkit (GATK) recommendations (Van der Auwera et al., Current Protocols in Bioinformatics, 2013, picard-1.45, GenomeAnalysisTK-2.2-16). Variant calling was performed in parallel using 3 variant callers: GenomeAnalysisTK-2.2-16, Samtools-0.1.18 and MuTect-1.1.4 (McKenna et al., Genome Res, 2010; Li et al., Oxford J, 2009 Aug; Cibulskis et al., Nature, 2013). Annovar-v2012-10-23 with cosmic-v64 and dbsnp-v137 were used for the annotation and RefSeq for the structural annotation. For GATK and Samtools, single nucleotide variants (SNVs) with a quality under 30, a depth of coverage under 6 or with less than 2 reads supporting the variant were filter out. MuTect with parameters following GATK and Samtools thresholds have been used to filter our irrelevant variants. .SNVs within and around exons of coding genes overlapping splice sites.. Then,variants reported in more than 1% of the population in the 1000 genomes (1000gAprl_2012) or Exome Sequencing Project (ESP6500) have been discarded in order to filter polymorphisms. Finally, synonymous variants were filtered out. MuTect focuses on somatic by filtering with constitutional sample. Mpileup comparison between constitutional and somatic DNAs allowed us to focus also on tumor specific SNVs with GATK and Samtools. Finally, every SNV called by our pipeline and also supported in any constitutional samples were filtered our in order to prevent putative constitutional DNA coverage deficiency. Then we analyzed CNVs (copy number variants) with HMMcopy-v0.1.1 (Gavin et al., Genome Res, 2012) and control-FREEC-v6.7 (Boeva et al., Bioinformatics 2011) with a respective window of 2000bp and 1000 bp, and auto-correction of normal contamination of tumor samples for Control-FREEC. Finally we explored Structural variants (SVs) including deletions, inversions, tandem duplications and translocations using DELLY-v0.5.5 with standard parameters (Rausch et al., Oxford J, 2012). In tumors, at least 10 supporting reads were required to make a call and 5 supporting reads for the sample NB0175 with a coverage of only 40X (see table 2). To predict SVs in constitutional samples for subsequent somatic filtering, only 2 supporting reads were required in order not to miss one. To identify somatic events, all the SVs in each normal sample were first flanked by 500 bp in both directions and any SVs called in a tumor sample which was in the combined flanked regions of respective normal sample was removed (see graph 1). Deletions with more than 5 genes impacted or larger than 1Mb and inversions or tandem duplications covering more than 4 genes, were removed. We focused on exonic and splicing events for deletions, inversions, and tandem duplications. For translocation, we keep all SVs that occurred in intronic, exonic, 5'UTR, upstream or splicing regions. Bioinformatics detection of variations with Deep sequencing approach Once PE reads merged and adaptors trimmed by SeqPrep with default parameters, merged reads were aligned via the BWA (Li H. and Durbin R. 2009 PMID 19451168) allowing up to 1 differences in the 22-base-long seeds and reporting only unique alignments. Only reads having a mapping quality 20 or more have been further analysed. Variant calling software was not used, since we aimed to predict variations at low frequencies, observed in less than 1% of reads. Such variants require a custom approach. Using DepthOfCoverage functions of the Genome Analysis Toolkit (GATK) v2.13.2 (McKenna A, et al., 2010 Genome Research PMID: 20644199), we focused on high quality coverage of bases A, C, G and T at the targeted variant position. Depth of coverage of each base following a mapping quality higher than 20 and a base quality higher than 10 have been taken into account in order to focus only on high quality data. Aiming to determine the background level of variability at the studied regions, 10 control samples were included in the analysis. The same approach and filtering criteria have been applied as introduced above over the entire amplicons. In order to highlight variants, for each sample the frequencies of each bases at each amplicon position were then compared to those observed in the set of controls. Statistical analyses were performed with the R statistical software (http://www.R-project.org). Fisher’s exact two-sided tests with a Bonferroni correction were performed to compare percentages of bases between the data sets, i.e. for a given base between a case and the controls. Finally, significant variations were filtered-in once (i) a significant increase in the percentage of avariant base and (ii) a significant decrease in the percentage of it's reference base following our p.values criteria was observed (p.val < 0.05).

  Dataset EGAD00001001356

• Oncogenic Ectodomain deletion of FGFR1 caused by Breakage-Fusion-Bridges in Squamous Cell lung Cancer

  The identification of recurrent 8p11.23 amplifications including FGFR1 raised the hope of a treatable target in squamous cell lung cancer (SQLC). However, only a minority of patients with FGFR1-amplified tumors respond to single agent inhibitor therapy targeting FGFR. To understand the underlying mechanism of FGFR1 dependency, we performed whole genome and transcriptome sequencing of 25 FGFR1-amplified primary tumors with unknown response upon FGFR inhibition. In addition, we performed deep sequencing of 26 FGFR1-amplified samples whereof the response upon FGFR inhibition was known for 25 samples. In both cohorts we identified intra-chromosomal tail-to-tail breaks close to the FGFR1 transcription start site, being responsible for focal amplification of FGFR1. These specific breaks are caused by a Breakage-Fusion-Bridge-like (BFB-like) mechanism. Here, we associate these breaks with FGFR inhibitor sensitivity. Moreover, in some cases these breaks are located within the open reading frame of FGFR1, which leads to the expression of an ΔEC-FGFR1 transcript that lacks the ecto-domain. Overexpression of ΔEC-FGFR1 transforms Baf3 cells and lead to an FGFR1-dependent phenotype. Our results demonstrate that the truncation of the FGFR1 ectodomain is a frequent event in 8p11.23-amplified squamous cell lung cancer caused by tail-to-tail breaks. These breaks might be used as a predictive therapeutic marker to stratify patients for FGFR-inhibitor therapy.

  Study EGAS00001005059

• Detection of rare mutations, copy number variation, and DNA methylation in the same template DNA molecules

  The analysis of cell-free DNA (cfDNA) from plasma offers great promise for the earlier detection of cancer. At present, changes in DNA sequence, methylation, or copy number are the most sensitive ways to detect the presence of cancer. To further increase the sensitivity of such assays with limited amounts of sample, it would be useful to be able to evaluate the same template molecules for all these changes. Here we report an approach, called MethylSaferSeqS, that achieves this goal, and can be applied to any standard library preparation method suitable for massively parallel sequencing. The innovative step was to copy both strands of each DNA-barcoded molecule with a primer that allows the subsequent separation of the original strands (retaining their 5-methylcytosine residues) from the copied strands (in which the 5-methylcytosine residues are replaced with unmodified cytosine residues). The epigenetic and genetic alterations present in the DNA molecules can then be obtained from the original and copied strands, respectively. We applied this approach to plasma from 265 individuals, including 198 with cancers of the pancreas, ovary, lung and colon, and found the expected patterns of mutations, copy number alterations, and methylation. Furthermore, we could determine which original template DNA molecules were methylated and/or mutated. MethylSaferSeqS should be useful for addressing a variety of questions relating genetics and epigenetics in the future.

  Study EGAS00001006839

• Genetic and Phenotypic Analysis of Multiple Sclerosis in Hispanics

  We sought to identify MS susceptibility variants specific to the Major Histocompatibility Complex and to assess the effect of ancestral risk modification within 2652 Latinx and Hispanic individuals, as well as 2435 Black and African American individuals ascertained in collaboration with the Alliance for Research in Hispanic MS (ARHMS). All DNA samples were obtained through whole blood extraction and were genotyped on the MS Chip, an Illumina Infinium custom genotyping array that contains targeted and dense coverage of the extended MHC, specifically designed for imputation. Classical HLA alleles, SNPs, and amino acid residues were imputed simultaneously from genotyped SNPs using HLA-TAPAS and a multi-ethnic reference panel of 2504 samples from the 1000 Genomes Project. We have identified several novel susceptibility alleles which are rare in European populations including HLA-B*53:01, and we have identified an independent role for HLA-DRB1*15:01 and HLA-DQB1*06:02 on MS risk. We found a decrease in Native American ancestry in MS cases vs controls across the MHC, peaking near the previously identified MICB locus; and we have identified several susceptibility variants for which MS risk is modified by global or local ancestry. We are making data on all genotyped and imputed SNPs, classical HLA alleles, and amino acid residues available through dbGaP for the Hispanic and African American samples utilized in these analyses.

  Study phs003105

• Whole Exome Sequencing and RNA-Seq to Characterize the Somatic Breast Cancer Landscape Among Latinas in California

  Somatic mutational profiles are increasingly being used to guide treatment. However, there is a dearth of data from non-European populations. To address this gap, we conducted paired tumor/normal whole exome paired-end sequencing and tumor RNA-Seq on 146 breast tumors from 140 Hispanic/Latina (H/L) women with breast cancer. We characterized somatic mutations, copy number alterations, tumor subtypes [Prediction Analysis of Microarray 50 (PAM50)], Catalogue Of Somatic Mutations In Cancer (COSMIC) mutation signatures and expression profiles, and compared them to breast tumors from non-Hispanic White (White) women in The Cancer Genome Atlas (TCGA). Using MutSigCV, we found that the most commonly mutated genes, including TP53, RUNX1, PTEN, PIK3CA, MAP3K1, GATA3, CDH1 and CBFB, were similar to those in breast tumors from non-Hispanic white women. Novel findings were that outliers of BLTP2 (aka KIAA0100) over-expression and COSMIC mutational signature 16 were more common in tumors from H/L women than White women. We also identified that the known germline APOBEC3A/B deletion variant was significantly more common in H/L women and was associated with indigenous American ancestry. Data available through dbGaP will be whole-exome sequencing (WES) and RNA sequencing (RNA-Seq) data, as well as phenotypic data, including estrogen receptor, progesterone receptor, and HER2 receptor status, tumor histology, stage, and grade, and age at diagnosis.

  Study phs003218

• BRAIN Cell EncyclOpeDia of Transcribed Elements (BRAINcode)

  BRAINCODE: How Does the Human Genome Function in Specific Brain Neurons? The human brain comprises about 86 billion neurons whose function is central to human biology. How does the human genome program high performing neurons and neural networks in response to experience? What subprograms does the genome express in physiologically and morphologically distinct brain cells? The goal of the BRAIN Cell encyclOpeDia of transcribed Elements Consortium (BRAINcode) is to provide a map of gene expression - both protein-coding and non-coding - in specific cell types, not in culture, but in situ in brains of people. Going beyond traditional mRNA sequencing, polyadenylated and non-polyadenylated transcripts were ultra deeply sequenced using ribo-depleted RNA from neurons laser-captured from human post-mortem brains. Three prototypical neuron types, dopamine neurons, pyramidal neurons, and Betz cells, were prioritized because of their key biologic roles and differential vulnerability to important neurodegenerative diseases such as Parkinson's or Alzheimer's disease. Genetic variation between individuals was examined for correlation with differences in transcribed sequences to identify regions of the genome that influence whether, how, and how much a transcript is expressed in specific cell types in human brains. Our results indicate a vast universe of annotated and novel non-coding RNAs expressed in brain cells and suggest a more diverse and much more complex transcriptional architecture than previously imagined.

  Study phs001556

• Drug screening of patient-derived organoids from colorectal peritoneal metastases

  Colorectal cancer (CRC) patients with peritoneal metastases (CRPM) have limited treatment options and the lowest CRC survival rates. We trialed the possibility of organoid directed precision treatment for CRPM patients. CRPM organoids (peritonoids) isolated from patients underwent next-generation sequencing and medium-throughput drug panel testing ex vivo to identify specific drug sensitivities for each patient. We measured the utility of such a service including: success of peritonoid generation, time to cultivate peritonoids, reproducibility of the medium-throughput drug testing, and documented changes to clinical therapy as a result of the testing. Peritonoids were successfully generated and validated from 68% (19/28) of patients undergoing standard care. Genomic and drug profiling was completed within 8 weeks and a formal report ranking drug sensitivities was provided to the medical oncology team upon failure of standard care treatment. This resulted in a treatment change for 2 patients, one of whom had a partial response despite previously progressing on multiple rounds of standard care chemotherapy. The barrier to implementing this technology in Australia is the need for drug access and funding for off-label indications. In conclusion, our approach is feasible, reproducible and can guide novel therapeutic choices in this poor prognosis cohort, where new treatment options are urgently needed. This platform is relevant to many solid organ malignancies.

  Study phs002023

• eMERGE Geisinger eGenomic Medicine (GeM) - MyCode Project Controls

  A research cohort of adult Geisinger Clinic patients was enrolled from community-based primary care clinics of the Geisinger Health System. Patients were eligible for enrollment if they were a primary care patient of a Geisinger Clinic physician and were scheduled for a non-emergent clinic visit. All participants provided written informed consent and HIPAA authorization. Consenting patient agreed to provide blood samples for broad biomedical research use, and permission to access data in their Geisinger electronic medical record for research. The enrollment rate was 90% of patients approached. The demographics of the cohort approximate those of the Geisinger Clinic outpatient population. Research blood samples were collected during an outpatient clinical phlebotomy encounter. Research blood samples are coded and stored in a central biorepository. Samples are linkable to clinical data in a de-identified manner for research via an IRB-approved data broker process. For genomic analysis, DNA is extracted from EDTA-anticoagulated whole blood. For the initial eMERGE Geisinger eGenomic Medicine (GeM) genotyping project, a subset of 1,232 unique samples were genotyped using Illumina HumanOmniExpress-12 v1.0 arrays, and used as population controls for other Geisinger Clinic case cohorts, including abdominal aortic aneurysm and gastric bypass surgery cases. These samples were selected from a larger subset of approximately 6,000 MyCode DNA samples using a partial matching algorithm that included age, sex, and body mass index as variables.

  Study phs000381

• Genomic Analyses in Neoadjuvant Immunotherapy-Treated Head and Neck Cancers

  This is a phase 2 multicenter interventional trial involving the preoperative (neoadjuvant) administration of pembrolizumab, an anti-programmed cell death protein-1 (anti-PD-1) agent, to patients with pathologically confirmed, surgically resectable advanced head and neck squamous cell carcinoma (HNSCC). Participating centers were Dana-Farber/Harvard Cancer Center, Alvin J. Siteman Cancer Center, Washington University and Memorial Sloan Kettering Cancer Center. The trial (NCT02296684) consisted of two cohorts: Cohort 1 enrolled and treated patients received one dose of pembrolizumab prior to surgical resection with curative intent, and adjuvant pembrolizumab was planned for patients with high risk pathological features (positive margins and/or extranodal extension) in the surgical specimen. Cohort 2 enrolled and treated patients received 2 doses of pembrolizumab prior to surgical resection with curative intent. Adjuvant treatment was not dictated by the trial protocol for Cohort 2. The primary endpoint was rate of at least 50% pathological tumor response (pTR), defined as the presence of tumor necrosis/keratinous debris in the tumor bed following neoadjuvant pembrolizumab. pTR was graded as follows: pTR0 (<10%), pTR1 (10-49%) and pTR2 (≥50%). Baseline, pre-treatment tumor biopsy and post-neoadjuvant pembrolizumab surgically resected tumor specimens were dissociated, and CD3+CD45+ cells were analyzed for T cell receptor (TCR) expression and 5' gene expression analysis using single cell sequencing technologies for a subset of patients.

  Study phs002864

• Epigenomics of Human CD8 T cell Differentiation and Aging

  The efficacy of the adaptive immune response declines dramatically with age, but the cell-intrinsic mechanisms driving the changes characteristic of immune aging in humans remain poorly understood. One hallmark of immune aging is the loss of self-renewing naive cells and the accumulation of differentiated but dysfunctional cells within the CD8 T cell compartment. Using ATAC-seq, we first inferred the transcription factor binding activities that maintain the naive and central and effector memory CD8 T cell states in young adults. Integrating our results with RNA-seq, we determined that BATF, ETS1, Eomes, and Sp1 govern transcription networks associated with specific CD8 T cell subset properties, including activation and proliferative potential. Extending our analysis to aged humans, we found that the differences between memory and naive CD8 T cells were largely preserved across age, but that naive and central memory cells from older individuals exhibited a shift toward a more differentiated pattern of chromatin openness. Additionally, aged naive cells displayed a loss in chromatin openness at gene promoters, a phenomenon that appears to be due largely to a loss in binding by NRF1, leading to a marked drop-off in the ability of the naive cell to initiate transcription of mitochondrial genes. Our findings identify BATF- and NRF1-driven gene regulation as targets for delaying CD8 T cell aging and restoring T cell function.

  Study phs001187

• The Research Institute at Nationwide Children's Hospital Genetics of Congenital Heart Disease (CHD)

  Congenital heart disease (CHD) is the most common type of birth defect, but the cause for the majority of cardiac birth defects remains unknown. Numerous epidemiologic studies have demonstrated evidence that genetic factors likely play a contributory, if not causative, role in CHD. A few genes have been identified by us and other investigators using traditional genetic approaches, but these genes only account for a small portion of the non-syndromic CHDs. Therefore, we are now utilizing whole exome sequencing (WES) and whole genome sequencing (WGS), with the addition of more traditional genetic techniques such as chromosomal microarray or traditional linkage analysis, to identify genetic causes of familial and isolated CHD. With WES/WGS, we are able to sequence all of the genes of an individual and apply different data analysis techniques based on whether we are analyzing a multiplex family or a cohort of trios (mother, father and child with CHD) with a specific isolated CHD. Therefore, WES/WGS is a robust method for identification of novel genetic causes of CHD, which will have important diagnostic and therapeutic consequences for these children. Subjects provided blood samples, skin biopsy, or saliva samples for DNA extraction and creation of cell lines. Tissue remaining after surgical procedures was also collected. Medical records were searched for detailed phenotype data and subjects were followed longitudinally for outcome.

  Study phs002010

• Genomics of Hepatocellular Carcinoma

  Despite the growing incidence of hepatocellular carcinoma (HCC) due to increased hepatitis C virus infection and non-alcoholic steatohepatitis in Western populations, the genetic determinants of this cancer have yet to be established. A minority (15-20%) of HCC occurs in livers without cirrhosis or other chronic disease. Because the liver is functionally preserved in these patients and surgical resection of the tumor may be performed safely, and because of the scarcity of livers available for transplantation, non-cirrhotic patients undergo surgical resection for HCC treatment. These resected tumors present investigatory opportunities in two ways: First, they provide tumor specimen which have not been treated with chemotherapy and embolization. Second, the absence of cirrhosis may provide an unconfounded background from which to investigate the genetic tumorigenesis of HCC. In livers with cirrhosis, genetic instability is prevalent as assessed by assays for microsatellite instability, loss of heterozygosity and aberrant DNA methylation. In contrast, in livers without cirrhosis, the non-tumor tissue surrounding HCC have been found to have fewer genetic aberrations. It has been postulated that the pathologic process of cirrhosis may result in the accumulation of many "passenger"; as well as "driver" mutations, and that the examination of the HCC cancer genome may be facilitated in non-cirrhotic livers because of the absence of an accumulated background noise of mutations.

  Study phs001106

• 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

• 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

• 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

• MutWP1__CRUK_Grand_Challenge__normal_kidney_Nanoseq

  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. Non-tumour DNA is also studied for some cancer types. 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 EGAS00001005451

• Double mutant DNMT3A AML: a unique subtype

  Mutation of DNMT3A, encoding a de novo methyltransferase essential for cytosine methylation, is a common early event in clonal hematopoiesis (CH) and adult acute myeloid leukemia (AML). Spontaneous deamination of methylated cytosines damages DNA, which is repaired by the base excision repair (BER) enzymes MBD4 and TDG. Congenital MBD4-deficiency has been linked to early-onset CH and AML, and is marked by exceedingly high levels of DNA damage and mutation of DNMT3A. Strikingly, wildtype (WT) DNMT3A binds TDG, thereby potentiating its repair activity. Since TDG is the only remaining BER enzyme in MBD4-deficient AML patients capable of repairing methylation damage, we investigated whether mutant DNMT3A negatively affects the repair function of TDG. We found that, whereas WT DNMT3A stimulates TDG function, mutant DNMT3A impairs TDG-mediated repair of DNA damage in vitro. In light of this finding and to extrapolate our observations to the broader AML patient population, we investigate here the genetic profiles and survival outcomes of AML patients with single (SM) versus double mutant (DM) DNMT3A. DM DNMT3A AML patients show a characteristic driver mutation landscape and reduced overall survival when compared to SM DNMT3A AML patients. Importantly, whole-genome sequencing showed a trend for increased DNA damage in primary DM DNMT3A AML samples, especially when DNMT3A mutations are located at the DNMT3A-TDG interaction interface.

  Study EGAS00001007966

• 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

• MutWP1__CRUK_Grand_Challenge__matched_blood_Nanoseq

  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. Non-tumour DNA is also studied for some cancer types. 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 EGAS00001006774

• Papuan Y chromosome Diversity Panel

  Island Southeast Asia and Oceania host one of the world’s richest assemblages of human phenotypic, linguistic and cultural diversity. Despite this, the region’s male genetic lineages are globally among the last to remain unresolved. We compiled ~9.7 Mb of Y chromosome sequence from a diverse sample of over 380 men from this region, including 152 first reported here. The granularity of this dataset allows us to fully resolve and date the regional Y chromosome phylogeny. This new high-resolution tree confirms two main population bursts: multiple rapid diversifications following the region’s initial settlement ~50 kya, and extensive expansions <6 kya. Notably, ~40-25 kya the deep rooting local lineages of C-M130, M-P256, and S-B254 show almost no further branching events in Island Southeast Asia, New Guinea and Australia, matching a similar pause in diversification seen in maternal mitochondrial DNA lineages. The main local lineages start diversifying ~25 kya, at the time of the Last Glacial Maximum. This improved Y chromosome topology highlights localized events with important historical implications, including pre-Holocene contact between Mainland and Island Southeast Asia, potential interactions between Australia and the Papuan world, and a sustained period of diversification following the flooding of the ancient Sunda and Sahul continents as the insular landscape observed today formed. The high-resolution phylogeny of the Y chromosome presented here thus enables a detailed exploration of

  Study EGAS00001006025

• 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

• 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

• 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