23144 results for "ega"

in 61.24 milliseconds.

• Tubulointerstitial fibrosis is the histological hallmark of chronic kidney disease (CKD). Hypoxia and inflammation (i.e., interleukin (IL)-1β signalling) are independent mediators of tubulointerstitial fibrosis. However, the physiological response of human kidney tubular cells to IL-1β/IL-1RI signalling under the hypoxic conditions of CKD is poorly understood and remains a clinical imperative for therapeutic targeting. This study reports that hypoxia and IL-1β act in synergy to trigger cell cycle arrest/cellular senescence of ex vivo patient-derived primary proximal tubular epithelial cells (PTECs).

  Hypoxia and interleukin (IL)-1β are independent mediators of tubulointerstitial fibrosis, the histological hallmark of chronic kidney disease (CKD). Here, we examine how hypoxia and IL-1β act in synergy to augment maladaptive proximal tubular epithelial cell (PTEC) repair in human CKD. Ex vivo patient-derived PTECs were cultured under normoxic (21% O2) or hypoxic (1% O2) conditions in the absence or presence of IL-1β and examined for maladaptive repair signatures. Hypoxic PTECs incubated with IL-1β displayed a discrete transcriptomic profile distinct from PTECs cultured under hypoxia alone, IL-1β alone or under normoxia. Hypoxia+IL-1β-treated PTECs had 692 ‘unique’ differentially expressed genes (DEGs) compared to normoxic PTECs, with ‘cell cycle’ the most significantly enriched KEGG pathway based on ‘unique’ down-regulated DEGs (including CCNA2, CCNB1 and CCNB2). Hypoxia+IL-1β-treated PTECs displayed signatures of cellular senescence, with reduced proliferation, G2/M cell cycle arrest, elevated senescence-associated β-galactosidase (SA-β-gal) activity and increased production of pro-inflammatory/fibrotic senescence-associated secretory phenotype (SASP) factors compared to normoxic conditions. Treatment of Hypoxia+IL-1β-treated PTECs with either a senolytic drug combination, quercetin+dasatinib, or a type I IL-1 receptor (IL-1RI) neutralizing antibody attenuated senescent cell burden. Our data identify a mechanism whereby hypoxia in combination with IL-1β/IL-1RI signalling trigger PTEC senescence, providing novel therapeutic and diagnostic check-points for restoring tubular regeneration in human CKD.

  Study EGAS00001007904

• Identification of intronic sequences promoting natural skipping of ABCA4 exon 15 using long-read transcriptomics and midigene assays

  Stargardt disease is an inherited retinal disease that leads to progressive loss of vision, and is caused by biallelic variants in the ABCA4 gene. A significant proportion of ABCA4 variants of variable penetrance affect pre-mRNA splicing. Previously, we observed the exclusion of certain ABCA4 exons into the final transcript in the absence of any disease-causing variants. In this study, we explored this phenomenon that we coined natural exon skipping (NES) for the entire ABCA4 transcript, using targeted RT-PCR analysis as well as long-read transcriptome sequencing of control retina samples. For the most abundant NES event, which is skipping of ABCA4 exon 15, overlapping ABCA4 deletion constructs were employed to identify responsible cis-regulatory sequence motifs. Subsequently, antisense oligonucleotides (AONs) were designed and tested to reverse NES and restore fully-spliced ABCA4 mRNA. We identified intronic regions on each side of exon 15 that, when deleted, increased the inclusion of exon 15 in the final transcripts. Interestingly, these motifs appear to work in a synergistic manner. Further fine-mapping of these regions revealed several cis-regulatory motifs that may be responsible for inducing NES, on each side of the exon. AONs targeting these motifs unfortunately did not decrease the skipping of exon 15. In conclusion, we discovered an intriguing phenomenon that, in the absence of disease-causing variants, results in the skipping of one or more exons in the ABCA4 transcript. These NES events appear to be regulated by cis-regulatory motifs within the flanking introns. This work captured the extend of NES throughout ABCA4 transcript and the role of cis-regulatory motifs in pre-mRNA splicing.. The use of AONs targeting cis-regulatory motifs is an innovative strategy for splicing modulation, suitable for complementing future variant-specific splicing modulation therapies in ABCA4 retinopathies.

  Study EGAS50000000071

• Blueprint Data Access Committee.

  Access to the sequence and alignment and genotype level data produced by the Blueprint Consortium is controlled by the Blueprint Data Access Committee (DAC). Access to data will be granted to qualified investigators for appropriate use. Please email blueprint-dac@ebi.ac.uk to request the Data Access Application and Data Access Agreement.

  Dac EGAC00001000135

• Integrating molecular imaging and transcriptomic profiling in advanced HER2-positive breast cancer receiving trastuzumab emtansine (T-DM1): an analysis of the ZEPHIR clinical trial

  Purpose: The ZEPHIR clinical trial evaluated the role of [89Zr]trastuzumab-PET/CT (HER2-PET/CT) and 2-[18F]fluoro-2-deoxy-D-glucose PET/CT ([18F]FDG-PET/CT) in predicting outcomes in patients with advanced HER2-positive breast cancer treated with trastuzumab emtansine (T-DM1). Here, we combined molecular/metabolic imaging and transcriptomic data to investigate the biological processes associated with [89Zr]trastuzumab and [18F]FDG uptake, and to dissect the mechanisms involved in T-DM1 resistance. Experimental design: RNA was extracted from metastasis biopsies obtained at study inclusion in the ZEPHIR trial. HER2-PET/CT and [18F]FDG-PET/CT imaging data of biopsied lesions were integrated with RNA sequencing data. Lesions were compared based on the level of [89Zr]trastuzumab uptake (dichotomized into high/low classes) as well as on the presence/absence of metabolic response, defined comparing baseline and on-treatment [18F]FDG-PET/CT. Results: We analyzed matched transcriptomic and molecular/metabolic imaging data for 24 metastases (one lesion/patient). Pathways involved in extracellular matrix (ECM) organization, glycosylphosphatidylinositol synthesis, and collagen formation were enriched in lesions presenting low [89Zr]trastuzumab uptake. [18F]FDG uptake at baseline correlated with proliferation and immune-related processes. Hypoxia and ECM-related processes were enriched in lesions that showed no metabolic response to T-DM1. Interestingly, immune-related processes (e.g., Fc gamma receptor activation, complement cascade) were associated with high [89Zr]trastuzumab uptake and metabolic response. Conclusions: To our knowledge, this study represents the first correlative analysis between [89Zr]trastuzumab tumor uptake and gene expression profiling in humans. Our findings point toward a role of ECM in impairing [89Zr]trastuzumab tumor uptake and T-DM1 metabolic response in advanced HER2-positive breast cancer, and highlight the potential of molecular imaging to depict tumor microenvironment features.

  Study EGAS50000000470

• Genomic Profiling of Pediatric Tumors by Cell Free DNA Sequencing

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

  Study EGAS50000000072

• Integrative analysis reveals a macrophage-predominant, immunosuppressive immune microenvironment and subtype-specific therapeutic vulnerabilities in advanced salivary gland cancer

  Background Salivary gland cancers (SGC) are rare and heterogeneous malignant tumors. Advanced SGC lack established treatment options and show poor response to immunotherapy. Here, an integrative multi-omics analysis in a large cohort of advanced SGC revealed insights into the tumor immune microenvironment (TIM) and distinct mechanisms of immune evasion. Methods A total of 104 patients with recurrent/metastatic SGC from the DKTK MASTER program were included in this study. Whole-exome or whole-genome sequencing and RNA-sequencing was performed on fresh frozen tumor tissue. The tumor immune microenvironment was analyzed using CIBERSORT deconvolution analysis and immune gene expression scores in bulk RNA-sequencing data. Single-nuclei sequencing and immunohistochemistry analyses were performed in selected samples. Results were validated in bulk RNA-sequencing data of a previously published independent dataset. Results Bulk transcriptome analysis revealed an immune-deserted TIM in the majority of advanced SGC samples. Immune exclusion was most prominent in adenoid cystic carcinoma (ACC) subgroup 1 exhibiting a downregulation of the antigen processing machinery. Only a small subset of advanced SGC, including few adenoid cystic carcinoma, exhibited T-cell inflammation, which was correlated with tumor mutational burden in Non-ACC samples. Subtype specific expression of immune checkpoints as well as cancer testis antigens were identified with prominent expression of VTCN1 in luminal cells within ACC. Single-cell RNA-sequencing and bulk RNA-seq deconvolution analysis validated immune cell exclusion and revealed a TIM that was dominated by M2 macrophages across SGC subtypes. Among evaluable patients treated with immune checkpoint inhibitors, a high T-cell to macrophage ratio was associated with clinical benefit. Conclusions These data support biomarker-based development of immune-checkpoint inhibition and the development of novel immune-checkpoint inhibitors and cellular therapies in SGC.

  Study EGAS50000000809

• A prospective study identifies MisMatch Repair genes as candidate predisposing genes for uveal melanoma.

  Objective: The aim of this study is to identify new predisposing genes for uveal melanoma (UM), the most common intraocular malignancy in adults. Design: Prospective cohort study from Institut Curie. Participants: 381 unselected patients diagnosed with UM between July 2021 and February 2023 at Institut Curie. Methods: All patients received genetic counseling and consented to extended genetic testing. A panel of 122 genes was analyzed by targeted sequencing on germline DNA from these patients. Main outcomes measures: UM was diagnosed by ophthalmologists at Institut Curie, and a senior pathologist confirmed the diagnosis when tumor or biopsy was available. Classification of the germline genetic variants was performed according to the ACMG guidelines. Results: We identified at least one pathogenic variant (PV) in 79 patients. Among those variants, 20 were found in clinically relevant genes for the patient’s follow-up and genetic counseling. We found a significant enrichment of germline pathogenic variants in Mismatch Repair (MMR) MLH1, MSH2 and MSH6 genes in this UM cohort, suggesting that MMR germline alterations could predispose to UM. From this cohort, one tumor was available from a patient carrying a MLH1 germline mutation. The tumor was characterized by a monosomy 3, in line with a loss of heterozygosity (LOH) of MLH1, located on chromosome 3. Loss of expression of MLH1 and PMS2 was observed by immunohistochemistry, and MMR mutational signatures SBS6, ID1 and ID2 were identified from the whole genome sequencing of the tumor. These data suggest that MLH1 is implicated in the oncogenesis of this UM. Conclusions: This study is the largest prospective constitutional study performed on UM patients, and the first demonstrating that MMR germline alterations are significantly enriched among these patients. Further analyses will be needed to estimate the risk of developing UM for patients with Lynch syndrome

  Study EGAS50000000914

