690 results for ""PhD+2001"+"Professor+of+Bioinformatics"+2012+18000+citations+profile"

in 31.52 milliseconds.

• Whole Exome Sequencing of Calcitonin Producing Pancreatic Neuroendocrine Neoplasms (CT-pNENs) Indicates a Unique Molecular Signature

  Calcitonin-producing pancreatic neuroendocrine neoplasms (CT-pNENs), which can mimic clinical and laboratory features of medullary thyroid carcinoma, are an extremely rare clinical entity with approximately 60 cases reported worldwide. The molecular profile of CT-pNENs is not yet known. Whole-exome sequencing (WES) was performed on tumor and corresponding serum samples of five patients with increased calcitonin serum levels and histologically proven calcitonin-positive CT-pNENs. cBioPortal analysis and gene enrichment analysis with DAVID were performed to validate candidate genes. Immunohistochemistry was used to detect protein expression of MUC4 and MUC16 in CT-pNEN specimens. Common somatic mutations of pNENs like MEN1, ATRX and PIK3CA were identified in cases of CT-pNENs. New somatic SNVs in ATP4A, HES4, and CAV3 have not yet been described in CT-pNENs. Pathogenic germline mutations in FGFR4 and DPYD were found in three of five cases. Mutations of CALCA (calcitonin) and the corresponding receptor CALCAR were found in all five tumor samples, but none of them resulted in clinical or protein sequelae. All five tumor cases showed single nucleotide variations (SNVs) in MUC4, and four cases showed SNVs in MUC16, both of which were membrane-bound mucins. Immunohistochemistry showed protein expression of MUC4 and MUC16 in one case each, and the liver metastasis of a third case was double positive for MUC4 and MUC16. The homologous recombination deficiency (HRD) score of all tumors was low.CT-pNENs have a unique molecular signature compared to other pNEN subtypes, specifically involving the FGFR4, DPYD, MUC4, MUC16 and the KRT family genes. However, these WES data from only five cases must be interpreted with caution because causal interpretation of the mutational landscape in CT-pNENs is difficult. Further research is needed to explain differences in pathogenesis compared with other pNENs. In particular, multi-omics data such as RNASeq, methylation, and whole genome sequencing could be informative.

  Study phs003060

• Population Architecture using Genomics and Epidemiology (PAGE): Multiethnic Cohort (MEC)

  The Multiethnic Cohort (MEC) has established a large biorepository of blood and urine (N=67,000) and cryopreserved lymphocytes (N=15,000) linked to extensive, prospectively collected risk factors (e.g., diet, smoking, physical activity), biomarkers and clinical data for five racial/ethnic groups. This cohort study of over 215,000 men and women in Hawaii and California is unique in that it is population-based and includes large representations of older adults (45-75 yrs at baseline) for five US racial/ethnic groups (Japanese Americans, African Americans, European Americans, Latinos and Native Hawaiians) at varying risks of chronic diseases. Within the PAGE investigation, the MEC proposes to study: 1) diseases for which we have DNA available for large numbers of cases and controls (breast, prostate, and colorectal cancer, diabetes, and obesity); 2) important cancers that are less common (e.g., lung, pancreas, endometrial cancers, NHL) but for which we propose to pool our data with other funded groups; 3) common traits that are risk factors for these diseases (e.g., body mass index/weight, waist-to-hip ratio, height) and 4) relevant disease-associated biomarkers (e.g., fasting insulin and lipids, steroid hormones). The specific aims are: 1) To determine the population-based epidemiologic profile (allele frequency, main effect, heterogeneity by disease characteristics) of putative causal variants in the five racial/ethnic groups in the MEC; 2) for variants displaying effect heterogeneity across ethnic/racial groups, we will utilize differences in LD to identify a more complete spectrum of associated variants at these loci; 3) investigate gene x gene and gene x environment interactions to identify modifiers; 4) examine the associations of putative causal variants with already measured intermediate phenotypes (e.g., plasma insulin, lipids, steroid hormones); and 5) for variants that do not fall within known genes, start to investigate their relationships with gene expression and epigenetic patterns in small genomic studies. This study is part of the Population Architecture using Genomics and Epidemiology (PAGE) study phs000356.

  Study phs000220

• The Sea Islands Genetic Network (SIGNET): Identifying genetic contributors to diabetes and dyslipidemia in African Americans

  Recent genome-wide association studies (GWAS) have successfully identified genetic variants that influence diabetes risk in European populations, however most do not have a major impact on diabetes risk in populations of African descent. The African American (AA) population from the Sea Islands of coastal South Carolina and Georgia has high rates of type 2 diabetes, low levels of admixture, and in general, consume a diet rich in saturated fats. We postulate that this unique combination of ancestral and environmental factors results in a more consistent penetrance of diabetes risk alleles, as well as enrichment of risk alleles of African origin. The existing DNA samples and rich phenotypic data from the Sea Island Families Project comprise a unique resource for genetic studies of type 2 diabetes and related metabolic traits such as dyslipidemia. Our central hypothesis is that the increased risk for T2DM in AA compared with European American (EA) is due, in part, to susceptibility alleles of African origin, and that these alleles can be identified using a GWAS. The Specific Aims are to: 1) Identify genetic risk factors for type 2 diabetes utilizing DNA samples and data from the Sea Island Families Project, Reasons for Geographic and Racial Differences in Stroke (REGARDS) Study recruited from SC, GA, NC, and AL; and a GWAS approach; 2) Identify genetic contributors to lipoprotein subclasses in African Americans using the lipoprotein subclass profile (particle size and concentration for multiple subclasses of VLDL, LDL, and HDL) assessed by NMR at LipoScience, Inc., and the GWAS data from Aim 1. The rationale for this project is that identification and validation of novel pathophysiological pathways and informed selection of candidate genes for diabetes risk will inform development of new, targeted prevention and treatment strategies in this underserved, high risk population.

  Study phs000433

• Molecular analysis of diffuse cerebellar gliomas

  Among diffuse gliomas, oligodendrogliomas show relatively better prognosis, respond well to radiotherapy and chemotherapy, and seldom progress to very aggressive tumors. To elucidate the genetic and epigenetic background for such behavior and tumor evolution during tumor relapse, we comparatively analyzed 12 pairs of primary and recurrent oligodendrogliomas with 1p/19q-codeletion. Initial treatment for these patients was mostly chemotherapy alone. Temozolomide was used for 3, and procarbazine, nimustine and vincristine (PAV chemotherapy) were used for 7 patients. World Health Organization histological grade at recurrence was mostly stable; it was increased in 2, the same in 9, and decreased in 1 cases. Whole-exome sequencing demonstrated that the rate of shared mutation between the primary and recurrent tumors was relatively low, ranging from 3.2-57.9% (average, 33.3%), indicating a branched evolutionary pattern. The trunk alterations that existed throughout the course were restricted to IDH1 mutation, 1p/19q-codeletion, and TERT promoter mutation, and mutation of the known candidate tumor suppressor genes CIC and FUBP1 were not consistently observed between primary and recurrent tumors. Multiple sampling from different regions within a tumor showed marked intratumoral heterogeneity. Notably, in general, the number of mutations was not significantly different after recurrence, remaining under 100, and no hypermutator phenotype was observed. FUBP1 mutation, loss of chr. 9p21, and TCF12 mutation were among a few recurrent de novo alterations that were found at recurrence, indicating that these events were clonally selected at recurrence but were not enough to enhance malignancy. Genome-wide methylation status, measured by Illumina 450 K arrays, was stable between recurrence and the primary tumor. In summary, although oligodendroglioma displays marked mutational heterogeneity, histological malignant transformation accompanying events such as considerable increase in mutation number and epigenetic profile change were not observed at recurrence, indicating that noticeable temporal and spatial genetic heterogeneity in oligodendrogliomas does not result in rapid tumor progression.

  Study JGAS000106

• Local In Time Statistics for processual research

  Background: Functional genomics in a processual analysis cover the time-dependent changes in transcriptomics and epigenetics before diagnosis of a disease, reflecting the changes in both life style and disease processes. The aim of this paper is to explore the dynamic, time-dependent mechanisms of the metastatic processes, using blood transcriptomics and including time in a continuous manner. For achieving this goal we develop new statistical methods based on statistics that are local in time. Methods: The new statistical method, Local In Time Statistics (LITS), is based on calculating statistics in moving windows and randomization. The method has been tested for the analysis of a dataset that collectively provides information on the blood transcriptome up to eight years before breast cancer diagnosis. The dataset from the NOWAC Post-genome Cohort consists of 467 case-control pairs matched on birth year and time of blood sampling. The data for a pair is the difference in log2 gene expression between the case and control. The stratified analyses are based on important biological differences like metastatic versus non-metastatic cancer, and the mode of cancer detection, i.e. screening detected versus clinically detected cancers. The dataset was used for examining whether the gene expression profile varies between cases and controls, with time, or between cases with and without metastases. Results: The null hypotheses of no differences between cases and controls, no time-dependent changes, and no differences between different strata were all rejected. For screening detected cancers the probability of correct prediction of metastasis status was best in year 1 before diagnosis compared to year 3 and 4 before diagnosis for clinically detected cancers. The predictor was not very sensitive to the number of genes included.Conclusions: Using a new statistical method, LITS, we have demonstrated time-dependent changes of the blood transcriptome up to eight years before breast cancer diagnosis.

  Study EGAS00001002520

• Liver biopsy derived induced pluripotent stem cells provide an unlimited supply for the generation of patient-specific hepatocyte-like cells

  Background: Hepatocyte-like cells (HLCs) differentiated from induced pluripotent stem cells (iPSCs) have emerged as a promising cell culture model to study metabolism, biotransformation, viral infections and genetic liver diseases. However the derivation of iPSCs is restricted by ethical and procedural constraints. Furthermore, incomplete HLC differentiation remains a major challenge. iPSC may carry-on a tissue of origin dependent expression memory influencing iPSC differentiation into different cell types. If liver derived iPSCs (Li-iPSCs) maintain a liver specific gene expression profile is not known. It is also not known, if Li-iPSCs would allow the generation of more fully differentiated HLCs.Methods: Primary liver cells (PLCs) were expanded from liver needle biopsies and reprogrammed into Li-iPSCs using a non-integrative Sendai virus-based system. Li-iPSCs were differentiated into HLCs using established differentiation protocols. The HLC phenotype was characterized at the protein, functional and transcriptional level. RNA sequencing data were generated from the originating liver biopsies, the Li-iPSCs, fibroblast derived iPSCs, and differentiated HLCs, and used to characterize and compare their transcriptome profiles.Results: PLCs were obtained from small pieces of liver biopsies and could be reprogrammed into Li-iPSCs. Li-iPSCs indeed retain a liver specific transcriptional footprint. Li-iPSCs can be propagated to provide an unlimited supply of cells for differentiation into Li-HLCs. Similar to HLCs derived from fibroblasts, Li-HLCs could not be fully differentiated into hepatocytes. Relative to the originating liver, Li-HLCs showed lower expression of liver specific transcription factors and increased expression of genes involved in the differentiation of other tissues.Conclusions: PLCs and Li-iPSCs obtained from small pieces of human needle liver biopsies constitute a novel unlimited source for the production of HLCs. Despite the preservation of a liver specific gene expression footprint in Li-iPSCs, the generation of fully differentiated hepatocytes cannot be achieved with the current differentiation protocols.

  Study EGAS00001002676

• Characterization of genetic intratumor heterogeneity in colorectal cancer and matching patient-derived spheroid cultures.

  Patient-derived in vitro cultures of colorectal cancer (CRC) may help guide treatment strategies prior to patient treatment. However, most previous studies have been performed on a single biopsy per tumor. The purpose of this study was to analyze multiple spatially distinct biopsies from CRCs and see how well intratumor heterogeneity (ITH) was recapitulated in matching patient-derived spheroids. Three to five biopsies were collected from six CRC tumors. Each biopsy was split in two; one half was used for spheroid culturing, the other half was used for direct DNA and RNA purification. For two of the patients, lymph node metastases were analyzed. Somatic mutations were called from whole exome sequencing data. Each tumor contained multiple mutations shared across all biopsies and spheroids, including major CRC drivers such as APC, KRAS, and TP53. At the same time, all tumors exhibited ITH on both mutation and copy number level. The concordance between biopsies and spheroids ranged between 40-70% for coding mutations. For three of the patients the biopsy and spheroid from matching areas clustered together, meaning that the spheroid resembled the area of origin more than the other areas. However, all biopsies and spheroids contained private mutations. Thereby, multiple cultures from spatially distinct sites of the tumor increase the insight into the genetic profile of the entire tumor. RNA sequencing based molecular subtyping was stable and seemed independent of genetic ITH. In conclusion, all examined CRC tumors contained genetic ITH on mutation and copy number level. Spheroid cultures partly reflected this ITH, meaning that multiple cultures from distinct sites of the tumor improved the representation of the different genetic subclones in the primary tumor. This should be taken in to account when establishing patient-derived models for drug screening.

  Study EGAS00001002684

