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

in 42.08 milliseconds.

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

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

  Study EGAS00001003586

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

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

  Study EGAS00001005679

• Targeted sequencing of candidate regions on chromosome 22q predisposing to multiple schwannomas

  Constitutional LZTR1 or SMARCB1 pathogenic variants have been found in ~86% of familial and ~40% of sporadic schwannomatosis cases. Hence, we performed massively parallel sequencing of the entire LZTR1, SMARCB1 and NF2 genomic loci in 35 individuals with schwannomas negative for constitutional first-hit pathogenic variants in the LZTR1/SMARCB1/NF2 coding sequences, however with 22q deletion and/or a different NF2 pathogenic variant in each tumor, including six cases with only one tumor available. Furthermore, we verified whether any other LZTR1/SMARCB1/NF2 (likely) pathogenic variant(s) could be found in 16 cases carrying a SMARCB1 constitutional variant in the 3’ untranslated region (3’-UTR) c.*17C>T, c.*70C>T or c.*82C>T. As no additional variants were found, functional studies were performed to clarify the effect of these 3’-UTR variants on the transcript. The 3’-UTR variants c.*17C>T and c.*82C>T showed pathogenicity by negatively affecting SMARCB1 transcript level. Two novel deep intronic SMARCB1 variants, c.500+883T>G and c.500+887G>A, resulting in out-of-frame missplicing of intron 4, were identified in two unrelated individuals. Further resequencing of the entire repeat - masked genomics sequences of chromosome 22q in individuals negative for pathogenic variants in the SMARCB1/LZTR1/NF2 coding- and non-coding regions revealed five potential schwannomatosis-predisposing candidate genes, i.e. MYO18B, NEFH, SGSM1, SGSM3 and SBF1, pending further verification.

  Study EGAS00001005680

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

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

  Study EGAS00001005737

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

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

  Study EGAS00001007204

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

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

  Study EGAS00001005678

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

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

  Study EGAS00001005235

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

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

  Study EGAS00001005677

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

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

  Study EGAS00001006118

• Prostate cancer ancestral genomic disparity

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

  Study EGAS00001006425

• Whole transcriptome profiling of liquid biopsies from tumour xenografted mouse models: validation cohort

  While cell-free DNA (cfDNA) in liquid biopsies is widely being used and investigated, free circulating RNA (extracellular RNA, exRNA) has the potential to improve therapy response monitoring and cancer detection due to its dynamic nature. A fundamental open question hampering the initiation of large-scale liquid biopsy collections for tumour exRNA analysis is that it remains unclear in which blood subcompartment tumour-derived exRNAs primarily reside. We set out to develop a host-xenograft deconvolution framework, exRNAxeno, with cDNA mapping strategies to either a combined human-mouse reference genome or both species genomes in parallel, that can be applied to exRNA sequencing data from liquid biopsies of human xenograft mouse models, enabling to distinguish (human) tumoural RNA from (murine) host RNA, and as such to specifically analyse tumour-derived exRNA. Subsequently, the preferred exRNAxeno combination pipeline was applied to total exRNA sequencing data from blood platelets and three plasma fractions from a breast cancer patient-derived and neuroblastoma cell line-derived xenograft mouse model. We show that tumoural exRNA concentrations are not determined by plasma platelet levels, while host exRNA concentrations increase with platelet content. Furthermore, a large variability in exRNA levels and gene content across individual mice is observed. In general, the tumoural gene detectability in plasma is correlated with the RNA expression levels in the tumour tissue or cancer cell line. Our results unravel new aspects of tumour-derived exRNA biology in xenograft mouse models.

  Study EGAS00001006582

