Search Results - EGA European Genome-Phenome Archive

ctDNA monitoring using patient-specific sequencing and integration of variant reads - Lung cohort

Circulating tumor-derived DNA (ctDNA) can be used to monitor cancer dynamics noninvasively. Detection of ctDNA can be challenging in patients with low-volume or residual disease, where plasma contains very few tumor-derived DNA fragments. We show that sensitivity for ctDNA detection in plasma can be improved by analyzing hundreds to thousands of mutations that are first identified by tumor genotyping. We describe the INtegration of VAriant Reads (INVAR) pipeline, which combines custom error-suppression methods and signal-enrichment approaches based on biological features of ctDNA. With this approach, the detection limit in each sample can be estimated independently based on the number of informative reads sequenced across multiple patient-specific loci. We applied INVAR to custom hybrid-capture sequencing data from 176 plasma samples from 105 patients with melanoma, lung, renal, glioma, and breast cancer across both early and advanced disease. By integrating signal across a median of >105 informative reads, ctDNA was routinely quantified to 1 mutant molecule per 100,000, and in some cases with high tumor mutation burden and/or plasma input material, to individual parts per million. This resulted in median Area Under the Curve (AUC) values of 0.98 in advanced cancers, and 0.80 in early stage and challenging settings for ctDNA detection. We generalized this method to whole-exome and whole-genome sequencing, showing that the INVAR may be applied without requiring personalized sequencing panels, so long as a tumor mutation list is available. As tumor sequencing becomes increasingly performed, such methods for personalized cancer monitoring may enhance the sensitivity of cancer liquid biopsies.

Study EGAS00001004447

Deep single-cell RNA sequencing data for 12346 T cells from tumour, adjacent normal tissue and peripheral blood of treatment-naive NSCLC patients

Cancer immunotherapies have shown sustained clinical responses in treating non-small cell lung cancer (NSCLC), but efficacy varies between patients and is believed to depend in part on the amount and properties of tumor infiltrating lymphocytes (TILs). To comprehensively depict and dissect the baseline landscape of the composition, lineage and functional states of TILs in lung cancer, here we generated deep single-cell RNA sequencing data for 12,346 T cells from the primary tumour, adjacent normal tissues and peripheral blood of 14 treatment-naive NSCLC patients. Combined expression and TCR-based lineage tracking revealed a significant proportion of effector T cells with common origins and similar functional states across peripheral blood and tumours pointing towards a highly migratory nature of these T cells. We also observed tumour-infiltrating CD8+ T cells undergoing extensive clonal expansion and exhaustion, with two clusters of cells exhibiting states preceding exhaustion. Survival analysis on independent datasets suggested that high ratio of pre-exhausted to exhausted T cells is associated with better prognosis of lung adenocarcinoma (LUAD). In addition, we observed further heterogeneity within the tumour regulatory T cells (Tregs), characterized by the bimodal distribution of TNFRSF9, an activation marker for antigen-specific Tregs. The gene signature of this group of activated tumour Tregs, which included IL1R2, correlated with poor prognosis in LUAD. The T cell clusters revealed by our single cell analyses provide a new approach for patient stratification, and the accompanying compendium of data will help the research community to gain further insight into the functional states and dynamics of T cell responses in lung cancer.

Study EGAS00001002430

Biomarker Screen for Efficacy of Oncolytic Virotherapy in Patient-derived Pancreatic Cancer Cultures

Pancreatic ductal adenocarcinoma (PDAC) is a tumour entity with unmet medical need. To assess the therapeutic potential of oncolytic virotherapy (OVT) against PDAC, different oncolytic viruses (OVs) are currently investigated in clinical trials. However, systematic comparisons of these different OVs in terms of efficacy against PDAC and biomarkers predicting therapeutic response are lacking. We screened fourteen patient-derived PDAC cultures for their sensitivity to five clinically relevant OVs, namely serotype 5 adenovirus Ad5-hTERT, herpes virus T-VEC, measles vaccine strain MV-NIS, reovirus jin-3, and protoparvovirus H 1PV using live cell analysis, quantification of viral genome/gene expression, cell viability as well as cytotoxicity assays. Sensitivity was correlated with transcriptome profiles to identify potential predictive biomarkers for response to OVT. Patient-derived PDAC cultures showed individual response patterns to OV treatment. Twelve of fourteen cultures were responsive to at least one OV, with no single OV proving superior or inferior across all cultures. Known host factors for distinct viruses were retrieved as potential biomarkers. Compared to the basal-like molecular subtype, the quasi-mesenchymal subtype of PDAC was more sensitive to H-1PV, jin-3, and T-VEC. Expression of viral entry receptors did not correlate with sensitivity to OV treatment. Rather, cellular pathways controlling immunological, metabolic, and proliferative signalling appeared to determine outcome. For instance, high baseline expression of interferon-stimulated genes (ISGs) correlated with relative resistance to oncolytic measles virus, whereas cGAS expression was associated with exceptional response. Combination treatment of MV-NIS with a cGAS inhibitor improved tumour cell killing in several PDAC cultures.

Study EGAS00001007001

Whole-genome Sequencing Suggests Mechanisms for 22q11.2 deletion-associated Parkinson’s disease