• Somatic copy number alteration and fragmentation analysis in circulating tumor DNA for cancer screening and treatment monitoring in colorectal cancer patients

  Background: Analysis of circulating free DNA (cfDNA) is a promising tool for personalized management of colorectal cancer (CRC) patients. Untargeted cfDNA analysis using whole-genome sequencing (WGS) does not need a priori knowledge of the patient´s mutation profile. Methods: Here we established LIquid biopsy Fragmentation, Epigenetic signature and Copy Number Alteration analysis (LIFE-CNA) using WGS with ~6x coverage for detection of circulating tumor DNA (ctDNA) in CRC patients as a marker for CRC detection and monitoring. Results: We describe the analytical validity and a clinical proof-of-concept of LIFE-CNA using a total of 259 plasma samples collected from 50 patients with stage I-IV CRC and 61 healthy controls. To reliably distinguish CRC patients from healthy controls, we determined cutoffs for the detection of ctDNA based on global and regional cfDNA fragmentation patterns, transcriptionally active chromatin sites, and SCNAs. We further combined global and regional fragmentation pattern into a machine learning classifier to accurately predict ctDNA for cancer detection. By following individual patients throughout their course of disease, we show that LIFE-CNA enables the reliable prediction of response or resistance to treatment up to 3.5 months before commonly used CEA. Conclusion: In summary, we developed and validated a sensitive and cost-effective method for untargeted ctDNA detection at diagnosis as well as for treatment monitoring of all CRC patients based on genetic as well as non-genetic tumor-specific cfDNA features. Thus, once sensitivity and specificity have been externally validated, LIFE-CNA has the potential to be implemented into clinical practice. To the best of our knowledge, this is the first study to consider multiple genetic and non-genetic cfDNA features in combination with ML classifiers and to evaluate their potential in both cancer detection and treatment monitoring.

  Study EGAS00001006490

• Separation of chronic myeloid leukemia stem cells from normal hematopoietic stem cells at single-cell resolution

  Chronic myeloid leukemia (CML) is a clonal disease, initiated and sustained by a subpopulation of CML stem cells (CML-SC) containing the BCR::ABL1 oncogene. Genetic markers that enrich for CML-SC have been identified, including CD26 being most specific. Genome sequencing analysis remain yet to be performed at single cell level in order to determine differences and heterogeneity in the genetic profile of CML-SC versus normal hematopoietic stem cells (HSC). Here we show that CD26 and CD33 separate CML-SC from normal HSC in freshly isolated bone marrow samples from newly diagnosed CML patients. DNA sequencing confirms that single CML-SC derived colonies are positive for BCR::ABL1 and harbor a recurrent set of patient-specific clonal somatic mutations across different CML-SC derived colonies. In contrast, marker-negative (CD26-CD33-) HSC from the same patients are negative for the BCR::ABL1 oncogene and do not carry recurrent, clonal mutation across different colonies. These results support the concept that CD26/CD33 expression separates genetically defined, clonal mutant CML-SC from normal HSC that carry individual sets of random mutations in CML patients. We provide whole exome DNA sequencing data for 23 clinical samples of CML-SC plus two buccal swipes derived normal samples. CML samples are comprised of 4 to 8 replicates from two patients, at diagnosis and after treatment. Single CML stem cells before treatment and single non-transformed HSC at remission were selected from bone marrow samples by FACS, according to genetic markers CD33+CD26+ at diagnosis and CD33+CD26-/CD33-CD26- at remission. WES libraries of colony forming assays derived CML-SC and HSC populations were prepared using Agilent SureSelect Human All Exon V6 kit and sequenced running 150 cycles (2x 75bp paired-end) on an Illumina NextSeq 500 platform. Following GATK best practices for somatic mutation calling, novel CML variants and random mutation events during remission phase were obtained.

  Study EGAS00001006904

• The Melbourne Urological Research Alliance (MURAL) Collection of Patient-Derived Models of Prostate Cancer

  1. Low-Dose Carboplatin Modifies the Tumor Microenvironment to Augment CAR T Cell Efficacy in Human Prostate Cancer Models.Single-cell RNA sequencing (scRNA-seq) was performed to analyze the transcriptional modifications in the tumor microenvironment of prostate cancer patient-derived xenografts (PDX). The analyzed PDX tumor was PDX-287R from the MURAL collection of PDXs (PMID: 34413304). Mice harboring PDX tumors were treated with the carboplatin chemotherapy, and tumors were harvested 3 weeks later. RNA sequencing (RNA-Seq) was used to define immune and non-immune cell populations and anti-tumor innate immune signaling within the tumor microenvironment.To investigate the response of tumor cells to carboplatin treatment in more detail, RNA-Seq was performed on FACs-isolated tumor cells from PDX-287R at one week post-treatment. Gene set enrichment analysis showed an enrichment for apoptotic pathways, with significantly increased expression of the pro-apoptotic genes FAS, BAX, BBC3, IFI6, and JUN. Previously, carboplatin has been shown to activate the cGAS/STING pathway, and here we also found an increase in STING signaling, with an enrichment for the cytosolic DNA sensing pathway and increased expression of STING1, STAT1 and STAT2. There was also a significant increase in the expression of genes involved in T cell chemotaxis, including CXCL10, CXCL11, and CCL20. Collectively, carboplatin-induced cell death in tumor cells likely initiated a pro-inflammatory phenotype in prostate cancer PDXs.2. Co-Targeting BET, CREBBP, and EP300 Inhibits Neuroendocrine Signaling in Androgen Receptor-Null Prostate Cancer.The landscape of castration-resistant prostate cancer (CRPC) is characterized by a multitude of phenotypes, among which neuroendocrine disease holds a prominent position, each displaying unique responses to therapeutic interventions. The efficacy of inhibitors targeting BET and CREBBP/EP300 in prostate cancer treatment is widely acknowledged, primarily due to their ability to attenuate androgen receptor (AR) signaling. However, the effectiveness of these inhibitors in prostate cancers devoid of AR remains uncertain.This investigation sought to elucidate the role of BRD4, CREBBP, and EP300, which are co-expressed in both AR-positive and AR-null prostate cancer. Our study revealed that a compound targeting all three proteins, NEO2734, exhibited efficacy in suppressing the growth of organoids from both AR-positive and AR-null tumors, as evidenced by alterations in viability, size, and cellular composition. Furthermore, NEO2734 treatment consistently led to transcriptional downregulation of cell cycle pathways across various models. In the context of neuroendocrine prostate cancer (NEPC), NEO2734 demonstrated notable efficacy. Treatment with NEO2734 resulted in decreased expression of ASCL1 and other neuroendocrine markers (including SYP) , accompanied by a reduction in tumor growth in vivo. These findings collectively underscore the potential of epigenome-targeted inhibitors in impeding the growth of neuroendocrine prostate cancer, by disrupting lineage regulators in a phenotype-dependent manner. Consistent with the NEPC organoids and explants, NKX3-1 was upregulated in NEO2734-treated PDXs.In this study, we examined the impact of NEO2734 treatment on 6 organoids (treated for 24 hrs), 3 explants treated with JQ1 (for 48 hrs) and 2 PDXs (at two time points) using bulk RNA sequencing (RNA-Seq). These results are available through this dbGaP submission. The implications of these results warrant further exploration and development of compounds possessing similar activities for clinical applications. The ability of NEO2734 to mitigate the growth of both AR-positive and AR-null prostate cancer underscores its potential as a versatile therapeutic agent. Moreover, its efficacy in targeting neuroendocrine phenotypes offers promising avenues for addressing a subtype of CRPC that poses significant clinical challenges. In conclusion, this study sheds light on the intricate interplay between epigenetic regulators and prostate cancer phenotypes. The observed effects of NEO2734 on growth inhibition and phenotype-specific disruption highlight the importance of continued research into compounds with similar epigenome-targeting mechanisms for advancing precision medicine approaches in the treatment of prostate cancer.3. Defining Focal Neuroendocrine Differentiation as a Transcriptionally Distinct Form of Prostate Cancer Pathology Characterized by the Expression of Androgen ReceptorsThe study aimed to understand the transcriptional differences among different neuroendocrine pathologies of prostate cancer. The study hypothesized that the diverse neuroendocrine pathology observed in prostate might be attributable to differences in transcriptional level that remain hidden in bulk RNA-Seq data. To investigate, we employed scRNA-seq to generate gene expression profiles of 18,632 individual tumor cells from 9 PDX models representing five distinct neuroendocrine pathologies of prostate cancer. Our analysis identified 3-8 transcriptionally distinct sub-populations per PDX. The expression of key oncogenic signaling pathways and master regulator activity varied across in neuroendocrine pathologies, with each type of pathology displaying a unique set of transcriptional sub-populations. We discovered that, like the amphicrine pathology, focal neuroendocrine differentiation (focal NED) cells maintain androgen receptor (AR) signaling. However, the expression profiles of focal NED cells differed from the other pathologies, indicating it should be considered as a distinct type of neuroendocrine pathology in prostate cancer. Analysis of copy number alterations suggested little clonal divergence between focal NED cells and neighboring adenocarcinoma cells, emphasizing the transcriptional distinctiveness of focal NED.Overall, the findings suggest potential differences in treatment approaches for tumors from prostate cancer patients showing neuroendocrine pathology, especially focal NED.4. Reprogramming of Androgen Receptor Activity in Castration-Resistant Prostate Cancer is Shaped by Truncated Variants.Under the influence of treatments that target the androgen receptor (AR) axis, prostate cancer cells undergo various changes to the AR gene. This includes the production of truncated AR variants that lack the ligand binding domain and are constitutively active. However, the exact role of these truncated variants in shaping the AR signaling axis and their importance in therapy resistance is not yet fully understood. By examining the AR cistrome in a group of patient-derived prostate cancer models with different mechanisms of castration resistance, we observed diverse changes in AR activity in metastatic prostate cancer, including castration-resistant models. Specifically, we identified a unique subset of tumors characterized by the expression of ARv567es, a variant resulting from structural rearrangements of the AR gene. Tumors positive for ARv567es showed a distinct AR binding profile and epigenome compared to those without ARv567es, which was linked to an altered transcriptional profile. The presence of ARv567es is associated with poor responses to treatments, such as castration and bipolar androgen therapy (BAT), regardless of the presence of full-length AR. This involves disruption of the autoregulatory loop that controls AR gene transcription and reduced transcriptional responses to treatment. Overall, this study demonstrates that ARv567es can lead to transcriptional reprogramming in advanced prostate cancer and drive therapy resistance. 5. The MURAL Collection of Prostate Cancer Patient-Derived Xenografts Enables Discovery Through Preclinical Models of Uro-Oncology.We created a comprehensive PDX resource for prostate cancer to facilitate the rapid and systematic evaluation of new therapies. This collection includes 59 tumors from 30 patients, gathered between 2012 and 2020, aligning with the availability of abiraterone and enzalutamide in the clinic. The PDXs cover the clinical, pathological, and genomic spectrum of prostate cancer, ranging from treatment-naïve primary tumors to castration-resistant metastases. We characterized heterogeneity in adenocarcinoma and neuroendocrine phenotypes using bulk and single-cell RNA sequencing. Additionally, organoids were cultured from PDXs, enhancing our preclinical study capabilities. Employing a 1x1x1 design, we quickly identified tumors with responses to combination therapies. To manage the distribution of PDXs, we established the Melbourne Urological Research Alliance (MURAL). This PDX collection significantly enhances our ability to test and prioritize effective treatments for future clinical trials in prostate cancer.6. Prostate Cancer Associated Fibroblasts Have Distinct Morpho-Mechanical Features That Are Associated With Patient Outcome.The tumor stroma, including cancer-associated fibroblasts (CAFs), has a crucial role in tumor progression. This study focused on the morphological and mechanical properties of CAFs and normal prostatic fibroblasts. The findings revealed that CAFs had distinct morphological and mechanical characteristics compared to normal fibroblasts, such as differences in nuclear size, shape, F-actin arrangement, cellular volumes, and elasticity. Additionally, the study correlated these biophysical properties with transcriptomic data, identifying pathways and cellular components linked to the observed changes. Overall, high-throughput assessments of biophysical properties in prostate cancer cells and stromal components may serve as predictive tools for patient outcomes and help identify new therapeutic approaches targeting the tumor stroma.

  Study phs003369