• deep-learning-powered tissue deconvolution for cfDNA

  Plasma cell-free DNA (cfDNA) is a noninvasive biomarker for cell death of all organs. Deciphering the tissue origin of cfDNA can reveal abnormal cell death because of diseases, which has great clinical potential in disease detection and monitoring. Despite the great promise, the sensitive and accurate quantification of tissue-derived cfDNA remains challenging to existing methods due to the limited characterization of tissue methylation and the reliance on unsupervised methods. To fully exploit the clinical potential of tissue-derived cfDNA, here we present one of the largest comprehensive and high-resolution methylation atlas based on 521 noncancer tissue samples spanning 29 major types of human tissues. We systematically identified fragment-level tissue-specific methylation patterns and extensively validated them in orthogonal datasets. Based on the rich tissue methylation atlas, we develop the first supervised tissue deconvolution approach, a deep-learning-powered model, cfSort, for sensitive and accurate tissue deconvolution in cfDNA. On the benchmarking data, cfSort showed superior sensitivity and accuracy compared to the existing methods. We further demonstrated the clinical utilities of cfSort with two potential applications: aiding disease diagnosis and monitoring treatment side effects. The tissue-derived cfDNA fraction estimated from cfSort reflected the clinical outcomes of the patients. In summary, the tissue methylation atlas and cfSort enhanced the performance of tissue deconvolution in cfDNA, thus facilitating cfDNA-based disease detection and longitudinal treatment monitoring.

  Study EGAS00001007213

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

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

  Study EGAS00001007385

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

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

  Study EGAS00001007555

• WGS and WES data of invasive micropapillary carcinoma of breast

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

  Dataset EGAD00001008450

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

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

  Dataset EGAD00001009512

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

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

  Dataset EGAD00001010881

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

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

  Study EGAS00001005738

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

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

  Study EGAS00001005747

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

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

  Study EGAS00001006879

• H3Africa - Stroke Investigative Research and Education Networks

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

  Study EGAS00001007331

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

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

  Study EGAS00001004739

• Characterization of individual foci of multicentric/multifocal breast cancer using targeted next generation sequencing

  Multifocality or multicentricity in breast cancer may be defined as the presence of two or more tumor foci within a single quadrant of the breast or within different quadrants of the same breast, respectively. This original classification of the breast cancer as multicentric or multifocal was based on the assumption that cancers arising in the same quadrant were more likely to arise from the same ductal structures than those occurring in separate areas of the breast. The problem with these definitions is that the ?quadrants? of the breast are arbitrary external designations, as no internal boundaries do exist. This project will therefore focus both on synchronous multifocal and multicentric tumors. The incidence of multifocal and multicentric breast cancers was reported to be between 13 and 75% depending on the definition used, the extent of the pathologic sampling of the breast and whether in situ disease is considered evidence of multicentricity (1). Although this incidence is variable, those figures show that it is a frequent phenomenon. Multiple (multifocal/multicentric) breast carcinomas, especially when occurring in the same breast, represent a real challenge for both pathologists and clinicians in terms of identifying the cellular origin and the best therapeutic management of the cancer. Multifocality or multicentricity has been associated with a number of more aggressive features including an increased rate of regional lymph node metastases and adverse patient outcome when compared with unifocal tumors (2-3), and a possible increased risk of local recurrence following breast conserving surgery (4). For the moment, the literature is divided on whether there is a corresponding impact on survival outcomes. Today, the current convention to stage and to treat multifocal and multicentric tumors is the classical tumor-node-metastasis (TNM) staging guidelines with which tumor size is assessed by the largest tumor focus without taking other foci of disease into consideration. If some papers, as the recent one from Lynch and colleagues, support the current staging convention (3), others, however, as Boyages et al. suggested that aggregate size and not the size of the largest lesion should be considered in order to refine the prognostic assessment of those tumors (5). On the top of that, the question whether multifocal/multicentric carcinomas are due to the spread of a single carcinoma throughout the breast or is due to multiple carcinomas arising simultaneously has been a matter of debate. Some studies suggested that multifocal breast cancer may result from either intramammary spread from a single primary tumor or multiple synchronous primary tumors; whereas others suggest that multiple breast carcinomas always arise from the same clone (6-8). Recently, Pietri and colleagues analyzed the biological characterization of a series of 113 multifocal/multicentric breast cancers (8) which were diagnosed over a 5-year period. The expression of estrogen (ER) and progesterone (PgR) receptors, Ki-67 proliferative index, expression of HER2 and tumor grading were prospectively determined in each tumor focus, and mismatches among foci were recorded. Mismatches in ER status were present in 5 (4.4%) cases and PgR in 18 (15.9%) cases. Mismatches in tumor grading were present in 21 cases (18.6%), proliferative index (Ki-67) in 17 (15%) cases and HER2 status in 11 (9.7%) cases. Interestingly, this heterogeneity among foci has led to 14 (12.4%) patients receiving different adjuvant treatments compared with what would have been indicated if we had only taken into account the biologic status of the primary tumor. This study therefore showed that differences in biological characteristics of multifocal/multicentric lesions play a crucial role in the adjuvant treatment decision making process. In this study, we will concentrate on a larger series of patients with multifocal invasive ductal breast cancer lesions. We aim at: 1. Evaluating the incidence of multifocality according to the different breast cancer molecular subtypes (ER-/HER2-, HER2+, ER+/HER2-). 2. Evaluating the incidence of multifocality in patients with hereditary breast cancer disease (presence of germline BRCA1 or BRCA2 mutations). Moreover, we would like to investigate if multifocal lesions with BRCA1 or BRCA2 mutations exhibit a characteristic combination of substitution mutation signatures and a distinctive profile of deletions as demonstrated recently by Nik-Zainal and colleagues (9). 3. Correlating multifocality with clinical information in order to define its influence on patients? survival (DFS and OS). 4. Carrying high coverage targeted gene sequencing of driver cancer genes and genes whose mutation is of therapeutic importance in order to compare clinically-relevant genetic differences between several multifocal breast cancer lesions. 5. Evaluating the impact of the distance between the different lesions on the clinical outcome but also on the genetic differences. 6. Comparing gene expression patterns between several multifocal breast cancer lesions and correlate them with the results of the targeted genes screen. 7. Characterizing the genomic and transcriptomic status of cancer related genes in metastatic lesions (local recurrence, positive lymph node or distant metastatic sites) from the same multifocal invasive ductal breast cancer patients in order to evaluate the consequence of genomic and transcriptomic heterogeneity of multifocal lesions on metastatic lesions. Multiple (multifocal/multicentric) breast carcinomas, especially when occurring in the same breast, represent a real challenge for both pathologists and clinicians in terms of identifying the cellular origin and the best therapeutic choice. This project has the potential to identify genetic/transcriptomic differences existing between several lesions constituting multifocal breast cancers, which in the routine clinical practice are usually considered to be homogeneous among them. We foresee validating significant results in a larger series of patients and this, in turn, could have a remarkable impact on the treatment and clinical management of multifocal breast cancers. Indeed, we hope to provide some evidence whether or not each focus matters in multifocal and multicentric breast cancer to define the adequate therapeutic approach, especially in the context of targeted therapies. The work to be done at Sanger will be target gene screen pooling of 1400 samples.

  Dataset EGAD00001000624