Objectives: To investigate disease risk mechanisms of early-onset Parkinson’s disease (PD) associated with the recurrent 22q11.2 deletion, a genetic risk factor for early-onset PD. Methods: In a proof-of-principle study, we used whole-genome sequencing (WGS) to investigate sequence variants in nine adults with 22q11.2DS, three with neuropathologically confirmed early-onset PD and six without PD. Adopting an approach used recently to study schizophrenia in 22q11.2DS, here we tested candidate gene-sets relevant to PD. Results: No mutations common to the cases with PD were found in the intact 22q11.2 region. While all were negative for rare mutations in a gene-set comprising PD disease-causing and risk genes, another candidate gene-set of 1000 genes functionally relevant to PD presented a nominally significant (P=0.03) enrichment of rare putatively damaging missense variants in the PD cases. Polygenic score results, based on common variants associated with PD risk, were non-significantly greater in those with PD. Conclusions: The results of this first-ever pilot study of WGS in PD suggest that the cumulative burden of genome-wide sequence variants may contribute to expression of early-onset PD in the presence of threshold-lowering dosage effects of a 22q11.2 deletion. We found no evidence that expression of PD in 22q11.2DS is mediated by a recessive locus on the intact 22q11.2 chromosome or mutations in known PD genes. These findings offer initial evidence of the potential effects of multiple within-individual rare variants on the expression of PD and the utility of next generation sequencing for studying the etiology of PD.

Study EGAS00001002275

Multi-Layered Molecular Profiling Informs the Diagnosis and Targeted Therapy of Desmoplastic Small Round Cell Tumor

Desmoplastic small round cell tumor (DSRCT) is an ultra-rare soft tissue sarcoma with dismal prognosis and limited therapeutic options. Here, we identified individual therapeutic targets in a cohort of 30 DSRCT patients by use of multilayered molecular profiling (genomic, transcriptomic, and (phospho-)proteomic analyses) within the DKFZ NCT DKTK MASTER program. 28 of the 30 patients received at least one personalized treatment recommendation. In line with the observation that DSRCTs are an entity with a low mutational load, most of the treatment recommendations were based on target gene or -protein expression. 13 patients received a recommended therapy, resulting in four partial responses, four stable diseases, and five progressive diseases. Eight of 13 patients received the multi tyrosine kinase inhibitor pazopanib, with mixed response. Based on high ERBB2 expression levels, two patients received the antibody drug conjugate trastuzumab deruxtecan (T-DXd), showed exceptional, long-lasting responses and are currently (8/2024) still on treatment and progression-free. We detected high ERBB2 protein expression in all DSRCT samples with proteome data (n = 9), which suggests that ERBB2 may be a particularly useful therapeutic target for DSRCT patients. For one of the two exceptional responders we did not detect a substantial ERBB2/HER2 overexpression, which may indicate that T-DXd could also be efficacious in DSRCT patients with low levels of HER2, as also described for breast cancer patients. In conclusion, multi-omic molecular profiling revealed personalized therapeutic targets in patients with DSRCT with HER2 targeting, with T-DXd as a promising therapeutic option for this ultra-rare and difficult-to-treat tumor entity.

Study EGAS00001007934

Whole genome sequencing of tumor samples from advanced pre-treated NSCLC patients recruited in a phase I/II single-arm trial utilising CVA21, an oncolytic coxsackie virus, in combination with pembrolizumab.

Acquired or de novo resistance to immune checkpoint inhibitors occurs in the majority of advanced non-small cell lung cancers (NSCLC). There is an unmet need to improve outcomes for patients with this condition. Oncolytic viruses represent an attractive treatment approach due to their dual activity in inducing tumor cell lysis directly and potentially augmenting anti-tumor immunity. Here, we present the safety, efficacy and translational findings from a phase I/II single-arm trial utilising CVA21, an oncolytic coxsackie virus, in combination with pembrolizumab in patients with advanced pre-treated NSCLC. We performed paired biopsies pre- and post-treatment in 10 patients who received intravenous CVA21 and pembrolizumab, 8 of whom had prior treatment with immune checkpoint inhibitor (ICI) therapy. Whole genome sequencing and spatial proteomics were performed to comprehensively characterise the response to CVA21. Combination CVA21/pembrolizumab (anti PD-1) therapy was well tolerated with no serious treatment-related adverse events. Partial responses were seen in two patients with prior acquired anti-PD-1 resistance and disease stabilisation in 6 patients, giving a clinical benefit rate of 80%. High baseline tumor mutational burden and PD-L1 expression were observed in patients with better response to treatment. Interestingly, an increase in antigen presentation and CD8+ T cell infiltration was observed on treatment compared with baseline in patients with better progression-free survival. This study demonstrates the potential of CVA21 to modulate the immunogenicity of tumor cells and remodelling the tumor microenvironment, providing novel insights for patient selection for trials involving novel immunotherapeutic approaches.

Study EGAS00001008258

Sequencing of paediatric High Grade Gliomas and DIPG