• Genomic profile of sporadic multiple meningiomas

  Molecular mechanisms underlying the formation of sporadic MMs and clonal formation etiology of these tumors are poorly understood. In this study, we worked on a cohort of patients with spatially separated MMs without prior radiation exposure or a family history who underwent surgical resection of at least two meningiomas. Unbiased, comprehensive next generation sequencing was performed and relevant clinical data was analyzed. We identified that both NF2-loss and non-NF2 driven MMs can form due to monoclonal expansion and those tumors can acquire inter-tumoral heterogeneity through branched evolution. Grade I and II meningiomas can occur in the same patient. Thus, the molecular make-up and clinical behavior of one tumor in MMs, cannot reliably lend insight into that of the others and suggests the clinical management strategy for MMs should be tailored individually.

  Study EGAS00001005700

• Somatic mutation in edited cholangiocyte organoids

  Primary human cells cultured in organoid format have great promise as potential regenerative cellular therapies. However, their immunogenicity and mutagenic profile remain unresolved, impeding effective long-term translation to the clinic. In this study we report, for the first time, the generation of human leukocyte antigen (HLA)-I and HLA-II knock-out human expandable primary cholangiocyte organoids (PCOs) using CRISPR-Cas9 as a potential ‘universal’ low-immunogenic therapy for bile duct disorders. HLA-edited PCOs (ePCOs) displayed the same phenotypic and functional characteristics as parental unedited PCOs. Despite minimal off-target edits, duplex sequencing approaches demonstrated that ePCOs and PCOs acquire mutations in culture at similar rates, but without evident selection for cancer-driver mutations. ePCOs induced reduced T cell-mediated immunity and donor-dependent NK cell cytotoxicity in vitro and evaded cytotoxic responses with increased graft survival in humanized mice in vivo. Our findings have important implications for assessment of safety and immunogenicity of primary cell-derived organoid cellular therapies.

  Dataset EGAD00001015455

• Somatic mutation in edited cholangiocyte organoids

  Primary human cells cultured in organoid format have great promise as potential regenerative cellular therapies. However, their immunogenicity and mutagenic profile remain unresolved, impeding effective long-term translation to the clinic. In this study we report, for the first time, the generation of human leukocyte antigen (HLA)-I and HLA-II knock-out human expandable primary cholangiocyte organoids (PCOs) using CRISPR-Cas9 as a potential ‘universal’ low-immunogenic therapy for bile duct disorders. HLA-edited PCOs (ePCOs) displayed the same phenotypic and functional characteristics as parental unedited PCOs. Despite minimal off-target edits, duplex sequencing approaches demonstrated that ePCOs and PCOs acquire mutations in culture at similar rates, but without evident selection for cancer-driver mutations. ePCOs induced reduced T cell-mediated immunity and donor-dependent NK cell cytotoxicity in vitro and evaded cytotoxic responses with increased graft survival in humanized mice in vivo. Our findings have important implications for assessment of safety and immunogenicity of primary cell-derived organoid cellular therapies.

  Dataset EGAD00001015456

• Somatic mutation in edited cholangiocyte organoids

  Primary human cells cultured in organoid format have great promise as potential regenerative cellular therapies. However, their immunogenicity and mutagenic profile remain unresolved, impeding effective long-term translation to the clinic. In this study we report, for the first time, the generation of human leukocyte antigen (HLA)-I and HLA-II knock-out human expandable primary cholangiocyte organoids (PCOs) using CRISPR-Cas9 as a potential ‘universal’ low-immunogenic therapy for bile duct disorders. HLA-edited PCOs (ePCOs) displayed the same phenotypic and functional characteristics as parental unedited PCOs. Despite minimal off-target edits, duplex sequencing approaches demonstrated that ePCOs and PCOs acquire mutations in culture at similar rates, but without evident selection for cancer-driver mutations. ePCOs induced reduced T cell-mediated immunity and donor-dependent NK cell cytotoxicity in vitro and evaded cytotoxic responses with increased graft survival in humanized mice in vivo. Our findings have important implications for assessment of safety and immunogenicity of primary cell-derived organoid cellular therapies.

  Dataset EGAD00001015457

• DERMATLAS: Hidradenoma papilliferum_RNAseq

  Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible.For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previously, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles.

  Dataset EGAD00001015480

• DERMATLAS: Hidradenoma papilliferum_WES

  Virtually no other tumour type is associated with so many different forms as skin cancer. Histologically, tumours of the skin may arise from epithelium, including epidermis, hair follicle, sebaceous or sweat gland, melanocytes, dermal-associated mesenchymal structures or tissue resident immune cells, making for a diversity of clinical presentations. Importantly, many skin tumour types have an extremely poor prognosis. Many skin tumour subtypes have never undergone molecular profiling, or if they have, targeted sequencing has been used and the number of cases analyzed has been so limited that deriving firm conclusions about the profile of driver genes, DNA mutational signatures and germline alleles has not been possible.For 70 key skin tumour subtypes defined by the World Health Organization, that have not undergone extensive genetic analysis previously, we propose to perform whole exome and transcriptome sequencing, from a range of body sites, to build a genomic atlas of dermatological tumours, including detailed maps of SNVs, copy number alterations, genome-wide methylation and expression profiles.

  Dataset EGAD00001015481

• TRanscriptomic ANalySis of left ventriCulaR gene Expression (TRANSCRibE)

  Myocardial ischemia occurs when there is a mismatch between coronary oxygen delivery and metabolic requirements of the myocardium, which may be clinically manifested during angina, coronary angioplasty or cardiopulmonary bypass (CPB). Myocardial ischemia may lead to a spectrum of myocardial stunning, hibernating myocardium, and ultimately cell death if the ischemic insult is severe. In the human heart, irreversible damage begins after approximately 20 to 40 minutes of oxygen deprivation. Observed molecular and cellular changes of myocardial ischemia are characteristic of an inflammatory response, but the exact mechanisms that underlie this pathological process are unclear and may not be full emulated by animal models of ischemia or infarction. Thus, we felt it valuable to investigate a human ischemia model. During cardiac surgery, CPB with aortic cross-clamping (AoXC) and cardioplegic arrest is associated with excellent clinical outcomes and suitable operative conditions. However, despite the use of cardioprotective strategies, AoXc during CPB is accompanied by a variable, yet obligate ischemic period lasting 1 to 3 hours, resulting in hypoxia, metabolic substrate depletion, reperfusion injury, apoptosis, and necrosis. Cardiac specific biomarkers of ischemia and infarction, including troponin, are elevated even after routine coronary artery bypass graft surgery and correlate with the duration of ischemia from AoXc.This process of CPB provides us with the ability to examine the transcriptional profile before and after an expected, consistent, and reproducible human ischemic event, albeit induced by cold cardioplegic arrest and not coronary occlusion. In addition, the absence of reperfusion in this time period allows us to examine the transcriptomic response to intermittent ischemia, without having to account for the perturbations of reperfusion injury. Although various animal models have been used to examine the effects of ischemia on cardiac function, no human data exist which examine the early transcriptomic response to a left ventricular (LV) ischemic insult. We therefore characterized the effect of cold cardioplegia induced acute ischemia on the transcriptional profile of the LV by performing whole transcriptome next-generation RNA-sequencing (RNA-seq) in patients undergoing cardiac surgery by sampling human LV tissue prior to, and after, the obligate ischemia during AoXC. We hypothesized that the cold cardioplegia induced ischemic injury will dramatically alter transcription in the human myocardium, and that we would identify genes and pathways, which will identify interventional targets for pharmacological therapy. Methods:We have collected left ventricle tissue samples and blood sample from patients undergoing heart surgery. We obtained punch biopsies (~3-5μg total RNA content) from the site of a routinely placed surgical vent in the anterolateral apical left ventricular wall of patients undergoing elective aortic valve replacement surgery with cardiopulmonary bypass. After an average of 79 minutes of aortic cross-clamping with intermittent cold blood cardioplegia for myocardial protection every 20 minutes, a second biopsy was obtained in the same manner. Tissue samples were immediately placed in RNAlater® (Ambion, ThermoFisher Scientific, Waltham, MA), and after 48 hours at +4°C were stored at -80°C until RNA extraction. Total RNA was isolated with Trizol and RNA quality was assessed using the Agilent Bioanalyzer 2100 (Agilent, Santa Clara, CA). Libraries were prepared by poly(A) mRNA isolation and reverse transcription Polymerase Chain Reaction (RT-PCR), then sequenced on the Illumina HiSeq2000 or HiSeq2500 (Illumina, San Diego, CA). As samples were analyzed at different times, different read lengths were employed, initially using single-end reads (n=20) and then transitioning to paired end reads (n=216), ranging from 36 - 100 base pairs. Raw sequencing files were processed using Sickle, Skewer, and STAR software, and aligned to GrCh37 or UCSC Hg19. DNA was isolated from whole blood using standard methods. SNP genotyping was performed using the Illumina Omni2.5Exome-8 BeadChip array with additional exome content (Illumina, San Diego, CA) chip, version 1.1. We first phased and imputed 93 subjects using a phasing tool called SHAPEIT and an imputation tool called MINIMAC, with 1000 Genomes phase 1 version 3 for the reference panel. We then phased and imputed 26 more subjects using SHAPEIT, an imputation tool called IMPUTE2, and 1000 Genomes phase 3 version 5.

  Study phs001679

• Identifying transcriptional programs underlying anti-EGFR small molecule response and resistance with TraCe-seq

  Genetic and non-genetic heterogeneity within cancer cell populations represents a major challenge to anti-cancer therapies. We currently lack robust methods to determine how pre-existing and adaptive features affect cellular responses to therapies. Here, by conducting clonal fitness mapping and transcriptional characterization using expressed barcodes and single-cell RNA-sequencing, we have developed TraCe-seq, a method that captures at clonal resolution the origin, fate, and differential early adaptive transcriptional programs of cells in a complex population in response to distinct treatments. We used TraCe-seq to benchmark how next-generation dual EGFR inhibitors-degraders compare to standard EGFR kinase inhibitors in EGFR-mutant lung cancer cells. To QC the TraCe-seq strategy, single-cell RNA-seq libraries were generated from a variety of human cancer cell lines transduced with the TraCe-seq library to validate the TraCe-seq strategy. Specifically, 5 different cell lines (PC9, MCF-10A, MDA-MB-231, NCI-H358, and NCI-H1373) were each transduced with a unique TraCe-seq barcode. The transduced cells were selected with puromycin only, dissociated to single cell suspensions, and then mixed together. The complex mixture of the 5 cell lines was profiled by 10X scRNA-seq. Furthermore, transduced NCI-H1373 cells were sorted by FACS to enrich for the top 50% of eGFP positive cells, and sorted cells were cultured briefly and used to construct scRNA-seq libraries and profiled by 10x scRNA-seq. To carry out the full TraCe-seq experiment, ~600 PC9 cells carrying unique TraCe-seq barcodes were expanded over 12 doublings to establish the barcoded population. A subset of the barcoded PC9 population was used to generate scRNA-seq libraries and profiled by 10x scRNA-seq prior to treatment to establish a baseline transcription profile for each barcoded clone. The rest of the cells were then treated for four days with 1 µM erlotinib, 1 µM GNE-069, or 1 µM GNE-104 respectively. scRNA-seq libraries were then generated form the treated cells and profiled by 10x scRNA-seq.

  Study EGAS00001005405

• Efficient prediction of a spatial transcriptomics profile better characterizes breast cancer tissue sections without costly experimentation

  Spatial transcriptomics is an emerging technology requiring costly reagents and considerable skills, limiting the identification of transcriptional markers related to histology. Here, we show that predicted spatial gene-expressions in unmeasured regions and tissues can enhance biologists��� histological interpretations. We developed the Deep learning model for Spatial gene Clusters and Expression, DeepSpaCE and confirmed its performance using the spatial-transcriptome profiles and immunohistochemistry images of consecutive human breast cancer tissue sections. For example, the predicted expression patterns of SPARC, an invasion marker, highlighted a small tumor-invasion region that is difficult to identify using raw data of spatial transcriptome alone because of a lack of measurements. We further developed semi-supervised DeepSpaCE using unlabeled histology images and increased the imputation accuracy of consecutive sections, enhancing applicability for a small sample size. Our method enables users to derive hidden histological characters via spatial transcriptome and gene annotations, leading to accelerated biological discoveries without additional experiments.

  Study JGAS000290

• Single-cell Kinnex sequencing of Alzheimer's disease isoform profile

  Alzheimer’s disease (AD) is the most common cause of dementia, yet its molecular causes remain incompletely understood. RNA diversity in part arising from dysregulated splicing may contribute to AD pathogenesis; however, the ability to interrogate the resulting full-length isoforms in single brain cells has been technologically limited. Here we report the use of PacBio Kinnex long-read sequencing combined with 10X Genomics single-cell preparations to identify both annotated and unannotated RNA isoforms. Eight AD and seven non-diseased post-mortem human brains yielded ~70,000 single nuclei showing diverse, differentially expressed, and switched transcripts in multiple genes. Cell-type-specific isoform expression and variants with intra-exonic junctions associated with reverse transcriptase-mediated somatic gene recombination were also detected. Novel isoforms, including CHI3L1 and SEPTIN4, were altered in AD. Kinnex sequencing of RNA isoforms from single nuclei detected vast isoform diversity amongst brain cell types, representing an under-explored element in AD and other brain disorders.

  Study EGAS50000001476

