22928 results for "ega"

in 69.07 milliseconds.

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

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

  Dataset EGAD00001000666

• RBMX functional retrocopy safeguards brain development

  Retrotransposition has generated thousands of intronless gene copies in mammalian genomes, yet their contribution to brain development and evolution remains largely unexplored. We investigated the role for RBMX retrocopy in shaping neurodevelopment and modulating disease. First, we reported 8 distinct pathogenic variants in RBMX, an X-linked splicing regulator, in 11 males presenting with intellectual disability, microcephaly, and cortical malformations. Then, we combined transcriptomic profiling, protein–protein and protein–RNA interaction studies both in human cellular models and mouse embryonic cortices to assess the functional redundancy between RBMX and its retrocopy RBMXL1. Finally, we use mouse genetics to dissect RBMX function and its compensation by RBMXL1 in corticogenesis. We show that RBMX pathogenic variants disrupt cortical development through both loss- and gain-of-function mechanisms. Surprisingly, despite severe phenotypes in humans, Rbmx-deficient mice display only mild cortical abnormalities, a discrepancy likely due to compensation by Rbmxl1, a retrocopy that arose independently in mice and humans. We demonstrate that RBMX and RBMXL1 share protein and RNA partners and act redundantly in brain development, with RBMXL1 buffering the impact of RBMX deficiency. Together, these findings establish that RBMXL1 functions as a functional paralog of RBMX safeguarding neurodevelopment by buffering deleterious variation. More broadly, these results identify retrocopies as active contributors to neurodevelopmental robustness and suggest that functional retrocopies may have facilitated the evolutionary diversification of the mammalian brain.

  Study EGAS50000001650

• Prognostic relevance of microenvironmental factors CD163 and CD8 combined with EZH2 and chromosome 18 gain in a validation cohort of follicular lymphoma patients of the Lunenburg Lymphoma Biomarker Consortium

  Background: In follicular lymphoma (FL), studies addressing the prognostic value of microenvironment-related immunohistochemical (IHC) markers and tumor cell-related genetic markers have yielded conflicting results, precluding implementation in practice. Therefore, the Lunenburg Lymphoma Biomarker Consortium (LLBC) performed a validation study for published markers. Methods: To maximize sensitivity, an end-of-spectrum design was applied for 122 uniformly immunochemotherapy-treated FL patients retrieved from international trials and registries; early failure (EF): progression or lymphoma-related death <2 years versus long remission: response duration of >5 years. IHC staining for T-cells and macrophages was performed on tissue microarrays from initial biopsy and scored with a validated computer-assisted protocol. Shallow whole-genome and deep targeted sequencing was performed on the same samples. Results: 96/122 cases with complete molecular and immunohistochemical data were included in the analysis. EZH2 wild-type (p=0.006), gain of chromosome 18 (p=0.002), low percentages of CD8+ cells (p=0.011) and CD163+ areas (p=0.038) were associated with EF. No significant differences in other markers were observed, thereby refuting previous claims on their prognostic significance. Conclusion: Using an optimized study design, this LLBC study validates wild-type EZH2 status, gain of chromosome 18, low percentages of CD8+ cells and CD163+ area as predictors of EF to immunochemotherapy in FL.

  Study EGAS00001002049

• The double-hit signature identifies double-hit diffuse large B-cell lymphoma with genetic events cryptic to FISH​

  High-grade B-cell lymphomas with MYC and BCL2 and/or BCL6 rearrangements (HGBL-DH/THs) include a group of diffuse large B-cell lymphomas (DLBCLs) with inferior outcomes after standard chemoimmunotherapy. We recently described a gene expression signature that identifies 27% of germinal center B-cell DLBCLs (GCB-DLBCLs) as having a double-hit-like expression pattern (DHITsig) and inferior outcomes; however, only half of these cases have both MYC and BCL2 translocations identifiable using standard breakapart fluorescence in situ hybridization (FISH). Here, 20 DHITsig+ GCB-DLBCLs apparently lacking MYC and/or BCL2 rearrangements underwent whole-genome sequencing. This revealed 6 tumors with MYC or BCL2 rearrangements that were cryptic to breakapart FISH. Copy-number analysis identified 3 tumors with MYC and 6 tumors with MIR17HG gains or amplifications, both of which may contribute to dysregulation of MYC and its downstream pathways. Focal deletions of the PVT1 promoter were observed exclusively among DHITsig+ tumors lacking MYC translocations; this may also contribute to MYC overexpression. These results highlight that FISH fails to identify all HGBL-DH/THs, while revealing a range of other genetic mechanisms potentially underlying MYC dysregulation in DHITsig+ DLBCL, suggesting that gene expression profiling is more sensitive for identifying the biology underlying poor outcomes in GCB-DLBCL.

  Study EGAS00001004285

• Bulk-tissue RNA-sequencing paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from post-mortem neuropathologically-confirmed control individuals.

  Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

  Study EGAS00001006380

• The Dutch Microbiome Project Batch 1

  The Dutch Microbiome Project (DMP) data includes shotgun metagenomic sequencing of faecal samples 8,208 Dutch individuals. Paired-end sequencing was performed using Illumina HiSeq 2000 platform. Data is archived in two batches to facilitate easier data access and upload to EGA. Batch 1 of DMP includes 4396 samples.

  Dataset EGAD00001007027

• The patterns and dynamics of genomic instability in metastatic pancreatic cancer

  Pancreatic cancer is an aggressive malignancy with a five-year mortality of 97–98%, usually due to widespread metastatic disease. Previous studies indicate that this disease has a complex genomic landscape, with frequent copy number changes and point mutations, but genomic rearrangements have not been characterized in detail. Despite the clinical importance of metastasis, there remain fundamental questions about the clonal structures of metastatic tumours, including phylogenetic relationships among metastases, the scale of ongoing parallel evolution in metastatic and primary sites, and how the tumour disseminates. Here we harness advances in DNA sequencing to annotate genomic rearrangements in 13 patients with pancreatic cancer and explore clonal relationships among metastases. We find that pancreatic cancer acquires rearrangements indicative of telomere dysfunction and abnormal cell-cycle control, namely dysregulated G1-to-S-phase transition with intact G2–M checkpoint. These initiate amplification of cancer genes and occur predominantly in early cancer development rather than the later stages of the disease. Genomic instability frequently persists after cancer dissemination, resulting in ongoing, parallel and even convergent evolution among different metastases. We find evidence that there is genetic heterogeneity among metastasis-initiating cells, that seeding metastasis may require driver mutations beyond those required for primary tumours, and that phylogenetic trees across metastases show organ-specific branches. These data attest to the richness of genetic variation in cancer, brought about by the tandem forces of genomic instability and evolutionary selection.

  Study EGAS00000000064

• Nanopore sequencing enables allelic phasing of FLG loss-of-function variants, intragenic copy number variation and methylation status in atopic dermatitis and ichthyosis vulgaris

  Loss-of-function (LoF) variants in the FLG gene are associated with ichthyosis vulgaris (IV) and atopic dermatitis (AD). IV and AD patients with two FLG LoF variants are reported to have a more severe disease course. However, subsets of compound heterozygous patients display phenotypic heterogeneity despite a prediction of a severe phenotype from their genotype. The underlying genetic mechanism of this phenotypic variability is unknown. One possibility is that two LoF FLG variants could be located in cis (on the same allele), thus altering the predicted gene dosage. To investigate the allelic contribution of compound heterozygous FLG LoF mutations to IV phenotypic variability, we employ long reads from Oxford Nanopore Technologies, utilising two methodologies to phase the entire coding region of FLG into parental alleles: an amplicon-sequencing approach and an adaptive-sampling strategy. We investigated the allelic phasing of FLG LoF variants and intragenic-repeat copy-number variations (CNV) of 21 subjects with known compound heterozygous genotypes and mixed IV severity. Our results suggest that IV severity is not linked to the presence of cis variants in our IV patient cohort. Nanopore reads enabled detection and accurate phasing of CNVs and LoF variants into parental alleles, providing a framework for future in-depth genetic analysis. These results highlight the usefulness of utilising long reads for interrogating the FLG gene and other highly repetitive genes for genetic stratification related to clinical phenotypes.

  Study EGAS50000000166

• The tissue origin of highly elevated cell-free DNA in patients with and without cancer

  Forty years ago, the first measurements of cell free DNA (cfDNA) from cancer patients revealed that its total concentration was often elevated compared to healthy controls. The basis for this increase, and in particular the source of the extra cfDNA in such patients, has never been determined. To address this issue, we assessed methylation patterns of cfDNA from 57 individuals, 34 of whom had unequivocally elevated cfDNA concentrations (i.e., >15 ng/mL of plasma). We found that the excess cfDNA in cancer patients did not come from either neoplastic cells or surrounding cells from the tumor's tissue of origin. Instead, 80% of the cfDNA was derived from neutrophils, B cells or T cells, with neutrophils accounting for ~2/3. B cells contributed more than T cells despite the great excess of T cells over B cells in the blood. To see if this result was generalizable to other clinical situations in which elevated cfDNA concentrations have been described, we evaluated nine patients before and one day following surgery. After surgical resection, cfDNA concentration increased by an average of 16.1-fold, as expected. The majority of this increase (55%) was contributed by leukocytes, though damaged hepatocytes contributed a substantial minority (31% of increase). These data suggest that the elevated concentrations of cfDNA in cancer patients, as with patients after surgery, arises largely from dying leukocytes—not from the neoplastic cells or the tissue of origin of the cancer.

  Study EGAS00001005400

• METABRIC miRNA landscape

  MicroRNAs (miRNAs) are key players in cancer pathogenesis1. They display differential expression across breast cancer subtypes, and play oncogenic and tumor-suppressive roles2-6. A systematic study of the global miRNA landscape in breast cancer and the factors shaping it has not been undertaken. Here we report the results of profiling miRNA expression in 1,302 breast tumors with detailed clinical annotation and long-term follow-up, and for which matching genomic and mRNA expression data were available7. This provides a comprehensive view of the quantity, distribution and variation of the miRNA population and dissects how much genomic, transcriptional and post-transcriptional events contribute to miRNA expression architecture. The key clinical parameters and cellular pathways related to the miRNA landscape are identified, exemplified by context-dependent interaction with cell adhesion and Wnt signaling. Strikingly, in breast cancer miRNAs appear to act as modulators of mRNA-mRNA interactions rather than molecular switches. We demonstrate an important modulatory role for miRNAs in the biology of breast tumors devoid of somatic copy number aberrations, a common subtype where immune response is prominent. Remarkably, prognostic miRNA signatures derived in this subtype are consistently prognostic across several other subtypes and validate in external cohorts. These findings represent a new framework to study the biology of miRNA in human cancer.

  Study EGAS00000000122

• Machine learning to detect the SINEs of cancer

  We previously described an approach called RealSeqS to evaluate aneuploidy in plasma cell-free DNA (cfDNA) through the amplification of ~350,000 repeated elements with a single primer. We hypothesized that an unbiased evaluation of the large amount of sequencing data obtained with RealSeqS might reveal other differences between plasma samples from patients with and without cancer. This hypothesis was tested through the development of a novel machine-learning approach called Alu Profile Learning Using Sequencing (A-PLUS) and its application to samples from 5108 individuals, 2037 with cancer and the remainder without cancer. Samples from cancer patients and controls were pre-specified into four cohorts used for: 1) model training, 2) analyte integration and threshold determination, 3) validation, and 4) reproducibility. A-PLUS alone provided a sensitivity of 40.5% across 11 different cancer types in the Validation Cohort, at a specificity of 98.5%. Combining A-PLUS with aneuploidy and 8 common protein biomarkers detected 51% of 1167 cancers at 98.9% specificity. We found that part of the power of A-PLUS could be ascribed to a single feature – the global reduction of AluS sub-family elements in the circulating DNA of cancer patients. We confirmed this reduction through the analysis of another independent dataset obtained with a very different approach (whole genome sequencing). The evaluation of Alu elements therefore has the potential to enhance the performance of several methods designed for the earlier detection of cancer.

  Study EGAS00001007169

• Neoadjuvant immune checkpoint blockade in women with mismatch repair deficient endometrial cancer

  Neoadjuvant immune checkpoint blockade (ICB) has shown unprecedented activity in mismatch repair deficient (MMRd) colorectal cancers, but its effectiveness in MMRd endometrial cancer (EC) remains unknown. In this investigator-driven, phase I, feasibility study (NCT04262089), 10 women with MMRd EC of any grade, planned for primary surgery, received two cycles of neoadjuvant pembrolizumab (200 mg IV) every three weeks. A pathologic response (primary objective) was observed in 5/10 patients, with 2 patients showing a major pathologic response. No patient achieved a complete pathologic response. A partial radiologic response (secondary objective) was observed in 3/10 patients, 5/10 patients had stable disease and 2/10 patients were non-evaluable on Magnetic Resonance Imaging (MRI). All patients completed treatment without severe toxicity (exploratory objective). At median duration of follow-up of 22.5 months, two non-responders experienced disease recurrence. In-depth analysis of the loco-regional and systemic immune response (predefined exploratory objective) showed that monoclonal T cell expansion significantly correlated with treatment response. Tumour-draining lymph nodes displayed clonal overlap with intra-tumoural T cell expansion. All pre-specified endpoints, efficacy in terms of pathologic response as primary endpoint, radiologic response as secondary outcome and safety and tolerability as exploratory endpoint, were reached. Neoadjuvant ICB with pembrolizumab proved safe and induced pathologic, radiologic, and immunologic responses in MMRd EC, warranting further exploration of extended neoadjuvant treatment.

  Study EGAS50000000483

• Biopsy-Derived Organoids in Personalised Early Breast Cancer Care: Challenges of Tumour Purity and Normal Cell Overgrowth Cap Their Practical Utility

  The ability to establish organoids composed exclusively of tumour rather than healthy cells is essential for their implementation into clinical practice. Organoids have in recent years emerged as a powerful tool to expand patient material in culture and generate modifiable 3D models derived from humans or animal models. For translational research, they enable the creation of model systems for an ever-increasing number of cell types and diseases. And in personalised medicine, they potentially allow for functional drug testing with high predictive power in certain settings. We found that using biopsy material from untreated, early-stage primary breast cancer patients poses significant challenges for consistently culturing tumour cells as organoids. Specifically, we observed frequent outgrowth of genetically normal, non-cancerous epithelial cells. We analysed >100 biopsy samples from early-stage breast cancer and present our large collection of >70 organoid lines. We also show methods of assessing successful tumour cell culture in a time and cost-efficient manner, proving a high rate (>85%) of normal cell overgrowth in early-stage breast cancer organoids. Finally, we show a number of successful attempts to culture cancer organoids from mastectomy-derived tissue of advanced, metastatic breast cancer. We conclude that the usefulness of organoids from early breast cancer for translational research and personalised medicine, especially guidance of adjuvant or post-surgical maintenance therapy, is strongly limited by the low success rate of culturing cancerous cells under organoid conditions.

  Study EGAS50000000605

• Single cell transcriptomics in expanded Tregs of APS-1 patients

  Aim of study: Exploring single cell transcriptomics in expanded Tregs of with autoimmune polyendocrine syndrome type 1 (APS-1) patients (global gene expression, immune panel, TCR). N= 9 APS-1 patients and 9 healthy controls Protocol 1. Regulatory T cells (Tregs) were isolated from EDTA-blood, activated with anti-CD3/CD28 and IL2, expanded in cell culture for 14 days and stored at -150℃ in AB serum or FBS + 5% DMSO until use. 2. Cells were thawed, run through a Live/Dead column and carefully counted. 10000 Tregs were used as input for the 10x protocol. 3. cDNA was generated and amplified using Single Cell VDJ 5’ Gel Beads, Chromium Next GEM Chip K Single cell Kit, Next GEM Single Cell 5’ GEM Kit v2 and 5’ Feature Barcode Kit (all from 10x genomics). 4. Gene expression (GEX) libraries were generated by using the Library Construction Kit from 10x genomics. 5. T-cell receptor (TCR) libraries were generated using the Single Cell Human TCR Amplification Kit and the Library Construction Kit from 10x genomics. 6. 10x Genomics Target Hybridization Kit and Human Immunology Panel was used with GEX libraries. Only 8 patient samples and 8 control samples were used for the immune panel. 7. All libraries have a unique sample index (Dual Index TT Set A). 8. All samples were sequenced on Illumina NovaSeq 6000 on a NovaSeq S2 flow cell.

  Study EGAS50000000181

• Multi-omics bulk and single-cell profiling of epithelioid sarcoma

  Epithelioid sarcoma (EpS) is an aggressive sarcoma, characterized by the loss of SMARCB1 expression. EpS is traditionally classified as distal or proximal according to clinicopathological features, but its molecular characteristics remain largely unknown. To establish an EpS molecular classification and uncover determinants of inter- and intra-patient heterogeneity in EpS, we used multi-omics profiling and integrated the genomic, transcriptional and methylome landscapes with single-cell RNA sequencing on fresh samples as well as spatial transcriptomics. We identified two molecular subtypes of EpS: “distal-like” and “proximal-like”, which were distinct from the histological subtypes. Distal-like tumors include all morphological distal subtypes and a subset of proximal and hybrid EpS. They express a specific molecular single-cell derived epithelial-to-mesenchymal transition signature, which associates with improved patient survival and includes Desmoglein 2 (DSG2), which we identify as a potential routine immunohistochemical biomarker to improve diagnosis. Distal-like EpS also display increased peri-tumoral CD8+ T cell infiltrates and specific tumor-immune cell interactions. Conversely, proximal-like EpS harbor a higher inter-tumoral molecular heterogeneity with some cases resembling other SMARCB1-deficient tumors by DNA methylation profiling, and increased intra-tumoral pro-tumoral macrophage infiltrates. Our study introduces a novel molecular classification of EpS with prognostic value, associated diagnostic biomarkers and distinct tumor immune microenvironments, which extends beyond traditional clinicopathological classification. This paves the way for precision medicine-based therapeutic strategies for patients with EpS.

  Study EGAS50000000973

• Decoding Human Endogenous Retrovirus Expression in Liver Metastatic Colorectal Cancers: Implications for Diagnosis and Prognosis

  The association between Human Endogenous RetroVirus (HERV) expression and cancer can be considered to both reflect the epigenetic deregulation of cancer and facilitate oncogenesis through its promoter/enhancer action on genes. More than 50% of colorectal cancers develop liver metastases and HERV has never been studied in this setting. We collected 400 RNAseq samples from more than 200 patients with primary and liver metastases from published and unpublished databases. We corroborated our findings with tandem mass spectrometry from primary colorectal cancer and single-cell RNAseq from liver metastatic colorectal cancer (lmCRC). We identified global stability in HERV expression between liver metastases and primary colorectal cancers, suggesting an early oncogenic footprint. Therefore, we proposed a list of seventeen HERV loci for lmCRC characterization with ERV-derived peptides but no expression in normal tissues. Among them, we discovered six prognostic markers after surgical resection of lmCRC. Four of them, HERVH_Xp22.32a, HERVH_20p11.23b, HERVH_13q33.3, HERVH_13q31.3, showed poor prognostic value when expressed (log-rank p-value 0.028, 0.0083, 9e-4, 0.05, respectively). HERVH_Xp22.2c (log-rank p-value 0.032) and HERVH_8q21.3b (only in multivariate models) were associated with better prognosis when expressed. These markers showed a cumulative effect on survival when expressed. Some HERVs were associated with decreased cytotoxic immune cells and most of them correlated with cell cycle pathways. Our findings provide new insights into lmCRC classification and suggest novel prognostic biomarkers and potential therapeutic targets.

  Study EGAS50000000307

