22929 results for "ega"

in 38.93 milliseconds.

• Circulating tumor DNA dynamics reveal KRAS G12C mutation heterogeneity and response to treatment with the KRAS G12C inhibitor divarasib in solid tumors

  Purpose: To inform prognosis, treatment response, disease biology, and KRAS G12C mutation heterogeneity, we conducted exploratory circulating tumor DNA (ctDNA) profiling on 134 patients with solid tumors harboring a KRAS G12C mutation treated with single-agent divarasib (GDC-6036) in a phase 1 study. Experimental design: Plasma samples were collected for serial ctDNA profiling at baseline (Cycle 1 Day 1 prior to treatment) and multiple on-treatment time points (Cycle 1 Day 15 and Cycle 3 Day 1). Results: KRAS G12C ctDNA was detectable from plasma samples in 72.9% (43/59) and 92.6% (50/54) of patients with non-small cell lung cancer and colorectal cancer, respectively, the majority of whom were eligible for study participation based on a local test detecting the KRAS G12C mutation in tumor tissue. Baseline ctDNA tumor fraction was associated with tumor type, disease burden, and metastatic sites. A decline in ctDNA level was observed as early as Cycle 1 Day 15. Serial assessment showed a decline in ctDNA tumor fraction associated with response and progression-free survival. Except for a few cases of KRAS G12C sub-clonality, on-treatment changes in KRAS G12C variant allele frequency mirrored changes in the overall ctDNA tumor fraction. Conclusion: Across tumor types, the KRAS G12C mutation likely represents a truncal mutation in the majority of patients. Rapid and deep decline in ctDNA tumor fraction was observed in patients responding to divarasib treatment. Early on-treatment dynamics of ctDNA were associated with patient outcomes and tumor response to divarasib treatment.

  Study EGAS50000000315

• Lactate metabolism in cancer stem cell fate regulation

  Tumors arise as a result of uncontrolled cell proliferation driven by mutations in the signaling pathways that regulate the balance between cell division and differentiation. Despite of that, in the tumor, cancer stem cells seem to retain their capability to differentiate, which can lead to differential chemotherapy response. We previously showed that intestinal stem cells and differentiated cells display specific metabolic features and that they interact through metabolites. Moreover, we found that altering mitochondrial metabolism is sufficient to lead and define stem cell differentiation. In cancer, cellular metabolism is highly derailed. However, how metabolic changes in the cancer cell contribute to cell fate decisions remains unknown. Recently, the crosstalk between metabolism and epigenetics has emerged as one of the players in tumor development. Here we aim to address the impact of the altered metabolism of cancer cells on their epigenome as a regulatory axis of cell identity and fate. We use human-derived tumor organoids where we introduce genetic encoded fluorescent-based reporters to visualize cell identity and metabolism at single cell level. We combine 4D live imaging with a recently developed machine learning method to re-construct cell-lineages, trace stem cell state and metabolic changes during tumor organoid development. We find that, mimicking the tumor metabolic microenvironment (Warburg effect-driven low-glucose high-lactate) leads to changes cellular metabolism and histone acetylation. Strikingly, these alterations provoke profound changes in tumor development by increasing the population of cancer stem cells and their proliferation and differentiation dynamics.

  Study EGAS50000000063

• Tumor Evolution Analysis Uncovered Immune-Escape Related Mutations in Relapsed Diffuse Large B-Cell Lymphoma

  More than 60% of patients with primary diffuse large B-cell lymphoma (DLBCL) achieve durable remissions upon the chemoimmunotherapy with R-CHOP. However, the outcome of the relapse patients is extremely poor. To characterize the molecular features that underlie treatment resistance, we performed whole-exome sequencing (n=44) and deep targeted sequencing (n=27) of 44 diagnosis-relapse DLBCL pairs. We found that the overall landscape of genetic alterations was concordant between diagnosis and relapse DLBCL. By a deep comparison of the longitudinal sequencings, we identified relapse-enriched mutations that involved in evasion of immune surveillance (CD58, TNFRSF14, PKD1, FAS), resistance to DNA damaging agents (ATM, TP53), and critical tumor suppressor gene (ZFHX3). Knockdown of ATM and ZFHX3 resulted in increased drug resistance in both germinal center B cell like (GCB) and activated B cell like (ABC) DLBCL cell lines in vitro, while the knockdown of TNFRSF14 and FAS conferred resistance specifically in GCB cell lines. Furthermore, we summarized two evolutionary patterns of DLBCL relapse as the “durable retention” model with stable clone structure, and the “selective sweep” model which exhibited significant clonal switch. Of note, genetic alterations leading to attenuated immune responses usually presented as trunk mutations in the “durable retention” model (7/12) while were newly acquired or expanded upon relapse in the “selective sweep” model (9/16). Overall, our study demonstrated that relapse DLBCL derives from distinct clonal evolutionary patterns, and highlighted immune escape as an important contributor of tumor relapse.

  Study EGAS50000000032

• Surgery in combination with immune checkpoint therapy as an effective treatment for patients with metastatic cancer.

  Abstract Surgical removal of primary tumors was shown to reverse tumor-mediated immune suppression in pre-clinical models with metastatic disease. However, how cytoreductive surgery in the metastatic setting modulates the immune responses in patients, especially in the context of immune checkpoint therapy (ICT)-containing treatments is not understood. Here, we report the first prospective, non-comparative clinical trial (N=104) using three different ICT-containing strategies plus cytoreductive or “debulking” surgery to remove the primary tumor-bearing kidney or a metastatic lesion as a treatment for patients (N=43) with metastatic clear cell renal cell carcinoma (mccRCC). For those patients (N=61) who were not candidates for cytoreductive surgery, a biopsy was obtained instead for correlative biological studies. Our data demonstrated that the combination of ICT with cytoreductive surgery was safe and feasible in patients with mccRCC. The 2-year overall survival was 84% with a median OS of 54.7 months for patients who received ICT containing regimens plus surgery. Immune-monitoring studies with co-detection by indexing (CODEX) identified distinct tumor spatial conformation of cellular subsets as a novel and improved predictor of response to ICT. Importantly, single-cell RNA-sequencing data demonstrated that surgical removal of the tumor increased antigen-presenting dendritic cell population with a concurrent reduction in KDM6B expressing immune-suppressive myeloid cells in the peripheral blood. Together, this study highlighted the feasibility of combining ICT with cytoreductive surgery in a metastatic setting and demonstrated potential enhancement of immune responses following ICT plus cytoreductive surgery in patients with metastatic disease.

  Study EGAS00001005667