• A Single Cell Transcriptomic Analysis of Human Neocortical Development

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

  Study phs001836

• SLCO1B1 Variants and Methotrexate Clearance

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

  Study phs000426

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

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

  Study phs002990

• Common Variation in Candidate Genes in the Diabetes Prevention Program

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

  Study phs000681

• Role of Tobacco Smoke in Clear Cell Renal Cell Carcinoma

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

  Study phs002083

• Strabismus, CCDD and other anomalies

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

  Study phs000478

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

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

  Study phs001982

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

  The NHLBI "Grand Opportunity" Exome Sequencing Project (GO-ESP), a signature project of the NHLBI Recovery Act investment, was designed to identify genetic variants in coding regions (exons) of the human genome (the "exome") that are associated with heart, lung and blood diseases. These and related diseases that are of high impact to public health and individuals from diverse racial and ethnic groups will be studied. These data may help researchers understand the causes of disease, contributing to better ways to prevent, diagnose, and treat diseases, as well as determine whether to tailor prevention and treatments to specific populations. This could lead to more effective treatments and reduce the likelihood of side effects. GO-ESP is comprised of five collaborative components: 3 cohort consortia - HeartGO, LungGO, and WHISP - and 2 sequencing centers - BroadGO and SeattleGO. Large epidemiological studies have demonstrated a significant heritable component in atrial fibrillation (AF), especially the Lone forms, suggesting a monogenic syndrome. Although substantial genetic contribution has been made to the etiology of AF, the specific genes have not yet been identified. The familial form of this disease remains poorly characterized and largely undetermined. Here we seek to identify, characterize and determine the natural course of AF in our clinical practice. We identified four large multi-generation families (FAF 1-4). In FAF 1-2, most family members have symptomatic paroxysmal Atrial Fibrillation (AF) and were adequately treated with a combination of rate and rhythm therapies. By contrast, the AF substrate in FAF 3 and 4 was resistant to anti-arrhythmic drugs and ablation therapies.

  Study phs000362

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

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

  Study phs001853

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

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

  Study phs000846