Paediatric high grade glioma (pHGG) and diffuse intrinsic pontine glioma (DIPG) are a phenotypically diverse collection of tumours which molecular profiling studies have refined into age- and location-based subgroups driven by unique genetic and epigenetic alterations. Despite this, the rarity of these tumours means individual studies are underpowered to investigate subgroup-specific changes and to identify low frequency events recurrently targeting specific pathways. We have collated genomic data from 143 unpublished cases up to 30 years of age along with those from 18 previously published datasets in an integrated analysis comprising 1067 unique patients across all anatomical compartments of the CNS, with histone mutation status available for 893 cases (n=67 H3.3G34R/V, n=316 H3.3K27M, n=68 H3.1/H3.2K27M). Distinct clinicopathological and molecular subgroups are defined by these histone mutations, with several novel co-segregating mutations identified through these aggregated analyses, including loss of FBXW7 in H3.3G34R/V, TOP3A rearrangements in H3.3K27M, and BCOR mutations in H3.1K27M. Histone wild-type subgroups can be further refined by the presence of key oncogenic events (PDGFRA or EGFR amplification/mutations) or methylation profiles which molecularly more closely resemble pleomorphic xanthoastrocytoma or other low grade gliomas. The burden of genomic aberrations decreases with age, highlighting the infant population as biologically and clinically distinct. Across the whole cohort, we identify novel previously unrecognised pathway dysregulation in small subsets of tumours (e.g. splicing, WNT, immune response) which make up the umbrella classification of paediatric diffusely infiltrating gliomas. The integrated dataset helps to further define the molecular diversity of the disease, opening up novel avenues for biological study and providing a basis for functionally defined future treatment stratification.

Study EGAS00001002314

Suppressors and activators of JAK-STAT signaling at diagnosis and relapse of acute lymphoblastic leukemia in Down Syndrome

Children with Down Syndrome (DS) are prone to development of high risk B-cell precursor ALL (DS-ALL) that differs genetically from most sporadic pediatric ALLs. Increased expression of CRLF2, the receptor to thymic stromal lymphopoietin (TSLP) characterizes about half of DS-ALLs and also a subgroup of sporadic “Philadelphia-like ALLs”. To understand the pathogenesis of relapsed DS-ALL we performed integrative genomic analysis of 25 matched diagnosis remission and relapse DS-ALLs. We found that the CRLF2 rearrangements are early events during DS-ALL evolution and generally stable between diagnoses and relapse. Secondary activating signaling events in JAK-STAT/RAS pathway were ubiquitous but highly redundant between diagnosis and relapse, suggesting that signaling is essential but that no specific mutations are “relapse driving”. We further found that activated JAK2 may be naturally suppressed in 25% of CRLF2 positive DS-ALL by loss-of-function aberrations in USP9X, a deubuitinase previously shown to stabilize the activated phosphorylated JAK2. Interrogation of large ALL genomic databases extended our findings up to 25% of CRLF2 positive, “Ph like” ALLs. Pharmacological or genetic inhibition of USP9X as well as treatment with low-dose Ruxolitinib enhanced the survival of pre-B ALL cells overexpressing mutated JAK2. Thus, somehow counterintuitive, we found that suppression of JAK-STAT “hyper-signaling” may be beneficial to leukemic B cell precursors. This and the reduction of JAK mutated clones at relapse suggest that the therapeutic effect of JAK specific inhibitors may be limited. Rather, combined signaling inhibitors or direct targeting of the TSLP receptor may be a useful therapeutic strategy for DS-ALL.

Study EGAS00001002410

Genomic diversity of the African-descent Makranis of Pakistan

From the eighth century onward, the Indian Ocean was the scene of extensive trade of sub-Saharan African slaves, via sea routes controlled by Muslim Arab and Swahili traders. Several populations in present-day Pakistan and India are thought to be the descendants of such slaves, yet their genetic ancestries and history of admixture and natural selection remain largely undefined. Here we studied the genome-wide diversity of the African-descent Makranis, residing on the Arabian Sea coast of Pakistan, together with that of four neighboring Pakistani populations, to investigate the genetic legacy, population dynamics and tempo of the Indian Ocean slave trade. We show that the Makranis result from an admixture event between local Baluch tribes and Bantu-speaking populations from eastern/southeastern Africa, which we dated to ~300 years ago during the Omani Empire domination. Levels of inbreeding, measured through cumulative runs of homozygosity, were found to be similar across Pakistani populations, suggesting that the admixed Makranis rapidly adopted the traditional practice of endogamous marriages. Finally, we searched for signatures of post-admixture selection at several traits evolving under positive selection, including skin color, lactase persistence and resistance to malaria. We demonstrate that the African-specific Duffy-null blood group — believed to confer resistance against Plasmodium vivax infection — was recently introduced to Pakistan through the slave trade and evolved under positive selection in this vivax malaria-endemic region. Collectively, our study reconstructs the genetic and adaptive history of a neglected episode of the African Diaspora, and illustrates the impact of recent admixture on the diffusion of adaptive traits across human populations.

Study EGAS00001002558

Ultra-Fast Patient-Derived Xenografts Identify Functional and Spatial Tumour Heterogeneities that Drive Therapeutic Resistance