• Impact of Glycolysis Inhibition on the Epigenome of Synovial Fluid T Cells in Juvenile Idiopathic Arthritis

  Juvenile Idiopathic Arthritis (JIA) encompasses a group of childhood-onset autoimmune conditions characterized by joint inflammation and mononuclear cell infiltration, especially of activated CD4+ memory/effector T (Tmem/Teff) cells in synovial fluid (SF). JIA CD4+ T cells exhibit a distinct inflammation-associated epigenomic profile, though the underlying molecular mechanisms remain unclear. This study investigates the role of glycolysis in shaping the epigenomic landscape of JIA SF CD4+ T cells, with a focus on H3K27ac histone modifications. Using an ex vivo model of antigen stimulation, CD4+ T cells were isolated from peripheral blood (PB) of healthy controls (HC) and synovial fluid from JIA patients. Cells were activated in vitro with anti-CD3/CD28 for 24 hours, and changes in H3K27ac were assessed by chromatin immunoprecipitation followed by sequencing (ChIP-seq). To explore the specific role of glycolysis, cells were also treated with 2-deoxyglucose (2DG), an inhibitor of glucose-6-phosphate production. This inhibition impacted both H3K27ac acetylation patterns and gene expression associated with the activated epigenomic regions. Additionally, HC CD4+ T cells were activated in vitro in the presence of 30% SF or 30% plasma from JIA patients. RNA sequencing (RNA-seq) was performed to examine SF-induced transcriptional changes, providing insight into the metabolic dysregulation and epigenomic alterations observed in JIA SF CD4+ T cells. This dataset includes ChIP-seq and RNA-seq data that collectively advance our understanding of how glycolysis and the SF environment contribute to the epigenetic and transcriptional regulation in JIA CD4+ T cells.

  Study EGAS50000000808

• Circulating Tumor DNA Analysis in ERBB2-Amplified Colorectal Cancer: Biomarker Analysis of the MyPathway Trial

  Combination of two HER2-directed antibodies, pertuzumab and trastuzumab (P+T), has antitumor activity in HER2-positive colorectal cancer (CRC). Although liquid biopsies are increasingly being used in clinical oncology, the association between tumor and circulating tumor DNA (ctDNA) ERBB2 status and ctDNA monitoring for early response and resistance are unknown. Eighty-five patients with ERBB2-amplified and/or overexpressed CRC were treated with P+T in the MyPathway trial; 42 had ctDNA testing at Cycle 1 Day 1 (C1D1), and 38 had longitudinal plasma tested for ctDNA. We analyzed the ctDNA vs tissue ERBB2 concordance, genomic co-alterations, and ctDNA dynamics and association with response. Forty-one (98%) of 42 patients had genomic alterations detected in ctDNA at C1D1, and 29 (69%) had ERBB2 amplification in ctDNA. There was a strong correlation between ERBB2 copy number on next-generation sequencing in tissue and C1D1 ERBB2 ctDNA copy number. Thirty-seven percent achieved a molecular response by C3D1 on P+T, which was associated with prolonged progression-free survival and overall survival (OS). CDKN2A and KRAS mutations were associated with shorter OS, and a trend was seen with PIK3CA mutations. Several emerging co-alterations were identified in ctDNA at progression, including in the MAPK and PI3K pathways and other tyrosine receptor kinases. ctDNA can detect ERBB2 amplification in many, but not all, patients with ERBB2 amplification detected in tumor samples. ctDNA molecular response was associated with better survival, and ctDNA co-alterations may offer insights into mechanisms of intrinsic and acquired resistance.

  Study EGAS50000000916

• Nanopore sequencing of blood, saliva, buccal mucosa samples of patients with inherited retinal disorders

  X-linked genetic disorders typically affect females less severely than males due to the presence of a second X chromosome not carrying the deleterious variant. However, the phenotypic expression in females is highly variable, which may be explained by an allelic skew in X chromosome inactivation. Accurate measurement of X inactivation skew is crucial to understand and predict disease phenotype in carrier females, with prediction especially relevant for degenerative conditions. We propose a novel approach using nanopore sequencing to quantify skewed X inactivation accurately. By phasing sequence variants and methylation patterns, this single assay reveals the disease variant, X inactivation skew, its directionality, and is applicable to all patients and X-linked variants. Enrichment of X-chromosome reads through adaptive sampling enhances cost-efficiency. Our study includes a cohort of 16 X-linked variant carrier females affected by two X-linked inherited retinal diseases: choroideremia and RPGR-associated retinitis pigmentosa. As retinal DNA cannot be readily obtained, we instead determine the skew from peripheral samples (blood, saliva and buccal mucosa), and correlate it to phenotypic outcomes. This revealed a strong correlation between X inactivation skew and disease presentation, confirming the value in performing this assay and its potential as a way to prioritise patients for early intervention, such as gene therapy currently in clinical trials for these conditions. Our method of assessing skewed X inactivation is applicable to all long-read genomic datasets, providing insights into disease risk and severity and aiding in the development of individualised strategies for X-linked variant carrier females.

  Study EGAS50000000440

• PCR-free shallow whole genome sequencing for chromosomal copy number detection from plasma of cancer patients is an efficient alternative to the conventional PCR-based approach

  Somatic copy number alterations (SCNAs) can be detected in tumor cell-free DNA (cfDNA) by shallow whole genome sequencing (sWGS). Polymerase chain reaction (PCR) is typically included in library preparations, but a PCR-free method could be an effective alternative for research and diagnostics. Peripheral blood samples (n=22) were collected in EDTA-containing tubes from non-small cell lung cancer patients (n=10) and healthy donors (n=12). PCR and PCR-free library preparations underwent sWGS and the following metrics were compared: mapping quality, proportion of duplicate reads, genome coverage, sequencing read length and copy number profiles. The PCR-free method was further evaluated on archived plasma samples (n=33) of varying pre-analytical quality collected in EDTA or lithium heparin tubes. Archived samples were from acute myeloid leukaemia (n=15) and B-cell lymphoma (n=18) patients. The percentage of unique reads was significantly higher for the 22 samples analysed by PCR-free (95.7%) compared to PCR (84%). PCR-free libraries from 33 archived samples produced data of a comparable quality, with 95.4% unique reads. All other evaluated metrics were highly comparable for PCR and PCR-free library preparations. We conclude that SCNA detection can be performed by PCR-free sWGS on cfDNA from cancer patients and is feasible for diagnostic and archival plasma collected in EDTA- or lithium heparin-containing tubes. Our results pave the way for an automated cfDNA workflow for cancer plasma samples.

  Study EGAS00001004692