• The Dynamic Immune Behavior of Primary and Metastatic Ovarian Carcinoma

  Patients with high-grade serous ovarian carcinoma are usually diagnosed at an advanced stage and often develop resistance to conventional therapy. Epithelial ovarian tumors have been shown to be rich in CD8+ lymphocytes, which infiltrate both primary and omental sites. In recent years, several prominent studies have highlighted a positive association between the presence of these immune cells and a good prognosis; however, clinical trials of immune checkpoint inhibitors have so far shown unsatisfactory results in ovarian cancer. In this study, we used an innovative combination of single-cell sequencing and spatial transcriptomics to identify the molecular mechanisms that lead to immunosuppression in primary and metastatic HGSC. Primary tumors consistently showed a more active immune microenvironment than did omental tumors. In addition, while several immune cell subtypes were present in all samples, their relative abundance differed between primary and omental tumors. Finally, we found that primary tumors were mostly populated by dysfunctional CD4 and CD8 T cells in later stages of differentiation, while omental tumors were mostly populated by naïve immune cells; this, in turn, was correlated with expression changes in the interferon α and γ pathways in epithelial cells, showing that cross-communication between the epithelial and immune compartments is important for immune suppression in ovarian cancer. These findings can impact the timing and choice of immune-modulating drugs and should be considered when planning treatment for patients with high-grade serous ovarian carcinoma.

  Study EGAS50000000038

• Cancer-associated fibroblasts promote drug resistance in ALK-driven lung adenocarcinoma cells by upregulating lipid biosynthesis

  Targeted therapy interventions using tyrosine kinase inhibitors (TKIs) provide encouraging treatment responses in ALK-rearranged lung adenocarcinomas, yet resistances occur almost inevitably. Apart from tumor cell-intrinsic resistance mechanisms, accumulating evidence supports a role of cancer-associated fibroblasts (CAFs) in affecting the therapeutic vulnerability of lung cancer cells. Here, we aimed to investigate underlying molecular networks shaping the therapeutic susceptibility of ALK-driven lung adenocarcinoma cells via tumor microenvironmental cues using three-dimensional (3D) spheroid co-culture settings. We show that CAFs promote therapy resistance of lung tumor cells against ALK inhibition by reducing apoptotic cell death and increasing cell proliferation. Using single-cell RNA-sequencing analysis, we show that genes involved in lipogenesis constitute the major transcriptional difference between TKI-treated homo- and heterotypic lung tumor spheroids. CAF-conditioned medium and CAF-secreted factors HGF and NRG1 were both able to promote resistance of 3D-cultured ALK-rearranged lung tumor cells via AKT signaling, which was accompanied by enhanced de novo lipogenesis and supression of lipid peroxidation. Notably, simultaneous targeting of ALK and SREBP-1 was able to overcome the established CAF-driven lipid metabolic-supportive niche of TKI-resistant lung tumor spheroids. Our findings highlight a crucial role of CAFs in mediating ALK-TKI resistance via lipid metabolic reprogramming and suggest new ways to overcome resistance towards molecular directed drugs by targeting vulnerabilities downstream of oncogenic signaling.

  Study EGAS50000000135

• Diagnostic Genomic Analysis is Prognostic in AYA ALL Patients Treated on a MRD-Stratified Paediatric Protocol

  Background: Over the last decade, next generation sequencing has enabled classification of multiple new recurrent genomic drivers of acute lymphoblastic leukemia (ALL). We provide a comprehensive genomic survey of adolescent and young adult (AYA) patients treated on the Australasian Leukaemia and Lymphoma Group (ALLG) ALL06 study. ALL06 clinical trial described the efficacy and tolerability of AYA ALL patients aged 15-39 years treated on a single minimal residual disease (MRD)-stratified platform. In this study we have taken the ALL06 cohort and retrospectively performed comprehensive genomic analyses assessing response of AYA patients with various genomic subtypes treated with a paediatric approach. Methods: Genomic risk classification of B- and T-ALL patients (n=57) enrolled to the study was based on mRNA-Sequencing, MLPA, immunophenotyping and cytogenetics. Results: Using this approach, 36/40 (90%) B- and 13/17 (76.5%) T-ALL patients were classified according to genomic risk. We have linked outcome data with each patient’s genomic risk and found a strong correlation between adverse genomic risk and MRD at the end of consolidation, translating to inferior overall and relapse free survival. Patients with adverse risk genomics who achieved negative MRD status had improved responses compared to those with persistent MRD. Patients with standard risk genomics had excellent responses regardless of their MRD status. Conclusion: Our data argue strongly for incorporation of a genomic risk classification into future ALL treatment paradigms at the time of diagnosis.

  Study EGAS50000000752

• Ultra-long sequencing for contiguous haplotype resolution of the human immunoglobulin heavy chain locus

  Genetic diversity within the human immunoglobulin heavy chain (IGH) locus influences the expressed antibody repertoire and susceptibility to infectious and autoimmune diseases. However, repetitive sequences and complex structural variation pose significant challenges for large-scale characterization. Here, we introduce a method that combines Oxford Nanopore ultra-long sequencing and adaptive sampling with a bioinformatic pipeline to produce haplotype-resolved, annotated IGH assemblies. Notably, our strategy overcomes prior limitations in phasing resolution, enabling single-contig haplotype assemblies that span the entire IGH locus. We applied this method to four individuals and validated the accuracy of the IGH assemblies using Pacific Biosciences HiFi reads, demonstrating near-complete sequence congruence, detecting only indels with a high degree of confidence. Moreover, when applied to the reference material HG002, our pipeline revealed no base differences and a limited number of indels compared with the Telomere-to-Telomere genome benchmark across the IGH region. Importantly, in the four individuals, our approach uncovered 30 novel alleles and previously uncharacterized large structural variants, including a 120 kb duplication spanning IGHE to IGHA1 and an expanded seven-copy IGHV3-23 gene haplotype. These findings underscore the power of our method to resolve the full complexity of the IGH locus and uncover previously unrecognized variants that may affect immune function and disease susceptibility. Thus, our method provides a strong basis for future immunological research and translational applications.

  Study EGAS50000001042

• Ultra-sensitive tumor-informed ctDNA monitoring of treatment response in advanced esophagogastric cancer patients

  To explore whether ultra-sensitive circulating tumor DNA (ctDNA) profiling enables earlier prediction of treatment response and detection of disease progression, we applied NeXT Personal, an ultra-sensitive bespoke tumor-informed liquid biopsy platform, to profile tumor samples from the KeyLargo study, a phase II trial in which metastatic esophagogastric cancer (mEGC) patients received capecitabine, oxaliplatin, and pembrolizumab. A total of 24 patients were evaluated, and all were ctDNA-positive at baseline. ctDNA levels varied from 406,067 down to 1.5 parts per million (PPM) with a median limit of detection of 2.03 PPM. ctDNA dynamics were highly correlated with changes in tumor size (ρ = 0.59, p = 7.3×10-9). Lack of early molecular response (50% or greater decrease in ctDNA levels at first available time point after 30 days, C2D1 or C4D1) was associated with worse overall survival (OS) (HR 6.6, 95% CI 1.8-24.1, p = 0.005) and progression-free survival (PFS) (HR 15.4, 95% CI 2.7-87.0, p = 0.002). Lack of molecular clearance of ctDNA was associated with worse OS (HR 6.9, 95% CI 1.5-30.8, p = 0.012) and PFS (HR 19.2, 95% CI 2.4-152.8, p = 0.005). Molecular progression (ctDNA increase) preceded imaging-derived progression by a median lead time of 65 days. These results suggest that ultra-sensitive liquid biopsy approaches could improve treatment decision-making for mEGC patients receiving chemotherapy and immunotherapy.

  Study EGAS50000001224

• Delineating pediatric brain tumor progression using single-nuclei sequencing

  Pediatric brain cancers are the most common solid tumors in children with varying survival rates depending on molecular features of cancer cells and associated microenvironment. Several studies have used single-cell RNA sequencing to define transcriptional programs active in cancer cells and highlighted supportive roles of non-cancerous cells embedded in these tumors. Although application of single-cell RNA sequencing technology has expanded our understanding of intratumoral heterogeneity and underlying transcriptional mechanisms in brain tumors at diagnosis, mechanisms of disease remains elusive. Therefore, the purpose of this proposal is to characterize cellular and transcriptional states of malignant and non-malignant compartments at primary diagnosis and at relapse. In this study, we performed single-nuclei sequencing of 42 patient samples taken at various disease stages (initial diagnosis, progressive disease, and relapsed disease) from 20 patients with low-grade glioma or Atypical Teratoid Rhabdoid Tumor (ATRT). Bulk whole genome and RNA-seq data for these samples are available from the Children’s Brain Tumor Network and the Gabriella Miller Kids First Data Resource Center. Findings from our study will help define the mechanistic shifts of tumor and non-cancerous cells over course of disease progression. Collectively with the single-cell Pediatric Cancer Atlas funded, this project lays the foundation of an important data set focused on lethal relapsed cancers that can continue to be expanded to identify cancer-intrinsic and extrinsic microenvironmental factors supporting recurrence of pediatric brain tumors.

  Study EGAS50000001288

• Gut metagenome associations with extensive digital health data in a volunteer-based EstMB cohort second timepoint data

  As part of the Estonian Biobank, we established the Estonian Microbiome Cohort which includes stool, oral and plasma samples from 2,509 participants and is supplemented with multi-omic measurements, questionnaires, and regular linkages to national electronic health records. Out of the 2,509 Estonian Microbiome Project participants, a sub-cohort of over 300 individuals provided an additional stool sample after a median follow-up period of 4.4 years. All participants of the EstMB cohort gave informed consent for the data and samples to be used for scientific purposes. In the second time point, the participants were instructed to send the samples by post immediately after taking the sample and time their sample collection to avoid shipping on the weekends. In the study, by utilising retrospective medication usage data from the electronic health records and the Estonian microbiome cohort shotgun metagenomics data set (n = 2,509), we systematically evaluate the long-term effects of antibiotics and human-targeted medications on the gut microbiome. We show that past usage of medications is associated with the gut microbiome. For example, the effects of antibiotics, psycholeptics, antidepressants, proton pump inhibitors, and beta-blockers are detectable several years after use. Furthermore, by analysing a subcohort (n >300 individuals) from the second time point, we show that similar changes in the gut microbiome occur after treatment initiation or discontinuation, possibly indicating causal effects.

  Study EGAS50000001181

• Genomic landscapes of endometrioid and mucinous ovarian cancers and morphologically similar tumor types

  While endometrioid and mucinous ovarian cancers represent nearly a quarter of ovarian cancers, their molecular characteristics and pathologic origins are poorly understood. To characterize the genomic and epigenomic alterations of these ovarian cancer subtypes and evaluate links to morphologically similar tumors from other sites, we performed sequence, copy number, mutation signature, rearrangement, and methylation analyses from tumor samples and matched normal tissues of 133 patients. These analyses included patients with ovarian endometrioid (n=44), ovarian mucinous (n=43), uterine endometrial (n=15), and gastrointestinal mucinous cancers (n=31). In addition to identifying genes previously known to be involved in these tumors, we identified alterations in RAD51C, NOTCH4, SMARCA1/4 and JAK1 in ovarian endometrioid, ESR1 in uterine endometrioid, and SMARCA4 in ovarian mucinous cancers. Whole genome sequencing revealed rearrangements involving PTEN, NF1, and NF2 in ovarian endometrioid tumors, and NF1 and MED1 in ovarian mucinous cancers. The number of alterations and affected genes in ovarian and uterine endometrioid cancers were similar as were genome-wide methylation profiles, suggesting that these tumors may have a similar tissue of origin. By contrast, methylation patterns in ovarian mucinous cancers were more similar to other gastrointestinal mucinous cancers, including mucinous cancers of the stomach, colon, and pancreas. These analyses provide insights into the genomic landscapes and origins of mucinous and endometrioid ovarian carcinomas, providing new avenues for early clinical intervention and management of patients with these cancers.

  Study EGAS50000000523

• Systematic Growth Factor Profiling Platform for 3D Tumor Models Reveals EstradiolResponsive Cellular Mechanisms of Immunotherapy Resistance

  Epithelial ovarian cancer (EOC) exhibits significant molecular heterogeneity, which is thought to be influenced by the tumor microenvironment (TME). However, the specific factors driving this heterogeneity remain poorly understood. In this study, we integrated growth factoromics and genomic data to identify TME factors contributing to EOC heterogeneity and immunosuppression. Our analysis revealed distinct subgroups of EOC samples with differential responsiveness to specific growth factors, particularly estradiol and R-spondin. Notably, samples with high responsiveness to estradiol and Wnt (C1 cluster) were enriched in the Mesenchymal (MES) molecular subtype, while samples with high responsiveness to R-spondin (C2 cluster) were associated with the Differentiated (DIF) and Immunoreactive (IMR) subtypes. We discovered that the G protein-coupled estrogen receptor (GPER), rather than the classical estrogen receptor alpha (ERα), was the key mediator of estrogen responsiveness and a determinant of molecular subtype classification. Single-cell sequencing analysis identified subsets of estrogen-responsive malignant cell (MAL.PDCD5) and cancer-associated fibroblast (FB.TNFSF10), which were associated with immunosuppressive patterns in the TCGA cohort. Importantly, the interaction between these cell populations significantly increased upon estrogen treatment, suggesting a potential mechanism by which estrogen signaling may contribute to the development of an immunosuppressive TME. Our findings highlight the critical role of estrogen signaling in shaping the molecular heterogeneity and immunosuppressive TME of EOC and provide novel insights into potential therapeutic targets for overcoming immunosuppression in this malignancy.

  Study EGAS50000000422

• Tumor-reactive heterotypic CD8 T cell clusters from clinical samples

  Emerging evidence suggests a correlation between CD8+ T cell-tumor cell proximity and anti-tumor immune response. However, it remains unknown whether these cells exist as functional clusters that can be isolated from clinical samples. By conventional and imaging flow cytometry, we show here that from 21/21 human melanoma metastases, we could isolate heterotypic clusters, comprising CD8+ T cells interacting with one or more tumor cells and/or antigen-presenting cells (APCs). Single-cell RNA-sequencing revealed that T cells from clusters were enriched for gene signatures associated with tumor reactivity and exhaustion. Clustered T cells exhibited increased TCR clonality indicative of expansion, while TCR-matched T cells showed more exhaustion and co-modulation when conjugated to APCs than to tumor cells. T cells that were expanded from clusters ex vivo exerted on average 9-fold increased killing activity towards autologous melanomas, which was accompanied by enhanced cytokine production. Upon adoptive cell transfer into mice, T cells from clusters showed improved patient-derived melanoma control, which was associated with increased T cell infiltration and activation. Together, these results demonstrate that tumor-reactive CD8+ T cells are enriched in functional clusters with tumor cells and/or APCs and that they can be isolated and expanded from clinical samples. Typically excluded by single cell gating in flow cytometry, these distinct heterotypic T cell clusters serve as a valuable source to decipher functional tumor-immune cell interactions, while they may also be therapeutically explored.

  Study EGAS50000000785

• The circulating cell-free DNA landscape in sepsis is dominated by impaired liver clearance

  Circulating cell-free DNA (cfDNA) is a promising molecular biomarker. However, its utility in severe infection remains poorly understood. Here, we isolated cfDNA from sepsis patients and controls, demonstrating a 41-fold increase in the amount of cfDNA in circulation during disease. We used sequencing to reconstruct cfDNA methylomes, fragmentation profiles, and nucleosome footprints across 58 samples. We observed no difference in cfDNA composition between patients and controls, challenging the idea that cfDNA increases due to higher immune cell death during sepsis. Instead, we suggest that liver dysfunction prevents efficient clearance of cfDNA during disease. This was supported by fragmentation and end-motif patterns, both of which showed evidence of cfDNA being exposed to circulating nucleases for a prolonged period, proportionally to the extent of liver dysfunction. Variation in cfDNA of megakaryocyte-erythroid progenitor origin was also a significant contributor in sepsis, increasing over time. Moreover, we showed that cfDNA retains nucleosome footprints with cell type-specific gene activity information. We developed a novel approach to study nucleosome phasing that successfully recovers tissue-specific signatures. By combining this with single-cell data, we demonstrated that sepsis patients with liver dysfunction have higher amounts of cfDNA derived from Kupffer cells and the liver parenchyma. In conclusion, we present the first high-throughput multi-modal study of cfDNA during sepsis, which will serve as a reference point for future studies on the role of this biomarker in critical illness.

  Study EGAS50000001033

• Relapse CHL study

  Classical Hodgkin Lymphoma (CHL) is characterized by a complex tumor microenvironment (TME) that supports disease progression. While immune cell recruitment by Hodgkin and Reed–Sternberg (HRS) cells is well-documented, the role of non-malignant B cells in relapse remains unclear. Using single-cell RNA sequencing (scRNA-seq) on paired diagnostic and relapsed CHL samples, we identified distinct shifts in B cell populations, particularly an enrichment of naïve B cells and a reduction of memory B cells in early-relapse compared to late-relapse and newly diagnosed CHL. The enrichment of naïve B cells in early relapse biopsies was confirmed in independent validation cohorts using scRNA-seq and immunohistochemistry. Notably, naïve B cells in early-relapse samples exhibited high expression of LGALS9, an immunosuppressive molecule known to interact with T cell immunoglobulin and mucin domain 3 (TIM-3) on regulatory T cells (Tregs). Cell-cell interaction analysis revealed the importance of interactions between LGALS9+ naïve B cells and TIM-3+ Tregs in the early-relapse setting. Spatial analysis by imaging mass cytometry confirmed close proximity of LGALS9+ naïve B cells with TIM-3⁺ CD4⁺ T cells and HRS cells, pointing to their role in shaping an immunosuppressive niche. Our findings highlight a previously unrecognized population of LGALS9⁺ naïve B cells with immunoregulatory potential in early-relapse CHL and provide new insights into the spatial and transcriptional architecture of the relapsed TME in CHL.

  Study EGAS00001008222