• Emirati Genome Project subset of 43,608 WGS samples for population-scale variant discovery and allele frequency mapping.

  Here, we present a comprehensive genomic characterization of a cohort of 43,608 Emirati genomes sequenced as part of the Emirati Genome Program (EGP). This study identified more than 421 million single nucleotide variants and indels and more than 600 million copy number- and structural variants. Small variants had 756 million molecular effects annotated. Of the identified variants, 69.04% were observed with more than one allele. Of 7.7 million polymorphic variants having an allele frequency (AF) of more than 5% in EGP, 1,348 have a predicted deleterious effect on a protein. Characterization of global variation reveals that EGP represents a genetic continuum encompassing the range of African, Asian, and European populations. It is best described by two Arabian, an Eurasian, and an African component, with the predominant Arabian component linked to mitochondrial haplogroups J and T, which are commonly attributed to the Middle East. Various aneuploidies of sex chromosomes were detected in 93 individuals overall and of chromosome 21 in 41 individuals. The median inbreeding coefficient and cumulative runs of homozygosity (ROHs) lengths were increased due to extensive consanguinity, being largest in the groups with Arabian main ancestry components and higher than reported for Qatar. Families were identified based on genetic relatedness and classified into 264 families with unrelated parents and 247 families with third- and fourth-degree consanguineous parents. Representative consanguineous pedigrees of families in EGP were outlined. Cumulative ROHs were affected by the main ancestry component and significantly increased in offspring of consanguineous parents, with a pronounced difference between 3rd and 4th-degree relatedness. Investigation of cumulative AFs of variants causing Mendelian diseases highlighted genes related to alpha- and beta-thalassemia (HBB, HBA2). It showed a high burden of variants causing severe recessive diseases, metabolic and retinal disorders, and hearing loss. In summary, EGP represents a landmark effort in characterizing the genetic diversity of the Emirati population, leveraging the largest Middle Eastern cohort reported to date.

  Study EGAS50000001071

• Single-Gene vs. Panel Sequencing in Advanced HR+/HER2− Breast Cancer

  Background: Hormone receptor-positive (HR+)/HER2-negative (HER2−) breast cancer is the most common subtype of breast cancer, with biomarker-driven therapies improving outcomes. Liquid biopsies, which analyze circulating tumor DNA (ctDNA), are increasingly essential in clinical management, offering a minimally invasive approach to assess somatic alterations and guide therapy decisions. However, the choice of assay directly affects clinical utility and global access to validated testing platforms remains inconsistent. To address these challenges, this study aimed to compare the performance of single-gene and panel-based sequencing for detecting PIK3CA mutations in HR+/HER2− breast cancer and evaluate the added value of broader genomic profiling in guiding treatment decisions. Patients and methods: We conducted a prospective, multicenter study comparing our previously established SiMSen-Seq (SSS) assay, a single-gene method targeting PIK3CA hotspot mutations, with the AVENIO ctDNA Expanded assay, enriching for 77 clinically relevant cancer genes. Additionally, tumor fractions were estimated with mFAST-SeqS by assessing chromosomal aneuploidy. We analyzed 161 plasma samples collected from 146 patients with HR+/HER2− advanced breast cancer before starting a new line of palliative therapy. Results: High concordance (92.6%) was observed between SSS and AVENIO for PIK3CA mutation detection, with accurate variant allele frequency (VAF) quantification confirmed by Bland-Altman analysis and linear mixed-effects modeling. AVENIO identified additional PIK3CA mutations in 19% of samples as well as other actionable alterations, including ESR1 mutations (17.5%) and PI3K pathway alterations (40.6%). Panel sequencing improved the interpretation of negative liquid biopsy results by integrating tumor fraction and VAF metrics. Conclusions: Panel sequencing assays such as AVENIO provide broader genomic coverage, enhance the detection of actionable mutations and refine the interpretation of negative liquid biopsy findings. These findings support the integration of comprehensive ctDNA profiling into routine clinical practice for HR+/HER2− advanced breast cancer while highlighting the need for improved assay accessibility and expanded genomic targets.

  Study EGAS00001008200

• CDK4/6 inhibition in advanced chordoma: final results of the NCT PMO-1601

  Background. Chordoma is a rare bone tumor with no approved systemic therapy. Inactivation of the tumor suppressor p16 (encoded by CDKN2A) leads to aberrant cyclin D-CDK4/6-RB pathway activity, which can be inhibited by palbociclib in chordoma cell lines. Patients and Methods. We conducted a phase II single-arm, open-labeled trial on palbociclib in adult patients with advanced chordomas with p16 loss (by immunohistochemistry) or CDKN2A loss (by genomic analysis) and CDK4/6 and RB1 presence (by immunohistochemistry or RNA-Sequencing). Patients received 125mg palbociclib once daily for 21 days in a 28-day cycle. The study used a Simon optimal 2-stage design. The primary endpoint was disease control rate after six cycles. Response was evaluated according to Response Evaluation Criteria in Solid Tumors (RECIST) version 1.1. Secondary endpoints included progression-free survival, overall survival, and biomarker analysis. The study was considered positive if 25% of patients reached the primary endpoint. Results. Between December 2017 and April 2023, 28 patients (median age 60 years, range 31-84) were enrolled. After a median follow-up of 28 months, the disease control rate was 39%, with 11 out of 28 patients achieving stable disease. No objective responses were observed. The median progression-free survival was 5.6 months, and the median overall survival was 24.6 months. Treatment was well tolerated without new safety signals. CDKN2A deletions were the most common genomic alterations underlying p16 loss and were accompanied by MTAP co-deletion. There was no correlation between responder phenotypes by immunohistochemistry and outcome. RNA-based analysis showed heterogeneous pathways deregulation at study entry suggesting underlying primary resistance mechanisms Conclusion. This trial met its primary endpoint with a 39% disease control rate, though antitumor activity was modest. Molecularly tailored combination therapies should be considered to improve efficacy. (ClinicalTrials.gov number, NCT03110744).

  Study EGAS00001007985

• A circulating biomarker for severity of facioscapulohumeral muscular dystrophy

  Facioscapulohumeral muscular dystrophy (FSHD) is a prevalent, incurable skeletal myopathy. Clinical trials for FSHD are hindered by heterogeneous biomarkers poorly associated with clinical severity, requiring invasive muscle biopsy. Macroscopically FSHD presents with slow fatty replacement of muscle, rapidly accelerated by inflammation. Mis-expression of the transcription factor DUX4 is currently accepted to underlie FSHD pathogenesis and mechanisms including PAX7 target gene repression have been proposed. Here we perform RNA-sequencing on MRI guided inflamed (n=24) and isogenic non-inflamed (n=24) muscle biopsies from the same clinically characterised FSHD patients, alongside isogenic peripheral blood mononucleated cells from a subset of patients (n=13), and unaffected controls. We employ multivariate models to evaluate the clinical associations of 5 published FSHD transcriptomic biomarkers. We demonstrate that PAX7 target gene repression can discriminate control, inflamed and non-inflamed FSHD muscle independently of age and sex (p<0.013), while the discriminatory power of DUX4 target genes is limited to distinguishing FSHD muscle from control. Importantly, the level of PAX7 target gene repression in non-inflamed muscle associates with clinical assessments of FSHD severity (p=0.04). DUX4 target gene biomarkers in FSHD muscle show associations with lower limb fat fraction and D4Z4 array length, but not clinical assessment. Lastly, PAX7 target gene repression in FSHD muscle correlates with the level in isogenic peripheral blood mononucleated cells (p=0.002). A refined PAX7 target gene biomarker comprising 143/601 PAX7 target genes computed in peripheral blood (the FSHD muscle-blood biomarker) associates with clinical severity in FSHD patients (p<0.036). Our new circulating biomarker validates as a classifier of clinical severity in an independent data set of 54 FSHD patient blood samples, with improved power in older patients (p=0.03). In summary, we present the minimally invasive FSHD muscle-blood biomarker of FSHD clinical severity valid in patient muscle and blood, with potential utility in routine disease monitoring and clinical trials.

  Study EGAS00001007350

• First genome-wide association study in an Australian Aboriginal population provides insights into genetic risk factors for Body Mass Index and Type 2 Diabetes

  A body mass index (BMI) >22kg/m2 is a risk factor for type 2 diabetes (T2D) in Aboriginal Australians. To identify loci associated with BMI and T2D we undertook a genome-wide association study using 1,355,133 quality-controlled single nucleotide polymorphisms (SNPs) genotyped (Illumina 2.5M Duo Beadchip) in 402 individuals in extended pedigrees from a Western Australian Aboriginal community. Imputation using the thousand genomes (1000G) reference panel extended the analysis to 6,724,284 post quality-control autosomal SNPs. No associations achieved genome-wide significance, commonly accepted as P<5x10-8. Nevertheless, genes/pathways in common with other ethnicities were identified despite the arrival of Aboriginal people in Australia >45,000 years ago. The top hit (rs10868204 Pgenotyped=1.50x10-6; rs11140653 Pimputed_1000G=2.90x10-7) for BMI lies 5’ of NTRK2, the type 2 neurotrophic tyrosine kinase receptor for brain-derived neurotrophic factor (BDNF) that regulates energy balance downstream of melanocortin-4 receptor (MC4R). PIK3C2G (rs12816270 Pgenotyped=8.06x10-6; rs10841048 Pimputed_1000G=6.28x10-7) was associated with BMI, but not with T2D as reported elsewhere. BMI also associated with CNTNAP2 (rs6960319 Pgenotyped=4.65x10-5; rs13225016 Pimputed_1000G=6.57x10-5), previously identified as the strongest gene-by-environment interaction for BMI in African-Americans. The top hit (rs11240074 Pgenotyped=5.59x10-6, Pimputed_1000G=5.73x10-6) for T2D lies 5’ of BCL9 that, along with TCFL2, promotes beta-catenin’s transcriptional activity in the WNT signaling pathway. Additional hits occurred in genes affecting pancreatic (KCNJ6, KCNA1) and/or GABA (GABRR1, KCNA1) functions. Notable associations observed for genes previously identified at genome-wide significance in other populations included MC4R (Pgenotyped=4.49x10-4) for BMI and IGF2BP2 Pimputed_1000G=2.55x10-6) for T2D. Our results may provide novel functional leads in understanding disease pathogenesis in this Australian Aboriginal population.

  Study EGAS00001001004

