1905 results

in 22.31 milliseconds.

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

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

  Study EGAS00001005679

• Defining the metastasome in colorectal cancer: Implications for hypotheses on metastasis evolution and personalized therapy (HIPO-032)

  Personalized cancer therapy aims at individual genetic changes characterizing primary cancers. Still, however, metastasis is the major clinically unmet problem, also due to an incomplete understanding of its molecular evolution.This study is the first to define the metastasis-specific whole genome landscape of human colorectal primary vs. matched metastatic lesions.Protein coding, ncRNA genes and 3’UTRs harbored about a third each of single nucleotide variations (SNVs)/indels, 19% of them being metastasis-specific. We found novel mutational hills among ncRNAs and 3’UTRs, copy number changes able to explain changes in microRNA expression in metastases, and novel metastasis-specific mutations in the 3’UTRs on important cancer signalling genes. Some metastatic lesions showed targetable mutations not detected in the primaries. Analysis of mutual exclusivity patterns supported a model that has progressive alterations in important colorectal cancer genes and revealed new partners. Metastatic lesions were enriched in SNVs, specific structural variations, and alterations in genes affecting cell adhesion and extracellular matrix interaction.Furthermore, a hepatic stellate activation cascade was enriched in metastases, suggesting genetic programs for site-specific colonization. Allele frequency distribution and mutational signature analysis suggests a common ancestor clone to both the primary tumor and the metastasis, with additional mutagenesis ongoing in both sites independently. Taken together, our results significantly add to hypothesis generation on metastasis evolution, and suggest novel metastasis-specific genomic changes which would be missed by current personalized therapy concepts. (HIPO-032)

  Study EGAS00001002717

• Partner-independent fusion gene detection by multiplexed CRISPR/Cas9 enrichment and long-read Nanopore sequencing

  Fusion genes are hallmarks of various cancer types and important determinants for diagnosis, prognosis and treatment possibilities. The promiscuity of fusion genes with respect to partner choice and exact breakpoint-positions restricts their detection in the diagnostic setting, even for known and recurrent fusion gene configurations. To accurately identify these gene fusions in an unbiased manner, we developed FUDGE: a FUsion gene Detection assay from Gene Enrichment. FUDGE couples target-selected and strand-specific CRISPR/Cas9 activity for enrichment and detection of fusion gene drivers (e.g. BRAF, EWSR1, KMT2A/MLL) - without prior knowledge of fusion partner or breakpoint-location - to long-read Nanopore sequencing. FUDGE encompasses a dedicated bioinformatics approach (NanoFG) to detect fusion genes from Nanopore sequencing data. Our strategy is flexible with respect to target choice and enables multiplexed enrichment for simultaneous analysis of several genes in multiple samples in a single sequencing run. We observe on average a 508-fold on-target enrichment and identify fusion breakpoints at nucleotide resolution - all within two days. We demonstrate that FUDGE effectively identifies fusion genes in cancer cell lines, tumor samples and on whole genome amplified DNA irrespective of partner gene or breakpoint-position in 100% of cases. Furthermore, we show that FUDGE is superior to routine diagnostic methods for fusion gene detection. In summary, we have developed a rapid and versatile fusion gene detection assay, providing an unparalleled opportunity for pan-cancer detection of fusion genes in routine diagnostics.

  Study EGAS00001003964

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

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

  Study EGAS00001005677