Treatment of patients with advanced cancers increasingly relies on agents targeting specific molecular or cellular aberrations. These expensive therapies often fail because of pre-existing drug resistance. Intra-tumoural molecular and micro-environmental heterogeneity resistance can create resistance in one or more subclones of a tumour, leading to spatio-functional heterogeneity. We created a rapid approach for quantifying tumour functional phenotypes based on tumour-implantation into the chorioallantoic membrane of chick embryos. We apply this technique to 43 separate biopsies from 6 patients with renal cell carcinoma (RCC) with 36 replicates per biopsy, achieving a 93.6% engraftment rate (1449/1548). Challenging these models with sunitinib, an anti-angiogenic therapy frequently used in RCC allowed us to quantify spatio-functional heterogeneity in individual patients. Prediction of de novo drug resistance in advanced cancer patients would facilitate precision medicine by tailoring systemic therapies and improve clinical trial design outcome measures. We describe a patient-derived xenograft platform where fresh patient tumour samples are implanted into the chorioallantoic membrane (CAM) of chick embryos at high engraftment efficiencies, permitting large scale “tumour avatar” studies to guide the selection of drugs and to predict clinical drug resistance outcomes in renal cell carcinoma (RCC) patients. Applying this “tumour avatar” model in a co-clinical trial on RCC patients, we observe intra-tumoural functional heterogeneity with sunitinib treatment, revealing the presence of drug-resistant clones prior to the initiation of therapy. Multi-regional exome Matched DNA sequencing of the primary tumour and metastases revealed DNA mutational signatures associated with sunitinib resistance at baseline and after treatment. These These studies demonstrate a rapid and affordable method of findings confirm the existence of genetic and surveying functional tumour heterogeneities heterogeneity.

Study EGAS00001002627

The light chain IgLV3-21 defines a new poor prognostic subgroup in Chronic Lymphocytic Leukemia: results from a multicenter study

Purpose: Unmutated (UM) immunoglobulin heavy chain region (IgHV) status or IgHV3-21 gene usage is associated with poor prognosis in chronic lymphocytic leukemia (CLL). Interestingly, IgHV3-21 is often co-expressed with light chain IgLV3-21 which is potentially able to trigger cell-autonomous BCR-mediated signaling. However this light chain has never been characterized independently of the heavy chain IgHV3-21.Experimental Design: we performed total RNA sequencing to characterize IgLV3-21 patients (n=32) and investigated IgLV3-21 prognostic impact in terms of treatment-free (TFS) and overall (OS) survival in 3 other independent cohorts for a total of 813 patients where IgLV3-21 presence was tested by real-time PCR and confirmed by Sanger sequencing.Results: Using total RNA sequencing to characterize 32 patients with high-risk CLL, we found a high frequency (28%) of IgLV3-21 rearrangements. Gene set enrichment analysis revealed that these patients present higher expression of genes responsible for ribosome biogenesis and translation initiation (P<0.0001) and MYC target genes (P=0.0003). Patients with IgLV3-21 rearrangements displayed a significantly shorter TFS and OS (P<0.05), in particular in IgHV mutated patients. In each of the 3 independent validation cohorts (n=813), we showed that IgLV3-21 rearrangements -similar to UM IgHV status- conferred poor prognosis compared to mutated IgHV (P<0.0001). Importantly, we confirmed in multivariate analysis that this was independent of IgHV mutational status or subset #2 stereotyped receptor (P<0.0001). Conclusions: We therefore demonstrate for the first time that a light chain can impact on CLL prognosis and that IgLV3-21 light chain usage defines a new poor prognostic subgroup in CLL.

Study EGAS00001002894

Exploiting immune cell receptor information to quantify index switching in single cell transcriptome sequencing experiment

By offering high sequencing speed and ultra-high-throughput at a low price, Illumina next-generation sequencing platforms has been widely adopted in recent years. However, an experiment with multiplexed library could be at risk of molecular recombination, known as “index switching”, which causes certain amount of reads to be assigned to an incorrect library. It is reported that a new advance, exclusion amplification (ExAmp) in conjunction with the patterned flow cell technology introduced on HiSeq 3000/HiSeq 4000/HiSeq X sequencing systems, potentially suffers from a higher rate of index switching than conventional bridge amplification. We took advantage of the diverse but highly cell-unique expression of immune cell receptors to quantify index switching on single cell RNA-seq data that were sequenced on HiSeq 3500 and HiSeq 4000. By utilizing the unique T-cell receptor (TCR) expression, we could quantify the spread-of-signal from many different wells (n=51 from total three batches) due to index switching. We used TraCer to reconstruct full-length TCR from all samples, and then used Kallisto to quantify TCR gene expression. We found index switching in all three batches of experiments investigated. The median percentage of incorrectly detected markers was estimated to be 4.2% (interquartile range (IQR): 2.0%-8.7%). We did not detect any consistent pattern of certain indices to be more prone for switching than others, suggesting that index switching is a stochastic process. We confirm that index switching is a problem that affects all samples run in multiplexed libraries on Illumina HiSeq 3500 and HiSeq 4000 platforms.

Study EGAS00001002911

A genome scan for genes underlying adult body size differences between Central African hunter-gatherers and farmers.