• Evolution of the African pygmy phenotype

  The evolutionary history of the human ‘pygmy’ phenotype (small body size), a characteristic of African and Southeast Asian rainforest hunter-gatherers, is largely unknown. Here we use a genome-wide admixture mapping analysis to identify 16 genomic regions that are significantly associated with the pygmy phenotype in the Batwa, a rainforest hunter-gatherer population from Uganda (East Central Africa). The identified genomic regions have multiple attributes that provide supporting evidence of genuine association with the pygmy phenotype, including enrichments for SNPs previously associated with stature variation in Europeans and for genes with growth hormone receptor and regulation functions. To test adaptive evolutionary hypotheses we computed the haplotype-based iHS statistic and the level of population differentiation (FST) between the Batwa and their agricultural neighbors, the Bakiga, for each SNP. Both |iHS| and FST values were significantly higher for SNPs within the Batwa pygmy phenotype-associated regions than the remainder of the genome, a signature of polygenic adaptation. In contrast, when we expanded our analysis to include Baka rainforest hunter-gatherers from Cameroon and Gabon (West Central Africa) and Nzebi and Nzime neighboring agriculturalists, we did not observe elevated |iHS| or FST values in these genomic regions. Together, these results suggest adaptive and at least partially convergent origins of the pygmy phenotype even within Africa, supporting the hypothesis that small body size confers a selective advantage for tropical rainforest hunter-gatherers but raising questions about the antiquity of this behavior.

  Study EGAS00001000908

• Evolutionary origin and methylation status of human intronic CpG islands that are not present in mouse

  Imprinting of the human RB1 gene is due to the presence of a differentially methylated CpG island (CGI) in intron 2, which is part of a retrogene derived from the PPP1R26 gene on chromosome 9. The murine Rb1 gene does not have this retrogene and is not imprinted. We have investigated whether the RB1/Rb1 locus is unique with respect to these differences. For this we have compared the CGIs from human and mouse by in silico analyses and found out that there are at least 2033 human intronic CGIs that are not present in the mouse. Among these CGIs, 104 show sequences similarities elsewhere in the human genome, which suggests that they arose by retrotransposition. We could narrow down the time points when most of these CGIs appeared during evolution. Their methylation status was analysed in two monocyte methylome datasets from whole genome bisulfite sequencing and in 18 published methylomes. Five CGIs appear to be differentially methylated, which suggests that they are subject of genomic imprinting. They are located in introns of the following genes: RB1, ASRGL1, PARP11, PDXDC1 and MYO1D. Four of these CGIs are parts of retrogenes. Our study supports the notion that imprinting builds on host defence mechanisms by which the genome protects itself against retrotransposed and foreign DNA and that the epigenetic fate of the retrotransposed DNA depends on selective forces at the integration site.

  Study EGAS00001000719

• UK10K_RARE_HYPERCHOL

  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.Familial Hypercholesterolemia is a condition where the affected person has a consistently high level of LDL which can lead to early clogging of the coronary arteries. All patients selected for this study will have been found not to carry the common APOB and PCSK9 mutations, and to have no detectable LDLR mutations by testing for 18 common mutations. For further information with regard to this cohort please contact Steve Humphries(steve.humphries@ucl.ac.uk).

  Study EGAS00001000129

• Integration of genomics and histology reveals diagnosis and effective therapy of refractory cancer of unknown primary with PDL1 amplification (H021)

  Identification of the tissue of origin in cancer of unknown primary (CUP) poses a diagnostic challenge and is critical for directing site-specific therapy. Currently, clinical decision making in patients with CUP primarily relies on histopathology and clinical features. Comprehensive molecular profiling has the potential to contribute to diagnostic categorization and, most importantly, guide CUP therapy through identification of actionable lesions. We here report the case of an advanced-stage malignancy initially mimicking poorly differentiated soft-tissue sarcoma that did not respond to multi-agent chemotherapy. Molecular profiling within a clinical whole-exome and transcriptome sequencing program revealed a heterozygous, highly amplified KRAS G12S mutation, compound-heterozygous TP53 mutation/deletion, high mutational load, and focal amplification of chromosomes 9p (including PDL1 [CD274] and JAK2) and 10p (including GATA3). Integrated analysis of molecular data and conventional histopathology suggested a diagnosis of triple-negative breast cancer (TNBC) and provided a rationale for immune checkpoint inhibitor therapy with pembrolizumab, which resulted in rapid clinical improvement and a lasting partial remission. Analysis of 157 TNBC samples from The Cancer Genome Atlas revealed focal, high-level PDL1 amplification coinciding with excessive PDL1 mRNA expression in 24% of cases. Collectively, these results illustrate the diagnostic utility of multidimensional tumor profiling in cases with non-descript histology and immune phenotype, demonstrate the predictive power of genomic PDL1 amplification for immune checkpoint inhibition, and suggest a targeted therapeutic strategy in chromosome 9p24.1/PDL1-amplified cancers.H021

  Study EGAS00001001846

• Negligible impact on missing heritability of autoimmune-locus rare coding-region variants

  Genome wide association studies (GWAS) have identified common variants of modest effect size at hundreds of loci for common autoimmune diseases - however a substantial fraction of heritability remains unexplained, and to which rare variants may contribute. To discover rare variants and test them for association with a phenotype, the majority of studies re-sequence a small initial sample size and then genotype the discovered variants in a larger sample set. This approach will fail to analyse a large fraction of the rare variants present in the entire sample set. Here we perform simultaneous amplicon sequencing-based variant discovery and genotyping for coding exons of 25 GWAS risk genes in 41,911 white-European origin UK subjects comprising 24,892 subjects with six autoimmune disease phenotypes and 17,019 controls, and show that rare coding-region variants at known loci play a negligible role in common autoimmune disease susceptibility. These results do not support the rare variant - synthetic genome-wide association hypothesis. Many known autoimmune disease risk loci contain multiple independent common and low frequency variant association signals, and so genes in these loci are a priori stronger candidates to harbor rare coding-region variants than other genes. Our data suggest that the missing heritability for common autoimmune diseases may not lie in the rare coding-region variant portion of the allelic spectrum, but perhaps as others have proposed in very many common variant signals of weak effect.

  Study EGAS00001000476

• Genome of the Netherlands

  The Genome of the Netherlands (GoNL) Project characterizes DNA sequence variation, common and rare, for SNVs, short insertions and deletions (indels) and larger deletions in 769 individuals of Dutch ancestry selected from five biobanks under the auspices of the Dutch hub of the Biobanking and Biomolecular Research Infrastructure (BBMRI-NL). The samples come from a representative sample of 250 trio-families from all provinces in the Netherlands. The parent-offspring trios include adult individuals ranging in age from 19 to 87 years (mean=53 years; SD=16 years) from birth cohorts 1910-1994. Sequencing was done on blood-derived DNA from uncultured cells and accomplished coverage was 14-15x. Samples where contributed by LifeLines (http://lifelines.nl/lifelines-research/general), The Leiden Longevity Study (http://www.healthy-ageing.nl; http://www.langleven.net), The Netherlands Twin Registry (NTR: http://www.tweelingenregister. org), The Rotterdam studies, (http://www.erasmus-epidemiology.nl/rotterdamstudy) and the Genetic Research in Isolated Populations program (http://www.epib.nl/research/geneticepi/research.html#gip). The sequencing was carried out in collaboration with the Beijing Institute for Genomics (BGI). The analysis was done by a consortium lead by UMCG, LUMC, Erasmus MC, VU university and UMCU, see http://www.nlgenome.nl. Funding for the project was provided by the Netherlands Organization for Scientific Research under award number 184021007, dated July 9, 2009 and made available as a Rainbow Project of the Biobanking and Biomolecular Research Infrastructure Netherlands (BBMRI-NL).

  Study EGAS00001000644

• Identification_of_synthetic_lethal_genes_by_CRISPR_Cas9_library

  Sequencing is performed on PCR products from the amplification of gRNA segments introduced into cell lines (some are commercial, some are derived by collaborators, as described in HMDMC 13/076). Cellular DNA or RNA are not sequenced. This data is defined as "Other Biological/Biochemical Assay Data" under the WTSI datasharing guidelines and consequently will be released at the time of publication and/or through faculty websites. SEQUENCING DETAILS: in October 2016 3 DNA samples will be submitted for PCR amplification and purification through CGP pipeline using Thompson's custom primer set and PCR settings and sequencing with Hiseq using a custom protocol (with 3 forward reads to read sgRNA1, sgRNA2 and the index, optimized with Illumina Bespoke Sanger team); read length 50bp. The following set of samples will be submitted in later 2016 and ealry 2017 as cell pellet to be processed through CGP Cas9 pipeline (Thompson's oligos and protocols) and sequences on Hiseq (same custom sequencing protocol described above) We are utilizing a custom CRISPR/Cas9 library to identify synthetic lethal gene pairs. Cells will be transduced with LV expressing Cas9, then infected with LV expressign the library of 42'000 sgRNA pairs. Cell pellet will be collected from the cell population at different timepoints. The samples will be processed by CGP pipeline customized for the specific protocol and then sequenced to identify the 19bp of the sgRNA and quantify their representation.

  Study EGAS00001002117

• Genetic and epigenetic characterization of adenoid cystic carcinoma

  Translocation events are frequent in cancer and may create chimeric fusions or regulatory rearrangements that drive oncogene overexpression. Although regulatory rearrangements are increasingly recognized in hematopoietic and even solid tumors, the underlying mechanisms remain obscure. Here we identify super-enhancer translocations that drive overexpression of the oncogenic transcription factor MYB as a unifying theme in adenoid cystic carcinoma (ACC). Whole genome sequencing data and chromatin state maps for 13 primary tumors and xenografts reveal distinct chromosomal rearrangements that juxtapose super-enhancer clusters to the MYB locus. Chromosome conformation capture confirms that the translocated enhancers physically interact with the MYB promoter. Remarkably, the MYB product itself binds to the translocated enhancers, thus creating a positive feedback loop that sustains its own expression. MYB also binds a large number of active enhancers that drive opposing regulatory programs in alternate cell lineages in ACC. MYB cooperates with the transcription factor TP63 in the myoepithelial component of the tumors, but promotes a Notch program in the luminal epithelial component. Bromodomain inhibitors, which disrupt enhancer function, slow tumor growth in ACC primagraft models in vivo, but are ineffective against high grade tumors that harbor activating mutations in the Notch pathway. Thus, our study identifies super-enhancer translocations as a unifying feature of ACC, and provides insight into the mechanism by which sustained MYB overexpression drives alternate cell fates in this disease.Processed ChIP-seq data is available on GEO under accession number GSE76465.

  Study EGAS00001001457

• Tissue-specific mutation accumulation in human adult stem cells during life

  Gradual accumulation of mutations in human adult stem cells during life is associated with various age-related diseases, including cancer. The number of stem cell divisions throughout life is believed to be a major determinant for mutation accumulation and could explain the extreme variation of cancer incidence across different organs. Yet, mutation patterns and rates of healthy adult stem cells remain unknown. Here, we determined genome-wide mutation patterns in primary adult stem cells of the small intestine, colon and liver of human donors with ages ranging from 3 to 87 years. We find that the number of mutations increases linearly with age up to several thousand mutations per cell at 87 years of age, while mutation spectra remain constant throughout life. Small intestine and colon stem cells have a 2-fold higher mutation rate per year compared with liver stem cells. These differences could be exclusively attributed to the mutagenic action of spontaneous deamination of cytosine residues and may reflect the high stem cell division rate in these tissues. The genomic distribution of somatic mutations is non-random and predominantly associated with DNA replication dynamics in the small intestine and colon, and with transcription in the liver. These results indicate that a stable balance between various mutagenic and DNA repair processes is maintained throughout life and that the activity of these processes in adult stem cells varies between tissues.

  Study EGAS00001001682

• Plasma-Seq of patients with metastatic prostate cancer

  Study 1 In this study we analysed patients with metastatic prostate cancer to scan their tumor genomes noninvasively in plasma DNA. We wanted to make whole-genome analysis from plasma DNA amenable to clinical routine applications and developed an approach based on a benchtop high-throughput platform, i.e. Illuminas MiSeq instrument. We performed whole-genome sequencing from plasma at a shallow sequencing depth to establish a genome-wide copy number profile of the tumor at low costs within 2 days. The genome-wide profiling in the plasma of our patients revealed multiple copy number aberrations including those previously reported in prostate tumors, such as losses in 8p and gains in 8q. High-level copy number gains in the AR locus were observed in patients with CRPC but not with CSPC disease. We identified the TMPRSS2-ERG rearrangement associated 3-Mbp deletion on chromosome 21 and found corresponding fusion plasma fragments in these cases. In an index case multiregional sequencing of the primary tumor identified different copy number changes in each sector, suggesting multifocal disease. Our plasma analyses of this index case, performed 13 years after resection of the primary tumor, revealed novel chromosomal rearrangements, which were stable in serial plasma analyses over a 9 months period, which is consistent with the presence of one metastatic clone. Our approach provides specific genomic signatures within 2 days which may therefore serve as “liquid biopsy”.

  Study EGAS00001000451

• RB1 Gene Inactivation by Chromothripsis in Human Retinoblastoma

  Retinoblastoma is a rare childhood cancer of the developing retina. Most retinoblastomas initiate with biallelic inactivation of the RB1 gene through diverse mechanisms including point mutations, nucleotide insertions, deletions, loss of heterozygosity and promoter hypermethylation. Recently, a novel mechanism of retinoblastoma initiation was proposed. Gallie and colleagues discovered that a small proportion of retinoblastomas lack RB1 mutations and had MYCN amplification [1]. In this study, we identified recurrent chromosomal, regional and focal genomic lesions in 94 primary retinoblastomas with their matched normal DNA using SNP 6.0 chips. We also analyzed the RB1 gene mutations and compared the mechanism of RB1 inactivation to the recurrent copy number variations in the retinoblastoma genome. In addition to the previously described focal amplification of MYCN and deletions in RB1 and BCOR, we also identified recurrent focal amplification of OTX2, a transcription factor required for retinal photoreceptor development. We identified 10 retinoblastomas in our cohort that lacked RB1 point mutations or indels. We performed whole genome sequencing on those 10 tumors and their corresponding germline DNA. In one of the tumors, the RB1 gene was unaltered, the MYCN gene was amplified and RB1 protein was expressed in the nuclei of the tumor cells. In addition, several tumors had complex patterns of structural variations and we identified 3 tumors with chromothripsis at the RB1 locus. This is the first report of chromothripsis as a mechanism for RB1 gene inactivation in cancer.

  Study EGAS00001000598

• UK10K COHORT TWINSUK

  The UK10K project proposes 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 directly analyse 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 TwinsUK samples will be part of the cohort study and will undergo whole genome sequencing. For further information with regard to this cohort please contact Brent Richards (brent.richards@kcl.ac.uk) or Nicole Soranzo (ns6@sanger.ac.uk).

  Study EGAS00001000108

• UK10K NEURO ASD GALLAGHER

  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 is an Irish sample set of individuals with ASD (approximately 50% with comorbid intellectual disability). Individuals have been diagnosed with ADI/ ADOS, measures of cognition/ adaptive function. They represent a more severe, narrowly defined cohort of ASD subjects. Family histories are available for some with measures of broader phenotype. For further information on this cohort please contact Nadia Bolshakova (bolshakn@tcd.ie).

  Study EGAS00001000112

• UK10K NEURO ASD BIONED

  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 BioNED (Biomarkers for Childhood onset neuropsychiatric disorders) study has been carrying out detailed phenotypic assessments evaluating children with an autism spectrum disorder. These assessments included ADI-R, ADOS, neuropsychology, EEG etc. There are 56 DNA samples from this study (25 extracted from blood). For further information with regard to this cohort please contact Patrick Bolton (patrick.bolton@kcl.ac.uk).

  Study EGAS00001000111

• Whole genome study of Hurthle cell thyroid carcinoma

  Oncocytic thyroid carcinoma, also known as Hürthle cell thyroid carcinoma, accounts for only a small percentage of all thyroid cancers. However, this malignancy often presents at an advanced stage and poses unique challenges to patients and clinicians. Surgical resection of the tumor accompanied in some cases by radioactive iodine treatment, radiation and chemotherapy are the established modes of therapy. Knowledge of the perturbed oncogenic pathways can provide better understanding of the mechanism of disease and thus opportunities for more effective clinical management. Initially, two oncocytic thyroid carcinomas and their matched normal tissues were profiled using whole genome sequencing. Subsequently, 72 oncocytic thyroid carcinomas, one cell line and 5 Hürthle cell adenomas were examined by targeted sequencing for the presence of mutations in multiple endocrine neoplasia I (MEN1) gene. We report the identification of MEN1 loss-of-function mutations in approximately 7% of patients diagnosed with oncocytic thyroid carcinoma. Whole genome sequence data also revealed large regions of copy number variation encompassing nearly the entire genomes of these tumours. Menin, a ubiquitously expressed nuclear protein, is a well-characterized tumor suppressor whose loss is the cause of multiple endocrine neoplasia type 1 syndrome. Menin is involved in several major cellular pathways such as regulation of transcription, control of cell cycle, apoptosis and DNA damage repair pathways. Mutations of this gene in a subset of Hürthle cell tumors point to a potential role for this protein and its associated pathways in thyroid tumorigenesis.

  Study EGAS00001000940

• Genomic Alterations in Gingivo-buccal Cancer: ICGC-India Project_YR01

  As a part of the ICGC, India has undertaken genomic studies on gingivobuccal cancer of the oral cavity, which is the most prevalent form of cancer among men in India. There are various known environmental (life-style) correlates of this cancer, the most important of which are tobacco chewing and HPV infection. Paired DNA samples – isolated from the tumor tissue and from the blood of fifty patients – have been analyzed to catalog germline and somatic mutations. Association between each observed genomic alteration and exposure to environmental risk factors is being explored.Detailed clinical characterization of the patients, collection of data on demographic and environmental exposures, and isolation of DNA samples from blood and tumor tissues collected from each patient are being done at the Advanced Centre for Research, Treatment and Education on Cancer, Mumbai. Exome capture and deep resequencing are being performed at the National Institute of Biomedical Genomics, Kalyani. We have analyzed the exomes of paired blood and tumor DNA samples using two orthogonal platforms – exome capture by Roche-Nimblegen SeqCap Ez kit followed by sequencing on Roche-454 GS FLX, and exome capture by Illumina TrueSeq Exome Enrichment kit followed by sequencing on Illumina HiSeq-2000. Each exome is being sequenced at a mean depth of 30x. In addition, we have used Illumina Omni Quad SNP-chips to generate genotype data at 1.14 million genomic markers. The average concordance between SNP-chip and sequence data is 99.8%.

  Study EGAS00001000249