• Exploration of CNV’s and SNV’s in cancers with well-known genetic rearrangements: Identification of additional genetic changes in rearrangements-driven cancer

  Specific cancers are well known to have specific gene rearrangements. Representative cancer is chronic myeloid leukemia (CML) which is famous for BCR/ABL gene rearrangements. By the way, also genetic rearrangements are not infrequently found in acute leukemia, and those include CBFB/MYH11 rearrangements, PML/RARA rearrangements, AML/ETO rearrangements, and BCR/ABL rearrangements. The situation is same is for lung adenocarcinoma of non-smokers and prostate cancers, and ALK and RET gene rearrangements are recently found. However, not all these rearrangements are sufficient to develop cancer, and it is suggested that some genetic aberrations are required for the development and progression of cancer in addition to these genetic rearrangements. In this study, we have a plan to explore the existence of additional genetic changes including single nucleotide variants (SNV’s) and copy number variations (CNV’s) in cancers with specific genetic rearrangements so as to define additional genetic changes that confer full oncogenic potential to cancer cells. So as to look at SNV’s and small indels, we are going to use whole exome sequencing (WES) data. And as for CNV’s, we are going to use array-CGH data. Also, to check the functionality of these genetic changes, we are going to use whole transcriptome data (WTS). In addition, in samples with whole genome sequencing (WGS ) data, we are going to explore genetic changes in intronic regions also. In fact, we sequenced 10 AML samples which are known to have specific genetic rearrangements. Analysis of these WGS, WES and WTS data are under the way. If we have a chance to access WES, array-CGH and WTS data of other cancer samples with specific genetic rearrangements, we would be able to compare additional genetic changes required for oncogenesis in cancer cells with various genetic rearrangements. We hope this result will enlighten our understanding of cancer genetics further.

  Study EGAS00001000673

• Breast_Heterogeneity_Validation

  Validation of subs and indels called in breast whole genome studies

  Study EGAS00001001972

• Recurrent DNMT3B gene rearrangements are associated with unfavorable outcome in dicentric (9;20)-positive pediatric BCP-ALL

  The dicentric chromosome dic(9;20)(p11∼13;q11), occurs in 2% of pediatric acute lymphoblastic leukemia (ALL) cases. However, the prognostic significance of the dic(9;20) alterations and the mechanism in which these alterations drive leukemogenesis remains elusive. Our study describes the demographic, clinical, prognostic, and molecular characteristics of the dic(9;20)-positive ALL in a cohort of 57 pediatric BCP-ALL patients. Targeted RNA-seq and Array CGH analysis unraveled heterogeneous breakpoints on chromosomes 9 and 20 and frequent deletions of the IKZF1 gene. Together with the intrinsic susceptibility for deletions in CDKN2A, CDKN2B, and PAX5 genes in the dic(9;20)-positive ALL, these patients are eligible for treatment escalation in the AIEOP-BFM ALL 2017 clinical study. Survival analysis revealed poor 5-year event-free survival (pEFS) of 69% in dic(9;20)-positive ALL cases, compared to BCP-ALL patients enrolled in modern treatment studies. Strikingly, we identified six dic(9;20)-positive cases with rearrangements involving the DNMT3B gene, and poor 5-year pEFS of 25% (SE = 20%), compared to 79% (SE = 9%) for the remaining cohort (P = 0.011). Despite the loss of the methyltransferase domain of the DNMT3B gene in most cases with DNMT3B rearrangements, we did not observe major changes in the global methylation patterns. Differential gene expression analysis of DNMT3B-rearranged ALL unraveled a dysregulation of genes involved in hematopoiesis, lymphocyte maturation, and the PI3K/AKT signaling pathway. After confirmation in an independent cohort, DNMT3B-rearranged ALL cases should be considered high-risk for relapse and treated accordingly. We provide clinical data sets of Array CGH, targeted RNA-seq, total RNA-seq, whole genome bisulfite (WGBS) and whole genome DNA sequencing (WGS) obtained from bone marrow (BM) or peripheral blood (PB) mononuclear cells of 57 pediatric patients with dicentric chromosome dic(9;20) positive Acute lymphocytic leukemia (ALL), from which in 6 cases DNMT3B gene rearrangement was identified. This data is complemented by total RNA-seq and WGBS of samples from 4 additional ALL patients with a t(12;21) translocation and ETV6-RUNX1 gene fusion. DNA was isolated from BM or PB lymphocytes using Qiagen QIAamp DNA Blood Midi Kit to perform i) Array CGH of 58 dic(9;20) positive samples by hybridizing 500ng DNA using a Agilent 400K SurePrint G3 Custom CGH Human Genome Microarray (e-Array design 84704) ii) WGBS of 6 DNMT3B rearrangement positive samples and 4 ETV6-RUNX1 positive samples using Tecan TrueMethyl oxBS-Seq module for library preparation and Illumina NovaSeq 6000 platform to run 2x151 cycles iii) WGS of DNMT3B rearrangement positive samples using Illumina Lotus DNA Library Prep Kit followed by sequencing running 2x160 cycles on an Illumina NovaSeq 6000 platform. RNA was isolated from PB lymphocytes using the PerkinElmer Chemagic 360 instrument, followed by i) targeted RNA-seq of 56 dic(9;20) positive samples prepared using Illumina TruSight RNA Pan-Cancer Panel and sequenced on an Illumina MiSeq platform running 2x75 cycles ii) total RNA-seq of 6 DNMT3B rearrangement positive samples and 4 ETV6-RUNX1 positive samples utilizing TruSeq Stranded Total RNA Library Prep Gold kit and running 2x100 cycles on an Illumina NovaSeq 6000 platform.

  Study EGAS00001007383