• Nala GSI GSAv3 PGx Study

  Background: Microarrays are a well-established and widely adopted technology capable of interrogating hundreds of thousands of loci across the human genome. Combined with imputation to cover common variants not included in the chip design, they offer a cost-effective solution for large-scale genetic studies. Beyond research applications, this technology can be applied for testing pharmacogenomics, nutrigenetics, and complex disease risk prediction. However, establishing clinical reporting workflows requires a thorough evaluation of the assay's performance, which is achieved through validation studies. In this study, we performed pre-clinical validation of a genetic testing workflow based on the Illumina Global Screening Array for 25 pharmacogenomic-related genes. Methods: To evaluate the accuracy of our workflow, we conducted multiple pre-clinical validation studies. Here, we present the results of accuracy and precision assessments, involving a total of 73 cell lines. These assessments encompass reference materials from the Genome-In-A-Bottle (GIAB), the Genetic Testing Reference Material Coordination Program (GeT-RM) projects, as well as additional samples from the 1000 Genomes project (1KGP). We conducted an accuracy assessment of genotype calls for target loci in each indication against established truth sets. Results: In our per-sample analysis, we observed a mean analytical sensitivity of 99.39% and specificity 99.98%. We further assessed the accuracy of star-allele calls by relying on established diplotypes in the GeT-RM catalogue or calls made based on 1KGP genotyping. On average, we detected a diplotype concordance rate of 96.47% across 14 pharmacogenomic-related genes with star allele-calls. Lastly, we evaluated the reproducibility of our findings across replicates and observed 99.48% diplotype and 100 % phenotype inter-run concordance. Conclusion: Our comprehensive validation study demonstrates the robustness and reliability of the developed workflow, supporting its readiness for further development for applied testing.

  Study EGAS00001007710

• mRNA capture sequencing and RT-qPCR for the detection of pathognomonic, novel and secondary fusion transcripts in formalin-fixed paraffin-embedded tissue: a sarcoma showcase

  RNA sequencing is widely used to study gene abundance patterns, but less frequently to assess structural alterations such as fusion transcripts. In this study, we assess the performance of messenger RNA (mRNA) capture sequencing to identify fusion transcripts in formalin-fixed paraffin-embedded (FFPE) tissue of different types of sarcoma, followed by orthogonal confirmation using RT-qPCR. To validate our workflow, we used the TruSight RNA Pan-Cancer Panel to analyze a first cohort of 6 positive control tumors for which the diagnostic workup demonstrated the presence of a specific chromosomal rearrangement. Fusion transcript calling by FusionCatcher confirmed all known aberrations and enabled the identification of both fusion gene partners and breakpoints. Next, whole transcriptome TruSeq RNA Exome sequencing was applied to a second cohort of 17 patients, diagnosed with fusion gene-negative alveolar rhabdomyosarcoma (ARMS) or undifferentiated round cell sarcoma (URCS), for whom fluorescence in situ hybridization (FISH) did not identify the classical pathognomonic rearrangements. Remarkably, for 6 out of 17 patients, mRNA capture sequencing readily detected a pathognomonic fusion transcript, i.e. PAX3-FOXO1 in 2 ARMS patients, and EWSR1-FLI1, EWSR1-ERG or EWSR1-NFATC2 in 4 URCS patients. These data show that RNA capture sequencing may enhance the detection rate of pathognomonic fusion genes in sarcoma. For the 11 patients with no known pathognomonic fusion transcript, 11 newly identified fusion transcripts were confirmed by RT-qPCR, including COPS3-TOM1L2, NCOA1-DTNB, WWTR1-LINC01986, PLAA-MOB3B, AP1B1-CHEK2 and BRD4-LEUTX fusion transcripts in ARMS patients. Additionally, recurrently detected secondary fusion transcripts in two patients diagnosed with EWSR1-NFATC2-positive sarcoma were confirmed (COPS4-TBC1D9, PICALM-SYTL2, SMG6-VPS53 and UBE2F-ALS2). Our results provide new insights into the underlying genetic causes of these malignancies.

  Study EGAS00001005202

• In this study, blood-brain barrier (BBB)-forming brain endothelial-like cells were generated from apolipoprotein E gene allele E4 (APOE4, high AD risk) and allele E3 (APOE3, lower AD risk) carrying patient-derived induced pluripotent stem cells (iPSCs). Cells were subsequently exposed to focused ultrasound and microbubbles (FUS+MB) to induce BBB opening and their transcriptome analysed. RNA sequencing (RNA-seq) results demonstrated minimal changes in the gene expression following FUS+MB suggesting safety of FUS+MB application in the clinical setting.

  The blood-brain barrier (BBB) is a major impediment to therapeutic intracranial drug delivery for the treatment of neurodegenerative diseases, including Alzheimer’s disease (AD). Focused ultrasound applied together with microbubbles (FUS+MB) is a novel technique to transiently open the BBB and increase drug delivery. Evidence suggests that FUS+MB is safe, however the effects of FUS+MB on human BBB cells, especially in relation to AD, remain sparsely investigated. Here we developed a human sporadic AD BBB cell platform to investigate the effects of FUS+MB on BBB cells and screen for the delivery of two potentially therapeutic AD antibodies. We generated BBB cells from apolipoprotein E gene allele E4 (APOE4, high AD risk) and allele E3 (APOE3, lower AD risk) carrying patient-derived induced pluripotent stem cells (iPSCs). We utilized this robust and reproducible platform to demonstrate increased delivery of potentially therapeutic antibodies across the BBB, including an analogue of Aducanumab (AduhelmTM) following FUS+MB treatment. Our results also demonstrate the safety of FUS+MB indicated by minimal changes in the cell transcriptome as well as little or no changes in cell viability and inflammatory responses within first 24 h post FUS+MB, with potentially important implications in a clinical setting. Finally, we report a more physiologically relevant hydrogel-based 2.5D BBB model as a key development for FUS+MB-mediated drug delivery screening, with potentially higher translational utility. Our results provide an important advancement for identifying FUS+MB-deliverable drugs and screening for cell- and patient-specific effects of FUS+MB, accelerating its use as a therapy in AD.

  Study EGAS00001005944

• Genomic profiling of paediatric high grade gliomas from the HERBY clinical trial

  The HERBY trial was a phase-II open-label, randomised, multicentre trial evaluating bevacizumab in patients with newly-diagnosed non-brainstem high grade glioma (HGG) between the ages of 3-18yrs. 121 patients were randomised with 1yr event-free survival as the primary end-point. Confirmation of HGG diagnosis by reference pathology was mandatory before randomisation, followed by review with five independent expert neuropathologists. We collected specimens from 89 patients consenting to translational research, and performed Sanger sequencing for H3F3A, Illumina 450K BeadArray methylation profiling, and whole exome sequencing and RNA sequencing where possible. 7/89 patients (8%) harboured H3F3A G34R/V mutations, whilst 24/89 (27%) harboured H3F3A K27M, the latter reflecting a high proportion of the novel entity recognised in the 2016 WHO classification of diffuse midline glioma with H3K27M mutation. Both histone mutations conferred a significantly shorter progression-free (p=0.0104) and overall survival (p=0.00159). 450K methylation subtyping additionally identified a number of subgroups in histone wild-type cases. These included infrequent (4/86, 4.6%) cases with IDH1 mutation in older children, further categorised into astrocytic ATRX-mutant (n=3) or 1p19q co-deleted with TERT promoter mutation (n=1). 9/74, (12%) of cases were classified as biologically resembling pleomorphic xanthoastrocytoma (PXA) by methylation profiling, with 5/9 harbouring BRAF V600E mutations and epithelioid histology, and 3/9 NF1 mutation and giant cell features. Three cases had methylation profiles more closely resembling low grade gliomas, though were morphologically high grade, and harboured MAPK dysregulation in two cases, with the third part of an additional infant (<3yrs) cohort. Finally, four cases had a mutational burden several orders of magnitude higher than the rest, with (2197-5332) somatic coding variants/sample reflecting a hypermutator phenotype. These data provide an important insight into the wide biological diversity of ‘HGG’ found in the paediatric age group and allow for a profound refinement of clinical trial interpretation.

  Study EGAS00001002328

• High-throughput telomere length measurement at nucleotide resolution using the PacBio high fidelity sequencing platform

  Telomeres are specialized nucleoprotein structures at the ends of linear chromosomes that are essential for the stability as well as the complete replication of the human genome. The progressive shortening of steady-state telomere length in normal human somatic cells is a promising biomarker for age-associated diseases. However, there remain substantial challenges in quantifying telomere length in clinics and research laboratories due to the lack of accurate and high-throughput telomere length measurement method. Here, we describe a novel workflow to capture and analyze telomeres accurately at the single-nucleotide resolution. By using newly designed telobaits that are complementary to the single-stranded G-rich 3’ telomeric overhangs, we can specifically and efficiently capture full-length telomeres and subject them to single-molecule real-time (SMRT) sequencing. We show that the telomere length can be measured at nucleotide resolution using the PacBio high fidelity (HiFi) sequencing platform, which provides the necessary accuracy and throughput to uncover the steady-state telomere length distribution in human culture cell lines as well as clinical patient samples. This new method is a valuable tool not only in population-wide epidemiology studies but also for the longitudinal assessment of telomere attrition in individuals. Our results also reveal the extreme heterogeneity of telomeric variant sequences (TVSs) that is dispersed throughout the telomere repeat region. Importantly, the unique distribution pattern and sequence of TVSs in individuals may be used as a fingerprint for sample identification. In addition, the presence of TVSs disrupts the continuity of the canonical (5’-TTAGGG-3’)n telomere repeats, which could affect the binding of shelterin complexes at the chromosomal ends. Therefore, TVSs may act as a form of genetic polymorphisms that impairs telomere protection. Together with the absolute telomere length in each individual, TVSs may have profound implications in human aging and diseases.

  Study EGAS00001006103