• Illumina HumanCoreExome genotyping data from the TEENAGE (TEENs of Attica: Genes and Environment) study.

  The TEENAGE study is a cross-sectional study comprising healthy Greek adolescents. The study was conducted according to the guidelines specified in the Declaration of Helsinki and all procedures involving human subjects were approved by the Institutional Review Board of Harokopio University and the Greek Ministry of Education, Lifelong Learning and Religious Affairs. The TEENAGE study target population comprised adolescent students aged 13 to 15 years attending the first three classes of public secondary schools located in the wider Athens area of Attica. The regions across Attica from which students were recruited were physically similar and all urban in nature. Sampling of regions and schools selected to participate was random. Students were enrolled from the wider urban region to ensure that they were coming from families with varying educational background and socio-economic status. Moreover, in order to exclude any vast economic disparities, all participating schools were exclusively public schools, which are free of charge. A total of twenty-four schools participated in the study out of 285 junior secondary public schools in Attica. Anthropometric, medical and dietary information was obtained. DNA samples were genotyped using Illumina HumanCoreExome BeadChips (Illumina, San Diego, CA, USA) at the Wellcome Trust Sanger Institute, Hinxton, UK. Genotypes were called using GenCall (Illumina Genome Studio) followed by zCall and quality control (QC) was performed in two stages (pre- and post-zCall).

  Study EGAS00001001001

• Efficacy of JAK/STAT pathway inhibition in murine xenograft models of early T-cell precursor (ETP) acute lymphoblastic leukemia

  Early T-cell precursor (ETP) acute lymphoblastic leukemia (ALL) is a recently described subtype of T-ALL characterized by a unique immunophenotype and genomic profile as well as a high rate of induction failure. Frequent mutations in cytokine receptor and JAK/STAT signaling pathways led us to hypothesize that ETP-ALL is dependent on JAK/STAT signaling. Here we demonstrate aberrant activation of the JAK/STAT pathway in ETP-ALL blasts relative to non-ETP T-ALL. Moreover, ETP-ALL showed hyperactivation of STAT5 in response to IL7, an effect that was abrogated by the JAK1/2 inhibitor ruxolitinib. In vivo, ruxolitinib displayed activity in 6/6 patient-derived murine xenograft models of ETP-ALL, with profound single-agent efficacy in 5 models. Ruxolitinib treatment decreased peripheral blast counts relative to pre-treatment levels and compared to control (P<0.01) in 5/6 ETP-ALL xenografts, with marked reduction in mean splenic blast counts (P<0.01) in 6/6 samples. Surprisingly, both JAK/STAT pathway activation and ruxolitinib efficacy were independent of the presence of JAK/STAT pathway mutations, raising the possibility that the therapeutic potential of ruxolitinib in ETP-ALL extends beyond those cases with JAK mutations. These findings establish the preclinical in vivo efficacy of ruxolitinib in ETP-ALL, a biologically distinct subtype for which novel therapies are needed.

  Study EGAS00001001146

• The clonal and mutational evolution spectrum of primary triple negative breast cancers

  Primary triple negative breast cancers (TNBC) represent approximately 16% of all breast cancers and are a tumour type defined by exclusion,for which comprehensive landscapes of somatic mutation have not been determined. Here we show in 104 early TNBC cases, that at the time of diagnosis these cancers exhibit a wide and continuous spectrum of genomic evolution, with some exhibiting only a handful of somatic aberrations in a few pathways, whereas others contain hundreds of somatic events and multiple pathways implicated. Integration with matched whole transcriptome sequence data revealed that only 36% of mutations are expressed. By examining single nucleotide variant (SNV) allelic abundance derived from deep re-sequencing (median > 20,000 fold) measurements in 2414 somatic mutations, we determine for the first time in an epithelial tumour, the relative abundance of clonal genotypes among cases in the population. We show that TNBC vary widely and continuously in their clonal frequencies at the time of diagnosis, with basal subtype TNBC exhibiting more variation than non-basal TNBC. Although p53 and somatic mutations appear clonally dominant compared with other pathways, in some tumours their clonal frequencies are incompatible with founder status. Mutations in cytoskeletal and cell shape/motility proteins occurred at lower clonal frequencies, suggesting they occurred later during tumour progression. Taken together our results show that future attempts to dissect the biology and therapeutic responses of TNBC will require the determination of individual tumour clonal genotypes.

  Study EGAS00001000132

• 33 patients with Monoclonal Gammopathy of Undetermined Significance (MGUS) had whole exome sequencing performed. Cells were sorted by FACSAria using CD138, CD19 and CD56 to obtain a pure abnormal plasma cell population. Generated somatic variants were compared to a previous study of 463 multiple myeloma patients.

  Monoclonal gammopathy of undetermined significance (MGUS) is a premalignant precursor of multiple myeloma (MM) with a 1% risk of progression per year. Although targeted analyses have shown the presence of specific genetic abnormalities such as IGH translocations, RB1 deletion, 1q gain, hyperdiploidy or RAS genes mutations, little is known about molecular mechanism of malignant transformation. We have performed whole exome sequencing together with SNP array analysis in 33 flow-cytometry separated abnormal PC samples of MGUS patients to describe somatic gene mutations and chromosome changes at the genome-wide level. Non-synonymous mutations (NS-SNVs) and copy number alterations (CNAs) were present in 97.0% and in 63.6% of cases, respectively. Importantly, the number of somatic mutations was significantly lower in MGUS compared to MM (p<10-4) and we have identified 6 myeloma significantly mutated genes which are KRAS, NRAS, DIS3, HIST1H1E, EGR1 and LTB in the MGUS dataset. We also found a positive correlation with increasing chromosome changes and somatic mutations. IGH translocations were present in 27.3% of cases comprising t(4;14), t(11;14), t(14;16) or t(14;20) and were in a similar frequency to MM, which corresponded with primary lesion hypothesis. Data from this study showed MGUS is a genetically comprehensive disease, however overall genetic instability is significantly lower compared to MM.

  Study EGAS00001001658

• Combining transcription factor binding affinities with open chromatin data for accurate gene expression prediction

  The binding and contribution of transcription factors (TF) to cell specific gene expression is often deduced from open-chromatin measurements to avoid cost and labour intensive TF ChIP-seq assays.It is important to develop reliable and fast computational methods for accurate TF binding prediction in open-chromatin regions (OCRs). Here, we report a novel segmentation-based method, TEPIC, to predict TF binding by combining sets of OCRs with position weight matrices.TEPIC can be applied to various open-chromatin data, e.g. DNaseI-seq and NOMe-seq, using either peaks or footprints as input.In addition to open-chromatin data, also Histone-Marks (HMs) can be used in TEPIC to identify candidate TF binding sites.TEPIC computes TF affinities and uses open-chromatin/HM signal intensity as quantitative measures of TF binding strength.Using machine learning techniques, we show that incorporating low affinity binding sites improves our ability to explain gene expression variability compared to the standard presence/absence classification of binding sites.Further, we show that both footprints and peaks capture essential TF binding events and lead to a good prediction performance.In our application, gene-based scores computed by TEPIC with one open-chromatin assay nearly reach the quality of several TF ChIP-seq datasets.Finally, we show that these scores correctly predict known transcriptional regulators as illustrated by the application to novel DNaseI-seq and NOMe-seq data for primary human hepatocytes and CD4+ T-cells, respectively.

  Study EGAS00001002073

• Transcriptomic profiles of chronic lymphocytic leukemia before and after frontline therapy: 5-year results from the randomized CLL14 study

  CLL14 (NCT02242942) is a multinational, open-label Phase III study designed to compare venetoclax-obinutuzumab (Ven-Obi) versus chlorambucil-obinutuzumab (Clb-Obi) in previously untreated patients with CLL and co-existing conditions. Between Aug 7, 2015, and Aug 4, 2016, 432 patients were enrolled and randomly assigned to receive either venetoclax plus obinutuzumab (Ven-Obi, n=216) or chlorambucil plus obinutuzumab (Clb-Obi, n=216). In both study arms, patients received fixed-duration treatment of overall 12 cycles, each cycle 28 days long, and entered subsequent post-treatment surveillance. Peripheral blood samples were collected before the start of therapy, as well as every 3-6 months during and after therapy, including at the time of disease progression. The primary endpoint of the study was progression-free survival, secondary endpoints included overall response rate, minimal residual disease response, and overall survival. The primary read-out was conducted after a median follow-up of 28.1 months, when the primary endpoint was met with 30 primary end-point events (disease progression or death) in the Ven-Obi arm and 77 in the Clb-Obi arm (hazard ratio, 0.35; 95% confidence interval [CI], 0.23 to 0.53;P<0.001). Follow-up is ongoing until 9 years after the last patient enrollment. In an exploratory analysis, bulk RNAseq was conducted with CD19-enriched blood from patients before start of therapy (418 patients) and at first relapse (44 patients).

  Study EGAS00001006596

• T cell receptor repertoire sequencing reveals chemotherapy-driven clonal expansion in colorectal liver metastases

  Colorectal liver metastasis (CLM) is a leading cause of colorectal cancer mortality, and the response to immune checkpoint inhibition (ICI) in microsatellite stable CRC has been disappointing. Administration of cytotoxic chemotherapy may cause increased density of tumour infiltrating T cells, which has been associated with improved response to ICI. This study aimed to quantify and characterize T cell infiltration in CLM using T cell receptor (TCR) repertoire sequencing. Eighty-five resected CLM from patients included in the Oslo CoMet study were subjected to TCR repertoire sequencing. Thirty-five and 15 patients had received neoadjuvant chemotherapy (NACT) within a short or long interval, respectively, prior to resection, while 35 patients had not been exposed to NACT. T cell fractions were calculated, repertoire clonality was analysed based on Hill evenness curves, and TCR sequential convergence was assessed using network analysis. Increased T cell fractions (10.6% vs 6.3%) were detected in CLM exposed to NACT within a short interval prior to resection, while modestly increased clonality was observed in NACT exposed tumours independently of the timing of NACT administration and surgery. While private clones made up >90% of detected clones, network connectivity analysis revealed that public clones contributed the majority of TCR sequence convergence. TCR repertoire sequencing can be used to quantify T cell infiltration and clonality in clinical samples. This study provides evidence to support chemotherapy-driven T cell clonal expansion in CLM in a clinical context.

  Study EGAS00001007136

• Assessing mitochondrial bioenergetics in coronary artery disease: A translational left ventricular tissue study in humans (The AMBITION study).

  Severe or recurrent myocardial ischemia can lead to long-term left ventricular (LV) dysfunction and heart failure in patients with coronary artery disease (CAD). In this study, we assessed for chronic metabolomic and transcriptomic adaptations within LV myocardium of patients with evidence of preoperative inducible ischemia undergoing coronary artery bypass grafting (CABG). During surgery, paired transmural LV biopsies were acquired on the beating heart from a region of inducible ischemia (determined by stress perfusion cardiovascular magnetic resonance), and a remote LV segment. From 33 patients, 63 LV biopsies were acquired on the beating heart during CABG, comparing to analysis of LV samples from 11 donor hearts. The global myocardial ATP/ADP ratio was reduced in CAD patients as compared to donor LV tissue with increased expression of oxidative phosphorylation (OXPHOS) genes encoding the electron transport chain complexes across multiple cell types. Paired analyses of biopsies obtained during CABG from LV segments with or without inducible ischemia revealed no significant difference in the ATP/ADP ratio, broader metabolic profile or expression of ventricular cardiomyocyte genes implicated in OXPHOS. Finally, differential metabolite analysis suggested dysregulation of several intermediates in patients with impaired LVEF; this notably included the tricyclic acid cycle metabolite, succinate. Overall, our results suggest that viable human myocardium in patients with stable CAD has global alterations in bioenergetic and transcriptional profile without large regional differences between areas with or without inducible ischemia.

  Study EGAS00001007351

• Whole exome sequencing data for paired non-serous endometrial and ovarian carcinomas from 27 patients with concurrent tumours.

  Introduction: Concurrent non-serous endometrial and ovarian tumours are often managed clinically as two separate primary tumours, but almost all exhibit evidence of a genomic relationship. Methodology: This study investigates the extent of relatedness using whole exome sequencing, which was performed on paired non-serous endometrial and ovarian carcinomas from 27 patients with concurrent tumours in a cohort with detailed clinicopathological annotation. Four whole exome sequencing-derived parameters (mutation, mutational burden, mutational signatures and mutant allele tumour heterogeneity scores) were used to develop a novel unsupervised model for the assessment of genomic similarity to infer genomic relatedness of paired tumours. Results: This novel model demonstrated genomic relatedness across all four parameters in all tumour pairs. Mutations in PTEN, ARID1A, CTNNB1, KMT2D and PIK3CA occurred most frequently and 24 of 27 (89 %) tumour pairs shared identical mutations in at least one of these genes, with all pairs sharing mutations in a number of other genes. Ovarian endometriosis, CTNNB1 exon 3 mutation, and progression and death from disease were present across the similarity ranking. Mismatch repair deficiency was associated with less genomically similar pairs. Discussion: Although there was diversity across the cohort, the presence of genomic similarity in all paired tumours supports the hypothesis that concurrent non-serous endometrial and ovarian carcinomas are genomically related and may have arisen from a common origin. Keywords: Co-existing; Concurrent; Endometrial carcinoma; Genomics; Mutation; Ovarian carcinoma; Paired tumours; Tumour mutational burden.

  Study EGAS00001008259

• Molecular Attributes Underlying Central Nervous System and Systemic Relapse in Diffuse Large B-cell Lymphoma

  Central nervous system (CNS) relapse in diffuse large B-cell lymphoma remains a significant conundrum. While clinical risk scores and biomarkers allow the identification of patients with higher risk of CNS relapse, their predictive ability is low. Herein, we report findings from a comprehensive study examining gene expression changes, mutations and clonal evolution. Firstly, considering diagnostic samples from patients with CNS relapse, systemic relapse or no relapse, we found a striking enrichment of gene expression changes reflecting protein biosynthesis in patients who experienced CNS or systemic relapse. Furthermore, diagnostic specimens of relapsing cases were characterized by low immune infiltrates. The most discriminative gene feature in our analysis was SLAMF1, high expression of which was associated with favourable outcomes. Next, we assessed whether specific gene mutations were associated with subsequent CNS or systemic relapse. Activated B-cell-like genes such as PIM1, CD79B and PRDM1 were more frequently mutated in cases with CNS relapse, whereas TP53, BCL2 and FOXO1 mutations were associated with systemic relapse. Lastly, we performed tumor evolution analysis in five patients, identifying significant clonal divergence that accompanied CNS relapse. In summary, we identified both common biology underlying both CNS and systemic relapse, independent of cell-of-origin (COO), and specific COO-enriched gene mutations that conferred a higher risk of either CNS or systemic relapse. Moreover, we identified a pattern of clonal evolution that can explain difficulties in both predicting and treating CNS relapse.

  Study EGAS00001004057

• Preclinical Pediatric Molecular Analysis for Therapy Choice (MATCH)

  Purpose: Pediatric solid tumors arise from diverse tissues during development and exhibit a wide range of molecular, cellular and genetic features. This diversity, combined with the low incidence of pediatric cancer makes it increasingly difficult to personalize therapy for individual patients based on the unique features of their tumors. Therefore, well-credentialed preclinical models that capture the diversity and heterogeneity of pediatric solid tumors are essential for identifying molecular targeted therapeutics for precision medicine. Experimental Design: Here, we report 281 orthotopic patient derived xenografts (O-PDXs) from 224 patients representing 24 different types of pediatric solid tumors. We have performed genomic characterization of the O-PDXs and compared them to their corresponding patient tumors. To demonstrate the feasibility and utility of using such a diverse collection of O-PDXs in preclinical studies, we performed a preclinical pediatric precision medicine trial based on the NCI-COG Pediatric MATCH trial enrollment criteria. We also tested molecular targeted therapy for a novel oncogenic fusion recently reported in pediatric melanoma and precision drug delivery using nano-liposomal irinotecan. Results: In the preclinical match study, we identified one specific single agent response to LY3023414 in an osteosarcoma model with a TSC2 variant of unknown significance. Combinations of nano-liposomal irinotecan were effective in treating both Ewing sarcoma and rhabdomyosarcoma xenografts. Conclusion: Our studies demonstrate the value of large, well-credentialed preclinical models for future precision medicine in pediatric oncology using single agents, drug combinations and novel drug formulations.

  Study EGAS00001008011

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

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

  Study EGAS00001005737

• Massively parallel nanowell-based single-cell gene expression profiling

  Technological advances have enabled transcriptome characterization of cell types at the single-cell level providing new biological insights. New methods that enable simple yet high-throughput single-cell expression profiling are highly desirable. Here we report a novel nanowell-based single-cell RNA sequencing system, ICELL8, which enables processing of thousands of cells per sample. The system employs a 5184-nanowell-containing microchip to capture ~1300 single cells and process them. Each nanowell contains preprinted oligonucleotides encoding poly-d(T), a unique well barcode, and a unique molecular identifier. The ICELL8 system uses imaging software to identify nanowells containing viable single cells and only wells with single cells are processed into sequencing libraries. Here, we report the performance and utility of ICELL8 using samples of increasing complexity from cultured cells to mouse solid tissue samples. Our assessment of the system to discriminate between mixed human and mouse cells showed that ICELL8 has a low cell multiplet rate (< 3%) and low cross-cell contamination. We characterized single-cell transcriptomes of more than a thousand cultured human and mouse cells as well as 468 mouse pancreatic islets cells. We were able to identify distinct cell types in pancreatic islets, including alpha, beta, delta and gamma cells. Overall, ICELL8 provides efficient and cost-effective single-cell expression profiling of thousands of cells, allowing researchers to decipher single-cell transcriptomes within complex biological samples.

  Study EGAS00001002320

• Platelet RNAseq data for SLFN14 K219N patients

  Pathogenic missense variants in SLFN14, which encodes an RNA endoribonuclease protein that regulates rRNA degradation, are known to cause inherited thrombocytopenia with impaired platelet aggregation and ATP secretion. Despite rather mild laboratory defects, these patients display an obvious bleeding phenotype. The function of SLFN14 in megakaryocyte (MK) and platelet biology is unknown. This study aims to characterize the platelet transcriptome in patients with a SLFN14 K219N variant and model the disease in the immortalized megakaryocyte cell line imMKCL. Total platelet RNA was sequenced for two patients and 19 healthy controls. Differential gene expression analysis yielded a total of 2999 and 2888 significantly (|log2FC|>1, FDR<0.05) up- and downregulated genes, respectively. Remarkably, these downregulated genes were not enriched for any biological pathway while upregulated genes were enriched for pathways involved in (mitochondrial) translation and transcription with a significant upregulation of 134 ribosomal protein genes (RPG). Heterozygous and homozygous SLFN14 K219N imMKCL showed a defect in MK differentiation that aggravated during cell passaging of undifferentiated cells and in proplatelet formation. SLFN14 defective platelets and MK showed signs of rRNA degradation while this was absent in undifferentiated imMKCL cells. Upregulation of RPG through increased mTOR signaling in SLFN14 K219N MK seems to be a compensatory response for rRNA degradation. Indeed, mTOR inhibition with rapamycin resulted in further enhanced rRNA degradation in SLFN14 K219N MK. Taken together, our study indicates dysregulation of mTOR-regulated ribosomal biogenesis as the disease mechanism for SLFN14-related thrombocytopenia.

  Study EGAS00001006339