• Genomic consequences of aberrant DNA repair mechanisms stratify ovarian cancer histotypes

  Aberrant DNA repair processes are a hallmark of human tumours. How these deficiencies variously impact the genomes of ovarian cancer between and within histotypes remains unknown. We studied the whole genome pointmutation and structural variation patterns of 133 tumours (59 high grade serous (HGSC), 35 clear cell (CCOC), 29 endometrioid (ENOC), and 10 adult granulosa cell tumours (GCT)) as a substrate for class discovery in ovarian cancer. Ab-initio clustering of integrated point mutation and structural variation signatures revealed seven novel subgroups, comprising between and within histotype stratification. Prevalence of fold-back inversions (FBI) co-localised with high level amplification events divided HGSC into two prognostically significant subgroups. This finding was recapitulated in two additional independent cohorts (total n=576 cases), and transcended gene based mutation status and gene expression as prognostically relevant features of HGSC. CCOC were divided into one group harbouring a signature reflective of active genome editing via APOBEC enzymes (26%) and a complementary group with an age-relatedmutational signature. ENOC were divided intomicrosatellite instable (MSI) cases (28%) with a distinctmismatch repairmutation signature and an outlying mutation rate, with the remainder of cases distributed amongst the remaining six groups. We suggest the FBI, MSI, and APOBEC groups represent biologic strata that could direct treatment strategies. Taken together, our work establishes the efficacy of the somatic genome as a biomarker to stratify ovarian cancers, simultaneously identifying patients that may benefit from emerging therapeutics and distinct subgroups within classic ovarian cancer histotypes.

  Study EGAS00001002390

• Analysis of the Elements Involved in the Enrichment of a Panel of Genomic Regions by Nanopore Sequencing Using Adaptive Sampling

  DNA sequencing using nanopore technologies with the affordable MinION device is useful for the identification and characterization of structural variants, long haplotypes, sequencing of repetitive regions and identification of epigenetic modifications. The main limitation of this approach is the low coverage obtained, which might be avoided by adaptive sampling, a computationally controlled method of enrichment for targeted genomic regions. This study dissects the factors involved in the enrichment by adaptive sampling of a panel of 18 human genome regions containing 20 genes implicated in breast cancer sequenced in 16 patients with familial breast cancer negative for NGS screening. An average coverage of 2.0x was obtained for the whole genome and 5.1x for selected regions. Sequencing time was the main factor improving coverage. The selection of long reads (>1 Kb) did not improve the enrichment. The length of the selected region, which in our study ranged from 126 to 565 Kb, did not play a significant role in enrichment. However, the region containing PMS2 showed significantly lower coverage, which could be explained by the high number of PMS2 pseudogenes (N = 14), which were also enriched. Our study shows new evidence of enrichment obtained by adaptive sampling in a panel of genomic regions and shows parameters, the relevance of sequencing time and the role of pseudogenes, that improve the enrichment yield with no library reloading or GPU use, data useful for a more efficient application of this procedure in future studies.

  Study EGAS00001007375