• Coding and regulatory somatic profiling of triple-negative breast cancer in Sub-Saharan African patients

  The burden of triple-negative breast cancer (TNBC) may be shaped by genetic factors, particularly inherited and somatic mutation profiles. However, data on this topic remain limited, especially for the African continent, where a higher TNBC incidence is observed. In the age of precision medicine, cataloguing TNBC diversity in African patients becomes imperative. We performed whole exome sequencing, including untranslated regions, on 30 samples from Angola and Cape Verde, which allowed to ascertain on potential regulatory mutations in TNBC for the first time. A high somatic burden was observed for the African cohort, with 86% of variants being so far unreported. Recurring to predictive functional algorithms, 17% of the somatic single nucleotide variants were predicted to be deleterious at the protein level, and 20% overlapped with candidate cis-regulatory elements controlling gene expression. Several of these somatic functionally-impactful mutations and copy number variation (mainly in 1q, 8q, 6 and 10p) occur in known BC- and all cancer-driver genes, enriched for several cancer mechanisms, including response to radiation and related DNA repair mechanisms. TP53 is the top of these known BC-driver genes, but our results identified possible novel TNBC driver genes that may play a main role in the African context, as TTN, CEACAM7, DEFB132, COPZ2 and GAS1. These findings emphasize the need to expand cancer omics screenings across the African continent, the region of the globe with highest genomic diversity, accelerating the discovery of new somatic mutations and cancer-related pathways.

  Study EGAS50000001013

• Whole-genome and transcriptome versus panel sequencing in precision oncology: A translational-clinical comparison

  Precision oncology offers new cancer treatment options. While the use of whole-genome (WGS) and transcriptome (TS) sequencing remains often confined to specific clinical programs, molecular analysis with gene panels is more common clinical practice. Methods: Twenty patients with rare or advanced tumors sequenced by WGS/TS within the DKFZ/NCT/DKTK-MASTER program from 2015 to 2020 (MASTER 1) were additionally sequenced with the TSO500 DNA and TST170 RNA Assay by Illumina. Original Molecular Tumor Board (MTB) therapy recommendations (TR) were juxtaposed with TR from NGS panel data analysis (comparison 1). To mitigate temporal differences between the two sequencings, WGS/TS was bioinformatically reanalyzed (MASTER 2), and new TR were again compared (comparison 2). Results: In comparison 1, 47.1% of the MASTER 1’s TR and 62.8% of the panel’s TR were identical, as were 45.9% of the TR of MASTER 2 and 54.9% of the panel in comparison 2. In both MASTER 1 and MASTER 2, 26.5% and 36.1% of TR, respectively, relied on Biomarkers (BM) which were not covered by the panel. Eight of the ten MTB-derived therapy implementations were also supported in the panel, the two remaining ones being based on BM absent from the panel. Conclusions: About half of the TR in WGS/TS and panel sequencing were identical. Around a third of TR issued by WGS/TS were based on BM which are not covered by the panel, leading to therapy implementations. The presented data underscores the clinical value of comprehensive molecular analyses in individuals with advanced cancers.

  Study EGAS50000000431

• Mutational signatures in normal tissues of patients treated with clinical stage G-quadruplex binder CX5461

  Drug induced photoreactions following exposure to UVA have been documented for several therapeutic compound classes, including psoralens and fluoroquinolones. CX5461, a quinolone-derived small molecule currently in phase 1 studies for cancer treatment, has been shown to be a clinically significant UV photosensitizer. Here we examined the mutagenic spectrum of CX5461 in normal tissues of patients undergoing treatment, by comparing deep whole genome sequencing (WGS) from skin biopsies taken from non-sun exposed areas pre and post drug exposure of patients in the Canadian Cancer Clinical Trials Group (CCTG) IND.231 trial (NCT02719977). We show that a low burden of SNVs (<100) and indels (<40) mutations was detected in post CX5461 skin biopsies, with a pattern resembling background or sequencing artefact signatures. To further characterize the mutational pattern of CX5461 photoreactivity, we compared mutational profiles of human retinal pigment epithelial cells exposed to UVA, CX5461, or a combination of both UV and CX5461. Treatment with CX5461 or UVA alone, resulted in a low SNV burden. In contrast, human cells exposed to UVA and CX5461 had a two to four orders higher mutational burden, with a distinct mutational pattern characterized by T>A, T>C mutations, dissimilar to that of UV alone or CX5461 alone. Our results suggest that human exposure to CX5461 alone is not associated with high levels of single base mutagenesis. However, CX5461 is a very potent UV photosensitizer, underscoring the need for adequate UV screening and controls when studying the effects of photosensitizing small molecules on single base DNA variation.

  Study EGAS50000001145

• Transcriptome profiling of megakaryocytes and platelets: application to GP9- and IKZF5-related thrombocytopenia

  Platelets are anucleate cells produced in the bone marrow and derived from large progenitor cells called megakaryocytes (MKs). Platelets receive RNA transcripts from their progenitorial MKs during thrombopoiesis. However, the correspondence between platelet and MK transcriptomes is poorly understood, particularly in the context of germline mutations that cause platelet formation defects or thrombocytopenia. We have studied the effects of two such mutations on MK and platelet transcriptomes. We generated immortalized MK cell line (imMKCL)-based models of Bernard-Soulier syndrome and IKZF5-related thrombocytopenia. MKs derived from imMKCLs with either a homozygous deletion of GP9 (GP9-/-) or a heterozygous Y121F variant in IKZF5 (IKZF5WT/Y121F) exhibited reduced proplatelet formation (reductions of 96% and 57%, respectively). Platelets from patients with either GP9-/- or IKZF5WT/Y121F genotypes had broad transcriptomic dysregulation, suggesting that (pro)platelet formation defects due to mutations in glycoprotein receptor and transcription factor genes such as GP9 and IKZF5 affect the trafficking of transcripts from MKs into platelets. RNA-seq data from MKs at four stages of differentiation revealed widespread but distinct changes in expression over time between the GP9-/- and the IKZF5WT/Y121F genotypes. Dysregulated genes in GP9-/- MKs were enriched for RNA metabolism and actin/tubulin folding pathways while those in IKZF5WT/Y121F MKs were enriched for cell cycle pathways. Most of these genes were also dysregulated in the platelets of patients with the corresponding diseases. Our results suggest that patients with inherited forms of thrombocytopenia present with a transcript trafficking defect during platelet formation.

  Study EGAS50000001276