• Local In Time Statistics for processual research

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

  Study EGAS00001002520

• RNAseq data to study PRPF6 regulated splice forms in colon cancer cell lines

  Alternative splicing plays critical roles in differentiation, development, and cancer (Pettigrew et al., 2008; Chen and Manley, 2009). The recent identification of specific spliceosome inhibitors has generated interest in the therapeutic potential of targeting this cellular process (van Alphen et al., 2009). Using an integrated genomic approach, we have identified PRPF6, an RNA binding component of the pre-mRNA spliceosome, as an essential driver of oncogenesis in colon cancer. Importantly, PRPF6 is both amplified and overexpressed in colon cancer, and only colon cancer cells with high PRPF6 levels are sensitive to its loss. Our data clearly point to an important role for PRPF6 in colon cancer growth and suggest that a better understanding of its role in alternative splicing in colon cancer is warranted. To determine the specific alternative splice forms that PRPF6 regulates in colon cancer, we plan three experiments: 1. The first involves knocking down expression of PRPF6 in two different cancer cell lines with 3 different siRNAs, and then completing RNA-seq to determine the gene expression changes that occur relative to a non-targeting control siRNA. Because of the role for PRPF6 in pre-mRNA splicing, we especially want to quantify the changes in splice-specific forms of all genes genome-wide to identify genes whose splicing is altered upon PRPF6 knockdown. 2. The second involves immunoprecipitating PRPF6 from two different cancer cell lines and isolating any RNA that is bound to PRPF6, since PRPF6 is an RNA-binding protein. We then want to carry out RNA-seq to identify which RNA molecules co-immunoprecipitated with PRPF6. This will help us determine possible functions for PRPF6 in regulating colon cancer growth. 3. The third involves overexpressing PRPF6 in cell lines and then carrying out RNA-seq to identify any changes in splice-specific gene expression. This will allow us to determine whether increased PRPF6 expression is sufficient to drive alternative splicing changes.

  Study EGAS00001000761

• Assessing the impact of low frequency coding variants on disease risk using the Exomechip

  Following several rounds of Genome Wide Association (GWA) scans and subsequent meta-analyses of results multiple risk loci have been identified for many common diseases. However, in most instances the implicated common variants identified so far explain only a modest fraction of the genetic risk. In parallel to our efforts to identify the causative variants in the known loci it is necessary to assess the full spectrum of sequence variants for disease risk and in particular low frequency and rare variants that have not been tested in a comprehensive way so far. To this end an international effort by investigators who have performed whole exome sequencing in circa 12,000 individuals have assembled a set of ~250,000 exonic variants of low frequency – defined as seen in at least two studies and a minimum of three individuals (non-sysnonymous) or two individuals for SNPs altering splice sites / stop codons. The SNP content from the exome studies was complemented with additional interesting SNP sets totalling 25,000 markers (GWA tag SNPs, grid of common variants, HLA, Mitochondrial, Y-chromosome etc) and were used to generate a custom iSELECT array, the exome chip. Adequately powered association studies to test low frequency variants for association to disease risk are still very expensive if conducted by whole exome sequencing. The exome chip provides a cost efficient way to undertake such an experiment and the Wellcome Trust Cases Control Consortium will be applying this approach to eight diseases, type 1 and type 2 diabetes, coronary artery disease, hypertension, multiple sclerosis, rheumatoid arthritis, bipolar, and ankylosing spondylitis. The first objective will be to generate a large set of at least 10,000 common UK controls as analysis of lower frequency variants will require large sample sizes. This includes ~6000 samples from the 1958 Birth Cohort.

  Study EGAS00001000584

• Whole Genome Sequencing of Asian Lung Cancers: Second Hand Smoke is Not Responsible for Higher Incidence of Lung Cancer Among Asian Never-Smokers

  Background Asian nonsmoking populations have a higher incidence of lung cancer compared to their European counterparts. There is a long-standing hypothesis that the increase of lung cancer in Asian never-smokers is due to environmental factors such as second-hand smoke. Results We analyzed whole-genome sequencing of 30 Asian lung cancers. Unsupervised clustering of mutational signatures separated the patients into two categories: i) all of the never-smokers, and ii) only smokers or ex-smokers. Half of the ex-smokers / smokers were in the never-smoker-like cluster. The overall somatic variant profiles of Asian lung cancers were similar to that of European origin with G.C>T.A being predominant in smokers. We found EGFR and TP53 are the most frequently mutated genes with mutations in 50% and 27% of individuals, respectively. Among these Asian never-smokers, 71% had an EGFR mutation compared to 20% of the smokers in the smoking cluster. Other frequently mutated genes include RYR2, SATB2, C1orf88, FERMT1 and CTNNB1. Somatic alterations occurred in WNT signaling pathway genes, suggesting this pathway plays a role in both Western and Asian lung cancers. Conclusions Asian never-smokers have lung cancer signatures distinct from the smoker signature and their mutation profiles are similar to that of European never-smokers. The profiles of Asian and European smokers are also similar. This suggests that the same mutational mechanisms underlie the etiology for both ethnic groups, and the high incidence of lung cancer in Asian never-smokers might not be due to second hand smoke or other carcinogens that cause oxidative DNA damage. Half of the ex-smokers/smokers have a molecular phenotype similar to never-smokers; of which, 50% had EGFR mutations. This suggests that routine EGFR testing is warranted in the Asian population, regardless of smoking status.

  Study EGAS00001000621

• Probabilistic modeling of personalized drug combinations from integrated chemical screen and molecular data in sarcoma

  Background: Cancer patients with advanced disease exhaust available clinical regimens and lack actionable genomic medicine results, leaving a large patient population without effective treatments options when their disease inevitably progresses. To address the unmet clinical need for evidence-based therapy assignment when standard clinical approaches have failed, we have developed a probabilistic computational modeling approach which integrates sequencing data with functional assay data to develop patient-specific combination cancer treatments. Methods: Tissue taken from a murine model of alveolar rhabdomyosarcoma was used to perform single agent drug screening and DNA/RNA sequencing experiments; results integrated via our computational modeling approach identified a synergistic personalized two-drug combination. Cells derived from the tumor were allografted into mouse models and used to validate the personalized two-drug combination. Computational modeling of single agent drug screening and RNA sequencing of multiple heterogenous sites from a single patient’s epithelioid sarcoma identified a personalized two-drug combination effective across all tumor regions. The heterogeneity-consensus combination was validated in a xenograft model derived from the patient’s primary tumor. Cell cultures derived from human and canine undifferentiated pleomorphic sarcoma were assayed by drug screen; computational modeling identified a resistance-abrogating two-drug combination common to both cell cultures. This combination was validated in vitro via a cell regrowth assay. Results Our computational modeling approach addresses three major challenges in personalized cancer therapy: synergistic drug combination predictions (validated in vitro and in vivo in a genetically engineered murine cancer model), identification of unifying therapeutic targets to overcome intra-tumor heterogeneity (validated in vivo in a human cancer xenograft), and mitigation of cancer cell resistance and rewiring mechanisms (validated in vitro in a human and canine cancer model). Conclusions These proof-of-concept studies support the use of an integrative functional approach to personalized combination therapy prediction for the population of high-risk cancer patients lacking viable clinical options and without actionable DNA sequencing-based therapy.

  Study EGAS00001003564

• McGill Epigenomics Mapping Centre

  The McGill Epigenomics Mapping Centre (EMC) and Data Coordination Centre (EDCC) were established in 2012 at the McGill University and Genome Quebec Innovation Centre (Montreal, Canada) to support large-scale human epigenome mapping for a broad spectrum of cell types and diseases. Previous epigenetic studies have provided profiles of specific chromatin marks in relation to basic biological processes, and future studies of molecular mechanisms must build upon these proofs-of-concept by integrating sequence-based variation with multiple levels of epigenetic and transcriptional regulation across the genome of human tissues and animal disease models. This project leverages the high-throughput sequencing infrastructure and expertise in genomics and transcriptomics at the McGill Innovation Centre to carry out this research. Controlled access to the data by collaborators and the greater scientific community is gained via a portal, which takes advantage of Compute Canada high-performance computing cluster resources to manage the large volume of data associated from the generation of reference epigenome maps. McGill University is one of two Canadian mapping and data coordination centres, the other based at the Michael Smith Genome Sciences Centre in Vancouver, British Columbia. The generation of comprehensive epigenome maps at McGill University is part of a larger international effort that is coordinated by the International Human Epigenome Consortium (IHEC), whose overall long-term objective is to determine the extent to which the epigenome shapes human populations over generations and in response to the environment. This project is funded under the Canadian Epigenetics, Environment, and Health Research Consortium (CEEHRC) by the Canadian Institutes of Health Research and by Genome Quebec, with additional support from Genome Canada. The computing and networking infrastructure, and part of the software development, are provided by Compute Canada and CANARIE.

  Study EGAS00001000995

• Polymorphisms in the mitochondrial genome are associated with bullous pemphigoid in Germans

  Bullous pemphigoid (BP) is the most prevalent autoimmune skin blistering disease and is characterized by the generation of autoantibodies against the hemidesmosomal proteins BP180 (type XVII collagen) and BP230. Most intriguingly, BP is distinct from other autoimmune diseases because it predominantly affects elderly individuals above the age of 75 years, raising the question why autoantibodies and the clinical lesions of BP emerges mostly in this later stage of life, even in individuals harboring known putative BP-associated germline gene variants. The mitochondrial genome (mtDNA) is a potential candidate to provide additional insights into the BP etiology; however, the mtDNA has not been extensively explored to date. Therefore, we sequenced the whole mtDNA of German BP patients (n=180) and age- and sex-matched healthy controls (n=188) using next generation sequencing (NGS) technology, followed by the replication study using Sanger sequencing of an additional independent BP (n=89) and control cohort (n=104). While the BP and control groups showed comparable mitochondrial haplogroup distributions, the haplogroup T exhibited a tendency of higher frequency in BP patients suffering from neurodegenerative diseases (ND) compared to BP patients without ND (p= 0.1448, Fisher’s exact test). A total of four single nucleotide polymorphisms (SNPs) in the mtDNA, namely, m.16263T>C, m.11914G>A, m.16051A>G, and m.16162A>G, were found to be significantly associated with BP based on the meta-analysis of our NGS data and the Sanger sequencing data (p=0.0017, p=0.0132, p=0.0129, and p=0.0076, respectively, Peto’s test). In summary, our study is the first to interrogate the whole mtDNA in BP patients and controls and to implicate multiple novel mtDNA variants in disease susceptibility. Studies using larger cohorts and more diverse populations are warranted to explore the functional consequences of the mtDNA variants identified in this study on immune and skin cells to understand their contributions to BP pathology.

  Study EGAS00001003932