• Clinical relevance of somatic copy-number alterations in plasma circulating tumor DNA from advanced EGFR-mutated NSCLC patients treated with osimertinib

  Background: Somatic copy-number alterations (SCNAs) are associated with drug resistance in patients with advanced epidermal growth factor receptor (EGFR)-mutated non-small cell lung cancer (NSCLC) treated with EGFR-tyrosine kinase inhibitors (TKIs). The purpose of our study was to assess the clinical relevance and development of genome-wide SCNAs in plasma circulating tumor DNA (ctDNA) under osimertinib therapy.Methods: We included 43 patients with advanced EGFR T790M-positive lung adenocarcinoma who were treated with osimertinib after progression under previous EGFR-TKI therapy. We performed genomic profiling of ctDNA in plasma samples from each patient obtained pre-osimertinib and after patients developed resistance to osimertinib. SCNAs were detected by shallow whole-genome plasma sequencing and EGFR mutations were assessed by droplet digital PCR.Results: SCNAs in resistance-related genes (rrSCNAs) were detected in 10 out of 39 (25.6%) evaluable patients before start of osimertinib. The presence of rrSCNAs in plasma before initiation of osimertinib therapy was associated with a lower response rate to osimertinib (50% versus 86%, p=0.02) and was an independent predictor for shorter progression-free survival (adjusted HR 3.26, 95% CI 1.37-7.75, p=0.008) and overall survival (adjusted HR 2.41, 95% CI 1.09-5.34, p=0.03).Conclusions: Genomic profiling of plasma ctDNA is clinically relevant and affects the efficacy and clinical outcome of osimertinib. Our approach allows comprehensive assessment of SCNAs in plasma samples of lung adenocarcinoma patients and may help to guide genotype-specific therapeutic strategies in the future.

  Study EGAS00001004539

• A whole genome analysis of single fetal human stem cells from the liver and the intestine

  Mutations of embryonic and fetal origin have the potential to affect a large proportion of adult cells and may alter cancer predisposition or lead to genetic disease syndromes. We have recently shown that human adult-stem cells progressively acquire approximately 40 novel tissue-specific mutations per year throughout postnatal life. Prenatal mutation rates are as yet unknown. Here we determined genome-wide mutation patterns of single stem cells in human development by sequencing of clonally expanded intestinal and liver organoid cultures of 2nd trimester human foetuses. Our results show that mutation rates in fetal stem cells are significantly higher than in adult stem cells.

  Study EGAS00001002886

• A Comprehensive Catalogue of Somatic Mutations from a Human Cancer Genome

  We have carried out complete sequencing of the genome of the human male maligant melanoma cell line COLO-829 using the Illumina Genome Analyzer II. We generated a sequencing library with a median insert size of ~200 bp following random fragmentation and gel fractionation of the genomic DNA. We sequenced 75 bases from both ends of these templates to cover the COLO-829 genome to an average depth of more than 40x. We have carried out purity-filtering (PF) to remove mixed reads, where two or more different template molecules are close enough on the surface of the flow-cell to form a mixed or overlapping cluster. No other filtering of the data has been carried out prior to submission. We have also submitted sequence data for a lymphoblastoid control cell line COLO-829BL from the same individual.

  Study EGAS00000000052

• Massive Genomic Rearrangment Acquired in a Single Catastrophic Event During Cancer Development

  Cancer is driven by mutation. Using Agilent exome hybridisation capture and Illumina GA massively parallel sequencing technology, we aim to sequence ~1600 microRNAs plus the protein coding genome of 25 matched human renal cancer samples. Bespoke algorithms are being developed to identify the somatically acquired point mutations, insertions and deletions in these samples. This project will give unprecedented insights into mutational processes, cellular repair pathways and gene networks associated with renal cancer development.Agilent whole exome hybridisation capture will be performed on genomic DNA derived from 25 renal cancers and matched normal DNA from the same patients. Three lanes of Illumina GA sequencing will be performed on the resulting 50 exome libraries and mapped to build 37 of the human reference genome to facilitate the identification of novel cancer genes.

  Study EGAS00000000029