• Cisplatin increases sensitivity to FGFR inhibition in patient-derived xenograft models of lung squamous cell carcinoma

  Lung squamous cell carcinoma (SqCC) is a molecularly complex and genomically unstable disease. No targeted therapy is currently approved for lung SqCC, although potential oncogenic drivers of SqCC have been identified, including amplification of the fibroblast growth factor receptor 1 (FGFR1). Reports from a recently completed clinical trial indicate low response rates in patients treated with FGFR tyrosine kinase inhibitors, suggesting inadequacy of FGFR1 amplification as a biomarker of response, or the need for combination treatment. We aimed to develop accurate models of lung SqCC and determine improved targeted therapies for these tumors. We show that detection of FGFR1 mRNA by RNA in situ hybridisation is a better predictor of response to FGFR inhibition than FGFR1 gene amplification using clinically relevant patient-derived xenograft (PDX) models of lung SqCC. FGFR1-overexpressing tumors were observed in all histological subtypes of non-small cell lung cancers (NSCLC) as assessed on a tissue microarray, indicating a broader range of tumors that may respond to FGFR inhibitors. In FGFR1-overexpressing PDX tumors, we observed increased differentiation and reduced proliferation following FGFR inhibition. Combination therapy with cisplatin was able to increase tumor cell death, and dramatically prolonged animal survival compared to single agent treatment. Our data suggest that FGFR tyrosine kinase inhibitors can benefit NSCLC patients with FGFR1-overexpressing tumors and provides a rationale for clinical trials combining cisplatin with FGFR inhibitors.

  Study EGAS00001002423

• Spatio-temporal Profiling of a Rhabdoid Tumor Case Study

  Background Immune checkpoint inhibition (ICI) has revolutionized oncology, offering extended survival and long-term remission in previously incurable cancers. While highly effective in tumors with high mutational burden, lowly mutated cancers, including pediatric malignancies, present low response rate and limited predictive biomarkers. Patients and methods We present a framework for the identification and validation of tumor-reactive T cells as a biomarker to quantify ICI efficacy and as candidates for a personalized TCR-T cell therapy. Therefore, we profiled a pediatric malignant rhabdoid tumor patient with complete remission after ICI therapy using deep single-cell T cell receptor (TCR) repertoire sequencing of the tumor microenvironment (TME) and the peripheral blood. Results Tracking T cell dynamics longitudinally from the tumor to cells in circulation over a time course of 12 months revealed a systemic response and durable clonal expansion of tumor-resident and ICI-induced TCR clonotypes. We functionally validated tumor reactivity of TCRs identified from the TME and the blood by co-culturing patient-derived tumor cells with TCR-engineered autologous T cells. Here, we observed unexpectedly high frequencies of tumor-reactive TCR clonotypes in the TME and confirmed T cell dynamics in the blood post-ICI to predict tumor-reactivity. Conclusion These findings strongly support spatio-temporal tracking of T cell activity in response to ICI to inform therapy efficacy and to serve as a source of tumor-reactive TCRs for personalized TCR-T designs.

  Study EGAS00001008174

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

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

  Study EGAS00001005679

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

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

  Study EGAS00001005680

• The evolutionary history of human colitis-associated colorectal cancer

  Inflammatory bowel disease confers an increased lifetime risk of developing colorectal cancer (CRC), and colitis-associated CRC (CA-CRC) is molecularly distinct from sporadic CRC (S-CRC). Here we have dissected the evolutionary history of CA-CRC using multi-region sequencing. Exome sequencing was performed on fresh-frozen multiple regions of carcinoma, adjacent non-cancerous mucosa and blood from 12 CA-CRC patients (n=55 exomes), and key variants were validated with orthogonal methods. Genome-wide copy number profiling was performed using SNP arrays and low-pass whole genome sequencing on archival non-dysplastic mucosa (n=9), low-grade dysplasia (LGD, n=30), high-grade dysplasia (HGD, n=13), mixed LGD/HGD (n=7) and CA-CRC (n=19). Phylogenetic trees were reconstructed, and evolutionary analysis used to reveal the temporal sequence of events leading to CA-CRC. 10/12 tumors were microsatellite stable with a median mutation burden of 3.0 single nucleotide alterations (SNAs) per Mb, around 20% higher than S-CRC (2.5 SNAs/Mb), and consistent with elevated aging-associated mutational processes. Non-dysplastic mucosa had considerable mutation burden (median 47 SNAs), including mutations shared with the neighboring CA-CRC, indicating a pre-cancer mutational field. CA-CRCs were often near-triploid (42%) or near-tetraploid (21%) and phylogenetic analysis revealed that copy number alterations (CNAs) began to accrue in non-dysplastic bowel, but the LGD/HGD transition often involved a punctuated catastrophic CNA increase.

  Study EGAS00001003028

• ICU_transcriptomics__Assessing_the_role_of_the_host_immune_response_in_patients_with_ventilator_associated_pneumonia

  The aim of the proposed project is to provide transcriptomic sequencing data analysis of cells isolated from bronchioalveolar lavage and blood samples from critically ill patients with pneumonia on an intensive care unit in order to investigate the host contribution to disease. The lung microbiome of these patients will be analysed in parallel by deep sequencing and 16S sequencing. These data will provide important information to the development of strategies aimed at improving patient outcomes, reducing the development of antimicrobial resistance and facilitating infection control measures in hospital. Specifically, the proposed work will be a proof of principle study of 100 patients. We will use transcriptomic data to define the host contribution to the symptoms of these critically ill patients diagnosed with pneumonia. The data will be correlated with cytokine/chemokine concentrations, routine microbiology analysis and pathogen genomic data from the same samples to determine the ecology of the lungs of these patients. We will look for associations between clinical features, host response and pathogen genomics in order to identify new diagnostic markers. This study should provide important information to improve the speed and accuracy of diagnosis - resulting in more targeted and successful treatment of critically ill patients in the ICU. 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/

  Study EGAS00001003074

• Characterization of HCV-specific CD4 T cells during DAA-Therapy

  Background: Chronic HCV-infection is characterized by a severe impairment of HCV-specific CD4 T cell help that is driven by chronic antigen stimulation. We aimed to study the fate of HCV-specific CD4 T cells after viral elimination. Methods: HCV-specific CD4 T cell responses were longitudinally analyzed using MHC class II tetramer-technology, multicolor flow cytometry and RNA sequencing in a cohort of chronically HCV-infected patients undergoing therapy with direct-acting antivirals. In addition, HCV-specific neutralizing antibodies and CXCL13 levels were analyzed. Results: We observed that the frequency of HCV-specific CD4 T cells increased within two weeks after initiation of DAA therapy. Multicolor flow cytometry revealed a downregulation of exhaustion and activation markers and an upregulation of memory-associated markers. While cells with a Th1 phenotype were the predominant subset at baseline, cells with phenotypic and transcriptional characteristics of follicular T helper cells increasingly shaped the circulating HCV-specific CD4 T cell repertoire, suggesting antigen-independent survival of this subset. These changes were accompanied by a decline of HCV-specific neutralizing antibodies and the germinal center activity. Conclusion: We identified a population of HCV-specific CD4 T cells with a follicular T helper cell signature that is maintained after therapy-induced elimination of persistent infection and may constitute an important target population for vaccination efforts to prevent re-infection and immunotherapeutic approaches for persistent viral infections.

  Study EGAS00001003950

• Whole genome sequencing with linked reads of pediatric glioblastoma samples

  Pediatric glioblastoma (pGBM) is a lethal cancer with no effective therapies. To understand mechanisms of tumor evolution in this cancer, we performed whole genome sequencing with linked reads on longitudinally resected pGBM samples. Our analyses showed that all diagnostic and recurrent samples were collections of genetically diverse subclones. Clonal composition rapidly evolved at recurrence, with less than 8% of non-synonymous single nucleotide variants being shared in diagnostic-recurrent pairs. In order to track the origins of the mutational events we observed in pGBM, we generated whole genome datasets for two patients and their parents. These trios showed that genetic variants could be (i) somatic, (ii) inherited from a healthy parent, or (iii) arose de novo in the germlines of pGBM patients. Analysis of variant allele frequencies supported a model of tumor growth involving slow-cycling cancer stem cells that give rise to fast-proliferating progenitor-like cells and to non-dividing cells. Interestingly, radiation and anti-mitotic chemotherapeutics did not increase overall tumor burden upon recurrence. These findings support an important role for slow-cycling stem cell populations in contributing to recurrences, since slow-cycling cell populations are expected to be less prone to genotoxic stress induced by these treatments and to accumulate few mutations. Our results highlight the need for new targeted treatments that account for the complex functional hierarchies and genomic heterogeneity of pGBM.

  Study EGAS00001003432

• Probabilistic cell type assignment of single-cell transcriptomic data reveals spatiotemporal microenvironment dynamics in human cancers

  Single-cell RNA sequencing (scRNA-seq) has transformed biomedical research, enabling decomposition of complex tissues into disaggregated, functionally distinct cell types. For many applications, investigators wish to identify cell types with known marker genes. Typically, such cell type assignments are performed through unsupervised clustering followed by manual annotation based on these marker genes, or via "mapping" procedures to existing data. However, the manual interpretation required in the former case scales poorly to large datasets, which are also often prone to batch effects, while existing data for purified cell types must be available for the latter. Furthermore, unsupervised clustering can be error-prone, leading to under- and over- clustering of the cell types of interest. To overcome these issues we present CellAssign, a probabilistic model that leverages prior knowledge of cell type marker genes to annotate scRNA-seq data into pre-defined and de novo cell types. CellAssign automates the process of assigning cells in a highly scalable manner across large datasets while simultaneously controlling for batch and patient effects. We demonstrate the analytical advantages of CellAssign through extensive simulations and exemplify real-world utility to profile the spatial dynamics of high-grade serous ovarian cancer and the temporal dynamics of follicular lymphoma. Our analysis reveals subclonal malignant phenotypes and points towards an evolutionary interplay between immune and cancer cell populations with cancer cells escaping immune recognition.

  Study EGAS00001003452

• Bioinformatic Methods and Bridging of Assay Results for Reliable Tumor Mutational Burden Assessment in Non-Small Cell Lung Cancer

  Tumor mutational burden (TMB) has emerged as a predictive biomarker of response to immune checkpoint inhibitors. Standardization of TMB measurement is essential for implementing diagnostic tools to guide treatment. Here we evaluate bioinformatic TMB analysis by whole exome sequencing (WES) in formalin-fixed, paraffin-embedded samples. In CheckMate 026, TMB was retrospectively assessed in 312 patients with non-small cell lung cancer (58% of the intent-to-treat population) who received first-line nivolumab treatment or chemotherapy. We examined the sensitivity of TMB assessment to bioinformatic filtering methods and assessed concordance between TMB data derived by WES and the FoundationOne CDx™ assay. TMB scores comprising synonymous, indel, frameshift, and nonsense mutations (all mutations) were 3.1-fold higher than data including missense mutations only, but values were highly correlated (Spearman’s r = 0.99). Scores from CheckMate 026 samples including missense mutations only were similar to those generated from data in The Cancer Genome Atlas, but those including all mutations were generally higher. Using databases for germline subtraction (instead of matched controls) showed a trend for race-dependent increases in TMB scores. Parameter variation can therefore impact TMB calculations, highlighting the need for standardization. Encouragingly, WES and FoundationOne CDx outputs were highly correlated (Spearman’s r = 0.90) and differences could be accounted for by empirical calibration, suggesting that reliable TMB assessment across assays, platforms and centers is achievable.

  Study EGAS00001003661

• Whole exome sequencing of an invasive diffuse intrinsic pontine glioma sampled from different sites

  The lateral ventricle subventricular zone (SVZ) is a frequent site of high-grade glioma (HGG) spread and tumor invasion of the SVZ is a predictor of a worse clinical prognosis. A range of HGG types invade the stem cell niche of the SVZ, including both adult glioblastoma and pediatric high-grade gliomas such as diffuse intrinsic pontine glioma (DIPG). The cellular and molecular mechanisms mediating this frequent invasion of the SVZ are poorly understood. Here we demonstrate that neural precursor cell (NPC):glioma cell communication underpins the propensity of glioma to colonize the SVZ. SVZ NPCs secrete chemoattractant signals toward which glioma cells home. Biochemical and proteomic analyses of factors secreted by SVZ NPCs revealed a discrete list of candidate proteins. Necessity and sufficiency testing implicated the neurite outgrowth-promoting factor pleiotrophin, along with three required binding partners (secreted protein acidic and rich in cysteine (SPARC), SPARC-like protein 1, and heat shock protein 90B), as key mediators of this chemoattractant effect. Pleiotrophin protein expression is strongly enriched in the SVZ, and knockdown of pleiotrophin expression starkly reduced glioma invasion of the SVZ in the adult murine brain. Pleiotrophin, in complex with the three binding partners, activated the Rho/ROCK pathway in DIPG cells, and inhibition of ROCK resulted in decreased DIPG invasion toward SVZ neural precursor cell-secreted factors. These findings demonstrate a pathogenic role for neural precursor cell – glioma interactions and potential therapeutic targets to limit glioma invasion.

  Study EGAS00001002326

• Identification of a Fumarate-Hydratase Deficient-like RCC Subtype with Convergent Pathway Alterations

  The development of molecularly targeted therapies builds upon the identification of tumor clades with shared molecular abnormalities. Hereditary leiomyomatosis and renal cell carcinoma (HLRCC) is a familial syndrome resulting from germline mutation of the two-hit tumor suppressor gene, Fumarate Hydratase (FH). HLRCC renal tumors present often in young adults, exhibit unique histological features and are particularly aggressive. HLRCC is uncommon and poorly understood, and whether there are sporadic renal tumors with similar characteristics, is unclear. Notably, FH can also be inactivated somatically in RCC, and should be included among FH-deficient (FHD) kidney tumors. We identified a tumor clade without FH mutations but with similar morphological and molecular characteristics as well as aggressive clinical behavior, which we name FHDL (FHD-like). We show that FHD/FHDL tumors represent a unique subtype of CIMP (CpG Island Methylator Phenotype) with the highest levels of overall genome methylation and CDKN2A silencing. FHD/FHDL tumors are characterized by convergent activation of the Hippo pathway as determined by frequent and mutually exclusive mutations in NF2 and PTPN14. AKR1B10, an NRF2 target, was upregulated in FHD/FHDL tumors as well as GATA3, which suggests a shared origin in the renal medulla and is consistent with topological studies. There is loss of 5-hydroxymethylcytosine in both FHD and FHDL tumors. These data define a group of particularly aggressive renal cancers with shared molecular features and potential targets for therapeutic intervention.

  Study EGAS00001002646

• Pseudodiastrophic dysplasia expands the known phenotypic spectrum of defects in proteoglycan biosynthesis

  Pseudodiastrophic dysplasia (PDD) is a severe skeletal dysplasia associated with prenatal manifestation and early lethality. Clinically, PDD is classified as a ‘dysplasia with multiple joint dislocations’ however, the molecular aetiology of the disorder is currently unknown. In order to identify the genetic defects underlying PDD, we performed whole exome sequencing (WES) on 3 patients from 2 unrelated families, clinically diagnosed with PDD. WES resulted in the identification of bi-allelic variants in the established skeletal dysplasia genes, B3GAT3 (Family 1) and CANT1 (Family 2). Mutations in these genes have previously been reported to cause ‘multiple joint dislocations, short stature, and craniofacial dysmorphism with or without congenital heart defects’ (B3GAT3) and Desbuquois dysplasia 1 (CANT1); disorders in the same nosological group as PDD. Follow-up of the B3GAT3 variants demonstrated significantly reduced B3GAT3/GlcAT-I expression. Downstream in vitro functional analysis revealed abolished biosynthesis of glycosaminoglycan side chains on proteoglycans. Functional evaluation of the CANT1 variant showed impaired nucleotidase activity, which results in inhibition of glycosaminoglycan synthesis through accumulation of uridine diphosphate. For the families described in this study, the PDD phenotype was caused by mutations in the known skeletal dysplasia genes B3GAT3 and CANT1, demonstrating the advantage of genomic analyses in delineating the molecular diagnosis of skeletal dysplasias. This finding expands the phenotypic spectrum of B3GAT3- and CANT1- related skeletal dysplasias to include PDD, and highlights the significant phenotypic overlap of conditions within the proteoglycan biosynthesis pathway.

  Study EGAS00001004145

• Successful BRAF/MEK-inhibition in a young patient with BRAF V600E-mutated extrapancreatic acinar cell carcinoma (HIPO-021)

  Pancreatic acinar cell carcinoma (PAC) is a rare disease with a poor prognosis. Treatment options for metastatic PAC are limited and often follow chemotherapeutic regimens for pancreatic ductal adenocarcinoma. Although recurrent genomic alterations, such as BRAF fusions and defects in genes involved in homologous recombination DNA repair, have been described in PAC, data on the clinical efficacy of molecularly guided, targeted treatment are scarce. Here we describe the case of a 27-yr-old patient with BRAF V600E-mutated PAC who was successfully treated with a combination of BRAF and MEK inhibitors. The patient presented to our clinic with abdominal pain and weight loss. Imaging showed extensive retroperitoneal disease as well as mediastinal lymphadenopathy. Because of elevated alpha-fetoprotein (AFP) levels and inconclusive histologic findings, a germ cell tumor was suspected; however, PEI chemotherapy was unsuccessful. A repeat biopsy yielded the diagnosis of PAC and treatment with FOLFIRINOX was initiated. Comprehensive molecular profiling within the MASTER (Molecularly Aided Stratification for Tumor Eradication Research) precision oncology program revealed a somatic BRAF V600E mutation and a germline PALB2 stop-gain mutation. Therapy was therefore switched to BRAF/MEK inhibition, resulting in almost complete remission and disease control for 12 mo and a remarkable improvement in the patient's general condition. These results indicate that BRAF alterations are a valid therapeutic target in PAC that should be routinely assessed in this patient population.

  Study EGAS00001004282

• H3.3G34R/V-transformed interneuron progenitors co-opt PDGFRA for gliomagenesis

  Histone H3.3 glycine 34 to arginine/valine (H3.3G34R/V) mutations occur in deadly hemispheric high-grade gliomas. These tumors show exquisite regional and temporal specificity, suggesting a developmental context permissive to the effects of G34R/V mutations. Here we present the molecular landscape of G34R/V gliomas (n=85) and show that 50% bear activating mutations in PDGFRA, with strong selection pressure for PDGFRAMUT clones at recurrence. We show that G34R/V tumors arise in interneuron progenitors of the foetal ventral forebrain expressing GSX2 and the DLX family of homeobox transcription factors, where terminal neuronal differentiation is impaired through aberrant G34R/V-mediated H3K27me3. Frequent co-occurrence of G34R/V & PDGFRAMUT is facilitated in this interneuron lineage-of-origin as PDGFRA forms an aberrant chromatin loop with the adjacent GSX2, hijacking its active chromatin conformation. At the single-cell level, G34R/V tumours entirely lack oligodendroglial transcriptional programs prominent in other glioma entities, and instead harbour dual neuronal and astroglial compartments. CRISPR-removal of H3.3G34R/V does not impact tumorigenicity suggesting this mutation becomes dispensable, while PDGFRAMUT are potently oncogenic regardless of G34 mutation. Collectively, our results suggest that G34R/V gliomas arise in foetal interneuron progenitors unable to undergo terminal differentiation, enabling co-option of PDGFRA through inappropriate expression and activating mutations to promote gliogenesis and oncogenicity. Reliance on PDGFRA for oncogenesis may be of therapeutic opportunity in G34R/V glioma.

  Study EGAS00001004301

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

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

  Study EGAS00001005678