• Stand Up To Cancer East Coast Prostate Cancer Research Group

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

  Study phs000915

• Molecular Evolution of Cancer

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

  Study phs001255

• Breast Cancer Family Registry

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

  Study phs002835

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

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

  Study phs000239

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

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

  Study phs001165

• Center for Craniofacial and Dental Genetics: Exome Sequencing of Orofacial Cleft Families

  Orofacial Clefts (OFCs) are genetically complex structural birth defects caused by genetic factors, environmental exposures, and their interactions. OFCs are the most common craniofacial anomalies in humans. Affect about 1 in 800 newborns, OFCs are also one of the most common structural birth defects worldwide. OFCs can impose a substantial personal cost resulting from the associated emotional and financial burden. On average a child undergoes 6 surgeries, spends 30 days in hospital, receives 5 years of orthodontic treatment, and participates in ongoing speech therapy, leading to an estimated total lifetime treatment cost of about $200,000. Individuals born with an OFC have higher infant mortality (especially in developing countries with limited access to care), higher mortality rates at all other stages of life, increased incidence of mental health problems, and higher risk for other disorders (notably including breast, brain, and colon cancers). Prior genetic studies have focused on the contribution of common genetic variants through genome-wide linkage and association studies. These studies have revealed a significant etiologic heterogeneity. One promising approach to dissect the etiology of OFCs into focus on subclinical phenotypic features present in the non-cleft family members of affected individuals. These subclinical features are subtle phenotypic features that may represent mild manifestations of the same risk genes; thus, their inclusion in family-based studies may clarify the inheritance patterns in OFC families. In this project, we conducted whole exome sequencing in multiplex families with both overt OFCs and subclinical phenotypes to identify high-impact genes causing OFCs.

  Study phs001675

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

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

  Study phs000482

• National Heart Lung and Blood Institute Exome sequencing in SCID

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

  Study phs000479

• Transposable Elements in FTLD-TDP and ALS-TDP

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

  Study phs001889

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

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

  Study phs002031

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

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

  Study phs003448

• ILyAD (Indolent Lymphoma And vitamin D)

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

  Study phs003503

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

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

  Study phs001486

• Single Cell and Spatial Transcriptomic Profiling of Haemophilus ducreyi Infection

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

  Study phs003754

• Single Chromatin Fiber Profiling in the Human Brain

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

  Study phs003771

• Multiple Myeloma Clinical Targeted Sequencing of Patient Samples

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

  Study phs003908

• dbGaP Collection: Open Translational Science in Schizophrenia (OPTICS)

  Schizophrenia is a chronic, severe, disabling brain disorder that affects approximately 1% of the population worldwide. Epidemiologic studies have clearly demonstrated that genetics play a strong role in etiology, but the inheritance is very complex. Innovative analytic approaches and creative ways of combining disparate data sets will be necessary for making breakthroughs in identifying causal pathways and ultimately new drug targets. This dbGaP Collection consists of all genetic studies of schizophrenia available in dbGaP that have been consented for general research use. The goal is to facilitate identification of datasets with related scientific content in order to expedite the application process and ascertainment of datasets of interest for increased scientific discovery. The Open Translational Science in Schizophrenia (OPTICS) Project, was launched by Janssen Research & Development, LLC, part of the Janssen Pharmaceutical Companies of Johnson & Johnson, with a group of leading research organizations including the National Institutes of Health, Yale University School of Medicine, Rutgers University and the Harvard T.H. Chan School of Public Health, to create a new forum for collaborative analysis of Janssen's schizophrenia clinical trial data and other publicly available data about schizophrenia with the goal of creating new models for conducting research. The OPTICS project is one part of a larger effort at Janssen and other Johnson & Johnson research companies to share clinical trials data to enhance public health and advance science and medicine. Qualified investigators and physicians may apply for access to anonymized clinical trials data from Janssen, for more information please visit https://sites.google.com/site/opticsschizophrenia/.

  Study phs000887