• G&T-seq: parallel sequencing of the genomes and transcriptomes of single cells

  The ability to simultaneously sequence the genome and transcriptome of the same single cell offersa powerful means to dissect functional genetic heterogeneity at the cellular level. Here we describeG&T-seq, a method for separating and sequencing genomic DNA and full-length mRNA from singlecells. By applying G&T-seq to over 220 single cells we reveal cellular properties that cannot beinferred from DNA or RNA sequencing alone, including associations between DNA copy numbervariation and gene expression dosage. We further demonstrate the detection of coding interchromosomalfusions and single nucleotide variants in both the genomes and transcriptomes ofindividual cells. G&T-seq enables the study of genotype-phenotype associations in single cells, andthe investigation of DNA cell lineage trees of healthy and diseased tissues with transcriptome inferredcell types and states.

  Study EGAS00001001204

• An integrative model of pathway convergence in genetically heterogeneous blast crisis chronic myeloid leukemia (CML)

  Targeted therapies against the BCR-ABL1 kinase have revolutionized treatment of chronic phase (CP) chronic myeloid leukemia (CML). In contrast, management of blast crisis (BC) CML remains challenging because BC cells acquire complex molecular alterations that confer stemness features to progenitor populations and resistance to BCR-ABL1 tyrosine kinase inhibitors. Comprehensive models of BC transformation have proved elusive because of the rarity and genetic heterogeneity of BC, but are important for developing biomarkers predicting BC progression and effective therapies. To better understand BC, we performed an integrated multiomics analysis of 74 CP and BC samples using whole-genome and exome sequencing, transcriptome and methylome profiling, and chromatin immunoprecipitation followed by high-throughput sequencing. Employing pathway-based analysis, we found the BC genome was significantly enriched for mutations affecting components of the polycomb repressive complex (PRC) pathway. While transcriptomically, BC progenitors were enriched and depleted for PRC1- and PRC2-related gene sets respectively. By integrating our data sets, we determined that BC progenitors undergo PRC-driven epigenetic reprogramming toward a convergent transcriptomic state. Specifically, PRC2 directs BC DNA hypermethylation, which in turn silences key genes involved in myeloid differentiation and tumor suppressor function via so-called epigenetic switching, whereas PRC1 represses an overlapping and distinct set of genes, including novel BC tumor suppressors. On the basis of these observations, we developed an integrated model of BC that facilitated the identification of combinatorial therapies capable of reversing BC reprogramming (decitabine+PRC1 inhibitors), novel PRC-silenced tumor suppressor genes (NR4A2), and gene expression signatures predictive of disease progression and drug resistance in CP.

  Study EGAS00001001751

• Identification of the dismal subtype of B-ALL with dysregulation of CDX2 and UBTF

  PURPOSE: Advance in the knowledge of genomic basis of B-cell acute lymphoblastic leukemia (B-ALL) has changed the treatment approaches and diagnostic assays and improved their outcome. Despite the recent development of next generation sequencing approaches, there are still some cases that their drivers are unknown. We aimed to define and describe the characteristics of an additional B-ALL subtype by integrating genomic, transcriptomic, and epigenomic approach. METHODS: More than 2,000 B-ALL cases of RNA-seq data from from Eastern Cooperative Oncology Group (ECOG) and the American College of Radiology Imaging Network (ACRIN) (n=764), the Children's Oncology Group (COG), Munich Leukemia Laboratory (MLL) (n=280), and other several collaborators were analyzed to identify the new subtype. Whole genome sequence was performed to detect mutations, stractural variants, and copy number alterations. The tridimentional analysis of chromatin was performed with HiChIP. RESULTS: We identified the new high-risk B-ALL subtype “CDX2/UBTF” that exhibited unique gene expression profiles, which is enriched in female (77.3%), adolescents and young adults (52.6%), and relapse cohort (3.3%). CD10 negativity and IgM positivity are hallmarks of this subtype with NTRK3 expression that can be an effective marker at diagnostic approach and a potential targeted therapy with TRK and mulkinase inbitors, larotrectinib and entrectinib. Genetically, two different alterations define CDX2/UBTF. One is UBTF-ATXN7L3 fusions caused by deletion of 17q21.31, and the other is extopic CDX2 expression through the enhancer hijacking mechanism induced by deletion of 13q12.2. Other genomic features include gain of 1q and PAX5 rearrangement (PAX5-ZCCHC7), might induce upregulation of histone cluster genes and PAX5. CONCLUSION: We described the novel B-ALL subtype “CDX2/UBTF” that has unique clinical and genomic characteristics.

  Study EGAS00001005863