• A no-stop mutation in MAGEB4 is a possible cause of rare X-linked azoospermia and oligozoospermia in a consanguineous Turkish family

  Purpose: The purpose of this study was to identify mutations that cause non-syndromic male infertility using whole exome sequencing of family cases.Methods: We recruited a consanguineous Turkish family comprising nine siblings with male triplets; two of the triplets were infertile as well as one younger infertile brother. Whole exome sequencing (WES) performed on two azoospermic brothers, identified a mutation in the MAGEB4 gene which was confirmed via Sanger sequencing and then screened for on control groups and unrelated infertile subjects. The effect of the mutation on mRNA and protein levels were tested after in vitro cell transfection. Structural features of MAGEB4 were predicted throughout the conserved MAGE domain.Results: The novel single base substitution (c.1041A>T) in the X-linked MAGEB4 gene was identified as a no-stop mutation. The mutation is predicted to add 24 amino acids to the C-terminus of MAGEB4. Our functional studies were unable to detect any effect either on mRNA stability, intracellular localization of the protein or the ability to homo / heterodimerize with other MAGE proteins. We thus hypothesize that these additional amino acids may affect the proper protein interactions with MAGEB4 partners. Conclusion: The whole exome analysis of a consanguineous Turkish family revealed MAGEB4 as a possible new X-linked cause of inherited male infertility. This study provides the first clue to the physiological function of a MAGE protein.

  Study EGAS00001002272

• Performance comparison of three DNA extraction kits on human whole-exome data from formalin-fixed paraffin-embedded normal and tumor samples

  We generated 42 human whole-exome sequencing data sets from fresh-frozen (FF) and FFPE samples. These samples include normal and tumor tissues from two different organs (liver and colon), that we extracted with three different FFPE extraction kits (QIAamp DNA FFPE Tissue kit and GeneRead DNA FFPE kit from Qiagen, Maxwell\textsuperscript{TM} RSC DNA FFPE Kit from Promega). Variant calling analysis shows a very high rate of concordance between matched FF / FFPE pairs and equivalent performance for the three kits we analyzed. We find a significant variation in the difference of total number of variants called between FF and FFPE samples for the three different FFPE DNA extraction kits. Coverage analysis shows that FFPE samples have less good indicators than FF samples, yet the coverage quality remains above accepted thresholds. We detect limited but significant variations in coverage indicator values between the three FFPE extraction kits. Globally, the GeneRead and QIAamp kits have better variant calling and coverage indicators than the Maxwell kit on the samples used in this study, although this kit performs better on some indicators and has advantages in terms of practical usage. Taken together, our results confirm the potential of FFPE samples analysis for clinical genomic studies, but also indicate that the choice of a FFPE DNA extraction kit should be done with careful testing and analysis beforehand in order to maximize the accuracy of the results.

  Study EGAS00001002631

• Genome-wide analysis of genetic risk factors for rheumatic heart disease in Aboriginal Australians provides support for pathogenic molecular mimicry

  Background. Rheumatic heart disease (RHD) following Group A Streptococcus (GAS) infections is heritable and prevalent in Indigenous populations. Molecular mimicry between human and GAS proteins triggers pro-inflammatory cardiac valve-reactive T-cells. Methods. Genome-wide genetic analysis was undertaken in 1263 Aboriginal Australians (398 RHD cases; 865 controls). Single nucleotide polymorphisms (SNPs) were genotyped using Illumina HumanCoreExome BeadChips. Direct typing and imputation was used to fine-map the human leukocyte antigen (HLA) region. Epitope binding affinities were mapped for human cross-reactive GAS proteins, including M5 and M6. Results. The strongest genetic association was intronic to HLA-DQA1 (rs9272622; P=1.86x10-7). Conditional analyses showed rs9272622 and/or DQA1*AA16 account for the HLA signal. HLA-DQA1*0101_DQB1*0503 (OR 1.44, 95%CI 1.09-1.90, P=9.56x10-3) and HLA-DQA1*0103_DQB1*0601 (OR 1.27, 95%CI 1.07-1.52, P=7.15x10-3) were risk haplotypes; HLA_DQA1*0301-DQB1*0402 (OR 0.30, 95%CI 0.14-0.65, P=2.36x10-3) was protective. Human myosin cross-reactive N-terminal and B repeat epitopes of GAS M5/M6 bind with higher affinity to DQA1/DQB1 alpha/beta dimers for the two risk haplotypes than the protective haplotype. Conclusions. Variation at HLA_DQA1-DQB1 is the major genetic risk factor for RHD in Aboriginal Australians studied here. Cross-reactive epitopes bind with higher affinity to alpha/beta dimers formed by risk haplotypes, supporting molecular mimicry as the key mechanism of RHD pathogenesis.

  Study EGAS00001002678

• Germline variants collaborate with somatic mutations and initiate and/or drive disease in primary myelodysplastic syndrome (MDS) and therapy-related myeloid neoplasms (t-MN)

  Therapy-related myeloid neoplasms (T-MN) are poorly characterized secondary hematological malignancies following chemotherapy/radiotherapy exposure. We compared the clinical and mutational characteristics of T-MN (n=129) and primary myelodysplastic syndrome (P-MDS, n=108) patients. Although the somatic mutation frequency was similar between T-MN and P-MDS patients (93% in both groups), the pattern was distinct. TP53 mutations were more frequent in T-MN (29.5% versus 7%), while spliceosomal complex mutations were more common in P-MDS (56.5% versus 25.6%).In contrast to P-MDS, the ring sideroblasts (RS) phenotype was not associated with better survival in T-MN, most probably due to genetic association with TP53 mutations. SF3B1 was mutated in 96% of P-MDS with ≥15% RS, but in only 32% T-MN. TP53 mutations were detected in 92% T-MN with ≥15% RS and SF3B1 wild-type cases. Interestingly, T-MN and P-MDS patients with “Very low” or “Low” Revised International Prognostic Scoring System (IPSS-R) showed similar biological and clinical characteristics. In a Cox-regression analysis, TP53 mutation was a poor prognostic factor in T-MN, independent of IPSS-R cytogenetic, disease-modifying therapy and NRAS mutation.Our data have direct implications for T-MN management and provide evidence that, in addition to conventional disease parameters, mutational analysis should be incorporated in T-MN risk stratification.

  Study EGAS00001003547

• Whole Tumor spatial heterogeneity in metastatic melanoma

  Heterogeneous inter-tumoral responses, sustained periods of apparent clinical benefit despite lack of objective response to various therapies, and even spontaneous remission are well known within a subpopulation of advanced melanoma patients. The molecular and cellular dynamics facilitating long-term survival remain poorly defined, particularly in the current era ofexposure to multiple potentially active therapies. We studied an exceptional case of long-term survival in a patient with non-responding metastatic melanoma in order to characterize the clonal and microenvironmental factors.We performed immunogenomic analysis of 3 metachronous tumors, and of 67 intratumor sub-regions of one metastatic lesion, usingwhole exome sequencing, RNA-sequencing, immunohistochemistry, and T cell receptor sequencing. Longitudinal analysesidentified mutations in severalknownresistance genes affecting distinctmetastases. Intra-tumoral immune signatures revealed spatially-distinct pockets of immune activation and suppression throughout the metastasis. A specific T cell VBeta CDR3 rearrangement was identified as both dominant and recurrent not only across multiple spatial points within a single tumor mass, but also across metachronous tumors spanning the patient’s disease course. Immunophenotyping of this T cell population with single-cell sequencingsuggestedrepeatedT-cell priming events leading to the persistence of both activated and exhausted T cells bearing this same TCR-Beta at any given time.Our findings highlight genomic and immune heterogeneityin metastatic disease, extending to the intra-tumoral leveland provide a rationale for integrative spatiotemporal analysis of metastatic lesions in the context of immune checkpoint blockade.

  Study EGAS00001003292

• The spatio-temporal evolution of lymph node spread in early breast cancer

  Purpose: The most significant prognostic factor in early breast cancer is lymph node involvement. This stage between localised and systemic disease is key to understanding breast cancer progression, however our knowledge of the evolution of lymph node malignant invasion remains limited, as most currently available data derive from primary tumours. Experimental design: In 11 treatment-naïve node positive early breast cancer patients without clinical evidence of distant metastasis, we investigated lymph node evolution using spatial multi-region sequencing (n=78 samples) of primary and lymph node deposits and genomic profiling of matched longitudinal circulating tumour DNA (ctDNA). Results: Linear evolution from primary to lymph node was rare (1/11) whereas the majority of cases displayed either early divergence between primary and nodes (4/11), or no detectable divergence (6/11) where both primary and nodal cells belonged to a single recent expansion of a metastatic clone. Divergence of metastatic subclones was driven in part by APOBEC. Longitudinal ctDNA samples from 2 subjects taken peri-operatively reflected the two major evolutionary patterns and demonstrate that private mutations can be detected even from early metastatic nodal deposits. Moreover, node removal resulted in disappearance of private lymph node mutations in ctDNA.Conclusions: This study sheds new light on a crucial evolutionary step in the natural history of breast cancer, demonstrating early establishment of axillary lymph node metastasis in a substantial proportion of patients.

  Study EGAS00001002947

• GCAT | Genomes for life: cohort study of the genomes of Catalonia

  The GCAT Study have recruited 20 000 participants aged 40–65 years. Participants who agreed to take part in the study completed a self-administered computer-driven questionnaire, and underwent blood pressure, cardiac frequency and anthropometry measurements. For each participant, blood plasma, blood serum and white blood cells are collected at baseline. A total of 5459 genomic profiles have been characterised by comprehensive genotyping. Genome-wide genotypes have been generated using Illumina Infinium SNP-bead array technology. We chose the Multi-Ethnic Global (MEGAEX, V.2) consortium array, a multipurpose, multiethnic genotyping array with two million selected markers (including previously described germline mutations, insertions-deletions (InDels) and SNPs).We have strictly followed the standard manufacturer recommended automated protocol for the Infinium HTS Assay scanned with a HiScan confocal scanner (Illumina, San Diego, California, USA). Genome Studio V.2011.1 has been used for raw data analysis. Genotyping was performed at the Genomics and Bioinformatics Unit of the PMPPC Institute for Health Science Research Germans Trias i Pujol, in Badalona, Spain. Future plans The first follow-up study started in December 2017 and will end by March 2018. Residences of all subjects will be geocoded during the following year. Several genomic analyses are ongoing, and metabolomic and genomic integrations will be performed to identify underlying genetic variants, as well as environmental factors that influence metabolites. http://dx.doi.org/10.1136/bmjopen-2017-018324

  Study EGAS00001003018

• Open_Targets_020_Epigenomes_of_Cell_Lines

  There is currently a drive to establish cell based assay systems of greater human biological and disease relevance through the use of well characterised transformed cell lines, primary cells and complex cellular models (e.g. co-culture, 3D models). However, although the field is gaining valuable experience in running more non-standard & complex cell assays for target validation and compound pharmacology studies, there is the lack of a systematic approach to determine if this expansion in cell assay models is reflected in increased human biological and disease relevance. The increasing wealth of publically available transcriptomic, and epigenome (ENCODE and Epigenome Roadmap) data represents an ideal reference mechanism for determining the relationship between cell types used for target & compound studies to primary human cells and tissues from both healthy volunteers & patients. The CTTV020 epigenomes of cell line project aims to generate epigenetic and transcriptomic profiles of cell lines and compare these with existing and newly generated reference data sets from human tissue and cell types. The aim is to identify assay systems which will provide greater confidence in translating target biology and compound pharmacology to patients. Multiple cell types commonly used within research have been grouped according to biology. Examples include erythroid, lung epithelial, hepatocyte cell types and immortalised models of monocyte / macrophage biology. 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 EGAS00001003138

• Cistrome-partitioning reveals convergence of somatic mutations and risk-variants on master transcription regulators in primary prostate tumors

  Thousands of noncoding somatic Single Nucleotide Variants (SNVs) of unknown function are reported in tumors. Partitioning the genome according to cistromes, reveals the enrichment of somatic SNVs in prostate tumors as opposed to adjacent normal tissue cistromes of master transcription regulators, including AR, FOXA1 and HOXB13. This parallels enrichment of prostate cancer genetic predispositions over these transcription regulators’ tumor cistromes, exemplified at the 8q24 locus harboring both risk-variants and somatic SNVs in cis-regulatory elements, upregulating MYC expression and altering the binding of transcription regulators to DNA. However, Massively-Parallel Reporter Assays reveal that few SNVs can alter the transactivation potential of individual CREs. Instead, SNVs accumulate, similarly to inherited riskvariants, in cistromes of master transcription regulators required for prostate cancer development. Difficulties in inferring the biological significance of noncoding mutations have limited their inclusion in precision genomics medicine pipelines. Most attempts to delineate a role for noncoding mutations relied on detecting evidence for positive selection within individual CREs, such as reported for the TERT gene promoter. By considering the enrichment of noncoding mutations in cistromes as opposed to individual CREs, we reveal their specificity towards master transcription regulators that promote prostate cancer development, a feature shared with inherited risk-variants. Overall, our work provides a blueprint for the functional interpretation of noncoding mutations in genomic tests relying on defining cis-regulatory units according to cistrome-partitioning to identify cancer driver transcription regulators.

  Study EGAS00001003928

• Neoadjuvant combination PD-L1 plus CTLA-4 blockade in patients with cisplatin-ineligible operable urothelial carcinoma

  Immune checkpoint therapy (ICT) is being tested in the neoadjuvant setting for patients with localized urothelial carcinoma (UC), with one study reporting data in cisplatin-ineligible patients who received anti PD-L1 monotherapy. The study reported that patients with bulky tumors, a known high-risk feature defined as greater than clinical T2 disease, had fewer responses, with pathologic complete response (pCR) rate of 17%. Here, we report on the first pilot combination neoadjuvant trial (NCT02812420) with anti-PD-L1 (durvalumab) plus anti-CTLA-4 (tremelimumab) in cisplatin-ineligible patients, with all tumors identified as having high risk features (N=28). Primary endpoint was safety and we observed 6 of 28 patients (21%) with grade ≥3 immune-related adverse events, consisting of asymptomatic laboratory abnormalities (N=4), hepatitis and colitis (N=2). We also observed pCR of 37.5% and downstaging to pT1 or less in 58% of patients who completed surgery (N=24). In summary, we provide initial safety, efficacy and biomarker data with neoadjuvant combination anti-PD-L1 plus anti-CTLA-4, which warrants further development for patients with localized UC, especially cisplatin-ineligible patients with high-risk features who do not currently have an established standard-of-care neoadjuvant treatment. Publication: Gao et al. Neoadjuvant PD-L1 plus CTLA-4 blockade in patients with cisplatin-ineligible operable high-risk urothelial carcinoma. Nature Medicine volume 26, pages1845–1851(2020)

  Study EGAS00001004074

• Exploiting evolutionary steering in cancer therapy

  Drug resistance, mediated by intra-tumour heterogeneity and clonal evolution, is arguably the biggest problem in cancer therapy today. However, evolving resistance to one drug may come at the cost of increased sensitivity to another due to so-called evolutionary trade-offs. This weakness can be exploited in the clinic using an approach called ‘evolutionary herding’ that aims at controlling the tumour cell population to delay or prevent resistance. However, model systems able to recapitulate cancer evolution experimentally are lacking and current in vitro techniques based on small populations, re-plating and escalating dose, are unsuitable to develop evolutionary herding strategies. We present a novel methodology for evolutionary herding in vitro and ex vivo using patient-derived organoids. Our approach is based on a combination of single-cell barcoding, very large populations of 108-109 cells grown without re-plating, realistic high drug doses, time-course monitoring of cancer clones, and mathematical modelling of tumour evolutionary dynamics. We demonstrate evolutionary herding in non-small cell lung cancer in vitro and in patient-derived colorectal cancer organoids (PDO). We show that herding causes controlled evolutionary bottlenecks that lead to collateral sensitivity. Through genomic analysis, we were also able to determine the mechanisms that drive such sensitivity. Our approach allows modelling evolutionary trade-offs experimentally to test patient-specific evolutionary herding strategies that can be translated into the clinic to control treatment resistance.

  Study EGAS00001003200

• Genome-wide study of the effect of blood collection tubes on the cell-free DNA methylome

  Background The methylation pattern of cfDNA, isolated from liquid biopsies, is gaining substantial interest for diagnosis and monitoring of diseases. We have evaluated the impact of type of blood collection tube and time delay between blood draw and plasma preparation on bisulfite-based cfDNA methylation profiling. Methods 15 tubes of blood were drawn from three healthy volunteer subjects (BD Vacutainer K2E EDTA spray tubes, Streck Cell-Free DNA BCT tubes, PAXgene Blood ccfDNA tubes, Roche Cell-Free DNA Collection tubes and Biomatrica LBgard blood tubes in triplicate). Samples were either immediately processed or stored at room temperature for 24 or 72 hours before plasma preparation. DNA fragment size was evaluated by capillary electrophoresis. Reduced representation bisulfite sequencing was performed on the cell-free DNA isolated from these plasma samples. We evaluated the impact of blood tube and time delay on several quality control metrics. Results All preservation tubes performed similar on the quality metrics that were evaluated. Furthermore, a considerable increase in cfDNA concentration and the fraction of it derived from NK cells was observed after a 72-hour time delay in EDTA tubes. Conclusion The methylation pattern of cfDNA is robust and reproducible in between the different preservation tubes. EDTA tubes processed as soon as possible, preferably within 24 hours, are the most cost effective. If immediate processing is not possible, preservation tubes are valid alternatives.

  Study EGAS00001004271

• Long-read trio sequencing of unsolved patients with intellectual disability

  Long-read sequencing (LRS) has the potential to comprehensively identify all medically-relevant genome variation, including variation commonly missed by short read sequencing (SRS) approaches. To determine this potential, we performed LRS up to 40x coverage using the Pacific Biosciences Sequel instrument for 5 patient-parent trios with intellectual disability, whose aetiology remained unresolved after SRS exomes and genomes. A first systematic assessment of LRS coverage showed that ~35Mb of the genome was only accessible by LRS and not SRS. Structural variant (SV) calling yielded on average 28,292 SV calls per individual, totalling 12.9 MB of sequence. Trio-based analyses showed concordance for up to 95% of these SV calls across the genome, and 80% of the LRS SV calls were not identified by SRS. De novo mutation analysis did not identify any de novo SVs, confirming these are rare events. Because of high sequence coverage, we were also able to call single nucleotide substitutions. On average, we identified 3 million substitutions per genome, with a Mendelian inheritance concordance of up to 97%. Of these, approximately 100,000 were located in the ~35Mb of the genome that was only captured by LRS. Moreover, these variants affected the coding sequence of 64 genes, including 32 known to cause Mendelian disorders. Although no clinical diagnoses were made in these 5 patients, our data show the potential added value of LRS compared to SRS for identifying medically relevant genome variation.

  Study EGAS00001004319