• Differential expression profile of gluten-specific T cells identified by single-cell RNA-seq

  Gluten-specific CD4+ T cells drive the pathogenesis of celiac disease and circulating gluten-specific T cells can be identified by staining with HLA-DQ:gluten tetramers. In this first single-cell RNA-seq study of tetramer-sorted T cells from untreated celiac disease patients blood, we found that gluten-specific T cells showed distinct transcriptomic profiles consistent with activated effector memory T cells that shared features with Th1 and follicular helper T cells. Compared to non-specific cells, gluten-specific T cells showed differential expression of several genes involved in T-cell receptor signaling, translational processes, apoptosis, fatty acid transport, and redox potentials. Many of the gluten-specific T cells studied shared T-cell receptor with each other, indicating that circulating gluten-specific T cells belong to a limited number of clones. Moreover, the transcriptional profiles of cells that shared the same clonal origin were transcriptionally more similar compared with between clonally unrelated gluten-specific cells.

  Study EGAS00001005517

• The clinicopathologic spectrum and genomic landscape of de-/trans-differentiated melanoma

  De- and trans-differentiation is a rare and only poorly understood phenomenon in cutaneous melanoma. To study this disease more comprehensively we have retrieved 11 primary cutaneous melanomas from our pathology archives showing biphasic features characterized by a conventional melanoma and additional areas of de-/trans-differentiation as defined by a lack of immunohistochemical expression of all conventional melanocytic markers (S-100 protein, SOX10, Melan-A and HMB-45). The clinical, histologic and immunohistochemical findings were recorded and follow-up was obtained. The patients were mostly elderly (median: 81 years; range: 42-86 years) without significant gender predilection, and the sun-exposed skin of the head and neck area was most commonly affected. The tumors were deeply invasive with a mean tumor thickness of 7 mm (range: 4-80 mm). The dedifferentiated component showed atypical fibroxanthoma-like features in the majority (7), while additional rhabdomyosarcomatous and epithelial transdifferentiation was noted histologically and/or immunohistochemically in two tumors each. The background conventional melanoma component was of desmoplastic (4), superficial spreading (3), nodular (2), lentigo maligna (1) or spindle cell (1) types. For the 7 patients with available follow-up data (median follow-up period of 25 months; range: 8-36 months), 2 died from their disease and 3 developed metastases. Next-generation sequencing of the cohort revealed somatic mutation of established melanoma drivers including mainly NF1 mutations in the conventional component (5 cases), which were also detected in the corresponding de-/trans-differentiated components. In summary, the diagnosis of de-/trans-differentiated melanoma is challenging and depends on the morphologic identification of the conventional melanoma component. Molecular analysis is diagnostically helpful as the mutated gene profile is shared between the conventional and de-/trans-differentiated components. Importantly, de-/trans-differentiation does not appear to confer a more aggressive behavior.

  Dataset EGAD00001007034

• The clinicopathologic spectrum and genomic landscape of de-/trans-differentiated melanoma

  De- and trans-differentiation is a rare and only poorly understood phenomenon in cutaneous melanoma. To study this disease more comprehensively we have retrieved 11 primary cutaneous melanomas from our pathology archives showing biphasic features characterized by a conventional melanoma and additional areas of de-/trans-differentiation as defined by a lack of immunohistochemical expression of all conventional melanocytic markers (S-100 protein, SOX10, Melan-A and HMB-45). The clinical, histologic and immunohistochemical findings were recorded and follow-up was obtained. The patients were mostly elderly (median: 81 years; range: 42-86 years) without significant gender predilection, and the sun-exposed skin of the head and neck area was most commonly affected. The tumors were deeply invasive with a mean tumor thickness of 7 mm (range: 4-80 mm). The dedifferentiated component showed atypical fibroxanthoma-like features in the majority (7), while additional rhabdomyosarcomatous and epithelial transdifferentiation was noted histologically and/or immunohistochemically in two tumors each. The background conventional melanoma component was of desmoplastic (4), superficial spreading (3), nodular (2), lentigo maligna (1) or spindle cell (1) types. For the 7 patients with available follow-up data (median follow-up period of 25 months; range: 8-36 months), 2 died from their disease and 3 developed metastases. Next-generation sequencing of the cohort revealed somatic mutation of established melanoma drivers including mainly NF1 mutations in the conventional component (5 cases), which were also detected in the corresponding de-/trans-differentiated components. In summary, the diagnosis of de-/trans-differentiated melanoma is challenging and depends on the morphologic identification of the conventional melanoma component. Molecular analysis is diagnostically helpful as the mutated gene profile is shared between the conventional and de-/trans-differentiated components. Importantly, de-/trans-differentiation does not appear to confer a more aggressive behavior.

  Dataset EGAD00001007033

• FASTQ files of total RNA-Seq data from the POPS SGA (Small for Gestational Age) samples

  All the samples were obtained from the Pregnancy Outcome Prediction–a prospective cohort study of nulliparous women attending the Rosie Hospital, Cambridge (UK) for their dating ultrasound scan between January 14, 2008, and July 31, 2012. Ethical approval for the study was given by the Cambridgeshire 2 Research Ethics Committee (reference number 07/H0308/163) and all participants provided written informed consent. Cases of preeclampsia (PET) were defined on the basis of the 2013 ACOG criteria and cases of small for gestational age (SGA)infants were confined to severe SGA, i.e. a customized birth weight <5th percentile. Chorionic villi from the corresponding placentas (free from decidua, visible infarction, calcification, hematoma, or damage) were collected and processed within 30 minutes of separation from the uterus. After repeated washes in chilled phosphate buffered saline, the samples were placed in RNA later (Applied Biosystems) and stored at -80°C. Total placental RNA was extracted using mirVana Isolation Kit (Ambion). For each placenta, approximately 5 mg of tissue were homogenized in the Lysis/Binding solution for 20 sec at 6 m/s using a bead beater (FastPrep24) and Lysing Matrix D Tubes (MP Biomedicals). The samples were then spun at 13,000 rpm for 5 min at 4°C and the supernatants recovered. Afterwards, the manufacturer's instructions were followed. Immediately after the RNA extraction, placental RNA samples were DNase-treated using DNA-free DNA Removal Kit (Ambion), aliquoted, and stored in -80°C. Quantity and quality of the RNA samples were assessed using the Agilent 2100 Bioanalyzer, the Agilent RNA 6000 Nano Kit (Agilent Technologies), and Qubit fluorometer. Libraries were prepared starting with 300-500 ng of good quality total RNA (RIN ≥7.5) using the TruSeq Stranded Total RNA Library Prep Kit with Ribo-Zero Human/Mouse/Rat (Illumina), according to the manufacturer's instructions. The kit contains 96 uniquely indexed adapter combinations in order to allow pooling of multiple samples prior to sequencing. After determining their size (with the Agilent 2100 Bioanalyzer and the Agilent High Sensitivity DNA Kit by Agilent Technologies) and concentration (by qPCR with the KAPA Illumina ABI Prism Library Quantification Kit, Kapa Biosystems), libraries have been pooled and sequenced (single-end, 125 bp) using a Single End V4 Cluster Kit and an Illumina HiSeq2500 or HiSeq4000 instrument.

  Dataset EGAD00001003507

• FASTQ files of total RNA-Seq data from the POPS PET (pre-eclamptic) samples

  All the samples were obtained from the Pregnancy Outcome Prediction–a prospective cohort study of nulliparous women attending the Rosie Hospital, Cambridge (UK) for their dating ultrasound scan between January 14, 2008, and July 31, 2012. Ethical approval for the study was given by the Cambridgeshire 2 Research Ethics Committee (reference number 07/H0308/163) and all participants provided written informed consent. Cases of preeclampsia (PET) were defined on the basis of the 2013 ACOG criteria and cases of small for gestational age (SGA)infants were confined to severe SGA, i.e. a customized birth weight <5th percentile. Chorionic villi from the corresponding placentas (free from decidua, visible infarction, calcification, hematoma, or damage) were collected and processed within 30 minutes of separation from the uterus. After repeated washes in chilled phosphate buffered saline, the samples were placed in RNA later (Applied Biosystems) and stored at -80°C. Total placental RNA was extracted using mirVana Isolation Kit (Ambion). For each placenta, approximately 5 mg of tissue were homogenized in the Lysis/Binding solution for 20 sec at 6 m/s using a bead beater (FastPrep24) and Lysing Matrix D Tubes (MP Biomedicals). The samples were then spun at 13,000 rpm for 5 min at 4°C and the supernatants recovered. Afterwards, the manufacturer's instructions were followed. Immediately after the RNA extraction, placental RNA samples were DNase-treated using DNA-free DNA Removal Kit (Ambion), aliquoted, and stored in -80°C. Quantity and quality of the RNA samples were assessed using the Agilent 2100 Bioanalyzer, the Agilent RNA 6000 Nano Kit (Agilent Technologies), and Qubit fluorometer. Libraries were prepared starting with 300-500 ng of good quality total RNA (RIN ≥7.5) using the TruSeq Stranded Total RNA Library Prep Kit with Ribo-Zero Human/Mouse/Rat (Illumina), according to the manufacturer's instructions. The kit contains 96 uniquely indexed adapter combinations in order to allow pooling of multiple samples prior to sequencing. After determining their size (with the Agilent 2100 Bioanalyzer and the Agilent High Sensitivity DNA Kit by Agilent Technologies) and concentration (by qPCR with the KAPA Illumina ABI Prism Library Quantification Kit, Kapa Biosystems), libraries have been pooled and sequenced (single-end, 125 bp) using a Single End V4 Cluster Kit and an Illumina HiSeq2500 or HiSeq4000 instrument.

  Dataset EGAD00001003508

• Institut Curie Neuroblastoma Whole Genome Sequencing Diagnosis Relapse

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

  Study EGAS00001001184

• The genomic complexity of early T-cell progenitor acute lymphoblastic leukemia

  CREST The accurate identification of structural variations using whole-genome DNA sequencing data generated by next-generation sequencing technology is extremely difficult. To address this challenge, we have developed CREST, an algorithm that uses sequencing reads with partial alignments to the reference human genome (so-called soft-clipped reads) to directly map the breakpoints of somatic structural variations. We applied CREST to paired tumor/normal whole genome sequencing data from five cases of T-lineage acute lymphoblastic leukemia (T-ALL). A total of 110 somatic structural variants were identified, >80% of which were validated by genomic PCR and Sanger sequencing. The validated structural variants included 31 inter-chromosomal translocations, 19 intra-chromosomal translocations, one inversion, 22 deletions and 16 insertions. A comparison of the results generated with CREST to those obtained using the traditional paired-end discordant mapping methods demonstrate CREST to have a much higher sensitivity and specificity. In addition, application of CREST to publicly available whole-genome sequencing data from the human melanoma cancer cell line COLO-829 demonstrated the identification of 50 novel structural variations not detected using the standard methods, 20 of which were selected for validation with a 90% success rate. These data demonstrate that direct mapping of soft-clipped reads offers an improved method for detecting structural variants at the nucleotide level of resolution. T-ALL Early T-cell precursor acute lymphoblastic leukaemia (ETP ALL) is an aggressive malignancy of unknown genetic basis. We performed whole-genome sequencing of 12 ETP ALL cases and assessed the frequency of the identified somatic mutations in 94 T-cell acute lymphoblastic leukaemia cases. ETP ALL was characterized by activating mutations in genes regulating cytokine receptor and RAS signalling (67% of cases; NRAS, KRAS, FLT3, IL7R, JAK3, JAK1, SH2B3 and BRAF), inactivating lesions disrupting haematopoietic development (58%; GATA3, ETV6, RUNX1, IKZF1 and EP300) and histone-modifying genes (48%; EZH2, EED, SUZ12, SETD2 and EP300). We also identified new targets of recurrent mutation including DNM2, ECT2L and RELN. The mutational spectrum is similar to myeloid tumours, and moreover, the global transcriptional profile of ETP ALL was similar to that of normal and myeloid leukaemia haematopoietic stem cells. These findings suggest that addition of myeloid-directed therapies might improve the poor outcome of ETP ALL.

  Study phs000340