The evolutionary and biological bases of the Central African pygmy phenotype, a characteristic of rainforest huntergatherers defined by reduced body size compared with neighboring farmers, remain largely unknown. Here, we perform a joint investigation in Central African hunter-gatherers and farmers of adult standing height, sitting height, leg length, and body mass index (BMI), considering 358 hunter-gatherers and 169 farmers with genotypes for 153,798 SNPs. In addition to reduced standing heights, hunter-gatherers have shorter sitting heights and leg lengths and higher sitting/standing height ratios than farmers and lower BMI for males. Standing height, sitting height, and leg length are strongly correlated with inferred levels of farmer genetic ancestry, whereas BMI is only weakly correlated, perhaps reflecting greater contributions of non-genetic factors to body weight than to height. Single- and multi-marker association tests identify one region and eight genes associated with hunter-gatherer/farmer status, and 24 genes associated with the height-related traits. Many of these genes have putative functions consistent with roles in determining their associated traits and the pygmy phenotype, and they include three associated with standing height in non-Africans (PRKG1, DSCAM, MAGI2). We find evidence that European height-associated SNPs or variants in linkage disequilibrium with them contribute to standing- and sitting-height determination in Central Africans, but not to the differential status of hunter-gatherers and farmers. These findings provide new insights into the biological basis of the pygmy phenotype, and they highlight the potential of cross-population studies for exploring the genetic basis of phenotypes that vary naturally across populations. Pemberton et al. Human Genetics 2018 https://doi.org/10.1007/s00439-018-1902-3

Study EGAS00001002975

Evolutionary dynamics of residual disease in human glioblastoma

Glioblastoma is the most common and aggressive adult brain malignancy against which conventional surgery and chemoradiation provide limited benefit. Even when a good treatment response is obtained, recurrence inevitably occurs either locally (c.80%) or distally (c.20%), driven by cancer clones that are often genomically distinct from those in the primary tumour. Glioblastoma cells display a characteristic infiltrative phenotype, invading the surrounding tissue and often spreading across the whole brain. We hypothesised that cancer cells responsible for relapse reside in two compartments of residual disease that are left behind after treatment: the infiltrated normal brain parenchyma, and the sub-ventricular zone (SVZ). In this model, residual disease subclones diverge early during glioblastoma evolution, may remain dormant in the normal parenchyma and are responsible for recurrence. To test our hypothesis, we performed whole-exome sequencing of 69 multi- region samples collected using fluorescence-guided resection from 11 patients, including the infiltrating tumour margin (M) and the SVZ for each patient, as well as matched blood. We used a phylogenomic approach to dissect the spatio-temporal evolution of each tumour and unveil the relation between residual disease and the main tumour mass. We also analysed two patients with paired primary-recurrence samples with matched residual disease. Our results suggest that the infiltrative phenotype is an early evolutionary trait of glioblastoma, giving rise to the tumour mass detected at surgery. After treatment, the same infiltrative clone may seed the growth of a recurrent tumour, thus representing the ‘missing link’ between the primary tumour and recurrent disease. These results are consistent with recognised clinical phenotypic behaviour and suggest that more specific therapeutic targeting of cells in the infiltrated brain parenchyma may improve patient’s outcome.

Study EGAS00001003043

Genomic landscape of oral cancers (Complete Genomics WGS)

Human papillomavirus (HPV) is a necessary but insufficient cause of a subset of oral squamous cell carcinomas (OSCC) that is increasing markedly in frequency. To identify contributory, secondary genetic alterations in these cancers, we used comprehensive genomics methods to compare 149 HPV-positive and 335 HPV-negative OSCC tumor/normal pairs. Different behavioral risk factors underlying the two OSCC types were reflected in distinctive genomic mutational signatures. In HPV-positive OSCC, the signatures of APOBEC cytosine deaminase editing, associated with anti-viral immunity, were strongly linked to overall mutational burden. By contrast, in HPV-negative OSCC, T>C substitutions in the sequence context 5’-ATN-3’ correlated with tobacco exposure. Universal expression of HPV E6*1 and E7 oncogenes was a sine qua non of HPV-positive OSCC. Significant enrichment of somatic mutations was confirmed or newly identified in PIK3CA, KMT2D, FGFR3, FBXW7, DDX3X, PTEN, TRAF3, RB1, CYLD, RIPK4, ZNF750, EP300, CASZ1, TAF5, RBL1, IFNGR1 and NFKBIA. Of these, many affect host pathways already targeted by HPV oncoproteins, including the p53 and pRB pathways, or disrupt host defenses against viral infections, including interferon and nuclear factor-B signaling. Frequent copy number changes were associated with concordant changes in gene expression. Chr. 11q (including CCND1) and 14q (including DICER1 and AKT1) were recurrently lost in HPV-positive OSCC, in contrast to their gains in HPV-negative OSCC. High-ranking variant allele fractions implicated ZNF750, PIK3CA and EP300 mutations as candidate driver events in HPV-positive cancers. We conclude that virus-host interactions cooperatively shape the unique genetic features of these cancers, distinguishing them from their HPV-negative counterparts.

Study EGAS00001002393

Indonesian Genome Diversity Project