• Persistent STAG2 mutation in recurrent pediatric glioblastoma

  Similar to their adult counterparts, the prognosis for pediatric patients with high-grade gliomas remains poor. At time of recurrence, treatment options are limited and remain without consensus. This report describes the genetic findings, obtained from whole exome sequencing of a pediatric patient with glioblastoma who underwent multiple surgical resections and treatment with standard chemoradiation, as well as a novel recombinant poliovirus vaccine therapy. Strikingly, despite the variety of treatments, there was persistence of a tumor clone, characterized by a deleterious STAG2 mutation, whose deficiency in pre-clinical studies can cause aneuploidy and aberrant mitotic progression, but remains understudied in the clinical setting. There was near elimination of an EGFR mutated and amplified tumor clone after gross total resection, standard chemoradiation, and poliovirus therapy, followed by the emergence of a persistently STAG2 mutated clone, with rare mutations in PTPN11 and BRAF, the latter composed of a novel deleterious mutation previously not reported in pediatric glioblastoma (p.D594G). This was accompanied by a mutation signature shift towards one characterized by increased DNA damage repair defects, consistent with the known underlying STAG2 deficiency. As such, this case represents a novel report following the clinical and genetic progression of a STAG2 mutated glioblastoma, including treatment with a novel and emerging immunotherapy. Although STAG2 deficiency comprises only a small subset of gliomas, this case adds clinical evidence to existing pre-clinical data supporting a role for STAG2 mutations in tumor formation and resistance to standard therapies.

  Study EGAS00001004340

• The Medical Genome Reference Bank: a whole genome data resource of 4,000 healthy elderly individuals.

  Allele frequency data from human reference populations is of increasing value for filtering and assignment of pathogenicity to genetic variants. Aged and healthy populations are more likely to be selectively depleted of pathogenic alleles, and therefore particularly suitable as a reference populations for the major diseases of clinical and public health importance. However, reference studies of the healthy elderly have remained under represented in human genetics. We have developed the Medical Genome Reference Bank (MGRB), a large scale comprehensive whole genome dataset of confirmed healthy elderly individuals, to provide a publicly accessible resource for health related research, and for clinical genetics. It also represents a useful resource for studying the genetics of healthy aging. The MGRB comprises 4,000 healthy, older individuals with no reported history of cancer, cardiovascular disease or dementia, recruited from two Australian community based cohorts. DNA derived from blood samples will be subject to whole genome sequencing. The MGRB will measure genome wide genetic variation in 4,000 individuals, mostly of European decent, aged 60 to 95 years (mean age 75 years). The MGRB has committed to a policy of data sharing, employing a hierarchical data management system to maintain participant privacy and confidentiality, whilst maximizing research and clinical usage of the database. The MGRB will represent a dataset of international significance, broadly accessible to the clinical and genetic research community.

  Study EGAS00001003511

• Initial whole genome sequencing of plasma cell neoplasms in First Responders exposed to the World Trade Center attack of September 11, 2001

  The World Trade Center (WTC) attack of September 11, 2001 created an unprecedented environmental exposure to known and suspected carcinogens. High incidence of multiple myeloma (MM) and precursor conditions has been reported among first responders to the WTC disaster. To expand on our prior screening studies, and to characterize the genomic impact of the exposure to known and potential carcinogens in the WTC debris, we were motivated to perform whole genome sequencing (WGS) of WTC first responders and recovery workers who developed a plasma cell disorder after the attack. PATIENTS AND METHODS: We performed WGS of 9 CD138-positive bone marrow mononuclear samples from patients who were diagnosed with plasma cell disorders after the WTC disaster. RESULTS: No significant differences were observed in comparing the post-WTC driver and mutational signatures landscape with 110 previously published WGS from 56 patients with MM and the CoMMpass WGS cohort (n=752). Leveraging constant activity of the single base substitution mutational signatures 1 and 5 over time, we estimated that tumor-initiating chromosomal gains were windowed to both pre- and post-WTC exposure. CONCLUSIONS: Although limitations in sample size preclude any definitive conclusions, our findings suggest that the observed increased incidence of plasma cell neoplasms in this population is due to complex and heterogeneous effects of the WTC exposure that may have initiated or contributed to progression of malignancy.

  Study EGAS00001004467

• Targets of MEK inhibition in DIPG

  The survival of children with diffuse intrinsic pontine glioma (DIPG) remains dismal, with new treatments desperately needed. In the era of precision medicine, targeted therapies represent an exciting treatment opportunity, yet resistance can rapidly emerge, playing an important role in treatment failure. In a prospective biopsy-stratified clinical trial, we combined detailed molecular profiling (methylation BeadArray, exome, RNAseq, phospho-proteomics) linked to drug screening in newly-established patient-derived models of DIPG in vitro and in vivo. We identified a high degree of in vitro sensitivity to the MEK inhibitor trametinib in samples which harboured genetic alterations targeting the MAPK pathway, including the non-canonical BRAF_G469V mutation, and those affecting PIK3R1 and NF1. However, treatment of PDX models and of a patient with trametinib at relapse failed to elicit a significant response. We generated trametinib-resistant clones in the BRAF_G469V model through continuous drug exposure, and identified acquired mutations in MEK1/2 (MEK1_K57N, MEK1_I141S and MEK2_I115N) with sustained pathway up-regulation. These cells showed the hallmarks of mesenchymal transition, and expression signatures overlapping with inherently trametinib-insensitive primary patient-derived cells that predicted an observed sensitivity to dasatinib. Combinations of trametinib with dasatinib and the downstream ERK inhibitor ulixertinib showed highly synergistic effects in vitro. These data highlight the MAPK pathway as a therapeutic target in DIPG, and show the importance of parallel resistance modelling and rational combinatorial treatments likely to be required for meaningful clinical translation.

  Study EGAS00001004495