• Circulating Genomic Determinants of Treatment Failure in Hodgkin Lymphoma

  In this study, we show that the plasma representation of mutations exceeds the bulk tumor representation in most cases, making classic Hodgkin Lymphoma (cHL) particularly amenable to noninvasive profiling. Leveraging single-cell transcriptional profiles of cHL tumors, we demonstrate hazard ratios (HRs) circulating tumor DNA (ctDNA) shedding to be shaped by DNASE1L3, whose increased tumor microenvironment-derived expression drives high ctDNA concentrations. Using this insight, we comprehensively profile 366 patients, revealing two distinct cHL genomic subtypes with characteristic clinical and prognostic correlates, as well as distinct transcriptional and immunological profiles. Furthermore, we identify a novel class of truncating IL4R-mutations that are dependent on IL13 signaling and therapeutically targetable with IL4R blocking antibodies. Finally, we demonstrate the clinical value of pre- and on-treatment ctDNA levels for longitudinally refining cHL risk prediction, and for detection of radiographically occult minimal residual disease. Collectively, these results support the utility of noninvasive strategies for genotyping and dynamic monitoring of cHL as well as capturing molecularly distinct subtypes with diagnostic, prognostic, and therapeutic potential.Methods utilized include a) targeted genotyping of baseline plasma samples, b) whole exome sequencing (WES) of baseline plasma samples, c) epigenetic expression inference from cell-free DNA-sequencing (EPIC-Seq) to infer gene expression from cell-free DNA (cfDNA) fragmentation profiles, d) ctDNA tracking in on-treatment samples, e) RNA-Sequencing, and f) Phased variant Enrichment and Detection Sequencing (PhasED-Seq).All patients were treated at cancer centers across Europe and North America between 2011 and 2020. The following 3 cancer centers, Stanford (CA, USA), Bellinzona [Confoederatio Helvetica (CH) or Switzerland], and Leuven (Belgium), were included. The 6 clinical trial collectives included in this study are the Phase III trial AHOD1331 trial (NCT02166463) and the Phase II cHOD17 study (NCT03755804) for enrolled pediatric patients, the Phase III AHL2011 (NCT01358747), the Phase II BREACH trial (NCT02292979), a pilot study evaluating concurrent Pembrolizumab with AVD (Doxorubicin, also known as Adriamycin, Vinblastine, and Dacarbazine; NCT03331341), and a Phase II study evaluating PVAB (Prednisone, Vinblastine, Doxorubicin and Bendamustine; NCT02414568) in older patients (61 years and older). In addition, lymph node specimens from a biorepository at St. Jude Children's Research Hospital (TN, USA) were included.

  Study phs003435

• Clinical Proteomic Tumor Analysis Consortium (CPTAC) Proteogenomic Confirmatory Study of Breast, Colon, Lung, and Ovarian Tumors

  Recently, significant progress has been made in characterizing and sequencing the genomic alterations in statistically robust numbers of samples from several types of cancer. For example, The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC) and other similar efforts are identifying genomic alterations associated with specific cancers (e.g., copy number aberrations, rearrangements, point mutations, epigenomic changes, etc.). The availability of these multi-dimensional data to the scientific community sets the stage for the development of new molecularly targeted cancer interventions. Understanding the comprehensive functional changes in cancer proteomes arising from genomic alterations and other factors is the next logical step in the development of high-value candidate protein biomarkers. Hence, proteomics can greatly advance the understanding of molecular mechanisms of disease pathology via the analysis of changes in protein expression, their modifications and variations, as well as protein-protein interaction, signaling pathways and networks responsible for cellular functions such as apoptosis and oncogenesis. Realizing this great potential, the NCI launched the second phase of the CPTC initiative in September 2011. Renamed the Clinical Proteomic Tumor Analysis Consortium, CPTAC is beginning to leverage its analytical outputs from Phase I to define cancer proteomes on genomically-characterized biospecimens. The purpose of this integrative approach is to provide the broad scientific community with knowledge that links genotype to proteotype and ultimately phenotype. The data contained in this dataset are derived from samples designed to confirm CPTAC findings from the TCGA samples. These confirmatory samples contain breast, ovarian, colon, and lung tumors collected via a protocol optimized for proteomics. Specifically, ischemic time of the sample was controlled and restricted to less than 30 minutes. ACGT, Inc. produced whole exome, mRNAseq, and miRNAseq for these samples. Corresponding proteomic data are available at: https://cptac-data-portal.georgetown.edu/cptacPublic/ The study design was to profile colon, breast, ovarian, and lung tumors both genomically and proteomically. Germline DNA was obtained from blood. Normal control samples for proteomics varied by organ site: adjacent colon tissue for colon cases, contralateral breast tissue for some breast cases, and Fallopian tube fimbria for some ovarian cases. Lung cases had no normal control for proteomic analysis. All cancer samples were derived from primary and untreated tumors.

  Study phs000892

• Institut Curie Neuroblastoma Whole Genome Sequencing Diagnosis Relapse

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

  Dataset EGAD00001001356

• A Clinical Trial of Pembrolizumab in Patients with Hepatitis B Virus-related Hepatocellular Carcinoma, with Parallel Study on Baseline and Serial Change in the Immune Environment

  This is a single-arm Phase II trial of pembrolizumab in patients with hepatitis B virus-related hepatocellular carcinoma with parallel study on baseline and serial change in the immune environment. Subjects should have a confirmed diagnosis of HCC (in accordance with the AASLD guideline) and confirmed chronic infection with hepatitis B virus as defined by positivity for HBsAg. Antiviral therapy for HBV must be given for at least 12 weeks and HBV viral load must be less than 100 IU/mL prior to first dose of study drug. They must have disease not amenable to a curative treatment approach or loco-ablation. Subject must be fit and agreeable with baseline and post-treatment biopsy of tumor. Subjects must have at least one measurable lesion per Response Evaluation Criteria in Solid Tumors (RECIST) 1.1, an Eastern Cooperative Oncology Group (ECOG) performance score of 0 or 1 and adequate organ functions. 30 subjects will be enrolled to receive pembrolizumab 200 mg IV every 3 weeks(Q3W). Pre-treatment and on-treatment biopsy after 2 cycles of Pembrolizumab will be preformed. Treatment will be stopped when progression of disease or intolerable toxicity occurs. The primary objectives of this trial are to study the efficacy and safety of pembrolizumab in patients with HBV-related HCC and to study the serial change in RNA expression of immune-related gene panel in post-treatment biopsy tissue. The secondary objectives of this trial are to study the serial change in cytokine profile between pre-treatment and post-treatment samples, to study the PD-L1 immunohistochemical (IHC) expression in tumor sample at baseline and post-treatment tissue samples and to study the presence of tumor infiltrating lymphocytes in the baseline and post-treatment tumor samples. The exploratory objective of this trial is to evaluate the possibility of using baseline and the serial change in RNA expression of immune-related gene panel or PD-L1/2 IHC to predict treatment response.

  Dataset EGAD50000001872

• The cellular immune and airway epithelial profile throughout childhood and in response to COVID-19 at multi-omic single cell level.

  Single-cell data gene expression data set (5’Chromium 10X) of healthy paediatric volunteers, and paediatric and adult COVID-19 patients. Gene expression was determined from samples of nasal, tracheal and bronchial brushings and blood (PBMCs). In addition to gene expression, PBMC’s were assayed by CITE-seq. A subset of samples have VDJ sequencing data for T cell receptors (TCR) and B cell receptors (BCR).

  Dataset EGAD00001007718

• RNA-sequencing of human skin samples obtained from SSc patients and healthy controls.

  Patient information SSc patients were recruited at the Department of Rheumatology of the Leiden University Medical Center (Leiden, The Netherlands). All patients met the American Rheumatism Association classification criteria for SSc (Subcommittee for scleroderma criteria 1980), and were classified according to LeRoy and Medsger criteria as either limited or diffuse cutaneous disease (LeRoy EC, Black C, Fleischmajer R, Jablonska S, Krieg T, Medsger TA Jr, Rowell N 1988). Institutional review board approval and written informed consent was obtained before patients entered this study. Two 4 mm skin biopsies were taken from a standardized location on the most proximal part of the lower arm, distal from the elbow. In 10 patients the skin biopsy came from a clinically affected area and in 4 patients the skin was locally unaffected. One sample was used for RNA sequencing and one sample was used for immunohistochemistry. Skin biopsies from healthy individuals were commercially sourced (Tissue Solutions, UK) and collected from donors undergoing skin resection surgery and after informed consent. To match the healthy skin with patients as much as possible, skin biopsies from healthy controls were also taken from a similar position (the under-arm (for 4 controls) and leg (for 2 controls)). Healthy skin donors were selected to match the age and sex of the SSc patient cohort. Biopsies from patients and controls were equally treated and were both stored at -80°C until RNA isolation was performed. RNA from frozen skin biopsies was isolated using RNeasy kit from fibrous tissue (Qiagen, the Netherlands). RNA quantity was determined by using SimplyNano 2000 and quality was assessed on Tapestation (Agilent, the Netherlands). All samples included in the study had a RIN score above 7.0. Transcriptome characterisation and analysis RNA sequencing was performed using polyA selection and a stranded protocol using Ion Torrent next generation sequencing technology (Service XS, The Netherlands). The Ion PI Template OT2 200 Kit v3 and Ion PI Sequencing 200 Kit v3 were used according to the manufacturer’s instructions. 20 samples were run on 11 PI chips. PI chip analyses, base calling and quality checks were performed using the Torrent Server Suite. An average of 42 million 100 bp reads was generated per sample. Following quality control, reads were aligned to the human genome (Homo sapiens GRh38.78) using Bowtie2 and STAR (Dobin et al. 2013; Langmead and Salzberg 2012). Reads were first aligned with STAR. For the unmapped reads from STAR, a second alignment step was performed using bowtie2 (local very sensitive options)

  Dataset EGAD00001003832

• Environmental Influences on Child Health Outcomes (ECHO) - PATHWAYS - Conditions Affecting Neurocognitive Development and Learning in Early Childhood (CANDLE)

  Pregnant women are exposed daily to multiple chemical and non-chemical stressors, including air pollutants, phthalates, and psychosocial stress. The Developmental Origins of Health and Disease model (DOHaD) states that exposure to these stressors in pregnancy affects fetal development in a manner that impacts offspring health across the life course. Yet epidemiologic data in these areas is limited, particularly with U.S. samples. It is also poorly understood whether these exposures prenatally affect childhood neurodevelopment and airway health in an independent or combined manner, and likely moderators of effects, such as the sex of the child, are infrequently addressed. We propose to unify three diverse extant pregnancy cohorts: TIDES (N=717, 2010-12), GAPPS (N=1133, 2012-16), and CANDLE (N=1385, 2007-11), for a combined sample size of 3235 mother-child dyads in the PATHWAYS study. PATHWAYS will investigate how chemical (air pollutants and phthalates) and non-chemical (psychosocial stress) exposures during pregnancy are related to placental gene expression (transcriptome) and childhood neurodevelopment and airway health (at ages 4-6, 8-9, and 10-11 years). Each cohort has rich resources of prenatal data and banked specimens (urine, blood, and placenta) that will be harmonized for the PATHWAYS study and will contribute to the ECHO consortium. We will develop a national model with high spatiotemporal resolution of key air pollutants and assess urinary markers of maternal exposure. A composite measure will capture multilevel maternal psychosocial stress across pregnancy, and urinary phthalate and blood stress hormone [Corticotropin-releasing hormone (CRH)] levels in the second and third trimesters will provide individual assessment of those exposures in potential critical periods. We will characterize the placental transcriptome using RNA sequencing (RNA-Seq) and will assess neurodevelopment and airway health prospectively into middle childhood. PATHWAYS will examine how these prenatal exposures are related to the placental transcriptome and child health outcomes in main effect and interactive models, with emphasis on sex-specific associations. For both neurodevelopment and airway health, we propose to measure both phenotypic precursors of health outcomes (i.e. fluid cognition, lung function growth), which yield dimensional tests of proposed associations, as well as clinically meaningful and policy relevant outcomes (i.e. asthma, mental health). Our study is powered to assess interactive effects of chemical and non-chemical stressors and will be the first study to characterize how prenatal environmental exposures relate to placental transcriptome pathways in relation to childhood health outcomes. This represents a significant scientific advance in testing DOHaD hypotheses. Major contributions to the ECHO Consortium include: 1) the development of a state of the art national model of air pollution, 2) a large, diverse pregnancy cohort with extensive biorepositories and extant prenatal and postnatal biomarkers, placental transcriptome, psychosocial and environmental data, and 3) an experienced, interdisciplinary team that will contribute meaningfully to the ECHO program of work.The sequencing data in this submission pertains only to the 794 subjects from the CANDLE study for whom placental tissue transcriptomics was performed.

  Study phs003619