Indonesia, the world’s fourth largest country by population, comprises an archipelago of about 900 permanently inhabited islands in tropical Asia, and hosts an astonishing array of human diversity that remains largely underrepresented in modern biological surveys. The region is unique for its key role in both the early and recent evolution of Homo outside of Africa. It has one of the first traces of anatomically modern humans in Eurasia, and was the focus of extensive interaction between archaic hominins and humans – both in the likely co-existence of H. floresiensis with modern humans, and in the introgression of Denisovans into the ancestors of its present inhabitants. More recently, Indonesia was a center of the spread of Neolithic culture by Austronesian speakers. Advancing maritime technologies allowed farming populations to treat this region as a springboard to reach oceanic islands as remote as Madagascar, Hawaii, Easter Island and New Zealand. Here we present complete human genomes from healthy individuals from multiple key islands in Indonesia and Papua. The genomes were sequenced to high depth (30x). The purpose of this dataset is to address different questions of human evolution in Island Southeast Asia and beyond, including early settlement of the area by anatomically modern humans, population structure, genetic admixture and adaptation, and contacts with archaic hominins. We employed a sampling strategy to capture a broad range of diversity across the region incorporating an average of 12 samples from each of 14 islands spanning Indonesia and Papua, from Sumatra in the west to New Britain in the east. Genomes from this dataset form part of an ongoing project to describe and understand human evolutionary history in Island Southeast Asia and Eurasia in general.

Study EGAS00001003054

Histone chaperone CHAF1A promotes proliferation and tumorigenicity in gastric cancer and impacts prognosis via context-depedent regulation of gene expression

The histone chaperone CHAF1A is associated with some tumors, but its mechanisms in tumor biology remain elusive. Gastric cancer (GC) has significant heterogeneity among patients, and novel prognostic markers and therapeutic targets may improve its outcomes. The role of CHAF1A in GC is unknown. In this study, we demonstrated that CHAF1A was overexpressed in GC tissues accompanied by low expression of wild-type P53 and high expression of the proliferation marker MKI67. In vitro, CHAF1A knockdown inhibited cell proliferation, induced cell cycle arrest and promoted apoptosis, while CHAF1A overexpression had contrary effects. In vivo, CHAF1A knockdown inhibited tumorigenicity of GC cells. A microarray assay with CHAF1A inhibition and a PathScan signaling antibody array with CHAF1A overexpression showed that CHAF1A inhibited the P53 pathway, activated stress response and induced glycolytic metabolism. In 665 patients, the expression level of CHAF1A protein was an independent predictor for overall survival and disease-free survival in non-cardia GC and also had survival associations in patients with small or perineural invasion-negative tumors. Tumor mRNA analyses based on next-generation sequencing in patients indicated that CHAF1A promoted the Warburg effect by depressing oxidative phosphorylation and increasing glycolysis. Furthermore, the heterogeneous clinical significance of CHAF1A among patient subgroups was associated with its differential regulation of gene expression involved with glycolysis and cell proliferation and survival, indicating that the role of CHAF1A is depended on patient characteristics. These results reveal critical roles and clinical values for CHAF1A in GC. CHAF1A may have context-dependent effects in GC, which increases understandings for GC heterogeneity.

Study EGAS00001003064

Exome sequencing in bipolar disorder families

Bipolar disorder (BD) is a major psychiatric disorder affecting around 1% of the global population. BD is characterized by recurrent manic and depressive episodes, and has an estimated heritability of around 70% Research has identified the first BD susceptibility genes. However, the underlying pathways and regulatory networks remain largely unknown. Since research has shown that the cumulative impact of common alleles with small effect appears to explain only around 25-38% of the phenotypic variance for BD, rare variants of high penetrance may also contribute to BD risk. The present study investigated the role of rare, nonsynonymous, and potentially functional variants via whole exome sequencing in 15 BD cases from two large, multiply affected families from Cuba. The high prevalence of BD in these pedigrees renders them promising in terms of the identification of genetic risk variants with large effect sizes. In addition, SNP array data were used to calculate polygenic risk scores for affected and unaffected family members. After correction for multiple testing, polygenic risk scores for common, BD-associated genetic variants were not significantly higher in BD cases than in healthy relatives. Exome sequencing identified a total of 17 rare and potentially damaging variants in 17 genes. The identified variants were shared by all investigated BD cases in the respective pedigree. The most promising variant was located in the gene SERPING1 (p.L349F), which has been reported as a genome-wide significant risk gene for schizophrenia. A further promising variant was located in CAPN2, which has been implicated in synaptic plasticity. The present data suggest novel candidate genes for BD susceptibility, and may facilitate the discovery of disease-relevant pathways and regulatory networks.

Study EGAS00001003085

Resolving the Full Spectrum of Human Genome Variation using Linked-Reads

