13380 results for ""60+mice"+"6+weeks"+diet+study+2014+2017+DOI"

in 50.01 milliseconds.

• Somatic mutation and selection at epidemiological scale - TwinsUK_ExomeNanoSeq_Buccal

  Exome NanoSeq in buccal swab samples obtained from patients with a wide range of clinical phenotypes. As we age, many tissues become colonised by microscopic clones carrying somatic driver mutations. Some of these clones represent a first step towards cancer whereas others may contribute to ageing and other diseases. However, our understanding of the clonal landscapes of human tissues, and their impact on cancer risk, ageing and disease, remains limited due to the challenge of detecting somatic mutations present in small numbers of cells. Here, we introduce a new version of nanorate sequencing (NanoSeq), a duplex sequencing method with error rates of less than 5 per billion base pairs, which is compatible with whole-exome and targeted gene sequencing. Deep sequencing of polyclonal samples with single-molecule sensitivity enables the simultaneous detection of mutations in large numbers of clones, yielding accurate somatic mutation rates, mutational signatures and driver mutation frequencies in any tissue. Applying targeted NanoSeq to 1,042 non-invasive samples of oral epithelium and 371 samples of blood from a twin cohort, we found an unprecedentedly rich landscape of selection, with 46 genes under positive selection driving clonal expansions in the oral epithelium, over 62,000 driver mutations, and evidence of negative selection in some genes. The high number of positively selected mutations in multiple genes provides high-resolution maps of selection across coding and non-coding sites, a form of in vivo saturation mutagenesis. Multivariate regression models enable mutational epidemiology studies on how carcinogenic exposures and cancer risk factors, such as age, tobacco or alcohol, alter the acquisition and selection of somatic mutations. Accurate single-molecule sequencing has the potential to unveil the polyclonal landscape of any tissue, providing a powerful tool to study early carcinogenesis, cancer prevention and the role of somatic mutations in ageing and disease.

  Dataset EGAD00001015620

• Somatic mutation and selection at epidemiological scale - TwinsUK_TargetedEMSeq_Blood

  Development of a targeted methylation assay to determine the cell-type composition of a sample. As we age, many tissues become colonised by microscopic clones carrying somatic driver mutations. Some of these clones represent a first step towards cancer whereas others may contribute to ageing and other diseases. However, our understanding of the clonal landscapes of human tissues, and their impact on cancer risk, ageing and disease, remains limited due to the challenge of detecting somatic mutations present in small numbers of cells. Here, we introduce a new version of nanorate sequencing (NanoSeq), a duplex sequencing method with error rates of less than 5 per billion base pairs, which is compatible with whole-exome and targeted gene sequencing. Deep sequencing of polyclonal samples with single-molecule sensitivity enables the simultaneous detection of mutations in large numbers of clones, yielding accurate somatic mutation rates, mutational signatures and driver mutation frequencies in any tissue. Applying targeted NanoSeq to 1,042 non-invasive samples of oral epithelium and 371 samples of blood from a twin cohort, we found an unprecedentedly rich landscape of selection, with 46 genes under positive selection driving clonal expansions in the oral epithelium, over 62,000 driver mutations, and evidence of negative selection in some genes. The high number of positively selected mutations in multiple genes provides high-resolution maps of selection across coding and non-coding sites, a form of in vivo saturation mutagenesis. Multivariate regression models enable mutational epidemiology studies on how carcinogenic exposures and cancer risk factors, such as age, tobacco or alcohol, alter the acquisition and selection of somatic mutations. Accurate single-molecule sequencing has the potential to unveil the polyclonal landscape of any tissue, providing a powerful tool to study early carcinogenesis, cancer prevention and the role of somatic mutations in ageing and disease.

  Dataset EGAD00001015623