• UK10K COHORT ALSPAC

  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 Avon Longitudinal Study of Parents and Children (ALSPAC) is a two-generation prospective study. Pregnant women living in one of three health districts in the former county of Avon with an expected delivery date between April 1991 and December 1992 were eligible to be enrolled in the study, and this formed the initial point of contact for the development of a large, family based resource. Information has been collected on the children and the mothers through retrieval of biological materials (e.g. antenatal blood samples, placentas), biological sampling (e.g. collection of cord blood, umbilical cord, milk teeth, hair, toenails, blood and urine), self-administered questionnaires, data extraction from medical notes, linkage to routine information systems and at repeat research clinics.

  Study EGAS00001000090

• Genetic_Heterogeneity_of_the_familial_gastric_neuroendocrine_tumors

  Gastric neuroendocrine tumors (gNETs) occur with an estimated frequency of 2 per 100,000 in the general population. Type I gastric neuroendocrine tumors (NETs) represent the 75% of gNTEs and arise from gastric enterochromaffin-like (ECL) cells. They have late age of onset and usually benigh course. Classically, hypergastrinemia in patients who have autoimmune atrophic gastritis, causes hyperplasia of gastric ECL cells that progresses into type I gastric NETs and parietal cell (PC) destruction. The genetic bases in families with this disease are unknown. We performed an exome sequencing study of an atypical aggressive familial gNETs case (with early age onset, nodal infiltrations and gastric adenocarcinomas) that followed a recessive model. We identified a deleterious mutation in homozygosis in the ATP4A gene, which encodes the proton pump responsible for acid secretion by gastric parietal cells. This mutation lead to achlorhydria first, and hypergastrinemia and gNET developing as consequence (Calvete et al. 2014). Recently, two more families with gNETs, classical clinical traits and recessive model have been studies by WES but we didn't find any mutation in the ATP4a gene. However, putative mutations affecting genes that contribute to the development and the integrity of PC have been found suggesting that genetic alterations associated to this disorder target to a unique cell type (parietal cells). In order to cinfirm this hypothesis, it is necessary the search for new genes implicated in the gNETs, more familial cases are needed to be studied. We have identified four more new familial gNETs cases. Here, we propose their study by WES. The first family is formed by thress siblings with gNETs. The other families include two siblings with gNETs. 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 EGAS00001002273

• Genome_Diversity_in_Africa_Project__Uganda

  Globally, human populations show structured genetic diversity as a result of geographical dispersion, selection and drift. Understanding this genetic variation can provide insights into the evolutionary processes that shape both human adaptation and variation in disease. Populations from SSA have the highest levels of genetic diversity. This characteristic, in addition to historical genetic admixture, can lead to complexities in the design of studies assessing the genetic determinants of disease and human variation. However, such studies of African populations are also likely to provide new opportunities to discover novel disease susceptibility loci and variants and refine gene–disease association signals. A systematic assessment of genetic diversity within SSA would facilitate genomic epidemiological studies in the region. The Genome Diversity in Africa Project (GDAP) aims to produce a comprehensive catalogue of human genetic variation in SSA, including single nucleotide polymorphisms (SNPs), structural variants, and haplotypes. This resource will make a substantial contribution to understanding patterns of genetic diversity within and among populations in SSA, as well as providing a global resource to help design, implement and interpret genomic studies in SSA populations and studies comprising globally diverse populations, complementing existing genomic resources. Specifically, we plan to carry out high depth whole genome sequencing of up to 2000 individuals across Africa (25 individuals from each ethnolinguistic group). Our scientific objectives are to: 1) develop a resource that provides a comprehensive catalogue of genetic variation in populations from SSA accessible to the global scientific community; 2) characterise population genetic diversity, structure, gene flow and admixture across SSA; 3) develop a cost-efficient, next-generation genotype array for diverse populations across SSA; and 4) facilitate whole genome-sequencing association studies of complex traits and diseases by developing a reference panel for imputation and resource for enhancing fine-mapping disease susceptibility loci. These scientific objectives will be supported by cross-cutting operational activities, including network and management of the consortium, research ethics, and research capacity building in statistical genetics and bioinformatics

  Study EGAS00001002523

• Next-Generation Sequencing of RNA and DNA Isolated from Paired Fresh-Frozen and Formalin-Fixed Paraffin-Embedded Samples of Human Cancer and Normal Tissue

  Formalin-fixed, paraffin-embedded (FFPE) tissues are an invaluable resource for clinical research. However, nucleic acids extracted from FFPE tissues are fragmented and chemically modified making them challenging to use in molecular studies. We analysed 23 fresh-frozen (FF), 35 FFPE and 38 paired FF/FFPE specimens, representing six different human tissue types (bladder, prostate and colon carcinoma; liver and colon normal tissue; reactive tonsil) in order to examine the potential use of FFPE samples in next-generation sequencing (NGS) based retrospective and prospective clinical studies. Two methods for DNA and three methods for RNA extraction from FFPE tissues were compared and were found to affect nucleic acid quantity and quality. DNA and RNA from selected FFPE and paired FF/FFPE specimens were used for exome and transcriptome analysis. Preparations of DNA Exome-Seq libraries was more challenging (29.5 % success) than that of RNA-Seq libraries, presumably because of modifications to FFPE tissue-derived DNA. Libraries could still be prepared from RNA isolated from two-decade old FFPE tissues. Data were analysed using the CLC Bio Genomics Workbench and revealed systematic differences between FF and FFPE tissue-derived nucleic acid libraries. In spite of this, pairwise analysis of DNA Exome-Seq data showed concordance for 70-80 % of variants in FF and FFPE samples stored for fewer than three years. RNA-Seq data showed high correlation of expression profiles in FF/FFPE pairs (Pearson Correlations of 0.90 +/- 0.05), irrespective of storage time (up to 244 months) and tissue type. A common set of 1,494 genes was identified with expression profiles that were significantly different between paired FF and FFPE samples irrespective of tissue type. Our results are promising and suggest that NGS can be used to study FFPE specimens in both prospective and retrospective archive-based studies in which FF specimens are not available.

  Study EGAS00001000737

• Dynamics of multiple resistance mechanisms in plasma DNA and their clinical implications for NSCLC patients receiving EGFR-targeted therapies

  Background:The development of resistance mechanisms to targeted therapies is complicated by underlying tumour heterogeneity. When multiple resistance mechanisms co-exist, identifying the dominant driver(s) is critical for selection of treatment. Patients and methods:We studied the relative dynamics of multiple oncogenic drivers and resistance mechanisms in 50 patients with EGFR-mutant non-small cell lung cancer (NSCLC), during treatment with gefitinib and hydroxychloroquine (NCT00809237). We performed digital PCR and targeted deep sequencing of cell-free DNA from longitudinal plasma samples from all patients, and shallow whole genome sequencing of serial samples from three patients who underwent histological transformation to small-cell lung cancer (SCLC). ResultsEGFR activating mutations were assessed in tumour samples and accurately identified in plasma samples of 95% of patients (41/43). We identified additional mutations including EGFR T790M (31/50, 62%), TP53 (23/50, 46%), PIK3CA (7/50, 14%), and PTEN (4/50, 8%), and tracked their relative levels in plasma. Patients with both TP53 and EGFR mutations detected in plasma pre-treatment tended to have worse overall survival than those where only EGFR was detected. The resistance-conferring EGFR T790M mutation was identified in 62% of the patients who progressed. Changes in relative levels of T790M and activating mutations in EGFR corresponded to clinical responses in two patients who were given sequential treatments. Patients who progressed without T790M detected in plasma had worse PFS during TKI continuation, and developed alternative resistance mechanisms, that could be identified in plasma DNA, including SCLC-associated copy number changes and TP53 mutations, that tracked with progression and response to subsequent therapies. Conclusion:Longitudinal analysis of multiple oncogenic drivers in ctDNA may help identifying dominant mechanisms of resistance to EGFR-targeted therapies, and highlight the importance of monitoring genetic events that are not direct drug targets but provide non-invasive molecular information that may guide clinical management.

  Study EGAS00001002908

• Population Structure and Genetic Diversity in Argentinean populations

  Argentina was one of the last places peopled when humans entered the Americas. Its earliest archeological sites are dated from 13,000 years BP, and throughout its history the region has exhibited broad cultural diversity, involving both hunter-gatherer and agro-pastoral societies, different linguistic families, the Southern expansion of the Inca Empire, the arrival of the Spanish conquistadors and slave trade from Africa. More recently, Argentina has had the highest proportion of influx of transoceanic migration between 1850 and 1940 of any Latin American country. Furthermore, the geographical distribution of immigrants was oriented towards the Pampa region of central-eastern Argentina and the Northeast, so different patterns of mixed ancestry throughout the country are expected. In this present work we aim to describe the history and structure of the Argentinean population. Thus far, we have collected DNA samples from 2904 participants from 15 dispersed regions in Argentina, from the Andean populations in the Northwest to the river lands of the Northeast. Our preliminary data includes genome-scale Illumina Exome Array 250K genotypes, plus next-generation sequencing data from uniformly spaced genomic regions targeted with a genotype-by- sequencing strategy (which covers 1.5% of the genome), for 375 participants from 9 of the 15 populations. We will identify proportions and patterns of European, African and Indigenous American ancestry in the Argentinean cohort via local ancestry inference. The first series of analyses will focus on genomic tracts of Indigenous American ancestry for fine-scale examination of sub-continental structure by comparing to a panel of 493 samples from 52 Indigenous American populations from throughout the Americas. Furthermore, this work will serve as a test of different low-cost strategies (exome genotype array vs. genotype-by-sequencing) for accurate inference of population structure in extant populations.

  Study EGAS00001001663