Large-scale population based analyses coupled with advances in technology have demonstrated that the human genome is more diverse than originally thought. Standard short-read approaches, used primarily due to accuracy, throughput and costs, fail to give a complete picture of a genome. They struggle to identify large, balanced structural events, cannot access repetitive regions of the genome and fail to resolve the human genome into its two haplotypes. Here we describe an approach that retains long range information while harnessing the power of short reads. Starting from only 1ng of DNA, we produce barcoded short read libraries. The use of novel informatic approaches allows for the barcoded short reads to be associated with the long molecules of origin producing a novel datatype known as 'Linked-Reads'. This approach allows for simultaneous detection of small and large variants from a single Linked-Read library. We have previously demonstrated the utility of whole genome Linked-Reads (lrWGS) for performing diploid, de novo assembly of individual genomes. In this manuscript, we show the utility of reference based analysis using a single Linked-Read library for full spectrum genome analysis. We demonstrate the ability of Linked-Reads to reconstruct megabase scale haplotypes and to recover parts of the genome that are typically inaccessible to short reads, including phenotypically important genes such as STRC, SMN1 and SMN2. We demonstrate the ability of both lrWGS and Linked-Read Whole Exome Sequencing (lrWES) to identify complex structural variations, including balanced events, single exon deletions, and single exon duplications. The data presented here show that Linked-Reads provide a scalable approach for comprehensive genome analysis that is not possible using short reads alone.

Study EGAS00001003121

RNA seq on 19 samples of Radiation-Induced Meningiomas

Purpose: Majority of pediatric cancers require the irradiation of the central nervous system (CNS), and as more patients survive into adulthood from improved oncological therapy the sequelae of brain radiation are increasing in prevalence. Radiation-induced meningiomas (RIMs), one such secondary effect, demonstrate a clinically more aggressive behaviour than sporadic meningiomas (SMs). We aimed to describe the genomic mutational landscape of RIMs.Methods: We analyzed a principal cohort of 18 RIMs, with 31 RIMs overall, from patients who received childhood radiation therapy and 30 SMs, as a comparator population. We performed a multiplatform integrative genomic analysis; including methylation, whole exome and RNA sequencing. Results: RIMs exhibited a five-fold increase in copy number alterations, commonly the loss of chromosome 1p (17/18 RIMs) and 22q (17/18 RIMs), which was significantly more than observed in sporadic meningiomas. Furthermore, RNA sequencing data revealed an NF2 gene fusion event in 35.3% of RIMs. In all 6 cases, there was a complete NF2 exon spliced into a complete exon of a reciprocal gene, suggesting that the breakpoints of genomic rearrangement are intronic. All tumours with the NF2 fusion also possessed monosomy of chromosome 22q, rendering the cells with homozygous disruption of NF2. Clinically, RIMs with the NF2 fusion exhibited ill-defined borders and a tendency to develop in anatomic frontal location. Also, targeted sequencing panel confirmed that RIMs had fewer nonsynonymous NF2 mutations (6.5% vs. 30% in SM) and absence of mutations in TRAF7, SMO, KLF4, PIK3CA and AKT1, genes traditionally involved in SMs.Conclusion: Our study demonstrates that RIMs have distinct genomic drivers of oncogenesis as compared to SMs, specifically NF2 inactivation through fusion event. Radiation therapy possibly triggers genomic structural rearrangements through error-prone repair of double-stranded DNA breaks.

Study EGAS00001002318

1 Intratumoral genetic heterogeneity and clonal evolution following neoadjuvant chemoradiotherapy (nCRT) in locally advanced rectal tumors.

Purpose: Neoadjuvant chemoradiotherapy (nCRT) is the standard treatment forpatients with locally advanced rectal cancer. However, little is known about the effect of nCRT in the mutational landscape and intratumoral genetic heterogeneity (ITGH) of rectal tumors. Moreover, response to nCRT vary substantially among patients and the genetic determinants of treatment response are not fully understood. Here, we sought to determine the impact of nCRT in the mutational landscape and ITGH of rectal tumors and to identify putative drivers of therapeutic resistance. Patients and Methods: We analyzed whole-exome sequencing (WES) and clinical data from 79 primary non-treated rectal cancers obtained from The Cancer Genome Atlas (TCGA). We also compared the mutational landscape of a matched set of peripheral blood cells, irradiated tumor-adjacent colonic mucosa and pre and post-treatment tumor samples to determine the iatrogenic effects of nCRT and its impact in tissue genetic heterogeneity. Finally, we performed WES of 7 matched pre and post-treatment tumor samples to examine how nCRT affects the clonal architecture of rectal tumors. Results: We show that rectal tumors exhibit remarkable ITGH, which increases with disease progression. We also show that nCRT, per se, does not introduce novel somatic mutations, nor alters the clonal architecture of the irradiated colonic mucosa. Instead, we show that it acts as a potent selective pressure, favoring the expansion of tumor cell subpopulations more prone to resist therapy and increasing ITGH in the residual tumor. Conclusion: Our results reveal the highly heterogeneous and dynamic clonal architecture of rectal tumors undergoing nCRT. Our results also uncover putative genetic determinants of response to nCRT and a possible role for ITGH in predicting treatment outcome.

Study EGAS00001003250

Enhanced detection of circulating tumor DNA by fragment size analysis