• Environmental Influences on Child Health Outcomes (ECHO) - PATHWAYS - Global Alliance to Prevent Prematurity and Stillbirth (GAPPS)

  Pregnant women are exposed daily to multiple chemical and non-chemical stressors, including air pollutants, phthalates, and psychosocial stress. The Developmental Origins of Health and Disease model (DOHaD) states that exposure to these stressors in pregnancy affects fetal development in a manner that impacts offspring health across the life course. Yet epidemiologic data in these areas is limited, particularly with U.S. samples. It is also poorly understood whether these exposures prenatally affect childhood neurodevelopment and airway health in an independent or combined manner, and likely moderators of effects, such as the sex of the child, are infrequently addressed. We propose to unify three diverse extant pregnancy cohorts: TIDES (N=717, 2010-12), GAPPS (N=1133, 2012-16), and CANDLE (N=1385, 2007-11), for a combined sample size of 3235 mother-child dyads in the PATHWAYS study. PATHWAYS will investigate how chemical (air pollutants and phthalates) and non-chemical (psychosocial stress) exposures during pregnancy are related to placental gene expression (transcriptome) and childhood neurodevelopment and airway health (at ages 4-6, 8-9, and 10-11 years). Each cohort has rich resources of prenatal data and banked specimens (urine, blood, and placenta) that will be harmonized for the PATHWAYS study and will contribute to the ECHO consortium. We will develop a national model with high spatiotemporal resolution of key air pollutants and assess urinary markers of maternal exposure. A composite measure will capture multilevel maternal psychosocial stress across pregnancy, and urinary phthalate and blood stress hormone [Corticotropin-Releasing Hormone (CRH)] levels in the second and third trimesters will provide individual assessment of those exposures in potential critical periods. We will characterize the placental transcriptome using RNA sequencing (RNA-Seq) and will assess neurodevelopment and airway health prospectively into middle childhood. PATHWAYS will examine how these prenatal exposures are related to the placental transcriptome and child health outcomes in main effect and interactive models, with emphasis on sex-specific associations. For both neurodevelopment and airway health, we propose to measure both phenotypic precursors of health outcomes (i.e., fluid cognition, lung function growth), which yield dimensional tests of proposed associations, as well as clinically meaningful and policy relevant outcomes (i.e., asthma, mental health). Our study is powered to assess interactive effects of chemical and non-chemical stressors and will be the first study to characterize how prenatal environmental exposures relate to placental transcriptome pathways in relation to childhood health outcomes. This represents a significant scientific advance in testing DOHaD hypotheses. Major contributions to the ECHO Consortium include: 1) the development of a state of the art national model of air pollution, 2) a large, diverse pregnancy cohort with extensive biorepositories and extant prenatal and postnatal biomarkers, placental transcriptome, psychosocial and environmental data, and 3) an experienced, interdisciplinary team that will contribute meaningfully to the ECHO program of work.The sequencing data in this submission pertain only to the 408 subjects from the GAPPS study for whom placental tissue transcriptomics was performed.

  Study phs003620

• DNA methylomes of monozygotic twins clinically discordant for multiple sclerosis

  Multiple sclerosis (MS) is an inflammatory demyelinating disease of the central nervous system. Although genetic susceptibility is important, a modest concordance rate for MS in monozygotic (MZ) twins suggests that interaction with other risk factors is required to develop clinical symptoms. In this study, we examined whether DNA methylation differences contribute to the discordant clinical manifestation of MS in MZ twins, and studied the impact of MS treatments on the DNA methylome. Genome-wide DNA methylation profiles of peripheral blood mononuclear cells (PBMCs) of 45 MZ twins clinically discordant for MS were generated using the Illumina HumanMethylationEPIC array. Repetitive element methylation and selected differentially methylated positions (DMPs) were validated using targeted deep bisulfite sequencing (TDBS). In addition, we performed whole genome bisulfite sequencing (WGBS) to profile CD4+ memory T-cells of a subset of four MS discordant MZ twins. Our results show that interferon treatment causes robust DNA methylation changes and several epigenetic biomarkers for interferon treatment response were identified. However, overall the PBMC-based methylomes of the MS discordant MZ twins were highly similar, since large systematic methylation differences (>5%) were absent in the data. This suggests that previously reported large methylation changes are probably caused by genetic rather than epigenetic differences. In addition, our data does not support the hypothesis that the observed maternal parent-of-origin effect in MS is due to genomic imprinting errors, and no evidence was found that copy number variations explain the discordant phenotype in these MZ twins. Although not genome-wide significant, a couple of DMPs associated with the MS phenotype were identified and successfully technically replicated using TDBS, including a differentially methylated region (DMR) in the promoter of the transmembrane protein encoding gene TMEM232. Another MS-associated DMP, located in an enhancer in the gene body of the transcription factor ZBTB16, was also associated with medium-term glucocorticoid treatment history. WGBS analysis confirmed this DMP as a promising epigenetic biomarker for glucocorticoid treatment. In conclusion, this epigenome-wide association study (EWAS) in genetically identical twins identified a DMR in the TMEM232 promoter as a candidate loci associated with the clinical manifestation of MS. In addition, epigenetic biomarkers for MS treatments were identified, revealing that not only short-term, but also medium-term treatment effects are detectable in immune cells, which should be taken into account in future EWAS designs.

  Study EGAS00001003147

• The subclonal architecture of metastatic breast cancer: Results from a prospective community-based rapid autopsy program 'CASCADE'

  Background: Understanding the cancer genome is seen as a key step in improving outcomes for cancer patients. Genomic assays are emerging as a possible avenue to personalised medicine in breast cancer. The majority of work in this area has targeted primary tumours however, and very few studies have performed comprehensive profiling of advanced disease. Evolution of the cancer genome during the natural history of breast cancer is largely unknown, as is the profile of disease at death. We sought to study in detail these aspects of advanced breast cancers that have resulted in lethal disease. Methods and Findings: Three patients with ER-positive, HER2-negative breast cancer and one patient with triple negative breast cancer underwent rapid autopsy as part of an institutional prospective community-based rapid autopsy program. Cases represented a range of management problems in breast cancer, including late relapse after early stage disease; de novo metastatic disease; discordant disease response and disease refractory to treatment. Between 5 and 12 metastatic sites were collected at autopsy together with available primary tumours and longitudinal metastatic biopsies taken during life. Samples underwent paired tumour-normal whole exome sequencing and single nucleotide polymorphism arrays. Subclonal architectures were inferred by jointly analysing all samples from each patient. Mutations were validated using high depth amplicon sequencing.Between cases, there were significant differences in mutational burden, driver mutations, mutational processes and copy number variation. Within each case, we found dramatic heterogeneity in subclonal structure from primary to metastatic disease and between metastatic sites, such that no single lesion captured the breadth of disease. Metastatic cross seeding was found in each case and treatment drove subclonal diversification. Subclones displayed parallel evolution of treatment resistance in some cases, and apparent augmentation of key oncogenic drivers as an alternative resistance mechanism. We also observed the key role of mutational processes in subclonal evolution.Limitations of this study include the potential for bias introduced by joint analysis of formalin fixed archival specimens with fresh specimens, and the difficulties in resolving subclones with whole exome sequencing. Other alterations that could define subclones such as structural variants or epigenetic modifications were not assessed. Conclusions: This study highlights the variety of mechanisms that shape the genome of metastatic breast cancer, and the value of studying advanced disease in detail. Treatment drives significant genomic heterogeneity in breast cancers which has implications for disease monitoring and treatment selection in the personalised medicine paradigm.

  Study EGAS00001002153

• Environmental Arsenic and Diabetes Mellitus (Chihuahua Cohort)

  The Chihuahua Cohort was established as part of a larger NIEHS-funded project that used laboratory models (tissue culture and mice) and human subjects to examine the association between exposure to arsenic and diabetes, including identification of mechanisms underlying this association. All procedures involving human subjects were approved by institutional review boards at the University of North Carolina at Chapel Hill and Cinvestav-IPN (Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Mexico City, Mexico). All participants provided signed informed consent. A total of 1,160 adults (≥ 18 years old) with a minimum 5-year uninterrupted residency in the Chihuahua study area were recruited between 2008 and 2012. Samples of drinking water were obtained from the participants' households. Data on residency, occupation, drinking-water sources and use, smoking, alcohol consumption, and medical history were recorded by a study questionnaire. Spot urine and fasting venous blood were collected during medical examination that included an oral glucose tolerance test (OGTT) with blood drawn 2 hr after a 75-g glucose dose. Plasma from both fasting and 2-hr blood samples was collected. Body mass index (BMI) of each participant was calculated from body weight and height. BMI cutoffs of ≥ 25.0, ≥ 30, and < 18.5 kg/m2 were used to define overweight, obese, and underweight individuals, respectively. Participants' waist circumference was also measured. Blood pressure was assessed using a manual sphygmomanometer. Three measurements were taken at intervals of at least 1 min, with a 5-min rest before obtaining the first reading; the mean of the last two measurements was used. The concentrations of arsenic in drinking water and arsenic species in urine were measured by hydride generation-atomic absorption spectrometry. The limit of detection for As in water as well as for As species in urine was 0.01 µg As/L. Creatinine concentration in urine was determined using an established colorimetric assay. DNA was isolated from blood for genotyping analysis.The following measures were taken: glucose concentration in fasting blood (FBG) and in blood collected 2 hours into OGTT (2HBG); concentrations of triglycerides, HDL and total cholesterol (TC) in fasting plasma; and LDL was calculated using the Friedewald equation. Diabetes was classified by FPG ≥ 126 mg/dL, 2HPG ≥ 200 mg/dL, or self-reported diabetes diagnosis or medication use (WHO/International Diabetes Federation 2006). Prediabetes was defined as the absence of diabetes with FPG ≥ 110 mg/ dL or 2HPG ≥ 140 mg/dL. Elevated fasting levels of each lipid were defined as plasma TG ≥ 150 mg/dL, TC ≥ 200 mg/dL, and LDL ≥ 130 mg/dL. Fasting HDL < 40 mg/dL in men and < 50 mg/dL in women were designated as low. Hypertension was defined by systolic blood pressure (SBP) > 140 mm Hg, diastolic blood pressure (DBP) > 90 mmHg, or self-reported use of anti-hypertensive medication. Genotyping analysis, which will be available in version 2 of this study, was carried out in a subset of ~600 subjects with focus on enzymes involved in the metabolism of arsenic, including AS3MT, GSTs.All data collected during the medical exams or generated by the above-described chemical and biochemical analysis of water, urine and blood/plasma, as well as results of the genotyping analysis are included in dbGaP. Data generated by metabolomics and epigenetic analyses in small subsets of the Chihuahua cohort are not included, but are available in published articles.

  Study phs002139

• Disease Variant Landscape of a Large Multiethnic Population of Moyamoya Patients by Exome Sequencing

  We report a significant advance in Moyamoya disease (MMD) research through our unique access to a large MMD patient population combined with high depth exome sequencing and bioinformatics. We confirm a reported RNF213 founder mutation (FM), provide evidence for its ethnic specificity and identify novel variants in other genes that associate with Caucasian and non-RNF213 FM MMD. This work has a broader impact in vascular research by highlighting novel genetic aspects of cerebrovascular disease. Moyamoya Disease (MMD) is a rare cerebrovascular disorder characterized by unilateral or bilateral progressive stenosis or occlusion of the internal carotid artery, with frequent involvement of the anterior cerebral artery and the middle cerebral artery (MCA). Fragile collateral vessels subsequently develop and are particularly prone to hemorrhage. MMD patients are usually diagnosed angiographically following presentation of developmental delays, seizures, migraines, ischemia, and/or hemorrhage. The most common approach to re-establish normal blood flow is bypass surgery. While some progress has been made in identifying disease-associated genes, the underlying disease biology is not well understood. Moyamoya disease (MMD) is a rare cerebrovascular disease. Our unique access to a large MMD patient population combined with high depth exome sequencing and bioinformatics has led to a significant advance in the field. We confirm a reported RNF213 founder mutation (FM), provide evidence for its ethnic specificity and identify novel variants in other genes that associate with Caucasian and non-RNF213 FM MMD. This work has a broader impact in vascular research by highlighting novel genetic aspects of cerebrovascular disease. The 125 ethnically diverse, unrelated MMD patients were matched based on sex and broad ethnic category to 125 controls. Control DNA from the 1000 Genomes project was used to design a control Library by Personalis. 125 controls were selected from this Personalis Control DNA Library for our study. (The 1000 Genomes Project Consortium 2012). Genomic DNA was extracted using the Gentra Puregene kit (Qiagen, Valencia, CA). Libraries were prepared from approximately 3 µg of high quality genomic DNA (50-200 ng/µl) using Illumina TruSeq Genomic DNA High throughput Sample Prep Kits (Illumina, San Diego, CA) and exome enrichment (targeting 62Mb) was accomplished using the TruSeq Exome Target Enrichment kit (Illumina, San Diego, CA), all according to manufacturer's protocols. Target enrichment validation was confirmed by determining the concentration of the library by PicoGreen-based quantitation. Library yields ranged from 100-1000ng of DNA, a portion of which was run on the Bioanalyzer HS DNA chip (Agilent, Santa Clara, CA), with an average size of 300-550nt for DNA fragments. Exome Sequencing. Sequencing was performed using Illumina Hiseq2000 or HiSeq2500 sequencers with single lane, paired-end 2X100bp reads. DNA fragments were generated and amplified using Clonal Single Molecule Array technology (Illumina, San Diego, CA). The sequences were determined using the Clonal Single Molecule Array and Sequencing-by-Synthesis using Illumina's instrumentation and Reversible Terminator Chemistry. Each sequencing lane interrogated the DNA sequences of a pool of 3 individual sample libraries each carrying a unique index. Sequencing reads of at least 2x100bp in length for a total of approximately 8 Gb of sequence data were generated for each sample.

  Study phs001700