• Distinct evolutionary trajectories of primary high grade serous ovarian cancers revealed through spatial mutational profiling

  High-grade serous ovarian cancer (HGSC) is characterized by poor outcome, often attributed to emergence of treatment-resistant sub-clones. We sought to measure the degree of genomic diversity within primary, untreated HGSC to examine the natural state of tumor evolution prior to therapy. We performed exome sequencing, copy number analysis, targeted amplicon deep sequencing and gene expression profiling on thirty-one spatially and temporally separated HGSC tumor specimens (six patients) including ovarian masses, distant metastases, and fallopian tube lesions. We found widespread intra-tumoral variation in mutation, copy number, and gene expression profiles, with key driver alterations in genes present in only a subset of samples (e.g. PIK3CA, CTNNB1, NF1). On average, only 51.5% of mutations were present in every sample of a given case (range: 10.2% to 91.4%), with TP53 as the only somatic mutation consistently present in all samples. Complex segmental aneuploidies, such as whole genome doubling, were present in a subset of samples from the same individual, with divergent copy number changes segregating independently of point mutation acquisition. Reconstruction of evolutionary histories showed one patient with mixed HGSC and endometrioid histology with common etiologic origin in the fallopian tube and subsequent selection of different driver mutations in the histologically distinct samples. In this patient, we observed mixed cell populations in the early fallopian tube lesion, indicating diversity arises at early stages of tumorigenesis. Our results reveal that HGSC exhibit highly individual evolutionary trajectories and diverse genomic tapestries prior to therapy, exposing an essential biological characteristic to inform future design of personalized therapeutic solutions and investigation of drug resistance mechanisms.

  Study EGAS00001000547

• Whole_genome_sequencing_of_100_DDD_trios_with_suspected_noncoding_causal_mutations

  We would like to whole-genome sequence 100 trios (300 individuals) from the DDD project in which we suspect the proband has a pathogenic noncoding mutation, to 30X coverage. These have been carefully selected to maximise the chances of finding pathogenic noncoding mutations. Specifically, we have selected probands who have no likely causal mutation within the coding region but who have a heterozygous loss-of-function mutation inherited from one parent in a known recessive developmental disorder gene, and whose phenotype is a good match for the typical phenotype caused by recessive variants in that gene. We suspect that these probands may have a heterozygous noncoding mutation on the other allele which reduces expression of the gene so that, together, these two mutations are effectively equivalent to a homozygous loss-of-function mutation. About a quarter of the probands have affected siblings and, if those siblings have been exome sequenced, we have checked that they are identical-by-descent in the region of the gene. We propose to search in the promoter region, intronic and annotated regulatory regions around these genes of interest for rare variants inherited from the other parent (i.e. not the one who transmitted the loss-of-function variant), or de novo variants. If we find such variants, and demonstrate an enrichment of this compound heterozygous LoF noncoding mechanism above that expected by chance, we intend to follow up specific candidates by RNAseq and potentially CRISPR/Cas9 experiments. 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 please see http://www.sanger.ac.uk/datasharing/

  Study EGAS00001002452