Existing methods to improve detection of circulating tumor DNA (ctDNA) have focused on sensitivity for detecting genomic alterations but have rarely considered the biological properties of plasma cell-free DNA (cfDNA). We hypothesized that differences in fragment lengths of circulating DNA could be exploited to enhance sensitivity for detecting the presence of ctDNA and for non-invasive genomic analysis of cancer. We surveyed ctDNA fragment sizes in 344 plasma samples from 200 cancer patients using low-pass whole-genome sequencing (0.4×). To establish the size distribution of mutant ctDNA, tumor-guided personalized deep sequencing was performed in 19 patients. We detected enrichment of ctDNA in fragment sizes between 90–150 bp, and developed methods for in vitro and in silico size selection of these fragments. Selecting fragments between 90–150 bp improved detection of tumor DNA, with more than 2-fold median enrichment in &gt;95% of cases, and more than 4-fold enrichment in &gt; 10% of cases. Analysis of size-selected cfDNA identified clinically actionable mutations and copy number alterations that were otherwise not detected. Identification of plasma samples from patients with advanced cancer was improved by predictive models integrating fragment length and copy number analysis of cfDNA, with AUC&gt; 0.99 compared to AUC &lt; 0.80 without fragmentation features. Increased identification of cfDNA from patients with glioma, renal, and pancreatic cancer was achieved with AUC &gt; 0.91, compared to AUC &lt; 0.5 without fragmentation features. Fragment size analysis and selective sequencing of specific fragment sizes can boost ctDNA detection. This could reduce the required depth of cfDNA sequencing for clinical applications, earlier diagnosis and study of tumor biology.

Study EGAS00001003258

Genetic characterization of B-cell prolymphocytic leukemia: a hierarchical prognostic model involving MYC and TP53 abnormalities - WXS

B-PLL is a rare disease, whose molecular pathology is largely unknown. We report here the cytogenetic and molecular findings in a large series of 34 B-PLL. Karyotype (K) was complex (≥3 abnormalities) in 73%, and highly complex (HCK≥5) in 45%. The most frequent chromosomal aberrations were: translocation targeting the MYC gene [t(MYC)] (62%), 17p deletion including TP53 gene (38%), trisomy 18/18q (30%), 13q14 deletion (29%), trisomy 3 (24%), trisomy 12 (24%) and 8p deletion (23%). Whole-Exome Sequencing performed in 16 patients revealed recurrent mutations in TP53 (6/16, 38%), MYD88 (n=4), BCOR (n=4), MYC (n=3), SF3B1 (n=3), FAT1 (n=3), SETD2 (n=2), CHD2 (n=2), CXCR4 (n=2) and BCLAF1 (n=2). The main group of patients (21/34, 62%) had a t(MYC) associated with a higher percentage of prolymphocytes (p=0.03), CD38 expression (p<0.001), lower K complexity (p=0.0004), mutations in MYC and in genes involved in RNA metabolism and chromatin remodeling. Principal component analysis of gene expression data showed that patients with t(MYC) clustered together. A second group with MYC gain (5/34, 15%), was associated with HCK≥5 (p=0,01) and trisomy 3 (p=0,008). Altogether, 26/34 patients (76%) had a MYC activation, translocation or gain, that were mutually exclusive. We identified 3 distinct cytogenetic prognostic groups (p=0.0006): lower risk: absence of MYC activation (median not reached); intermediate risk: MYC activation without del17p (125 months); high risk: MYC activation + del17p (11 months). Our results show that cytogenetic analysis is a useful diagnostic tool in B-PLL that improves prognostic stratification.

Study EGAS00001003275

Genomic landscape of oral cancers (Illumina RNA-Seq)

Human papillomavirus (HPV) is a necessary but insufficient cause of a subset of oral squamous cell carcinomas (OSCC) that is increasing markedly in frequency. To identify contributory, secondary genetic alterations in these cancers, we used comprehensive genomics methods to compare 149 HPV-positive and 335 HPV-negative OSCC tumor/normal pairs. Different behavioral risk factors underlying the two OSCC types were reflected in distinctive genomic mutational signatures. In HPV-positive OSCC, the signatures of APOBEC cytosine deaminase editing, associated with anti-viral immunity, were strongly linked to overall mutational burden. By contrast, in HPV-negative OSCC, T>C substitutions in the sequence context 5’-ATN-3’ correlated with tobacco exposure. Universal expression of HPV E6*1 and E7 oncogenes was a sine qua non of HPV-positive OSCC. Significant enrichment of somatic mutations was confirmed or newly identified in PIK3CA, KMT2D, FGFR3, FBXW7, DDX3X, PTEN, TRAF3, RB1, CYLD, RIPK4, ZNF750, EP300, CASZ1, TAF5, RBL1, IFNGR1 and NFKBIA. Of these, many affect host pathways already targeted by HPV oncoproteins, including the p53 and pRB pathways, or disrupt host defenses against viral infections, including interferon and nuclear factor-B signaling. Frequent copy number changes were associated with concordant changes in gene expression. Chr. 11q (including CCND1) and 14q (including DICER1 and AKT1) were recurrently lost in HPV-positive OSCC, in contrast to their gains in HPV-negative OSCC. High-ranking variant allele fractions implicated ZNF750, PIK3CA and EP300 mutations as candidate driver events in HPV-positive cancers. We conclude that virus-host interactions cooperatively shape the unique genetic features of these cancers, distinguishing them from their HPV-negative counterparts.

Study EGAS00001003237

22929 results for "ega"

in 50.58 milliseconds.