• Metagenomic Analysis of the Structure and Function of the Human Gut Microbiota in Crohn's Disease

  Inflammatory bowel diseases (IBD), such as Crohn's disease, are chronic, immunologically mediated disorders that have severe medical consequences. The current hypothesis is that these diseases are due to an overly aggressive immune response to a subset of commensal enteric bacteria. Studies to date on IBD have suggested that the disorder may be caused by a combination of bacteria and host susceptibility; however the etiologies of these diseases remain an enigma. In this application, we propose to develop and demonstrate the ability to profile Crohn's disease at an unprecedented molecular level by elucidation of specific biomarkers (bacterial strains, genes, or proteins) that correlate to disease symptoms. To achieve this goal, we will employ a multidisciplinary approach based on metagenomic and metaproteomic molecular tools to elucidate the composition of the commensal microbiota in monozygotic twins that are either healthy or exhibit Crohn's disease (for concordant, both are diseased; for discordant, one is healthy and one is diseased). The central hypotheses of this proposal are (1) that specific members and/or functional activities of the gastrointestinal (GI) microbiota differ in patients with Crohn's disease as compared to healthy individuals, and (2) that it will be possible to elucidate microbial signatures which correlate with the occurrence and progression of this disease by integration of data obtained from 16S rRNA-based molecular fingerprinting, metagenomics, and metaproteomics approaches. To address these hypotheses, three specific aims are proposed: 1) Obtain data on community gene content (metagenome) in a subset of healthy twins and twins with Crohn's Disease to assess potential differences in the metabolic capabilities of the gut microbiota associated with CD, 2) Obtain data on community protein content (metaproteome) in a subset of healthy twins and twins with Crohn's Disease to assess the state of expressed proteins associated with CD, 3) Apply various statistical clustering and classification methods to correlate/associate microbial community composition, gene and protein content with patient metadata, including metabolite profiles and clinical phenotype. The ultimate goal of these efforts is to identify novel biomarkers for non-invasive diagnostics of CD and to eventually identify drug targets (i.e. bacterial strains) for cure or suppression of disease symptoms. PUBLIC HEALTH RELEVANCE: This study aims to unravel the contribution of the bacteria that normally inhabit the human gastrointestinal tract to Crohn's disease by using a multidisciplinary approach to study changes in the structure and function of gut microbial communities in three sets of patient cohorts who have Crohn's disease. These results will be compared with those obtained from the study of healthy individuals and have the potential to identify new biomarkers of disease severity, location, and progression.

  Study phs000257

• The Role of CDX2 in Controlling ABCB1 Expression and Chemosensitivity in Human Colon Cancer

  This study reports the results of RNA sequencing (RNA-seq) experiments performed on two human colon cancer cell lines (HT-29, DLD-1) that were genetically engineered to achieve either the constitutive expression (in HT-29 cells) or the selective inactivation (in DLD-1 cells) of the transcription factor CDX2, a master regulator of embryonic gut development and a key regulator of multi-lineage differentiation in adult colon epithelial cells. Our previous studies had identified lack of CDX2 expression as a prognostic biomarker of poor prognosis and a predictive biomarker of benefit from adjuvant chemotherapy in Stage-II and Stage-III colon cancer patients (Dalerba et al., NEJM, 374:211-222, 2016; PMID: 26789870). The molecular basis of these associations, however, remained poorly understood. The present study was designed to identify which genes are under the transcriptional control of CDX2, in order to clarify their role in shaping colon cancer resistance to chemotherapy. Of the two lines included in this study, the first (HT-29) is CDX2-negative at baseline, and was infected with a lentivirus vector encoding for the cDNA of the human CDX2 gene under the control of a constitutive promoter (CMV). The second line (DLD-1) is CDX2-positive at baseline, and was infected with a series of lentivirus vectors encoding for CRISPR/Cas9 constructs designed to selectively inactivate the CDX2 gene. In both cases, infected cells were purified by fluorescence-activated cell sorting (FACS), based on the differential expression of a green fluorescent reporter (EGFP) that was also encoded by the lentivirus vectors, and transcribed in tandem with either the CDX2 cDNA (HT-29) or CRISPR/Cas9 gRNAs targeting CDX2 (DLD-1). The experiments were designed to include 3 experimental replicates per condition, and to compare the transcriptional profile of genetically engineered cells with that of autologous controls (i.e., HT-29 cells infected with an "empty" version of the same lentivirus vector used to over-express CDX2, DLD-1 cells engineered with a "non-targeting" CRISPR/Cas9 gRNA). In the specific case of DLD-1 cells engineered with CRIPSR-Cas9 constructs targeting CDX2, infected cells were sorted based on the differential expression of EGFP and then plated as single-cells, in order to sub-culture three isogenic clones, all of which were genetically sequenced to confirm bi-allelic knock-out of the CDX2 gene, and confirmed to lack CDX2 protein expression by Western blot. The results of our experiments demonstrate that, in human colon epithelial cells, CDX2 controls the expression of the ABCB1 efflux pump, by positively regulating the transcription of the corresponding gene.

  Study phs002903

• Action to Control Cardiovascular Risk in Diabetes (ACCORD) Clinical Trial

  The Action to Control Cardiovascular Risk in Diabetes (ACCORD) clinical trial was a randomized, multicenter, double 2 x 2 factorial design study involving 10,251 middle-aged and older participants with type-2 diabetes who are at high risk for CVD events because of existing CVD or additional risk factors. The purpose of ACCORD was to determine if intensive glycemic control, intensive lipid management and intensive blood pressure control could prevent major cardiovascular events (myocardial infarction, stroke or cardiovascular death) in adults with type 2 diabetes mellitus. Secondary hypotheses included treatment differences in other cardiovascular outcomes, total mortality, microvascular outcomes, health-related quality of life and cost-effectiveness. The ACCORD trial failed to show a beneficial effect of intensive blood pressure or lipid therapy on the primary outcome, and intensive glycemia management actually increased mortality. The hypothesis underlying this ancillary study is that the failure of ACCORD to achieve its goal of reducing cardiovascular risk in diabetic patients through intensive management of hyperglycemia, dyslipidemia, and hypertension may be the result of variation in drug response due to genetic variation between individual participants. Benefits of intensive therapy may accrue to subsets of subjects with specific genetic variants predisposing to efficacious responses to particular therapeutic regimens, and harms may accrue to those with other variants predisposing to poor efficacy or adverse events. Identifying these variants could lead to a precision medicine approach to treating diabetes where each patient's genetic profile could identify the most efficacious treatment regimen with the lowest likelihood of adverse events. To test this hypothesis, a genome-wide genetic analysis was undertaken, incorporating both common variants distributed across the genome and rare variants targeted to exonic regions. Associations of genetic variants with short term responses to individual medicines as well as long term outcomes were investigated. The dataset is composed of genetic data from the ~6100 participants who agreed to participate in the ACCORD optional genetic studies and who allowed broad investigator access to their samples and the data derived from those samples, and from whom a DNA sample of sufficient quality was obtained. While a total of 8514 participants consented to the optional genetics studies, not all consented to broad investigator access, and those who did not are not included in this dataset, although they were also genotyped. Access to these additional genotypes can only be obtained by direct collaboration with the investigators of this study. Phenotype data used in the association analyses are derived from the ACCORD public release clinical data set, which has been made available through BioLINCC.

  Study phs001411

• A Clone's Genomic Stability as a Biomarker of Its DNA-Damage Resilience

  A main problem in the treatment of advanced cancers, including gastric cancers and glioblastoma, is the incertitude at which we predict how individual patients will respond to DNA-damaging agents, especially on the long run. Knowing the mechanism behind a patient's response, or the lack thereof, will help us depart from the oversimplified “more-is-better” and “one-size-fits-all” principles according to which DNA-damaging agents are administered. This will improve clinical outcome by allowing us to pinpoint those who would respond better and longer to lower doses than to higher doses of DNA-damaging agents. Under the assumption that the success of DNA-damaging therapy increases with the proliferation rate of a relatively homogeneous tumor population, there was little reason to assume anything other than monotonic dose-response relations. With the recent paradigm shift that most cancers are in fact DNA mosaic products of ongoing evolution, comes the urgency to reconsider these fundamental principles behind DNA-damaging therapy administration. As the developers of one of the first DNA deconvolution methods and with access to technologies to profile the transcriptomes of up to 10,000 cells simultaneously, we are equipped to embark on first personalized dose-finding strategies for DNA-damaging therapies. We will test the potential of the very long-term legacy that DNA-damage entails on a cell “genomic instability” as new biomarker of DNA-damage response. Our preliminary studies showed that, for most cancer types, DNA-damaging agents change a clone's genomic instability and that clones succumb to a limit in the amount of genomic instability they can tolerate. In particular, our results showed that patients with intermediate genomic instability have a very poor outcome and that this relation is only evident among treatment-naive patients, but not among patients treated with DNA-damaging agents. Further they show that we can measure genomic instability per clone and that clones with extreme genomic instability typically don't grow large. Our hypothesis that genomic instability, rather than proliferation rate, determines how sensitive a tumor is to DNA-damaging agents on the long-term, is founded on two unexpected findings. Patients with extremely high genomic instability per tumor clone have an exceptionally good outcome. Aim 1 will integrate exome- and single cell RNA-Seq (scRNA-Seq) data to characterize clones and to measure how much genomic instability they can tolerate. Low genomic instability is associated with reduced benefit from DNA-damaging agents. Aim 2 will use comet assays and treatment history to quantify DNA damage per clone, relating it to the clones' ability to tolerate DNA damage and to changes in the genomic instability of therapy-surviving clones.

  Study phs003762

• Somatic Mutation Profile by Next Generation Sequencing in HER2+ Breast Cancer

  The HER pathway is the driving force behind 30% of human breast cancers. It is important to understand how targeted therapies block different cellular pathways, and mechanisms of escape from this blockage. Therapies directed at HER2 establish a successful treatment paradigm, but de novo and acquired resistance exist. The HER signaling system is a complex network with four receptors and eleven ligands, a phosphorylation signaling cascade, and many transcription factors, all complicated by both positive and negative feedback circuits. Analysis of genomes, exomes and transcriptomes by next generation sequencing is aimed at uncovering the genetic factors responsible for patient responses to HER2-directed therapies. We are sequencing HER2-overexpressing cancers, in order to identify potential somatic changes that may better select patients who will benefit from therapy, to determine new targets that may overcome resistance, and to improve outcomes with known current HER2-targeted therapies. Whole exome capture sequencing will determine somatic mutation profiles in HER2-overexpressing tumors, to comprehensively characterize the somatic alterations, with the goal of identifying those patients most likely to respond, as well as discovering new targets that may overcome resistance to HER2-directed therapy. We augment the whole exome data with RNA-seq data to determine expression levels of somatic mutations we discover.

  Study phs000770

• Sex chromosome aneuploidies give rise to changes in the circular RNA profile

  Purpose The landscape of circular RNAs (circRNAs), an important class of non-coding RNAs that regulate gene expression, has never been described in human disorders of sex chromosome aneuploidies. We profiled circRNAs in Turner syndrome females (45,X;TS) and Klinefelter syndrome males (47,XXY; KS) to investigate how circRNAs respond to a missing or an extra X chromosome. Methods Samples of blood, muscle and fat were collected from individuals with TS (n = 33) and KS (n = 22) and from male (n = 16) and female (n = 44) controls. CircRNAs were identified using a combination of circRNA identification pipelines (CIRI2, CIRCexplorer2 and circRNA_finder). Results Differential expression of circRNAs was observed throughout the genome in TS and KS, in all tissues. The host-genes from which several of these circRNAs were derived, were associated with known phenotypic traits. Furthermore, several differentially expressed circRNAs had the potential to capture micro RNAs that targeted protein-coding genes with altered expression in TS and KS. Conclusion Sex chromosome aneuploidies introduce pervasive changes in the circRNA transcriptome, demonstrating that the genomic changes in these syndromes are more complex than hitherto thought. CircRNAs may help explain some of the genomic and phenotypic traits observed in these syndromes.

  Study EGAS00001006404