• Complex genotype-phenotype relationships shape the response to treatment of Down Syndrome Childhood Acute Lymphoblastic Leukaemia

  Extensive genetic and epigenetic variegation has been demonstrated in many malignancies. Importantly, their interplay has the potential to contribute to disease progression and treatment resistance. To shed light on the complex relationships between these different sources of intra-tumour heterogeneity, we explored their relative contributions to the evolutionary dynamics of Acute Lymphoblastic Leukaemia (ALL) in children with Down syndrome, which has particularly poor prognosis. We quantified the tumour propagating potential of genetically distinct sub-clones using serial transplantation assays and SNP-arrays. While most leukaemias were characterized by a single dominant subclone, others were highly heterogeneous. Importantly, we provide clear and direct evidence that genotypes and phenotypes with functional relevance to leukemic progression and treatment resistance can co-segregate within the disease. Hence, individual genetic lesions can be restricted to well-defined cell immunophenotypes, corresponding to different stages of the leukemic differentiation hierarchy and varied proliferation potentials. As a result of this difference in fitness, which can be accurately quantified via competitive transplantation assays, matching diagnostic, post-treatment, and relapse leukaemias can be dominated by different genotypes, including pre-leukemic clones persisting throughout the disease progression and treatment. Intriguingly, plasticity also appears to be a temporally defined property that can segregate with genotype. These results suggest that Down Syndrome ALL should be viewed as a complex matrix of cells exhibiting genetic and epigenetic heterogeneity that foster extensive clonal evolution and competition. Therapeutic intervention reshapes this ‘eco-system’ and may provide the right conditions for the preferential expansion of selected compartments and subsequently relapse.

  Study EGAS50000001287

• Heterogeneity of IKZF1 genomic alterations and risk of relapse in childhood B-cell precursor acute lymphoblastic leukemia

  Purpose Genomic alterations of IKZF1 are common and associated with adverse clinical features in B-ALL, but the relationship between the type of IKZF1 alteration, disease subtype and outcome are incompletely understood. Patients and Methods We used transcriptome and genomic sequencing, and SNP microarray data to determine leukemia subtype and genomic alterations for 688 pediatric patients with B-ALL enrolled into St. Jude Total Therapy 15 and 16 studies. Results IKZF1 alterations were identified in 115 (16.7%) patients, commonly in BCR::ABL1 (78%) and CRLF2-rearranged, BCR::ABL1-like B-ALL (70%), and were associated with five-year cumulative incidence of relapse (CIR) of 14.8 ± 3.3% compared to 5.0 ± 0.9% for patients without any IKZF1 alteration (P < 0.0001). IKZF1 deletions of exon 4-7 (P = 0.0002), genomic IKZF1plus with any IKZF1 deletion (P = 0.006) or with focal IKZF1 deletion (P = 0.0007), and unfavorable genomic subtypes (P < 0.005) were independently adversely prognostic. Patients with both IKZF1 exon 4-7 deletion and unfavorable genomic subtype had a significantly increased risk of relapse (HR = 58.3; 95% CI, 11.9 – 285.4, P < 0.0001) whereas IKZF1 sequence mutations were not independently predictive of outcome. Conclusions Genomic IKZF1plus with any IKZF1 deletion, IKZF1 deletion of exon 4-7, and unfavorable subtype confer increased risk of relapse. The combination of unfavorable genomic subtype and IKZF1 deletion of exon 4-7 identified patients at greatest risk of relapse despite MRD-directed therapy. The type of IKZF1 alteration together with subtype are informative for risk stratification and predict response in patients with B-ALL.

  Study EGAS50000000106

• WES analysis of a mixed cohort of pituitary tumors

  BACKGROUND: Sporadic pituitary adenomas are the second most common primary brain tumor, but their genetics remains incompletely understood. We performed deep DNA sequencing of 42 pituitary macroadenomas, including hormonally active and atypical adenomas. METHODS: Pituitary adenomas representative of a spectrum of histopathologic subtypes, primary and recurrent tumors, and MIB-1 proliferative indices were profiled. Whole-exome sequencing was performed on tumors and paired blood followed by analysis of somatic copy number alteration and mutation. RESULTS: On the whole, the tumors harbored a low rate of mutations, similar to other benign brain tumors. However, almost one-third (29%) exhibited somatic copy-number alterations across large fractions of the genome (up to 99%). Despite such widespread genomic disruption, these tumors had few focal events, which is atypical among highly-disrupted cancers. The other 71% of tumors formed a distinct molecular class, with somatic copy-number alterations involving less than 6% of the genome. Among the highly-disrupted group, 92% were functional or atypical null-cell adenomas, whereas 87% of the less-disrupted group were clinically non-functional adenomas. The only significantly recurrent events were arm-level losses of chromosomes 1, 2, 11, and 18. No significantly recurrent mutations or insertion/deletion events were identified, suggesting no genes are altered by exomic mutations in &gt;15% of pituitary tumors. CONCLUSIONS: Genomic characterization of a large series of sporadic pituitary adenomas reveals two classes of tumor, based on genomic disruption from arm-level copy-number alterations, which associate with hormonal and histologic subtypes.

  Study EGAS00001001714