• Ancestry and somatic profile predict acral melanoma origin and prognosis – WXS

  Acral melanoma, which is not ultraviolet (UV)-associated, is the most common type of melanoma in several low- and middle-income countries including Mexico. Latin American samples are significantly underrepresented in global cancer genomics studies, which directly affects patients in these regions as it is known that cancer risk and incidence may be influenced by ancestry and environmental exposures. To address this, we characterise the genome and transcriptome of 123 acral melanoma tumours from 92 Mexican patients, a population notable because of its genetic admixture. Compared with other studies of melanoma, we found fewer frequent mutations in classical driver genes such as BRAF, NRAS or NF1. While most patients had predominantly Amerindian genetic ancestry, those with higher European ancestry had increased frequency of BRAF mutations and a lower median number of structural variants. The tumours with activating BRAF mutations have a transcriptional profile more similar to cutaneous non-volar melanocytes, suggesting that acral melanomas in these patients may arise from a distinct cell of origin compared to other tumours arising in these locations. KIT mutations were found in a subset of these tumours, and quadruple wild-type samples (non BRAF/NRAS/NF1/KIT) differed from mutated samples in their structural genomic profile and overall and recurrence-free survival patterns. Transcriptional profiling defined three expression clusters; these characteristics were associated with recurrence-free and overall survival. We highlight potential novel low-frequency drivers, such as PTPRJ, NF2 and RDH5. Our study enhances knowledge of this understudied disease and underscores the importance of including samples from diverse ancestries in cancer genomics studies.

  Dataset EGAD00001015755

• Phylogenetic reconstruction of adult blood cancer reveals early origins and lifelong evolution - WGS

  Mutations in cancer-associated genes drive tumour outgrowth. However, the timing of driver mutations and dynamics of clonal expansion that lead to human cancers are largely unknown. We used 580,133 somatic mutations from whole-genome sequencing of 1013 clonal haematopoietic colonies to reconstruct the phylogeny of haematopoiesis, from embryogenesis to clinical disease, in 12 patients with myeloproliferative neoplasms which are blood cancers more common in older age. JAK2V617F, the pathognomonic mutation driving the majority of these cancers, was acquired in utero or childhood, with upper estimates of age of acquisition from 33 weeks gestation to 10.8 years, in all 5 patients in whom JAK2V617F was either the only or the first driver event.  Driver mutations associated with age-related clonal haematopoiesis occurred prior to or following JAK2V617F,  as independent clonal expansions in JAK2V617F-mutated patients, and as large clonal expansions in JAK2V617F-unmutated patients . These mutations were also acquired in utero or childhood, with DNMT3A mutations occurring by 8 weeks of gestation to 7.6 years across 4 patients, and PPM1D mutation occurring by age 5.8yrs in a patient with MPN lacking phenotypic driver mutations. Sequential driver mutation acquisition was common, separated by decades across life, and often outcompeted ancestral clones. The mean latency between JAK2V617F acquisition and clinical presentation was 30 years (range 11-54 years). Rates of clonal expansion were inferred from phylogenetic trees and varied substantially (3% to 190% expansion/year), were affected by additional driver mutations, and were predictive of latency to clinical presentation. Driver mutations and rates of expansion would have been detectable in blood one to four decades before clinical presentation. This study reveals how driver mutation acquisition early in life with life-long growth and evolution underlie adult myeloproliferative neoplasms, providing opportunities for early detection and intervention and a new paradigm for cancer development.

  Dataset EGAD00001007714

• Immune variation leads to diverse outcomes in human malaria (2020-01-15)

  Falciparum malaria is clinically heterogeneous and yet in most cases the risk of life-threatening disease dramatically declines after the first few infections of life because children rapidly acquire disease tolerance (resistance to severe malaria without improved control of parasite burden). Identifying the factors that determine clinical outcome in a malaria-naive host is therefore paramount to reduce malaria mortality. However, the relative contribution of disease-causing variants of the Plasmodium var gene family versus pathogenic inflammatory cytokine cascades remains fiercely debated - we sought to reconcile these conflicting arguments by studying their interaction in vivo. To this end, two human challenge models were used to reveal the parasite-host interactions that underpin variation in falciparum malaria. To capture the diversity of human immune responses, each individual was analysed independently by tracking dynamic changes in their whole blood transcriptome through time. And to uncover evidence of preferential expansion of disease-causing variants, var gene expression was tracked in vivo from the start to end of infection. In this way, we could show that group A var genes are always expressed upon liver egress but in a minority population that does not increase over 10-days of blood cycling; there is no selection of disease-causing variants in the naive host. In fact, parasites do not respond in any way to differences or changes in host environment. On the other hand, host-intrinsic variation determines the intensity of inflammation and progression to clinical malaria. And furthermore, regulation of the interferon signaling network controls host fate. These data emphasise the role of human immune decision-making in shaping course & outcome of infection. This data is part of a pre-publication release. For information on the proper use of pre-publication data shared by the Wellcome Trust Sanger Institute (including details of any publication moratoria), please see http://www.sanger.ac.uk/datasharing/ . This dataset contains all the data available for this study on 2020-01-15.

  Dataset EGAD00001005790