• Genomic Profiling of Pediatric Tumors by Cell Free DNA Sequencing

  BACKGROUND Analysis of cell-free DNA (cfDNA) from peripheral blood has been studied extensively for detection and monitoring therapy response in adult cancers, but utility in pediatric tumors is unclear. METHODS We studied 233 children with hematologic, solid and brain tumors. cfDNA was extracted from 1 mL of plasma collected at diagnosis. Circulating tumor DNA (ctDNA) in cfDNA was quantified and diverse classes of somatic genomic variants were detected using PeCan-Seq. PeCan-Seq findings were compared with matched tumor whole-genome and whole-exome sequencing data. Minimal residual disease (MRD) status was evaluated by cfDNA PeCan-Seq in 14 patients with B-cell acute lymphoblastic leukemia. RESULTS Plasma cfDNA yields at diagnosis were higher for hematologic malignancies than for solid and brain tumors. Patients with hematologic malignancy were 100% ctDNA-positive, the median ctDNA fraction in cfDNA was 0.77. We detected 97% variants that were present in tumors and covered by the capture panel, and 50 subclonal variants that were present in cfDNA only. Our approach enabled detection of sequence and structural variants in B-ALL ctDNA at varying levels of MRD. Nineteen of 38 patients with solid tumors were ctDNA-positive, with high prevalence in Ewing Sarcoma, liver tumors and osteosarcoma. One of 19 brain tumor patients was ct-DNA positive, indicating lack of utility for plasma-derived cfDNA in this disease. Notably, cfDNA analysis of a single plasma sample identified genomic features in ctDNA derived from B-ALL, concurrent primary neuroblastoma, and a later neuroblastoma relapse, indicating that genomic variations from multiple tumors could be detected simultaneously. CONCLUSIONS Pediatric patients with hematologic malignancy had an abundance of plasma ctDNA at diagnosis. Plasma ctDNA levels varied between different solid tumor types, more studies are needed to determine utility. Panel-based cfDNA sequencing could serve as a non-invasive approach for pediatric cancer diagnosis and disease monitoring.

  Study EGAS50000000072

• One-step generation of tumor models by base editor multiplexing in adult stem cell-derived organoids

  Recent optimization of CRISPR/Cas9-mediated genome engineering has resulted in the development of base editors that can efficiently mediate C>T and A>G transitions. Combining these genome engineering tools with human adult stem cell (ASC)-derived organoid technology holds promise for disease modeling. Here, we demonstrate the application of base editors for the generation of complex tumor models in human ASC-derived hepatocyte, endometrial and intestinal organoids. First, using conventional and evolved Cas9-variants, we show efficacy of both cytosine and adenine base editors and use them to model four hot-spot point mutations in CTNNB1 in hepatocyte organoids. Next, we apply C>T base editors in endometrial organoids to insert nonsense mutations in PTEN and demonstrate tumorigenicity even in the heterozygous state. Furthermore, we use cytosine base editors for simultaneous oncogene activation (PIK3CA) and tumor-suppressor inactivation (APC and TP53). To increase the flexibility of base editor multiplexing, we then combine SpCas9 and SaCas9 base editors for simultaneous C>T and A>G editing at individual target sites. Finally, we show the power of base editor multiplexing by modeling colorectal tumorigenesis in a single step by simultaneously transfecting sgRNA’s targeting four cancer genes. Seven clonal organoid lines and one bulk wild-type control sample were paired-end whole-genome sequenced using the Illumina Novaseq 6000 system. We sequenced four clonal intestinal organoid lines harbouring engineered TP53 and FBXW7 mutations as well as three lines targeted for oncogenic APC/TP53/PIK3CA/SMAD4 mutations. This WGS showed, as previously reported, a genome-wide increase in C>T mutations due to C>T base editor off-target activity, which is not enriched in predicted off-target regions based on the sgRNA sequences. Furthermore, we confirmed the absence of editing-induced driver mutations and lack of off-target mutational hotspots created by the genomic engineering.

  Study EGAS00001006886