• Whole Genome Sequencing of clonal expansions of single healthy somatic cells (human, female)

  HSC73_clone: Bone marrow mononuclear cells from the healthy 73 years old female were thawed and labeled with Alexa-Fluor 488-conjugated anti-CD34 (581, Biolegend), Alexa-Fluor 700-conjugated anti-CD38 (HIT2, eBioscience), a cocktail of APC-conjugated lineage antibodies consisting of anti-CD4 (RPA-T4), anti-CD8 (RPA-T8), anti-CD11b (ICRF44), anti-CD20 (2H7), anti-CD56 (B159, all BD Biosciences), anti-CD14 (61D3), anti-CD19 (HIB19) and anti-CD235a (HIR2, all eBiocience) and 1 micro-gram/ml propidium iodide (Sigma). Using a BD FACSAria cell sorter, single Lin-CD34+CD38-PI- cells were individually sorted into low-adhesion 96-well tissue culture plates (Corning) containing 100micro-litre of StemSpan Serum-Free Expansion Medium (Stemcell technologies) supplemented with 100ng/ml of human SCF and FLT-3L, 50ng/ml of human TPO, 20ng/ml of human IL-3, IL-6 and G-CSF (all cytokines from Peprotech) and 50U/ml of penicillin and 50μg/ml of streptomycin (Sigma). Cells were incubated at 37 degrees C in a humidified atmosphere with 5% CO2 in air. After 5 days in culture, another 100micro litres of cytokine-containing medium were added. 13 days after seeding, clones B6 and G2 had expanded to approx. 105 cells and were selected for whole genome sequencing (2x101bp, paired-end, Illumina HiSeq2500) after tagmentation-based library preparation (see Extended Experimental Procedures) for clone B6 and standard library preparation for clone G2. For germline-control ~106 unsorted bone marrow mononuclear cells from the same donor were used for sequencing. An average of 30-fold sequence coverage for each the clones and the matching control were obtained.L4clone: A progenitor cell clone was raised from a peripheral blood sample of the 39 year old healthy female. Frozen peripheral blood mononuclear cells (PBMCs) were isolated from 2 ml heparinised peripheral blood via Ficoll Paque density centrifugation. A methylcellulose assay was performed as described earlier (Weisse et al., 2012). In brief, non-adherent mononuclear cells were incubated in the presence of the recombinant human cytokines IL-3, IL-5 and GM-CSF (R&D systems) over 14 days to induce colony formation. Colonies were detected under an inverted light microscope, and plucked by a pipette when colonies had approximately 10,000 cells/CFU. Each colony was washed three times in PBS and finally frozen as a cell pellet in -80 degrees C. Genomic DNA was isolated using the QIAamp DNA micro kit according to the instructions of the manufacturer (Qiagen, Hilden, Germany). Whole genome sequencing (2x101bp, paired-end, Illumina HiSeq2500) was performed for colony 4 after tagmentation-based library preparation and resulted in 15-fold sequence coverage for each the colony and the matching whole blood.

  Dataset EGAD00001000666

• Coagulation and Fibrinolysis in a Pediatric Insulin Titration Trial

  Summary of the Design and Aims: The Coagulation and Fibrinolysis in a Pediatric Insulin Titration Trial (CAF-PINT) is an ancillary study to the HALF PINT (NCT01565941) randomized, controlled trial that was originally designed to study changes in inflammatory and thrombosis pathways in the setting of an interventional trial. The parent HALF PINT study was designed to study the impact of tight glycemic control (TGC) using an explicit insulin titration algorithm to safely achieve a normal glucose in target range of 80-110 mg/dL versus 150-180mg/dL on clinical outcomes among critically ill hyperglycemic children with Heart and Lung Failure. The CAF PINT study collected blood samples pre randomization and at 2 and 4 days after randomization. Blood samples collected in EDTA tubes were centrifuged at local sites to separate the plasma from cell pellets and frozen prior to shipment. In addition, 0.5 mL of whole blood collected in PAXgene tubes modified for use in pediatrics frozen 24 hours after collection prior to shipment.Population Information: Setting: Twenty-one academic children's hospitals (HALF-PINT) and 1 academic children's hospital (IIT-SBPP). Participants: Critically ill hyperglycemic children requiring mechanical ventilation/vasoactive support (n=293).Interventions and Exposures: Patients in HALF PINT were randomized to achieve a low target glucose (80-110 mg/dL) vs. a higher target (150-180 mg/dL). The primary trial was stopped early due to low likelihood of achieving a statistically significant difference in the primary outcomes. Molecular Technologies Employed: Plasma biomarkers were assayed on thawed plasma using a Luminex panel. RNA was extracted from PAXgene Blood RNA Kits modified for pediatric use. Extracted RNA was sequenced with the Novaseq S4 system (Illumina) to generate 2 x 150 base paired end reads to a target depth of 50 million read-pairs per sample.Principal Findings of the Study: We tested for heterogeneity of treatment effects according to baseline blood inflammatory proteins using logistic regression with individual biomarkers as interaction terms and Cox regression analysis using biomarker-derived latent classes. Children with higher inflammatory proteins including TNFR-1, IL-6, IL-8, and IL-10 had lower mortality when treated with insulin to a target glucose of 80-110 mg/dL as compared to 150-180 mg/dL (interaction p<0.05). Causal forest estimates ranged from a 40% mortality reduction to a 15% mortality increase with TGC, depending on each patient’s individual biomarker profile. Latent class categorization strongly interacted with TGC with respect to mortality (Cox regression interaction p=0.005). Specifically, patients with higher inflammation benefited from TGC (6% mortality with TGC vs. 48% mortality, p=0.004) whereas hypoinflamed patients showed a trend towards harm (14% mortality with TGC vs. 7% mortality, p=0.055).Data available through dbGaP: Transcriptomic (Gene Counts) and phenotype data

  Study phs003016

• Unraveling the transcriptomic profile of utero-tubal lavage fluid of ovarian cancer patients

  The majority of ovarian cancers are diagnosed at an advanced stage and have a high mortality rate. Current screening strategies fail to improve prognosis because markers of early stage disease are lacking. This medical need justifies the search for biomarkers enabling early detection of ovarian cancer. In this study, we explore the extracellular transcriptome of utero-tubal lavage fluid obtained from 26 ovarian cancer patients and 48 controls using messenger RNA (mRNA) capture and small RNA sequencing. Because of its contact with the fallopian tube and the ovaries, collection of utero-tubal lavage might be the preferred non-invasive method to identify localized lesions compared with analyses based on peripheral blood samples. We observed an enrichment of ovarian and fallopian tube specific messenger RNAs in utero-tubal lavage fluid compared to other human biofluids. Over 300 mRNAs and 41 miRNAs were upregulated in ovarian cancer samples compared with controls. Upregulated genes were enriched for genes involved in cell cycle activation and proliferation, hinting at a tumor derived signal. In conclusion, we provide proof-of-principle that mRNA capture sequencing of utero-tubal lavage fluid is technically feasible, and that the extracellular transcriptome of utero-tubal lavage should be further explored in larger cohorts to assess the diagnostic value of the biomarkers identified in this study.

  Study EGAS00001005498

• Acute Respiratory Distress Network Early Versus Delayed Enteral Feeding to Treat People with Acute Lung Injury or Acute Respiratory Distress Syndrome (ARDSNet EDEN-BioLINCC)

  Data Access NOTE: Please refer to the “Authorized Access” section below for information about how access to the data from this accession differs from many other dbGaP accessions. Biospecimens: Access to Biospecimens is through the NHLBI Biologic Specimen and Data Repository Information Coordinating Center (BioLINCC). Biospecimens from ARDSNet-EDEN include bronchial lavage, plasma, DNA, and urine. Please note that use of biospecimens in genetic and non-genetic research is subject to a tiered consent. Objectives: To determine if initial lower-volume trophic enteral feeding would increase ventilator-free days (VFDs) and decrease gastrointestinal intolerances compared with initial full enteral feeding. Background: Mechanically ventilated patients cannot eat normally and if not fed for long periods become malnourished. Because malnutrition is associated with poor outcomes in critically ill patients, artificial nutrition is often provided, especially in those with acute lung injury (ALI) and with expected longer duration of mechanical ventilation. When feasible, enteral nutrition targeting full caloric needs has been advocated over parenteral nutrition. However, feeding intolerance and common care practices often serve as practical barriers to reaching recommended goals. Although confounded by indication and severity of illness, several observational studies have shown improved clinical outcomes, including fewer infections, shorter duration of mechanical ventilation, and lower mortality for patients receiving a higher percentage of calculated caloric needs. Nonetheless, the best timing, formulation, and amount of enteral nutrition remain unknown. Participants: A total of 1,000 participants were enrolled. Design: Participants were randomized, stratified by site and presence of shock at enrollment, to receive either trophic or full enteral feeding for the first 6 days of mechanical ventilation. The initial 272 patients were also simultaneously randomized to a separate trial (the OMEGA study) comparing a nutritional supplement containing omega-3 fatty acids and antioxidants with an isocaloric, isovolemic control in a 2 × 2 factorial design.The designated feeding strategy was initiated within 6 hours of randomization and continued until death, extubation, or day 6. The care of mechanically ventilated patients still receiving enteral feedings after day 6 was managed according to the full feeding strategy in both groups. In extubated patients who then required reintubation, enteral nutrition was restarted and managed according to the study protocol.In the full-feeding group, enteral nutrition was initiated at 25 mL/h and advanced to goal rates as quickly as possible. Gastric residual volumes were checked every 6 hours while enteral feeding was increased. Patients randomized to the initial trophic-feeding group had enteral nutrition initiated at 10 mL/h (10-20 kcal/h) for the first 272 patients who also received the omega-3 or control supplement (240 mL volume per day). After the data and safety monitoring board stopped the OMEGA portion of the factorial design, the initial trophic feeding rate was changed to 20 kcal/h to approximate the calories that had been delivered in the OMEGA study. Gastric Residual Volumes (GRVs) were checked every 12 hours during trophic feeding. In patients randomized to trophic feeding, enteral nutrition was advanced to full-energy feeding rates following the same protocol used for the full-feeding group if they were still receiving mechanical ventilation at 144 hours. Conclusions: There was no difference between groups with regard to the primary end point, VFDs to day 28. There also were no differences in 60-day mortality, organ failure−free days, ICU-free days, or the incidence of infection between groups. Similarly, there were no differences between groups in VFDs or survival when analyzed by body mass index category or when subsets of patients with shock or more severe lung injury (acute respiratory distress syndrome) were examined.JAMA. 2012 Feb 22;307(8):795-803

  Study phs004168

• Korea Epigenome Project(KEP), Korea National Research Institute of Health(KNIH)

  The KEP was launched in 2012 to produce 50 epigenomic sets on Korean chronic diseases related target cells with the participation on Internationsl Human Epigenome Consortium(IHEC) of which focus is constructing 1000 reference epigenome map on 250 human cell type. The objective of the KEP was to provide a comprehensive epigenomic map of Korean chronic disease-related target tissues. Single cells were isolated from the pancreas, fat and kidney. Islet, ductal, acinar, and beta cells from the pancreas were purified into a single homogenous cell type using flow-assisted cell sorting technique. Adipocytes and pre-adipocytes were purified from fat, using Percoll gradient centrifugation technique. Mesangial, distal tubule, proximal tubule, podocytes, and collecting duct cells from the kidney were collected using mechanical sieving and immunostaining methods. Each cell type was derived from normal tissue or those with diabetes, obesity, or chronic kidney disease. Sixty reference epigenome datasets of pancreatic, fat, kidney cells were produced. Cell types such as the islets, and beta from the pancreas, adipocytes, pre-adipocytes from fat tissue, podocyte, mesangial, tubules from kidney were used for primary target to produce a Korean metabolic reference epigenome. And the bisulfite converted-whole genome sequencing, Infinium 450k DNA methylation bead array, mRNA, miRNA-Sequencing were performed to produce a comprehensive epigenomic map. The aim of roadmap reports is preparation of Long-term and Mid-term Technology Roadmaps for Korean Human Epigenome Project towards a competitive execution of International Epigenome Project and successful expansion of the research outcomes. Roadmap reports reviewed on recent technological development of DNA methylation and histone modification analysis, emerging new trends in epigenome medicine, plans for hardware and software analysis systems aimed to epigenome analysis, medical expert opinion on target cells and isolation scheme, bioinformatics expert opinion on analytical tools, new ideas for a putative epigenome studies, trends in single cell epigenome analysis. In conclusion, roadmap report suggested several projects for the future KNIH epigenome research in collaboration with epigenome research group. These are including refocusing the goals of KNIH-driven research in identification of disease mechanism, development of new studies on GWAS-EWA, higher priority on the epigenomic studies of human primary cells, establishment of collaborative relationship with tissue banking centers, establishment of mapping centers and data analysis portals, establishment of Korean Epigenome Research Consortium, establishment of steering committee governing both research and funding, launching of additional funding for epigenomic researches.

  Study EGAS00001001774