• Genome-Wide Association Study of aspirin-induced PUD in a UK cohort

  Aspirin is a widely prescribed non-steroidal anti-inflammatory drug used as an effective analgesic and at low doses as anti-platelet drug for the prevention of cardiovascular diseases. However, aspirin is associated with a significant incidence of upper gastrointestinal (GI) complications, notably peptic ulceration. This study aims to identify genetic risk factors associated with the pathogenesis of aspirin-induced upper GI ulceration.A cohort of 248 patients with aspirin-induced, endoscopiclly confirmed GI ulceration and 491 controls (ulcer, without aspirin or aspirin without ulcer) were genotyped for 1.524,956 SNPs using the Illumina Omni 2.5 SNP array. A replication cohort (206 NSAID-induced ulcer cases [of which 87 were aspirin] and 309 controls) was subsequently genotyped for identified SNP association signals using Sequenom iPLEX Massarry. Logistic regression analysis of identified two SNP signal with suggestive association with aspirin-induced GI ulceration. A SNP within the EYA1 locus (rs12678747) was found to be associated with aspirin GI ulceration in both the discovery cohort (p=3.79x10-7, OR=2.06 ([1.55-2.73]) and, to a lesser extent, in the replication cohort (p=.0005, OR=1.71[1.17-2.50]. Combined analysis of both cohorts yielded a genome wide significant association (p=1.46x10-8). The rs12678747 genotype was found to be associated with significantly higher EYA1 transcript in ulcer biopsy tissue (p&lt;0.05).Data suggest that EYA1 may represent a risk locus for aspirin-induced GI ulceration. Though it is thought that EYA1 may play a key role in mediating cell death/ survival

  Study EGAS00001002052

• Microsatellite data show recent demographic expansions in sedentary but not in nomadic human populations in Africa and Eurasia.

  The transition from hunting and gathering to plant and animal domestication was one of the most important cultural and technological revolutions in human history. According to archeologists and paleoanthropologists, this transition triggered major demographic expansions. However, few genetic studies have found traces of Neolithic expansions in the current repartition of genetic polymorphism, pointing rather toward Paleolithic expansions. Here, we used microsatellite autosomal data to investigate the past demographic history of 87 African and Eurasian human populations with contrasted lifestyles (nomadic hunter-gatherers, semi-nomadic herders and sedentary farmers). Likely due to the combination of a higher mutation rate and the possibility to analyze several loci as independent replicates of the coalescent process, the analysis of microsatellite data allowed us to infer more recent expansions than previous genetic studies, potentially resulting from the Neolithic transition. Despite the variability in their location and environment, we found consistent expansions for all sedentary farmers, while we inferred constant population sizes for all hunter-gatherers and most herders, which could result from constraints linked to a nomadic or semi-nomadic lifestyle and/or competition for land between herders and farmers. As an exception, we inferred expansions for Central Asian herders. This might be linked with the arid environment of this area, which may have been more favorable to nomadic herders than to sedentary farmers. Alternatively, current Central Asian herders may descent from populations who have first experienced a transition from hunter-gathering to sedentary agropastoralism, and then a second transition to nomadic herding.

  Study EGAS00001000652

• Exome sequencing of uterine leiomyosarcomas

  Uterine leiomyosarcomas (ULMSs) are aggressive smooth muscle tumors associated with poor clinical outcome. Despite previous cytogenetic and molecular studies, their molecular background has remained elusive. To examine somatic variation in ULMS, we performed exome sequencing on 19 tumors. Altogether 43 genes were mutated in at least two ULMSs. Most frequently mutated genes included tumor protein P53 (TP53; 6/19; 33%), alpha thalassemia/mental retardation syndrome X-linked (ATRX; 5/19; 26%), and mediator complex subunit 12 (MED12; 4/19; 21%). Unlike ATRX mutations, both TP53 and MED12 alterations have repeatedly been associated with ULMSs. All the observed ATRX alterations were either nonsense or frameshift mutations. ATRX protein levels were reliably analyzed by immunohistochemistry in altogether 44 ULMSs, and the majority of tumors (23/44; 52%) showed clearly reduced expression. Loss of ATRX expression has been associated with alternative lengthening of telomeres (ALT), and thus the telomere length was analyzed with telomere-specific fluorescence in situ hybridization. The ALT phenotype was confirmed in all ULMSs showing diminished ATRX expression. Exome data also revealed one nonsense mutation in death-domain associated protein (DAXX), another gene previously associated with ALT, and the tumor showed ALT positivity. Aberrant expression of both TP53 and ATRX were associated with poor overall survival. In conclusion, exome sequencing revealed that TP53, ATRX, and MED12 are frequently mutated in ULMSs. ALT phenotype was commonly seen in tumors, indicating that ATR inhibitors, which were recently suggested as possible new drugs for ATRX-deficient tumors, could provide a potential novel therapeutic option for ULMS.

  Study EGAS00001001612

• Recurrent epimutations activate gene body promoters in primary glioblastoma

  Aberrant DNA hypomethylation may play an important role in the growth rate of glioblastoma (GBM), but the functional impact on transcription remains poorly understood. We assayed the GBM methylome with MeDIP-seq and MRE-seq, adjusting for copy number differences, in a small set of non-glioma CpG island methylator phenotype (non-G-CIMP) primary tumors. Recurrent hypomethylated loci were enriched within a region of chromosome 5p15 that is specified as a cancer amplicon and also encompasses TERT, encoding telomerase reverse transcriptase, which plays a critical role in tumorigenesis. Overall, 76 gene body promoters were recurrently hypomethylated, including TERT and the oncogenes GLI3 and TP73. Recurring hypomethylation also affected previously unannotated alternative promoters, and luciferase reporter assays for three of four of these promoters confirmed strong promoter activity in GBM cells. Histone H3 lysine 4 trimethylation (H3K4me3) ChIP-seq on tissue from the GBMs uncovered peaks that coincide precisely with tumor-specific decrease of DNA methylation at 200 loci, 133 of which are in gene bodies. Detailed investigation of TP73 and TERT gene body hypomethylation demonstrated increased expression of corresponding alternate transcripts, which in TP73 encodes a truncated p73 protein with oncogenic function and in TERT encodes a putative reverse transcriptase-null protein. Our findings suggest that recurring gene body promoter hypomethylation events, along with histone H3K4 trimethylation, alter the transcriptional landscape of GBM through the activation of a limited number of normally silenced promoters within gene bodies, in at least one case leading to expression of an oncogenic protein.

  Study EGAS00001000685

• DNA methylation and the adverse metabolic outcomes of adiposity

  Overweight and obesity affect ~1.5 billion people worldwide, and are major risk factors for type-2 diabetes (T2D), cardiovascular disease and related metabolic and inflammatory disturbances.1,2 Although the mechanisms linking adiposity to its clinical sequelae are poorly understood, recent studies suggest that adiposity may influence DNA methylation,3-6 a key regulator of gene expression and molecular phenotype.7 Here we use epigenome-wide association to show that body mass index (BMI, a key measure of adiposity) is associated with widespread changes in DNA methylation (187 genetic loci at P&lt;1x10-7, range P=9.2x10-8 to 6.0x10-46; N=10,261 samples). Genetic association analyses demonstrate that the alterations in DNA methylation are predominantly the consequence of adiposity, rather than the cause. We find the methylation loci are enriched for functional genomic features in multiple tissues (P&lt;0.05), and show that sentinel methylation markers identify gene expression signatures at 38 loci (P&lt;9.0x10-6, range P=5.5x10-6 to 6.1x10-35, N=1,785 samples). The methylation loci identified highlight genes involved in lipid and lipoprotein metabolism, substrate transport, and inflammatory pathways. Finally, we show that the disturbances in DNA methylation predict future type-2 diabetes (relative risk per 1SD increase in Methylation Risk Score: 2.3 [2.07-2.56]; P=1.1x10-54). Our results provide new insights into the biologic pathways influenced by adiposity, and may enable development of new strategies for prediction and prevention of type-2 diabetes and other adverse clinical consequences of obesity.

  Study EGAS00001001922

• Methylation differences in trisomy 21 using monozygotic twins

  We have hypothesized that methylation differences induced by trisomy 21 (T21) contribute to the phenotypic characteristics and heterogeneity in T21. In order to determine the methylation differences in T21 without the interference of the interindividual genomic variation, we have used fetal skin fibroblasts from monozygotic (MZ) twins discordant for T21. We also used skin fibroblasts from MZ twins concordant for T21, normal MZ twins without T21, and unrelated normal and T21 individuals. We applied Reduced Representation Bisulfite Sequencing (RRBS) to generate genome wide nucleotide resolution of DNA methylation based on high throughput sequencing between each pair of samples. RRBS revealed differentially methylated promoter regions (DMRs) in MZ twins discordant for T21 that have also been observed in unrelated normal and T21 individuals. The identified DMRs are enriched for genes involved in embryonic organ morphogenesis. These DMRs are maintained in iPS cells generated from this twin pair and are correlated with the gene expression changes. We have also observed an increase in DNA methylation level in the T21 methylome compared to the normal euploid methylome. This observation is concordant with the up regulation of DNA methyltransferase enzymes (DNMT3B and DNMT3L) and down regulation of DNA demethylation enzymes (TET2 and TET3) in the iPSC of the T21 versus normal twin. Additionally, two sets of T21 MZ twins discordant for heart defect highlighted DMRs in genes that are associated with heart development. In conclusion the study of DNA methylation differences in MZ twins discordant for genomic abnormalities is a promising approach to understand the molecular pathophysiology of aneuploidies.

  Study EGAS00001001051

• RB Loss in Resistant EGFR Mutant Lung Adenocarcinomas that Transform to Small Cell Lung Cancer

  Tyrosine kinase inhibitors (TKIs) are effective treatments for non-small cell lung cancers (NSCLCs) with epidermal growth factor receptor (EGFR) mutations. However, they do not lead to cures, and, on average, relapse occurs typically occurs afterafter an average of one year of continuous treatment. In a subset of patients, aA fundamental histological transformation from NSCLC to small cell lung cancer (SCLC) is observed in a subset of the resistant cancers, but the molecular changes associated with this transformation remain unknown. Analysis of a cohort of tumor samples and cell lines derived from resistant EGFR mutant patients with SCLC transformation revealed that RB is lost in 100% of these SCLC transformed cases, but rarely in those that remain NSCLC. Global changes in gene expression, including increased neuroendocrine marker expression and absence of EGFR expression, are observed in resistant EGFR mutant cancers that transformed to SCLC. Consistent Further, with their similarities to classical SCLC at the genetic and gene expression levelincreased neuroendocrine marker and decreased EGFR expression as well as greater sensitivity to BCL2 family inhibition are observed in resistant SCLC transformed cancers compared to resistant NSCLCs. , cell lines derived from resistant EGFR mutant SCLC biopsies are substantially more sensitive to Bcl-2 family inhibition with ABT-263 compared to those derived from resistant EGFR mutant NSCLCs. Together, these findings suggest that despite developing initially as EGFR mutant adenocarcinomas, this subset of resistant cancers ultimately adopt many of the molecular and phenotypic characteristics of classical SCLC.

  Study EGAS00001001102

• High-throughput Detection of Clinically Relevant Mutations in Archived Tumor Samples By Multiplexed PCR and Next Generation Sequencing

  PURPOSE: Tailoring cancer treatment to tumor molecular characteristics promises to make personalized medicine a reality. However, reliable genetic profiling of archived clinical specimens has been hindered by limited sensitivity and high false positive rates. Here, we describe a novel methodology, MMP-seq, which enables sensitive and specific high-throughput, high content genetic profiling in archived clinical samples. EXPERIMENTAL DESIGN: We first validated MMP-seq’s technical performance in 66 cancer cell lines and a Latin square cross-dilution of known somatic mutations. We next characterized the performance of MMP-seq in 17 formalin-fixed paraffin embedded (FFPE) clinical samples using matched fresh frozen (FF) tissue from the same tumors as benchmarks. To demonstrate the potential clinical utility of our methodology, we profiled FFPE tumor samples from 73 endometrial cancer patients. RESULTS: We demonstrated that MMP-seq enabled rapid and simultaneous profiling of a panel of 88 cancer genes in 48 samples, and detected variants at frequencies as low as 0.4%. We identified DNA degradation and deamination as the main error sources and developed practical and robust strategies for mitigating these issues, and dramatically reduced the false positive rate. Applying MMP-seq to a cohort of endometrial tumor samples identified extensive, potentially actionable alterations in the PI3K and RAS pathways, including novel PIK3R1 hotspot mutations that may disrupt negative regulation of PIK3CA. CONCLUSIONS: MMP-seq provides a robust solution for comprehensive, reliable and high- throughput genetic profiling of clinical tumor samples, paving the way for the incorporation of genomic-based testing into clinical investigation and practice.

  Study EGAS00001000674

• Exome Sequencing Of 75 Individuals From Multiply Affected Coeliac Families

  Coeliac disease is a highly heritable common autoimmune disease involving chronic small intestinal inflammation in response to dietary wheat. The human leukocyte antigen (HLA) region, and 40 newer regions identified by genome wide association studies (GWAS) and dense fine mapping, account for ~40% of the disease heritability. We hypothesized that in pedigrees with multiple individuals with coeliac disease rare [minor allele frequency (MAF) &lt;0.5%] mutations of larger effect size (odds ratios of ~ 2 – 5) might exist. We sequenced the exomes of 75 coeliac individuals of European ancestry from 55 multiply affected families. We selected interesting variants and genes for further follow up using a combination of: an assessment of shared variants between related subjects, a model-free linkage test, and gene burden tests for multiple, potentially causal, variants. We next performed highly multiplexed amplicon resequencing of all RefSeq exons from 24 candidate genes selected on the basis of the exome sequencing data in 2,248 unrelated coeliac cases and 2,230 controls. 1,335 variants with a 99.9% genotyping call rate were observed in 4,478 samples, of which 939 were present in coding regions of 24 genes (Ti/Tv 2.99). 91.7% of coding variants were rare (MAF &lt;0.5%) and 60% were novel. Gene burden tests performed on rare functional variants identified no significant associations (P&lt;1x10-3) in the resequenced candidate genes. Our strategy of sequencing multiply affected families with deep follow up of candidate genes has not identified any new coeliac disease risk mutations.

  Study EGAS00001001093

• UK10K NEURO FSZNK

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. This Finnish schizophrenia sample set has been collected from a population cohort using national registers. The entire sample collection consists of 2756 individuals from 458 families of whom 931 are diagnosed with schizophrenia spectrum disorder. Families outside Kuusamo (n=288) all had at least two affected siblings. All diagnoses are based on DSM-IV and for a large fraction of cases there is cognitive data.For further details/descriptions with regard to this data set please contact Tiina Paunio (tiina.paunio@thl.fi)

  Study EGAS00001000119

• UK10K NEURO ASD SKUSE

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches. This sample set of UK origin consists of clinically identified subjects with Autism Spectrum Disorders, mostly without intellectual disability (ie. Verbal IQs &gt;70). The subjects represent children and adults with Autism, Asperger syndrome or Atypical Autism, identified according to standardized research criteria (ADI-algorithm, ADOS). A minority has identified comorbid neurodevelopmental disorders (e.g. ADHD). Family histories are available, with measures of broader phenotype in first-degree relatives.For further information on this cohort please contact David Skuse (d.skuse@ucl.ac.uk).

  Study EGAS00001000114

• UK10K_NEURO_MUIR

  In the UK10K project we propose a series of complementary genetic approaches to find new low frequency/rare variants contributing to disease phenotypes. These will be based on obtaining the genome wide sequence of 4000 samples from the TwinsUK and ALSPAC cohorts (at 6x sequence coverage), and the exome sequence (protein coding regions and related conserved sequence) of 6000 samples selected for extreme phenotypes. Our studies will focus primarily on cardiovascular-related quantitative traits, obesity and related metabolic traits, neurodevelopmental disorders and a limited number of extreme clinical phenotypes that will provide proof-of-concept for future familial trait sequencing. We will analyse directly quantitative traits in the cohorts and the selected traits in the extreme samples, and also use imputation down to 0.1% allele frequency to extend the analyses to further sample sets with genome wide genotype data. In each case we will investigate indels and larger structural variants as well as SNPs, and use statistical methods that combine rare variants in a locus or pathway as well as single-variant approaches.The sample selection consists of subjects with schizophrenia (SZ), autism, or other psychoses all with mental retardation (learning disability). The samples were initially collected under the leadership of Walter J Muir (deceased), now with Prof. Blackwood, Dr McKechanie and Prof McIntosh as custodians. These subjects represent the intersection of severe forms of neurodevelopmental disorders, appear to have a higher rate of familiality of SZ than typical, and are likely to have more serious and penetrant forms of mutations.For further information on this cohort please contact Andrew McIntosh (andrew.mcintosh@ed.ac.uk).

  Study EGAS00001000122

• These are from Korean HCC samples with exome sequencing

  Hepatic resection is the most curative treatment option for early-stage hepatocellular carcinoma, but is associated with a high recurrence rate, which exceeds 50% at 5 years after surgery. Understanding the genetic basis of hepatocellular carcinoma at surgically curable stages may enable the identification of new molecular biomarkers that accurately identify patients in need of additional early therapeutic interventions. Whole exome sequencing and copy number analysis was performed on 231 hepatocellular carcinomas (72% with hepatitis B viral infection) that were classified as early-stage hepatocellular carcinomas, candidates for surgical resection. Recurrent mutations were validated by Sanger sequencing. Unsupervised genomic analyses identified an association between specific genetic aberrations and postoperative clinical outcomes. Recurrent somatic mutations were identified in 9 genes, including TP53, CTNNB1, AXIN1, RPS6KA3, and RB1. Recurrent homozygous deletions in FAM123A, RB1, and CDKN2A, and high-copy amplifications in MYC, RSPO2, CCND1, and FGF19 were detected. Pathway analyses of these genes revealed aberrations in the p53, Wnt, PIK3/Ras, cell cycle, and chromatin remodelling pathways. RB1 mutations were significantly associated with cancer-specific and recurrence-free survival after resection (p = 0.016 and p = 0.001, respectively). FGF19 amplifications, known to activate Wnt signalling, were mutually exclusive with CTNNB1 and AXIN1 mutations, and significantly associated with cirrhosis (p = 0.017). RB1 mutations can be used as a prognostic molecular biomarker for resectable hepatocellular carcinoma. Further study is required to investigate the potential role of FGF19 amplification in driving hepatocarcinogenesis in patients with liver cirrhosis and to investigate the potential of anti-FGF19 treatment in these patients.

  Study EGAS00001000604

• Chromatin immunoprecipitation linked to next-generation whole genome sequencing (ChIP-Seq) for H3K36me3 in paediatric high grade glioma cell lines KKNS4 and SF188 with and without a G34V mutation in H3F3A

  Glioblastomas of children and young adults have a median survival of only 12-15months and are clinically and biologically distinct from histologically similar cancers in older adults1. They are defined by highly specific mutations in the gene encoding the histone H3.3 variant H3F3A2, occurring either at or close to key residues marked by methylation for regulation of transcription – K27 and G34. We performed chromatin immunoprecipitation linked to next-generation whole genome sequencing (ChIP-Seq) for H3K36me3 in order to test the hypothesis that, rather than total H3K36me3, the G34V mutation may instead result in differential binding of the trimethyl mark throughout the genome. Compared with H3F3A wild-type SF188 paediatric GBM cells, H3K36me3 was found to be significantly differentially bound in KNS42 cells at 5130 distinct regions of the genome corresponding to 156 genes. Concurrent whole genome DNA sequencing demonstrated that this was not confounded by copy number alterations. The cerebral hemispheric-specific G34 mutation drives a distinct expression signature through differential genomic binding of the K36 trimethylation mark (H3K36me3). The transcriptional program induced recapitulates that of the developing forebrain, and involves numerous markers of stem cell maintenance, cell fate decisions and self-renewal. Critically, H3F3A G34 mutations cause profound upregulation of MYCN, a potent oncogene which is causative of glioblastomas when expressed in the correct developmental context. This driving aberration is selectively targetable in this patient population by inhibiting kinases responsible for stabilisation of the protein.

  Study EGAS00001001437

• Spatial and temporal evolution of distal 10q deletion, a prognostically unfavorable event in diffuse low-grade gliomas

  Background: The disease course of patients with diffuse low-grade glioma is notoriously unpredictable. Temporal and spatially distinct samples may provide insight into the evolution of clinically relevant copy number aberrations (CNAs). The purpose of this study is to identify CNAs that are indicative of aggressive tumor behaviour and can thereby complement the prognostically favorable 1p/19q co-deletion. Results: Genome-wide, 50 base pair single-end, sequencing was performed to detect CNAs in a clinically well-characterized cohort of 98 formalin-fixed paraffin-embedded low-grade gliomas. CNAs are correlated with overall survival as an endpoint. Seventy-five additional samples from spatially distinct regions and paired recurrent tumors of the discovery cohort were analysed to interrogate the intratumoral heterogeneity and spatial evolution. Loss of 10q25.2-qter is a frequent subclonal event and significantly correlates with an unfavorable prognosis. A significant correlation is furthermore observed in a validation set of 126 and confirmation set of 184 patients. Loss of 10q25.2-qter arises in a longitudinal manner in paired recurrent tumor specimens, whereas the prognostically favorable 1p/ 19q co-deletion is the only CNA that is stable across spatial regions and recurrent tumors. Conclusions: CNAs in low-grade gliomas display extensive intratumoral heterogeneity. Distal loss of 10q is a late onset event and a marker for reduced overall survival in low-grade glioma patients. Intratumoral heterogeneity and higher frequencies of distal 10q loss in recurrences suggest this event is involved in outgrowth to the recurrent tumor.

  Study EGAS00001000643

• Study on the consequences of prenatal famine exposure on DNA methylation.

  The potential of the epigenome to undergo environmentally-induced changes during early development is central to its postulated role in human disease. Here, we report a genome-scale view of DNA methylation differences after early gestational malnutrition that was caused by the Dutch Hunger Winter, a severe famine at the end of World War II. We adopted an annotation guided analysis of reduced-representation bisulfite sequencing data generated in 24 middle-aged individuals who had been exposed to famine in utero and 24 unexposed sibling controls. We observed differential DNA methylation primarily at genomic regions with regulatory potential, including enhancers active during early development. In-depth annotation of 181 prenatally-induced differentially methylated regions (P-DMRs) showed that they were predominantly located in gene bodies and co-occurred with histone marks denoting an active state. In line with the early gestational timing of the exposure, the P-DMRs mapped to genes that were enriched for differential expression patterns during blastocyst development and early organogenesis. Validation and further explorative analyses of 6 P-DMRs mapping to SMAD7, CDH23, INSR, RFTN1, CPT1A and KLF13 highlighted the critical role of gestational timing, indicated that differential methylation extended along pathways involved in growth, development and metabolism, and suggested associations with birth weight and serum LDL cholesterol in adulthood and confirmed an enhancer function for the INSR and CPT1A P-DMRs. Our results point toward a link between prenatal malnutrition and epigenetic modulation of growth and lipid metabolism that may underlie the adverse metabolic phenotype of exposed individuals in later life.

  Study EGAS00001000668

• Dilgom_Exome

  The Finrisk sample sets are part of the National FINRISK Study. It is a large population survey on risk factors of chronic, noncommunicable diseases. The survey is carried out since 1972 every five years using independent, random and representative population samples from different parts of Finland. The main results from the previous FINRISK 2007 survey are published.The National FINRISK Study Survey was carried out in 5 areas in Finland and 2000 inhabitants aged 25-75 years were invited to participate in each year. Among findings were that Finns continue to gain weight.Data from FINRISK surveys are used for many different research projects and for national health monitoring needs. The recent research activities deal, in addition to cardiovascular diseases and the classical risk factors, also with e.g. asthma and allergy, alcohol, socioeconomic factors and genetic epidemiology.The FINRISK study is part of the MORGAM Project (MONICA Risk, Genetics, Archiving and Monigraph), sponsored by the EU and MDECODE (Molecular Diversity and Epidemiology of Common Disease) program coordinated by the University of Michigan.The exome sequencing study will be part of the Dilgom study, which is a part of the larger Finrisk population based health study performed in Finland. It consists of 5000 individuals with a prospective aspect of metabolic traits. The cohort has been extensively phenotyped for their cardiovascular and metabolic status. So far, we have performed a 660K Illumina GWAS and a full genome wide expression study from peripheral blood cells of 500 individuals. The cohort has also been in total genotyped by the cardiometabochip

  Study EGAS00001000086

• Quantification of chromosomal copy number aberrations by shallow whole-genome sequencing

  Detection of DNA copy number aberrations by shallow whole-genome sequencing (WGS) faces many challenges including lack of completion and errors in the human reference genome, repetitive sequences, polymorphisms, variable sample quality, and biases in the sequencing procedures.Formalin-fixed paraffin-embedded (FFPE) archival material, the analysis of which is important for studies of cancer, presents particular analytical difficulties due to degradation of the DNA and frequent lack of matched reference samples. We present a robust, cost-effective WGS method for DNA copy number analysis that addresses these challenges more successfully than currently available procedures. In practice very useful profiles can be obtained with 0.1 fold genome coverage. We improve on previous methods by; first, implementing a combined correction for sequence mappability and GC content, and second, applying this procedure to sequence data from the 1000 Genomes Project in order to develop a blacklist of problematic genome regions. A small subset of these blacklisted regions were previously identified by ENCODE, but the vast majority are novel unappreciated problematic regions. Our procedures are implemented in a pipeline called QDNAseq. We have analyzed over 1,000 samples, most of which were obtained from the fixed tissue archives of over 25 institutions.We demonstrate that for most samples our sequencing and analysis procedures yield genome profiles with noise levels near the statistical limit imposed by read counting. The described procedures also provide better correction of artifacts introduced by low DNA quality than prior approaches, and better copy number data than high-resolution microarrays at substantially lower cost.

  Study EGAS00001000642

• Prediction and quantification of splice events from RNA-seq data

  Analysis of splice variants from short read RNA-seq data remains a challenging problem. Here we present a novel method for the genome-guided prediction and quantification of splice events from RNA-seq data, which enables the analysis of unannotated and complex splice events. Splice junctions and exons are predicted from reads mapped to a reference genome and are assembled into a genome-wide splice graph. Splice events are identified recursively from the graph and are quantified locally based on reads extending across the start or end of each splice variant. We assess prediction accuracy based on simulated and real RNA-seq data, and illustrate how different read aligners (GSNAP, HISAT2, STAR, TopHat2) affect prediction results. We validate our approach for quantification based on simulated data, and compare local estimates of relative splice variant usage with those from other methods (MISO, Cufflinks) based on simulated and real RNA-seq data. In a proof-of-concept study of splice variants in 16 normal human tissues (Illumina Body Map 2.0) we identify 249 internal exons that belong to known genes but are not related to annotated exons. Using independent RNA samples from 14 matched normal human tissues, we validate 9/9 of these exons by RT-PCR and 216/249 by paired-end RNA-seq (2 x 250 bp). These results indicate that de novo prediction of splice variants remains beneficial even in well-studied systems. An implementation of our method is freely available as an R/Bioconductor package SGSeq.

  Study EGAS00001001026