• RRBS sequencing of 7 tumour regions and a normal sample from a single TRACERx patient.

  Background: Epigenetic heterogeneity within a tumour can play an important role in tumour evolution and the emergence of resistance to treatment. It is increasingly recognised that the study of DNA methylation (DNAm) patterns along the genome - so-called `epialleles' - offers greater insight into epigenetic dynamics than conventional analyses which examine DNAm marks individually. Results: We have developed a Bayesian model to infer which epialleles are present in multiple regions of the same tumour. We apply our method to reduced representation bisulfite sequencing (RRBS) data from multiple regions of one lung cancer tumour and a matched normal sample. The model borrows information from all tumour regions to leverage greater statistical power. The total number of epialleles, the epiallele DNAm patterns, and a noise hyperparameter are all automatically inferred from the data. Uncertainty as to which epiallele an observed sequencing read originated from is explicitly incorporated by marginalising over the appropriate posterior densities. The degree to which tumour samples are contaminated with normal tissue can be estimated and corrected for. By tracing the distribution of epialleles throughout the tumour we can infer the phylogenetic history of the tumour, identify epialleles that differ between normal and cancer tissue, and define a measure of global epigenetic disorder. Conclusions: Detection and comparison of epialleles within multiple tumour regions enables phylogenetic analyses, identification of differentially expressed epialleles, and provides a measure of epigenetic heterogeneity.

  Study EGAS00001002484

• Renal_habitat_WXS

  This project aims at correlating genetic diversity (as well as transcriptomic, metabolic and immune diversity) observed in multiple targeted tumour biopsies and circulating tumour DNA with heterogeneity on multiparametric MRI. To this end, 20 patients with RCC (15 patients with localised/locally advanced disease and 5 with mRCC) will undergo multiparametric MRI including metabolic assessment by 13C-MRSI as well as measurements of tumour perfusion, oxygenation and water diffusion. Based on unique combinations of these imaging parameters, tumour habitats will be defined through cluster mapping. Intra-operative biopsy and immediate freezing of tumour cores from different imaging habitats will be performed. Multiregional ex-vivo tumour biopsies will be spatially registered to the MR data by using 3D printed tumour moulds generated from high-resolution volumetric sequences. These moulds allow slicing tumour specimens in identical planes to the MR acquisition and increase the accuracy of tissue sampling. The biopsies will then undergo molecular profiling such as metabolomic and RNA sequencing (including mRNA immune profiling), as well as DNA sequencing to assess driver and clonal mutations in different tumour habitats and guide design of appropriate cfDNA probes. The targeted cfDNA capture panel will be employed to determine the ability of cfDNA to capture the full genetic diversity observed in the primary tumour and to relate cfDNA to various imaging habitats in both renal and pancreatic cancer. Data integration will be undertaken between imaging habitats and molecular profiles for each region.

  Study EGAS00001003703

• Copy number profiling of putative cancer stem from pleural effusion aspirates from breast cancer patients

  A common symptom during late stage breast cancer disease is pleural effusion, which is related to poor prognosis. In pleural effusions malignant cells can be detected, which indicates a metastatic spread from the primary tumor site. Beside the usefulness for investigation of the metastatic process, it has been a source to isolate cells with putative cancer stem cell (CSC) properties. For this study, pleural effusion aspirate from 25 metastatic breast cancer patients were processed to establish in vitro culture. Patient-derived aspirates were cultured under sphere forming conditions and isolated primary cultures were sorted for cancer stem cell populations ALDH1+ and CD44+CD24-/low and sphere forming efficiency of CSC and non-CSC subpopulation was determined. Unsorted cells had a higher potential to form spheres compared to CSC subpopulations. Subsequently, DNA was isolated from pleural effusions before and after cell sorting, to compare their copy number profiles generated by low-coverage whole genome sequencing to DNA from FFPE primary tumors or bone metastasis. Unbalanced profiles could be detected in four pleural samples before sorting, in two ALDH1+ sorted samples, in one CD44+CD24-/low and in all control samples. Additionally, spheroid culture from one pleural sample could be passaged on and was characterized. Therefore, pleural effusions are an appropriate media for enrichment of putative CSC. However, low-coverage sequencing is difficult due to low quantities of abnormal cells compared to normal cells.

  Study EGAS00001002343

• ATRX mutant neuroblastoma is sensitive to EZH2 inhibition via modulation of neuronal differentiation.

  Mutations and structural alterations of the SWI/SNF-like chromatin remodeler ATRX have been reported at high frequency in a number of adult and pediatric tumors1. However, the consequences of ATRX (Alpha Thalassemia/Mental Retardation, X-linked) mutations in cancer and their underlying epigenetic sensitivities remain ill defined. Particularly intriguing are the large N-terminal deletions of ATRX in neuroblastoma that generate in-frame fusion (IFF) proteins1–3 devoid of key chromatin interaction domains. Here we demonstrate that neuroblastoma cells harbouring ATRX IFFs have distinct gene expression programs compared to neuroblastoma cells that are wild type for ATRX. This is due in part to H3K27me3-mediated silencing of REST (RE1 Silencing Transcription Factor) target genes involved in neuronal differentiation. In turn, we find that ATRX IFF cells display exquisite sensitivity to EZH2 inhibition in both adherent and tumorsphere conditions, due in part to derepression of neurogenesis genes, including REST targets. Examination of the epigenomic landscape of a pediatric neuroblastoma tumor harboring an ATRX IFF revealed that H3K27me3 occupies a subset of REST target genes that are transcriptionally silenced, and are sensitive to EZH2 inhibition in our cell-based assays. Thus, our study greatly advances our understanding of indolent neuroblastoma and identifies EZH2 inhibition as a potential therapy for patients with ATRX mutant disease. Further, these studies may be applicable to other ATRX mutant pediatric malignancies, particularly those carrying similar structural alterations of ATRX.

  Study EGAS00001002507

• Mutational landscape and patterns of clonal evolution in relapsed pediatric acute lymphoblastic leukemia

  Relapse of acute lymphoblastic leukemia (ALL) remains a leading cause of childhood death. Prior studies have shown clonal mutations at relapse often arise from relapse-fated subclones that exist at diagnosis. However, the genomic landscape, evolutionary trajectories and mutational mechanisms driving relapse are incompletely understood. In an analysis of 92 cases of relapsed childhood ALL, incorporating multimodal DNA and RNA sequencing, deep digital mutational tracking and xenografting to formally define clonal structure, we identify 50 significant targets of mutation with distinct patterns of mutational acquisition or enrichment. CREBBP, NOTCH1, and Ras signaling mutations rose from diagnosis subclones, whereas variants in NCOR2, USH2A and NT5C2 were exclusively observed at relapse. Evolutionary modeling and xenografting demonstrated that relapse-fated clones were minor (50%), major (27%) or multiclonal (18%) at diagnosis. Putative second leukemias, including those with lineage shift, were shown to most commonly represent relapse from an ancestral clone rather than a truly independent second primary leukemia. A subset of leukemias prone to repeated relapse exhibited hypermutation driven by at least three distinct mutational processes, resulting in heightened neoepitope burden and potential vulnerability to immunotherapy. Finally, relapse-driving sequence mutations were detected prior to relapse using deep digital PCR at levels comparable to orthogonal approaches to monitor levels of measurable residual disease. These results provide a genomic framework to anticipate and circumvent relapse by earlier detection and targeting of relapse-fated clones.

  Study EGAS00001003975

• The Southern African Human Genome Programme

  The Southern African Human Genome Programme (SAHGP) is a national and regional initiative that aims to unlock the unique genetic character of southern African populations. Its vision is to improve quality of life through understanding human genetic diversity. It is timely to advance the SAHGP initiative since global genomic initiatives are increasingly pointing to populations in Africa, more specifically the San and Khoe from southern Africa, as harbouring the most ancient genetic signatures among the living populations of the world. Studies on African populations have already generated a wealth of information upon which the SAHGP can build. This information will be used to promote and support genomic research programmes both nationally and regionally to address critical questions that would benefit the people of our region. The SAHGP aims to make a significant contribution to the understanding of DNA variation among southern Africans and how this impacts on the health of the people of this region. Aims of the SAHGP: To develop capacity for genomic research in southern Africa; To establish sustainable resources for genomic research in the region; and To translate knowledge and information into improvements in human health in the region. The objective is to derive the maximum benefit from the SAHGP, and that the data should be widely shared with the scientific community for the benefit of the people of the sub-continent and beyond.

  Study EGAS00001002639

• 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

• DNA hypermethylation and differential gene expression associated with Klinefelter syndrome

  The molecular basis for the phenotypic traits and morbidity in Klinefelter syndrome (KS) are not clarified. As DNA methylation affect gene expression and thereby play a role in disease susceptibility, we performed genome-wide DNA methylation profiling of leucocytes from peripheral blood samples from 67 KS patients, 67 male controls and 33 female controls, in addition to genome-wide RNA-sequencing profiling in a subset of 9 KS patients, 9 control males and 13 female controls. Characterization of the methylome as well as the transcriptome of both coding and non-coding genes identified a unique epigenetic and genetic landscape of both autosomal chromosomes as well as the X chromosome in KS. A subset of genes show significant correlation between methylation values and expression values. Gene set enrichment analysis of differentially methylated positions yielded terms associated with well-known comorbidities seen in KS. In addition, differentially expressed genes revealed enrichment for genes involved in the immune system, wnt-signaling pathway and neuron development. Based on our data we point towards many new candidate genes (AKAP17A, AMOT, APOB, DACT1, DDX58, DOCK7, EIF2S3, FIGNL1, G3BP1, HENMT1, HERC5, IFI44, IFI44L, IFIT1, IFIT3, ISG15, KANK1, LGALS1, NSD1, PEX10, PLSCR1, RSAD2, SHROOM2, SLC25A6, SPEG, SPON2, TXLNG), which may be implicated in the phenotype and further point towards non-coding genes (RP13-216E22.4, RP13-36G14.4, G087825, G088512), which may be involved in X chromosome inactivation in KS and in the regulation of escape genes.

  Study EGAS00001002797

• Atypical 3q26/MECOM rearrangements genocopy inv(3)/t(3;3) in acute myeloid leukemia

  Acute myeloid leukemia (AML) with inv(3)/t(3;3)(q21q26) is a distinct WHO recognized entity, characterized by its aggressive course and poor prognosis. In this subtype of AML, the translocation of a GATA2 enhancer (3q21) to MECOM (3q26) results in overexpression of the MECOM isoform EVI1 and monoallelic expression of GATA2 from the unaffected allele. The full-length MECOM transcript, MDS1-EVI1, is not expressed as the result of the 3q26 rearrangement. Besides the classical inv(3)/t(3;3), a number of other 3q26/MECOM rearrangements with poor treatment response have been reported in AML. Here we demonstrate, in a group of 33 AML patients with atypical 3q26 rearrangements, MECOM involvement with EVI1 overexpression, but no or low MDS1-EVI1 levels. Moreover, the 3q26 translocations in these AML patients often involve super-enhancers of genes active in myeloid development (e.g. CD164, PROM1, CDK6 or MYC). In more than 50% of these cases allele specific GATA2 expression was observed, either by copy number loss or by an unexplained allelic imbalance. Altogether, atypical 3q26 recapitulate the main leukemic mechanism of inv(3)/t(3;3) AML, namely EVI1 overexpression driven by enhancer hijacking, absent MDS1-EVI1 expression and potential GATA2 involvement. Therefore, we conclude that both atypical 3q26/MECOM and inv(3)/t(3;3) can be classified as a single entity of 3q26-rearranged AMLs. Routine analyses determining MECOM rearrangements, EVI1 and MDS1-EVI1 expression are required to recognize 3q-rearranged AML cases.

  Study EGAS00001004325

• WGS and WES data of invasive micropapillary carcinoma of breast

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

  Dataset EGAD00001008450

• 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

• Memory-like HCV-specific CD8+ T cells retain a chronic scar after cure of chronic HCV infection

  Chronic hepatitis C virus (HCV) infection is associated with CD8+ T-cell exhaustion characterized by limited effector functions and thus compromised anti-viral activity. Exhausted HCV-specific CD8+ T cells are comprised of memory-like and terminally exhausted CD8+ T-cell subsets. So far, little is not known about the molecular profile and fate of these cells after elimination of chronic antigen stimulation by direct acting antiviral therapy (DAA). Here, we report an antigen-driven molecular core signature underlying exhausted CD8+ T-cell subset heterogeneity in chronic viral infection with a progenitor/progeny relationship of memory-like and terminally exhausted HCV-specific CD8+ T cells via an intermediate stage. Furthermore, transcriptional profiling reveals that the memory-like cells remain after DAA-mediated cure while terminally exhausted HCV-specific CD8+ T-cell subsets are lost. Thus, the memory polarization of the overall HCV-specific CD8+ T-cell response after cure does not result from re-differentiation of exhausted T cells. Consequently, antigen elimination has little impact on the exhausted core signature of memory-like CD8+ T cells that remains clearly different from bona fide T-cell memory. These results identify a molecular signature of T-cell exhaustion that is imprinted like a chronic scar in HCV-specific CD8+ T cells even after HCV cure, highlighting the requirement of re-programming to elicit full effector potential of exhausted T cells.

  Study EGAS00001004538

• Identification of the mutational consequences of precancerous liver disease (including alcohol abuse) on the genomes of human adult stem cells.

  Stem cells are believed to be the cells-of-origin cancer. Nevertheless, their genome remains remarkably stable during life. To increase insight into how these cells can drive tumorigenesis, we set out to determine the genomic integrity of human adult stem cells from patients with liver disease associated with increased cancer risk. To this end, we characterized genome-wide mutation profiles of single ASCs from livers of patients with alcoholic liver disease, non-alcoholic liver disease and primary sclerosing cholangitis. The patients included in our study received a liver transplant due to their liver disease. Strikingly, we found that ASCs in these livers maintain a stable genome and have comparable mutation rates as healthy human liver ASCs. We did observe more nonsense and nonsynonymous mutations in COSMIC cancer genes that are hit by dominant mutations in tumors in comparison to healthy ASCs (n.s.), including two ASCs in alcoholic liver disease patients with a nonsense mutation in PTPRK. Next, we sequenced 5 biopsies of an alcoholic tumour to identify the mutational patterns in the tumour-initiating ASC. We found that the mutational patterns are, again, similar to healthy liver, but the mutational load is slightly higher. This suggests that the increased cancer risk in the liver disease patients is not caused by an increased mutational load, but rather that the precancerous liver environment promotes outgrowth of cells that would normally be selected against.

  Study EGAS00001002983

• Mutational Signature and Transcriptomic Classification Analyses as the Decisive Diagnostic Tools for a Cancer of Unknown Primary with Neuroendocrine Differentiation

  PURPOSE: Cancer of unknown primary (CUP) is a group of metastatic tumors in which the standard diagnostic work-up fails to identify the site of origin of the tumor. The potential impact of precision oncology on this group of patients is large since their tumors might have actionable driver mutations that can provide treatment options otherwise not available for patients with these fatal cancers. This study investigated if comprehensive genomics analyses could inform on the origin of the tumor. PATIENT AND METHODS: Here we describe a patient whose tumor was misdiagnosed at least three times. Next-generation sequencing, a PDX mouse model and bioinformatics was used to identify an actionable mutation, predict resistance development to the targeted therapy, and to correctly diagnose the origin of the tumor. The Cancer Genome Atlas was used to benchmark the bioinformatics workflow. RESULTS: Despite the lack of a known primary tumor site and the absence of diagnostic immunohistochemical markers, the origin of the patient's tumor was established using the novel bioinformatics workflow. This included a mutational signature analysis of the sequenced metastases and comparison of their transcriptomic profiles to a pan-cancer panel of tumors from The Cancer Genome Atlas. We further discuss the strengths and limitations of the latter approaches in the context of three potentially incorrectly diagnosed TCGA lung tumors. CONCLUSION: Comprehensive genomics analyses could inform on the origin of tumors in patients suffering from CUP.

  Study EGAS00001003026

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

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

  Study EGAS00001005235

• DNA methylation analysis on PBL obtained from male patients with UC-PSC, UC and HC

  Primary sclerosing cholangitis (PSC) is a chronic inflammatory liver disease affecting the intra- and extrahepatic bile ducts, and is strongly associated with ulcerative colitis (UC). In this study, we explored the peripheral blood DNA methylome and its immune cell composition in patients with PSC-UC, UC, and healthy controls (HC) with the aim to develop a predictive assay in distinguishing patients with PSC-UC from those with UC alone. The peripheral blood DNA methylome of male patients with PSC and concomitant UC, UC and HCs was profiled using the Illumina HumanMethylation Infinium EPIC BeadChip (850K) array. Differentially methylated CpG position (DMP) and region (DMR) analyses were performed alongside gradient boosting classification analyses to discern PSC-UC from UC patients. As observed differences in the DNA methylome could be the result of differences in cellular populations, we additionally employed mass cytometry to characterize the immune cell compositions. Genome wide methylation analysis revealed no differentially methylated positions between PSC-UC and UC patients. Nonetheless, using gradient boosting we were capable of discerning PSC-UC from UC with an area under the receiver operator curve (AUROC) of 0.80. Four CpG sites annotated to the NINJ2 gene were found to strongly contribute to the predictive performance. While mass cytometry analyses corroborated the largely similar blood cell composition among patients with PSC-UC, UC and HC, a higher abundance of myeloid cells was observed in UC compared to PSC-UC patients.

  Study EGAS00001005832

• Predictive value of chromosome 18q11.2-q12.1 loss for benefit from bevacizumab in metastatic colorectal cancer; a post-hoc analysis of the randomized controlled trial AGITG-MAX

  Patients and methods: We analyzed tumor samples of 256 patients enrolled in the AGITG-MAX 3-arm randomized controlled trial (RCT) who were randomly assigned to receive either capecitabine, or capecitabine and bevacizumab with- or without mitomycin. Chromosome 18q11.2-q12.1 copy number status was measured and correlated with PFS according to a predefined analysis plan with marker-treatment interaction as the primary end-point.Results: Chromosome 18q11.2-q12.1 losses were detected in 71% of patients (181/256) characteristic for mCRC. Consistent with the non-randomized study, significant PFS benefit of bevacizumab was observed in chromosome 18q11.2-q12.1 loss patients (P=0.009), and not in 18q no-loss patients (P=0.67). Although significance for marker-treatment interaction was not reached (Pinteraction = 0.28), hazard ratio and 95% confidence interval (HRinteraction 0.72, CI 0.39 - 1.32) show striking overlap with the non-randomized study cohorts (HRinteraction 0.41, CI 0.32 - 0.8) supported by a non-significant Cochrane chi-squared (P=0.11) for heterogeneity. Conclusions: PFS characteristics for marker positive- and negative patients randomized for treatment with- or without bevacizumab were consistent with those of the non-randomized discovery study, though the number of AGITG-MAX samples fell short to reach significance for marker-treatment interaction. Post-hoc analysis of the AGITG-MAX RCT provides supportive evidence for chromosome 18q11.2-q12.1 as a clinically relevant predictive marker for mCRC patients.

  Study EGAS00001005453

• A novel orthotopic patient-derived xenograft model of radiation-induced glioma following medulloblastoma

  Radiation-induced glioma (RIG) is a highly aggressive brain cancer arising as a consequence of radiation therapy. We report a case of RIG that arose in the brain stem following treatment for paediatric medulloblastoma, and the development and characterisation of a matched orthotopic patient-derived xenograft (PDX) model (TK-RIG915). Patient and PDX tumours were analysed using DNA methylation profiling, whole genome sequencing (WGS) and RNA sequencing. While initially thought to be a diffuse intrinsic pontine glioma (DIPG) based on disease location, results from methylation profiling and WGS were not consistent with this diagnosis. Furthermore, clustering analyses based on RNA expression suggested the tumours were distinct from primary DIPG. Additional gene expression analysis demonstrated concordance with a published RIG expression profile. Multiple genetic alterations that enhance PI3K/AKT and Ras/Raf/MEK/ERK signalling were discovered in TK-RIG915 including an activating mutation in PIK3CA, upregulation of PDGFRA and AKT2, inactivating mutations in NF1, and a gain-of-function mutation in PTPN11. Additionally, deletion of CDKN2A/B, increased IDH1 expression, and decreased ARID1A expression were observed. Detection of phospho-S6, -4EBP1 and -ERK via immunohistochemistry confirmed PI3K pathway and ERK activation. Here, we report one of the first PDX models for RIG, which recapitulates the patient disease and is molecularly distinct from primary brain stem glioma. Genetic interrogation of this model has enabled the identification of potential therapeutic vulnerabilities in this currently incurable disease.

  Study EGAS00001004709

• 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

• A Dose Escalation Study of Efmarodocokin Alfa (UTTR1147A) in Healthy Volunteers and Patients with Ulcerative Colitis

  Background: The interleukin-22 cytokine (IL-22) has demonstrated efficacy in nonclinical colitis models with a non-immunosuppressive mechanism-of-action. Efmarodocokin alfa (UTTR1147A) is a fusion protein agonist that links IL-22 to the crystallizable fragment (Fc) of human immunoglobulin (IgG4) for improved pharmacokinetic characteristics, but with a mutation to minimize Fc effector functions. Methods: This randomized, phase 1b study evaluated the safety, tolerability, pharmacokinetics, and pharmacodynamics of repeat intravenous dosing of efmarodocokin alfa in healthy volunteers (HVs; n=32) and patients with ulcerative colitis (UC; n=24) at 30–90 µg/kg doses given biweekly (Q2W) or monthly (Q4W) (6:2 active:placebo per cohort). Results: The most common adverse events (AEs) were on-target dermatological effects (dry skin, erythema, and pruritus) that were reversible. Dose-limiting non-serious dermatological AEs (severe dry skin, erythema, exfoliation, and discomfort) were seen in two HVs and one patient at 90 ug/kg Q2W. Pharmacokinetics were generally dose-proportional across the dose levels tested, but patients demonstrated lower drug exposures relative to HVs at the same dose. IL-22 serum biomarkers and IL-22 responsive genes in colon biopsies were induced with active treatment. Patients demonstrated changes in microbiota composition following active treatment, thereby reversing baseline dysbiosis. Clinical response was observed in 7/18 active- and 1/6 placebo-treated patients; clinical remission was observed in 5/18 active- and 0/6 placebo-treated patients.

  Study EGAS00001006172

• Single_cell_characterization_of_T_cell_lymphoma_

  Lymphomas are types of blood cancer derived from lymphocytes, and can be classified into two main categories, Hodgkin's lymphomas (HL) and non-Hodgkin lymphomas (NHL). T-cell lymphoma is a type of NHL that develops from T lymphocytes and accounts for 10-15% of total NHL cases. As a rare disease, T-cell lymphoma is functionally, pathologically, and clinically complex and heterogeneous. The 2016 revision of the WHO classification has listed 29 subtypes of mature T/NK-cell neoplasms, among which over 25 are T lymphocyte in origin. These subtypes of T-cell lymphomas have different biological behaviours and clinical prognosis. Clinically, T-cell lymphoma ranges from indolent (slow-growing) to aggressive (fast-growing and spreading) disease. According to the guideline booklet provided by the Lymphoma Research Foundation (https://lymphoma.org/aboutlymphoma/nhl/), indolent T-cell lymphoma subtypes include adult T-cell leukaemia/lymphoma (ATLL), cutaneous T-cell lymphoma (CTCL), and mycosis fungoides (MF), while the aggressive subtypes include anaplastic large cell lymphoma (ALCL), angioimmunoblastic T-cell lymphoma (AITL), peripheral T-cell lymphoma (PTCL), and T lymphoblastic lymphoma (TLBL). PTCL and ALCL represent the most common subtypes of aggressive T-cell lymphoma whereas TLBL is a rare and aggressive subtype mostly diagnosed in children. In this study, we will use single cell RNA sequencing and single cell multiomic methods three aggressive T-cell lymphoma subtypes (ALCL, PTCL and TLBL) and try to gain a deeper insight into the genotypic and phenotypic features of these T-cell lymphomas.

  Study EGAS00001005750

• Genomic stratification and liquid biopsy in a rare adrenocortical carcinoma (ACC) case, with dual lung metastases

  Adrenocortical carcinoma is a rare malignancy with a poor prognosis and few treatment options. Molecular characterization of this cancer remains limited. We present a case of an adrenocortical carcinoma (ACC) in a 37-yr-old female, with dual lung metastases identified 1 yr following commencement of adjuvant mitotane therapy. As standard therapeutic regimens are often unsuccessful in ACC, we undertook a comprehensive genomic study into this case to identify treatment options and monitor disease progress. We performed targeted and whole-genome sequencing of germline, primary tumor, and both metastatic tumors from this patient and monitored recurrence over 2 yr using liquid biopsy for ctDNA and steroid hormone measurements. Sequencing revealed the primary and metastatic tumors were hyperhaploid, with extensive loss of heterozygosity but few structural rearrangements. Loss-of-function mutations were identified in MSH2, TP53, RB1, and PTEN, resulting in tumors with mismatch repair signatures and microsatellite instability. At the cellular level, tumors were populated by mitochondria-rich oncocytes. Longitudinal ctDNA mutation and hormone profiles were unable to detect micrometastatic disease, consistent with clinical indicators of disease remission. The molecular signatures in our ACC case suggested immunotherapy in the event of disease progression; however, the patient remains free of cancer. The extensive molecular analysis presented here could be applied to other rare and/or poorly stratified cancers, to identify novel or repurpose existing therapeutic options, thereby broadly improving diagnoses, treatments, and prognoses.

  Study EGAS00001003497

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

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

  Study EGAS00001006118

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

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

  Study EGAS00001003586

• Prostate cancer ancestral genomic disparity

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

  Study EGAS00001006425

• Single-cell ATAC-seq analysis for COVID19 patients

  While SARS-CoV-2 infection causes mild respiratory disease in most individuals, a small group of patients develops severe COVID-19. Dysfunctional innate immune responses have been identified to contribute to differences in COVID-19 severity, but the key regulators are still unknown. Here, we present an integrative single-cell epigenetics, transcriptomic, and genetics analysis of peripheral blood mononuclear cells from hospitalized and convalescent COVID-19 patients. In classical monocytes, we identified 41.3% of significantly up-regulated genes in hospitalized COVID-19 patients potentially induced by differential chromatin accessibility. Sub-clustering and motif-enrichment analyses of monocytes reveal disease condition-specific regulation by transcription factors, such as C/EBPs and SPI1, and their targets, including a long-noncoding RNA LUCAT1, which further regulates interferon responses and is associated with the need for oxygen supply of COVID-19 patients. The interaction between C/EBPs and LUCAT1 was validated through loss-of-function experiments. Finally, we investigated genetic risk variants that exhibit allele-specific open chromatin (ASoC) in promoters/enhancers of COVID-19 patients. Integrating our data with publicly available expression quantitative trait loci and chromosomal interactions indicates that ASoC SNP rs6800484-C is associated with lower expression of CCR2, which may contributeto higher viral loads in lungs and higher risk of COVID-19 hospitalization. Altogether, our study highlights the diverse genetic and epigenetic regulators that contribute to the innate immune responses of different COVID-19 patients.

  Study EGAS00001006559

• Single-cell RNA-seq analysis for COVID19 patients

  While SARS-CoV-2 infection causes mild respiratory disease in most individuals, a small group of patients develops severe COVID-19. Dysfunctional innate immune responses have been identified to contribute to differences in COVID-19 severity, but the key regulators are still unknown. Here, we present an integrative single-cell epigenetics, transcriptomic, and genetics analysis of peripheral blood mononuclear cells from hospitalized and convalescent COVID-19 patients. In classical monocytes, we identified 41.3% of significantly up-regulated genes in hospitalized COVID-19 patients potentially induced by differential chromatin accessibility. Sub-clustering and motif-enrichment analyses of monocytes reveal disease condition-specific regulation by transcription factors, such as C/EBPs and SPI1, and their targets, including a long-noncoding RNA LUCAT1, which further regulates interferon responses and is associated with the need for oxygen supply of COVID-19 patients. The interaction between C/EBPs and LUCAT1 was validated through loss-of-function experiments. Finally, we investigated genetic risk variants that exhibit allele-specific open chromatin (ASoC) in promoters/enhancers of COVID-19 patients. Integrating our data with publicly available expression quantitative trait loci and chromosomal interactions indicates that ASoC SNP rs6800484-C is associated with lower expression of CCR2, which may contributeto higher viral loads in lungs and higher risk of COVID-19 hospitalization. Altogether, our study highlights the diverse genetic and epigenetic regulators that contribute to the innate immune responses of different COVID-19 patients.

  Study EGAS00001006560

• Direct Comparative Analysis of 10X Genomics Chromium and Smart-seq2

  Single cell RNA sequencing (scRNA-seq) is widely used for profiling transcriptomes of individual cells. The droplet-based 10X Genomics Chromium (10X) approach and the plate-based Smart-seq2 full-length method are two frequently-used scRNA-seq platforms, yet there are only a few thorough and systematic comparisons of their advantages and limitations. Here, by directly comparing the scRNA-seq data by the two platforms from the same samples of CD45- cells, we systematically evaluated their features using a wide spectrum of analysis. Smart-seq2 detected more genes in a cell, especially low abundance transcripts as well as alternatively spliced transcripts, but captured higher proportion of mitochondrial genes. The composite of Smart-seq2 data also resembled bulk RNA-seq data better. For 10X-based data, we observed higher noise for mRNA in the low expression level. Despite the poly(A) enrichment, approximately 10-30% of all detected transcripts by both platforms were from non-coding genes, with lncRNA accounting for a higher proportion in 10X. 10X-based data displayed more severe dropout problem, especially for genes with lower expression levels. However, 10X-data can better detect rare cell types given its ability to cover a large number of cells. In addition, each platform detected different sets of differentially expressed genes between cell clusters, indicating the complementary nature of these technologies. Our comprehensive benchmark analysis offers the basis for selecting the optimal scRNA-seq strategy based on the objectives of each study.

  Study EGAS00001003973

• Whole-genome sequencing of normal Singaporean volunteers

  Sleep is associated with various health outcomes. Despite their growing adoption, the potential for consumer wearables to contribute sleep metrics to sleep-related biomedical research remains largely uncharacterized. Here we analyzed sleep tracking data, along with questionnaire responses and multi-modal phenotypic data generated from 482 normal volunteers. First, we compared wearable-derived and self-reported sleep metrics, particularly total sleep time (TST) and sleep efficiency (SE). We then identified demographic, socioeconomic and lifestyle factors associated with wearable-derived TST. Among others, male gender (β = -15.539, 95% confidence interval [CI] = -26.245 - -4.832, p = 0.005), older age (β = -0.493, CI = -0.941 - -0.044, p = 0.032) and manual labor (β = -26.856, CI = -49.715 - -3.997, p = 0.022) were associated with reduced TST, whereas alcohol consumption (β = 19.247, CI = 8.008 - 30.486, p = 8.54x10-04) was associated with increased TST. Multi-modal phenotypic data analysis showed that wearable-derived TST and SE were associated with various cardiovascular disease risk markers such as body mass index (β = -0.006, CI = -0.00 - -0.000, p = 0.040) and waist circumference (β = 0.001, CI = -0.015-0.017, p = 0.893), whereas self-reported measures were not. We then showed that insufficient sleep was associated with premature telomere attrition, using telomere length estimated using whole-genome sequencing (β = 1.275, CI = 0.187 - 2.363, p = 0.023) and quantitative PCR (β = 0.001, CI = 8.318x10-05 - 0.001, p = 0.028). Our study highlights the potential for sleep metrics from consumer wearables to provide novel insights into data generated from population cohort studies.

  Study EGAS00001004007

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

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

  Study EGAS00001006582