• Copy number profiling of primary samples and cell lines of retinoblastoma

  While RB1 loss initiates retinoblastoma development, additional somatic copy number alterations (SCNAs) can drive tumor progression. Although SCNAs have been identified with good concordance between studies at a cytoband resolution, accurate identification of single genes for all recurrent SCNAs is still challenging. This study presents a comprehensive meta-analysis of genome-wide SCNAs integrated with gene expression profiling data, narrowing down the list of plausible retinoblastoma driver genes. We performed SCNA profiling of 45 primary retinoblastoma samples and 8 retinoblastoma cell lines by high-resolution microarrays. We combined our data with genomic, clinical and histopathological data of ten published genome-wide SCNA studies, which strongly enhanced the power of our analyses (N=310). Comprehensive recurrence analysis of SCNAs in all studies integrated with gene expression data allowed us to reduce candidate gene lists for 1q, 2p, 6p, 7q and 13q to a limited gene set. Besides the well-established driver genes RB1 (13q-loss) and MYCN (2p-gain) we identified CRB1 and NEK7 (1q-gain), SOX4 (6p-gain) and NUP205 (7q-gain) as novel retinoblastoma driver candidates. Depending on the sample subset and algorithms used, alternative candidates were identified including MIR181 (1q-gain) and DEK (6p gain). Remarkably, our study showed that copy number gains rarely exceeded change of one copy, even in pure tumor samples with 100% homozygosity at the RB1 locus (N=34), which is indicative for intra-tumor heterogeneity. In addition, profound between-tumor variability was observed that was associated with age at diagnosis and differentiation grades. Since focal alterations at commonly altered chromosome regions were rare except for 2p24.3 (MYCN), further functional validation of the oncogenic potential of the described candidate genes is now required. For further investigations, our study provides a refined and revised set of candidate retinoblastoma driver genes.

  Study EGAS00001001715

• NEOPREDICT-Lung: Neoadjuvant Immunotherapy for Resectable NSCLC

  Links via: DOI: https://doi.org/10.1038/s41591-024-02965-0 PubMed: https://pubmed.ncbi.nlm.nih.gov/38689060/ Journal: https://www.nature.com/articles/s41591-024-02965-0 Abstract Antibodies targeting the immune checkpoint molecules PD-1, PD-L1 and CTLA-4 administered alone or in combination with chemotherapy are standard of care in most patients with metastatic non-small cell lung cancers. When given prior to curative surgery, tumor responses and improved event-free survival are achieved. New antibody combinations may be more efficacious and tolerable. In an ongoing, open-label phase II study, 60 biomarker-unselected, treatment-naïve patients with resectable non-small cell lung cancer were randomized to receive two preoperative doses of nivolumab (anti-PD-1) with or without relatlimab (anti-LAG-3) antibody therapy. The primary study endpoint was feasibility of surgery within 43 days, which was met by all patients. Curative resection was achieved in 95% of patients. Secondary endpoints included pathological and radiographic response rates, pathologically complete resection rates, disease-free and overall survival rates, and safety. Major pathological (≤10% viable tumor cells) and objective radiographic responses were achieved in 27% and 10% (nivolumab), and 30% and 27% (nivolumab and relatlimab) of patients, respectively. In 100% (nivolumab) and 90% (nivolumab and relatlimab) of patients tumors and lymph nodes were pathologically completely resected. With 12 months median duration of follow-up, disease-free survival and overall survival rates at 12 months were 89% and 93% (nivolumab), and 93% and 100% (nivolumab and relatlimab). Both treatments were safe with grade ≥ 3 treatment-emergent adverse events reported in 10% and 13% of patients per study arm. Exploratory analyses provided insights into biological processes triggered by preoperative immunotherapy. This study establishes the feasibility and safety of dual targeting of PD-1 and LAG-3 prior to lung cancer surgery.

  Study EGAS00001007753

• The G2 gene expression signature and MYC overexpression are independent poor prognostic factors in childhood high-grade osteosarcoma

  Purpose In high-grade osteosarcoma, prognostic factors at diagnosis are insufficient for stratifying patients into relevant subgroups. Recently, a transcriptomic study developed the G1/G2 gene expression signature, in which the G2 signature was associated with unfavorable survival. An orthogonal study identified MYC amplification as an unfavorable prognostic factor using targeted next-generation sequencing. The purpose of this study was to validate the independent prognostic value and to investigate the combined prognostic value of the G1/G2 signature with, MYC amplification and/or MYC expression for survival prediction. Material and methods This study included pediatric and adolescent patients with high-grade osteosarcoma. RNA sequencing was performed in 48 patients. Whole-exome sequencing was performed in 40 patients. Gene expression signature scores, MYC amplification (defined as > 7 copies) and MYC expression levels were calculated. Multivariable Cox proportional hazards analysis was performed for event-free survival (EFS) (primary endpoint) and overall survival (OS) (secondary endpoint). Results In the full cohort, the 3-year event-free survival rate was 37%. In multivariable Cox regression analysis with metastatic disease stage (n=21, 44%) as covariate, the G2 signature and MYC expression were independently associated with worse outcomes in terms of EFS (Hazard ratio (HR) 3.32, 95% confidence interval (CI) 1.34-8.21 and HR 3.38, 95%-CI 1.71-6.66, respectively) and OS (HR 4.07, 95%-CI 1.19-13.9 and HR 2.88, 95%-CI 1.22-6.76, respectively). MYC amplification was not associated with EFS or OS in univariable analysis (HR 1.88, 95% -CI 0.74-4.77 and HR 0.79, 95%-CI 0.21-3.05, respectively). Conclusion The G2 gene expression signature and MYC expression were independently associated with unfavorable outcomes in a pediatric cohort of patients with high-grade osteosarcoma. The combined prognostic value warrants further prospective validation and could potentially serve as a stratification marker for future osteosarcoma treatment protocols.

  Study EGAS00001008073

• Whole exome and RNA sequencing of organoid samples derived from TRACERx patients

  Immune checkpoint blockade in non-small cell lung cancer (NSCLC) rarely leads to complete tumour remission, suggesting that tumours harbour subpopulations of cancer cells that can escape immune predation. However, while tumours consist of multiple genetically distinct clones, the extent to which separate tumour subclones differ intrinsically in their functional capacity for immune evasion remains largely unexplored. This has been challenging to address due to an inability to isolate and propagate individual tumour subclones from human cancers. Here, we leverage the multi-region TRACERx lung cancer evolution study to generate a patient-derived platform that allows the functional analysis of immune escape at single clone resolution. We generated libraries of >20 separate NSCLC organoid lines per patient by isolating individual (clonal) organoids established from multiple spatially separated tumour regions. Each organoid subline was co-cultured with autologous tumour infiltrating lymphocytes (TIL) to evaluate their capacity to elicit a T cell response. Organoid lines derived from separate tumour regions, or from single clones within individual regions, differ vastly in their capacity to elicit a T cell response, showing the coexistence of intrinsically immune evading (‘hot’) and non-evading (‘cold’) clones in the same primary tumour. Organoid lines escaping T cell recognition represent genetically and transcriptionally unique subclones with a distinct evolutionary history, showing that evolutionary sculpting of subclonal diversity can have functional consequences for anti-tumour immunity. Notably, the identification of immune-evasive subclones would not have been possible using computational approaches alone, highlighting the unique value of functional co-culture platforms. This patient-derived functional co-culture platform reveals that tumour evolution can give rise to distinct cancer clones with intrinsic differences in immune evasion capacity and provides an approach to prospectively identify, isolate and characterise immune evading subclones from patients with cancer.

  Study EGAS00001008092

• Raw bulk TCRseq data fromTIL and non-TIL expanded TRACERx samples in the manuscript Bulk TCRseq from TRACERx samples from 'Subclonal immune evasion in non-small cell lung cancer'

  Immune checkpoint blockade in non-small cell lung cancer (NSCLC) rarely leads to complete tumour remission, suggesting that tumours harbour subpopulations of cancer cells that can escape immune predation. However, while tumours consist of multiple genetically distinct clones, the extent to which separate tumour subclones differ intrinsically in their functional capacity for immune evasion remains largely unexplored. This has been challenging to address due to an inability to isolate and propagate individual tumour subclones from human cancers. Here, we leverage the multi-region TRACERx lung cancer evolution study to generate a patient-derived platform that allows the functional analysis of immune escape at single clone resolution. We generated libraries of >20 separate NSCLC organoid lines per patient by isolating individual (clonal) organoids established from multiple spatially separated tumour regions. Each organoid subline was co-cultured with autologous tumour infiltrating lymphocytes (TIL) to evaluate their capacity to elicit a T cell response. Organoid lines derived from separate tumour regions, or from single clones within individual regions, differ vastly in their capacity to elicit a T cell response, showing the coexistence of intrinsically immune evading (‘cold) and non-evading (‘hot) clones in the same primary tumour. Organoid lines escaping T cell recognition represent genetically and transcriptionally unique subclones with a distinct evolutionary history, showing that evolutionary sculpting of subclonal diversity can have functional consequences for anti-tumour immunity. Notably, the identification of immune-evasive subclones would not have been possible using computational approaches alone, highlighting the unique value of functional co-culture platforms. This patient-derived functional co-culture platform reveals that tumour evolution can give rise to distinct cancer clones with intrinsic differences in immune evasion capacity and provides an approach to prospectively identify, isolate and characterise immune evading subclones from patients with cancer.

  Study EGAS00001008118