• Somatic mutation and selection at epidemiological scale - TwinsUK_TargetedEMSeq_Buccal

  Development of a targeted methylation assay to determine the cell-type composition of a sample. As we age, many tissues become colonised by microscopic clones carrying somatic driver mutations. Some of these clones represent a first step towards cancer whereas others may contribute to ageing and other diseases. However, our understanding of the clonal landscapes of human tissues, and their impact on cancer risk, ageing and disease, remains limited due to the challenge of detecting somatic mutations present in small numbers of cells. Here, we introduce a new version of nanorate sequencing (NanoSeq), a duplex sequencing method with error rates of less than 5 per billion base pairs, which is compatible with whole-exome and targeted gene sequencing. Deep sequencing of polyclonal samples with single-molecule sensitivity enables the simultaneous detection of mutations in large numbers of clones, yielding accurate somatic mutation rates, mutational signatures and driver mutation frequencies in any tissue. Applying targeted NanoSeq to 1,042 non-invasive samples of oral epithelium and 371 samples of blood from a twin cohort, we found an unprecedentedly rich landscape of selection, with 46 genes under positive selection driving clonal expansions in the oral epithelium, over 62,000 driver mutations, and evidence of negative selection in some genes. The high number of positively selected mutations in multiple genes provides high-resolution maps of selection across coding and non-coding sites, a form of in vivo saturation mutagenesis. Multivariate regression models enable mutational epidemiology studies on how carcinogenic exposures and cancer risk factors, such as age, tobacco or alcohol, alter the acquisition and selection of somatic mutations. Accurate single-molecule sequencing has the potential to unveil the polyclonal landscape of any tissue, providing a powerful tool to study early carcinogenesis, cancer prevention and the role of somatic mutations in ageing and disease.

  Dataset EGAD00001015622

• 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