• Whole genome characterisation of lung cancer organoids and tissue

  Organoids are self-organizing 3D structures grown from stem cells that recapitulate essential aspects of organ structure and function. Here we describe a method to establish long-term expanding human airway organoids from broncho-alveolar biopsies or lavage material. The pseudostratified airway organoid epithelium consists of basal cells, functional multi-ciliated cells, mucus-producing goblet cells, and CC10-secreting club cells. Airway organoids derived from cystic fibrosis (CF) patients allow assessment of CFTR function in an organoid swelling assay. Organoid culture conditions also allow gene editing as well as the derivation of various types of lung cancer organoids. Respiratory syncytial virus (RSV) infection recapitulated central disease features and dramatically increases organoid cell motility, found to be driven by the non-structural viral NS2 protein. We conclude that human airway organoids represent versatile models for the in vitro study of hereditary, malignant, and infectious pulmonary disease.

  Study EGAS00001002899

• Epigenomic data of HEMa_LP

  To profile the epigenomic landscape of Human Epidermal Melanocytes , adult, lightly pigmented(HEMa-LP), whole genome bisulphite sequencing (WGBS) and ChIP-seq of histone marks data were generated.

  Study EGAS00001004016

• A small cell lung cancer genome reports complex tobacco exposure signatures

  Cancer is driven by mutation. Worldwide, tobacco smoking is the major lifestyle exposure that causes cancer, exerting carcinogenicity through 60 chemicals that bind and mutate DNA. Using massively parallel sequencing technology, we sequenced a small cell lung cancer cell line, NCI-H209, to explore the mutational burden associated with tobacco smoking. 22,910 somatic substitutions were identified, including 132 in coding exons. Multiple mutation signatures testify to the cocktail of carcinogens in tobacco smoke and their proclivities for particular bases and surrounding sequence context. Effects of transcription-coupled repair and a second, more general expression-linked repair pathway were evident. We identified a tandem duplication that duplicates exons 3-8 of CHD7 in-frame, and another two lines carrying PVT1-CHD7 fusion genes, suggesting that CHD7 may be recurrently rearranged in this disease. These findings illustrate the potential for next-generation sequencing to provide unprecedented insights into mutational processes, cellular repair pathways and gene networks associated with cancer.

  Study EGAS00000000051

• Genome sequencing identifies splice-disrupting variants in childhood heart disease

  We used genome sequencing (GS) to identify cardiac specific splice disrupting variants that were subsequently validated in a subset of participants. In turn, we developed a heart-specific model for canonical and non-canonical splice variants, which was applied to patients with congenital heart disease (CHD). These included patients with two of the most common forms of cyanotic CHD i.e. tetralogy of Fallot (TOF) and dextro-transposition of the great arteries (TGA). In addition to identifying canonical splice-disrupting variants in known CHD-related genes in 1% cases, this approach identified putatively damaging non-canonical splice-disrupting variants in 11% of isolated CHD, with deeply intronic variants representing 53% of non-canonical splice-disrupting variants in CHD genes. GS was critical for the identification of variants that would not be captured by routine clinical genetic tests including exome sequencing, while cardiac RNA-Seq allowed for high specificity in the interpretation of splice-disrupting effects in the heart.

  Study EGAS50000000586

• Multiplexing cortical brain organoids for the longitudinal dissection of developmental traits at single cell resolution

  Dissecting human neurobiology at high resolution and with mechanistic precision requires a major leap in scalability, given the need for experimental designs to multiple individuals and, prospectively, population cohorts. To lay the foundation for this, we have implemented and benchmarked complementary strategies to multiplex brain organoids by pooling cells from different pluripotent stem cell lines (PSC) either during organoid generation (mosaic models) or before single cell-RNAseq library preparation (downstream multiplexing). We have also developed a new computational method, SCanSNP, and a consensus call to deconvolve cell identities, overcoming current criticalities in doublets prediction and low quality cells identification. We validated both multiplexing methods for charting neurodevelopmental trajectories at high resolution, thus linking specificity to genetic variation. Finally, we modelled their scalability across different multiplexing combinations and showed that mosaic organoids represent an enabling method for high-throughput settings. Together, this multiplexing suite of experimental and computational methods provides a highly scalable resource for disease modelling.

  Study EGAS50000000698