• Successful immune checkpoint blockade in a patient with advanced stage microsatellite unstable biliary tract cancer (H021)

  Cancers acquire multiple somatic mutations that can lead to the generation of immunogenic mutation-induced neoantigens. These neoantigens can be recognized by the host’s immune system. However, continuous stimulation of immune cells against tumor antigens can lead to immune cell exhaustion, which allows uncontrolled outgrowth of tumor cells. Recently, immune checkpoint inhibitors have emerged as a novel approach to overcome immune cell exhaustion and re-activate anti-tumor immune responses. In particular, antibodies blocking the exhaustion-mediating programmed death receptor (PD-1)/PD-L1 pathway have shown clinical efficacy. The effects were particularly pronounced in tumors with DNA mismatch repair deficiency and a high mutational load, which typically occur in the colon and endometrium. Here, we report on a 24-year old woman diagnosed with extrahepatic cholangiocarcinoma who showed strong and durable response to the immune checkpoint inhibitor pembrolizumab, although treatment was initiated at an advanced stage of disease. The patient’s tumor displayed DNA mismatch repair deficiency and microsatellite instability (MSI), but lacked other features commonly discussed as predictors of response towards checkpoint blockade, such as PD-L1 expression or dense infiltration with cytotoxic T cells. Notably, high levels of HLA class I and II antigen expression were detected in the tumor, suggesting a potential causal relation between functionality of the tumor’s antigen presentation machinery and the success of immune checkpoint blockade. We suggest determining MSI status in combination with HLA class I and II antigen expression in tumors potentially eligible for immune checkpoint blockade even in the absence of conventional markers predictive for anti-PD-1/PD-L1 therapy and in entities not commonly linked to MSI phenotype. Further studies are required to determine the value of these markers for predicting the success of immune checkpoint blockade. (H021)

  Study EGAS00001002441

• Exome Sequencing to Define the Landscape of Plasma Cells in Systemic Light chain Amyloidosis

  Systemic light chain amyloidosis (AL) is characterized by the deposition of immunoglobulin light chains as amyloid fibrils in different organs, where they form toxic protein aggregates. The underlying disease is a plasma cell disorder, but limited whole exome data are available. We report the findings of an exome sequencing study in AL to define a plasma cell signature and compare this to monoclonal gammopathy of undefined significance (MGUS) and myeloma (MM). Twenty-four samples from unselected newly diagnosed untreated AL patients were analysed. CD138+ cells were isolated from bone marrow cells using MACSorting (Miltenyi Biotech, Bisley, UK). Cells were lysed in RLT+ buffer and DNA/RNA extracted using the AllPrep kit (Qiagen, Manchester, UK). Peripheral blood was isolated from patients, white blood cells purified by Ficoll-Pacque and DNA extracted using the QIAamp DNA mini kit (Qiagen). There were a median of 47 acquired nonsynonymous variants (range 4-253) per patient. Thirty percent of patients yielded a mutation in a myeloma driver gene. These mutations were either clonal or sub-clonal. There was evidence of MAPK activation with NRAS mutations, NF-κB activation and DNA repair pathway alterations. Exome data were used to determine the cytogenetic groups of AL samples and identified hyperdiploidy in 33% of cases and t(11;14) in 33% of cases. All translocations were generated via class switch recombination and breakpoints on chromosome 11 fell into the expected regions. A novel translocation, unlikely to be driving disease progression was identified involving RCC1. Finally, further analysis of the APCS gene encoding for the SAP protein did not show any biased genotype in the AL population in comparison to the general population. We conclude that exome sequencing identifies a genetic make-up of AL that is similar to other plasma cell disorders such as MM and MGUS. This not only includes copy number abnormalities and translocations but also nonsynonymous mutations.

  Study EGAS00001001418

• H3Africa - Genomic and Environmental Risk Factors for Cardiometabolic Disease in Africans

  The long-term vision of the AWI-Gen Collaborative Centre (CC) is to build sustainable capacity in Africa for research that leads to an understanding of the interplay between genetic, epigenetic and environmental risk factors for obesity and related cardiometabolic diseases (CMD) in sub-Saharan Africa. The CC will be consolidated under the auspices of the University of the Witwatersrand (Wits) and the International Network for the Demographic Evaluation of Populations and Their Health in Developing Countries (INDEPTH). It will capitalize on the unique strengths of existing longitudinal cohorts, including the urban Soweto and rural Agincourt studies in South Africa (Wits based), and the well established INDEPTH demographic health and surveillance centers in Kenya, Ghana, Burkina Faso and South Africa. The centers offer established infrastructure, trained fieldworkers, long-standing community engagement, and detailed longitudinal phenotypic data, focusing on obesity and cardiometabolic health. Key strengths are harmonized phenotyping across sites, building on strong existing cohorts, and representation of the geographic and social variability of African populations. We aim to: 1. Build sustainable infrastructure (biobanks and laboratories) and capabilities (well characterized population cohorts, genotyping and bioinformatics) for genomic research on the African continent; 2. Understand the genomic architecture of sub-Saharan populations from west, east and south Africa to guide genomic studies (genome sequencing and high throughput SNP and CNV arrays using unrelated individuals and family trios to improve the accuracy of haplotype analyses) and; 3. Investigate the independent and synergistic genomic contributions to body fat distribution (BMI, hip/waist circumference, subcutaneous and visceral fat) in these populations considering the relevant environmental and social contexts (rural/urban communities, quickly transitioning obesity prevalence, differential HIV, TB, and malaria infection histories). We will investigate the effect of obesity and fat distribution on the risk for CMD in the longitudinal cohorts.

  Study EGAS00001002482

• UK10K_NEURO_ASD_FI

  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.These samples are a subset of a nationwide collection of Finnish autism spectrum disorder (ASD) samples. The samples have been collected from Central Hospitals across Finland in collaboration with the University of Helsinki. The samples consist of 93 individuals with a diagnosis of autistic disorder or Asperger syndrome from 36 families with at least two affected individuals. Of these individuals, 16 can be genealogically connected to form two large pedigrees originating from Central Finland, suggesting possible genetic risk factors shared identical by descent within the pedigrees. All diagnoses are based on ICD-10 and DSM-IV diagnostic criteria for ASDs. Additional phenotypic data is available for a subset of the individuals.For further information with regard to this cohort please contact Aarno Palotie (Aarno.palotie@helsinki.fi).

  Study EGAS00001000110

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

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

  Study EGAS00001002676

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

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

  Study EGAS00001002684

• Intra-tumor heterogeneity of localized lung adenocarcinomas defined by multi-region sequencing

  Background: Intra-tumor heterogeneity (ITH) may have profound impacts on biopsy strategy, characterization of actionable targets, treatment planning, and drug resistance. ITH of lung cancers has not been well studied. Methods: We performed multi-region whole exome sequencing (average sequencing depth 277x) on 48 tumor regions from 11 resected lung adenocarcinomas (8 stage I, 2 stage II and one stage III tumors) followed by deep sequencing validation (average sequencing depth 863x). Results: All tumors showed clear evidence of spatial ITH. On average, 76% of all mutations and 20/21 known cancer gene mutations were present in all regions of individual tumors (trunk mutations). Copy number analysis did not show substantial difference in large-scale chromosome aberrations. A canonical KRAS p.G12C mutation was detected in only 1 of 4 tumor regions (branch mutations) at exome depth but detected in all four tumor regions (trunk mutations) at increased sequencing depth. In addition, many heterogeneous branch mutations defined by exome sequencing were detectable in all regions of individual tumors with increasing sequencing depth. Significant differences in mutational spectrum were observed in 6 tumors indicating that specific mutational processes were likely operative at different times during development of these tumors. With a median follow up of 21 months post surgery, 3 patients have relapsed and all 3 patients presented with significantly larger fractions of subclonal mutations in their primary tumors. Conclusions: Although multi-region sampling is still required to fully assess ITH complexity, single biopsy analysis at appropriate depth might be sufficient to identify the majority of known cancer gene mutations in localized lung adenocarcinomas. Larger subclonal mutation fractions may be associated with increased likelihood of postsurgical relapse in patients with localized lung adenocarcinomas. Studies on larger cohort, ideally with comprehensive clinical annotation and repeat biopsy at relapse are needed to fully understand the clinical impact of ITH.

  Study EGAS00001000930

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

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

  Study EGAS00001000673

• Whole exome sequencing of germline DNA was performed and subsequent polymorphisms in genes known and putatively involved in the innate immune response to fungi were identified

  Disease can be conceptualized as the result of interactions between infecting microbe and holobiont, the combination of a host and its microbial communities. It is likely that genomic variation in the host, infecting microbe, and commensal microbiota are key determinants of infectious disease clinical outcomes. However, until recently, simultaneous, multiomic investigation of infecting microbe and holobiont components has rarely been explored. Herein, we characterized the infecting microbe, host, micro- and mycobiomes leading up to infection onset in a leukemia patient that developed invasive mucormycosis. We discovered that the patient was infected with a strain of the recently described Mucor velutinosus species which we determined was hypervirulent in a Drosophila challenge model and has a predisposition for skin dissemination. After completing the infecting M. velutinosus genome and genomes from four other Mucor species, comparative pathogenomics was performed and assisted in identifying 66 M. velutinosus-specific putatively secreted proteins, including multiple novel secreted aspartyl proteinases which may contribute to the unique clinical presentation of skin dissemination. Whole exome sequencing of the patient revealed multiple non-synonymous polymorphisms in genes critical to control of fungal proliferation, such as TLR6 and PTX3. Moreover, the patient had a non-synonymous polymorphism in the NOD2 gene and a missense mutation in FUT2, which have been linked to microbial dysbiosis and microbiome diversity maintenance during physiologic stress, respectively. In concert with host genetic polymorphism data, the micro- and mycobiome analyses revealed that the infection developed amid a dysbiotic microbiome with low α-diversity, dominated by staphylococci. Additionally, longitudinal mycobiome data showed that M. velutinosus DNA was detectable in oral samples preceding disease onset. Our genome-level study of the host-infecting microbe-commensal triad extends the concept of personalized genomic medicine to the holobiont-infecting microbe interface thereby offering novel opportunities for using synergistic genetic methods to increase understanding of infectious diseases pathogenesis and clinical outcomes.

  Study EGAS00001001542