• Prevention and Early Treatment of Acute Lung Injury Network - Reevaluation of Systemic Early Neuromuscular Blockade (PETAL ROSE-BioLINCC)

  Data Access NOTE: Please refer to the "Authorized Access" section below regarding accessing data through the BioData Catalyst ecosystem. The data from this accession is not available for download through dbGaP. Available Data: The data available for request now include Long Term Outcome data.Biospecimens: Access to Biospecimens is through the NHLBI Biologic Specimen and Data Repository Information Coordinating Center (BioLINCC). Biospecimens from PETAL ROSE include Plasma, DNA, Whole Blood, and Urine. Please note that use of biospecimens in genetic research is subject to a tiered consent.Objectives: To determine the efficacy and safety of early neuromuscular blockade with concomitant heavy sedation as compared with a strategy of usual care with lighter sedation targets in patients with moderate-to-severe ARDS.Background: It has been well established that the approaches used for the application of mechanical ventilation in patients with acute respiratory distress syndrome (ARDS) can affect survival and outcomes after discharge from the intensive care unit (ICU). A large, multicenter trial conducted a decade before this study reported that the early administration of a 48-hour infusion of neuromuscular blockade in patients with moderate-to-severe ARDS (defined by a ratio of the partial pressure of arterial oxygen [Pao2] to the fraction of inspired oxygen [Fio2] of Participants: There were 1006 participants.Design: PETAL-ROSE was a multicenter, unblinded, randomized trial of patients with moderate-to-severe ARDS. Participants were randomly assigned in a 1:1 ratio to receive 48 hours of continuous neuromuscular blockade with concomitant deep sedation (intervention group) or to receive usual care without routine neuromuscular blockade and with lighter sedation targets (control group).Patients in the intervention group who were not under deep sedation at baseline were deeply sedated within 4 hours after randomization. Subsequently, patients in this group received an intravenous bolus of 15 mg of cisatracurium, followed by a continuous infusion of 37.5 mg per hour for 48 hours. After the 48-hour trial intervention period, decisions regarding further use of neuromuscular blockade, including the choice of agent, were left to the discretion of the treating clinician. Neuromuscular blockade could be stopped early if the patient met the criteria for freedom from mechanical ventilation (Fio2 ≤0.40 and PEEP ≤8 cm of water) for at least 12 hours. All patients were treated with a strategy of low tidal volume ventilation within 2 hours after randomization and a high PEEP strategy for up to 5 days after randomization. Assessors who were unaware of the group assignment interviewed surviving patients or their proxies at 3, 6, and 12 months after randomization. The primary end point was in-hospital death from any cause at 90 days (in-hospital was defined as the time in the trial hospital plus transfer to another hospital, including the time in long-term acute care facilities). Conclusions: After the second interim analysis, the decision to stop the trial for futility was made independently by the data and safety monitoring board.In critically ill patients identified shortly after the diagnosis of moderate-to-severe ARDS, the addition of early continuous neuromuscular blockade with concomitant deep sedation did not result in lower mortality than a usual-care approach to mechanical ventilation that included lighter sedation targets.

  Study phs003878

• Submitter Portal

  Submitter Portal The Submitter Portal does not support array submissions! Please refer to the array documentationAny new objects registered in the Submitter Portal will not be available for use in an array submission For further information please check our Submission FAQs, submission quickguide as well as submission terms! The metadata submission process can be difficult and time-consuming. For this reason, the EGA has developed the Submitter Portal, a tool that was created to offer a simplified and user-friendly method of registering metadata. Our portal provides features that are intended to make the input of metadata easier, ensuring that your data is registered correctly and effectively. Our page is divided into logical sections and includes a helpful video instruction to further assist you as you complete the submission process. With the Submitter Portal, you can rest assured that the submission of your data is in good hands. Submitter Portal Index Previous steps Access to Submitter Portal Start a Submission Take The Tour Tutorial Previous steps Create your EGA account To submit data, you first need to create your EGA user. Then, once your account has been verified by the Help Desk team, you will be able to request a submitter role and send the EGA Data Processing Agreement (DPA). In the case you already have an EGA account, or an ega-box (submission account), you can skip this step. Upload your files Please note that all your files must be encrypted using the Crypt4GH tool before upload them.As soon as you are assigned with a submitter role and added your public key to your profile, you will be able to connect to the EGA inbox and upload your files. Understand the EGA metadata schema It's crucial to comprehend the EGA metadata schema, a set of rules that specify how data is organised, described, and shared inside the EGA, in order to get the most out of this resource. Learn all about EGA metadata schema! Register your DAC and policy There are two objects in the EGA metadata schema that are registered in a separate portal. All DACs and policy objects are registered using the DAC Portal, a tool developed by EGA to help data controllers manage their data stored at the EGA. You can find relevant information in the DAC Portal Guide. Access to Submitter Portal 1. Request Submitter Role Once you have created your EGA user, request a Submitter Role and complete relevant information such as submission type and size. 2. Signature of the Data Processing Agreement The Data Processing Agreement (DPA) signed by CRG and EMBL-EBI (Joint Data Processors for the EGA) can be downloaded to the submitter for its signature. Further information related to the DPA can be in the Privacy Notice section. The DPA should be signed by a legal representative of the submitter with power of attorney/appointment duly authorised to sign contracts on behalf of the institution. The submitter will return a fully-executed copy of the DPA to us. Please note that due to the high volume of submissions to EGA the DPA is non-negotiable. Once you have completed all the information and attached the DPA signed, send your Submitter request and the EGA Helpdesk Team will validate it and send a confirmation to start submitting your metadata. Start a Submission The steps below will show you how to start a submission and are accompanied with a tutorial and the Take the tour of Submitter Portal. 1. Login Go to the Submitter Portal interface and log in using your EGA account credentials. 2. Create a submission The Submitter Portal interface will show you any already existing open submission. To start a new submission click on “Create a Submission”.In this phase you have the option to add any collaborator(s) as well as their permissions (“Read only” or “Read & Write”). If you want to invite all the collaborators from a previous submission, you can add them by selecting said submission.In the event that the submission has been closed, you will need to reopen it. To do this, you will need to modify one of the objects in the submission, for example, copy and paste the title. Once the submission is open, you can add a collaborator. Further information can be found in the Submission Tour 3. Create Study Click on add Study and complete the form with the relevant information. The fields with an asterisk (*) are mandatory. Further information can be found in the Submission Tour4. Create Sample Sample registration can be done either individually or in batch by using a template. Batch upload: download the template and complete it with information about the sequencing samples. Upload the samples.csv and the sample will be created.Add single sample: complete the mandatory fields with the metadata information of your sample.  Further information can be found in the Submission Tour5. Create Experiment Register an experiment with ​information about the sequencing methods, protocols and machines. Experiments generate the connection between samples and study. Further information can be found in the Submission Tour6. Create Run Samples, experiments and files are linked through runs. To register a run, you will be asked to select the appropriate experiment and the file format. Run registration can be done either individually or in batch by using a template. Batch upload: download the template and complete it with the sample - file linkage. Upload the run.csv and the runs will register.As the INBOX is a directory, bear in mind to include the file path (/). Example: /filename.c4gh. Add single run: select the file format, sample and experiment from the drop box. For re-using submitted samples, please use the EGA accession ID (EGANXXX). Please note that files must be uploaded in your inbox before linking your files and samples. Further information can be found in the Submission Tour7. Create Analysis Analysis should be only used for BAM/BAI or CRAM/CRAI pair, VCF and phenotype linkage to samples. The analysis is an EGA specific metadata object that links samples to files. This object also stores some metadata about your experiments, such as the experiment type, genome reference, or the platform used. **If only BAM or CRAM alignment files are submitted but not the original unaligned FASTQ files, then please make sure that the BAM or CRAM files also contain the unaligned reads. This is critical to enable primary re-analysis and re-alignment of the dataset using new tools or future genome assemblies.** Start by selecting the samples to be linked to the file, and then populate the required attribute fields. Please be aware of the existence of mandatory fields that must be populated.When populating the chromosome field (mandatory), please select the chromosomes from the checkbox. Please take in consideration that EGA allows for the re-use of any registered metadata. Therefore, a previously registered study or samples can be referenced for the analysis in the current data submission.Further information can be found in the Submission Tour 8. Create Dataset Dataset contains the collection of runs/analysis data files subjected to controlled access by being associated with a Policy. Complete the mandatory fields such as title, description, type and policy. Please, take in consideration that DAC and Policy must be registered in advance in the DAC Portal. Similar to other metadata objects, you can re-use previous runs or analyses, using their EGA accessions IDs (runs - EGAR, analyses - EGAZ). Further information can be found in the Submission Tour 9. Finalise submission Once the files are ingested please delete them from your SFTP INBOX before finalising your submission. This will prevent the files from reappearing in your Submitter Portal Inbox. Once you have completed your submission, click on the Finalise button, and schedule an estimated release date. Afterwards, your submission will be sent to Helpdesk for approval. Please note that once the Helpdesk has approved your submission, you will receive the identifiers for your data. See the Frequently Asked Questions (FAQ) for more information. 10. Update closed submission If you want to append more file to a closed submission, you will need to re-open the submission where the dataset was registered. Find the Submission: go to the Submitter Portal Menu -> Datasets section and find the submission id (EGA0XXX).  Reopen the submission where the dataset is located: to reopen a submission, you must edit any object within the submission. For example, you can copy and paste the title. Take The Tour The Submitter Portal Take The Tour is a visual guidance of all the steps described previously about how to start a submission. Tutorial In this tutorial we will show you how to use the Submitter Portal to register your metadata. Please note the video focuses on the run-based submission (for raw files - fastq - and aligned data - BAM/CRAM) and analysis-based submission (for your BAM/BAI pairs, variation -VCF - and phenotype files).

  Documentation submission/metadata/submission/sequencing-phenotype/submitter-portal

• 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