• Bulk-RNA Sequencing of high-grade pancreatic and non-pancreatic Neuroendocrine Neoplasms

  Therapeutic decisions in oncology depend on a precise pathological classification of individual neoplasms. Recent years have seen an intensification of research activities aimed at the extraction of clinically relevant information from patient-derived 'omics' data based on Machine-Learning models. However, a comprehensive training of Machine-Learning models requires sufficiently large numbers of training samples, which are usually not available for rare cancer types. The problem is worsened when individual tissues segregate into different cancer subtypes, as their discrimination would require even more training samples. Methods: Here, we report on a new data-augmentation technique to support the training of Machine-Learning models on ‘omics’ data from pancreatic neuroendocrine neoplasms (panNEN). PanNENs display all properties described above: Only about 2-3% of all pancreatic neoplasms are neuroendocrine and they fall into different subtypes with distinctly different prognosis, which makes the precise classification of such samples both difficult and important for therapy decisions. The approach reconstructs a given transcriptome based on healthy pancreatic cell type signatures and creates a Machine-Learning model that integrates the observed reconstruction error and predicted cell type proportions as features. Results: A benchmark of the deconvolution model predictions with the ground-truth found that the model could efficiently predict the sample grading, disease-related patient survival time, and differentiate between different subtypes in four panNEN and one mixed panNEN and non-pancreatic NENs dataset. We compared the predictive performance of the deconvolution-trained model to that of a model trained directly on the transcriptomic data, under inclusion of the Ki-67 classification gold standard biomarker, and found the performances to be comparable. Conclusion: Our approach serves as a data-augmentation technique to facilitate the training of Machine Learning models for rare cancer types via utilization of data from healthy origin. Additionally, the deconvolution results were clinically interpretable and further research can render the deconvolution approach an effective complementary asset for the clinical classification of neoplasms.

  Study EGAS00001004861

• Molecular Tumor Board to Inform on Personalized Medicine for a Man with Advanced Prostate Cancer

  Purpose: Molecular profiling of cancer is increasingly common as part of routine care in oncology, and germline and somatic profiling may provide insights and actionable targets for men with metastatic prostate cancer. However, all reported cases are of deidentified individuals without full medical and genomic data available in the public domain. Patient and Methods: We present a case of whole genome tumor and germline sequencing in a patient with advanced prostate cancer, who has agreed to make his genomic and clinical data publicly available in an accompanying personal testimonial. Results: We describe an 84-year-old Caucasian male with a Gleason 10 oligometastastic hormone-sensitive prostate cancer. Whole-genome sequencing provided insights into his tumor’s underlying mutational processes and the development of an SPOP mutation. It also revealed an androgen-receptor dependency of his cancer which was reflected in his durable response to radiation and hormonal therapy. Potentially actionable genomic lesions in the tumor were identified through a personalized medicine approach for potential future therapy, but at the moment, he remains in remission, illustrating the hormonal sensitivity of his SPOP-driven prostate cancer. We also placed this patient in the context of a large prostate-cancer cohort from the PCAWG (Pan-cancer Analysis of Whole Genomes) group. In this comparison, the patient's cancer appears typical in terms of the number and type of somatic mutations, but it has a somewhat larger contribution from the mutational process associated with aging. Conclusion: We combined the expertise of medical oncology and genomics approaches to develop a molecular tumor board to integrate the care and study of this patient, who continues to have an outstanding response to his combined modality treatment. This identifiable case potentially helps overcome barriers to clinical and genomic data sharing.

  Study EGAS00001004648

• Genetic analysis of HLA and immune escape genes in Diffuse Large B-cell Lymphoma

  Fifty percent of diffuse large B-cell lymphoma (DLBCL) cases lack cell-surface expression of the class-I major histocompatibility complex (MHC-I), thereby escaping immune recognition by cytotoxic T cells. In order to comprehensively identify the mechanisms involved in MHC-I loss, we first performed immunophenotypic analysis of both MHC class-I and -II in 657 cases across the spectrum of B-NHL, revealing that loss of MHC-I, but not MHC-II, is preferentially restricted to DLBCL. We then used whole exome and targeted deep-sequencing to examine genes involved in MHC-I expression in 74 DLBCL samples representative of MHC-I positive and negative cases. We show here that somatic biallelic or monoallelic inactivation of B2M and/or HLA-I is present in 80.9% (34/42) of MHC-I negative tumors. Furthermore, 68.8% (22/32) of MHC-I positive DLBCLs also harbored monoallelic HLA-I genetic alterations (MHC-I positivemono) that lead to allelic imbalance, suggesting allele-specific inactivation. Both MHC-I negative and MHC-I positivemono cases showed a significantly higher mutational burden as well as a higher number of inferred neo-antigens, suggesting co-selection of HLA-I loss and sustained neo-antigen production. Interestingly, the analysis of > 500.000 individuals in two databases revealed homozygosity of germline HLA-I genes in 26-38% of DLBCL patients, a frequency significantly higher than that observed in any other cancer type. In mice, germinal-center B cells lacking HLA-I expression did not progress to lymphoma and were counterselected in the context of oncogene-driven lymphomagenesis, suggesting that additional events are needed to license immune evasion. These results suggest a multi-step process of HLA-I loss including both germ-line and somatic events in DLBCL development, and have direct implications for the pathogenesis and immunotherapeutic targeting of this disease.

  Study EGAS00001005054

• Same-day genomic and epigenomic diagnosis of brain tumors using realtime nanopore sequencing

  Molecular classification of cancer has entered clinical routine to inform diagnosis, prognosis and treatment decisions. At the same time, new tumor entities have been identified that cannot be defined histologically. For central nervous systems tumors, the current World Health Organization classification explicitly demands molecular testing, e.g. for 1p/19q- codeletion or IDH mutations, to make an integrated histomolecular diagnosis. However, a plethora of sophisticated technologies is currently needed to assess different genomic and epigenomic alterations and turnaround times are in the range of weeks, which makes standardized and widespread implementation difficult and hinders timely decision making. Here, we explored the potential of a pocket-size nanopore sequencing device for multimodal and rapid molecular diagnostics of cancer. Low-pass whole genome sequencing was used to simultaneously generate copy number (CN) and methylation profiles from native tumor DNA in the same sequencing run. Single nucleotide variants in IDH1, IDH2, TP53, H3F3A and the TERT promoter region were identified using deep amplicon sequencing. Nanopore sequencing yielded ~0.1X genome coverage within six hours and resulting CN and epigenetic profiles correlated well with matched microarray data. Diagnostically relevant alterations, such as 1p/19q codeletion, and focal amplifications could be recapitulated. Using ad hoc random forests, we could perform supervised pan-cancer classification to distinguish gliomas, medulloblastomas and brain metastases of different primary sites. Single nucleotide variants in IDH1, IDH2 and H3F3A were identified using deep amplicon sequencing within minutes of sequencing. Detection of TP53 and TERT promoter mutations shows that sequencing of entire genes and GC-rich regions is feasible. Nanopore sequencing allows same-day detection of structural variants, point mutations and methylation profiling using a single device with negligible capital cost. It outperforms hybridization-based and current sequencing technologies with respect to time-to diagnosis and required laboratory equipment and expertise, aiming to make precision medicine possible for every cancer patient, even in resource restricted settings.

  Study EGAS00001002213

• Deep single-cell RNA sequencing data for 11138 T cells from tumour, adjacent normal tissue and peripheral blood of treatment-naive CRC patients. The DATA ACCESS AGREEMENT is provided at https://github.com/zhangyybio/single-T-cell-data-access. Applicants can request access to the data by directly downloading it or by sending an email to cancerpku@pku.edu.cn. The process that is used to approve an application includes verifying the institution, participants and research purposes of the application. In general this process will take about two weeks. In principal, any scientific research program complying with the laws and bioethic regulation policies of China will be approved.

  T cells are central players in cancer immunotherapy1, yet some of their fundamental properties such as development and migration within tumours remain elusive. The enormous T cell receptor (TCR) repertoire, required for recognising foreign and self-antigens2,3, could serve as lineage tags to track these T cells in tumours4. Here, we obtained transcriptomes of 11,138 single T cells from 12 colorectal cancer (CRC) patients and developed STARTRAC (Single T-cell Analysis by Rna-seq and Tcr TRACking) indices to quantitatively analyse dynamic relationships among 20 identified T cell subsets with distinct functions and clonalities. While both CD8+ effector and “exhausted” T cells exhibited high clonal expansion, they were independently connected with tumour-resident CD8+ effector memory cells, implicating a TCR-based fate decision. Of the CD4+ T cells, the majority of tumour-infiltrating Tregs showed clonal exclusivity, whereas certain Treg clones were developmentally linked to multiple TH clones. Notably, we identified two IFNG+ TH1-like clusters in tumours, the GZMK+ TEM and CXCL13+ TH1-like clusters, which were associated with distinct IFN-γ-regulating transcription factors, EOMES/RUNX3 and BHLHE40, respectively. Only CXCL13+BHLHE40+ TH1-like cells were preferentially enriched in tumours of microsatellite-instable (MSI) patients, which might explain their favourable response rates to immune-checkpoint blockade. Furthermore, we found IGFLR1 to be highly expressed in both CXCL13+BHLHE40+ TH1-like and CD8+ exhausted T cells and possessed co-stimulatory functions. Our integrated STARTRAC analyses provides a powerful avenue to comprehensively dissect the T cell properties in CRC, which could shed new insights into the dynamic relationships of T cells in other cancers

  Study EGAS00001002791

• WES_of_adult_intellectual_disabilities_with_co_morbid_psychiatric_disorders

  ost adults with intellectual disabilities (ID) do not undergo genetic diagnostic investigation as part of their clinical care and have 'missed the boat' with regard to the WES and WGS genetic testing that is now being provided for children with ID. There is a dramatically increased risk of psychistric disorders in adults with ID, e.g. the risk of psychoses is 10X higher than in the general population. It remains an open question as to how much of adult ID is genetic in origin and how similar the genetic forms of adult ID are to those being diagnosed in children, in part due to survivor bias. There is also the opportunity to identify adults with treatable forms of ID, of which over 80 have been described, thus improving their clinical management. Furthermore, analysis of medical records of adults with genetic forms of ID can help to characterise the 'natural history' of individual disorders, resulting in more accurate prognoses for diagnosed children and identifying opportunities for improved management and possibly therapeutic intervention (e.g. optimal anti-epileptic therapy). Here we propose to exome sequence (to ~50X coverage) 200 adults with ID and co-morbid psychiatric disorders. This cohort has previously been assayed with chromosomal microarrays (Wolfe et al 2017 EJHG, 25, 66–72) identifying a diagnostic yield of ~11% which is comparable to the CNV diagnostic yield in various child ID cohorts (10-15%). The authors observed no substantive biases in diagnostic yield between different psychiatric diagnostic classes. The WES data will be analysed using the diagnostic workflows developed in the DDD study to ensure comparability between child and adult ID datasets. This study is intended as a pilot study to demonstrate the value of WES in adults with ID. 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 EGAS00001002962